*HEADING shear failure criterion with Mises plasticity *RESTART,WRITE,NUMBER INTERVAL=10 *PREPRINT,MODEL=YES *NODE,NSET=NALL *NODE,NSET=NALL 1,0,0,0 2,0,1,0 3,1,0,0 4,1,1,0 5,2,0,0 6,2,1,0 7,3,0,0 8,3,1,0 9,4,0,0 10,4,1,0 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,10,0,0 20,11,0,0 21,11,1,0 22,10,1,0 23,10,0,0.1 24,11,0,0.1 25,11,1,0.1 26,10,1,0.1 27,15,0,0 28,16,0,0 29,16,1,0 30,15,1,0 31,20,0,0 32,21,0,0 33,21,1,0 34,20,1,0 35,25,0,0 36,26,0,0 37,26,1,0 38,25,1,0 39,30,0,0 40,31,0,0 41,31,1,0 42,30,1,0 43,35,0,0 44,36,0,0 45,36,1,0 46,35,1,0 47,40,0,0 48,41,0,0 49,41,1,0 50,40,1,0 51,45,0,0 52,47,0,0 53,47,2,0 54,45,2,0 61,55,0,0 62,56,0,0 63,56,1,0 64,55,1,0 65,55,0,1 66,56,0,1 67,56,1,1 68,55,1,1 77,60,0,0 78,61,0,0 79,61,1,0 80,60,1,0 ********TRUSS ELEMENTS******************* *ELEMENT,TYPE=T2D2,ELSET=TRUSS 1,1,2 *ELEMENT,TYPE=T3D2,ELSET=TRUSS 2,3,4 *NSET,NSET=TRUSS_BOTTOM 1,3 *NSET,NSET=TRUSS_TOP 2,4 *NSET,NSET=TRUSS_ALL TRUSS_BOTTOM,TRUSS_TOP *******BEAM AND AXIAL SHELL************** *ELEMENT,TYPE=B31,ELSET=BEAM 3,5,6 *ELEMENT,TYPE=B31,ELSET=BEAM 4,7,8 *ELEMENT,TYPE=SAX1,ELSET=SHELL 5,9,10 *NSET,NSET=BEAM_BOTTOM 5,7,9 *NSET,NSET=BEAM_TOP 6,8,10 *NSET,NSET=BEAM_ALL BEAM_BOTTOM,BEAM_TOP ******3D SOLID CONTINUUM AND CONTINUUM SHELL************ *ELEMENT,TYPE=C3D8R,ELSET=SOLID3D 6,11,12,13,14,15,16,17,18 *ELEMENT,TYPE=C3D8,ELSET=SOLID3D 66,61,62,63,64,65,66,67,68 *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 61,62,65,66 *NSET,NSET=SOLID_TOP 13,14,17,18,21,22,25,26 63,64,67,68 *NSET,NSET=SOLID_LEFT 11,14,15,18,19,22,26,23 61,64,65,68, *NSET,NSET=SOLID_RIGHT 12,13,17,16,24,20,21,25 62,63,67,66, *NSET,NSET=SOLID_FRONT 15,16,17,18,23,24,25,26 65,66,67,68, *NSET,NSET=SOLID_BACK 11,12,13,14,19,20,21,22 61,62,63,64, *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=S4,ELSET=SHELL 78,77,78,79,80 *ELEMENT,TYPE=S4R,ELSET=SHELL 9,31,32,33,34 *NSET,NSET=SHELL_BOTTOM 27,28,31,32 77,78 *NSET,NSET=SHELL_TOP 30,29,34,33 80,79 *NSET,NSET=SHELL_LEFT 27,30,31,34 77,80 *NSET,NSET=SHELL_RIGHT 28,29,32,33 78,79 *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 TRUSS,BEAM,SOLID,MEMBRANE,SHELL *ELSET,ELSET=QA_TEST_ALLELEMENTS EALL ************************************** *BOUNDARY TRUSS_BOTTOM,1,2 TRUSS_ALL,3,3 BEAM_BOTTOM,1,2 BEAM_ALL,3,6 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 *************************************** *SOLID SECTION,MATERIAL=MAT_MISES,ELSET=TRUSS,CONTROLS=CONT 0.1 *BEAM SECTION,MATERIAL=MAT_MISES,ELSET=BEAM,SECTION=CIRC 0.1 *SHELL SECTION,MATERIAL=MAT_MISES,ELSET=SHELL,CONTROLS=CONT 0.1 *SOLID SECTION,MATERIAL=MAT_MISES,ELSET=SOLID2D,CONTROLS=CONT 1.0 *SOLID SECTION,MATERIAL=MAT_MISES,ELSET=SOLID3D,CONTROLS=CONT *MEMBRANE SECTION,MATERIAL=MAT_MISES,ELSET=MEMBRANE,CONTROLS=CONT 0.1 *SECTION CONTROLS,ELEMENT DELETION=YES,MAX DEGRADATION=0.75,NAME=CONT *************************************** *MATERIAL,NAME=MAT_MISES *Density 2700 *Elastic 7.00E+010, 0.33 *Plastic, Hardening=Isotropic, Rate=0 1.45E+08, 0.00E+00 1.51E+08, 3.44E-03 1.60E+08, 8.59E-03 1.73E+08, 1.63E-02 1.89E+08, 2.79E-02 2.08E+08, 4.53E-02 2.31E+08, 7.14E-02 2.54E+08, 1.11E-01 2.76E+08, 1.69E-01 2.97E+08, 2.57E-01 3.24E+08, 3.90E-01 3.61E+08, 5.88E-01 4.17E+08, 8.85E-01 5.00E+08, 1.33E+00 6.25E+08, 2.00E+00 6.25E+08, 2.10E+00 *Plastic, Hardening=Isotropic, Rate=1.00e-02 1.46E+08, 0.00E+00 1.53E+08, 3.44E-03 1.62E+08, 8.59E-03 1.74E+08, 1.63E-02 1.90E+08, 2.79E-02 2.10E+08, 4.53E-02 2.32E+08, 7.14E-02 2.55E+08, 1.11E-01 2.77E+08, 1.69E-01 2.99E+08, 2.57E-01 3.25E+08, 3.90E-01 3.62E+08, 5.88E-01 4.18E+08, 8.85E-01 5.01E+08, 1.33E+00 6.26E+08, 2.00E+00 6.26E+08, 2.10E+00 *Plastic, Hardening=Isotropic, Rate=1.00e-01 1.49E+08, 0.00E+00 1.56E+08, 3.44E-03 1.65E+08, 8.59E-03 1.77E+08, 1.63E-02 1.93E+08, 2.79E-02 2.13E+08, 4.53E-02 2.35E+08, 7.14E-02 2.58E+08, 1.11E-01 2.80E+08, 1.69E-01 3.02E+08, 2.57E-01 3.28E+08, 3.90E-01 3.65E+08, 5.88E-01 4.21E+08, 8.85E-01 5.04E+08, 1.33E+00 6.29E+08, 2.00E+00 6.29E+08, 2.10E+00 *Plastic, Hardening=Isotropic, Rate=1.00e+00 1.54E+08, 0.00E+00 1.60E+08, 3.44E-03 1.69E+08, 8.59E-03 1.82E+08, 1.63E-02 1.98E+08, 2.79E-02 2.18E+08, 4.53E-02 2.40E+08, 7.14E-02 2.63E+08, 1.11E-01 2.85E+08, 1.69E-01 3.06E+08, 2.57E-01 3.32E+08, 3.90E-01 3.69E+08, 5.88E-01 4.24E+08, 8.85E-01 5.05E+08, 1.33E+00 6.26E+08, 2.00E+00 6.26E+08, 2.10E+00 *Plastic, Hardening=Isotropic, Rate=1.00e+01 1.61E+08, 0.00E+00 1.68E+08, 3.44E-03 1.77E+08, 8.59E-03 1.89E+08, 1.63E-02 2.05E+08, 2.79E-02 2.25E+08, 4.53E-02 2.47E+08, 7.14E-02 2.69E+08, 1.11E-01 2.90E+08, 1.69E-01 3.10E+08, 2.57E-01 3.33E+08, 3.90E-01 3.65E+08, 5.88E-01 4.09E+08, 8.85E-01 4.70E+08, 1.33E+00 5.46E+08, 2.00E+00 5.46E+08, 2.10E+00 *Plastic, Hardening=Isotropic, Rate=2.50e+01 1.65E+08, 0.00E+00 1.72E+08, 3.44E-03 1.81E+08, 8.59E-03 1.93E+08, 1.63E-02 2.08E+08, 2.79E-02 2.28E+08, 4.53E-02 2.49E+08, 7.14E-02 2.70E+08, 1.11E-01 2.89E+08, 1.69E-01 3.05E+08, 2.57E-01 3.22E+08, 3.90E-01 3.43E+08, 5.88E-01 3.69E+08, 8.85E-01 3.97E+08, 1.33E+00 4.18E+08, 2.00E+00 4.18E+08, 2.10E+00 *Plastic, Hardening=Isotropic, Rate=1.00e+02 1.73E+08, 0.00E+00 1.79E+08, 3.44E-03 1.88E+08, 8.59E-03 2.00E+08, 1.63E-02 2.15E+08, 2.79E-02 2.34E+08, 4.53E-02 2.55E+08, 7.14E-02 2.74E+08, 1.11E-01 2.90E+08, 1.69E-01 3.03E+08, 2.57E-01 3.13E+08, 3.90E-01 3.24E+08, 5.88E-01 3.34E+08, 8.85E-01 3.39E+08, 1.33E+00 3.31E+08, 2.00E+00 3.31E+08, 2.10E+00 *Plastic, Hardening=Isotropic, Rate=1.00e+03 1.85E+08, 0.00E+00 1.91E+08, 3.44E-03 2.00E+08, 8.59E-03 2.12E+08, 1.63E-02 2.27E+08, 2.79E-02 2.46E+08, 4.53E-02 2.66E+08, 7.14E-02 2.85E+08, 1.11E-01 3.00E+08, 1.69E-01 3.11E+08, 2.57E-01 3.20E+08, 3.90E-01 3.27E+08, 5.88E-01 3.33E+08, 8.85E-01 3.33E+08, 1.33E+00 3.19E+08, 2.00E+00 3.19E+08, 2.10E+00 *Plastic, Hardening=Isotropic, Rate=1.00e+04 1.85E+08, 0.00E+00 1.91E+08, 3.44E-03 2.00E+08, 8.59E-03 2.12E+08, 1.63E-02 2.27E+08, 2.79E-02 2.46E+08, 4.53E-02 2.66E+08, 7.14E-02 2.85E+08, 1.11E-01 3.00E+08, 1.69E-01 3.11E+08, 2.57E-01 3.20E+08, 3.90E-01 3.27E+08, 5.88E-01 3.33E+08, 8.85E-01 3.33E+08, 1.33E+00 3.19E+08, 2.00E+00 3.19E+08, 2.10E+00 *Damage Initiation, Criterion = Shear, KS=0.3 0.5929, -10., 0.001 0.5929, 1.4236, 0.001 0.6123, 1.4625, 0.001 0.6364, 1.5013, 0.001 0.6653, 1.5401, 0.001 0.6993, 1.5789, 0.001 0.7386, 1.6177, 0.001 0.7836, 1.6566, 0.001 0.8344, 1.6954, 0.001 0.8917, 1.7342, 0.001 0.9557, 1.773, 0.001 1.027, 1.8118, 0.001 1.106, 1.8506, 0.001 1.1935, 1.8895, 0.001 1.2901, 1.9283, 0.001 1.3965, 1.9671, 0.001 1.5135, 2.0059, 0.001 1.6421, 2.0447, 0.001 1.7831, 2.0835, 0.001 1.9377, 2.1224, 0.001 2.107, 2.1612, 0.001 2.2923, 2.2, 0.001 2.2923, 10., 0.001 0.5361, -10., 250. 0.5361, 1.4236, 250. 0.5831, 1.4625, 250. 0.6313, 1.5013, 250. 0.6808, 1.5401, 250. 0.7317, 1.5789, 250. 0.7842, 1.6177, 250. 0.8382, 1.6566, 250. 0.8941, 1.6954, 250. 0.9517, 1.7342, 250. 1.0114, 1.773, 250. 1.0732, 1.8118, 250. 1.1372, 1.8506, 250. 1.2035, 1.8895, 250. 1.2724, 1.9283, 250. 1.3439, 1.9671, 250. 1.4182, 2.0059, 250. 1.4955, 2.0447, 250. 1.5759, 2.0835, 250. 1.6595, 2.1224, 250. 1.7467, 2.1612, 250. 1.8374, 2.2, 250. 1.8374, 10., 250. *DAMAGE EVOLUTION,TYPE=ENERGY,SOFTENING=EXPONENTIAL 1E8 *AMPLITUDE,NAME=AMP1,DEFINITION=DECAY -1,1,0,-0.05 *STEP,nlgeom=yes,inc=1000 *static 0.001,0.05,,0.001 *BOUNDARY,TYPE=DISPLACEMENT,AMP=AMP1 TRUSS_TOP,2,2,2 TRUSS_TOP,1,1,0 BEAM_TOP,2,2,2 BEAM_TOP,1,1,0 SOLID_TOP,2,2,2 SOLID_RIGHT,1,1,2.0 SHELL_TOP,2,2,2 SHELL_RIGHT,1,1,2.0 PLANE_TOP,2,2,2 PLANE_RIGHT,1,1,2.0 MEMBRANE_TOP,2,2,4 MEMBRANE_RIGHT,1,1,4.0 ** ** *OUTPUT,FIELD,freq=10 *ELEMENT OUTPUT,ELSET=QA_TEST_ALLELEMENTS STATUS,SDEG,DMICRT,SHRRATIO PEEQ,S *NODE OUTPUT,NSET=NALL U, *OUTPUT,HISTORY *ELEMENT OUTPUT,ELSET=EALL STATUS,SDEG,DMICRT,SHRRATIO PEEQ,LEP,SP,MISES,PRESS *END STEP