程序如下，熱源形狀參數(shù)abc和焊接電流、電壓、焊接速度以及關(guān)鍵點(diǎn)的坐標(biāo)和時(shí)間可以根據(jù)你的條件修改。

      SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,
     1                 TEMP,PRESS,SNAME)
C
      INCLUDE 'ABA_PARAM.INC'

      DIMENSION COORDS(3),FLUX(2),TIME(2)
      CHARACTER*80 SNAME
          
          x=COORDS(1)
      y=COORDS(2)
      z=COORDS(3)
          PI=3.1415926
          
                    
      !焊接工藝參數(shù)，使用同一熱源
      wu=60.0
      wi=1.0
      effi=1.0
      v=0.1
      q=wu*wi*effi
          
      !焊接熱源形狀參數(shù)(球熱源)
      a=0.0005
          b=0.0005
          c=0.0005


      !各個(gè)關(guān)鍵點(diǎn)坐標(biāo)及焊槍經(jīng)過(guò)時(shí)間(相對(duì)總時(shí)間累計(jì)shijian)
      x1=-0.00255
      y1=0.00095
      z1=0.0002
          t1=0.0
      
      x2=0.00245
      y2=0.00095
      z2=0.0002
          t2=0.05
      


      !熱源分布
          heat=6*sqrt(3.0)*q/(a*b*c*PI*sqrt(PI))
        
        !第一段直線
        if(TIME(2).lt.t2) then      
          d=v*TIME(2)
        !熱源中心在總體坐標(biāo)系的坐標(biāo)
          xd=x1+d
          yd=y1
          zd=z1
        !熱源局部坐標(biāo)
          xx=x-xd
          yy=y-yd
          zz=z-zd
         
           shape=exp(-3.0*xx**2/a**2-3.0*yy**2/b**2-3.0*zz**2/c**2)
            


         
      
        end if
            
C     JLTYP＝1，表示為體熱源
      JLTYP=1                                                                                                               
      FLUX(1)=heat*shape
      RETURN
      END，