A unified Abaqus implementation of the phase field fracture method using only a user material subroutine

Abstract

We present a simple and robust implementation of the phase field fracture method in Abaqus. Unlike previous works, only a user material (UMAT) subroutine is used. This is achieved by exploiting the analogy between the phase field balance equation and heat transfer, which avoids the need for a user element mesh and enables taking advantage of Abaqus' in-built features. A unified theoretical framework and its implementation are presented, suitable for any arbitrary choice of crack density function and fracture driving force. Specifically, the framework is exemplified with the so-called AT1, AT2 and phase field-cohesive zone models (PF-CZM). Both staggered and monolithic solution schemes are handled. We demonstrate the potential and robustness of this new implementation by addressing several paradigmatic 2D and 3D boundary value problems. The numerical examples show how the current implementation can be used to reproduce numerical and experimental results from the literature, and efficiently capture advanced features such as complex crack trajectories, crack nucleation from arbitrary sites and contact problems. The code developed can be downloaded from www.empaneda.com/codes.

References in corpus (9)

• A phase field formulation for hydrogen assisted cracking
• Phase field modeling of quasi-static and dynamic crack propagation: COMSOL implementation and case studies
• A novel framework to model the fatigue behavior of brittle materials based on a variational phase-field approach
• Phase field modelling of crack propagation in functionally graded materials
• A phase field formulation for dissolution-driven stress corrosion cracking
• A phase field model for elastic-gradient-plastic solids undergoing hydrogen embrittlement
• An assessment of phase field fracture: crack initiation and growth
• Phase field modelling of fracture and fatigue in Shape Memory Alloys
• Phase-field Fracture Modelling of Thin Monolithic and Laminated Glass Plates under Quasi-static Bending

Cited by in corpus (4)

• A modified combined active-set Newton method for solving phase-field fracture into the monolithic limit
• A finite element implementation of a large deformation gradient-damage theory for fracture with Abaqus user material subroutines
• A generalised phase field model for fatigue crack growth in elastic-plastic solids with an efficient monolithic solver
• A smoothed particle hydrodynamics approach for phase field modeling of brittle fracture