A gradient-based non-local GTN model: Explicit finite element simulation of ductile damage and fracture

Non-local models have over the years been established as an effective approach to solve the pathological mesh dependency problem observed in finite element simulation of strain-softening materials. This paper presents the formulation, implementation, and application of a gradient -based non-local Gurson-Tvergaard-Needleman (GTN) model for explicit finite element analysis. The porosity is taken as the non-local variable where the increment in porosity is averaged over the volume using an implicit gradient model. The gradient model is implemented in Abaqus/Explicit by utilising the coupled thermal-mechanical solver, which proves to be both a simple and computationally efficient approach. Due to the use of an explicit integration scheme, a transient term is introduced to the partial differential equation of the gradient formulation. The non-local GTN model is compared to the local counterpart for increasing mesh refinements using a plane strain shear band specimen, a plane strain tension specimen, and a plane strain compact tension specimen. The proposed approach can remedy the pathological mesh dependency problem. For the plane strain tension specimen, it is shown that the non -local GTN model will preserve the fracture mode (slant versus cup-cup fracture mode) when the mesh is refined. The non-local GTN model is also able to predict the same fracture mode as observed in ductile tearing experiments. However, non-local averaging can over-smooth the fields and exclude the slant fracture mode from occurring if the material length Lc is too large.

Automated summary

This paper presents a gradient-based non-local GTN model to solve the pathological mesh dependency problem in finite element simulation of strain-softening materials. The porosity is averaged using an implicit gradient model, and the gradient model is implemented in Abaqus/Explicit. The proposed approach can effectively remedy the mesh dependency problem and accurately predict the fracture mode of the material.