用ANSYS做桥墩与水的耦合作用下桥墩的自振频率，代码如下，大家看看问题在哪

ywl_5234

求解圆形空心桥墩与水耦合作用时候桥墩的自振频率

fini
/clear
/prep7
cylind,5,6,0,60
cylind,6,10,0,60
cylind,10,20,0,60
et,1,solid45
mp,ex,1,3.45e10
mp,nuxy,1,0.2
mp,dens,1,2500!2.5e-6
et,2,fluid30
et,3,fluid130
r,3,20
et,4,fluid30
keyopt,4,2,1

mp,dens,2,1030
mp,sonc,2,1460
mp,dens,3,1030
mp,sonc,3,1460
mp,dens,4,1030
mp,sonc,4,1460
!!!对桥墩的网格化
type,1
real,1
mat,1
MSHKEY,0
MSHAPE,1,3D !
csys,1
vsel,s,volu,,1
vmesh,all
!!!!!!!!!!!
!!!对桥墩外围与桥墩相接触的水进行网格化
type,2
real,2
mat,2
MSHKEY,0
MSHAPE,1,3D !
csys,1
vsel,s,volu,,2
vmesh,all
!对桥墩外围与桥墩不接触的水进行网格化
type,4
real,4
mat,4
MSHKEY,0
MSHAPE,1,3D !
csys,1
vsel,s,volu,,3
vmesh,all
! 指定无限吸收边界
csys,1
nsel,s,loc,x,20
type,3
real,3
mat,3
esurf  
! 标识流固交接面
nsel,s,loc,x,6
esel,s,type,,2
sf,all,fsi,1
allsel,all
/solu
nsel,s,loc,x,5,6
nsel,r,loc,z,0
d,all,all
antype,modal
modopt,damp,10
mxpand,10,,,yes
allsel,all
solve
finish
/post1
set,list
运行的时候出现错误：
Node no. 2620 of the Fluid130 Acoustic Absorbing element 13065 should lie on a SPhere of radius 20;instead lies at a radius of 54.156 .Ensure that the radius input through the element real constant matches with that of the absorbing boundary.
大家看看问题出在什么地方。

下面是这个圆形空心墩在空气中的自振频率计算的代码，计算的结果和在水中计算的结果完全一样，大家看看是怎么回事啊？fini
/clear
/prep7
cylind,5,6,0,60
et,1,solid45
mp,ex,1,3.45e10
mp,nuxy,1,0.2
mp,dens,1,2500!2.5e-6
SMRTSIZE,2
vmesh,all
fini
/solu
antype,modal
modopt,lanb,20
nsel,s,loc,z,0
d,all,all
mxpand,20
allsel,all
solve
finish
/post1
set,list

桥墩在水中和在空气中计算出来的频率都是一样的，如下：
*****  INDEX OF DATA SETS ON RESULTS FILE  *****
   SET   TIME/FREQ    LOAD STEP   SUBSTEP  CUMULATIVE
     1  2.1532             1         1         1
     2  2.1558             1         2         2
     3  10.069             1         3         3
     4  10.806             1         4         4
     5  10.813             1         5         5
     6  15.497             1         6         6
     7  24.587             1         7         7
     8  24.617             1         8         8
     9  26.141             1         9         9
    10  26.433             1        10        10
    11  27.790             1        11        11
    12  28.065             1        12        12
    13  30.237             1        13        13
    14  31.776             1        14        14
    15  31.873             1        15        15
    16  38.720             1        16        16
    17  38.792             1        17        17
    18  39.230             1        18        18
    19  39.295             1        19        19
    20  46.345             1        20        20

各位大侠能否看看是怎么回事呢？

[ 本帖最后由 ywl_5234 于 2009-4-21 07:51 编辑 ]