一个涉及设备内两种气体传热分析的问题，能量传递过程包括流体热辐射、对流传热等

iamsisi

大家帮我看看这个命令流，为什么感觉边界条件除了温度，其他的都没有加载上？
问题出在哪里？而且好像如果我用GUI的话，边界条件就加载上了。
r0=0.025
r1=0.045
r2=0.2032
r3=0.2112
r4=2.0
!rad1=0.35
!rad2=0.5
!rad3=0.8
SAVE


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!model
/PREP7  
PCIRC,r1,r0,0,360
PCIRC,r2,r1,0,360
PCIRC,r3,r2,0,360
PCIRC,r4,r3,0,360
allsel,all
aglue,all
lglue,all
numcmp,area
numcmp,line
save
/pnum,line,1
/pnum,area,1
/replot

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!material
et,1,fluid141
keyopt,1,1,2
mp,dens,2,2791
mp,kxx,2,202.4
mp,c,2,871

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!material attribute
allsel,all
asel,s,area,,1,,,1
aplot
aatt,2,,1,0,
ASUM, DEFAULT     
*GET,s1,AREA,1,AREA
allsel,all
asel,s,area,,3,,,1
aplot
aatt,2,,1,0,
ASUM, DEFAULT     
*GET,s2,AREA,3,AREA
allsel,all
/replot
SAVE

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!mesh
esize,0,40
mshape,1,2d
mshkey,0
asel,s,area,,4,,,1
amesh,all

allsel,all
esize,0,80
mshape,1,2d
mshkey,0
asel,s,area,,3,,,
asel,a,area,,2,,,
asel,a,area,,1,,,
amesh,all
SAVE

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!load
allsel,all
/replot

lsel,s,,,17,20,1,1
lplot                                                !temperature
nsll,s,1
d,all,temp,298
save


allsel,all                          !velocity
LPLOT   
lsel,s,,,all
nsll,s,1
d,all,vx,0
d,all,vy,0
save

allsel,all                          !source
asel,s,,,         1
aplot
bfa,all,hgen,18.966/s1
save

allsel,all
asel,s,,,         3
aplot
bfa,all,hgen,12.532/s2
SAVE

allsel,all                          !surface-surface radiation
ASEL,S, , ,       2
aplot
ESLA,
nsle,
lsel,s,,,1,4,1
nsll,s,1
sf,all,rdsf,0.35,-2
lsel,s,,,9,12,1
nsll,s,1
sf,all,rdsf,0.5,-2
save

allsel,all                          !ambient radiation
asel,s,,,         4
esla,s
nsle,s
lsel,s,,,13,16,1
nsll,s,1
sf,all,rad,0.8,298
save

stef,5.67e-8                            !radiation parameters                                                       

! Set Stefan-Boltzmann constant
toffset,0                            ! Set temperature offset to zero
radopt,.5,0.00001               

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!fluid property
allsel,all
SAVE
FLDATA1,SOLU,TRAN,0
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TEMP,1
FLDATA1,SOLU,TURB,1
FLDATA1,SOLU,COMP,0
FLDATA1,SOLU,VOF,0  
FLDATA1,SOLU,SFTS,0
FLDATA1,SOLU,IVSH,0
FLDATA1,SOLU,SWRL,0
FLDATA1,SOLU,SPEC,0
FLDATA1,SOLU,ALE,0  
FLDATA1,SOLU,RDSF,1


/COM,,Steady State Analysis,0   
FLDATA2,ITER,EXEC,1000,
FLDATA2,ITER,OVER,50,   
FLDATA2,ITER,APPE,50,   
FLDATA3,TERM,VX,0.01,   
FLDATA3,TERM,VY,0.01,   
FLDATA3,TERM,VZ,0.01,   
FLDATA3,TERM,PRES,1e-008,   
FLDATA3,TERM,TEMP,1e-008,   
FLDATA3,TERM,ENKE,0.01,
FLDATA3,TERM,ENDS,0.01,
FLDATA5,OUTP,SUMF,10,   


FLDA,PROT,DENS,CMIX         ! Fluid density
FLDA,PROT,VISC,CMIX         ! Fluid viscosity
FLDA,PROT,COND,CMIX         ! Fluid conductivity
FLDA,PROT,SPHT,CMIX         ! Fluid specific heat
flda,vary,dens,t            !all properties MUST be variable
flda,vary,visc,t
flda,vary,cond,t
flda,vary,spht,t

ACEL,0,9.81,0,  

FLDATA18,METH,PRES,1
FLDATA19,TDMA,PRES,100,
FLDATA18,METH,TEMP,6
FLDATA22,MAXI,TEMP,1000,
FLDATA20,SRCH,TEMP,2,   
FLDATA20B,PBCG,FILL,6,  
FLDATA21,CONV,TEMP,1e-012,

KEYOPT,1,1,2,   
MSDATA,2,8314.29,   


MSSPEC,1,air,29,1,  
MSQUAD,1,0,0,   
MSNOMF,1,0.5,   
MSPROP,1,DENS,gas,1.166,298,1.01325e5,0,  
MSVARY,1,DENS,1
MSPROP,1,VISC,gas,1.72e-5,273,80.469,0,  
MSVARY,1,VISC,1
MSPROP,1,COND,gas,0.0244,273,180.67,0,  
MSVARY,1,COND,1
MSPROP,1,SPHT,CONSTANT,1005,   

MSSPEC,2,SF6,140.07,1,  
MSQUAD,2,0,0,   
MSNOMF,2,0.5,   
MSPROP,2,DENS,gas,6.14,298,1.2367e5,0,  
MSVARY,2,DENS,1
MSPROP,2,VISC,gas,1.45e-5,275.11,267.344,0,  
MSVARY,2,VISC,1
MSPROP,2,COND,gas,0.0127,273,264.3,0,  
MSVARY,2,COND,1
MSPROP,2,SPHT,CONSTANT,688.96,  
save

allsel,all
ASEL,S, , ,       2
aplot
ESLA,S
NSLE,S
ic,all,sf6,1.0
ic,all,air,0.0

ASEL,S, , ,       4
aplot
ESLA,S
NSLE,S
ic,all,sf6,0.0
ic,all,air,1.0

allsel,all
save
/solu
solve

iamsisi

其中图是借用了以前fossile兄的图，做了标注和修改，我的模型跟fossile兄的有相似之处，希望不要见怪
参数及边界条件如下：
外壳
外半径： 211.2mm
厚度： 8mm，
内半径为：203.2

导体
外半径： 45mm
内半径： 25m

大空间
取半径1.2m

外壳
外表面外部辐射系数：0.8
参考温度：25℃

内表面surface-to-surface辐射系数
0.5

导体
外表面surface-to-surface辐射系数
0.3

气体
空气（理想空气） SF6气体


其中外部空气的压强为0.1MPa，内部SF6气体的压强为0.45MPa


外壳和导体材料为铝

导体生热率(导体面积为a)：18.966/a w/m3

重力y轴方向：9.81

其中
图中红色部分为导体

黄色部分为外壳（图仅是示意，不代表实际情况）

iamsisi

能关注一下我的帖子，给我一些指导，谢谢各位。