各向同性与各向异性三相材料接头应力奇异性研究

研究了各向同性与各向异性三相材料接头的应力奇性指数,通过引入奇异点附近区域位移场渐近展开的典型项,将各向同性与各向异性组合材料接头的控制方程和径向边界条件转化为变系数常微分方程的特征值问题;再利用插值矩阵法求解所建立的特征方程,得到接头端部的应力奇性指数和特征角函数。对由两个各向异性材料和一个各向同性材料以任意楔形角组成的三相接头结构的奇异性进行了研究,并比较了它们的应力奇性指数。计算结果表明:对于粘结接头,各向同性材料刚度越大应力奇异性越强;对于剥离接头,各向同性材料楔形角或材料刚度越大,第一阶应力奇异性越弱。计算结果与已有文献的结果对比吻合良好,证明了本文方法的有效性。     

层状压电陶瓷致动器中力电耦合场奇异性的数值分析

首先推导了不同压电材料界面裂纹尖端处的扇形区域内包含基本方程、裂纹面D-P边界条件和交界面处边界条件的弱形式。通过假设力电耦合位移场(位移和电势)与到裂纹尖端距离的(λ 1)次方成正比,可以得到一个分析压电材料裂纹尖端处力电耦合场奇异性的特殊的一维有限元列式。该一维有限元列式只需对扇形区域在角度方向上离散,最后的总体方程为一个关于λ的二次特征根方程。探讨了层状压电陶瓷致动器中可能出现奇异力电耦合场的部位的裂纹面边界条件及交界面处边界条件,进而将该一维有限元法进行推广,用于研究了这些部位的力电耦合场的奇异性。通过数值算例与相应的精确解的比较表明该方法是正确的,而且仅用很少单元就可以得到非常精确的结果。     

PMN-PT层/弹性基底结构中声表面波特性分析

研究了PMN-PT 压电层/弹性(金刚石) 基底结构中表面波的传播特性,压电层表面是机械自由的,电学边界条件分为电学开路和电学短路,压电层与基底之间采用理想连接. 得到了满足控制方程和边界条件的电弹场以及弹性波在结构中传播时的频散方程,通过数值算例分析了压电材料PMN-PT 的极化方向对弹性波频散曲线和机电耦合系数的影响,以及不同极化方向时弹性位移和电势随结构深度方向的变化,结果可为PMN-PT 压电材料在高频声表面波器件中的应用提供有价值的理论参考.      

压电材料裂纹顶端条状电饱和区模型的力学分析

在线性压电本构方程框架下,对裂纹顶端条状电饱和区模型进行了严格的数学分析．完整地考虑了各向异性力电耦合效应．建立了电饱和区尺寸与外加电场的依赖关系．证实了当裂纹垂直极化轴时,压电材料的断裂应力随着外加正电场的增加而减小,随着外加负电场的增加而增加．当裂纹平行于极化轴时,与极化轴平行的外加电场对断裂应力无影响     

压电材料裂纹顶端条状电饱和区模型的力学分析

在线性压电本构方程框架下,对裂纹顶端条状电饱和区模型进行了严格的数学分析．完整地考虑了各向异性力电耦合效应．建立了电饱和区尺寸与外加电场的依赖关系．证实了当裂纹垂直极化轴时,压电材料的断裂应力随着外加正电场的增加而减小,随着外加负电场的增加而增加．当裂纹平行于极化轴时,与极化轴平行的外加电场对断裂应力无影响     

正交异性多材料三维接头应力奇异性分析

不同材料组成的接头结构在实际工程中是很常见的,这些应力集中较高的接头常常是工程结构断裂失效所处的最关键位置。基于正交异性多材料接头端部物理场的幂级数渐近展开假设,由弹性力学理论导出了关于正交异性多材料接头端部应力奇性指数的特征微分方程组,并将接头端部的力学交界和边界条件表达为奇性指数和特征角函数的组合,从而将正交异性多材料接头端部应力奇性指数的计算转化为相应边界条件下常微分方程组特征值的求解问题,运用插值矩阵法求解可以获得多材料接头端部若干阶应力奇性指数和相应的特征角函数。数值计算结果与现有结果对比表明了该方法的有效性和具有较高的计算精度;同时计算结果还表明多材料接头结构发生断裂失效的主要原因是位移场特征角函数的一阶导函数在不同材料粘结界面发生了突变。     

各向同性材料接头和界面相交裂纹应力奇异性特征分析

提出了用插值矩阵法分析各向同性材料接头以及与界面相交的平面裂纹应力奇异性。基于接头和裂纹端部附近区域位移场渐近展开,将位移场的渐近展开式的典型项代入线弹性力学基本方程,得到关于平面内各向同性材料接头以及与两相材料界面相交裂纹应力奇异性指数的一组非线性常微分方程的特征值问题,运用插值矩阵法求解,获得了两相材料平面接头端部应力奇异性指数以及与界面以任意角相交的裂纹尖端的应力奇异性指数随裂纹角的变化规律,数值计算结果与已有结果比较表明,本文方法具有很高的精度和效率。     

功能梯度压电圆板自由振动问题的三维精确分析

本文对周边为广义刚性滑动和广义简支两种边界条件下的功能梯度压电材料圆板自由振动问题进行分析。根据轴对称横观各向同性压电材料基本方程，并利用有限Hankel变换得到了功能梯度压电材料圆板的状态空间方程。假设材料的机械和电学性质均沿板厚方向按统一的指数函数形式梯度分布，从而获得了周边为广义刚性滑动和广义弹性简支两种边界条件下功能梯度压电圆板自由振动问题的三维精确频率方程，该方程是一个关于自由振动频率的超越方程，通过求解该超越方程可得到在不同板厚以及不同的材料性质梯度变化情况下的圆板自由振动频率值，结果表明在相同的材料性质梯度变化情况下频率均随着板厚增加而增大，而在相同的板厚情况下频率则随材料性质梯度变化指数的增大而减小的结论。     

压电材料反平面切口奇性指数分析

为理解压电材料反平面切口尖端奇异状态,提出了一种切口奇性特征分析法.基于切口根部位移场幂级数渐近展开假设,从应力平衡方程和电荷守恒条件出发,导出了关于压电材料反平面切口奇性指数的特征微分方程组,并将切口的力电学边界条件以及界面协调条件表达为奇性指数和特征角函数的组合.从而,压电材料切口反平面奇性指数的计算被转化为相应边界条件下常微分方程组特征值的求解问题,采用插值矩阵法可以计算出各阶奇性指数和相应的特征角函数.该法既适合裂纹奇性分析,也可用于单、双材料切口的奇性计算,并避免了用迭代法求解超越方程的不足.因而具有适应性强的特点.计算发现,压电材料反平面切口存在两个奇性指数,切口的奇异性程度随着开角的增大而增强.     

轴对称压电-压磁圆柱夹层结构的圣维南端部效应

研究了半无限长轴对称压电-压磁夹层结构的圆柱体圣维南端部效应的衰减问题。圆柱的端部承受自平衡磁电弹载荷;圆柱的内外表面为机械自由表面,但承受不同的电磁边界条件,即电学短路或电学开路及磁学短路或磁学开路边界条件。基于横贯各向同性压电或压磁材料在轴对称圆柱坐标系下的本构方程,推导了关于衰减率的特征方程并求得问题的数值解。结果表明,边界条件、内外径之比、材料厚度比对结构的衰减率都有显著的影响。     

基于非协调元特征法的奇异性问题分析

提出了一个基于位移的、分析平面尖劈尖端奇性应力场和位移场问题的非协调FE特征分析法．该方法与过去原有求解裂纹尖端近似场的有限元特征分析方法导出公式的出发点不同，并且采用的单元形式为非协调元，尖劈尖端邻域内的位移场假定没有采用奇异变换技术，运用该方法处理了若干尖劈和接头的算例，所有的计算结果表明，该方法较原有方法使用的单元少而且精度高，具有应用广泛性。     

基于新型裂尖杂交元的压电材料断裂力学研究

提出了一种裂尖邻域杂交元模型，将其与标准杂交应力元结合来求解压电材料裂纹尖 端的奇性电弹场和断裂参数的数值解．裂纹尖端杂交元的建立步骤为：1) 利用高次内插有限元特征法求解特征问题，得到反映裂尖奇异性电弹场状况的特 征值和特征角分布函数；2) 利用广义Hellinger-Reissner变分泛函以及特征问题的解来建立裂尖邻域杂交元模型．该 方法求解电弹场时，摒弃了传统有限元方法中裂尖奇异性场需要借助解析解的做法，也避免 了单纯有限元方法中需要在裂尖端部进行高密度单元划分．采用PZT5板中心裂纹问题 作为考核例，数值结果显示了良好的精确性．作为进一步应用，求解了含中心界面裂纹 的PZT4-PZT5两相压电材料的应力强度因子和电位移强度因子．所有的算例都考虑 了3种裂纹面电边界条件．     

2D dynamic analysis of cracks and interface cracks in piezoelectric composites using the SBFEM

The dynamic stress and electric displacement intensity factors of impermeable cracks in homogeneous piezoelectric materials and interface cracks in piezoelectric bimaterials are evaluated by extending the scaled boundary finite element method (SBFEM). In this method, a piezoelectric plate is divided into polygons. Each polygon is treated as a scaled boundary finite element subdomain. Only the boundaries of the subdomains need to be discretized with line elements. The dynamic properties of a subdomain are represented by the high order stiffness and mass matrices obtained from a continued fraction solution, which is able to represent the high frequency response with only 3–4 terms per wavelength. The semi-analytical solutions model singular stress and electric displacement fields in the vicinity of crack tips accurately and efficiently. The dynamic stress and electric displacement intensity factors are evaluated directly from the scaled boundary finite element solutions. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required. Numerical examples are presented to verify the proposed technique with the analytical solutions and the results from the literature. The present results highlight the accuracy, simplicity and efficiency of the proposed technique.     

电弹性体力学中的偏场方法及其应用

主要综述了当前有关叠加于偏场之 上的电弹性体小位移问题的求解方法,即偏场方法.首先介绍了作为偏场方法理论基 础的非线性电弹性力学理论,接着总结了偏场方法的研究进展以及 受偏场作用后,电弹性梁、板、壳结构的分析方法,随后综述了偏场方法的诸多应用:其中, 包括在薄壁电弹性结构屈曲分析中的应用、在记时与通信压电谐振器和基于频率漂移原理所 设计的声波传感器的频率稳定性分析方面的应用、在非线性电弹性材料系数的测定以及偏场 作用下电致伸缩陶瓷特征的分析等方面的应用.最后给出了该领域当前和未来的一些可能的 研究课题.全文参考文献166篇.     

The effective electroelastic property of piezoelectric media with parallel dielectric cracks

Existing studies on the coupled electroelastic behaviour of cracked piezoelectric media have been based mostly on the electrically impermeable and permeable crack models. The current paper presents a study of the effective electroelastic property of piezoelectric media weakened by parallel cracks using a dielectric crack model with the electric boundary condition along the crack surfaces being governed by the opening displacement. The theoretical formulation is obtained using the dilute model of distributed cracks and the solution of a single dielectric crack problem. It is observed that the effective electroelastic property of cracked piezoelectric media is nonlinear and sensitive to loading conditions. Different modes of crack deformation are predicted and discussed. Attention is paid to the transition between electrically permeable and impermeable crack models.     

On the singularities of the thermo-electro-elastic fields near the apex of a piezoelectric bonded wedge

Based on the generalized Lekhnitskii formulation and Mellin transform, the thermo-electro-elastic fields of a piezoelectric bonded wedge are investigated in this paper. From the potential theory in a wedge-shaped region, a general form of the temperature change is proposed as a particular solution in the generalized Lekhnitskii formulation. The emphasis is on the singular behavior near the apex of the piezoelectric bonded wedge, including singularity orders and angular functions, which can be computed numerically. The interface between two materials can be either perfectly bonded, namely type A, so that the continuity of electric displacements holds, or a thin electrode, namely type B, so that the electric potential is grounded. Case studies of PZT-5H/PZT-4 and graphite-epoxy/PZT-4 bonded wedges reveal that, in most cases, the type B continuity condition has more severe singularities than type A due to the mixed boundary point of the electrostatics at the apex of the wedge. The results of this study show that the reduction or disappearance of singularity orders is possible through the appropriate selection of poling/fiber orientations and wedge angles.     

压电材料反平面界面裂纹的杂交元分析

基于奇异性电弹场数值特征解开发了一种新型反平面界面裂纹尖端单元。将新型单元与四节点压电P-S单元组装,求解从绝缘到导通的任意电边界条件下,压电结构反平面界面裂纹尖端电弹场的数值解。考察了层厚、载荷类型和裂纹面间电边界条件等对反平面界面裂纹尖端断裂参数的影响。算例证明新型单元能使P-S单元数显著降低,计算结果更为精确。     

Analysis of bonded anisotropic wedges with interface crack under anti-plane shear loading

The antiplane stress analysis of two anisotropic finite wedges with arbitrary radii and apex angles that are bonded together along a common edge is investigated. The wedge radial boundaries can be subjected to displacement-displacement boundary condi- tions, and the circular boundary of the wedge is free from any traction. The new finite complex transforms are employed to solve the problem. These finite complex transforms have complex analogies to both kinds of standard finite Mellin transforms. The traction free condition on the crack faces is expressed as a singular integral equation by using the exact analytical method. The explicit terms for the strength of singularity are extracted, showing the dependence of the order of the stress singularity on the wedge angle, material constants, and boundary conditions. A numerical method is used for solving the resul- tant singular integral equations. The displacement boundary condition may be a general term of the Taylor series expansion for the displacement prescribed on the radial edge of the wedge. Thus, the analysis of every kind of displacement boundary conditions can be obtained by the achieved results from the foregoing general displacement boundary condition. The obtained stress intensity factors （SIFs） at the crack tips are plotted and compared with those obtained by the finite element analysis （FEA）.     

Damage analysis of thermopiezoelectric properties: Part I - crack tip singularities

This is Part I of the work on a two-dimensional analysis of thermal and electric fields of a thermopiezoelectric solid damaged by cracks. It deals with finding the singular crack tip behavior for the temperature, heat flow, displacements, electric potential, stresses and electric displacements. By application of Fourier transformations and the extended Stroh formalism, the problem is reduced to a pair of dual integral equations for the temperature field with the aid of an auxiliary function. The electroelastic field is governed by another pair of dual integral equations. The inverse square root singularity is found for the heat flow field while the logarithmic singularity prevailed for the electroelastic field regardless of whether the crack lies in a homogeneous piezoelectric solid or at an interface of two dissimilar piezoelectric materials. Results are given for the energy release rate and a finite length crack oriented at an arbitrarily angle with reference to the external disturbances. Part II of this paper considers the modelling of a piezoelectric material containing microcracks. A representative cracked area element is used to obtain the effective conductivity and electroelastic modulus. Numerical results are given for a peizoelectric Bati O₃ ceramic with cracks.     

Dynamic fracture behavior of a cracked piezoelectric half space under anti-plane mechanical and in-plane electric impact

Summary  The dynamic response of a cracked piezoelectric half-space under anti-plane mechanical and in-plane electric impacting loads is investigated in the present paper. In the study, the crack is assumed parallel to the free surface of the half-space. Laplace and Fourier transforms are used to reduce the mixed boundary value problems to Cauchy-type singular integral equations in the Laplace transform domain, which are solved numerically. Then, a numerical Laplace inversion is performed and the dynamic stress and electric displacement factors are obtained as functions of time and geometry parameters. The dynamic energy release rate is derived for piezoelectric materials in terms of the electroelastic intensities and is displayed graphically. Received 5 January 2000; accepted for publication 28 June 2000