Dispersion of Love wave in an isotropic layer sandwiched between orthotropic and prestressed inhomogeneous half-spaces

An in-depth study has been carried out for the dispersion of Love waves in an isotropic elastic layer sandwiched between orthotropic and prestressed inhomogeneous elastic half-spaces. The inhomogeneities in density and rigidity of the lower half-space are space dependent and an arbitrary function of depth. Simple mathematical techniques are used to obtain dispersion relation for Love wave propagation in an isotropic layer. An extensive analysis is carried out through numerical computation to explore the effect of inhomogeneity and initial stress the lower half on the phase velocity of the Love waves. The numerical analysis of dispersion equation manifests that the phase velocity of the Love wave increases with the increase of stress parameter. The results further indicate that the inhomogeneity of the half space affect the wave velocity significantly. These results can be useful to study geophysical prospecting and understanding the cause and estimation of damage due to earthquakes.


INTRODUCTION
The deformation at any point of the medium is useful to analyze the deformation field around mining tremors and drilling into the crust of the earth.It may also find application in various engineering problems, crystal physics and solid-earth geophysics regarding deformation of an anisotropic solid.In fact, study of surface waves in non-homogeneous and layered media has been of central interest to theoretical and experimental seismologists.Our Earth is a spherical and layered solid under high initial stress.Due to variation of temperature, gravitating pull, atmosphere, slow process of creep and pressure due to crustal layer, the critical initial stresses are stored in the layer of the Earth.At the present time the usefulness of dislocation theory in seismology is restricted by the absence of detailed knowledge of either the tectonic stress which drives the system or the stress which resists slip on the fault plane and by the absence of detailed observations of deformation preceding, accompanying of dislocation theory to seismology lie in the mathematical theory but rather in the basis mechanics of faulting.The stresses which exist in an elastic body even though the external forces are absent are termed as prestresses.These stresses might exert significant effect on the elastic waves produced by earthquakes.The propagation of Love waves in a nonhomogeneous elastic media is of considerable importance in earth-quake engineering and seismology on account of occurrence of non-homogeneities in the earth crust, as the earth is made up of different layers.The mathematical expression provides the bridge between modelling results and field applications.
Surface waves are very important in the study of earthquake engineering, geophysics and geodynamics.Love waves cause more destruction to the structure than that of the body waves due to its slower attenuation of the energy.The supplement of surface wave analysis and other wave propagation problems to anisotropic elastic materials has been the subject of many studies.Many authors have discussed Love wave propagation by considering various irregularities, inhomogeneities and boundaries of the Earth.Love (1944) and Ewing et al. (1957) proposed the propagation of waves in transversely isotropic medium.Chatopadhyay (1975) discussed Love waves due to irregularity in the thickness of the non-homogeneous crystal layer.Deresiewicz (1962) studied the propagation of Love waves in a homogeneous crust overlying an inhomogeneous substratum.Bhattacharya (1969) examined the Love waves in intermediate heterogeneous layer placed between isotropic elastic half-spaces.Midya (2004) discussed Love waves in micropolar homogeneous elastic media.Manna et al. (2013) discussed propagation of Love wave in hetrogeneous elastic half-space and piezoelectric layer.Du et al. (2008) studied the effect of initial stress on the propagation of piezoelectric layered structures loaded with viscous liquid.Liu and Wang (2005) studied Love waves in functionally graded layered piezoelectric structure.Chakraborty and Dey (1982) discussed the propagation of Love waves in water saturated soil underlain by heterogeneous elastic medium.Ke et al. (2005) discussed Love waves in nonhomogeneous fluid saturated porous layered half-space.Kundu et al. (2013) discussed propagation of Love wave in porous rigid layer kept over prestressed half space.Chattaraj et al. (2013) discussed Love wave propagation in irregular prestressed anisotropic porous stratum.Ghorai et al. (2010) showed the effect of rigid boundary on the propagation of Love wave in porous layer placed over an elastic half-space.Kadian and Singh (2010) studied the influence of size of barrier on Love wave reflection.Ahmed and Abd-Dahab (2010) studied the effect of initial stress on Love waves in an orthotropic Granular layer.Gupta et al. (2013) proposed a mathematical model to study Love wave propagation in homogeneous and initially stressed hetrogeneous half-spaces.Presently, Madan et al. (2014) investigated propagation of Love waves in saturated porous anisotropic layer.Kakar and Gupta (2014) studied Love waves in an intermediate heterogeneous layer lying in between homogeneous and inhomogeneous isotropic elastic half-spaces.More recently, Kundu et al. (2014a;2014b) have examined Love wave propagation in fiber-reinforced media.
The present paper deals with the study of propagation of Love wave in a sandwiched layer lying between orthotropic and inhomogeneous half spaces.Five different cases have been studied for propagation of Love waves in a layer.The dispersion equations of Love waves under assumed conditions have been derived.Also numerical computation of dispersion equation has been performed to show the effect of initial stresses and inhomogeneity parameters on the propagation of Love waves.It has been found that initial stress parameter, rigidity parameter and density parameter of the lower half-space affect the phase velocity of Love waves.

FORMULATION OF THE PROBLEM
We have considered an isotropic and homogeneous layer of thickness H (denoted as 2 M ) sandwiched between two orthotropic and prestressed inhomogeneous (denoted as 1 M and 3 M ) halfspaces (as shown in Fig. 1).Let 2  and 2  be the rigidity and density of the intermediate layer and rigidity and density in the upper half-space are 1  and 1  .The origin has been taken at the lower interface, Love wave propagates toward x -axis, while the positive z -axis toward the interior of the lower half space.The rigidity and density of the lower half are space dependent and an arbitrary function of depth i.e.


 . Here  is the inhomogeneous parameter of lower half-space and having dimension that are inverse of length.Here  is inhomogeneous parameter of lower half-space and having dimension that are inverse of length.The upper portion of superficial half-space corresponds to the free surface with zero relative density and rigidity as

SOLUTION OF THE PROBLEM
3.1 Solution for the upper half-space Equation of motion for upper half-space in the absence of body forces can be written as (Love, 1911)  e and zz e are defined by Using Love wave conditions 1) and ( 2), the equation of motion for the upper orthotropic half-space becomes and stress-strain relations reduces to To solve Eq. ( 4) we take the following substitution and k is wave number.Using Eq. ( 6) in Eq. ( 4), we get Latin American Journal of Solids and Structures 12 (2015) 1934-1949 where Therefore, the solution for the upper orthotropic half-space is given by 1 A e exp ( ) where  is arbitrary constant.

Solution for the lower half-space
Equation of motion for lower half-space under initial stress P acting along x -axis can be written as (Love, 1911) w are the components of the displacement vector and 3  is the density of the lower half-space.Here, x  , y  and z  are the rotational components in the lower half-space, which are defined by , Eq. ( 10) can be reduced to Latin American Journal of Solids and Structures 12 (2015) 1934-1949 The stress-strain relations are The inhomogeneity of rigidity and density of the lower half-space are Now, substituting the inhomogeneity of rigidity from Eq. ( 14) in Eq. ( 13), we have The equation of motion ( 12) with the help of equations ( 14) and ( 15) can be written as To solve Eq. ( 16) we take the following substitution Using Eq. ( 17) in Eq. ( 16), we get where c is phase velocity and Latin American Journal of Solids and Structures 12 (2015) 1934-1949 Introducing the non-dimensional quantities in Eq. ( 19), we get, Eq. ( 20) is the well known Whittaker's equation (Whittaker and Watson, 1990).
The solution Eq. ( 20) is given by where  and 1  are arbitrary constants and The solution of Eq. ( 22) is given by Eq. ( 23) is the displacement for the Love wave in the half space.Now, expanding Eq. ( 23) up to linear term, we have  w are the components of the displacement vector and 2  is the density of the intermediate layer.
The stress displacement relation for isotropic media is where  and 2  are known as Lame's constants for homogeneous media, ij  is Kronecker delta and where, i u are the components of displacement vector and can be defined as , the stress-strain relations are Using Eq. ( 28), the Eq. ( 25) can be written as where 2 2 2 c    .To solve Eq. ( 29) we take the following substitution where kc   , c is phase velocity and k is wave number.Using Eq. (30) in Eq. ( 29), we get where Therefore, the solution for the intermediate layer is given by where C and D are arbitrary constants.

BOUNDARY CONDITIONS
Both displacement and stress components are continuous at

DISPERSION RELATIONS
Applying 2 nd boundary conditions in equations ( 2), (28) and equations ( 9), (33), we have Now, applying 1 st boundary conditions in equations ( 15), (28) and equations ( 24), (33), we have Eliminating A, B, C and D from equations ( 34) to (37), the dispersion relation for Love waves can be calculated as On solving further above equation, we get On solving further equation ( 38), we get

which implies
Latin American Journal of Solids and Structures 12 (2015) 1934-1949 Eq. ( 39) is dispersion relation of Love waves in an intermediate isotropic vertical layer placed in between orthotropic and prestressed inhomogeneous half-spaces.

Special cases
Case 1 If Eq. ( 40) is dispersion relation of Love waves in an intermediate isotropic vertical layer placed in between homogeneous and prestressed inhomogeneous half-spaces.
Case 2 If 0   , the Eq. ( 39) reduces to Eq. ( 41) is dispersion relation of Love waves in an intermediate isotropic vertical layer placed in between orthotropic and prestressed homogeneous half-spaces.Case 3 If the Eq. ( 39) reduces to Eq. ( 42) is dispersion relation of Love waves in an intermediate isotropic vertical layer placed in between homogeneous and prestressed homogeneous half-spaces.
Case 4 If 0   , 0 P  , the Eq. ( 39) reduces to Eq. ( 43) is dispersion relation of Love waves in an intermediate isotropic vertical layer placed in between orthotropic and homogeneous half-spaces.
Case 5 If Also, on neglecting the lower half space, Eq. ( 44) reduces to Eq. ( 45) and Eq. ( 46) are classic Love wave dispersion relation; hence it validates our solution for Love waves in an intermediate isotropic vertical layer placed in between orthotropic and prestressed inhomogeneous half-spaces.

NUMERICAL CALCULATIONS AND DISCUSSION
To show the effect of inhomogeneity parameters and initial stress parameters of lower half-space on Love wave propagation in intermediate layer, we take following parameters Gubbins (1990).
(i) Material Parameters for upper half-space.
(ii) Material Parameters for intermediate layer.
(iii) Material Parameters for lower half-space.
We have plotted dimensionless phase velocity   0.5 present in the lower half-space.The values of stress parameters for curves have been taken as 0.0, 0.3, 0.5 and 0.7, respectively.It is observed from these curves that as the stress parameters in the half-space increases, the velocity of Love wave increases.Figure 4 shows the effect of initial stress parameters on dimensionless phase velocity

CONCLUSIONS
In this paper, an analytical approach is used to investigate the propagation of Love wave in a homogeneous isotropic layer of finite thickness between orthotropic and inhomogeneous half spaces.

Figure 1 :
Figure 1: Geometry of the problem.
     and zz  are the incremental stress components, 3 u , 3 v and 3


are the Whittaker functions.Now considering the condition V( )0 z  as z   i.e. ( ) 0 s  as s   in Eq. (21), the exact solution becomes z   , therefore the geometry of the problem leads to the following conditions:1 st Boundary conditions (i) At the interface, 0 z the continuity of the displacement along the x direction requires that2 3 v v , where 2 v and 3 v are the displacement components, along the y direction only, in the intermediate layer and lower half-space respectively.(ii) At the interface, 0 z  the continuity of the stress requires that  stress component.(iii) Also, stability conditions leads to 3 0 v  as z   .

2
nd Boundary conditions (i) At the interface, z   , the upper boundary plane is not free surface, the continuity of the displacement along the x direction requires that 1 2 v v  , where 1 v is the displacement component in the upper half-space along the y direction only.(ii)At the interface, z   , the continuity of the stress requires that  the relevant is stress component.(iii) Also, stability conditions leads to 1 0 v  as z   .

2cc
against dimensionless wave number kH for Eq.(39) using MATLAB software.The effects of initial stress parameters   2 P  and inhomogeneity parameters k  on Love wave propagation have been shown in Figs.2-4. Figure 2 is plotted for dimensionless phase velocity2 c c in intermediate layer against dimensionless wave number kH of Love wave for different values of inhomogeneity parameter k  and in the presence of constant initial stress parameter   2 P   0.5 present in the lower half-space.It is clear from this figure, the phase velocity increases with increase of inhomogeneity parameters k  .Figure 3 represents the variation of dimensionless phase velocity 2 c c in intermediate layer against dimensionless wave number kH of Love wave for different values of initial stress parameter   2 P  and in the presence of constant inhomogeneity parameter k

2
c c in intermediate layer against dimensionless wave number kH of Love wave in the presence of constant inhomogeneity parameter for homogeneous media.From above numerical analysis, the following observations are made: i.In entire figures, dimensionless phase velocity 2 c c of Love waves in intermediate layer decreases with increase of dimensionless wave number kH .ii.The dimensionless phase velocity 2 c c of Love wave in intermediate layer shows remarkable change with inhomogeneity k  and stress parameters   2 P  .iii.It is observed as the depth increases the velocity of Love wave in intermediate layer decreases.iv.The phase velocity 2 c c of Love wave in intermediate layer decreases with the decrease of initial stress   2 P  of lower half-space.v.The phase velocity 2 c c of Love wave in intermediate layer increases with the increase of inhomogeneity parameter k  of lower half-space.