finish$/clear$/prep7
et,1,solid45
et,2,fluid30,,1 !定义流体单元fluid30,表示不与结构接触的流体单元
et,3,fluid30,,0 !定义流体单元fluid30,keyopt(2)设置为0,表示流体单元与结构接触
mp,ex,1,2.1e11
mp,prxy,1,0.3
mp,dens,1,7850
mp,dens,2,1000
mp,sonc,2,1460
blc4,-10,-2.75,20,5.5,373
blc4,-60,-16.5,120,33,373
vptn,1,2 !分割体后,内部棱体1为结构域,外部棱体3为流体域
vsel,s,,,1
vatt,1,,1, !体1附上结构材料属性
vsel,s,,,3
vatt,2,,2, !体3附上流体材料属性
allsel
vglue,1,3
lesize,5,,,10
lesize,6,,,5
lesize,9,,,30
lesize,19,,,40
lesize,17,,,40
lesize,18,,,20
lesize,20,,,20
mshape,0,3D
mshkey,1
vsweep,all,13,14,0 !扫掠流体域划分流体单元
asel,s,,,3,6,1, !选择流体结构接触面
nsla,s,1 !选择这些面上的所有节点
esln,s,0, !选择这些节点相邻的单元
esel,u,type,,1 !去除1号单元(剩下的是2号单元)
emodif,all,type,3 !将接触面上的2号单元强制修改为3号单元
!sf,all,fsi,1
asel,s,,,3,6,1
sfa,all,,fsi
asel,s,,,1
nsla,s,1
d,all,all !悬臂梁底部全部约束
asel,s,,,13
nsla,s,0
d,all,uz,0 !流体单元z向位移约束
asel,s,,,9,12,1 !流体域其他面给定约束pres为0
nsla,s,1
d,all,pres,0.0
asel,s,,,14
nsla,s,1
d,all,pres,0.0
allsel
finish
/solu
antype,modal
modopt,unsym,10,,,
mxpand,10,1,1000
solve
|