Reflection of plane waves from the free surface of a hard sphere-filled elastic metacomposite

Abstract

We use an effective medium model to study the problem of reflection of plane waves from the free surface of a half-space occupied by an elastic particulate metacomposite. This problem has received little attention in the recent literature despite its significance from both practical and theoretical points of view. Classical formulas for the reflection angles and amplitudes of the reflected waves for a homogeneous elastic half-space with no wave attenuation are extended to a particulate metacomposite half-space with wave attenuation. We also include a detailed discussion concerning the reflected plane shear wave and surface compressional wave in the case of an incident shear wave propagating at an incident angle smaller than the critical angle. The efficiency and accuracy of the model are demonstrated via detailed comparisons between the predicted phase velocity and attenuation coefficient of plane waves in an (infinite) entire space and the corresponding results available in the literature. The implications of our results on the reflection of plane waves from the free surface of a hard sphere-filled elastic metacomposite are discussed. We mention that a quantitative validation of our results cannot be made here as a result of the lack of availability of established data in the existing literature.

Keywords

Reflection of plane waves particulate composite free surface

1. Introduction

Studying the dynamic characteristics of different elastic solids is a key in a number of practical applications including structural vibrations in engineering, ultrasonic material inspection, seismology, and numerous other domains. Typically, such solids can be described effectively using classical dynamic equations from linear isotropic elasticity [1]. However, certain materials with complex microstructure (e.g., certain composite materials [2,3], metacomposites [4 –7], and granular materials [8,9]) exhibit a distinctively characteristic reaction to an applied load. Numerous studies have examined the phenomenon of wave reflection from free surfaces across various materials, considering different boundaries and initial conditions. For example, Singh and Singh [10] studied the reflection of quasi-P and quasi-S waves at the free surface of a fiber-reinforced, anisotropic elastic half-space media and examined how the amplitude ratios change with incident wave angle. Their results showed that reinforcement significantly impacts both the amplitude ratios and the critical angle of the wave. Graff and Pao [11] investigated the plane wave reflection at a stress-free surface of an elastic half-space solid in the presence of couple stresses. They illustrated that by considering couple stresses, three different kinds of waves are reflected from a stress-free surface when the incident waves are either plane pressure or plane shear waves. The additional reflected wave, referred to as the SS wave (represents surface disturbance), quickly decays from the surface for any viable value of the couple stress material constant. Dey and Addy [12] performed an analytical study of plane wave reflection at a free surface of an initially stressed medium. In this study, several results are presented to illustrate how the reflected P and SV waves change with different initial stresses. Guha and Singh [13] examined how the initial stresses affected the reflection properties of plane waves at the free surface of a half-space rotating piezo thermoelastic fiber-reinforced composite. They employed three thermoelasticity theories, namely Green–Lindsay, Lord–Shulman, and classical dynamical coupled theories, to investigate the issue. They demonstrated how the amplitude ratios of reflected waves depend on the angle of incidence, the characteristics of the rotating medium, and different levels of initial stresses. Sidhu and Singh [14] examined the reflection of plane waves at the free surface of an initially stressed elastic medium directly without using potentials. In this research, they demonstrated that the velocities of quasi-P and quasi-S waves vary based on the propagation angle. In addition, the model of the initially stressed medium employed in this particular study is more comprehensive than the one used by Dey and Addy [12]. Singh and Yadav [15] examined how plane waves reflect from the stress-free surface of an electrically conductive and rotating fiber-reinforced elastic solid half-space in the presence of a magnetic field. They numerically calculated the energy ratios of reflected waves and the velocities of plane waves for a specific material. In addition, they showcased the impacts of anisotropy, rotation, and magnetic fields on the velocities and energy ratios of waves. Deswal et al. [16] studied plane wave reflection at the boundary surface of a magneto-thermoelastic medium, considering the impacts of changing thermal conductivity and mass diffusivity. They presented the phase velocity, energy ratio, and amplitude ratio of waves for a specific material. Their findings indicate that the results vary based on the angle of incidence, frequency, and the thermoelastic attributes of the half-space. Yadav [17] investigated the reflection of four types of quasi-planar waves from a thermally insulated stress-free surface of a rotating orthotropic magneto-thermoelastic half-space medium. They calculated energy ratio and reflection coefficients of reflected waves for specific materials and illustrated the influence of rotation angular frequency, diffusion, angle of incidence, and magnetic parameters on the energy ratio and reflection coefficients. Singh and Guha [18] employed the classical dynamic couple theory, Green–Lindsay theory, and Lord–Shulman theory to analyze reflection of plane waves from a free surface of a piezothermoelastic fiber-reinforced composite half-space solid. They illustrated the existence of four different types of waves in the two-dimensional model of the composite medium and computed the amplitude ratio of reflected waves by employing the corresponding boundary condition. Selim [19] studied how the initial stresses influence the reflection coefficients of plane waves in a dissipative medium. They employed Biot’s incremental deformation theory to obtain reflection coefficients of reflected waves. Their research demonstrates that the presence of compressive initial stresses amplifies the speed of P waves while decreasing that of SV waves. Singh et al. [20] examined the propagation of time-harmonic plane waves in an infinite nonlocal elastic solid material containing voids. They demonstrated the impact of voids and nonlocality parameters on the dispersion and attenuation of compressional and transverse waves in the medium. Moreover, they explored the reflection phenomenon of incident coupled compressional waves from the stress-free boundary surface of a nonlocal elastic solid half-space containing voids. By employing suitable boundary conditions, they derived formulas for various reflection coefficients and their corresponding energy ratios. In addition, they visually depicted the impact of nonlocality on reflection coefficients through graphical representations. Furthermore, several studies [21,22] have investigated wave reflection from rough surfaces and examined the effects of multiple scattering on acoustic reflection from rough surfaces. The aforementioned review clearly indicates the intensive and continuing research dedicated to the study of wave reflection phenomena from a free surface in different kinds of materials. However, none of these studies have addressed the phenomenon of wave reflection in particulate elastic metacomposites which represent a new type of composite material displaying remarkable dynamic characteristics, including wave attenuation and vibration isolation. To the best of our knowledge, the study of the problem concerning the reflection of plane waves from the free surface of a particulate composite half-space remains absent from the existing literature. The main contribution of the present research in contrast to other related papers is to study the reflection of plane waves (both compressional and shear) when they are incident at an angle $θ$ on the free boundary surface of a hard sphere-reinforced random composite half-space. Particulate composites form an important category within the group of elastic metacomposites, wherein the elastic modulus and mass density of the embedded spherical particles are notably greater than that of the compliant elastic matrix. A notable advantage of the present model lies in its provision of relatively simple formulas, eliminating the necessity for complex numerical calculations.

Our paper is organized as follows. In section 2, we present the governing equations of wave reflection at the boundary of a metacomposite for two different incident waves. Subsequently, in section 3, we first verify the effectiveness and precision of the wave propagation model presented here by comparing its predicted outcomes to existing, known data in the literature. Next, we use the model to make a comprehensive analysis of the reflection behavior of two different incident waves from the free surface of metacomposites. In section 4, we offer an overview of our findings together with some concluding remarks.

2. Wave reflection on the free surface of a metacomposite half-space

An effective medium model (details given in Appendix 1) is used to investigate the wave reflection of an incident compressional or shear plane wave on the free surface of a hard sphere-filled metacomposite half-space under the conditions of plane strain. In comparison with existing related models, the present model enjoys mathematical simplicity although restricted to heavy hard sphere-filled elastic composites in which the elastic moduli and mass density of embedded spheres are much greater than that of the softer and lighter matrix (such as those listed in Tables 1 and 2).

Table 1.

Characteristic parameters of particulate metacomposite, taken from Sabina and Willis [23].

	Epon828z	Lead
Density ( $ρ$ ) [ $\frac{gr}{c m^{3}}$ ]	$1.202$	$11.3$
P-wave speed ( $c_{p}$ ) [ $\frac{mm}{μ s}$ ]		$2.21$
S-wave speed ( $c_{s}$ ) [ $\frac{mm}{μ s}$ ]	$1.2$	$0.86$
Poisson ratio ( $ν$ ) [ - ]	$0.37$	—
Radius ( $a$ ) [ $μ m$ ]	—	$660$

Table 2.

Characteristic parameters of particulate metacomposite, taken from Luppé et al [24].

	Epoxy	Tungsten carbide
Density ( $ρ$ ) [ $\frac{gr}{c m^{3}}$ ]	$1.142$	$14.95$
Shear modulus ( $G$ )[ $GPa$ ]	$1.6 - 0.06 i$	$273$
Poisson ratio ( $ν$ ) [ - ]	$0.37$	−
Radius ( $a$ ) [ $μ m$ ]	−	$198.5$

We first establish a fixed rectangular coordinate system ${x_{i}} (i = 1, 2, 3)$ . Accordingly, the corresponding displacement field (u₁, u₂, u₃) of a plane wave can be described as having u₁ and u₃ dependent on x₁, x₃, and time t, with u₂ set to zero. Consider the reflection of an incident plane wave on the free surface of a particulate elastic composite occupying the lower half-space x₃ ≥ 0 (as illustrated in Figure 1) where the incident wave with a known frequency and amplitude is advancing toward the free boundary of the half-space in the indicated direction.

Figure 1.

The reflected compressional and shear plane waves of an incident (compressional or shear) plane wave.

At the free boundary of the half-space, the imposed boundary condition requires the traction to be zero, resulting in both the normal stress and the two components of shear stress on the boundary also being zero:

{σ_{13} |}_{x_{3} = 0} = 0, {σ_{23} |}_{x_{3} = 0} = 0, {σ_{33} |}_{x_{3} = 0} = 0 .

(1)

In consideration of the two-dimensional motion, the stress component is zero. In addition, the two remaining boundary conditions provide two requirements that the displacement field must satisfy:

{σ_{13} |}_{x_{3} = 0} = {μ (\frac{\partial u_{1}}{\partial x_{3}} + \frac{\partial u_{3}}{\partial x_{1}}) |}_{x_{3} = 0} = 0,

(2)

{σ_{33} |}_{x_{3} = 0} = {λ \frac{\partial u_{1}}{\partial x_{1}} + (λ + 2 μ) \frac{\partial u_{3}}{\partial x_{3}} |}_{x_{3} = 0} = 0 .

(3)

here, $μ$ and $λ$ are the effective Lame constants of the isotropic metacomposite. It is obvious that the boundary conditions cannot be satisfied if we assume that the incident compressional wave results solely in a reflected compressional wave. The satisfaction of the boundary conditions is only possible when we consider, in addition, the presence of a reflected shear wave.

2.1. An incident compressional plane wave

We first consider the reflected compressional and shear waves of an incident compressional plane wave as shown in Figure 1. For an incident compressional wave shown in Figure 1 with the displacement in the direction ( $x'$ ) of wave propagation given by (the real part of):

u_{I} = I e^{i (k_{P} x' - ω t)},

(4)

where $x' = x_{1} \cos θ - x_{3} \sin θ$ ; $I$ and $k_{p}$ are the (complex-valued) amplitude and wave number of the incident wave; $ω$ , $i$ , and $t$ are, respectively, frequency, imaginary unit, and time. Thus, the displacement components of the incident compressional wave in the $x_{1}$ and $x_{3}$ directions are expressed as:

{(u_{1})}_{I} = I \cos θ e^{i (k_{P} x_{1} \cos θ - k_{P} x_{3} \sin θ - ω t)},

(5)

{(u_{3})}_{I} = - I \sin θ e^{i (k_{P} x_{1} \cos θ - k_{P} x_{3} \sin θ - ω t)} .

(6)

Therefore, the total displacements of the incident compressional wave and the reflected compressional and shear plane waves in the $x_{1}$ and $x_{3}$ directions are written as:

\begin{matrix} u_{1} = {(u_{1})}_{I} + {(u_{1})}_{P} + {(u_{1})}_{S} \\ = I \cos θ e^{i (k_{P} x_{1} \cos θ - k_{P} x_{3} \sin θ - ω t)} \\ + P \cos θ_{P} e^{i (k_{P} x_{1} \cos θ_{p} + k_{P} x_{3} \sin θ_{p} - ω t)} \\ - S \sin θ_{S} e^{i (k_{S} x_{1} \cos θ_{S} + k_{S} x_{3} \sin θ_{S} - ω t)}, \end{matrix}

(7)

\begin{matrix} u_{3} = {(u_{3})}_{I} + {(u_{3})}_{P} + {(u_{3})}_{S} \\ = - I \sin θ e^{i (k_{P} x_{1} \cos θ - k_{P} x_{3} \sin θ - ω t)} \\ + P \sin θ_{P} e^{i (k_{P} x_{1} \cos θ_{p} + k_{P} x_{3} \sin θ_{p} - ω t)} \\ + S \cos θ_{S} e^{i (k_{S} x_{1} \cos θ_{S} + k_{S} x_{3} \sin θ_{S} - ω t)} . \end{matrix}

(8)

where ${(u_{1})}_{I}$ , ${(u_{1})}_{P}$ , and ${(u_{1})}_{S}$ denote, respectively, the displacements of incident compressional, reflected compressional, and reflected shear waves along the $x_{1}$ axis; ${(u_{3})}_{I}$ , ${(u_{3})}_{P}$ , and ${(u_{3})}_{S}$ , respectively, are the displacements of incident compressional, reflected compressional, and reflected shear waves along the $x_{3}$ axis; P and S represent the amplitudes of the reflected compressional and shear waves, respectively; and $θ$ , $θ_{P}$ , and $θ_{S}$ denote, respectively, the angles of incident compressional, reflected compressional, and reflected shear wave propagation with respect to the $x_{1}$ axis, as shown in Figure 1. Here, unlike the real-valued wave numbers of plane waves in a homogeneous elastic body with no embedded particles, $k_{P}$ and $k_{S}$ represent the (complex-valued) wave numbers of compressional and shear plane waves in the particulate metacomposite with wave attenuation (see Equations (52) and (53) in Appendix 1). It is noteworthy to mention that the attenuation coefficients for compressional and shear waves, defined in Equations (57) and (59), correspond to the imaginary parts of the respective complex-valued wave numbers. Of paramount importance in this study is that, as seen from Equations (52) and (53), the ratio (k_P/k_S) is a positive constant determined by the effective elastic moduli of the particulate metacomposite.

By substituting the displacement field (Equations (7) and (8)) into the given boundary condition (Equations(2) and (3)), we obtain:

\begin{matrix} 2 k_{P} \sin θ \cos θ I e^{i (k_{P} x_{1} \cos θ - ω t)} \\ - 2 k_{P} \sin θ_{P} \cos θ_{P} P e^{i (k_{P} x_{1} \cos θ_{P} - ω t)} \\ + k_{S} (\sin^{2} θ_{S} - \cos^{2} θ_{S}) S e^{i (k_{S} x_{1} \cos θ_{S} - ω t)} = 0, \end{matrix}

(9)

\begin{matrix} - (λ \cos^{2} θ + (λ + 2 μ) \sin^{2} θ) k_{P} I e^{i (k_{P} x_{1} \cos θ - ω t)} \\ - (λ \cos^{2} θ_{P} + (λ + 2 μ) \sin^{2} θ_{P}) k_{P} P e^{i (k_{P} x_{1} \cos θ_{P} - ω t)} \\ - 2 μ k_{S} \sin θ_{S} \cos θ_{S} S e^{i (k_{S} x_{1} \cos θ_{S} - ω t)} = 0 . \end{matrix}

(10)

The exponential terms in these equations are linearly independent functions of $x_{1}$ . Consequently, the equations yield only trivial solutions for the complex amplitudes P and S unless the coefficients of $x_{1}$ are identical:

\begin{matrix} k_{P} \cos θ = k_{P} \cos θ_{P} = k_{S} \cos θ_{S}, \\ \to θ = θ_{P}, \\ \to \cos θ_{S} = \frac{k_{P}}{k_{S}} \cos θ . \end{matrix}

(11)

where $k_{P}$ and $k_{S}$ represent the (complex-valued) wave numbers of compressional and shear plane waves, respectively, in the particulate metacomposite while the ratio (k_P/k_S) is a positive constant. So, the exponential terms cancel from Equations (9) and (10), and we have:

1 - \frac{P}{I} + \frac{(\sin^{2} θ_{S} - \cos^{2} θ_{S})}{2 \sin θ \cos θ} (\frac{k_{S}}{k_{P}}) \frac{S}{I} = 0,

(12)

1 + \frac{P}{I} + \frac{2 μ \sin θ_{S} \cos θ_{S}}{(λ \cos^{2} θ + (λ + 2 μ) \sin^{2} θ)} (\frac{k_{S}}{k_{P}}) \frac{S}{I} = 0 .

(13)

Thus, the reflection angles $(θ_{P}, θ_{S})$ and the amplitudes $(P, S)$ of the reflected waves are determined by Equations (11)–(13). Here, it is emphasized that, unlike plane waves in an isotropic homogeneous elastic body with real wave numbers and zero wave attenuation, the wave numbers $(k_{P}, k_{S})$ of plane waves in a hard sphere-filled metacomposite are generally complex-valued with their non-zero imaginary parts representing wave attenuation. To the best of our knowledge, the above results (Equations (11)–(13)) remain absent in the existing literature dealing with the reflection of plane waves on the free surface of a particulate elastic composite with wave attenuation.

2.2. An incident shear wave

We now consider a plane shear wave with known frequency and amplitude incident upon the free boundary of the half-space of a particulate metacomposite in the indicated direction $θ$ . It is asserted that the incident shear wave gives rise to a reflected shear wave with the propagation direction $θ_{S}$ and a reflected compressional wave with the propagation direction $θ_{P}$ , as illustrated in Figure 1. Let the displacement of the incident shear wave be expressed as (the real part of):

u_{I} = I e^{i (k_{s} x ″ - ω t)},

(14)

where $x ″ = x_{1} \cos θ - x_{3} \sin θ$ ; $I$ and $k_{s}$ are the (complex-valued) amplitude and wave number of the incident shear wave; $ω$ , $i$ , and $t$ are, respectively, frequency, imaginary unit, and time. Since the displacement direction ( $x ″$ ) of the incident shear wave is perpendicular to its propagation direction, the displacements of the incident shear wave in the $x_{1}$ and $x_{3}$ directions are expressed as:

{(u_{1})}_{I} = I \sin θ e^{i (k_{S} x_{1} \cos θ - k_{S} x_{3} \sin θ - ω t)},

(15)

{(u_{3})}_{I} = I \cos θ e^{i (k_{S} x_{1} \cos θ - k_{S} x_{3} \sin θ - ω t)} .

(16)

Therefore, the total displacements of the incident shear wave and the reflected compressional and shear waves in the $x_{1}$ and $x_{3}$ directions are given by:

\begin{matrix} u_{1} = {(u_{1})}_{P} + {(u_{1})}_{I} + {(u_{1})}_{S} \\ = P \cos θ_{P} e^{i (k_{P} x_{1} \cos θ_{P} + k_{P} x_{3} \sin θ_{P} - ω t)} \\ + I \sin θ e^{i (k_{S} x_{1} \cos θ - k_{S} x_{3} \sin θ - ω t)} \\ - S \sin θ_{S} e^{i (k_{S} x_{1} \cos θ_{S} + k_{S} x_{3} \sin θ_{S} - ω t)}, \end{matrix}

(17)

\begin{matrix} u_{3} = {(u_{3})}_{P} + {(u_{3})}_{I} + {(u_{3})}_{S} \\ = P \sin θ_{P} e^{i (k_{P} x_{1} \cos θ_{P} + k_{P} x_{3} \sin θ_{P} - ω t)} \\ + I \cos θ e^{i (k_{S} x_{1} \cos θ - k_{S} x_{3} \sin θ - ω t)} \\ + S \cos θ_{S} e^{i (k_{S} x_{1} \cos θ_{S} + k_{S} x_{3} \sin θ_{S} - ω t)} . \end{matrix}

(18)

By substituting the displacement field into the given boundary condition (Equations (2) and (3)), we obtain:

\begin{matrix} - 2 k_{P} \sin θ_{P} \cos θ_{P} P e^{i (k_{P} x_{1} \cos θ_{P} - ω t)} \\ + k_{S} (\sin^{2} θ - \cos^{2} θ) I e^{i (k_{S} x_{1} \cos θ - ω t)} \\ + k_{S} (\sin^{2} θ_{S} - \cos^{2} θ_{S}) S e^{i (k_{S} x_{1} \cos θ_{S} - ω t)} = 0, \end{matrix}

(19)

\begin{matrix} - (λ \cos^{2} θ_{P} + (λ + 2 μ) \sin^{2} θ_{P}) k_{P} P e^{i (k_{P} x_{1} \cos θ_{P} - ω t)} \\ + 2 μ k_{S} \sin θ \cos θ I e^{i (k_{S} x_{1} \cos θ - ω t)} \\ - 2 μ k_{S} \sin θ_{S} \cos θ_{S} S e^{i (k_{S} x_{1} \cos θ_{S} - ω t)} = 0 . \end{matrix}

(20)

The exponential terms within Equations (19) and (20) are linearly independent functions of $x_{1}$ . Consequently, the equations yield only trivial solutions for the complex amplitudes P and S unless the coefficients of $x_{1}$ are identical.

\begin{matrix} k_{S} \cos θ = k_{P} \cos θ_{P} = k_{S} \cos θ_{S}, \\ \to θ = θ_{S}, \\ \to \cos θ_{P} = \frac{k_{S}}{k_{P}} \cos θ . \end{matrix}

(21)

For the incident angle $θ$ such that $\frac{k_{S}}{k_{P}} \cos θ_{S} > 1$ , we cannot solve the equation $\cos θ_{P} = \frac{k_{S}}{k_{P}} \cos θ_{S}$ for the direction of propagation of the reflected compressional wave. As the propagation direction of the incident shear wave decreases to a specific angle, the wavelength of the shear wave along the boundary aligns with the wavelength of the reflected compressional wave. This specific angle is termed the critical angle $θ_{C}$ . To study the wave reflection of a shear wave, we need to find a separate solution when $θ$ is less than the critical angle.

2.2.1. Incident angle $θ$ greater than the critical angle $θ_{C}$

When the angle $θ$ is greater than the critical angle, the analysis of the reflection of a shear wave is similar to the analysis of the reflection of a compressional wave. So, we have:

1 + \frac{S}{I} - \frac{2 \sin θ_{P} \cos θ_{P}}{\sin^{2} θ - \cos^{2} θ} (\frac{k_{P}}{k_{S}}) \frac{P}{I} = 0,

(22)

1 - \frac{S}{I} - \frac{(λ \cos^{2} θ_{P} + (λ + 2 μ) \sin^{2} θ_{P})}{2 μ \sin θ \cos θ} (\frac{k_{P}}{k_{S}}) \frac{P}{I} = 0 .

(23)

Thus, for an incident shear wave at an incident angle larger than the critical angle, the reflection angles $(θ_{P}, θ_{S})$ and the amplitudes $(P, S)$ of the reflected waves are determined by Equations (21)–(23). We again mention here that the above results (Equations (21)–(23) remain absent in the existing literature concerning the reflection of plane waves on the free surface of a particulate elastic composite with wave attenuation.

2.2.2. Incident angle $θ$ less than the critical angle $θ_{C}$

As illustrated in Figure 2, when the angle $θ$ is less than the critical angle, in addition to the reflected shear wave, there is a reflected surface wave of the form $u (x_{1}, x_{3}, t) = u_{0} e^{- h x_{3} + i (k x_{1} - wt)}$ that propagates in the $x_{1}$ direction and attenuates in the $x_{3}$ direction. Also, $k$ is the wave number and $h$ is the decay index (with positive real part) in the $x_{3}$ direction with $u_{0} = (d, 0, b)$ .

Figure 2.

The reflected shear plane wave and the surface wave of an incident shear plane wave with an incident angle less than the critical angle.

To compute the decay index $h$ of the surface wave, we utilize the modified version of the classical elastodynamic equations (Equations (46)–(47) in Appendix 1). Substituting $u (x_{1}, x_{3}, t) = u_{0} e^{- h x_{3} + i (k x_{1} - wt)}$ and $u_{mass} (x_{1}, x_{3}, t) = u_{0 mass} e^{- h x_{3} + i (k x_{1} - wt)}$ with u ₀ = ((u₁)_0mass, 0, (u₃)_0mass) into Equations (46), (47), (49), and (50), and eliminating $e^{- h x_{3}} e^{i (k x_{1} - ω t)}$ from both sides of the equations, we obtain:

- ω^{2} {(u_{1})}_{0 mass} = \frac{G}{ρ (1 - 2 ν)} (- k^{2} d - hikb) + \frac{G}{ρ} (- k^{2} d + h^{2} d),

(24)

{(u_{1})}_{0 mass} - η i ω {(u_{1})}_{0 mass} = d - η i ω d - \frac{ρ_{S}}{β} (- ω^{2} {(u_{1})}_{0 mass} + \frac{ρ_{m}}{ρ} (1 - δ) ω^{2} d),

(25)

- ω^{2} {(u_{3})}_{0 mass} = \frac{G}{ρ (1 - 2 ν)} (- ikhd + h^{2} b) + \frac{G}{ρ} (- k^{2} b + h^{2} b),

(26)

{(u_{3})}_{0 mass} - η i ω {(u_{3})}_{0 mass} = b - η i ω b - \frac{ρ_{S}}{β} (- ω^{2} {(u_{3})}_{0 mass} + \frac{ρ_{m}}{ρ} (1 - δ) ω^{2} b) .

(27)

in which $ρ$ , $ρ_{s}$ , and $ρ_{m}$ are, respectively, the mass density of the composite, stiff spheres, and matrix phase. Also, $δ$ and $β$ denote the volume fraction of embedded spheres and spring constant (per unit volume of the composite) for the composite. Furthermore, $G$ , $ν$ , and $η$ stands for shear modulus, Poisson ratio, and viscous coefficient of composite. Here, $u_{mass}$ is the displacement field of the mass center of a representative unit cell. From Equations (25) and (27), we can compute ${(u_{1})}_{0 mass} and {(u_{3})}_{0 mass}$ as follows:

{(u_{1})}_{0 mass} = \frac{1 - η i ω - \frac{ρ_{S}}{β} \frac{ρ_{m}}{ρ} (1 - δ) ω^{2}}{1 - η i ω - \frac{ρ_{S}}{β} ω^{2}} d = \frac{R + iJ}{{(1 - \frac{ρ_{S}}{β} ω^{2})}^{2} + η^{2} ω^{2}} d,

(28)

{(u_{3})}_{0 mass} = \frac{1 - η i ω - \frac{ρ_{S}}{β} \frac{ρ_{m}}{ρ} (1 - δ) ω^{2}}{1 - η i ω - \frac{ρ_{S}}{β} ω^{2}} b = \frac{R + iJ}{{(1 - \frac{ρ_{S}}{β} ω^{2})}^{2} + η^{2} ω^{2}} b .

(29)

here, R and J represent the dimensionless real numbers as defined in Equation (54). By substituting Equations (28) and (29) into Equations (24) and (26), we obtain:

- ω^{2} \frac{R + i J}{{(1 - \frac{ρ_{S}}{β} ω^{2})}^{2} + η^{2} ω^{2}} d = \frac{G}{ρ (1 - 2 ν)} (- k^{2} d + h i k b) + \frac{G}{ρ} (- k^{2} d + h^{2} d),

(30)

- ω^{2} \frac{R + iJ}{{(1 - \frac{ρ_{S}}{β} ω^{2})}^{2} + η^{2} ω^{2}} b = \frac{G}{ρ (1 - 2 ν)} (- ikhd + h^{2} b) + \frac{G}{ρ} (- k^{2} b + h^{2} b) .

(31)

The conditions for the existence of non-zero solution $(d, b)$ give a quadratic equation for $h^{2}$ , and the admissible roots of $h^{2}$ (with positive real part) are determined by:

h_{1}^{2} = k^{2} - k_{P}^{2},

(32)

h_{2}^{2} = k^{2} - k_{S}^{2} .

(33)

Therefore, with the $d - b$ relation $kb = i h_{1} d$ and $kd = - i h_{2} b$ , the general form of the displacement field for the surface wave are as follows:

u_{1} (x_{1}, x_{3}, t) = (k Q_{1} e^{- h_{1} x_{3}} + i h_{2} Q_{2} e^{- h_{2} x_{3}}) e^{i (k x_{1} - ω t)},

(34)

u_{3} (x_{1}, x_{3}, t) = (i h_{1} Q_{1} e^{- h_{1} x_{3}} - k Q_{2} e^{- h_{2} x_{3}}) e^{i (k x_{1} - ω t)} .

(35)

where $Q_{1}$ and $Q_{2}$ are two independent arbitrary complex constants (to be determined by the surface conditions). Specifically, when the incident shear wave has a small incident angle (less than the critical angle), the traction-free surface conditions (Equations (2) and (3)) are required $k = k_{S} (\cos θ)$ . Consequently, from Equations (32) and (33), we have:

\begin{matrix} h_{1}^{2} = k^{2} - k_{P}^{2} = (\cos^{2} θ - \frac{c_{2}^{2}}{c_{1}^{2}}) k_{S}^{2}, \\ \frac{c_{2}^{2}}{c_{1}^{2}} = \frac{(1 - 2 ν)}{2 (1 - ν)}, \end{matrix}

(36)

h_{2}^{2} = k^{2} - k_{S}^{2} = (\cos^{2} θ - 1) k_{S}^{2} .

(37)

where $c_{1}$ and $c_{2}$ are the effective compressional and shear wave speeds of the metacomposite half-space. Therefore, given that $(\cos^{2} θ - \frac{c_{2}^{2}}{c_{1}^{2}}) > 0$ and $(\cos^{2} θ - 1) < 0$ , the reflected surface wave has only one branch of $h_{1}$ with a positive real part. This is because the $h_{2}$ branch does not decay with $x_{3}$ . Consequently, the displacement field for the reflected surface wave with $k = k_{S} (\cos θ)$ propagating along the $x_{1}$ -direction is expressed as:

u_{1} (x_{1}, x_{3}, t) = (k Q_{1} e^{- h_{1} x_{3}}) e^{i (k_{S} \cos θ x_{1} - ω t)} = P e^{- h x_{3}} e^{i (k_{S} \cos θ x_{1} - ω t)},

(38)

u_{3} (x_{1}, x_{3}, t) = (i h_{1} Q_{1} e^{- h_{1} x_{3}}) e^{i (k_{S} \cos θ x_{1} - ω t)} = i \frac{hP}{k} e^{- h x_{3}} e^{i (k_{S} \cos θ x_{1} - ω t)} .

(39)

where $P$ is the amplitude of the surface wave. As a result, the displacement field is considered as follows:

\begin{matrix} u_{1} = {(u_{1})}_{P} + {(u_{1})}_{I} + {(u_{1})}_{S} \\ = P e^{- h x_{3}} e^{i (k_{S} \cos θ x_{1} - ω t)} \\ + I \sin θ e^{i (k_{S} x_{1} \cos θ - k_{S} x_{3} \sin θ - ω t)} \\ - S \sin θ_{S} e^{i (k_{S} x_{1} \cos θ_{S} + k_{S} x_{3} \sin θ_{S} - ω t)}, \end{matrix}

(40)

\begin{matrix} u_{3} = {(u_{3})}_{P} + {(u_{3})}_{I} + {(u_{3})}_{S} \\ = i \frac{hP}{k} e^{- h x_{3}} e^{i (k_{S} \cos θ x_{1} - ω t)} \\ + I \cos θ e^{i (k_{S} x_{1} \cos θ - k_{S} x_{3} \sin θ - ω t)} \\ + S \cos θ_{S} e^{i (k_{S} x_{1} \cos θ_{S} + k_{S} x_{3} \sin θ_{S} - ω t)} . \end{matrix}

(41)

By substituting the displacement field into the provided boundary condition equations (Equations (2) and (3)), we derive:

\begin{matrix} (- 2 ihP) e^{i (k_{S} \cos θ x_{1} - ω t)} \\ + k_{S} (\sin^{2} θ - \cos^{2} θ) I e^{i (k_{S} x_{1} \cos θ - ω t)} \\ + k_{S} (\sin^{2} θ_{S} - \cos^{2} θ_{S}) S e^{i (k_{S} x_{1} \cos θ_{S} - ω t)} = 0, \end{matrix}

(42)

\begin{matrix} (λ (\frac{h^{2} P}{k_{S} \cos θ} - k_{S} \cos θ P) + 2 μ (\frac{h^{2} P}{k_{S} \cos θ})) e^{i (k_{S} \cos θ x_{1} - ω t)} \\ + 2 μ k_{S} \sin θ \cos θ I e^{i (k_{S} x_{1} \cos θ - ω t)} \\ - 2 μ k_{S} \sin θ_{S} \cos θ_{S} S e^{i (k_{S} x_{1} \cos θ_{S} - ω t)} = 0 . \end{matrix}

(43)

These equations have only trivial solutions for the complex amplitudes P and S unless $θ = θ_{S}$ . Thus, the exponential terms cancel from Equations (42) and (43), and we have:

\frac{(- 2 i \sqrt{(\cos^{2} θ - \frac{c_{2}^{2}}{c_{1}^{2}})})}{(\sin^{2} θ - \cos^{2} θ)} \frac{P}{I} + 1 + \frac{S}{I} = 0,

(44)

\frac{- \frac{1}{\cos θ} + 2 \cos θ}{2 \sin θ \cos θ} \frac{P}{I} + 1 - \frac{S}{I} = 0 .

(45)

Thus, for an incident shear wave at an incident angle smaller than the critical angle, the amplitudes $(P, S)$ of the reflected surface compressional wave and shear plane wave are determined by Equations (44) and (45). It is again worth noting that the formulas in Equations (44) and (45) remain absent in the existing literature concerning the reflection of plane waves on the free surface of a particulate elastic composite with wave attenuation.

3. Illustrated examples

Various examples are presented in this section to evaluate the performance of the effective medium model used in the present work. The characteristic parameters of the particulate metacomposites used in this analysis are listed in Tables 1 and 2.

3.1. Validation of the model with known data on plane waves in an infinite body

In this section, we examine the validation of the present model for wave propagation in heavy hard sphere-filled random composites. To this end, since numerical results on the reflection of plane waves from the free surface of a particulate composite are not yet available in the existing literature, we compare the phase velocity and attenuation of plane waves in a whole (infinite) space predicted by Equations (56)–(59) of the present model with the known data provided by the numerical models in two commonly cited or most recent papers [23,24]. For this purpose, Figure 3(a) illustrates the phase velocity diagram of a compressional wave, normalized to the matrix’s longitudinal velocity, plotted against the normalized frequency of a metacomposite. This composite consists of an epoxy matrix (EPON 828Z) embedded with lead inclusions. The material characteristics of this metacomposite are represented in Table 1 and are taken from Sabina and Willis [23]. Also, Figure 3(b) plots the normalized attenuation of compressional wave versus the normalized frequency of the abovementioned metacomposite. In Figure 3, $c_{p}$ , $c_{pm}$ , $ω$ , $a$ , and $α$ denote the phase velocity of the compressional wave of the composite, phase velocity of the compressional wave of the surrounding matrix, frequency, radius of the lead sphere, and attenuation of the wave, respectively. Moreover, the volume fraction of lead in the composite is $5 %$ in Figure 3, and the damping coefficient of this elastic metacomposite (due to wave radiation rather than viscous effect) is determined by Equation (55) (see Ilinskii, et al. [25], Moon and Mow [26]), where $ε = 0.7448$ is a positive number on the order of unity.

Figure 3.

Plots of dimensionless compressional wave (a) phase velocity and (b) attenuation against the normalized frequency of the metacomposite as mentioned in Table 1.

Moreover, the wave propagation behavior of the shear wave for the two models is compared in Figure 4. Figure 4(a) presents the phase velocity diagram of a shear wave, normalized to the longitudinal velocity of the matrix, plotted against the normalized frequency of the abovementioned metacomposite in Table 1. Furthermore, Figure 4(b) illustrates the normalized attenuation of the shear wave versus the normalized frequency of the same material. Also, the volume fraction of lead in the composite is $5 %$ in Figure 4, and the radiation damping coefficient of this metacomposite is determined by Equation (55), where $ε = 0.7448$ is a positive number on the order of unity.

Figure 4.

Plots of dimensionless shear wave (a) phase velocity and (b) attenuation against the normalized frequency of the metacomposite as mentioned in Table 1.

According to both diagrams in Figures 3 and 4, the current wave propagation model is generally in reasonable agreement with known data reported by numerical simulation [23] for both compressional and shear waves, except the higher-frequency range in Figure 4(b) beyond the local resonance frequency peak. In view of the fact that even different numerical models often give quite different results for exactly the same examples (e.g., Kim et al. [6], Sabina and Willis [23]), the agreement between the predicted results of Equation (56)–(59) of the present model with known data is considered to be quite good.

To demonstrate the utility of the current wave propagation model, two additional examples are presented here. Figure 5(a) depicts the phase velocity diagram of a compressional wave plotted against the frequency of a metacomposite, composed of an epoxy matrix embedded with tungsten carbide inclusions. The material characteristics of this metacomposite are detailed in Table 2 and sourced from Luppé et al [24]. In addition, Figure 5(b) illustrates the attenuation of the compressional wave versus the frequency of the aforementioned metacomposite. In Figure 5, $c_{p}$ , $f$ , and $α$ represent, respectively, the phase velocity of the compressional wave of the composite, frequency, and attenuation of the wave. Furthermore, the volume fraction of lead in the composite is given as $5 %$ in Figure 5, and the radiation damping coefficient of this metacomposite is determined by Equation (55), where $ε = 0.66$ is a positive number on the order of unity.

Figure 5.

(a) phase velocity and (b) attenuation of a compressional wave of the metacomposite as mentioned in Table 2.

Furthermore, a comparison of the wave propagation characteristics of the shear wave for the two models is depicted in Figure 6. In Figure 6(a), the phase velocity diagram of a shear wave is presented and plotted against the frequency of the metacomposite as mentioned in Table 2. In addition, Figure 6(b) illustrates the attenuation of the shear wave versus the frequency of the same material. Moreover, the volume fraction of lead in the composite is represented in Figure 6; the radiation damping coefficient of this metacomposite is determined by Equations (55), where $ε = 0.6945$ is a positive number on the order of unity.

Figure 6.

(a) phase velocity and (b) attenuation of a shear wave of the metacomposite as mentioned in Table 2.

Based on the information depicted in both Figures 5 and 6, the current wave propagation model generally aligns reasonably well with the established data obtained through numerical simulations [24], for both compressional and shear waves.

3.2. Examples for reflection of an incident compressional wave

In this section, we consider two examples concerning the wave reflection of both compressional and shear incident waves at the free boundary of a half-space metacomposite. Initially, the analysis focuses on the wave reflection of incident compressional waves for two distinct types of metacomposites. Subsequently, a similar analysis is conducted for incident shear waves for two different metacomposites. Figure 7 depicts plots of the amplitude ratio of the reflected compressional and shear waves resulting from an incident compressional wave at the free surface of the metacomposite as mentioned in Table 1, as a function of $θ$ . According to Figure 7, the reflected shear wave vanishes at $θ = 0 °$ and at $θ = 90 °$ . The amplitude of the compressional and shear waves reaches its maximum at approximately $θ = 27 °$ and $θ = 41 °$ , respectively.

Figure 7.

Amplitude ratio of the reflected shear and compressional wave of the metacomposite as mentioned in Table 1 for an incident compressional wave with the incident angle of θ.

Similar to the previous example, Figure 8 presents the amplitude ratio of reflected waves as a function of for a different metacomposite as detailed in Table 2. As indicated in Figure 8, the reflected shear wave vanishes at $θ = 0 °$ and at $θ = 90 °$ . The amplitude of the compressional and shear waves peaks at around $θ = 23 °$ and $θ = 41 °$ , respectively.

Figure 8.

Amplitude ratio of the reflected shear and compressional wave of the metacomposite as mentioned in Table 2 for an incident compressional wave with the incident angle of θ.

3.3. Examples for reflection of an incident shear wave

In this section, a comparable analysis to the previous section is carried out for incident shear waves in two different metacomposites. However, as previously noted, the wave reflection analysis differs for incident angles greater or less than the critical angle. Initially, the wave reflection of two types of metacomposites is examined for incident angles greater than the critical angle. Subsequently, for the same metacomposites, a wave reflection analysis is presented for incident angles less than the critical angle. Figure 9 illustrates plots of the amplitude ratio of the reflected compressional and shear waves resulting from an incident shear wave at the free surface of the metacomposite as detailed in Table 1, as a function of $θ$ , where $θ$ is greater than the critical angle. According to Figure 9, the critical angle for this composite is approximately $θ = 63 °$ , and the amplitude of the reflected shear wave matches that of the incident shear wave.

Figure 9.

Amplitude ratio of the reflected shear and compressional wave of the metacomposite as mentioned in Table 1 for an incident shear wave with incident angle θ greater than the critical angle.

Similar to the previous example, Figure 10 depicts the amplitude ratio of reflected waves as a function of $θ$ , where $θ$ is greater than the critical angle for a different metacomposite as detailed in Table 2. According to Figure 10, the amplitude of the reflected shear wave matches that of the incident shear wave at the critical angle, which is approximately $θ = 59 °$ for the chosen composite. In addition, it is worth noting that the maximum amplitude of the reflected compressional wave occurs at the critical angle.

Figure 10.

Amplitude ratio of the reflected shear and compressional wave of the metacomposite as mentioned in Table 2 for an incident shear wave with incident angle θ greater than the critical angle.

As previously mentioned, when the angle of incidence for the shear wave is less than the critical angle, there is not only a reflected shear wave but also a reflected surface wave that propagates in the $x_{1}$ direction and attenuates in the $x_{3}$ direction. In contrast to the previous example, the amplitude ratio here is a complex number. To illustrate their variations with respect to $θ$ , Figure 11 plots both the real and imaginary parts of the amplitude ratio for both reflected waves from the free boundary of the metacomposite as mentioned in Table 1.

Figure 11.

Amplitude ratio of the reflected shear and surface wave of the metacomposite as mentioned in Table 1 for an incident shear wave with incident angle θ less than the critical angle.

Similar to the previous example, incident shear waves with an angle less than the critical angle result in both a reflected surface wave and a shear wave. Importantly, the amplitude ratio of these reflected waves is represented by complex numbers. Figure 12 illustrates both the real and imaginary parts of the amplitude ratio for these reflected waves reflecting from the free boundary of the metacomposite as detailed in Table 2. Unfortunately, a comparison with known data for the results predicted by the present model on the reflection of plane waves from the free surface of a particulate composite cannot be made here due to the lack of available data in the existing literature

Figure 12.

Amplitude ratio of the reflected shear and surface wave of the metacomposite as mentioned in Table 2 for an incident shear wave with incident angle θ less than the critical angle.

4. Summary and conclusion

In this paper, we focus on the reflection of an incident plane compressional or shear wave on the free surface of a hard sphere-filled metacomposite half-space. The main findings are as follows:

To the best of our knowledge, for the first time in the literature, the present work offers explicit formulas and numerical examples for the reflection of plane waves on the free surface of a hard sphere-filled metacomposite half-space.

The efficiency and accuracy of the present effective medium model for particulate elastic composites are demonstrated by comparing the predicted phase velocity and attenuation coefficient of plane waves in an entire (infinite) space with established known data in the literature.

Explicit formulas are presented for the reflection angles and amplitudes of the reflected waves for an incident compressional or shear plane wave for a hard sphere-filled metacomposite half-space with wave attenuation.

In addition, we provide an in-depth analysis of the reflected plane shear wave and the surface compressional wave when an incident shear wave propagates at an angle less than the critical angle.

Notably, unlike reflected plane waves in a homogeneous elastic body characterized by real wave numbers and constant amplitudes, the wave numbers of reflected plane and surface waves in a hard sphere-filled metacomposite are generally complex-valued with the amplitudes decaying with the propagating distance.

Footnotes

Appendix 1 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 a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada (grant no. RGPIN-2023-03227 Schiavo).

ORCID iDs

Chongqing Ru

Peter Schiavone

References

Achenbach

. Wave propagation in elastic solids. Amsterdam: Elsevier, 2012.

Datta

Shah

Bratton

, et al. Wave propagation in laminated composite plates. J Acoust Soc Am 1988; 83(6): 2020–2026.

Maio

Memmolo

Ricci

, et al. Ultrasonic wave propagation in composite laminates by numerical simulation. Compos Struct 2015; 121: 64–74.

Duranteau

Valier-Brasier

Conoir

J-M

, et al. Random acoustic metamaterial with a subwavelength dipolar resonance. J Acoust Soc Am 2016; 139(6): 3341–3352.

Kanaun

Levin

Propagation of longitudinal elastic waves in composites with a random set of spherical inclusions (effective field approach). Arch Appl Mech 2007; 77: 627–651.

Kim

Lee

. Dispersion of elastic waves in random particulate composites. J Acoust Soc Am 1995; 97(3): 1380–1388.

Tayari

Franklin

Tinel

, et al. Effective longitudinal wave in an elastic matrix containing coated or uncoated spheres: influence of the second order approximation including multiple scattering. J Sound Vib 2022; 541: 117325.

O’Donovan

Ibraim

O’sullivan

, et al. Micromechanics of seismic wave propagation in granular materials. Granul Matter 2016; 18: 1–18.

Tang

Yang

. Wave propagation in granular material: what is the role of particle shape. J Mech Phys Solids 2021; 157: 104605.

10.

Singh

. Reflection of plane waves at the free surface of a fibre-reinforced elastic half-space. Sadhana 2004; 29: 249–257.

11.

Graff

Pao

Y-H

. The effects of couple-stresses on the propagation and reflection of plane waves in an elastic half-space. J Sound Vib 1967; 6(2): 217–229.

12.

Dey

Addy

. Reflection of plane waves under initial stresses at a free surface. Int J Non-lin Mech 1977; 12(6): 371–381.

13.

Guha

Singh

. Effects of initial stresses on reflection phenomenon of plane waves at the free surface of a rotating piezothermoelastic fiber-reinforced composite half-space. Int J Mech Sci 2020; 181: 105766.

14.

Sidhu

Singh

. Reflection of P and SV waves at the free surface of a prestressed elastic half-space. J Acoust Soc Am 1984; 76(2): 594–598.

15.

Singh

Yadav

. Reflection of plane waves from a free surface of a rotating fibre-reinforced elastic solid half-space with magnetic field. Int J Appl Math Mech 2013; 9(9): 75–91.

16.

Deswal

Jangra

Punia

. Reflection of plane waves at the free surface of a magneto-thermoelastic medium with variable thermal conductivity and variable mass diffusivity. Waves Random Complex Media 2022; 32: 1–23.

17.

Yadav

. Reflection of plane waves from the free surface of a rotating orthotropic magneto-thermoelastic solid half-space with diffusion. J Therm Stress 2020; 44(1): 86–106.

18.

Singh

Guha

. Reflection of plane waves from the surface of a piezothermoelastic fiber-reinforced composite half-space. Mech Adv Mater Struct 2021; 28(22): 2370–2382.

19.

Selim

. Reflection of plane waves at free surface of an initially stressed dissipative medium. InTech 2009; 2: 1.

20.

Singh

Kaur

Tomar

. Waves in nonlocal elastic solid with voids. J Elast 2017; 128(1): 85–114.

21.

Biot

. Generalized boundary condition for multiple scatter in acoustic reflection. J Acoust Soc Am 1968; 44(6): 1616–1622.

22.

Wenzel

. Smoothed boundary conditions for randomly rough surfaces. J Math Phys 1974; 15(3): 317–323.

23.

Sabina

Willis

. A simple self-consistent analysis of wave propagation in particulate composites. Wave Motion 1988; 10(2): 127–142.

24.

Luppé

Conoir

J-M

Valier-Brasier

. Longitudinal and transverse coherent waves in media containing randomly distributed spheres. Wave Motion 2022; 115: 103082.

25.

Ilinskii

Meegan

Zabolotskaya

, et al. Gas bubble and solid sphere motion in elastic media in response to acoustic radiation force. J Acoust Soc Am 2005; 117(4): 2338–2346.

26.

Moon

Mow

. Wave propagation in a composite material containing dispersed rigid spherical inclusions. RAND, 1970, https://www.rand.org/pubs/research_memoranda/RM6139.html

27.

. A simple model for elastic wave propagation in hard sphere-filled random composites. J Acoust Soc Am 2022; 152(3): 1595–1604.

Reflection of plane waves from the free surface of a hard sphere-filled elastic metacomposite

Abstract

Keywords

1. Introduction

2. Wave reflection on the free surface of a metacomposite half-space

2.1. An incident compressional plane wave

2.2. An incident shear wave

2.2.1. Incident angle θ greater than the critical angle θ C

2.2.2. Incident angle θ less than the critical angle θ C

3. Illustrated examples

3.1. Validation of the model with known data on plane waves in an infinite body

3.2. Examples for reflection of an incident compressional wave

3.3. Examples for reflection of an incident shear wave

4. Summary and conclusion

Footnotes

Appendix 1

Declaration of conflicting interests

Funding

ORCID iDs

References

2.2.1. Incident angle $θ$ greater than the critical angle $θ_{C}$

2.2.2. Incident angle $θ$ less than the critical angle $θ_{C}$