Analytical and numerical approaches to piezoelectric bimorph

We propose an efficient and accurate approach to piezoelectric bimorph based on a refined expansion of the elastic displacement and electric potential. The field approximation of the through-the-thickness variation accounts for a shear correction and a layerwise modelling for the electric potential. A particular attention is devoted to the boundary conditions on the bottom and top faces of the plate as well as to the interface continuity conditions for the electromechanical variables. The continuity condition on the electric potential imposes some restrictions on the approximation of the electric potential. Moreover, the continuity condition on the normal component of the electric induction at the bimorph interface is ensured by a Lagrange multiplier. The equations of the piezoelectric bimorph are obtained by using variational formulation involving the appropriate boundary and continuity conditions.A selection of numerical illustrations is presented for the series and parallel piezoelectric bimorphs simply supported under cylindrical bending conditions. Two types of electromechanical load are considered (i) a surface density of force applied on the top face and (ii) an electric potential applied on the bottom and top faces of the bimorph. The results thus obtained are compared to those provided by finite element computations performed for the full 3D model and by a simplified model without shear effect. At last, the problem of piezoelectric bimorph vibration is also examined for both closed and open circuit conditions. Excellent predictions with low error estimates of the local (profile) and global responses as well as resonant frequencies are observed. The comparisons assess of the effectiveness of the present approach to piezoelectric bimorph.     

Nonlinear dynamic response of piezoelastic laminated plates considering interfacial damage effects

This paper presents a nonlinear model for cross-ply piezoelastic laminated plates containing the damage effect of the intralayer materials and interlaminar interfaces. The model is based on the general six-degrees-of-freedom plate theory, the discontinuity of displacement, and electric potential on the interfaces are depicted by three shape functions, which are formulated according to solutions about three equilibrium equations and conservation of charge. By using the Hamilton variation principle, the three-dimensional nonlinear dynamic equations of piezoelastic laminated plates with damage are presented. Then using the finite difference method and the Newmark scheme, an analytical solution is presented. In numerical examples, the effects of different damage models, damage evolution, amplitude and frequency of electric loads on the nonlinear dynamic response of piezoelectric laminated plate with interfacial imperfections are investigated.     

带剪切型压电激励器的板的高阶剪切变形理论

利用压电材料固有的正，逆压电效应可以对结构变形和振动进行控制。与外加电场与极化方向平行于板厚度的压电材料的拉伸作动机制相比，外加电场与极化方向垂直的压电材料的剪切作动机制可以在作动器内产生较小的应力，从而降低作动器边界产生分层破坏的危险。本文对于压电材料的剪切作动机制进行研究，应用三阶剪切变形理论建立带剪切型压电激励器的智能层合板模型。采用哈密顿原理导出带剪切型压电激励器的层合板的控制方程。采用空间法得到了各种边界条件组合条件下板的解析解。数值算例对一三层板采用高阶和一阶剪切变形理论进行计算，结果表明两种理论所得的变形曲线很相似。但对于厚度剪切型激励器而言，由于激励器是引起板的剪切变形，而高阶剪切变形理论比一阶剪切变形理论能更好地反映结构的剪切应变能，因此高阶剪切变形理论可以提供板变形的更为精确的解。因此，对于厚度剪切型激励器，剪切变形理论的选取对于板变形结果的好坏有重要的作用。     

曲线加筋Kirchhoff-Mindlin板自由振动分析

相比传统加筋板,曲线加筋板能够更充分地发挥材料力学性能.在加筋板力学分析中,厚板通常采用Reissner-Mindlin理论,然而当板厚较薄时易出现剪切自锁,离散的Kirchhoff-Mindlin理论采用假设剪切应变场可避免该问题.针对曲线加筋Kirchhoff-Mindlin板自由振动分析,采用离散的Kirchhoff-Mindlin三角形单元和Timoshenko曲梁单元分别模拟板和加强筋,根据板的位移插值函数及筋板交界面的位移协调条件,建立基于板单元位移自由度的有限元方程.为了验证方法的有效性和准确性,采用直线加筋薄板、曲线加筋薄板和厚板3种模型进行算例研究,通过收敛性和精度分析来选择合理的有限元网格密度.直线加筋薄板前20阶固有频率均与文献结果吻合良好;曲线加筋板算例中,本文方法满足收敛条件的板单元数目为2469,Nastran模型板单元数目为6243;本文所得曲线加筋板固有频率与Nastran计算结果最大误差为3.4%.研究结果表明,本文方法无需筋板单元共节点,可使用较少的有限元网格数量,并能够保证计算精度;在离散Kirchhoff-Mindlin三角形板单元基础上构造Timoshenko梁单元可同时适用于曲线加筋薄板与厚板自由振动分析.     

Analytical solution for the cylindrical bending vibration of piezoelectric composite plates

An analytical solution for the cylindrical bending vibrations of linear piezoelectric laminated plates is obtained by extending the Stroh formalism to the generalized plane strain vibrations of piezoelectric materials. The laminated plate consists of homogeneous elastic or piezoelectric laminae of arbitrary thickness and width. Fourier basis functions for the mechanical displacements and electric potential that identically satisfy the equations of motion and the charge equation of electrostatics are used to solve boundary value problems via the superposition principle. The coefficients in the infinite series solution are determined from the boundary conditions at the edges and continuity conditions at the interfaces between laminae, which are satisfied in the sense of Fourier series. The formulation admits different boundary conditions at the edges of the laminated piezoelectric composite plate. Results for laminated elastic plates with either distributed or segmented piezoelectric actuators are presented for different sets of boundary conditions at the edges.     

A semi-analytical solution for static and dynamic analysis of plates with piezoelectric patches

A modified mixed variational principle for piezoelectric materials is established and the state-vector equation of piezoelectric plates is deduced directly from the principle. Then the exact solution of the state-vector equation is simply given, and based on the semi-analytical solution of the state-vector equation, a realistic mathematical model is proposed for static analysis of a hybrid laminate and dynamic analysis of a clamped aluminum plate with piezoelectric patches. Both the plate and patches are considered as two three-dimensional piezoelectric bodies, but the same linear quadrilateral element is used to discretize the plate and patches. This method accounts for the compatibility of generalized displacements and generalized stresses on the interface between the plate and patches, and the transverse shear deformation and the rotary inertia of the plate and patches are also considered in the global algebraic equation system. Meanwhile, there is no restriction on the thickness of plate and patches. The model can be also modified to achieve a semi-analytical solution for the transient responses to dynamic loadings and the vibration control of laminated plate with piezoelectric patches or piezoelectric stiffeners.     

Finite element analysis of the piezoelectric vibrations of quartz plate resonators with higher-order plate theory

A finite element formulation of the piezoelectric vibrations of quartz resonators based on Mindlin plate theory is derived. The higher-order plate theory is employed for the development of a collection of successively higher-order plate elements which can be effective for a broad frequency range including the fundamental and overtone modes of thickness-shear vibrations. The presence of electrodes is also considered for their mechanical effects.The mechanical displacements and electric potential are combined into a generalized displacement field, and the subsequent derivations are carried out with all the generalized equations. Through the standard finite element procedure, the vibration frequency, the vibration mode shapes and the electric potential distribution are obtained. The frequency spectra are compared with some well-known experimental results with good agreement.Our previous experience with finite element analysis of high-frequency quartz plate vibrations leads us to believe that memory and computing time will always remain as key issues despite the advances in computers. Hence, the use of sparse matrix techniques, efficient eigenvalue solvers, and other reduction procedures are explored.     

Static analysis of simply supported functionally graded and layered magneto-electro-elastic plates

In this article, static analysis of functionally graded, anisotropic and linear magneto-electro-elastic plates have been carried out by semi-analytical finite element method. A series solution is assumed in the plane of the plate and finite element procedure is adopted across the thickness of the plate such a way that the three-dimensional character of the solution is preserved. The finite element model is derived based on constitutive equation of piezomagnetic material accounting for coupling between elasticity, electric and magnetic effect. The present finite element is modeled with displacement components, electric potential and magnetic potential as nodal degree of freedom. The other fields are calculated by post-computation through constitutive equation. The functionally graded material is assumed to be exponential in the thickness direction. The numerical results obtained by the present model are in good agreement with available functionally graded three-dimensional exact benchmark solutions given by Pan and Han [Pan, E., Han, F., in press. Green’s function for transversely isotropic piezoelectric functionally graded multilayered half spaces. Int. J. Solids Struct.]. Numerical study includes the influence of the different exponential factor, magneto-electro-elastic properties and effect of mechanical and electric type of loading on induced magneto-electro-elastic fields. In addition further study has been carried out on non-homogeneous transversely isotropic FGM magneto-electro-elastic plate available in the literature [Chen, W.Q., Lee, K.Y., Ding, H.J., 2005. On free vibration of non-homogeneous transversely isotropic magneto-electro-elastic plates].     

Nonlinear transient isogeometric analysis of smart piezoelectric functionally graded material plates based on generalized shear deformation theory under thermo-electro-mechanical loads

We present a generalized shear deformation theory in combination with isogeometric (IGA) approach for nonlinear transient analysis of smart piezoelectric functionally graded material (FGM) plates. The nonlinear transient formulation for plates is formed in the total Lagrange approach based on the von Kármán strains, which includes thermo-piezoelectric effects, and solved by Newmark time integration scheme. The electric potential through the thickness of each piezoelectric layer is assumed to be linear. The material properties vary through the thickness of FGM according to the rule of mixture and the Mori–Tanaka schemes. Various numerical examples are presented to demonstrate the effectiveness of the proposed method.     

A two-dimensional closed-form solution for the free-vibrations analysis of piezoelectric sandwich plates

This work presents a two-dimensional (2D) closed-form solution for the free-vibrations analysis of simply-supported piezoelectric sandwich plates. It has the originality to consider all components of the electric field and displacement, thus satisfying exactly the electric equilibrium equation. Besides, the formulation considers full layerwise first-order shear deformation theory and through-thickness quadratic electric potential. Its independent mechanical and electric variables are decomposed using Fourier series expansions, then substituted in the derived mechanical and electric 2D equations of motion. The resulting eigenvalue system is then condensed so that only nine mechanical unknowns are retained. After its validation on single- and three-layer piezoelectric, and hybrid sandwich plates, the present approach was then used to analyze thickness modes of a square sandwich plate with piezoceramic faces and elastic cross-ply composite core. It was found that only the first three thickness modes are global, thus can be modeled by the mixed equivalent single-layer/layerwise approach, often retained in the literature; the remaining higher thickness modes being characteristic of sandwich behavior; i.e., dominated by the deformations of either the core or the faces. These results, together with presented through-thickness variations of the mechanical and electric variables clearly recommend full layerwise modeling. Several numerical results are provided for future reference for validation of 2D approximate analytical or numerical approaches; in particular, of 2D piezoelectric adaptive finite elements.     

安装蜂窝板动力学特性分析及主动控制试验研究

通过对蜂窝芯的等效化处理,建立了ANSYS的壳单元模型,并作有限元模态分析,然后与其试验模态分析结果比较,有限元分析结果和试验结果基本一致。通过ANSYS的PSD随机功率谱分析和模态分析所得振型确定了贴片位置,对蜂窝板进行了压电智能结构振动主动控制试验研究,得到了较好的振动抑制效果。分析结果为仪器安装蜂窝板的设计和实现智能结构控制提供了重要参考依据。     

A refined finite element method for bending analysis of laminated plates integrated with piezoelectric fiber-reinforced composite actuators

This research presents a finite element formulation based on four-variable refined plate theory for bending analysis of cross-ply and angle-ply laminated composite plates integrated with a piezoelectric fiber-reinforced composite actuator under electromechanical loading. The four-variable refined plate theory is a simple and efficient higher-order shear deformation theory, which predicts parabolic variation of transverse shear stresses across the plate thickness and satisfies zero traction conditions on the plate free surfaces. The weak form of governing equations is derived using the principle of minimum potential energy, and a 4-node non-conforming rectangular plate element with 8 degrees of freedom per node is introduced for discretizing the domain. Several benchmark problems are solved by the developed MATLAB code and the obtained results are compared with those from exact and other numerical solutions, showing good agreement.     

Finite element analysis of smart functionally graded plates

This paper deals with the derivation of a finite element model for the static analysis of functionally graded (FG) plates integrated with a layer of piezoelectric fiber reinforced composite (PFRC) material. The layer of PFRC material acts as the distributed actuator of the FG plates. The Young’s modulus of the FG plate is assumed to vary exponentially along the thickness of the plate while the Poisson’s ratio is assumed to be constant over the domain of the plate. The finite element model has been verified with the exact solutions for both thick and thin plates. Emphasis has been placed on investigating the effect of variation of piezoelectric fiber angle in the PFRC layer on its actuating capability of the FG plates. The finite element solutions also revealed that the activated PFRC layer is more effective in controlling the deformations of the FG plates when the layer is attached to the surface of the FG plate with minimum stiffness than when it is attached to the surface of the same with maximum stiffness.     

Topology optimization of piezoelectric nanostructures

We present an extended finite element formulation for piezoelectric nanobeams and nanoplates that is coupled with topology optimization to study the energy harvesting potential of piezoelectric nanostructures. The finite element model for the nanoplates is based on the Kirchoff plate model, with a linear through the thickness distribution of electric potential. Based on the topology optimization, the largest enhancements in energy harvesting are found for closed circuit boundary conditions, though significant gains are also found for open circuit boundary conditions. Most interestingly, our results demonstrate the competition between surface elasticity, which reduces the energy conversion efficiency, and surface piezoelectricity, which enhances the energy conversion efficiency, in governing the energy harvesting potential of piezoelectric nanostructures.     

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

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

基于高阶位移模式的压电层合板动力有限元模型

建立了含压电片层合板的有限元动力学模型。以位于压电层上下表面处的电场强度和层间电压为未知量,给出了三次函数的电势分布模式,采用Reddy的高阶剪切理论描述板的位移场,假设板厚度方向的正应力为零给出了减缩的本构方程,采用有限元方法,基于Hamilton原理导出结构的动力学方程,然后用静态缩聚的方法压缩掉电场自由度和次要的位移自由度。最后用四边形矩形单元求解了一对称铺层和非对称铺层悬臂板的固有频率,并与ANSYS结果对比验证了本文模型的精确性。     

Constitutive computational modelling foundation of piezoelectronic microstructures and application to high-frequency microchip DSAW resonators

This paper establishes a piezoelectric constitutive computational approach based on generalized eigenvalue and multivariable finite element solutions with potential applications to accurate and effective analysis of layered piezoelectric microstructures of arbitrary geometries and different anisotropic materials, to ease the limitation of current computer capacity in analyzing large-scale high-frequency disturbed surface acoustic waves (DSAW) by mounted electrodes in piezoelectric devices such as microchip SAW resonators. A new incompatible generalized hybrid/mixed element GQM5 is also proposed for improving predictions of the piezoelectric surface mount thermal stresses that are shear-dominated. The (generalized) plane strain constitutive model is numerically verified for piezoelectric finite element computation. With the help of computational piezoelectricity (electro-mechanics) for general layered structures with metal electrodes and anisotropic piezoelectric substrates, some new interesting, reliable and fundamental constitutive finite element results are obtained for high-frequency piezoelectric and mechanical SAW propagations and can be used for further applications. The ST-cut FEA results agree quite well with available exact and lab solutions for free surface case. The project supported by SRF for ROCS, SEM of China, the past Rutgers Univer-Seiko Epson Joint Fund and Zhejiang Provincial NSF     

压电复合材料层板柱面弯曲的精化理论

本研究旨在建立精确的压电复合材料层板理论。位移场和电势场采用近似表达，其沿板厚的分布通过构造高精度的位移分布函数和电势分布函数来描述。这两个函数由三雏弹性平衡方程和静电平衡方程的特解来导出，从而满足复杂的力电耦合关系和各类连续条件，保证了本文理论的高精度。本文理论仅涉及4个位移和电势变量，且不随层数的增加而增多，较之变量随层数而增多的分层理论简单得多，平衡方程形式简单；也便于发展成有限元等数值模型。通过与三维精确解比较，算例显示了本文理论的高精度和有效性。     

Hamiltonian parametric element and semi-analytical solution for smart laminated plates

Based on the Hellinger-Reissner (H-R) mixed variational principle for piezoelectric material, a unified 4-node Hamiltonian isoparametric element of anisotropy piezoelectric material is established. A new semi-analytical solution for the natural vibration of smart laminated plates and the transient response of the laminated cantilever with piezoelectric patch is presented. The major steps of mathematical model are as follows: the piezoelectric layer and host layer of laminated plate are considered as unattached three-dimensional bodies and discretized by the Hamiltonian isoparametric elements. The control equation of whole structure is derived by considering the compatibility of generalized displacements and generalized stresses on the interface between layers. There is no restriction for the side-face geometrical boundaries, the thickness and the number of layers of plate by the use of the present isoparametric element. Present method has wide application area.     

Vibration analysis of FG annular sector in moderately thick plates with two piezoelectric layers

Free vibration of functionally graded(FG) annular sector plates embedded with two piezoelectric layers is studied with a generalized differential quadrature(GDQ)method. Based on the first-order shear deformation(FSD) plate theory and Hamilton's principle with parameters satisfying Maxwell's electrostatics equation in the piezoelectric layers, governing equations of motion are developed. Both open and closed circuit(shortly connected) boundary conditions on the piezoelectric surfaces, which are respective conditions for sensors and actuators, are accounted for. It is observed that the open circuit condition gives higher natural frequencies than a shortly connected condition. For the simulation of the potential electric function in piezoelectric layers, a sinusoidal function in the transverse direction is considered. It is assumed that properties of the FG material(FGM) change continuously through the thickness according to a power distribution law.The fast rate convergence and accuracy of the GDQ method with a small number of grid points are demonstrated through some numerical examples. With various combinations of free, clamped, and simply supported boundary conditions, the effects of the thicknesses of piezoelectric layers and host plate, power law index of FGMs, and plate geometrical parameters(e.g., angle and radii of annular sector) on the in-plane and out-of-plane natural frequencies for different FG and piezoelectric materials are also studied. Results can be used to predict the behaviors of FG and piezoelectric materials in mechanical systems.