abaqus 高斯热源子程序怎么加载在圆柱体上？

我想问一下在abaqus中，如何写高斯热源的子程序加载在圆柱体上，如图所示。请问我的子程序有什么问题吗？

      SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,JLTYP,

     1                 TEMP,PRESS,SNAME)

C

      INCLUDE 'ABA_PARAM.INC'

      parameter(one=1.d0)

      DIMENSION COORDS(3),FLUX(2),TIME(2)

      CHARACTER*80 SNAME

C     Q,电弧有效热功率W

C     v,焊枪移动速度m/s

C     Rh,加热斑点半径，95%的热量落在以Rh为半径的面积内

C d,当前时刻焊接斑点中心跟焊接初始位置的距离

Q=400

v=0.005

Rh=0.002

d=v*TIME(2)

x=COORDS(1)

y=COORDS(2)

z=COORDS(3)

C     焊接板厚度为0.005m,焊枪移动从坐标0,0开始，沿着x方向移动

x0=0

y0=0

C     

PI=3.1415

C

R=sqrt((x-x0-d)**2+(y-y0)**2)

C     JLTYP＝0，表示为面热源

JLTYP=0

IF (KSTEP. eq. one) THEN

        FLUX(1)=3*Q/(PI*Rh**2)*exp(-3*R**2/Rh**2)

      ENDIF

      RETURN

      END

捕获.png