*HEADING
 ductile criterion with drucker-prager plasticity
*RESTART,WRITE,NUMBER INTERVAL=10
*PREPRINT,MODEL=YES
*NODE,NSET=NALL
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,7,0,0
20,8,0,0
21,8,1,0
22,7,1,0
23,7,0,0.1
24,8,0,0.1
25,8,1,0.1
26,7,1,0.1
27,9,0,0
28,10,0,0
29,10,1,0
30,9,1,0
31,11,0,0
32,12,0,0
33,12,1,0
34,11,1,0
35,13,0,0
36,14,0,0
37,14,1,0
38,13,1,0
39,15,0,0
40,16,0,0
41,16,1,0
42,15,1,0
43,17,0,0
44,18,0,0
45,18,1,0
46,17,1,0
47,19,0,0
48,20,0,0
49,20,1,0
50,19,1,0
51,21,0,0
52,23,0,0
53,23,2,0
54,21,2,0
******3D SOLID CONTINUUM AND CONTINUUM SHELL************
*ELEMENT,TYPE=C3D8R,ELSET=SOLID3D
6,11,12,13,14,15,16,17,18
*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
*NSET,NSET=SOLID_TOP
13,14,17,18,21,22,25,26
*NSET,NSET=SOLID_LEFT
11,14,15,18,19,22,26,23
*NSET,NSET=SOLID_RIGHT
12,13,17,16,24,20,21,25
*NSET,NSET=SOLID_FRONT
15,16,17,18,23,24,25,26
*NSET,NSET=SOLID_BACK
11,12,13,14,19,20,21,22
*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=S4R,ELSET=SHELL
9,31,32,33,34
*NSET,NSET=SHELL_BOTTOM
27,28,31,32
*NSET,NSET=SHELL_TOP
30,29,34,33
*NSET,NSET=SHELL_LEFT
27,30,31,34
*NSET,NSET=SHELL_RIGHT
28,29,32,33
*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
SOLID,MEMBRANE,SHELL
*ELSET,ELSET=QA_TEST_ALLELEMENTS
EALL
**************************************
*BOUNDARY
SOLID_BOTTOM,2,2
SOLID_LEFT,1,1
SOLID_BACK,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
 1.0, 0.0, 0.0, 0.0, 1.0, 0.0
 3, 0.0
*SHELL SECTION,MATERIAL=MAT_DP,ELSET=SHELL,CONTROLS=CONT,ORIENT=RECT
0.1
*SOLID SECTION,MATERIAL=MAT_DP,ELSET=SOLID2D,CONTROLS=CONT,ORIENT=RECT
1.0
*SOLID SECTION,MATERIAL=MAT_DP,ELSET=SOLID3D,CONTROLS=CONT,ORIENT=RECT
*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
40.,0.8,20.
*DRUCKER PRAGER HARDENING
90E6,0.0
*DAMAGE INITIATION,CRITERION=MaxPS
50.0e+7
*DAMAGE INITIATION,CRITERION=DUCTILE
3E-2
*DAMAGE EVOLUTION,TYPE=DISPLACEMENT,SOFTENING=EXPONENTIAL
 3E-2,3
*AMPLITUDE,NAME=AMP1,DEFINITION=DECAY
-1,1,0,-5
*STEP,nlgeom=yes,inc=1000
*static
 0.01,0.5,,0.01
*BOUNDARY,TYPE=DISPLACEMENT,AMP=AMP1
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
*NODE OUTPUT,NSET=NALL
 U,
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=EALL
 DMICRT,
 PEEQ,
 MISES,
 PRESS,
 TRIAX
 SDEG
*END STEP