*HEADING
 HC ductile criterion with hill plasticity
 Input file name: damage_ductile_hill_hc_std.inp
*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_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
 1.0, 0.0, 0.0, 0.0, 1.0, 0.0
 3, 0.0
*SHELL SECTION,MATERIAL=MAT_HILL,ELSET=SHELL,CONTROLS=CONT,ORIENT=RECT
0.1
*SOLID SECTION,MATERIAL=MAT_HILL,ELSET=SOLID2D,CONTROLS=CONT,ORIENT=RECT
1.0
*SOLID SECTION,MATERIAL=MAT_HILL,ELSET=SOLID3D,CONTROLS=CONT,ORIENT=RECT
*MEMBRANE SECTION,MATERIAL=MAT_HILL,ELSET=MEMBRANE,CONTROLS=CONT,ORIENT=RECT
0.1
*SECTION CONTROLS,ELEMENT DELETION=YES,MAX DEGRADATION=0.75,NAME=CONT
***************************************
*MATERIAL,NAME=MAT_HILL
*DENSITY
8.96E3
*ELASTIC
124.E9,0.34
*PLASTIC,HARDENING=ISOTROPIC
  90.0000E6,   0.0000
 142.7256E6,   0.0040
 155.3643E6,   0.0080
 164.1189E6,   0.0120
 171.0327E6,   0.0160
 176.8365E6,   0.0200
 187.4375E6,   0.0290
 195.9536E6,   0.0380
 203.1705E6,   0.0470
 209.4871E6,   0.0560
 215.1371E6,   0.0650
 220.2701E6,   0.0740
 224.9886E6,   0.0830
 229.3661E6,   0.0920
 233.4572E6,   0.1010
 237.3040E6,   0.1100
 240.9393E6,   0.1190
 244.3896E6,   0.1280
 247.6762E6,   0.1370
 250.8171E6,   0.1460
 253.8271E6,   0.1550
 256.7187E6,   0.1640
 259.5029E6,   0.1730
 262.1888E6,   0.1820
 264.7846E6,   0.1910
 267.2973E6,   0.2000
*POTENTIAL
1.5,1.2,1.0,1.0,1.0,1.0
*DAMAGE INITIATION,CRITERION=HC
1.0,0.1,0.01,0.1
*DAMAGE EVOLUTION,TYPE=ENERGY,SOFTENING=LINEAR
 1E6
*AMPLITUDE,NAME=AMP1,DEFINITION=DECAY
-1,1,0,-5
*STEP,nlgeom=yes,inc=10000
*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