Extension of the GISSMO fracture model for thin-walled structures under combined tensile and bending loads

Ductile fracture prediction for thin-walled structures requires computationally efficient simulation tools able to approximately represent the key effects that occur on the sub-thickness scale. Unfortunately, classical shell elements, typically used to model deformation and failure of thin components, are inherently in a state of plane stress and therefore unable to capture the through-thickness stress distribution. This is consequential considering recent studies that demonstrated the differences between fracture in bending vs. in-plane tension. A laterally constrained thin metal plate under in-plane tension is likely to fracture under significantly lower plastic strain than the same plate under bending (with fracture initiating on the tensile side), even though both these conditions are examples of plane strain tension. Under in-plane tension fracture is preceded by a through-thickness neck, and therefore higher stress triaxiality than in the case of plane strain bending, where necking is absent. To account for these differences, we propose an extension of the GISSMO model, which relies on a simple fracture criterion based on stress-dependent fracture strain defined by the user (i.e. fracture locus). The fracture strain is typically determined experimentally using a combination of in-plane tensile tests under varying degree of lateral constraint and shear tests. The proposed extension involves defining a separate fracture locus for bending, also determined experimentally using bending tests. Fracture occurs when the equivalent plastic strain reaches its critical level represented by interpolation between the two bounding cases, i.e. bending and in-plane tension, with a bending index Ω used as an interpolation parameter.