The effect of temperature-dependent physical properties and fractional thermoelasticity on nonlocal nanobeams

doi:10.15406/oajmtp.2018.01.00009

Open Access Journal of

eISSN: 2641-9335

Mathematical and Theoretical Physics

Forum Article Volume 1 Issue 2

The effect of temperature-dependent physical properties and fractional thermoelasticity on nonlocal nanobeams

Ahmed E Abouelregal

Verify Captcha

Regret for the inconvenience: we are taking measures to prevent fraudulent form submissions by extractors and page crawlers. Please type the correct Captcha word to see email ID.

Department of Mathematics, Mansoura University, Egypt

Correspondence: Ahmed E Abouelregal, Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516 Egypt

Received: November 10, 2017 | Published: April 4, 2018

Citation: Abouelregal AE. The effect of temperature – dependent physical properties and fractional thermoelasticity on nonlocal nanobeams. Open Acc J Math Theor Phy. 2018;1(2):46-55 DOI: 10.15406/oajmtp.2018.01.00009

Download PDF

Abstract

There are many papers are presented in the field of nanoengineering and nanotechnology. In this work, a new model of heat equation with fractional derivatives and relaxation time is introduced to investigate the thermo elastic vibrations of nonlocal nanobeams. In addition, the thermal conductivity and the modulus of elasticity are taken as a linear function of environmental temperature and the body is due to a ramp−type heating. The basic equations of the considered problem are formulated and solved using Laplace transform along with its numerical inversion. The distributions of the temperature, deflection as well as other variable fields are numerically obtained and illustrated graphically. The effects of the nanoscale, fractional order derivative, ramping-time and relaxation time on the considered variables are concerned and discussed in detail. The introduced model can be reduced to the corresponding classical coupled and Lord and Shulman theories for different values of phase lags and fractional order parameters.

Keywords: fractional heat equation, nanobeams, temperature−dependent properties, ramp−type heating, fractional order.

Introduction

In recent years, many interesting models have been established by using fractional calculus to study the physical manners, particularly in the area of diffusion, viscoelasticity, heat conduction, mechanics of solids and biological systems. The use of fractional system derivatives and integrals leads to the formulation of some physical problems that are more economical and useful than the classical approach.

In the second half of the nineteenth century, Caputo,¹ Caputo and Mainardi² found an agreement between experimental and theoretical results when using fractional derivatives to describe viscoelastic materials. Povstenko³ proposed a fractional time derivative heat conduction equation for a quasi−static uncoupled theory of thermo elasticity. He also investigated in⁴ a generalizations of the heat conduction and Fourier law by using space and time fractional derivatives. Jiang & Xu⁵ obtained a fractional heat equation with a time fractional derivative in the general orthogonal curvilinear and other coordinate system. Samko et al. in their book⁶ provided an excellent historical review of this area under discussion. Formulae establishing relations between the two types of fractional derivatives of Riemann−Leuville and Caputo are discussed in detail in.^7,8 It should be confirmed that if attention is given, the results obtained using the Caputo construction can be modified to the Riemann−Liouville version.

The theory of couple thermo elasticity is extended by Lord & Shulman,⁹ Green and Lindsay¹⁰ by including the thermal relaxation time in constitutive relations. Green & Naghdi¹¹ proposed a new generalized thermo elasticity theory by including the thermal−displacement gradient among the independent constitutive variables. An important feature of this theory, which is not present in other thermo elasticity theories, is that it does not accommodate dissipation of thermal energy.

Recently, a completely new fractional order thermo elasticity model involves one thermal relaxation parameter was introduced by Sherief et al.¹² Also, Ezzat¹³ established a model of fractional heat conduction equation by using the new Taylor series expansion of time−fractional order developed by Jumarie.¹⁴ In addition, authors in¹⁵ presented a two generalized models of heat conduction with fractional order for an isotropic nonhomogeneous thermo elastic solid.

Modern structural elements are often exposed to temperature changes of such magnitude that their physical properties are not fixed even in an approximate sense. The thermal and mechanical properties of materials vary with temperature, so that the temperature dependence of material properties must be taken into consideration in the thermal stress analysis of these elements.^16,17 At high temperatures, material properties such as the thermal conductivity, Poisson’s ratio, the modulus of elasticity and the thermal expansion coefficients are no longer keep constant.¹⁸ In numerous of papers,^19–24 fractional calculus is often applied in vibration studies of rods, beams, plates as well as in other fields of mechanics.

Micro−scale mechanical resonators have high sensitivity as well as fast response and are widely used as sensors and modulators. Micro and nano−mechanical resonators have attracted considerable attention recently due to their many important technological applications. Accurate analysis of various effects on the characteristics of resonators, such as resonant frequencies and quality factors, is crucial for designing high−performance components. Many authors have studied the vibration and heat transfer process of beams. Some of them are found in the references.²⁵⁻²⁷

The nonlocal theory of elasticity was used to study applications in nano–mechanics including lattice dispersion of elastic waves, wave propagation in composites, dislocation mechanics, fracture mechanics, surface tension fluids, etc. Of all the nanostructures, the mechanical behaviors of nanotubes and nanobeams have been most widely investigated. The models of the nonlocal beams expected increasing attention in the earlier few years. In 1972, Eringen introduced the theory of nonlocal continuum mechanics,²⁸ in an effort to deal with the small–scale structure problems. The theories of nonlocal continuum consider the state of stress at a point as a function of the states of strain of all points in the body while the classical continuum mechanics assumes the state of stress at a certain point uniquely depends on the state of strain at that same point. Solutions from various problems support this theory.²⁸⁻³⁵

The current manuscript is an effort to study a thermo elastic problem of a nanobeam loaded thermally by ramp–type heating. Also, a new fractional heat conduction equation model is introduced and the modulus of elasticity is considered to be a linear function of reference temperature. Laplace transform analytical technique is used to calculate the vibration of displacements and temperature. The variations along the axial direction and the through–the–thickness distributions of all fields are investigated. Some comparisons have been also shown graphically to estimate the effects of the small–scale, fractional order and ramping time parameters on all the studied fields.

Mathematical model for nonlocal thermo elasticity with fractional heat conduction

According to the nonlocal elasticity theory of Eringen, the nonlocal differential constitutive equations for a homogenous thermoelastic materials is²⁸⁻³⁰

$(1 - ξ \nabla^{2}) σ_{i j} = τ_{i j}$ (1)

where

$σ_{i j}$ and

$τ_{i j}$ are the nonlocal and local stress tensors, respectively,

$\nabla^{2}$ is Laplacian operator and

$ξ$ is the nonlocal parameter. One may see that when the internal characteristic length is neglected, i.e., the particles of a medium are considered to be continuously distributed,

$ξ$ is zero, and Eq. (1) reduces to the constitutive equation of classical local thermo elasticity. The classical thermo elasticity is based on the principles of the classical theory of heat conductivity, specifically on the classical Fourier law, which relates the heat flux vector

$q$ to the temperature gradient as follows: $q = - K \nabla θ$ (2)

where $K$ is the thermal conductivity of a solid, $θ = T - T_{0}$ is the excess temperature distribution, in which $T_{0}$ is the environmental temperature. Equation (4) together with the energy equation yields the heat conduction equation or the parabolic heat conduction equation and is diffusive with the notion of infinite speed of propagation of thermal disturbances.

$ρ C_{E} \frac{\partial θ}{\partial t} + γ T_{0} \frac{\partial}{\partial t} (div (u)) = - div (q) + Q$ (3)

where $ρ$ is the density, $C_{E}$ is the specific heat, $u$ is the displacement vector, $γ = \frac{Ε α_{τ}}{1 - 2 ν} = α_{τ} (3 λ + 2 μ)$ , being Lame's constants, $α_{τ}$ being the coefficient of linear thermal expansion, $E$ is Young modulus, $ν$ is Poisson's ratio, and $Q$ is the intensity of heat source. The Riemann‒Liouville fractional integral is introduced as a natural generalization of the convolution type integral:^36,37

$I^{α} f (t) = \int_{0}^{t} \frac{{(t - τ)}^{α - 1}}{Γ (α)} f (τ) d τ, (0 < α < 1) .$ (4)

Where $Γ (α)$ is the Gamma function and $f (t)$ is a Lebesgue integrable continuous function satisfies

$\lim_{α \to 1} \frac{d^{α}}{d t^{α}} f (t) = f' (t)$ (5)

Using the Taylor series expansion of time fractional order as in,¹³ a new non–Fourier model was constructed by¹⁴ in the form

$(1 + \frac{τ_{0}^{α}}{α!} \frac{\partial^{α}}{\partial t^{α}}) q = - K \nabla θ$ (6)

Taking divergence of both sides of Eq. (6), we get

$(1 + \frac{τ_{0}^{α}}{α!} \frac{\partial^{α}}{\partial t^{α}}) (div (q)) = - K \nabla^{2} θ$ (7)

From Eqs. (3) and (7), we can get the fractional ordered generalized heat conduction equation as

$K \nabla^{2} θ = (1 + \frac{τ_{0}^{α}}{α!} \frac{\partial^{α}}{\partial t^{α}}) [ρ C_{E} \frac{\partial θ}{\partial t} + γ T_{0} \frac{\partial}{\partial t} (div (u)) - Q] .$ (8)

The introduced new model of thermo elasticity using the methodology of fractional calculus with wide range

$0 < α \leq 1$ covering two cases of conductivity,

$0 < α < 1$ corresponds to weak conductivity,

$α = 1$ for normal conductivity. For the Lord and Shulman (LS) theory

$α = 1$ and for the calssical coupled theory of thermoelasticity (CTE)

$α = 0$ and

$τ_{0} = 0$ . Constitutive equations:

$τ_{i j} = 2 μ e_{i j} + λ e_{i j} - γ θ δ_{i j}$ (9)

Equation of motion:

$σ_{j i, j} + F_{i} = ρ {\ddot{u}}_{i}$ (10)

Equations (1) and (8) describe the nonlocal thermo elasticity theory with fractional order derivative. It can be seen that the corresponding local thermoelasticity with fractional order is recovered by putting

$ξ = 0$ in equation (1). It is known that the values of physical properties, such as the moduli of elasticity and thermal expansion coefficients, cannot be assumed constant over a wide range of high–temperature applications. Consider a thermoelastic body of material having temperature–dependent properties on the form^16,38,39

$(E, γ, K) = (E_{0}, γ_{0}, K_{0}) F (θ)$ (11)

where

$E_{0}, γ_{0}$ and

$K_{0}$ are considered to be constants;

$F (θ)$ is given in a nondimensional function of temperature. In the case of temperature-independent modulus of elasticity,

$F (θ) = 1$ and

$(E, γ, K) = (E_{0}, γ_{0}, K_{0}) .$

The variation of some or all of the these mechanical and thermal properties with temperature can be approximated by linear, exponential, quadratic laws or any other law appropriately approximating the experimental data. We will consider that16,38,39

$F (θ) = e^{- β θ} ≅ 1 - β θ$ (12)

where

$β$ is called the empirical material constant measured in

$1 / K$ . Equation (12) introduces an empirical model that usually works for a limited range of temperatures.

Since only the inﬁnitesimal temperature deviations from reference temperature are considered and for linearity of the governing partial differential equations of the problem, we have to take into account the condition $| T - T_{0} | / T_{0} ≪ 1$ ,which give us the approximating function of $F (θ)$ to be in the form.^40,41

$F (θ) ≅ F (T_{0}) = 1 - β T_{0}$ (13)

Formulation of the problem and modeling of beam structures

We consider a thermo elastic thin nanobeam initially at temperature $T_{0}$ such that x axis is drawn along the axial direction of the beam and $y$ , $z$ axes correspond to the width and thickness, respectively (Figure 1). The small flexural deflections of the nanobeam with dimensions of length $L$ ,width $b$ and thickness $h$ are considered.

The displacement components are given by

$u = - z \frac{\partial w}{\partial x}, v = 0, w (x, y, z, t) = w (x, t),$ (14)

where w is the lateral deflection.

Figure 1 Schematic diagram for the nanobeam.

For a one–dimensional problem, the differential form of the constitutive Eq. (9) after using Eqs. (1) and (14) can be expressed as:^32,33

$σ_{x} - ξ \frac{\partial^{2} σ_{x}}{\partial x^{2}} = - E_{0} F (T_{0}) [\frac{\partial^{2} w}{\partial x^{2}} + α_{T} θ]$ (15)

where

$σ_{x}$ is the nonlocal axial stress, and

$α_{T} = α_{t} / (1 - 2 ν)$ .For transversely vibration of nanobeams, the equilibrium conditions of Euler–Bernoulli theory can be written as

$\frac{\partial^{2} M}{\partial x^{2}} = ρ A \frac{\partial^{2} w}{\partial t^{2}}$ (16)

where

$A = b h$ is the cross–section area.

With aid of Eq. (15), the flexure moment is given by

$M (x, t) - ξ \frac{\partial^{2} M}{\partial x^{2}} = - I E_{0} F (T_{0}) [\frac{\partial^{2} w}{\partial x^{2}} + α_{T} M_{T}]$ (17)

Where

$I = b h^{3} / 12$ is the inertia moment of the cross–section,

$I E_{0}$ is the flexural rigidity of the beam and

$M_{T}$ is the thermal moment,

$M_{T} = \frac{12}{h^{3}} \int_{- h / 2}^{h / 2} θ (x, z, t) z d z$ (18)

Substituting Eq. (17) into Eq. (16), one can get the motion equation of the nanobeam as

$\frac{\partial^{4} w}{\partial x^{4}} + \frac{ρ}{I E_{0} F (T_{0})} (\frac{\partial^{2} w}{\partial t^{2}} - ξ \frac{\partial^{4} w}{\partial t^{2} \partial x^{2}}) + α_{t} \frac{\partial^{2} M_{T}}{\partial x^{2}} = 0$ (19)

Also, it is exactly seen that the flexure moment of the nonlocal nanobeams is given by

$M (x, t) = ξ A ρ \frac{\partial^{2} w}{\partial t^{2}} - I E_{0} F (T_{0}) [\frac{\partial^{2} w}{\partial x^{2}} + α_{t} M_{T}]$ (20)

Substituting the Euler–Bernoulli assumption and Eq. (14) into Eq. (8), gives the geralized heat conduction equation with fractional order derivative without the heat source

$(Q = 0)$ , as

$\frac{\partial^{2} θ}{\partial x^{2}} + \frac{\partial^{2} θ}{\partial z^{2}} = (1 + \frac{τ_{0}^{α}}{α!} \frac{\partial^{α}}{\partial t^{α}}) [\frac{1}{k} \frac{\partial θ}{\partial t} - \frac{γ_{0} T_{0}}{K_{0}} z \frac{\partial}{\partial t} (\frac{\partial^{2} w}{\partial x^{2}})]$ (21)

Solution of the problem

There is no heat flow across the upper and lower surfaces of the nanobeam (thermally insulated), so that $\frac{\partial θ}{\partial z}$ should be vanish at the upper and lower surfaces of the nanobeam $z = \pm h / 2$ . For a very nanobeam, assuming that the increment temperature varies in a sinusoidal form along the thickness direction. That is

$θ (x, z, t) = Θ (x, t) \sin (\frac{π z}{h})$ (22)

Substituting Eq. (22) into Eq. (19), one can get the motion equation of the nanobeams as

$\frac{\partial^{4} w}{\partial x^{4}} + \frac{ρ A}{E_{0} F (T_{0})} (\frac{\partial^{2} w}{\partial t^{2}} - ξ \frac{\partial^{4} w}{\partial t^{2} \partial x^{2}}) + \frac{24 α_{t}}{h π^{2}} \frac{\partial^{2} Θ}{\partial x^{2}} = 0$ (23)

Also, the flexure moment can be determined from Eqs. (20) and (22) as $M (x, t) = ξ A ρ \frac{\partial^{2} w}{\partial t^{2}} - I E_{0} F (T_{0}) [\frac{\partial^{2} w}{\partial x^{2}} + \frac{24 α_{t}}{h π^{2}} Θ]$ (24)

Now, multiplying Eq. (21) by means of

$12 z / h^{3}$ and integrating it with respect to

$z$ through the beam thickness from

$- h / 2$ to

$h / 2$ , yields

$\frac{\partial^{2} Θ}{\partial x^{2}} - \frac{π^{2}}{h^{2}} \frac{\partial^{2} Θ}{\partial z^{2}} = (δ + \frac{τ_{0}^{α}}{α!} \frac{\partial^{α}}{\partial t^{α}}) [\frac{1}{k} \frac{\partial Θ}{\partial t} - \frac{γ_{0} T_{0} π^{2} h}{24 K_{0}} \frac{\partial}{\partial t} (\frac{\partial^{2} w}{\partial x^{2}})]$ (25)

Now, for simplicity we will use the following non-dimensional variables:

$\begin{array}{l} {x^{'}, z^{'}, u^{'}, w^{'}, L^{'}, h^{'}} = c_{0} η_{0} {x, z, u, w, L, h}, {t^{'}, τ'_{0}, ξ^{'}} = c_{0}^{2} η_{0} {t, τ_{0}, η_{0} ξ}, \\ Θ^{'} = \frac{Θ}{T_{0}}, {σ^{'}}_{x} = \frac{σ_{x}}{E_{0}}, M^{'} = \frac{M}{c_{0} η_{0} I E_{0} F (T_{0})}, c_{0} = \sqrt{\frac{E_{0}}{ρ}} . \end{array}$ (26)

So, the basic equations in nondimensional forms are simplified as (dropping the primes for convenience)

$\begin{array}{l} \frac{\partial^{4} w}{\partial x^{4}} + A_{1} (\frac{\partial^{2} w}{\partial t^{2}} - ξ \frac{\partial^{4} w}{\partial t^{2} \partial x^{2}}) = - A_{2} \frac{\partial^{2} Θ}{\partial x^{2}} \\ (\frac{\partial^{2}}{\partial x^{2}} - A_{4}) Θ = (1 + \frac{τ_{0}^{α}}{α!} \frac{\partial^{α}}{\partial t^{α}}) [\frac{\partial Θ}{\partial t} - A_{5} \frac{\partial}{\partial t} (\frac{\partial^{2} w}{\partial x^{2}})] \end{array}$ (27)

$M (x, t) = A_{1} (ξ \frac{\partial^{2} w}{\partial t^{2}} - \frac{\partial^{2} w}{\partial x^{2}}) - A_{2} Θ$ (28)

where

$A_{1} = \frac{12}{F (T_{0}) h^{2}}, A_{2} = \frac{24 T_{0} α_{t}}{π^{2} h}, A_{3} = ξ A_{1}, A_{4} = \frac{π^{2}}{h^{2}}, A_{5} = \frac{γ_{0} k π^{2} h}{24 K_{0}}$

Initial and boundary conditions

To solve the problem, the initial and boundary conditions must be taken into consideration. The homogeneous initial conditions are taken as

$Θ (x, 0) = \frac{\partial Θ (x, 0)}{\partial t} = 0 = w (x, 0) = \frac{\partial w (x, 0)}{\partial t}$ (29)

We will assume that the two ends of the nanobeam are clamped i.e.

$w (0, t) = w (L, t) = 0 = \frac{\partial^{2} w (0, t)}{\partial x^{2}} = \frac{\partial^{2} w (L, t)}{\partial x^{2}}$ (30)

Also, we consider the nanobeam is loaded thermally by ramp-type heating, which give

$Θ (x, t) = Θ_{0} {\begin{matrix} 0, t \leq 0, \\ \frac{t}{t_{0}}, 0 \leq t \leq t_{0} \\ 1, t > 0 \end{matrix},$ (31)

Where t₀ is a non–negative constant called ramp–type parameter and

$Θ_{0}$ is a constant. In addition, the temperature at the end boundary should satisfy the following relation

$\frac{\partial Θ}{\partial x} = 0 on x = L$ (32)

Solution of the problem in the Laplace transform domain

The closed form solution of the governing and constitutive equations can be possible by adapting the Laplace transformation method. Taking the Laplace transform defined by the relation

$\bar{f} (x, t) = \int_{0}^{\infty} f (x, t) e^{- s t} d t$ (33)

to both sides of Eqs. (27) and (28) and using the homogeneous initial conditions (29), one gets the field equations in the Laplace transform space as

$\begin{array}{l} (\frac{d^{4}}{d x^{4}} - A_{3} s^{2} \frac{d^{2}}{d x^{2}} + A_{1} s^{2}) \bar{w} = - A_{2} \frac{d^{2} \bar{Θ}}{d x^{2}}, \\ (\frac{\partial^{2}}{\partial x^{2}} - B_{1}) \bar{Θ} = - B_{2} \frac{\partial^{2} \bar{w}}{\partial x^{2}}, \end{array}$ (34)

$\bar{M} (x, t) = A_{1} (ξ s^{2} \bar{w} - \frac{d^{2} \bar{w}}{d x^{2}}) - A_{2} \bar{Θ}$ (35)

where

$B_{1} = s (δ + \frac{s^{α} τ_{0}^{α}}{α!}) + A_{4}, B_{2} = s (δ + \frac{s^{α} τ_{0}^{α}}{α!}) A_{5}$ Elimination

$\bar{Θ}$ or<

$\bar{w}$ from Eqs. (34), one obtains:

$(D^{6} - A D^{4} + B D^{2} - C) {\bar{Θ}, \bar{w}} (x) = 0,$ (36)

where the coefficients

$A$ ,

$B$ and

$C$ are given by

$A = A_{3} s^{2} + A_{2} B_{2} + B_{1}, B = A_{1} s^{2} + B_{1} A_{3} s^{2}, C = A_{1} B_{1} s^{2}, D = \frac{d}{d x}$ (37)

Equation (36) can be moderated to

$(D^{2} - m_{1}^{2}) (D^{2} - m_{2}^{2}) (D^{2} - m_{3}^{2}) {\bar{Θ}, \bar{w}} (x) = 0,$ (38)

where

$m_{n}^{2},$

$n = 1, 2, 3, 4$ are roots of

$m^{6} - A m^{4} + B m^{2} - C = 0,$ (39)

The solution of the governing equations (39) in the Laplace transformation domain can be represented as

${\bar{w}, \bar{Θ}} (x) = \sum_{n = 1}^{3} ({1, β_{n}} C_{n} e^{- m_{n} x} + {1, β_{n + 3}} C_{n + 3} e^{m_{n} x}) .$ (40) Where the compatibility between these two equations and Eq. (34), gives

$β_{n} = - \frac{m_{n}^{2} B_{2}}{m_{n}^{2} - B_{1}}$ (41)

where

$C_{n}$ and

$β_{n}$ are parameters depending on

$s$ .

The axial displacement after using Eq. (40) takes the form $\bar{u} (x) = - z \frac{d \bar{w}}{d x} = z \sum_{n = 1}^{3} m_{n} (C_{n} e^{- m_{n} x} - C_{n + 3} e^{m_{n} x}) .$ (42)

Substituting the expressions of $\bar{w}$ and $\bar{Θ}$ from (40) into (36), we get at the solution for the bending moment $\bar{M}$ as follows:

$\bar{M} (x) = \sum_{n = 1}^{3} (A_{3} s^{2} - m_{n}^{2} - A_{2} β_{n}) (C_{n} e^{- m_{n} x} C_{n + 3} e^{m_{n} x}) .$ (43)

In addition, the strain will be

$\bar{e} (x) = \frac{d \bar{u}}{d x} = - z \sum_{n = 1}^{3} m_{n}^{2} (C_{n} e^{- m_{n} x} + C_{n + 1} e^{m_{n} x}) .$ (44)

After using Laplace transform, the boundary conditions (30)‒(32) take the forms

$\begin{array}{l} \bar{w} (0, s) = \bar{w} (L, s) = 0, \\ \frac{\partial^{2} \bar{w} (0, s)}{\partial x^{2}} = \frac{\partial^{2} \bar{w} (L, s)}{\partial x^{2}} = 0, \\ \frac{\partial \bar{Θ} (0, s)}{\partial x} = Θ_{0} (\frac{1 - e^{- s t_{0}}}{s^{2} t_{0}}) = \bar{G} (s), \\ \frac{\partial \bar{Θ} (L, s)}{\partial x} = 0 \end{array}$ (45)

Substituting Eq. (40) into the above boundary conditions, one obtains six linear equations;

$\begin{array}{l} \sum_{n = 1}^{3} (C_{n} + C_{n + 1}) = 0, \\ \sum_{n = 1}^{3} (C_{n} e^{- m_{n} L} + C_{n + 1} e^{m_{n} L}) = 0, \end{array}$ (46)

$\begin{array}{l} \sum_{n = 1}^{3} m_{n}^{2} (C_{n} + C_{n + 1}) = 0, \\ \sum_{n = 1}^{3} m_{n}^{2} (C_{n} e^{- m_{n} L} + C_{n + 1} e^{m_{n} L}) = 0, \end{array}$ (47)

$\begin{array}{l} \sum_{n = 1}^{3} m_{n} (β_{n} C_{n} - β_{n + 1} C_{n + 1}) = - G (s), \\ \sum_{n = 1}^{3} m_{n} (β_{n} C_{n} e^{- m_{n} L} - β_{n + 1} C_{n + 1} e^{m_{n} L}) = 0, \end{array}$ (48)

The solution of the above system of linear equations gives the unknown parameters,. In order to determine the studied fields in the physiacl domain, the Riemann-sum approximation method is used to obtain the numerical results. The details of these methods can be found in Honig & Hirdes.33

Special case

The following special cases can be obtained from the system of Eqs. (1), (8) and (9):

The equations of a coupled theory of nonlocal thermoelasticity with temperature-dependent mechanical properties result from (1), (8), in the limiting caseand (11) by lettingand putting.
The equations of a coupled theory of local thermoelasticity with temperature-dependent mechanical properties result from (1), (8), in the limiting caseand (11) by lettingand putting.
The equations of a coupled theory of nonlocal thermoelasticity with temperature-independent mechanical properties result from (1), (8) and (9), in the limiting caseand (11) by letting,and putting.
The equations of a coupled theory of local thermoelasticity with temperature-independent mechanical properties result from (1), (8) and (9), in the limiting caseby lettingand putting.
The equations of a generalized theory of nonlocal thermoelasticity with temperature-dependent mechanical properties without fractional derivatives are obtained from Eqs. (1), (8) and (9) by letting the fractional parameterand putting.
The equations of a generalized theory of nonlocal thermoelasticity with temperature-independent mechanical properties without fractional derivatives are obtained from Eqs. (1), (8) and (9) by letting the fractional parameterand puttingand.
The equations of the generalized nonlocal thermoelasticity without energy dissipation with temperature-dependent mechanical properties without fractional derivatives result from Eqs. (1), (8) and (9) by letting,and.

Numerical results

In the present work, the thermoelastic coupling effect is analyzed by considering a beam made of silicon (Table 1).

Material properties/ Material	Metal (Aluminum)
Thermal conductivity $K$ $(W m^{(- 1)} K^{(- 1)})$	156
Young' modulus $E$ $(G P a)$	169
Density $ρ ({Kgm}^{- 3})$	2330
Thermal expansion $α_{t}$ $(K^{- 1})$	$2.59 \times 10^{- 6}$
Thermal diffusivity $k (m^{2} s^{- 1})$	$84.18 \times 10^{- 6}$
Poisson's ratio $ν$	0.22
Specific heat $C_{E} (J / kgK)$	713

Table 1 Mechanical and thermoelastic properties parameter of the graded nanobeam

Also, the length−to−thickness ratio of the nanobeam is fixed asand the dimensional parameters carried out in the numerical simulations areand. For all numerical calculations Mathematica programming Language has been used. The numerical calculations of the flexure momentthermodynamical temperature, displacement, lateral vibrationand strainhave been considered for various values of the nonlocal parameterfractional orderthickness of nanobeam and periodic frequency. The results are investigated graphically in Figures 2−17 in the wide range ofat different positions of. Some plots consider the present quantitis throug the length of the beam and others take into account both the length and thickness directions.

Figure 2 The transverse deflection with different fractional order parameter .

Figure 3 The temperature with different fractional order parameter .

Figure 4 The displacement with different fractional order parameter .

Figure 5 The flexure moment with different fractional order parameter .

Figure 6 The transverse deflection with different nonlocal parameter .

!--Figure start here-->

Figure 7 The temperature with different nonlocal parameter .

Figure 8 The displacement with different nonlocal parameter .

Figure 9 The flexure moment with different nonlocal parameter .

Figure 10 The transverse deflection with different ramping time .

Figure 11 The temperature with different ramping time .

Figure 12 The displacement with different ramping time .

Figure 13 The flexure moment with different ramping time .

Figure 14 The transverse deflection in the case of temperature−independent properties.

Figure 15 The temperature in the case of temperature−independent properties.

Figure 16 The displacement in the case of temperature−independent properties.

Figure 17 The flexure moment in the case of temperature−independent properties.

From the results we obtained it is clear that, the field quantities such as the deflection, temperature, displacement, and bending moment distributions depend not only on the timeand space coordinate, but also depend on the fractional parameter, nonlocal parameter ramping time parameterand effect of reference temperature. Numerical calculations and graphs have been divided into four cases. In this work, the obtained results have been compared with the results obtained previously in.^22,23,35,4

The effect of the fractional order parameter

In the first case, we represent the dimensionless lateral vibration, thermodynamic temperature, displacement, and bending moment at different fractional order parameter to stand on the effect of this parameter on all the studied fields (Figures 2−5). The computations were carried out for wide range of , for different values of the parameter with wide rangewhen,, and.The different values of the parameter,describes two types of conductivity (weak conductivity,,normal conductivity,), respectively. From the different figures it is observed that the nature of variations of all the field variables for fractional order parameter is significantly different. The difference is more prominent for higher values of fractional derivative. This observation is consistent with previous results in.²²We also observes the following important facts from Figures 3−6.

The lateral vibrationdistribution has large change against to fractional order parameter. Also, the lateral vibration increases with the increase of the parameter.
From figure 1, lateral vibration atis zero as shown, which agrees with the boundary condition prescribed. This coincided with the mechanical boundary condition that the two ends of the nanobeam are clamped.
We observed that at the lateral vibrationreaches a maximum value at some distance from the surface of the nanobeam and then it decreases with the increase of distance.
The temperature distribution have small differences for different values of.
Figure 3 indicate that thetemperature has a maximum value at boundaryand decreases smoothly and finally tends to zero after some distance (which satisfy the boundary conditions).
Figure 3 shows that the displacementincreases as It can be observed that the fractional order parameterhas a great effects on the displacement distribution.
From Figure 4 it is observed that the bending momentincreases as the distanceincreases and then it decreases with the increase of distance. Figure 5 shows that the bending moment larg dependent on the variation of the fractional order parameter.
According to these results, we have to construct a new classification to all the materials according to its fractional parameter, where this parameter becomes new indicator of its ability to conduct the thermal energy; confirms the results obtained in.²²

The effect of the nonlocal parameter

Figures 6−9 describe the temperature, displacement, lateral vibration and bending moment field distributions for different values of nonlocal parameter. In this case, we notice that when the nonlocal parametervanishingindicates the old situation (local model of elasticity) while other values indicate the nonlocal theories of elasticity and thermo elasticity. The figures 6−9 show that this parameter has significant effect on all the fields. The waves reach the steady state depending on the value of the nonlocal parameter. The conclusion remarks from the Figures can be shortened as follows:

The lateral vibrationdecrease when the value of increases while the values of the bending momentincreases when the value of
It can be seen from Figure 7 that the temperaturesmall dependent on the variation of nonlocal parameter .
From Figure 5, the values of the displacementstart increasing with the nonlocal parameter in the range,thereafter increasing to maximum amplitudes in the range. This shows the difference between the local generalized thermoelasticity and the nonlocal generalized thermoelasticity models.
It can be seen from Figure 9 that the bending moment larg dependent on the variation of nonlocal parameter.
As in,^23,35,40 the field variables have pronounced sensitive to the variation of the nonlocal parameter.

The effect of the ramping time parameter

This case illustrating how the field quantities vary with the different values of the ramping time parameterwith constant,, and. The numerical results are obtained and presented graphically in Figures 10−13. We can see the significant effect of the ramping time parameteron all the studied fields. Also, we can conclude that:

The increasing in the value of the ramping time parameter causes decreasing in the values of thelateral vibration, temperatureand bending momentwhich is very obvious in the peek points of the curves.
Also, the values of the displacementstart decreasing with the ramping time parameter in the range, thereafter increasing to maximum amplitudes in the range.
It can be observed that the ramping time parameterhas great effects on the distribution of field quantities. This is also consistent with what has been obtained in.⁴⁰

The effect of the temperature−independent modulus

This case presents the effect of temperature−dependent properties (reference temperature) on the field variables for fixed values of the ramping time parameterand the nonlocal parameteraccording to the equation. For a comparison of the results, the lateral vibration, the temperature, the displacement, and the bending moment of the nano−beam are shown in Figures 14−17. We found that, the increasing in the value of causes increasing in the values of the lateral vibration, temperature, displacement fields which is very obvious in the peek points of the curves. When the modulus of elasticity is a linear function of reference temperatureand in the case of a temperature−independent modulu. On the other hand, high sensitivity is observed for the effect of the reference temperature on all of the studied variables. This also agrees with what Zenkour has reached in.⁴⁰

Conclusion

The thermal and mechanical properties of materials vary with temperature, so that the temperature dependent on material properties must be taken into consideration in the thermal stress analysis of these elements. In this work, a new model of nonlocal generalized thermoelasticity based upon Euler–Bernoulli theory for nanobeams is constructed. The vibration characteristics of the deflection, temperature, displacement and bending moment of nanobeam subjected to ramp–type heating are investigated. The effects of the nonlocal parameter, fractional order parameterand the ramping time parameteron the field variables are investigated. Numerical technique based on the Laplace transformation has been used. The effects of the nonlocal parameter, reference temperature and the ramping time parameter on all the studied field quantities have been shown and presented graphically.

According to the results shown in all figures, it is find that:

The nonlocal parameterhas significant effects on all the studied fields.
The effect of the fractional parameter on all the studied fields is very significant.
The thermo elastic stresses, displacement and temperature have a strong dependency on the ramping time parameter.
The parameterand the dimension cause significant changes on all the studied fields and the phenomenon of finite speeds of propagation are manifested in all figures.
New classification of the materials must be constructed according to the fractional parameter which describes the ability of the material to conduct the heat.
This study is very important for micro scale problems because in these cases the material parameters are temperature dependent.

$λ, μ$	Lam´e’s constants	$K$	thermal conductivity
$α_{t}$	thermal expansion coefficient	$α$	Fractional order parameter
$γ$	coupling parameter	$q_{i}$	components of the heat flows vector
$T_{0}$	environmental temperature	$δ_{i j}$	Kronecker's delta function
$θ = T - T_{0}$	temperature increment	$u_{i}$	displacement components
$T$	absolute temperature	$F_{i}$	body force components
$C_{E}$	specific heat	$Q$	heat source
$e$	cubical dilatation	$τ_{0}$	relaxation time
$σ_{i j}$	nonlocal stress tensor	$h$	nanobeam thickness
$e_{i j}$	strain tensor	$ρ$	material density
$L$	nanobeam length	$b$	nanobeam width
$A = b h$	cross-section area	$o x y z$	Cartesian coordinate
$τ_{i j}$	local stress tensor	$\nabla^{2}$	Laplacian operator
$ξ$	Nonlocal parameter	$E$	Young’s modulus
$Γ (α)$	Gamma function	$t$	the time
$Γ (α)$	Gamma function
$ν$	Poisson's ratio