*HEADING
 fld criterion with mises 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,4,0,0
10,5,0,0
11,5,1,0
12,4,1,0
13,8,0,0
14,9,0,0
15,9,1,0
16,8,1,0
17,12,0,0
18,13,0,0
19,13,1,0
20,12,1,0
21,16,0,0
22,17,0,0
23,17,1,0
24,16,1,0
25,20,0,0
26,22,0,0
27,22,2,0
28,20,2,0
29,27,0,0
30,28,0,0
31,28,1,0
32,27,1,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=S4,ELSET=SHELL
22,29,30,31,32
*ELEMENT,TYPE=S4R,ELSET=SHELL
3,13,14,15,16
*NSET,NSET=SHELL_BOTTOM
9,10,13,14
29,30
*NSET,NSET=SHELL_TOP
12,11,16,15
32,31
*NSET,NSET=SHELL_LEFT
9,12,13,16
29,32
*NSET,NSET=SHELL_RIGHT
10,11,14,15
30,31
*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_MISES,ELSET=SHELL,CONTROLS=CONT
0.1
*SOLID SECTION,MATERIAL=MAT_MISES,ELSET=SOLID,CONTROLS=CONT
1.0
*MEMBRANE SECTION,MATERIAL=MAT_MISES,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_MISES
*Density
7800
*Elastic	
2.10E+011,	0.3
*Plastic, Hardening=Isotropic, Rate=0
7.76E+08,	0,	
8.09E+08,	0.01,	
8.29E+08,	2.00E-02,	
8.42E+08,	3.00E-02,	
8.66E+08,	6.00E-02,	
8.83E+08,	1.00E-01,	
8.95E+08,	1.50E-01,	
9.10E+08,	2.50E-01,	
9.22E+08,	4.00E-01,	
9.53E+08,	2.00E+00,	
*Plastic, Hardening=Isotropic,Rate=0.001
7.91E+08,	0.00E+00,	
8.24E+08,	1.00E-02,	
8.46E+08,	2.00E-02,	
8.63E+08,	3.00E-02,	
8.99E+08,	6.00E-02,	
9.31E+08,	1.00E-01,	
9.58E+08,	1.50E-01,	
9.95E+08,	2.50E-01,	
1.03E+09,	4.00E-01,	
1.17E+09,	2.00E+00,	
*Plastic, Hardening=Isotropic, Rate=0.01
7.99E+08,	0.00E+00,	
8.31E+08,	1.00E-02,	
8.55E+08,	2.00E-02,	
8.74E+08,	3.00E-02,	
9.16E+08,	6.00E-02,	
9.55E+08,	1.00E-01,	
9.89E+08,	1.50E-01,	
1.04E+09,	2.50E-01,	
1.09E+09,	4.00E-01,	
1.28E+09,	2.00E+00,	
*Plastic, Hardening=Isotropic, Rate=0.1
8.05E+08,	0.00E+00,	
8.38E+08,	1.00E-02,	
8.63E+08,	2.00E-02,	
8.84E+08,	3.00E-02,	
9.33E+08,	6.00E-02,	
9.78E+08,	1.00E-01,	
1.02E+09,	1.50E-01,	
1.08E+09,	2.50E-01,	
1.14E+09,	4.00E-01,	
1.39E+09,	2.00E+00,	
*Plastic, Hardening=Isotropic, Rate=1.
8.08E+08,	0.00E+00,	
8.42E+08,	1.00E-02,	
8.69E+08,	2.00E-02,	
8.93E+08,	3.00E-02,	
9.46E+08,	6.00E-02,	
9.98E+08,	1.00E-01,	
1.05E+09,	1.50E-01,	
1.12E+09,	2.50E-01,	
1.19E+09,	4.00E-01,	
1.49E+09,	2.00E+00,	
*Plastic, Hardening=Isotropic, Rate=10.
8.10E+08,	0.00E+00,	
8.46E+08,	1.00E-02,	
8.76E+08,	2.00E-02,	
9.01E+08,	3.00E-02,	
9.60E+08,	6.00E-02,	
1.02E+09,	1.00E-01,	
1.07E+09,	1.50E-01,	
1.15E+09,	2.50E-01,	
1.24E+09,	4.00E-01,	
1.60E+09,	2.00E+00,	
*Plastic, Hardening=Isotropic, Rate=100.
8.12E+08,	0.00E+00,	
8.50E+08,	1.00E-02,	
8.82E+08,	2.00E-02,	
9.09E+08,	3.00E-02,	
9.74E+08,	6.00E-02,	
1.04E+09,	1.00E-01,	
1.10E+09,	1.50E-01,	
1.19E+09,	2.50E-01,	
1.28E+09,	4.00E-01,	
1.70E+09,	2.00E+00,	
*Plastic, Hardening=Isotropic, Rate=1000.
8.15E+08,	0.00E+00,	
8.55E+08,	1.00E-02,	
8.88E+08,	2.00E-02,	
9.17E+08,	3.00E-02,	
9.87E+08,	6.00E-02,	
1.06E+09,	1.00E-01,	
1.13E+09,	1.50E-01,	
1.23E+09,	2.50E-01,	
1.33E+09,	4.00E-01,	
1.81E+09,	2.00E+00,	
*DAMAGE INITIATION,CRITERION=FLD
0.2756,		-0.1253
0.2371,		-0.0863
0.2056,		-0.0561
0.1867,		-0.034
0.1674,		-0.0152
0.1485,		0.
0.1578,		0.0158
0.1813,		0.0363
0.2166,		0.065
0.2567,		0.1027
0.2918,		0.1459
0.3266,		0.196
0.3531,		0.2472
0.3688,		0.295
0.3643,		0.3279
0.3531,		0.3531
*DAMAGE EVOLUTION,TYPE=ENERGY,SOFTENING=EXPONENTIAL
1E9
*AMPLITUDE,NAME=AMP1,DEFINITION=DECAY
-1,1,0,-0.05
*amplitude,name=amp2
 0,0,0.05,0.2
*STEP,nlgeom=yes,inc=1000
*static
 0.001,0.05,,0.002
*BOUNDARY,TYPE=DISPLACEMENT,AMP=AMP2
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,freq=10
*ELEMENT OUTPUT,ELSET=QA_TEST_ALLELEMENTS
STATUS,SDEG,DMICRT
PEEQ,S
*NODE OUTPUT,NSET=NALL
U,
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=EALL
STATUS,SDEG,DMICRT
PEEQ,LEP,SP,MISES,PRESS
*END STEP