Next Article in Journal
Fast Prototyping Microfluidics: Integrating Droplet Digital Lamp for Absolute Quantification of Cancer Biomarkers
Next Article in Special Issue
Improved Multilayered (Bi,Sc)O3-(Pb,Ti)O3 Piezoelectric Energy Harvesters Based on Impedance Matching Technique
Previous Article in Journal
An Array Switching Strategy for Direction of Arrival Estimation with Coprime Linear Array in the Presence of Mutual Coupling
Previous Article in Special Issue
A Tri-Stable Piezoelectric Vibration Energy Harvester for Composite Shape Beam: Nonlinear Modeling and Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical and Experimental Investigations of a Pseudo-Magnetic Levitation System for Energy Harvesting

Department of Applied Mechanics, Lublin University of Technology, 20-618 Lublin, Poland
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(6), 1623; https://doi.org/10.3390/s20061623
Submission received: 6 February 2020 / Revised: 4 March 2020 / Accepted: 10 March 2020 / Published: 14 March 2020
(This article belongs to the Special Issue Energy Harvesting Sensors)

Abstract

:
The paper presents an analytical, numerical and experimental analysis of the special designed system for energy harvesting. The harvester system consists of two identical magnets rigidly mounted to the tube’s end. Between them, a third magnet is free to magnetically levitate (pseudo-levitate) due to the proper magnet polarity. The behaviour of the harvester is significantly complicated by a electromechanical coupling. It causes resonance curves to have a distorted shape and a new solution from which the recovered energy is higher is observed. The Harmonic Balance Method (HBM) is used to approximately describe the response and stability of the mechanical and electrical systems. The analytical results are verified by a numerical path following (continuation) method and experiment test with use of a shaker. The influence of harvester parameters on the system response and energy recovery near a main resonance is studied in detail.

Graphical Abstract

1. Introduction

1.1. Energy Harvesting

Today, more than 40% of all energy consumption is in the form of electrical energy, which is expected to grow to 60% by 2040 [1]. Among the technologies that support the growth of clean energy is the energy harvesting (EH) process technology [2]. The EH is defined as a process where small amounts of energy are recovered from ambient energy and transformed into an electrical energy used immediately (or stored for later usage). An ambient energy is defined as many different forms energy: thermal (e.g., radiation, solar), chemical (e.g., nuclear), electrical, magnetic (e.g., magnetisation), mechanical (e.g., potential, kinetic, elastic, fluid), triboelectric [3,4] and more. The typical energy harvester device comprises one or more transducer mechanisms, power conditioning module, and energy storage (Figure 1).
The transducer mechanism converts energy from one energy type to a another, usually electricity. The power conditioning module is necessary because the output of the transducer can be intermittent or strongly disturbed, and at the wrong frequency, voltage and current to directly drive the device [5]. It rectifies the power and regulates the voltage value. The energy storage unit is needed to the energy accumulate. Usually a rechargeable battery, capacitor, or supercapacitor is used. The power harvested is large enough to power small electronic devices such as sensors [6], structural health monitoring [7,8,9,10] or wireless transmission [11]. The wireless power transfer has been proposed as a new method energy harvesting method. Liu et al. [12] proposed a 5G-based green broadband communication system with simultaneous wireless information and power transfer. A joint optimization unit has been proposed to improve the system performance. In [13], simultaneous cooperative spectrum sensing and energy harvesting model for multi-channel cognitive radio is presented. The obtained results show that this model can achieve higher throughput compared with the traditional. Moreover, the harvested energy can compensate the sensing energy loss and guarantee enough transmission power.

1.2. Pseudo-Magnetic Levitation Harvesters

The electromagnetic induction is the most cost effective type of the transducer mechanism [14]. They are based on Faraday’s law of induction—that a change in the magnetic flux of a circuit will result in the induction of an electromotive force. The electromagnetic harvesters are often modelled as a linear devices, then only perform well over a narrow bandwidth close to their natural frequency [15,16,17]. Increasing the damping of the harvester by modifying the electromechanical coupling will extend the bandwidth over which it functions effectively [17,18].
The effective improvement in the operating bandwidth can be achieved by introducing the non-linearity, tuneable resonance, and bi-stability [19,20,21]. The interesting type of the strong nonlinear electromagnetic harvesters use a pseudo-magnetic levitation effect [22,23,24] for energy recovery. Note that magnetic levitation always occurs with a help of a mechanical constraint for stability. The Earnshaw’s theorem proves that it is not possible to achieve magnetic levitation using any combination of the fixed magnets and electric charges. Therefore, the proper term should be the “pseudo-magnetic levitation”. However, in many papers, one can find the term magnetic levitation (maglev) [10,25,26,27] instead of pseudo-magnetic levitation. The pseudo-maglev levitation harvesters are characterized by simplicity of their construction, lack of springs and dampers (not a physical spring that is easily worn out), therefore the time period of usage can be a very long. Moreover, the harvester can recover energy under the external excitation in a low and high frequency ranges. As mentioned earlier, the pseudo-magnetic levitation harvesters are strongly nonlinear systems, therefore are difficult to analyse. Usually, the Finite Element Method has been used to solve the differential equations that govern the dynamics of these systems. Then the magnetic forces, magnetic field distributions are usually easy modelled [25,28]. Mann and Sims [26] presented a design for electromagnetic energy harvesting from the nonlinear oscillations of magnetic levitation. They showed the nonlinear response behaviour of the harvester to be strongly dependent on the damping level. Additionally, they showed interesting aspect that the governing equation for relative displacement reduces to the form of Duffing’s equation. This means that the potential energy function of the maglev oscillator may be expressed as: monostable hardening, monostable softening and bistable. The exemplary potential functions of the Duffing oscillator in Figure 2 are shown. The linear oscillator and nonlinear oscillator which exhibits hardening and softening nonlinearity is monostable. However, for some parameter configurations the oscillator becomes nonlinear bistable, having new stable equilibria.
Their second paper [27] presents analytical solutions for the linear response behaviour of the energy maglev harvester. Their analysis include the inductance coil, but they assumed small value of the electromechanical coupling and the resistance level. They applied the method of undetermined coefficients for analytical analysis. Berdy et al. [29] investigated the power output of a magnetic levitation vibration energy harvester on human participants while they walk and run on a treadmill. The measurement results show that the variation in power generated is relatively significant due to the variation in walking and running gait styles as well as the angle of attachment of the device.
The relevant studies that report the major achievements in the magnetic levitation harvesters in paper [30] are presented. The authors proposed four selection criteria for the harvester comparison: (1) two or more magnets (and one or more coils), (2) including at least one levitating magnet (and one or more fixed magnets), (3) executes axial motions of the levitating magnet(s) and (4) architectures with mono-stable electromagnetic-induction configurations. Moreover, the authors show a comparison of our harvester with other single levitating devices.
This paper presents complex analysis of the single pseudo-levitating harvester which is solved analytically by the HBM, numerically by the continuation method (CM) and verified by experimental shaker tests. In the literature, the numerical and experimental studies of the maglev harvesters were used. We propose the approximate analytical method to find an analytical solutions, and show analytical influence of the system parameters on energy harvesting. Additionally, the stability analysis of obtained analytical solutions is performed. All analytical results show a good compliance with the numerical research and an experiment.

2. Energy Harvesting System Design

2.1. Magnetic Levitation Architectures

The analyzed single pseudo-levitating magnet harvester plays role of the pendulum tuned mass damper. It is mounted on a special laboratory harvester–absorber system, which allows for simultaneous vibration reduction and energy recovery. More information about this can be found in [31,32].
The single levitating harvester in Figure 3a is presented. It consists of the cylindrical non-magnetic tube with the two fixed ring permanent magnets (Figure 3b).
The ring magnets have 20 mm in diameter and 10 mm in height. The levitating (moving) cylindrical permanent neodymium magnet of 20 mm in diameter and 30 mm in height is placed inside the tube between the fixed ring magnets. It motion is limited by the repulsive force exerted by the opposing magnetic field of magnets placed. The teflon surface was applied around the moving magnet to reduce friction between the moving magnet and the inner tube’s surface. Additionally, special air holes and gaps to air compression reduction are made. To prevent the magnet’s impact special rubber bumpers are installed. The induction coil is wrapped around outside of the tube, and is wound from 140 μ m diameter Copper wire. The coil has 12,740 wires, the inductance of 1.46 H and the own resistance of 1.15 k Ω . The separation distance between the moving and the bottom fixed magnet was 0.041 m. It can be changed by the screw system. This allows a modification in the magnetic suspension parameters. The levitating magnet moving axially through the center of a coil will induce a voltage across the coil terminals. The practical application of the magnet coil system can be shaker flashlights. If the flashlights vigorously shaken back and forth, then a magnet to move through a coil causing charge to the battery [33].

2.2. Mathematical Model

The dynamics of the pseudo-maglev harvester can be modeled as the nonlinear spring-mass-damper mechanical system (Figure 4a) with an external base excitation y ( t ) , connected to the electrical circuit (Figure 4b). The circuit consists of the coil with the inductance (L) and the own ( R C ) and load ( R L ) resistances. The sum of R L and R C the total resistance (R) is called. When the magnet moves inside the coil, then electromagnetic induction of voltage U ( t ) and current i ( t ) flow in the electrical circuit occurs. The electrical circuit dissipates the produced energy across a load resistor (the recovered current is converted into heat across the load resistor).
Defining a new coordinate z ( t ) = x ( t ) y ( t ) which means the relative displacement between the vibrating structure and the moving magnet, the equations of motion can be written:
m z ¨ ( t ) + c z ˙ ( t ) + k z ( t ) + k 1 z 3 ( t ) + α i ( t ) + m g = m Q ω 2 c o s ( ω t ) ,
L i ˙ ( t ) + ( R L + R C ) i ( t ) = α z ˙ ( t ) .
The pseudo-maglev suspension is modelled as a magnetic spring with the stiffness components k and k 1 (Duffing type) and the linear dashpot described by the coefficient c. A cubic nonlinearity comes from the magnetic restoring forces between the magnets [26,27,32]. The parameter m is the magnet mass, while g is the gravitational acceleration. The excitation is assumed to be harmonic y ( t ) = Q c o s ( ω t ) . The electrical and mechanical systems are coupled by the electromechanical coupling coefficient α . It characterizes how the induced voltage is related to the velocity of the magnet. The terms α i ( t ) and α z ˙ ( t ) are called electrodynamic Lorentz force ( e d f ) and the well known electromotive force ( e m f ), respectively. The fundamental relations for electromechanical systems is that e m f = e d f [26].
In the literature, different models of α can be found [25,34,35]. Generally, this parameter depends on the coil design and is a function of the magnet’s position versus the coil [18,34]. However in paper [32] authors show that fixed value of the coupling coefficient can be accepted, provided that it is properly chosen. They showed comparison between classical (fixed value) and novel (nonlinear polynomial) electromechanical coupling modelling. Moreover, in the literature the coil inductance is assumed to be very small and finally is neglected. This leads that the electrical current flow equation is reduced [36] and the electromechanical coupling coefficient is part of damping called electromagnetic.

2.3. Shaker Test

The laboratory experiments have been performed using the electromagnetic shaker system T I R A v i b 50101 which in Figure 5 is shown. The shaker was controlled by LMS.TestLab software which provided the sinusoidal input excitation signal during vibration tests and reproduced environmental conditions over a required frequency band. The two acceleration sensors for measuring the magnet’s response and to control of the excitation have been applied. Measurements of the shaker acceleration were obtained by mounting a sensor to the shaker surface. Measurements of the centre moving magnet acceleration were obtained by a second accelerometer mounted on the special element glued to the magnet.
The vibration signals are recorded in LMS SCADAS and the data is processed with LMS.Test.Lab software to obtain the acceleration and displacement spectres. The special designed harvester module consists of the DSP MicroDAQ module with multi-core application processor OMAP L137, the amplifiers, and the load resistance module which allowed to modify the resistance. The experimental results have been formed by an upward and downward frequency sweep.

3. Harmonic Balance Method (HBM)

3.1. Stability Analysis

Because the equations of motion (1) and (2) include nonlinear terms it is difficult to find their strictly correct solutions. Therefore, in the neighbourhood of main resonance the approximate solutions by the harmonic balance method (HBM) are sought. Generally, the HBM is an approximate analytical method for the study of non-linear systems described by ODE equations. The main essence is to replace the nonlinear parts the by specially-constructed functions to find approximate solutions of the nonlinear systems. To verify the accuracy of the HBM, numerical simulations have been used. The numerical simulation based on the numerical continuation technique in A u t o 07 p software [37]. This method gives a deeper understanding of the solution behaviour: stability, multiplicity and bifurcations.
The pseudo-maglev harvester response has a predominant basic frequency ω under the external periodic excitation. Therefore, the first approximate solutions are assumed as follows:
z ( t ) = z 0 ( t ) + A ( t ) s i n ( ω t ) + B ( t ) c o s ( ω t ) ,
i ( t ) = C ( t ) s i n ( ω t ) + D ( t ) c o s ( ω t ) ,
where A ( t ) , B ( t ) , C ( t ) , D ( t ) are unknown amplitudes of the harvester and amplitudes of the induced current. The amplitude z 0 is responsible for the vibration centre shift of the levitating magnet.
After substituting Equations (3) and (4) into Equations (1) and (2), and then balancing the coefficients of the corresponding s i n e and c o s i n e terms we get a set of first-order approximate differential modulation equations. For better clarity, assumed notations A ( t ) A , B ( t ) B , C ( t ) C , D ( t ) D and z 0 ( t ) z 0 .
z ˙ 0 c + m g + z 0 ( k + 3 2 A 2 k 1 + 3 2 B 2 k 1 + k 1 z 0 2 ) = 0 ,
c ( B ˙ + A ω ) + m ( 2 A ˙ ω B ω 2 ) + B k + 3 4 A 2 B k 1 + 3 4 B 3 k 1 m Q ω 2 + 3 B k 1 z 0 2 + D α = 0 ,
c ( A ˙ B ω ) m ( 2 B ˙ ω + A ω 2 ) + A k + 3 4 A 3 k 1 + 3 4 A B 2 k 1 + 3 A k 1 z 0 2 + C α = 0 ,
L ( D ˙ + C ω ) α ( B ˙ + A ω ) + D R = 0 ,
L ( C ˙ D ω ) α ( A ˙ B ω ) + C R = 0 .
Stability analysis of the harmonic solutions is carried out by using the approximate Equations (5)–(9). Determining derivatives A ˙ , B ˙ , C ˙ , D ˙ and z ˙ 0 , we get the so-called amplitude modulation equations which can be written:
z ˙ 0 = 3 k 1 A 2 z 0 + 3 k 1 B 2 z 0 + 2 k 1 z 0 3 + 2 k z 0 + 2 g m 2 c ,
A ˙ = 6 k 1 A 2 B m ω 3 k 1 A B 2 c 24 k 1 B z 0 2 m ω + 4 B c 2 ω + 8 B m 2 ω 3 8 k B m ω 6 k 1 B 3 m ω 4 ( c 2 + 4 m 2 ω 2 ) 3 k 1 A 3 c 12 k 1 A z 0 2 c 4 A c m ω 2 4 k A c 4 C α c + 8 Q m 2 ω 3 8 D α m ω 4 ( c 2 + 4 m 2 ω 2 ) ,
B ˙ = 6 k 1 A 3 m ω 3 k 1 A 2 B c + 6 k 1 A B 2 m ω + 24 k 1 A z 0 2 m ω 4 A c 2 ω 8 A m 2 ω 3 + 8 k A m ω 4 ( c 2 + 4 m 2 ω 2 ) 3 k 1 B 3 c 12 k 1 B z 0 2 c 4 B c m ω 2 4 k B c + 4 Q c m ω 2 4 D α c + 8 C α m ω 4 ( c 2 + 4 m 2 ω 2 ) ,
C ˙ = 3 k 1 A 3 α c 6 k 1 A 2 B α m ω 3 k 1 A B 2 α c 12 k 1 A z 0 2 α c 4 A α c m ω 2 4 k A α c 6 k 1 B 3 α m ω 4 L ( c 2 + 4 m 2 ω 2 ) 24 k 1 B z 0 2 α m ω 8 B α m 2 ω 3 8 k B α m ω 4 C α 2 c 8 D α 2 m ω + 8 Q α m 2 ω 3 + 4 D L c 2 ω 4 L ( c 2 + 4 m 2 ω 2 ) 4 C R c 2 + 16 D L m 2 ω 3 16 C R m 2 ω 2 4 L ( c 2 + 4 m 2 ω 2 ) ,
D ˙ = 6 k 1 A 3 α m ω 3 k 1 A 2 B α c + 6 k 1 A B 2 α m ω 4 k B α c + 8 A α m 2 ω 3 + 8 k A α m ω 3 k 1 B 3 α c 4 L ( c 2 + 4 m 2 ω 2 ) 24 k 1 A z 0 2 α m ω 12 k 1 B z 0 2 α c 4 B α c m ω 2 4 D α 2 c + 8 C α 2 m ω + 4 Q α c m ω 2 4 C L c 2 ω 4 L ( c 2 + 4 m 2 ω 2 ) 4 D R c 2 16 C L m 2 ω 3 16 D R m 2 ω 2 4 L ( c 2 + 4 m 2 ω 2 ) .
To determined stability of the obtained solutions is based on the Jacobian eigenvalues analysis, where the Jacobian takes the form
z ˙ 0 z 0 λ z ˙ 0 A z ˙ 0 B z ˙ 0 C z ˙ 0 D A ˙ z 0 A ˙ A λ A ˙ B A ˙ C A ˙ D B ˙ z 0 B ˙ A B ˙ B λ B ˙ C B ˙ D C ˙ z 0 C ˙ A C ˙ B C ˙ C λ C ˙ D D ˙ z 0 D ˙ A D ˙ B D ˙ C D ˙ D λ = 0 .
The stability of the obtained solutions depends on the roots of λ . If all roots have the negative real part then the solution is stable.

3.2. Approximate Analytical Solutions

For a steady state, amplitudes are constant thus the first-order derivatives are equal to zero
z ˙ 0 = 0 , A ˙ = 0 , B ˙ = 0 , C ˙ = 0 , D ˙ = 0 .
Introducing Equation (16) into Equations (5)–(9) and after few mathematical manipulations we obtain the following algebraic equations:
m g + z 0 ( k + 3 2 A 2 k 1 + 3 2 B 2 k 1 + k 1 z 0 2 ) = 0 ,
B k + 3 4 A 2 B k 1 + 3 4 B 3 k 1 m Q ω 2 + A c ω B m ω 2 + 3 B k 1 z 0 2 + D α = 0 ,
A k + 3 4 A 3 k 1 + 3 4 A B 2 k 1 B c ω A m ω 2 + 3 A k 1 z 0 2 + C α = 0 ,
D R + C L ω A ω α = 0 ,
C R D L ω + B ω α = 0 .
Note that Equations (20)–(21) are linear. Equations (17)–(21) have been solved with help of the Mathematica software. After a few mathematical manipulations, the amplitudes C and D are:
C = ω ( B R + A L ω ) α R 2 + L 2 ω 2 ,
D = ω ( A R + B L ω ) α R 2 + L 2 ω 2 .
Then Equations (22) and (23) are squared and added, where z 1 2 = A 2 + B 2 . This leads to the fact that current i ( i 2 = C 2 + D 2 ) and the recovered power P = R i 2 responses are written in the simple form:
i = ω α z 1 R 2 + L 2 ω 2 , P = ω 2 α 2 z 1 2 R R 2 + L 2 ω 2 .
Introducing amplitudes C and D (Equations (22) and (23)) into Equations (18) and (19), we get the amplitudes A and B:
A = m Q ω 3 ( c R α 2 R 2 + L 2 ω 2 ) ( c R 2 ω + c L 2 ω 3 + R ω α 2 ) 2 ( R 2 + L 2 ω 2 ) 2 + ( k + 3 4 k 1 ( z 1 2 + 4 z 0 2 ) + ω 2 ( m + L α 2 R 2 + L 2 ω 2 ) ) 2 ,
B = m Q ω 2 ( k + 3 4 k 1 ( z 1 2 + 4 z 0 2 ) + ω 2 ( m + L α 2 R 2 + L 2 ω 2 ) ) ( c R 2 ω + c L 2 ω 3 + R ω α 2 ) 2 ( R 2 + L 2 ω 2 ) 2 + ( k + 3 4 k 1 ( z 1 2 + 4 z 0 2 ) + ω 2 ( m + L α 2 R 2 + L 2 ω 2 ) ) 2 .
Finally, the sixth-order non-linear polynomial equation describing the magnet’s amplitude and the magnet vibration centre shift are obtained:
[ 9 k 1 2 ( R 2 + L 2 ω 2 ) ] z 1 6 + 24 k 1 ( R 2 + L 2 ω 2 ) [ k m ω 2 + 3 k 1 z 0 2 + L ω 2 α 2 R 2 + L 2 ω 2 ] z 1 4 + + 16 ( R 2 + L 2 ω 2 ) [ k 2 + c 2 ω 2 2 k m ω 2 + m 2 ω 4 + 6 k k 1 z 0 2 6 k 1 m ω 2 z 0 2 + 9 k 1 2 z 0 4 ] z 1 2 + + 16 ω 2 α 2 [ 2 k L + 2 c R 2 L m ω 2 + 6 k 1 L z 0 2 + α 2 ] z 1 2 16 m 2 Q 2 ω 6 ( R 2 + L 2 ω 2 ) = 0 ,
m g + z 0 ( k + 3 2 k 1 z 1 2 + k 1 z 0 2 ) = 0 .
The frequency–amplitude harvester response is given by solving of Equations (27) and (28).

4. Methods

In order to check the correctness of the HBM method, the obtained results were compared with the numerical continuation method (CM) and the experiment. The CM is a very good numerical method suited for tracing one-dimensional manifolds, curves (called branches) of solutions. It allows to detect the bifurcation points, new solution branches and calculate the stable and unstable solutions. The theoretical analysis is based on the parameters from the experimental laboratory rig (Figure 5). Some of these are readily measurable: m = 0.09 kg, L = 1.46 H, R L = 0 ÷ 10 k Ω , R C = 1.15 k Ω and Q = 0.6 mm. The damping coefficient is estimated from the logarithmic decrement method equals c = 0.07 Ns/m. In order to determine the magnetic suspension parameters, the force–displacement relationship was conducted with use of Shimadzu tensile tests. Of course, the linear stiffness depends on the distance between the fixed and moving magnets [34]. In our case, these parameters were k = 640 N/m and k 1 = 460 kN/m 3 . The value of k 1 depends on geometry and size of the all magnets. The distance between the moving magnet (with sensor) and the shaker was 0.15 m. It has been estimated based on our own tests and suggestions from [38].

5. Results and Discussion

5.1. Single Pseudo-Levitating Magnet Vibration Centre Shift

The magnet-coil configuration influences the electromechanical coupling [32]. Therefore, the magnet vibration centre must be taken into consideration. The detail analysis of the Equation (28) shows that the magnet vibration centre z 0 depends on the magnet oscillation amplitude z 1 .
In Figure 6a, we can see that the magnet amplitude z 1 influences the vibration center shift. At high enough oscillation amplitudes, the moving magnet vibration centre shift is practically independent from the linear stiffness (Figure 6a).
Moreover, it can be seen that for z 1 = 0 we observe the static displacement of the moving magnet. Of course, an increase of the linear stiffness component k reduces the static displacement of magnet.
Figure 6b depicts the magnet amplitude influencing the vibration centre shift. The analysis was done for three different nonlinear stiffness: k 1 = 100 kN/m 3 (red line), k 1 = 460 kN/m 3 (black line) and k 1 = 1000 kN/m 3 (blue line). For the low magnet oscillation, the influence of the nonlinearity on the vibration centre is weak. An increase of the oscillation amplitude causes that the vibration centre is reduced.

5.2. Analytical Results

In order to show the influence of the amplitude excitation (Q), a series of analyses have been performed. The exemplary resonance curves of the moving magnet and the recovered power versus the frequency ω in Figure 7 are shown. The blue line corresponds with the case where Q = 0.1 mm, the red line to Q = 0.5 mm, the green line to Q = 0.8 mm and the black line to Q = 1.1 mm.
The most effective energy harvesting occurs in the resonance conditions. The resonance peak is located close to the frequency ω 88 rad/s. Of course, a greater magnet oscillation amplitude Q causes an increase in the magnet response and leads to higher level of the recovered power. For example, an increase of the excitation amplitude about two times (from 0.6 mm to 1.1 mm) causes a two times increase in the magnet vibration and a four times increase in the recovered power. Note that all resonance curves are stable and have the linear behaviour.
The frequency response curves of the pseudo-maglev harvester under the total resistance R in Figure 8 are shown. The analytical study has been done for the resistances: 0.5 k Ω (blue line), 1 k Ω (red line), 2 k Ω (green line) and 5 k Ω (black line). Upon analysing and comparing the obtained results from both diagrams, we can conclude that the resistance strongly influences the magnet dynamics and energy harvesting.
An increase of R causes the magnet oscillation to become higher. This effect can be explained by analysis of Equations (1) and (2), and neglecting the inductance (L = 0 H). Then, e d f force (Equation (1)) can be written as e d f = α 2 z ˙ R . This means that the higher resistance decreases the electrical damping level ( e d f force). However, the EH under R looks quite different. A low value of R causes the resonance peak not to be observed. An increase of R causes higher magnet oscillation but reduces the effective bandwidth of the power (Figure 8b).
The parameter L (inductance) characterizes behaviour of the coil. It is defined in terms of that opposing e m f force or its generated magnetic flux and the corresponding electric current. The coil inductance depends on the geometry of the current path as well as the magnetic permeability of nearby materials. As mentioned earlier, much of the literature suggests that the coil inductance does not influence the frequency response as well as energy harvesting. This is true only for the low inductance (usually typical in real practice). Our analytical results (Figure 9a,b) show that the inductance lower than 1H can be neglected.
A main question seems to be when the L can be ignored. In this aim, the relationship between the moving magnet response (Figure 10a) and the recovered current (Figure 10b) versus the coil inductance has been plotted. For the analysis, three different resistance levels: R = 1 k Ω (red line), R = 2 k Ω (black line) and R = 5 k Ω (blue line) have been selected. Assuming, that the neglect of L causes a 2% error in the magnet’s oscillation, then the maximal coil inductance can be: 1.28 H for R = 1 k Ω , 1.69 H for R = 2 k Ω and 2.57 H for R = 5 k Ω (Figure 10a). Whereas, for a 2% error in the recovered current, the inductance can range: 1.84 H for R = 1 k Ω , 2.07 H for R = 2 k Ω and 2.25 H for R = 5 k Ω (Figure 10b). The analysis clearly shows that the coil inductance can be ignored, but this depends on the total resistance level.
The last analysed parameter is the electromechanical coupling coefficient ( α ). This parameter characterizes how the induced current in the coil is related to the velocity of the moving magnet [32]. As shown in the amplitude-frequency, plot this parameter is a crucial from the dynamics (Figure 11a) and energy harvesting (Figure 11b) point of view. We can see that the electromechanical coupling coefficient strongly affects the magnet amplitude responses. If the electromechanical coefficient increases, then the magnet’s amplitude and the recovered energy are significantly reduced. This means that α can be treated as the electrical damping coefficient (increases e d f force in Equation (1)).
Additionally, the jump amplitude phenomena (well-known as the foldover effect) close to the main resonance is observed. For higher values of α the jump phenomena gradually disappears. The jump phenomena causes coexistence of two stable and one unstable solutions and is positive for EH. One of the responses is characterized by the higher energy input (from top branch). For example, for α = 12 V/A the maximal recovered power from top branch equals about P = 2.7 W while from the bottom branch is P = 0.03 W (Figure 11b). Moreover, the foldover effect causes a broadening effective frequency bandwidth. For the coupling value of α = 20 V/A the effective frequency bandwidth is located ω 82 ÷ 92 rad/s, but for α = 12 V/A the region became wider to ω 85 ÷ 125 rad/s. Of course, the higher EH is strongly related with increasing of the magnet oscillation (Figure 11a).

5.3. Experimental Verification

A series of experimental tests for the two different resistances have been performed: R = 2 k Ω and R = 10 k Ω . The exemplary resonance curves of the moving magnet response ( z 1 ) and the recovered current (i) in Figure 12 and Figure 13 are shown respectively. The experiment tests are marked by the red plus signs. The blue line means the numerical stable solution, while the blue circle points denote the stable solution obtained by the HBM.
The HBM and CM results are compliance in the total frequency range. This means that the first approximation of the assumed solutions is sufficient. The experimental results are similar to the analytical and numerical. A few small differences probably come from the magnet-tube friction. The maximal recovered current has been obtained close to the resonance peak ω = 88 rad/s equals i = 5.8 mA (what correspond power P = 0.067 W).
Modification in the coupling coefficient and the resistance values cause the resonance curves to have a distorted shape and the foldover effect to appear (Figure 13a,b).
The blue lines mean the stable numerical solution, the blue circle points denote the stable solution obtained by the HBM, the red lines mean the unstable numerical solution, and the red point show the unstable analytical solution. The experimental verification marked again by the red plus signs. The analytical and numerical results are in compliance. The foldover effect exhibits two stable and one unstable solutions has been confirmed experimentally. The numerical recovered power from the top branch equals P = 3.6 W, while about P = 0.036 W from the bottom branch.
The exemplary experimental and numerical recovered time histories in Figure 14 are shown. We can see a current flow in the electrical circuit for two different electrical parameters configurations.
The red color means the recovered current from the pseudo-maglev harvester, while the black line means numerical result. The numerical time histories are clear periodic signals, while the experimental signals are slightly disturbed probably by the nonlinear friction between the levitating magnet and tube, the model simplification and parameters identification. However, both signals have similar periods and amplitudes.

6. Conclusions

The paper presents theoretical and experimental analysis of the single pseudo-levitating harvester. The HBM method is successfully applied to obtain the analytical responses and calculation of the stability solutions. The analytical, numerical and experimental results are in compliance already for the first assumed approximation solutions. Moreover, the influence of the electrical parameters on energy harvesting is investigated in detail.
The obtained results show that the magnet vibration centre shift depends on the magnet’s oscillation amplitude. This is very a important conclusion because the magnet shift relative to the coil changes the electromechanical coupling value.
The frequency response strongly depends on the resistance and electromechanical coupling coefficient. An increase in the resistance causes the higher oscillation of the levitating magnet. Moreover, the coupling coefficient strongly influences the resonance curve shape. In low values of α , the strong nonlinear behaviour is observed with two stable periodic solutions. The obtained results confirm that the coil inductance can be neglected; however, this depends on the resistance level. For the analyzed parameters, the coil inductance can be neglected if is smaller then 1.28 H. The maximum recovered current from numerical tests was 3.6 W, while from the experiment it was about 1 W.
The obtained results show that for the tested device the recovered power reaches few watts. It is expected that in larger structures, the harvested energy will get higher amounts. The effectiveness of the systems is relatively small—about 2% of input. However, the proposed device is a part of the pendulum absorber/harvester system. Therefore, the recovered current mainly will come from the vibration mitigation.
A more detailed experimental investigation and the shaping of the electromechanical coupling by the special designed magnet-separators stacks and the energy leakage problem will be carried out in a future work.

Author Contributions

Conceptualization, methodology, software, experimental and analytical investigations, writing—original draft and final preparation, numerical simulation, HBM and data visualization K.K.; stability analysis and experiment A.M. All authors have read and agree to the published version of the manuscript.

Funding

The project/research was financed in the framework of the project Lublin University of Technology—Regional Excellence Initiative, funded by the Polish Ministry of Science and Higher Education (contract no. 030/RID/2018/19).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
HBMharmonic balance method
EHenergy harvesting
maglevmagnetic levitation
CMcontinuation method
emfelectromotive force
edfelectrodynamic (Lorentz) force

References

  1. European Center for Power Electronics (ECPE). Position Paper on Energy Efficiency—The Role of Power Electronics. March 2007. Available online: https://www.ecpe.org (accessed on 20 January 2020).
  2. Kilner, J.A.; Skinner, S.J.; Irvine, S.J.C.; Edwards, P.P. (Eds.) Functional Materials for Sustainable Energy Applications; Woodhead Publishing: Cambridge, UK, 2012. [Google Scholar]
  3. Slabov, V.; Kopyl, S.; Dos Santos, M.P.S.; Kholkin, A.L. Natural and eco-friendly materials for triboelectric energy harvesting. Nano-Micro Lett. 2020, 12, 1–18. [Google Scholar] [CrossRef] [Green Version]
  4. Tarelho, J.P.G.; Dos Santos, M.P.S.; 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]
  5. Spreemann, D.; Manoli, Y. (Eds.) Electromagnetic Vibration Energy Harvesting Devices: Architectures, Designs, Modeling and Optimization, Advanced Microelectronics; Springer: New York, NY, USA, 2012. [Google Scholar]
  6. Shaikh, F.K.; Zeadally, S. Energy harvesting in wireless sensor networks: A comprehensive review. Renew. Sustain. Energy Rev. 2016, 55, 1041–1054. [Google Scholar] [CrossRef]
  7. Qing, X.; Li, W.; Wang, Y.; Sun, H. Piezoelectric transducer-based structural health monitoring for aircraft in applications. Sensors 2019, 19, 545. [Google Scholar] [CrossRef] [PubMed]
  8. Korosak, Z.; Suhadolnik, N.; Pletersek, A. The implementation of a low power environmental monitoring and soil moisture measurement system based on UHF RFID. Sensors 2019, 19, 5527. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Bernardo, R.; Rodrigues, A.; Dos Santos, M.P.S.; Carneiro, P.; Lopes, A.; Amaral, J.S.; Amaral, V.S.; Morais, R. Novel magnetic stimulation methodology for low-current implantable medical devices. Med. Eng. Phys. 2019, 73, 77–84. [Google Scholar] [CrossRef] [PubMed]
  10. Dos Santos, M.P.; Ferreira, J.A.F.; Simões, J.A.O.; Pascoal, R.; Torrao, J.; Xue, X.; Furlani, E.P. Magnetic levitation-based electromagnetic energy harvesting: A semi-analytical non-linear model for energy transduction. Sci. Rep. 2016, 6, 18579. [Google Scholar] [CrossRef]
  11. Kaleta, J.; Mech, R.; Wiewiorski, P. Energy harvester based on magnetomechanical effect as a power source for multi-node wireless network. Small-scale energy harvesting. In A Guide to Small-Scale Energy Harvesting Techniques; Lallart, M., Ed.; IntechOpen: London, UK, 2012. [Google Scholar]
  12. Liu, X.; Zhang, X.; Jia, M.; Fan, L.; Lu, W.; Zhai, Z. 5G-based green broadband communication system design with simultaneous wireless information and power transfer. Phys. Commun. 2018, 28, 130–137. [Google Scholar] [CrossRef]
  13. Liu, X.; Li, F.; Na, Z. Optimal resource allocation in simultaneous cooperative spectrum sensing and energy harvesting for multichannel cognitive radio. IEEE Access 2017, 5, 3801–3812. [Google Scholar] [CrossRef]
  14. Cepnik, C.; Lausecker, R.; Wallrabe, U. Review on electrodynamic energy harvestersa classification approach. Micromachines 2013, 5, 168–196. [Google Scholar] [CrossRef] [Green Version]
  15. Williams, C.B.; Yates, R.B. Analysis of a micro-electric generator for microsystems. Sens. Actuators A Phys. 1996, 52, 8–11. [Google Scholar] [CrossRef]
  16. Mitcheson, P.D.; Green, T.C.; Yeatman, E.M.; Holmes, A.S. Architectures for vibration-driven micropower generators. J. Miecroelectromech. Syst. 2004, 13, 429–440. [Google Scholar] [CrossRef] [Green Version]
  17. Stephen, N.G. On energy harvesting from ambient vibration. J. Sound Vib. 2006, 293, 409–425. [Google Scholar] [CrossRef] [Green Version]
  18. Beeby, S.P.; Tudor, M.J.; White, N.M. Energy harvesting vibration sources for microsystems applications. Meas. Sci. Technol. 2006, 17, R175. [Google Scholar] [CrossRef]
  19. Jiang, W.A.; Chen, L.Q. Energy harvesting of monostable Duffing oscillator under Gaussian white noise excitation. Mech. Res. Commun. 2013, 53, 85–91. [Google Scholar] [CrossRef]
  20. Zhou, S.; Cao, J.; Inman, D.J.; Lin, J.; Liu, S.; Wang, Z. Broadband tri-stable energy harvester: Modelling and experiment verification. Appl. Energy 2014, 133, 33–39. [Google Scholar] [CrossRef]
  21. Jiang, W.A.; Chen, L.Q. A piezoelectric energy harvester based on internal resonance. Acta Mech. Sin. 2015, 31, 223–228. [Google Scholar]
  22. Sun, K.; Liu, G.Q.; Xu, X.Y. A piezoelectric energy harvester based on internal resonance. Appl. Mech. Mater. 2011, 128–129, 923–927. [Google Scholar] [CrossRef]
  23. Saha, C.; O’Donnell, T.; Wang, N.; McCloskey, P. Electromagnetic generator for harvesting energy from human motion. Sens. Actuators A Phys. 2008, 147, 248–253. [Google Scholar] [CrossRef]
  24. Dallago, E.; Marchesi, M.; Venchi, G. Analytical model of a vibrating electromagnetic harvester considering nonlinear effects. IEEE Trans. Power Electron. 2010, 25, 1989–1997. [Google Scholar] [CrossRef]
  25. Bernal, A.; Garcia, L. The modelling of an electromagnetic energy harvesting architecture. Appl. Math. Model. 2012, 36, 4728–4741. [Google Scholar] [CrossRef]
  26. Mann, B.; Sims, N. Energy harvesting from the nonlinear oscillations of magnetic levitation. J. Sound Vib. 2009, 319, 515–530. [Google Scholar] [CrossRef] [Green Version]
  27. Mann, B.; Sims, N. On the performance and resonant frequency of electromagnetic induction energy harvesters. J. Sound Vib. 2010, 329, 1348–1361. [Google Scholar] [CrossRef] [Green Version]
  28. Dos Santos, M.P.S.; Coutinho, J.; Marote, A.; Sousa, B.; Ramos, A.; Ferreira, J.A.F.; Bernardo, R.; Rodrigues, A.; Marques, A.T.; Da Cruz e Silva, O.A.B.; et al. Capacitive technologies for highly controlled and personalized electrical stimulation by implantable biomedical systems. Sci. Rep. 2019, 9, 5001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Berdy, D.F.; Valentino, D.J.; Peroulis, D. Kinetic energy harvesting from human walking and running using a magnetic levitation energy harvester. Sens. Actuators A Phys. 2015, 222, 262–271. [Google Scholar] [CrossRef]
  30. Carneiro, P.; Dos Santos, M.P.S.; Rodrigues, A.; Ferreira, J.A.F.; Simoes, J.A.O.; Marques, A.T.; Kholkin, A.L. Electromagnetic energy harvesting using magnetic levitation architectures: A review. Sens. Actuators A Phys. 2020, 260, 114191. [Google Scholar] [CrossRef]
  31. Kecik, K. Assessment of energy harvesting and vibration mitigation of a pendulum dynamic absorber. Mech. Syst. Signal Process. 2018, 106, 198–209. [Google Scholar] [CrossRef]
  32. Kecik, K.; Mitura, A.; Lenci, S.; Warminski, J. Energy harvesting from a magnetic levitation system. Int. J. Nonliner Mech. 2017, 94, 200–206. [Google Scholar] [CrossRef]
  33. Gieras, J.; Kucharski, A.; Piechowski, J. Performance characteristics of a shake flashlight. In Proceedings of the 2017 IEEE International Electric Machines and Drives Conference (IEMDC), Miami, FL, USA, 21–24 May 2017. [Google Scholar]
  34. Kecik, K. Architecture and optimization of a low-frequency maglev energy harvester. Int. J. Struct. Stab. Dyn. 2019, 8, 1–13. [Google Scholar] [CrossRef]
  35. Munaz, A.; Lee, B.; Chung, G. A study of an electromagnetic energy harvester using multi-pole magnet. Sens. Actuators A Phys. 2013, 201, 134–140. [Google Scholar] [CrossRef]
  36. Yurchenko, D.; Alevras, P. Parametric pendulum based wave energy converter. Mech. Syst. Signal Process 2018, 99, 504–515. [Google Scholar] [CrossRef] [Green Version]
  37. Doedel, E.; Oldeman, B. Auto-07p: Continuation and Bifurcation Software; Concordia Univ. Canada: Montreal, QC, Canada, 2012; Available online: http://indy.cs.concordia.ca/auto/ (accessed on 20 January 2020).
  38. Medina, N.; Vicente, J.; Robles, J. Magnetic fields influence on sensors with electrical output under sinusoidal excitations. Acta IMEKO 2016, 5, 51–58. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Typical energy harvesting process.
Figure 1. Typical energy harvesting process.
Sensors 20 01623 g001
Figure 2. The potential energy function of the Duffing oscillator.
Figure 2. The potential energy function of the Duffing oscillator.
Sensors 20 01623 g002
Figure 3. Photo (a) and scheme of the single pseudo-levitating harvester (b). This system was originally designed and applied as the pendulum tuned mass damper (see paper [31]).
Figure 3. Photo (a) and scheme of the single pseudo-levitating harvester (b). This system was originally designed and applied as the pendulum tuned mass damper (see paper [31]).
Sensors 20 01623 g003
Figure 4. Physical model of pseudo-maglev harvester (a) and the electrical circuit comprising of the coil and the external load resistor (b). The magnet oscillation causes the current induction in the electrical circuit.
Figure 4. Physical model of pseudo-maglev harvester (a) and the electrical circuit comprising of the coil and the external load resistor (b). The magnet oscillation causes the current induction in the electrical circuit.
Sensors 20 01623 g004
Figure 5. Measuring stand: shaker TIRA vib 50101 with the harvester (a), the LMS SCADAS module (b) with the analogue TIRA amplifiers (c), the load resistance module (d), the harvester DSP module (e), the harvester amplifiers (f) and the PC with the TestLab software (g).
Figure 5. Measuring stand: shaker TIRA vib 50101 with the harvester (a), the LMS SCADAS module (b) with the analogue TIRA amplifiers (c), the load resistance module (d), the harvester DSP module (e), the harvester amplifiers (f) and the PC with the TestLab software (g).
Sensors 20 01623 g005
Figure 6. An influence of the linear (a) and nonlinear (b) stiffness components on the magnet vibration centre. The diagrams have been obtained for parameters k 1 = 460 kN/m 3 and k 1 for k = 640 N/m from the Equation (28).
Figure 6. An influence of the linear (a) and nonlinear (b) stiffness components on the magnet vibration centre. The diagrams have been obtained for parameters k 1 = 460 kN/m 3 and k 1 for k = 640 N/m from the Equation (28).
Sensors 20 01623 g006
Figure 7. Analytical HBM results. The moving magnet’s resonance curves (a) and the recovered power (b) under the external periodic excitation. The results have been calculated for the resistance R = 2.0 k Ω . The results have obtained from Equations (27) and (28). All resonance responses are stable.
Figure 7. Analytical HBM results. The moving magnet’s resonance curves (a) and the recovered power (b) under the external periodic excitation. The results have been calculated for the resistance R = 2.0 k Ω . The results have obtained from Equations (27) and (28). All resonance responses are stable.
Sensors 20 01623 g007
Figure 8. Analytical HBM results. Magnet’s resonance curves (a) and recovered power (b) under the resistance influence. The results have been calculated for amplitude Q = 0.6 mm and have obtained from Equations (27) and (28). All resonance curves are stable.
Figure 8. Analytical HBM results. Magnet’s resonance curves (a) and recovered power (b) under the resistance influence. The results have been calculated for amplitude Q = 0.6 mm and have obtained from Equations (27) and (28). All resonance curves are stable.
Sensors 20 01623 g008
Figure 9. Analytical HBM results. Magnet’s resonance curves (a) and the recovered power (b) under the coil inductance influence. The results have been calculated for amplitude Q = 0.6 mm and the resistance R = 2 k Ω . The results have been obtained from Equations (27) and (28). All resonance curves are stable.
Figure 9. Analytical HBM results. Magnet’s resonance curves (a) and the recovered power (b) under the coil inductance influence. The results have been calculated for amplitude Q = 0.6 mm and the resistance R = 2 k Ω . The results have been obtained from Equations (27) and (28). All resonance curves are stable.
Sensors 20 01623 g009
Figure 10. Influence of the coil inductance L on the magnet’s response (a) and the recovered current (b). The results have been obtained for frequency ω = 88 rad/s.
Figure 10. Influence of the coil inductance L on the magnet’s response (a) and the recovered current (b). The results have been obtained for frequency ω = 88 rad/s.
Sensors 20 01623 g010
Figure 11. Analytical HBM results. Magnet’s resonance curves (a) and the recovered power (b) under the electromechanical coupling coefficient. The results have been calculated for Q = 0.6 mm and R = 2 k Ω . The results have been obtained from Equations (27) and (28). The dashed line means unstable solution obtained from HBM.
Figure 11. Analytical HBM results. Magnet’s resonance curves (a) and the recovered power (b) under the electromechanical coupling coefficient. The results have been calculated for Q = 0.6 mm and R = 2 k Ω . The results have been obtained from Equations (27) and (28). The dashed line means unstable solution obtained from HBM.
Sensors 20 01623 g011
Figure 12. Magnet’s resonance curve (a) and the recovered current (b) for R = 2 k Ω and α = 70 V/A. The marked black points denote the HBM results, the blue lines mean numerical results and the plus signs are experimental results.
Figure 12. Magnet’s resonance curve (a) and the recovered current (b) for R = 2 k Ω and α = 70 V/A. The marked black points denote the HBM results, the blue lines mean numerical results and the plus signs are experimental results.
Sensors 20 01623 g012
Figure 13. Magnet’s resonance curve (a) and the recovered current (b) for R = 10 k Ω and α = 25 Vs/m. The marked black points denote the HBM results, the blue and red lines mean numerical results and the plus signs are experimental results. The foldover effect toward high frequencies is observed. This is commonly named hardening.
Figure 13. Magnet’s resonance curve (a) and the recovered current (b) for R = 10 k Ω and α = 25 Vs/m. The marked black points denote the HBM results, the blue and red lines mean numerical results and the plus signs are experimental results. The foldover effect toward high frequencies is observed. This is commonly named hardening.
Sensors 20 01623 g013
Figure 14. Experimental and numerical recovered currents obtained for parameters: α = 70 V/A, R = 2 k Ω (a) and α = 25 V/A, R = 10 k Ω (b). The frequency and amplitude of excitation were ω = 88 rad/s and Q = 0.6 m.
Figure 14. Experimental and numerical recovered currents obtained for parameters: α = 70 V/A, R = 2 k Ω (a) and α = 25 V/A, R = 10 k Ω (b). The frequency and amplitude of excitation were ω = 88 rad/s and Q = 0.6 m.
Sensors 20 01623 g014

Share and Cite

MDPI and ACS Style

Kecik, K.; Mitura, A. Theoretical and Experimental Investigations of a Pseudo-Magnetic Levitation System for Energy Harvesting. Sensors 2020, 20, 1623. https://doi.org/10.3390/s20061623

AMA Style

Kecik K, Mitura A. Theoretical and Experimental Investigations of a Pseudo-Magnetic Levitation System for Energy Harvesting. Sensors. 2020; 20(6):1623. https://doi.org/10.3390/s20061623

Chicago/Turabian Style

Kecik, Krzysztof, and Andrzej Mitura. 2020. "Theoretical and Experimental Investigations of a Pseudo-Magnetic Levitation System for Energy Harvesting" Sensors 20, no. 6: 1623. https://doi.org/10.3390/s20061623

APA Style

Kecik, K., & Mitura, A. (2020). Theoretical and Experimental Investigations of a Pseudo-Magnetic Levitation System for Energy Harvesting. Sensors, 20(6), 1623. https://doi.org/10.3390/s20061623

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