我想問(wèn)一下在abaqus中，如何寫高斯熱源的子程序加載在圓柱體上，如圖所示。請(qǐng)問(wèn)我的子程序有什么問(wèn)題嗎？

      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,焊槍移動(dòng)速度m/s

C     Rh,加熱斑點(diǎn)半徑，95%的熱量落在以Rh為半徑的面積內(nèi)

C d,當(dāng)前時(shí)刻焊接斑點(diǎn)中心跟焊接初始位置的距離

Q=400

v=0.005

Rh=0.002

d=v*TIME(2)

x=COORDS(1)

y=COORDS(2)

z=COORDS(3)

C     焊接板厚度為0.005m,焊槍移動(dòng)從坐標(biāo)0,0開始，沿著x方向移動(dòng)

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