Next Article in Journal
Modeling Temperature-Dependent Vibration Damping in C/SiC Fiber-Reinforced Ceramic-Matrix Composites
Next Article in Special Issue
The Multidirectional Auxeticity and Negative Linear Compressibility of a 3D Mechanical Metamaterial
Previous Article in Journal
Material Analysis of Steel Fibre Reinforced High-Strength Concrete in Terms of Flexural Behaviour. Experimental and Numerical Investigation
Previous Article in Special Issue
Planar Mechanical Metamaterials with Embedded Permanent Magnets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybridized Love Waves in a Guiding Layer Supporting an Array of Plates with Decorative Endings

1
IMSIA, ENSTA ParisTech, 828 Bd des Maréchaux, 91732 Palaiseau, France
2
Institut Langevin, CNRS, ESPCI ParisTech, 1 rue Jussieu, 75005 Paris, France
3
LAUM, CNRS UMR 6613, Le Mans Université, avenue Olivier Messiaen, 72085 Le Mans, France
4
Aix Marseille Univ., CNRS, Centrale Marseille, Institut Fresnel, 13013 Marseille, France
*
Author to whom correspondence should be addressed.
Materials 2020, 13(7), 1632; https://doi.org/10.3390/ma13071632
Submission received: 11 November 2019 / Revised: 20 December 2019 / Accepted: 27 December 2019 / Published: 1 April 2020
(This article belongs to the Special Issue Advances in Mechanical Metamaterials)

Abstract

:
This study follows from Maurel et al., Phys. Rev. B 98, 134311 (2018), where we reported on direct numerical observations of out-of-plane shear surface waves propagating along an array of plates atop a guiding layer, as a model for a forest of trees. We derived closed form dispersion relations using the homogenization procedure and investigated the effect of heterogeneities at the top of the plates (the foliage of trees). Here, we extend the study to the derivation of a homogenized model accounting for heterogeneities at both endings of the plates. The derivation is presented in the time domain, which allows for an energetic analysis of the effective problem. The effect of these heterogeneous endings on the properties of the surface waves is inspected for hard heterogeneities. It is shown that top heterogeneities affect the resonances of the plates, hence modifying the cut-off frequencies of a wave mathematically similar to the so-called Spoof Plasmon Polariton (SPP) wave, while the bottom heterogeneities affect the behavior of the layer, hence modifying the dispersion relation of the Love waves. The complete system simply mixes these two ingredients, resulting in hybrid surface waves accurately described by our model.

1. Introduction

The problem of waves propagating in an elastic half-space supporting an array of beams or plates is well known in seismology, where the site–city interaction aims at understanding the interaction of seismic waves with a set of buildings. Starting with the seminal work of Housner [1] (see also [2]), the site–city interaction has been intensively studied numerically [3,4,5] and analytically [6,7,8,9,10,11]. In this context, seismic shields, or metabarriers, have been considered using resonators buried in soil [12,13,14,15] or arrays of trees with a gradient in their heights [16,17,18]. More generally, this configuration is the elastic analog of a corrugated interface able to support surface waves, studied in acoustics [19] and in electromagnetism [20,21], where they are known as Spoof Plasmon Polaritons (SPPs). SPPs play a fundament role in the extraordinary transmission of long wavelength electromagnetic waves through metallic gratings [22,23] and have been studied intensively in the past twenty years for their potential applications in subwavelength optics, data storage, light generation, microscopy, and bio-photonics; see, e.g., [24]. Such similarities between surface waves in electromagnetism and elastodynamics fuel research in seismic metamaterials [25], as they lead to simplified models that see behind the tree that hides the forest [26].
To describe classical SPPs, the homogenization of a stratified medium is an easy and efficient tool [27,28]; the analysis is valid in the low frequency regime, namely owing to the existence of a small parameter measuring the ratio of the array spacing to the typical wavelength, and it provides, at the dominant order, the dispersion relation of SPPs. Thanks to the mathematical analogy between the problem in electromagnetism and in elasticity, this approach was applied in [18] accounting for the presence of a guiding soil layer underlain by an elastic half-space. Simple dispersion relations have been obtained from the effective model for the resulting spoof Love waves, so-called because of the characteristics they share with classical Love waves (surface waves supported by the layer on its own) and SPPs. Next, to account for the presence of heterogeneities (a foliage) at the top of the plates (trees), a hybrid model was used where the homogenization was performed locally (near the top of the plates) at the second order.
The present study generalizes and complements this study following two ways: (i) from a physical point of view, we include the effect of heterogeneities at the bottom of the plates (Figure 1), and (ii) from a technical point of view, we derive the full model at second order. This produces a significant increase in the accuracy of the theoretical prediction: in the reported examples, the model at order two is accurate up to a 1–2% error margin, while the model previously used in [18], at order one, would be accurate up to 10–30%. The second order model (see Equations (2) and (3)) provides a one-dimensional problem along the z-direction with a succession of homogeneous layers: the substrate occupying a half-space, the guiding layer, and an effective anisotropic layer replacing the region of the plates (see Figure 2). The effect of the heterogeneities at the bottom is encapsulated in transmission conditions, which tell us that the displacement and the normal stress are not continuous; this holds for plates without ending heterogeneities, a fact that was disregarded in [18]. The effect of the heterogeneities at the top is encapsulated in a boundary condition that differs from the usual stress free condition, as in [18]. We recover that for most of the frequencies, the plates do not interact efficiently with the layer; in the present case, it results that the surface wave resembles that of the layer only, hence a wave of the Love type. However, the resonances of the plates produce cut-off frequencies around which the dispersion relations are deeply affected. For simple plates, this can already produce drastic modifications in the dispersion relations (hybridization of the Love branches, avoided crossings at the cut-off frequencies of the plates). When heterogeneities at the endings of the plates are accounted for, additional changes happen. The heterogeneities at the bottom of the plates modify the behavior of the layer on its own, resulting in modified Love waves. The heterogeneities at the top of the plates modify the resonances of the plates, hence the cut-off frequencies. These two simple ingredients allow us to interpret qualitatively the various dispersion relations obtained in the configuration of the plates decorated at both ends. Next, the dispersion relations are accurately recovered by our homogenized model.
The paper is organized as follow. Section 2 summarizes the main results of the analysis: the effective model, Equations (2) and (3), and the resulting equation of energy conservation, Equation (10). The full derivation of the effective model is detailed in the Appendix A and Appendix B. In Section 3, we inspect the characteristics of waves guided by an array of decorated plates. The dispersion relations of these waves are exhibited numerically and compared to the closed forms provided by the effective model, Equations (21)–(23). The heterogeneities have the form of an additional thin hard layer at the bottom of the plates and a thin hard cap on the top. These simple shapes of heterogeneities allow us to discuss the Love waves modified by the bottom heterogeneity only and the resonances of a plate modified by the top heterogeneity.
Throughout the paper, we use the following notations:
-
Material properties: mass density ρ and shear modulus μ , with subscripts “S” for the substrate, “L” for the guiding Layer, “P” for the Plates, and “b,t” for the heterogeneities at the bottom and at the top of the plates.
-
Geometrical parameters: the layer has a total height H L = h L + h b with h b occupied by the heterogeneities. The array of plates is periodic with spacing , with plate thickness φ P h P and total height H P = h P + h t ( h t occupied by the heterogeneities). The heterogeneities at the bottom and top of the plates have surfaces S b = φ b h b and S t = φ t h t .

2. Summary of the Main Results

In the actual problem, the Navier equations for shear waves simplify to a wave equation for the antiplane displacement u = u y ( x , t ) and the stress vector σ ( x , t ) , of the form [29]:
σ ( x , t ) = μ ( x ) u ( x , t ) , ρ ( x ) 2 u t 2 ( x , t ) = d i v σ ( x , t ) ,
with x = ( x , z ) and t the time. The mass density ρ ( x ) and the shear modulus μ ( x ) are piecewise constant in the different materials, substrate/layer/plate/heterogeneities; see Figure 1. At each boundary between the elastic materials, the continuity of the displacement u and of the normal stress σ · n holds (with n the local normal vector). Eventually, at the boundaries separating elastic media and air, the stress-free boundary condition σ · n = 0 applies. In this section, we present the effective model deduced from the asymptotic analysis developed in Appendix A.

2.1. Effective Model

In the effective homogenized model, the regions of the substrate z ( , H L ) and of the guiding layer z ( H L , h b ) are kept as in the actual problem, while the region of the plates z ( 0 , h P ) is replaced by an equivalent homogeneous region of the same height. In this region, the medium is highly anisotropic, with propagation being allowed in the vertical direction z only; this calculation follows from [18] and applies almost identically in the acoustic case for arrays of Helmholtz resonators [30]. The boundary condition at the top of the effective medium, z = h P , is a condition of the Robin type for the normal stress. The transmission conditions at the bottom of the effective medium apply across the actual region of the heterogeneity, and they involve four parameters depending on the geometry of the heterogeneity and of the plates. Specifically, the homogenized model reads as:
for z ( , H L ) , σ = μ S u , ρ S 2 u t 2 = div σ , for z ( H L , h b ) , σ = μ L u , ρ L 2 u t 2 = div σ , for z ( 0 , h P ) , σ = μ P φ P 0 0 0 1 u , ρ P φ P 2 u t 2 = div σ ,
along with the continuity of u and σ · n at z = H L and the dynamic effective conditions:
across   the   region   ( h b , 0 ) , u = b μ L σ z ¯ + l b u ¯ x , σ z = l b σ z ¯ x μ L L b 2 u ¯ x 2 + h b ρ ^ b 2 u ¯ t 2 , at   the   top   of   the   plates   z = h P , σ z ( x , h P , t ) = L t σ z z ( x , h P , t ) .
The transmission conditions involve u = u ( x , 0 , t ) u ( x , h b , t ) and u ¯ = 1 2 u ( x , 0 , t ) + u ( x , h b , t ) , being the jump of u across the bottom heterogeneity and its mean value, respectively, and the same for σ z .
Among the five effective parameters ( b , l b , L b , ρ ^ b , L t ) entering in the effective conditions, two are known explicitly, while three are defined by elementary problems on ( V 1 , V 2 ) that satisfy static problems set in non-dimensional coordinate χ = ( χ , ζ ) = ( x / , z / ) in the vicinity of z = 0 (see Figure A3 in Appendix A.3). These problems read as:
div μ μ L V 1 = 0 , lim ζ V 1 = e z , lim ζ + V 1 = μ L φ P μ P e z , div μ μ L V 2 + χ = 0 , lim ζ V 2 = 0 , lim ζ + V 2 = e x ,
with V 1 , μ V 1 · n continuous at each interface between two elastic media and V 1 · n = 0 at the boundaries with the air and V 1 and μ V 1 one periodic with respect to χ for ζ < 0 (the same for V 2 and μ V 2 + χ ). Then, we have:
Effective parameters in ( 3 ) b = lim ζ + V 1 μ L φ P μ P ζ + h b , l b = lim ζ + ( V 2 + χ ) , L b = Y P μ P μ L χ ( V 2 + χ ) d χ + Y b μ μ L χ V 2 d χ + h b μ L φ b μ b + ( 1 φ b ) μ L , ρ ^ b = φ b ρ b + ( 1 φ b ) ρ L , L t = h t ρ t φ t ρ P φ P .
It is worth noting that the homogenized problem is set in a domain where the regions ( h b , 0 ) and ( h P , H P ) occupied by the heterogeneities have disappeared. It should be possible to extend the anisotropic region to ( 0 , H P ) as done in [18]; this would lead to a different, but as accurate effective model, with slightly different boundary condition at z = H P (specifically, a different value of L t ). However, this is not suitable from an energetic point of view (see Section 2.2). Similarly, the transmission conditions involve jumps of the displacement and of the normal stress across a non-zero interface. It should be possible to express the transmission conditions across a zero thickness interface located say at z = 0 . Again, this would lead to a different and as accurate effective model, with slightly different transmission conditions (with different values of b and L b ); again, our choice guaranties good properties of the energy in the effective problem.

2.2. Effective Energy

The solution ( u , σ ) of the homogeneous problem is expected to approximate the, say numerical, solution ( u num , σ num ) of the actual problem. Hence, we expect that the actual elastic energy is also correctly approximated in the effective problem. In the actual problem, the elastic energy simply reads as [29]:
E num = 1 2 D num 1 μ | σ num | 2 + ρ u num t 2 d x .
We shall now interrogate the equation of energy conservation in the homogenized problem where the effective boundary and jump conditions in (3) make additional energies appear. These terms appear primarily as fluxes within the bounded region D (see Figure 3), but they can be written as the time derivative of effective energies supported by the surface γ at the top of the plates and across the heterogeneities at the bottom of the plates ( Γ ± ).
By simple manipulation of the equations in (2), the equation of energy conservation in the homogenized problem is found to be of the form:
d d t E S + E L + E P + Φ = 0 ,
with:
E S , L = 1 2 D S , L | σ | 2 μ S , L + ρ S , L u t 2 d x , E P = 1 2 D P | σ z | 2 μ P φ P + ρ P φ P u t 2 d x ,
and:
Φ = D u t σ · n d l
(here, Φ is a line integral). The flux Φ has a contribution on Σ and two contributions that do not cancel even if the region D is bounded, that is if Σ is associated with Neumann or Dirichlet boundary conditions. Specifically, Φ = Φ Σ + Φ b + Φ t with:
Φ b , t = d d t E b , t and , E b = 1 2 Γ μ L L b u ¯ x 2 + h b ρ ^ b u ¯ t 2 + b μ L σ z ¯ 2 d x , E t = 1 2 γ ρ P φ P L t u t 2 d x ,
where n is the normal interior and D S , L , P the parts of D occupied by the substrate, the layer, and the plates, respectively. We have used that σ z = L t ρ P φ P t t u on γ from (2) and (3). We also have that Φ b = Γ t u ¯ σ z + t u σ z ¯ d x ; hence, Φ b = Γ l b x σ z ¯ μ L L b x x u ¯ + h b ρ ^ b t t u ¯ t u ¯ + b μ L t σ z ¯ + l b x t u ¯ σ z ¯ d x . The two terms in l b cancel after integration by parts of one of them, and we integrate also by parts the term in L b . It is worth noting that the integrations by parts make boundary terms (b.t.) appear. These terms can be interpreted as concentrated forces at the ending points of Γ ± along x; they are disregarded in the present study. Next, E b , t in (10) are energies since they are definite positive quadratic forms. Indeed, L t > 0 from (5) and ρ ^ b > 0 from (5), and it is shown in Appendix B that b and L b are positive as well. It is also worth mentioning that choosing a different position for γ would produce a different and possibly negative value of L t . Similarly, expressing the transmission conditions across a zero thickness interface would produce a possibly negative value of L b . Discussions on the effective energies can be found in [31,32].
We further stress that the homogenized problem is set on D , which differs from D num ; the regions D b for z ( h b , 0 ) and D t for z ( h P , H P ) are missing. Intuitively, we expect that the effective energies E b , t represent the elastic energies in D b and D t in the actual problem; specifically, we expect that:
E b , t 1 2 D b , t 1 μ | σ num | 2 + ρ u num t 2 d x .
We shall illustrate in Section 3.4 that these intuitive relations are indeed legitimate.

3. Hybrid Love Waves in a Guiding Layer Supporting Decorated Plates

In this section, we inspect the ability of the effective problem (2) and (3) to reproduce the scattering properties of an actual array. We consider the geometry of Figure 4: = 1 in arbitrary unit length, φ P = φ t = 0.5 and φ b = 1 . The total heights H P = h P + h t = 12 , H L = h L + h b = 8 are fixed. When the heterogeneities are considered, we set h t = 1 (hence, h P = 11 ) and/or h b = 1 (hence, h L = 7 ). We give in the tables below the material properties and the values of the effective parameters entering in the effective conditions (3).
We consider the time-harmonic regime with a time dependence e i ω t , which is omitted in the following, and inspect the solution of a scattering problem for a wave coming from z = with a wavenumber β along x, resulting in a reflected wave with a complex reflection coefficient R. This scattering problem allows us to cover the case of an incoming propagating wave, with | R | = 1 for β ω / c S , and the case of guided waves, when | R | = for β > ω / c S . The actual problem has to be solved numerically, and this was done using classical multimodal calculations.
In the rest of this section, we shall use for β the component of the wavenumber along x and make use of the following quantities:
k P = ω c P , γ L = ω 2 c L 2 β 2 , γ S = ω 2 c S 2 β 2 ,
( c a = μ a / ρ a for a = P, L, S).

3.1. Two Reference Solutions

To begin with, we establish two families of reference solutions that will be useful to analyze our problem. The first is that of Love waves supported by a guiding layer on the top of a substrate with c L < c S , which can be affected by the presence of the bottom heterogeneities. The second family is that of the Spoof Plasmon Polaritons (SPPs) in the plates, which can be affected by the presence of heterogeneities at the bottom of the plates.

3.1.1. Love Waves and Modified Love Waves

If we remove the array (Figure 5), the problem is reduced to a guiding layer sandwiched between air and the semi-infinite substrate (classical Love wave), and its modified version when a thin hard layer is added. The exact solutions of these problems are easily obtained. For classical Love waves, the solution of the scattering problem reads as:
u ( x , z ) = e i β x × A cos γ L z , z ( h L , 0 ) , e i γ S ( z + h L ) + R Love e i γ S ( z + h L ) , z ( , h L ) ,
and using the continuity of the displacement and of the normal stress provides ( A , R Love ) , in particular:
R Love = tan γ L h L i Y tan γ L h L + i Y , with Y = μ S γ S μ L γ L .
We recover the usual dispersion relation of Love waves for Y imaginary ( γ S imaginary with a positive imaginary part) and | R Love | = , which guaranties a family of Love wave dispersion branches in ω / c S < β < ω / c L ; see Figure 5.
If we add a layer of thickness h b in the guiding layer, the exact solution reads as:
u ( x , z ) = e i β x × A cos γ b z , z ( h b , 0 ) , B cos γ L z + C sin γ L z , z ( h L , h b ) , e i γ S ( z + h L ) + R Love b e i γ S ( z + h L ) , z ( , h L ) .
Again, applying the continuity of the displacement and of the normal stress at z = h L , h b provides ( A , B , C , R Love b ) and, in particular:
R Love b = tan γ L h L + Θ b i Y tan γ L h L + Θ b + i Y , Θ b = tan 1 μ b γ b μ L γ L tan γ b h b ,
where we have defined γ b = ω 2 c b 2 β 2 , c b = μ b / ρ b . Surface waves in this configuration have a dispersion relation tan γ L h L + Θ b + i Y = 0 , which can differ significantly from the dispersion relation of the classical Love waves; see Figure 5.

3.1.2. SPPs and Modified SPPs

The dispersion relation of spoof plasmons was derived using approximate methods [19,21] including classical homogenization [27,28]. However, the asymptotes at the cutoff frequencies can be straightforwardly calculated since they correspond to resonances of the plates associated with Dirichlet–Neumann boundary conditions at the bottom-top of the plate. For the classical SPPs, with k P = ω / c P , the solution simply reads as u ( x , z ) = A cos k P ( z H P ) e i β x for z ( 0 , H P ) , where we have anticipated the Neumann boundary condition (stress-free condition) at the top of the plate. At resonance, the Dirichlet boundary condition applies at z = 0 (the plate is clamped to the layer or to the substrate), resulting in the resonance frequency defining the asymptotes for:
Asymptotes of SPPs : cot an k P H P = 0 , hence ω n SPP = ( 2 n + 1 ) π 2 c P H P .
When the plate is terminated by a cap of the same thickness φ P and height h t (with k t = ω / c t ), the solution reads as u ( x , z ) = B cos k P z + C sin k P z e i β x for z ( 0 , h P ) and u ( x , z ) = A cos k t ( z H P ) e i β x for z ( h P , H P ) . Still at resonance, the Dirichlet boundary condition at z = 0 imposes B = 0 ; then applying the continuity of the displacement and normal stress at z = h P provides two relations on ( A , C ) , which are compatible if:
Asymptotes of modified SPPs : cot an k P h P = μ t k t μ P k P tan k t h t k P L t ,
and the last equality holds in the case where k t h t 1 with L t defined in (5) (and L t = 0 for h t = 0 ). In our geometry, with L t = 10 , H P = 12 , and c P = 240 , the first three asymptotes of the classical SPPs are obtained for ω n SPP / ( 2 π ) = 3 , 15 , 25 . In the presence of the caps of thickness h t = 1 (hence, h P = 11 ), solving the implicit relation cot an X = 10 11 X , with X = ω h P / c P , provides the first three modified asymptotes at ω n SPP , t / ( 2 π ) = 3.1 , 12.0 , 22.4 ; see Figure 6.

3.2. Dispersion Relation of Hybridized Love Waves

We shall see that the properties of the hybridized Love waves can be understood in light of the two ingredients studied in the previous section. On the one hand, for weak coupling between the plates and the guiding layer, our surface waves resemble Love waves, which are affected by the presence of heterogeneities. On the other hand, their dispersion relation is deeply modified in the vicinity of the SPPs asymptotes at cut-off frequencies (the coupling is maximum at those frequencies); the cut-off frequencies are dictated by the characteristics of the plates, hence sensitive to the presence of heterogeneities at their tops.
For the geometry of Figure 4, the solution of the scattering problem reads as:
u ( x , z ) = e i β x × A cos k P ( z h P ) + k P L t sin k P ( z h P ) , z ( 0 , h P ) , B cos γ L ( z + h b ) + C sin γ L ( z + h b ) , z ( H L , h b ) , e i γ S ( z + H L ) + R e i γ S ( z + H L ) , z ( , H L ) .
We have accounted for the boundary condition σ z = L t z σ z at z = h P (from (3)). Next, accounting for the continuities of the displacement and of the normal stress at x = H L and for the effective transmission conditions in (3) between x = h b and x = 0 provides the four relations needed to deduce ( A , B , C , R ) . This leaves us with:
R = D * ( ω , β ) D ( ω , β ) ,
where:
D ( ω , β ) = φ P μ P k P μ L γ L tan k P h P + k P L t + C b 1 k P L t tan k P h P 1 i Y tan γ L h L + 1 k P L t tan k P h P φ P μ P μ L k P b tan k P h P + k P L t tan γ L h L + i Y ,
and:
Y = μ S γ S μ L γ L , C b = h b ρ ^ b ω 2 μ L L b β 2 μ L γ L .
In (20), D * is deduced from D by substituting i Y by i Y . It results that for waves propagating in the substrate ( γ L and γ S real for β < ω / c S < ω / c L ), D * is the complex conjugate of D, and | R | = 1 , as expected. Next, surface waves correspond to γ S imaginary with a positive imaginary part, and | R | = ; hence, we deduce that:
Dispersion relation of guided waves : D ( ω , β ) = 0 , i γ S real > 0 .
It is worth noting that we recover the exact dispersion relations of the classical and modified Love waves for H P = 0 , whence h P = h t = 0 and L t = 0 in (5). We still have to determine C b in (22), and to do so, we have to adapt the elementary problem for V 2 to find L b in (5). This can be done easily by setting σ z 0 ( x , 0 , t ) = 0 in (A30) and replacing the limit to + by a boundary condition on χ = 0 ; it results that the limit at + of V 2 in (4) is replaced by V 2 = e x at ζ = 0 , and the integral over Y P cancels in (5); the problem is simpler, but still non-trivial. However, in the case where the heterogeneity in the bottom is a thin layer ( φ b = 1 ), the integral over Y b cancels as well by periodicity, and L b = h b μ b / μ L . Thus, we get D ( ω , β ) = C b 1 i Y tan γ L h L + tan γ L h L + i Y , with C b = μ b γ b 2 h b μ L γ L tan Θ b . Expectedly, we recover the exact dispersion relation of Love waves for h b = 0 ( C b = 0 ) and that of the modified Love waves in the limit of small h b .

3.3. Validation of the Homogenized Solution

Figure 7 shows the main results of the present study; we report the dispersion relations in four cases from the undecorated plates to the plates decorated at both endings; see the Table 3. They are visible by means of a large (diverging) reflection coefficient in R num computed in the direct numerics and from the explicit homogenized R value in (20)–(22).
The exact dispersion relations of the classical, (14), and modified, (16), Love waves are given for comparison. As previously said, the interaction of the plates with the layer is weak except in the vicinity of the cut-off frequencies. However, for the relatively tall plates that we have considered, these cut-off frequencies are sufficiently close to each other to modify the dispersion relation of our guided waves deeply.
This is already visible for the undecorated plates (Case 1); the guided waves tend to resemble the classical Love waves, but they experience several hybridizations at the cut-off frequencies of the SPPs in (17), accompanied by avoided crossings (sometimes not so pronounced). In the presence of the caps on the top of the plates, the same scenario is observed, with now cut-off frequencies given by the modified SPPs in (18). Eventually, Cases 3 and 4 with the thin hard bottom layer reproduce the same sequence as Cases 1 and 2 with the guided waves, which tend to resemble the modified Love waves.
The ability of our homogenized solution to reproduce the actual dispersion relations accurately is excellent, less than 2% on average in the reported ranges of ω and β , once | R | -values larger than 10 have been saturated. This is particularly visible in the zooms of Figure 8 in the vicinity of avoided crossings for the decorated plates (AC1 and AC2 in Figure 7). In particular, we stress that in the absence of heterogeneity, we have a stress-free condition since L t = 0 in (5), but we do not have the continuity of the displacement and normal stress across z = 0 since b and L b do not vanish (see Table 1). In comparison, the homogenized model (see Table 2 for homogenized coefficients) at the leading order provides the usual stress-free condition and continuity relations regardless of the presence of heterogeneities. Thus, it misses the effects of the heterogeneities and provides the same prediction for Cases 1 to 4. The resulting error is significantly higher, about 10% for Cases 1 and 2, and about 30% for Cases 3 and 4.
Figure 7. Dispersion relations of guided waves in four configurations of plates from direct numerics and from the homogenized solution (20)–(22); the dispersion relations are revealed by large | R | values. The exact reference dispersion relations of classical and modified Love waves are given in the top panel for comparison. Dotted lines are a guide for the eye showing light lines for Love waves and asymptotes for SPPs. AC1 and AC2 are avoided crossings magnified in Figure 8. Details of homogenized coefficients and geometrical parameters are given in Table 2 and Table 3.
Figure 7. Dispersion relations of guided waves in four configurations of plates from direct numerics and from the homogenized solution (20)–(22); the dispersion relations are revealed by large | R | values. The exact reference dispersion relations of classical and modified Love waves are given in the top panel for comparison. Dotted lines are a guide for the eye showing light lines for Love waves and asymptotes for SPPs. AC1 and AC2 are avoided crossings magnified in Figure 8. Details of homogenized coefficients and geometrical parameters are given in Table 2 and Table 3.
Materials 13 01632 g007
Eventually, we report in Figure 9 examples of the displacement fields of the guided wave at ω / ( 2 π ) = 24 for the arrays of undecorated plates and of decorated plates. Both in the numerics and in the homogenized problem (19), the whole solutions were divided by R num and R respectively in order to produce an evanescent wave in the substrate of the form e i β x e | γ S | ( x + H L ) that has the unitary amplitude at z = H L . This allows for quantitative comparison between the direct numerics and the homogenized solution without any adjustable parameter, see Table 2. The agreement in the strengths of the resonances and in the repartition of the amplitudes in the substrate, in the layer, and in the plates is again excellent. At the reported frequency, the wave is evanescent in the substrate as soon as β > 0.15 , and it becomes evanescent in the layer for β > 0.75 ; this is visible for the undecorated plates for the guided waves with β 2 = 0.77 , which is supported by the array only.

3.4. Energies in the Actual/Homogenized Problems

In this section, we inspect the intuitive relations announced in (11), namely that the effective energies E b , t coincide with the elastic energies stored in the regions of the heterogeneities.
In the actual problem, we define D num = { x ( 0 , 1 ) , z ( H * , 0 ) } { x ( 0 , φ P ) , z ( 0 , H P ) } . The energies in the actual problem are the usual elastic energies, which read, in the succession of regions, substrate, layer, bottom, plate, and top, as:
E num , A = 1 2 D A μ num , A | u num | 2 + ρ num , A ω 2 | u num | 2 d x ,
where A = S, L, b, P, t and D S = { x ( 0 , 1 ) , z ( H * , h L ) } , D L = { x ( 0 , 1 ) , z ( h L , h b ) } , D b = { x ( 0 , 1 ) , z ( h b , 0 ) } , D P = { x ( 0 , φ P ) , z ( 0 , h P ) } , D t = { x ( 0 , φ P ) , z ( h P , H P ) } .
In the effective problem, the energies are obtained explicitly owing to the solution in (19), which provides the fields in the substrate, layer, and effective region accounting for the plates. We denote u ( x , z ) = f S , P , L ( z ) e i β x in (19); hence, f S ( z ) = e i γ S ( z + H L ) + R e i γ S ( z + H L ) , f L = B cos γ L ( z + h b ) + C sin γ L ( z + h b ) , and f P ( z ) = A cos k P ( z h P ) + k P L t sin k P ( z h P ) , R given by (20)–(22), and:
A = 2 i Y D ( ω , β ) cos k P h P cos γ L h L , B = 2 i Y D ( ω , β ) cos γ L h L 1 k P L t tan k P h P φ P μ P μ L k P b ( tan k P h P + k P L t ) , C = 2 i Y D ( ω , β ) cos γ L h L φ P μ P k P μ L γ L tan k P h P + k P L t + C b 1 k P L t tan k P h P .
It follows that the effective energies in (8) read as:
E S = H * h L μ S | f S | 2 + ρ S ω 2 + μ S β 2 | f S | 2 d z , E L = h L h b μ L | f L | 2 + ρ L ω 2 + μ L β 2 | f L | 2 d z , E P = 0 h P μ P φ P | f P | 2 + ρ P φ P ω 2 | f P | 2 d z ,
and for the last integral, we accounted for the effective stress–displacement relations in (2). Next, from (10) along with (3), it is easy to see that:
E b = 2 μ L L b β 2 | f P ( 0 ) | 2 + h b ρ ^ b ω 2 | f P ( 0 ) | 2 + b μ L ( μ P φ P k P ) 2 | f P ( 0 ) | 2 , E t = 2 ρ P φ P L t ω 2 | f P ( h P ) | 2 .
We computed the energies in the actual problem, (24), and in the homogenized problem, (26) and (27), for an incident propagating wave ( γ S real in (19)). The real part of the reflection coefficient R is reported in the left panel of Figure 10. In the case of weak coupling with the array of plates, R R Love b , (16); hence, R 1 , except in the vicinity of the resonances of the layer (diverging tan γ L h L + Θ b ) where it goes to −1; see light grey arrows. Next, strong coupling with the array occurs at the resonance of the plates, resulting in R 1 ; see the dark grey arrows. The resulting repartition of the energies is plotted against the frequency for a wave at incidence 45 (we normalized the energies to the total energy). For ω / ( 2 π ) ( 030 ) , three resonances of the plates and two resonances of the layer take place. Expectedly, the energy in the plates E P is small except at the resonance of the plates where almost all the energy is shared in the plates and their top heterogeneities. This is particularly visible at the first resonance where 35% of the total energy is supported by the heterogeneities. Symmetrically, at the resonances of the layer, most of the energy is supported by the layer ( E L ) and the bottom heterogeneities ( E t ). The ability of the homogenized solution to reproduce the solution in the substrate, in the layer, and in the plates is recovered in the energies with error margins of 0.1%, 0.5%, and 4%, respectively, in the reported case. More remarkably, the effective energies E t and E b accurately reproduce the variations of the actual elastic energies, with error margins of 0.7%, which legitimizes the intuitive relations (11).
Eventually, we report in Figure 11 the surface density of energy e computed numerically ( e = 1 2 μ num , A | u num | 2 + ρ num , A ω 2 | u num | 2 for A = S, L, b, P, t) and that of the effective problem (which varies with z only from (26)); it is worth noting that the energies E t , b do not give rise to surface density since they are defined along lines. We recover the observations of Figure 10: at a resonance of the layer ( ω / ( 2 π ) = 19 ) , almost all the energy is stored in the layers; at a resonance of the plates ( ω / ( 2 π ) = 19 ) , it is stored in the plates; and otherwise, it is equally distributed.

4. Concluding Remarks

We studied the problem of wave propagation in a geometry that combined two resonators, a soft layer in a substrate and an array of plates. In particular, we focused on the ability of thin heterogeneities at the endings of the plates to impact on the response of the system significantly. This was done thanks to asymptotic homogenization accounting for the boundary effects to be captured at the endings of the plates (e.g. foliages and roots for a model of trees). Such an analysis provides a simple effective model in which the region of the plates are replaced by a homogeneous highly anisotropic region and the effects of the heterogeneities were encapsulated in effective dynamic conditions. It was shown that the resulting effective model accurately predicts the dispersion relation of surface waves far beyond the quasi-static limit. These anti-plane shear waves share common features with Love waves in geophysics and surface plasmon polaritons (SPPs) in photonics; the dispersion relation of such hybridized Love waves was obtained in a closed form that allowed us to discriminate the role of the layer and that of the plates. Besides, we showed that the presence of heterogeneities at the decorative endings of the plates may affect significantly the characteristics of the surface waves. We also showed that the variations of the actual elastic energies (in the different regions) are accurately reproduced by the effective energies identified in the actual problem. In particular, the contributions of the effective interface and of the effective surface correspond to the actual energies stored in the thin regions containing the heterogeneities, and they disappear in the effective problem.
Our approach is useful for at least two reasons. On the one hand, it provides a simpler problem for which explicit solutions are available. This was illustrated in the present study where the model was shown to be very accurate up to frequencies corresponding to a ratio of the wavelength to the spacing close to one (typically in the layer and in the plates), where the long wavelength homogenization failed. Next, numerical resolution in the time domain may become intractable due to the separation of the scales associated with the typical wavelength, the array spacing, and the possible thinner scales in the heterogeneities. Eventually, the interest in dealing with effective problems was exemplified for scalar waves in a 2D geometry; it would be all the more evident for polarized elastic waves in 3D geometries. We finally note that our approach is well adapted to handle substrates with a gradient in elastic properties such as granular media [33,34].

Author Contributions

All authors contributed to the analysis and discussed the results. All authors contributed to the final version of the manuscript. All authors read and agreed to the published version of the manuscript.

Acknowledgments

We acknowledge J.-J. Marigo for fruitful discussion.SG wishes to thank the Department of Mathematics at Imperial College London for a visiting position in the group of R.V. Craster in 2018–2019.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Asymptotic Analysis

In the asymptotic analysis that we shall conduct, we use the macroscopic (usual) coordinate x = ( x , z ) , and we introduce the so-called microscopic coordinate χ = ( χ , ζ ) , defined by:
χ = x ,
where = ε is a small positive parameter (and small means small compared to the typical, or maximal, wavenumber imposed by the source of order unity). The analysis is firstly conducted in the region of the plates far from their endings. Afterwards specific analyses are conducted in the vicinity of plate endings to account for the boundary layer effects.

Appendix A.1. The Homogenized Wave Equation

In the region of the plates, the wavefields vary over long distances in the two directions x and z; these long distance variations are accounted for by x . Next, short distance variations of the fields occur within a single plate, which are accounted for thanks to the additional coordinate χ ; note that ζ is not needed since the plates are invariant along z (Figure A1).
Figure A1. Homogenization in the bulk of the array. In the x = ( x , z ) coordinate, the array has a spacing = ε . The rescaling in the horizontal χ = x / ε Y P coordinate is shown in the inset; Y P is a one-dimensional domain.
Figure A1. Homogenization in the bulk of the array. In the x = ( x , z ) coordinate, the array has a spacing = ε . The rescaling in the horizontal χ = x / ε Y P coordinate is shown in the inset; Y P is a one-dimensional domain.
Materials 13 01632 g0a1
Thus, the differential operator reads as:
e x ε χ + x ,
which will be used in the plates where:
σ = μ P u , ρ P 2 u t 2 = div σ ,
apply, owing to the expansions:
u = u 0 ( x , χ , t ) + ε u 1 ( x , χ , t ) + , σ = σ 0 ( x , χ , t ) + ε σ 1 ( x , χ , t ) + ,
with χ Y P . During the homogenization procedure, the coordinate χ aims to disappear, and we shall see that the homogenized wave equation in (2) involves the effective fields u n and σ n defined by:
u n ( x , t ) 1 φ P Y P u n ( x , χ , t ) d χ , σ n ( x , t ) Y P σ n ( x , χ , t ) d χ .
Doing so, we anticipate the macroscopic equilibrium of the forces by implicitly extending the stress by zero in Y \ Y P , with Y = { χ ( 1 / 2 , 1 / 2 ) } .
We start the analysis at the leading order in 1 / ε , with u 0 χ = σ x 0 χ = 0 ; hence u 0 ( x , t ) and σ x 0 ( x , t ) do not depend on χ . It follows that σ x 0 is constant in Y P , and as it vanishes at the boundaries with air at χ = ± φ P / 2 , it is zero everywhere in Y P . Thus, we have:
σ x 0 ( x , t ) = 0 , u 0 ( x , t ) .
At the order ε 0 and accounting for (A5), we have:
σ z 0 ( x , χ , t ) = μ P u 0 z ( x , t ) , ρ P 2 u 0 t 2 ( x , t ) = σ z 0 z ( x , t ) + σ x 1 χ ( x , χ , t ) , in Y P .
It follows that σ z 0 ( x , t ) does not depend on χ and σ x 1 ( x , t ) = 0 (as for σ x 0 ), and thus:
σ z 0 ( x , t ) = μ P φ P u 0 z ( x , t ) , ρ P φ P 2 u 0 t 2 ( x , t ) = σ z 0 z ( x , t ) ,
by simply using that:
σ z 0 ( x , t ) = φ P σ z 0 ( x , t ) , σ x 1 = 0 .
We now move on to the second order. Starting with (A5), hence 0 = σ x 0 = μ P x u 0 ( x , t ) + χ u 1 ( x , χ , t ) , we deduce the displacement u 1 of the form:
u 1 ( x , χ , t ) = χ u 0 x ( x , t ) + u 1 ( x , t ) , in Y P ,
with the origin of χ such that χ = 0 . It follows that:
σ z 1 ( x , χ , t ) = μ P u 1 z ( x , χ , t ) = μ P χ 2 u 0 z x ( x , t ) + u 1 z ( x , t ) ,
which, after integration over Y P and thanks to χ = 0 , leaves us with:
σ z 1 ( x , t ) = μ P φ P u 1 z ( x , t ) .
Eventually, the equation of equilibrium at order one reads as:
ρ P 2 u 1 t 2 = d i v x σ 1 + σ x 2 χ , in Y P ,
which after integration over Y P gives:
ρ P φ P 2 u 1 t 2 = σ z 1 z ( x , t ) .
The effective wave equations at the leading order (A7) Â and at order one (A11)–(A13) have the same forms; hence, up to O ( ε 2 ) , ( u , σ ) satisfies the effective equation announced in (2).

Appendix A.2. The Boundary Condition at the Top of the Plates

The homogenized wave equation derived in the previous section has to be supplied with a boundary condition at the top of the plates. To derive this condition, we have to analyze the near (or evanescent) field excited in the vicinity of z = h P (Figure A2). In this region of small extent along z, the long distance variations of the macroscopic fields occur along x only (across the plates). Next, to describe the short distance variations of the evanescent field, we use the microscopic coordinates χ and:
ζ = z h P ε .
Accordingly, we consider the following asymptotic expansions:
u = v 0 ( x , χ , t ) + ε v 1 ( x , χ , t ) + , σ = τ 0 ( x , χ , t ) + ε τ 1 ( x , χ , t ) + ,
with χ Y = Y t Y P where Y t is the bounded region containing the heterogeneity (of vertical extent h t / and | Y t | = φ t h t / ), Y P = Y P × ( , 0 ) is the unbounded region of the plate since in rescaled coordinate χ , and the bottom of the plate has been sent to .
Figure A2. Elementary cell at the top of the plate in the χ = ( χ , ζ ) coordinate, in the unbounded two-dimensional region Y = Y P + Y t ; the region of the plate is Y P = Y P × ( , 0 ) , and the heterogeneity Y t is bounded in Y × ( 0 , h t ) , with | Y t | = φ t h t / .
Figure A2. Elementary cell at the top of the plate in the χ = ( χ , ζ ) coordinate, in the unbounded two-dimensional region Y = Y P + Y t ; the region of the plate is Y P = Y P × ( , 0 ) , and the heterogeneity Y t is bounded in Y × ( 0 , h t ) , with | Y t | = φ t h t / .
Materials 13 01632 g0a2
In Y , the fields ( u , σ ) satisfy:
σ = μ ( χ ) u , ρ ( χ ) 2 u t 2 = div σ ,
with ( μ ( χ ) , ρ ( χ ) ) varying within Y t depending on the characteristics of the heterogeneity and being equal to ( μ P , ρ P ) in Y P . The above system is complemented by a condition of zero normal stress at the boundaries in contact with air and conditions of continuity of displacement and normal stress at the interfaces between two elastic media. Eventually, boundary conditions are missing when ζ ; these boundary conditions are obtained by imposing that the fields in (A15) match those defined in (A3), which hold far from the top of the plate. This is written in an intermediate region where z h P and ζ . Using that z = h P + ε ζ in (A3) and re-expanding in Taylor expansions for small ε , we get the so-called matching conditions at each order. At the first and second orders and accounting for the fact that u 0 ( x , t ) and σ 0 ( x , t ) do not depend on χ from (A5) and (A8), we get:
u 0 ( x , h P , t ) = lim ζ v 0 ( x , χ , t ) , σ 0 ( x , h P , t ) = lim ζ τ 0 ( x , χ , t ) , u 1 ( x , h P , χ , t ) = lim ζ v 1 ( x , χ , t ) ζ u 0 z ( x , h P , t ) , σ 1 ( x , h P , χ , t ) = lim ζ τ 1 ( x , χ , t ) ζ σ 0 z ( x , h P , t ) .
We can now start the analysis, using e x x + 1 ε χ in (A16) along with (A15). At the leading order, we have d i v χ τ 0 = 0 , which after integration over Y and accounting for the boundary conditions and for the matching condition on σ 0 leaves us with:
0 = lim ζ Y P τ z 0 ( x , χ , t ) d χ = σ z 0 ( x , h P , t ) .
At the leading order, the effective boundary condition is the usual stress-free condition, regardless of the presence of the heterogeneity. Hence, we move to the next order to get the boundary condition on σ z 1 , and to do so, we have to determine τ 0 . We start with χ v 0 = 0 ; hence, v 0 ( x , t ) is independent of χ ; from the matching condition (A17) on u 0 , we get that:
v 0 ( x , t ) = u 0 ( x , h P , t ) .
This allows us to define the problem satisfied by ( v 1 , τ 0 ) in χ coordinate, which reads as:
d i v χ τ 0 = 0 , τ 0 ( x , χ , t ) = μ ( χ ) u 0 x ( x , h P , t ) e x + χ v 1 ( x , χ , t ) , in Y , v 1 , τ 0 · n , continuous   at   the   interfaces   between   two   elastic   media , τ 0 · n = 0 , at   the   boundaries   in   contact   with   air , lim ζ τ 0 = 0 .
For the limit ζ , we used in the matching condition (A17) for τ 0 that σ 0 ( x , h P , t ) = 0 from (A5) and (A18). It is easy to check that the system (A20) has an explicit solution, which reads as:
τ 0 = 0 , v 1 ( x , χ , t ) = χ u 0 x ( x , h P , t ) + v ^ 1 ( x , t ) ,
where v ^ 1 ( x , t ) does not need to be specified, but it appears since v 1 is defined in (A20) up to a function of ( x , t ) . Owing to the above results, the equation of equilibrium at the order ε 0 in (A16), specifically ρ 2 v 0 t 2 = d i v χ τ 1 + τ x 0 x = 0 , simplifies because of (A19) and (A21). After integration over Y and using that τ 1 · n is either zero (on the boundaries with the air) or continuous, we get:
2 u 0 t 2 ( x , h P , t ) Y ρ ( χ ) d χ = Y d i v χ τ 1 d χ = lim ζ Y P τ z 1 ( x , χ , t ) d χ .
Integrating over Y P the matching condition for σ z 1 in (A17), along with z σ z 0 = ρ P φ P t t u 0 from (A7), we find that (A22) can be written as:
σ z 1 ( x , h P , t ) = σ z 0 z ( x , h P , t ) lim ζ ζ + Y ρ ( χ ) ρ P φ P d χ .
The integral of the mass density over Y is Y ρ ( χ ) d χ = Y t ρ ( χ ) d χ lim ζ ζ ρ P φ P , from which:
σ z 1 ( x , h P , t ) = Y t ρ ( χ ) ρ P φ P d χ σ z 0 z ( x , h P , t ) .
The above expression is valid if one considers a heterogeneity with varying mass density, and in the case of uniform mass density ρ ( χ ) = ρ t , it simplifies to:
σ z 1 ( x , h P , t ) = h t ρ t φ t ρ P φ P σ z 0 z ( x , h P , t ) .
Eventually, making use of σ z 0 ( x , h P , t ) = 0 in (A18) and of (A25), we get that σ z 0 + ε σ z 1 (with ε = ); hence σ z up to O ( ε 2 ) satisfies the boundary condition announced in (3).

Appendix A.3. Jump Conditions

The analysis of the problem in the vicinity of the heterogeneity at the bottom of the plates (Figure A3) is similar to that conducted in the previous section. The elementary cell in which the analysis is conducted is unbounded for ζ ± since in the microscopic coordinate, the top of the plate and the lower interface of the layer were set to infinity. We use the same expansions as in (A15), and for simplicity, we keep the same notations, with χ = ( χ , ζ ) and:
χ = x ε , ζ = z ε .
In the present case, the elementary cell involves a part of ζ < 0 where we impose that ( v n , τ n ) in (A15) are periodic with respect to χ Y .
Figure A3. Elementary cell at the bottom of a single plate in the χ = ( χ , ζ ) coordinate. Y = Y P Y b with Y P = { χ Y P × ( 0 , + ) } , and Y b = { χ Y × ( , 0 ) } ; φ b h b / is the non-dimensional surface of the heterogeneity of vertical extent h b / .
Figure A3. Elementary cell at the bottom of a single plate in the χ = ( χ , ζ ) coordinate. Y = Y P Y b with Y P = { χ Y P × ( 0 , + ) } , and Y b = { χ Y × ( , 0 ) } ; φ b h b / is the non-dimensional surface of the heterogeneity of vertical extent h b / .
Materials 13 01632 g0a3
As in the previous section, the matching conditions tell us that the solution in Y when ζ ± matches the solution valid far from the bottom of the plates z ( h b , 0 ) . When ζ + , it matches (A3), which holds in the plates; when ζ , it matches the solution in the layer. As the layer is a homogeneous region, the expansion of the solution ( u , σ ) is trivial: of the same form as in (A3) with all the terms ( u n , σ n ) being a function of ( x , t ) only. We thus get that:
u 0 ( x , 0 ± , t ) = lim ζ ± v 0 ( x , χ , t ) , σ 0 ( x , 0 ± , t ) = lim ζ ± τ 0 ( x , χ , t ) , u 1 ( x , 0 , t ) = lim ζ v 1 ( x , χ , t ) ζ u 0 z ( x , 0 , t ) , u 1 ( x , 0 + , χ , t ) = lim ζ + v 1 ( x , χ , t ) ζ u 0 z ( x , 0 + , t ) , σ 1 ( x , 0 , t ) = lim ζ τ 1 ( x , χ , t ) ζ σ 0 z ( x , 0 , t ) , σ 1 ( x , 0 + , χ , t ) = lim ζ + τ 1 ( x , χ , t ) ζ σ 0 z ( x , 0 + , t ) .
For the limit of ( v n , τ n ) ( n = 0 , 1 ) when ζ , we used that all the ( u n , σ n ) do not depend on χ in the layer. For the limit of ( v n , τ n ) ( n = 0 , 1 ) when ζ + , we used that ( u 0 , σ 0 ) do not depend on χ in the plates from (A5) and (A8), but ( u 1 , σ 1 ) do.
We can now start the analysis, using as in the previous section, the differential operator e x x + 1 ε χ and the expansions (A15) in:
σ = μ ( χ ) u , ρ ( χ ) 2 u t 2 = div σ ,
with ( μ ( χ ) , ρ ( χ ) ) varying within Y b = { χ Y × ( , 0 ) } and being equal to ( μ P , ρ P ) in Y P = { χ Y P × ( 0 , + ) } . As previously, the leading order starts with χ v 0 = 0 ; hence: v 0 is independent of χ , and from (A27), v 0 ( x , t ) = u 0 ( x , 0 ± , t ) . Next, integrating the relation d i v χ τ 0 = 0 over Y along with the (A27) leaves us with:
v 0 ( x , t ) = u 0 ( x , 0 , t ) , σ z 0 0 = u 0 0 = 0 ,
where we have defined w 0 = w ( x , 0 + , t ) w ( x , 0 , t ) at this stage. At the leading order, the continuity of the displacement and of the normal stress apply regardless of the presence of the heterogeneities. Thus, we move to the next order, and as in the previous section, we need to define the problem on ( v 1 , τ 0 ) , which reads as:
d i v χ τ 0 = 0 , τ 0 = μ ( χ ) χ v 1 + u 0 x ( x , 0 , t ) e x , in Y , v 1 , τ 0 · n , continuous   at   the   interfaces   between   two   elastic   media , τ 0 · n = 0 , a t χ = ± φ P / 2 , in Y P , v 1 , τ 0 , one   periodic   with   respect   to   χ in Y b , lim ζ χ v 1 = σ z 0 ( x , 0 , t ) μ L e z , lim ζ + χ v 1 = σ z 0 ( x , 0 , t ) μ P φ P e z u 0 x ( x , 0 , t ) e x .
To find the above limits for ζ ± , we used the matching conditions for σ 0 in (A27) along with σ 0 ( x , 0 , t ) = σ z 0 ( x , 0 , t ) e z + μ L x u 0 ( x , 0 , t ) e x (since in the layer, σ z 0 = σ z 0 and σ x 0 = μ L x u 0 ) and along with σ 0 ( x , 0 + , t ) = σ z 0 ( x , 0 , t ) e z / φ P from (A5) and (A8).
The system (A30) is the counterpart of the system (A20) that we obtained at the top of the plate, but now, the solution is not trivial; hence, (A20) has to be solved numerically. However, instead of solving (A30) for a given scattering problem, that is for given external loadings ( x u 0 ( x , 0 , t ) , σ z 0 ( x , 0 , t ) ) , we shall use that (A30) is linear with respect to those loadings. Specifically, we set:
v 1 = 1 μ L σ z 0 ( x , 0 , t ) V 1 ( χ ) + u 0 x ( x , 0 , t ) V 2 ( χ ) + v ^ 1 ( x , t ) , τ 0 = μ ( χ ) μ L σ z 0 ( x , 0 , t ) χ V 1 ( χ ) + μ ( χ ) u 0 x ( x , 0 , t ) χ V 2 ( χ ) + χ ,
and it is sufficient that ( V 1 , V 2 ) satisfy (4) to ensure that v 1 satisfies (A30). The advantage is obvious: as ( V 1 , V 2 ) satisfy static problems, they can be computed once and for all, independently of the scattering problem that will be considered afterwards. We shall see that these elementary solutions provide the coefficients defined in (5).
From (4), V 1 has a linear behavior in ζ for ζ ± , and V 2 is linear in χ when ζ + . As V 1 , V 2 in (4) are defined up to a constant, we set the constant at zero for ζ and define:
lim ζ ( V 1 ζ ) = 0 , lim ζ V 2 = 0 , lim ζ + V 1 μ L φ P μ P ζ = α 1 , lim ζ + ( V 2 + χ ) = α 2 .
The jump of u 1 is obtained using the matching conditions (A27) for u 1 , with v 1 in (A31) along with (A32). Using in addition that σ z 0 = μ P φ P z , which holds at z = 0 + from (A7), and σ z 0 = σ z 0 = μ L z u 0 , which holds at z = 0 , we easily get that:
u 1 ( x , 0 , χ , t ) = v ^ 1 ( x , t ) , and   u 1 ( x , 0 + , χ , t ) = α 1 μ L σ z 0 ( x , 0 , t ) + ( α 2 χ ) u 0 x ( x , 0 , t ) + v ^ 1 ( x , t )
(incidentally, we recover that u 1 is linear w.r.t. χ as in (A9)). After integration over Y P (with χ = 0 ), we get:
u 1 0 = α 1 μ L σ z 0 ( x , 0 , t ) + α 2 u 0 x ( x , 0 , t ) .
To get the jump on σ z 1 , we use the equation of equilibrium in (A16) at order zero integrated over Y , specifically:
2 u 0 t 2 ( x , 0 , t ) Y ρ ( χ ) d χ = Y d i v χ τ 1 + τ x 0 x d χ
where we use that v 0 = u 0 ( x , 0 , t ) from (A29). We use (A27) along with (A31), and we account for the fact that Y b χ V 2 d χ and Y P χ ( V 2 + χ ) d χ are bounded and that Y P d χ = φ P lim ζ + ζ and Y b d χ = lim ζ + ζ . We also use that ρ P t t u 0 = z σ z 0 in the plate (from (A7) and (A8), which holds for z = 0 + ) and that ρ L t t u 0 μ L x x u 0 z σ z 0 = 0 in the layer (from (1), which holds for z = 0 ). Eventually, for constant ρ = ρ b and μ = μ b in the heterogeneity, we use that Y b ρ ρ L d χ = φ b h b / ( ρ b ρ L ) , Y b μ μ L d χ = φ b h b / ( μ b μ L ) . This leaves us with the jump in σ z 1 of the form:
σ z 1 0 = β 1 σ z 0 x ( x , 0 , t ) μ L β 2 2 u 0 x 2 ( x , 0 , t ) + φ b h b ( ρ b ρ L ) 2 u 0 t 2 ( x , 0 , t ) ,
with:
β 1 = Y μ μ L V 1 χ d χ , β 2 = Y P μ μ L ( V 2 + χ ) χ d χ + Y b μ μ L V 2 χ d χ + μ b μ L 1 φ b h b .
We shall express the jump conditions obtained in (A29) and in (A34)–(A36) in a different form, but equivalent up to O ( ε 2 ) . Specifically, we want the jump between h b and zero (we shall comment on this choice later on). With = ε and h b = O ( ) , we can use the Taylor expansion of w 0 ( x , 0 , t ) = w 0 ( x , h b , t ) + h b z w 0 ( x , 0 , t ) + O ( ε 2 ) for w 0 = u 0 , σ z 0 to get the jump of w = w 0 + ε w 1 defined as w = w ( x , 0 + , t ) w ( x , h b , t ) (as in (3)). From (A29) along with (A34) Â and (A36), we get:
u 0 + ε u 1 = α 1 + h b μ L σ z 0 ( x , 0 , t ) + α 2 u 0 x ( x , 0 , t ) + O ( 2 ) ,
where we used that σ z 0 = σ z 0 = μ L z u 0 . Replicating this for the jump in the normal stress, we get:
σ z 0 + ε σ z 1 = β 1 σ z 0 x ( x , 0 , t ) μ L h b + β 2 2 u 0 x 2 ( x , 0 , t ) + h b φ b ρ b + ( 1 φ b ) ρ L 2 u 0 t 2 ( x , 0 , t ) + O ( 2 ) ,
where we used that z σ z 0 = z σ z 0 = μ L x x u 0 + ρ L t t u 0 (and it will be shown in Appendix B.1 that β 1 = α 2 ). Eventually, we define w ¯ = 1 2 w ( x , h b , t ) + w ( x , 0 + , t ) ; with w 0 ( x , 0 , t ) = w ¯ + O ( 2 ) in the right hand side terms, the jumps in (A38) and (A39) are equivalent to that written in (3) omitting O ( 2 ) , and with:
b = α 1 + h b , l b = α 2 , L b = h b + β 2 ,
according to (5) along with (A32) and (A37). We made a choice on the expression of the jump conditions, and we shall see that this guaranties that the energy in the effective problem is a definite positive quadratic form. Note that we also made a choice on the position of the effective boundary condition on the top of the plate; it is possible to choose a different position in the vicinity of the order of of that one; however, as already stressed, if all the resulting effective problems are equivalent up to O ( 2 ) , all of them do not guaranty proper energy of the effective problem; see, e.g., [31,32].

Appendix B. Properties of the Effective Parameters

Here, we shall prove three properties that have been used in the previous Appendix: (1) In (A39), we used that α 2 + β 1 = 0 . Furthermore, the effective energy E eff in (10) is a definite positive quadratic form if: (2) b > 0 and (3) L b > 0 , with ( b , L b ) defined in (A40).

Appendix B.1. α2 + β1 = 0

The parameters ( α 2 , β 1 ) are defined by α 2 = lim ζ + ( V 2 + χ ) (see (A32)) and β 1 = Y μ μ L V 1 χ d χ (see (A37)), with V 1 and V 2 satisfying the elementary problems: (4) (note that b = α 2 in (5)). We start with:
0 = Y V 2 div μ μ L V 1 d χ = Y μ μ L V 1 · V 2 d χ + α 2 ,
where we used that all the boundary terms μ V 2 V 1 · n on Y vanish except at ζ + where Y P μ μ L V 2 V 1 · e z d χ = Y P 1 φ P ( α 2 χ ) d χ = α 2 since μ = μ P and χ = 0 . Indeed, V 2 and μ V 1 · n are continuous at the interface between two elastic media; next, μ V 1 · n = 0 at the boundaries in contact with air; eventually, for ζ , V 2 vanishes. Next, considering:
0 = Y V 1 div μ μ L V 2 + χ d χ = Y μ μ L V 1 · V 2 d χ β 1 ,
and here, all the boundary terms vanish since μ ( V 2 + χ ) · n is continuous or vanish and for ζ ± ( V 2 + χ ) · ( ± e z ) = 0 . It follows that:
α 2 = β 1 = Y d χ μ μ L V 1 · V 2 .

Appendix B.2. ℓb > 0

To show that b > 0 or equivalently α 1 + h b / 0 from (A40), we rely on the variational formulation of the elementary problem on square integrable field V = V 1 H with H ( χ , ζ < 0 ) = ζ and H ( χ , ζ > 0 ) = μ L φ P μ P ζ . Thus, V satisfies div ( μ / μ L ) ( V + H ) = 0 with lim ζ ± V = 0 , from (4), and ( V + H ) , ( μ / μ L ) ( V + H ) · n are continuous at the interfaces between two elastic media and vanish at the boundaries in contact with air.
We now define the set V of admissible fields V ˜ being continuous with V ˜ 0 for | ζ | + . Next, we introduce the energy defined over V by:
V ˜ V , E ( V ˜ ) = Y μ μ L 1 2 | V ˜ | 2 + V ˜ · H d χ α ( V ˜ ) , with α ( V ˜ ) = 1 φ P Y P V ˜ ( χ , + ) d χ Y V ˜ ( χ , ) d χ .
One can show using standard arguments of calculus of variations that the minimizer V of E defined by:
V = arg min V ˜ V E ( V ˜ ) ,
corresponds to V = V 2 H with α ( V ) = α 2 since, by the definition of V and from (A32), lim ζ V = 0 and lim ζ + V = α 1 . To conclude, we need the expression of E ( V ) . We multiply the relation div ( μ / μ L ) ( V + H ) = 0 by V, then by H, and integrate by parts. We get that Y μ μ L | V | 2 d χ + Y μ μ L V · H d χ = α 1 and, after straightforward calculations, that Y d χ μ μ L V · H = φ b h b 1 μ b μ L . It follows that, from (A44), we have:
E ( V ) = 1 2 φ b h b 1 μ b μ L α 1 .
Thus, by bounding the energy E ( V ) , we shall get a lower bound for α 1 , hence for b = α 1 + h b . To do so, we chose a test field V ˜ being piecewise linear along ζ , with V ˜ ( χ ) = f ( ζ ) and f ( ζ ) = 0 , b ( ζ + h b / ) / h b , b for ζ ( , h b ) , ( h b , 0 ) , and ( 0 , + ) , respectively. At this stage, b is a free parameter that we shall fix to minimize E ( V ˜ ) . It is easy to see that:
E ( V ˜ ) = h b 1 2 φ b μ b μ L + 1 φ b b h b 2 + φ b μ b μ L 1 b h b ,
whose minimum with respect to b is obtained for b * = ( h b / ) φ b ( μ b / μ L 1 ) / ( φ b μ b / μ L + 1 φ b ) and provides:
min E ( V ˜ ) = h b 2 φ b 2 ( μ b / μ L 1 ) 2 φ b μ b / μ L + 1 φ b
It is now sufficient to use (A45) along with (A46) and (A48) to find that:
b h b φ b μ b μ L + 1 φ b > 0 .

Appendix B.3. Lb0

To obtain a lower bound on L b , we rely on the variational formulation of the elementary problem on V 2 . First, we define the set T of admissible fields τ ˜ such that:
τ ˜ · n , continuous   at   the   interfaces   between   two   elastic   media , τ ˜ · n = 0 , a t χ = ± φ P / 2 , in Y P , τ ˜ , one   periodic   with   respect   to   χ in Y b , lim ζ τ ˜ = e x , lim ζ + τ ˜ = 0 .
Next, we introduce the complementary energy defined over T :
τ ˜ T , E * ( τ ˜ ) = Y b μ μ L τ ˜ μ μ L e x 2 d χ + Y P μ L μ | τ ˜ | 2 2 τ ˜ · e x d χ .
One can show from standard arguments of the calculus of variations that the minimizer τ of E * defined by:
τ = arg min τ ˜ T E * ( τ ˜ ) ,
corresponds to the stress field associated with the solution V 2 , that is to say:
τ = μ μ L ( χ V 2 + e x ) .
Now, multiplying by V 2 the equilibrium equation (4) associated with V 2 , integrating by parts, and using (A51) and (A53), we find the relation:
L b h b φ b μ b μ L + ( 1 φ b ) = Y P μ μ L ( V 2 + χ ) χ d χ + Y b μ μ L V 2 χ d χ = E * ( τ ) .
Thus, by bounding the complementary energy of the solution E * ( τ ) , we shall get a bound on the left hand side in (A54), hence on L b given by (5). To do so, we chose the piecewise constant test field τ ˜ such that τ ( χ ) = g ( ζ ) e x with g ( ζ ) = 1 , a, and 0 for ζ ( , h b / ) , ( h b / , 0 ) , and ( 0 , + ) , respectively. At this stage, a is a free parameter that we shall fix to minimize E * ( τ ˜ ) . Calculating the energy of such a test field gives:
E * ( τ ˜ ) = a μ b μ L 2 μ L μ b φ b h b + ( a 1 ) 2 h b ( 1 φ b ) .
Now, minimizing the energy with respect to a gives the optimum a * = φ b μ L μ b + 1 φ b 1 . By injecting a * in (A55) and using (A54), we finally obtain:
L b h b φ b μ L μ b + 1 φ b > 0 .

References

  1. Housner, G.W. Interaction of building and ground during an earthquake. Bull. Seismol. Soc. Am. 1957, 47, 179–186. [Google Scholar]
  2. Bielak, J. Dynamic behaviour of structures with embedded foundations. Earthq. Eng. Struct. Dyn. 1975, 3, 259–274. [Google Scholar] [CrossRef]
  3. Clouteau, D.; Aubry, D. Modification of the ground motion in dense urban areas. J. Comput. Acoust. 2001, 9, 1659–1675. [Google Scholar] [CrossRef]
  4. Kham, M.; Semblat, J.-F.; Bard, P.Y.; Dangla, P. Seismic site–city interaction: Main governing phenomena through simplified numerical models. Bull. Seismol. Soc. Am. 2006, 96, 1934–1951. [Google Scholar] [CrossRef]
  5. Guéguen, P.; Colombi, A. 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. 2016, 106, 2855–2864. [Google Scholar] [CrossRef]
  6. Guéguen, P.; Bard, P.Y.; Chavez-Garcia, F.J. Site-city seismic interaction in Mexico city–like environments: An analytical study. Bull. Seismol. Soc. Am. 2002, 92, 794–811. [Google Scholar] [CrossRef]
  7. Boutin, C.; Roussillon, P. Assessment of the urbanization effect on seismic response. Bull. Seismol. Soc. Am. 2004, 94, 251–268. [Google Scholar] [CrossRef]
  8. Ghergu, M.; Ionescu, I.R. Structure-soil-structure coupling in seismic excitation and city effect. Int. J. Eng. Sci. 2009, 47, 342–354. [Google Scholar] [CrossRef]
  9. Schwan, L.; Boutin, C.; Padrón, L.A.; Dietz, M.S.; Bard, P.Y.; Taylor, C. Site-city interaction: Theoretical, numerical and experimental crossed-analysis. Geophys. J. Int. 2016, 205, 1006–1031. [Google Scholar] [CrossRef] [Green Version]
  10. Brûlé, S.; Enoch, S.; Guenneau, S. Structured soils under dynamic loading: The metamaterials in Geotechnics. Rev. Fr. Geotech. 2017, 151, 4. [Google Scholar] [CrossRef]
  11. Ungureanu, B.; Guenneau, S.; Achaoui, Y.; Diatta, A.; Farhat, M.; Hutridurga, H.; Craster, R.V.; Enoch, S.; Brûlé, S. The influence of building interactions on seismic and elastic body waves. EPJ Appl. Metamat. 2019, 6, 18. [Google Scholar] [CrossRef]
  12. Cacciola, P.; Tombari, A. Vibrating barrier: A novel device for the passive control of structures under ground motion. Proc. R. Soc. A 2015, 471, 20150075. [Google Scholar] [CrossRef]
  13. Dertimanis, V.K.; Antoniadis, I.A.; Chatzi, E.N. Feasibility Analysis on the Attenuation of Strong Ground Motions Using Finite Periodic Lattices of Mass-in-Mass Barriers. J. Eng. Mech. 2016, 142, 04016060. [Google Scholar] [CrossRef]
  14. Palermo, A.; Krodel, S.; Marzani, A.; Daraio, C. Engineered metabarrier as shield from seismic surface waves. Sci. Rep. 2016, 6, 39356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Palermo, A.; Marzani, A. Control of Love waves by resonant metasurfaces. Sci. Rep. 2018, 8, 7234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Colombi, A.; Colquitt, D.J.; Roux, P.; Guenneau, S.; Craster, R.V. A seismic metamaterial: The resonant metawedge. Sci. Rep. 2016, 6, 27717. [Google Scholar] [CrossRef] [Green Version]
  17. Colquitt, D.J.; Colombi, A.; Craster, R.V.; Roux, P.; Guenneau, S. Seismic metasurfaces: Sub-wavelength resonators and Rayleigh wave interaction. J. Mech. Phys. Solids 2017, 99, 379–393. [Google Scholar] [CrossRef]
  18. Maurel, A.; Marigo, J.-J.; Pham, K.; Guenneau, S. Conversion of Love waves in a forest of trees. Phys. Rev. B 2018, 98, 134311. [Google Scholar] [CrossRef] [Green Version]
  19. Kerders, L.; Allard, J.F.; Lauriks, W. Ultrasonic surface waves above rectangular-groove gratings. J. Acoust. Soc. Am. 1998, 103, 2730–2733. [Google Scholar] [CrossRef]
  20. Hurd, R.A. The propagation of an electromagnetic wave along an infinite corrugated surface. Can. J. Phys. 1954, 32, 727–734. [Google Scholar] [CrossRef]
  21. Pendry, J.B.; Martin-Moreno, L.; Garcia-Vidal, F. Mimicking surface plasmons with structured surfaces. Science 2004, 305, 847–848. [Google Scholar] [CrossRef]
  22. McPhedran, R.C.; Botten, L.C.; Bliek, P.; Deleuil, R.; Maystre, D. Inductive grids in the region of diffraction anomalies: Theory, experiment and applications. IEEE Trans. Microw. Theory Tech. 1980, 28, 1119–1125. [Google Scholar]
  23. Ebbesen, T.W.; Lezec, H.J.; Ghaemi, H.F.; Thio, T.; Wolff, P. Extraordinary optical transmission through sub-wavelength hole arrays. Nature 1998, 391, 667–669. [Google Scholar] [CrossRef]
  24. Barnes, W.L.; Dereux, A.; Ebbesen, T.W. Surface plasmon subwavelength optics. Nature 2003, 424, 824–830. [Google Scholar] [CrossRef] [PubMed]
  25. Brûlé, S.; Enoch, S.; Guenneau, S. Role of nanophotonics in the birth of seismic megastructures. Nanophotonics 2019, 8, 1591–1600. [Google Scholar] [CrossRef] [Green Version]
  26. Lott, M.; Roux, P.; Garambois, S.; Guéguen, P.; Colombi, A. Evidence of metamaterial physics at the geophysics scale: The METAFORET experiment. Geophys. J. Int. 2019, 220, 1330–1339. [Google Scholar] [CrossRef]
  27. Mercier, J.F.; Cordero, M.L.; Félix, S.; Ourir, A.; Maurel, A. Classical homogenization to analyse the dispersion relations of spoof plasmons with geometrical and compositional effects. Proc. R. Soc. A 2015, 471, 20150472. [Google Scholar] [CrossRef]
  28. Norris, A.N.; Su, X. Enhanced acoustic transmission through a slanted grating. C. R. Mécanique 2015, 343, 622–634. [Google Scholar] [CrossRef]
  29. Achenbach, J.D. Wave Propagation in Elastic Solids; North-Holland: Amsterdam, The Netherlands, 1973. [Google Scholar]
  30. Maurel, A.; Marigo, J.-J.; Mercier, J.-F.; Pham, K. Modelling resonant arrays of the Helmholtz type in the time domain. Proc. R. Soc. A 2018, 474, 20170894. [Google Scholar] [CrossRef] [Green Version]
  31. Marigo, J.-J.; Maurel, A.; Pham, K.; Sbitti, A. Effective dynamic properties of a row of elastic inclusions: The case of scalar shear waves. J. Elast. 2017, 128, 265–289. [Google Scholar] [CrossRef] [Green Version]
  32. Maurel, A.; Marigo, J.-J.; Pham, K. Effective boundary condition for the reflection of shear waves at the periodic rough boundary of an elastic body. Vietnam J. Mech. 2018, 40, 303–323. [Google Scholar] [CrossRef]
  33. Bodet, L.; Dhemaied, A.; Martin, R.; Mourgues, R.; Rejiba, F.; Tournat, V. Small-scale physical modeling of seismic-wave propagation using unconsolidated granular media. Geophysics 2014, 79, T323. [Google Scholar] [CrossRef]
  34. Palermo, A.; Krödel, S.; Matlack, K.H.; Zaccherini, R.; Dertimanis, V.K.; Chatzi, E.N.; Marzani, A.; Daraio, C. Hybridization of guided surface acoustic modes in unconsolidated granular media by a resonant metasurface. Phys. Rev. Appl. 2018, 9, 054026. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Periodic array of plates decorated at their endings with spacing = 1 , height h P , and thickness φ P ; the substrate occupying a half-space is surmounted by a guiding layer of thickness h L able to support Love waves. The insets show a zoom on the two endings with heterogeneity surfaces S t = φ t h t and S b = φ b h b .
Figure 1. Periodic array of plates decorated at their endings with spacing = 1 , height h P , and thickness φ P ; the substrate occupying a half-space is surmounted by a guiding layer of thickness h L able to support Love waves. The insets show a zoom on the two endings with heterogeneity surfaces S t = φ t h t and S b = φ b h b .
Materials 13 01632 g001
Figure 2. Configuration of the effective problem (2) and (3): The region of the plates has been replaced by a homogeneous medium; the effective boundary condition and transmission conditions encapsulate the effects of the heterogeneities at the decorative endings of the plates.
Figure 2. Configuration of the effective problem (2) and (3): The region of the plates has been replaced by a homogeneous medium; the effective boundary condition and transmission conditions encapsulate the effects of the heterogeneities at the decorative endings of the plates.
Materials 13 01632 g002
Figure 3. Domain D where the energy is conserved in the absence of incoming/outcoming fluxes through Σ . The effective boundary condition on γ and jump conditions between Γ ± in (3) result in additional effective energies E t , b in (10).
Figure 3. Domain D where the energy is conserved in the absence of incoming/outcoming fluxes through Σ . The effective boundary condition on γ and jump conditions between Γ ± in (3) result in additional effective energies E t , b in (10).
Materials 13 01632 g003
Figure 4. Configuration of the array. The total thickness H P = h P + h b = 12 of the array and the total thickness H L = h L + h t = 8 of the layer are kept constant; = 1 and φ b = 1 , φ t = φ P = 0.5 . When the heterogeneities are considered, h b = h t = 1 .
Figure 4. Configuration of the array. The total thickness H P = h P + h b = 12 of the array and the total thickness H L = h L + h t = 8 of the layer are kept constant; = 1 and φ b = 1 , φ t = φ P = 0.5 . When the heterogeneities are considered, h b = h t = 1 .
Materials 13 01632 g004
Figure 5. Reference solutions: Dispersion relations of Love waves and modified Love waves (with a bottom layer of thickness h b ). In the presence of a thin layer b atop the guiding layer, the dispersion relation is modified ( Θ b in (16)) resulting in different shapes of the Love dispersion branches. SPP, Spoof Plasmon Polariton.
Figure 5. Reference solutions: Dispersion relations of Love waves and modified Love waves (with a bottom layer of thickness h b ). In the presence of a thin layer b atop the guiding layer, the dispersion relation is modified ( Θ b in (16)) resulting in different shapes of the Love dispersion branches. SPP, Spoof Plasmon Polariton.
Materials 13 01632 g005
Figure 6. Asymptotes of the SPPs (at ω n SPP / ( 2 π ) = 3 , 15 , 25 ) and modified SPPs (at ω n SPP , t / ( 2 π ) = 3.1 , 12.0 , 22.4 ).
Figure 6. Asymptotes of the SPPs (at ω n SPP / ( 2 π ) = 3 , 15 , 25 ) and modified SPPs (at ω n SPP , t / ( 2 π ) = 3.1 , 12.0 , 22.4 ).
Materials 13 01632 g006
Figure 8. Magnified views of the two avoided crossings AC1 and AC2 from Figure 7 (Case 4).
Figure 8. Magnified views of the two avoided crossings AC1 and AC2 from Figure 7 (Case 4).
Materials 13 01632 g008
Figure 9. Displacement fields corresponding to the two branches of guided waves ( β 1 and β 2 ) at ω / ( 2 π ) = 24 for the undecorated plates (top) and decorated plates (bottom). On each panel, the field from direct numerics is plotted for x < 0 , and the homogenized solution from (19) is plotted for x > 0 . In both cases, the displacement at z = H L is unitary, which allows for a quantitative comparison without any tuning parameter.
Figure 9. Displacement fields corresponding to the two branches of guided waves ( β 1 and β 2 ) at ω / ( 2 π ) = 24 for the undecorated plates (top) and decorated plates (bottom). On each panel, the field from direct numerics is plotted for x < 0 , and the homogenized solution from (19) is plotted for x > 0 . In both cases, the displacement at z = H L is unitary, which allows for a quantitative comparison without any tuning parameter.
Materials 13 01632 g009
Figure 10. Left: Real part of the reflexion coefficient R ( 1 , 1 ) in colorscale against β and ω c S β . The dashed white line corresponds to an incident propagating wave at oblique incidence with β = ω c S sin 45 . Right: Repartition of the energies in the bottom and top heterogeneities (upper panel) and in the substrate, layer, and plate (normalized with the total energy); see the lower panel. Open symbols are obtained from direct numerics, Equations (24), and plain lines from the homogenized problem, Equations (26) and (27).
Figure 10. Left: Real part of the reflexion coefficient R ( 1 , 1 ) in colorscale against β and ω c S β . The dashed white line corresponds to an incident propagating wave at oblique incidence with β = ω c S sin 45 . Right: Repartition of the energies in the bottom and top heterogeneities (upper panel) and in the substrate, layer, and plate (normalized with the total energy); see the lower panel. Open symbols are obtained from direct numerics, Equations (24), and plain lines from the homogenized problem, Equations (26) and (27).
Materials 13 01632 g010
Figure 11. Surface densities of energy in the actual problem (computed numerically) and given by the homogenized solution at ω / ( 2 π ) = 19 (corresponding to the resonance of the Love type with R = 1 , the energy is stored in the layer and in the bottom layer), at ω / ( 2 π ) = 22.7 (resonance with R = 1 of the SPP type; the energy is stored in the plate), and at ω / ( 2 π ) , which is a standard case ( R 1 , the energy is spread).
Figure 11. Surface densities of energy in the actual problem (computed numerically) and given by the homogenized solution at ω / ( 2 π ) = 19 (corresponding to the resonance of the Love type with R = 1 , the energy is stored in the layer and in the bottom layer), at ω / ( 2 π ) = 22.7 (resonance with R = 1 of the SPP type; the energy is stored in the plate), and at ω / ( 2 π ) , which is a standard case ( R 1 , the energy is spread).
Materials 13 01632 g011
Table 1. Elastic material properties (in arbitrary units).
Table 1. Elastic material properties (in arbitrary units).
SubstrateLayerPlateBottomTop
μ S = 2000 μ L = 72 μ P = 14.4 μ b = 1600 μ t = 1600
ρ S = 2000 ρ L = 1800 ρ P = 250 ρ b = 2500 ρ t = 2500
c S = 1000 c L = 200 c P = 240 c b = 800 c t = 800
Table 2. Homogenized coefficients entering in the effective conditions (3) in (5).
Table 2. Homogenized coefficients entering in the effective conditions (3) in (5).
Coeff. in (5) L t ρ ^ b b L b l b
h t = 0 02500 h b = 0 0.13240.01200
h t = 1 102500 h b = 1 0.051122.23570
Table 3. Geometries of the four configurations, whose dispersion relations are given in Figure 7.
Table 3. Geometries of the four configurations, whose dispersion relations are given in Figure 7.
h L h b H L h P h t H P
Love808000
Modified Love718000
Case 180812012
Case 280811112
Case 371812012
Case 471811112

Share and Cite

MDPI and ACS Style

Pham, K.; Maurel, A.; Félix, S.; Guenneau, S. Hybridized Love Waves in a Guiding Layer Supporting an Array of Plates with Decorative Endings. Materials 2020, 13, 1632. https://doi.org/10.3390/ma13071632

AMA Style

Pham K, Maurel A, Félix S, Guenneau S. Hybridized Love Waves in a Guiding Layer Supporting an Array of Plates with Decorative Endings. Materials. 2020; 13(7):1632. https://doi.org/10.3390/ma13071632

Chicago/Turabian Style

Pham, Kim, Agnès Maurel, Simon Félix, and Sébastien Guenneau. 2020. "Hybridized Love Waves in a Guiding Layer Supporting an Array of Plates with Decorative Endings" Materials 13, no. 7: 1632. https://doi.org/10.3390/ma13071632

APA Style

Pham, K., Maurel, A., Félix, S., & Guenneau, S. (2020). Hybridized Love Waves in a Guiding Layer Supporting an Array of Plates with Decorative Endings. Materials, 13(7), 1632. https://doi.org/10.3390/ma13071632

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop