流固耦合热分析实例

startear

这是我做的流回耦合热分析的例子，因为是生手，不知道结果好不好，请高手指点指点。

/prep7  

et,1,fluid142   
et,2,solid70
et,3,fluid141   
et,4,plane55
keyopt,1,4,1
keyopt,3,4,1

/view,1,1,0,0
cyl4,0,0,30,0,40,,  !建立几何形状
cyl4,0,0,0,0,80,,   
aovlap,1,2  
aplot   
agen,2,1   !流固耦合界面设置要求(个人理解)
numcmp,all
allsel  

alist
lsel,s,,,1,4
lsel,a,,,13,16  
lesize,all,4
lsel,s,,,5,8
lsel,a,,,17,20  
lesize,all,4
lsel,s,,,9,12   
lesize,all,8

allsel  

asel,s,,,1  
aatt,,,4
asel,inve   
aatt,,,3

allsel
mshape,1
mshkey,0
amesh,all   
  
allsel  
/pnum,area,1

type,2  
mat,2   
esize,,2
asel,s,,,1  
vext,all,,,0,0,8   

type,1  
mat,1   
esize,,2
asel,s,,,3,4
vext,all,,,0,0,8

esize,,4
asel,s,,,2,4
vext,all,,,0,0,-32  
allsel  

aclear,all  
etdel,3
etdel,4
numcmp,all

vsel,s,,,1  
nslv,s,1
ic,all,temp,293

sfa,5,,hflux,50000

esel,s,type,,1
nsle,s,1
ic,all,temp,293  !不知道如何设置流体的初始温度

allsel
asel,s,,,6,13   
asel,a,,,1  
nsla,s,1
sf,all,fsin,1     !设置流固界面

asel,s,,,15,18  
asel,a,,,24,27  
asel,a,,,2  
nsla,s,1
sf,all,fsin,1   

asel,s,,,19
asel,a,,,14
nsla,s,1
d,all,pres,0

nsel,s,loc,z,-32
d,all,ux,0  
d,all,uy,0  
d,all,uz,0  
d,all,vx,0  
d,all,vy,0  
d,all,vz,1.0

csys,1  
nsel,s,loc,x,80
d,all,uz,0  
d,all,vz,0  
allsel

/solu

mp,c,2,1.36
mp,dens,2,7800  
mp,kxx,2,400


fldata,solu,flow,1  
fldata,solu,ale,1   
fldata,solu,turb,1  
fldata,solu,tran,1  
fldata,solu,temp,1  
fldata,iter,exec,50
fldata,temp,nomi,293
fldata,temp,bulk,293
toffst,0

fldata,prot,dens,constant   
fldata,nomi,dens,800
fldata,prot,visc,constant   
fldata,nomi,visc,0.2
fldata,prot,cond,constant   
fldata,nomi,cond,10
fldata,prot,spht,constant   
fldata,nomi,spht,0.95  

fsan,on    !流固耦合求解设置
fsco,all,0.1
fsre,all,0.5
fsou,last   
fsit,5  
fsti,2.5  
fsdt,0.1   
fsor,solid  
fstr,solid  
fstr,fluid  
fldata,time,step,0.1  
deltim,0.1
kbc,0   

solve   
fini

/post1
file,,rfl
set,last
plnsol,temp

startear

我怎么传图上来呀

startear

chendan

楼主你选择的单元怎么二维跟三维混在一起啊？你这么做是什么意思啊？

mbrslydddd

哈哈，好，先给你加个分，不过能不能把图贴上来呢？

startear

这段时间没有上网了，试着回答chendan的问题：
采用二维单元是面积网格划分的需要，因为建立流固耦合界面的要求面上的节点对应（这是我从出错结果中理解的，也不知道对否），所以我先对面积分网然后再拉伸成实体，这样就可以做到节点的一一对应了。但在求解这前已经将这些单元删除了，所以应该不是二维和三维单元混合使用吧

startear

版主呀，我不知道怎样才能把图贴上来呀，挺郁闷的
你教一下吧

songbl123

startear，如果流固耦合界面的面上的节点不对应，会出现什么样的结果，我按照你这个例子做了一个简单的流固耦合的例子，就是没有把耦合界面的面上的节点对应，结果固体的温度好像没有传到流体上，你觉得是不是节点没有对应导致的？

startear

不好意思，前阵子忙着做实验很少上网。
这个问题我也不是很清楚，说实在我也是刚学ANSYS，有很多问题也搞不清楚

songbl123

哪位大侠能帮帮我呀
如果固体结构很不规则，无法做到边界上流固节点对应怎么办？