如何应用ABAQUS软件利用面体复合热源分别加载到试板的上下表面

本人应用ABAQUS软件做焊接温度场分析时，想通过dflux子程序在0.005米厚的试板的上表面加载一个面热源，然后加载一个开口位于下表面的体热源，但貌似面热源一直没有加载上去，结果中并没有面热源加载的效果。求大侠帮忙。 SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,
1 TEMP,PRESS,SNAME)
C
INCLUDE 'ABA_PARAM.INC'
parameter(one=1.d0)
DIMENSION COORDS(3),FLUX(2),TIME(2)
CHARACTER*80 SNAME
q=4000*0.6
v=0.005
d=v*TIME(2)
x=COORDS(1)
y=COORDS(2)
z=COORDS(3)
PI=3.1415
r0=3.4e-4
f=0.5
x0=0.00
y0=0.00
z0=-0.0001
if (JSTEP. eq. one) then
C JLTYP＝1，表示为体热源
JLTYP=1
h=0.011
e=2.71828
if(COORDS(3).gt.-0.0001) then
heat1=9*q/(PI*r0**2*h*(1-1/(e**3)))
zz=log(h/(z-z0))
shape1=exp(-9*((x-x0)**2+(y-y0-d)**2)/(r0**2*zz))
FLUX1=f*heat1*shape1
endif
C JLTYP＝0，表示为面热源
JLTYP=0
x1=0.00
y1=0.00
if(COORDS(3).ge.0.005) then
heat2=3*q/(PI*r0**2)
shape2=exp(-3*((x-x1)**2+(y-y1-d)**2)/r0**2)
FLUX2=(1-f)*heat2*shape2
endif
FLUX(1)=FLUX1+FLUX2
endif
RETURN
END