求解空气流体与周围岩体之间的热量交换问题

zyaiyyf

最近想做空气流体与周围岩体之间的热交换问题，想知道空气对周围岩体的热影响范围，但是不知道为什么算不出影响范围，求解释，命令流如下：
finish
/clear,start
/config,nproc,2
/filename,fsi,1
/title,FSI fen xi
/prep7
c=-6
et,1,fluid142       !流体三维单元
et,2,solid70        !热体三维单元
et,3,fluid141       !流体二维单元
et,4,plane55      !热学平面
keyopt,1,4,1
keyopt,3,4,1
mp,dens,1,1800
mp,c,1,1.01
mp,kxx,1,1.19
k,1,0,0
k,2,-20,20
k,3,-20,0
k,4,-20,-20
k,5,0,-20
k,6,0,20

l,2,3
l,3,4
l,4,5
l,5,1
l,1,6
l,6,2
l,1,3
circle,1,4.5
circle,1,5
circle,1,5.25

alls
lcsl,all
nsel,all
lsel,all
nummrg,all
numcmp,all

lsel,s,,,5
lsel,a,,,9
lsel,a,,,13
lsel,a,,,8
lsel,a,,,12
lsel,a,,,16
ldele,all
ksel,s,,,11,15,4
ksel,a,,,7
kdele,all
nsel,all
lsel,all
nummrg,all
numcmp,all
alls
ldele,12,18,6,1
lcomb,5,6
lcomb,11,13
alls
nummrg,all
numcmp,all

lsel,s,,,5,10,5
lesize,all,,,10
al,5,10
mshape,1
mshkey,0
aatt,,,3
amesh,1
alls
lsel,s,,,6,7
lesize,all,,,5
alls
lsel,s,,,14,15
lesize,all,,,1

alls
al,14,5,15,7,6
mshape,1
mshkey,0
aatt,1,,4
amesh,2

alls
lsel,s,,,16,18
lesize,all,,,1
alls
al,8,16,6,17
al,17,7,18,9
mshape,1
mshkey,0
aatt,1,,4
amesh,all

lsel,s,,,1,4
lsel,a,,,11,13
lesize,all,,,15
alls
al,1,11,8,13,4
al,2,11,9,12,3
aatt,1,,4
mshape,1
mshkey,0
amesh,all
alls
nummrg,all,,,,low       !先合并 在压缩
numcmp,all
arsymm,x,all           
nummrg,all,,,,low
numcmp,all
alls
agen,2,1,7,6
aclear,1,7,6
adel,1,7,6,1
agen,2,2,8,6
aclear,2,8,6
adel,2,8,6,1     !存在多余的线
type,1
asel,s,type,,3
aplot
esiz,,6
vext,all,,,0,0,c
asel,s,type,,4
aplot
type,2
vext,all,,,0,0,c
alls
aclear,all
etdel,3
etdel,4

!流体部分的左右边界
asel,s,,,17
asel,a,,,8
aplot
nsla,s,1
d,all,vx,0 !设置流体的三个方向速度为0,y\x方向位移为0
d,all,vy,0
d,all,vz,0
!d,all,uy,0
!d,all,ux,0
!流固界面的标记
!流体部分
alls
asel,s,,,8,17,9   
aplot
nsla,s,1
sf,all,fsin,1
!固体部分
alls
asel,s,,,23,43,20   
aplot
nsla,s,1
sf,all,fsin,1
!空气进口设置
asel,s,,,13,14
aplot
nsla,s,1
!d,all,ux,0 !入口处三个方向位移为0,只有z方向的速度
!d,all,uy,0
!d,all,uz,0
d,all,vx,0
d,all,vy,0
d,all,vz,-5
d,all,temp,263
!空气出口设置
asel,s,,,2,16,14
aplot
nsla,s,1
d,all,pres,0
!固体边界的位移
alls
asel,s,,,23,43,20   
aplot
nsla,s,1
!d,all,ux,0
!d,all,uy,0
!d,all,uz,0
d,all,temp,280
sf,all,conv,5,280
!围岩初始温度
alls
vsel,s,type,,2
vplot
nslv,s,1
ic,all,temp,280

alls
fini
/solu
fldata1,solu,temp,1
flda,solu,flow,1
flda,solu,ale,1
flda,solu,turb,1
!flda,solu,tran,1
flda,iter,exec,50

flda,prot,dens,air-si
flda,prot,visc,air-si
flda,prot,cond,air-si
flda,prot,spht,air-si
ACEL,,9.81
FLDATA8,NOMI,DENS,1060   
FLDATA8,NOMI,VISC,0.0424
FLDATA8,NOMI,COND,1
FLDATA8,NOMI,SPHT,4

fsan,on
!FSI analysis on
fsco,all,1.0e-1 !FSI convergence criteria
fsre,all,0.5 !FSI relaxation
fsou,1 !FSI output frequency
fsit,5 !Number of stagger loops
fsti,5 !FSI end time
fsdt,0.01 !FSI delta time
fsor,fluid !FSI order of solution
fstr,solid, !Transient fluid for FSI analysis
antype,trans
fstr,fluid, !Transient structural for FSI analysis
deltim,0.01
flda,time,step,0.01
save
solve
save
fini