A finite element implementation of a large deformation gradient-damage theory for fracture with Abaqus user material subroutines
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.