Issue 
EPJ Appl. Metamat.
Volume 6, 2019
Metamaterials'2018 – Microwave, mechanical, and acoustic metamaterials



Article Number  18  
Number of page(s)  12  
DOI  https://doi.org/10.1051/epjam/2019015  
Published online  20 September 2019 
https://doi.org/10.1051/epjam/2019015
Review
The influence of building interactions on seismic and elastic body waves
^{1}
Imperial College London, Department of Mathematics, London SW7 2AZ, UK
^{2}
Aix Marseille Univ, CNRS, Centrale Marseille, Institut Fresnel, 13013 Marseille, France
^{3}
FEMTOST, CNRS, Université de Bourgogne FrancheComté, 25044 Besançon, France
^{4}
Division of Computer, Electrical, and Mathematical Sciences and Engineering, King Abdullah University of Science and Technology (KAUST), Thuwal 2395569100, Saudi Arabia
^{5}
Indian Institute of Technology Bombay, Department of Mathematics, SW7 2AZ Mumbai, India
^{*} email: b.ungureanu@imperial.ac.uk
Received:
2
April
2019
Accepted:
23
July
2019
Published online: 20 September 2019
We outline some recent research advances on the control of elastic waves in thin and thick plates, that have occurred since the large scale experiment [S. Brûlé, Phys. Rev. Lett. 112, 133901 (2014)] that demonstrated significant interaction of surface seismic waves with holes structuring sedimentary soils at the meter scale. We further investigate the seismic wave trajectories of compressional body waves in soils structured with buildings. A significant substitution of soils by inclusions, acting as foundations, raises the question of the effective dynamic properties of these structured soils. Buildings, in the case of perfect elastic conditions for both soil and buildings, are shown to interact and strongly influence elastic body waves; such sitecity seismic interactions were pointed out in [Guéguen et al., Bull. Seismol. Soc. Am. 92, 794–811 (2002)], and we investigate a variety of scenarios to illustrate the variety of behaviours possible.
Key words: Seismic metamaterial / sitecity interaction / elastic cloaking / homogenisation / earthquake engineering
© B. Ungureanu et al., published by EDP Sciences, 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction: sitecity interaction
We assume that smart cities of the future will increasingly include more buried structures, especially in regions with large variations of temperature with the motivation being, for instance, to improve energy savings. This important substitution of soil by rigid or void elements will inevitably raise the question of the real dynamic response of the sedimentary soils at the free surface. In this article we explore what the theoretical tools borrowed from transformation optics can bring to the design of structured soils such as for a large scale photonic elastic crystal experimentally tested in [1], and we envision what a metacity cloak might achieve. Allied to the discussion of soil modification is the dynamic influence of buildings upon the elastic body waves and upon each other.
Seventeen years ago, a group of research scientists [4] made the bold claim that a megalopolis such as Mexicocity could have tremendous interactions with seismic waves, and deeply impact the registered ambient seismic noise. At that time, this research grouping might not have been fully aware of the deep consequences their work would have a decade afterwards, following the rise of metamaterials [6] for wave control. We show in Figure 1A the water front of Singapore, with its skyline, and for a physicist well versed in locally resonant structures such as metamaterials, highstorey buildings that can conceivably couple to waves with a wavelength much larger than their crosssectional size. Actually, one can view such a cluster of buildings as a scaledup forest of ‘concrete trees’. Some plate experiments with forests of metallic rods performed a few years ago by the group of Roux at University of Grenoble, have indeed shown significant filtering effects for Lamb waves [7], and these experiments motivated experiments on filtering effects of Rayleigh waves in a forest of trees [8].
Then, with the viewpoint that buildings are potential resonators it becomes clear that the Rayleigh waves could hybridize with the buildings [9] and behave in a way conceptually similar to spoof plasmon polaritons in optical metasurfaces [3]. A collection of buildings such as these at the waterfront of Singapore can be idealised using an elastic rainbow [2], as one can note that there is a gradient in the height of buildings in Figure 1A, and COMSOL simulations then predict socalled rainbow trapping in Figure 1B, where some of the buildings (the third, fifth and ninth, reading from smallest to tallest) can be seen to resonate with the seismic signal and other ones do not. This is reminiscent of the experimental and numerical analysis of the impact of three closely located towers in the city of Grenoble on local seismicity, and its perception by inhabitants, during an earthquake of small magnitude [10]. This study supports the hypothesis that buildings in urban areas need not, or indeed should not, be considered as standalone constructions with perfect embedding; they are usually studied as being independent dynamically. Our simulation in Figure 1B of an elastic rainbow effect has been investigated previously for forests of trees of an increasing height [2], and in the present case it could prove uncomfortable or even dangerous for inhabitants of the resonant buildings.
It is even possible that such a large scale metawedge [2] could convert some surface seismic wave activity into downward propagating waves, and thus shield other buildings within the megalopolis. The concept of the elastic rainbow is borrowed from the optics community [3], and this suggests intimate links between waves in optics and geophysics. For instance, during the summer of 2012, the civil engineering team of Stéphane Brûlé at the Menard company (Vinci Group) performed the first large scale experiments on a seismic shield [1] near the city of Grenoble, making use of some analogies with photonic crystals within which light is disallowed to propagate over certain frequency ranges, and the same team subsequently tested a flat seismic lens via negative refraction [11] near the city of Lyon. These two experiments confirmed that concepts from the physics of photonic crystals can be translated to the control of Rayleigh waves in structured sedimentary soils. Actually, analogies between optics and geophysics have been further extended to encompass locally resonant structures such as metamaterials going beyond the control of electromagnetic waves thanks to low frequency stop bands [6]. For instance, inertial resonators buried in the soil have been proposed to achieve shielding effects for low frequency seismic waves [12].
Fig. 1
Sitecity interaction and the rainbow effect: photo (A) of the waterfront in Singapore (courtesy of B. Ungureanu) and COMSOL simulation (B) of a compressional seismic wave, generated by a pulsating point force (of typical frequency 50 Hz) oriented as shown by black arrow, propagating in a concrete plate 40 m in thickness. This body Pwave interacts with a cluster of concrete buildings 20 m × 20 m in crosssection of increasing heights (from 10 to 100 m, by increments of 10 m) with centre to centre spacing of 40 m. The cluster is reminiscent of an elastic metawedge [2] inspired by trapped rainbow storage of light in graded metamaterials [3]. In light of the seminal work by Guéguen et al. [4], we know that the high storey buildings in the background in (A) can be viewed as a metacity which could display metamateriallike locally resonant features when a seismic wave propagates through it [5]. An elastic rainbow [2] effect could occur, which would filter the seismic wave frequencies in resonance with some buildings, as can be seen in the third, fifth and ninth buildings of the cluster in (B), reading from the smallest to the tallest building. 
2 Cloaking in optics and geophysics
With all the aforementioned works, it will not be a surprise that the theory of transformation optics introduced by Pendry et al. [13] and conformal optics introduced by Leonhardt [14], have had such a tremendous impact on the elastic wave community in the past ten years. Following Milton et al. [15], transformed elastodynamic equations have been studied for body elastic waves [16–18] and surface elastic waves in thin plates [19–22].
We show in Figure 2 the power of some wave analogies drawn between optics, hydrodynamics [23,24] and geophysics, that led to the concept of an invisibility cloak (Fig. 2C), working for microwaves (Fig. 2A), and that can serve as a basis for a design of a seismic cloak (Fig. 2B). The ideal goal of a physicist might be to render buildings invisible to seismic waves, but it is not a simple task, as civil engineers know. The design of a mechanical cloak as experimentally achieved by the group of Wegener [25] remains a challenge in the dynamic regime, so controlling inplane elastic wave trajectories is still an open problem. Of course, some great progress has been made on the control of outofplane, flexural, and Rayleighlike, wave trajectories in thin and thick plates [19–32], but there are other sorts of deleterious, surface Love, and body, seismic waves, that cannot be satisfactorily controlled yet, at least from a theoretical standpoint [33]. One might argue that at least controlling Rayleigh waves [34] would greatly reduce the damaging effect of an earthquake in sedimentary soils, as any building located inside the seismic zone would be in theory untouched by Rayleigh waves. There are other issues that need thought, for instance, there is a deep connection between cloaking and the theory of inverse problems, in the context of partial differential equations [35], that suggests a seismic cloak might just have the opposite effect as what one wishes to achieve. Indeed, from the standpoint of inverse problems [35], a perfect cloak should make the reconstruction of material properties based on boundary measurements impossible. Intuitively, an aspect of cloaking is that waves cannot penetrate the invisibility region within the cloak, so one expects that this can have interesting applications in wave protection. However, some fundamental theorems in inverse problems [36], show that while one could not reconstruct material properties within the cloak and its core, some trapped modes might occur (associated with properties of the socalled Dirichlet to Neumann map). Such trapped modes in the context of a seismic cloak might be even more devastating than the incoming seismic wave would have been in the absence of cloak! We will show a range of these features in the next section.
Fig. 2
Cloaking and wave analogies: water waves, [23] and also microwaves and sound waves [24], experiments have all been performed with the same geometrical cloak. The metamaterial cloak works over a finite interval of wave wavelengths that do not exceed the cloak's diameter (20 cm) and for which the crosssection of its structural elements can be considered small enough (not exceeding one third of the wavelength); (A) microwave experiment (courtesy of R. Abdeddaim and E. Georget, Institut Fresnel) showing both reduced field in the centre of the cloak and cloak's transparency for a transverse electric wave at 3.5 GHz [24]; (B) artistic view of a seismic cloak (Institute Fresnel/INSIS/CNRS); (C) photo of the fabricated aluminium cloak (courtesy of S. Enoch). 
3 Transformation seismology for the design of metacity cloaks
We now investigate the possibility of using interacting buildings and waves to try and create a design of a metacity that could have the ability to selfprotect itself or, at least, the centre of the city against earthquakes. We proceed as follows: first we consider a geometry like in Figure 2C with concentric layers of pillars. From our knowledge of transformation optics, we know that such a metamaterial corresponds to transformed anisotropic parameters that create an exclusion zone inside which an object is made invisible (known as a hole in electromagnetic space). Such a cloak, designed for electromagnetic waves, tells us that a scaledup cloak of increasing anisotropy towards its centre might well behave in a similar manner for elastic waves, see e.g. for a seismic cloak designed with Maxwellfisheyes or Luneburg lenses [37]: an exclusion zone could be created and a building placed inside this zone would be protected from elastic waves. Then we approximate the elastic parameters of these layers using an engineering algorithm that amounts to placing graded buildings of increasing height when one moves from the outer concentric ring of the cloak towards the inner ring, see Figure 3. Different designs have been tested and our aim here is to show the range of possible behaviours.
Such a plate design should work also for Rayleigh waves (and possibly Love waves), but in the next section we show numerical results for inplane coupled shear and pressure elastic waves, which according to transformation elastodynamic theory [16,18] require some twisting modes (elastic chirality associated with an effective elasticity tensor without the minor symmetry).
We have performed COMSOL computations for two types of cloaks: the first type consists of concrete buildings on a concrete (E = 25e9 [Pa], ρ = 2300 [kg/m^{3}], ν = 0.2) plate, Figure 4 and (left panels in Figs. 5–7) and the second type consists of steel (E = 200e9 [Pa], ρ = 7870 [kg/m^{3}], ν = 0.29) buildings in a concrete plate (right panels in Figs. 5–7). In Figure 4, we start with the most challenging case of a seismic source placed in the intense nearfield region of a potential cloak. Clearly, the wavefront behind the cloak is disturbed in panel (A), by comparison with benchmark for a body Pwave in homogeneous plate in panel (B), but upon inspection of panels (C) and (D), the field amplitude inside the centre of the cloak is clearly reduced. Such an effect is similar to the shielding experiment in [1], although in the present case the shield is an annulus rather than a slab. We shall see later that when the source is placed in the inside of the cloak, the wave pattern outside the cloak is nearly isotropic which is a hallmark of invisibility (see Fig. 8). So, the situation is more complex than it seems.
Let us now consider a seismic source far away from the cloak, it is now a plane body wave incident upon the cloak as in Figures 5–7. One can see that while almost perfect invisibility is achieved for the first type of cloak (concrete resonators) at 10 Hz (the seismic field is nearly unperturbed outside the cloak, Fig. 5A), the seismic wave field inside its core has a wave pattern reminiscent of what one would expect for a concentrator (the distance between wavefronts is decreased inside the core, see Fig. 5B). However, the second type of cloak (steel resonators) demonstrates good protection features in Figure 5B, compare also Figure 5C and D, but is as neutral as the cloak in (A) for its surroundings. When we increase the frequency, one can see in Figure 6 that cloaking is preserved to certain extent at 20 Hz, and the concentrating effect in (A), (C) disappears, whereas the protection in (B) is maintained. However, at 40 Hz, one can see in Figure 7 that both cloaks are no longer neutral for their environment. Actually, the second type of cloak in Figure 7B acts as a seismic wave shield and therefore buildings in front of it would be even more damaged, by the reflection, than in the absence of the cloak, but buildings placed after the cloak are in a shadow, relatively safe, zone, see Figure 7C which provides another perspective on the wave phenomenon.
These findings are a follow up of the proposal in [5] for shielding a place against seismic waves, by appropriately designed buildings of a city, used as aboveground resonators.
To show that the first type of cloak (concrete resonators) has indeed good invisibility features, we make a numerical experiment of a source placed inside its core at 10 Hz, see Figure 8A, C that shows that the seismic wavefield outside the cloak is almost isotropic, which is a hallmark of cloaking [21]. Obviously, when the frequency gets higher, see Figure 8B for a source at 50 Hz, homogenization breaks down and the seismic wavefield starts sensing the structural elements within the cloak, which now behaves as many secondary sources and a complex wave pattern outside the cloak with raylike wave trajectories in radial directions, ensues. One finally note from inspection of Figure 8C, D that the source placed inside the invisibility region leads to an enhanced seismic wave amplitude inside this cavitylike region, and to nearly suppressed seismic wave emission outside the cloak.
Fig. 3
Design of a large scale elastodynamic cloak in a plate with 10 concentric layers of pillars simulating buildings with height varying linearly from 10 m to 100 m, of square cross section 20 m × 20 m, with centre to centre spacing of 40 m along the radial direction. The plate has thickness of 40 m and the cloak has diameter of 1000 m, with a protected area in the centre which is 200 m in diameter and buildings have parameters of concrete (E = 25e9 [Pa], ρ = 2300 [kg/m^{3}], ν = 0.2) with 20 m of interstitial distance in between them. 
Fig. 4
Scattering of the metacity cloak (consisting of concrete buildings) by a point source at frequency 10 Hz (A), (C) and benchmark in concrete plate (B), (D); the seismic source is situated at a distance of 700 m from the centre of the cloak along a radial direction marked by red segment in (A), where the arc length is the Euclidean distance; Insets in (C) and (D) are the 3D views of (A) and (B), respectively. Blue curve in (C) and (D) shows profile of magnitude of total displacement field computed along the red segment in (A) and (B), respectively. Dark grey region in (C) marks the protected region inside the cloak with around 40 percent of reduction of the magnitude of total displacement in comparison with same region in (D). Light grey region in (C) marks the cloak itself. 
Fig. 5
Scattering of the metacity cloak by a plane wave source at 10 Hz; (A) concrete buildings on a concrete plate (E = 25e9 [Pa], ρ = 2300 [kg/m^{3}], ν = 0.2); (B) steel buildings (E = 200e9 [Pa], ρ = 7870 [kg/m^{3}], ν = 0.29) on a concrete plate. Color scale is as in Figure 4. One notes the antagonistic behavior of the field inside the core of the cloaks in (A), (C) and (B), (D). 
Fig. 6
Scattering of the metacity cloak by a plane wave source at 20 Hz; (A), (C) concrete buildings on a concrete plate (E = 25e9 [Pa], ρ = 2300 [kg/m^{3}], ν = 0.2); (B), (D) steel buildings (E = 200e9 [Pa], ρ = 7870 [kg/m^{3}], ν = 0.29) on a concrete plate. Color scale is as in Figure 5. 
Fig. 7
Scattering of the metacity cloak by a plane wave source at 40 Hz with different materials in (A) and (B): (A) concrete buildings on a concrete plate (E = 25e9 [Pa], ρ = 2300 [kg/m^{3}], ν = 0.2); (B) consists of steel (E = 200e9 [Pa], ρ = 7870 [kg/m^{3}], ν = 0.29) buildings on a concrete plate (right panel). Lower panel (C) is a 3D view of (B) Profile of magnitude of total displacement field integrated along the red segment in (A) and (B) suggests reduction of total displacement by up to 10% in the core of seismic cloak in (A) and reduction of total displacement by up to 80% in the core of seismic cloak as well as behind it in (B). 
Fig. 8
Cloaking of a seismic source located inside the core of the metacity cloak at 10 Hz (A), (C). These are concrete buildings on a concrete plate. One notices in (A), (C) that the field outside the cloak has reduced amplitude and is almost isotropic, which is a hallmark of cloaking a sensor in optics. The benchmark of a 10 Hz seismic source in a homogeneous plate in (D) makes clear the cavity like effect of the invisibility region in (A), (C). 
4 Open questions, new ideas, and future trends
We have presented some intriguing results that show buildings do interact with each other and the wavefield and can, under some circumstances, cloak a large scale elastic metamaterial. On the one hand, a perfect cloak would not prevent boundary measurements to be made for complete reconstruction of the elastic properties of a stiff object. In fact, depending upon the type of elastic cloaking theory − preserving the symmetries of the elasticity tensor at the cost of dealing with a transformed equation of the Willis type [15], or preserving the structure of the elasticity equation at the cost of breaking the symmetries of the elasticity tensor [20] – leads to different theoretical bounds on elastic cloaking [5]. This could have important consequences in the effectiveness of a seismic protection as the design of seismic cloaks relies upon some structured media that should approximate to certain extent the transformed parameters [11,25–27]. We believe that our two cloaks' (either concrete or steel resonators) can be further improved by better understanding how one can homogenize elastic plates with resonant buildings [38,39]. This is a delicate task which could be achieved by some reverse engineering algorithm using the transformed parameters in Willis's equation (which has fully symmetric rank 4, 3 and 2) tensors and doing same for the Cosserat case.
Our numerical simulations have shown that Young's modulus, Poisson's ratio and density play a prominent role in the way the seismic wave couples with the elements of the cloak. When the plate and its structural elements have same elastic parameters, we achieve good cloaking quality at low frequencies (10 Hz and 20 Hz), but protection is rather poor (at 10 Hz, we just achieve the opposite effect, of a field concentration inside the centre of the cloak). When the plate and its structural elements have markedly different elastic parameters, we achieve poor cloaking quality but good protection. Thus, changing the Young modulus leads to a shielding effect instead of cloaking as it can be seen in Figure 7C, which is very common in reality, buildings having very different properties compared with the alluvial soils, propagation medium of the seismic waves.
In fact, civil engineers compute the eigenmodes of each structure they wish to erect in order to avoid the interaction of the seismic wave with the structure they want to protect. This is achieved by choosing the characteristics of these structures in a way that they will absorb the minimum of the energy transported by the elastic wave so as to avoid the resonance.
Our proposal is that we should simply do the opposite to design structures in a smart way, in such a way that they will absorb exactly what we need to change the properties of these surface Rayleigh waves and so to obtain zones with zero seismic impact − the cloaked zones.
Another possibility is to convert and redirect the surface Rayleigh waves thanks to the transformation elasticity theory, or by filtering the surface waves frequencies thanks to the elastic rainbow effect [3].
Literature in geophysics and earthquake engineering is vast but ideas supporting the importance of sitecity interaction [4,10,40–42] are starting to be taken into consideration. It will take more time to convince the earthquake engineering community that ideas borrowed from the metamaterials community, notably regarding invisibility and cloaking [43], can help in the design of new cities − metacities.
In conclusion, the advantages of cloaking versus that of shielding protection is a nontrivial topic, which has hardly been addressed in the existing literature. The field of seismic metamaterials is still in its infancy and further theoretical and experimental results are required to answer the level of seismic protection one could achieve with structured soils in civil engineering. One could, for instance, envision that a skyscraper like that in Figure 9 may act as a directive antenna when an earthquake hits the megalopolis. This skyscraper might not only vibrate tremendously, but it might also act as a secondary seismic source which would modify the dynamic response of peripheral buildings. The way that architects could design cities of the future, with higher and higher buildings, might actually have deep consequences in terms of local seismicity and this study, and other previous ones based on geophysical measurements in situ, shows that one might wish to design a structure considering its close neighbours.
Fig. 9
Photo of a skyscraper in Bangkok (A) which stands out in the background (courtesy of B. Ungureanu). This building, modelled by a concrete pillar of height 300 m and square cross section 60 m × 60 m, atop a stressfree plate of thickness 40 m and infinite extent in the horizontal plane (modelled with elastic PMLs) should interact with an incoming seismic wave at 10 Hz according to a COMSOL simulation (B), with two smaller buildings placed in front being strongly impacted by the back reflection from the taller one. (C) and (D) simulations show the difference of the propagation of body waves considering 3 identical buildings or 10 buildings with height varying linearly from 10 m to 100 m, of square cross section 20 m × 20 m, with centre to centre spacing of 40 m along the propagation path, which behave like a wave guide path. These simulations exemplify the richness of sitecity interactions. 
Acknowledgments
Richard V. Craster acknowledges support of CNRS through the Unité Mixte Internationale Abraham de Moivre. Bogdan Ungureanu acknowledges funding of European Union (MARIE SKŁODOWSKACURIE ACTIONS project Acronym/Full Title: METAQUAKENG − Metamaterials in Earthquake Engineering − MSCA IF − H2020).
References
 S. Brûlé, E.H. Javelaud, S. Enoch, S. Guenneau, Experiments on seismic metamaterials: molding surface waves, Phys. Rev. Lett. 112, 133901 (2014) [CrossRef] [Google Scholar]
 A. Colombi, D. Colquitt, P. Roux, S. Guenneau, R. Craster, A seismic metamaterial: the resonant metawedge, Sci. Rep. 6, 27717 (2016) [CrossRef] [Google Scholar]
 K.L. Tsakmakidis, A.D. Boardman, O. Hess, Trapped rainbow storage of light in metamaterials, Nature 450 , 397 (2007) [CrossRef] [PubMed] [Google Scholar]
 P. Guéguen, P.Y. Bard, F.J. ChavezGarcia, Sitecity seismic interaction in Mexico city like environments: an analytical study, Bull. Seismol. Soc. Am. 92, 794–811 (2002) [CrossRef] [Google Scholar]
 S. Brûlé, B. Ungureanu, Y. Achaoui, A. Diatta, R. Aznavourian, T. Antonakakis, R.V. Craster, S. Enoch, S. Guenneau, Metamateriallike transformed urbanism, Innov. Infrastruct. Solut. 2, 20 (2017) [CrossRef] [Google Scholar]
 M. Kadic, T. Bückmann, R. Schittny, M. Wegener, Metamaterials beyond electromagnetism, Rep. Prog Phys. 76, 126501 (2013) [CrossRef] [Google Scholar]
 P. Roux, M. Rupin, F. Lemoult, G. Lerosey, A Colombi, R. V. Craster, New Trends Toward LocallyResonant Metamaterials at the Mesoscopic Scale, in: R.V. Craster, S. Guenneau, W.A. Kuperman, E.G. Williams (Eds.), World Scientific Handbook of Metamaterials and Plasmonics, 16, Elastic, Acoustic and Seismic Metamaterials, Vol. 2, 2017, pp. 251–299 [CrossRef] [Google Scholar]
 P. Roux, D. Bindi, T. Boxberger, A. Colombi, F. Cotton, I. Douste‐Bacque, S. Garambois, P. Guéguen, G. Hillers, D. Hollis, T. Lecocq, I. Pondaven, Toward seismic metamaterials: the METAFORET project, Seismol. Res. Lett. 89, 582–593 (2018) [CrossRef] [Google Scholar]
 A. Colombi, P. Roux, S. Guenneau, P. Guéguen, R. Craster, Forests as a natural seismic metamaterial: Rayleigh wave bandgaps induced by local resonances, Sci. Rep. 6, 19238 (2016) [CrossRef] [Google Scholar]
 P. Guéguen, A. Colombi, Experimental and numerical evidence of the clustering effect of structures on their response during an earthquake: a case study of three identical towers in the city of Grenoble, France, Bull. Seismol. Soc. Am. 106, 2855–2864 (2016) [CrossRef] [Google Scholar]
 S. Brûlé, E.H. Javelaud, S. Enoch, S. Guenneau, Flat lens effect on seismic waves propagation in the subsoil, Sci. Rep. 7, 18066 (2017) [CrossRef] [Google Scholar]
 Y. Achaoui, B. Ungureanu, S. Enoch, S. Brûlé, S. Guenneau, Seismic waves damping with arrays of inertial resonators, Extreme Mech. Lett. 8, 30–37 (2016) [CrossRef] [Google Scholar]
 J.B. Pendry, D. Schurig, D.R. Smith, Controlling electromagnetic fields, Science 312, 1780 (2006) [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 U. Leonhardt, Optical conformal mapping, Science 312, 1777 (2006) [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 G.W. Milton, M. Briane, J.R. Willis, On cloaking for elasticity and physical equations with a transformation invariant form, New J. Phys. 8, 248 (2006) [CrossRef] [Google Scholar]
 M. Brun, S. Guenneau, A.B. Movchan, Achieving control of inplane elastic waves, Appl. Phys. Lett., 94, 061903 (2009) [CrossRef] [Google Scholar]
 A. Norris, A.L. Shuvalov, Elastic cloaking theory, Wave Motion, 48, 525 (2011) [Google Scholar]
 A. Diatta, S. Guenneau, Controlling solid elastic waves with spherical cloaks, Appl. Phys. Lett. 05, 021901 (2014) [CrossRef] [Google Scholar]
 M. Farhat, S. Guenneau, S. Enoch, Ultrabroadband elastic cloaking in thin plates, Phys. Rev. Lett. 103, 024301 (2009) [CrossRef] [Google Scholar]
 N. Stenger, M. Wilhelm, M. Wegener, Experiments on elastic cloaking in thin plates, Phys. Rev. Lett. 108, 014301 (2012) [CrossRef] [Google Scholar]
 A. Darabi, A. Zareei, M.R. Alam, M.J. Leamy, Experimental demonstration of an ultrabroadband nonlinear cloak for flexural waves, Phys. Rev. Lett. 121, 174301 (2018) [CrossRef] [Google Scholar]
 D. Misseroni, D.J. Colquitt, A.B. Movchan, N.V. Movchan, I.S. Jones, Cymatics for the cloaking of flexural vibrations in a structured plate, Sci. Rep. 6, 23929 (2016) [CrossRef] [Google Scholar]
 M. Farhat, S. Enoch, S. Guenneau, A.B. Movchan, Broadband cylindrical acoustic cloak for linear surface waves in a fluid, Phys. Rev. Lett. 101, 134501 (2008) [CrossRef] [PubMed] [Google Scholar]
 J. Xu, X. Jiang, N. Fang, E. Georget, R. Abdeddaim, J.M. Geffrin, M. Farhat, P. Sabouroux, S. Enoch, S. Guenneau, Molding acoustic, electromagnetic and water waves with a single cloak, Sci. Rep. 5, 10678 (2015) [CrossRef] [Google Scholar]
 T. Bückman, M. Kadic, R. Schittny, M. Wegener, Mechanical cloak design by lattice transformation, Proc. Nat. Acad. Sci. 112, 4930–4934 (2015) [CrossRef] [Google Scholar]
 A. Colombi, P. Roux, S. Guenneau, M. Rupin, Directional cloaking of flexural waves in a plate with a locally resonant metamaterial, J. Acoust. Soc. Am. 137, 1783–1789 (2015) [CrossRef] [Google Scholar]
 A. Diatta, Y. Achaoui, S. Brûlé, S. Enoch, S. Guenneau, Control of Rayleighlike waves in thick plate Willis metamaterials, AIP Adv. 6, 121707 (2016) [CrossRef] [Google Scholar]
 A. Darabi, A. Zareei, M.R. Alam, M.J. Leamy, Experimental demonstration of an ultrabroadband nonlinear cloak for flexural waves, Phys. Rev. Lett. 121, 174301 (2018) [CrossRef] [Google Scholar]
 S. Brûlé, S. Enoch, S. Guenneau, R.V. Craster, Seismic metamaterials: controlling surface Rayleigh waves using analogies with electromagnetic metamaterials, Handbook of metamaterials (World Scientific, Singapore, 2018), Vol. 2, Chap. 7 [Google Scholar]
 A. Climente, D. Torrent, J. SánchezDehesa, Analysis of flexural wave cloaks, AIP Adv. 6, 121704, (2016) [CrossRef] [Google Scholar]
 T. Abbas, H. Ammari, G. Hu, A. Wahab, J.C. Ye, Twodimensional elastic scattering coefficients and enhancement of nearly elastic cloaking, J. Elast. 128, 203–243 (2017) [CrossRef] [Google Scholar]
 L. Pomot, S. Bourgeois, C. Payan, M. Remillieux, S. Guenneau, On form invariance of the KirchhoffLove plate equation, arXiv:1901.00067 (2019) [Google Scholar]
 R.V. Craster, A. Diatta, S. Guenneau, H. Hutridurga, Some results on nearcloaking for linear elasticity, arXiv:1803.01360 (2018) [Google Scholar]
 L. Rayleigh, On waves propagated along the plane surface of an elastic solid, Proc. Lond. Math. Soc. 17, 4–11 (1885) [CrossRef] [Google Scholar]
 R. Kohn, H. Shen, M. Vogelius, M. Weinstein, Cloaking via change of variables in electric impedance tomography, Inverse Probl. 24, 015016 (2008) [CrossRef] [Google Scholar]
 A. Greenleaf, Y. Kurylev, M. Lassas, G. Uhlmann, Cloaking Devices, Electromagnetic Wormholes, and transformation optics, SIAM Rev. 51, 3 (2009) [CrossRef] [Google Scholar]
 A. Colombi, S. Guenneau, P. Roux, R.V. Craster, Transformation seismology: composite soil lenses for steering surface elastic Rayleigh waves, Sci. Rep. 6, 25320 (2016) [CrossRef] [Google Scholar]
 C. Chesnais, C. Boutin, S. Hans, Wave propagation and nonlocal effects in periodic frame materials: Generalized continuum mechanics, Math. Mech. Sol. 20, 929–958 (2015) [CrossRef] [Google Scholar]
 R.V. Craster, J. Kaplunov, J. Postnova, Highfrequency asymptotics, homogenisation and localisation for lattices, Quart. J. Mech. Appl. Math. 63, 497–519 (2010) [CrossRef] [Google Scholar]
 A. Kham, J.F. Semblat, P.Y. Bard, P. Dangla, Seismic sitecity interaction: main governing phenomena through simplified numerical models, Bull. Seismol. Soc. Am. 96, 1934–1951 (2006) [CrossRef] [Google Scholar]
 J.F. Semblat, A. Pecker, Waves and vibrations in soils: earthquakes, traffic, shocks, construction works (IUSS Press, Austria, 2009) [Google Scholar]
 S. Brûlé, F. Cuira, Pratique de l'interaction solstructure sous seisme − Applications aux fondations et aux soutenements, AFNOR ISBN: 9782124656004, 2018 [Google Scholar]
 R. Fleury, F. Monticone, A. Alu, Invisibility and cloaking, origins, present, and future perspectives, Phys. Rev. Appl. 4, 037001 (2015) [CrossRef] [Google Scholar]
Cite this article as: Bogdan Ungureanu, Sebastien Guenneau, Younes Achaoui, Andre Diatta, Mohamed Farhat, Harsha Hutridurga, Richard V. Craster, Stefan Enoch, Stephane Brûlé, The influence of building interactions on seismic and elastic body waves, EPJ Appl. Metamat. 6, 18 (2019)
All Figures
Fig. 1
Sitecity interaction and the rainbow effect: photo (A) of the waterfront in Singapore (courtesy of B. Ungureanu) and COMSOL simulation (B) of a compressional seismic wave, generated by a pulsating point force (of typical frequency 50 Hz) oriented as shown by black arrow, propagating in a concrete plate 40 m in thickness. This body Pwave interacts with a cluster of concrete buildings 20 m × 20 m in crosssection of increasing heights (from 10 to 100 m, by increments of 10 m) with centre to centre spacing of 40 m. The cluster is reminiscent of an elastic metawedge [2] inspired by trapped rainbow storage of light in graded metamaterials [3]. In light of the seminal work by Guéguen et al. [4], we know that the high storey buildings in the background in (A) can be viewed as a metacity which could display metamateriallike locally resonant features when a seismic wave propagates through it [5]. An elastic rainbow [2] effect could occur, which would filter the seismic wave frequencies in resonance with some buildings, as can be seen in the third, fifth and ninth buildings of the cluster in (B), reading from the smallest to the tallest building. 

In the text 
Fig. 2
Cloaking and wave analogies: water waves, [23] and also microwaves and sound waves [24], experiments have all been performed with the same geometrical cloak. The metamaterial cloak works over a finite interval of wave wavelengths that do not exceed the cloak's diameter (20 cm) and for which the crosssection of its structural elements can be considered small enough (not exceeding one third of the wavelength); (A) microwave experiment (courtesy of R. Abdeddaim and E. Georget, Institut Fresnel) showing both reduced field in the centre of the cloak and cloak's transparency for a transverse electric wave at 3.5 GHz [24]; (B) artistic view of a seismic cloak (Institute Fresnel/INSIS/CNRS); (C) photo of the fabricated aluminium cloak (courtesy of S. Enoch). 

In the text 
Fig. 3
Design of a large scale elastodynamic cloak in a plate with 10 concentric layers of pillars simulating buildings with height varying linearly from 10 m to 100 m, of square cross section 20 m × 20 m, with centre to centre spacing of 40 m along the radial direction. The plate has thickness of 40 m and the cloak has diameter of 1000 m, with a protected area in the centre which is 200 m in diameter and buildings have parameters of concrete (E = 25e9 [Pa], ρ = 2300 [kg/m^{3}], ν = 0.2) with 20 m of interstitial distance in between them. 

In the text 
Fig. 4
Scattering of the metacity cloak (consisting of concrete buildings) by a point source at frequency 10 Hz (A), (C) and benchmark in concrete plate (B), (D); the seismic source is situated at a distance of 700 m from the centre of the cloak along a radial direction marked by red segment in (A), where the arc length is the Euclidean distance; Insets in (C) and (D) are the 3D views of (A) and (B), respectively. Blue curve in (C) and (D) shows profile of magnitude of total displacement field computed along the red segment in (A) and (B), respectively. Dark grey region in (C) marks the protected region inside the cloak with around 40 percent of reduction of the magnitude of total displacement in comparison with same region in (D). Light grey region in (C) marks the cloak itself. 

In the text 
Fig. 5
Scattering of the metacity cloak by a plane wave source at 10 Hz; (A) concrete buildings on a concrete plate (E = 25e9 [Pa], ρ = 2300 [kg/m^{3}], ν = 0.2); (B) steel buildings (E = 200e9 [Pa], ρ = 7870 [kg/m^{3}], ν = 0.29) on a concrete plate. Color scale is as in Figure 4. One notes the antagonistic behavior of the field inside the core of the cloaks in (A), (C) and (B), (D). 

In the text 
Fig. 6
Scattering of the metacity cloak by a plane wave source at 20 Hz; (A), (C) concrete buildings on a concrete plate (E = 25e9 [Pa], ρ = 2300 [kg/m^{3}], ν = 0.2); (B), (D) steel buildings (E = 200e9 [Pa], ρ = 7870 [kg/m^{3}], ν = 0.29) on a concrete plate. Color scale is as in Figure 5. 

In the text 
Fig. 7
Scattering of the metacity cloak by a plane wave source at 40 Hz with different materials in (A) and (B): (A) concrete buildings on a concrete plate (E = 25e9 [Pa], ρ = 2300 [kg/m^{3}], ν = 0.2); (B) consists of steel (E = 200e9 [Pa], ρ = 7870 [kg/m^{3}], ν = 0.29) buildings on a concrete plate (right panel). Lower panel (C) is a 3D view of (B) Profile of magnitude of total displacement field integrated along the red segment in (A) and (B) suggests reduction of total displacement by up to 10% in the core of seismic cloak in (A) and reduction of total displacement by up to 80% in the core of seismic cloak as well as behind it in (B). 

In the text 
Fig. 8
Cloaking of a seismic source located inside the core of the metacity cloak at 10 Hz (A), (C). These are concrete buildings on a concrete plate. One notices in (A), (C) that the field outside the cloak has reduced amplitude and is almost isotropic, which is a hallmark of cloaking a sensor in optics. The benchmark of a 10 Hz seismic source in a homogeneous plate in (D) makes clear the cavity like effect of the invisibility region in (A), (C). 

In the text 
Fig. 9
Photo of a skyscraper in Bangkok (A) which stands out in the background (courtesy of B. Ungureanu). This building, modelled by a concrete pillar of height 300 m and square cross section 60 m × 60 m, atop a stressfree plate of thickness 40 m and infinite extent in the horizontal plane (modelled with elastic PMLs) should interact with an incoming seismic wave at 10 Hz according to a COMSOL simulation (B), with two smaller buildings placed in front being strongly impacted by the back reflection from the taller one. (C) and (D) simulations show the difference of the propagation of body waves considering 3 identical buildings or 10 buildings with height varying linearly from 10 m to 100 m, of square cross section 20 m × 20 m, with centre to centre spacing of 40 m along the propagation path, which behave like a wave guide path. These simulations exemplify the richness of sitecity interactions. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.