An asymptotically compatible treatment of traction loading in linearly elastic peridynamic fracture
Introduction
Peridynamics provides a description of continuum mechanics in terms of integral operators rather than classical differential operators [1], [2], [3], [4], [5], [6], [7]. These nonlocal models are defined in terms of a lengthscale , referred to as a horizon, which denotes the extent of nonlocal interaction. The nonlocal viewpoint allows a natural description of processes requiring reduced regularity in the relevant solution, such as fracture mechanics [8], [9]. An important feature of such models is that when classical continuum models still apply, they revert back to classical continuum models as . Discretizations which preserve this limit under refinement are termed asymptotically compatible (AC) [10], and there has been significant work in recent years toward establishing such discretizations — for an incomplete list see [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20]. Broadly, strategies either involve adopting traditional finite element shape functions and carefully performing geometric calculations to integrate over relevant horizon/element subdomains, or adopt a strong-form meshfree discretization where particles are associated with abstract measure. The former is more amenable to mathematical analysis due to a better variational setting, while the latter is simple to implement and generally faster [21], [22]. In this paper we pursue the meshfree viewpoint.
For fracture mechanics problems one often refines both and at the same rate under so-called M-convergence, , for [23]. In this setting, one obtains banded stiffness matrices allowing scalable implementations. Typically in the literature a scheme is termed AC if it recovers the solution in both the finite and M-convergence limit — in this work we abuse the definition slightly and only require the M-convergence case for asymptotic compatibility as the relevant limit for problems with a corresponding local limit. This AC property is only one necessary ingredient in achieving a convergent simulation, and our recent work focused upon establishing convergence in this setting for boundary value problems [18], [19]. To achieve similar convergence for problems involving fracture, one must also consider the interplay between consistency of quadrature for discrete operators and the imposition of traction loads as fracture surfaces open up and evolve [24]. For peridynamic fracture problems where the free surface evolves implicitly via the breaking of bonds [17], [25], one lacks an explicit boundary representation over the course of a simulation. In addition to providing challenges regarding accurate imposition of traction loads, the breaking of bonds also renders higher-order numerical quadrature inaccurate, as consistent AC quadrature weights are typically derived in the absence of damage.
Our goal is to provide a comprehensive treatment of fracture, nonlocal quadrature, and traction loading which is able to perform more accurate state-based peridynamic fracture simulations free of spurious surface effects. In particular, when no fracture occurs and therefore the classical continuum theory applies, the formulation should preserve the AC limit under M-convergence. When fracture occurs, the formulation should be able to capture the material damage and the evolving fracture surfaces via bond breaking. This practically means that one is able to incorporate all of the necessary ingredients to perform non-trivial simulations of fracture mechanics while maintaining a scalable implementation and guaranteeing convergence. Such a capability is elusive in the peridynamic literature; while peridynamics has been shown to provide a powerful modeling platform for a broad range of applications [26], [27], the development of efficient discretizations with rigorous underpinnings has lagged behind until the last few years.
The challenge in incorporating traction loading into a peridynamic framework stems from the fact that, in contrast to local mechanics, peridynamic boundary conditions must be defined on a finite volume region outside the surface [9], [20], [28]. Theoretical and numerical challenges arise in how to mathematically impose nonhomogeneous Neumann boundary conditions properly in the nonlocal model. In peridynamic models, careless imposition of traction loads leads to a smaller effective material stiffness close to the boundary, since the integral on those material points is over a smaller region. Therefore, an unphysical strain energy concentration is induced, leading in turn to an artificial softening of the material near the boundary. Such undesirable phenomena are referred to in the literature as a “surface” or “skin” effects [29], [30]. We propose a novel treatment of nonlocal traction-type boundary conditions which avoid the surface effect by designing a loading aimed to recover the corresponding local traction boundary condition as . The approach requires no explicit representation of the boundary, imposing the traction volumetrically using the same information that would be available during a traditional meshfree bond-based peridynamics simulation. Although the Neumann-constrained nonlocal problem and its AC limit were investigated in nonlocal diffusion models [18], [19], [20], [28], [31], [32], to the authors’ best knowledge, the development of AC peridynamic formulations with traction-type boundary conditions remains restricted to weak formulations, simple traction loadings and/or simple geometries. Several modeling and numerical approaches have been proposed to correct the surface effect [7], [33], [34], [35], [36], [37], [38], [39] but mostly restricted to free surfaces. For nonzero loadings, the tractions are often applied as prescribed body forces through a layer of finite thickness at the material boundary [27], [33], [40], as a surface integral through a weak form [41], or by modifying the nonlocal operator through eigenvalues analysis [42]. Therefore, developing an AC meshfree discretization method for peridynamics which is capable to handle nonhomogeneous traction loadings on complex boundaries is critical for the general practice of peridynamics in realistic engineering applications.
We consider the linear peridynamic solid (LPS) model [43] as a prototypical state-based model appropriate for brittle fracture. The LPS model may be interpreted as a nonlocal generalization of the mixed form of linear elasticity, evolving both displacements and a dilatation. We will show that consistent treatment evolving traction loading will require a modification to the definition of dilatation to guarantee consistency in the presence of fractures; conceptually this corresponds to the fact that dilatation is a kinematic variable without associated boundary conditions, and should be estimated consistently independently of whether a fracture is occurring in the vicinity of a given point. Based on the modified nonlocal dilatation, we further propose a new nonlocal generalization of classical traction loads in the LPS model. Particularly, we convert the local traction loads to a correction term in the momentum balance equation, which provides an estimate for the nonlocal interactions of each material point with points outside the domain. Based on this traction-type boundary condition, a meshfree formulation is developed for the LPS model based on the optimization-based quadrature rule [17], which preserves the AC limit under M-convergence and naturally represents the evolving free surfaces in dynamic fracture problems. We note that asymptotic compatibility is not well-defined for dynamic fracture, as there is no known corresponding local theory for peridynamics with bond breaking.1 However, our modified LPS formulation preserves the AC limit for the linear elastic model with traction loading on the evolving fracture surfaces. This fact, together with the consistent discretization introduced here, provides the opportunity for efficient and accurate peridynamic fracture simulations.
We remark that the paper is organized to first establish the rigorous mathematical underpinnings of the approach, while the second half focuses on a more engineering-oriented exploration of its application. Readers with more applied interests may skip many of the proofs in the work without issue. The work is organized as follows. We recall first the linear peridynamic solid (LPS) model definition in Section 2. In Section 3.1, we introduce a novel approach to apply classical traction loads on the LPS model. After establishing the continuous limits of the scheme, we next pursue a consistent discretization. In Section 4 we introduce meshfree quadrature which preserves asymptotic compatibility in the limit, and establish the discrete scheme for boundary value problems in the absence of fracture. We proceed to investigate a number of two-dimensional statics problems with analytic solutions for the local limit in Section 5. These test cases include: linear patch tests (Section 5.1); manufactured local limits to illustrate asymptotic convergence rates (Section 5.2); homogeneous materials with free-surfaces or non-zero traction loading on curvilinear surfaces (Section 5.3); composite materials with internal interfaces (Section 5.4). In Section 6, we further extend the proposed formulation to handle dynamic brittle fracture, and provide preliminary validation results by comparing our numerical results with available numerical simulations and experimental measurements on three benchmark problems. Section 7 summarizes our findings and discusses future research.
Section snippets
A linear state-based peridynamic model
We consider the state-based linear peridynamic solid (LPS) model in a body occupying the domain , or . Let be the nonlocal dilatation, generalizing the local divergence of displacement, and denote a positive radial function with compactly supported on the -ball . The momentum balance and nonlocal dilatation are then given by the following,
Neumann and mixed-type constraint problems
We now consider a state-based peridynamic problem with general mixed boundary conditions: and . Here and are both 1D curves. We denote the regions near the boundary as Note that to apply the nonlocal Dirichlet-type boundary condition, is required in a layer with non-zero volume outside , while the proposed traction load is applied as a Neumann boundary condition on the sharp
Optimization-based meshfree quadrature rules
In this section, we introduce a strong-form particle discretization of the state-based peridynamics introduced above. Discretizing the whole interaction region by a collection of points , we aim to solve for the displacement and the nonlocal dilatation on all . We first characterize the distribution of collocation points as follows. Recall the definitions [52] of fill distance and separation distance
Numerical verification and asymptotic compatibility
In this section we numerically verify the approach by investigating accuracy when recovering analytic solutions in the M-convergence limit with mixed boundary conditions. We consider: linear patch tests, smooth manufactured solutions, analytical solutions to curvilinear surface loading problems, and analytical solutions to linearly elastic composites. For each case, we consider various combinations of Dirichlet and traction-type boundary conditions, exploring also the effect of reduced
Fracture dynamics for brittle fracture experiments
The previous sections have established the ability of the scheme to recover local solutions of boundary value problems in elasticity governed by traction loadings and ensured that the breaking bonds treatment does not impair the AC convergence of the quadrature treatment. Of course, the main appeal of peridynamic discretizations is to handle fracture problems, and we devote the remainder of the paper to demonstrating how the scheme prescribed previously adapts to practical engineering settings,
Conclusion and future work
Peridynamics presents a flexible framework for modeling fracture mechanics. In particular, bond-based fracture models admit a sharp representation of fracture surfaces while avoiding the loss of mass associated with damage models and element death [73]. This flexibility comes with a cost however, as the free-surface introduced during fracture compounds traditional challenges in peridynamic models related to nonlocal boundary conditions. This work has presented a complete workflow demonstrating
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgments
Sandia National Laboratories is a multi-mission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC., a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. This paper describes objective technical results and analysis. Any subjective views or opinions that might be expressed in the paper do not necessarily represent the views of the U.S.
References (73)
Reformulation of elasticity theory for discontinuities and long-range forces
J. Mech. Phys. Solids
(2000)- et al.
Implementing peridynamics within a molecular dynamics code
Comput. Phys. Comm.
(2008) - et al.
A reproducing kernel enhanced approach for peridynamic solutions
Comput. Methods Appl. Mech. Engrg.
(2018) - et al.
Convergence studies in meshfree peridynamic simulations
Comput. Math. Appl.
(2016) - et al.
An asymptotically compatible meshfree quadrature rule for nonlocal problems with applications to peridynamics
Comput. Methods Appl. Mech. Engrg.
(2019) - et al.
An asymptotically compatible formulation for local-to-nonlocal coupling problems without overlapping regions
Comput. Methods Appl. Mech. Engrg.
(2020) - et al.
Nonlocal diffusion and peridynamic models with Neumann type constraints and their numerical approximations
Appl. Math. Comput.
(2017) - et al.
A meshfree method based on the peridynamic model of solid mechanics
Comput. Struct.
(2005) - et al.
Characteristics of dynamic brittle fracture captured with peridynamics
Eng. Fract. Mech.
(2011) - et al.
Integral approximations to classical diffusion and smoothed particle hydrodynamics
Comput. Methods Appl. Mech. Engrg.
(2015)
Peridynamics via finite element analysis
Finite Elem. Anal. Des.
Weak form of peridynamics for nonlocal essential and natural boundary conditions
Comput. Methods Appl. Mech. Engrg.
An improved SPH method: Towards higher order convergence
J. Comput. Phys.
Minimal positive stencils in meshfree finite difference methods for the Poisson equation
Comput. Methods Appl. Mech. Engrg.
A meshfree method based on the peridynamic model of solid mechanics
Comput. Struct.
An energy-based coupling approach to nonlocal interface problems
Comput. & Fluids
Numerical simulation of crack curving and branching in brittle materials under dynamic loads using the extended non-ordinary state-based peridynamics
Eur. J. Mech. A Solids
Visualization of fracture progression in peridynamics
Comput. Graph.
A partitioned coupling framework for peridynamics and classical theory: analysis and simulations
Comput. Methods Appl. Mech. Engrg.
Dynamic fracture of soda-lime glass: A full-field optical investigation of crack initiation, propagation and branching
J. Mech. Phys. Solids
Some computational capabilities for nonlinear finite element analysis
Nucl. Eng. Des.
Peridynamics as an upscaling of molecular dynamics
Multiscale Model. Simul.
A Continuum Theory with Long-Range Forces for Solids
Analysis and numerical approximation of an integro-differential equation modeling non-local effects in linear elasticity
Math. Mech. Solids
Mathematical analysis for the peridynamic nonlocal continuum theory
ESAIM Math. Model. Numer. Anal.
Handbook of Peridynamic Modeling
Nonlocal integral formulations of plasticity and damage: survey of progress
J. Eng. Mech.
A nonlocal vector calculus, nonlocal volume-constrained problems, and nonlocal balance laws
Math. Models Methods Appl. Sci.
Asymptotically compatible schemes and applications to robust discretization of nonlocal models
SIAM J. Numer. Anal.
Numerical methods for nonlocal and fractional models
Asymptotically compatible reproducing kernel collocation and meshfree integration for nonlocal diffusion
Generalized reproducing kernel peridynamics: unification of local and non-local meshfree methods, non-local derivative operations, and an arbitrary-order state-based peridynamic formulation
Comput. Part. Mech.
Local limits and asymptotically compatible discretizations
An asymptotically compatible approach for Neumann-type boundary condition on nonlocal problems
ESAIM Math. Model. Numer. Anal.
A meshfree unification: reproducing kernel peridynamics
Comput. Mech.
Convergence, adaptive refinement, and scaling in 1D peridynamics
Internat. J. Numer. Methods Engrg.
Cited by (25)
Phase-field modeling of anisotropic crack propagation based on higher-order nonlocal operator theory
2024, International Journal of Solids and StructuresAn explicit updated Lagrangian fragile points method for dynamic fracture simulations
2024, Engineering Fracture MechanicsHigher-order nonlocal operator theory for phase-field modeling of ductile fracture in elasto-plastic materials
2023, Computer Methods in Applied Mechanics and EngineeringMeshless numerical solution for nonlocal integral differentiation equation with application in peridynamics
2022, Engineering Analysis with Boundary ElementsAn asymptotically compatible probabilistic collocation method for randomly heterogeneous nonlocal problems
2022, Journal of Computational PhysicsA meshfree peridynamic model for brittle fracture in randomly heterogeneous materials
2022, Computer Methods in Applied Mechanics and EngineeringCitation Excerpt :The LPS model may be interpreted as a nonlocal generalization of the mixed form of linear elasticity, evolving both displacements and a dilatation. We begin with a review of the deterministic LPS model for heterogeneous materials [34] in Section 2.1, then extend the formulation to the stochastic LPS problem with random parameters in Section 2.2. Finally, we discuss the treatment of material fracture, including the damage criteria and the handling of free surfaces created by the evolving fractures, in Section 2.3.