三角形热源 磨削子程序 abaqus

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

    1 JLTYP,TEMP,PRESS,SNAME)

C

     INCLUDE 'ABA_PARAM.INC'

C

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

     CHARACTER*80 SNAME

     qm = 3.2321*10000

     lc = 2

     dx = (4000/60)*time(1)

     distance=coords(1)-dx

     if(distance>=0.0.and.distance<=2)then

     Flux(1) = 2*qm*distance/2

     else

     Flux(1)=0

     end if

     RETURN

     END