Efficiency of energy conversion for a piezoelectric power harvesting system
Y C Shu and I C Lien
Institute of Applied Mechanics, National Taiwan University, Taipei 106, Taiwan, Republic of China
E-mail: yichung@spring.iam.ntu.edu.tw
Received 22 June 2006, in final form 7 September 2006
Published 29 September 2006
Online at stacks.iop.org/JMM/16/2429
Abstract
This paper studies the energy conversion efficiency for a rectified piezoelectric power harvester. An analytical model is proposed, and an expression of efficiency is derived under steady-state operation. In addition, the relationship among the conversion efficiency, electrically induced damping and ac–dc power output is established explicitly. It is shown that the optimization criteria are different depending on the relative strength of the coupling. For the weak electromechanical coupling system, the optimal power transfer is attained when the efficiency and induced damping achieve their maximum values. This result is consistent with that observed in the recent literature. However, a new finding shows that they are not simultaneously maximized in the strongly coupled electromechanical system.
1. Introduction
Advances in low-power electronic design and fabrication have opened the possibility of self-powered microsensors and communication nodes [1]. At the same time, the need to power remote systems or embedded devices independently has motivated many research efforts harvesting electric energy from various ambient sources, including solar power, thermal gradients and vibration [2]. Among these energy scavenging sources, mechanical vibration is a potential power source that is abundant enough to be of use, is easily accessible through microelectromechanical systems (MEMS) technology for conversion to electric energy, and is ubiquitous in applications ranging from small household appliances to large infrastructures [3, 4].
Piezoelectric vibration-to-electricity converters have received much attention as transducers, since they have high electromechanical coupling, require no external voltage source and are particularly attractive for use in MEMS [5–8]. As a result, piezoelectric materials for scavenging energy from ambient vibration sources have recently seen a dramatic rise in use for power harvesting. This includes the use of resonant piezoelectric-based structures of cantilever beam configuration [9–15] as well as plate (membrane) configuration [16–19]. Other harvesting schemes include the use of long strips of piezoelectric polymers in ocean or river-water flows [20, 21], the use of piezoelectric ‘cymbal’ transducers [22, 23], and the use of piezoelectric windmill for generating electric power from wind energy [24].
Jeon et al [5] have successfully developed the first MEMSbased micro-scale power generator using a {3-3} mode of PZT transducer. A 170 µm × 260 µm PZT beam has been fabricated, and a maximum dc voltage of 3 V across the load 10.1 M has been observed. In addition, the energy density of the power generator has been estimated at around 0.74 mWh cm−2, which compares favourably to the use of lithium ion batteries. Roundy et al [3] subsequently created prototyes of thin PZT structures with a target volume power density of . Recently, duToit et al [25] provided in-depth design principles for MEMS-scale piezoelectric energy harvesters and proposed a prototype of 30 µW cm−3 from low-level vibration. Related works on the modelling of miniaturized piezoelectric power harvesting devices can be found in [26–28].
As the testing, characterization and fabrication of MEMSscale energy harvesters are not always available compared to the similar tasks in bulk power harvesters, a normalization scheme is particularly useful for comparing the performance of micro-scale power harvesters. One good method uses the parameter of efficiency of mechanical to electric energy conversion. Umeda et al [29, 30] have studied the efficiency of mechanical impact energy to electric energy using a piezoelectric vibrator. Goldfarb and Jones [31] subsequently investigated the efficiency of the piezoelectric material in a stack configuration for converting mechanical harmonic excitation into electric energy. Roundy [32] provided an expression for effectiveness that can be used to compare various approaches and designs for vibration-based energy harvesting devices. Recently, in contrast to efforts where the conversion efficiency was examined numerically [29], Richards et al [33] derived an analytic formula to predict the energy conversion efficiency of piezoelectric energy harvesters in the case of ac power output. Since the electronic load requires a stabilized dc voltage while a vibrating piezoelectric element generates an ac voltage, the desired output needs to be rectified, filtered and regulated to ensure the electric compatibility. Thus, we investigate the conversion efficiency for a rectified piezoelectric power harvesting system based on the analysis proposed by Shu and Lien [34] in section 2. We show that the conversion efficiency is dependent on the frequency ratio, the normalized resistance and, in particular, the ratio of the electromechanical coupling coefficient to the mechanical damping ratio. In general, the conversion efficiency can be improved with a larger coupling coefficient and smaller damping. Recently, Cho et al [35, 36] performed a series of experiments and proposed a set of design guidelines for the performance optimization of micromachined piezoelectric membrane generators by enhancing the electromechanical coupling coefficient.

flowchart
graph TD A["F(t)"] --> B["M"] B --> C["u(t)"] C --> D["K"] D --> E["ηm"] E --> F["Piezoelectric element characterized by Θ&Cp"] G["Energy storage system"] --> H["Regulator"] H --> I["Ce"] I --> J["Ground"]
Figure 1. An equivalent model for a piezoelectric vibration energy harvesting system.
When an energy harvester is applied to a system, energy is removed from the vibrating structure and supplied to the desired electronic components, resulting in additional damping of the structure [37, 38]. Because the efficiency is defined as the ratio of the time-averaged power dissipated across the load to that done by the external force, electrically induced damping can be defined explicitly, and its connection to the conversion efficiency is established in section 3. It is shown that the load to maximize the conversion efficiency is the same as that to maximize the induced damping. However, the extraction of harvested power may not be simultaneously optimized. It is demonstrated in section 4 that optimization criteria vary according to the relative strength of the coupling. The conclusions are drawn in section 5.
Finally, our result can be applied to the investigation of charging a battery from a vibrating piezoelectric harvester. It is shown that direct charging will result in a low efficiency of energy conversion, since the equivalent impedance of the battery may not match that of the optimal load in most situations. Ottman et al [39, 40] have developed an adaptive electric circuit to optimize the energy transfer from the piezoelectric element to the stored device. It is based on the principle of load impedance adaptation by tuning the load impedance to obtain a higher power flow. Related work based on the synchronous electric charge extraction can be found in [41–44].
2. A piezoelectric power harvesting model
Consider an energy conversion device which includes a vibrating piezoelectric structure together with an energy storage system. If the modal density of such a device is widely separated and the structure is vibrating at around its resonance frequency, we may model the power generator as a mass+spring+damper+piezo structure, as schematically shown in figure 1 [33, 39, 44]. It consists of a piezoelectric element coupled to a mechanical structure. In this approach, a forcing function is applied to the system and an effective mass M is bounded on a spring of effective stiffness on a damper of coefficient , and on a piezoelectric element characterized by effective piezoelectric coefficient and capacitance These effective coefficients are dependent on the material constants and the design of energy harvesters and can be derived using the standard modal analysis [25, 38, 45, 46].
Let u be the displacement of the mass and be the voltage across the piezoelectric element. The governing equations of the piezoelectric vibrator can be described by [25, 38, 43, 44]
where I (t) is the current flowing into the specified circuit. Since most applications of piezoelectric materials for power generation involve the use of periodic straining of piezoelectric elements, the vibrating generator is assumed to be driven at around resonance by the harmonic excitation
where is the constant magnitude and w (in radians per second) is the angular frequency of vibration.
The power generator considered here is connected to a storage circuit system, as illustrated in figure 1. Since the electrochemical battery needs a stabilized dc voltage while a vibrating piezoelectric element generates an ac voltage, this requires a suitable circuit to ensure the electric compatibility. Typically an ac–dc rectifier followed by a filtering capacitance is added to smooth the dc voltage, as shown in figure 1. A controller placed between the rectifier output and the battery is included to regulate the output voltage. Figure 2 is a simplified energy harvesting circuit commonly adopted for design analysis. It can be used to estimate an upper bound of the real power that the piezoelectric generator is able to deliver at a given excitation. Note that the regulation circuit and battery are replaced with an equivalent resistor R and is the rectified voltage across it.

text_image
Piezoelectric vibrator I(t) Ce R Vc(t)
Figure 2. A typical ac–dc harvesting circuit.

line chart
| Time (t) | U(t) | Vp(t) |
|---|---|---|
| t0 | 0 | -Vc |
| t1 | -u0 | -Vc |
| t2 | 0 | -Vc |
| t3 | 0 | -Vc |
| t4 | -u0 | -Vc |
| t5 | 0 | -Vc |
Figure 3. Typical waveforms of displacement u(t) and piezo voltage for an ac–dc power harvesting circuit.
The common approach to having the stable output dc voltage is to assume that the filter capacitor e is large enough so that the rectified voltage Vc is essentially constant [39]. Specifically, , where Vc(t)- and are the average and ripple of , respectively. This average is independent of provided that the time constant is much larger than the oscillating period of the generator [43]. The magnitude of , however, depends on and is negligible for large . Under this hypothesis, , and therefore in the following, we use instead of -, to represent the average of for notation simplicity.
The rectifying bridge shown in figure 2 is assumed to be perfect here. Thus, it is open circuited if the piezo voltage is smaller than the rectified voltage . As a result, the current flowing into the circuit vanishes, and this implies varies proportionally with respect to , as seen from . On the other hand, when reaches , the bridge conducts and the piezo voltage is kept equal to the rectified voltage, , Finally, the conduction in the rectifier diodes is blocked again when the absolute value of the piezo voltage starts decreasing. Typical waveforms of u(t) and satisfying these properties are schematically shown in figure 3.
To solve connected to an ac–dc circuit shown in figure 2, we first determine the relation between the average value of the rectified voltage and displacement magnitude. From figure 3, the steady-state solution of u(t) is assumed to take the following form:
with being the constant magnitude. Let be the period of vibration, and and be two time instants such that the displacement u undergoes from the minimum −u0 to the maximum as illustrated in figure . Assume that during the semi-period from ti to . It follows that . Note that voltage and for during which the rectifier conducts. This gives
since the average current flowing through the capacitance is zero, i.e., at the steady-state operation. The integration of (2) from time to is therefore
or
Note that (5) is identical with that derived by [39, 43, 44].
We next need to find out to determine . There are two approaches to estimating it in the case of ac–dc power harvesting system in the recent literature [39, 43, 44]. The first one models the piezoelectric device as the current source in parallel with its internal electrode capacitance , 11, 39]. It is based on the assumption that the internal current source of the generator is independent of the external load impedance. This is equivalent to assuming that the coupling is very weak and the term can be dropped from (1). On the other hand, if the coupling is not so weak, Guyomar et al [43] and Lefeuvre et al [44] have provided another approach to estimate the displacement magnitude by assuming that the external forcing function and the velocity of the mass are in phase. Recently, Shu and Lien [34] have proposed a new method for determining u0 without the uncoupled and in-phase assumptions. They have shown that this new estimation is more accurate than the other two. We here briefly outline the steps of derivation of since some of them are required to derive the efficiency of energy conversion.
Consider the balance of energy. Let (1) be multiplied by u(t ) ˙ and (2) be multiplied by . Integration of the addition of these two equations from time to gives the equation of the energy balance
where
Note that (3) and (4) are used to derive (7).
Right now there are two equations (5) and (6) and three unknowns and θ. We need the third one to solve them. Differentiating (1) with respect to time t and using (2), we have
Integrating (8) with respect to time t from to and using (4) provides the third equation
Thus, the unknown variable θ can be eliminated from (7) and (9). This gives
As the magnitude of displacement is related with the rectified voltage by (5), the above equation (10) can be further simplified to find The result is
The average harvested power can also be obtained once is determined since
To summarize, the normalized displacement rectified voltage and average harvested power can be expressed by
where several non-dimensionless variables are introduced by
Above is the alternative electromechanical coupling coefficient1, the mechanical damping ratio, the natural frequency of the short circuit, and r the normalized frequency and electric resistance. Note that there are two resonances for the system since the piezoelectric structure exhibits both short circuit and open circuit stiffness. They are defined by
where and are the frequency ratios of the short circuit and the open circuit, respectively. Note that the shift in device natural frequency is pronounced if the coupling factor is large. Besides, (12), (13) and (14) are validated both numerically and experimentally, and are also compared well with other existing estimates in [34].
Since (14) is expressed in terms of non-dimensionless parameters, it can be used as the power normalization scheme to compare the performance and efficiency of the devices relatively. This is particularly useful in the design of a MEMSscale power generator since the testing, characterization and fabrication of micro-scale energy harvesters are not always available compared to the similar tasks in bulk power harvesters. Further, in most vibration-based power harvesting systems the vibration source is due to the periodic excitation of some base. This gives , where A is the magnitude of acceleration of the exciting base. Therefore, the harvested average power per unit mass is described by
This shows that the harvested average power per unit mass depends on the input vibration characteristics (frequency ratio and acceleration , the normalized electric resistance the short circuit resonance , the mechanical damping ratio and the overall electromechanical coupling coefficient of the system. Thus, the scheme to optimize the power either by tuning the electric resistance, selecting suitable operation points, or adjusting the coupling coefficient by optimal structural design can be guided completely by (14).
3. Conversion efficiency and electrically induced damping
The efficiency of mechanical to electrical energy conversion is a fundamental parameter in order to compare energy harvesters of various sizes and with different vibration inputs. If the generator is excited by a harmonic forcing function as in (3), the energy conversion efficiency is commonly defined as the time-averaged power ratio by
where · dt denotes the average over time [26, 33]. Above is the time-averaged power dissipated across the load resistor R and is the time-averaged power done by the external force. The balance of energy in (6) gives
where is the time-averaged power dissipated due to the structural damping. Further, the use of (7) provides an expression of efficiency of energy conversion under the steadystate operation. Indeed, it is given by
due to (5). It can also be written in terms of the nondimensionless parameters defined in (15) by
It is clear from (19) that the conversion efficiency depends on the normalized resistance r, the applied frequency ratio and, in particular, the relative magnitudes of the electromechanical coupling coefficient and the mechanical damping ratio . However, one has to be cautious when applying (19) to the study of conversion efficiency far below resonance since the dielectric loss term is not included in the present model. We refer to [29] which has studied the effect of dielectric loss on the efficiency of mechanical impact energy transformed into electric energy.
Williams and Yates [48] have proposed a model to study the conversion of the kinetic energy to electric power without specifying the mechanism by which the conversion process takes place. It is based on the idea that the conversion of energy from the oscillating mass to electricity is similar to a linear damper in the conventional mass-spring system. According to their model, the total damping ratio of the system can be decomposed as
where is the electrically induced damping ratio due to the removal of mechanical energy from the vibrating system. Using the efficiency derived in (19), we can determine the induced damping ζe for an ac–dc piezoelectric power harvesting system. Indeed, the efficiency of the energy conversion can be re-defined by
if the effect of the electric system on the mechanical system is proportional to the velocity of the oscillating mass. Thus, from (19) and (20), the induced damping added to the system can be found to be
Note that is small at a small load resistance since only a slight fraction of energy is removed from the system. It is also small at a larger electric load since the circuit behaves like an open-circuit condition, preventing the generated charges from flowing out of the piezoelectric elements.
Another quantity known as the loss factor is commonly used for comparing the damping capacity of a vibrating system. It is defined as the ratio of the energy dissipated per radian and the total strain energy, and is related to the total damping ratio by
for small values of damping [49]. Thus, the loss factor added to the system due to the energy dissipated across the load resistor is therefore
We now turn to the study of the optimal efficiency of energy conversion which is important in the design of an energy harvester. From (19), the normalized load to maximize efficiency for fixed , can be obtained according to
This gives
Besides, it can be shown from (21) that the electrically induced damping ratio evaluated at has also achieved the maximum value; i.e.,
This in turn gives the maximum value of the electrically induced loss factor
Lesieutre et al [37] have studied the induced electric damping associated with a piezoelectric energy harvesting system and derived an expression for the maximum loss factor. It is
for small values of which is defined by
Obviously, (25) and (26) are almost identical provided that the coupling coefficient is small and the applied frequency ratio . Note that Lesieutre et al [37] have derived an expression of the maximum loss factor by assuming that it occurs at the condition of the optimal power transfer. Their argument is in generally true for most common cases except the situation where the ratio is large. We will discuss it in the next section.
4. Discussion
The shift in device natural frequency from to is pronounced if the electromechanical coupling coefficient is large, as seen from (16). It occurs often in devices whose the piezoelectric element’s contribution to the overall structural stiffness is significant [25]. In addition, this frequency shift could be large at the piezoelectric micro power generator operated utilizing the {3-3} mode since the piezoelectric effect is further enhanced in the longitudinal mode [5, 25]. The consequence of pronounced shift in natural frequency from to leads to completely distinct optimization schemes for maximum power extraction. Shu and Lien [34] have studied the ac–dc power output for a rectified piezoelectricbased harvesting device. They have shown that the average harvested power has two identical peaks evaluated at two different electric loads at the respective operating points provided that k e while it has only one peak otherwise. ζm Thus, the study of relation among the conversion efficiency, electrically induced damping and power transfer has to be classified according to the relative magnitudes of the coupling coefficient and the mechanical damping ratio.
4.1. Weak electromechanical coupling
In the following, we take and which are typical parameters for piezoelectric power generators designed as the cantilever-beam type and operated in the {3-1} mode.
The conversion efficiency, normalized displacement and harvested power against the frequency ratio are plotted in figure for various normalized resistances. In addition, they are also plotted versus normalized resistance, with varying frequency ratios in figures . Consider figures and (d) first. The conversion efficiency is small and around 18% for very small load at , and is increasing as the resistance becomes large. According to (23), it achieves its maximum value
at r eff 2 decreases as the load exceeds for ≈ 1. The conversion efficiency then . Besides, the overlapping of curves in figure shows that the efficiency of energy conversion is not sensitive with respect to the frequency ratio in the case of the weak electromechanical coupling.
We next consider figures 4(b) and . Each curve of displacement has a peak around the resonance whose value depends on the resistance. From (21) the electrically induced damping is small at both small and large electric loads while it attains its maximum value at , which in turn also gives the maximum conversion efficiency of the generator. a result, the peak of displacement descends until to the point which corresponds to the maximum induced damping as well as efficiency. When r moves past , the induced damping decreases, resulting in the elevation of the peak of displacement. Figure illustrates the dependence of displacement on resistance at the different operating points. The displacement becomes monotonically decreasing if the device is excited in the vicinity of its short circuit resonance while it becomes monotonically increasing at the open circuit resonance. If the frequency is operated around the middle point of and , the curve has a local minimum at the turning point close to where the induced damping is maximized.
Note that our theoretical prediction on the behaviour of displacement against frequency ratio for various electric loads, as illustrated in figure , qualitatively agrees well with the experimental results observed by Lesieutre et al (see figure 6 in [37]).
Each curve of harvested power has a maximum around the resonance whose value depends on the resistance, as demonstrated in figure . In contrast with displacement, the peak of power and the electrically induced damping simultaneously ascend as the electric load increases. It attains the maximum value at corresponding to the local minimum of displacement as well as the maximum of efficiency. After that, the peak of power descends when the load exceeds . Finally, the optimal load to maximize the harvested power is illustrated in figure with varying frequency ratios. The peaks of power are slightly higher for frequencies operated around than those operated around . However, unlike the conversion efficiency in figure which is insensitive to the frequency ratio, a small derivation from the optimal load results in a more significant drop in power operated around than that operated around .
4.2. Strong electromechanical coupling
We use and here to demonstrate the effect of strong coupling on the relation among the efficiency, damping and the power transfer.
The conversion efficiency, normalized displacement and power versus frequency ratio are plotted, with varying load resistances in figures . Besides, they are also plotted against normalized resistance in figures for a variety of frequency ratios. Unlike the case of weak coupling, the average harvested power has two identical optimal peaks evaluated at different resistances and frequencies; i.e.,
where
Obviously, ≈ ≈ as illustrated 1 as demonstrated in figure . In contrast with figure , the peak of harvested power decreases while the electrically induced damping increases as the load exceeds . The envelope of these peaks

line chart
| Frequency Ratio | r:0.1 | r:0.25 | r:0.5 | r:π/2 | r:5 | r:10 | r:20 |
|---|---|---|---|---|---|---|---|
| 0.9 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 0.95 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.0 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.05 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.1 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.15 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.2 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.25 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.3 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.35 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.4 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.45 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.5 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.55 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.6 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.65 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.7 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.75 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.8 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.85 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.9 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 1.95 | 0.3 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.2 |
| 2.0 | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| - | - | - | - | - | - | - | - |
| … | … … | … … | … … | … … | … | … | … |

line chart
| Normalized Resistance | Ω:1.016 | Ω:1.024 | Ω:1.032 | Ω:1.039 | Ω:1.047 | Ω:1.055 | Ω:1.062 |
|---|---|---|---|---|---|---|---|
| 0 | 0.45 | 0.45 | 0.45 | 0.45 | 0.45 | 0.45 | 0.45 |
| 4 | 0.38 | 0.38 | 0.38 | 0.38 | 0.38 | 0.38 | 0.38 |
| 8 | 0.32 | 0.32 | 0.32 | 0.32 | 0.32 | 0.32 | 0.32 |
| 12 | 0.26 | 0.26 | 0.26 | 0.26 | 0.26 | 0.26 | 0.26 |
| 16 | 0.21 | 0.21 | 0.21 | 0.21 | 0.21 | 0.21 | 0.21 |
| 20 | 0.18 | 0.18 | 0.18 | 0.18 | 0.18 | 0.18 | 0.18 |

line chart
| Frequency Ratio | r:0.1 | r:0.25 | r:0.5 | r:π/2 | r:5 | r:10 | r:20 |
|---|---|---|---|---|---|---|---|
| 0.9 | ~3.5 | ~3.8 | ~4.0 | ~3.6 | ~3.7 | ~3.4 | ~3.2 |
| 0.95 | ~6.0 | ~6.5 | ~7.0 | ~6.2 | ~6.8 | ~6.0 | ~5.5 |
| 1.0 | ~10.0 | ~11.0 | ~12.0 | ~10.5 | ~11.5 | ~10.0 | ~9.5 |
| 1.05 | ~13.0 | ~12.5 | ~11.5 | ~10.0 | ~9.5 | ~9.0 | ~8.5 |
| 1.1 | ~10.0 | ~9.5 | ~8.5 | ~7.5 | ~6.5 | ~6.0 | ~5.5 |
| 1.15 | ~6.0 | ~5.5 | ~4.5 | ~3.8 | ~3.2 | ~3.0 | ~2.8 |
| 1.2 | ~3.5 | ~3.2 | ~2.8 | ~2.5 | ~2.2 | ~2.0 | ~1.8 |

line chart
| Normalized Resistance | Ω:1.016 | Ω:1.024 | Ω:1.032 | Ω:1.039 | Ω:1.047 | Ω:1.055 | Ω:1.062 |
|---|---|---|---|---|---|---|---|
| 0 | 14.0 | 13.5 | 13.0 | 12.5 | 12.0 | 11.5 | 11.0 |
| 4 | 8.5 | 9.0 | 9.5 | 9.8 | 9.7 | 9.6 | 9.4 |
| 8 | 7.8 | 9.5 | 10.0 | 10.2 | 10.1 | 10.0 | 9.8 |
| 12 | 7.5 | 10.0 | 10.5 | 10.7 | 10.6 | 10.5 | 10.2 |
| 16 | 7.3 | 10.5 | 11.0 | 11.2 | 11.1 | 11.0 | 10.7 |
| 20 | 7.2 | 11.0 | 11.5 | 11.7 | 11.6 | 11.5 | 11.2 |

line chart
| Frequency Ratio | r:0.1 | r:0.25 | r:0.5 | r:π/2 | r:5 | r:10 | r:20 |
|---|---|---|---|---|---|---|---|
| 0.9 | 0.25 | 0.25 | 0.25 | 0.25 | 0.25 | 0.25 | 0.25 |
| 0.95 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 |
| 1.0 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 |
| 1.05 | 1.75 | 1.75 | 1.75 | 1.75 | 1.75 | 1.75 | 1.75 |
| 1.1 | 1.25 | 1.25 | 1.25 | 1.25 | 1.25 | 1.25 | 1.25 |
| 1.15 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 |
| 1.2 | 0.25 | 0.25 | 0.25 | 0.25 | 0.25 | 0.25 | 0.25 |
(c)

line chart
| Normalized Resistance | Ω:1.016 | Ω:1.024 | Ω:1.032 | Ω:1.039 | Ω:1.047 | Ω:1.055 | Ω:1.062 |
|---|---|---|---|---|---|---|---|
| 0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 |
| 2 | 1.8 | 1.9 | 2.0 | 2.1 | 2.2 | 2.3 | 2.4 |
| 4 | 1.5 | 1.7 | 1.8 | 1.9 | 2.0 | 2.1 | 2.2 |
| 6 | 1.2 | 1.4 | 1.5 | 1.6 | 1.7 | 1.8 | 1.9 |
| 8 | 0.9 | 1.1 | 1.2 | 1.3 | 1.4 | 1.5 | 1.6 |
| 10 | 0.7 | 0.9 | 1.0 | 1.1 | 1.2 | 1.3 | 1.4 |
| 12 | 0.5 | 0.7 | 0.8 | 0.9 | 1.0 | 1.1 | 1.2 |
| 14 | 0.4 | 0.6 | 0.7 | 0.8 | 0.9 | 1.0 | 1.1 |
| 16 | 0.3 | 0.5 | 0.6 | 0.7 | 0.8 | 0.9 | 1.0 |
| 18 | 0.25 | 0.45 | 0.55 | 0.65 | 0.75 | 0.85 | 0.95 |
| 20 | 0.2 | 0.4 | 0.5 | 0.6 | 0.7 | 0.8 | 0.9 |
Figure 4. Weak electromechanical coupling demonstrated using and (a), (b) and (c) are the conversion efficiency, displacement and power against frequency ratio for various resistances, while (d), (e) and (f ) are those against resistance with varying frequency ratios.
has a local minimum which is closely related to the maximum induced damping When the load further increases, the peak of power then rises to the second maximum, closely corresponding to as described by (29). Figure 5(f ) shows the dependence of the harvested power versus the normalized resistance for a variety of frequency ratios. Switching between these two peaks, corresponding to t, is attained by varying the applied frequency from to . It indicates the importance of selecting the correct operating point in the case of strong coupling.
The dependence of conversion efficiency on the frequency as well as on the load resistance in the strong coupling case is qualitatively similar to the case of weak coupling. But the magnitude of conversion efficiency is much higher here

line chart
| Frequency Ratio | r:0.067 | r:0.15 | r:0.3 | r:1.058 | r:4 | r:8 | r:18.27 |
|---|---|---|---|---|---|---|---|
| 0.8 | 0.9 | 0.85 | 0.8 | 0.75 | 0.7 | 0.65 | 0.5 |
| 1.0 | 0.85 | 0.8 | 0.75 | 0.7 | 0.65 | 0.6 | 0.45 |
| 1.2 | 0.8 | 0.75 | 0.7 | 0.65 | 0.6 | 0.55 | 0.4 |
| 1.4 | 0.75 | 0.7 | 0.65 | 0.6 | 0.55 | 0.5 | 0.35 |
| 1.6 | 0.7 | 0.65 | 0.6 | 0.55 | 0.5 | 0.45 | 0.3 |
(a)

line chart
| Normalized Resistance | Ω:1.021 | Ω:1.140 | Ω:1.212 | Ω:1.265 | Ω:1.304 | Ω:1.342 | Ω:1.394 |
|---|---|---|---|---|---|---|---|
| 0 | 0.85 | 0.85 | 0.85 | 0.85 | 0.85 | 0.85 | 0.85 |
| 5 | 0.78 | 0.77 | 0.76 | 0.75 | 0.74 | 0.73 | 0.72 |
| 10 | 0.72 | 0.70 | 0.69 | 0.68 | 0.67 | 0.66 | 0.65 |
| 15 | 0.68 | 0.66 | 0.65 | 0.64 | 0.63 | 0.62 | 0.61 |
| 20 | 0.64 | 0.62 | 0.61 | 0.60 | 0.59 | 0.58 | 0.57 |
| 25 | 0.60 | 0.58 | 0.57 | 0.56 | 0.55 | 0.54 | 0.53 |
(d)

line chart
| Frequency Ratio | r:0.067 | r:0.15 | r:0.3 | r:1.058 | r:4 | r:8 | r:18.27 |
|---|---|---|---|---|---|---|---|
| 0.8 | ~2.0 | ~1.5 | ~1.0 | ~0.8 | ~0.7 | ~0.6 | ~0.5 |
| 1.0 | ~9.5 | ~6.5 | ~4.5 | ~3.0 | ~2.5 | ~2.0 | ~1.5 |
| 1.2 | ~3.0 | ~2.5 | ~2.0 | ~1.5 | ~1.0 | ~0.8 | ~0.6 |
| 1.4 | ~1.0 | ~0.8 | ~0.6 | ~0.5 | ~0.4 | ~0.3 | ~0.2 |
| 1.6 | ~0.5 | ~0.4 | ~0.3 | ~0.2 | ~0.1 | ~0.1 | ~0.1 |
(b)

line chart
| Normalized Resistance | Q:1.021 | Q:1.140 | Q:1.212 | Q:1.265 | Q:1.304 | Q:1.342 | Q:1.394 |
|---|---|---|---|---|---|---|---|
| 0 | 10.0 | 3.8 | 2.5 | 2.0 | 1.5 | 1.2 | 1.0 |
| 5 | 2.5 | 3.5 | 3.0 | 2.8 | 2.2 | 2.0 | 1.8 |
| 10 | 3.0 | 4.0 | 3.5 | 3.2 | 2.8 | 2.5 | 2.3 |
| 15 | 3.5 | 4.5 | 4.0 | 3.8 | 3.2 | 3.0 | 2.8 |
| 20 | 4.0 | 5.0 | 4.5 | 4.2 | 3.8 | 3.5 | 3.2 |
| 25 | 4.5 | 5.5 | 5.0 | 4.8 | 4.2 | 4.0 | 3.8 |
(e)

line chart
| Frequency Ratio | r:0.067 | r:0.15 | r:0.3 | r:0.58 | r:4 | r:8 | r:18.27 |
|---|---|---|---|---|---|---|---|
| 0.8 | ~0.1 | ~0.1 | ~0.1 | ~0.1 | ~0.1 | ~0.1 | ~0.1 |
| 1.0 | ~2.5 | ~2.2 | ~1.8 | ~1.5 | ~1.0 | ~0.5 | ~0.3 |
| 1.2 | ~0.5 | ~1.5 | ~1.2 | ~1.0 | ~0.5 | ~0.3 | ~0.2 |
| 1.4 | ~0.2 | ~2.5 | ~2.2 | ~2.0 | ~2.5 | ~2.0 | ~1.8 |
| 1.6 | ~0.1 | ~0.2 | ~0.2 | ~0.2 | ~0.2 | ~0.2 | ~0.2 |
(c)

line chart
| Normalized Resistance | Ω:1.021 | Ω:1.140 | Ω:1.212 | Ω:1.265 | Ω:1.304 | Ω:1.342 | Ω:1.394 |
|---|---|---|---|---|---|---|---|
| 0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 |
| 5 | 1.5 | 1.7 | 1.8 | 1.9 | 2.0 | 2.1 | 2.2 |
| 10 | 1.0 | 1.2 | 1.4 | 1.6 | 1.8 | 2.0 | 2.3 |
| 15 | 0.5 | 0.7 | 0.9 | 1.1 | 1.3 | 1.5 | 2.4 |
| 20 | 0.2 | 0.4 | 0.6 | 0.8 | 1.0 | 1.2 | 2.5 |
| 25 | 0.1 | 0.2 | 0.4 | 0.5 | 0.7 | 0.9 | 2.4 |
( f )
Figure 5. Strong electromechanical coupling demonstrated using and , (b) and (c) are the conversion efficiency, displacement and power against frequency ratio for various resistances, while (d), (e) and (f ) are those against resistance with varying frequency ratios.
because of the relatively large ratio k 2e . Our current example points out ≈ 90% which is much larger than that (46% in (27)) in the weak coupling case. Besides, the condition to maximize the conversion efficiency does not lead to the maximum power transfer in the strong electromechanical coupling; instead, it is closely related to the local minimum of the power envelope. The conversion efficiency takes on the smaller value around 50% at the optimal power generation, as shown in figures and (f ).
Finally, we consider figures 5(b) and (e). The effect of the induced electric damping is much pronounced in the strong coupling case than that in the case of weak coupling since, from (21), is proportional to . As a result, at the maximum induced damping, the envelope of the peaks of displacement has a local minimum whose value is much smaller than that of the first peak. Figure 5(b) demonstrates that the peak of displacement decreases approximately from 10 to 2 at the maximum induced damping, while the weak coupling case shows a moderate descent in figure 4(b). Unlike the power, the second maximum of these peaks due to the reduction of at large r is smaller than the first one since the system has more damping at the open-circuit condition. Finally, in contrast with the results obtained in the weak coupling case, the two pronounced peaks of displacement evaluated around and correspond to the maximum extraction of power rather than the local minimum point of the envelope of displacement by comparing figures 5(b) and (c) with figures 4(b) and (c).
5. Conclusion
This paper establishes the relation among the energy conversion efficiency, electrically induced damping and power transfer for a rectified piezoelectric power harvester. An analytical model is proposed, and an exact formula for the conversion efficiency is derived under the steady-state operation. It is shown that the efficiency depends on the normalized resistance the frequency ratio and, in particular, the relative magnitudes of the electromechanical coupling coefficient and the mechanical damping ratio . In general, high energy conversion efficiency can be achieved with large k2eζ ; the improvement of the coupling coefficient m for micromachined piezoelectric membrane harvesters has been recently investigated by Cho et al [35, 36].
The induced damping added to the system due to the removal of mechanical energy from the vibrating structure is obtained based on the derived formula of conversion efficiency. It is shown that the maximum conversion efficiency corresponds to the maximum induced electric damping as well as the optimal power transfer in the case of weak electromechanical coupling. This result is consistent with that observed by Lesieutre et al [37].
However, unlike [37], a new finding shows that the optimal electric load maximizing the conversion efficiency and induced electric damping is very different from that maximizing the harvested power in strongly coupled electromechanical systems. This gives completely distinct optimization schemes, since the harvested power has two identical peaks evaluated at two different electric loads at the respective operating points, provided that k 2e , while it has only one peak otherwise.
Acknowledgments
We are glad to acknowledge the Ministry of Economic Affairs for support under grant no 94-EC-17-A-05-S1-017 (WHAM-BioS).
References
[1] Chandrakasan A, Amirtharajah R, Goodman J and Rabiner W 1998 Trends in low power digital signal processing Int. Symp. Circuits Syst. 4 604–7
[2] Roundy S, Steingart D, Frechette L, Wright P and Rabaey J 2004 Power sources for wireless sensor networks Lecture Notes Comput. Sci. 2920 1–17
[3] Roundy S, Leland E S, Baker J, Carleton E, Reilly E, Lai E, Otis B, Rabaey J M, Wright P K and Sundararajan V 2005 Improving power output for vibration-based energy scavengers IEEE Pervasive Comput. 4 28–36
[4] Sodano H A, Inman D J and Park G 2004 A review of power harvesting from vibration using piezoelectric materials Shock Vib. Dig. 36 197–205
[5] Jeon Y B, Sood R, Jeong J H and Kim S G 2005 MEMS power generator with transverse mode thin film PZT Sensors Actuators A 122 16–22
[6] Poulin G, Sarraute E and Costa F 2004 Generation of electric energy for portable devices: comparative study of an electromagnetic and a piezoelectric system Sensors Actuators A 116 461–71
[7] Roundy S, Wright P K and Rabaey J 2003 A study of low level vibrations as a power source for wireless sensor nodes Comput. Commun. 26 1131–44
[8] Sodano H A, Inman D J and Park G 2005 Generation and storage of electricity from power harvesting devices J. Intell. Mater. Syst. Struct. 16 67–75
[9] Elvin N G, Elvin A A and Spector M 2001 A self-powered mechanical strain energy sensor Smart Mater. Struct. 10 293–9
[10] Elvin N G, Lajnef N and Elvin A A 2006 Feasibility of structural monitoring with vibration powered sensors Smart Mater. Struct. 15 977–86
[11] Ng T G and Liao W H 2005 Sensitivity analysis and energy harvesting for a self-powered piezoelectric sensor J. Intell. Mater. Syst. Struct. 16 785–97
[12] Roundy S and Wright P K 2004 A piezoelectric vibration based generator for wireless electronics Smart Mater. Struct. 13 1131–42
[13] Sodano H A, Inman D J and Park G 2005 Comparison of piezoelectric energy harvesting devices for recharging batteries J. Intell. Mater. Syst. Struct. 16 799–807
[14] Sodano H A, Lloyd J and Inman D J 2006 An experimental comparison between several active composite actuators for power generation Smart Mater. Struct. 15 1211–6
[15] White N M, Glynne-Jones P and Beeby S P 2001 A novel thick-film piezoelectric micro-generator Smart Mater. Struct. 10 850–2
[16] Ericka M, Vasic D, Costa F, Poulin G and Tliba S 2005 Energy harvesting from vibration using a piezoelectric membrane J. Physique IV 128 187–93
[17] Kim S, Clark W W and Wang Q M 2005 Piezoelectric energy harvesting with a clamped circular plate: analysis J. Intell. Mater. Syst. Struct. 16 847–54
[18] Kim S, Clark W W and Wang Q M 2005 Piezoelectric energy harvesting with a clamped circular plate: experimental study J. Intell. Mater. Syst. Struct. 16 855–63
[19] Whalen S, Thompson M, Bahr D, Richards C and Richards R 2003 Design, fabrication and testing of the micro heat engine Sensors Actuators A 104 290–8
[20] Allen J J and Smits A J 2001 Energy harvesting EEL J. Fluids Struct. 15 629–40
[21] Taylor G W, Burns J R, Kammann S M, Powers W B and Welsh T R 2001 The energy harvesting Eel: a small subsurface ocean/river power generator IEEE J. Ocean. Eng. 26 539–47
[22] Kim H W, Batra A, Priya S, Uchino K, Markley D, Newnham R E and Hofmann H F 2004 Energy harvesting using a piezoelectric ‘Cymbal’ transducer in dynamic environment Japan. J. Appl. Phys. 43 6178–83
[23] Kim H W, Priya S, Uchino K and Newnham R E 2005 Piezoelectric energy harvesting under high pre-stressed cyclic vibrations J. Electroceram. 15 27–34
[24] Priya S, Chen C T, Fye D and Zahnd J 2005 Piezoelectric windmill: a novel solution to remote sensing Japan. J. Appl. Phys. 44 L104–7
[25] duToit N E, Wardle B L and Kim S G 2005 Design considerations for MEMS-scale piezoelectric mechanical vibration energy harvesters Integr. Ferroelectr. 71 121–60
[26] F Lu, Lee H P and Lim S P 2004 Modeling and analysis of micro piezoelectric power generators for micro-electromechanical-systems applications Smart Mater. Struct. 13 57–63
[27] Ramsay M J and Clark W W 2001 Piezoelectric energy harvesting for bio MEMS applications Proc. SPIE 4332 429–38
[28] Horowitz S B, Sheplak M, Cattafesta L N III and Nishida T 2006 A MEMS acoustic energy harvester J. Micromech. Microeng. 16 S174–81
[29] Umeda M, Nakamura K and Ueha S 1996 Analysis of the transformation of mechanical impact energy to electric energy using piezoelectric vibrator Japan. J. Appl. Phys. 35 3267–73
[30] Umeda M, Nakamura K and Ueha S 1997 Energy storage characteristics of a piezo-generator using impact induced vibration Japan. J. Appl. Phys. 36 3146–51
[31] Goldfarb M and Jones L D 1999 On the efficiency of electric power generation with piezoelectric ceramic Trans. ASME J. Dynamic Syst. Meas. Control 121 566–71
[32] Roundy S 2005 On the effectiveness of vibration-based energy harvesting J. Intell. Mater. Syst. Struct. 16 809–23
[33] Richards C D, Anderson M J, Bahr D F and Richards R F 2004 Efficiency of energy conversion for devices containing a piezoelectric component J. Micromech. Microeng. 14 717–21
[34] Shu Y C and Lien I C 2006 Analysis of power output for piezoelectric energy harvesting systems Smart Mater. Struct. 15 1499–512
[35] Cho J, Anderson M, Richards R, Bahr D and Richards C 2005 Optimization of electromechanical coupling for a thin-film PZT membrane: I. Modeling J. Micromech. Microeng. 15 1797–803
[36] Cho J, Anderson M, Richards R, Bahr D and Richards C 2005 Optimization of electromechanical coupling for a thin-film PZT membrane: II. Experiment J. Micromech. Microeng. 15 1804–9
[37] Lesieutre G A, Ottman G K and Hofmann H F 2004 Damping as a result of piezoelectric energy harvesting J. Sound Vib. 269 991–1001
[38] Sodano H A, Park G and Inman D J 2004 Estimation of electric charge output for piezoelectric energy harvesting J. Strain 40 49–58
[39] Ottman G K, Hofmann H F, Bhatt A C and Lesieutre G A 2002 Adaptive piezoelectric energy harvesting circuit for wireless remote power supply IEEE Trans. Power Electron. 17 669–76
[40] Ottman G K, Hofmann H F and Lesieutre G A 2003 Optimized piezoelectric energy harvesting circuit using step-down converter in discontinuous conduction mode IEEE Trans. Power Electron. 18 696–703
[41] Badel A, Benayad A, Lefeuvre E, Lebrun L, Richard C and Guyomar D 2006 Single crystals and nonlinear process for outstanding vibration-powered electrical generators IEEE Trans. Ultrason. Ferroelectr. Freq. Control 53 673–84
[42] Badel A, Guyomar D, Lefeuvre E and Richard C 2005 Efficiency enhancement of a piezoelectric energy harvesting device in pulsed operation by synchronous charge inversion J. Intell. Mater. Syst. Struct. 16 889–901
[43] Guyomar D, Badel A and Lefeuvre E 2005 Toward energy harvesting using active materials and conversion improvement by nonlinear processing IEEE Trans. Ultrason. Ferroelectr. Freq. Control 52 584–95
[44] Lefeuvre E, Badel A, Richard C and Guyomar D 2005 Piezoelectric energy harvesting device optimization by synchronous electric charge extraction J. Intell. Mater. Syst. Struct. 16 865–76
[45] Hagood N W, Chung W H and Flotow A V 1990 Modelling of piezoelectric actuator dynamics for active structural control J. Intell. Mater. Syst. Struct. 1 327–54
[46] Wang Q M and Cross L E 1999 Constitutive equations of symmetrical triple layer piezoelectric benders IEEE Trans. Ultrason. Ferroelectr. Freq. Control 46 1343–51
[47] Lesieutre G A and Davis C L 1997 Can a coupling coefficient of a piezoelectric device be higher than those of its active material? J. Intell. Mater. Syst. Struct. 8 859–67
[48] Williams C B and Yates R B 1996 Analysis of a micro-electric generator for microsystems Sensors Actuators A 52 8–11
[49] Rao S S 1990 Mechanical Vibrations (Reading, MA: Addison-Wesley)