# Effect of Initial Stress on the Propagation Characteristics of Waves in Fiber-Reinforced Transversely Isotropic Thermoelastic Material under an Inviscid Liquid Layer

Authors: Rajneesh Kumar; Sanjeev Ahuja; S. K. Garg
Journal: Journal of Thermodynamics
 (2014)

---

## Abstract

The present investigation deals with the propagation of waves in fiber-reinforced transversely isotropic thermoelastic solid half space with initial stresses under a layer of inviscid liquid. The secular equation for surface equation in compact form is derived after developing the mathematical model. The phase velocity and attenuation coefficients of plane waves are studied numerically for a particular model. Effects of initial stress and thickness of the layer on the phase velocity, attenuation coefficient, and specific loss of energy are predicted graphically in the certain model. A particular case of Rayleigh wave has been discussed and the dispersion curves of the phase velocity and attenuation coefficients have also been presented graphically. Some other particular cases are also deduced from the present investigation.


---

## Body

## 1. Introduction

Fiber-reinforced materials are widely used in engineering structures due to their superiority over the structural materials in applications requiring high strength and stiffness in lightweight components. Consequently, characterization of their mechanical behavior is of particular importance for structural design using these materials. Fibers are assumed an inherent material property, rather than some form of inclusion in models as Spencer [1]. In the case of an elastic solid reinforced by a series of parallel fibers it is usual to assume transverse isotropy.The idea of continuous self-reinforcement at every point of an elastic solid was introduced by Belfield et al. [2]. The characteristic property of reinforced concrete member is that its components, namely, concrete and steel, act together as a single anisotropic unit as main long as they remain in the elastic condition; that is, the two components are bound together so that there can be no relative displacement between them. The dynamical interaction between the thermal and mechanical fields in solids has great practical applications in modern aeronautics, astronautics, nuclear reactors, and high energy particle accelerators.The analysis of stress and deformation of fiber-reinforced composite materials has been an important subject of solid mechanics for last three decades. Pipkin [3] did pioneer works on the subject. Sengupta and Nath [4] discussed the problem of surface waves in fiber-reinforced anisotropic elastic media.Lord and Shulman [5] introduced a theory of generalized thermoelasticity with one relaxation time for an isotropic body. The theory was extended for anisotropic body by Dhaliwal and Sherief [6]. In this theory, a modified law of heat conduction including both the heat flux and its time derivatives replaces the conventional Fourier’s Law. The heat equation associated with this theory is hyperbolic and hence eliminates the paradox of infinite speeds of propagation inherent in both coupled and uncoupled theories of thermoelasticity. Erdem [7] derived heat conduction equation for a composite rigid material containing an arbitrary distribution of fibers. Recently, Kumar and Gupta [8] discussed the Wave motion in an anisotropic fiber-reinforced thermoelastic solid.Surface waves carry a lot of information about the Earth’s crust and dispersion analysis of surface waves is concerned with the phase velocity and wave number. Further the velocity of seismic waves depends upon the elastic parameters of the medium through which they travel. Therefore, many investigators have studied the surface wave propagation by considering different kinds of models. Chadwick and Seet [9] and, Singh and Sharma [10] have discussed the propagation of plane harmonic waves in transversely isotropic thermoelastic materials. Singh [11] studied a problem on wave propagation in an anisotropic generalized thermoelastic solid and obtained a cubic equation, which gives the dimensional velocities of various plane waves.Chadwick [12] has discussed the propagation of surface waves in homogeneous thermoelastic media. Chadwick and Windle [13] studied the effect of heat conduction on the propagation of Rayleigh waves in the semi-infinite media (i) when the surface is maintained at constant temperature and (ii) when the surface is thermally insulated. Sharma and Singh [14] have studied thermoelastic surface waves in a transversely isotropic half-space.The study of wave propagation in a generalized thermoelastic media with additional parameters like prestress, porosity, viscosity, microstructure, temperature, and other parameters provide vital information about existence of new or modified waves. Such information may be useful for experimental seismologists in correcting earthquake estimation.Prestressed materials have various applications, for example, in oil and geophysical industry, NDT in prestressed materials, use of rubber composites in automotive, and aerospace and defense industries (often in prestressed states) and in the study of biological tissues (lung, tendon, etc.) which are all nonlinear prestressed viscoelastic composites.Montanaro [15] investigated the isotropic linear thermoelasticity with hydrostatic initial stress. Wang and Slattery [16] formulated the thermoelastic equations without energy dissipation for a body which has previously received a large deformation and is at nonuniform temperature. Ieşan [17] presented a theory of Cosserat thermoelastic solids with initial stresses. Ames and Straughan [18] derived the continuous dependence results for initially prestressed thermoelastic bodies.Some theorems in the generalized theory of thermoelasticity for prestressed bodies were studied by Wang et al. [19]. Marin and Marinescu [20] studied the asymptotic partition of total energy for the solutions of the mixed initial boundary value problem within the context of the thermoelasticity of initially stressed bodies. Kalinchuk et al. [21] studied the problem of steady-state harmonic oscillations for a nonhomogeneous thermoelastic prestressed medium.Othman and Song [22] discussed the reflection of plane waves from a thermoelastic elastic solid half-space under hydrostatic initial stress without energy dissipation. Singh [23] studied wave propagation in an initially stressed transversely isotropic thermoelastic half-space. However, no attempt has been made to discuss the effect of hydrostatic initial stress on propagation of waves in the layer of fiber-reinforced transversely isotropic thermoelastic material. Othman et al. [24] investigated the influences of fractional order, hydrostatic initial stress, and gravity field on the plane waves in a fiber-reinforced isotropic thermoelastic medium using normal mode analysis. Abd-Alla et al. [25] studied the effects of rotation and gravity field on surface waves in fiber-reinforced thermoelastic media under four theories. Othman and Atwa [26] studied the effect of rotation, gravity, and reinforcement on the total deformation and mutual interaction of the body in a fiber-reinforced thermoelastic solid using normal mode analysis in the context of Green-Naghdi theory of thermoelasitcity. Abbas [27] studied the transient phenomena in the magnetothermoelastic model in the context of the Lord and Shulman theory in a perfectly conducting medium. Abbas [28] constructed the equations for generalized thermoelasticity of an unbounded fiber-reinforced anisotropic medium with a circular hole in the context of Green and Naghdi (GN) theory. Recently, Abbas and Zenkour [29] investigated two-temperature generalized thermoelasticity with two relaxation times for an infinite fiber-reinforced anisotropic plate containing a circular cavity using a finite element method.In the present paper, the governing equations of thermoelastcity of fiber-reinforced transversely isotropic solid half space with initial stresses under a layer of inviscid liquid are formulated and solved analytically in two dimensions at uniform temperature. The phase velocity and attenuation coefficients of plane waves are studied numerically for a particular model. Specific loss of energy is obtained and discussed numerically for a particular model to study the effects of initial stress. The study of present investigation has applications in civil engineering and geophysics.


## 2. Basic Equations

The basic equations in the dynamic theory of the homogeneous thermally conducting fiber-reinforced medium with an initial hydrostatic stress without body forces and heat sources are given by Lord and Shulman [5] and Abbas and Othman [30] as
(1)σij,j-Pωij,j=ρu¨i,i,j=1,2,3
and heat conduction equation is given by
(2)KijT,ij=(∂∂t+τ0∂2∂t2)(Toβiju˙i,j+ρCeT˙),i,j=1,2,3.The constitutive relations for thermally conducting transversely isotropic, fiber-reinforced linearly elastic media are(3)σij=λekkδij+2μTeij+α(akamekmδij+aiajekk)+2(μL-μT)(aiakekj+ajakeki)+β(akamekmaiaj)-βijTδij,i,j,k,m=1,2,3,
where
(4)eij=12(ui,j+uj,i),ωij=12(uj,i-ui,j),i,j=1,2,3,
and λ,μT are elastic parameters, α,β,(μL-μT) are reinforced anisotropic elastic parameters, ρ is the mass density, σij are components of stress tensor, ui are the displacement components, eij are components of infinitesimal strain, T is the temperature change of a material particle, To is the reference uniform temperature of the body, kij are coefficients of thermal conductivity, βij are thermal elastic coupling tensor, Ce is the specific heat at constant strain, δij is the Kronecker delta, and P is the initial Pressure. The comma notation is used for spatial derivatives and superimposed dot represents time differentiation. aj are components ofa, all referred to Cartesian coordinate. The vectora may be a function of position. We choosea so that its components are (1, 0, 0).The equation which governs the motion in terms of velocity potentialψ of the homogeneous inviscid liquid is given by  Ewing et al. [31]:(5)∂2ψ∂x12+∂2ψ∂x22=1(αl)2∂2ψ∂t2,
where αl=λl/ρl is the velocity of sound liquid.The displacement componentsu1l,u2l and pressure pl in the medium M2 are given by
(6)u1l=∂ψ∂x1,u2l=∂ψ∂x2,pl=-ρl∂ψ∂t.


## 3. Formulation of the Problem

We consider a homogeneous, thermally conducting transversely isotropic fiber-reinforced medium (M1) with an initial hydrostatic stress lying under a uniform homogeneous inviscid liquid layer medium (M2) of thickness H. We take the rectangular Cartesian coordinate system Ox1x2x3 at any point on the plane horizontal surface and x1-axis in the direction of wave propagation and x2-axis pointing vertically downward into the half-space so that all particles on a line parallel to x3-axis are equally displaced; therefore, all the field quantities will be independent of x3-coordinates. Hence the liquid layer medium M2 occupies the region -H<x2<0 and x2>0 is occupied by the half-space medium M1. The plane x2=-H represents the free surface of the liquid and x2=0 is taken as interface between M1 and M2 medium (Figure 1).Figure 1
Geometry of the problem.The displacement components for mediumM1 are taken as
(7)u→=(u1,u2,0).With the help of (3), (4), and (7) and (1)-(2), take the form
(8)C11∂2u1∂x12+(C12+C0-P2)∂2u2∂x1∂x2+(C0+P2)∂2u1∂x22-β11∂T∂x1=ρ∂2u1∂t2,(9)C22∂2u2∂x22+(C12+C0-P2)∂2u1∂x1∂x2+(C0+P2)∂2u2∂x12-β22∂T∂x2=ρ∂2u2∂t2,(10)∂2T∂x12+k-∂2T∂x22-ρCek11(∂T∂t+τo∂2T∂t2)=ε[(∂2u1∂x1∂t+τo∂3u1∂x1∂t2)+β-(∂2u2∂x2∂t+τo∂3u2∂x2∂t2)],
where
(11)β11=(C11+C12)α11+C12α22,β22=(C12+C23-C55)α11+C22α22,C12=λ+α,C11=λ+2α+4μL-2μT+β,C22=C33=λ+2μT,C44=C66=2μL,C55=2μT,C23=C33-C55,Co=C442,ε=Toβ11k11,k-=k22k11,β-=β22β11,
and λ,α,β,μL,μT are material constants, α11,α22 are components of linear thermal expansion, and τo is thermal relaxation time.To facilitate the solution, the following dimensionless quantities are introduced:(12)(x1′,x2′)=w*v1(x1,x2),(u1′,u2′)=ρv1w*β11To(u1,u2),tij′=tijβ11To,T′=TTo,t′=w*t,h′=v1w*h,τo′=w*τo,ψ′=w*ψv12,(u1l′,u2l′)=w*v1(u1l,u2l),pl′=plλl,t22l′=t22lλl,
where w*=CeC11/k11,v12=C11/ρ and w* is the characteristic frequency of the medium.


## 4. Solution of the Problem

We assume the solutions of the form(13)(u1,u2,T,ψ)=(1,W,R,S)Uexp⁡{iξ(x1+mx2-ct)},
where c=ω/ξ is the nondimensional phase velocity, ω is the frequency, and ξ is the wave number. m is still unknown parameter. 1,W,R,S are, respectively, the amplitude ratios of u1,u2,T, and u2l with respect to u1. Substituting the values of u1,u2  and  T from (13) in (8), (9), and (10), we obtain
(14)[mC11(C12+Co-P2)]W+(-1iξ)R+[1-c2+1C11(Co+P2)m2]=0,[C22C11m2+1C11(Co+P2)-c2]W+(-mβ-iξ)R+[1C11(C12+Co-P2)m]=0,[ε1β-(-mc+mc2τoiξ)]W+[1+k-m2+(ciξ-τoc2)]R+ε1[c-τoc2iξ]=0,
where
(15)ε1=Toβ112ρw*k11.The system of (14) has a nontrivial solution if the determinant of the coefficients [1,W,R]T vanishes, which yields to the following polynomial characteristic equation:
(16)m6+d21m4+d31m2+d41=0,
where
(17)d1=[k-(L2g1)],d2=[L2g1+(β-)2ε1FL2-g1FL2+k-(L1g1+L2g2-A12)+L2g1+(β-)2ε1FL2],d3=[L2g2+(β-)2ε1FL1+A12(F-1)-Fg1(L1+ε1)-Fg2L2+β-FA1(1-ε1)+k-(L1g2)+L1g1+L2g2+(β-)2ε1FL1+A12(F-1)],d4=[-g2F(L1+ε1)+L1g2],d21=d2d1,d31=d3d1,d41=d4d1,F=-c+τoc2iξiξ,A=mA1,A1=1C11(C12+Co-P2),G=g1m2+g2,g1=C22C11,g2=1C11(Co+P2)-c2,L=L1+L2m2,L1=1-c2,L2=1C11(Co+P2).The characteristic equation (16) is cubic in m2 and hence possesses three roots (mp2,p=1,2,3). Therefore, there exist three types of quasi-waves in transversely isotropic elastic half-space, namely, quasi-longitudinal waves (QP), quasi-transverse waves (QSV), and quasi-thermal waves (QT).The formal expression for displacements and volume fraction field, satisfying the radiation condition thatRe(mp)≥0, can be written as
(18)(u1,u2,T)=∑p=13Gp(1,n1p,n2p)exp⁡{iξ(x1+impx2-ct)},
where
(19)n1p=Δ1Δ,n2p=Δ2Δ,Δ1=(-β-mε1c(1-τociξ)iξ-AY),Δ2=(Gε1c(1-τociξ)-AZ),Δ=(GY+β-mZiξ),Y=1+k-m2+(ciξ-τoc2),Z=ε1β-mc(-1+cτoiξ)
and G1,G2,G3 are arbitrary constants.Substituting the values ofψ from (13) in (5) and after simplification, we obtain
(20)ψ=[G4sin(ξm4x2)+G5cos⁡(ξm4x2)]exp⁡(iξ(x1-ct)),
where
(21)m4=±((v1cαl)2-1)1/2
and G4,G5 are arbitrary constants.The normal displacement and normal stress for mediumM2, with the help of (5), (6), (12), and (20), are given by
(22)u2l=m4ξ[G4cos⁡(ξm4x2)-G5sin(ξm4x2)]×exp⁡(iξ(x1-ct)),t22l=-ξ2(1+m42)[G4sin(ξm4x2)+G5cos⁡(ξm4x2)]×exp⁡(iξ(x1-ct)).Also from (3) we have
(23)t22=C12u1,1+C22u2,2-β22T,t21=Co(u1,2+u2,1).


## 5. Boundary Conditions

The appropriate boundary conditions are as follows:(i)
vanishing of normal stress component at the free surface of the liquid layer,x2=-H; that is,
(24)(t22l)M2=-pl=0,(ii)
continuity of normal stress components at the interfacex2=0; that is,
(25)(t22)M1=-pl,(iii)
normal displacement components at the interfacex2=0; that is,
(26)(u2l)M2=(u˙2)M1,(iv)
vanishing shear stress for mediumM1 at the interface x2=0; that is,
(27)(t21)M1=0,(v)
thermal boundary conditions(28)∂T∂x2+hT=0,whereh⟶0 corresponds Insulated boundary condition, h⟶∞ corresponds Isolated boundary condition.


## 6. Derivations of the Secular Equations

Substituting the value ofu2,T,u2l,t22l from (18) and (22) in (24)–(28) and with the aid of (7) and (12), after simplification, we obtain
(29)c2=(αlv1)2[(tan-1(T11))2+ξ2H2ξ2H2],
where
(30)T11=C11m4(E11D1-E12D2+E13D3)icρlv12(n11D1-n12D2+n13D3),E1p=iξ(C12C11)+(C22C11)n1p(i2ξmp)-(β22β11)n2p,E2p=ξ(n1pi-mp),E3p=n2p(h-ξmp),(p=1,2,3),D1=E22E33-E32E23,D2=E21E33-E31E23,D3=E21E32-E31E22.Equation (29) is surface waves secular equation of fiber-reinforced transversely isotropic thermoelastic half-space under a uniform homogeneous inviscid liquid layer. The secular equation has complete information about the phase velocity, wave number, and attenuation coefficient of surface waves in such medium.


### 6.1. Particular Cases

(i)
If we takeh→0 in (29), we obtain the frequency equation for the thermally insulated boundary.
(ii)
If we takeh→∞ in (29), we obtain the frequency equation for the thermally isolated boundary.
(iii)
TakingP→0 in secular equation (29), yield the frequency equation for fiber-reinforced transversely isotropic thermoelastic half space without initial stress under a uniform homogeneous inviscid liquid layer.
(iv)
ForμT=μL=μ,  K11=K22=K, β11=β22=β=(3λ+2μ)αt, (29) reduces to the frequency equation in a homogeneous Isotropic half space under inviscid liquid.
(v)
τ
o
→
0 in (10), we obtained the frequency equation in fiber-reinforced coupled thermoelastic half space under inviscid liquid layer.


## 6.1. Particular Cases

(i)
If we takeh→0 in (29), we obtain the frequency equation for the thermally insulated boundary.
(ii)
If we takeh→∞ in (29), we obtain the frequency equation for the thermally isolated boundary.
(iii)
TakingP→0 in secular equation (29), yield the frequency equation for fiber-reinforced transversely isotropic thermoelastic half space without initial stress under a uniform homogeneous inviscid liquid layer.
(iv)
ForμT=μL=μ,  K11=K22=K, β11=β22=β=(3λ+2μ)αt, (29) reduces to the frequency equation in a homogeneous Isotropic half space under inviscid liquid.
(v)
τ
o
→
0 in (10), we obtained the frequency equation in fiber-reinforced coupled thermoelastic half space under inviscid liquid layer.


## 7. Rayleigh Wave

(1) TakingH→0, the frequency equation (29) reduces to
(31)E11(E22E33-E23E32)-E12(E21E33-E31E23)+E13(E21E32-E31E22)=0.Equation (31) is the frequency equation for Rayleigh surface wave propagation in a fiber-reinforced transversely isotropic thermoelastic half-space with initial stress. In general wave number and hence phase velocity of waves is complex quantity therefore the waves are attenuated in space. If we write
(32)c-1=v-1+iω-1F
so that, ξ=K+iF, where K=ω/v, v and F are real. Also the roots of characteristic equation (16) are in general complex and hence we assume that mp=pp+iqp, so that the exponent in the plane wave solutions (18) for the half-space becomes
(33)iK(x1-mpIx2-vt)-K(FKx1+mpRx2),
where
(34)mpR=pp-qpFK,mpI=qp+ppFK,(p=1,2,3).This shows thatv is the propagation velocity and F is the attenuation coefficient of wave.The equation (18) can be rewritten as
(35)u1=∑p=13Gpexp⁡(-Fx1-λpRx2)×exp⁡[i(K(x1-vt)-λpIx2)],u2=∑p=13n1pGpexp⁡(-Fx1-λpRx2)×exp⁡[i(K(x1-vt)-λpIx2)],T=∑p=13n2pGpexp⁡(-Fx1-λpRx2)×exp⁡[i(K(x1-vt)-λpIx2)]
with λp=K(mpR+impI)=λpR+iλpI(say), (p=1,2,3). Moreover it is clear that
(36)|λpR|2-|λpI|2=K2[(mpR)2-(mpI)2].


### 7.1. Normal Stress

Normal stress at the interfacex2=0 is given by
(37)t22=(C12iξ∑p=13Gp-C22ξ∑p=13mpGpn1p-β22∑p=13Gpn2p)×exp⁡(iξ(x1-vt)),(p=1,2,3).


### 7.2. Specific Loss

The Specific loss is the ratio of energy(ΔW) dissipated in taking a specimen through a stress cycle, to the elastic energy (W) stored in the specimen when the strain is maximum. The Specific loss is the most direct method of defining internal friction for a material. For a sinusoidal plane wave of small amplitude, Kolsky [32] shows that the specific loss ΔW/W equals 4π times the absolute value of imaginary part of ξ to the real part of ξ; that is,
(38)ΔWW=4π|Im⁡(ξ)Re(ξ)|=4π|vGω|=4π|GF|.


## 7.1. Normal Stress

Normal stress at the interfacex2=0 is given by
(37)t22=(C12iξ∑p=13Gp-C22ξ∑p=13mpGpn1p-β22∑p=13Gpn2p)×exp⁡(iξ(x1-vt)),(p=1,2,3).


## 7.2. Specific Loss

The Specific loss is the ratio of energy(ΔW) dissipated in taking a specimen through a stress cycle, to the elastic energy (W) stored in the specimen when the strain is maximum. The Specific loss is the most direct method of defining internal friction for a material. For a sinusoidal plane wave of small amplitude, Kolsky [32] shows that the specific loss ΔW/W equals 4π times the absolute value of imaginary part of ξ to the real part of ξ; that is,
(38)ΔWW=4π|Im⁡(ξ)Re(ξ)|=4π|vGω|=4π|GF|.


## 8. Numerical Results and Discussion

For numerical computations, we take the following values of the relevant parameters for generalized fiber-reinforced transversely isotropic thermoelastic solid:(39)λ=5.66×1010N/m2,ρ=2660Kg/m3,μT=2.46×1010N/m2,μL=5.66×1010N/m2,α=-1.28×1010N/m2,β=220.90×1010N/m2,K11=0.0921×103Jm-1deg-1s-1,P=100K22=0.0963×103Jm-1deg-1s-1,α1=0.017×104deg-1,α2=0.015×104deg-1,To=293K,τo=0.05s,h=1,H=1,Ce=0.787×103JKg-1deg-1,ω=2s-1.The physical constants for water are given by Ewing et al. [31]:(40)ρl=1.0×103Kg/m3,λl=2.14×103N/m2.For comparison with the generalized fiber-reinforced isotropic thermoelastic solid following Ezzat [33], we take the following values of relevant parameters for the case of copper material as
(41)λ=7.76×1010N/m2,μT=μL=μ=3.86×1010N/m2,K11=K22=K=386NK-1s-1,β11=β22=β=(3λ+2μ)αt,αt=1.78×10-5K-1α1=1.005deg-1,α2=1.001deg-1,Ce=383.1m2K-1s-1,P=0,H=1,ρ=8954Kg/m3,τo=0.1s,h=1.Using the above values of parameters, the dispersion curves of nondimensional phase velocity, attenuation coefficient and specific heat are shown graphically with respect to nondimensional wave number. Here, we consider abbreviations FTTIWIS, FTTIIS, FTIWIS, and FTIIS for fiber-reinforced transversely isotropic thermoelastic material without initial stress, fiber-reinforced transversely isotropic thermoelastic material with initial stress, fiber-reinforced isotropic thermoelastic material without initial stress, and fiber-reinforced isotropic thermoelastic material with initial stress respectively andH refers to the thickness of the layer.Figure2 depicts the variations of phase velocity with respect to the wave number. The behavior of phase velocity for FTTIWIS, FTTIIS, FTIWIS, and FTIIS is similar but the corresponding values are different in magnitude. The values of phase velocity decrease for lower wave number and become stable for higher wave number in all the cases. Figure 3 shows the variation of attenuation coefficient with respect to the wave number. For FTTIWIS and FTTIIS the values of attenuation coefficient increase monotonically up to wave number (ξ) = 2 and after that it becomes constant for whole range of wave number. On the other hand, the values of attenuation coefficients are consistent within the whole range of wave number for FTIWIS and FTIIS.Figure 2
Variation of phase velocity w.r.t. wave number (with and without initial stress).Figure 3
Variation of attenuation coefficient w.r.t. wave number (with and without initial stress).Figure4 indicates the variation of phase velocity with respect to the wave number for different values of thickness of layer. Here, we consider the particular case of fiber-reinforced transversely isotropic thermoelastic material and analyses of the behavior of curves with respect to initial stress. It is observed that for small value of wave number the graph of phase velocity shows opposite behavior w.r.t. thickness, that is, on increasing the values of thickness of layer, the values of corresponding phase velocity decreases, and it became steady for large values of wave number.Figure 4
Variation of phase velocity w.r.t. wave number (with and without initial stress for different value ofH).Figure5 shows the variation of attenuation coefficient with respect to the wave number for different value of thickness of the layer. Here also we considered the particular case of fiber-reinforced transversely isotropic thermoelastic material and analyses of the behavior of curves with respect to initial stress. It is perceived that on decreasing the value of H the curve shows consistency for small value of wave number and on increasing the value of H the curve shows uniformity for large value of ξ.Figure 5
Variation of attenuation coeff. w.r.t. wave number (with and without initial stress for different value ofH).Figure6 depicts the variations of specific loss verses wave number. For FTIWIS and FTIIS the curve shows sharp decrease in the value of specific loss for the range 0≤ξ≤1 and become parallel for ξ>1. It is observed that graphs corresponding to FTTIWIS and FTTIIS follow the similar pattern for whole range of ξ but the values are different in magnitude.Figure 6
Variation of specific loss w.r.t. wave number (with and without initial stress).


### 8.1. Rayleigh Wave

Figure7 shows the variations of phase velocity of Rayleigh wave verses wave number for FTTIWIS, FTTIIS, FTIWIS, and FTIIS. The value of phase velocity decreases for small values of ξ and become stable for higher value of ξ. It is analyzed that in all the cases behavior of curves are almost the same within the range 2≤ξ≤5. Figure 8 shows the variations of attenuation coefficient of Rayleigh wave verses wave number for FTTIWIS, FTTIIS, FTIWIS, and FTIIS. For FTIWIS and FTIIS the graph shows a sudden fall for ξ=0.5 and then increases monotonically up to ξ=2 and becomes constant for the remaining range of ξ. A large difference between the behavior of curves for the transversely isotropic and isotropic case is also perceived.Figure 7
Variation of phase velocity w.r.t. wave number for Rayleigh wave (with and without initial stress).Figure 8
Variation of attenuation coeff. w.r.t. wave number for Rayleigh wave (with and without initial stress).


## 8.1. Rayleigh Wave

Figure7 shows the variations of phase velocity of Rayleigh wave verses wave number for FTTIWIS, FTTIIS, FTIWIS, and FTIIS. The value of phase velocity decreases for small values of ξ and become stable for higher value of ξ. It is analyzed that in all the cases behavior of curves are almost the same within the range 2≤ξ≤5. Figure 8 shows the variations of attenuation coefficient of Rayleigh wave verses wave number for FTTIWIS, FTTIIS, FTIWIS, and FTIIS. For FTIWIS and FTIIS the graph shows a sudden fall for ξ=0.5 and then increases monotonically up to ξ=2 and becomes constant for the remaining range of ξ. A large difference between the behavior of curves for the transversely isotropic and isotropic case is also perceived.Figure 7
Variation of phase velocity w.r.t. wave number for Rayleigh wave (with and without initial stress).Figure 8
Variation of attenuation coeff. w.r.t. wave number for Rayleigh wave (with and without initial stress).


## 9. Conclusions

The propagation of surface waves in an initially stressed fiber-reinforced transversely isotropic thermoelastic material under an inviscid liquid layer has been studied. After developing the formal solutions, the secular equation for surface wave propagation is derived. Some special cases of frequency equation are also discussed. The attenuation coefficients and specific loss (liquid layer and half-space) have also been studied and computed for FTTIWIS, FTTIIS, FTIWIS, and FTIIS. The trend of phase velocity is the same in all these four cases for surface wave but varies for attenuation coefficient. From numerical analysis, it is shown that the phase speeds and attenuation coefficients for Rayleigh wave are also affected significantly due to the presence of initial stress parameters.



---

*Source: 134276-2014-08-26.xml*