scispace - formally typeset
Search or ask a question

Showing papers in "Journal of Elasticity in 2009"


ReportDOI
TL;DR: In this paper, the response of a state-based peridynamic material model for a small deformation superposed on a large deformation was investigated for the case of a crack.
Abstract: 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 Frechet 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 Poincare’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.

245 citations


Journal ArticleDOI
TL;DR: In this article, the Francfort-Marigo model was modified to incorporate the idea of less brittle, "deviatoric-type fracture" and apply to materials such as confined stone.
Abstract: Recently, Francfort and Marigo (J. Mech. Phys. Solids 46, 1319–1342, 1998) have proposed a novel approach to fracture mechanics based upon the global minimization of a Griffith-like functional, composed of a bulk and a surface energy term. Later on the same authors, together with Bourdin, introduced (in J. Mech. Phys. Solids 48, 797–826, 2000) a variational approximation (in the sense of Γ-convergence) of such functional, essentially for computational purposes. Here, we utilize this new variational approach to show how it might be altered to incorporate the idea of less brittle, “deviatoric-type fracture” and apply to materials such as confined stone. To do so, we modify the original formulation of Francfort and Marigo, in particular its approximation of Bourdin, Francfort and Marigo, to only allow for discontinuities in the deviatoric part of the strain. We apply such modified model to gain insight on the deterioration and cracking in the ashlar masonry work of the French Pantheon, which are so repetitious and particular to be a distinguishable symptom of ongoing damage. Numerical experiments are performed and the results compared to those obtained using the original Francfort-Marigo model and to actual crack patterns from the Pantheon. The modified formulation allows one to reproduce fracture paths surprisingly similar to that observed in situ, to sort out the possible causes of damage, and to confirm, with a quantitative analysis, the main structural deficiencies in the French monument. This practical example enhances the importance of this promising new theory based in the mathematical sciences.

134 citations


Journal ArticleDOI
TL;DR: In this article, a physically-based approach to non-local elasticity theory is introduced, which is formulated by reverting the continuum to an ensemble of interacting volume elements, where interactions between adjacent elements are classical contact forces while long-range interactions between non-adjacent elements are modelled as distance-decaying central body forces.
Abstract: In this paper the physically-based approach to non-local elasticity theory is introduced. It is formulated by reverting the continuum to an ensemble of interacting volume elements. Interactions between adjacent elements are classical contact forces while long-range interactions between non-adjacent elements are modelled as distance-decaying central body forces. The latter are proportional to the relative displacements rather than to the strain field as in the Eringen model and subsequent developments. At the limit the displacement field is found to be governed by an integro-differential equation, solved by a simple discretization procedure suggested by the underlying mechanical model itself, with corresponding static boundary conditions enforced in a quite simple form. It is then shown that the constitutive law of the proposed model coalesces with the Eringen constitutive law for an unbounded domain under suitable assumptions, whereas it remains substantially different for a bounded domain. Thermodynamic consistency of the model also has been investigated in detail and some numerical applications are presented for different parameters and different functional forms for the decay of the long range forces. For simplicity, the problem is formulated for a 1D continuum while the general formulation for a 3D elastic solid has been reported in the appendix.

131 citations


Journal ArticleDOI
TL;DR: In this paper, a methodology for investigating the dynamics of heterogeneous media using the nonlocal continuum model given by the peridynamic formulation is presented. But the approach presented here provides the ability to model the macroscopic dynamics while at the same time resolving the dynamics at the length scales of the microstructure.
Abstract: A methodology is presented for investigating the dynamics of heterogeneous media using the nonlocal continuum model given by the peridynamic formulation. The approach presented here provides the ability to model the macroscopic dynamics while at the same time resolving the dynamics at the length scales of the microstructure. Central to the methodology is a novel two-scale evolution equation. The rescaled solution of this equation is shown to provide a strong approximation to the actual deformation inside the peridynamic material. The two scale evolution can be split into a microscopic component tracking the dynamics at the length scale of the heterogeneities and a macroscopic component tracking the volume averaged (homogenized) dynamics. The interplay between the microscopic and macroscopic dynamics is given by a coupled system of evolution equations. The equations show that the forces generated by the homogenized deformation inside the medium are related to the homogenized deformation through a history dependent constitutive relation.

94 citations


Journal ArticleDOI
TL;DR: In this article, the authors derived necessary and sufficient conditions for the strong ellipticity condition of anisotropic elastic materials, which holds if and only if a second order tensor function is positive definite for any unit vectors.
Abstract: In this paper, we derive necessary and sufficient conditions for the strong ellipticity condition of anisotropic elastic materials. We first observe that the strong ellipticity condition holds if and only if a second order tensor function is positive definite for any unit vectors. Then we further link this condition to the rank-one positive definiteness of three second-order tensors, three fourth-order tensors and a sixth-order tensor. In particular, we consider conditions of strong ellipticity of the rhombic classes, for which we need to check the copositivity of three second-order tensors and the positive definiteness of a sixth-order tensor. A direct method is presented to verify our conditions.

85 citations


Journal ArticleDOI
TL;DR: In this article, an axisymmetrical FGM hollow cylinder or disk with arbitrarily varying material properties is considered, and a simple and efficient approach is suggested, which reduces the associated problem to solving a Fredholm integral equation.
Abstract: The elastic analysis of a pressurized functionally graded material (FGM) annulus or tube is made in this paper. Different from existing studies, this study deals with an axisymmetrical FGM hollow cylinder or disk with arbitrarily varying material properties. A simple and efficient approach is suggested, which reduces the associated problem to solving a Fredholm integral equation. The resulting equation is approximately solved by expanding the solution as series of Legendre polynomials. The stresses and displacements can be represented in terms of the solution to the equation. For radius-dependent Young’s modulus, numerical results of the distribution of the radial and circumferential stresses are presented graphically. Our results indicate that change in the gradient of the FGM tube does not produce a substantial variation of the radial stress, but strongly affects the distribution of the hoop stress. In particular, the hoop stress may reach its maximum at an internal position or at the outer surface when the tube is internally pressurized. The results obtained are helpful in designing FGM cylindrical vessels to prevent failure.

74 citations


Journal ArticleDOI
TL;DR: In this article, a closed-form criterion for the onset of cavitation in compressible, isotropic, hyperelastic solids subjected to non-symmetric loading conditions is derived.
Abstract: In this work, we derive a closed-form criterion for the onset of cavitation in compressible, isotropic, hyperelastic solids subjected to non-symmetric loading conditions. The criterion is based on the solution of a boundary value problem where a hyperelastic solid, which is infinite in extent and contains a single vacuous inhomogeneity, is subjected to uniform displacement boundary conditions. By making use of the “linear-comparison” variational procedure of Lopez-Pamies and Ponte Castaneda (J. Mech. Phys. Solids 54:807–830, 2006), we solve this problem approximately and generate variational estimates for the critical stretches applied on the boundary at which the cavity suddenly starts growing. The accuracy of the proposed analytical result is assessed by comparisons with exact solutions available from the literature for radially symmetric cavitation, as well as with finite element simulations. In addition, applications are presented for a variety of materials of practical and theoretical interest, including the harmonic, Blatz-Ko, and compressible Neo-Hookean materials.

49 citations


Journal ArticleDOI
TL;DR: In this article, the authors investigated conservation and balance laws in the framework of linear elastodynamics considering the strain energy density depending on the gradients of the displacement up to the third order.
Abstract: In the present paper we investigate conservation and balance laws in the framework of linear elastodynamics considering the strain energy density depending on the gradients of the displacement up to the third order, as originally proposed by Mindlin (Int. J. Solids Struct. 1, 417–438, 1965). The conservation and balance laws that correspond to the symmetries of translation, rotation, scaling and addition of solutions are derived using Noether’s theorem. Also, the formulas of the dynamical J,L and M-integrals are presented for the problem under study. Moreover, the balance law of addition of solutions gives rise to explore the dynamical reciprocal theorem as well as the restrictions under which it is valid.

37 citations


Journal ArticleDOI
TL;DR: In this paper, an extensive global treatment of radial motions of compressible nonlinearly viscoelastic cylindrical and spherical shells under time-dependent pressures is presented, where the shells are described by a geometrically exact 2-dimensional theory in which the shells suffer thickness strains as well as the standard stretching of their base surfaces, and their motions are governed by fourth-order systems of semilinear ordinary differential equations.
Abstract: This paper contains an extensive global treatment of radial motions of compressible nonlinearly viscoelastic cylindrical and spherical shells under time-dependent pressures. It furnishes a variety of conditions on a general class of material properties and on the pressure terms ensuring that there are solutions existing for all times, there are unbounded globally defined solutions, there are solutions that blow up in finite time, and there are solutions having the same period as that of the pressure terms. The shells are described by a geometrically exact 2-dimensional theory in which the shells suffer thickness strains as well as the standard stretching of their base surfaces. Consequently their motions are governed by fourth-order systems of semilinear ordinary differential equations. This work shows that there are major qualitative differences between the nonlinear dynamical behaviors of cylindrical and spherical shells.

34 citations


Journal ArticleDOI
TL;DR: This work obtains membrane theory directly from the local equations and boundary conditions of the three-dimensional theory of elastic surfaces.
Abstract: Membrane theory may be regarded as a special case of the Cosserat theory of elastic surfaces, or, alternatively, derived from three-dimensional elasticity theory via asymptotic or variational methods. Here we obtain membrane theory directly from the local equations and boundary conditions of the three-dimensional theory.

33 citations


Journal ArticleDOI
Duvan Henao1
TL;DR: In this paper, it was shown that energy minimizers for nonlinear elasticity in which cavitation is allowed only at a finite number of prescribed flaw points can be obtained, in the limit as e→0, by introducing micro-voids of radius e in the domain at the prescribed locations and minimizing the energy without allowing for cavitation.
Abstract: We prove that energy minimizers for nonlinear elasticity in which cavitation is allowed only at a finite number of prescribed flaw points can be obtained, in the limit as e→0, by introducing micro-voids of radius e in the domain at the prescribed locations and minimizing the energy without allowing for cavitation. This extends the result by Sivaloganathan, Spector, and Tilakraj (SIAM J. Appl. Math. 66:736–757, 2006) to the case of multiple cavities, and constitutes a first step towards the numerical simulation of cavitation (in the nonradially-symmetric case).

Journal ArticleDOI
TL;DR: In this paper, the problem of finding the transversely isotropic tensor closest to a given elasticity tensor with respect to the Frobenius norm was considered and a method for finding the optimal orientation of the coordinate system was formulated.
Abstract: We consider the problem of finding the transversely isotropic elasticity tensor closest to a given elasticity tensor with respect to the Frobenius norm. A similar problem was considered by other authors and solved analytically assuming a fixed orientation of the natural coordinate system of the transversely isotropic tensor. In this paper we formulate a method for finding the optimal orientation of the coordinate system—the one that produces the shortest distance. The optimization problem reduces to finding the absolute maximum of a homogeneous eighth-degree polynomial on a two-dimensional sphere. This formulation allows us a convenient visualization of local extrema, and enables us to find the closest transversely isotropic tensor numerically.

Journal ArticleDOI
TL;DR: In this article, the macroscopic response of porous Neo-Hookean solids with random microstructures was derived for an exact and closed-form solution to a Hamilton-Jacobi equation.
Abstract: Making use of the particulate microgeometries of Idiart (J. Mech. Phys. Solids 56:2599–2617, 2008), we derive an exact and closed-form result for the macroscopic response of porous Neo-Hookean solids with random microstructures. The stored-energy function is a solution to a Hamilton-Jacobi equation with initial porosity and macroscopic deformation gradient playing the roles of time and space. The main theoretical and practical aspects of the result are discussed.

Journal ArticleDOI
TL;DR: In this paper, the eigenvalues of both the given tensor and the projected tensor are used to estimate the distance of a given elasticity tensor to other tensors exhibiting particular material symmetries.
Abstract: Using orthogonal projections, we investigate distance of a given elasticity tensor to classes of elasticity tensors exhibiting particular material symmetries. These projections depend on the orientation of the elasticity tensor; hence the distance is obtained as the minimization of corresponding expressions with respect to the action of the orthogonal group. These expressions are stated in terms of the eigenvalues of both the given tensor and the projected one. The process of minimization is facilitated by the fact that, as we prove, the traces of the corresponding Voigt and dilatation tensors are invariant under these orthogonal projections. For isotropy, cubic symmetry and transverse isotropy, we formulate algorithms to find both the orientation and the eigenvalues of the elasticity tensor endowed with a particular symmetry and closest to the given elasticity tensor.

Journal ArticleDOI
TL;DR: In this paper, the authors consider a body whose stored energy is quadratic with respect to the Green-St. Venant strain and show that suitable minimizing sequences converge to the unique minimizer of the energy functional of linear elasticity.
Abstract: We consider a residually-stressed, uniform hyperelastic body whose stored energy is quadratic with respect to the Green–St. Venant strain. We show that, in the limit of vanishing loads, suitable minimizing sequences converge to the unique minimizer of the energy functional of linear elasticity. We also deduce the standard stress-strain relations for linear elasticity with residual stress.

Journal ArticleDOI
TL;DR: In this article, a method for the solution of the axisymmetric boundary value problem for a finite elastic cylinder with assigned stress and/or displacements acting on the ends and side was proposed.
Abstract: In this paper, we propose a method for the solution of the axisymmetric boundary value problem for a finite elastic cylinder with assigned stress and/or displacements acting on the ends and side. The technique utilizes the Love representation, which allows for reduction of the solution of the elastic problem to the search for a biharmonic function on a cylindrical domain. In the solution method suggested here, we write the Love function with a Bessel expansion and analyze in detail the conditions under which it is possible to differentiate the expansion term by term. We show that this is possible only for a restricted class of elastic solutions. In the general case, we introduce two new auxiliary functions of the z-coordinate. In this way, we obtain the general form of the axisymmetric biharmonic function, which is discussed in relation to certain specific boundary conditions applied on the side and ends of the cylinder. We obtain an exact explicit solution of practical interest for a cylinder with free ends and assigned displacements applied to the side.

Journal ArticleDOI
TL;DR: In this paper, Liu et al. showed that the classical entropy flux relation need not be valid in general, by considering a transversely isotropic elastic material body, and they proved that the results are identical to the classical results based on the Clausius-Duhem inequality (Coleman and Noll in Arch. Mech. Anal. 90:259-270, 2008).
Abstract: Recently (Liu in J. Elast. 90:259-270, 2008) thermodynamic theory of elastic (and viscoelastic) material bodies has been analyzed based on the general entropy inequality. It is proved that for isotropic elastic materials, the results are identical to the classical results based on the Clausius-Duhem inequality (Coleman and Noll in Arch. Ration. Mech. Anal. 13:167-178, 1963), for which one of the basic assumptions is that the entropy flux is defined as the heat flux divided by the absolute temperature. For anisotropic elastic materials in general, this classical entropy flux relation has not been proved in the new thermodynamic theory. In this note, as a supplement of the theory presented in (Liu in J. Elast. 90:259-270, 2008), it will be proved that the classical entropy flux relation need not be valid in general, by considering a transversely isotropic elastic material body.

Journal ArticleDOI
TL;DR: In this article, a state-space approach for exact analysis of axisymmetric deformation and stress distribution in a circular cylindrical body of transversely isotropic material is developed.
Abstract: A state-space approach for exact analysis of axisymmetric deformation and stress distribution in a circular cylindrical body of transversely isotropic material is developed. By means of Hamiltonian variational formulation via Legendre’s transformation, the basic equations in cylindrical coordinates are formulated into a state-space framework in which the unknown state vector comprises the displacements and associated stress components as the dual variables and the system matrix possesses the symplectic characteristics of a Hamiltonian system. Upon delineating the symplecticity of the formulation, a viable solution approach using eigenfunction expansion is developed. For illustration, an exact analysis of a finite thick-walled circular cylinder under internal and external pressures is presented, with emphasis on the end effects.

Journal ArticleDOI
TL;DR: In this paper, a large class of steady-state or transient nonlinear boundary conditions which provide asymptotic kinematics of Kirchhoff-love type were determined in the framework of small strains.
Abstract: Most of the derivations of the mechanical behavior of a plate as the limit behavior of a three-dimensional solid whose thickness tends to zero deal with stationary homogeneous linear boundary conditions on the lateral boundary. Here, in the framework of small strains, we rigorously determine a large class of steady-state or transient nonlinear boundary conditions which provide asymptotic kinematics of Kirchhoff-Love type.

Journal ArticleDOI
TL;DR: In this article, the geometry of the weave pattern affects the symmetry properties of the elastic energy of the surface of a woven fabric, and the relation of these invariants with the curvature and twist of the fibers is discussed.
Abstract: A woven fabric can be described as a surface made of two families of fibers: in this work we study how the geometry of the weave pattern affects the symmetry properties of the elastic energy of the surface Four basic symmetry classes of weave patterns are possible, depending on the angle between the fibers and their material properties The properties of the pattern determine the material symmetry group of the network, under which the elastic energy is invariant We derive representations for the energy of a woven fabric that are invariant under the symmetry group of the network, and discuss the relation of these invariants with the curvature and twist of the fibers

Journal ArticleDOI
TL;DR: In this paper, the authors use the theory of Lie groups/algebras to discuss the symmetries of crystals with uniform distributions of defects, and use the Lie group theory to describe the distribution of defects.
Abstract: I use the theory of Lie groups/algebras to discuss the symmetries of crystals with uniform distributions of defects

Journal ArticleDOI
TL;DR: In this paper, the authors studied the deformation of an anisotropic elliptic inclusion embedded in an infinite dissimilar matrix subject to a uniform loading at infinity, where the surface traction is continuous across the interface while the displacement is discontinuous.
Abstract: In this paper we study the two-dimensional deformation of an anisotropic elliptic inclusion embedded in an infinite dissimilar anisotropic matrix subject to a uniform loading at infinity. The interface is assumed to be imperfectly bonded. The surface traction is continuous across the interface while the displacement is discontinuous. The interface function that relates the surface traction and the displacement discontinuity across the interface is a tensor function, not a scalar function as employed by most work in the literature. We choose the interface function such that the stress inside the elliptic inclusion is uniform. Explicit solution for the inclusion and the matrix is presented. The materials in the inclusion and in the matrix are general anisotropic elastic materials so that the antiplane and inplane displacements are coupled regardless of the applied loading at infinity.

Journal ArticleDOI
TL;DR: In this article, the strong Eshelby conjecture holds for simply-connected anisotropic inclusion problems under plane or anti-plane deformation, where the interfaces can be either perfect or dislocation-like.
Abstract: Based on the Stroh formalism for anisotropic elasticity and the complex variable function method, we prove in this paper that the strong Eshelby conjecture holds for simply-connected anisotropic inclusion problems under plane or anti-plane deformation. The interfaces can be either perfect or dislocation-like. For these inclusion problems, if the induced stress field inside the inclusion is uniform for a single uniform eigenstrain, the inclusion is of the elliptic shape. Thanks to the generality of the proof method, we obtain also alternative proofs of the strong Eshelby conjecture for isotropic inclusion problems, which are given in the Appendix.

Journal ArticleDOI
TL;DR: In this article, the authors proposed a solution of the Saint-Venant's Problem for a slender compound piezoelectric beam presented in this paper generalizes the recent solution by the authors and E. Harash (J. Appl. Mech. 11:1-10, 2007) for a homogeneous beam and the solution for a compound elastic beam developed by O. Rand.
Abstract: The solution of the Saint-Venant’s Problem for a slender compound piezoelectric beam presented in this paper generalizes the recent solution by the authors and E. Harash (J. Appl. Mech. 11:1–10, 2007) for a homogeneous piezoelectric beam and the solution for a compound elastic beam developed by O. Rand and the first author (Analytical Methods in Anisotropic Elasticity with Symbolic Computational Tools, Birkhauser, Boston, 2005). Justification for this approximation emerges from the St. Venant’s Principle. The stress, strain and (electrical) displacement components (“solution hypothesis”) are presented as a set of initially assumed expressions involving twelve tip loading parameters, six unknown weight coefficients, and three pairs of torsion/bending functions of two variables. Each pair of functions satisfies the so-called coupled non-homogeneous Neumann problem (CNNP) in the cross-sectional domain. The work develops concepts of the torsion/bending functions, the torsional rigidity and piezoelectric shear center, the tip coupling matrix, for a compound piezoelectric beam. Examples of exact and approximate solutions for rectangular laminated beams made of transtropic materials are presented.

Journal ArticleDOI
TL;DR: In this paper, a one-dimensional variational problem for an anisotropic, partially inhomogeneous, residually stressed, rectangular thin-walled beam is derived, by Γ-convergence, from the three dimensional theory of linear elasticity with residual stress.
Abstract: A one-dimensional variational problem for an anisotropic, partially inhomogeneous, residually stressed, rectangular thin-walled beam is derived, by Γ-convergence, from the three-dimensional theory of linear elasticity with residual stress

Journal ArticleDOI
TL;DR: A theoretical framework for the description of the interaction between diffusion, mechanics, and degradation in elastic solids is developed in this article to avoid complications that obscure the essential features of these interactions, and a particular specialization of the general theory is selected and a numerical implementation based on the finite element method, a backward Euler time-stepping scheming, and an operator-splitting algorithm is described.
Abstract: A theoretical framework for the description of the interaction between diffusion, mechanics, and degradation in elastic solids is developed To avoid complications that obscure the essential features of these interactions, we work within a one-dimensional setting A particular specialization of the general theory is selected and a numerical implementation based on the finite-element method, a backward Euler time-stepping scheming, and an operator-splitting algorithm is described An application involving the time-independent end-loading of a notched cylindrical bar is used to illustrate the ability of the theory to describe some essential features of solute-assisted degradation

Journal ArticleDOI
TL;DR: In this paper, it was shown that for homogeneous bodies, the linear elastic solution is by no means a rigid rotation (in a particular case it is an “explosion”) while the residual stress may even become infinite.
Abstract: Material frame indifference implies that the solution in non-linear elasticity theory for a connected body rigidly rotated at its border is a rigid, stress-free, deformation. If the same problem is considered within linear elasticity theory, considered as an approximation to the true elastic situation, one should expect that if the angle of rotation is small, the body still undergoes a rigid deformation while the corresponding stress, though not zero, remains consistently small. Here, we show that this is true, in general, only for homogeneous bodies. Counterexamples of inhomogeneous bodies are presented for which, whatever small the angle of rotation is, the linear elastic solution is by no means a rigid rotation (in a particular case it is an “explosion”) while the stress may even become infinite. If the same examples are re-interpreted as problems in an elasticity theory based upon genuinely linear constitutive relations which retain their validity also for finite deformations, it is shown that they would deliver constraint reaction forces that are not in equilibrium in the actual, deformed, state. This furnishes another characterization of the impossibility of an exact linear constitutive theory for elastic solids with zero residual stress.

Journal ArticleDOI
TL;DR: In this article, a static, purely flexural mechanical analysis is presented for a Kirchhoff solid circular plate, deflected by a transverse central force, and clamped along two antipodal arcs, the remaining part of the boundary being free.
Abstract: A static, purely flexural mechanical analysis is presented for a Kirchhoff solid circular plate, deflected by a transverse central force, and clamped along two antipodal arcs, the remaining part of the boundary being free. By adopting an integral formulation, the contact reaction is assumed to be formed by four equal concentrated forces acting at the support extremities, accompanied by two distributed moments with radial and circumferential axis, respectively. This plate problem is rephrased in terms of a complex-valued Hilbert singular integral equation of the second kind, whose solution is obtained in analytical, integral form. A design chart is presented that reports the plate central deflection as a function of the angular width of the plate supports.

Journal ArticleDOI
TL;DR: In this paper, the authors proposed to remove the singular stresses at the vertices of reentrant corners with traditional, stress-free, boundary conditions by recognizing cohesive stress action at sharp corners, and also introducing consistent cohesive laws ahead of corners.
Abstract: Elastic stress singularities occur at the vertices of reentrant corners with traditional, stress-free, boundary conditions. By recognizing cohesive stress action at sharp corners, and also introducing consistent cohesive laws ahead of corners, these singular stresses can be removed. This was earlier shown asymptotically for plane strain states: Here it is shown asymptotically for antiplane shear states.

Journal ArticleDOI
TL;DR: In this article, the authors show that the approach by Lur'e may have no solution and derive new Papkovich-Neuber harmonic functions for the ellipsoidal cavity problem from the solutions obtained by the Eshelby equivalent inclusion method.
Abstract: Lur’e (Three-dimensional Problem of the Theory of Elasticity. Interscience, New York, 1964, §6.9) presented an approach to solve the problem of an ellipsoidal cavity in a linear, elastic and isotropic medium loaded by uniform principal stresses at infinity. In this paper we show that the approach by Lur’e may have no solution. Derivation mistakes are first pointed out in his (6.9.22), (6.9.23), (6.9.30) and (6.9.31). With the correct expressions, we then prove the coefficient matrix in his (6.9.32) to be singular. Therefore constants A,A4,A5 may have no solution. The problem lies in the harmonic functions chosen by Lur’e for the Papkovich-Neuber solution. From the solutions obtained by the Eshelby equivalent inclusion method, the present paper derives new Papkovich-Neuber harmonic functions for the ellipsoidal cavity problem.