有限长压电层合简支板自由振动的三维精确解

基于三维弹性理论和压电理论，导出了有限长矩形压电层合简支板的动力学方程及相应的边界条件，给出了一种求解压电层合板自由振动三维精确解的方法；分析了正、逆向压电效应对层合板振动频率的影响．本文所述的方法和结果对于求解其他三维动态问题，验证、比较其他简化模型、有限元计算结果以及工程应用都有指导意义．     

基于遗传算法和梯度算法压电层合结构的最优形状控制

基于压电层合结构的有限元方程,运用ANSYS/APDL语言,编制了力-电耦合有限元分析程序(MPFEMP).以该程序为计算基础,采用遗传算法和一阶梯度优化算法,以压电片尺寸为设计变量,以压电层合梁和板的预期位移或最小重量为目标函数,给定初始变量和适应度函数,通过循环迭代MPFEMP计算程序,研究了多点控制的压电层合梁板结构的形状最优控制.结论对比分析证明了两种优化方法分析压电层合结构的有效性,同时,对复杂多层智能结构的最优形状控制和主动控制研究具有一定的参考价值.     

采用压电机敏元件进行结构振动控制Ⅲ：控制系统设计与实验研究

采用作者在上篇导出的压电耦合体动力学模型，给出了压电主动阻尼控制系统的设计方法；导出了压电耦合梁系统的作动方程和检测方程的显式表达式。以此为基础，以简单梁为对象，对压电检测器和作动器的性能、粘结层的影响、压电主动阻尼控制及压电主、被动阻尼双控制进行了实验研究     

层合厚板混合状态Hamiltonian动力元及其半解析解

文中在时间方向采用Ｌａｐｌａｃｅ变换，给出了层合厚板动力学问题混合状态Ｈａｍｉｌｔｏｎ正则方程及其半解析法．该方法在层板平面内采用通常的有限元离散，而沿板厚方向采用状态控制方程给出解析解答．在层与层之间采用迁移矩阵法，给出层合板上下表面力学量之间的关系式．利用打靶法得到响应在象空间的一般解，然后再利用拉氏逆变换的数值解求出层合板的瞬时位移场和应力场．     

压电复合材料层合梁的分岔、混沌动力学与控制

研究了简支压电复合材料层合梁在轴向、横向载荷共同作用下的非线性动力学、分岔和混沌动力学响应. 基于vonKarman理论和Reddy高阶剪切变形理论,推导出了压电复合层合梁的动力学方程. 利用Galerkin法离散偏微分方程,得到两个自由度非线性控制方程,并且利用多尺度法得到了平均方程. 基于平均方程,研究了压电层合梁系统的动态分岔,分析了系统各种参数对倍周期分岔的影响及变化规律. 结果表明,压电复合材料层合梁周期运动的稳定性和混沌运动对外激励的变化非常敏感,通过控制压电激励,可以控制压电复合材料层合梁的振动,保持系统的稳定性,即控制系统产生倍周期分岔解,从而阻止系统通过倍周期分岔进入混沌运动,并给出了控制分岔图.      

主动约束阻尼板压电层电压拓扑优化研究

建立主动约束层阻尼板有限元模型,以结构模态阻尼比最大化为目标函数,压电层总电能消耗为约束条件,压电层单元控制电压为设计变量,对主动约束层阻尼板压电层电压进行了拓扑优化,获得了压电层电压最优拓扑分布。通过引入虚拟设计变量,将压电层电压控制不连续问题转化为连续问题。考虑实际工程应用的需要,采用指数函数对电压中间变量进行惩罚。在灵敏度分析基础上,采用移动渐进线(MMA)法,求解了主动约束层阻尼板电压拓扑优化问题。数值算例证实了电压拓扑优化模型以及数值求解方法的有效性。     

基于应变梯度有限元的单层石墨烯振动研究

建立了单层石墨烯等效非局部薄板的一种新的有限元模型,并运用有限元法分析不同边界条件下单层石墨烯振动的小尺度效应。给出了基于弹性应变梯度理论下Kirchhoff板的振动方程。发展了一种4节点24自由度的板单元,用于离散化求解考虑微纳结构尺度效应的高阶微分方程。在研究四边简支板振动时,考虑应变梯度的非局部弹性有限元数值计算结果与理论分析结果相一致。用有限元方法研究了不同尺寸、振动波长、振动模态阶数、边界条件类型以及非局部参数的单层石墨烯振动。     

含压电片层合壳的有限元分析与控制仿真

本文首先针对含压电材料的一般结构，推导了按位移和电位移求解时的混合变分原理，在此基础上，通过对面电荷取变分，直接得到了含按电压驱动压电片层合壳的有限元方程，最后给出了利用压电片进行静变形和振动控制的仿真算例，并进行了分析。结果表明，本文建立的有限元方程准确、可靠。     

功能梯度压电圆板轴对称自由振动问题精确解

将功能梯度压电圆板的位移变量和电势变量写为分离变量的形式，由压电动力学平衡方程导出以位移、电势及其一阶导数为状态变量的状态方程，考虑周边固支接地的边界条件，导出了求解功能梯度压电圆板自振频率精确解的方程。将方程退化至一般的非梯度纯弹性圆板的形式，求解其自振频率，得到的结果与相应的理论解完全吻合，从而验证了本文方法的正确性。更进一步地对梯度函数沿板厚以指数形式变化的功能梯度压电圆板的自振频率进行了计算，并得到了梯度化对板自振频率的影响规律。     

压电结构的热弹性比拟建模方法

通过比较逆压电效应本构方程和热弹性本构方程,在线性范围内建立逆压电效应和热弹性效应的比拟关系:将驱动电压比拟为温度载荷,压电应变系数比拟为热膨胀系数,用热弹性有限元方法分析逆压电效应,并采用热弹性比拟方法求解基于长度伸缩逆压电效应和厚度剪切逆压电效应的驱动器问题,给出了平面结构“热-机-电”耦合问题的热弹性比拟求解方法。热弹性比拟方法为快速建立复杂压电结构的有限元模型提供了一种有效的途径,可缩短复杂压电主动结构的设计周期,降低设计成本。最后给出数值算例说明文中方法的有效性。     

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.     

横观各向同性压电材料中裂纹问题的边界元法

采用Somigiliana公式给出了三维横观各向同性压电材料中的非渗漏裂纹问题的一般解和超奇异积分方程,其中未知函数为裂纹面上的位移间断和电势间断．在此基础上,使用有限部积分和边界元结合的方法,建立了超奇异积分方程的数值求解方法,并给出了一些典型数值算例的应力强度因子和电位移强度因子的数值结果,结果令人满意．     

PZT-4紧凑拉伸试样的断裂分析

基于线性压电材料的复势理论,通过解析分析,导出了一种分析有限压电板裂纹问题的解析数值方法. 首先,计算了含中心裂纹有限板的断裂参数,与Woo和Wang的解析数值法(Int J Fract, 1993, 62: 203$\sim$218)相比较,表明该方法具有很高的精度和很好的计算效率. 随后,采用该方法和有限元法计算了PZT-4紧凑拉伸试样在绝缘裂纹面边界条件下断裂时的断裂参数,发现各断裂参数的临界值分散性很大,不能作为压电材料的单参数断裂准则. 进而,针对试样真实的裂隙形状,采用有限元法计算了裂隙尖端的应力、电位移场,比较了裂隙内介质的介电性能对裂隙尖端场的影响,计算了带微裂纹的真实裂隙模型的断裂参数并进行了理论分析.      

Basic solutions of a 3-D rectangular limited-permeable crack or two 3-D rectangular limited-permeable cracks in?piezoelectric materials

The solutions of a 3-D rectangular limited-permeable crack or two 3-D rectangular limited-permeable cracks in piezoelectric materials were given by using the generalized Almansi’s theorem and the Schmidt method. At the same time, the electric permittivity of the air inside the rectangular crack was considered. The problem was formulated through Fourier transform as three pairs of dual integral equations, in which the unknown variables are the displacement jumps across the crack surfaces. To solve the dual integral equations, the displacement jumps across the crack surfaces were directly expanded as a series of Jacobi polynomials. Finally, the effects of the electric permittivity of the air inside the rectangular crack,the shape of the rectangular crack and the distance between two rectangular cracks on the stress and electric displacement intensity factors in piezoelectric materials were analyzed.     

A 3-D rectangular permeable crack or two 3-D rectangular permeable cracks in a piezoelectric material

The solutions of a 3-D rectangular permeable crack and two 3-D rectangular permeable cracks in a piezoelectric material were investigated by using the generalized Almansi’s theorem and the Schmidt method. The problem was formulated through Fourier transform into three pairs of dual integral equations, in which the unknown variables are the displacement jumps across the crack surfaces. To solve the dual integral equations, the displacement jumps across the crack surfaces were directly expanded as a series of Jacobi polynomials. Finally, the effects of the shape of the rectangular crack and the distance between two rectangular cracks on the stress and electric displacement intensity factors in a piezoelectric material were analyzed. These results can be used for the strength and the coupling effect evaluation of cracked piezoelectric materials.     

An effective boundary element method for analysis of crack problems in a plane elastic plate

A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.     

内部压力作用下矩形板中源于椭圆孔的分支裂纹应力强度因子的一种数值分析

应用一种边界元方法来研究内部压力作用下矩形板中源于椭圆孔的分支裂纹。该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者最近提出的裂尖位移不连续单元构成。在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界。本数值结果进一步证实这种数值方法对计算有限大板中复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示裂纹体几何对应力强度因子的影响。     

刚度微分法计算压电材料平面断裂问题

把计算应变能释放率的刚度微分法推广到压电材料平面断裂问题．在此基础上，利用压电材料平面断裂问题的有限元数值解作为真实场，用Ｓｏｓａ的平面问题裂端渐近解作为辅助场，由推广的交互Ｍ积分法求得了应力强度因子ＫＩ，ＫＩＩ和电位移强度因子ＫＩＶ．算例表明，计算结果与理论解符合得很好     

Green function and its application for a piezoelectric plate with various openings

Summary For a two-dimensional piezoelectric plate subjected to mechanical and electric load, a Green function satisfying traction free and exact electric boundary conditions along a hole is developed using Lekhnitskii's formalism and the technique of conformal mapping. The critical points for the mapping function used is investigated numerically, and the study indicates that the transformation of a polygonal hole in a piezoelectric plate into a unit circle is nonsingle-valued. A simple approach is presented to treat such a situation. Based on the Green function developed in this paper, a system of singular integral equations for the unknown dislocation defined on crack faces is presented to study the interaction between cracks and holes. Numerical results are presented to elucidate the effects of crack orientation on stress and electric displacement (SED) intensity factors and to illustrate the application of the proposed formulation. Received 4 June 1998, accepted for publication 15 July 1988     

THE BASIC SOLUTION OF A 3-D RECTANGULAR PERMEABLE CRACK IN A PIEZOELECTRIC/PIEZOMAGNETIC COMPOSITE MATERIAL

The solution of a 3-D rectangular permeable crack in a piezoelectric/piezomagnetic composite material was investigated by using the generalized Almansi’s theorem and the Schmidt method.The problem was formulated through Fourier transform into three pairs of dual integral equations,in which the unknown variables are the displacement jumps across the crack surfaces.To solve the dual integral equations,the displacement jumps across the crack surfaces were directly expanded as a series of Jacobi polynomials.Finally,the relations between the electric filed,the magnetic flux field and the stress field near the crack edges were obtained and the efects of the shape of the rectangular crack on the stress,the electric displacement and magnetic flux intensity factors in a piezoelectric/piezomagnetic composite material were analyzed.