Shock and Vibration
Volume 2016 (2016), Article ID 9614842, 13 pages
http://dx.doi.org/10.1155/2016/9614842

Research Article

A Broadband Vibration-Based Energy Harvester Using an Array of Piezoelectric Beams Connected by Springs

V. Meruane and K. Pichara

Department of Mechanical Engineering, Universidad de Chile, Beauchef 851, 9370456 Santiago, Chile

Received 2 June 2015; Revised 25 August 2015; Accepted 27 August 2015

Academic Editor: Carlo Trigona

Copyright © 2016 V. Meruane and K. Pichara. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Piezoelectric cantilevered beams have been widely used as vibration-based energy harvesters. Nevertheless, these devices have a narrow frequency band and if the excitation is slightly different there is a significant drop in the level of power generated. To handle this problem, the present investigation proposes the use of an array of piezoelectric cantilevered beams connected by springs as a broadband vibration-based energy harvester. The equations for the voltage and power output of the system are derived based on the analytical solution of the piezoelectric cantilevered energy harvester with Euler-Bernoulli beam assumptions. To study the advantages and disadvantages of the proposed system, the results are compared with those of an array of disconnected beams (with no springs). The analytical model is validated with experimental measurements of three bimorph beams with and without springs. The results show that connecting the array of beams with springs allows increasing the frequency band of operation and increasing the amount of power generated.

1. Introduction

The idea of harvesting energy from the environment resulting in self-powered monitoring systems is fascinating and has attracted the attention of many researchers. This is demonstrated by the large numbers of investigations in the area [1–8]. Examples of remote sensors that could benefit from energy harvesting are weather monitoring, sea surface water monitoring, tactical military surveillance, monitoring of equipment and structures, oceanography, and volcano and seismic monitoring to name a few. For most of these applications, renewable power sources are needed for long-term remote monitoring. Some of the different power sources available are thermoelectric energy, mechanical vibrations, hydroelectric power, wind energy, solar power, and ambient radio frequency.

Mechanical vibrations are an attractive source due to their high availability in certain environments. The transduction mechanisms that can be used to convert ambient vibration into electrical energy are electromagnetic, electrostatic, and piezoelectric. Among these mechanisms, piezoelectric transduction has received great attention in the last years due to its ability to directly convert strain energy into usable electric energy [1, 9]. In piezoelectric energy harvesting the mechanical energy of vibrations is transformed by a piezoelectric material to electrical energy. When a piezoelectric material is deformed, the central molecules in the crystal become polarized and form a dipole. If the dipoles are arranged, then two surfaces of the material become positively and negatively charged. This property can be exploited to transform mechanical strain into electrical energy. Figure 1 presents a scheme of two cantilever piezoelectric power generators. The generators are submitted to base excitation, which causes the beams to vibrate, inducing an alternating voltage. The voltage generated by the piezoelectric elements can be used to power a device directly or to charge a rechargeable battery. In the case of a unimorph beam, a layer of piezoceramic material is attached to a metal layer (substructure), whereas in a bimorph configuration a metal layer is sandwiched between two piezoelectric layers. In both cases, a proof mass is attached at the end to increase the strain.

(a)

(b)

Figure 1: Cantilever piezoelectric energy configurations: (a) unimorph, (b) bimorph (parallel connection).

The main problem of conventional energy harvesters is their narrow frequency band of operation. Conventional energy harvesters, as the ones shown in Figure 1, have an operating frequency equal to their resonance frequency, and if the excitation is slightly different there is a significant drop in the level of power generated. This problem has led to the search for broadband energy harvesters. Twiefel and Westermann [10] present a survey of broadband techniques for vibration-based energy harvesting. These techniques are classified as follows:

(i)Linear generators: arrays, combined modes.(ii)Nonlinear generators: bi/multiple stable, vibroimpact, resonant tuning.(iii)Advanced electronic networks: switching networks, synthetic impedance.

According to the researchers there is no universal solution; depending on the application, one technique may be better than the others. For example, linear generators are more efficient for stochastic excitations whereas nonlinear techniques are more suitable for harmonic excitations or very low excitation frequencies.

In the case of linear generators, recent research has been focused on developing systems that exhibit multiple resonances that are closely spaced. The simplest solution is to use multiple cantilevers. The cantilevers are built with different geometries (large, tip mass or width) so each beam has a slightly different resonant frequency, thus increasing the bandwidth of the full system. A first design of a cantilever array was presented by Shahruz [11, 12]. He concluded that if the dimensions of the beams and masses are chosen appropriately, then the system performs as a broadband harvester. However, the maximal frequency band is limited and is independent of the beams dimensions. This approach has been studied extensively since then [13–15].

Another method is to use structures with two or more degrees of freedom, resulting in two or more resonant frequencies that are closely spaced [16–19]. Z. Yang and J. Yang [18] analyzed the power generated by the flexural vibrations of two piezoelectric beams connected by springs; they showed that this structure can be used as a broadband energy harvester. Nevertheless, the model presented by the authors is restricted to two beams and the extension to more than two beams is not trivial. Springs have also been used in MEMS-based energy harvesters, where it has been demonstrated that the use of nonlinear springs increases the generated power [20–22].

The research community related to piezoelectric energy harvesting is very wide and includes researchers from mechanical, electrical, materials, and civil engineering. This has led to the development of many different models for the same problem: a piezoelectric cantilevered beam. Most of these models use simplifications, such as oversimplifying of the piezoelectric coupling as viscous damping [23], using of the static deflection of the beam in a dynamic model [24, 25], no inclusion of the piezoelectric effect in the mechanical model [26], or approximate solutions as the Rayleigh-Ritz approximation [27, 28]. Erturk and Inman [29] provide a comprehensive discussion of the different simplifications found in literature to model piezoelectric beams, highlighting the common indiscretions and providing corrections and necessary clarifications for researchers from different engineering fields. The best solution can be obtained by solving the continuous electromechanical model analytically. Different approaches for solving the continuous model can be found in literature [30–33].

In this paper we present a general approach to characterize multiple piezoelectric beams connected with springs. The analytical solution of the cantilevered piezoelectric energy harvester with Euler-Bernoulli beam assumptions proposed by Erturk and Inman [33] is extended to the case of an array of piezoelectric beams connected by springs. The proposed configuration is attractive since it combines the properties of beam arrays and of systems with multiple degrees of freedom. Furthermore, the results demonstrate that connecting an array of beams with springs allows increasing the frequency band of operation and increasing the amount of power generated when compared with an array of disconnected beams.

The remainder of this work is organized as follows. Section 2 presents the structure under investigation and the methodology to compute its modes shapes and natural frequencies. Section 3 provides the derivation of the electromechanical equations of a bimorph cantilever beam. In Section 4 the equations of the output voltage and power of an array of bimorph cantilever beams connected by springs are derived. Section 5 presents analytical and experimental results for different combinations of tip masses and springs constants, comparing the results with those of an array of disconnected beams (with no springs). Finally, conclusions and forthcoming work are presented in Section 6.

2. Structure and Vibration Modes

The structure under investigation is an array of piezoelectric beams connected by springs, as illustrated in Figure 2(a). All beams are clamped at the left-end to a vertical wall. The right-end of each beam is connected to a concentrated tip mass and to its neighbors beams by springs. The inertia of the concentrated masses is neglected. All beams have identical dimensions and the only differences are the tip masses.

(a) Full system

(b) Scheme of the th beam

Figure 2: Structure consisting in an array of piezoelectric beams connected by springs.

To derive the vibration modes of this system, it is necessary to study the boundary conditions of each beam separately. Looking at the th beam, as shown in Figure 2(b), and assuming a separation-of-variables solution, the transverse displacement of the th beam can be expressed as . The function represents the mode shape and depends on the boundary conditions. In an Euler-Bernoulli beam, the general solution for is

The boundary conditions at the left-end are

which implies that

Therefore the mode shape can be written as

At the right-end, the bending moment and shear force must satisfy the following conditions:

where is Young’s modulus, is the area moment of inertia of the cross section, is the th tip mass, and is the stiffness of the th spring. Replacing that and , being the natural frequency, yields

where is the length of the beam. Replacing the definition of from (4), the following equations are obtained:

with

Equation (7) can be obtained for each of the beams in the structure. Particular cases are the first and last beams. In the case of the first beam there is no previous neighbor; thus the following equation is obtained:

whereas for the last beam there is no next neighbor and the resulting equations are

Coupling all equations yields a system of equations, which can be written as the following matrix equation:

where is a coefficient matrix and . This equation can have nonzero solution for only if the determinant of the coefficient matrix vanishes (i.e., the coefficient matrix is singular). Setting the determinant of equal to zero yields the characteristic equation that is satisfied by an infinite number of choices for , denoted by . The vector of unknowns is calculated by replacing in the coefficient matrix and then computing its null space. For each value of a different vector is obtained.

Given and , the mode shapes are given by

The undamped natural frequency associated with the th mode is given by

where is the mass per unit length. As an example, Figure 3 illustrates the vibration modes of three elastically connected beams. In this example, all beams have a mass tip of , and all springs have a stiffness of .

(a) Mode 1

(b) Mode 2

(c) Mode 3

(d) Mode 4

(e) Mode 5

(f) Mode 6

(g) Mode 7

(h) Mode 8

(i) Mode 9

Figure 3: First nine vibration modes of three elastically connected beams with identical tip masses and spring constants ().

3. Electromechanical Model of a Bimorph Piezoelectric Beam

3.1. Mechanical Equation

Let us consider the bimorph cantilever beam shown in Figure 4. The bimorph beam consists in two layers of piezoelectric material (PZT) bonded to a substructure layer. Assuming an Euler-Bernoulli beam, the equation of motion can be written as [33]

where is the transverse deflection of the beam relative to the base, is the base displacement, is the internal bending moment, is the equivalent coefficient of strain rate damping, is the viscous air damping, is the area moment of inertia of the cross section, is the mass per unit length, is the tip mass, and is the Dirac delta function.

Figure 4: Scheme of the bimorph piezoelectric beam.

The internal moment is given by the integral of the stress distribution over the cross section of the beam shown in

where is the width of the beam, is the thickness of each PZT layer, is the thickness of the substructure (see Figure 4), and is the stress. Superscripts and stand for PZT and substructure, respectively. The stress distribution across the section is given by the strain-stress relations of the substructure and PZT layers, which are given by

where is the strain, is Young’s modulus, is the piezoelectric constant, and is the electric field. The bending strain can be expressed in terms of the radius of curvature as

Since the piezoelectric layers are connected in parallel they have the same voltage, but the instantaneous electric fields are in opposite directions. The electric field can be written in terms of the voltage across the PZT, , and the thickness of the PZT layer, , by the following equation:

Replacing (16) to (19) in (15), we obtain

Equation (20) can be expressed as

where is the length of the beam and the Heaviside function. It has been assumed that the PZT layer covers the entire length of the beam. The constant represents the bending stiffness of the beam and is given by

and the term represents the coupling between the strain and the voltage and is given by

Substituting (21) in (14), we obtain the mechanical equation of motion with electrical coupling:

3.2. Electrical Equation

To derive the electrical equation let us start with the piezoelectric constitutive relation given in

where is the electric displacement and is the permittivity at constant stress. Replacing the definition of from (17) yields

with . The strain evaluated at the middle of the PZT layer is given by

where is the distance of the center of the PZT layer to the neural axis. Thus,

The electric charge is obtained by integrating the electric displacement over the area of the beam. The resulting expression is

Therefore the current generated by one PZT layer is

Since both layers are connected in parallel, the current generated by the bimorph beam is given by

If the beam is connected directly to a resistive load, the voltage is given by . Substituting this relation in (31) yields

The electrical equation with mechanical coupling is obtained by rearranging (32):

4. Voltage and Power Output of an Array of Piezoelectric Beams Connected by Springs

Equations (24) and (33) are the electromechanical equations for a single bimorph beam under transverse vibrations. Let us consider the case of multiple beams connected by springs with a parallel electrical connection, as shown in Figure 5. In that case the mechanical equation of motion is given by

where is the voltage of the th beam. Since the beams are electrically connected in parallel they have the same voltage; thus . It has been assumed that all the beams are identical except for the tip mass. It should be noted that the springs are included as boundary conditions. Therefore, they affect the vibration modes, but not the equation of motion.

Figure 5: A parallel electrical connection of an array of bimorph beams connected by springs.

The solution for the relative motion of the system can be represented by a sum of its mode shapes as

where and are the th mode shape and modal coordinate, respectively. The th mode shape can be defined by segments as

where is the th mode shape of the th beam, is the length of one beam, and is the number of beams.

Replacing (35), the mechanical equation of motion can be written as

with

Assuming a harmonic base motion, , and a harmonic response, and , (37) becomes

Rearranging (39), the amplitude of the temporal term, , can be written as follows:

For simplicity, (40) can be expressed as

with

The electrical equation needs to be derived for each beam individually. The relative motion for the th beam is given by

Replacing (43) in (30) yields

If the beams are connected directly to a resistive load, the voltage is given by ; substituting this condition equation (44) becomes

Replacing the harmonic response condition in (45) yields

Rearranging (46), we obtain

with

Replacing from (41) in (47) the following equation is obtained:

Rearranging (49) the voltage amplitude can be written as shown in

The voltage FRF is defined as the ratio of the voltage output to the base acceleration, . The instantaneous power output can be obtained from the relation . Therefore, the FRF of the power output is simply the square of the voltage FRF divided by .

5. Results

5.1. Analytical

This section presents the results of the power FRF for different number of beams and different combinations of tip masses and spring constants in the array. The geometric, material, and electromechanical parameters of the beams are given in Table 1. These are the parameters of the experimental beams used in the next section.

Table 1: Parameters of the bimorph beams.

Length of the beam, (mm)
48.2

Width of the beam, (mm)
10

Thickness of the substructure, (mm)
0.19

Thickness of the PZT, (mm)
0.24

Young’s modulus of the substructure, (GPa)
40.4

Young’s modulus of the PZT, (GPa)
55

Mass density of the substructure, (kg/m³)
1519

Mass density of the PZT, (kg/m³)
8140

Piezoelectric constant, (pm/V)
−180

Permittivity, (nF/m)
13.28

Strain rate damping, (kgm⁴/s)
1.16 × 10⁻⁶

Viscous air damping, (kg/s)
0.4

The voltage is computed by (50) and the power FRF is computed as . The value of the external load resistance is defined for each case as the optimum resistor that maximizes the power output. To study the advantages and disadvantages of the proposed system, the results are compared with those obtained by an array of disconnected beams (with no springs).

First, the case of identical beams is investigated; the results are shown in Figure 6(a). It is obtained that, in the case of identical beams, the only modes that are excited are the modes where all beams move in phase, such as modes 1, 4, and 7 in Figure 3. In this case, it is also found that the natural frequencies do not depend on the spring stiffness because all beams are in phase and the springs do not deform. Furthermore, the results are exactly the same as the ones obtained with disconnected beams, as illustrated in Figure 6(b).

(a) Connected

(b) Disconnected

Figure 6: Power FRF with an increasing number of identical bimorph beams, and .

In the case of beams with different tip masses, the results of the FRF power do depend on the springs connecting the beams, as shown in Figure 7. If the beams are disconnected and new beams with different tip masses are added to the array, then the frequency band of the device is increased obtaining a similar power amplitude for all resonant frequencies. On the other hand, if the beams are connected by spring and new beams are included in the array, then the frequency band is increased, but the amplitude related to each resonant frequency varies.

(a) Connected

(b) Disconnected

Figure 7: Power FRF with an increasing number of bimorph beams with different tip masses. Two beams: , , ; three beams: , , , ; six beams: , , , , , , .

For an array of piezoelectric beams connected by springs, the frequency band of operation and the power amplitude associated with each resonant frequency depend on the combination of tip masses and springs constants. This is more evident in Figures 8 and 9, where different combinations of springs and tip masses are studied for the case of two and three beams. In general, increasing the tip mass of the beams decreases some of the resonant frequencies and increasing the spring stiffness increases the frequency band.

(a) Varying tip masses ()

(b) Varying spring constants ()

Figure 8: Power FRF of two bimorph beams with different combinations of tip masses and spring constants.

(a) Varying tip masses ()

(b) Varying spring constants (

Figure 9: Power FRF of three bimorph beams with different combinations of tip masses and spring constants.

To compare the performances of an array of beams connected by springs and one of the disconnected beams, a case of ten beams with both situations is plotted in Figure 10. It is possible to see that the array with connected beams has a slightly wider frequency band, which increases with an increment of the springs constants. Furthermore, by comparing the areas below both curves in the frequency range 50–200 Hz, the array with connected beams provides an area between 20% and 36% larger than the array with disconnected beams. Thus, by connecting the beams, the performance of the array can be improved. It should be noted that the previous results are obtained with arbitrary values for the tip masses and spring constants. Therefore, this performance can be further improved by a proper design.

(a)

(b)

(c)

(d)

Figure 10: Power FRF with ten bimorph beams: , , , , , , , , , .

5.2. Experimental

To validate the proposed model a preliminary experimental campaign was performed. Figure 11 presents the experimental setup. It consists of a signal generator connected to a power amplifier and to an electrodynamic shaker, which produces a controlled base vibration. An array of three piezoelectric bimorph beams is mounted over the shaker; Figure 12 shows this array. The cantilever bimorph beams have a tip mass at the free-end and can be connected with springs between them. Three experiments were performed. In the first case, no springs were used and each beam was connected independently to a  kΩ resistive load; in this case the voltage generated by each beam is monitored individually. In the second case, no springs were used and the beams were electrically connected in parallel to a resistive load of  kΩ. In the last case, the beams were connected with springs between them and they were electrically connected in parallel to a resistive load of  kΩ.

Figure 11: Experimental setup.

Figure 12: Experimental bimorph beams.

The properties of the bimorph beams are the ones given in Table 1; the tip masses are , , and the spring constants are .

During the three experiments the system is subjected to a harmonic base displacement with frequencies in the range 20–70 Hz. The results for the first experiment are presented in Figure 13(a) and are compared with the ones obtained by an analytical model in Figure 13(b). Both, analytical and experimental, curves are very similar. The only difference is that the experimental power increases from beams 1 to 3, and the analytical model predicted the opposite behavior. This unexpected behavior seems to be related to a resonance of the supporting structure; in the three experimental curves, it is possible to see an increment of the power near 60 Hz, which can be explained by a resonance of the structure.

(a) Experimental

(b) Analytical

Figure 13: Power FRF for the three beams with no springs and independent electrical connections to  kΩ resistive loads: (a) experimental measurement, (b) analytical results.

Figure 14 presents the experimental and analytical results for the second and third experiments. It is possible to see a good agreement between the experimental and analytical curves, although, similarly to the results in Figure 13(a), there is an increment of the power near 60 Hz that confirms a structural resonance.

(a) Experimental

(b) Analytical

Figure 14: Power FRF for the three beams with and without springs and electrical connection in parallel to a resistive load of  kΩ: (a) experimental measurement, (b) analytical results.

Comparing the areas below both experimental curves in Figure 13(b) in the frequency range 20–70 Hz, the array with connected beams provides an area 5% larger than the array with disconnected beams, which confirms that the addition of springs can increment the generated power. These results can be improved by using softer springs, which have lower energy losses, and by solving the structural resonance problem.

6. Conclusions

This paper investigated the use of an array of bimorph beams connected by springs as a broadband vibration-based energy harvester. The equations for the voltage and power output of the full system have been deduced based on the analytical solution of the cantilevered piezoelectric energy harvester with Euler-Bernoulli beam assumptions. To study the advantages and disadvantages of the proposed system, the results are compared with those of an array of disconnected beams (with no springs). The analytical model is validated with experimental measurements of three bimorph beams with and without springs.

Both analytical and experimental results show that connecting the array of beams with springs allows increasing the frequency band of operation and increasing the amount of power generated. Thus, the maximal frequency band for a certain array of beams can be extended by connecting them by springs, and at the same time the output power is increased. The performance of this system can be optimized for specific applications by a proper selection of the tip masses and spring constants.

The preliminary experimental results are consistent with the analytical ones, but further research is necessary to improve the experimental setup and to investigate different experimental conditions, such as variations of the springs stiffness.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

S. R. Anton and H. A. Sodano, “A review of power harvesting using piezoelectric materials (2003–2006),” Smart Materials and Structures, vol. 16, no. 3, article R1, 2007.

S. Chalasani and J. M. Conrad, “A survey of energy harvesting sources for embedded systems,” in Proceedings of the IEEE SoutheastCon, pp. 442–447, April 2008.

K. A. Cook-Chennault, N. Thambi, and A. M. Sastry, “Powering MEMS portable devices—a review of non-regenerative and regenerative power supply systems with special emphasis on piezoelectric energy harvesting systems,” Smart Materials and Structures, vol. 17, no. 4, Article ID 043001, 2008.

A. Dewan, S. U. Ay, M. N. Karim, and H. Beyenal, “Alternative power sources for remote sensors: a review,” Journal of Power Sources, vol. 245, pp. 129–143, 2014.

A. Harb, “Energy harvesting: state-of-the-art,” Renewable Energy, vol. 36, no. 10, pp. 2641–2654, 2011.

S. Saadon and O. Sidek, “A review of vibration-based MEMS piezoelectric energy harvesters,” Energy Conversion and Management, vol. 52, no. 1, pp. 500–504, 2011.

H. A. Sodano, D. J. Inman, and G. Park, “A review of power harvesting from vibration using piezoelectric materials,” Shock and Vibration Digest, vol. 36, no. 3, pp. 197–205, 2004.

R. J. M. Vullers, R. van Schaijk, I. Doms, C. Van Hoof, and R. Mertens, “Micropower energy harvesting,” Solid-State Electronics, vol. 53, no. 7, pp. 684–693, 2009.

R. Caliò, U. B. Rongala, D. Camboni et al., “Piezoelectric energy harvesting solutions,” Sensors, vol. 14, no. 3, pp. 4755–4790, 2014.

J. Twiefel and H. Westermann, “Survey on broadband techniques for vibration energy harvesting,” Journal of Intelligent Material Systems and Structures, vol. 24, no. 11, pp. 1291–1302, 2013.

S. M. Shahruz, “Design of mechanical band-pass filters for energy scavenging,” Journal of Sound and Vibration, vol. 292, no. 3–5, pp. 987–998, 2006.

S. M. Shahruz, “Limits of performance of mechanical band-pass filters used in energy scavenging,” Journal of Sound and Vibration, vol. 293, no. 1-2, pp. 449–461, 2006.

S. M. Shahruz, “Design of mechanical band-pass filters for energy scavenging: multi-degree-of-freedom models,” Journal of Vibration and Control, vol. 14, no. 5, pp. 753–768, 2008.

J.-Q. Liu, H.-B. Fang, Z.-Y. Xu et al., “A MEMS-based piezoelectric power generator array for vibration energy harvesting,” Microelectronics Journal, vol. 39, no. 5, pp. 802–806, 2008.

H. C. Lin, P. H. Wu, I. C. Lien, and Y. C. Shu, “Analysis of an array of piezo-electric energy harvesters connected in series,” Smart Materials and Structures, vol. 22, no. 9, Article ID 094026, 2013.

S.-J. Jang, E. Rustighi, M. J. Brennan, Y. P. Lee, and H.-J. Jung, “Design of a 2DOF vibrational energy harvesting device,” Journal of Intelligent Material Systems and Structures, vol. 22, no. 5, pp. 443–448, 2011.

L. Tang and Y. Yang, “A multiple-degree-of-freedom piezoelectric energy harvesting model,” Journal of Intelligent Material Systems and Structures, vol. 23, no. 14, pp. 1631–1647, 2012.

Z. Yang and J. Yang, “Connected vibrating piezoelectric bimorph beams as a wide-band piezoelectric power harvester,” Journal of Intelligent Material Systems and Structures, vol. 20, no. 5, pp. 569–574, 2009.

W. Zhou, G. R. Penamalli, and L. Zuo, “An efficient vibration energy harvester with a multi-mode dynamic magnifier,” Smart Materials and Structures, vol. 21, no. 1, Article ID 015014, 2012.

D. S. Nguyen, E. Halvorsen, G. U. Jensen, and A. Vogl, “Fabrication and characterization of a wideband MEMS energy harvester utilizing nonlinear springs,” Journal of Micromechanics and Microengineering, vol. 20, no. 12, Article ID 125009, 2010.

S. D. Nguyen and E. Halvorsen, “Nonlinear springs for bandwidth-tolerant vibration energy harvesting,” Journal of Microelectromechanical Systems, vol. 20, no. 6, pp. 1225–1227, 2011.

S. Boisseau, G. Despesse, and B. A. Seddik, “Nonlinear h-shaped springs to improve efficiency of vibration energy harvesters,” Journal of Applied Mechanics, vol. 80, no. 6, Article ID 061013, 2013.

S. Priya, “Advances in energy harvesting using low profile piezoelectric transducers,” Journal of Electroceramics, vol. 19, no. 1, pp. 165–182, 2007.

S. Roundy, P. K. Wright, and J. Rabaey, “A study of low level vibrations as a power source for wireless sensor nodes,” Computer Communications, vol. 26, no. 11, pp. 1131–1144, 2003.

S. Roundy and P. K. Wright, “A piezoelectric vibration based generator for wireless electronics,” Smart Materials and Structures, vol. 13, no. 5, pp. 1131–1142, 2004.

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, 2004.

N. E. Dutoit, B. L. Wardle, and S.-G. Kim, “Design considerations for mems-scale piezoelectric mechanical vibration energy harvesters,” Integrated Ferroelectrics, vol. 71, no. 1, pp. 121–160, 2005.

S.-N. Chen, G.-J. Wang, and M.-C. Chien, “Analytical modeling of piezoelectric vibration-induced micro power generator,” Mechatronics, vol. 16, no. 7, pp. 379–387, 2006.

A. Erturk and D. J. Inman, “Issues in mathematical modeling of piezoelectric energy harvesters,” Smart Materials and Structures, vol. 17, no. 6, Article ID 065016, 2008.

F. Lu, H. P. Lee, and S. P. Lim, “Modeling and analysis of micro piezoelectric power generators for micro-electromechanical-systems applications,” Smart Materials and Structures, vol. 13, no. 1, pp. 57–63, 2004.

S. Jiang, X. Li, S. Guo, Y. Hu, J. Yang, and Q. Jiang, “Performance of a piezoelectric bimorph for scavenging vibration energy,” Smart Materials and Structures, vol. 14, no. 4, pp. 769–774, 2005.

A. Erturk and D. J. Inman, “A distributed parameter electromechanical model for cantilevered piezoelectric energy harvesters,” Journal of Vibration and Acoustics, vol. 130, no. 4, Article ID 041002, 2008.

A. Erturk and D. J. Inman, “An experimentally validated bimorph cantilever model for piezoelectric energy harvesting from base excitations,” Smart Materials and Structures, vol. 18, no. 2, Article ID 025009, 2009.