面热源子程序
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