这段时间在研究单纯利用ansys中的单元进行流固耦合,发现这方面的资料很少而且很零散,例子也很少。我看了很多资料琢磨了半天,做了下面这个例子。这个例子利用shell181与fluid80单元进行直接耦合自由度的方式进行计算。没有用173和170来做接触。大家讨论下,看看我这样做有没有什么问题。大家多交流争取能讨论出来一个正确常用的例子,供后来者用。需要注意fluid80单元适用于没有进出口的流固耦合分析,同时其需要建在z坐标以下,自由液面放置在z=0的位置。
fini
/cle
/prep7
!***********单元及材料建立
ET,1,shell181
ET,2,FLUID80
KEYOPT,2,2,1
MP,EX,1,2E11
MP,PRXY,1,0.3
MP,DENS,7.85E3
R,1,0.006
MP,EX,2,3E8
MP,DENS,2,1E3
MP,VISC,2,0.001
!************建立罐体
BLOCK,0,1,0,1,0,-5
alls
VDELE,1
ASEL,S,LOC,Z,0
ADELE,ALL
ALLS
esize,0.1
!***********划分网格
AATT,1,,1
ASEL,U,LOC,Z,0
AMESH,ALL
ALLS
!***********建立流体域
BLOCK,0,1,0,1,0,-5
!**********划分网格
esize,0.1
vatt,2,,2
vsweep,all
ALLS
fini
/SOL
!********静力学求解
ANTYPE,0
!**********施加重力载荷
ACEL,1,-1,9.8,
ALLS
!*********耦合垂直罐体方向节点自由度
NSEL,S,LOC,X,0
NSEL,A,LOC,X,1
CPINTF,UX,0.0001,
ALLS
NSEL,S,LOC,Y,0
NSEL,A,LOC,Y,1
CPINTF,UY,0.0001,
ALLS
NSEL,S,LOC,Z,-5
CPINTF,UZ,0.0001,
ALLS
!************约束罐体底面
ESEL,S,TYPE,,1
ALLSEL,BELOW,ELEM
NSEL,S,LOC,Z,-5
D,all,ALL,
ALLS
SOLVE
FINISH
附结果图如下:
![利用SHELL181,FLUID80单元进行流固耦合计算]() |