ANSYS温度成命令流问题

xiaoqiang1

有哪位大神帮忙看一看下面的命令流哪有问题，时间很紧 ，在此跪谢了！！！！
命令流如下：
FINISH
/FILNAME,TEMP-FIELD
/PREP7
ET,1,SOLID70
ET,2,SURF152
MPTEMP,1,20,200,400,800,1200
MPDATA,KXX,1,,51.8,48.6,43.9,25.9,39.7
MPDATA,KXX,2,,14.1,18,20.8,26.3,31.9
MPDATA,C,1,,461,523,607,954,661
MPDATA,C,2,,461,523,607,644,669
MP,DENS,1,7850
MP,DENS,2,7930
K,1,0,0,0
K,2,0,2,0
K,3,2,0,0
K,4,4,0,0
K,5,20,0,0,
K,6,20,-10,0
K,7,0,-10,0
K,8,0,-2,0
K,9,0,0,80
L,1,2
L,1,3
L,3,4
L,4,5
L,5,6
L,6,7
L,7,8
L,1,8
L,1,9
LARC,2,3,1,2
A,1,2,3
A,1,3,4,5,6,7,8

LSEL,S,,,2,4
LCCAT,ALL     !连接线，以便VDRAG成体
ALLSEL
LCCAT,7,8
ALLSEL
LESIZE,2,,,8
LESIZE,3,,,4
LSEL,S,,,4,5
LESIZE,ALL,,,8
ALLSEL
LESIZE,6,,,20
LSEL,S,,,7,8
LESIZE,ALL,,,4
ALLSEL
LSEL,S,,,1,2
LSEL,A,,,10
LESIZE,ALL,,,8
ALLSEL
LESIZE,9,1
MSHK,1
AMESH,ALL
LDELE,11,12
MSHAPE,0,3D
MSHK,1
MAT,1
TYPE,1
VDRAG,2,,,,,,9,,,,,,
MSHAPE,0,3D
MSHK,1
MAT,2
TYPE,1
VDRAG,1,,,,,,9,,,,,,
V=1
L=80
NUM_TIME=L/v+1
NUM_Z=L/v+1
*DIM,FLUX2,TABLE,NUM_TIME,NUM_Z,9,TIME,Z,X
*DO,K,1,9,1
*DO,I,1,NUM_TIME,1
FLUX2(I,0,K)=(I-1)
*ENDDO
*DO,J,1,NUM_Z,1
FLUX2(0,J,K)=(J-1)
*ENDDO
Q=720
KAPA=1/4
qmax=((Q*KAPA)/3.14159)
PERCENT=0.98
RADIUS=ABS(SQRT(1/KAPA*LOG(1/(1-PERCENT))))  
*DO,ANTIME,1,NUM_TIME,1

*DO,N,0,8,1
ZDIST=(N*V)
*IF,(ANTIME-N),GT,0,THEN
*IF,SQRT((ZDIST)**2+((RADIUS/8)*(K-1))**2),LE,RADIUS,THEN
GAUSS=qmax/EXP(KAPA*((ZDIST)**2+((RADIUS/8)*(K-1))**2))
FLUX2(ANTIME,ANTIME-N,K)=GAUSS*10e5
*ELSE
*ENDIF
*ENDIF
*IF,(ANTIME+N),LE,NUM_TIME,THEN
*IF,SQRT((ZDIST)**2+((RADIUS/8)*(K-1))**2),LE,RADIUS,THEN
GAUSS=qmax/EXP(KAPA*((ZDIST)**2+((RADIUS/8)*(K-1))**2))
FLUX2(ANTIME,ANTIME+N,K)=GAUSS*10e5

*ELSE
*ENDIF
*ENDIF
*ENDDO
*ENDDO

FLUX2(0,0,K)=(K-1)*(RADIUS/8)
*ENDDO
/SOLU
ANTYPE,TRANS  !瞬态分析
timint,on  !时间积分效果，如果为OFF则进行稳态分析
autots,off  ! 关闭时间步长
kbc,1 !0表示渐变,1阶跃
TUNIF,20 !所有节点初始温度
asel,s,,,2,6,1
ASEL,A,,,10
sfa,all,,conv,6,20
ALLSEL
ASEL,S,,,7
DA,7,TEMP,0
allsel
nropt,full !牛顿拉普森法
ESEL,S,MAT,,2        
EKILL,ALL
ALLSEL
*do,tm,1e-6,80,1
/VIEW,1,1,1,1
AUTOTS,on !开启时间步长
TIME,tm !指定载荷步结束时间
DELTIM,,,,1 !定义时间步长
KBC,0 !0表示渐变,1阶跃
*IF,TM,GT,1e-6,THEN
ZALIVE=TM*V+1
NSEL,S,LOC,Z,(ZALIVE-1)
ESLN
NSEL,S,LOC,Z,(ZALIVE-2)
ESLN,U
ESEL,R,MAT,,2
EALIVE,ALL
ALLSEL
*ELSE
NSEL,S,LOC,Z,0
ESLN
ESEL,R,MAT,,2
EALIVE,ALL
ALLSEL
*ENDIF
TSRES,ERASE
ASEL,S,,,3,5,1
SFA,all,1,HFLUX,%FLUX2%
SOLVE
ALLSEL

/PSF,hflux,,3
ESEL,S,LIVE
eplot
/replot
*enddo