Abstract and Applied Analysis

Volume 2012 (2012), Article ID 290461, 14 pages

http://dx.doi.org/10.1155/2012/290461

## Structural-Electrical-Coupled Formulation for the Free Vibration of a Piezoelectric-Laminated Plate Using the Analytical Arbitrary Quadrilateral *p* Element

Department of Civil and Architectural Engineering, City University of Hong Kong, Hong Kong

Received 7 November 2011; Revised 13 February 2012; Accepted 14 February 2012

Academic Editor: Svatoslav Staněk

Copyright © 2012 Y. Y. Lee et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

An analytical quadrilateral *p* element is developed for solving the free vibrations of piezoelectric-laminated plates. The formulations of the displacement and strain fields are based on first-order shear deformation plate theory. The coupling effect between the electrical and stress fields is also considered. The Legendre orthogonal polynomials are used as the element interpolation functions, and the analytical integration technique is adopted. It is found that the present *p* element method gives high numerical precision results, fast and monotonic convergence rate. In the numerical cases, the effects of the number of hierarchical terms and mesh size on the convergence rate are investigated. Examples of square plates with different displacement and potential boundary conditions are studied. In the comparisons, the solutions of the present element are in good agreement with those obtained from other classical and finite element methods.

#### 1. Introduction

Various beam/plate problems have been hot research topics for decades (e.g., [1–15]), and many of them are solved using the finite element method (e.g., [16–21]). Beams/plates embedded with piezoelectric materials are widely used in various intelligent material systems for active control because of the electromechanical-coupling properties. For cases of simple geometric shapes, the two-dimensional (2D) and three-dimensional (3D) analytical methods that were developed by Benjeddou and Deü [22], Correia et al. [23], Heyliger and Saravanos [24], and Ding and Chen [25] are effective in analyzing the free vibrations of piezoelectric-laminated plates and provide clear physical meanings in the analytical and symbolical formulations. For cases of complex geometries, boundary conditions, and loadings, the finite element method (FEM) is more powerful and versatile than classical solutions. Most of the classical continuum solutions of plates have been limited to single- or two-mode approximations. This is due to the difficulties in obtaining the general multiple mode governing equations using Galerkin’s approach, especially for plates with complex boundary conditions. Hence, Saravanos et al., Lam et al., He et al., and Loja et al. [26–30] developed the finite-element and finite-strip methods for complicated piezoelectric-laminated plates. There are two main kinds of theories used for finite-element formulations. One is the classical plate theory (CPT), and the other one is the shear deformation theory, which branches out into first-order shear deformation theory (FSDT) and higher order shear deformation theory (HSDT). In general, theories that consider the shear deformation effect are more accurate for the analysis of piezoelectric-laminated plates. Thus, one such theory is used in this paper. The performance of an FEM model can be improved by a finer mesh (*h*-version) or refined integrated piezoelectric sensor and actuator interpolation (*p*-version). Zienkiewicz and Taylor [31] mentioned that the convergence rates of *p*-version elements would be more rapid than those of *h*-version elements for the same degrees of freedom (DOFs). The *p*-version elements that were developed by Houmat [32] and Woo et al. [33] were demonstrated to have fast monotonic convergence for cases of triangular plates and thick skew plates. According to the free-vibration analysis of skew Mindlin plates that Leung and Zhu [34] conducted, if the *p*-version approach is employed, then analytical integration is preferred in the procedures of setting up the element formulation; otherwise, the numerical integration errors may adversely affect the accuracy, and monotonic convergence of the natural frequencies cannot be guaranteed. In the axial free vibration analyses of beams, Leung and Chan [35] adopted trigonometric interpolation functions and the analytical integration technique to improve the accuracy. The analytical *p* element approach has been further applied to various vibrating plate problems [36–39] such as membrane vibration and thick plates and laminated plates. Although this kind of element can be applied to the analysis of plates with triangular and other complicated shapes, it is tedious to set up the mesh and assemble the elements when compared with quadrilateral elements.

Based on FSDT, the quadrilateral hierarchical element is developed for the free vibrations of piezoelectric-laminated plates in this paper. Legendre orthogonal polynomials are adopted in the interpolation functions to improve the accuracy. Comparisons with published results and case studies show the good accuracy and efficiency of the present quadrilateral *p* element.

#### 2. Theoretical Formulation

##### 2.1. Analytical Model

Figure 1 shows a laminated plate and the coordinate system defined at the midplane. The fiber direction is indicated by an angle , which is the positive-rotation angle of the principal material axes from the arbitrary axes. The elasticity modulus for a layer parallel to the fibers is , and perpendicular to the fibers is . The layers can be piezoelectric or of composite materials. It is assumed that the potential of the piezoelectric layer varies linearly through its thickness.

##### 2.2. Displacements and Strains

According to FSDT [40], the displacement fields of the plate are expressed as: where , and are the displacement components along the , , and axes, respectively, , , and are the associated midplane displacements, and , are the normal rotations about the and axes, respectively.

The linear strain-displacement relations are given by

The electrical potential is assumed to be linear through the thickness in each piezoelectric layer. Hence, the potential of the *k*th layer can be expressed as
where are the electrical potentials of the bottom and top surfaces of the *k*th layer, respectively, and and are the -coordinates corresponding to the bottom and top surfaces.

Therefore, the electric field in the layer can be written as:

In two-dimensional problems, the transverse displacement can be interpolated using the following equation (2.5): where ; and are Legendre orthogonal polynomials [40] (which are also used in the interpolations of displacements , , , ).

When , and 2, (2.5) represents a typical FEM interpolation. The hierarchical shape functions, when or , lead to zero displacement at the corner nodes. Additional DOFs appear along the four edges and in the interior of the element. The DOFs at the four corner nodes are represented by and , and the DOFs along the four edges are represented by or . Finally, the DOFs in the interior are represented by and .

##### 2.3. Constitutive Relationships

The constitutive relationships for the lamina-oriented arbitrarily, taking into account the piezoelectric effects are given as follows [25]: where is the stress matrix, is the strain vector, is the electrical displacement vector, is the electric field vector, is the transformed elastic stiffness coefficient matrix, the coefficients of which are explicitly given by Reddy [41], is the permittivity matrix, and is the transformed piezoelectric stress coefficient matrix, which is given by where are the piezoelectric stress coefficients in the directions parallel and perpendicular to the fibers, and is the fiber direction.

##### 2.4. Finite Element Model

The coordinate systems that are used to define an arbitrary quadrilateral plate element are shown in Figures 1 and 2 depict the Cartesian plane coordinate system of the element and the mapped square plane region. The Jacobian matrix is expressed in terms of the Cartesian coordinates at the four corner nodes: where , , , , , and .

The determinant of the Jacobian matrix is , and

The displacement fields , and and the electric potential are interpolated using the Legendre orthogonal polynomials as
where is the vector of generalized DOFs, is the nodal electric potential in the bottom surface of the *k*th layer, and is the identity matrix. Substituting the expressions of the displacement fields into (2.2) and (2.4) gives
where , is the total number of layers.

The potential energies of the element include the elastic strain energy, the piezoelectric energy, and the electrical energy, which are given as

The kinetic energy of the element is given by

The Hamilton principle is adopted in the derivation of the equation of the eigenvalue problem, which is as follows: where In the procedure of coordinate mapping, the derivative in and should be replaced by Note that analytic integration is adopted in (2.16). Hence, the analytic integrations of the stiffness and mass matrices are in the form of , where , , and are constants.

#### 3. Numerical Results and Discussions

In this section, case studies are carried out to show the performance of the quadrilateral *p* element. The elastic, piezoelectric, and dielectric material properties that are adopted in the case studies are given in Table 1. There are two sets of electric boundary conditions for the outer surface of the piezoelectric layers: a closed circuit condition , (i.e., the potential is grounded) and an open circuit condition , (i.e., the electric displacement is zero). Unless otherwise stated, the number of the hierarchical term *p* is 4, and the mesh in Figure 3(c) is used.

##### 3.1. Convergence Study of a Square Plate

The free vibration of a simply supported hybrid sandwich square plate is investigated here. The accuracy of the present element can be seen in this convergence study. The lamina configuration is PZT-4/GE/PZT-4. The thicknesses of the layers are 0.1 h, 0.8 h, 0.1 h, respectively. The close circuit condition is considered. The plate is meshed into one, two, or four quadrilateral elements, as shown in Figure 3. The fundamental frequency parameter is shown in Table 2 for different numbers of hierarchical terms. It can be observed that the monotonic convergence rate is very fast with respect to the number of hierarchical terms, and the results are in good agreement with those obtained from the 2D analytical method [22] and the 3D state space method [42]. It should be noted that the electrical potentials in [22, 42] are assumed to be quadratic through the plate thickness, whereas they are linear in the present method. Thus, there is a small difference between the results.

##### 3.2. Square Plates with Different Displacement Boundary Conditions

The first ten frequencies of a three-ply piezoelectric-laminated square plate (G-1195N/Ti-6A1-4V/G-1195N, 0.1 mm/0.48 mm/0.1 mm) are computed for three boundary conditions: simply supported (SSSS), fully clamped (CCCC), and cantilevered (CFFF). The electric boundary condition is a closed circuit. He et al. [28] derived finite element solutions using a traditional laminated plate element and CPT. A comparison of the results of our element and those of the traditional elements is shown in Table 3. Besides, a comparison for a static case is shown in Figure 4. The transverse normal stress, transverse displacement, and in-plane stress results of the sensor plate from the proposed method reasonably agree with those from [8].

##### 3.3. Square and Trapezoidal Plates with Different Potential Boundary Conditions

To study the validity of the present element for different potential boundary conditions, a five-ply simply supported plate, which is made of three plies of G/E with the laminate sequence (0°, 90°, 0°)and two surface-bonded PZT-4 piezoelectric layers, is considered. The thickness of both piezoelectric layers is 0.1 h, and all G/E layers have the same thickness. The unit density is adopted for fair comparison. The solutions of the present method are presented in Table 4 and are compared with the 3D exact solution. Besides, Figure 5 shows the trapezoidal plates with PZT-5H layers perfectly bonded on both the upper and lower surfaces having various boundary conditions (the number of Fourier terms used are ). The thicknesses of the trapezoidal plates and PZT-5H layers are 10 mm and 1 mm, respectively. Table 5 shows the natural frequencies of the aluminum trapezoidal plates with PZT-5H layers.

#### 4. Conclusion

A quadrilateral *p* element with analytical integration for the free vibration of piezoelectric-laminated composite plates is presented. The Legendre orthogonal polynomials are used as the element interpolation functions, and the analytical integration technique is adopted. The monotonic convergence rate of the present element is very fast with respect to the number of hierarchical terms. Comparisons between the solutions that are obtained from our method and those of other methods show that the element is accurate and efficient for free-vibration analyses of piezoelectric-laminated plates.

#### Acknowledgment

The research described herein was supported by Grant number CityU SRG7008100.

#### References

- A. Borisovich and J. Janczewska, “Stable and unstable bifurcation in the von Kármán problem for a circular plate,”
*Abstract and Applied Analysis*, vol. 2005, no. 8, pp. 889–899, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. L. Santos, J. Ferreira, and C. A. Raposo, “Existence and uniform decay for a nonlinear beam equation with nonlinearity of Kirchhoff type in domains with moving boundary,”
*Abstract and Applied Analysis*, vol. 2005, no. 8, pp. 901–919, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. A. Esquivel-Avila, “Dynamic analysis of a nonlinear Timoshenko equation,”
*Abstract and Applied Analysis*, Article ID 724815, 36 pages, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - Y. X. Li and H. Yang, “An existence and uniqueness result for a bending beam equation without growth restriction,”
*Abstract and Applied Analysis*, vol. 2010, Article ID 694590, 9 pages, 2010. View at Publisher · View at Google Scholar - M. Grasselli, V. Pata, and G. Prouse, “Longtime behavior of a viscoelastic Timoshenko beam,”
*Discrete and Continuous Dynamical Systems A*, vol. 10, no. 1-2, pp. 337–348, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. E. Muñoz Rivera and R. Racke, “Global stability for damped Timoshenko systems,”
*Discrete and Continuous Dynamical Systems A*, vol. 9, no. 6, pp. 1625–1639, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - E. Carrera, “An improved reissner-mindlin-type model for the electromechanical analysis of multilayered plates including piezo-layers,”
*Journal of Intelligent Material Systems and Structures*, vol. 8, no. 3, pp. 232–248, 1997. View at Publisher · View at Google Scholar · View at Scopus - E. Carrera and M. Boscolo, “Classical and mixed finite elements for static and dynamic analysis of piezoelectric plates,”
*International Journal for Numerical Methods in Engineering*, vol. 70, no. 10, pp. 1135–1181, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - E. Carrera, S. Brischetto, and P. Nali, “Variational statements and computational models for multifield problems and multilayered structures,”
*Mechanics of Advanced Materials and Structures*, vol. 15, no. 3-4, pp. 182–198, 2008. View at Publisher · View at Google Scholar · View at Scopus - E. Carrera and P. Nali, “Mixed piezoelectric plate elements with direct evaluation of transverse electric displacement,”
*International Journal for Numerical Methods in Engineering*, vol. 80, no. 4, pp. 403–424, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - E. Carrera and P. Nali, “Multilayered plate elements for the analysis of multifield problems,”
*Finite Elements in Analysis and Design*, vol. 46, no. 9, pp. 732–742, 2010. View at Publisher · View at Google Scholar - E. Carrera, S. Brischetto, and P. Nali,
*Plates and Shells for Smart Structures*, John Wiley and Sons, New York, NY, USA, 2011. - E. Carrera, S. Brischetto, and M. Cinefra, “Variable kinematics and advanced variational statements for free vibrations analysis of piezoelectric plates and shells,”
*Computer Modeling in Engineering and Sciences*, vol. 65, no. 3, pp. 259–341, 2010. View at Google Scholar · View at Zentralblatt MATH · View at Scopus - S. Brischetto and E. Carrera, “Coupled thermo-mechanical analysis of one-layered and multilayered plates,”
*Composite Structures*, vol. 92, no. 8, pp. 1793–1812, 2010. View at Publisher · View at Google Scholar · View at Scopus - E. Carrera, A. Büttner, and P. Nali, “Mixed elements for the analysis of anisotropic multilayered piezoelectric plates,”
*Journal of Intelligent Material Systems and Structures*, vol. 21, no. 7, pp. 701–717, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. Caboussat and R. Glowinski, “Numerical solution of a variational problem arising in stress analysis: the vector case,”
*Discrete and Continuous Dynamical Systems A*, vol. 27, no. 4, pp. 1447–1472, 2010. View at Publisher · View at Google Scholar - C. Thai-Hoang, N. Nguyen-Thanh, H. Nguyen-Xuan, and T. Rabczuk, “An alternative alpha finite element method with discrete shear gap technique for analysis of laminated composite plates,”
*Applied Mathematics and Computation*, vol. 217, no. 17, pp. 7324–7348, 2011. View at Publisher · View at Google Scholar - R. Sharma, R. Bhargava, and I. V. Singh, “Combined effect of magnetic field and heat absorption on unsteady free convection and heat transfer flow in a micropolar fluid past a semi-infinite moving plate with viscous dissipation using element free Galerkin method,”
*Applied Mathematics and Computation*, vol. 217, no. 1, pp. 308–321, 2010. View at Publisher · View at Google Scholar - M. Sajid, R. Mahmood, and T. Hayat, “Finite element solution for flow of a third grade fluid past a horizontal porous plate with partial slip,”
*Computers and Mathematics with Applications*, vol. 56, no. 5, pp. 1236–1244, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - Z. Wang and B. Hu, “Research of combined hybrid method applied in the Reissner-Mindlin plate model,”
*Applied Mathematics and Computation*, vol. 182, no. 1, pp. 49–66, 2006. View at Publisher · View at Google Scholar - C.-S. Chen, C.-P. Fung, and R.-D. Chien, “A further study on nonlinear vibration of initially stressed plates,”
*Applied Mathematics and Computation*, vol. 172, no. 1, pp. 349–367, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Benjeddou and J. F. Deü, “A two-dimensional closed-from solution for the free-vibrations analysis of piezoelectric sandwich plates,”
*International Journal of Solids and Structures*, vol. 39, no. 6, pp. 1463–1486, 2002. View at Publisher · View at Google Scholar · View at Scopus - V. M. Franco Correia, M. A. Aguiar Gomes, A. Suleman, C. M. Mota Soares, and C. A. Mota Soares, “Modelling and design of adaptive composite structures,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 185, no. 2–4, pp. 325–346, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - P. Heyliger and D. A. Saravanos, “Exact free-vibration analysis of laminated plates with embedded piezoelectric layers,”
*Journal of the Acoustical Society of America*, vol. 98, no. 3, pp. 1547–1557, 1995. View at Publisher · View at Google Scholar · View at Scopus - H. J. Ding and W. Q. Chen,
*Three Dimensional Problems of Piezoelectricity*, Nova Science Publishers, New York, NY, USA, 2001. - D. A. Saravanos, P. R. Heyliger, and D. A. Hopkins, “Layerwise mechanics and finite element for the dynamic analysis of piezoelectric composite plates,”
*International Journal of Solids and Structures*, vol. 34, no. 3, pp. 359–378, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - K. Y. Lam, X. Q. Peng, G. R. Liu, and J. N. Reddy, “A finite-element model for piezoelectric composite laminates,”
*Smart Materials and Structures*, vol. 6, no. 5, pp. 583–591, 1997. View at Publisher · View at Google Scholar · View at Scopus - X. Q. He, T. Y. Ng, S. Sivashanker, and K. M. Liew, “Active control of FGM plates with integrated piezoelectric sensors and actuators,”
*International Journal of Solids and Structures*, vol. 38, no. 9, pp. 1641–1655, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - M. A. Ramos Loja, C. M. Mota Soares, and C. A. Mota Soares, “Higher-order B-spline finite strip model for laminated adaptive structures,”
*Composite Structures*, vol. 52, no. 3-4, pp. 419–427, 2001. View at Publisher · View at Google Scholar · View at Scopus - M. A. R. Loja, J. Infante Barbosa, C. M. Mota Soares, and C. A. Mota Soares, “Analysis of piezolaminated plates by the spline finite strip method,”
*Computers and Structures*, vol. 79, no. 26–28, pp. 2321–2333, 2001. View at Publisher · View at Google Scholar · View at Scopus - O. C. Zienkiewicz and R. L. Taylor,
*The Finite Element Method*, McGraw-Hill, New York, NY, USA, 4th edition, 1989. - A. Houmat, “Triangular Fourier
*p*element for the analysis of membrane vibrations,”*Journal of Sound and Vibration*, vol. 230, no. 1, pp. 31–43, 2000. View at Publisher · View at Google Scholar · View at Scopus - K. S. Woo, C. H. Hong, P. K. Basu, and C. G. Seo, “Free vibration of skew Mindlin plates by p-version of F.E.M,”
*Journal of Sound and Vibration*, vol. 268, no. 4, pp. 637–656, 2003. View at Publisher · View at Google Scholar · View at Scopus - A. Y. T. Leung and B. Zhu, “Comments on ‘Free vibration of skew Mindlin plates by
*ρ*-version of F.E.M.’,”*Journal of Sound and Vibration*, vol. 278, no. 3, pp. 699–703, 2004. View at Publisher · View at Google Scholar · View at Scopus - A. Y. T. Leung and J. K. W. Chan, “Fourier
*p*element for the analysis of beams and plates,”*Journal of Sound and Vibration*, vol. 212, no. 1, pp. 179–185, 1998. View at Publisher · View at Google Scholar · View at Scopus - A. Y. T. Leung, B. Zhu, J. Zheng, and H. Yang, “Analytic trapezoidal Fourier
*p*element for vibrating plane problems,”*Journal of Sound and Vibration*, vol. 271, no. 1-2, pp. 67–81, 2004. View at Publisher · View at Google Scholar · View at Scopus - A. Y. T. Leung, B. Zhu, J. Zheng, and H. Yang, “A trapezoidal Fourier
*p*element for membrane vibrations,”*Thin-Walled Structures*, vol. 41, no. 5, pp. 479–491, 2003. View at Publisher · View at Google Scholar · View at Scopus - A. Y. T. Leung and B. Zhu, “Transverse vibration of thick polygonal plates using analytically integrated trapezoidal Fourier
*p*element,”*Computers and Structures*, vol. 82, no. 2-3, pp. 109–119, 2004. View at Publisher · View at Google Scholar · View at Scopus - A. Y. T. Leung, C. Xiao, B. Zhu, and S. Yuan, “Free vibration of laminated composite plates subjected to in-plane stresses using trapezoidal
*p*element,”*Composite Structures*, vol. 68, no. 2, pp. 167–175, 2005. View at Publisher · View at Google Scholar · View at Scopus - A. Houmat, “Hierarchical finite element analysis of the vibration of membranes,”
*Journal of Sound and Vibration*, vol. 201, no. 4, pp. 465–472, 1997. View at Publisher · View at Google Scholar · View at Scopus - J. N. Reddy,
*Mechanics of Laminated Composite Plates: Theory and Analysis*, CRC Press, New York, NY, USA, 1996. - H. J. Ding, W. Q. Chen, and R. Q. Xu, “New state space formulations for transversely isotropic piezoelasticity with application,”
*Mechanics Research Communications*, vol. 27, no. 3, pp. 319–326, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus