Product: Abaqus/Standard
User subroutine DFLUX:
can be used to define a nonuniform distributed flux as a function of position, time, temperature, element number, integration point number, etc. in a heat transfer or mass diffusion analysis;
will be called at each flux integration point for each element-based or surface-based (heat transfer only) nonuniform distributed flux definition in the analysis;
ignores any amplitude references that may appear with the associated nonuniform distributed flux definition; and
uses the nodes as flux integration points for first-order heat transfer, first-order coupled temperature-displacement, first-order coupled thermal-electrical-structural, and mass diffusion elements.
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
user coding to define FLUX(1) and FLUX(2)
RETURN
ENDFLUX(1)
Magnitude of flux flowing into the model at this point. In heat transfer cases the units are JT–1L–2 for surface fluxes and JT–1L–3 for body flux. In transient heat transfer cases where a non-default amplitude is used to vary the applied fluxes, the time average flux over the time increment must be defined rather than the value at the end of the time increment. In mass diffusion cases the units are PLT–1 for surface fluxes and PT–1 for body flux.
FLUX(1) will be passed into the routine as the magnitude of the flux specified as part of the element-based or surface-based flux definition. If the magnitude is not defined, FLUX(1) will be passed in as zero.
This flux is not available for output purposes.
FLUX(2)
In heat transfer cases: 
, the rate of change of the flux with respect to the temperature at this point. The units are JT–1L–2
–1 for surface fluxes and JT–1L–3–1 for body flux.
In mass diffusion cases: 
, the rate of change of the flux with respect to the mass concentration at this point. The units are LT–1 for surface fluxes and T–1 for body flux.
The convergence rate during the solution of the nonlinear equations in an increment is improved by defining this value, especially when the flux is a strong function of temperature in heat transfer analysis or concentration in mass diffusion analysis.
SOL
Estimated value of the solution variable (temperature in a heat transfer analysis or concentration in a mass diffusion analysis) at this time at this point.
KSTEP
Step number.
KINC
Increment number.
TIME(1)
Current value of step time (defined only in transient analysis).
TIME(2)
Current value of total time (defined only in transient analysis).
NOEL
Element number.
NPT
Integration point number in the element or on the element's surface. The integration scheme depends on whether this is a surface or a body flux.
COORDS
An array containing the coordinates of this point. These are the current coordinates if geometric nonlinearity is accounted for during the step (see “Defining an analysis,” Section 6.1.2 of the Abaqus Analysis User's Guide); otherwise, the array contains the original coordinates of the point.
JLTYP
Identifies the flux type for which this call to DFLUX is being made. The flux type may be a body flux, a surface-based flux, or an element-based surface flux. For element-based surface fluxes, this variable identifies the element face for which this call to DFLUX is being made. This information is useful when several different nonuniform distributed fluxes are being imposed on an element at the same time. See Part VI, “Elements,” of the Abaqus Analysis User's Guide for element face identification. The key is as follows:
TEMP
Current value of temperature at this integration point (defined only for a mass diffusion analysis). Temperature for a heat transfer analysis is passed in as variable SOL.
PRESS
Current value of the equivalent pressure stress at this integration point (defined only for a mass diffusion analysis).
SNAME
Surface name for a surface-based flux definition (JLTYP=0). For a body flux or an element-based surface flux the surface name is passed in as blank.