ELECTROELASTIC INTENSIFICATION NEAR ANTI-PLANE CRACK IN A FUNCTIONALLY GRADIENT PIEZOELECTRIC CERAMIC STRIP

Following the theory of linear piezoelectricity, we consider the electro-elastic problems of a finite crack in a functionally gradient piezoelectric ceramic strip. By the use of Fourier transforms we reduce the problem to solving two pairs of dual integral equations. The solution to the dual integral equations is then expressed in terms ofa Fredholm integral equation of the second kind. Numerical calculations are carried out for piezoelectric ceramics. The electric field intensity factors and the energy release rate are shown graphically, and the electroelastic interactions are illustrated.     

COUPLED ANALYSIS FOR THE HARVESTING STRUCTURE AND THE MODULATING CIRCUIT IN A PIEZOELECTRIC BIMORPH ENERGY HARVESTER

The authors analyze a piezoelectric energy harvester as an electro-mechanically coupled system. The energy harvester consists of a piezoelectric bimorph with a concentrated mass attached at one end, called the harvesting structure, an electric circuit for energy storage, and a rectifier that converts the AC output of the harvesting structure into a DC input for the storage circuit. The piezoelectric bimorph is assumed to be driven into flexural vibration by an ambient acoustic source to convert the mechanical energies into electric energies. The analysis indicates that the performance of this harvester, measured by the power density, is characterized by three important non-dimensional parameters, i.e., the non-dimensional inductance of the storage circuit, the non-dimensional aspect ratio （length/thickness） and the non-dimensional end mass of the harvesting structure. The numerical results show that： （1） the power density can be optimized by varying the non-dimensional inductance for each fixed non-dimensional aspect ratio with a fixed non-dimensional end mass; and （2） for a fixed non-dimensional inductance, the power density is maximized if the non-dimensional aspect ratio and the non-dimensional end mass are so chosen that the harvesting structure, consisting of both the piezoelectric bimorph and the end mass attached, resonates at the frequency of the ambient acoustic source.     

ANALYSIS OF CRACK-TIP SINGULARITIES FOR AN INTERFACIAL PERMEABLE CRACK IN METAL/PIEZOELECTRIC BIMATERIALS

By modeling metal as a special piezoelectric material with extremely small piezoelec- tricity and extremely large permittivity,we have obtained the analytical solutions for an interfacial permeable crack in metal/piezoelectric bimaterials by means of the generalized Stroh formalism. The analysis shows that the stress fields near a permeable interfacial crack tip are usually with three types of singularities:r~(-1/2 iε)and r~(-1/2).Further numerical calculation on the oscillatory indexεare given for 28 types of metal/piezoelectric bimaterials combined by seven commercial piezoelectric materials: PZT-4,BaTiO_3,PZT-5H,PZT-6B,PZT-7A,P-7 and PZT-PIC 151 and four metals:copper,silver,lead and aluminum,respectively.The explicit expressions of the crack tip energy release rate(ERR)and the crack tip generalized stress intensity factors(GSIF)are obtained.It is found that both the ERR and GSIF are independent of the electric displacement loading,although they seriously depends on the mechanical loadings.     

A novel type of transverse surface wave propagating in a layered structure consisting of a piezoelectric layer attached to an elastic half-space

The existence and propagation of transverse surface waves in piezoelectric coupled solids is investigated, in which perfect bonding between a metal/dielectric substrate and a piezoelectric layer of finite-thickness is assumed. Dis- persion equations relating phase velocity to material con- stants for the existence of various modes are obtained in a simple mathematical form for a piezoelectric material of class 6mm. It is discovered and proved by numerical examples in this paper that a novel Bleustein-Gulyaev （B-G） type of transverse surface wave can exist in such piezoelectric cou- pled solid media when the bulk-shear-wave velocity in the substrate is less than that in the piezoelectric layer but greater than the corresponding B-G wave velocity in the same pie- zoelectric material with an electroded surface. Such a wave does not exist in such layered structures in the absence of pie- zoelectricity. The mode shapes for displacement and electric potential in the piezoelectric layer are obtained and discussed theoretically. The study extends the regime of transverse sur- face waves and may lead to potential applications to surface acoustic wave devices.     

CONSERVATION LAWS IN FINITE MICROCRACKING BRITTLE SOLIDS

This paper addresses the conservation laws in finite brittle solids with microcracks. The discussion is limited to the 2-D cases. First, after considering the combination of the Pseudo-Traction Method and the indirect Boundary Element Method, a versatile method for solving multi-crack interacting problems in finite plane solids is proposed, by which the fracture parameters （SIF and path-independent integrals） can be calculated with a desirable accuracy. Second, with the aid of the method proposed, the roles the conservation laws play in the fracture analysis for finite microcracking solids are studied. It is concluded that the conservation laws do play important roles in not only the fracture analysis but also the analysis of damage and stability for the finite microcracking system. Finally, the physical interpretation of the M-integral is discussed further. An explicit relation between the M-integral and the crack face area, i.e., M = GS, has been discovered using the analytical method, which can shed some light on the Damage Mechanics issues from a different perspective.     

Piezoelectric generator based on torsional modes for power harvesting from angular vibrations

Torsional vibration of a circular piezoelectric shell of polarized ceramics mounted on a rotationally vibrating base is analyzed. The shell is properly electroded and connected to a circuit such that an electric output is generated. The structure analyzed represents a piezoelectric generator for converting mechanical energy from angular vibrations to electrical energy. Analytical expressions and numerical results for the output voltage, current, power, efficiency and power density are given.     

Nonlinearity analysis of piezoelectric micromachined ultrasonic transducers based on couple stress theory

This paper studies the static deformation behavior of a piezoelectric micromachined ultrasonic transducer (PMUT) actuated by a strong external electric field. The transducer membrane consists of a piezoelectric layer, a passive layer and two electrode layers. The nonlinearities of the piezoelectric layer caused by electrostriction under a strong electric field are analyzed. Because the thickness of the transducer membrane is on the microscale, the size dependence of the deformation behavior is evaluated using the couple stress theory. The results show that the optimal ratio of the top electrode diameter and the membrane diameter is around 0.674. It is also found that this optimal value does not depend on any other parameters if the thicknesses of the two electrodes are negligible compared with those of the piezoelectric and passive layers. In addition, the nonlinearities of the piezoelectric layer will become stronger along with the increase of the electric field, which means that softening of the membrane stiffness occurs when a strong external electric field is applied. Meanwhile, the optimal thickness ratio for the passive layer and the piezoelectric layer is not equal to 1.0 which is usually adopted by previous researchers. Because there exists size dependence of membrane deformation, the optimal value of this thickness ratio needs to be greater than 1.0 on the microscale.     

ELECTROELASTIC INTENSIFICATION NEAR ANTI-PLANE CRACK IN A FUNCTIONALLY GRADIENT PIEZOELECTRIC CERAMIC STRIP

Following the theory of linear piezoelectricity,we consider the electro-elastic prob-lems of a finite crack in a functionally gradient piezoelectric ceramic strip.By the use of Fouriertransforms we reduce the problem to solving two pairs of dual integral equations.The solution tothe dual integral equations is then expressed in terms of a Fredholm integral equation of the secondkind.Numerical calculations are carried out for piezoelectric ceramics.The electric field intensityfactors and the energy release rate are shown graphically,and the electroelastic interactions areillustrated.     

A continuum damage model for piezoelectric materials

In this paper, a constitutive model is proposed for piezoelectric material solids containing distributed cracks. The model is formulated in a framework of continuum damage mechanics using second rank tensors as internal variables. The Helrnhotlz free energy of piezoelectric mate- rials with damage is then expressed as a polynomial including the transformed strains, the electric field vector and the tensorial damage variables by using the integrity bases restricted by the initial orthotropic symmetry of the material. By using the Talreja＇s tensor valued internal state damage variables as well as the Helrnhotlz free energy of the piezoelectric material, the constitutive relations of piezoelectric materials with damage are derived. The model is applied to a special case of piezoelectric plate with transverse matrix cracks. With the Kirchhoff hypothesis of plate, the free vibration equations of the piezoelectric rectangular plate considering damage is established. By using Galerkin method, the equations are solved. Numerical results show the effect of the damage on the free vibration of the piezoelectric plate under the close-circuit condition, and the present results are compared with those of the three-dimensional theory.     

压电材料中的Bueckner功共轭积分及和J积分M积分的关系

研究横观各向同性压电材料中裂纹问题,提出了Bueckner功共轭积分在这类材料中的表达式：并通过引出两类辅助的应力-位移-电位移-电势场,证明功共轭积分和这类材料中的J积分和M积分仍然存在简单的两倍关系由此,各类在脆性材料断裂问题中已广泛应用的权函数方法可顺理成章地推广到压电材料的研究中来．这对独立地确定电位移强度因子和经典的I、II型应力强度因子提供了有力的数学上的工具．进而通过计算机械应变能释放率对压电材料中裂纹的稳定做出判断．     

压电材料中三维I型断裂力学分析

讨论了不可导通情况下三维横观各向刚性压电材料中受拉伸和电载荷作用的平片裂纹Ⅰ型断裂力学问题．使用自限部分概念，从二维线性压电理论出发，严格得到了一组以裂纹面位移间断和电势间断为未知变量的超奇异积分方程组；应用二维超奇异积分的主部分析法，从理论上分析得到了裂纹前沿应力和电势奇性指数以及应力和电位移奇性场，从而找到了以裂纹面位移间断和电势间断表示的应力和电位移强度因子、能量释放率表达式；为所得到的超奇异积分方程组建立了数值法，并用此计算了若干典型的平片裂纹问题，数值结果令人满意．     

The weight function for cracks in piezoelectrics

The weight function in fracture mechanics is the stress intensity factor at the tip of a crack in an elastic material due to a point load at an arbitrary location in the body containing the crack. For a piezoelectric material, this definition is extended to include the effect of point charges and the presence of an electric displacement intensity factor at the tip of the crack. Thus, the weight function permits the calculation of the crack tip intensity factors for an arbitrary distribution of applied loads and imposed electric charges. In this paper, the weight function for calculating the stress and electric displacement intensity factors for cracks in piezoelectric materials is formulated from Maxwell relationships among the energy release rate, the physical displacements and the electric potential as dependent variables and the applied loads and electric charges as independent variables. These Maxwell relationships arise as a result of an electric enthalpy for the body that can be formulated in terms of the applied loads and imposed electric charges. An electric enthalpy for a body containing an electrically impermeable crack can then be stated that accounts for the presence of loads and charges for a problem that has been solved previously plus the loads and charges associated with an unsolved problem for which the stress and electric displacement intensity factors are to be found. Differentiation of the electric enthalpy twice with respect to the applied loads (or imposed charges) and with respect to the crack length gives rise to Maxwell relationships for the derivative of the crack tip energy release rate with respect to the applied loads (or imposed charges) of the unsolved problem equal to the derivative of the physical displacements (or the electric potential) of the solved problem with respect to the crack length. The Irwin relationship for the crack tip energy release rate in terms of the crack tip intensity factors then allows the intensity factors for the unsolved problem to be formulated, thereby giving the desired weight function. The results are used to derive the weight function for an electrically impermeable Griffith crack in an infinite piezoelectric body, thereby giving the stress intensity factors and the electric displacement intensity factor due to a point load and a point charge anywhere in an infinite piezoelectric body. The use of the weight function to compute the electric displacement factor for an electrically permeable crack is then presented. Explicit results based on a previous analysis are given for a Griffith crack in an infinite body of PZT-5H poled orthogonally to the crack surfaces.     

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

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

Numerical solution of polarization saturation/dielectric breakdown model in 2D finite piezoelectric media

The polarization saturation (PS) model [Gao, H., Barnett, D.M., 1996. An invariance property of local energy release rates in a strip saturation model of piezoelectric fracture. Int. J. Fract. 79, R25–R29; Gao, H., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids 45, 491–510], and the dielectric breakdown (DB) model [Zhang, T.Y., Zhao, M.H., Cao, C.F., 2005. The strip dielectric breakdown model. Int. J. Fract. 132, 311–327] explain very well some experimental observations of fracture of piezoelectric ceramics. In this paper, the nonlinear hybrid extended displacement discontinuity-fundamental solution method (NLHEDD-FSM) is presented for numerical analysis of both the PS and DB models of two-dimensional (2D) finite piezoelectric media under impermeable and semi-permeable electric boundary conditions. In this NLHEDD-FSM, the solution is expressed approximately by a linear combination of fundamental solutions of the governing equations, which includes the extended point force fundamental solutions with sources placed at chosen points outside the domain of the problem under consideration, and the extended Crouch fundamental solutions with extended displacement discontinuities placed on the crack and the electric yielding zone. The coefficients of the fundamental solutions are determined by letting the approximated solution satisfy certain conditions on the boundary of the domain, on the crack face and the electric yielding zone. The zero electric displacement intensity factor in the PS model or the zero electric field strength intensity factor in the DB model at the outer tips of the electric yielding zone is used as a supplementary condition to determine the size of the electric yielding zone. Iteration approaches are adopted in the NLHEDD-FSM. The electric yielding zone is determined, and the extended intensity factors and the local J-integral are calculated for center cracks in piezoelectric strips. The effects of finite domain size, saturation property and different electric boundary conditions, as well as different models on the electric yielding zone and the local J-integral, are studied.     

Invariant formulation of the electromechanical enthalpy function of transversely isotropic piezoelectric materials

Summary Theoretical and numerical aspects of the formulation of electromechanically coupled, transversely isotropic solids are discussed within the framework of the invariant theory. The main goal is the representation of the governing constitutive equations for reversible material behaviour based on an anisotropic electromechanical enthalpy function, which automatically fulfills the requirements of material symmetry. The introduction of a preferred direction in the argument list of the enthalpy function allows the construction of isotropic tensor functions, which reflect the inherent geometrical and physical symmetries of the polarized medium. After presenting the general framework, we consider two important model problems within this setting: i) the linear piezoelectric solid; and ii) the nonlinear electrostriction. A parameter identification of the invariant- and the common coordinate-dependent formulation is performed for both cases. The tensor generators for the stresses, electric displacements and the moduli are derived in detail, and some representative numerical examples are presented.We thank Dipl.-Ing. H. Romanowski for his support and helpful remarks.     

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     

压电复合材料中正n边形孔边裂纹的反平面问题研究

本文基于Cauchy积分理论和Schwarz-Christoffel(SC)变换技术,针对压电复合材料中带一条裂纹的正n边形孔口缺陷的反平面断裂力学进行了探究．假设满足电不可通边界条件,利用Cauchy积分公式和留数定理,获得了任意正n边形裂尖处应力和电位移两个场强度因子以及全能量释放率的封闭形式的显式解．当正n边形边数取定时,所得解可退化为已有结果,以此验证方法的有效性．并通过数值算例,对比分析了n=3, n=4, n=5三种特殊情形对应的等效场强度因子和无量纲能量释放率与孔口边长、裂纹长度和受到的力、电载荷之间的曲线图．数值结果显示：正n边形孔洞的尺寸和裂纹长度均会促进裂纹扩展,且前者的影响更显著一些;正n边形边的数量增加会阻止裂纹的扩展;在电不可通边界条件下,机械载荷对裂纹的扩展始终贡献显著,电场对断裂行为的影响取决于机械载荷．本研究结果具有一般性,适用于任意正n边形孔边裂纹问题的求解,为压电复合材料元器件的优化设计和断裂特性分析提供了新思路．     

Thermo-electro-elastic transient responses in piezoelectric hollow structures

The paper presents an analytical method to solve thermo-electro-elastic transient response in piezoelectric hollow structures subjected to arbitrary thermal shock, sudden mechanical load and electric excitation. Volterra integral equation of the second kind caused by interaction between elastic deformation and electric field is solved by using an interpolation method. Thus, the exact expressions for the transient responses of displacement, stresses, electric displacement and electric potential in the piezoelectric hollow structures are obtained by means of Hankel transform, Laplace transform, and their inverse transforms. In Section 2, based on spherical coordinates, the governing equation of thermo-electro-elastic transient responses in a piezoelectric hollow sphere is found and the associated numerical results are carried out. In Section 3, based on cylindrical coordinates, the governing equation of thermo-electro-elastic transient responses in a non-homogeneous piezoelectric hollow cylinder is found and the corresponding numerical results are carried out. The results carried out may be used as a reference to solve other transient coupled problems of thermo-electro-elasticity in piezoelectric structures.