*HEADING
 : FLUID LINK TEST FOR SS ANALYSIS WITH F2D2
 (2-1x1x1 BLOCKS, PNEUMATIC FLUID)
*RESTART, WRITE, FREQ=1
*NODE, NSET=FLUID_CAVITY-1_NODES
*ELEMENT, TYPE=F2D2, ELSET=FLUID_CAVITY-1
*ELEMENT, TYPE=F2D2, ELSET=FLUID_CAVITY-2
*ELEMENT, TYPE=FLINK, ELSET=FLUID_LINK-1
**
**              MASS ELEMENTS
**
*ELEMENT, TYPE=MASS, ELSET=MASS1
*ELEMENT, TYPE=MASS, ELSET=MASS2
*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
**
**              SPRINGS
**
*ELEMENT, TYPE=SPRING1, ELSET=SPRINGY1
*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
*NSET,NSET=out_pcav
*NSET,NSET=out_ppor
**
**              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,CVOL
*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