*HEADING
 AIRSPRING -- AXISYMMETRIC MODEL
*NODE
*NODE, NSET = RIGID
*NSET, NSET = CAVITY
*NSET, NSET=ALL, GENERATE
*NGEN, LINE = C
**
** RUBBER SPRING
**
*ELEMENT, TYPE = MAX1
*ELGEN, ELSET = SPRING
*MEMBRANE SECTION, ELSET = SPRING, MATERIAL=RUBBER
5.0E-3,
*REBAR LAYER
RBAR, 3.5E-3, 3.5E-3, 0, STEEL, 8.66,1
*MATERIAL, NAME = RUBBER
*HYPERELASTIC, POLYNOMIAL, N = 1
3.2E6, 0.8E6
*DENSITY
1000.,
*MATERIAL, NAME = STEEL
*ELASTIC
2.5086E9,0.3
**
** RIGID BODY/SURFACE
**
*NODE,NSET=RIGNODE
*NGEN,NSET=RIGNODE,LINE=C
*ELEMENT,TYPE=RAX2,ELSET=RIGELEM
*ELGEN,ELSET=RIGELEM
*RIGID BODY, REF NODE=1000, ELSET=RIGELEM
**
** FLUID ELEMENTS
**
*ELEMENT, TYPE = FAX2, ELSET = FLUID
*ELGEN, ELSET = FLUID
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16
*FLUID PROPERTY, REF NODE=500, TYPE=PNEUMATIC,ELSET=FLUID, AMBIENT=101.36E+3
*FLUID DENSITY, PRESSURE=0.0, TEMPERATURE=27.0
1.774,
**
*BOUNDARY
1,1,2
4,1,2
4,6,6
6,1,1
RIGID,1,2
RIGID,6,6
*INITIAL CONDITIONS, TYPE=STRESS, REBAR, UNBALANCED STRESS=STEP
SPRING,,70.
*NSET,NSET=RIG
*SURFACE,TYPE=ELEMENT,NAME=ASURF
ECON,SPOS
*SURFACE,TYPE=ELEMENT,NAME=BSURF
RIGELEM,SNEG
**
** CONTACT DEFINITION
**
*ELSET,ELSET=ECON,GENERATE
*CONTACT PAIR,INTERACTION=SMOOTH
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
************** HISTORY DEFINITION **********************
**
** STEP 0 - EQUILIBRATE MEMBRANE ELEMENTS INITIAL STRESS
**
*STEP, NLGEOM
*STATIC
*BOUNDARY
ALL,1,6
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ENDSTEP
**
** STEP 1 - INFLATION STEP (5 ATMOSPHERES)
**
*STEP, NLGEOM
*STATIC
0.1,1.0
*MONITOR, NODE=26, DOF=1
*PRINT, CONTACT = YES
*BOUNDARY, OP=NEW
CAVITY,8,8, 5.066E5
1,1,2
4,1,2
4,6,6
6,1,1
RIGID,1,2
RIGID,6,6
*OUTPUT,FIELD
*NODE OUTPUT
U,RF
CF,
*ELEMENT OUTPUT
S,
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF
*EL PRINT, FREQUENCY=99
*NODE PRINT, NSET=CAVITY, FREQUENCY=99
PCAV, CVOL
*NODE PRINT, NSET=RIGID, FREQUENCY=99
U, RF
*NODE FILE, NSET=CAVITY
PCAV, CVOL
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=CAVITY
PCAV,CVOL
*NODE OUTPUT, NSET=RIGID
U, RF
CF,
*END STEP
**
** STEP 2 - SEAL CAVITY, RELEASE GRIP ON
**          RIGID SURFACE,
**          AND LOAD RIGID SURFACE
**
*STEP, NLGEOM
*STATIC
1.0, 1.0
*MONITOR, NODE=1000, DOF=2
*BOUNDARY, OP=NEW
1,1,2
4,1,2
4,6,6
6,1,1
RIGID,1
RIGID,6
*CLOAD
RIGID,2,-150.E3
*END STEP
**
** STEP 3 - PERTURBATION STEP
**          2 LOAD CASES
*STEP, PERTURBATION
*STATIC
*BOUNDARY, OP=NEW
1,1,2
4,1,2
4,6,6
6,1,1
RIGID,1
RIGID,6
** LOAD CASE 1
** UNIT AXIAL DISPLACEMENT WITH CAVITY PRESSURE FREE
*LOAD CASE, NAME="PRESSURE FREE"
*BOUNDARY, OP=NEW
RIGID, 2, 2, 1.E-3
*END LOAD CASE
** LOAD CASE 2
** UNIT LOAD WITH CAVITY PRESS FIXED
*LOAD CASE, NAME="PRESSURE FIXED"
*BOUNDARY, OP=NEW
CAVITY, 8, 8, 0
RIGID,  2, 2, 1.E-3
*END LOAD CASE
*END STEP
**
** STEP 4 - NONLINEAR LOAD DEFLECTION CURVE
**
*STEP, NLGEOM
*STATIC
0.1,1.0
*CLOAD
RIGID,2,-240.E3
*END STEP