*HEADING
 AIRSPRING -- 180 DEG MODEL 
*RESTART,WRITE,NUMBER INTERVAL=10
*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
*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
*RIGID BODY,REFNODE=1,PINNSET=N1M1
*RIGID BODY,REFNODE=11000,ELSET=RIGELEM,
TIENSET=RIG
*NSET, NSET=QA_TEST_CAVITY
CAVITY,
*NSET, NSET=QA_TEST_RIGID
RIGID,
*************************************************************
**General contact definition
*Contact
*Contact Inclusions, All Element Based
*Contact Property Assignment
 , , SMOOTH
*SURFACE INTERACTION, NAME = SMOOTH
*************************************************************
*STEP
*DYNAMIC,EXPLICIT
,1.E-2
*BOUNDARY,AMPLITUDE=RAMP
CAVITY,8,8, 5.066E5
*NSET,NSET=RIG
N40,N41
** OUTPUT REQUESTS
**
*FILE OUTPUT,NUMBER INTERVAL=2
*NODE FILE, NSET=CAVITY
PCAV, CVOL
*NODE FILE, NSET=RIGID
U, 
RF,
*****
*OUTPUT, HISTORY,TIME INTERVAL=1.E-5
*NODE OUTPUT, NSET=QA_TEST_CAVITY
PCAV, CVOL
*OUTPUT, HIST, FREQ=9999
*ENERGY OUTPUT, VAR=PRESELECT
****** 
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=4
*NODE OUTPUT
U,
***OUTPUT,HISTORY,TIME INTERVAL=1.E-5
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLAE,ALLSE,ETOTAL,DT,DTFB,ALLVD
*ELSET,ELSET=TOP
39,2139,1039
***ELEMENT OUTPUT,ELSET=TOP
**S,E
***ELEMENT OUTPUT,ELSET=TOP,REBAR=PLSREBAR
**S,E,RBANG,RBROT
***NODE OUTPUT,NSET=CAVITY
**PCAV,CVOL
***NODE OUTPUT,NSET=RIG
**U,RF
***NODE OUTPUT,NSET=RIGID
**U,RF
*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