*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, 
*ELEMENT, TYPE = SFM3D3, ELSET = FLUID
6001,  1,  4, 104
6040, 40, 41, 140
*ELGEN, ELSET = FLUID
6001, 22, 100, 100
6040, 22, 100, 100
*SURFACE SECTION, ELSET=FLUID
*ELSET, ELSET=SPRINGSURF
    2,    4,    6,    8,   10,   12,   14,   16,   18,   20,   22,   26,   27,   28,   29,   30
   31,   32,   33,   34,   35,   36,   37,   38,   39,  102,  104,  106,  108,  110,  112,  114
  116,  118,  120,  122,  126,  127,  128,  129,  130,  131,  132,  133,  134,  135,  136,  137
  138,  139,  202,  204,  206,  208,  210,  212,  214,  216,  218,  220,  222,  226,  227,  228
  229,  230,  231,  232,  233,  234,  235,  236,  237,  238,  239,  302,  304,  306,  308,  310
  312,  314,  316,  318,  320,  322,  326,  327,  328,  329,  330,  331,  332,  333,  334,  335
  336,  337,  338,  339,  402,  404,  406,  408,  410,  412,  414,  416,  418,  420,  422,  426
  427,  428,  429,  430,  431,  432,  433,  434,  435,  436,  437,  438,  439,  502,  504,  506
  508,  510,  512,  514,  516,  518,  520,  522,  526,  527,  528,  529,  530,  531,  532,  533
  534,  535,  536,  537,  538,  539,  602,  604,  606,  608,  610,  612,  614,  616,  618,  620
  622,  626,  627,  628,  629,  630,  631,  632,  633,  634,  635,  636,  637,  638,  639,  702
  704,  706,  708,  710,  712,  714,  716,  718,  720,  722,  726,  727,  728,  729,  730,  731
  732,  733,  734,  735,  736,  737,  738,  739,  802,  804,  806,  808,  810,  812,  814,  816
  818,  820,  822,  826,  827,  828,  829,  830,  831,  832,  833,  834,  835,  836,  837,  838
  839,  902,  904,  906,  908,  910,  912,  914,  916,  918,  920,  922,  926,  927,  928,  929
  930,  931,  932,  933,  934,  935,  936,  937,  938,  939, 1002, 1004, 1006, 1008, 1010, 1012
 1014, 1016, 1018, 1020, 1022, 1026, 1027, 1028, 1029, 1030, 1031, 1032, 1033, 1034, 1035, 1036
 1037, 1038, 1039, 1102, 1104, 1106, 1108, 1110, 1112, 1114, 1116, 1118, 1120, 1122, 1126, 1127
 1128, 1129, 1130, 1131, 1132, 1133, 1134, 1135, 1136, 1137, 1138, 1139, 1202, 1204, 1206, 1208
 1210, 1212, 1214, 1216, 1218, 1220, 1222, 1226, 1227, 1228, 1229, 1230, 1231, 1232, 1233, 1234
 1235, 1236, 1237, 1238, 1239, 1302, 1304, 1306, 1308, 1310, 1312, 1314, 1316, 1318, 1320, 1322
 1326, 1327, 1328, 1329, 1330, 1331, 1332, 1333, 1334, 1335, 1336, 1337, 1338, 1339, 1402, 1404
 1406, 1408, 1410, 1412, 1414, 1416, 1418, 1420, 1422, 1426, 1427, 1428, 1429, 1430, 1431, 1432
 1433, 1434, 1435, 1436, 1437, 1438, 1439, 1502, 1504, 1506, 1508, 1510, 1512, 1514, 1516, 1518
 1520, 1522, 1526, 1527, 1528, 1529, 1530, 1531, 1532, 1533, 1534, 1535, 1536, 1537, 1538, 1539
 1602, 1604, 1606, 1608, 1610, 1612, 1614, 1616, 1618, 1620, 1622, 1626, 1627, 1628, 1629, 1630
 1631, 1632, 1633, 1634, 1635, 1636, 1637, 1638, 1639, 1702, 1704, 1706, 1708, 1710, 1712, 1714
 1716, 1718, 1720, 1722, 1726, 1727, 1728, 1729, 1730, 1731, 1732, 1733, 1734, 1735, 1736, 1737
 1738, 1739, 1802, 1804, 1806, 1808, 1810, 1812, 1814, 1816, 1818, 1820, 1822, 1826, 1827, 1828
 1829, 1830, 1831, 1832, 1833, 1834, 1835, 1836, 1837, 1838, 1839, 1902, 1904, 1906, 1908, 1910
 1912, 1914, 1916, 1918, 1920, 1922, 1926, 1927, 1928, 1929, 1930, 1931, 1932, 1933, 1934, 1935
 1936, 1937, 1938, 1939, 2002, 2004, 2006, 2008, 2010, 2012, 2014, 2016, 2018, 2020, 2022, 2026
 2027, 2028, 2029, 2030, 2031, 2032, 2033, 2034, 2035, 2036, 2037, 2038, 2039, 2102, 2104, 2106
 2108, 2110, 2112, 2114, 2116, 2118, 2120, 2122, 2126, 2127, 2128, 2129, 2130, 2131, 2132, 2133
 2134, 2135, 2136, 2137, 2138, 2139, 6001, 6040, 6101, 6140, 6201, 6240, 6301, 6340, 6401, 6440
 6501, 6540, 6601, 6640, 6701, 6740, 6801, 6840, 6901, 6940, 7001, 7040, 7101, 7140, 7201, 7240
 7301, 7340, 7401, 7440, 7501, 7540, 7601, 7640, 7701, 7740, 7801, 7840, 7901, 7940, 8001, 8040
 8101, 8140
*SURFACE, TYPE=ELEMENT, NAME=CAVSURF
SPRINGSURF, SPOS
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16, UNIVERSAL GAS CONSTANT=8.314
**
**FLUID CAVITY DEFINITON
**
*FLUID CAVITY, NAME=SPRING, BEHAVIOR=GAS, REF NODE=50000, SURFACE=CAVSURF, AMBIENT PRESSURE=101.36E+3

*FLUID BEHAVIOR, NAME=GAS
*MOLECULAR WEIGHT
0.044
*CAPACITY, TYPE=TABULAR
30
**
** 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_CAVITY
CAVITY,
*NSET, NSET=QA_TEST_RIGID
RIGID,
*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_RIGID
UT,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=4
*NODE OUTPUT
U,
*OUTPUT, HIST
*NODE OUTPUT, NSET=QA_TEST_CAVITY
PCAV, CVOL
*ENERGY OUTPUT, VAR=PRESELECT
*NODE OUTPUT, NSET=RIGID
U, RF
***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