【讨论】帮助看看这个辐射求解

mwei

试着将help的一个辐射算例4.8. Example of a 2D Radiation Analysis Using the Radiosity Method ，用辐射矩阵方法做了一下，结果对不上，百思不得其解
下面是我的命令流：
请大家帮我看看，哪的问题？？？
谢谢！
/PREP7   
CYL4,0,0,.5,0,.25,180     ! Circular annulus 1
CYL4,0.2,0,1,0,.75,180    ! Circular annulus 2
ET,1,PLANE55  
et,2,link32
et,3,matrix50,1
MP,KXX,1,.1  
MP,KXX,2,.1  
MP,KXX,3,.1  
!mp,emis,2,.9
!mp,emis,3,.7
ESIZE,0.05,   
AMESH,ALL   
!定义辐射面
LSEL,S,LINE,,1   
nsll,s,1
mat,2
type,2
esurf
LSEL,S,LINE,,7   
nsll,s,1
mat,3
type,2
esurf
n,1000,1,-10 !define a space node
  
!生成辐射矩阵
/aux12
esel,s,type,,2
nsle,s,full
geom,2d
stef,0.119E-10  
emis,2,.9
emis,3,.7
vtype,0
space,1000
write,123
alls
!使用辐射矩阵进行热分析
/prep7
type,3
se,123
esel,s,type,,2
edele,all
alls
/solu
toffst,460
d,1000,temp,70
DL,ALL, ,TEMP,1500,1    ! Temperature on inner annulus
LSEL,S,LINE,,5   
DL,ALL, ,TEMP,100,1     ! Temperature on outer annulus
ALLSEL
  
solve
fini
/post1
plnsol,temp

mwei

按上述命令流运行得到的温度分布，跟不考虑辐射得到的温度分布一样。好像我处理的辐射根本没有起作用！不解！

mwei

今天检查了一下昨天的命令流发现有个地方犯了一个弱智的错误，更正后运行结果和帮助文件的结果基本吻合，每个节点的温度误差基本都不超过0.5度。
以下是我更正后的命令流：
若有错误，敬请指出！
谢谢！
/PREP7   
CYL4,0,0,.5,0,.25,180     ! Circular annulus 1
CYL4,0.2,0,1,0,.75,180    ! Circular annulus 2
ET,1,PLANE55  
et,2,link32
et,3,matrix50,1
MP,KXX,1,.1
MP,KXX,2,.1
MP,KXX,3,.1
ESIZE,0.05,   
AMESH,ALL   
!定义辐射面
LSEL,S,LINE,,1
nsll,s,1   
type,2
mat,2  
esurf
LSEL,s,LINE,,7   
nsll,s,1
type,2
mat,3
esurf
n,1000,1,-1
  
!生成辐射矩阵
/aux12
esel,s,type,,2
nsle,s,full
geom,1
emis,2,0.9
emis,3,0.7
stef,0.119E-10  
vtype,0
space,1000
mprint,1
write,481
alls
!使用辐射矩阵进行热分析
/prep7
type,3
se,481
esel,s,type,,2
edele,all
alls
/solu
toffst,460
d,1000,temp,70
LSEL,S,LINE,,3
DL,ALL, ,TEMP,1500,1    ! Temperature on inner annulus
LSEL,S,LINE,,5   
DL,ALL, ,TEMP,100,1     ! Temperature on outer annulus
ALLSEL
  
solve
fini
/post1
plnsol,temp
ASEL,S,AREA,,1
NSLA,S,1
PRNSOL,TEMP
FINISH