---
title: Localizing Gradient Damage Method (LGDM)
url: https://www.emergentmind.com/topics/localizing-gradient-damage-method-lgdm
type: topic
---

# Localizing Gradient Damage Method (LGDM)

Searching arXiv for recent papers on the Localizing Gradient Damage Method and closely related gradient-damage formulations.
Localizing Gradient Damage Method (LGDM) denotes, in the literature cited here, a class of gradient-regularized damage formulations designed to suppress pathological mesh dependence while permitting damage to localize into a finite-width band. Across these formulations, the central idea is the introduction of a nonlocal field—directly the damage variable, or an auxiliary micromorphic quantity such as nonlocal chain stretch, nonlocal equivalent strain, or regularized damage invariants—together with an internal length scale that controls the width of the damage process zone. A distinctive strand of recent work, particularly for elastomers, frames LGDM as a thermodynamically consistent alternative to phase-field fracture in which fracture localization and diffuse network scission are both represented and fracture toughness may emerge as an output rather than an input [2502.02822]. Other works use the same label for phase-field-like gradient damage, modified non-local damage laws with fixed-width bands, mixed finite-element regularizations, and micromorphic tensorial damage extensions [2506.24099] [2301.06503] [2311.15918].

## 1. Terminology and conceptual scope

The nomenclature surrounding LGDM is not uniform. In the elastomer formulation of Mousavi et al., the method is a stretch-based gradient-enhanced damage model with a local damage variable \(d(X)\in[0,1]\), a perturbed-Lagrangian pressure \(p(X)\), and a nonlocal chain-stretch measure \(\bar\lambda(X)\). Damage is not assigned its own evolution PDE; instead, it is obtained from a constitutive damage function driven by the historic nonlocal stretch, while \(\bar\lambda\) satisfies a micromorphic balance [2502.02822]. In the modified non-local damage model of Saji et al., LGDM is realized through a nonlocal equivalent strain \(\bar\varepsilon\) governed by a Helmholtz-type equation, while the local damage \(d\) is updated through a chosen damage law \(d=d(\bar\varepsilon)\) [2506.24099].

In other papers summarized under the same label, LGDM is closer to gradient-damage or phase-field practice. Alkhoury et al. formulate a regularized total energy in terms of displacement \(\mathbf u\) and scalar damage \(d\), with degradation \(g(d)=(1-d)^2\), fracture toughness \(G_c\), and internal length \(l_0\), leading to a canonical Helmholtz-type damage PDE [2507.02702]. Riesselmann and Balzani employ a mixed finite-element formulation for finite strains with displacement \(u\), damage \(\alpha\), and a Lagrange multiplier \(\lambda\) enforcing irreversibility exactly, while Holthusen et al. regularize anisotropic damage by introducing nonlocal invariants \(\bar d_i\) associated with a second-order damage tensor \(D\) [2211.07964] [2311.15918].

This variation in usage suggests that LGDM is best understood not as a single constitutive model, but as a methodological family centered on localization-capable gradient regularization.

## 2. Variational structure and governing fields

The formulations collected under LGDM are variationally structured, but they differ in which field is regularized. In the elastomer model of Mousavi et al., the reference configuration is \(\Omega_0\subset\mathbb R^3\), with deformation map \(x=X+u(X)\), deformation gradient \(F=\partial x/\partial X\), Jacobian \(J=\det F\), and right Cauchy–Green tensor \(C=F^T F\). The local elastic energy is compressible neo-Hookean,
\[
\psi(F)=\frac{\mu}{2}\bigl[I_1(C)-3-2\ln J\bigr],
\]
and near-incompressibility is enforced by the perturbed-Lagrangian term \(-p(J-1)-p^2/(2\kappa)\) with \(\kappa\gg \mu\). The nonlocal contribution introduces a stiffness \(h_{nl}>0\), an internal length \(\ell>0\), and a damage-weighted relaxation function \(g(d)\) [2502.02822].

In that framework, damage is introduced through a statistical-mechanics-inspired conjugate force \(f_d(d,\bar\lambda)\). A critical chain stretch \(\lambda_{cr}>1\), steepness \(\gamma\ge 0\), maximum damage parameter \(c\in(0,1]\), and a history operator \(\mathcal H(\bar\lambda)\) define a constitutive law such that no damage occurs for \(\bar\lambda<\lambda_{cr}\), whereas for \(\bar\lambda\ge \lambda_{cr}\), \(d\) increases with the maximum historic nonlocal stretch. Solving \(f_d=0\) yields the damage function \(d=d(\bar\lambda)\), so the only genuine additional PDE is the one for \(\bar\lambda\) [2502.02822].

Saji et al. formulate the free-energy density as
\[
\psi\bigl(\varepsilon(u),\bar\varepsilon,d,\nabla\bar\varepsilon\bigr)
=
\frac12\,g(d)\,\varepsilon:\mathbf C:\varepsilon
+\frac12\,h\Bigl[(f_r(d)\varepsilon-\bar\varepsilon)^2+c\,|\nabla\bar\varepsilon|^2\Bigr],
\]
with
\[
g(d)=(1-d)+f_a(d,d^{k-1}),
\qquad
c=\ell_c^2/2.
\]
Here \(\bar\varepsilon\) carries the regularization, the stress is degraded by \(g(d)\), and the nonlocal field is driven by \(f_r(d)\varepsilon\) rather than by \(\varepsilon\) alone [2506.24099].

By contrast, the phase-field-like version summarized by Alkhoury et al. uses the regularized total energy
\[
\Pi(\mathbf u,d)
=
\int_{\Omega_0} g(d)\,\psi_0(F)\,dV_0
+\int_{\Omega_0} G_c\Bigl(\frac{1}{2l_0}d^2+\frac{l_0}{2}|\nabla d|^2\Bigr)\,dV_0
-\langle \text{external work}\rangle,
\]
which yields mechanical equilibrium together with the damage balance
\[
-2(1-d)\psi_0(F)+\frac{G_c}{l_0}d-G_c l_0 \Delta d=0.
\]
Irreversibility is imposed by replacing \(\psi_0(F)\) with the history field \(H(t)=\max_{s\le t}\psi_0(F(s))\) [2507.02702].

Micromorphic anisotropic extensions retain the same general logic. Holthusen et al. introduce a symmetric damage tensor \(D(x)\), a tuple of local invariants \(d_i(D)\), and their nonlocal counterparts \(\bar d_i(x)\), with micromorphic energy
\[
\psi^{\mathrm{micr}}
=
\frac12\sum_{i=1}^{n_d} H_i\bigl(d_i(D)-\bar d_i\bigr)^2
+
\frac12\sum_{i=1}^{n_d} A_i |\nabla \bar d_i|^2.
\]
The regularized variables may be the full six independent tensor components, three trace invariants, or a reduced volumetric-deviatoric pair [2311.15918].

## 3. Localization mechanism and suppression of damage-band broadening

The principal theoretical motivation for LGDM is the avoidance of spurious localization pathologies. Classical local damage models are mesh dependent after loss of ellipticity. Conventional gradient-damage formulations remove that mesh dependence by introducing an internal length, but several papers emphasize that they may still exhibit unrealistic widening of the damage band as damage progresses [2506.24099]. Mousavi et al. make the same criticism of classical gradient-enhanced damage models for elastomers, noting that they often fail to capture the cascade from damage to fracture and instead produce damage-zone broadening [2502.02822].

In the stretch-based elastomer formulation, two specific modeling features are introduced to overcome this deficiency. First, the local chain-stretch driver \(\lambda_{ch}=\sqrt{I_1/3}\) is bounded by replacing it with
\[
\lambda_{ch}^{\max}-\langle \lambda_{ch}^{\max}-\lambda_{ch}\rangle,
\]
so that once \(\lambda_{ch}\ge \lambda_{ch}^{\max}\), the local driver saturates. Second, nonlocal coupling is relaxed inside highly damaged regions through
\[
g(d)=(1-d)^m,\qquad m>0.
\]
According to the authors, these two features permit a physically meaningful diffuse damage field together with simultaneous fracture localization [2502.02822].

The modified non-local damage model achieves a related objective through two different thermodynamic conditions. The auxiliary degradation term
\[
f_a(d,d^{k-1})=\frac12 d^2-\frac12(d^{k-1})^2
\]
is added so that \(\partial g/\partial d\to 0\) as \(d\to 1\), and hence the thermodynamic damage driving force \(Y=-\partial\psi/\partial d\) vanishes at full damage. In parallel, the forcing term in the nonlocal equation is attenuated with
\[
f_r(d)=1-d^n,\qquad n\gg 1,
\]
so that \(f_r(d)\varepsilon\to 0\) as \(d\to 1\). The stated consequence is a non-growing constant-width damage band [2506.24099].

Earlier analytical localization studies provide a useful backdrop. In a one-dimensional rate-independent regularized damage model, the internal length \(\ell\) controls the width of the localized zone, and linearization around the onset of damage gives a critical wavelength
\[
\lambda_c=\frac{2\pi \ell}{\sqrt2},
\qquad
\frac{L_{d,2}}{2}=\frac{\pi \ell}{\sqrt2}.
\]
Modified regularization with a variable length \(\ell(\alpha)=\ell_0/(1-\alpha)^p\) yields a constant-zone-width response for \(p\approx 0.5\), whereas \(p>0.5\) produces expansion and \(p<0.5\) contraction [1412.5539]. This suggests that a central issue in LGDM is not only regularization per se, but the precise asymptotic behavior of the nonlocal coupling as the damage variable approaches its terminal state.

## 4. Irreversibility, fracture measures, and the distinction from phase-field fracture

Irreversibility is enforced in several ways. In the stretch-based elastomer model, the constitutive structure is governed by \( \dot d\ge 0 \) together with the Clausius–Planck condition \( \partial \Psi/\partial d \le 0 \); because \(d\) appears only algebraically, irreversibility is handled pointwise through the history-dependent constitutive relation \(d=d[\bar\lambda]\) [2502.02822]. In phase-field-like variants, a history field \(H(t)\) stores the maximum previous energy density and enters the damage PDE directly [2507.02702]. In the mixed finite-element formulation of Riesselmann and Balzani, irreversibility is written as a KKT system,
\[
\alpha-\bar\alpha\ge 0,\qquad \lambda\le 0,\qquad \lambda(\alpha-\bar\alpha)=0,
\]
with \(\bar\alpha=\max\{\alpha_n,\alpha\}\) and exact enforcement by a Lagrange multiplier [2211.07964].

A major conceptual distinction concerns fracture toughness. In standard phase-field formulations, \(G_c\) is a prescribed material input appearing explicitly in the crack-density part of the free energy [2507.02702]. Mousavi et al. emphasize that their LGDM does not prescribe an input \(G_c\). Instead, one loads a cracked or notched specimen in displacement control, computes the domain \(J\)-integral, and identifies the plateau in
\[
G=J(D)
\]
as the predicted fracture toughness. In their \(1\times 1\) plane-strain panel with a pre-damage zone, the plateau gives \(G_c\approx 0.104\) [2502.02822].

That predicted value was then used as the input toughness in a phase-field benchmark with the same geometry, \(\mu\), \(\kappa\), and \(\ell\). The reported force-displacement curves, \(J\)-integral versus displacement curves, and damage-field versus phase-field contours all matched to within a few percent, while LGDM retained two stated differences: \(\ell\) had a micromechanical interpretation as a chain-statistical nonlocal length, and the damage field \(d(X)\) followed chain-stretch physics rather than serving as a crack interpolant [2502.02822].

The comparison is not uniformly favorable to all gradient-enhanced damage variants. In an earlier elastomer study, the measured energy release rate during crack propagation in a phase-field model with artificial viscosity did not comply with the imposed critical energy release rate and showed non-monotonic behavior. The stretch-based GED alternative again produced \(J\) as an output, but was reported to suffer from numerical issues to which strain-based GED formulations are susceptible, including damage-zone broadening and high sensitivity of \(J\) to artificial viscosity \(\eta\) [2408.05162]. This is one of the clearest controversies in the area: some damage-based GED formulations were said to be identical to phase-field methods and to inherit their limitations, whereas other LGDM variants were explicitly designed to avoid that equivalence [2502.02822].

## 5. Discretization strategies and computational implementations

LGDM has been implemented in several computational settings, and the implementation strategy depends strongly on which field is regularized.

| Implementation line | Core discretization | Notable claim |
|---|---|---|
| Stretch-based elastomer LGDM | Taylor–Hood type \(P_2/P_1/P_1\) for \((u,p,\bar\lambda)\) | Fracture toughness extracted from \(J\)-integral plateau |
| Abaqus UMAT + UMATHT | Damage PDE mapped to heat equation | Applicable in large deformation |
| Mixed FE with Lagrange multiplier | \(P_2\)-\(u\), bubble-enriched \(P_1\)-\(\alpha\), \(P_0\)-\(\lambda\) | No penalty or stabilization |
| FFT-based nonlocal ductile failure | FFT-Galerkin mechanics + Krylov Helmholtz solve | Millions of DOFs in 3D RVEs |
| JAX-FEM differentiable GED | HEX8 with automatic differentiation | Mesh-independent, data-driven extension |
| Vectorized MATLAB LGDM | Sparse vectorized assembly | Significant speedup on consumer hardware |

In the elastomer formulation of Mousavi et al., the finite-element spaces are \(u\in[H^1]^3\), \(p\in L^2\), and \(\bar\lambda\in H^1\), discretized with vector-valued quadratics for \(u\), linears for \(p\), and linears for \(\bar\lambda\). A staggered quasi-static algorithm is used: first solve the mechanics block for \((u^k,p^k)\) via Newton–SNES, then solve the nonlocal PDE for \(\bar\lambda^k\) with bound constraints using vinewtonssls, then update \(d^k=d[\bar\lambda^k]\) algebraically until \(\|\bar\lambda^k-\bar\lambda^{k-1}\|_\infty<\text{tol}\). Near incompressibility is enforced through \(\kappa/\mu\approx 10^3\Rightarrow \nu\approx 0.499\) [2502.02822].

Alkhoury et al. translate the large-deformation damage PDE into the spatial frame to exploit Abaqus’s heat-equation solver. With the identification \(\theta\equiv d\), the transient heat equation is matched term by term via
\[
\rho C = J^{-1},\qquad
K=l_0^2 J^{-1}B,\qquad
r=2(1-d)J^{-1}H-J^{-1}d.
\]
UMATHT prescribes \(DUDT=C_{\mathrm{eff}}=1/J\) and \(FLUX=-(l_0^2 J^{-1}B)\,DTEMDX\), while UMAT uses the RPL mechanism with
\[
RPL=2(1-d)J^{-1}H-J^{-1}d,\qquad
DRPLDT=-2J^{-1}H-J^{-1}.
\]
The implementation was validated against an independent UEL code, with coincidence to within numerical precision in the reported decoupled tests [2507.02702].

Riesselmann and Balzani propose a mixed finite element for gradient damage in finite strains based on a Lagrange multiplier and a bubble-enriched damage space. The damage field uses \(P_1\)-Lagrange interpolation on vertices plus one quartic bubble at the element center, and the multiplier is element-wise constant. Because the multiplier and bubble degree of freedom are local to each element, both are eliminated by static condensation, leaving a global system with only the \(P_2\)-displacement and vertex-based \(P_1\)-damage unknowns. The formulation is reported to require no penalty or stabilization [2211.07964].

For large heterogeneous RVEs, an FFT framework solves the fully coupled problem in staggered form: a conventional mechanical problem by FFT-Galerkin under mixed macroscopic loading control, and a generalized Helmholtz equation for each nonlocal variable by preconditioned conjugate gradient. In three-dimensional particle-reinforced composites, \(128^3\) runs with roughly \(2\times 10^6\) voxels completed in a few hours on a moderate multicore workstation, and localized-band thickness scaled linearly with the regularization length \(\ell_j\) [2103.04770].

Two implementation studies highlight opposite ends of the software spectrum. A vectorized MATLAB code reported runtime reductions from \(145.3\) s to \(9.67\) s for a one-dimensional bar with \(1000\) elements, from \(3.40\) s to \(1.19\) s for a \(100\times100\) two-dimensional SEN problem, and from \(12.93\) s to \(5.67\) s for a \(100\times100\times5\) three-dimensional SEN benchmark, relative to non-vectorized baselines on a Ryzen 5 5600H laptop [2301.06503]. At the other end, a differentiable JAX-FEM framework combines gradient-enhanced damage with physics-augmented neural networks; the nonlocal field \(\phi\) obeys a PDE of the form
\[
-\mathrm{div}(g_c l^2 \nabla \phi)+g_c(\phi-\kappa)=0,
\]
and convergence studies showed virtually identical damage evolution across three meshes in uniform tension and notch-tip damage \(d\approx 0.79\) across meshes of \(241\), \(625\), and \(2350\) elements in a notched-plate simulation [2604.03411].

## 6. Applications, limitations, and evolving research directions

Applications of LGDM span quasi-brittle elastomers, metallic fracture, heterogeneous ductile failure, and anisotropic damage at finite strains. In elastomers, the method is used to represent large deformation, near incompressibility, diffuse chain scission, and crack propagation. Mousavi et al. explicitly compare LGDM with phase-field fracture and argue that both can deliver nearly identical macroscopic toughness and load response, while LGDM additionally carries a built-in model of diffuse network scission [2502.02822]. A later elastomer study extends this perspective by relating the gradient-enhanced damage zone to a fractocohesive length
\[
\ell_{fc}=G_c/W,
\]
defined as the width of the region where irreversible bond scission and energy dissipation occur; the reported implication is that the full damage-zone width is typically of order \(2\!-\!3\ell\) [2509.00313].

For metallic fracture, the gradient-damage framework is coupled to mechanism-based strain-gradient plasticity, introducing separate length scales for fracture regularization \(\ell_f\) and plastic hardening \(\ell_p\). The reported effect of \(\ell_p\) is twofold: it elevates crack-tip stresses and facilitates brittle-style failure for sharp defects, but delays localization and increases peak load for blunt defects [2108.04908]. In multiphase microstructures, phase-specific characteristic lengths \(\ell_j(x)\) allow nonlocal variables to diffuse only within their host phase, producing mesh-objective stress-strain responses where local Gurson or Lemaitre damage models diverge strongly with grid size [2103.04770]. For anisotropic damage, a reduced two-degree-of-freedom volumetric-deviatoric regularization reproduces the full six-degree-of-freedom response in post-peak softening and damage pattern, whereas a three-invariant alternative yields slightly wider damage bands and higher dissipation [2311.15918].

Several limitations recur. One is the longstanding tendency of classical gradient-damage models toward damage-zone broadening, which motivates the localizing modifications of both Mousavi et al. and Saji et al. [2502.02822] [2506.24099]. Another is numerical sensitivity under unstable crack growth, especially when artificial viscosity is introduced; in the elastomer comparison of phase-field and GED, the reported effect on \(J\) was particularly pronounced for GED [2408.05162]. A further source of confusion is conceptual rather than numerical: some “damage-based GED” formulations are described as identical to phase-field fracture and therefore as inheriting phase-field limitations, while other authors reserve LGDM for formulations that regularize auxiliary micromorphic drivers and derive damage algebraically [2502.02822].

Taken together, these works position LGDM as a broad research program rather than a single closed formalism. Its unifying objective is stable, mesh-independent softening with a finite process-zone width; its differentiating questions concern which internal variable should be regularized, how localization should be arrested at full damage, whether \(G_c\) is an input or an emergent output, and how closely the regularized damage field should be tied to micromechanical physics.

Source: https://www.emergentmind.com/topics/localizing-gradient-damage-method-lgdm