请帮看看程序哪里不对（激光加热新热源项）？

dhf19822

主要是/solution后面加载部分，模型节点通过循环施加热生成率载荷（hegn），其不仅与位置有关，还与激光作用时间及计算时间有关，但是加载求解后，模型无温升变化，仍为37度，不知问题出在什么地方，载荷设置有什么问题：设的载荷步时间间隔，与子步时间一样，会不会有问题？循环求解多载荷步！请高手们指点下，急！谢谢～不胜感激～！




finish
/clear
/filnam,tissue
/units,si
/prep7
ua=70
us=1.8e4
ut=ua+us
ref=0.05 !反射率
fai=1e5  !表面功率密度
tao=1e-3 !脉冲作用时间
w=5e-4
l=1e-3
size=1e-5 ！网格尺寸
r0=3e-4  !光斑半径
et,1,plane55,,,1
mp,kxx,1,0.42 !unit
mp,c,1,3.35e3
mp,dens,1,1.09e3
rectng,0,w,-l,0
aesize,1,size
mshkey,0
mshape,0,2d
mat,1
amesh,1
finish
/solu
antype,trans
toffst,273
allsel,all
ic,all,temp,37
kbc,1
timint,on
autots,on
outres,all,all
t2=5e-3
tt=tao+t2            !j激光周期
del=2e-4                     ！载荷子步时间
*do,t,del,tao,del
time,t                          ！载荷步（载荷步间隔与子步一样，即一个载荷步一个子步）
sfl,3,,conv,10,20
deltime,del,del,del  
*do,j,0,l/2*size  ! 节点循环施加载荷（x,y,t）
y=2*size*j
*do,i,0,w/2*size
x=2*size*i
*if,x,le,r0,and,t,le,tao,then                !mod(t,tt),le,mod(tao,tt),then
bf,node(x,y,0),hgen,ua*(1-ref)*fai*exp(-0.5*x**2/r0**2*exp(-us*y)-ut*y)*exp(-4*(t-tao)**2/tao**2) ！热生成率（hegn）
*else
bf,node(x,y,0),hgen,0
*endif
*enddo
*enddo
allsel
solve
*enddo
finish

hnpsd

好像没有这样加载热源的吧,直接把数学公式放上去,ansys不支持
你可以采取表的形式加载,多看看以前的人怎么做的