Nonlocal thermoelastic vibrations for variable thermal conductivity nanobeams due to harmonically varying heat
Ashraf M. Zenkour1 , Ahmed E. Abouelregal2
1Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
1Department of Mathematics, Faculty of Science, Kafrelsheikh University, Kafr El-Sheikh 33516, Egypt
2Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
2Department of Mathematics, College of Science and Arts, University of Aljouf, El-Qurayat, Saudi Arabia
Journal of Vibroengineering, Vol. 16, Issue 8, 2014, p. 3665-3678.
Received 27 August 2014; received in revised form 23 October 2014; accepted 5 November 2014; published 30 December 2014
This article constructs a new model of nonlocal thermoelasticity beam theory with phase-lags considering the thermal conductivity to be variable. A nanobeam subjected to a harmonically varying heat is considered. The nonlocal theories of coupled thermoelasticity and generalized thermoelasticity with one relaxation time can be extracted as limited and special cases of the present model. The effects of the variable thermal conductivity parameter, the nonlocal parameter, the phase-lags and the angular frequency of thermal vibration on the lateral vibration, the temperature, the displacement, and the bending moment of the nanobeam are investigated.
Keywords: nonlocal thermoelasticity beam theory, phase-lags, harmonically varying heat, nanobeam, variable thermal conductivity.
Nanotechnology has wide applications, which include the fabrication of engineering structures at nanoscale enabling production of materials with revolution properties and devices with enhanced functionality. One of these structures is nanobeams, which have been used widely in systems and devices such as atomic force microscope, nanowires, nano-actuators and nano-probes (Chakraverty ). The thermoelastic solutions may be completely incorrect with the ignoring of small-scale effects in sensitive nano-designing fields. Since mechanical behavior of nanostructures could not be predicted accurately by classical continuum theories. So, nonlocal continuum theories have been proposed to include the small scale effect in vibration and buckling of nanobeams.
In 1972, Eringen  proposed a theory, called nonlocal continuum mechanics, in an effort to deal with the small-scale structure problems. In the classical (local) continuum elasticity, the material particles are assumed to be continuously distributed and the stress tensor at a reference point is uniquely determined by the strain mechanics is based on the constitutive functionals being functional of the past deformation histories of all material points of the body concerned. The small length scale effects are counted by incorporating the internal characteristic length such as the length of the C-C bond into the constitutive relationship. Solutions from various problems support this theory (Eringen [2, 3]; Eringen and Edelen ). For examples, the dispersion curves by the nonlocal model are in excellent agreement with those by the Born-Karman theory of lattice dynamics; the dislocation core and cohesive stress predicted by the nonlocal theory are close to those known in the physics of solids (Eringen ; Eringen and Edelen ). Zhang et al.  and Lu et al.  determined the values of scaling effect parameter experimentally and they used these values in vibration and buckling problems of nanostructures. In addition, the first author and his colleagues [7-10] have recently discussed the behavior of nanobeam in the context of a nonlocal thermoelasticity theory.
The present paper investigates the coupling effect (Lord and Shulman ; Green and Lindsay ; Green and Naghdi ) between temperature and strain rate in nanobeams with variable thermal conductivity (Berman ) based on the nonlocal dual-phase-lag model. The solution for the generalized thermoelastic vibration of a nanobeam subjected to time-dependent temperature is developed. The Laplace transform method is used to determine the lateral vibration, temperature, displacement, and bending moment of the nanobeam. The effects of nonlocal parameter, phase-lags, angular frequency of thermal vibration and the variability thermal conductivity parameter are investigated.
2. Mathematical modeling and basic equations
Let us consider a homogenous isotropic thermoelastic solid in Cartesian coordinate systems initially un-deformed and at uniform temperature The modified of classical thermoelasticity model is given by Tzou theory (Tzou [15-17]) in which the Fourier law is replaced by an approximation of the equation:
where denotes the heat flux vector, denotes the thermal conductivity, is the temperature increment of the resonator, in which is the environmental temperature, denotes the temperature distribution, and represent the delay times. The delay time is said to be the phase-lag of temperature gradient and is the phase-lag of heat flux. The above equation may be approximated as:
where . The heat conduction equation corresponding to the dual-phase-lag model proposed by Tzou [15-17] in this case takes the form:
where denotes the specific heat per unit mass at constant strain, is the volumetric strain, in which is the thermal expansion coefficient, is Young’s modulus and is Poisson’s ratio, is the material density of the medium, and is the heat source.
The general nonlocal constitutive equations for a nanobeam can be written as :
where denotes the nonlocal stress tensor, denotes the classical (Cauchy) or local stress tensor and is the Laplacian. The parameter is the internal characteristic length, is the external characteristic length (e.g., crack length, wave-length), and is a constant appropriate to each material.
Eq. (3) describes the non-local coupled dynamical thermoelasticity theory, the generalized thermoelasticity theories proposed by Lord and Shulman , and dual-phase-lag model (Tzou [15-17]) for different sets of values of phase-lags parameters and as follows:
• The nonlocal coupled thermoelasticity (CTE) theory results from Eq. (3) by taking .
• The nonlocal generalized theory with one relaxation time (LS theory) can be obtained when 0, ( is the relaxation time).
• The nonlocal generalized theory with two dual-phase-lags (DPL) theory is given by setting .
It can be seen that the corresponding local thermoelasticity is recovered by putting in the above equations.
3. Formulation of the problem
Let us consider small flexural deflections of a thin elastic nanobeam with dimensions of length , width and thickness , as shown in Fig. 1. We define the -axis along the axis of the beam and the - and -axes correspond to the width and thickness, respectively.
Fig. 1. Schematic diagram of the nanobeam
In equilibrium, the beam is unstrained, unstressed and at temperature everywhere. The displacements according to Euler-Bernoulli beam theory are given by:
where is the lateral deflection. It is to be noted that, the nonlocal theory assumes that stress at a point depends not only on the strain at that point but also on strains at all other points of the body. So, the one-dimensional constitutive equation gives the uniaxial tensile stress only, according to the differential form of the nonlocal constitutive relation proposed by Eringen [2, 3], as:
where is the nonlocal parameter. The flexure moment of the cross-section is given by:
where is the flexural rigidity of the beam in which denotes the inertia moment of the cross-section and is the thermal moment:
The equation of motion is given by:
where denotes the cross-section area. Substituting Eq. (7) into Eq. (9), the equation of motion becomes:
In case of the nonlocal theory, the bending moment can be written as:
It is exactly seen that when the parameter is set identically to zero, the local Euler-Bernoulli beam theory will be obtained.
4. Temperature-dependent thermal conductivity
Generally, the assumption that the solid body is thermosensitivity (the thermal properties of the material vary with temperature) leads to a nonlinear heat conduction problem. The exact solution of such problem can be found by assuming the material to be (simply nonlinear) meaning that the thermal conductivity and specific heat depend on the temperature, but thermal diffusivity is assumed be constant.
The thermal conductivity is assumed to vary linearly with temperature according to :
where is the thermal conductivity at ambient temperature and is the slope of the thermal conductivity-temperature curve divided by the intercept . Now, we will consider the Kirchhoff transformation :
where is a new function expressing the heat conduction. By substituting Eq. (12) in Eq. (13), we get:
From Eq. (13), it follows that:
After substituting Eqs. (15) and (16) into Eq. (3), the new form of the general heat equation of solids with temperature-dependent thermal conductivity is obtained by the form:
where . Substituting the Euler-Bernoulli assumption given in Eq. (5) into Eq. (17) gives the thermal conduction equation for the beam without the heat source , as:
For a nanobeam, assuming that the temperature increment varies sinusoidally along the thickness direction. That is:
Using Eqs. (8), (12) and (15), Eq. (10) will be in the form:
where denotes the speed of propagation of isothermal elastic waves. Since such that , then for linearity, Eq. (20) will be:
Now, multiplying Eq. (18) by means of and integrating it with respect to through the beam thickness from to , yields:
The preceding governing equations can be put in non-dimensional forms by using the following dimensionless parameters:
So, the governing equations, and the bending and constitutive equations in non-dimensional forms are simplified as (dropping the primes for convenience):
5. Initial and boundary conditions
In order to solve the problem, both the initial and boundary conditions should be considered. The initial conditions of the problem are taken as:
These conditions are supplemented by considering the two ends of the nanobeam are subjected to the following boundary conditions:
Let us also consider the micro-beam is loaded thermally by a harmonically varying incidents into the surface of the nanobeam :
where is the angular frequency of thermal vibration ( for a thermal shock problem) and is constant. Using Eqs. (13) and (19), then one gets:
Assuming that the boundary is thermally insulated, this means that the following relation will be satisfied:
Using Eqs. (14) and (19), then one gets:
6. Solution in the Laplace transform domain
The closed form solution of the governing and constitutive equations may be possible by adapting the Laplace transformation method. Applying the Laplace transform to Eqs. (24)-(26), one gets the field equations in the Laplace transform space as:
where an over bar symbol denotes its Laplace transform, denotes the Laplace transform parameter and:
Elimination or from Eqs. (34) and (35) gives the following differential equation for both or :
where the coefficients , and are given by:
Introducing into Eq. (38), one gets:
where and , and are the roots of the characteristic equation:
These roots are given by:
The solution of the governing equations in the Laplace transformation domain can be represented as:
where and are parameters depending on . The compatibility between Eqs. (35) and (44), gives:
Then, the axial displacement after using Eqs. (5) and (44) takes the form:
In addition, the strain will be:
After using Laplace transform, the boundary conditions take the forms:
Substituting Eqs. (44) into the above boundary conditions, one obtains six linear equations in the matrix form as:
After solving the above equations, we have the values of the constants , 1, 2,..., 6 whose solution complete the solution of the problem in the Laplace transform domain. Hence, we obtain the expressions for the non-dimensional lateral vibration, displacement, the stress and other physical quantities of the medium.
After obtaining , the temperature increment can be obtained by solving Eq. (14) and using Eq. (19) to give:
Then, the bending moment given in Eq. (36) in the Laplace domain with the aid of the Eqs. (44) and (45), become:
7. Inversion of the Laplace transforms
In order to invert the Laplace transform in the above equations, we adopt a numerical inversion method based on a Fourier series expansion. Durbin  derived the approximation formula:
where being the Laplace transform. It is known that and tend to when tends to infinity. To apply efficiently each of the above mentioned procedures, one would have first to find, in each case, the value of after which decreases monotonically to 0. This is virtually impossible for the very complicated . So, the real and imaginary part of are evaluated together through a complex, single precision arithmetic subroutine, but are converted into double precision constants for the summation up to NSUM; the results are then turned back to single precision expressions. Thus, one avoids time and storage consuming systematic double precision computation; NSUM can be determined by the convergence criterion:
It should be noted that a good choice of the free parameters and is not only important for the accuracy of the results, but also for the application of the Korrecktur method and the methods for the acceleration of convergence. We found that to 10 give good results for NSUM ranging from 50 to 5000. The values of all parameters in Eq. (39) are defined as and in this paper (Durbin ).
Fig. 2. Distribution of the field quantities through the axial direction for different theories of thermoelasticity
a) The deflection
b) The temperature
c) The radial displacement
d) The bending moment
8. Numerical results
In terms of the Riemann-sum approximation defined in Eq. (52), numerical Laplace inversion is performed to obtain the non-dimensional lateral vibration, temperature, displacement and bending moment in the nanobeam. In the present work, the thermoelastic coupling effect is analyzed by considering a beam made of silicon. The material parameters are given as:
The aspect ratios of the beam are fixed as and . The figures were prepared by using the non-dimensional variables which are defined in Eq. (14) for a wide range of beam length when , and 0.12 sec.
Numerical calculations are carried out for four cases. In the first case the non-dimensional lateral vibration, temperature, displacement and bending moment with different phase-lags are investigated when the pulse-width and nonlocal parameter remain constants 5, . The graphs in Fig. 2(a)-(d) represent the curves predicted by three different theories of thermoelasticity obtained as special cases of the present dual-phase-lag model. The computations are performed for various values of the parameters and to obtain the coupled theory (CTE) (), the Lord and Shulman (LS) theory and and the generalized theory of thermoelasticity proposed by Tzou (DPL) ( and ).
Fig. 3. Distribution of the field quantities through the axial direction for various thermal conductivity parameters
a) The deflection
b) The temperature
c) The radial displacement
d) The bending moment
Fig. 2(a) depicts the distribution of the lateral vibration which always begins and ends at the zero values (i.e. vanishes) and satisfies the boundary condition at . In the context of three theories, the values of starts increasing to maximum amplitudes in the range thereafter it moves in the direction of wave propagation. The behavior of DPL model may be differ than those of other two theories.
It is observed from Fig. 2(b) that the temperature decreases as the distance increases to move in the direction of wave propagation. The temperature of DPL model may be similar than those of other models.
Fig. 2(c) shows that the values of the displacement start decreasing in the range thereafter increasing to maximum amplitudes in the range . The displacement moves directly in the direction of wave propagation in the range .
In Fig. 2(d), the bending moment is zero when at the boundary of the beam and it attains a stationary maximum value at some positions.
It can be observed that the phase-lags parameters have great effects on the distribution of field quantities. The mechanical distributions indicate that the wave propagates as a wave with finite velocity in medium. The fact that in thermoelasticity theories (DPL, LS, CTE), the waves propagate with finite speeds is evident in all these figures. The behavior of the three theories is generally quite similar.
Fig. 4. Distribution of the field quantities through the axial direction for various angular frequencies of thermal vibration
a) The deflection
b) The temperature
c) The radial displacement
d) The bending moment
In second case, we consider three different values the variability thermal conductivity parameter to discuss the thermosensitivity (the thermal properties of the material vary with temperature). We take –1, –0.5 for variable thermal conductivity and when thermal conductivity is independent of temperature. In this case the angular frequency of thermal vibration and phase-lag of the heat flux and the phase-lag of temperature gradient . From Fig. 3(a)-(d), the parameter has significant effects on all the fields. From Fig. 3(a), 3(c) and 3(d), the lateral vibration , the temperature and bending moment decreases as the decreases. From Fig. 3(b), it can be found that the displacement decreases as the parameter value decreases in the range and increases in the range . The medium close to the beam surface suffers from tensile stress, which becomes larger with the time passing. We have noticed from the figures that, the variability thermal conductivity parameter has a significant effect on all the fields, which add an importance to our consideration about the thermal conductivity to be variable. It is clear that the maximum values of occur near the surface of the beam and its magnitude decreases with the increase of .
Fig. 5. Distribution of the field quantities through the axial direction for various nonlocal parameters
a) The deflection
b) The temperature
c) The radial displacement
d) The bending moment
In the third case, three different values of the angular frequency of thermal vibration were considered. For thermal shock problem, we put and for harmonically heat is set to . In this case the variability thermal conductivity parameter is fixed to –0.5. From Fig. 4 (a)-(d), the values of the lateral vibration, temperature, displacement and bending moment decreases as decreases. These figures illustrate that, the angular frequency parameter has a significant effects on all the fields.
Next, the DPL model is used to investigate the non-dimensional lateral vibration, temperature and displacement for different values of the nonlocal parameter –0.5, . The numerical results are obtained and presented graphically in Fig. 5(a)-(d). It is seen that the effect of the nonlocal parameter on all the studied fields is highly significant.
The last case is to investigate the non-dimensional lateral vibration, temperature and displacement for various values of the nonlocal parameter 10–7, (0, 1, 2) when the pulse width remains constant It can be seen that the corresponding local thermoelasticity model is recovered by putting . The variations of the field variables with various values of the nonlocal parameter are depicted in Fig. 5(a)-(d). Generally, the field variables are very sensitive to the variation of the nonlocal parameter.
In this paper, a new model of nonlocal generalized thermoelasticity with dual-phase-lags for the Euler-Bernoulli nanobeam under variable thermal conductivity is constructed. The vibration characteristics of the deflection, temperature, the displacement and the bending moment of nanobeam subjected to harmincally varying heat heating are investigated. The effects of the nonlocal parameter and the angular frequency of thermal vibration on the field variables are investigated.
Numerical simulation results yield the following conclusions:
• The variability thermal conductivity parameter has significant effects on the speed of the wave propagation of all the studied fields.
• The dependence of the thermal conductivity on the temperature has a significant effect on thermal and mechanical interactions.
• According to the results shown in all figures, we find that the nonlocal parameter has significant effects on all the field quantities. On the other hand the thermoelastic stresses, displacement and temperature have a strong dependency on the angular frequency of thermal vibration parameter.
• The dual-phase-lag model of thermoelasticity predicts a finite speed of wave propagation that made the generalized theorem of thermoelasticity more consistent with the physical properties of the material. Also, it was found that the phase lags of the heat flux and temperature gradient and play a significant role on the behavior of all field variables.
• The Lord and Shulman (LS) theory and the classical thermoelasticity theory (CTE) are obtained as special cases of the present model.
• The paper also concludes the governing equation of motion for a nonlocal nanobeam can be formed by replacing the bending moment term in the classical equation of motion with an effective nonlocal bending moment as presented herewith. It is found that the effect of the nonlocal parameter is very pronounced. Finally, this work is expected to be useful to design and analyze the wave propagation properties of nanoscale devices.
This work was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah. The authors, therefore, acknowledge with thanks the DSR technical and financial support.
- Chakraverty L. B. S. Free vibration of nonhomogeneous Timoshenko nanobeams. Meccanica, Vol. 49, 2014, p. 51-67. [Search CrossRef]
- Eringen A. C. Nonlocal polar elastic continua. International Journal of Engineering Scince, Vol. 10, 1972, p. 1-16. [Search CrossRef]
- Eringen A. C. On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. Journal of Applied Physics, Vol. 54, 1983, p. 4703-4710. [Search CrossRef]
- Eringen A. C., Edelen D. G. B. On nonlocal elasticity. International Journal of Engineering Scince, Vol. 10, 1972, p. 233-248. [Search CrossRef]
- Zhang Y. Q., Liu G. R., Xie X. Y. Free transverse vibrations of double-walled carbon nanotubes using a theory of nonlocal elasticity. Physical Review B, Vol. 71, 2005, p. 195404. [Search CrossRef]
- Lu P., Lee H. P., Lu C., Zhang P. Q. Dynamic properties of ﬂexural beams using a nonlocal elasticity model. Journal of Applied Physics, Vol. 99, 2006, p. 1-7. [Search CrossRef]
- Zenkour A. M., Abouelregal A. E. The effect of two temperatures on a functionally graded nanobeam induced by a sinusoidal pulse heating. Structural Engineering and Mechanics, an International Journal, Vol. 51, 2014, p. 199-214. [Search CrossRef]
- Zenkour A. M., Abouelregal A. E. Effect of harmonically varying heat on FG nanobeams in the context of a nonlocal two-temperature thermoelasticity theory. European Journal of Computational Mechanics, Vol. 23, 2014, p. 1-14. [Search CrossRef]
- Zenkour A. M., Abouelregal A. E.,Alnefaie K. A., Abu-Hamdeh N. H., Aifantis E. A refined nonlocal thermoelasticity theory for the vibration of nanobeams induced by ramp-type heating. Applied Mathematics and Computation, Vol. 248, 2014, p. 169-183. [Search CrossRef]
- Zenkour A. M., Abouelregal A. E. Vibration of FG nanobeams induced by sinusoidal pulse heating via a nonlocal thermoelastic model. Acta Mechanica, Vol. 225, Issue 12, 2014. [Search CrossRef]
- Lord H. W., Shulman Y. A generalized dynamical theory of thermoelasticity. Journal of Mechanics and Physics of Solid, Vol. 15, 1967, p. 299-309. [Search CrossRef]
- Green A. E., Lindsay K. A. Thermoelasticity. Journal of Elasticity, Vol. 2, 1972, p. 1-7. [Search CrossRef]
- Green A. E., Naghdi P. M. Thermoelasticity without energy dissipation. Journal of Elasticity, Vol. 31, 1993, p. 189-209. [Search CrossRef]
- Berman R. The thermal conductivity of dielectric solids at low temperatures. Advances on Physics, Vol. 2, 1953, p. 103-140. [Search CrossRef]
- Tzou D. Y. A unified approach for heat conduction from macro- to micro-scales. Journal of Heat Transfer, Vol. 117, 1995, p. 8-16. [Search CrossRef]
- Tzou D. Y. Macro-to-Microscale Heat Transfer: the Lagging Behavior. Washington, DC, Taylor & Francis, 1996. [Search CrossRef]
- Tzou D. Y. Experimental support for the Lagging behavior in heat propagation. Journal of Thermophysics and Heat Transfer, Vol. 9, 1995, p. 686-693. [Search CrossRef]
- Durbin F. Numerical inversion of Laplace transforms: an efficient improvement to Dubner and Abate's method. The Computer Journal, Vol. 17, 1974, p. 371-376. [Search CrossRef]