A deep energy method formulation for phase-field modelling of ductile fracture
Abstract
This work extends the Deep Energy Method to phase-field modelling of ductile fracture in elasto-plastic solids. Its key contribution is an incremental DEM-based formulation that couples history-dependent J2 plasticity and plasticity-dependent fracture toughness with separate energy-optimised displacement and phase-field representations. The formulation combines the second-order Ambrosio-Tortorelli phase-field fracture model with small-strain J2 plasticity and linear isotropic hardening. A volumetric-deviatoric split defines the elastic crack-driving energy, while accumulated equivalent plastic strain governs the degradation of local fracture toughness. Plastic internal variables are updated locally through radial return. The reduced energy is minimised alternately with respect to the displacement and phase fields, represented by a ReZero residual neural network and spatially localised Gaussian radial basis functions with trainable centres, respectively. This split addresses the diffuse phase-field profiles observed in preliminary calculations with a shared global network. Alternating optimisation is combined with transfer learning to reuse trainable parameters and optimiser states across loading increments, while the converged constitutive history is carried forward separately at quadrature points. The numerical results demonstrate close agreement with finite-element reference solutions in predicting force-displacement responses, energy evolution, plastic deformation and crack propagation. The framework captures the main features of ductile fracture, although accuracy becomes more sensitive to loading increments near abrupt fracture events.
Related
- Formation of grain boundaries in ductile single crystals under plane-strain simple shear: a block-coordinate finite element method
- Phase-resolved coastal wave modelling using a non-hydrostatic multilayer formulation
- A new formulation of Discontinuous Galerkin method for interface capturing approaches
- Gluonic nucleon energy correlators and fracture functions for Color Glass Condensate
- Modeling of variable angle tow composites using a variable separation method based on a unified formulation
BibTeX
@misc{plungė2026deep,
title = {{A deep energy method formulation for phase-field modelling of ductile fracture}},
author = {Nojus Plungė and Peter Brommer and Rachel S. Edwards and Emmanouil G. Kakouris},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.19094},
eprint = {2610.19094},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.19094},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - A deep energy method formulation for phase-field modelling of ductile fracture
AU - Nojus Plungė
AU - Peter Brommer
AU - Rachel S. Edwards
AU - Emmanouil G. Kakouris
PY - 2026
DA - 2026/10/13
PB - arXiv
DO - 10.48550/arxiv.2610.19094
UR - https://arxiv.org/abs/2610.19094
AB - This work extends the Deep Energy Method to phase-field modelling of ductile fracture in elasto-plastic solids. Its key contribution is an incremental DEM-based formulation that couples history-dependent J2 plasticity and plasticity-dependent fracture toughness with separate energy-optimised displacement and phase-field representations. The formulation combines the second-order Ambrosio-Tortorelli phase-field fracture model with small-strain J2 plasticity and linear isotropic hardening. A volumetric-deviatoric split defines the elastic crack-driving energy, while accumulated equivalent plastic strain governs the degradation of local fracture toughness. Plastic internal variables are updated locally through radial return. The reduced energy is minimised alternately with respect to the displacement and phase fields, represented by a ReZero residual neural network and spatially localised Gaussian radial basis functions with trainable centres, respectively. This split addresses the diffuse phase-field profiles observed in preliminary calculations with a shared global network. Alternating optimisation is combined with transfer learning to reuse trainable parameters and optimiser states across loading increments, while the converged constitutive history is carried forward separately at quadrature points. The numerical results demonstrate close agreement with finite-element reference solutions in predicting force-displacement responses, energy evolution, plastic deformation and crack propagation. The framework captures the main features of ductile fracture, although accuracy becomes more sensitive to loading increments near abrupt fracture events.
ER -