1. Introduction
Earthquakes are one of the most hazardous natural disasters in the world, which release seismic energy through seismic waves and pose a huge hazard to buildings [
1]. The seismic waves are combined with surface waves and body waves; Rayleigh waves are vertically elliptical surface waves traveling along the ground surface, which has been demonstrated to pose potential damage to buildings. Simultaneously, the surface wave barriers periodically arranged around the building can protect buildings by attenuating the effects of Rayleigh surface waves.
In recent years, scholars have conducted numerous research on seismic metamaterial surface waves barriers using phononic crystals [
2,
3,
4,
5,
6]. Surface waves barriers can be classified into two types based on the mechanisms of bandgap generation: Bragg scattering [
4,
7] and local resonance [
8,
9,
10,
11,
12,
13]. Bragg scattering type barriers can generate one wide bandgap, effectively attenuating most frequencies within 0–20 Hz. Muhammad et al. [
14] designed two 2D columnar barriers with lengths of 12.5 m and calculated the bandgap for both single and multilayered soils. The results showed that the barriers can attenuate surface waves in both types of soils, offering insights for soil conditions in practical implementation. Witarto et al. [
15] analyzed a 3D square cube buried in soil and discovered that a square cube with a period constant of 1 m can generate an elastic waves bandgap from 18.44 Hz to 27.09 Hz. Liu et al. [
16] designed a continuous wall that was partially buried underground, which can attenuate elastic waves in the ranges of 6.2–8.8 Hz and 15.1–20 Hz. The attenuation effect of barriers with a positive gradient was better than that of a periodic arrangement, which offers a new method of barrier arrangement. Wang et al. [
17] conducted numerical simulations and scaling tests on both the petal-shaped-filled barriers and ordinary square-filled barriers. The results showed that the petal-shaped-filled barriers can generate multiple bandgaps in the range of 0–20 Hz and their effect is superior to that of the ordinary square-filled barriers. Numerous studies [
18,
19,
20,
21,
22,
23,
24,
25,
26,
27,
28,
29] demonstrated that the Bragg scattering barrier effectively attenuates Rayleigh surface waves. However, the sizes of the most Bragg scattering barriers are too large to apply and implement in practical engineering.
Many studies have been conducted on seismic surface wave resonators based on the theory of local resonance. These resonators were much smaller in size compared to Bragg scattering barriers, making them more beneficial for engineering implementation. Liu et al. [
30] proposed a local resonance metamaterial composed of rubber wrapped around a lead ball and arranged in squares. This design achieves the goal of controlling a large wavelength using a small size, providing a new idea for the design of seismic surface waves barriers. A resonant meta-wedge [
31] was designed to convert surface waves into body waves and to reflect Rayleigh surface waves, and an analytical model was developed to accurately predict the rainbow effect of earthquakes. Some scholars [
32,
33,
34,
35] calculated the bandgap of local resonant resonators using the finite element method and concluded that these resonators were capable of generating a low-frequency bandgap, which was helpful for attenuating Rayleigh waves. Some references [
36,
37,
38] theoretically deduced the attenuation properties of resonators and numerically proved that the theoretical model can achieve an efficient low-frequency bandgap. Yan et al. [
39] designed a cubic resonator foundation made of three materials and demonstrated that it can effectively attenuate seismic surface waves by theoretical analysis and field experiments. Palermo et al. [
40] designed a 2D super-barrier and showed that it can convert surface waves into body waves to achieve the attenuation of Rayleigh waves. The finite element method showed that the metabarrier was able to attenuate surface waves by converting surface waves into body waves. Additionally, a scaling experiment was conducted to verify that the resonator was indeed able to attenuate surface waves.
In recent years, machine learning has been developed and has achieved certain results in various fields. Wang et al. [
41] combined the comparative learning pretraining method with the self-supervised pretraining mode to learn the internal features of the image, and the prediction results were remarkable compared with the traditional methods. Purohit and Dave [
42] compared six algorithms and suggested that combining self-supervised learning and weakly-supervised learning to form a self-supervised learning model can better analyze labeled and unlabeled datasets and produce more accurate image segmentation models. Gaur et al. [
43] used a neural network to calculate the deformation of elastic–plastic structures and derived the stress–strain relationship, which is consistent with the linear theory and shows the accuracy and effectiveness of the method. Liu and Lu [
44] used particle swarm optimization and 10-fold cross-validation to analyze the parameters and their predictions matched the experimental results, indicating that the computational framework is useful for parametric analysis of materials. Koy and Colak [
45] provided better prediction of data with short-term variations by combining artificial neural network and vector autoregression. A review of the literature on machine learning revealed that it can be applied to the analysis of dispersion curves of resonators. This enables the estimation of the amplitude of vibrations generated by the resonator and the soil under different conditions. Furthermore, machine learning can be used to optimize the parameter design of the resonator and to analyze the elastic dynamics. These will be studied in detail later.
As mentioned above, theoretical calculations indicate that the local resonance resonator can generate ultra-low-frequency elastic waves bandgap with a smaller size compared to the Bragg scattering barrier, which makes it more favorable for engineering applications. However, this research on local resonance type structures is insufficient, especially the single form, such as only two to three materials wrapped into a spherical structure or arranged into a layered structure. Therefore, a new resonator with a three-vibrator is proposed in this study, which can absorb the energy of elastic waves by the internal vibrators of the resonator to attenuate the energy of the elastic waves. Then, fewer studies have been conducted on seismic surface waves in the range of 0–20 Hz. In this study, we focus on the resonator for the most destructive Rayleigh waves, which can be attenuated effectively from the range of 0 Hz to 20 Hz. Therefore, the resonance effectively attenuates the 0–20 Hz seismic surface waves, which extends the research on bandgap generation by local resonance theory.
This study is organized as follows.
Section 2 proposes a theoretical model of a resonator with a three-vibrator based on the local resonance theory and establishes the dynamic equilibrium equation.
Section 3 introduces the theoretical calculation method of the finite element model.
Section 4 calculates the bandgap of the novel resonator and analyzes the mechanism of bandgap. In
Section 5, the transmission spectrum is calculated and the displacement modes are analyzed for the mechanism of bandgap. In
Section 6, time history analysis is used to verify the bandgap and attenuation effect of the resonator on seismic surface waves.
Section 7 analyzes the effect of resonator size parameters on the bandgap range.
Section 8 analyzes the effects of soil and material damping on attenuation characteristics for Rayleigh waves. The main conclusions are presented in
Section 9.
2. Theoretical Modeling and Development of Equations of Motion
A theoretical model of resonator with three vibrators is proposed based on the local resonance theory [
30,
31,
34], as shown in
Figure 1. The model consists of an outer frame and three internal vibrators. The internal vibrators and the outer frame are connected by springs and the resonators are also connected by springs. Moreover, the springs between the resonators represent the propagation of vibration displacement, while the springs between the outer frame and the internal vibrators represent the vibrational displacement of the internal vibrators. The vibrations generated by the internal vibrators are used to absorb the energy of the elastic waves, thus achieving the purpose of attenuating the elastic waves. The parameters of the resonator are shown in
Table 1. The materials are all linearly elastic and undamped.
The springs are used to connect the outer frame of the resonator to the inner vibrators. When the elastic waves propagate, the energy carried by the elastic waves is transferred to the outer frame m0 of the resonator through the spring k0, and then to the internal vibrators m1, m2, and m3 through the springs k1, k2, and k3, respectively, so the internal vibrators can vibrate. The internal vibrators absorb part of the energy of the elastic waves, thus attenuating the elastic waves. On the other hand, damping is an intrinsic property of the material that allows the absorbed energy to be dissipated.
Three theoretical methods can be used to calculate the bandgaps of structure, including plane wave expansion (PWE) [
4,
46], multiple scattering theory (MST) [
10,
47,
48,
49,
50], and mass-spring theory [
37,
38,
51,
52,
53,
54]. The mass-spring theory is the most appropriate method according to the characteristic of the resonators. Thus, as shown in
Figure 1, the dynamic equilibrium equations of the resonator can be expressed as follows:
where
m0 is mass of the outer frame.
m1,
m2, and
m3 are the internal vibrators, respectively.
k0 is the spring stiffness between the resonators.
k1,
k2, and
k3 are the spring stiffness between the resonator and the internal vibrators, respectively.
is the displacement of the
jth resonator.
is the displacement of the
ith vibrator in the
jth resonator.
and
are the acceleration of the
jth resonator with respect to time and that of the
ith vibrator in the
jth resonator with respect to time. The displacement function is shown as
where
u0 and
un are the displacement of the outer frame and the
nth internal vibrator, respectively.
U0 and
Un are amplification of the outer frame and the
nth internal vibrator, respectively. ω is the circular frequency,
t is the time variable, and i
2 = −1. The periodic boundary conditions of the resonator are shown as
where
and
are the displacement of the (
j − 1)th and (
j + 1)th resonator.
k is the wave number and
a is the lattice constant. Therefore, the dispersion equation of resonator can be calculated by Equations (1)–(3) and shown as
Equation (4) can be simplified as
where
is the fourth-order matrix of the wave vector.
is the fourth-order matrix of the displacement.
is a fourth-order unit matrix. Therefore, the dispersion curve of resonator can be calculated by Equation (4) and
Table 1 and shown in
Figure 2. It can be seen from
Figure 2 that the theoretical model of the resonator can generate a low-frequency wide bandgap from 2.7 Hz to 4.5 Hz. It can be obtained that the resonator produces one low-frequency wide bandgap when the mass and the spring stiffness of the internal vibrator are the same. Therefore, the theoretical model of the resonator is useful for the attenuation of elastic waves, which provides a theoretical reference for the establishment of the practical model of the resonator.
3. Numerical Models
The finite element model of the resonator is established according to the theoretical model of multi-vibrator resonators. Based on local resonance theory, we select the internal vibration method to design a 2D resonator comprising a copper outer frame and three copper internal vibrators connected by rubber bars. Some previous studies [
14,
17,
20,
21,
35] have shown that the hollow square columns can produce a wider low-frequency bandgap, so a hollow square column is selected for the outer frame of the resonator in this study.
Figure 3 shows the arrangement and dimensions of the resonator. Huang and Shi [
23] pointed out that cylindrical barriers arranged in a hexagonal shape have a larger filling ratio compared to a square arrangement, which is more advantageous for a low-frequency wide bandgap. However, the literature [
23] focuses on solid cylindrical barriers, which generate a bandgap by Bragg scattering, and the wide bandgap is a directional bandgap that attenuates elastic waves in one direction only. The resonator proposed in this paper generates bandgaps by the mechanism of local resonance and the bandgaps are omnidirectional bandgaps that can attenuate elastic waves in all directions. Therefore, the attenuation effect can be achieved by using an arrangement perpendicular to the building, as shown in
Figure 3a. The square arrangement is chosen for simplicity of calculation, and the hexagonal arrangement of the resonator will be calculated in a subsequent study.
The lattice constant is denoted by
a, while
a1 and
a2 represent the outer and inner length of the square frame, respectively.
hb represents the depth of the resonator buried in the soil. The spacing between the three internal vibrators and the internal vibrator radius are represented by
b and
r. Simultaneously, the side length of the rubber connecting the internal vibrators to the outer frame is represented by
a3. The height of the unit cell is denoted by
H. It has been demonstrated that when the height of the soil is greater than 10 times of the lattice constant, it has no significant effect on the bandgap created by the structure [
55]. The dimensions and material parameters of the barrier are presented in
Table 2 and
Table 3, respectively.
A surface waves resonator structure is created by arranging a unit cell periodically along the
x–y direction. Firstly, the periodically arranged resonator can be simplified to a square unit cell based on phononic crystal and periodicity theory, and the first Brillouin zone is shown in
Figure 3d. Then, the eigenfrequency of the unit cell can be calculated by scanned the path Γ-Χ-Μ-Γ in the first Brillouin zone in the eigenfrequency module of the finite element software COMSOL Multiphysics 5.6.
The material is defined as linearly elastic and undamped, and soil properties are consistent with Biot’s theory and are undamped. The governing equation used in the calculation of the eigenfrequencies can be written as
where
is the position vector,
u is the displacement vector,
ρ is the mass density, and λ and μ are the elastic constants. According to the Floquet–Bloch theorem, in the present periodic system, the solution of Equation (6), can be written as
where
uk(
r) is modulation function with the same periodicity as the elementary unit cell. The periodic boundary condition can be set around the unit cell and written as
where
a is the constant vector, which gives the repetition of a periodic structure. Combining the governing equation of Equation (6) and the periodic boundary condition of Equation (8), the dispersion relation of an infinite periodic system can be transferred into an eigenvalue problem of an elementary unit cell. The dispersion equation is an implicit function of wave vector
k and angular frequency ω in mathematics,
where
K and
M are the stiffness and mass matrices of the unit cell, respectively.
To eliminate the effect of structural dimensions, the relative bandgap can be used and written as
where
fu and
fl are the upper and lower bandgap frequencies, respectively.
4. Bandgaps of the Resonator
Rayleigh waves are a type of seismic surface waves, which carry most of the energy of the generated elastic waves within a frequency range of 0–20 Hz. The resonator designed focuses on attenuating Rayleigh waves within the frequency bandgap. To establish the unit cell model, the periodic boundary conditions are set around the unit cell to simulate the periodically arranged resonator. The fixed boundary conditions can be set at the bottom of the unit cell to simulate the soil in a semi-infinite space. Then, the eigenfrequency of the structure can be obtained by scanning the path Γ-Χ-Μ-Γ in the first Brillouin zone. However, the eigenfrequency diagram includes both body waves and surface waves. The previous research shows that the surface waves can be distinguished using the energy identification method [
56]. The energy identification method and the model presented in the Pu et al. [
56] was used to verify the accuracy of the finite element method. The model and the bandgaps are plotted in
Figure 4a [
56]. The periodic constant of the unit, the radius of pile, and the parameter about depth are taken as
a = 0.8 m,
r = 0.3 m, and
h = 20
a, respectively. The model comprises foam columns that are completely buried in the soil. The Young’s modulus
E, mass density ρ, and the Poisson ratio υ of foam and soil are taken as 37 MPa, 20 MPa, 60 kg/m
3, 1800 kg/m
3, 0.32, and 0.35, respectively. The two materials are also assumed to be linearly elastic and undamped. The COMSOL is used for the eigenfrequencies calculation and the energy distribution parameter identification method is used for surface wave identification. It can be seen from
Figure 4b that the present results are in agreement with the previous paper, which proves that the calculation method used in this study is correct.
Figure 5 shows the bandgaps for the multi-vibrator resonator. The first five bandgaps cover the frequency of the Rayleigh waves, and the ranges are 4.68–6.01 Hz, 6.03–9.31 Hz, 9.48–10.90 Hz, 10.94–18.72 Hz, and 18.79–22.01 Hz, respectively. All the bandgaps are omnidirectional and can attenuate Rayleigh waves in all directions.
The finite element model is equivalent to the theoretical model and its dispersion curve is calculated to verify the accuracy of the finite element modeling. Given that the theoretical model is two-dimensional, the cross-section of the finite element model is taken as the two-dimensional equivalent section. The masses of the internal vibrators of the resonator are the volume times the density, where the volume is the cross-sectional area times the unit length, as follows:
The mass of the outer frame of the resonator is the volume times the density, where the volume is the cross-sectional area times the unit length, as follows:
The spring stiffnesses between the resonators and between the resonator and its internal vibrators are provided by that of the soil and rubber, respectively. The stiffness formula is as follows:
where
En is the Young’s modulus of soil or rubber.
An is the area of the resonator in contact with the soil or rubber in the plane normal to the direction of spring action.
Ln is the length of soil or rubber equivalent to the spring. The equivalent parameters are shown in
Table 4. The equivalent model dispersion curves can be obtained by
Table 4 and Equation (4), as shown in
Figure 6. The bandgap of the equivalent model is from 4.36 Hz to 5.99 Hz in
Figure 6, while the first bandgap of the finite element model is from 4.68 Hz to 6.01 Hz. The errors in the upper and lower boundaries are minimal, amounting to 0.3% and 6.8%, respectively. The reason for the errors may be the different dimensions of the finite element model and the theoretical model. In the finite element model, the vibration displacements of the resonator may be distributed in the
x,
y, and
z directions, while that of the resonator can only be generated along the axis of the spring. Nevertheless, the two-dimensional theoretical resonator model proposed in this study can effectively analyze the lower boundary of the bandgap.
The special points are selected at the upper and lower boundary of bandgaps to investigate the mechanism of the bandgap generated by the resonators.
Figure 5 shows the points of the resonator for A1, A2, A3, A4, and A5, respectively. Additionally,
Figure 7 displays the vibration modal diagrams. As illustrated in
Figure 7, the vibrations of surface waves modes are primarily concentrated in the internal vibrators of the resonator. In
Figure 7a, the vibrational displacements generated by the elastic waves are concentrated in the second and third vibrators, and the second vibrator triggers the greatest displacement. Furthermore, the second vibrator produces a horizontal displacement along the
y-axis, while the third vibrator produces a horizontal displacement in the opposite direction. As illustrated in
Figure 7b, the vibrational displacement occurs primarily in the first vibrator, which rotates clockwise along the
y–z plane. The maximum displacement is observed on the four edges of the rubber bar. In
Figure 7c, the horizontal vibrational displacement is concentrated in the second vibrator along the
x-axis, which causes stretching on one side of the rubber bar and compression on the other side. Small displacements are observed in the third vibrator, which are in the opposite direction to the second vibrator. The maximum displacement occurs on the copper ball and on the compressed side of the rubber bar.
Figure 7d illustrates that a rotation is produced in three dimensions in the first and second vibrators. In addition to horizontal displacements in the
x- and
y-axis directions, rotations in the
y–z plane were also observed. The maximum displacement occurs on the edge of the rubber bar. As illustrated in
Figure 7e, the third vibrator produced two rotational displacements in opposite directions in the
y–z plane. The maximum displacement occurs on the four edges of the rubber bar. The energy in the elastic waves is consumed by the vibrations that are generated by the internal vibrators, which can attenuate the elastic waves and protect the building situated behind the resonator.
5. Transmission Spectrum Analysis
In
Section 4, the bandgaps of the resonator are calculated at infinite periods. However, it is not possible to set up an infinite period resonator in practical engineering. Therefore, this section will calculate the transmission spectrum to verify the attenuation effect of the proposed resonator on elastic waves.
The transmission spectrum function is defined as
where
u1 and
u2 are the displacements at the pickup positions without and with barrier, respectively.
Figure 8 shows the finite element model of the resonator periodically arranged in soil. To simulate the attenuation effect of the resonator on elastic waves in a semi-infinite space, a perfectly matched layer (PML) is set at the boundary of the model to absorb the elastic waves, which can reduce the reflection of elastic waves and improve the accuracy of the results. A distance of 16
a is set between the excitation and the resonator to convert the unit line load into surface waves. To determine the attenuation of elastic waves by the soil, the transmission spectrum analysis with no barrier can be calculated. A unit line load is applied along the negative direction of the
z-axis at the excitation position and the displacement is collected at the pickup position. Then, the transmission spectrum diagram is calculated using Equation (14).
Figure 9 shows the transmission spectrum diagram of the 10-cycle resonator within the frequency of the Rayleigh surface waves. It can be seen that the resonator effectively attenuates the elastic waves in the bandgap. The maximum attenuation values reached −19.85 dB and −21.67 dB at 5.14 Hz in the first bandgap, and 10.66 Hz in the third bandgap, respectively.
To further investigate the mechanism of attenuation of elastic waves by the resonator, the total displacements at 1 Hz and 5.1 Hz were analyzed and are compared in
Figure 10. As shown in
Figure 10a,b, it can be seen that the resonator cannot attenuate elastic waves with the frequency of 1 Hz, which means that the propagation path of the elastic waves remains unaltered and the energy carried in the elastic waves is not reduced.
In
Figure 10c,d, the elastic waves for 5.1 Hz are greatly attenuated by the resonator. As shown in
Figure 10c, the elastic waves emitting from the excitation point are gradually converted into body and surface waves. In
Figure 10d, the resonator blocks the surface waves and changes the surface waves into body waves transmitting to the deep soil, resulting in minimal vibration of the soil surface behind the resonator. The vibrations of soil in front of the resonator generated by the elastic waves are intensified. Some of the energy of the elastic waves is converted into body waves, while the remainder is concentrated in the internal vibrators of the resonator.
Figure 10e illustrates the vibration of the internal vibrators of the resonator when the excitation frequency of the surface waves is 5.1 Hz. It can be found that the internal vibrators produce a substantial vibration, which consumes the energy of the elastic waves.
6. Time History Analysis
Time history analysis of artificial and seismic waves is used to further verify the attenuation of the resonator in this section. The Hanning modulation window function is selected for the artificial waves, and is shown as follows:
where
Amax is the maximum amplitude of the artificial waves,
fc is the center frequency,
C is the periodic number, and
t is the duration.
Figure 11 shows the numerical model of the resonator arranged periodically along the
x direction. To simulate the semi-infinite space, a low-reflection boundary condition is set around the model, which is used to reduce the reflection of elastic waves at the boundary. The displacement of the artificial waves or seismic waves will be applied at the excitation position, and the acceleration amplitude will be collected at the pickup position. The center frequency of the artificial waves is selected as 5.1 Hz in the bandgap range and 1 Hz out of the frequency bandgap, respectively, as shown in
Figure 12. The acceleration time history of the El Centro waves is shown in
Figure 13.
Figure 14 shows the acceleration response of the resonator under artificial waves with main frequencies of 1 Hz and 5.1 Hz. It can be seen that the resonator can effectively attenuate the artificial waves when the excitation frequency of 5.1 Hz falls into the frequency bandgap, which means that the acceleration of the artificial waves with the frequency of 5.1 Hz can be reduced by 65%, as shown in
Figure 14b, whereas the resonator almost cannot mitigate the artificial waves when the central frequency of 1 Hz falls outside the frequency bandgap, as shown in
Figure 14a.
Figure 15 illustrates the acceleration time history and the response spectrum of the resonator under the El Centro seismic wave.
Figure 15a demonstrates that the resonator significantly attenuates the El Centro seismic wave during the time period of 0–50 s. The resonator reduces the peak of the seismic waves to 0.39 with an attenuation rate of 61%. In
Figure 15b, the resonator significantly attenuates the amplitude of the El Centro seismic wave in the bandgap with an attenuation rate of 65.75%. These results demonstrate that the resonator can effectively attenuate seismic waves within the frequency bandgap of the resonators.
In
Figure 15a, the maximum acceleration of the El Centro seismic wave occurs at 2.75 s with the resonator, while it is only 0.00958 at 2.75 s with the 10-cycle resonator. In order to further analyze the attenuation of the El Centro seismic wave by the 10-cycle resonator, the total displacements of the seismic waves with and without the resonator are shown in
Figure 16 at 2.75 s. As shown in
Figure 16a, the surface of the soil body at the pickup point produced a large ground displacement without the barrier. The upper parts of the soil body produced many concentrated displacements, with the largest concentrated displacement occurring in the upper part of the soil body at the pickup point. In
Figure 16b, it can be seen that the surface displacement of the soil body behind the resonator is significantly attenuated and the displacement of the soil body below the resonator is also attenuated. Therefore, it further indicates that the resonator not only effectively attenuates the surface waves of seismic waves, but also converts part of the seismic waves into body waves and also proves the conclusions drawn from the previous section.
7. Effect of Resonator Size
The size of the resonator is the key factor affecting the bandgaps. Based on the mechanism of bandgap generation in resonators as well as the previous research results [
7,
10,
21,
23,
35,
36,
37,
48], this study analyzes the effects of the outer frame size and the radius of the internal vibrator of the resonator on the bandgap, respectively. The outer frame size of the resonator is related to the other dimensions of the resonator structure. In order to ensure that the other important parameters of the structure remain unchanged (the size of the internal vibrator and the thickness of the frame), the length of the rubber is increased while increasing the size of the outer frame, and the other parameters of the resonator are shown in
Table 2. The bandgap is calculated in the range of 0–20 Hz for frame sizes
a1 of 1.0 m, 1.1 m, 1.2 m, and 1.3 m. Frequency dispersion curves are shown in
Figure 17. In calculating the effect of the radius of the internal vibrator on the bandgap range, other parameters are kept constant, as shown in
Table 2, and the bandgap is calculated in the range of 0–20 Hz when the radius
r of the internal vibrator is 0.05 m, 0.10 m, 0.15 m, and 0.20 m, respectively, and frequency dispersion curves are shown in
Figure 18.
In
Figure 17, it can be seen that when the outer frame size
a1 is 1 m, the first five bandgaps can cover the frequency range of the Rayleigh wave. As the frame size
a1 increases, the number of bandgap increases to cover the frequency range of the Rayleigh wave, and the number of bandgaps are 6, 7, and 8 when
a1 is 1.1 m, 1.2 m, and 1.3 m, respectively. As the outer frame size
a1 increases from 1 m to 1.3 m, the upper and lower boundaries of the first five bandgaps are decreasing, and the lower boundary of the first bandgap is 3.4249 Hz when
a1 is 1.3 m, which is 1.2546 Hz lower compared to 4.6795 Hz when
a1 is 1 m; the widths of the first, second, fourth, and fifth bandgaps are decreasing, while the third bandgap has a slight increase in width. The larger the outer frame size
a1, the more the lower boundary of the bandgap decreases. When increasing the outer frame size, the length of the rubber connecting the internal vibrator to the outer frame increases accordingly, and the increase in rubber length helps to increase the vibration displacement of the internal vibrator as well as consume the energy of the elastic wave. Therefore, as the size of the outer frame of the resonator and the length of the rubber increase, the lower boundary of the bandgap generated by the resonator gradually decreases, which contributes to the attenuation of the low-frequency elastic waves and enables the generation of more bandgaps.
In
Figure 18, it can be concluded that the lower boundary of the bandgap decreases as the radius of the internal vibrators increase. The lower boundary of the bandgap is 7.0796 Hz when the radius of the internal vibrators
r = 0.05 m and decreases to 4.6795 Hz when the radius of the internal vibrators
r = 0.20 m. From the previous analysis of the vibration diagrams of the resonator, it can be seen that the bandgap of the resonator is produced due to the absorption of the energy of the elastic wave by the internal vibrators. The larger the radius of the vibrators, the more elastic wave energy is absorbed by the internal vibrators, so the lower the bandgap that can be generated.
Based on the analysis of the parameters of the resonator, it is clear that the larger the resonator, the lower the frequency bandgap, and the larger the radius of the internal vibrators, the lower the frequency bandgap.
8. Effect of Damping
As mentioned above, soil properties are consistent with Biot’ s theory and undamped in the previous sections. While the actual soil has damping characteristics, the influence of the soil damping is considered in this section. The Rayleigh damping is used in this study, and the equation for Rayleigh damping is written as follows:
where
is mass damping coefficient, and
is stiffness damping coefficient. Rayleigh damping is affected by the frequency range and damping ratio of the structure. In this section, the frequency range of the structure is set as 0.1–20 Hz and can be written as
where
f1 and
f2 are the lower and upper boundaries of the structural frequency, respectively. The angular frequencies can be calculated as follows:
where
and
are the lower and upper boundaries of the angular frequency for structure, respectively. Rayleigh damping coefficients are written as
where
is damping ratio.
The soil damping ratio is from 0.01 to 0.08, with a general value of 0.05 taken in this study. The mass and stiffness damping coefficients of Rayleigh damping are calculated to be 0.0625 s
−1 and 0.0008 s, respectively. The time history analysis of the El Centro seismic wave with soil damping is calculated and the results are shown in
Figure 19. It can be found in
Figure 19a that the soil damping significantly attenuates the El Centro seismic wave with a maximum acceleration of 0.2286, and the attenuation efficiency of the resonator considering soil damping is improved by 32.67% compared to the case of the resonator without soil damping.
Figure 19b illustrates that soil damping attenuates frequencies within 1.96–20 Hz, with a maximum amplitude of 0.477. However, soil damping does not attenuate frequencies from 0 Hz to 1.96 Hz.
Simultaneously, the material damping of the rubber and copper is also considered.
Figure 20 shows the acceleration responses of the resonator considering the soil and material damping. It can be found that the attenuation by the resonator considering both soil and material damping is almost the same as when only considering soil damping, which indicates that material damping has a negligible effect on the attenuation performance of the resonator.
9. Conclusions
In this study, a low-frequency bandgap resonator with three internal vibrators was investigated using theoretical and finite element methods, respectively. Firstly, a resonator with three internal vibrators was proposed based on the mass spring theory. Secondly, the attention effect of the resonator was studied through transmission spectrum and time history analysis. Finally, the effects of soil and material damping of the resonator on the mitigation of the Raleigh waves were discussed and compared. The main conclusions were obtained as follows.
(1) A theoretical model of the resonator was proposed based on the local resonance theory, and the mass-spring theory was used to establish the dynamic equilibrium equation. The results showed that the resonator can generate a wide low-frequency bandgap.
(2) The resonator has omnidirectional bandgaps with a first five bandgaps (4.68–22.01 Hz) covering the frequency of Rayleigh surface waves. It was indicated that the mechanism of bandgap of the resonator is a local resonance phenomenon through vibration modal analysis.
(3) It was demonstrated from transmission spectrum of the 10-cycle resonator that the resonator can effectively attenuate elastic waves within the bandgap. The reason is that the resonator converts surface waves into body waves and the internal vibrators consume the energy of elastic waves.
(4) The time history analysis results show that the resonator can decrease the accelerations of artificial waves or seismic waves by 61% when the main frequency is in the bandgap of the resonator. The amplitude of the El Centro seismic wave in the bandgap was attenuated by 65.75%.
(5) The effects of soil and material damping of the resonator on the attenuation of seismic waves were, respectively, considered and calculated. The results showed that the attenuation rate of the acceleration was 32.67% when the numerical model considered the soil damping, while the material damping of the resonator had little effect in attenuating seismic waves.
(6) In this study, the resonators could attenuate seismic surface waves in the main frequency range of 4.68 Hz to 22.01 Hz and it was found that a 10-cycle resonator could effectively attenuate El Centro seismic waves. Moreover, the higher the number of resonators, the better the attenuation effect. Therefore, the number of resonators can be flexibly arranged according to the earthquake magnitude, building strength, and site conditions in practical engineering applications.
(7) In this study, a resonator with some bandgaps from 4.68 Hz to 22.01 Hz was proposed, but the study is not perfect and there are still many deficiencies, which will be further improved. Firstly, a structure that is not only small in size but also has a lower frequency bandgap to attenuate seismic waves with ultra-low frequency as the main frequency can be developed. Then, scaled or full-scale tests can be conducted to further demonstrate the actual attenuation of elastic waves by the resonator structure. Finally, specific protected building structures will be added to the finite element model calculations to analyze the vibration displacements of the building with and without resonator protection and to further analyze the protective effect of the resonator on the building.
(8) The computational method of machine learning can be applied to the calculation of seismic surface wave attenuation by resonators, such as analyzing the dispersion curves of the resonators, estimating the amplitude of vibrations generated by the resonators and the ground under different conditions, optimizing the parameter design of the resonators, and calculating the theoretical analysis of the resonators in terms of elastic dynamics. These machine learning aspects will be studied and analyzed later.