A finite element implementation of a large deformation gradient-damage theory for fracture with Abaqus user material subroutines

TL;DR

This paper presents a detailed finite element implementation of a large deformation gradient-damage fracture model in Abaqus, addressing challenges in modeling damage in materials undergoing large strains, with validation across various material behaviors.

Contribution

It introduces a robust Abaqus UMAT/UMATHT implementation of a gradient-damage fracture theory for large deformations, including detailed code and validation benchmarks.

Findings

Successful implementation across different material models

Robustness demonstrated through benchmark problems

Applicable to large deformation and rate-dependent plasticity

Abstract

Recent advancements in computations have enabled the application of various modeling approaches to predict fracture and failure, such as the gradient-damage (phasefield) method. Several existing studies have leveraged the heat equation solver in Abaqus to model gradient-damage, due to its mathematical resemblance to the heat equation. Particular care is required when extending the approach to large deformation scenarios due to differences in the referential and spatial configurations, especially since the heat equation in Abaqus is solved in the spatial configuration, whereas most gradient-damage frameworks are formulated in the referential configuration. This work provides a pedagogic view of an appropriate Abaqus implementation of a gradient-damage theory for fracture in materials undergoing large deformation using Abaqus UMAT and UMATHT user subroutines. Key benchmark problems from the literature are used to demonstrate the robustness of our implementation across various materials exhibiting different constitutive behaviors, such as non-linear elasticity, linear elasticity, and large deformation rate-dependent plasticity, ensuring its applicability regardless of the specific material constitutive choice. The details of the implementation, along with the codes, which are a direct outcome of this work, are also provided.