# Models for 31-Mode PVDF Energy Harvester for Wearable Applications

Authors: Jingjing Zhao; Zheng You
Journal: The Scientific World Journal
 (2014)

---

## Abstract

Currently, wearable electronics are increasingly widely used, leading to an increasing need of portable power supply. As a clean and renewable power source, piezoelectric energy harvester can transfer mechanical energy into electric energy directly, and the energy harvester based on polyvinylidene difluoride (PVDF) operating in 31-mode is appropriate to harvest energy from human motion. This paper established a series of theoretical models to predict the performance of 31-mode PVDF energy harvester. Among them, the energy storage one can predict the collected energy accurately during the operation of the harvester. Based on theoretical study and experiments investigation, two approaches to improve the energy harvesting performance have been found. Furthermore, experiment results demonstrate the high accuracies of the models, which are better than 95%.


---

## Body

## 1. Introduction

As a power supply, an electrochemical battery is necessary to wearable electronics, which is one focus of wearable electronics. But two obvious disadvantages of these batteries, limited energy storage capability and potential damage to environment, have limited further development and applications of wearable electronics. Many efforts have been focused on wearable energy harvesters, which can harvest the mechanical energy dissipated from human motion to supply renewable and clean energy [1]. There have been several concepts of harvesting energy from human motion with varied mechanisms to be studied [2–13], including electromagnetic [5–8], electrostatic [9], thermoelectric [10], and piezoelectric harvesters [11–13]. Among them, piezoelectric energy harvester can convert mechanical energy into electric power directly, resulting in a more compact and simpler structure. Lead zirconate titanate (PZT) and PVDF are two common functional materials for piezoelectric energy harvesters. Compared with rigid, brittle, and heavy PZT, PVDF has considerable flexibility, good stability, and easiness to handle and shape [4]. Considering that human motion is characterized by high amplitude and low frequency (less than 5 Hz typically [14, 15]), the advantages presented above make PVDF more appropriate for wearable applications where flexibility is necessary, compared with PZT. PVDF has been designed for wearable energy harvesters implemented in shoes [16–19], bags [20, 21], and clothing [22, 23]. Almost all the reported PVDF energy harvesters are based on 31-mode, because 31-mode can easily yield larger strain by smaller forces in comparison with 33-mode.As discussed above, 31-mode PVDF energy harvester for wearable applications has been studied widely. But there is still a lack within the general theoretical models for design, making it hard to analyze and optimize the energy harvester. So this paper focuses on a series of models, including basic model, resistive load model, capacitive load model, and energy storage model, which can give theoretical support to this kind of harvester. Design and analysis of the device with accuracy and facility can be realized by using the models. Experiments show that the models were accurate with errors less than 5%. Furthermore, with the help of the models and experimental results, it has been found that two approaches can improve the performance of energy harvesting.


## 2. Models for 31-Mode PVDF Energy Harvester

### 2.1. PVDF Energy Harvester Unit

A PVDF energy harvester unit studied in this paper is shown in Figure1(a), which is made of a layer of PVDF film. A cartesian coordinate system is established to identify the directions within the PVDF unit. 1-axis and 3-axis are mechanical stress direction and polarization direction during PVDF fabrication, respectively. Electrodes are placed on the surfaces perpendicular to 3-axis and U
O is the output voltage between the two electrodes. l, w, and h are length, width, and thickness of the PVDF unit. A
3 is the area of one 3-axis surface or one electrode, and A
1 is the area of cross-section along 1-axis. The PVDF energy harvester unit can be used individually, or several units are stacked and connected in series or in parallel for higher output voltage or current. Figure 2(a) shows the typical equivalent circuit used to present the operation of the unit. Q
F is an equivalent charge generator representing the piezoelectric charges. R
P and C
P are the resistance and capacitance of the PVDF unit, and U
O is the output voltage. R
P is extremely large with little current flowing, therefore considering R
P infinite. Thus, the equivalent can be simplified as shown in Figure 2(b).(a) A PVDF energy harvester unit, (b) 31-mode, and (c) 33-mode.
(a)
(b)
(c)(a) The equivalent circuit for a PVDF energy harvester unit, (b) The ideal equivalent circuit.
(a)
(b)31-mode and 33-mode are two common modes in which piezoelectric energy harvesters operate generally. A force is applied along 1-axis then charges accumulate on the 3-axis electrodes when the harvester unit works in 31-mode, and the applied force direction parallels 3-axis then charges generated on the 3-axis electrodes in 33-mode, as shown in Figures1(b) and 1(c). Benefiting from A
1 which is much smaller than A
3, the PVDF unit operated in 31-mode can gain a far larger strain comparing to the operation in 33-mode, when the applied forces have the same amplitude. Therefore, more charges can be produced in 31-mode. Furthermore, the forces along 1-axis are very common in wearable applications, which can be produced by limbs motion [23]. Based on the analysis above, the models for PVDF energy harvester operated in 31-mode are more valuable and this paper focuses on them.The force along 1-axis produced by human motion (F
1) is an alternating force with a low frequency less than 5 Hz. The abstract model of F
1 is shown in Figure 3, and during one cycle the force can be divided into four stages: rising, roof, falling, and pause. F
n is the amplitude of the Nth cycle, which is a key parameter for calculating the quantity of stored energy or charges. Basic model, resistive model, capacitive load model, and energy storage model of the PVDF unit are established step by step as below.Figure 3
The abstract model ofF
1 which can be divided into four stages during one cycle.


### 2.2. Basic Model

The basic model describes the piezoelectric performance of a PVDF unit operated in 31-mode. And the piezoelectric constitutive equations for 31-mode are(1)
S
1
=
s
11
E
T
1
+
d
31
E
3
,
D
3
=
d
31
T
1
+
ε
33
T
E
3
,
where S
1 is the strain, s
11
E is the elastic compliance under constant electric field, T
1 is the stress, d
31 is the piezoelectric constant, E
3 is electric field, D
3 is the electric displacement, and ε
33
T is the permittivity under constant stress. F
1 is produced by human motion and dynamic, leading to inertial force and damping force in PVDF unit. Considering that the frequency of the human motion is low enough and the weight of PVDF unit is light, the values of inertial and damping forces can be neglected. The constitutive equations can be simplified by the following relationships [20]:
(2)
D
3
(
t
)
=
Q
3
(
t
)
A
3
=
Q
3
(
t
)
l
w
,
E
3
(
t
)
=
U
O
(
t
)
h
,
T
1
(
t
)
=
F
1
(
t
)
A
1
=
F
1
(
t
)
w
h
,
S
1
(
t
)
=
y
(
t
)
l
,
where Q
3 is the free charges collected on the electrodes (equal to the charges flowing to the external circuit), U
O is the output voltage, F
1 is the dynamic force applied on PVDF along 1-axis, and y is the displacement along 1-axis. The piezoelectric constitutive equations can be rewritten as
(3)
F
1
(
t
)
=
F
11
(
t
)
+
F
12
(
t
)
,
Q
F
(
t
)
=
d
31
l
h
F
1
(
t
)
=
Q
3
(
t
)
+
C
P
U
O
(
t
)
,
F
11
(
t
)
=
Y
A
1
l
y
(
t
)
,
F
12
(
t
)
=
d
31
Y
w
U
O
(
t
)
,
C
P
=
ε
33
A
3
h
,
where Y is the PVDF elastic modulus, Q
F is the generated charges, C
P is the capacitance of PVDF unit, F
11 is the deformation force, and F
12 is the force caused by the secondary converse piezoelectric effect. F
12 reaches maximum when the PVDF unit is open without no load (Q
3
=
0), and the maximum value is far less than F
1, as expressed in (4). Hence, F
12 can be neglected and (3) can be simplified as (5)
(4)
F
12
(
t
)
=
d
31
2
Y
ε
33
F
1
(
t
)
≪
F
1
(
t
)
,
(5)
F
1
(
t
)
=
Y
A
1
l
y
(
t
)
,
Q
F
(
t
)
=
d
31
l
h
F
1
(
t
)
=
Q
3
(
t
)
+
C
P
U
O
(
t
)
.
The model, expressed by (5), describes the piezoelectric performance of PVDF unit operated in 31-mode, and it also acts as the basic model in order to give the fundamental supports to other models discussed below.


### 2.3. Resistive Load Model

The schematic diagram of the equivalent circuit is transferred as shown in Figure4, when the PVDF unit is terminated with a resistive load R
L. R
W presents the wire resistance which is much lower than R
L typically and can be ignored. The applied force is defined in (6), consisting of a sinusoidal force component and a preload force component for convenience of experimental verifications. The amplitude of the preload force is greater than that of the sinusoidal one in order to keep the film unit tensional. Equation (5) can be rewritten as (7). The output voltage is given in (8) and the power output can be calculated by U
O
2
/
R
L. The resistive load model offers an approach to test the piezoelectric performance of PVDF unit operating in 31-mode by measuring the output voltage, and it can further evaluate the generating capacity by calculating the output power
(6)
F
(
t
)
=
A
+
B
sin
(
ω
t
)
A
>
B
>
0
,
(7)
F
1
(
t
)
=
Y
A
1
l
y
(
t
)
=
B
sin
(
ω
t
)
,
d
31
l
h
F
1
(
t
)
=
Q
3
(
t
)
+
C
P
R
L
Q
˙
3
(
t
)
,
(8)
U
O
(
t
)
=
1
(
1
/
ω
C
P
R
L
)
2
+
1
×
d
31
B
ε
33
w
sin
[
ω
t
+
arctan
⁡
(
1
ω
C
P
R
L
)
]
.Figure 4
PVDF unit is connected with a resistive load.


### 2.4. Capacitive Load Model

Figure5 shows the equivalent circuit of the capacitive load model in which the PVDF unit is connected with a capacitor to store the generated charges. The applied force is composed of rising and roof stages expressed in Figure 3 in order to keep the output voltage direction unchanged. The applied force is shown in Figure 6, where T
R is the force rise time, F is the amplitude value, C
L is the storage capacitor, and C
P is the capacitance of PVDF unit.Figure 5
PVDF unit with capacitive load.Figure 6
The applied force is composed of rising and roof stages.By using Kirchhoff’s voltage law and Kirchhoff’s current law, (9) can be given as
(9)
∫
0
t
i
L
(
t
)
d
t
+
∫
0
t
i
P
(
t
)
d
t
=
Q
F
(
t
)
,
∫
0
t
i
P
(
t
)
d
t
C
P
=
∫
0
t
i
L
(
t
)
d
t
C
L
+
R
W
i
L
(
t
)
+
Q
L
(
0
)
C
L
,
Q
L
(
t
)
=
∫
0
t
i
L
(
t
)
d
t
+
Q
L
(
0
)
,
Q
F
(
t
)
=
d
31
l
h
F
1
(
t
)
,
nnnnnnnnnhhnhhnnn
0
≤
t
≤
T
R
,
where R
W is the wire resistance, i
P and i
L are branch currents, Q
L is the charges stored in C
L, and Q
L(0) is the initial value. Equation (10) is obtained by Laplace transform
(10)
Q
L
(
s
)
=
C
L
Q
F
(
s
)
s
C
L
C
P
R
W
+
(
C
L
+
C
P
)
+
C
L
2
C
P
R
W
U
L
(
0
)
s
C
L
C
P
R
W
+
(
C
L
+
C
P
)
+
C
L
2
U
L
(
0
)
s
2
C
L
C
P
R
W
+
s
(
C
L
+
C
P
)
.
The inverse Laplace transform of (10) is (11). Q
L
1 is the zero input response and Q
L
2 is the zero state response. Consider
(11)
Q
L
(
t
)
=
Q
L
1
(
t
)
+
Q
L
2
(
t
)
(
0
≤
t
≤
T
R
)
,
Q
L
1
(
t
)
=
C
L
Q
F
(
t
)
C
L
+
C
P
-
C
L
e
-
t
/
τ
∫
0
t
e
t
/
τ
d
Q
F
(
t
)
,
Q
L
2
(
t
)
=
e
-
t
/
τ
Q
L
(
0
)
+
C
L
Q
L
(
0
)
C
L
+
C
P
(
1
-
e
-
t
/
τ
)
,
τ
=
C
L
C
P
R
W
C
L
+
C
P
,
where τ is the time constant. Typically, the wire resistance R
W approaches zero, and C
L is much greater than C
P. Hence, τ is approximately equal to zero and can be neglected when comparing to T
R. Equation (11) can be simplified to
(12)
Q
L
(
t
)
=
Q
L
1
(
t
)
+
Q
L
2
(
t
)
,
Q
L
1
(
t
)
≈
C
L
Q
F
(
t
)
C
L
+
C
P
,
Q
L
2
(
t
)
≈
C
L
Q
L
(
0
)
C
L
+
C
P
,
nnnhh
0
≤
t
≤
T
R
.
When the applied force F
1 reaches the amplitude F, the charges on the C
P and C
L can be given by (13) according to the law of conservation of charge:
(13)
Q
L
(
t
)
=
C
L
C
L
+
C
P
[
Q
F
(
T
R
)
+
Q
L
(
0
)
]
,
Q
P
(
t
)
=
C
P
C
L
+
C
P
[
Q
F
(
T
R
)
+
Q
L
(
0
)
]
,
Q
F
(
T
R
)
=
d
31
l
h
F
.
The numerical relationship is established between the quantity of stored charge and the amplitude of applied force, which gives supports for the energy storage model.


### 2.5. Energy Storage Model

This model can predict the quantities of the collected charges and stored energy and help to design and optimize a 31-mode PVDF energy harvester. The energy storage circuit with a PVDF unit is shown in Figure7(a). The circuit is widely used for piezoelectric harvester, consisted of a bridge rectifier and a storage capacitor. The rectifier alternates the generated current, and the charges accumulate on the storage capacitor C
L. When two diodes are conducting, the equivalent circuit is presented in Figure 7(b), which is similar to the circuit discussed in Section 2.3. U
D is the diode forward voltage, and R
D is the diode forward resistance. The time constant τ is defined as
(14)
τ
=
C
L
C
P
(
R
W
+
2
R
D
)
C
L
+
C
P
.(a) The circuit for energy storage and (b) the equivalent circuit when bridge rectifier conducting.
(a)
(b)The increase of the charges stored inC
L during one cycle is determined by the force amplitude of every cycle according to the analysis in Section 2.4. Therefore, using the abstract force in Figure 3 as the applied force to analyze the charging process is reasonable. The force period is far greater than τ; therefore the delay between Q
L
(
t
) and Q
F
(
t
) can be neglected. The equations obtained in Section 2.4 can be used here. The four stages in the Nth cycle are analyzed below, and the charging process of C
L is shown in Figure 8.Figure 8
The charging process ofC
L during one cycle.Rising Stage (
t
n
,
1~t
n
,
2
). Rising stage can be divided into two parts by t
D which is the first rectifier turn-on time. The voltage on C
P is lower than rectifier threshold voltage and the generated charges accumulate in C
P before t
D. After t
D, the current i
L flows through the rectifier and charges are stored in C
P and C
L. At the moment t
n
,
2, the force stops rising and reaches the amplitude F
n of the Nth cycle. This stage can be presented by the following equations, where U
P is the voltage across C
P, and U
L is the voltage across C
L:
(15)
t
n
1
≤
t
≤
t
D


:
{
U
L
(
t
)
=
U
L
(
t
n
1
)
U
P
(
t
)
=
U
P
(
t
n
1
)
+
Q
F
n
(
t
)
C
P
,
t
=
t
D


:
{
U
P
(
t
D
)
=
U
L
(
t
n
1
)
+
2
U
D
Q
F
n
(
t
D
)
=
[
U
L
(
t
n
1
)
+
2
U
D
-
U
P
(
t
n
1
)
]
C
P
,
t
D
≤
t
≤
t
n
2


:
{
U
L
(
t
)
=
U
L
(
t
n
1
)
+
Q
F
n
(
t
)
-
Q
F
n
(
t
D
)
C
P
+
C
L
U
P
(
t
)
=
U
L
(
t
n
1
)
+
2
U
D
+
Q
F
n
(
t
)
-
Q
F
n
(
t
D
)
C
P
+
C
L
,
t
=
t
n
2


:
{
Q
L
(
t
n
2
)
=
C
L
[
Q
F
n
(
t
n
2
)
+
Q
L
(
t
n
1
)
+
Q
P
(
t
n
1
)
]
C
P
+
C
L
Q
L
(
t
n
2
)
-
2
C
L
C
P
U
D
C
P
+
C
L
Q
P
(
t
n
2
)
=
C
P
[
Q
F
n
(
t
n
2
)
+
Q
L
(
t
n
1
)
+
Q
P
(
t
n
1
)
]
C
P
+
C
L
Q
L
(
t
n
2
)
+
2
C
L
C
P
U
D
C
P
+
C
L
Q
F
n
(
t
n
2
)
=
d
31
l
h
F
n
.Roof Stage (
t
n
,
2~t
n
,
3
). During roof stage, the applied force keeps unchanged, so no charge is generated and the rectifier turns off. The leakage current of capacitors is quite low and the stage’s duration is short. It can be concluded that no charge on C
P and C
L leaks away, so equations can be given as
(16)
t
n
2
≤
t
≤
t
n
3


:
{
Q
L
(
t
)
=
Q
L
(
t
n
2
)
Q
P
(
t
)
=
Q
P
(
t
n
2
)
.
Falling Stage (
t
n
,
3~t
n
,
4
). The input force decreases to zero from F
n. This stage is similar to the rising stage and can be divided into two parts by t
D
′, the rectifier turn-on time during this stage. By referring to the analysis of rising stage, the charging process can be presented by the following equations:
(17)
t
n
3
≤
t
≤
t
D
′


:
{
U
L
(
t
)
=
U
L
(
t
n
3
)
U
P
(
t
)
=
U
P
(
t
n
3
)
+
Q
F
n
(
t
)
C
P
,
t
=
t
D
′


:
{
U
P
(
t
D
′
)
=
-
[
U
L
(
t
n
3
)
+
2
U
D
]
Q
F
n
(
t
D
′
)
=
-
[
U
L
(
t
n
3
)
+
2
U
D
+
U
P
(
t
n
3
)
]
C
P
,
t
n
3
≤
t
≤
t
n
4


:
{
U
L
(
t
)
=
U
L
(
t
n
3
)
-
Q
F
n
(
t
)
-
Q
F
n
(
t
D
′
)
C
P
+
C
L
U
P
(
t
)
=
-
[
U
L
(
t
n
3
)
+
2
U
D
-
Q
F
n
(
t
)
-
Q
F
n
(
t
D
′
)
C
P
+
C
L
]
,
t
=
t
n
4


:
{
Q
L
(
t
n
4
)
=
C
L
[
Q
L
(
t
n
3
)
-
Q
F
n
(
t
n
4
)
-
Q
P
(
t
n
3
)
]
C
P
+
C
L
Q
L
(
t
n
4
)
-
2
C
L
C
P
U
D
C
P
+
C
L
Q
P
(
t
n
4
)
=
-
C
P
[
Q
L
(
t
n
3
)
-
Q
F
n
(
t
n
4
)
-
Q
P
(
t
n
3
)
]
C
P
+
C
L
Q
L
(
t
n
4
)
-
2
C
L
C
P
U
D
C
P
+
C
L
Q
F
n
(
t
n
4
)
=
-
d
31
l
h
F
n
.Pause Stage (
t
n
,
4~t
n
+
1,1
). No force is applied on PVDF until the (
N
+
1
)th cycle begins. The charges on C
P and C
L remain constant for the same reasons presented in stay stage. Consider
(18)
t
n
4
≤
t
≤
t
(
n
+
1
)
1


:
{
Q
L
(
t
)
=
Q
L
(
t
n
4
)
Q
P
(
t
)
=
Q
P
(
t
n
4
)
.
As discussed above, the charges are generated and collected on C
L during rising and falling stages, and the applied force keeps constant with no charge generated in roof and pause stages. Recursion formulae of Q
L and Q
P can be derived and be presented by
(19)
Q
L
,
0
=
Q
L
(
N
)
,
Q
P
,
0
=
Q
P
(
N
)
,
Q
L
,
1
=
t
1
(
Q
F
+
Q
L
,
0
+
Q
P
,
0
)
-
2
t
1
U
D
C
P
,
Q
P
,
1
=
t
2
(
Q
F
+
Q
L
,
0
+
Q
P
,
0
)
+
2
t
1
U
D
C
P
,
Q
L
,
2
=
t
1
(
Q
F
+
Q
L
,
1
-
Q
P
,
1
)
-
2
t
1
U
D
C
P
,
Q
P
,
2
=
-
t
2
(
Q
F
+
Q
L
,
1
-
Q
P
,
1
)
-
2
t
2
U
D
C
L
,
t
1
=
C
L
(
C
L
+
C
P
)
,
t
2
=
C
P
(
C
L
+
C
P
)
,
Q
F
=
F
n
d
31
l
h
,
Q
L
(
N
+
1
)
=
Q
L
,
2
,
Q
P
(
N
+
1
)
=
Q
P
,
2
,
where N is the cycle number, Q
L
(
N
) is the quantity of the charges on C
L at the end of the Nth cycle, Q
P
(
N
) is the quantity of charges on C
P at the end of the Nth cycle, F
n is the force amplitude in the Nth cycle, and t
1 and t
2 are two defined parameters for simplification. Q
L
,
0 and Q
P
,
0 are the initial electric quantities of one cycle, Q
L
,
1 and Q
P
,
1 are the electric quantities when the applied force reaches F
n, and Q
L
,
2 and Q
P
,
2 are the electric quantities when the force returns to its initial value. If the force amplitude in each cycle can be acquired by test or calculation, the recursion formulae can be utilized to estimate the charges collected on storage capacitor. In some cases, the force amplitude in each cycle is a constant value F and C
P as well as C
L stores no charge initially; then (19) can be rewritten below as
(20)
Q
L
(
N
+
1
)
=
A
1
B
+
A
+
-
1
A
+
n
-
A
1
B
+
A
+
-
1
+
B
1
,
Q
P
(
N
+
1
)
=
A
2
B
+
A
+
-
1
A
+
n
-
A
2
B
+
A
+
-
1
+
B
2
,
(21)
A
1
=
(
t
1
-
t
2
)
t
1
,
A
2
=
(
t
2
-
t
1
)
t
2
,
B
1
=
-
4
t
1
2
U
D
C
P
-
2
t
1
U
D
C
P
+
2
t
1
2
Q
F
,
B
2
=
4
t
1
t
2
U
D
C
P
-
2
t
1
U
D
C
P
-
2
t
1
t
2
Q
F
,
A
+
=
A
1
+
A
2
=
(
t
1
-
t
2
)
2
,
B
+
=
B
1
+
B
2
=
2
Q
F
(
t
1
-
t
2
)
t
1
-
8
t
1
2
U
D
C
P
,
Q
F
=
F
d
31
l
h
.
It is known that C
L is much greater than C
P, and U
D is not more than 0.7 V; thus t
1 is approximately 1 and t
2 is close to 0. And the approximate expression of (20) is given below as
(22)
Q
L
(
N
+
1
)
≈
2
Q
F
+
2
Q
F
1
-
A
+
N
1
-
A
+
,
Q
P
(
N
+
1
)
≈
0
.
According to (20) and (22), it is shown that C
P and C
L are main factors that influence the charge collection, while U
D has little influence on Q
L, so it is not an important parameter that needs selecting carefully. Hence, researchers can focus on key parameters by analyzing the energy storage model. Besides, Q
L is approximately proportional to Q
F while Q
F is proportional to the applied force amplitude; therefore Q
L is approximately proportional to the force amplitude. In reality, C
L is far larger than C
P; therefore most of the generated charges are transferred to C
L. The charge increment per cycle declines as the cycle number increases, as shown in (23). The charge transfer efficiency η
Q
L is defined in (24). The energy stored on C
L is given in (25), and (26) presents the energy increment per cycle. Consider
(23)
Δ
Q
L
(
N
)
=
Q
L
(
N
)
-
Q
L
(
N
-
1
)
=
{
A
1
B
+
A
+
N
-
1
N
≠
1
B
1
≈
2
Q
F
N
=
1
,
(24)
η
Q
L
(
N
)
=
Δ
Q
L
(
N
)
2
Q
F
,
(25)
E
L
(
N
)
=
Q
L
(
N
)
2
2
C
L
,
(26)
Δ
E
L
(
N
)
=
E
L
(
N
)
-
E
L
(
N
-
1
)
.
The energy storage model is based on the energy storage circuit shown in Figure 7. The recursion formulae of stored charges are presented by (19), which can predict the charge quantity when the amplitude of applied force in each cycle is known. If the applied force keeps constant, the energy storage model can be presented by (20). Several parameters shown in (23)~(25) can evaluate the energy harvesting performance.


### 2.6. Two Approaches to Improve Energy Harvesting Performance

On the basis of the energy storage model, two approaches to improve energy harvesting performance are derived. The first approach is that the voltage on storage capacitor should to be kept low for high energy storage efficiency. Several parameters of the energy storage model are assigned to detail for convenience to clear analysis and easy demonstration. ConsiderC
L
=
112 μF, C
P
=
133 nF, Q
F
=
43.3 μC, and U
D
=
0.53 V; then Q
L, η
Q
L, E
L, and Δ
E
L are simulated in Figure 9. Apparently, η
Q
L decreases exponentially as the cycle number N increases. Δ
E
L first increases and then decreases with an increasing N. It is because that as the voltage on C
L increases, more charges are dissipated on C
P. The stored energy on C
L should be moved to another energy storage device, like a battery, or consumed by load timely before Δ
E
L begins to decrease or η
Q
L is lower than a certain value (85% for instance). After transferring the stored energy, the voltage on C
L will decrease and C
L can store energy with a high efficiency again. In conclusion, reducing the voltage across the storage capacitor is a viable approach for high energy harvesting performance.(a) The charges and energy stored on storage capacitor versus the cycle number. (b) The charges transfer efficiency and energy increment per cycle versus the cycle number.
(a)
(b)The second approach is to improve the generated chargesQ
F. The electrical connection configuration of PVDF units affects Q
F considerably, and three common configurations are studied. The first configuration is one h thick PVDF unit, and the other two configurations are n layers of h
/
n thick PVDF units that are stacked and wired in parallel and series. In the three cases, PVDF units have the same length, width, and the total thickness. The generated charges of the three cases can be given by (27) and the PVDF equivalent capacitance is expressed by (28), where S means single PVDF unit, n
P means n PVDF units in parallel, and n
S means n PVDF units in series. Obviously, thin PVDF units with a parallel connection can generate most charges, proving that connecting thin PVDF units in parallel is an effective approach to improving the energy harvesting performance. Consider
(27)
Q
F
,
S
=
d
31
l
h
F
,
Q
F
,
n
P
=
n
d
31
l
h
F
=
n
Q
F
,
S
,
Q
F
,
n
S
=
d
31
l
h
F
=
Q
F
,
S
,
(28)
C
P
,
S
=
ε
33
A
3
h
,
C
P
,
n
P
=
n
ε
33
A
3
h
/
n
=
n
2
C
P
,
S
,
C
P
,
n
S
=
1
n
ε
33
A
3
h
/
n
=
C
P
,
S
.


## 2.1. PVDF Energy Harvester Unit

A PVDF energy harvester unit studied in this paper is shown in Figure1(a), which is made of a layer of PVDF film. A cartesian coordinate system is established to identify the directions within the PVDF unit. 1-axis and 3-axis are mechanical stress direction and polarization direction during PVDF fabrication, respectively. Electrodes are placed on the surfaces perpendicular to 3-axis and U
O is the output voltage between the two electrodes. l, w, and h are length, width, and thickness of the PVDF unit. A
3 is the area of one 3-axis surface or one electrode, and A
1 is the area of cross-section along 1-axis. The PVDF energy harvester unit can be used individually, or several units are stacked and connected in series or in parallel for higher output voltage or current. Figure 2(a) shows the typical equivalent circuit used to present the operation of the unit. Q
F is an equivalent charge generator representing the piezoelectric charges. R
P and C
P are the resistance and capacitance of the PVDF unit, and U
O is the output voltage. R
P is extremely large with little current flowing, therefore considering R
P infinite. Thus, the equivalent can be simplified as shown in Figure 2(b).(a) A PVDF energy harvester unit, (b) 31-mode, and (c) 33-mode.
(a)
(b)
(c)(a) The equivalent circuit for a PVDF energy harvester unit, (b) The ideal equivalent circuit.
(a)
(b)31-mode and 33-mode are two common modes in which piezoelectric energy harvesters operate generally. A force is applied along 1-axis then charges accumulate on the 3-axis electrodes when the harvester unit works in 31-mode, and the applied force direction parallels 3-axis then charges generated on the 3-axis electrodes in 33-mode, as shown in Figures1(b) and 1(c). Benefiting from A
1 which is much smaller than A
3, the PVDF unit operated in 31-mode can gain a far larger strain comparing to the operation in 33-mode, when the applied forces have the same amplitude. Therefore, more charges can be produced in 31-mode. Furthermore, the forces along 1-axis are very common in wearable applications, which can be produced by limbs motion [23]. Based on the analysis above, the models for PVDF energy harvester operated in 31-mode are more valuable and this paper focuses on them.The force along 1-axis produced by human motion (F
1) is an alternating force with a low frequency less than 5 Hz. The abstract model of F
1 is shown in Figure 3, and during one cycle the force can be divided into four stages: rising, roof, falling, and pause. F
n is the amplitude of the Nth cycle, which is a key parameter for calculating the quantity of stored energy or charges. Basic model, resistive model, capacitive load model, and energy storage model of the PVDF unit are established step by step as below.Figure 3
The abstract model ofF
1 which can be divided into four stages during one cycle.


## 2.2. Basic Model

The basic model describes the piezoelectric performance of a PVDF unit operated in 31-mode. And the piezoelectric constitutive equations for 31-mode are(1)
S
1
=
s
11
E
T
1
+
d
31
E
3
,
D
3
=
d
31
T
1
+
ε
33
T
E
3
,
where S
1 is the strain, s
11
E is the elastic compliance under constant electric field, T
1 is the stress, d
31 is the piezoelectric constant, E
3 is electric field, D
3 is the electric displacement, and ε
33
T is the permittivity under constant stress. F
1 is produced by human motion and dynamic, leading to inertial force and damping force in PVDF unit. Considering that the frequency of the human motion is low enough and the weight of PVDF unit is light, the values of inertial and damping forces can be neglected. The constitutive equations can be simplified by the following relationships [20]:
(2)
D
3
(
t
)
=
Q
3
(
t
)
A
3
=
Q
3
(
t
)
l
w
,
E
3
(
t
)
=
U
O
(
t
)
h
,
T
1
(
t
)
=
F
1
(
t
)
A
1
=
F
1
(
t
)
w
h
,
S
1
(
t
)
=
y
(
t
)
l
,
where Q
3 is the free charges collected on the electrodes (equal to the charges flowing to the external circuit), U
O is the output voltage, F
1 is the dynamic force applied on PVDF along 1-axis, and y is the displacement along 1-axis. The piezoelectric constitutive equations can be rewritten as
(3)
F
1
(
t
)
=
F
11
(
t
)
+
F
12
(
t
)
,
Q
F
(
t
)
=
d
31
l
h
F
1
(
t
)
=
Q
3
(
t
)
+
C
P
U
O
(
t
)
,
F
11
(
t
)
=
Y
A
1
l
y
(
t
)
,
F
12
(
t
)
=
d
31
Y
w
U
O
(
t
)
,
C
P
=
ε
33
A
3
h
,
where Y is the PVDF elastic modulus, Q
F is the generated charges, C
P is the capacitance of PVDF unit, F
11 is the deformation force, and F
12 is the force caused by the secondary converse piezoelectric effect. F
12 reaches maximum when the PVDF unit is open without no load (Q
3
=
0), and the maximum value is far less than F
1, as expressed in (4). Hence, F
12 can be neglected and (3) can be simplified as (5)
(4)
F
12
(
t
)
=
d
31
2
Y
ε
33
F
1
(
t
)
≪
F
1
(
t
)
,
(5)
F
1
(
t
)
=
Y
A
1
l
y
(
t
)
,
Q
F
(
t
)
=
d
31
l
h
F
1
(
t
)
=
Q
3
(
t
)
+
C
P
U
O
(
t
)
.
The model, expressed by (5), describes the piezoelectric performance of PVDF unit operated in 31-mode, and it also acts as the basic model in order to give the fundamental supports to other models discussed below.


## 2.3. Resistive Load Model

The schematic diagram of the equivalent circuit is transferred as shown in Figure4, when the PVDF unit is terminated with a resistive load R
L. R
W presents the wire resistance which is much lower than R
L typically and can be ignored. The applied force is defined in (6), consisting of a sinusoidal force component and a preload force component for convenience of experimental verifications. The amplitude of the preload force is greater than that of the sinusoidal one in order to keep the film unit tensional. Equation (5) can be rewritten as (7). The output voltage is given in (8) and the power output can be calculated by U
O
2
/
R
L. The resistive load model offers an approach to test the piezoelectric performance of PVDF unit operating in 31-mode by measuring the output voltage, and it can further evaluate the generating capacity by calculating the output power
(6)
F
(
t
)
=
A
+
B
sin
(
ω
t
)
A
>
B
>
0
,
(7)
F
1
(
t
)
=
Y
A
1
l
y
(
t
)
=
B
sin
(
ω
t
)
,
d
31
l
h
F
1
(
t
)
=
Q
3
(
t
)
+
C
P
R
L
Q
˙
3
(
t
)
,
(8)
U
O
(
t
)
=
1
(
1
/
ω
C
P
R
L
)
2
+
1
×
d
31
B
ε
33
w
sin
[
ω
t
+
arctan
⁡
(
1
ω
C
P
R
L
)
]
.Figure 4
PVDF unit is connected with a resistive load.


## 2.4. Capacitive Load Model

Figure5 shows the equivalent circuit of the capacitive load model in which the PVDF unit is connected with a capacitor to store the generated charges. The applied force is composed of rising and roof stages expressed in Figure 3 in order to keep the output voltage direction unchanged. The applied force is shown in Figure 6, where T
R is the force rise time, F is the amplitude value, C
L is the storage capacitor, and C
P is the capacitance of PVDF unit.Figure 5
PVDF unit with capacitive load.Figure 6
The applied force is composed of rising and roof stages.By using Kirchhoff’s voltage law and Kirchhoff’s current law, (9) can be given as
(9)
∫
0
t
i
L
(
t
)
d
t
+
∫
0
t
i
P
(
t
)
d
t
=
Q
F
(
t
)
,
∫
0
t
i
P
(
t
)
d
t
C
P
=
∫
0
t
i
L
(
t
)
d
t
C
L
+
R
W
i
L
(
t
)
+
Q
L
(
0
)
C
L
,
Q
L
(
t
)
=
∫
0
t
i
L
(
t
)
d
t
+
Q
L
(
0
)
,
Q
F
(
t
)
=
d
31
l
h
F
1
(
t
)
,
nnnnnnnnnhhnhhnnn
0
≤
t
≤
T
R
,
where R
W is the wire resistance, i
P and i
L are branch currents, Q
L is the charges stored in C
L, and Q
L(0) is the initial value. Equation (10) is obtained by Laplace transform
(10)
Q
L
(
s
)
=
C
L
Q
F
(
s
)
s
C
L
C
P
R
W
+
(
C
L
+
C
P
)
+
C
L
2
C
P
R
W
U
L
(
0
)
s
C
L
C
P
R
W
+
(
C
L
+
C
P
)
+
C
L
2
U
L
(
0
)
s
2
C
L
C
P
R
W
+
s
(
C
L
+
C
P
)
.
The inverse Laplace transform of (10) is (11). Q
L
1 is the zero input response and Q
L
2 is the zero state response. Consider
(11)
Q
L
(
t
)
=
Q
L
1
(
t
)
+
Q
L
2
(
t
)
(
0
≤
t
≤
T
R
)
,
Q
L
1
(
t
)
=
C
L
Q
F
(
t
)
C
L
+
C
P
-
C
L
e
-
t
/
τ
∫
0
t
e
t
/
τ
d
Q
F
(
t
)
,
Q
L
2
(
t
)
=
e
-
t
/
τ
Q
L
(
0
)
+
C
L
Q
L
(
0
)
C
L
+
C
P
(
1
-
e
-
t
/
τ
)
,
τ
=
C
L
C
P
R
W
C
L
+
C
P
,
where τ is the time constant. Typically, the wire resistance R
W approaches zero, and C
L is much greater than C
P. Hence, τ is approximately equal to zero and can be neglected when comparing to T
R. Equation (11) can be simplified to
(12)
Q
L
(
t
)
=
Q
L
1
(
t
)
+
Q
L
2
(
t
)
,
Q
L
1
(
t
)
≈
C
L
Q
F
(
t
)
C
L
+
C
P
,
Q
L
2
(
t
)
≈
C
L
Q
L
(
0
)
C
L
+
C
P
,
nnnhh
0
≤
t
≤
T
R
.
When the applied force F
1 reaches the amplitude F, the charges on the C
P and C
L can be given by (13) according to the law of conservation of charge:
(13)
Q
L
(
t
)
=
C
L
C
L
+
C
P
[
Q
F
(
T
R
)
+
Q
L
(
0
)
]
,
Q
P
(
t
)
=
C
P
C
L
+
C
P
[
Q
F
(
T
R
)
+
Q
L
(
0
)
]
,
Q
F
(
T
R
)
=
d
31
l
h
F
.
The numerical relationship is established between the quantity of stored charge and the amplitude of applied force, which gives supports for the energy storage model.


## 2.5. Energy Storage Model

This model can predict the quantities of the collected charges and stored energy and help to design and optimize a 31-mode PVDF energy harvester. The energy storage circuit with a PVDF unit is shown in Figure7(a). The circuit is widely used for piezoelectric harvester, consisted of a bridge rectifier and a storage capacitor. The rectifier alternates the generated current, and the charges accumulate on the storage capacitor C
L. When two diodes are conducting, the equivalent circuit is presented in Figure 7(b), which is similar to the circuit discussed in Section 2.3. U
D is the diode forward voltage, and R
D is the diode forward resistance. The time constant τ is defined as
(14)
τ
=
C
L
C
P
(
R
W
+
2
R
D
)
C
L
+
C
P
.(a) The circuit for energy storage and (b) the equivalent circuit when bridge rectifier conducting.
(a)
(b)The increase of the charges stored inC
L during one cycle is determined by the force amplitude of every cycle according to the analysis in Section 2.4. Therefore, using the abstract force in Figure 3 as the applied force to analyze the charging process is reasonable. The force period is far greater than τ; therefore the delay between Q
L
(
t
) and Q
F
(
t
) can be neglected. The equations obtained in Section 2.4 can be used here. The four stages in the Nth cycle are analyzed below, and the charging process of C
L is shown in Figure 8.Figure 8
The charging process ofC
L during one cycle.Rising Stage (
t
n
,
1~t
n
,
2
). Rising stage can be divided into two parts by t
D which is the first rectifier turn-on time. The voltage on C
P is lower than rectifier threshold voltage and the generated charges accumulate in C
P before t
D. After t
D, the current i
L flows through the rectifier and charges are stored in C
P and C
L. At the moment t
n
,
2, the force stops rising and reaches the amplitude F
n of the Nth cycle. This stage can be presented by the following equations, where U
P is the voltage across C
P, and U
L is the voltage across C
L:
(15)
t
n
1
≤
t
≤
t
D


:
{
U
L
(
t
)
=
U
L
(
t
n
1
)
U
P
(
t
)
=
U
P
(
t
n
1
)
+
Q
F
n
(
t
)
C
P
,
t
=
t
D


:
{
U
P
(
t
D
)
=
U
L
(
t
n
1
)
+
2
U
D
Q
F
n
(
t
D
)
=
[
U
L
(
t
n
1
)
+
2
U
D
-
U
P
(
t
n
1
)
]
C
P
,
t
D
≤
t
≤
t
n
2


:
{
U
L
(
t
)
=
U
L
(
t
n
1
)
+
Q
F
n
(
t
)
-
Q
F
n
(
t
D
)
C
P
+
C
L
U
P
(
t
)
=
U
L
(
t
n
1
)
+
2
U
D
+
Q
F
n
(
t
)
-
Q
F
n
(
t
D
)
C
P
+
C
L
,
t
=
t
n
2


:
{
Q
L
(
t
n
2
)
=
C
L
[
Q
F
n
(
t
n
2
)
+
Q
L
(
t
n
1
)
+
Q
P
(
t
n
1
)
]
C
P
+
C
L
Q
L
(
t
n
2
)
-
2
C
L
C
P
U
D
C
P
+
C
L
Q
P
(
t
n
2
)
=
C
P
[
Q
F
n
(
t
n
2
)
+
Q
L
(
t
n
1
)
+
Q
P
(
t
n
1
)
]
C
P
+
C
L
Q
L
(
t
n
2
)
+
2
C
L
C
P
U
D
C
P
+
C
L
Q
F
n
(
t
n
2
)
=
d
31
l
h
F
n
.Roof Stage (
t
n
,
2~t
n
,
3
). During roof stage, the applied force keeps unchanged, so no charge is generated and the rectifier turns off. The leakage current of capacitors is quite low and the stage’s duration is short. It can be concluded that no charge on C
P and C
L leaks away, so equations can be given as
(16)
t
n
2
≤
t
≤
t
n
3


:
{
Q
L
(
t
)
=
Q
L
(
t
n
2
)
Q
P
(
t
)
=
Q
P
(
t
n
2
)
.
Falling Stage (
t
n
,
3~t
n
,
4
). The input force decreases to zero from F
n. This stage is similar to the rising stage and can be divided into two parts by t
D
′, the rectifier turn-on time during this stage. By referring to the analysis of rising stage, the charging process can be presented by the following equations:
(17)
t
n
3
≤
t
≤
t
D
′


:
{
U
L
(
t
)
=
U
L
(
t
n
3
)
U
P
(
t
)
=
U
P
(
t
n
3
)
+
Q
F
n
(
t
)
C
P
,
t
=
t
D
′


:
{
U
P
(
t
D
′
)
=
-
[
U
L
(
t
n
3
)
+
2
U
D
]
Q
F
n
(
t
D
′
)
=
-
[
U
L
(
t
n
3
)
+
2
U
D
+
U
P
(
t
n
3
)
]
C
P
,
t
n
3
≤
t
≤
t
n
4


:
{
U
L
(
t
)
=
U
L
(
t
n
3
)
-
Q
F
n
(
t
)
-
Q
F
n
(
t
D
′
)
C
P
+
C
L
U
P
(
t
)
=
-
[
U
L
(
t
n
3
)
+
2
U
D
-
Q
F
n
(
t
)
-
Q
F
n
(
t
D
′
)
C
P
+
C
L
]
,
t
=
t
n
4


:
{
Q
L
(
t
n
4
)
=
C
L
[
Q
L
(
t
n
3
)
-
Q
F
n
(
t
n
4
)
-
Q
P
(
t
n
3
)
]
C
P
+
C
L
Q
L
(
t
n
4
)
-
2
C
L
C
P
U
D
C
P
+
C
L
Q
P
(
t
n
4
)
=
-
C
P
[
Q
L
(
t
n
3
)
-
Q
F
n
(
t
n
4
)
-
Q
P
(
t
n
3
)
]
C
P
+
C
L
Q
L
(
t
n
4
)
-
2
C
L
C
P
U
D
C
P
+
C
L
Q
F
n
(
t
n
4
)
=
-
d
31
l
h
F
n
.Pause Stage (
t
n
,
4~t
n
+
1,1
). No force is applied on PVDF until the (
N
+
1
)th cycle begins. The charges on C
P and C
L remain constant for the same reasons presented in stay stage. Consider
(18)
t
n
4
≤
t
≤
t
(
n
+
1
)
1


:
{
Q
L
(
t
)
=
Q
L
(
t
n
4
)
Q
P
(
t
)
=
Q
P
(
t
n
4
)
.
As discussed above, the charges are generated and collected on C
L during rising and falling stages, and the applied force keeps constant with no charge generated in roof and pause stages. Recursion formulae of Q
L and Q
P can be derived and be presented by
(19)
Q
L
,
0
=
Q
L
(
N
)
,
Q
P
,
0
=
Q
P
(
N
)
,
Q
L
,
1
=
t
1
(
Q
F
+
Q
L
,
0
+
Q
P
,
0
)
-
2
t
1
U
D
C
P
,
Q
P
,
1
=
t
2
(
Q
F
+
Q
L
,
0
+
Q
P
,
0
)
+
2
t
1
U
D
C
P
,
Q
L
,
2
=
t
1
(
Q
F
+
Q
L
,
1
-
Q
P
,
1
)
-
2
t
1
U
D
C
P
,
Q
P
,
2
=
-
t
2
(
Q
F
+
Q
L
,
1
-
Q
P
,
1
)
-
2
t
2
U
D
C
L
,
t
1
=
C
L
(
C
L
+
C
P
)
,
t
2
=
C
P
(
C
L
+
C
P
)
,
Q
F
=
F
n
d
31
l
h
,
Q
L
(
N
+
1
)
=
Q
L
,
2
,
Q
P
(
N
+
1
)
=
Q
P
,
2
,
where N is the cycle number, Q
L
(
N
) is the quantity of the charges on C
L at the end of the Nth cycle, Q
P
(
N
) is the quantity of charges on C
P at the end of the Nth cycle, F
n is the force amplitude in the Nth cycle, and t
1 and t
2 are two defined parameters for simplification. Q
L
,
0 and Q
P
,
0 are the initial electric quantities of one cycle, Q
L
,
1 and Q
P
,
1 are the electric quantities when the applied force reaches F
n, and Q
L
,
2 and Q
P
,
2 are the electric quantities when the force returns to its initial value. If the force amplitude in each cycle can be acquired by test or calculation, the recursion formulae can be utilized to estimate the charges collected on storage capacitor. In some cases, the force amplitude in each cycle is a constant value F and C
P as well as C
L stores no charge initially; then (19) can be rewritten below as
(20)
Q
L
(
N
+
1
)
=
A
1
B
+
A
+
-
1
A
+
n
-
A
1
B
+
A
+
-
1
+
B
1
,
Q
P
(
N
+
1
)
=
A
2
B
+
A
+
-
1
A
+
n
-
A
2
B
+
A
+
-
1
+
B
2
,
(21)
A
1
=
(
t
1
-
t
2
)
t
1
,
A
2
=
(
t
2
-
t
1
)
t
2
,
B
1
=
-
4
t
1
2
U
D
C
P
-
2
t
1
U
D
C
P
+
2
t
1
2
Q
F
,
B
2
=
4
t
1
t
2
U
D
C
P
-
2
t
1
U
D
C
P
-
2
t
1
t
2
Q
F
,
A
+
=
A
1
+
A
2
=
(
t
1
-
t
2
)
2
,
B
+
=
B
1
+
B
2
=
2
Q
F
(
t
1
-
t
2
)
t
1
-
8
t
1
2
U
D
C
P
,
Q
F
=
F
d
31
l
h
.
It is known that C
L is much greater than C
P, and U
D is not more than 0.7 V; thus t
1 is approximately 1 and t
2 is close to 0. And the approximate expression of (20) is given below as
(22)
Q
L
(
N
+
1
)
≈
2
Q
F
+
2
Q
F
1
-
A
+
N
1
-
A
+
,
Q
P
(
N
+
1
)
≈
0
.
According to (20) and (22), it is shown that C
P and C
L are main factors that influence the charge collection, while U
D has little influence on Q
L, so it is not an important parameter that needs selecting carefully. Hence, researchers can focus on key parameters by analyzing the energy storage model. Besides, Q
L is approximately proportional to Q
F while Q
F is proportional to the applied force amplitude; therefore Q
L is approximately proportional to the force amplitude. In reality, C
L is far larger than C
P; therefore most of the generated charges are transferred to C
L. The charge increment per cycle declines as the cycle number increases, as shown in (23). The charge transfer efficiency η
Q
L is defined in (24). The energy stored on C
L is given in (25), and (26) presents the energy increment per cycle. Consider
(23)
Δ
Q
L
(
N
)
=
Q
L
(
N
)
-
Q
L
(
N
-
1
)
=
{
A
1
B
+
A
+
N
-
1
N
≠
1
B
1
≈
2
Q
F
N
=
1
,
(24)
η
Q
L
(
N
)
=
Δ
Q
L
(
N
)
2
Q
F
,
(25)
E
L
(
N
)
=
Q
L
(
N
)
2
2
C
L
,
(26)
Δ
E
L
(
N
)
=
E
L
(
N
)
-
E
L
(
N
-
1
)
.
The energy storage model is based on the energy storage circuit shown in Figure 7. The recursion formulae of stored charges are presented by (19), which can predict the charge quantity when the amplitude of applied force in each cycle is known. If the applied force keeps constant, the energy storage model can be presented by (20). Several parameters shown in (23)~(25) can evaluate the energy harvesting performance.


## 2.6. Two Approaches to Improve Energy Harvesting Performance

On the basis of the energy storage model, two approaches to improve energy harvesting performance are derived. The first approach is that the voltage on storage capacitor should to be kept low for high energy storage efficiency. Several parameters of the energy storage model are assigned to detail for convenience to clear analysis and easy demonstration. ConsiderC
L
=
112 μF, C
P
=
133 nF, Q
F
=
43.3 μC, and U
D
=
0.53 V; then Q
L, η
Q
L, E
L, and Δ
E
L are simulated in Figure 9. Apparently, η
Q
L decreases exponentially as the cycle number N increases. Δ
E
L first increases and then decreases with an increasing N. It is because that as the voltage on C
L increases, more charges are dissipated on C
P. The stored energy on C
L should be moved to another energy storage device, like a battery, or consumed by load timely before Δ
E
L begins to decrease or η
Q
L is lower than a certain value (85% for instance). After transferring the stored energy, the voltage on C
L will decrease and C
L can store energy with a high efficiency again. In conclusion, reducing the voltage across the storage capacitor is a viable approach for high energy harvesting performance.(a) The charges and energy stored on storage capacitor versus the cycle number. (b) The charges transfer efficiency and energy increment per cycle versus the cycle number.
(a)
(b)The second approach is to improve the generated chargesQ
F. The electrical connection configuration of PVDF units affects Q
F considerably, and three common configurations are studied. The first configuration is one h thick PVDF unit, and the other two configurations are n layers of h
/
n thick PVDF units that are stacked and wired in parallel and series. In the three cases, PVDF units have the same length, width, and the total thickness. The generated charges of the three cases can be given by (27) and the PVDF equivalent capacitance is expressed by (28), where S means single PVDF unit, n
P means n PVDF units in parallel, and n
S means n PVDF units in series. Obviously, thin PVDF units with a parallel connection can generate most charges, proving that connecting thin PVDF units in parallel is an effective approach to improving the energy harvesting performance. Consider
(27)
Q
F
,
S
=
d
31
l
h
F
,
Q
F
,
n
P
=
n
d
31
l
h
F
=
n
Q
F
,
S
,
Q
F
,
n
S
=
d
31
l
h
F
=
Q
F
,
S
,
(28)
C
P
,
S
=
ε
33
A
3
h
,
C
P
,
n
P
=
n
ε
33
A
3
h
/
n
=
n
2
C
P
,
S
,
C
P
,
n
S
=
1
n
ε
33
A
3
h
/
n
=
C
P
,
S
.


## 3. Experiments

Two experiments were performed to verify the models and approaches established above. The testing system for model verification is shown in Figure10. The PVDF samples consisted of several PVDF energy harvester units which were wired in series or parallel. The PVDF units were 70 mm × 20 mm in size, with sputter-deposited aluminum layer electrodes of 30 mm × 20 mm. The mechanical and electrical properties of the PVDF films made into the units were listed in Table 1. The sample capacitance C
P was measured by a digital LCR meter (Agilent E4980, at 20 Hz). The characteristics of the PVDF samples are listed in Table 2. A material testing system machine (MTS, INSTRON-8874) was run by load control to apply the forces defined in (6) on the PVDF samples. The output voltage of each sample was measured by an oscilloscope (DSOX2024A).Table 1
Properties of the PVDF materials.
Material property
Symbol
Value
Relative permittivity
ε
r
9.5 ± 1
Resistivity
ρ
r
1011 Ω·m
Mass density
ρ
1.78 g/cm3
Piezoelectric constant
d
31
17 × 10−12 C/N
Elastic modulus
Y
2500 MPa
Planar Poisson’s ratio
υ
0.3
Yield strength
σ
s
45~55 MPaTable 2
PVDF samples for verifying the models.
Sample number
C
P (nF)
Thickness (μm)
Connection configuration
1
0.127
500
Single unit
2
0.129
3
0.125
4
0.492
100
Single unit
5
0.517
6
0.510
7
0.493
100
Single unit
8
0.515
9
0.272
100 × 2
2 units in series
10
0.287
11
1.095
100 × 2
2 units in parallel
12
1.038
13
1.093
14
0.184
100 × 3
3 units in series
15
0.182
16
1.480
100 × 3
3 units in parallel
17
1.533Figure 10
Testing system.


## 4. Results and Discussions

To test the basic model and resistive load model, the experiments were performed by the MTS machine with samples 1~6. The force frequency ranged from 1 Hz to 5 Hz because of human motion frequency of about 1~5 Hz. Samples 1~3 were 500μm thick and the applied force was (
100
+
80
sin
ω
t
)
N, and samples 4~6 were 100 μm thick with the applied force of (
40
+
20
sin
ω
t
)
N. During testing, the output of a sample was terminated into the oscilloscope to measure the output voltage. The internal resistance of the oscilloscope was 10 MΩ, serving as a resistive load, and the output voltage amplitudes were simulated by (8). Figure 11 compares the experimental and simulated results of the output voltage amplitudes of samples 1~6. The six samples were tested at 1~5 Hz, and all the results of simulation were in good agreement with the ones from the experiment measurements with errors less than 5%. According to (8), (29) can be given. It shows that samples 1~3 can produce approximately the same voltage as that samples 4~6 can produce, which is also demonstrated in Figure 11. In conclusion, the resistive load model is capable of predicting the output voltages of different samples accurately, indirectly proving that the basic model is accurate. Consider
(29)
U
Amplitude
,
Sample
1
~
3
U
Amplitude
,
Sample
4
~
6
≈
C
P
,
Sample
4
~
6
C
P
,
Sample
1
~
3
·
B
Sample
1
~
3
B
Sample
4
~
6
=
4
C
P
,
Sample
4
~
6
C
P
,
Sample
1
~
3
≈
1
.
The results of PVDF samples 7~17 were presented in Figure 12 to verify the capacitive load model and energy storage model. Samples 7~17 consisted of 1~3 layers of 100 μm PVDF energy harvester units and the units were wired in series or parallel. The energy storage circuit shown in Figure 7 was used to store the generated charges on energy storage C
L of 11.3 μF. The samples were operated by the applied force (
50
+
30
sin
4
π
t
)
N. After one sample kept working for 1 min (or 120 cycles), the voltage on the storage capacitor U
L was measured to calculate the stored energy and charges. The simulated values of the voltage on C
L were calculated by (20) with Q
F predicted by (27). The samples had three connection configurations. Q
F for samples of one unit and several units in parallel was 0.153 μC, Q
F for samples of two units in series was 0.0756 μC, and Q
F for samples of three units in series was 0.051 μC. In Figure 12, the simulated and experimental results of the sample were presented, and the errors were listed. There was great consistency between the experimental data and theoretical data with errors of 3%, which demonstrated that the capacitive load model and the energy storage model were accurate enough to describe PVDF energy harvester units. Obviously, the samples of PVDF units wired in parallel generated more charges and energy than the samples of PVDF units in series did. For example, the samples of two units in parallel (sample 9 and sample 10) produced about twice as many charges as the samples of two units in series (samples 11~13) did. The results were consistent with the theoretical investigation derived from (27).Figure 11
Experimental and simulated results of the output voltage amplitude of samples 1~6.Figure 12
The experimental and simulated data ofU
L after the PVDF samples 7~17 had been operated by the applied force (
50
+
30
sin
4
π
t
)
N for 1 min.


## 5. Conclusions

This paper establishes a series of models for 31-mode PVDF energy harvester for wearable applications. Basic model presents the operation of PVDF with 31-mode. Resistive load model can test the performance of 31-mode PVDF energy harvester and predict harvested energy quantity. Capacitive load model describes the charge distribution. And energy storage model provides a way to calculate the stored charges and energy. Experiment results demonstrate the accuracies of the series of models that are better than 95% for every case tested. In order to improve and optimize the harvester, two approaches have been found through theoretical analysis and experimental investigation. Firstly, it is necessary to timely reduce the voltage across the energy storage device (a capacitor for instance). Second approach is that thin PVDF energy harvester units wired in parallel can generate more charges and improve the collected energy. The models and the approaches can serve as general rules for design and analysis of 31-mode PVDF energy harvester for wearable applications.



---

*Source: 893496-2014-07-09.xml*