线状加热温度场分布!!
我做的水火弯板的温度场模拟,和焊接的做法相似,不用单元生死,但是温度云图是这么显示的,很奇怪的形状,恳请版主及各位好手帮我分析下原因。此外,加热区内的节点温度一直保持在20度左右,不知是什么问题所致。附上命令流和温度场云图。……
/solu
antype,trans !定义求解类型
tunif,20 !给所有节点分配一个温度
kbc,1 !阶跃载荷
timint,on !打开瞬态效应
autots,on !自动时间步长
alls
nsel,s,loc,y,0 !在y=0处施加对称条件
dsym,symm,y
esel,s,type,,3 !表面对流的施加
nsel,s,node,,all
sf,all,conv,%conve%,20
alls
pi=3.1415926
a_C4H4=0.000143503 !C4H4的燃烧热值
q_C4H4=104600 !C4H4的流量
eta=0.325 !热效率
r0=38e-3 !热源有效半径m
v=4e-3 !热源移动速度m/s
qmax=3*a_C4H4*q_C4H4*eta/(r0*r0*pi) !热流密度最大值
tm=len1/v !总加热时间
tinc=esize1/v !单步时间长度
allsel
asel,s,,,2 !选择加热区的上表面152单元
nsla,s,0
esln,s,0
esel,u,type,,3
cm,surf1,elem
*get,nemax,elem,,num,max
*get,nemin,elem,,num,min
max_time=tm/tinc !循环步数
*do,j,1,max_time,1
cmsel,s,surf1
*do,i,nemin,nemax,1
*if,esel(i),eq,1,then !如果该单元在上面的集合里,则
xsy=centrx(i) !读取该单元的中心坐标
ysy=centry(i)
zsy=centrz(i)
a=v*tinc*j
b=0
c=thic
di=abs(sqrt((xsy-a)*(xsy-a)+(ysy-b)*(ysy-b)+(zsy-c)*(zsy-c)))
!该单元中心距离热源中心
*if,di,le,r0,then !该单元中心在加热半径范围内
qr=qmax*EXP(-3*di**2/(r0**2))!该单元中心处的热流大小
sfe,i,6,hflux,,qr !在每个单元所指定的面上施加热源载荷
*endif
*endif
*enddo
allsel,all
time,j*tinc !载荷步结束的时间值
nsubst,1 !子步数和时间步长
solve
sfedele,all,all,hflux
*enddo !全选求解
save
finish
E\:wenduchang.png http://ww4.sinaimg.cn/large/654e2f0ejw1dmebzo2ts6j.jpg 自己顶!不要沉!!! 问题已解决。 请教楼主解决的办法 sanfengw 发表于 2011-10-25 10:08 static/image/common/back.gif
请教楼主解决的办法
图为热流分布,而一般我们只要查看节点的温度分布就好了。是我没搞清状况。
至于节点的温度为什么会这么低,原因是定义燃气的热值时少了三个零,应该为104600000。
问题都是概念不清和马虎所致。希望你不要像我这样。 huqunbo0504 发表于 2011-10-25 15:44 static/image/common/back.gif
图为热流分布,而一般我们只要查看节点的温度分布就好了。是我没搞清状况。
至于节点的温度为什么会这么 ...
谢谢回答,通过交流让我也学到了不少 sanfengw 发表于 2011-10-27 09:06 static/image/common/back.gif
谢谢回答,通过交流让我也学到了不少
请问,方便把这个问题以及完整的命令流发上来么,我想学学你这种现状加热的方法
页:
[1]