*HEADING
AIRSPRING -- 180 DEG MODEL 
*NODE, NSET = FPLANE
1,2.E-2,0.,0.
4,2.E-2,0.,20.E-2,0.,0.,1.0    
6,0.,0.,20.E-2,0.,0.,1.0    
26,0.,0.,40.E-2,0.,0.,-1.0 
32,6.293E-2,0.,37.771E-2,-0.58431,0.,-0.81153 
38,10.698E-2,0.,36.211E-2,0.,0.,-1.0  
39,11.35E-2,0.,36.211E-2,0.,0.,-1.0 
40,12.E-2,0.,36.211E-2,0.,0.,-1.0 
41,12.E-2,0.,0.       
*NGEN, LINE = C, NSET = FPLANE
6, 26, 2, , 0., 0., 30.E-2, 0., 1., 0.
26, 32, 1, , 0., 0., 30.E-2    
32, 38, 1, , 10.698E-2, 0., 43.211E-2 
*NCOPY,OLDSET=FPLANE,CHANGE NUMBER=100,SHIFT,
MULTIPLE=22,NEWSET=NALL
0., 0., 0.
0., 0., 0., 1., 0., 0., 8.18181818
*NODE, NSET = RIGID
11000, 12.E-2, 0. , 0.
*NODE, NSET = CAVITY
50000, 9.E-2 , 0. , 0.
**
** RUBBER SPRING
**
*ELEMENT, TYPE = S4R
2,   4,  6, 106, 104
26, 26, 27, 127, 126
*ELGEN, ELSET = SPRING1
2,  11, 2, 2, 22, 100, 100
*ELGEN, ELSET = SPRING2
26, 14, 1, 1, 22, 100, 100
*ELSET,ELSET=SPRING
SPRING1,SPRING2
*SHELL SECTION, MATERIAL = RUBBER, ELSET = SPRING
6.0E-3, 
*MATERIAL, NAME = RUBBER
*HYPERELASTIC, POLYNOMIAL, N = 1
3.2E6, 0.8E6
*DENSITY
1000., 
*REBAR,ELEMENT=SHELL,MATERIAL=STEEL,GEOMETRY=SKEW,
NAME=PLSREBAR
SPRING, 1.E-6, 3.5E-3, 0.,  18.0
*REBAR,ELEMENT=SHELL,MATERIAL=STEEL,GEOMETRY=SKEW,
NAME=MNSREBAR
SPRING, 1.E-6, 3.5E-3, 0., -18.0
*MATERIAL, NAME = STEEL
*ELASTIC
2.1E11, 0.3
*DENSITY
7.8E3, 
**
** FLUID ELEMENTS
**
*ELEMENT, TYPE = F3D3, ELSET = FLUID
6001,  1,  4, 104
6040, 40, 41, 140
*ELGEN, ELSET = FLUID
6001, 22, 100, 100
6040, 22, 100, 100
*ELEMENT, TYPE = F3D4, ELSET = FLUID
6002,  4,  6, 106, 104
6026, 26, 27, 127, 126
*ELGEN, ELSET = FLUID
6002, 11, 2, 2, 22, 100, 100
6026, 14, 1, 1, 22, 100, 100
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16
*FLUID PROPERTY, REF NODE=50000, TYPE=PNEUMATIC,
ELSET=FLUID, 
AMBIENT=101.36E+3
*FLUID DENSITY, PRESSURE=0.0, TEMPERATURE=27.0
1.774, 
**
** RIGID BODY/SURFACE
**
*NODE, NSET = RIGID
11000, 12.E-2, 0., 0.
*NODE,NSET=RIGNODE
11001, 12.E-2,  0.,  36.511E-2
11021, 6.5E-2,  0.,  42.011E-2
12000, 12.E-2,  0.,  42.011E-2
*NGEN,NSET=RIGNODE,LINE=C
11001, 11021, 1, 12000
*NCOPY,OLDSET=RIGNODE,CHANGE NUMBER=100,SHIFT,
MULTIPLE=22,NEWSET=NALL
0., 0., 0.
0., 0., 0., 1., 0., 0., 8.18181818
*ELEMENT,TYPE=R3D4,ELSET=RIGELEM
11001, 11001,11101,11102,11002
*ELEMENT, TYPE=R3D3, ELSET=RIGELEM
112001, 11000, 11101, 11001
112002, 11000, 11201, 11101
112003, 11000, 11301, 11201
112004, 11000, 11401, 11301
112005, 11000, 11501, 11401
112006, 11000, 11601, 11501
112007, 11000, 11701, 11601
112008, 11000, 11801, 11701
112009, 11000, 11901, 11801
112010, 11000, 12001, 11901
112011, 11000, 12101, 12001
112012, 11000, 12201, 12101
112013, 11000, 12301, 12201
112014, 11000, 12401, 12301
112015, 11000, 12501, 12401
112016, 11000, 12601, 12501
112017, 11000, 12701, 12601
112018, 11000, 12801, 12701
112019, 11000, 12901, 12801
112020, 11000, 13001, 12901
112021, 11000, 13101, 13001
112022, 11000, 13201, 13101
*ELGEN,ELSET=RIGELEM
11001, 22,100,100,  20,1,1
*ELEMENT,TYPE=MASS,ELSET=MASS_RIGID
11000, 11000
*MASS,ELSET=MASS_RIGID
.5, 
**
** KINEMATIC CONSTRAINTS
**
*NSET, NSET = N1M1, GENERATE
101, 2201, 100
*NSET, NSET = N4, GENERATE
4, 2204, 100
*NSET, NSET = N6, GENERATE
6, 2206, 100
*NSET, NSET = N40, GENERATE
40, 2240, 100
*NSET, NSET = N41, GENERATE
41, 2241, 100
*NSET, NSET = YPLN1, GENERATE
4,  26, 2
26, 39, 1
*NSET, NSET = YPLN2, GENERATE
2204, 2226, 2
2226, 2239, 1
*BOUNDARY
1,1,6
N4,1,6
N6,2,3
YPLN1, 2
YPLN1, 4
YPLN1, 6
YPLN2, 2
YPLN2, 4
YPLN2, 6
RIGID,1,6
**
** STEP 1 - INFLATION STEP (5 ATMOSPHERES)
**
*AMPLITUDE,NAME=RAMP,DEFINITION=SMOOTH STEP
0.,0.,  1.E-2,1.0
*AMPLITUDE,NAME=RAMPSMOOTH,DEFINITION=SMOOTH STEP
0.,0.,  2.E-2,1.0
*SURFACE,TYPE=ELEMENT, NAME = ASURF
ASURF,SNEG
*SURFACE,TYPE=ELEMENT, NAME = RSURF
RIGELEM,SPOS
*RIGID BODY,REFNODE=1,PINNSET=N1M1
*RIGID BODY,REFNODE=11000,ELSET=RIGELEM,TIENSET=RIG
*NSET, NSET=QA_TEST
2216,2218,2220,2222,
*STEP
*DYNAMIC,EXPLICIT
,1.E-2
*BOUNDARY,AMPLITUDE=RAMP
CAVITY,8,8, 5.066E5
**
** CONTACT DEFINITION
**
*NSET,NSET=RIG
N40,N41
*ELSET,ELSET=ASURF,GEN
26,2126,100
27,2127,100
28,2128,100
29,2129,100
30,2130,100
31,2131,100
32,2132,100
33,2133,100
34,2134,100
35,2135,100
36,2136,100
37,2137,100
38,2138,100
*CONTACT PAIR, INTERACTION = SMOOTH
ASURF, RSURF
*SURFACE INTERACTION, NAME = SMOOTH
**
** OUTPUT REQUESTS
**
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=2
*NODE OUTPUT, NSET=QA_TEST
U,
*OUTPUT, HIST
*ENERGY OUTPUT, VAR=PRESELECT
*NODE OUTPUT, NSET=RIGID
U, RF
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=4
*NODE OUTPUT
U,
*ELSET,ELSET=TOP
39,2139,1039
*END STEP
**
** STEP 2 - SEAL CAVITY, LOAD FURTHER
**
*STEP
*DYNAMIC,EXPLICIT
,2.E-2
*BOUNDARY, OP=NEW,AMPLITUDE=RAMPSMOOTH
1,1,6
N4,1,6
N6,2,3
YPLN1, 2
YPLN1, 4
YPLN1, 6
YPLN2, 2
YPLN2, 4
YPLN2, 6
RIGID,1,2
RIGID,4,6
RIGID,3,3, 2.E-2
*END STEP