 Open Access
 Total Downloads : 7
 Authors : Sourav Debasis, Sapan K.Samal, Ranjan Chattaraj
 Paper ID : IJERTCONV3IS25008
 Volume & Issue : NCRAEEE – 2015 (Volume 3 – Issue 25)
 Published (First Online): 30072018
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Love Waves Propagation in a TransverseIsotropic Fluid Saturated Porous Layer Lying Over a SelfReinforced Half Space
Sourav Debasis1, Sapan K.Samal2 ,Ranjan Chattaraj3
1 Department of Mathematics,
Gandhi Institute For Technology, Gangapada, BBSR752054,India.
2 Department of Mathematics,
School of Applied Sciences, KIIT University, Bhubaneswar751024,Odisha, India.
3 Department of Applied Mathematics,
Birla Institute of Technology, Mesra, Deoghar Campus814142, India.
Abstract: This paper reveals propagation of Love waves in a fluid saturated anisotropic porous layer lying over a self reinforced half space. Transverse isotropic fluid flow resistivity is taken into consideration for porous layer. The equations of motion have been formulated separately for different media under suitable boundary conditions at the interface of porous layer and selfreinforced half space. We obtain the dispersion relation for assuming model and the result shows there is attenuation in propagation of Love waves which arises due to the resistivity, porosity and reinforced parameters. The effect of porosity, resistivity and reinforced parameter in phase velocity is studied and shown graphically. Some particular cases are also derived and compared with results obtained earlier.
Keywords: Love waves, Porous Layer, Transverse Isotropic, Selfreinforced, Attenuation, Phase velocity.

INTRODUCTION
Studies of propagation of elastic waves in layered media have long been a research subject because of its practical importance in geophysics, earthquake engineering and exploration of oil and underground water. A good amount of information about the propagation of seismic waves in layered media is available in the well known book of Ewing et al.(1957). Konczak Z.(1989) discussed the propagation of Love waves in a fluid saturated porous anisotropic layer and he showed that attenuation in the propagation of surface waves occurs due to porosity of the layer. The propagation of Love wave in layered anisotropic media has been discussed by Kuznetsov S V.(2006), Yongqiang Guo et al.(2008) and others. Recently the study of wave propagation in porous medium has gained prime interest.
The earth contains fluid saturated porous rocks on or below its surface in the form of sandstone and other sediments permeated by groundwater or oil. Water saturated ocean sediments may also be considered as porous material. Biot M.A.(1956a,b) developed the theory of plane wave propagation in fluid saturated porous media in his two classical papers. Later Deresiewicz H.(1961), studied the effect of boundaries on the propagation of waves in liquid filled porous solid and derived the dispersion equation for the Love waves in a porous solid. Practically, the saturated porous materials are anisotropic due to bed, compaction and the presence of aligned micro cracks. Anisotropy may
have significant effects on wave characteristic in layered media. Ghorai et al.(2010) studied the effect of gravity on propagation of Love wave in a porous layer under rigid boundary. Chattaraj et al(2013) discussed the dispersion of Love wave propagating in irregular anisotropic porous stratum under initial stress. Again Chattaraj et al(2013) shows the dispersion of torsional surface waves in an anisotropic layer over the porous half space under gravity. Pradhan et al (2002) discuss the dispersion of shear waves in a fluid saturated elastic plate, taking resistivity of fluid and obtained some useful result. Samal et al discussed surface wave propagation in fiberreinforced anisotropic elastic layer between liquid saturated porous halfspace ad uniform liquid layer. Wang et al. (1998) studied the propagation of Love wave in a transversely isotropic fluid saturated porous layered half space. Later Ke et al. (2006) studied the same problem with a linearly varying elastic property.
Maradudinad et al. (1976) studied the attenuation of Rayleigh surface waves by surface roughness. Billy M.de et al. (1987) analyzed the attenuation measurements of an ultrasonic Rayleigh wave propagation along rough surfaces.
Fiber reinforced materials are widely used in engineering structures due to its high strengthtoweight and stiffnesstoweight ratio, which makes them ideally suited for use in weight sensitive structures. As the mechanical behavior of many fiberreinforced composite materials is adequately modeled by the theory of linear elasticity for transversely isotropic materials in the preferred direction, coinciding with the fiber direction, the fiber reinforced composites are assumed to exhibit anisotropy. Also, some hard and soft rocks or minerals beneath the earth surface exhibit fiber reinforcement property. Pipkin and Rogers(1971) developed the plane strain theory of finite deformation for fiber reinforced materials. Spencer(1972) introduced the constitutive equation for a fiber reinforced linearly anisotropic elastic medium with respect to preferred direction. The concept of continuous reinforcement at every point of an elastic solid was given by Belfield et al. (1983). Hashing and Rosen(1964) gave the elastic moduli for fiber reinforced materials. Chattopadhyay and Michel(2006) studied a model for spherical second harmonic SH wave propagation in selfreinforced linearly elastic media. A study of torsional waves in a fiber reinforced composite material has been done by Bao et al. (2006). Chattopadhyay et al.(2009) discussed the propagation of Torsional waves in semi infinite fiber reinforced cylindrical rod. Chattopadhyay et al. (2012) studied the torsional surface waves in fiber reinforced layer lying over inhomogeneous half space with linearly varying rigidity and density and showed that the presence of reinforcement increases the phase velocity.
In the present paper, we discuss the propagation of Love waves in a fluid saturated anisotropic porous layer lying over a reinforced half space. The transverse isotropic fluid flow resistance in porous layer is considered in the porous medium which first of all taken by Konczak Z.(1989). Further, very few have considered the resistivity in fluid saturated porous layer over a reinforced halfspace due to complexity in getting solution. The equations of motion have been formulated for both the media under suitable boundary conditions at the interface. Following Biot (1956a, b) we derive the dispersion equation for love waves. In this problem we observed that the attenuation of propagation of Love waves occurred due to the effect of anisotropy, porosity and resistivity of the medium. Love waves for a particular model has been discussed and shown graphically. It is observed that the phase velocity gradually decreases with increase in wave number for a fixed value of attenuation parameters. The attenuation also decreases with the increase of frequency.

FORMULATION
We consider a fluid saturated anisotropic porous layer of thickness h lying over a selfreinforced halfspace. The free surface of the porous layer is assumed to be traction free. A coordinate system ( x, y, z ) has been taken with
origin at the interface and positive axial ( z) direction is directed downward. The wave is assumed to propagate along x direction (Fig1).
Anisotropic porous layer
h O x
Fiber reinforced halfspace
z
Fig.1: Geometry of the problem

Dynamics of Porous layer
The equations of motion for the fluidsaturated porous layer without body forces are (1962)
ij,j = 11u i + 12U i – bij(U j – u j) (1)
,i = 12u i + 22U i – bij(U j – u j) (2)
Where aij are stress components of solid skeletn, a = fp is the reduced pressure of the fluid( p is the pressure of the fluid and f is the porosity factor of the medium), ui and Uj are displacement components of solid and fluid in porous media, p11, p12, and p22 are dynamics coefficients arise due to the inertia effect of fluid flow and are related to the mass density of solid ps and fluid pf .
b11
0
0
[bij] = 0 b11
0
0
0
b33
b11
0
0
[bij] = 0 b11
0
0
0
b33
The components of the flow resistance tensor b , for the transverseisotropy are
(3)
Using the conventional Love waves conditions
u1 = 0, u2 = u2(x, z, t), u3 = 0 U1 = 0, U2 = U2(x, z, t), U3 = 0
Equation (1) and (2) reduces to
N a2u2 + G a2u2 = a2 (p
u + p
U ) – b
a (U
– u ) (4)
ax2
az2
at2
11 2
12 2
11 at 2 2
0 = a2 (p
u + p
U ) + b
a (U
u ) (5)
at2
11 2
12 2
11 at 2 2
By eliminating U2 from these equations, we obtain
a2u2
au2 2
N a2u2 + G a2u2 – p
a2u2 – b
au2 – p12 at2 b11 at
= 0 (6)
ax2
az2
11 at2
11 at
p a2u2+b au2
22 at2 11 at
For the waves changing harmonically, we take u2 = 0(z)ei(kxwt) (7)
2
2
The equation (6) takes the form ( a2
az
+ 2)< = 0 (8)
1
1
1
1
Where 2 is a complex quantity and is defined by
= k C2 F + i C2 R – N = k( , and ( = C2 F + i C2 R – N
(9)
C
C
C
C
C
C
C
C
G
G
G
G
1 2 2 2 2
G G G G
F = F( ) = w2 (C1 222+1) 22 , R= R( ) = w2 (C1+ 22) 22 ,c =
– 2 ,
C
C
C
C
G
G
G
G
2 (1+2 222)C1
2 (1+2 22)C1 1
p2 pw C1
11 22 12
kl = pkl, k, l = 1,2. c2 = G/p , p = p11 – ( 12) , =
, d =
G
This can be drawn that
p22
b11
22
as f 0 i. e. d 1, the porous layer becomes non – porous solid as f 1 i. e.d 0, the porous layer becomes fluid
as 0 < d < 1, the layer is poro – elastic
The solution of (8) is u2 = (A1 eik( + B1 eik()ei(kxwt) (10)

Dynamic of Fiber reinforced material:
According to Spencer A.J.M. (1972), the Cauchy stress tensor for self reinforced elastic medium is given by
3
3
*
*
* 3
* *
ij
k 1
ekk ij
2 Teij
3
3
3
3
k 1 m1
akamekm ij
k 1
aia jekk
(11)
2
a
a e* a
a e*
a a e* a a
3
3
L T
i
3
3
k 1
k kj
j
k 1
k ki
3
3
3
3
k 1 m1
k m km i j
where e* (u* u* ) / 2 , u* is the components of displacement, a a T
is the preferred fiber direction
ij i, j j ,i j i
( i, j 1,2,3 which corresponds to x, y, z ). The comma before an index represents (partial) space differentiation.
Here * are the components of stress, e* components of infinitesimal strain, a the components of a and all refer
ij ij i
to rectangular Cartesian coordinates. The coefficients , T , L ,,
are elastic constant with the dimension of
stress. T
and
L are the transverse shear and longitudinal shear modulus in the preferred direction respectively.
and are the specific stress components to take into account different layers for concrete part of the composite material.
The equation of motion without body force are given by
P a2V + 2iR
a2V
+ Q a2V = a2V
(12)
1 ax2
1 axaz
1 az2 T at2
Where 1 = 1 + (Âµ – 1) a12 , R1 = (Âµ – 1) a1a3 , 1 = 1 + (Âµ – 1)a2
ÂµT ÂµT
ÂµT 3
The solution of (12) is v = A2 ekm1z ei(kxwt) , (13)
where m = i (R
C
C
+ R2 + ( – C2) ) , c2 = ÂµT , p is the density of self reinforced medium
1 Q1 1
1 1 1 2
2
2 p

Boundary Conditions are
a23 = 0 at z = h
a23 = r23 at z = 0 (14)
u2 = v at z = 0
Using these boundary conditions, finally we obtain the dispersion equation
tan(k(h) = m1D , where D = ÂµL + (ÂµL – ÂµT)a2 (15)
G( 3
Now taking k = k1 + ia , where a is the attenuation coefficient, we obtain the real part of the dispersion equation(Using Mathematica software) in the form
2 4
2 4
1
C
C
C
C
G
G
k h C2 F – N + C2 R2 cos 1
C2 F – N +i C2 R) – ah
1 2 2
G G
2 4
2 4
1
2 Arg( 2 2
G G
C
C
C
C
G
G
C2 F – N + C2 R2 sin 1 C2 F – N +i C2 R) = 1 Ar (1 – D1 ) Ar (1 +
C2 G C2
2 Arg( C2
G C2 2 2 2
C
C
C
C
G G G G
C FN+i C R G
2 G 2
G G
D1
C2FN+i C2R G
) ] (16)
C
C
C
C
2 G 2
G G
Where D1 = m1D.

Particular Cases:
Case1: In case the resistivity of fluid b11 = 0 i.e. F=1 and R=0 , then the dispersion equation takes the form
(Âµ +(Âµ Âµ )a2) R + R2+Q P C2
C
C
tan k h C2 – N =
T 3 1
1 1 1 2
2
(17)
C
C
2
2
G
G
1 2
G Q1G C N
C
C
2 G G
The term containing the attenuation term in equation (16) becomes zero and this is the dispersion equation when love wave propagates through a fluid saturated porous layer without fluid resistance rest upon a reinforced half space. When reinforcement parameter is zero i.e. ÂµL = ÂµT = Âµ, R1 = 0, 1 = 1 = 1, then Half space becomes an isotropic medium. The dispersion equation becomes
Âµ 1C2
C
C
C
C
tan k h C2
2
2
2
G
G
– N = 2
(18)
1 2
G G C N
C
C
2 G G
This is similar to the dispersion equation obtained for the prorogation of Love wave in a fluid saturated porous layer lying over an isotropic halfspace.
Case2: If the upper layer becomes isotropic i.e. N=G=Âµ and c1 = cG = jÂµ/p , then dispersion equation takes the standard form as
1C2
C
C
tan k h C2
2
– 1 = 2
(19)
C
C
1 2
2
2
1 C21
C1
This is the dispersion equation of Love waves in a isotropic layer lying over a isotropic halfspace.

Numerical Results
To see the effect of porosity, resistivity of fluid and reinforcement parameters we have taken a particular model. The model consists of a fiber reinforced half space lying below a kerosene sandstone layer.
For this model we calculate the phase velocity ratio i.e. c/cG which is a function of the real part of nondimensional wave number h .
The material constants for a reinforced medium has been considered as
ÂµL = 5.66 x 109N/m2, ÂµT = 2.46 x 109N/m2, p = 7800 k /m2
Moreover, we take < = rr/3 so that a1 = sin (1[), a2 = cos (1[) and other data are specified on the body of the
3 3
graphs.
The elastic parameters for kerosene saturated porous sandstone are taken as Yew and Jogi (1976)
N = 2, p = 1.926137 x 103 k /m2,p
= 0.002137 x 103k /m2,
G 1 12
p13 = 0.215337 x 103 k /m2
Figure 2 shows the attenuation of Love waves decreases as the frequency increases. As d increases, the attenuation also increases. The difference is quite distinct for smallest frequency.
Figure 3 shows the variation of phase velocity with respect to wave number of different frequencies with attenuation ah = 0. It is observed that as the wave number increases the magnitude phase velocity increases. After a certain wave number the phase velocity becomes constant.
Figure 4 also depicts the variation of phase velocity with respect to wave number for different frequencies with attenuation h = 0.2 . It is observed that the magnitude of phase velocity is more in comparison to figure 3. But as the wave number increases the phase velocity decreases like as figure 1.
Figure 5 shows a variation of phase velocity with respect to porosity d for attenuation factor ah = 1,2,3. But it is interesting to note that the phase velocity of Love waves remains same when d=1 i.e. the porous layer becomes solid.
Figure 6 presents the phase velocity variation with respect to wave number for the elastic constant ratio N = 2,3,4
G
of porous medium. As N increases the phase velocity increases. But the increase in wave numbers also increases
G
phase velocity.

Conclusions
There is a significant effect of attenuation, porosity and anisotropy simultaneously in the propagation of Love waves in a transverse isotropic fluid saturated porous layer rest upon a selfreinforced half space. The attenuation is more for low frequency and less for high frequency. The porosity factor d also affects the attenuation. As d increases the magnitude of attenuation also increases. The increase in wave number decreases the phase velocity. It is also seen that when the layer becomes a perfect solid the phase velocity remains same for all attenuation factors. The phase velocity increases as the elastic constant ratio of porous medium increases.

BIBLIOGRAPHY

Bao, L., Sakurai, M. and Nakazawa, M. (2006). A study on the torsional wave of fiber reinforced composite materials. Journal of Composite Materials, 40(4): 338391.

Belfield, A.J., Rogers, T.G. and Spencer, A.J.M.,(1983). Stress in elastic plates reinforced by fibers lying in concentric circles. Journal of Mechanics and Physics of Solids.
31: 2554.

Billy M. de, G Quentin and E Baron (1987) Attenuation measurements of an ultrasonic Rayleigh wave propagation along rough surfaces, Journal of Applied Physics, 61: 21402145.

Biot, M.A., (1956.a). Theory of propagation of elastic waves in a fluidsaturated porous solid: I. Lowfrequency range. ,Journal of Acoustical Society of America, 28(2): 168178.

Biot, M.A.,(1956.b). Theory of propagation of elastic waves in a fluidsaturated porous solid: II. Highfrequency range., Journal of Acoustical Society of America, 28(2):178191.

Biot M.A.,(1962) Mechanics of deformation and acoustic propagation in porous medium, Journal of Applied Physics, 33:14821498

Chattopadhyay, A. and Michel, V., (2006). A model for spherical SHwave propagation in self reinforced linearly elastic media. Archive of Applied Mechanics. 75,23: 113124.

Chattopadhyay, A., Gupta, S., Samal, S.K. and Sharma, V.K.,(2009). Torsional Wave in Self reinforced Medium. International Journal of Geomechanics, 9(1):913.

Chattopadhyay, A., Gupta, S., Sahu, S.A. and Singh, A.K.,(2012). Torsional surface waves in a SelfReinforced medium over a heterogeneous half space, International Journal of Geomechanics, 12(2): 193197.

Chattaraj Ranjan, Samal K. Sapan, and N.C.Mahanti(2013), Dispersion of Love wave propagating in irregular Aisotropic Porous Stratum under initial stress, International Journal of Geomechanics, 13(4): 402408.

Chattaraj Ranjan, Samal K. Sapan, and N.C.Mahanti(2013), Dispersion of Torsional surface waves in anisotropic layer over porous half space under gravity, ZAMM Journal of Applied Mathematics and Mechanics, 19

Deresiewicz H.,(1961),The effect of boundaries on wave propagation in liquidfilled porous solid: II Love waves in a porous layer, Bulletin Seismological Society of America, 5159

Ewing, W.M., Jardetzky, W.S. and Press, F.,(1957). Elastic waves in layered Media. Mc. Graw Hill , New York, Ghorai

P. , Samal S. K. Mahanti N.C.,(2010), Love waves in a fluidsaturated Porous layer under a rigid boundary and lying over an elastic halfspace under gravity, Applied Mathematical Modeling 34:18731883


Hashin, Z. and Rosen, B.W., (1964). The Elastic Modulli of FiberReinforced Materials. Journal of Applied Mechanics, 21: 233242.

Ke, L.L., Wang, Y.S. and Zhang, Z.M., (2006), Love waves in an inhomogeneous fluid saturated porous layered half space with linearly varying properties; Soil Dynamics and Earthquake Engineering; 26: 574581.

Konczak Z.,(1989).The Propagation of Love Waves in a FluidSaturated Porous Anisotropic Layer, Acta Mechanica 79: 155168.

Kuznetsov S. V.,(2006), Love waves in layered anisotropic media, Journal of Applied Mathematics and Mechanics,70: 116127

Maradudinad A. A., Mills D. L.(1976) Attenuation of Rayleigh surface waves by a surface roughness. Annals of Physics, 100: 262309

Pipkin, A.C. and Rogers, T.G., 1971. Plane deformations of incompressible fiber reinforced materials. Journal of Applied Mechanics, 38: 634640.

Pradhan A., Samal S. K.,Mahanti N.C.,(2002), Shear waves in a fluid saturated elastic plate, Sadhana, 27(6), 595 604.

Samal S.K., and Chattaraj R.(2011), Surface wave propagation in fiberreinforced anisotropic elastic layer between liquid saturated porous half space and uniform liquid layer, Acta Geophysica,59(3),470482

Spencer, T.G.,(1972).Deformation of fiberreinforced material, Oxford University Press, London.

Wang Y. S. , Zang Z.M.. (1998),Propagation of Love waves in a transversely isotropic fluid saturated porous layered halfspace, Journal of Acoustical Society of America, 103(2):695701

Yew, C.H., ad Jogi, P.N. (1976), Study of wave motions in fluidsaturated porous rocks, Journal of Acoustical Society of America, 60(1): 28

Yongqiang Guo .,Chen Weiqiu( 2008), On free wave propagation in anisotropic layered media, Acta Mechnica solida sinica, 21(6): 500506