Next Article in Journal
Maximizing the Minimal Satisfaction—Characterizations of Two Proportional Values
Next Article in Special Issue
The Effects of Variable Thermal Conductivity in Semiconductor Materials Photogenerated by a Focused Thermal Shock
Previous Article in Journal
Construction of Solitary Two-Wave Solutions for a New Two-Mode Version of the Zakharov-Kuznetsov Equation
Previous Article in Special Issue
Kane’s Method-Based Simulation and Modeling Robots with Elastic Elements, Using Finite Element Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Size-Dependent Thermoelastic Vibrations of Nanobeams Subjected to Harmonic Excitation and Rectified Sine Wave Heating

by
Ahmed E. Abouelregal
1,2 and
Marin Marin
3,*
1
Department of Mathematics, College of Science and Arts, Jouf University, Al-Qurayyat 72388, Saudi Arabia
2
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
3
Department of Mathematics and Computer Science, Transilvania University of Brasov, 500036 Brasov, Romania
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(7), 1128; https://doi.org/10.3390/math8071128
Submission received: 15 June 2020 / Revised: 6 July 2020 / Accepted: 8 July 2020 / Published: 10 July 2020
(This article belongs to the Special Issue Applied Mathematics and Solid Mechanics)

Abstract

:
In this article, a nonlocal thermoelastic model that illustrates the vibrations of nanobeams is introduced. Based on the nonlocal elasticity theory proposed by Eringen and generalized thermoelasticity, the equations that govern the nonlocal nanobeams are derived. The structure of the nanobeam is under a harmonic external force and temperature change in the form of rectified sine wave heating. The nonlocal model includes the nonlocal parameter (length-scale) that can have the effect of the small-scale. Utilizing the technique of Laplace transform, the analytical expressions for the studied fields are reached. The effects of angular frequency and nonlocal parameters, as well as the external excitation on the response of the nanobeam are carefully examined. It is found that length-scale and external force have significant effects on the variation of the distributions of the physical variables. Some of the obtained numerical results are compared with the known literature, in which they are well proven. It is hoped that the obtained results will be valuable in micro/nano electro-mechanical systems, especially in the manufacture and design of actuators and electro-elastic sensors.

1. Introduction

In recent decades, due to the rapid advancement in engineering technology and stringent training requirements, dynamics and stability and their control over mechanical vibration have gradually transformed into a fundamental and indivisible branch of study in applied mechanics and related engineering. Mechanical vibration is considered in many conditions and circumstances as a useful phenomenon, used in many areas and can also better serve people’s lives. Vibration analysis is fundamental for developments, as well as structural and mechanical system design. These data support us to predict the performance of the structure under various external loads and to design a control system, which is used to analyze the development vibrations in a cantilever beam instead of the structure itself [1].
The field of micro-electro-mechanical systems (MEMSs) is fast becoming involved in many resistance and correspondence applications. Modern technologies have been created to manufacture a variety of MEMS gadgets to meet the demand for many precision industries. MEMSs consist of elastic mechanical parts, such as micro-bridges, cantilevers, and microscopic films of various geometrical sizes and engineering forms that often contain loads [2]. It is imperative for MEMS designers to understand the mechanical properties of elastic micro-components, considering the ultimate goal of predicting the amount of deflection from the applied load and the additional method of avoiding cracking, performance development, and increasing the shelf life of MEMS devices [3].
Because of its applicability to a variety of problems, MEMSs have been applied in the fields of engineering applications and mathematical physics, including numerical methods and theoretical studies [4,5,6,7,8,9,10,11]. Advanced apparatuses created on the basis of resonators, modern science, microscale switches, telephones, mirrors, and pumps are cases of this approach [12,13].
Both investigations and atomic reproduction calculations have established an important dimension influence in the properties of mechanical materials when the sizes of these engineering structures are very small. As a result, the scale effect plays a principal role in the dynamic and static behavior of materials and micro/nanostructures and cannot be slighted. It is also known that classical mechanics do not represent such dimensional influences in small-scale and nano-scale structures.
The fundamental distinction between the classical theory of elasticity and the nonlocal Eringen theory [14] depends on the meaning of stress. In the local elasticity theory, stress at a point is a function of strain only at that point, while in the nonlocal elasticity theory, the stress at any point is an element of the strains at all points in the continuum. In the theory of local resilience, stress at a point is a function of motion at that point, while in non-local elasticity, stress at any point is an element of strains at all points in the continuum. Thus, the theory of nonlocal elasticity includes information about long-range powers around atoms, and therefore the length of the internal scale is taken into account. However, the nonlocal theory has been connected to different fields of physics, containing the lattice distribution of the elastic and thermoelastic waves, the dislocation mechanism [15], etc.
The dynamical interaction of solid and strong materials that are subjected to load movement is one of strong building and engineering areas, for instance, sea industrial, structural designing, earthquake design, and tribal science. For example, ground motion and stresses are prompted in submerged soils by rapidly moving vehicular loads or surface influence waves due to explosives. Many researchers have considered the dynamic interaction of beams and rods under moving forces. To consider the thermoelastic behavior of the structure caused by a harmonic external force, it is important to obtain the temperature distribution first [16,17,18,19,20,21,22,23,24].
All the studies mentioned above do not involve varying temperatures. Elastic bodies, as well as microbeams, are often in actual engineering applications under a variable temperature environment, so the variable temperature must be taken into account. The generalized thermoelasticity theories consider the influence of the coupling between the temperature and the rate of the strain. In addition, the resulting associated equations are of the hyperbolic type. Consequently, the discrepancy with respect to the infinite velocity of the spread of thermal waves in the classical coupled theory is disregarded. One of the most famous theories of generalized thermoelasticity that we will use in this work is Lord and Shulman's theory [25] that includes the time of thermal relaxation.
The free transversal vibrations of a complex system of coupled nanobeams attracted much attention in the scientific community. Few researchers have linked the theory of thermoelasticity to mechanical vibrations within nanomaterials. Many of them take the effect of the temperature in the form of a thermal load only, similar to the initial stress, and impose the thermal distribution in a specific way. They obtain the temperature distribution by solving the steady-state heat conduction equation with boundary conditions on the lower and upper surfaces of the nanobeam across the thickness.
According to the best knowledge of the authors, this is only one of a few studies that combine the system coupled with elastic nanobeams and generalized thermoelasticity (non-Fourier conduction equation). In addition to what was mentioned, in this paper, we obtained analytical solutions for different physical distributions using Laplace transforms without the need to impose the form of solutions, whether a harmonic solution or otherwise. In addition, the proposed model and methodology introduced in this work can be applied to study and describe the thermodynamic behavior of axial nanomaterial systems, such as nanoplates or nanoscale rods with thermoelastic properties.
In this investigation, the thermoelastic nonlocal theory is applied to the Euler Bernoulli beam problem subjected to a dispersed harmonic excitation load per unit length. The non-Fourier conduction equation depends on the thermal relaxation times [25] that are applied. The Laplace transform procedure is utilized as part of the deduction. The effects due to the harmonic external load, nonlocal and angular frequency parameters are represented graphically and are investigated. The current model can be used in micro/nano-electro-mechanical applications, such as mass flow sensors, accelerometers, relay switches, frequency filters and resonators. The vibration of nanobeams is a significant topic for the study of nanotechnology, as it relates to the optical and electronic properties of the nanobeams.

2. Theoretical Problem Formulations

The schematic representation of the considered system is illustrated in Figure 1, showing a hinged–hinged nanobeam of length L ( 0 x L ) , width b ( b / 2 y b / 2 ) and thickness h ( h / 2   z h / 2 ) and Young’s modulus E , assuming that the Euler–Bernoulli beam theory is employed for modeling the nanobeams and the cross section of the nanobeam is uniform along the entire length. Hence, the displacements of the beam can be written as
u = z w x ,   v = 0 ,   w = w ( x , t ) ,  
where w ( x , t ) denotes the transverse displacement (lateral deflection) of the nanobeam.
Based on Eringen’s nonlocal elasticity theory [14,15], the constitutive relation for a one-dimensional problem, after using Equation (11), can be written as
σ x ξ 2 σ x x 2 = E ( z 2 w x 2 + α T θ ) ,
where σ x is the nonlocal axial stress, θ = T T 0 is the resonator temperature change, T is the distribution of temperature and T 0 denotes the environmental temperature, α T = α t / ( 1 2 ν ) , ν is Poisson’s ratio and α t is the linear thermal expansion.
In Equation (2), the nonlocal parameter is ξ = ( e 0 a ) 2 , where a is the internal characteristic length and e 0 is a suitable parameter that can be determined by the experiment. When the characteristic parameter a is neglected, then ξ = 0 , and Equation (2) reduces to the classical constitutive relation (local elasticity).
The bending moment M ( x , t ) is given by
M ( x , t ) = h / 2 h / 2 z σ x dz .
Substitution of Equation (2) into Equation (3) results in
M ( x , t ) ξ 2 M x 2 = EI ( 2 w x 2 + α T M T )
where I = bh 3 / 12 and M T is the thermal moment defined by
M T = 12 h 3 h / 2 h / 2 θ ( x , z , t ) zdz
When the beam is due to a distributed load q ( x , t ) , the transverse motion equation is as follows [24]:
2 M x 2 = q ( x , t ) + ρ A 2 w t 2 .
where ρ is the density of the material and A is the cross section of the nanobeam.
Substitution of Equation (4) into Equation (6) results in
M ( x , t ) = ξ ( ρ A 2 w t 2 q ) EI ( 2 w x 2 + α T M T ) .
Substituting Equation (7) into Equation (6), Equation (6) can be rewritten in the form
[ 4 x 4 + ρ A EI 2 t 2 ( 1 ξ 2 x 2 ) ] w 1 EI ( 1 ξ 2 x 2 ) q + α T 2 M T x 2 = 0 .
The non-Fourier heat transfer equation, considering the entropy balance proposed by Lord and Shulman [25], which includes the heat flux in addition to its time derivative, is given as
( K θ , i ) , i + ( 1 + τ 0 t ) ( ρ Q ) = ( 1 + τ 0 t ) t ( ρ C E θ + γ T 0 e ) .
where K is the coefficient of thermal conductivity, Q is the heat source, C E denotes the specific heat at constant strain, τ 0 is the relaxation time characteristic according to the Lord and Shulman theory and e = u x is the volumetric strain. Substituting Equation (1) into the heat equation, Equation (9), when   Q = 0 , is given by
K ( 2 x 2 + 2 z 2 ) θ = ( 1 + τ 0 t ) t ( ρ C E θ γ T 0 z 2 w x 2 )

3. Solution of the Problem

We assume that the increment of temperature change is in terms of a sine function (sinusoidal variation) as
θ ( x , z , t ) = Θ ( x , t ) sin ( π h z )
Using relation (11) in the governing Equations (7), (8) and (10), we obtain
[ 4 x 4 + ρ A EI 2 t 2 ( 1 ξ 2 x 2 ) ] w 1 EI ( 1 ξ 2 x 2 ) q + 24 α T π 2 h 2 Θ x 2 = 0 ,
M ( x , t ) = ξ ( ρ A 2 w t 2 q ) EI ( 2 w x 2 + 24 T 0 α T π 2 h Θ ) ,
( 1 + τ θ t ) ( 2 Θ x 2 π 2 h 2 Θ ) = ( 1 + τ q t ) t ( ρ C E K Θ γ T 0 π 2 h 24 K 2 w x 2 ) .
Introducing the following non-dimensional quantities:
{ x , w , u , z , L , h , b }   =   η c { x , w , u , z , L , h , b } , { t , τ q , τ θ }   =   η c 2 { t , τ q , τ θ } , Θ   =   1 T 0 Θ , ξ   =   η 2 c 2 ξ , M   =   1 η cEI M ,   q   =   A EI q , c 2   =   E ρ , η   =   ρ C E K
Substitution of (15) into (12)–(14) results in (dropping the primes for convenience)
[ 4 x 4 + 12 h 2 2 t 2 ( 1 ξ 2 x 2 ) ] w ( 1 ξ 2 x 2 ) q + 24 T 0 α T π 2 h 2 Θ x 2 = 0 ,
2 Θ x 2 π 2 h 2 Θ = ( 1 + τ 0 t ) t ( Θ γ π 2 h 24 K η 2 w x 2 ) ,
M ( x , t ) = 12 ξ h 2 2 w t 2 ξ q ( x ) 2 w x 2 24 T 0 α T π 2 h Θ
In general, the harmonic excitation q ( x , t ) has the form of a sine or cosine function of a single frequency. In this problem, the dynamic load q ( x , t ) can be considered as
q ( x , t ) = F 0 cos ( Ω t ) ,
where F 0 is the magnitude of forcing excitation and Ω is the frequency of the external excitations ( Ω = 0 for the uniformly distributed load).

4. Boundary Conditions

Let us also consider what the boundary conditions are assumed to be, as in Table 1
In Table 1, the function f ( x , t ) is a varying rectified sine wave function which is described mathematically as
f ( x , t ) = Θ 0 | sin ( ω t ) | ,
where Θ 0 is a constant and ω is the frequency of the rectified sine wave.

5. Solution in the Transformed Space

In this problem, the initial conditions are assumed to be
w ( x , t ) | t = 0 = w ( x , t ) t | t = 0 = 0 , Θ ( x , t ) | t = 0 = Θ ( x , t ) t | t = 0 = 0 .
Employing the technique of Laplace transform to Equations (16)–(18), we obtain
[ d 4 dx 4 A 3 s 2 d 2 dx 2 + A 1 s 2 ] w ¯ = A 2 d 2 Θ ¯ dx 2 + g ¯ ( s ) ,
( d 2 dx 2 B 1 ) Θ ¯ = B 2 d 2 w ¯ dx 2 ,
M ¯ ( x , s ) = ( d 2 dx 2 A 3 s 2 ) w ¯ A 2 Θ ¯ ξ g ¯ ( s ) ,
where
A 1   =   12 h 2 , A 2   =   24 T 0 α T π 2 h , A 3   =   ξ A 1 , A 4   =   π 2 h 2 , A 5   =   γ π 2 h 24 K η , B 1   =   A 4   +   s ( 1 + τ 0 s ) , B 2   =   s ( 1 + τ 0 s ) A 5 , g ¯ ( s )   =   sF 0 Ω 2 + s 2
From Equations (22) and (23), we can obtain
[ d 6 dx 6 A d 4 dx 4 + B d 2 dx 2 C ] w ¯ = 0 ,
where
A = B 1 + A 2 B 2 + A 3 s 2 ,   B = s 2 ( A 1 + A 3 B 1 ) ,   C = A 1 B 1 s 2 .
The general solution to Equation (26) can be presented as:
w ¯ ( x , s ) = j = 1 3 ( L j sinh ( m j x ) + M j cosh ( m j x ) ) ,
where L j and M j are the integral parameters and ± m j , ( j = 1 , 2 , 3 ) satisfy the equation
m 6 Am 4 + Bm 2 C = 0 .
From Equations (26) and (25), we can obtain
Θ ¯ ( x , s ) = 1 A 2 B 1 [ d 4 w ¯ dx 4 ( A 2 B 2 + A 3 s 2 ) d 2 w ¯ dx 2 + A 1 s 2 w ¯ g ¯ ( s ) ] .
The general solutions of Equation (10) with the help of Equation (30) can be simplified as
Θ ¯ ( x , s ) = j = 1 3 H j ( L j sinh ( m j x ) + M j cosh ( m j x ) ) + H 4 ,
where
H j = 1 A 2 B 1 [ m j 4 ( A 2 B 2 + A 3 s 2 ) m j 2 + A 1 s 2 ] ,   j = 1 , 2 , 3 ,   H 4 = g ¯ ( s ) A 2 B 1 .
Introducing the solutions of w ¯ and Θ ¯ into Equation (24) yields
M ¯ ( x , s ) = j = 1 3 ( m j 2 A 3 s 2 + A 2 H j ) ( L j sinh ( m j x ) + M j cosh ( m j x ) ) A 2 H 4 g ¯ ( s )
The axial displacement u ¯ is obtained using Equations (28) and (1) as
u ¯ = z w ¯ x = z j = 1 3 m j ( L j cosh ( m j x ) + M j sinh ( m j x ) ) .
Using the Laplace transform domain, the boundary conditions given in Table 1 may be written as
w ¯ ( x , s ) | x = 0 , L = 0 , d 2 w ¯ ( x , s ) dx 2 | x = 0 , L = 0 ,
Θ ¯ ( x , s ) | x = 0 = Θ 0 ( ω cosh ( π s 2 ω ) ω 2 + s 2 ) = G ¯ ( s ) ,
d Θ ¯ dx | x = L = 0 .
Substituting (28) and (30) into Equations (35)–(37), one obtains
j = 1 3 M j = 0
j = 1 3 ( L j sinh ( m j L ) + M j cosh ( m j L ) ) = 0
j = 1 3 m j 2 M j = 0
j = 1 3 m j 2 ( L j sinh ( m j L ) + M j cosh ( m j L ) ) = 0
j = 1 3 H j M j + H 4 = G ¯ ( s )
j = 1 3 H j ( L j m j cosh ( m j x ) + M j m j sinh ( m j x ) ) = 0
It is difficult to obtain direct Laplace transform inversion for the complicated solutions of the transformed studied fields. Thus, the physical solutions are obtained numerically using the Fourier expansion technique [26]. In this method, any function g ¯ ( x , s ) can be transformed to that in the domain of time g ( x , t ) by the relation:
g ( x , t ) = e ct t ( 1 2 g ¯ ( x , c ) + Re n = 1 N f g ¯ ( x , c + in π t ) ( 1 ) n ) ,
where N f is a finite number and the parameter c satisfies the relation ct 4.7 [27].

6. Numerical Results

The defined Tzou approximation technique [28] is employed to achieve the numerical results of the studied fields. The silicon (Si) material is selected for numerical calculation, and the physical constants at T 0 = 293   K are [23]:
E = 169   GPa ,   ρ = 2330   kg / m 3 ,   C E = 713   J / ( kg   K ) , α T = 2.59 × 10 6 ( 1 / K ) ,   ν = 0.22 ,   K = 156   W / ( mK ) .
In the calculation analysis, we take L / h = 10 and b / h = 0.5 and the parameter ξ ¯ ( ξ ¯ = 10 2 ξ ) is also considered. The effects of several effective parameters on the conduct of the deflection of the nanobeam are explored. The numerical results are plotted in Figure 2, Figure 3 and Figure 4, which explain the variations of the distributions of lateral vibration, displacement, temperature and moment, concerning different parameters.
Figure 2a–d illustrate the variety of the dimensionless studied fields of the system versus the distance x when L = 1 , z = h / 3 and t = 0.1 . The figures clearly show the influence of the nanoscale parameter ξ ¯ on the dynamic deflection and other fields of the nanobeam. When taking the nonlocal parameter as ξ ¯ = 0 , the governing equations of the nonlocal theory are reduced to the original local theory. As shown in Figure 2a, when ξ ¯ is increased from 0 to 3.0 , the dynamic response of the deflection w changes from a softening to a stiffness type behavior. It is also detected that the deflection vibration of the nanobeam corresponds to the midpoint and meets the boundary condition at the two ends x = 0 , L . Figure 2b demonstrates the temperature profile of the thermoelastic nanobeam for different values of the nanoscale parameter ξ ¯ . From the figure, it is found that that temperature θ decreased over time, which means the mechanical energy of the nanobeam is wasted in the form of thermal energy. Additionally, by increasing the values of ξ ¯ , we note a decrease in the temperature profile.
From Figure 2c, it is observed that an increase in the value of parameter ξ ¯ results in an increase in the dynamic displacement u profile. The bending moment M increases as the distance x increases, as is shown in Figure 2d. It can be seen from Figure 2d that the bending moment M increases with an increase in the nonlocal parameter ξ ¯ . In previous research [23,24], they found that all the studied fields of the nanobeam clearly depend on the nanoscale parameter ξ ¯ . Now, our calculation demonstrates that this phenomenon is also valid for various values of the parameter ξ ¯ .
In Figure 3a–d the variation of the studied fields of the beam with respect to the magnitude of the forcing excitation F 0 is shown. In the following case, we take the values of the parameters ξ ¯ , ω and Ω , respectively, as { 1 , 3 , 5 } [26]. One can see from these figures that the effect of the forcing excitation F 0 has pronounced influences on the distribution of all studied fields. Additionally, it can be seen that the increase in the value of the forcing excitation F 0 causes an increase in the thermodynamic temperature values, displacement and the lateral fields, which is evident in the peak points of the profiles. This result is consistent with the results of [18,20]. It is clear from Figure 3d that the bending moment M decrease when the forcing excitation is increased. The value of F 0 = 0 indicates no forcing excitation, while other values indicate the harmonic forcing excitation.
The variation in the dimensionless field quantities with the external excitation frequency Ω are shown in Figure 4a–d. The value of Ω = 0 indicates the static forcing excitation, while other values indicate the harmonic forcing excitation. It is apparent from these figures that, when the excitation frequency coefficient is increased, the amplitude of the studied fields increases. The data would seem to suggest that increasing the amount of angular frequency causes a decrease in both the static equilibrium position and the amplitude of the response. When an external harmonic excitation is applied to the nanobeam, the field quantities are more sensitive to the excitation frequency Ω . The results obtained in our research are found to be in good agreement with [18]. In Harrington et al.’s research [29], they found that the temperature changes with the resonance frequency of the external excitation.
The influence of the angular frequency of the rectified sine wave ω of the varying rectified sine wave heating on the vibrations of the studied fields of the nanobeam under harmonic external excitation is presented in Figure 5a–d. The small scale and other parameters ξ ¯ , F 0 and Ω are assumed to be constant in this case. The numerical results of the studied fields are displayed for various angular frequency values ω = 1 , 2 , 3 . It is observed that the parameters of the angular frequency ω affect the characteristics of the nanobeam significantly. From Figure 5a–d, it is found that the amplitude of temperature, deflection and axial displacement, as well as the bending moment, decrease with increasing time.

7. Conclusions

In the current investigation, the thermoelastic vibration of nanobeams under the effect of a harmonic external force and rectified sine wave heating is discussed. Using the nonlocal elasticity theory and non-Fourier heat conduction model, basic equations are derived. The effect of different parameters, such as the nonlocal parameter ξ ¯ , the magnitude of forcing excitation F 0 , the angular frequency of thermal vibration ω and the frequency of the external excitation Ω on the studied fields of the nanobeam was additionally examined. From the obtained results, we found that:
-
The magnitude of forcing excitation, the nonlocal parameter, the frequency of the external excitation and the frequency of the rectified sine wave heating field have a considerable influence on the response of the system behavior.
-
When an external harmonic excitation is applied to the nanobeam, the field quantities are more sensitive to excitation frequencies.
-
The increase in the value of the forcing excitation causes an increase in the thermodynamic temperature values, displacement and in the lateral fields, which is evident in the peak points of the profiles.
-
The results in this study may find applications and requests in the development and design of resonators due to thermal environmental loading.
-
The vibration of the heat response of nanobeams can vary by changing the size of the external force and thermal loads without having to change any other engineering and physical parameters of the nanobeams.
-
The temperature dependence of the frequency of the resonance of resonator devices can be utilized to design and create precision thermometers that are not at all suitable under high sensitivity to rectified sine wave heating and harmonic loads.

Author Contributions

Writing, review and editing, A.E.A. and M.M. Both authors have read and agreed to the present version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mohammedi, S.; Hassanirad, A. Applied and Theoretical Cantilever Beam Free Vibration Analysis. World Acad. Sci. Eng. Tech. 2012, 61, 1619–1622. [Google Scholar]
  2. Younis, M.I. MEMS Linear and Non-Linear Statics and Dynamics; Springer: New York, NY, USA, 2011. [Google Scholar]
  3. Allameh, S.M. An introduction to mechanical-properties-related issues in MEMS structures. J. Mater. Sci. 2003, 38, 4115–4123. [Google Scholar] [CrossRef]
  4. Abouelregal, A.E. Response of thermoelastic microbeams to a periodic external transverse excitation based on MCS theory. Microsyst. Technol. 2018, 24, 1925–1933. [Google Scholar] [CrossRef]
  5. Anh, N.D.; Hai, N.Q.; Hieu, D.V. The equivalent linearization method with a weighted averaging for analyzing of nonlinear vibrating systems. Lat. Am. J. Solids Struct. 2017, 14, 1723–1740. [Google Scholar] [CrossRef] [Green Version]
  6. Sidhardh, S.; Ray, M.C. Size-Dependent elastic response in functionally graded microbeams considering generalized first strain gradient elasticity. Q. J. Mech. Appl. Math 2019, 72, 273–304. [Google Scholar] [CrossRef]
  7. Abouelregal, A.E.; Zenkour, A.M. Thermoviscoelastic response of an axially loaded beam under laser excitation and resting on Winkler’s foundation. Multidiscip. Model. Mater. Struct. 2019, 15, 1238–1254. [Google Scholar] [CrossRef]
  8. Gholipour, A.; Ghayesh, M.H.; Zhang, Y. A comparison between elastic and viscoelastic asymmetric dynamics of elastically supported AFG beams. Vibration 2020, 3, 3–17. [Google Scholar] [CrossRef] [Green Version]
  9. Abouelregal, A.E.; Zenkour, A.M. Thermoelastic problem of an axially moving microbeam subjected to an external transverse excitation. J. Theor. Appl. Mech. 2015, 53, 167–178. [Google Scholar] [CrossRef] [Green Version]
  10. Farokhi, H.; Ghayesh, M.H.; Gholipour, A. Dynamics of functionally graded micro-cantilevers. Int. J. Eng. Sci. 2017, 115, 117–130. [Google Scholar] [CrossRef]
  11. Ghayesh, M.H. Mechanics of viscoelastic functionally graded microcantilevers. Eur. J. Mech. A Solids 2019, 73, 492–499. [Google Scholar] [CrossRef]
  12. Rezazadeh, G.; Khatami, F.; Tahmasebi, A. Investigation of the torsion and bending effects on static stability of electrostatic torsional micro-mirrors. Microsyst. Technol. 2007, 13, 715–722. [Google Scholar] [CrossRef]
  13. Chen, J.Y.; Hsu, Y.C.; Lee, S.S.; Mukherjee, T.; Fedder, G.K. Modeling and simulation of a condenser microphone. Sens. Actuators A Phys. 2008, 145, 224–230. [Google Scholar] [CrossRef]
  14. Eringen, A.C. Nonlocal Continuum Field Theories; Springer Verlag: New York, NY, USA, 2002. [Google Scholar]
  15. Eringen, A.C. On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. J. Appl. Phys. 1983, 54, 4703–4710. [Google Scholar] [CrossRef]
  16. Marin, M.; Craciun, E.M.; Pop, N. Considerations on mixed initial-boundary value problems for micropolar porous bodies. Dyn. Syst. Appl. 2016, 25, 175–196. [Google Scholar]
  17. Riaz, A.; Ellahi, R.; Bhatti, M.M.; Marin, M. Study of heat and mass transfer in the Eyring—Powell model of fluid propagating peristaltically through a rectangular compliant channel. Heat. Transf. Res. 2019, 50, 1539–1560. [Google Scholar] [CrossRef]
  18. Bhatti, M.M.; Ellahi, R.; Zeeshan, A.; Marin, M.; Ijaz, N. Numerical study of heat transfer and Hall current impact on peristaltic propulsion of particle-fluid suspension with compliant wall properties. Mod. Phys. Lett. B 2019, 33, 1950439. [Google Scholar] [CrossRef]
  19. Warminska, A.; Manoach, E.; Warminski, J. Vibrations of a composite beam under thermal and mechanical loadings. Procedia Eng. 2016, 144, 959–966. [Google Scholar] [CrossRef] [Green Version]
  20. Warminska, A.; Warminski, J.; Manoach, E. Nonlinear dynamics of a reduced multimodal Timoshenko beam subjected to thermal and mechanical loadings. Meccanica 2014, 49, 1775–1793. [Google Scholar] [CrossRef] [Green Version]
  21. Alcheikh, N.; Tella, S.A.; Younis, M.I. Adjustable static and dynamic actuation of clamped-guided beams using electrothermal axial loads. Sens. Actuators A Phys. 2018, 273, 19–29. [Google Scholar] [CrossRef] [Green Version]
  22. Chang, T.; Fu, Y.M.; Dai, H.L. Nonlinear dynamic analysis of fiber metal laminated beams subjected to moving loads in thermal environment. Compos. Struct. 2016, 140, 410–416. [Google Scholar]
  23. Abouelregal, A.E.; Zenkour, A.M. Thermoelastic response of nanobeam resonators subjected to exponential decaying time varying load. J. Theor. Appl. Mech. 2017, 55, 937–948. [Google Scholar] [CrossRef] [Green Version]
  24. Zenkour, A.M.; Abouelregal, A.E. Decaying temperature and dynamic response of a thermoelastic nanobeam to a moving load. Adv. Comput. Des. 2018, 3, 1–16. [Google Scholar]
  25. Lord, H.W.; Shulman, Y. A generalized dynamical theory of thermoelasticity. J. Mech. Phys. Solids 1967, 15, 299–309. [Google Scholar] [CrossRef]
  26. Hanig, G.; Hirdes, U. A method for the numerical inversion of Laplace transform. J. Comput. Appl. Math 1984, 10, 113–132. [Google Scholar] [CrossRef] [Green Version]
  27. Zhang, Y.Q.; Liu, G.R.; Xie, X.Y. Free transverse vibrations of double-walled carbon nanotubes using a theory of nonlocal elasticity. Phys. Rev. B 2005, 71, 195–404. [Google Scholar] [CrossRef]
  28. Tzou, D.Y. Macro-To Microscale Heat Transfer: The Lagging Behavior; John Wiley & Sons: Hoboken, NJ, USA, 2014. [Google Scholar]
  29. Harrington, D.A.; Mohanty, P.; Roukes, M.L. Energy dissipation in suspended micromechanical resonators at low temperatures. Phys. B Condens. Matter 2000, 284, 2145–2146. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram for the nanobeam.
Figure 1. Schematic diagram for the nanobeam.
Mathematics 08 01128 g001
Figure 2. (a): The deflection w for different values of the nonlocal parameter ξ ¯ . (b): The temperature θ for different values of the nonlocal parameter ξ ¯ . (c): The displacement u . for different values of the nonlocal parameter ξ ¯ . (d): The distribution of bending moment M for different values of the nonlocal parameter ξ ¯ .
Figure 2. (a): The deflection w for different values of the nonlocal parameter ξ ¯ . (b): The temperature θ for different values of the nonlocal parameter ξ ¯ . (c): The displacement u . for different values of the nonlocal parameter ξ ¯ . (d): The distribution of bending moment M for different values of the nonlocal parameter ξ ¯ .
Mathematics 08 01128 g002aMathematics 08 01128 g002b
Figure 3. (a) The deflection w for different magnitudes of the forcing excitation F0. (b): The temperature θ for different magnitudes of the forcing excitation F0 (c): The displacement u for different magnitudes of the forcing excitation F0. (d): The bending moment M for different magnitudes of the forcing excitation F 0
Figure 3. (a) The deflection w for different magnitudes of the forcing excitation F0. (b): The temperature θ for different magnitudes of the forcing excitation F0 (c): The displacement u for different magnitudes of the forcing excitation F0. (d): The bending moment M for different magnitudes of the forcing excitation F 0
Mathematics 08 01128 g003
Figure 4. (a) The distribution of deflection w for different frequencies of the external excitation Ω (b): The distribution of temperature θ for different frequencies of the external excitation Ω. (c): The distribution of displacement u for different frequencies of the external excitation Ω . (d): The distribution of bending moment M for different frequencies of the external excitation Ω .
Figure 4. (a) The distribution of deflection w for different frequencies of the external excitation Ω (b): The distribution of temperature θ for different frequencies of the external excitation Ω. (c): The distribution of displacement u for different frequencies of the external excitation Ω . (d): The distribution of bending moment M for different frequencies of the external excitation Ω .
Mathematics 08 01128 g004
Figure 5. (a) The deflection w for different thermal vibration frequencies ω . (b) The temperature θ for different thermal vibration frequencies ω . (c) The displacement u for different thermal vibration frequencies ω . (d) The bending moment M for different thermal vibration frequencies ω .
Figure 5. (a) The deflection w for different thermal vibration frequencies ω . (b) The temperature θ for different thermal vibration frequencies ω . (c) The displacement u for different thermal vibration frequencies ω . (d) The bending moment M for different thermal vibration frequencies ω .
Mathematics 08 01128 g005
Table 1. The boundary conditions of the problem.
Table 1. The boundary conditions of the problem.
xMechanical Boundary ConditionsThermal Boundary Conditions
x = 0 w ( x , t ) = 0 , 2 w ( x , t ) x 2 = 0 θ = θ ( x , z , t ) = sin ( pz ) f ( x , t )
x = L w ( x , t ) = 0 , 2 w ( x , t ) x 2 = 0 K θ x = 0

Share and Cite

MDPI and ACS Style

Abouelregal, A.E.; Marin, M. The Size-Dependent Thermoelastic Vibrations of Nanobeams Subjected to Harmonic Excitation and Rectified Sine Wave Heating. Mathematics 2020, 8, 1128. https://doi.org/10.3390/math8071128

AMA Style

Abouelregal AE, Marin M. The Size-Dependent Thermoelastic Vibrations of Nanobeams Subjected to Harmonic Excitation and Rectified Sine Wave Heating. Mathematics. 2020; 8(7):1128. https://doi.org/10.3390/math8071128

Chicago/Turabian Style

Abouelregal, Ahmed E., and Marin Marin. 2020. "The Size-Dependent Thermoelastic Vibrations of Nanobeams Subjected to Harmonic Excitation and Rectified Sine Wave Heating" Mathematics 8, no. 7: 1128. https://doi.org/10.3390/math8071128

APA Style

Abouelregal, A. E., & Marin, M. (2020). The Size-Dependent Thermoelastic Vibrations of Nanobeams Subjected to Harmonic Excitation and Rectified Sine Wave Heating. Mathematics, 8(7), 1128. https://doi.org/10.3390/math8071128

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop