请高人帮忙看看热分析问题出在哪？
我要做的是一个厚为5mm的薄壁筒包在圆柱型金属材料上，在薄壁筒上施加功率20kw的热源对圆柱型金属进行加热，圆柱型金属材料参数随温度变化，空气对流系数取10w/m2，求半小时后金属材料的温度分布和温度随时间的变化曲线图.我将问题简化一维来处理，ansys算到50-60秒就死机没有响应了！求各位高人帮忙看看问题出在哪？命令流：
finish
/clear
/filenam,electric heat
/title,electric heat of a solid cylinder barrel
keyw,pr_therm,1
/prep7
smrt,off
riw=.03 ！inner of solid cylinder barrel
row=.065 ！outer of solid cylinder barrel
roc=.07 ！outer of solid wall
t=.001 ！model thickness
Hair=10
Tair=20
et,1,55,,,1
et,2,55,,,1
et,3,55,,,1
emunit,mks
mp,dens,1,7.85e-3
mptemp
mptemp,1,20,100,200,300,400,500
mptemp,7,600,700,800,900,1000
mpdata,kxx,1,1,47.68,43.53,40.44,38.13,36.02,34.16
mpdata,kxx,1,7,31.98,28.66,26.49,25.92,24.02
mptemp
mptemp,1,20,100,200,300,400,500
mptemp,7,600,700,720,742,760
mpdata,c,1,1,472,480,498,524,560,615
mpdata,c,1,7,700,854,916,986,1064
/pnum,area,1
/pnum,line,1
/pnum,kp,1
rectng,0,riw,0,t
rectng,riw,row,0,t
rectng,row,roc,0,t
aglue,all
numcmp,area
asel,s,area,,1
aatt,1,1,1
asel,s,area,,2
aatt,1,1,2
asel,s,area,,3
aatt,1,1,3
asel,all
lsel,s,loc,x,riw/2
lesize,all,,,15
lsel,s,loc,y,t/2
lesize,all,,,1
lsel,all
mshape,0,2d
mshkey,1
amesh,1
lsel,s,loc,y,0
lsel,a,loc,y,t
lsel,u,loc,x,riw/2
lesize,all,.001
lsel,all
mshape,0,2d
mshkey,1
amesh,all
/rgb,index,100,100,100,0
/rgb,index,80,80,80,13
/rgb,index,60,60,60,14
/rgb,index,0,0,0,15
/replot
eplot
/solu
antype,trans
timint,on
autots,on
outres,,all
time,1800
!nsubst,360
deltim,5
toffst,273
!bfunif,temp,20
allsel
tunif,20
allsel
kbc,1
lsel,s,,,10
nsll,s,1
sf,all,conv,Hair,Tair
esel,s,type,,3
nsle,s,1
bf,all,hgen,3.932e7
allsel
eqslv,jcg
solve
finish
/post26
/axlab,x,TIME,(sec)
/axlab,y,TEMPERATURE
/gthk,axis,3
/gthk,curve,3
/replot
/xrange,0,1800
! /yrange,20,
nsol,2,1,temp,,tempcl
nsol,3,2,temp,,tempsurf
nsol,4,6,temp,,tempi
plvar,2,3,4
finish 没人回？哭了！
页:
[1]