*HEADING
 AIRSPRING -- 180 DEG MODEL 
*Preprint, Contact=Yes, Model=Yes, History=Yes
*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, 
*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
*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, FIELD, TIME MARKS=YES, NUMBER INTER=2
*NODE OUTPUT, NSET=QA_TEST_RIGID
U,
*OUTPUT, HIST, TIME INTERVAL=1e-3
*NODE OUTPUT, NSET=QA_TEST_CAVITY
PCAV, CVOL
*ENERGY OUTPUT, VAR=PRESELECT
****** 
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=10
*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