流固耦合实例（附命令流）

zj300600

跟大家一块分享一个流固耦合的实例！

finish
/clear
/PREP7
et,1,solid65 !实体单元
et,2,fluid30,,1!!!流体单元
et,3,fluid30!!!!流固耦合单元
!定义材料
mp,ex,1,3.45e7
mp,dens,1,2.6
mp,prxy,1,0.1667
mp,dens,2,1000
mp,sonc,2,1400
!建立桥墩
!建立外墩
k,1,3,-3.3
k,2,3,3.3
k,3,-3,3.3
k,4,-3,-3.3
k,5,3.7,-4.0,-70
k,6,3.7,4.0,-70
k,7,-3.7,4.0,-70
k,8,-3.7,-4.0,-70
v,1,2,3,4,5,6,7,8
*GET,kpM,kp,0,NUM,MAX
!建立内墩
K,kpM+1,2.1,2.4
K,kpm+2,2.1,-2.4
k,kpm+3,-2.1,-2.4
k,kpm+4,-2.1,2.4
K,kpM+5,2.8,3.1,-70
K,kpm+6,2.8,-3.1,-70
k,kpm+7,-2.8,-3.1,-70
k,kpm+8,-2.8,3.1,-70
v,kpm+1,kpm+2,kpm+3,kpm+4,kpm+5,kpm+6,kpm+7,kpm+8
vsbv,1,2
!建立流体单元
!建立实体流体单元
*GET,kpM,kp,0,NUM,MAX
K,kpM+1,31.4,-34.4,-14
K,kpm+2,31.4,34.4,-14
k,kpm+3,-31.4,34.4,-14
k,kpm+4,-31.4,-34.4,-14
K,kpM+5,37,-40,-70
K,kpm+6,37,40,-70
k,kpm+7,-37,40,-70
k,kpm+8,-37,-40,-70
v,kpm+1,kpm+2,kpm+3,kpm+4,kpm+5,kpm+6,kpm+7,kpm+8
!建立内部流体
*GET,kpM,kp,0,NUM,MAX
K,kpM+1,3.14,3.44,-14
K,kpm+2,3.14,-3.44,-14
k,kpm+3,-3.14,-3.44,-14
k,kpm+4,-3.14,3.44,-14
K,kpM+5,3.7,4.0,-70
K,kpm+6,3.7,-4.0,-70
k,kpm+7,-3.7,-4.0,-70
k,kpm+8,-3.7,4.0,-70
v,kpm+1,kpm+2,kpm+3,kpm+4,kpm+5,kpm+6,kpm+7,kpm+8
vsbv,1,2
!划分网格
wplane,,0,2.4,0,-2.1,2.4,0,-2.8,3.1,-70
vsbw,all
wplane,,0,-2.4,0,2.1,-2.4,0,2.8,-3.1,-70
vsbw,all
wplane,,-2.1,0,0,-2.1,-2.4,0,-2.8,-3.1,-70
vsbw,all
wplane,,2.1,0,0,2.1,2.4,0,2.8,3.1,-70
vsbw,all
wplane,,0,3.44,-14,-3.14,3.44,-14,-3.7,4,-70
vsbw,all
wplane,,0,-3.44,-14,3.14,-3.44,-14,3.7,-4,-70
vsbw,all
wplane,,-3.14,0,-14,-3.14,-3.44,-14,-3.7,-4,-70
vsbw,all
wplane,,3.14,0,-14,3.14,3.44,-14,3.7,4,-70
vsbw,all
wplane,,3.14,3.44,-14,-3.14,3.44,-14,-3.14,-3.44,-14
vsbw,all
vsel,s,loc,x,-3.7,3.7
vsel,r,loc,y,-4,4
vatt,1,,1
vsel,inve
vatt,2,,2
vsel,all
esize,2
vmesh,all
!流体无限远处
allsel
asel,s,loc,x,-37,-31.4
asel,a,loc,x,31.4,37
asel,a,loc,y,-40,-34.4
asel,a,loc,y,34.4,40
nsla,s,1
d,all,pres,0
!墩的底部
allsel
asel,s,loc,z,-70
asel,r,loc,x,-3.7,3.7
asel,r,loc,y,-4,4
nsla,s,1
d,all,all
!定义流固耦合截面
allsel
asel,s,,,24,34,10
asel,a,,,49,59,10
asel,a,,,141,144,3
asel,a,,,159,169,10
asel,a,,,128,152,24
asel,a,,,33,60,27
asel,a,,,73,77,4
asel,a,,,100,150,50
asel,a,,,91,139,48
asel,a,,,86,96,10
asel,a,,,157,164,7
asel,a,,,170,176,6
nsla,s
esln,s
esel,r,type,,2
emodif,all,type,3
sf,all,fsi,1
finish
!!!!!!!!!!
!!!求解
/solu
acel,,9.8
antype,0
allsel
solve
finish
/solu
antype,modal
modopt,unsym,100
mxpand,100,,,yes
allsel
solve