Nonlinear Dynamics of an Unsymmetric Cross-Ply Square Composite Laminated Plate for Vibration Energy Harvesting

: The nonlinear behaviors and energy harvesting of an unsymmetric cross-ply square composite laminated plate with a piezoelectric patch is presented. The unsymmetric cross-ply square composite laminated plate has two stable equilibrium positions by applying thermal stress, thus having snap-through with larger amplitude between the two stable equilibrium positions relative to the general laminated plate. Based on the von-Karman large deformation theory, the nonlinear electromechanical coupling equations of motion of the unsymmetric composite laminated plate with a piezoelectric patch are derived by using Hamilton’s principle. The inﬂuence of the base excitation amplitude on nonlinear behaviors and energy harvesting are investigated. For different base excitation amplitudes, the motions of the system demonstrate periodic motion, quasi-periodic motion, chaotic motion and snap-through, and two single-well chaotic attractors and a two-well chaos attractor coexist. Moreover, the power generation efﬁciency is optimal when the excitation amplitude is in a certain range due to its own unique nonlinear characteristics. The unsymmetric cross-ply square composite laminated plate subjected to thermal stress can actually be called a kind of bistable composite shell structure that has a broad application prospect in combination with morphing aircraft, large deployable antenna and solar panel, which are very likely to have nonlinear vibration.


Introduction
Vibration phenomenon and its energy are ubiquitous in nature and life. There are various forms of vibration in various mechanical equipment, vehicles, ocean waves, human activities and so on. How to effectively use vibration has become the focus of attention. Energy harvesting has been widely concerned because of its advantages and great engineering application value. Energy-harvesting devices can extract energy from mechanical vibration [1,2], human and animal movement [3,4], wind or air flow [5,6], ocean current [7] and road and bridge vibration [8]. Ideally, the vibration energy harvesting can replace chemical battery [9] and can also supply power for wireless sensor network equipment [10]. An electromagnetic energy harvester for freight railcars was presented to enable the smart technologies, of which design, modeling, in-lab and onboard field-tests were implemented [2]. A previously introduced piezoelectric rotational harvester was taken, relying on beam plucking for the energy conversion, to the next step by testing the device during a halfmarathon race [3]. The feasibility of employing these vibrational energy harvesters was considered on flying vertebrates for the purpose of powering a bio-logging device [4]. A dynamic model and behavior of a novel double-beam piezo-magneto-elastic nonlinear wind energy harvester (DBPME-WEH) is investigated, which is a double-beam structure and contains the magnet-induced bistable nonlinearity to enhance the performance of the galloping-based vibration energy harvesting [5]. An experimental study of wind energy harvesting by self-sustained oscillations (flutter) of a flexible piezoelectric membrane fixed in a novel orientation called the "inverted flag" is described [6]. Buoy-type ocean-wave energy converters were designed to exhibit resonant responses when subject to excitation by ocean waves. A novel excitation scheme was proposed which has the potential to improve the energy-harvesting capabilities of these converters [7]. With a nonlinear X-shaped structure connected with piezoelectric harvesters through two types of special mounting configurations (horizontal and vertical cases), novel coupled vibration energy-harvesting systems were purposely constructed and investigated for exploiting nonlinearity and structural coupling effect in vibration energy harvesting [8]. A wide variety of nanomaterials were being explored for technologies that can extract energy from the environment to generate electrical power for sensors and other devices [9]. A ME/EM composite VEH was proposed that can simultaneously obtain a high voltage, large current, high power and wide bandwidth compared to a traditional single ME VEH or single EM VEH [10].
Energy-harvesting devices use a transducer to convert mechanical vibrations to electrical energy. The main methods of transduction include electromagnetic [11], electrostatic [12][13][14], triboelectric [15] and piezoelectric [16,17]. The model of energy harvesting in the environment by electromagnetic conversion has been established maturely, and it has been implemented in many vibration energy harvesters [18]. Amirtharajah and chandrakasan [19] proposed an electromagnetic vibration energy-harvesting system, which is mainly composed of permanent magnet, coil, shell, spring and lumped mass. Mann and Owens [20] proposed an electromagnetic induction energy harvester to transform the vibration energy and use it for electric loads. The energy harvester consists of a movable suspended magnet and several static external magnets. However, the electromagnetic power generation device is vulnerable to the electromagnetic interference of the surrounding environment, and it is not convenient for the miniaturization design. The electrostatic conversion device is composed of electret and variable capacitors, which change the capacity of the variable capacitance diode under vibration excitation so as to realize energy conversion. Its disadvantage is that it needs additional power supply to provide the initial voltage drive [21]. The theoretical support of piezoelectric conversion is the positive piezoelectric effect of piezoelectric materials. Piezoelectric energy harvesting has the advantages of a simple structure, easy system integration, high power density, low cost and miniaturization and has attracted the attention of domestic and foreign scholars. Erturk and Inman [22] gave the analytical expressions of the voltage output and excitation frequency of the piezoelectric cantilever beam energy-harvesting system based on the Euler-Bernoulli beam and verified the reliability of the theoretical analysis through experiments. Anton and Sodano [23] designed energy scavenging devices to capture the ambient energy surrounding the electronics and convert it into usable electrical energy. Ferrari et al. [24] designed a multifrequency vibration energy harvester by using three cantilever beams that have the same size and different mass. The system has three resonance frequency points, but when the external excitation frequency deviates from three resonance frequencies, the output power of the system will also be greatly reduced.
A key design index for energy-harvesting devices is the effective operating frequencies bandwidth at which the energy can be extracted. A linear vibration collector [25,26] is only suitable to work near the natural frequencies; when the vibration energy spectrum of environment is wide, the energy efficiency of linear vibration collector will be very low. In order to enhance the harvesting performance and to broaden the harvesting bandwidth, the research objectives of the vibration energy harvester develop to a wide-frequency band, nonlinear, multi-steady state (two or several steady states), multi-degree of freedom and so on [27,28]. When the vibration frequency of the external environment does not match the natural frequency of the designed energy harvester, the energy conversion rate will be very low. Therefore, a frequency modulation energy-harvesting system is proposed through adjusting the natural frequency of the harvester to make it the same as the external excitation frequency [29,30]. However, in most of the cases, the external vibration frequencies in the actual environment are broadband or change with time. In order to improve the energy-harvesting efficiency, it is found that introducing nonlinearity technology into energy harvesting can effectively improve the frequency bandwidth and response. Daqaq et al. [31] studied the nonlinear monostable energy harvester with soft spring characteristics through an analytical analysis and experimental research, obtained the expressions of the output voltage and power under a monostable state by using the multi-scale method and analyzed the influence of the impedance and coupling coefficient on the energy harvester. Yildirim et al. [32] manufactured and researched a parametrically excited cantilever beam in contact with mechanical stoppers. Yuan et al. [33] investigated the nonlinear behaviors of a vibration-based energy harvester. The harvester consists of a circular composite plate with the clamped boundary, a mass and two steel rings. Jin et al. [34] designed an energyharvesting device to realize the double target of reducing the vibration amplitude and harvesting extremely low ambient vibrations. The device consists of a nonlinear energy sink (NES) and a nonlinear energy harvester with cascaded essential nonlinearities. A more complex nonlinear multi-stable design is an important measure to solve the energy harvesting of low-frequency, ultra-low-frequency and low-intensity vibrations to improve the output power and to expand the harvesting frequency bandwidth. At present, according to the formation mechanisms, the bistable energy-harvesting devices can be divided into three types: a bistable piezoelectric energy-harvesting structure of magnetic repulsion [35,36], a bistable piezoelectric energy-harvesting structure of magnetic suction [37] and a bistable piezoelectric energy-harvesting structure of compression buckling [38]. Since then, the emergence of more complex multi-stable energy capture structures further improves the efficiency of energy harvesting [39,40]. Shukla and Pradyumna [41] explored the use of curvilinear fiber in a bimorph energy harvester and investigated the effect on tailoring the resonance frequency and extracted power. Piezoelectric energy harvesting and piezoelectric driving are similar, except for their ultimate purpose. Gohari et al. [42] utilized inclined piezoelectric actuators to drive the twisting deformation of smart laminated cantilever composite plates/beams; proposed a novel explicit analytical solution for obtaining twisting deformation and optimal shape control and developed a series of simple, accurate and robust finite element (FE) analysis models and realistic electromechanical-coupled FE procedures for verification. Gohari et al. [43] developed an analytical solution for the electromechanical flexural response of smart laminated piezoelectric composite rectangular plates encompassing flexible-spring boundary conditions at two opposite edges that can potentially be applied to real-life structural systems, such as smart floors and bridges, and can be used to analyze the flexural deformation response.
Different from traditional energy-harvesting devices that contain only one stable equilibrium position, this study focuses on nonlinear and bistable characteristics. The novelty of this work lies in that the unsymmetric cross-ply square composite laminated plate has two stable equilibrium positions by applying thermal stress, thus having snapthrough with a larger amplitude between the two stable equilibrium positions relative to the general laminated plate.
In this paper, in order to improve the power generation efficiency, a very efficient nonlinear energy harvesting model is designed. An unsymmetric cross-ply square composite laminated plate and a piezoelectric patch are the main components of the energy-harvesting device. Under the experimental conditions, the center of the plate is connected with the exciter through a supporting bar; therefore, it is assumed that the center of the plate is fixed and four edges are kept free, and a base excitation is given at the center. The residual thermal stress is taken into account to make the unsymmetric composite laminated plate have two stable equilibrium positions. A nonlinear dynamic model of an unsymmetric cross-ply square composite laminated plate with a piezoelectric patch is established. The nonlinear dynamics of the unsymmetric cross-ply square composite laminated plate tend to behave as snap-through and chaotic motions as two single-well chaotic attractors, and as two-well chaos attractors coexist, two-well chaotic motion with large amplitude occurs.
By adjusting the excitation amplitude and frequency, this kind of two-well chaotic motion with a large amplitude will occur in a very wide range of parameters. The nonlinear characteristics of the unsymmetric cross-ply square composite laminated plate can cause snap-through and chaotic motions, which greatly enhance the power generation capacity and greatly improve the power generation efficiency of the energy harvester. When the excitation amplitude is in a specific scope, the motions of the plate around the double-well attractor, the energy harvesting can obtain the highest power generation efficiency and the strongest power generation capacity. Therefore, the research of energy-harvesting devices has great theoretical significance and practical value.
The novelty of this work is that the unsymmetric composite laminated plate has two stable equilibrium positions by applying thermal stress, thus having snap-through with a larger amplitude between the two stable equilibrium positions relative to the general laminated plate.
The unsymmetric cross-ply square composite laminated plate subjected to thermal stress can actually be called a kind of bistable composite shell structure that has a broad application prospect in combination with morphing aircraft, large deployable antenna and a solar panel, which are very likely to have nonlinear vibrations.

Equation of Motion
The unsymmetric cross-ply square composite laminated plate with a piezoelectric patch is depicted in Figure 1. The structure consists of an unsymmetric cross-ply square composite laminated plate and a piezoelectric patch; the piezoelectric patch is pasted on the top surface of the plate, and the ply angle of the unsymmetric cross-ply square composite laminated plate is [0 , as can be seen in Figure 1. A cartesian coordinate oxyz is located in the center of the unsymmetric cross-ply square composite laminated plate. The edge lengths of the unsymmetric cross-ply square composite laminated plate in the x and y directions are 2L x and 2L y , respectively, and the thickness of the unsymmetric cross-ply square composite laminated plate is 2H. The edge lengths of the piezoelectric patch in the x and y directions are L 1 and L 2 , respectively, and the thickness of the piezoelectric patch is h E . Three variables u, v and w are used to denote the displacements of the plate in the x, y and z directions. Under the experimental conditions, the center of the plate is connected with the exciter through a supporting bar. Therefore, it is assumed that the center of the plate is fixed and the four edges are kept free, and a base excitation is given at the center. The unsymmetric cross-ply square composite laminated plate is subjected to a base excitation Y. Figure 2 shows the two stable equilibrium positions of the structure subjected to thermal stress.

Analysis of Unsymmetric Composite Laminated Plate
The displacement field according to the classical deformation plate theory for an arbitrary point is assumed to be where u 0 , v 0 and w 0 denote the longitudinal and transverse displacements of the mid-plane, and z is the transverse coordinate of the plate.

Analysis of Unsymmetric Composite Laminated Plate
The displacement field according to the classical deformation plate theory for an arbitrary point is assumed to be where u0, v0 and w0 denote the longitudinal and transverse displacements of the midplane, and z is the transverse coordinate of the plate.  Based on the von-Karman large deformation theory, the strain-displacement relations can be obtained: Based on the von-Karman large deformation theory, the strain-displacement relations can be obtained: where {ε 0 } represent the membrane strains, and {k} denotes the bending curvature. The stress-strain relationship of the kth layer is given as follows: where α xx , α yy and α xy are the coefficients of the linear expansion. ∆T is the temperature change. Q ij is the material elastic constant, and the expression is shown as follows: in which E 11 and E 22 are the modulus of elasticity, v 12 and v 21 are the Poisson's ratios, G 12 is the shear modulus and all these coefficients belong to the unsymmetric cross-ply square composite laminated plate.
The stress resultants of the unsymmetric cross-ply square composite laminated plate are represented as follows: where A ij are the extensional stiffnesses, D ij are the bending stiffnesses and B ij are the bendingextensional coupling stiffnesses, which are defined in terms of the lamina stiffnesses Q ij as The thermal force and moment resultants are given as

Analysis of Piezoelectric Patch
The stress-strain relationship of the piezoelectric patch is given as follows: where e ij are the piezoelectric stress constants. ∈ ii is the dielectric constant. c E ij are the material elastic constants, and the expression is shown as follows: (15) in which E E 11 and E E 22 are the modulus of elasticity, v E 12 and v E 21 are the Poisson's ratios, G E 12 is the shear modulus and all these coefficients belong to the piezoelectric patch. The stress resultants of the piezoelectric patch are represented as follows: where A E ij are called extensional stiffnesses, D E ij are called the bending stiffnesses and B E ij are called the bending-extensional coupling stiffnesses, which are defined in terms of the lamina stiffnesses c E ij as

Equation of Motion of the Unsymmetric Composite Laminated Plate with the Piezoelectric Patch
The Hamilton's principle will be used in the following analysis to derive the nonlinear governing equations of motions of the unsymmetric cross-ply square composite laminated plate with the piezoelectric patch. The mathematical statement of the Hamilton's principle is given as where T total , U total and δW total are the kinetic energy, the strain energy and the virtual work done of the unsymmetric composite laminated plate with a piezoelectric patch. According to Equations (6)-(10), the potential energy of the unsymmetric cross-ply square composite laminated plate is where S lam and V lam are the area and volume of the unsymmetric cross-ply square composite laminated plate, respectively. n lam is the layer number of the plate.
According to Equations (14)- (16), the strain energy of the piezoelectric patch is where S pie and V pie are the area and volume of the piezoelectric patch, respectively. n pie is the number of piezoelectric patches. The electric potential energy of the piezoelectric patch is The potential energy of the piezoelectric patch is The total potential energy of the unsymmetric cross-ply square composite laminated plate with the piezoelectric patch is The total kinetic energy of the unsymmetric cross-ply square composite laminated plate with the piezoelectric patch is where ρ lam and ρ pie are the density of the unsymmetric cross-ply square composite laminated plate and the piezoelectric patch, respectively. The internal electrical work is given: where E z is the electric field applied in the z direction, and D z is the electric displacement in the z direction. The work of the external force is where c is the damping coefficient. Q(t) is the electric charge quantity and the expression is as follows: The relationship between the current and voltage is as follows: where R is the external electrical resistance. The total potential energy of the unsymmetric cross-ply square composite laminated plate with the piezoelectric patch is According to Reference [44], when the unsymmetric cross-ply square composite laminated plate plus the piezoelectric patch is clamped in the center, the static displacements are described by the assumptions: (1) u s is odd in the x and even in the y directions.
(2) v s is even in the x and odd in the y directions.
(3) w s is even in the x and y directions, and it vanishes at the origin. Therefore, according to Reference [44], the static displacements can be expressed by It is noticed that Equations (31)-(33) only describes the static displacements; in order to describe the dynamic displacements, we introduced the time-related functions u n,m−n (t), v n,m−n (t) and w n,m−n (t) to the replace the functions u n,m−n , v n,m−n and w n,m−n that are independent of the time in Equations (31)- (33). Therefore, the dynamic displacements are considered by the following shape functions.
The dynamic displacements are considered by the following shape functions: It should be pointed out that u n,m−n , v n,m−n and w n,m−n can represent the static curvature of the unsymmetric composite laminated plate, and u n,m−n (t), v n,m−n (t) and w n,m−n (t) can represent the time-varying curvature of the unsymmetric cross-ply square composite laminated plate.
In order to derive the dimensionless equations, the transformations of the variables and parameters are given as follows: Thus, Equations (49)-(51) can be rewritten in dimensionless forms: .. .. .
where the coefficients V 1 , V 2 , V 3 , V 4 , V 5 and V 6 are associated with the output voltage V.

Analysis of Nonlinear Behaviors and Energy Harvesting
In this section, the parameters of the unsymmetric cross-ply square composite laminated plate with the piezoelectric patch are set to be E 11  To present numerical results of the energy harvesting of the unsymmetric cross-ply square composite laminated plate with the piezoelectric patch, the sixth-order Runge-Kutta algorithm is utilized to numerically analyze Equations (53)-(55). For convenience in this study, the overbars in Equations (53)-(55) are dropped in the following analysis. The paper focuses on research of the influence of the foundation excitation amplitude f on energy harvesting. The bifurcation diagrams, waveforms and phase portraits are presented to demonstrate energy harvesting.
The bifurcation diagrams are obtained in Figures 3 and 4. In the figure, f is the external excitation amplitude, and w is the transverse displacement which is the adopted position (0.15, 0.15) of the unsymmetric composite laminated plate, and V is the voltage output. It can be seen from Figures 3 and 4, with the excitation amplitude f increasing, the motions of the unsymmetric composite laminated plate behave as a periodic motion, quasi-periodic motion and chaotic motion. The periodic motion behaves as a vibration, with a small amplitude around one of the two single-well chaotic attractors. The quasiperiodic motion behaves as the oscillation, with a large amplitude around one of the two single-well chaotic attractors. The chaotic motion behaves as an oscillation, with a large amplitude around the two-well chaos attractor. It should be pointed out that the two attractors actually correspond to two equilibrium positions. That is because the unsymmetric composite laminated plate demonstrates two equilibrium positions when the thermal stress is considered, as shown in Equations (12)-(13).
From Figure 3, it also can be found that the unsymmetric cross-ply square composite laminated plate with a piezoelectric patch has two stable positions. When the base excitation amplitude is in a certain range, a chaotic motion with a large amplitude between the two positions will occur. When the base excitation amplitude is out of a certain range, a periodic motion or chaotic motion with a small amplitude around one of the two positions will occur. Compared with Figure 3, the output voltage in Figure 4 also has two stable positions. When the base excitation amplitude is in a certain range, a chaotic motion with a large amplitude occurs, and the corresponding voltage output is relatively large. When the base excitation amplitude is out of a certain range, the voltage output is small. Hence, in order to improve the efficiency of the power generation, appropriate excitation amplitudes should be taken.
Moreover, in order to analyze the detailed motion status and energy harvesting of the unsymmetric cross-ply square composite laminated plate with a piezoelectric patch under different external excitation amplitudes, the dynamic essence and energy harvesting of the above motion region (e.g., Figure 3) is examined via the time domain response, phase portraits, voltage of the harvester and the Poincare map, as shown in Figures 5-13.
By analyzing the phase portrait and time history diagram in Figures 5-13, two singlewell chaotic attractors and a two-well chaos attractor can be observed. In the figures, it is found that the motions behave as a periodic motion around one of the two single-well chaotic attractors, a quasi-periodic motion around one of the two single-well chaotic attractors, a chaotic motion around one of the two single-well chaotic attractors, a chaotic motion around the two-well chaos attractor, a quasi-periodic motion around another one of the two single-well chaotic attractors and a periodic motion around another one of the two single-well chaotic attractors.   From Figures 5-13, the influence of the excitation amplitude on the voltage of the harvester can be clearly identified. In Figure 5, because the base excitation amplitude is too small, the voltage output is weak, which is a small disturbance around the first stable position. Compared with Figure 5, the voltage output is significantly enhanced due to the increasing base excitation amplitude in Figure 6, which is a disturbance with a large amplitude near the first stable position. When the base excitation amplitude reaches a critical value, snap-through with sharp increasing amplitudes occurs in the process of voltage output, as shown in Figure 7. When the base excitation amplitude is in a certain range, constant snap-through with a large amplitude occurs in the process of voltage output, which represents the extremely powerful power generation capacity, as shown in Figures  8-10. From Figures 8-10, it can be found that the output voltage oscillates around two stable positions with a large amplitude, which means a significantly higher generation efficiency. When the base excitation amplitude is beyond the scope, which is shown in Figure 4, the output voltage oscillates near the second stable position with sharp decreas- By analyzing the phase portrait and time history diagram in Figures 5-13, two singlewell chaotic attractors and a two-well chaos attractor can be observed. In the figures, it is found that the motions behave as a periodic motion around one of the two singlewell chaotic attractors, a quasi-periodic motion around one of the two single-well chaotic attractors, a chaotic motion around one of the two single-well chaotic attractors, a chaotic motion around the two-well chaos attractor, a quasi-periodic motion around another one of the two single-well chaotic attractors and a periodic motion around another one of the two single-well chaotic attractors.
pearing accompanying chaotic motion guarantees the highest power generation efficiency and the strongest power generation capacity.
The unsymmetric cross-ply square composite laminated plate subjected to thermal stress can actually be called a kind of bistable composite shell structure. That is to say, this kind of bistable energy-harvesting device has a broad application prospect in combination with morphing aircraft, large deployable antenna and a solar panel, which are very likely to have nonlinear vibrations.   From Figures 5-13, the influence of the excitation amplitude on the voltage of the harvester can be clearly identified. In Figure 5, because the base excitation amplitude is too small, the voltage output is weak, which is a small disturbance around the first stable position. Compared with Figure 5, the voltage output is significantly enhanced due to the increasing base excitation amplitude in Figure 6, which is a disturbance with a large amplitude near the first stable position. When the base excitation amplitude reaches a critical value, snap-through with sharp increasing amplitudes occurs in the process of voltage output, as shown in Figure 7. When the base excitation amplitude is in a certain range, constant snap-through with a large amplitude occurs in the process of voltage output, which represents the extremely powerful power generation capacity, as shown in Figures 8-10. From Figures 8-10, it can be found that the output voltage oscillates around two stable positions with a large amplitude, which means a significantly higher generation efficiency. When the base excitation amplitude is beyond the scope, which is shown in Figure 4, the output voltage oscillates near the second stable position with sharp decreasing amplitudes, which is shown in Figure 11. Figure 11 shows almost the same generation efficiency as Figure 6. Similar to Figure 5, when the base excitation amplitude is too large, the voltage output is again weak, which is a small disturbance around the second stable position in Figures 12 and 13.
By comparing the figures of each group, it can be seen that, when the motions go around single-well attractor, the amplitude of the response is relatively small, while when the motions go around a double-well attractor, the amplitude of the response is quite large. In order to improve the efficiency of power generation, we try to make the motions of the system happen around the double-well attractor.

Conclusions
The nonlinear behaviors and energy harvesting of an unsymmetric cross-ply square composite laminated plate with a piezoelectric patch was studied. The nonlinear electromechanical coupling equation of motion was derived. The influence of the external excitation amplitude on nonlinear behaviors and energy harvesting was analyzed. The main conclusions can be drawn as follows: (1) Compared with the general composite laminated plate, the unsymmetric crossply square composite laminated plates subjected to thermal stress have two equilibrium Based on the above analysis, it can be concluded that, with the change of the base excitation amplitude, the complex nonlinear dynamics of the system behave as a periodic motion, quasi-periodic motion and chaotic motion. When the base excitation amplitude is in a specific scope, snap-through with a large amplitude appears accompanying a chaotic motion. It is based on these complex nonlinear dynamic behaviors that we can determine how to make the energy-harvesting device have the highest power generation efficiency and the strongest power generation capacity. By adjusting the base excitation amplitude, the response of the system is made in a range with large amplitudes to make sure the occurrence of snap-through and chaotic motion. Snap-through with a large amplitude appearing accompanying chaotic motion guarantees the highest power generation efficiency and the strongest power generation capacity.
The unsymmetric cross-ply square composite laminated plate subjected to thermal stress can actually be called a kind of bistable composite shell structure. That is to say, this kind of bistable energy-harvesting device has a broad application prospect in combination with morphing aircraft, large deployable antenna and a solar panel, which are very likely to have nonlinear vibrations.

Conclusions
The nonlinear behaviors and energy harvesting of an unsymmetric cross-ply square composite laminated plate with a piezoelectric patch was studied. The nonlinear electromechanical coupling equation of motion was derived. The influence of the external excitation amplitude on nonlinear behaviors and energy harvesting was analyzed. The main conclusions can be drawn as follows: (1) Compared with the general composite laminated plate, the unsymmetric crossply square composite laminated plates subjected to thermal stress have two equilibrium positions, which can lead to snap-through with a large amplitude, thus greatly improving the efficiency of power generation and greatly enhancing the power generation capacity.
(2) With the change of the base excitation amplitude, the motions of the unsymmetric cross-ply square composite laminated plate demonstrate periodic motion, quasi-periodic motion and chaotic motion.
(3) When the unsymmetric cross-ply square composite laminated plate is subjected to external excitation, two single-well chaotic attractors and a two-well chaos attractor may coexist. When the excitation amplitude is in a specific scope, the motions of the plate around the double-well attractor (meaning the snap-through occurs), at this moment, the unsymmetric cross-ply square composite laminated plate with piezoelectric patch can obtain the highest power generation efficiency and the strongest power generation capacity.
(4) The unique nonlinear characteristics of the designed structure make its own excellent power generation efficiency. Moreover, the power generation efficiency is optimal when the excitation amplitude is in a certain range.
(5) The unsymmetric cross-ply square composite laminated plate subjected to thermal stress can actually be called a kind of bistable composite shell structure that has two stable equilibrium positions. Due to snap-through, bistable composite shell structures are admittedly ideal materials for energy-harvesting devices.
Actually, the conclusions were mainly about the bi-stability of the unsymmetric composite laminated plate subjected to thermal stress, which had two stable equilibrium positions. The reason why this kind of structure is important for energy harvesting devices is actually the snap-through between the two stable equilibrium positions.