THE REMARKABLE NATURE OF RADIALLY SYMMETRIC DEFORMATION OF ANISOTROPIC PIEZOELECTRIC INCLUSION

The present paper deals with spherically symmetric deformation of an inclusion- matrix problem, which consists of an infinite isotropic matrix and a spherically uniform anisotropic piezoelectric inclusion. The interface between the two phases is supposed to be perfect and the system is subjected to uniform loadings at infinity. Exact solutions are obtained for solid spherical piezoelectric inclusion and isotropic matrix. When the system is subjected to a remote traction, analytical results show that remarkable nature exists in the spherical inclusion. It is demonstrated that an infinite stress appears at the center of the inclusion. Furthermore, a cavitation may occur at the center of the inclusion when the system is subjected to uniform tension, while a black hole may be formed at the center of the inclusion when the applied traction is uniform pressure. The appearance of different remarkable nature depends only on one non-dimensional material parameter and the type of the remote traction, while is independent of the magnitude of the traction.     

PLANE STRAIN PROBLEM OF PIEZOELECTRIC SOLID WITH ELLIPTIC INCLUSION

By using the complex variables function theory, a plane strain electro-elastic analysis was performed on a transversely isotropic piezoelectric material containing an elliptic elastic inclusion, which is subjected to a uniform stress field and a uniform electric displacement loads at infinity. Based on the present finite element results and some related theoretical solutions, an acceptable conjecture was found that the stress field is constant inside the elastic inclusion. The stress field solutions in the piezoelectric matrix and the elastic inclusion were obtained in the form of complex potentials based on the impermeable electric boundary conditions.     

Stress and electric fields in an infinite electrostrictive plate with an elliptic inhomogeneity

The stress and electric fields in electrostrictive materials under general electric loading at infinity are obtained in this paper. It is shown that the pseudo total stresses are continuous in the whole body. The elliptic inhomogeneity problem is first discussed in this paper and its solution is also given. The results show that the stress in the inhomogeneity is not uniform which is different from the solution of Eshelby theory for elastic materials. When the inhomogeneity and matrix have the same dielectric permittivity or the matrix is a non-electrostrictive material, the stress field is uniform in the inhomogeneity. The form of stress function is simple when the inhomogeneity degenerates to a circle.     

On the anti-plane shear of an elliptic nano inhomogeneity

The elastic field of an elliptic nano inhomogeneity embedded in an infinite matrix under anti-plane shear is studied with the complex variable method. The interface stress effects of the nano inhomogeneity are accounted for with the Gurtin–Murdoch model. The conformal mapping method is then applied to solve the formulated boundary value problem. The obtained numerical results are compared with the existing closed form solutions for a circular nano inhomogeneity and a traditional elliptic inhomogeneity under anti-plane. It shows that the proposed semi-analytic method is effective and accurate. The stress fields inside the inhomogeneity and matrix are then systematically studied for different interfacial and geometrical parameters. It is found that the stress field inside the elliptic nano inhomogeneity is no longer uniform due to the interface effects. The shear stress distributions inside the inhomogeneity and matrix are size dependent when the size of the inhomogeneity is on the order of nanometers. The numerical results also show that the interface effects are highly influenced by the local curvature of the interface. The elastic field around an elliptic nano hole is also investigated in this paper. It is found that the traction free boundary condition breaks down at the elliptic nano hole surface. As the aspect ratio of the elliptic hole increases, it can be seen as a Mode-III blunt crack. Even for long blunt cracks, the surface effects can still be significant around the blunt crack tip. Finally, the equivalence between the uniform eigenstrain inside the inhomogeneity and the remote loading is discussed.     

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.     

A cracked sliding interface between anisotropic bimaterials

A general solution for the stresses and displacements of a cracked sliding interface between anisotropic bimaterials subjected to uniform tensile stress at infinity is given by using the Stroh’s formulation. Horizontal and vertical opening displacements on the interface, stress intensity factors, and energy release rate are expressed in real form, which are valid for any kind of anisotropic materials including the degenerate materials such as isotropic materials. It is observed that stresses exhibit the traditional inverse square root singularities near the crack tips, and the vertical opening displacement and energy release rate are intimately related to a real parameter λ determined by the elastic constants of the anisotropic bimaterials.     

应用边界积分法求圆形夹杂问题的解析解

边界元方法作为一种数值方法, 在各种科学工程问题中得到了广泛的应用.本文参考了边界元法的求解思路, 从Somigliana等式出发, 利用格林函数性质,得到了一种边界积分法, 使之可以用来寻求弹性问题的解析解.此边界积分法也可以从Betti互易定理得到. 应用此新方法, 求解了圆形夹杂问题.首先设定夹杂与基体之间完美连接, 将界面处的位移与应力按照傅里叶级数展开,根据问题的对称性与三角函数的正交性来简化假设, 减少待定系数的个数.其次选择合适的试函数(试函数满足位移单值条件以及无体力的线弹性力学问题的控制方程),应用边界积分法, 求得界面处的位移与应力的值. 然后再求解域内位移与应力.得到了问题的精确解析解, 当夹杂弹性模量为零或趋向于无穷大时,退化为圆孔或刚性夹杂问题的解析解. 求解过程表明,若问题的求解区域包含无穷远处时, 所取的试函数应满足无穷远处的边界条件.若求解区域包含坐标原点, 试函数在原点处位移与应力应是有限的.结果表明了此方法的有效性.      

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.     

Stress field of a coated arbitrary shape inclusion

This paper presents an effective method for the plane problem of a coated inclusion of arbitrary shape embedded in an isotropic matrix subjected to uniform stresses at infinity. Based on the complex variable method combined with the expansion of Faber series and Laurent series, the complex potentials in the matrix, the coating and the arbitrary shape inclusion are given in the form of series with unknown coefficients. The stress and displacement continuous conditions on the interfaces are then used to produce a set of linear equations containing all the coefficients. Through solving these linear equations, the complex potentials are finally obtained in the three phases. Additionally, numerical results are presented and graphically shown to investigate the influence of inclusion geometry and coating on the stress distribution along the interfaces for the cases of a coated elliptic, square and triangle inclusions, respectively. It is found that the coating has little effects on the interface stress for a hard inclusion, while it impacts greatly for a soft inclusion. Especially, it is also found that the stresses show the nature of intense fluctuations near the corner of the triangle inclusion, since the inclusion in this case is similar to a wedge.     

General treatment of mode III interfacial crack problems in piezoelectric materials

Summary  The anti-plane problem of N collinear interfacial cracks between dissimilar transversely isotropic piezoelectric media, which are subjected to piecewise uniform out-of-plane mechanical loading combined with in-plane electric loading at infinity, and also a line loading at an arbitrary point, is addressed by using the complex function method. In comparison with other relevant works, the present study has two features: one is that the analysis is based on the permeable crack model, i.e. the cracks are considered as permeable thin slits, and, thus, both the normal component of electric displacement and the tangential component of electric field are assumed to be continuous across these slits. The other feature is that explicit closed-form solutions are given not only in piezoelectric media, but also inside cracks when the media are subjected to the most general loading. It is shown that the singularities of electric displacement and electric field in the media are always dependent on that of stress for the general case of loading, and all the singularities of field variables are independent of the applied uniform electric loads at infinity. For the interfacial cracks the electric field is square-root singular at the crack tips and shows jumps across the interface, while the normal component of the electric field is linearly variable inside the crack, but the tangential component is square-root singular. However, for a homogeneous medium with collinear cracks, the electric field is always nonsingular in the medium while the electric displacement exhibits square-root singularity. Moreover, in this case, the electric field inside any crack is equal to a constant when uniform loads are applied at infinity. Received 22 November 1999; accepted for publication 20 July 2000     

Closed form solution and numerical analysis for Eshelby’s elliptic inclusion in plane elasticity

This paper presents a closed form solution and numerical analysis for Es- helby＇s elliptic inclusion in an infinite plate. The complex variable method and the confor- real mapping technique are used. The continuity conditions for the traction and displace- ment along the interface in the physical plane are reduced to the similar conditions along the unit circle of the mapping plane. The properties of the complex potentials defined in the finite elliptic region are analyzed. From the continuity conditions, one can separate and obtain the relevant complex potentials defined in the inclusion and the matrix. From the obtained complex potentials, the dependence of the real strains and stresses in the inclusion from the assumed eigenstrains is evaluated. In addition, the stress distribution on the interface along the matrix side is evaluated. The results are obtained in the paper for the first time.     

Formulation and implementation of a singular anisotropic finite element

The formulation and implementation of a singular finite element for analyzing homogeneous anistropic materials is presented in this paper. Lekhnitskii's stress function method is used to formulate the boundary value problem with the stress function expressed as a Laurent series. The development of the element stiffness matrix and the method of integrating the element to conventional displacement based finite element programs is shown. The stiffness matrix generation is based on a least squates collocation technique to satisfy displacement continuity boundary conditions at the element interface. Implementation of the element is demonstrated for cracked anisotropic materials subjected to inplane loading. Center cracked, on and off-axis coupons under tensile loading are analyzed using the element. It is shown that the stress distributions and intensity factors compare well with those obtained using other methods.     

面向剪切下各向异性三相椭圆夹杂中均匀应力

论证了只要合适选择中间界面层的弹性常数，各向异性线弹性固体在远场均匀反平面剪切应力作用下三相椭圆夹杂内椭圆上仍存在均匀应力场。讨论了内外两椭圆除过其中心相同外无其它任何几何限制条件。所给出的数值算例显示出该结论的正确性。该方法为纤维增强复合材料的设计提供了一条新途径。     

压电体椭圆孔边的力学分析

基于复变函数的方法,以PZT-4材料为例,分别采用精确电边界条件和非导通电边界条件进行了远场均匀载荷作用下的横观各向同性压电体椭圆孔的力学分析并与相关结果进行对比。结果表明当椭圆孔退化为圆孔时,无论在远场作用力载荷或电载荷,两种电边界条件下的结果均能完全吻合。随着椭圆孔的愈加尖锐化,非导通电边界条件逐渐不能适用。     

2维各向异性弹性力学中的Stroh公式和某些不变量

本文简要综述了各向异性弹性体2维变形的Stroh公式。该公式数学上优美而应用上有力。它的最近进展使我们将已有的复数形式的解转变成实数形式,进而可去解决一些迄今尚未解出的问题。结果,发现了许多有趣的而在物理上难以预料的现象。例如,对于界面包含x_3轴的各向异性双介质体,如果线力和线位错都沿x_3轴,那末包含x_3轴的任意径向平面上的面拉力将不随径向平面的选择而变化。而且,在任何固定点上的应力与界面的方位无关。要揭示的其他一些不变性质是:双介质材料中的虚拟力,界面裂纹上的应力奇异,界面裂纹表面上的振荡,等等;另外,在具有一个有限裂纹,或一个椭圆孔,或一个刚性椭圆核的无穷大均匀介质中的某些应力和位移也是不变量。这些不变量的大部分尚未有它们的物理解释。在实数形式的解中,Barnett-Lothe张量常常出现。在绕x_3轴转动时,这些张量的分量的某些组合是不变的。这些组合将在本文中推导并列出。最后我们给出了某些各向异性弹性材料,它们的Barnett-Lothe张量是不变量,即这些张量的所有分量在绕x_3轴转动时是不变的。     

考虑晶界滑错的金属多晶体响应的自洽方法

本文首先建立含有三种介质（各向异性基体、各向异性夹杂，界面层）的平面应变夹杂模型，将基体和夹杂位移场展开为多项式级数，假设界面层很薄，运用变分原理得出这一问题的近似解。将上述夹杂问题的解和ＨＩＬＬ自洽方法相结合，给出了考虑晶界滑错效应的金属多晶体弹塑性响应。     

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.     

压电复合材料中的Eshelby夹杂问题

通过采用解析延拓和共形映射技术，获得了压电复合材料中有关Eshelby夹杂几个典型问题的精确弹性解答，即横观各向同性压电介质中任意形状的Eshelby夹杂与圆柱异相夹杂间相互作用；一般各向异性压电介质中任意形状的Eshelby夹杂与双压电材料所形成界面的相互作用．成功求解这些问题的关健在于构造一个辅助函数．与Ru所采用的方法不同，所引入的辅助函数在无穷远点不存在极点，从而使得所展开的分析更加自然合理．分析结果清楚地揭示出Eshelby夹杂的存在对压电复合材料机电耦合响应将产生不容被忽视的影响．很典型的一个例于是当一个Eshelby椭圆夹杂与圆柱异相夹杂相互作用时，每个夹杂体内部的应力场和电场都将是不均匀的；另一个例于是位于界面附近的Eshelby夹杂有可能是界面发生损伤的一个重要原因．     

Contact zone assessment for a fast growing interface crack in an anisotropic bimaterial

Summary A plane strain problem for a crack with a frictionless contact zone at the leading crack tip expanding stationary along the interface of two anisotropic half-spaces with a subsonic speed under the action of various loadings is considered. The cases of finite and infinite-length interface cracks under the action of a moving concentrated loading at its faces are considered. A finite-length crack for a uniform mixed-mode loading at infinity is considered as well. The associated combined Dirichlet-Riemann boundary value problems are formulated and solved exactly for all above-mentioned cases. The expressions for stresses and the derivatives of the displacement jumps at the interface are presented in a closed analytical form for an arbitrary contact zone length. Transcendental equations are obtained for the determination of the real contact zone length, and the associated closed form asymptotic formulas are found for small values of this parameter. It is found that independently of the types of the crack and loading, an increase of the crack tip speed leads to an increase of the real contact zone length and the correspondent stress intensity factor. The latter increase significantly for an interface crack tip speed approaching the Ragleigh wave speed.     

Elastic field of an isotropic matrix with a nanoscale elliptical inhomogeneity

In traditional continuum mechanics, the effect of surface energy is ignored as it is small compared to the bulk energy. For nanoscale materials and structures, however, the surface effects become significant due to the high surface/volume ratio. In this paper, two-dimensional elastic field of a nanoscale elliptical inhomogeneity embedded in an infinite matrix under arbitrary remote loading and a uniform eigenstrain in the inhomogeneity is investigated. The Gurtin–Murdoch surface/interface elasticity model is applied to take into account the surface/interface stress effects. By using the complex variable technique of Muskhelishvili, the analytic potential functions are obtained in the form of an infinite series. Selected numerical results are presented to study the size-dependency of the elastic field and the effects of surface elastic moduli and residual surface stress. It is found that the elastic field of an elliptic inhomogeneity under uniform eigenstrain is no longer uniform when the interfacial stress effects are taken into account.