Peridynamic modeling of membranes and fibers

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.     

Linearized Theory of Peridynamic States

A state-based peridynamic material model describes internal forces acting on a point in terms of the collective deformation of all the material within a neighborhood of the point. In this paper, the response of a state-based peridynamic material is investigated for a small deformation superposed on a large deformation. The appropriate notion of a small deformation restricts the relative displacement between points, but it does not involve the deformation gradient (which would be undefined on a crack). The material properties that govern the linearized material response are expressed in terms of a new quantity called the modulus state. This determines the force in each bond resulting from an incremental deformation of itself or of other bonds. Conditions are derived for a linearized material model to be elastic, objective, and to satisfy balance of angular momentum. If the material is elastic, then the modulus state is obtainable from the second Fréchet derivative of the strain energy density function. The equation of equilibrium with a linearized material model is a linear Fredholm integral equation of the second kind. An analogue of Poincaré’s theorem is proved that applies to the infinite dimensional space of all peridynamic vector states, providing a condition similar to irrotationality in vector calculus.     

Peridynamic modeling of pitting corrosion damage

In this paper we introduce a peridynamic model for the evolution of damage from pitting corrosion capable of capturing subsurface damage. We model the anodic reaction in corrosion processes (in which electroplating is negligible) as an effective peridynamic diffusion process in the electrolyte/solid system coupled with a phase-change mechanism that allows for autonomous evolution of the moving interface. In order to simulate creation of subsurface damage, we introduce a corrosion damage model based on a stochastic relationship that connects the concentration in the metal to the damage of peridynamic mechanical-bonds that are superposed onto diffusion-bonds. We study convergence of this formulation for diffusion-dominated stage. The model leads to formation of a subsurface damage layer, seen in experiments. We validate results against experiments on pit growth rate and polarization data for pitting corrosion. We extend the 1D model to the 2D and 3D, and introduce a new damage-dependent corrosion model to account for broken mechanical bonds that enhance the corrosion rate. This coupled model can predict the pit shape and damage profile in materials with microstructural heterogeneities, such as defects, interfaces, inclusions, and grain boundaries.     

Peridynamic States and Constitutive Modeling

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.      

Nonlocal Constrained Value Problems for a Linear Peridynamic Navier Equation

In this paper, we carry out further mathematical studies of nonlocal constrained value problems for a peridynamic Navier equation derived from linear state-based peridynamic models. Given the nonlocal interactions effected in the model, constraints on the solution over a volume of nonzero measure are natural conditions to impose. We generalize previous well-posedness results that were formulated for very special kernels of nonlocal interactions. We also give a more rigorous treatment to the convergence of solutions to nonlocal peridynamic models to the solution of the conventional Navier equation of linear elasticity as the horizon parameter goes to zero. The results are valid for arbitrary Poisson ratio, which is a characteristic of the state-based peridynamic model.     

Peridynamic beams: A non-ordinary,state-based model

This paper develops a new peridynamic state based model to represent the bending of an Euler–Bernoulli beam. This model is non-ordinary and derived from the concept of a rotational spring between bonds. While multiple peridynamic material models capture the behavior of solid materials, this is the first 1D state based peridynamic model to resist bending. For sufficiently homogeneous and differentiable displacements, the model is shown to be equivalent to Eringen’s nonlocal elasticity. As the peridynamic horizon approaches 0, it reduces to the classical Euler–Bernoulli beam equations. Simple test cases demonstrate the model’s performance.     

A Peridynamic Model of Fracture Mechanics with Bond-Breaking

Wave propagation in elastic dielectrics with flexoelectricity, micro-inertia and strain gradient elasticity is investigated in this paper. Dispersion phenomenon, which does not exist in classical elastic dielectric theory, is observed in the flexoelectric microstructured solids. Analytical solutions for the phase velocity \(C_{p}\), group velocity \(C_{g}\) and their ratio \(\gamma = C_{g} / C_{p}\) are calculated for the case of harmonic decomposition. The magnitudes of the phase velocity and group velocity changed with the increasing of the wave number, while they are constant in the classical elastic dielectric theory. It is shown that the flexoelectricity, micro-inertia and microstructural effects are significant to predict the real behavior of longitudinal wave propagating in flexoelectric microstructured solids. Microstructural effects are not sufficient for dealing with realistic dispersion curves in flexoelectric solids, the micro-inertia and flexoelectricity are needed to obtain a physically acceptable value of the phase and group velocities.     

Analysis of the Volume-Constrained Peridynamic Navier Equation of Linear Elasticity

Well-posedness results for the state-based peridynamic nonlocal continuum model of solid mechanics are established with the help of a nonlocal vector calculus. The peridynamic strain energy density for an elastic constitutively linear anisotropic heterogeneous solid is expressed in terms of the field operators of that calculus, after which a variational principle for the equilibrium state is defined. The peridynamic Navier equilibrium equation is then derived as the first-order necessary conditions and are shown to reduce, for the case of homogeneous materials, to the classical Navier equation as the extent of nonlocal interactions vanishes. Then, for certain peridynamic constitutive relations, the peridynamic energy space is shown to be equivalent to the space of square-integrable functions; this result leads to well-posedness results for volume-constrained problems of both the Dirichlet and Neumann types. Using standard results, well-posedness is also established for the time-dependent peridynamic equation of motion.     

复合材料层合板冲击损伤近场动力学模型与分析

近场动力学(简称PD)理论通过域内积分建立物质基本运动方程。不同于传统理论中通过微分建立运动方程的方法,该理论对场函数没有连续性的要求,因而适合求解各类不连续问题。基于此,本文建立了正交各向异性单层板PD理论模型,进而引入单层板层间作用,发展了正交各向异性层合板PD模型及其损伤模型,并模拟了各向同性与各向异性层合板冲击损伤;通过对比分析,对模型的有效性进行了验证。     

Peridynamic modelling of impact damage in three-point bending beam with offset notch

The nonlocal peridynamic theory has been proven to be a promising method for the material failure and damage analyses in solid mechanics.Based upon the integrodifferential equations,peridynamics enables predicting the complex fracture phenomena such as spontaneous crack nucleation and crack branching,curving,and arrest.In this paper,the bond-based peridynamic approach is used to study the impact damage in a beam with an offset notch,which is widely used to investigate the mixed I-II crack propagation in brittle materials.The predictions from the peridynamic analysis agree well with available experimental observations.The numerical results show that the dynamic fracture behaviors of the beam under the impact load,such as crack initiation,curving,and branching,rely on the location of the offset notch and the impact speed of the drop hammer.     

Peridynamic modeling and simulation of coupled thermomechanical removal of ice from frozen structures

Meccanica - In this work, a bond-based peridynamic de-icing model has been developed to simulate the thermo-mechanical ice removal process of frozen structures. In the proposed numerical method,...     

Saint-Venant Torsion of a Two-Phase Circumferentially Symmetric Compound Bar

This work is concerned with the Saint-Venant torsion problem of a two-phase circumferentially symmetric compound prismatic bar. By generalizing a method originally proposed by Packham and Shail [12], we demonstrate that for a particular two-phase configuration, simply or multiply connected, which is invariant with phase interchange, the solutions can be constructed from solutions of two analogous problems with constant material properties. An effective shear modulus is derived in analytic form, which is approximately the harmonic mean of the component shear moduli. We also show that the effective torsional shear modulus is homogeneous for arbitrary configurations. This revised version was published online in August 2006 with corrections to the Cover Date.     

Deformation field of a misfitting inclusion

J.D. Eshelby (1957, 1959) has calculated the deformation field associated with an ellipsoidal inclusion in a state of homogeneous strain within an infinite matrix. Since most real precipitates occur with facets, the strain within such an inclusion is not uniform. Thus, plate precipitates of θ′ in Al-Cu and η in Al-Au have coherent broad faces with mismatches of 1.34 and 4.95 % respect- ively and semicoherent or disordered interfaces at the edges with residual mismatches of about ?4.3 and ?1.00% normal to the broad faces. The deformation field in the matrix around such precipitates has been calculated using Kelvin's (1848) result for the stress field due to a point force. The calculations show the existence of high stresses near the edges of the precipitates where they have an appreciable misfit. Unlike the case of an ellipsoidal inclusion, the stress fields of these precipitates have dilatational components which can affect the diffusion of solute atoms to them and, thus, the kinetics of interface migration. The behavior of alloys containing these precipitates indicates that the moduli of the precipitates are somewhat greater than those of the matrices. The present calculations, based on the assumption that the two moduli are the same, underestimate the actual deformation field in the matrix. In real systems, therefore, the effects of the deformation field on misfit dislocation nucleation and kinetics of interface migration are likely to be somewhat greater in general.     

Thermoelastoplastic Deformation of a Multilayer Ball

The problem of centrally symmetric deformation of a multilayer elastoplastic ball in the process of successive accretion of preheated layers to its outer surface is considered in the framework of small elastoplastic deformations. The problems of residual stress formation in the elastoplastic ball with an inclusion and a cavity are solved under various mechanical boundary conditions on the inner surface and for prescribed thermal compression distributions. The graphs of residual stress and displacement fields are constructed.     

多阶变截面杆的弹性稳定计算

本文研究单阶和多阶截面杆在集中和分布轴向荷载作用下的弹性稳定计算方法,首次给出了其整体稳定和局部稳定的解析解.