编辑代码

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

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

  U=16.5
  AI=60.
  v=250./60
  yita=0.75
  R0=5.

  qm=yita*U*AI*1000./3.14/R0/R0

  dx=v*time(1)

  rr=((COORDS(1)-dx)**2+COORDS(2)**2

  FLUX(1)=3*qm*EXP(-(3*rr)/R0**2)

  RETURN
  END