请教版主和各位高手！

sh_0320

/CLEAR,NOSTART
/FILNAM,sample
/TITLE,RESIDUAL STRESS ANALYSIS
/UNITS,SI
/CONFIG,NRES,10000
/PREP7
ET,1,SOLID70
MP,DENS,1,7850
MP,C,1,438
MP,KXX,1,53

/VIEW,1,1,1,1
BLOCK,0,0.04,0,0.04,0,0.004
SAVE
LESIZE,1,,,80
LESIZE,4,,,80
LESIZE,9,,,4
MSHKEY,1
VATT,1,1,1
VMESH,1
/PNUM,DEFA
EPLOT
FINISH

/SOLU
LSIZE=0.0005
V=0.04 !速度
TINC=LSIZE/V !单步时间长度
KK=1 !能量集中系数
P=1400
R=0.002
r=0.01
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+(V*(K-1)*TINC)
     DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
     *IF,DISTANCE,LE,R,THEN
     FLUX2(I,J,K)=(2*P*0.4)/(3.14*R*R)/EXP(2*kk*DISTANCE**2/(0.002**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.5

*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
  !*
  !*
  SFA,P51X,1,HFLUX,%FLUX2%
  /STATUS,SOLU
  SOLVE
  OUTRES,ALL,ALL,
  /PSF,HFLUX,,2
   
  /REPLOT
      
*ENDDO

OUTRES,ALL,ALL,

FINISH

上面是激光沿斜线扫描的温度场，算出的温度有问题，只有50－60度。
在同样参数下，直线扫描温度场很正常，最高温度为900度左右。