Energy-norm positive/negative projection scheme for phase-field fracture with anisotropic elasticity
Authors:
- Shank KulkarniSKShank Kulkarni
- Mohammad SafiMSMohammad Safi
- Timothy J. TrusterTTTimothy J. Truster
- Khalid AlshibliKAKhalid Alshibli
Preprint
- 6 September 2026
- preprint
- Published by Elsevier in SSRN Electronic Journal
Abstract
Phase-field theory has proven to be an efficient approach for predicting the damage and fracture of engineering materials. It can simulate complex crack patterns without previously defined crack paths. In this work, a novel framework was developed for quantifying the fracture behavior of anisotropic brittle materials. This novel approach includes a phase-field formulation for fracture with (i) distinct contributions arising from anisotropy in the elastic moduli and (ii) a positive/negative projection (PNP) scheme in the complementary-energy norm for anisotropic elasticity, so that the crack-driving energy is consistently associated with tensile deformation while suppressing unphysical damage growth under compression. The present work considers anisotropy in the elastic constitutive response, while the fracture surface energy is assumed isotropic. The resulting coupled displacement–damage problem is implemented in Abaqus/Standard through a user element subroutine (UEL), while a user material subroutine (UMAT) is used to transfer stresses, strains, and damage variables for post-processing. A monolithic solution strategy is adopted to solve the displacement and phase-field equations simultaneously. The implementation is verified using one-element tests under tensile, compressive, and shear loading, patch tests with distorted meshes, and several quasi-static fracture benchmarks, including single-edge-notched tension, compression, shear, and three-point bending problems. Parametric studies are also performed to examine the effects of mesh size, regularization length scale, and critical energy release rate on the predicted fracture response. For anisotropic trigonal materials, the material-point simulations capture the orientation-dependent elastic and damage responses, while the notched fracture simulations demonstrate orientation-dependent crack trajectories arising from elastic anisotropy.Keywords
This publication has 43 references indexed in Scilit:
- 3D implementation of push-out test in ABAQUS using the phase-field methodJournal of Mechanical Science and Technology, 2023
- A phase-field fracture model for brittle anisotropic materialsComputational Mechanics, 2022
- Brittle anisotropic fracture propagation in quartz sandstone: insights from phase-field simulationsComputational Geosciences, 2020
- Multi-phase-field modeling of anisotropic crack propagation for polycrystalline materialsComputational Mechanics, 2017
- Linear and nonlinear solvers for variational phase-field models of brittle fractureInternational Journal for Numerical Methods in Engineering, 2016
- Thermodynamically consistent phase‐field models of fracture: Variational principles and multi‐field FE implementationsInternational Journal for Numerical Methods in Engineering, 2010
- The Variational Approach to FractureJournal of Elasticity, 2008
- A finite element method for crack growth without remeshingInternational Journal for Numerical Methods in Engineering, 1999
- Reproducing kernel particle methodsInternational Journal for Numerical Methods in Fluids, 1995
- Element‐free Galerkin methodsInternational Journal for Numerical Methods in Engineering, 1994