*HEADING
 fld criterion with drucker-prager plasticity
*RESTART,WRITE,NUMBER INTERVAL=10
*PREPRINT,MODEL=YES
*NODE,NSET=NALL
1,0,0,0
2,1,0,0
3,1,1,0
4,0,1,0
5,0,0,0.1
6,1,0,0.1
7,1,1,0.1
8,0,1,0.1
9,2,0,0
10,3,0,0
11,3,1,0
12,2,1,0
13,4,0,0
14,5,0,0
15,5,1,0
16,4,1,0
17,6,0,0
18,7,0,0
19,7,1,0
20,6,1,0
21,8,0,0
22,9,0,0
23,9,1,0
24,8,1,0
25,10,0,0
26,12,0,0
27,12,2,0
28,10,2,0
*******CONTINUUM SHELL**********
*ELEMENT,TYPE=SC8R,ELSET=SHELL
1,1,2,3,4,5,6,7,8
*NSET,NSET=CONTINUUM_BOTTOM
1,2,5,6
*NSET,NSET=CONTINUUM_TOP
3,4,8,7
*NSET,NSET=CONTINUUM_LEFT
1,4,5,8
*NSET,NSET=CONTINUUM_RIGHT
2,3,6,7
*NSET,NSET=CONTINUUM_BACK
1,2,3,4
*NSET,NSET=CONTINUUM_FRONT
5,6,7,8
*NSET,NSET=CONTINUUM_ALL
CONTINUUM_BOTTOM,CONTINUUM_TOP,CONTINUUM_LEFT,CONTINUUM_RIGHT
CONTINUUM_BACK,CONTINUUM_FRONT
*******SHELL*****************
*ELEMENT,TYPE=S4R,ELSET=SHELL
2,9,10,11,12
*ELEMENT,TYPE=S4RS,ELSET=SHELL
3,13,14,15,16
*NSET,NSET=SHELL_BOTTOM
9,10,13,14
*NSET,NSET=SHELL_TOP
12,11,16,15
*NSET,NSET=SHELL_LEFT
9,12,13,16
*NSET,NSET=SHELL_RIGHT
10,11,14,15
*NSET,NSET=SHELL_ALL
SHELL_BOTTOM,SHELL_TOP,SHELL_LEFT,SHELL_RIGHT
************CPS4R AND MEMBRANE***************
*ELEMENT,TYPE=CPS4R,ELSET=SOLID
4,17,18,19,20
*ELEMENT,TYPE=M3D4R,ELSET=MEMBRANE
5,21,22,23,24
*ELEMENT,TYPE=M3D4,ELSET=MEMBRANE
6,25,26,27,28
*NSET,NSET=PLANE_BOTTOM
17,18,21,22,25,26
*NSET,NSET=PLANE_LEFT
17,20,21,24,25,28
*NSET,NSET=PLANE_TOP
19,20,23,24
*NSET,NSET=MEMBRANE_TOP
27,28
*NSET,NSET=PLANE_RIGHT
18,19,22,23
*NSET,NSET=MEMBRANE_RIGHT
26,27
*NSET,NSET=PLANE_ALL
PLANE_BOTTOM,PLANE_TOP,PLANE_LEFT,PLANE_RIGHT,
MEMBRANE_RIGHT,MEMBRANE_TOP
***************BOUNDARY************************
*BOUNDARY
CONTINUUM_BOTTOM,2,2
CONTINUUM_LEFT,1,1
CONTINUUM_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
***************************************
*SHELL SECTION,MATERIAL=MAT_DP,ELSET=SHELL,CONTROLS=CONT
0.1
*SOLID SECTION,MATERIAL=MAT_DP,ELSET=SOLID,CONTROLS=CONT
1.0
*MEMBRANE SECTION,MATERIAL=MAT_DP,ELSET=MEMBRANE,CONTROLS=CONT
0.1
*SECTION CONTROLS,ELEMENT DELETION=YES,MAX DEGRADATION=0.75,NAME=CONT
*ELSET,ELSET=EALL
SOLID,SHELL,MEMBRANE
*ELSET,ELSET=QA_TEST_ALLELEMENTS
EALL
***************************************
*MATERIAL,NAME=MAT_DP
*DENSITY
8.96E3
*ELASTIC
124.E9,0.34
*DRUCKER PRAGER
40.,0.8,20.
*DRUCKER PRAGER HARDENING
90e6,0.0
*DAMAGE INITIATION,CRITERION=FLD
3E-2
*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
CONTINUUM_TOP,2,2,1
CONTINUUM_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.001
*ELEMENT OUTPUT,ELSET=EALL
STATUS,SDEG,DMICRT
PEEQ,LEP,SP,MISES,PRESS
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIME INTERVAL=0.25
*END STEP