Abstract
This article provides a method to predict the steadystate response of a piezoelectric cantilever under the excitation of Gaussian colored noise which is a better approximation to the ambient random excitation than the white noise. In this paper, a cantilever energy harvester model, which is fullypaved with piezoelectric ceramic and excited by Gaussian colored noise, is built by the Lagrange’s equations and is nondimensionalized. Then the dynamic equation and piezoelectric equation are transformed into a pair of stochastic differential equations about transient equivalent amplitude and transient phase. After that, the stochastic averaging method is applied to simplify the stochastic differential equations into an Itô type equation about the equivalent amplitude, the drift coefficient and the diffusion coefficient are obtained. This method provides another approach which is different from the stochastic linearization method and moment method. On this basis, some crucial outputs such as the steadystate probability density function (PDF) of the equivalent amplitude, the beam’s displacement, the transient output voltage, the joint probability density function of displacement and velocity as well as the mean square value of output voltage are all obtained. With different noise intensity values and delay factors, the Monte Carlo numerical simulation has verified the correctness of the theoretical prediction.
1. Introduction
Piezoelectric cantilever model which can transfer the ambient energy into electrical energy has attracted many researchers’ attention in recent decades. The efficiency of the transforming is affected by the piezoelectric cantilever’s response characters. At the beginning of this century, Sndano et.al. [1] and Erturk et.al. [2] designed a type of monostablestate piezoelectric cantilever harvester and theoretically predicted the vibrational voltage output. Furthermore, they verified the correctness of this model experimentally. In fact, this type of model only performs large vibration and large voltage output near the fundamental frequency. In other words, when the environmental excitation frequency is far from the fundamental frequency, the voltage output will decrease significantly. After that, scholars extended the piezoelectric harvester model into double well [36], triple well [79], and even tetradic well [10] models with the help of tip magnet and ancillary repulsive magnets. These structures can improve the efficiency of collecting energy. Scholars also found that potential well as shallow as possible and potential barrier as low as possible will widen the working frequency range. Based on these multiwell models, many dynamical characters such as bifurcation, chaos [8], and randomly chaos [11] are studied sufficiently.
Besides the nonlinear forces induced by repulsive magnets, the nonlinearity of this type of model is caused by two main reasons. The first one is the nonlinear strainstress relations of the piezoelectric ceramic; the second one is the geometrical nonlinearity of the cantilever’s deformation. Tim Usher et al. [12], Stanton et.al. [13], Angela Triplett et al. [14] and Guo et al. [15] found out that the nonlinear strainstress relations of piezoelectric ceramic will lead to complex nonlinear dynamical responses even if the cantilever model only has small amplitudes (i.e. the model is under the linear deformation hypothesis). Based on the nonlinear deformation hypothesis, Peyman Firoozy et al. [16] built a new piezoelectric cantilever model, in which he took both the geometrical nonlinearity and the longitudinal inertia nonlinearity into account. Subsequently, the model was quantitatively studied. Some researchers enriched the piezoelectric cantilever study with different academic fields such as aerodynamic and composite material. C. B. Fevzi [17] developed the piezoelectric cantilever with the help of the aerodynamic wing and B. F. C. B. Fevzi [18] and Cakmak et al. [19] deeply investigated the composite beams with different lamination angles.
However, these researches mentioned above are all based on the sinusoidal excitation assumption. Due to the fact that most of the environmental excitations are random, it is very important to extend the research field from deterministic cases into random vibrations. Some researchers used the numerical simulation method to study the monostablestate model [20] and multistablestate [21, 22] piezoelectric cantilever models which were excited by random excitations. However, more researches about the theoretical prediction of the responses of the randomly excited piezoelectric cantilever harvester model are in need. Daqaq et.al applied the moment method to theoretically studied the mean square values of displacement as well as the output voltage of a monostable^{}[23]^{}piezoelectric harvester and a bistable piezoelectric [24] harvester. These harvesters are all excited by the Gaussian white noise as well as Gaussian colored noise. One striking conclusion was obtained that if the model was excited by white noise the nonlinear stiffness term had no influence on the output voltage mean square value. Based on the minimum meansquare error criterion, a kind of equivalent linearization method was proposed by Wenan Jiang et.al. [25] to investigate the output voltage of a piezoelectric harvester excited by white noise. Numerical simulation verified the effectiveness of this method. Compared to the moment method and the equivalent linearization method, the stochastic averaging method can give more detailed predictions on the output values when the nonlinearity of the model cannot be neglected. Wenan Jiang, Liqun Chen et.al. [26] used the standard stochastic averaging method to study the piezoelectric harvester which was externally or parametrically excited by white noise. The influences of different parameters in the model on the output voltage were studied as well. Subsequently, Wenan Jiang, Liqun Chen et.al. [27] developed an improved stochastic averaging method, by which they found that the quadratic nonlinearity and cubic nonlinearity would improve the mean square value of the output voltage and electrical power. These scholars have made great contribution to the model research about the white noise excitation.
However, we must pay attention to the importance of the colored noise excitations. As is well known, Gaussian colored noise, which can be obtained by a zeromean Gaussian white noise passing through a linear filter $\dot{\eta}=\lambda \eta +\lambda \xi \left(t\right)$, is a better approximation of the natural ambient noise than the white noise, and has the correlation function in the form of:
where $\lambda $ denotes the correlation time parameter and $\xi \left(t\right)$ denotes a Gaussian white noise with 0 mean and intensity of 2$D$. The expression of Eq. (1) can be also seen in reference [28]. If $D=$1, $\lambda =\mathrm{}$0.1, 0.2, Fig. 1 shows the diagram of the correlation function $R.$ As is well known, the autocorrelation function of the white noise is $R\left(\tau \right)=\mathrm{E}\left[\xi \right(t\left)\xi \right(t+\tau \left)\right]=D\delta \left(\tau \right)$, the $\delta \left(\tau \right)$ is the Dirac delta function. If the parameter $\lambda $ in Eq. (1) is infinitely large, the colored noise will decrease to be a Gaussian with noise. G. Falsone and I. Elishakoff [29] developed a stochastic linearization technique for colored noise excited Duffing oscillator. After that, Bo Li et.al. [30] and Gen Ge. et al. [31] extended the stochastic averaging method into wider fields such as nonlinear oscillator with stiffness nonlinearity and inertia nonlinearity.
The detailed response characters of colored noise excited piezoelectric cantilever model, such as stationary probability density function (PDF) of displacement, PDF of voltage, joint PDF of displacement and velocity and mean square values of voltage, haven’t been studied yet. Thus, a method on stationary responses of piezoelectric cantilever with colored noise excitation is proposed in this paper.
The research subject is a uniform cantilever with two layers of piezoelectric ceramic on the top as well as the bottom. The linear strainstress relation as well as the nonlinear deformation is taken into account. Firstly, the dynamic model is built by the Lagrange equation and is transformed into dimensionless form. Secondly, a pair of differential equations of amplitude and phase is obtained by coordinate transformation. After that, the pair of equations is simplified to be an Itô type stochastic differential equation. Thirdly, the stationary probability density functions about the amplitude, the displacement and the output voltage, the joint probability density function of the displacement and velocity, as well as the mean square value of voltage, are all obtained by solving the averaged FokkerPlankKolmogorov (FPK) equation. At last, the influence of the noise density and the noise delay time on these PDFs are also studied. The Monte Carlo simulations verifies this procedure.
Fig. 1The correlation function of the colored noise
2. The modeling
As is shown in Fig. 2, the piezoelectric cantilever is characterized as: the thickness of the piezoelectric ceramic ${t}_{p}$, the length of the base beam ${L}_{b}$, the width of the cantilever. The vector 1 denotes the longitudinal direction, while vector 3 denotes the transversal direction. $z\left(t\right)$ is the pedestal excitation in the form of colored noise.
Fig. 2Cantilever beam with piezoelectric sheets fully covered under excitation of base noise
The strainstress relation of the piezoelectric ceramic is expressed as [1]:
where, ${T}_{1}^{p}$ denotes the 1direction stress of ceramic, ${S}_{1}^{p}$^{}denotes the 1direction strain of ceramic, ${C}_{11}^{E}$^{}denotes Young’s modulus of elasticity of piezoelectric ceramic, ${E}_{3}$ denotes 3direction electric field intensity, ${D}_{3}$ denotes 3direction potential shift.
Although every mode can be excited by the colored noise, the first mode’s amplitude predominates those of the lower modes. Thus, only the first mode is studied. It is assumed that the point away from the base has the transverse displacement $w\left(s,t\right)=\phi \left(s\right)q\left(t\right)$, where $\phi \left(s\right)$ is the first order mode function, and $q\left(t\right)$ denotes the modal coordinate. The piezoelectric cantilever stores up the kinetic energy, the elastic potential energy and the electric energy. The kinetic energy of the base layer and that of one piece of piezoelectric ceramic are as follows:
In Eq. (4) and Eq. (5), ${\rho}_{b}$ denotes the density of the base layer, ${\rho}_{p}$ denotes the density of the ceramic, $\dot{w}+\dot{z}$ denotes the transverse velocity.
The strain of the point away from the neutral layer with distance y is expressed as: ${S}_{1}^{p}=yw\text{'}\text{'}$.
The total elastic energy is:
where ${v}_{p}$ denotes the volume of the piezoelectric ceramic, ${I}_{p}=b{\int}_{\frac{{t}_{b}}{2}}^{\frac{{t}_{b}}{2}+{t}_{p}}{y}^{2}dy$ is the inertia moment to the neutral layer, ${{I}_{p}}^{\text{'}}=b{\int}_{\frac{{t}_{b}}{2}}^{\frac{{t}_{b}}{2}+{t}_{p}}{y}^{4}dy$ is the fourthorder moment to the neutral layer, $d=\frac{{t}_{b}+{t}_{p}}{2}$ denotes the distance from the center of the ceramic layer to the neutral layer. Thus, one has:
$+{{C}_{111}}^{E}{{I}_{p}}^{\text{'}}{\int}_{0}^{L}{{w}^{\text{'}\text{'}}}^{4}ds\frac{u\left(t\right)}{2}{e}_{31}bd{\int}_{0}^{L}{{w}^{\text{'}\text{'}}}^{2}ds.$
The work done by the electric force is:
The Lagrange function is defined as $L=TU+{W}_{e}$. Considering the damping dissipation energy $\psi ={\int}_{0}^{L}\frac{1}{2}C{\dot{w}}^{2}ds$, the Lagrange equation is written as follows:
Integrating ${D}_{3}$ on the ceramic surface aero, the electric quantity $Q\left(t\right)$ passing through ${R}_{L}$ is:
Taking derivative to $Q\left(t\right)$ with respect to time $t$, and considering the Ohm law, the electric current $I\left(t\right)$ satisfies:
where, the electric capacity is ${C}_{p}=\frac{{{\epsilon}_{33}}^{S}b{L}_{p}}{2{t}_{p}}$.
The piezoelectric function is:
Introducing the character length and the character time scale, the dimensionless transformation is as follows:
The dynamic equation and the piezoelectric equation are obtained as:
where parameters are:
${h}_{4}=\frac{bd{e}_{31}}{L}{\int}_{0}^{1}\frac{{\partial}^{2}\phi}{{L}^{2}\partial {\zeta}^{2}}d\zeta ,{h}_{5}=\frac{\rho A{L}^{2}}{{T}^{2}}{\int}_{0}^{1}\phi d\zeta ,{c}^{\mathrm{*}}=\frac{C{L}^{2}}{T}{\int}_{0}^{1}{\phi}^{2}d\zeta ,$
$\beta =\frac{bd{e}_{31}}{{C}_{p}}{\int}_{0}^{1}\frac{{\partial}^{2}\phi}{\partial {\zeta}^{2}}d\zeta ,\gamma =\frac{T}{{C}_{p}{R}_{L}}.$
By defining $\mu =\frac{{c}^{\mathrm{*}}}{{h}_{1}}$, ${{\omega}_{0}}^{2}=\frac{{h}_{2}}{{h}_{1}}$, $\alpha =\frac{{h}_{3}}{{h}_{1}}$, $\vartheta =\frac{{h}_{4}}{{h}_{1}}$, $\eta =\frac{{h}_{5}}{{h}_{1}}\ddot{f}$, Eq. (13) and Eq. (14) can be transformed as:
Parameters in Eq. (15) and Eq. (16) are as follows: ${\omega}_{0}$ denotes the fundamental frequency, $\alpha $ is the cubic nonlinearity coefficient, $\mu $ is the damping ratio, $\beta $ and $\vartheta $ are electromechanical coupling coefficient in piezoelectric equation and that in dynamic equation, $\eta \left(t\right)$ denotes the colored noise.
3. Stochastic averaging
Standard stochastic averaging method is an effective tool in dealing with quasilinear oscillators subject to wide band noise. When the nonlinear stiffness term is weak enough, it has weak influence on the responsive frequency. The model we use in this paper has weak nonlinearity when the deformation is small.
For convenience, the dimensionless time $\tau $ is still expressed as $t$:
where $\theta ={\omega}_{0}t+\phi $.
Substituting Eq. (17) into Eq. (16), one has:
Then the voltage can be expressed as:
where $B=\sqrt{{{B}_{1}}^{2}+{{B}_{2}}^{2}}=\frac{A\beta {\omega}_{0}}{\sqrt{{\gamma}^{2}+{{\omega}_{0}}^{2}}}$, ${\phi}_{u}=\phi +\mathrm{\Delta}\phi $, $\mathrm{t}\mathrm{a}\mathrm{n}\phi =\frac{{B}_{1}}{{B}_{2}}=\frac{\gamma}{{\omega}_{0}}$.
The presumption of the transformation is that the nonlinear term has weak influence on the frequency. Subsequently the responsive frequency is still assumed to be ${\omega}_{0}$. The correctness of this assumption will be verified in the numerical simulation section.
Taking derivative to $x=A\mathrm{c}\mathrm{o}\mathrm{s}\theta $, $y=A{\omega}_{0}\mathrm{s}\mathrm{i}\mathrm{n}\theta $ with respect to time $t$ yields:
Solving the two equations, one can obtain:
By now, a pair of stochastic differential equation has been solved out which can be expressed in the standard form:
where:
An Itô type stochastic differential equation of the limiting diffusion process can be given as:
where $m\left(A\right)$ denotes the drift coefficient, ${\sigma}^{2}\left(A\right)$ denotes the diffusion coefficient, $W\left(t\right)$ is standard unit Brownian motion (see reference [29, 30]). The operator $\u27e8\u2022\u27e9=\frac{1}{2\pi}{\int}_{0}^{2\pi}\u2022d\theta $ does averaging with respect to time, $R\left(\tau \right)$ denotes the noise correlation function, 2D is the noise strength, $\lambda $ is the time delay coefficient:
where $E\u27e8\u2022\u27e9$ denotes the mathematical expectation, $\delta \left(\tau \right)$ is the Diracdelta function.
For application, people pay more attention to the amplitude than to the phase. Using the Eq. (26), the drift coefficient $m\left(A\right)$, and the diffusion coefficient ${\sigma}^{2}\left(A\right)$ are given as:
Thus the Eq. (22) can be averaged to be Eq. (30):
4. Stationary responses
The averaged FokkerPlankKolmogorov (FPK) equation associated with Itô equation (30) is of the form:
where $p(A,t)$ is the transition probability density of displacement amplitude $A$. Under assumption that zero probability flows at the two boundaries $A=0$ and $A=\mathrm{\infty}$, the stationary solution of FPK equation $p\left(A\right)$ for Eq. (31) is of the form:
where constant $N$ is the normalization constant.
Substituting Eq. (28), Eq. (29) into Eq. (32) yields the stationary PDF of the equivalent amplitude:
Considering the relation $B=\frac{A\beta {\omega}_{0}}{\sqrt{{\gamma}^{2}+{{\omega}_{0}}^{2}}}$ one can deduce the PDF of the output voltage as:
The joint stationary probability density of displacement$x$and velocity $y=\dot{x}$ can be further obtained from $p\left(A\right)$ as follows:
Subsequently, the PDF of displacement $p\left(x\right)$ can be gained by integrating $y$ form –∞ to ∞ in Eq. (35):
On the other hand, by integrating $x$ form –∞ to ∞ in Eq. (35), the PDF of velocity $p\left(y\right)$ can also be obtained.
Considering the relation $p\left(u\right)={\left.\frac{1}{2\pi \omega B}p\left(B\right)\right}_{B\to \sqrt{{u}^{2}}}$, PDF of the output voltage $p\left(u\right)$ is:
By Eq. (18), one knows $u=\frac{\beta {{\omega}_{0}}^{2}}{{\gamma}^{2}+{{\omega}_{0}}^{2}}x\frac{\beta \gamma}{{\gamma}^{2}+{{\omega}_{0}}^{2}}y$, thus the mean square value of voltage is:
From Eq. (33) to Eq. (38), one has obtained all theoretical predictions in need to investigate the efficiency of collecting the electric energy. In the nest section, necessary numerical simulations will be carried out to verify these equations.
5. Numerical simulation
The Monte Carlo method is applied to verify the effectiveness of the above Eq. (33), Eq. (35), Eq. (35), Eq. (37) and Eq. (38). The parameters are set as: ${\omega}_{0}=$ 1, $\mu =$0.05, $\alpha =$ 1, $\beta =$ 0.01, $\gamma =$0.1, $\vartheta =$ 0.05. The parameter D is chosen as $D=$0.001 and $D=$0.002, separately. As is well known, the larger the correlation time parameter is, the narrower the noise band becomes. The time delay coefficient is chosen as $\lambda =$0.5 and $\lambda =$5 for comparison. Totally 3,000 sets of colored noise are imposed to the Eq. (15) and Eq. (16). Each set of noise includes 30,000 numbers. The time step is set as $\mathrm{\Delta}t=\text{0.005}$. The numerical solutions could be obtained by an order2 stochastic RungeKutta algorithm [31]. The last 10,000 steps are kept as the stationary responses. For each step the transient amplitude $A={x}^{2}+\frac{{y}^{2}}{{{\omega}_{0}}^{2}}$ is calculated out. Finally, do statistic on total stationary $n=$3,000×10,000 values of $A$, $x$, and $u$, and the stationary probability density functions of the amplitude $p\left(A\right)$, of the displacement $p\left(x\right)$, of the voltage $p\left(u\right)$ are shown in Fig. 3Fig. 5.
Fig. 3The analytical solution and numerical solution of stationary PDF of amplitude (solid lines: Eq. (33) dots: Monte Carlo simulation on Eq. (15) and Eq. (16))
Fig. 4The analytical solution and numerical solution of stationary PDF of displacement (solid lines: Eq. (36); dots: Monte Carlo simulation on Eq. (15) and Eq. (16))
From Fig. 3Fig. 5, one can easily find that larger noise density will induce wider spread of the amplitude, displacement and voltage. And, with the same noise density, smaller time delay coefficient $\lambda $ i.e. wider noise band, will clearly lead to wider spread of these three variables. Also, it is obvious that the numerical simulation coincides with the theoretical prediction.
Subsequently, the joint PDF of displacement $x$ and velocity $y=\dot{x}$ are numerically simulated. We take 40×40 grids in a range $\left\{0.5\le x\le 0.5,0.5\le y\le 0.5\right\}$ with a space gap: $\mathrm{\Delta}x=\mathrm{\Delta}y=\text{0.025}$. Then we count the numbers ${n}_{ij}$ ($i$, $j=$140) of the $n=$4,000×10,000 stationary $\left(x,y\right)$ dots which fall within each grid. Finally, the joint PDF $p\left(x,y\right)$ value of each grid can be calculated out by $p\left({x}_{i},{y}_{i}\right)=\frac{{n}_{ij}}{n\times \mathrm{\Delta}x\times \mathrm{\Delta}y}$. The joint PDF are illustrated in Fig. 6.
Fig. 5The analytical solution and numerical solution of stationary PDF of voltage (solid lines: Eq. (37); dots: Monte Carlo simulation on Eq. (15) and Eq. (16))
Fig. 6The joint stationary PDF of displacement x and velocity y: a), b) for D= 0.001, λ = 0.5, c), d) for D= 0.001, λ = 5; e, f) for D= 0.002, λ = 0.5; a), c), e) are given by Eq. (35); b), d), f) are obtained by Monte Carlo simulation)
a)
b)
c)
d)
e)
f)
Obviously, Fig. 6 gives a full view of the probability density function of the displacement and velocity. If one do integral to PDF $p\left(x,y\right)$ with respect to variable y from $\mathrm{\infty}$ to $+\mathrm{\infty}$, the joint PDF will degrade to the PDF of displacement $p\left(x\right)$ as shown in Fig. 4. Similarly, the larger the noise density $D$ is and the smaller the time delay coefficient is, the lower peak value of the Fig. 6 becomes.
Fig. 7Mean square values of voltage with different coefficients: a) for D= 0.002, μ= 0.05; b) for D= 0.002, λ= 0.5; The solid line is the theoretical prediction value obtained by Eq. (38), and the hollow circle is the numerical result obtained by Monte Carlo simulation)
a)
b)
Mean square value of the output voltage $E\left({u}^{2}\right)$ is an important index in estimating the piezoelectric cantilever’s behavior. One can see the output $E\left({u}^{2}\right)$ is basically proportional to the time delay factor $\lambda $ from 0.3 to 0.9. As is well known, the larger the time delay factor $\lambda $ is, the wider the noise band becomes. There comes the conclusion that the wider band noise will lead to larger voltage output even if the noise has the same noise density. Also, it is not hard to understand that when the damping increases the output will decrease in Fig. 7(b). This finding is consistent with the conclusion given by Daqaq [23].
6. Conclusions
A piezoelectric cantilever model subject to colored noise is studied and the stationary responses are theoretically predicted by applying the stochastic averaging method. Generally, this is an effective method which enriches the methods in random vibrational field of the piezoelectric cantilever. Some conclusions are summarized as follows:
1) The Cubic nonlinear coefficient $\alpha $ has small influence on the responses. In the theoretical prediction equations of amplitude, displacement, voltage, (see Eq. (3439)), the coefficient $\alpha $ was eliminated during the stochastic averaging procedure. But the good fitness between the equations and the numerical simulation shows this coefficient $\alpha $ has few impacts on these responses. This phenomenon coincides with the reference [24].
2) Stronger noise intensity will not only lead to larger vibration amplitude and displacement but also lead to lager output mean square voltage. Time delay factor also has a significant influence on the output mean square voltage. With the same noise density level, the larger the time delay factor is (i.e. the wider the noise band is), the larger the voltage output becomes.
References

H. A. Sodano, G. Park, and D. J. Inman, “Estimation of electric charge output for piezoelectric energy harvesting,” Strain, Vol. 40, No. 2, pp. 49–58, May 2004, https://doi.org/10.1111/j.14751305.2004.00120.x

A. Erturk and D. J. Inman, “An experimentally validated bimorph cantilever model for piezoelectric energy harvesting from cantilevered beams,” Smart Materials and Structures, Vol. 18, No. 2, pp. 025009–18, Feb. 2009, https://doi.org/10.1088/09641726/18/2/025009

S. Sun and S. Q. Cao, “Dynamic modeling and analysis of bistable piezoelectric cantilever beam power generation system,” (in Chinese), Chinese Journal of Physics, Vol. 61, No. 21, 2012.

W. Yang and S. Towfighian, “A hybrid nonlinear vibration energy harvester,” Mechanical Systems and Signal Processing, Vol. 90, pp. 317–333, Jun. 2017, https://doi.org/10.1016/j.ymssp.2016.12.032

K. Guo, S. Cao, and S. Wang, “Numerical and experimental studies on nonlinear dynamics and performance of a bistable piezoelectric cantilever generator,” Shock and Vibration, Vol. 2015, pp. 1–14, 2015, https://doi.org/10.1155/2015/692731

M. H. Yao, Y. B. Li, and W. Zhang, “Nonlinear dynamics of longitudinal auxiliary magnetic bistable piezoelectric cantilever beam,” (in Chinese), Journal of Beijing University of Technology, Vol. 41, No. 11, pp. 1756–1760, 2015.

J. Cao, S. Zhou, W. Wang, and J. Lin, “Influence of potential well depth on nonlinear tristable energy harvesting,” Applied Physics Letters, Vol. 106, No. 17, p. 173903, Apr. 2015, https://doi.org/10.1063/1.4919532

G. T. Oumbé Tékam, C. A. Kitio Kwuimy, and P. Woafo, “Analysis of tristable energy harvesting system having fractional order viscoelastic material,” Chaos: An Interdisciplinary Journal of Nonlinear Science, Vol. 25, No. 1, p. 013112, Jan. 2015, https://doi.org/10.1063/1.4905276

H.X. Zou et al., “A broadband compressivemode vibration energy harvester enhanced by magnetic force intervention approach,” Applied Physics Letters, Vol. 110, No. 16, p. 163904, Apr. 2017, https://doi.org/10.1063/1.4981256

C. Wang, Q. Zhang, and W. Wang, “Wideband quinstable energy harvesting via combined nonlinearity,” AIP Advances, Vol. 7, No. 4, p. 045314, Apr. 2017, https://doi.org/10.1063/1.4982730

S. Sun and S. Q. Cao, “Response analysis of bistable piezoelectric power generation system excited by white noise,” (in Chinese), Piezoelectric and AcoustoOptic, Vol. 37, No. 6, pp. 969–972, 2015.

T. Usher and A. Sim, “Nonlinear dynamics of piezoelectric high displacement actuators in cantilever mode,” Journal of Applied Physics, Vol. 98, No. 6, p. 064102, Sep. 2005, https://doi.org/10.1063/1.2041844

S. C. Stanton, A. Erturk, B. P. Mann, and D. J. Inman, “Nonlinear piezoelectricity in electroelastic energy harvesters: Modeling and experimental identification,” Journal of Applied Physics, Vol. 108, No. 7, p. 074903, Oct. 2010, https://doi.org/10.1063/1.3486519

A. Triplett and D. D. Quinn, “The effect of nonlinear piezoelectric coupling on vibrationbased energy harvesting,” Journal of Intelligent Material Systems and Structures, Vol. 20, No. 16, pp. 1959–1967, Nov. 2009, https://doi.org/10.1177/1045389x09343218

K. K. Guo and S. Q. Cao, “Analysis of the main resonance response of a piezoelectric cantilever beam considering material nonlinearity,” (in Chinese), Vibration and Shock, Vol. 33, No. 19, pp. 8–16, 2014.

P. Firoozy, S. E. Khadem, and S. M. Pourkiaee, “Broadband energy harvesting using nonlinear vibrations of a magnetopiezoelastic cantilever beam,” International Journal of Engineering Science, Vol. 111, pp. 113–133, Feb. 2017, https://doi.org/10.1016/j.ijengsci.2016.11.006

L. Tang, Y. Yang, and C.K. Soh, “Improving functionality of vibration energy harvesters using magnets,” Journal of Intelligent Material Systems and Structures, Vol. 23, No. 13, pp. 1433–1449, Sep. 2012, https://doi.org/10.1177/1045389x12443016

F. C. Bolat, “An experimental analysis and parametric simulation of vibrationbased piezoaeroelastic energy harvesting using an aerodynamic wing profile,” Arabian Journal for Science and Engineering, Vol. 45, No. 7, pp. 5207–5214, Jul. 2020, https://doi.org/10.1007/s13369020043731

F. Cakmak Bolat, S. Basaran, and M. Kara, “Investigation of energy harvesting in composite beams with different lamination angles under dynamic effects,” Composite Structures, Vol. 270, p. 114056, Aug. 2021, https://doi.org/10.1016/j.compstruct.2021.114056

F. Cakmak Bolat, S. Basaran, and S. Sivrioglu, “Piezoelectric and electromagnetic hybrid energy harvesting with lowfrequency vibrations of an aerodynamic profile under the air effect,” Mechanical Systems and Signal Processing, Vol. 133, p. 106246, Aug. 2019, https://doi.org/10.1016/j.ymssp.2019.106246

J.T. Lin and B. Alphenaar, “Enhancement of energy harvested from a random vibration source by magnetic coupling of a piezoelectric cantilever,” Journal of Intelligent Material Systems and Structures, Vol. 21, No. 13, pp. 1337–1341, Sep. 2010, https://doi.org/10.1177/1045389x09355662

M. Ferrari, V. Ferrari, M. Guizzetti, B. Andò, S. Baglio, and C. Trigona, “Improved energy harvesting from wideband vibrations by nonlinear piezoelectric converters,” Sensors and Actuators A: Physical, Vol. 162, No. 2, pp. 425–431, Aug. 2010, https://doi.org/10.1016/j.sna.2010.05.022

M. F. Daqaq, “Response of unimodal duffingtype harvesters to random forced excitations,” Journal of Sound and Vibration, Vol. 329, No. 18, pp. 3621–3631, Aug. 2010, https://doi.org/10.1016/j.jsv.2010.04.002

M. F. Daqaq, “Transduction of a bistable inductive generator driven by white and exponentially correlated Gaussian noise,” Journal of Sound and Vibration, Vol. 330, No. 11, pp. 2554–2564, May 2011, https://doi.org/10.1016/j.jsv.2010.12.005

W.A. Jiang and L.Q. Chen, “An equivalent linearization technique for nonlinear piezoelectric energy harvesters under Gaussian white noise,” Communications in Nonlinear Science and Numerical Simulation, Vol. 19, No. 8, pp. 2897–2904, Aug. 2014, https://doi.org/10.1016/j.cnsns.2013.12.037

W.A. Jiang and L.Q. Chen, “Stochastic averaging of energy harvesting systems,” International Journal of NonLinear Mechanics, Vol. 85, pp. 174–187, Oct. 2016, https://doi.org/10.1016/j.ijnonlinmec.2016.07.002

W.A. Jiang and L.Q. Chen, “Stochastic averaging based on generalized harmonic functions for energy harvesting systems,” Journal of Sound and Vibration, Vol. 377, pp. 264–283, Sep. 2016, https://doi.org/10.1016/j.jsv.2016.05.012

Y. Xu, R. Gu, H. Zhang, W. Xu, and J. Duan, “Stochastic bifurcations in a bistable DuffingVan der Pol oscillator with colored noise,” Physical Review E, Vol. 83, No. 5, p. 056215, May 2011, https://doi.org/10.1103/physreve.83.056215

B. Li, K. Hu, G. Jin, Y. Song, and G. Ge, “Response of cantilever model with inertia nonlinearity under transverse basal Gaussian colored noise excitation,” Mathematical Problems in Engineering, Vol. 2021, pp. 1–9, Jan. 2021, https://doi.org/10.1155/2021/1823596

G. Ge, Z. Li, Q. Gao, and J. Duan, “A stochastic averaging method on the strongly nonlinear DuffingRayleigh oscillator under Gaussian colored noise excitation,” Journal of Vibroengineering, Vol. 18, No. 7, pp. 4766–4775, Nov. 2016, https://doi.org/10.21595/jve.2016.17011

R. L. Honeycutt, “Stochastic RungeKutta algorithms. II. Colored noise,” Physical Review A, Vol. 45, No. 2, pp. 604–610, Jan. 1992, https://doi.org/10.1103/physreva.45.604
About this article
The authors gratefully acknowledge the support of the Program for Innovative Research Team in University of Tianjin (No. TD 135037,6002030152010107), the Natural Science Foundation of China (NSFC) through Grant Nos. 11402186, the Tianjin Research Program of Application Foundation and Advanced Technology through Grant No. 14JCQNJC05600.