Next Article in Journal
FEM Simulation of the Riveting Process and Structural Analysis of Low-Carbon Steel Tubular Rivets Fracture
Previous Article in Journal
Development of VUMAT and VUHARD Subroutines for Simulating the Dynamic Mechanical Properties of Additively Manufactured Parts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Low-Frequency Bandgaps of the Lightweight Single-Phase Acoustic Metamaterials with Locally Resonant Archimedean Spirals

1
Department of Mechanics and Tianjin Key Laboratory of Nonlinear Dynamics and Control, Tianjin University, Tianjin 300350, China
2
Key Laboratory of Aeroacoustics and Dynamics, Aircraft Strength Research Institute, Xi’an 710065, China
3
School of Civil Engineering, Henan Polytechnic University, Jiaozuo 454003, China
*
Authors to whom correspondence should be addressed.
Materials 2022, 15(1), 373; https://doi.org/10.3390/ma15010373
Submission received: 1 December 2021 / Revised: 21 December 2021 / Accepted: 31 December 2021 / Published: 5 January 2022

Abstract

:
In order to achieve the dual needs of single-phase vibration reduction and lightweight, a square honeycomb acoustic metamaterials with local resonant Archimedean spirals (SHAMLRAS) is proposed. The independent geometry parameters of SHAMLRAS structures are acquired by changing the spiral control equation. The mechanism of low-frequency bandgap generation and the directional attenuation mechanism of in-plane elastic waves are both explored through mode shapes, dispersion surfaces, and group velocities. Meanwhile, the effect of the spiral arrangement and the adjustment of the equation parameters on the width and position of the low-frequency bandgap are discussed separately. In addition, a rational period design of the SHAMLRAS plate structure is used to analyze the filtering performance with transmission loss experiments and numerical simulations. The results show that the design of acoustic metamaterials with multiple Archimedean spirals has good local resonance properties, and forms multiple low-frequency bandgaps below 500 Hz by reasonable parameter control. The spectrograms calculated from the excitation and response data of acceleration sensors are found to be in good agreement with the band structure. The work provides effective design ideas and a low-cost solution for low-frequency noise and vibration control in the aeronautic and astronautic industries.

1. Introduction

Lightweight periodic honeycomb materials (LPHM) are a typical kind of acoustic metamaterials [1,2,3,4,5,6]. With the rational design of sub-wavelength units to control elastic wave propagation, these structures can be widely used in vibration isolation [7,8], noise insulators/absorbers [9,10], negative refraction [11,12], and acoustic cloaking [13,14]. As the Bragg scattering bandgap requires the lattice size to be comparable to the corresponding wavelength [15,16,17], the reduction of low-frequency vibrations and noise below 500 Hz remains a challenging assignment for LPHM design [18,19,20] when used in engineering applications in aeronautics [21] and railways [22].
Introducing a local resonant mass block is currently one of the popular ways to achieve a low-frequency bandgap below 500 Hz for LPHM. Liu et al. [23] used a simple cubic lattice form to constitute a local resonance unit using a high-density mass block wrapped in soft rubber material. The artificial periodic structure is formed by periodically arranged local resonance units in an elastic medium, which successfully exploits the local resonance effect of elastic waves to achieve a bandgap of around 400 Hz in a 20 mm cubic lattice. In contrast, traditional local resonance structures typically allow for a very narrow resonance bandgap at low frequencies [24]. For the purpose of obtaining a wide resonant bandgap, Dong et al. [25] proposed a new multiple bandgap phononic crystals capable of generating lots of flat bandgaps in the low-frequency band, however, the flat bandgaps generated by this unique local resonant phononic crystal are still all above 500 Hz. To address the above problem, Li et al. [26] proposed a phonon crystal structure in the form of a cylinder with a periodic arrangement attached to a thin plate. This cylinder structure attached to a thin plate corresponds to the formation of a “solid-state Helmholtz resonator” and succeeds in generating a large number of broadband gaps in the range 256 Hz–855 Hz. In addition, Ning et al. [27] similarly designed a tunable metamaterial consisting of a frame structure, airbag, and a counterweight by introducing a local resonant mass block. They modulated the designed acoustic metamaterials by the gauge pressure and gas temperature inside the airbag, effectively attenuating wave propagation in the 13 Hz to 90 Hz range. Although the low-frequency bandgap of LPHM below 500 Hz can be obtained by introducing local resonance mass blocks in the above-mentioned references, the lightweight efficiency of LPHM has severely deteriorated. On the other hand, the discontinuous distribution of materials with different properties in the complex spatial structure can bring great challenges to the fabrication [28,29].
To endow the LPHM with the simultaneous advantages of lightweight efficiency and low-frequency bandgaps, the single-phase lightweight periodic honeycomb materials with subwavelength bandgaps have attracted much attention in the research community [18,30,31,32,33,34]. Chen et al. [35,36] designed a star-assisted metamaterial with a low-frequency bandgap and double-negative characteristics in a certain frequency range using single-phase materials, which solved the problem that conventional double-negative acoustic metamaterials are difficult to be applied due to their complex structure and multi-phase material composition. The novel lightweight bidirectional re-entrant lattice metamaterials were proposed by Ren et al. [37], which forms a wide bandgap of about 2 kHz in the range of 2.7 kHz to 4.7 kHz. The single-phase acoustic metamaterials with periodically arranged diverging star-shaped cells resulting in a “low frequency” bandgap from 1.44 kHz to 1.56 kHz were proposed by Kumar and Pal [38]. In the above-mentioned, single-phase LPHM, although the simultaneous advantages of lightweight efficiency and low-frequency bandgaps are realized, the low-frequency bandgaps are all larger than 500 Hz.
Aiming at the challenging problem of vibration and noise reduction below 500 Hz for the single-phase LPHM, a new kind of local resonance single-phase lightweight periodic honeycomb materials, the square honeycomb acoustic metamaterials with locally resonant Archimedean spirals (SHAMLRAS) are proposed in this paper. The SHAMLRAS consists of Archimedean spirals of the same material combined in a square honeycomb structure. These acoustic metamaterials with the introduction of Archimedean spirals have special resonance characteristics compared to existing single-phase acoustic metamaterials and have the capacity to form multiple bandgaps below 500 Hz (down to approximately 184.5 Hz) within a lattice size of 40 mm. The research has the potential to provide an invaluable guide to the engineering suppression of low-frequency noise and vibration.
The paper is organized into six sections including the introduction above. The second part describes the geometry of SHAMLRAS and illustrates the tools for wave propagation studies. The low-frequency bandgap characteristics and directional attenuation properties of elastic waves in the SHAMLRAS periodic structure are thoroughly discussed in Section 3. For more flexibility in obtaining the low-frequency bandgap width we need, the effect of the material parameters, the arrangement of spirals, and the variables of parameter control equations on the bandgap width and position being investigated in Section 4. The transmission loss experiments and COMSOL verification of the low-frequency bandgap are demonstrated in Section 5. Finally, the major achievements of this work are presented in Section 6.

2. Mechanical Design and FE Modeling of the SHAMLRAS

2.1. Mechanical Design

In this study, a low-frequency single-phase acoustic metamaterials consisting of multiple spirals is proposed, based on a square honeycomb structure and a circular array of four spirals. The structural composition and lattice arrangement of acoustic metamaterials are shown in Figure 1a–c, with a lattice constant of L = 40 mm and a ligament thickness of p = 1 mm. We can effectively reduce the bandgap frequency by tuning the parameters of the spiral line equation.
The spiral line control equation is derived from the Archimedean spiral line and is derived as follows:
r = α + β θ
where θ is the pole angle, α is the pole diameter when θ = 0, i.e., the initial radius, and β is the increase (or decrease) of the pole diameter per 1 rad of rotation.
Suppose R2 is the inner diameter, n is the turn number, and d is the circle distance, then this polar equation can be written in the following form:
{ r = ( R 2 + n d t ) θ = n 2 π t ,   t [ 0 , 1 ]
Therefore, the corresponding Cartesian coordinate equation for the Archimedean spiral plane is shown below, which is the parametric control equation for the spiral in this paper.
{ x ( t ) = ( R 2 + n d t ) cos ( n 2 π t ) y ( t ) = ( R 2 + n d t ) sin ( n 2 π t ) ,   t [ 0 , 1 ]

2.2. FE Modeling for the Free Wave Propagation

The unit cell and Brillouin zone of periodic lattices are shown in Figure 1d. In the unit, e i ( i = 1 , 2 ) is the basic lattice vector, which can be expressed by orthogonal Cartesian basic vector and lattice constant as:
e 1 = L i e 2 = L j
A reciprocal lattice space is defined according to the direct lattice space. In general terms, the relationship between the basis vectors of the direct and reciprocal lattices is satisfied as follows:
e i e j * = 2 π δ i j
where e i denote the basis vectors of the direct lattice, e j * denote the basis of reciprocal lattice, and the subscripts i and j take the integer value 1 and 2. The δ i j is the Kronecker delta function, and the expressions are as follows:
δ i j = { 1 , i = j 0 , i j
The coordinate positions of the lattice points in the reciprocal lattice can be represented by the reciprocal lattice vector G , which is a linear combination of the reciprocal lattice basis vectors:
G = r o + n 1 e 1 + n 2 e 2
where n 1 and n 2 are integers, and r o is the displacement of Point O.
Additionally, the reciprocal lattice vectors of the SHAMLRAS can be expressed as:
e 1 * = ( 2 π L , 0 ) e 2 * = ( 0 , 2 π L )
Thus, the entire SHAMLRAS periodic structure can be constructed by shifting the unit cell along the basic lattice vector ( e 1 , e 2 ) in a two-dimensional space. The Brillouin zone of the basic lattice can also be obtained, as shown in Figure 1d, in which the black shadow region is the irreducible Brillouin zone (IBZ). The coordinates of the boundary points of the IBZ are shown in Table 1.
According to Bloch’s theorem, the part u ( r ) of the eigenwave field volume associated with the spatial position r can be expressed in the form of a spatial plane wave
u ( r ) = U k ( r ) e i k · r
where r = ( x , y , z ) is the position vector, k = ( k x , k y , k z ) is the wave vector in the first Brillouin zone, i is the imaginary unit, and U k ( r ) is eigenwave amplitude.
The boundary displacement of the periodic structure can be controlled by Bloch boundary conditions (Floquet periodicity), so there is
u ( r + R ) = e i k · r u ( r )
where r is the position vector of the node on the boundary, and R is the basis vector of the lattice.
Determine the form function according to the general finite element procedure and establish the stiffness and mass matrices within the unit to obtain the generalized eigenequations of the unit:
( k s ω 2 M s ) u ( ν ) = f
where k s is the unit stiffness matrix, M s is the mass matrix, u ( ν ) and f are the vectors of generalized nodal displacements and forces.
The SHAMLRAS periodic structure can be reduced to a series of cells using Equation (10) for the whole calculation. Two edges in x- and y-directions are selected as source boundaries in COMSOL Multiphysics 5.5, and the two edges selected for the destination boundary correspond to the two edges of the source boundary. The Floquet periodicity conditions at the corresponding boundaries of the periodicity cell are expressed as:
u destination = e i k f ( r destination r source ) u source
where u is a vector of dependent variables, and the vector k f represents the spatial periodicity of the excitation.
The wave vector k is parametrically swept along the path O A B O according to the first Brillouin zone of the SHAMLRAS structure given in Figure 1d. With Equations (11) and (12), the structural eigenfrequencies can be found for a given k . Thus, the dispersion curve is finally plotted with the wave vector k as the horizontal coordinate and the eigenfrequency as the vertical coordinate to obtain the energy band diagram of the SHAMLRAS, and the forbidden band between the dispersion curves is the band structure. Moreover, the eigenfrequencies on the dispersion branches all correspond to the mode shape the structure has. The complete surface ω = ω ( k ) is the dispersion surface mentioned in Section 3.2. In this case, the corresponding dispersion surfaces of the different dispersion branches represent a class of modes at the corresponding eigenfrequencies. Hence, the n-th ( n = 1 , 2 , 3 , ) dispersion surface can also be called the n-th ( n = 1 , 2 , 3 , ) mode. In this case, there are as many orders of eigenfrequency as there are dispersion surfaces.

3. Bandgap and Directional Propagation Characteristics of Elastic Waves for the SHAMLRAS

3.1. Bandgap Characteristics

In this section, the in-plane dynamical properties of SHAMLRAS based on Figure 1 are investigated. The geometrical parameters of the SHMLRAS and the parameters of the photosensitive resin used are shown in Table 2.
Based on the material and geometrical parameters provided, the calculated bandgap structure of the SHMLRAS in the low-frequency range is shown in Figure 2a, where the pink shaded area is the bandgap in the calculated frequency range. For ease of discussion and analysis, a detailed schematic of the first bandgap is shown in Figure 2b. Clearly, there are a total of five complete bandgaps in the 0–500 Hz range. These five bandgaps are located between the third and fourth dispersion curves, between the sixth and seventh dispersion curves, between the eighth and ninth dispersion curves, between the tenth and eleventh dispersion curves, and between the twelfth and thirteenth dispersion curves respectively. The corresponding frequency ranges are [184.42, 184.69] Hz, [237.18, 249.80] Hz, [252.98, 272.13] Hz, [339.70, 368.00] Hz, and [368.60, 425.75] Hz respectively. For single-phase cellular metamaterials, this bandgap result is very promising and shows that lightweight SHMLRAS has excellent elastic wave attenuation properties in the low-frequency region as well as favorable prospects for engineering applications.
The application of point group theory, which is used extensively in physics and molecular chemistry, to analyze and predict the electromagnetic response of complex structures is a highly rewarding endeavor [39]. However, the mode analysis approach used in these works [38,40] more easily helps us to understand the mechanism of bandgap generation. The vibrational modes at Points O, A, and B are analyzed below for the upper (green circles) and lower (red circles) dispersion curves corresponding to the former three bandgaps in Figure 2a. The size and direction of the arrows in the mode shape diagrams represent the displacement distance and direction of the structure relative to its original position. Figure 3 shows the mode shapes at Points O, A, and B for the upper and lower boundaries of the first bandgap respectively. It can be seen that the displacement and direction of the mode shapes at the upper and lower boundaries of the fourth (1st upper) and third (1st lower) dispersion curves corresponding to the Point O are very similar in the second and fourth quadrants. In the first and third quadrants, however, the directions of the modes are opposite. At Point A, we can notice a clear exchange of energy in the structural vibrations, with the spiral vibrations in the second and fourth quadrants becoming the spiral vibrations in the first and third quadrants. In the A B direction, attenuation occurs in the vibration modes corresponding to the fourth and third dispersion curves. Moreover, the 1st up vibration in the second and fourth quadrants decays significantly, which leads to the spiral vibration in the first three quadrants. Comparing Figure 2b with the previous discussion, we can see that the energy exchange occurring at Point A is the main reason for the creation of the first bandgap.
Figure 4a–f shows the mode shape diagrams corresponding to the sixth (2nd lower) and ninth (3rd upper) dispersion branches at the Points O, A, and B. The marked arrows show that there is a clear difference in the direction of the spiral vibrations in the different quadrants, i.e., a greater degree of local resonance. For example, at Point B, the sixth dispersion branch tends to vibrate downwards in the first quadrant, leftwards in the second quadrant, upwards in the third quadrant, and to the right in the fourth quadrant. The ninth dispersion branch tends to vibrate in the lower left, lower right, upper right, and upper left in the four quadrants respectively. The effect of the elastic wave on the structure is weakened or completely counteracted by the effect of the local resonance spirals on the structure, i.e., there is a positive correlation between the local resonance stratification of the spiral and the attenuation performance of the elastic wave.
The local resonant mode of the unit structure is excited as the frequency approaches the natural frequency of the resonant. At this point, the elastic waves in the structure will be strongly coupled with the structural local resonant mode. Furthermore, the energy is localized due to the constant exchange into the resonant unit and the elastic wave will not propagate further. In the band structure, this is manifested by the truncation of the energy band starting at Point O by the resonant straight band, resulting in the formation of the third and fourth bandgaps. The vibration of the spirals on the seventh (2nd upper) dispersion branch is essentially constant in the third quadrant, as shown in Figure 4g–i. In the fourth quadrant, however, only the direction of vibration does not change. Meanwhile, the direction of vibration of the spiral in the second quadrant shows a similar attenuation of amplitude after reversal as in the third and fourth quadrants. On the O A B direction, the eighth (3rd lower) dispersion branch shows marked attenuation of the spiral amplitude in the first three quadrants. The spiral vibration only shows a reversal of direction in the second quadrant, except in the fourth quadrant which remains essentially unchanged.
Comparing Figure 3 and Figure 4, the consumption of elastic wave energy by the structure in the calculated frequency range originates from the local resonance of the spiral. Furthermore, the degree of local resonance of the spiral increases, the better performance of designed structure for elastic wave attenuation.

3.2. Dispersion Surfaces

In this paper, the dispersion surfaces of the SHAMLRAS periodic structure are calculated using the finite element method in the solid mechanics interface of COMSOL. The four edges of the square honeycomb metamaterials cell are divided by taking mesh = 40, i.e., the four edges of the first Brillouin zone in Figure 1d are divided into 1600 points. Once we have swept the wave vector k for all the points inside and on the boundary of the first Brillouin zone, the complete eigenfrequency surface ω = ω ( k ) corresponding to the different dispersion branches can be calculated using COMSOL. In this case, a complete eigenfrequency surface ω = ω ( k ) is also a mode. The corresponding dispersion surfaces of the different dispersion branches represent a class of modes at the corresponding eigenfrequencies. Hence, the n-th ( n = 1 , 2 , 3 , ) dispersion surface can also be called the n-th ( n = 1 , 2 , 3 , ) mode. Here we select modes near the first three bandgaps (third to ninth dispersion branches) for further discussion and analysis.
The dispersion surface method can be used to show the elastic wave propagation properties in the periodic structure in a three-dimensional and visual way. Figure 5a shows the dispersion relations of the seven dispersion surfaces of the third to ninth dispersion branches of the SHAMLRAS periodic structure in three-dimensional space, where the x- and y-axis coordinates represent the coordinate of the first Brillouin zone point and the z-axis represents the frequency f. The comparative relationship between the SHAMLRAS dispersion relation and the bandgap structure in the horizontal view is shown in Figure 5b. There is a good correspondence between the dispersion curves of the SHAMLRAS structure and the dispersion curves of the bandgap structure, as well as between the bandgap positions of the two diagrams. Simultaneously, it proves the correctness of the dispersion relation calculated by the above method.

3.3. Directional Propagation Property of Elastic Waves

We can obtain the ISO-frequency contours corresponding to the modes by projecting the dispersion surfaces in the k x and k y plane in Figure 5a and extracting the eigenfrequencies according to different division intervals.
The group velocity of two-dimensional periodic structures along the x- and y-directions for a given frequency can be written as:
c g x = a x ω k x ,   c g y = a y ω k y
where a x and a y denote the lattice constants in the x- and y-directions respectively.
The group velocity is defined as the gradient of an ISO-frequency curve when the dispersion surface is described by a two-dimensional contour. At the same time, the direction of the outer normal at each point on the contour is the direction of the group velocity at that point, which represents the direction of energy propagation of the vibration at that frequency. Therefore, the group velocity can be used to express the speed of energy propagation as well as the direction and magnitude of the elastic wave propagation. With the calculation of the gradient at each point on the contour, we can obtain the direction of propagation and the region of propagation for a given frequency of the vibration. The values of c g x and c g y are used as the x- and y-coordinates to obtain the group velocity for a specific frequency.
With increasing frequencies, the SHAMLRAS structure vibration propagation is mostly concentrated in the k x and k y directions, as shown in Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12. At the same time, the vibration propagation generally behaves as follows: it is enhanced, then weakened, and then enhanced again. In particular, the sixth mode case has the weakest vibration propagation in the k x and k y directions. On the other hand, we can observe from Figure 6a, Figure 7a, Figure 8a, Figure 9a, Figure 10a, Figure 11a and Figure 12a that there is weak anisotropy in the k x and k y directions for the fifth and sixth modes at the lower frequencies of the dispersive surface, whereas there are strong anisotropy in this direction for the fourth, seventh and eighth modes. Furthermore, only the fifth and ninth modes have strong anisotropy in the direction with the diagonal as the axis of symmetry.
To be more precise in evaluating the propagation of elastic waves on different dispersion surfaces, the ISO-frequency contours for the three frequency cases are marked in (b) of Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 with heavy black solid lines, and the group velocities corresponding to these ISO-frequency contours are shown in (c,d) of Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 respectively. The group velocity in Figure 6c–e shows a distribution of points in the 0° to 360° direction. The points on the outside of the group velocity are more densely distributed in the k x and k y directions as the frequency increases, however, the group velocity points in this direction disappear with the frequency increases to f = 184.2 Hz. In contrast, with increasing frequency f the group velocity distribution is more densely distributed in the region near the diagonal as the axis of symmetry.
The clusters of velocities in the k x and k y directions increase remarkably with increasing frequency, and the points in the other directions decrease significantly, according to (c,d) of Figure 9, Figure 10 and Figure 12. However, a highly visible concentration of group velocity points in the diagonal direction of the graph is observed in Figure 6c,d as well as in Figure 11c,d when the frequencies are located in the fourth and eighth modes, and the degree of the clusters increases noticeably with increasing frequency. That demonstrates that there is a strong energy aggregation at this location, and it also indicates that the propagation of energy in the k x and k y direction is much weaker in this direction compared to the diagonal direction due to the small group velocity distribution. Consequently, we can consider that elastic waves form vibration blind zones in the k x and k y directions, and this characteristic is very useful for vibration isolation needs in specific directions of engineering.

4. Influence of the SHAMLRAS Bandgap

4.1. Influence of Spiral Geometry on the Bandgap

From the analysis in Section 2.1, the structure of SHAMLRAS changes markedly when the inner diameter (R2), turn number (n), circle distance (d), the tart parameters (t0), and the End parameters (t1) of parameter t are changed in Equation (5). The changes in these parameters are further influenced by the bandgap characteristics of the SHAMLRAS structure and the dispersion relations. To be more flexible in dealing with different frequencies of elastic waves in engineering damping applications, the influence of the spiral geometry parameters on the bandgap width and position is discussed below.
The influence of the circle distance d and the inner diameter R2 on the bandgap width of the SHAMLRAS structure are shown in Figure 13a–d respectively. With increasing circle distance d in the range of 1000 Hz, the second bandgap width is found to show an increase followed by a decrease. Although the bandgap width of the SHAMLRAS structure decreases as the circle distance d increases, the frequency of each bandgap also decreases remarkably. Furthermore, a comparison of Figure 13a,b shows that the first bandgap is generated between the third and fourth dispersion curves when d = 2.75, 3.25, 3.75, 4.25. The second to fifth bandgaps then corresponds to between the fifth and sixth dispersion curves, between the seventh and eighth dispersion curves, between the ninth and tenth dispersion curves, and between the eleventh and twelfth dispersion curves respectively. However, in the case of d = 2.25, no bandgap can be observed to be created between the third and fourth dispersion curves, with the first bandgap located between the sixth and seventh dispersion curves. As the inner diameter R2 is varied, the width of the first three bandgaps is in Figure 13c is negatively correlated with the change in R2, while only the width of the fifth bandgap is positively correlated with the change in R2, with only the fourth bandgap showing an increase followed by a decrease. Moreover, the position of the dispersion curves for the first five bandgaps of the SHAMLRAS structure has not changed by comparing Figure 13a,b. From the above discussion, there is a clear indication that we should increase the value of the circle distance d as well as R2 in the SHAMLRAS structure within a limited dimension when we need a low-frequency bandgap.
The dependence of the bandgap width on turn number n, starting value t0, and end value t1 is shown in Figure 14a–f respectively. All the shapes of the SHAMLRAS structures are noted to change more markedly with the variation of parameters. The variation of these parameters also leads to changes in the position of the bandgap of the SHAMLRAS structure between the different dispersion curves, especially in the low-frequency bandgap region. There is a narrower bandgap that is consistent between the second and third dispersion curves as the values of these three parameters gradually approach the structural parameters selected in Table 2. As the selected parameter change to the structural parameters shown in Table 2, the bandgap tends to decrease at all frequency positions, in particular, the frequency values between the second and third dispersion curves overlap, and the bandgap at this position disappears. We also observe an interesting phenomenon: there is a striking similarity in the effect of the turn number n in Figure 14a,b and the end value t1 in Figure 14e,f on the SHAMLRAS geometry and band structure. This suggests a strong one-to-one correspondence between the turn number n and the end value t1 for the SHAMLRAS structure. On the other hand, variations of the three parameters in Figure 14 have a relatively high effect on the bandgap width of the SHAMLRAS structure. Although the variation in bandgap width at this point does not have the same continuity of variation as the bandgap width shown in Figure 13, the position of the dispersion curve is the same for each bandgap, except for the structural parameters in Table 2.
In the above case, the dispersion curves of the second to fifth bandgaps of the different SHAMLRAS structures are located between the third and fourth dispersion curves, between the sixth and seventh dispersion curves, between the eighth and ninth dispersion curves, and between the tenth and eleventh dispersion curves respectively. This indicates that the relationship between the position of the bandgap and dispersion curves is relatively stable for the range of parameters chosen.

4.2. Influence of Spiral Arrangement on the Bandgap

The effect on bandgap width and position of spirals in a square honeycomb structure is explored here. As shown in Figure 15, the bandgap width distribution is calculated for each case by varying the arrangement of the spirals in different quadrants. When the spiral arrangement is varied in accordance with Figure 15, there are decreases to disappear of the bandgap between the third and fourth dispersion curves, the bandgap between the sixth and seventh dispersion curves, and the bandgap between the eighth and ninth dispersion curves when the spiral arrangement is transformed in the first and third quadrants. The different arrangement of spirals demonstrates that only the first three bandgaps of structural parameters shown in Table 2 are more significantly affected.

4.3. Influence of Material Parameters on the Bandgap

The influence of the material parameters Young’s modulus E, Poisson’s ratio ν, and density ρ (only one of these parameters is changed at a time) on the width and position of the bandgap in the numerical simulation are illustrated in Figure 16. The frequency position of each bandgap rises gradually with increasing Young’s modulus E and Poisson’s ratio ν as shown in Figure 16a,c. Simultaneously, the width of the five bandgaps shown in the diagram increases. In contrast, both the width and position of the bandgap in Figure 16e show a negative correlation with the variation in density ρ. Additionally, it can be observed that higher frequency bandgaps are more sensitive to the change of material parameters. In other words, the higher the bandgap frequency, the faster the bandgap width increases (or decreases) with the change of material parameters. On the other hand, the relative positions of each bandgap have not changed in the dispersion curves. It is demonstrated that changes in Young’s modulus E, Poisson’s ratio ν, and density ρ during the simulations do not open a new bandgap (or disappear the existing bandgap) in the band structure of SHAMLRAS. As the analysis above shows, the smaller the values of Young’s modulus E and Poisson’s ratio ρ of the material, the lower the frequency of the bandgap. At the same time, it is easier to obtain bandgaps at lower frequencies with higher values of the density ρ.

5. Filtering Properties of the Finite Periodic Structure of the SHAMLRAS

The dispersion relations and free wave propagation properties in the case of SHAMLRAS infinity structures have been analyzed in the previous sections, but they may not be adaptable to the demands of flexible load-bearing in engineering. In this section, the performance of SHAMLRAS for vibration isolation finite structures is investigated from both experimental tests and numerical simulations. The geometry and loading environment in Figure 17a were used for both the transmission loss experiments and the COMSOL simulations. To evaluate the vibration isolation performance of the SHAMLRAS structure, the frequency response function (FRF) was calculated with the left-hand panels of 5 mm thickness as the excitation input and the right-hand panels of 2 mm thickness as the response receiver. The geometric and material parameters shown in Table 2 were used as the basis for obtaining the specimens required in the experiments through 3D printing, and the experiment was carried out based on the transmission loss experimental setup shown in Figure 17b. After the printed specimen has been fixed to the vertical vibration table through the 5 mm thick prefabricated base, the excitation and response data from the experimental test are calculated using an accelerometer to obtain the spectrum in Figure 17c. We used acceleration as the input signal for excitation in both experimental and simulations, and the calculated spectra were placed in Figure 17c for comparative analysis.
Five decay peaks are corresponding to the band structure that we observed both in experimental and in COMSOL simulations, and their frequencies and ranges match well with the calculated bandgaps in Figure 17c. At the same time, the spectrogram shows a significant amplification at approximately 120 Hz. This amplification phenomenon is attributed to the resonance in this frequency region of the SHAMLRAS periodic structure used in the experiments and simulations. It is shown that as the sweep frequency increases, the vibration frequency of the vertical vibrating table gradually coincides with the eigenfrequency of the SHAMLRAS periodic structure. At this point, the structure deformed significantly and is represented by a resonance peak in the spectrogram.
In addition, the behavior of the elastic wave propagation is shown in Figure 17d for excitation frequencies of f = 220 Hz (outside the bandgap range) as well as for frequencies f = 240 Hz (inside the bandgap range). The elastic wave is significantly suppressed when the excitation frequency f = 240 Hz is in the bandgap range, which corresponds well to the bandgap shown in Figure 2a.

6. Conclusions

In this paper, lightweight single-phase acoustic metamaterials with low-frequency bandgap properties are designed by combining a square honeycomb structure with multiple Archimedean spirals.
The existent Archimedean spirals are shown to open multiple complete bandgaps below 500 Hz with Bloch’s theorem and finite element analysis. The vibrational modes are discussed for the dispersion curves near the first three bandgaps at the boundary points of the IBZ. The generation of low-frequency bandgaps was found to be related to the degree of local resonance of spirals. Characteristics of the directional propagation of elastic waves in SHAMLRAS periodic structures and the attenuation properties of in-plane vibrations are analyzed using dispersion surfaces and group velocities at specific frequencies. Optimizing the width and position of the low-frequency bandgap to enhance wave attenuation can be achieved by adjusting the material parameters, the spiral arrangement, and the parameters of the control equations. In this case, the spirals in circular arrays, the increase of the circle distance d and the inner diameter R2 within a limited dimension are helpful in obtaining a lower frequency bandgap. In the end, the spectrum obtained through transmission loss experiments and COMSOL simulations are used to demonstrate and verify the vibration isolation performance of SHAMLRAS with limited dimensions. That proves the great potential of the SHAMLRAS structure in achieving low-frequency noise and vibration control using single-phase materials.

Author Contributions

Conceptualization, H.G., Q.Y. and X.L.; formal analysis, H.G., Q.Y., X.L. and L.W.; investigation, H.G., Y.Z. and Y.S.; methodology, H.G., X.L., Y.Z., Y.S., Q.D., L.W., J.X. and H.Y.; resources, Y.S., Q.D. and J.X.; writing—original draft, H.G. and X.L.; writing—review & editing, Q.Y., Y.Z., Y.S., Q.D., X.L. and H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The project was supported by the National Natural Science Foundation of China (Grant No. 12072222, 12132010, 12021002, 11991032), the State Key Laboratory of Mechanical Behavior and System Safety of Traffic Engineering Structures (Grant No. SKLTESKF1901), the Aeronautical Science Foundation of China (Grant No. ASFC-201915048001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors appreciate Yingli Li and Gengwang Yan of Central South University for the helpful and insightful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, L.; Liu, H.-T. 3D compression–torsion cubic mechanical metamaterial with double inclined rods. Extrem. Mech. Lett. 2020, 37, 100706. [Google Scholar] [CrossRef]
  2. Fok, L.; Ambati, M.; Zhang, X. Acoustic Metamaterials. MRS Bull. 2008, 33, 931–934. [Google Scholar] [CrossRef]
  3. Cummer, S.A.; Christensen, J.; Alù, A. Controlling sound with acoustic metamaterials. Nat. Rev. Mater. 2016, 1, 16001. [Google Scholar] [CrossRef] [Green Version]
  4. Ye, M.; Gao, L.; Li, H. A design framework for gradually stiffer mechanical metamaterial induced by negative Poisson’s ratio property. Mater. Des. 2020, 192, 108751. [Google Scholar] [CrossRef]
  5. Sun, Y.; Pugno, N.M. In plane stiffness of multifunctional hierarchical honeycombs with negative Poisson’s ratio sub-structures. Compos. Struct. 2013, 106, 681–689. [Google Scholar] [CrossRef] [Green Version]
  6. Sun, Y.; Wang, B.; Pugno, N.; Wang, B.; Ding, Q. In-plane stiffness of the anisotropic multifunctional hierarchical honeycombs. Compos. Struct. 2015, 131, 616–624. [Google Scholar] [CrossRef]
  7. Hu, G.; Tang, L.; Das, R.; Gao, S.; Liu, H. Acoustic metamaterials with coupled local resonators for broadband vibration suppression. AIP Adv. 2017, 7, 025211. [Google Scholar] [CrossRef] [Green Version]
  8. Koh, C.Y.; Jorba, D.A.; Thomas, E.L. Phononic Metamaterials for Vibration Isolation and Focusing of Elastic Waves. U.S. Patent 8,833,510, 16 September 2014. [Google Scholar]
  9. Wang, Y.-F.; Wang, Y.-S.; Wang, L. Two-dimensional ternary locally resonant phononic crystals with a comblike coating. J. Phys. D Appl. Phys. 2013, 47, 015502. [Google Scholar] [CrossRef]
  10. Chiang, T.-Y.; Wu, L.-Y.; Tsai, C.-N.; Chen, L.-W. A multilayered acoustic hyperlens with acoustic metamaterials. Appl. Phys. A 2011, 103, 355–359. [Google Scholar] [CrossRef]
  11. Oh, J.H.; Seung, H.M.; Kim, Y.Y. Doubly negative isotropic elastic metamaterial for sub-wavelength focusing: Design and realization. J. Sound Vib. 2017, 410, 169–186. [Google Scholar] [CrossRef]
  12. Lu, M.-H.; Zhang, C.; Feng, L.; Zhao, J.; Chen, Y.-F.; Mao, Y.-W.; Zi, J.; Zhu, Y.-Y.; Zhu, S.-N.; Ming, N.-B. Negative birefraction of acoustic waves in a sonic crystal. Nat. Mater. 2007, 6, 744–748. [Google Scholar] [CrossRef]
  13. Cummer, S.A.; Schurig, D. One path to acoustic cloaking. New J. Phys. 2007, 9, 45. [Google Scholar] [CrossRef]
  14. Zigoneanu, L.; Popa, B.-I.; Cummer, S.A. Three-dimensional broadband omnidirectional acoustic ground cloak. Nat. Mater. 2014, 13, 352–355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kushwaha, M.S.; Halevi, P.; Dobrzynski, L.; Djafari-Rouhani, B. Acoustic band structure of periodic elastic composites. Phys. Rev. Lett. 1993, 71, 2022–2025. [Google Scholar] [CrossRef] [PubMed]
  16. Sigalas, M.M.; Economou, E.N. Elastic and acoustic wave band structure. J. Sound Vib. 1992, 158, 377–382. [Google Scholar] [CrossRef]
  17. Brillouin, L. Wave Propagation in Periodic Structures: Electric Filters and Crystal Lattices; Dover Publications: Mineola, NY, USA, 2003; Volume 2. [Google Scholar]
  18. Javid, F.; Wang, P.; Shanian, A.; Bertoldi, K. Architected Materials with Ultra-Low Porosity for Vibration Control. Adv. Mater. 2016, 28, 5943–5948. [Google Scholar] [CrossRef]
  19. Norton, M.P.; Karczub, D.G. Fundamentals of Noise and Vibration Analysis for Engineers; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  20. Zhang, H.; Wen, J.; Xiao, Y.; Wang, G.; Wen, X. Sound transmission loss of metamaterial thin plates with periodic subwavelength arrays of shunted piezoelectric patches. J. Sound Vib. 2015, 343, 104–120. [Google Scholar] [CrossRef]
  21. Vaziri, A.; Hutchinson, J.W. Metal sandwich plates subject to intense air shocks. Int. J. Solids Struct. 2007, 44, 2021–2035. [Google Scholar] [CrossRef] [Green Version]
  22. Lestari, W.; Qiao, P. Dynamic Characteristics and Effective Stiffness Properties of Honeycomb Composite Sandwich Structures for Highway Bridge Applications. J. Compos. Constr. 2006, 10, 148–160. [Google Scholar] [CrossRef]
  23. Liu, Z.; Zhang, X.; Mao, Y.; Zhu, Y.Y.; Yang, Z.; Chan, C.T.; Sheng, P. Locally Resonant Sonic Materials. Science 2000, 289, 1734–1736. [Google Scholar] [CrossRef] [PubMed]
  24. Cai, L.; Xiaoyun, H.; Xisen, W. Band-structure results for elastic waves interpreted with multiple-scattering theory. Phys. Rev. B 2006, 74, 153101. [Google Scholar] [CrossRef]
  25. Dong, Y.; Yao, H.; Du, J.; Zhao, J.; Ding, C. Research on bandgap property of a novel small size multi-band phononic crystal. Phys. Lett. A 2019, 383, 283–288. [Google Scholar] [CrossRef]
  26. Li, S.; Dou, Y.; Chen, T.; Xu, J.; Li, B.; Zhang, F. Designing a broad locally-resonant bandgap in a phononic crystals. Phys. Lett. A 2019, 383, 1371–1377. [Google Scholar] [CrossRef]
  27. Ning, S.; Yan, Z.; Chu, D.; Jiang, H.; Liu, Z.; Zhuang, Z. Ultralow-frequency tunable acoustic metamaterials through tuning gauge pressure and gas temperature. Extrem. Mech. Lett. 2021, 44, 101218. [Google Scholar] [CrossRef]
  28. Bandyopadhyay, A.; Heer, B. Additive manufacturing of multi-material structures. Mater. Sci. Eng. R Rep. 2018, 129, 1–16. [Google Scholar] [CrossRef]
  29. Chen, D.; Zheng, X. Multi-material Additive Manufacturing of Metamaterials with Giant, Tailorable Negative Poisson’s Ratios. Sci. Rep. 2018, 8, 9139. [Google Scholar] [CrossRef]
  30. Abueidda, D.W.; Jasiuk, I.; Sobh, N.A. Acoustic band gaps and elastic stiffness of PMMA cellular solids based on triply periodic minimal surfaces. Mater. Des. 2018, 145, 20–27. [Google Scholar] [CrossRef]
  31. Warmuth, F.; Wormser, M.; Körner, C. Single phase 3D phononic band gap material. Sci. Rep. 2017, 7, 3843. [Google Scholar] [CrossRef] [Green Version]
  32. Gao, P.; Climente, A.; Sánchez-Dehesa, J.; Wu, L. Single-phase metamaterial plates for broadband vibration suppression at low frequencies. J. Sound Vib. 2019, 444, 108–126. [Google Scholar] [CrossRef]
  33. Fei, X.; Jin, L.; Zhang, X.; Li, X.; Lu, M. Three-dimensional anti-chiral auxetic metamaterial with tunable phononic bandgap. Appl. Phys. Lett. 2020, 116, 021902. [Google Scholar] [CrossRef]
  34. Huang, Y.; Li, J.; Chen, W.; Bao, R. Tunable bandgaps in soft phononic plates with spring-mass-like resonators. Int. J. Mech. Sci. 2019, 151, 300–313. [Google Scholar] [CrossRef]
  35. Chen, M.; Xu, W.; Liu, Y.; Yan, K.; Jiang, H.; Wang, Y. Band gap and double-negative properties of a star-structured sonic metamaterial. Appl. Acoust. 2018, 139, 235–242. [Google Scholar] [CrossRef] [Green Version]
  36. Chen, M.; Jiang, H.; Zhang, H.; Li, D.; Wang, Y. Design of an acoustic superlens using single-phase metamaterials with a star-shaped lattice structure. Sci. Rep. 2018, 8, 1861. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Ren, F.; Wang, L.; Liu, H. Low frequency and broadband vibration attenuation of a novel lightweight bidirectional re-entrant lattice metamaterial. Mater. Lett. 2021, 299, 130133. [Google Scholar] [CrossRef]
  38. Kumar, N.; Pal, S. Low frequency and wide band gap metamaterial with divergent shaped star units: Numerical and experimental investigations. Appl. Phys. Lett. 2019, 115, 254101. [Google Scholar] [CrossRef]
  39. Isik, O.; Esselle, K.P. Analysis of spiral metamaterials by use of group theory. Metamaterials 2009, 3, 33–43. [Google Scholar] [CrossRef]
  40. Qi, D.; Yu, H.; Hu, W.; He, C.; Wu, W.; Ma, Y. Bandgap and wave attenuation mechanisms of innovative reentrant and anti-chiral hybrid auxetic metastructure. Extrem. Mech. Lett. 2019, 28, 58–68. [Google Scholar] [CrossRef]
Figure 1. (a) Description of the square honeycomb acoustic metamaterials with locally resonant Archimedean spirals (SHAMLRAS); (b) the unit cell of SHAMLRAS; (c) the circular array unit of Archimedes spiral; (d) the first Brillouin zone and irreducible Brillouin zone (IBZ, black shadow region).
Figure 1. (a) Description of the square honeycomb acoustic metamaterials with locally resonant Archimedean spirals (SHAMLRAS); (b) the unit cell of SHAMLRAS; (c) the circular array unit of Archimedes spiral; (d) the first Brillouin zone and irreducible Brillouin zone (IBZ, black shadow region).
Materials 15 00373 g001
Figure 2. (a) Band structure of the SHMLRAS and (b) detailed schematic of the first bandgap.
Figure 2. (a) Band structure of the SHMLRAS and (b) detailed schematic of the first bandgap.
Materials 15 00373 g002
Figure 3. Mode shapes of the proposed structure: (ac) represent the points (O, A, B) of the IBZ on the fourth branches the lower edge of the first bandgap; (df) represent the points (O, A, B) of the IBZ on the third branches the upper edge of the first bandgap; the arrows represent the magnitude and direction of the displacement of the mode. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.).
Figure 3. Mode shapes of the proposed structure: (ac) represent the points (O, A, B) of the IBZ on the fourth branches the lower edge of the first bandgap; (df) represent the points (O, A, B) of the IBZ on the third branches the upper edge of the first bandgap; the arrows represent the magnitude and direction of the displacement of the mode. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.).
Materials 15 00373 g003
Figure 4. Mode shapes of the proposed structure: (ac) represent the points (O, A, B) of the IBZ on the sixth branches the lower edge of the second bandgap; (df) represent the points (O, A, B) of the IBZ on the ninth branches the upper edge of the third bandgap; (gi) represent the points (O, A, B) of the IBZ on the seventh branches the upper edge of the second bandgap; (jl) represent the points (O, A, B) of the IBZ on the eighth branches the lower edge of the third bandgap; the arrows represent the magnitude and direction of the displacement of the mode.
Figure 4. Mode shapes of the proposed structure: (ac) represent the points (O, A, B) of the IBZ on the sixth branches the lower edge of the second bandgap; (df) represent the points (O, A, B) of the IBZ on the ninth branches the upper edge of the third bandgap; (gi) represent the points (O, A, B) of the IBZ on the seventh branches the upper edge of the second bandgap; (jl) represent the points (O, A, B) of the IBZ on the eighth branches the lower edge of the third bandgap; the arrows represent the magnitude and direction of the displacement of the mode.
Materials 15 00373 g004
Figure 5. (a) Dispersion relations for SHAMLRAS structures in isometric perspective; (b) comparison of dispersion relations and band structure in horizontal view.
Figure 5. (a) Dispersion relations for SHAMLRAS structures in isometric perspective; (b) comparison of dispersion relations and band structure in horizontal view.
Materials 15 00373 g005
Figure 6. (a) Third mode; (b) ISO-frequency contours (the thick black solid lines are the three selected frequencies to be analyzed and similar depictions can be seen in Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12); (ce) group velocity of third wave mode at 183.3 Hz, 183.6 Hz, and 184.2 Hz.
Figure 6. (a) Third mode; (b) ISO-frequency contours (the thick black solid lines are the three selected frequencies to be analyzed and similar depictions can be seen in Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12); (ce) group velocity of third wave mode at 183.3 Hz, 183.6 Hz, and 184.2 Hz.
Materials 15 00373 g006
Figure 7. (a) Fourth mode; (b) ISO-frequency contours; (ce) group velocity of fourth wave mode at 185.5 Hz, 186.4 Hz, and 187.8 Hz.
Figure 7. (a) Fourth mode; (b) ISO-frequency contours; (ce) group velocity of fourth wave mode at 185.5 Hz, 186.4 Hz, and 187.8 Hz.
Materials 15 00373 g007
Figure 8. (a) Fifth mode; (b) ISO-frequency contours; (ce) group velocity of fifth wave mode at 198.2 Hz, 208.2 Hz, and 223.2 Hz.
Figure 8. (a) Fifth mode; (b) ISO-frequency contours; (ce) group velocity of fifth wave mode at 198.2 Hz, 208.2 Hz, and 223.2 Hz.
Materials 15 00373 g008
Figure 9. (a) Sixth mode; (b) ISO-frequency contours; (ce) group velocity of sixth wave mode at 211.2 Hz, 231.3 Hz, and 233.9 Hz.
Figure 9. (a) Sixth mode; (b) ISO-frequency contours; (ce) group velocity of sixth wave mode at 211.2 Hz, 231.3 Hz, and 233.9 Hz.
Materials 15 00373 g009
Figure 10. (a) Seventh mode; (b) ISO-frequency contours; (ce) group velocity of seventh wave mode at 250.6 Hz, 251.3 Hz, and 252.1 Hz.
Figure 10. (a) Seventh mode; (b) ISO-frequency contours; (ce) group velocity of seventh wave mode at 250.6 Hz, 251.3 Hz, and 252.1 Hz.
Materials 15 00373 g010
Figure 11. (a) Eighth mode; (b) ISO-frequency contours; (ce) group velocity of eighth wave mode at 252.1 Hz, 252.4 Hz, and 252.8 Hz.
Figure 11. (a) Eighth mode; (b) ISO-frequency contours; (ce) group velocity of eighth wave mode at 252.1 Hz, 252.4 Hz, and 252.8 Hz.
Materials 15 00373 g011
Figure 12. (a) Ninth mode; (b) ISO-frequency contours; (ce) group velocity of ninth wave mode at 277.6 Hz, 294.3 Hz, and 311 Hz.
Figure 12. (a) Ninth mode; (b) ISO-frequency contours; (ce) group velocity of ninth wave mode at 277.6 Hz, 294.3 Hz, and 311 Hz.
Materials 15 00373 g012
Figure 13. The dependence of the bandgap width on (a,b) the circle distance d, and (c,d) the inner diameter R2.
Figure 13. The dependence of the bandgap width on (a,b) the circle distance d, and (c,d) the inner diameter R2.
Materials 15 00373 g013
Figure 14. The dependence of the bandgap width on (a,b) the turn number n, (c,d) the start value t0, and (e,f) the end value t1.
Figure 14. The dependence of the bandgap width on (a,b) the turn number n, (c,d) the start value t0, and (e,f) the end value t1.
Materials 15 00373 g014
Figure 15. The dependence of the bandgap width on the spiral (a,b) transforming quadrants, where the number “original” represents the bandgap distribution for the parameter shown in Table 2, the number “I” represents a transformation of the spiral arrangement in the first quadrant, the numbers “I&II” represents a transformation of the spiral arrangement in the first and second quadrants, and the numbers “I&III” represents a transformation of the spiral arrangement in the first and third quadrants.
Figure 15. The dependence of the bandgap width on the spiral (a,b) transforming quadrants, where the number “original” represents the bandgap distribution for the parameter shown in Table 2, the number “I” represents a transformation of the spiral arrangement in the first quadrant, the numbers “I&II” represents a transformation of the spiral arrangement in the first and second quadrants, and the numbers “I&III” represents a transformation of the spiral arrangement in the first and third quadrants.
Materials 15 00373 g015
Figure 16. The dependence of the bandgap width on (a,b) Young’s modulus E, (c,d) the Poisson’s ratio ν, and (e,f) the density ρ.
Figure 16. The dependence of the bandgap width on (a,b) Young’s modulus E, (c,d) the Poisson’s ratio ν, and (e,f) the density ρ.
Materials 15 00373 g016
Figure 17. (a) Geometry and loading conditions of the SHAMLRAS periodic structure; (b) experimental setup for transmission loss in SHAMLRAS periodic structures; (c) comparative spectrum diagram obtained from transmission loss experiment and simulation; (d) deformation of SHAMLRAS periodic structures under harmonic loading frequency at f = 220 Hz (outside the bandgap region) and f = 240 Hz (inside the bandgap region).
Figure 17. (a) Geometry and loading conditions of the SHAMLRAS periodic structure; (b) experimental setup for transmission loss in SHAMLRAS periodic structures; (c) comparative spectrum diagram obtained from transmission loss experiment and simulation; (d) deformation of SHAMLRAS periodic structures under harmonic loading frequency at f = 220 Hz (outside the bandgap region) and f = 240 Hz (inside the bandgap region).
Materials 15 00373 g017
Table 1. Coordinates of the boundary points of the IBZ.
Table 1. Coordinates of the boundary points of the IBZ.
Boundary PointsCartesian BasisReciprocal Basis
O(0, 0)(0, 0)
A ( π / L , 0)(1, 0)
B ( π / L ,   π / L )(1, 1)
Table 2. Geometry and material parameters of the structure.
Table 2. Geometry and material parameters of the structure.
Lattice ParametersMaterial Parameters of Photosensitive Resin
Radius of tangent circleR1 = 9.5 mmYoung’s modulus2.642 GPa
Inner diameterR2 = 2.0 mmDensity1150 kg/m3
Circle distanced = 4.25 mmPoisson’s ratio0.42
Ligament thicknessp = 1.0 mm
Turn numbern = 2
Start valuet0 = 0
End valuet1 = 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gao, H.; Yan, Q.; Liu, X.; Zhang, Y.; Sun, Y.; Ding, Q.; Wang, L.; Xu, J.; Yan, H. Low-Frequency Bandgaps of the Lightweight Single-Phase Acoustic Metamaterials with Locally Resonant Archimedean Spirals. Materials 2022, 15, 373. https://doi.org/10.3390/ma15010373

AMA Style

Gao H, Yan Q, Liu X, Zhang Y, Sun Y, Ding Q, Wang L, Xu J, Yan H. Low-Frequency Bandgaps of the Lightweight Single-Phase Acoustic Metamaterials with Locally Resonant Archimedean Spirals. Materials. 2022; 15(1):373. https://doi.org/10.3390/ma15010373

Chicago/Turabian Style

Gao, Haoqiang, Qun Yan, Xusheng Liu, Ying Zhang, Yongtao Sun, Qian Ding, Liang Wang, Jinxin Xu, and Hao Yan. 2022. "Low-Frequency Bandgaps of the Lightweight Single-Phase Acoustic Metamaterials with Locally Resonant Archimedean Spirals" Materials 15, no. 1: 373. https://doi.org/10.3390/ma15010373

APA Style

Gao, H., Yan, Q., Liu, X., Zhang, Y., Sun, Y., Ding, Q., Wang, L., Xu, J., & Yan, H. (2022). Low-Frequency Bandgaps of the Lightweight Single-Phase Acoustic Metamaterials with Locally Resonant Archimedean Spirals. Materials, 15(1), 373. https://doi.org/10.3390/ma15010373

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop