*HEADING
 ductile failure with Mises plasticity
*RESTART,WRITE,NUMBER INTERVAL=10
*PREPRINT,MODEL=YES
*NODE,NSET=NALL
1,0,0,0
2,0,1,0
3,1,0,0
4,1,1,0
5,2,0,0
6,2,1,0
7,3,0,0
8,3,1,0
9,4,0,0
10,4,1,0
11,5,0,0
12,6,0,0
13,6,1,0
14,5,1,0
15,5,0,1
16,6,0,1
17,6,1,1
18,5,1,1
19,9,0,0
20,10,0,0
21,10,1,0
22,9,1,0
23,9,0,0.1
24,10,0,0.1
25,10,1,0.1
26,9,1,0.1
27,13,0,0
28,14,0,0
29,14,1,0
30,13,1,0
31,17,0,0
32,18,0,0
33,18,1,0
34,17,1,0
35,21,0,0
36,22,0,0
37,22,1,0
38,21,1,0
39,25,0,0
40,26,0,0
41,26,1,0
42,25,1,0
43,29,0,0
44,30,0,0
45,30,1,0
46,29,1,0
47,33,0,0
48,34,0,0
49,34,1,0
50,33,1,0
51,37,0,0
52,38,0,0
53,38,1,0
54,37,1,0
61,41,0,0
62,42,0,0
63,42,1,0
64,41,1,0
65,41,0,1
66,42,0,1
67,42,1,1
68,41,1,1
77,45,0,0
78,46,0,0
79,46,1,0
80,45,1,0
**
101, 0,2,0
102, 0,3,0
103, 1,2,0
104, 1,3,0
105, 2,2,0
106, 2,3,0
107, 3,2,0
108, 3,3,0
109, 4,2,0
110, 4,3,0
111, 5,2,0
112, 6,2,0
113, 6,3,0
114, 5,3,0
115, 5,2,1
116, 6,2,1
117, 6,3,1
118, 5,3,1
119, 9,2,0
120,10,2,0
121,10,3,0
122, 9,3,0
123, 9,2,0.1
124,10,2,0.1
125,10,3,0.1
126, 9,3,0.1
127,13,2,0
128,14,2,0
129,14,3,0
130,13,3,0
131,17,2,0
132,18,2,0
133,18,3,0
134,17,3,0
135,21,2,0
136,22,2,0
137,22,3,0
138,21,3,0
139,25,2,0
140,26,2,0
141,26,3,0
142,25,3,0
143,29,2,0
144,30,2,0
145,30,3,0
146,29,3,0
147,33,2,0
148,34,2,0
149,34,3,0
150,33,3,0
151,37,2,0
152,38,2,0
153,38,3,0
154,37,3,0
161,41,2,0
162,42,2,0
163,42,3,0
164,41,3,0
165,41,2,1
166,42,2,1
167,42,3,1
168,41,3,1
177,45,2,0
178,46,2,0
179,46,3,0
180,45,3,0
********TRUSS ELEMENTS*******************
*ELEMENT,TYPE=T2D2,ELSET=TRUSS
1,1,2
*ELEMENT,TYPE=T3D2,ELSET=TRUSS
2,3,4
*ELEMENT,TYPE=T2D2,ELSET=TRUSS_REF
101,101,102
*ELEMENT,TYPE=T3D2,ELSET=TRUSS_REF
102,103,104
**
*NSET,NSET=TRUSS_BOTTOM
1,3,
101,103
*NSET,NSET=TRUSS_TOP
2,4,
102,104
*NSET,NSET=TRUSS_ALL
TRUSS_BOTTOM,TRUSS_TOP
*******BEAM AND AXIAL SHELL**************
*ELEMENT,TYPE=B31,ELSET=BEAM
3,5,6
*ELEMENT,TYPE=PIPE31,ELSET=BEAM
4,7,8
*ELEMENT,TYPE=SAX1,ELSET=SHELL
5,9,10
*ELEMENT,TYPE=B31,ELSET=BEAM_REF
103,105,106
*ELEMENT,TYPE=PIPE31,ELSET=BEAM_REF
104,107,108
*ELEMENT,TYPE=SAX1,ELSET=SHELL_REF
105,109,110
**
*NSET,NSET=BEAM_BOTTOM
5,7,9,
105,107,109
*NSET,NSET=BEAM_TOP
6,8,10,
106,108,110
*NSET,NSET=BEAM_ALL
BEAM_BOTTOM,BEAM_TOP
******3D SOLID CONTINUUM AND CONTINUUM SHELL************
*ELEMENT,TYPE=C3D8R,ELSET=SOLID
6,11,12,13,14,15,16,17,18
*ELEMENT,TYPE=C3D8,ELSET=SOLID
66,61,62,63,64,65,66,67,68
*ELEMENT,TYPE=SC8R,ELSET=SHELL
7,19,20,21,22,23,24,25,26
*ELEMENT,TYPE=C3D8R,ELSET=SOLID_REF
106,111,112,113,114,115,116,117,118
*ELEMENT,TYPE=C3D8,ELSET=SOLID_REF
166,161,162,163,164,165,166,167,168
*ELEMENT,TYPE=SC8R,ELSET=SHELL_REF
107,119,120,121,122,123,124,125,126
**
*NSET,NSET=SOLID_BOTTOM
11,12,15,16,19,20,23,24,
61,62,65,66,
111,112,115,116,119,120,123,124,
161,162,165,166
*NSET,NSET=SOLID_TOP
13,14,17,18,21,22,25,26,
63,64,67,68,
113,114,117,118,121,122,125,126,
163,164,167,168
*NSET,NSET=SOLID_LEFT
11,14,15,18,19,22,26,23,
61,64,65,68,
111,114,115,118,119,122,126,123,
161,164,165,168
*NSET,NSET=SOLID_RIGHT
12,13,17,16,24,20,21,25,
62,63,67,66,
112,113,117,116,124,120,121,125,
162,163,167,166
*NSET,NSET=SOLID_FRONT
15,16,17,18,23,24,25,26,
65,66,67,68,
115,116,117,118,123,124,125,126,
165,166,167,168
*NSET,NSET=SOLID_BACK
11,12,13,14,19,20,21,22,
61,62,63,64,
111,112,113,114,119,120,121,122,
161,162,163,164
*NSET,NSET=SOLID_ALL
SOLID_BOTTOM,SOLID_TOP,SOLID_LEFT,SOLID_RIGHT,SOLID_FRONT,SOLID_BACK
******CONVENTIONAL SHELL***********
*ELEMENT,TYPE=S4R,ELSET=SHELL
8,27,28,29,30
*ELEMENT,TYPE=S4,ELSET=SHELL
78,77,78,79,80
*ELEMENT,TYPE=S4RS,ELSET=SHELL
9,31,32,33,34
*ELEMENT,TYPE=S4R,ELSET=SHELL_REF
108,127,128,129,130
*ELEMENT,TYPE=S4,ELSET=SHELL_REF
178,177,178,179,180
*ELEMENT,TYPE=S4RS,ELSET=SHELL_REF
109,131,132,133,134
**
*NSET,NSET=SHELL_BOTTOM
27,28,31,32,
77,78,
127,128,131,132,
177,178
*NSET,NSET=SHELL_TOP
30,29,34,33,
80,79,
130,129,134,133,
180,179
*NSET,NSET=SHELL_LEFT
27,30,31,34,
77,80,
127,130,131,134,
177,180
*NSET,NSET=SHELL_RIGHT
28,29,32,33,
78,79,
128,129,132,133,
178,179
*NSET,NSET=SHELL_ALL
SHELL_BOTTOM,SHELL_TOP,SHELL_LEFT,SHELL_RIGHT
*****PLANE STRESS,PLANE STRAIN,AXIAL SYMMETRIC,MEMBRANE****
*ELEMENT,TYPE=CPS4R,ELSET=SOLID_PS
10,35,36,37,38
*ELEMENT,TYPE=CPE4R,ELSET=SOLID_PE
11,39,40,41,42
*ELEMENT,TYPE=CAX4R,ELSET=SOLID_AX
12,43,44,45,46
*ELEMENT,TYPE=M3D4R,ELSET=MEMBRANE
13,47,48,49,50
*ELEMENT,TYPE=M3D4,ELSET=MEMBRANE
14,51,52,53,54
*ELEMENT,TYPE=CPS4R,ELSET=SOLID_REF_PS
110,135,136,137,138
*ELEMENT,TYPE=CPE4R,ELSET=SOLID_REF_PE
111,139,140,141,142
*ELEMENT,TYPE=CAX4R,ELSET=SOLID_REF_AX
112,143,144,145,146
*ELEMENT,TYPE=M3D4R,ELSET=MEMBRANE_REF
113,147,148,149,150
*ELEMENT,TYPE=M3D4,ELSET=MEMBRANE_REF
114,151,152,153,154
**
*NSET,NSET=PLANE_BOTTOM
35,36,39,40,43,44,47,48,
51,52,
135,136,139,140,143,144,147,148,
151,152
*NSET,NSET=PLANE_TOP
38,37,42,41,45,46,50,49,
54,53,
138,137,142,141,145,146,150,149,
154,153
*NSET,NSET=PLANE_LEFT
35,38,39,42,43,46,47,50,
51,54,
135,138,139,142,143,146,147,150,
151,154
*NSET,NSET=PLANE_RIGHT
36,37,40,41,44,45,48,49,
52,53,
136,137,140,141,144,145,148,149,
152,153
*NSET,NSET=PLANE_ALL
PLANE_BOTTOM,PLANE_TOP,PLANE_LEFT,PLANE_RIGHT
*ELSET,ELSET=EALL
TRUSS,BEAM,SOLID,
SOLID_PS,
SOLID_PE,
SOLID_AX,
SOLID,MEMBRANE,SHELL,
TRUSS_REF,BEAM_REF,SOLID_REF,
SOLID_REF_PS,
SOLID_REF_PE,
SOLID_REF_AX,
MEMBRANE_REF,SHELL_REF
*ELSET,ELSET=QA_TEST_ALLELEMENTS
EALL
**************************************
*BOUNDARY
TRUSS_BOTTOM,1,2
TRUSS_ALL,3,3
BEAM_BOTTOM,1,2
BEAM_ALL,3,6
SOLID_BOTTOM,2,2
SOLID_LEFT,1,1
SOLID_ALL,3,3
SHELL_BOTTOM,2,2
SHELL_LEFT,1,1
SHELL_ALL,3,6
PLANE_BOTTOM,2,2
PLANE_LEFT,1,1
PLANE_ALL,3,3
***************************************
*SOLID SECTION,MATERIAL=MAT_MISES,ELSET=TRUSS,CONTROLS=CONT
0.1
*BEAM SECTION,MATERIAL=MAT_MISES,ELSET=BEAM,SECTION=PIPE
0.1, 0.01
*SHELL SECTION,MATERIAL=MAT_MISES,ELSET=SHELL,CONTROLS=CONT
0.1
*SOLID SECTION,MATERIAL=MAT_MISES,ELSET=SOLID,CONTROLS=CONT
1.0
*SOLID SECTION,MATERIAL=MAT_MISES,ELSET=SOLID_PS,CONTROLS=CONT
1.0
*SOLID SECTION,MATERIAL=MAT_MISES,ELSET=SOLID_PE,CONTROLS=CONT
1.0
*SOLID SECTION,MATERIAL=MAT_MISES,ELSET=SOLID_AX,CONTROLS=CONT
1.0
*MEMBRANE SECTION,MATERIAL=MAT_MISES,ELSET=MEMBRANE,CONTROLS=CONT
0.1
*SOLID SECTION,MATERIAL=MAT_MISES_REF,ELSET=TRUSS_REF,CONTROLS=CONT
0.1
*BEAM SECTION,MATERIAL=MAT_MISES_REF,ELSET=BEAM_REF,SECTION=PIPE
0.1, 0.01
*SHELL SECTION,MATERIAL=MAT_MISES_REF,ELSET=SHELL_REF,CONTROLS=CONT
0.1
*SOLID SECTION,MATERIAL=MAT_MISES_REF,ELSET=SOLID_REF,CONTROLS=CONT
1.0
*SOLID SECTION,MATERIAL=MAT_MISES_REF,ELSET=SOLID_REF_PS,CONTROLS=CONT
1.0
*SOLID SECTION,MATERIAL=MAT_MISES_REF,ELSET=SOLID_REF_PE,CONTROLS=CONT
1.0
*SOLID SECTION,MATERIAL=MAT_MISES_REF,ELSET=SOLID_REF_AX,CONTROLS=CONT
1.0
*MEMBRANE SECTION,MATERIAL=MAT_MISES_REF,ELSET=MEMBRANE_REF,CONTROLS=CONT
0.1
*SECTION CONTROLS,ELEMENT DELETION=YES,NAME=CONT
***************************************
*MATERIAL,NAME=MAT_MISES
*Density
2700
*Elastic, dependencies = 1
7.00E+010,	0.33,,  0.0
7.00E+004,      0.33,,  1.0
*Plastic, dependencies = 1
1.45E+08,	0.00E+00,,   0.0
1.51E+08,	3.44E-03,,   0.0
1.60E+08,	8.59E-03,,   0.0
1.73E+08,	1.63E-02,,   0.0
1.89E+08,	2.79E-02,,   0.0
2.08E+08,	4.53E-02,,   0.0
2.31E+08,	7.14E-02,,   0.0
2.54E+08,	1.11E-01,,   0.0
2.76E+08,	1.69E-01,,   0.0
2.97E+08,	2.57E-01,,   0.0
3.24E+08,	3.90E-01,,   0.0
3.61E+08,	5.88E-01,,   0.0
4.17E+08,	8.85E-01,,   0.0
5.00E+08,	1.33E+00,,   0.0
6.25E+08,	2.00E+00,,   0.0
6.25E+08,	2.10E+00,,   0.0
1.45E+02,	0.00E+00,,   1.0
1.45E+02,	2.00E+00,,   1.0
*User defined field, properties=3
** eqpsCrit,  failDisp,  Dmax
     0.52  ,    0.1  ,  1.0
*Depvar, Delete=3
 3
 1, PEEQ  ,"Equivalent Plastic Strain"
 2, DAMAGE,"Damage Variable"
 3, STATUS,"Material Point Status"
**
*MATERIAL,NAME=MAT_MISES_REF
*Density
2700
*Elastic
7.00E+010,	0.33,
*Plastic
1.45E+08,	0.00E+00
1.51E+08,	3.44E-03
1.60E+08,	8.59E-03
1.73E+08,	1.63E-02
1.89E+08,	2.79E-02
2.08E+08,	4.53E-02
2.31E+08,	7.14E-02
2.54E+08,	1.11E-01
2.76E+08,	1.69E-01
2.97E+08,	2.57E-01
3.24E+08,	3.90E-01
3.61E+08,	5.88E-01
4.17E+08,	8.85E-01
5.00E+08,	1.33E+00
6.25E+08,	2.00E+00
6.25E+08,	2.10E+00
*DAMAGE INITIATION, CRITERION=DUCTILE
 0.52
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT
 0.1
**
*AMPLITUDE,NAME=AMP1,DEFINITION=DECAY
-1,1,0,-0.05
*STEP
*DYNAMIC,EXPLICIT
,0.05
*BOUNDARY,TYPE=DISPLACEMENT,AMP=AMP1
TRUSS_TOP,2,2,2
TRUSS_TOP,1,1,0
BEAM_TOP,2,2,2
BEAM_TOP,1,1,0
SOLID_TOP,2,2,1
SOLID_RIGHT,1,1,1
SHELL_TOP,2,2,1
SHELL_RIGHT,1,1,1
PLANE_TOP,2,2,1
PLANE_RIGHT,1,1,1
**
**
*OUTPUT,FIELD,NUMBER INTERVAL=20,TIME MARKS=YES
*ELEMENT OUTPUT,ELSET=QA_TEST_ALLELEMENTS
STATUS,SDV
PEEQ,S
*NODE OUTPUT,NSET=NALL
U,
*OUTPUT,HISTORY,TIME INTERVAL=0.0001
*ELEMENT OUTPUT,ELSET=EALL
STATUS,SDV,FV
PEEQ,MISES
SDEG, DMICRT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIME INTERVAL=0.025
*END STEP