哪位大哥帮我看看我的这个问题(有APDL)
小弟要做个稳态热辐射分析的课题,要求温度场的分布。只有一个热流密度载荷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左右。
出现如图错误。
1
[ 本帖最后由 huntfora 于 2007-5-14 15:58 编辑 ]
回复 #1 huntfora 的帖子
没人理我啊?有人懂吗? 温度怎么越来越高? 把热流定义成负值看一下; 另外,你不需要另外加辐射点吗?你不防换一种求解辐射的方法, 我用rdsf这样的方法计算辐射总是不成功. 可尝试在模型外另加一辐射点, 对辐射点进行网格划分,在求解中制定环境稳定(D命令).
回复 #3 weichen 的帖子
我只是找了个差不多的例子,改了下才会我的这个的。对辐射不太会。希望多讨论。
QQ:105095353,
楼上的也留下号吧,谢谢
回复 #3 weichen 的帖子
我只是找了个差不多的例子,改了下才会我的这个的。对辐射不太会。希望多讨论。
QQ:105095353,
楼上的也留下号吧,谢谢 你可以参考一下<ANSYS8.0热分析教程与实力解析>中关于辐射问题的几个例子.最后的例题是我所提到的计算辐射问题的方法.
非常可惜我很落伍,没有QQ
回复 #6 weichen 的帖子
非常谢谢!但这书我已经买了,发现写得不怎么样,还有很多错误。
但看帮助又不太懂,希望有人可以多交流交流。
回复 #7 huntfora 的帖子
有人懂吗?帮帮忙啊谢谢了
回复 #1 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 编辑 ] 自己想了好久,可能结果不对称是网格不具有周向对称性。
我改了下网格,网格是对称了。
但在模型的四周的面上(底面没加)结果还是不太理想。
下面是我的网格和结果。
回复 #10 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
讨论讨论,
谢谢! 重新建模了,解决周向不对称的问题了。
但还有瞬态和稳态的问题,不知道到底有什么区别。
我是用的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
哪位运行下我的命令啊
谢谢 下面两个图是我的结果,温度是对称了。但是用瞬态求出来的,而我只要稳态就可以了。另外,最高温度应该是2000多度,现在才200多一点。
页:
[1]