那位高手帮帮我看下这个程序那里错了。谢谢

z471584973

/CLEAR,NOSTART
/FILNAM,2211
/TITLE,RESIDUAL STRESS ANALYSIS
/UNITS,SI
/CONFIG,NRES,10000
/PREP7
ET,1,SOLID70
MP,DENS,1,5250
MP,C,1,650
MP,KXX,1,5

/VIEW,1,1,1,1
BLOCK,0,0.04,0.02,0.04,0,0.004
block,0,0.04,0,0.02,0,0.004
esize,0.001
vsel,all
vmesh,all
allsel
!建立接触单元
ET,2,170
ET,3,174
KEYOPT,3,9,0
R,3,
RMORE,
RMORE,,150
RMORE,0
KEYOPT,3,1,2
! Generate the target surface
vsel,s,volu,,1,1
aslv,s
ASEL,S,loc,y,0.02,0.02,0
CM,_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,S,0
ESURF,ALL
allsel,all
! Generate the contact surface
vsel,s,volu,,2,2
aslv,s
ASEL,S,loc,y,0.02,0.02,0
CM,_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF,ALL
RMORE,,150
allsel,all
!接触完毕

FINISH

/SOLU
LSIZE=0.001
V=0.4 !速度
TINC=LSIZE/V !单步时间长度
KK=1 !能量集中系数
Qmax=1800000
R=0.0005
MAX_TIME=41
MAX_X=41
MAX_Y=41
   
*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME
*DO,K,1,MAX_TIME,1
  *DO,I,1,MAX_X,1
  FLUX2(I,0,K)=(I-1)*LSIZE
  *ENDDO
   
  *DO,J,1,MAX_Y,1
  FLUX2(0,J,K)=(J-1)*LSIZE
  *ENDDO
   
   *DO,I,1,MAX_X,1
    *DO,J,1,MAX_Y,1
     XCENTER=V*(K-1)*TINC
     YCENTER=0.02
     DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
     *IF,DISTANCE,LE,R,THEN
     FLUX2(I,J,K)=Qmax/EXP(KK*DISTANCE**2/(3*r**2))
     *ELSE
     FLUX2(I,J,K)=0
     *ENDIF
    *ENDDO
   *ENDDO

  FLUX2(0,0,K)=(K-1)*TINC
*ENDDO

!开始计算
  ANTYPE,TRANS
  TUNIF,20 !初始温度20度
  KBC,1
   
  *SET,DELT,1
   TIMINT,ON
   AUTOTS,OFF
   
   TT=0.1

*DO,TM,1E-6,TT,TINC
  /VIEW,1,1,1,1
  /ANG,1
  /REP,FAST
   
  FLST,2,1,5,ORDE,1
  FITEM,2,2
   
  /GO
  !*
  TIME,TM
  AUTOTS,0
  DELTIM, , , ,1
  KBC,0
  !*
  TSRES,ERASE
  !FLST,2,1,5,ORDE,1
  !FITEM,2,2
  /GO
asel,s,loc,z,0.004,0.004
  SFA,all,1,HFLUX,%FLUX2%
  /STATUS,SOLU
  SOLVE
allsel,all
  OUTRES,ALL,ALL,
  /PSF,HFLUX,,2
   
  /REPLOT
      
*ENDDO

OUTRES,ALL,ALL,

FINISH