The vibration analysis of a nanobeam due to a ramp-type heating under Moore-Gibson-Thompson theory of thermoelasticity

. Micro/nanobeams find extensive appreciation due to their employment as resonators. The current study investigates a thermoelastic homogeneous isotropic nanobeam exposed to ramp-type heating as a thermal loading. The governing equations of the nanobeam have been deduced in the context of the Moore-Gibson-Thompson thermoelasticity theory. The Laplace transform technique has been applied, while its inversions have been calculated by using Tzou numerical technique with approximation. When a rectangular thermoelastic nanobeam composed of silicon nitride is easily supported, the numerical findings have been confirmed and presented in the figures. It is worth noting that the ramp-time heat parameter has a substantial impact on all of the functions investigated. The parameter that distinguishes the Moore-Gibson-Thompson model from the standard thermoelastic model has a minor influence on temperature increment but has a considerable effect on lateral deflection, strain, and stress distributions.


Introduction
Nanomechanical resonators have been of great significance. For instance, in communications systems and mechanical electrometers or magnetometers, resonators made of silicone crystal material typically with semi-conductive techniques for the development of MHz or even GHz core frequencies may be used for the identification, force, pressure or acceleration sensitive loads [1][2][3]. Researchers are interested in understanding how the physical and mechanical properties of resonant parameters affect them. Sun et al. [4] test thermoelastic damping (TED) in the microbeam resonators based on generalized thermoelasticity theory with one relaxation time of Lord and Shulman (LS) model [5]. To solve combined thermoelastic equations, he used finite Fourier sinus transformation methods and Laplace methods, which extracted the expressions of fluidity and heat. The deflect vibrates in a simulated state mode when the thermal moment vibrates in the same mode as the deflection at first. In the form of the updated theory of couple tension, Hosseini [6] studied TED and microbeam resonator dynamical behaviour using the theory of thermoelasticity. For the solution of thermoelastic variations and the representation of the beam deflection and thermoelectric moment, the Laplace Transformation technology was adopted. However, the results of the modified pair stress theory on deflection and thermal moment are not evident in their work.
Kumar recently investigated and expanded on the role of the Sun and its impacts on deflection and thermal moment in microbeam resonators, taking into consideration the three-phase of thermoelasticity phase-lags ideas (TPL) [7]. We also address several experiments on thermoelastic beam vibration, as references in the literature indicate [8][9][10][11][12][13].
Given the foregoing description of beam resonator thermomechanical testing, this study intends to investigate the thermal and dynamic behaviour of simply by using MGT thermoelasticity theory, as recently mentioned by Quintanilla [14]; this is the Moore-Gibson-Thompson (MGT) equation. In this paper, the heat conduction equation is described. In reality, this modified heat conduction equation is reached when the relaxation parameter is implemented in the Green-Naghdi type-III model (GN-III) [15,16]. The coupled thermoelastic problem's governing equations are derived based on the MGT, LS, and GN-III models [16]. Through comparing the findings achieved under the MGT model with those of LS and GN-III models, vibration responses of the deflection and the thermal time are analyzed [16].
In the current work, ramp-type heating as a thermal loading will be applied on thermoelastic isotropic nanobeam. A simplified form of thermoelasticity theory will be used for Moore-Gibson-Thompson. The transform technique for Laplace will be used and its inversions will be computed using the approximation numerical Tzou procedure. The numerical results for a rectangular nanobeam made of silicon nitride will be applied to validate the results.

Basic equations
A homogeneous thermal isotropic, thermal viscoelastic solid in the Cartesian coordination category Kelvin-Voigt was considered at a uniform temperature and initially undeformed. The fundamental differential motion and heat conduction equations for generalized thermoelasticity are based on the thermal control legislation of Non-Fourier.
Without body forces and heat sources, the displacement components and absolute temperature are given by [17].
-If * ≠ 0 and ≠ 0, we obtain the Moore-Gibson-Thompson (MGT) model [14]. The deformation-displacement relations are in the form: where , = , , , = 3 + 2 , is the coefficient of linear thermal expansion, is the density, is the thermal relaxation time, , Lamè's parameter in the usual case, is the Kronecker delta symbol, is the thermal conductivity, * is the conductivity rate parameter, and is the specific heat at constant strain, and for stable solution > * [14,19].

Problem formulation
We consider a thin thermoelastic nanobeam of length ℓ 0 ≤ ≤ ℓ , width − /2 ≤ ≤ /2 , and thickness ℎ −ℎ/2 ≤ ≤ ℎ/2 with small flexural deflections. The In the state of equilibrium, the nanobeam is unstrained, unstressed, and without any mechanical damping. Moreover, the reference temperature is everywhere [6]. In this work, the well-known Euler-Bernoulli equation [20] is adopted, therefore any plane cross-section that is at first perpendicular to the axis of the beam stays flat and perpendicular during bending to the neutral surface as in Fig. 1.

Formulation of the problem by using Laplace transform
The Laplace transform will be applied for Eq. (21) and Eq. (22), which is defined by the following formula: The Laplace transform inversions could be obtained by using the following formula: where is an imaginary number unit and Re is the real part. For faster convergence, many numerical experiments have concluded that the value satisfies the relation ≈ 4.7 Tzou [22,23]. Hence: where = + + ⁄ and = + . We can re-write Eq. (26) and Eq. (27) to be in the forms: where = ⁄ . After elimination, we get: where = + , = , and = . The general solutions of Eq. (30) and Eq. (31) will take the forms: where , , and are unknown parameters to be determined, and the roots of the characteristic equation are ± , ± , and ± which satisfy the following characteristic equation: To get the parameters, we will substitute from Eq. (33) and Eq. (34) into the Eq. (30). Hence, we have: which gives: To calculate the parameters , = 1, 2, 3, we must use the boundary conditions. We consider that the first end of the beam = 0 is thermally loaded by a function and simply supported as follows: where is constant. The second end of the beam = ℓ is simply supported, and zero temperature increment: Applying the Laplace transform, we get: Then we obtain the system of linear equations shown below: A mathematical software MAPLE 2020 has been used to solve the above system and obtain the constants , , .
Hence, we get: The lateral deflection is: which complete the solution in the Laplace transform domain.

Numerical Results
For the thermal loading, we considered the ramp-type heating, which takes the form [22]: By applying Laplace transform on the above equation, we get: The figures will be prepared by using the non-dimensional variables for nanobeam as ℓ = 0.5, = 1.0, = ℎ/4, and = 1.0.

Discussion
Figs. 2-5 represent the temperature increment, lateral deflection, strain, and stress distributions for the three studied models: Lord-Shulman (LS), Green-Naghdi type-III (GN-III), and the newest one Moore-Gibson-Thompson (MGT) when = .  Fig. 2 shows that the three models offer the same behaviour and values for the temperature increment, where the three curves are very close from the beginning up to the end. Fig. 3 shows that the parameters and have significant effects on the lateral deflection distributions. The absolute values of the lateral deflection take the following order: Fig. 4 shows that the parameters and have significant effects on the strain distributions. The absolute values of the strain take the following order: Fig. 5 shows that the parameters and have significant effects on the stress distributions. The absolute values of the strain take the following order:  We noted that an increase in the value of the ramp-time heat parameter leads to a decrease in the temperature increment, lateral deflection, strain, and stress of the nanobeam.
Thus, the ramp-time heat parameter has significant effects on all the studied functions and the thermomechanical waves' propagation through the nanobeam. Fig. 8. The strain for different models when < Fig. 9. The stress for different models when <

Conclusions
The current study looks at the temperature increment, lateral deflection, strain, and stress of a simply supported nanobeam resonator heated using a ramp. In terms of dimensionless temperature increment, lateral deflection, strain, and stress functions, the governing equations for three distinct coupled thermoelasticity theories (LS, GN-III, and MGT) of nanobeams have been developed. The following are the major achievements of this work: 1) The temperature increment under the three studied models is approximately the same.
2) The absolute values of the vibration, strain, and stress in the context of the GN model are greater than its values under the LS model, and they have smaller values under the MGT model.
3) As the uniform load on the upper surface of the nanobeam rises, so do the deflection amplitudes and thermal momentum. 4) As the ramp-time heat parameter is increased, the temperature rise, lateral deflection, strain, and stress decrease.
5) The relaxation parameter has a considerable effect on the vibration of the lateral deflection, strain, and stress distributions. 6) In future work, this work could be extended to include other types of thermal and mechanical loading and different types of nanobeams such as circular nanobeams.