# A Novel Elastic Metamaterial with Multiple Resonators for Vibration Suppression

Authors: Saman Ahmadi Nooraldinvand; Hamid M. Sedighi; Amin Yaghootian
Journal: Advances in Condensed Matter Physics
 (2021)

---

## Abstract

In this paper, two models of elastic metamaterial containing one and two resonators are proposed to obtain the bandgaps with the aim of providing broadband vibration suppression. The model with one DOF is built by assembling several unite cells in which each unite cell consists of a rectangular frame as the base structure and a rack-and-pinion mechanism that is joined to the frame with a linear spring on both sides. In the second model with two DOF, a small mass is added while its center is attached to the center of the pinion on one side and the other side is connected to the rectangular frame via a linear spring. In the first mechanism, the pinion is considered as the single resonator, and in the 2DOF model, on the other hand, the pinion and small mass acted as multiple resonators. By obtaining the governing equations of motion for a single cell in each model, the dynamic behavior of two metastructures is thoroughly investigated. Therefore, the equations of motion for the two models are written in matrix form, and then, the dispersion relations are presented to analyze the influences of system parameters on the bandgaps’ starting/ending frequencies. Finally, two models are successfully compared and then numerically simulated via MATLAB-SIMULINK and MSC-ADAMS software. With the aid of closed-form expressions for starting/ending frequencies, the correlation between the system parameters and bandgap intervals can be readily recognized.


---

## Body

## 1. Introduction

Metamaterials are referred to the types of advanced materials that are synthetically made including small substructures which generally behave like a continuous material. The frequency bands at which acoustic and elastic waves cannot propagate are called bandgap, which is the most prominent feature of metamaterials. The propagation of waves of different wavelengths is controlled by the low/high-frequency bandgaps generated by the metamaterials. Given this unique property, metamaterials can be used in the field of acoustic insulation, filtration of waves, oscillations reduction, and sonic transmission [1–3]. The cell body of an acoustic metamaterial consists of a basic structure in which one or more locally vibrating components are located. Interactions between the base environment and local vibrators enable interesting physical phenomena such as bandwidth, negative effective density, negative modulus of elasticity, directional filtering and wave conduction, increased dissipation, and vibroacoustic attenuation [4, 5]. Basically, two types of local resonators are used: (1) translational resonator and (2) rotary resonator, which cause negative effective density and negative effective elastic modulus, respectively [6]. Huang et al. [7] demonstrated the consequence of using different equivalent models to represent a lattice system consisting of mass-in-mass units. Then they studied the dispersive wave propagation and compared to various equivalent models. They found that, if the classical elastic continuum is utilized to represent the original mass-in-mass structure, the effective mass density will be frequency dependent and may become negative for frequencies near the resonance frequency of the internal mass. In the meantime, if a multi-displacement microstructure continuum model is utilized, the dispersive behavior of wave propagation and the band gap structure can be adequately described.Active and passive periodic structures are made up of similar substructures or unit cells that are joined to each other in a similar way. The periodicity of these structures exhibits unique dynamic properties because they can prevent wave propagation and act as a mechanical filter [8]. Wave propagation is not possible in infinite periodic structures at particular frequency intervals. Photonic bandgaps are frequency bands in which acoustic/elastic waves cannot be transmitted. Unlike infinite periodic structures, wave transmission or vibration occurs sparingly if finite periodic structures are excited within the band gaps, which can be examined for frequency response function diagrams to obtain this wave or vibration transmission [9]. Propagation or passbands allow waves along these structures which propagate in certain frequency ranges without attenuation; and stopbands, attenuation bands, or bandgaps are the frequency bands in which waves can be attenuated [10]. In periodic structures, two common methods for generating phononic bandgaps are Bragg scattering and local resonances [9]. The band structure and wave reduction mentioned above are due to these two types of mechanisms that are clearly different from each other [11]. In the high-frequency mode, Bragg bandgaps can be used to filter the waves, but in the low-frequency mode, this is not possible [12]. On the other hand, at much lower frequencies than Bragg scattering, bandgaps can be created by local resonance [9]. Control of resonant-type bandgaps due to local resonance (LR) mechanism is easily possible, unlike Bragg-type bandgaps. In order to adjust the bandgap properties, the properties of the locally resonant structure can be changed. Chen and Chien [12] proposed a mass grid system, and before this new idea was applied to the field of elastic, it was presented in the field of acoustic.Lazarov and Jensen [10] conducted their study in both linear and nonlinear modes to propagate the wave in one-dimensional chains in which nonlinear oscillators are joined. They found that in the nonlinear state, the bandgap could shift depending on the amplitude and degree of local nonlinearities. Also, the results of their research showed that in the linear mode, the location of the gap is around the resonant frequency which allows the gap to be generated at a lower frequency domain. Casalotti et al. [13] examined an Euler-Bernoulli beam with an arrangement of nonlinear mass-spring substructures attached to investigate the ability of the nonlinear metamaterial beam to absorb multistate vibrations. Beam-connected substructures that act as local resonators or vibration absorbers created bandgaps. At first, they identified the frequency stop bands of the unite cell, and then in order to investigate the multifrequency stopband, the frequency response for the nonlinear metamaterial beam was extracted. They indicated that, for the metamaterial beam, the oscillations related to the lowest three vibration modes were significantly reduced by properly adjusting the constituent parameters of local absorbers. Zhou et al. [14] presented two models of an acoustic metamaterial beam in which local resonators were attached owning to both flexible bandgaps and high static stiffness. Initially, an equivalent model by mass-spring-beam of acoustic metamaterial beam with variable cross section was considered and analyzed analytically. The dispersion relationship of acoustic metamaterial beam was extracted, and the effect of different control parameters on bandgaps was investigated. Then, they simulated and examined a two-dimensional finite element model of acoustic metamaterial beam for the sake of validation in COMSOL Multiphysics. Their results showed that there is good compatibility between the analytical model and the two-dimensional model. Finally, a three-dimensional model of acoustic metamaterial beam was presented and discussed in two forms of equal and unequal thickness to measure the wave attenuation. Huang and Sun [15] studied a multiresonator mass-in-mass grid system and obtained its dispersion curves and bandgap. It was shown that bandgap can be displaced by changing the magnitude of the internal masses and the spring constant. They proposed a monoatomic model equivalent to the original system and found that the effective mass at bandgaps frequencies is negative. Finally, they introduced a microstructure continuous model that can obtain the dispersive behavior and bandgap structure of the original system. According to the study conducted by Sun et al. [16], it was shown that an acoustic absorber based on metamaterial can be thought of as a uniform isotropic beam with a number of very small mass-spring subsystems at different points. They also discussed how to create bandgap by subsystems, the negative effective mass, and stiffness. Their results showed that common mechanical vibration absorbers are the basis of the actual working mechanism of metamaterial beam. Shear forces and bending moments were then created to hold the beam straight and prevent wave propagation. An important and in-depth topical review in the field of active metamaterials and metadevices was performed by Xiao et al. [17]. They reported the advancement of active metadevices and metamaterials ranging from microwave to visible wavelengths, including milestones as well as the state of the art, and finally presented the future prospects together with several emerging tuning strategies and materials.Most studies of metamaterial models are based on a damper-spring-damper combination. However, this paper uses a new combination of rack-and-pinion mechanism and mass to achieve wider bandgaps and investigate the dynamic behavior of the new combination. In the current research, two models of elastic metamaterials are proposed to obtain broadband bandgaps for vibration suppression. The models are made from a combination of a rack-and-pinion mechanism connected to a concentrated mass. The governing equations of both models with one and two degrees of freedom are solved by MATLAB-SIMULINK and dispersion curves and frequency response diagrams are also drawn. Then, the effect of control parameters on the starting and ending frequencies of the bandgaps is mathematically investigated. Finally, both models are modeled and simulated in MSC-ADAMS software, and the desired graphs are numerically extracted. There is a satisfactory agreement between the simulated results and those obtained by numerical algorithms.


## 2. Models: 1DOF and 2DOF

In this section, two models consisting of one- and two-degree-of-freedom elastic metamaterials (EMM) are introduced and discussed. The models include rectangular cells containing a rack-and-pinion mechanism and a concentrated mass. The presence of a rack prevents slipping between the pinion and the cell which converts the transfer motion to a pure rolling motion. The aim of this study is to analytically obtain the bandgaps of the proposed models in order to reduce or eliminate the mechanical vibrations. For both one- and two-degree-of-freedom models, the mathematical models are presented and the staring/ending frequencies of bandgaps are mathematically extracted with close-form expressions.


### 2.1. One DOF Model of EMM

#### 2.1.1. One DOF Model: Mathematical Formulation

The proposed one-degree-of-freedom model of elastic metamaterials is shown in Figure1. As shown, the main body of this model consists of rectangular frames with a rack-and-pinion inside. Rectangular frames are then connected by linear springs. On the other hand, the pinion is connected to the walls of the rectangular frame by two linear springs.Figure 1
A schematic of the metamaterial arrangement of the 1DOF model.The parameters defined for the problem are shown in Figure2, which represents the unit cell of the one-degree-of-freedom model. Therefore, using Newton’s second law, the governing equations of a unit cell are derived as follows:(1)∑Fx=m1x¨n⇒−k1xn−xn−1−k1xn−xn+1−m2x¨n+rθ¨=m1x¨n,(2)∑MO=IOθ¨⇒−2k2rθr−m2x¨nr=IOθ¨,where m1 and m2 are frame and pinion masses, respectively, randθ are the radius and the rotational displacement of pinion, respectively, xn represents the translational displacement of the nth cell, and IO denotes the mass moment of inertia of the pinion with respect to the point o.Figure 2
A schematic of a unite cell of the 1DOF model.Now, Equations (1) and (2) are written in a matrix form as follows:(3)m1+m2m2rm2rIOx¨nθ¨+2k1002k2r2xnθ+−k1xn−1+xn+10=00,and the dimensionless parameters used in this model are defined by(4)M=m1m2,I⌢=IOm2r2,K⌢=k1k2,ω⌢n=k2m2,rx¯n=xn,τ=ωt,ω˜=ωω⌢n,where ω⌢n is the local natural frequency, M,I⌢,K⌢,x¯n,τ,andω˜ are the dimensionless mass, moment of inertia, stiffness, displacement, time, and frequency, respectively. By replacing these dimensionless parameters in Equation (3), the new form of governing equations can be written as follows:(5)M+1ω˜2ω˜2ω˜2I⌢ω˜2d2x¯ndτ2d2θdτ2+2K⌢002x¯nθ+−K⌢x¯n−1+x¯n+10=00.


#### 2.1.2. One DOF Model: Dispersion Properties

In this section, the solution corresponding to the harmonic wave for the 1DOF model is defined for extracting the dispersion relations and the bandgap frequencies:(6)x¯nθ=x¯n∘θ∘eiqa.n+τ+c.c,x¯n±1=x¯n±1∘eiqa.n+τ+c.c,where x¯n∘and θ∘ are the steady-state amplitude of the solution, qa is the phase factor, and nrepresents the periodic number. Moreover, a.n indicates the distance between nthcell from the origin (support) [18–22]. Substituting Equations (6) into (5), the new form of the governing equation in its matrix form is as follows:(7)M+1ω˜2ω˜2ω˜2I⌢ω˜2x¯n∘−1eiqa.n+τ+c.cθ∘−1eiqa.n+τ+c.c+2K⌢002x¯n∘eiqa.n+τ+c.cθ∘eiqa.n+τ+c.c+−K⌢x¯n∘eiqan−1eiτ+x¯n∘eiqan+1eiτ+c.c0=00.Now, by factoring and arranging the different expressions in Equation (7), the final form of the equations of motion for the considered system is given by(8)−M+1ω˜2+2K⌢−K⌢2cosqa−ω˜2−ω˜2−I⌢ω˜2+2x¯n∘θ∘eiqa.n+τ+c.c=00.The dispersion relation for the 1DOF model can be obtained from Equation (8). For a nontrivial solution, the determinant of the coefficient matrix must be equal to zero. Therefore, the dispersion relation is extracted as(9)I⌢1+M−1ω˜4+2−M−1−I⌢K⌢1+cosqaω˜2+4K⌢1−cosqa=0.By solving Equation (9) for ω˜, the relation for the acoustic and optical frequencies of the EMM can be achieved as described in Equations (10) and (11). Then, by employing the obtained closed-form expressions, the dispersion curves for the 1DOF model are shown in Figure 3.(10)ω˜Acoustic=Δ−ΦΨ,(11)ω˜Optical=Δ+ΦΨ,where(12)Δ=M+I⌢K⌢−I⌢K⌢cosqa+1,Ψ=I⌢+I⌢M−1,Φ=I⌢2K⌢2cos2qa+−2I⌢2K⌢2+4I⌢K⌢−4K⌢cosqa+1+4K⌢+2M−2I⌢K⌢−I⌢K⌢M.Figure 3
Dispersion curves for the 1DOF model.Equations (10) and (11) express the acoustic and optical frequency branches, respectively. Low-frequency and high-frequency dispersive curves are then presented by ω˜Acoustic and ω˜Optical, respectively [1]. In the above equations, at the point of qa=±π, the slope of the dispersion curve, namely, the velocity of the traveling wave, is equal to zero. This means the wave cannot propagate, similar to what happened at the bandgaps. This area is the first Brillouin zone’s boundary. Hence, substituting qa=π and qa=0 into Equations (10) and (11), and then solving these expressions, the starting and ending frequencies of the bandgaps for the 1DOF model are extracted [1, 20].


### 2.2. Two DOF Model of EMM

#### 2.2.1. Two DOF Model: Mathematical Formulation

In this section, a model of two degrees of freedom of EMM is presented (see Figure4). This model is an extended model of the first model discussed in the previous section. The unite cell of the model is made of rectangular frames with a rack-and-pinion inside it connected to a concentrated mass. As shown in Figure 5, the pinion and concentrated mass are joined to the rectangular frame by a linear spring with stiffness of k2. In addition, the rectangular frames are connected to each other by a linear spring with stiffness of k1 and create an integrated metastructure. Therefore, the governing equations of a unit cell can be described as follows:(13)∑Fx=m1x¨n⇒−k1xn−xn−1−k1xn−xn+1−m2x¨n+rθ¨−m3x¨n+x¨=m1x¨n,(14)∑MO=IOθ¨⇒−2k2rθr+k2xr−m2x¨nr=IOθ¨,(15)∑Fx=m3x¨⇒−m3x¨n−2k2x+k2rθ=m3x¨,where m3is the concentrated mass, and x denotes the translational displacement of concentrated mass.Figure 4
A schematic of the metamaterial arrangement of the 2DOF model.Figure 5
A schematic of a unite cell of 2DOF model.Now, Equations (13) to (15) are written in matrix form as follows:(16)m1+m2+m3m2rm3m2rIO0m30m3x¨nθ¨x¨+2k10002k2r2−k2r0−k2r2k2xnθx+−k1xn−1+xn+100=000.To express the problem in the dimensionless form, the following nondimensional parameters are defined for the 2DOF model:(17)m¯1=m2+m3m1,m¯2=m2m1,I¯=IOm2r2,K¯=k2k1,ωn1=k1m1ωn2=2k2m3,rx¯n=xn,rx¯=x,τ=ωt,ω¯1=ωωn1,ω¯2=ωωn2,where ωn1 and ωn2 stand for the natural local frequencies of the first and second modes, respectively, m¯1 and m¯2 are defined as the ratio of the total mass of the pinion and the concentrated mass to the mass of the rectangular frame and pinion-to-frame mass ratio, respectively, I¯,K¯,x¯n,x¯, and τ are the dimensionless moment of inertia, stiffness, displacement for cell, displacement for concentrated mass, and dimensionless time, and ω¯1andω¯2 denote the first and second dimensionless frequencies, respectively. Using the introduced dimensionless parameters, Equation (16) can be rewritten in a dimensionless form as follows:(18)1+m¯1ω¯12m¯2ω¯12m¯1−m¯2ω¯12ω¯12I¯ω¯120ω¯220ω¯22d2x¯ndτ2d2θdτ2d2x¯dτ2+20002K¯m¯2−K¯m¯20−121x¯nθx¯+−x¯n−1+x¯n+100=000.


#### 2.2.2. Two DOF Model: Dispersion Properties

In this section, the solutions corresponding to the harmonic wave for the 2DOF model are presented as follows for extracting the dispersion relations as well as the bandgap frequencies:(19)x¯nθx¯=x¯n∘θ∘x¯∘eiqa.n+τ+c.c,x¯n±1=x¯n±1∘eiqa.n+τ+c.c.Similar to the previous section,x¯n∘, θ∘, and x¯∘represent the steady-state amplitude of the solutions. Substituting Equations (19) into (18), the new forms of the governing equations using matrix notation are given byBy simplifying Equation (1), the final form of the governing equations yields(20)1+m¯1ω¯12m¯2ω¯12m¯1−m¯2ω¯12ω¯12I¯ω¯120ω¯220ω¯22x¯n∘−1eiqa.n+τ+c.cθ∘−1eiqa.n+τ+c.cx¯∘−1eiqa.n+τ+c.c+20002K¯m¯2−K¯m¯20−121x¯n∘eiqa.n+τ+c.cθ∘eiqa.n+τ+c.cx¯∘eiqa.n+τ+c.c+−x¯n∘eiqan−1eiτ+x¯n∘eiqan+1eiτ+c.c00=000.(21)−1+m¯1ω¯12+2−2cosqa−m¯2ω¯12−m¯1−m¯2ω¯12−ω¯12−I¯ω¯12+2K¯m¯2−K¯m¯2−ω¯22−121−ω¯22x¯n∘θ∘x¯∘eiqa.n+τ+c.c=000.Dispersion relation for the 2DOF model can be then obtained by manipulating Equation (21). For a nonzero solution, the determinant of the coefficient matrix should be set to zero. Therefore, the dispersion relation is extracted as follows:(22)Aω6+Bωn12+Cωn22ω4+Dωn12ωn22+Eωn14ω2+F=0,where(23)A=2m¯2m¯2I¯−1−I¯,B=−2K¯2+m¯2+4m¯2I¯1−cosqa,C=m¯2m¯2+2I¯+m¯1m¯21+2I¯,D=3K¯1+m¯1−4m¯2I¯1−cosqa,E=8K¯1−cosqa,F=−6K¯1−cosqaωn14ωn22.Similar to the explanations given for the one-degree-of-freedom model, Equation (22) shows the relationship between ω and qa. By plotting this equation, the dispersion curves are then available. By replacing qa=π in Equation (22), the expressions for the starting frequencies can be achieved as follows:(24)Fs1=H3A−12CB+B,Fs2=H3A+CB+B.Now, by substituting qa = 0 in Equation (22), the correlations between the system parameters and ending frequencies are given by(25)Fe1=K−L2A,Fe2=K+L2A,where(26)A=2m¯22I¯+2m¯22+2m¯2I¯,B=E+E+D−F2−C31/2+D−F1/3,C=13AH23A+G,D=GH6A2,E=H327A3,F=6Kωn14ωn22A,G=16Kωn14+3K+3Km¯1−8m¯2I¯ωn12ωn22,H=m¯2ωn22m¯2+m¯1+2I¯+2m¯1I¯+ωn12−4K¯−2K¯m¯2+8m¯2I¯,I=m¯12m¯222I¯+12+2m¯1m¯23+4m¯1m¯22I¯2I¯+1+2I¯m¯2+m¯222,J=4K¯m¯1m¯224I¯+5+8K¯m¯1m¯2I¯−1+16K¯m¯22I¯+1+8K¯m¯2I¯−4K¯m¯23,K=H−8m¯2I¯ωn12,L=16K¯m¯224+m¯2+1ωn14+Iωn24+Jωn12ωn22.


## 2.1. One DOF Model of EMM

### 2.1.1. One DOF Model: Mathematical Formulation

The proposed one-degree-of-freedom model of elastic metamaterials is shown in Figure1. As shown, the main body of this model consists of rectangular frames with a rack-and-pinion inside. Rectangular frames are then connected by linear springs. On the other hand, the pinion is connected to the walls of the rectangular frame by two linear springs.Figure 1
A schematic of the metamaterial arrangement of the 1DOF model.The parameters defined for the problem are shown in Figure2, which represents the unit cell of the one-degree-of-freedom model. Therefore, using Newton’s second law, the governing equations of a unit cell are derived as follows:(1)∑Fx=m1x¨n⇒−k1xn−xn−1−k1xn−xn+1−m2x¨n+rθ¨=m1x¨n,(2)∑MO=IOθ¨⇒−2k2rθr−m2x¨nr=IOθ¨,where m1 and m2 are frame and pinion masses, respectively, randθ are the radius and the rotational displacement of pinion, respectively, xn represents the translational displacement of the nth cell, and IO denotes the mass moment of inertia of the pinion with respect to the point o.Figure 2
A schematic of a unite cell of the 1DOF model.Now, Equations (1) and (2) are written in a matrix form as follows:(3)m1+m2m2rm2rIOx¨nθ¨+2k1002k2r2xnθ+−k1xn−1+xn+10=00,and the dimensionless parameters used in this model are defined by(4)M=m1m2,I⌢=IOm2r2,K⌢=k1k2,ω⌢n=k2m2,rx¯n=xn,τ=ωt,ω˜=ωω⌢n,where ω⌢n is the local natural frequency, M,I⌢,K⌢,x¯n,τ,andω˜ are the dimensionless mass, moment of inertia, stiffness, displacement, time, and frequency, respectively. By replacing these dimensionless parameters in Equation (3), the new form of governing equations can be written as follows:(5)M+1ω˜2ω˜2ω˜2I⌢ω˜2d2x¯ndτ2d2θdτ2+2K⌢002x¯nθ+−K⌢x¯n−1+x¯n+10=00.


### 2.1.2. One DOF Model: Dispersion Properties

In this section, the solution corresponding to the harmonic wave for the 1DOF model is defined for extracting the dispersion relations and the bandgap frequencies:(6)x¯nθ=x¯n∘θ∘eiqa.n+τ+c.c,x¯n±1=x¯n±1∘eiqa.n+τ+c.c,where x¯n∘and θ∘ are the steady-state amplitude of the solution, qa is the phase factor, and nrepresents the periodic number. Moreover, a.n indicates the distance between nthcell from the origin (support) [18–22]. Substituting Equations (6) into (5), the new form of the governing equation in its matrix form is as follows:(7)M+1ω˜2ω˜2ω˜2I⌢ω˜2x¯n∘−1eiqa.n+τ+c.cθ∘−1eiqa.n+τ+c.c+2K⌢002x¯n∘eiqa.n+τ+c.cθ∘eiqa.n+τ+c.c+−K⌢x¯n∘eiqan−1eiτ+x¯n∘eiqan+1eiτ+c.c0=00.Now, by factoring and arranging the different expressions in Equation (7), the final form of the equations of motion for the considered system is given by(8)−M+1ω˜2+2K⌢−K⌢2cosqa−ω˜2−ω˜2−I⌢ω˜2+2x¯n∘θ∘eiqa.n+τ+c.c=00.The dispersion relation for the 1DOF model can be obtained from Equation (8). For a nontrivial solution, the determinant of the coefficient matrix must be equal to zero. Therefore, the dispersion relation is extracted as(9)I⌢1+M−1ω˜4+2−M−1−I⌢K⌢1+cosqaω˜2+4K⌢1−cosqa=0.By solving Equation (9) for ω˜, the relation for the acoustic and optical frequencies of the EMM can be achieved as described in Equations (10) and (11). Then, by employing the obtained closed-form expressions, the dispersion curves for the 1DOF model are shown in Figure 3.(10)ω˜Acoustic=Δ−ΦΨ,(11)ω˜Optical=Δ+ΦΨ,where(12)Δ=M+I⌢K⌢−I⌢K⌢cosqa+1,Ψ=I⌢+I⌢M−1,Φ=I⌢2K⌢2cos2qa+−2I⌢2K⌢2+4I⌢K⌢−4K⌢cosqa+1+4K⌢+2M−2I⌢K⌢−I⌢K⌢M.Figure 3
Dispersion curves for the 1DOF model.Equations (10) and (11) express the acoustic and optical frequency branches, respectively. Low-frequency and high-frequency dispersive curves are then presented by ω˜Acoustic and ω˜Optical, respectively [1]. In the above equations, at the point of qa=±π, the slope of the dispersion curve, namely, the velocity of the traveling wave, is equal to zero. This means the wave cannot propagate, similar to what happened at the bandgaps. This area is the first Brillouin zone’s boundary. Hence, substituting qa=π and qa=0 into Equations (10) and (11), and then solving these expressions, the starting and ending frequencies of the bandgaps for the 1DOF model are extracted [1, 20].


## 2.1.1. One DOF Model: Mathematical Formulation

The proposed one-degree-of-freedom model of elastic metamaterials is shown in Figure1. As shown, the main body of this model consists of rectangular frames with a rack-and-pinion inside. Rectangular frames are then connected by linear springs. On the other hand, the pinion is connected to the walls of the rectangular frame by two linear springs.Figure 1
A schematic of the metamaterial arrangement of the 1DOF model.The parameters defined for the problem are shown in Figure2, which represents the unit cell of the one-degree-of-freedom model. Therefore, using Newton’s second law, the governing equations of a unit cell are derived as follows:(1)∑Fx=m1x¨n⇒−k1xn−xn−1−k1xn−xn+1−m2x¨n+rθ¨=m1x¨n,(2)∑MO=IOθ¨⇒−2k2rθr−m2x¨nr=IOθ¨,where m1 and m2 are frame and pinion masses, respectively, randθ are the radius and the rotational displacement of pinion, respectively, xn represents the translational displacement of the nth cell, and IO denotes the mass moment of inertia of the pinion with respect to the point o.Figure 2
A schematic of a unite cell of the 1DOF model.Now, Equations (1) and (2) are written in a matrix form as follows:(3)m1+m2m2rm2rIOx¨nθ¨+2k1002k2r2xnθ+−k1xn−1+xn+10=00,and the dimensionless parameters used in this model are defined by(4)M=m1m2,I⌢=IOm2r2,K⌢=k1k2,ω⌢n=k2m2,rx¯n=xn,τ=ωt,ω˜=ωω⌢n,where ω⌢n is the local natural frequency, M,I⌢,K⌢,x¯n,τ,andω˜ are the dimensionless mass, moment of inertia, stiffness, displacement, time, and frequency, respectively. By replacing these dimensionless parameters in Equation (3), the new form of governing equations can be written as follows:(5)M+1ω˜2ω˜2ω˜2I⌢ω˜2d2x¯ndτ2d2θdτ2+2K⌢002x¯nθ+−K⌢x¯n−1+x¯n+10=00.


## 2.1.2. One DOF Model: Dispersion Properties

In this section, the solution corresponding to the harmonic wave for the 1DOF model is defined for extracting the dispersion relations and the bandgap frequencies:(6)x¯nθ=x¯n∘θ∘eiqa.n+τ+c.c,x¯n±1=x¯n±1∘eiqa.n+τ+c.c,where x¯n∘and θ∘ are the steady-state amplitude of the solution, qa is the phase factor, and nrepresents the periodic number. Moreover, a.n indicates the distance between nthcell from the origin (support) [18–22]. Substituting Equations (6) into (5), the new form of the governing equation in its matrix form is as follows:(7)M+1ω˜2ω˜2ω˜2I⌢ω˜2x¯n∘−1eiqa.n+τ+c.cθ∘−1eiqa.n+τ+c.c+2K⌢002x¯n∘eiqa.n+τ+c.cθ∘eiqa.n+τ+c.c+−K⌢x¯n∘eiqan−1eiτ+x¯n∘eiqan+1eiτ+c.c0=00.Now, by factoring and arranging the different expressions in Equation (7), the final form of the equations of motion for the considered system is given by(8)−M+1ω˜2+2K⌢−K⌢2cosqa−ω˜2−ω˜2−I⌢ω˜2+2x¯n∘θ∘eiqa.n+τ+c.c=00.The dispersion relation for the 1DOF model can be obtained from Equation (8). For a nontrivial solution, the determinant of the coefficient matrix must be equal to zero. Therefore, the dispersion relation is extracted as(9)I⌢1+M−1ω˜4+2−M−1−I⌢K⌢1+cosqaω˜2+4K⌢1−cosqa=0.By solving Equation (9) for ω˜, the relation for the acoustic and optical frequencies of the EMM can be achieved as described in Equations (10) and (11). Then, by employing the obtained closed-form expressions, the dispersion curves for the 1DOF model are shown in Figure 3.(10)ω˜Acoustic=Δ−ΦΨ,(11)ω˜Optical=Δ+ΦΨ,where(12)Δ=M+I⌢K⌢−I⌢K⌢cosqa+1,Ψ=I⌢+I⌢M−1,Φ=I⌢2K⌢2cos2qa+−2I⌢2K⌢2+4I⌢K⌢−4K⌢cosqa+1+4K⌢+2M−2I⌢K⌢−I⌢K⌢M.Figure 3
Dispersion curves for the 1DOF model.Equations (10) and (11) express the acoustic and optical frequency branches, respectively. Low-frequency and high-frequency dispersive curves are then presented by ω˜Acoustic and ω˜Optical, respectively [1]. In the above equations, at the point of qa=±π, the slope of the dispersion curve, namely, the velocity of the traveling wave, is equal to zero. This means the wave cannot propagate, similar to what happened at the bandgaps. This area is the first Brillouin zone’s boundary. Hence, substituting qa=π and qa=0 into Equations (10) and (11), and then solving these expressions, the starting and ending frequencies of the bandgaps for the 1DOF model are extracted [1, 20].


## 2.2. Two DOF Model of EMM

### 2.2.1. Two DOF Model: Mathematical Formulation

In this section, a model of two degrees of freedom of EMM is presented (see Figure4). This model is an extended model of the first model discussed in the previous section. The unite cell of the model is made of rectangular frames with a rack-and-pinion inside it connected to a concentrated mass. As shown in Figure 5, the pinion and concentrated mass are joined to the rectangular frame by a linear spring with stiffness of k2. In addition, the rectangular frames are connected to each other by a linear spring with stiffness of k1 and create an integrated metastructure. Therefore, the governing equations of a unit cell can be described as follows:(13)∑Fx=m1x¨n⇒−k1xn−xn−1−k1xn−xn+1−m2x¨n+rθ¨−m3x¨n+x¨=m1x¨n,(14)∑MO=IOθ¨⇒−2k2rθr+k2xr−m2x¨nr=IOθ¨,(15)∑Fx=m3x¨⇒−m3x¨n−2k2x+k2rθ=m3x¨,where m3is the concentrated mass, and x denotes the translational displacement of concentrated mass.Figure 4
A schematic of the metamaterial arrangement of the 2DOF model.Figure 5
A schematic of a unite cell of 2DOF model.Now, Equations (13) to (15) are written in matrix form as follows:(16)m1+m2+m3m2rm3m2rIO0m30m3x¨nθ¨x¨+2k10002k2r2−k2r0−k2r2k2xnθx+−k1xn−1+xn+100=000.To express the problem in the dimensionless form, the following nondimensional parameters are defined for the 2DOF model:(17)m¯1=m2+m3m1,m¯2=m2m1,I¯=IOm2r2,K¯=k2k1,ωn1=k1m1ωn2=2k2m3,rx¯n=xn,rx¯=x,τ=ωt,ω¯1=ωωn1,ω¯2=ωωn2,where ωn1 and ωn2 stand for the natural local frequencies of the first and second modes, respectively, m¯1 and m¯2 are defined as the ratio of the total mass of the pinion and the concentrated mass to the mass of the rectangular frame and pinion-to-frame mass ratio, respectively, I¯,K¯,x¯n,x¯, and τ are the dimensionless moment of inertia, stiffness, displacement for cell, displacement for concentrated mass, and dimensionless time, and ω¯1andω¯2 denote the first and second dimensionless frequencies, respectively. Using the introduced dimensionless parameters, Equation (16) can be rewritten in a dimensionless form as follows:(18)1+m¯1ω¯12m¯2ω¯12m¯1−m¯2ω¯12ω¯12I¯ω¯120ω¯220ω¯22d2x¯ndτ2d2θdτ2d2x¯dτ2+20002K¯m¯2−K¯m¯20−121x¯nθx¯+−x¯n−1+x¯n+100=000.


### 2.2.2. Two DOF Model: Dispersion Properties

In this section, the solutions corresponding to the harmonic wave for the 2DOF model are presented as follows for extracting the dispersion relations as well as the bandgap frequencies:(19)x¯nθx¯=x¯n∘θ∘x¯∘eiqa.n+τ+c.c,x¯n±1=x¯n±1∘eiqa.n+τ+c.c.Similar to the previous section,x¯n∘, θ∘, and x¯∘represent the steady-state amplitude of the solutions. Substituting Equations (19) into (18), the new forms of the governing equations using matrix notation are given byBy simplifying Equation (1), the final form of the governing equations yields(20)1+m¯1ω¯12m¯2ω¯12m¯1−m¯2ω¯12ω¯12I¯ω¯120ω¯220ω¯22x¯n∘−1eiqa.n+τ+c.cθ∘−1eiqa.n+τ+c.cx¯∘−1eiqa.n+τ+c.c+20002K¯m¯2−K¯m¯20−121x¯n∘eiqa.n+τ+c.cθ∘eiqa.n+τ+c.cx¯∘eiqa.n+τ+c.c+−x¯n∘eiqan−1eiτ+x¯n∘eiqan+1eiτ+c.c00=000.(21)−1+m¯1ω¯12+2−2cosqa−m¯2ω¯12−m¯1−m¯2ω¯12−ω¯12−I¯ω¯12+2K¯m¯2−K¯m¯2−ω¯22−121−ω¯22x¯n∘θ∘x¯∘eiqa.n+τ+c.c=000.Dispersion relation for the 2DOF model can be then obtained by manipulating Equation (21). For a nonzero solution, the determinant of the coefficient matrix should be set to zero. Therefore, the dispersion relation is extracted as follows:(22)Aω6+Bωn12+Cωn22ω4+Dωn12ωn22+Eωn14ω2+F=0,where(23)A=2m¯2m¯2I¯−1−I¯,B=−2K¯2+m¯2+4m¯2I¯1−cosqa,C=m¯2m¯2+2I¯+m¯1m¯21+2I¯,D=3K¯1+m¯1−4m¯2I¯1−cosqa,E=8K¯1−cosqa,F=−6K¯1−cosqaωn14ωn22.Similar to the explanations given for the one-degree-of-freedom model, Equation (22) shows the relationship between ω and qa. By plotting this equation, the dispersion curves are then available. By replacing qa=π in Equation (22), the expressions for the starting frequencies can be achieved as follows:(24)Fs1=H3A−12CB+B,Fs2=H3A+CB+B.Now, by substituting qa = 0 in Equation (22), the correlations between the system parameters and ending frequencies are given by(25)Fe1=K−L2A,Fe2=K+L2A,where(26)A=2m¯22I¯+2m¯22+2m¯2I¯,B=E+E+D−F2−C31/2+D−F1/3,C=13AH23A+G,D=GH6A2,E=H327A3,F=6Kωn14ωn22A,G=16Kωn14+3K+3Km¯1−8m¯2I¯ωn12ωn22,H=m¯2ωn22m¯2+m¯1+2I¯+2m¯1I¯+ωn12−4K¯−2K¯m¯2+8m¯2I¯,I=m¯12m¯222I¯+12+2m¯1m¯23+4m¯1m¯22I¯2I¯+1+2I¯m¯2+m¯222,J=4K¯m¯1m¯224I¯+5+8K¯m¯1m¯2I¯−1+16K¯m¯22I¯+1+8K¯m¯2I¯−4K¯m¯23,K=H−8m¯2I¯ωn12,L=16K¯m¯224+m¯2+1ωn14+Iωn24+Jωn12ωn22.


## 2.2.1. Two DOF Model: Mathematical Formulation

In this section, a model of two degrees of freedom of EMM is presented (see Figure4). This model is an extended model of the first model discussed in the previous section. The unite cell of the model is made of rectangular frames with a rack-and-pinion inside it connected to a concentrated mass. As shown in Figure 5, the pinion and concentrated mass are joined to the rectangular frame by a linear spring with stiffness of k2. In addition, the rectangular frames are connected to each other by a linear spring with stiffness of k1 and create an integrated metastructure. Therefore, the governing equations of a unit cell can be described as follows:(13)∑Fx=m1x¨n⇒−k1xn−xn−1−k1xn−xn+1−m2x¨n+rθ¨−m3x¨n+x¨=m1x¨n,(14)∑MO=IOθ¨⇒−2k2rθr+k2xr−m2x¨nr=IOθ¨,(15)∑Fx=m3x¨⇒−m3x¨n−2k2x+k2rθ=m3x¨,where m3is the concentrated mass, and x denotes the translational displacement of concentrated mass.Figure 4
A schematic of the metamaterial arrangement of the 2DOF model.Figure 5
A schematic of a unite cell of 2DOF model.Now, Equations (13) to (15) are written in matrix form as follows:(16)m1+m2+m3m2rm3m2rIO0m30m3x¨nθ¨x¨+2k10002k2r2−k2r0−k2r2k2xnθx+−k1xn−1+xn+100=000.To express the problem in the dimensionless form, the following nondimensional parameters are defined for the 2DOF model:(17)m¯1=m2+m3m1,m¯2=m2m1,I¯=IOm2r2,K¯=k2k1,ωn1=k1m1ωn2=2k2m3,rx¯n=xn,rx¯=x,τ=ωt,ω¯1=ωωn1,ω¯2=ωωn2,where ωn1 and ωn2 stand for the natural local frequencies of the first and second modes, respectively, m¯1 and m¯2 are defined as the ratio of the total mass of the pinion and the concentrated mass to the mass of the rectangular frame and pinion-to-frame mass ratio, respectively, I¯,K¯,x¯n,x¯, and τ are the dimensionless moment of inertia, stiffness, displacement for cell, displacement for concentrated mass, and dimensionless time, and ω¯1andω¯2 denote the first and second dimensionless frequencies, respectively. Using the introduced dimensionless parameters, Equation (16) can be rewritten in a dimensionless form as follows:(18)1+m¯1ω¯12m¯2ω¯12m¯1−m¯2ω¯12ω¯12I¯ω¯120ω¯220ω¯22d2x¯ndτ2d2θdτ2d2x¯dτ2+20002K¯m¯2−K¯m¯20−121x¯nθx¯+−x¯n−1+x¯n+100=000.


## 2.2.2. Two DOF Model: Dispersion Properties

In this section, the solutions corresponding to the harmonic wave for the 2DOF model are presented as follows for extracting the dispersion relations as well as the bandgap frequencies:(19)x¯nθx¯=x¯n∘θ∘x¯∘eiqa.n+τ+c.c,x¯n±1=x¯n±1∘eiqa.n+τ+c.c.Similar to the previous section,x¯n∘, θ∘, and x¯∘represent the steady-state amplitude of the solutions. Substituting Equations (19) into (18), the new forms of the governing equations using matrix notation are given byBy simplifying Equation (1), the final form of the governing equations yields(20)1+m¯1ω¯12m¯2ω¯12m¯1−m¯2ω¯12ω¯12I¯ω¯120ω¯220ω¯22x¯n∘−1eiqa.n+τ+c.cθ∘−1eiqa.n+τ+c.cx¯∘−1eiqa.n+τ+c.c+20002K¯m¯2−K¯m¯20−121x¯n∘eiqa.n+τ+c.cθ∘eiqa.n+τ+c.cx¯∘eiqa.n+τ+c.c+−x¯n∘eiqan−1eiτ+x¯n∘eiqan+1eiτ+c.c00=000.(21)−1+m¯1ω¯12+2−2cosqa−m¯2ω¯12−m¯1−m¯2ω¯12−ω¯12−I¯ω¯12+2K¯m¯2−K¯m¯2−ω¯22−121−ω¯22x¯n∘θ∘x¯∘eiqa.n+τ+c.c=000.Dispersion relation for the 2DOF model can be then obtained by manipulating Equation (21). For a nonzero solution, the determinant of the coefficient matrix should be set to zero. Therefore, the dispersion relation is extracted as follows:(22)Aω6+Bωn12+Cωn22ω4+Dωn12ωn22+Eωn14ω2+F=0,where(23)A=2m¯2m¯2I¯−1−I¯,B=−2K¯2+m¯2+4m¯2I¯1−cosqa,C=m¯2m¯2+2I¯+m¯1m¯21+2I¯,D=3K¯1+m¯1−4m¯2I¯1−cosqa,E=8K¯1−cosqa,F=−6K¯1−cosqaωn14ωn22.Similar to the explanations given for the one-degree-of-freedom model, Equation (22) shows the relationship between ω and qa. By plotting this equation, the dispersion curves are then available. By replacing qa=π in Equation (22), the expressions for the starting frequencies can be achieved as follows:(24)Fs1=H3A−12CB+B,Fs2=H3A+CB+B.Now, by substituting qa = 0 in Equation (22), the correlations between the system parameters and ending frequencies are given by(25)Fe1=K−L2A,Fe2=K+L2A,where(26)A=2m¯22I¯+2m¯22+2m¯2I¯,B=E+E+D−F2−C31/2+D−F1/3,C=13AH23A+G,D=GH6A2,E=H327A3,F=6Kωn14ωn22A,G=16Kωn14+3K+3Km¯1−8m¯2I¯ωn12ωn22,H=m¯2ωn22m¯2+m¯1+2I¯+2m¯1I¯+ωn12−4K¯−2K¯m¯2+8m¯2I¯,I=m¯12m¯222I¯+12+2m¯1m¯23+4m¯1m¯22I¯2I¯+1+2I¯m¯2+m¯222,J=4K¯m¯1m¯224I¯+5+8K¯m¯1m¯2I¯−1+16K¯m¯22I¯+1+8K¯m¯2I¯−4K¯m¯23,K=H−8m¯2I¯ωn12,L=16K¯m¯224+m¯2+1ωn14+Iωn24+Jωn12ωn22.


## 3. The Integrity of One and Two DOF Models

To verify the results of this study extracted by mathematical modeling, both elastic metamaterials are modeled and simulated in MSC-ADAMS software under the assumption of rigid dynamics, as shown in Figures6–9. Furthermore, the frequency response analysis (FRA) is performed for both models and frequency diagrams are then plotted. In the following section, the transient specifications of the system are expressed by parameter Π and the ratio of the fifth cell displacement to the displacement of the first cell and then compared with numerical solutions.Figure 6
A schematic of the EMM arrangement of 1DOF model-MSC-ADAMS software.Figure 7
A schematic of a unite cell of 1DOF model-MSC-ADAMS software.Figure 8
A schematic of the EMM arrangement of 2DOF model-MSC-ADAMS software.Figure 9
A schematic of a unite cell of 2DOF model-MSC-ADAMS software.For numerical simulations, MATLAB/Simulink software is utilized, and the governing equations of proposed models are solved numerically. The block diagrams drawn by MATLAB/Simulink models are displayed in Figure10. It is worth mentioning that each section in the diagram contains the equations governing the dynamics of a single cell.Figure 10
Block diagrams plotted in MATLAB/Simulink for the numerical solution of governing equations. (a) 1DOF model. (b) 2DOF model, array of cells (top) and a unit cell (bottom).
(a)(b)For both models, the frequency response diagrams obtained from numerical solutions and simulations are plotted in Figures11 and 12. As indicated, both analyses provide the same intervals for system bandgaps, demonstrating the satisfactory agreement between the results of both procedures which indicates the integrity of our mathematical modeling as well as numerical calculations.Figure 11
FRA for 1DOF model, comparison of numerical solution (MATLAB/Simulink) with simulation (MSC-ADAMS).Figure 12
FRA for 2DOF model, comparison of numerical solution (MATLAB/Simulink) with simulation (MSC-ADAMS).


## 4. Results and Discussion

After verifying the soundness of the numerical solutions, the dispersion characteristics and the effect of different system parameters on the bandgap intervals are comprehensively studied. In this section, it is aimed to examine the bandgaps of the introduced models and demonstrate in which conditions they can be extended or transmitted.


### 4.1. One DOF Model of EMM

#### 4.1.1. Dispersion Curve for the 1DOF Model

Figure2 shows the dispersion curves of the one-degree-of-freedom model, which shows the acoustic and optical wave modes as well as the bandgap interval. The values of the selected parameters are M=1.56,I⌢=1.5,K⌢=0.35,andω⌢n=8.9. Considering these parameters, the bandgap of the one-degree-of-freedom model is placed between the frequencies of 6.15 to 11.94 (Hz), which is the same as the frequency range obtained from Figure 11.


#### 4.1.2. The Effect of System Parameters in the 1DOF Model

The surfaces drawn in Figure13 exhibit the effect of different dimensionless parameters on the starting and ending frequencies of their corresponding bandgaps. The top and bottom surfaces indicate the starting and ending frequencies, respectively, and the distance between them indicates the bandgap interval of the 1DOF model. As can be seen in Figure 13, there is a direct relationship between the bandgap and the dimensionless stiffness, meaning that as the dimensionless stiffness increases, the bandgap is also enhanced. On the other hand, by increasing the dimensionless mass, the bandgap decreases. It is deduced that the stiffness parameter has a positive effect on the expansion of the bandgap and the effect of the mass parameter is to weaken the bandgap interval. Moreover, according to the illustrated results, in general, the dimensionless moment of inertia directly affects the bandgap range, which means that as the dimensionless moment of inertia increases, the bandgap of the model is also enlarged. Although the bandgap interval slightly decreases at the lower values, when the dimensionless moment of inertia takes the higher values, this interval is then enhanced satisfactorily.Figure 13
The effect of changing (a)K⌢ and M and (b) K⌢ and I⌢on the starting and ending frequencies of the bandgap for 1DOF model.
(a)(b)


### 4.2. 2DOF Model of EMM

#### 4.2.1. Dispersion Curve for the 2DOF Model

In order to display the dispersion curves of the modified model with multiple resonators, one can solve Equation (22) for ω and extract the desired curves. The values of the selected parameters are m¯1=1.04,m¯2=0.64,I¯=1.5,andK¯=2.86,ωn1=4.2rad/s,ωn2=15.9rad/s. Considering the above parameters, the bandgap starting frequencies (Fs1and Fs2) are calculated by substituting qa=±π in the acoustic and optical wave modes, respectively. The resulting plots are shown in Figure 14 by Ds(1) and Ds(2), respectively. In addition, the bandgap ending frequencies (Fe1and Fe2) are also extracted by replacing qa=0 in two optical roots of Equation (22), respectively. The plotted results are illustrated in Figure 14 by Ds(2) and Ds(3), respectively. As indicated in Figure 14, the difference between the first starting frequency and the first ending frequency demonstrates the first bandgap, and the difference between the second starting frequency and the second ending frequency represents the second bandgap. According to the illustrated results in Figure 12, one can find that there is a small area between the first and second bandgaps. This consequence is also demonstrated in Figure 14 with a slight difference between the highest and lowest points of Ds(2). The slight difference leads to the flattening of the Ds(2) curve.Figure 14
Dispersion curves for the 2DOF model.


#### 4.2.2. Comparison of Dispersion Curves for Two Models

As shown in Figure15, by adding the concentrated mass and modifying the 1DOF model to a 2DOF model, two bandgaps are created in the dynamics of the new model, which allows us to reduce or suppress the vibrations of the modified system in more frequency ranges. From the comparison between the dispersion curves of both models, one can conclude that the first bandgap in the two-degree-of-freedom model starts and ends at lower frequencies than that of the one-degree-of-freedom model.Figure 15
Comparison of dispersion curves in two models.


#### 4.2.3. The Effect of System Parameters in the 2DOF Model

Figure16 displays well that when the parameter m¯1 shifts upward, the first starting frequency also increases and the first ending frequency reduces, which satisfactorily expands the first bandgap of the modified model. On the other hand, by increasing the first bandgap, the second one decreases and the starting and ending frequencies eventually reach a constant value.Figure 16
The effect ofm¯1 parameter on bandgap for the 2DOF model.With increasing the parameterK¯, the first bandgap can be extended to cover a larger frequency range, as shown in Figure 17, although the first starting frequency will eventually approach a constant value. However, the change in the second bandgap of the model does not show the same trend, eventually converges to an almost constant value, and no longer changes significantly.Figure 17
The effect ofK¯ parameter on bandgap for the 2DOF model.For the sake of extending the first and second bandgaps of the modified model, the dimensionless moment of inertiaI¯ has to be increased. As indicated in Figure 18, by continuing the process of increasing the dimensionless moment of inertia, the first ending and the second starting frequencies will not change and take a constant value.Figure 18
The effect ofI¯ parameter on bandgap for the 2DOF model.As the first and second local frequencies increase, the same behavior cannot be observed in the first bandgap of this model, although it decreases with small values and follows an increasing trend at higher values (see Figures19 and 20). In general, it can be inferred that the first bandgap is expanded. Finally, it is found that with the increase of the first local frequency, the expansion of the second bandgap is not uniform, although, in general, it becomes larger. Moreover, it should be pointed out that when the second local frequency increases, the second bandgap of the modified model is reduced.Figure 19
The effect ofωn1 parameter on bandgap for the 2DOF model.Figure 20
The effect ofωn2 parameter on bandgap for the 2DOF model.


## 4.1. One DOF Model of EMM

### 4.1.1. Dispersion Curve for the 1DOF Model

Figure2 shows the dispersion curves of the one-degree-of-freedom model, which shows the acoustic and optical wave modes as well as the bandgap interval. The values of the selected parameters are M=1.56,I⌢=1.5,K⌢=0.35,andω⌢n=8.9. Considering these parameters, the bandgap of the one-degree-of-freedom model is placed between the frequencies of 6.15 to 11.94 (Hz), which is the same as the frequency range obtained from Figure 11.


### 4.1.2. The Effect of System Parameters in the 1DOF Model

The surfaces drawn in Figure13 exhibit the effect of different dimensionless parameters on the starting and ending frequencies of their corresponding bandgaps. The top and bottom surfaces indicate the starting and ending frequencies, respectively, and the distance between them indicates the bandgap interval of the 1DOF model. As can be seen in Figure 13, there is a direct relationship between the bandgap and the dimensionless stiffness, meaning that as the dimensionless stiffness increases, the bandgap is also enhanced. On the other hand, by increasing the dimensionless mass, the bandgap decreases. It is deduced that the stiffness parameter has a positive effect on the expansion of the bandgap and the effect of the mass parameter is to weaken the bandgap interval. Moreover, according to the illustrated results, in general, the dimensionless moment of inertia directly affects the bandgap range, which means that as the dimensionless moment of inertia increases, the bandgap of the model is also enlarged. Although the bandgap interval slightly decreases at the lower values, when the dimensionless moment of inertia takes the higher values, this interval is then enhanced satisfactorily.Figure 13
The effect of changing (a)K⌢ and M and (b) K⌢ and I⌢on the starting and ending frequencies of the bandgap for 1DOF model.
(a)(b)


## 4.1.1. Dispersion Curve for the 1DOF Model

Figure2 shows the dispersion curves of the one-degree-of-freedom model, which shows the acoustic and optical wave modes as well as the bandgap interval. The values of the selected parameters are M=1.56,I⌢=1.5,K⌢=0.35,andω⌢n=8.9. Considering these parameters, the bandgap of the one-degree-of-freedom model is placed between the frequencies of 6.15 to 11.94 (Hz), which is the same as the frequency range obtained from Figure 11.


## 4.1.2. The Effect of System Parameters in the 1DOF Model

The surfaces drawn in Figure13 exhibit the effect of different dimensionless parameters on the starting and ending frequencies of their corresponding bandgaps. The top and bottom surfaces indicate the starting and ending frequencies, respectively, and the distance between them indicates the bandgap interval of the 1DOF model. As can be seen in Figure 13, there is a direct relationship between the bandgap and the dimensionless stiffness, meaning that as the dimensionless stiffness increases, the bandgap is also enhanced. On the other hand, by increasing the dimensionless mass, the bandgap decreases. It is deduced that the stiffness parameter has a positive effect on the expansion of the bandgap and the effect of the mass parameter is to weaken the bandgap interval. Moreover, according to the illustrated results, in general, the dimensionless moment of inertia directly affects the bandgap range, which means that as the dimensionless moment of inertia increases, the bandgap of the model is also enlarged. Although the bandgap interval slightly decreases at the lower values, when the dimensionless moment of inertia takes the higher values, this interval is then enhanced satisfactorily.Figure 13
The effect of changing (a)K⌢ and M and (b) K⌢ and I⌢on the starting and ending frequencies of the bandgap for 1DOF model.
(a)(b)


## 4.2. 2DOF Model of EMM

### 4.2.1. Dispersion Curve for the 2DOF Model

In order to display the dispersion curves of the modified model with multiple resonators, one can solve Equation (22) for ω and extract the desired curves. The values of the selected parameters are m¯1=1.04,m¯2=0.64,I¯=1.5,andK¯=2.86,ωn1=4.2rad/s,ωn2=15.9rad/s. Considering the above parameters, the bandgap starting frequencies (Fs1and Fs2) are calculated by substituting qa=±π in the acoustic and optical wave modes, respectively. The resulting plots are shown in Figure 14 by Ds(1) and Ds(2), respectively. In addition, the bandgap ending frequencies (Fe1and Fe2) are also extracted by replacing qa=0 in two optical roots of Equation (22), respectively. The plotted results are illustrated in Figure 14 by Ds(2) and Ds(3), respectively. As indicated in Figure 14, the difference between the first starting frequency and the first ending frequency demonstrates the first bandgap, and the difference between the second starting frequency and the second ending frequency represents the second bandgap. According to the illustrated results in Figure 12, one can find that there is a small area between the first and second bandgaps. This consequence is also demonstrated in Figure 14 with a slight difference between the highest and lowest points of Ds(2). The slight difference leads to the flattening of the Ds(2) curve.Figure 14
Dispersion curves for the 2DOF model.


### 4.2.2. Comparison of Dispersion Curves for Two Models

As shown in Figure15, by adding the concentrated mass and modifying the 1DOF model to a 2DOF model, two bandgaps are created in the dynamics of the new model, which allows us to reduce or suppress the vibrations of the modified system in more frequency ranges. From the comparison between the dispersion curves of both models, one can conclude that the first bandgap in the two-degree-of-freedom model starts and ends at lower frequencies than that of the one-degree-of-freedom model.Figure 15
Comparison of dispersion curves in two models.


### 4.2.3. The Effect of System Parameters in the 2DOF Model

Figure16 displays well that when the parameter m¯1 shifts upward, the first starting frequency also increases and the first ending frequency reduces, which satisfactorily expands the first bandgap of the modified model. On the other hand, by increasing the first bandgap, the second one decreases and the starting and ending frequencies eventually reach a constant value.Figure 16
The effect ofm¯1 parameter on bandgap for the 2DOF model.With increasing the parameterK¯, the first bandgap can be extended to cover a larger frequency range, as shown in Figure 17, although the first starting frequency will eventually approach a constant value. However, the change in the second bandgap of the model does not show the same trend, eventually converges to an almost constant value, and no longer changes significantly.Figure 17
The effect ofK¯ parameter on bandgap for the 2DOF model.For the sake of extending the first and second bandgaps of the modified model, the dimensionless moment of inertiaI¯ has to be increased. As indicated in Figure 18, by continuing the process of increasing the dimensionless moment of inertia, the first ending and the second starting frequencies will not change and take a constant value.Figure 18
The effect ofI¯ parameter on bandgap for the 2DOF model.As the first and second local frequencies increase, the same behavior cannot be observed in the first bandgap of this model, although it decreases with small values and follows an increasing trend at higher values (see Figures19 and 20). In general, it can be inferred that the first bandgap is expanded. Finally, it is found that with the increase of the first local frequency, the expansion of the second bandgap is not uniform, although, in general, it becomes larger. Moreover, it should be pointed out that when the second local frequency increases, the second bandgap of the modified model is reduced.Figure 19
The effect ofωn1 parameter on bandgap for the 2DOF model.Figure 20
The effect ofωn2 parameter on bandgap for the 2DOF model.


## 4.2.1. Dispersion Curve for the 2DOF Model

In order to display the dispersion curves of the modified model with multiple resonators, one can solve Equation (22) for ω and extract the desired curves. The values of the selected parameters are m¯1=1.04,m¯2=0.64,I¯=1.5,andK¯=2.86,ωn1=4.2rad/s,ωn2=15.9rad/s. Considering the above parameters, the bandgap starting frequencies (Fs1and Fs2) are calculated by substituting qa=±π in the acoustic and optical wave modes, respectively. The resulting plots are shown in Figure 14 by Ds(1) and Ds(2), respectively. In addition, the bandgap ending frequencies (Fe1and Fe2) are also extracted by replacing qa=0 in two optical roots of Equation (22), respectively. The plotted results are illustrated in Figure 14 by Ds(2) and Ds(3), respectively. As indicated in Figure 14, the difference between the first starting frequency and the first ending frequency demonstrates the first bandgap, and the difference between the second starting frequency and the second ending frequency represents the second bandgap. According to the illustrated results in Figure 12, one can find that there is a small area between the first and second bandgaps. This consequence is also demonstrated in Figure 14 with a slight difference between the highest and lowest points of Ds(2). The slight difference leads to the flattening of the Ds(2) curve.Figure 14
Dispersion curves for the 2DOF model.


## 4.2.2. Comparison of Dispersion Curves for Two Models

As shown in Figure15, by adding the concentrated mass and modifying the 1DOF model to a 2DOF model, two bandgaps are created in the dynamics of the new model, which allows us to reduce or suppress the vibrations of the modified system in more frequency ranges. From the comparison between the dispersion curves of both models, one can conclude that the first bandgap in the two-degree-of-freedom model starts and ends at lower frequencies than that of the one-degree-of-freedom model.Figure 15
Comparison of dispersion curves in two models.


## 4.2.3. The Effect of System Parameters in the 2DOF Model

Figure16 displays well that when the parameter m¯1 shifts upward, the first starting frequency also increases and the first ending frequency reduces, which satisfactorily expands the first bandgap of the modified model. On the other hand, by increasing the first bandgap, the second one decreases and the starting and ending frequencies eventually reach a constant value.Figure 16
The effect ofm¯1 parameter on bandgap for the 2DOF model.With increasing the parameterK¯, the first bandgap can be extended to cover a larger frequency range, as shown in Figure 17, although the first starting frequency will eventually approach a constant value. However, the change in the second bandgap of the model does not show the same trend, eventually converges to an almost constant value, and no longer changes significantly.Figure 17
The effect ofK¯ parameter on bandgap for the 2DOF model.For the sake of extending the first and second bandgaps of the modified model, the dimensionless moment of inertiaI¯ has to be increased. As indicated in Figure 18, by continuing the process of increasing the dimensionless moment of inertia, the first ending and the second starting frequencies will not change and take a constant value.Figure 18
The effect ofI¯ parameter on bandgap for the 2DOF model.As the first and second local frequencies increase, the same behavior cannot be observed in the first bandgap of this model, although it decreases with small values and follows an increasing trend at higher values (see Figures19 and 20). In general, it can be inferred that the first bandgap is expanded. Finally, it is found that with the increase of the first local frequency, the expansion of the second bandgap is not uniform, although, in general, it becomes larger. Moreover, it should be pointed out that when the second local frequency increases, the second bandgap of the modified model is reduced.Figure 19
The effect ofωn1 parameter on bandgap for the 2DOF model.Figure 20
The effect ofωn2 parameter on bandgap for the 2DOF model.


## 5. Concluding Remarks

In this paper, at first, a new model of elastic metamaterials (1DOF model) involving a rack-and-pinion mechanism was presented. Then, by adding a concentrated mass, the model was modified by converting to a novel metastructure with multiple resonators in which more broadband bandgaps were produced. Different system parameters allow us to shift or expand the frequency ranges at which the system vibrations are reduced or eliminated by changing them. The influences of different system parameters on the bandgap intervals of the considered metamaterial are summarized as follows:(i)
By adding the concentrated mass and converting the model from 1DOF to 2DOF, two bandgaps are satisfactorily produced in the dynamics of the system(ii)
In the 1DOF model, in order to achieve the wider bandgap, in general, the dimensionless parametersI⌢ and K⌢should be increased or the dimensionless mass parameter M should be reduced(iii)
In the 2DOF model, the width of the first bandgap is enhanced by increasing the parametersK¯, I¯, m¯1, and ωn1 and decreasing ωn2(iv)
In the 2DOF model, the width of the second bandgap is improved by increasingI¯, ωn1, and ωn2 or decreasing the parameter m¯1. No significant effect was demonstrated by varying the parameter K¯



---

*Source: 3914210-2021-06-02.xml*