【讨论】关于辐射计算中空间节点的位置问题！

soaringbird

这个问题，我也想了好久，
做的时候也是很郁闷，刚开始想法跟你一样
  
后来我想是这样：除了对外辐射的那个表面外呢，书中并没有添加边界条件，那么意味着是绝热边条，所以，就那个表面温度低
  
你说呢？是这样吧？！
：）

mwei

这时我用radiosity求解的该问题
结果似乎可以理解
fini
/clear
/prep7   
et,1,solid70
et,2,plane55
  
mp,kxx,1,49.8   
mp,c,1,465   
mp,dens,1,7840   
mp,emis,2,1  
  
rpr4,6,0,0,0.03,,0.3
esize,0.03,
type,1
mat,1
mshkey,0
mshape,1,3d
vmesh,1
  
sfa,all,,rdsf,1,1
/psf,rdsf,emis,3
/psf,rdsf,encl,3
  
/aux12
stef,5.67e-8
toffst,273
radopt,0.1,,0
spctmp,1,20
hemiopt,10  
v2dopt,0  
vfcalc,flo,vf
  
/solu
antype,trans
autots,on  
tunif,800
kbc,1  
deltim,25   
outres,,all  
outpr,nsol,last  
time,1200   
alls
solve   
/post1
set,last
plnsol,temp
  
欢迎大家讨论，究竟那本书上出现的问题是什么原因造成的，我也对那个结果无法理解

mwei

我用ansys提供的几种方法进行了分别求解并进行对照：
1、radiosity method
fini
/clear
/prep7   
et,1,solid70
  
mp,kxx,1,49.8   
mp,c,1,465   
mp,dens,1,7840   
mp,emis,1,1  
  
rpr4,6,0,0,0.03,,0.3
esize,0.03,
type,1
mat,1
mshkey,0
mshape,1,3d
vmesh,1
  
sfa,all,,rdsf,1,1
/psf,rdsf,emis,3
/psf,rdsf,encl,3
  
/aux12
stef,5.67e-8
toffst,273
radopt,0.1,,0
spctmp,1,20
hemiopt,10  
v2dopt,0  
vfcalc,flo,vf
  
/solu
antype,trans
autots,on  
tunif,800
kbc,1  
deltim,25   
outres,,all  
outpr,nsol,last  
time,1200   
alls
solve   
/post1
set,last
plnsol,temp
  
2、radiation matrix method
fini
/clear
/prep7   
et,1,solid70
et,2,shell57
et,3,matrix50,1
  
mp,kxx,1,49.8   
mp,c,1,465   
mp,dens,1,7840   
mp,emis,2,1
  
rpr4,6,0,0,0.03,,0.3
esize,0.03,
type,1
mat,1
mshkey,0
mshape,1,3d
vmesh,1
  
type,2
mat,2
asel,all
nsla,s,1
esurf
  
n,500,0,0,0.35
  
/aux12
esel,s,type,,2
nsle,s
nsel,a,node,,500
geom,0
emis,2,1
stef,5.67e-8
vtype,0  
space,500
write,radm
alls
fini
  
/prep7
type,3
se,radm
esel,s,type,,2
edele,all
d,500,temp,20
fini
   
/solu  
alls   
antype,trans
toffst,273   
tunif,800   
autots,on   
kbc,1   
deltim,25   
outres,,all  
outpr,nsol,last  
time,1200   
solve   
  
/post1
set,last
plnsol,temp
  
3、surf152 method
fini
/clear
/prep7   
et,1,solid70
et,2,surf152,,,,1,1
  
keyopt,2,9,1
r,2,1,5.67e-8   
mp,kxx,1,49.8   
mp,c,1,465   
mp,dens,1,7840   
mp,emis,2,1  
  
rpr4,6,0,0,0.03,,0.3
esize,0.03,
type,1
mat,1
mshkey,0
mshape,1,3d
vmesh,1
  
type,2   
real,2   
mat,2   
n,500,0,0,0.35
esurf,500   
finish  
   
/solu  
alls   
antype,trans
toffst,273   
d,500,temp,20   
tunif,800   
autots,on   
kbc,1   
deltim,25   
outres,,all  
outpr,nsol,last  
time,1200   
solve   
/post1
set,last
plnsol,temp
  
三种方法的得到的温度场分布是一致的，温度场是对称的，两端低，中间高，最高温度均为310.34
那本书上的结果感觉应该是有问题的

mwei

不应该的，后面两个我运行过，可以的
你把后两个的运行结果也贴出来吧

clarino

代劳一下a，运行mwei第三种方法结果

clarino

运行mwei第二种方法结果，ai

nwpuyl

楼主很早就给我发了E－MAIl，可是由于时间问题，很抱歉昨天我才看到
mwei自己做的结果很好，但是大家对那本书的错误却没提出来，不知解决没有？
另外mwei的radoisity方法结果差异比较大哦！
我觉得市面上买的书，很难有质量比较高的，我从图书馆借了这本看了一下，笔误非常多......所以
《ansys工程应用教程热与电磁学篇》第106--112页的例子缺陷如下：
    1、建模繁琐，这是产生此例错误的根本原因－－rpr4一个命令就可以完成作者的rpr4和voffst“组合”
    2、网格划分不精确，导致对称结构对称边界模型的结果却不对称
    3、作者对实际模型和有限元模型的差异认识不清－－增加plane55只是为了方便建模，却把网格带入计算过程
修改如下即可：
/prep7   
et,1,solid70
et,2,surf22,1,,,1,1  
et,3,plane55
keyopt,2,9,1
r,2,1,5.67e-8   
mp,kxx,1,49.8   
mp,c,1,465   
mp,dens,1,7840   
mp,emis,2,1  
type,3   
rpr4,6,0,0,0.03,
esize,,3
mshkey,0          !使网格规则，结果才会对称
mshape,0,2d
amesh,1  
type,1   
esize,,10   
voffst,1,0.3
aclear,1           !清除辅助建模产生的无用网格，这是错误的关键
type,2   
real,2   
mat,2   
n,500,0,0,-0.05  
esurf,500   
finish   
/solu   
antype,trans
toffst,273   
d,500,temp,20   
tunif,800   
autots,on   
kbc,1   
deltim,25   
outres,,all  
outpr,nsol,last  
time,1200   
solve
/post1
esel,s,type,,1   !加上本句命令，才能清楚的看见温度梯度分布
plnsol,temp

huntfora

我把代码运行了一下，得出的结果怎么不和楼上的不一样啊。
是不是代码 的问题？

huntfora

有没人讨论下这个题啊？

xieyongrenren

呵呵 这样的

huntfora

为什么我把求解改成稳态时，结果却是初始温度的20度啊。
这个问题里和时间也没关系吧。

原帖由 nwpuyl 于 2004-6-21 09:45 发表
楼主很早就给我发了E－MAIl，可是由于时间问题，很抱歉昨天我才看到
mwei自己做的结果很好，但是大家对那本书的错误却没提出来，不知解决没有？
另外mwei的radoisity方法结果差异比较大哦！
我觉得市面上 ...