A Metamaterial-Inspired Structure for Simultaneous Vibration Attenuation and Energy Harvesting

In this article, a magnetomechanical metamaterial structure capable of simultaneous vibration attenuation and energy harvesting is presented. )e structure consists of periodically arranged local resonators combining cantilever beams and permanent magnetcoil systems. A prototype of the metamaterial dual-function structure is fabricated, and models are developed. Results show good agreement between model simulation and experiment. Two frequency bandgaps are measured: 205–257Hz and 587–639Hz. Within these bandgaps, vibrations are completely attenuated.)e level of vibration attenuation in the first bandgap is substantially larger than the level of vibration attenuation observed in the second bandgap. Mode shapes suggest that bending deformations experienced by the local resonators in the second bandgap are less than the deformations experienced in the first bandgap, and most vibrational energy is localized within the first bandgap where the fundamental resonant frequency is located, i.e., 224Hz.)e ability of the fabricated metamaterial structure to harvest electric power in these bandgaps is examined. Results show that vibration attenuation and energy harvesting characteristics of the metamaterial structure are coupled. Stronger vibration attenuation within the first bandgap has led to enhanced energy harvesting capabilities within this bandgap. Power measurements at optimum load resistance of 15Ω reveal that maximum power generated within the first bandgap reaches 5.2 μW at 245Hz. Compared with state-of-the-art, the metamaterial structure presented here shows a significant improvement in electric power generation, at considerably lower load resistance, while maintaining the ability to attenuate undesired vibrations within the frequency bandgap.


Introduction
Elastic mechanical metamaterials have shown extraordinary properties including bandgap phenomenon [1][2][3]. Consequently, elastic metamaterials are gaining unprecedented attention in advanced engineering applications such as vibration attenuation and wave filtering [4][5][6] and energy harvesting [7,8]. Typically, the overall mechanical metamaterial structure is man-made and consists of periodically arranged cavities housing locally resonating masses that are connected to the cavities by elastic elements, i.e., springs [9]. e bandgap phenomenon is the result of the dynamic behavior of the periodically arranged local resonators and Bragg's scattering in periodic structures. e overall dynamics of the metamaterial structure is determined by the locally resonating masses. When the mechanical metamaterial structure is subject to excitations from external source, the kinetic energy contained in these oscillations is trapped in the local resonators. is results in the unique capability of generating frequency bandgaps within which oscillations are prevented from passing through. e generated bandgaps are formed at frequency matching the fundamental resonant frequency of the local resonators that are acting as local absorbers within the metamaterial structure [2,10].
To this end, research interests in mechanical metamaterials have grown in a few directions including: vibration attenuation and wave filtering, vibration energy harvesting, and simultaneous vibration attenuation and energy harvesting. For instance, Matlack et al. presented an elastic metastructure for broadband vibration absorption [5]. e structure was additively manufactured using fused deposition-modeling technique and then assembled manually in order to insert metallic cubes inside the printed structure. Moreover, a metamaterial structure embedded with hierarchically organized local resonators was fabricated and characterized by Xu et al. [11]. Both experimental and numerical methods were used to fully characterize the fabricated structure. Results from this study demonstrated the ability of the fabricated structure to attenuate vibrations within two bandgaps in frequency regions of 144-188 Hz and 480-810 Hz. In a similar fashion, Zhu et al. built a chiral lattice-based metamaterial beam that is embedded with local resonators to achieve broadband vibration attenuation in the frequency range of, approximately, 210-700 Hz [12]. e use of metamaterial structures to generate electric power via energy harvesting has also been investigated recently. A thorough review of metamaterials-based energy harvesters was presented by Chen et al. [7]. For example, a metamaterial energy harvester employing piezoelectrics was fabricated [13]. e harvester was shown to produce peak power in the order to 100 µW-1 mW at 55 kHz and 30-80 kΩ. Concurrent vibration attenuation and energy harvesting metamaterials have been recently, also, demonstrated. Chen et al. built a metamaterial beam embedded with periodic cells that consist of membrane-split-ring resonators [14]. PVDF piezoelectric patch was attached to the membranes to convert the kinetic energy captured by the resonators into useful electric power. e fabricated structure exhibited simultaneous vibration attenuation and energy harvesting capabilities. Excellent and moderate vibration attenuation characteristics were reported in two bandgaps, i.e., 340-426 Hz and 460-500 Hz, respectively. e fabricated structure recovered approximately 0.5 µW across 200 kΩ at 348 Hz. Power was only measured and reported within the first bandgap, i.e., at 348 Hz. Also, a recent study by Li et al. reported a mechanical metamaterial structure for simultaneous vibration isolation and energy harvesting [15]. e structure consisted of beams with PVDF piezoelectric transducers to extract electric energy from the resonating cells. e metamaterial structure was able to attenuate vibrations and generate electric power, simultaneously, within 146-171 Hz frequency bandgap. Specifically, a peak power of 0.05 µW was measured across 1 MΩ within the frequency bandgap. e current state-of-the-art literature reveals that a growing number of researchers have recently attempted to use elastic mechanical metamaterials for simultaneous vibration attenuation and energy harvesting. e work in this area is still new and rapidly growing [16]. erefore, the work we present in this article is focused on investigating a magnetomechanical-based metamaterial structure for simultaneous vibration attenuation and energy harvesting. e presented structure consists of a combination of cantilever beams and permanent magnet-coil systems in order to achieve two main functions, simultaneously. ese functions are vibration attenuation and energy harvesting. Compared with other efforts and available literature, the metamaterial structure we present in this work is shown to generate, significantly, more electric power while maintaining its ability to attenuate undesired vibrations. Additionally, the presented structure requires lower load resistance to achieve optimum power transfer, compared with what has been presented in the literature [13][14][15]. is is an important feature of the presented structure because, typically, electronic circuits for sensors require an input current in the order of 10 mA [17]. As the load resistance increases the deliverable electric current decreases. It is desirable to achieve maximum power transfer at the lowest load resistance values [17]. e rest of the article is organized as follows: Section 2 presents the aspects of the proposed design and fabrication methodology of the metamaterial structure. Section 3 deals with modeling of the metamaterial structure. Section 4 is focused on characterization tests and experimental methods. Findings and results are presented and discussed in Section 5. Figure 1 displays the design and concept of the metamaterial vibration attenuation energy harvesting system presented in this work. e metamaterial structure, shown in Figure 1(a), consists of local resonators (unit cells). Each unit cell, shown in Figure 1(b), contains four angled cantilevers with fixed-free ends. A permanent magnet is placed in the center of each unit cell. e tip mass at the free end of each cantilever is made of copper coils wrapped around the free end of the cantilever as shown in Figure 1(c). e coils serve two purposes: first, the coils form additional mass to lower the resonant frequency of the vibrating cantilevers, thus lowering the frequency bandgap of the metamaterial structure; second, the magnet-coil system is used to extract electric energy from kinetic energy of the resonating cantilevers.

Design and Fabrication
When the metamaterial structure is subject to external vibrations, the vibrational energy is trapped into the local resonators, causing the cantilevers to resonate.
at is, when the driving frequency of vibration matches the resonant frequency of the local resonators, the vibration energy is transferred into kinetic energy using the resonating cantilevers, and therefore, this energy is localized inside the resonators. is results in frequency bandgaps. Consequently, these undesired vibrations are blocked from traveling across the ends of the metamaterial structure.
us, the first function of the proposed metamaterial vibration attenuation energy harvesting system is achieved. Concurrently, by integrating the coils and permanent magnets into the design of the local resonators, as shown in Figure 1(c), the kinetic energy of the resonating cantilevers is converted into useful electric power. Subsequently, voltage is induced in the coils wrapped around the free end of the cantilevers. erefore, the second function of the proposed system is achieved by converting the kinetic energy contained in these vibrations into useful electric power.

Theory
eoretical models of the proposed metamaterial vibration attenuation energy harvesting structure are developed next. First, COMSOL FEM software is used to investigate frequency bandgaps and modes of vibrations. Second, a theoretical model is developed using first principle and used to investigate wave attenuation capabilities of the presented metamaterial structure. Figure 2 shows the geometries and dimensions of the cantilever, i.e., the building block of the metamaterial structure presented in this work. First, a 3D CAD drawing of the unit cell was prepared using SolidWorks software, and the predetermined geometries and dimensions are shown in Figure 2. e CAD drawing was then imported into COMSOL finite-element multiphysics software to analyze the dynamic response of the metamaterial structure and obtain the dispersion diagram that depicts the relationship between the wave number, K, and frequencies for an infinite number of periodically arranged unit cells. at is, for an array of unit cells, at frequency intervals that meet Bragg's scattering condition, destructive interference occurs.

COMSOL FE Model.
is yields the desired frequency bandgap. e presence of the local resonators in the structure further enhances the frequency bandgap because, at the resonant frequency of the locally resonating structures, the kinetic energy contained in the external vibrations is transferred to the local resonators, causing them to resonate while creating a lower frequency bandgap in the dispersion diagram of the metamaterial structure.
In the COMSOL model, a fine mesh coupled with the choice of material properties (given in Table 1) and geometries (shown in Figure 2) is used to simulate the dispersion curve and obtain frequency bandgaps. e model assumes no damping effects and solves the governing equation of the metamaterial structure given by where ω and K d are frequency and dynamic stiffness matrix of the unit cell, respectively, K u and M u are stiffness and mass matrices of the unit cell, respectively. Here, q L , q R , and q int are displacement vector of the left side of the unit cell, displacement vector of the right side of the unit cell, and displacement vector of the internal nodes of the unit cell, respectively. Similarly, F L and F R are force vectors on the left and right sides of the unit cell, respectively. Upon finite element discretization, the unit cell contained 11,647 elements with a maximum size of 1.81 mm and a minimum size of 0.0775 mm.
To solve the governing equation of the system, standard Bloch-Floquet theory is imposed on the displacement, q, and force, F, vectors of each node, thus taking care of the periodicity condition of the unit cells given by where l u � 50 mm is the length of the unit cell as shown in Figure 2. Next, the described eigenvalue problem is solved for a given wave number in order to obtain the corresponding Eigen frequencies. A band structure is obtained by sweeping the wave vector along the irreducible Brillouin zone. is results in dispersion curves defining frequency bandgaps where oscillations are blocked and prevented from passing through. e tolerance was set to 1 × 10 −6 .

Analytical Model.
e metamaterial vibration attenuation energy harvester structure is represented as Euler-Bernoulli beam with periodic arrangement of local resonators attached to the beam. is is shown in Figure 3. Each local resonator consists of a cantilever with linear stiffness coefficient, k, and mass, m. In this work, transfer matrix method is used in modeling the wave propagation across the beam, assuming small displacement of the beam [18]. e general equation of motion of the Euler-Bernoulli beam is given as where y(x, t) is the dynamic deflection of the neutral axis at point x, E is the elastic modulus of the beam, I is the area moment of inertia, ρ is the density of the beam, and A is the cross-sectional area. In this work, it is assumed that where c is the flexural wave number given as For the n th unit cell, the mode shape function can be rewritten as Y n x′ � Q n cos cx′ + B n sin cx′ + C n cosh cx′ + D n sinh cx′ , Here, Q n , B n , C n , and D n are the unknown amplitudes describing the mode shape function [19].

Analysis of the Local Resonator (LR).
e equation of the n th local resonator in Figure 3 is described as where F(x n , t) is the force of the n th resonator developed as a result of interaction with the beam at the contact point as shown in Figure 3 and G n is the vertical displacement of the mass of the resonator, m, during the external excitation. e force, F(x n , t), is given by where the spring stiffness coefficient, k, of the local resonator cantilever is given as k � EWH 3 /4L 3 , where W and H are the width and height of the cantilever, and L is the overall length of the cantilever. Assuming the solution for the vertical displacement of the mass resonator of the form: where V n is the displacement of the local resonator (LR), and combining equations (7)-(9) yields the following equation: Substituting G n and y(x n , t) into equation (10) gives Rearranging equation (11), the vertical displacement, V n , can be written as Moreover, equation (8) can be rewritten as and substituting equation (12) into equation (13) yields e model for the four unconnected resonators (r � 1, 2, 3, 4) shown in Figure 3, at the interface with the beam can be derived as [20] m € G n, r � k r y x n , t − G n, r (t) , r � 1, 2, 3, 4, with the amplitude of displacement of the local resonators, V n , given by Consequently, the force acting on the beam due to these four resonators can be obtained from equations (13) and (14) and is given by: Next, the continuity conditions at the point where the resonators meet the beam, as shown in Figure 3, are imposed between the n th unit cell and the (n − 1) th cell. ese include displacement, slope, bending moment, and shear force and, respectively, given by: Resonator cantilever Shock and Vibration 5 Substituting equations (6) and (17) into equation (18) yields where ψ n � (Q n , B n , C n, D n ) T . Here, U and M are given by respectively, and From Bloch-Floquet's equation of period structures in the x direction, we have ψ n � e izl u ψ n−1 .
Further, by letting M −1 U � S and substituting equation (19) into equation (22) yield Equation (23) can be further rearranged into a standard eigenvalue problem described by where I is a [4 × 4] identity matrix.

Experimental Methods
A prototype of the metamaterial vibration attenuation energy harvesting structure was built and characterized experimentally. e fabrication process involved two main steps. First, the CAD design of the metamaterial structure was prepared using SolidWorks software and then converted into STL (STereoLithography) format in order to perform 3D printing of the prototype. e Ultimaker 3 Fused Deposition Modeling (3D FDM) printer was used to build the metamaterial structure using PETG (polyethylene terephthalate) filament. Second, 35 AWG enameled copper wires were winded around the free end of the cantilevers and neodymium iron boron (NdFeB) permanent solid magnets were manually integrated into the 3D printed metamaterial structure. Properties and dimensions of the fully assembled metamaterial vibration attenuation energy harvesting structure are given in Table 1. Figure 4 shows the fully fabricated and assembled metamaterial vibration attenuation energy harvester structure prototype. e experiment apparatus used to characterize the response of the metamaterial vibration attenuation energy harvester prototype is shown in Figure 5. e apparatus consists of a shaker table (VT-500, SENTEK DYNAMICS), power amplifier (LA-800, SENTEK DYNAMICS), vibration controller (S81B-P02, SENTEK DYNAMICS), accelerometers (PCB333B30 model, PCB Piezotronics), decade box (Global Specialties RDB-10), and PC. Characterization tests were focused on, simultaneously, measuring vibration transmissibility and the amount of electric power recovered from a unit cell. e electromagnetic shaker table was used to excite the metamaterial structure harmonically at predetermined frequencies and acceleration levels. An accelerometer was fixed, on top of the shaker table, and used to monitor the input waveform and acceleration level delivered to the metamaterial structure. When performing these experiments, the lower end of the metamaterial vibration attenuation energy harvester was secured on top of the shaker table as shown in Figure 5. A fixture and clamps were used to constrain and fix the top end of the metamaterial structure (not shown in Figure 5). A second accelerometer was fixed on the top end of the metamaterial structure and used to measure the vibration transmissibility of the metamaterial structure. Both input and output signals from the top and bottom accelerometers were measured and stored on the PC to analyze the vibration transmissibility. at is, the vibration transmissibility was measured by tracking the level of acceleration transmitted from the metamaterial structure relative to the input acceleration from the shaker table. Concurrently, the output voltage from a unit cell was measured across a load resistance and used to estimate the amount of recovered electric power.

Results and Discussion
First, we start by investigating wave propagation through the metamaterial structure. Using equation (24), for every selected frequency value, ω, the solutions of the wave vector, z,  are obtained. Figure 6 shows the results obtained from solving the eigenvalue problem given in equation (24). Here, a stopband exists when z is imaginary, while a pass-band exists when z is real [19]. As shown in Figure 6, the stopbands of the wave are marked in gray. Results reveal that there exist two main stopbands in the frequency range 1-1000 Hz. ese bands are occurring in the range of 121-310 Hz (lower frequency range) and 528-674 Hz (higher frequency range), respectively. ese two bands are the result of the beam itself and the local resonators (LR) [21]. at is, a higher frequency stopband is created as a result of vibration phase differences-in-phase and out-ofphase-between adjacent unit cells. A low-frequency stopband is typically generated by the LR stopband due to the attenuation effect of the resonators.
Next, to fully understand the behavior of the presented metamaterial vibration attenuation energy harvesting structure, we further analyze the structure using COMSOL FEM model described in Section 3.1 and experiment described in Section 4. e frequency bandgaps of the metamaterial vibration attenuation energy harvesting structure shown in Figure 7 are analyzed using COMSOL FEM. Periodic boundary conditions were imposed, and geometric and material properties were set similar to those measured and reported in Table 1. Figure 7 shows the dispersion curve of the metamaterial structure obtained using COMSOL model described in Section 3. Figure 7 suggests that there exists two major bandgaps within the frequency range of 100-1000 Hz. e first and second bandgaps are 218-247 Hz and 589-780 Hz, respectively.
To further explore these bandgaps, Figure 8 shows the corresponding mode shapes obtained using COMSOL at the lower and upper edges of the first and second bandgaps, i.e., 218, 247, 589, and 780 Hz, respectively. Clearly, within these bandgaps, vibrations are attenuated and trapped within the locally resonating cantilevers. Also, the fundamental resonant frequency of the unit cell (combined cantilever beams, tip coils, and central magnet) was estimated using COMSOL at approximately 224 Hz. e fundamental resonant frequency falls within the first bandgap of the metamaterial structure, i.e., 218-247 Hz. As shown in Figure 8, at the lower and upper edges of the first bandgap, the deformations of the unit cell are localized within the cantilevers. e vibrational energy is localized within the frequency bandgap and converted into kinetic energy of the resonating cantilevers. Also, Figure 8 suggests that maximum bending deformations occur within the first bandgap. erefore, it is expected that most of energy harvesting will occur within the first bandgap. is observation agrees with the vibration energy harvesting measurements that were performed in this work and discussed later in this article. Since most mechanical vibrations in nature occur at frequencies lower than 300 Hz [15], in this work, energy harvesting is targeted within this bandgap, i.e., 218-247 Hz.
Next, vibration transmissibility of the fabricated metamaterial structure was measured using experiment apparatus shown in Figure 5. In these experiments, a sinusoidal vibration sweep was performed in the range of 50-1000 Hz at an acceleration level of 0.5g m.s −2 . Figure 9 shows the measured vibration transmissibility of the metamaterial structure versus frequency. Clearly, the results from this experiment demonstrate two frequency bandgaps (marked in red). Results from COMSOL model simulations and analytical model (equation (24)) are displayed in Table 2. Results show a good agreement between models' predictions and measured data for both bandgaps. However, some discrepancy between models' predictions and measured data is evident. For example, the measured bandwidth of the first bandgap (205-257 Hz) is wider than the bandwidth of the COMSOL simulation-based first bandgap. A similar observation was made and reported by Chen et al. [14]. is may be attributed to some of the approximations made in the model simulations including absence of damping [14] or other experimental variations such as insignificant distortion in the metamaterial structure during the experiment. Similarly, the first bandgap obtained from equation (24) is a little wider than the bandgap predicted by COMSOL simulations and experiment as shown in Table 2. is is likely because of some of the simplifications and approximations made in the model (equation (24)) including absence of damping. Also, the analytical model does not factor in some of the design complexities including the flange attachment at the contact point with the beam. e model also assumes  small displacements of the transverse beam and ignores the mass contribution from the transverse beam. Nonetheless, both models predict similar behaviors and show similar trends compared with measured data. Furthermore, results from experiment show that the level of vibration attenuation, represented as the dip in the transmissibility curve in Figure 9, occurring in the first bandgap is substantially larger than the level of vibration attenuation observed in the second bandgap. is behavior is expected since, as discussed earlier and shown in Figure 8, the second bandgap experiences less deformations compared with the first bandgap, and most vibrational energy is localized within the first bandgap where the fundamental resonant frequency is located, i.e., 224 Hz. e ability of the fabricated metamaterial structure to scavenge electric power while concurrently attenuating undesired vibrations is examined next. Figure 10 shows electric power harvested from one cantilever located in the  bottom central unit cell of the metamaterial structure. In this experiment, the metamaterial structure was driven harmonically at 0.5 g·m.s −2 and fixed frequency (selected within the first bandgap, i.e., 223 Hz) while the load resistance was swept in the range of 1-1000 Ω (Figure 10(a)). e same experiment was repeated at a fixed frequency selected from the second bandgap, i.e., 615 Hz, as shown in Figure 10(b). Results from these experiments confirm that most of the vibrational energy is localized within the first bandgap, causing the cantilevers to locally resonate. Consequently, within the first bandgap, kinetic energy in these resonating cantilevers is effectively converted into electric power through the coils at the free end of these cantilevers. Maximum power outputs of, approximately, 2.5 µW and 0.6 nW, were measured at 223 Hz (first bandgap) and 615 Hz (second bandgap), respectively. us, results demonstrate that both vibration transmissibility ( Figure 9) and energy harvesting ( Figure 10) characteristics of the fabricated metamaterial structure are more prominent in the first frequency bandgap. Also, Figure 10 reveals that the maximum power output occurs at optimum load resistance of 15 Ω.
It is also worth noting that the improvement in energy harvesting within the first bandgap (205-257 Hz) is accompanied by stronger vibration attenuation. at is, vibration attenuation and energy harvesting characteristics of the metamaterial structure are coupled. Stronger vibration attenuation leads to enhanced vibration energy harvesting capabilities. is is further demonstrated in Figure 11 which shows the output electric power versus frequency. In this experiment, the metamaterial structure was swept harmonically at 0.5g m.s −2 in the frequency range of 50-500 Hz. e load resistance was fixed at the optimum value, i.e., 15 Ω, in order to maximize power generation. It is clear that, within the bandgap (205 Hz-257 Hz), the large dip in the transmissibility curve (Figure 9) corresponds to the maximum power generation ( Figure 11). Also, Figure 11 reveals that the maximum power generated within the frequency bandgap can reach up to, approximately, 5.2 µW at 245 Hz.
Next, performance metrics of the presented metamaterial vibration attenuation energy harvesting structure are compared to the state-of-the-art work reported in the recent literature. is is shown in Table 3. e fabricated metamaterial vibration attenuation energy harvesting structure performed very well compared with other efforts recently reported in the literature. For example, the metamaterial structure was able to attenuate undesired vibrations at the, typically, targeted frequencies, i.e., less than 300 Hz. Concurrently, the presented metamaterial structure produced significantly more electric power compared with what has been reported in the literature. Table 3 also shows another advantage of the work presented in this article. e optimum load resistance to maximize electric power was significantly lower compared with what has been reported in the literature. is low impedance guarantees higher electric current that is desirable to operate electronic circuits for sensors, as discussed earlier in this article. Although our initial prototyping efforts produced a metamaterial structure with comparable metrics to what has been reported in the literature, it should be noted that optimization of this metamaterial dual-purpose structure is beyond the scope of work reported in this article. e work reported in this article is focused on design, proof-of-concept experiments, and characterization tests of the metamaterial structure rather than optimization. Several design improvements can be made to improve the overall performance of the fabricated metamaterial dual-function structure.
is includes optimization of geometries and dimensions of the cantilever beams, tilt angle of cantilever beams, number of unit cells, and size and location of central magnets, etc. is will require developing a detailed and sophisticated coupled model of the magnetomechanical metamaterial vibration attenuation energy harvesting structure and, therefore, will be the focus of future work.

Conclusions
e work presented in this article is focused on the design and characterization of a magnetomechanical-based metamaterial structure for simultaneous vibration attenuation and energy harvesting. e structure consists of periodically arranged unit cells that are made of a combination of cantilever beams and permanent magnet-coil systems. Each unit cell contains four angled cantilevers with fixed-free ends. A permanent magnet is placed in the center of each unit cell. e tip mass at the free end of each cantilever is made of copper coils that serve two purposes: lowering the frequency bandgap of the metamaterial structure and extracting electric energy from kinetic energy of the resonating cantilevers as they move near the central permanent magnet. A prototype of the metamaterial dual-function structure has been fabricated using additive manufacturing and then assembled manually. Models of the structure have been developed and validated against the experiment.
Results show good agreement between model simulations and experiment. Two frequency bandgaps, i.e., 205-257 Hz and 587-639 Hz, respectively, have been measured. Within these two bandgaps, vibrations are completely attenuated. Results suggest that the level of vibration attenuation in the first bandgap is substantially larger than the level of vibration attenuation observed in the second bandgap. Mode shapes obtained using COMSOL FEM model reveal that the second bandgap is accompanied by less bending deformations compared with the first bandgap and most vibrational energy is localized within the first bandgap where the fundamental resonant frequency is located, i.e., 224 Hz. Moreover, the ability of the fabricated metamaterial structure to harvest electric power while synchronously attenuating undesired vibrations has been inspected. Electric power outputs of, approximately, 2.5 µW and 0.6 nW were measured at 223 Hz (first bandgap) and 615 Hz (second bandgap) at optimum load resistance of 15 Ω. ese results confirm that most of the vibrational energy is localized within the first bandgap, causing the cantilevers to locally resonate and kinetic energy in these resonating cantilevers to be, effectively, converted into useful electric power through the coils at the free ends of these cantilevers. Results also show that vibration attenuation and energy harvesting characteristics of the metamaterial structure are coupled. Stronger vibration attenuation has led to enhanced vibration energy harvesting capabilities. Frequency sweep measurements at optimum load resistance of 15 Ω reveal that maximum power generated within the first frequency bandgap reaches, approximately, 5.2 µW at 245 Hz.
Compared with the state-of-the-art and reported literature, the metamaterial structure we present in this work has shown an order of magnitude improvement in electric power generation, at substantially lower optimum load resistance, while maintaining the ability to attenuate undesired  vibrations in a frequency bandgap of 205-257 Hz. Future work will focus on optimization of the metamaterial structure to maximize energy harvesting while maintaining vibration attenuation capabilities.

Data Availability
Data reported in this article are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this article.

Authors' Contributions
Winner Anigbogu was involved in the data collection, formal analysis, software, validation, investigation, visualization, and formal analysis. Hamzeh Bardaweel was responsible for the conceptualization, funding acquisition, methodology, project administration, resources, software supervision, validation, visualization, and writing, reviewing, and editing of the original draft.