Reflection of Plane Waves from Free Surface of an Initially Stressed Rotating Orthotropic Dissipative Solid Half-Space ()
1. Introduction
The Earth is considered as an elastic body with various additional parameters, e.g. porosity, initial stress, viscosity, dissipation, temperature, voids, diffusion, etc. Initial stresses in a medium are caused by various reasons such as creep, gravity, external forces, difference in temperatures, etc. The reflection of plane waves at free surface, interface and layers is important in estimating the correct arrival times of plane waves from the source. Various researchers studied the reflection and transmission problems at free surface, interfaces and in layered media [1- 12]. The study of the reflection of plane waves in the presence of initial stresses and dissipation finds significant applications in various engineering fields. Following Biot [13] theory of incremental deformation, Selim [14] studied the reflection of plane waves at a free surface of an initially stressed dissipative medium. Recently, Singh and Arora [15] studied the reflection of plane waves from a free surface of an initially stressed transversely isotropic dissipative medium.
In the present paper, we studied the problem on reflection of plane waves at a stress-free surface of an initially stressed rotating orthotropic dissipative solid half-space. The reflection coefficients of reflected waves are computed numerically to observe the effects of initial stress and rotation.
2. Formulation of the Problem and Solution
We consider an initially stressed orthotropic half-space rotating about y-axis with. Following Biot [13] and Schoenberg and Censor [16], the basic dynamical equations of motion in x-z plane for an infinite, initially stressed and rotating medium, in the absence of external body forces are,
(1)
(2)
where r is the density, is rotational component, sij (i, j = 1, 3) are incremental stress components, u and w arethe displacement components.
Following Biot [13], the stress-strain relations are:
(3)
where Cij are the incremental elastic coefficients.
For dissipative medium, elastic coefficients are replaced by the complex constants:
(4)
where, , are real. Following Fung [17], the stress and strain components in dissipative medium are,
(5)
where (i, j = 1, 3) and being the angular frequency.
With the help of Equations (4) and (5), the Equation (3) becomes,
(6)
The displacement vector is given by where (n) assigns an arbitrary direction of propagation of waves,
is the unit displacement vector and
is the phase factor, in which is the unit propagation vector, cn is the velocity of propagation, , and kn is corresponding wave number, which is related to the angular frequency by. The displacement components u(n) and w(n) are written as
(7)
Making use of Equations (6) and (7) into the Equations (1) and (2), we obtain a system of two homogeneous equations, which has non-trivial solution if
(8)
where,
, ,
The rootscorrespond to quasi-P (qP) waves and quasi-SV (qSV) waves respectively. These two roots give the square of velocities of propagation as well as damping. Real parts of the right hand sides correspond to phase velocities and the respective imaginary parts correspond to damping velocities of qP and qSV waves, respectively. It is observed that both and depend on initial stresses, rotation, damping and direction of propagation. In the absence of initial stresses, rotation and damping, the above analysis corresponds to the case of orthotropic elastic solid.
3. Reflection of Plane Waves from Free Surface
We consider an initially stressed rotating orthotropic dissipative half-space occupying the region z > 0 (Figure 1). In this section, we shall drive the closed form expressions for the reflection coefficients for incident qP or qSV waves.
The displacement components of incident and reflected waves are as,
(9)
where,
(10)
Here, subscripts 1, 2, 3 and 4 correspond to incident qP wave, incident qSV wave, reflected qP wave and reflected qSV wave, respectively.
In the x-z plane, the displacement and stress components due to the incident qP wave are written as
(11)
where,
In the x-z plane, the displacement and stress components due to the incident qSV wave are written as
(12)
In the x-z plane, the displacement and stress components due to the reflected qP wave are written as
(13)
In the x-z plane, the displacement and stress components due to the reflected qSV wave are written as
(14)
The boundary conditions required to be satisfied at the free surface z = 0,
(15)
The above boundary conditions are written as
(16)
The Equations (11) to (14) will satisfy the boundary conditions (16), if the following Snell’s law holds
(17)
with the following relations A1δ1 + A2δ2 + A3δ3 + A4δ4 = 0, (18)
A1δ5 + A2δ6 + A3δ7 + A4δ8 = 0, (19)
where
(20)
and.
For incident qP wave (A2 = 0), we obtain from equations (18) and (19),
(21)
For incident qSV wave (A1 = 0), we obtain from equations (18) and (19),
(22)
For isotropic case, B11 = λ + 2μ + P, B13 = λ + P, B33 = λ + 2μ, Q = μ, P = –S11, Ω = 0, then the above theoretical derivations reduce to those obtained by Selim [14]
4. Numerical Example
For numerical purpose, a particular example of the material (Zinc) is chosen with the following physical constants,
, ,
, ,
, ,
, ,
.
From Equations (21) and (22), the reflection coefficients of reflected qP and qSV waves are computed for the incident qP and qSV waves. The numerical values of the reflection coefficients of reflected qP and qSV waves are shown graphically in Figures 2 and 3 for incident qP wave and in Figures 4 and 5 for incident qSV wave.
In Figure 2, from comparison of solid line with dashed lines, it is observed that the reflection coefficients of qP and qSV waves change due to the presence of initial stresses at each angle of incidence of qP wave except
Figure 2. Effect of initial stresses on the reflection coefficients of qP and qSV waves for incidence of qP wave.
Figure 3. Effect of rotation parameter on the reflection coefficients of qP and qSV waves for incidence of qP wave.
Figure 4. Effect of initial stresses on the reflection coefficients of qP and qSV waves for incidence of qSV wave.
Figure 5. Effect of dissipation on the reflection coefficients of qP and qSV waves for incidence of qSV wave.
grazing incidence. The effect of initial stresses is observed maximum in the range 45˚ < e1 < 90˚.
From Figure 3, it is observed that the reflection coefficients of qP and qSV waves change due to the presence of rotation in the medium at each angle of incidence of qP wave except grazing incidence.
In Figure 4, from comparison of solid line with dashed lines, it is observed that the reflection coefficients of qP and qSV waves change due to the presence of initial stresses in the medium at each angle of incidence of qSV wave except grazing incidence.
From Figure 5, it is observed that the reflection coefficients of qP and qSV waves change due to the presence of dissipation in the medium at each angle of incidence of qSV wave.
5. Conclusions
The reflection from the stress-free surface of an initially stressed rotating orthotropic dissipative medium is considered. The expressions for the reflection coefficients of reflected qP and qSV waves are obtained in closed form for the incidence of qP and qSV waves. For a particular material, these coefficients are computed and depicted graphically against the angle of incidence for different values of initial stress and rotation parameters. From the figures, it observed that 1) the initial stresses affect significantly the reflection coefficients of reflected qP and qSV waves; 2) the rotation parameter also affects significantly the reflection coefficients of qP and qSV waves; 3) Reflection coefficients are also affected due to the presence of dissipation.