热分析

高手们请看看这个贴，我在做激光表面淬火，下面是操作命令流，用的是ansys14，为什么总是出现

这个提示，然后ansys就不响应了。
求高手解释！
/prep7 !进入前处理
/units,si !设置单位制为国际单位制
*set,x0,0 !赋值，类似于C++中的“＝”，下同
*set,y0,0
*set,rb,0.002 !定义激光光斑半径
*set,Q,8000 !定义激光的功率密度
*set,PI,3.1415926
*set,v,0.05 !定义激光的移动速度
*set,tim,0
QC=Q/(PI* rb*rb)
et,1,solid70 ! 定义单元类型
mp,dens,1,7833 !定义密度
mp,c,1,465 !定义比热
mp,kxx,1,54 !定义热传导系数
block,0,0.052,0,0.025,-0.025,0 !建长方体模型
/view,1,1,1,1 !转变视图方向为ISO视角
esize,0.002 !单元大小
mshkey,0 !网格设置
mshape,1,3D
vmesh,all !划分网格
nsel,s,loc,x,0,0.05 !选择激光光斑区域
nsel,r,loc,y,0,0.005
nsel,r,loc,z,-0.01,0
nrefine,all,3 !对激光光斑作用的区域进行网格细化
save
x_center=x0+rb !设置光斑的初始位置
y_center=y0
然后通过gui方式调用函数编辑器定义要加载的热流密度函数,函数的形式为
QC*exp(-2*(({X}-x_center)^2+({Y}-y_center)^2)/(rb^2)),保存后退出,利用函数加载器
加载这个函数,在弹出的对话框中定义Table Parameter Name为
laser,然后将函数式中的参数分别定义为QC=QC,x_center=x_center,
y_center=y_center,rb=rb
光斑的移动采用循环编制,每一次移动一个光斑的位置.在加载区域的选择上,由于圆形区域
不好选择,所以采用圆形区域的外接正方形区域来近似.
再输入第二部分命令流
/sol
antype,4,new !定义分析类型
trnopt,full
toffst,0 !设置采用的温度单位为K
tunif,300 !设置初始温度为300K
!---------------------------
*do,i,1,13
tim=tim+0.04
time,tim !设置载荷步作用时间
x_center=x0+(i-1)*2*rb !设置当前作用区域的中心点坐标 (每次x向前移动一个直径的距离）
y_center=y0
*set,laser(4,0,1),x_center !改变定义的函数中响应的中心点坐标值
nsel,s,loc,x,x_center-0.002,x_center+0.002
nsel,r,loc,y,y_center,y_center+0.002
nsel,r,loc,z,0
sf,all,hflux,%laser% !在选择的节点上面施加定义好的热流密度函数
allsel,all
nsubst,5
solve
*enddo