1. Introduction

Horizontally polarized shear surface waves (SH-waves) propagating in multilayered plates resemble Love waves ^[1] in polarization, but differ in absence of a contacting half-space (substrate), and, hence excluding necessity to impose Sommerfeld’s emission condition:

(1.1)

where is the displacement field in the substrate; is the coordinate along depth of the substrate, and is the unit normal to the plane boundary of the substrate; see Figure 1.

As will be shown later, absence of condition (1.1) results in a completely different behavior of the SH-waves in layered plates comparing to Love waves. For example, it is known an existence inequality ^[2] for a genuine Love wave propagating in an isotropic traction-free layer contacting with isotropic substrate:

(1.2)

where denotes speed of the corresponding shear bulk wave propagating in , and polarized in direction. Violating inequality (1.2) prevents Love waves to exist. As will be shown later, the SH-waves in two-layered plates exist at any physically admissible properties of homogeneous layers and at traction-free, clamped, or mixed boundary conditions (one outer surface is traction free, and another is clamped).

The main method used for constructing analytic and numerical solutions for SH and Love waves, is based on a combination of a complex formalism ^{[3, 4]} and the modified transfer matrix (MTM) method ^{[5, 6]}. The latter method along with the multiprecision computations allowed us to avoid numerical instability of the solution arising at layers of large thickness and high frequencies of the waves. This is known as the “large -problem’, where is the layer thickness and is frequency; see ^{[7, 8, 9, 10]}. We should also mention that herein we do not consider the exceptional waves corresponding to appearing Jordan blocks in the Christoffel equation; see ^[5].

Other methods for analyzing propagation of SH and Love waves in connection to geophysical and geotechnical applications, are presented in recent publications ^{[11, 12, 13, 14]}.

2. Basic Notations

All the regarded layers and a substrate are assumed to be homogeneous, anisotropic and linearly hyperelastic. Equations of motion for homogeneous anisotropic elastic medium can be written in the form:

(2.1)

where is the material density, and is the elasticity tensor assumed to be positive definite:

(2.2)

Remark 2.1. a) The other assumption concerns symmetry of the elasticity tensor. It will be assumed that all the regarded materials possess planes of elastic symmetry coinciding with the sagittal plane , where vector is the polarization vector of the SH-wave. This is achieved by the elasticity tensor belonging to the monoclinic system, and the latter is equivalent to vanishing all of the decomposable components of the tensor having odd number of entries of the vector in the orthogonal basis in generated by the vector and any two orthogonal vectors belonging to the sagittal plane.

b) It will be shown later that assuming monoclinic symmetry provides a sufficient condition for the surface tractions acting on any plane to be collinear with vector .

Following ^[5], we will seek a horizontally polarized shear wave in a layer or substrate in the form:

(2.3)

where coordinate is as defined in (1.1); is the unknown scalar complex-valued function; the exponential multiplier in (2.3) corresponds to propagation of the plane wave front along direction with the phase speed ; is the wave number.

Remark 2.2. The displacement field defined by (2.3) is generally complex. In reality, either real or imaginary part of the right-hand side of (2.3) represents physical displacement field that will be implicitly assumed in the subsequent analysis. However, retaining complex expressions for the displacement field, will allow us to describe situations with the phase shift in a more convenient manner.

Substituting representation (2.3) into Eq. (2.1) and taking into account Remark 2.1.a, yields the following differential equation:

(2.4)

Characteristic equation for the differential equation (2.4), known also as the Christoffel equation, has the form:

(2.5)

Left-hand side of Eq. (2.5) represents a polynomial of degree 2 with respect to the Christoffel parameter . Thus, for the monoclinic elastic symmetry only two partial waves form the regarded SH-wave in a layer, and one partial wave with forms Love wave in the substrate; see ^[5].

The following lemma flows out from solving the Cauchy problem for Eq. (2.4):

Lemma 2.1. A necessary and sufficient condition for the real-analytic solution of Eq.(2.4), to be a non-zero function, is a simultaneous non-vanishing and its first derivative at some .

Remark 2.3. a) For an orthotropic medium and the SH-wave propagating in a direction of the principle elasticity, Eq.(2.5) is simplified:

(2.6)

The solution for Eq.(2.6) is:

(2.7)

For the considered case, the general solution of Eq. (2.4) can be represented in the form:

(2.8)

where is the positive root in (2.7) . For a Love wave in a substrate with attenuation condition (1.1) the considered solution becomes

(2.9)

in representation (2.9) we chose the positive root for .

b) Supposing roots of the Christoffel equation (2.5) are multiple, we arrive at necessity to modify the solution of Eq. (2.4) by placing a logarithmic term (this corresponds to appearing one Jordan block in Eq. (2.4) that is reduced to a system of the first order). However, in the following analysis the case of multiple roots will not be studied; see ^{[5, 6]} for the solutions related to multiple roots.

c) Representation (2.3) allows us to express surface tractions acting on any plane in the form:

(2.10)

Assumption of monoclinic symmetry ensures surface tractions (2.10) to be collinear with vector . For an orthotropic material with axes of elastic symmetry coinciding with vectors and expression (2.10) reduces to

(2.11)

As before, vector is necessary collinear with vector . Expression (2.11) shows that on a boundary plane of a traction-free orthotropic plate the following condition must be satisfied: .

d) There is a major difference between bulk and surface SH-waves. While both have the same polarization defined by vector , the corresponding stress components are different; see Figure 2. This concerns all kinds of surface SH-waves including the genuine Love waves satisfying attenuation condition (1.1) and the so-called leakage Love waves propagating with speed . The latter wave has the same speed and polarization as the corresponding bulk wave, does not attenuate in a substrate, but it is not a bulk wave in a substrate due to different stress components it generates. Indeed, for a surface SH-wave propagating in an orthotropic medium the stress vector is defined by Eq. (2.11), while for a bulk SH-wave, the stress vector is:

(2.12)

3. Energy of SH and Love Waves

3.1. Specific Kinetic and Elastic (potential) Energy

Herein, we derive expressions for specific kinetic and elastic (potential) energy of the SH-waves. Taking into account representation (2.3), the specific kinetic energy can be defined by:

(3.1)

where the following relation between the phase speed and frequency is used:

(3.2)

Equations (2.4), (3.1), and (3.2) allow us to represent the specific kinetic energy in the form:

(3.3)

Another useful expression flows out from (3.1) and (3.2):

(3.4)

Similarly, the specific elastic energy can be defined by:

(3.5)

Remark 3.1. a) In view of Remark 2.2, expressions (3.1) and (3.5) coincide with the corresponding expressions for kinetic and elastic specific energy, obtained without using complex displacement fields.

b) Analysis of expressions (2.1), (2.3), (2.6), and (2.8) reveals that for the regarded waves , due to presence of the generally non-constant function . At the same time, for bulk waves , and hence from (2.1) and (2.3) we arrive at ; see also ^{[2, 3]} for discussions.

Proposition 3.1. a) If at some finite value of the phase speed the corresponding frequency vanishes, then both specific kinetic and elastic energies vanish also.

b) The specific kinetic energy vanishes at any , on a plane , if function vanishes at this .

c) The specific elastic energy does not vanish at any finite value of the phase speed and any non-vanishing frequency .

Proofs of conditions a) and b) are obvious. Proof of condition c) follows from the positive definite condition for the elasticity tensor, Lemma 2.1, and expressions (3.2), (3.5).

3.2. Group Speed

The vector-valued group speed is defined by ^{[2, 3]}:

(3.6)

where denotes gradient with respect to the independent spatial variable . For the subsequent analysis the scalar group speed will also be needed:

(3.7)

Now, combining (3.3), (3.4), and (3.7) yields:

(3.8)

where as before, stands for the phase speed.

Proposition 3.2. a) At any physically admissible properties of a medium and any SH-wave propagating with the finite phase speed , the corresponding group speed is delimited from zero.

b) If at , where takes some finite value, then .

Proof a) flows out from observation that the radicand in (3.8) is strictly positive due to (2.2) and Lemma 2.1. Proof b) is obvious.

3.3. Ray Speed

The vector-valued ray speed can be defined by (see ^[3]):

(3.9)

where is the flux of elastic energy:

(3.10)

The corresponding scalar ray speed is:

(3.11)

Substituting (2.3) into (3.11) and exploiting (3.3), (3.5), yields:

(3.12)

Proposition 3.3. a) At any physically admissible properties of a medium and any SH-wave propagating with the finite phase speed , the corresponding ray speed is delimited from zero.

b) If at , where is finite, then .

c) A necessary and sufficient condition for , is as follows:

(3.13)

Proofs a) and b) are analogous to the proof of Proposition 3.2. Proof c) follows directly from (3.8), (3.12), with account of (3.1), (3.5).

4. SH-waves in a Single-layered Orthotropic Plate

Hence it will be assumed that vectors , and coincide with the axes of elastic symmetry of an orthotropic medium.

Remark 4.1. It can be shown (see ^[6]) that regardless of boundary conditions and at imaginary roots of Eq. (2.6), no SH-wave can propagate in directions of elastic symmetry of an orthotropic single-layered plate. Thus, the following inequality

(4.1)

naturally arising from (2.7), delivers a necessary condition for existing surface SH-wave. Thus, for the regarded plate all surface SH-waves are necessary supersonic, since the radicand in the right-hand side of (4.1) defines speed of the corresponding shear bulk wave . In this section we assume condition (4.1) to hold.

4.1. Traction-free Plate

Herein we consider a single-layered plate with the traction-free boundary conditions:

(4.2)

where is the thickness of the plate (we choose origin of coordinates at the median plane).

For such a plate, finding function from (2.6), (4.2), yields:

(4.3)

where is defined by (2.7).

Proposition 4.1. a) On planes , where

(4.4)

both the displacement field and specific kinetic energy vanish. That is equivalent to existence of the internal immovable layers under propagating SH-wave on a traction-free plate.

b) At any finite phase speed satisfying inequality (4.1), there are no waves propagating at vanishing frequency (both phase speed and frequency are delimited from zero).

Proof a) follows from considering zeroes of the function, defined by (4.3). Proof b) follows from analyzing expressions (4.3), (3.2). It reveals that no non-trivial solutions exist at .

4.2. Clamped Plate

For a single-layered plate with clamped outer surfaces, boundary conditions are:

(4.5)

Finding function from Eq. (2.4) and satisfying boundary conditions (4.5), yields:

(4.6)

Similarly to Proposition 4.1, we have

Proposition 4.2. a) On planes , where

(4.7)

both the displacement field and specific kinetic energy vanish. That is equivalent to existence of the internal immovable layers under propagating surface SH-wave on a clamped plate.

b) At any finite phase speed satisfying inequality (4.1), there are no waves propagating at vanishing frequency (both phase speed and frequency are delimited from zero).

4.3. Plate with Mixed Boundary Conditions

Herein we consider a plate with traction-free upper and clamped bottom surface:

(4.8)

Direct analysis reveals that function satisfying homogeneous boundary conditions (4.8) takes the form:

(4.9)

Proposition 4.3. a) On planes , where

(4.10)

both the displacement field and specific kinetic energy vanish. That is equivalent to existence of immovable layers under propagating surface SH-wave on a clamped plate.

b) At any finite phase speed satisfying inequality (4.1), there are no waves propagating at vanishing frequency (both phase speed and frequency are delimited from zero).

4.4. Concluding Remarks

For monoclinic single-layered plates and all the considered boundary conditions (i) the admissible speed interval is transonic:

(4.11)

(ii) at any phase speed satisfying (4.11) there are immovable longitudinal layers, and (iii) there are no anomalous SH-waves corresponding to vanishing frequency.

Below, the dispersion relations will be plotted in terms of the phase speed and frequency , or the phase speed and wave number , the latter is connected with and by (3.2). The typical dispersion curves for a homogeneous traction-free plate are presented in Figure 3.

These curves relate to a plate with . The lower mode dispersion curves for this plate but with different boundary conditions are plotted in Figure 4.

It should also be noted that a lower mode wave is usually the easiest to launch, and it appears to be the most informative due to its dispersion. The presented data reveal also that for a given plate its boundary conditions can be distinguished easily by analyzing lower mode dispersion curves.

5. Love Waves

In this section we assume inequality (1.2) to hold ensuring existence of Love wave. Different reference coordinate systems for the layer (with the origin at the median plane) and for the substrate (with the origin at the interface) will be chosen.

5.1. Analytical Solutions

For Love wave we consider the following boundary conditions: (i) the traction-free boundary conditions at the outer surface of a layer, (ii) contact type conditions at the interface, and (iii) the attenuation condition (1.1) at :

(5.1)

In Eqs. (5.1) and further index 1 is referred to the layer, and 2 to the substrate.

Applying the Modified Transfer Matrix (MTM) method ^[5], functions, which define the corresponding displacement fields, can be represented in the form:

(5.2)

The implicit secular equation for the wave number can be written in the form ^[5]:

(5.3)

Proposition 5.1. a) Suppose that

(5.4)

then on planes (in a layer), where

(5.5)

(5.6)

and satisfies Eq. (5.3), both the displacement field and specific kinetic energy vanish.

b) No Love waves propagate, when

(5.7)

c) At the phase speed , there is a lower mode leakage Love wave propagating with vanishing wave number.

Proof a) flows out from analyzing expression (5.2) for function . Values for defined by (5.5) and (5.6), are zeroes of this function.

Proof b) follows from an observation that at the phase speed satisfying the first inequality in (5.7), both and are imaginary due to (2.7) and (1.2). But, Eq. (5.3) at imaginary reduces to

(5.8)

The latter equation has no positive solutions with respect to the wave number . Now, at the phase speed satisfying the second inequality in (5.7), we see that the exponent term in (5.2)₂ does not satisfy the attenuating condition at .

To prove c) we need to consider Eq. (5.3) at small :

(5.9)

Supposing that , which implys , and letting , which gives a wave corresponding to the lower branch, we arrive at the following condition:

(5.10)

uniformly on any compact in . The latter proves condition c).

Remark 5.1. a) However, leakage Love wave introduced in Proposition 5.1.c, does not coincide with the horizonally polarized bulk wave in the substrate. That is because leakage Love wave travelling with the same speed as the shear bulk wave in the substrate (or better speaking, travelling with the phase speed ), leads to different stress components; see Figure 2. Moreover, as expresion (2.11) shows this leakage wave is a stress-free wave.

b) Another interesting observation concerns amplitude distribiution along depth of a layer and substrate. According to representations (2.3) and (5.2), the distribution in both layer and substrate tends to a constant as .

5.2. Numerical Example

Herein we consider an orthotropic layer lying on an orthotropic substrate. The outer plane (Figure 1) of the layer is traction-free. The layer and substrate are characterized by the following parameters:

(5.11)

The corresponding dispersion curves are presented in Figure 5.

These data demonstrate presence of a lower mode leakage wave ^[5] propagating with the phase speed and having vanishing frequency.

6. Two-layered Orthotropic Plates

It is assumed that (i) both layers are orthotropic with axes of elastic symmetry coincident with vectors , and ; and (ii) the corresponding shear bulk waves differ:

(6.1)

Remark 6.1. a) If inequality (6.1) violates, then the two-layered plate appears as a single-layered, with respect to the considered SH-wave.

b) In this section and further we assume that each layer has its own reference coordinate system with the origin lying in the median plane of a layer.

6.1. Traction-free Plate

Boundary conditions for a traction-free plate are:

(6.2)

where lower indicies are referred to the corresponding layers.

Applying the MTM method ^{[5, 6]}, functions which define the displacement field in the corresponding layeres, can be represented in the form:

(6.3)

at the wave number satisfying the following equation ^[6]:

(6.4)

Proposition 6.1. a) Suppose that

(6.5)

where is the bulk wave speed in the corresponding layer, then on planes where

(6.6)

and

(6.7)

and satisfies Eq. (6.4)), both the displacement field and specific kinetic energy vanish in a layer with the minimal bulk wave speed .

b) Suppose that

(6.8)

(the phase speed is transonic in both layers), then on planes where satisfies Eq. (6.6) and satisfies Eq. (6.7), the displacement field and specific kinetic energy vanish in both layers.

c) At the phase speed , where

(6.9)

there is a lower mode SH-wave propagating with vanishing wave number .

Proofs a) and b) flow out from expression (6.3) for functions . Values for defined by (6.6) and (6.7), are zeroes of these functions.

To prove c) we need to consider Eq. (6.4) at small :

(6.10)

Equating to zero the coefficient at in the left-hand side of Eq. (6.10), we arrive at the solution for the phase speed given by (6.9).

Definition 6.1. A lower mode surface SH-wave propagating with the phase speed and with vanishing wave number will be called the anomalous surface wave.

Proposition 6.1 along with expressions (3.1), (3.2), (3.5) ensure:

Corollary 6.1. a) Both kinetic and potential energy vanish at the anomalous surface wave.

b) At sufficiently small and real in both planes there can be no planes with vanishing displacement field (all SH-waves in the vicinity of the anomalous SH-wave do not have planes at which the displacement field vanishes).

Remark 6.2. a) Proposition 6.1 shows that planes with vanishing kinetic energy arise only if the phase speed becomes transonic for the corresponding layer.

b) Direct analysis reveals that the wave speed satisfies the inequalities:

(6.11)

c) The considered SH-waves in the vicinty of the anomalous SH-wave resemble solitones, since their wave number tends to zero, as .

d) Taking into account Eqs. (3.1), (3.5), we have for the lower branch of the SH-waves at :

(6.12)

In obtaining (6.12), restriction of at is accounted (herein denotes the corresponding Hörmander’s norm). Thus, conditions (6.12) ensure minimum energy required to initiate lower branch of the SH-wave at . Moreover, similarly to the leakage Love wave discussed in Remark 5.1, the lower branch of the SH-waves at is a stress-free wave, and the amplitude distribution in both layers tends to a constant at .

6.2. Clamped Plate

Boundary conditions for a clamped plate are:

(6.13)

Application of the Modified Transfer Matrix (MTM) method gives functions in the form:

(6.14)

The wave number in representation (6.14) satisfies the following equation ^[6]:

(6.15)

Similarly to the preceding case, we have

Proposition 6.2. a) Suppose that

(6.16)

then on planes where

(6.17)

and

(6.18)

and satisfies Eq. (6.4)), both the displacement field and specific kinetic energy vanish in a layer with the minimal bulk wave speed .

b) Suppose that

(6.19)

(the phase speed is transonic for both layers), then on planes where

(6.20)

the displacement field and specific kinetic energy vanish.

c) At the phase speed , where satisfies Eq. (6.9), there is a lower mode SH-wave propagating with vanishing wave number .

6.3. Plate with Mixed Boundary Conditions

Boundary conditions for the considered plate are:

(6.21)

Application of the Modified Transfer Matrix (MTM) method gives functions in the form:

(6.22)

For the considered case, the wave number satisfies the following equation ^[6]:

(6.23)

Similarly to the preceding cases, we have

Proposition 6.3. a) Suppose that

(6.24)

then on planes where satisfies Eqs. (6.17), (6.18), both the displacement field and specific kinetic energy vanish in a layer with the minimal bulk wave speed .

b) Suppose that

(6.25)

(the phase speed is transonic for both layers), then on planes where satisfies Eqs (6.20), the displacement field and specific kinetic energy vanish in both layers.

c) At the phase speed , where satisfies Eq. (6.9), there is a lower mode SH-wave propagating with vanishing wave number .

6.4. Numerical Examples

Thus, for a two-layered monoclinic plate with all the considered boundary conditions (i) the admissible speed interval is partly transonic

(6.26)

(ii) at any phase speed satisfying (6.26) there are immovable longitudinal planes, and (iii) there are anomalous lower-mode SH-waves at .

Typical dispersion curves (in terms of the frequency and phase speed ) for a two-layered traction-free plate with

(6.27)

are presented in Figure 6.

Existence of the anomalous SH-wave at () propagating in a two-layered plate, resembles a lower mode leakage Love wave propagating with vanishing frequency at the speed

(6.28)

see Remark 5.1 and Figure 5.

We also analyzed influence of the depth increase of a bottom layer in a two-layered plate on behavior of the lower branch of dispersion curves; see Figure 7.

The presented data reveal that while ratio increases from 1 to 25, the lower mode dispersion curves of the SH-wave are getting closer to the corresponding curve for Love wave; compare with Figure 5. However, even at the highest computed ratio there is a principle possibility to differentiate a SH-wave from Love wave, since there limiting speeds remain distinguishable.

7. Multilayered Orthotropic Plate

Herein, we consider multilayered traction-free plates consisting of alternating layers (up to 31) of the unit thickness: . It is assumed that (i) all the layers are orthotropic with the axes of elastic symmetry coincident with vectors , and ; and (ii) the corresponding shear bulk waves of the contacting layers differ. Mechanical properties of the layers needed for the subsequent analysis are as follows:

(7.1)

Since for a plate with more than 2 orthotropic layers no explicit analytical solution for the dispersion relations is known, a numerical algorithm exploiting the modified transfer matrix method will be adopted.

7.1. MTM Method

Applying the MTM method to a more general case of the -layered plate with orthotropic layers of arbityrary thickness, we arrive at the implicit secular equation in the form ^[6]:

(7.2)

In Eq. (7.2) vectors and represent the traction-free boundary conditions at outer planes, and are the transfer matrices:

(7.3)

where .

Equation (7.2) is the transcendental equation with respect to the phase speed and wave number . Applying relation (3.2), allows us to reformulate Eq. (7.2) in terms of the phase speed and frequency . Direct evaluation of the determinant of the transfer matrix reveals that . Thus, the determinants do not vanish at , and the transfer matrices in (7.2) are correctly defined at any (real and positive) phase speed satisfying the condition

(7.4)

Condition (7.4) is assumed to hold below.

Remark 7.1. Expressions (7.3) correspond to orthotropic layers. When layers with monoclinic symmetry need to be considered, a more general expression accounting possible appearing of multiple roots in the Christoffel equation should be used; see ^{[5, 6]}.

7.2. Numerical Results

Analysis of dispersion curves for SH-waves propagating in -layered plates with different number of alternating layers (7.1) revealed that all the considered plates have similar lower branches ending at some finite phase speeds , and the latter depend on number of layers; see Figure 8.

All these curves resemble one obtained for a two-layered plate; see Figure 6.

The following plot in Figure 9 demonstrates the first 50 dispersion curves in a 31-layered traction-free plate.

The most interesting observable features in Figure 9 are: (i) existence of the anomalous wave speed ; (ii) fusion of the dispersion curves at their origin; and (iii) bundled structure of these curves.

Remark 7.2. It should be noted that while properties (i) and (ii) are common for all the considered plates (7.1) at any , property (iii) is becoming stably observable at a relatively large number of layers ().

Another plot in Figure 10 represents variation of the lower branch of the dispersion curves in a 31-layered traction-free plate due to the thickness variation of the 15-th layer (that is the middle layer).

These data shows that even relatively small 10% variation of the thickness of the middle layer is detectable by analyzing behavior of the lower branch in the vicinity of . Similar results were obtained at variation of the shear bulk wave speed of the middle layer. Thus, analyzing lower branches of dispersion curves provides a robust and nondestructive method for evaluating physical and geometrical properties of the internal layers.

Remark 7.3. To achieve the desired computational precision and avoid possible numerical instability, the main computations were implemented with mantissas having over 500 decimal digits. That was achieved by applying Bailey’s numerical package for multiprecision computations.

Acknowledgements

Author thanks Dr. D.H. Bailey for the Multiprecision Package, the Russian Foundation for Basic Research (Grant 04-01-00781) and the Russian Academy of Sciences Program OEMMPU No.12 for partial financial support.

References

[1]	A.E.H. Love, Some Problems of Geodynamics. Cambridge University Press, London (1911).
	In article
[2]	E. Dieulesaint and D. Royer, Elastic Waves in Solids. Wiley, N.Y. (1980).
	In article
[3]	T.C.T. Ting and D.M. Barnett, Classifications of surface waves in anisotropic elastic materials. Wave Motion 26 (1997) 207-218.
	In article		CrossRef
[4]	S.V. Kuznetsov, Subsonic Lamb waves in anisotropic plates. Quart. Appl. Math. 60 (2002) 577-587.
	In article
[5]	S.V. Kuznetsov, Love waves in stratified monoclinic media. Quart. Appl. Math. 62 (2004) 749-766.
	In article
[6]	S.V. Kuznetsov, SH-waves in multilayered plates. Quart. Appl. Math. 64 (2006) 153-165.
	In article
[7]	S.V. Kuznetsov, Love waves in non-destructive diagnostics of layered composites, Acoustical Physics, 56 (2010) 877-892.
	In article		CrossRef
[8]	M.J.S. Lowe, Matrix techniques for modeling ultrasonic waves in multilayered media. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control 42 (1995) 525-542.
	In article		CrossRef
[9]	D. Lévesque and L. Piché, A robust transfer matrix formulation for the ultrasonic response of multilayered absorbing media. J. Acoust. Soc. 92 (1992) 452-467.
	In article		CrossRef
[10]	M. Castaings and B. Hosten, Transfer matrix of multilayered absorbing and anisotropic media. Measurments and simulations of ultrasonic wave propagation through composite materials. J. Acoust. Soc. Am. 94 (1993) 1488-1495.
	In article		CrossRef
[11]	P. Michaels and V. Gottumukkula, Theory of viscoelastic Love waves and their potential application to near-surface sensing to permeability. In: Advances in near-surface seismology and ground-penetrating radar. Geophysical Developments Series. (2010) 263-278.
	In article		CrossRef
[12]	D. Restrepo, J.D. Gomez, and J.D. Jaramillo, SH wave number Green’s function for a layered, elastic half-space. Part I: Theory and dynamic canyon response by the discrete wave number boundary element method. Pure Appl. Geophys. (2014) 1-14.
	In article
[13]	M. Behm and R. Snieder, Love waves from local traffic noise interferometry. The Leading Edge. 32 (2013) 628-632.
	In article		CrossRef
[14]	J. Xia, X. Yin, and Y. Xu, Feasibility of determining Q of near-surface materials from Love waves. J. Appl. Geophysics. 95 (2013) 47-52.
	In article		CrossRef