*HEADING
*PREPRINT,MODEL=NO,HISTORY=NO
*****************
** TEST 1: CAX4H
*****************
*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=CAX4H,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: CPE4H
*******************
**
*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=CPE4H,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 3: CGAX4H
******************
**
*NODE,NSET=NALL
401, 7.,  0.
402, 8.,  0.
403, 8.,  1.
404, 7.,  1.
411, 7.,  2.
412, 8.,  2.
413, 8.,  3.
414, 7.,  3.
*NSET,NSET=LEFT
401,404,411,414
*NSET,NSET=BOT
401,402,411,412
*NSET,NSET=RIGHT
402,403,412,413
*ELEMENT,TYPE=CGAX4H,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 4: CAXA4H1
******************
**
*NODE,NSET=NALL
701, 10.,  0.
702, 11.,  0.
703, 11.,  1.
704, 10.,  1.
711, 10.,  2.
712, 11.,  2.
713, 11.,  3.
714, 10.,  3.
*NSET,NSET=LEFT
701,704,711,714
10701,10704,10711,10714
*NSET,NSET=BOT
701,702,711,712
10701,10702,10711,10712
*NSET,NSET=RIGHT
702,703,712,713
10702,10703,10712,10713
*ELEMENT,TYPE=CAXA4H1,ELSET=TEST7,OFFSET=10000
701,701,702,703,704
711,711,712,713,714
*ELSET,ELSET=TEST
711, 
*ELSET,ELSET=REFER
701, 
*SOLID SECTION,ELSET=TEST7,MATERIAL=ELAS
***************************
** TEST 5: CAX4H 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
*NSET,NSET=LEFT
801,804,811,814
*NSET,NSET=BOT
801,802,811,812
*NSET,NSET=RIGHT
802,803,812,813
*ELEMENT,TYPE=CAX4H,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,REBAR,10.000000,1
*EMBEDDED ELEMENT,HOST ELSET=CAX4H_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,REBAR,10.000000,1
*EMBEDDED ELEMENT,HOST ELSET=CAX4H_TWO
1812,
*INITIAL CONDITIONS,TYPE=STRESS,REBAR
1811, REB8A, -1.0
1812, REB8B, -1.0
*ELSET,ELSET=TEST
811,1812 
*ELSET,ELSET=REFER
801,1811 
*ELSET,ELSET=CAX4H_ONE
801
*ELSET,ELSET=CAX4H_TWO
811
*SOLID SECTION,ELSET=TEST8,MATERIAL=ELAS
*****************
** TEST 6: CPE4RH
*******************
**
*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=CPE4RH,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 7: CGAX4RH
******************
**
*NODE,NSET=NALL
1101, 7.,  5.
1102, 8.,  5.
1103, 8.,  6.
1104, 7.,  6.
1111, 7.,  7.
1112, 8.,  7.
1113, 8.,  8.
1114, 7.,  8.
*NSET,NSET=LEFT
1101,1104,1111,1114
*NSET,NSET=BOT
1101,1102,1111,1112
*NSET,NSET=RIGHT
1102,1103,1112,1113
*ELEMENT,TYPE=CGAX4RH,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, 
*******************
** TEST 8: CAXAH
******************
**
*NODE,NSET=NALL
1401, 10.,  5.
1402, 11.,  5.
1403, 11.,  6.
1404, 10.,  6.
1411, 10.,  7.
1412, 11.,  7.
1413, 11.,  8.
1414, 10.,  8.
*NSET,NSET=LEFT
1401,1404,1411,1414
11401,11404,11411,11414
*NSET,NSET=BOT
1401,1402,1411,1412
11401,11402,11411,11412
*NSET,NSET=RIGHT
1402,1403,1412,1413
11402,11403,11412,11413
*ELEMENT,TYPE=CAXA4H1,ELSET=TEST14,OFFSET=10000
1401,1401,1402,1403,1404
1411,1411,1412,1413,1414
*ELSET,ELSET=TEST
1411, 
*ELSET,ELSET=REFER
1401, 
*SOLID SECTION,ELSET=TEST14,MATERIAL=ELAS
**************************
** TEST 9: S4R WITH REBAR
**************************
**
*NODE,NSET=NALL
1501, 0.,  10.
1502, 1.,  10.
1503, 1.,  11.
1504, 0.,  11.
1511, 0.,  12.
1512, 1.,  12.
1513, 1.,  13.
1514, 0.,  13.
*NSET,NSET=LEFT
1501,1504,1511,1514
*NSET,NSET=BOT
1501,1502,1511,1512
*NSET,NSET=RIGHT
1502,1503,1512,1513
*ELEMENT,TYPE=S4R,ELSET=TEST15
1501,1501,1502,1503,1504
1511,1511,1512,1513,1514
*ELSET,ELSET=TEST
1511, 
*ELSET,ELSET=REFER
1501, 
*SHELL SECTION,ELSET=TEST15,MATERIAL=ELAS
1., 
*REBAR LAYER
REB15,0.5,1.0,0,REBAR,30,1 
**************************
** TEST 10: M3D4R WITH REBAR
**************************
**
*NODE,NSET=NALL
1601, 4.,  10.
1602, 5.,  10.
1603, 5.,  11.
1604, 4.,  11.
1611, 4.,  12.
1612, 5.,  12.
1613, 5.,  13.
1614, 4.,  13.
*NSET,NSET=LEFT
1601,1604,1611,1614
*NSET,NSET=BOT
1601,1602,1611,1612
*NSET,NSET=RIGHT
1602,1603,1612,1613
*ELEMENT,TYPE=M3D4R,ELSET=TEST16
1601,1601,1602,1603,1604
1611,1611,1612,1613,1614
*ELSET,ELSET=TEST
1611, 
*ELSET,ELSET=REFER
1601, 
*MEMBRANE SECTION,ELSET=TEST16,MATERIAL=ELAS
1.,
*REBAR LAYER
REB16,0.5,1.0,0,REBAR,30,1  
**
** REMAINDER OF DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*MATERIAL,NAME=REBAR
*ELASTIC
1000., .25
*MATERIAL,NAME=ELAS
*HYPERELASTIC
80., 20.
*INITIAL CONDITIONS,TYPE=STRESS
 EALL,1.0
**
**
*STEP,,NLGEOM, INC=100
 step #1: REMOVE TEST ELEMENTS
*STATIC
0.1, 1.0
*BOUNDARY
NALL,3,6
LEFT,1
BOT,2
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*EL PRINT,POSITION=CENTROID,FREQ=99
1,
S,
STH,
E,
*EL PRINT,REBAR,FREQ=99
S,E,RBANG,RBROT
*EL FILE,REBAR,FREQ=99
S,E,RBANG
*EL FILE,REBAR,FREQ=99
RBROT,
*EL FILE,POSITION=CENTROID,FREQ=99
S,E
*NODE PRINT,NSET=LEFT,FREQ=99
RF   , 
*NODE FILE,FREQ=99
RF   , 
*END STEP
**
**
*STEP,NLGEOM, INC=100
 STEP #2: ADD TEST ELEMENT BACK IN
*STATIC
.1,1.
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST, 
*BOUNDARY
RIGHT,1,1,.2
*END STEP
**
**
*STEP,NLGEOM, INC=100
 STEP #3: REMOVE TEST ELEMENT
*STATIC
.1,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
NALL,3,6
LEFT,1
BOT,2
RIGHT,1
*EL FILE,REBAR,FREQ=99
S,E,RBANG
*EL FILE,POSITION=CENTROID,FREQ=99
S,E
*END STEP
**
**
*STEP,NLGEOM, INC=100
 STEP #4: STRETCH ELEMENTS FURTHER WHILE ELEMENTS ARE ABSENT
*STATIC
.1,1.
*BOUNDARY
RIGHT,1,1,.3
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST, 
*EL FILE,REBAR,FREQ=99
S,E,RBANG
*EL FILE,REBAR,FREQ=99
RBROT,
*EL FILE,POSITION=CENTROID,FREQ=99
S,E
*END STEP
**