这是用流固耦合方式求解
- /prep7
- et,1,fluid142
- et,2,solid70
- et,3,fluid141
- et,4,plane55
- ! 定义了五种单元类型cyl4,0,0,0,0,10 !建立几何形状
- cyl4,0,0,0,0,80 ! 画了一个圆面和一个圆环
- /pnum,line,1 显示实体编号
- /pnum,area,1
- /replot
- aovlap,1,2 !面搭接
- aplot !画面
- numcmp,all !压缩编号
- allsellsel,s,,,1,8 选择线
- lesize,all,6 对于所选线段进行网格划分设置asel,s,,,1 选择面
- aatt,,,3 设置所选面的单元类型为(用于未划分网格的面)
- asel,inve 选择除以外的所有面
- aatt,,,4mshape,1 使用三角形网格
- mshkey,0 使用自由网格划分
- alls
- amesh,allagen,2,1 !通过复制原有的圆面,产生一个新的重合的圆面
- aclear,1 ! 清除网格
- adel,1,,,1 !删除未划分的面积区域allseltype,2 !为随后生成的单元设置单元号和材料号
- esize,,10 !设置线分段数
- asel,s,,,2
- vext,all,,,0,0,80 !对于选择的面拉伸生成实体type,1
- esize,,10
- asel,s,,,3
- vext,all,,,0,0,80allsel
- aclear,all !删除所有面上的节点和网格
- etdel,3 !删除单元类型
- etdel,4
- numcmp,all
- asel,s,,,13,16
- nsla,s,1
- d,all,vx,0
- d,all,vy,0
- d,all,vz,0 !设置边界的无滑移边界条件
- SF,ALL,FSIN,1 !设置流固中流体的界面asel,s,,,12
- nsla,s,1
- d,all,temp,310 !设置入口温度
- d,all,vx,0.2 !设置入口速度
- d,all,vy,0
- d,all,vz,0
- asel,s,,,3
- nsla,s,1
- d,all,pres,0 !设置出口压力asel,s,,,8,11
- nsla,s,1
- SF,ALL,FSIN,1 !设置流固中固体的界面asel,s,,,1
- nsla,s,1
- d,all,temp,293
- asel,s,,,4,7
- nsla,s,1
- sf,all,conv,3,293
- allselfini
- /solutoffst,0
- mp,kxx,1,0.67 定义固体材料FLDA,SOLU,TEMP,T
- FLDA2,ITER,EXEC,10
- FLDATA8,NOMI,DENS,1060
- FLDATA8,NOMI,VISC,0.0424
- FLDATA8,NOMI,COND,1
- FLDATA8,NOMI,SPHT,4
- FLDATA34,mir,MOME,0.5fsan,on !流固耦合求解设置
- fsor,solid
- fstr,solid,stat
- fstr,fluid,stat
- fsin,CONS
- fsit,1
- fsdt,0.1
- fsco,all,0.1
- fsre,all,0.5
- fsou,last
- fldata,time,step,0.1
- deltim,0.1SAVESOLVE
复制代码
|