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

     

1 JLTYP,TEMP,PRESS,SNAME)

      

INCLUDE 'ABA_PARAM.INC'

      

parameter(one=1.d0)

      

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

      

CHARACTER*80 SNAME

      

C     熱源功率（熱源瞬時給焊件的熱能）Q,熱效率n；熱源半徑Rh，單位m ;熱源穿透深度80微米，0.00008m

      

Q=170

      

n=1

      

Rh=0.0001

      

PI=3.1415

      

h=0.00008

     

 

C     移動速度v=0.05m/s

      

v=0.05

      

x=COORDS(1)

     

y=COORDS(2)

      

z=COORDS(3)

 

C     定義移動熱源起始點。

      

x0=-0.001

      

y0=0

      

z0=0

      

d=v*TIME(1)

C     JLTYP＝1，表示為體熱源      

      

JLTYP=1 

      

C     計算熱源熱流密度，FLUX,單位J/(m^2*s)

      

R=sqrt((x-x0-d)**2+(z-z0)**2)

      

heat=3*Q*n/(h*PI*Rh**2)

      

shape1=exp(-3*R**2/Rh**2)

     

if (JSTEP.eq.one) then

      

FLUX(1)=heat*shape1

      

    

      

RETURN

      

END