Jtam-A4.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 57, 1, pp. 155-166, Warsaw 2019 DOI: 10.15632/jtam-pl.57.1.155 FORCED VIBRATIONS OF A THERMOELASTIC DOUBLE POROUS MICROBEAM SUBJECTED TO A MOVING LOAD Rajneesh Kumar Kurukshetra University, Department of Mathematics, Kurukshetra, Haryana, India Richa Vohra H.P. University, Department of Mathematics and Statistics, Shimla, HP, India e-mail: richavhr88@gmail.com The present paper deals with forced vibrations of a homogeneous, isotropic thermoelastic double porous microbeam subjected to moving load, in context of Lord-Shulman theory of thermoelasticity with one relaxation time. The Laplace transform has been applied to obtain expressions for the axial displacement, lateral deflection, volume fraction field and temperature distribution. A numerical inversion technique has been used to recover the resulting quantities in the physical domain. Effects of velocity and time parameters are shown graphically by plotting axial displacement, lateral deflection, volume fraction field and temperature distribution against distance. Some particular cases are also deduced. Keywords: double porosity, thermoelasticity, Lord-Shulman theory,microbeam,moving load 1. Introduction Recently, dynamical analysis of engineering structures subjected tomoving loads has gained gre- at importance.Vehicle-bridge interactions are a vast area of interest in themoving load problem. Advances in transport technology and automobile engineering have resulted in high speeds and heaviness of vehicles and other moving bodies. As a result, corresponding structures have been subjected to vibration and dynamic stressmuch higher than ever before. The engineering struc- tures with moving loads often come out in buildings, bridges, railways and cranes. Beam type structures are widely used inmany fields like civil, mechanical and aerospace engineering.Many researchers have investigated dynamical behavior of beams on elastic foundations subjected to moving loads, especially in railway engineering. Themodern trend towards higher speeds in the railways has further intensified the research in order to accurately predict the vibration behavior of railway tracks. Pores or fractures can be observed in engineering structures due to reasons like erosion, corrosion, fatigue or accidents which affect the dynamic behavior of the entire structure to a considerable extent.This leads to thedevelopment of thedoubleporositymodelwhichhas its ap- plications in geophysics, rockmechanics andmany branches of engineering like civil engineering, chemical engineering and the petroleum industry. Biot (1941) proposed amodel for porousme- dia with single porosity. Later on Barenblatt et al. (1960) introduced amodel for porousmedia with a double porosity structure. The double porosity model consists of two coexisting degrees of porosity in which one corresponds to the porousmatrix and the other to the fissurematrix. Nunziato and Cowin (1979) developed a nonlinear theory of an elastic material with voids. Later, Cowin and Nunziato (1983) developed a theory of linear elastic materials with voids for mathematical study of the mechanical behavior of porous solids. In this theory, the skeletal materials are elastic, and interstices are void ofmaterial, hence an additional degree of freedom, 156 R. Kumar, R. Vohra the volume fraction of the void, is added. Iesan andQuintanilla (2014) derived a theory of ther- moelastic solids with a double porosity structure by using the theory developed by Nunziato and Cowin (1979). Darcy’s law was not used in developing that theory. So far, not much work has been done on the theory of thermoelasticity with the double porosity based on the model proposed by Iesan and Quintanilla (2014). Recent investigations have been started in the the- ory of thermoelasticity with double porosity which has a significant application in continuum mechanics. Kumar et al. (2015) applied the state space approach to a boundary value problem for thermoelastic materials with double porosity. The dynamic behavior of different isotropic structures subjected to moving loads has been investigated by many researchers. Olsson (1991) studied the dynamic problem of a simply sup- ported beam subjected to a constant force moving at a constant speed. The linear dynamic response of a simply supported uniform beam under a moving load of constant magnitude and velocitywas investigated byMichaltsos et al. (1996).Rao (2000) studied thedynamic responseof amulti-span Euler-Bernoulli beamdue tomoving loads.Mehri et al. (2009) presented the linear dynamic response of uniform beams with different boundary conditions under a moving load based on the Euler-Bernoulli beam theory. Sharma andGrover (2011) analysed a thermoelastic vibrations in micro-/nano-scale beam resonators with the presence of voids. Kargarnovin et al. (2012) studied the dynamic response of a delaminated composite beam under the action of a moving oscillatory mass. Esen (2015) investigated the transverse and longitudinal vibrations of a thin plate which carried a load moving along an arbitrary trajectory with variable velocity. Kumar (2016) studied the response of a thermoelastic beam due to the thermal source in the modified couple stress theory. Kaghazian et al. (2017) investigated free vibrations of a piezo- electric nanobeam using nonlocal elasticity theory. Zenkour (2017) studied the thermoelastic response of a microbeam embedded in visco-Pasternak’s medium based on GN-III model. In the present work, forced vibrations of a homogeneous, isotropic thermoelastic double porous microbeam, subjected to a moving load in the context of Lord-Shulman theory of ther- moelasticity has been investigated. The Laplace transform has been applied to find expressions for axial displacement, lateral deflection, volume fraction fields and temperature distribution. The resulting quantities are obtained in the physical domain by using a numerical inversion technique. Variations of the axial displacement, lateral deflection, volume fraction field and tem- perature distribution against the axial distance are depicted graphically to show the effect of the velocity parameter. Some particular cases have also been deduced. 2. Basic equations Following Iesan andQuintanilla (2014) as well as Lord and Shulman (1967); the field equations and the constitutive relation for a homogeneous isotropic thermoelastic material with a double porosity structure in the absence of body forces, extrinsic equilibrated body forces and heat sources can be written as µ∇2ui +(λ+µ)uj,ji + bϕ,i +dψ,i −βT,i = ρüi α∇2ϕ+ b1∇ 2ψ− bur,r −α1ϕ−α3ψ+γ1T = κ1ϕ̈ b1∇ 2ϕ+γ∇2ψ−dur,r −α3ϕ−α2ψ+γ2T = κ2ψ̈ (2.1) and ( 1+τ0 ∂ ∂t ) (βT0u̇j,j +γ1T0ϕ̇+γ2T0ψ̇+ρC ∗Ṫ)= K∗∇2T (2.2) tij = λerrδij +2µeij + bϕδij +dψδij −βTδij (2.3) Forced vibrations of a thermoelastic double porous microbeam... 157 where λ and µ are Lame’s constants, ρ is mass density, β = (3λ+2µ)αt, αt is linear thermal expansion, C∗ is specific heat at a constant strain, ui are displacement components, tij is the stress tensor, κ1 and κ2 are coefficients of equilibrated inertia, ϕ is the volume fraction field corresponding to pores, and ψ is the volume fraction field corresponding to fissures, K∗ is the coefficient of thermal conductivity, τ0 is the thermal relaxation time, κ1 and κ2 are coefficients of equilibrated inertia, and b,d,b1,γ,γ1,γ2 are constitutive coefficients, δij is Kronecker’s delta, T is the temperature changemeasured form the absolute temperature T0 (T0 6=0), a superposed dot represents differentiation with respect to time variable t. 3. Formulation of the problem We consider a homogeneous, isotropic thermoelastic double porous microbeam having dimen- sions: length L (0 ¬ x ¬ L), width a (−a/2 ¬ y ¬ a/2) and thickness h (−h/2 ¬ z ¬ h/2) in a Cartesian coordinate sytem Oxyz as shown in Fig. 1. The microbeam undergoes bending vibrations of a small amplitude about the x-axis such that the deflection is consistent with the linear Euler-Bernoulli theory. Therefore, the displacements can be written as u1 = u =−z ∂w ∂x u2 =0 u3 = w(x,t) (3.1) where w is the lateral deflection and u is the axial displacement. Fig. 1. Geometry of the beam The equation of motion for forced vibrations of the beam can be written as ∂2M ∂x2 +ρA ∂2w ∂t2 = F(x,t) (3.2) where A = ah is the cross-section area, M is the flexural moment of cross section of the micro- beam and F(x,t) is the appliedmoving load. By substituting Eqs. (2.3) and (3.1) into Eq.(3.2), we obtain the equation ofmotion for forcedvibrations of anEuler-Bernoulli thermoelastic double porousmicrobeam subjected to moving load as (λ+2µ)I ∂4w ∂x4 +ρA ∂2w ∂t2 − ∂2Mϕ ∂x2 − ∂2Mψ ∂x2 + ∂2MT ∂x2 = F(x,t) (3.3) where F = F0δ(x−vt) is the appliedmoving load, v is its velocity, δ is theDirac delta function, t is time in seconds, I = ah3/12 is the moment of inertia of the cross-section and Mϕ, Mψ are the volume fraction field moments, and MT is the thermal moment of the beam given by Mϕ = b h/2 ∫ −h/2 aϕz dz Mψ = d h/2 ∫ −h/2 aψz dz MT = β h/2 ∫ −h/2 aTz dz (3.4) 158 R. Kumar, R. Vohra Equations (2.1)2,3 and (2.2) with the help of Eq. (3.1) can be written as α (∂2ϕ ∂x2 + ∂2ϕ ∂z2 ) + b1 (∂2ψ ∂x2 + ∂2ψ ∂z2 ) + bz ∂2w ∂x2 −α1ϕ−α3ψ+γ1T = κ1 ∂2ϕ ∂t2 b1 (∂2ϕ ∂x2 + ∂2ϕ ∂z2 ) +γ (∂2ψ ∂x2 + ∂2ψ ∂z2 ) +dz ∂2w ∂x2 −α3ϕ−α2ψ+γ2T = κ2 ∂2ψ ∂t2 K∗ (∂2T ∂x2 + ∂2T ∂z2 ) = ( 1+ τ0 ∂ ∂t )[ −βT0z ∂ ∂t (∂2w ∂x2 ) +γ1T0ϕ̇+γ2T0ψ̇+ρC ∗Ṫ ] (3.5) 4. Solution of the problem For the present microbeam, we assume that there is no flow of heat and the volume fraction fields across the surfaces (z =±h/2) so that ∂T/∂z = ∂ϕ/∂z = ∂ψ/∂z =0 at z =±h/2. For a very thin beam, assuming that the volume fraction fields and temperature increment in terms of sin(πz/h) function along the thickness direction one obtains ϕ(x,z,t) = Φ(x,t)sin πz h ψ(x,z,t) = Ψ(x,t)sin πz h T(x,z,t)= Θ(x,t)sin πz h (4.1) Introducing non-dimensional variables as x′ = 1 L x u′ = 1 L u t′x = tx E Φ′ = L α Φ Ψ ′ = L α Ψ Θ′ = β E Θ t′ = c1 L t τ ′0 = c1 L τ0 F ′ 0 = L2 ah2(λ+2µ) F0 (4.2) where c21 =(λ+2µ)/ρ and E = µ(3λ+2µ)/(λ+µ) is Young’s modulus. Making use of Eqs. (4.1) in Eq. (3.3) andwith the aid of Eqs. (4.2), yields (after suppressing primes) ∂4w ∂x4 +a1 ∂2w ∂t2 −a2 ∂2Φ ∂x2 −a3 ∂2Ψ ∂x2 +a4 ∂2Θ ∂x2 = F0δ(x−vt) (4.3) On multiplying Eqs. (3.5) by z and integrating them with respect to z from −h/2 to h/2 and after using Eq.(4.2), we obtain (suppressing primes for convenience) a5 ∂2Φ ∂x2 −a6Φ+a7 ∂2Ψ ∂x2 −a8Ψ +a9 ∂2w ∂x2 −a10Φ−a11Ψ +a12Θ = ∂2Φ ∂t2 a13 ∂2Φ ∂x2 −a14Φ+a15 ∂2Ψ ∂x2 −a16Ψ +a17 ∂2w ∂x2 −a18Φ−a19Ψ +a20Θ = ∂2Ψ ∂t2 (4.4) and ∂2Θ ∂x2 −a21Θ = ( 1+ τ0 ∂ ∂t )[ a22 ∂ ∂t (∂2w ∂x2 ) +a23 ∂Φ ∂t +a24 ∂Ψ ∂t +a25 ∂Θ ∂t ] (4.5) Forced vibrations of a thermoelastic double porous microbeam... 159 where a1 = ρc21L I(λ+2µ) a2 = 2bα Iπ2(λ+2µ)L a3 = 2dα Iπ2(λ+2µ)L a4 = 2E Iπ2(λ+2µ) a5 = α κ1c 2 1 a6 = απ2L2 κ1c 2 1h 2 a7 = b1 κ1c 2 1 a8 = b1π 2L2 κ1c 2 1h 2 a9 = bhπ2L2 24ακ1c 2 1 a10 = α1L 2 κ1c 2 1 a11 = α3L 2 κ1c 2 1 a12 = γ1EL 3 αβκ1c 2 1 a13 = b1 κ2c 2 1 a14 = b1π 2L2 κ2c 2 1h 2 a15 = γ κ2c 2 1 a16 = γπ2L2 κ2c 2 1h 2 a17 = dhπ2L2 24ακ2c 2 1 a18 = α3L 2 κ2c 2 1 a19 = α2L 2 κ2c 2 1 a20 = γ2EL 3 αβκ2c 2 1 a21 = π2L2 h2 a22 =− β2T0hc1π 2 24EK∗ a23 = αβT0γ1c1 EK∗ a24 = αβT0γ2c1 EK∗ a25 = ρC∗c1L K∗ The initial conditions of the problem are assumed to be homogeneous and are taken as w(x,t) ∣ ∣ ∣ t=0 = ∂w(x,t) ∂t ∣ ∣ ∣ t=0 = Φ(x,t) ∣ ∣ ∣ t=0 = ∂Φ(x,t) ∂t ∣ ∣ ∣ t=0 = Ψ(x,t) ∣ ∣ ∣ t=0 = ∂Ψ(x,t) ∂t ∣ ∣ ∣ ∣ ∣ t=0 =0 Θ(x,t) ∣ ∣ ∣ t=0 = ∂Θ(x,t) ∂t ∣ ∣ ∣ ∣ ∣ t=0 =0 (4.6) These initial conditions are supplemented by considering that the two ends of the microbe- am are clamped and remain at zero increment of the volume fraction fields and temperature. Mathematically, it can be written as w(x,t) ∣ ∣ ∣ x=0,L = ∂w(x,t) ∂x ∣ ∣ ∣ ∣ ∣ x=0,L =0 Φ(x,t) ∣ ∣ ∣ x=0,L =0 Ψ(x,t) ∣ ∣ ∣ x=0,L =0 Θ(x,t) ∣ ∣ ∣ x=0,L =0 (4.7) 5. Solution in the Laplace transform domain Applying the Laplace transform defined by f(s)= L[f(t)] = ∞ ∫ 0 f(t)e−st dt (5.1) to Eqs. (4.3)-(4.5) under initial conditions (4.6), and after some simplifications, we obtain ( d10 dx10 +B1 d8 dx8 +B2 d6 dx6 +B3 d4 dx4 +B4 d2 dx2 +B5 ) (w,Φ,Ψ,Θ)= (f1,f2,f3,f4)e − s v x (5.2) where B1, B2, B3, B4, B5, f1, f2, f3, f4, are given in Appendix I. 160 R. Kumar, R. Vohra The solution to system of Eqs. (5.2), in the Laplace transform domain, can be written as w = H1e − s v x + 5 ∑ i=1 (Die −mix +Di+5e mix) Φ = H2e − s v x + 5 ∑ i=1 g1i(Die −mix +Di+5e mix) Ψ = H3e − s v x + 5 ∑ i=1 g2i(Die −mix +Di+5e mix) Θ = H4e − s v x + 5 ∑ i=1 g3i(Die −mix +Di+5e mix) (5.3) where Hi = fiν 10 s10+B1s8ν2+B2s6ν4+B3s4ν6+B4s2ν8+B5ν10 i =1,2,3,4 and g1i, g2i, g3i; (i =1,2,3,4,5) are given in Appendix II. Here±mi, i =1,2, . . . ,5 are the roots of the characteristic equation m10+B1m 8+B2m 6+B3m 4+B4m 2+B5 =0 Therefore, the corresponding expressions for the axial displacement in the Laplace transform domain can be written as u =−z dw dx =−z [ 5 ∑ i=1 (−miDie −mix +miDi+5e mix)− s v H1e − s v x ] (5.4) Boundary conditions (4.7) in the Laplace transform domain take the form as w(x,s) ∣ ∣ ∣ x=0,L = dw(x,s) ∂x ∣ ∣ ∣ ∣ ∣ x=0,L =0 Φ(x,s) ∣ ∣ ∣ x=0,L =0 Ψ(x,s) ∣ ∣ ∣ x=0,L =0 Θ(x,s) ∣ ∣ ∣ x=0,L =0 (5.5) By substituting Eqs. (5.3) into boundary conditions (5.5), we obtain a system of ten linear equations in the matrix form as AD=E (5.6) where A=                   1 1 1 1 1 1 1 1 1 1 b1 b2 b3 b4 b5 a1 a2 a3 a4 a5 −m1 −m2 −m3 −m4 −m5 m1 m2 m3 m4 m5 −m1b1 −m2b2 −λ3b3 −m4b4 −m5b5 m1a1 m2a2 m3a3 m4a4 m5a5 g11 g12 g13 g14 g15 g11 g12 g13 g14 g15 g11b1 g12b2 g13b3 g14b4 g15b5 g11a1 g12a2 g13a3 g14a4 g15a5 g21 g22 g23 g24 g25 g21 g22 g23 g24 g25 g21b1 g22b2 g23b3 g24b4 g25b5 g21a1 g22a2 g23a3 g24a4 g25a5 g31 g32 g33 g34 g35 g31 g32 g33 g34 g35 g31b1 g32b2 g33b3 g34b4 g35b5 g31a1 g32a2 g33a3 g34a4 g35a5                   ai =e miL bi =e −miL i =1,2,3,4,5 D= [ D1 D2 D3 D4 D5 D6 D7 D8 D9 D10 ]T E= [ E1 E2 E3 E4 E5 E6 E7 E8 E9 E10 ]T Forced vibrations of a thermoelastic double porous microbeam... 161 and E1 =−H1 E2 =−H1e − s v L E3 = s v H1 E4 = s v H1e − s v L E5 =−H2 E6 =−H2e − s v L E7 =−H3 E8 =−H3e − s v L E9 =−H4 E10 =−H4e − s v L By solving the above system of equations (5.6), we obtain the values of unknownparameters Di, i =1,2, . . . ,10. This completes the solution of the problem in the Laplace transform domain. In order to determine the axial displacement, lateral deflection, volume fraction field and temperature distribution in the physical domain, we will adopt a numerical inversion method given by Honig and Hirdes (1984). In this method, the Laplace domain f(s) can be inverted to the time domain f(t) as f(t)= 1 t1 exp(Ωt) [1 2 f(Ω)+Re N ∑ k=1 f ( Ω + ikπ t1 ) exp (ikπt t1 )] 0 < t1 < 2t where Re is the real part and i is the imaginary unit. The value of N is chosen sufficiently large and it represents the number in terms of the truncated Fourier series such that f(t)= exp(Ωt)Re [ f ( Ω + iNπ t1 ) exp (iNπt t1 )] ¬ ε1 where ε1 is a prescribed small positive number. Also, the value of Ω should satisfy the relation Ωt ≃ 4.7 for faster convergence, Tzou (1996). Particular cases i) If τ0 =0 in equations (5.6), it yields correspondingexpressions for the thermoelastic double porousmicrobeam in context of the coupled theory (CT) of thermoelasticity. ii) If b1 = α3 = γ = α2 = γ2 = d → 0 in equations (5.6), we obtain the corresponding expressions for the thermoelastic single porousmicrobeam (thermoelasticmicrobeamwith voids). 6. Numerical results and discussion Numerical computations have been done for a copper like material microbeam. The material parameters are taken as inKumar et al. (2015): λ =7.76·1010Nm−2, C∗ =3.831·103m2s−2K−1, µ = 3.86 · 1010Nm−2, K∗ = 3.86 · 103Ns−1K−1, T0 = 298K, ρ = 8.954 · 10 3Kgm−3, αt =1.78·10 −5K−1,α2 =2.4·10 10Nm−2,α3 =2.5·10 10Nm−2, γ =1.1·10−5N,α =1.3·10−5N, γ2 = 0.219 ·10 5Nm−2, κ1 = 0.1456 ·10 −12Nm−2s2, b = 0.9 ·1010Nm−2, α1 = 2.3 ·10 10Nm−2, κ2 =0.1546 ·10 −12Nm−2s2, τ0 =0.01s. The aspect ratio of the beam is fixed as L/h = 10, a/h = 0.5, z = h/6. When h is varied, L and a change accordingly with h. For the microscale beam, we take the range of the beam length L =(1−10) ·10−6m. The plots are prepared by using the dimensionless variables for a wide range of the beam length when, unless otherwise stated, L =4.0, a = h/2 and z = h/6. The software MATLAB has been used to find the values of axial displacement, lateral de- flection, volume fraction field and temperature distribution. Variations of these quantities with respect to the axial distance have been shown in Figs. 2-5 to indicate the effects of velocity and time parameters. In Figs. 2 and 3, solid line, small and big dashes lines correspond to the values of velocity v = 1.0, 2.0 and 3.0, respectively with the fixed value of time t = 0.15, whereas in 162 R. Kumar, R. Vohra Figs. 4 and 5, solid line, small and big dashes lines correspond to the values of time t = 0.15, 0.175 and 0.2, respectively with the fixed value of velocity v =1.0. Effect of the velocity parameter InFig. 2a, it is noticed that the value of axial displacement u initially decreases in 0 < x < 1, then increases in 1¬ x < 2.8 and again decreases slowly and steadily in the remaining region. It is also found that the magnitude of u decreases with an increase in the value of velocity v near the source application point while the trend gets reversed asmoving away from the source. Figure 2b depicts that the lateral deflection w decreases in 0 < x < 1, increases in 1¬ x ¬ 2.2 and then becomes stationary as x ­ 2.2. The amplitude of variation is higher near the point of application of the source while as moving away from the source, the values become almost stationary for all the values of velocity v. Figure 3a shows that the volume fraction field ϕ initially decreases sharply in 0 < x < 1, then increases abruptly in 1 ¬ x < 3.5 and then decreases slowly and steadily in the remaining region. Also, the magnitude of ϕ decreases as v increases near the source application point while the trend gets reversed away from the source. From Fig. 3b, it is clear that the value of temperature distribution T initially increases in the region 0 < x < 1 and decreases monotonically as x ­ 1. It is also evident that the magnitude of T increases with an increase in velocity v. Fig. 2. (a) Axial displacement u and (b) lateral deflection w versus axial distance x Fig. 3. (a) Volume fraction field ϕ and (a) temperature distribution T versus axial distance x Forced vibrations of a thermoelastic double porous microbeam... 163 Effect of the time parameter InFig. 4a, it is noticed that the value of axial displacement u initially decreases in 0 < x < 1, then increases in 1¬ x < 2.8andagain decreases in the remaining region.Also, themagnitudeof u increases with an increase in the value of time t. Figure 4b depicts that the lateral deflection w decreases in 0 < x < 1and then increases afterwards asx ­ 1. It is found that as time t increases, the magnitude of w decreases. Figure 5a shows that the volume fraction field ϕ is oscillatory in nature. The value of ϕ initially decreases in the range 0 < x < 1, then increases in 1¬ x < 3.5 and thendecreases in 3.5¬ x < 5.3 andagain starts increasing slowly asx ­ 5.3.Themagnitude value of ϕ also increases with an increase in the value of time t. From Fig. 5b, it is clear that the value of temperature distribution T initially increases in the region 0 < x < 1 and decreases monotonically as x ­ 1. It is also found that the magnitude of T decreases with an increase in the value of time t. Fig. 4. (a) Axial displacement u and (b) lateral deflection w versus axial distance x Fig. 5. (a) Volume fraction field ϕ and (b) temperature distribution T versus axial distance x 164 R. Kumar, R. Vohra 7. Conclusions In the present work, forced vibrations of an Euler-Bernoulli thermoelastic double porousmicro- beam, in context of Lord-Shulman theory of thermoelasticity, subjected tomoving load has been investigated. Effects of velocity and time parameters are shown graphically on axial displace- ment, lateral deflection, volume fraction field and temperature distribution. All field quantities are observed to be very sensitive towards the velocity as well as time parameters. • It is observed that the amplitude of variation is higher near the application point of the source while the values become almost stationary as moving away from the source due to the effect of velocity. The values of axial displacement u and volume fraction field ϕ decreasewith an increase in the velocity v near the application point of the source,whereas an opposite trend of variation is noticed away from the source. The values of lateral deflection w also decrease as the velocity parameter increases near the point of application of the source and becomes stationary as moving away from the source. Themagnitude of temperature distribution T gets greater with an increase in the velocity parameter. • Due to the effect of time parameter, the values of axial displacement u and volume frac- tion field ϕ increase with an increase in the value of time t, whereas an opposite trend and behavior of variation is observed in the case of lateral deflection w and temperature distribution T , i.e. the values of w and T decrease as there is an increase in the value of time t. This type of study is useful due to its physical application inmany fields of engineering like civil, mechanical, aerospace and industrial sectors. The results obtained in this investigation should prove to be beneficial for researchers working on the theory of thermoelasticity with the double porosity structure.The introduction of the double porous parameter to the thermoelastic medium represents a more realistic model for further studies. Appendix I a26 = s(1+ τ0s) a27 =−s(1+τ0s)a22 a28 =−s(1+ τ0s)a23 a29 =−s(1+ τ0s)a24 a30 =−[a21+s(1+ τ0s)a25] n1 =−(a6+a10+s 2) n2 =−(a8+a11) n3 =−(a18+a14) n4 =−(a16+a19+s 2) r1 = a5a15−a7a13 r2 = a5(a15a30+n4)−a13n2+n1a15−a7(a13a30+n3) r3 = n1(a15a30+n4)+a5(n4a30−a20a29)−a7(n3a30−a20a28)−n2(a13a30+n3) +a12(a13a29−a15a28) r4 = n1(n4a30−a20a29)+a12(n3a29−n4a28)+n2(a20a28−n3a30) r5 = a9a15−a7a17 r6 = a9(a15a30+n4)−a7(a17a30−a20a27)−n2a17−a12a15a27 r7 = a9(n4a30−a20a29)+a12(a17a29−n4a27)−n2(a17a30−a20a27) r8 = a9a13−a5a17 r9 = a9(a13a30+n3)−n1a17−a5(a17a30−a20a27)−a12a10a27 r10 = a9(n3a30−a20a28)−n1(a30a17−a27a20)+a12(a17a28−n3a27) r11 = a27(a5a12+a7a13) r12 = a9(a13a29−a15a28)+a5(n4a27−a17a29)+a27(n1a15+g3a13)−a7(a17a28−n3a27) r13 = a9(n3a29−n4a28)−n1(a17a29−n4a27)−n2(a17a28−n3a27) Forced vibrations of a thermoelastic double porous microbeam... 165 B1 = r2+a2r5−a3r8−a4r11 r1 B2 = a1r1s 2+a2r6−a3r9−a4r12+r3 r1 B3 = a1r3s 2+a2r7−a3r10−a4r13+r4 r1 B4 = a1r3s 2 r1 B5 = a1r4s 2 r1 f1 = 1 v7 F0(r1s 6+r2s 4v2+r3s 2v4+r4v 6) f2 =− 1 v7 F0(r5s 6+r6s 4v2+r7v 4s2) f3 = 1 v7 F0(r8s 6+r9s 4v2+r10v 4s2) f4 =− 1 v7 F0(r11s 6+r12s 4v2+r13v 4s2) Appendix II g1i =− r5m 6 i +r6m 4 i +r7m 2 i r1m 6 i +r2m 4 i +r3m 2 i +r4 g2i = r8m 6 i +r9m 4 i +r10m 2 i r1m 6 i +r2m 4 i +r3m 2 i +r4 g3i =− r11m 6 i +r12m 4 i +r13m 2 i r1m 6 i +r2m 4 i +r3m 2 i +r4 i =1,2, . . . ,5 References 1. Barenblatt G.I., Zheltov I.P., Kochina I.N., 1960, Basic concept in the theory of seepage of homogeneous liquids in fissured rocks (strata), Journal of Applied Mathematics and Mechanics, 24, 1286-1303 2. Biot M.A., 1941, General theory of three-dimensional consolidation, Journal of Applied Physics, 12, 155-164 3. Chang T.P., Liu Y.N., 1996, Dynamic finite element analysis of a nonlinear beam subjected to a moving load, International Journal of Solids and Structures, 33, 12, 1673-1688 4. Cowin S.C., Nunziato J.W., 1983, Linear elasticmaterialswith voids, Journal of Elasticity, 13, 125-147 5. Esen I., 2015, A new FEM procedure for transverse and longitudinal vibration analysis of thin rectangular plates subjected to a variable velocitymoving load along an arbitrary trajectory,Latin American Journal of Solids and Structures, 12, 808-830 6. Honig G., Hirdes U., 1984, A method for the numerical inversion of the Laplace transforms, Journal of Computational and Applied Mathematics, 10, 113-132 7. Iesan D., Quintanilla R., 2014, On a theory of thermoelastic materials with a double porosity structure, Journal of Thermal Stresses, 37, 1017-1036 8. Kaghazian A., Hajnayeb A., Foruzande H., 2017, Free vibration analysis of a piezoelectric nanobeam using nonlocal elasticity theory, Structural Engineering and Mechanics, 61, 5, 617-624 9. KargarnovinM.H.,AhmadianM.T.,TalookolaeiR.A.J., 2012,Dynamics of a delaminated Timoshenko beam subjected to a moving oscillatory mass,Mechanics Based Design of Structures and Machines, 40, 2, 218-240 10. Khalili N., 2003, Coupling effects in double porositymedia with deformablematrix,Geophysical Research Letters, 30, 22, 2153, DOI: 10.1029/2003GL018544 11. Kumar R., 2016,Response of thermoelastic beamdue to thermal source inmodified couple stress theory,Computational Methods in Science and Technology, 22, 2, 95-101 12. Kumar R., Vohra R., Gorla M.G., 2015, State space approach to boundary value problem for thermoelastic material with double porosity, Applied Mathematics and Computation, 271, 1038-1052 166 R. Kumar, R. Vohra 13. Lord H., Shulman Y., 1967, A generalized dynamical theory of thermoelasticity, Journal of Mechanics and Physics of Solids, 15, 299-309 14. Mehri B., Davar A., Rahmani O., 2009, Dynamic Green function solution of beams under a moving load with different boundary conditions, Scientia Iranica, 16, 3, 273-279 15. MichaltsosG., SopianopouloD.,KounadisA.N., 1996,The effect ofmovingmass and other parameters on the dynamic response of a simply supported beam, Journal of Sound and Vibration, 191, 357-362 16. Nunziato J.W., Cowin S.C., 1979, A nonlinear theory of elastic materials with voids, Archive of Rational Mechanics and Analysis, 72, 175-201 17. Olsson M., 1991, On the fundamental moving load problem, Journal of Sound and Vibration, 145, 2, 299-307 18. Rao G.W., 2000, Linear dynamics of an elastic beam under moving loads, Journal of Vibration and Acoustics, 122, 3, 281-289 19. Sharma J.N.,GroverD., 2011,ThermoelasticVibrations inmicro-/nano-scalebeam resonators with voids, Journal of Sound and Vibration, 330, 2964-2977 20. Sherief H., Saleh H., 2005, A half space problem in the theory of generalized thermoelastic diffusion, International Journal of Solids and Structures, 42, 4484-4493 21. Sun Y., Fang D, Saka M , Soh A.K., 2008, Laser-induced vibrations of micro-beams under different boundary conditions, International Journal of Solids and Structures, 45, 1993-2013 22. Tzou D., 1996,Macro-to-Micro Heat transfer, Taylor & Francis,Washington DC 23. Zenkour A.M., 2017, Thermoelastic response of a microbeam embedded in visco-Pasternak’s medium based onGN-III model, Journal of Thermal Stresses, 40, 2, 198-210 Manuscript received May 9, 2017; accepted for print August 23, 2018