The Rayleigh surface wave is studied at a stress-free thermally insulated surface of an isotropic, linear, and homogeneous thermoelastic solid half-space with microtemperatures. The governing equations of the thermoelastic medium with microtemperatures are solved for surface wave solutions. The particular solutions in the half-space are applied to the required boundary conditions at stress-free thermally insulated surface to obtain the frequency equation of the Rayleigh wave. Some special cases are also derived. The non-dimensional speed of Rayleigh wave is computed numerically and presented graphically to reveal the dependence on the frequency and microtemperature constants.
1. Introduction
The theory of materials with microstructures has been a subject of intensive study in the literature since E. Cosserat and F. Cosserat [1]. The microtemperature and/or microdeformation of the nanoparticles could be considered very important in future technologies. The thermoelasticity with microtemperatures considers the microstructure of the body, in which each microelement possesses a microtemperature. The theory of thermodynamics for elastic material with innerstructures was developed by Grot [2] according to which the molecules possess microtemperatures along with macrodeformation of the body. The experimental data for the silicone rubber containing spherical aluminum particles and for human blood presented by Říha [3] conform closely to the predicted theoretical model of thermoelasticity with microtemperatures.
Iesan and Quintanilla [4] developed the linear theory for elastic materials with inner structure whose particles, in addition to the classical displacement and temperature fields, possess microtemperatures. Ieşan [5] proposed a theory of thermoelastic bodies with microstructures and microtemperatures, where the microelements of the material possess microtemperatures and can undergo microrotation, microstretch and translation. Ieşan and Quintanilla [6] presented the theory of thermoelastic bodies with inner structure and microtemperatures, which permits the transmission of heat as thermal waves at finite speed. Several papers based on the theory of thermoelasticity with microtemperatures have been published such as Iesan [7], Svanadze [8, 9], Casas and Quintanilla [10], Scalia and Svanadze [11, 12], Aoudai [13], Iesan [14], Scalia et al. [15], Ieşan and Scalia [16], Yang and Huang [17], Quintanilla [18], Svanadze and Tracinà [19], Chirita et al. [20], and Steeb et al. [21].
Due to increasing interest in nanomaterials, the significance of microtemperature and/or microdeformation of the nanoparticles cannot be ignored. The studies related to wave propagation in the theory of thermoelastic materials with microtemperature may be important in future technologies. The theory of thermoelasticity with microtemperatures (Iesan and Quintanilla [4]) is applied to study the Rayleigh wave at the thermally insulated stress-free surface of an isotropic, homogeneous thermoelastic solid half-space with microtemperature. The frequency equation of the Rayleigh wave is obtained. The dependence of numerical values of the speed of the Rayleigh wave on material parameters, frequency, and microtemperature constants is shown graphically for a particular material of the model.
2. Basic Equations
Following Iesan and Quintanilla [4], the constitutive relations for homogeneous and isotropic thermoelastic medium with microtemperatures are
(1)σij=λδijerr+2μeij-βTδij,qi=KT,i+K1wi,qij=-K4wr,rδij-K5wi,j-K6wj,i,Qi=(K1-K2)wi+(K-K3)T,i,ρη*=βerr+aT,ρϵi=-bwi,
where
(2)eij=12(ui,j+uj,i),
and λ, μ, β, a, b, K and Kj(j=1,2,…,6) are constitutive coefficients. σij are the components of the stress tensor. eij are the components of the strain tensor. ρ is the reference mass density of the medium. η* is entropy per unit mass. ϵi are the components of the first moment of energy vector. qij are the components of the first heat flux moment vector. Qi are the components of the mean heat flux vector. qi are components of the heat flux vector. ui are the components of the displacement vector u→. wi are the components of the microtemperature vector w→. T=Θ-T0, where Θ is the temperature at time t. T0 is the temperature of the medium in its natural state and assumed to be such that |T/T0|<<1. A comma in the subscript denotes the spatial derivative and δij is the Kronecker delta.
Following Iesan and Quintanilla [4], the constitutive equations (1) combined with the reduced Clausius-Duhem inequality in context of the linear theory of thermoelasticity with microtemperature imply the following inequalities:
(3)3K4+K5+K6≥0,K6+K5≥0,K6-K5≥0,K≥0,(K1+T0K3)2-4T0KK2≤0.
Following Iesan and Quintanilla [4], the fundamental system of field equations of the linear theory of thermoelasticity with microtemperatures
the equations of motion
(4)σji,j+ρfi=ρu¨i,
the balance energy
(5)ρT0η*˙=qi,i+ρS,
the first moment of energy
(6)ρϵ˙i=qji,j+qi-Qi+ρMi,
where fi are the components of the body force vector, Mi are the components of the first heat source moment vector, and S is the heat supply. Superposed dot represents the temporal derivative and other symbols are described previously.
Using (1) and (2) in (4) to (6), the following system of linear partial differential equations is obtained:
(7)μui,jj+(λ+μ)uj,ij-βT,i+ρfi=ρu¨i,KT,ii-βT0u˙i,i+K1wi,i+ρS=aT0T˙,K6wi,jj+(K4+K5)wj,ji-K3T,i-K2wi-ρMi=bw˙i.
The field equations (7) in term- of displacement, macro- and microtemperatures for a linear homogeneous elastic solid in the absence of body force, heat source, and first-heat source moment vector are written in the following form:
(8)μ∇2u→+(λ+2μ)∇(∇·u→)-β∇T=ρu→¨,K∇2T-βT0∇·u→˙+K1(∇·w→)-aT0T˙=0,K6∇2w→+(K4+K5)∇(∇·w→)-K3∇T-K2w→-bw→˙=0.
3. Analytical 2D Solution
We consider a homogeneous and isotropic thermoelastic medium of an infinite extent with Cartesian coordinates system (x,y,z), which is previously at uniform temperature. The origin is taken on the plane surface and the z-axis is taken normally into the medium (z≥0). The surface z=0 is assumed stress free and thermally insulated. The present study is restricted to the plane strain parallel to the x-z plane, with the displacement vector u→=(u1,0,u3). Introducing the scalar potentials ϕ and ξ, and vector potential ψ→ through Helmholtz representation of a vector field as
(9)u→=∇ϕ+∇×ψ→,∇·ψ→=0,w→=∇ξ.
Inserting (9) in (8), we obtain
(10)μ∇2ψ→=ρψ→¨,(11)(λ+2μ)∇2ϕ-βT=ρϕ¨,(12)K∇2T-βT0∇2ϕ˙-aT0T˙+K1∇2ξ=0,(13)(K4+K5+K6)∇2ξ-K2ξ-bξ˙-K3T=0.
For thermoelastic surface waves in the half-space propagating in x-direction, the potential functions ϕ, T, and ξ are taken in the following form:
(14){ϕ,T,ξ,ψ}={ϕ^(z),T^(z),ξ^(z),ψ^(z)}exp{ι(ηx-χt)},
where χ2=η2c2, η is the wave number, c is the phase velocity, and ψ=(ψ→)y. Substituting (14) in (11) to (13) and eliminating ϕ^, T^, and ξ^, we obtain the following auxiliary equation:
(15)D6-η2A0D4+η4A1D2-η6A2=0,
where D=d/dz and
(16)A0=3-(a1+a2+a3+a4+a5)c2,A1=3-2(a1+a2+a3+a4+a5)c2A1=+(a1a2+a2a3+a2a4+a3a4+a4a5)c4,A2=1-(a1+a2+a3+a4+a5)c2A2=+(a1a2+a2a3+a2a4+a3a4+a4a5)c4-(a2a3a4)c6,a1=ββ0Kμ0χ2,a2=K8K7χ2,a3=a0Kχ2,a4=ρμ0,a5=K1K3KK7χ2,a0=iχaT0,β0=iχβT0,K7=K4+K5+K6,K8=ibχ-K2.
Taking into account (15) and keeping in mind that ϕ^,T^,ξ^→0 as z→∞ for surface waves, the solutions ϕ, T, and ξ are written as
(17)ϕ=[Aexp(-ηλ1z)+Bexp(-ηλ2z)ϕ=+Cexp(-ηλ3z)]exp{ι(ηx-χt)},T=[p1Aexp(-ηλ1z)+p2Bexp(-ηλ2z)T=+p3Cexp(-ηλ3z)]exp{ι(ηx-χt)},ξ=[q1Aexp(-ηλ1z)+q2Bexp(-ηλ2z)ξ=+q3Cexp(-ηλ3z)]exp{ι(ηx-χt)},
where
(18)λ12+λ22+λ32=A0,λ12λ22+λ22λ32+λ12λ32=A1,λ12λ22λ32=A2,pi=η2{ρc2-(λ+2μ)(1-λi2)}β,(i=1,2,3),qi=K3piK8-η2(1-λi2)K7,(i=1,2,3).
Substituting (14) in to (10) and keeping in mind that ψ^→0 as z→∞ for surface waves, we obtain the following solution:
(19)ψ=Eexp[-ηλ4z+ι(ηx-χt)],
where
(20)λ42=1-c2c22,c22=μρ.
4. Derivation of Frequency Equation
The mechanical and thermal conditions at the thermally insulated surface z=0 are as follows:
vanishing of the normal stress component
(21)σzz=0,
vanishing of the tangential stress component
(22)σzx=0,
vanishing of the normal heat flux component
(23)qz=0,
vanishing of normal first heat flux moment vector component
(24)qzz=0,
where
(25)σzz=λ(∂2ϕ∂x2+∂2ϕ∂z2)+2μ(∂2ψ∂x∂z)+2μ∂2ϕ∂z2-βT,σzx=μ[2∂2ϕ∂x∂z-∂2ψ∂z2+∂2ψ∂x2],qz=K∂T∂z+K1∂ξ∂z,qzz=K4∂2ξ∂x2+K7∂2ξ∂z2.
Using the solutions (17) and (19) for ϕ, T, ξ, and ψ in (21) to (24) and eliminating A, B, C, and E, the following equation is obtained:
(26)X1Y1+X2Y2+X3Y3=0,
where
(27)X1=(1-λ42)s1+4μλ1λ4,X2=(1-λ42)s2+4μλ2λ4,X3=(1-λ42)s3+4μλ3λ4,Y1=λ2q3(Kp2+K1q2)(K4-K7λ32)Y1=-λ3q2(Kp3+K1q3)(K4-K7λ22),Y2=λ3q1(Kp3+K1q3)(K4-K7λ12)Y2=-λ1q3(Kp1+K1q1)(K4-K7λ32),Y3=λ1q2(Kp1+K1q1)(K4-K7λ22)Y3=-λ2q1(Kp2+K1q2)(K4-K7λ12),si=λ(1-λi2)-2μλi2+βpiη2,(i=1,2,3).
Equation (26) is the frequency equation of Rayleigh wave in thermoelastic medium with microtemperature.
5. Special Cases5.1. Isotropic Thermoelastic Case
In the case where microtemperature is absent, that is, when K1=K2=K3=K4=K5=K6=b=0, the frequency equation (26) is reduced to
(28)(1-λ42)(λ2*p2*s1*-λ1*p1*s2*)+4μλ1*λ2*λ4(p2*-p1*)=0,
where
(29)λ1*2+λ2*2=2-c2(a1+a3+a4),λ1*2λ2*2=1-c2(a1+a3+a4)+a3a4c4,pi*=η2{ρc2-(λ+2μ)(1-λi*2)}β,(i=1,2),si*=λ(1-λi*2)-2μλi*2+βpi*η2,(i=1,2).
Equation (28) is the frequency equation of Rayleigh wave in an isotropic thermoelastic solid.
5.2. Isotropic Elastic Case
In the case where thermal parameters are neglected, the frequency equation (28) is reduced to
(30)(2-c2c22)2=41-c2c121-c2c22,
which is the frequency equation of Rayleigh wave for an isotropic elastic case.
6. Numerical Example
The speed of propagation of Rayleigh wave is computed for the following physical constants of the model: λ=2.17×1010 N·m−2, μ=3.278×1010 N·m−2, ρ=1.74×103 kg·m−3, aT0=1.8×106 J·m−3 deg−1, β=2.68×106 N·m−2 deg−1, K=1.7×102 W·m−1 deg−1, K1=2×1010 W·m−1, K2=0.1×1010 W·m−1, K3=0.4×1010 W·m−1, K4=0.3×1010 W·m−1, K5=0.5×1010 W·m−1, K6=0.7×1010 W·m−1, b=1.3849×1010 N, T0=298K, and χ=3.58×1011 s−1.
The non-dimensional speed ρc2/μ of Rayleigh wave is computed and plotted against frequency for the range 2 Hz ≤χ≤30 Hz. It increases very sharply with the increase of frequency as shown in Figure 1. The non-dimensional speed ρc2/μ of Rayleigh wave is also computed for certain ranges of microtemperature constants. The non-dimensional speed ρc2/μ of Rayleigh wave remains almost constant at different values of microtemperature constants K1 and K3. The non-dimensional speed ρc2/μ of Rayleigh wave decreases very sharply with the increase of K2 in the range 0.01≤K2≤1W·m-1. This variation of non-dimensional speed against K2 is shown in Figure 2. The variation of the non-dimensional speed ρc2/μ of Rayleigh wave is similar against microtemperature constants K4, K5, and K6. It first decreases very sharply to its minimum value and thereafter it increases with the increase of K4, K5 and K6 values in the range 0≤K4, K5, or K6≤1W·m-3 as shown in Figures 3, 4, and 5. The non-dimensional speed of Rayleigh wave significantly depends on the frequency and microtemperature constants, as evident from Figures 1 to 5.
Variation of the nondimensional speed (ρc2/μ)1/2 of Rayleigh wave versus frequency χ.
Variation of the nondimensional speed (ρc2/μ)1/2 of Rayleigh wave versus microtemperature constant K2.
Variation of the nondimensional speed (ρc2/μ)1/2 of Rayleigh wave versus microtemperature constant K4.
Variation of the nondimensional speed (ρc2/μ)1/2 of Rayleigh wave versus microtemperature constant K5.
Variation of the nondimensional speed (ρc2/μ)1/2 of Rayleigh wave versus microtemperature constant K6.
7. Conclusion
The appropriate solutions of all the governing equations of thermoelastic medium with microtemperatures are applied at the boundary conditions at a thermally insulated free surface of a half-space to obtain the frequency equation of Rayleigh wave. From the frequency equation of Rayleigh wave, it is observed that the phase speed of Rayleigh wave depends on various material parameters including the microtemperature parameters. The dependence of numerical values of non-dimensional speed of propagation on the frequency and microtemperature parameters is shown graphically for a particular material representing the model. The problem though is theoretical but it can provide useful information for experimental researchers working in the field of geophysics and earthquake engineering and seismologist working in the field of mining tremors and drilling into the Earth crust. The study on wave propagation phenomenon in thermoelasticity with microtemperature is at its early stage. Recently, Steeb et al. [21] introduced the plane waves in such material. The present paper studied the propagation of Rayleigh wave in thermoelastic half-space with microtemperature. Based on theoretical results obtained by Steeb et al. [21] and in this paper, it is quiet early to predict possible specific applications of this phenomenon. The rapid advancement of MEMS/NEMS technology needs design and fabrication of microstructures. The possible applications of such studies may be in development of microtemperature sensors.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
CosseratE.CosseratF.1909Paris, FranceA. Hermann et FilsGrotR. A.Thermodynamics of a continuum with microstructure1969788018142-s2.0-0014550115ŘíhaP.On the microcontinuum model of heat conduction in materials with inner structure19761465295352-s2.0-0016892943IesanD.QuintanillaR.On a theory of thermoelasticity with microtemperatures20002331992152-s2.0-003435040210.1080/014957300280407IeşanD.Thermoelasticity of bodies with microstructure and microtemperatures20074425-26864886622-s2.0-3544894960310.1016/j.ijsolstr.2007.06.027IeşanD.QuintanillaR.On thermoelastic bodies with inner structure and microtemperatures2009354112232-s2.0-5934910435010.1016/j.jmaa.2008.12.017IesanD.On a theory of micromorphic elastic solids with microtemperatures20012487377522-s2.0-003561034010.1080/014957301300324882SvanadzeM.Boundary value problems of the theory of thermoelasticity with
microtemperatures20033188189SvanadzeM.Fundamental solutions of the equations of the theory of thermoelasticity with microtemperatures20042721511702-s2.0-1242332730CasasP. S.QuintanillaR.Exponential stability in thermoelasticity with microtemperatures2005431-233472-s2.0-2294448998610.1016/j.ijengsci.2004.09.004SvanadzeM.TracinàR.Representations of solutions in the theory of thermoelasticity with microtemperatures for microstretch solids20113421611782-s2.0-7865140736710.1080/01495739.2010.511946ScaliaA.SvanadzeM.Potential method in the linear theory of thermoelasticity with microtemperatures20093210102410422-s2.0-7034945163110.1080/01495730903103069AouadiM.Some theorems in the isotropic theory of microstretch thermoelasticity with microtemperatures20083176496622-s2.0-4584908391510.1080/01495730801981772IeşanD.Singular surfaces in the theory of thermoelasticity with microtemperatures20093212127912922-s2.0-7044941783210.1080/01495730903310631ScaliaA.SvanadzeM.TracinàR.Basic theorems in the equilibrium theory of thermoelasticity with microtemperatures20103387217532-s2.0-7795494321910.1080/01495739.2010.482348IeşanD.ScaliaA.Plane deformation of elastic bodies with microtemperatures20103776176212-s2.0-7864933156310.1016/j.mechrescom.2010.09.005YangL.HuangL.Propagation of singularities for Cauchy problems of semilinear thermoelastic systems with microtemperatures201028346156292-s2.0-7795249388210.1002/mana.200610821QuintanillaR.On growth and continuous dependence in thermoelasticity with microtemperatures20113499119222-s2.0-8005234385910.1080/01495739.2011.586278SvanadzeM.TracinàR.Representations of solutions in the theory of thermoelasticity with microtemperatures for microstretch solids20113421611782-s2.0-7865140736710.1080/01495739.2010.511946ChiritaS.CiarlettaM.D'ApiceC.On the theory of thermoelasticity with microtemperatures2013397349361SteebH.SinghJ.TomarS. K.Time harmonic waves in thermoelastic material with microtemperatures201348818