崩溃ing，谁能告诉我，简单的热处理模拟到底哪里出错了

samlin123

大家好，我在做一个简单的热处理模拟，底面为铝基板，尺寸为半径13mm，厚度1.63mm，铝基板上一个LED，LED尺寸为半径2mm，厚度2.55mm，LED上加载1W热源，周围空气对流系数10W/M2，我做了以下ansys分析：
finish
/clear
/filname,LED module
/title,LED module
/prep7
keyw,pr_thermal,1
et,1,solid70
mp,kxx,1,159!定义铝导热系数,单位 W/M.C
mp,dens,1,2700!定义铝密度,单位 KG/M3
mp,c,1,963!定义铝比热,单位 J/Kg.C
cylind,0.013,0,0,0.00163!定义基板尺寸
cylind,0.004,0,0,0.00418!v=32!定义热源尺寸
block,-0.001,0.001,-0.001,0.001,0,0.00418!分割热源，以便网格划分
vovlap,all
vglue,all
numcmp,all
/pnum,volu,1
/pnum,area,1
vplot
vplot
allsel
esize,0.0002
vsweep,1
vsweep,2
allsel
esize,0.0004
vsweep,4
vsweep,5
allsel
esize,0.002
vsweep,3
finish
/solu
antype,4
TRNOPT,FULL
LUMPM,0
NSUBST,5,8,3
OUTRES,ERASE
OUTRES,ALL,ALL  
AUTOTS,1
TIME,5000  
allsel
TUNIF,18,
allsel
asel,s,loc,z,0!选择底面
nsla,s,1
sf,all,CONV,10,18!应用底面，对流系数10W/M2
allsel
asel,s,loc,z,0.00163
asel,a,,,1,2,1
asel,u,,,3
asel,u,,,19
nsla,s,1!选择正面及侧面
sf,all,CONV,10,18!对流系数10W/M2
allsel
vsel,s,,,2
vsel,s,,,5
nslv,s,1!选择体2及体5，体2＋体5为发热源
bf,all,hgen,2.6e7!应用发热率，发热量为1W，发热率我以为1W/（小圆柱体体积）×85％计
allsel
solve
finish
/post1
plnsol,temp
finish




计算结果竟然打到300度，明显不对啊，实际温度正常只有50度左右，为什么误差这么大，请各位指点下，不胜感激
偶用mm，mg为长度及重量单位，计算结果也是一样，
命令如下：
finish
/clear
/filname,LED module
/title,LED module
/prep7
keyw,pr_thermal,1
et,1,solid70
mp,kxx,1,0.159!定义铝导热系数,单位 W/mm.C
mp,dens,1,2.700!定义铝密度,单位 mg/mm3
mp,c,1,0.000963!定义铝比热,单位 J/mg.C
cylind,13,0,0,1.63!定义基板尺寸
cylind,4,0,0,4.18!v=32定义热源尺寸
block,-1,1,-1,1,0,4.18!分割热源，以便网格划分
vovlap,all
vglue,all
numcmp,all
/pnum,volu,1
/pnum,area,1
vplot
vplot
allsel
esize,0.2
vsweep,1
vsweep,2
allsel
esize,0.4
vsweep,4
vsweep,5
allsel
esize,2
vsweep,3
finish
/solu
antype,4
TRNOPT,FULL
LUMPM,0
NSUBST,5,8,3
OUTRES,ERASE
OUTRES,ALL,ALL  
AUTOTS,1
TIME,5000  
allsel
TUNIF,18,
allsel
asel,s,loc,z,0!选择底面
nsla,s,1
sf,all,CONV,0.000010,18!应用底面，对流系数0.000010W/mm2
allsel
asel,s,loc,z,1.63
asel,a,,,1,2,1
asel,u,,,3
asel,u,,,19
nsla,s,1!选择正面及侧面
sf,all,CONV,0.00001,18!应用底面，对流系数0.000010W/mm2
allsel
vsel,s,,,2
vsel,s,,,5
nslv,s,1!选择体2及体5，体2＋体5为发热源
bf,all,hgen,0.026!应用发热率，发热量为1W，发热率我以为1W/（小圆柱体体积）×85％计
allsel
solve
finish
/post1
plnsol,temp
finish





不明白啊，不明白，请高手来解答下

[ 本帖最后由 samlin123 于 2009-1-13 20:10 编辑 ]

samlin123

高手来解答一下啊，真郁闷啊