『求助』流体分析设置问题

xlpe

我用ansys分析一个发热棒在一个管子中的温度变化，管子埋在土壤中，管子和发热棒中间为空气。     我先给发热棒加一个负载进行稳态分析（不考虑发热棒与管子的辐射，只考虑管子内空气的对流），然后改变这个负载在稳态基础上进行瞬态分析，为什么瞬态分析第一步的初始条件不是以稳态结果开始的？是不是我的设置有问题，请高手指导下，命令如下。

/FILNAME,fare
/TITLE,TEMPERATURE DISTRIBUTION IN A SHORT SOLID CYLINDER-2D SOLUTION
KEYW,PR_CFD,1
/PREP7
toffst,273
ET,1,FLUID141
PCIRC,0,0.075,0,360               !以工作平面原点为圆心生成圆形区域，半径7.5cm，排管内径
PCIRC,0,0.02,0,360              !发热棒半径2cm
CYL4,0,0,0.075,0,0.085,360        !圆心坐标，排管内外径，7.5-8.5cm
RECTNG,-2,2,-2,0.5
AOVLAP,ALL
NUMMERG,ALL
NUMCMP,ALL
ALLSEL,ALL
MAT,2        !发热棒材料为2
Smrt,1
Mshape,1,2D      !三角形单元;
Mshkey,0         !自由剖分
AMESH,1
MAT,1        !排管空气材料为1
Smrt,2
Mshape,1,2D      !三角形单元;
Mshkey,0         !自由剖分
AMESH,2
MAT,3        !排管材料为3
Smrt,3
Mshape,1,2D      !三角形单元;
Mshkey,0         !自由剖分
AMESH,3
MAT,4        !土壤为4
Smrt,4
AMESH,4
MP,DENS,2,2740
MP,DENS,3,2300
MP,DENS,4,2050
MP,KXX,2,50
MP,KXX,3,0.15
MP,KXX,4,1.66
MP,C,2,896
MP,C,3,653
MP,C,4,1840
FINISH

/SOLU
Ic,all,temp,19
LSEL,S,,1,8    !选择排管中空气
NSLL,S,1
D,ALL,VX,0
D,ALL,VY,0             !加载速度0
ALLSEL
lsel,s,,,14,16,2
SFL,all,HFLUX,0
DL,13,,TEMP,19          !下边界加载温度19
ALLSEL
SFL,15,CONV,2.476,,39   !上边界对流换热系数2.476，温度39
ALLSEL                                                     
ASEL,S,,,1     !选择发热棒
bfa,ALL,HGEN,20000    !面上加载体积载荷，热生成率
ALLSEL

FLDATA7,PROT,DENS,AIR-SI    !空气属性国际制
FLDATA7,PROT,VISC,AIR-SI
FLDATA7,PROT,COND,AIR-SI
FLDATA7,PROT,SPHT,AIR-SI
FLDATA13,VARY,DENS,1    !密度变化  
FLDA,SOLU,TEMP,T        !激活求解温度方程
FLDA,SOLU,FLOW,F        !关闭流体方程
FLDA,RELX,TEMP,1.0      !流体物性常值，温度松弛因子1，在一个全局迭代步内即可求解温度方程，
                        !不必解流体方程，两次全局迭代导致第二次结果有小的变化。
FLDA,ITER,EXEC,1        !确定稳态分析的迭代和输出，迭代次数1
FLDA,METH,TEMP,3        !求解病态传热问题选择预条件共轭残差法（PCR）
SAVE
SOLVE
FLDATA1,SOLU,FLOW,T     !激活流体方程
FLDATA1,SOLU,VOF,0      !
FLDATA1,SOLU,TEMP,T     !激活温度方程
FLDATA1,SOLU,TURB,T     !激活紊流模型，默认关
FLDATA1,SOLU,TRAN,F     !关闭暂态分析，默认关
FLDATA1,SOLU,COMP,F     !不可压缩
FLDATA1,SOLU,IVSH,F     !关闭不可压缩流体的粘性发热选项
FLDATA1,SOLU,SWRL,F     !关闭求解轴对称旋流选项
FLDATA1,SOLU,SPEC,F     !关闭多组份疏运选项
FLDATA18,METH,PRES,1     !选取压力的TDMA方法
FLDATA19,TDMA,PRES,100   !
FLDA,METH,TEMP,6         !选取能量方程的PBCGM
FLDA,CONV,TEMP,1E-20     !控制温度的收敛监测
FLDA,RELX,TEMP,0.3       !流体物性常值，温度松弛因子0.3
FLDA,RELX,PRES,0.2       !流体物性常值，压力松弛因子0.2
FLDA,ITER,EXEC,100      !确定稳态分析的迭代和输出，迭代次数500
FLDA,ITER,OVER,50       !结果文件覆盖频率，每隔50次迭代
FLDA,OUTP,SUMF,50       !输出总结的频率，每隔50次迭代
ACEL,0,9.81,0    !线形加速度x0，y9.81，z0
SAVE
SOLVE              !稳态求解

/solu
   !改变加载，计算5秒的加热过程
asel,s,,,1
bfa,all,hgen,50000
allsel,all
FLDATA1,SOLU,TRAN,1  
FLDATA4,TIME,STEP,0.2   !时间步长
!FLDATA4,TIME,ISTEP,20  !第一步长，只在程序自定义步长时
FLDATA4,TIME,NUMB,50    !步数，缺省为10
FLDATA4,TIME,TEND,5   !暂态结束时间5秒
FLDATA4,TIME,GLOB,10    !每步迭代次数，缺省为20
FLDAta4,time,TEMP,1E-20     !控制温度的收敛监测，缺省10-6
fldata4a,step,appe,1  !输出控制，每隔1子步
SOLVE
FINISH  

/POST1
set,last
PLNSOL,TEMP

/post26
/AXLAB,X,TIME,(sec)       !指定xy轴的标签
/AXLAB,Y,TEMP
/GTHK,AXIS,3              !修改坐标线宽度、比率3为默认值的1.5倍
/GTHK,CURVE,3    !曲线宽度
/COLOR,CURVE,MRED,1   !曲线颜色
NSEL,S,LOC,X,0
NSEL,R,LOC,Y,0.02      
*GET,QA1,NODE,,NUM,MAX
NSOL,2,QA1,TEMP
plvar,2
prva,2
FINISH