ABAQUS Analysis User's Manual  

25.2.23 UFLUID
User subroutine to define fluid density and fluid compliance for hydrostatic fluid elements.

Product: ABAQUS/Standard  

References

Overview

User subroutine UFLUID:

  • is called for each cavity for which a user-defined fluid constitutive model is being specified;

  • is called for every fluid element (“Hydrostatic fluid elements,” Section 18.8.1) and for every fluid link element (“Fluid link elements,” Section 18.8.3) connected to a cavity reference node;

  • requires that the fluid density, , and the fluid pressure compliance, , be defined;

  • requires that the fluid temperature compliance, , be defined if the routine is to be used in a linear perturbation step and the fluid is subjected to a temperature excursion; and

  • ignores any data specified for the fluid constitutive model outside the user subroutine.

Density and fluid mass

At the start of the analysis (prior to the first iteration) the density calculated in user subroutine UFLUID (for the initial pressure, , and temperature, ) is used to calculate the fluid mass from the initial cavity volume. During the analysis the expected cavity volume is calculated from the fluid mass and the density.

User subroutine interface

      SUBROUTINE UFLUID(RHO,CP,CT,PNEWDT,ENER,PRESS,DPRESS,PRESSI,
     1 TEMP,DTEMP,TEMPI,TIME,DTIME,KSTEP,KINC,NONUM,FLNAME,LFLAG)
C
      INCLUDE 'ABA_PARAM.INC'
C
      CHARACTER*80 FLNAME
      DIMENSION TIME(2)

      user coding to define RHO, CP, and CT

      RETURN
      END

Variables to be defined

RHO

Fluid density, , at the end of the increment.

CP

Fluid pressure compliance, , at the end of the increment. For a linear perturbation step this is the base state compliance. Fluid pressure compliance is defined as

where is the fluid cavity pressure.

CT

Fluid temperature compliance, . This variable is needed only if a fluid temperature excursion occurs in a linear perturbation step and is the base state compliance. Fluid temperature compliance is defined as

where is the fluid cavity temperature.

Variables that can be updated

PNEWDT

Ratio of suggested new time increment to the time increment being used (DTIME, see below). This variable allows you to provide input to the automatic time incrementation algorithms in ABAQUS/Standard (if automatic time incrementation is chosen).

PNEWDT is set to a large value before each call to UFLUID.

If PNEWDT is redefined to be less than 1.0, ABAQUS/Standard must abandon the time increment and attempt it again with a smaller time increment. The suggested new time increment provided to the automatic time integration algorithms is PNEWDT × DTIME, where the PNEWDT used is the minimum value for all calls to user subroutines that allow redefinition of PNEWDT for this iteration.

If PNEWDT is given a value that is greater than 1.0 for all calls to user subroutines for this iteration and the increment converges in this iteration, ABAQUS/Standard may increase the time increment. The suggested new time increment provided to the automatic time integration algorithms is PNEWDT × DTIME, where the PNEWDT used is the minimum value for all calls to user subroutines for this iteration.

If automatic time incrementation is not selected in the analysis procedure, values of PNEWDT that are greater than 1.0 will be ignored and values of PNEWDT that are less than 1.0 will cause the job to terminate.

ENER

Energy per unit mass stored in the fluid. This variable is used for energy output only and has no effect on the solution.

Variables passed in for information

PRESS

Fluid cavity pressure at the end of the increment. For a linear perturbation step this is the base state pressure.

DPRESS

Fluid cavity pressure increment. For a linear perturbation step this value is zero.

PRESSI

Fluid cavity pressure at the beginning of the analysis.

TEMP

Fluid cavity temperature at the end of the increment. For a linear perturbation step this is the base state temperature.

DTEMP

Fluid cavity temperature increment. For a linear perturbation step this value is zero.

TEMPI

Fluid cavity temperature at the beginning of the analysis.

TIME(1)

Current value of step time at the start of the increment.

TIME(2)

Current value of total time at the start of the increment.

DTIME

Time increment.

KSTEP

Step number.

KINC

Increment number.

NONUM

Cavity reference node number.

FLNAME

User-specified fluid property name, left justified.

LFLAG

Linear perturbation flag for the step. If this is a linear perturbation step, LFLAG=1. For a general analysis step LFLAG=0.

Example