PROPAGATION OF WAVES AT AN IMPERFECTLY

Similar documents
Reflection of quasi-p and quasi-sv waves at the free and rigid boundaries of a fibre-reinforced medium

Reflection and refraction of thermoelastic plane waves at an interface between two thermoelastic media without energy dissipation

Chapter I. Introduction

International Journal of Pure and Applied Sciences and Technology

Propagation of Rayleigh Wave in Two Temperature Dual Phase Lag Thermoelasticity

Time Harmonic Inclined Load in Micropolar Thermoelastic Medium Possesing Cubic Symmetry with One Relaxation Time

LASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE

ON THE FLEXURAL AND EXTENSIONAL THERMOELASTIC WAVES IN ORTHOTROPIC PLATES WITH TWO THERMAL RELAXATION TIMES

Reflection of plane micropolar viscoelastic waves at a loosely bonded solid solid interface

Reflection of SV- Waves from the Free Surface of a. Magneto-Thermoelastic Isotropic Elastic. Half-Space under Initial Stress

DAMPING OF GENERALIZED THERMO ELASTIC WAVES IN A HOMOGENEOUS ISOTROPIC PLATE

The effect of rigidity on torsional vibrations in a two layered poroelastic cylinder

Dr. Parveen Lata Department of Basic and Applied Sciences, Punjabi University, Patiala, Punjab, India.

Research Article Propagation of Plane Waves in a Thermally Conducting Mixture

Receiver. Johana Brokešová Charles University in Prague

Shear Stresses and Displacement for Strike-slip Dislocation in an Orthotropic Elastic Half-space with Rigid Surface

Research Article Dispersion of Love Waves in a Composite Layer Resting on Monoclinic Half-Space

Stoneley Waves at the Boundary Surface of Modified Couple Stress Generalized Thermoelastic with Mass Diffusion

Mathematical Model for Thermal Shock Problem of a Generalized Thermoelastic Layered Composite Material with Variable Thermal Conductivity

REFLECTION AND TRANSMISSION OF PLANE WAVES AT AN INTERFACE BETWEEN ELASTIC AND MICROPOLAR THERMOELASTIC DIFFUSION MEDIA

THE REFLECTION PHENOMENA OF SV-WAVES IN A GENERALIZED THERMOELASTIC MEDIUM

EFFECT OF COUPLE-STRESS ON THE REFLECTION AND TRANSMISSION OF PLANE WAVES AT AN INTERFACE

Internal Heat Source in Temperature Rate Dependent Thermoelastic Medium with Hydrostatic Initial Stress

EFFECT OF INITIAL STRESS ON THE REFECTION OF MAGNETO-ELECTRO-THERMO-ELASTIC WAVES FROM AN ISOTROPIC ELASTIC HALF-SPACE

Part 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA

Iranian Journal of Mathematical Sciences and Informatics Vol.2, No.2 (2007), pp 1-16

On propagation of Love waves in an infinite transversely isotropic poroelastic layer

Rajneesh Kumar, 1 S. C. Rajvanshi, 2,3 and Mandeep Kaur 2,3. 1. Introduction

Wave propagation in temperature rate dependent thermoelasticity with two temperatures

Waves propagation in an arbitrary direction in heat conducting orthotropic elastic composites

REFLECTIONOFPLANEWAVESFROMAFREESURFACEOF A GENERALIZED MAGNETO-THERMOELASTIC SOLID HALF-SPACE WITH DIFFUSION

Propagation of Plane Waves in a Generalized Thermo-magneto-electro-elastic Medium

Reflection of Plane Waves from a Rotating Magneto Thermoelastic Medium with Two Temperature and Initial Srtress Under Three Theories

EFFECT OF DISTINCT CONDUCTIVE AND THERMODYNAMIC TEMPERATURES ON THE REFLECTION OF PLANE WAVES IN MICROPOLAR ELASTIC HALF-SPACE

Waves in Linear Optical Media

SECOND ORDER ASYMPTOTIC BOUNDARY CONDITIONS FOR

Generalized Thermoelasticity Plane Waves in Rotating Media with Thermal Relaxation under the Temperature Dependent Properties

' ' ' ' ). The midplane of the plate is chosen to FREE WAVE PROPAGATION IN PLATES OF GENERAL ANISOTROPIC MEDIA

Mathematical modelling of Stoneley wave in a transversely isotropic thermoelastic media

where d is the vibration direction of the displacement and c is the wave velocity. For a fixed time t,

Plane waves in a rotating generalized thermo-elastic solid with voids

c. Better work with components of slowness vector s (or wave + k z 2 = k 2 = (ωs) 2 = ω 2 /c 2. k=(kx,k z )

Propagation of Plane Waves in Micro-stretch Elastic Solid in Special Case

Reflection of Plane Waves from Electro-magneto-thermoelastic Half-space with a Dual-Phase-Lag Model

Available online at ScienceDirect. Procedia Engineering 144 (2016 )

Thermoelastic Interactions without Energy Dissipation Due to Inclined Load

Influence of Irregularity and Rigidity on the Propagation of Torsional Wave

ISSN: X (p); (e)

Propagation and Reflection of Plane Waves in a Rotating Magneto Elastic Fibre Reinforced Semi Space with Surface Stress

Research Article Reflection of Plane Waves in Generalized Thermoelastic Half Space under the Action of Uniform Magnetic Field

Lawrence Berkeley National Laboratory

Decomposition of Waveguides Propagating in Piezoelectric Crystals subject to Initial Fields

Wave Propagation in Uniaxial Media. Reflection and Transmission at Interfaces

Research Article Fundamental Solution in the Theory of Thermomicrostretch Elastic Diffusive Solids

Reflection, refraction, and transmission of SH waves at a micropolar layer separating two elastic media

Available online at ScienceDirect. Procedia Engineering 144 (2016 )

Wave Propagation in Anisotropic, Anelastic, Porous and Electromagnetic Media

COMPARISON OF OPTICAL AND ELASTIC BREWSTER S ANGLES TO PROVIDE INVUITIVE INSIGHT INTO PROPAGATION OF P- AND S-WAVES. Robert H.

21. Fractional order magneto-thermoelasticity in a rotating media with one relaxation time

Static Deformation Due To Long Tensile Fault Embedded in an Isotropic Half-Space in Welded Contact with an Orthotropic Half-Space

Graduate School of Engineering, Kyoto University, Kyoto daigaku-katsura, Nishikyo-ku, Kyoto, Japan.

Relationships between the velocities and the elastic constants of an anisotropic solid possessing orthorhombic symmetry

Module 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method

Numerical Modeling for Different Types of Fractures

Effects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses

Application of fractional order theory of thermoelasticity to a 1D problem for a cylindrical cavity

Linearized AVO in viscoelastic media Shahpoor Moradi,Kristopher A. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada

Effective boundary condition method and Rayleigh waves in orthotropic half-spaces coated by a thin layer with sliding contact

Electromagnetic Waves Across Interfaces

Modeling of Variable Lamé s Modulii for a FGM Generalized Thermoelastic Half Space

Forced Vibration of the Pre-Stressed and Imperfectly Bonded Bi-Layered Plate Strip Resting on a Rigid Foundation

Conversion coefficients at a liquid/solid interface

INTERFACE WAVES ALONG FRACTURES IN TRANSVERSELY ISOTROPIC MEDIA

Cellular solid structures with unbounded thermal expansion. Roderic Lakes. Journal of Materials Science Letters, 15, (1996).

Available online Journal of Scientific and Engineering Research, 2016, 3(6): Research Article

Module III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina

Horizontally polarized shear waves in stratified anisotropic (monoclinic) media

Introduction to Seismology Spring 2008

Unphysical negative values of the anelastic SH plane wave energybased transmission coefficient

Chap. 1 Fundamental Concepts

Radiation pattern in homogeneous and transversely isotropic attenuating media

16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity

Effect of Thermal Stress and Magnetic Field on Propagation of Transverse Wave in an Anisotropic Incompressible Dissipative Initially Stressed Plate

Waveform inversion and time-reversal imaging in attenuative TI media

3D and Planar Constitutive Relations

Effect of Rotation and Initial Magnetic Field in Fibre-Reinforced Anisotropic Elastic Media

Far-field radiation from seismic sources in 2D attenuative anisotropic media

PART A. CONSTITUTIVE EQUATIONS OF MATERIALS

Rayleigh waves of arbitrary profile in anisotropic media

I. INTRODUCTION J. Acoust. Soc. Am. 109 (4), April /2001/109(4)/1398/5/$ Acoustical Society of America 1398

Mode-I crack in a two-dimensional fibre-reinforced generalized thermoelastic problem

Solutions: Homework 7

DISPLACEMENTS AND STRESSES IN AN ANISOTROPIC MEDIUM DUE TO NON-UNIFORM SLIP ALONG A VERY LONG STRIKE-SLIP FAULT

Multilayer Reflectivity

Phase-shift modelling for HTI media

COMPLETE SOLUTIONS OF A COUPLED SYSTEM OF PARTIAL DIFFERENTIAL EQUATIONS ARISING IN THERMOELASTICITY * D. S. CHANDRASEKHARAIAH

Dynamic Response of Heterogeneity and Reinforcement on the Propagation of Torsional Surface Waves

An acoustic wave equation for orthorhombic anisotropy

Electric field dependent sound velocity change in Ba 1 x Ca x TiO 3 ferroelectric perovskites

Chapter 2 - Macromechanical Analysis of a Lamina. Exercise Set. 2.1 The number of independent elastic constants in three dimensions are: 2.

Transcription:

Journal of Theoretical and Applied Mechanics, Sofia, 2011, vol. 41, No. 3, pp. 77 92 PROPAGATION OF WAVES AT AN IMPERFECTLY BONDED INTERFACE BETWEEN TWO MONOCLINIC THERMOELASTIC HALF-SPACES Joginder Singh Department of Mathematics, Indo Global College of Engineering, Abhipur, Mohali, India Baljeet Singh Department of Mathematics, Post Graduate Government College, Sector-11, Chandigarh-160011, India e-mail:bsinghgc11@gmail.com Praveen Ailawalia Department of Mathematics, RIMT Institute of Engineering and Technology, Mandi Govindgarh, Punjab, India [Received 09 February 2011. Accepted 14 November 2011] Abstract. An imperfectly bonded interface between two monoclinic thermoelastic half-spaces is chosen to study the reflection and transmission of plane waves in context of generalized thermoelasticity. Six relations between amplitudes of incident, reflected and transmitted quasi-p (qp) waves, quasi thermal (qt) waves and quasi-sv (qsv ) waves are obtained. Some particular cases are obtained which agree with earlier well established results. A procedure for computing reflected and refracted angles is derived for a given incident wave, where the angles of reflection are found not equal to the angles of incident waves in a monoclinic thermoelastic medium. Key words: generalized thermoelasticity, imperfect boundary, plane waves, reflection and transmission, monoclinic anisotropic. 1. Introduction Keith and Crampin [1] investigated that three types of body waves with mutually orthogonal particle motion can propagate in an anisotropic elastic solid medium. In general, the particle motion is neither purely longitudinal nor * Corresponding author e-mail: bsinghgc11@gmail.com

78 Joginder Singh, Baljeet Singh, Praveen Ailawalia purely transverse. Due to this, the three types of body waves in an anisotropic medium are referred to as qp, qsv and qsh, rather than as P, SV and SH, the symbols used for propagation in an isotropic medium. A monoclinic medium has one plane of elastic symmetry. SH motion is decoupled from the P SV motion for wave propagation in the plane of symmetry. It is neither purely longitudinal nor purely transverse motion in the case of P SV waves, while the particle motion of SH waves is purely transverse. Chattopadhyay and Choudhury [2] discussed the reflection of qp waves at the plane free boundary of amonoclinic half-space. Inasubsequentpaper, Chattopadhyaet al [3] studied the reflection of qsv waves at a plane free boundary of a monoclinic half-space. Inthese papers, they assumedthatqp waves are purelylongitudinal andqsv waves are purely transverse. Singh [4] pointed out the errors in these papers and SinghandKhurana [5] restudiedthe reflection ofp and SV waves at the free surface of a monoclinic half-space. Biot [6] developed the coupled theory of thermoelasticity, where heat equation being parabolic predicts an infinite speed of propagation for thermal wave, which is not physically relevant. Lord and Shulman [7] introduced the theory of generalized thermoelasticity with one relaxation time by postulating a new law of heat conduction to replace the classical Fourier s law. This law contains the heat flux vector as well as its time derivative. In addition, it contains a new constant which acts as a relaxation time. The heat equation of this theory is of wave type, ensuring finite speeds of propagation for thermal and elastic waves. The remaining governing equations for this theory, namely the equation of motion and the constitutive relation remain the same as those for the coupled and uncoupledtheories. Dhaliwal and Sherief [8] extended the Lord and Shulman (L-S) theory for an anisotropic media. The second generalization was developed by Green and Lindsay [9] and is known as G-L theory. This theory contains two constants that act as relaxation times. Many authors contribute towards the wave propagation in isotropic and anisotropic medium with the application of generalized thermoelasticity, for example, Sinha and Sinha [10], Sharma [11], Kumar and Singh [12] and Singh [13 15]. The boundary between the solids may behave as slip, perfect or neither, depending on the properties of the layers, as demonstrated by Rokhlin and Marom [16], and its state significantly affects elastic wave reflection phenomena. Imperfect bonding means that the stress components, heat flux are continuous and the small displacement field and temperature are not continuous. The small vector difference in the displacement is assumed to depend linearly on the traction vector and in temperature it is assumed to depend linearly on the heat flux. In the present paper, the reflection and transmis-

Propagation of waves at an imperfectly bonded interface... 79 sion of plane waves is studied at an imperfectly bonded interface between two dissimilar thermoelastic half-spaces of monoclinic type in the context of gerneralized thermoelasticity. Relations between amplitudes of incident, reflected and transmitted quasi-p (qp) waves, quasi thermal (qt) waves and quasi- SV (qsv ) waves are obtained, which are reduced to earlier well established relations for some limiting cases. A procedure for computing reflected and refracted angles is also derived for a given incident angle of a striking wave. 2. Nomenclature c phase velocity c ij isothermal elasticities C e specific heat at constant strain e ij components of strain tensor k wave number t time τ 0,τ 1 relaxation times K 2,K 3 thermal conductivities T absolute temperature T 0 the initial uniform temperature α 2,α 3 coefficients of linear thermal expansion ρ density σ ij components of stress tensor ω circular frequency Θ = T T 0, Θ T 1 β 2 = (c 12 + c 22 )α 2 + c 23 α 3, β 3 = 2c 23 α 2 + c 33 α 3. 3. Formulation of the problem and solution Considerahomogeneous, anisotropic, generalized thermoalstic medium of monoclinic type at a uniform temperature T 0. The origin is taken on the thermally insulated and stress-free plane surface and z-axis is directed normally into the half-space, which is representedby z >0. Following Lord andshulman [7] and Green ( and Lindsay ) [9] theories, the governing field equations in the y z plane x 0 for generalized monoclinic thermoelasticity, in absence of body forces and heat sources, are: 2 v (1) c 22 y 2 + c 2 v 44 z 2 + c 2 w 24 y 2 + c 2 w 34 z 2 + 2c 2 v 24 y z

80 Joginder Singh, Baljeet Singh, Praveen Ailawalia ( ) + (c 23 + c 44 ) 2 w y z β Θ 2 Θ + τ 1 = ρ 2 v y t t 2, 2 v c 24 y 2 + c 2 v 34 z 2 + c 2 w 44 y 2 + c 2 w 33 z 2 + 2c 2 w 34 y z (2) + (c 23 + c 44 ) 2 v y z β 3 z ( ) Θ Θ + τ 1 = ρ 2 w t t 2, 2 Θ K 2 y 2 + K 2 [ ( Θ 2 ) 3 z 2 T v 0 β 2 y t + τ 0Ω 3 v y t 2 ( 2 )] ( w (3) +β 3 z t + τ 0Ω 3 w Θ z t 2 = ρc e t + τ 2 ) Θ 0 t 2. The use of symbol Ω and τ 1 in equation (3) makes these equations possible for L-S and G-L theories of generalized thermoelasticity. For the L-S theory, τ 1 = 0, Ω = 1 and for the G-L theory τ 1 > 0, Ω = 0. The thermal relations τ 0 and τ 1 satisfy the inequality τ 1 > τ 0 > 0 for the G-L theory only. For plane waves of circular frequency ω, wave number k and phase velocity c, incident at the free boundary z = 0 at an angle θ with the z-axis, we assume: (4) (v, w, Θ) = (A, B, C)exp (ιp 1 ), where A, B, C are the amplitude factors and: (5) P 1 = ωt k (y sin θ z cos θ), is the phase factor. For waves reflected at z = 0, we assume: (6) (v, w, Θ) = (A, B, C)exp (ιp 2 ), where: (7) P 2 = ωt k (y sin θ + z cos θ). (8) Making use of Eq. (4) or (6) in Eqs. (1) to (3), we obtain: ( D1 ρc 2) ( ι A + L 1 B τ k) 1 β 2 sin θ C = 0,

Propagation of waves at an imperfectly bonded interface... 81 (9) L 1 A + ( D 2 ρc 2) ( ι B ± τ k) 1 β 3 cos θ C = 0, ( ι ) (D3 (10) τt 0 c 2 β 2 sin θ A τt 0 c 2 β 3 cos θ B ρc 2 ) τ C e C = 0, k Here, D 1 (θ) = c 22 sin 2 θ + c 44 cos 2 θ 2c 24 sinθ cos θ, D 2 (θ) = c 44 sin 2 θ + c 33 cos 2 θ 2c 34 sinθ cos θ, D 3 (θ) = K 3 cos 2 θ + K 2 sin 2 θ, L 1 (θ) = c 24 sin 2 θ + c 34 cos 2 θ (c 23 + c 44 ) sinθ cos θ, ( ι ( ι τ = τ 0, τ = τ 0 Ω, τ ω) ω) 1 = (1 + ιωτ 1 ), where, the upper sign corresponds to the incident waves and the lower sign corresponds to the reflected waves. Equations (8) (10) in A, B, C admit a non-trivial solution only if the determinant of their coefficients vanishes, i.e.: (11) A 0 ζ 3 + A 1 ζ 2 + A 2 ζ + A 3 = 0 where A 0 = τ, A 1 = [ D 1 τ + D 2 τ + D 4 + εη ( β2 cos 2 θ + sin 2 θ )], A 2 = ( D 1 D 2 L 2 ) 1 τ + D 1 D 4 + D 2 D 4 + εη ( β2 D 1 cos 2 θ + D 2 sin 2 θ ) ± 2L 1 εη β sin θ cos θ, ( A 3 = D 4 D1 D 2 L 2 1), ( ) ( ζ = ρc 2 D3 β 2, D 4 =, ε = 2 T 0 C e ρc e v1 2 ), η = ττ 1 v 2 1, v 2 1 = c 22 ρ, β = ( β3 Following Singh [13], the cubic equation (11) with complex coefficients has three roots ζ 1, ζ 2, ζ 3 corresponding to phase velocities of three plane waves, namely quasi-p (qp), quasi-thermal (qt) and quasi-sv (qsv ) waves, respectively. If we write c 1 j = Vj 1 iω 1 q j, j = 1,2,3, then V j and q j are the speeds of propagation and the attenuation coefficients of qp, qt and qsv waves, respectively. β 2 ).

82 Joginder Singh, Baljeet Singh, Praveen Ailawalia 4. Reflection and transmission Consider a homogeneous, monoclinic generalized thermoelastic halfspace occupying the region z > 0 in imperfectly bonded contact with another homogeneous, monoclinic generalized thermoelastic half-space occupying the region z < 0. Incident qp or qsv or qt waves at the interface will generate reflected qp,qsv,qt waves inthe half-space z >0and transmitted qp,qsv,qt waves in the half-space z < 0 as shown in Fig. 1. Fig. 1. Geometry of the problem showing incident, reflected and transmitted waves at an interface between two dissimilar thermoelastic half-spaces of monoclinic type The required boundary conditions at an imperfectly bonded interface z = 0 are (Lavrentyev and Rokhlin [17]): (i) σ 33 = K n [w w ], (ii) σ 23 = K t [v v ], (12) (iii) K 3 T z = K c [T T ], (iv) σ 33 = σ 33, (v) σ 23 = σ 23, (vi) K 3 T z = K T 3 z,

Propagation of waves at an imperfectly bonded interface... 83 where, K n, K t are the normal and transverse stiffness coefficients. The symbols with primes correspond to the upper half-space z < 0. The appropriate components of the displacement vector and the temperature field which satisfy the boundary condition (12) at z = 0 are as follows: For half-space z > 0, (13a) v = 6 A j e ip j, w = j=1 For half-space z < 0, (13b) v = 9 A j e ip j, w = j=7 6 B j e ip j, T = j=1 9 B j e ip j, T = j=7 6 C j e ip j. j=1 9 C j e ip j, j=7 where [ P j = ω t (y sin e ] j z cos e j ), (j = 1,2,3,7,8,9); c j [ (14) P j = ω t (y sine ] j + z cos e j ), (j = 4,5,6), c j (15) A j = F j B j, C j = F j B j, For j = 1, 2, 3, 7, 8, 9, ( ) D 2j ρc 2 j β 2 sin θ + L 1j β 3 cos θ F j = ( ), D 1j ρc 2 j β 3 cos θ + L 1j β 2 sin θ F j = ( i k )( ) L 2 1j (D 1j ρc 2 j D 2j ρc 2 j ) [ ( ) ], τ L 1j β 2 sin θ + D 1j ρc 2 j β 3 cos θ where, D 2j = c 44 sin 2 e j + c 33 cos 2 e j 2c 34 sine j cos e j, D 1j = c 22 sin 2 e j + c 44 cos 2 e j 2c 24 sine j cos e j, L 1j = c 24 sin 2 e j + c 34 cos 2 e j (c 23 + c 44 ) sine j cos e j.

84 Joginder Singh, Baljeet Singh, Praveen Ailawalia For j = 4, 5, 6, ( ) D 2j ρc 2 j β 2 sin θ L 1j β 3 cos θ F j = ( ), D 1j ρc 2 j β 3 cos θ L 1j β 2 sin θ ( )( ) D 1j ρc 2 Fj = j D 2j ρc 2 j L 2 1j ) [( ) ], τ D 1j ρc 2 j β 3 cos θ β 2 sin θ L 1j where, ( i k D 2j = c 44 sin 2 e j + c 33 cos 2 e j + 2c 34 sine j cos e j, D 1j = c 22 sin 2 e j + c 44 cos 2 e j + 2c 24 sine j cos e j, L 1j = c 24 sin 2 e j + c 34 cos 2 e j + (c 23 + c 44 ) sine j cos e j. Here, the subscriptj =1 corresponds to incident qp waves, j =2 corresponds to incident qt waves, j = 3 corresponds to incident qsv waves, j = 4 corresponds to reflected qp waves, j = 5 corresponds to reflected qt waves, j = 6 corresponds to reflected qsv waves, j = 7 corresponds to transmitted qp waves, j = 8 corresponds to transmitted qt waves and j = 9 corresponds to transmitted qsv. The boundary conditions at z = 0 given by equation (12) must be satisfied for all values of y, then we have: (16) P 1 (y, 0, t) = P 2 (y, 0, t) = P 3 (y, 0, t) = P 4 (y, 0, t) = P 5 (y, 0, t) which can be written as: (17) = P 6 (y, 0, t) = P 7 (y, 0, t) = P 8 (y, 0, t) = P 9 (y, 0, t) sin e 1 c 1 (e 1 ) = sin e 2 c 2 (e 2 ) = sin e 3 c 3 (e 3 ) = sin e 4 c 4 (e 4 ) = sin e 5 c 5 (e 5 ) From equations (12) to (17), we have = sine 6 c 6 (e 6 ) = sin e 7 c 7 (e 7 ) = sin e 8 c 8 (e 8 ) = sin e 9 c 9 (e 9 ) (18) a 1 B 1 + a 2 B 2 + a 3 B 3 + a 4 B 4 + a 5 B 5 + a 6 B 6 + a 7 B 7 + a 8 B 8 + a 9 B 9 = 0, b 1 B 1 + b 2 B 2 + b 3 B 3 + b 4 B 4 + b 5 B 5 + b 6 B 6 + b 7 B 7 + b 8 B 8 + b 9 B 9 = 0, d 1 B 1 + d 2 B 2 + d 3 B 3 + d 4 B 4 + d 5 B 5 + d 6 B 6 + d 7 B 7 + d 8 B 8 + d 9 B 9 = 0, f 1 B 1 + f 2 B 2 + f 3 B 3 + f 4 B 4 + f 5 B 5 + f 6 B 6 + f 7 B 7 + f 8 B 8 + f 9 B 9 = 0, g 1 B 1 + g 2 B 2 + g 3 B 3 + g 4 B 4 + g 5 B 5 + g 6 B 6 + g 7 B 7 + g 8 B 8 + g 9 B 9 = 0, h 1 B 1 + h 2 B 2 + h 3 B 3 + h 4 B 4 + h 5 B 5 + h 6 B 6 + h 7 B 7 + h 8 B 8 + h 9 B 9 = 0,

Propagation of waves at an imperfectly bonded interface... 85 where the expressions for a l, b l, d l, f l, g l and h l are given in Appendix. The system of equations (18) can be reduced for different boundaries by taking: (i) (ii) K n 0, K t, K c, for Normal stiffness boundary, K n, K t 0, K c, for Transverse stiffness boundary, (iii) K n, K t, K c 0, for Thermal contact conductance, (iv) K n, K t 0, K c, for Slip boundary, (v) K n, K t, K c, for Welded contact. 5. Particular cases (a) If we take, c 24 = c 34 = c 24 = c 34 = 0, the system of equations (18) reduces for imperfectly bounded interface between two dissimilar orthotropic thermoelastic half-spaces as discussed by Kumar and Singh [12], where we have found minor but significant error in the coefficients of their equations (28) to (33). (b) If we take, K 1 = K 3 = K 1 = K 3 = 0, β 1 = β 3 = β 1 = β 3 = 0, the system of equations (18) reduces to for imperfectly bounded interface between two dissimilar monoclinic elastic half-spaces. If furthermore, we take K n, K t, the above results agree with Khurana and Singh [18]. (c) If we take, c 23 = c 33 2c 44 c 23 = c 33 2c 44, c 24 = c 34 = c 24 = c 34 = 0, the system of equations (18) reduces for imperfectly bonded interface between two dissimilar transversely isotropic thermoelastic half-spaces. (d) If we take, c 22 = c 33 = λ + 2µ, c 13 = c 23 = c 12 = λ, c 44 = c 55 = c 66 = µ, c 14 = c 24 = c 34 = c 56 = 0, β 2 = β 3 = β = (3λ + 2µ)α t, α 2 = α 3 = α t, K 2 = K 3 = K, c 22 = c 23 = λ + 2µ, c 13 = c 23 = c 12 = λ, c 44 = c 55 = c 66 = µ, c 14 = c 24 = c 34 = c 56 = 0, β 2 = β 3 = β = (3λ + 2µ )α t, α 2 = α 3 = α t, K 2 = K 3 = K, the system of equation (18) reduces for imperfectly bondedinterface between two dissimilar isotropic thermoelastic half-spaces. If we take further K n, K t, K c the above case agrees with Sinha and Elsibai [19, 20]. 6. Procedure for computing reflected and refracted angles The reflection coefficients depend on the velocities c i (e i ), i = 1, 2, 3, 4, 5, 6, which are the functions of reflected and refracted angles. They are

86 Joginder Singh, Baljeet Singh, Praveen Ailawalia supposed to be known for an incident qp wave, e 1 and, therefore c 1 (e 1 ). One has to compute e 4, e 5 and e 6 for given e 1. The velocities c 4 (e 4 ), c 5 (e 5 ) and c 6 (e 6 ) can then be computed from explicit algebraic formulae. The procedure is given below for computing e 4, e 5 and e 6 for given e 1 in the case of incident qp wave, for given e 2 in the case of incident qt wave, and for given e 3 in the case of incident qsv wave. Equation (11) can be written as (19) A 0 c 6 + A 1 c 4 + A 2 c 2 + A 3 = 0. Here, A 0 = ρ 3 τ, A 1 = ρ 2 [Uτ + Zτ + W C e + εη ( β2 p 2 3 + p2 2) ], [ (UZ A 2 = ρ V 2 ) τ + UW C e A 3 = W C e ( UZ V 2 ), where, + ZW ] + εη ( β2 Up 2 3 C + ) Zp2 2 2V η βp2 p 3, e U = c 22 p 2 2 + c 44p 2 3 + 2c 24p 2 p 3, Z = c 44 p 2 2 + c 33p 2 3 + 2c 34p 2 p 3, V = c 44 p 2 2 + c 34p 2 3 + (c 23 + c 24 )p 2 p 3, W = K 3 p 2 3 + K 2p 2 2. We define dimensionless apparent velocity c by the relation: (20) c = c a /β = c/p 2 β, where, β = (c 44 /ρ) 1 2, ρc 2 = p 2 2c 44 c 2. With the help of relation (20), the equation (19) becomes [ (21) c 6 (Ū + Z)τ + W ] c ( β 2 p 2 4 + 1) τ + εη C e c 2 44 [ + (Ū Z W V 2 )τ + (Ū + Z) + εη C e c 2 ( β 2 Ūp 2 + Z) 2εη β V 44 c 2 p 44 ] c 2 τ W C e (Ū Z V 2 ) 1 τ = 0.

Propagation of waves at an imperfectly bonded interface... 87 Here, (22) Ū = c 22 + p 2 + 2 c 24 p, Z = 1 + c 33 p 2 + 2 c 34 p, V = c 24 + c 34 p 2 + ( c 23 + 1)p, W = K3 p 2 + K 2, where, c ij = c ij c 44, K3 = K 3 c 44, K2 = K 2 c 44, p = p 3 p 2. For incident qp waves, p = cot e 1 ; for incident qt waves, p = cot e 2 ; for incident qsv waves, p = cot e 3 ; for reflected qp waves, p = cot e 4 ; for reflected qt waves, p = cot e 5 ; for reflected qsv waves, p = cot e 6. For a given p, equation (21) can be solved for c 3, the three roots corresponding to qp, qt and qsv waves. However, for a given c, equation (21) is a six degree equation in p, corresponding to incident qp, qt and qsv, reflected qp, qt and qsv. The positive roots corresponding to the reflected waves and the negative roots corresponding to the incident waves. On inserting the expressions for Ū, Z, V and W from equation (22) into equation (21), we obtain the six degree equation in p as: (23) G 0 p 6 + G 1 p 5 + G 2 p 4 + G 3 p 3 + G 4 p 2 + G 5 p + G 6 = 0. The expressions for G 0, G 1, G 2, G 3, G 4, G 5 and G 6 are given in Appendix. If we neglect thermal effects, the equation (23) reduces to the equation (48) obtained by Singh and Khurana [5]. Now, we define q = 1, The equation p (23) transforms into (24) G 6 q 6 + G 5 q 5 + G 4 q 4 + G 3 q 3 + G 2 q 2 + G 1 q + G 0 = 0. Equation will possess three positive roots for angles of incidence, for which all three reflected qp, qsv and qt waves exist, the smaller positive root (say q 6 ) corresponding to reflected qt and the root (q 5 ) corresponding to reflected qsv and the larger positive root (q 4 ) corresponding to reflected qp. Further, e 4 = tan 1 (q 4 ), e 5 = tan 1 (q 5 ), e 6 = tan 1 (q 6 ). For an isotropic thermoelastic medium, if we take: c 14 = c 24 = c 34 = 0, c 11 = c 22 = c 33 = λ + 2µ, c 44 = c 55 = c 66 = µ, c 12 = c 13 = c 23 = λ, K2 = K 3 = K, β 2 = β 3 = β

88 Joginder Singh, Baljeet Singh, Praveen Ailawalia Then, G 0 = K C e τ γ, G 1 = 0, { G 2 = γ + K εη (1 + γ) + C e τ µτ { G 4 = 1 + γ + K C e τ + εη µτ G 5 = 0, } c 2 + } c 4 + ( G 6 = c 6 1 + γ + K C e τ + εη ) ( µτ c 4 + γ + K C e τ ( 3γ), G 3 = 0, { 2γ + 2 K 2εη C e τ (1 + γ) + µτ K ) εη C e τ (1 + γ) + µτ c 2 } c 2 3 K C e τ γ, K C e τ γ, and, the equation (23) reduces to (25) sγ ( p 2 c 2 + 1 ) ( p 2 α 1 c 2 γ + 1 )(p 2 β 1 c 2 ) γ + 1 = 0, where α 1 = r + r 2 4sγ 2s, β 1 = r r 2 4sγ, r = γ + K 2s In this case, the Snell s law becomes: C e τ + εη K µτ, s = C e τ. (26) sin e 1 v qp = sin e 2 v qsv = sin e 3 v qt. Therefore, the roots of equation (25) are given by p 2 = c 2 1 = cot 2 e 3 corresponding to qsv waves, p 2 = c2 γ α 1 1 = cot 2 e 1, corresponding to qp wave and p 2 = c2 γ β 1 1 = cot 2 e 2 corresponding to qt wave. Thus, we may choose q = 1 p : (27) q 1 = tan e 1, q 2 = tan e 2, q 3 = tan e 3, q 4 = tan e 1, q 5 = tan e 2, q 6 = tan e 3,... as the six roots of the equation (24) for an isotropic medium. This choice will act as a guiding factor in computing the angles of reflection of qp, qt and qsv

Propagation of waves at an imperfectly bonded interface... 89 waves in a monoclinic thermoelastic medium. For an orthotropic medium, it can be shown that g 1 = g 3 = g 5 = 0 in equation (24). Therefore, equation (24) reduces to a cubic equation in q 2. Thus we may choose q 1 = q 4, q 2 = q 5, q 3 = q 6, so that the angle of reflection of qp, qt and qsv waves is equal to the angle of incidence of qp, qt and qsv waves, respectively. This is not true for a monoclinic thermoelastic medium. A similar analysis can be derived to compute the angles of refraction of qp, qsv and qt waves in the upper monoclinic half-space for given incident angle. 7. Conclusion A problem of reflection and transmission is considered at an imperfectly bonded interface between two dissimilar thermoelastic solid half-spaces of monoclinic type, where we have obtained the six relations between amplitudes of incident, reflected and transmitted waves. These relations are verified theoretically for some limiting cases. A procedure for computing the reflected and refracted angles for a given incident angle of a striking wave is explained in detail. It is concluded that the angle of particular reflected wave is not equal to the angle of incidence of that wave in a monoclinic medium. Numerical verification of present theoretical results can be attempted by following the procedure given in the previous section. REFERENCES [1] Keith, C. M., S. Crampin. Seismic Body Waves in Anisotropic Media: Reflection and Refraction at a Plane Interface. Geophys. J. R. Astr. Soc, 49 (1977), 181 208. [2] Chattopadhyay, A., S. Choudhary. The Reflection Phenomenon of P- waves in a Medium of Monoclinic Type. Int. J. Engg. Sci, 33 (1995), 195 207. [3] Chattopadhyay, A., S. Saha, M. Chakraborty. The Reflection of SVwaves in a Monoclinic Medium. Indian J. pure Appl. Math., 27 (1996), 1029 1042. [4] Singh, S. J. Comments on The Reflection Phenomenon of P-waves in a Medium ofmonoclinic Type by ChattopadhyayandChoudhary. (Int. J. Engg. Sci., 33 (1995), 195 207), Int. J. Engg. Sci., 37 (1999), 407 410.

90 Joginder Singh, Baljeet Singh, Praveen Ailawalia [5] Singh, S., S. Khurana. Reflection of P and SV Waves at the Free Surface of a Monoclinic Elastic Half- Space, Proc. Indian Acad. Sci. (Earth Planet. Sci.), 111, 2002, 401 412. [6] ] Biot, M. Thermoelasticity and Irreversible Thermodynamics. J. Appl. Phys., 27 (1956), 249 253. [7] Lord, H. W., Y. Shulman. A Generalized Dynamical Theory of Thermoelasticity. J. Mech. Phys. Solids, 15 (1967), 299 309. [8] Dhaliwal, R. S., H. H. Sherief. Generalized Thermoelasticity for Anisotropic Media. Quart. Appl. Math., 33 (1980) 1 8. [9] Green, A. E., K. A. Lindsay. Thermoelasticity. J. Elasticity, 2 (1972), 1 7. [10] Sinha, A. N., S. B. Sinha. Reflection of Thermoelastic Waves at a Solid Halfspace With Thermal Relaxation. J. Phys. Earth, 22 (1974), 237 244. [11] Sharma, J. N. Reflection of Thermoelastic Waves from the Stress-free Insulated Boundary of an Anisotropic Half-space. Indian J. pure Appl. Math., 19 (1988), 294 304. [12] Kumar, R., M. Singh. Reflection/ Transmission of Plane Waves at an Imperfectly Bonded Interface of Two Orthotropic Generalized Thermoelastic Half- Spaces. Material Sci. Engg. A, 472 (2008), 83 96. [13] Singh, B. Dispersion Relations in a Generalized Monoclinic Thermoelastic Solid half-space. Journal of Technical Physics, 47 (2006), 119 129. [14] Singh, B. Reflection of Plane Waves at the Free Surface of a Monoclinic Thermoelastic Solid Half-space. European J. Mech. A/Solids, 29 (2010), 911 916. [15] Singh, B. Wave Propagation in an Initially Stressed Transversely Isotropic Thermoelastic Solid Half-space. Applied Mathematics and Computation, 217 (2010), 705 715. [16] Rokhlin, S.I., D. Marom. Study of Adhesive Bonds using Low-frequency Obliquely Incident Ultrasonic Waves. J. Acoust. Soc. Am., 80 (1986), 585 590. [17] Lavrentyev, A.I., S.I. Rokhlin. Ultrasonic Spectroscopy of Imperfect ContactInterfaces Betweena Layerand TwoSolids. J. Acoust. Soc. Am., 103 (1998), 657 664. [18] Singh, S. J., S. Khurana. Reflection and Transmission of P and SV Waves atthe Interface BetweenTwo Monoclinic Elastic Half-spaces, Proc. Natl. Acad. Sci. India 71(A), 2001, 305 319. [19] Sinha, S. B., K. A. Elsibai. Reflection of Thermoelastic Waves at a Solid Half-space with two Relaxation Times. J. Thermal Stresses, 19 (1996), 749 762. [20] Sinha, S. B., K. A. Elsibai. Reflection and Refraction of Thermoelastic Waves atan Interface of Two Semi-infinite Media with two RelaxationTimes. J. Thermal Stresses, 20 (1997), 129-145.

Propagation of waves at an imperfectly bonded interface... 91 Appendix The expressions for a l, b l, d l, f l, g l and h l are a l = ik n, (l = 1,...,6) [ ] a l = ik n c ω sine l 23F l + c ω cose l ω cose l 33 + c 34 F l c ω sine l 34 + iβ c l c l c l c 3τ 1F l, l b l = ik t F l, (l = 1,...,6) [ ] b l = ik t F l c ω sine l 24 F l + c ω cose l 34 + c ω cose l 44 F l c ω sine l 44, c l c l c l c l [ d l = ik c Fl, (l = 1,...,6), d l = ik c Fl + K 3 ω cos e l c l ω sine l ω cose l f l = c 23 F l + c 33 c l c l ω sine l ω cose l f l = c 23 F l c 33 c l c l F l ω cose l ω sine l + c 34 F l c 34 + iβ 3 τ 1 Fl c l c, l ω cose l ω sine l c 34 F l c 34 + iβ 3 τ 1 Fl, c l c l (l = 7,...,9) (l = 7,...,9) ], (l = 7,...,9) (l = 1,...,3) (l = 4,...,6) f l = c ω sine l 23 F l c ω cose l 33 c ω cose l 34 F l + c ω sine l 34 iβ 3 c l c l c l c τ 1 F l, l (l = 7,...,9) ω sine l ω cose l g l = c 24 F l + c 34 c l c l ω sine l ω cose l g l = c 24 F l c 34 c l c l ω cose l ω sine l + c 44 F l c 44, (l = 1,...,3) c l c l ω cose l ω sine l c 44 F l c 44, (l = 4,...,6) c l c l g l = c ω sine l 24 F l c ω cose l 34 c ω cose l 44 F l + c ω sine l 44, (l = 7,...,9) c l c l c l c l h l = K 3 F l ω cose l, (l = 1,...,3), h l = K 3 F ω cose l l, (l = 4,...,6) c l c l h l = K 3F l ω cos e l c l, (l = 7,...,9)

92 Joginder Singh, Baljeet Singh, Praveen Ailawalia τ 1 = (1 + iωt 1), τ 1 = (1 + iωt 1 ) The expressions for G 0, G 1, G 2, G 3, G 4, G 5 and G 6 are G 0 = K 3 ) 2 ( c C e τ 34 c 33, G1 = 2 K 3 C e τ ( c 34 c 23 c 24 c 33 ) ( G 2 = c 33 c 2 34 + K 3 C e τ (1 + c 33) + εη ) c 44 τ β 2 c 2 K 3 ( c22 c C e τ 33 2 c 24 c 34 c 2 23 2 c ) K2 23 ( c33 C e τ c 2 ) 34 [{ G 3 = 2 c 24 c 33 c 34 c 23 + K } 2 3 C e τ ( c εη β 24 + c 34 ) + c 44 τ c εη β 24 c 44 τ c 34 K ] 3 C e τ ( c 22 c 34 c 24 c 23 ) K 2 C e τ ( c 24 c 33 c 34 c 23 ) { G 4 = (1 + c 33 ) + K 3 C e τ + εη } c 44 τ β 2 G 5 = 2 { c 4 + c 22 c 33 + 1 + 2 c 24 c 34 (1 + c 23 ) 2 + K 3 C e τ ( c 22+1) + K 2 2 C e τ (1+ c εη β 33)+ c 44 τ + c 22 2 εη β c 44 τ (1+ c 23) + εη } c 44 τ c 33 c 2 + K 3 ) ( ) 2 K2 ( c C e τ 24 c 22 + C e τ (1 + c 23 ) 2 c 22 c 33 1 + 2 c 24 c 34 [ { ( c 24 + c 34 ) c 4 + c 22 c 34 c 24 c 23 + K 2 + K 2 G 6 = c 6 ] C e τ ( c 24 c 23 c 22 c 34 ) { (1 + c 22 ) + K 2 C e τ + εη c 44 τ + K 2 C e τ ( c2 24 c 22) C e τ ( c 24 + c 34 ) + εη c 2 } c 44 τ ( c 34 c 24 β) } { c 4 + ( c 22 c 2 24)+ K 2 C e τ ( c 22 + 1)+ εη } c 44 τ c 2 c 2