#### Abstract

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, and Hamilton’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 type multipulse 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.

#### 1. 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 examine the buckling bifurcation and sling-shot buckling of active piezo-laminated plates. Lee et al. [5] employed third-order 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 fiber-reinforced 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 four-edge simply supported laminated composite piezoelectric rectangular plate in the case of the 1 : 2 internal resonances.

The global bifurcations and chaotic dynamics of high-dimensional 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 high-dimensional 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 two-degree-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 high-dimensional 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 four-dimensional 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.

#### 2. 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 , and the displacement of is free to move at the edge of . 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 and 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 and the -direction at with the form of and , respectively. Transverse excitation loaded to the composite laminated piezoelectric rectangular plate is expressed as . The dynamic electrical loading is represented by .

Considering Reddy’s third-order shear deformation description of the displacement field, we havewhere () are the deflection of a point on the middle surface, 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 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 , , , , and . It is desirable that we select an appropriate mode function to satisfy the boundary condition. Thus, we can rewrite , , , , and in the following forms:

By means of the Galerkin method, substituting (6a), (6b), (6c), (6d), (6e) into (4a), (4b), (4c), (4d), (4e), integrating, and neglecting all inertia terms in (4a), (4b), (4d), and (4e), we obtain the expressions of , , , , , , , and via and 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 and 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 .

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 four-dimensional averaged equations in the Cartesian form are obtained as follows:

#### 3. Computation of Normal Form

In order to assist the analysis of the Shilnikov type multipulse orbits and chaotic dynamics of the laminated composite piezoelectric rectangular plate, it is necessary to reduce the averaged equation (12a), (12b), (12c), (12d) to a simpler normal form. It is found that there are and symmetries in the averaged equation (12a), (12b), (12c), (12d) without the parameters. Therefore, these symmetries are also held in normal form.

We take into account the excitation amplitude as a perturbation parameter. Amplitude can be considered as an unfolding parameter when the Shilnikov type multipulse orbits are investigated. Obviously, when we do not consider the perturbation parameter, (12a), (12b), (12c), (12d) becomewhere .

Executing the Maple program given by Zhang et al. [24], the nonlinear transformation used here is given as follows:

Substituting (14a), (14b), (14c), (14d) into (13a), (13b), (13c), (13d) yields a simpler 3rd-order normal form with the parameters for averaged equation (12a), (12b), (12c), (12d) as follows: where the coefficients are , , , and , respectively.

Further, let Substituting (16) into (15a), (15b), (15c), (15d) yields

In order to get the unfolding of (17a), (17b), (17c), (17d) a linear transformation is introduced:

Substituting (18) into (17a), (17b), (17c), (17d) and canceling nonlinear terms including the parameter yield the unfolding as follows: where , , and .

The scale transformations to be introduced into (19a), (19b), (19c), (19d) are

Then, normal form (19a), (19b), (19c), (19d) can be rewritten in the form with the perturbationswhere the Hamiltonian function is of the form and , , , and are the perturbation terms induced by the dissipative effects

#### 4. Heteroclinic Bifurcations of Unperturbed System

In this section, we focus on studying the nonlinear dynamics characteristic of the unperturbed system. When , 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 of (21a), (21b), (21c), (21d) as a parameter since . Consider the first two decoupled equations of (21a), (21b), (21c), (21d),

Since , (24a), (24b) can exhibit the heteroclinic bifurcations. It is obvious from (24a), (24b) that when , the only solution to (24a), (24b) is the trivial zero solution, , which is the saddle point. On the curve defined by , that is, the trivial zero solution bifurcates into three solutions through a pitchfork bifurcation, which are given by and , respectively, where

From the Jacobian matrix evaluated at the nonzero solutions, it can be found that the singular point 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 and (25) becomes such that for all , (24a), (24b) have two hyperbolic saddle points, , which are connected by a pair of heteroclinic orbits, ; that is, . Thus, in the full four-dimensional phase space, the set defined by is a two-dimensional invariant manifold.

From the results obtained by Feng et al. [9–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 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 yieldswhere

From the results obtained by Feng et al. [9–11], it is known that if , = constant is called a periodic orbit, and if , = constant is known as a circle of the singular points. Any value of at which 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

**(a)**

**(b)**

The physical interpretation of the phase shift is the phase difference between the two end points of the orbit. In the subspace , 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 . 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 and , (24a), (24b) can be rewritten as

Set ; (35a), (35b) is a system with the Hamiltonian function

When , there is a heteroclinic loop which consists of the two hyperbolic saddle points and a pair of heteroclinic orbits . 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 . Therefore, the phase shift may be expressed as

#### 5. Existence of Multipulse Orbits

After obtaining detailed information on the nonlinear dynamic characteristics of the subspace for the unperturbed system from (21a), (21b), (21c), (21d), the next step is to examine the effects of small perturbation terms on the unperturbed system from (21a), (21b), (21c), (21d). The extended Melnikov method developed by Kovačič et al. [15–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).

##### 5.1. 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–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) yieldsIt 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) yieldswhere .

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 and , we can know the stabilities of these singular points. Therefore, it is known that the singular point is a center point. The singular point 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).

**(a)**

**(b)**

It is found that for the sufficiently small parameter , the singular point remains a hyperbolic singular point of the saddle stability type. For small perturbations, the singular point 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 , and substituting in (46) into (45), the estimate of the basin of the attractor for 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 three-pulse jumping orbit which depicts the formation mechanism of the multipulse orbits.

##### 5.2. 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 high-dimensional nonlinear systems in engineering applications.

The extended Melnikov method was first presented by Kovačič et al. [15–17], which is an extension of the global perturbation method developed by Feng et al. [9–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 four-dimensional, near-integrable Hamilton systems or higher-dimensional, 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–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 with the certain combination of parameters , and , which we name the -pulse Melnikov function*.*

It is important to obtain the detailed expression of the -pulse Melnikov function. We compute the 1-pulse Melnikov function based on the formula obtained by Kovačič et al. [15–17] at the resonant case . The 1-pulse Melnikov function coincides with the standard Melnikov function . The 1-pulse Melnikov function on both heteroclinic manifolds and is given as follows:

Based on the results given by Kovačič et al. [15–17], it is known that the -pulse Melnikov function is defined as where for and .

It is noticed that the angle function is the complex formula where appears as the argument of a logarithm. When resonance occurs, the periodic orbit corresponding to the value degenerates into a circle of equilibria. Under this case, there exists = 0, . Based on the expression obtained by Kovačič et al. [15–17], the -pulse Melnikov function can be written as follows:

If we set and , (52) can be rewritten as follows:

Based on Proposition 3.1 given by Kovačič et al. [15–17], the nonfolding condition is always satisfied in the resonant case. We obtain the following two conditions:

The main aim of the following analysis focuses on identifying simple zeroes of the -pulse Melnikov function. Define a set that contains all such simple zeroes to be The -pulse Melnikov function has two simple zeroes in the interval

##### 5.3. Geometric Structure of Multipulse Orbits

Based on the aforementioned analysis, we obtain the following conclusions. When the parameters of , , , , , and satisfy condition (54), the -pulse Melnikov function (53) has simple zeroes at and . For or , when the -pulse Melnikov function has no simple zeroes, the stable and unstable manifolds and intersect transversely along a symmetric pair of the two-dimensional, -pulse surfaces . This signifies that the presence of the Shilnikov type -pulse orbits leads to chaotic dynamics in the sense of the Smale horseshoes for the nonlinear motion for the laminated composite piezoelectric rectangular plate. In the phase space of the unperturbed system from (21a), (21b), (21c), (21d), this symmetric pair of the two-dimensional, -pulse surfaces breaks down smoothly onto a pair of limiting -pulse surfaces, , parametrized by (37a), (37b), and (38) with , , and an arbitrary . The sign in (49) is determined by the sign of the corresponding -pulse Melnikov function .

From the discussion given by Kovačič et al. [15–17], it is easily found that for ( or ), the values of the -pulse Melnikov functions are not zero for all , and all have the same sign. It is known that this sign is negative for and positive for . Therefore, the -pulse heteroclinic surfaces and indeed exist, and the limiting -pulse surfaces and also exist when . Since the regions enclosed by the stable and unstable heteroclinic manifolds and are both convex, and the normal vector is known to point out of these manifolds, it demonstrates that the orbits forming each of the surfaces are parametrized by (37a), (37b), and (38) with the alternating signs, and the orbits forming each of the surfaces are parametrized by (37a), (37b), and (38) with the same signs.

For the parameter , there exist orbit segments on the annulus , where the end points of the segments are and , respectively. The trajectories of (44a), (44b) on the segments travel from the end points to in forward time. Therefore, the end points and are, respectively, referred to as the departure and landing points of the heteroclinic jumping . In addition, the line transversely intersects the segments at the end point for , while the line transversely intersects the segments at the end point when . For all , the difference in the coordinates of two end points and is zero, namely, For each , one of the heteroclinic orbits represented by and contained in the limiting surfaces at the value , connects two intersection points and . Therefore, a heteroclinic orbit on the limiting surfaces connects the certain point on the annulus to the end point on the segment . It is also known that a heteroclinic orbit on the limiting surfaces