Next Article in Journal
Numerical Simulation and Flow Field Analysis of Porous Water Jet Nozzle Based on Fluent
Previous Article in Journal
Classification of Rock Mass Quality in Underground Rock Engineering with Incomplete Data Using XGBoost Model and Zebra Optimization Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Vibrational Energy Harvesting Sensor Based on Linear and Rotational Electromechanical Effects

Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, CA 94550, USA
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(16), 7072; https://doi.org/10.3390/app14167072
Submission received: 11 July 2024 / Revised: 31 July 2024 / Accepted: 6 August 2024 / Published: 12 August 2024
(This article belongs to the Topic Advanced Energy Harvesting Technology)

Abstract

:
In this investigation, a magnetically coupled double-spring design is presented for harvesting low-level non-stationary random vibrational energy. The sensor relies on multimodal coupling between the translation and rotation of a two-spring magnet and coil system to widen the harvesting bandwidth. Energy methods are used to develop a model to characterize the electromechanical response of the system, the solution of which is obtained using stochastic techniques based on a particle swarm algorithm. This approach provides an efficient method to estimate system parameters that otherwise are difficult or impossible to determine with independent measurements. The experimental results demonstrate agreement with the theoretical predictions over a limited bandwidth. The sensor can effectively harvest non-stationary vibration energy down to 10−4 g within a limited bandwidth of 130–150 Hz . The sensor prototype has an operational volume of 2.6 c m 3 with a calculated power density of 0.2 W / c m 3 . The sensor’s small size results in a coupling efficiency of approximately 6% across the tested bandwidth.

1. Introduction

Self-powered sensors and networks have become a tangible reality, thanks to the ability to harness ambient energy from the surroundings to power electronic devices. This advancement has a profound impact on the realization of intelligent and adaptable sensor networks. As microelectronics and sensing technologies continue to advance, the power requirements for these devices have reduced to the microwatt regime. At this level, the available ambient energy in the environment becomes an attractive resource to power these systems. Efficient energy harvesting technologies will greatly diminish the reliance on short-lived, battery-powered approaches and eliminate the need for direct power supply to the devices [1,2,3]. This paves the way for the development of truly autonomous wireless sensor networks that can operate indefinitely, independent of routine maintenance and infrastructure requirements. These networks will have diverse applications in structural health monitoring, security perimeter tamper detection, autonomous systems, and seamless integration with existing Internet of Things (IoT) technologies [4,5].
A key challenge for persistent sensing is power supply. Even when power requirements are in the microwatt range, they need to be continuous. Many research efforts are focused on developing sensor technologies that leverage ambient energy sources. For systems operating on machinery or structures with well-defined and repeatable dynamic characteristics (e.g., rotating machinery), the energy harvesting mechanism is designed to target specific vibrational frequencies [6]. This approach achieves optimal performance when the harvesting design resonates at or near one of the fundamental operating frequencies of the monitored system. Even small levels of vibrational energy can contribute significantly to harvested power when the harvesting system operates near a resonance. However, linear systems face performance degradation as the operational frequency deviates from the harvester’s resonance design point [7,8,9]. Single-resonance VEH devices operating in a linear regime typically have an FWHM resonance that is 2–3 Hz . Nonlinear effects have shown an improvement by a factor 5, or 10–15 Hz [7,8,9].
One approach to introduce nonlinearities is through opposing magnetic fields [10,11]. Muskat provides a thorough review of many electromagnetic vibrational energy harvesting (VEH) designs [12]. In particular, designs utilizing multiple masses and magnets to generate favorable harvesting performance have been the focus of numerous investigations [4,13,14,15]. These designs are based on modifying the electromagnetic coupling of the magnet and the coil with the introduction of additional magnets to generate nonlinear forcing functions that lead to local minima in the potential energy function. The local minima generate a multi-stable vibrational response that leads to improved bandwidth and conversion capabilities. The approach is not without challenges—large potential wells can “trap” the system, reducing the high-energy interwell oscillations. To mitigate these effects, several investigators have developed tri-stable approaches resulting in lower potential interwell boundaries, resulting in lower excitation energies to achieve the interwell boost for improved energy harvesting [16,17].
In the following sections, we explore the use of nonlinear coupled multi-modal behavior to enhance bandwidth and increase energy conversion efficiency. Our focus is on developing simplified ordinary differential equation (ODE) approximations to capture the electromechanical behavior. Experimental measurements and comparisons to the ODE model are used to validate the effectiveness of the sensor design. The nonlinear ODE system is also used to explore parameter values for analyzing performance and “tuning” of the device. Here, model parameters can be fit to match empirical data.
The device is described in Section 2. The state-space ODE theory, derived in Section 3, is based upon a Lagrangian formulation. Experience has shown that when fitting the model parameters to match its output to data, the cost function has multiple local minima. Thus, a particle swarm optimizer was used to fit the model parameters. This is described in Section 4. Section 5 describes the experimental design and testing method. The results and harvesting performance are presented in Section 6. The results are discussed in Section 7 with conclusions presented in Section 8.

2. Device Description

The vibrational energy harvesting (VEH) device is pictured in Figure 1. An idealized diagram of the sensor is shown in Figure 2. The prototype sensor uses dual cylindrical magnets each attached to its own planar spring. Threaded rods are utilized to adjust the separation distance of the magnets and their angular orientation. One magnet is at the center of a cylindrically wrapped copper coil. The magnets are oriented with similar poles counterposed. The repulsive magnetic coupling introduces an additional magnetic force in addition to the back electromotive force (EMF) resulting from the coil/magnet interaction. The effect introduces additional nonlinear force behavior to the response of the VEH. Both magnets experience the same base excitation. Another characteristic of the design which results in rotational motion can be attributed to the intentional asymmetries caused by the threaded rods adjustments and unintentional asymmetries due to device imperfections [17].
A result of the asymmetrical repulsive force is to induce rotational motion between the upper magnet and the lower magnet-coil assembly, which adds an additional degree of freedom to the system. It was observed that the planar springs do not operate in a simple vertical (up/down) motion but rather experience out-of-plane tilting as the opposing magnetic forces push the magnets away from their vertical orientations. Consider holding two bar magnets, one in each hand, arranged with identical poles facing each other. As the magnets are moved closer to each other, the opposing force will not only push them apart along their axes but also deflect them with a force felt in the hands. Thus, the planar springs are modeled as bars with springs on their ends as opposed to a single vertically restoring spring. The center of rotation is not necessarily the center of the bar.
In deriving the theory, attention must be paid to the sign of the tilt angles, θ 1 and θ 2 . Following the diagram of Figure 2, θ 1 is negative. Thus, sin θ 1 is negative and sin θ 1 is positive, making x 1 l 11 sin θ 1 a positive displacement and x 1 + l 12 sin θ 1 a negative displacement. The diagram shows θ 1 is positive. Thus, sin θ 2 is positive and sin θ 2 is negative, making x 2 l 21 sin θ 2 a negative displacement and x 2 + l 22 sin θ 2 a positive displacement.
Referring to Figure 2, the planar spring radii are 10 m m . Thus, the bar lengths are 10 m m . The center of rotation is located at l 11 from the edge of spring 1 and l 21 from the edge of spring 2. By definition, l 12 = 10 l 11 and l 22 = 10 l 21 .
The measured linear planar spring constants are listed in Table 1. However, the spring constants, k 11 , k 12 , k 21 and k 22 , which permit both vertical and tilt motion, are unknown. Additionally, the center of rotation as given by l 11 , l 12 , l 21 and l 22 , the damping coefficients, c 1 and c 2 , and the moments of inertia, J 1 and J 2 , are also unknown.
The unknown parameters are identified by fitting the empirical output voltage spectrum to that of a lumped parameter nonlinear ordinary differential equation (ODE) model of the device. The following section develops the ODE model. Section 4 describes how the unknown parameters are fit.

3. Magnetically Coupled Double-Spring Design

Applying an energy approach [18,19,20], a Lagrangian is defined as
L x 1 , x 2 , x ˙ 1 , x ˙ 2 , θ 1 , θ 2 , θ ˙ 2 , θ ˙ 1 , q , q ˙ = K . E . P . E . Mechanical energy + W m Magnetic energy ,
where the kinetic energy is,
K . E . = 1 2 m 1 x ˙ 1 2 + 1 2 m 2 x ˙ 2 2 + 1 2 J 1 θ ˙ 1 2 + 1 2 J 2 θ ˙ 2 2 ,
the potential energy is,
P . E . = 1 2 k 11 x 1 l 11 sin θ 1 2 + 1 2 k 12 x 1 + l 12 sin θ 1 2 + 1 2 k 21 x 2 l 21 sin θ 2 2 + 1 2 k 22 x 2 + l 22 sin θ 2 2 + U 12 ( r 12 ) ,
where θ n , n = 1 , 2 is positive when the right-hand side of the bar rotates counter-clockwise.
The magnetic energy is [21]
W m = N B 12 ( r 12 ) l x 2 x 1 q ˙ + 1 2 L c q ˙ 2 .
The magnetic potential energy at a magnet located at r 2 due to a magnet at r 1 is [10]
U 12 ( r 12 ) = m 2 · B 12 ,
where the magnetic flux density at point r 2 due to a source at r 1 is
B 12 ( r 12 ) = μ 0 4 π 2 m 1 · r 12 | r 12 | 3
and the gradient, 2 , is with respect to r 2 and r 12 r 2 r 1 , and μ 0 4 × 10 7 π H / m is the permeability of free space. The physical parameters are as follows:
m n are the masses of the magnets for n = 1 , 2 ,
x n are the vertical displacements of the magnets for n = 1 , 2 ,
k n 1 , k n 2 are the linear spring constants for n = 1 , 2 ,
l n 1 , l n 2 are the spring lengths about the center of rotation for n = 1 , 2 ,
J n are the moments of inertia for n = 1 , 2 ,
m n = M n V n are the magnetic moments of the magnets for n = 1 , 2 ,
M n are the magnetizations of the magnets for n = 1 , 2 ,
V n are the volumes of the magnets for n = 1 , 2 . Note: V 1 = V 2 V ,
qis the induced charge,
L c is the coil inductance,
Nis the number of coil turns,
Ais the coil area,
lis the coil width,
B 12 ( r 12 ) is the magnetic flux density at r 2 from a magnet at r 1 ,
B 12 ( r 12 ) is B 12 , and 
U 12 ( 12 ) is the magnetic potential energy between the two magnets.
Refer to Appendix A for the physical units.
Substituting Equations (2)–(4) into Equation (1), the full Lagrangian is obtained,
L x 1 , x 2 , x ˙ 1 , x ˙ 2 , θ 1 , θ 2 , θ ˙ 1 , θ ˙ 2 , q , q ˙ = 1 2 m 1 x ˙ 1 2 + 1 2 m 2 x ˙ 2 2 + 1 2 J 1 θ ˙ 1 2 + 1 2 J 2 θ ˙ 2 2 1 2 k 11 x 1 l 11 sin θ 1 2 1 2 k 12 x 1 + l 12 sin θ 1 2 1 2 k 21 x 2 l 21 sin θ 2 2 1 2 k 22 x 2 + l 22 sin θ 2 2 U 12 ( r 12 ) + N B 12 ( r 12 ) l x 2 x 1 q ˙ + 1 2 L c q ˙ 2 .
The electromechanical dissipation is
D x 1 , x 2 , x ˙ 1 , x ˙ 2 , q , q ˙ = 1 2 c 1 x ˙ 1 2 + 1 2 c 2 x ˙ 2 2 + 1 2 R l + R c q ˙ 2 ,
where
c n are the damping constants for n = 1 , 2 ,
R c is the coil resistance, and 
R l is the load resistance.
The mechanical Lagrangian equations are
t L x ˙ n + D x ˙ n L x n = m n y ¨ ,
t L θ ˙ n + D θ ˙ n L θ n = 0 ,
where n = 1 , 2 and y ¨ is the input acceleration. The electrical Lagrangian equation is
t L q ˙ + D q ˙ L q = E 0 No external voltage .
Solving Equations (9)–(11) for the highest-order derivatives in each physical value and noting the current is the time derivative of the charge,
i = q ˙ ,
yields
m 1 x ¨ 1 = k 11 + k 12 x 1 c 1 x ˙ 1 + k 11 l 11 k 12 l 12 sin θ 1 + 3 μ 0 V 2 2 M x 2 M y 2 4 π | x 1 x 2 | 5 x 1 x 2 l N B 12 2 π i + m 1 y ¨
J 1 θ ¨ 1 = k 11 l 11 k 12 l 12 x 1 cos θ 1 k 11 l 11 2 + k 12 l 12 2 cos θ 1 sin θ 1
m 2 x ¨ 2 = k 21 + k 22 x 2 c 2 x ˙ 2 + k 21 l 21 k 22 l 22 sin θ 2 3 μ 0 V 2 2 M x 2 M y 2 4 π | x 1 x 2 | 5 x 1 x 2 + l N B 12 2 π i + m 2 y ¨
J 2 θ ¨ 2 = k 21 l 21 k 22 l 22 x 2 cos θ 2 k 21 l 21 2 + k 22 l 22 2 cos θ 2 sin θ 2
i ˙ = 1 L c 2 l N B 12 x ˙ 2 x ˙ 1 R c + R l i .
Note, the magnet volumes are assumed identical; thus, V V 1 = V 2 . The derivation details are presented in Appendix B.
Define the state vector as,
x ( t ) = x 1 ( t ) x ˙ 1 ( t ) x 2 ( t ) x ˙ 2 ( t ) θ 1 ( t ) θ ˙ 1 ( t ) θ 2 ( t ) θ ˙ 2 ( t ) i ( t ) ,
and express Equation (13) through (17) as a state-space system of ODEs,
d d t x ( t ) = A linear x ( t ) + A nonlinear ( x ( t ) ) + u ( t ) ,
where
A linear = 0 1 0 0 0 0 0 0 0 1 m 1 ( k 11 + k 12 ) c 1 m 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 m 2 ( k 21 + k 22 ) c 2 m 2 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 R l + R c L c ,
A nonlinear = 0 1 m 1 k 11 l 11 k 12 l 12 sin θ 1 + 3 μ 0 V 2 2 M x 2 M y 2 4 π x 1 x 2 5 x 1 x 2 l N B 12 2 π i 0 1 m 2 k 21 l 21 k 22 l 22 sin θ 2 3 μ 0 V 2 2 M x 2 M y 2 4 π x 1 x 2 5 x 1 x 2 + l N B 12 2 π i 0 1 J 1 k 11 l 11 k 12 l 12 x 1 cos θ 1 k 11 l 11 2 + k 12 l 12 2 cos θ 1 sin θ 1 0 1 J 2 k 21 l 21 k 22 l 22 x 2 cos θ 2 k 21 l 21 2 + k 22 l 22 2 cos θ 2 sin θ 2 2 N B 12 l L c x ˙ 2 x ˙ 1
and
u ( t ) = 0 1 0 1 0 0 0 0 0 y ¨ ( t ) .
Note, as derived in Appendix B,
B 12 = μ 0 V 4 π 4 M x 2 + M y 2 x 1 x 2 6 ,
Following Hooke’s law,
F = k x ,
one may identify the linear spring constant and the nonlinear spring function within Equations (20) and (21),
F 1 = k 11 + k 12 Linear   spring   constant x 1 + 3 μ 0 V 2 2 M x 2 M y 2 4 π | x 1 x 2 | 5 Nonlinear spring function x 1 x 2 ,
with a similar expression for x 2 . The highly nonlinear nature of the spring function results in a system which is sensitive to driving functions which are orders of magnitude less than the acceleration due to gravity, g. The driving acceleration used in the current model is of the order 10 2 g [22].
To understand the nonlinear nature of the system, consider an order of magnitude analysis of Equation (25):
x n 10 5   m ,
k n , 1 + k n , 2 10 2   N   m 1 ,
μ 0 10 6   H   m 1 ,
V 10 8   m 3 ,
2 M x 2 M y 2 10 9   A   m 1 ,
x 1 x 2 10 5   m .
So, the force order of magnitude is
F 10 3 Linear   ±   10 7 Nonlinear N .
The conclusion is the nonlinearity is ten times stronger than the linear components.

Static Condition

The system of equations described by Equation (19) through (22) represent a boundary value problem which is defined for all t > 0 . However, the initial conditions are based on an equilibrium condition between the static forces from the magnetic fields and the forces due to the static deflection of the supporting springs. Hence, when there is zero input, there must be zero output. Thus, the state-space model of Equation (19) must be identically zero for all time when ( t ) 0 :
A linear x ( t ) + A nonlinear ( x ( t ) ) = 0 .
This is the static condition. Consider holding two magnets, one in each hand, orienting them with identical poles facing each other and bringing one’s hands together. The force of the same pole arrangement will push the hands apart in distance and angle. The same occurs in the static case: The springs will displace vertically and rotate to avoid alignment of the magnets.
It is not possible to solve analytically for the static physical values of ( 0 ) , particularly the initial displacements and angles, because of the nonlinear nature of the system. Rather, the static values are determined by minimizing,
f ( x ( 0 ) ) = A linear x ( 0 ) + A nonlinear ( x ( 0 ) )
over x ( 0 ) , where the linear and rotational velocities were constrained to be zero. The static values are listed in Table 1 and Figure 3.
It is suspected there are multiple local minima in the static condition objective function. Certain combinations of the known and guessed parameters of Table 1 resulted in static angles in excess of tens of degrees. These static solutions were rejected and the parameters changed before rerunning the static optimization. The initial guesses were not altered before rerunning the optimization.
Some “reasonable” and “acceptable” static values resulted in simulations with “collisions” that are ODE solutions where the magnets crossed paths. These cases were rejected. Again, the initial guesses were not altered before rerunning the static solver and the ODE.

4. Estimation of the Effective Model Parameters

Table 1 lists the model physical parameters divided into four categories,
  • Known Those parameters whose values may be empirically measured.
  • Guesses Unknown parameters whose values were guessed based upon expertise and approximate empirical measurement.
  • Particle Swarm Optimized Fit Values Parameters whose values were determined by a particle swarm optimization which fit the model spectrum to the empirical spectrum, that minimized the difference between the two spectra. A particle swarm optimizer was selected because it was hypothesized that the nonlinear nature of the problem resulted in an error measure with multiple local minima.
  • Unconstrained Derivative-Free Optimized Static Values The device static parameters which cannot be solved analytically due to the nonlinear nature of the problem and had to be solved via an optimization method.
Identifying the unknown parameters and verifying and validating the lumped parameter nonlinear ordinary differential equation (ODE) model involved utilizing experimental measurements obtained from the vibration harvesting device prototype. The measurements served as a basis to fine-tune and validate the predictions made by the model. By comparing the experimental results with the numerical simulations, adjustments and optimizations can be made to improve the accuracy of the model’s predictions.
To accomplish the verification, a particle swarm optimization (PSO) algorithm [23,24] was employed. The PSO algorithm was utilized to optimize and fit the parameters of the theoretical model to align with the experimental results. Through an iterative process, the algorithm searches for the best parameter values that minimize the discrepancy between the model predictions and the experimental data. The cost function was the root mean square error between the empirical and modeled output voltage spectra:
f ( p ) = 1 N n = 1 N V e ( n ) α V m ( n ; p ) 2 ,
where V e ( n ) is the empirical spectrum, V m ( n ; ) is the model spectrum, N is the number of sample points and is the vector of parameters to be optimized:
p = k 11 k 12 k 21 k 22 J 1 J 2 α .
α is a positive amplitude used to account for unknown scaling between the empirical and model spectra.
The effectiveness of this calibration approach hinges upon the assumption that the theoretical model adequately captures all the relevant physics associated with the system under investigation. It is crucial for the model to approximate the behavior of the system with an acceptable degree of accuracy. By validating the model against experimental data, an assessment of whether this assumption holds true and determination of the reliability of the calibrated model can be made. The results obtained from the PSO optimization indicate that the fitted parameters of the model are not significantly different from the independently measured physical parameters, such as the inertial masses and spring stiffness. This suggests that the calibrated model successfully captures the essential characteristics of the system. The agreement between the model parameters and the independently measured values further reinforces the confidence in the accuracy and reliability of the numerical model.

5. Experimental Design and Testing

A computer controlled and calibrated experimental vibration testbed was developed for analyzing the prototype sensor designs. Figure 4 depicts the testbed and its fundamental components. The input excitation for the sensor is generated by an electrodynamic shaker, designed to operate within a bandwidth ranging from 1 Hz to 5 k Hz . The computer (not shown for clarity) interfaces with the lock-in amplifier, digital oscilloscope and a laser vibrometer. The lock-in amplifier provides a sinusoidal drive signal to the shaker which is controlled by a MATLAB-GPIB [25] interface on the computer. The generated drive signal can be swept across the experimental bandwidth with very fine resolution. The drive signal from the lock-in amplifier can be amplified 40 dB as needed to boost the output acceleration of the shaker and drive the VEH sensor into a nonlinear operational regime. A 20 c m aluminum angle bar is bolted to the shaker on one end. The VEH sensor and the reference accelerometer are mounted on the free end of the aluminum bar. The primary need for the aluminum bar is to minimize the effects of the strong magnetic field from the electrodynamic shaker on the VEH sensor prototypes. Measurements with a laser vibrometer on the free and clamped ends of the stiff aluminum beam over the frequency range (1 Hz through 200 Hz ) indicate very little variation for modest acceleration levels ( y ¨ ( t ) < 0.5 g). The experimental measurement setup incorporates three independent methods for characterizing the sensor response. The first method utilizes an optical laser vibrometer which can be focused on the inertial mass, the sensor case of the VEH sensor or on the aluminum beam. This non-contact measurement enables precise velocity measurement of the inertial mass without influencing the dynamic response of the sensor. Output from the laser vibrometer is recorded onto the scope and saved via internal storage to a thumb drive. The nature of this sensor design makes this measurement challenging as the inertial mass is designed to wobble, which reduces the signal to noise of the vibrometer and ultimately limits its utility for this situation. The second measurement is performed by an accelerometer mounted near the VEH sensor location on the rigid aluminum beam. This accelerometer captures the input acceleration, y ¨ ( t ) , experienced by the sensor base. The accelerometer provides an accurate standard reference measurement of the magnitude and transient behavior (e.g., frequency content) of the input signal y ¨ ( t ) to the VEH sensor. Also, since the accelerometer is calibrated, it is also utilized to determine the relative g-force on the sensor for different input voltages to the shaker. Finally, the third measurement is the actual output of the VEH sensor. Here, either a transient response can be recorded on the oscilloscope or a forced response function can be measured using the in-phase and quadrature capabilities of the lock-in amplifier [26]. The lock-in amplifier provides an alternative approach to measuring the response of a system by not relying upon calculating the instantaneous spectrum in a short time window, as is the case for spectrum analyzers. The nature of the phased locked loop in the lock-in amplifier allows for very accurate narrow-band measurements of the order of a few Hertz to be recorded over the bandwidth of the sweep. In addition, the time constants and filter settings associated with the dwell time can be adjusted to improve the sensitivity.
The current output from the VEH sensor is also directed into an original equipment manufacturer (OEM) energy harvesting circuit, EH300 [27]. EH300 harvesting modules are designed for storing power from weak and intermittent voltage signals generated by the VEH sensors. The modules provide solutions for long-term power storage applications for powering diagnostic sensors and communication requirements.
The characteristic input for this investigation is a small-amplitude, long-duration random noise signal. A signal of this type is what might be experienced by the sensor attached to a structure responding to non-correlated ambient vibrations. Here, the emphasis is to test the VEH sensor’s bandwidth and performance to low-level input stimuli in a controlled and repeatable manner. Harvesting performance of the current VEH sensor prototype is investigated in the next section.

6. Empirical and Model Analysis of the VEH Sensor

In this section, the performance of the VEH sensor is investigated using spectral techniques over a limited bandwidth. The goal is to determine an estimate of the harvester’s sensitivity ( V / g ) and power conversion efficiency to low-level vibration input. The first approach analyzes the response spectrum of the sensor utilizing the narrowband capabilities of a lock-in amplifier and a swept sinusoidal drive signal. The second approach utilizes a random noise signal to measure the harvesting performance when the sensor is exposed to stochastic non-stationary vibration input. Here, the noise signal is defined to cover the same bandwidth as determined for the harvesting sensitivity. A calibrated accelerometer is used as a reference to maintain similar input characteristics for the two methods.

6.1. Swept Sine Measurements of VEH Sensor and Optimization of Model

A range of sinusoidal input voltages are generated using the internal frequency generator of the lock-in amplifier (Stanford Research, SRS856, Stanford Research Systems, Sunnyvale, CA, USA). The signal is swept across the estimated operational bandwidth of the prototype. The drive voltage to the shaker is adjusted using a reference accelerometer (PCB Piezotronics, model 352A24, Depew, NY, USA). The accelerometer is mounted on the end of the beam and provides a measurement of the input acceleration experienced by the harvesting sensor. Using the calibrated accelerometer, the drive voltage from the power amplifier was adjusted to create an approximate input acceleration from the shaker of y ¨ ( t ) = 10 2 g across the range of 1 Hz to 200 Hz . The controlled response of the shaker (BK vibration exciter 4809, Marlborough, MA, USA) across this frequency range and drive level is flat. A laser vibrometer is used to verify the shaker response. Additional vibrometer measurements detected no variation in acceleration magnitude across the length of the aluminum bar either, thus verifying the shaker and the attached aluminum bar are reproducing the drive signal characteristics from either the lock-in amplifier or the signal generator.
The goal is to determine an estimate of the harvester’s coupling sensitivity ( V / g ) and power conversion efficiency for low-level acceleration input. Figure 5 shows the empirical voltage spectrum from the VEH prototype measured by the SRS 856. The load resistance is 50 Ω and the input drive acceleration from the shaker is 8 × 10−3 g. The spectrum shows one large peak at 138 Hz and a secondary smaller one at near 140 Hz . The larger peak at 138 Hz has a usable bandwidth of approximately 20 Hz . Outside this range, the sensor coupling sensitivity of the prototype is insufficient for harvesting; this will be elaborated upon in the next section. This measurement can now be utilized to improve the analytical model. Since the input acceleration is constant, the root mean square (rms) sensitivity of the sensor is estimated as (rms(Vout)/10−2 g); so, between 130 Hz and 150 Hz , the rms coupling sensitivity is 0.58 m V / g , with a peak sensitivity of 0.87 m V / g at 138 Hz .
The optimization of the model requires tuning the parameters of the model to fit a given measured output. The parameters utilized in the model described in Equations (19)–(22) are listed in Table 1. Some of these parameters can be easily measured while others cannot. The initial fabrication data of the components (i.e., mass, magnetic strength, dimensions, estimated spring stiffness, etc.) are measured using standard laboratory techniques. The spring stiffness is estimated using a measurement of the first resonances of each of the magnet spring sub-assemblies using the laser vibrometer. The magnetic field strength of each magnet was measured using a EXTECH AC/DC Magnetic Meter, Nashua, NH, USA. These techniques will introduce errors when coupled with the complexity of the equations of motion. It is found through trial and error that the stiffness estimates for the spherical springs and the rotational inertial terms have the strongest effect on the response of the model for the low-level input acceleration. The optimization step is utilized in conjunction with the parameters that have the least uncertainty to estimate the parameters that have higher uncertainty. With this approach, the model is “tuned” to fit the empirical results. For each optimization iteration, the model of Equation (19) is run for 30 s of model time using MATLAB’s ODE45 [25] solver (this is also used to solve the static case of Equation (34)). The model output is the state vector of Equation (18). The time series of interest is the current, i ( t ) , which when multiplied by the 50 Ω load resistance yields the harvested voltage, v ( t ) . The model voltage spectrum, V m ( n ; ) , is estimated using the Welch [28] method with a window of 8.2 s with a 50% overlap. Here, the empirical spectrum, V e ( n ) , is taken from Figure 5.
Figure 6 shows the optimized model output using the parameters of Table 1 along with the harvested voltage spectrum from the VEH prototype over the frequency range 20−200 Hz . The main spectral peak matches the empirical spectrum and follows the trend on the left-hand side of the peak. The 4 degree-of-freedom lumped parameter model (two linear displacements and two rotations from each magnet) follows the empirical spectrum well after parameter fitting. The root mean square error is 1.4 × 10 7 .
As a measure of similarity, the normalized root mean squared error was chosen. It is defined as
s = 100 1 V e V m V e V ¯ e ,
where V ¯ e is the mean of the empirical spectrum. The similarity is 81.7% (an identical match is 100%).
The similarity measure provides confidence the model parameters have been tuned such that the model is accurately predicting the response of the VEH prototype. The empirical data and the fitted model show that the sensor has the highest potential harvesting performance over a 20 Hz bandwidth centered at approximately 140 Hz . In the next section, the response of the sensor to low-level random noise is investigated for harvesting applications.

6.2. Harvesting from Low-Level Non-Stationary Signals

Long-term storage and monitoring application require successful harvesting of potentially low-level and infrequent sources of vibrational energy. In this second analysis, the energy harvesting capability of the VEH sensor is investigated utilizing a random acceleration spectrum generated from the white noise feature of the Keysight 3352B arbitrary waveform generator, Santa Rosa, CA, USA (AWG). For these tests, the AWG drives the shaker directly and the accelerometer is used to set the input acceleration magnitude of y ¨ RMS 10 2 g, as determined from the reference accelerometer. The targeted noise spectrum is assumed white, centered at 130 Hz , with a 75 Hz bandwidth.
Transient accelerometer data and the voltage output from the VEH harvester are recorded over a 7 min time window using a digital oscilloscope. The sampling rate for the oscilloscope is set at 0.2 m s , which resulted in a data record of approximately 20 k samples . The frequency spectra are calculated from the transient data using Fourier methods [28]. Figure 7 shows the acceleration input at the base of the harvesting sensor and the output spectra of the VEH harvester measured with a 50 Ω load termination. The black line is the input acceleration spectrum measured by the reference accelerometer and converted to g   m   s 2 . The red line is the VEH sensor voltage output spectrum corresponding to an input random acceleration of y ¨ RMS 10 2 g for comparison.
Initially, it was thought that the random acceleration input spectrum would be flat across the bandwidth and the VEH sensor would be driven by a stochastically uniform acceleration input. The calculated input acceleration spectrum is found to be biased towards a peak rms amplitude around 75 Hz , with an approximate 2 dB/decade on either side. Unfortunately, a large portion of the energy in this spectrum falls outside the targeted range of the VEH; however, it still provides a test of the sensor’s off nominal harvesting capability.
The highlighted region of the graph is the useful harvesting bandwidth ( 130 Hz < VEH frequency < 150 Hz ) of the sensor with the EH300 harvesting circuit. Essentially, this region represents an estimate based on the EH300 circuit’s power harvesting capabilities defined as 4 V at 200 n A or 0.8 μ W [27]. Across the 50 Ω load resistance, this translates into a minimum power of 0.04 μ W across the 20 Hz operational bandwidth of the highlighted region and is marked by the red dashed line relative to the voltage axis. A smaller region exists for the second peak; however, this is not considered for this analysis.
An estimate of the harvesting efficiency of the VEH sensor based on the random acceleration is shown in Figure 8. The harvesting efficiency is defined as η = P / P in , where P is the average output power from the harvesting circuit, calculated as P = V out 2 / R l Δ f , and P in is the mechanical input power, calculated as P in = Y ¨ 2 m Δ f . Here, m are the magnet masses and Δf is the frequency band [29].
The ratio between the output voltage and the input acceleration spectrum are plotted as a function of frequency across the harvesting bandwidth, essentially normalizing the VEH sensor output sensitivity by dividing out the noise signal response. The similarity between the Figure 5 and Figure 8 is evident in the bandwidth (130−150 Hz ). In this region, the sensor has an average rms harvesting efficiency of 6.3%, with a peak value of 9.8% at 138 Hz . Here, the impedance match between the coil resistance, 60 Ω , and the load resistance, 50 Ω , supports that the approximation is close to optimal in terms of impedance matching with the load resistance. In a third set of measurements, the harvesting performance of the VEH prototype is recorded while connected to the EH300 harvesting circuit. The harvested voltage across the 50 Ω load resistor is connected to the input buffering circuit of the EH300. The oscilloscope is now used to record the voltage gain across the circuit’s storage capacitor, and as before, the sensor is exposed to the same band limited random acceleration signal described previously. The total harvested voltage from the VEH sensor in the EH300 is recorded for approximately 7 min. The sampling rate for the oscilloscope was set at 0.2 m s , which resulted in a data record of approximately 20 k samples . This is the largest window that could be held in the oscilloscope’s window before the first samples are lost. The sampling rate ensures a Nyquist limit of 5 k Hz , which is a factor of 25 times larger than the highest frequency associated with the noise spectrum. Figure 9 shows a linear fit of the data results in a voltage harvesting rate of approximately 10 m V / h for random acceleration noise y ¨ rms 10 2 g. The spikes in the temporal voltage signal may be generated in part by electronic noise present in the unshielded circuit. Another consideration is that the random acceleration magnitudes may be introducing the fluctuations in the voltage in the storage capacitor. Finally, it is possible to estimate the VEH sensor’s power density by normalizing the output power with the volume, Vol mag , of the effective moving elements (i.e., the volume of the two coil and magnet assemblies [29]). This can be calculated as P = V out 2 / R l Δ f / Vol mag , with a resulting power density of 0.2 μ W / c m 3 . With this information, the prototype sensor can be compared to existing results published in the literature [30].

7. Discussion

A prototype of a nonlinear broadband vibrational energy harvester (VEH) was modeled, designed, and built. The experimental performance was investigated over an input range of 1 Hz to 200 Hz with low-level random vibrations y ¨ ( t ) g . Over this range, the useful harvesting bandwidth is determined to be due to the sensitivity of the harvesting sensor and the minimum detectable voltage level from the EH300 harvesting circuit. The tested prototype sensor demonstrated the capability of harvesting usable energy from random ambient vibrations of the order of 10 3 g over the operating range of 130 Hz < f < 150 Hz . This underscores the device’s potential to harness energy from minute mechanical disturbances. While the bandwidth is limited to 20 Hz , further testing may improve the performance. The small size of the inertial masses plays a large role in the device’s low efficiency; however, this might be improved with further optimization. In the context of longer duration exposure, the sensor does generate energy which can slowly be stored. The EH300 harvesting circuit is responsible for storing the accumulated energy over time. It is observed that the VEH had a charging performance of 10 m V   h 1 , based on the measured voltage across the storage capacitor (1000 n F ) in the EH300 circuit. Based on this performance estimate, the current prototype would take approximately 40 h to fully cycle the EH300 to an operational state suitable to support typical IoT microprocessor applications. This long charging window is due to the misalignment of the input acceleration spectrum and the VEH sensor’s peak sensitivity. An order of magnitude improvement would be realized if the sensor’s design point matched the input acceleration spectrum. The EH300 provided a convenient method for energy capture, with similar approaches given in the literature [31]; however, an investigation of other potential harvesting circuits was not conducted. The EH300 harvesting circuit can also charge different external capacitors with different capacitances, providing the designer additional flexibility and capabilities [32].
Equations of motion for the coupled electro-mechanical system (ODE) model were developed and used to predict the harvesting performance. With the fitted parameters, the model output voltage spectrum shows good correlation with the model over the predicted resonances, with the overall trend between 20 Hz and 200 Hz showing a similarity measure of 81.7%. The similarity in the region of interest, 130 Hz to 150 Hz , is 73.8%. The decrease in similarity over the 130−150 Hz range is an artifact of fitting the model to the data over the entire spectrum, 20−200 Hz . Because of the fitting process, the model is biased towards the entire measured spectra, so it is not surprising that the fitted model deviates from the experiment in this reduced range. The similarity would improve by refitting the entire model over the smaller bandwidth; however, this was not done in this investigation. The comparison between the model’s performance and the experimental results provides interesting insights into the behavior of the sensor design. Notably, the presence of asymmetries linked to magnetic fields and height disparities within the mounting hardware introduces a rotational motion component into the sensor response, a factor that the experiment reveals but the model only partially captures. These findings highlight the complexities of the system and emphasize the necessity of refining the model to encompass the nuanced 3D effects and accurately represent the future experimental comparisons. The design also possesses the capability to be mechanically adjusted (tuned) to accommodate different tilt angles and magnetic separation distances. The adjusting process is simple but a bit tedious. With the sensor on the test stand, the hex nuts can be loosened and moved up or down on each threaded rod using a small non-magnetic wrench. The adjustable angular travel is ∼ ± 3 degrees and the vertical travel is ∼ ± 4 m m . The repelling forces of the two magnets also contribute to the angular tilt and the overall separation between the coil and the upper magnet. Adjustments can be used in conjunction with the model parameters to improve harvesting and modify performance; however, this was not done in this investigation.

8. Conclusions

Overall, the current sensor was shown to be sensitive at capturing low-level random noise centered near its design spectrum. Further tuning by adjusting the tilt and separation of the magnets may improve the performance. Vibration harvesting sensors targeted to harvest low-level ambient vibration energy face many design challenges. Targeting the correct operational spectrum such that the sensor dynamics take full advantage of the ambient vibrational field can be aided by models that are tuned to the target spectrum. Once tuned, the model then indicates which parameters are more effective at controlling the performance. Future investigations will focus on developing sensor designs and operational parameters that optimize the harvesting process.

Author Contributions

Hardware design and construction, K.A.F.; theory development and analysis, K.A.F. and S.K.L.; modeling and simulation, S.K.L. All authors have read and agreed to the published version of the manuscript.

Funding

The Department of Energy Office of Weapon Technology Manufacturing and Maturation under Contract DE-AC52-07NA27344.

Institutional Review Board Statement

This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. Document release number: LLNL-JRNL-858500.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are not publicly available until they have been reviewed and released by Lawrence Livermore National Laboratory.

Acknowledgments

The authors thank David Obenhauf and Bert Vollmer for their support.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Units

Table A1. Summary of parameter units.
Table A1. Summary of parameter units.
VariablePropertyUnit
x ( t ) Vertical displacement m
θ ( t ) Rotation angle rad
y ¨ ( t ) Input acceleration m   s 2
mMass k g
k s Linear spring constant k g   s 2 = N   m 1
cDamping coefficient k g   s 1 = N   s   m 1
JMoment of inertia k g m 2 / rad
BMagnetic flux density T = V   s   m 2 = N / ( A m )
U 12 Potential energy between two magnets N m
NNumber of coil turnsUnitless
lCoil length m
μ Permeability H   m 1 = N   A 2
M Magnetization A   m 1
m M V Magnetic moment (magnetization × magnet volume) A   m 2
VMagnet volume m 3
RResistance Ω = m 2 k g / ( A 2   s 3 )
LInductance H = k g   m 2 / ( A s ) 2
CCapacitance F = C   V 1 = s 4 A 2 / ( k g   m 2 )
iCurrent A

Appendix B. Derivation Notes

The magnet coordinate system is
r n = y n x n z ,
for n = 1 , 2 . x n is the vertical displacement, as shown in Figure 2, and y n is the horizontal displacement. Note: y n is not to be confused with y, which is the input acceleration. There is no displacement into the page, thus z 0 .
The magnetic flux density of a magnet located at r 1 with magnetic moment m 1 at an observation point r 2 is [10],
B 12 ( r 12 ) = μ 0 4 π 2 m 1 · r 12 | r 12 | 3 , = μ 0 M y V 4 π | x 2 x 1 | 3 μ 0 M x V 2 π | x 2 x 1 | 3 0 ,
where the gradient 2 is with respect to r 2 and r 12 r 2 r 1 . Note,
| B 12 | = B 12 = μ 0 V 4 π 4 M x 2 + M y 2 x 1 x 2 6 .
The potential energy between magnets located at r n , n = 1 , 2 , is
U 12 ( r 12 ) = m 2 · B 12 ( r 12 ) , = μ 0 V 2 2 M x 2 M y 2 4 π | x 1 x 2 | 3 .
The Lagrangian derivation of the ODE model requires the derivative of the potential energy with respect to x n , n = 1 , 2 :
x 1 U 12 = 3 μ 0 V 2 2 M x 2 M y 2 4 π | x 1 x 2 | 4 d d x | x | x x 1 x 2 , = 3 μ 0 V 2 2 M x 2 M y 2 4 π | x 1 x 2 | 4 x 1 x 2 | x 1 x 2 | , = 3 μ 0 V 2 2 M x 2 M y 2 4 π | x 1 x 2 | 5 x 1 x 2 ,
x 2 U 12 = 3 μ 0 V 2 2 M x 2 M y 2 4 π x 1 x 2 5 x 1 x 2 .

References

  1. Harne, R.L.; Wang, K.W. A Review of the Recent Research on Vibration Energy Harvesting via Bistable Systems. Smart Mater. Struct. 2013, 22, 023001. [Google Scholar] [CrossRef]
  2. Lee, J.; Yuen, S.; Li, W.; Leong, P. Development of an AA size energy transducer with micro resonators. In Proceedings of the 2003 IEEE International Symposium on Circuits and Systems (ISCAS), Bangkok, Thailand, 25–28 May 2003; Volume 4. [Google Scholar]
  3. Seeman, M.; Sanders, S.; Rabaey, J. An ultra-low-power power management IC for wireless sensor nodes. In Proceedings of the 2007 IEEE Custom Integrated Circuits Conference, San Jose, CA, USA, 16–19 September 2007; pp. 925–931. [Google Scholar]
  4. Jia, Y. Review of nonlinear vibration energy harvesting: Duffing, bistability, parametric, stochastic and others. J. Intell. Mater. Syst. Struct. 2020, 31, 921–944. [Google Scholar] [CrossRef]
  5. Waterbury, A. Vibration Harvesting Using Electromagnetic Transduction. Ph.D. Thesis, University of California, Berkeley, CA, USA, 2011. [Google Scholar]
  6. Chai, H.D. Electromechanical Motion Devices; Prentice-Hall PTR: Hoboken, NJ, USA, 1988. [Google Scholar]
  7. Junxiang, J.; Liu, S.; Feng, L.; Zhao, D. A Review of Piezoelectric Vibration Energy Harvesting with Magnetic Coupling Based on Different Structural Characteristics. Micromachines 2021, 12, 436. [Google Scholar] [CrossRef]
  8. Stanton, S.C.; McGehee, C.C.; Mann, B.P. Reversible hysteresis for broadband magnetopiezoelastic energy harvesting. Appl. Phys. Lett. 2010, 95, 174103. [Google Scholar] [CrossRef]
  9. Zhang, H.; Corr, L.R.; Ma, T. Issues in vibration energy harvesting. J. Sound Vib. 2018, 421, 79–90. [Google Scholar] [CrossRef]
  10. Furlani, E.P. Permanent Magnet and Electromechanical Devices: Materials, Analysis and Applications; Academic Press: Cambridge, MA, USA, 2001. [Google Scholar]
  11. Carneiro, P.; dos Santos, M.P.; Rodrigues, A.; Ferreira, J.A.; Simões, J.A.; Marques, A.T.; Kholkin, A.L. Electromagnetic energy harvesting using magnetic levitation architectures: A review. Appl. Energy 2020, 260, 114191. [Google Scholar] [CrossRef]
  12. Muscat, A.; Bhattacharya, S.; Zhu, Y. Electromagnetic Vibrational Energy Harvesters: A Review. Sensors 2022, 22, 5555. [Google Scholar] [CrossRef] [PubMed]
  13. Beeby, S.P.; Torah, R.N.; Tudor, M.J.; Glynne-Jones, P.; O’Donnell, T.; Saha, C.R.; Roy, S. A micro electromagnetic generator for vibration energy harvesting. J. Micromech. Microeng. 2007, 17, 1257. [Google Scholar] [CrossRef]
  14. Mann, B.; Sims, N. Energy harvesting from the nonlinear oscillations of magnetic levitation. J. Sound Vib. 2009, 319, 515–530. [Google Scholar] [CrossRef]
  15. Nico, V.; Boco, E.; Frizzell, R.; Punch, J. A high figure of merit vibrational energy harvester for low frequency applications. Appl. Phys. Lett. 2016, 108, 013902. [Google Scholar] [CrossRef]
  16. Zhou, S.; Cao, J.; Inman, D.J.; Lin, J.; Liu, S.; Wang, Z. Broadband tristable energy harvester: Modeling and experiment verification. Appl. Energy 2014, 133, 33–39. [Google Scholar] [CrossRef]
  17. Zhou, S.; Zuo, L. Nonlinear dynamic analysis of asymmetric tri-stable energy harvesters for enhanced energy harvesting. Commun. Nonlinear Sci. Numer. Simul. 2018, 61, 271–284. [Google Scholar] [CrossRef]
  18. Goldstein, H. Classical Mechanics; Addison-Wesley: Boston, MA, USA, 1980. [Google Scholar]
  19. Landau, L.D.; Lifshits, E.M. Mechanics, 3rd ed.; Butterworth-Heinemann: Oxford, UK, 1976. [Google Scholar]
  20. Preumont, A. Mechatronics: Dynamics of Electromechanical and Piezoelectric Systems; Solid Mechanics and Its Applications; Springer: Dordrecht, The Netherlands, 2006. [Google Scholar]
  21. Walker, J. Halliday & Resnick Fundamentals of Physics, 10th ed.; Wiley: Hoboken, NJ, USA, 2014. [Google Scholar]
  22. Neri, I.; Travasso, F.; Mincigrucci, R.; Orfei, H.V.F.; Gammaitoni, L. A real vibration database for kineticenergy harvesting application. J. Intell. Mater. Syst. Struct 2012, 23, 2095–2101. [Google Scholar] [CrossRef]
  23. Harrison, K.R.; Engelbrecht, A.P.; Ombuki-Berman, B. Self-adaptive particle swarm optimization: A review and analysis of convergence. Swarm Intell. 2018, 12, 187–226. [Google Scholar] [CrossRef]
  24. Pedersen, M.E.H. Good Parameters for Particle Swarm Optimization; Technical Report; Hvass Laboratories: Luxembourg, 2010. [Google Scholar]
  25. MATLAB; The MathWorks, Inc.: Natick, MA, USA, 2024.
  26. Elhady, A.; Alghamdi, M.S.; Abdel-Rahman, E. Experimental construction of force- and frequency-response curves of nonlinear resonators. Chaos 2023, 33, 063102. [Google Scholar] [CrossRef]
  27. EH300; Advanced Linear Devices, Inc.: Sunnyvale, CA, USA, 2024.
  28. Oppenheim, A.V.; Schafer, R.W.; Buck, J.R. Discrete-Time Signal Processing, 2nd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 1999; ISBN 0137549202. [Google Scholar]
  29. Roundy, S. On the Effectiveness of Vibration-based Energy Harvesting. J. Intell. Mater. Syst. Struct. 2005, 16, 809–823. [Google Scholar] [CrossRef]
  30. Brenes, A.; Morel, A.; Gibus, D.; Yoo, C.S.; Gasnier, P.; Lefeuvre, A.B.E. Large-bandwidth piezoelectric energy harvesting with frequency-tuning synchronized electric charge extraction. Sens. Actuators A Phys. 2020, 302, 111759. [Google Scholar] [CrossRef]
  31. Rodriguez, J.C.; Nico, V.; Punch, J. A Vibration Energy Harvester and Power Management Solution for Battery-Free Operation of Wireless Sensor Nodes. Sensors 2019, 19, 3776. [Google Scholar] [CrossRef] [PubMed]
  32. Advanced Linear Devices EH300. Available online: https://www.mouser.com/ProductDetail/Advanced-Linear-Devices/EH300?qs=NoieiqYAh1AkX2jUqjbnpw%3D%3D (accessed on 5 August 2024).
Figure 1. (Left): The dual magnet device. Threaded rods are utilized to adjust the separation distance of the magnets and their angular orientation (Right): The planar spring facilitates out-of-plane displacement and rotation of the supported upper magnet and the lower magnet/coil assembly.
Figure 1. (Left): The dual magnet device. Threaded rods are utilized to adjust the separation distance of the magnets and their angular orientation (Right): The planar spring facilitates out-of-plane displacement and rotation of the supported upper magnet and the lower magnet/coil assembly.
Applsci 14 07072 g001
Figure 2. Lumped parameter one-dimensional model with ODE parameters of the device shown in Figure 1. The top magnet is represented by the mass m 1 . The combination bottom magnet and coil is represented by mass m 2 .
Figure 2. Lumped parameter one-dimensional model with ODE parameters of the device shown in Figure 1. The top magnet is represented by the mass m 1 . The combination bottom magnet and coil is represented by mass m 2 .
Applsci 14 07072 g002
Figure 3. The counterposed magnetic poles cause the magnets to deflect and separate. The static values are determined by numerical optimization.
Figure 3. The counterposed magnetic poles cause the magnets to deflect and separate. The static values are determined by numerical optimization.
Applsci 14 07072 g003
Figure 4. Diagram of the experimental testing setup. The system can simultaneously measure the coil current, i ( t ) , inertial mass velocity, x ˙ ( t ) , input acceleration, y ¨ ( t ) and the harvesting performance (volts per time). The input excitation can be amplified to push operation of the VEH sensor into a nonlinear regime as needed. The computer-controlled lock-in amplifier is utilized as both a signal source and as a tool to record the response of the sensor.
Figure 4. Diagram of the experimental testing setup. The system can simultaneously measure the coil current, i ( t ) , inertial mass velocity, x ˙ ( t ) , input acceleration, y ¨ ( t ) and the harvesting performance (volts per time). The input excitation can be amplified to push operation of the VEH sensor into a nonlinear regime as needed. The computer-controlled lock-in amplifier is utilized as both a signal source and as a tool to record the response of the sensor.
Applsci 14 07072 g004
Figure 5. VEH coupling sensitivity as measured across a 50 Ω load resistance for a constant input acceleration of 8 × 10 3 g over 1−200 Hz . The spectrum shows one large peak at 138 Hz and a secondary smaller one at near 140 Hz . The larger peak at 138 Hz has a usable bandwidth of approximately 20 Hz between 130 and 150 Hz (highlighted region). The effective rms coupling sensitivity in this region is 0.58 m V / g with a peak value of 0.87 m V / g at 138 Hz . Outside this range, the sensor coupling sensitivity of the prototype is not considered for harvesting.
Figure 5. VEH coupling sensitivity as measured across a 50 Ω load resistance for a constant input acceleration of 8 × 10 3 g over 1−200 Hz . The spectrum shows one large peak at 138 Hz and a secondary smaller one at near 140 Hz . The larger peak at 138 Hz has a usable bandwidth of approximately 20 Hz between 130 and 150 Hz (highlighted region). The effective rms coupling sensitivity in this region is 0.58 m V / g with a peak value of 0.87 m V / g at 138 Hz . Outside this range, the sensor coupling sensitivity of the prototype is not considered for harvesting.
Applsci 14 07072 g005
Figure 6. Results of a 30 s simulated run with a 20−200 Hz input chirp of amplitude 8 × 10 3 g. (a) The normalized root mean square similarity measure is 81.7%. (b) Detail of the 130−150 Hz useful region of interest where the similarity measure is 78.3%. The decrease in similarity over the region of interest is an artifact of fitting the model and the data over the entire spectrum (20−200 Hz ). The similarity would improve by refitting the entire model over the smaller bandwidth; however, this was not done in this investigation.
Figure 6. Results of a 30 s simulated run with a 20−200 Hz input chirp of amplitude 8 × 10 3 g. (a) The normalized root mean square similarity measure is 81.7%. (b) Detail of the 130−150 Hz useful region of interest where the similarity measure is 78.3%. The decrease in similarity over the region of interest is an artifact of fitting the model and the data over the entire spectrum (20−200 Hz ). The similarity would improve by refitting the entire model over the smaller bandwidth; however, this was not done in this investigation.
Applsci 14 07072 g006
Figure 7. VEH sensor response to stochastic, non-stationary vibration input. Black line: Input acceleration spectrum in g ( m   s 2 ) of the noise source. Red line: VEH sensor voltage power output spectrum. The output voltage is measured with a 50 Ω load termination. The highlighted region of the graph is the useful sensor harvesting bandwidth, 130−150 Hz . The red dashed horizontal line is the lowest harvestable power attainable by the EH300 circuit, which is estimated to be 0.04 μ W (see text).
Figure 7. VEH sensor response to stochastic, non-stationary vibration input. Black line: Input acceleration spectrum in g ( m   s 2 ) of the noise source. Red line: VEH sensor voltage power output spectrum. The output voltage is measured with a 50 Ω load termination. The highlighted region of the graph is the useful sensor harvesting bandwidth, 130−150 Hz . The red dashed horizontal line is the lowest harvestable power attainable by the EH300 circuit, which is estimated to be 0.04 μ W (see text).
Applsci 14 07072 g007
Figure 8. VEH sensor harvesting efficiency defined in terms of the ratio of the powers: η = P out / P in . Here, P out is defined as the output voltage power spectrum per Hz across the 50 Ω load resister (the red line of Figure 7). P in represents the mechanical power spectrum per Hz of the combined mass of the coil and magnets (the black line of Figure 7) as a function of frequency over the 130−150 Hz region of interest. The rms harvesting efficiency for this region is η rms = 6.3 %, with a peak harvesting efficiency of η pk = 9.7 % at 138 Hz .
Figure 8. VEH sensor harvesting efficiency defined in terms of the ratio of the powers: η = P out / P in . Here, P out is defined as the output voltage power spectrum per Hz across the 50 Ω load resister (the red line of Figure 7). P in represents the mechanical power spectrum per Hz of the combined mass of the coil and magnets (the black line of Figure 7) as a function of frequency over the 130−150 Hz region of interest. The rms harvesting efficiency for this region is η rms = 6.3 %, with a peak harvesting efficiency of η pk = 9.7 % at 138 Hz .
Applsci 14 07072 g008
Figure 9. VEH sensor charging performance over a seven−minute charging window using the EH300 interface circuit. Blue line: Raw DC voltage across the output terminals from the EH300. Red line: Fitted data using least squares approximation. The measured harvesting performance is approximately 10 m V / h , as harvested from the random acceleration input acceleration power spectra shown in Figure 7.
Figure 9. VEH sensor charging performance over a seven−minute charging window using the EH300 interface circuit. Blue line: Raw DC voltage across the output terminals from the EH300. Red line: Fitted data using least squares approximation. The measured harvesting performance is approximately 10 m V / h , as harvested from the random acceleration input acceleration power spectra shown in Figure 7.
Applsci 14 07072 g009
Table 1. Physicalparameters used in the simulation which approximately match those of the experiment.
Table 1. Physicalparameters used in the simulation which approximately match those of the experiment.
ParameterValueUnit
Known
Magnet masses, m n ,   n = 1 , 2 1.52, 0.453 g
Linear spring constants, k s , n ,   n = 1 , 2 649, 160 k g   s 2
Planar spring radius, l10 m m
Magnetic flux density peak amplitude, B1 T
Damping, c n ,   n = 1 , 2 5 × 10−2, 5 × 10−2 k g   s 1
Number of coil turns, N650Unitless
Coil area, A25 π m m 2
Coil resistance, R c 60 Ω
Coil inductance, L c 1.3 H
Load resistance, R l 50 Ω
Guesses
Magnet radii, ρ 1 m m
Magnet height, h5 m m
Magnet volume, V n π ρ 2 h = 15.7 m m 3
Magnetizations, M n x , M n y , M n z , n = 1 , 2 5 × 10−4, 10−4, 10−4 A   m 1
Spring radii, l 11 , l 12 , l 21 , l 22 2, 8, 3, 7 m m
Damping, c 1 , c 2 0.05, 0.05 N   s   m 1
Particle Swarm Optimized Fit Values
Linear spring constants, k 11 , k 12 559, 461 k g   s 2
Linear spring constants, k 21 , k 22 130, 161 k g   s 2
Moments of inertia, J 1 , J 2 2.2 × 10−2, 4.8 × 10−2 k g   m 2
Unconstrained Derivative-Free Optimized Static Values
x 1 ( 0 ) , x 2 ( 0 ) 0.05, −24.14 μ m
θ 1 ( 0 ) , θ 2 ( 0 ) 0.23, 0.08degrees
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lehman, S.K.; Fisher, K.A. A Vibrational Energy Harvesting Sensor Based on Linear and Rotational Electromechanical Effects. Appl. Sci. 2024, 14, 7072. https://doi.org/10.3390/app14167072

AMA Style

Lehman SK, Fisher KA. A Vibrational Energy Harvesting Sensor Based on Linear and Rotational Electromechanical Effects. Applied Sciences. 2024; 14(16):7072. https://doi.org/10.3390/app14167072

Chicago/Turabian Style

Lehman, Sean K., and Karl A. Fisher. 2024. "A Vibrational Energy Harvesting Sensor Based on Linear and Rotational Electromechanical Effects" Applied Sciences 14, no. 16: 7072. https://doi.org/10.3390/app14167072

APA Style

Lehman, S. K., & Fisher, K. A. (2024). A Vibrational Energy Harvesting Sensor Based on Linear and Rotational Electromechanical Effects. Applied Sciences, 14(16), 7072. https://doi.org/10.3390/app14167072

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