Vibration Energy Harvesting from Plates by Means of Piezoelectric Dynamic Vibration Absorbers

Preprint

Article

Vibration Energy Harvesting from Plates by Means of Piezoelectric Dynamic Vibration Absorbers

Altmetrics

Downloads

Views

Comments

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

15 December 2023

Posted:

24 December 2023

You are already at the latest version

Alerts

Abstract

In this paper the possibility of harvesting energy from the vibrations of a plate is analyzed. The harvester takes the form of a cantilever dynamic vibration absorber equipped with a piezoelectric layer and tuned by means of a tip mass to the first mode of vibration of the plate. A mathematical model of the coupled system composed of the plate and the harvester is presented. The validity of the proposed harvester is proved by means of simulations carried out with the modal expansions approach. Simulation results highlighting the effect of harvester tuning and location are presented as well. Then the validity of the harvester is confirmed by experimental tests carried out both with a concentrated impulsive load and with a distributed pressure load. Simulations and experimental tests show that the cantilever piezoelectric dynamic vibration absorber generates open circuit voltage an order of magnitude larger than the one generated by the same piezoelectric layer when directly bonded to the plate surface.

Keywords:

Subject: Engineering - Mechanical Engineering

1. Introduction

The present era, which is characterized by IoT and by Industry 4.0 requires a lot of sensors to monitor the industrial, urban and natural environments. The use of wirings to feed the sensors is sometimes unfeasible, e.g. this is the case of remote sensors installed in buoys in deep sea or in tall buildings and antennas. Nowadays remote sensors nodes are fed by means of batteries, which are not eco-friendly and require complex and expensive replacement operations. Fortunately, the energy consumption of sensor nodes is constantly decreasing, and new technical solutions are possible. Vibration energy harvesting is a great opportunity, since remote sensors can be fed by using existing sources of energy. Vibration energy harvesting can be performed exploiting different physical phenomena to convert mechanical energy into electrical energy. There are piezoelectric [1], electromagnetic [2], electrostatic [3] and triboelectric [4] harvesters. Piezoelectric harvesters usually take the form of thin layers that either are bonded to a vibrating surface, or are bonded to a cantilever beam to make a cantilever harvester. Piezoelectric harvesters generate high voltage, can be integrated with the vibrating structure and have small mass and encumbrance. For these reasons they are adopted to feed sensors and other small electronic equipment in vehicles, industrial machinery, civil structures and biomedical devices.

The vibrating plate is one of the most common vibrating structures and it is a typical part of machines, domestic appliances, vehicles and aircrafts [5,6]. Vibrating plates can be found in buildings (panels of facades) and in other civil structures as well [7,8]. Hence, vibration energy harvesting from plates is a very present research topic. Plates are sometimes excited by concentrated loads, but there are many cases in which this structure is excited by the distributed loads caused by wind, traffic or by water flows. The simplest piezoelectric harvester for a vibrating plate is a piezoelectric layer directly bonded to the plate surface. It works in the bending mode (the 3-1 mode [9]) and exploits the strain caused by plate bending. The strain of the plate reaches the maximum value at the natural frequencies of the modes of vibrations that are excited by concentrated or distributed loads. A possible alternative, in order to better exploit the resonance phenomenon, is a cantilever harvester with tip mass tuned to the natural frequency of the most exited mode of vibration of the plate. This cantilever harvester behaves as a dynamic vibration absorber (DVA) that is able to receive a large amount of energy from the plate. On the one hand plate vibrations are strongly reduced [10], on the other hand the cantilever harvester experiences large amplitude vibrations.

DVA with piezoelectric layers have been studied in recent years both with the aim of suppressing the vibrations of the host structure and with the aim of scavenging energy. Ali and Adhikari in 2013 [11] proposed an energy harvesting DVA and analytically studied the coupling of this device with a 1 degree of freedom (DOF) vibrating system. A lumped parameter 1 DOF model of the piezoelectric harvester was adopted. An analytical expression of the natural frequencies of the coupled system was given, the presence of two resonance peaks was demostrated and the influence of the piezoelectric layer on the damping of the DVA was highlighted. Abdelmoula ed alii in 2017 [12] studied the coupling of a cantilever harvester with a 1 DOF vibrating system. The cantilever harvester was modeled according to the modal expansion approach like in [13]. A strategy was proposed to minimize the vibrations of the host structure and to maximize the harvested power. Experimental results dealing with the effect of a cantilever piezoelectric dynamic absorber mounted on a vibrating rigid beam were presented by Sulaiman et alii in 2021 [14]. Rezaei et alii in 2022 [15] extended the previous analyses considering the coupling between a flexible beam and a cantilever harvester. Both the cantilever harvester and the flexible beam were modeled adopting the modal expansion approach. The number of modes of vibration of the host structure and of the cantilever harvester needed to accurately model the system’s dynamics was investigated. Results showed that considering only the first mode of the host structure and of the cantilever harvester is enough to capture the main features of the system. In 2012 Rajarathinam and Ali [16] improved the analysis of a 2 DOF system composed of a host vibrating structure and an energy harvesting DVA considering both random excitation and system uncertainties. Results showed that uncertainties in the natural frequency lead to a decrease in the harvested power. Recently, Rezaei et alii [17] presented experimental results dealing with a flexible beam equipped with a cantilever harvester and excited by a shaker. They developed a detailed integrated mode taking into account shaker dynamics. A good agreement between experimental and analytical results was obtained.

The effect of a piezoelectric DVA on a vibrating plate was not addressed in the above mentioned researches. Hence the first aim of this paper is to study vibration energy harvesting from a vibrating plate by means of cantilever piezoelectric dynamic vibration absorber (CPDVA). The second aim of this paper is to study the performance of a CPDVA installed on a plate exited by distributed loads that are very common in actual structures. The third aim of this paper is to make a comparison between the performances of a simple piezoelectric layer and a CPDVA mounted on the same plate.

The paper is organized as follows. In section 2 the mathematical model of a vibrating plate with a CPDVA is developed, plate dynamics are modeled with the modal expansion approach, whereas the harvester is modeled with a lumped parameter approach that takes into account the first mode of vibration. In section 3 some calculated results are presented showing the effect of harvester tuning and positioning on the generated voltage. The experimental setup is described in section 4. A prototype CPDVA was built using a commercial cantilever harvester M2814P2C-01 made by Smart Material GmbH. This device was installed in a clamped-clamped-free-free plate excited by concentrated and distributed loads. Section 5 deals with experimental results and focuses on the comparison between the voltage generated by the CPDVA and the one generate by a simple piezoelectric layer directly bonded to the plate. Finally, conclusions are drawn in section 6.

2. Mathematical Model

In this study a rectangular aluminium plate clamped on two opposite sides is considered. The system and the main dimensions are represented in Figure 1. The plate lays on the

x - y

plane, and vibrations occur in the z direction. The cantilever piezoelectric dynamic vibration absorber (CPDVA) consists of a cantilever beam with its axis parallel to the plate and fixed to the plate through one of its extremities. A tip mass is attached to the free end of the cantilever for tuning. Concentrated and distributed loads in the z direction are considered, since they excite the vibrations of the plate in the transversal direction. The vibrations are then transferred to the CPDVA which harvests energy through the piezoelectric patch attached to the cantilever.

2.1. Plate model

The equation of the forced vibrations of a plate is:

ρ \ddot{w} (x, y, t) + C_{p l a t e} \dot{w} (x, y, t) + D (\frac{\partial^{4} w (x, y, t)}{\partial x^{4}} + 2 \frac{\partial^{4} w (x, y, t)}{\partial x^{2} \partial y^{2}} + \frac{\partial^{4} w (x, y, t)}{\partial y^{4}}) = F (x, y, t)

(1)

where w is the transversal displacement of any point of the plate, which is a function of time and space,

F (x, y, t)

is a force distribution,

ρ

is the surface density of the plate, D is the flexural rigidity of the plate, and

C_{p l a t e}

is a damping coefficient; a proportional damping is assumed. The plate has homogeneous thickness h, and the flexural rigidity can be calculated from the modulus of elasticity E and from the Poisson’s ratio

ν

D = E h^{3} / (12 (1 - ν^{2}))

Vibration modes of the plate can be found through the associated homogeneous equation of (1). According to the Rayleigh method, when the i-th mode of the plate is considered, which has m nodal lines in the x direction and n nodal lines in the y direction, the transversal displacement of each point of the rectangular plate can be expressed by the function:

Φ_{i} (x, y) \equiv Φ_{m n} (x, y) = X_{m} (x) Y_{n} (y)

(2)

X_{m} (x)

and

Y_{n} (y)

are the fundamental mode shapes along the x and y direction respectively and depend on the boundary conditions of the plate. For the rectangular plate clamped at

x = 0

and

x = a

the

X_{m} (x)

, functions are:

X_{m} (x) = \{\begin{matrix} cos (γ_{1} \cdot (\frac{x}{a} - \frac{1}{2})) + \frac{sin (γ_{1} / 2)}{sinh (γ_{1} / 2)} cosh (γ_{1} \cdot (\frac{x}{a} - \frac{1}{2})), & m = 2, 4, 6, \dots \\ sin (γ_{2} \cdot (\frac{x}{a} - \frac{1}{2})) - \frac{sin (γ_{2} / 2)}{sinh (γ_{2} / 2)} sinh (γ_{2} \cdot (\frac{x}{a} - \frac{1}{2})), & m = 3, 5, 7, \dots \end{matrix}

(3)

where

γ_{1}

and

γ_{2}

are respectively the roots of:

tan (\frac{γ_{1}}{2}) + tanh (\frac{γ_{1}}{2}) = 0

(4)

tan (\frac{γ_{2}}{2}) - tanh (\frac{γ_{2}}{2}) = 0

(5)

The functions

Y_{n} (y)

are defined as:

Y_{n} (y) = \{\begin{matrix} 1, & n = 0 \\ 1 - \frac{2 y}{b}, & n = 1 \\ cos (γ_{1} \cdot (\frac{y}{b} - \frac{1}{2})) - \frac{sin (γ_{1} / 2)}{sinh (γ_{1} / 2)} cosh (γ_{1} \cdot (\frac{y}{b} - \frac{1}{2})), & n = 2, 4, 6, \dots \\ sin (γ_{2} \cdot (\frac{y}{b} - \frac{1}{2})) + \frac{sin (γ_{2} / 2)}{sinh (γ_{2} / 2)} sinh (γ_{2} \cdot (\frac{y}{b} - \frac{1}{2})), & n = 3, 5, 7, \dots \end{matrix}

(6)

The following orthogonality conditions hold for plate modes:

\int_{0}^{a} \int_{0}^{b} ρ Φ_{j} (x, y) Φ_{i} (x, y) d x d y = m_{i} δ_{i j}

(7)

\int_{0}^{a} \int_{0}^{b} Φ_{j} (x, y) \nabla^{4} Φ_{i} (x, y) d x d y = ω_{i}^{2} δ_{i j}

(8)

where

δ_{i j}

is equal to 1 if

i = j

and 0 otherwise;

m_{i}

and

ω_{i}

are the modal mass and the natural angular frequency of the mode i respectively; and the operator

\nabla^{4}

is defined as

\nabla^{4} = \frac{\partial^{4}}{\partial x^{4}} + 2 \frac{\partial^{4}}{\partial x^{2} \partial y^{2}} + \frac{\partial^{4}}{\partial y^{4}}

The parameters of the plate considered in this study are:

a = 1 m

b = 0.6 m

h = 3 \cdot 10^{- 3} m

E = 6.9 \cdot 10^{10} P a

, and

ν = 0.33

. From the plate parameters, natural frequency and shape of each mode can be determined (for the calculation of natural frequencies refer to [18]). In this paper, only the first six modes of the plate will be used to describe plate dynamics; the natural frequencies and the corresponding m-n combinations are reported in Table 1, and Figure 2 shows the modal shapes.

2.2. Harvester lumped element model

The CPDVA consists of a cantilever beam partially covered by a piezoelectric layer. The cantilever is fixed to the plate at coordinates (

x_{a}, y_{a}

), and a tip mass is fixed on the cantilever tip. Figure 3 represents the harvester and its geometrical parameters.

The harvester is modeled through a single-degree-of-freedom (SDOF) approach, the variable that describes the configuration of the harvester is the tip displacement with respect to the harvester base (

u_{r} (t)

). The equation of motion of the harvester subject to a base excitation is:

M {\ddot{u}}_{r} (t) + C {\dot{u}}_{r} (t) + K u_{r} (t) = - α M \ddot{w} (x_{a}, y_{a}, t)

(9)

where M, C, and K are the lumped mass, damping, and stiffness of the harvester respectively;

α

is a corrective factor for the base excitation in SDOF models ([1,9]).

The lumped parameters are calculated following the approach used and validated in [19]. First of all, the displacement of the cantilever along the beam axis is expressed as the product between a shape function

g (ξ)

and the relative tip displacement:

u (ξ, t) = g (ξ) u_{r} (t)

(10)

where

ξ

is the coordinate along the beam axis,

ξ

is zero at the clamp location and it increases moving to the tip. The static deformed shape of a cantilever with a concentrated load on its tip is assumed to approximate the shape of the first mode of the harvester. Therefore,

g (ξ)

is defined as:

g (ξ) = \frac{3 ξ^{2} L - ξ^{3}}{2 L^{3}}

(11)

where L is the length of the cantilever. Moreover, with this assumption, the ratio between the tip rotation (

φ

) and the tip displacement is:

\frac{φ}{u_{r}} = \frac{3}{2 L}

(12)

Using the Lagrangian approach, the following expressions for the lumped parameters of the harvester are found:

M = \frac{33}{140} (m_{s} + m_{p}) + M_{t} {(1 + \frac{3 w_{t}}{4 L})}^{2} + \frac{M_{t}}{8 L} (h_{t}^{2} + w_{t}^{2})

(13)

K = \frac{3 E_{h} I_{h}}{L^{3}} + \frac{θ^{2}}{C_{p}}

(14)

C = 2 ζ \sqrt{K M}

(15)

where

m_{s}

m_{p}

are the masses of the structural and of the piezoelectric layer of the cantilever respectively;

M_{t}

w_{t}

h_{t}

are the mass, the width, and the height of the tip mass;

E_{h} I_{h}

is the equivalent bending stiffness of the beam;

C_{p}

is the capacitance of the piezoelectric layer;

θ

is the electro-mechanical coupling coefficient; and

ζ

is the damping coefficient of the harvester. It is worth noting that stiffness K depends on the backward piezoelectric effect as well.

In open circuit conditions, the electro-mechanical coupling coefficient is defined by the following equation ([13]):

θ = - e_{3, 1} b_{p} h_{p c} \int_{ξ_{1}}^{ξ_{2}} \frac{d^{2} g (ξ)}{d ξ^{2}} d ξ

(16)

where

e_{3, 1}

is the piezoelectric constant of the material;

b_{p}

is the width of the piezoelectric layer;

h_{p c}

is the distance between the neutral axis of the composite beam and the center of the piezoelectric layer; and

ξ_{1}

ξ_{2}

are the coordinates that identify the extremities of the piezoelectric layer along the beam axis, see Figure 3.

Finally, the open circuit voltage produced by the harvester is expressed as ([19]):

v_{o c} (t) = \frac{θ}{C_{p}} u_{r} (t)

(17)

2.3. Coupled equations

The aim of this section is the calculation of frequency response functions (FRFs) between applied loads, generated voltage, and plate displacement, since they make possible the calculation of the response of the system in the presence of harmonic, periodic, transient and random loads. To generalize the analysis the following loads acting on the plate are considered: a uniform distribuited load

q_{z} (t)

, a concentrated force

F_{c} (t)

applied in the point of coordinates (

x_{c}, y_{c}

), and the forcing term due to the CPDVA

F_{a} (t)

. The uniform distribuited load can represent both the action of a fluid in contact with the plate and an inertial load resulting from the acceleration of the two clamped sides of the plate (

{\ddot{z}}_{f} (t)

) due to environmental vibrations. In this case:

q_{z} (t) = - ρ {\ddot{z}}_{f} (t)

(18)

The concentrated force is a reference load condition that simulates the effect of impulsive tests carried out by means of hammers for modal testing.

The force due to the CPDVA is the sum of two terms. The first term is the inertial force related to the oscillation of the harvester base which has mass

M_{a}

and is rigidly fixed to the plate. The second term is related to the oscillation of the cantilever with tip mass. The second term can be obtained by isolating the damping and the stiffness terms in Equation (9). Therefore, the following expression is found:

F_{a} (t) = - M ({\ddot{u}}_{r} (t) + (α + \frac{M_{a}}{M}) \ddot{w} (x_{a}, y_{a}, t))

(19)

In order to obtain the equations of motion of the system composed by the plate and the CPDVA, the distributed and the concentrated forces are inserted in (1), and the resulting equation is combined with (9):

\{\begin{matrix} ρ \ddot{w} (x, y, t) + C_{p l a t e} \dot{w} (x, y, t) + D \nabla^{4} w (x, y, t) = q_{z} (t) + F_{c} (t) δ (x - x_{c}) δ (y - y_{c}) + F_{a} (t) δ (x - x_{a}) δ (y - y_{a}) \\ M {\ddot{u}}_{r} (t) + C {\dot{u}}_{r} (t) + K u_{r} (t) = - α M \ddot{w} (x_{a}, y_{a}, t) \end{matrix}

(20)

where

δ

is the Dirac delta function.

The transversal displacement of the plate is expressed as a linear combination of the modes of the plate:

w (x, y, t) = \sum_{i = 1}^{\infty} Φ_{i} (x, y) η_{i} (t)

(21)

in which

η_{i} (t)

is the modal coordinate of mode i. In practical simulations, the summation in (21) is truncated taking into account only the first N modes of the plate in the frequency band of interest. Substituting (19) and (21) in (20) gives:

\{\begin{matrix} ρ \sum_{i = 1}^{N} Φ_{i} (x, y) {\ddot{η}}_{i} (t) + C_{p l a t e} \sum_{i = 1}^{N} Φ_{i} (x, y) {\dot{η}}_{i} (t) + D \nabla^{4} Φ_{i} (x, y) η_{i} (t) = q_{z} (t) + F_{c} (t) δ (x - x_{c}) δ (y - y_{c}) - \\ M ({\ddot{u}}_{r} (t) + (α + \frac{M_{a}}{M}) \sum_{i = 1}^{N} Φ_{i} (x_{a}, y_{a}) {\ddot{η}}_{i} (t)) δ (x - x_{a}) δ (y - y_{a}) \\ M {\ddot{u}}_{r} (t) + C {\dot{u}}_{r} (t) + K u_{r} (t) = - α M \sum_{i = 1}^{N} Φ_{i} (x_{a}, y_{a}) {\ddot{η}}_{i} (t) \end{matrix}

(22)

Then, the orthogonality conditions (7) and (8) are exploited by multiplying (22) for the generic mode

Φ_{i} (x, y)

and integrating on the whole plate. In this way, the following system of

N + 1

differential equations in modal coordinates is obtained:

\{\begin{matrix} m_{i} {\ddot{η}}_{i} (t) + 2 m_{i} ζ_{i} ω_{i} {\dot{η}}_{i} (t) + m_{i} ω_{i}^{2} η_{i} (t) = \int_{0}^{a} \int_{0}^{b} Φ_{i} (x, y) q_{z} (t) d x d y + Φ_{i} (x_{c}, y_{c}) F_{c} (t) - \\ M Φ_{i} (x_{a}, y_{a}) ({\ddot{u}}_{r} (t) + (α + \frac{M_{a}}{M}) \sum_{i = 1}^{N} Φ_{i} (x_{a}, y_{a}) {\ddot{η}}_{i} (t)) \\ M {\ddot{u}}_{r} (t) + C {\dot{u}}_{r} (t) + K u_{r} (t) = - α M \sum_{i = 1}^{N} Φ_{i} (x_{a}, y_{a}) {\ddot{η}}_{i} (t) \end{matrix}

(23)

where the first equation holds for

i = 1, 2, \dots, N

For the computation of the frequency response functions (FRFs), the following harmonic solutions are assumed:

η_{i} (t) = η_{i, 0} \cdot e^{i ω t}

(24)

u_{r} (t) = u_{r, 0} \cdot e^{i ω t}

(25)

q_{z} (t) = q_{z, 0} \cdot e^{i ω t}

(26)

F_{c} (t) = F_{c, 0} \cdot e^{i ω t}

(27)

where i is the imaginary unit.

Introducing (24),(25), (26), and(27) in (23), this system of linear equations holds:

\{\begin{matrix} - ω^{2} m_{i} η_{i, 0} + 2 i ω m_{i} ζ_{i} ω_{i} η_{i, 0} + m_{i} ω_{i}^{2} η_{i, 0} (t) = \int_{0}^{a} \int_{0}^{b} Φ_{i} (x, y) q_{z, 0} d x d y + Φ_{i} (x_{c}, y_{c}) F_{c, 0} + \\ ω^{2} M Φ_{i} (x_{a}, y_{a}) (u_{r, 0} + (α + \frac{M_{a}}{M}) \sum_{i = 1}^{N} Φ_{i} (x_{a}, y_{a}) η_{i, 0}) \\ - ω^{2} M u_{r, 0} + i ω C u_{r, 0} + K u_{r, 0} = ω^{2} α M \sum_{i = 1}^{N} Φ_{i} (x_{a}, y_{a}) η_{i, 0} \end{matrix}

(28)

where the first equation holds for

i = 1, 2, \dots, N

System (28) is rearranged letting the unknown

η_{1, 0}

η_{2, 0}

, ...,

η_{N, 0}

, and

u_{r, 0}

at the left side and moving the assigned forcing terms to the right side. Moreover, for brevity, the following symbols are adopted:

\begin{matrix} Φ_{i}^{a} = Φ_{i} (x_{a}, y_{a}) \\ Φ_{i}^{c} = Φ_{i} (x_{c}, y_{c}) \\ A_{i} = - m_{i} ω^{2} + 2 i m_{i} ζ_{i} ω_{i} ω + m_{i} ω_{i}^{2} \\ B_{i} = M ω^{2} Φ_{i}^{a} (α + \frac{M_{a}}{M}) \\ P_{i} = \int_{0}^{a} \int_{0}^{b} Φ_{i} (x, y) d x d y \\ D = - M ω^{2} + i C ω + K \end{matrix}

(29)

Therefore:

\{\begin{matrix} (A_{i} - Φ_{i}^{a} B_{i}) η_{i, 0} - Φ_{i}^{a} \sum_{j = 1}^{n} B_{j} η_{j, 0} (1 - δ_{i j}) - B_{i} \frac{M}{α M + M_{a}} u_{r, 0} = P_{i} q_{z, 0} + Φ_{i}^{c} F_{c, 0} \\ - α M ω^{2} \sum_{i = 1}^{n} Φ_{i}^{a} η_{i, 0} + D u_{r, 0} = 0 \end{matrix}

(30)

where the first equation holds for

i = 1, 2, \dots, N

Considering the uniform distributed load amplitude

q_{z, 0}

or the concentrated force amplitude

F_{c, 0}

as inputs, Equation (30) can be solved for a range of forcing frequencies

ω

. The solution can be computed in Matlab using standard routines, and amplitudes

η_{1, 0}

η_{2, 0}

, ...,

η_{N, 0}

, and

u_{r, 0}

can be obtained at each frequency. The amplitude of the plate displacement

w_{0}

can be calculated from the modal coordinates by means of:

w_{0} (x, y, ω) = \sum_{i = 1}^{N} Φ_{i} (x, y) η_{i, 0} (ω)

(31)

Therefore FRFs between the plate displacement (and acceleration) in the generic point

x_{P}, y_{P}

and the loads can be calculated using the following definitions:

F R F_{w_{P}, q_{z}} (ω) = \frac{w_{0} (x_{P}, y_{P}, ω)}{q_{z, 0}} F R F_{{\ddot{w}}_{P}, q_{z}} (ω) = - ω^{2} \frac{w_{0} (x_{P}, y_{P}, ω)}{q_{z, 0}}

(32)

F R F_{w_{P}, F_{c}} (ω) = \frac{w_{0} (x_{P}, y_{P}, ω)}{F_{c, 0}} F R F_{{\ddot{w}}_{P}, F_{c}} (ω) = - ω^{2} \frac{w_{0} (x_{P}, y_{P}, ω)}{F_{c, 0}}

(33)

where (32) are used for the uniform distributed load, and (33) are used for the concentrated load.

The amplitude of the open circuit (OC) voltage can be retrieved exploiting Equation (17); the following expression holds:

v_{o c, 0} = \frac{θ}{C_{p}} u_{r, 0}

(34)

Finally, the FRFs between OC voltage and applied loads are:

F R F_{v, q_{z}} (ω) = \frac{v_{o c, 0} (ω)}{q_{z, 0}}

(35)

F R F_{v, F_{c}} (ω) = \frac{v_{o c, 0} (ω)}{F_{c, 0}}

(36)

for the case with distributed and concentrated load respectively.

3. Calculated results

The first mode of vibration of a plate has the largest probability of being stimulated by the environmental sources of excitation (wind, traffic, rain) for two reasons, the former related to frequency, the latter related to the modal shape.

Typically, environmental sources are characterized by power spectral densities (PSDs) showing the maximum amplitudes in the low frequency range. Turbulence of grazing flows used to collect energy form artificial piezoelectric grass [20] shows relevant PSD amplitudes below

50 H z

. Wind turbulence PSDs used for the design of civil structures show the maximum amplitude below

1 H z

[21,22,23]. The PSDs of traffic induced vibrations show their main peaks below

20 H z

[24].

The first mode of vibration of a plate clamped on two opposite sides has the simplest modal shape (Figure 2). Higher order modes show an increasing number of nodal lines which define different zones of the plate that vibrate in phase opposition. A spatially uniform pressure fluctuation on the surface of the plate can excite the first mode of vibration but is not able to transfer much energy to higher order modes with different areas that vibrate in phase opposition. To highlight this concept the FRF between the acceleration of the point at the center of the plate and the distributed load

q_{z} (t)

was calculated by means of the mathematical model of Section 2. The FRF magnitude, which is plotted in Figure 4, shows that only the 1st and the 6th mode of vibration are excited by the distributed load. The other modes are not excited at all, because there is a complete cancellation between the pressure loads applied in the areas of the plate the move in opposition. Owing to the dependence on frequency squared (Equation (32)) the acceleration FRF shows the highest peak for the 6th mode. If the sizes of the plate have the order of magnitude of

1 m

, they are much smaller than the typical length scale of wind turbulence that excites buildings [23], therefore pressure can be considered spatially uniform. For pressure fluctuations related to turbulence having smaller length scales the possible excitation of higher order modes has to be evaluated using more sophisticated methods like the acceptance integral [25], nevertheless, also in this case the first mode is the most excited [25]. In many practical conditions the panels can be excited by vibrations coming from the supports [26], also in this case the lower order modes are the most excited.

For the above-mentioned reasons this research focuses on vibration energy harvesting from the first mode of vibration of the plate and the CPDVA is tuned to the natural frequency of this mode and is located at the center of the plate, which is an anti-node for the first mode of vibration. The mathematical model of the CPDVA is developed starting from the properties of a commercial cantilever harvester M2814P2C-01 built by Smart Material. This harvester has a piezoelectric layer made of Macro Fiber Composite [27] and has a structural layer made of glass fiber composite (FR4). Its properties are summarized in Table 2. The transfer function of the M2814P2C-01 was experimentally evaluated with the impulsive method [28] and showed a resonance peak at

146.69 H z

. Therefore, the addition of a large tip mass is needed to tune the harvester to the first natural frequency of the plate (

16.5 H z

A small clamp is needed to fix the cantilever base to the plate. If the clamp is made of polymeric material, its mass is about

0.018 k g

. Before tuning the harvester to the plate the effect of the mass of the clamp on the natural frequency of the first plate mode was evaluated by means of the mathematical model of Section 2 setting

M = 0 k g

. The clamp mass has a very small effect on the natural frequency of the first mode that moves from

16.5 H z

16.4 H z

Harvester tuning is then carried out by means of Equation (9) setting the term at the right hand side to zero (free vibrations analysis). The undamped natural frequency is given by:

ω_{n} = \sqrt{K / M}

(37)

In which K and M are given by equations (14) and (13) respectively. Since all the terms in the expressions of K and M are known except

M_{t}

, the undamped natural frequency (

ω_{n}

) is set equal to the natural frequency of the first mode of the plate, and the tip mass that tunes the harvester is found: the calculated value is

M_{t} = 0.051 k g

The calculations reported in this section aim to show the validity of the CPDVA and the importance of tuning and location of the device. The focus is on voltage generation. Figure 5 makes a comparison between the acceleration FRF at the center of the bare plate and the acceleration FRF of the same point when the CPDVA is tuned to the first mode of the plate and is mounted at the center of the plate. To calculate the FRFs a distributed load with constant amplitude and variable frequency (

q_{z} (t)

) is assumed to excite the plate.

As in the other cases considered in the scientific literature (e.g. beams [15]) the CPDVA eliminates the original resonance peak of the plate and generates a couple of new peaks: the first at lower frequency, the second at higher frequency. The amplitudes of the new acceleration peaks are significantly lower than the one of the original resonance peak of the bare plate. Figure 6 makes a comparison between the OC voltage FRF of the CPDVA tuned to the first mode of the plate and mounted at the center of the plate and the OC voltage FRF of a piezoelectric patch directly bonded under the center of the plate. The first mode of the plate corresponds to the bending of a beam in the x-z plane, hence, when the first mode is excited, the curvature at the center of the plate is large and the electro-mechanical coupling between the patch and the plate is large. Indeed, for the 1st mode, (16) can be generalized for the patch bonded on the plate as:

θ = - e_{3, 1} b_{p} h_{p c p} \int_{x_{1}}^{x_{2}} \frac{\partial^{2} Φ_{1}}{\partial x^{2}} d x

(38)

where

h_{p c p}

is the distance between the center of the piezoelectric layer and the neutral axis of the plate, and

x_{1}

x_{2}

are the coordinates of the extremities of the patch. This piezoelectric patch has the same properties and sizes of the piezoelectric layer included in the M2814P2C-01 harvester. Figure 6 highlights that the CPDVA generates much more voltage than the patch, because not only the two peaks of CPDVA are higher than the peak of the patch, but also the bandwidth around resonances are larger. The latter effect can be very useful in the presence of broad-band random excitation.

In the previous simulations the CPDVA was mounted at the center of the plate (

x_{a} = a / 2 y_{a} = b / 2

), where the first mode of vibration shows the largest amplitudes (see Figure 2). In the simulations reported in Figure 7 the CPDVA is moved from the best location along the longitudinal direction (x). It is worth noting that, since the first mode of the plate is symmetric, leftwards and rightwards displacements have the same effect on the coupling between the CPDVA and this mode of the plate. An harmonic distributed load with constant amplitude and variable frequency (

q_{z} (t)

) is assumed to excite the plate.

The FRF between acceleration of the point of the plate just above the the CPDVA and the distributed load shows that the coupling between the CPDVA and the plate decreases when the CPDVA is moved from the the center of the plate i.e. when

x_{a} < a / 2

. Indeed the peaks of the FRF become closer and their amplitudes decrease. Eventually, when the CPDVA is located near the nodal line the two peaks tend to merge into a unique peak.

The FRF between generated OC voltage and distributed load shows that the bandwidth including the two resonance peaks decreases when the CPDVA is moved from the central position. This is a negative effect when the plate is excited by broad-band random pressure fluctuations. The effect of CPDVA position on the peaks of the OC voltage FRF is more complex. If the lateral displacement is small (e.g.

x_{a} = a / 4

) the amplitudes of the two peaks decrease. The maximum amplitude is reached when the two peaks merge into a unique peak (

x_{a} = a / 10

). Eventually, peak amplitude decreases when the CPDVA is close to the nodal line of the plate mode. This behavior can be exploited if the plate is excited by a harmonic load having a well defined frequency.

In Figure 8 the CPDVA is moved from the center of the plate in the lateral direction (y). This displacement does not affect the coupling with the first mode of vibration of the plate and with the other modes with

n = 0

(see 1 and Figure 2), but it alters the coupling with the modes with

n \neq 0

, e.g. the second mode. Hence, this effect can be found only when the CPDVA is coupled with a plate. The FRF between the acceleration of the plate point above the CPDVA and the distributed load shows small modifications in the position and height of the two main resonance peaks, but, when the CPDVA is moved towards the anti-node of the second mode (

y_{a} = 0

), a minor resonance peak at the frequency of the second peak appears. This phenomenon happens because when

y_{a}

decreases there is a coupling between the CPDVA and the second mode as well. Since this coupling alters the shape of the second mode, the uniform distributed load is able to cause a small excitation of the second mode. The excitation of the second mode slightly alters the effect of the CPDVA on the first mode and the OC voltage FRF shows small variations in the position and height of the main resonance peaks. For

x_{a} = a / 2, y_{a} = 0

the second peak becomes the higher.

Figure 9 deals with the effect of CPDVA tuning, in this case the harvester is always at the center of the plate, which is excited by a harmonic distributed load (

q_{z} (t)

The maximum effect on plate acceleration takes place when the CPDVA is tuned to the natural frequency of the first mode of the plate. This behavior is typical of many vibrating systems [15]. A shift in the tuning frequency

Δ f = 3 H z

(about 20%) is enough to make appear a large resonance peak again. The frequency of this resonance is close to the one of the original resonance of the plate and depends on CPDVA tuning.

Figure 9 shows that all the OC voltage FRFs intersect at the frequency of the original mode of the plate (

16.5 H z

). The exact tuning leads to the most symmetric behavior with two peaks of the OC voltage having nearly the same amplitude and the minimum at the natural frequency of the original plate mode.

If the tuning frequency is shifted from the value corresponding to the natural frequency of the first mode, the maximum peak of voltage FRF always coincides with the maximum peak of the acceleration FRF. When the tuning frequency decreases the first peak moves to lower frequencies and decreases its amplitude, whereas the second peak moves towards the natural frequency of the plate, increases and eventually decreases when

Δ f = - 6 H z

. In a simply supported beam forced by a concentrated load and equipped with a cantilever DVA a similar behavior was found, when the value of the tuning mass was increased (decreased tuning frequency) [15]. When the tuning frequency increases the second peak moves to higher frequencies, whereas the first moves towards the original natural frequency of the plate, increases and eventually decreases when

Δ f = - 6 H z

. The behavior is not symmetric, since positive and negative frequency shifts have different effects. The behavior depicted in Figure 9 suggests that, when the distributed load has a broad-band spectrum, the best tuning coincides with the exact tuning to the natural frequency of the plate mode. When a narrow-band distributed load excites the plate, the largest amount of energy can be harvested when the harvester is tuned to a slightly higher frequency.

4. Experimental equipment

The experimental equipment was developed with the aim of evaluating the performance of the harvesting system in different load conditions. To make the clamps, the two shorter sides of the aluminium plate were fixed to a very stiff steel base by means of screws (Figure 10). The CPDVA used in experiments was built using a commercial cantilever harvester (PEH M2814P2C-01) made by Smart Materials GmbH. The base to attach the harvester to the plate was 3D printed. Table 2 reports the parameters of the harvester and Figure 11 shows the CPDVA attached to the backside of the plate. The CPDVA base was bonded to the center of the plate (

x = a / 2

y = b / 2

) by means of structural adhesive. In order to verify the advantages of the CPDVA, experiments were also performed with a piezoelectric patch in place of the CPDVA. In particular, the same piezoelectric patch used in the PEH M2814P2C-01 cantilever harvester was bonded under the center of the plate as shown in Figure 12.

During tests the OC voltage generated by the piezoelectric layer and the acceleration of the center of the plate were measured. The latter measurement was carried out using a small size PCB 352C22 accelerometer (sensitivity

10.48 m V / g

). Since the mass of this accelerometer is less than 1 gram, it does not affect the dynamics of the plate.

Two load cases were considered in the framework of this research. In the first the center of the plate was excited by means of an hammer for modal testing (PCB 086C03). Since in this case the measurement of the exciting force was possible and reliable, the FRFs of the system were calculated. It is worth noting that many simulations showed that there are small differences between the FRFs obtained with uniform pressure load and the FRFs obtained with concentrated force acting at the center of the plate.

In the second load case a fan was installed over the plate at a distance of

1.2 m

, see Figure 10. The plate was excited by the pressure generated by the air jet impinging on the plate surface. The fan diameter is

0.148 m

and the maximum flow rate is

595 m^{3} / h

, which corresponds to a Reynolds number

R e \approx 10^{5}

. Experimental tests dealing with turbulent jets impinging on flat surfaces were presented in [29,30,31]. Results show axial-symmetric distributions of both mean velocities and fluctuating components, which reach the maximum values at small distances from jet axis, that is at the center of the plate. The pressure fluctuations caused by the abovementioned pressure distributions can excite the first mode of the plate, but cannot excite the modes with odd number of nodal lines.

The voltage generated by the CPDVA and the signals of the sensors were acquired through a NI9230 board, typically, 4096 samples were collected with a sampling frequency of

2048 H z

. Then, digital signals were analyzed with the software NI Signal Express to obtain the FRFs between voltage, acceleration, and force and PSDs of OC voltage.

5. Experimental results

5.1. Tests with impulsive excitation

A preliminary experimental modal analysis of the clamped plate was carried out to check the natural frequencies and the modal shapes. A good agreement with analytical results was found. Since the experimental natural frequency of the first mode resulted slightly lower than than the analytical value, the tip mas was increased (

M_{t} = 0.058 k g

) .

For the load case with impulsive excitation, the plate equipped with the CPDVA and the accelerometer was excited through the hammer for modal testing, the test was repeated for seven times. Figure 13 shows the impulsive force applied by the hammer and the measured acceleration and voltage. As shown by experimental data, applying this type of excitation gives a good repeatability of the results. The presence of beats in the signal of the voltage generated by the CPDVA highlights the presence of two close natural frequencies in the system.

The measured data was elaborated to retrieve the FRF between the acceleration at the center of the plate and the concentrated load applied in the center of the plate and the FRF between the OC voltage and the concentrated load applied in the center of the plate. The magnitudes of these FRFs are shown in Figure 14 together with the magnitudes of the FRFs predicted by the mathematical model.

As expected, the original resonance peak of the plate (

15.4 H z

) is eliminated by the CPDVA and a couple of new peaks is generated; the first at lower frequency (

14.0 H z

), the second at higher frequency (

16.9 H z

). The presence of the two peaks in the voltage FRF enables energy harvesting in a broad frequency band; this is useful in the case of random excitation. Curves of different tests are very similar. There is a good agreement between the experimental FRFs and the FRFs predicted by the mathematical model with concentrated load excitation (Section 2). In particular, the model predict with good accuracy the frequency and the height of the FRF peaks, which represent the largest contributes to voltage generation.

In order to evaluate the performance of the CPDVA, the tests with hammer excitation were repeated substituting the CPDVA with the piezoelectric patch bonded to the plate. Measured acceleration and voltage generated by the simple patch are shown in Figure 15, and the amplitudes of the FRFs obtained from these measures are shown in Figure 16.

The analysis of time domain data shows that the maximum voltage generated by the patch is about 3 times smaller than the maximum voltage generated by the CPDVA and very quickly reaches values lower than

\pm 1 V

. The experimental FRFs show the highest peaks at the first natural frequency of the plate. The height of the peak of the voltage FRF obtained with the patch is one order of magnitude smaller than the peaks obtained with the CPDVA.

In Figure 16, a small peak corresponding to the second natural frequency of the plate can be observed. Indeed, the patch was bonded with an offset in the y-direction with respect to the center of the plate. Moreover, the hammer hit was manually generated and it was not exactly centered and aligned with the perpendicular of the plate. Thus, it was able to weakly excite the second mode of the plate.

5.2. Tests with air excitation

The fan depicted in Figure 10 was used to excite the plate by means of random pressure oscillations. Several tests were carried out with air velocity ranging from

4.1

9.0 m / s

. In the first series of tests the piezoelectric patch was bonded under the center of the plate, whereas in the second series of tests the CPDVA was fixed under the center of the plate. The OC voltage generated by the piezoelectric layer was analyzed in the frequency domain by means of the NI Signal Express to obtain the power spectral density (PSD). In every test condition the test was repeated seven times, a good repeatability was found. Figure 17 refers to the voltage generated by the patch. The voltage PSD shows a peak at the natural frequency of the first mode of the plate. The height of this peak strongly increases with air velocity and reaches the maximum value (

0.02 V^{2} / H z

) when air velocity is

9.0 m / s

. It is worth noting that at the maximum air velocity a minor peak at the natural frequency of the second mode of the plate appears as well. Since the air flow is able to strongly excite only the first mode of vibration of the plate, pressure fluctuations generated by the air flow can be considered spatially uniform on the plate surface and dominated by low frequency components, as happens in the excitation of buildings by wind.

Figure 18 shows the PSD of OC voltage generated by the CPDVA. Also in this case the generated voltage strongly increases with air velocity, but amplitudes are much larger than in the previous case with peaks ranging form

0.5

1 V^{2} / H z

when air velocity reaches the maximum value. The effect of the CPDVA is clearly shown by the presence of the two peaks around the frequency of the first mode of the plate, like in the impulsive tests, the first peak is the higher.

6. Conclusions

In many conditions relevant to industrial and civil applications the response of a vibrating plate is dominated by the first mode of vibration, having the lowest natural frequency and the simplest modal shape. Direct energy harvesting from plate vibrations by means of a piezoelectric patch bonded to the plate is not very efficient, since the plate is stiff an the strain is small. The CPDVA transfers the vibration energy of the plate to the small cantilever of the harvester which has a less stiff section (made of piezoelectric and plastic materials), therefore the generation of voltage is strongly enhanced. This result is confirmed both by simulations and by experimental tests carried out with impulsive excitation and air excitation. Simulations showed that the harvester’s performance is affected by location and tuning, nevertheless the system is rather robust and small variations in tuning and location do not have a dramatic effect on generated voltage. In the presence of narrow-band excitation a shift of the tuning frequency can improve performance. Future studies will deal with CPDVA optimization in the presence of narrow- and broad-band excitations.

Author Contributions

Conceptualization, M.T., A.P., and A.D.; methodology, M.T., A.P., and A.D.; software, M.T. and A.P.; validation, M.T., A.P., and A.D.; formal analysis, M.T., A.P., and A.D.; investigation, M.T., A.P., and A.D.; writing—original draft preparation, M.T., A.P., and A.D.; writing—review and editing, M.T., A.P., and A.D.; supervision, A.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

Erturk, A.; Inman, D.J. On mechanical modeling of cantilevered piezoelectric vibration energy harvesters. Journal of intelligent material systems and structures 2008, 19, 1311–1325. [Google Scholar] [CrossRef]
Xie, Z.; Teng, L.; Wang, H.; Liu, Y.; Fu, M.; Liang, J. A Self-Powered Synchronous Switch Energy Extraction Circuit for Electromagnetic Energy Harvesting Enhancement. IEEE Transactions on Power Electronics 2023. [Google Scholar] [CrossRef]
Chen, L.; Ma, Y.; Hou, C.; Su, X.; Li, H. Modeling and analysis of dual modules cantilever-based electrostatic energy harvester with stoppers. Applied Mathematical Modelling 2023, 116, 350–371. [Google Scholar] [CrossRef]
Ye, C.; Liu, D.; Chen, P.; Cao, L.N.; Li, X.; Jiang, T.; Wang, Z.L. An Integrated Solar Panel with a Triboelectric Nanogenerator Array for Synergistic Harvesting of Raindrop and Solar Energy. Advanced Materials 2023, 35, 2209713. [Google Scholar] [CrossRef]
Warburton, G. The vibration of rectangular plates. Proceedings of the Institution of Mechanical Engineers 1954, 168, 371–384. [Google Scholar] [CrossRef]
Tommasino, D.; Moro, F.; Zumalde, E.; Kunzmann, J.; Doria, A. An Analytical–Numerical Method for Simulating the Performance of Piezoelectric Harvesters Mounted on Wing Slats. Actuators. MDPI, 2023, Vol. 12, p. 29. [CrossRef]
Bedon, C.; Fasan, M.; Amadio, C. Vibration Analysis and Dynamic Characterization of Structural Glass Elements with Different Restraints Based on Operational Modal Analysis. Buildings 2019, 9. [Google Scholar] [CrossRef]
Pipitone, G.; Barone, G.; Palmeri, A. Optimal design of double-skin façades as vibration absorbers. Structural Control and Health Monitoring 2018, 25. [Google Scholar] [CrossRef]
Priya, S.; Inman, D.J. Energy harvesting technologies; Springer, 2009; Volume 21. [Google Scholar]
Thorby, D. Structural Dynamics and Vibration in Practice; Elsevier Ltd., 2008. [Google Scholar]
Ali, S.F.; Adhikari, S. Energy harvesting dynamic vibration absorbers. Journal of Applied Mechanics 2013, 80, 041004. [Google Scholar] [CrossRef]
Abdelmoula, H.; Dai, H.; Abdelkefi, A.; Wang, L. Control of base-excited dynamical systems through piezoelectric energy harvesting absorber. Smart Materials and Structures 2017, 26, 095013. [Google Scholar] [CrossRef]
Erturk, A.; Inman, D.J. An experimentally validated bimorph cantilever model for piezoelectric energy harvesting from base excitations. Smart materials and structures 2009, 18, 025009. [Google Scholar] [CrossRef]
Sulaiman, L.H.; Anuar, M.A.; Zamri, A. An empirical study on the effectiveness of energy harvesting from dynamic vibration absorber. ARPN Journal of Engineering and Applied Sciences 2021, 16, 2674–2683. [Google Scholar]
Rezaei, M.; Talebitooti, R.; Liao, W.H. Concurrent energy harvesting and vibration suppression utilizing PZT-based dynamic vibration absorber. Archive of Applied Mechanics 2022, 1–20. [Google Scholar] [CrossRef]
Rajarathinam, M.; Ali, S.F. Parametric uncertainty and random excitation in energy harvesting dynamic vibration absorber. ASCE-ASME Journal of Risk and Uncertainty in Engineering Systems, Part B: Mechanical Engineering 2021, 7, 010905. [Google Scholar] [CrossRef]
Rezaei, M.; Talebitooti, R.; Liao, W.H.; Friswell, M.I. Integrating PZT layer with tuned mass damper for simultaneous vibration suppression and energy harvesting considering exciter dynamics: An analytical and experimental study. Journal of Sound and Vibration 2023, 546, 117413. [Google Scholar] [CrossRef]
Leissa, A.W. Vibration of plates; Scientific and Technical Information Division, National Aeronautics, 1969; Volume 160. [Google Scholar]
Tommasino, D.; Moro, F.; de Pablo Corona, E.; Vandi, L.; Baietta, A.; Pracucci, A.; Doria, A. Optimization of a Piezoelectric Wind-Excited Cantilever for Energy Harvesting from Facades. Advances in Italian Mechanism Science; Niola, V., Gasparetto, A., Quaglia, G., Carbone, G., Eds.; Springer International Publishing: Cham, 2022; pp. 848–856. [Google Scholar] [CrossRef]
Hobeck, J.D.; Inman, D.J. Artificial piezoelectric grass for energy harvesting from turbulence-induced vibration. Smart Materials and Structures 2012, 21, 105024. [Google Scholar] [CrossRef]
Eurocode 1: Actions on structures - Part 1-4: General actions - Wind actions. Standard, European committee for standardization, Brussels, BE, 2010.
Tao, T.; Wang, H.; Wu, T. Comparative Study of the Wind Characteristics of a Strong Wind Event Based on Stationary and Nonstationary Models. Journal of Structural Engineering 2017, 143, 04016230. [Google Scholar] [CrossRef]
Roncallo, L.; Gimondo, M.; Tubino, F. Dynamic Response of Slender Vertical Structures Subjected to Thunderstorm Outflows. Applied Sciences 2023, 13. [Google Scholar] [CrossRef]
Hao, H.; Ang, T.C. Analytical Modeling of Traffic-Induced Ground Vibrations. Journal of Engineering Mechanics 1998, 124, 921–928. [Google Scholar] [CrossRef]
Au-Yang, M.K. Joint and Cross Acceptances for Cross-Flow-Induced Vibration—Part I: Theoretical and Finite Element Formulations. Journal of Pressure Vessel Technology 2000, 122, 349–354. [Google Scholar] [CrossRef]
Tommasino, D.; Moro, F.; Bernay, B.; De Lumley Woodyear, T.; de Pablo Corona, E.; Doria, A. Vibration Energy Harvesting by Means of Piezoelectric Patches: Application to Aircrafts. Sensors 2022, 22. [Google Scholar] [CrossRef]
Kuang, Y.; Zhu, M. Evaluation and validation of equivalent properties of macro fibre composites for piezoelectric transducer modelling. Composites Part B: Engineering 2019, 158, 189–197. [Google Scholar] [CrossRef]
Tommasino, D.; Tonan, M.; Moro, F.; Doria, A. Identification of the Piezoelectric Properties of Materials From Impulsive Tests on Cantilever Harvesters. International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, 2023, Vol. Volume 12: 35th Conference on Mechanical Vibration and Sound (VIB). International Design Engineering Technical Conferences and Computers and Information in Engineering Conference V012T12A005. [CrossRef]
Tummers, M.J.; Jacobse, J.; Voorbrood, S.G. Turbulent flow in the near field of a round impinging jet. International Journal of Heat and Mass Transfer 2011, 54, 4939–4948. [Google Scholar] [CrossRef]
Hassan, M.E.; Assoum, H.; Sobolìk, V.; Vétel, J.; Abed-Meraïm, K.; Garon, A.; Sakout, A. Experimental investigation of the wall shear stress and the vortex dynamics in a circular impinging jet. Experiments in Fluids 2012. [Google Scholar] [CrossRef]
Shiyong, Y.; Yong, G.; Nan, J.; Junjie, L. An experimental study of a turbulent jet impinging on a flat surface. International Journal of Heat and Mass Transfer 2015, 83, 820–832. [Google Scholar] [CrossRef]

Figure 1. Energy harvesting system (top). Front and lateral views of the plate (bottom).

Figure 2. Calculated modes of vibration of the plate.

Figure 3. Scheme of the CPDVA (left), and harvester cross section (right).

Figure 4. The FRF (magnitude) between the acceleration at the center of the plate and the uniform distributed load.

Figure 5. Comparison between the magnitude of calculated FRF of the acceleration at the center of the plate and the uniform distributed load. Bare plate with patch and plate with CPDVA.

Figure 6. Comparison between the magnitude of calculated FRF of the generated OC voltage and the uniform distributed load. Bare plate with patch and plate with CPDVA.

Figure 7. FRF (magnitude) of the acceleration at the center of the plate (top) and FRF of the OC voltage generated by the CPDVA (bottom). Comparison considering different positions of the harvester along x-axis (and the same position along y-axis) and the uniform distributed load.

Figure 8. FRF (magnitude) of the acceleration at the center of the plate (top) and FRF of the OC voltage generated by the CPDVA (bottom). Comparison considering different positions of the harvester along y-axis (and the same position along x-axis) and the uniform distributed load.

Figure 9. FRF (magnitude) of the acceleration at the center of the plate (top) and FRF of the OC voltage generated by the CPDVA (bottom). Comparison considering different tuning frequencies of the harvester and the uniform distributed load.

Figure 10. Clamped plate and fan for the excitation of the plate through the air jet.

Figure 11. Picture of the CPDVA. The tip mass is on the left hand side and the base glued to the plate is on the right hand side.

Figure 12. Piezoelectric patch directly bonded to the plate.

Figure 13. Experimental results obtained by means of an impulse applied at the center of the plate witch CPDVA: impulse (left), acceleration at the center of the plate (center) and OC voltage generated by the CPDVA (right).

Figure 14. Magnitude of the FRF between the acceleration at the center of the plate and the impulse (left) and magnitude of the FRF between the the OC voltage generated by the CPDVA and the impulse (right). Experimental and predicted FRF for the impulse applied at the center of the plate.

Figure 15. Experimental results obtained by means of an impulse applied at the center of the plate with the piezoelectric patch directly bonded to the plate: impulse (left), acceleration at the center of the plate (center) and OC voltage generated by the patch (right).

Figure 16. Magnitude of the FRF between the acceleration at the center of the plate and the impulse (left) and the magnitude of the FRF between the OC voltage generated by the piezoelectric patch bonded to the plate and the impulse (right). Experimental and predicted FRF for the impulse applied at the center of the plate.

Figure 17. PSD of the OC voltage generated by the patch directly bonded on the plate excited by the air jet. Experimental results are reported for three values of air velocity.

Figure 18. PSD of the OC voltage generated by the CPDVA excited by the air jet. Experimental results are reported for three values of air velocity.

Table 1. Calculated natural frequency and m-n combination for the first six plate modes.

mode number (i)	m	n	frequency [Hz]
1	2	0	16.49
2	2	1	23.89
3	3	0	45.45
4	3	1	56.40
5	2	2	61.27
6	4	0	89.08

Table 2. Parameters of the M2814P2C-01 cantilever harvester.

Structural layer length	$70 m m$
Structural layer thickness	$1 m m$
Structural layer width	$25 m m$
Structural layer Elastic modulus	$24 G P a$
Structural layer density	$1920 k g m^{- 3}$
Piezoelectric patch length	$28 m m$
Piezoelectric patch thickness	$0.18 m m$
Piezoelectric patch width	$14 m m$
Piezoelectric patch Elastic modulus	$30.24 G P a$
Piezoelectric patch density	$7400 k g m^{- 3}$
Piezoelectric patch capacitance	$30.8 n F$
Piezoelectric constant $e_{31}$	$5.02 p C m^{-} 2$

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Vibration Energy Harvesting from Plates by Means of Piezoelectric Dynamic Vibration Absorbers

Abstract

1. Introduction

2. Mathematical Model

2.1. Plate model

2.2. Harvester lumped element model

2.3. Coupled equations

3. Calculated results

4. Experimental equipment

5. Experimental results

5.1. Tests with impulsive excitation

5.2. Tests with air excitation

6. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Conflicts of Interest

References

MDPI Initiatives

Important Links

Subscribe