*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 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 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.
*ELSET,ELSET=TEST
1812,811 
*ELSET,ELSET=REFER
1811,801 
*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
******************
** 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
*******************
** 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
**
**
**
** REMAINDER OF DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*INITIAL CONDITIONS,TYPE=STRESS
EALL,-100., .0, 0., 
*MATERIAL,NAME=ELAS
*ELASTIC
12857,0.28
*CAP PLASTICITY
173.20508,30.0,0.61858957,0.027,0.01,1.0
*CAP HARDENING
** PB, EVOLP
 213.0, 0.00
 222.0, 0.01
 242.0, 0.02
 282.0, 0.03
 362.0, 0.04
 522.0, 0.05
 842.0, 0.06
1482.0, 0.07                                              
2762.0, 0.08
**
*DENSITY
1.5E-4,
**
**
*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,
PEQC,
 
*NODE PRINT,FREQ=99
RF   , 
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT
S,E,PE,PEQC,RD,VVF
*ELEMENT OUTPUT, REBAR
S,
E,
*NODE OUTPUT
RF
*EL FILE, REBAR
S,
E,
*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: STRETCH ELEMENT
*STATIC
0.1,1.
*BOUNDARY
RIGHT,1,1,.1
*END STEP