*HEADING
 ductile failure with Mises plasticity
*RESTART,WRITE,NUMBER INTERVAL=10
*PREPRINT,MODEL=YES
*NODE,NSET=NALL
********TRUSS ELEMENTS*******************
*ELEMENT,TYPE=T2D2,ELSET=TRUSS
*ELEMENT,TYPE=T3D2,ELSET=TRUSS
*NSET,NSET=TRUSS_BOTTOM
*NSET,NSET=TRUSS_TOP
*NSET,NSET=TRUSS_ALL
TRUSS_BOTTOM,TRUSS_TOP
*******BEAM AND AXIAL SHELL**************
*ELEMENT,TYPE=B31,ELSET=BEAM
*ELEMENT,TYPE=B31,ELSET=BEAM
*ELEMENT,TYPE=SAX1,ELSET=SHELL
*NSET,NSET=BEAM_BOTTOM
*NSET,NSET=BEAM_TOP
*NSET,NSET=BEAM_ALL
BEAM_BOTTOM,BEAM_TOP
******3D SOLID CONTINUUM AND CONTINUUM SHELL************
*ELEMENT,TYPE=C3D8R,ELSET=SOLID3D
*ELEMENT,TYPE=C3D8,ELSET=SOLID3D
*ELEMENT,TYPE=SC8R,ELSET=SHELL
*NSET,NSET=SOLID_BOTTOM
*NSET,NSET=SOLID_TOP
*NSET,NSET=SOLID_LEFT
*NSET,NSET=SOLID_RIGHT
*NSET,NSET=SOLID_FRONT
*NSET,NSET=SOLID_BACK
*NSET,NSET=SOLID_ALL
SOLID_BOTTOM,SOLID_TOP,SOLID_LEFT,SOLID_RIGHT,SOLID_FRONT,SOLID_BACK
******CONVENTIONAL SHELL***********
*ELEMENT,TYPE=S4R,ELSET=SHELL
*ELEMENT,TYPE=S4,ELSET=SHELL
*ELEMENT,TYPE=S4R,ELSET=SHELL
*NSET,NSET=SHELL_BOTTOM
*NSET,NSET=SHELL_TOP
*NSET,NSET=SHELL_LEFT
*NSET,NSET=SHELL_RIGHT
*NSET,NSET=SHELL_ALL
SHELL_BOTTOM,SHELL_TOP,SHELL_LEFT,SHELL_RIGHT
*****PLANE STRESS,PLANE STRAIN,AXIAL SYMMETRIC,MEMBRANE****
*ELEMENT,TYPE=CPS4R,ELSET=SOLID2D
*ELEMENT,TYPE=CPE4R,ELSET=SOLID2D
*ELEMENT,TYPE=CAX4R,ELSET=SOLID3D
*ELEMENT,TYPE=M3D4R,ELSET=MEMBRANE
*ELEMENT,TYPE=M3D4,ELSET=MEMBRANE
*NSET,NSET=PLANE_BOTTOM
*NSET,NSET=PLANE_TOP
*NSET,NSET=MEMBRANE_TOP
*NSET,NSET=PLANE_LEFT
*NSET,NSET=PLANE_RIGHT
*NSET,NSET=MEMBRANE_RIGHT
*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
*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
*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
*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
*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
*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
*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
*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
*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
*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
**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,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