Next Article in Journal
Research on Dynamic Risk Assessment and Active Defense Strategy of Active Distribution Network under Ice Weather
Next Article in Special Issue
Investigation and Enhancement of the Detectability of Flaws with a Coarse Measuring Grid and Air Coupled Ultrasound for NDT of Panel Materials Using the Re-Radiation Method
Previous Article in Journal
A Simplified Mathematical Model for the Analysis of Varying Compliance Vibrations of a Rolling Bearing
Previous Article in Special Issue
A Bio-Inspired Flapping Wing Rotor of Variant Frequency Driven by Ultrasonic Motor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling of Electromagnetic Energy Harvester with Rotational Pendulum Using Mechanical Vibrations to Scavenge Electrical Energy

by
Bartłomiej Ambrożkiewicz
*,
Grzegorz Litak
and
Piotr Wolszczak
Department of Automation, Faculty of Mechanical Engineering, Lublin University of Technology, Nadbystrzycka 36, 20-618 Lublin, Poland
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(2), 671; https://doi.org/10.3390/app10020671
Submission received: 22 December 2019 / Revised: 14 January 2020 / Accepted: 15 January 2020 / Published: 17 January 2020

Abstract

:

Featured Application

A concept of a non-linear electromagnetic energy harvesting (EH) system with a rotational pendulum was proposed. In the system, non-linearity was introduced through the interaction of additional stationary magnets placed into coils with a tip magnet on the pendulum by creating double potential well. In the system, both oscillations and rotations are possible depending on the system parameters including the strength of magnets, the excitation conditions including amplitude and frequency, and the initial conditions. Moreover, a housing for an EH prototype was made in a 3D printing technique in order to avoid undesirable electromagnetic coupling. Data obtained from a finite element method (FEM) model, real data, and data from an experiment were applied to reflect the system’s operation. The proposed EH device can be used in applications, where a kinematic longitude excitation with fairly small frequency is available.

Abstract

A concept of non-linear electromagnetic system with the rotational magnetic pendulum for energy harvesting from mechanical vibrations was presented. The system was stimulated by vertical excitation coming from a shaker. The main assumption of the system was the montage of additional regulated stationary magnets inside coils creating double potential well, and the system was made with a 3D printing technique in order to avoid a magnetic coupling with the housing. In validation process of the system, modelling of electromagnetic effects in different configurations of magnets positions was performed with the application of a finite element method (FEM) obtaining the value of magnetic force acting on the pendulum. A laboratory measurement circuit was built and an experiment was carried out. The voltage and power outputs were measured for different excitations in range of system operational frequencies found experimentally. The experimental results of the physical system with electrical circuit and numerical estimations of the magnetic field of a stationary magnet’s configuration were used to derive a mathematical model. The equation of motion for the rotational pendulum was used to prove the broadband frequency effect.

Graphical Abstract

1. Introduction

Through the last decade, research on a new field of activity, energy harvesting (EH), has gained recognition both in scientific and industrial environments. It denotes the effort to develop small distributed renewable sources of energy and the need to power up low-voltage autonomous electrical devices and batteries. Many EH devices have already found application in a few technical solutions i.e., automotive [1], wireless sensor networks (WSN) [2], health-monitoring devices [3], and civil engineering [4]. The main advantage of the aforementioned device is its substitutability for a battery as a source of energy. The energy is scavenged from the environment, examples of which can be solar energy, tides, gravitational force or mechanical vibrations [5]. Most of the currently designed systems are electromagnetic, piezoelectric and electrostatic [5,6,7,8,9,10,11].
From the application point of view, research aimed at increasing the density of recovered energy is important for systems generating small portions of energy in a fairly wide frequency band [12,13,14]. That is why energy harvesting systems with non-linear behavior are researched [15,16,17,18,19,20]. The main advantages of non-linearities introduced to the system are the possibilities to modify the resonance behavior [21,22], to obtain good off-resonance performance [23,24,25] or to attain multi-resonance [26].
Our proposed concept is based on an electromagnetic system with a rotational pendulum. Energy harvesters with a rotational pendulum as a resonator can evolve in time with a large amplitude. Earlier works on the system’s development started from theoretical analysis of its dynamics and studies on its response by different excitation conditions [27,28,29,30,31,32]. Next, several systems were built, optimized and tested for different applications, primarily energy harvesting from human motion [33,34,35,36,37] and automotives [38,39]. In our system we focused on introducing double potential well by mounting adjustable additional magnets inside coils expecting non-linear behavior brought by magnetic interactions. The magnet’s position and its repelling alignment will enable pendulum’s rotation only by specific excitation frequency which was found experimentally. To obtain a pure system’s response, the corresponding parts were made with a 3D printing technique. Some improvements connected with design had to be introduced after performing tests. The system was validated with mathematical model and variables used in the pendulum’s equation of motion were based on simulation and experiment. Voltage and power outputs were measured with various excitation frequencies.

2. Materials and Methods

For the implementation of a specific research task, the idea of electromagnetic system with a rotational pendulum for energy harvesting from mechanical vibrations has been proposed. A prototype of the device was made for carrying out a physical experiment and a model for carrying out numerical simulations was developed. The experiment’s methodology started with the performance of the finite element method (FEM) model, where the value of the magnetic force acting on a pendulum was found. The physical experiment provided information on operational frequency range, acceleration amplitude acting on the system during test and rotational velocity of a pendulum. Values found in the experiment were used in formulation equation of motion. Other values were calculated from system’s geometry, such as moment of inertia and distance from pivot to the center of gravity. The system’s operation was validated with the adjusted mathematical model.

2.1. Prototype and Physical Experiment

The prototype of the energy harvester with a rotational pendulum is based on electromagnetic induction phenomena. To rid it of magnetic coupling, the system’s housing and pendulum’s housing are made of ABS (Acrylonitrile Butadiene Styrene) in 3D printing technique (Figure 1). Many steps are taken to miniaturize the EH systems, because in case of electromagnetic systems this is hard to do. Our goal was not to exceed 10 cm of system’s design in any direction, what was succeeded. The system consists of 12 coils connected in series mounted symmetrically on both sides (6 on one side) of the energy harvester. By such a connection we wanted to obtain the biggest possible voltage/power response from the system. The system’s design provide to apply versatile coils connection, exemplary connecting 6 coils mounted on one wall of the system. Inside 4 of them, additional adjustable permanent magnets are installed, two of them on each side of the system (Figure 2). They will have an influence on electromagnetic induction of the system and movement of the pendulum by creating double potential well. Because the initial design of the pendulum was too light, additional inertial masses had to be mounted on the system’s sides to increase its moment of inertia. The total weight of entire system was equal to ca. 900 g. In the system, both mechanical and electrical damping will occur. Mechanical damping was in form of friction inside bearings, and an electrical one from additional resistors. The first one was only considered later in the formulation of the equation of motion.
In the prototype presented above it is possible to adjust a few factors such as the position of additional magnets, additional inertial masses of pendulum and the resistance load attached to the electric circuit. The aforementioned variables will have an influence on output results: voltage, power, magnetic force and magnetic induction in the system. Change of magnetic force and electromagnetic induction were simulated in Finite Element Method Magnetics (FEMM) software. The principle of calculations will be explained in next sections. Because we used the strongest available magnets made of neodymium N37 for the magnetic forces acting between them, in the experiment, results were only obtained for the furthest distance between them. Then, the additional load was regulated and output voltage and power were measured. Such measurements were performed for different frequencies of excitation to obtain chart output voltage/power in the frequency of excitation domain.
The scheme of measuring setup is presented in Figure 3. The measurement system consisted of the following equipment: signal generator—KZ 1406 (ZOPAN, Warsaw, Poland), signal amplifier—BAA 120 (TIRA GmbH, Schalkau, Germany), vibration exciter—S 513 (TIRA, Schalkau, Germany), oscilloscope—MDO3024 (Tektronix, Beaverton, OR, USA), accelerometer—AltIMU-10 v3 (Pololu, Las Vegas, NV, USA). Entire measuring system starts with signal generator, where periodic sinusoidal signal was generated with specific frequency matched to the operational frequency of the energy harvester, where the rotational movement of the pendulum could be sustained. Next, the signal was processed with a signal amplifier, and such a signal was steering the shaker’s movement exciting the energy harvester. Signals of output voltage and output power were transferred to the oscilloscope. Different values of output voltage and power were obtained by additional resistance in the form of a potentiometer. On the energy harvester housing, an accelerometer was mounted, the measured amplitude of acceleration acting on the system was equal to 0.55 g, by input current equal to 1.1 A, and input voltage of 4.0 V displayed on the amplifier. The value of acceleration amplitude was taken into the mathematical model. Below in Table 1 and Table 2, the most important characteristics of used vibration exciter and accelerometer are given.

2.2. Numerical Simulations

The numerical calculations of electromagnetic induction and magnetic force values were performed using the finite element method. Namely, the FEMM software (Ph.D. David Meeker, Waltham, MA, USA) was used, which consider magnetic, electrostatic, heat flow and current flow problems. The first simulation task was the change of magnetic induction value in the angular and linear position change of the additional magnet inside the coils and mounted on the pendulum. Then, we received results in the form of a 3D surface presenting the change of electromagnetic phenomena in the system. The estimation was performed for a coil, an additional permanent magnet and the permanent tip magnet mounted on the rotating pendulum keeping the geometrical and magnetic features of applied magnets (NdFeB 37) and coils (28 AWG); also, the repulsive orientation of magnets was taken into account. The aim of the simulation is to check how the change of additional magnet position influence on magnetic induction arose in the coil [40,41,42]. In the Figure 4 the layout of the magnets and coil and exemplary results in FEMM software are presented.
In total, created surface will consist of 900 results, changing the position of magnet inside coil by 0.5 mm and position of magnet on pendulum by 2.5°. The resultant value of electromagnetic induction is taken from the coil in the x-axis. The magnetic induction results surface is presented in Figure 5 and Figure 6.
The position of the additional magnet was limited by the geometry of the plastic cylindrical channels on the sides of the system. Value discrepancies of magnetic induction are the result of magnetic interactions between magnets. It must be emphasized that only two pairs of parallel coils with additional magnets were installed. In other coils without additional magnets, the value of electromagnetic induction would be smaller.
The second simulation task using the FEMM software was to find the value of magnetic force between acting magnets and use such value for modelling of mechanical system. The manual of FEMM software predicts the application of co-energy of entire energetic system for two different close positions [43]. The magnetic force equation has the following form:
F m a g   =   δ W c δ l ,
where δWc is variation of co-energy value of entire energetic system for two different close positions, Fmag is the magnetic force acting on the pendulum, δl is the variation of the additional magnet position.
To compare the magnetic forces’ values in the system for different distances between additional magnet and the pendulum, Figure 7 shows the curves for x = 1 mm and x = 13 mm. The differences between the curves in Figure 7 are caused by the magnetic interactions. The shape with three local maxima results from the rectangular shape of cross-section of the magnets in 2D simulation. The amplitude and distribution of local maxima changes with the distance between the magnets. In Figure 7b, the local maxima disappear due to the decrease in amplitude of the magnetic force for the 13 mm distance. Two-dimensional simulation provided very important information concerning change of magnetic force depending on the distance between magnets and the obtained values can be used in the equation of motion.
The calculated resultant value of magnetic force acting on the pendulum was applied to the equation of motion of the system. The aforementioned magnetic force comprised in some way the damping of the pendulum’s movement. That is why its value is so important in modelling. Because of the appearance of magnetic forces that were too strong by close relative alignment of the magnets, in the experiment and in the mathematical model only the long distance between the magnets was considered.

3. Experimental Results

Using the measuring system described in the previous section, measurements of output voltage and output power were performed in the domain of variable operating frequency [44,45]. Also measurements were conducted in open loop (directly, without additional resistance). First tests showed that the movement of the pendulum is supported by frequencies from 10 Hz to 15 Hz. Experimental tests gave important information about which environment the system can succeed in. The disadvantage can be short broadband of operational frequency, but it is possible to broaden it by application of heavier masses on additional inertial pendulums. Output signals were displayed on the oscilloscope and its results were stored. Then, visualization and proper calculations were performed in Matlab software. Measured values of output voltage were stored with a sampling rate of 100,000 S/s per channel. For the tests, seven different values of additional resistance were considered mainly close to the value of internal resistance of 12 coils in the system connected in series which value was 21.5 Ω. The most relevant characteristics of performed tests are collected in Table 3.
After collecting the raw signal of output voltage, we calculated its root mean square value (2). After such calculations output power was calculated (3). All the results were visualized in the operating frequency domain (Figure 8).
U R M S   =   1 N k   =   0 N 1 u ( k ) 2
P   =   U R M S 2 R
An increase of the additional resistance value causes the increase of the voltage RMS value from the excited system. However, it can be expected that the mentioned value would not be as great and reach almost the limit of 30.3 Ω. Cross-checking measurements must be performed for higher values of additional resistance or by wider range of operating frequency. Unfortunately, design constraints hindered it. The results of RMS voltage value in open loop show its linear change with operational frequency. Tests with the open loop are not in practical usage, but they can show the limits of the system. The above presented results were only performed for the frequency range, where rotational movement was supported. To expand it, pendulum’s moment of inertia must be changed. Moreover, for the optimization of the system’s power efficiency, extended experiments can be performed in wider range of additional resistance [46].

4. Mechanical Model of the Mechanical Resonator and a Broadband Effect

Research on small portions of mechanical energy conversion into electrical energy was started from mechanical vibrations. Results obtained were adjusted to needs of specific energy receiver paying special attention for output power, resonance frequency or energy conversion. Continuous development favors usage of new materials and systems managing conversion and flow of energy to the receiver. Nowadays, in the beginning of work on the design of new device, mathematical model is derived. Few papers have been dedicated only to mathematical modeling of different types of EH and provide crucial and useful information on this topic [47,48,49].
The mechanical model of the energy harvester set with a rotational pendulum can be represented by the equation of motion [50,51,52]. Into the equation of motion (4) for a pendulum the equation of motion correcting factor was introduced in order to reflect real situation from the experiment. To keep the rotational movement of the pendulum, the entire system with a vibration exciter was tilted, and for this angle φ is added to the equation of motion. We did not consider the equation for the magnetic circuit. In our simplified model we took into account all the variables having an impact on pendulum’s movement and they are contained in the following simplified equation of motion:
I φ ¨ + c φ ¨ + m ω 2 A h cos ( ω t ) sin φ + m g h sin φ + F m a g s i n ( 2 φ + 1.5 ) = 0 ,
where: I—moment of inertia of pendulum, c—damping coefficient, h—distance from center of gravity to pivot axis, m—mass of pendulum and rotational additional masses, g—gravitation, ω—angular velocity of pendulum, A—amplitude of vertical excitation, Fmag—magnetic force amplitude acting on pendulum, φ—tilt angle of pendulum.
The phase shift of 1.5 rad present in the magnetic force term denotes the tilting of the pendulum (Figure 2). Here, we simplified the angular dependence of the magnetic force using sinusoidal dependence.
Some selected features of the Hamiltonian version of the considered model together with the separatrices which defines the regions of two basic solutions: oscillations and rotations are presented in Figure 9a–c. Note that the double well potential (Figure 9c) is generated by the interaction between the pendulum tip permanent magnet and additional permanent magnets distributed in the system as presented in Figure 2. In the studied system the magnetic forces are much stronger than gravity, which leads to a small splitting of the separatrices (see the differences of red and black lines in Figure 9b). The numerical results of the frequency sweeps are presented in Figure 10.
Figure 10 shows the frequency sweep results for pendulum response in its angular velocity. Our interest was focused on the angular velocity and the corresponding variability that are directly proportional to the electromotive forces in the electrical circuit. We used nodal initial conditions defined by the vertical position of the pendulum for initial angular displacement and two different initial conditions for angular velocity: zero and fairly large 80 rad/s. Two different excitation amplitudes were used, A = 0.002 and A = 0.05 m.
In Figure 10a we show the results of standard deviation and mean value for nodal initial velocity and small excitation. We noticed that the pendulum oscillates with small amplitude with visible resonance at the frequency of 6.5 Hz. This frequency corresponds to the natural frequency of the mechanical resonator. Simultaneously, the mean angular velocity is zero. In Figure 10b we present the results of calculations repeated for large initial velocity and we observe another response in angular velocity. It remains finite up to frequency f of about 24 Hz and collapses to the solution presented in Figure 10a. During the fairly large angular velocity response of 54 rad/s, the variability of this velocity was fairly small. Note that these results were obtained for the excitation amplitude which can be generated on the shaker. To complete the story, we enlarged the amplitude of excitation and obtained richer responses. In Figure 10b, we can observe a number of higher resonances which multiplies f = 6.5 Hz. This is a typical behavior for a non-linear system. Non-linearity based higher order resonances manifest with large enough excitation. On the other hand higher value of initial velocity drives the system to rotations. Interestingly, we observe a wide spectrum of rotational solutions separated by rotational instabilities leading to oscillatory vibrations. In addition, various rotational solutions appear. For small frequency f < 6 Hz we observe a similar dynamical response to the situation presented in Figure 10a. The output velocity is stable with changing frequency raging about 54 rad/s. This solution is related to overcoming the potential barrier (Figure 9c). Referring to Figure 9a,b it is easy to find that a velocity of 54 Hz is sufficient to overcome the potential barrier and obtain a rotational solution. This solution resembles the type of vibration coherence resonance [20], where excitation frequency is less important than the amplitude important to generate a large orbit motion with a minimum of excitation energy necessary to overcome the potential barrier. For higher frequency intervals such as [9 Hz, 15 Hz] and [35 Hz, 45 Hz], the response angular velocity is proportional to excitation frequency. Interestingly, for the first mentioned interval, proportionality is 1:1 while for the next it is 1:1/3. This is evidence that the solution has a sub-harmonic 3 property. In a summary of the results of Figure 10, one can conclude that the system has a great capability toward frequency broad band responses. In spite of the strong model simplifications, the simulations shows the large frequency region of rotational solutions. Such sub-harmonic solutions can be very useful for energy harvesting as proved by [17,53].

5. Conclusions

This paper presents an electromagnetic energy harvester with a rotational pendulum. The entire work presents a logical sequence from the system design, to the application of the finite element method to check magnetic interactions in the system, through testing and mathematical modelling. Performance of the FEM model showed magnetic interactions, which appear in the system and what are the expected magnetic forces by usage of co-energy. Results of calculated magnetic forces served to perform a simplified model of the rotational pendulum reflecting its kinematics in tests. Also, the application of the mentioned method let to us obtaining the surface of electromagnetic induction in domain of magnets’ positions. Owing to this, three specific areas of interest were found depending on the position of additional magnets in coils. By the conducted tests, information about output voltage RMS and output power from the system was obtained. An experiment was performed with different values of operational frequencies and for different values of additional resistance, presenting in such a way a map of results. In the last stage, a simplified mechanical model of rotational pendulum was derived. To formulate the pendulum’s equation of motion, an element from the experiment was added in the form of acting magnetic forces in the system and some phase shift was needed to start the pendulum’s permanent rotational movement. Finally, simulations checked the application potential of the proposed device in the context of frequency broadness, which is important for variable ambient working conditions. However, one should remember that the large orbit (more energetic) is only one possible solution as a non-linear system is known for multiple solutions. To check the availability of the larger orbit solution, one needs to perform more studies including basins of attraction approach [54]. The system can be still optimized by changing its moment of inertia, the strength of magnets, improving the electrical circuit, and by carrying out the tests in wider range of frequencies. The next step in research on the system will be the creation of multiple potential wells, by increasing number of coils with magnets and studying multi-stable phenomena as in the oscillating energy harvesting systems [18,55].

Author Contributions

B.A.: methodology, validation, formal analysis, writing—original draft preparation, visualization, G.L.: methodology, investigation, software, formal analysis, supervision, P.W.: formal analysis, investigation, data curation, writing—review and editing, project administration, funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

Applsci 10 00671 i001This publication was supported by the program of the Polish Ministry of Science and Higher Education under the project DIALOG 0019/DLG/2019/10 in the years 2019–2021.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maurya, D.; Kumar, P.; Khaleghian, S.; Sriramdas, R.; Kang, M.G.; Kishore, R.A.; Kumar, V.; Song, H.C.; Park, J.M.; Taheri, S.; et al. Energy harvesting and strain sensing in smart tire for next generation autonomous vehicles. Appl. Energy 2018, 232, 312–322. [Google Scholar] [CrossRef]
  2. Shaikh, K.; Zeadally, S. Energy harvesting in wireless sensor networks: A comprehensive review. Renew. Sustain. Energy Rev. 2016, 55, 1041–1054. [Google Scholar] [CrossRef]
  3. Dong, L.; Wen, C.; Liu, Y.; Closson, A.B.; Han, X.; Escobar, G.L.; Oglesby, M.; Feldman, M.; Chen, Z.; Zhang, J.X.J. Piezoelectric Buckled Beam Array on a Pacemaker Lead for Energy Harvesting. Adv. Mater. Technol. 2019, 4. [Google Scholar] [CrossRef] [Green Version]
  4. Balguvhar, S.; Bhalla, S. Evaluation of power extraction circuits on piezo-transducers operating under low-frequency vibration-induced strains in bridges. Strain 2019, 55, e12303. [Google Scholar] [CrossRef]
  5. Priya, S.; Inman, D.J. Energy Harvesting Technologies; Springer: Blacksburg, VA, USA, 2009. [Google Scholar]
  6. Kaźmierski, T.J.; Beeby, S. Energy Harvesting Systems, Principles, Modelling and Applications; Springer: New York, NY, USA, 2011. [Google Scholar]
  7. 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 2019, 300, 111634. [Google Scholar] [CrossRef]
  8. Tarelho, J.P.G.; Soares dos Santos, M.P.; Ferreira, J.A.F.; Ramos, A.; Kopyl, S.; Kim, S.O.; Hong, S.; Kholkin, A. Graphene-based materials and structures for energy harvesting with fluids—A review. Mater. Today 2018, 21, 1019–1041. [Google Scholar] [CrossRef]
  9. Dai, F.; Pan, D. Piezoelectric Energy Harvesting Based on Bi-Stable Composite Laminate. In Energy Harvesting; Manyala, R., Ed.; IntechOpen: London, UK, 2018; pp. 63–82. [Google Scholar]
  10. Koszewnik, A. Analytical modeling and experimental validation of an energy harvesting system for the smart plate with an integrated piezo-harvester. Sensors 2019, 19, 812. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Spreeman, D.; Manoli, Y. Electromagnetic Vibration Energy Harvesting Devices: Architectures, Design, Modelling and Optimization; Springer: Berlin, Germany, 2012. [Google Scholar]
  12. Erturk, A.; Hoffmann, J.; Inman, D.J. A piezomagnetoelastic structure for broadband vibration energy harvesting. Appl. Phys. Lett. 2009, 94, 254102. [Google Scholar] [CrossRef]
  13. Cottone, F.; Vocca, H.; Gammaitoni, L. Nonlinear energy harvesting. Phys. Rev. Lett. 2009, 102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Litak, G.; Friswell, M.I.; Adhikari, S. Magnetopiezoelastic energy harvesting driven by random excitations. Appl. Phys. Lett. 2010, 96, 214103. [Google Scholar] [CrossRef]
  15. Ando, B.; Baglio, S.; Marletta, V.; Pergolizzi, E.; Ferrari, V.; Ferrari, M.; Bulsara, A.R. Nonlinear Snap-Through-Buckling Devices for Energy Harvesting from Vibrations. Sensors 2015, 319, 409–413. [Google Scholar]
  16. Ferrari, M.; Ferrari, V.; Guizzetti, M.; Ando, B.; Baglio, S.; Trigona, C. Improved energy harvesting from wideband vibrations by nonlinear piezoelectric converters. Sens. Actuators A 2010, 162, 425–431. [Google Scholar] [CrossRef] [Green Version]
  17. Friswell, M.I.; Ali, S.F.; Adhikari, S.; Lees, A.W.; Bilgen, O.; Litak, G. Nonlinear piezoelectric vibration energy harvesting from a vertical cantilever beam with tip mass. J. Intell. Mater. Syst. Struct. 2012, 23, 1505–1521. [Google Scholar] [CrossRef]
  18. Huang, D.; Zhou, S.; Litak, G. Theoretical analysis of multi-stable energy harvesters with high-order stiffness terms. Commun. Nonlinear Sci. Numer. Simul. 2019, 69, 270–286. [Google Scholar] [CrossRef]
  19. Malaji, P.V.; Faruque Ali, S. Analysis and experiment of magneto-mechanically coupled harvesters. Mech. Syst. Signal Process. 2018, 108, 304–316. [Google Scholar] [CrossRef]
  20. Kumar, R.; Gupta, S.; Faruque Ali, S. Energy harvesting from chaos in base excited double pendulum. Mech. Syst. Signal Process. 2019, 124, 49–64. [Google Scholar] [CrossRef]
  21. Daqaq, M.F.; Masana, R.; Erturk, A.; Quinn, D.D. On the role of nonlinearities in vibratory energy harvesting: A critical review and discussion. Appl. Mech. Rev. 2014, 66. [Google Scholar] [CrossRef]
  22. Coccolo, M.; Litak, G.; Seoane, J.M.; Sanjuan, M.A.F. Optimizing the electrical power in an energy harvesting system. Int. J. Bifurc. Chaos 2015, 25, 1550171. [Google Scholar] [CrossRef]
  23. Łygas, K.; Wolszczak, P.; Litak, G.; Stączek, P. Complex response of an oscillating vertical cantilever with clearance. Meccanica 2019, 54, 1689–1702. [Google Scholar] [CrossRef] [Green Version]
  24. Wolszczak, P.; Litak, G.; Łygas, K. Analysis of dynamics of a vertical cantilever in rotary coupling to the moving frame with movement limiters. In Proceedings of the International Conference on Structural Nonlinear Dynamics and Diagnosis, Tangier, Morocco, 25–27 June 2018; Volume 241. [Google Scholar]
  25. Łygas, K.; Wolszczak, P.; Stączek, P.; Litak, G. Energy harvesting from an oscillating vertical piezoelectric cantilever with clearance. In Energy Harvesting for Wireless Sensor Networks; Kanoun, O., Ed.; De Gruyter Oldenbourg: Chemnitz, Germany, 2018; pp. 125–135. [Google Scholar]
  26. Izadgoshasb, I.; Lim, Y.Y.; Padilla, R.V.; Sedighi, M.; Novak, J.P. Performance Enhancement of a Multiresonant Piezoelectric Energy Harvester for Low Frequency Vibrations. Energies 2019, 12, 2770. [Google Scholar] [CrossRef] [Green Version]
  27. De Paula, A.S.; Savi, M.A.; Wiercigroch, M.; Pavlovskaia, E. Chaos Control Methods Applied to Avoid Bifurcations in Pendulum Dynamics; IUTAM Bookseries; IUTAM: Aberdeen, UK, 2013; Volume 32, pp. 387–395. [Google Scholar]
  28. Najdecka, A.; Kapitaniak, T.; Wiercigroch, M. Synchronous rotational motion of parametric pendulums. Int. J. Non-Linear Mech. 2013, 70, 84–94. [Google Scholar] [CrossRef]
  29. Dostal, L.; Pick, M.A. Theoretical and experimental study of a pendulum excited by random loads. Eur. J. Appl. Math. 2019, 30, 912–927. [Google Scholar] [CrossRef]
  30. Teh, S.H.; Woo, K.C.; Demrdash, H. Rotating orbits of pendulum in stochastic excitation. J. Theor. Appl. Mech. 2016, 54, 717–730. [Google Scholar] [CrossRef] [Green Version]
  31. Teh, S.H.; Woo, K.C.; Demrdash, H. Rotations of pendulum when its pivot oscillates chaotically. J. Comput. Nonlinear Dyn. 2018, 13. [Google Scholar] [CrossRef]
  32. Ding, W.; Cao, H.; Zhang, B.; Wang, K. A low frequency tunable miniature inertial pendulum energy harvester. J. Appl. Phys. 2018, 124, 164506. [Google Scholar] [CrossRef]
  33. Ylli, K.; Hoffmann, D.; Willmann, A.; Folkmer, B.; Manoli, Y. Investigation of Pendulum Structures for Rotational Energy Harvesting from Human Motion. J. Phys. Conf. Ser. 2015, 660. [Google Scholar] [CrossRef]
  34. Halim, M.A.; Rantz, R.; Zhang, Q.; Gu, L.; Yang, K.; Roundy, S. Electromagnetic energy harvesting from swing-arm motion using rotational eccentric mass structure. In Proceedings of the 2017 19th International Conference on Solid-State Sensors, Actuators and Microsystems, Kaohsiung, Taiwan, 18–22 June 2017; pp. 1863–1866. [Google Scholar]
  35. Halim, M.A.; Rantz, R.; Zhang, Q.; Gu, L.; Yang, K.; Roundy, S. An electromagnetic rotational energy harvester using sprung eccentric rotor driven by pseudo-walking motion. Appl. Energy 2018, 217, 66–74. [Google Scholar] [CrossRef]
  36. Hou, C.; Chen, T.; Li, Y.; Huang, M.; Shi, Q.; Liu, H.; Sun, L.; Lee, C. A rotational pendulum based electromagnetic/triboelectric hybrid-generator for ultra-low-frequency vibrations aiming at human motion and blue energy applications. Nano Energy 2019, 63, 103871. [Google Scholar] [CrossRef]
  37. Izadgoshasb, I.; Lim, Y.Y.; Tang, L.; Padilla, R.V.; Tang, Z.S.; Sedighi, M. Improving efficiency of piezoelectric based energy harvesting from human motions using double pendulum system. Energy Convers. Manag. 2019, 184, 559–570. [Google Scholar] [CrossRef]
  38. Payne, O.R.; Moss, S.D. Rotational energy harvesting from a drive shaft for structural health monitoring. In Proceedings of the 8th Australian Congress on Applied Mechanics, Melbourne, Australia, 23–26 November 2014; pp. 834–841. [Google Scholar]
  39. Dotti, F.E.; Reguera, F.; Machado, S.P. Rotations of the parametric pendulum excited by a reciprocating motion with a view on energy harvesting. Lect. Notes Mech. Eng. 2018, 6, 385–397. [Google Scholar]
  40. Rubes, O.; Smilek, J.; Hadas, Z. Development of vibration energy harvester fabricated by rapid prototyping technology. In Proceedings of the 16th International Conference on Mechatronics Mechatronika, Brno, Czech Republic, 3–5 December 2014. [Google Scholar]
  41. De Vlieger, J.; Vanoost, D.; Peuteman, J.; Verbrugge, S.; De Gersem, H.; Pissoort, D. Design, simulation and use of energy harvester based on permanent magnet synchronous generator. In Proceedings of the 2016 IEEE International Power Electronics and Motion Control Conference (PEMC), Varna, Bulgaria, 25–28 September 2016; pp. 557–562. [Google Scholar]
  42. Zhang, Y.; Cao, J.; Liao, W.H.; Zhao, L.; Lin, J. Theoretical modelling and experimental verification of circular Halbach electromagnetic energy harvesters for performance enhancement. Smart Mater. Struct. 2018, 27. [Google Scholar] [CrossRef]
  43. Moghaddam, R.R.; Sandulescu, P.; Tenca, P.; Zanuso, G. Efficient and accurate fully coenergy-based inductance-less FEM model of electrical machines for modern control and identification purposes. In Proceedings of the 2016 XXII International Conference on Electrical Machines (ICEM), Lausanne, Switzerland, 4–7 September 2016; pp. 299–305. [Google Scholar]
  44. Rysak, A.; Gregorczyk, M.; Chwełatiuk, K.; Gąska, D. Study of the High-Amplitude Solutions in the System of Magnetic Sliding Oscillator with many degrees of freedom. In Dynamical Systems Theory and Applications; Springer Proceedings in Mathematics & Statistics; Springer: Łódź, Poland, 2018; Volume 248, pp. 295–310. [Google Scholar]
  45. Kęcik, K. Assessment of energy harvesting and vibration mitigation of pendulum dynamic absorber. Mech. Syst. Signal Process. 2018, 106, 198–209. [Google Scholar] [CrossRef]
  46. Koszewnik, A.; Grześ, P.; Walendziuk, W. Mechanical and electrical impedance matching in a piezoelectric beam for Energy Harvesting. Eur. Phys. J. Spec. Top. 2015, 224, 2719–2731. [Google Scholar] [CrossRef]
  47. Li, Z.; Zuo, L.; Luhrs, G.; Lin, L.; Qin, Y.X. Electromagnetic energy-harvesting shock absorbers: Design, modelling, and road tests. IEEE Trans. Veh. Technol. 2013, 62, 1065–1074. [Google Scholar] [CrossRef]
  48. Kim, H.S.; Kim, J.H.; Kim, J. A review of piezoelectric energy harvesting based on vibration. Int. J. Precis. Eng. Manuf. 2011, 12, 1129–1141. [Google Scholar] [CrossRef]
  49. Masana, R.; Daqaq, M.F. Electromechanical modeling and nonlinear analysis of axially loaded energy harvesters. J. Vib. Acoust. 2011, 133. [Google Scholar] [CrossRef]
  50. Marszal, M.; Witkowski, B.; Jankowski, K.; Perlikowski, P.; Kapitaniak, T. Energy harvesting from pendulum oscillations. Int. J. Non-Linear Mech. 2017, 94, 251–256. [Google Scholar] [CrossRef]
  51. Marszal, M.; Jankowski, K.; Perlikowski, P.; Kapitaniak, T. Bifurcations of Oscillatory and Rotational Solutions of Double Pendulum with Parametric Vertical Excitation. Math. Probl. Eng. 2014. [Google Scholar] [CrossRef] [Green Version]
  52. Witkowski, B.; Perlikowski, P.; Prasad, A.; Kapitaniak, T. The dynamics of co- and counter rotating spherical pendula. Eur. Phys. J. Spec. Top. 2014, 223, 707–720. [Google Scholar] [CrossRef] [Green Version]
  53. Huguet, T.; Badel, A.; Druet, O.; Lallart, M. Drastic bandwidth enhancement of bistable energy harvesters: Study of subharmonic behaviors and their stability robustness. Appl. Energy 2018, 226, 607–617. [Google Scholar] [CrossRef]
  54. Syta, A.; Litak, G.; Friswell, M.I.; Adhikari, S. Multiple solutions and corresponding power output of a nonlinear bistable piezoelectric energy harvester. Eur. Phys. J. 2016, 89, 99. [Google Scholar] [CrossRef] [Green Version]
  55. Zhou, S.; Junyi, C.; Litak, G.; Jing, L. Numerical analysis and experimental verification of broadband tristable energy harvesters. Tm-Tech. Mess. 2018, 85, 521–532. [Google Scholar] [CrossRef]
Figure 1. (a) Solid model and real system of the electromagnetic energy harvester: (1) coil, (2) ball bearing, (3) longitudinal, cylindrical free space in form channels for additional magnets, (b) system made using 3D printing technique, (c) cross-section of system: (4) main pendulum with magnet, (5) exemplary additional magnet mounted inside coil, (d) system with additional inertial masses: (6) additional external inertial masses of the pendulum.
Figure 1. (a) Solid model and real system of the electromagnetic energy harvester: (1) coil, (2) ball bearing, (3) longitudinal, cylindrical free space in form channels for additional magnets, (b) system made using 3D printing technique, (c) cross-section of system: (4) main pendulum with magnet, (5) exemplary additional magnet mounted inside coil, (d) system with additional inertial masses: (6) additional external inertial masses of the pendulum.
Applsci 10 00671 g001
Figure 2. Principle of operation and alignment of coils and additional magnets. Interior view of the system in a plane perpendicular to the axis of rotation of the pendulum. The system could be fixed to the moving frame with an additional constant tilting angle.
Figure 2. Principle of operation and alignment of coils and additional magnets. Interior view of the system in a plane perpendicular to the axis of rotation of the pendulum. The system could be fixed to the moving frame with an additional constant tilting angle.
Applsci 10 00671 g002
Figure 3. Scheme of measuring setup.
Figure 3. Scheme of measuring setup.
Applsci 10 00671 g003
Figure 4. (a) Magnets and coil layout in Finite Element Method Magnetics (FEMM) pre-processor, (b) Exemplary results in FEMM solver.
Figure 4. (a) Magnets and coil layout in Finite Element Method Magnetics (FEMM) pre-processor, (b) Exemplary results in FEMM solver.
Applsci 10 00671 g004
Figure 5. Magnetic induction results in x-axis in additional magnet position and angular magnet position domain.
Figure 5. Magnetic induction results in x-axis in additional magnet position and angular magnet position domain.
Applsci 10 00671 g005
Figure 6. Magnetic induction in x-axis in additional magnet position domain. Three specific areas of interests were found: 1—additional magnet fully outside coil, 2—additional magnet partially inside in and outside coil, 3—additional magnet fully inside coil.
Figure 6. Magnetic induction in x-axis in additional magnet position domain. Three specific areas of interests were found: 1—additional magnet fully outside coil, 2—additional magnet partially inside in and outside coil, 3—additional magnet fully inside coil.
Applsci 10 00671 g006
Figure 7. Comparison of total magnetic forces values in the system for: (a) x = 1 mm (b) x = 13 mm, x is the distance between the additional magnet and the pendulum. Note that, the secondary peaks visible in Figure 7a originated from the finite size of magnet dipoles.
Figure 7. Comparison of total magnetic forces values in the system for: (a) x = 1 mm (b) x = 13 mm, x is the distance between the additional magnet and the pendulum. Note that, the secondary peaks visible in Figure 7a originated from the finite size of magnet dipoles.
Applsci 10 00671 g007
Figure 8. Results of (a) root mean square (RMS) voltage in the operating frequency domain for different values of additional resistance, (b) output power in the operating frequency domain for different values of additional resistance, (c) RMS voltage in the operating frequency domain for open loop, (d) output power in operating frequency domain for open loop.
Figure 8. Results of (a) root mean square (RMS) voltage in the operating frequency domain for different values of additional resistance, (b) output power in the operating frequency domain for different values of additional resistance, (c) RMS voltage in the operating frequency domain for open loop, (d) output power in operating frequency domain for open loop.
Applsci 10 00671 g008
Figure 9. Selected model features for the mechanical resonator: (a) The total energy for the Hamiltonian system (without damping and denoted in colors [J] (see the color bar); (b) Possible solutions: oscillation (blue), rotation (green) and the separatrices—orbital lines crossing the saddle points (black and red); (c) Double well potential defined with in the pendulum angular variable [−π, π]. The system parameters are indicated in Table 4.
Figure 9. Selected model features for the mechanical resonator: (a) The total energy for the Hamiltonian system (without damping and denoted in colors [J] (see the color bar); (b) Possible solutions: oscillation (blue), rotation (green) and the separatrices—orbital lines crossing the saddle points (black and red); (c) Double well potential defined with in the pendulum angular variable [−π, π]. The system parameters are indicated in Table 4.
Applsci 10 00671 g009
Figure 10. The frequency sweep results for a mechanical resonator. Ωstd and Ωmean denote variability (standard deviation) and mean value (averaged in the time interval [100 s, 200 s]) of the rotational velocity versus excitation frequency. Here the amplitude of excitation was fixed at 2 mm (a,b) and 5 cm (c,d), while the initial conditions: angle φ(t = 0) = 0, velocity φ ˙ (t = 0) = 0 (for a,c) and 80 (for b,d).
Figure 10. The frequency sweep results for a mechanical resonator. Ωstd and Ωmean denote variability (standard deviation) and mean value (averaged in the time interval [100 s, 200 s]) of the rotational velocity versus excitation frequency. Here the amplitude of excitation was fixed at 2 mm (a,b) and 5 cm (c,d), while the initial conditions: angle φ(t = 0) = 0, velocity φ ˙ (t = 0) = 0 (for a,c) and 80 (for b,d).
Applsci 10 00671 g010
Table 1. Characteristics of the vibration exciter TIRA S 513.
Table 1. Characteristics of the vibration exciter TIRA S 513.
FeatureValue
Rated peak force [N]100
Frequency range [Hz]2–7000
Max. rated travel [mm] Pk-Pk13
Max. velocity [m/s] sine1.5
Max. acceleration [g] sine45
Table 2. Characteristics of the accelerometer AltIMU-10 v3.
Table 2. Characteristics of the accelerometer AltIMU-10 v3.
FeatureValue
Operating voltage [V]2.5–5.0
Supply current [mA]6
Output formatI2C
Sensitivity range [g]±2, ±4, ±6, ±8, ±10
Table 3. Characteristics of the test.
Table 3. Characteristics of the test.
Test FeatureInput Data
Operational frequencies [Hz]10; 11; 12; 13; 14; 15
Internal resistance of the system [Ω]21.5
Additional resistance [Ω]3.3; 6.5; 12.3; 17; 21.5; 25.7; 30.3
Sampling rate per channel [S/s]100,000
Table 4. Model system parameters used for numerical simulations (as written in Equation (4)).
Table 4. Model system parameters used for numerical simulations (as written in Equation (4)).
System ParametersValues [Unit]
h0.03325 [m]
c0.000003 [N s/(kg m2)]
g9.81 [m/s2]
m0.00186 [kg]
Fm0.03 [N]
I0.0000377 [kg m2]
A0.002 or 0.05 [m]

Share and Cite

MDPI and ACS Style

Ambrożkiewicz, B.; Litak, G.; Wolszczak, P. Modelling of Electromagnetic Energy Harvester with Rotational Pendulum Using Mechanical Vibrations to Scavenge Electrical Energy. Appl. Sci. 2020, 10, 671. https://doi.org/10.3390/app10020671

AMA Style

Ambrożkiewicz B, Litak G, Wolszczak P. Modelling of Electromagnetic Energy Harvester with Rotational Pendulum Using Mechanical Vibrations to Scavenge Electrical Energy. Applied Sciences. 2020; 10(2):671. https://doi.org/10.3390/app10020671

Chicago/Turabian Style

Ambrożkiewicz, Bartłomiej, Grzegorz Litak, and Piotr Wolszczak. 2020. "Modelling of Electromagnetic Energy Harvester with Rotational Pendulum Using Mechanical Vibrations to Scavenge Electrical Energy" Applied Sciences 10, no. 2: 671. https://doi.org/10.3390/app10020671

APA Style

Ambrożkiewicz, B., Litak, G., & Wolszczak, P. (2020). Modelling of Electromagnetic Energy Harvester with Rotational Pendulum Using Mechanical Vibrations to Scavenge Electrical Energy. Applied Sciences, 10(2), 671. https://doi.org/10.3390/app10020671

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