An efficient peridynamic model for damage and fracture analysis: the PD-GL framework with Gauss–Legendre–Lebedev quadrature

The computational cost of peridynamics has historically limited its use to small-scale problems, despite its conceptual appeal for fracture modeling and failure analysis. This work introduces a revised peridynamic integration paradigm that enables high-fidelity simulations using substantially fewer integration points while maintaining accuracy. The key observation is that conventional peridynamic volume integration can be reformulated as two separable components, radial and angular, allowing independent numerical treatment of each. Building on this decoupling, we propose a peridynamic integration framework that couples Gauss– quadrature in the radial direction with Lebedev spherical quadrature for angular integration (PD-GL). This reformulation yields three main advances: (i) systematic reduction of integration points via structured projection algorithms, (ii) removal of repeated neighbor-search overhead through spatial hashing, and (iii) retention of high accuracy even for small horizons δ = 2 Δ x , where conventional peridynamics typically deteriorates. Across a range of horizon sizes ( δ = 2 Δ x to 4 Δ x ), PD-GL achieves more than a sevenfold speedup relative to standard implementations, with the advantage increasing as the horizon grows. By integrating multiple established acceleration strategies into a single unified scheme, PD-GL addresses the primary computational bottlenecks of peridynamic simulations and facilitates practical, high-fidelity fracture and failure modeling in engineering applications.