*HEADING
 POROUS PLASTICITY WITH GURSON'S MODEL
*RESTART,TIMEMARKS=NO,WRITE,NUMBER INTERVAL=20
*NODE,NSET=N01
 1, 0. , 0.
 2, 1. , 0.
 3, 1. , 1.
 4, 0. , 1.
 5, 0. , 0. , 1.0
 6, 1. , 0. , 1.0
 7, 1. , 1. , 1.0
 8, 0. , 1. , 1.0
*NCOPY,OLD SET=N01,CHANGE NUMBER=10,SHIFT,NEW SET=NR1
 1.5, 0.0, 0.0
 0.0,
*NCOPY,OLD SET=N01,CHANGE NUMBER=20,SHIFT,NEW SET=NR1
 3.0, 0.0, 0.0
 0.0,
*NCOPY,OLD SET=N01,CHANGE NUMBER=30,SHIFT,NEW SET=NR1
 4.5, 0.0, 0.0
 0.0,
*NSET,NSET=NR1
 N01,
*NCOPY,OLD SET=NR1,CHANGE NUMBER=40,SHIFT,NEW SET=NR2
 0.0, 1.5, 0.0
 0.0,
*NCOPY,OLD SET=NR2,CHANGE NUMBER=40,SHIFT,NEW SET=NR3
 0.0, 1.5, 0.0
 0.0,
*NCOPY,OLD SET=NR3,CHANGE NUMBER=40,SHIFT,NEW SET=NR4
 0.0, 1.5, 0.0
 0.0,
*NSET,NSET=R1TOP
  3,  4,  7,  8
 13, 14, 17, 18
 23, 24, 25, 28
 33, 34, 35, 38
*NSET,NSET=R2TOP
 43, 44, 47, 48
 53, 54, 57, 58
 63, 64, 65, 68
 73, 74, 75, 78
*NSET,NSET=R3TOP
 83, 84, 87, 88
 93, 94, 97, 98
103,104,105,108
113,114,115,118
*NSET,NSET=R1BOT
  1,  2,  5,  6
 11, 12, 15, 16
 21, 22, 25, 26
 31, 32, 35, 36
*NSET,NSET=R2BOT
 41, 42, 45, 46
 51, 52, 55, 56
 61, 62, 65, 66
 71, 72, 75, 76
*NSET,NSET=R3BOT
 81, 82, 85, 86
 91, 92, 95, 96
101,102,105,106
111,112,115,116
*NSET,NSET=R4BOT1
121,122,125,126
141,142,145,146
*NSET,NSET=R4BOT2
131,132,135,136
151,152,155,156
*NSET,NSET=R2LHS,GENERATE 
 41, 71, 10
 44, 74, 10
 45, 75, 10
 48, 78, 10
*NSET,NSET=R34LHS,GENERATE
 81, 151, 10
 84, 154, 10
 85, 155, 10
 88, 158, 10
*NSET,NSET=R2RHS,GENERATE
 42, 72, 10
 43, 73, 10
 47, 77, 10
 46, 76, 10
*NSET,NSET=R34RHS,GENERATE
 82, 152, 10
 83, 153, 10
 87, 157, 10
 86, 156, 10
*NSET,NSET=BOTTOM,GENERATE
 1, 151, 10
 2, 152, 10
 5, 155, 10
 6, 156, 10
*NSET,NSET=BACK,GENERATE
 1, 121, 40
 2, 122, 40
 3, 123, 40
 4, 124, 40
*NSET,NSET=FRONT,GENERATE
 5, 125, 40
 6, 126, 40
 7, 127, 40
 8, 128, 40
*ELEMENT,TYPE=C3D8R,ELSET=C3D
  1,    1,   2,   3,   4,   5,   6,   7,   8
  2,   11,  12,  13,  14,  15,  16,  17,  18
  5,   41,  42,  43,  44,  45,  46,  47,  48
  6,   51,  52,  53,  54,  55,  56,  57,  58
  9,   81,  82,  83,  84,  85,  86,  87,  88
 10,   91,  92,  93,  94,  95,  96,  97,  98
 13,  121, 122, 123, 124, 125, 126, 127, 128
 14,  131, 132, 133, 134, 135, 136, 137, 138
*ELEMENT,TYPE=CPE4R,ELSET=CPE
  3,   21,  22,  23,  24
  7,   61,  62,  63,  64
 11,  101, 102, 103, 104
 15,  141, 142, 143, 144
*ELEMENT,TYPE=CPS4R,ELSET=CPS
  4,   31,  32,  33,  34
  8,   71,  72,  73,  74
 12,  111, 112, 113, 114
 16,  151, 152, 153, 154
*ELSET,ELSET=EALL,GEN
 1, 16, 1
*ELSET,ELSET=R123,GEN
 1, 12, 1
*ELSET,ELSET=R4,GEN
 13, 16, 1
***
*ELSET,ELSET=SOLID2D_R123
3,4,7,8,11,12
*ELSET,ELSET=SOLID3D_R123
1,2,5,6,9,10
*ELSET,ELSET=SOLID2D_R4
15,16
*ELSET,ELSET=SOLID3D_R4
13,14
***
*SOLID SECTION,MATERIAL=GURSON1,ELSET=SOLID2D_R123
*SOLID SECTION,MATERIAL=GURSON1,ELSET=SOLID3D_R123
*SOLID SECTION,MATERIAL=GURSON2,ELSET=SOLID2D_R4
*SOLID SECTION,MATERIAL=GURSON2,ELSET=SOLID3D_R4
***
*MATERIAL,NAME=GURSON1
*DENSITY
 0.001, 
*ELASTIC
 300., 0.3
*PLASTIC
 1.0, 
*POROUS METAL PLASTICITY, RELATIVE DENSITY=0.95
1. , 1. , 1.
*MATERIAL,NAME=GURSON2
*DENSITY
 0.001, 
*ELASTIC
 300., 0.3
*PLASTIC
 1.00, 0.0
 2.00, 0.5
*POROUS METAL PLASTICITY, RELATIVE DENSITY=0.95
1.5 , 1. , 2.25
*VOID NUCLEATION
0.3 , 0.1 , 0.04
*BOUNDARY
 R1BOT, 2
 R1TOP, 2
 FRONT, 3
 BACK,  3
*AMPLITUDE,NAME=RAMP
 0.0, 0.0, 1.0, 1.0
**
**
*ELSET,ELSET=QA_TEST_ALLELEMS
R123, R4
**
*NSET, NSET=QA_TEST_ALLNODES
NR1, R1TOP, R2TOP, R3TOP
R1BOT, R2BOT, R3BOT, R4BOT1
R4BOT2, R2LHS, R34LHS, R2RHS
R34RHS, BOTTOM, BACK, FRONT
**
*STEP
*DYNAMIC,EXPLICIT
 ,1.0 
*BOUNDARY,TYPE=VEL,AMPLITUDE=RAMP
 R1BOT ,1,,-0.05
 R1TOP ,1,, 0.05
 R2BOT ,2,, 0.05 
 R2TOP ,2,,-0.05 
 R2LHS ,1,, 0.05
 R2RHS ,1,,-0.05
 R3BOT ,2,,-0.05
 R4BOT1,2,, -0.05 
 R4BOT2,2,, 0.05
 R3TOP ,2,, 0.05
 R34LHS,1,,-0.05
 R34RHS,1,, 0.05
**
*FILE OUTPUT,NUMBER INTERVAL=3,TIME MARKS=YES
*EL FILE
 MISES, 
 PRESS, 
 PEEQ, 
 PE,
 ELEN,
**
*OUTPUT, FIELD, NUMBER INTERVAL=3, TIMEMARKS=YES
*ELEMENT OUTPUT, ELSET=QA_TEST_ALLELEMS
 MISES, 
 PRESS, 
 PEEQ, 
 PE,
 ELEN,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.333333
***OUTPUT,HISTORY,VARIABLE=PRESELECT
*END STEP