A time-delayed feedback control is applied to a nonlinear piezoelectric energy harvesting system excited by additive and multiplicative Gaussian white noises to improve its energy harvesting performance. An equivalent decoupling system can be obtained by using a variable transformation. Based on the standard stochastic averaging method, the Fokker–Plank–Kolmogorov equation and the stationary probability density functions of the amplitude, displacement, and velocity of the harvester are obtained, respectively. In addition, the approximate expressions of mean square electric voltage and the mean extracted output power are derived. Finally, the paper explores the influences of parameters on the mean square electric voltage. The results show that noise intensity, time delay, feedback strength, time constant ratio, and coupling coefficients have great influences on the mean square electric voltage. The accuracy of the theoretical method is verified by the Monte Carlo simulation.
National Natural Science Foundation of China1167223211772254119722881. Introduction
Due to the development of wireless sensors towards low energy consumption and miniaturization, conventional power supply devices are no longer adapted to technological innovation [1, 2]. Nowadays, how to improve the efficiency of an environmental vibration energy harvesting system is becoming a hot research field.
In the recent years, vibration energy harvesters based on the piezoelectric effect have become one of the most effective energy harvesting methods due to their advantages of simple structure, no electromagnetic interference, easy manufacture, and miniaturization [3, 4]. As the first design scheme adopted, a linear vibration energy harvester relies on the resonance phenomenon and has a narrow effective frequency band. Only when the excitation frequency matches its basic natural frequency can it work effectively. A slight perturbation of the basic frequency of a linear harvester to the excitation frequency causes a sharp drop in the energy output [5–7]. Therefore, the linear harvesters based on the resonance principle are not suitable for obtaining energy from broadband and stochastic excitation sources, both of which are typical excitation types in the real world. To extend the harvester’s bandwidth, some recent solutions call for utilizing energy harvesters with stiffness-type nonlinearities. [8–13]. Cottone et al. [14] proposed that the performance of the nonlinear oscillators under stochastic excitation is superior to the linear oscillators through numerical and experimental results. McInnes et al. [15] established the theoretical model of an energy harvesting system based on stochastic resonance and found that its performance was significantly improved. Huang et al. [16] analyzed the dynamical characteristics of the nonlinear vibration energy harvester with an uncertain parameter. Erturk et al. [17, 18] investigated the stochastic response of a bistable energy harvester under harmonic excitation through numerical calculation and experiment. Daqaq [19] explored the stochastic response of the bistable energy harvesting system driven by white noise and exponential correlated Gaussian noise. However, when the bistable energy harvesters are driven by low level excitation, it will not be able to conduct high-energy interwell oscillation. Zhou et al. [20–23] designed a tristable energy harvester, and their results demonstrated that compared to the bistable energy harvester, tri-stable energy harvesters have a smaller high-energy motion threshold and can efficiently collect energy over a larger low-frequency range. Panyam and Daqaq [24] explored the responses of the tristable energy harvester to harmonic excitation, and the results showed that the tristable energy harvesters have higher output power than the bistable energy harvesters. Huang et al. [25] investigated the resonance mechanism of nonlinear vibrational multistable energy harvesters under narrow-band stochastic parametric excitations.
In fact, vibration energy harvesters are inevitably affected by environmental excitation. Therefore, the researchers investigated the performance of the nonlinear energy harvesting system under stochastic excitation [26–28]. Sun et al. [26] motivated to study stochastic dynamics in a fractional-order system subjected to Gaussian noise. Jin et al. [27] established the semianalytical solution of random response of the nonlinear vibration energy harvester under Gaussian white noise excitation by using generalized harmonic transformation and equivalent nonlinearization. Xiao and Jin [29] explored the energy harvesting of a monostable duffing type harvester with piezoelectric coupling under correlated multiplicative and additive white noise.
Some recent studies have shown that deliberately introducing delay into a continuous-time control dynamical system yields favorable results. In distributed generation systems, the introduction of time delays of a communication network is to ensure the precise voltage restoration and optimal load sharing. [30–32]. Belhaq and Hamdi [33] investigated a delayed van der Pol oscillator with time-varying delay amplitude coupled to an electromagnetic energy harvesting device. Zhang et al. [34] studied the fractional-delay system in the presence of both the colored noise and delayed feedback. Xiao et al. [35, 36] found that the coupling delay provides a strategy to control the oscillating dynamics in fractional-order oscillators. The robust distributed control algorithms consider time-delay communications and noise communication disturbances of a communication network for distributed generation systems.
In general, due to the time required for signal propagation, each electronic equipment inevitably has time delay. However, time delays in both controllers and actuators may deteriorate the control performance or even cause the instability of the piezoelectric energy harvesting system. The appropriate choices of feedback strength and time delays may improve the performance of dynamical systems. However, there are still rarely nonlinear analyses into the energy harvesting enhancement with time-delayed feedback control. Guo et al. [37] proposed a time-delayed feedback control to improve collection performance of the multiple attractors wind-induced vibration energy harvester system. The results show that time-delayed feedback is a powerful tool for achieving a wide range of operating regimes and enhancing amplitude-frequency characteristics and controlling the stochastic or deterministic dynamics of nonlinear systems. However, time-delayed feedback control of the nonlinear piezoelectric energy harvesting system [38, 39] may influence the harvesting performance and the power generation under different excitation conditions, but it has not yet been explored. Therefore, it is necessary to deeply analyze the influences of time-delayed feedback control on the energy harvesting system.
In this paper, the aimed system is a nonlinear piezoelectric energy harvesting system containing time-delayed feedback control in the presence of additive and multiplicative Gaussian white noise. We devote ourselves to studying the stochastic responses induced by time-delayed feedback in the energy harvesting system. It has been found that not only the noise time, constant ratio, and coupling coefficients but also the time delay and the feedback strength can increase the harvested energy from vibrations. The main work and the results are as follows. The piezoelectric energy harvesting system and its equivalent system are given in Section 2. The stationary probability density functions (PDFs) of displacement and velocity are calculated by using the stochastic average method, and the approximate expressions of the mean square electric voltage (MSEV) are derived in Section 3. In addition, the influences of parameters on the MSEV are discussed in detail. Finally, conclusions are drawn in Section 5.
2. Model and the Equivalent System2.1. Model
A simplified model of a piezoelectric energy harvester is shown in Figure 1 [38, 39]. Through theoretical analysis and experimental measurement, the control electromechanical model of the piezoelectric energy harvesting system is deduced:(1)mX¯¨+cX¯˙+δ1X¯+δ3X¯3+θV¯˙=−mX¯¨b,(2)CpV¯˙+V¯R=θX¯˙,where X¯ is the displacement of mass m and c, δ1, and δ3 are the damping coefficient and the linear and nonlinear stiffness coefficients, respectively. θ is the electromechanical coupling coefficient, Cp is the piezoelectric capacitance, and V is the voltage measured across the equivalent resistive load R.
A simplified model of a piezoelectric energy harvester.
For simplicity, a set of nondimensionless parameters is introduced as(3)X=X¯lc,t=ω1τ,ω1=δ1m,λ=δ3lc2δ1,V=CplcθV¯,ξ=c2δ1m,κ=θ2Cpk1,α=1RCpω1.
The nonlinear piezoelectric energy harvesting system is dimensionless as follows:(4)X¨+2ξX˙+X+λX3+κV=−X¨b,(5)V˙+αV=X˙.
On the other hand, time-delayed feedback control Ft=uX˙t−τ−X˙t is proposed [40]. Meanwhile, the energy harvesting system is inevitably affected by the external environment noise, and the additive and multiplicative Gaussian white noise excitation can more accurately replace the harmonic excitation as the stochastic noise source.
Therefore, the corresponding model with time-delayed feedback control under additive and multiplicative Gaussian white noise is(6)X¨+2ξX˙+X+λX3+κV=ξ1t+Xξ2t+uX˙τ−X˙,(7)V˙+αV=X˙,where 2ξ is the nondimensional linear damping coefficient. κ is the piezoelectric coupling coefficient in the mechanical equation. λ is the nondimensional nonlinear stiffness coefficient. α is time constant ratio. u is the feedback strength, and τ is the time delay. ξ1t and ξ2t denote Gaussian white noise which have the following properties:(8)ξit=0,Rijτ=ξitξjt+τ=2Dijδτ,i=j,0,i≠j,i,j=1,2,where D11 and D22 are the Gaussian white noise intensities of ξ1t and ξ2t, respectively.
2.2. The Equivalent System
Supposing that the noise excitation intensity, damping coefficient, and cubic nonlinear term coefficient are small, the responses of equation (6) can be regarded as stochastic disturbance of periodic solutions about a conservative linear system. Therefore, assume that the solution of system (5) has the following form [41, 42]:(9)Xt=AtΦ,X˙t=−Atω0sinΦ,where Φt=ω0t+θt, At is the amplitude, and Φt is the phase angle.
Substituting equation (9) into equation (6) yields [43](10)Vt=Cte−αt+Aω0α2+ω02ω0cosΦ−αsinΦ.
The stationary part of equation (10) is(11)Vt=Aω0α2+ω02ω0cosΦ−αsinΦ=ω02α2+ω02Xt+αα2+ω02X˙t.
As(12)ω0cosΦ−αsinΦ=α2+ω02cosΦ+tan−1αω0,
the amplitude of the stationary voltage can be obtained as(13)L=Aω0α2+ω02.
Substitution of equations (11) and (14) into equation (6) yields a modified decoupling equation:(15)X¨+SX˙+ω2X+λX3=ξ1t+Xξ2t,where ω2=1+κω02/α2+ω02−uω0sinω0τ and S=2ξ+κα/α2+ω02−ucosω0τ+u.
3. Theoretical Analysis
Using the stochastic average method, suppose the responses of equation (15) can be expressed as follows [44]:(16)Xt=AtcosΦ,X˙t=−AtωsinΦ,Φt=ωt+Θt.
Substituting equation (16) into equation (15) yields(17)A˙=1ω−AωSsinΦ+λA3cos3ΦsinΦ−1ωξ1tsinΦ−1ωξ2tAcosΦsinΦ,(18)Θ˙=cosΦAω−AωSsinΦ+λA3cos3Φ−1Aωξ1tsinΦ−1Aωξ2tAcosΦsinΦ.
By adopting the standard stochastic averaging method, the stochastic differential equation (in the Ito sense) for the instantaneous amplitude At is determined by(19)dA=mAdt+σAdWt,where Wtis an independent normalized Wiener process. A,Θ can be approximated as the two-dimensional stochastic diffusion process. The drift coefficient and diffusion coefficient of the Fokker–Plank–Kolmogorov (FPK) equation can be obtained as follows:(20)mA=−ASsin2Φ+λω2A3cos3ΦsinΦ+1Aω2D11cos2Φ−I11ωcosΦsinΦ+D22A2cos4Φ−A2cos3ΦsinΦI222ωΦ,(21)σ2A=2D11ω2sin2Φ+2D22A2cos2Φω2sin2ΦΦ,where Φ=limT⟶∞1/T∫0Tdt=1/2π∫02πdΦ. Sω and Iω have the following form:(22)Sω=1π∫−∞0cosωtRτdτ,Iω=1π∫−∞0sinωtRτdτ.
Obviously, the amplitude A does not depend on the phase Θ. In conclusion, the FPK equation of amplitude A yields(23)∂p∂t=−∂∂AmAp+12∂2∂A2σ2Ap.
In the stationary sense, equation (23) can be solved to yield the stationary PDFs for amplitude(24)pA=Cσ2Aexp∫0Amuσ2udu,in which C is a normalization factor.
The joint PDF of the amplitude and phase can be retrieved as(25)pA,Θ=pApΘ=12πpA.
According to the transformation from A,Θ to the X,X˙, we can obtain the expression of the joint PDF of X and X˙ as follows:(26)pX,X˙=1ω1ApA,Θ=12πω1ApA.
The corresponding marginal stationary PDFs of displacement X and velocity X˙ are shown in equations (27) and (28), respectively.(27)pX=∫−∞+∞pX,X˙dX˙,(28)pX˙=∫−∞+∞pX,X˙dX.
Thus, the expected value of the MSEV and the mean extracted output power are(29)V2t=ω02α2+ω02A2=ω02α2+ω02∫0+∞A2pAdA,(30)P=αkV2t.
4. Numerical Simulation
In this section, we mainly focus on the influences of the parameters on the stochastic responses of the piezoelectric energy harvesting system. To verify the validity of the analysis results, equations (6) and (7) will be directly simulated by Monte Carlo simulations. Let α=0.05, λ=0.05, and the initial conditions are set as the static equilibrium position X0=0,X˙0=0,V0=0. The Monte Carlo simulations are selected as the reference to evaluate the accuracy of the analytical results.
4.1. The Influences of Parameters on PDFs
When τ=0.5, k=0.3, and u=0.05, the stationary PDFs of amplitude, displacement, and velocity under different noise intensities are investigated, as in Figure 2. Specifically, when Gaussian white noise intensity increases continuously, the stationary PDFs peak values decrease and shift to the right, which indicates that the system leads to a large response and gives out a high harvested energy. It can be concluded that the larger excitation intensity plays an important role in the stationary responses of the nonlinear piezoelectric energy harvester. The results are consistent with those of numerical simulation, which can prove the validity of theoretical analysis.
The stationary PDFs under different noise intensities: (a) amplitude, (b) displacement, and (c) velocity. Theoretical analysis results (-) and Monte Carlo simulation results (∗).
While fixing parameters D11=0.02, D22=0.02, and τ=0.5, the stationary PDFs of the amplitude, displacement, and the velocity for different feedback strengths are demonstrated in Figure 3. As the feedback strength increases from negative to positive, the shape of the stationary PDFs become steeper gradually, which results in the probability of the system stabilizing at the equilibrium position increasing. It can be clearly explained that under this set of parameters, negative feedback intensity is more conducive to the energy harvester. According to Figures 2 and 3, we can conclude that the results obtained by the theoretical method are consistent well with those from Monte Carlo simulations.
The stationary PDFs under feedback strength. (a) Amplitude, (b) displacement, and (c) velocity. Theoretical analysis results (-) and Monte Carlo simulation results (∗).
4.2. The Influences of Parameters on the MSEV
In order to completely clarify the influences of the parameters on the performance of the piezoelectric energy harvesters, noise intensity, feedback strength, and time delay are successively introduced to investigate their influences on MSEVs in the following analysis. By taking u=−0.05, k=0.5, and τ=0.5, the effects of additive and multiplicative Gaussian white noise intensity on the MSEV are illustrated in Figures 4(a) and 4(b), respectively. It is found that with the increase of the additive and multiplicative noise intensity, MSEVs increase monotonically. The results show that additive noise plays a leading role in energy harvesting. The energy extracted from multiplicative noise by the piezoelectric energy harvester is relatively small, but it still makes sense when the noise intensity is large. Therefore, appropriately increasing the additive and multiplicative excitation intensity is a feasible method to increase the output voltage.
The influences of noise intensity on the MSEV. Theoretical analysis results (solid round) and Monte Carlo simulation results (∗).
The variation of the MSEV with the time delay and feedback strength is illustrated in Figure 5. It is seen that for any time delay τ, the MSEV V2 increases monotonically with the decrease of the feedback strength u. Nevertheless, the influences of time delay on the MSEV are nonmonotone.
The influences of time delay and feedback strength on the MSEV.
To illustrate the influences of the feedback strength on MSEV, Figure 6(a) shows the changing of the MSEV with the feedback strength u=0.01 and u=−0.01. For D11=0.02, D22=0.02, and k=0.5, the variation of the MSEV V2 over the time delay τ is investigated. The results show that MSEV of the energy harvesting system without delay feedback system can be stabilized at 0.51. When u=−0.01, the time delay τ reaches the critical value of 3.14, and the MSEV of the system reaches the maximum and then decreases monotonically. However, when the feedback strength u=0.01, the curve of the MSEV changing with time delay τ seems to be completely opposite. For τ<3.14, the MSEV decreases monotonically as the time delay increases. When τ>3.14, the MSEV increases with the continuous increase of time delay. Furthermore, negative feedback strength is more favorable to the energy harvesting system than positive feedback strength. Therefore, the harvesting efficiency of the energy harvesting system can be improved by selecting the optimal control parameters u and τ.
The influences of time delay and feedback strength on the MSEV. Theoretical analysis results (solid line) and Monte Carlo simulation results (∗).
By taking D11=0.02, D22=0.02, and τ=1, the dependence of the MSEV V2on feedback strength u and time constant ratio α is depicted in Figure 7(a). Specifically, the MSEV decreases rapidly as the time constant ratio α increases for a negative feedback strength u=−0.02. However, the MSEV decreases slowly with the increase of the time constant ratio for a positive feedback strength u=0.02. As the time constant ratio continues to increase, both the positive and negative feedback strengths keep the MSEV in a small range. Through the abovementioned analysis, a smaller time constant ratio and the negative feedback control will make the energy harvesting system induce larger MSEV.
The influences of time delay and feedback strength on the MSEV. Theoretical analysis results (solid line) and Monte Carlo simulation results (∗).
The influences of the MSEV V2on feedback strength u and the piezoelectric coupling coefficient κ are shown in Figure 8. The MSEV almost decreases proportionally with the increases of κ at any feedback strength u in Figure 8(b). This means a smaller κ and the negative feedback control will increase the production of output power.
The influences of piezoelectric coupling coefficient and feedback strength on the MSEV. Theoretical analysis results (solid line) and Monte Carlo simulation results (∗).
5. Conclusions
In this paper, the time-delayed feedback control is utilized to improve the efficiency of a nonlinear piezoelectric energy harvesting system under the influences of the additive and multiplicative Gaussian white noise excitation. The standard stochastic averaging method has been applied to derive the stationary PDFs of amplitude, displacement, and velocity. Then, the expression of MSEV is deduced via the approximate relationship between voltage and mechanical states. Furthermore, the influences of the parameters on the MSEV are investigated. The results show that the efficiency of the harvester can be improved by selecting the optimal time delay and negative feedback strength when the other parameters unaltered. Moreover, the increases of noise intensity will lead to the increases of MSEV. Besides, we find that time constant ratio and coupling coefficients will reduce the MSEV. The analysis results are almost consistent with the Monte Carlo simulation results of the original system, which verifies the rationality and validity of the theoretical method.
Compared with the nonlinear monostable, the multistable piezoelectric energy harvesting systems have attracted much attention due to the activation of the large-orbit interwell snap-through mechanism. Therefore, we will focus on the bistable energy harvesting system with the time-delayed feedback control next steps.
Data Availability
The data used to support the findings of the study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest concerning the publication of this article.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grant nos. 11672232, 11772254, and 11972288).
LiangC.AiJ.ZuoL.Design, fabrication, simulation and testing of an ocean wave energy converter with mechanical motion rectifier201713619020010.1016/j.oceaneng.2017.03.0242-s2.0-85015653570KuangY.ZhuM.Characterisation of a knee-joint energy harvester powering a wireless communication sensing node201625505501310.1088/0964-1726/25/5/0550132-s2.0-84968677725ElvinN.ErturkA.2013Berlin, GermanySpringerErturkA.InmanD. J.2011Hoboken, NJ, USAJohn Wiley & SonsRoundyS.WrightP. K.A piezoelectric vibration based generator for wireless electronics20041351131114210.1088/0964-1726/13/5/0182-s2.0-5744241231DutoitN. E.WardleB. L.KimS. G.Design considerations for memsscale piezoelectric mechanical vibration energy harvesters200571112116010.1080/105845805909645742-s2.0-33644749219ErturkA.InmanD. J.A distributed parameter electromechanical model for cantilevered piezoelectric energy harvesters2004130404100210.1115/1.28904022-s2.0-44349192514MannB. P.SimsN. D.Energy harvesting from the nonlinear oscillations of magnetic levitation20093191-251553010.1016/j.jsv.2008.06.0112-s2.0-56749117734MannB. P.OwensB. A.Investigations of a nonlinear energy harvester with a bistable potential well201032991215122610.1016/j.jsv.2009.11.0342-s2.0-74149090793MbongT. D.SieweM. S.TchawouaC.Controllable parametric excitation effect on linear and nonlinear vibrational resonances in the dynamics of a buckled beam20185437738810.1016/j.cnsns.2017.06.0192-s2.0-85020752594MallickD.AmannA.RoyS.Surfing the high energy output branch of nonlinear energy harvesters20161171919770110.1103/physrevlett.117.1977012-s2.0-84994591094ChenL. Q.JiangW. A.PanyamM.DaqaqM. F.A broadband internally resonant vibratory energy harvester2016138606100710.1115/1.40342532-s2.0-84984653325SebaldG.KuwanoH.GuyomarD.DucharneB.Experimental duffing oscillator for broadband piezoelectric energy harvesting2011201010200110.1088/0964-1726/20/10/1020012-s2.0-80053506652CottoneF.VoccaH.GammaitoniL.Nonlinear energy harvesting2009102808060110.1103/physrevlett.102.0806012-s2.0-61449128189McInnesC. R.GormanD. G.CartmellM. P.Enhanced vibrational energy harvesting using nonlinear stochastic resonance20083184-565566210.1016/j.jsv.2008.07.0172-s2.0-55349089888HuangD. M.ZhouS. X.HanQ.LitakG.Response analysis of the nonlinear vibration energy harvester with an uncertain parameter201916239340710.1177/1464419319893211ErturkA.HoffmannJ.InmanD. J.A piezomagnetoelastic structure for broadband vibration energy harvesting2009942525410210.1063/1.31598152-s2.0-67649482443ErturkA.InmanD. J.Broadband piezoelectric power generation on high-energy orbits of the bistable Duffing oscillator with electromechanical coupling2011330102339235310.1016/j.jsv.2010.11.0182-s2.0-79952438592DaqaqM. F.Transduction of a bistable inductive generator driven by white and exponentially correlated Gaussian noise2011330112554256410.1016/j.jsv.2010.12.0052-s2.0-79952539208ZhouS.CaoJ.InmanD. J.LinJ.LiuS.WangZ.Broadband tristable energy harvester: modeling and experiment verification2014133333910.1016/j.apenergy.2014.07.0772-s2.0-84905638500ZhouS.ZuoL.Nonlinear dynamic analysis of asymmetric tristable energy harvesters for enhanced energy harvesting20186127128410.1016/j.cnsns.2018.02.0172-s2.0-85042374468CaoJ.ZhouS.WangW.LinJ.Influence of potential well depth on nonlinear tristable energy harvesting20151061717390310.1063/1.49195322-s2.0-84929486139ZhouS.CaoJ.InmanD. J.LinJ.LiD.Harmonic balance analysis of nonlinear tristable energy harvesters for performance enhancement201637322323510.1016/j.jsv.2016.03.0172-s2.0-84979468143PanyamM.DaqaqM. F.Characterizing the effective bandwidth of tri-stable energy harvesters201738633635810.1016/j.jsv.2016.09.0222-s2.0-84992724029HuangD. M.ZhouS. X.YangZ. C.Resonance mechanism of nonlinear vibrational multistable energy harvesters under narrow-band stochastic parametric excitations2019201920105014310.1155/2019/1050143SunZ.DangP.XuW.Detecting and measuring stochastic resonance in fractional-order systems via statistical complexity2019125344010.1016/j.chaos.2019.05.0152-s2.0-85065859009JinX.WangY.XuM.HuangZ.Semi-analytical solution of random response for nonlinear vibration energy harvesters201534026728210.1016/j.jsv.2014.11.0432-s2.0-84921284169SuM.XuW.ZhangY.Theoretical analysis of piezoelectric energy harvesting system with impact under random excitation202011910332210.1016/j.ijnonlinmec.2019.103322XiaoS.JinY.Response analysis of the piezoelectric energy harvester under correlated white noise20179032069208210.1007/s11071-017-3784-72-s2.0-85029713919LaiJ. G.LuX. Q.YuX. H.AntonelloM.HongH.Distributed voltage regulation for cyber-physical microgrids with coupling delays and slow switching topologies201950110011010.1109/TSMC.2019.2924612LaiJ. G.LuX. Q.AntonelloM.LiuG. P.Stochastic distributed pinning control for Co-Multi-Inverter networks with A virtual leader201910.1109/TCSII.2019.2950764LaiJ. G.LuX. Q.YuX. H.AntonelloM.Stochastic distributed secondary control for ac microgrids via event-triggered communication20201142746275910.1109/TSG.2020.2966691BelhaqM.HamdiM.Energy harvesting from quasi-periodic vibrations20168642193220510.1007/s11071-016-2668-62-s2.0-84975748613ZhangJ.SunZ.YangX.XuW.Controlling bifurcations in fractional-delay systems with colored noise20182811185013710.1142/s02181274185013772-s2.0-85055655813XiaoR.SunZ.YangX.XuW.Amplitude death islands in globally delay-coupled fractional-order oscillators20199532093210210.1007/s11071-018-4678-z2-s2.0-85057747000XiaoR.SunZ.YangX.XuW.Emergence of death islands in fractional-order oscillators via delayed coupling20196916817510.1016/j.cnsns.2018.09.0152-s2.0-85054007356GuoQ.SunZ. K.ZhangY.XuW.Time-delayed feedback control in the multiple attractors wind-induced vibration energy harvesting system2019201911797382310.1155/2019/79738232-s2.0-85060184567XuM.JinX.WangY.HuangZ.Stochastic averaging for nonlinear vibration energy harvesting system20147821451145910.1007/s11071-014-1527-62-s2.0-84911007606DaqaqM. F.On intentional introduction of stiffness nonlinearities for energy harvesting under white Gaussian excitations20126931063107910.1007/s11071-012-0327-02-s2.0-84863962204PyragasK.Continuous control of chaos by self-controlling feedback1992170642142810.1016/0375-9601(92)90745-82-s2.0-34250315551ZhuW. Q.HuangZ. L.SuzukiY.Response and stability of strongly non-linear oscillators under wide-band random excitation20013681235125010.1016/s0020-7462(00)00093-72-s2.0-0035422261RobertsJ. B.SpanosP. D.Stochastic averaging: an approximate method of solving random vibration problems198621211113410.1016/0020-7462(86)90025-92-s2.0-0022563269JiangW.-A.ChenL.-Q.Stochastic averaging of energy harvesting systems20168517418710.1016/j.ijnonlinmec.2016.07.0022-s2.0-84978696283ZhuW. Q.HuangZ. L.YangY. G.Stochastic averaging of quasi-integrable Hamiltonian systems199764497598410.1115/1.27890092-s2.0-0000761056