*HEADING
 KEYWORD-PARAMETER COMBINATION TESTED: CRUSHABLE FOAM
*PREPRINT,MODEL=NO,HISTORY=NO
*****************
** TEST 1: CAX4
*****************
*NODE,NSET=NALL
1,  1.,  0.
2,  2.,  0.
3,  2.,  1.
4,  1.,  1.
11, 1.,  2.
12, 2.,  2.
13, 2.,  3.
14, 1.,  3.
*NSET,NSET=LEFT
1,4,11,14
*NSET,NSET=BOT
1,2,11,12
*NSET,NSET=RIGHT
2,3,12,13
*ELEMENT,TYPE=CAX4,ELSET=TEST1
1,1,2,3,4
11,11,12,13,14
*ELSET,ELSET=TEST
11, 
*ELSET,ELSET=REFER
1, 
*SOLID SECTION,ELSET=TEST1,MATERIAL=ELAS
*****************
** TEST 2: CPE4
*******************
**
*NODE,NSET=NALL
101,  4.,  0.
102,  5.,  0.
103,  5.,  1.
104,  4.,  1.
111,  4.,  2.
112,  5.,  2.
113,  5.,  3.
114,  4.,  3.
*NSET,NSET=LEFT
101,104,111,114
*NSET,NSET=BOT
101,102,111,112
*NSET,NSET=RIGHT
102,103,112,113
*ELEMENT,TYPE=CPE4,ELSET=TEST2
101,101,102,103,104
111,111,112,113,114
*ELSET,ELSET=TEST
111, 
*ELSET,ELSET=REFER
101, 
*SOLID SECTION,ELSET=TEST2,MATERIAL=ELAS
******************
** TEST 4: CGAX4
******************
**
*NODE,NSET=NALL
401, 10.,  0.
402, 11.,  0.
403, 11.,  1.
404, 10.,  1.
411, 10.,  2.
412, 11.,  2.
413, 11.,  3.
414, 10.,  3.
*NSET,NSET=LEFT
401,404,411,414
*NSET,NSET=BOT
401,402,411,412
*NSET,NSET=RIGHT
402,403,412,413
*ELEMENT,TYPE=CGAX4,ELSET=TEST4
401,401,402,403,404
411,411,412,413,414
*ELSET,ELSET=TEST
411, 
*ELSET,ELSET=REFER
401, 
*SOLID SECTION,ELSET=TEST4,MATERIAL=ELAS
***************************
** TEST 8: CAX4 WITH REBAR
***************************
*NODE,NSET=NALL
801,  1.,  5.
802,  2.,  5.
803,  2.,  6.
804,  1.,  6.
811,  1.,  7.
812,  2.,  7.
813,  2.,  8.
814,  1.,  8.
1814,1.3,5.
1815,2.,5.2
1816,1.3,7.
1817,2.,7.2
*ELEMENT,TYPE=CAX4,ELSET=TEST8
801, 801,802,803,804
811, 811,812,813,814
**
**
*ELEMENT,TYPE=SFMAX1,ELSET=SFMAX1-1
1811,1814,1815
*SURFACE SECTION,ELSET=SFMAX1-1
*REBAR LAYER,GEOMETRY=ANGULAR
REB8A,1.000000,19.188136,0.000000,ELASLIN,10.000000,1
*EMBEDDED ELEMENT,HOST ELSET=CAX4_ONE
1811,
*ELEMENT,TYPE=SFMAX1,ELSET=SFMAX1-2
1812,1816,1817
*SURFACE SECTION,ELSET=SFMAX1-2
*REBAR LAYER,GEOMETRY=ANGULAR
REB8B,1.000000,19.188136,0.000000,ELASLIN,10.000000,1
*EMBEDDED ELEMENT,HOST ELSET=CAX4_TWO
1812,
**
**
*NSET,NSET=LEFT
801,804,811,814
*NSET,NSET=BOT
801,802,811,812
*NSET,NSET=RIGHT
802,803,812,813
*INITIAL CONDITIONS,TYPE=STRESS,REBAR
1811, REB8A, -1.0
1812, REB8B, -1.0
*ELSET,ELSET=CAX4_ONE
801
*ELSET,ELSET=CAX4_TWO
811
*ELSET,ELSET=TEST
811, 1812
*ELSET,ELSET=REFER
801, 1811
*SOLID SECTION,ELSET=TEST8,MATERIAL=ELAS
*****************
** TEST 9: CPE4R
*******************
**
*NODE,NSET=NALL
901,  4.,  5.
902,  5.,  5.
903,  5.,  6.
904,  4.,  6.
911,  4.,  7.
912,  5.,  7.
913,  5.,  8.
914,  4.,  8.
*NSET,NSET=LEFT
901,904,911,914
*NSET,NSET=BOT
901,902,911,912
*NSET,NSET=RIGHT
902,903,912,913
*ELEMENT,TYPE=CPE4R,ELSET=TEST9
901,901,902,903,904
911,911,912,913,914
*ELSET,ELSET=TEST
911, 
*ELSET,ELSET=REFER
901, 
*SOLID SECTION,ELSET=TEST9,MATERIAL=ELAS
*HOURGLASS STIFFNESS
48.0, 
******************
** TEST 11: CGAX4R
******************
**
*NODE,NSET=NALL
1101, 10.,  5.
1102, 11.,  5.
1103, 11.,  6.
1104, 10.,  6.
1111, 10.,  7.
1112, 11.,  7.
1113, 11.,  8.
1114, 10.,  8.
*NSET,NSET=LEFT
1101,1104,1111,1114
*NSET,NSET=BOT
1101,1102,1111,1112
*NSET,NSET=RIGHT
1102,1103,1112,1113
*ELEMENT,TYPE=CGAX4R,ELSET=TEST11
1101,1101,1102,1103,1104
1111,1111,1112,1113,1114
*ELSET,ELSET=TEST
1111, 
*ELSET,ELSET=REFER
1101, 
*SOLID SECTION,ELSET=TEST11,MATERIAL=ELAS
*HOUR GLASS STIFFNESS
48.0, 
**
** REMAINDER OF DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
 TEST, REFER
*MATERIAL,NAME=ELAS
*ELASTIC
 3.E6, 0.
*CRUSHABLE FOAM
 1.1, 0.1
*CRUSHABLE FOAM HARDENING
 220000.0,  0.0
 246512.1,  0.1
 272938.0,  0.2
 299019.2,  0.3
 324549.1,  0.4
 349345.3,  0.5
 373259.8,  0.6
 396170.9,  0.7
 418011.1,  0.8
 438721.0,  0.9
 458268.2,  1.0
 493843.0,  1.2
 524836.3,  1.4
 551528.4,  1.6
 574307.3,  1.8
 593594.8,  2.0
 629359.6,  2.5
 651994.7,  3.0
 683338.7,  5.0
 688332.2, 10.0
**
*MATERIAL,NAME=ELASLIN
*ELASTIC
 0.6666667E6, 0.0
**
*INITIAL CONDITIONS,TYPE=STRESS
 EALL,-1.0
**
**
*STEP,,NLGEOM, UNSYMM=YES, INC=100
 step #1: REMOVE TEST ELEMENT
*STATIC
 0.1, 1.0
*BOUNDARY
NALL,3,6
LEFT,1
BOT,2
RIGHT,1
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*EL PRINT,POSITION=CENTROID,FREQ=99
1,
S,
E,
PE,
*EL PRINT, REBAR
S,
E,
*EL FILE,POSITION=CENTROID,FREQ=99
S,E,PE
*EL FILE, REBAR
S,
E,
*NODE PRINT,FREQ=99
RF   , 
*NODE FILE,FREQ=99
RF   , 
*END STEP
**
**
*STEP,NLGEOM,UNSYMM=YES, INC=100
 STEP #2: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
NALL,3,6
LEFT,1
BOT,2
RIGHT,1
*END STEP
**
**
*STEP,NLGEOM,UNSYMM=YES,INC=100
 STEP #3: COMPRESS ELEMENT
*STATIC
0.01,1.
*BOUNDARY
RIGHT,1,1,-.1
*END STEP