我这个焊接命令流的错误在那里？

wangsheng040

大家好！我最近做焊接，但是这个命令流有错误，不知道错误在那里？哪个高人帮我看看！谢谢
finish
/clear
/title,weld analysis by "ELEMENT BIRTH AND DEATH"
/prep7
/unit,si
et,1,13,4
et,2,13,4

MPTEMP,1,20,500,800,1200,1500                              !温度
MPDATA,EX,1,1,2.09E11,1.72E11,1.33E11,0.84E11,0.45E11      !左边板的弹性模量
MPDATA,EX,2,1,2.16E11,1.78E11,1.42E11,1.04E11,0.68E11      !焊缝的弹性模量
MPDATA,EX,3,1,2.12E11,1.75E11,1.39E11,1.07E11,0.83E11      !右边板的弹性模量
TB,BKIN,1,5                                                !以表格的形式输入经典双线性随动强化参数
TBTEMP,20,1                                                !左板材料特性
TBDATA,1,315E6,0.209E11
TBTEMP,500,2
TBDATA,1,212E6,0.172E11
TBTEMP,800,3
TBDATA,1,153E6,0.133E11
TBTEMP,1200,4
TBDATA,1,93E6,0.084E11
TBTEMP,1500,5
TBDATA,1,33E6,0.045E11

TB,BKIN,2,5                                                 ! 焊缝材料特性
TBTEMP,20,1
TBDATA,1,314E6,0.216E11
TBTEMP,500,2
TBDATA,1,203E6,0.178E11
TBTEMP,800,3
TBDATA,1,153E6,0.142E11
TBTEMP,1200,4
TBDATA,1,63E6,0.104E11
TBTEMP,1500,5
TBDATA,1,13E6,0.068E11

TB,BKIN,3,5                                                 !右板的材料特性
TBTEMP,20,1
TBDATA,1,330E6,0.212E11
TBTEMP,500,2
TBDATA,1,213E6,0.175E11
TBTEMP,800,3
TBDATA,1,153E6,0.139E11
TBTEMP,1200,4
TBDATA,1,73E6,0.107E11
TBTEMP,1500,5
TBDATA,1,13E6,0.083E11

MP,DENS,1,7850                                             !材料的密度
MP,DENS,2,7770
MP,DENS,3,7860
MP,ALPX,1,1.23E-5                                          !材料的线膨胀系数
MP,ALPX,2,1.38E-5
MP,ALPX,3,1.48E-5
MP,NUXY,1,0.30                                             
MP,NUXY,2,0.29
MP,NUXY,3,0.29
MP,KXX,1,33                                                !材料热导率
MP,KXX,2,32
MP,KXX,3,34
MP,C,1,561                                                 ! 材料的比热容
MP,C,2,683
MP,C,3,983
MP,MURX,1,1                                                !材料磁场相对系数
MP,MURX,2,1
MP,MURX,3,1
MP,REFT,1,20                                               !材料的参考温度
MP,REFT,2,1500
MP,REFT,3,20

!建模和网格化
*AFUN,DEG
CSYS,0
WIDTH=1
HEIGHT=0.3
EX=0.01*TAN(60)
AR=0.3*TAN(30)
DF=AR+EX
CF=DF/COS(60)

K,1,0,0,0
K,2,WIDTH/2-0.01,0,0
K,3,WIDTH,0,0
K,4,WIDTH,HEIGHT,0
K,5,WIDTH/2+AR,HEIGHT,0
K,6,WIDTH/2-AR,HEIGHT,0
K,7,0,HEIGHT,0
K,8,WIDTH/2+0.01,0,0
K,9,WIDTH/2-0.01,0.02,0
K,10,WIDTH/2+0.01,0.02,0

N,1,WIDTH/2+AR,HEIGHT,0
N,2,WIDTH,HEIGHT,0
N,3,WIDTH/2,0,0

LARC,5,6,2,CF
A,1,2,9,6,7
A,2,9,6,5,10,8
A,8,3,4,5,10
/TRIAD,OFF                      !关闭坐标的三角符号显示
/PLOTS,INFO,0                   !关闭图形屏幕上的信息显示
/DEVICES,VECTOR,1               !采用向量显示

ESIZE,0.03
TYPE,2
MAT,2
AMESH,2
ESIZE,0.03
TYPE,1
MAT,1
AMESH,1
MAT,3
AMESH,3
NUMMRG,ALL,,,,LOW
NUMCMP,ALL
EPLOT
/IMAGE,SAVE,VMESH,BMP
/DEVICES,VECTOR,0
FINISH

/SOLU
ANTYPE,4
TRNOPT,FULL
NSEL,ALL
*GET,MINX,NODE,,MNLOC,X
NSEL,S,LOC,X,MINX
D,ALL,UX,0
*GET,MINY,NODE,,MNLOC,Y
NSEL,R,LOC,Y,MINY
D,ALL,UY,0
NSEL,ALL
*GET,MINX,NODE,,MNLOC,X
NSEL,S,LOC,X,MINX
D,ALL,TEMP,20
NSEL,ALL
*GET,MAXX,NODE,,MXLOC,X
NSEL,S,LOC,X,MAXX
D,ALL,TEMP,20

*SET,nna,2
ESEL,ALL
*GET,EMAX,ELEM,,NUM,MAX
ASEL,S,AREA,,nna
ESLA
*GET,NSE,ELEM,,COUNT
*DIM,NE,,NSE
*DIM,N1,,NSE
*DIM,N2,,NSE
*DIM,NEX,,NSE
*DIM,NEY,,NSE
*DIM,NEORDER,,NSE

II=0
*DO,I,1,EMAX
*IF,ESEL(I),EQ,1,THEN
II=II+1
NE(II)=I
*ENDIF
*ENDDO

*DO,I,1,NSE
*GET,NEY(I),ELEM,NE(I),CENT,Y
*GET,NEX(I),ELEM,NE(I),CENT,X
*ENDDO

D=0
*DO,I,1,NSE
*IF,NEY(I),GE,HEIGHT,THEN
d=d+1
N2(d)=NE(I)
*ENDIF
*ENDDO


!ESEL,NONE
!*DO,I,1,d
!*IF,N2(I),NE,0,THEN
!ESEL,A,ELEM,,N2(I)
!*ENDIF
!*ENDDO
!EPLOT

MINE=0
*DIM,NE2,,d
*DIM,NEX2,,d
*DIM,NEY2,,d
*DO,I1,1,d
ESEL,U,ELEM,,MINE
*GET,NSE2,ELEM,,COUNT
II=0
*DO,I1,1,EMAX
*IF,ESEL(I),EQ,1,THEN
II=II+1
NE2(II)=I
*ENDIF
*ENDDO
*DO,I,1,NSE2
*GET,NEY2(I),ELEM,NE2(I),CENT,Y
*GET,NEX2(I),ELEM,NE2(I),CENT,X
*ENDDO

MINY=1E20
MINX=1E20

*DO,I,1,NSE2
*IF,NEY2(I),LT,MINY,THEN
MINY=NEY2(I)
MINX=NEX2(I)
MINE=NE2(I)
*ELSE
*IF,NEY2(I),EQ,MINY,THEN
*IF,NEX2(I),LT,MINT,THEN
MINY=NEY2(I)
MINX=NEX2(I)
MINE=NE2(I)
*ENDIF
*ENDIF
*ENDIF
*ENDDO
NEORDER(I1)=MINE
*ENDDO


MAX_TEM=1500
dt1=1e-3
dt=5
t=0
esel,all
eplot
/auto,1
/replot
*do,i,1,nse
ekill,neorder(i)
esel,s,live
eplot
*enddo
allsel,all
outres,all,all
ic,all,temp,20
kbc,1
timint,0,struct
timint,1,therm
TIMINT,0,MAG
TINTP,0.005,,,1,0.5,0.2

NSUB1=2
NSUB2=40
*DO,I,1,NSE
EALIVE,NEORDER(I)
ESEL,S,LIVE
EPLOT
ESEL,ALL
t=t+Dt1
TIME,t
NSUBST,1
*DO,J,1,4
D,NELEM(NEORDER(I),J),TEMP,MAX_TEM
*ENDDO
SOLVE
t=t+Dt1
TIME,t
SOLVE
*DO,J,1,4
DDELE,NELEM(NEORDER(I),J),TEMP
*ENDDO
t=t+Dt-2*Dt1
TIME,t
NSUBST,NSUB1
SOLVE
*ENDDO
t=t+50000
TIME,t
NSUBST,NSUB2
SOLVE
FINISH
!节点379的温度随时间变化的曲线
/POST26
NSOL,2,379,TEMP,,TEMP379
/GROPT,AXNSC,1.5,
PLVAR,2
/IMAGE,SAVE,TEMP379,BMP
FINISH
/POST1
/SEG,DELE
/CONT,1,15,0,1200E6/16,1200E6
/DSCALE,1,1.0
AVPRIN,0,0
AVRES,1
/SEG,MULTI,STRESS,0.1
ESEL,ALL
*DO,I,1,NSE
ESEL,U,ELEM,,NEORDER(I)
*ENDDO
*DO,I,1,NSE
ESEL,A,ELEM,,NEORDER(I)
SET,(I-1)*3+1,1
PLNSOL,S,EQV
*DO,J,1,NSUB1
SET,(I-1)*3+3,J
PLNSOL,S,EQV
*ENDDO
*ENDDO
*DO,I,1,NSUB2
SET,(NSE-1)*3+4,I
PLNSOL,S,EQV
*ENDDO
/SEG,OFF,STRESS,0.1
/IMAGE,SAVE,STRESS,BMP   

/POST1
/SEG,DELE
/CONT,1,15,0,1500/16,1500
/DSCALE,1,1.0
AVPRIN,0,0
AVRES,1
/SEG,MUNTI,TEMP,0.1
ESEL,ALL
*DO,I,1,NSE
ESEL,U,ELEM,,NEORDER(I)
*ENDDO
*DO,I,1,NSE
ESEL,A,ELEM,,NEORDER(I)
SET,(I-1)*3+1,1
PLNSOL,TEMP
*DO,J,1,NSUB1
SET,(I-1)*3+3,J
PLNSOL,TEMP
*ENDDO
*ENDDO