哪位大哥帮我看看我的这个问题（有APDL）

huntfora

小弟要做个稳态热辐射分析的课题，要求温度场的分布。
只有一个热流密度载荷2.5MW/M2。密度6000kg/m3比热610J/(kg*k)。
模型尺寸总长为100mm,圆柱半径为37.5mm.
发现设为瞬态和稳态时结果相差很大，不知为什么。哪位帮忙看看，谢谢。
命令如下：
fini
/clear
/PREP7
!定义单位制
/UNITS,SI

!modeling
bom_r=0.0375                   !导弹半径
bom_l=0.0635                    !导弹长度
CYL4, , ,bom_r, , , ,bom_l

wpro,,-30.000000,   
wpro,,-30.000000,   
wpro,,-30.000000,   

SPHERE,bom_r, ,0,180,

!select element
et,1,solid70

!material properity
mp,kxx,1,60
mp,c,1,610
mp,dens,1,6000
mp,emis,1,1  

!mesh
esize,,5
mshkey,0
mshape,1,3d
vmesh,all

sfa,all,,rdsf,1,1
/psf,rdsf,emis,3
/psf,rdsf,encl,3

/aux12
stef,5.67e-8
toffst,273
radopt,0.1,0.0001,0         
spctmp,1,20
hemiopt,10  
v2dopt,0, , 0,200
vfcalc,hemic,vf   

/solu
antype,trans
!当设为稳态时，结果和瞬态相差很大。
autots,on
aSEL,s,area,,3,6
SFA,all,1,HFLUX,2.5e6
allsel
tunif,20     !uniform temp
kbc,1        !stepped loading
deltim,0.5
outres,all    !Controls the solution data written to the database.
outpr,all,last  
time,1
alls
solve   
/post1
set,last
plnsol,temp

瞬态分析的结果图片太大，上传不了。当设为稳态时，出现错误，实际最高温度应该为2600左右。
出现如图错误。
[local]1[/local]

[ 本帖最后由 huntfora 于 2007-5-14 15:58 编辑 ]

huntfora

没人理我啊？
有人懂吗？

weichen

温度怎么越来越高? 把热流定义成负值看一下; 另外,你不需要另外加辐射点吗?你不防换一种求解辐射的方法, 我用rdsf这样的方法计算辐射总是不成功. 可尝试在模型外另加一辐射点, 对辐射点进行网格划分,在求解中制定环境稳定(D命令).

huntfora

我只是找了个差不多的例子，改了下才会我的这个的。
对辐射不太会。希望多讨论。
QQ：105095353，
楼上的也留下号吧，谢谢

huntfora

我只是找了个差不多的例子，改了下才会我的这个的。
对辐射不太会。希望多讨论。
QQ：105095353，
楼上的也留下号吧，谢谢

weichen

你可以参考一下<ANSYS8.0热分析教程与实力解析>中关于辐射问题的几个例子.最后的例题是我所提到的计算辐射问题的方法.
非常可惜我很落伍,没有QQ

huntfora

非常谢谢！
但这书我已经买了，发现写得不怎么样，还有很多错误。
但看帮助又不太懂，希望有人可以多交流交流。

huntfora

有人懂吗？帮帮忙啊
谢谢了

huntfora

fini
/clear
keyw,pr_therm,1
/prep7   
et,1,solid70
et,2,surf22,1,,,1,1  
keyopt,2,9,1
r,2,1,5.67e-8   
mp,kxx,1,60
mp,c,1,610   
mp,dens,1,6000
mp,emis,2,0.8
type,1  

!modeling
bom_r=0.0375              
bom_l=0.0635                  
CYL4, , ,bom_r, , , ,bom_l

wpro,,-30.000000,   
wpro,,-30.000000,   
wpro,,-30.000000,   

SPHERE,bom_r, ,0,180,
vadd,all

esize,,5
mshkey,0         
mshape,1,3d
vmesh,all

vclear,1           
type,2   
real,2   
mat,2   
n,1000,0,0,-0.05  
esurf,1000   
finish   

/solu   
antype,trans
！这里还是瞬态，我的要求是稳态的
csys,1
aSEL,s,area,,3,10
SFA,all,1,HFLUX,2.5e6       ！这里怎么样加成周向对称的载荷？
allsel
toffst,273   
d,1000,temp,20   
tunif,20  
autots,on   
kbc,1   
deltim,0.5
outres,,all  
outpr,nsol,last  
time,1
solve
/post1
esel,s,type,,1   
plnsol,temp

这是另一种方法的命令，哪位运行下我的命令，谢谢了
改正了些错误
还有两个错误：
1、热流密度应该是同向对称的，故温度场也应是周向对称的
2、现在上面的设置是瞬态，但我的结果要求是稳态的。
下面两张图片是我的结果。希望有人多讨论讨论，
谢谢


[ 本帖最后由 huntfora 于 2007-5-19 10:39 编辑 ]

huntfora

自己想了好久，可能结果不对称是网格不具有周向对称性。
我改了下网格，网格是对称了。
但在模型的四周的面上（底面没加）结果还是不太理想。
下面是我的网格和结果。

huntfora

下面是命令：
fini
/clear
/PREP7  
et,1,solid70
et,2,surf22,1,,,1,1  
keyopt,2,9,1
r,2,1,5.67e-8   
mp,kxx,1,60
mp,c,1,610   
mp,dens,1,6000
mp,emis,2,0.8
type,1  

!*
*SET,bom_r,0.0375   
*SET,bom_l,0.0635   
RECTNG,0,bom_r,0,-bom_l,
PCIRC, ,bom_r,0,90,
K,20,bom_r/5,-0.01-bom_l,,
K,21,bom_r/5,0.01+bom_r,,
LSTR,20,21
ASBL,ALL,8
NUMMRG,KP, , , ,LOW
LESIZE,2, , ,10, , , , ,1
LESIZE,14, , ,10, , , , ,1
LESIZE,10, , ,8,0.2, , , ,1
LESIZE,12, , ,8,0.2, , , ,1
LESIZE,11, , ,8, , , , ,1
LESIZE,9, , ,8,0.2, , , ,1
AMESH,4
AMESH,3
ALLSEL,ALL  
TYPE,   1   
EXTOPT,ESIZE,6,0,   
EXTOPT,ACLEAR,1
!*  
EXTOPT,ATTR,1,1,0   
ESYS,0  
!*
FLST,2,4,5,ORDE,2   
FITEM,2,3   
FITEM,2,-6  
FLST,8,2,3  
FITEM,8,1   
FITEM,8,6   
VROTAT,P51X, , , , , ,P51X, ,360,4,  
csys,1
VSEL,S,LOC,X,0,bom_r/5

VSEL,S, , ,       3
VSEL,A, , ,       4
VSEL,A, , ,       7
VSEL,A, , ,       8
VSEL,A, , ,       11
VSEL,A, , ,       12
VSEL,A, , ,       15
VSEL,A, , ,       16
VMESH,ALL

NUMMRG,NODE
NUMMRG,KP
NUMMRG,ELEMENT

type,2   
real,2   
mat,2   
n,10000,0,0,-0.05  
esurf,10000
finish   

/solu   
antype,trans
aSEL,s,area,,1,18
aSEL,a,area,,16,29
aSEL,a,area,,42,52
aSEL,a,area,,31,44
aSEL,a,area,,10,48
aSEL,a,area,,23,36
SFA,all,1,HFLUX,2.5e6
allsel
toffst,273   
d,1000,temp,20   
tunif,20  
autots,on   
kbc,1   
deltim,0.5
outres,,all  
outpr,nsol,last  
time,1
solve
/post1
esel,s,type,,1   !加上本句命令，才能清楚的看见温度梯度分布
plnsol,temp
讨论讨论，
谢谢！

huntfora

重新建模了，解决周向不对称的问题了。
但还有瞬态和稳态的问题，不知道到底有什么区别。
我是用的vrotat来生成对称网格的。
命令好下：
fini
/clear
keyw,pr_therm,1
/prep7   
et,1,solid70
et,2,surf152
keyopt,2,5,1
keyopt,2,9,1
r,2,1,5.67e-8   
mp,kxx,1,60
mp,c,1,610   
mp,dens,1,6000
mp,emis,2,0.8
type,1  

!modeling
m_r=0.0375                   !导弹半径
m_l=0.0635                    !导弹长度
recting,0,m_r,0,m_l
CYL4,0,m_l, , ,m_r,90
asel,s,,,1,2,1
aadd,all


lsel,s,,,2,4,2
lesize,all,,,10
lsel,s,,,1,9,8
lesize,all,,,5
!lsel,s,,,6,7,1
!lesize,all,,,5
lsel,s,,,8
lesize,all,,,8
MSHKEY,0          !使网格规则，结果才会对称
mshape,0,2d
amesh,all          !面划网格

ALLSEL
type,1
esize,,5
vrotat,all,,,,,,1,6,360
aclear,all

type,2   
real,2   
mat,2   
n,10000,0,0,-0.05  
esurf,10000
finish   

/solu   
antype,trans
csys,1
aSEL,s,area,,1,4,3
aSEL,a,area,,6,8,2
aSEL,a,area,,10,12,2
aSEL,a,area,,14,16,2
SFA,all,1,HFLUX,2.5e6
allsel
toffst,273   
d,10000,temp,20   
tunif,20  
autots,on   
kbc,1   
deltim,0.5
outres,,all  
outpr,nsol,last  
time,1
solve
/post1
esel,s,type,,1   !加上本句命令，才能清楚的看见温度梯度分布
plnsol,temp
哪位运行下我的命令啊
谢谢

huntfora

下面两个图是我的结果，温度是对称了。但是用瞬态求出来的，而我只要稳态就可以了。另外，最高温度应该是2000多度，现在才200多一点。