Elastodynamic response of an orthotropic layered elastic half-space to time-harmonic loading

Abstract

The elastodynamics of an orthotropic half-space coated by a thin orthotropic layer is theoretically investigated in this article. We newly propose explicit expressions of free Rayleigh waves in a layered half-space that are dependent on only one unknown constant representing amplitude. The main contribution is on deriving, in a simple manner, the theoretical predictions of far-field Rayleigh wave motion arising from time-harmonic loads using elastodynamic reciprocity theorems. These are the very first closed-form exact solutions found for the forced motion of Rayleigh waves in a layered half-space of orthotropic materials. To demonstrate the theoretical results, computation of Rayleigh wave motion in a jointed rock, including a layer of quartz-schist and a half-space of soil, is considered. We present the phase and group dispersion curves superimposed with the amplitude spectra that provide useful information on wave modes, frequencies, and displacement amplitudes. The inclusion of the amplitude spectra in the dispersion curves is a significant improvement over other dispersion curves currently available in the literature. The analytical predictions are compared with numerical results found by finite element analysis, and they show excellent agreement for the cases of a uniform distributed load and a varying distributed load both applied over a strip on the layer surface. The calculations obtained in the current study could generally be very useful for applications in seismology and materials characterization of coated structures.

Keywords

Rayleigh waves layered half-space orthotropic materials reciprocity theorem time-harmonic loads

Introduction

The study of elastic waves in layered half-space structures is of great importance in seismology and geophysics. Seismologists and geophysicists are interested in understanding and predicting the behavior of seismic waves, including P, S, and Rayleigh surface waves, under different earthquake loadings. A numerical analysis of stress wave propagation through a joint rock mass, including three orthogonal sets, was explored.¹ It used orthotropic continua to replace the discontinuous media and transferred the problem into an equivalent continuous model. The dynamic response agreed with the in-situ records and the response obtained by discrete modeling. Rayleigh waves in a self-reinforced layer over an inhomogeneous half-space were studied for geophysical and civil engineering applications.² In particular, they discussed the influence of inhomogeneity parameters and layer thickness on the Rayleigh wave phase velocity. In another study, Gupta and Ahmed³ obtained analytical solutions of Rayleigh waves in a layered structure consisting of an anisotropic layer and a sandy medium. Moreover, Rayleigh waves are commonly employed to inspect and evaluate mechanical and civil engineering structures such as fiber-reinforced polymer-retrofitted concrete structures⁴ and composite coating structures.^5,6

Fundamentals of propagation of Rayleigh surface waves in a layered half-space are widely available in the literature.^7–9 The topic is also discussed in detail in numerous research papers. In particular, Achenbach and Keshava¹⁰ found the dispersion curves of Rayleigh waves in a layered half-space and Tiersten¹¹ explored the effect of thin films on the propagation of guided waves. The approximate characteristic equations for Rayleigh waves in an orthotropic layered half-space with a sliding contact at the interface were reported by Vinh et al.¹² Elastodynamics of a layered half-space including a thin, soft layer with a clamped upper face overlying a half-space was recently considered.¹³ The authors derived the non-traditional boundary conditions along the substrate surface combined with the effect of the layer using a long-wave high-frequency procedure.

The problem of Rayleigh wave motions in layered half-spaces induced by forces is rather complicated and thus there are few published works. The conventional approach is based on integral transform techniques, commonly used to solve wave motions in single solids such as a thin plate and a half-space,^7,14 becomes extremely difficult to use in solving the dynamic responses of Rayleigh waves in an orthotropic layered half-space. Xu and Ma¹⁵ used the direct stiffness method and Fourier transform to investigate the responses of a multilayered half-space to a spatially periodic moving load. They derived the displacement and traction components along layer interfaces for both the internal soil layers and the soil half-space. However, the soil materials were all considered homogeneous and isotropic. Based on the dual vector form and the Hankel transform, elastodynamic responses of a multi-layered transversely isotropic piezoelectric medium under time-harmonic loadings were studied.¹⁶ To obtain the piezoelectric field, however, a numerical scheme is required to evaluate the complex integrals that involve products of polynomial and Bessel functions. Recently, Bratov et al.¹⁷ proposed an asymptotic hyperbolic-elliptic expression for Rayleigh waves in a layered half-space subjected to a vertical load.

In recent years, another method has been introduced based on the reciprocity theorem, i.e. a relation between displacements, traction components and forces for two different loading states of a body, to obtain the closed-form solutions of guided wave motions.¹⁸ The approach is simpler than the conventional approach using integral transform techniques¹⁹ and applicable to inhomogeneous solids,^20,21 multilayered structures^22,23 and anisotropic materials.²⁴ The validity of the reciprocity-based approach was verified for Rayleigh waves in a half-space; the results found by the reciprocity application are identical to the ones by the integral transform methods.^19,24,25

One of the main contributions of this article is to derive the closed-form solutions of Rayleigh wave motions in an orthotropic layered half-space under time-harmonic loads using reciprocity theorems. The reciprocity approach was successfully employed to calculate Rayleigh wave fields in a layered half-space where both the layer and the half-space were homogeneous and isotropic solids.²² Preliminary results on the progress involved in the current work have also been reported by the authors.²⁶ The rest of the paper is organized into four sections. The main purpose of Section 2 is to propose explicit expressions of Rayleigh waves propagating in an orthotropic layered half-space. Section 3 discusses the procedure to obtain the closed-form solutions of Rayleigh wave motion due to various time-harmonic loads by application of reciprocity. Section 4 presents several calculations of Rayleigh wave motions in a jointed rock including a quartz-schist layer and a soil half-space.

Moreover, we present the phase and group dispersion curves superimposed with the amplitude spectra that provide useful information on wave modes, frequencies, and displacement amplitudes. The analytical predictions are compared with numerical results found by finite element analysis, and they show excellent agreement for several loading conditions. The calculations obtained in the current study could generally be beneficial for applications in seismology and materials characterization of coated structures. Major conclusions drawn from the present work are given in Section 5.

Propagation of free Rayleigh waves in an orthotropic layered half-space

We propose here the explicit solutions of Rayleigh surface waves traveling in an orthotropic half-space coated by an orthotropic thin layer. These solutions are also important for the computation of Rayleigh wave fields in the layered half-space due to the application of time-harmonic sources using the reciprocity theorem in the following section. Let us begin by considering a layered structure including a layer $Ω$ of thickness h perfectly bonded with a half-space $\hat{Ω}$ along $z = 0$ of the coordinate system $(x, z)$ as shown in Figure 1. Both the layer and the half-space are of homogeneous orthotropic elastic solids with the properties indicated by the stiffness matrices C and $\hat{C}$ , respectively. The governing equation for elastic waves traveling in an orthotropic medium can be written as⁷

σ_{i j, i} = ρ {\ddot{u}}_{j} i, j = x, z

(1)

where

σ_{i j}

represent stresses,

u_{j}

indicate displacements, and

ρ

denotes mass density. Equation (1) describes the relationship between stresses and displacements without the presence of force terms. This equation is therefore used to illustrate the propagation of free Rayleigh waves.

Figure 1.

Free Rayleigh waves in a layered half-space.

In general, guided waves that propagate in a thin layer can be separated into six partial waves, including three incident waves $P^{+}, S V^{+}, S H^{+}$ and three reflected waves $P^{-}, S V^{-}, S H^{-}$ .²⁷ Meanwhile, there are only three partial waves in a half-space since there is no reflection back to the surface of the half-space.⁷ In this study, we are interested in a two-dimensional plane strain problem so that the $S H$ waves are ignored. The displacement fields of the layer may then be expressed as a superposition of four partial waves

u_{x} = (A_{1} e^{i k α_{1} z} + A_{2} e^{i k α_{2} z} + A_{3} e^{i k α_{3} z} + A_{4} e^{i k α_{4} z}) e^{i k (x - c t)}

(2)

u_{z} = (W_{1} A_{1} e^{i k α_{1} z} + W_{2} A_{2} e^{i k α_{2} z} + W_{3} A_{3} e^{i k α_{3} z} + W_{4} A_{4} e^{i k α_{4} z}) e^{i k (x - c t)} .

(3)

Two partial waves are involved in the half-space, thus

{\hat{u}}_{x} = ({\hat{A}}_{1} e^{- k {\hat{α}}_{1} z} + {\hat{A}}_{2} e^{- k {\hat{α}}_{2} z}) e^{i k (x - c t)}

(4)

{\hat{u}}_{z} = i ({\hat{W}}_{1} {\hat{A}}_{1} e^{- k {\hat{α}}_{1} z} + {\hat{W}}_{2} {\hat{A}}_{2} e^{- k {\hat{α}}_{2} z}) e^{i k (x - c t)}

(5)

where k is wavenumber in the x -direction, i is the imaginary unit, c denotes phase velocity while

A_{j} (j = 1, 2, 3, 4)

and

{\hat{A}}_{j} (j = 1, 2)

are unknown amplitudes of the partial waves in the layer and the half-space, respectively. The expressions of dimensionless quantities

W_{j}, {\hat{W}}_{j}

and parameters

α_{j}, {\hat{α}}_{j}

which represent the ratios of the wavenumbers between z -direction and x -direction are given in the Appendix.

Note that the corresponding stress components $σ_{x x}, σ_{x z}, σ_{z z}$ can be easily calculated using Hooke's law. By assuming a perfectly bonded boundary condition at the interface of the layer and the half-space, we may write

σ_{z z} = 0; σ_{x z} = 0 a t z = - h

u_{x} = {\hat{u}}_{x}; u_{z} = {\hat{u}}_{z}; σ_{z z} = {\hat{σ}}_{z z}; σ_{x z} = {\hat{σ}}_{x z} a t z = 0

(6)

Substituting the displacements and stresses into the boundary conditions yields,

T A = 0

(7)

where

A = [A_{1} A_{2} A_{3} A_{4} {\hat{A}}_{1} {\hat{A}}_{2}]^{T}

and

T

is a six-by-six matrix whose expression is provided in the Appendix. For nontrivial solutions of coefficients in

A

, the determinant of matrix

T

must vanish, i.e.

det (T) = 0

, leading to a relation called the characteristic equation of Rayleigh waves in an orthotropic layered half-space. Solving this equation, which is generally complicated and requires a numerical procedure, results in dispersion curves shown in Figure 3.

Figure 3.

Dispersion curves of Rayleigh waves in a rock mass model.

Note that, in Eq. (7), there are only five independent equations with six unknowns $A_{1}$ , $A_{2}$ , $A_{3}$ , $A_{4}$ , ${\hat{A}}_{1}$ and ${\hat{A}}_{2}$ . We find a general solution in the form of

A_{1} = A d_{1}; A_{2} = A d_{2}; A_{3} = A d_{3}; A_{4} = A d_{4}; {\hat{A}}_{1} = A {\hat{d}}_{1}; {\hat{A}}_{2} = A {\hat{d}}_{2}

(8)

where A represents the only relative amplitude while

d_{1}

d_{2}

d_{3}

d_{4}

{\hat{d}}_{1}

and

{\hat{d}}_{2}

are dimensionless coefficients depending on material properties of the layered half-space (see the Appendix).

For the layer, the displacement and stress fields can now be rewritten as

u_{x} = A U_{x} (z) e^{i k (x - c t)}, u_{z} = A U_{z} (z) e^{i k (x - c t)}

(9)

σ_{x x} = i k C_{55} A T_{x x} (z) e^{i k (x - c t)}, σ_{x z} = i k C_{55} A T_{x z} (z) e^{i k (x - c t)}

(10)

where

\begin{aligned} U_{x} (z) = d_{1} e^{i k α_{1} z} + d_{2} e^{i k α_{2} z} + d_{3} e^{- i k α_{1} z} + d_{4} e^{- i k α_{2} z}, \\ U_{z} (z) = W_{1} d_{1} e^{i k α_{1} z} + W_{2} d_{2} e^{i k α_{2} z} - W_{1} d_{3} e^{- i k α_{1} z} - W_{2} d_{4} e^{- i k α_{2} z} \end{aligned}

(11)

\begin{aligned} T_{x x} (z) = \frac{D_{31}}{C_{55}} d_{1} e^{i k α_{1} z} + \frac{D_{32}}{C_{55}} d_{2} e^{i k α_{2} z} + \frac{D_{31}}{C_{55}} d_{3} e^{- i k α_{1} z} + \frac{D_{32}}{C_{55}} d_{4} e^{- i k α_{2} z}, \\ T_{x z} (z) = (α_{1} + W_{1}) d_{1} e^{i k α_{1} z} + (α_{2} + W_{2}) d_{2} e^{i k α_{2} z} - (α_{1} + W_{1}) d_{3} e^{- i k α_{1} z} - (α_{2} + W_{2}) d_{4} e^{- i k α_{2} z} . \end{aligned}

(12)

For the half-space,

{\hat{u}}_{x} = A {\hat{U}}_{x} (z) e^{i k (x - c t)}, {\hat{u}}_{z} = i A {\hat{U}}_{z} (z) e^{i k (x - c t)}

(13)

{\hat{σ}}_{x x} = i k {\hat{C}}_{55} A {\hat{T}}_{x x} (z) e^{i k (x - c t)}, {\hat{σ}}_{x z} = - k {\hat{C}}_{55} A {\hat{T}}_{x z} (z) e^{i k (x - c t)}

(14)

where

{\hat{U}}_{x} (z) = {\hat{d}}_{1} e^{- k {\hat{α}}_{1} z} + {\hat{d}}_{2} e^{- k {\hat{α}}_{2} z}, {\hat{U}}_{z} (z) = {\hat{W}}_{1} {\hat{d}}_{1} e^{- k {\hat{α}}_{1} z} + {\hat{W}}_{2} {\hat{d}}_{2} e^{- k {\hat{α}}_{2} z}

(15)

\begin{aligned} {\hat{T}}_{x x} (z) = \frac{{\hat{D}}_{31}}{{\hat{C}}_{55}} {\hat{d}}_{1} e^{- k {\hat{α}}_{1} z} + \frac{{\hat{D}}_{32}}{{\hat{C}}_{55}} {\hat{d}}_{2} e^{- k {\hat{α}}_{2} z}, \\ {\hat{T}}_{x z} (z) = ({\hat{α}}_{1} + {\hat{W}}_{1}) {\hat{d}}_{1} e^{- k {\hat{α}}_{1} z} + ({\hat{α}}_{2} + {\hat{W}}_{2}) {\hat{d}}_{2} e^{- k {\hat{α}}_{2} z} . \end{aligned}

(16)

Here, the terms

C_{i j}

denote components of material stiffness matrix and the quantities

D_{i j}

can be found in the Appendix. The Rayleigh wave fields are now expressed via only one unknown amplitude A instead of six, as in Eqs. (2)–(5). This amplitude depends on the loading and will be computed using the reciprocity theorem in the following section.

Reciprocity application for the calculation of Rayleigh wave motions

This section aims to derive closed-form solutions of Rayleigh waves in an orthotropic layered half-space subjected to a time-harmonic source. It will be shown that the amplitudes of generated Rayleigh waves are obtained, in a simple manner, using reciprocity relations between two loading states of an elastic body expressed as¹⁸

\int_{Ω}^{} (f_{j}^{A} u_{j}^{A} - f_{j}^{B} u_{j}^{B}) d Ω + \int_{\hat{Ω}}^{} ({\hat{f}}_{j}^{A} {\hat{u}}_{j}^{A} - {\hat{f}}_{j}^{B} {\hat{u}}_{j}^{B}) d \hat{Ω} = \int_{S}^{} (σ_{i j}^{B} u_{j}^{A} - σ_{i j}^{A} u_{j}^{B}) n_{i} d S + \int_{\hat{S}}^{} ({\hat{σ}}_{i j}^{B} {\hat{u}}_{j}^{A} - {\hat{σ}}_{i j}^{A} {\hat{u}}_{j}^{B}) {\hat{n}}_{i} d \hat{S}

(17)

Here, S and

\hat{S}

are external boundaries of the layer and the half-space without the interface, respectively, n and

\hat{n}

are normal vectors, and f represents body force, see Figure 2. The superscripts

A

and

B

stand for state

A

(actual state) and state

B

(virtual state), respectively. For this problem, State

A

is the Rayleigh wave motion due to the loading, while state

B

is generally chosen as a free Rayleigh wave propagating in the layered structure.

Figure 2.

Layered half-space under a time-harmonic load.

We first consider a vertical load applied at $(x_{0}, z_{0})$ in the form of Delta function as

f_{z}^{A} = P δ (z - z_{0}) δ (x - x_{0}) e^{- i k c t}

(18)

Under the loading, both body waves and Rayleigh waves are generated. The body waves, however, rapidly attenuate after several wavelengths²⁸ making the Rayleigh waves more dominant. In the far-field, the motion of Rayleigh waves propagating in the positive x -direction can be expressed as a summation of wave modes as

u_{x}^{m} = \sum_{m = 0}^{\infty} A_{m}^{P +} U_{x}^{m} (z) e^{i k_{m} (x - c_{m} t)}, u_{z}^{m} = \sum_{m = 0}^{\infty} A_{m}^{P +} U_{z}^{m} (z) e^{i k_{m} (x - c_{m} t)}

(19)

σ_{x x}^{m} = \sum_{m = 0}^{\infty} i k_{m} C_{55} A_{m}^{P +} T_{x x}^{m} (z) e^{i k_{m} (x - c_{m} t)}, σ_{x z}^{m} = \sum_{m = 0}^{\infty} i k_{m} C_{55} A_{m}^{P +} T_{x z}^{m} (z) e^{i k_{m} (x - c_{m} t)}

(20)

and

{\hat{u}}_{x}^{m} = \sum_{m = 0}^{\infty} A_{m}^{P +} {\hat{U}}_{x}^{m} (z) e^{i k_{m} (x - c_{m} t)}, {\hat{u}}_{z}^{m} = \sum_{m = 0}^{\infty} i A_{m}^{P +} {\hat{U}}_{z}^{m} (z) e^{i k_{m} (x - c_{m} t)}

(21)

{\hat{σ}}_{x x}^{m} = \sum_{m = 0}^{\infty} i k_{m} {\hat{C}}_{55} A_{m}^{P +} {\hat{T}}_{x x}^{m} (z) e^{i k_{m} (x - c_{m} t)}, {\hat{σ}}_{x z}^{m} = \sum_{m = 0}^{\infty} - k_{m} {\hat{C}}_{55} A_{m}^{P +} {\hat{T}}_{x z}^{m} (z) e^{i k_{m} (x - c_{m} t)}

(22)

where

A_{m}^{P +}

and

k_{m}

are the amplitude and wavenumber of mode m, respectively, while

U_{x}^{m}

U_{z}^{m}

T_{x x}^{m}

T_{x z}^{m}

and

{\hat{U}}_{x}^{m}

{\hat{U}}_{z}^{m}

{\hat{T}}_{x x}^{m}

{\hat{T}}_{x z}^{m}

are the functions which have been defined previously in Eqs. (11)–(12) and Eqs. (15)–(16).

We choose virtual state $B$ as a free Rayleigh wave of mode n traveling in negative x -direction written as

u_{x}^{n} = - B_{n} U_{x}^{n} (z) e^{- i k_{n} (x + c_{n} t)}, u_{z}^{n} = B_{n} U_{z}^{n} (z) e^{- i k_{n} (x + c_{n} t)}

(23)

σ_{x x}^{n} = i k_{n} C_{55} B_{n} T_{x x}^{n} (z) e^{- i k_{n} (x + c_{n} t)}, σ_{x z}^{n} = - i k_{n} C_{55} B_{n} T_{x z}^{n} (z) e^{- i k_{n} (x + c_{n} t)}

(24)

and

{\hat{u}}_{x}^{n} = - B_{n} {\hat{U}}_{x}^{n} (z) e^{- i k_{n} (x + c_{n} t)}, {\hat{u}}_{z}^{n} = i B_{n} {\hat{U}}_{z}^{n} (z) e^{- i k_{n} (x + c_{n} t)}

(25)

{\hat{σ}}_{x x}^{n} = i k_{n} {\hat{C}}_{55} B_{n} {\hat{T}}_{x x}^{n} (z) e^{- i k_{n} (x + c_{n} t)}, σ_{x z}^{n} = k_{n} C_{55} B_{n} T_{x z}^{n} (z) e^{- i k_{n} (x + c_{n} t)}

(26)

where

B_{n}

imply the amplitudes of free Rayleigh mode n.

Substituting these expressions into Eq. (17) yields

P B_{n} U_{z}^{n} (z_{0}) e^{- i k_{n} (x_{0} + c_{n} t)} e^{- i k_{n} c_{n} t} = \sum_{m = 0}^{\infty} (\int_{- h}^{0} F_{A B}^{m n} |_{x = b} d z + \int_{0}^{\infty} {\hat{F}}_{A B}^{m n} |_{x = b} d z)

(27)

where

F_{A B}^{m n} = σ_{x x}^{B n} u_{x}^{A m} + σ_{x z}^{B n} u_{z}^{A m} - σ_{x x}^{A m} u_{x}^{B n} - σ_{x z}^{A m} u_{z}^{B n}

(28)

{\hat{F}}_{A B}^{m n} = {\hat{σ}}_{x x}^{B n} {\hat{u}}_{x}^{A m} + {\hat{σ}}_{x z}^{B n} {\hat{u}}_{z}^{A m} - {\hat{σ}}_{x x}^{A m} {\hat{u}}_{x}^{B n} - {\hat{σ}}_{x z}^{A m} {\hat{u}}_{z}^{B n}

(29)

Here the integral along

x = a

is neglected due to the counter-propagating wave property.¹⁹ After some manipulations, Eq. (27) can be simplified as

P U_{z}^{n} (z_{0}) e^{- i k_{n} x_{0}} = \sum_{m = 0}^{\infty} A_{m}^{P +} e^{i (k_{m} - k_{n}) b} (C_{55} I_{m n} + {\hat{C}}_{55} {\hat{I}}_{m n})

(30)

where

I_{m n} = \int_{- h}^{0} i [k_{n} (T_{x x}^{n} (z) U_{x}^{m} (z) - T_{x z}^{n} (z) U_{z}^{m} (z)) + k_{m} (T_{x x}^{m} (z) U_{x}^{n} (z) - T_{x z}^{m} (z) U_{z}^{n} (z))] d z

(31)

{\hat{I}}_{m n} = \int_{0}^{\infty} i [k_{n} ({\hat{T}}_{x x}^{n} (z) {\hat{U}}_{x}^{m} (z) - {\hat{T}}_{x z}^{n} (z) {\hat{U}}_{z}^{m} (z)) + k_{m} ({\hat{T}}_{x x}^{m} (z) {\hat{U}}_{x}^{n} (z) - {\hat{T}}_{x z}^{m} (z) {\hat{U}}_{z}^{n} (z))] d z .

(32)

The integrals

I_{m n}

and

{\hat{I}}_{m n}

only yield non-zero values in the case of

m = n

because of the orthogonality relation, which has been meticulously discussed by Phan et al.²² The amplitude of Rayleigh waves is then obtained as

A_{m}^{P +} = \frac{P U_{z}^{n} (z_{0}) e^{- i k_{n} x_{0}}}{2 (C_{55} I_{n n} + {\hat{C}}_{55} {\hat{I}}_{n n})}

(33)

where

I_{n n} = i k_{n} \int_{- h}^{0} (T_{x x}^{n} (z) U_{x}^{n} (z) - T_{x z}^{n} (z) U_{z}^{n} (z)) d z

(34)

{\hat{I}}_{n n} = i k_{n} \int_{0}^{\infty} ({\hat{T}}_{x x}^{n} (z) {\hat{U}}_{x}^{n} (z) - {\hat{T}}_{x z}^{n} (z) {\hat{U}}_{z}^{n} (z)) d z .

(35)

If we choose state

B

in the positive x -direction, the amplitude is calculated as

A_{m}^{P -} = \frac{P U_{z}^{n} (z_{0}) e^{i k_{n} x_{0}}}{2 (C_{55} I_{n n} + {\hat{C}}_{55} {\hat{I}}_{n n})} .

(36)

Consider now Rayleigh wave motion due to a horizontal load of the form

f_{z}^{A} = Q δ (z - z_{0}) δ (x - x_{0}) e^{- i k c t}

(37)

Following a similar procedure as illustrated using the vertical load, we find the amplitudes of Rayleigh waves as

A_{m}^{Q +} = \frac{- Q U_{x}^{n} (z_{0}) e^{- i k_{n} x_{0}}}{2 (C_{55} I_{n n} + {\hat{C}}_{55} {\hat{I}}_{n n})}, A_{m}^{Q -} = \frac{Q U_{x}^{n} (z_{0}) e^{i k_{n} x_{0}}}{2 (C_{55} I_{n n} + {\hat{C}}_{55} {\hat{I}}_{n n})} .

(38)

Based on the calculations of Rayleigh wave motion due to a time-harmonic load, we may derive solutions of Rayleigh waves generated by a distribution of loadings. In the case of a uniformly distributed vertical load applied on a strip from

x = - a

x = a

, amplitudes of Rayleigh waves in the positive and negative x -directions are computed as

A_{m}^{U +} = \frac{2 A_{m}^{P +} \sin (k a)}{k}, A_{m}^{U -} = \frac{2 A_{m}^{P -} \sin (k a)}{k} .

(39)

In case of varying vertical loadings of the form

\bar{p} (x_{0}) = P e^{i k x_{0}}

on a strip from

x = - a

x = a

, the amplitudes in the positive and negative

x -

directions are

A_{m}^{V +} = 2 a A_{m}^{P +}, A_{m}^{V -} = \frac{A_{m}^{P -} \sin (2 k a)}{k} .

(40)

Results and discussions

The analytical solutions derived in the earlier sections are illustrated here by computing the Rayleigh wave motion in a rock mass modeled as a quartz-schist layer overlaid a soil half-space. Material properties of the layer and the half-space are given in Table 1.¹ The layer thickness is chosen as $h = 1000 m$ . In order to obtain the phase and group velocity dispersion curves of Rayleigh waves in the joint structure via Eq. (7), a computer code is built based on the numerical root-searching algorithms and extrapolation schemes.²⁹ The displacement amplitudes of the Rayleigh waves generated by vertical and horizontal sources of unit magnitude applied on the free surface are calculated using our proposed method. These results are then integrated into the dispersion curves to highlight the amplitude spectra of different wave modes.

Table 1.

Material properties of the quartz-schist layer and the soil half-space.

Material	$C_{11}$ $(G P a)$	$C_{13}$ $(G P a)$	$C_{33}$ $(G P a)$	$C_{55}$ $(G P a)$	$ρ$ $(k g / m^{3})$
Quartz-schist	$13.55$	$3.00$	$12.36$	$3.75$	$2660$
Soil	$201.92$	$86.54$	$201.92$	$57.69$	$2000$

The phase and group velocity dispersion curves superimposed by the displacement amplitudes are plotted in Figure 3. Here, Figure 3(a) and Figure 3(c) show the horizontal-amplitude spectra while Figure 3(b) and Figure 3(d) show the vertical-amplitude spectra. There are upper and lower bounds for the phase velocities in the layered half-space shown in Figure 3(a) and Figure 3(b). In particular, as the thickness of the quartz-schist layer becomes small, the phase velocity of the first Rayleigh mode approaches the value of the Rayleigh wave velocity of the soil half-space $({\hat{c}}_{R})$ . The transverse wave velocity of the soil half-space $({\hat{c}}_{T})$ defines the upper bound. This threshold velocity is a consequence of the existing condition of Rayleigh waves in the half-space that satisfies $0 < ρ c^{2} < min ({\hat{C}}_{11}, {\hat{C}}_{55})$ .³⁰ The Rayleigh wave phase velocity of the quartz-schist layer $(c_{R})$ serves as the lower bound. When the dimensionless quantity $f h$ increases to infinity, all dispersion curves essentially become asymptotic to the lower limit.

With the presence of the amplitude spectra, the curves produced in this study are considerably superior to other dispersion diagrams available in the literature. They facilitate the analysis of seismic waves, such as source localization or earthquake prediction, based on the observed wave amplitudes. In ultrasonic nondestructive evaluation, the superimposed dispersion curves could offer instant optimal selections of wave modes and frequencies, increasing the signal-to-noise ratio and simplifying signal-processing stages.

Analytical predictions of Rayleigh wave fields presented above are now verified using a finite element simulation carried out by COMSOL Multiphysics software. In this benchmark test, we compute the amplitudes of the Rayleigh waves generated by both a uniform distributed load and a varying distributed load applied over a strip of width $2 a$ on the surface of the $1000 m$ quartz-schist layer. By the use of a $0.1 H z$ vertical distributed load, the lowest mode of the Rayleigh waves is triggered.

The horizontal displacement amplitudes $(u_{1})$ and the vertical displacement amplitudes $(u_{2})$ of Rayleigh waves are shown in Figure 4 for the uniform distributed loading and Figure 5 for the varying distributed loading. It can be seen in Figure 4 that both displacement amplitudes $u_{1}, u_{2}$ first increase and then decrease with increase of the width $2 a$ because of the presence of the term $\sin k a$ in Eq. (39). For this type of loading, the displacement amplitudes in the positive and negative directions are the same because this is a symmetric problem. On the other hand, the amplitudes of the Rayleigh surface waves due to the varying load in the positive and negative directions are different in Figure 5. Here, the displacement amplitudes in the positive direction are linearly proportional to the strip width where the loads are applied. This can be explained by Eq. 40. The amplitudes in the negative direction, which depend on the term $\sin 2 k a$ , increase at low values of $k a$ and then decrease with the increase of the loading width $2 a$ .

Figure 4.

Uniform strip load of unit magnitude. Numerical results: $u_{2}$ – circle sign; $u_{1}$ – asterisk sign. Analytical results: solid lines.

Figure 5.

Varying strip load of unit magnitude. Numerical results: $u_{2}^{+}$ – circle sign; $u_{1}^{+}$ – asterisk sign; $u_{2}^{+}$ – multiply sign; $u_{1}^{-}$ – plus sign. Analytical results: solid and dash lines.

In both cases of loading conditions, comparisons between the analytical prediction and numerical results are in excellent agreement. Small differences appearing in the comparisons are accounted for by the limitation of the numerical approach in dealing with the problems of multi-mode and dispersive guided waves. Due to the closed-form nature of solutions, analytical calculations are performed swiftly in these testing cases, while numerical simulation takes significant time to reproduce these results. Thus, the proposed method could be appropriate for solving inverse problems such as material characterization or defect quantification, which typically necessitate a large amount of data or a fast forward-problem solver. Furthermore, the most intriguing aspect of our approach is that reciprocity calculations can be feasible for individual modes at high frequencies. Thus, the method can perform effectively in highly dispersive regions with numerous overlapping modes, which could be challenging for numerical simulations.

Conclusions

The motion of Rayleigh waves in an orthotropic half-space coated by a thin orthotropic layer due to the application of time-harmonic loadings is discussed in this article. The explicit expressions of Rayleigh surface waves propagating in a layered half-space and that depend on only one unknown amplitude have been introduced. The main contribution is the derivation of the analytical solutions of far-field Rayleigh waves using elastodynamic reciprocity theorems. The proposed approach has shown an advantage for obtaining the Rayleigh wave fields in a simple manner and in closed-form exact solutions. Computation of Rayleigh wave motion in a jointed rock mass, including a quartz-schist layer and a soil half-space, has been demonstrated for the cases of a uniform distributed load and a varying distributed load both applied over a strip on the surface of the quartz-schist layer. In particular, we have presented the dispersion curves superimposed with the amplitude spectra which could offer useful information on Rayleigh wave modes, frequencies, and displacement amplitudes. The comparisons between the analytical predictions and numerical results by finite element analysis have shown excellent agreement. With the closed-form property and applicability to complex structures, the expressions and results reported in this research might be beneficial for application in seismology and materials characterization of coated half-space. Further studies on the propagation of Rayleigh waves in multilayered half-space can be conducted using the reciprocity theorem. These results can serve as fundamental solutions for scattering problems of Rayleigh waves due to a cavity in structures.

Footnotes

Acknowledgments

This research is funded by Graduate University of Science and Technology under grant number GUST.STS.ĐT2020-CH01; and Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 107.02-2019.21.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Vietnam National Foundation for Science and Technology Development (NAFOSTED), Graduate University of Science and Technology, (grant number 107.02-2019.21, GUST.STS.ÐT2020-CH01).

Author biography

Ductho Le received his master's degree in Solid Mechanics from VNU University of Science, Hanoi. His research interest includes bone quantitative ultrasound, structural health monitoring, and material characterization.

Yen Nguyen received his PhD in Engineering Physics from Hanoi University of Science and Technology, Vietnam in 2018. His research focuses on Molecular Dynamic Simulation and Physical Optics. He is currently a senior lecturer & researcher at Institute of Theoretical and Applied Research, Duy Tan University.

Hoang Ngoc Quy graduated with a master's degree in Mechanical Engineering from VNU University of Engineering and Technology, Hanoi. His research focuses on solid-body mechanics, vibrations of sound, and ultrasonic waves.

Lawrence H. Le received his PhD degree in Earth Physics and MBA degree in Finance and Technology Commercialization from the University of Alberta, Edmonton, AB, Canada, in 1991 and 1999, respectively. He held a Natural Sciences and Engineering Research Council of Canada (NSERC) postdoctoral fellowship at the Schlumberger-Doll Research Laboratory, Ridgefield, CT, USA. He is currently a Clinical Professor leading the graduate program in the Department of Radiology and Diagnostic Imaging and a Senior Medical Physicist with Alberta Health Services, Edmonton, AB. He is also a Senior Visiting Scholar with the Center for Biomedical Engineering, Fudan University, Shanghai, China.

Haidang Phan received his PhD from Pusan National University, South Korea, majoring in Mechanical Engineering in 2013. He has made original and significant contributions to both theoretical and applied studies of guided waves in aerospace composites and in cortical and/or cancellous bones for characterization of their material properties. Currently working as a senior lecturer at the University of Engineering and Technology, Vietnam National University, Dr. Phan published two monographs, several book chapters, along with more than thirty journal papers, and is the principal investigator of several national-funded projects.

Appendix

Expressions of $W_{j}$ and ${\hat{W}}_{j}$ :

W_{j} = \frac{ρ c^{2} - C_{11} - C_{55} α_{j}^{2}}{(C_{13} + C_{55}) α_{j}}, {\hat{W}}_{j} = \frac{{\hat{C}}_{11} - {\hat{C}}_{55} {\hat{α}}_{j}^{2} - \hat{ρ} c^{2}}{({\hat{C}}_{13} + {\hat{C}}_{55}) {\hat{α}}_{j}}, j = 1, 2

where

C_{i j}

and

{\hat{C}}_{i j}

are components of the stiffness matrices C and

\hat{C}

, respectively.

Expressions of $α_{j}$ and ${\hat{α}}_{j}$ :

α_{1} = - α_{3} = \sqrt{\frac{S + \sqrt{S^{2} - 4 P}}{2}}, α_{2} = - α_{4} = \sqrt{\frac{S - \sqrt{S^{2} - 4 P}}{2}}

{\hat{α}}_{1}^{2} + {\hat{α}}_{2}^{2} = \hat{S}, {\hat{α}}_{1}^{2} {\hat{α}}_{2}^{2} = \hat{P}

where

S = \frac{{(C_{13} + C_{55})}^{2} - (C_{11} - ρ c^{2}) C_{33} - (C_{55} - ρ c^{2}) C_{55}}{C_{33} C_{55}}, P = \frac{(C_{11} - ρ c^{2}) (C_{55} - ρ c^{2})}{C_{33} C_{55}}

\hat{S} = \frac{({\hat{C}}_{11} - \hat{ρ} c^{2}) {\hat{C}}_{33} + ({\hat{C}}_{55} - \hat{ρ} c^{2}) {\hat{C}}_{55} - {({\hat{C}}_{13} + {\hat{C}}_{55})}^{2}}{{\hat{C}}_{33} {\hat{C}}_{55}}, \hat{P} = \frac{({\hat{C}}_{11} - \hat{ρ} c^{2}) ({\hat{C}}_{55} - \hat{ρ} c^{2})}{{\hat{C}}_{33} {\hat{C}}_{55}}

Matrix

T

T = [\begin{matrix} D_{11} e^{- i k α_{1} h} & D_{12} e^{- i k α_{2} h} & D_{11} e^{i k α_{1} h} & D_{12} e^{i k α_{2} h} & 0 & 0 \\ D_{21} e^{- i k α_{1} h} & D_{22} e^{- i k α_{2} h} & - D_{21} e^{i k α_{1} h} & - D_{22} e^{i k α_{2} h} & 0 & 0 \\ 1 & 1 & 1 & 1 & - 1 & - 1 \\ W_{1} & W_{2} & - W & - W_{2} & - i {\hat{W}}_{1} & - i {\hat{W}}_{2} \\ D_{11} & D_{12} & D_{11} & D_{12} & - {\hat{D}}_{11} & - {\hat{D}}_{12} \\ D_{21} & D_{22} & - D_{21} & - D_{22} & - i {\hat{D}}_{21} & - i {\hat{D}}_{22} \end{matrix}]

where

D_{1 j} = C_{13} + C_{33} α_{j} W_{j}, D_{2 j} = C_{55} (α_{j} + W_{j}), D_{3 j} = C_{11} + C_{13} α_{j} W_{j}, j = 1, 2

{\hat{D}}_{1 j} = {\hat{C}}_{13} - {\hat{C}}_{33} {\hat{α}}_{j} {\hat{W}}_{j}, {\hat{D}}_{2 j} = {\hat{C}}_{55} ({\hat{α}}_{j} + {\hat{W}}_{j}), {\hat{D}}_{3 j} = {\hat{C}}_{11} - {\hat{C}}_{13} {\hat{α}}_{j} {\hat{W}}_{j}, j = 1, 2

Expressions of $d_{1}$ , $d_{2}$ , $d_{3}$ , $d_{4}$ , ${\hat{d}}_{1}, {\hat{d}}_{2}$ :

d_{1} = H_{11} {\hat{d}}_{1} + H_{12} {\hat{d}}_{2}, d_{2} = H_{21} {\hat{d}}_{1} + H_{22} {\hat{d}}_{2}

d_{3} = H_{31} {\hat{d}}_{1} + H_{32} {\hat{d}}_{2}, d_{4} = H_{41} {\hat{d}}_{1} + H_{42} {\hat{d}}_{2}

{\hat{d}}_{1} = \frac{D_{11}}{s_{1}} (H_{12} + s_{1}^{2} H_{32}) + \frac{D_{12}}{s_{2}} (H_{22} + s_{2}^{2} H_{42})

{\hat{d}}_{2} = - \frac{D_{11}}{s_{1}} (H_{11} + s_{1}^{2} H_{31}) - \frac{D_{12}}{s_{2}} (H_{21} + s_{2}^{2} H_{41})

where

s_{1} = e^{i k α_{1} h}, s_{2} = e^{i k α_{2} h}

H_{11} = \frac{{\hat{D}}_{11} - D_{12}}{2 (D_{11} - D_{12})} + i \frac{{\hat{D}}_{21} W_{2} - D_{22} {\hat{W}}_{1}}{2 (D_{21} W_{2} - D_{22} W_{1})}

H_{12} = \frac{{\hat{D}}_{12} - D_{12}}{2 (D_{11} - D_{12})} + i \frac{{\hat{D}}_{22} W_{2} - D_{22} {\hat{W}}_{2}}{2 (D_{21} W_{2} - D_{22} W_{1})}

H_{21} = \frac{D_{11} - {\hat{D}}_{11}}{2 (D_{11} - D_{12})} + i \frac{D_{21} {\hat{W}}_{1} - {\hat{D}}_{21} W_{1}}{2 (D_{21} W_{2} - D_{22} W_{1})}

H_{22} = \frac{D_{11} - {\hat{D}}_{12}}{2 (D_{11} - D_{12})} + i \frac{D_{21} {\hat{W}}_{2} - {\hat{D}}_{22} W_{1}}{2 (D_{21} W_{2} - D_{22} W_{1})}

H_{31} = \frac{{\hat{D}}_{11} - D_{12}}{2 (D_{11} - D_{12})} - i \frac{{\hat{D}}_{21} W_{2} - D_{22} {\hat{W}}_{1}}{2 (D_{21} W_{2} - D_{22} W_{1})}

H_{32} = \frac{{\hat{D}}_{12} - D_{12}}{2 (D_{11} - D_{12})} - i \frac{{\hat{D}}_{22} W_{2} - D_{22} {\hat{W}}_{2}}{2 (D_{21} W_{2} - D_{22} W_{1})}

H_{41} = \frac{D_{11} - {\hat{D}}_{11}}{2 (D_{11} - D_{12})} - i \frac{D_{21} {\hat{W}}_{1} - {\hat{D}}_{21} W_{1}}{2 (D_{21} W_{2} - D_{22} W_{1})}

H_{42} = \frac{D_{11} - {\hat{D}}_{12}}{2 (D_{11} - D_{12})} - i \frac{D_{21} {\hat{W}}_{2} - {\hat{D}}_{22} W_{1}}{2 (D_{21} W_{2} - D_{22} W_{1})}

References

Babanouri

Fattahi

. Evaluating orthotropic continuum analysis of stress wave propagation through a jointed rock mass. Bull Eng Geol Environ 2018; 77: 725–733.

Gupta

Ahmed

. On Rayleigh waves in self-reinforced layer embedded over an incompressible half-space with varying rigidity and density. Procedia Eng 2017; 173: 1021–1028.

Pal

Kumar

Bose

. Propagation of Rayleigh waves in anisotropic layer overlying a semi-infinite sandy medium. Ain Shams Eng J 2015; 6: 621–627.

Mohseni

C-T

. Rayleigh Wave propagation and scattering characteristics at debondings in fibre-reinforced polymer-retrofitted concrete structures. Struct Health Monit 2018; 18: 303–317.

Ebrahiminejad

Mardanshahi

Kazemirad

. Nondestructive evaluation of coated structures using lamb wave propagation. Appl Acoust 2022; 185: 108378.

Abo-Dahab SM, Kilany AA, Allam MNM, Mohamed RA, Rida SZ. Influence of several fields on Rayleigh waves propagation in a fiber-reinforced orthotropic half-space material under four thermoelastic models. Waves in Random and Complex Media. 2022; 32(5): 2197–2220.

Achenbach

. Wave Propagation in Elastic Solids. Amsterdam: North-Holland Publishing Company, 1973.

Rose

. Ultrasonic Guided Waves in Solid Media. Cambridge: Cambridge University Press, 2014.

Nayfeh

. Wave Propagation in Layered Anisotropic Media: With Applications to Composites. vol 39. Amsterdam: Elsevier; 1995.

10.

Achenbach

Keshava

. Free waves in a plate supported by a semi-infinite Continuum. J Appl Mech 1967; 34: 397–404.

11.

Tiersten

. Elastic surface waves guided by thin films. J Appl Phys 1969; 40: 770–789.

12.

Vinh

Ngoc Anh

. Rayleigh Waves in an orthotropic half-space coated by a thin orthotropic layer with sliding contact. Int J Eng Sci 2014; 75: 154–164.

13.

Kaplunov

Prikazchikov

Sultanova

. Rayleigh-type waves on a coated elastic half-space with a clamped surface. Philos Trans R Soc A: Math Phys Eng Sci 2019; 377: 20190111.

14.

Vasudevan

Mal

. Response of an elastic plate to localized transient sources. J Appl Mech 1985; 52: 356–362.

15.

. Dynamic response of the multilayered half-space medium due to the spatially periodic harmonic moving load. Soil Dyn Earthquake Engn 2022; 157: 107246.

16.

Liu

Zhang

Lin

, et al. Elastodynamic response for the multi-layered transversely isotropic piezoelectric solid subjected to time-harmonic loads. Compos Struct 2016; 153: 755–772.

17.

Bratov

Kaplunov

Lapatsin

, et al. Elastodynamics of a coated half-space under a sliding contact. Math Mech Solids 2022; 27: 1480–1493.

18.

Achenbach

. Reciprocity in Elastodynamics. Cambridge: Cambridge University Press, 2003.

19.

Phan

Cho

Achenbach

. Validity of the reciprocity approach for determination of surface wave motion. Ultrasonics 2013; 53: 665–671.

20.

Kulkarni

Achenbach

. Application of the reciprocity theorem to determine line-load-generated surface waves on an inhomogeneous transversely isotropic half-space. Wave Motion 2008; 45: 350–360.

21.

Balogun

Achenbach

. Surface waves generated by a line load on a half-space with depth-dependent properties. Wave Motion 2013; 50: 1063–1072.

22.

Phan

Cho

, et al. A closed-form solution to propagation of guided waves in a layered half-space under a time-harmonic load: an application of elastodynamic reciprocity. Ultrasonics 2019; 96: 40–47.

23.

Phan

Bui

Nguyen

HTL

, et al. Computation of interface wave motions by reciprocity considerations. Wave Motion 2018; 79: 10–22.

24.

Dao

Ngo

Phan

, et al. Rayleigh Wave motions in an orthotropic half-space under time-harmonic loadings: a theoretical study. Appl Math Model 2020; 87: 171–179.

25.

Phan

Cho

Achenbach

. Verification of surface wave solutions obtained by the reciprocity theorem. Ultrasonics 2014; 54: 1891–1894.

26.

Dao

, et al. Guided wave propagation in a layered half-space structure of anisotropic materials. 2020 5th international conference on green technology and sustainable development (GTSD). 27-28 Nov. 2020 2020:398–404. doi: https://doi.org/10.1109/GTSD50082.2020.9303064.

27.

Solie

Auld

. Elastic waves in free anisotropic plates. J Acoust Soc Am 1973; 54: 50–65.

28.

Phan

Cho

Achenbach

. Application of the reciprocity theorem to scattering of surface waves by a cavity. Int J Solids Struct 2013; 50: 4080–4088.

29.

Lowe

MJS

. Matrix techniques for modeling ultrasonic waves in multilayered media. IEEE Trans Ultrason Ferroelectr Freq Control 1995; 42: 525–542.

30.

Vinh

Ogden

. Formulas for the Rayleigh wave speed in orthotropic elastic solids. Arch Mech 2004; 56: 247–265.