Forum Article Volume 1 Issue 2
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
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.
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,1 Caputo and Mainardi2 found an agreement between experimental and theoretical results when using fractional derivatives to describe viscoelastic materials. Povstenko3 proposed a fractional time derivative heat conduction equation for a quasi−static uncoupled theory of thermo elasticity. He also investigated in4 a generalizations of the heat conduction and Fourier law by using space and time fractional derivatives. Jiang & Xu5 obtained a fractional heat equation with a time fractional derivative in the general orthogonal curvilinear and other coordinate system. Samko et al. in their book6 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,9 Green and Lindsay10 by including the thermal relaxation time in constitutive relations. Green & Naghdi11 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.12 Also, Ezzat13 established a model of fractional heat conduction equation by using the new Taylor series expansion of time−fractional order developed by Jumarie.14 In addition, authors in15 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.18 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.25−27
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,28 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.28−35
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.
According to the nonlocal elasticity theory of Eringen, the nonlocal differential constitutive equations for a homogenous thermoelastic materials is28−30
(1−ξ∇2)σij=τij(1−ξ∇2)σij=τij (1)
whereσijσij andτijτij are the nonlocal and local stress tensors, respectively,∇2∇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 vectorqq to the temperature gradient as follows: q=−K∇θq=−K∇θ (2)
whereKK is the thermal conductivity of a solid,θ=T−T0θ=T−T0 is the excess temperature distribution, in whichT0T0 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.
ρCE∂θ∂t+γT0∂∂t(div(u))=−div(q)+QρCE∂θ∂t+γT0∂∂t(div(u))=−div(q)+Q (3)
where ρρ is the density,CECE is the specific heat,uu is the displacement vector,γ=Εατ1-2ν=ατ(3λ+2μ)γ=Eατ1−2ν=ατ(3λ+2μ) , being Lame's constants,ατατbeing the coefficient of linear thermal expansion,EE is Young modulus,νν is Poisson's ratio, andQQ 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)=∫t0(t−τ)α−1Γ(α)f(τ)dτ, (0<α<1).Iαf(t)=∫t0(t−τ)α−1Γ(α)f(τ)dτ, (0<α<1). (4)
WhereΓ(α)Γ(α) is the Gamma function andf(t)f(t) is a Lebesgue integrable continuous function satisfies
limα→1dαdtαf(t)=f'(t)limα→1dαdtαf(t)=f'(t) (5)
Using the Taylor series expansion of time fractional order as in,13 a new non–Fourier model was constructed by14 in the form
(1+τα0α!∂α∂tα)q=−K∇θ(1+τα0α!∂α∂tα)q=−K∇θ (6)
Taking divergence of both sides of Eq. (6), we get(1+τα0α!∂α∂tα)(div(q))=−K∇2θ(1+τα0α!∂α∂tα)(div(q))=−K∇2θ (7)
From Eqs. (3) and (7), we can get the fractional ordered generalized heat conduction equation asK∇2θ=(1+τα0α!∂α∂tα)[ρCE∂θ∂t+γT0∂∂t(div(u))−Q].K∇2θ=(1+τα0α!∂α∂tα)[ρCE∂θ∂t+γT0∂∂t(div(u))−Q]. (8)
The introduced new model of thermo elasticity using the methodology of fractional calculus with wide range0<α≤10<α≤1 covering two cases of conductivity,0<α<10<α<1 corresponds to weak conductivity,α=1α=1 for normal conductivity. For the Lord and Shulman (LS) theoryα=1α=1 and for the calssical coupled theory of thermoelasticity (CTE) α=0α=0 andτ0=0τ0=0 . Constitutive equations:τij=2μeij+λeij−γθδijτij=2μeij+λeij−γθδij (9)
Equation of motion:σji,j+Fi=ρ¨uiσji,j+Fi=ρ¨ui (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ξ=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 form16,38,39(E,γ,K)=(E0,γ0,K0)F(θ)(E,γ,K)=(E0,γ0,K0)F(θ) (11)
whereE0,γ0E0,γ0 andK0K0 are considered to be constants;F(θ)F(θ) is given in a nondimensional function of temperature. In the case of temperature-independent modulus of elasticity, F(θ)=1F(θ)=1 and(E,γ,K)=(E0,γ0,K0).(E,γ,K)=(E0,γ0,K0).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−βθF(θ)=e−βθ≅1−βθ (12)
whereββ is called the empirical material constant measured in1/K1/K . Equation (12) introduces an empirical model that usually works for a limited range of temperatures.Since only the infinitesimal 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−T0|/T0≪1|T−T0|/T0≪1 ,which give us the approximating function of F(θ)F(θ) to be in the form.40,41
F(θ)≅F(T0)=1−βT0F(θ)≅F(T0)=1−βT0 (13)
We consider a thermo elastic thin nanobeam initially at temperatureT0T0 such that x axis is drawn along the axial direction of the beam andyy ,zz axes correspond to the width and thickness, respectively (Figure 1). The small flexural deflections of the nanobeam with dimensions of length LL ,widthbb and thicknesshh are considered.
The displacement components are given by
u=−z∂w∂x, v=0, w(x,y,z,t)=w(x,t),u=−z∂w∂x, v=0, w(x,y,z,t)=w(x,t), (14)
where w is the lateral deflection.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−ξ∂2σx∂x2=−E0F(T0)[∂2w∂x2+αTθ]σx−ξ∂2σx∂x2=−E0F(T0)[∂2w∂x2+αTθ] (15)
whereσxσx is the nonlocal axial stress, andαT=αt/(1−2ν)αT=αt/(1−2ν) .For transversely vibration of nanobeams, the equilibrium conditions of Euler–Bernoulli theory can be written as∂2M∂x2=ρA∂2w∂t2∂2M∂x2=ρA∂2w∂t2 (16)
whereA=bhA=bh is the cross–section area.With aid of Eq. (15), the flexure moment is given by
M(x,t)−ξ∂2M∂x2=−IE0F(T0)[∂2w∂x2+αTMT]M(x,t)−ξ∂2M∂x2=−IE0F(T0)[∂2w∂x2+αTMT] (17)
WhereI=bh3/12I=bh3/12 is the inertia moment of the cross–section, IE0IE0 is the flexural rigidity of the beam and MTMT is the thermal moment,MT=12h3∫h/2−h/2θ(x,z,t)zdzMT=12h3∫h/2−h/2θ(x,z,t)zdz (18)
Substituting Eq. (17) into Eq. (16), one can get the motion equation of the nanobeam as∂4w∂x4+ρIE0F(T0)(∂2w∂t2−ξ∂4w∂t2∂x2)+αt∂2MT∂x2=0∂4w∂x4+ρIE0F(T0)(∂2w∂t2−ξ∂4w∂t2∂x2)+αt∂2MT∂x2=0 (19)
Also, it is exactly seen that the flexure moment of the nonlocal nanobeams is given byM(x,t)=ξAρ∂2w∂t2−IE0F(T0)[∂2w∂x2+αtMT]M(x,t)=ξAρ∂2w∂t2−IE0F(T0)[∂2w∂x2+αtMT] (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)(Q=0) , as ∂2θ∂x2+∂2θ∂z2=(1+τα0α!∂α∂tα)[1k∂θ∂t−γ0T0K0z∂∂t(∂2w∂x2)]∂2θ∂x2+∂2θ∂z2=(1+τα0α!∂α∂tα)[1k∂θ∂t−γ0T0K0z∂∂t(∂2w∂x2)] (21)There is no heat flow across the upper and lower surfaces of the nanobeam (thermally insulated), so that∂θ∂z∂θ∂z should be vanish at the upper and lower surfaces of the nanobeamz=±h/2z=±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(πzh)θ(x,z,t)=Θ(x,t)sin(πzh) (22)
Substituting Eq. (22) into Eq. (19), one can get the motion equation of the nanobeams as
∂4w∂x4+ρAE0F(T0)(∂2w∂t2−ξ∂4w∂t2∂x2)+24αthπ2∂2Θ∂x2=0∂4w∂x4+ρAE0F(T0)(∂2w∂t2−ξ∂4w∂t2∂x2)+24αthπ2∂2Θ∂x2=0 (23)
Also, the flexure moment can be determined from Eqs. (20) and (22) as M(x,t)=ξAρ∂2w∂t2−IE0F(T0)[∂2w∂x2+24αthπ2Θ]M(x,t)=ξAρ∂2w∂t2−IE0F(T0)[∂2w∂x2+24αthπ2Θ] (24)
Now, multiplying Eq. (21) by means of 12z/h312z/h3 and integrating it with respect to zz through the beam thickness from −h/2−h/2 to h/2h/2 , yields∂2Θ∂x2−π2h2∂2Θ∂z2=(δ+τα0α!∂α∂tα)[1k∂Θ∂t−γ0T0π2h24K0∂∂t(∂2w∂x2)]∂2Θ∂x2−π2h2∂2Θ∂z2=(δ+τα0α!∂α∂tα)[1k∂Θ∂t−γ0T0π2h24K0∂∂t(∂2w∂x2)] (25)
Now, for simplicity we will use the following non-dimensional variables:{x′,z′,u′,w′,L′,h′}=c0η0{x,z,u,w,L,h}, {t′,τ'0,ξ′}=c20η0{t,τ0,η0ξ}, Θ′=ΘT0, σ′x=σxE0, M′=Mc0η0IE0F(T0), c0=√E0ρ. (26)
So, the basic equations in nondimensional forms are simplified as (dropping the primes for convenience)
∂4w∂x4+A1(∂2w∂t2−ξ∂4w∂t2∂x2)=−A2∂2Θ∂x2(∂2∂x2−A4)Θ=(1+τα0α!∂α∂tα)[∂Θ∂t−A5∂∂t(∂2w∂x2)] (27)
M(x,t)=A1(ξ∂2w∂t2−∂2w∂x2)−A2Θ (28)
whereA1=12F(T0)h2,A2=24T0αtπ2h, A3= ξA1,A4=π2h2,A5=γ0kπ2h24K0
To solve the problem, the initial and boundary conditions must be taken into consideration. The homogeneous initial conditions are taken as
Θ(x,0)=∂Θ(x,0)∂t=0=w(x,0)=∂w(x,0)∂t (29)
We will assume that the two ends of the nanobeam are clamped i.e.w(0,t)=w(L,t)=0=∂2w(0,t)∂x2=∂2w(L,t)∂x2 (30)
Also, we consider the nanobeam is loaded thermally by ramp-type heating, which giveΘ(x,t)=Θ0{0, t≤0,tt0, 0≤t≤t01, t>0, (31)
Where t0 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∂Θ∂x=0 on x=L (32)
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
ˉf(x,t)=∫∞0f(x,t)e−stdt (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(d4dx4−A3s2d2dx2+A1s2)ˉw=−A2d2ˉΘdx2, (∂2∂x2−B1)ˉΘ=−B2∂2ˉw∂x2, (34)
ˉM(x,t)=A1(ξs2ˉw−d2ˉwdx2)−A2ˉΘ (35)
where B1=s(δ+sατα0α!)+A4, B2=s(δ+sατα0α!)A5 Elimination ˉΘ or< ˉw from Eqs. (34), one obtains:(D6−AD4+BD2−C){ˉΘ,ˉw}(x)=0, (36)
where the coefficients A , B and C are given byA=A3s2+A2B2+B1, B=A1s2+B1A3s2, C=A1B1s2,D=ddx (37)
Equation (36) can be moderated to (D2−m21)(D2−m22)(D2−m23){ˉΘ,ˉw}(x)=0, (38) wherem2n, n=1,2,3,4 are roots ofm6−Am4+Bm2−C=0, (39)
The solution of the governing equations (39) in the Laplace transformation domain can be represented as {ˉw,ˉΘ}(x)=∑3n=1({1,βn}Cne−mnx+{1,βn+3}Cn+3emnx). (40) Where the compatibility between these two equations and Eq. (34), givesβn=−m2nB2m2n−B1 (41)
where Cn and βn are parameters depending on s .The axial displacement after using Eq. (40) takes the form ˉu(x)=−zdˉwdx=z∑3n=1mn(Cne−mnx−Cn+3emnx). (42)
Substituting the expressions of ˉw and ˉΘ from (40) into (36), we get at the solution for the bending moment ˉM as follows:
ˉM(x)=∑3n=1(A3s2−m2n−A2βn)(Cne−mnxCn+3emnx). (43)
In addition, the strain will be
ˉe(x)=dˉudx=−z∑3n=1m2n(Cne−mnx+Cn+1emnx). (44)
After using Laplace transform, the boundary conditions (30)‒(32) take the formsˉw(0,s)=ˉw(L,s)=0, ∂2ˉw(0,s)∂x2=∂2ˉw(L,s)∂x2=0,∂ˉΘ(0,s)∂x=Θ0(1−e−st0s2t0)=ˉG(s), ∂ˉΘ(L,s)∂x=0 (45)
Substituting Eq. (40) into the above boundary conditions, one obtains six linear equations;∑3n=1(Cn+Cn+1)=0,∑3n=1(Cne−mnL+Cn+1emnL)=0, (46)
∑3n=1m2n(Cn+Cn+1)=0,∑3n=1m2n(Cne−mnL+Cn+1emnL)=0, (47)
∑3n=1mn(βnCn−βn+1Cn+1)=−G(s),∑3n=1mn(βnCne−mnL−βn+1Cn+1emnL)=0, (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):
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 (Wm(−1)K(−1)) |
156 |
Young' modulus E (GPa) |
169 |
Density ρ(Kgm−3) |
2330 |
Thermal expansion αt (K−1) |
2.59×10−6 |
Thermal diffusivity k (m2s−1) |
84.18×10−6 |
Poisson's ratio ν |
0.22 |
Specific heat CE (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 start here-->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.22We also observes the following important facts from Figures 3−6.
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 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 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.40
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:
λ,μ | Lam´e’s constants |
K | thermal conductivity |
αt | thermal expansion coefficient |
α |
Fractional order parameter |
γ |
coupling parameter |
qi |
components of the heat flows vector |
T0 |
environmental temperature |
δij |
Kronecker's delta function |
θ=T−T0 |
temperature increment |
ui |
displacement components |
T |
absolute temperature |
Fi |
body force components |
CE |
specific heat |
Q |
heat source |
e |
cubical dilatation |
τ0 |
relaxation time |
σij |
nonlocal stress tensor |
h |
nanobeam thickness |
eij |
strain tensor |
ρ |
material density |
L |
nanobeam length |
b |
nanobeam width |
A=bh |
cross-section area |
oxyz |
Cartesian coordinate |
τij |
local stress tensor |
∇2 |
Laplacian operator |
ξ |
Nonlocal parameter |
E |
Young’s modulus |
Γ(α) |
Gamma function |
t |
the time |
Γ(α) |
Gamma function |
|
|
ν |
Poisson's ratio |
|
|
None
The author declares that there is no conflict of interest.
©2018 Abouelregal. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.