*HEADING
 : UFLUID TEST USING F2D2 
(1x1x1 BLOCK, USER DEFINED FLUID)
*RESTART, WRITE, FREQ=99
*NODE, NSET=CAVITY-1_NODES
  1, 1., 1.
  2, 2., 1.
  3, 2., 2.
  4, 1., 2.
*ELEMENT, TYPE=F2D2, ELSET=CAVITY-1
1, 2, 3
2, 3, 4
**
**              HYDROSTATIC FLUID PROPERTY
**
*FLUID PROPERTY, ELSET=CAVITY-1, REFNODE=1, TYPE=USER, NAME=CAVITY-GAS
**
**              SPRINGS
**
*ELEMENT, TYPE=SPRING1, ELSET=SPRINGY1
21, 4
*SPRING, ELSET=SPRINGY1
2, 
1.E-8,
*BOUNDARY
2,  1, 2, 0.0
4,  1, 1, 0.0
*INITIAL CONDITIONS, TYPE=TEMPERATURE
1, 200.
*EQUATION
2, 
3, 1, 1.0, 2, 1, -1.0
2, 
3, 2, 1.0, 4, 2, -1.0
**
**              STEP 1: APPLY LOAD AT T=200
**
*STEP, NLGEOM
*STATIC
.2, 1.
*CLOAD
4, 2, -100.
*NODE PRINT
U,NT11,PCAV,CVOL
*NODE FILE,FREQ=99
U,NT,PCAV,CVOL
*OUTPUT,FIELD,FREQ=99
*NODE OUTPUT
U,NT,
*ENERGY PRINT
*ENERGY FILE
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*END STEP
**
**              STEP 2: INCREASE FLUID TEMPERATURE TO 300
**
*STEP, NLGEOM
*STATIC
.5, 1.
*TEMPERATURE
1, 300.
*END STEP
**
**              STEP 3: ADD 5.0 UNITS OF MASS AT T=300
**
*STEP, NLGEOM
*STATIC
.5, 1.
*FLUID FLUX
1, 5.0
*END STEP
**
**              STEP 4: REMOVE 5.0 UNITS OF MASS AT T=300
**
*STEP, NLGEOM
*STATIC
.5, 1.
*FLUID FLUX
1, -5.0
*END STEP
**
**              STEP 5: DECREASE FLUID TEMPERATURE TO 200
**
*STEP, NLGEOM
*STATIC
.5, 1.
*TEMPERATURE
1, 200.
*FLUID FLUX,OP=NEW
*END STEP