Abstract
Taking into account the gravitational potential energy of the piezoelectric energy harvester, the size effect, and the rotary inertia of tip magnet, a more accurate distributed parametric electromechanical coupling equation of tristable cantilever piezoelectric energy harvester is established by using the generalized Hamilton variational principle. The effects of magnet spacing, the mass of tip magnet, the thickness ratio of piezoelectric layer and substrate, and the load resistance and piezoelectric material on the performance of piezoelectric energy capture system are studied by using multiscale method. The results show that the potential well depth can be changed by reasonably adjusting the magnet spacing, so as to improve the energy capture efficiency of the system. Increasing the mass of tip magnet can enhance the output power and frequency bandwidth of the interwell motion. When the thickness of the piezoelectric beam remains unchanged, the optimal load impedance of the system increases along with the increase of thickness ratio of piezoelectric layer and substrate. Compared with the traditional model, which neglects the system gravitational potential energy, the eccentricity, and the rotary inertia of the tip magnet, the calculation results of the frequency bandwidth and the peak power of the modified model have significantly increased.
1. Introduction
The extensive application of wireless sensors in building structural health monitoring, environmental monitoring, and military fields has put forward higher requirements for the cleanliness and sustainability of the power supply. The piezoelectric energy harvester (PEH) can collect vibration energy in the surrounding environment and convert it into sustainable and clean electricity. Due to its simple structure and compact and easy integration, it has great potential to implement wireless sensor nodes and has received extensive attention in recent years [1–4]. The primary structure of the PEH initially is a linear resonator that mainly consists of a piezoelectric cantilever beam with a tip mass. Such linear PEH has a narrow frequency bandwidth that can achieve a peak value only when the external excitation frequency is in agreement with its resonant frequency [5–7]. This makes such linear PEH difficult to meet the requirements of the practical application [8, 9]. To make it more applicable, a variety of complex energy capture systems with active or adaptive techniques have been developed. Among these techniques, nonlinear techniques have been proved to be more suitable for energy harvesting from ambient vibrations in practical applications due to the less sensitivity to variations in the excited frequency of the realistic operational conditions [10,11].
The nonlinearity of the PEHs caused by the magnetic force can be generally classified into three main categories, namely, monostable [12, 13], bistable [14–17], and tristable [18–20] PEHs. Bistable piezoelectric energy harvester (BPEH) has two symmetric potential wells and can oscillate between these two potential wells under low-frequency external excitations, which can greatly enhance the power outputs. Stanton et al. [21] established the distribution parameter model of the cantilever BPEH. Based on the harmonic balance method, the analytical expression of the steady-state response of the energy capture system is obtained and its working characteristics are analyzed. Tang et al. [22] investigated the voltage output of a bistable cantilever piezoelectric energy harvesting system under different excitations and analyzed the influence of magnet spacing on the system response. Kim et al. [23] proposed an electromechanical coupling equation for a hysteresis reversible cantilever-type magnetoelastic BPEH. Wang et al. [24] designed an asymmetric bistable piezoelectric energy harvester composed of a piezoelectric cantilever beam and rotatable magnets and further proposed a method for compensating asymmetric potential wells by appropriately changing the magnet's declination angle. These investigations have shown that the working bandwidth and output power of the BPEH can be greatly improved after oscillating into an interwell motion. However, when the external excitation is low, it is difficult for the BPEH to get rid of the confinement of the potential well and move into the large interwell motion.
In order to overcome this challenge and to produce high energy output at low excitation intensity, tristable piezoelectric energy harvesters (TPEHs) are proposed. Zhu et al. [25] established a TPEH based on the magnetic bistable piezoelectric cantilever beam model, when the angle and distance between two fixed magnets are appropriate, the potential energy function of the energy capture system presents a three-steady state and the potential barrier decreases. Numerical simulation and experimental results show that the TPEH can generate higher output power than the bistable piezoelectric energy harvester in a wider frequency range under lower excitation intensity. Li et al. [26] analyzed a tristable energy harvester which can obtain a high harvesting efficiency at low-frequency base excitation by tristable coherence resonance. Panyam et al. [27] theoretically and experimentally studied the influence of the electromechanical coupling coefficient on the effective bandwidth of the TPEH. Kim et al. [28] proposed potential energy diagrams and found that the tristability of the TPEH is initiated by a new pitchfork bifurcation or a degenerate pitchfork bifurcation that leads to a pair of saddle-node bifurcations. Leng et al. [29] obtained a more precise magnetic force model for the TPEH using equivalent magnetizing current method and the calculation results of magnetic force are in good agreement with experimental data. Yang et al. [18,19] investigated the energy harvesting of the tristable hybrid vibration energy harvester by using geometric nonlinearity to tune the resonant frequency.
In this paper, considering the gravitational potential energy of the system, the eccentricity, and rotary inertia of the tip magnet, a more accurate nonlinear three-steady-state distribution parameter model of the piezoelectric cantilever energy harvester is established based on the generalized Hamilton variational principle, and the analytic solution of the motion equation of the system is obtained by using the multiscale method. In the magnetic dipole model, the influence of eccentricity of the tip magnet is taken into account and the magnetic potential energy expressed by polynomial is obtained. The effects of the distance between the magnets, the mass of the tip magnet, the load resistance, the thickness ratio of the piezoelectric layer to the substrate layer, and the piezoelectric materials on the performance of the energy capture system were studied.
2. Theoretical Model of the TPEH
Figure 1 shows the configuration of the TPEH with magnetic coupling considered in this paper. It mainly consists of a piezoelectric cantilever beam and three magnets (denoted as A, B, and C). The piezoelectric cantilever beam is fixed at the base and is composed of a metal substrate and a pair of piezoelectric layers (PZTs). Two identical PZTs, poled oppositely in the thickness direction, are tightly bonded on the upper and lower surfaces of the substrate. The two PZTs are electrically connected in series with an equivalent load resistance (R) represented as a low power electronic device. Magnet A (named tip magnet) is attached at the tip of the cantilever beam and two external magnets (B and C) are fixed at the frame. The horizontal gap between the tip magnet and magnets B and C is dh. The vertical gap between magnet B and magnet C is 2dv. Here, l and b denote the length and width of the piezoelectric cantilever beam, respectively; hs and tp denote the thickness of the metal substrate and the PZTs, respectively; e is the eccentricity of the tip magnet.

represents the vibration displacement of the base, and s is the coordinate along the neutral axis of the piezoelectric cantilever beam. represents the displacement of the beam at s relative to its fixed end. The modeling of the PEH is based on the linear Euler-Bernoulli beam theory and the linear constitutive equations of the piezoelectric beam are assumed as follows:where Y is Young’s modulus, subscript/superscript p and s represent the piezoelectric layers and substrate layer, S1 and T1 are the strain and the stress of the piezoelectric cantilever beam, respectively. 1 and 3 indicate x and y directions, D3 is the electric displacement, d31 is the piezoelectric constant, and is the dielectric constant. E3 is the electric field.
The strain generated in the piezoelectric cantilever beam can be expressed as
The Lagrange function of the piezoelectric energy harvesting system is as follows:where Tk is the kinetic energy, Ue is the potential energy, We is the electrical energy, Ug is the gravitational potential energy, and Um is the magnetic potential energy. Tk and We are as follows:where , in which and are the density of the piezoelectric layers and substrate layer, respectively. Mt and J are mass and the rotary inertia of the tip mass, respectively. is the permittivity at constant strain.
Ue is expressed aswhere , .
Ug can be expressed as
According to the magnetic dipole model [25] and considering the eccentricity of the tip magnet, the magnetic potential energy can be given bywhere is the magnetic permeability constant. As shown in Figure 1, the horizontal displacement can be evaluated by , where is the rotation angle of the tip magnet A. and are the vector directed from the magnetic moment source of magnet B and C to that of magnet C, respectively.
, and represent the magnetic moment vectors of magnet A, B, and C, respectively, where MA (MB or MC) and VA (VB or VC) are the magnetization intensity and material volume of the magnet A (B or C), respectively.
Using the Galerkin approach, the displacement is assumed as the following form:where and are the R-order mode shape function and the generalized mode coordinates of the cantilever beam, respectively.
The modal shape function satisfies the following orthogonal relations:where is the Kronecker delta, represents the resonance frequency of the r-th mode, and is the eigenvalue. The calculation of the oscillation function and eigenvalues is described in the literature [6].
Substituting equation (9) into (8) and considering only the 1st order mode, Taylor's expansion of Um at can be written aswhere, ,
Expressions for the coefficients and are shown in the appendix.
Considering only the 1st order mode, equation (2) is substituted into the following Lagrangian variational equation using equation (9).where is the generalized dissipative force, denotes the 1st order resonance frequency, is the damping ratio, and represents the generalized output charge. , where R denotes the load impedance. The electromechanical coupling equations of the piezoelectric beam energy harvesting system can be obtained by using equation (13):where , , , , and .
Assume that the excitation acceleration is , where denotes the amplitude of the excitation, denotes the circular frequency, and Cp denotes the capacitance. Introducing the dimensionless parameters shown in equation (17), the electromechanical coupling equations (14) and (16) can be rewritten as equations (18) and (19) in the dimensionless form.where , , , , , , and , .
3. Multiscale Method Analysis
The total potential energy function of the system can be expressed as
Figure 2 shows the potential energy curve of the system at different magnet spacing. For the traditional model (G0 = 0), when , , and , the two symmetric potential wells associated with the equilibria at are separated by the potential well associated with the trivial equilibrium, xs = 0. In this scenario, the potential energy function is symmetrically tristable. With the increase of dv, the middle potential well depth increases while the outer potential well depth decreases when dh remains constant. When the dv is constant, the middle and outer middle potential well depths decrease with the increase of dh. When the modified model (G0 ≠ 0) is used, the potential energy curve will no longer be symmetric, but the variation pattern of potential well depth with magnet spacing is consistent with that of the conventional model.

Introducing a small perturbation parameter , a new independent time variable Tn can be expressed as
The derivative with respect to is
The displacement and output voltage response of the system can be respectively expressed as
The constant parameters such as nonlinear term coefficients, electromechanical coupling coefficients, and the excitation forcing in equation (18) are scaled so that the viscous damping effect appears at the same order of the perturbation problem. That is, we let
Substituting equation (24) into equation (18), we getwhere represents the linearized oscillation frequency within the middle potential well.
To express the nearness of the excitation frequency to the linearized resonance frequency within the middle potential well, we let
Substituting equations (22)–(26) into equation (25) and truncating at order ε and separating the terms of ε and equating them to zero, we obtain
The solutions of equations (27a) and (27b) can be written aswhere A(T1) is a complex-valued function, is the complex conjugate of A, and cc is the complex conjugate of the preceding term. Substituting equations (29) and (30) into equation (28a) and eliminating the secular terms, we get
A can be expressed in the polar form , where a represents the steady-state displacement response amplitude and θ is the argument. Introducing , substituting it into equation (29) and then separating the real and imaginary parts, we obtainwhere , , , and .
For the steady-state response, let the time derivatives in equations (32a) and (32b) be equal to zero, and square and add the resulting equations to obtainwhere the steady-state displacement response amplitude a can be obtained by equation (33). The steady-state solutions for the displacement and the voltage can then be expressed in the following form:where is the steady-state voltage response amplitude , and .
The steady-state solutions of output power can be written as follows:
The stability of the solution can be determined according to the Routh–Hurwitz method.
4. Results and Discussion
In this section, we discuss the magnet spacing, the mass of tip magnet, thickness ratio of piezoelectric layer and substrate, load resistance, and piezoelectric material on the performance of TPEH for two different computational models. The geometric and material properties are as follows: , , , , , , , , . PZT-5 A: , , , . PZT-5H: , , , and . Figures 3 and 4 depict variations of displacement and output power versus excited frequency for different magnet spacing, respectively, in which Mt = 14.9 g, R = 300 kΩ. Figures 3 and 4 show that the displacement and output power amplitude of the interwell motion increase while the bandwidth decreases as dh increases for a fixed dv. The amplitude of the intrawell motion displacement of the conventional model decreases as dh increases, while the amplitude of the intrawell motion displacement of the modified model is not sensitive to changes in dh. With the dv increasing, the displacement and output power amplitude of the interwell motion decreases, while the bandwidth and intrawell motion displacement amplitude increases in the case of the fixed dh. The trend of interwell power varies with the distance between magnets being the same for the conventional model and the modified model, but the peak power and interwell bandwidth of the two models differ significantly. Taking = 7.5 mm, dh = 20 mm and = 7.5 mm, dh = 21 mm as examples, when the modified model is adopted, the corresponding peak power of the two cases is increased by 10.4% and 10.8%, respectively, compared with that of the traditional model, and the corresponding interwell bandwidth is increased by 25.8% and 27.3%, respectively, compared with that of the traditional model.


Figures 5 and 6 give the displacement and output power frequency response curves for different tip magnet masses Mt when = 7.5 mm, dh = 20 mm, and R = 300 kΩ. Figure 5 shows that increasing the tip magnet mass Mt can significantly increase the peak displacement and interwell bandwidth. When the excited frequency is low, the displacement amplitude of the interwell motion increases with tip magnet mass Mt increasing for the conventional model case while there is no significant change in interwell motion displacement amplitude for the modified model. As can be seen from Figure 6, the interwell motion output power amplitude and peak power for both computational model cases increase significantly with Mt increasing. The peak power calculated by the modified model is 0.13 W and 0.23 W, when Mt = 9.9 g and Mt = 14.9 g, respectively, which is 21.2% and 7.4% higher than the peak power of the conventional model, and the interwell motion bandwidth of the modified model is 38.1% and 25.5% higher than that of the conventional model, respectively.


Figures 7 and 8 show the steady-state amplitude response curves of the interwell motion output voltage with the variation of the excited acceleration amplitude for different tip magnet masses Mt when excited frequency ω = 0.5 and ω = 1. Figures 7 and 8 show that there is an excited acceleration threshold, and when the TPEH is excited by an excited amplitude greater than this threshold, it can cross the potential barrier and enter the high energy orbit. The excited acceleration threshold of interwell motion decreases with the increase of Mt when excited frequency ω = 0.5 and ω = 1. Compared to the conventional model, the modified model has a lower excitation acceleration threshold for interwell motion. In addition, the larger the excitation frequency, the larger the excitation acceleration amplitude required for the TPEH to generate interwell motion.


The variation of output power amplitude with load resistance for the two calculation models at low excited frequency is studied when = 7.5 mm, dh = 20 mm, and Mt = 14.9 g. Figure 9 shows that the output power amplitude at each excitation frequency tends to increase at the beginning and decrease afterwards with the increase of load resistance at each excited frequency. Each excitation frequency corresponds to an optimal load resistance to maximize the amplitude of output power of the TPEH, and the optimal load resistance decreases with the increase of excitation frequency. At lower excited frequencies (in the range of 0.6 to 1.1), the optimal load resistance of the modified model is higher than that of the conventional model, but the corresponding peak power is significantly lower than that of the conventional model.

(a)

(b)
The effect of load resistance R on the peak displacement of the cantilever beam for two different piezoelectric materials (PZT-5A, PZT-5H) at = 7.5 mm, dh = 20 mm, and Mt= 14.9 g is shown in Figure 10. The results show that the load resistance R has little effect on the peak displacement of the system under the short circuit or open circuit conditions. The peak displacements of the two piezoelectric materials are very similar to each other. The reason is that they have similar elastic modulus and density, the natural frequencies of short circuit and open circuit are close, and the flexibility under short circuit or open circuit resonance excitation is also very close.

Figure 11 depicts the peak power variation curve with load resistance R for the cantilever beam of two piezoelectric materials (PZT-5A, PZT-5H) when = 7.5 mm, dh = 20 mm, Mt = 14.9 g. The results show that the peak power of the system increases sharply with the increase of the load resistance (the first peak power maximum occurs, and the corresponding local optimal load resistance is Ropt1), then decreases slightly and continues to increase (the second peak power maximum occurs, and the corresponding local optimal load resistance is Ropt2), and then decreases gradually. Both local optimum load resistances of the PZT-5A cantilever beam are larger than those of the PZT-5H cantilever beam, but the maximum values of the system peak power with both piezoelectric materials are very close, due to the fact that the magnitude of the peak power depends on the flexibility of the piezoelectric cantilever beam. The peak power maximum values of the modified model for both piezoelectric materials are larger than those of the conventional model, and the corresponding local optimal load resistance Ropt1 is larger than that of the conventional model, while Ropt2 is smaller than that of the conventional model in both cases. Taking the PZT-5A cantilever beam as an example, the Ropt1 = 59K obtained by using the modified mode is 10.1% higher than that of the traditional model, and the Ropt2 = 281K is 21.6% lower than that of the traditional model. The maximum peak power of the modified model is 7% higher than that of the traditional model.

Figure 12 shows the variation curve of peak power of PZT-5A and PZT-5H cantilever beams with thickness ratio TP/HS of piezoelectric layer and substrate under the two calculation models when the thickness of piezoelectric beam (2tp + hS) is 0.6 mm. It can be seen from Figure 12 that the peak power of PZT-5A and PZT-5H cantilever beams under both calculation models show a trend of increasing sharply and then leveling off with the increase of tp/hs. Under the traditional model, when tp/hs < 0.5, the peak power of PZT-5H cantilever beam is significantly lower than that of PZT-5A cantilever beam; when tp/hs is between 0.5 and 1.1, the peak power of PZT-5H cantilever beam is slightly larger than that of PZT-5A cantilever beam. After that, with the continuous increase of tp/hs, the peak power of PZT-5H cantilever beam is lower than that of PZT-5A cantilever beam again. Under different piezoelectric materials, when tp/hs is small, the peak power obtained by the two models is close to each other, while when tp/hs exceeds a certain value, the peak power obtained by the modified model is obviously higher than that obtained by the traditional model.

The optimal load resistance curves of PZT-5A and PZT-5H cantilever beams with two excitation frequencies (ω = 0.5, ω = 0.5) for dv = 7.5 mm, dh = 20 mm, Mt = 14.9 g, and piezoelectric beam thickness 2tp + hs = 0.6 mm are given in Figure 13. It shows that the optimal load resistance of PZT-5A and PZT-5H cantilever beams under both calculation models decreases with increasing excited frequency and increases with increasing tp/hs, and the growth rate of optimal load resistance slows down with increasing tp/hs. The optimal load resistance of both calculated models is greater for the PZT-5A cantilever beam compared to the PZT-5H cantilever beam for the same excited frequency and tp/hs. Under the same conditions, the optimal load resistance of the modified model is larger than that of the conventional model.

5. Conclusions
In this paper, based on the generalized Hamilton variational principle, considering the gravity potential energy of the energy harvester, the size effect, and the rotary inertia of tip magnet, a more accurate electromechanical coupling equation of the nonlinear tristable piezoelectric cantilever energy harvester is obtained, and the analytical solution of the equation is obtained by using the multiscale method. The effects of magnet spacing, the mass of tip magnet, thickness ratio of piezoelectric layer and substrate, load resistance, and piezoelectric material on the performance of piezoelectric energy capture system under the modified and conventional models were investigated and the following conclusions were obtained:(1)Increasing dh or decreasing dv can improve the interwell displacement and output power amplitude. The interwell bandwidth decreases with the increase of dh and increases with the increase of .(2)Increasing Mt can significantly increase the peak power and bandwidth of the system, but the excited acceleration threshold to generate interwell motion decreases with increasing Mt.(3)The maximum peak power of two piezoelectric cantilever beams (PZT-5A and PZT-5H) is very close. When tp/hs is very small or too large, the peak power of PZT-5A cantilever beam is slightly larger than that of PZT-5H cantilever beam, while when tp/hs is moderate, the peak power of PZT-5H cantilever beam is significantly larger than that of PZT-5A cantilever beam. The optimal load resistance of both piezoelectric material cantilever beams increases with tp/hs increasing. Under the same excited frequency and tp/hs, the optimal load resistance of PZT-5A cantilever beam is larger.(4)The effect on the optimal load resistance of the system is significant after introducing the gravitational potential energy of the TPEH, the eccentricity, and the rotational inertia of the tip magnet into the conventional calculation model. In addition, compared with the traditional model, the peak power and interwell motion bandwidth of the modified model are significantly increased. Therefore, it is of great significance to establish a modified model to improve the accuracy of predicting the energy harvesting efficiency of the TPEH.
Appendix
Data Availability
Requests for the data used to support the findings of this study will be considered by the corresponding author.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was funded by the Doctoral Startup Foundation of Anhui Jianzhu University (2020QDZ07, 2018QD28, and 2017QD05), the Natural Science Research Project of Higher Education Institutions in Anhui Province (KJ2019A0747 and KJ2019A0748), the Anhui Provincial Natural Science Foundation (grant no. 2008085QE245), the Project of Science and Technology Plan of Department of Housing and Urban-Rural Development of Anhui Province (2020-YF15 and 2021-YF22), the ideological and political demonstration course for the Quality Engineering Project of Anhui Province (2020szsfkc0298), the ideological and political demonstration course of Anhui Jianzhu University (2020szkc05), and the research project of Anhui Agricultural University (rc412009).