Thermal vibration in rotating nanobeams with temperature-dependent due to exposure to laser irradiation

Effective classical representations of heterogeneous systems fail to have an effect on the overall response of components on the spatial scale of heterogeneity. This effect may be critical if the effective continuum subjects' scale differs from the material's microstructure scale and then leads to size-dependent effects and other deviations from conventional theories. This paper is concerned with the thermoelastic behavior of rotating nanoscale beams subjected to thermal loading under mechanical thermal loads based on the non-local strain gradient theory (NSGT). Also, a new mathematical model and governing equations were constructed within the framework of the extended thermoelastic theory with phase delay (DPL) and the Euler-Bernoulli beam theory. In contrast to many problems, it was taken into account that the thermal conductivity and specific heat of the material are variable and linearly dependent on temperature change. A specific operator has been entered to convert the nonlinear heat equation into a linear one. Using the Laplace transform method, the considered problem is solved and the expressions of the studied field variables are obtained. The numerical findings demonstrate that a variety of variables, such as temperature change, Coriolis force due to rotation, angular velocity, material properties, and nonlocal length scale parameters, have a significant influence on the mechanical and thermal waves.

Generalized thermoelasticity theories have been developed to solve the infinite speed heat propagation problems predicted by the classical thermoelastic theory. By substituting a unique law of thermal conductivity for the standard Fourier law, Lord and Shulman [42] proposed one of the modified generalized thermoelasticity theories with one relaxation time. Recently, the Tzou model is one of the most recent models that has gained fame in the field of thermoelasticity [43][44][45]. In order to address the problems of the infinite speed of thermal waves, which are anticipated by conventional thermoelastic theory, several proposals have been offered in addition to the above-mentioned models [46][47][48][49].
Many studies focusing on the dynamics, bending, statics, and vibrational properties of nanobeams have been presented, with few studies focusing on the thermal elastic behavior of small-sized rotating beams with variable physical properties. In this paper, motivated by the need for improvement, the dynamic transverse vibration properties and thermal vibration analyzes of rotating nanobeams are studied based on the proposal of a new methodology based on nonlocal strain gradient theory (NSGT) that captures the small-scale effects. Also, the equations for the governing system were derived based on the Euler-Bernoulli theory and generalized thermoelasticity with two-phase lag (DPL). Contrary to many other concerns, the thermal conductivity and specific heat of the material have been considered variable and linearly dependent on changes in temperature.
The proposed model is able to study rotating nanobeams under thermal loads and includes Coriolis effects caused by rotation. To the authors' knowledge, for the first time, a rotational model of a thermal nanobeam has been reported whose properties depend on the temperature change. This study will also contribute to a full understanding of the behavior of rotating nanobeams, which have variable physical properties. The Laplace transform method was used to solve constitutive equations and partial differential governing equations, and is an important addition to this article. The Laplace transforms are inverted using a numerical approach based on Fourier series expansions. According to the numerical results, the mechanical and thermal wave responses of the rotating nanobeam are greatly influenced by many parameters such as temperature change, angular velocity, nonlocality parameters, and thermal conductivity change. The results of this study may be useful for future research and accurate design of nanomachines, such as nanoparticle bearings and nanogears, etc.
A micro/nanoscale size influence, for example, is represented by a stiffer elastic response to external stresses. This has been found in the plastic deformation of metals and polymers. In terms of the size effect on elasticity, Lam et al. [54] discovered a rise in the bending rigidities of epoxy microbeams. When the beam thicknesses were reduced from 120 to 20 µm, the bending rigidities were roughly 2.4 times higher than expected by conventional theory. Similarly, additional studies have revealed an apparent rise in Young's modulus without reference to higher-order models. In the absence of strain gradients (for instance, in uniaxial tensile testing), Lam et al. [54] demonstrated that the elastic response of epoxy is independent of sample thickness, which is supported by strain gradient theories.

The dual-phase-lag thermoelastic mode (DPL)
The classical Fourier's law and energy equation for a homogeneous, isotropic, thermoelastic are as follows: Combining (1) and (2) yields the classical heat conduction equation as Tzou [43][44][45] devised the dual-phase-lag (DPL) model that ignores instantaneous temperatureenergy interactions. The DPL model, in particular, makes heat transmission in thermoelastic materials accurate and easy. The modified Fourier law, which has been superseded by Tzou [43][44][45], can be expressed as The proposed new generalized dual-phase-lag thermoelastic model can be obtained by combining (13) and (15) as In Eqs (12)- (16), ⃗ is the heat flux vector, is the heat source, is the thermal conductivity, is the specific heat, is the material density, = div ⃗⃗, is the phase lag of heat flux, and is the temperature gradient's phase lag. Equation (16) is simplified to the constitutive relation of the LS model by setting = 0 and omitting the influence of the term 2 . Equation (7) can be reduced to the classical Fourier law when = = 0.

Problem formulation
As shown in Figure 1, a nanobeam with length , width and thickness ℎ , internal crosssectional area , and bending stiffness are the components of the system under discussion. Under the Euler-Bernoulli principle, every plane cross-section is perpendicular to the nanobeam pivot in the beginning and perpendicular to the unbiased surface. The displacements are given by [53] where is the deflection. Equation (17) can be used to simplify Eq (9) as where is the nonlocal axial stress, = /(1 − 2 ) and = ( 0 ) 2 . The bending moment ( , ) is given by: Substituting Eq (18) into Eq (19), we get where = ℎ 3 /12 and The transverse motion equation can be written as [51] We suppose that the nanobeam rotates about an axis parallel to the -axis with an angular velocity Ω centered at a small distance from the first edge of the nanobeam. The centrifugal tensional force ( ) is introduced because of rotation. In this case, the equation of transverse motion (22) can be written as [27,28] The axial force ( ) due to centrifugal stiffening at a distance from the origin ( Figure 1) is given by [27,28,39] ( ) = ∫ Ω 2 ( + ) .

Thermal properties of materials
The conductivity is considered a linear function of the temperature variation as [23,39,40] = ( ) = 0 (1 + 1 ), where 0 denotes the thermal conductivity at temperature 0 and 1 is the temperature-dependent thermal conductivity fluctuation.
By substituting Eq (29) into (28), we obtain To linearize the governing equation (30), we define a variable as follows [39,40] = ∫ After inserting Eq (29) into Eq (31) and integrating, we have By differentiating the relationship (31) times in terms of distances and once in terms of time, we get Then, the heat conduction equation can be expressed as

Sinusoidal solution
To solve the problem, we take the temperature change solution as a sinusoidal function as Presenting Eq (35) in Eqs (26), (27) and (34) leads to The following non-dimensional variables are provided for convenience We can get the following after introducing dimensionless quantities (39) into the Eqs (36)-(38) (primes are omitted) The maximum axial force ( ) due to centrifugal stiffening at the root ( = 0) is given by [38,39] The motion Eq (40) can therefore be described as The bending moment ( , ) in Eq (42) can also be expressed as

Application
We assume that the nanobeam fulfills the following boundary conditions: i) Mechanical boundary conditions [55] ( , )| =0, = 0, ii) Thermal boundary conditions: In Eq (58), is the time length of a laser pulse and 0 the laser intensity. Introducing Eqs (32) and (33) into Eqs (58) and (59), we get In the Laplace transform domain, the boundary conditions (57)-(61) can be written as When the previously mentioned conditions are applied, then we have By solving the aforementioned system equations, the unknown parameters , ( = 1,2. . ,8) may be computed.

Laplace transform inversion
The Fourier series expansion, as well as the Honig and Hirdes technique [56], will be briefly introduced in this paper. In this method, all functions ̅ ( , ) in the Laplace field can be converted to the time field ̅ ( , ) by using the following relation: where the parameter 1 denotes the time interval and Re stands for the real part of the complex function. By summing a specific number of , Equation (80) can now be calculated numerically. As a result, and must be adjusted to improve accuracy [57,58].

Results and discussion
In this work, numerical results are provided by focusing on physical parameters and constants in SI units, which are used in formulating the physical variables under study, and by using silicon as a material. The following physical values are used [39]: We divided graphical representations of several properties, such as size-dependent effects ( and ), parameter change of thermal conductivity 1 , and rotation Ω on the investigated thermalphysical fields, into three groups in order to analyze their effects.
The current work was prompted by a strong fit between the conclusions of nonlocal strain gradient theory and evidence from experimental research and molecular dynamics simulations. The investigation of wave propagation in the classical theory was established using the NSG theory and the DPL heat transfer theory, and the calculation results were compared with the theoretical and experimental results [59,60]. The nonlocal strain gradient elasticity model, unlike the conventional model, can demonstrate good agreement with experimental results.

Effect of nonlocal and material length scale coefficients
Nonlocal gradient theory (NSGT) has been used to capture the size-dependent effect of the nanobeam, as well as other specific cases such as classical theory (CET), nonlocal theory (NET), and strain regression theory (SGT). A stress-gradient parameter is also included to account for the stress gradient effect, as well as a strain-gradient parameter to account for the strain gradient effect. If = 0, the introduced nonlocal gradient theory (NSGT) model may be reduced to a nonlocal elastic model (NET), and if = 0, to a strain gradient model (SGT). • It is clear that the non-local parameter and the stress gradient modulus have a great influence on the interactions within the nanobeam.
• The dimensionless values of deflection and displacement grow in proportion to the nonlocal parameter values. From Figures 2 and 4, it is apparent that the magnitudes of NSGT, SGT and NET models for deflections and displacement have greater values than those achieved using conventional continuum theory (CET).
• When the parameter > , the NSGT model yields lower results than the nonlocal elastic model (NET). These facts suggest the nanobeam has a stiffness-reduction impact when the length scale parameter is lower than the non-local parameter, and when the parameter of the material length scale is higher than the non-local parameter, the nanobeam exerts a stiffness-reduction effect. Thus, it can be said that the non-local parameter has the effect of stiffness-softening while the gradient modulus has the effect of stiffness-hardening. Similar results were also found in several papers, such as in [60][61][62] on the dynamic response of beam-type structures based on NSGT.
• From the figures, it is clear that the difference between the NSGT and SGT models is less pronounced at lower non-local parameter values but becomes more pronounced at larger non-local parameter values. This pattern applies to both the CET and the NET. Indeed, regarding the reduction of the results due to the raising of the non-local parameter values, one can mention that the strain increases while stress is reduced in Eringen's differential formulation. This phenomenon causes a softening of the effect. Thus, an increase in the non-local parameter decreases the aberrations within the nanobeam. Compared with the published results, the collected results show good agreement, for example, with those of Jena et al. [22] and Zeng et al. [23].
• Also, when the parameter is smaller than the parameter ; the nonlocal effect predominates, which leads to a stiffness-softening effect of the nanobeam. When the parameter is greater than the parameter , the effect of the stress gradient takes precedence, causing the nanobeam to stiffen. Figure 2 shows that the deflection starts and terminates at zero values and fulfills the boundary conditions at = 0 and = . It is further shown that the aberrations recorded by the Euler-Bernoulli beam theory are usually slower than those predicted by non-local beam theories, indicating that the nanobeam is attenuated and softened when the size-dependent is taken into account.
• From Figure 3, we notice that the parameters and have a very weak influence on the distribution of temperature change . Thermal waves are continuous and smooth phenomena that reach a steady state depending on the phase lag factors, which means that heat transfers easily from one particle to another, causing the temperature to drop even more. This phenomenon is in contrast to the traditional theory of thermal elasticity, which predicts that the rate of propagation of thermal perturbation is infinite.
• Increases in temperature and thermal load cause a drop in beam stiffness, which reduces thermal and mechanical waves.
• As seen in Figure 4 it is evident that the magnitude of the displacement in the NSGT, SGT and NET models is less than that obtained using standard contact theory (CET).
• In contrast to the behavior shown in Figures 2 and 3, the bending moment behavior as shown in Figures 5. It is clear that the sizes of the NSGT, SGT, and NET models of the bending moment have lower values than those achieved using the conventional contact theory (CET). In other words, the bending moment decreases with the height of the nonlocal modulus.    • Finally, the classical theory, which does not take into account the size dependence and can yield results suitable for studying narrow beams, is well known. However, in non-local models, the effect of size-dependent mechanical behavior becomes significant and cannot be ignored. From theoretical calculations, we may conclude that the nanobeam can display either a hardening softening or a hardening-stiffness effect, depending on the relative amplitude of the non-local factor and the lengthscale factor of the material.
• This may explain why the stiffness-hardening effect of nanobeams is frequently found in experimental studies as in [61], but the stiffness-softening effect of nanostructures is more common as discussed in [62].

The effect of the variability of thermal conductivity
Thermal conductivity is an important property of a material that is often taken for granted. On the other hand, many experimental and theoretical investigations have shown that thermal conductivity is closely related to temperature change [63,64]. Depending on the general one-dimensional thermoplastics problem, Xiong and Guo [65] emphasized the effect of changing temperaturedependent properties on field variables. As described earlier, the thermal conductivity is assumed to be a linear function of the temperature change (see Eq (25)). A material with a high is considered an excellent conductor of heat, while a material with a low is considered a good thermal insulator. As a result, this parameter has a significant impact on both the working conditions in deep mines and the ability of subterranean vents to store thermal energy. It is also related to geothermal energy generation and radioactive waste disposal.
In this subsection, the effect of the variability parameter of thermal conductivity 1 on the nondimensional deflection, temperature spread, displacement, and bending moment is investigated. When the coefficient of thermal conductivity is dependent on heat, there will be As can be seen from Figure 6, the magnitudes of decrease as the parameter 1 increases. As illustrated in Figure 7, as the distance increases, the temperature drops, driving wave propagation. The temperature distribution is increased when the parameter 1 is reduced. Physically, increasing the parameter 1 improves the heat transfer process and raises the local temperature of the beam. It is also noted that with an increase in the parameter 1 , the temperature rises. The results were compared with those presented in [66] and it was found that there was agreement in the results.   Figure 8 shows that the displacement increases when 0.0 ≤ ≤ 0.1 for the greatest amplification and drops when 0.1 ≤ ≤ 0.6. In the last range 0.6 ≤ ≤ 1 of wave propagation, the displacement is traveling straight. The variability parameter 1 has a significant impact on the displacement. In Figure 9, the increase in the parameter 1 is designed to enhance the distribution of the bending moment .
The results demonstrate that the change in thermal conductivity has an effect that should not be ignored [67]. The mechanical and thermal distributions of the nanobeam reveal that the wave propagates as a velocity-limited wave in the medium and depends on the coefficient of change of thermal conductivity 1 [34]. We may deduce from the analysis that the thermoelastic behavior is affected by variations in thermal material characteristics, and these impacts are mostly focused on the distribution's peak values.

Effects of the rotation speed
As mentioned earlier, spinning nanobeams play an important role in many nano-devices, and properly modeling and studying them is a very challenging task. In this section, a study was conducted to investigate the effect of rotational speed Ω on the functions of the physical fields. Numerical calculations were performed when the values of the other effective parameters were constant. This work will use non-local gradient theory (NSGT) to elucidate the vibrational behavior of nanobeams for size-dependent beams, as there has been no precise inquiry into the nanobeams' rotation. It is found that the speed of rotation, elastic media, and non-local scaling factors greatly influence the bending vibration of the system. Figures 10, 11, 12, and 13 show the variance of the investigated fields for three distinct angular velocity values (Ω = 0, 0.1, 0.3). The rotation coefficient is zero (Ω = 0) in the absence of rotation, which is a special case of the current work. Figure 10 depicts the rotation Ω that affects the deflection . In the presence and absence of rotation (centrifugal force), this parameter was shown to have a significant impact on the deflection and variance in results. It is noticeable that increasing the angular velocity Ω reduces the deflection of the non-local transverse nanobeam. These findings and behavior are similar to those seen in [68,69]. In Figure 11, the effect of angular velocity on temperature change is examined and displayed. We note from the figure that the effect of rotation on the temperature change is very slight. The results of the preceding literature are compatible with those of the author, such as the equivalent results achieved in [37][38][39].
The effect of the angular velocity of rotation Ω on the variance of the displacement is shown in Figure 12. It is found that the speed of rotation Ω most prominently affects the curves reflecting the displacement field. The figure also displays that the displacement distribution decreases in certain periods with increasing rotation and rises in other periods along the axis of the beam. Figure  13 shows the changes in bending moment for different values of rotational velocity Ω. The graph indicates that rotation greatly affects the moment curves and that as the angular velocity Ω values decrease, the bending moment amplitude increases. Ultimately, the amplitude of the studied fields changes as the non-dimensional angular velocity increases. In other words, the speed of rotation Ω is most prominently responsible for the stiffening effect of the centrifugal force. The functional features of rotating structures with diverging effects of external stimulation can be improved based on understanding the influence of rotational speed in hybrid nano-generators to achieve better shaft motion.

Conclusions
In the current work, a nonlocal strain gradient beam model is constructed to examine the vibration response of the thermoelastic nanobeam, including its rotational effect. The Euler-Bernoulli beam theory coupled with the generalized thermoelastic model with phase delays and nonlocal straingradient theory, is used to construct the governing equations for the rotating nanobeam. It was taken into consideration that the thermal conductivity coefficient and the specific heat of the nanobeam are temperature dependent. The effect of nonlocal, material length scale parameters, rotation, and thermal conductivity variability on thermal and mechanical wave scattering is investigated. The following important outcomes are achieved: • The non-local gradient elasticity model showed results that were different from the results of the traditional model, and this is in agreement with the previous literature.
• Deflections and displacement are reduced as the non-local parameters increase. The strain gradient effect is a priority when the non-local parameter is greater than the length scale parameter of the material, which results in the annealing of the nanobeam.
• The difference between NSGT and SGT is less pronounced at lower non-local parameter values, but becomes more pronounced at larger non-local parameter values. This pattern applies to both the CET and the NET.
• The behavior of different fields is affected by changes in the thermal properties of materials, and these effects appear mostly in the peak values of the distributions.
• Depending on the relative amplitude of the nonlocal factor and the material length-scale parameter, the nanobeam may show the effect of stiffness-softening or stiffness-hardening.
• The fluctuation of thermal conductivity and its dependence on temperature change have a significant impact on different physical fields.
• The angular velocity of rotation is most prominently responsible for the strengthening effect of the centrifugal force.
• The effect of rotation on the temperature change is very slight. The fundamental effect of rotation on the behavior of different physical distributions can also be inferred from the obtained results and should be taken into account in production and design procedures.