Bond-based peridynamic modelling of singular and nonsingular crack-tip fields

A static meshfree implementation of the bond-based peridynamics formulation for linearly elastic solids is applied to the study of the transition from local to nonlocal behavior of the stress and displacement fields in the vicinity of a crack front and other sources of stress concentration. The long-range nature of the interactions between material points that is intrinsic to and can be modulated within peridynamics enables the smooth transition from the square-root singular stress fields predicted by the classical (local) linear theory of elasticity, to the nonsingular fields associated with nonlocal theories. The accuracy of the peridynamics scheme and the transition from local to nonlocal behavior, which are dictated by the lattice spacing and micromodulus function, are assessed by performing an analysis of the boundary layer that surrounds the front of a two dimensional crack subjected to mode-I loading and of a cracked plate subjected to far-field tension.     

On displacement methods in planar anisotropic elasticity

Two displacement formulation methods are presented for problems of planar anisotropic elasticity. The first displacement method is based on solving the two governing partial differential equations simultaneously/ This method is a recapitulation of the orignal work of Eshelby, Read and Shockley [7] on generalized plane deformations of anisotropic elastic materials in the context of planar anisotropic elasticity.The second displacement method is based on solving the two governing equations separately. This formulation introduces a displacement function, which satisfies a fourth-order partial differential equation that is identical in the form to the one given by Lekhnitskii [6] for monoclinic materials using a stress function. Moreover, this method parallels the traditional Airy stress function method and thus the Lekhnitskii method for pure plane problems. Both the new approach and the Airy stress function method start with the equilibrium equations and use the same extended version of Green's theorem (Chou and Pagano [13], p. 114; Gao [11]) to derive the expressions for stress or displacement components in terms of a potential (stress or displacement) function (see also Gao [10, 11]). It is therefore anticipated that the displacement function involved in this new method could also be evaluated from measured data, as was done by Lin and Rowlands [17] to determine the Airy stress function experimentally.The two different displacement methods lead to two general solutions for problems of planar anisotropic elasticity. Although the two solutions differ in expressions, both of the depend on the complex roots of the same characteristic equation. Furthermore, this characteristic equation is identical to that obtained by Lekhnitskii [6] using a stress formulation. It is therefore concluded that the two displacement methods and Lekhnitskii's stress method are all equivalent for problems of planar anisotropic elasticity (see Gao and Rowlands [8] for detailed discussions).     

Analytical solution to a mixed boundary value elastic problem of a roller-guided panel of laminated composite

This study presents an analytical solution to elastic field in a roller-guided panel of symmetric cross-ply laminated composite material. The mixed boundary value two-dimensional plane stress elasticity problem is formulated in terms of a single displacement potential function. This reduces the problem to the solution of a single fourth order partial differential equation of equilibrium as the other equilibrium equation is satisfied automatically. The solution is obtained in terms of an infinite Fourier series. To present some numerical results, a panel of glass/epoxy laminated composite is considered and different components of stress and displacement at different sections of the panel are presented graphically. To justify the present analytical solution, it is compared with the finite element solution obtained by using the commercial software ANSYS. It is found that the two solutions agree well with each other. This ensures that the formulation developed in this study based on the displacement potential approach can be used to obtain analytical solution of an elastic field in structural elements of laminated composite under any mode of boundary conditions prescribed in terms of either stress, displacement or any combination of these.     

A peridynamic implementation of crystal plasticity

This paper presents the first application of peridynamics theory for crystal plasticity simulations. A state-based theory of peridynamics is used (Silling et al., 2007) where the forces in the bonds between particles are computed from stress tensors obtained from crystal plasticity. The stress tensor at a particle, in turn, is computed from strains calculated by tracking the motion of surrounding particles. We have developed a quasi-static implementation of the peridynamics theory. The code employs an implicit iterative solution procedure similar to a non-linear finite element implementation. Peridynamics results are compared with crystal plasticity finite element (CPFE) analysis for the problem of plane strain compression of a planar polycrystal. The stress, strain field distribution and the texture formation predicted by CPFE and peridynamics were found to compare well. One particular feature of peridynamics is its ability to model fine shear bands that occur naturally in deforming polycrystalline aggregates. Peridynamics simulations are used to study the origin and evolution of these shear bands as a function of strain and slip geometry.     

The G2 constant displacement discontinuity method – Part I: Solution of plane crack problems

A new constant displacement discontinuity element was presented in a previous paper applied initially for the numerical solution of either isolated straight cracks or for co-linear cracks of the three fundamental deformation modes I, II and III due to the special form of the solution. It was based on the strain-gradient elasticity theory in its simplest possible Grade-2 variant. The assumption of the G2 expression for the stresses has resulted to a better average stress value at the mid-point of the straight displacement discontinuity compared to the classical elasticity solution. This new element gave considerably better predictions of the stress intensity factors compared to the constant displacement discontinuity element and the linear displacement discontinuity element. Moreover, it preserved the simplicity and hence the high speed of computations. In this Part I, the solution for this element is extended for the analysis of cracks of arbitrary shape in an infinite plane isotropic elastic body and it is validated against three known analytical solutions.     

Two-dimensional complete rational analysis of functionally graded beams within symplectic framework

Exact solutions for generally supported functionally graded plane beams are given within the framework of symplectic elasticity. The Young’s modulus is assumed to exponentially vary along the longitudinal direction while the Poisson’s ratio remains constant. The state equation with a shift-Hamiltonian operator matrix has been established in the previous work, which is limited to the Saint-Venant solution. Here, a complete rational analysis of the displacement and stress distributions in the beam is presented by exploring the eigensolutions that are usually covered up by the Saint-Venant principle. These solutions play a significant role in the local behavior of materials that is usually ignored in the conventional elasticity methods but possibly crucial to the material/structure failures. The analysis makes full use of the symplectic orthogonality of the eigensolutions. Two illustrative examples are presented to compare the displacement and stress results with those for homogenous materials, demonstrating the effects of material inhomogeneity.     

On the displacement potential solution of plane problems of structural mechanics with mixed boundary conditions

The present paper describes the advancement of displacement potential approach in relation to solution of plane problems of structural mechanics with mixed mode of boundary conditions. Both the conditions of the plane stress and the plane strain are considered for analyzing the displacement and stress fields of the structural problem. Using the finite difference technique based on the present displacement potential approach for the case of the plane stress and the plane strain conditions, firstly an elastic cantilever beam subjected to a pure shear at its tip is solved and these two solutions (plane stress and plane strain) are compared with Timoshenko and Goodier cantilever beam bending solutions (Theory of elasticity, 2nd edn. McGraw-Hill, New York, 1951); secondly the above-mentioned displacement potential approach for the case of the plane stress and the plane strain conditions are applied to solve a one-end fixed square plate subjected to a combined loading at its tip. Effects of plane stress and plane strain on the elastic field of the plate are discussed in a comparative fashion. Limitations of Timoshenko and Goodier cantilever beam bending solutions (Theory of elasticity, 2nd edn. McGraw-Hill, New York, 1951) over the displacement potential approach for the case of the plane stress and the plane strain conditions are not only discussed but also the superiority of the present displacement potential approach for the case of the plane stress and the plane strain conditions are reflected in the present research work.     

Axisymmetric contact problem of cubic quasicrystalline materials

The axisymmetric elasticity theory of cubic quasicrystal was developed in Ref. [1]. The axisymmetric elasticity problem of cubic quasicrystal is reduced to a single higher-order partial differential equation by introducing a displacement function, based on which, the exact analytic solutions for the elastic field of an axisymmetric contact problem of cubic quasicrystalline materials are obtained for universal contact stress or contact displacement. The result shows that if the contact stress has order −1/2 singularity on the edge of the contact domain, the contact displacement is a constant in the contact domain. Conversely, if the contact displacement is a constant, the contact stress must have order −1/2 singularity on the edge of the contact domain. Project supported by the National Natural Science Foundation of China (No. 19972011).     

Static and Dynamic Green’s Functions in Peridynamics

We derive the static and dynamic Green’s functions for one-, two- and three-dimensional infinite domains within the formalism of peridynamics, making use of Fourier transforms and Laplace transforms. Noting that the one-dimensional and three-dimensional cases have been previously studied by other researchers, in this paper, we develop a method to obtain convergent solutions from the divergent integrals, so that the Green’s functions can be uniformly expressed as conventional solutions plus Dirac functions, and convergent nonlocal integrals. Thus, the Green’s functions for the two-dimensional domain are newly obtained, and those for the one and three dimensions are expressed in forms different from the previous expressions in the literature. We also prove that the peridynamic Green’s functions always degenerate into the corresponding classical counterparts of linear elasticity as the nonlocal length tends to zero. The static solutions for a single point load and the dynamic solutions for a time-dependent point load are analyzed. It is analytically shown that for static loading, the nonlocal effect is limited to the neighborhood of the loading point, and the displacement field far away from the loading point approaches the classical solution. For dynamic loading, due to peridynamic nonlinear dispersion relations, the propagation of waves given by the peridynamic solutions is dispersive. The Green’s functions may be used to solve other more complicated problems, and applied to systems that have long-range interactions between material points.     

圆柱型各向异性弹性力学平面问题

本文对圆柱型各向异性弹性力学平面问题的基本方程进行了改写。在此基础上，导出了应力函数Ｇ和位移函数φ，它们满足相同的控制方程，比文〔１〕的应力函数Ｆ的控制方程要简单，便于求得特解，并有Ｆ＝ｒＧ的关系。还对若干经典问题进行了求解。     

Non-local theory solution for in-plane shear of through crack

A non-local theory of elasticity is applied to obtain the plane strain stress and displacement field for a through crack under in-plane shear by using Schmidt's method. Unlike the classical elasticity solution, a lattice parameter enters into the problem that make the stresses finite at crack tip. Both the angular variations of the circumferential stress and strain energy density function are examined to associate their stationary value with locations of possible fracture initiation. The former criterion predicted a crack initiation angle of 54° from the plane of shear for the non-local solution as compared with about 75° for the classical elasticity solution. The latter criterion based on energy density yields a crack initiation angle of 80° for a Poisson's ratio of 0.28. This is much closer to the value that is predicted by the classical crack tips solution of elasticity.     

Solution of the plane stress problems of strain-hardening materials described by power-law using the complex pseudo-stress function

In the present paper,the compatibility equation for the plane stress problems of power-law materials is transformed into a biharmonic equation by introducing the so-calledcomplex pseudo-stress function,which makes it possible to solve the elastic-plastic planestress problems of strain hardening materials described by power-law using the complexvariable function method like that in the linear elasticity theory.By using this generalmethod,the close-formed analytical solutions for the stress,strain and displacementcomponents of the plane stress problems’of power-law materials is deduced in the paper,which can also be used to solve the elasto-plastic plane stress problems of strain-hardeningmaterials other than that described by power-law.As an example,the problem of a power-law material infinite plate containing a circular hole under uniaxial tension is solved byusing this method,the results of which are compared with those of a known asymptoticanalytical solution obtained by the perturbation method.     

An Exact Derivation from Three-Dimensions of the Linear Equations for the Bending of an Elastically Monoclinic Plate and Associated Three-Dimensional Solutions

The exact linear three-dimensional equations for a elastically monoclinic (13 constant) plate of constant thickness are reduced without approximation to a single 4th order differential equation for a thickness-weighted normal displacement plus two auxiliary equations for weighted thickness integrals of a stress function and the normal strain. The 4th order equation is of the same form as in classical (Kirchhoff) theory except the unknown is not the midsurface normal displacement. Assuming a solution of these plate equations, we construct so-called modified Saint-Venant solutions—“modified” because they involve non-zero body and surface loads. That is, solutions of the exact three-dimensional elasticity equations that exhibit no boundary layers and that are subject to a special set of body and surface loads that leave the analogous plate loads arbitrary.     

Elasticity solutions for plane anisotropic functionally graded beams

This paper considers the plane stress problem of generally anisotropic beams with elastic compliance parameters being arbitrary functions of the thickness coordinate. Firstly, the partial differential equation, which is satisfied by the Airy stress function for the plane problem of anisotropic functionally graded materials and involves the effect of body force, is derived. Secondly, a unified method is developed to obtain the stress function. The analytical expressions of axial force, bending moment, shear force and displacements are then deduced through integration. Thirdly, the stress function is employed to solve problems of anisotropic functionally graded plane beams, with the integral constants completely determined from boundary conditions. A series of elasticity solutions are thus obtained, including the solution for beams under tension and pure bending, the solution for cantilever beams subjected to shear force applied at the free end, the solution for cantilever beams or simply supported beams subjected to uniform load, the solution for fixed–fixed beams subjected to uniform load, and the one for beams subjected to body force, etc. These solutions can be easily degenerated into the elasticity solutions for homogeneous beams. Some of them are absolutely new to literature, and some coincide with the available solutions. It is also found that there are certain errors in several available solutions. A numerical example is finally presented to show the effect of material inhomogeneity on the elastic field in a functionally graded anisotropic cantilever beam.     

Crack problem under shear loading in cubic quasicrystal

IntroductionQuasicrystalasanewstructureofsolidmatter[1,2 ]bringsprofoundnewideastothetraditionalcondensedmatterphysicsandencouragesconsiderabletheoreticalandexperimentalstudiesonthephysicalandmechanicalpropertiesofthematerial,includingtheelasticitytheoryofthequasicrystal,manyvaluableresultsweregiven[3～ 5 ].Defectsinthematerialwereobservedsoonafterthediscoveryofthequasicrystal[6 ,7].Cracksareonetypeofdefects,theirexistencegreatlyinfluencesthephysicalandmechanicalpropertiesofthequasicrystalinem…     

An efficient method for computing local quantities of interest in elasticity based on finite element error estimation

We present an efficient finite element method for computing the engineering quantities of interest that are linear functionals of displacement in elasticity based on a posteriori error estimate. The accuracy of quantities is greatly improved by adding the approximate cross inner product of errors in the primal and dual problems, which is calculated with an inexpensive gradient recovery type error estimate, to the quantities obtained from the finite element solution. With less CPU time, the accuracy of the improved quantities obtained with the proposed method on the coarse finite element mesh is similar to that of the quantities obtained from the finite element solutions on the finer mesh. Three quantities related to the local displacement, local stress and stress intensity factor are computed with the proposed method to verify its efficiency.     

The scattering of harmonic elastic anti-plane shear waves by two collinear cracks in anisotropic material plane by using the non-local theory

In this paper, the dynamic behavior of two collinear cracks in the anisotropic elasticity material plane subjected to the harmonic anti-plane shear waves is investigated by use of the nonlocal theory. To overcome the mathematical difficulties, a one-dimensional nonlocal kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress field near the crack tips. By use of the Fourier transform, the problem can be solved with the help of a pair of triple integral equations, in which the unknown variable is the displacement on the crack surfaces. To solve the triple integral equations, the displacement on the crack surfaces is expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present near crack tips. The nonlocal elasticity solutions yield a finite hoop stress at the crack tips, thus allowing us to using the maximum stress as a fracture criterion. The magnitude of the finite stress field not only depends on the crack length but also on the frequency of the incident waves and the lattice parameter of the materials.     

Three-dimensional elasticity solution for bending of transversely isotropic functionally graded plates

This paper presents a three-dimensional elasticity solution for a simply supported, transversely isotropic functionally graded plate subjected to transverse loading, with Young’s moduli and the shear modulus varying exponentially through the thickness and Poisson’s ratios being constant. The approach makes use of the recently developed displacement functions for inhomogeneous transversely isotropic media. Dependence of stress and displacement fields in the plate on the inhomogeneity ratio, geometry and degree of anisotropy is examined and discussed. The developed three-dimensional solution for transversely isotropic functionally graded plate is validated through comparison with the available three-dimensional solutions for isotropic functionally graded plates, as well as the classical and higher-order plate theories.     

A unified treatment of axisymmetric adhesive contact problems using the harmonic potential function method

A unified treatment of axisymmetric adhesive contact problems is provided using the harmonic potential function method for axisymmetric elasticity problems advanced by Green, Keer, Barber and others. The harmonic function adopted in the current analysis is the one that was introduced by Jin et al. (2008) to solve an external crack problem. It is demonstrated that the harmonic potential function method offers a simpler and more consistent way to treat non-adhesive and adhesive contact problems. By using this method and the principle of superposition, a general solution is derived for the adhesive contact problem involving an axisymmetric rigid punch of arbitrary shape and an adhesive interaction force distribution of any profile. This solution provides analytical expressions for all non-zero displacement and stress components on the contact surface, unlike existing ones. In addition, the newly derived solution is able to link existing solutions/models for axisymmetric non-adhesive and adhesive contact problems and to reveal the connections and differences among these solutions/models individually obtained using different methods at various times. Specifically, it is shown that Sneddon’s solution for the axisymmetric punch problem, Boussinesq’s solution for the flat-ended cylindrical punch problem, the Hertz solution for the spherical punch problem, the JKR model, the DMT model, the M-D model, and the M-D-n model can all be explicitly recovered by the current general solution.