仿真3D打印過程溫度變化，但一直出問題，請各位大神幫忙看一下
我這是仿真的3D打印過程溫度變化的命令流。我是現(xiàn)將所有單元都?xì)⑺溃缓笠粋€一個的激活，激活一個單元后，給這個單元加載進行分析，分析完之后將該單元 卸載，然后給下一單元加載再進行分析，一直重復(fù)此過程，直到所有單元分析完成。但在生死單元循環(huán)中，即*DO循環(huán)中。那個RESTART命令提示必須加一 個FINISH命令才能開始RESTART,但是加了之后，仍然不對。請各位幫忙看一下我的命令流哪里有問題，謝謝啦
/BATCH
FINISH
/CLEAR,START
/CWD,G:\Exercise
/FILNAME,FDMTEMPERTURE,0
/TITLE,ANSYS OF FDM TEMPERTURE
*SET,L1,1E-2
*SET,W1,1.2E-3
*SET,H1,0.3E-3 !設(shè)置模型尺寸為10mm*1.2mm*0.3mm
!設(shè)置網(wǎng)格尺寸
*SET,L2,L1/25
*SET,W2,0.4E-3
*SET,H2,H1/3 !設(shè)置單元尺寸為0.4mm*0.4mm*0.1mm（噴嘴直徑為0.4，分層厚度0.1mm）
ARG1=L1/L2*W1/W2*H1/H2 !計算單元總數(shù)量
/CONFIG,NRES,10000
!*************************************
!進入前處理
!*************************************
/PREP7
ET,1,SOLID70
MP,DENS,1,1150 !設(shè)置密度
MPTEMP,1,50,100,150,200 !指定四個溫度
MPDATA,C,1,1,1470,1490,1710,2020 !比熱容
MPDATA,KXX,1,1,3.0E-2,2.8E-2,2.9E-2,3.3E-2 !導(dǎo)熱系數(shù)
MPDATA,EX,1,1,3.5E5,2.48E9,1.68E9,0.5E9 !彈性模量
MPDATA,NUXY,1,1,0.38,0.39,0.40,0.41 !泊松比
MPDATA,ALPX,1,1,8.51E-5,8.42E-5,8.40E-5,8.38E-5 !熱膨脹系數(shù)
BLOCK,0,L1,0,W1,0,H1 !建立成型件的20mm*1mm*0.3mm的模型
!網(wǎng)格劃分
/PNUM,AREA,1 !顯示面編號
VATT,1,,0 !體單元類型solid70
/PNUM,LINE,1
MSHAPE,0,3D !單元形狀六面體
MSHKEY,1 !映射網(wǎng)格劃分
LESIZE,1,W2
LESIZE,4,L2
LESIZE,9,H2 !單元尺寸大小
VMESH,1
EPLOT
/VIEW,,-1,-2,-3
/IMAGE,SAVE,VMESH,BMP
FINISH
!**************************************
!進入載荷加載和求解
!*************************************
/SOLU
ANTYPE,TRANSIENT,NEW
AR21=0 !每一個載荷步加載時間是0.02s，單元激活時間間隔也是0.02s，掃描速度20mm/s
AR22=0.02
TIME,AR21
AUTOTS,ON !自動時間步選項（非線性分析打開加速收斂）
NROPT,FULL,,ON ! 牛頓-拉普森選項
TRNOPT,FULL !瞬態(tài)分析選項，指定最大分析模態(tài)數(shù)
LUMPM,OFF !使用一致質(zhì)量矩陣
KBC,1
DELTIM,AR21/4,AR21/4,AR21/2,0 !載荷子步設(shè)定
TIMINT,ON
MAX_TEMP=220 !設(shè)定溫度值220攝氏度
NSEL,S,LOC,X,0,L1,L2
SF,ALL,CONV,72,25
NSEL,S,LOC,Y,0,W1,W2
SF,ALL,CONV,72,25
NSEL,S,LOC,Z,0,H1,H2
SF,ALL,CONV,72,25 !初始化所有節(jié)點，加載對流
/REPLOT,FAST
OUTRES,ALL,ALL
OUTPR,NSOL,ALL
!生死單元循環(huán)命令
ESEL,S,ELEM,,1,ARG1,1,0
EKILL,ALL
ALLSEL,ALL
*DO,AR25,1,H1/H2,1
*DO,AR23,1,L1/L2,1
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0 !依次選擇下一個單元激活
EALIVE,ALL
*DO,J,1,8,1
D,NELEM((AR25-1)*W1/W2*L1/L2+AR23,J),TEMP,MAX_TEMP !給每個單元的八個節(jié)點施加恒定溫度
*ENDDO
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0
NSLE,S,ALL !選擇單元中的所有節(jié)點
SF,ALL,CONV,72,25
ALLSEL,ALL
SOLVE
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0
*DO,J,1,8,1
DDELE,NELEM((AR25-1)*W1/W2*L1/L2+AR23,J),TEMP !刪除每個單元的八個節(jié)點施加恒定溫度
*ENDDO
ALLSEL,ALL
AR22=AR22+AR21
!*IF,AR22,GT,AR21*ARG1,*EXIT
FINISH
ANTYPE,TRANSIENT,RESTART
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0
*DO,J,1,8,1
DDELE,NELEM((AR25-1)*W1/W2*L1/L2+AR23,J),TEMP !刪除每個單元的八個節(jié)點施加恒定溫度
*ENDDO
TIME,AR22
*ENDDO
*DO,AR26,L1/L2,1,-1
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+L1/L2+AR26,,,0
EALIVE,ALL
*ENDDO
*DO,AR27,L1/L2,1,-1
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0
EALIVE,ALL
*DO,J,1,8,1
D,NELEM((AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,J),TEMP,MAX_TEMP !給每個單元的八個節(jié)點施加恒定溫度
*ENDDO
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0
NSLE,S,ALL !選擇單元中的所有節(jié)點
SF,ALL,CONV,72,25
ALLSEL,ALL
SOLVE
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0
*DO,J,1,8,1
DDELE,NELEM((AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,J),TEMP !刪除每個單元的八個節(jié)點施加恒定溫度
*ENDDO
ALLSEL,ALL
AR22=AR22+AR21
!*IF,AR22,GT,AR21*ARG1,*EXIT
FINISH
ANTYPE,TRANSIENT,RESTART
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0
*DO,J,1,8,1
DDELE,NELEM((AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,J),TEMP !刪除每個單元的八個節(jié)點施加恒定溫度
*ENDDO
TIME,AR22
*ENDDO
*ENDDO
SAVE
/POST1
PLNSOL,TEMP,,1
ANDATA,0.02,1,2,0,0,1,0,0 !生成某一范圍內(nèi)的結(jié)果數(shù)據(jù)的順序梯度線動畫
FINISH
！請各位大神幫幫忙呀
怎么解決了嘛 ·？ 這邊