Published: 15 May 2014

Stochastic bifurcation characteristics of cantilevered piezoelectric energy harvester

Jia Xu1
Ming-Yi Luan2
Zhi-Wen Zhu3
Kang-Kang Guo4
1Tianjin Key Laboratory of Nonlinear Dynamics and Chaos Control, 92 Weijin Road, Tianjin, 300072, P. R. China
1, 2, 3, 4School of Mechanical Engineering, Tianjin University, 92 Weijin Road, Tianjin, 300072, P. R. China
Corresponding Author:
Kang-Kang Guo
Views 39
Reads 16
Downloads 1165

Abstract

Stochastic bifurcation characteristics of cantilevered piezoelectric energy harvester were studied in this paper. Von de Pol differencial item was introduced to interpret the hysteretic phenomena of piezoelectric ceramics, and then the nonlinear dynamic model of piezoelectric cantilever beam subjected to axial stochastic excitation was developed. The stochastic stability of the system was analyzed, and the steady-state probability density function and the joint probability density function of the dynamic response of the system were obtained, and then the conditions of stochastic Hopf bifurcation were analyzed. Numerical simulation shows that stochastic Hopf bifurcation appears when bifurcation parameter varies, which can increase vibration amplitude of cantilever beam system and improve the efficiency of piezoelectric energy harvester. Finally, the theoretical and numerical results were proved by experiments. The results of this paper are helpful to application of cantilevered piezoelectric energy harvester in engineering fields.

1. Introduction

Piezoelectric ceramics is a kind of smart material. It can be used to convert mechanical energy into electrical energy, which is known as piezoelectric effect. Based on this effect, piezoelectric energy harvester can be designed to gather vibration energy of structures. Compared with other kinds of power generation systems, piezoelectric energy harvester has many advantages, such as small size, high electro-mechanical conversion efficiency, long service life, and low cost, which cause it be applied as green energy widely.

Many scholars studied piezoelectric energy harvester. DuToit designed MEMS-scale piezoelectric mechanical vibration energy harvesters firstly [1]. Erturk developed the mechanical model of cantilevered piezoelectric vibration energy harvesters [2]. Priya proposed the criterion for material selection in design of bulk piezoelectric energy harvesters [3]. Liao studied parameters optimization and power characteristics of piezoelectric energy harvesters with RC circuit [4]. Although many advances were obtained, the modeling problem limits the application of cantilever piezoelectric energy harvester in industry fields. In order to optimize piezoelectric energy harvester effectively, it is necessary to build a model in high accuracy to describe the nonlinear characteristics of piezoelectric energy harvester.

Due to the hysteretic characteristics of piezoelectric ceramics, most of the piezoelectric models were shown as equations with subsection function or double integral function, which were hard to be analyzed in theory [5-9]. Usually, research results could only be obtained by numerical or experiment method [10, 11]. In this paper, hysteretic nonlinear theory was introduced to develop a new kind of continuous piezoelectric model, and then stochastic bifurcation characteristics of cantilevered piezoelectric energy harvester were analyzed.

2. Hysteretic nonlinear model of piezoelectric energy harvester

The voltage-displacement curve of piezoelectric ceramics was shown in Fig. 1. Obviously, there is hysteretic nonlinearity in piezoelectric ceramics.

Fig. 1The displacement-voltage curves of piezoelectric ceramics

The displacement-voltage curves of piezoelectric ceramics

In this paper, Von del Pol hysteretic model was introduced to describe the hysteretic nonlinear characteristics of piezoelectric ceramics. The initial Von del Pol hysteretic model describes hysteretic loop which is symmetrical about the initial point (0, 0). It can be shown as follows:

1
y=fx=f0x+a1-xb2x˙,

where f0x is skeleton curve of hysteretic loop and usually expressed in polynomial function, a and b are coefficients which determine the difference between the skeleton curve and the real curve. The essence of Von del Pol item a1-xb2x˙ are two parabolic lines which are symmetrical about the original point (0, 0).

Supposing the displacement-voltage curves of piezoelectric ceramics is symmetrical about the point G(V0,x0), it can be shown as follows:

2
V-V0=b1x-x0+b2x-x03+b31-x-x0b42x˙,

where V is voltage, x is displacement, bi(i=1, 2, 3, 4) are coefficients, skeleton curve is chosen as f0(x)=b1x+b2x3.

b4=x0 since the loading curve have the same value as the unloading curve when x=0, and σ0-b1x0-b2x03=0 because the initial voltage of piezoelectric ceramics is zero. Thus, Eq. (2) can be rewritten as follows:

3
V=a1x+a2x2+a3x3+a4x-a5x2x˙,

where a1=b1+3b2x02, a2=-3b2x0, a3=b2, a4=2b3b4, a5=b3b42.

The structure of cantilevered piezoelectric energy harvester was shown in Fig. 2. It can be regarded as cantilever composite beam. The mechanical model of cantilevered piezoelectric energy harvester was shown in Fig. 3, where the thickness of adhesive layer was ignored.

According to Hamilton's principle, we obtained:

4
t1t2[δ(T-U+We)+δW]dt=0,

where T is kinetic energy of structure, U is potential energy of structure, We is electric energy of piezoelectric ceramics, W is work done by external force:

5
T=12Vs-ρs-u˙2dVs-+12Vp-ρp-u˙2dVp-,
6
U=12Vs-σs-εs-dVs-+12Vp-σp-εp-dVp-,
7
We=12Vp-E3D3dVp-,
8
δW=0LδuFdx+δul,tF+δϕq-,

where ρi(i=s-, p-) is density of substrate and PZT piezoelectric ceramics, where the subscript s- is substrate, p- is PZT piezoelectric ceramics; Vi is volume, σi is stress, εi is strain; u=u(x,t) is deflection of composite cantilever beam, E3 is electric field intensity, D3 is electrostrictive displacement, F=eς(t) is axial stochastic excitation, e is intensity of stochastic excitation, ς(t) is Gauss white nosie whose mean is zero and intensity is 2D, ϕ is electric potential, q- is electric charge.

Fig. 2The structure of cantilevered piezoelectric energy harvester

The structure of cantilevered piezoelectric energy harvester

Fig. 3The mechanical model of cantilevered piezoelectric energy harvester

The mechanical model of cantilevered piezoelectric energy harvester

The dynamic equation of piezoelectric cantilever beam subjected to axial stochastic excitation can be obtained as follows:

9
x¨+2ζmx˙+kx-ϑ2V+ϑ2α1xV+12ϑ2α2V2-α3x2+α4x3=exςt,

where ζm is damping, k is stiffness, ϑ and αi(i= 1, 2, 3, 4) are coefficients.

Substituting Eq. (3) into Eq. (9), we obtained:

10
x¨+A1xx˙+A2xx˙2-B1x+B2x2+B3x3+B4x4+B5x5+B6x6=exςt,

where:

11
A1x=2ζm-ϑ2a4x-a5x2+ϑ2α1a4x-a5x2x+ϑ2α2a1a4x-a5x2x
+ϑ2α2a2a4x-a5x2x2+ϑ2α2a3a4x-a5x2x3,
12
A2x=12ϑ2α2a4x-a5x22,
13
B1=k+ϑ2a1,
14
B2=ϑ2α1a1-a2+12α2a12-α3,
15
B3=ϑ2α1a2-a3+a1a2+α4,
16
B4=ϑ2α1+12α2a22+α2a1a3,
17
B5=ϑ2α2a2a3,
18
B6=12ϑ2α2a32.

3. Stochastic stability analysis of system

Let x=q, x˙=p, Eq. (10) can also be shown as follows:

19
q˙=p,p˙=-A1qp-A2qp2+B1q-B2q2-B3q3-B4q4-B5q5-B6q6+eqςt.

The Hamiltonian function of Eq. (19) can be shown as follows:

20
H=12p2-12B1q2+13B2q3+14B3q4+15B4q5+16B5q6+17B6q7.

According to the quasi-nonintegrable Hamiltonian system theory, the Hamiltonian function H(t) converges weakly in probability to an one-dimensional Ito diffusion process. The averaged Ito equation about the Hamiltonian function can be shown as follows:

21
dH=mHdt+σHdBt,

where B(t) is standard Wiener process, m(H) and σ(H) are drift and diffusion coefficients of Ito stochastic process, which can be obtained in stochastic averaging method as follows:

22
mH=-ζmHR-14ϑ2α1a4R2-18ϑ2a5+14ϑ2α2a1a4R3
+116ϑ2α2a2a5-a3a4R5+2e2DHR,
23
σ2H=22e2DHR2,

where R is the solution of the following equation:

24
-12B1R2+13B2R3+14B3R4+15B4R5+16B5R6+17B6R7=H.

Then the associated largest Lyapunov exponent of the system is:

25
λ=limt1tlnH1/2=2ϑ2α1a4-e2D4(k+ϑ2a1).

Now the local stochastic stability of the system can be discussed as follows:

1) The trivial solution H=0 is locally asymptotic stable if and only if λ< 0, which means the vibration amplitude of the system will tend to zero;

2) The trivial solution H=0 is locally asymptotic unstable if and only if λ> 0, which means the vibration amplitude of the system will tend to be large;

3) Bifurcation should appear near the trivial solution H=0 if and only if λ= 0, which means the vibration amplitude of the system will jump between the small value and the large value.

The largest Lyapunov exponent can only estimate the local stability. In this paper, the boundary classification method was used to analyze the global stability of the trivial solution of the system. Generally, the boundaries of diffusion process are singular, and the boundary classification is often determined by diffusion exponent, drift exponent and character value [12].

When H0, we obtained that: mHOH,σ2HOH2,αl=2,βl=1,cl=2ϑ2α1a4De2, where αl is diffusion exponent, βl is drift exponent, cl is character value, l is left boundary. Thus, the left boundary H=0 belongs to the first kind of singular boundary.

According to the classification for singular boundary, we obtained:

1) The left boundary H=0 is repulsively natural if cl> 1;

2) The left boundary H=0 is strictly natural if cl= 1;

3) The left boundary H=0 is attractively natural if cl< 1.

When H, we obtained that: m(H)=O(H8/7),σ2(H)=O(H9/7),αr=97,βr=87, where r is the right boundary. Thus, the right boundary H= belongs to the first kind of singular boundary. Thus, the right boundary H= is an entrance boundary.

The necessary and sufficient conditions for globally asymptotic stability of the trivial solution require that the left undary be attractively natural and the right boundary be entrance. Thus, the trivial solution H=0 is globally asymptotically stable if and only if cl< 1, which means 2ϑ2α1a4<e2D. The influence of the character value to the stability was shown in Fig. 4.

Fig. 4The influence of the character value to the stability

The influence of the character value to the stability

4. Stochastic bifurcation and simulation

The averaged FPK equation of Eq. (19) is:

26
ft=-HmHf+122σ2HfH2,

where f is probability density.

Thus, the stationary probability density function of the system is:

27
fH=A-Hηexp22De22-ζm-18ϑ2a5-14ϑ2α2a1a4H-12-14ϑ2α1a4H-1-116ϑ2α2a2a5-a3a4H1/2,

where A- is a normalization constant, η=2ϑ2α1a4-2De2De2.

The result of numerical simulation were shown in Fig. 5-8, where k=0.5, D=0.5, l=1, ξm=0.05, M=40, E=1011, A=10-4, I=10-11, a1=3.8, a2=–0.27, a3=0.014.

Fig. 5The steady-state probability density of the system when a4=2 and a5=4.5

The steady-state probability density of the system when a4=2 and a5=4.5

Fig. 6The joint probability density of the system when a4=2 and a5=4.5

The joint probability density of the system when a4=2 and a5=4.5

Fig. 7The steady-state probability density of the system when a4=3.5 and a5=4

The steady-state probability density of the system when a4=3.5 and a5=4

Fig. 8The joint probability density of the system when a4=3.5 and a5=4

The joint probability density of the system  when a4=3.5 and a5=4

From Fig. 5-8, we can see that:

1) The hysteretic nonlinear damping coefficients ai (i=4, 5) can induce stochastic Hopf bifurcation of the system. From Fig. 6 and Fig. 8, we can obviously see that there are two limit cycles in the stationary probability density, which means that there are two vibration amplitudes whose probability are both very high. Jumping phenomena between the two vibration amplitudes will appear when the conditions are changed;

2) The stationary probability density of the response of the system can be changed through adjusting the parameters ai (i=4, 5). It means that different PZT piezoelectric ceramics materials will cause different vibration amplitudes of the system since the parameters ai (i=4, 5) are determined by PZT piezoelectric ceramics. It provide a way to improve the efficiency of the energy harvester since the stationary probability density of the big vibration amplitude can be increased by choosing appropriate PZT piezoelectric ceramics.

The experimental results of cantilevered piezoelectric energy harvester were shown in Fig. 9 and Fig. 10, where the vibration amplitude was shown as output voltage of sensor. We can see that system paremeters ai (i=4, 5) can influence the vibration of the system, and jumping phenomena between the two vibration amplitudes appears when the conditions are changed.

Fig. 9The response of cantilevered piezoelectric energy harvester when a4=3.5 and a5=4

The response of cantilevered piezoelectric energy harvester when a4=3.5 and a5=4

Fig. 10The response of cantilevered piezoelectric energy harvester when a4= 2 and a5= 4.5

The response of cantilevered piezoelectric energy harvester when a4= 2 and a5= 4.5

5. Conclusions

Stochastic bifurcation characteristics of cantilevered piezoelectric energy harvester had been studied in this paper. Von de Pol differencial item was introduced to interpret the hysteretic phenomena of piezoelectric ceramics, and then the nonlinear dynamic model of piezoelectric cantilever beam subjected to axial stochastic excitation was developed. The stochastic stability of the system was analyzed, and the steady-state probability density function and the joint probability density function of the dynamic response of the system were obtained, and then the conditions of stochastic Hopf bifurcation were analyzed. Numerical simulation shows that stochastic Hopf bifurcation appears when bifurcation parameter varies, which can increase vibration amplitude of cantilever beam system and improve the efficiency of piezoelectric energy harvester. Finally, the theoretical and numerical results were proved by experiments. The results of this paper are helpful to application of cantilevered piezoelectric energy harvester in engineering fields.

References

  • DuToit N. E., Wardle B. L., Kim S. G. Design considerations for MEMS-scale piezoelectric mechanical vibration energy harvesters. Integrated Ferroelectrics, Vol. 71, Issue 1, 2005, p. 121-160.
  • Erturk S. A., Inman D. J. On mechanical modeling of cantilevered piezoelectric vibration energy harvesters. Journal of Intelligent Material Systems and Structures, Vol. 19, Issue 11, 2008, p. 1311-1325.
  • Priya P., Shashank M. Criterion for material selection in design of bulk piezoelectric energy harvesters. IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control, Vol. 57, Issue 12, 2010, p. 2610-2612.
  • Liao M. R. Y. B., SodanoH. A. Optimal parameters and power characteristics of piezoelectric energy harvesters with an RC circuit. Smart Materials and Structures, Vol. 18, Issue 4, 2009, p. 045011.
  • Hagood N. W., Chung W. H., Flotow A. V. Modelling of piezoelectric actuator dynamics for active structural control. Journal of Intelligent Material Systems and Structures, Vol. 1, Issue 3, 1990, p. 327-354.
  • Akella P., Chen X., Cheng W. Modeling and control of smart structures with bonded piezoelectric sensors and actuators. Smart Materials and Structures, Vol. 3, Issue 3, 1994, p. 344-353.
  • Zhu Z. W., Han J. X., Xu J. Modeling of piezoelectric actuator based on hysteretic nonlinear theory. Sensor Letters, Vol. 9, Issue 5, 2011, p. 2052-2057.
  • Goldfarb M., Celanovic N. Modeling piezoelectric stack actuators for control of micromanipulation. IEEE Control Systems Magazine, Vol. 17, Issue 3, 1997, p. 69-79.
  • Shen M. H. A new modeling technique for piezoelectrically actuated beams. Computers and Structures, Vol. 57, Issue 3, 1995, p. 361-363.
  • DuToit N. E., Wardle B. L. Experimental verification of models for microfabricated piezoelectric vibration energy harvesters. AIAA Journal, Vol. 45, Issue 5, 2007, p. 1126-1137.
  • Carlos D. M. J., Erturk A., Inman D. J. An electromechanical finite element model for piezoelectric energy harvester plates. Journal of Sound and Vibration, Vol. 327, Issue 1, 2009, p. 9-25.
  • Zhu W. Q., Huang Z. L. Stochastic Hopf bifurcation of quasi non-integrable Hamiltonian systems. International Journal of Non-Linear Mechanics, Vol. 34, Issue 3, 1999, p. 437-447.

About this article

Received
11 October 2013
Accepted
13 December 2013
Published
15 May 2014
Keywords
piezoelectric energy harvester
hysteretic loop
stochastic bifurcation
Acknowledgements

The authors gratefully acknowledge the support of Natural Science Foundation of China (NSFC) through grant No. 11272229 and 11302144, the Ph. D. Programs Foundation of Ministry of Education of China through grant No. 20120032120006, and Tianjin Research Program of Application Foundation and Advanced Technology through grant No. 13JCYBJC17900.