# Auxetic-like metamaterials as novel earthquake protections

Bogdan Ungureanu, Younes Achaoui, Stefan Enoch, Stéphane Brulé, Sébastien Guenneau

► To cite this version:

Bogdan Ungureanu, Younes Achaoui, Stefan Enoch, Stéphane Brulé, Sébastien Guenneau. Auxetic-like metamaterials as novel earthquake protections. EPJ Applied Metamaterials, 2016, 2, pp.17. 10.1051/epjam/2016001. hal-01283651

HAL Id: hal-01283651

https://hal.science/hal-01283651

Submitted on 8 Nov 2018

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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Auxetic-like metamaterials as novel earthquake protections

Bogdan Ungureanu¹,²,ª, Younes Achaoui²,ª, Stefan Enoch², Stéphane Brûlé³, and Sébastien Guenneau²,\*

¹ Faculty of Civil Engineering and Building Services Technical University “Gheorghe Asachi” of Iasi, 43, Dimitrie Mangeron Blvd., 700050 Iasi, Romania

² Aix-Marseille Université, CNRS, Centrale Marseille, Institut Fresnel UMR7249, 13013 Marseille, France

³ Dynamic Soil Laboratory, Ménard, 91620 Nozay, France

**Abstract** – We propose that wave propagation through a class of mechanical metamaterials opens unprecedented avenues in seismic wave protection based on spectral properties of auxetic-like metamaterials. The elastic parameters of these metamaterials like the bulk and shear moduli, the mass density, and even the Poisson ratio, can exhibit negative values in elastic stop bands. We show here that the propagation of seismic waves with frequencies ranging from 1 Hz to 40 Hz can be influenced by a decimeter scale version of auxetic-like metamaterials buried in the soil, with the combined effects of impedance mismatch, local resonances and Bragg stop bands. More precisely, we numerically examine and illustrate the markedly different behaviors between the propagation of seismic waves through a homogeneous isotropic elastic medium (concrete) and an auxetic-like metamaterial plate consisting of 4³ cells (40 m × 40 m × 40 m), utilized here as a foundation of a building one would like to protect from seismic site effects. This novel class of seismic metamaterials opens band gaps at frequencies compatible with seismic waves when they are designed appropriately, what makes them interesting candidates for seismic isolation structures.

**Key words:** Stop bands, Auxetics, Mechanical metamaterials, Seismic waves.

## 1. Introduction

Auxetics, a term coined by Ken Evans, are composites that become thicker perpendicular to the applied force when stretched, what leads to negative Poisson ratios [1, 2]. Mechanical metamaterials [3, 4] are periodic structures, counterpart of electromagnetic metamaterials, consisting of materials with high contrast in mechanical properties, which have been studied at small scales ranging from micrometers to centimeters. Control of surface seismic waves has been experimentally demonstrated in soils structured at the metric scale [5] that can be viewed as analogs of phononic crystals with holes [6]. In order to achieve control of seismic waves in the subwavelength regime, it seems natural to look at large scale analogs of mechanical [3, 4] and acoustic [7–16] metamaterials, which are periodic structures that can manipulate acoustic (e.g. pressure) and elastic (e.g. surface Lamb, Rayleigh, or bulk coupled shear and pressure) waves. The Veselago-Pendry flat lens via negative refraction is perhaps the most famous paradigm of electromagnetic [17–19], acoustic [20] and platonic [21] crystals. An interesting application that arises from the

periodic distribution of boreholes or inclusions embedded in a soil is as aforementioned the seismic wave shielding for Rayleigh wave frequencies within stop bands [5] (frequency intervals where, under certain conditions, propagation of elastic waves, or some of their polarizations, is forbidden), but also the flat seismic lens [22] for Rayleigh waves. Interestingly, forest of trees can also serve as seismic metamaterials for Rayleigh waves [23]. Such experiments on control of Rayleigh waves can find some applications in protection of urban infrastructures against earthquakes in soft sedimentary soils [24, 25].

However, one would like to design seismic metamaterials that have the ability to create band gaps not only for surface Rayleigh waves [5, 22, 23] but also for all other elastic wave polarizations: within certain frequency ranges known as complete stop bands, an incoming mechanical wave would be completely reflected by the structure, whether the seismic wave signal propagates near the air-soil interface (Rayleigh or Love waves) or within the soil (coupled shear and pressure mechanical waves). Notably, it is an interesting challenge to achieve the base isolation of a structure by applying the concept of metamaterials in civil engineering. In what follows, we would like to demonstrate that base isolation with auxetic-like

<footnote>ª Both authors contributed equally.</footnote>

<footnote>\*e-mail: sebastien.guenneau@fresnel.fr</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Illustration of the basic bow-tie element (known to lead to negative to Poisson's ratios in the static limit [4]), which is designed here to achieve tunable elastic stop bands for seismic waves (the periodic cell is 10 m wide).</figcaption></figure>

metamaterials could be an effective way to improve the seismic response of a building due to an earthquake in contrast to traditional seismic design methods which aimed to increase the strength of the structural parts of the building.

The Poisson ratio $v$ is a positive parameter for many common isotropic elastic natural materials, with values ranging from 0 (e.g. cork) to 0.5 (e.g. rubber), which means that they expand laterally upon pushing on them axially. Man-made auxetics were introduced when Roderik Lakes studied three-dimensional foams with $v < 0$ [1] in the mid eighties. More recently, there has been a keen interest in rationally designing auxetic mechanical metamaterials, which can additionally be intentionally anisotropic. However, auxetics need not be necessarily engineered, e.g. living bone tissue is a natural anisotropic auxetic material [2]. All fabricated rationally designed auxetic metamaterials are based on a few basic motifs [3]. One of the simplest motifs is the bow-tie element, which we shall consider in our study. Upon pushing (resp. compressing) along the vertical $z$ direction, a usual (i.e. not auxetic) material contracts (resp. expands) along the horizontal $x$ and $y$ directions. For a negative $v$, a compression along the vertical axis $z$ leads to a contraction of the shown bow-tie motif in Figure 1 along the horizontal $x$ and $y$ direction [4].

In Figure 1, the sign and magnitude of $v$ can be controlled via the angle $\alpha$ [4]. A zero Poisson ratio $v$ is expected to occur near $\alpha = \pi/2$. The Poisson ratio effectively changes if the modulus of the strain is not small compared to unity, i.e., if one leaves the linear mechanical regime and enters the nonlinear regime, in which the change in the angle $\alpha$ is no longer negligible [3]. This basic bow-tie motif can be assembled into two-dimensional model systems and into three-dimensional (anisotropic) mechanical metamaterials [4].

The angle $\alpha$ and hence the Poisson ratio $v$ can be adjusted in the fabrication process, as was demonstrated by the group of Martin Wegener at Karlsruhe Institute for Technology [4]. The quoted Poisson ratios $v$ measured by this group (more precisely, the components $v_{ij}$ of the Poisson matrix since the

composite medium has an anisotropic elastic response) can be tailored to positive, near zero or even negative values. In the present paper, we propose to use such auxetic-like composites in order to design elastic stop band metamaterials for seismic waves propagating in sedimentary soils at hertz frequencies [5].

## 2. Stop band properties of auxetic-like metamaterials

We select concrete structures that should be fairly simple to fabricate and with elastic and geometric parameters that permit seismic waves molding. For the calculation of the band structure, the finite element method was used (Comsol Multiphysics). In the computations of Figure 2, only one basic cell was considered with quasi-periodic boundary conditions enforced along each direction, in accordance with the Floquet-Bloch theorem. One usually starts with a vanishing Bloch vector $k = (0, 0, 0)$, so that the resulting frequency spectrum consists of a well-defined countable set of eigenfrequencies associated with periodic eigenfields (stationary waves) sitting on local extrema of the band structure (sometimes at band gap edges). Changing the $k$ vector value used in the quasi-periodic boundary conditions, one finds the eigenfrequencies and associated eigenfields that describe the overall band structure. In our calculations, we considered the irreducible Brillouin zone $\Gamma XMR$, where $\Gamma = (0, 0, 0)$, $X = (2\pi/a, 0, 0)$, $M = (2\pi/a, 2\pi/a, 0)$ and $R = (2\pi/a, 2\pi/a, 2\pi/a)$, with $a$ the array pitch. We realize that one should in principle describe all points within the irreducible zone (which is the volume of a tetrahedron), but it is customary to only look at its edges (invoking some symmetry reasons) for the sake of saving computational time.

The dimensions used for the computations of band structures associated with the bow-tie motifs shown in Figures 2a–2d are cells $10 \times 10 \times 10 \text{ m}^3$ (same length along the $x$, $y$ and $z$ directions). Each cell is connected to its neighboring cells within the lattice via 6 bars are $1 \times 1 \text{ m}^2$ in cross-section. The cross-section of each bar constituting the bow-tie element in the basic cell is also $1 \times 1 \text{ m}^2$.

The four cases studied in this work consist of bars made from concrete and the only thing that changes between the cases is the angle: $\alpha = 0^\circ$; $\alpha = 30^\circ$; $\alpha = 45^\circ$; $\alpha = 60^\circ$; Note that we use the density $\rho = 2.3 \text{ g/cm}^3$, the longitudinal (compres sional) $c_l = 3.475 \text{ km/s}$ and the transverse (shear) $c_s = 2.022 \text{ km/s}$ components of elastic wave velocity for con crete. This corresponds to a Young modulus $E = 30 \text{ GPa}$ and a Poisson ratio $v = 0.3$. Assuming that the lattice constant $a = 10 \text{ m}$, frequencies are expected to be in the order of a few Hertz.

The first structure that was numerically examined is for an angle $\alpha = 0^\circ$ and is shown in Figure 2a. In this case one can observe on the dispersion diagram that there are no stop bands in the frequency range 0–60 Hz. Stop bands actually occur at frequencies less relevant for seismic wave protection, i.e., above 50 Hz [5].

The second structure (for $\alpha = 30^\circ$, see Figure 2b), which is in fact the first bow-tie auxetic cell, gives rise to a tiny total band gap in a frequency range of interest for seismic wave protection, around 40 Hz. A second, larger, stop band lies
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Band diagrams of auxetic-like metamaterials generated with elementary cells (10 m × 10 m × 10 m) with bars (concrete) making an angle (a) α = 0°; (b) α = 30°; (c) α = 45°; (d) α = 60° where the horizontal axis is the normalized k wavevector describing the edges of the irreducible Brillouin zone ΓXMR and the vertical axis is the frequency (Hz). Green curves correspond to shear and pressure waves propagating in homogeneous isotropic bulk (concrete). Shaded regions mark the location of complete (dark gray) and partial (light gray) stop bands.</figcaption></figure>

above the frequency range of interest. This design could be used as an alternative to models of seismic protection proposed in references [27, 28].

The third structure (for α = 45°, see Figure 2c), displays a wide total stop band that opens up around 40 Hz, from 33 Hz to 47 Hz (14 Hz of bandwidth, so with a relative bandwidth of 35%).

By comparison, the fourth structure (for α = 60°, see Figure 2d) displays wide complete in-plane stop bands, which shrink along the MR and RΓ: This is easily understood by the fact that the vertical bars now meet at the center of the cell, making homogeneous vertical bars along which waves can travel at any frequency.

It then transpires from panels (a)–(d) that the auxetic-like structure with an angle of α = 45° has the most promising band structure. Also worth noticing is the markedly different anisotropy of auxetic-like metamaterials in panels (a)–(d), wherein the second (resp. third) acoustic band displays decreasing (resp. increasing) group velocity with increasing angle α. This might find some applications in seismic wave control such as cloaking.

### 3. Interpretation of stop bands with quasi-static and dynamic effective properties

There are interesting connections between models of metamaterials using low-frequency high-contrast and high-frequency low-contrast homogenization theories (see for instance Chapter 1 in Ref. [29] where simple scalar wave equations are studied), in that the two routes lead to frequency dependent material parameters. Low-frequency high-contrast homogenization of the Navier equation is known to lead to a

frequency dependent rank-2 tensor of effective density [30], but in what follows we would like to give some elements of proof for similar effective behavior of Young's modulus and Poisson's ratio in the high-frequency regime. Let us start with some well-known facts of homogenization theory in the low-frequency low-contrast case. According to Christensen and Lo [31], and Nemat-Nasser and Willis [32] the effective density $\langle\rho\rangle$ and Poisson's ratio squared $\langle v^2\rangle$ can be approximated as follows in the case of a cubic array of ellipsoidal inclusions

$$ \langle v^2 \rangle = \frac{\omega^2 a^2 \langle \rho \rangle}{\langle C_{1111} \rangle}, \quad \langle \rho \rangle = \rho^{(1)} f_1 + \rho^{(2)} f_2, \quad (1) $$

$$ \langle C_{1111} \rangle = \bar{C}_{1111}^{(1)} f_1 + \bar{C}_{1111}^{(2)} f_2, \\ f_1 = 1 - f_2, \quad f_2 = \frac{\pi b_1 b_2 b_3}{6a^3}, \quad (2) $$

where $\rho^{(1)}$ and $\rho^{(2)}$ are the densities of media 1 (bulk medium) and 2 (ellipsoidal inclusion), $\bar{C}_{1111}^{(1)}$ and $\bar{C}_{1111}^{(2)}$ are the first component of the elasticity tensor of media 1 and 2, respectively. Moreover, $f_1, f_2$ denotes the volume fraction of media 1 and 2, a is the side length of the periodic cell, and $b_1, b_2, b_3$ are the semi-axes of the ellipsoidal inclusion. If we use these formulas to get a first estimate of the effective property of the structured medium of Figure 1 in the quasi-static regime (leading order approximation), assuming that medium 1 is the volume occupied by the bars and medium 2 is the outer (air) medium, we cannot find any unusual behavior, as testified by the lowest branches in the dispersion diagrams in Figures 2a–2d. The fact that the pressure (third lowest) band on one hand and the shear (first and second lowest) bands on
the other hand get increasingly different wavespeeds is due to the fact that in the quasi-static limit the auxetic-type metamaterial behaves effectively like a fluid.

However, the long wavelength limit breaks down at higher frequencies (or when the contrast gets bigger), where one needs to invoke dynamic effective properties that take into account the fine structure (here, bow-tie geometry) of the structured media, as described in references [32–36]. Notably, in references [35, 36] it is clearly established that both the effective elasticity tensor $C_{\text{eff}}$ and the effective density $\rho_{\text{eff}}$ are frequency dependent, and [36] exhibits simultaneously negative values of entries of $C_{\text{eff}}$ and $\rho_{\text{eff}}$ in stop bands. This suggests negative effective Poisson's ratios might occur in many periodic structures in the stop bands, not necessarily with auxetic-type media such as first introduced by Lakes [1], but also in many phononic crystals and mechanical metamaterials.

Indeed, if we assume that we are near resonances, then homogenization requires dynamic effective parameters that were first identified by Auriault in 1983 in the context of diffusion processes in high-contrast periodic media [37] using multiple scale expansion techniques, with subsequent extension to elastodynamics [30, 38]. At the turn of the millennium, the mathematician Zikhov proposed a rigorous mathematical framework using two-scale convergence techniques in weak forms of scalar partial differential equations with fast oscillating parameters [39], that can be extended to elastodynamic equations [40]. Similar models of acoustic metamaterials exist using asymptotic analysis of Bessel functions [8], multi-structures [9], periodic unfolding techniques [10]. The link between the Poisson ratio and the appearance of stop bands in auxetic metamaterials, was pointed out in reference [41].

In our case, we consider a one-phase medium with stress-free boundary conditions like in reference [41] but we carry out a more in-depth analysis by noticing that near resonances, modes get localized inside the bars of the structured medium in a way similar to what they would within a two-phase high-contrast medium ($\bar{C}_{1111}^{(2)}/\bar{C}_{1111}^{(1)} = O(\eta^2)$), where the stiff phase $\bar{C}_{1111}^{(1)}$ occupies the volume of the bars and the soft phase $\bar{C}_{1111}^{(2)}$ is the vacuum between the bars. This configuration of a composite with connected high-rigidity solid and soft inclusions has been studied by Auriault and Boutin [30], using the usual ansatz in the displacement field $u = u^0(x, y) + \eta u^1(x, y) + \eta^2 u^2(x, y) + \dots$ where $u^i(x, y) = u^i(x_1, x_2, x_3, y_1, y_2, y_3)$ are periodic fields in the microscopic variable $y$. They found that the effective medium is described by the following homogenized Navier equation:

$$ \mathrm{div}_x(\sigma_{\mathrm{eff}}^x(u_{\mathrm{eff}})) + \omega^2 \rho_{\mathrm{eff}}^x(\omega)u_{\mathrm{eff}} = 0, \quad (3) $$

where $u_{\text{eff}}(x) = \int u^0(x, y)dy := \iiint (u_1^0, u_2^0, u_3^0)(x_1, x_2, x_3, y_1, y_2, y_3)dy_1dy_2dy_3$ is the homogenized displacement field, $\sigma_{\text{eff}}^x$ and $\rho_{\text{eff}}^x$ are the rank-2 effective stress and density tensors, respectively (both of which only depend upon the macroscopic variable $x$).

The effective (anisotropic) density is given by the following Drude-like expression:

$$ \rho_{\text{eff}}(\omega) = \langle\rho\rangle I + \rho \sum_{i=1}^{\infty} \frac{\langle\phi^i\rangle \otimes \langle\phi^i\rangle}{\langle\int \phi^i \phi^i dy\rangle} \frac{1}{\omega_i^2/\omega^2 - 1}, \quad (4) $$

with $\phi^i(y) = (\phi_1^i, \phi_2^i, \phi_3^i)(y_1, y_2, y_3)$ a vector valued eigenfunction of the eigenvalue problem

$$ \begin{aligned} \mathrm{div}_y(C^y : \varepsilon^y(\phi^i)) &= -\lambda_i \phi^i \text{ in the soft phase and } \phi^i \\ &= 0 \text{ on its boundary,} \end{aligned} \quad (5) $$

which is set on the microscopic periodic cell, which from renormalization is the unit cube $[0, 1]^3$. It is well known that for this spectral problem there is a countable set of eigenvalues $0 \le \lambda_1 \le \lambda_2 \le \lambda_3 \le \dots < +\infty$ with associated eigenfunctions $\phi_i^i$, $i = 1, 2, 3, \dots$ generating an orthogonal basis in the space of finite energy functions. When the frequency squared $\omega^2$ gets close to $\omega_i^2 = \lambda_i/\rho$ (4) takes negative values, what can be used as an interpretation of stop bands as seen in Figures 2b–2d. A negative effective mass density could be notably invoked in stop bands of [42–45], but this would come short to a complete explanation of band structures.

Besides from the frequency dependent effective density in equation (3), the effective stress tensor is related to an effective elasticity tensor through $\sigma_{\text{eff}} = C_{\text{eff}} : \varepsilon^x$ where the tensor of deformation has the form $\varepsilon_{ij}^x = \frac{1}{E} \left[ \sigma_{ij}(1 + v) - v\delta_{ij} \sum_{k=1}^3 \sigma_{kk} \right]$

with $\delta_{ij}$ the Kronecker symbol, E the Young modulus and v the Poisson ratio. The effective elasticity tensor is computed from an annex problem of elasto-static type set on the unit cell excluding the soft phase by noticing that the stiff phase moves like a rigid body (translation) at the leading order of elastic field displacement. The effective elasticity tensor essentially contains the anisotropic features of the stiff phase, and it does not seem to have much to do with appearance of stop bands i.e. it is not frequency dependent. However, as recognized in reference [36], there is no uniqueness in the description of the homogenized medium at high frequencies, and alternative effective parameters can be achieved by assuming other types of limits in phase contrasts when one leaves the quasi-static regime. This is exactly what has been found by Craster and coauthors in 2010 [33] when introduced the concept of high-frequency homogenization (HFH). In their seminal paper, these authors take the usual ansatz for the displacement field, but instead of considering a high contrast, they rescale the frequency and they are led to frequency dependent effective parameters for both stress and density tensors in the homogenized Navier equations [35]. Regarding HFH of thin plates, Antonakakis and Craster find frequency dependent effective rigidity and density in the Kirchoff-Love equations [34], which is consistent with the work of Torrent et al. that demonstrates negative effective Young modulus, density and Poisson ratio using scattering matrix asymptotics in thin periodic plates with soft inclusions [46]. We thus claim that stop bands in Figures 2b–2d can be interpreted in terms of negative effective density, Young's modulus and Poisson ratio, depending upon the homogenization approach used. Classical work on auxetic materials [47–52] and mechanical metamaterials [53] as well as on homogenization [54–57] do not seem to have looked into the dynamic effective properties of auxetic-type materials, so we hope our section will foster theoretical and experimental studies in this direction.
<figure><img src="image_4.png" /><figcaption>**Figure 3.** (a) Band diagram for a concrete building (20 m × 20 m × 80 m) atop a concrete homogeneous plate (40 m × 40 m × 40 m); Floquet-Bloch conditions are set on vertical sides of the plate and stress-free conditions hold elsewhere; Flat bands correspond to eigenmodes of the building around 1 Hz, 9 Hz, 17 Hz, etc. that couple to the flexural band (lowest band at $\Gamma$). (b) Representative eigenmode of (a) at 17 Hz. (c) Eigenmode with suppressed building's vibration around 17 Hz when the plate is structured with auxetic metamaterial (4³ elementary cells with bars making an angle as in Figure 2c).</figcaption></figure>

**4. Quantification of energy loss through auxetic-like metamaterials**

As we discussed in the previous section, there are further dynamic effects induced by the effective elasticity tensor $C_{\text{eff}}$ that lead to other negative effective parameters, including the bulk modulus and Poisson ratio, as shown by Torrent et al. in the context of a thin plate theory [46] that seems particularly well suited to the configuration of an auxetic-like metamaterial plate shown in Figure 3. The detailed asymptotic analysis required for the identification of the frequency dependent effective elasticity tensor goes beyond the scope of the present paper but, one can clearly see the dramatic changes in the higher part of the dispersion diagrams in Figures 2a–2d when the angle between bars changes, which is caused by a dynamic effective tensor changing radically of nature. Besides from that, there are marked local resonances in Figure 3a that can be interpreted with a negative effective Poisson ratio according to the reference [46]. One should note in passing that the fact that there is a dynamic negative effective density does not necessarily mean that there is a negative effective Poisson ratio in 3D periodic structures, but we find that for plates the dynamic effective density, rigidity and Poisson ratio have the form of equation (4), which is consistent with [46], provided that the wave wavelength is much larger than the plate thickness. Furthermore, one should note that auxetic-like metamaterials exhibit slow modes in Figure 2 (dispersion of pressure and shear waves propagating in the homogeneous concrete bulk are marked by green lines), and this can be interpreted in terms of vanishing effective density, Young's modulus and Poisson's ratio in light of HFH.

Let us now focus our attention on slow (localized) modes responsible for devastating buildings' resonances. In order to

carry out their band structure analysis, we consider a model for a doubly periodic array of bars 80m high, with a square cross-section of 20 m × 20 m (that stand for buildings) lying atop a homogeneous concrete plate (40 m × 40 m × 40 m). We enforce Floquet-Bloch conditions on the vertical sides of the plate, with stress-free conditions elsewhere. The elastic band structure is shown in Figure 3a. Note that the periodicity of the cell is irrelevant in the frequency range of interest since the band folding occurs far above 60 Hz (about 200 Hz for compressional waves). Eigenfrequencies of the buildings appear as flat bands on the dispersion diagram. One can also see that the first band on Figure 3a has a quadratic behavior around the $\Gamma$ point, unlike in Figure 2, which is due to the fact that the buildings lie atop a 40 m deep concrete plate, so the corresponding elastic wave is a Lamb (flexural type) surface wave.

We display in Figure 3b a representative eigenmode at 17 Hz of a building which corresponds to the fifth flat band in Figure 3a. This vibration is completely suppressed when the plate is structured with an auxetic-like metamaterial in Figure 3c. The structured plate consists of a supercell made of 4³ elementary cells with bars making an angle $\alpha = 45^\circ$ like in Figure 2c with Floquet-Bloch conditions on vertical sides and stress-free conditions elsewhere. Upon inspection of the band structure corresponding to the configuration in Figure 3c – which is not shown as it has a large number of bands – we notice a stop band that opens around 17 Hz: this tiny stop band that is marked by two closely located local extrema in the transmission (black) curve near 20 Hz in Figure 4a, can be interpreted in terms of a negative effective density, rigidity and Poisson's ratio: indeed, the building's vibration creates a monopolar resonance (a defect mode) within the auxetic supercell array, so that the Drude-like effective formula (4) applies
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Quantification of energy loss (a) with red (resp. black) curves representing transmission through a building atop a concrete plate (resp. a concrete plate structured with an auxetic metamaterial); Magenta curve gives the total displacement field stored inside the building above an auxetic-like metamaterial normalized with respect to the same total displacement field above the homogeneous plate. Representative out-of-plane displacement without (b1), (c1) and (d1) and with (b2), (c2) and (d2) auxetic-like metamaterial illustrate how seismic protection works at 1 Hz (b1), (b2) with the impedance mismatch between the homogeneous plate and the auxetic-like metamaterial; at 9 Hz (c1), (c2) with a local resonance; at 37 Hz (d1), (d2) with the Bragg band gap.</figcaption></figure>

to all of these parameters, see reference [46] for similar effective properties of a periodic plate.

We observe in Figure 4, some low transmission from 0 to 15 Hz, which can be attributed to the liquid-like effective properties of the plate in the quasi-static limit (see Sect. 3). Indeed, Figure 2c displays anisotropic bands with small velocities compared to a homogeneous isotropic plate, and the same holds true from the band structure associated with the supercell consisting of $4^3$ cells like in Figure 2c. Physically speaking, by using auxetic-like elements we soften the medium underneath the building, which leads to a strong impedance mismatch between the auxetic-like building's foundation and the homogeneous plate medium. This results in a total reflection over a broad range of frequencies, according to the transmission loss of Figure 4, that completely suppress the vibrations of the building at 1 Hz and 9 Hz (marked by the first two local resonances in Figure 3a). The suppression of the building's resonance at 37 Hz is more conventional, since this falls within a Bragg stop band according to Figure 2c.

Let us now quantify the amount of elastic energy stored within the building and the transmission loss through this building, with and without the auxetic-like metamaterial. To do so, we apply periodic boundary conditions on the transverse sides of an elongated homogeneous concrete plate comprising a building (see Figures 4b1–4d2) so as to prevent appearance of transverse beam modes. A line source is applied at the left edge of the plate in order to simulate an incoming earthquake. The red (resp. black) curve in Figure 4a shows the transmission through the plate with a building on its own (resp. with the auxetic-like metamaterial). One can see that regardless of the frequency of the incoming wave, transmission through the building lying atop the auxetic-like metamaterial foundation is at least 10 dB lower than without the metamaterial. One should note that any building located behind the auxetic-like metamaterial would be protected. In order to quantify the level of protection of the building lying atop the auxetic-like metamaterial, we compute (magenta curve) the ratio of the norm of the total displacement field stored inside the

building above the plate by the norm of the total displacement field stored inside the building above the plate structured with auxetic-like metamaterial. The building's protection which is mainly achieved via wave velocity impedance mismatch is clearly demonstrated throughout the frequency range 0–100 Hz, with a slightly more pronounced protection in the range of stop band frequencies from 33–47 Hz.

Some typical examples of out-of-plane (vertical) displacement field are shown at 1 Hz with b1 (resp. without b2) auxetic-like metamaterial, at 9 Hz with c1 (resp. without c2) auxetic-like metamaterial and at 37 Hz with d1 (resp. without d2) auxetic-like metamaterial. The periodic foundation is clearly different from the traditional base isolation in which it causes a fundamental frequency shift in the structure, thus reducing its response and generating a frequency gap. Implementing auxetic-like cells in the foundations of sensitive buildings would significantly affect their static and dynamic responses during earthquakes.

**5. Concluding remarks and perspectives on seismic metamaterials**

The main conclusion of our study is that auxetic-like materials who have been widely studied for their special properties linked to a negative Poisson ratio $v$, also have interesting stop band properties that can be used in the context of seismic wave protection. The stop bands can be associated with frequency dependent mass density $\rho$, but also bulk modulus [12, 58], shear modulus, in a way similar to what was unveiled in acoustic metamaterials. Nonetheless, it is also possible to achieve frequency dependent Poisson's ratios in stop bands of certain mechanical metamaterials [36, 46], what is less well known. Theories of low-frequency high-contrast homogenization [37–40] and high-frequency homogenization [33–35] exist that allow to investigate such dynamic effective properties of auxetic metamaterials. Using arrays of auxetic cells for seismic structural protection is a new way to perform elastic band gaps
in order to prevent seismic wave propagation over specific frequency ranges, but we also discovered that the very strong impedance mismatch between elastic wave velocity within the homogeneous bulk of concrete and that within the auxetic metamaterial enables suppression of wave transmission over a very large frequency range (0–100 Hz), which opens interesting avenues in earthquake protection. We stress that we have performed our study in elastic conditions, but we are well aware that soils display visco-elastic properties causing frequency dependent damping effects in soft sediments [59].

One could also envisage to combine auxetic-like seismic metamaterials within the soil with other types of protections such as forests of trees [23] to further widen the range of stop band frequencies for Rayleigh waves. It might be interesting to also investigate thick pillars above the soil, since previous numerical and experimental studies on small scale photonic crystals composed of nickel pillars grown on a lithium niobate substrate [14, 15] have shown complete stop bands at low frequencies which are also robust versus disorder. Such seismic metamaterials would offer alternative ways to protect buildings, which necessitate specific modal analysis [24, 59] by civil engineers in order to avoid disastrous seismic site effects.

*Acknowledgements.* Y. Achaoui, B. Ungureanu and S. Guenneau acknowledge funding from European Research Council (ERC Grant ANAMORPHISM).

## References

1. R. Lakes, Foam structures with a negative Poisson's ratio, Science 235 (1987) 1038-1040.

2. R. Shahar, P. Zaslansky, M. Barak, A.A. Friesem, J.D. Currey, S. Weiner, Anisotropic Poisson's ratio and compression modulus of cortical bone determined by speckle interferometry, J. Biomech 40 (2007) 252-264.

3. M. Kadic, T. Bückmann, R. Schittny, M. Wegener, Metamaterials beyond electromagnetism, Rep. Prog. Phys. 76 (2013) 126501.

4. T. Bückmann, N. Stenger, M. Kadic, J. Kaschke, A. Frölich, T. Kennerknecht, C. Eberl, M. Thiel, M. Wegener, Tailored 3D mechanical metamaterials made by dip-in direct-laser-writing optical lithography, Adv. Mater. 24 (2012) 2710.

5. S. Brülé, E. Javelaud, S. Enoch, S. Guenneau, Experiments on seismic metamaterials: molding surface waves, Phys. Rev. Lett. 112 (2014) 133901.

6. S. Benchabane, O. Gaiffe, G. Ulliac, R. Salut, Y. Achaoui, V. Laude, Observation of surface-guided waves in holey hypersonic phononic crystal, Appl. Phys. Lett. 98 (2011) 171908.

7. Z. Liu, X. Zhang, Y. Mao, Y.Y. Zhu, Z. Yang, C.T. Chan, P. Sheng, Locally resonant sonic materials, Science 289 (2000) 1734.

8. J. Li, C.T. Chan, Double-negative acoustic metamaterials, Phys. Rev. E 70 (2004) 055602.

9. A.B. Movchan, S. Guenneau, Split-ring resonators and localized modes, Phys. Rev. B 70 (2004) 125116.

10. A. Ávila, G. Griso, B. Miara, Bandes phoniques interdites en élasticité linéarisée, C.R. Acad. Sci. Paris: Ser. I 340 (2005) 933-938.

11. Z. Liu, C.T. Chan, P. Sheng, Analytic model of phononic crystals with local resonance, Phys. Rev. B 71 (2005) 014103.

12. N. Fang, D.J. Xi, J.Y. Xu, M. Ambrati, W. Spriturvanich, C. Sun, X. Zhang, Ultrasonic metamaterials with negative modulus, Nat. Mater. 5 (2006) 452.

13. G.W. Milton, J.R. Willis, On modifications of Newton's second law and linear continuum elastodynamics, Proc. R. Soc. A 463 (2007) 855-880.

14. A. Khelif, Y. Achaoui, S. Benchabane, V. Laude, B. Aoubiza, Locally resonant surface acoustic wave band gaps in a two-dimensional phononic crystal of pillars on a surface, Phys. Rev. B 81 (2010) 214303.

15. Y. Achaoui, V. Laude, S. Benchabane, A. Khelif, Local resonances in phononic crystals and in random arrangements of pillars on a surface, J. Appl. Phys. 114 (2013) 104503.

16. A. Colombi, P. Roux, M. Rupin, Sub-wavelength energy trapping of elastic waves in a meta-material, J. Acoust. Soc. Am. 136 (2014) EL192-EL196.

17. V.G. Veselago, The electrodynamics of substances with simultaneously negative values of $ε$ and $μ$, Sov. Phys. Usp. 10 (1968) 509-514.

18. J.B. Pendry, Negative refraction makes a perfect lens, Phys. Rev. Lett. 85 (2000) 3966-3969.

19. R.A. Shelby, D.R. Smith, S. Shultz, Experimental verification of a negative index of refraction, Science 292 (2001) 77-79.

20. S. Yang, J.H. Page, L. Zhengyou, M.L. Cowan, C.T. Chan, P. Sheng, Focusing of sound in a 3D phononic crystal, Phys. Rev. Lett. 93 (2004) 024301.

21. M. Dubois, E. Bossy, S. Enoch, S. Guenneau, G. Lerosey, P. Sebbah, Time-driven super-oscillations with negative refraction, Phys. Rev. Lett. 114 (2013) 013902.

22. S. Brulé, S. Enoch, S. Guenneau, Flat seismic lens, arXiv:1602.04492, 2016 (http://arxiv.org/abs/1602.04492).

23. A. Colombi, P. Roux, S. Guenneau, P. Gueguen, R.V. Craster, Scientific Reports 6 (2016) 19238, DOI: 10.1038/srep19238

24. S. Brulé, E. Javelaud, M. Marchand, Chimneys health monitoring during a nearby heavy dynamic compaction site, in Journées Nationales de Géotechnique et de Géologie de l'Ingénieur JNGG2012, 4-6 July 2012, Bordeaux, France, Tome II, pp, 919-926, 2012.

25. T. Gmur, Dynamique des structures, analyse modale numérique, Presses Polytechniques et Universitaires Romandes, Lausanne, Switzerland, 2008.

26. S.H. Kim, M.P. Das, Seismic waveguide of metamaterials, Mod. Phys. Lett. B 26 (2012) 1250105.

27. Z. Shi, Z. Cheng, H. Xiang, Seismic isolation foundations with effective attenuation zones, Soil Dyn. Earthquake Eng. 57 (2014) 143-151.

28. S. Krodel, N. Thome, C. Daraio, Wide band-gap seismic metastructures, Ex. Mech. Letters 4 (2015), 111-117, DOI: 10.1016/j.eml.2015.05.004.

29. R.V. Craster, S. Guenneau, Acoustic metamaterials, Springer Verlag, London, 2012.

30. J.L. Auriault, C. Boutin, Long wavelength inner-resonance cut-off frequencies in elastic composite materials, Int. J. Solids Struct. 49 (2012) 3269-3281.

31. R.M. Christensen, K.H. Lo, Solution for effective shear properties in three phase sphere and cylinder models, J. Mech. Phys. Solids 27 (1979) 315-330.
32. S. Nemat-Nasser, J.R. Willis, Homogenization of periodic elastic composites and locally resonant sonic materials, Phys. Rev. B. 83 (2011) 104103.

33. R.V. Craster, J. Kaplunov, A.V. Pichugin, High frequency homogenization for periodic media, Proc. R. Soc. Lond. A 466 (2010) 2341–2362.

34. T. Antonakakis, R.V. Craster, High frequency asymptotics for microstructured thin elastic plates and platonics, Proc. R. Soc. Lond. A 468 (2012) 1408–1427.

35. T. Antonakakis, R.V. Craster, S. Guenneau, Homogenisation for elastic photonic crystals and metamaterials, J. Mech. Phys. Solids 17 (2014) 84.

36. A. Srivastava, Elastic metamaterials and dynamic homogenization: a review, Int. J. Smart & Nano Mat. 6 (2015) 41–60.

37. J.-L. Auriault, Effective macroscopic description for heat conduction in periodic composites, Int. J. Heat Mass Transf. 26 (1983) 861–869.

38. J.-L. Auriault, G. Bonnet, Dynamique des composites élastiques périodiques, Arch. Mech. 37 (1985) 269–284.

39. V.V. Zhikov, On an extension of the method of two-scale convergence and its applications, Sbor. Math. 191 (2000) 973–1014.

40. V.P. Smyshlyaev, Propagation and localization of elastic waves in highly anisotropic periodic composites via two-scale homogenization, Mech. Mater. 41 (2009) 434–447.

41. X.X. Su, Y.F. Wang, Y.S. Wang, Effects of Poisson's ratio on the band gaps and defect states in two-dimensional vacuum: solid porous phononic crystals, Ultrasonics 52 (2012) 225–265.

42. O. Sigmund, J.S. Jensen, Systematic design of phononic band-gap materials and structures by topology optimization, Philos. Trans. R. Soc. Lond. Ser. A 361 (2003) 1001–1019.

43. Y. Achaoui, B. Ungureanu, S. Enoch, S. Brûlé, S. Guenneau, Seismic waves damping with arrays of inertial resonators, arXiv:1512.06078, 2015.

44. P. Wang, F. Casadei, S. Shan, JC Weaver, K Bertoldi, Harnessing buckling to design tunable locally resonant acoustic metamaterials, Phys. Rev. Lett. 113 (2014) 014301.

45. Y. Xiao, J. Wen, W. Wen, Longitudinal wave band gaps in metamaterial-based elastic rods containing multi-degree-of-freedom resonators, New J. Phys. 14 (2012) 033042.

46. D. Torrent, Y. Pennec, B. Djafari-Rouhani, Effective medium theory for elastic metamaterials in thin elastic plates, Phys. Rev. B 90 (2014) 104110.

47. J.N. Grima, L. Mizzi, K.M. Azzopardi, R. Gatt, Auxetic perforated mechanical metamaterials with randomly oriented cuts, Adv. Mater. 28 (2015) 385–389, DOI: 10.1002/adma.201503653

48. J.N. Grima, R. Caruana-Gauci, K.W. Wojciechowski, R. Gatt, Smart metamaterials with tunable auxetic and other properties, Smart Mater. Struc. 22 (2013) 084016.

49. J.N. Grima, R. Caruana-Gauci, Mechanical metamaterials: materials that push back, Nat. Mater. 11 (2012) 565–566.

50. Z.G. Nicolau, A.E. Motter, Mechanical metamaterials with negative compressibility transitions, Nat. Mater. 11 (2012) 608–613.

51. S. Babaee, J. Shim, J.C. Weaver, E.R. Chen, N. Patel, K. Bertoldi, 3D soft metamaterials with negative Poisson's ratio, Adv. Mater. 25 (2013) 5044–5049.

52. R. Lakes, Advances in negative Poisson's ratio materials, Adv. Mater. 5 (1993) 293–296.

53. J. Christensen, M. Kadic, O. Kraft, M. Wegener, Vibrant times for mechanical metamaterials, MRS Communications 5 (2015) 453–462.

54. P.S. Theocaris, G.E. Stavroulakis, P.D. Panagiotopoulos, Negative Poisson's ratio in composites with star-shaped inclusions: a numerical homogenization approach, Arch. Appl. Mech. 67 (1997) 274–286.

55. P.S. Theocaris, G.E. Stavroulakis, The homogenization method for the study of Poisson's ratio in fiber composites, Arch. Appl. Mech. 68 (1998) 281–295.

56. R.F. Almgren, An isotropic three-dimensional structure with Poisson's ratio = -1, J. Electricity 15 (1985) 427–430.

57. G.W. Milton, The theory of composites, Cambridge University Press, 2002.

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

59. R. Reitherman, Earthquakes and Engineers: An International History, ASCE Press, Reston, VA, 2012.