*HEADING
 : FLUID LINK TEST FOR SS ANALYSIS WITH F2D2
 (2-1x1x1 BLOCKS, PNEUMATIC FLUID)
*RESTART, WRITE, FREQ=1
*NODE, NSET=FLUID_CAVITY-1_NODES
  1, 1., 1.
  2, 2., 1.
  3, 2., 2.
  4, 1., 2.
 11, 4., 1.
 12, 5., 1.
 13, 5., 2.
 14, 4., 2.
*ELEMENT, TYPE=F2D2, ELSET=FLUID_CAVITY-1
 1, 2, 3
 2, 3, 4
*ELEMENT, TYPE=F2D2, ELSET=FLUID_CAVITY-2
 3, 12, 13
 4, 13, 14
*ELEMENT, TYPE=FLINK, ELSET=FLUID_LINK-1
 5, 1, 11
**
**              MASS ELEMENTS
**
*ELEMENT, TYPE=MASS, ELSET=MASS1
 6, 4
*ELEMENT, TYPE=MASS, ELSET=MASS2
 7, 14
*MASS,ELSET=MASS1
 1.0, 
*MASS,ELSET=MASS2
 1.0E-6,
**
**              HYDROSTATIC FLUID PROPERTY
**
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-460.
*FLUID PROPERTY, ELSET=FLUID_CAVITY-1, REFNODE=1, TYPE=PNEUMATIC,AMBIENT=14.7
*FLUID DENSITY, PRESSURE=0., TEMPERATURE=200.
 10., 
*FLUID PROPERTY, ELSET=FLUID_CAVITY-2, REFNODE=11, TYPE=PNEUMATIC,AMBIENT=14.7
*FLUID DENSITY, PRESSURE=0., TEMPERATURE=200.
 10., 
**
**              FLUID LINK PROPERTY
**
*FLUID LINK, ELSET=FLUID_LINK-1
 10., 
**
**              SPRINGS
**
*ELEMENT, TYPE=SPRING1, ELSET=SPRINGY1
 21, 4
 41, 14
*SPRING, ELSET=SPRINGY1
 2, 
 1.E-8,
*BOUNDARY
 2,  1, 2, 0.0
 4,  1, 1, 0.0
 4,  3, 3, 0.0
 12, 1, 2, 0.0
 14, 1, 1, 0.0
 14, 3, 3, 0.0
*INITIAL CONDITIONS, TYPE=TEMPERATURE
 1, 200.
 11, 200.
*EQUATION
 2, 
 3, 1, 1.0
 2, 1,-1.0
 2, 
 3, 2, 1.0
 4, 2,-1.0
 2, 
 13, 1, 1.0, 12, 1, -1.0
 2, 
 13, 2, 1.0, 14, 2, -1.0
*NSET,NSET=out_pcav
1,
*NSET,NSET=out_ppor
1,11
**
**              STEP 1: STEADY STATE DYNAMICS
**                      WITH NO PRELOAD
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
 0.1,
*CLOAD
4, 2, -10.
*ENERGYPRINT
*NODE PRINT
 U,PCAV,PPOR,CVOL
 PU,
 V, A
*NODE FILE
 U,PU
*NODE FILE,nset=out_pcav
PCAV,
*NODE FILE,nset=out_ppor
PPOR,CVOL
*EL PRINT, ELSET=FLUID_LINK-1
 MFL,PHMFL
 MFLT,PHMFT
*EL FILE, ELSET=FLUID_LINK-1
 MFL,PHMFL
 MFLT,PHMFT
*OUTPUT, HISTORY
*ENERGYOUTPUT
*OUTPUT,FIELD
*NODE OUTPUT
U,POR,V,A
*ELEMENTOUTPUT,ELSET=FLUID_LINK-1
MFL,MFLT
*END STEP
**
**              STEP 1B: STEADY STATE DYNAMICS
**                       WITH NO PRELOAD 
**                       (NO IMAGINARY PART)
**
*STEP
*STEADY STATE DYNAMICS,DIRECT=REAL ONLY
 0.1,
*CLOAD
 4, 2, -10.
*ENERGYPRINT
*NODE PRINT
 U,PCAV,PPOR,CVOL
 PU,
 V, A
*NODE FILE
 U,PU
*NODE FILE,nset=out_pcav
PCAV,
*NODE FILE,nset=out_ppor
PPOR,CVOL
*EL PRINT, ELSET=FLUID_LINK-1
 MFL,PHMFL
 MFLT,PHMFT
*EL FILE, ELSET=FLUID_LINK-1
 MFL,PHMFL
 MFLT,PHMFT
*OUTPUT, HISTORY
*ENERGYOUTPUT
*END STEP
**
**              STEP 2: PRELOAD
**
*STEP, NLGEOM
*STATIC
 .2, 1.
*CLOAD
 4, 2, -10.
 14, 2, -10.
*ENERGYPRINT,FREQ=99
*NODE PRINT,FREQ=99
 U,PCAV,CVOL
*NODE FILE,FREQ=99
 U,
*NODE FILE,nset=out_pcav,FREQ=99
PCAV,
*NODE FILE,nset=out_ppor,FREQ=99
CVOL
*ENERGY PRINT
*OUTPUT, HISTORY
*ENERGYOUTPUT
*END STEP
**
**              STEP 3: STEADY STATE DYNAMICS
**                      NO PRELOAD
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
 0.1,
*CLOAD
 4, 2, -10.
*ENERGYPRINT
*NODE PRINT
 U,PCAV,PPOR,CVOL
 PU,
 V, A
*NODE FILE
 U,PU
*NODE FILE,nset=out_pcav
PCAV,
*NODE FILE,nset=out_ppor
PPOR,CVOL
*EL PRINT, ELSET=FLUID_LINK-1
 MFL,PHMFL
 MFLT,PHMFT
*EL FILE, ELSET=FLUID_LINK-1
 MFL,PHMFL
 MFLT,PHMFT
*OUTPUT, HISTORY
*ENERGYOUTPUT
*END STEP
**
**              STEP 3B: STEADY STATE DYNAMICS
**                       WITH NO PRELOAD 
**                       (NO IMAGINARY PART)
**
*STEP
*STEADY STATE DYNAMICS,DIRECT=REAL ONLY
 0.1,
*CLOAD
 4, 2, -10.
*ENERGYPRINT
*NODE PRINT
 U,PCAV,PPOR,CVOL
 PU,
 V, A
*NODE FILE
 U,PU
*NODE FILE,nset=out_pcav
PCAV,
*NODE FILE,nset=out_ppor
PPOR,CVOL
*EL PRINT, ELSET=FLUID_LINK-1
 MFL,PHMFL
 MFLT,PHMFT
*EL FILE, ELSET=FLUID_LINK-1
 MFL,PHMFL
 MFLT,PHMFT
*OUTPUT, HISTORY
*ENERGYOUTPUT
*END STEP