Publications / SAND Report

Strong Local-Nonlocal Coupling for Integrated Fracture Modeling

Littlewood, David J.; Silling, Stewart A.; Mitchell, John A.; Seleson, Pablo D.; Bond, Stephen D.; Parks, Michael L.; Turner, Daniel Z.; Burnett, Damon J.; Ostien, Jakob O.; Gunzburger, Max G.

Peridynamics, a nonlocal extension of continuum mechanics, is unique in its ability to capture pervasive material failure. Its use in the majority of system-level analyses carried out at Sandia, however, is severely limited, due in large part to computational expense and the challenge posed by the imposition of nonlocal boundary conditions. Combined analyses in which peridynamics is em- ployed only in regions susceptible to material failure are therefore highly desirable, yet available coupling strategies have remained severely limited. This report is a summary of the Laboratory Directed Research and Development (LDRD) project "Strong Local-Nonlocal Coupling for Inte- grated Fracture Modeling," completed within the Computing and Information Sciences (CIS) In- vestment Area at Sandia National Laboratories. A number of challenges inherent to coupling local and nonlocal models are addressed. A primary result is the extension of peridynamics to facilitate a variable nonlocal length scale. This approach, termed the peridynamic partial stress, can greatly reduce the mathematical incompatibility between local and nonlocal equations through reduction of the peridynamic horizon in the vicinity of a model interface. A second result is the formulation of a blending-based coupling approach that may be applied either as the primary coupling strategy, or in combination with the peridynamic partial stress. This blending-based approach is distinct from general blending methods, such as the Arlequin approach, in that it is specific to the coupling of peridynamics and classical continuum mechanics. Facilitating the coupling of peridynamics and classical continuum mechanics has also required innovations aimed directly at peridynamic models. Specifically, the properties of peridynamic constitutive models near domain boundaries and shortcomings in available discretization strategies have been addressed. The results are a class of position-aware peridynamic constitutive laws for dramatically improved consistency at domain boundaries, and an enhancement to the meshfree discretization applied to peridynamic models that removes irregularities at the limit of the nonlocal length scale and dramatically improves conver- gence behavior. Finally, a novel approach for modeling ductile failure has been developed, moti- vated by the desire to apply coupled local-nonlocal models to a wide variety of materials, including ductile metals, which have received minimal attention in the peridynamic literature. Software im- plementation of the partial-stress coupling strategy, the position-aware peridynamic constitutive models, and the strategies for improving the convergence behavior of peridynamic models was completed within the Peridigm and Albany codes, developed at Sandia National Laboratories and made publicly available under the open-source 3-clause BSD license.