『求助』流体分析设置问题
本帖最后由 xlpe 于 2009-5-10 18:01 编辑我用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
页:
[1]