Elastic Wave Properties of an Adaptive Electromechanical Metamaterial Beam

An adaptive metamaterial beam with a couple of electromechanical resonators in each unit cell is proposed in this study to open Bragg bandgap or locally resonant bandgap for flexural wave attenuation in Euler–Bernoulli beams. 0e electromechanical resonator is composed of a piezoelectric layer with segmented electrodes and shunt circuits, which affect the dynamic equivalent stiffness. It is illustrated that there is only a Bragg bandgap when the circuits of the two adjacent resonators are approximate to a short circuit or open circuit, and the locally resonant bandgap will be generated in the pure inductance circuits when the resonant frequencies are different in general. 0e locally resonant bandgap can be broadened by adding more resonators into the unit cell with the resonant frequencies of the shunting circuits satisfying a proper ratio.


Introduction
e research of acoustic/elastic metamaterials has lately attracted growing attention [1][2][3]. In particular, locally resonant metamaterials yield bandgap formation at wavelengths much larger than the lattice size of the substructure, and acoustic metamaterial has shown promising advantages in enabling vibration absorption and sound attenuation. In 2000, Liu et al. firstly used hard metal shperes embedded in a soft rubber as mechanical resonators to demonstrate the locally resonant bandgap properties.
is seminal work opened up a new way for the study of acoustic metamaterials.
A variety of metamaterials with purely mechanical locally resonators have been researched using the plane wave expansion method, transfer matrix method, or finite element method [5][6][7][8][9][10]. However, the resonant bandwidth yielded by elastic/mechanical metamaterials is narrow, which seriously limits their use in noise or vibration absorption. Moreover, actively controlling the position and width of the bandgap frequency region is very difficult in practice, if not impossible, for passive metamaterials.
One of the most pronounced challenges in elastic metamaterial development is the ability to tune their performance in an adaptive manner without requiring physical structural modifications. e piezoelectric shunting technique was originally proposed by Forward [11] and further investigated by Hagood and von Flotow [12]. ese electromechanical resonators have the benefit of significantly lower mass requirements than purely mechanical massspring resonators. e use of piezoelectric elements shunted to resonating circuits has lately been explored for locally resonant bandgap formation.
ese structures with resonating circuits can exhibit an electromechanical locally resonant bandgap.
Recently, Sugino et al. [13,14] derived the bound frequencies of the locally resonant bandgap for a piezoelectric bimorph beam, finding that the purely piezoelectric locally resonant bandgap size mainly depends on the electromechanical coupling coefficient. Xu et al. [15] established a semianalytical model of adaptive piezoelectric materials and proposed a new method to widen the bandgap by adding geometric cavities, which can increase the electromechanical coupling coefficient and effectively broaden the bandwidth.
A one-dimensional piezoelectric stack phononic crystal with a shunt circuit was studied by Chen et al. [16] and Lia et al. [17] using the transfer matrix method. eoretical and numerical results show that the shunt circuits can tune the regulation of structure flexural wave propagation. In addition, a hybrid adaptive metamaterial [18,19] with mechanical and electromechanical resonators is proposed to control the flexural waves and reveals that the phenomena for the mechanical and electromechanical locally resonant bandgaps are caused by the negative dynamic mass and negative dynamic stiffness, respectively. Two types of locally resonant bandgaps can be combined into a single, larger bandgap for increased bandwidth, by adjusting the natural frequency of the mechanical and electromechanical resonators.
In order to be suitable for more structural applications, two-dimensional electromechanical metamaterials with segmented electrodes and specified boundary conditions have been studied. Recently, to control the wave propagation in a thin plate, two-dimensional metamaterials with periodic arrays of shunted piezoelectric patches have been employed in Zhang et al. [20] by using the finite element method and Bloch theorem.
Kherraz et al. [21] demonstrated numerically and experimentally the opening of a locally resonant bandgap in an active phononic crystal and the application of periodic electroboundary conditions to control the dispersion of elastic waves in plates. By using Hamilton's principle, Sugino et al. [22] established a fully coupled two-dimensional governing equation and boundary conditions of the piezoelectric plate and found that the effective dynamic stiffness of the system depends on the admittance of the shunt circuit and electromechanical coupling in the system. By using the plane wave expansion method and the modal expansion method, the local resonant bandgap and wave attenuation have been obtained, and the finite element method was used to verify the analytical results. Metamaterials with piezoelectric shunting technology can also be used for active vibration control [23,24] and piezoelectric energy harvesting [25][26][27][28].
Inspired by Xiao et al. [29] and Zhou et al. [10], more mechanical resonators are added to the unit cell to broaden the bandgaps, especially for very low target frequencies. In this paper, a couple of electromechanical resonators containing a piezoelectric layer with segmented electrodes and shunt circuit in one-unit cell are proposed to open Bragg bandgap or locally resonant bandgap for flexural wave attenuation in Euler-Bernoulli beams. Band structures in both pure resistance and inductance circuits have been studied by using the transfer matrix method and the discretization method will be applied as numerical validation. e effects of the circuit parameters, connections, and length ratio on the band structures of this electromechanical metamaterial beam are investigated. e analytical and numerical results demonstrate that a couple of electromechanical resonators in each unit cell is practicable to open Bragg bandgap or locally resonant bandgap. Our design is of great significance for the practical applications of low-frequency vibration or wave attenuations in engineering. Figure 1(a) shows a simplified structural modal of our designed adaptive metamaterial beam, which is comprised of a host beam with periodically surfacebonded piezoelectric layers. Each of the piezoelectric layers with opposite polarizations in the transverse direction is made of two rectangular cross-section piezoelectric plates. e structure of each unit cell containing two subcells named subcell 1 and subcell 2 is illustrated in Figure 1(b). e symbols l 1 and l 2 represent the length of subcell 1 and subcell 2, respectively, h s and h p are the thicknesses of the host beam and piezoelectric layers, and v j1 , v j2 and Y j1 , Y j2 are the voltage and the external load admittance across the two subcells of the jth electrode pair. e host beam and the piezoelectric layer share the same width b. e circuits of the subcells 1 and 2 have the same connections and different impedances, which causes different dynamic effective modulus and different natural frequencies. It will be found that the periodic arranged piezoelectric layers and alternate impedances give rise to various phenomena of band structures.

Model and Method of Calculation
In the following two subsections, the transfer matrix method and modes discretization method are used, respectively, to infinite and finite modes to validate the results of each other.

Propagation of Flexural Waves in Electromechanical Metamaterial Beam by Using the Transfer Matrix Method.
e mechanical interface boundary effect between the host beam and the piezoelectric layers is ignored since the bonding is assumed to be perfect. e constitutive equations of the piezoelectric beam are where S and T are the normal stress and normal strain in the x-direction and D and E 3 denote the electric displacement and electric field in z-direction, respectively. e symbol s E 11 is the elastic compliance at the constant electric field, d 31 is the piezoelectric strain constant, and ε T 33 is the permittivity component at constant stress. Here, for the modeling of the piezoelectric metamaterial, we adopt the equivalent Young's modulus representation [12] where Z j1 (ω) and Z j2 (ω) are the impedance of the two subcells, depending on the type of the shunt circuit and the value of the circuit element. e impedance Z j and the admittance Y j mentioned in Figure 1 are reciprocal of each other Y j � 1/Z j . In pure resistance/inductance circuit, the impedance can be expressed as R j /iωL j . e circuit can be approximate to open/short circuit when the impedance is infinite or goes to zero indefinitely. e short-circuit Young's modulus E p , electromechanical coupling coefficient k 31 , the permittivity component at constant strain ε S 33 , and the capacitance of the subcell 1 and subcell 2C p,j1 and 2C p,j2 are, respectively, In order to study the complex band structure of electromechanical metamaterial beams, the traditional transfer matrix method is used. Based on the Euler-Bernoulli beam theory, the governing equation of the transverse displacement for the host beam segment can be written as where In equations (5)-(7), w(x, t) denotes the transverse displacement of the host beam at position x and time t. Symbols E s , ρ s , I s , and A s are Young's modulus, mass density, area moment of inertia, and area of cross section of the host beam, respectively. Symbols E p , ρ p , I p , and A p are Young's modulus under constant electric field, mass density, area moment of inertia, and area of cross-section of the piezoelectric plate, respectively.
To characterize the wave-attenuation properties of the metamaterial beam, the harmonic solution to equation (5) can be written as where W(x) is the mode shape function and ω is the frequency of vibration. e general solution of subcell 1 and subcell 2 for the jth unit cell can be written as where A j , B j , C j , and D j are unknown constant and Shock and Vibration 3 Consider the continuity conditions of the displacement, rotation angle, bending moment, and shear force at the connecting point between the two subcells: By substituting equations (9) and (10) into equation (12), we can express the continuity conditions in the matrix form as where and In equations (15) and (16), Similarly, applying the continuity conditions of the displacement, rotation angle, bending moment, and shear force between the two consecutive unit cells, we have Each unit has the same material structure parameters and circuit parameters and the equivalent Young's modulus E p,j1 � E p,(j+1)1 and E p,j2 � E p,(j+1)2 . By substituting equations (9) and (10) into equation (19), we express the continuity conditions in the matrix form as where Combining equations (13) and (20) and eliminating ψ j2 yield e transfer matrix between the two cells is then defined as Further, by employing the Floquet-Bloch theorem, for a one-dimensional periodic structure, the relationship between ψ j+1 and ψ j is where q is the wave vector along the x-direction. By inserting equation (24) into equation (22), a standard eigenvalue problem is formed where I is the 4 × 4 identity matrix. e dispersion relations of the metamaterial beam with electromechanical resonators can be obtained by solving the above dispersion equation (25). It can be seen that the dispersion equation is related to the equivalent Young's modulus mentioned in equations (2) and (3), which is affected by the circuit elements, including resistance, inductance, and inherent capacitance. e effects of short circuit and open circuit will be taken into account in further discussions. Based on the given frequencies, the complex wavenumbers are determined, from which the wave-propagation properties are derived. e propagation of a flexural wave in the beam is shown bandgap or passband when the solution of q is imaginary or real, correspondingly.

Modes Discretization of the Electromechanical Metamaterial Beam with a Finite Length.
For the purpose of the practical design, it is necessary to investigate the frequency response of the metamaterial beam with a finite number of electromechanical resonators. A uniform Euler-Bernoulli beam with piezoelectric layers is considered. e metamaterial beam with transverse base excitation w b (t) is of length L, and 2S electromechanical resonators are periodically placed onto the host beam at the distance of l 1 and l 2 and l 1 and l 2 satisfying l 1 + l 2 � 2l. Without loss of generality, we assume the host beam and the resonators are all undamped and the length of the two subcells is the same, l 1 � l 2 � l. e governing equation for the electromechanical metamaterial beam is then e governing equations for the shunt circuit of the two subcells are where w(x, t) is the transverse displacement of the beam at position x and time t and v j1 and Y j1 are, respectively, the voltage and the external load admittance across the subcell 1 of the jth electrode pair. v j2 and Y j2 are parameters of subcell 2. H(x) is the Heaviside function, and f(x, t) is the external forcing. Furthermore, EI and m are the short-circuit flexural rigidity and mass per length of the beam motioned in equation (1), ϑ is the term associated with electromechanical coupling coefficient, and C p, j1 and C p, j2 are the internal piezoelectric capacitance across the subcell 1 and subcell 2 of jth electrode pair, given by Shock and Vibration 5 We assume a modal expansion [30] for the transverse displacement of the beam w(x, t) where functions ϕ r (x) are the normalized mode shape functions of the plain beam (i.e., without resonators) at short circuit, η r (t) are the generalized coordinates to be confirmed, and N is the discretized number. e modal functions of the host beam satisfy the orthogonality conditions From equations (27) and (28), we obtain where V j1 and V j2 are the voltage amplitudes of the two subcells. Substituting equations (32) and (33) into equation (26), multiplying by ϕ n (x), and integrating over the beam length from 0 to L and then using the orthogonal conditions, the discretized governing equation can be obtained where By solving equation (34), the values of generalized coordinates can be calculated; substituting them into equation (30), the relative deflection amplitude can be obtained. e frequency response functions of the metamaterial beam are defined and calculated as In summary, we proposed two methods to treat the dispersions: the transfer matrix method based on the state vectors and the continuity conditions between the two subcells and the units and the discretization method based on the vibrational point of view. e band structures of the proposed periodic structures will be derived by the transfer matrix method and the discretization method will be applied as numerical validation. Hence, equations (25) and (36) will be used in the following sections.

Dispersion Relations of the Adaptive Metamaterial Beam
In this section, the complex band structures of the metamaterial beam with various values of the shunt circuit elements containing resistance and inductance will be discussed based on the transfer matrix method. e parameters used in this study have been listed in Table 1.

Bragg Bandgap in Pure Resistance
Circuit. In a pure resistance circuit, there are only the resistance elements and the inherent capacitance of the piezoelectric layers in the circuits which cannot generate resonant frequency. It is shown when the resistance values of the two circuits differ greatly and can be approximated as an open circuit and short circuit, the structure can only generate Bragg bandgap, and the band characteristics are most obvious. Figure 2 shows the first two-order Bragg bandgaps by R 1 � 10 −6 Ω and R 2 � 10 6 Ω. e Bragg bandgaps in pure inductance circuit mainly reflected in the first-and second-order bandgaps and the bandwidths are generally narrow. e impedance here rests with the resistance values of the external circuit. It has little effect on the band frequency region and bandwidth of the Bragg bandgap except the amplitude of the propagation constant.

Dispersion Relations in Pure Inductance
Circuit. e band structure caused by the inductance circuit is different from that caused by the resistance circuit. In a pure inductance circuit, the circuit contains an inductance element, which forms an LC oscillation circuit with the inherent capacitance of the piezoelectric layer. e impedances mentioned in equations (2) and (3) are iω e1 L 1 and iω e2 L 2 , where ω e1 and ω e2 are the oscillation frequencies of the two circuits. Due to the existence of the oscillation frequency, the metamaterial beam containing a pure inductance circuit can produce both locally resonant bandgap and Bragg bandgap. Moreover, the resonant bandgap frequencies are generally lower than Bragg bandgap frequencies. 6 Shock and Vibration

Bragg Bandgap.
In a pure inductance circuit, the addition of inductance will also cause the impedance mismatch between the two subcells. In addition, the oscillation frequency generated by the circuits will also affect the Bragg bandgap and the locally resonant bandgap. Similar to the case of the pure resistance circuit, the Bragg bandgap appearing under the condition of the shunt circuit connected to L 1 is equivalent to an open circuit, while the shunt circuit connected to L 2 is equivalent to a short circuit. Figures 2(a) and 2(b) show the same phenomenon that the first two Bragg bandgaps with L 1 approach infinity and L 2 approach zero. e locally resonant bandgap generated by L 1 is infinitesimal which can be ignored and that generated by L 2 is infinite, and the system is mainly manifested as Bragg type. Moreover, the imaginary part value of the propagation constant is also relatively small.

Locally Resonant Bandgap.
In general, short and open circuits are not considered, the resonant frequency generated by the inductance circuit is limited, and the impedance difference between the subcells is not significant. In this case, the electromechanical metamaterial beam mainly exhibits resonant bandgap in a pure inductance circuit. Next, the effect of the inductance value on locally resonant bandgap is discussed.
Considering the particular case L 1 � L 2 , the impendence and the natural frequency of the two subcells are exactly the same, and the model can be reduced to traditional electromechanical metamaterials. e corresponding local resonant bandgap is shown in Figure 3. It is noted that the bandwidth is narrower and the frequency range is lower than that of the Bragg bandgap. e bandwidth here seems to be relatively narrow because the electromechanical coupling of PZT-5H piezoelectric material is relatively weak. If the piezoelectric material is replaced by that with higher electromechanical coupling coefficient, the bandwidth here will be increased. Figure 4 shows the band structures of the electromechanical metamaterials beam for different values of the inductance. It can be observed that two separate local resonant bandgaps appear due to the inequality of the inductance values. As we all know, the larger the inductance, the smaller the natural frequency of the corresponding circuit. So the resonant bandgap caused by L 1 has a higher frequency and a wider frequency region than that of L 2 .
In order to predict the influence of inductance on locally resonant bandgap, the assumption of L 1 � 10 H is given and the bandwidth with respect to various values of L 2 is discussed. According to Figure 5(a), the position and width of the bandgap generated by subcell 1 are determined with bound frequencies ω from 113.46 Hz to 117.03 Hz, and the bandwidth and frequency range are generated by subcell 2 decrease as L 2 increases. When L 2 increases to a certain value, the local resonant bandgap generated by subcell 1 and the bandgap generated by subcell 2 will be infinitely close to  each other. erefore, determining appropriate inductance or effective modulus can be employed to achieve a single lower and wider bandgap. As shown in Figure 5(b), it is theoretically possible to perfectly combine the two bandgaps by setting L 2 � 10 H and L 1 � 9.3573 H, and the frequency ranges are exactly adjacent, creating a single and wider bandgap. For the structure proposed in this paper, wider and lower bandgap structures can be obtained. In general, when the structure contains multiple subunits, the bandgap will be extremely wide and low.

e Influence of length Ratio on Dispersion Relations.
In this section, the effect of the length of the two subcells on band structure will be discussed. Firstly, the coefficient η is introduced which can be written as η � l 1 /(l 1 + l 2 ), where l 1 and l 2 are the length of the subcell 1 and subcell 2.

Bragg Bandgap in Pure Inductance
Circuit. Based on the results given in Sections 3.1 and 3.2, the Bragg bandgap is most obvious when the circuits connected to the two subcells are open and short, respectively. Here, the resistance values of the two subcells circuits are kept at constant R 1 � 10 6 Ω and R 2 � 10 −6 Ω, respectively. Figure 6 shows the Bragg bandgaps with various η. It should be noted that, with the increase of the length ratio η, the location of the first Bragg bandgap moves toward a higher frequency and the width of the Bragg bandgap varies to wider. e capacitance of the piezoelectric material described by equations (36) and (37) affected by the length of the two subcells can change the impendence and Bragg bandgap. Coincidentally, the phenomenon in the pure inductance circuits is similar to the pure resistance circuits.

Locally Resonant Bandgap in Pure Inductance Circuit.
e capacitance of the two subcells is affected by the length. For the given inductance L 1 � L 2 � 10 H, Figure 7 shows the locally resonant bandgaps of the electromechanical metamaterial beam for various length ratio η. Comparing the four plots, it is found that the bandgap generated by subcell 1 with a higher frequency region will move toward lower frequency and the width will vary to wider with the increase of η. On the contrary, the bandgap generated by subcell 2 with lower frequency region will move toward a higher frequency and the width will vary to narrower. In addition, the locations of the two bandgaps are getting closer to each other. e two bandgaps are coincident consistent with the traditional piezoelectric metamaterials as shown in Figure 7(d) when η � 0.5. e bandwidth and frequency range of the locally resonant bandgap mainly depend on the value of η when the inductance value is constant. Next, the relationship between the bandwidth of two locally resonant bandgaps and η is analyzed. e bandwidth of the two subcells versus the length ratio η is depicted in Figure 8(a). It can be observed that the bandgap generated by subcell 1 is shifted from highfrequency range to a lower one as η increase, and the bandwidth becomes wider. In contrast, the position and width of the bandgap generated by subcell 2 are opposite to the bandgap generated by subcell 1. It also can be seen that there are two values near both sides of η � 0.5 which can widen the bandgap. As shown in Figure 8(b), there are similar phenomena to L 1 � L 2 � 10 H, but no symmetry due to the unequal of the inductance values. In addition, there are two length ratios near the location where the two frequency ranges overlap which can tune the two bandgaps to be adjacent, creating a single, wider bandgap.
According to Sugino et al. [13], the starting and ending frequency of the local resonant bandgap with electromechanical resonators can be given by where α is the electromechanical coupling coefficient, which measures how strongly coupled the piezoelectric material is to the host beam. en it can be written as Based on equation (38), one can easily find that the electromechanical coupling coefficient is affected by the geometric and structural parameters of the metamaterial beam. To obtain a wider locally resonant bandgap, the frequencies of the electromechanical resonators should be selected as ω s1 � ω e2 , and the bandgap can be widened at η � 0.328654, where the ending frequency of the bandgap generated by subcell 2 and the starting frequency of the bandgap generated by subcell 1 coincide as seen in Figure 9.

Numerical Results
e frequency response functions in an elastic metamaterial beam can examine the wave properties effectively. However, for a finite length of the electromechanical metamaterial beam, the number of unit cells and the order of truncation will always affect the efficiency of wave attenuation. In order to explore the influence of various circuit elements on wavepropagation characteristics, the numerical simulations of wave attenuation with twenty unit cells and eight orders truncation are conducted based on equation (36). In the calculations, the geometric and structural parameters are also the same as listed in Table 1, and the length of metamaterial beam and each subcell is L � 1 m and l � 25 mm.

Bragg Bandgap in
Open-Short Circuit. In pure resistance circuit or pure inductance circuit, the governing equations (27) and (28) for shunting circuits of the two subcells have the same form except the admittances Y j , where 1/R 1 and 1/R 2 are for resistance circuit and 1/iω e1 L 1 and 1/iω e2 L 2 are for inductance circuit. As can be seen from Section 3.2, the Bragg bandgaps in the pure resistance circuit or the pure inductance circuit are uniform. e frequency responses of the electromechanical metamaterial beam for open-short circuits are shown in Figure 10. It  shows that the first Bragg bandgap shown in Figure 2 and the current frequency responses are in good agreement with each other. Figure 11(a) shows the frequency responses of the traditional metamaterial beam when the piezoelectric layers are connected with pure inductance circuits by setting L 1 � L 2 � 10 H. Here, geometric and structural parameters are the same as those used in Figure 10. As illustrated in Figure 11(a), a single narrow bandgap is observed at frequencies between 112.71 Hz and 117.10 Hz. is transmittance dip and the bandgap frequencies predicted by the analytical dispersion relations shown in Section 3 are in good agreement with each other. Figure 11(b) shows the frequency responses of the electromechanical metamaterial beam with L 1 � 10 H and L 2 � 5 H, and it is clear from the figure that there are two narrow separate bandgaps whose frequencies are between 113.22 Hz and 117.10 Hz and the other frequencies are between 159.52 Hz and 165.81 Hz. In order to achieve the desired performance demonstrated in Figure 5(b), the inductance of the subcell should be modified slightly. Specifically, for the pure inductance circuit, L 2 is selected as 9.3573 H. In this case, the ending frequency of the bandgap generated by subcell 2 and the starting frequency of the bandgap generated by subcell 1 coincide. As shown in Figure 11(c), the frequency responses of the electromechanical metamaterial beam with L 1 � 10 H and L 2 � 9.3573 H demonstrate a single and wider bandgap from 113.11 Hz to 121.79 Hz. It should be noted that a transmission peak exists inside the widen bandgap, which is caused by the coupling between the two circuits of the adjacent subcells. In the future, geometric optimizations are still in demanding in attenuation performance.

Locally Resonant Bandgap in Pure Inductance Circuit.
As Bragg bandgap can only be obtained when the circuits of the two adjacent resonators are approximate to short circuit and open circuit, the existing formula (i.e., equation (37)) cannot predict the Bragg bandgap. Table 2   structure analysis, transmittance, and the existing formula (i.e., equation (37)). e parameters for calculating the results of Table 2 can be obtained from Table 1. e length ratio is taken as 0.5, and the electromechanical coupling coefficient of the two subcells α is 0.0642. It can be seen from Table 2 that the starting frequency of the bandgap obtained by using the band structure analysis is basically the same as that obtained by the existing formula (i.e., equation (37)), and the ending frequency of the bandgap obtained by using the band structure analysis is slightly smaller than that obtained by using the existing formula. Due to the difference in the computational accuracy between the band structure analysis and the existing formula, the bandgap obtained is slightly different. In addition, when considering twenty unit cells and eighth order truncations, the transmittance obtained by equation (36) has a high accuracy, which can be used to verify well in the band structure. In general, the band structure analysis, transmittance, and the existing formula (i.e., equation (37)) are in good agreement with each other in bandgap prediction.

Conclusions
In analogy with the mechanical metamaterial, electromechanical metamaterial made from elastic substrates with piezoelectric layers shunted to shunting circuits also exhibits bandgaps. However, the bandgap formation in electromechanical metamaterial is associated with a frequency-dependent modal stiffness and the frequency-dependent modal mass in the mechanical metamaterial. Here, an adaptive metamaterial beam with a couple of electromechanical resonators in each unit cell is proposed to open multiple locally resonant bandgaps or Bragg bandgaps. We presented a detailed analytical and numerical investigation of the metamaterial beam with electromechanical resonators. e numerical results show that the Bragg bandgap of the system is the most obvious when the pure resistance circuit of two subcells in a unit can be approximated as open circuit and short circuit, respectively. Similarly, when the pure inductance circuits of the two subcells are approximate to open and short circuit respectively, the system also presents mainly Bragg bandgap characteristics. Otherwise, multiple separate locally resonant bandgaps can be created. e frequency responses as numerical validation are in good agreement with the analytical bandgaps. It is concluded that the length ratio η can also affect the band structures dramatically by varying the capacitance of the subcell. Further, the local resonant bandgap can be broadened, creating a single, wider local resonant bandgap, by adding two resonators to the unit cell with the resonant frequencies of the shunting circuits satisfying a proper ratio. In general, when the structure contains multiple subunits, the bandgap will be extremely wide and low. e system proposed in this paper opens up important perspectives for the development of adaptive vibration or wave-attenuation devices for broadband and low-frequency applications.

Data Availability
All the data used in the study are available. Previously reported data were used to support this study and are available.
ese prior studies are cited at relevant places within the text as references.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.