Propagation of SH Waves in a Double Non-Homogeneous Crustal Layers of Finite Depth Lying Over an Homogeneous Half-Space 1

The present paper studies the Propagation of SH waves in a double non-homogeneous crustal layers lying over an isotropic homogeneous half-space, where upper layer ((i.e. rigidity and density varying trigonometrically with depth) and intermediate layer (i.e. rigidity and density varying parabolically with depth). The wave velocity equation has been obtained. Closed form solutions have been derived separately for the displacements in two non-homogeneous crustal layers and lower half-space. The dispersion curves are depicted by means of graphs for different values of non-homogeneity parameters and thickness ratio for layers.

rocket propellants, polymeric materials and growing demand of engineering and industrial applications.
The propagation of surface waves in elastic media is of considerable importance in earth-quake engineering and seismology on account of occurrence of stratification in the earth crust, as the earth is made up of different layers.As a result, the theory of surface waves has been developed by Stoneley [1924], Bullen [1965], Ewing et al. [1957], Hunters [1970] and Jeffreys [1970].
Many results of theoretical and experimental studies revealed that a real earth is considerably more complicated than the models presented earlier.This has led to a need for more realistic representation as a medium through which seismic waves propagate.The wave propagation in crystalline media plays a very interesting role in geophysics and also in ultrasonic and signal processing.Monoclinic medium is an example of such medium, keeping in the mind the fact that the nonhomogeneity characteristic is one of the most generalized elastic conditions inside the earth, many authors have studied the propagation of different waves in different media with non-homogeneity.Sezawa [1927] studied the dispersion of elastic waves propagated on curved surfaces.The transmission of elastic waves through a stratified solid medium was studied by Thomson [1950].Haskell [1953] studied the dispersion of surface waves in multilayered media.Biot [1965] studied the influence of gravity on Rayleigh waves, assuming the force of gravity to create a type of initial stress of hydrostatic nature and the medium to be incompressible.
Propagation of Love waves in a non-homogeneous stratum of finite depth sandwiched between two semi infinite isotropic media has been studied earlier by Sinha [1967].Wave propagation in a thin two-layered laminated medium with couple under initial stress was studied by Roy [1984].Datta [1986] studied the effect of gravity on Rayleigh wave propagation in a homogeneous, isotropic elastic solid medium.Effects of irregularities on the propagation of guided SH waves has been studied by Chattopadhyay et al [1983].Goda [1992] studied the effect of non-homogeneity and anisotropy on Stoneley waves.Gupta et al. [2012] investigated the influence of linearly varying density and rigidity on torsional surface waves in inhomogeneous crustal layer.
Some of the recent notable works on propagation of seismic waves in various media with different geometries were done by Chattopadhyay et al [2009Chattopadhyay et al [ , 2010Chattopadhyay et al [ , 2010]].
Recently Sethi et al. [2013] investigated the surface waves in homogeneous visco-elastic media of higher order under the influence of surface stresses.
In the present problem, we have considered the propagation of SH waves in a double nonhomogeneous crustal layers lying over an isotropic homogeneous half-space.The expression for displacement and dispersion equation is found in the closed form.The dispersion curves are depicted by means of graphs for different values of non-homogeneity parameters and thickness ratio.The influence of non-homogeneity parameters, wave number and thickness ratio of the layers on the dimensionless phase velocity has been studied.
Here, the non-homogeneity in the upper layer can be taken as Also, the non-homogeneity in intermediate layer can be tak- nz n and the lower half-space is homogeneous, where α, n are all constants and having the dimensions that are inverse of length.

FORMULATION OF THE PROBLEM
Let the density and rigidity of layers (medium-I) and (medium-II) are     and   ,   respectively and   ,   are the density and rigidity of homogeneous half-space (medium-III).
The thicknesses of the first and second layers are   and   respectively.Let us consider x-axis along the direction of the wave propagation and z-axis is taken vertically downwards.Let u, v, w are displacements along x, y, z axis respectively.Let (u1,v1,w1),(u2,v2,w2),(u3,v3,w3) are the displacements for the layers (medium-I) and (medium-II) and half-space (medium-III) respectively.
For SH wave propagation we have, The equation governing the propagation of Love wave in homogeneous isotropic elastic medium in the absence of body forces are Also, Hooke's law for isotropic medium, where λ, µ are lame's constants and Δ be the dilatation.

SOLUTION FOR UPPER LAYER (MEDIUM-I)
In view of equations ( 1), ( 3) and ( 4), equation of motion (2) gives the non-vanishing equations of motion for propagation of SH wave in upper layer (medium-II) as For a wave propagating in x-direction, equation ( 5) reduces to, where c K   be the velocity of SH waves.
in equation ( 7), to eliminate dW dz , we get We take variations in rigidity and density as Here α > 0 for non-homogeneity and α = 0 for homogeneity Introducing (9) in equation ( 8), we obtain Thus solution for equation ( 10) can be taken as where A and B are arbitrary constants.

Thus displacement component for non-homogenous layer (medium-II) is given by
Latin American Journal of Solids and Structures 13 (2016) 2628-2642

SOLUTION FOR INTERMEDIATE LAYER (MEDIUM-II)
In view of equations ( 1), ( 3) and ( 4), equation of motion (2) gives the non-vanishing equations of motion for propagation of SH wave in upper layer (medium-II) as For a wave propagating in x-direction, in equation ( 14), to eliminate dV dz , we get We take variations in rigidity and density as Here n > 0 for non-homogeneity and n= 0 for homogeneity Introducing (16) in equation ( 15), we get Thus solution for equation ( 17) can be taken as where C and D are arbitrary constants.
Thus displacement component for non-homogeneous intermediate layer (medium-II) is given by

SOLUTION FOR SEMI INFINITE HALF SPACE (MEDIUM-III)
In view of equations ( 1), ( 3) and ( 4), equation of motion (2) gives the non-vanishing equations of motion for propagation of SH wave in upper layer (medium-II) as For a wave propagating in x-direction, equation ( 20) reduces to, We take variations in rigidity and density as Introducing (23) in equation ( 22), we get Thus solution for equation ( 24) can be taken as where E and F are arbitrary constants.Thus displacement component for homogenous half-space (medium-III) is given by u x z t Ee e pz i Kx wt 3 ( , , ) .
1 2 where which gives the period equation for propagation of SH waves in a non-homogeneous isotropic two layers lying over homogeneous isotropic semi-infinite medium. where, .
which gives the wave velocity equation for propagation of SH waves in a homogeneous isotropic layer lying between two initially stressed non-homogeneous isotropic semi-infinite medium.

Case (II):
When upper layer is non-homogeneous i.e. α≠0, and intermediate layer is homogeneous, i.e. n=0, dispersion equation (32) takes the form, where Latin American Journal of Solids and Structures 13 (2016) 2628-2642 which gives the dispersion relation for propagation of SH waves in a non-homogeneous isotropic layer lying between initially stressed Gibson half-space and initially stressed non-homogeneous isotropic semi-infinite medium.Case (III): When both layers are homogeneous i.e. α=0, n=0, dispersion equation ( 32) reduces to, where, which gives the wave velocity equation for propagation of SH waves in a homogeneous isotropic layer lying over homogeneous isotropic semi-infinite medium.

Case (IV):
In the absence of upper layer i.e.H1→0, n≠0, and H2≠0, dispersion equation ( 32) which gives the wave velocity equation for propagation of SH waves in a non-homogeneous isotropic layer lying over homogeneous isotropic semi-infinite medium.

Case (V):
In the absence of upper layer i.e.H1→0, n=0, and H2≠0, dispersion equation ( 32 which gives the wave velocity equation for propagation of SH waves in a homogeneous isotropic layer lying over homogeneous isotropic semi-infinite medium.

Case (VI):
In the absence of intermediate layer i.e.H2→0, α≠0, and H1≠0, dispersion equation (32) reduces to, which gives the wave velocity equation for propagation of SH waves in a non-homogeneous isotropic upper layer lying over homogeneous isotropic semi-infinite medium.

Case(VII):
In the absence of intermediate layer i.e.H2→0, α=0, and H1≠0, dispersion equation (32) reduces to, which gives the wave velocity equation for propagation of SH waves in a homogeneous upper isotropic layer lying over homogeneous isotropic semi-infinite medium.

NUMERICAL COMPUTATIONS AND DISCUSSION
To study the effect of various dispersion non-homogeneities and thickness ratio on the propagation of SH wave propagating in double non-homogeneous layers lying over an isotropic homogeneous semi-infinite media, where upper layer (i.e.rigidity and density varying trigonometrically with depth) and intermediate layer (i.e.rigidity and density varying parabolically with depth), phase velocity is calculated numerically for equation (32), we take the following data: For non-homogeneous crustal layer (Medium-I) (Gubbins and Tierstein [1990]) For non-homogeneous crustal layer (Medium-II) (Gubbins and Tierstein [1990])    Figure 2 shows the effect of non-homogeneity parameter α/2K involved in the rigidity and density of the upper crustal layer.Following observations and effects are obtained under the above considered values.
(1a) For a particular dimensionless wave number KH and fixed value of non-homogeneity in rigidity and density of intermediate layer and thickness ratio i.e. n/K=1.0 and H2/H1=1.0, the dimensionless phase velocity (c/β1) 2 of SH waves decreases, as the value of α/2K increases from 1.0 to 2.0.(1b) For various values of α/2K and particular values of n/K and H2/H1=1.0, the phase velocity (c/β1) 2 increases as the wave number decreases in all curves 1-3.(1c) Curve 1(for α/2K=2.0) is steeper than the curve 2 (for α/2K =1.5), which in turn is steeper than curves 3 (for α/2K =1.0), that reveals that the dimensionless non-homogeneity factor α/2K has a prominent effect on SH wave propagation.
Figure 3 shows the effect of non-homogeneity parameter n/K involved in the rigidity and density of the non-homogeneous intermediate layer.Following observations and effects are obtained under the above considered values.
(2a) For a particular dimensionless wave number KH and fixed value of non-homogeneity of upper layer i.e. α/2K = 1.0, and fixed thickness ratio i.e.H2/H1=1.0, the dimensionless phase velocity (c/β1) 2 of SH waves increases, as the value of n/K increases from 0.1 to 1.0.(2b) For various values of n/K and particular values of α/2K and H2/H1, the phase velocity (c/β1) 2 increases as the wave number decreases in all curves 1-4.(2c) Curve 1(for n/K=0.0) is steeper than the curve 2 (for n/K=0.1),which in turn is steeper than curve 3, 4 (for n/K=0.6,n/K=1.0),that reveals that the dimensionless non-homogeneity factor n/K has a prominent effect on SH wave propagation.(2d) All the four curves are coinciding as wave number approaching 0.7 and dimensionless phase velocity (c/β1) 2 approaching 0.85.
Figure 4 shows the effect of thickness ratio H2/H1 on dimensionless phase velocity (c/β1) 2 .Following observations and effects are obtained under the above considered values.
(3a) For a particular dimensionless wave number KH and fixed value of non-homogeneity of upper and intermediate layers i.e. α/2K=1.0 and n/K = 1.0, the dimensionless phase velocity (c/β1) 2 of SH waves decreases, as the value of H2/H1 increases from 1.0 to 2.0.(3b) For various values of H2/H1and particular values of n/K and α/2K, the phase velocity (c/β1) 2 increases as the wave number decreases in all curves 1-3.

CONCLUSIONS
Here, we has been studied the propagation of SH waves in double non-homogeneous crustal layers lying over an isotropic homogeneous half space.Closed form solutions have been derived separately for the displacements in non-homogeneous crustal layers and lower half-space.By using the asymptotic expansion of Whittaker's function we have derived the wave velocity equation for the SH waves in compact form.Dimensionless phase velocity (c/β1) 2 is calculated numerically.The effect of various dimensionless elastic parameters, thickness ratio and non-homogeneity factors on the dimensionless phase velocity (c/β1) 2 have been shown graphically.We make the following observations 1.For a particular dimensionless wave number KH and fixed value of non-homogeneity in rigidity and density of intermediate layer and thickness ratio i.e. n/K and H2/H1, the dimensionless phase velocity (c/β1) 2 of SH waves decreases, as the value of α/2K increases.2. For various values of α/2K and particular values of n/K and H2/H1, the phase velocity (c/β1) 2 increases as the wave number decreases.3.For various values of n/K and particular values of α/2K and H2/H1, the phase velocity (c/β1) 2 increases as the wave number decreases.4. For a particular dimensionless wave number KH and fixed value of non-homogeneity of upper layer i.e. α/2K, and fixed thickness ratio i.e.H2/H1, the dimensionless phase velocity (c/β1) 2 of SH waves increases, as the value of n/K increases.5.For a particular dimensionless wave number KH and fixed value of non-homogeneity of upper and intermediate layers i.e. α/2K and n/K = 1.0, the dimensionless phase velocity (c/β1) 2 of SH waves decreases, as the value of H2/H1 increases.6.For various values of H2/H1and particular values of n/K and α/2K, the phase velocity (c/β1) 2 increases as the wave number decreases.The study of the shear wave propagation in layered media has some possible applications in the field of earthquake engineering, seismology, civil engineering and geophysics.In geophysics, they are helpful not only in the exploration of internal structure of earth which consist of many crustal layers but also in the exploration of valuable materials buried inside the earth like minerals, metals, hydrocarbons and petroleum, etc.In earthquake engineering and seismology, Shear waves are very useful in predicting earthquakes (in-terms of finding their epicentre) in dynamic response of soils and man-made structures.SH-waves cause more destruction to the structure than the body waves due to slower attenuation of the energy.So this study may be useful for understanding the cause and estimation of damage due to earthquakes.Hence, the present theoretical study may be useful for the practical application of shear waves in the non-homogeneous layered media.

Figure 1 :
Figure 1: Geometry of the problem.