Peridynamics (PD), similarly to some other nonlocal continuum theories, exhibits truncated interaction horizons near free surfaces, cracks, and voids in bounded domains. This loss of neighbors causes artificial surface effects and inconsistencies in the derivative and energy operators that persist under discretization refinement, limiting PD accu- racy and robustness. First, the analysis and numerical isolation of the effects ofhorizon truncation on these operators are carried out, showing how it induces surface softening in bond-based and stiffening/softening in state-based PD formulations. Then, a purely PD-based surface correction is proposed to restore full-horizon behavior while preserv- ing the original horizon radius and bond topology. Each node is assigned a scalar influence weight, and every bond is scaled by the average of the endpoint weights to maintain symmetry. The nodal influence weights are computed in a preprocessing optimization step enforcing agreement between truncated- and full-horizon derivative and energy operators. Unlike existing approaches, the proposed method does not rely on ghost particles, variable horizons, or reference solutions from classical continuum mechanics (CCM). Benchmarks in one, two, and three dimensions, including a dynamic brittle fracture test, show that the optimized nodal influence weights reduce boundary-induced surface effects, improve energy consistency, andaccelerate convergence of PD solutions. Analytical evaluation of the calibration targets over full-horizons also improves interior volume integration accuracy. The proposed scheme is sim- ple to implement in existing PD codes, incurs only a modest preprocessing cost, performs robustly for flat and curved boundaries on both regular and irregular point sets, produces reliable crack paths, and applies to both bond-based and ordinary state-based PD formulations. An open-source Python implementation, perifit, is provided.
Surface and boundary corrections in peridynamics using optimized nodal influence weights / Shojaei, A., Ongaro, G., Hermann, A., Silling, S.A., Trovalusci, P., Cyron, C.J.. - In: INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING. - ISSN 0029-5981. - 127:14(2026). [10.1002/nme.70385]
Surface and boundary corrections in peridynamics using optimized nodal influence weights
Ongaro, Greta
;Trovalusci, Patrizia;
2026
Abstract
Peridynamics (PD), similarly to some other nonlocal continuum theories, exhibits truncated interaction horizons near free surfaces, cracks, and voids in bounded domains. This loss of neighbors causes artificial surface effects and inconsistencies in the derivative and energy operators that persist under discretization refinement, limiting PD accu- racy and robustness. First, the analysis and numerical isolation of the effects ofhorizon truncation on these operators are carried out, showing how it induces surface softening in bond-based and stiffening/softening in state-based PD formulations. Then, a purely PD-based surface correction is proposed to restore full-horizon behavior while preserv- ing the original horizon radius and bond topology. Each node is assigned a scalar influence weight, and every bond is scaled by the average of the endpoint weights to maintain symmetry. The nodal influence weights are computed in a preprocessing optimization step enforcing agreement between truncated- and full-horizon derivative and energy operators. Unlike existing approaches, the proposed method does not rely on ghost particles, variable horizons, or reference solutions from classical continuum mechanics (CCM). Benchmarks in one, two, and three dimensions, including a dynamic brittle fracture test, show that the optimized nodal influence weights reduce boundary-induced surface effects, improve energy consistency, andaccelerate convergence of PD solutions. Analytical evaluation of the calibration targets over full-horizons also improves interior volume integration accuracy. The proposed scheme is sim- ple to implement in existing PD codes, incurs only a modest preprocessing cost, performs robustly for flat and curved boundaries on both regular and irregular point sets, produces reliable crack paths, and applies to both bond-based and ordinary state-based PD formulations. An open-source Python implementation, perifit, is provided.I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.


