含刚性线夹杂及裂纹的各向异性压电材料耦合场分析

采用各向异性弹性力学中Ｓｔｒｏｈ方法对含刚性线夹杂及裂纹的无限大各向异性压电材料耦合的弹性场和电场进行了分析。并得到夹杂和基体界面间耦合场的实型显函表达式及夹杂尖端的１／２阶奇异性。     

含界面层圆形夹杂中螺旋位错的弹性分析

研究了复合材料含界面层圆形夹杂内部的一个螺旋位错在夹杂，界面层与基体材料中产生的弹性干涉，将复变函数的分区亚纯函数理论与柯西型积分，罗朗级数相结合，求出了各分区复势的解析关系，化为一个关于界面层复势的函数方程，用显示式表达了问题的结果，揭示了界面层参数对位错干涉能与位错力的影响规律。该解析方法较经典级数方法未知量大量减少，表达式更加简洁，结果的特殊情形包含了若干已有成果。     

含椭圆夹杂正交各向异性体的界面应力研究

本文研究了远场作用反平面载荷时含椭圆夹杂正交各向异性体的界面应力分布规律．利用解析函数边值问题理论和共形映射技术,推导了反平面载荷下含椭圆夹杂正交异性体的精确解,获得了夹杂和基体内应力场的闭合解,并通过有限元结果验证了本文解析解的有效性．研究表明：基体材料主方向弹性模量比C55/C44和夹杂形状比 对界面应力影响显著;基体材料主方向模量比C55/C44对界面应力的影响受夹杂/基体模量比Cf/C44的限制．     

滑动界面球形夹杂对平面压缩波的散射

非理想粘结界面对多相材料力学性能的影响日益受到重视。本文研究了无限各向同性基体中的滑介面球形单夹杂对平面压缩的散射问题。夹杂与基体间的界面为非理想粘结界面,在剪应力的作用下将出现界面两侧相对滑移。假定界面相对滑动位移与界面剪应力成正比,在这种线弹簧型滑动界面条件下,通过波函数的级数展开法,获得了夹杂在基体中反射波和折射波以及应力场的解析表达式,并讨论了界面自由滑动和刚性夹杂等特例。     

无限介质中单个纳米涂层夹杂反平面问题的两个模型

采用零厚度界面模型和界面层模型研究了无限介质中单个纳米涂层圆柱形夹杂的反平面问题,利用复变函数方法获得了两种模型夹杂、涂层和基体内应力场的封闭解析解．研究表明,当界面层模型中的界面相厚度趋于零时,界面层模型可以解析地退化为零厚度界面模型．数值算例分析了界面模量不同取值时应力场的分布和应力的尺度依赖性．本文结果丰富了对纳米夹杂力学行为的认识,并可为直接采用零厚度界面模型有困难的纳米夹杂问题的研究提供有价值的参考．     

弹性滑错界面椭圆夹杂应力场

本文采用弹性弹簧来模拟夹杂和基体之间界面存在的缺陷．应用Ｐａｐｋｏｖｉｃｈ Ｎｅｕｂｅｒ位移函数法给出了带有缺陷界面无穷大基体中含椭圆形夹杂问题的级数解，计算结果显示了滑错夹杂应力场的一些特点．     

功能梯度界面相对周期分布纤维增强复合材料反平面剪切的影响

基于复变函数理论和边界配点法,探索了功能梯度界面相在周期均匀分布纤维增强复合材料反平面剪切问题中所起的作用.由于纤维在复合材料基体中的周期分布是均匀的,将其简化成含一功能梯度界面相夹杂的方形单胞.采用分层均匀化方法,将功能梯度界面相离散成K层界面层.当K足够大时,每个界面层可视为匀质材料,同时计算得到的复合材料宏观性能...     

动载下复合材料基体裂纹与夹杂相互作用的数值模拟

采用ABAQUS软件及粘聚裂纹模型对动态拉伸载荷下复合材料中垂直于基体-夹杂界面的基体裂纹与夹杂的相互作用进行了数值模拟.结果表明:在给定的界面强度下,当加载率(或应变率)低于某一临界值(临界应变率)时裂纹将沿基体-夹杂界面扩展,当高于该临界值时裂纹可以穿过界面在夹杂中扩展.此外,随着界面强度的提高,临界应变率降低;当界面强度超过一定值后,裂纹扩展方向将不受外部载荷影响,裂纹将沿自相似方向扩展;当界面强度低于该临界值时,裂纹将沿界面扩展;并且临界应变率随夹杂尺寸的增大而降低,即小夹杂更难于破坏.上述结果可为前人的混凝土动态实验数据提供合理的细观理论解释.     

压电材料椭圆夹杂界面开裂问题的电弹性耦合解

本文研究了在反平面剪切和面内电场的共同作用下，压电材料椭圆夹杂的界面开裂问题，假定夹杂是刚性的导体，采用复变函数保角变换和级数展开方法，可确定压电材料基体的复势表达式，进而求得夹杂界面开裂的电弹性耦合的能量释放率。     

含界面效应纳米尺度圆环形涂层中螺型位错分析

研究了纳米尺度圆环形涂层(界面层)中螺型位错与圆形夹杂以及无限大基体材料的干涉效应.涂层与夹杂的界面和涂层与基体的界面均考虑界面应力效应.运用复势方法,获得了三个区域复势函数的解析解答.利用求得的应力场和Peach-Koehler公式,得到了作用在螺型位错上位错力的精确表达式.主要讨论了界面应力对涂层(界面层)中螺型位错运动和平衡稳定的影响规律.结果表明,界面应力对界面附近位错的运动有大的影响,由于界面应力的存在,可以改变涂层内位错与夹杂/基体干涉的引斥规律,并使位错在涂层内部产生三个稳定或非稳定的平衡点.考虑界面效应后,有一个额外的排斥力或吸引力作用在位错上,使原有的位错力增加或减小.     

Uniform Stress Inside an Anisotropic Elliptic Inclusion with Imperfect Interface Bonding

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. T.C.T. Ting is Professor Emeritus of University of Illinois at Chicago and Consulting Professor of Stanford University.     

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.     

非理想界面对三相复合材料应力场的影响

对非理想界面的三相复合材料,提出了计算弹性应力场的微观力学模型,在适当的简化假设下,对带界相的颗粒增强和纤维增强复合材料,得到了应力场的计算公式。以剪切载荷为例给出了数值例子。给出的数值结果表明非理想界面对三相复合材料应力场的影响。     

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.     

Recursive geometric integrators for wave propagation in a functionally graded multilayered elastic medium

The differential equations governing transfer and stiffness matrices and acoustic impedance for a functionally graded generally anisotropic magneto-electro-elastic medium have been obtained. It is shown that the transfer matrix satisfies a linear 1st order matrix differential equation, while the stiffness matrix satisfies a nonlinear Riccati equation. For a thin nonhomogeneous layer, approximate solutions with different levels of accuracy have been formulated in the form of a transfer matrix using a geometrical integration in the form of a Magnus expansion. This integration method preserves qualitative features of the exact solution of the differential equation, in particular energy conservation. The wave propagation solution for a thick layer or a multilayered structure of inhomogeneous layers is obtained recursively from the thin layer solutions. Since the transfer matrix solution becomes computationally unstable with increase of frequency or layer thickness, we reformulate the solution in the form of a stable stiffness-matrix solution which is obtained from the relation of the stiffness matrices to the transfer matrices. Using an efficient recursive algorithm, the stiffness matrices of the thin nonhomogeneous layer are combined to obtain the total stiffness matrix for an arbitrary functionally graded multilayered system. It is shown that the round-off error for the stiffness-matrix recursive algorithm is higher than that for the transfer matrices. To optimize the recursive procedure, a computationally stable hybrid method is proposed which first starts the recursive computation with the transfer matrices and then, as the thickness increases, transits to the stiffness matrix recursive algorithm. Numerical results show this solution to be stable and efficient. As an application example, we calculate the surface wave velocity dispersion for a functionally graded coating on a semispace.     

Analysis of a Partially Debonded Conducting Rigid Elliptical Inclusion in a Piezoelectric Matrix

ＩｎｔｒｏｄｕｃｔｉｏｎＤｕｅｔｏｔｈｅｉｒｉｎｔｒｉｎｓｉｃｅｌｅｃｔｒｏｍｅｃｈａｎｉｃａｌｃｏｕｐｌｉｎｇｐｒｏｐｅｒｔｉｅｓ,ｐｉｅｚｏｅｌｅｃｔｒｉｃｃｅｒａｍｉｃｓｈａｖｅｂｅｅｎｅｘｔｅｎｓｉｖｅｌｙｕｓｅｄｉｎｄｅｓｉｇｎｏｆｖａｒｉｏｕｓｅｌｅｃｔｒｏｎｉｃａｎｄｅｌｅｃｔｒｏｍｅｃｈａｎｉｃａｌｄｅｖｉｃｅｓｓｕｃｈａｓｓｅｎｓｏｒｓａｎｄａｃｔｕａｔｏｒｓ.Ｉｎｒｅｃｅｎｔｙｅａｒｓ,ｍｅｃｈａｎｉｃａｌａｎａｌｙｓｉｓｏｆｄｉｓｌｏｃａｔｉｏｎｓ ,ｃｒａｃｋｓ,ｃａｖｉｔｉｅ…     

Three-dimensional solution of smart laminated anisotropic circular cylindrical shells with imperfect bonding

This paper presents an exact solution for a simply-supported and laminated anisotropic cylindrical shell strip with imperfect bonding at the off-axis elastic layer interfaces and with attached anisotropic piezoelectric actuator and sensor subjected to transverse loading. In this research, the imperfect interface conditions are described in terms of linear relations between the interface tractions in the normal and tangential directions, and the respective discontinuities in displacements. The solution for an elastic (or piezoelectric) layer of the smart laminated cylindrical shell strip is obtained in terms of the six-dimensional (or eight-dimensional) pseudo-Stroh formalism, solution for multilayered system is then derived based on the transfer matrix method. Finally, a numerical example is presented to demonstrate the effect of imperfect interface on the static response of the smart laminated cylindrical shell. The derived solutions can serve as benchmark results to assess various approximate shell theories and numerical methods.     

A solution to the thermo-elastic interface crack branching in dissimilar anisotropic bi-material media

An interface crack or delamination may often branch out of the interface in a laminated composite due to thermal stresses developing around the delamination/crack tip when the media is exposed to heat flow induced by environmental events such as a sudden short-duration fire. In this paper, the thermo-elastic problem of interface crack branching in dissimilar anisotropic bi-media is studied by using the theory of Stroh’s dislocation formalism, extended to thermo-elasticity in matrix notation. Based on the complex variable method and the analytical continuation principle, the thermo-elastic interface crack/delamination problem is examined and a general solution in compact form is derived for dissimilar anisotropic bi-media. A set of Green’s functions is proposed for the dislocations (conventional dislocation and thermal dislocation/heat vortex) in anisotropic bi-media. These functions may be more suitable than those which have appeared in the literature on addressing thermo-elastic interface crack branching in dissimilar anisotropic bi-materials. Using the contour integral method, a closed form solution to the interaction between the dislocations and the interface crack is obtained. Within the scope of linear fracture mechanics, the thermo-elastic problem of interface crack branching is then solved by modelling the branched portion as a continuous distribution of dislocations. The influence of thermal loading and thermal properties on the branching behavior is examined, and criteria for predicting interface crack branching are suggested, based on the extensive numerical results from the study of various cases.     

General coupled solution of anisotropic piezoelectric materials with an elliptic inclusion

In this investigation, the Stroh formalism is used to develop a general solution for an infinite, anisotropic piezoelectric medium with an elliptic inclusion. The coupled elastic and electric fields both inside the inclusion and on the interface of the inclusion and matrix are given. The project supported by the National Natural Science Foundation of China     

层状横观各向同性地基上异性矩形薄板的弯曲解析解

选用更具广泛性的层状横观各向同性弹性地基模型,来分析四边自由各向异性矩形地基板的弯曲解析解。先基于直角坐标下横观各向同性体的静力胡海昌通解,借助双重傅里叶变换及矩阵传递法,获得层状横观各向同性地基的静力位移场和应力场;然后将异性薄板的弯曲控制方程,与基于层状横观各向同性弹性地基的位移解建立的板与地基变形协调方程相结合,先按对称性分解,再用三角级数法,得出层状横观各向同性弹性地基上四边自由各向异性矩形薄板的弯曲解析解,包括地基反力、板的挠度及板的内力的解析表达式。克服了数值法的弊端,取消了对地基反力的假设,且避免了矩阵指数函数的计算;同时考虑了地基的层状性及板和地基的各向异性,从而得到板的内力及地基反力更切实际的分布规律。算例结果与文献的有限元结果吻合良好,证明本文方法是切实可行的。