Search

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.

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

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.

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

The peridynamic model of solid mechanics is a nonlocal theory containing a length scale. It is based on direct interactions between points in a continuum separated from each other by a finite distance. The maximum interaction distance provides a length scale for the material model. This paper addresses the question of whether the peridynamic model for an elastic material reproduces the classical local model as this length scale goes to zero. We show that if the motion, constitutive model, and any nonhomogeneities are sufficiently smooth, then the peridynamic stress tensor converges in this limit to a Piola-Kirchhoff stress tensor that is a function only of the local deformation gradient tensor, as in the classical theory. This limiting Piola-Kirchhoff stress tensor field is differentiable, and its divergence represents the force density due to internal forces. The limiting, or collapsed, stress-strain model satisfies the conditions in the classical theory for angular momentum balance, isotropy, objectivity, and hyperelasticity, provided the original peridynamic constitutive model satisfies the appropriate conditions.

The peridynamic model is a framework for continuum mechanics based on the idea that pairs of particles exert forces on each other across a finite distance. The equation of motion in the peridynamic model is an integro-differential equation. In this paper, a notion of a peridynamic stress tensor derived from nonlocal interactions is defined. At any point in the body, this stress tensor is obtained from the forces within peridynamic bonds that geometrically go through the point. The peridynamic equation of motion can be expressed in terms of this stress tensor, and the result is formally identical to the Cauchy equation of motion in the classical model, even though the classical model is a local theory. We also establish that this stress tensor field is unique in a certain function space compatible with finite element approximations.

We developed an abrasion model that predicts mass loss and change in nose shape for steel projectiles that penetrate concrete targets. Mass loss data from four sets of experiments with two ogive-nose projectile geometries and concrete targets with limestone and quartz aggregates were used to develop the abrasion model. We plotted post-test mass loss versus initial kinetic energy and found a nearly linear dependence for striking velocities to approximately 1000 m/s. With this linear relationship, we derived a mathematical model that was implemented into the Sandia-developed, Eulerian hydrocode CTH. Predictions from CTH agreed well with experimental observations.

A generalization of the original peridynamic framework for solid mechanics is proposed. This generalization permits the response of a material at a point to depend collectively on the deformation of all bonds connected to the point. This extends the types of material response that can be reproduced by peridynamic theory to include an explicit dependence on such collectively determined quantities as volume change or shear angle. To accomplish this generalization, a mathematical object called a deformation state is defined, a function that maps any bond onto its image under the deformation. A similar object called a force state is defined, which contains the forces within bonds of all lengths and orientation. The relation between the deformation state and force state is the constitutive model for the material. In addition to providing a more general capability for reproducing material response, the new framework provides a means to incorporate a constitutive model from the conventional theory of solid mechanics directly into a peridynamic model. It also allows the condition of plastic incompressibility to be enforced in a peridynamic material model for permanent deformation analogous to conventional plasticity theory.

The peridynamic model described in Silling (Silling, S.A., 1998. Reformation of Elasticity Theory for Discontinuous and Long-Range Forces, SAND98-2176. Sandia National Laboratories, Albuquerque, NM), being a central-force model, is limited to modeling materials with a Poisson's ratio of 1/4. In this paper, the peridynamic model is generalized by adding pairwise peridynamic moments to simulate linear elastic materials with varying Poisson's ratios. The new model is called the “micropolar peridynamic model”. The micropolar peridynamic model is placed within a finite element context to enable efficacious application of boundary conditions and efficient computational solutions using an implicit, rather than an explicit solution algorithm. The implicit solution algorithm is suitable for quasistatic simulation of damage and cracking in concrete structures. With this new model, very simple tensile damage mechanisms at the micro structural (peridynamic) level are sufficient to explain a great deal of the microcracking (damage) and fracture mechanics observed in concrete structures. The new implementation appears to be computationally efficient.

An alternative theory of solid mechanics, known as the peridynamic theory, formulates problems in terms of integral equations rather than partial differential equations. This theory assumes that particles in a continuum interact with each other across a finite distance, as in molecular dynamics. Damage is incorporated in the theory at the level of these two-particle interactions, so localization and fracture occur as a natural outgrowth of the equation of motion and constitutive models. A numerical method for solving dynamic problems within the peridynamic theory is described. Accuracy and numerical stability are discussed. Examples illustrate the properties of the method for modeling brittle dynamic crack growth.

The peridynamic theory of continuum mechanics allows damage, fracture, and long-range forces to be treated as natural components of the deformation of a material. In this paper, the peridynamic approach is applied to small thickness two- and one-dimensional structures. For membranes, a constitutive model is described appropriate for rubbery sheets that can form cracks. This model is used to perform numerical simulations of the stretching and dynamic tearing of membranes. A similar approach is applied to one-dimensional string like structures that undergrow stretching, bending, and failure. Long-range forces similar to van der Waals interactions at the nanoscale influence the equilibrium configurations of these structures, how they deform, and possibly self-assembly.

Some materials may naturally form discontinuities such as cracks as a result of deformation. As an aid to the modeling of such materials, a new framework for the basic equations of continuum mechanics, called the "peridynamic" formulation, is proposed. The propagation of linear stress waves in the new theory is discussed, and wave dispersion relations are derived. Material stability and its connection with wave propagation is investigated. It is demonstrated by an example that the reformulated approach permits the solution of fracture problems using the same equations either on or off the crack surface or crack tip. This is an advantage for modeling problems in which the location of a crack is not known in advance.

The deformation of an infinite bar subjected to a self-equilibrated load distribution is investigated using the peridynamic formulation of elasticity theory. The peridynamic theory differs from the classical theory and other nonlocal theories in that it does not involve spatial derivatives of the displacement field. The bar problem is formulated as a linear Fredholm integral equation and solved using Fourier transform methods. The solution is shown to exhibit, in general, features that are not found in the classical result. Among these are decaying oscillations in the displacement field and progressively weakening discontinuities that propagate outside of the loading region. These features, when present, are guaranteed to decay provided that the wave speeds are real. This leads to a one-dimensional version of St. Venant's principle for peridynamic materials that ensures the increasing smoothness of the displacement field remotely from the loading region. The peridynamic result converges to the classical result in the limit of short-range forces. An example gives the solution to the concentrated load problem, and hence provides the Green's function for general loading problems.