The energy-release rate and “self-force” of dynamically expanding spherical and plane inclusion boundaries with dilatational eigenstrain

In the context of the linear theory of elasticity with eigenstrains, the radiated field including inertia effects of a spherical inclusion with dilatational eigenstrain radially expanding is obtained on the basis of the dynamic Green's function, and one of the half-space inclusion boundary (with dilatational eigenstrain) moving from rest in general subsonic motion is obtained by a limiting process from the spherically expanding inclusion as the radius tends to infinity while the eigenstrain remains constrained, and this is the minimum energy solution. The global energy-release rate required to move the plane inclusion boundary and to create an incremental region of eigenstrain is defined analogously to the one for moving cracks and dislocations and represents the mechanical rate of work needed to be provide for the expansion of the inclusion. The calculated value, which is the “self-force” of the expanding inclusion, has a static component plus a dynamic one depending only on the current value of the velocity, while in the case of the spherical boundary, there is an additional contribution accounting for the jump in the strain at the farthest part at the back of the inclusion having the time to reach the front boundary, thus making the dynamic “self-force” history dependent.     

可压超弹性-塑性材料中的空穴生成

本文研究了材料的弹塑性性质对球体中空穴生成问题的影响，材料的弹性用一种可压超弹性材料的本构关系来描述，材料的塑性用满足材料的不可压条件和Tresca屈服条件的理想塑性材料的本构关系来描述。这类超弹性．塑性材料中可以发生空穴的生成现象，得到了在表面拉伸作用下球体中空穴生成时空穴半径与临界拉伸之间的关系式和临界拉伸。球体的变形可分为弹-塑性变形阶段和完全塑性变形阶段，球体中心首先形成塑性变形区域，并有空穴的突然生成；塑性变形区域能够快速增长，并且使球体很快进入完全塑性变形阶段；空穴在弹-塑性变形阶段迅速增长，但进入完全塑性变形阶段后增长较慢。同时给出了不同变形阶段球体中的应力分布。数值结果表明材料的塑性性质对材料中的空穴生成有明显的影响。     

Dynamic void nucleation and growth in solids: A self-consistent statistical theory

We present a framework for a self-consistent theory of spall fracture in ductile materials, based on the dynamics of void nucleation and growth. The constitutive model for the material is divided into elastic and “plastic” parts, where the elastic part represents the volumetric response of a porous elastic material, and the “plastic” part is generated by a collection of representative volume elements (RVEs) of incompressible material. Each RVE is a thick-walled spherical shell, whose average porosity is the same as that of the surrounding porous continuum, thus simulating void interaction through the resulting lowered resistance to further void growth. All voids nucleate and grow according to the appropriate dynamics for a thick-walled sphere made of incompressible material. The macroscopic spherical stress in the material drives the response in all volume elements, which have a distribution of critical stresses for void nucleation, and the statistically weighted sum of the void volumes of all RVEs generates the global porosity. Thus, macroscopic pressure, porosity, and a distribution of growing microscopic voids are fully coupled dynamically. An example is given for a rate-independent, perfectly plastic material. The dynamics of void growth gives rise to a rate effect in the macroscopic material even though the parent material is rate independent.     

Stress state of a transversely isotropic spherical shell with a rigid circular inclusion

Analytical expressions are derived for the stresses near a rigid circular inclusion in a transversely isotropic shallow spherical shell under uniform pressure. The form of solution depends on the range of the transverse shear compliance parameter. The influence of the relative radius of a rigid inclusion and transverse shear compliance on stress concentration is analyzed __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 12, pp. 67–73, December 2005.     

直角平面区域内固定圆形刚性夹杂问题的Green函数解

利用复变函数法、多极坐标移动技术研究了直角平面区域内含有固定圆形夹杂时的反平面问题Green函数解.首先构造出不含夹杂的完整直角平面区域内满足边界应力条件的入射位移场;其次,建立直角平面区域内固定圆形夹杂对该入射场产生的满足直角边界应力自由条件的散射波解,并由叠加原理得到介质内的总波场.最后利用夹杂边界处的位移条件确定出散射波解中的未知系数,最终得到问题的Green函数解,还通过算例讨论了夹杂边界处的径向应力和环向应力随不同波数、角度和不同夹杂位置及不同点源位置的变化情况.算例结果表明了该文方法的有效实用性.     

Onset of residual stress fields near solitary spherical inclusions in a perfectly elastoplastic medium

When computing residual stresses in deformable solids, one has to use the theory of elastoplastic solids, because the final level and distribution of residual stresses is determined exactly by the accumulated reversible strains. In turn, to compute the elastic strains, one needs to determine the displacement field. The problem of determining displacements in statically determinate problems of the theory of perfect elastoplastic solids was considered for the first time in [1, 2]. The techniques proposed there permitted solving the problem of finding the residual stresses near a cylindrical cavity in a perfectly elastoplastic medium [3]. It was shown that secondary plastic flow [4] may arise in the unloading processes, which significantly redistributes the final residual stresses. In the present paper, we consider the loading and unloading problems for a ball with a rigid or elastic spherical inclusion. We study the onset of secondary plastic flow under unloading and compute the residual stresses. Thus, we model the onset of the residual stress field near a more rigid inhomogeneity. The case of a softer inhomogeneity was essentially considered in [3], where the onset of the residual stress field near a continuity flaw was studied.     

Rolling contact of a rigid sphere/sliding of a spherical indenter upon a viscoelastic half-space containing an ellipsoidal inhomogeneity

In this paper, the frictionless rolling contact problem between a rigid sphere and a viscoelastic half-space containing one elastic inhomogeneity is solved. The problem is equivalent to the frictionless sliding of a spherical tip over a viscoelastic body. The inhomogeneity may be of spherical or ellipsoidal shape, the later being of any orientation relatively to the contact surface. The model presented here is three dimensional and based on semi-analytical methods. In order to take into account the viscoelastic aspect of the problem, contact equations are discretized in the spatial and temporal dimensions. The frictionless rolling of the sphere, assumed rigid here for the sake of simplicity, is taken into account by translating the subsurface viscoelastic fields related to the contact problem. Eshelby's formalism is applied at each step of the temporal discretization to account for the effect of the inhomogeneity on the contact pressure distribution, subsurface stresses, rolling friction and the resulting torque. A Conjugate Gradient Method and the Fast Fourier Transforms are used to reduce the computation cost. The model is validated by a finite element model of a rigid sphere rolling upon a homogeneous vciscoelastic half-space, as well as through comparison with reference solutions from the literature. A parametric analysis of the effect of elastic properties and geometrical features of the inhomogeneity is performed. Transient and steady-state solutions are obtained. Numerical results about the contact pressure distribution, the deformed surface geometry, the apparent friction coefficient as well as subsurface stresses are presented, with or without heterogeneous inclusion.     

A coated rigid elliptical inclusion loaded by a couple in the presence of uniform interfacial and hoop stresses

We consider a confocally coated rigid elliptical inclusion, loaded by a couple and introduced into a remote uniform stress field. We show that uniform interfacial and hoop stresses along the inclusion–coating interface can be achieved when the two remote normal stresses and the remote shear stress each satisfy certain conditions. Our analysis indicates that: (i) the uniform interfacial tangential stress depends only on the area of the inclusion and the moment of the couple; (ii) the rigid-body rotation of the rigid inclusion depends only on the area of the inclusion, the coating thickness, the shear moduli of the composite and the moment of the couple; (iii) for given remote normal stresses and material parameters, the coating thickness and the aspect ratio of the inclusion are required to satisfy a particular relationship; (iv) for prescribed remote shear stress, moment and given material parameters, the coating thickness, the size and aspect ratio of the inclusion are also related. Finally, a harmonic rigid inclusion emerges as a special case if the coating and the matrix have identical elastic properties.     

Convective cooling/heating induced thermal stresses in a fluid saturated porous medium undergoing local thermal non-equilibrium

Local thermal non-equilibrium (LTNE) may have profound effects on the pore pressure and thermal stresses in fluid saturated porous media under transient thermal loads. This work investigates the temperature, pore pressure, and thermal stress distributions in a porous medium subjected to convective cooling/heating on its boundary. The LTNE thermo-poroelasticity equations are solved by means of Laplace transform for two fundamental problems in petroleum engineering and nuclear waste storage applications, i.e., an infinite porous medium containing a cylindrical hole or a spherical cavity subjected to symmetrical thermo-mechanical loads on the cavity boundary. Numerical examples are presented to examine the effects of LTNE under convective cooling/heating conditions on the temperature, pore pressure and thermal stresses around the cavities. The results show that the LTNE effects become more pronounced when the convective heat transfer boundary conditions are employed. For the cylindrical hole problem of a sandstone formation, the thermally induced pore pressure and the magnitude of thermal stresses are significantly higher than the corresponding values in the classical poroelasticity, which is particularly true under convective cooling with moderate Biot numbers. For the spherical cavity problem of a clay medium, the LTNE effect may become significant depending on the boundary conditions employed in the classical theory.     

An extension of generalized plane stress for problems with out-of-plane restraint

The standard concept of generalized plane stress is extended to obtain a new mathematical model for studying the effect of local out-of-plane displacement restraint on the in-plane stresses and displacements in thin plates. It is pointed out how this model may be used by the photoelastician, whose otherwise plane-stress experiment introduces an unavoidable out-of-plane restraint condition in the model, to obtain some estimate of the deviation to be expected between the results of his experiment and the actual plane-stress solution of the problem. In this way, the model may be applied to aid in the interpretation of a large class of two-dimensional photoelastic analyses involving the determination of stresses near rigid inclusions and rigid boundaries. The extended model is then applied to the problem of an annular disk subjected to thermal shrinkage and completely restrained at its outer boundary. In view of the simplicity of the model, the predicted radial and circumferential stress distributions agree remarkably well with existing photoelastic data. In contrast, results obtained from standard generalized plane-stress theory, which cannot account for the out-of-plane displacement restraint at the outer boundary, show substantial deviation from experimental values, especially near the restrained boundary.     

A rigid triangular inclusion partially bonded in an elastic matrix

The two-dimensional problem of a rigid rounded-off angle triangular inclusion partially bonded in an infinite elastic plate is studied. The unbonded part of the inclusion boundary forms an interfacial crack. Based on the complex variable method for curvilinear boundaries, the problem is reduced to a non-homogeneous Hilbert problem and the stress and displacement fields in the plate are obtained in closed form. Special attention is paid in the investigation of the stress field in the vicinity of the crack tip. It is found that the stresses present an oscillatory singularity and the general equations for the local stresses are derived. The singular stress field is coupled with the maximum circumferential stress and the minimum strain energy density criteria to study the fracture characteristics of the composite plate. Results are given for the complex stress intensity factors, the local stresses, the crack extension angles and the critical applied loads for unstable crack growth from its more vulnerable tip or two types of interfacial cracks along the inclusion boundary.     

Slow Motion of a Porous Sphere Translating Along the Axis of a Circular Cylindrical Pore Subject to a Stress Jump Condition

The coupled flow problem of an incompressible axisymmetrical quasisteady motion of a porous sphere translating in a viscous fluid along the axis of a circular cylindrical pore is discussed using a combined analytical–numerical technique. At the fluid–porous interface, the stress jump boundary condition for the tangential stress along with continuity of normal stress and velocity components are employed. The flow through the porous particle is governed by the Brinkman model and the flow in the outside porous region is governed by Stokes equations. A general solution for the field equations in the clear region is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are satisfied first at the cylindrical pore wall by the Fourier transforms and then on the surface of the porous particle by a collocation method. The collocation solutions for the normalized hydrodynamic drag force exerted by the clear fluid on the porous particle is calculated with good convergence for various values of the ratio of radii of the porous sphere and pore, the stress jump coefficient, and a coefficient that is proportional to the permeability. The shape effect of the cylindrical pore on the axial translation of the porous sphere is compared with that of the particle in a spherical cavity; it found that the porous particle in a circular cylindrical pore in general attains a lower hydrodynamic drag than in a spherical envelope.     

On Structured Surfaces with Defects: Geometry,Strain Incompatibility,Stress Field,and Natural Shapes

In this paper we study the stress and deformation fields generated by nonlinear inclusions with finite eigenstrains in anisotropic solids. In particular, we consider finite eigenstrains in transversely isotropic spherical balls and orthotropic cylindrical bars made of both compressible and incompressible solids. We show that the stress field in a spherical inclusion with uniform pure dilatational eigenstrain in a spherical ball made of an incompressible transversely isotropic solid such that the material preferred direction is radial at any point is uniform and hydrostatic. Similarly, the stress in a cylindrical inclusion contained in an incompressible orthotropic cylindrical bar is uniform hydrostatic if the radial and circumferential eigenstrains are equal and the axial stretch is equal to a value determined by the axial eigenstrain. We also prove that for a compressible isotropic spherical ball and a cylindrical bar containing a spherical and a cylindrical inclusion, respectively, with uniform eigenstrains the stress in the inclusion is uniform (and hydrostatic for the spherical inclusion) if the radial and circumferential eigenstrains are equal. For compressible transversely isotropic and orthotropic solids, we show that the stress field in an inclusion with uniform eigenstrain is not uniform, in general. Nevertheless, in some special cases the material can be designed in order to maintain a uniform stress field in the inclusion. As particular examples to investigate such special cases, we consider compressible Mooney-Rivlin and Blatz-Ko reinforced models and find analytical expressions for the stress field in the inclusion.     

A STUDY OF THE CATASTROPHE AND THE CAVITATION FOR A SPHERICAL CAVITY IN HOOKE’S MATERIAL WITH 1/2 POISSON’S RATIO

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ｔｈｅｒｅｓｅａｒｃｈｅｓｏｎｃａｖｉｔａｔｉｏｎａｎｄｃａｔａｓｔｒｏｐｈｅｏｆａｃａｖｉｔｙｈａｖｅｓｕｐｐｌｉｅｄａｎｅｗｍｅｔｈｏｄｆｏｒｉｎｖｅｓｔｉｇａｔｉｎｇｔｈｅｍｅｃｈａｎｉｃｓｏ...     

DECAGONAL QU ASICRYSTALLINE ELLIPTICAL INCLUSIONS UNDER THERMOMECHANICAL LOADING

In this work, an elegant method is proposed to derive the thermoelastic field in- duced by thermomechanical loadings in a decagonal quasicrystalline composite composed of an infinite matrix reinforced by an elliptical inclusion. The thermomechanical loadings include a uniform temperature change, remote uniform in-plane heat fluxes and remote uniform in-plane stresses. The corresponding boundary value problem is ultimately reduced to the solution of two independent sets of four coupled linear algebraic equations, each of which involves four complex constants characterizing the internal stress field. The solution demonstrates that a uniform tem- perature change and remote uniform stresses will induce an internal uniform stress field, and that uniform heat fluxes will result in a linearly distributed internal stress field within the elliptical inclusion. The induced uniform rigid body rotation within the inclusion is given explicitly.     

Spherical ceramic pebbles subjected to multiple non-concentrated surface loads

This paper presents an analytical solution for the stress distributions within spherical ceramic pebbles subjected to multiple surface loads along different directions. The method of solution employs a displacement approach together with the Fourier associated Legendre expansion for piecewise boundary loads. The solution corresponds to spherically isotropic elastic spheres. The classical solution for isotropic spheres subjected diametral point loads is recovered as a special case of our solution. For the isotropic pebbles under consideration, stresses within spheres are numerically evaluated. The results show that the number of loads does have significant influence on the maximum tensile stress inside the sphere. Moreover, the applicability of solutions using the series expansion method for stresses near surface load areas is also examined. The stresses evaluated with large enough number of terms agree quite well with those derived from FEM simulations, except around the edge of circle load area.     

Cavity formation at the center of a sphere composed of two compressible hyper-elastic materials

IntroductionHorgan[1] reviewedthecavitatedbifurcationproblemforhyper_elasticmaterials,includinginhomogeneousandanisotropicmaterialsaswellashomogeneousandisotropicmaterials .Forincompressiblematerials,HorganandPence[2 ,3 ] examinedtheeffectofmaterialinhomogeneityontheformationandgrowthofvoidandobtainedananalyticsolutionofthecavitatedbifurcationproblemforasolidspherecomposedoftwoneo_Hookeanmaterials.Thebifurcationmayoccurnotonlytotherightbutalsototheleftforthecomposedsphere .Thestabilitiesofth…     

The Cosserat Subspace ũ(−1) for Bodies of Spherical Geometry

We construct the orthogonal bases of the Cosserat eigenvectors ũ⁽⁻¹⁾ for the first boundary value problem of an elastic solid sphere and an infinite elastic space containing a spherical rigid inclusion. These orthogonal bases are expressed in terms of the Jacobi and Legendre polynomials. An example of a nonharmonic heat source shows the convergence of the sequence of the eigenvectors ũ⁽⁻¹⁾. This revised version was published online in August 2006 with corrections to the Cover Date.     

Uncoupled magnetoelastic problem for a ferromagnetic body with a spherical cavity

An uncoupled stress problem for an unbounded elastic soft ferromagnetic body with a spherical cavity in a magnetic field uniform at infinity is solved. The stresses, displacements, and magnetic quantities in the body are determined. The features of stress distribution over the body and its boundary surface are studied __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 10, pp. 42–48, October 2007.     

Bifurcation phenomena in the rigid inclusion power law matrix composite sphere

Bifurcation of interface separation related to cavity nucleation is analyzed for a radially loaded composite sphere consisting of a rigid inclusion separated from a power law matrix by a uniform, non-linear cohesive zone. Equations for the spherically symmetric and non-symmetric problems are obtained from a hyperelastic finite strain theory by a limiting process that preserves non-linear matrix and interface response at infinitesimal strain. A complete solution to the symmetric problem is presented including bifurcation load, stresses, and evolution of elasto-plastic boundary and interface separation. An analysis of non-symmetric bifurcation, under symmetric conditions of geometry and loading, yields the bifurcation load and first non-symmetric mode shape associated with rigid inclusion displacement. An energy analysis is carried out for both symmetric and non-symmetric problems in order to assess stability of spherically symmetric states to spherically symmetric and non-symmetric “rigid body mode” perturbations.Results are provided for an interface force law that captures interface failure in normal mode and linear response in shear mode. For the symmetric problem, (i) there are threshold parameter values above which bifurcation will generally not occur, (ii) threshold values below which there do not exist equilibria in the post bifurcation regime, (iii) bifurcation occurs after attainment of the maximum interface strength. For the non-symmetric problem, (i) bifurcation always occurs, although it can be delayed by interfacial shear, (ii) for the smooth interface, non-symmetric bifurcation occurs after attainment of the maximum interface strength and always precedes symmetric bifurcation.