# Effects of Rotation and Gravity Field on Surface Waves in Fibre-Reinforced Thermoelastic Media under Four Theories

Authors: A. M. Abd-Alla; S. M. Abo-Dahab; A. Al-Mullise
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

Estimation is done to investigate the gravitational and rotational parameters effects on surface waves in fibre-reinforced thermoelastic media. The theory of generalized surface waves has been firstly developed and then it has been employed to investigate particular cases of waves, namely, Stoneley waves, Rayleigh waves, and Love waves. The analytical expressions for surface waves velocity and attenuation coefficient are obtained in the physical domain by using the harmonic vibrations and four thermoelastic theories. The wave velocity equations have been obtained in different cases. The numerical results are given for equation of coupled thermoelastic theory (C-T), Lord-Shulman theory (L-S), Green-Lindsay theory (G-L), and the linearized (G-N) theory of type II. Comparison was made with the results obtained in the presence and absence of gravity, rotation, and parameters for fibre-reinforced of the material media. The results obtained are displayed by graphs to clear the phenomena physical meaning. The results indicate that the effect of gravity, rotation, relaxation times, and parameters of fibre-reinforced of the material medium is very pronounced.


---

## Body

## 1. Introduction

A reinforced concrete member will be designed for all conditions of stress that may occur accordance with the principle of mechanics. Fibre-reinforced composites are used in a variety of structures due to their low weight and high strength. The characteristic property of a reinforced composite is that its components act together as single anisotropic units as long as they remain in the elastic condition. Investigation of interaction between magnetic field, stress, and strain in a thermoelastic solid is very important due to its many applications in the field of geophysics, plasma physics, and related topics, especially in the nuclear field, where the extremely high temperature and temperature gradients, as well as the magnetic fields originating inside nuclear reactors, influence their design. Recently, more attention has been studied the dynamical problem of propagation of surface waves in a homogeneous and non-homogeneous elastic and thermoplastic media because of its utilitarian aspects in earthquake, engineering, and seismology on account of the occurrence of nonhomogeneities in the Earth’s crust, as the Earth is made up of different layers. Abd-Alla et al. [1] investigated propagation of Rayleigh waves in generalized magnetothermoelastic orthotropic material under initial stress and gravity field. Thermal stresses in an infinite circular cylinder have been investigated by Abd-Alla et al. [2]. Abd-Alla and Mahmoud [3] studied the magnetothermoelastic problem in rotational nonhomogeneous orthotropic hollow cylinder under the hyperbolic heat conduction model. Abd-Alla and Ahmed [4] investigated the Rayleigh waves in an orthotropic thermoelastic medium under gravity field and initial stress. Rayleigh waves in a magnetoelastic half-space of orthotropic material under the influence of initial stress and gravity field have been investigated by Abd-Alla et al. [5]. Elnaggar and Abd-Alla [6] studied Rayleigh waves in magneto-thermomicroelastic half-space under initial stress. Abd-Alla and Ahmed [7] discussed Rayleigh waves in an orthotropic thermoelastic medium under gravity field and initial stress. Effect of rotation and initial stress on an infinite generalized magnetothermoelastic diffusion body with a spherical cavity has been investigated by Abd-Alla and Abo-Dahab [8]. Wu and Chai [9] studied propagation of surface waves in anisotropic solids: theoretical calculation and experiment. Wu and Liu [10] investigated the measurement of anisotropic elastic constants of fibre-reinforced composite plate by using ultrasonic bulk wave and laser generated Lamb wave. The group velocity variation of Lamb wave in Fibre-reinforced composite plate studied by Rhee et al. [11]. Fu and Zhang [12] investigated the continuum-mechanical modeling of kink-band formation in fibre-reinforced composites. Espinosa et al. [13] discussed the modeling impact induced delamination of woven Fibre-reinforced composites with contact/cohesive laws. Wave propagation in materials reinforced with bi-directional Fibres presented by Weitsman and Benveniste [14]. Weitsman [15] introduced the wave propagation and energy scattering in materials reinforced by inextensible Fibres. Dai and Wang [16] considered the stress wave propagation in piezoelectric Fibre-reinforced laminated composites subjected to thermal shock. Ohyoshi [17] studied the propagation of Rayleigh waves along an obliquely cut surface in a directional Fibre-reinforced composite. Rogerson [18] investigated the penetration of impact waves in a six-ply fibre composite laminate. Weitsrian [19] studied the reflection of harmonic waves in fibre-reinforced materials. Huang et al. [20] investigated the effect of fibre-matrix interphase on wave propagation along, and scattering from, multilayered fibres in composites. Fu and Zhang [21] discussed the continuum-mechanical modeling of kink-band formation in fibre-reinforced composites. B. Singh and S. J. Singh [22] investigated the reflection of plane waves at the free surface of a fibre-reinforced elastic half-space. Sengupta and Nath [23] studied the surface waves in fibre-reinforced anisotropic elastic media. Samal and Chattaraj [24] studied the surface wave propagation in Fibre-reinforced anisotropic elastic layer between liquid saturated porous half-space and uniform liquid layer. Chatopadhyay et al. [25] investigated the reflection of quasi-P and quasi-SV waves at the free and rigid boundaries of a fibre-reinforced medium. Singh [26] discussed the wave propagation in an incompressible transversely isotropic fibre-reinforced elastic media. Singh [27] studied the wave propagation in thermally conducting linear fibre-reinforced composite materials. Abd-Alla et al. [28] studied the thermal stresses in a nonhomogeneous orthotropic elastic multilayered cylinder.Recently, Othman and Lotfy [29] investigated the effect of magnetic field and rotation of the 2D problem of a fibre-reinforced thermoelastic under three theories with influence of gravity. Abouelregal and Abo-Dahab [30] discuss the dual phase lag model on magnetothermoelasticity infinite nonhomogeneous solid having a spherical cavity. Abd-Alla and Abo-Dahab [31] investigated the effect of magnetic field on poroelastic bone model for internal remodeling. Abd-Alla et al. [32, 33] studied the Rayleigh waves in generalized magneto thermoviscoelastic granular medium under the influence of rotation, gravity field, and initial stress and magnetothermoelastic half-space of a homogeneous orthotropic material under the effect of the rotation, initial stress and gravity field. Abo-Dahab [34] investigated propagation of Stoneley waves in magnetothermoelastic materials with voids and two relaxation times.The present investigation is devoted to study the propagation of surface waves in fibre-reinforced anisotropic thermoelastic medium subjected to gravity field and relaxation times leading to particular cases such as Rayleigh waves, Love waves and Stoneley waves. The surface wave velocity and attenuation coefficient are obtained in the physical domain by using the harmonic vibrations. The effects of the gravity, relaxation times, and parameters for fibre-reinforced of the material medium on surface waves are studied simultaneously. The analytical expressions for surface waves velocity and attenuation coefficient are represented graphically. Numerical results for the surface waves velocity and attenuation coefficient are given and illustrated graphically in the presence and absence of the gravity field and rotation.


## 2. Formulation of the Problem

We consider homogeneous thermally conducting transversely two fibre-reinforced media. LetM1 and M2 be two fibre-reinforced elastic thermoelastic semi-infinite solid media. They are perfectly welded in contact to prevent any relative motion or sliding before and after the disturbances and that the continuity of displacement, stress, and so for the hold good across the common boundary surface. Further the mechanical properties of M1 are different from those of M2. These media extend to an infinite great distance from the origin and are separated by a plane horizontal boundary and M2 is to be taken above M1. Let Oxyz be a set of orthogonal Cartesian coordinates, let O be any point of the plane boundary, and let Oy points vertically downward to the medium M1. We consider the possibility of a type of wave travelling in the direction Ox in such a manner that the disturbance is largely confined to the neighborhood of the boundary and, at any instant, all particles in any line parallel to y-axis have equal displacements. These two assumptions conclude that the wave is a surface wave and all partial derivatives with respect to z are zero. Further let us assume that u, v are the components of displacements at any point (x,y,z) at any time t. It is also assumed that gravitational field produces a hydrostatic initial stress that is produced by a slow process of creep where the shearing stresses tend to become smaller or vanish after a long period of time.The equilibrium equation of the initial stress is in the form(1)∂τ∂x=0,∂τ∂y+ρg=0.The equations and constitutive relations for such medium in the absence of body forces and heat sources are(2)(λ+2α+4μL-2μT+β)∂2u∂x2+(λ+α+μL)∂2v∂x∂y(λ(λ+μL∂2u∂y2-β1(1+ϑ0∂∂t)∂T∂x-ρg∂v∂x=ρ(∂2u∂t2-2Ω∂v∂t-Ω2u),(3)(λ+2μT)∂2v∂y2+(λ+α+μL)∂2u∂x∂y+μL∂2v∂x2(λ(λ-β2(1+ϑ0∂∂t)∂T∂y+ρg∂u∂x=ρ(∂2v∂t2+2Ω∂u∂t-Ω2v),(4)k1*(n*+t1∂∂t)∂2T∂x2+k2*(n*+t1∂∂t)∂2T∂y2=ρC*(n1∂∂t+τ0∂2∂t2)T+T0(n1∂∂t+τ0n0∂2∂t2)(λ×(β1∂u∂x+β2∂v∂y).


## 3. Special Cases of Thermoelasticity Theory

We get that the forms show the behavior under gravity field; the theories are as follows.(1) The equation of the coupled thermoelasticity (C-T theory) is obtained when(5)n*=n1=1,t1=τ0=θ0=0,τ0=0.(2) Consider Lord-Shulman theory.For the Lord-Shulman (L-S) theory(6)n*=n1=n0=1,t1=θ0=0,τ0>0.(3) Consider the Green-Lindsay theory.For the Green-Lindsay (G-L) theory(7)n*=n1=1,n0=t1=0,θ0≥τ0>0.(4) Consider the Green-Nagdhi theory.For the Green-Nagdhi theory (G-N) of type 2,(8)n*>0,n1=0,n0=1,t1=ϑ0=0,τ0=1.For plane strain deformation in the x-y plane, the components of stress take the form
(9)τ11=(λ+2α+4μL-2μT+β)∂u∂x+(λ+α)∂v∂y-β1(1+ϑ0∂∂t)T,(10)τ12=μL∂v∂x+μL∂u∂y,(11)τ21=μL∂v∂x+μL∂u∂y,(12)τ22=(λ+α)∂u∂x+(λ+2μT)∂v∂y-β2(1+ϑ0∂∂t)T,
where
(13)β1=(2λ+3α+4μL-2μT+β)α1+(λ+α)α2,β1=(2λ+α)α1+(λ+2μT)α2.
By Helmholtz’s theorem [35], the displacement vector u can be written in the form(14a)u→=gradφ+curlψ→,
where the scalar φ and the vector ψ→ represent irrotational and rotational parts of the displacement u→. It is possible to take only one component of the vector ψ→ to be nonzero, as
(14b)ψ→=(0,0,ψ).
From (14a) and (14b) we obtain
(14c)u=∂φ∂x+∂Ψ∂y,v=∂φ∂y-∂Ψ∂x.Substituting from (14c) into (2) and (3) we get
(15)(λ+2α+4μL-2μT+β)ρ∂2φ∂x2+(λ+α+2μL)ρ∂2φ∂y2(λ-β1(1+ϑ0∂∂t)T+g∂Ψ∂x=∂2φ∂t2+2Ω∂Ψ∂t-Ω2φ,(16)(α+3μL-2μT+β)ρ∂2Ψ∂x2+μLρ∂2Ψ∂y2-g∂φ∂x=∂2Ψ∂t2-2Ω∂φ∂t-Ω2Ψ,(17)(λ+α+2μL)ρ(∂2φ∂x2)+(λ+2μT)ρ(∂2φ∂y2)-β2ρ(1+ϑ0∂∂t)T+g∂Ψ∂x=∂2φ∂t2+2Ω∂Ψ∂t-Ω2φ,(18)-μLρ(∂2Ψ∂x2)+(α+μL-2μT)ρ(∂2Ψ∂y2)+g∂φ∂x=-∂2Ψ∂t2+2Ω∂φ∂t+Ω2Ψ.
From (16) and (18) we get
(19)∂2Ψ∂x2+f1∂2Ψ∂y2=0,
where
(20)f1=(α+2μL-2μT)(α+2μL-2μT+β).Also from (15) and (17) we get
(21)∂2φ∂x2+f2∂2φ∂y2-f3(1+ϑ0∂∂t)T+f4∂Ψ∂x=f5(∂2φ∂t2+2Ω∂Ψ∂t-Ω2φ),
where
(22)f2=(2λ+α+2μL-2μT)(2λ+3α+6μL-2μT+β),f3=(β2+β1)(2λ+3α+6μL-2μT+β),f4=2g2λ+3α+6μL-2μT+β,f5=2(2λ+3α+6μL-2μT+β).The solution of (19) and (21) has the form
(23)φ=F(y)eiω(x-ct),Ψ=G(y)eiω(x-ct),T=H(y)eiω(x-ct).
Substituting from (23) in to (19) we get
(24)G′′-m2G=0.
The solution for (24) as the form
(25)G=C1emy+C2e-my.
Equation (24) must have exponential solutions in order that G willdescribe surface waves; they must become vanishingly small as y→∞:
(26)G=C2e-my,
where m2=ω2/f1.Substituting from (23) in to (21) we get
(27)H=f6F′′+f7F+f8C2e-my,
where
(28)f6=f2f3(1-iωcϑ0),f7=ω2(c2f5-1)+f5Ω2f3(1-iωcϑ0),f8=iω(f4-2cf5Ω)f3(1-iωcϑ0).From (26) and (27), we can obtain the temperature T and the potential function Ψ as follows:
(29)Ψ=C2e-myeiω(x-ct),T=[f6F′′+f7F+f8C2e-my]eiω(x-ct).
Substituting from (29) and (14c) into (4) for T, Ψ, u, and v we get
(30)F′′′′+aF′′+bF=dC2e-my,
where
(31)a=f10f9,b=f11f9,d=f12f9f9=(n*-iωct1)k2*f6,f10=-ω2k1*f6(n*-iωct1)+(n*-iωct1)k2*f7-ρC*(-iωcn1-ω2c2τ0)f6-T0β2(-iωcn1-ω2c2τ0n0),f11=-ω2k1*(n*-iωct1)f7-ρC*(-iωcn1-ω2c2τ0)f7+ω2T0β1(-iωcn1-ω2c2τ0n0),f12=-ω2k1*(n*-iωct1)f8+k2*(n*-iωct1)m2f8-ρC*(-iωcn1-ω2c2τ0)f8-T0(-iωcn1-ω2c2τ0n0)×(-iωmβ1+iωmβ2).
The solution of (30) is
(32)F=C4e-f14y+C6e-f16y+MC2e-my,
where
(33)f14=-aa2-4b2,f16=-a2+12a2-4b,M=dm4+am2-bm.After finding the values ofF we get the following values:
(34)φ=[C4e-f14y+C6e-f16y+MC2e-my]eiω(x-ct),(35)Ψ=C2e-myeiω(x-ct),(36)T={(f6f142+f7)C4e-f14y+(f6f162+f7)C6e-f′16y+[(m2f6+f7)M+f8]C2e-my{(f6f142+f7)C4e-f14y+(f6f162+f7)C6e-f′16y}eiω(x-ct).
Substituting from (34), (35), and (36) into (14c), we get the components of displacement as
(37)u={iωC4e-f14y+iωC6e-f16yj+(iωM-m)C2e-my{iωC4e-f14y+iωC6e-f16y{iωC4e-f14y+iωC6e-f16y}eiω(x-ct),(38a)v={-f14C4e-f14y-f16C6e-f16yh-(mM+iω)C2e-my{-f14C4e-f14y-f16C6e-f16y}eiω(x-ct).
Substituting from (36), (37), and (38a) in to (9)–(12) we obtain
(38b)τ11={λ(f142-ω2)+α(f142-2ω2)-4μLω2+2μTω2-βω2-β1(1-iωϑ0)(f6+f7)}×C4e-f14yeiω(x-ct)+{λ(f162-ω2)+α(f162-2ω2)×C4e-f14yeiω(x-ct)+++-4μLω2+2μTω2-βω2×C4e-f14yeiω(x-ct)+++-β1(1-iωϑ0)(f6f162+f7)}×C6e-f16yeiω(x-ct)+{λ(m2-ω2)+α[(m2-2ω2)M-iωm]-(4μL-2μT+β)(ω2M+iωm)-β1(1-iωϑ0)[(m2f6+f7)M+f8]}×C2e-myeiω(x-ct),(38c)τ22={λ(f142-ω2)-αω2+2μTf142-β2(1-iωcϑ0)×(f6f142+f7)}C4e-f14yeiω(x-ct)+{λ(f162-ω2)-αω2+2μTf162-β2(1-iωcϑ0)(f6f162+f7)}×C6e-f16yeiω(x-ct)+{λ(m2-ω2)M-α(ω2M+iωm)+2μTm(mM+iω)-β2(1-iωcϑ0)[(m2f6+f7)M+f8]}×C2e-myeiω(x-ct)(38d)τ12={[-2iωμLf14]e-f14yC4+[-2iωμLf16]C6e-f16y+[μL(-2iωmM+ω2-m2)]C2e-my}eiω(x-ct).


## 4. Boundary Conditions

The stresses components, the displacement components, and temperature at the boundary surface between the mediaM1 and M2 must be continuous at all times and positions.Consider the following:(39)τ22=τ22′aty=0,τ12=τ12′aty=0,u=u′aty=0,v=v′aty=0,∂T∂y+hT=∂T′∂y+h′T′aty=0.Finally, eliminating the constantsC2, C4, C6, C2′, C4′, and C6′ from (10), (12), (36), (37), and (38a), (38b), (38c), and (38d), we get
(40)det(aij)=0,i,j=1,2,3,
where
(41)a11=λ(f142-ω2)-αω2+2μTf142-β2(1-iωcϑ0)(f6f142+f7),a12=λ(f162-ω2)-αω2+2μTf162-β2(1-iωcϑ0)(f6f162+f7),a13=λ(m2-ω2)M-α(ω2M+iωm)+2μTm(mM+iω)-β2(1-iωcϑ0)×[(m2f6+f7)M+f8],a14=-[λ′(f14′2-ω2)-α′ω2+2μ′Tf14′2-β2′(1-iωcϑ0)(f6′f14′2+f7′)],a15=[λ′(f16′2-ω2)-α′ω2+2μ′Tf16′2-β2′(1-iωcϑ0)(f6′f16′2+f7′)],a16=-[λ′(m′2-ω2)M′-α′(ω2M′+iωm′)+2μ′Tm′(m′M′+iω)-β2′(1-iωcϑ0)×[(m′2f6′+f7′)M′+f8′][λ′(m′2-ω2)M′-α′(ω2M′+iωm′)[λ′(m′2-ω2)M′-α′(ω2M′+iωm′)],a21=-2iωμLf14,a22=-2iωμLf16,a23=μL(-2iωmM+ω2-m2),a24=2iωμL′f14′,a25=2iωμL′f16′,a26=-[μL′(-2iωm′M′+ω2-m′2)],a31=iω,a32=iω,a33=iωM-m,a34=-iω,a35=-iω,a36=-(iωM′-m′),a41=-f14,a42=-f16,a43=-(mM+iω),a44=f14′,a45=f16′,a46=m′M+iω,a51=(f6f142+f7)(h-f14),a52=(f6f162+f7)(h-f16),a53=[(m2f6+f7)M+f8](1-m),a54=-[(f6′f14′2+f7′)(h′-f14′)],a55=-[(f6′f16′′2+f7′)(h′-f16′)],a56=-[(m′2f6′+f7′)M′+f8′](1-m′),a61=(f6f142+f7),a62=(f6f162+f7),a63=[(m2f6+f7)M+f8],a64=-(f6′f14′2+f7′),a65=-(f6′f16′2+f7′),a66=-[(m′2f6′+f7′)M′+f8′].


## 5. Stoneley Waves

It is the generalized form of Rayleigh waves in which we assume that the waves are propagated along the common boundary of two semi-infinite mediaM1 and M2. Therefore (40) determines the wave velocity equation for Stoneley waves in anisotropic fibre-reinforced solid thermoelastic media under the influence of gravity. Clearly, from (40), it follows that wave velocity of the Stoneley waves depends upon the parameters for fibre-reinforced of the material medium, gravity, rotation, and the densities of both mediums, since the wave velocity equation (40) for Stoneley waves under the present circumstances depends on the particular value of ω and creates a dispersion of a general wave form. Further, (40), of course, is in complete agreement with the corresponding classical result, when the effect of gravity, rotation, and parameters for fibre-reinforcement is ignored.


## 6. Rayleigh Waves

To investigate the possibility of Rayleigh waves in anisotropic fibre-reinforced thermoelastic media, we replace mediumM2 by vacuum, in the preceding problem, since the boundary(y=0)is adjacent to vacuum. So the boundary condition in this case may be expressed as
(42)τ12=0aty=0,v=0aty=0,∂T∂y+hT=0aty=0.
Eliminating the constants C2, C4, and C6from (10), (12), (34), (35), (36), and (37), we get
(43)det(bi,j)=0,i,j=1,2,3,
where
(44)b11=-2iωμLf14,b12=-2iωμLf16,b13=μL(-2iωmM+ω2+m2),b21=f14,b22=f16,b23=(mM+iω),b31=(h-f14)(f6f142+f7),b32=(h-f16)(f6f162+f7),b33=(h-m)[(m2f6+f7)M+f8].


## 7. Love Waves

To investigate the possibility of love waves in a fibre-reinforced thermoelastic solid media, we replace mediumM2 that is obtained by two horizontal plane surfaces at a distance H-apart, while M1 remains infinite. For medium M1, the displacement component ν remains the same as in general case given by (38a), (38b), (38c), and (38d).The boundary condition of Love wave is as follows:(45)τ12=τ12′,aty=0,v=v′,aty=0,T=T′,aty=0,τ12′=0,aty=-H,v′=0aty=-H,T′=0aty=-H.Eliminating the constantsC2, C4, C6, C2′, C4′, and C6′ from (10), (12), (34), (35), (36), and (37), we get
(46)det(ci,j)=0,i,j=1,2,3,(47)c11=-2iωμLf14,c12=-2iωμLf16,c13=μL(-2iωmM+ω2+m2),c14=2iωμL′f14′,c15=2iωμL′f16′,c16=-μL′(-2iωm′M′+ω2+m′2),c21=f14,c22=f16,c23=mM+iω,c24=-f14′c25=-f16′,c26=-(m′M′+iω),c31=f6f142+f7,c32=f6f162+f7,c33=(m2f6+f7)M+f8,c34=-(f6′f14′2+f7′),c35=-(f6′f16′2+f′7),c36=-[(m′2f6′+f7′)M′+f8′],c41=c42=c43=0,c44=2iωμL′f14′ef14′H,c45=2iωμL′f16′ef16′H,c46=-μL′(-2iωm′M′+ω2+m′2)em′H,a51=0,c52=0,c53=0,c54=f14′ef14′H,c55=f16′ef16′H,c56=(m′M′+iω)em′H,c61=0,c62=0,c63=0,c64=(f6′f14′2+f7′)ef14′H,c65=(f6′f16′2+f7′)ef16′H,c66=[(m′2f6′+f7′)M′+f8′]em′H.


## 8. Numerical Results and Discussion

With a view to illustrating the analytical procedure presented earlier, we now consider a numerical example for which computational results are given. The results depict the variation of Stoneley waves velocity and attenuation coefficient, Stoneley secular equation, Rayleigh secular equation, Rayleigh waves velocity and attenuation coefficient, Love waves velocity and attenuation coefficient, and Love secular equation in the context of four theories. To study the influence of gravity field and rotation with reinforcement on wave propagation, we now present some numerical results for the physical constants:(48)ρ=2660kg/m3,λ=5.65×1010Nm-2,μT=2.46×1010Nm-2,μL=5.66×1010Nm-2,C*=0.787×103JKg-1deg-1,α=-1.28×1010Nm-2,β=220.90×1010Nm-2,K1*=0.0921×103Jm-1deg-1s-1,K2*=0.0963×103Jm-1deg-1s-1,α1=0.017×104deg-1,α2=0.015×104deg-1,T0=2930K.The output is plotted in Figures1–30.Figure 1
Variation of the Stoneley waves velocity and attenuation coefficient (CT) model with respect to the gravity.Figure 2
Variation of the Stoneley waves velocity and attenuation coefficient (LS) model with respect to the gravity.Figure 3
Variation of the Stoneley waves velocity and attenuation coefficient (GL) model with respect to the gravity.Figure 4
Variation of the Stoneley waves velocity and attenuation coefficient (GN type (II)) model with respect to the gravity.Figure 5
Variation of the Stoneley waves velocity and attenuation coefficient (CT) model with respect to the rotation.Figure 6
Variation of the Stoneley waves velocity and attenuation coefficient (LS) model with respect to the rotation.Figure 7
Variation of the Stoneley waves velocity and attenuation coefficient (GL) model with respect to the rotation.Figure 8
Variation of the Stoneley waves velocity and attenuation coefficient (GN type (II)) model with respect to the rotation.Variation of the Stoneley secular equation for the (CT, LS, GL, and GN (type (II)) models with respect to the phase velocity.
(a)
(b)
(c)
(d)Figure 10
Variation of the Stoneley waves velocity and attenuation coefficient (LS) model with respect to the phase velocity.Figure 11
Variation of the Stoneley waves velocity and attenuation coefficient (GL) model with respect to the phase velocity.Figure 12
Variation of the Stoneley waves velocity and attenuation coefficient (GN type (II)) model with respect to the phase velocity.Figure 13
Variation of the Stoneley waves velocity and attenuation coefficient (GN type (II)) model with respect to the phase velocity.Figure 14
Variation of the Rayleigh waves velocity and attenuation coefficient (CT) model with respect to the gravity.Figure 15
Variation of the Rayleigh waves velocity and attenuation coefficient (LS) model with respect to the gravity.Figure 16
Variation of the Rayleigh waves velocity and attenuation coefficient (GL) model with respect to the gravity.Figure 17
Variation of the Rayleigh waves velocity and attenuation coefficient (GN type (II)) model with respect to the gravity.Figure 18
Variation of the Rayleigh waves velocity and attenuation coefficient CT-model with respect to the rotation.Figure 19
Variation of the Rayleigh waves velocity and attenuation coefficient LS-model with respect to the rotation.Figure 20
Variation of the Rayleigh waves velocity and attenuation coefficient GL-model with respect to the rotation.Figure 21
Variation of the Rayleigh waves velocity and attenuation coefficient (GN type (II)) model with respect to the rotation.Variation of the Rayleigh secular equation for the (CT, LS, GL, and GN (type (II)) models with respect to the phase velocity.
(a)
(b)
(c)
(d)Figure 23
Variation of the Rayleigh waves velocity and attenuation coefficient (LS) model with respect to the phase velocity.Figure 24
Variation of the Rayleigh waves velocity and attenuation coefficient (GL) model with respect to the phase velocity.Figure 25
Variation of the Rayleigh waves velocity and attenuation coefficient (GN type (II)) model with respect to the phase velocity.Figure 26
Variation of the Rayleigh waves velocity and attenuation coefficient (GN type (II)) model with respect to the phase velocity.Variation of the Love secular equation for the (CT, LS, GL, and GN (type (II)) models with respect to the phase velocity.
(a)
(b)
(c)
(d)Figure 28
Variation of the Love waves velocity and attenuation coefficient (LS) model with respect to the phase velocity.Figure 29
Variation of the Love waves velocity and attenuation coefficient (GL) model with respect to the phase velocity.Figure 30
Variation of the Love waves velocity and attenuation coefficient (GN type (II)) model with respect to the phase velocity.Figures1, 2, 3, 4, 5, 6, 7, and 8 show the variation of the Stoneley waves and attenuation coefficient value with respect to gravity field g and rotation Ω, respectively, has oscillatory behavior in the whole range of gravity g and rotation for C-T theory, L-S theory, G-L theory, and G-N theory. These figures indicate that the medium along g and Ω undergoes compressive deformation due to the thermal shock. In both figureures, it is clear that the Stoneley wave velocity and attenuation coefficient have a nonzero value only in a bounded region of space.Figure9 shows the variation of the Stoneley secular equation with respect to phase velocity c, that has an oscillatory behavior of C-T theory, L-S theory, and G-L theory for different values of relaxation time τ0 in the whole range of the phase velocity c, while in G-N theory it decreases with increasing the constant n* and it increases with increasing phase velocity c. This is due to the fact that the thermal waves in the coupled theory travel with an infinite speed of propagation as opposed to a finite speed in the generalized case.Figure10 displays the variation of the Stoneley wave velocity and attenuation coefficient with respect to phase velocity c that has oscillatory behavior in the whole range of c for L-S theory, which changes from the positive to the negative gradually.Figure11 shows the variations of the Stoneley waves velocity and attenuation coefficient with respect to phase velocity c that has oscillatory behavior in the whole range of c for (G-L theory) that increases with increasing relaxation time τ0, θ0.Figure12 shows the variations of the Stoneley waves velocity and attenuation coefficient with respect to phase velocity c for G-N theory in the whole range of c for different values of the constant n* which changes from the positive to the negative gradually. In both figures, it is clear that the Stoneley wave velocity and attenuation coefficient have a nonzero value only in a bounded region of space. It is observed that the Stoneley wave velocity decreases with increasing the phase velocity, while the attenuation coefficient increases with increasing the constant  n*.Figure13 shows the variations of the Stoneley waves velocity and attenuation coefficient with respect to phase velocity c for G-N theory in the whole range of c for different values of constant n*. In both figureures, it is clear that the Stoneley wave velocity and attenuation coefficient have a nonzero value only in a bounded region of space. It is observed that the Stoneley wave velocity increases with increasing the phase velocity, while it decreases with increasing constant n* and attenuation coefficient increases with increasing the phase velocity, while it increases with increasing constant n*.Figures14, 15, 16, 17, 18, 19, 20, and 21 show the variations of the Rayleigh velocity and attenuation coefficient with respect to gravity field g and rotation Ω, respectively, which has an oscillatory behavior in the whole range of g and Ω for C-T theory, L-S theory, G-L theory, and G-N theory. These figureures indicate that the medium of attenuation coefficient along g undergoes compressive deformation due to the thermal shock, except (C-T theory) it is a tension deformation. In both figureures, it is clear that the Rayleigh wave velocity and attenuation coefficient have a nonzero value only in a bounded region of space. Figure 22 shows the variations of the Rayleigh wave secular with respect to phase velocity c, it is increases with increasing the phase velocity for C-T theory, L-S theory, G-L theory, and G-N theory, while it decreases with increasing constant n* for G-N theory. In both figureures, it is clear that the Rayleigh waves secular equation has a nonzero value only in a bounded region of space.Figures23, 24, and 25 show the variations of the Rayleigh wave velocity and attenuation coefficient with respect to phase velocity c, it has oscillatory behavior in the whole range of c for different values of relaxation times τ0, θ0, and n*, respectively, for L-S theory, G-L theory, and G-N theory, which changes from the positive to the negative gradually. In both figureures, it is clear that the Rayleigh wave velocity and attenuation coefficient have a nonzero value only in a bounded region of space.Figure26 shows the variations of the Rayleigh wave velocity and attenuation coefficient with respect to phase velocity c for G-N theory, the Rayleigh wave velocity increases with increasing phase velocity, while it decreases with increasing the constant n*, and the attenuation coefficient decreases with increasing phase velocity and the constant n*; it is clear that the Rayleigh wave velocity and attenuation coefficient have a nonzero value only in a bounded region of space.Figure27 shows the variation of the Love waves secular equation with respect to phase velocity c, which increases with increasing phase velocity for C-T theory, L-S theory, G-L theory, and G-N theory, while it decreases with increasing the constant n*; it is clear that Love waves secular equation has a nonzero value only in a bounded region of space.Figures28, 29, and 30 show the variations of the Love wave velocity and attenuation coefficient with respect to phase velocity c, it has oscillatory behavior in the whole range of c for different values of relaxation times τ0, θ0, and  n*, respectively, for L-S theory, G-L theory, and G-N theory, which changes from the positive to the negative gradually. In both figures, it is clear that Love wave velocity and attenuation coefficient have a nonzero value only in a bounded region of space.


## 9. Conclusion

Due to the complicated nature of the governing equations of the generalized thermoelasticity Fibre-reinforced theory, the work done in this field is unfortunately limited in number. The method used in this study provides a quite succesful in dealing with such problems. This method gives exact solutions in the generalized thermoelastic medium without any assumed restrictions on the actual physical quantities that appear in the governing equations of the problem considered. Important phenomena are observed in all these computations.(i)
It was found that for large values of time they give close results. The solutions obtained in the context of generalized thermoelasticity theory, however, exhibit the behavior of speed of surface wave propagation.(ii)
By comparing Figures1–30, it was found that the surface wave velocity has the same behavior in both media. But with the passage of rotation, relaxation times, and gravity, numerical values of surface wave velocity in the thermoelastic medium are large in comparison due to the influences of gravity and relaxation times in the elastic medium.(iii)
Special cases are considered as Rayleigh waves, Love wave, and Stoneley surface waves in anisotropic generalized thermoelastic medium, as well as in the isotropic case.(iv)
The results presented in this paper should prove useful for researchers in material science and designers of new materials.(v)
Study of the phenomenon of rotation and gravity is also used to improve the conditions of oil extractions.



---

*Source: 562369-2013-12-26.xml*