Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 425317, 9 pages

http://dx.doi.org/10.1155/2014/425317

Research Article

An Efficient Spectral Element Model with Electric DOFs for the Static and Dynamic Analysis of a Piezoelectric Bimorph

Xingjian Dong, Zhike Peng, Wenming Zhang, HongXing Hua, and Guang Meng

Institute of Vibration Shock & Noise, State Key Laboratory of Mechanical System and Vibration, Shanghai Jiao Tong University, Shanghai 200240, China

Received 12 February 2014; Accepted 26 March 2014; Published 24 April 2014

Academic Editor: Weichao Sun

Copyright © 2014 Xingjian Dong 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 efficient spectral element (SE) model for static and dynamic analysis of a piezoelectric bimorph is proposed. It combines an equivalent single layer (ESL) model for the mechanical displacement field with a sublayer approximation for the electric potential. The 2D Gauss-Lobatto-Legendre (GLL) shape functions are used to discretize the displacements and then the governing equation of motion is derived following the standard SE method procedure. It is shown numerically that the present SE model can well predict both the global and local responses such as mechanical displacements, natural frequencies, and the electric potentials across the bimorph thickness. In the case of bimorph sensor application, it is revealed that the distribution of the induced electric potential across the thickness does not affect the global natural frequencies much. Furthermore, the effects of the order of Legendre polynomial and the mesh size on the convergence rate are investigated. Comparison of the present results for a bimorph sensor with those from 3D finite element (FE) simulations establishes that the present SE model is accurate, robust, and computationally efficient.

1. Introduction

Piezoelectric materials, especially lead zirconate titanate (PZT), can function either as actuator or as sensor for their inherent coupling electromechanical character. Due to its wide ranging applications in electroacoustic transducers, medical devices, microrobot, and atomic force microscope (AFM) cantilevers, the study of smart structures consisting of PZT sensors and PZT actuators has drawn considerable attention of many researchers in the fields of active vibration control, noise attenuation, and damage detection. The most popular simple PZT sensor or actuator consists usually of a slab of piezoelectric ceramic such that the PZT layer expands or contracts mainly in its length direction. However, the motion of a single layer, as well as the induced electric potential, is extremely small. To achieve practically meaningful actuation or sensing capabilities in PZT devices, a piezoelectric bimorph consisting of two PZT layers is commonly used for the reason that it can produce flexural deformation significantly larger than the length or thickness deformation of the individual layers [1–6]. However, the applications of the piezoelectric bimorph require the development of admissible approaches entailing capabilities to predict the global responses of the bimorph structure, such as the deflection and natural frequencies. Additionally, the approaches should address the local responses, such as the through-the-thickness variation of the electric potentials.

In the past three decades, a large variety of models have been developed to predict the static and dynamic responses of piezoelectric bimorph structures under all kinds of electromechanical loads with emphasis on approximating the mechanical displacement and electric potential correctly. The classification of the various models is mainly based on the kinematic assumption for approximating the through-the-thickness variation of the electromechanical state variables and representation method of the piezoelectric layers [1]. The accurate responses of the piezoelectric bimorph structures can be obtained by solving the 3-dimensional (3D) coupled field equations with exact satisfaction of the mechanical and electric boundary conditions [7]. However, the 3D analytical solutions do not provide the results for more general case of loading and complicated boundary conditions. Thus, there is a genuine requirement of improved numerical techniques such as finite element (FE) method. The 3D FE modeling of piezoelectric bimorph structures, which is now available in many commercial FE software programs such as ABAQUS and ANSYS, will result in systems with a large number of degrees of freedom (DOFs). To overcome the computational inefficiency associated with the 3D FE models, the equivalent single layer (ESL) model applied to both mechanical and electrical unknowns has been proposed. There are two main kinds of theories used for ESL model. One is the classical laminated plate theory (CLPT) [8, 9], and the other is the shear deformation theory, which branches out into first-order shear deformation theory (FSDT) [10, 11] and higher-order shear deformation theory (HSDT) [12, 13]. The ESL model is simple and capable of predicting the global responses of the bimorph, but it does not account for the nonlinear distribution of the electric potential across piezoelectric layers as observed from 3D solutions [4]. This shortcoming inspired the researchers to develop layer-wise theory [4, 14–16] or the sublayer theory [2, 7, 17–19] for approximation of the electric potential. In the latter case, the piezoelectric layer is divided into appropriate number of sublayers along the thickness direction and a linear through-the-thickness electric potential distribution for each sublayer is assumed. It is expected that the actual distribution of the electric potential across the bimorph thickness can be approached with more sublayers adopted. Therefore, a sensible strategy is to use an ESL model for the mechanical variables and a layer-wise theory or a sublayer theory for the electric variables, respectively [20–22].

The FE model which considers a global ESL approximation for the mechanical field variables and a layer-wise theory or a sublayer theory for the electric potential will inevitably result in a large number of DOFs for practical dynamical problems. Recently, the spectral element (SE) method, which combines the geometric flexibility of the FE method with the accuracy of the global pseudospectral method, is widely used in many research areas related to seismology, fluid dynamics, and acoustics [23–25]. More recently, the SE method was used to simulate wave propagations for the purpose of damage detection in structures [26]. In fact, the SE method and FE method are closely related and built on the same ideas. The main difference between them is that SE method uses orthogonal polynomials, such as Legendre and Chebyshev polynomials, in the approximation functions; therefore, the mass matrix is diagonal, which is a very significant advantage over conventional FE method especially for dynamic analysis. Moreover, to have an accurate simulation with conventional FE method, a mesh with large number of nodes and elements is inevitably needed. The SE method, in which the polynomial order is increased and the mesh size is decreased, can be used to overcome this problem. However, it appears that the use of SE method for problems of static and dynamic analysis of piezoelectric bimorph has not been widely reported in the literatures so far.

The objective of the present work is to develop an efficient and accurate electromechanical coupled 2D SE model for static and dynamic analysis of a piezoelectric bimorph using combination of ESL for the displacements and a sublayer model for the electric potentials. The remainder of this paper is outlined as follows. The approximations for the displacements and the electric potentials are given in Section 2; the SE model with electric DOFs is then derived utilizing Legendre orthogonal polynomials in the interpolation function. In Section 3, the global responses, such as the deflection and natural frequencies, and the local responses, such as the electric potential across the thickness, are presented. The results are also compared to those provided by the 3D FE simulations. The convergence of the present SE model and the robustness of algorithm with respect to the order of the Legendre polynomial as well as mesh size are then investigated. At last the closing remarks and discussion of the results are given in Section 4.

2. Mathematical Formulation

2.1. Constitutive Relationships, Displacement, and Strain

A piezoelectric bimorph made of two PZT-4 piezoelectric layers (as shown in Figure 1) undergoing a surface density of normal force and electric potential applied to the top and bottom surfaces is considered in this work. The length and width of the bimorph are 25 mm and 12.5 mm, respectively. Both piezoelectric layers have the same thickness and are poled in the -direction. The reference plane is chosen to be the middle plane of the bimorph, and the-axis is defined as the direction normal to the plane according to the right-hand rule. Generally, the linear constitutive equations of piezoelectric materials, including the converse and direct piezoelectric effects, can be written as

where and represent stress vector and strain vector, respectively. is the electric field vector, is the electric displacement vector, is the elastic coefficient matrix, is the dielectric coefficient matrix, and is the piezoelectric stress coefficient matrix.

Figure 1: Geometry of a piezoelectric bimorph.

The FSDT is based on the constant transverse shear strain assumption, which leads to the displacement field [27]

where , , and are the displacements in the -, -, and -directions, respectively. , , and are the in-plane and transverse displacements of a point on the middle plane, respectively. and - denote the rotations of a transverse normal about the - and -axes, respectively. We define

where is the displacement vector and is a generalized displacement vector. Then (2) can be expressed in matrix form as

where

The relationship between the strains and the displacements can be written as follows:

where is the derivation operator defined as

2.2. Spectral Element Discretization of a Plate

As in the classical FE method, the bimorph is firstly meshed to a number of nonoverlapping rectangular elements. Each rectangular element, denoted by , is then mapped to a reference element, denoted by , using an invertible local mapping. The discretization procedure is illustrated in Figure 2. Subsequently, a set of nodes, denoted by , are defined in the local coordinate system of the reference element as roots of the following polynomial expression:

where is the th order Legendre polynomial. In fact, the nodes are the 2D Gauss-Lobatto-Legendre (GLL) points. In this way the nodes of the element can be specified in the local coordinate system of the element, as shown in Figure 2. The 2D shape function built on the specified node can be written as

where is the 1D shape function of order at the 1D GLL points which can be defined as [26]

The 2D shape functions are used for interpolating both the element coordinates and the element displacements. Consequently, coordinates and within each may be uniquely related to and upon the invertible mapping

where and denote the coordinates of and , respectively, of the element nodes . The generalized displacements , , , , and over a reference element are discretized by the 2D shape functions as

where , , , , and are the nodal values of the generalized displacements. The discrete element nodal displacement vector is expressed as

where is the displacement vector of the node

Figure 2: Discretization of a plate and an example of spectral element.

Substituting (12) into (4) yields

where is the matrix of displacement shape function which can be written as

with

where is a identity matrix. Substituting (15) into (6) yields

where is strain-displacement matrix which can be written as

2.3. A Sublayer Model for the Electric Potentials

In order to accurately model the through-the-thickness distribution of the electric potential, each layer of the piezoelectric bimorph is subdivided mathematically into thinner sublayers. As shown in Figure 1, the sublayers are numbered in top-to-bottom order. The -coordinates of the top and bottom surfaces of the th sublayer are denoted by and , respectively. It is assumed that in each sublayer the electric potential has a linear variation across the thickness such that

where is the shape function of the electric potential and is a column matrix composed of the electric potentials at the top and the bottom surfaces of the th sublayer, which can be expressed as

In this way the electric potential is approximated as piecewise linear across the thickness and it is expected that the actual nonlinear electric potential field of the piezoelectric bimorph can be approached with an appropriate number of sublayers. Considering that the piezoelectric bimorph is discretized using 2D mesh, each sublayer is also discretized using the same mesh to keep the compatibility. Consequently, an element potential vector which is of the following form is then introduced in the spectral plate finite element:

The surface potential of the sublayer, , is assumed to be constant over the element and are treated as elemental degrees of freedom (DOFs), as illustrated in Figure 2. Furthermore, the top and bottom surfaces of the piezoelectric layers are always coated with metallic coatings of zero thickness and the potentials on the electrodes should be taken as independent of ,. Thus the present approach combines an ESL theory for the displacement field and a piecewise linear approximation for the electric potential. Under the quasi electrostatic approximation, the electric field and the electric potential in each sublayer have the following relationship:

where is the electric field of the th sublayer and is the electric field-potential matrix, given by

2.4. Governing Equations

Hamilton’s variational principle is adopted in the derivation of the elementary governing equation of motion

where denotes the elementary mass matrix; is mechanical stiffness matrix; and are the piezoelectric coupling matrices; is the dielectric permittivity matrix; and denote the nodal external force vector and the nodal externally applied charge vector, respectively:

where is the mass density, is the surface force vector, and is the surface charge density vector. is the Jacobian matrix of the mapping (11) which is defined by

The GLL integration rule is then used to calculate the characteristic matrices and the nodal force vector in (25) at the elemental level [26].

In this study, the interface between the two PZT layers is grounded. The bimorph is supposed to be used as a sensor in closed circuit, in which the top and bottom surfaces are grounded and a uniform pressure load is applied to the upper surface. By applying the electric boundary conditions, the DOFs for the electric potentials are condensed out such that (25) is finally of the form

where is the mechanical stiffness matrix induced by the electromechanical coupling of the piezoelectric materials and is the mechanical forces induced by the applied voltages of piezoelectric actuators [2]. The electric potential is then recovered by the inverse process of the aforementioned condensation. Assembling all elementary equations, one can have a global dynamic system equation

where , , , , and are the assembled counterparts of matrices , , , , and ; is the global nodal displacement vector. Since the electric potential DOFs for the sublayers have been condensed out, this approach will not result in an excessive number of potential field variables. For the purpose of static analysis, the governing equations of motion in (29) reduces to

For the dynamic frequency analysis problems, the corresponding eigenvalue problem is

where is the natural frequency.

3. Numerical Results and Discussion

In this section, case studies are carried out to demonstrate the efficiency and accuracy of the present model in estimating both the global responses, such as the deflection and natural frequencies, and the local responses, such as the through-the-thickness variation of the electric potential. A simply supported rectangular piezoelectric bimorph shown in Figure 1, which was analyzed by Fernandes and Pouget [1], is considered here. The material constants of PZT-4 are given as

Different slenderness ratios such as , , and , which represent the thick, moderately thick, and thin bimorph plate, respectively, are considered. Unless otherwise stated, the order of Legendre polynomial is chosen as , and the mesh in Figure 2 is used in this work.

3.1. Static Responses of the Bimorph Sensor

The bimorph is supposed to be used as a sensor in closed circuit. To achieve practically meaningful sensor capabilities and guarantee that the piezoelectric material behaves linearly, a uniform pressure load of is applied to the top surface. Equation (30) is utilized to calculate the static responses of the bimorph sensor. The numerical results for the deflection and the electric potential are given with the following dimensionless units:

where is the stiffness constant of elastic coefficient matrix shown in (32) and the amplification factor is taken as . For the purpose of comparison, a coupled 3D analysis is carried out using 20-noded hexahedral 3D piezoelectric elements (C3D20RE) with a mesh size of in ABAQUS and the results from the coupled 3D FE analysis are taken as accurate.

The through-the-thickness variations of both the deflection and the electric potential at the centre of the bimorph are collected in Figure 3 for the slenderness ratio . It can be observed from Figure 3(a) that the deflection estimated by the present method adopting different number of sublayers is constant through the thickness and it is a good approximation of the nonlinear distribution described by the coupled 3D FE analysis. The present model based on FSDT which assumes uniform deflection across the thickness cannot predict the nonlinear variation of through the thickness. Moreover, it is noticed that by using more sublayers the deflections tend to be smaller. According to (25), a higher stiffness would be obtained with more constraints on the electric potentials across the thickness. The electric potentials induced by the plate deformation of the bimorph through the direct piezoelectric effects are shown in Figure 3(b). It is observed that the distribution of the electric potential across the bimorph thickness predicted by the present model with more than 2 sublayers is in good agreement with the nonlinear distribution predicted by the coupled 3D FE analysis. Furthermore, it is expected that with more sublayers adopted the nonlinear distribution of the electric potential across the bimorph thickness can be accurately approached without introducing any higher-order electric potential assumptions. However, the conventional linear through-the-thickness electric potential model [28] would be inaccurate to predict the local electric potential response of the piezoelectric bimorph for the case of sensors.

(a)

(b)

Figure 3: Bimorph of under pressure load. (a) Dimensionless deflection. (b) Dimensionless electric potential. 3D FE analysis (full line), present model with (triangles) and present model with (small circles).

Comparisons of the numerical results from the present model with those from the 3D FE analysis are presented in Table 1 for three typical slenderness ratios. The maximum estimating errors for the deflection are 7.29% for , 5.19 for , and 2.61% for . With 10 sublayers adopted, the discrepancy for the maximum values of the electric potential at the plate center is only 0.89% for the thin plate and 4.20% for the thick plate. Consequently, it can be concluded that the present model can provide a good approximation to the global responses such as the deflections. Furthermore, it should be noted that the local responses such as the induced electric potentials of the piezoelectric bimorph sensor can also be well predicted by the present model.

Table 1: Static responses of the piezoelectric bimorph sensor.

Approaches

Error (%)

Error (%)

5
3D FE
−178.93
—
−0.1547
—

Present ()
−169.73
5.14
−0.1710
10.53

Present ()
−165.87
7.29
−0.1612
4.20

10
3D FE
−2487.93
—
−0.6193
—

Present ()
−2374.38
4.56
−0.6839
10.43

Present ()
−2358.72
5.19
−0.6474
4.53

50
3D FE

—
−15.65
—

Present ()

1.85
−16.42
4.92

Present ()

2.61
−15.79
0.89

3.2. Natural Frequencies of the Piezoelectric Bimorph

We then propose the prediction of natural frequencies of the piezoelectric bimorph plate for closed circuit condition on the top and bottom surfaces of the structure for the typical slenderness ratio . The order of Legendre polynomial is chosen as and the mesh in Figure 4(a) is adopted. The natural frequencies of the bimorph plate predicted by the present model are shown in Table 2 in comparison to the results provided by the 3D FE analysis. A rether good agreement between the present model and the 3D FE analysis is observed (maximum error is 5.70%), which indicates that the distribution of the induced electric potential across the thickness does not affect the global natural frequencies much. It can then be safely concluded that the number of sublayers adopted in the present model will not contribute to the accuracy of the global natural frequencies. By using the present SE model with 2 sublayers, the first six modes shapes of the bimorph plate are shown in Figure 5. It is observed that modes 1, 3, and 6 are flexural modes, modes 2 and 5 are torsional modes, and mode 4 is an in-plane shear mode. Therefore, it should be noted that the vibration characteristics of the bimorph plate can be well addressed by the present model. The transversally isotropic piezoelectric material cannot affect the in-plane shear mode, so that the natural frequency of the fourth mode remains the same when using different number of sublayers for the PZT layer, as can be seen in Table 2.

Table 2: Natural frequencies for the bimorph plate in close circuit (Hz).

Modes
1
2
3
4
5
6

3D FE
5989
14938
23276
26091
34035
49112

Present ()
6191
14168
23862
27578
32854
50158

Error (%)
3.37
5.15
2.52
5.70
3.47
2.13

Present ()
6187
14179
23938
27578
32776
50668

Error (%)
3.31
5.08
2.84
5.70
3.70
3.17

(a)

(b)

(c)

Figure 4: The 3 meshes adopted for the bimorph plate.

(a) Mode 1

(b) Mode 2

(c) Mode 1

(d) Mode 2

(e) Mode 1

(f) Mode 2

Figure 5: Mode shapes of the first six modes of the bimorph plate for .

Finally, the convergence of the present SE model is investigated. The bimorph plate is meshed into 2, 3, or 8 quadrilateral elements, as shown in Figures 4(b), 4(c), and 2, respectively. For this case, 2 sublayers are adopted. The first 2 natural frequencies are shown in Table 3 for different order of Legendre polynomials. It can be observed that the monotonic convergence rate of the present method is very fast with respect to the order of the Legendre polynomial and the algorithm is robust with respect to mesh size.

Table 3: Natural frequencies for different order of Legendre polynomials (Hz).

Mesh
2 quadrilaterals
3 quadrilaterals
8 quadrilaterals

Mode
1
2
1
2
1
2

6438
14731
6292
14403
6199
14183

6241
14286
6213
14183
6196
14171

6201
14180
6192
14171
6191
14168

6193
14176
6191
14166
6191
14166

4. Conclusions

An efficient SE model based on the combination of an ESL approach for the mechanical displacement and a sublayer approximation for the electric potential is presented for the static and dynamic analysis of a piezoelectric bimorph. 2D GLL shape functions are used to discretize the displacements. It is observed that the monotonic convergence rate of the developed SE model is very fast with respect to the order of Legendre polynomials. The capability of the present model for prediction of global and local responses such as mechanical displacements, natural frequencies, and the electric potentials across the thickness for static and dynamic processes has been verified by good agreement in numerical solutions with a 3D FE model. The advantage of the present model is that, with appropriate number of sublayers adopted, the nonlinear distribution of the electric potential across the bimorph thickness can be accurately predicted even for rather thick bimorph sensors without introducing any higher-order electric potential assumptions. It is shown that the number of sublayers adopted in the present model will not contribute to the accuracy of the global natural frequencies. Consequently, it can be further concluded that the distribution of the induced electric potential across the thickness of the bimorph sensor does not affect the global natural frequencies much.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The research was supported by National Science Fund for Distinguished Young Scholars (Grant no. 11125209) and Natural Science Foundation of China (Grants nos. 11322215 and 10702039).

References

A. Fernandes and J. Pouget, “Analytical and numerical approaches to piezoelectric bimorph,” International Journal of Solids and Structures, vol. 40, no. 17, pp. 4331–4352, 2003.

S. Y. Wang, “A finite element model for the static and dynamic analysis of a piezoelectric bimorph,” International Journal of Solids and Structures, vol. 41, no. 15, pp. 4075–4096, 2004.

A. Fernandes and J. Pouget, “An accurate modelling of piezoelectric multi-layer plates,” European Journal of Mechanics A, vol. 21, no. 4, pp. 629–651, 2002.

R. P. Khandelwal, A. Chakrabarti, and P. Bhargava, “An efficient hybrid plate model for accurate analysis of smart composite laminates,” Journal of Intelligent Material Systems and Structures, vol. 24, no. 16, pp. 1927–1950, 2013.

A. Fernandes and J. Pouget, “Accurate modelling of piezoelectric plates: single-layered plate,” Archive of Applied Mechanics, vol. 71, no. 8, pp. 509–524, 2001.

J. Fei, Y. Fang, and C. Yan, “The Comparative study of vibration control of flexible structure using smart materials,” Mathematical Problems in Engineering, vol. 2010, Article ID 768256, 13 pages, 2010.

S. Kapuria, “An efficient coupled theory for multilayered beams with embedded piezoelectric sensory and active layers,” International Journal of Solids and Structures, vol. 38, no. 50-51, pp. 9179–9199, 2001.

J. M. Simões Moita, I. F. P. Correia, C. M. Mota Soares, and C. A. Mota Soares, “Active control of adaptive laminated structures with bonded piezoelectric sensors and actuators,” Computers & Structures, vol. 82, no. 17–19, pp. 1349–1358, 2004.

G. R. Liu, X. Q. Peng, K. Y. Lam, and J. Tani, “Vibration control simulation of laminated composite plates with integrated piezoelectrics,” Journal of Sound and Vibration, vol. 220, no. 5, pp. 827–846, 1999.

S. B. Kerur and A. Ghosh, “Active vibration control of composite plate using AFC actuator and PVDF sensor,” International Journal of Structural Stability and Dynamics, vol. 11, no. 2, pp. 237–255, 2011.

J.-P. Jiang and D.-X. Li, “Robust vibration control for smart solar array structure,” Journal of Vibration and Control, vol. 17, no. 4, pp. 505–515, 2011.

X. Q. Peng, K. Y. Lam, and G. R. Liu, “Active vibration control of composite beams with piezoelectrics: a finite element model with third order theory,” Journal of Sound and Vibration, vol. 209, no. 4, pp. 635–649, 1998.

M. Y. Yasin, N. Ahmad, and M. N. Alam, “Finite element analysis of actively controlled smart plate with patched actuators and sensors,” Latin American Journal of Solids and Structures, vol. 7, no. 3, pp. 227–247, 2010.

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.

S. Kapuria and S. D. Kulkarni, “Static electromechanical response of smart composite/sandwich plates using an efficient finite element with physical and electric nodes,” International Journal of Mechanical Sciences, vol. 51, no. 1, pp. 1–20, 2009.

O. Polit and I. Bruant, “Electric potential approximations for an eight node plate finite element,” Computers and Structures, vol. 84, no. 22-23, pp. 1480–1493, 2006.

P. Bisegna and G. Caruso, “Evaluation of higher-order theories of piezoelectric plates in bending and in stretching,” International Journal of Solids and Structures, vol. 38, no. 48-49, pp. 8805–8830, 2001.

D. A. Saravanos, “Mixed laminate theory and finite element for smart piezoelectric composite shell structures,” AIAA Journal, vol. 35, no. 8, pp. 1327–1333, 1997.

S. Kapuria, “A coupled zig-zag third-order theory for piezoelectric hybrid cross-ply plates,” ASME Journal of Applied Mechanics, vol. 71, no. 5, pp. 604–614, 2004.

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.

M. Krommer and H. Irschik, “Reissner-Mindlin-type plate theory including the direct piezoelectric and the pyroelectric effect,” Acta Mechanica, vol. 141, no. 1, pp. 51–69, 2000.

D. A. Saravanos, “Damped vibration of composite plates with passive piezoelectric-resistor elements,” Journal of Sound and Vibration, vol. 221, no. 5, pp. 867–885, 1999.

G. Seriani, “3-D large-scale wave propagation modeling by spectral element method on Cray T3E multiprocessor,” Computer Methods in Applied Mechanics and Engineering, vol. 164, no. 1-2, pp. 235–247, 1998.

D. Komatitsch, C. Barnes, and J. Tromp, “Simulation of anisotropic wave propagation based upon a spectral element method,” Geophysics, vol. 65, no. 4, pp. 1251–1260, 2000.

D. Komatitsch and J. Tromp, “Introduction to the spectral element method for three-dimensional seismic wave propagation,” Geophysical Journal International, vol. 139, no. 3, pp. 806–822, 1999.

P. Kudela, A. Zak, M. Krawczuk, and W. Ostachowicz, “Modelling of wave propagation in composite plates using the time domain spectral element method,” Journal of Sound and Vibration, vol. 302, no. 4-5, pp. 728–745, 2007.

J. N. Reddy, Mechanics of Laminated Composite Plates and Shells: Theory and Analysis, CRC Press, Boca Raton, Fla, USA, 2003.

P. Phung-Van, T. Nguyen-Thoi, T. Le-Dinh, and H. Nguyen-Xuan, “Static and free vibration analyses and dynamic control of composite plates integrated with piezoelectric sensors and actuators by the cell-based smoothed discrete shear gap method (CS-FEM-DSG3),” Smart Materials and Structures, vol. 22, no. 9, pp. 1–17, 2013.