Internally Pressurized Radially Polarized Piezoelectric Cylinders

The piezoelectric phenomenon has been exploited in science and engineering for decades. Recent advances in smart structures technology have lead to a resurgence of interest in piezoelectricity, and in particular, in the solution of fundamental boundary-value problems. In this paper, we develop an analytic solution to the axisymmetric problem of an infinitely long, radially polarized, radially orthotropic piezoelectric hollow circular cylinder. The cylinder is subjected to uniform internal pressure, or a constant potential difference between its inner and outer surfaces, or both. An analytic solution to the governing equilibrium equations (a coupled system of second-order ordinary differential equations) is obtained. On application of the boundary conditions, the problem is reduced to solving a system of linear algebraic equations. The stress distributions in the cylinder are obtained numerically for two typical piezoceramics of technological interest, namely PZT-4 and BaTiO₃. It is shown that the hoop stresses in a cylinder composed of these materials can be made virtually uniform throughout the cross-section by applying an appropriate set of boundary conditions. This revised version was published online in June 2006 with corrections to the Cover Date.     

Anti-plane time-harmonic Green's functions for a circular inhomogeneity in piezoelectric medium with a spring- or membrane-type interface

The present work focuses on the two-dimensional anti-plane time-harmonic dynamic Green's functions for a circular inhomogeneity immersed in an infinite matrix with an imperfect interface, where both the inhomogeneity and matrix are assumed to be piezoelectric and transversely isotropic. Two types of imperfect interface, the spring-type interface with electromechanical coupling and the membrane-type interface, are considered. The former type is often used to model the electromechanical damage along the interface while the latter is largely employed to simulate surface/interface effect of nano-sized inhomogeneity. By using the Bessel function expansions, explicit solutions for the electromechanical fields induced by a time-harmonic anti-plane line force and line charge located in an unbounded matrix as well as the circular inhomogeneity are respectively derived. The present solutions can recover the anti-plane Green's functions for some special cases, such as the dynamic or quasi-static Green's functions of piezoelectricity with perfect interface as well as the dynamic or quasi-static Green's functions of pure elasticity with imperfect interface. For detailed discussions, selected analytical results are presented. Dependence of the electromechanical fields on circular frequency as well as interface properties is examined. The size effect related to interfacial imperfection is also discussed.     

Saint-Venant end effects for plane deformations of linear piezoelectric solids

The recent developments in smart structures technology have stimulated renewed interest in the fundamental theory and applications of linear piezoelectricity. In this paper, we investigate the decay of Saint-Venant end effects for plane deformations of a piezoelectric semi-infinite strip. First of all, we develop the theory of plane deformations for a general anisotropic linear piezoelectric solid. Just as in the mechanical case, not all linear homogeneous anisotropic piezoelectric cylindrical solids will sustain a non-trivial state of plane deformation. The governing system of four second-order partial differential equations for the two in-plane displacements and electric potential are overdetermined in general. Sufficient conditions on the elastic and piezoelectric constants are established that do allow for a state of plane deformation. The resulting traction boundary-value problem with prescribed surface charge is an oblique derivative boundary-value problem for a coupled elliptic system of three second-order partial differential equations. The special case of a piezoelectric material transversely isotropic about the poling axis is then considered. Thus the results are valid for the hexagonal crystal class 6mm. The geometry is then specialized to be a two-dimensional semi-infinite strip and the poling axis is the axis transverse to the longitudinal direction. We consider such a strip with sides traction-free, subject to zero surface charge and self-equilibrated conditions at the end and with tractions and surface charge assumed to decay to zero as the axial variable tends to infinity. A formulation of the problem in terms of an Airy-type stress function and an induction function is adopted. The governing partial differential equations are a coupled system of a fourth and third-order equation for these two functions. On seeking solutions that exponentially decay in the axial direction one obtains an eigenvalue problem for a coupled system of fourth and second-order ordinary differential equations. This problem is the piezoelectric analog of the well-known eigenvalue problem arising in the case of an anisotropic elastic strip. It is shown that the problem can be uncoupled to an eigenvalue problem for a single sixth-order ordinary differential equation with complex eigenvalues characterized as roots of transcendental equations governing symmetric and anti-symmetric deformations and electric fields. Assuming completeness of the eigenfunctions, the rate of decay of end effects is then given by the real part of the eigenvalue with smallest positive real part. Numerical results are given for PZT-5H, PZT-5, PZT-4 and Ceramic-B. It is shown that end effects for plane deformations of these piezoceramics penetrate further into the strip than their counterparts for purely elastic isotropic materials.     

Analysis of non-axisymmetric wave propagation in a homogeneous piezoelectric solid circular cylinder of transversely isotropic material

A study concerning the propagation of free non-axisymmetric waves in a homogeneous piezoelectric cylinder of transversely isotropic material with axial polarization is carried out on the basis of the linear theory of elasticity and linear electro-mechanical coupling. The solution of the three dimensional equations of motion and quasi-electrostatic equation is given in terms of seven mechanical and three electric potentials. The characteristic equations are obtained by the application of the mechanical and two types of electric boundary conditions at the surface of the piezoelectric cylinder. A novel method of displaying dispersion curves is described in the paper and the resulting dispersion curves are presented for propagating and evanescent waves for PZT-4 and PZT-7A piezoelectric ceramics for circumferential wave numbers m = 1, 2, and 3. It is observed that the dispersion curves are sensitive to the type of the imposed boundary conditions as well as to the measure of the electro-mechanical coupling of the material.     

The scattering of electroelastic waves by an ellipsoidal inclusion in piezoelectric medium

The scattering problem for a single ellipsoidal piezoelectric inclusion embedded in piezoelectric medium is investigated. Based on the polarization method, the extended displacements are expressed in terms of integral equations, whose kernels are obtained from the Green’s functions of homogenous matrix. In this paper, the 3D dynamic Green’s functions are derived by means of the Radon transform technique. To illustrate the use of the equations, scattering by a piezoelectric, ellipsoidal inhomogeneity in a piezoelectric medium is considered in the low frequency and an asymptotic formula for this scattering cross-section is obtained. Numerical results of the scattering cross-sections are carried out for a spheroidal BaTiO₃-inclusion in a PZT-5H-matrix.     

Stress singularities near the apex of a cylindrically polarized piezoelectric wedge

Summary In this paper, the stress singularities for a cylindrically polarized piezoelectric wedge are investigated. The characteristic equations are derived analytically by using the extended Lekhnitskii formulation. The piezoelectric material (PZT-4) is polarized in the radial, circular or axial direction, respectively. Similar to the rectilinearly polarized piezoelectric problem, the inplane and antiplane stress fields are uncoupled. The results show the variations of the singularity orders with the changes of the wedge angle, material constants, polarized direction, and the boundary conditions.     

The scattering of P-waves by a piezoelectric particle with FGPM interfacial layers in a polymer matrix

Propagation of P-wave in an unbounded elastic polymer medium which contains a set of nested concentric spherical piezoelectric inhomogeneities is formulated. The polymer matrix is made of Epoxy and is isotropic; each phase of the inhomogeneity is made of a different piezoelectric material and is radially polarized and has spherical isotropy. Note that the individual phases are homogeneous, and all interfaces are perfectly bonded. The scattered displacement and electric potentials in the matrix are expressed in terms of spherical wave vector functions and Legendre functions, respectively. The transmitted displacement and electric potentials within each phase of the piezoelectric particle are expressed in terms of Legendre functions. The equations of motion and electrostatics in each phase of the piezoelectric inhomogeneity lead to a system of coupled second order differential equations, which is solved using the generalized Frobenius series. The present theory is extended to the case where the core of the inhomogeneity is made of PZT-4 and its coating is made of functionally graded piezoelectric material (FGPM) whose microstructural composition varies smoothly from PZT-4 at the core–coating interface to Epoxy at the coating–matrix interface. The effects of different types of variation in the electro-mechanical properties of FGPM on scattering cross-section and other electro-mechanical fields are addressed. The present theory is valid for arbitrary coating thickness, and arbitrary frequencies.     

Green’s functions for transversely isotropic piezoelectric functionally graded multilayered half spaces

By virtue of two systems of vector functions and the propagator matrix method, Green’s functions for transversely isotropic, piezoelectric functionally graded (exponentially graded in the vertical direction), and multilayered half spaces are derived. It is observed that the homogeneous solution and propagator matrices for each functionally graded layer in the transformed domain are independent of the choice of the two systems of vector functions. For a point force and point charge density applied at any location of the functionally graded half space, Green’s functions are expressed in terms of one-dimensional infinite integrals. To carry out the numerical integral involving Bessel functions, an adaptive Gauss quadrature approach is introduced and modified. The piezoelectric functionally graded Green’s functions include those in the corresponding elastic functionally graded media as special results with the latter being also unavailable in the literature. Two piezoelectric functionally graded half-space models are analyzed numerically: one is a functionally graded PZT-4 half space, and the other a coated functionally graded PZT-4 layer over a homogeneous BaTiO₃ half space. The effects of different exponential factors on Green’s function components are clearly demonstrated, which could be useful in the design and manufacturing of piezoelectric functionally graded structures.     

An analytical study of vibration in functionally graded piezoelectric nanoplates: nonlocal strain gradient theory

In this paper, we analytically study vibration of functionally graded piezoelectric(FGP) nanoplates based on the nonlocal strain gradient theory. The top and bottom surfaces of the nanoplate are made of PZT-5 H and PZT-4, respectively. We employ Hamilton's principle and derive the governing differential equations. Then, we use Navier's solution to obtain the natural frequencies of the FGP nanoplate. In the first step, we compare our results with the obtained results for the piezoelectric nanoplates in the previous studies. In the second step, we neglect the piezoelectric effect and compare our results with those obtained for the functionally graded(FG) nanoplates. Finally, the effects of the FG power index, the nonlocal parameter, the aspect ratio, and the lengthto-thickness ratio, and the nanoplate shape on natural frequencies are investigated.     

Magnetoelectric Green's functions and their application to the inclusion and inhomogeneity problems

Explicit expressions of magnetoelectric Green's functions are obtained for a transversely isotropic medium exhibiting coupling between the static electric and magnetic fields utilizing the contour integral representation. Four Green's functions exist which represent the coupled static electric and magnetic response to a unit point electric or magnetic charge. The Green's functions are applied to analyze the inclusion and inhomogeneity problems in an infinite magnetoelectric medium, and explicit, closed form expressions are obtained for the Eshelby type tensors. The magnetoelectric Eshelby's tensors can be readily used in the solution of numerous problems in the mechanics and physics of magnetoelectric solids.     

Penny shaped crack in a piezoelectric cylinder surrounded by an elastic medium subjected to combined in-plane mechanical and electrical loads

The fracture problem of a penny shaped crack in a piezoelectric ceramic cylinder surrounded by an infinite elastic medium under in-plane normal mechanical and electrical loads is considered with the electric continuous boundary conditions on the crack surface. By using the potential theory and Hankel transform, a system of dual integral equations is obtained, and expressed to a Fredholm integral equation of the second kind. The mechanical and electrical field equations and all sorts of field intensity factors of mode I are obtained, and the numerical values of various field intensity factors for PZT-6B piezoelectric ceramic surrounded by several different elastic media are graphically shown for a uniform load and a ring-shaped load, respectively. And the effects of the size of the piezoelectric cylinder and the elastic material properties on various field intensity factors are obtained.     

考虑挠曲电效应的压电纳米薄板力-电-热耦合特性研究

摘要：挠曲电效应是由应变梯度引起的,与尺度相关的力电耦合效应。基于Kirchhoff板假设和挠曲电理论,本文推导了温度和电压作用下的压电薄板力-电-热耦合微分控制方程,定量分析了微分控制方程中非线性项的影响,并针对四周固支压电薄板采用Ritz法求解,数值计算了压电薄板的弯曲和振动行为。在研究温度和挠曲电效应对薄板耦合特性和力学行为的影响时,本文分别考虑了材料系数不随温度变化和随温度线性变化两种情况。以PZT-5H为例,我们讨论了挠曲电和温度对压电薄板的横向位移和固有频率的影响。研究结果表明挠曲电效应对压电纳米薄板的力学行为影响很大,且具有明显的尺寸效应。此外,薄板对温度变化非常敏感。因此,可通过挠曲电效应和温度来调控压电纳米薄板的多场耦合特性和力学行为,进而优化基于压电薄板的NEMS/MEMS中传感器、作动器等电子器件的性能。     

考虑挠曲电效应的压电纳米薄板力-电-热耦合特性研究

摘要：挠曲电效应是由应变梯度引起的,与尺度相关的力电耦合效应。基于Kirchhoff板假设和挠曲电理论,本文推导了温度和电压作用下的压电薄板力-电-热耦合微分控制方程,定量分析了微分控制方程中非线性项的影响,并针对四周固支压电薄板采用Ritz法求解,数值计算了压电薄板的弯曲和振动行为。在研究温度和挠曲电效应对薄板耦合特性和力学行为的影响时,本文分别考虑了材料系数不随温度变化和随温度线性变化两种情况。以PZT-5H为例,我们讨论了挠曲电和温度对压电薄板的横向位移和固有频率的影响。研究结果表明挠曲电效应对压电纳米薄板的力学行为影响很大,且具有明显的尺寸效应。此外,薄板对温度变化非常敏感。因此,可通过挠曲电效应和温度来调控压电纳米薄板的多场耦合特性和力学行为,进而优化基于压电薄板的NEMS/MEMS中传感器、作动器等电子器件的性能。     

Electroelastic field of a piezoelectric annular finite cylinder

This paper presents a theoretical study of a piezoelectric annular cylinder under axisymmteric electromechanical loading. The piezoelectric material is assumed to be transversely isotropic and the general solutions of the governing equations are obtained in terms of a Fourier–Bessel series containing Bessel functions of the first and second kind. The boundary-value problems for vertical pressure and an electric charge loading applied to the ends of an annular cylinder are solved by expanding the applied loading in terms of a Fourier–Bessel series. Selected numerical results for the electroelastic field of an annular cylinder are presented for different aspect ratios of a cylinder and material properties.     

On the Green's Functions and Boundary Integral Formulation of Elastic Planes with Cutouts

ABSTRACT

The problem of an infinite elastic plane that contains a hole of arbitrary shape and is subjected to a concentrated unit load is considered. The Green's function (influence function) for the problem is formulated by means of two complex potential functions. This is accomplished by mapping the region that is exterior to the hole onto a unit circle. A class of closed contour hole shapes is analyzed. Green's functions for an elliptical hole and a class of triangular holes are determined. Green's functions for a class of rectangular holes are also discussed. In order to determine stress and displacement fields for the finite plane problem, Green's function is employed and an indirect boundary integral equation is formulated, with the integrand of the integral equation incorporating the effect of the hole. The contour of the hole is no longer considered a part of the boundary and only the contour of the region that is exterior to the hole is subdivided into boundary elements. Examples for elliptical and triangular holes are solved.     

考虑几何与压电非线性行为压电层合轴对称圆板的精确解

研究可移简支及夹支边界条件下,轴对称压电层合圆板在强电场和机械荷载联合作用下的非线性变形.考虑电致伸缩的非线性压电效应及几何非线性的影响,导出轴对称压电层合圆板的控制方程.通过调整坐标轴的位置对控制方程进行简化,得到关于挠度和径向力的4阶非线性控制方程.再通过简单的积分并引入无量刚变量将控制方程等价地化为2阶非线性耦合微分方程组.利用幂级数法得到可移简支及夹支边界条件下强电场和均布荷载共同作用时的挠度、径向力及径向位移的幂级数精确解.通过对双、单压电晶片执行器的数值计算及分析,得到电场、外载对于位移、径向力的影响关系.     

Penny-shaped crack in a piezoelectric cylinder under electro-mechanical loads

Summary The fracture behavior of a penny-shaped crack in a axisymmetrical piezoelectric ceramic cylinder of finite radius under mechanical and electrical loads is analyzed under electric continuous boundary conditions on the crack surface. The potential theory and Hankel transform are used to obtain a system of dual integral equations, which is then expressed as a Fredholm integral equation. Singular mechanical and electrical fields and field intensity factors of mode I are obtained. The numerical values of various field intensity factors for PZT-6B piezoelectric ceramics are graphically shown for a uniform load and a ring-shaped load, respectively. The effects of the radius of the cylinder on the field intensity factors are investigated. This work was supported by the Korea Research Foundation Grant (KRF-2001-041-E00057).     

轴对称横观各向同性热弹性圆柱的分解

为了得到精确的应力场、位移场、温度场,将扭转圆轴的精化理论研究方法推广到轴对称横观各向同性热弹性圆柱。利用Bessel函数以及轴对称横观各向同性热弹性圆柱的通解,给出了轴对称横观各向同性热弹性圆柱的分解定理。根据柱面齐次边界条件获得了精确的精化方程,精化方程可以分解为一阶方程、超越方程、温度方程,从而将横观各向同性热弹性圆柱的轴对称问题分解为轴向拉压问题、超越问题、热-应力耦合问题。超越部分对应端部自平衡情况,可以清晰地了解到端部应力分布对内部应力场的影响,热-应力耦合部分对应无外加应力场时圆柱内部因温度变化引起的热应力。     

Dynamic behavior of piezoelectric/magnetostrictive composite hollow cylinder

The dynamic behavior of a multilayered, perfectly bonded piezoelectric/magnetostrictive composite hollow cylinder under radial deformation is investigated. The superposition method, the state space method as well as the separation of variables method are elegantly integrated in the solution approach. The governing equations are finally transformed into two Volterra integral equations of the second kind with respect to two functions of time. The elastic, electric and magnetic fields are finally obtained according to solving the integral equations. Free vibrations and transient responses are demonstrated by numerical experiments.     

压电空心圆柱中波的传播

应用三维压电弹性体轴对称模型对压电空心圆柱中波的传播进行了研究. 发现在 圆柱中是否具有压电性质会对波的传播带来显著的差异. 当波长趋向于零时，在压电圆柱中 拟P波的波速渐进趋向于横观各向同性弹性体的准P波波速，而非压电圆柱中拟P波的波速 渐进趋向于一维杆模型中的P波波速；在压电圆柱中拟SV波存在驻波现象. 圆柱中的SH波 同电场无关，所以在压电圆柱和非压电圆柱中SH波具有相同的频散曲线. 应用积分变 换方法将圆柱的控制方程同其侧边界条件相结合，得到了一组动力学方程. 针对具体的侧边 界条件，得到相对应的波导条件和频散方程. 在此基础上通过数值计算模拟了圆柱受到端部 应力脉冲激励后的瞬态响应. 同时讨论了空心圆柱的半径比对轴向波传播的影响.