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

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

压电螺型位错和共线界面刚性线夹杂的干涉效应

研究了压电材料中压电螺型位错和共线界面导电刚性线夹杂的电弹干涉效应．运用复变函数解析延拓技术与奇性主部分析方法,获得了该问题的一般解答．作为算例,求出了界面含一条刚性线夹杂时两种压电介质区域广义应力函数的封闭形式解．导出了作用在位错上的像力和刚性线夹杂表面剪应力和电位移的解析表达式．讨论了界面刚性线长度,两种材料的剪切模量比和压电系数比对位错力和刚性线表面剪应力的影响规律．为进一步研究该类问题提供了一个基本解。     

一种颗粒复合材料物理性能计算模型与工程应用

结合傅里叶描述子法与格林函数法提出了任意形状夹杂物的Eshelby张量计算方法,基于分步夹杂法建立了一种具有任意形状颗粒的复合材料物理性能预测理论模型。通过傅里叶描述子法对扫描电镜(SEM)图像中颗粒形貌进行数学表征,构建由多个细观结构参数组合的几何模型;利用格林函数求解夹杂物Eshelby张量,考虑颗粒形态、大小和组分性能,研究细观结构对颗粒复合材料有效性能的影响。在此基础上,设计了一款颗粒复合材料物理性能预测软件。以铜铬合金为研究对象,分析重构的细观模型与预测的结果,并与有限元计算结果进行对比,验证了该计算模型的有效性和可靠性,以及在工程计算应用中的可行性。     

含界面刚性线夹杂与广义压电位错双压电介质的电弹性分析

研究了双压电材料中广义压电位错与分布于界面的刚性共线线夹杂相互作用问题.基于线性压电理论Stroh框架,相应的混合边值问题,可化为常见的Hilbert问题.求解Hilben问题,得到存在界面刚性线夹杂与位错时,压电体内所有场变量的显式表达.给出了由于界面和刚性线夹杂的存在,作用于广义压电位错上的广义Peach-Koehler镜像力.对均匀压电材料这一特殊情况,给出了数值算例,讨论了位错对刚性线夹杂端部场强度因子的干涉和它们之间的相互作用.结果可作为求解界面刚性线夹杂与微裂纹交互作用问题的Green函数,也可作为边界元方法的核函数.     

压电螺位错与含界面裂纹圆形涂层夹杂的干涉

研究位于基体或夹杂中任意点的压电螺型位错与含界面裂纹圆形涂层夹杂的电弹耦合干 涉问题. 运用复变函数方法，获得了基体，涂层和夹杂中复势函数的一般解答. 典型例 子给出了界面含有一条裂纹时，复势函数的精确级数形式解. 基于已获得的复势函数和广 义Peach-Koehler公式，计算了作用在位错上的像力. 讨论了裂纹几何条件，涂层厚度和材 料特性对位错平衡位置的影响规律. 结果表明，界面裂纹对涂层夹杂附近的位错运动有很大 的影响效应，含界面裂纹涂层夹杂对位错的捕获能力强于完整粘结情况；并发现界面裂纹长 度和涂层材料常数达到某一个临界值时可以改变像力的方向. 解答的特殊情形包含了以 往文献的几个结果.     

基于界面层模型的双周期纳米夹杂复合材料反平面问题研究

利用复变函数方法并结合双准周期Riemann边值问题理论，获得了含双周期分布非均匀相（夹杂/界面层）的复合材料在远场均匀反平面应力下弹性场的全场解答．该解答可用于对纳米夹杂复合材料的应力进行分析，结合平均场理论也用于预测纳米夹杂复合材料的有效性能．计算结果表明：当夹杂尺度在纳米量级时，应力和有效反平面剪切模量具有明显的尺度依赖性，并且随着夹杂尺寸的增加，趋近于不考虑界面效应时的结果；界面层厚度和性能对应力和有效反平面剪切模量明显变化时所对应的夹杂尺度范围和趋近于无界面效应结果的快慢有显著影响；当界面厚度足够薄时，界面层模型可用于模拟零厚度界面情况．     

压电材料反平面应变状态的任意形状夹杂问题

应用复函数的Ｆａｂｅｒ级数展开方法，分析了含任意形状夹杂的压电材料反平面应变问题，给出了问题的复势函数解。利用这个解，具体讨论了椭圆形夹杂及其极限（几何方面与物理方面）问题。并给出了三角形、正方形夹杂的近似结果。其特例结果与早期工作一致     

等参三角形热夹杂的构造及位移场数值计算

在材料制备和机械设计中,局部温升是造成材料失效和故障形成的重要因素之一.依照微观力学中,采用热夹杂模型可以定量深入地揭示与局部温升所关联的力学机理.在过往的研究中,受均匀热本征应变的夹杂模型广受关注;而相关非均匀分布的热本征应变问题,因其理论推导复杂而研究不多.论文首先给出在平面无限域中,受线性分布热本征应变作用的多边形夹杂的位移场解析解.基于格林函数法和围道积分,推导边界线单元的位移响应封闭解,该解通过叠加可直接给出线性热本征应变作用下的任意多边形夹杂的解析表达式.受到有限元分析中等参单元思想的启发,论文进一步将这种“等参元”方法扩展至求解Eshelby夹杂问题中.在该研究中,三角形单元的本征应变插值公式与位置坐标变换式均使用了相同的形函数与节点参数,因而所构建的单元模型称为等参三角形夹杂模型.论文方法可便捷地用于处理受任何分布热本征应变的任意形状二维Eshelby夹杂问题.相较于传统的有限元分析,论文所构建的数值求解方案实施方便且优势明显：只需在夹杂域上进行三角形网格剖分、而无需在无限的基体域上划分网格,因而可以极大地提高前处理便捷性及计算效率.此外,论文所给出的多边形夹杂解析解,...     

多圆形刚性夹杂的反平面问题

构造任意分布且相互影响的多个圆形刚性夹杂模型的复应力函数，采用复变函数方法，达到满足各个夹杂的边界条件，利用坐标变换和围线积分将求解方程组化为线性代数方程组，推导出了圆形刚性夹杂任意分布的界面应力解析表达式，算例对多夹杂与单夹杂两种模型的界面应力最大值进行了对比，同时还给出了界面应力最大值随夹杂间距的变化规律，求出了刚性夹杂的合理间距。本文发展的分析方法为研究夹杂材料的细观机理探索了一条有效的分析途径。     

压电圆柱形夹杂对邻近裂纹三维应力场的影响

研究了压电复合材料薄板中压电圆柱形夹杂与邻近宏观钝裂纹间的相互作用。重点分析了外加电场，裂尖与压电圆柱形夹杂间韧带长度对裂尖三维应力场的影响。计算结果表明：在不同的外加电场作用下，压电体不仅能改变裂尖张开应力的大小，还能改变其分布。所得结果对进一步探讨线弹性介质中裂纹的启裂控制有参考价值。     

On the interaction between a dislocation and a circular inhomogeneity with imperfect interface in antiplane shear

The solution of appropriate elasticity problems involving the interaction between inclusions and dislocations plays a fundamental role in many practical and theoretical applications, namely, it increases the understanding of material defects thereby providing valuable insight into the mechanical behavior of composite materials.Although the problem of a three-phase circular inclusion interacting with a dislocation in antiplane shear has been presented [Xiao and Chen, Mech. Mater. 32 (2000) 485], the analysis is limited to the classical perfect bonding condition. The current paper considers the solution for a homogeneous circular inclusion interacting with a dislocation under thermal loadings in antiplane shear. The bonding along the inhomogeneity–matrix interface is considered to be imperfect with the assumption that the interface imperfections are constant. It is found that when the inhomogeneity is soft, regardless of the level of interface imperfection, the inhomogeneity will always attract the dislocation. As a result, no equilibrium positions are available. Alternatively, when the inhomogeneity is hard, an unstable equilibrium position is found which depends on the imperfect interface condition and the shear moduli ratio μ₂/μ₁.     

The effective properties of piezocomposites, part I: Single inclusion problem

The problem of a piezoelectric ellipsoidal inclusion in an infinite nonpiezoelectric matrix is very important in engineering. In this paper, it is solved via Green's function technique. The closed-form solutions of the electroelastic Eshelby's tensors for this kind of problem are obtained. The electroelastic Eshelby's tensors can be expressed by the Eshelby's tensors of the perfectly elastic inclusion problem and the perfectly dielectric inclusion problem. Since the closed-form solutions of the Eshelby's tensors of the perfectly elastic inclusion problem and the perfectly dielectric inclusion problem can be given by theory of elasticity and electrodynamics, respectively, the electroelastic Eshelby's tensors can be obtained conveniently. Using these results, the closed-form solutions of the constraint elastic fields and the constraint electric fields inside the piezoelectric ellipsoidal inclusion are also obtained. These expressions can be readily utilized in solutions of numerous problems in the micromechanics of piezoelectric solids, such as the deformation and energy analysis, damage evolution and fracture of the piezoelectric materials. The project supported by the National Natural Science Foundation of China     

ANALYSIS OF A SCREW DISLOCATION INSIDE A CIRCULAR INCLUSION WITH INTERFACIAL CRACKS IN PIEZOELECTRIC COMPOSITES

IntroductionPiezoelectric materials have potentials for use in many modern devices and compositestructures. The presence of various defects, such as inclusions, holes, dislocations andcracks, can greatly influence their characteristics and coupled behavio…     

On a semi-infinite crack penetrating a piezoelectric circular inhomogeneity with a viscous interface

We investigate a semi-infinite crack penetrating a piezoelectric circular inhomogeneity bonded to an infinite piezoelectric matrix through a linear viscous interface. The tip of the crack is at the center of the circular inhomogeneity. By means of the complex variable and conformal mapping methods, exact closed-form solutions in terms of elementary functions are derived for the following three loading cases: (i) nominal Mode-III stress and electric displacement intensity factors at infinity; (ii) a piezoelectric screw dislocation located in the unbounded matrix; and (iii) a piezoelectric screw dislocation located in the inhomogeneity. The time-dependent electroelastic field in the cracked composite system is obtained. Particularly the time-dependent stress and electric displacement intensity factors at the crack tip, jumps in the displacement and electric potential across the crack surfaces, displacement jump across the viscous interface, and image force acting on the piezoelectric screw dislocation are all derived. It is found that the value of the relaxation (or characteristic) time for this cracked composite system is just twice as that for the same fibrous composite system without crack. Finally, we extend the methods to the more general scenario where a semi-infinite wedge crack is within the inhomogeneity/matrix composite system with a viscous interface.     

Debonded arc-shaped interface conducting rigid line inclusions in piezoelectric composites

We consider an arc-shaped conducting rigid line inclusion located at the interface between a circular piezoelectric inhomogeneity and an unbounded piezoelectric matrix subjected to remote uniform anti-plane shear stresses and in-plane electric fields. Moreover, one side of the rigid line inclusion has become fully debonded from the matrix or the inhomogeneity leading to the formation of an insulating crack. After the introduction of two sectionally holomorphic vector functions, the problem is reduced to a vector Riemann–Hilbert problem, which can be decoupled sequentially by repeated application of the orthogonality relations between the eigenvectors for two corresponding generalized eigenvalue problems.     

双周期圆截面纤维复合材料平面问题的解析法

结合双准周期Riemann边值问题理论与Eshelby等效夹杂原理，为双周期圆截面纤维复合材 料平面问题发展了一个实用有效的解析方法，获得了问题的全场级数解并与有限元结果进行 了比较. 该方法为非均匀材料的力学性质分析和复合材料等新材料的微结构设计提供了 一个有效的计算工具，也可用来评估有限元等数值与近似方法的精度.     

Thermoelectric field for a coated hole of arbitrary shape in a nonlinearly coupled thermoelectric material

We study the thermoelectric field for an electrically and thermally insulated coated hole of arbitrary shape embedded in an infinite nonlinearly coupled thermoelectric material subject to uniform remote electric current density and uniform remote energy flux. A conformal mapping function for the coating and matrix is introduced, which simultaneously maps the hole boundary and the coating-matrix interface onto two concentric circles in the image plane. Using analytic continuation, we derive a general solution in terms of two auxiliary functions. The general solution satisfies the insulating conditions along the hole boundary and all of the continuity conditions across the perfect coating-matrix interface. Once the two auxiliary functions have been obtained in the elementary-form, the four original analytic functions in the coating and matrix characterizing the thermoelectric fields are completely and explicitly determined. The design of a neutral coated circular hole that does not disturb the prescribed thermoelectric field in the thermoelectric matrix is achieved when the relative thickness parameter and the two mismatch parameters satisfy a simple condition. Finally, the neutrality of a coated circular thermoelectric inhomogeneity is also accomplished.     

Electro-elastic interaction between a piezoelectric screw dislocation and circular interfacial rigid lines

The electro-elastic interaction between a piezoelectric screw dislocation located either outside or inside inhomogeneity and circular interfacial rigid lines under anti-plane mechanical and in-plane electrical loads in linear piezoelectric materials is dealt with in the framework of linear elastic theory. Using Riemann–Schwarz’s symmetry principle integrated with the analysis of singularity of complex functions, the general solution of this problem is presented in this paper. For a special example, the closed form solutions for electro-elastic fields in matrix and inhomogeneity regions are derived explicitly when interface containing single rigid line. Applying perturbation technique, perturbation stress and electric displacement fields are obtained. The image force acting on piezoelectric screw dislocation is calculated by using the generalized Peach–Koehler formula. As a result, numerical analysis and discussion show that soft inhomogeneity can repel screw dislocation in piezoelectric material due to their intrinsic electro-mechanical coupling behavior and the influence of interfacial rigid line upon the image force is profound. When the radian of circular rigid line reaches extensive magnitude, the presence of interfacial rigid line can change the interaction mechanism.     

Analysis of a partially debonded elliptic inhomogeneity in piezoelectric materials

A generalized solution was obtained for the partially debonded elliptic inhomogeneity problem in piezoelectric materials under antiplane shear and inplane electric loading using the complex variable method. It was assumed that the interfacial debonding induced an electrically impermeable crack at the interface. The principle of conformal transformation and analytical continuation were employed to reduce the formulation into two Riemann-Hilbert problems. This enabled the determination of the complex potentials in the inhomogeneity and the matrix by means of series of expressions. The resulting solution was then used to obtain the electroeiastic fields and the energy release rate involving the debonding at the inhomogeneity-matrix interface. The validity and versatility of the current general solution have been demonstrated through some specific examples such as the problems of perfectly bonded elliptic inhomogeneity , totally debonded elliptic inhomogeneity, partially debonded rigid and conducting elliptic inhomogeneity, and partially debonded circular inhomogeneity.     

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.