*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: 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 4: 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 5: 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,ELAS1,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,ELAS1,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, UNBALANCED STRESS=STEP
1811, REB8A, 60.
1812, REB8B, 60.
*ELSET,ELSET=TEST
811,1812 
*ELSET,ELSET=REFER
801,1811 
*SOLID SECTION,ELSET=TEST8,MATERIAL=ELAS
*****************
** TEST 6: 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 7: 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 8: CAXA41
******************
**
*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, UNBALANCED STRESS=STEP
EALL,-10., .0, 0.,
*MATERIAL,NAME=ELAS1
*ELASTIC
3000.,0.3
*DENSITY
1.E-4, 
*PLASTIC
10.0,0.0
17.41,0.05
27.27,0.50
29.99,0.80
*MATERIAL,NAME=ELAS
*ELASTIC
3000., .3
*DENSITY
1.E-4, 
*PLASTIC
10.0,0.0
17.41,0.05
27.27,0.50
29.99,0.80
*POROUS METAL PLASTICITY, RELATIVE DENSITY=0.95
1.5,1.0,2.25
*VOID NUCLEATION
0.3,0.1,0.04
**
**
*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,
VVF, 
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT
S,E,PE
VVF, 
*NODE OUTPUT
RF, 
*NODE PRINT,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: STRETCH ELEMENT
*STATIC
0.1,1.
*BOUNDARY
RIGHT,1,1,.1
*END STEP
*STEP,PERT
 STEP #4: STRETCH ELEMENT, LINEAR PERT.
*STATIC
0.1,1.
*LOAD CASE,name=stretch-1
*BOUNDARY
RIGHT,1,1,.1
*END LOAD CASE
*LOAD CASE,name=stretch-2
*BOUNDARY
RIGHT,1,1,.2
*END LOAD CASE
*LOAD CASE,name=stretch-3
*BOUNDARY
RIGHT,1,1,.3
*END LOAD CASE
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,E
*NODE OUTPUT
U,RF
*END STEP