急！激光表面加热温度不升，反而降！！

dhf19822

对一圆柱表面，用激光加热，建模建一截面，采用plane55轴对称，表面加一热流，边界考虑与空气的对流，对流系数为1，空气温度为25度，圆柱截面所有节点起始温度也为25度。如图，在y=0的线上加热流，即图中的上表面！温度分布如图所示，上面的温度不升，反降！大侠们慷慨些，指点下！谢谢！
finish
/clear
/filnam,tissue
/units,si
/prep7
et,1,plane55
keyopt,1,3,1
mp,kxx,1,0.2 !unit
mp,c,1,2e3
mp,dens,1,1.2e3
mp,kxx,2,0.45
mp,c,2,3.5e3
mp,dens,2,1.075e3
mp,kxx,3,0.486
mp,c,3,3.67e3
mp,dens,3,1.06e3
rectng,0,0.01,-5e-3,0
rectng,0,0.01,-0.01,-5e-3
rectng,0,0.01,-0.02,-0.01
aglue,all
numcmp,all
aesize,1,5e-4
mshkey,1
mshape,0,2d
mat,1
amesh,1
aesize,2,1e-3
mshkey,1
mshape,0,2d
mat,2
amesh,2
aesize,3,1e-3
mshkey,1
mshape,0,2d
mat,3
amesh,3
finish
/solu
h=0.3        !加热效率
r=5e-3        !高斯热源半径
pi=3.1415926
q=100              
qmax=3*q/(pi*r*r)
size1=5e-4
max_x=21
*dim,laser,table,max_x,1,1,x,y,laser1
*do,i,1,max_x,1
laser(i,0,1)=(i-1)*size1
xx=(i-1)*size1
*if,laser(i,0,1),le,r,then
laser(i,1,1)=qmax*exp(-3*xx**2/r**2)
*else
laser(i,1,1)=0
*endif
*enddo
antype,trans
toffst,273
allsel,all
ic,all,temp,25
nsel,s,node,loc,y,0,0
sf,all,hflux,%laser%
lsel,s,line,loc,x,0.01,0.01
lsel,a,line,loc,y,-0.02,-0.02
sfl,all,conv,1,25
time,5
kbc,1
timint,on
deltime,2e-1,2e-1,8e-1
autots,on
outres,all,all
solve
finish
/post1

dhf19822

大哥，大姐们，慷慨些，来指点下吧！

dhf19822

怒了，斑竹答一下吧！ :L

Strosailor

sfl,all,conv,1,25
有问题。。你仔细查查这个语句的用法。。似乎有些问题

dhf19822

好的，谢了！是有问题，应该sfl,线号，conv,1,25.情况改善不小。
但是改过之后，还是有点问题，模型初始温度25度，周围空气温度25度，加对流边界条件与激光热流密度，但是算出来之后，有些节点温度23，22度左右，竟然比初始温度25度低，不知什么原因？不是因为网格划分的有点疏吧？
不吝赐教！谢谢！

dhf19822

边界条件哪里出了问题？大侠们点点啊！antype,trans
toffst,273
allsel,all
ic,all,temp,25  ！这条命令没错吧？所有节点施加初始温度25度！
nsel,s,node,loc,y,0,0
sf,all,hflux,%laser%
lsel,s,line,loc,x,0.01,0.01
lsel,a,line,loc,y,-0.02,-0.02
sfl,all,conv,1,25
但为什么有些节点算出来温度竟然比25度低呢？

csuwumin

我下载下来看看，我现在也在搞这个激光成型。多指教啊

zwj_12

sfl加载conv的一些参数我也搞不清楚是怎么设置的？