*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,39,0,0
53,39,2,0
54,37,2,0
61,42,0,0
62,43,0,0
63,43,1,0
64,42,1,0
65,42,0,1
66,43,0,1
67,43,1,1
68,42,1,1
77,45,0,0
78,46,0,0
79,46,1,0
80,45,1,0
********TRUSS ELEMENTS*******************
*ELEMENT,TYPE=T2D2,ELSET=TRUSS
1,1,2
*ELEMENT,TYPE=T3D2,ELSET=TRUSS
2,3,4
*NSET,NSET=TRUSS_BOTTOM
1,3
*NSET,NSET=TRUSS_TOP
2,4
*NSET,NSET=TRUSS_ALL
TRUSS_BOTTOM,TRUSS_TOP
*******BEAM AND AXIAL SHELL**************
*ELEMENT,TYPE=B31,ELSET=BEAM
3,5,6
*ELEMENT,TYPE=B31,ELSET=BEAM
4,7,8
*ELEMENT,TYPE=SAX1,ELSET=SHELL
5,9,10
*NSET,NSET=BEAM_BOTTOM
5,7,9
*NSET,NSET=BEAM_TOP
6,8,10
*NSET,NSET=BEAM_ALL
BEAM_BOTTOM,BEAM_TOP
******3D SOLID CONTINUUM AND CONTINUUM SHELL************
*ELEMENT,TYPE=C3D8R,ELSET=SOLID3D
6,11,12,13,14,15,16,17,18
*ELEMENT,TYPE=C3D8,ELSET=SOLID3D
66,61,62,63,64,65,66,67,68
*ELEMENT,TYPE=SC8R,ELSET=SHELL
7,19,20,21,22,23,24,25,26
*NSET,NSET=SOLID_BOTTOM
11,12,15,16,19,20,23,24
61,62,65,66
*NSET,NSET=SOLID_TOP
13,14,17,18,21,22,25,26
63,64,67,68
*NSET,NSET=SOLID_LEFT
11,14,15,18,19,22,26,23
61,64,65,68,
*NSET,NSET=SOLID_RIGHT
12,13,17,16,24,20,21,25
62,63,67,66,
*NSET,NSET=SOLID_FRONT
15,16,17,18,23,24,25,26
65,66,67,68,
*NSET,NSET=SOLID_BACK
11,12,13,14,19,20,21,22
61,62,63,64,
*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
*NSET,NSET=SHELL_BOTTOM
27,28,31,32
77,78
*NSET,NSET=SHELL_TOP
30,29,34,33
80,79
*NSET,NSET=SHELL_LEFT
27,30,31,34
77,80
*NSET,NSET=SHELL_RIGHT
28,29,32,33
78,79
*NSET,NSET=SHELL_ALL
SHELL_BOTTOM,SHELL_TOP,SHELL_LEFT,SHELL_RIGHT
*****PLANE STRESS,PLANE STRAIN,AXIAL SYMMETRIC,MEMBRANE****
*ELEMENT,TYPE=CPS4R,ELSET=SOLID2D
10,35,36,37,38
*ELEMENT,TYPE=CPE4R,ELSET=SOLID2D
11,39,40,41,42
*ELEMENT,TYPE=CAX4R,ELSET=SOLID3D
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
*NSET,NSET=PLANE_BOTTOM
35,36,39,40,43,44,47,48,
51,52
*NSET,NSET=PLANE_TOP
38,37,42,41,45,46,50,49
*NSET,NSET=MEMBRANE_TOP
54,53
*NSET,NSET=PLANE_LEFT
35,38,39,42,43,46,47,50
51,54
*NSET,NSET=PLANE_RIGHT
36,37,40,41,44,45,48,49
*NSET,NSET=MEMBRANE_RIGHT
52,53
*NSET,NSET=PLANE_ALL
PLANE_BOTTOM,PLANE_TOP,PLANE_LEFT,PLANE_RIGHT,
MEMBRANE_RIGHT,MEMBRANE_TOP
*ELSET,ELSET=SOLID
SOLID2D,SOLID3D
*ELSET,ELSET=EALL
TRUSS,BEAM,SOLID,MEMBRANE,SHELL
*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=CIRC
0.1
*SHELL SECTION,MATERIAL=MAT_MISES,ELSET=SHELL,CONTROLS=CONT
0.1
*SOLID SECTION,MATERIAL=MAT_MISES,ELSET=SOLID2D,CONTROLS=CONT
1.0
*SOLID SECTION,MATERIAL=MAT_MISES,ELSET=SOLID3D,CONTROLS=CONT
1.0
*MEMBRANE SECTION,MATERIAL=MAT_MISES,ELSET=MEMBRANE,CONTROLS=CONT
0.1
*SECTION CONTROLS,ELEMENT DELETION=YES,MAX DEGRADATION=0.75,NAME=CONT
***************************************
*MATERIAL,NAME=MAT_MISES
*Density
2700
*Elastic
7.00E+010,	0.33
*Plastic, Hardening=Isotropic, Rate=0
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
*Plastic, Hardening=Isotropic, Rate=1.00e-02
1.46E+08,	0.00E+00
1.53E+08,	3.44E-03
1.62E+08,	8.59E-03
1.74E+08,	1.63E-02
1.90E+08,	2.79E-02
2.10E+08,	4.53E-02
2.32E+08,	7.14E-02
2.55E+08,	1.11E-01
2.77E+08,	1.69E-01
2.99E+08,	2.57E-01
3.25E+08,	3.90E-01
3.62E+08,	5.88E-01
4.18E+08,	8.85E-01
5.01E+08,	1.33E+00
6.26E+08,	2.00E+00
6.26E+08,	2.10E+00
*Plastic, Hardening=Isotropic, Rate=1.00e-01
1.49E+08,	0.00E+00
1.56E+08,	3.44E-03
1.65E+08,	8.59E-03
1.77E+08,	1.63E-02
1.93E+08,	2.79E-02
2.13E+08,	4.53E-02
2.35E+08,	7.14E-02
2.58E+08,	1.11E-01
2.80E+08,	1.69E-01
3.02E+08,	2.57E-01
3.28E+08,	3.90E-01
3.65E+08,	5.88E-01
4.21E+08,	8.85E-01
5.04E+08,	1.33E+00
6.29E+08,	2.00E+00
6.29E+08,	2.10E+00
*Plastic, Hardening=Isotropic, Rate=1.00e+00
1.54E+08,	0.00E+00
1.60E+08,	3.44E-03
1.69E+08,	8.59E-03
1.82E+08,	1.63E-02
1.98E+08,	2.79E-02
2.18E+08,	4.53E-02
2.40E+08,	7.14E-02
2.63E+08,	1.11E-01
2.85E+08,	1.69E-01
3.06E+08,	2.57E-01
3.32E+08,	3.90E-01
3.69E+08,	5.88E-01
4.24E+08,	8.85E-01
5.05E+08,	1.33E+00
6.26E+08,	2.00E+00
6.26E+08,	2.10E+00
*Plastic, Hardening=Isotropic, Rate=1.00e+01
1.61E+08,	0.00E+00
1.68E+08,	3.44E-03
1.77E+08,	8.59E-03
1.89E+08,	1.63E-02
2.05E+08,	2.79E-02
2.25E+08,	4.53E-02
2.47E+08,	7.14E-02
2.69E+08,	1.11E-01
2.90E+08,	1.69E-01
3.10E+08,	2.57E-01
3.33E+08,	3.90E-01
3.65E+08,	5.88E-01
4.09E+08,	8.85E-01
4.70E+08,	1.33E+00
5.46E+08,	2.00E+00
5.46E+08,	2.10E+00
*Plastic, Hardening=Isotropic, Rate=2.50e+01
1.65E+08,	0.00E+00
1.72E+08,	3.44E-03
1.81E+08,	8.59E-03
1.93E+08,	1.63E-02
2.08E+08,	2.79E-02
2.28E+08,	4.53E-02
2.49E+08,	7.14E-02
2.70E+08,	1.11E-01
2.89E+08,	1.69E-01
3.05E+08,	2.57E-01
3.22E+08,	3.90E-01
3.43E+08,	5.88E-01
3.69E+08,	8.85E-01
3.97E+08,	1.33E+00
4.18E+08,	2.00E+00
4.18E+08,	2.10E+00
*Plastic, Hardening=Isotropic, Rate=1.00e+02
1.73E+08,	0.00E+00
1.79E+08,	3.44E-03
1.88E+08,	8.59E-03
2.00E+08,	1.63E-02
2.15E+08,	2.79E-02
2.34E+08,	4.53E-02
2.55E+08,	7.14E-02
2.74E+08,	1.11E-01
2.90E+08,	1.69E-01
3.03E+08,	2.57E-01
3.13E+08,	3.90E-01
3.24E+08,	5.88E-01
3.34E+08,	8.85E-01
3.39E+08,	1.33E+00
3.31E+08,	2.00E+00
3.31E+08,	2.10E+00
*Plastic, Hardening=Isotropic, Rate=1.00e+03
1.85E+08,	0.00E+00
1.91E+08,	3.44E-03
2.00E+08,	8.59E-03
2.12E+08,	1.63E-02
2.27E+08,	2.79E-02
2.46E+08,	4.53E-02
2.66E+08,	7.14E-02
2.85E+08,	1.11E-01
3.00E+08,	1.69E-01
3.11E+08,	2.57E-01
3.20E+08,	3.90E-01
3.27E+08,	5.88E-01
3.33E+08,	8.85E-01
3.33E+08,	1.33E+00
3.19E+08,	2.00E+00
3.19E+08,	2.10E+00
*Plastic, Hardening=Isotropic, Rate=1.00e+04
1.85E+08,	0.00E+00
1.91E+08,	3.44E-03
2.00E+08,	8.59E-03
2.12E+08,	1.63E-02
2.27E+08,	2.79E-02
2.46E+08,	4.53E-02
2.66E+08,	7.14E-02
2.85E+08,	1.11E-01
3.00E+08,	1.69E-01
3.11E+08,	2.57E-01
3.20E+08,	3.90E-01
3.27E+08,	5.88E-01
3.33E+08,	8.85E-01
3.33E+08,	1.33E+00
3.19E+08,	2.00E+00
3.19E+08,	2.10E+00
*Damage Initiation, Criterion = Ductile 
2.08E+02,	-3.33,		0.001
207.9706,	-0.333,		0.001
118.7949,	-0.267,		0.001
67.8571,	-0.2,		0.001
38.7612,	-0.133,		0.001
22.1419,	-0.0667,	0.001
12.6496,	0,		0.001
7.229,		0.0667,		0.001
4.1354,		0.133,		0.001
2.3728,		0.2,		0.001
1.3739,		0.267,		0.001
0.8173,		0.333,		0.001
0.5237,		0.4,		0.001
0.3988,		0.467,		0.001
0.4021,		0.533,		0.001
0.535,		0.6,		0.001
0.84,		0.667,		0.001
0.5008,		0.730,		0.001
0.187,		0.851,		0.001
0.0457,		1.02,		0.001
0.0076,		1.24,		0.001
0.0001,		1.51,		0.001
0.0001,		3.33,		0.001
465.088,	-3.33,		250.
465.088,	-0.333,		250.
250.1916,	-0.267,		250.
134.5893,	-0.2,		250.
72.4019,	-0.133,		250.
38.9487,	-0.0667,	250.
20.9533,	0.,		250.
11.2736,	0.0667,		250.
6.068,		0.133,		250.
3.2707,		0.2,		250.
1.7715,		0.267,		250.
0.9753,		0.333,		250.
0.5662,		0.4,		250.
0.3817,		0.467,		250.
0.3488,		0.533,		250.
0.4542,		0.6,		250.
0.74,		0.667,		250.
0.4193,		0.730,		250.
0.1421,		0.851,		250.
0.0302,		1.02,		250.
0.0042,		1.24,		250.
0.0001,		1.51,		250.
0.0001,		3.33,		250.
*DAMAGE EVOLUTION,TYPE=ENERGY,SOFTENING=EXPONENTIAL
1E8
*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
MEMBRANE_TOP,2,2,2
MEMBRANE_RIGHT,1,1,2
**
**
*OUTPUT,FIELD,NUMBER INTERVAL=20,TIME MARKS=YES
*ELEMENT OUTPUT,ELSET=QA_TEST_ALLELEMENTS
STATUS,SDEG,DMICRT
PEEQ,S
*NODE OUTPUT,NSET=NALL
U,
*OUTPUT,HISTORY,TIME INTERVAL=0.0001
*ELEMENT OUTPUT,ELSET=EALL
STATUS,SDEG,DMICRT
PEEQ,LEP,SP,MISES,PRESS,TRIAX
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIME INTERVAL=0.025
*END STEP