Open access peer-reviewed chapter

Energy Harvesting Prediction from Piezoelectric Materials with a Dynamic System Model

By José Carlos de Carvalho Pereira

Reviewed: February 15th 2021Published: March 23rd 2021

DOI: 10.5772/intechopen.96626

Downloaded: 155


Piezoelectric vibration energy harvesting has been investigated for different applications due to the amount of wasted vibration from dynamic systems. In the case of piezoelectric materials, this energy lost to the environment can be recovered through the vibration of energy harvesting devices, which convert mechanical vibration into useful electrical energy. In this context, this chapter aims to present the mechanical/electrical coupling on a simple dynamic system model in which a linear piezoelectric material model is incorporated. For this purpose, a mechanical/electrical element of a piezoelectric disk is developed and integrated into a lumped mass, viscous damping, and spring assembling, similar to a quarter car suspension system. Equations of motion for this dynamic system in the time domain can be solved using the finite element method. The recovered electric power and energy density for PZT (Lead Zirconate Titanate) from the wasted vibration can be predicted considering that the road roughness is introduced as an input mode.


  • harvesting energy
  • wasted vibration
  • dynamic system
  • linear piezoelectric model
  • PZT

1. Introduction

Energy use is widely discussed nowadays, as energy conversion and management. In this way, new sources of energy are required to be investigated. Thus, one of the energy sources that can be used is from vibration systems, which can be subject to different excitations. This wasted energy to the environment can be recovered through vibration energy harvesting devices, which convert mechanical vibration into useful electrical energy in a way that low power devices may utilize.

Piezoelectric materials are known to have the electro/mechanical coupling effect. This property has a large range of applications in engineering. Currently, they are extensively used as sensors and actuators in vibration control systems. As a sensor, it can monitor the vibrations when bonded to a flexible structure. As an actuator, it can control the vibration level by introducing a restored force or by adding damping to the system.

In the context of recovered energy from mechanical vibrations based on the conversion of piezoelectric harvesting devices and its application on powering electronic devices, this subject has received the attention from various researchers [1, 2, 3]. An example of this type of recovered energy is the suspension system vibration for use of the vehicle itself, such as an energetic source for an active and semi-active suspension.

On a typical road, vehicles suffer accelerations due to its roughness, which excite undesired vibration. Some recently conducted reviews mentioned the potential of recovering a few hundred watts for a passenger car driven in experimental tests as well as some mathematical models [4, 5]. One of the ways to convert the mechanical energy from the vehicle suspension to electric energy is through piezoelectric materials [6]. Therefore, the objective of this chapter is to present the coupling between a piezoelectric element and a dynamic system in the context of predicting the recovered electric power and energy density for piezoelectric materials, especially the PZT (Lead Zirconate Titanate).


2. Piezoelectric material modelling

Mathematical models for predicting the harvesting energy in piezoelectric materials submitted to axial loads consider its geometric properties, diameter Dpand thickness hp, and its mechanical and electrical properties, Young’s modulus is czzE, the piezoelectric constant is ezzand the dielectric constant is ϵzzS, as shown in Figure 1. Points 1 and 2 represent the two faces of the piezoelectric disk, and wand Vare the mechanical displacement and electric potential, respectively, at these two points.

Figure 1.

The piezoelectric material model.

The electromechanical coupling effect of the piezoelectric material can be described using a set of basic equations as given in the IEEE Standard on Piezoelectricity [7]:


Where σzis the normal stress, Dzis the electric flux density, both in direction z. As seen in the above equations, the electro/mechanical coupling occurs due to the piezoelectric constant ezz.

The mechanical strain energy Umand the electric energy Ueof the piezoelectric material are written as [7]:


The electrical power can be calculated as the partial derivative of the electric energy, presented by Eq. (2), in respect to time:


3. Dynamic system modelling

This chapter aims to predict the wasted energy from vibration systems that can be further transformed into electrical energy. A typical vibration system can be described as mass, spring and damper elements, and that can represent a suspension system assembly.

The suspension system is an assembly of suspension arms or linkages, springs and shock absorbers that connect the wheels to the vehicle’s chassis in order to isolate passengers from vibrations due to bumps and roughness of the road. Furthermore, it must maintain the contact of the wheels with the road to ensure drivability. Thus, the suspension system is the mechanical system where the stability and handling of the vehicle, besides energy harvesting, must be equilibrated.

Mathematical models were initially developed for vertical vehicle performance, and the one-dimensional quarter car model is the simplest from the frequently used suspension system [8]. It is composed of the sprung mass ms, which represents ¼ of the vehicle’s body and the unsprung mass mu, which represents the wheel assembling mass. Both are considered rigid bodies. Its displacements are wsand wu, respectively and both are vertically aligned. There are some studies that include a third degree in the system to describe road roughness excitation wr. In this case, only the bounce input mode, or the vertical displacement can be implemented. Other elements of the suspension system are included, such as tire stiffness kt, suspension stiffness ks, and viscous damping cs. All vertical displacements are a function of the independent variable tthat represents the time. Figure 2 illustrates this 1D quarter car model, which could represent both the front and rear of the vehicle.

Figure 2.

1D quarter car model.

The expressions of kinetical energy from the masses, strain energy and dissipation function from the shock absorber, and the virtual work from the road roughness excitation are written as:


4. The suspension system and piezoelectric disk coupling

The conversion of the mechanical energy from the vehicle suspension to electric energy through piezoelectric materials can be predicted by piezoelectric disk coupling illustrated in Figure 1 and the suspension system illustrated in Figure 2. Since the conversion of mechanical energy to electrical energy in this system is produced by compressive efforts, the piezoelectric disk is located between the shock absorber system composed by stiffness ksand viscous damping cson the bottom side, along with the sprung mass mson the upper side, as illustrated in Figure 3. As stated previously, vertical displacements wand now the electric potential V, are all a function of the independent variable t.

Figure 3.

The suspension system and piezoelectric coupling model.


5. The vertical displacement and electric potential approach

Understanding physical problems can be accomplished when the numerical simulation of equations within the variables that describe them are represented. In the case of the electric energy prediction from the vehicle suspension system, the variables are vertical displacement wand electric potential Vof different points, or nodes, and the domain within the model is considered valid.

Some numerical methods are often used to understand physical problems, and among them, the most widely used in engineering is the finite element method. Within this method, the variables in the equations that describe the physical problem must be approximated to the polynomial functions of these variables. For the application of this method in solving problems, all elements that compose it must be represented by matrices, called elementary matrices, which are the result of the adopted polynomial functions. For further details on this method, the following references [9, 10] are recommended.

If the linear approximation function of vertical displacement wand electric potential Vover thickness hpin direction z of the piezoelectric disk are considered, stiffness elementary matrices of the piezoelectric finite element can be obtained by applying the Lagrange equations [11] over the energy expressions presented by Eq. (2).


Where Km, Km_eland Kelare mechanical elementary stiffness, electromechanical coupling elementary stiffness and electric elementary stiffness matrices, respectively, and Spis the cross-sectional area of the piezoelectric disk.

Moreover, if the same linear approximation function of the vertical displacement wover stiffness ksand viscous damping csof the shock absorber in direction z are considered, the suspension system’s differential equation of motion, illustrated in Figure 2, can be obtained by applying the Lagrange equations [11] over the energy expressions presented by Eq. (4).


Where the elementary matrices and the vectors are:


Using the technique of assembling the elementary matrices of the finite element method, the suspension system’s differential equation of motion and piezoelectric disk coupled model, as illustrated in Figure 3, is as shown in Eq. (8):


Eq. (9) can be manipulated and substituted for Eq. (8). Thus, the final equation of motion as a function of only mechanical variables is:


Eq. (10) can be solved in the time domain with an integration method, such as the Newmark Method [12]. The mechanical force is due to road roughness as well as the tire characteristics of the wheel as shown in Eq. (7). These data are used to apply the condition in each time step in solving Eq. (10), in which all variables ws=w2, w1, and wuand their time derivatives are obtained.

The response in the time domain in respect to vertical displacements is obtained by using Eq. (10) and subsequently, the electric potential is obtained as:


The electrical energy can be calculated as presented by Eq. (2). The development of this equation follows:


Thus, the expression of the electrical energy due to the piezoelectric disk on this suspension system is as:


In addition, the electrical power can be calculated as presented by Eq. (3) and its development is as follows:


Thus, the expression of the electrical power due to the piezoelectric disk on this suspension system is as:


Where ẇ1, ẇ2, V̇1and V̇2are time derivatives of the displacements and electric potential of nodes 1and 2of the piezoelectric disk.


6. Application example

For a simpler demonstration of the potential for harvesting energy in a suspension system by means of piezoelectric material, a MATLAB® code to obtain the results was developed and presented below.

The dimensions of the piezoelectric disk are diameter Dp=0.065mand thickness hp=0.025m. The vehicle data and the properties of piezoelectric material PZT-5H are shown in Tables 1 and 2, as indicated in references [13, 14], respectively.

Property descriptionValue
Body and wheel mass
¼ Vehicle sprung mass (mu)362.5 kg
Unsprung mass (ms)39 kg
Tire stiffness (kt)200 kN/m
Shock absorber suspension
Suspension stiffness (ks)30 kN/m
Suspension damping (cs)4 kN s/m

Table 1.

Vehicle data.

Piezoelectric constant – ezz[C/m2]23.30
Dielectric constant – ϵzzS[F/m]1.30 x 10−8
Density – ρ [kg/m3]7500
Young’s modulus – czzE[GPa]23.0

Table 2.

Properties of piezoelectric material.

The road excitation is as depicted in [6]. The excitation is the bounce input mode vertical position of the road wrtwith class D road (poor) and a driving speed for a vehicle of 20 m/s.

In Figure 4, the instantaneous electric energy density for piezoelectric material PZT-5H calculated by Eq. (13), with a simulation time of 30 s and a fixed step time of 0.00001 s.

Figure 4.

Electric energy density response for PZT-5H.

Figure 5 exhibits the overall instantaneous electric power density response for piezoelectric material PZT-5H calculated by Eq. (15).

Figure 5.

Electric power density response for PZT-5H.


7. Conclusions

This chapter proposes a coupled suspension system and piezoelectric model to predict the potential of harvested electric power in vehicle suspension systems. The performance of piezoelectric material PZT-5H was investigated, in respect to harvesting energy based on energy density and electric power density.

The approach presented in this chapter is a way to simulate the electric power generated in vehicle suspension systems established by piezoelectric harvesting. Nonetheless, these results would need to be compared with experimental results to demonstrate the validity of the proposed model.

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

José Carlos de Carvalho Pereira (March 23rd 2021). Energy Harvesting Prediction from Piezoelectric Materials with a Dynamic System Model, Piezoelectric Actuators - Principles, Design, Experiments and Applications, Hu Huang and Jianping Li, IntechOpen, DOI: 10.5772/intechopen.96626. Available from:

chapter statistics

155total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Design, Characterisation and Prospect of Piezoelectric Microfluidic Technology

By Song Chen, Zhonghua Zhang, Junwu Kan, Jianping Li and Jianming Wen

Related Book

First chapter

Introductory Chapter: Properties and Processing of Metallic Glasses

By Hu Huang

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us