SEARCH · Search NASA
Results for “peridynamic”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The Evolution of the Peridynamics Co-Authorship Network
We report peridynamics is a relatively new field in continuum mechanics that has developed over the past 20 years. This paper studies the evolution of collaborations in the field of peridynamics since its inception using social network analysis. For this purpose, we construct a network for each year from 2000 to 2019 describing co-authorship between scientists in peridynamics. In these networks, each node represents a scientist and each link connects two co-authoring scientists with a link weight representing the frequency and strength of their collaboration; each network as a whole can be thought of as a graph representation of the peridynamics community for the given year. By constructing a network for each year, we are able to analyze the evolution of the network in time and discuss the implications of this evolution for the peridynamics community. Our study demonstrates that the peridynamics community has been growing exponentially in size in recent years. Centrality metrics are also used to identify the most collaborative scientists in the community. Moreover, we compute link recommendations based on both elevating a scientist’s position in the network with respect to certain centrality metrics or closing structural holes in the network identified with persistent homology. We further extend the analysis to higher-order networks whose nodes represent groups of scientists in the community and whose links connect collaborating groups. In some sense, our work studies the past, present, and future of the peridynamics community.
PDMATLAB2D: A Peridynamics MATLAB Two-dimensional Code
PDMATLAB2D is a meshfree peridynamics implementation in MATLAB suitable for simulation of two-dimensional fracture problems. The purpose of this code is twofold. First, it provides an entry-level peridynamics computational tool for educational and training purposes. Second, it serves as an accessible and easily modifiable computational tool for peridynamics researchers who would like to adapt the code for a multitude of peridynamics simulation scenarios. The current version of the code implements a bond-based brittle elastic peridynamic model and a critical stretch criterion for bond breaking. However, the code is designed to be extendable for other peridynamic models and computational features. In this paper, we provide an overview of the code structure and functions with illustrative examples. In conclusion, due to the integrated computation and postprocessing MATLAB capabilities, PDMATLAB2D can serve as an effective testbed for testing new constitutive models and advanced numerical features for peridynamics computations.
Modeling Metallic Fuel using Peridynamics
Based on available modeling and simulation capabilities of peridynamics module in MOOSE framework for oxide fuel, the overall goal of this project is to further develop the peridynamics capabilities for modeling metallic fuel. It includes two major tasks: 1) develop validated scheme to handle the shape tensor singularity due to insufficient active neighbors of a material particle in the peridynamic correspondence model for fracture problems, and 2) develop failure modeling scheme including failure criterion for metallic fuels. Before the peridynamics can be applied to model metallic fuel, the formulation instability of the peridynamic correspondence model should be addressed. The PI first worked on developing new stabilization method to improve the performance of the peridynamic correspondence model and reduce the possibility of getting a singular shape tensor while applying the model for fracture problems. The new stabilization scheme uses bond-associated weight function rather than bond-associated horizon. Compared to bondassociated horizon stabilized method, this new stabilization scheme has better performance with improved prediction accuracy and reduced free surface effect. Using this newly developed stabilization, materials models from BISON can be directly used in peridynamics for metallic fuels, such as fission rate and burnup dependent creep and swell models. Publication of this work is under preparation.
Peridynamics and surrogate modeling of pressure-driven well stimulation
In this work we use the peridynamics theory of solid mechanics to simulate fracture in an annular rock domain subject to an in-situ stress and create surrogate models that predict the area of the resulting cracks. Peridynamics is a non-local formulation of continuum mechanics that naturally accommodates material discontinuities. Furthermore, unlike other fracture modeling techniques there is no need to provide information about the crack path. We utilize the peridynamics code Peridigm and take a two-stage approach to fracture modeling. First an implicit solve is performed to compute the in-situ stress state. We then execute an explicit solve where a pressure loading designed to emulate fluid-driven hydraulic fracture is applied at the borehole and transmitted to the pre-stressed rock. We present results from polynomial and single and multi-level Gaussian process surrogate models constructed from a sampling study of the peridynamics model. The surrogates predict crack area given a measure of the in-situ stress anisotropy and rise time and amplitude of the pressure loading. These surrogates take a minuscule fraction of peridynamics model's running time to evaluate and are a step towards enabling advanced optimization and uncertainty quantification workflows that require many model evaluations.
On the peridynamic effective force state and multiphase constitutive correspondence principle
This article concerns modeling unsaturated deformable porous media as an equivalent single-phase and single-force state peridynamic material through the effective force state. The balance equations of linear momentum and mass of unsaturated porous media are presented by defining relevant peridynamic states. The energy balance of unsaturated porous media is utilized to derive the effective force state for the solid skeleton that is an energy conjugate to the nonlocal deformation state of the solid, and the suction force state. Through an energy equivalence, a multiphase constitutive correspondence principle is built between classical unsaturated poromechanics and peridynamic unsaturated poromechanics. The multiphase correspondence principle provides a means to incorporate advanced constitutive models in classical unsaturated porous theory directly into unsaturated peridynamic poromechanics. Finally, numerical simulations of localized failure in unsaturated porous media under different matric suctions are presented to demonstrate the feasibility of modeling the mechanical behavior of such three-phase materials as an equivalent single-phase peridynamic material through the effective force state concept.
Validation Testing of a Peridynamic Impact Damage Model Using NASA's Micro-Particle Gun
Through a collaborative effort between the Virginia Commonwealth University and Raytheon, a peridynamic model for sand impact damage has been developed1-3. Model development has focused on simulating impacts of sand particles on ZnS traveling at velocities consistent with aircraft take-off and landing speeds. The model reproduces common features of impact damage including pit and radial cracks, and, under some conditions, lateral cracks. This study focuses on a preliminary validation exercise in which simulation results from the peridynamic model are compared to a limited experimental data set generated by NASA's recently developed micro-particle gun (MPG). The MPG facility measures the dimensions and incoming and rebound velocities of the impact particles. It also links each particle to a specific impact site and its associated damage. In this validation exercise parameters of the peridynamic model are adjusted to fit the experimentally observed pit diameter, average length of radial cracks and rebound velocities for 4 impacts of 300 μm glass beads on ZnS. Results indicate that a reasonable fit of these impact characteristics can be obtained by suitable adjustment of the peridynamic input parameters, demonstrating that the MPG can be used effectively as a validation tool for impact modeling and that the peridynamic sand impact model described herein possesses not only a qualitative but also a quantitative ability to simulate sand impact events.
Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures
Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.
Peridynamic Models for Random Media Found by Coarse Graining
Using coarse graining, the upscaled mechanical properties of a solid with small scale heterogeneities are derived. The method maps internal forces at the small scale onto peridynamic bond forces in the coarse grained mesh. These upscaled bond forces are used to calibrate a peridynamic material model with position-dependent parameters. These parameters incorporate mesoscale variations in the statistics of the small scale system. The upscaled peridynamic model can have a much coarser discretization than the original small scale model, allowing larger scale simulations to be performed efficiently. The convergence properties of the method are investigated for representative random microstructures. In conclusion, a bond breakage criterion for the upscaled peridynamic material model is also demonstrated.
An asymptotically compatible treatment of traction loading in linearly elastic peridynamic fracture
Meshfree discretizations of state-based peridynamic models are attractive due to their ability to naturally describe fracture of general materials. However, two factors conspire to prevent meshfree discretizations of state-based peridynamics from converging to corresponding local solutions as resolution is increased: quadrature error prevents an accurate prediction of bulk mechanics, and the lack of an explicit boundary representation presents challenges when applying traction loads. Herein, we develop a reformulation of the linear peridynamic solid (LPS) model to address these shortcomings, using improved meshfree quadrature, a reformulation of the nonlocal dilatation, and a consistent handling of the nonlocal traction condition to construct a model with rigorous accuracy guarantees. In particular, these improvements are designed to enforce discrete consistency in the presence of evolving fractures, whose a priori unknown location render consistent treatment difficult. In the absence of fracture, when a corresponding classical continuum mechanics model exists, our improvements provide asymptotically compatible convergence to corresponding local solutions, eliminating surface effects and issues with traction loading which have historically plagued peridynamic discretizations. When fracture occurs, our formulation automatically provides a sharp representation of the fracture surface by breaking bonds, avoiding the loss of mass. We provide rigorous error analysis and demonstrate convergence for a number of benchmarks, including manufactured solutions, free-surface, nonhomogeneous traction loading, and composite material problems. Finally, we validate simulations of brittle fracture against a recent experiment of dynamic crack branching in soda-lime glass, providing evidence that the scheme yields accurate predictions for practical engineering problems.
Bond-Based Peridynamics: a Tale of Two Poisson’s Ratios
This paper explores the restrictions imposed by bond-based peridynamics, particularly with respect to plane strain and plane stress models. We begin with a review of the derivations in Gerstle et al. (2005) wherein for isotropic materials a Poisson’s ratio restriction of $\frac{1}{4}$ for plane strain and $\frac{1}{3}$ for plane stress is deduced. Next, we show Cauchy’s relations are an intrinsic limitation of bond-based peridynamics and specialize these restrictions to plane strain and plane stress models. This generalizes the results from Gerstle et al. (2005) and demonstrates that the Poisson’s ratio restrictions described in Gerstle et al. (2005) are merely a consequence of Cauchy’s relations for isotropic materials. We conclude with a discussion of the validity of peridynamic plane strain and plane stress models formulated from two-dimensional bond-based peridynamic models.
An FEM-Based Peridynamic Model for Failure Analysis of Unidirectional Fiber-Reinforced Laminates
To predict the mixed damage modes of unidirectional fiber-reinforced polymer (FRP) laminates under dynamic loading, an FEM-based peridynamic model is introduced in this paper. Based on its geometric structure and material composition, a long fiber lamina is considered a transversely isotropic medium as a result of homogenization at the meso-scale. The laminated structure is modeled by stacking surface mesh layers with arbitrary fiber angles along the thickness direction. The peridynamic bonds between Gauss points connect the separated elements. These bonds are classified as inner-layer bonds and inter-layer bonds. To represent the anisotropy of a laminate, the micro-elastic modulus of the inner-layer and inter-layer bonds is calculated from the anisotropic engineering material constants separately. To capture complex failure behaviors of laminate structures, an empirical damage model is proposed for the tension/ compression breakage of peridynamic bonds. This damage model can control the in-plane and delamination failure process. Finally, benchmark tests are conducted to validate the elastic response of laminates under dynamic loading. In terms of damage analysis, the proposed model can capture the complex damage modes and resistive force of laminate structures.
Asymptotically compatible reproducing kernel collocation and meshfree integration for the peridynamic Navier equation
Here, we study reproducing kernel (RK) collocation method for peridynamic Navier equation. In the first part, we apply a linear RK approximation to both displacement and dilatation, and then back-substitute dilatation and solve the peridynamic Navier equation in a pure displacement form. The RK collocation scheme converges to the nonlocal limit for a fixed nonlocal interaction length and also to the local limit as nonlocal interactions vanish. The stability is shown by comparing the collocation scheme with the standard Galerkin scheme using Fourier analysis. In the second part, we apply the RK collocation to the quasi-discrete peridynamic Navier equation and show its convergence to the correct local limit when the ratio between the nonlocal length scale and the discretization parameter is fixed. The analysis is carried out on a special family of rectilinear Cartesian grids for the RK collocation method with a designated kernel with finite support. We assume the Lamé parameters satisfy λ ≥ μ to avoid extra assumptions on the nonlocal kernel. Finally, numerical experiments are conducted to validate the theoretical results.
Overall equilibrium in the coupling of peridynamics and classical continuum mechanics
Coupling peridynamics based computational tools with those using classical continuum mechanics can be very beneficial, because it can provide a means to generate a computational method that combines the efficiency of classical continuum mechanics with the capability to simulate crack propagation, typical of peridynamics. This paper presents an overlooked issue in this type of coupled computational methods: the lack of overall equilibrium. This can be the case even if the coupling strategy satisfies the usual numerical tests involving rigid body motions as well as uniform and linear strain distributions. We focus our investigation on the lack of overall equilibrium in an approach to couple peridynamics and classical continuum mechanics recently proposed by the authors. In our examples, the magnitude of the out-of-balance forces is a fraction of a percent of the applied forces, but it cannot be assumed to be a numerical round-off error. We show analytically and numerically that the main reason for the existence of out-of-balance forces is a lack of balance between the local and nonlocal tractions at the coupling interface. This usually results from the presence of high-order derivatives of displacements in the coupling zone.
Crack nucleation at forging flaws studied by non-local peridynamics simulations
We present a computational study and framework that allows us to study and understand the crack nucleation process from forging flaws. Forging flaws may be present in large steel rotor components commonly used for rotating power generation equipment including gas turbines, electrical generators, and steam turbines. The service life of these components is often limited by crack nucleation and subsequent growth from such forging flaws, which frequently exhibit themselves as non-metallic oxide inclusions. The fatigue crack growth process can be described by established engineering fracture mechanics methods. However, the initial crack nucleation process from a forging flaw is challenging for traditional engineering methods to quantify as it depends on the details of the flaw, including flaw morphology. We adopt the peridynamics method to describe and study this crack nucleation process. For a specific industrial gas turbine rotor steel, we present how we integrate and fit commonly known base material property data such as elastic properties, yield strength, and S-N curves, as well as fatigue crack growth data into a peridynamic model. The obtained model is then utilized in a series of high-performance two-dimensional peridynamic simulations to study the crack nucleation process from forging flaws for ambient and elevated temperatures in a rectangular simulation cell specimen. The simulations reveal an initial local nucleation at multiple small oxide inclusions followed by micro-crack propagation, arrest, coalescence, and eventual emergence of a dominant micro-crack that governs the crack nucleation process. The dependence on temperature and density of oxide inclusions of both the details of the microscopic processes and cycles to crack nucleation is also observed. The results are compared with fatigue experiments performed with specimens containing forging flaws of the same rotor steel.
A data-driven peridynamic continuum model for upscaling molecular dynamics
Nonlocal models, including peridynamics, often use integral operators that embed lengthscales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computational limitations in the length and time scales it can address. To combine the strengths of both methods and to obtain a coarse-grained, homogenized continuum model that efficiently and accurately captures materials’ behavior, we propose a learning framework to extract, from MD data, an optimal Linear Peridynamic Solid (LPS) model as a surrogate for MD displacements. To maximize the accuracy of the learnt model we allow the peridynamic influence function to be partially negative, while preserving the well-posedness of the resulting model. To achieve this, we provide sufficient well-posedness conditions for discretized LPS models with sign-changing influence functions and develop a constrained optimization algorithm that minimizes the equation residual while enforcing such solvability conditions. This framework guarantees that the resulting model is mathematically well-posed, physically consistent, and that it generalizes well to settings that are different from the ones used during training. We illustrate the efficacy of the proposed approach with several numerical tests for single layer graphene. Our two-dimensional tests show the robustness of the proposed algorithm on validation data sets that include thermal noise, different domain shapes and external loadings, and discretizations substantially different from the ones used for training.
Peridynamic Model for Single-Layer Graphene Obtained from Coarse Grained Bond Forces
An ordinary state-based peridynamic material model is proposed for single sheet graphene. The model is calibrated using coarse grained molecular dynamics simulations. The coarse graining method allows the dependence of bond force on bond length to be determined, including the horizon. The peridynamic model allows the horizon to be rescaled, providing a multiscale capability and allowing for substantial reductions in computational cost compared with molecular dynamics. The calibrated peridynamic model is compared to experimental data on the deflection and perforation of a graphene monolayer by an atomic force microscope probe.
Determination of Ballistic Limit for IM7/8552 Using Peridynamics
Significant testing is required to design and certify primary aircraft structure subject to High Energy Dynamic Impact (HEDI) events; current work under the NASA Advanced Composites Consortium (ACC) HEDI Project seeks to determine the state-of-the-art of dynamic fracture simulations for composite structures in these events. This paper discusses one of four Progressive Damage Analysis (PDA) methods selected for this project: peridynamics, through EMU implementation. A brief discussion of peridynamic theory is provided, followed by an outline of ballistic impact testing performed for model development and assessment. Detailed modeling approach and test-analysis correlation for a single open test case are presented, followed by the results of a series of blind predictions made prior to testing and test-analysis correlation performed with measured NASA test results. Specifically, we present simulation results for the ballistic limit (V50) of IM7/8552 composite panels ballistically tested with an impactor representative of a high-velocity fan-blade-out condition. In particular, force and displacement history and the damage state determined analytically are compared to measured results. Ultimately, peridynamics has the ability to predict damage patterns, impact force and deflections during a high energy dynamic impact event on composite panels of different layups using two different types of impactors. Blind predictions were promising and increased confidence in the model for impact simulation. There are open questions regarding the fidelity of the test fixture idealization in regards to stiffness and damping which will need to be addressed in future work.