Modeling and Chaotic Dynamics of the Laminated Composite Piezoelectric Rectangular Plate

This paper investigates the multipulse heteroclinic bifurcations and chaotic dynamics of a laminated composite piezoelectric rectangular plate by using an extended Melnikov method in the resonant case. According to the von Karman type equations, Reddy’s third-order shear deformation plate theory, andHamilton’s principle, the equations of motion are derived for the laminated composite piezoelectric rectangular plate with combined parametric excitations and transverse excitation. The method of multiple scales and Galerkin’s approach are applied to the partial differential governing equation. Then, the four-dimensional averaged equation is obtained for the case of 1 : 3 internal resonance and primary parametric resonance. The extended Melnikov method is used to study the Shilnikov typemultipulse heteroclinic bifurcations and chaotic dynamics of the laminated composite piezoelectric rectangular plate. The necessary conditions of the existence for the Shilnikov type multipulse chaotic dynamics are analytically obtained. From the investigation, the geometric structure of the multipulse orbits is described in the four-dimensional phase space. Numerical simulations show that the Shilnikov type multipulse chaotic motions can occur. To sum up, both theoretical and numerical studies suggest that chaos for the Smale horseshoe sense in motion exists for the laminated composite piezoelectric rectangular plate.


Introduction
The need for high-speed, light-weight, and energy-saving structures in the aerospace and aviation industry has led to the composite materials instead of traditional materials.Additional requirements for multifunctionality, active vibration, shape control, vibration suppression, and acoustic control have made the development of smart and intelligent structures.A piezoelectric composite laminate is composed of piezoelectric layers which are embedded in laminated composite structures or are boned on the surface of structures.The direct and converse piezoelectric effects are used to suppress the transient vibration and to control the deformation, shape, and buckling of the structures.Such lightweight flexible structures generate large deformations, geometrical nonlinearity, and structural instability when piezoelectric composite laminates are subjected to the coupling between the mechanical and electrical loads.Therefore, it is necessary to study geometrically nonlinear effects on dynamic characteristics of structures in order to accurately design and effectively control vibrations of piezoelectric composite laminate structures.It is very important to investigate the large amplitude nonlinear vibrations of smart structures with piezoelectric materials in order to achieve and predict the desired performance of the systems.
Recently, the studies on dynamics of composite structures with piezoelectric materials have made some progress.Tzou et al. [1] used spatially distributed orthogonal piezoelectric actuators to perform the distributed structural control of elastic shell.They utilized a gain factor and a spatially distributed mode actuator function to describe modal feedback functions.Purekar et al. [2] presented phased array filters with piezoelectric sensors to detect damage in isotropic plates and adopted wave propagation to describe plate dynamics.Ishihara and Noda [3] took into account the effect of transverse shear to analyze the dynamic behavior of the laminate composed of fiber-reinforced laminae and piezoelectric layers constituting a symmetric cross-ply laminate rectangular plate with simply supported edges.Oh [4] considered snap-through thermopiezo-elastic behaviors to 2 Mathematical Problems in Engineering examine the buckling bifurcation and sling-shot buckling of active piezo-laminated plates.Lee et al. [5] employed thirdorder shear deformation theory and nonlinear finite element to canvass deflection suppression characteristics of laminated composite shell structures with smart material laminae.Panda and Ray [6] exploited the first-order shear deformation theory and the three-dimensional finite element method to delve into the open-loop and closed-loop nonlinear dynamics of functionally graded plates with the piezoelectric fiberreinforced composite material under the thermal environment.Dumir et al. [7] used the extended Hamilton's principle to derive the coupled nonlinear equations of motion and the boundary conditions for buckling and vibration of symmetrically laminated hybrid angle-ply piezoelectric panels under in-plane electrothermomechanical loading.Yao and Zhang [8] employed the third-order shear deformation plate theory to explore the bifurcations and chaotic dynamics of the fouredge simply supported laminated composite piezoelectric rectangular plate in the case of the 1 : 2 internal resonances.
The global bifurcations and chaotic dynamics of highdimensional nonlinear systems have been at the forefront of nonlinear dynamics for the past two decades.There are two ways of solutions on Shilnikov type chaotic dynamics of high-dimensional nonlinear systems.One is Shilnikov type single-pulse chaotic dynamics and the other is Shilnikov type multipulse chaotic dynamics.Most researchers focused on Shilnikov type single-pulse chaotic dynamics of highdimensional nonlinear systems.Much research in this field has concentrated on Shilnikov type single-pulse chaotic dynamics of thin plate structures.Feng and Sethna [9] utilized the global perturbation method to study the global bifurcations and chaotic dynamics of the thin plate under parametric excitation and obtained the conditions in which the Shilnikov type homoclinic orbits and chaos can occur.Tien et al. [10] applied the Melnikov method to investigate the global bifurcation and chaos for the Smale horseshoe sense of a two-degree-of-freedom shallow arch subjected to simple harmonic excitation for the case of 1 : 2 internal resonance.Malhotra and Sri Namachchivaya [11] employed the averaging method and Melnikov technique to canvass the local, global bifurcations and chaotic motions of a twodegree-of-freedom shallow arch subjected to simple harmonic excitation for the case of 1 : 1 internal resonance.The global bifurcations and chaotic dynamics were investigated by Zhang [12] for the simply supported rectangular thin plates subjected to the parametrical-external excitation and the parametrical excitation.Yeo and Lee [13] made use of the global perturbation technique to examine the global dynamics of an imperfect circular plate for the case of 1 : 1 internal resonance and obtained the criteria for chaotic motions of homoclinic orbits and heteroclinic orbits.Yu and Chen [14] adopted the global perturbation method to explore the global bifurcations of a simply supported rectangular metallic plate subjected to a transverse harmonic excitation for the case of 1 : 1 internal resonance.
While most of studies are on the Shilnikov type single-pulse global bifurcations and chaotic dynamics of high-dimensional nonlinear systems, there are researchers investigating the Shilnikov type multipulse homoclinic and heteroclinic bifurcations and chaotic dynamics.So far, there are two theories of the Shilnikov type multipulse chaotic dynamics.One is the extended Melnikov method and the other theory is the energy phase method.Much achievement is made in the former theory of high-dimensional nonlinear systems.In 1996, Kovačič and Wettergren [15] used a modified Melnikov method to investigate the existence of the multipulse jumping of homoclinic orbits and chaotic dynamics in resonantly forced coupled pendula.Furthermore, Kaper and Kovačič [16] studied the existence of several classes of the multibump orbits homoclinic to resonance bands for completely integral Hamiltonian systems subjected to small amplitude Hamiltonian and damped perturbations.Camassa et al. [17] presented a new Melnikov method which is called the extended Melnikov method to explore the multipulse jumping of homoclinic and heteroclinic orbits in a class of perturbed Hamiltonian systems.Until recently, Zhang and Yao [18] introduced the extended Melnikov method to the engineering field.They came up with a simplification of the extended Melnikov method in the resonant case and utilized it to analyze the Shilnikov type multipulse homoclinic bifurcations and chaotic dynamics for the nonlinear nonplanar oscillations of the cantilever beam.
The study on the second theory of the Shilnikov type multipulse chaotic dynamics was stated by Haller and Wiggins [19].They presented the energy phase method to investigate the existence of the multipulse jumping homoclinic and heteroclinic orbits in perturbed Hamiltonian systems.Up to now, few researchers have made use of the energy phase method to study the Shilnikov type multipulse homoclinic and heteroclinic bifurcations and chaotic dynamics of highdimensional nonlinear systems in engineering applications.Malhotra et al. [20] used the energy-phase method to investigate multipulse homoclinic orbits and chaotic dynamics for the motion of flexible spinning discs.Yu and Chen [21] made use of the energy-phase method to examine the Shilnikov type multipulse homoclinic orbits of a harmonically excited circular plate.
This paper focuses on the Shilnikov type multipulse orbits and chaotic dynamics for a simply supported laminated composite piezoelectric rectangular plate under combined parametric excitations and transverse load.Based on the von Karman type equations and Reddy's third-order shear deformation plate theory, Hamilton's principle is employed to obtain the governing nonlinear equations of the laminated composite piezoelectric rectangular plate with combined parametric excitation and transverse load.We apply Galerkin's approach and the method of multiple scales to the partial differential governing equations to obtain the fourdimensional averaged equation for the case of 1 : 3 internal resonance and primary parametric resonance.From the averaged equation, the theory of normal form is used to find the explicit formulas of normal form.We study the heteroclinic bifurcations of the unperturbed system and the characteristic of the hyperbolic dynamics of the dissipative system, respectively.Finally, we employ the extended Melnikov method to analyze the Shilnikov type multipulse orbits and chaotic dynamics in the laminated composite piezoelectric plate.In this paper, the extended Melnikov function can be simplified in the resonant case and does not depend on the perturbation parameter.We have used the extended Melnikov method to investigate heteroclinic bifurcations and multipulse chaotic dynamics of the laminated composite piezoelectric plate under the case of 1 : 3 internal resonances.The analysis indicates that there exist the Shilnikov type multipulse jumping orbits in the perturbed phase space for the averaged equations.We present the geometric structure of the multipulse orbits in the four-dimensional phase space.The results from numerical simulation also show that the chaotic motion can occur in the motion of the laminated composite piezoelectric plate, which verifies the analytical prediction.The Shilnikov type multipulse orbits are discovered from the results of numerical simulation.In summary, both theoretical and numerical studies demonstrate that chaos for the Smale horseshoe sense in the motion exists.This paper demonstrates how to employ the extended Melnikov method to analyze the Shilnikov type multipulse heteroclinic bifurcations and chaotic dynamics of high-dimensional nonlinear systems in engineering applications.
The laminated composite piezoelectric rectangular plates are widely applied in space stations, satellite solar panels, sensors, and actuators for the active control of structures and so on.In this paper, we have investigated the multipulse global bifurcations and chaotic dynamics of a laminated composite piezoelectric rectangular plate by using an extended Melnikov method and numerical simulations in detail.We have understood nonlinear vibration characteristics of a laminated composite piezoelectric rectangular plate.Our theoretical results can be used to solve some engineering problems.Since these smart structures are generally light weight and relatively large structural flexibility, laminated composite piezoelectric rectangular plates can induce large vibration deformation during the rapid deployment.In order to eliminate or suppress large vibration and chaotic motion, theoretical results can help optimize the design of the structural parameters of laminated composite piezoelectric rectangular plates.Therefore, the theoretical studies on the multipulse global bifurcations and chaotic dynamics of laminated composite piezoelectric rectangular plates play a very important role in applications in aerospace and mechanical engineering.

Equations of Motion and Perturbation Analysis
We consider a four-edge simply supported laminated composite piezoelectric rectangular plate, where the length, the width, and the thickness are denoted by , , and ℎ, respectively.The laminated composite piezoelectric rectangular plate is subjected to in-plane excitation, transverse excitation, and piezoelectric excitation, as shown in Figure 1.We consider the laminated composite piezoelectric rectangular plate as regular symmetric cross-ply laminates with  layers with respect to principal material coordinates alternatively oriented at 0 ∘ and 90 ∘ to the laminated coordinate axes.Some of layers are made of the piezoelectric materials as actuators, and the other layers are made of fiber-reinforced composite materials.It is assumed that different layers of the symmetric cross-ply composite laminated piezoelectric rectangular plate are perfectly clung to each other, and piezoelectric actuator layers are embedded in the plate.The fiber direction of odd-numbered layers is the -direction of the laminate.The fiber direction of even-numbered layers is the -direction of the laminate.Simply supported plate with immovable edges satisfies the symmetry requirement that eliminates the coupling between bending and extension.However, the displacement of  is free to move at the edge of  = 0, and the displacement of  is free to move at the edge of  = 0. Therefore, the membrane stress is smaller and there exists the coupling between bending and extension.A Cartesian coordinate system  is located in the middle surface of the composite laminated piezoelectric rectangular plate.Assume that (, V, ) and ( 0 , V 0 ,  0 ) describe the displacements of an arbitrary point and a point in the middle surface of the composite laminated piezoelectric rectangular plate in the , , and  directions, respectively.It is also assumed that in-plane excitations of the composite laminated piezoelectric rectangular plate are loaded along the -direction at  = 0 and the -direction at  = 0 with the form of  0 +   cos Ω 1  and  1 +  cos Ω 2 , respectively.Transverse excitation loaded to the composite laminated piezoelectric rectangular plate is expressed as  =  3 cos Ω 3 .The dynamic electrical loading is represented by   =   cos Ω 4 .Considering Reddy's third-order shear deformation description of the displacement field, we have  (, , , ) =  0 (, , ) +   (, , ) V (, , , ) = V 0 (, , ) +   (, , ) (, , , ) =  0 (, , ) , where ( 0 , V 0 ,  0 ) are the deflection of a point on the middle surface, (, V, ) are the displacement components along the (, , ) coordinate directions, and   and   represent the rotation components of normal to the middle surface about the  and  axes, respectively.The nonlinear strain-displacement relations are assumed to have the following form: ( Stress constitutive relations are presented as follows: where   and   denote the mechanical stresses and strains in extended vector notation,    represents the elastic stiffness tensor,   stands for the electric field vector, and   is the piezoelectric tensor. According to Hamilton's principle, the nonlinear governing equations of motion in terms of generalized displacements ( 0 , V 0 ,  0 ,   ,   ) for the composite laminated piezoelectric rectangular plate are given in the previous studies as follows [8]: The simply supported boundary conditions of the composite laminated piezoelectric rectangular plate can be represented as follows [8,22]: The boundary condition (5f) includes the influence of the in-plane load.We consider complicated nonlinear dynamics of the composite laminated piezoelectric rectangular plate in the first two modes of  0 , V 0 ,  0 ,   , and   .It is desirable that we select an appropriate mode function to satisfy the boundary condition.Thus, we can rewrite  0 , V 0 ,  0 ,   , and   in the following forms: By means of the Galerkin method, substituting (6a), (6b), (6c), (6d), (6e) into (4a), (4b), (4c), (4d), (4e), integrating, Mathematical Problems in Engineering and neglecting all inertia terms in (4a), (4b), (4d), and (4e), we obtain the expressions of  01 ,  02 , V 1 , V 2 ,  1 ,  2 ,  3 , and  4 via  1 and  2 as follows: where the coefficients presented in (7a), (7b), (7c), (7d), (7e), (7f) can be found in the previous studies [8].
In order to obtain the dimensionless governing equations of motion, we introduce the transformations of the variables and parameters For simplicity, we drop the overbar in the following analysis.Substituting (5a), (5b), (5c), (5d), (5e), (5f)-( 8) into (4c) and applying the Galerkin procedure, we obtain the governing equations of motion of the composite laminated piezoelectric rectangular plate for the dimensionless as follows: where the coefficients presented in (9a), (9b) are given in the previous studies [8].
The above equations include the cubic terms, in-plane excitation, transverse excitation, and piezoelectric excitation.Equation (9a), (9b) can describe the nonlinear transverse oscillations of the composite laminated piezoelectric rectangular plate.We only study the case of primary parametric resonance and 1 : 3 internal resonances.In this resonant case, there are the following resonant relations: where  1 and  2 are two detuning parameters.The method of multiple scales [23] is employed to (9a), (9b) to find the uniform solutions in the following form: where  0 = ,  1 = .Substituting ( 10) and (11a), (11b) into (9a), (9b) and balancing the coefficients of corresponding powers of  on the left-hand and right-hand sides of equations, the fourdimensional averaged equations in the Cartesian form are obtained as follows: (12d)

Heteroclinic Bifurcations of Unperturbed System
In this section, we focus on studying the nonlinear dynamics characteristic of the unperturbed system.When  = 0, it can be known that system from (21a), (21b), (21c), (21d) is an uncoupled two-degree-of-freedom nonlinear system.The variable  appears in the subspace ( 1 ,  2 ) of (21a), (21b), (21c), (21d) as a parameter since İ = 0. Consider the first two decoupled equations of (21a), (21b), (21c), (21d), Since  1 > 0, (24a), (24b) can exhibit the heteroclinic bifurcations.It is obvious from (24a), (24b) that when  −  6  2 < 0, the only solution to (24a), (24b) is the trivial zero solution, ( 1 ,  2 ) = (0, 0), which is the saddle point.On the curve defined by  =  6  2 , that is, the trivial zero solution bifurcates into three solutions through a pitchfork bifurcation, which are given by  0 = (0, 0) and  ± () = (, 0), respectively, where From the Jacobian matrix evaluated at the nonzero solutions, it can be found that the singular point  0 is the center point and the singular points  ± () are saddle points.It is observed that  and  actually represent the amplitude and phase of vibrations.Therefore, we assume that  ≥ 0 and (25) becomes such that for all  ∈ [ 1 ,  2 ], (24a), (24b) have two hyperbolic saddle points,  ± (), which are connected by a pair of heteroclinic orbits,  ℎ ± ( 1 , ); that is, lim  1 → ±∞  ℎ ± ( 1 , ) =  ± ().Thus, in the full four-dimensional phase space, the set defined by From the results obtained by Feng et al. [9][10][11], it is known that two-dimensional invariant manifold  is normally hyperbolic.The two-dimensional normally hyperbolic invariant manifold  has the three-dimensional stable and unstable manifolds represented as   () and   (), respectively.The existence of the heteroclinic orbit of (24a), (24b) to  ± () = (, 0) indicates that   () and   () intersect nontransversally along a three-dimensional manifold denoted by Γ, which can be written as We analyze the dynamics of the unperturbed system of (21a), (21b), (21c), (21d) restricted to .Considering the unperturbed system of (21a), (21b), (21c), (21d) restricted to  yields where From the results obtained by Feng et al. [9][10][11], it is known that if   ( ± (), ) ̸ = 0,  = constant is called a periodic orbit, and if   ( ± , ) = 0,  = constant is known as a circle of the singular points.Any value of  ∈ [ 1 ,  2 ] at which   ( ± , ) = 0 is a resonant value  and these singular points are resonant singular points.We denote a resonant value by   such that Then, we obtain The geometric structure of the stable and unstable manifolds of  in the full four-dimensional phase space for the unperturbed system of (21a), (21b), (21c), (21d) is given in Figure 2. Since  represents the phase of the oscillations, when  =   , the phase shift Δ of oscillations is defined by The physical interpretation of the phase shift is the phase difference between the two end points of the orbit.In the subspace ( 1 ,  2 ), there exists a pair of heteroclinic orbits connecting to saddle points.Therefore, the homoclinic orbit in the subspace (, ) is, in fact, a heteroclinic connecting in the full four-dimensional space ( 1 ,  2 , , ).The phase shift denotes the difference of the value  when a trajectory leaves and returns to the basin of attraction of .We will use the phase shift in subsequent analysis to obtain the condition for the existence of the Shilnikov type multipulse orbit.The phase shift will be calculated later in the heteroclinic orbit analysis.
We consider the heteroclinic orbits of (24a), (24b).Let  1 =  −  6  2 and  3 =  2 , (24a), (24b) can be rewritten as Set  = 0; (35a), (35b) is a system with the Hamiltonian function When  = 0, there is a heteroclinic loop Γ 0 which consists of the two hyperbolic saddle points  ± () and a pair of heteroclinic orbits  ± ( 1 ).In order to calculate the phase shift and the extended Melnikov function, it is necessary to obtain the equations of a pair of heteroclinic orbits, which are given as follows: We turn our attention to the computation of the phase shift.Substituting the first equation of (37a), (37b) into the fourth equation of the unperturbed system of (21a), (21b), (21c), (21d) and integrating yield where At  =   , there is   ≡ 0. Therefore, the phase shift may be expressed as

Existence of Multipulse Orbits
After obtaining detailed information on the nonlinear dynamic characteristics of the subspace ( 1 ,  2 ) for the unperturbed system from (21a), (21b), (21c), (21d), the next step is to examine the effects of small perturbation terms (0 <  ≪ 1) on the unperturbed system from (21a), (21b), (21c), (21d).The extended Melnikov method developed by Kovačič et al. [15][16][17] is utilized to discover the existence of the multipulse orbits and chaotic dynamics of the nonlinear vibration for the laminated composite piezoelectric rectangular plate.We start by studying the influence of such small perturbations on the manifold .The objective of the research is to identify the parameter regions where the existence of the multipulse orbits is possible in the perturbed phase space.The main aim is to verify whether these parameters satisfy the transversality condition of multipulse chaotic dynamics.It will be shown that these multipulse orbits can occur in the Hamilton system with dissipative perturbations if the parameters meet the transversality condition.The existence of such multipulse orbits provides a robust mechanism for the existence of the complicated dynamics in the perturbed system.In this section, the emphasis is put on the application aspects of the extended Melnikov method to (21a), (21b), (21c), (21d).

Dissipative Perturbations of the Homoclinic Loop.
We analyze dynamics of the perturbed system and the influence of small perturbations on .Based on the analysis by Kovačič et al. [15][16][17], we know that  along with its stable and unstable manifolds is invariant under small, sufficiently differentiable perturbations.It is noticed that  ± () in (24a), (24b) maintains the characteristic of the hyperbolic singular point under small perturbations, in particular,  →   .Therefore, we obtain Considering the last two equations of (21a), (21b), (21c), (21d) yields It is known from the above analysis that the last two equations of (21a), (21b), (21c), (21d) are of a pair of pure imaginary eigenvalues.Therefore, the resonance can occur in (41a), (41b).Also introduce the scale transformations Substituting the above transformations into (41a), (41b) yields where  =  2 6 −  2  1 .When  = 0, (43a), (43b) become The unperturbed system from (44a), (44b) is a Hamilton system with the function The singular points of (44a), (44b) are given as Based on the characteristic equations evaluated at the two singular points  0 and  0 , we can know the stabilities of these singular points.Therefore, it is known that the singular point  0 is a center point.The singular point  0 is a saddle which is connected to itself by a homoclinic orbit.The phase portrait of system for (44a), (44b) is shown in Figure 3(a).
It is found that for the sufficiently small parameter , the singular point  0 remains a hyperbolic singular point   of the saddle stability type.For small perturbations, the singular point  0 becomes a hyperbolic sink   .The phase portrait of the perturbed system from (43a), (43b) is depicted in Figure 3(b).
Using the function (45), at ℎ = 0, and substituting   in (46) into (45), the estimate of the basin of the attractor for  min is obtained as Define an annulus   near  =   as where  is a constant and is sufficiently large so that the unperturbed orbit is enclosed within the annulus.It is noticed that the three-dimensional stable and unstable manifolds of   , denoted as   (  ) and   (  ), are the subsets of the manifolds   (  ) and   (  ), respectively.We will indicate that for the perturbed system, the saddle focus   on   has the multipulse orbits which come out of the annulus   and can return to the annulus in the full four-dimensional space.These orbits, which are asymptotic to some invariant manifolds in the slow manifold   , leave and enter a small neighborhood of   multiple times and finally return and approach an invariant set in   asymptotically, as shown in Figure 4.In Figure 4, this is an example of the threepulse jumping orbit which depicts the formation mechanism of the multipulse orbits.

The 𝑘-Pulse Melnikov Function.
Most researchers focused on Shilnikov type single-pulse chaotic dynamics of the high-dimensional nonlinear systems from the thin plate structures in the past.There exist multipulse chaotic dynamics in the practical engineering systems.The extended Melnikov method is a kind of theory which can be used to investigate the multipulse jumping orbits in the high-dimensional nonlinear systems.Since the theory on multipulse chaotic dynamics is very esoteric and abstract, it is difficult to be extended to solve the engineering problems.Up to now, few researchers have made use of the extended Melnikov method to study the Shilnikov type multipulse homoclinic and heteroclinic bifurcations and chaotic dynamics of highdimensional nonlinear systems in engineering applications.
The extended Melnikov method was first presented by Kovačič et al. [15][16][17], which is an extension of the global perturbation method developed by Feng et al. [9][10][11].Camassa et al. [17] gave the detailed procedure of mathematical proof on the extended Melnikov method, which unifies several disjoint perturbation theoretical methods.This method can be also utilized to detect the Shilnikov type multipulse homoclinic or heteroclinic orbits to the slow manifolds of fourdimensional, near-integrable Hamilton systems or higherdimensional, nonlinear systems.The extended Melnikov function is different from the usual Melnikov function and describes slow dynamics of the multipulse orbits on the hyperbolic manifold.
The key of the extended Melnikov method is how to calculate the extended Melnikov method.The extended Melnikov function is computed by a recursion procedure from the usual 1-pulse Melnikov function and depends on the small perturbation parameter  through a logarithmic function which calculates the asymptotic in the particularly delicate small  limit.In this paper, the extended Melnikov function can be simplified in the resonant case and does not depend on the perturbation parameter.We have used the extended Melnikov method to investigate heteroclinic bifurcations and multipulse chaotic dynamics of the laminated composite piezoelectric rectangular plate.We use the extended Melnikov method described by Kovačič et al. [15][16][17] to find the Shilnikov type multipulse orbits for nonlinear vibration for the laminated composite piezoelectric rectangular plate.We search for the multipulse excursions to find the nondegenerate zeroes of the extended Melnikov function   (, ,  0 ,  2 ) with the certain combination of parameters , ,  0 , and  2 , which we name the -pulse Melnikov function.
We obtain a countable infinity of the singular heteroclinic jumping orbits as follows.Each orbit starts along one branch of the manifold ( 0 ) of the saddle  0 on the annulus .Then, the singular heteroclinic jumping orbit departs from the annulus , goes along one of the singular -pulse orbits Γ  , and lands back at a point on the separatrix that connects the saddle  0 to itself on the annulus .After traveling along the separatrix for a while, the singular heteroclinic jumping orbit takes off again along the singular -pulse orbit Γ  and continues such process.Eventually, the singular heteroclinic jumping orbit lands back on the separatrix.
Therefore, it is concluded that the multipulse orbits of (21a), (21b), (21c), (21d) consist of several portions of the slow time scale on the hyperbolic manifold   and many fast time scale heteroclinic pulses leaving from the manifold   , and these multipulse heteroclinic orbits form a consecutive and recurrence process.

Numerical Results of Chaotic Motions
Based on the above qualitative analysis for the multipulse orbits and chaotic dynamics of the laminated composite piezoelectric rectangular plate, the conditions of the chaotic motion under the sense of the Smale horses are obtained.The heteroclinic bifurcations of (12a), (12b), (12c), (12d) appear when  1 > 0. Therefore, the above theoretical analysis is focused on the situation which there exist heteroclinic bifurcations in (12a), (12b), (12c), (12d).The parameter  1 is the combination of the parameters  6 ,  7 , and  6 , where  1 = (9 6  7 )/2 6 .In this section, we have only performed numerical simulations of the multipulse chaotic motions of the laminated composite piezoelectric rectangular plate under heteroclinic bifurcations in order to further verify the theoretical analysis.Consequently, the parameters  6 ,  7 , and  6 are chosen to satisfy  1 > 0.

Mathematical Problems in Engineering
We choose the averaged equation (12a), (12b), (12c), (12d) to conduct numerical simulations.A numerical approach through the computer software Matlab is utilized to explore the existence of the Shilnikov type multipulse chaotic motions in the laminated composite piezoelectric rectangular plate.Based on the above qualitative analysis, it is found that the damping coefficients  1 ,  2 and transverse excitation  2 play an important role in the multipulse chaotic motions of the laminated composite piezoelectric rectangular plate.In addition, the parameters  2 and  3 are related to the inplane excitation in the -direction and the in-plane excitation in the -direction, respectively.The parameter  4 is the piezoelectric excitation which reflects the characteristics of the piezoelectric material.Hence, the parameters  1 ,  2 ,  4 , and  2 are selected as the controlling parameters to discover the law for complicated nonlinear dynamics of the laminated composite piezoelectric rectangular plate.
We begin to draw bifurcation diagrams of the parameters  2 ,  1 ,  2 , and  4 .Bifurcation diagrams describe the vibration law of the modal displacements  1 and  3 , respectively, when the parameters  2 ,  1 ,  2 , and  4 change in a certain region.We draw bifurcation diagrams according to the rules of the Runge-Kutta algorithm and the Poincaré map theory.For the periodic motions, Poincaré map is of several separate points.For a chaotic motion, the Poincaré map consists of a number of points on the limited Poincaré section.Therefore, it can be observed that chaotic motion and periodic motion of nonlinear system appear from bifurcation diagrams.The chaotic and periodic responses can be identified by several conventional criteria such as phase portraits and Poincaré map.Based on the response law of bifurcation diagrams, phase portraits and Poincaré map are utilized to further verify the existence of the chaotic motions and the periodic motions.In order to compare the influence of these parameters  2 ,  1 ,  2 , and  4 on nonlinear vibration in the laminated composite piezoelectric rectangular plate, we choose the same initial conditions to carry out numerical simulation.7 that the excitation  2 is an important parameter that influences on the nonlinear dynamic responses of the laminated composite piezoelectric rectangular plate.Figure 7 shows that the chaotic motion of the laminated composite piezoelectric rectangular plate appears first, followed by a periodic motion of that.With the increase of excitation  2 , Figure 7 presents the following law: chaotic motion → multi-period motion.
We study the impact of the damping parameter on the nonlinear dynamic responses of the laminated composite piezoelectric rectangular plate.Figure 8 is the bifurcation diagram of the laminated composite piezoelectric rectangular plate with the damping coefficient  1 .The figure demonstrates that system is beginning to enter into the region of the chaotic motion then appears the periodic motion window and finally comes into the region of the chaotic motion again as the damping coefficient  1 varies in the interval  1 = 0.01 ∼ 0.7.Other parameters and initial conditions are the same as those in Figure 7 when excitation is chosen as  2 = 82.7.Figures 8(a  the planes ( 1 ,  1 ) and ( 3 ,  1 ), respectively, as well as the impact of the damping coefficient  1 on the system.Figure 9 portraysthe bifurcation diagram for the laminated composite piezoelectric rectangular plate when the inplane excitation  2 in the -direction varies in the interval  2 = 2∼65.Other parameters and initial conditions remain the same as those in Figure 8 when the damping coefficient  1 is selected as  1 = 0.2.Figures 9(a) and 9(b) display the bifurcation diagram on the plane ( 1 ,  2 ) and ( 3 ,  2 ), respectively.Figure 9 presents that the beginning movement of the system is the periodic motion; then the system appears the chaotic motion.With the increase of the excitation  2 , Figure 9 shows the following evolution law: periodic motion → chaotic motion.Based on the above bifurcation diagram, the excitations  2 ,  2 ,  4 , and the damping coefficient  1 are selected as specific values in order to find the multipulse chaotic motions of the laminated composite piezoelectric rectangular plate.Figure 11 indicates existence of the multipulse chaotic motion of the laminated composite piezoelectric rectangular plate when the excitation  2 is 82.7.In this case, the chosen parameters and initial conditions are the same as those in Figure 7. Figures 11(a) and 11(b) are the three-dimensional phase portrait in the space ( 1 ,  2 ,  3 ) and the Poincaré map on the plane ( 1 ,  2 ), respectively.Figure 11 shows that the excitation  2 has a noticeable effect on the existence of the multipulse chaotic motions on the laminated composite piezoelectric rectangular plate.
Besides the excitations  2 ,  2 ,  4 and the damping coefficient  1 , the multipulse chaotic motions of the laminated composite piezoelectric rectangular plate also depend on other parameters.Figure 12 is obtained when the parameters and initial conditions are chosen as  1 = 14.37,  2 = 11.42, 1 = 0.  and 12, it is found that there are differences in the phase portrait and the Poincaré map, respectively.From the three-dimensional phase portrait in Figure 12, we can see that there exists obvious multipulse jumping phenomenon.The three-dimensional phase portrait is composed of the four regions.The different regions are connected by the multipulse orbits.
In the following numerical simulations, several different sets of parameters and initial conditions are given in order to investigate the different shapes of the multipulse chaotic motion.Figure 13 demonstrates the multipulse chaotic response in the laminated composite piezoelectric rectangular plate for  2 = 92.38.Some parameters and initial conditions are chosen as  1 = 3.61,  2 = 3.13,  5 = −15.01, 8 = 4.09,  10 = −0.01, 20 = −0.09, 30 = −0.05, 40 = −0.05.In this case, other parameters are the same as those in Figure 7. From Figure 13, we can see that there is another shape for the multipulse chaotic motion.It is found that the shapes of these two phenomena depicted in Figures 12 and 13 are completely different.From the threedimensional phase portrait in Figure 13, it is found that multipulse jumping phenomenon is more prominent.

Conclusions
In this paper, the nonlinear vibrations of the laminated composite piezoelectric rectangular plate are studied by applying the theories of the global bifurcations and chaotic dynamics for high-dimensional nonlinear systems.The multipulse heteroclinic orbits and chaotic dynamics are investigated using the extended Melnikov method for the case where the averaged equations have one nonsemisimple double zero and a pair of pure imaginary eigenvalues.The extended Melnikov method can be applied to study the Shilnikov type multipulse heteroclinic bifurcations and chaotic dynamics of highdimensional nonlinear systems in engineering applications.Analysis of the multipulse heteroclinic orbits in the laminated composite piezoelectric rectangular plate demonstrates that such an analysis is a typical singular perturbation problem in which there are two different time scales.Dynamics on the hyperbolic manifold   are of the slow time scale and

Figure 1 :
Figure 1: The model of a laminated composite piezoelectric rectangular plate is given.

Figure 2 :
Figure 2: The geometric structure of manifolds ,   (), and   () is given in the full four-dimensional phase space.

Figure 3 :
Figure 3: Dynamics on the normally hyperbolic manifold is described; (a) the unperturbed case; (b) the perturbed case.

Figure 4 :
Figure 4: The Shilnikov type three-pulse orbits are obtained.

Figure 5 :
Figure 5: The 3-bump orbit with the single-pulse is depicted.

Figure 10 Figure 9 :Figure 10 :
Figure 10 indicates the bifurcation diagram for the laminated composite piezoelectric rectangular plate when the piezoelectric excitation  4 varies from  4 = 2 to  4 = 120.Other parameters and initial conditions remain the same as those in Figure 9 when the in-plane excitation  2 is chosen as  2 = 23.Figures 10(a) and 10(b) demonstrate the bifurcation diagram on the planes ( 1 ,  4 ) and ( 3 ,  4 ), respectively.It is observed from Figure 10 that the piezoelectric excitation  4 has a significant influence on the complicated nonlinear dynamic behaviors of the laminated composite piezoelectric rectangular plate.As the piezoelectric excitation  4 increases, Figure 10 reveals the following law: chaotic motion → multiperiod motion → one-period motion → multiperiod motion → chaotic motion.