Search

The peridynamic theory of solid mechanics is applied to modeling the deformation and fracture of micrometer-sized particles made of organic crystalline material. A new peridynamic material model is proposed to reproduce the elastic–plastic response, creep, and fracture that are observed in experiments. The model is implemented in a three-dimensional, meshless Lagrangian simulation code. In the small deformation, elastic regime, the model agrees well with classical Hertzian contact analysis for a sphere compressed between rigid plates. Under higher load, material and geometrical nonlinearity is predicted, leading to fracture. The material parameters for the energetic material CL-20 are evaluated from nanoindentation test data on the cyclic compression and failure of micrometer-sized grains.

The propagation of a wave pulse due to low-speed impact on a one-dimensional, heterogeneous bar is studied. Due to the dispersive character of the medium, the pulse attenuates as it propagates. This attenuation is studied over propagation distances that are much longer than the size of the microstructure. A homogenized peridynamic material model can be calibrated to reproduce the attenuation and spreading of the wave. The calibration consists of matching the dispersion curve for the heterogeneous material near the limit of long wavelengths. It is demonstrated that the peridynamic method reproduces the attenuation of wave pulses predicted by an exact microstructural model over large propagation distances.

The dynamic behavior of an elastic peridynamic material with a nonconvex bond potential is studied. In spite of the material’s inherently unstable nature, initial value problems can be solved using essentially the same techniques as with conventional materials, both analytically and numerically. In a suitably constructed material model, small perturbations grow exponentially over time until the material fails. The time for this growth is computed explicitly for a stretching bar that passes from the stable to the unstable phase of the material model. This time to failure represents an incubation time for the nucleation of a crack. The finiteness of the failure time in effect creates a rate dependence in the failure properties of the material. Thus, the unstable nature of the elastic material leads to a rate effect even though it does not contain any terms that explicitly include a strain rate dependence.

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.

The effect of spatial nonlocality on the decay of waves in a dissipative material is investigated. The propagation and decay of waves in a one-dimensional, viscoelastic peridynamic medium is analyzed. Both the elastic and damping terms in the material model are nonlocal. Waves produced by a source with constant amplitude applied at one end of a semi-infinite bar decay exponentially with distance from the source. The model predicts a cutoff frequency that is influenced by the nonlocal parameters. A method for computing the attenuation coefficient explicitly as a function of material properties and source frequency is presented. The theoretical results are compared with direct numerical simulations in the time domain. The relationship between the attenuation coefficient and the group velocity is derived. It is shown that in the limit of long waves (or small peridynamic horizon), Stokes’ law of sound attenuation is recovered.

We show that machine learning can improve the accuracy of simulations of stress waves in one-dimensional composite materials. We propose a data-driven technique to learn nonlocal constitutive laws for stress wave propagation models. The method is an optimization-based technique in which the nonlocal kernel function is approximated via Bernstein polynomials. The kernel, including both its functional form and parameters, is derived so that when used in a nonlocal solver, it generates solutions that closely match high-fidelity data. The optimal kernel therefore acts as a homogenized nonlocal continuum model that accurately reproduces wave motion in a smaller-scale, more detailed model that can include multiple materials. We apply this technique to wave propagation within a heterogeneous bar with a periodic microstructure. Several one-dimensional numerical tests illustrate the accuracy of our algorithm. The optimal kernel is demonstrated to reproduce high-fidelity data for a composite material in applications that are substantially different from the problems used as training data.

Conventional engineering methods oftentimes have challenges in the quantification of crack nucleation processes from manufacturing defects that are relevant for engineering component lifing. We present peridynamic simulation framework for the description of the crack nucleation process from cluster of non-metallic inclusions in aluminum alloys. Our non-local simulation framework characterizes crack nucleation process as multiple micro-crack nucleation events from individual inclusions, and eventually one micro-crack dominates. We define individual stages of the crack nucleation process, i.e., nucleation, micro-crack, technical, and crack initiation, that allow a quantification and meta model development of both the individual stages and the entire crack nucleation process.

Most previous development of the peridynamic theory has assumed a Lagrangian formulation, in which the material model refers to an undeformed reference configuration. In the present work, an Eulerian form of material modeling is developed, in which bond forces depend only on the positions of material points in the deformed configuration. The formulation is consistent with the thermodynamic form of the peridynamic model and is derivable from a suitable expression for the free energy of a material. It is shown that the resulting formulation of peridynamic material models can be used to simulate strong shock waves and fluid response in which very large deformations make the Lagrangian form unsuitable. The Eulerian capability is demonstrated in numerical simulations of ejecta from a wavy free surface on a metal subjected to strong shock wave loading. The Eulerian and Lagrangian contributions to bond force can be combined in a single material model, allowing strength and fracture under tensile or shear loading to be modeled consistently with high compressive stresses. This capability is demonstrated in numerical simulation of bird strike against an aircraft, in which both tensile fracture and high pressure response are important.

Peridynamic correspondence material models provide a way to combine a material model from the local theory with the inherent capabilities of peridynamics to model long-range forces and fracture. However, correspondence models in a typical particle discretization suffer from zero-energy mode instability. These instabilities are shown here to be an aspect of material stability. A stability condition is derived for state-based materials starting from the requirement of potential energy minimization. It is shown that all correspondence materials fail this stability condition due to zero-energy deformation modes of the family. To eliminate these modes, a term is added to the correspondence strain energy density that resists deviations from a uniform deformation. The resulting material model satisfies the stability condition while effectively leaving the stress tensor unchanged. Computational examples demonstrate the effectiveness of the modified material model in avoiding zero-energy mode instability in a peridynamic particle code.

The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. The computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.

The propagation of large amplitude nonlinear waves in a peridynamic solid is analyzed. With an elastic material model that hardens in compression, sufficiently large wave pulses propagate as solitary waves whose velocity can far exceed the linear wave speed. In spite of their large velocity and amplitude, these waves leave the material they pass through with no net change in velocity and stress. They are nondissipative and nondispersive, and they travel unchanged over large distances. An approximate solution for solitary waves is derived that reproduces the main features of these waves observed in computational simulations. It is demonstrated by numerical studies that the waves interact only weakly with each other when they collide. Wavetrains composed of many non-interacting solitary waves are found to form and propagate under certain boundary and initial conditions.

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. 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.

A notion of material homogeneity is proposed for peridynamic bodies with vari- able horizon but constant bulk properties. A relation is derived that scales the force state according to the position-dependent horizon while keeping the bulk properties un- changed. Using this scaling relation, if the horizon depends on position, artifacts called ghost forces may arise in a body under homogeneous deformation. These artifacts de- pend on the second derivative of horizon and can be reduced by use of a modified equilibrium equation using a new quantity called the partial stress . Bodies with piece- wise constant horizon can be modeled without ghost forces by using a technique called a splice between the regions. As a limiting case of zero horizon, both partial stress and splice techniques can be used to achieve local-nonlocal coupling. Computational examples, including dynamic fracture in a one-dimensional model with local-nonlocal coupling, illustrate the methods.

A position-aware linear solid (PALS) peridynamic constitutive model is proposed for isotropic elastic solids. The PALS model addresses problems that arise, in ordinary peridynamic material models such as the linear peridynamic solid (LPS), due to incomplete neighborhoods near the surface of a body. We improved model behavior in the vicinity of free surfaces through the application of two influence functions that correspond, respectively, to the volumetric and deviatoric parts of the deformation. Furthermore, the model is position-aware in that the influence functions vary over the body and reflect the proximity of each material point to free surfaces. Demonstration calculations on simple benchmark problems show a sharp reduction in error relative to the LPS model.

A simple demonstration of nonlocality in a heterogeneous material is presented. By analysis of the microscale deformation of a two-component layered medium, it is shown that nonlocal interactions necessarily appear in a homogenized model of the system. Explicit expressions for the nonlocal forces are determined. The way these nonlocal forces appear in various nonlocal elasticity theories is derived. The length scales that emerge involve the constituent material properties as well as their geometrical dimensions. A peridynamic material model for the smoothed displacement field is derived. It is demonstrated by comparison with experimental data that the incorporation of nonlocality in modeling dramatically improves the prediction of the stress concentration in an open hole tension test on a composite plate.

The peridynamic theory of solid mechanics provides a natural framework for modeling constitutive response and simulating dynamic crack propagation, pervasive damage, and fragmentation. In the case of a fragmenting body, the principal quantities of interest include the number of fragments, and the masses and velocities of the fragments. We present a method for identifying individual fragments in a peridynamic simulation. We restrict ourselves to the meshfree approach of Silling and Askari, in which nodal volumes are used to discretize the computational domain. Nodal volumes, which are connected by peridynamic bonds, may separate as a result of material damage and form groups that represent fragments. Nodes within each fragment have similar velocities and their collective motion resembles that of a rigid body. The identification of fragments is achieved through inspection of the peridynamic bonds, established at the onset of the simulation, and the evolving damage value associated with each bond. An iterative approach allows for the identification of isolated groups of nodal volumes by traversing the network of bonds present in a body. The process of identifying fragments may be carried out at specified times during the simulation, revealing the progression of damage and the creation of fragments. Incorporating the fragment identification algorithm directly within the simulation code avoids the need to write bond data to disk, which is often prohibitively expensive. Results are recorded using fragment identification numbers. The identification number for each fragment is stored at each node within the fragment and written to disk, allowing for any number of post-processing operations, for example the construction of cumulative distribution functions for quantities of interest. Care is taken with regard to very small clusters of isolated nodes, including individual nodes for which all bonds have failed. Small clusters of nodes may be treated as tiny fragments, or may be omitted from the fragment identification process. The fragment identification algorithm is demonstrated using the Sierra/SolidMechanics analysis code. It is applied to a simulation of pervasive damage resulting from a spherical projectile impacting a brittle disk, and to a simulation of fragmentation of an expanding ductile ring.

This work presents a rigorous derivation for the formulation of the J-integral in bond-based peridynamics using the crack infinitesimal virtual extension approach. We give a detailed description of an algorithm for computing this nonlocal version of the J-integral. We present convergence studies (m-convergence and δ-convergence) for two different geometries: a single edge-notch configuration and a double edge-notch sample. We compare the results with results based on the classical J-integral and obtained from FEM calculations that employ special elements near the crack tip. We identify the size of the nonlocal region for which the peridynamic J-integral value is near the classical FEM solutions. We discuss how the boundary conditions and the peridynamic “skin effect” may influence the peridynamic J-integral value. We also observe, computationally, the path-independence of the peridynamic J-integral.

This study presents an analysis approach based on a merger of the finite element method and the peridynamic theory. Its validity is established through qualitative and quantitative comparisons against the test results for a stiffened composite curved panel with a central slot under combined internal pressure and axial tension. The predicted initial and final failure loads, as well as the final failure modes, are in close agreement with the experimental observations. This approach demonstrates the capability of the PD approach to assess the durability of complex composite structures.

The design of an aerospace structure using an off-the-shelf composite would involve increasing the gauge thickness until all the design requirements are met. This can lead to an inefficient design, because excess margins will exist for all properties except the one that determines the gauge. The design of a material can be made practical by creating a hybrid composite consisting of two or more types of fibers or resins, each embellishing a particular trait or function to the material. This paper investigates both high- and low-energy hail impact against a toughened-epoxy, intermediate-modulus, carbon-fiber composite using both experimental and analytical means. The effect of introducing ply-level hybridization by substituting up to 20% of the plies with glass-reinforced plies is considered. It is found that delamination can be reduced by this hybridization, but the benefits are dependent on the impact energy and the test conditions. A computational model based on the peridynamic theory of solid mechanics ...

A condition for the emergence of a discontinuity in an elastic peridynamic body is pro- posed, resulting in a material stability condition for crack nucleation. The condition is derived by deter- mining whether a small discontinuity in displace- ment, superposed on a possibly large deformation, grows over time. Stability is shown to be determined by the sign of the eigenvalues of a tensor field that dependsonlyonthelinearizedmaterialproperties.This condition for nucleation of a discontinuity in dis- placement can be interpreted in terms of the dynamic stability of plane waves with very short wavelength. A numerical example illustrates that cracks in a peri- dynamic body form spontaneously as the body is loaded.

In this paper, we compare small deformations in an infinite linear elastic body due to the presence of point loads within the classical, local formulation to the corresponding deformations in the peridynamic, non-local formulation. Owing to the linearity of the problem, the response to a point load can be used to obtain the response to general body force loading functions by superposition. Using Laplace and Fourier transforms, we thus obtain an integral representation for the three-dimensional peridynamic solution with the help of Green’s functions. We illustrate this new theoretical result by dynamic and static examples in one and three dimensions. In addition to this main result, we also derive the non-local three-dimensional jump conditions, as well as the weak formulation of peridynamics together with the associated finite element discretization.

Peridynamics is a continuum reformulation of the standard theory of solid mechanics. Unlike the partial differential equations of the standard theory, the basic equations of peridynamics are applicable even when cracks and other singularities appear in the deformation field. The assumptions in the original peridynamic theory resulted in severe restrictions on the types of material response that could be modeled, including a limitation on the Poisson ratio. Recent theoretical developments have shown promise for overcoming these limitations, but have not previously incorporated rate dependence and have not been demonstrated in realistic applications. In this paper, a new method for implementing a rate-dependent plastic material within a peridynamic numerical model is proposed and demonstrated. The resulting material model implementation is fitted to rate-dependent test data on 6061-T6 aluminum alloy. It is shown that with this material model, the peridynamic method accurately reproduces the experimental results for Taylor impact tests over a wide range of impact velocities. The resulting model retains the advantages of the peridynamic formulation regarding discontinuities while allowing greater generality in material response than was previously possible. Copyright © 2009 John Wiley & Sons, Ltd.

In this paper, we develop a new non-ordinary state-based peridynamic method to solve transient dynamic solid mechanics problems. This new peridynamic method has advantages over the previously developed bond-based and ordinary state-based peridynamic methods in that its bonds are not restricted to central forces, nor is it restricted to a Poisson’s ratio of 1/4 as with the bond-based method. First, we obtain non-local nodal deformation gradients that are used to define nodal strain tensors. The deformation gradient tensors are used with the nodal strain tensors to obtain rate of deformation tensors in the deformed configuration. The polar decomposition of the deformation gradient tensors are then used to obtain the nodal rotation tensors which are used to rotate the rate of deformation tensors and previous Cauchy stress tensors into an unrotated configuration. These are then used with conventional Cauchy stress constitutive models in the unrotated state where the unrotated Cauchy stress rate is objective. We then obtain the unrotated Cauchy nodal stress tensors and rotate them back into the deformed configuration where they are used to define the forces in the nodal connecting bonds. As a first example we quasi-statically stretch a bar, hold it, and then rotate it ninety degrees to illustrate the methods finite rotation capabilities. Next, we verify our new method by comparing small strain results from a bar fixed at one end and subjected to an initial velocity gradient with results obtained from the corresponding one-dimensional small strain analytical solution. As a last example, we show the fracture capabilities of the method using both a notched and un-notched bar.

We introduce here adaptive refinement algorithms for the non-local method peridynamics, which was proposed in (J. Mech. Phys. Solids 2000; 48:175–209) as a reformulation of classical elasticity for discontinuities and long-range forces. We use scaling of the micromodulus and horizon and discuss the particular features of adaptivity in peridynamics for which multiscale modeling and grid refinement are closely connected. We discuss three types of numerical convergence for peridynamics and obtain uniform convergence to the classical solutions of static and dynamic elasticity problems in 1D in the limit of the horizon going to zero. Continuous micromoduli lead to optimal rates of convergence independent of the grid used, while discontinuous micromoduli produce optimal rates of convergence only for uniform grids. Examples for static and dynamic elasticity problems in 1D are shown. The relative error for the static and dynamic solutions obtained using adaptive refinement are significantly lower than those obtained using uniform refinement, for the same number of nodes. Copyright © 2008 John Wiley & Sons, Ltd.

The traditional methods for analyzing deformation in structures attempt to solve the partial differential equations of the classical theory of continuum mechanics. Yet these equations, because they require the partial derivatives of displacement to be known throughout the region modeled, are in some ways unsuitable for the modeling of discontinuities caused by damage, in which these derivatives fail to exist. As a means of avoiding this limitation, the peridynamic model of solid mechanics has been developed for applications involving discontinuities. The objective of this method is to treat crack and fracture as just another type of deformation, rather than as pathology that requires special mathematical treatment. The peridynamic theory is based on integral equations so there is no problem in applying the equations across discontinuities. The peridynamic method has been applied successfully to damage and failure analysis in composites. It predicts in detail the delamination and matrix damage process in c...