Next Article in Journal
Machine Learning Algorithm for Efficient Design of Separated Buffer Super-Junction IGBT
Next Article in Special Issue
Research on Output Characteristics of Microscale BST Laminate Structure Based on Mixed Finite Element Method
Previous Article in Journal
Multiband Signal Receiver by Using an Optical Bandpass Filter Integrated with a Photodetector on a Chip
Previous Article in Special Issue
A Pavement Piezoelectric Energy Harvester for Small Input Displacements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Objective Topology Optimization of a Broadband Piezoelectric Energy Harvester

1
Department of Engineering, Jade University of Applied Sciences, Friedrich-Paffrath-Str. 101, 26389 Wilhelmshaven, Germany
2
Institute for Electronic Appliances and Circuits, Faculty of Computer Science and Electrical Engineering, University of Rostock, Albert-Einstein-Str. 2, 18059 Rostock, Germany
3
Institute of Microstructure Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(2), 332; https://doi.org/10.3390/mi14020332
Submission received: 16 December 2022 / Revised: 20 January 2023 / Accepted: 20 January 2023 / Published: 27 January 2023

Abstract

:
In recent years, topology optimization has proved itself to be state of the art in the design of mechanical structures. At the same time, energy harvesting has gained a lot of attention in research and industry. In this work, we present a novel topology optimization of a multi-resonant piezoelectric energy-harvester device. The goal is to develop a broadband design that can generate constant power output over a range of frequencies, thus enabling reliable operation under changing environmental conditions. To achieve this goal, topology optimization is implemented with a combined-objective function, which tackles both the frequency requirement and the power-output characteristic. The optimization suggests a promising design, with satisfactory frequency characteristics.

1. Introduction

Energy harvesting devices convert free ambient energy into usable electrical energy. This energy can then be used as an additional power source for wireless autonomous devices, e.g., inside a wireless sensor network. These devices are often deployed under harsh environmental conditions, causing them to be not easily accessible [1]. There are different sources of ambient energy available for energy harvesting, such as solar, wind, and mechanical energy. In our works, we focus on energy harvesting from mechanical vibrations that occur due to the motion of rotatory machines, e.g., in an industrial setting or in automotive. Different techniques are applicable for vibrational energy harvesting: electromagnetic, electrostatic, magnetostrictive, and piezoelectric [2]. Of all these, piezoelectric energy harvesters (PEH) are the most common type of vibration energy harvesters, due to their high achievable power density, scalability, and simplistic geometry [3,4,5,6]. The simplest design consists of a flexible beam structure that carries piezoelectric material, which transforms vibration energy into usable electric power. Typically, a proof mass is added, in order to achieve the desired resonant frequency for a given application. However, very often the ”harvesting frequency” changes over time, due to, e.g., changes in environmental conditions (temperature and humidity) or due to the aging of the vibration source. In such case, a simple PEH fails, as the structure represents a single mass–spring–damper system and, therefore, only exhibits a single dominant resonant mode.
An adaptation or tuning of the PEH’s resonant frequency is considered one of many broadband energy-harvesting techniques applied to avoid a mismatch of excitation frequency and the PEH’s resonant frequency. Comprehensive reviews on the topic of broadband piezoelectric energy harvesting can be found in [7,8]. In this work, we consider a so-called multi-modal broadband energy-harvesting approach. Multi-modal energy harvesting requires a device that exhibits multiple resonant frequencies that, therefore, can cover a larger frequency bandwidth. In [9], a bi-directional, U-shaped structure with cross-connected beams was introduced that can capture vibration in 3D space. In [10,11], ”star”- and a ‘’pizza”-shaped arrays were proposed and, in [12], a trapezoidal geometry and optimal segmentation for higher modes were presented. Most of these recent designs are based on combining multiple simple beam PEHs with different resonant frequencies into an array, which is often referred to as an ”in parallel” configuration. In our case, we focus on serial designs, i.e., cascading simple beam PEHs instead of arranging them into an array [13]. In [14], we showed that cascading/serial designs significantly outperform array configurations and confirmed previously presented results from, e.g., [15,16].
We refer to our compact dual-frequency resonator as the folded-beam design. Its design goal is to have the first two resonant frequencies with similar power amplitudes in close proximity. Therefore, it can provide a broad usable frequency range for energy harvesting. This design approach is our key differentiation from other authors. A major challenge of our design is to account for unwanted interactions between the two modes, such as switching or merging. That is why manual design approaches become unfeasible and mathematical optimization approaches are to be favored.
The authors of [17] discuss various optimization algorithms with regard to their respective suitability for finding the suitable geometric dimensions of PEHs. More device-specific optimizations are presented in [18], where multiple topologies of a tip-excited single beam PEH are parametrically optimized over a power spectrum from 0 Hz to 800 Hz. In [19], the parameter optimization of a magnetic coupled PEH was performed, optimizing the spacing of two coupled resonators for different magnetic flux densities. In [20], stack PEH for automotive is optimized, with regard to parameters, such as the number of stack layers or height to cross-section ratio. For the folded-beam PEH, we implemented a two-stage global parametric design optimization and proposed novel geometries in [14,21,22]. Parameter optimizations, however, are strongly limited, as the topology of the final design is already determined by the initial design. For example, the final design cannot be star-shaped if the initial design was rectangular. Hence, in this work, we explore topology optimization (TO) to gain more design freedom.
The methodology of TO was introduced in [23] and has become very popular across all fields of engineering. The TO methodology is designed to find an optimal layout of a structure given the available design space, applied loads, support conditions, and constraints such as maximum weight. This is achieved by reformulating the classical optimization problem with a set of distributed functions on a fixed domain [24]. In most applications, this reformulation results in an optimal material distribution problem on a finite-element grid. The most common TO application is weight or compliance minimization. However, recent works show that TO is also considered for thermal applications [25], where a cold plate for a lithium-ion battery was optimized, or for optical applications [26], where the desired band gap characteristics for phononic crystals were obtained.
The first studies of TO on the eigenfrequency can be found in [27], where the support of a plate-like structure was optimized. Later, in [28,29], the modified optimality criteria and the mean frequency goal function were introduced, which significantly improved the performance of TO for dynamical eigenfrequency problems. In [30], we introduced a dynamic compliance-based optimization approach for resonant structures that uses the compliance goal function to have more control over the modal shape and for a more stable optimization process. As for applications, in [31], a 2D in-plane single mass MEMS gyroscopes were optimized using TO, which was later extended to a tuning fork resonator, i.e., a two mass configuration, in [32]. In [33], a composite material PEH was proposed as a result of a multi-material TO. More recently, the same authors proposed plate-type PEH designs [34], maximizing open-circuit voltage by topology optimizing the through-thickness of piezoelectric material distribution.
In this work, we present the TO of a multi-modal PEH and propose a novel folded-beam design that fulfills its unique frequency requirements. The TO procedure implements a multi-objective goal function to tackle both the frequency and the power amplitude requirements. The achieved optimal design is characterized in a simulation. The manuscript is structured as follows: the introduction of the folded-beam PEH is provided in Section 2; in Section 3, we introduce our TO approach and establish the optimization setup; and, in Section 4 we present the optimization results and discuss them in Section 5. Section 6 concludes this work and provides an outlook to future research.

2. Definition and Modeling of Multi-Modal Piezoelectric Energy Harvester

Piezoelectric materials can convert mechanical vibration into electrical energy using the piezoelectric effect. The corresponding governing equation is usually provided in coupled strain-charge form [35]:
S = s T + σ T E ,
D = σ T + ϵ E ,
where S R 6 × 1 and T R 6 × 1 are the stress and strain tensors, respectively, and D is the electric displacement vector. σ R 3 × 6 s R 6 × 6 denotes the compliance tensor measured with a constant electric field E and ϵ R 3 × 3 denotes the permittivity tensor at constant strain. σ is the piezoelectric tensor. As most piezoelectric materials are classified as ceramics, they are extremely brittle and require an elastic support structure for this designed use case. Therefore, a PEH is typically composed of an elastic mechanical resonator that carries piezoelectric material, which is attached to the areas that are exposed to mechanical stress or strain. The design of a PEH mostly focuses on the geometry of the mechanical resonator in order to achieve desired characteristics.

2.1. Finite Element Modeling of Piezoelectric Energy Harvester

The finite element discretization of the governing equation results in a second-order system of dynamic algebraic equations:
Σ = M m 0 0 0 = : M x ¨ m x ¨ el + D x ˙ m x ˙ el + K m K m , el K m , el T K el = : K x m x el = b 0 u y = C x m x el ,
where x m are the nodal displacements and x el the nodal electrical potentials. M m and K m denote the mechanical mass and stiffness matrices, respectively. K el is the dielectric matrix and K m , el describes the piezoelectric coupling. The input u acts on the system via the input vector b and C gathers the user-defined outputs of the system y . Finally, D = α M + β K defines the Rayleigh damping for the system.

2.2. Folded-Beam Resonator

The folded-beam design was introduced in [13]. Its geometry equals two cascading simple beam harvesters, with the second beam folded toward the fixed end of the first beam instead of further extending it. The initially proposed design is shown in Figure 1.
The design can be fabricated from a single sheet of steel. It includes piezoelectric films on both cantilever beams and permanent magnets at their free ends. These magnets can be used for bi-directional frequency tuning via external magnetic forces, which is, however, out of the scope of the current work. The folded-beam geometry is designed to have two fundamental modes, corresponding to the two respective beams, at specific frequencies in the range of from 50 Hz to 100 Hz.

3. Topology Optimization Procedure

In general, a TO problem is defined as follows:
min ρ g ( ρ )
subject to:
Σ : M ( ρ ) x ¨ + D ( ρ ) x ˙ + K ( ρ ) x = b u y = C x
e = 1 n ρ e V e V * 0 , ( ρ = [ ρ e ] e [ 1 , n ] , V * = α V 0 ) ,
0 ρ e 1 ,
where the design parameter ρ is a vector containing all the elemental pseudo-densities ρ e [ 0 , 1 ] . g is the objective function and Σ is the dynamical system resulting from finite element discretization introduced in (3). The (6) introduces a volume constraint for the optimization. The density-dependent system matrices are provided by the classical simple isotropic material model with penalization (SIMP model):
M ( ρ ) = e = 1 n ρ e M e , K ( ρ ) = e = 1 n ρ e p K e ,
with M e and K e , respectively, being the elemental mass and stiffness matrices of a full finite element.

3.1. Objective Function and Sensitivity Analysis

The main design goal of the folded-beam PEH is to have the two fundamental modes inside a pre-defined small frequency range. In this way, their response peaks overlap and allow for efficient energy harvesting over the entire frequency range (cf. Figure 2). However, the modes are also required to generate similar levels of power output, in order to realize this usable frequency range, as if one of the modes has a significantly higher power output, the second mode may not even be observable. These two goals are taken into account via a combined-objective function that consists of two parts. The first part g 1 is the mean frequency objective function first introduced in [29], which has already been involved in our previous works [36]. It is defined as:
g 1 : = 1 w 1 + w 2 w 1 ( λ 1 λ 1 , goal ) 2 + w 2 ( λ 2 λ 2 , goal ) 2 1 2 ,
where λ 1 , λ 2 and λ 1 , goal , λ 2 , goal are the first two eigenvalues of the generalized eigenvalue problem corresponding to model (3) and their respective target values. Note that each of the eigenfrequencies of the strucure f i are connected to λ i as follows: λ i = ( 2 π f i ) 2 . The weighting factors w 1 and w 2 are chosen as 1 / λ 1 , goal and 1 / λ 2 , goal , respectively. This nondimensionlizes the expression and ensures an equal convergence for all eigenvalues [29]. The second part of the objective function g 2 represents the relative difference of the electrical energy of the structure at both resonant modes p e l , 1 and p e l , 2 , i.e.:
g 2 : = 2 ( p e l , 1 p e l , 2 ) p e l , 1 + p e l , 2 2 .
The dielectric energy of the structure [37,38] can be computed as:
p e l , i = 1 2 ϕ i , el T K el ϕ i , el ,
where ϕ i , el contains the electrical degrees of freedom of the i-th eigenvector. The partial objective functions are combined in a weighted sum:
g = v 1 g 1 + v 2 g 2 ,
and, analogously, the sensitivity equals the weighted sum of both partial sensitivities:
g ρ e = v 1 g 1 ρ e + v 2 g 2 ρ e .
The weighting factors v 1 and v 2 are chosen such that both sensitivities are in the same order of magnitude. The partial sensitivities can be computed as:
g 1 ρ e = 1 ( w 1 + w 2 ) g 1 w 1 ( λ 1 λ 1 , goal ) λ 1 ρ e + w 2 ( λ 2 λ 2 , goal ) λ 2 ρ e
for the first part, where λ i ρ e is obtained by (cf. Appendix A):
λ i ρ e = ϕ i , e T K e ρ e λ i M e ρ e ϕ i , e .
The second part of the sensitivity is obtained using the quotient rule:
g 2 ρ e = 8 g 2 p e l , 1 ρ e p e l , 2 ρ e ( p e l , 1 + p e l , 2 ) ( p e l , 1 p e l , 2 ) p e l , 1 ρ e + p e l , 2 ρ e ( p e l , 1 + p e l , 2 ) 2 ,
where p e l , i ρ e are, again, computed via the adjoint method (cf. Appendix B):
p e l , i ρ e = μ i , m T K m ρ e λ i M m ρ e ϕ i , m ,
and μ i , m is the mechanical part of the adjoint state vector that solves the adjoint problem:
K m λ i M m K m , el K m , el T K el μ i , m μ i , el = 0 1 2 K el ϕ i , el .

3.2. Filters

Having elemental pseudo-density as its design variable, TO suffers from a number of issues. Firstly, the solution depends on the discretization of the design domain. Finer meshes often generate thin details that cannot be obtained by coarser meshes and would therefore converge to a completely different topology. Secondly, TO on a coarse mesh often converges to designs suffering from the well-known checkerboard problem [39,40]. These problems can be tackled by applying appropriate filter techniques [24,41], these filters themselves, however, introduce a third issue. The aforementioned filtering techniques regularize elemental density based on a weighted average of the elemental density of adjacent elements and, therefore, introduce blurred, i.e., gray areas with intermediate elemental densities, which are not manufacturable. We have experienced this issue extensively in our previous works [30,36]. In order to resolve this issue and to obtain manufacturable designs, we have implemented a Heaviside filter utilizing the threshold projection suggested in [42]. Figure 3 shows both filters and their effects on the same optimization process.
To summarize, we have implemented both the density filter and the threshold projection in our optimization to achieve a manufacturable design. First, we apply a density filter on the candidate design encoded in ρ , which is suggested by the optimizer. Subsequently, the density-filtered design ρ ˜ is projected onto a pre-defined threshold using the Heaviside filter, providing us with the final design ρ ˜ ^ . Both filters will be introduced in detail in the following subsections.

3.2.1. Density Filter

Let ρ i be the elemental pseudo-density of the i-th element computed by the optimizer. The density filter replaces this value with ρ ˜ i by performing a weighted average of its adjacent elemental densities:
ρ ˜ i = j N e , i w ( r j ) ρ j j N e , i w ( r j ) , w ( r j ) = R min | r i r j | ,
where N e , i is the set of adjacent elements in the neighborhood of element i and w ( r ) is the weight function based on a user-specified filter radius R min and the Euclidean distance between the elements i and j ( r i , r j denotes the coordinates of the center point of the respective element).
Accordingly, the sensitivity of the density filter (19) with reference to the design variable is provided by:
ρ ˜ i ρ j = w ( r j ) j N e , i w ( r j ) .
Equation (20) needs to be added to the sensitivity of the goal function according to the chain rule.

3.2.2. Threshold Projection

The threshold projection sets all the density values above a user-defined threshold to 1 and all the values below this threshold to 0. There are several versions of this filter. In [42], the filter is constructed as follows:
ρ ˜ ^ e = tanh ( β η ) + tanh ( β ( ρ ˜ e η ) ) tanh ( β η ) + tanh ( β ( 1 η ) ) ,
where η is the projection threshold and β is a projection parameter tuning the ”sharpness” of the projection (cf. Figure 4).
The corresponding elemental correction term for the sensitivity is provided by:
ρ ˜ ^ e ρ ˜ e = 1 β cosh 2 ( β ( ρ ˜ e η ) ) tanh ( β η ) + tanh ( β ( 1 η ) ) .

4. Optimization Results

The optimization is initialized with a design space that measures 80 × 78   mm 2 . The design space is discretized with 1 × 1 × 1   mm 3 linear cubic finite elements. From this design space, we exclude regions reserved for the application of commercially available piezoelectric patches as well as pre-defined gaps between the inner and outer beam, as well as between the inner beam and the fixed end. Figure 5 depicts the setting for the optimization. The design space is marked light gray. Both masses (M) are modeled as point masses, as we found this approximation to be sufficiently accurate for our application. Finally, the piezoelectric patches are modeled as a secondary layer on top of the design space using the material parameters provided in [43]. The material parameters are presented in Table 1.
For this contribution, the target values for the structure’s resonant frequencies are set to f 1 = 74 Hz and f 2 = 76 Hz . Furthermore, the penalty factor is set to p = 2 and the volume fraction to α = 0.7 . At the beginning of the optimization, β = 1 and this parameter is subsequently increased during optimization until β = 8 . The optimization is stopped either if the maximum change of elemental density is less than 1% or when the number of iterations reaches 100. The resulting geometry is presented in Figure 6. The structure’s resonant frequencies are f 1 = 73.59 Hz and f 2 = 74.48 Hz . The relative difference in dielectric energy between these resonant modes is g 2 = 15.03 % .
To achieve a manufacturable design, the final geometry is cut-off at ρ e = α (volume fraction), i.e., elements with ρ e > α are considered full and those with ρ e α are considered void. The boundary of the geometry is smoothed as the final step of this conversion. Figure 7 shows the mode shapes of the first two resonant modes of the structure after conversion. They value at f 1 = 74.573 Hz and f 2 = 79.054 Hz , respectively. The corresponding open-circuit voltage and maximum power outputs are presented in Figure 8 and Figure 9. For the computation of the electrical power output, we assumed optimally matched load resistances at each frequency. Note that both outer patches are connected in parallel.

5. Discussion

As presented in Section 4, the intended objectives are not achieved by the optimization. This is to be expected, as the optimization result can only partially fulfill both objective contributions ( g 1 and g 2 ). Furthermore, we observe a change in the resonant frequency f 2 after deriving a manufacturable geometry from the optimization result. However, both values are still acceptable for our intended application.
After analyzing the power output of the manufacturable geometry using a harmonic analysis (Figure 9), we observe a bigger relative difference in the power amplitudes at the two resonant frequencies compared to the optimization process. During the optimization, the structure’s dielectric energy has been used as a power potential indicator (as implemented in g 2 ). After optimization, the relative difference of the dielectric energy g 2 at both resonant frequencies was minimized to 15 % , whereas the relative difference in the power output values to 124%. There are several factors contributing to this discrepancy. Firstly, we use the eigenvectors for the evaluation of the dielectric energy of the system in (11). This is computationally efficient, as the mode shapes, i.e., the eigenvectors, are directly obtained during the computation of the resonant frequencies. By definition, the nodal displacements and potentials in the eigenvectors represent the mode shape and cannot be evaluated quantitatively. Hence, proper scaling of the eigenvectors would improve the accuracy of the indicator. Such an improvement can be implemented by a harmonic analysis within the optimization loop, which motivates the following discussion.
Our experimental observations revealed that the damping of a multi-resonant structure is mode-specific and generally dependent on frequency. It was not possible to determine a single damping ratio, which applies to the whole frequency range under investigation. Common practice is to adjust the damping model, such that the simulation results match the experimental findings. Due to the unavailability of a reliable extent of damping, we have to state that the additional computational effort for the calculating power output from a harmonic analysis is of limited use. Thereby, we favored the use of dielectric energies as obtained from the eigenvectors for power estimation. Only experimental characterization can provide mode-specific damping information and resulting power output.

6. Conclusions and Outlook

In this work, we performed a TO of a multi-resonant PEH. For the optimization, we proposed a novel combined-objective function with the following requirements: (1) the final structure should have its first two resonant frequencies at f 1 = 74 Hz and f 2 = 76 Hz and (2) the relative difference in the power output value corresponding to the two modes should be minimal. To obtain a manufacturable design, a regularization filter for the elemental material density has been implemented. The optimization has converged to a promising structure that fulfills the requirements. In post-processing, the optimal result is converted into a manufacturable geometry. The subsequent finite element simulation of this manufacturable design shows satisfactory results for the structure’s resonant frequencies. The corresponding power output amplitudes can be improved. However, these values may not be a true representation of the device’s real-world performance as, in the course of simulations, constant damping was assumed. This assumption is not realistic but is a common practice, since accurate frequency-dependent damping models can only be obtained experimentally. Therefore, only an experimental characterization can determine the actual power output level of the device.
As a next step, we will fabricate the optimized design obtained from this work and perform experimental characterizations for further validation of the design approach. Furthermore, we will implement a harmonic analysis of the power output values within the optimization loop for a more accurate estimation of the output power. The additional computational effort will be minimized by the state-of-the-art model order-reduction techniques.

Author Contributions

Conceptualization, S.H. and T.B.; inspiration and discussion during conceptualization, K.C.N., D.H. and J.G.K.; methodology, S.H.; software, S.H. and U.F.; simulation and data processing, S.H.; original draft preparation, S.H.; writing—review and editing, S.H., D.H. and T.B.; supervision, D.H. and T.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by Deutsche Forschungsgemeinschaft (DFG), the German Research Foundation—Project number 446995632.

Data Availability Statement

The data presented in this study are contained within the article. No additional data were created or analyzed in this study.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
PEHPiezoelectric Energy Harvester
TOTopology Optimization

Appendix A

The generalized eigenvalue of a Matrix pencil ( K , M ) is defined by:
K ϕ i = λ i M ϕ i ( K λ i M ) ϕ i = 0 .
Therefore, the derivative w.r.t to ρ e can be computed by:
( K λ i M ) ρ e ϕ i + ( K λ i M ) ϕ i ρ e = 0 .
By multiplying with ϕ i T from the left, we can eliminate the second part of the equation:
ϕ i T ( K λ i M ) ρ e ϕ i = 0
ϕ i T K ρ e ϕ i ϕ i T λ i ρ e M ϕ i ϕ i T λ i M ρ e ϕ i = 0 .
As ϕ i T M ϕ i = 1 by normalization, it follows:
λ i ρ e = ϕ i , e T K e ρ e λ i M e ρ e ϕ i , e .

Appendix B

The derivative of p e l can be computed as:
p e l ρ e = ρ e 1 2 ϕ e l T K e l ϕ e l = 1 2 ϕ e l T K e l ϕ e l ρ e + 1 2 ϕ e l T K e l ρ e ϕ e l .
K e l ρ e = 0 since dielectric properties do not change with the change of design space. ϕ e l ρ e can be computed using the adjoint method. The harmonic equilibrium equation from the System (3) for an arbitrary excitation frequency f reads:
K m λ M m K m , e l K m , e l T K e l = : K ¯ ϕ m ϕ e l = : ϕ = B u 0 ,
with λ = ( 2 π f ) 2 . Its derivative w.r.t ρ e is provided by:
K ¯ ρ e ϕ + K ¯ ϕ ρ e = 0 ϕ ρ e = K ¯ 1 K ¯ ρ e ϕ .
Now, let [ μ m , μ e l ] T solve the following adjoint system:
K ¯ μ m μ e l μ = 0 1 2 K e l ϕ e l .
Combining (A8) with (A9) and using the symmetry of K ¯ , (A6) can be rewritten as:
p e l ρ e = μ T K ¯ ρ e ϕ
From here, only:
p e l ρ e = μ e T K e ρ e λ M e ρ e ϕ e
remains, as the piezoelectric properties stored K m , e l is also independent from ρ e and μ e l , e = 0 and ϕ e l , e = 0 for all the design space elements (steel structure).

References

  1. Sah, D.K.; Amgoth, T. Renewable energy harvesting schemes in wireless sensor networks: A Survey. Inf. Fusion 2020, 63, 223–247. [Google Scholar] [CrossRef]
  2. Mohanty, A.; Parida, S.; Behera, R.K.; Roy, T. Vibration energy harvesting: A review. J. Adv. Dielectr. 2019, 9, 1930001. [Google Scholar] [CrossRef]
  3. Adendorff, K.; Naidoo, R.; Jimoh, A.A.; Okojie, D.E. A hybrid piezoelectric micro-power generator for use in low power applications. Renew. Sustain. Energy Rev. 2015, 49, 1136–1144. [Google Scholar]
  4. Bowen, C.; Kim, H.; Weaver, P.; Dunn, S. Piezoelectric and ferroelectric materials and structures for energy harvesting applications. Energy Environ. Sci. 2014, 7, 25–44. [Google Scholar] [CrossRef] [Green Version]
  5. Li, H.; Tian, C.; Deng, Z.D. Energy harvesting from low frequency applications using piezoelectric materials. Appl. Phys. Rev. 2014, 1, 041301. [Google Scholar] [CrossRef] [Green Version]
  6. Toprak, A.; Tigli, O. Piezoelectric energy harvesting: State-of-the-art and challenges. Appl. Phys. Rev. 2014, 1, 031104. [Google Scholar] [CrossRef]
  7. Liu, H.; Zhong, J.; Lee, C.; Lee, S.W.; Lin, L. A comprehensive review on piezoelectric energy harvesting technology: Materials, mechanisms, and applications. Appl. Phys. Rev. 2018, 5, 041306. [Google Scholar] [CrossRef] [Green Version]
  8. Tang, L.; Yang, Y.; Soh, C.K. Toward broadband vibration-based energy harvesting. J. Intell. Mater. Syst. Struct. 2010, 21, 1867–1897. [Google Scholar] [CrossRef]
  9. Qin, H.; Mo, S.; Jiang, X.; Shang, S.; Wang, P.; Liu, Y. Multimodal Multidirectional Piezoelectric Vibration Energy Harvester by U-Shaped Structure with Cross-Connected Beams. Micromachines 2022, 13, 396. [Google Scholar] [CrossRef] [PubMed]
  10. Caetano, V.J.; Savi, M.A. Multimodal pizza-shaped piezoelectric vibration-based energy harvesters. J. Intell. Mater. Syst. Struct. 2021, 32, 2505–2528. [Google Scholar] [CrossRef]
  11. Caetano, V.J.; Savi, M.A. Star-shaped piezoelectric mechanical energy harvesters for multidirectional sources. Int. J. Mech. Sci. 2022, 215, 106962. [Google Scholar] [CrossRef]
  12. Pertin, O.; Shrivas, P.; Guha, K.; Rao, K.S.; Iannacci, J. New and efficient design of multimode piezoelectric vibration energy harvester for MEMS application. Microsyst. Technol. 2021, 27, 3523–3531. [Google Scholar] [CrossRef]
  13. Bouhedma, S.; Zheng, Y.; Lange, F.; Hohlfeld, D. Magnetic frequency tuning of a multimodal vibration energy harvester. Sensors 2019, 19, 1149. [Google Scholar] [CrossRef] [Green Version]
  14. Bouhedma, S.; Hu, S.; Schütz, A.; Lange, F.; Bechtold, T.; Ouali, M.; Hohlfeld, D. Analysis and Characterization of Optimized Dual-Frequency Vibration Energy Harvesters for Low-Power Industrial Applications. Micromachines 2022, 13, 1078. [Google Scholar] [CrossRef]
  15. Wang, H.; Tang, L. Modeling and experiment of bistable two-degree-of-freedom energy harvester with magnetic coupling. Mech. Syst. Signal Process. 2017, 86, 29–39. [Google Scholar] [CrossRef]
  16. Xiao, H.; Wang, X.; John, S. A dimensionless analysis of a 2DOF piezoelectric vibration energy harvester. Mech. Syst. Signal Process. 2015, 58, 355–375. [Google Scholar] [CrossRef]
  17. Sarker, M.R.; Julai, S.; Sabri, M.F.M.; Said, S.M.; Islam, M.M.; Tahir, M. Review of piezoelectric energy harvesting system and application of optimization techniques to enhance the performance of the harvesting system. Sens. Actuators A Phys. 2019, 300, 111634. [Google Scholar] [CrossRef]
  18. Park, J.; Lee, S.; Kwak, B.M. Design optimization of piezoelectric energy harvester subject to tip excitation. J. Mech. Sci. Technol. 2012, 26, 137–143. [Google Scholar] [CrossRef]
  19. Rui, X.; Li, Y.; Liu, Y.; Zheng, X.; Zeng, Z. Experimental study and parameter optimization of a magnetic coupled piezoelectric energy harvester. Appl. Sci. 2018, 8, 2609. [Google Scholar] [CrossRef] [Green Version]
  20. Yang, H.; Wei, Y.; Zhang, W.; Ai, Y.; Ye, Z.; Wang, L. Development of piezoelectric energy harvester system through optimizing multiple structural parameters. Sensors 2021, 21, 2876. [Google Scholar] [CrossRef]
  21. Bouhedma, S.; Rao, Y.; Schütz, A.; Yuan, C.; Hu, S.; Lange, F.; Bechtold, T.; Hohlfeld, D. System-level model and simulation of a frequency-tunable vibration energy harvester. Micromachines 2020, 11, 91. [Google Scholar] [CrossRef] [Green Version]
  22. Hu, S.; Bouhedma, S.; Schütz, A.; Stindt, S.; Hohlfeld, D.; Bechtold, T. Design optimization of multi-resonant piezoelectric energy harvesters. Microelectron. Reliab. 2021, 120, 114114. [Google Scholar] [CrossRef]
  23. Bendsøe, M.P.; Kikuchi, N. Generating optimal topologies in structural design using a homogenization method. Comput. Methods Appl. Mech. Eng. 1988, 71, 197–224. [Google Scholar] [CrossRef]
  24. Bendsoe, M.P.; Sigmund, O. Topology Optimization: Theory, Methods, and Applications; Springer: Berlin/Heidelberg, Germany, 2003. [Google Scholar]
  25. Chen, F.; Wang, J.; Yang, X. Topology optimization design and numerical analysis on cold plates for lithium-ion battery thermal management. Int. J. Heat Mass Transf. 2022, 183, 122087. [Google Scholar] [CrossRef]
  26. Sharma, A.K.; Kosta, M.; Shmuel, G.; Amir, O. Gradient-based topology optimization of soft dielectrics as tunable phononic crystals. Compos. Struct. 2022, 280, 114846. [Google Scholar] [CrossRef]
  27. Díaaz, A.R.; Kikuchi, N. Solutions to shape and topology eigenvalue optimization problems using a homogenization method. Int. J. Numer. Methods Eng. 1992, 35, 1487–1502. [Google Scholar] [CrossRef]
  28. Ma, Z.D.; Kikuchi, N.; Hagiwara, I. Structural topology and shape optimization for a frequency response problem. Comput. Mech. 1993, 13, 157–174. [Google Scholar] [CrossRef]
  29. Ma, Z.D.; Cheng, H.C.; Kikuchi, N. Structural design for obtaining desired eigenfrequencies by using the topology and shape optimization method. Comput. Syst. Eng. 1994, 5, 77–89. [Google Scholar] [CrossRef] [Green Version]
  30. Hu, S.; Liu, Z.; Bechtold, T. Fast Topology Optimization for Resonating Structures via Generalized Incremental Frequency Method and Modal Superposition-Based Model Order Reduction. In Advances in Structural and Multidisciplinary Optimization, Proceedings of the 13th World Congress of Structural and Multidisciplinary Optimization, Beijing, China, 20–24 May 2019; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
  31. Giannini, D.; Braghin, F.; Aage, N. Topology optimization of 2D in-plane single mass MEMS gyroscopes. Struct. Multidiscip. Optim. 2020, 62, 2069–2089. [Google Scholar] [CrossRef]
  32. Giannini, D.; Aage, N.; Braghin, F. Topology optimization of MEMS resonators with target eigenfrequencies and modes. Eur. J. Mech.-A/Solids 2022, 91, 104352. [Google Scholar] [CrossRef]
  33. He, M.; Zhang, X.; dos Santos Fernandez, L.; Molter, A.; Xia, L.; Shi, T. Multi-material topology optimization of piezoelectric composite structures for energy harvesting. Compos. Struct. 2021, 265, 113783. [Google Scholar] [CrossRef]
  34. He, M.; He, M.; Zhang, X.; Xia, L. Topology Optimization of Piezoelectric Energy Harvesters for Enhanced Open-Circuit Voltage Subjected to Harmonic Excitations. Materials 2022, 15, 4423. [Google Scholar] [CrossRef]
  35. Ikeda, T. Fundamentals of Piezoelectricity/by Takuro Ikeda; Oxford University Press: Oxford, UK, 1990. [Google Scholar]
  36. Hu, S.; Fitzer, U.; Stindt, S.; Bechtold, T. Topology Optimization of a Folded Beam Piezoelectric Energy Harvester. IFAC-PapersOnLine 2022, 55, 379–384. [Google Scholar] [CrossRef]
  37. Lerch, R. Simulation of piezoelectric devices by two- and three-dimensional finite elements. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1990, 37, 233–247. [Google Scholar] [CrossRef]
  38. Zheng, B.; Chang, C.J.; Gea, H.C. Topology optimization of energy harvesting devices using piezoelectric materials. Struct. Multidiscip. Optim. 2009, 38, 17–23. [Google Scholar] [CrossRef]
  39. Diaz, A.; Sigmund, O. Checkerboard patterns in layout optimization. Struct. Optim. 1995, 10, 40–45. [Google Scholar] [CrossRef]
  40. Jog, C.S.; Haber, R.B. Stability of finite element models for distributed-parameter optimization and topology design. Comput. Methods Appl. Mech. Eng. 1996, 130, 203–226. [Google Scholar] [CrossRef]
  41. Sigmund, O.; Petersson, J. Numerical instabilities in topology optimization: A survey on procedures dealing with checkerboards, mesh-dependencies and local minima. Struct. Optim. 1998, 16, 68–75. [Google Scholar] [CrossRef]
  42. Wang, F.; Lazarov, B.S.; Sigmund, O. On projection methods, convergence and robust formulations in topology optimization. Struct. Multidiscip. Optim. 2011, 43, 767–784. [Google Scholar] [CrossRef]
  43. dos Santos Guimarães, C.; de Silva Bussamra, F.L.; Pommier-Budinger, V.; Hernandes, J.A. Structural Shape Control Using Macro Fiber Composite Piezoelectric Sensors and Actuators. Mec. Comput. 2010, 29, 8263–8279. [Google Scholar]
Figure 1. Folded-beam design introduced in [13]. The geometry is shown on the left. It is composed of a steel structure (gray), two tip masses (black), and three piezoelectric patches (mustard). The first two mode shapes of the structure are shown on the right.
Figure 1. Folded-beam design introduced in [13]. The geometry is shown on the left. It is composed of a steel structure (gray), two tip masses (black), and three piezoelectric patches (mustard). The first two mode shapes of the structure are shown on the right.
Micromachines 14 00332 g001
Figure 2. An ideal case for a frequency range of efficient operation, formed by two overlapping resonant modes (log scale).
Figure 2. An ideal case for a frequency range of efficient operation, formed by two overlapping resonant modes (log scale).
Micromachines 14 00332 g002
Figure 3. Filter techniques in TO. Left: A TO-obtained structure without any filter. Center: Same optimization with a density filter. Right: Same optimization with both filters active.
Figure 3. Filter techniques in TO. Left: A TO-obtained structure without any filter. Center: Same optimization with a density filter. Right: Same optimization with both filters active.
Micromachines 14 00332 g003
Figure 4. Threshold projection (21) for threshold η = 0.7 and different β values.
Figure 4. Threshold projection (21) for threshold η = 0.7 and different β values.
Micromachines 14 00332 g004
Figure 5. Design space of the problem. The areas beneath the piezo (-electric) patches and the masses M are excluded from the design space.
Figure 5. Design space of the problem. The areas beneath the piezo (-electric) patches and the masses M are excluded from the design space.
Micromachines 14 00332 g005
Figure 6. Optimized geometry. The gray level indicates the value of ρ e for each element.
Figure 6. Optimized geometry. The gray level indicates the value of ρ e for each element.
Micromachines 14 00332 g006
Figure 7. Mode shapes of the first two resonant modes of the optimal structure after the conversion to a manufacturable design value at f 1 = 74.573 Hz and f 2 = 79.054 Hz .
Figure 7. Mode shapes of the first two resonant modes of the optimal structure after the conversion to a manufacturable design value at f 1 = 74.573 Hz and f 2 = 79.054 Hz .
Micromachines 14 00332 g007
Figure 8. Open- circuit voltage output of the structure obtained through harmonic analysis. The structure is excited with 0.2 g at its fixed end. The damping is set to 0.008.
Figure 8. Open- circuit voltage output of the structure obtained through harmonic analysis. The structure is excited with 0.2 g at its fixed end. The damping is set to 0.008.
Micromachines 14 00332 g008
Figure 9. Maximum power output of the structure at optimal load value computed from values presented in Figure 8. Both outer patches are connected in parallel. The total power output is computed as the sum of both patches.
Figure 9. Maximum power output of the structure at optimal load value computed from values presented in Figure 8. Both outer patches are connected in parallel. The total power output is computed as the sum of both patches.
Micromachines 14 00332 g009
Table 1. Mechanical properties of commercial MFC piezoelectric patches [43].
Table 1. Mechanical properties of commercial MFC piezoelectric patches [43].
Density (kg/m 3 )Piezoelectric Constants (C/m 2 )
ρ 4700 e 31 −2.227
Young’s Modulus (GPa) e 32 −0.671
E 1 45.21 e 33 16.665
E 2 12.39 e 24 0.0258
E 3 40.44 e 15 13.668
Shear Modulus (GPa)Dielectric Relative Constants
G 12 6.03 ε 11 T / ε 0 1574.8
G 23 6.68 ε 22 T / ε 0 24.7
G 31 17.01 ε 33 T / ε 0 1528.7
Poisson’s Ratio
ν 12 0.39
ν 23 0.17
ν 13 0.44
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.

Share and Cite

MDPI and ACS Style

Hu, S.; Fitzer, U.; Nguyen, K.C.; Hohlfeld, D.; Korvink, J.G.; Bechtold, T. Multi-Objective Topology Optimization of a Broadband Piezoelectric Energy Harvester. Micromachines 2023, 14, 332. https://doi.org/10.3390/mi14020332

AMA Style

Hu S, Fitzer U, Nguyen KC, Hohlfeld D, Korvink JG, Bechtold T. Multi-Objective Topology Optimization of a Broadband Piezoelectric Energy Harvester. Micromachines. 2023; 14(2):332. https://doi.org/10.3390/mi14020332

Chicago/Turabian Style

Hu, Siyang, Ulrike Fitzer, Khai Chau Nguyen, Dennis Hohlfeld, Jan G. Korvink, and Tamara Bechtold. 2023. "Multi-Objective Topology Optimization of a Broadband Piezoelectric Energy Harvester" Micromachines 14, no. 2: 332. https://doi.org/10.3390/mi14020332

APA Style

Hu, S., Fitzer, U., Nguyen, K. C., Hohlfeld, D., Korvink, J. G., & Bechtold, T. (2023). Multi-Objective Topology Optimization of a Broadband Piezoelectric Energy Harvester. Micromachines, 14(2), 332. https://doi.org/10.3390/mi14020332

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop