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

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

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

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

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

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

横观各向同性介质中椭圆夹杂受非弹性剪切变形引起的弹性场的确定

本文求解了横观各向同性介质中椭圆夹杂内受非弹性剪切变形引起的弹性场。采用各向异性弹性力学平面问题的复变函数解法,结合保角变换,获得夹杂内应变能和基体内边界的应力分布和相应的应变能的表达式。进一步,根据最小应变能原理,获得表征夹杂平衡边界的两个特征剪切应变,从而得到了弹性场的解析解。通过应力转换关系,验证了应力解满足夹杂边界上法向正应力和剪应力的连续条件,表明了该解的正确性。本文解可用于复合材料断裂强度的分析中。     

磁电弹性复合材料热接触特性研究

磁电弹性(Magneto-electro-elastic, MEE)复合材料在接触热载荷作用下将产生复杂热应力、热电和热磁多物理场响应.基于半解析法建立MEE复合材料热接触模型,其中,推导了单位集中法向力、切向力、电荷、磁荷和温升载荷下各物理量的频率响应函数,引入热流-温升影响系数计算摩擦热产生的温升,并采用离散卷积快速傅里叶变换和共轭梯度法加速其计算过程.将模型计算结果与有限元仿真进行对比,验证模型有效性.进一步利用所提模型分析热接触过程中摩擦热对各物理场的影响规律,结果表明：滑动速度、摩擦系数和表面形貌改变将影响摩擦热分布,进而显著影响应力、电势和磁势的大小和分布;弹性场和电场与摩擦温升呈负相关,而磁场与其为正相关,并且耦合场对温升的敏感程度由高至低为弹性场、电场和磁场.     

各向异性介质中的椭圆夹杂问题

本文研究了含椭圆夹杂的弹性体在多项式荷载作用下的二维变形问题，获得了介质和夹杂中的弹性场，证明了夹杂内部的应力变场以荷载的同阶多项式形式出现，而介质中的弹性场也能用椭圆坐标ζ＾－１２的多项式形式表示出来，并在此基础上，以受剪力作用的含夹杂或空孔的悬臂梁为例，求解了梁中的应力扰动现象，并获得了夹杂或空孔周转的环向应力。     

电致伸缩圆柱体电场冲击响应研究

根据弹性力学轴对称平面应变问题的基本方程,采用有限Hankel变换及其逆变换辅以Laplace变换技术,得到了轴对称径向突加电场载荷条件下电致伸缩材料实心圆柱体的动态位移和应力响应的解析解．由于电冲击引起圆柱体内弹性波的传播,动态位移和应力随时间呈不同峰值的周期性变化．数值计算表明,随着半径增大,位移的响应相应增加,在圆柱表面达到静态位移数值的5倍以上;在圆柱表面附近,动态应力响应呈周期性拉压变化,最大幅度可达到静态应力的20倍左右．因此,在计算位移和应力场时,必须考虑电场冲击因素．     

基于数字散斑相关实验测量的材料构型力的计算方法

材料构型力学主要研究材料中的缺陷(夹杂、空穴、位错、裂纹、塑性区等)的构型(形状、尺寸和位置)改变时,所引起的系统自由能的变化。本研究将基于数字散斑相关技术,实验测量材料试件的位移场分布,随后通过材料构型力的定义式,计算求得弹塑性材料中缺陷构型力的分布。其方法概括如下:位移场通过数字图像相关技术测得;应变及位移梯度场利用三次样条拟合获得;线弹性材料应力通过简单线弹性本构方程获取,而塑性材料的表面应力场通过Ramberg-Osgood本构方程计算求得;弹塑性应变能密度分布则由应力-应变曲线数值积分获得。该方法对普通弹性材料或者弹塑性材料均适用,可以用于各种不同的缺陷及缺陷群的材料构型力测量。     

双周期圆柱形压电夹杂反平面问题的精确解及其应用

研究无限压电介质中双周期圆柱形压电夹杂的反平面问题.借鉴Eshelby等效夹杂原理,通过引入双周期非均匀本征应变和本征电场,构造了一个与原问题等价的均匀介质双周期本征应变和本征电场问题.利用双准周期Riemann边值问题理论,获得了夹杂内外严格的电弹性解.作为压电纤维复合材料的一个重要模型,预测了压电纤维复合材料的有效电弹性模量.     

轴承钢滚动接触疲劳亚表面夹杂处损伤分析

轴承钢在滚动接触疲劳(RCF)中失效的主要原因之一是亚表面白蚀区(WEA)的形成.本文中从塑性应变累积引起剪切局域化新的角度对WEA进行了研究.通过耦合晶体塑性和相场损伤理论建立了损伤演化本构模型,研究了非金属夹杂处塑性应变累积和损伤演化.研究表明,接触疲劳载荷引起的塑性应变局域化导致了剪切带的形成.剪切带的形貌、取向和应变与WEA的一致,表明WEA实际上是应变局域化的剪切带.晶体取向对WEA损伤的形成和发展有很大的影响,WEA仅在择优的晶体取向下形成.与软夹杂周围的剪切带和损伤演化不同,硬夹杂处的剪切带与夹杂相切,形成的4条剪切变形带将夹杂“包围”.剪切带内部处于高应变和低应力的状态,带中心处应变达到最大,随带宽两侧急剧减小,而中心处应力却最小,几乎为零,沿带宽两侧增大,这说明裂纹在剪切带内萌生和扩展.该结论阐明了裂纹和WEA形成的关系,即裂纹是在WEA形成过程中因应变不协调产生,而非裂纹先产生,裂纹上下表面相互摩擦导致WEA形成.     

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.     

Inclusion of arbitrary polygon with graded eigenstrain in an anisotropic piezoelectric full plane

In this paper, an exact closed-form solution for the Eshelby problem of a polygonal inclusion with a graded eigenstrain in an anisotropic piezoelectric full plane is presented. For this electromechanical coupling problem, by virtue of Green’s function solutions, the induced elastic and piezoelectric fields are first expressed in terms of line integrals on the boundary of the inclusion. Using the line-source Green’s function, the line integral is then carried out analytically for the linear eigenstrain case, with the final expression involving only elementary functions. Finally, the solution is applied to the semiconductor quantum wire (QWR) of square, triangle, circle and ellipse shapes within the GaAs (0 0 1) substrate. It is demonstrated that there exists significant difference between the induced field by the uniform eigenstrain and that by the linear eigenstrain. Since the misfit eigenstrain in most QWR structures is actually non-uniform, the present solution should be particularly appealing to nanoscale QWR structure analysis where strain and electric fields are coupled and are affected by the non-uniform misfit strain.     

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

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

The determination of the elastic field of a polygonal star shaped inclusion

Recently we found that the elastic field is uniform in a pentagonal star (five-pointed star inclusion) [1], and in a triangular inclusion [2], when an eigenstrain is distributed uniformly in these inclusions. This result is similar to the famous result of Eshelby (1957) that the elastic field is uniform in an ellipsoidal inclusion in an infinitely body when an eigenstrain is distributed uniformly in the ellipsoidal inclusion. We also found that for a Jewish star (Star of David or six points star) or a rectangular inclusion subjected to a uniform eigenstrain, the stress field is not uniform in these inclusions. These results also hold for two dimensional plane strain cases. Furthermore these analytical results are confirmed experimentally by photoelasticity method. In this paper, we investigate a more general inclusion of an m-pointed polygonal inclusion subjected to the uniform eigenstrain. We conclude that the stress field is uniform when m is odd number. This conclusion agrees with the speculation made by B. Boley after the author's talk at Shizuoka [2].     

Size dependent,non-uniform elastic field inside a nano-scale spherical inclusion due to interface stress

The primary objective of the present paper is to analyze the influence of interface stress on the elastic field within a nano-scale inclusion. Special attention is focused on the case of non-hydrostatic eigenstrain. From the viewpoint of practicality, it is assumed that the inclusion is spherically shaped and embedded into an infinite solid, within which an axisymmetric eigenstrain is prescribed. Following Goodier’s work, the elastic fields inside and outside the inclusion are obtained analytically. It is found that the presence of interface stress leads to conclusion that the elastic field in the inclusion is not only dependent on inclusion size but also on non-uniformity. The result is in strong contrast to Eshelby’s solution based on classical elasticity, and it is helpful in the understanding of relevant physical phenomena in nano-structured solids.     

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.     

平面任意形状等剪切模量异质夹杂问题的解析解

本文研究任意形状夹杂域在受到远端均匀荷载和均匀本征应变作用下的弹性场问题,其中基体和夹杂的材料不同但具有相同的剪切模量。利用等效理论将远端均匀荷载引起的扰动转化为等效均匀本征应变的作用,再利用K-M势函数表达扰动场问题的界面连续条件;借助于黎曼映射定理,用洛朗多项式将平面光滑闭合曲线外部区域映射到单位圆外部区域,借助柯西积分公式和Faber多项式求解了等剪切本征应变下夹杂和基体的K-M势函数的显式解析解,其中考虑了夹杂相对于基体的刚体位移。将得到的结果与相关文献的结果进行对比,表明了本论文的方法和结果是有效的和正确的。     

Piezoelectric composites with periodic multi-coated inhomogeneities

A new, robust homogenization scheme for determination of the effective properties of a periodic piezoelectric composite with general multi-coated inhomogeneities is developed. In this scheme the coating does not have to be thin, the shape and orientation of the inclusion and coatings do not have to be identical, their centers do not have to coincide, their properties do not have to remain uniform, and the microstructure can be with the 2D elliptic or the 3D ellipsoidal inclusions. The development starts from the local electromechanical equivalent inclusion principle through the introduction of the position-dependent equivalent eigenstrain and electric field. Then with a Fourier series expansion and a superposition procedure, the volume-averaged equivalent eigenstrain and electric field for each phase are obtained. The results in turn are used in an energy equivalent criterion to determine the effective properties of the composite. In this model the interphase interactions in each multi-coated particle and the long-range interactions between the periodically distributed particles are fully accounted for. To demonstrate its wide range of applicability, we applied it to examine the properties of several periodic composites: (i) piezoelectric PZT spherical particles in a polymer matrix, (ii) continuous glassy fibers with thin PZT coating in an epoxy matrix, (iii) spherical PZT particles coated by thick or functionally graded piezoelectric layer, (iv) spheroidal voids coated with a thick non-piezoelectric layer in a PZT matrix, and (v) spherical piezoelectric inhomogeneities with eccentric, non-uniform thickness coating. The calculated results reflect the complex nature of interplay between the properties of core, matrix, and coating, as well as whether the coating is uniform, functionally graded, or eccentric. The accuracy of this new scheme is checked against the double-inclusion and other micromechanics models, and good agreement is observed.