两边是固体,中间是流体,环境温度25°,瞬态求解出现温度低于环境温度的情况!大家帮忙看看哪里出错了,谢谢!
FINISH
/CLEAR
/FILNAME,GISHEAT
/TITLE,PHYSICAL EXPERIMENT
/UNITS,SI
/PREP7
TOFFST,273
R0=0.03
LEN0=0.3
LEN1=0.2
V=136.8 !电压
T0=12 !环境温度
Q=3E6 !发热陶瓷生热率
ET,1,FLUID141
MP,KXX,2,50.2 !外壳的导热系数
MP,DENS,2,7850 !外壳密度
MP,C,2,460 !外壳比热
MP,KXX,3,500 !外壳的导热系数
MP,C,3,400
MP,DENS,3,2700
RECTNG,0,R0,0,LEN0
RECTNG,R0,R0+LEN1,0,LEN0
RECTNG,R0+LEN1,R0+LEN1+R0,0,LEN0
AGLUE,ALL
NUMCMP,AREAS
NUMCMP,LINES
ALLSEL,ALL
ASEL,S,,,1
AATT,2,,1
ALLSEL,ALL
ASEL,S,,,3
AATT,3,,1
ASEL,S,,,2
AATT,1,,1
ALLSEL,ALL
LSEL,S,,,1,3,2
LESIZE,ALL,,,10,-5
ALLSEL,ALL
LSEL,S,,,7,8
LESIZE,ALL,,,30,-5
ALLSEL,ALL
LSEL,S,,,9,10
LESIZE,ALL,,,10,-5
ALLSEL,ALL
LSEL,S,,,2,4,2
LSEL,A,,,5,6
LESIZE,ALL,,,60
ALLSEL,ALL
MSHKEY,2
MSHAPE,0,2D
AMESH,ALL
FINISH
/SOLU
ALLSEL,ALL
ASEL,S,,,1
ESLA,S
BFE,ALL,HGEN,,Q
ALLSEL,ALL
ASEL,S,,,2
ESLA,S
NSEL,S,EXT
D,ALL,,,,,,VX,VY
ALLSEL,ALL
LSEL,S,,,6
NSLL,S,1
SF,ALL,CONV,3.59,T0 !对流换热
SF,ALL,RAD,0.85,T0 !对环境辐射
ALLSEL,ALL
LSEL,S,,,2,5,3
SFL,ALL,RDSF,0.85,,-1
/PSF,RAD,EMIS,1,1
/PSF,RAD,TAMP,1,1
/PSF,RDSF,EMSS,1,1
/PSF,RDSF,ENCL,1,1
ALLSEL,ALL
STEF,5.67E-8
RADOPT,0.5,0.0001
SPCTEMP,1,T0
HEMIOPT,10
TOFFST,273
ALLSEL,ALL
FLDATA1,SOLU,TRAN,1 !瞬态分析
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TEMP,1 !设置求解温度场
FLDATA1,SOLU,TURB,1 !设置求解湍流
FLDATA4,TIME,STEP,5
FLDATA4,TIME,NUMB,60
FLDATA4,TIME,TEND,300
FLDATA4,TIME,GLOB,50
FLDATA4,TIME,OVER,50
FLDATA4,TIME,APPE,50
FLDATA4,TIME,SUMF,50
FLDATA4,TIME,BC,1
!FLDATA4,TIME,METH,NEWM
FLDATA12,PROP,DENS,4
FLDATA13,VARY,DENS,1 !密度可变
FLDATA12,PROP,VISC,4
FLDATA13,VARY,VISC,1
FLDATA12,PROP,COND,4
FLDATA13,VARY,COND,1
FLDATA12,PROP,SPHT,4
FLDATA13,VARY,SPHT,0
FLDATA7,PROT,DENS,AIR-SI !密度
FLDATA8,NOMI,DENS,-1
FLDATA9,COF1,DENS,0
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,AIR-SI !粘性
FLDATA8,NOMI,VISC,-1
FLDATA9,COF1,VISC,0
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA7,PROT,COND,AIR-SI !导热系数
FLDATA8,NOMI,COND,-1
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,AIR-SI !比热
FLDATA8,NOMI,SPHT,-1
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
ACEL,0,9.8,0 !Y方向重力加速度
FLDATA15,PRES,REFE,101350
FLDATA16,BULK,BETA,1E15
FLDATA17,GAMM,COMP,1.4
FLDATA14,TEMP,NOMI,T0 !初始化属性时的温度
FLDATA14,TEMP,BULK,T0
FLDATA18,METH,PRES,3 !指定求解器
FLDATA18,METH,TEMP,6
ALLSEL,ALL
SOLVE
FINISH
/POST1
SET,LAST
PLNSOL,TEMP
FINISH |