面热源子程序
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=yitaUAI*1000./3.14/R0/R0
dx=v*time(1)
rr=((COORDS(1)-dx)2+COORDS(2)2
FLUX(1)=3qmEXP(-(3rr)/R0*2)
RETURN
END