yangyuebiao 发表于 2009-4-7 14:19:21

【求助】热辐射问题！（帮我看看APDL程序）

本帖最后由 yangyuebiao 于 2009-4-7 18:12 编辑

下面是我的命令流，但是计算出来的结果仍然是模型的原温度，不知道是什么原因，希望大家帮忙看看，谢谢了！

!****************************前处理部分************************************
/prep7
et,1,plane55                           !定义单元
et,2,surf151
KEYOPT,2,4,1
KEYOPT,2,5,1
KEYOPT,2,6,0
KEYOPT,2,7,1
KEYOPT,2,8,2
KEYOPT,2,9,1
mptemp,1,30,1000,1200,1509.6,1528.5,1600 !定义温度         
mpdata,dens,1,1,7800,7405.93,7324.68,7200,7061,7061
                                       !定义密度                           
mpdata,kxx,1,1,52,27,30,31.333,420,420   !定义钢的导热系数
mpdata,enth,1,1,158e6,5226.38e6,6271.41e6,7875e6,9320e6,9880e6
                                       !定义钢的热焓
MP,EMIS,1,0.85
MP,EMIS,2,1                                       
!*************几何模型建立及网格划分***************************************
cyl4,0.215,0.170,0,0,0.010,90,0          ! 定义圆角
rectng,0,0.215,0,0.170   
rectng,0,0.215,0.170,0.180
rectng,0.215,0.225,0,0.170
aadd,all
numcmp,area
numcmp,line
lcomb,2,6
lcomb,3,7
numcmp,line
lesize,5,,,43
lesize,2,,,45
lesize,3,,,34
lesize,4,,,36
lesize,1,,,4
lsel,all
mat,1
amesh,1
TYPE,2
REAL,2
MAT,2
*GET,NN,NODE,,num,max
LSEL,S,,,5
NSLL,S,1
LSEL,a,,,1
NSLL,a,1
LSEL,a,,,4
NSLL,a,1
N,NN+1,0.3,0.375,0
ESURF,NN+1
ALLSEL                                       
SAVE
!**************************************************************************
!**************************求解部分****************************************
/solu
antype,trans
TIME,500
AUTOTS,1
DELTIM,5, , ,1
KBC,0
ALLSEL
CNVTOL,TEMP,,0.05,2
bfunif,temp,1548.5
STEF,5.67E-8
TOFFST,273
ALLSEL
d,NN+1,temp,27
KBC,1
NSUBST,2
OUTRES,all,ALL
outpr,nsol,last
SOLVE
finish
http://home.simwe.com/attachment/200904/7/377883_1239097300f647.jpg

该平面以一定的初始温度置于空气中，其中2和3两条边是绝热的，其余的边与空气发生辐射传热，我想求解一段时间后该平面的温度场！

yangyuebiao 发表于 2009-4-7 16:15:02

我修改了一下，可是程序还是有问题啊，大家帮帮忙？

!****************************前处理部分************************************
/prep7
et,1,plane55                           !定义单元
et,2,surf151
KEYOPT,2,4,1
KEYOPT,2,5,1
KEYOPT,2,6,0
KEYOPT,2,7,1
KEYOPT,2,8,2
KEYOPT,2,9,1
mptemp,1,30,1000,1200,1509.6,1528.5,1600 !定义温度         
mpdata,dens,1,1,7800,7405.93,7324.68,7200,7061,7061
                                       !定义密度                           
mpdata,kxx,1,1,52,27,30,31.333,420,420   !定义钢的导热系数
mpdata,enth,1,1,158e6,5226.38e6,6271.41e6,7875e6,9320e6,9880e6
                                       !定义钢的热焓
MP,EMIS,1,0.85
MP,EMIS,2,1
                                       
!*************几何模型建立及网格划分***************************************
cyl4,0.215,0.170,0,0,0.010,90,0          ! 定义圆角
rectng,0,0.215,0,0.170   
rectng,0,0.215,0.170,0.180
rectng,0.215,0.225,0,0.170
aadd,all
numcmp,area
numcmp,line
lcomb,2,6
lcomb,3,7
numcmp,line
lesize,5,,,43
lesize,2,,,45
lesize,3,,,34
lesize,4,,,36
lesize,1,,,4
lsel,all
mat,1
amesh,1
TYPE,2
REAL,2
MAT,2
*GET,NN,NODE,,num,max
LSEL,S,,,5
NSLL,S,1
LSEL,a,,,1
NSLL,a,1
LSEL,a,,,4
NSLL,a,1
N,NN+1,0.3,0.375,0
ESURF,NN+1
ALLSEL                                       
SAVE
!**************************************************************************

!**************************************************************************
/AUX12
esel,s,type,,2
nsle,s,all
nsel,a,node,,NN+1
geom,1
EMIS,1,0.85
EMIS,2,1
STEF,5.67E-8
vtype,1
N,NN+1,0.28,0.0.35,0
space,NN+2
write,casting
ALLSEL
FINISH

/prep7
DOF,TEMP
et,3,matrix50
type,3
se,casting
ALLSEL
FINISH
!**************************求解部分****************************************
/SOLU
antype,trans
TIME,5000
AUTOTS,1
DELTIM,1, , ,1
KBC,0
ALLSEL

CNVTOL,TEMP,,0.05,2

bfunif,temp,100
TOFFST,273
ALLSEL
d,NN+1,temp,27
KBC,1
NSUBST,2
OUTRES,all,ALL
outpr,nsol,last
SOLVE
finish

yangyuebiao 发表于 2009-4-7 17:47:37

本帖最后由 yangyuebiao 于 2009-4-7 18:21 编辑

似乎是个很简单的辐射问题，可以就是做不出来啊
看来还自己不会的原因
希望大家多多指点！

yangyuebiao 发表于 2009-4-7 18:14:09

RE: 【求助】热辐射问题！（帮我看看APDL程序）

http://home.simwe.com/attachment/200904/7/377883_12390973072ksn.jpg
网格划分

yangyuebiao 发表于 2009-4-7 19:24:32

:'(怎么没有回帖啊
麻烦版主帮帮忙，急啊

yangyuebiao 发表于 2009-4-8 07:48:35

大家帮帮忙啊
页: [1]
查看完整版本: 【求助】热辐射问题！（帮我看看APDL程序）