*HEADING
 ductile failure with Mises plasticity; user defined element characteristic length
*RESTART,WRITE,NUMBER INTERVAL=10
*PREPRINT,MODEL=YES
*NODE,NSET=NALL
**
********TRUSS ELEMENTS*******************
*ELEMENT,TYPE=T2D2,ELSET=TRUSS
*ELEMENT,TYPE=T3D2,ELSET=TRUSS
*ELEMENT,TYPE=T2D2,ELSET=TRUSS_REF
*ELEMENT,TYPE=T3D2,ELSET=TRUSS_REF
**
*NSET,NSET=TRUSS_BOTTOM
*NSET,NSET=TRUSS_TOP
*NSET,NSET=TRUSS_ALL
TRUSS_BOTTOM,TRUSS_TOP
*******AXIAL SHELL**************
*ELEMENT,TYPE=SAX1,ELSET=SHELL
*ELEMENT,TYPE=SAX1,ELSET=SHELL_REF
**
*NSET,NSET=SAX1_BOTTOM
*NSET,NSET=SAX1_TOP
*NSET,NSET=SAX1_ALL
SAX1_BOTTOM,SAX1_TOP
******3D SOLID CONTINUUM AND CONTINUUM SHELL************
*ELEMENT,TYPE=C3D8R,ELSET=SOLID
*ELEMENT,TYPE=C3D8,ELSET=SOLID
*ELEMENT,TYPE=SC8R,ELSET=SHELL
*ELEMENT,TYPE=C3D8R,ELSET=SOLID_REF
*ELEMENT,TYPE=C3D8,ELSET=SOLID_REF
*ELEMENT,TYPE=SC8R,ELSET=SHELL_REF
**
*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=S4RS,ELSET=SHELL
*ELEMENT,TYPE=S4R,ELSET=SHELL_REF
*ELEMENT,TYPE=S4,ELSET=SHELL_REF
*ELEMENT,TYPE=S4RS,ELSET=SHELL_REF
**
*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=SOLID_PS
*ELEMENT,TYPE=CPE4R,ELSET=SOLID_PE
*ELEMENT,TYPE=CAX4R,ELSET=SOLID_AX
*ELEMENT,TYPE=M3D4R,ELSET=MEMBRANE
*ELEMENT,TYPE=M3D4,ELSET=MEMBRANE
*ELEMENT,TYPE=CPS4R,ELSET=SOLID_REF_PS
*ELEMENT,TYPE=CPE4R,ELSET=SOLID_REF_PE
*ELEMENT,TYPE=CAX4R,ELSET=SOLID_REF_AX
*ELEMENT,TYPE=M3D4R,ELSET=MEMBRANE_REF
*ELEMENT,TYPE=M3D4,ELSET=MEMBRANE_REF
**
*NSET,NSET=PLANE_BOTTOM
*NSET,NSET=PLANE_TOP
*NSET,NSET=PLANE_LEFT
*NSET,NSET=PLANE_RIGHT
*NSET,NSET=PLANE_ALL
PLANE_BOTTOM,PLANE_TOP,PLANE_LEFT,PLANE_RIGHT
*ELSET,ELSET=EALL
TRUSS,SOLID,
SOLID_PS,
SOLID_PE,
SOLID_AX,
SOLID,MEMBRANE,SHELL,
TRUSS_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
SAX1_BOTTOM,1,2
SAX1_ALL,6,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
*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
*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+005,      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
*User defined field, properties=3
** eqpsCrit,  failDisp,  Dmax
*Depvar, Delete=3
 1, PEEQ  ,"Equivalent Plastic Strain"
 2, DAMAGE,"Damage Variable"
 3, STATUS,"Material Point Status"
 4, CHARLENGTH,"Element Characteristic Length"
*characteristic length,definition=user,properties=1
**
*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
*DAMAGE INITIATION, CRITERION=DUCTILE
 0.52
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT
 0.1
*characteristic length,definition=user,properties=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
SAX1_TOP,2,2,2
SAX1_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