Next Article in Journal
Explicit Integrating Factor Runge–Kutta Method for the Extended Fisher–Kolmogorov Equation
Next Article in Special Issue
Mathematical and Computational Modelling in Mechanics of Materials and Structures
Previous Article in Journal
Fractional Hermite–Hadamard-Type Inequalities for Differentiable Preinvex Mappings and Applications to Modified Bessel and q-Digamma Functions
Previous Article in Special Issue
Finite Element Analysis of Hierarchical Metamaterial-Based Patterns for Generating High Expansion in Skin Grafting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Consideration of Rayleigh Waves on a Coated Orthorhombic Elastic Half-Space Reinforced Using an Elastic Winkler Foundation

Department of Mathematics and Statistics, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
Math. Comput. Appl. 2023, 28(6), 109; https://doi.org/10.3390/mca28060109
Submission received: 6 June 2023 / Revised: 12 November 2023 / Accepted: 13 November 2023 / Published: 15 November 2023

Abstract

:
This article derives approximate formulations for Rayleigh waves on a coated orthorhombic elastic half-space with a prescribed vertical load acting as an elastic Winkler foundation. In addition, perfect continuity conditions are imposed between the coating layer and the substrate, while suitable decaying conditions are slated along the infinite depth of the half-space. The effect of the thin layer is modeled using appropriate effective boundary conditions within the long-wave limit. By applying the Radon transform and using the perturbation method, the derived model successfully captures the physical characteristics of elastic surface waves in coated half-spaces. The model consists of a pesudo-static elliptic equation decaying over the interior of the half-space and a singularly perturbed hyperbolic equation with a pseudo-differential operator. The pseudo-differential equation gives the approximate dispersion of surface waves on the coated half-space structure and is analyzed numerically at the end.

1. Introduction

Rayleigh surface waves are a known type of seismic wave described by Lord Rayleigh [1] that propagates along the surface of elastic media like the Earth’s crust; some of the developments recorded with regards to Rayleigh surface waves can be found in references [2,3,4] and the references appended therein. These waves are formed due to the interaction between compressional (P-waves) and shear (S-waves) waves near the surface of the Earth. In addition, when an earthquake or any other source (like volcanic activity, explosions, and even human-made sources like traffic or construction) generates seismic waves, both P-waves and S-waves are produced. However, Rayleigh waves are slower than the P-waves and S-waves, but they have longer wavelengths and are dispersive in nonhomogenous media, upon which different frequencies travel at different speeds. Furthermore, these types of waves are comprehensively studied in seismology to understand the behavior of earthquakes, evaluate the structural integrity of buildings, and aid in the exploration of subsurface geological structures [5,6,7], among other applications.
Now that the Rayleigh surface wave has been found to have a vast relevance in the exploration of subsurface geological structures, we therefore further dissect coated elastic media [8,9,10] as a particular case of these structures. In fact, coated elastic media are structures that combine elastic properties with a protective coating or layer. Elastic media, such as elastomers or polymers, are known for their ability to deform under stress and return to their original shape when the stress is removed. More so, the imposition of the additional layer as a coating to an elastic medium can serve quite a lot of functions, including aesthetic considerations, surface modification, and protection for the underlying elastic material, to mention a few [11]. In addition, coated media are very useful in our daily activities, and are found to model numerous real-life applications. For instance, in the medical industry, coating enhances the biocompatibility of elastic material, reduces friction, and further provides a sterile barrier for medical devices like hand gloves, catheters, and bandages, to state but just a few [12]. Also, the huge relevance of coated elastic structures can equally be found in the design and modeling of coated fabrics, coated elastic bands, and coated cables/wires, among others. Please refer to references [12,13,14,15,16,17,18,19,20,21,22] for more information on the application of such structures amidst the influence of external forces and excitations.
In particular, as the present study aims to examine the dynamic characteristics of the propagation of Rayleigh waves on an orthorhombic-coated orthorhombic-elastic-loaded elastic half-space, it then becomes imperative to explore a little about orthorhombic material [23]. Generally, orthorhombic is a crystallographic term that is used to describe a specific type of crystal structure exhibited by certain materials [24]. In an orthorhombic crystal structure, the lattice is defined by three mutually perpendicular axes of unequal lengths with angles of 90 degrees between each axis. Furthermore, various materials can have an orthorhombic crystal structure, including, for instance, minerals and certain metals. Some examples of orthorhombic minerals include aragonite, azurite, and topaz. These minerals admit distinctive physical properties due to their crystal structure, such as optical properties and cleavage planes. In addition, we mention the notable orthorhombic crystal, titanium dioxide [25], that is formed naturally as the mineral rutile. In essence, the orthorhombic crystal structure is one of several possible arrangements in crystalline materials aside from monoclinic, cubic, hexagonal, and tetragonal materials, to mention a few [26].
In this regard, the theory of surface waves is concerned with the development of hyperbolic–elliptic asymptotic models that capture the contribution of surface waves to the overall dynamic response when surface tractions were first prescribed by Kaplunov and Kossovich [27] and Kaplunov et al. [28]. Within these formulations, the Rayleigh wave propagation is described using a hyperbolic equation along the surface (specifically, a forced wave equation), with decay into the interior governed by quasi-static elliptic equations. They are derived by perturbing the inhomogeneous dynamic equations in linear elasticity around the eigen-solution, corresponding to surface waves of arbitrary profiles that were formerly examined by Sobolev [29], Friedlander [30], and Chadwick [21], among others, for the plane strain case, and recently extended, by Kiselev and Parker [31], to the 3D setup. In addition, the approach in reference [28] was later extended to a coated isotropic elastic half-space by Dai et al. [32]. Moreover, this extension then leads to elegant explicit approximate solutions for the near-resonant regimes of a moving load on an elastic half-space (see Erbas et al. [33]; Kaplunov et al. [34]), and for examining the significance of flexural-seismic meta-surfaces (see Wootton et al. [35]). A more methodical clarification of the approach could be found in works by Ege et al. [36], Kaplunov and Prikazchikov ([37,38]), and Mubaraki and Almalki [8], among others. Later on, the approximate model of Kaplunov et al. [28] was extended to the orthorhombic elastic half-plane by Nobili and Prikazchikov [24], and to elastic half-space of arbitrary anisotropy by Fu et al. [39].
However, the current manuscript intends to make use of the asymptotic approximation method [32] to explicitly derive approximate equations of motions and the resulting dispersion relation, governing the propagation of Rayleigh waves on an orthorhombic-coated orthorhombic-elastic-loaded elastic half-space. Furthermore, the prescribed vertically loaded excitation under consideration is taken to be induced using the Winkler elastic foundation [40], as an extension case to the known work in the literature (see reference [24] and the references therewith); equally, one may read reference [41] on the refinement of the Winkler–Fuss elastic foundation. Further, suitable perfect interfacial continuity conditions are imposed between the layer and substrate of elastic half-space, while decaying boundary conditions are presumed along the depth of the half-space. Furthermore, it is our aim to derive an approximate model with the help of the long-wave limit approximation to exhaustively capture the dynamic characteristics of surface waves on the examining structure. Indeed, the propagation of Rayleigh waves on such media is presided over using a perturbed singular hyperbolic equation with a pseudo-differential operator; such an equation shall be acquired in this study, incorporating all the physical assumptions imposed. In addition, the study shall analyze the derived model with regard to some special cases of material constants in elasticity. Furthermore, the novelty of the present work is the generalization of various considerations (see references [8,24], for instance) in the case of a 3D orthorhombic-coated orthorhombic-elastic-loaded elastic half-space and is further supported by the Winkler elastic foundation. In fact, looking at the six (6) elastic constants posed by an orthorhombic material [23] is enough to figure out the generality, or rather the complexity, of the present consideration.

2. Formulation of the Problem

Let us consider a thin orthorhombic layer of thickness h coated orthorhombic elastic half-space, which occupies the domain 0 < x j < for j = 1 , 2 and x 3 0 , further subject to a prescribed surface loading, which is reinforced using the Winkler elastic foundation; see Figure 1 for a schematic vision of the coated structure.
The 3D equations of motion are followed by (see, e.g., Achenbach [42])
σ 11 , 1 + σ 12 , 2 + σ 13 , 3 = ρ u 1 , t t , σ 21 , 1 + σ 22 , 2 + σ 23 , 3 = ρ u 2 , t t , σ 31 , 1 + σ 32 , 2 + σ 33 , 3 = ρ u 3 , t t ,
with the comma ( , ) indicating differentiation with the corresponding variables, u n = u n ( x 1 , x 2 , x 3 , t ) , and n = 1 , 2 , 3 being the plane displacements for the coating − and the half-space (substrate) + layers. Further, ρ are the mass volume densities, and σ 1 n = σ n 1 , σ 2 n = σ n 2 , σ 3 n = σ n 3 are the symmetric stress components for the orthorhombic medium, which are defined as follows:
σ 11 = c 11 u 1 , 1 + c 12 u 2 , 2 + c 13 u 3 , 3 , σ 12 = c 66 u 1 , 2 + u 2 , 1 , σ 22 = c 12 u 1 , 1 + c 22 u 2 , 2 + c 23 u 3 , 3 , σ 13 = c 55 u 1 , 3 + u 3 , 1 , σ 33 = c 13 u 1 , 1 + c 23 u 2 , 2 + c 33 u 3 , 3 , σ 23 = c 44 u 2 , 3 + u 3 , 2 ,
c 11 , c 12 , c 22 , c 13 , c 23 , c 33 , c 44 , c 55 , and c 66 are elastic constants through Voigt notation, for the coating “−” and substrate “+” layers, respectively.
Further, upon inserting the constitutive equations for the stress-displacement relation expressed in (2) into (1), one gets the following explicit equations of motions for the governing coated half-space:
c 11 u 1 , 11 + c 66 u 1 , 22 + c 55 u 1 , 33 + c 12 + c 66 u 2 , 12 + c 13 + c 55 u 3 , 13 = ρ u 1 , t t , c 66 u 2 , 11 + c 22 u 2 , 22 + c 44 u 2 , 33 + c 12 + c 66 u 1 , 12 + c 23 + c 44 u 3 , 23 = ρ u 2 , t t , c 55 u 3 , 11 + c 44 u 3 , 22 + c 33 u 3 , 33 + c 13 + c 55 u 1 , 13 + c 23 + c 44 u 2 , 23 = ρ u 3 , t t .
Additionally, the impulsive boundary conditions are prescribed on the surface of the coating x 3 = h as follows:
σ i 3 = 0 , and σ 33 = P ,
where i = 1 , 2 and P = P x 1 , x 2 , t is the prescribed vertical load, which is presumed to be induced by the Winkler elastic foundation; that is, it takes the following expression [40]:
P = a u 3 ,
where u 3 is the displacement component of the coated layer, in which the load is exerted upon, while a is the stiffness of the reinforced Winkler elastic foundation.
However, the imposed perfect continuity conditions on the interface of the two layers, that is, at x 3 = 0 , take the following expression:
σ n 3 = σ n 3 + , u n = u n + , n = 1 , 2 , 3 ,
while the decay depth-wise boundary conditions are assumed to be x 3 as follows:
u n + 0 , n = 1 , 2 , 3 .
Hence, the given equations of motions expressed in (3) for the propagation of Rayleigh waves on a coated orthorhombic elastic half-space will be asymptotically examined. Indeed, the prescribed impulsive boundary conditions on the surface of the coated layer, coupled with the imposed perfect continuity conditions, as expressed in (4)–(7), will be utilized for the acquisition of the resulting approximate solution, as well as the approximate equations of motions.

3. Derivation of the Effective Boundary Conditions

This section derives the required effective boundary conditions for the acquisition of the optimal approximate solution as well the approximate equations of motions for the governing formulation. Thus, we start off this approximation by suppressing the significance of the thin coated layer at the interface, that is, at x 3 = 0 .
Here, we take into consideration the following dimensionless small parameter [34]:
K = k h 1 ,
where k is the wavenumber. In fact, a very small wavenumber implies that the propagation of waves happens with a long-wave, while a very small frequency implies the propagation of waves is with a low-frequency (which is not our case). Please refer to reference [35] for related studies on the propagation of waves within a low-frequency long-wave band. Further, we assume the continuity conditions at x 3 = 0 to be as follows:
u n = w n + ,
with w n + = w n + x 1 , x 2 , t for n = 1 , 2 , 3 are the displacement components on the surface of the substrate “+”.
Next, we introduce the following scaled variables:
ξ i = k x i , η = x 3 h , τ = k C t ,
with
u n * = k u n , w n * = k w n + , a * = 1 k K c 44 a , σ i j * = 1 c 44 σ i j , and σ n 3 * = 1 K c 44 σ n 3 ,
where ξ i for i = 1 , 2 , and η are the scaled dimensionless spatial variables; τ is the scaled dimensionless temporal variable in the coating; a * is the scaled dimensionless stiffness of the reinforced Winkler’s foundation; u n * and w n * are scaled dimensionless displacements; σ i j * and σ n 3 * are scaled dimensionless stresses, all for n = 1 , 2 , 3 , and i j = 1 , 2 ; and C is the speed in the coating defined by
C = c 44 ρ .
Indeed, the above scaling becomes imperative in order to restrain the complete dependence of the entire structure on the coating layer. Certainly, the coating layer is partially ignored, thereby utilizing its full relevance with regards to its prescribed boundary data. In this case, only the equations of motions in the substrate remain, with an infusion of the coating boundary conditions in both the substrate equations and the resulting new scaled boundary and interfacial data. Further, the equation of motions (1) and the constitutive relations expressed in (2) can then be re-expressed in terms of these new variables, as given below:
σ i i , ξ i * + σ j i , ξ j * + σ 3 i , η * = u i , τ τ * , σ 33 , η * + K σ i 3 , ξ i * + σ j 3 , ξ j * = u 3 , τ τ * ,
and
K σ i i * = 1 c 44 c i 3 u 3 , η * + K c i i u i , ξ i * + c i j u j , ξ j * , K 2 σ 33 * = 1 c 44 c 33 u 3 , η * + K c i 3 u i , ξ i * + c j 3 u j , ξ j * , σ i j * = c 66 c 44 u i , ξ j * + u j , ξ i * , K 2 σ 31 * = c 55 c 44 u 1 , η * + ϵ u 3 , ξ 1 * , K 2 σ 32 * = u 2 , η * + ϵ u 3 , ξ 2 * .
In addition, after utilizing the scaled new variables expressed above, the prescribed boundary conditions in (6) and (9) then take the following expression:
σ i 3 * = 0 , σ 33 * = a * u 3 * at η = 1 , and u n * = w n * , at η = 0 .
It is appropriate to express the related displacement and stress components in the following expansion form
u n * w n * σ m n * = u n ( 0 ) w n ( 0 ) σ m n ( 0 ) + K u n ( 1 ) w n ( 1 ) σ m n ( 1 ) + , m , n = 1 , 2 , 3 .
Therefore, upon using the above equation in (12) and (13), the following system is obtained at the leading order:
σ i i , ξ i ( 0 ) + σ j i , ξ j ( 0 ) + σ 3 i , η ( 0 ) = u i , τ τ ( 0 ) , σ i j ( 0 ) = c 66 c 44 u i , ξ j ( 0 ) + u j , ξ i ( 0 ) , σ 33 , η ( 0 ) = u 3 , τ τ ( 0 ) , u n , η * = 0 ,
while the corresponding boundary conditions from (14) take the following form:
σ i 3 ( 0 ) = a * u 3 * , σ 33 ( 0 ) = 0 , at η = 1 , and u n ( 0 ) = w n * , at η = 0 .
The leading order displacement components satisfying (16) 4 and (17) 2 are then obtained in the following form:
u n ( 0 ) = w n * .
From (16) 3 , (17) 1 , and (18), we obtain
σ 33 ( 0 ) = η + 1 w 3 , τ τ * a * w 3 * .
Moreover, at the next order O ( ϵ ) , (13) 2 and the boundary value problem (14) 2 lead to the acquisition of
u 3 ( 1 ) = 1 c 33 c i 3 u i , ξ i ( 0 ) + c j 3 u j , ξ j ( 0 ) ,
and
u n ( 1 ) = 0 , at η = 0 .
From (20) and the boundary conditions expressed in (21), one obtains
u 3 ( 1 ) = η c 33 c i 3 w i , ξ i * + c j 3 w j , ξ j * .
By substituting (18) and (22) into (13) 1 , we obtain
σ i i ( 0 ) = 1 c 44 c i i c i 3 c 33 2 w i , ξ i * + c i j c i 3 c j 3 c 33 w j , ξ j * .
Finally, we have, from (16) 1 , (16) 2 , (18), and (23), at the same time satisfying (17) 1 , the following:
σ i 3 ( 0 ) = 1 + η c 44 c 44 w i , τ τ * c 66 w i , ξ j ξ j * c i i c i 3 c 33 2 w i , ξ i ξ i * c i i c i 3 c 33 2 w i , ξ i ξ i * c 66 + c i j c i 3 c j 3 c 33 w j , ξ i ξ j * ,
In the original variables, the stress components at the interface x 3 = 0 can then be expressed from (19) and (24) as follows:
σ i 3 + = h ρ u i , t t + c 66 u i , j j + c i i c i 3 c 33 2 u i , i i + c 66 + c i j c i 3 c j 3 c 33 u j , i j + , σ 33 + = ρ h u 3 , t t + a u 3 + .
Note that, in the absence of the effect of the Winkler elastic foundation ( a = 0 ) , the conditions in (25) may obviously be affirmed to correspond to the results reported in reference [43].

4. Application of the Perturbation Technique

Now that the related effective boundary conditions are derived in (25), we then proceed to derive the resulting pseudo-differential equation for the transverse and longitudinal potentials of the governing half-space. Indeed, the model examination of the half-space ( x 3 0 ) involves the wave dynamic equations of motions expressed in (2), and subject to the derived boundary conditions in (25).
Furthermore, we consider the special case of material constants, that is, when
c 11 = c 22 , c 13 = c 23 , and c 44 = c 55 .
The above presumption, which is referred to as a pure mode, means that the displacement component is parallel everywhere to the anti-plane motion with no additional symmetry. Subsequently, it is appropriate to deploy the Radon integral transform, defined as follows [32]:
u l r x , r , x 3 , t = u l + x cos r y sin r , x sin r + y cos r , x 3 , t d y ,
where
x = x 1 cos r + x 2 sin r , y = x 1 sin r + x 2 cos r ,
and
u x r = u 1 r cos r + u 2 r sin r , u y r = u 1 r sin r + u 2 r cos r ,
with r 0 , 1 2 π .
Now, we set u y r = 0 , which means that the anti-plane dynamic motion is dissuaded by the presence of the elastic Winkler foundation. Moreover, the equations of motions expressed in (2) for the substrate + are rewritten in terms of the present transformation as follows:
c 11 + u x , x x r + β + u 3 , x 3 r + c 55 + u x , 33 r = ρ + u x , t t r , c 33 + u 3 , 33 r + β + u x , x 3 r + c 55 + u 3 , x x r = ρ + u 3 , t t r ,
and subject to
c 55 + u x , 3 r + u 3 , x r = h ρ u x , t t r δ u x , x x r , c 13 + u x , x r + c 33 + u 3 , 3 r = ρ h u 3 , t t r a u 3 r ,
where
β + = c 13 + + c 55 + , and δ = c 11 c 13 c 33 2 .
Now, let us introduce yet another scaling of the following format:
ξ = k x c R t , γ = k x 3 τ = k K c R t .
where c R is the speed of the Rayleigh wave.
Then, the transformed equations of motions expressed in (30) can now be rewritten in the latter new scaling as follows:
c 11 + ρ + c R 2 u x , ξ ξ r + c 55 + u x , γ γ r + β + u 3 , ξ γ r = ρ + c R 2 ϵ 2 u x , τ τ r 2 ϵ u x , ξ τ r , c 33 + u 3 , γ γ r + β + u x , ξ γ r + c 55 + ρ + c R 2 u 3 , ξ ξ r = ρ + c R 2 ϵ 2 u 3 , τ τ r 2 ϵ c R u 3 , ξ τ r .
Certainly, (34) can be rewritten in the form of a single partial differential equation of the fourth-order, contacted by u x r as follows:
C 1 u x , ξ ξ ξ ξ r + C 2 u x , ξ ξ γ γ r + C 3 u x , γ γ γ γ r + K D 1 u x , ξ ξ ξ τ r + D 2 u x , ξ γ γ τ r K 2 E 1 u x , ξ ξ τ τ r + E 2 u x , γ γ τ τ r K 3 F 1 u x , ξ τ τ τ r + K 4 F 2 u x , τ τ τ τ r = 0 ,
where the coefficients C 1 , C 2 , C 3 , D i ,   E i , and F i , for i = 1 , 2 , are specified as
C 1 = c 11 + ρ + c R 2 c 55 + ρ + c R 2 ,   C 2 = c 11 + c 33 + + c 55 + 2 β + 2 c 33 + + c 55 + ρ + c R 2 ,   C 3 = c 33 + c 55 + ,   D 1 = 2 ρ + c R 2 c 11 + + c 55 + 2 ρ + c R 2 ,   D 2 = 2 ρ + c R 2 c 33 + + c 55 + ,   E 1 = ρ + c R 2 c 11 + + c 55 + 6 ρ + c R 2 ,   E 2 = ρ + c R 2 c 33 + + c 55 + ,   F 1 = 4 ρ + 2 c R 4 ,   and   F 2 = ρ + 2 c R 4 .
Moreover, the boundary conditions (31) are then reformed at γ = 0 as follows:
u x , γ r + u 3 , ξ r = c 55 c 55 + K c R 2 c 0 2 δ c 55 U x , ξ ξ r + c R 2 c 0 2 K 3 U x , τ τ r 2 K 2 U x , ξ τ r , c 13 + u x , ξ r + c 33 + u 3 , γ r = c 55 c R 2 c 0 2 K u 3 , ξ ξ r 2 K 2 u 3 , ξ τ r + K 3 u 3 , τ τ r a k u 3 r ,
where c 0 = c 55 / ρ .
Thus, accordingly, let us now expand the displacement components u x r and u 3 r as asymptotic series as follows:
u x r = K 1 U x 0 ξ , γ , τ + U x 1 ξ , γ , τ + , u 3 r = K 1 U 3 0 ξ , γ , τ + U 3 1 ξ , γ , τ + .
Then, at the leading order, (34) 1 becomes
c 11 + ρ + c R 2 U x , ξ ξ 0 + c 55 + U x , γ γ 0 + β + U 3 , ξ γ 0 = 0 ,
while (35) gives
C 1 U x , ξ ξ ξ ξ 0 + C 2 U x , ξ ξ γ γ 0 + C 3 U x , γ γ γ γ 0 = 0 .
Undeniably, the obtained elliptic equation in (39) can alternatively be represented using an operator notation as follows:
Δ 1 Δ 2 U x 0 = 0 ,
with
Δ i = ξ ξ 2 + q i γ γ 2 , i = 1 , 2 ,
where q i for i = 1 , 2 , which is determined using
q i = C 2 + 1 i C 2 2 4 C 1 C 3 2 C 1 , i = 1 , 2 ,
where C 2 2 4 C 1 C 3 0 for i = 1 , 2 . Indeed, this restriction allows the assumption of only real quantities, that is, q i 0 . Therefore, the solution for (40) can be obtained with the help of a pair of plane harmonic functions as follows:
U x 0 = ϕ 0 ξ , q 1 γ , τ + ψ 0 ξ , q 2 γ , τ ,
Then, on inserting the solution (43) into (38), amidst exploiting the application of the Cauchy–Riemann identities for the function g ξ , q γ , shown as
g , γ = q H g , ξ , g , ξ = 1 q H g , γ and H H g = g ,
where H is Hilbert transform, then we arrive at
U 3 0 = α 1 H ϕ 0 ξ , q 1 γ , τ + α 2 H ψ 0 ξ , q 2 γ , τ ,
with
α i = ρ + c R 2 c 11 + + q i 2 c 55 + β + q i , i = 1 , 2 .
Implying (43) and (45) into leading boundary conditions (36), we deduce at the surface γ = 0 the following:
α 1 q 1 ϕ , ξ 0 + α 2 q 2 ψ , ξ 0 = 0 , c 13 + + c 33 + α 1 q 1 ϕ , ξ 0 + c 13 + + c 33 + α 2 q 2 ψ , ξ 0 = 0 .
Thus, the classical Rayleigh wave equation follows:
Det α 1 q 1 α 2 q 2 c 13 + + c 33 + α 1 q 1 c 13 + + c 33 + α 2 q 2 = 0 ,
having the equivalent expression
λ = α 1 q 1 α 2 q 2 = c 13 + + c 33 + α 1 q 1 c 13 + + c 33 + α 2 q 2 .
Then, the elastic potentials ψ ( 0 ) and ϕ ( 0 ) can easily be related to each other as follows:
ψ 0 = λ ϕ 0 at γ = 0 .
Therefore, the solution obtained in (45) may be expressed in terms of only one potential function ψ 0 or ϕ 0 as follows:
U 3 0 = α 1 α 2 λ H ϕ 0 ξ , 0 , τ = 1 λ λ α 2 α 1 H ψ 0 ξ , 0 , τ .
Furthermore, upon going further to the next order, (34) 1 and (35) then take the following expressions:
c 11 + ρ + c R 2 U x , ξ ξ 1 + c 55 + U x , γ γ 1 + β + U 3 , ξ γ 1 = 2 ρ + c R 2 U x , ξ τ 0 ,
and
C 3 Δ 1 Δ 2 U x 1 = 2 ρ + c R 2 c 11 + + c 55 + 2 ρ + c R 2 U x , ξ ξ ξ τ 0 + c 33 + + c 55 + U x , ξ γ γ τ 0 .
The general solutions for U x 1 and U 3 1 are obtained in a similar manner to reference [24] as follows:
U x 1 ξ , γ , τ = ϕ 1 ξ , q 1 γ , τ + ψ 1 ξ , q 2 γ , τ + γ 2 C 3 q 2 2 q 1 2 ϑ 1 q 1 ϕ ¯ , τ ( 0 ) ϑ 2 q 2 H ψ , τ ( 0 ) ,
U 3 , γ 1 ξ , γ , τ = α 1 q 1 ϕ , ξ 1 ξ , q 1 γ , τ + α 2 q 2 ψ , ξ 1 ξ , q 2 γ , τ 1 β + 2 ρ + c R 2 + ϑ 1 c 3333 + q 2 2 q 1 2 ϕ , τ 0 1 β + 2 ρ + c R 2 + ϑ 2 c 3333 + q 1 2 q 2 2 ψ , τ 0 + γ 2 C 3 q 2 2 q 1 2 ϑ 1 α 1 H ϕ , ξ τ 0 ϑ 2 α 2 H ψ , ξ τ 0 ,              
and
U 3 , ξ 1 ξ , γ , τ = α 1 H ϕ , ξ 1 + α 2 H ψ , ξ 1 1 q 1 β + 2 ρ + c R 2 + ϑ 1 2 c 55 + q 1 α 1 β + 2 C 2 q 1 q 2 2 q 1 2 ϕ ¯ , τ 0 1 q 2 β + 2 ρ + c R 2 + ϑ 2 2 c 55 + q 2 α 2 β + 2 C 2 q 2 q 1 2 q 2 2 H ψ , τ 0 + γ 2 C 2 q 2 2 q 1 2 ϑ 2 α 2 q 2 ψ , ξ τ 0 ϑ 1 α 1 q 1 ϕ , ξ τ 0 ,
where
ϑ i = 2 ρ + c R 2 c 11 + + c 55 + c 33 + + c 55 + q i 2 2 ρ + c R 2 .
At order O ( 1 ) , the boundary conditions (36) lead to
U x , γ 1 + U 3 , ξ 1 = c 55 c 55 + c R 2 c 0 2 δ c 55 U x , ξ ξ 1 ,   c 13 + U x , ξ 1 + c 33 + U 3 , γ 1 = c 55 c R 2 c 0 2 U 3 , ξ ξ 0 a k U 3 0 , at γ = 0 .
Now, upon putting the solutions acquired in (43), (45), and (54)–(56) into (58), and via the application of the Cauchy–Riemann relations earlier expressed in (44), one obtains
α 1 q 1 H ϕ , ξ 1 + α 2 q 2 H ψ , ξ 1 α 1 q 1 δ 11 H ϕ , τ 0 + δ 12 H ψ , τ 0 c 55 c 55 + c R 2 c 0 2 δ c 1313 ϕ , ξ ξ 0 + ψ , ξ ξ 0 = 0 , c 13 + + c 33 + α 1 q 1 ϕ , ξ 1 + c 13 + + c 33 + α 2 q 2 ψ , ξ 1 c 13 + + c 33 + α 1 q 1 δ 21 ϕ , τ 0 + δ 22 ψ , τ 0 c 55 c R 2 c 0 2 α 1 H ϕ , ξ ξ 0 + α 2 H ψ , ξ ξ 0 = a k α 1 α 2 λ H ϕ 0 , at γ = 0 ,
where
α 1 q 1 δ 1 i = 1 q i β + 2 ρ + c R 2 + ϑ j α j β + + β + 2 c 55 + q i 2 C 2 q i q i 2 q j 2 , and c 13 + + c 33 + α 1 q 1 δ 2 i = 1 β + 2 ρ + c 33 + c R 2 ϑ i q i 2 q j 2 .
Using (50) and implicit differentiation with respect to ξ , we arrive at
c 13 + + c 33 + α 2 q 2 c 13 + + c 33 + α 1 q 1 α 2 q 2 α 1 q 1 ψ , ξ ξ 1 δ 21 δ 11 λ δ 22 δ 12 ϕ , ξ τ 0 ρ c R 2 δ 1 λ c 55 + α 1 q 1 + ρ c R 2 α 1 λ α 2 c 13 + + c 33 + α 1 q 1 H ϕ , ξ ξ ξ 0 = a α 1 α 2 λ H ϕ , ξ 0 k c 13 + + c 33 + α 1 q 1 , at γ = 0 .
By simplifying this formula, we obtain
2 ϕ , ξ τ 0 + b q 1 ϕ , ξ ξ γ 0 = a α 1 α 2 λ ϕ , γ 0 k q 1 c 13 + + c 33 + α 1 q 1 , at γ = 0 ,
where the constant b inherits the properties of both coating and substrate, given explicitly by
b = 1 B ρ c R 2 δ 1 λ c 55 + α 1 q 1 + ρ c R 2 α 1 λ α 2 c 13 + + c 33 + α 1 q 1 ,
which takes both positive and negative values corresponding to the local minimum and maximum of the phase velocity equal to the Rayleigh wave speed, while the constant B contains properties of the substrate only takes the following form:
B = 1 2 δ 21 δ 11 λ δ 22 δ 12 .
Then, on re-expressing (62) in the form of the original variables x , x 3 , t , one obtains
ϕ , x x r 1 c R 2 ϕ , t t r + b h q 1 ϕ , x x 3 r = a α 1 α 2 λ ϕ , 3 0 B q 1 c 13 + + c 33 + α 1 q 1 , at x 3 = 0 .
In addition, the elliptic equation for the potentials ϕ r and ψ r are then found to be
ϕ , x x r + q 1 ϕ , 33 r = 0 , and ψ , x x r + q 2 ψ , 33 r = 0 .
Moreover, let us now introduce their respective inverse transforms as follows:
ψ 1 r = ψ r cos r , and ψ 2 r = ψ r sin r .
Hence, (65) and (66) may be reinterpreted by inverting the respective transforms as follows:
Δ ϕ + q 1 ϕ , 33 = 0 , and Δ ψ i + q 2 ψ i , 33 = 0 ,
which govern the decay along the depth of the half-space x 3 0 , the interior, and in addition to the boundary conditions at x 3 = 0 given below:
Δ ϕ 1 c R 2 ϕ , t t + b h q 1 Δ ϕ , 3 = a α 1 α 2 λ ϕ , 3 B q 1 c 13 + + c 33 + α 1 q 1 ,
and
Δ ψ i 1 c R 2 ψ i , t t + b h q 2 Δ ψ i , 3 = a λ α 2 α 1 ψ i , 3 λ B q 2 c 13 + + c 33 + α 1 q 1 .
Once again, let us say that Equation (68) can be presented in terms of a pseudo-differential equation on the surface x 3 = 0 of the coated structure as
Δ ϕ 1 c R 2 ϕ , t t b h Δ Δ ϕ = a α 1 α 2 λ Δ ϕ B c 13 + + c 33 + α 1 q 1 ,
where Δ is the pseudo-differential operator (for more details, see reference [32]).
Therefore, (70) reduces to the plane strain problem x 1 , x 3 , t ; that is, it takes the following form:
ϕ , 11 1 c R 2 ϕ , t t b h 11 2 ϕ , 11 = a α 1 α 2 λ 11 2 ϕ B c 13 + + c 33 + α 1 q 1 at x 3 = 0 .
Certainly, Equation (71) leads to the acquisition of the resulting approximate dispersion by using the solution ϕ ( x 1 , 0 , t ) = f ( 0 ) e i k ( x 1 c t ) , where c is the phase velocity, and f ( 0 ) is an arbitrary function, which further results in obtaining
c c R = 1 b K + a h α 2 λ α 1 K B c 13 + + c 33 + α 1 q 1 ,
with c R in the latter equation equally denoting the speed of the Rayleigh wave.

5. Model Verification

In the case of the absence of the coating layer and the Winkler elastic foundation, that is, when h = 0 and a = 0 , then (71) may be identical to the hyperbolic Equation (38) in reference [24]. Moreover, for an isotropic material case, the specified elastic constants c m n m = n = 1 , 2 , 3 for the orthorhombic elastic now take the following reduced form:
c 11 = c 33 = λ + 2 μ , c 13 = λ , c 55 = 2 μ ,
where λ and μ are the respective Lame’s elastic constants [42], then the equation expressed in (70) may be compared with Equation (5.3) of Dai et al. [32] while disregarding the effect of external loading as follows:
Δ ϕ 1 c R 2 ϕ , t t b 0 h Δ Δ ϕ = 0 ,
where
b 0 = μ 1 β R 2 2 B 0 ρ c R 2 μ α R + β R 4 β R 1 κ 0 2 ,
with
B 0 = β R α R 1 α R 2 + α R β R 1 β R 2 1 β R 4 , α R = 1 ρ + c R 2 λ + + 2 μ + , β R = 1 ρ + c R 2 μ + ,
and
μ = μ μ + κ 0 = λ + 2 μ μ .
Furthermore, the approximate dispersion relation expressed in (72) may then be reduced, leading to a similar approximate dispersion relation, as reported in reference [8], given as
c c R = 1 b 0 K b 1 ζ K ,
where
b 1 = μ α R 1 β R 2 2 B 0 , and ζ = h a μ .
In short, it is part of the novelty of the present work to categorically state that various considerations have been generalized by our examination. As an example, in the absence of the reinforcement induced by the elastic Winkler foundation, the results of the present study are reduced to those of Nobili and Prikazchikov, as reported in reference [24]. Further, when the isotropic material case is considered, the present study matches the results obtained by Dai et al. [32]. In the same vein, our finding coincides with the recent results of Mubaraki and Almalki [8] in the absence of the action of magnetic field force.
In this regard, we endeavour here to numerically analyze the significance of the imposed loading on the coated substrate, which was presided over by the elastic Winkler foundation. In light of this, the obtained approximate dispersion relation in (72) for orthorhombic-coated orthorhombic half-space, is simulated numerically, considering the combination of soft-stiff and stiff-soft materials, respectively. In fact, the cilicon (Si) material [44] is sought after as a soft material, which admits the following material data:
c 11 = 11.6 GPa , c 13 = 5.4 GPa , c 33 = 16.6 GPa , c 55 = 9.5 GPa , ρ = 2329 kgm 3 ,
while for the stiff material, aluminum nitride (AIN) material [45] is considered, which has the following physical data:
c 11 + = 34.5 GPa , c 13 + = 12.0 GPa , c 33 + = 39.5 GPa , c 55 + = 11.8 GPa , ρ + = 3260 kgm 3 .
Thus, we have portrayed in Figure 2 and Figure 3 the influence of the scaled dimensionless Winkler foundation parameter ζ + = h a c 55 + on scaled dimensionless phase speed c c R against the dimensionless wavenumber K. More precisely, Figure 2 and Figure 3 shows the variational significance of the Winkler foundation parameter on the dispersion of waves in a soft-coated stiff-substrate, and in stiff-coated soft-substrate structures, respectively. Notably, from Figure 2, it is noted that an increase in the Winkler foundation parameter lessens the acquired approximate dispersion relation through the phase speed versus the wavenumber curve. Moreover, one can also observe that the dispersion relation attains its maximum in the absence of the elastic Winkler foundation. Additionally, when a swap of materials is made between the material constants of the coating and that of the half-space substrate layers, as portrayed in Figure 3, that is, a media with stiff-coated soft-substrate, an opposite tread is realized. Thus, the choice of material or combination of materials is very important with regard to the vibration analysis and design of dissimilar single, coated, and multilayered structures.

6. Conclusions

The present study asymptotically derived the approximate equations of motions and dispersion relation governing the propagation of Rayleigh waves on a loaded orthorhombic-coated orthorhombic elastic half-space. More precisely, the prescribed vertically loaded excitation was presumed to be in favor of an elastic Winkler foundation. Indeed, perfect continuity conditions were imposed between the coated layer and elastic half-space. Certainly, the derived model was found to comprehensively capture the physical characteristics of elastic surface waves, where the propagation of Rayleigh waves on the governing media was described using a singularly perturbed hyperbolic equation, admitting a pseudo-differential operator. Furthermore, upon utilizing the long-wave limit approximation for elastic surface waves, the decay over the interior of the half-space was described using a pseudo-static elliptic equation, through the acquisition of appropriate effective boundary conditions; further, the significance of the thin coating layer on the dispersion of surface waves on the coated structure was equally examined. Finally, as the orthorhombic material happened to generalize several other materials of real-life relevance, the present study then serves as an interesting monograph for the examination of the dispersion of surface waves on coated media in the fields of linear elasticity and material science.

Funding

This research received no external funding.

Acknowledgments

The researcher would like to acknowledge Deanship of Scientific Research, Taif University for funding this work.

Conflicts of Interest

The author declares that he has no conflict of interest.

References

  1. Rayleigh, L. On waves propagated along the plane surface of an elastic solid. Proc. Lond. Math. Soc. 1885, 1, 4–11. [Google Scholar] [CrossRef]
  2. Barnett, D.M.; Lothe, J. Consideration of the existence of surface wave (Rayleigh wave) solutions in anisotropic elastic crystals. J. Phys. F Metal Phys. 1974, 4, 671. [Google Scholar] [CrossRef]
  3. Fu, Y.B.; Mielke, A. A new identity for the surface-impedance matrix and its application to the determination of surface-wave speeds. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 2002, 458, 2523–2543. [Google Scholar] [CrossRef]
  4. Destrade, M. Seismic Rayleigh waves on an exponentially graded, orthotropic half-space. Proc. R. Soc. A Math. Phys. Eng. Sci. 2007, 463, 495–502. [Google Scholar] [CrossRef]
  5. Palermo, A.; Krodel, S.; Marzani, A.; Daraio, C. Engineered metabarrier as shield from seismic surface waves. Sci. Rep. 2016, 6, 39356. [Google Scholar] [CrossRef] [PubMed]
  6. Cho, Y.S. Non-destructive testing of high strength concrete using spectral analysis of surface waves. NDT E Int. 2003, 36, 229–235. [Google Scholar] [CrossRef]
  7. Krylov, V.V. (Ed.) Noise and Vibration from High-Speed Trains; Thomas Telford: London, UK, 2001. [Google Scholar]
  8. Mubaraki, A.M.; Almalki, F.M. Surface waves on a coated homogeneous half-space under the effects of external forces. Symmetry 2022, 14, 2241. [Google Scholar] [CrossRef]
  9. Mubaraki, A.; Prikazchikov, D.; Kudaibergenov, A. Explicit model for surface waves on an elastic half-space coated by a thin vertically inhomogeneous layer. In DSTA 2019: Perspectives in Dynamical Systems I: Mechatronics and Life Sciences; Springer: Cham, Switzerland, 2019; pp. 267–275. [Google Scholar] [CrossRef]
  10. Althobaiti, S.; Mubaraki, A.; Nuruddeen, R.I.; Gomez-Aguilar, J.F. Wave propagation in an elastic coaxial hollow cylinder when exposed to thermal heating and external load. Results Phy. 2022, 38, 105582. [Google Scholar] [CrossRef]
  11. Tiainen, V.M. Amorphous carbon as a bio-mechanical coating-mechanical properties and biological applications. Diam. Relat. Mater. 2001, 10, 153–160. [Google Scholar] [CrossRef]
  12. Asif, M.; Nuruddeen, R.I.; Nawaz, R. Propagation of elastic waves in a magneto-elastic layer laying over a light Winkler foundation with rotation. Waves Random Complex Media 2023. [Google Scholar] [CrossRef]
  13. Knopoff, L. The interaction between elastic wave motions and magnetic field in electrical conductors. J. Geophys. Res. 1955, 60, 441–456. [Google Scholar] [CrossRef]
  14. Chadwick, P. Elastic wave propagation in a magnetic field. Aces IX Congr. Int. Mech. Appl. 1957, 7, 143–158. [Google Scholar]
  15. Kaliski, S.; Petykiewicz, J. Equation of motion coupled with the field of temperature in a magnetic field involving mechanical and electrical relaxation for anisotropic bodies. Proc. Vib. Probl. 1959, 4, 1. [Google Scholar]
  16. Abubakar, I. Magneto-elastic SH-type of motion. Pure Appl. Geophys. 1964, 59, 10–20. [Google Scholar] [CrossRef]
  17. Mubaraki, A.M.; Nuruddeen, R.I.; Gomez-Aguilar, J.F. Modelling the dispersion of waves on a loaded bi-elastic cylindrical tube with variable material constituents. Results Phys. 2023, 53, 106927. [Google Scholar] [CrossRef]
  18. Nawaz, R.; Asif, M.; Nuruddeen, R.I.; Alahmadi, H. Phase velocity analysis by multi-parametric variations in a highly heterogeneous sandwich plate structure embedded in the Pasternak foundations with viscoelastic interlayer. Mech. Based Des. Struct. Mach. 2023. [Google Scholar] [CrossRef]
  19. Nath, S.; Sengupta, P.R. Influence of gravity on propagation of waves in a medium in presence of a compressional source. Sadhana 1999, 24, 495–505. [Google Scholar] [CrossRef]
  20. Nuruddeen, R.I.; Nawaz, R.; Zaigham Zia, Q.M. Effects of thermal stress, magnetic field and rotation on the dispersion of elastic waves in an inhomogeneous five-layered plate with alternating components. Sci. Prog. 2020, 103, 0036850420940469. [Google Scholar] [CrossRef]
  21. Chadwick, P. Surface and interfacial waves of arbitrary form in isotropic elastic media. J. Elast. 1976, 6, 73–80. [Google Scholar] [CrossRef]
  22. Abo-Dahab, S.M.; Abd-Alla, A.M.; Khan, A. Rotational effect on Rayleigh, Love and Stoneley waves in non-homogeneous fibre-reinforced anisotropic general viscoelastic media of higher order. Struct. Eng. Mech. 2016, 58, 181–197. [Google Scholar] [CrossRef]
  23. Chadwick, P. The existence of pure surface modes in elastic materials with orthorhombic symmetry. J. Sound Vibr. 1976, 47, 39–52. [Google Scholar] [CrossRef]
  24. Nobili, A.; Prikazchikov, D.A. Explicit formulation for the Rayleigh wave field induced by surface stresses in an orthorhombic half-plane. Eur. J. Mech.-A 2018, 70, 86–94. [Google Scholar] [CrossRef]
  25. Siddiqui, H. Modification of Physical and Chemical Properties of Titanium Dioxide (TiO2) by Ion Implantation for Dye Sensitized Solar Cells. In Ion Beam Techniques and Applications; IntechOpen: Rijeka, Croatia, 2019. [Google Scholar] [CrossRef]
  26. Marzouki, R. Introductory Chapter: Crystalline Materials and Applications. In Synthesis Methods and Crystallization; IntechOpen: Rijeka, Croatia, 2020. [Google Scholar] [CrossRef]
  27. Kaplunov, Y.D.; Kossovich, L.Y. Asymptotic model of Rayleigh waves in the far-field zone in an elastic half-plane. In Doklady Physics; Nauka/Interperiodica: Noida, India, 2004; Volume 49, pp. 234–236. [Google Scholar]
  28. Kaplunov, J.; Zakharov, A.; Prikazchikov, D. Explicit models for elastic and piezoelastic surface waves. IMA J. Appl. Math. 2006, 71, 768–782. [Google Scholar] [CrossRef]
  29. Sobolev, S.L.; Frank, P.; von Mises, R. Some problems in wave propagation. In Differential and Integral Equations of Mathematical Physics; ONTI: Moscow, Russia, 1937; pp. 468–617. [Google Scholar]
  30. Friedlander, F.G. On the total reflection of plane waves. Q. J. Mech. Appl. Math. 1948, 1, 376–384. [Google Scholar] [CrossRef]
  31. Kiselev, A.P.; Parker, D.F. Omni-directional Rayleigh, Stoneley and Schölte waves with general time dependence. Proc. R. Soc. Math. Phys. Eng. 2010, 466, 2241–2258. [Google Scholar] [CrossRef]
  32. Dai, H.-H.; Kaplunov, J.; Prikachikov, D.A. long-wave model for the surface wave in a coated half-space. Proc. R. Soc. A 2010, 466, 3097–3116. [Google Scholar] [CrossRef]
  33. Erbaş, B.; Kaplunov, J.; Prikazchikov, D.A.; Şahin, O. The near-resonant regimes of a moving load in a three-dimensional problem for a coated elastic half-space. Math. Mech. Solids 2017, 22, 89–100. [Google Scholar] [CrossRef]
  34. Kaplunov, J.; Prikazchikov, D.A.; Erbaş, B.; Şahin, O. On a 3D moving load problem for an elastic half space. Wave Motion 2013, 50, 1229–1238. [Google Scholar] [CrossRef]
  35. Wootton, P.T.; Kaplunov, J.; Colquitt, D.J. An asymptotic hyperbolic–elliptic model for flexural-seismic metasurfaces. Proc. R. Soc. A 2019, 475, 20190079. [Google Scholar] [CrossRef]
  36. Ege, N.; Erbaş, B.; Prikazchikov, D.A. On the 3D Rayleigh wave field on an elastic half-space subject to tangential surface loads. ZAMM J. Appl. Math. Mech./Z. Angew. Math. Mech. 2015, 95, 1558–1565. [Google Scholar] [CrossRef]
  37. Kaplunov, J.; Prikazchikov, D.A. Explicit models for surface, interfacial and edge waves. In Dynamic Localization Phenomena in Elasticity, Acoustics and Electromagnetism; Springer: Berlin/Heidelberg, Germany, 2013; pp. 73–114. [Google Scholar]
  38. Kaplunov, J.; Prikazchikov, D.A. Asymptotic theory for Rayleigh and Rayleigh-type waves. Adv. Appl. Mech. 2017, 50, 1–106. [Google Scholar]
  39. Fu, Y.; Kaplunov, J.; Prikazchikov, D. Reduced model for the surface dynamics of a generally anisotropic elastic half-space. Proc. R. Soc. A 2020, 476, 20190590. [Google Scholar] [CrossRef]
  40. Erbas, B.; Kaplunov, J.; Nobili, A.; Kilic, G. Dispersion of elastic waves in a layer interacting with a Winkler foundation. J. Acoust. Soc. Am. 2018, 144, 2918–2925. [Google Scholar] [CrossRef] [PubMed]
  41. Kaplunov, J.; Prikazchikov, D.; Sultanova, L. Justification and refinement of Winkler-Fuss hypothesis. Z. Math. Phys. 2018, 69, 1–15. [Google Scholar] [CrossRef]
  42. Achenbach, J.D. Wave Propagation in Elastic Solids, Eight Impression; Elsevier: Amsterdam, The Netherlands, 1999. [Google Scholar]
  43. Vinh, P.C.; Linh, N.T.K. An approximate secular equation of Rayleigh waves propagating in an orthotropic elastic half-space coated by a thin orthotropic elastic layer. Wave Motion 2012, 49, 81–689. [Google Scholar] [CrossRef]
  44. Sotnikov, A.V.; Schmidt, H.; Weihnacht, M.; Smirnova, E.P.; Chemekova, T.Y.; Makarov, Y.N. Elastic and piezoelectric properties of AlN and LiAlO2 single crystals. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2010, 57, 808–811. [Google Scholar] [CrossRef] [PubMed]
  45. Tsubouchi, K.; Sugai, K.; Mikoshiba, N. AlN material constants evaluation and SAW properties on AlN/Al2O3 and AlN/Si. In Proceedings of the 1981 Ultrasonics Symposium, Winnipeg, MB, Canada, 10–12 September 1981; pp. 375–380. [Google Scholar]
Figure 1. A coated elastic orthorhombic half-space reinforced using the Winkler elastic foundation.
Figure 1. A coated elastic orthorhombic half-space reinforced using the Winkler elastic foundation.
Mca 28 00109 g001
Figure 2. Influence of the scaled dimensionless Winkler foundation parameter ζ + on scaled dimensionless phase speed c c R versus the dimensionless wave number K on a Si-coated AIN-substrate.
Figure 2. Influence of the scaled dimensionless Winkler foundation parameter ζ + on scaled dimensionless phase speed c c R versus the dimensionless wave number K on a Si-coated AIN-substrate.
Mca 28 00109 g002
Figure 3. Influence of the scaled dimensionless Winkler foundation parameter ζ + on scaled dimensionless phase speed c c R versus the dimensionless wavenumber K on an AIN-coated Si-substrate.
Figure 3. Influence of the scaled dimensionless Winkler foundation parameter ζ + on scaled dimensionless phase speed c c R versus the dimensionless wavenumber K on an AIN-coated Si-substrate.
Mca 28 00109 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mubaraki, A.M. Asymptotic Consideration of Rayleigh Waves on a Coated Orthorhombic Elastic Half-Space Reinforced Using an Elastic Winkler Foundation. Math. Comput. Appl. 2023, 28, 109. https://doi.org/10.3390/mca28060109

AMA Style

Mubaraki AM. Asymptotic Consideration of Rayleigh Waves on a Coated Orthorhombic Elastic Half-Space Reinforced Using an Elastic Winkler Foundation. Mathematical and Computational Applications. 2023; 28(6):109. https://doi.org/10.3390/mca28060109

Chicago/Turabian Style

Mubaraki, Ali M. 2023. "Asymptotic Consideration of Rayleigh Waves on a Coated Orthorhombic Elastic Half-Space Reinforced Using an Elastic Winkler Foundation" Mathematical and Computational Applications 28, no. 6: 109. https://doi.org/10.3390/mca28060109

APA Style

Mubaraki, A. M. (2023). Asymptotic Consideration of Rayleigh Waves on a Coated Orthorhombic Elastic Half-Space Reinforced Using an Elastic Winkler Foundation. Mathematical and Computational Applications, 28(6), 109. https://doi.org/10.3390/mca28060109

Article Metrics

Back to TopTop