求助：生死单元后，再次求解不合理

xiebao143

激光去除加工过程仿真，在第一次加载后，设定超过熔化部分的全去除（杀死或者改变材料属性），自然冷却，在加载，求解，去除，冷却。重启动重复这个过程，但是在冷却时结果不合理，也就是生死单元后，就出错了，无法冷却求解，期待您的帮助。程序如下：

/filename,killelment

/title, the simulate of laser cladding process
/prep7
BLOCK,0,100e-6,0,-50e-6,-50e-6,50e-6,

ET,1,SOLID70

ET,2,SURF152

keyopt,2,4,1

keyopt,2,5,0

keyopt,2,8,1  !hflux

toffst,273

mp,DENS,1,7840

mp,kxx,1,34

mp,c,1,500

Mshkey,1
Mshape,0,3d
type,1

Mat,1
esize,5e-6,
Vmesh,all
Eplot
asel,all
type,2

Mat,1
amesh,all
esel,s,type,,1

Asel,s,,,1,6

Nsla,s,1
sf,all,conv,800,25
allsel,all
/psf,hflux,,2,0

/SOLU  
antype,transient,new
TRNOPT,FULL
Kbc,1
tunif,25

neqit,25

!nlgeom,on
autots,on
timint,on,therm
TINTP,,,,1
Lnsrch,on

nropt,full
outres,all,all
Allsel
nsel,all

Aaa=50e-6
local,11,1, Aaa,0,0,0,-90,0
CSYS,11
esel,s,type,,2
nsel,r,loc,x, 0,10e-6
nsel,r,loc,y,0,180
nsel,r,loc,z,0
sf,,hflux,7e10

t=0
dt=5e-6
csys,0
t=t+dt
time,t
deltim,1e-6,0.5e-6,1e-6
Allsel
Solve

save

Finish

/post1

Set,last

etable,melty22,temp,

ESEL,S,ETAB,MELTY22,400

finish

/solu

*get,elemnum,elem,0,count
   *if,elemnum,gt,0,then
   ekill,all

esel,s,live

   nsle,s,1

   nsel,invert

   d,all,all,0
   *else

allsel
   *endif

Esel,s,live
!在这之前都没问题，求解后结果就出错了

sfdele,all,hflux！删除载荷，冷却

t=t+2*dt
time,t

Allsel
Solve

finish

!后面的就更无法运算

/solu

   resume

   antype,,rest！重复上述过程