New Green’s function for stress field and a note of its application in quantum-wire structures

The elastic field caused by the lattice mismatch between the quantum wires and the host matrix can be modeled by a corresponding two-dimensional hydrostatic inclusion subjected to plane strain conditions. The stresses in such a hydrostatic inclusion can be effectively calculated by employing the Green’s functions developed by Downes and Faux, which tend to be more efficient than the conventional method based on the Green’s function for the displacement field. In this study, Downes and Faux’s paper is extended to plane inclusions subjected to arbitrarily distributed eigenstrains: an explicit Green’s function solution, which evaluates the stress field due to the excitation of a point eigenstrain source in an infinite plane directly, is obtained in a closed-form. Here it is demonstrated that both the interior and exterior stress fields to an inclusion of any shape and with arbitrarily distributed eigenstrains are represented in a unified area integral form by employing the derived Green’s functions. In the case of uniform eigenstrain, the formulae may be simplified to contour integrals by Green’s theorem. However, special care is required when Green’s theorem is applied for the interior field. The proposed Green’s function is particularly advantageous in dealing numerically or analytically with the exterior stress field and the non-uniform eigenstrain. Two examples concerning circular inclusions are investigated. A linearly distributed eigenstrain is attempted in the first example, resulting in a linear interior stress field. The second example solves a circular thermal inclusion, where both the interior and exterior stress fields are obtained simultaneously.     

Eshelby's problem of non-elliptical inclusions

The Eshelby problem consists in determining the strain field of an infinite linearly elastic homogeneous medium due to a uniform eigenstrain prescribed over a subdomain, called inclusion, of the medium. The salient feature of Eshelby's solution for an ellipsoidal inclusion is that the strain tensor field inside the latter is uniform. This uniformity has the important consequence that the solution to the fundamental problem of determination of the strain field in an infinite linearly elastic homogeneous medium containing an embedded ellipsoidal inhomogeneity and subjected to remote uniform loading can be readily deduced from Eshelby's solution for an ellipsoidal inclusion upon imposing appropriate uniform eigenstrains. Based on this result, most of the existing micromechanics schemes dedicated to estimating the effective properties of inhomogeneous materials have been nevertheless applied to a number of materials of practical interest where inhomogeneities are in reality non-ellipsoidal. Aiming to examine the validity of the ellipsoidal approximation of inhomogeneities underlying various micromechanics schemes, we first derive a new boundary integral expression for calculating Eshelby's tensor field (ETF) in the context of two-dimensional isotropic elasticity. The simple and compact structure of the new boundary integral expression leads us to obtain the explicit expressions of ETF and its average for a wide variety of non-elliptical inclusions including arbitrary polygonal ones and those characterized by the finite Laurent series. In light of these new analytical results, we show that: (i) the elliptical approximation to the average of ETF is valid for a convex non-elliptical inclusion but becomes inacceptable for a non-convex non-elliptical inclusion; (ii) in general, the Eshelby tensor field inside a non-elliptical inclusion is quite non-uniform and cannot be replaced by its average; (iii) the substitution of the generalized Eshelby tensor involved in various micromechanics schemes by the average Eshelby tensor for non-elliptical inhomogeneities is in general inadmissible.     

具有光滑界面的椭球夹杂的轴对称弹性场

本文在旋转椭球坐标系下,利用Papkovich—Neuber位移通解求解了具有光滑界面椭球夹杂由于均匀的特征应变引起的轴对称弹性场,与理想界面不同,在夹杂与基体界面不能经受剪应力而可自由滑动的情况下,解答只能是无穷级数形式,因此文中给出了数值算例。     

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.     

The axisymmetrical elastic field of an ellipsoidal inclusion with a slipping interface

This paper is concerned with the axisymmetrical elastic fields caused by an ellipsoidal inclusion with a slipping interface which undergoes a uniform eigenstrain. The problem is solved under a revolving ellipsoidal coordinate with the aid of Papkovich-Neuber general dipacement formula. In contrast to the perfectly bonded interface, when the interface between the inclusion and the matrix cannot sustain shear stress, and is free to slip, the solution cannot be expressed in closed form and involves infinite series. Therefore, the results are illustrated by numerical examples.     

任意形状热夹杂位移场的三角形单元离散算法

平面夹杂模型在纤维增强型复合材料中有广泛应用.复合材料内部通常含有不规则形状夹杂,而夹杂物的存在能严重影响材料的机械力学性能,往往导致应力集中及裂纹萌生等失效先兆.先前关于多边形夹杂的研究大多数关注受均匀本征应变下的应力/应变解,而对位移的分析较少. 基于格林函数方法和围道积分,本文给出了平面热夹杂边界线单元的封闭解析解,可方便应用于受任意分布本征应变的任意形状平面热夹杂位移场的数值计算.当夹杂受均匀本征应变时, 只需将该夹杂边界进行一维离散,因而本文方法可直接得出受均匀分布热本征应变的任意多边形夹杂位移场的封闭解析解.当夹杂区域存在非均匀分布本征应变时,可将该区域划分为足够小的三角形单元进行数值计算. 众所周知,应力应变场在多边形夹杂顶点处具有奇异性,容易导致数值计算上的处理困难及相应的数值稳定性问题; 然而本文工作表明,在多边形顶点处位移场是连续有界的, 因而数值稳定性较好.本文算法可以便捷高效地通过计算机编程实现. 文中给出的验证算例,均体现了本文离散方法的高精度、以及计算编程的鲁棒性.      

The Eshelby property of sliding inclusions

The present paper solves the problem of sliding ellipsoidal/elliptical inclusion with the Eshelby property. Results show that the sliding ellipsoidal/elliptical inclusions can have uniform eigenstresses if the prescribed uniform eigenstrains fulfill certain prerequisites. Solutions and prerequisites are obtained for ellipsoidal and elliptical inclusions, respectively. It is shown that the eligible uniform eigenstrains inducing uniform eigenstresses can be shear or non-shear eigenstrains, depending on the geometric shape and the material constants of the inclusion. The study indicates that inclusions of degenerated form, like spheroids, spheres and circles, may also maintain uniform eigenstresses. At last, the corresponding discussion for the inhomogeneous sliding inclusion problem with both eigenstrain and remote loading is also given.     

弱界面异质球形夹杂本征应变问题的解析解

本文研究了具有线弹簧弱界面的异质球形夹杂的本征应变问题，所采用的线弹簧界面模型既能界面的切线方向滑动，又能考虑界面的法线方向张开，根据叠加原理、原问题的弹性场可分成三部分；二部分由真实均匀本征应变所引起，另一部分由等效的非均匀本征应变所引起，后一部分则由虚拟的Ｓｏｍｉｇｌｉａｎａ位错场所产生。本文求得了等效非均匀本征应变和虚拟位错场的Ｂｕｒｇｅｒ矢量的解析表达式，进而确定的问题的弹性场。     

The effective thermoelectroelastic properties of microinhomogeneous materials

The paper is concerned with composite materials which consist of a homogeneous matrix phase with a set of inclusions uniformly distributed in the matrix. The components of these materials are considered to be ideally elastic and exhibit piezoelectric properties. One of the variants of the self-consistent scheme, the Effective Field Method (EFM) is applied to calculate effective dielectric, piezoelectric and thermoelastic properties of such materials, taking into account the coupled electroelastic effects. At first the coupled thermoelectroelastic problem for a homogeneous medium with an isolated inclusion is solved. For an ellipsoidal inclusion and constant external field the solution of this problem is found in a closed analytic form. This solution is then used in the EFM to derive the effective thermoelectroelastic operator for the composite containing a random array of ellipsoidal inclusions. Explicit formulae for the electrothermoelastic constants are given for composites, reinforced by spheroidal inclusions.     

Effect of interface stresses on the elastic deformation of an elastic half-plane containing an elastic inclusion

The effect of the interface stresses is studied upon the size-dependent elastic deformation of an elastic half-plane having a cylindrical inclusion with distinct elastic properties. The elastic half-plane is subjected to either a uniaxial loading at infinity or a uniform non-shear eigenstrain in the inclusion. The straight edge of the half-plane is either traction-free, or rigid-slip, or motionless, which represents three practical situations of mechanical structures. Using two-dimensional Papkovich–Neuber potentials and the theory of surface/interface elasticity, the elastic field in the elastic half-plane is obtained. Comparable with classical result, the new formulation renders the significant effect of the interface stresses on the stress distribution in the half-plane when the radius of the inclusion is reduced to the nanometer scale. Numerical results show that the intensity of the influence depends on the surface/interface moduli, the stiffness ratio of the inclusion to the surrounding material, the boundary condition on the edge of the half-plane and the proximity of the inclusion to the edge.     

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.     

基于等效特征应变原理的复合材料有效模量细观力学分析

基于等效特征应变原理,提出了一种新的复合材料有效模量细观力学分析方法。首先,在等效特征应变原理基础上提出平均等效特征应变原理,它可用于解决有限体下任意形状(无论是凸或凹形)的单个夹杂或多个夹杂的弹性变形问题。其次,将平均等效特征应变原理与细观力学直接均匀法相结合,来分析确定复合材料的有效模量。最后利用复合材料纤维与基体的力学性能参数及纤维的体分比,借助MATLAB编程方法,预测其有效模量。通过将理论预测值与已有的的试验值、其它理论预测值进行对比,验证了新分析方法的合理性和分析精度。     

Thickness effect of a thin film on the stress field due to the eigenstrain of an ellipsoidal inclusion

Solutions of the stress field due to the eigenstrain of an ellipsoidal inclusion in the film/substrate half-space are obtained via the Fourier transforms and Stroh eigenrelation equations. Based on the acquired solutions, the effect of a thin film’s thickness on the stress field is investigated with two types of ellipsoidal inclusions considered. The results in this paper show that if the thickness of the thin film increases, its effect on the stress field will become weaker, and can even be neglected. In the end, a guide rule is introduced to simplify the calculation of similar problems in engineering.     

Some three-dimensional problems for an elastic medium with isolated rigid inclusions

Problems of determining stresses in isolated ellipsoidal rigid inclusions contained in an isotropic elastic space subjected to uniformly distributed loading at infinity are considered. A solution in a closed form is constructed for inclusions shaped as ellipsoids of revolution.     

An investigation of elastic field created by randomly distributed inclusions

In this paper, the elastic field created by randomly distributed inclusions is studied. The inclusions are considered to be randomly distributed in the material, and have random orientation and size. The random point field model is proposed to describe the randomness of inclusion position, orientation and size. As a special case, when phase transformation inclusions are uniformly distributed in the material, and have non-random orientation, the theory gives the same result as Mori and Tanaka (1973. Acta Metallurgica 21, 571). The elastic field created by randomly distributed dislocation loops is also considered in some detail, and it is found that the continuum theory of dislocation loops is applicable only when the size of the dislocation loop becomes infinitesimal.     

Bounds on the overall properties of composites with ellipsoidal inclusions

A new model is put forward to bound the effective elastic moduli of composites with ellipsoidal inclusions. In the present paper, transition layer for each ellipsoidal inclusion is introduced to make the trial displacement field for the upper bound and the trial stress field for the lower bound satisfy the continuous interface conditions which are absolutely necessary for the application of variational principles. According to the principles of minimum potential energy and minimum complementary energy, the upper and lower bounds on the effective elastic moduli of composites with ellipsoidal inclusions are rigorously derived. The effects of the distribution and geometric parameters of ellipsoidal inclusions on the bounds of the effective elastic moduli are analyzed in details. The present upper and lower bounds are still finite when the bulk and shear moduli of ellipsoidal inclusions tend to infinity and zero, respectively. It should be mentioned that the present method is simple and needs not calculate the complex integrals of multi-point correlation functions. Meanwhile, the present paper provides an entirely different way to bound the effective elastic moduli of composites with ellipsoidal inclusions, which can be developed to obtain a series of bounds by taking different trial displacement and stress fields.     

Harmonic elastic inclusions in the presence of point moment

We employ conformal mapping techniques to design harmonic elastic inclusions when the surrounding matrix is simultaneously subjected to remote uniform stresses and a point moment located at an arbitrary position in the matrix. Our analysis indicates that the uniform and hydrostatic stress field inside the inclusion as well as the constant hoop stress along the entire inclusion–matrix interface (on the matrix side) are independent of the action of the point moment. In contrast, the non-elliptical shape of the harmonic inclusion depends on both the remote uniform stresses and the point moment.     

A circular inclusion with a nonuniform interphase layer in anti-plane shear

An analytical solution is derived for the problem of a nonuniformly coated circular inclusion in an unbounded matrix under anti-plane deformations. The inclusion/interphase/matrix system is subject to (1) remote uniform shear and uniform eigenstrain imposed on the circular inclusion, and (2) a screw dislocation or a point force in the matrix. It is found that the varying interphase thickness will exert a significant influence on the nonuniform stress field within the circular inclusion, and on the direction and magnitude of the image force acting on a screw dislocation. In the course of development, it is found that the presence of certain coated inclusions, which are termed stealth, will not cause change of elastic energy in the body. The derived analytical solution for a screw dislocation is then employed as Green’s function to investigate a radial matrix crack interacting with the nonuniformly coated inclusion. The numerical results show that the varying interphase thickness will also affect the stress intensity factors.     

The stress field of a sliding inclusion

This paper discusses the stress fields when a spheroidal inclusion, free to slip along its interface, is subjected to a constant nonshear eigenstrain, and when an elastic body containing the inhomogeneity is under all-around tension or uniaxial tension at infinity. In each case the stress field in the inclusion or the inhomogeneity is not constant, contrary to Eshelby's solution. When sliding takes place, the stress increases locally compared with the perfect bonding case, but the elastic energy decreases due to the relaxation. The relative displacement (slip) along the interface is also evaluated.