ABAQUS热源程序.txt

         SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,
         $                TEMP,PRESS,SNAME)

         INCLUDE 'ABA_PARAM.INC'
         
         DIMENSION COORDS(3),FLUX(2),TIME(2)
         CHARACTER*80 SNAME

C        u,焊接电压
C        i,焊接电流
C        effi,焊接效率系数
C        q,焊接有效热功率W
C        v,焊接速度m/s
         u=20
         i=165
         effi=0.5
         v=0.005
         q=u*i*effi
         d=v*TIME(2)

         x=COORDS(1)
         y=COORDS(2)
         z=COORDS(3)

C        从坐标x0,y0,z0开始,沿着z方向移动(起始点坐标）
         x0=-0.05
         y0=0
         z0=0.003

C        a1,a2,b,c为双椭球的形状参数