Layered Metamaterial Beam Structures With Local Resonators for Vibration Attenuation: Model and Experiment

Abstract

Layered metamaterial beam structures are gaining attention in a variety of fields including vibration attenuation and energy harvesting. Exhaustive research on single-beam metamaterial vibration attenuation structures using local resonators exists in literature. Moreover, there are recent attempts at modelling double-layered beams with different kinds of constraints. The double-layered beam models in literature are limited to simple beams and not extended to metamaterials with local resonators. This article is focused on developing a design criterion and a modelling platform for layered metamaterial structures with multiple beams and local resonators for vibration isolation. The model is developed using Euler-Bernoulli beam equations, superposition of mode shapes and Galerkin methods. A prototype layered metamaterial structure is fabricated and characterized experimentally. The prototype consists of horizontal beams, local resonators forming unit cells, and vertical beams linkages. Each local resonator consists of cantilevers with tip masses. Results show good agreement between model and experiment. Two major bandgaps are observed at 190–410 Hz and 550–710 Hz. Results reveal that the low frequency bandgap can be further reduced through the design of the local resonators. Results also show that alternating the length of the local resonators causes a shift in the first frequency bandgap. An increase in the number of local resonators opens up extra frequency bandgaps at lower frequencies with the drawback of reducing the depth in vibration transmissibility. Moreover, the higher frequency bandgaps are mostly affected by the horizontal beams. An increase in the length of the horizontal beams, while the number and design of the local resonators are fixed, broadens the second frequency bandgap and shifts it to lower values.

1 Introduction

Recently, mechanical metamaterials have become the centerpiece of many research studies and engineering applications (Xia et al., 2020). This is mainly because of their unique properties including the presence of frequency bandgaps influenced by the local resonators (Wang et al., 2016; Jiang, He, 2017). This is particularly of interest for many engineering applications including wave guiding, vibration attenuation and, more recently, dual-purpose vibration suppression and energy scavenging (Matlack et al., 2016; Reichl and Inman, 2017; Casablanca et al., 2018). Frequency bandgaps created in these metamaterial structures originate from the ability of the local resonators to absorb kinetic energy from the oscillations passing through the metamaterial structure (Li et al., 2017). When the frequency of the external vibration source aligns with the resonant frequency of the local resonators, energy from these oscillations is transferred to the local resonators (Krödel et al., 2015; Li et al., 2017). Typically, these local resonators are either attached to the surface of the metamaterial structures or embedded inside (Fang et al., 2006; Chen and Wang, 2014; Krödel et al., 2015). Moreover, the Bragg’s scattering in the periodic structures leads to generation of frequency bandgaps (Kushwaha et al., 1993). The bandgaps generated by Bragg’s scattering are influenced by the wavelengths of the Bragg scattering mechanism, which is limited by the design of the beam of the mechanical metamaterial structure (Krödel et al., 2015). The addition of local resonators to the structure makes it possible to generate bandgaps at frequencies much lower than the Bragg’s frequencies (Liu et al., 2000; Liu et al., 2005). In the case of vibration attenuation, when the mechanical metamaterial is subject to vibrations by an external vibration source, it has the ability to suppress these oscillations, at low frequencies, that are closely aligned with the resonant frequency of the local resonators (Sugino et al., 2017).

Exhaustive literature has already looked at design criteria and characterization of metamaterial structures that are made of single beam structures with local resonators (Yu et al., 2006a; Yu et al., 2006b; Wang et al., 2013; Wang and Wang, 2013; Xiao et al., 2013; Zhou et al., 2017). For example, the work by Liu et al. was among the first few studies to look at metamaterial structures and local resonators (Liu et al., 2000). Their pioneer seminal work showed that localized resonant structures exhibit valuable characteristics such as negative elastic constants, leading to generation of low frequency bandgaps. The transfer matrix method was then used by Yu et al. for characterization of a single beam metamaterial structure with resonators that are made of rubber and copper rings (Yu et al., 2006a; Yu et al., 2006b). The results from their model simulations were validated experimentally. Xiao et al. built on the analytical model of the single beam metamaterial design developed by Yu et al. and showed that resonance and Bragg’s bandgaps exist as a result of local resonance effect and periodicity of unit cells (Xiao et al., 2013). Further studied the effects of different design parameters on the generated frequency bandgaps Liu et al. (2007). The results from their work revealed that, although a large local resonator generates wide attenuation, a large number of small local resonators spread across the beam gives enhanced attenuation spectra.

Other enhanced designs and studies of single beam-based metamaterial structures have also been developed (Chen et al., 2011; Xiao et al., 2012; Wang et al., 2013; Wang and Wang, 2013; Hu et al., 2021). For instance, a design of a beam with spring-mass local resonators sandwiched inside the beam was proposed by Chen et al. (2011). Results showed that the amplitude of waves was attenuated near the resonance frequency of the local resonators. Additionally, their results revealed that bandgaps were influenced mostly by the size of the internal mass rather than the spring constant of the resonator. To create more flexible bandgaps, designed an array of local resonators with varying resonant frequencies Xiao et al. (2012). Their work showed that, compared to a beam with a single array of local resonators, a beam with varying array of local resonators can generate broader bandgaps above and below Bragg’s bandgap. Used a combination of both vertical springs and oblique springs to investigate the effect of negative stiffness introduced by the oblique springs on single beam-based metamaterial structures Zhou et al. (2017). Their results showed that adding two oblique springs to the single beam-based metamaterial structure introduced negative stiffness which effectively produced attenuation at lower frequencies. More recently, studied the effect of nonlinearity that was purposefully introduced into the design of a single beam-based metamaterial structure Xia et al. (2020). This was done by adding bistable attachments to the metamaterial structure. Their work showed that nonlinear vibrations of the bistable resonator attachments resulted in wider frequency bandgap compared to the case when only linear local resonators were considered.

All the aforementioned studies have focused on single-beam based metamaterial structures with local resonators. On the other hand, few recent studies have looked at modelling layered beams. For example, Oniszczuk investigated two beams connected by a Winker elastic layer (Oniszczuk, 2000). The work used Bernoulli-Fourier method to show that the free vibration of a simply supported double-beam system is similar to a double-string system. Used both Fourier and Laplace approach to develop a closed form solution for double infinite Euler-Bernoulli beams with harmonic loads applied Li et al. (2019). Douglas et al. treated a system of three elastic-viscoelastic beams as two non-identical Euler-Bernoulli beams with a viscoelastic layer between them (Douglas and Yang, 1978). The method was limited to a specific fixed-free boundary condition and an external load applied only at the free end. Extended Euler-Bernoulli’s approach to arbitrary boundary conditions and applied loads Vu et al. (2000). Then applied a change of variable method in decoupling the set of fourth-order differential equations from the modal solution of double Euler-Bernoulli beams with moving loads Abu-Hilal (2006). Recently demonstrated the use of layered metamaterial for simultaneous vibration isolation and energy harvesting Li et al. (2017). Anigbogu et al. presented a unique metamaterial design using 3D printed layered beams to achieve vibration attenuation (Anigbogu and Bardaweel, 2020). Their results showed that layered metamaterial beams can achieve vibration attenuation at low frequencies that are commonly encountered in nature.

A key observation from these recent studies is that these efforts have mainly focused on studying and investigating either single beam-based metamaterial structures or double-beams with different support systems. None of these studies has tackled the issue of developing a general approach to modelling layered metamaterial beams with local resonators. The work presented in this article is motivated by the growing interest in layered metamaterial beams with local resonators. Recently, these layered metamaterial beams with local resonators have found numerous potential engineering applications including vibration attenuation and energy scavenging (Li et al., 2017; Anigbogu and Bardaweel, 2020). To the best knowledge of the authors, the work presented here is the first to tackle the issue of developing a modelling platform for layered metamaterial structures with multiple beams and local resonators. The work is focused on developing design criteria and modelling platform for layered metamaterial structures with multiple beams and local resonators. The model is developed using Galerkin method and superposition of mode shapes to solve an application of Euler-Bernoulli equation to layered-metamaterial beam structures. A prototype of the layered metamaterial structure with local resonators is fabricated and characterized experimentally. The model developed in this work is then validated against measured data. Additionally, in this work design criteria and the effects of core design parameters on the frequency bandgaps generated by the layered metamaterial structure are studied. Section 2 describes the theory and mathematical model developed for the layered metamaterial structure. Section 3 introduces the experimental methods, while Section 4 describes the main findings and discussions of these results. Conclusions and summary of the work are articulated in Section 5.

2 Theory and Model Development

The overall representation of the layered metamaterial vibration attenuation structure is shown in Figure 1A. The metamaterial structure is comprised of J-layered horizontal beams, the local resonators per unit cell and the vertical beams linking the horizontal beams together. The local resonators are made of cantilevers with tip masses attached to them as shown in Figure 1A. Each unit cell has four local resonators linked to the same point as shown in Figure 1B. For simplicity we represent the unit cell as one single local resonator force acting on the same point, in tandem with the approach adopted by Wang et al. in (Wang and Wang, 2013). The effects of the vertical beams are factored into the governing equations as spring forces and damping forces following the relative velocities of the beams as shown in Figure 2. The free body diagram of the local resonators within a unit cell is shown in Figure 3. An approximate model of the transmissibility of the layered metamaterial structure is obtained using Euler-Bernoulli’s equations and Galerkin method (Zhou et al., 2017).

FIGURE 1

FIGURE 2

FIGURE 3

In order to analyze the dynamic behavior of the layered metamaterial vibration attenuation structure when subject to external vibrations, the Frequency Response Function (FRF) approach is used to determine the transmissibility of the device. The Galerkin method is used to analyze the beams’ response to external vibrations. The length of each horizontal beam is , where α is the number of unit cells, and is the length of each unit cell in the design adopted in this work (Figure 1B). The transverse vibration of a simple beam with attachment and a vibration source at a point is given by the standard Euler Bernoulli equation (Zhou et al., 2017):where , and are the transverse displacement, Young’s modulus of elasticity and the cross-sectional moment of inertia of the beam, respectively. Moreover, and are the density of the material, and the cross-sectional area of the beam, respectively. The expressions and represent the force from the vibration source, the delta Dirac function and the force from the attached local resonators, respectively. The force, , from the vibration source only applies to the first beam. In Eq. 1, represents the position of the local resonator, while is the total number of local resonators on each horizontal beam. This Euler Bernoulli equation shown in Eq. 1 is used to analyze the displacement of each of the horizontal beams. Applying Galerkin method and the superposition of modes, the transverse displacement of horizontal beam, is given by (Zhou et al., 2017):

In Eq. 2, is the generalized displacement of the horizontal beam at the mode, and is the trial function of the horizontal beam which is taken to be the mode shape of the horizontal beam at the mode. The index is the total number of layered horizontal beams.

Next, we extend Eq. 1 to factor in the connecting vertical beams as shown in Figures 1, 2. This yields the set of equations representing each of the horizontal beams from the base beam to the top, respectively. A generalized form of this set of equations from the lowest horizontal beam, i.e., , to middle horizontal beam, i.e., , till the topmost beam, i.e., , for any number of horizontal beams constituting the layered metamaterial structure yields Eqs 3–5 below (Abu-Hilal, 2006):where and are the lowest horizontal beam and topmost horizontal beam, respectively, and the middle section beams are indexed with . In Eq. 3, is the force exerted by the vibration source, represented as a sinusoidal input. Also is the forcing frequency of vibration source, is the amplitude of the external vibration force and represents the stiffness of the vertical beams given by (Sharma and Tyagi, 2013):where , and represent the width, thickness and length of the vertical beams respectively. The damping of each vertical beam, , and the equivalent mass, , are given by (Waichi et al., 2010):

Grouping the impact of the vertical beams as forces acting at specific points on the horizontal beams, Eq. 3 can be rewritten in the following generalized forms:where and are the forces due to the vertical beams acting on each horizontal beam (from the lowest horizontal beam to horizontal beam) and the topmost horizontal beam, respectively, and represents the point of contacts on each horizontal beam. Using Eq. 2, these forces, which are due to the interlinking vertical beams, can be further developed into the following generalized set of equations:

Also, the forces from the attached local resonators and the internal damping force are given as:where , and are the internal damping force, the damping ratio and the mass of the local resonator, respectively. is the vertical displacement of the local resonator.

Substituting Eq. 2 into Eq. 5 yields the following set of equations, respectively:

Multiplying Eq. 9 by the trial function , integrating over the length of each of the horizontal beams, i.e., from 0 to and applying orthogonality of mode shapes, yields the following transformation:

Here, the transverse displacements, , in the forces , and in Eq. 10, represented in Eqs 6–8, are also affected by the integration and application of orthogonality condition in Eq. 10. Also, and are given as:

In Eq. 10 the impact of damping on the layered metamaterial structure is considered. Here, and are the horizontal beam’s damping coefficient and the damping ratio of the horizontal beam in mode, respectively. The mode functions of the fixed-fixed beams, , are assumed to be the same as the trial function, , and given as (Zhou et al., 2017):

Moreover, the wave number is given by:

Furthermore, using Newton’s second law, the movement of each local resonator on each horizontal beam (including the horizontal beam) can be described as:

Solving Eqs 10, 14 simultaneously using Runge-Kutta numerical method, the generalized displacements and the local resonator displacement are obtained. The obtained resonator displacement is then substituted into Eq. 2, in order to obtain the displacement of the horizontal beam. This Galerkin solution is truncated at N = 4. Transmissibility is taken as the output displacement at the center of the most upper horizontal beam, i.e. , relative to the input vibration at the center of the lowest horizontal beam, i.e., vibration source. This is obtained by substituting the generalized displacement result, , obtained from solving Eqs 10, 14 simultaneously, into Eq. 2. The transmissibility , is then given as:

3 Experimental Methods

In this section experimental methods are described. Mainly, a brief description of fabrication and assembly of the layered metamaterial vibration attenuation structure as well as experimental setup used for characterization are discussed. Experimental work was focused on measuring the transmissibility of the layered metamaterial vibration attenuation structure. Model simulations were then validated against these experimental measurements.

The layered metamaterial structure was designed using Solidworks CAD software, then changed into an STL (STereoLithography) format that is suitable for 3D printing. PETG (polyethylene terephthalate) was used as the printing filament in the Ultimaker 3 Fused Deposition Modeling (3D FDM) printer. The local resonators (cantilever beams) were designed in a way that allowed insertion of additional masses at their tips. The tip-masses were created separately and then slotted into the tips of the local resonators. These additional tip-masses were vital in lowering the resonant frequencies of local resonators to less than 300 Hz. Figure 4 shows a prototype of the layered metamaterial vibration attenuation structure fabricated in this work. For this work, a layered metamaterial beam structure with four horizontal beams was chosen. This choice was informed by the need for the layered metamaterial structure to fit on the 3D printer bed used for fabrication. The dimensions and specifications of the manufactured layered metamaterial structure are outlined in Table 1.

FIGURE 4

The set-up of the experiment is shown in Figure 5. First, an external vibration waveform was generated by the controller (S81B–P02, SENTEK DYNAMICS). This waveform was then passed through the power amplifier (LA-800, SENTEK DYNAMICS) for amplification, before reaching the shaker table (VT-500, SENTEK DYNAMICS). The electromagnetic shaker table was used to apply the external vibration input to the layered metamaterial structure that was fixed on top of the shaker table. The shaker table was operated at predetermined frequency range of 50–700 Hz, and acceleration level of 0.5 g [m s⁻²]. The sides of the layered metamaterial vibration attenuation structure were held with a fixture and clamps. Two accelerometers (PCB333B30 model, PCB Piezotronics) were used to capture the input signal from the shaker table and the output signal from the top of the layered metamaterial structure. The first accelerometer was placed on top of the shaker table to capture the input signal data, while the second accelerometer was positioned on top of the layered metamaterial structure to capture the output signal. The input and output signals were both monitored using Crystal Instrument’s Engineering Data Management (EDM) software on a PC. The transmissibility was analyzed by comparing the output acceleration from the top-most beam of the metamaterial structure to the input acceleration from the vibration source, i.e., shaker table.

FIGURE 5

4 Results and Discussion

The results from the FRF analysis were obtained from Eq. 15 as the Root Mean Square ratio between the output displacement (at the center of the upper beam) and the input displacement (from the vibration source) (Zhou et al., 2017). The parameters used in this study are shown in Table 1. For the purpose of model validation and design investigation, a layered metamaterial structure made of 3 unit cells per beam, three vertical columns between beams, and a total of four horizontal beams was manufactured (Figure 4). In this study modes beyond the fourth modes were truncated because higher modes are unnecessary as have been shown in literature (Li et al., 2017). For example, the 5th mode of the unit cell is 1,469 Hz, which is much higher than the frequency range of interest in the current study (below 700 Hz). When frequencies much higher than 700 Hz are considered, the value of N should be increased to capture those higher modes (Zhou et al., 2017).

4.1 Model Simulation

First, we start by examining the displacement amplitude of the topmost local resonator, i.e., cantilever beam, and the lowest local resonator. Displacement response was obtained using Eq. 14 and results are shown in Figure 6. Results reveal peak displacement at approximately 220 Hz. This peak corresponds to the first fundamental resonant frequency of the unit cell which was also determined and verified using COMSOL simulations and found at, approximately, 224 Hz (Anigbogu and Bardaweel, 2020). The COMSOL simulation of the unit-cell was done using Bloch-Floquet’s periodic boundary condition. The unit-cell was modelled as a fine mesh comprised of 11,647 elements with a maximum size of 1.81 mm and a minimum size of 0.0775 mm (Anigbogu and Bardaweel, 2020).

FIGURE 6

Next, the transmissibility of the metamaterial structure was estimated using Eq. 15 and the results are shown in Figure 7. Results reveal two main bandgaps of interests, i.e., below 700 Hz. Frequencies higher than 700 Hz were not considered in this work because most of mechanical vibrations in nature occur at frequencies well below the 700 Hz frequency limit adopted in this work (Li et al., 2017). The first bandgap expands from 190 to 410 Hz with its lowest transmissibility dip occurring at approximately 245 Hz. The second frequency bandgap takes place between 550 and 710 Hz with its lowest dip at approximately 600 Hz.

FIGURE 7

It is worth noting that the level of vibration attenuation, shown in Figure 7, observed within the first bandgap is significantly larger than the level of vibration attenuation observed in the second bandgap. This is likely because most of the vibrational energy is localized within the first bandgap. To further investigate this behavior, the input and output displacements of the layered metamaterial structure at frequency within the first bandgap, , are obtained using Eqs 2, 10. Results are shown in Figure 8A. Similarly, output and input displacements are obtained at higher frequency within the second bandgap, i.e., 600 Hz using Eq. 10 and Eq. 2 and results are shown in Figure 8B. Figures 8A,B reveals that, at those selected frequencies, the output displacement amplitudes are suppressed compared to the input displacement amplitudes. Nonetheless, near the fundamental resonant frequency of the cantilever, i.e., 224 Hz the vibrations are trapped in the local resonators. This localization of vibrational energy results in significant reduction in the displacement of the most-upper beam (output) compared to the bottom beam near the vibration source (input) as shown in Figure 8A. This observation explains the substantial dip in transmissibility within the first bandgap compared to the second bandgap as shown in Figure 7.

FIGURE 8

4.2 Experimental Validation

Next, model simulations are validated against experimentally measured data and the results are shown in Figure 9. Measured data shown in Figure 9 were obtained using experiment setup shown in Figure 5. Results show a good agreement between model simulations (Figure 7) and measured data (Figure 9). While results from both model simulations and experimental data exhibit alike characteristics and trends some discrepancy is present. For example, model simulations predict a second bandgap of 550–710 Hz compared to the measured one, i.e., 587–639 Hz. This mismatch between model simulations and measured data may be attributed to a few factors. First, the mismatch between model simulations and experimental data may be attributed to the protrusions and cuts that exist in the fabricated layered metamaterial structure and were not factored into the model simulations. The presence of some extruded cuts and holes in the 3D printed layered metamaterial structure, coupled with defects accrued from assembling the tip masses and the 3D printed layered metamaterial structure may contribute to the slight mismatch between model and experiment. Second, printing parameters such as temperature, extrusion rate and speed may affect the quality of 3D printed structures (Mueller and Shea, 2000; Pilipović et al., 2009; Singh, 2011; Mueller et al., 2015; Li et al., 2017; Reichl et al., 2021). These are factors that cannot be readily accounted for in the developed model. For example, Reichl et al., found that changes in 3D printing temperature may affect the damping properties of cantilever vibration absorbers and therefore likely to adjust their resonant frequencies (Reichl et al., 2021). In this work, similar to the approach taken by Xia et al. fixed modal damping ratio was assumed for simplicity (Xia et al., 2020). Additionally, the slight mismatch between model simulations and experiment may be attributed to the experiment setup itself. For instance, deviation between dynamic model predictions and experiment was reported by Berdy et al. (2014), Lee et al. (2016), Dong et al. (2017), Dhote et al. (2018) and it was attributed to the likely presence of a small tilt in the experimental rig which may lead to multi-direction vibrations and damping that are not accounted for in the dynamic model.

FIGURE 9

4.3 Further Studies on the Dynamic Behavior of the Layered Metamaterial Structure

Next, design criteria and effects of different design parameters on the response of the layered metamaterial vibration attenuation structure are investigated using the developed model. Figure 10 shows the effect of the number of local resonators (unit cells) per beam on the transmissibility of the layered metamaterial structure. It is worth noting that each unit cell has a length, and one local resonator. The nominal original horizontal beam is comprised of three-unit cells, implying three local resonators per beam. Therefore, an increase in the number of local resonators (unit cells) per horizontal beam from 3 to 5 yields an increase of the overall length of all the horizontal beams (L) from 150 to 250 mm (Figure 10A). Figure 10B shows that an increase in the number of local resonators (unit cells) per beam leads to new frequency bandgaps opening up at lower frequencies. For instance, increasing the number of resonators from 3 to 5 introduces a new bandgap at lower frequencies i.e., 98–128 Hz. The same pattern of an extra bandgap opening up at lower frequencies persists when the local resonators increased to 7 and then 10. This is likely because with an increase in the number of unit cells more resonant peaks are expected to appear at lower frequencies (Zhou et al., 2017). Also, an increase in the number of local resonators means that the length of the horizontal beams and the forces acting on these beams are changed while the geometry of the local resonators remained the same. Therefore, in Figure 10B, the new bandgaps can be also attributed to Bragg’s scattering that is resulting from intrinsic changes on the horizontal beam. It is also worth noting that the depth of the bandgap decreased when the number of local resonators increased. This reduction in the depth of the transmissibility is likely due to the fact that, with an increase in the number of local resonators, there is also a concurrent increase in the number of vertical beams holding the horizontal beams. The increase in the number of vertical beams makes the overall layered metamaterial vibration attenuation structure stiffer; resulting in smaller transmissibility depths as shown in Figure 10B. Therefore, a tradeoff between having more bandgaps at lower frequencies and reduction in the transmissibility dips exists.

FIGURE 10

To further investigate the impact of the local resonators on the first frequency bandgap, the lengths of the resonators are varied and results are shown in Figure 11. Results show that increasing the length of the local resonators causes the first bandgap to move to lower frequencies. Results also suggest that when the length of the local resonators is increased above the nominal length, , the high frequency bandgap, i.e., second frequency bandgap, stays almost unchanged. For instance, when the length of the local resonators is increased from to the first bandgap shifts from 118–250 Hz–93–200 Hz while the second frequency bandgap remains the same i.e., 440–660 Hz. This is likely because an increase in the length of the local resonator reduces the stiffness and subsequently its resonant frequency; causing the first bandgap to shift to lower values. The second bandgap remains unchanged because it is mostly affected by the dimensions of the horizontal beams which remained fixed. This is because higher frequency bandgaps are mostly affected by vibration phase differences between adjacent unit cells (Wang and Wang, 2013). When the dimensions of the horizontal beam are not changed the higher frequency bandgap is minimally impacted. The transmissibility dip of the first bandgap, when the length of the local resonator is changed to , falls within the region of the new resonant frequency, , of the local resonator.

FIGURE 11

Next, the impact of a change in the horizontal beams’ characteristics on the frequency bandgaps is further investigated. This is done through studying the effect of variations in the length of the unit cell beam, , while the number of resonators is kept the same, i.e. α = 3 and their dimensions unaltered (Figure 12A). Transmissibility results are shown in Figure 12B. Here, the total length of each of the four horizontal beams, , is increased from 150 mm (when ), to 180 mm (when ). A variation in the length of the unit cell’s horizontal beam from the original size, , results in noticeable changes in the second frequency bandgap (Figure 12B). When the length of the unit cell of the horizontal beams is changed from 60 to 65 mm, as shown in Figure 12B, the position of the first bandgap experiences slight changes. That is, as shown in Figure 12B, only the edges of the first bandgap moved slightly from 150–380 Hz to 138 HZ–360 Hz. Moreover, while altering the horizontal beams’ dimensions slightly affects the width of the first frequency bandgap, the depth of its transmissibility, which is primarily created by the local resonators, remains unchanged. On the other hand, a closer look at the second bandgap, one can notice that there are more prominent changes. Here, the second bandgap becomes broader and shifts to the left with an increase in the length of the horizontal beams. This is expected because increasing the length of the unit cell, , makes the total length of the horizontal beam longer and hence significantly impacts the vibration phase differences between adjacent unit cells. Thus, altering the horizontal beams’ dimensions has a strong impact on the second bandgap, i.e., the higher frequency bandgap.

FIGURE 12

5 Conclusion

References

• 1

  Abu-HilalM. "Dynamic Response of a Double Euler–Bernoulli Beam Due to a Moving Constant Load." J. Sound Vibration Vol. 297No. 3 (2006): ppp. 477–491. 10.1016/j.jsv.2006.03.050

  • CrossRef
  • Google Scholar

• 2

  AnigboguW.BardaweelH. (2020). A Metamaterial-Inspired Structure for Simultaneous Vibration Attenuation and Energy Harvesting. Shock and Vibration2020, 4063025. 10.1155/2020/4063025

  • CrossRef
  • Google Scholar

• 3

  BerdyD. F.ValentinoD. J.PeroulisD. "Design and Optimization of a Magnetically Sprung Block Magnet Vibration Energy Harvester." Sensors Actuators A: Phys. Vol. 218 (2014): ppp. 69–79. 10.1016/j.sna.2014.06.011

  • CrossRef
  • Google Scholar

• 4

  CasablancaO.VenturaG.GarescìF.AzzerboniB.ChiaiaB.ChiappiniM.et al (2018). Seismic Isolation of Buildings Using Composite Foundations Based on Metamaterials. J. Appl. Phys.123 (No. 17), 174903. 10.1063/1.5018005

  • CrossRef
  • Google Scholar

• 5

  ChenJ. S.SharmaB.SunC. T. "Dynamic Behaviour of Sandwich Structure Containing Spring-Mass Resonators." Compos. Structures Vol. 93No. 8 (2011): ppp. 2120, 2125. 10.1016/j.compstruct.2011.02.007

  • CrossRef
  • Google Scholar

• 6

  ChenY.WangL. (2014). Periodic Co-continuous Acoustic Metamaterials with Overlapping Locally Resonant and Bragg Band Gaps. Appl. Phys. Lett.105 (No.19), 191907. 10.1063/1.4902129

  • CrossRef
  • Google Scholar

• 7

  DhoteS.YangZ.BehdinanK.ZuJ. "Enhanced Broadband Multi-Mode Compliant Orthoplanar Spring Piezoelectric Vibration Energy Harvester Using Magnetic Force." Int. J. Mech. Sci. Vol. 135 (2018): ppp. 63–71. 10.1016/j.ijmecsci.2017.11.012

  • CrossRef
  • Google Scholar

• 8

  DongG.ZhangX.XieS.YanB.LuoY. "Simulated and Experimental Studies on a High-Static-Low-Dynamic Stiffness Isolator Using Magnetic Negative Stiffness Spring." Mech. Syst. Signal Process. Vol. 86 (2017): ppp. 188, 203. 10.1016/j.ymssp.2016.09.040

  • CrossRef
  • Google Scholar

• 9

  DouglasB. E.YangJ. C. S. "Transverse Compressional Damping in the Vibratory Response of Elastic-Viscoelastic-Elastic Beams." AIAA J. Vol. 16No. 9 (1978): ppp. 925–930. 10.2514/3.7595

  • CrossRef
  • Google Scholar

• 10

  FangN.XiD.XuJ.AmbatiM.SrituravanichW.SunC.et al "Ultrasonic Metamaterials with Negative Modulus." Nat. Mater Vol. 5No. 6 (2006): ppp. 452–456. 10.1038/nmat1644

  • CrossRef
  • Google Scholar

• 11

  HuG.C. M. AustinA.SorokinV.TangL. (2021). Metamaterial Beam with Graded Local Resonators for Broadband Vibration Suppression. Mech. Syst. Signal Process.146, 106982. 10.1016/j.ymssp.2020.106982

  • CrossRef
  • Google Scholar

• 12

  JiangT.HeQ. (2017). Dual-Directionally Tunable Metamaterial for Low-Frequency Vibration Isolation. Appl. Phys. Lett.110 (No.2), 21907. 10.1063/1.4974034

  • CrossRef
  • Google Scholar

• 13

  KrödelS.ThoméN.DaraioC. "Wide Band-Gap Seismic Metastructures." Extreme Mech. Lett. Vol. 4 (2015): ppp. 111–117. 10.1016/j.eml.2015.05.004

  • CrossRef
  • Google Scholar

• 14

  KushwahaM. S.HaleviP.DobrzynskiL.Djafari-RouhaniB. (1993). Acoustic Band Structure of Periodic Elastic Composites. Phys. Rev. Lett.71 (13), 2022–2025. 10.1103/physrevlett.71.2022

  • CrossRef
  • Google Scholar

• 15

  LeeC.StampD.KapaniaN. R.Mur-MirandaOscarJ. (2016). “Harvesting Vibration Energy Using Nonlinear Oscillations of an Electromagnetic Inductor,” in Proceedings of SPIE-The International Society for Optical Engineering. 10.1117/12.849895

  • CrossRef
  • Google Scholar

• 16

  LiY.BakerE.ReissmanT.SunC.LiuW. K. (2017). Design of Mechanical Metamaterials for Simultaneous Vibration Isolation and Energy Harvesting. Appl. Phys. Lett.111 (25), 251903. 10.1063/1.5008674

  • CrossRef
  • Google Scholar

• 17

  LiB.ChengY.ZhuZ.ZhangF. "A Closed-form Solution for a Double Infinite Euler-Bernoulli Beam on a Viscoelastic Foundation Subjected to Harmonic Line Load." Earthq. Eng. Eng. Vib. Vol. 18No. 1 (2019): ppp. 129, 140. 10.1007/s11803-019-0494-9

  • CrossRef
  • Google Scholar

• 18

  LiuZ.ZhangX.MaoY.ZhuY. Y.YangZ.ChanC. T.et al (2000).Locally Resonant Sonic Materials. Science, 289. 1734–1736. 10.1126/science.289.5485.1734

  • CrossRef
  • Google Scholar

• 19

  LiuY.YuD.LiLi.ZhaoH.WenJ.WenX. (2007). Design Guidelines for Flexural Wave Attenuation of Slender Beams with Local Resonators. Phys. Lett. A362 (5), 344–347. 10.1016/j.physleta.2006.10.056

  • CrossRef
  • Google Scholar

• 20

  LiuZ.ChanC. T.ShengP. (2005). Analytic Model of Phononic Crystals with Local Resonances. Phys. Rev. B71 (No.1), 14103. 10.1103/physrevb.71.014103

  • CrossRef
  • Google Scholar

• 21

  MatlackK. H.BauhoferA.KrödelS.PalermoA.DaraioC. "Composite 3D-Printed Metastructures for Low-Frequency and Broadband Vibration Absorption." Proc. Natl. Acad. Sci. USA Vol. 113No. 30 (2016): ppp. 8386, 8390. 10.1073/pnas.1600171113

  • CrossRef
  • Google Scholar

• 22

  MuellerJ.KimS. E.SheaK.DaraioC. (2015). Tensile Properties of Inkjet 3D Printed Parts: Critical Process Parameters and Their Efficient Analysis. ASME Computers and Information in Engineering. 10.1115/DETC2015-48024

  • CrossRef
  • Google Scholar

• 23

  MuellerJ.SheaK. (2000). “The Effect of Build Orientation on the Mechanical Properties in Inkjet 3D-Printing,” in Proceedings of the 26th Annual International Solid Freeform Fabrication (SFF) Symposium, Austin, TX.

  • Google Scholar

• 24

  OniszczukZ. "Free Transverse Vibrations Of Elastically Connected Simply Supported Double-Beam Complex System." J. Sound Vibration Vol. 232No. 2 (2000): ppp. 387–403. 10.1006/jsvi.1999.2744

  • CrossRef
  • Google Scholar

• 25

  PilipovićA.RaosP.ŠercerM. (2009). Experimental Analysis of Properties of Materials for Rapid Prototyping. Int. J. Adv. Manufacturing Tech.40, 105.

  • Google Scholar

• 26

  ReichlK. K.AvciO.InmanD. J. "Temperature Dependent Damping in Additively Manufactured Polymer Structures." J. Appl. Comput. Mech. Vol. 7(2021): ppp. 993–1008.

  • Google Scholar

• 27

  ReichlK. K.InmanD. J. "Lumped Mass Model of a 1D Metastructure for Vibration Suppression with No Additional Mass." J. Sound Vibration Vol. 403 (2017): ppp. 75–89. 10.1016/j.jsv.2017.05.026

  • CrossRef
  • Google Scholar

• 28

  SharmaT.TyagiA. "Review of Mechanical Modelling of Fixed-Fixed Beams in RF MEMS Switches.", 2013Third International Conference on Advanced Computing and Communication Technologies (ACCT), ppp. 211–214. 10.1109/acct.2013.53

  • CrossRef
  • Google Scholar

• 29

  SinghR. (2011). Process Capability Study of Polyjet Printing for Plastic Components. J. Mech. Sci. Technol.25 (No. 4), 1011–1015. 10.1007/s12206-011-0203-8

  • CrossRef
  • Google Scholar

• 30

  SuginoC.XiaY.LeadenhamS.RuzzeneM.ErturkA. "A General Theory for Bandgap Estimation in Locally Resonant Metastructures." J. Sound Vibration Vol. 406 (2017): ppp. 104, 123. 10.1016/j.jsv.2017.06.004

  • CrossRef
  • Google Scholar

• 31

  VuH. V.OrdóñezA. M.KarnoppB. H. "Vibration Of A Double-Beam System." J. Sound Vibration Vol. 229No. 4 (2000): ppp. 807, 822. 10.1006/jsvi.1999.2528

  • CrossRef
  • Google Scholar

• 32

  WaichiW.AzidA.MajlisB. "Formulation of Stiffness Constant and Effective Mass for a Folded Beam." Arch. Mech. Vol. 62 (2010): ppp. 405–418.

  • Google Scholar

• 33

  WangM. Y.WangX. (2013). Frequency Band Structure of Locally Resonant Periodic Flexural Beams Suspended with Force-Moment Resonators. J. Phys. D: Appl. Phys.46 (No. 25), 255502. 10.1088/0022-3727/46/25/255502

  • CrossRef
  • Google Scholar

• 34

  WangZ.ZhangP.ZhangY. (2013). Locally Resonant Band Gaps in Flexural Vibrations of a Timoshenko Beam with Periodically Attached Multioscillators. Math. Probl. Eng.2013, 146975. 10.1155/2013/146975

  • CrossRef
  • Google Scholar

• 35

  WangZ.ZhangQ.ZhangK.HuG. "Tunable Digital Metamaterial for Broadband Vibration Isolation at Low Frequency." Adv. Mater. Vol. 28No. 44 (2016): ppp. 9857, 9861. 10.1002/adma.201604009

  • CrossRef
  • Google Scholar

• 36

  XiaY.RuzzeneM.ErturkA. "Bistable Attachments for Wideband Nonlinear Vibration Attenuation in a Metamaterial Beam." Nonlinear Dyn. Vol. 102No. 3 (2020): ppp. 1285, 1296. 10.1007/s11071-020-06008-4

  • CrossRef
  • Google Scholar

• 37

  XiaoY.WenJ.WenX. "Broadband Locally Resonant Beams Containing Multiple Periodic Arrays of Attached Resonators." Phys. Lett. A Vol. 376No. 16 (2012): ppp. 1384, 1390. 10.1016/j.physleta.2012.02.059

  • CrossRef
  • Google Scholar

• 38

  XiaoY.WenJ.YuD.WenX. "Flexural Wave Propagation in Beams with Periodically Attached Vibration Absorbers: Band-Gap Behavior and Band Formation Mechanisms." J. Sound Vibration Vol. 332No. 4 (2013): ppp. 867–893. 10.1016/j.jsv.2012.09.035

  • CrossRef
  • Google Scholar

• 39

  YuD.LiuY.WangG.ZhaoH.QiuJ. (2006). Flexural Vibration Band Gaps in Timoshenko Beams with Locally Resonant Structures. J. Appl. Phys.100 (No. 12), 124901. 10.1063/1.2400803

  • CrossRef
  • Google Scholar

• 40

  YuD.LiuY.ZhaoH.WangG.QiuJ. (2006). Flexural Vibration Band Gaps in Euler-Bernoulli Beams with Locally Resonant Structures with Two Degrees of Freedom. Phys. Rev. B73 (No. 6), 64301. 10.1103/physrevb.73.064301

  • CrossRef
  • Google Scholar

• 41

  ZhouJ.WangK.XuD.OuyangH. (2017). Local Resonator with High-Static-Low-Dynamic Stiffness for Lowering Band Gaps of Flexural Wave in Beams. J. Appl. Phys.121 (No. 4), 44902. 10.1063/1.4974299

  • CrossRef
  • Google Scholar