# Quasi-static stop band with flexural metamaterial having zero rotational stiffness

Joo Hwan Oh, Badreddine Assouar

<figure><img src="image_1.png" /></figure>

To cite this version:

Joo Hwan Oh, Badreddine Assouar. Quasi-static stop band with flexural metamaterial having zero rotational stiffness. Scientific Reports, 2016, 6 (1), pp.33410. 10.1038/srep33410. hal-04953644

HAL Id: hal-04953644

https://hal.science/hal-04953644v1

Submitted on 18 Feb 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Quasi-static stop band with flexural metamaterial

having zero rotational stiffness

Joo Hwan Oh<sup>a,b,c</sup> and Badreddine Assouar<sup>a,b</sup>*

<sup>a</sup> University of Lorraine, Institut Jean Lamour, Boulevard des Aiguillettes, BP: 70239, Vandoeuvre-lès-Nancy 54506, France

<sup>b</sup> CNRS, Institut Jean Lamour, Vandoeuvre-lès-Nancy 54506, France

<sup>c</sup> Institute of Advanced Machine and Design, Seoul National University, 599 Gwanak-ro, Gwanak-gu, Seoul, 151-744, Korea

* Corresponding Author, Email: Badreddine.Assouar@univ-lorraine.fr; Tel: +33 3 83 68 49 05

Metamaterials realizing stop bands have attracted much attentions recently since they can break-through the well-known mass law. However, achieving the stop band at extremely low frequency has been still a big challenge in the fields of elastic metamaterials. In this paper, we propose a new metamaterial based on the idea of the zero rotational stiffness, to achieve extremely low frequency stop band for flexural elastic waves. Unlike the previous ways to achieve the stop band, we found that the zero rotational stiffness can provide a broad stop band at extremely low frequency, which starts from even almost zero frequency. To achieve the zero rotational stiffness, we propose a new elastic metamaterial consisting of blocks and links with the hinge connection. Analytic developments as well as numerical simulations evidence that this new metamaterial can exhibit extremely low and broad stop band, even at the quasi-static ranges. In addition, the metamaterial is shown to exhibit the negative group velocity at extremely low frequency ranges, as well as the quasi-static stop band, if it is properly designed.
Keywords: Elastic Metamaterial, Zero Rotational Stiffness, Low Frequency Stop band, Flexural Elastic Wave.

Introduction

Recently, researches on the artificial sub-wavelength acoustic or elastic structures, called acoustic/elastic metamaterials, have extended our knowledge beyond the current limit. By utilizing the intrinsic dynamic motions inside their unit cell structures, metamaterials have shown various interesting wave phenomena, such as the negative refraction or even the wave cloaking, which were impossible to be realized previously. Above all, the metamaterial can be designed to have the stop band that can break-through the well-known mass-law, which is very important characteristic in various vibration/acoustic devices. In fact, the break-through of the mass-law has been already shown in the phononic crystals¹,², which utilize the Bragg scattering to achieve a stop band originated from the crystal’s periodicity³,⁴. However, there has been a limitation that the stop bands in phononic crystals generally exist at too high frequencies to be applied in vibration or acoustic devices.

On the other hand, recent researches on the resonance-based metamaterials have realized a stop band at relatively low frequencies where the wavelength is much larger than the size of the unit cell. From the first realization of the resonance-based acoustic metamaterial by Ping Sheng et al.⁵, there have been a number of researches on the resonance-based metamaterials. For the in-plane elastic waves, Huang et al.⁶ developed an analytic theory to explain the formation of the negative density with the resonance-based metamaterials. Zhu et al.⁷ realized the anisotropic elastic metamaterial having different wave stop band for different directions. Hou et al.⁸ proposed the idea of the tunable elastic metamaterial having negative stiffness. Zhu et al.⁹,¹⁰ realized the chiral elastic metamaterials having broad stop band and double negativity, respectively. Romero-Garcia et al.¹¹ realized the resonance-based acoustic
metamaterial having various stop band. Recently, Oh *et al.*¹²,¹³ realized the elastic metamaterial whose negative density and stiffness can be independently realized and adjoin each other, respectively. For the out-of-plane elastic waves (flexural elastic waves), Oudich *et al.*¹⁴ realized the locally resonant metamaterial with the stubbed plates. Williams *et al.*¹⁵ developed the multiresonant metamaterials for the flexural waves. Assouar *et al.*¹⁶ proposed the double-sides stubbed plates for broad stop band. The latter extended the idea of the double-sides stubbed plates to the hybrid plates having both cylinders and holes to further enlarge the stop band¹⁷. In addition, there have been several researches to achieve the extremely low frequency stop band without the internal resonance phenomena. Li *et al.*¹⁸ utilized the concept of the acoustic coiling-up space to achieve perfect absorption and zero transmission at very low frequency without the internal resonance metamaterial. Oh *et al.*¹⁹ used the idea of the corrugated elastic structure to obtain broad stop band at very low frequency ranges, and extended the idea to achieve the far field sub-wavelength imaging lens.

Despite of the previous researches on the low frequency stop band metamaterials, challenge still remains on achieving much lower frequency stop band, especially the stop band covering almost quasi-static frequency ranges. In fact, there have been only few researches on the metamaterials that allow the quasi-static stop band. In acoustics, Lee *et al.*²⁰,²¹ proposed and realized the acoustic metamaterial having negative density and bulk modulus at the quasi-static frequencies, respectively. However, due to the physical difference between acoustics (based on the vector physics) and elasticity (based on the tensor physics), the acoustic metamaterials cannot be directly used to achieve the elastic quasi-static stop band. In elasticity, Yao *et al.*²² and Yu *et al.*²³ utilized the idea of the elastic foundation to realize the elastic quasi-static stop band. Nevertheless, the elastic metamaterials require the connection with the external fixed boundary to achieve the quasi-static wave stop band. For the best of the author's knowledge, there has been no research to obtain the elastic
metamaterial with the quasi-static stop band by its intrinsic dynamic characteristics only.

In this research, we proposed a new metamaterial to obtain the elastic quasi-static stop band for out-of-plane flexural elastic waves. The main idea in this research is called ‘the zero rotational stiffness,’ which is completely different from the previous metamaterials based on the internal resonance, the coiling up space or the Bragg scattering. Here, from the classical Timoshenko beam theory, we found that the broad quasi-static stop band metamaterial can be achieved if the metamaterial has a zero rotational stiffness. Based on this idea, we newly designed a new elastic metamaterial, consisting of the hinge connection based structure with blocks and links, to achieve the zero rotational stiffness. The analytic theory and the numerical simulations for the metamaterial showed that the metamaterial has broad quasi-static stop band from almost zero frequency. Moreover, with proper design, we found that the metamaterial can be tuned to exhibit negative group velocities in addition to the quasi-static stop band. Since extremely low frequency stop band can be obtained with this research, the findings in this paper are expected to be further applied in various vibration and flexural wave devices.

## Results

**Background theory: mass-spring system for general flexural waves**

Before showing the main idea of the zero rotational stiffness, it is worth to review the periodic mass-spring system for the general flexural wave. Consider a flexural wave propagating through a slender beam as in Fig. 1 (a). To analyze the wave dispersion relation of the flexural wave, this slender beam can be converted to the equivalent mass-spring system shown in Fig. 1 (b). It should be emphasized that each mass in Fig. 1 (b) are connected with two springs. In acoustics or in-plane elastic waves, only one spring connecting each mass is considered in the equivalent mass-spring system. However, as well known from the classical Timoshenko beam
theory, the flexural wave is governed by two types of motions, the material's vertical movement *w* and the rotation *θ*²⁴. Accordingly, in the equivalent mass-spring system of the flexural waves, both springs, the vertical shear spring *α* (which represents the shear stiffness of the beam) and the rotational spring *β* (which represents the bending stiffness of the beam)²⁵, should be considered between each mass to properly model the flexural waves.

Also, both the mass *m* and the momentum of inertia *I* are considered for each mass.

Following the detailed procedures shown in the Supplementary Material, the equation of motion of the equivalent mass-spring system can be written as

$$ -\omega^2 \begin{bmatrix} m w_n \\ I \theta_n \end{bmatrix} = \begin{bmatrix} A & B \\ \operatorname{conj}(B) & C \end{bmatrix} \begin{bmatrix} w_n \\ \theta_n \end{bmatrix}, \qquad (1a) $$

$$ A = \alpha(\exp(-ika) + \exp(ika) - 2), \qquad (1b) $$

$$ B = 0.5a^* \alpha(\exp(ika) - \exp(-ika)), \qquad (1c) $$

$$ C = \beta(\exp(-ika) + \exp(ika) - 2) - (0.5a^*)^2 \alpha(\exp(-ika) + \exp(ika) + 2). \qquad (1d) $$

Here, *ω* and *k* are the angular frequency and the wavevector, respectively.

In homogeneous continuous beam structure, equation (1) can be further reduced by using the fact that $a^* = a \ll 1$. In this case, $\exp(ika) + \exp(-ika) \sim 2-(ka)^2$ and $\exp(ika) - \exp(-ika) \sim 2ika$ can be assumed. As a result, the wave dispersion curve can be obtained as (see the Supplementary Material for the detailed procedures)

$$ \alpha \beta a^2 k^4 - a \rho \omega^2 (\alpha I_b + \beta A_b) k^2 + \rho A_b \omega^2 (\rho I_b \omega^2 - \alpha a) = 0. \qquad (2) $$

As can be seen in equation (2), unlike those in acoustics or in-plane elastic waves, the wave dispersion equation of the flexural wave has a fourth order equation. Fig. 1 (c) shows the plot of the dispersion curves shown in equation (2). Here, $A_b, I_b, \rho, a, \alpha$ and $\beta$ are considered to be 1. As can be seen here, there are two wave dispersion curves for general flexural waves. The first dispersion branch starts from the origin, i.e., the zero frequency. On
the other hand, the second dispersion branch starts to appear above the certain cutoff frequency.

**Main idea of the zero rotational stiffness** Based on the wave dispersion relation of general flexural wave, let us clarify our idea. The main idea in this work, the zero rotational stiffness, is to make the rotational spring $\beta$ to be zero. From equation (2), it can be clearly seen that if the rotational spring $\beta$ is zero, the 4th order term (the $k^4$ term) in the wave dispersion equation becomes zero and the wave dispersion equation is reduced to the second order equation as

$$-\omega^2(\alpha a \rho I_b k^2 + \alpha a \rho A_b - \rho^2 \omega^2 A_b I_b) = 0. \quad (3)$$

Accordingly, the wave dispersion equations for the zero rotational stiffness are

$$\omega = 0 \text{ (for all } k) \text{ or } k = \left( \frac{A_b}{I_b} - \frac{\rho\omega}{\alpha a} \omega^2 \right)^{1/2}. \quad (4)$$

Fig. 1 (d) shows the wave dispersion equation with $\beta=0$ while the other constants are same as in Fig. 1 (c). From Fig. 1 (d), one can clearly see that if $\beta=0$, the first branch becomes totally flat while the second branch is almost unaltered. Thus, one can obtain the quasi-static stop band below the cutoff frequency of the second branch.

As explained, the idea of the zero rotational stiffness is a totally different mechanism compared to the Bragg scattering and the local resonance phenomena that have been utilized to achieve the stop band previously. The previous ideas of the Bragg scattering and the local resonance phenomena are to create a stop band at a frequency range which was pass band naturally. However, in this research, we don't create a stop band; we make the first dispersion branch to be a flat branch so that we can utilize the existing cutoff frequency of the second branch as a stop band. Also, the idea of the zero rotational stiffness is valid for the flexural
wave only. In acoustic or in-plane elastic waves, if the bulk modulus or the corresponding
stiffness is zero, the stiffness term becomes zero for the whole frequency and the wave
motion, such as the wave dispersion equation cannot be defined. However, because the
flexural wave is governed by both the vertical shear spring and the rotational spring, although
the rotational spring is zero, each component in the stiffness matrix is not zero unless the
vertical shear spring is zero too. Thus, the wave dispersion equation can still be defined
although $\beta = 0$.

**Metamaterial design to achieve the zero rotational stiffness**

There can be various metamaterial designs to achieve the zero rotational stiffness. Here, among the various ways, we consider a hinge connection based structure to achieve the zero rotational stiffness. The metamaterial design consisting of the hinge connection is shown in Fig. 2. In the current metamaterial design, the most important part is the rigid cylinder that connects the blocks and the links. The rigid cylinder connecting the blocks and the links makes two structures move together along x, y and z directions. However, the rigid cylinder does not provide any constraints for their rotation with respect to the cylinder. Thus, if any moment is induced, the links and the blocks do not deform but rotate independently. Therefore, no moment can be transferred between the blocks and the links, i.e., the structure has the zero rotational stiffness.

**Analytic modeling of the metamaterial**

Based on the proposed metamaterial, let us focus on the analytic investigation of the metamaterial first. To carry out this analytic investigation, the metamaterial shown in Fig. 2 is replaced by the equivalent mass-spring system as in Fig. 3 (a). Here, we will focus on the linear case only when the motion of the blocks is small and thus the change of the angle of the links $\Phi$ is ignorable during the
motion. In the equivalent mass-spring system in Fig. 3 (b), it is assumed that the dominant deformation takes place at the links while the blocks exhibit almost zero deformation. This assumption is usually valid since the links generally have a low stiffness than the blocks. Also, the rotational springs are not used to describe the links since the links are connected with the blocks by the hinge connection and thus no moment can be transferred from the links to the blocks. As a result, the links are modeled as the longitudinal spring with the spring coefficient $\gamma$. From the detailed analytic procedures shown in the Supplementary Material, the wave dispersion equation can be derived as

$$-\omega^2 \begin{bmatrix} mw_n \\ I\theta_n \end{bmatrix} = \begin{bmatrix} A' & B' \\ conj(B') & C' \end{bmatrix} \begin{bmatrix} w_n \\ \theta_n \end{bmatrix}, \qquad (5a)$$

$$A' = 2\gamma \sin^2 \Phi(\exp(-ika) + \exp(ika) - 2), \qquad (5b)$$

$$B' = 2\gamma \sin \Phi(b \sin \Phi + d \cos \Phi)(\exp(ika) - \exp(-ika)), \qquad (5c)$$

$$C' = -2\gamma(b \sin \Phi + d \cos \Phi)^2 (\exp(-ika) + \exp(ika) + 2). \qquad (5d)$$

Comparing equation (5) with the general wave dispersion equation in equation (1), one can find that two equations are in fact identical to each other with the following equivalences;

$$\alpha = 2\gamma \sin^2 \Phi, \quad a^* = 2(b \sin \Phi + d \cos \Phi) / \sin \Phi, \quad \beta = 0. \qquad (6)$$

Here, one can clearly see that the metamaterial is equivalent to the general mass-spring system with $\beta = 0$, i.e., the proposed unit cell structure has indeed the zero rotational stiffness. Also, it can be seen that the angle of the links $\Phi$ should be non-zero value; otherwise, the stiffness matrix in equation (5) becomes zero for all frequencies and the wave phenomena inside the metamaterial cannot be defined.

**Numerical methodology for the metamaterial**

Now, the numerical methodology to evaluate the wave dispersion curves of the continuum unit cell structure is focused. Fig. 4 (a)
shows the actual continuum model used in the current numerical investigation. Here, aluminum (Young's modulus = 70 GPa, Poisson's ratio = 0.33 and density = 2700 kg/m³) is considered for the blocks and the polytetrafluoroethylene (PTFE, Young's modulus = 0.5 GPa, Poisson's ratio = 0.49 and density = 2000 kg/m³) is considered for the links. As used for general metamaterials, the Floquet-Bloch boundary condition is applied at the sides of the unit cell in Fig. 4 (a) and the eigenvalue problem is solved for various wavevectors *k* inside the 1st irreducible Brillouin zone²⁶. In the simulation, the structural module (linear analysis) of COMSOL 4.3b is used.

In the current numerical simulation, a special attention should be paid to properly model the hinge connection condition. To explain how the hinge connection is modeled, let us consider the hinge connection between only one block and link, as in Fig. 4 (b). Here, Boundary 1 indicates the contact boundary of the block with the rigid cylinder, and Boundary 2, that of the link with the rigid cylinder. To introduce the hinge connection condition, all displacements of the nodes at the Boundary 1 are expressed by the displacements and the rotations of the point located at the center of the Boundary 1, respectively. For instance, the displacement **U** = (u, v, w)ₙ of any node at the Boundary 1 can be expressed as

$$ \mathbf{U} = (u, v, w)_1 + (\theta_x, \theta_y, \theta_z)_1 \times \mathbf{r} \quad (7) $$

where (u, v, w)₁ and (θx, θy, θz)₁ are the displacements and the rotations of the center of the Boundary 1, as in Fig. 4 (b). Also, **r** is the vector which connects from the center to the node in the Boundary 1. In the same manner, all the displacements of the nodes at the Boundary 2 are expressed by (u, v, w)₂ and (θx, θy, θz)₂, which are defined at the center of the Boundary 2. After that, the following constraints are imposed between two center points of the Boundary 1 and 2 to achieve the hinge connection condition;

$$ (u, v, w)_2 = (u, v, w)_1 + (\theta_x, 0, \theta_z)_1 \times \mathbf{r}_{1-2}, \quad (\theta_x, 0, \theta_z)_2 = (\theta_x, 0, \theta_z)_1. \quad (8) $$
Here, $\mathbf{r}_{1-2}$ is the vector which connects from the center of the Boundary 1 to the center of the Boundary 2. As can be seen in equation (8), the displacements and the rotations of the center of the Boundary 1 are set to be identical to those of the center of the Boundary 2, except the rotation with respect to the y axis, $(\theta_y)_1$ and $(\theta_y)_2$. Therefore, the block and the link can freely rotate independently with respect to the y axis, i.e., $(\theta_y)_1$ and $(\theta_y)_2$ are perfectly independent to each other. By this way, one can model the hinge connection condition numerically.

**Validation of the analytic investigations** By using the numerical method, the validation of the analytic theory is made by comparing the wave dispersion curve obtained by the analytic mass-spring system with the numerically calculated wave dispersion curve. The detailed geometric parameter of the unit cell, and the mass and spring coefficients used for the analytic calculation, are given in the Supplementary Material.

Fig. 5 (a) shows the wave dispersion curves obtained by the analytic mass-spring system (equation (5)) and the numerical simulation. Note that in Fig. 5 (a), only the flexural wave mode's wave dispersion curves are plotted among various possible elastic wave modes. Good agreements can be found between two results, validating the analytic equations derived with the equivalent mass-spring system. Also, it can be clearly seen that very broad stop band is indeed formed at extremely low frequency. The quasi-static stop band is measured to be formed from 0.006 Hz to 425 Hz. This clearly shows that the proposed elastic metamaterial has extremely low frequency stop band, owing to the zero rotational stiffness.

**Mode shape investigation of the proposed metamaterial** For a better understanding of the underlying physics, the mode shapes for each dispersion branch are studied. Fig. 5 (b)
shows the mode shapes at the point A₁ and A₂ in Fig. 5 (a), respectively. From the left figure in Fig. 5 (b), it can be clearly seen that the links have almost same displacement, i.e., the links do not exhibit any deformation. Considering that the links act as the spring which dominates the overall stiffness of the unit cell, almost no deformation at links indicates that the mode shape involves almost zero stiffness. Thus, the frequency of the first dispersion branch should be zero, as expected from the analytic investigation. In fact, the flat dispersion branch at zero frequency indicates that the unit cell should have an instability under the zero frequency (static condition). If a metamaterial has the zero rotational stiffness, each block can freely rotate and thus the structure cannot sustain its own weight. However, in reality, the perfectly zero rotational stiffness cannot be realized due to the friction and thus the metamaterial can be stable. For more about this issue, see the Supplementary Material.

On the other hand, in the right figure in Fig. 5 (b), the deformation varies in the links, showing that there is high deformation induced in the links. Therefore, the mode shape involves nonzero stiffness and the second dispersion branch should be located at nonzero frequencies, as shown in Fig. 5 (a). Note that in the right figure in Fig. 5 (b), one can see that the links have also some internal bending motion. This bending motion provides additional moments and shear forces, which were not considered in the analytic model. Thus, in strict meaning, the metamaterial may not provide ‘purely zero’ rotational stiffness. However, the effects of these moments and shear forces are extremely small and ignorable; due to the hinge connection condition, the moments are not transferred to the block, and the shear forces only provide very small rotation of the block. This insignificance of the internal bending motion can be supported from the good agreements between the dispersion curves shown in Fig. 5 (a); although the internal bending motion is not considered, the analytic approach can precisely predict the actual wave dispersion curve of the continuum unit cell structure. This indicates that the internal bending motion is not significant in the current metamaterial,
especially in the frequency range of the interest.

**Wave transmission analysis of the proposed metamaterial**

To further validate the formation of the quasi-static stop band, the wave transmission is numerically calculated with the same program, the structural module (linear analysis) of COMSOL 4.3b. In the simulation, various metamaterial systems consisting of 1, 2, 3 and 4 unit cells, respectively, are considered. Fig. 6 (a) shows the continuum model consisting of 4 unit cells used for the wave transmission analysis. As in Fig. 6 (a), the z-directional force is imposed at the left end and the z-directional displacements, $w_1$ and $w_2$, are measured at the points $P_1$ and $P_2$, respectively. After that, the wave transmission is evaluated as $abs(w_2 / w_1)$ for various frequencies. Although this definition of the wave transmission does not provide the exact transmission coefficient, it can clearly show the existence of the stop band. For the numerical efficiency, only the half of the structure along the y direction is considered and the symmetric boundary condition along the y direction is applied as in Fig. 6 (a). Also, the frequencies over 100 Hz are considered for all simulations since there is a convergence problem at the lower frequencies due to the instability of the metamaterial.

Fig. 6 (b) shows the numerically calculated wave transmission for the metamaterial system consisting of various unit cells. As can be seen in the simulation results, very low wave transmissions are measured for all metamaterial system at the frequencies from 100 Hz to 425 Hz, validating the formation of the quasi-static stop band. Also, it can be seen that only few unit cells can provide very low wave transmission at the stop band frequency ranges. For instance, the wave transmission below 0.01 at 100 Hz can be achieved with only two unit cells. This suggests that the stop band in the metamaterial is not originated from the Bragg scattering. In addition, several peaks can be observed in the wave pass band from 425 Hz to
538 Hz. Considering that the number of the peaks are same as the number of the unit cells used in the metamaterial system, it can be seen that the peaks are originated from the Fabry-Perot resonance inside the metamaterial system.

To more clearly see why the stop band and the pass band are formed, the deformation characteristics of the metamaterial layer at 200 Hz (belongs to the stop band) and 500 Hz (belongs to the pass band) are shown in Figs. 6 (c, d). In Fig. 6 (c), it can be seen that the dominant motion of the blocks at 200 Hz is the rotation with respect to the y axis. Since the links do not have any stiffness for the rotation of the blocks, no wave energy can be transmitted through the metamaterial. However, at 500 Hz, the blocks exhibit not only the rotation but also the vertical movement, as can be seen in Fig. 6 (d). In this case, the links exhibit internal deformation according to the vertical movement of the blocks. Thus, the wave energies can be transmitted between the neighboring blocks, making the pass band at the frequency.

In fact, from the classical Timoshenko theory, it has been well-known that the flexural wave is dominantly governed by the rotation at low frequencies, while it is governed by the coupled motion of the rotation and the vertical movement at high frequencies. Since the metamaterial has the zero rotational stiffness, it cannot transfer the wave energy at low frequencies where the wave is dominated by the rotation. However, at high frequencies, the vertical movement begins to take place and the metamaterial can transfer the wave energy due to the vertical movement. This is how the quasi-static stop band is formed with the proposed elastic metamaterial.

**Expansion to the negative group velocity** In addition to the quasi-static stop band, we will show that the metamaterial can be tailored to achieve the negative group velocity. To clarify this point, equation (5) is re-visited. From the latter, the frequencies at $k=0$ and
$k = \pi / a$ can be analytically calculated as

$$ \text{at } k = 0: \omega_1^{k=0} = 0 \text{ (1st branch) or } \omega_2^{k=0} = (b \sin \Phi + d \cos \Phi)\sqrt{8\gamma / I} \text{ (2nd branch), (9a)} $$

$$ \text{at } k = \pi / a: \omega_1^{k=\pi/a} = 0 \text{ (1st branch) or } \omega_2^{k=\pi/a} = \sin \Phi \sqrt{8\gamma / m} \text{ (2nd branch). (9b)} $$

In the previously considered metamaterial, $\omega_2^{k=0}$ was smaller than $\omega_2^{k=\pi/a}$ and the second dispersion branch had the positive group velocity. However, if $\omega_2^{k=0}$ is larger than $\omega_2^{k=\pi/a}$, one can expect to have the negative group velocity at the second dispersion branch. Considering equation (S17), one can clearly see that $\omega_2^{k=\pi/a}$ decreases as $l$ increases since it will increase the total mass $m$. For $\omega_2^{k=0}$, on the other hand, increasing $l$ will increase $I$ but it also increases $b$. In fact, in the current metamaterial, $\omega_2^{k=0}$ increases as $l$ increases. Thus, to achieve $\omega_2^{k=0} > \omega_2^{k=\pi/a}$, we enlarge the length of the block as in Fig. 7 (a). Figs. 7 (b, c) show the wave dispersion curve and the mode shapes, respectively. As can be seen in Fig. 7 (b), the new metamaterial is shown to have the negative group velocity from 369 Hz to 465 Hz, in addition to the quasi-static stop band from 0 Hz to 369 Hz. This clearly supports that the proposed metamaterial can even realize negative group velocity.

Fig. 8 (a) shows the numerical model to calculate the wave transmission of the metamaterial consisting of the unit cell shown in Fig. 7 (a). Here, the detailed procedures are same as that used in Fig. 6 (a). The transmission plot shown in Fig. 8 (b) clearly supports the wave dispersion curve in Fig. 7 (b); the existence of the quasi-static stop band from 100 to 369 Hz and the pass band from 369 Hz to 465 Hz. The important point is that the formation of the negative group velocity does not affect the formation of the quasi-static stop band. In fact, from the deformation characteristics at 200 Hz (belongs to the stop band) and 400 Hz (belongs to the pass band) shown in Fig. 8 (c), it can be seen that regardless of the formation of the negative group velocity, the metamaterial exhibits similar dynamic behavior as the
previous metamaterial system shown in Fig. 6 (c). Thus, one can achieve the quasi-static stop band and the negative group velocity simultaneously, which is favorable in various metamaterial system design. The transition from the positive to the negative group velocity by adjusting the length of the block is shown in the Supplementary Material.

## Conclusions

In this paper, we proposed a new idea of the zero rotational stiffness to achieve extremely low frequency stop band, which may cover even the quasi-static frequency range. Unlike the previous ways to achieve the stop band, such as the Bragg scattering or the internal resonance phenomena, we made the first dispersion branch of the flexural wave mode as a flat curve at the zero frequency by using the idea of the zero rotational stiffness. In this way, we found that one can achieve very broad stop band from almost zero frequency to the cutoff frequency of the second dispersion branch. As a possible way to achieve the zero rotational stiffness, we developed the new metamaterial consisting of the blocks and the links. Here, the hinge connection was introduced between the blocks and the links to achieve the zero rotational stiffness. Analytic and numerical investigations showed that the developed elastic metamaterial can exhibit very broad stop band at extremely low frequency range. In addition, further investigation showed that the metamaterial can have the negative group velocity in addition to the quasi-static stop band. Considering the importance of the extremely low frequency stop band in vibration and acoustic devices, the results in this paper is expected to open a new way in various vibration metamaterial devices. The actual realization of the metamaterial with the structure stabilization will be a subject of next researches.

## Acknowledgements

This work was supported by the Center for Advanced Meta-Materials (CAMM) funded by the
Ministry of Science, ICT and Future Planning as Global Frontier Project" (CAMM-2014M3A6B3063711)

## Author Contributions

J. H. Oh proposed the idea and performed the analytic and numerical investigations. B. Assouar contributed to the analysis of the results and supervised the research. J. H. Oh and B. Assouar wrote the manuscript.

## Additional Information

**Competing Financial Interests** The authors declare no competing financial interests.

## References

1. Sigalas, M. & Economou, E. N. Band structure of elastic waves in two dimensional systems. *Solid State Commun.* **86**, 141-143 (1993).

2. Kushwaha, M. S., Halevi, P., Dobrzynski, L. & Djafari-Rouhani, B. Acoustic band structure of periodic elastic composites. *Phys. Rev. Lett.* **71**, 2022-2025 (1993).

3. Yablonovitch, E. Inhibited spontaneous emission in solid-state physics and electronics. *Phys. Rev. Lett.* **58**, 2059-2062 (1987).

4. Brillouin, L. *Wave propagation in periodic structures*. (Dover Publications, Inc., 1946).

5. Liu, Z., Zhang, X., Mao, Y., Zhu, Y. Y., Yang, Z., Chan, C. T. & Sheng, P. Locally resonant sonic materials, *Science* **8**, 1734-1736 (2000).

6. Huang, H. H., Sun, C. T. & Huang, G. L. On the negative effective mass density in acoustic metamaterials. *Int. J. Engng. Sci.* **47**, 610-617 (2009).

7. Zhu, R., Liu, X. N., Huang, G. L., Huang, H. H. & Sun, C. T. Microstructural design
and experimental validation of elastic metamaterial plates with anisotropic mass
density, *Phys. Rev. B* **86**, 144307 (2012).

8 Hou, Z. & Assouar, B. M. Tunable solid acoustic metamaterial with negative elastic modulus. *Appl. Phys. Lett.* **106**, 251901 (2015).

9 Zhu, R., Liu, X. N., Hu, G. K., Sun, C. T. & Huang, G. L. A chiral elastic metamaterial beam for broadband vibration suppression. *J. Sound Vib.* **333**, 2759-2773 (2014).

10 Zhu, R., Liu, X. N., Hu, G. K., Sun, C. T. & Huang, G. L. Negative refraction of elastic waves at the deep-subwavelength scale in a single-phase metamaterial. *Nat. Commun.* **5**, 5510 (2014).

11 Romero-Garcia, V., Krynkin, A., Garcia-Raffi, L. M., Umnova, O. & Sanchez-Perez, J. V. Multi-resonant scatterers in sonic crystals: Locally multi-resonant acoustic metamaterial, *J. Sound Vib.* **332**, 184–198 (2013).

12 Oh, J. H., Kwon, Y. E., Lee, H. J. & Kim, Y. Y. Elastic metamaterials for independent realization of negativity in density and stiffness. *Sci. Rep.* **6**, 23630 (2016).

13 Oh, J. H., Seung, H. M. & Kim, Y. Y. Adjoining of negative stiffness and negative density bands in an elastic metamaterial. *Appl. Phys. Lett.* **108**, 093501 (2016).

14 Oudich, M., Senesi, M., Assouar, B. M., Ruzenne, M., Sun, J.-H., Vincent, B., Hou, Z. & Wu, T.-T. Experimental evidence of locally resonant sonic band gap in two-dimensional phononic stubbed plates. *Phys. Rev. B* **84**, 165136 (2011).

15 Williams, E. G., Roux, P., Rupin, M. & Kuperman, W. A. Theory of multiresonant metamaterials for A₀ Lamb waves, *Phys. Rev. B* **91**, 104307 (2015).

16 Assouar, B. M. & Oudich, M. Enlargement of a locally resonant sonic band gap by using double-sides stubbed phononic plates, *Appl. Phys. Lett.* **100**, 123506 (2012).

17 Assouar, B. M., Sun, J.-H., Lin, F.-S. & Hsu, J.-C. Hybrid phononic crystal plates
for lowering and widening acoustic band gaps. *Ultrasonics* **54**, 2159–2164 (2014).

18 Li, Y & Assouar, B. M. Acoustic metasurface-based perfect absorber with deep subwavelength thickness. *Appl. Phys. Lett.* **108**, 063502 (2016).

19 Oh, J. H., Seung, H. M. & Kim, Y. Y. A truly hyperbolic elastic metamaterial lens. *Appl. Phys. Lett.* **104**, 073503 (2014).

20 Lee, S. H., Park, C. M., Seo, Y. M., Wang, Z. G. & Kim, C. K. Acoustic metamaterial with negative density. *Phys. Lett. A* **373**, 4464-4469 (2009).

21 Lee, S. H., Park, C. M., Seo, Y. M., Wang, Z. G. & Kim, C. K. Acoustic metamaterial with negative modulus. *J. Phys.: Condens. Matter* **21**, 175704 (2009).

22 Yao, S., Zhou, X. & Hu, G. Investigation of the negative-mass behaviors occurring below a cut-off frequency. *New J. Phys.* **12**, 103025 (2010).

23 Yu, D., Wen, J., Shen, H., Xiao, Y. & Wen, X. Propagation of flexural wave in periodic beam on elastic foundations. *Phys. Lett. A* **376**, 626-630 (2012).

24 Graff, K. F. *Wave Motion in Elastic Solids* (Dover, New York, 1991).

25 Meirovitch, L. *Fundamentals of Vibrations* (McGraw-Hill, 2001).

26 Langlet, P., Hladky-Hennion, A.-C. & Decarpigny, J.-N. Analysis of the propagation of plane acoustic waves in passive periodic materials using the finite element method. *J. Acoust. Soc. Am.* **98**, 2792-2800 (1995).

## Figure Legends

Fig. 1. (Color online) (a) Conversion of general beam under flexural wave motion to the mass-spring system, (b) general mass-spring system for flexural wave, (c) typical wave dispersion curve of flexural wave in general beam, (d) wave dispersion curve of flexural wave with the zero rotational stiffness.
Fig. 2. (Color online) The proposed unit cell structure having zero rotational stiffness.

Fig. 3. (Color online) (a) Conversion of the proposed metamaterial to the mass-spring system, (b) the equivalent mass-spring system of the proposed metamaterial.

Fig. 4. (Color online) (a) Continuum model of the proposed metamaterial’s unit cell, (b) the hinge connection condition imposed in the current numerical simulation.

Fig. 5. (Color online) (a) Plot of the analytically and numerically calculated wave dispersion curves, (b) the mode shapes of the unit cell at the point A₁ (left) and A₂ (right) in (a), respectively. Here, the color represents the total displacement.

Fig. 6. (Color online) (a) Continuum model of the metamaterial system to measure the wave transmission, (b) the measured wave transmission for various number of the unit cell, total displacement plot and deformed shape of the metamaterial system at (c) 200 Hz and (d) 500 Hz.

Fig. 7. (Color online) (a) Continuum model of the metamaterial with the negative group velocity, (b) plot of the analytically and numerically calculated wave dispersion curves, (c) the mode shapes of the unit cell at the point A₃ (left) and A₄ (right) in (a), respectively. Here, the color represents the total displacement.

Fig. 8. (Color online) (a) Continuum model of the metamaterial system consisting of the metamaterial in Fig. 7 (a) to measure the wave transmission, (b) the measured wave transmission for various number of the unit cell, total displacement plot and deformed
shape of the metamaterial system at (c) 200 Hz and (d) 400 Hz.

<figure><img src="image_3.png" /><figcaption>Fig. 1.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>Fig. 2.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>Fig. 3.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 4.</figcaption></figure>
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>Fig. 5.</figcaption></figure>
<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /><figcaption>Fig. 6.</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>Fig. 7.</figcaption></figure>
<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

<figure><img src="image_21.png" /><figcaption>Fig. 8.</figcaption></figure>