*HEADING
 shear criterion with druker-prager plasticity
*RESTART,WRITE,NUMBER INTERVAL=10
*PREPRINT,MODEL=YES
*NODE,NSET=NALL
******3D SOLID CONTINUUM AND CONTINUUM SHELL************
*ELEMENT,TYPE=C3D8R,ELSET=SOLID
*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=S4RS,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=SOLID
*ELEMENT,TYPE=CPE4R,ELSET=SOLID
*ELEMENT,TYPE=CAX4R,ELSET=SOLID
*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=EALL
SOLID,MEMBRANE,SHELL
*ELSET,ELSET=QA_TEST_ALLELEMENTS
EALL
**************************************
*BOUNDARY
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
***************************************
*ORIENTATION,NAME=RECT
*SHELL SECTION,MATERIAL=MAT_DP,ELSET=SHELL,CONTROLS=CONT,ORIENT=RECT
0.1
*SOLID SECTION,MATERIAL=MAT_DP,ELSET=SOLID,CONTROLS=CONT,ORIENT=RECT
1.0
*MEMBRANE SECTION,MATERIAL=MAT_DP,ELSET=MEMBRANE,CONTROLS=CONT,ORIENT=RECT
0.1
*SECTION CONTROLS,ELEMENT DELETION=YES,MAX DEGRADATION=0.75,NAME=CONT
***************************************
*MATERIAL,NAME=MAT_DP
*DENSITY
8.96E3
*ELASTIC
124.E9,0.34
*DRUCKER PRAGER
*DRUCKER PRAGER HARDENING
*DAMAGE INITIATION,CRITERION=SHEAR,KS=0.5
2E-2,0
4E-2,5
*DAMAGE EVOLUTION,TYPE=DISPLACEMENT,SOFTENING=EXPONENTIAL
3E-2,3
*AMPLITUDE,NAME=AMP1,DEFINITION=DECAY
-1,1,0,-5
*STEP
*DYNAMIC,EXPLICIT
,0.5
*BOUNDARY,TYPE=DISPLACEMENT,AMP=AMP1
SOLID_TOP,2,2,1
SOLID_RIGHT,1,1,0.5
SHELL_TOP,2,2,1
SHELL_RIGHT,1,1,0.5
PLANE_TOP,2,2,1
PLANE_RIGHT,1,1,0.5
MEMBRANE_TOP,2,2,2
MEMBRANE_RIGHT,1,1,1
**
**
*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.001
*ELEMENT OUTPUT,ELSET=EALL
STATUS,SDEG,DMICRT
PEEQ,LEP,SP,MISES,PRESS
*OUTPUT,HISTORY,VARIABLE=PRESELECTED,TIME INTERVAL=0.25
*END STEP