We study the analogy between the linear elasticity theory equations and classical Maxwell equation with chiral effects and we propose a new method of an earthquake-resistant design to support conventional aseismic designs using acoustic metamaterials. We suggest a simple and practical method to reduce the amplitude of a seismic wave exponentially. Our device is like an attenuator of a chiral seismic wave. Constructing a cylindrical shell-type waveguide that creates a stop-band for the chiral seismic wave, we can convert the wave into an evanescent wave for some frequency range without touching the building we want to protect.

Earthquakes are the result of sudden release of huge amount of energy in the Earth’s crust that produces seismic waves. There is a lot of observational evidence concerning rotational motions on the Earth’s surface excited by earthquakes, we can infer the existence of rotational motions in earthquake focal zones. It is observed that rotational motions (twist) in earthquake sources naturally excite rotational seismic waves, solitons and chiral waves [1-3].

The key idea of this article consists in the description of rotational seismic waves into a chiral regime [4,5] by considering chiral rotational seismic waves excited in earthquake sources in strong earthquakes. We can distinguish two kinds of rotational seismic waves: 1) rotational longitudinal waves, that is, PR waves; and 2) rotational shear waves, that is, SR waves. Rotational seismic waves propagate faster in solid rocks and much slower in fractured media along tectonic faults. It has been observed that these waves may have a form of rotational seismic solitons and that they can trigger earthquakes [

It is the collapse of structures that causes extensive damage and loss of life during earthquakes. Aseismic capabilities are highly relevant to public safety and a large amount of research has gone into establishing practical analysis and design methods for them.

The main principles underpinning the protection of structures from seismic waves are applied to waves of force that travel through the Earth, which is an elastic body. Travelling pressure waves are always generated with longitudinal polarization. This means that the ground is alternately compressed and dilated in the direction of propagation. Travelling shear waves are transverse, or perpendicular to the direction of propagation. One of the effects of the metamaterial is to de-couple these components of the seismic wave. By applying these principles, it is possible to design a cylindrical, seismic metatmaterial cloak to remedy the destructive effects of coupled pressure and shear waves that would travel through the same earthen geometrical plane as the structure. The components of the seismic waves should work together to disrupt the elastic earth, on which a building is anchored.

Recent development of metamaterial science opens a new direction to control the seismic waves. Farhat et al. proposed a design of a cloak to control bending waves propagated in isotropic heterogeneous thin plates [7-9].

Their cloaking of shear elastic waves pass smoothly into the material rather than reflecting or scattering at the material’s surface. In this paper we propose a method to control seismic waves by using a new class of materials called metamaterials. Metamaterials are artificially engineered materials which have a special property of negative refractive index.

We present here a solution that a metamaterial acts as an attenuator by converting the destructive seismic wave into an evanescent wave by making use of the imaginary velocity of stop-band of the wave.

There are many representations to express the scale of earthquakes. Among them the magnitude that comes from the amplitude of the seismic waves is most important. The common form of the magnitude, R, is the Richter-scale defined by comparing the two amplitudes in a logarithmic scale as, where is the maximum amplitude of the seismic wave and is the maximum amplitude of the background vibration and order of µm. We focus here on how to reduce the amplitude of the seismic wave by using the properties of chiral metamaterials. In Section 2 of this paper, we develop a theory for chiral seismic waves. Based on the similarity of elastodynamics and electrodynamics we discuss the region of elastic metamaterial with negative modulus and negative density that enable a seismic wave becomes evanescent wave. In Section 3, we discuss a scheme of seismic attenuator. The novelty of this proposal is that the attenuation is due to the combined effect of metamaterial and evanescence in S waves.

In connection with seismic events, the Elastodynamics is linked with Electrodynamics through the transverse waves. Seismic waves encountering interfaces that separate rocks of different elastic properties also undergo reflection, refraction, and scattering phenomena and chiral materials are noncentrosymmetric due to handedness in their microstructures. The elastic field behavior in a chiral medium is readily described herein using the governing equations and constitutive relations for noncentrosymmetric, isotropic micropolar materials. Accordingly, linearly polarized longitudinal waves and left and right circularly polarized transverse waves are eigenstates for elastic waves in the chiral medium. There are two classes of seismic body waves which travel through the interior of the Earth: P waves and S waves. P waves are the fastest seismic waves, and consequently, the first to arrive at any given location. Torsion waves, often called S waves, represent the spiraling motion of particles twisting between inner structures. P waves are usually the first to be recorded on a seismogram because they travel the fastest. S waves usually have more height, or amplitude, than P waves.

We will be started with consideration of an analogy between the linear elasticity theory equations and classical Maxwell equation with chiral effects. Chiral media have the property to split two degenerate modes [10-14]. For a chiral medium with an electric or magnetic resonance, a negative refractive band is created and there is a continuous transition between the negative and positive bands, as described in [

where the vector represents the vector of displacements is a elastic medium, are the velocities of propagation of the longitudinal (P waves) and transverse (S waves) waves correspondingly. In order to separate the wave processes explicitly we introduce the chiral seismic parameter which take into account the mode conversion from P to S waves with the time varation as

Here we identify so using the Born-Fedorov formalism

and

the equivalent Maxwell equations are

Solving these equations we have that the effective phase velocity of S-waves is where is the wavevector without and the wave equation is

Metamaterials with negative electromagnetic constitutive parameters give new propagation characteristics for electromagnetic waves [10-14]. Negative permittivity occurs naturally due to the plasma oscillation, but negative permeability became available only artificially [

Although from an electromagnetic point of view, chiral homogeneous material can be described by different specific equations, by analogy between elastodynamics and electrodynamics in this work we use the BornFedorov equations as the most suitable for applications of our interest [5-6]. Here the density of matter is equivalent to the electric permittivity, and the Lamé parameter is equivalent to the reciprocal of the magnetic permeability. The speed of transversal seismic S-waves without chirality is given by. The chiral effect given by the factor which modify the phase velocity .

On the other hand, acoustic waves are created by compressibility or elasticity of the medium. Young’s modulus, Y, is a one-dimensional compressibility defined by, where is the pressure or stress and l is the length. Shear modulus, , is a two-dimensional one for a surface wave defined by, where is the horizontal shift and h is the height of the object. Bulk modulus, , is a three-dimensional one for a body wave defined by

Seismic medium of Earth crust can be considered as an accumulation of infinite number of elastic plates. Although the seismic surface wave is not pure two dimensional, the velocity is mainly dependent upon the density, and shear modulus, , of the seismic medium. Seismic wave is a kind of acoustic wave and every acoustic wave propagates following two wave equations in principle. Assuming the plane wave time dependence, the pressure, p, and the velocity, of the S-wave in the two dimensions are expressed as the Newton’s law [

and the continuity equation

where is the Laplacian operator at the surface, p is the pressure, ω is the angular frequency of the wave, and is the velocity.

The Equations (4) and (5) generates the wave equation as

if we compare Equation (3) with Equation (6), we complete the analogy and connection between the chiral electromagnetism and the acoustic seismic waves. The velocity of the seismic wave is

The Equation (4) is the counterpart of the Ampere’s law and Equation (5) is of the Faraday’s law by the analogy of electromagnetism and acoustic mechanics.

If the shear modulus or the density becomes negative, the velocity becomes imaginary. Then, so does the refractive index, n, or the inverse of the velocity as, where is the background velocity. Since the wavevector becomes imaginary, too. Therefore, the imaginary wavevector makes the amplitude of the seismic wave become an evanescent wave. We call it an attenuator. Negative shear modulus of elastic media has been studied and realized very recently [16-18] (Helmholtz resonators with the negative bulk modulus). 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.

In

Therefore, the sound intensity decays exponentially at some resonant frequency ranges. In general the elastic material is described by three independent effective parameters of, and ρ [

From the formalism of electromagnetic response in metamaterials, effective electric permittivity and effecttive magnetic permeability show negative values at some specific frequency ranges around resonances [19,20]. The Helmholtz resonator is a realization of an electrical resonance circuit by mechanical correspondence. It is known as that the ferriplasmon frequency in metals or in an array of ring metal wires produces the magnetic permeability as [

where is the ferriplasma frequency and is a

loss by damping.

Considering the structural loss, the general form of the effective shear modulus, , is given similarly with the general form of the bulk modulus as [10-14]

where is the resonance frequency and F is a geometric factor [15,21,22]. The real and imaginary part of the is plotted in

which will be proportional to density, so

which may be negative too and in this case the phase velocity is negative near of resonance conditions (and are simultaneously negatives) or when .

We can build an attenuator or a barrier of a seismic wave by filling-up many resonators under the ground around the building that we want to protect. Then, the amplitude of the seismic wave that passed the waveguide is reduced exponentially by the imaginary wavevector at the frequency ranges of negative modulus. Mixing up many different kinds of resonators will cover many different corresponding frequency ranges of the seismic waves. The chiral factor can enhance this attenuation if we take with the appropriate sign [23-26].

If we assume that the plain seismic wave of wavelength λ propagates in x-direction, the amplitude of the wave reduces exponentially as

Let the initial seismic wave, that is, before entering the waveguide, have amplitude and magnitude, and final seismic wave, that is after leaving the waveguide, have amplitude and magnitude. Then, is written as

The amplitude of the seismic wave reduces exponenttially as passing the waveguide of metamaterials. We can rewrite Equation (9) with the definition of the magnitude in

Taking logarithms both sides of Equation (10), we obtain the width of the waveguide, , as

where

For example, if the refractive index is and, the wavelength of the surface wave is , we need to reduce. If we need That is the chiral factor can substantially reduce the volume of materials, with the significant reduction in engineering costs.

If the aseismic level of the building is and the width of the waveguide surrounding the building is about m, then the effective aseismic level of the building is increased to Therefore, a high refracttive index material can be obtained by the factor for a narrow attenuator. In civil engineering earthquake proofing methods must be practical, that is, with great versatility, simplicity and easy to construct.

The resonator should be easy to build. We show an example of a resonator in

The size of the cylinder can be estimated from the analogy between electric circuits and mechanical pipes.

An example of chiral metamaterials, is chiral Swiss Rolls, which are constructed by wrapping a conducting sheet around a cylindrical mandrel in a helical shape (as shown in

chirality, at least 50 times larger than other chiral structures previously reported in the literature. This property ensures the change of a wave polarization by 90˚ in less than a wavelength. Hence, chiral Swiss Rolls support a negative refractive band for one-wave polarization, which is predicted analytically and calculated numerically and therefore a macroscopic resonant magnetic behaviour is created. Chirality arises from the helical shape of the structure.

An example of the design of the meta-cylinder for the seismic frequency range is followings: the diameter of the hole is order of 0.1 m, the thickness of the cylinder is order of 0.1 m, and the volume inside is order of. Since the meta-cylinder is considerably smaller than the corresponding wavelengths, the array meta-cylinders behave as a homogenized medium. The chiral factor can reduce the phase velocity so the S-wave can be stopped.

The shape of the meta-cylinder is neither necessary to be circular nor to have 4 holes. It could be any form of a concrete box with several side holes. Cubic or hexagonal boxes would be fine. Various kinds of resonators may cover various kinds of resonance frequencies of the seismic waves. There happens an energy dissipation of the seismic waves inside of the waveguide and the absorbed energy will turn into sound and heat. They make the temperature of the waveguide increasing depending on the magnitude of energy that arrives at the waveguide [20,21].

A vertical view of the metamaterial barrier with many meta-cylinders are in

Seismic waves cannot pass through water. Then, we can imagine a water barrier of a big trench filled with water. However, it is not stable to stand a series of attacks of “foreshock → main shock → aftershock.” Because, if the outer part of the water trench is brought down by a seismic wave, we do not have enough time to rebuild it to prepare for the main shock and the after shock. Water cannot have a high refractive index n in Equation (11), too. Maintaining the depth of the water trench up to the wavelength of the seismic wave require huge amount of water and, therefore, not be practical.

We introduced a supportive method for aseismic design. with an earthquake proof barrier around the building to be protected. This barrier is a kind of waveguide that reduces exponentially the amplitude of the dangerous seismic waves. Then, if the structures are protected with ground rubber layers would be possible to divert the S waves which do not impact a building.

The cloak can use giant-sized metamaterials that dampen the most destructive surface waves that are coupled to pressure and shear-body waves. Using large rings of a flexible metamaterials spaced around building foundations and tuned to the wavelength of the waves to be diverted, the researchers claimed the most destructive surface waves can be diverted around entire buildings.

Controlling the width, the chirality and refractive index of the waveguide, we can upgrade the aseismic range of the building as needed in order to defend it, at will, without touching it. It could be a big advantage of the waveguide method. This method will be effective for isolated buildings because we need some areas to construct the aseismic shell. It may be applicable for social overhead capitals such as power plants, dams, airports, nuclear reactors, oil refining complexes, long-span bridges, express rail-roads, etc.