下面是一个焊接分析的代码，转载自网络，仅供学习研究参考，本站为对结果的正确性做验证。

/PREP7
/UNITS,SI

ET,1,SOLID70
MP,C,1,544.28
MP,KXX,1,15.06
MP,DENS,1,4500

BLOCK,0,30,0,15,0,3

ESIZE,1
VMESH,ALL
FINISH

/SOLU
!PARAMETER TABLE TO INPUT HEAT FLUX AS A MOVING SOURCE
!SPEED OF HEAT SOURCE IN X DIRECTION V=(200/60)M/S
!GAUSS VARIATION ASSUMED IN X & Y DIRECTION AS Q=QMAX/EXP(K*R**2)
!TOTAL HEAT INPUT Q=3.14*QMAX/K
!DIMENSIONING TABLE NAMED "FLUX2"
V=1
L=30
NUM_TIME=L/V+1
NUM_X=L/V+1

*DIM,FLUX2,TABLE,NUM_TIME,NUM_X,7,TIME,X,Y

*DO,K,1,7,1

*DO,I,1,NUM_TIME,1
FLUX2(I,0,K)=(I-1)
*ENDDO

*DO,J,1,NUM_X,1
FLUX2(0,J,K)=(J*V-V)
*ENDDO

Q=1000
KAPA=1/4
QMAX=((Q*KAPA)/3.14159)
PERCENT=0.98
RADIUS=ABS(SQRT(1/KAPA*LOG(1/(1-PERCENT))))

!N=[RADIUS/V]

*DO,ANTIME,1,NUM_TIME,1
*DO,N,0,15,1
XDIST=(N*V)
*IF,(ANTIME-N),GT,0,THEN
*IF,SQRT((XDIST)**2+((RADIUS/6)*(K-1))**2),LE,RADIUS,THEN
!
GAUSS=QMAX/EXP(KAPA*((XDIST)**2+((RADIUS/6)*(K-1))**2))
FLUX2(ANTIME,ANTIME-N,K)=GAUSS*10E5
!
*ELSE
*ENDIF
*ENDIF
!
*IF,(ANTIME+N),LE,NUM_TIME,THEN
*IF,SQRT((XDIST)**2+((RADIUS/6)*(K-1))**2),LE,RADIUS,THEN
!
GAUSS=QMAX/EXP(KAPA*((XDIST)**2+((RADIUS/6)*(K-1))**2))
FLUX2(ANTIME,ANTIME+N,K)=GAUSS*10E5
!
*ELSE
*ENDIF
*ENDIF
*ENDDO
*ENDDO
!
!INPUT Y DIRECTION DIVISIONS
!---------------------------
FLUX2(0,0,K)=(K-1)*(RADIUS/6)
*ENDDO

ANTYPE,TRANS
TUNIF,20
KBC,1

*SET,DELT,1
TIMINT,ON
AUTOTS,OFF

*DO,TM,1E-6,30,0.5
/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