线状加热温度场分布！！

huqunbo0504

我做的水火弯板的温度场模拟，和焊接的做法相似，不用单元生死，但是温度云图是这么显示的，很奇怪的形状，恳请版主及各位好手帮我分析下原因。此外，加热区内的节点温度一直保持在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

huqunbo0504

huqunbo0504

自己顶！不要沉！！！

huqunbo0504

问题已解决。

sanfengw

请教楼主解决的办法

huqunbo0504

sanfengw 发表于 2011-10-25 10:08
请教楼主解决的办法

图为热流分布，而一般我们只要查看节点的温度分布就好了。是我没搞清状况。
至于节点的温度为什么会这么低，原因是定义燃气的热值时少了三个零，应该为104600000。
问题都是概念不清和马虎所致。希望你不要像我这样。

sanfengw

huqunbo0504 发表于 2011-10-25 15:44
图为热流分布，而一般我们只要查看节点的温度分布就好了。是我没搞清状况。
至于节点的温度为什么会这么 ...

谢谢回答，通过交流让我也学到了不少

sanfengw

sanfengw 发表于 2011-10-27 09:06
谢谢回答，通过交流让我也学到了不少

请问，方便把这个问题以及完整的命令流发上来么，我想学学你这种现状加热的方法