*HEADING
*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 3: CPS4
*******************
**
*NODE,NSET=NALL
301,  7.,  0.
302,  8.,  0.
303,  8.,  1.
304,  7.,  1.
311,  7.,  2.
312,  8.,  2.
313,  8.,  3.
314,  7.,  3.
*NSET,NSET=LEFT
301,304,311,314
*NSET,NSET=BOT
301,302,311,312
*NSET,NSET=RIGHT
302,303,312,313
*ELEMENT,TYPE=CPS4,ELSET=TEST3
301,301,302,303,304
311,311,312,313,314
*ELSET,ELSET=TEST
311, 
*ELSET,ELSET=REFER
301, 
*SOLID SECTION,ELSET=TEST3,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 5: M3D4
******************
**
*NODE,NSET=NALL
501, 13.,  0.
502, 14.,  0.
503, 14.,  1.
504, 13.,  1.
511, 13.,  2.
512, 14.,  2.
513, 14.,  3.
514, 13.,  3.
*NSET,NSET=LEFT
501,504,511,514
*NSET,NSET=BOT
501,502,511,512
*NSET,NSET=RIGHT
502,503,512,513
*ELEMENT,TYPE=M3D4,ELSET=TEST5
501,501,502,503,504
511,511,512,513,514
*ELSET,ELSET=TEST
511, 
*ELSET,ELSET=REFER
501, 
*MEMBRANE SECTION,ELSET=TEST5,MATERIAL=ELAS
.25, 
******************
** TEST 6: S4R
******************
**
*NODE,NSET=NALL
601, 16.,  0.
602, 17.,  0.
603, 17.,  1.
604, 16.,  1.
611, 16.,  2.
612, 17.,  2.
613, 17.,  3.
614, 16.,  3.
*NSET,NSET=LEFT
601,604,611,614
*NSET,NSET=BOT
601,602,611,612
*NSET,NSET=RIGHT
602,603,612,613
*ELEMENT,TYPE=S4R,ELSET=TEST6
601,601,602,603,604
611,611,612,613,614
*ELSET,ELSET=TEST
611, 
*ELSET,ELSET=REFER
601, 
*SHELL SECTION,ELSET=TEST6,MATERIAL=ELAS
.25, 
*******************
** TEST 7: CAXA
******************
**
*NODE,NSET=NALL
701, 19.,  0.
702, 20.,  0.
703, 20.,  1.
704, 19.,  1.
711, 19.,  2.
712, 20.,  2.
713, 20.,  3.
714, 19.,  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=CAXA41,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 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=SFMAX1,ELSET=SFMAX1-1
1811,1814,1815
*SURFACE SECTION,ELSET=SFMAX1-1
*REBAR LAYER,GEOMETRY=ANGULAR
REB8A,1.000000,19.188136,0.000000,ELAS,10.000000,1
*EMBEDDED ELEMENT
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,ELAS,10.000000,1
*EMBEDDED ELEMENT
1812,
*NSET,NSET=LEFT
801,804,811,814
*NSET,NSET=BOT
801,802,811,812
*NSET,NSET=RIGHT
802,803,812,813
*ELEMENT,TYPE=CAX4,ELSET=TEST8
801, 801,802,803,804
811, 811,812,813,814
*INITIAL CONDITIONS,TYPE=STRESS,REBAR
1811, REB8A, 60.
1812, REB8B, 60.
*ELSET,ELSET=TEST
1812,811
*ELSET,ELSET=REFER
1811,801
*SOLID SECTION,ELSET=TEST8,MATERIAL=ELAS
*****************
** TEST 9: CPE4
*******************
**
*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=CPE4,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
*****************
** TEST 10: CPS4
*******************
**
*NODE,NSET=NALL
1001,  7.,  5.
1002,  8.,  5.
1003,  8.,  6.
1004,  7.,  6.
1011,  7.,  7.
1012,  8.,  7.
1013,  8.,  8.
1014,  7.,  8.
*NSET,NSET=LEFT
1001,1004,1011,1014
*NSET,NSET=BOT
1001,1002,1011,1012
*NSET,NSET=RIGHT
1002,1003,1012,1013
*ELEMENT,TYPE=CPE4,ELSET=TEST10
1001,1001,1002,1003,1004
1011,1011,1012,1013,1014
*ELSET,ELSET=TEST
1011, 
*ELSET,ELSET=REFER
1001, 
*SOLID SECTION,ELSET=TEST10,MATERIAL=ELAS
******************
** TEST 11: CGAX4
******************
**
*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=CGAX4,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
***************************
** TEST 12: M3D4 WITH REBAR
***************************
**
*NODE,NSET=NALL
1201, 13.,  5.
1202, 14.,  5.
1203, 14.,  6.
1204, 13.,  6.
1211, 13.,  7.
1212, 14.,  7.
1213, 14.,  8.
1214, 13.,  8.
*NODAL THICKNESS
1201, .125
1202, .175
1203, .175
1204, .125
1211, .125
1212, .175
1213, .175
1214, .125
*NSET,NSET=LEFT
1201,1204,1211,1214
*NSET,NSET=BOT
1201,1202,1211,1212
*NSET,NSET=RIGHT
1202,1203,1212,1213
*ELEMENT,TYPE=M3D4,ELSET=TEST12
1201,1201,1202,1203,1204
1211,1211,1212,1213,1214
*INITIAL CONDITIONS,TYPE=STRESS,REBAR
1201, REB12, 60.
1211, REB12, 60.
*ELSET,ELSET=TEST
1211, 
*ELSET,ELSET=REFER
1201, 
*MEMBRANE SECTION,ELSET=TEST12,MATERIAL=ELAS,NODAL THICKNESS
.25,
*REBAR LAYER
REB12,1,0.5,0,ELAS,30,1  
**************************
** TEST 13: S4R WITH REBAR 
**************************
**
*NODE,NSET=NALL
1301, 16.,  5.
1302, 17.,  5.
1303, 17.,  6.
1304, 16.,  6.
1311, 16.,  7.
1312, 17.,  7.
1313, 17.,  8.
1314, 16.,  8.
11314, 16.3,           5.
11315, 17.,          5.2
11316, 16.3,           7.
11317, 17.,          7.2
**ELEMENT,TYPE=SFMAX1,ELSET=SFMAX1-1
**11311,11314,11315
**SURFACE SECTION,ELSET=SFMAX1-1
**REBAR LAYER,GEOMETRY=ANGULAR
**REB9,1.000000,19.188136,0.000000,ELAS,10.000000,1
**EMBEDDED ELEMENT
**11311,
**ELEMENT,TYPE=SFMAX1,ELSET=SFMAX1-1
**11312,11316,11317
**SURFACE SECTION,ELSET=SFMAX1-1
**REBAR LAYER,GEOMETRY=ANGULAR
**REB8,1.000000,19.188136,0.000000,ELAS,10.000000,1
***EMBEDDED ELEMENT
**1812,
*NSET,NSET=LEFT
1301,1304,1311,1314
*NSET,NSET=BOT
1301,1302,1311,1312
*NSET,NSET=RIGHT
1302,1303,1312,1313
*ELEMENT,TYPE=S4R,ELSET=TEST13
1301,1301,1302,1303,1304
1311,1311,1312,1313,1314
*INITIAL CONDITIONS,TYPE=STRESS,REBAR
1301, REB13, 60.
1311, REB13, 60.
*ELSET,ELSET=TEST
1301, 
*ELSET,ELSET=REFER
1311, 
*SHELL SECTION,ELSET=TEST13,MATERIAL=ELAS
.25,
*REBAR LAYER
REB13,1,0.5,0,ELAS,30,1  
*******************
** TEST 14: CAXA
******************
**
*NODE,NSET=NALL
1401, 19.,  5.
1402, 20.,  5.
1403, 20.,  6.
1404, 19.,  6.
1411, 19.,  7.
1412, 20.,  7.
1413, 20.,  8.
1414, 19.,  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=CAXA41,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 15: S4 WITH REBAR 
**************************
**
*NODE,NSET=NALL
1501, 16.,  5.
1502, 17.,  5.
1503, 17.,  6.
1504, 16.,  6.
1511, 16.,  7.
1512, 17.,  7.
1513, 17.,  8.
1514, 16.,  8.
*NSET,NSET=LEFT
1501,1504,1511,1514
*NSET,NSET=BOT
1501,1502,1511,1512
*NSET,NSET=RIGHT
1502,1503,1512,1513
*ELEMENT,TYPE=S4,ELSET=TEST15
1501,1501,1502,1503,1504
1511,1511,1512,1513,1514
*INITIAL CONDITIONS,TYPE=STRESS,REBAR
1501, REB13, 60.
1511, REB13, 60.
*ELSET,ELSET=TEST
1511, 
*ELSET,ELSET=REFER
1501, 
*SHELL SECTION,ELSET=TEST15,MATERIAL=ELAS
.25,
*REBAR LAYER
REB15,1,0.5,0,ELAS,30,1  
**
**
**
** REMAINDER OF DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*INITIAL CONDITIONS,TYPE=STRESS
EALL,-100., .0, 0., 
*MATERIAL,NAME=ELAS
*ELASTIC
2000., .3
*DENSITY
1.E-4, 
**
**
*STEP,,NLGEOM
 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, 
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT
S,E
*ELEMENT OUTPUT,ELSET=TEST12,REBAR
S,E
*NODE OUTPUT
RF   , 
*NODE PRINT,FREQ=99
RF   , 
*EL FILE,REBAR
S,
E
*END STEP
**
**
*STEP,NLGEOM
 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
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT
S,E
*ELEMENT OUTPUT,ELSET=TEST12,REBAR
S,E
*END STEP
**
**
*STEP,NLGEOM
 STEP #3: STRETCH ELEMENT
*STATIC
1.,1.
*BOUNDARY
RIGHT,1,1,.1
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT
S,E
*ELEMENT OUTPUT,ELSET=TEST12,REBAR
S,E
*END STEP