*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=S4R,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
*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 = HC 
1.0,0.05,0.01,0.1
*DAMAGE EVOLUTION,TYPE=ENERGY,SOFTENING=EXPONENTIAL
 1E8
*AMPLITUDE,NAME=AMP1,DEFINITION=DECAY
-1,1,0,-0.05
*STEP,nlgeom=yes,inc=1000
*static
 0.001,0.05,,0.001
*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,freq=10
*ELEMENT OUTPUT,ELSET=QA_TEST_ALLELEMENTS
 DMICRT,
 PEEQ,
 S,
 TRIAX
*NODE OUTPUT,NSET=NALL
 U,
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=EALL
 DMICRT,
 PEEQ,
 MISES,
 PRESS,
 TRIAX,
 SDEG
*END STEP