*HEADING msfld criterion with mises plasticity, msfld definition. *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 *******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_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=MSFLD,DEFINITION=MSFLD 2.75978E-01, -4.54644E-01 2.40007E-01, -3.63981E-01 2.12554E-01, -2.72860E-01 1.98881E-01, -1.82110E-01 1.85146E-01, -9.08005E-02 1.71473E-01, 0.00000E+00 1.91985E-01, 1.00127E-01 2.33148E-01, 2.00221E-01 2.94889E-01, 3.00092E-01 3.70235E-01, 4.00078E-01 4.45732E-01, 5.00000E-01 5.28012E-01, 6.00122E-01 6.03406E-01, 7.00085E-01 6.65166E-01, 7.99892E-01 6.92519E-01, 9.00082E-01 7.06200E-01, 1.00000E+00 *DAMAGE EVOLUTION,TYPE=ENERGY,SOFTENING=EXPONENTIAL 1E9 *AMPLITUDE,NAME=AMP1,DEFINITION=DECAY -1,1,0,-0.05 *INITIAL CONDITIONS,TYPE=DAMAGE INITIATION,CRITERION=MSFLD SHELL,0.15 SOLID,0.15 MEMBRANE,0.15 *STEP *DYNAMIC,EXPLICIT ,0.01 *BOUNDARY,TYPE=DISPLACEMENT,AMP=AMP1 CONTINUUM_TOP,2,2,0 CONTINUUM_RIGHT,1,1,0 SHELL_TOP,2,2,0 SHELL_RIGHT,1,1,0 PLANE_TOP,2,2,0 PLANE_RIGHT,1,1,0 MEMBRANE_TOP,2,2,0 MEMBRANE_RIGHT,1,1,0 ** ** *OUTPUT,FIELD,NUMBER INTERVAL=20,TIME MARKS=YES *ELEMENT OUTPUT,ELSET=QA_TEST_ALLELEMENTS DMICRT S *NODE OUTPUT,NSET=NALL U, *OUTPUT,HISTORY,TIME INTERVAL=0.0001 *ELEMENT OUTPUT,ELSET=EALL DMICRT *OUTPUT,HISTORY,VARIABLE=PRESELECT,TIME INTERVAL=0.025 *END STEP