*HEADING
 : F3D4 FLUID ELEMENT TEST 
(1x1x1 BLOCK, HYDRAULIC, NO BULK MODULUS, THERMAL EXPANSION)
*RESTART, WRITE, FREQ=99
*NODE, NSET=CAVNOD
  1, 1., 1., 1.
  2, 2., 1., 1.
  3, 2., 2., 1.
  4, 1., 2., 1.
  5, 1., 1., 0.
  6, 2., 1., 0.
  7, 2., 2., 0.
  8, 1., 2., 0.
*ELEMENT,TYPE=SFM3D4,ELSET=STRUCTURE
 1, 2, 3, 7, 6
 2, 3, 4, 8, 7
 3, 6, 7, 8, 5
*SURFACE SECTION, ELSET=STRUCTURE
*SURFACE,TYPE=ELEMENT,NAME=CAV
 STRUCTURE, SPOS
*PHYSICAL CONSTANTS, ABSOLUTE ZERO = 0
*FLUID CAVITY, NAME=CAVITY, BEHAVIOR=FLUID, SURFACE=CAV, REFNODE=1
 1.0
*FLUID BEHAVIOR, NAME=FLUID
*FLUID DENSITY
10., 
*FLUID BULK MODULUS
4.E+5,   0.
1.E+5, 300.
*FLUID EXPANSION, ZERO=100.
1.E-5,   0.
4.E-5, 300.
*ELEMENT, TYPE=SPRING1, ELSET=SPRINGX1
21, 2
*SPRING, ELSET=SPRINGX1
1, 
400., 
*ELEMENT, TYPE=SPRING1, ELSET=SPRINGY1
42, 4
*SPRING, ELSET=SPRINGY1
2, 
1.E-6,
*BOUNDARY
2, 2, 3, 0.0
3, 3, 3, 0.0
4, 1, 1, 0.0
4, 3, 3, 0.0
5, 1, 3, 0.0
6, 2, 3, 0.0
7, 3, 3, 0.0
8, 1, 1, 0.0
8, 3, 3, 0.0
*EQUATION
2, 
3, 1, 1.0, 2, 1, -1.0
2, 
6, 1, 1.0, 2, 1, -1.0
2, 
7, 1, 1.0, 2, 1, -1.0
2, 
3, 2, 1.0, 4, 2, -1.0
2, 
7, 2, 1.0, 4, 2, -1.0
2, 
8, 2, 1.0, 4, 2, -1.0
*INITIAL CONDITIONS, TYPE=TEMPERATURE
1, 200.
**
**              STEP 1: APPLY LOAD AT T=200
**
*STEP, NLGEOM
*STATIC
.2, 1.
*CLOAD
4, 2, -600.
*NODE PRINT
U,NT11,PCAV,CVOL
*OUTPUT,FIELD,FREQ=99
*NODE OUTPUT
 U,NT,
*OUTPUT,HISTORY,FREQ=99
*NODE OUTPUT, NSET=CAVNOD
U,NT,PCAV,CVOL
*OUTPUT,HISTORY
*ENERGY OUTPUT
*ENERGY PRINT
*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