arXiv:1210.5589v2 [cond-mat.mtrl-sci] 1 Mar 2013 aeegh[,2.Tewvlntso ufc waves surface 100 of of wavelengths order The the of 2]. are [1, wavelength a rvlsoe hnbd ae n h mltdsde- amplitudes travel waves waves 1 epi- the Surface Surface about the depth. and the waves with and generation. exponentially body crease earth of than point the slower main of travel body the surface when is the generated center is at wave arrive surface waves are A waves Love(L) waves. and Rayleigh(R) surface waves and waves, Secondary(S) waves body and body are Primary(P) waves: seismic waves. of surface types and two are There wave. directly is that surface Earth’s the hypocenter. the the of epi- on above depth The point the the kilometers. and of is tens center fire’ exceed of not of ‘ring does mantle orig- hypocenters so-called the earthquakes the tectonic of of at parts majority inate upper The very [1]. crust the Earth the the at between or Most border mantle the begins. and at is rupture located energy earthquake are strain the hypocenters stored the and where released point the first is hypocen- focus The the or in waves. ter energy seismic creates of which amounts crust, huge Earth’s of release sudden the from eoe ehd rvd e a orss seismic resist to way de- new recently a However, provide thus earthquakes. methods are from veloped and safe waves not seismic the still control earthquakes. effective, cannot to response generally Humans passive are a comprise methods basically These but swing- vibration. ensure from or collapse to not ing do aimed themselves buildings are the engineering that geotechnical and tural property. most and the life to produce dura- hazards which long serious amplitudes, of causing of large frequencies destruction that low of than of and slowly are tions Those more decays waves. waves body surface of tude 30 esi aei ido nooeeu acoustic inhomogeneous of kind a is wave seismic A results disasters, largest nature’s of one earthquake, An ovninlerhukpofigmtoso struc- of methods earthquakeproofing Conventional Hz b eryblwteadbefeuny h ampli- The frequency. audible the below nearly , eateto hoeia hsc,RP,TeAsrla Na Australian The RSPE, Physics, Theoretical of Department ∼ a 3 iiino aieEgneig op ainlMrtm Un Maritime National Mokpo Engineering, Marine of Division km/sec ASnmes 10.j 16.j 91.30.-f 91.60.Lj, 81.05.Xj, numbers: PACS tested is zone barrier. shadow normal seismic a The with zone. creates gene the meta-barrier level o within The protection sound buildings point method. the the cloaking estimated from traditional differs We the mecha protection a area heat. of is and method This sound energy. redu into seismic which energy the waves, dissipating seismic i by becomes the wave nentially attenuates seismic meta-boxes the b of many velocity and of empty waves the protect, seismic huge to of want designing frequencies we resonance By the to metamaterials. sponding acoustic using zone edvlpdanwmto ferhukpofegneigto engineering earthquakeproof of method new a developed We rica esi hdwZn yAosi Metamaterials Acoustic by Zone Shadow Seismic Artificial n aygetywti h et of depth the within greatly vary and m n h rqece r below are frequencies the and agHo Kim Sang-Hoon Dtd ac ,2013) 4, March (Dated: a ∗ n uud .Das P. Mukunda and xrsino h form: the of expression sound. a into when converted level is sound of wave the simulation seismic of we calculation the computer and paper a zone all this shadow including the also In method but the build- metamaterials. metamaterials, report the the by only behind surrounded not buildings protects is method that This stop-band ing the wave. of velocity the wave imaginary of attenuated the of an use into attenuator making wave by an us- seismic zone as the shadow converting act by artificial Metamaterials an metamaterials. create sug- ing to We method [12–15]. a introduced gested previously are that method jeop- rather buildings. but surrounding problem, the the solve ardizes cloaked not the does behind This buildings how- region. they to waves, destructive as seismic still cloaked building are The ever, the plate. of the waves foundation through The pass the ring. from innermost in- away the refractive to bend low ring the outer the at from flexible dex most high the to a at pro- index is stiffest from refractive design arranged They Their rings concentric of plates. propagated plate waves thin [6–8]. heterogeneous bending al isotropic of in design et theoretical Farhat a posed by suggested method types two are seismic There attenuating and deflecting cloaking. waves. of methods: is, cloaking form metamaterials of a with design principle, wave- seismic a in The of that [3–5]. than di- length smaller with much potentially structures mensions are homogeneous Metamaterials effectively man-made metamaterials. acoustic using waves i eim[2.Sermodulus, Shear [12]. medium mic ueytodmninl eoiyi anydependent mainly density, is the velocity on dimensional, two purely property) h hs eoiyo l ehnclwvsflosthe follows waves mechanical all of velocity phase The cloaking extended an is method attenuated The cloaking traditional a is waves seismic Deflecting inlUiest,Cner,AT00,Australia 0200, ACT Canberra, University, tional esi hdwzn,poetn l the all protecting zone, shadow seismic a ovninlsimcdsg,including design, seismic conventional f ia ehdo ovrigteseismic the converting of method nical rigte rudtebidnsthat buildings the around them urying ycmue iuainadcompared and simulation computer by aiay h eabrircomposed meta-barrier The maginary. } e h mltd ftewv expo- wave the of amplitude the ces vriy op 3-2,R .Korea O. R. 530-729, Mokpo iversity, 1 oe ihafwsd-oe corre- side-holes few a with boxes / 2 raea rica esi shadow seismic artificial an create ae rmasimcwv.This wave. seismic a from rated lhuhsimcsraewvsaenot are waves surface seismic Although . ρ n ha modulus, shear and , b † v = { eatcproperty)/(inertial (elastic E sdfie by defined is , E fteseis- the of , 2

∆P = E∆x/h, where ∆x is the horizontal shift and h is the height of the object. The negative shear modulus causes an axisymmetric deformation under an opposite axisymmetric loading. The negative shear modulus of elastic media has been realized several times by Helmholtz resonators or its ap- plications [16–18]. The resonance of accumulated waves in the Helmholtz resonator reacts against the applied pressure at some specific frequency ranges. Then, the negative modulus is realized by passing the acoustic wave through an array of Helmholtz resonators called “meta- boxes.” Therefore, the acoustic intensity decays at some resonant frequency ranges. Considering the structural loss, the general form of the effective shear modulus of FIG. 1: Meta-boxes with 6 side holes and the positioning elastic materials, Eeff , is given similarly with the general form of the bulk modulus [9–12, 19, 20]. of the 4 meta-boxes. The empty spaces among the boxes for passage between the holes act as a large capacitor. 2 1 1 Fω E− = E− 1 o , (1) eff  − ω2 ω2 + iΓω  − o the imaginary wave-vector at the frequency ranges of the negative modulus. where ω is the resonance frequency, Γ is a loss by damp- o The size of the meta-boxes can be estimated by analogy ing, and F is a geometric factor [21]. The effective shear between electric circuits and mechanical pipes. A pipe or modulus has a negative value at 1 < ω/ω < √1+ F o tube with open ends corresponds to an inductor, and a when the loss is small. The negative range of the real closed end corresponds to a capacitor [12]. The resonance part is the stop-band of the wave. frequency of a Helmholtz resonator is [21–23] If the shear modulus becomes negative, then the ve- locity becomes imaginary, as does the refractive index, S ω = v, (3) n, or the inverse of the velocity as n = vo/v = vo ρ/E, o rl′V where vo is the background velocity. It has imaginaryp √ where S is the area of the cross-section, V is the volume, impedance of absorption: ρv = ρE. Because the wave- ′ vector, k = 2πn/λ, becomes imaginary, too, the ampli- and v is the background acoustic velocity. l is the effec- tive length which is given by l′ l +0.85d, where l is the tude of the seismic wave is attenuated. It functions as a ≃ “barrier” or “attenuator” to the seismic waves. length of the hole or thickness of the meta-box and d is The intensity of an earthquake or the earthquake mag- the diameter of the hole [23]. nitude, M, is defined by comparing the two amplitudes An example of the design of the meta-boxes for the in a logarithmic scale [1, 2]. seismic frequency range is plotted in Fig. 1. It is a 6- hole cylinder or rectangular box. The diameter of the A hole and the thickness of the cylinder is on the order of M = log , (2) 3 Ao 0.3m or 1 foot. The volume inside could be 10 100m . The greater the number of holes, the higher∼ the reso- where A is the maximum amplitude of the seismic wave nant frequency. The shape of the meta-box could be any and Ao is the maximum amplitude of the background form of a concrete box with a few side holes. Because vibration. Ao depends on the distance from the epicenter the meta-boxes are considerably smaller than the cor- and the order of micrometers. If A is recorded at 100km responding wavelengths, the array of the boxes behaves from the epicenter and Ao = 1µm, then the magnitude as a homogenized medium. Various types of resonators becomes Richter-scale or local magnitude ML. may cover various types of resonance frequencies of the Here we focus on reducing the amplitude of the seismic seismic waves. wave using the properties of the metamaterials. The en- A seismic wave is an elastic wave that can be approx- ergy of the seismic waves is dissipated inside of the meta- imated by a sine wave. Assuming a plane wave of wave- boxes, and the absorbed energy is converted into sound length λ propagates in x-direction, the amplitude of the and heat [22]. The seismic energy is converted into sound wave reduces exponentially [12]. the same way breathing pressure makes a resonant sound ikx −2π|n|x/λ in a musical instrument. We built an attenuator or an Ae = Ae . (4) earthquakeproof barrier to a seismic wave by filling up We can rewrite Eq. (4) with the definition of the magni- many resonators under the ground around the building tude in Eq. (2). that we want to protect. The amplitude of the seismic 2 wave that passes through the meta-boxes is reduced by 10Mi e− π|n|x/λ = 10Mf , (5) 3

FIG. 2: A design of ∆M = 3 meta-barrier made of many meta-boxes: w = 60m, L = 300m. The seismic shadow zone and its diffraction areas are shown.

where Mi and Mf are the seismic magnitude before and after entering the seismic barrier. From Eq. (5) we obtain the thickness of the seismic barrier ∆x [12].

ln 10 λ∆M 0.366λ ∆x = = ∆M, (6) 2π n n FIG. 3: Pressure distribution by (a)a meta-barrier in Fig. 2 | | | | and (b)a normal barrier. freq.= 10Hz. where ∆M = Mi Mf . Note that a high refractive index material is desirable− for the narrow seismic barrier. The depth of the seismic barrier should be at least the same as E = 1.0GP a for the background. The background the foundation of the building, but it is not necessary to could be rocks, water, or mixed. We applied ρ = be more than the wavelength of the surface waves. Mag- 2.0 103kg/m3 and E = 0.50GP a for the meta-barrier, nitude 3 or lower earthquakes are almost imperceptible and×E =0.50GP a with the− same density for the normal and rarely cause damage. On the other hand, magnitude barrier. Note that the normal barrier has the refractive 6 or higher ones are very destructive, and it is expensive index 2 and matches the impedance ρv with the back- to build buildings protected from these large earthquakes ground. We see from the Fig. 3(a) that the meta-barrier using conventional earthquakeproof design. The largest creates the seismic shadow zone and diffraction areas. On magnitude earthquakes on record were on the order of 9 the other hand, the normal barrier in Fig. 3(b) can not [1]. Most buildings in earthquake belts have a seismic de- stand from the seismic wave. sign of 3 M 6. Therefore, ∆M = 3 is the minimum Converting mechanical energy into acoustic energy cre- ≤ ≤ magnitude required to effectively protect the buildings. ates sound and heat. It is simple to estimate the sound We designed an artificial seismic shadow zone in Fig. level. The acoustic intensity or sound level in deci- 2 by constructing a meta-barrier which was composed of bels(dB) is defined by [22, 23] many meta-boxes. The other side of the meta-barrier of I the seismic wave is the seismic shadow zone protecting all β = 10 log , (7) of the buildings. The width of the barrier is on the order Iref of the seismic wavelength, but the length of the barrier should be much longer than the wavelength because there where Iref is the reference sound intensity at the thresh- old of hearing or I =1.0 10−12W/m2. The intensity are diffraction waves from the both ends of the barrier. ref × Then, we may need an order of hundreds of thousand is proportional to the square of its amplitude. Compar- meta-boxes. If the relative refractive index of the meta- ing Eq. (2) with (7), we can roughly estimate the sound box is n = 2 and the wavelength of the surface wave is level in decibels. | | λ = 100m, a width of ∆x 60m for ∆M = 3 resistance. A The shape of the barrier is≃ neither circular nor concentric β 20 log 20M, (8) ≃ Aref ≃ rings like in the traditional clocking method [6–8]. We have shown that water trench is not satisfactory for the where Aref is the reference amplitude and A is the root- attenuation of seismic wave [12]. mean-square of the amplitude of the seismic wave. Note The seismic shadow zone by a meta-barrier was dis- that M is the magnitude at the barrier and is different placed using COMSOL multiphysics and compared with from Richter’s scale ML. a normal barrier. We applied some typical material pa- The sound level in Eq. (8) can be rewritten with rameters. Density and modulus are ρ =1.0 103kg/m3, Richter’s scale M . The intensity of a three-dimensional × L 4 wave decreases with 1/r2, and that of a two-dimensional nuclear reactors, oil refineries , long-span bridges, express wave decreases with 1/r, where r is the distance from rail-roads, historical monuments, among others. the source. A seismic wave is a combination of these two. Then, the magnitude M is written as [1, 2] This research was supported by Basic Science Research Program through the National Research Foundation of M = a b log D, (9) Korea(NRF) funded by the Ministry of Education, Sci- − ence and Technology(2012-0001425). where D is the distance from the epicenter, and a, b are empirical positive constants. b is the power of the distance and order of 1.5. Because Richter’s scale ∗ Electronic address: [email protected] ML is defined at D = 100km, Eq. (9) is written as † M = a b log 100. From the relation between M and Electronic address: [email protected] L − [1] R. Villaverde, Fundamental Concepts of Earthquake En- ML, the sound level in Eq. (8) is expressed with ML gineering, Ch. 4, 5, App. (CRC, New York, 2009). instead of M. [2] A. Udias, Principles of Seismology, Ch. 15 (Cambridge, New York, 2010). 100 5 β 20 ML + b log , (10) [3] J. B. Pendry, Nature Phot. , 599 (2006). ≃  D  [4] D. R. Smith, J. B. Pendry, and M. C. K. Wiltshire, Sci- ence 305, 788 (2004). where D is measured in km. The dominant factor is [5] J. Li and C. T. Chan, Phys. Rev. E 70, 055602(R) (2004). not the distance D, but the magnitude ML. If ML = 6 [6] M. Farhat, S. Guenneau, S. Enoch and A. B. Movchan, and D = 100km, then β 120dB, which corresponds Phys. Rev. B 79, 033102 (2009). to a siren or rock concert.≃ However, the sound level is [7] M. Farhat, S. Guenneau, and S. Enoch, Phys. Rev. Lett. 103 actually lower than this because some part of the seismic , 024301 (2009). [8] M. Farhat, S. Guenneau, and S. Enoch, Phys. Rev. B 85, energy is converted into heat by friction. The frequency 020301(R) (2012). will be near infrasonic and barely audible. [9] N. Fang, X. Xi, J. Xu, M. Ambati, W. Srituravanich, C. We developed a different kind of cloaking method for Sun and X. Zhang, Nature Mater. 5, 452 (2006). the attenuation of seismic waves. The method produces [10] S. H. Lee, C. M. Park, Y. M. Seo, Z. G. Wang, and C. an artificial seismic shadow zone that abolishes the seis- K. Kim, J. Phys. Condens. Matter 21, 175704 (2009). mic wave by changing the seismic wave-vector to imagi- [11] S. H. Lee, C. M. Park, Y. M. Seo, Z. G. Wang, and C. 104 nary using meta-boxes. We designed a meta-barrier and K. Kim, Phys. Rev. Lett. , 054301 (2010). [12] S.-H. Kim and M. P. Das, Mod. Phys. Lett. B 16, tested it by a computer simulation. This method does not 1250105 (2012). just add another seismic system to a building, but rather [13] MIT, Technology Review, Feb. (2012), constructs an earthquakeproof barrier around the build- http://www.technologyreview.com/blog/arxiv/27566/ ing to be protected. This barrier is a kind of waveguide [14] A. Piore, Discover Magazine, 70-72, Jul./Aug. (2012), that exponentially reduces the amplitude of the danger- http://discovermagazine.com/2012/jul-aug/06-how-to-make-anything-disapp ous seismic waves. The seismic energy in the barrier is [15] S.-H. Kim, J. Acoust. Soc. Am. 131, 3292 (2012). 39 converted into a very low frequency sound. [16] Y. C. Wang, and R. S. Lakes, J. Comp. Mat. , 1645 This method differs from the traditional cloaking (2005). [17] Y. Wu, Y. Lai, and Z.-Q. Zhang, Phys. Rev. Lett. 107, method of deflecting seismic waves. The traditional 105506 (2011). method sacrifices the buildings behind the meta-regions. [18] X. Zhou, and G. Hu, Phys. Rev. B 79, 195109 (2009). In addition, the seismic barrier method has several ad- [19] Y. Cheng, J. Y. Xu, and X. J. Liu, Appl. Phys. Lett. 92, vantages. The building we want to protect does not have 051913 (2008). to be altered, and the seismic range of the building can [20] Y. Cheng, J. Y. Xu, and X. J. Liu, Phys. Rev. B 77, be upgraded or downgraded by adjusting the width of 045134 (2008). [21] C. Ding, L. Hao, and X. Zhao, J. Appl. Phys. 108, 074911 the barrier. This method will be effective for buildings (2010). that are not well-equipped with conventional seismic de- [22] L. E. Kinsler, A. R. Frey, A. B. Coppens, and J. V. signs. The biggest advantage is that all buildings behind Sanders, Fundamentals of Acoustics, 4th edn. Ch. 5. the barrier will be protected, which decreases the cost of (John Wiley & Sons, New York, 1999). seismic protection. This method is applicable for social [23] Beranek, L. L. Acoustics, Ch. 1, 3, 5. (McGraw-Hill, New overhead capitals, such as power plants, dams, airports, York, 1954).