Application of results on Eshelby tensor to the determination of effective poroelastic properties of anisotropic rocks-like composites

The present work is devoted to the determination of the macroscopic poroelastic properties of anisotropic elastic porous materials saturated by a fluid under pressure. It makes use of the theoretical results provided by Withers [Withers, P.J., 1989. The determination of the elastic field of an ellipsoidal inclusion in a transversely isotropic medium, and its relevance to composite materials. Philosophical Magazine A 59 (4), 759–781.] for the problem of an ellipsoidal inclusion embedded in a transversely isotropic elastic medium. The particular case of a spherical inclusion is very important for rock-like composites such as argillite and shales. The implementation of these results in a micromechanical theory of poroelasticity allows to quantify the effects of the solid matrix anisotropy and of pore space on the effective poromechanical properties. Closed form expressions of Biot tensor and of Biot modulus are presented as well as numerical applications for anisotropic shales.     

Elastodynamic wave scattering by an incompressible elastic inclusion

The problem of scattering of time-harmonic elastodynamic waves by an incompressible elastic inclusion is solved by means of the null field approach. The solution is obtained both directly and as a limit of the solution to the corresponding problem for a compressible inclusion. It is also demonstrated that the null field solution to the problem of scattering by a rigid movable scatterer can be obtained from a null field solution for the incompressible scatterer by taking the limit of infinite shear modulus. Some numerical results for spherical and spheroidal inclusions are given.     

A unified treatment of the elastic elliptical inclusion under antiplane shear

Summary A generalized and unified treatment is presented for the antiplane problem of an elastic elliptical inclusion undergoing uniform eigenstrains and subjected to arbitrary loading in the surrounding matrix. The general solution to the problem is obtained through the use of conformal mapping technique and Laurent series expansion of the associated complex potentials. The resulting elastic fields are derived explicitly in both transformed and physical planes for the inclusion and the surrounding matrix. These relations are universal in the sense of being independent of any particular loading as well as the geometry of the matrix. The complete field solutions are provided for an elliptical inclusion under uniform loading at inifinity, and for a screw dislocation interacting with the elastic elliptical inclusion.     

Stress concentrations in the particulate composite with transversely isotropic phases

The accurate series solution have been obtained of the elasticity theory problem for a transversely isotropic solid containing a finite or infinite periodic array of anisotropic spherical inclusions. The method of solution has been developed based on the multipole expansion technique. The basic idea of method consists in expansion the displacement vector into a series over the set of vectorial functions satisfying the governing equations of elastic equilibrium. The re-expansion formulae derived for these functions provide exact satisfaction of the interfacial boundary conditions. As a result, the primary spatial boundary-value problem is reduced to an infinite set of linear algebraic equations. The method has been applied systematically to solve for three models of composite, namely a single inclusion, a finite array of inclusions and an infinite periodic array of inclusions, respectively, embedded in a transversely isotropic solid. The numerical results are presented demonstrating that elastic properties mismatch, anisotropy degree, orientation of the anisotropy axes and interactions between the inclusions can produce significant local stress concentration and, thus, affect greatly the overall elastic behavior of composite.     

Cloaking of a circular cylindrical elastic inclusion from antiplane elastic waves and resonance effects

Cloaking of a circular cylindrical elastic inclusion embedded in a homogeneous linear isotropic elastic medium from antiplane elastic waves is studied. The transformation or change-of-variables method is used to determine the material properties of the cloak and the homogenization theory of composites is used to construct a multilayered cloak consisting of many bi-material cells. The large system of algebraic equations associated with this problem is solved by using the concept of multiple scattering with wave expansion coefficient matrices. Numerical results for cloaking of an elastic inclusion and a rigid inclusion are compared with the case of a cavity. It is found that while the cloaking patterns for the three cases are similar, the major difference is that standing waves are generated in the elastic inclusion and the multilayered cloak cannot prevent the motion inside the elastic inclusion, even though the cloak seems nearly perfect. Waves can penetrate into and cause vibrations inside the elastic inclusion, where the amplitude of standing waves depend on the material properties of the inclusion but are very much reduced when compared to the case when there is no cloak. For a prescribed mass density, the displacements inside the elastic cylinder decrease as the shear modulus increases. Moreover, the cloaking of the elastic inclusion over a range of wavenumbers is also investigated. There is significant low frequency scattering even if the cloak consists of a large number of layers. When the wavenumber increases, the multilayered cloak is not effective if the cloak consists of an insufficient number of layers. Resonance effects that occur in cloaking of elastic inclusions are also discussed.     

The rigid movable inclusion in elastostatics and elastodynamics

The problem of a single rigid movable inclusion is solved both for elastostatics and elastodynamic, using the null field approach. It is also shown how the solution can be obtained as the limit of the solution for an elastic inclusion. Numerical results for scattering by a superspheroidal inclusion are given.     

Scattering by an anisotropic circle

The scattering by a circle is considered when the outside medium is isotropic and the inside medium is anisotropic (orthotropic). The problem is a scalar one and is phrased as a scattering problem for elastic waves with polarization out of the plane of the circle (SH wave), but the solution is with minor modifications valid also for scattering of electromagnetic waves. The equation inside the circle is first transformed to polar coordinates and it then explicitly contains the azimuthal angle through trigonometric functions. Making an expansion in a trigonometric series in the azimuthal coordinate then gives a coupled system of ordinary differential equations in the radial coordinate that is solved by power series expansions. With the solution inside the circle complete the scattering problem is solved essentially as in the classical case. Some numerical examples are given showing the influence of anisotropy, and it is noted that the effects of anisotropy are generally strong except at low frequencies where the dominating scattering only depends on the mean stiffness and not on the degree of anisotropy.     

Scattering of harmonic waves by a disk-shaped rigid inclusion in a three-dimensional elastic matrix

We consider a three-dimensional problem on the interaction of harmonic waves with a thin rigid movable inclusion in an infinite elastic body. The problem is reduced to solving a system of two-dimensional boundary integral equations of Helmholtz potential type for the stress jump functions on the opposite surfaces of the inclusion. We propose a boundary element method for solving the integral equations on the basis of the regularization of their weakly singular kernels. Using the asymptotic relations between the amplitude-frequency characteristics of the wave farzone field and the obtained boundary stress jump functions, we determine the amplitudes of the shear plane wave scattering by a circular disk-shaped inclusion for various directions of the wave incident on the inclusion and for a broad range of wave numbers.     

Elliptical inclusion in an anisotropic plane: non-uniform interface effects

We study the plane deformation of an elastic composite system made up of an anisotropic elliptical inclusion and an anisotropic foreign matrix surrounding the inclusion. In order to capture the influence of interface energy on the local elastic field as the size of the inclusion approaches the nanoscale, we refer to the Gurtin-Murdoch model of interface elasticity to describe the inclusion-matrix interface as an imaginary and extremely stiff but zero-thickness layer of a finite stretching modulus. As opposed to isotropic cases in which the effects of interface elasticity are usually assumed to be uniform (described by a constant interface stretching modulus for the entire interface), the anisotropic case considered here necessitates non-uniform effects of interface elasticity (described by a non-constant interface stretching modulus), because the bulk surrounding the interface is anisotropic. To this end, we treat the interface stretching modulus of the anisotropic composite system as a variable on the interface curve depending on the specific tangential direction of the interface. We then devise a unified analytic procedure to determine the full stress field in the inclusion and matrix, which is applicable to the arbitrary orientation and aspect ratio of the inclusion, an arbitrarily variable interface modulus, and an arbitrary uniform external loading applied remotely. The non-uniform interface effects on the external loading-induced stress distribution near the interface are explored via a group of numerical examples. It is demonstrated that whether the nonuniformity of the interface effects has a significant effect on the stress field around the inclusion mainly depends on the direction of the external loading and the aspect ratio of the inclusion.     

A novel hybrid finite element analysis of inplane singular elastic field around inclusion corners in elastic media

This paper deals with the inplane singular elastic field problems of inclusion corners in elastic media by an ad hoc hybrid-stress finite element method. A one-dimensional finite element method-based eigenanalysis is first applied to determine the order of singularity and the angular dependence of the stress and displacement field, which reflects elastic behavior around an inclusion corner. These numerical eigensolutions are subsequently used to develop a super element that simulates the elastic behavior around the inclusion corner. The super element is finally incorporated with standard four-node hybrid-stress elements to constitute an ad hoc hybrid-stress finite element method for the analysis of local singular stress fields arising from inclusion corners. The singular stress field is expressed by generalized stress intensity factors defined at the inclusion corner. The ad hoc finite element method is used to investigate the problem of a single rectangular or diamond inclusion in isotropic materials under longitudinal tension. Comparison with available numerical results shows the present method is an efficient mesh reducer and yields accurate stress distribution in the near-field region. As applications, the present ad hoc finite element method is extended to discuss the inplane singular elastic field problems of a single rectangular or diamond inclusion in anisotropic materials and of two interacting rectangular inclusions in isotropic materials. In the numerical analysis, the generalized stress intensity factors at the inclusion corner are systematically calculated for various material type, stiffness ratio, shape and spacing position of one or two inclusions in a plate subjected to tension and shear loadings.     

Scattering far field solution of SH-wave by movable rigid cylindrical interface inclusion

The Green's function is used to solve the scattering far field solution of SH-wave by a movable rigid cylindrical interface inclusion in a linear elastic body. First, a suitable Green's function is developed, which is the fundamental displacement solution of an elastic half space with a movable rigid half-cylindrical inclusion impacted by out-of-plane harmonic line source loaded at any point of its horizontal surface. By using the Green's function, a series of Fredholm integral equations of the first kind which determine the scattering far field can be set up. Then the paper gives the expressions on the far field including the displacement mode of scattering wave and the solution of scattering cross-section. Finally, some examples and numerical results are discussed to analyze the influence of the combination of different media parameters on the answer of far field.     

直角平面内弹性圆夹杂对入射平面SH波的散射

利用复变函数法、多极坐标移动技术及傅立叶级数展开求解二维直角平面内圆形弹性夹杂对稳态入射平面SH波的散射问题。首先写出直角平面内不含夹杂时的入射波场和反射波场;其次建立直角平面内含夹杂时夹杂外的散射波解和夹杂内的驻波解,并利用叠加原理写出问题的总波场,借助夹杂边界处应力和位移的连续条件建立求解散射波解和驻波解中未知系数的无穷代数方程组并求解,通过算例具体讨论了直角平面水平边界点的位移幅度比和夹杂边界处径向应力集中系数随不同无量纲波数、入射角及圆孔位置的变化情况,结果表明了算法的有效实用性。     

各向异性介质的弹性波散射问题的边界元方法

本文引用加权残数法建立了各向异性介质内含任意形式异质夹杂时的散射问题的边界积分方程式,导出了相应的辐射条件,计算了内含圆柱体,椭圆柱体、界面裂纹情形下对SH 波的散射位移场、应力场以及散射横截面.数值结果表明本方法用于解答各向异性介质的弹性波散射问题具有良好的精度和应用前景.     

Stress State of a Transversely Isotropic Medium with an Arbitrarily Oriented Spheroidal Inclusion

The stress-concentration problem for an elastic transversely isotropic medium containing an arbitrarily oriented spheroidal inclusion (inhomogeneity) is solved. The stress state in the elastic space is represented as the superposition of the principal state and the perturbed state due to the inhomogeneity. The problem is solved using the equivalent-inclusion method, the triple Fourier transform in space variables, and the Fourier-transformed Green function for an infinite anisotropic medium. Double integrals over a finite domain are evaluated using the Gaussian quadrature formulas. In special cases, the results are compared with those obtained by other authors. The influence of the geometry and orientation of the inclusion and the elastic properties of the medium and inclusion on the stress concentration is studied__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 2, pp. 33–40, February 2005.     

The elastic field caused by a general ellipsoidal inclusion and the application to martensite formation 2

The elastic field throughout an ellipsoidal inclusion in an indefinitely-extended anisotropic material is investigated when an eigenstrain (a stress-free transformation strain) is periodically distributed throughout the inclusion. This is an extention of the results obtained by J.D. Eshelby (1961) for uniform eigenstrains and by R.J. Asaro and D.M. Barnett (1975) for polynomial eigenstrains. The solution is applied to the evaluation of elastic strain energies of a disc-shaped martensite with alternating twins and of a spherical precipitate with a banded structure. The significant amount of the elastic strain energies explains the necessity of the supercooling of austenite steel for the martensitic transformation to occur.     

Multi-Scale Asymptotic Expansion for a Small Inclusion in Elastic Media

The aim of this paper is to present an asymptotic expansion of the influence of a small inclusion of different stiffness in an elastic media. The applicative interest of this study is to provide tools which take into account this influence and correct the deformation without inclusion by additive terms that can be precalculated and which depend only on the shape of the inclusion. We treat two problems: an anti-plane linearized elasticity problem and a plane strain problem. On every expansion order we provide corrective terms modeling the influence of the inclusion using techniques of scaling and multi-scale asymptotic expansions. The resulting expansion is validated by comparing it to a test case obtained by solving the Poisson transmission problem in the case of an inclusion of circular shape using the separation of variables method. Proofs of existence and uniqueness of our fields on unbounded domains are also adapted to the bidimensional Poisson problem and the linear elasticity problem.     

On interfacial discontinuities in elastic composites

The paper contains a succinct analysis of interfacial discontinuities in anisotropic elastic solids. The results are combined with known results on the ellipsoidal inclusion problem to provide some general formulae for the determination of stress (and strain) concentration factors. Some explicit results are given for cavities in an infinite matrix under arbitrary uniform loading at infinity.     

二维直角平面内固定圆形夹杂对稳态入射反平面剪切波的散射

利用复变函数法、多极坐标及傅立叶级数展开技术求解了二维直角平面内固定圆形夹杂对稳态入射反平面剪切（shearing horizontal, SH）波的散射问题。首先构造出介质内不存在夹杂时的入射波场和反射波场,然后建立介质内存在夹杂时由夹杂边界产生的能够自动满足直角边应力自由条件的散射波解,从而利用叠加原理写出介质内的总波场。利用夹杂边界处位移条件和傅立叶级数展开方法列出求解散射波中未知系数的无穷代数方程组,在满足计算精度的前提下通过有限项截断,得到相应有限代数方程组的解,最后通过算例具体讨论了二维直角平面水平边界点的位移幅度比和相位随量纲一波数、入射波入射角及夹杂位置的不同而变化的情况,结果表明了算法的有效实用性。     

Elastic properties of inhomogeneous transversely isotropic rocks

The problem to determine the effective elastic moduli and velocities of elastic wave propagation in transversely isotropic solid containing aligned spheroidal inhomogeneities (solid grains, vugs and micro-cracks) has been solved using the self-consistent scheme known as effective medium approximation (EMA). Since a solution of so-called one-particle problem is a base for each self-consistent method, we solved this problem as a first step for spheroidal inhomogeneity in a transversely isotropic medium. In contrast to the known solution of this problem by Lin and Mura we obtained the expressions for the strain field inside inclusion in the explicit form (without quadratures). The obtained solution was used then in the symmetric variant of the EMA where each component of the system was considered as spheroid with its own aspect ratio. This approach was applied to simulate the properties of the rocks containing isolated pores and micro-cracks. For connected fluid-filled pores we used the anisotropic variant of the Gassmann theory. The results of the calculations, obtained for the effective elastic moduli, have been compared with the experimental data and theoretical simulations of the other authors. Unlike many other rock mechanics theories, EMA approximation gives correct elastic moduli values even in the nondilute concentration of inhomogeneities. The comparison of the experimental data for oriented crack system with the EMA predictions indicates their good correspondence.     

A boundary element method for calculating elastic waves scattered by axisymmetric inclusion

A Boundary Element Method (BEM) is described to compute the scattering of elastic waves by an axisymmetric inclusion in an infinite elastic medium. The boundary loads applied to the inclusion is expanded in terms of Fourier series in an infinite space. The boundary integral equation is solved in the general direction of the axisymmetric inclusion by BEM. The problem of the 3-D scattering of elastic waves is reduced to a 1-Done. According to the geometric features of the axisymmetric in clusion the ring shell elements are adopted in this method. A comparison is made with other BEM methods. The numerical results show this method can reduce the amount of calculation and enhance the speed of convergence. Supported by Foundation of Ph. D Program of State Education Commission of China