THERMAL EFFECT ON FREE VIBRATION AND BUCKLING FACTA UNIVERSITATIS Series: Mechanical Engineering Vol. 15, N o 1, 2017, pp. 45 - 62 DOI: 10.22190/FUME161115007S © 2017 by University of Niš, Serbia | Creative Commons Licence: CC BY-NC-ND Original scientific paper THERMAL EFFECT ON FREE VIBRATION AND BUCKLING OF A DOUBLE-MICROBEAM SYSTEM UDC 620.179.13+620.174]:534.1 Marija Stamenković Atanasov 1 , Danilo Karličić 1 , Predrag Kozić 2 , Goran Janevski 2 1 Mathematical Institute of the Serbian Academy of Sciences and Arts, Serbia 2 University of Niš, Department of Mechanical Engineering, Serbia Abstract. The paper investigates the problem of free vibration and buckling of an Euler- Bernoulli double-microbeam system (EBDMBS) under the compressive axial loading with a temperature change effect. The system is composed of two identical, parallel simply- supported microbeams which are continuously joined by the Pasternak’s elastic layer. Analytical expressions for the critical buckling load, critical buckling temperature, natural frequencies and frequencies of transverse vibration of the EBDMBS represented by the ratios are derived and validated by the results found in the literature. Also analytical expressions are obtained for various buckling states and vibration-phase of the EBDMBS. The temperature change effect is assumed to have an influence on both the microbeams. The length scale parameter, temperature change effect, critical buckling load, thickness/material parameter, Pasternak’s parameter and Poisson’s effect are discussed in detail. Also, as a clearer display of the thermo-mechanical response of EBDMBS, the paper introduces a critical scale load ratio of the modified and the local critical buckling loads in low- temperature environs. Numerical results show that the critical buckling temperatures for classical theories are always higher than the critical buckling temperature for MCST systems. Key Words: Thermal Effect, Double-microbeam System, Critical Buckling Load, Pasternak’s Parameter, Poisson’s Effect 1. INTRODUCTION Micro and nano structures became an object of interest in modern science and technology just after their invention. They possess important mechanical, electrical and thermal performances that are higher than conventional structural materials. Using micro/nano Received November 15, 2016 / Accepted February 23, 2017 Corresponding author: Goran Janevski University of Niš, Department of Mechanical Engineering, A. Medvedeva 14, 18000 Nis, Serbia E-mail: gocky.jane@gmail.com 46 M. STAMENKOVIŠ ATANASOV, D. KARLIŢIŠ, P. KOZIŠ, G. JANEVSKI structures in a high temperature environment leads to certain changes in stiffness. Recently, the vibration and buckling studies of beams with the microstructure effect have been increasingly present in the scientific community. Researchers are motivated to develop theories such as the modified couple stress theory (MCST) which contains the material length scale parameter; also, they are able to describe size effects on the micro-scale. On the other hand, the classical continuum mechanics theories neither contain any internal material length scale parameter nor are they able to describe these effects. The structural elements such as beams, plates, and membranes in the micro or nano length scale are frequently used as components in micro/nano electromechanical systems (MEMS/NEMS). With the rapid development of technology, functionally graded (FG) beams and plates are often used in MEMS/NEMS, such as the components in the shape of memory thin films alloy with a global thickness in the micro/nano scale, atomic force microscopes (AFMs), and electrically actuated MEMS devices [1–8]. As opposed to the strain theories which were introduced by Mindlin and Eshel [9], with five constants besides the Lamé constants, Lam et al. [1] presented a modified theory consisting of only three non-classical constants. Wang et al. [10] used the above theory to analyze the behavior of micro beams considering Euler–Bernoulli and Timoshenko beam theories. Analysis of bending and buckling of a thin beam were presented by Lazopoulos and Lazopoulos [11]. These results imply that the gradient coefficient has a significant effect on the buckling load while the surface effect of the energy is negligible. The modified couple stress theory has been used by many authors just as it has been mostly applied to micro beams. The non-local Bernoulli–Euler beam model was proposed by Peddieson et al. [12], using a constitutive equation after Eringen et al. [13] which contains two additional material constants. The non-local theories for the Bernoulli–Euler, Timoshenko, Reddy, and Levinson beams were developed by Reddy [14] in a unified way using the Hamilton principle and the non-local constitutive relation of Eringen et al. [13]. Park and Gao [18] used a modified couple stress theory with an Euler-Bernoulli formulation for the bending analysis of cantilever beams. Ma et al. [19] and Reddy et al. [17] developed a modified Timoshenko beam theory and investigation of the bending and free vibration of simply- supported beams using the Navier solution process and finite elements method. The buckling analysis of functionally graded micro beams based on the MCST is presented in Nateghi et al. [20]. The free vibration of a single-layered graphene sheet resting on an elastic matrix as a Pasternak foundation model explored by using the modified couple stress theory is presented by Bekir and Civalek [21]. M. Simsek and Reddy [22] developed a united higher order beam theory for a functionally graded (FG) microbeam embedded in an elastic Pasternak medium using the modified couple stress theory. In the paper of Hendou and Mohammadi [23], an Euler–Bernoulli model has been used for vibration analysis of micro-beams with large transverse deflection where thermoelastic damping is considered to be the main damping mechanism and displayed as imaginary stiffness into the equation of motion by evaluating the temperature profile as a function of lateral displacement. Free vibration and buckling of microbeams with the temperature change effect is presented by Ke et al. [24]. The finding that the thermal effect on the fundamental frequency and critical buckling load is very low when the thickness of the microbeam has a similar value to that of the material length scale parameter and that it becomes significant when the thickness of the microbeam becomes larger, was of particular interest, and the present paper considers what may happen with a double microbeam system. Scientists are trying to comprehend the vibration behavior of micro/nano structures with the effect of temperature changes. Primarily motivated by the last few studies [22, Thermal Effect on the Free Vibration and Buckling of a Double-Microbeam System 47 23, 24], this paper analyzes the free vibration and buckling behavior with the temperature change effect of the Euler-Bernoulli double-microbeam system (EBDMBS). To solve the higher-order main equations of the EBDMBS we use the Bernoulli–Fourier method. The system is composed of two identical and parallel, simply-supported beams which are continuously joined by the Pasternak’s elastic layer. It is assumed that the temperature change effect has an impact on both the microbeams. The length scale parameter, temperature change effect, critical buckling load, thickness/material parameter, Pasternak’s parameter and Poisson’s effect are discussed in detail. The paper presents the impact of different above-mentioned parameters on the natural frequency, frequency under compressive axial loading, critical buckling load and critical temperature of EBDMBS with thermal effect. Also, results for various buckling state and vibration-phase of the EBDMBS are obtained. The vibration phases include out-of-phase and in-phase modes of vibration. In order to verify the present study, a comparison of the thermal effect on the dimensionless natural frequency of the system for the first three modes with the results found in the literature is given in the tabular form. Because of the strong coupling between mechanical and electrical phenomenon in electromechanical microdevices, there is a growing need for results with temperature effect since they can give contribution to the making of modern microsensors. The ability of the MEMS device is precision and sensitivity without the need for any cumbersome electrical components. 2. FORMULATION On the basis of the MCST, we discuss the oscillatory system of two parallel Euler– Bernoulli microbeams which are continuously joined by the Pasternak elastic layer under the influence of axial loading including the temperature change effect (see Fig. 1).The Pasternak foundation assumes the presence of shear interaction among the spring elements which is achieved by connecting the ends of the springs to a beam that only undergoes transverse shear deformation, see [25]. The load–deflection relationship is obtained by taking into account the vertical equilibrium of a shear layer. The pressure–deflection relationship is given by ,2,1, 2 0 2 0     i x w Gkwp i i (1) where G is the modulus of a shear Pasternak's layer and K is the stiffness modulus of a Winkler elastic layer. Beams are continuously connected with the Winkler elastic layer which represents an idealized medium formed of system springs. Both the microbeams are rectangular and have the same length L, thickness h, width b. Fig. 1 Double-microbeam system coupled by the Pasternak's layer 48 M. STAMENKOVIŠ ATANASOV, D. KARLIŢIŠ, P. KOZIŠ, G. JANEVSKI The beams are simply supported at the ends and under the effect of the axial compressive load with the temperature change effect. 2.1. Introduction of the modified coupled stress theory The MCST was developed from the classical couple stress theory, which was well grounded by Mindlin [26], Mindlin and Tiersten [27], Toupin [28] and Koiter [29]. This theory, suggested by Yang et al.(2002), holds that energy density is a function of strain as well as curvature. According to the modified couple stress theory, Yang [30], Park and Gao [18], Ma et al. [19] and Reddy [17], strain energy Us in an isotropic linear elastic material occupying area Ω with a volume element V, and, can be written as 1 ( : : ) , 2 s U dV      m (2) where σ and ε are the Cauchy stress tensor and strain tensor, respectively, m is the deviatory part of the couple stress tensor, and χ is the symmetric part of the curvature tensor. These tensors correspond to the geometrical equations: 1 1 [ ( ) ] , [ ( ) ] , , 2 2 T T x y z                        u u (3) where  is the Nabla operator, u is the displacement vector. The rotation vector and constitutive equations are defined by , 1 curl 2  u (4)   .2, 2   l2tr  mI  (5) Where material length scale parameter l has the dimension of length which is mathematically the square of the ratio of the curvature modulus to the shear modulus and is physically regarded as a material property measuring the effect of couple stress, Mindlin [26], μ and λ are the Lamé constants that are given as: and . (1 )(1 2 ) 2(1 ) E            (6) In order to implement the linear constitutive relations presented in Eq. (6) the microbeam material should be made homogeneous, isotropic and linearly elastic. 2.2. Mathematical model of the double-microbeam system Based on the Euler-Bernoulli beam theory, axial displacements u(x,z,t) and transverse displacements of any point of the beam, w(x,z,t) are given by Reddy [14] as 0 0 ( , ) ( ) , ( ) 0 , ( ) ( , ), W x t u x, z,t z v x, z,t w x, z,t W x t x       (7) Thermal Effect on the Free Vibration and Buckling of a Double-Microbeam System 49 where W0(x,t) is the midplane displacement. From vector Eqs. (3)-(5) and displacements (7) it follows that 2 0 0 2 ( , ) ( , ) , 0 , , 0, xx yy zz xz yz xy y x z W x t W x t z xx                        (8) 2 0 2 ( , )1 , 0 , 2 xy yy zz xz yz xy W x t x                (9) 2 2 0 2 ( , ) , 0 , xy xx yy zz xz yz W x t m μl m m m m m x          (10) 2 0 2 ( , )(1 ) , 0 , (1 )(1 2 ) xx yy zz xz yz xy W x tE z T x                           (11) where E is Young’s modulus, α is the coefficient of thermal expansion, ν is Poisson's ratio and ΔT=TT0 is the temperature change with a respect to reference temperature T0 and assuming no shear strains are created by temperature change. In this study, the equilibrium equations are derived by the principle of total potential energy [15, 16]. From Eqs (2) and (8) - (11), the variation of strain energy in the double- beam system can be determined as 2 2 0 0 2 2 0 ( 2 ) , ( 1, 2), L i i s xx xx xy xy xi xyi W W U m dV M Y dx i x x                       (12) where Mxi and Yxyi are the stress resultant moments and couple moments for the first and second microbeam, respectively, defined as ., 2 0 2 2 2 0 2 x W lAdAmY x W DdAzM i xzi A xyxyi i xxi A xxxi ii         (13) The stiffness components in Eqs. (13) are defined as [22] 2(1 ) { , } {1, } , , ( 1, 2). (1 )(1 2 ) i i i xx xx i xz i A A E A D z dA A dA i             (14) Using the displacement field components given in Eq. (7), we obtain the variation of kinetic energy in the form ),2,1(,, 0 0 02 0 2 0      iAdAmdxW t W mK iii A ii L i i ie i  (15) where ρi is the mass density for the first and second microbeam. The first variation of the additional strain energy caused by the elastic medium is written by 01 01 02 02 01 02 01 02 0 ( ) ( ) , L ad W W W W U K W W W W G G dx x x x x                      (16) 50 M. STAMENKOVIŠ ATANASOV, D. KARLIŢIŠ, P. KOZIŠ, G. JANEVSKI where K and G are the spring constants of the Winkler and Pasternak elastic medium, respectively. The first variation of the work done by axial forces Fxi=Fmi+FT, (i=1,2) can be given as 0 0 0 ( ) , ( 1, 2), L i i ext mi T W W W F F dx i x x           (17) where FT=AxxαΔT is the axial force due to the influence of the temperature change and Fmi, (i=1,2) is the axial forces due to the mechanical loading for the first and second microbeams. The main equation and the boundary conditions can be derived by the Hamilton principles as follows 0 [ ( )] 0. t e s ad ext K U U W dt       (18) If we substitute the expressions for δUs, δKe, δUad and δWext from Eqs. (12), (15), (16) and (17) into Eq. (18) and after integrating by parts and then collecting the coefficients of δW01 and δW02, the equations of motion of the double microbeam system are obtained in the form ),()(: 02012 01 2 2 01 2 12 1 2 1 2 2 01 2 0101 WWK x W G x W FF x Y x M t W mW Tm xyx                 (19) ).()(: 02012 02 2 2 02 2 22 2 2 2 2 2 02 2 0202 WWK x W G x W FF x Y x M t W mW Tm xyx                 (20) The boundary and initial conditions of the double-microbeam system are assumed to be simply supported and considered as 2 2 0 0 0 0 2 2 (0, ) ( , ) 0, (0, ) ( , ) 0, ( 1, 2),i i i i W W W t W L t t L t i x x          (21) 0 0 ( , 0) ( ), ( , 0) ( ), ( 1, 2).i i i i W W x f x x g x i t      (22) 3. ANALYTICAL SOLUTION PROCEDURE For the sake of simplicity, we assume that the two parallel beams of the elastically connected double-beam system have the same bending stiffness EI1=EI2=EI and cross- sectional area A1=A2=A. Both microbeams have the same length L and same material characteristics ρ1=ρ2=ρ. The equations of motion can be expressed in the terms of displacements W01 and W02. By substituting Eqs. (14) into Eqs. (19) and (20) the main equations of EBDMBS in terms of the displacements are given below 4 2 2 2 2 01 01 01 01 0 1 01 024 2 2 2 ( ) ( ) ( ) 0, xx xz m T W W W W D A l m F F G K W W x t x x                 (23) Thermal Effect on the Free Vibration and Buckling of a Double-Microbeam System 51 4 2 2 2 2 02 02 02 02 0 2 01 024 2 2 2 ( ) ( ) ( ) 0. xx xz m T W W W W D A l m F F G K W W x t x x                 (24) In order to simplify the solving of Eqs. (23) and (24), we will introduce the following dimensionless parameters: ).2,1(,,,, ,,,,,, 2 0 0 0 0   i A A a A D d A F F A F F A G k A KL k m A L t L l l L x L W w xx xz xz xx xx xx xx mi mi xx T T xx p xx w xxi i  (25) Assuming time harmonic motion and using separation of variables, the solutions of Eqs. (23) and (24) with the main boundary conditions (21) can be written in the form 0 1 ( , ) ( ) ( ), ( ) sin( ), , ( 1, 2), i n in n n n n w X S X k k n i             (26) where Sin(τ) is the unknown time function, and Xn(ξ) is the known mode shape function for a simply supported single microbeam. Introducing the general solutions (26) into Eqs. (23) and (24) we obtain the following equations 2 4 2 2 1 0 1 1 2 [( ) ( ) ] 0 , n xx xz n m T n p n w n w n S d a l k F F k k k k S k S        (27) 2 4 2 2 2 0 2 2 1 [( ) ( ) ] 0 , n xx xz n m T n p n w n w n S d a l k F F k k k k S k S        (28) The solutions of Eqs. (27) and (28) are assumed in the following forms ,1,, 21  jeDSeCS nn j nn j nn  (29) where ωn marks the natural frequency of the double-microbeam system, and Cn and Dn present the amplitude coefficients of the two microbeams, respectively. By substituting Eqs. (29) into Eqs. (27) and (28) the determinant can be written from which the nontrivial solutions for constants Cn and Dn can be obtained only when this determinant of the coefficients vanishes. This gives the following frequency equations 2 4 2 2 2 0 1 2 2 4 2 2 2 4 2 2 0 1 0 2 2 [2( ) ( ) ( ) 2 2 ] [( ) ( ) ][( ) ( ) ] 0. n xx xz n m T n m T n p n w n xx xz n m T n p n w xx xz n m T n p n w w d a l k F F k F F k k k k d a l k F F k k k k d a l k F F k k k k k                         (30) Finally, when the bi-axial compression forces due to the mechanical loading 1 2m m F F = 0 are ignored, the natural frequency of the system is written by the formula 0 2 2 4 2 2 , 0 ( ) ( ) , nI II xx xz n T n p n w w d a l k F k k k k k      (31) where 0 nI  is the lower natural frequency and 0 nII  is the higher natural frequency of the EBDMBS. When the bi-axial load applied on the double-beam system reaches a certain critical 52 M. STAMENKOVIŠ ATANASOV, D. KARLIŢIŠ, P. KOZIŠ, G. JANEVSKI value, the double-beam system becomes unstable which means that the system begins to buckle. Introducing ωn=0 into Eq. (30), and substituting mechanical load ratio ,12 mm FF where FF m  1 and , 2 FF m  we obtain the equation for the critical buckling load as follows 2 4 2 2 2 2 4 0 1 14 1 {[( ) ](1 ) 4 } , 2 cr xx xz n T n p n w n n n F d a l k F k k k k k b k c k            (32) 2 4 2 2 2 1 0 [( ) ](1 ) , xx xz n T n p n w n b d a l k F k k k k k      (33) 2 4 2 2 2 2 1 0 [( ) ] . xx xz n T n p n w w c d a l k F k k k k k      (34) Based on equation (32), the critical buckling temperature of the EBDMBS for the bi- axial compression equal to zero 0 21  mm FF is of the form 2 2 0 1 ( ) [( ) ]. cr MCST xx xz n p xx T d a l k k A     (35) The illustrated analytical expressions for the natural frequency equation (31), critical buckling load equation (32) and the critical buckling temperature (35) are common equations for the EBDMBS with thermal influence. 3.1. Out-of-phase modes of vibration and buckling state A detailed analysis for different cases of phase modes of vibration and buckling state is shown in the paper of Murmu and Adhikari [33, 34, 35]. Fig. 2 Out-of-phase vibration of the double-microbeam system For the EBDMBS we can use a change in variables by considering W0i(x,t) as the relative displacement of the microbeam-1 with respect to the microbeam-2 , 02010 WWW  (36) then . 02001 WWW  (37) Subtracting Eq. (23) from Eq. (22) and using Eqs. (36) and (37) we obtain 4 2 2 2 2 0 0 0 0 0 04 2 2 2 ( ) ( ) 2 0, xx xz m T W W W W D A l m F F G KW x t x x                (38) Thermal Effect on the Free Vibration and Buckling of a Double-Microbeam System 53 4 2 2 2 2 02 02 02 02 0 04 2 2 2 ( ) ( ) . xx xz m T W W W W D A l m F F G KW x t x x               (39) In above Eqs. (38) and (39) for sake of simplicity we assume that Fm1=Fm2=Fm. If material length scale parameter l is ignored, the above equations become those of the classical Euler–Bernoulli beam theory. For the present out-of-phase analysis of the EBDMBS, we see simplicity in using Eq. (38). The general solution of Eq. (38) is written as ,1,)( 00  iexWW ti (40) where W0(x) is the corresponding deformation shape of the EBDMBS and ω is frequency. For vibration analysis we know that is Fm=0. By introducing Eq. (40) in Eq. (38) we get 4 2 2 20 0 0 04 2 ( ) ( ) ( ) ( ) ( 2 ) ( ) 0, xx xz T W x W x D A l F G m K W x x x            (41) or ,0)( )()( 032 0 2 24 0 4 1       xWA x xW A x xW A (42) where the coefficients are 2 2 1 2 0 3 ( ) , ( ) , ( 2 ) . xx xz T D A l A F G A m K A     (43) The general solution of Eq. (42) can be written as ,coshsinhcossin)( 242312110 xCxCxCxCxW  (44) where Ck, (k=1,2,3,4) can be determined from the boundary conditions (21) and 2 2 2 2 1 2 2 1 3 2 2 2 1 3 1 1 1 1 ( 4 ), ( 4 ). 2 2 A A A A A A A A A A        (45) Further, the solving of frequency for the out-of-phase vibration is presented in this section. The configuration of the EBDMBS with out-of-phase vibration mode (W01  W02  0) is shown in Fig. 2. By using the boundary conditions of simply-supported microbeam system from Eq. (21) yields C2=0 and C4=0. From that we can write 1 2 1 2 2 2 2 1 1 2 2 3 sin sinh 0 . ( ) sin ( ) sinh 0 xx xz xx xz L L C D A l L D A l L C                         (46) For the nontrivial solution of Eq. (46) the determinant is zero, it follows 2 2 2 2 1 2 2 1 2 sin [( ) sinh ( ) sinh ] 0. xx xz xx xz L D A l L D A l L        (47) From Eq. (47) the frequency equation is ,0sin 1  L (48) 54 M. STAMENKOVIŠ ATANASOV, D. KARLIŢIŠ, P. KOZIŠ, G. JANEVSKI and implies ...2,1, 1  nnL (49) Using Eq. (45) yields .0 3 2 12 4 11  AAA  (50) Using the dimensionless parameters (25) and Eq. (43), the natural frequency of the EBDMBS for out-of-phase vibration mode we get 2 4 2 0 ( ) ( ) 2 . n xx xz n T p n w d a l k F k k k      (51) Using dimensionless parameters (25) and Eq. (38) we get the expression of buckling load in out-of-phase sequence as 2 4 2 2 ( ) ( ) 2 , xx xz n T p n w n n d a l k F k k k F k      (52) 3.2. In-phase modes of vibration and buckling state The configuration of the EBDMBS with in-phase modes of vibration is shown in Fig. 3. The relative displacements between the two microbeams are absent (W01  W02 = 0). For the mentioned EBDMBS vibration we solve the Eq. (39). Fig. 3 Out-of-phase vibration of the double-microbeam system By applying the same procedure from the previous chapter the natural frequencies of the EBDMBS for in-phase vibration mode can be expressed as 2 4 2 0 ( ) ( ) . n xx xz n T p n d a l k F k k     (53) The microbeams are buckled in the same direction (synchronous), see Fig. 3.Using dimensionless parameters (25) and Eq. (39) we get the expression of buckling load in in- phase sequence as 2 2 ( ) ( ). n xx xz n T p F d a l k F k    (54) It is shown from Eqs. (53) and (54) that the in-phase vibration mode and buckling state of the EBDMBS is independent of the stiffness of the connecting springs while it is dependent on Pasternak's layer and temperature effect and hence the EBDMBS can be treated as a single microbeam. A similar analysis for nanobeam system is presented in the paper of Murmu and Adhikari [33, 35]. Thermal Effect on the Free Vibration and Buckling of a Double-Microbeam System 55 4. NUMERICAL RESULTS AND DISCUSSION In this section, we have illustrated a comparative study of the analytical results written in this paper and the results found in the literature. The microbeams of the system are made of epoxy with the following properties: ν=0.38, ρ=1220kg/m 3 , E=1.44GPa, l=17.6μm, α=54×10 -6 / ◦ C from [24]. The cross-section shape and length are kept the same by letting b/h=2 and L/h=10 respectively. Temperature and material length scale parameter effect on the two different cases of phase vibration modes and buckling state will be presented. A detailed parametric study is carried out by investigating the influence of different parameters on the natural frequency, frequency under the compressive axial loading, critical buckling load and critical temperature of the EBDMBS with thermal effect. 4.1. Temperature and material length scale parameter effect on the phase vibration modes and phase buckling of the EBDMBS The frequency results of the EBDMBS are presented in terms of the frequency parameters for out-of-phase in Eq. (51) and in phase vibration mode in Eq. (53).Variation in frequency parameter Ωn with material length scale parameter l0, for different phase vibration due to temperature change is shown in Fig. 4. For the Winkler and Pasternak parameter we used constant values of kw=10 and kp=0.1, while for the temperature change we used two different values ΔT=50 ◦ C and ΔT=100 ◦ C. It can be noticed from Fig. 4 that with increasing material length scale parameter l0, frequency parameter Ωn also increases for both considered cases of phase vibration. Frequency parameter Ωn decreases as the temperature effect increases for both considered cases of phase vibration. The buckling state results of the EBDMBS are presented in terms of the buckling parameters for out-of-phase in Eq. (52) and in phase buckling in Eq. (54). Material length scale parameter (l0) on buckling parameter Fn for the different phase buckling due to temperature change is shown in Fig. 5. Also, it can be seen that as material length scale parameter l0 increases, buckling parameter Fn also increases for both considered cases of phase buckling. Buckling parameter Fn decreases as the temperature effect increases for both considered cases of phase buckling. Fig. 4 Variation in frequency parameter Ωn with material length scale parameter (l0) for different phase vibrations due to temperature change 56 M. STAMENKOVIŠ ATANASOV, D. KARLIŢIŠ, P. KOZIŠ, G. JANEVSKI Fig. 5 Material length scale parameter (l0) on buckling parameter Fn for the different phase buckling due to temperature change 4.2. Thermal effect on the natural frequency of the EBDMBS It is commonly known that the lowest natural frequency and buckling load of systems of a larger number of coupled nano/micro structures correspond to the natural frequency and buckling load of one beam or plate, see Karliţiš et al. [31]. In this paper, the results for the lowest natural frequency and critical buckling load of the EBDMBS may be compared with those obtained for a microbeam one, presented in Ke et al. [24] and Ma et al. [19]. In order to confirm the present analytical method, Table 1 shows a comparison of thermal effect on the dimensionless natural frequency of the system from Eq. (31) for the first three modes with results Ke et al [24]. Perfect agreement between the present frequencies and those of Ke et al. [24] can be observed from Table 1. It is shown that the inclusion of the thermal effect decreases the frequencies of the microbeam one. It is seen that the effect of Pasternak parameter kp=0.01, for a greater mode, leads to an increase in natural frequencies. Table 1 Thermal effect on the dimensionless natural frequencies for the three modes of the microbeams with h/l=2 Mode ( )T C n=1 n=2 n=3 n=1 n=2 n=3 n=1 n=2 n=3 Ke et al.[24] Present study for 0 p k Present study for 01.0 p k 0 20 40 60 80 100 0.3478 0.3322 0.3159 0.2986 0.2804 0.2608 1.2890 1.2727 1.2562 1.2394 1.2225 1.2053 2.6277 2.6099 2.5920 2.5739 2.5558 2.5374 0.3582 0.3429 0.3271 0.3104 0.2927 0.2739 1.4328 1.4178 1.4027 1.3875 1.3719 1.3563 3.2238 3.2089 3.1939 3.1788 3.1637 3.1485 0.4764 0.4651 0.4535 0.4416 0.4294 0.4168 1.5645 1.5508 1.5370 1.5231 1.5090 1.4948 3.3588 3.3444 3.3300 3.3156 3.3011 3.2866 It can be seen from Fig. 6 that the natural frequency with Poisson’s ratio (i.e. ν=0.38), suggested by the present Euler-Bernoulli beam model is always higher than that by Poisson’s ratio ν=0. The similar results, merely for a Timoshenko beam, are presented by Ma et al. [19]. We can conclude that there is perfect agreement between the present frequencies and those of Ma et al. [19], when we ignore an effect of the elastic medium, i.e. kp=0 and kw=0. Thermal Effect on the Free Vibration and Buckling of a Double-Microbeam System 57 The temperature effect and the Pasternak’s parameter on the natural frequency can also be noticed in Fig. 6. The natural frequency decreases with temperature effect, in this case 100°C, for a given value of h/l. It is noticed that the inclusion of the constant values of Pasternak’s parameter (kp=0.01) increases the natural frequency of the EBDMBS. As a significant result, Fig. 6 shows that the increase in bending rigidity is suggested by the present model. Also important is that the difference between the natural frequency with the Poisson’s ratio and the one without it is important only when the beam thickness is too small. Fig. 6 The natural frequency of the EBDMBS varying with microbeam thickness, temperature effect and Pasternak’s parameter 4.3. The effect of the compression axial load and temperature effect on the EBDMBS To investigate the influence of the compressive axial loading on the natural frequencies of EBDMBS transverse vibration, we can compare the results of natural frequencies under the compressive axial loading and those without axial loading. , 2 4 2 2 222 , cbb IInI     (55) 2 4 2 2 2 2 0 2( ) ( ) ( ) 2 2 , xx xz n cr T n cr T n p n w b d a l k F F k F F k k k k         (56) 2 4 2 2 2 0 2 4 2 2 2 0 [( ) ( ) ] [( ) ( ) ] , xx xz n cr T n p n w xx xz n cr T n p n w w c d a l k F F k k k k d a l k F F k k k k k               and .1 cr m F F  (57) If we define .,        nII nII nI nI  (58) with vibration mode number n=1 the impact of the compressive axial loading on the natural frequencies of transverse vibration of the EBDMBS presented by ratios of ψ1 and 58 M. STAMENKOVIŠ ATANASOV, D. KARLIŢIŠ, P. KOZIŠ, G. JANEVSKI ψ2 are shown in Fig. 7. Fig. 7 shows that the ratios of frequencies ψ1 and ψ2 decrease with increasing axial compressive load . It can be noticed that the effect of the compressive axial loading on lower natural frequency ωnI is practically independent of axial compression ratio ϑ, whereas on higher natural frequency ωnII it is dependent on it. For the Winkler and Pasternak parameter we used constant values of kw=10 and kp=0.1, while for the temperature change ΔT=50 ◦ C. It can be noticed from Fig. 7(b) that as the axial compression ratio ϑ increases, the ratio of frequency ψ2 decreases. Also, it can be seen that the axial compression ratio on the ratio of frequency ψ2 is independent of axial compression ratio ϑ for small axial compressive load , while it is significant for a large axial compressive load. Fig. 8 shows the thermal effect and effect of the Pasternak parameter on critical buckling load Fcr for the EBDMBS as a function of axial load ratio ϑ . With the axial load ratio ϑ increase, the critical buckling load decreases. For a taken value of axial load ratio ϑ, the critical buckling load of the EBDMBS decreases with an increase in temperature change. As can be seen, for the higher value of the Pasternak parameter of kp=0.1, the critical buckling load has a noticeably higher value. Fig. 7 The thermal effect on the relationships between ratios ψ1 and ψ2 and dimensionless parameter  with increasing axial compressive load ratio ϑ Fig. 8 The thermal effect and effect of Pasternak parameter on the critical buckling load Fcr for the EBDMBS as a function of axial load ratio ϑ Thermal Effect on the Free Vibration and Buckling of a Double-Microbeam System 59 The critical scale load ratio of the modified and the local critical buckling loads at a low temperature environs is presented as . lcrClassica crMCST cr       (59) Fig. 9 Length scale parameter Tcr at low temperature environs In order to make a better illustration of the thermo-mechanical response of the EBDMBS, we introduced a scale parameter. Fig. 9 shows the influence of length scale parameter Tcr at the low temperature environs. The influence of nonlocal parameter at low temperature environs is shown in Karliţiš et al. [32]. It can be observed from Fig. 9 that this parameter increases for a length scale parameter increase. 5. CONCLUSIONS The thermal effect on the free vibration and buckling of the Euler-Bernoulli double- microbeam system is examined in this paper based on the modified couple stress theory. The system is composed of two identical, parallel, simply-supported beams which are continuously joined by the Pasternak’s elastic layer. The temperature change effect is assumed to have an influence on both microbeams. The higher-order main equations and boundary conditions are derived using the Hamilton principle. The separation of variables method (known as the Fourier method) is used for the main equations to obtain free vibration frequencies and critical buckling loads of the EBDMBS. The length scale parameter, temperature change effect, critical buckling load, thickness/material parameter, Pasternak’s parameter and Poisson’s effect are discussed in detail. Also, the effect of different mentioned parameters on the natural frequency, frequency under the compressive axial loading, critical buckling load and critical temperature of the EBDMBS with thermal effect are presented. Effect of the material length scale parameter and thermal effect on the two different cases of phase modes of vibration and buckling state are discussed. Based on the presented analysis we conclude that the in-phase vibration mode and buckling state of the EBDMBS is independent of the stiffness of the connecting springs while it is 60 M. STAMENKOVIŠ ATANASOV, D. KARLIŢIŠ, P. KOZIŠ, G. JANEVSKI dependent on Pasternak's layer and temperature effect and hence the EBDMBS can be treated as a single microbeam. In order to confirm the present study, we have shown in tabular form a comparison of thermal effect on the dimensionless natural frequency of the system for three modes with the results found in the literature. It is concluded that the presented results are in perfect agreement with the results observed in Ke et al. [24]. It is shown that the inclusion of the thermal effect decreases the frequencies of the microbeam one. Also, the effect of Pasternak parameter kp for a greater mode leads to an increase in natural frequencies, but including the temperature change, the frequency is decreased and leads to the decreased stiffness of the system. The numerical results obtained for the natural frequency with Poisson’s effect and suggested by the present Euler-Bernoulli beam model are always higher than those without Poisson’s effect. The thermal effect on the natural frequency is very low for the microbeam one of the EBDMBS and with a small ratio of h/l, while it is significant for the microbeam with a large ratio of h/l. The impact of the compressive axial loading on the natural frequencies of the EBDMBS transverse vibration leads to the following observations:  The temperature change effect has an impact on both microbeams.  The lower and higher natural frequency under the compressive axial loading decrease with increasing axial compressive load and also decrease with a temperature change increase. The reason for that is that the thermal effect leads to the reduction in stiffness and such a behavior leads to the softening of the materials of the EBDMBS.  The effect of the compressive axial loading on the lower natural frequency is almost independent of the axial compression ratio, whereas on the higher natural frequency it depends on it.  For a higher value of the Pasternak parameter, the critical buckling load has a higher value which decreases with increasing temperature change.  The critical buckling temperature for the presented systems is always lower than for the classical theories.  The critical scale load ratio of the modified and the local critical buckling loads at the low temperature environs increases with the increasing length scale parameter. All these observations can be useful for modern electromechanical systems. Physical views of this paper may be useful for the design and vibration analysis of microresonators and microsensors applications. We have shown that using the presented system with the temperature change leads to considerable changes in stiffness, i.e. the thermal effect leads to the reduction in stiffness and such a behavior leads to the softening of the materials of the EBDMBS. Acknowledgements: This research is supported by the research grant of the Serbian Ministry of Science and Environmental Protection under the numbers OI 174001 and OI 174011. Thermal Effect on the Free Vibration and Buckling of a Double-Microbeam System 61 REFERENCES 1. Lam, D. C. C., Yang, F., Chong, A. C. M., Wang, J., Tong, P., 2003, Experiments and theory in strain gradient elasticity, Journal of the Mechanics and Physics of Solids, 51(8), pp. 1477-1508. 2. Gallacher, B. J., Burdess, J. S., Harish, K. M., 2006, A control scheme for a MEMS electrostatic resonant gyroscope excited using combined parametric excitation and harmonic forcing, Journal of Micromechanics and Microengineering, 16(2), 320. 3. Kacem, N., Baguet, S., Hentz, S., Dufour, R., 2011, Computational and quasi-analytical models for non- linear vibrations of resonant MEMS and NEMS sensors, International Journal of Non-Linear Mechanics, 46(3), pp. 532-542. 4. Harish, K. M., Gallacher, B. J., Burdess, J. S., Neasham, J. A., 2009, Experimental investigation of parametric and externally forced motion in resonant MEMS sensors, Journal of Micromechanics and Microengineering, 19(1), 015021. 5. Magrab, E. B., 2012, Vibrations of Elastic Systems: With Applications to MEMS and NEMS, Vol. 184, Springer. 6. Ilic, B., Krylov, S., Bellan, L. M., Craighead, H. G., 2007, Dynamic characterization of nanoelectromechanical oscillators by atomic force microscopy, Journal of applied physics, 101(4), 044308. 7. Hasanyan, DJ., Batra RC., Harutyunyan S., 2008, Pull-in instabilities in functionally graded microthermoelectromechanical systems, J Thermal Stress, 31, pp1006–1021. 8. Rahaeifard, M., Kahrobaiyan, M. H., Ahmadian, M. T., 2009, Sensitivity analysis of atomic force microscope cantilever made of functionally graded materials, 3rd international conference on micro-and nanosystems (MNS3), San Diego (CA, USA), In: DETC 2009-86254. 9. Mindlin, R. D., Eshel, N. N., 1968, On first strain-gradient theories in linear elasticity, International Journal of Solids and Structures, 4(1), pp. 109-124. 10. Wang, B., Zhao, J., Zhou, S., 2010, A micro scale Timoshenko beam model based on strain gradient elasticity theory, European Journal of Mechanics-A/Solids, 29(4), pp. 591-599. 11. Lazopoulos, K. A., Lazopoulos, A. K., 2010, Bending and buckling of thin strain gradient elastic beams, European Journal of Mechanics-A/Solids, 29(5), pp. 837-843. 12. Peddieson, J., Buchanan, G. R., McNitt, R. P., 2003, Application of nonlocal continuum models to nanotechnology, International Journal of Engineering Science, 41(3), pp. 305-312. 13. Eringen, A. C., 1983, On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves, Journal of Applied Physics, 54(9), pp. 4703-4710. 14. Reddy, J. N., 2007, Nonlocal theories for buckling bending and vibration of nanobeams, International Journal of Engineering Science, 45, pp. 288–307. 15. Reddy, J. N., 2002, Energy principles and variational methods in applied mechanics, 2nd ed. New York: John Wiley and Sons. 16. Reddy, J. N., 2008, An introduction to continuum mechanics with applications, New York, Cambridge University Press. 17. Reddy, J. N., 2011, Microstructure-dependent couple stress theories of functionally graded beams, Journal of the Mechanics and Physics of Solids, 59(11), pp. 2382-2399. 18. Park, S. K., Gao, X. L., 2006, Bernoulli–Euler beam model based on a modified couple stress theory, Journal of Micromechanics and Microengineering, 16(11), 2355. 19. Ma, H. M., Gao, X. L., Reddy, J. N., 2008, A microstructure-dependent Timoshenko beam model based on a modified couple stress theory, Journal of the Mechanics and Physics of Solids, 56(12), pp. 3379-3391. 20. Nateghi, A., Salamat-talab, M., Rezapour, J., Daneshian, B., 2012, Size dependent buckling analysis of functionally graded micro beams based on modified couple stress theory, Applied Mathematical Modelling, 36(10), pp. 4971-4987. 21. Bekir, A., Civalek, Ö., 2011, Strain gradient elasticity and modified couple stress models for buckling analysis of axially loaded micro-scaled beams, International Journal of Engineering Science, 49(11), pp. 1268-1280. 22. Şimşek, M., Reddy, J. N., 2013, A unified higher order beam theory for buckling of a functionally graded microbeam embedded in elastic medium using modified couple stress theory, Composite Structures, 101, pp. 47-58. 23. Hendou, R. H., Mohammadi, A.K., 2014, Transient analysis of nonlinear Euler–Bernoulli micro-beam with thermoelastic damping, via nonlinear normal modes, Journal of Sound and Vibration, 333(23), pp. 6224-6236. 62 M. STAMENKOVIŠ ATANASOV, D. KARLIŢIŠ, P. KOZIŠ, G. JANEVSKI 24. Ke, L. L., Wang, Y. S., Wang, Z. D., 2011, Thermal effect on free vibration and buckling of size-dependent microbeams, Physica E: Low-dimensional Systems and Nanostructures, 43(7), pp. 1387-1393. 25. Dutta S. C., Roy R., 2002, A critical review on idealization and modeling for interaction among soil – foundation–structure system, Computers and Structures, 80(1), pp. 1579–159. 26. Mindlin, R. D., 1963, Influence of couple-stresses on stress concentrations. Experimental mechanics, 3(1), pp. 1-7. 27. Mindlin, R. D., Tiersten, H. F., 1962, Effects of couple-stresses in linear elasticity, Archive for Rational Mechanics and Analysis, 11(1), pp. 415-448. 28. Toupin, R. A., 1962, Elastic materials with couple-stresses, Archive for Rational Mechanics and Analysis, 11(1), pp. 385-414. 29. Koiter, W. T., 1964, Couple-stresses in the theory of elasticity: I and II, Proceedings of the Koninklijke Nederlandse Akademie van Wetenschappen, B67, pp. 17–44. 30. Yang, F. A. C. M., Chong, A. C. M., Lam, D. C. C., Tong, P., 2002, Couple stress based strain gradient theory for elasticity. International Journal of Solids and Structures, 39(10), pp. 2731-2743. 31. Karliţiš, D., Koziš, P., Pavloviš, R., 2014, Free transverse vibration of nonlocal viscoelastic orthotropic multi- nanoplate system (MNPS) embedded in a viscoelastic medium, Composite Structures, 115, pp. 89-99. 32. Karliţiš, D., Cajiš, M., Koziš, P., Pavloviš, I., 2015, Temperature effects on the vibration and stability behaviour of multi-layered graphene sheets embedded in an elastic medium, Composite Structures, 131, pp. 672-681. 33. Murmu, T., Adhikari, S., 2010, Nonlocal transverse vibration of double-nanobeam-systems, Journal of Applied Physics, 108(8), p. 083514. 34. Murmu, T., Adhikari, S., 2010, Nonlocal effects in the longitudinal vibration of double-nanorod systems, Physica E: Low-dimensional Systems and Nanostructures, 43(1), pp. 415-422. 35. Murmu, T., Adhikari, S., 2011, Axial instability of double-nanobeam-systems, Physics Letters A, 375(3), pp. 601-608.