# Convergence Analysis and Numerical Studies for Linearly Elastic Peridynamics with Dirichlet-Type Boundary Conditions

arXiv: https://arxiv.org/abs/2106.13878 · HTML: https://mlsanity.com/paper/2106.13878/

## Key facts

- Cited by 4 papers in the map (data through 14 Oct 2026).
- First posted on arXiv 25 Jun 2021.

## Built on by (most cited first, 4 of 4)

- [Peridynamic Neural Operators: A Data-Driven Nonlocal Constitutive Model for Complex Material Responses](https://mlsanity.com/paper/2401.06070/) (2024)
- [A data-driven peridynamic continuum model for upscaling molecular dynamics](https://mlsanity.com/paper/2108.04883/) (2021)
- [Heterogeneous Peridynamic Neural Operators: Discover Biotissue Constitutive Law and Microstructure From Digital Image Correlation Measurements](https://mlsanity.com/paper/2403.18597/) (2024)
- [Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures](https://mlsanity.com/paper/2301.04540/) (2023)

## Abstract (arXiv)

The nonlocal models of peridynamics have successfully predicted fractures and deformations for a variety of materials. In contrast to local mechanics, peridynamic boundary conditions must be defined on a finite volume region outside the body. Therefore, theoretical and numerical challenges arise in order to properly formulate Dirichlet-type nonlocal boundary conditions, while connecting them to the local counterparts. While a careless imposition of local boundary conditions leads to a smaller effective material stiffness close to the boundary and an artificial softening of the material, several strategies were proposed to avoid this unphysical surface effect. In this work, we study convergence of solutions to nonlocal state-based linear elastic model to their local counterparts as the interaction horizon vanishes, under different formulations and smoothness assumptions for nonlocal Dirichlet-type boundary conditions. Our results provide explicit rates of convergence that are sensitive to the compatibility of the nonlocal boundary data and the extension of the solution for the local model. In particular, under appropriate assumptions, constant extensions yield frac12 order convergence rates and linear extensions yield frac32 order convergence rates. With smooth extensions, these rates are improved to quadratic convergence. We illustrate the theory for any dimension d≥ 2 and numerically verify the convergence rates with a number of two dimensional benchmarks, including linear patch tests, manufactured solutions, and domains with curvilinear surfaces. Numerical results show a first order convergence for constant extensions and second order convergence for linear extensions, which suggests a possible room of improvement in the future convergence analysis.