1. Introduction
Many manufacturing processes are usually exposed to moving heat sources such as welding, cutting, surface laser treatment, and some tribological problems [
1,
2]. Two important applications of a laser pulse beam, i.e., laser texturing and laser piercing processes, are illustrated in
Figure 1. Thanks to its simplicity and effectiveness, “beam theory” is an essential tool in solid mechanics that enables the designers and engineers to evaluate the behavior of different structures and systems from macro- to micro-scale dimensions. Beam theories are commonly used in the pre-design stages because they offer a substantial understanding of the performance of structures. This is true even though more sophisticated theories for stress investigation of complex structural systems are now accessible. Several theories concerning beams have been developed, all of which are predicated on various assumptions and produce varying levels of precision. The Euler-Bernoulli beam model is widely recognized as one of the field’s most fundamental and significant concepts. Microscale beams use microscale beams, many microstructure devices, systems like sensors and actuators, microelectromechanical systems (MEMS), and very thin electronic gadgets [
3,
4]. Understanding the thermal behavior of such structural components in the manufacturing processes is critical for predicting the performance of such sensitive devices.
The response of beams that are on two-parameter elastic foundations, which have a variety of features in design and related industries, has been the subject of a significant amount of research in the past few years, and various theoretical works have been presented recently that have analyzed this behavior and attitude [
5]. Many different structures and systems can be described as beam structures that foundations support. An essential area of geotechnical engineering is the study of structural foundations that are supported by soil media. It is both necessary and practical to idealize the soils’ responses to external loads because the majority of naturally existing soils have complicated stress-strain properties.
Geotechnical assessments of beams that have been resting on elastic soil layers or pile behaviors caused by horizontal stresses, etc., have all benefited from using the Winkler model (1867). Based on “Winkler’s hypothesis”, which argues that the deflection at every point on the surface of an elastic continuum is proportional exclusively to the load that is being applied to the surface, this model was developed. An independent set of vertical linear springs is considered to compose the continuum in this mechanical model, which follows from the premise in question. Disadvantages arise, however, from the fact that the displacements that are in adjacent cells do not blend smoothly [
6]. By introducing shear interaction between the nearby Winkler spring elements, Pasternak (1954) revolutionized the field. By using the words “Pasternak model” to describe a generalized foundation model, Kerr (1964) viewed Pasternak’s mechanical model as a shear layer on the Winkler framework. The Pasternak model, referred to as PM from here on out, is applied in a similar way. For the past many decades, the approaches for obtaining the parameters have been addressed because the two-parameter models inevitably need both values for their use in practical applications. Despite not providing the parameters’ exact values, Pasternak himself recommended utilizing plate loading tests to evaluate the two parameters.
The interaction between the elastic substrate and the main structures and their influence on the system’s thermal behavior may be important in several engineering fields. Despite their significant structural design challenges, elastic-founded structures are frequently utilized in modern architecture, civil engineering, and electronic devices. Consequently, in recent years, numerous attempts have been conducted to explore elastic structures’ dynamic and thermal behaviors, where different foundation models were described in detail [
7,
8,
9,
10,
11,
12].
Any temperature variation in a physical system allows heat to move from one temperature range to another. Until the temperature of the device is steady, this kind of transmission occurs. Although, over the past two decades, scientists have paid considerable attention to Fourier’s thermal conduction equations for modeling the thermal analysis of structures, various phenomena continue to occur in everyday life requiring specific consideration and interpretation. For instance, the fundamental Fourier law appears to be wrong when extremely quick events and small scales are taken into account, and complicated theories are required to characterize the thermal conductivity of systems in science, accurately. The extended thermoelasticity hypothesis considers the effects of the interplay between temperature and strain speed. However, the following coupling formulas are hyperbolic.
Figure 1.
(
a) Schematic of laser texturing process using pulsed laser [
11] (
b) and the laser cutting process [
12].
Figure 1.
(
a) Schematic of laser texturing process using pulsed laser [
11] (
b) and the laser cutting process [
12].
Consequently, the paradox of the unlimited rate of heat waves in the conventional coupled model will be resolved. Lord and Shulman [
13] created a more advanced generalized heat conduction model that included a novel heat transfer law to overcome the traditional Fourier law flaw. The updated rule covers heat transport and its partial derivative which is connected to time. Three thermoelastic theories were developed by Green and Naghdi [
14,
15,
16] for homogeneous elastic materials. These theories for the coupled equations aim to transform type I to conventional heat conduction concepts, where the associated systems are linear. Thermal waves can propagate at a restricted rate in the linearized forms of types II and III. Several investigators [
17,
18,
19,
20,
21] have presented new generalized theories to improve the standard Fourier’s law by employing higher-order time derivatives.
The dynamic thermal behavior of an elastic medium that is exposed to heat sources has attracted more attention in the past decade. By developing an equation with two unique phase lags, one which is in the temperature gradient and the other which is in the heat flux vector, Tzou [
22] devised a two-stage heat transfer law in order to determine the amount of time that passes before the reactions occur due to the effect of the microstructure throughout the duration. The theory of three-phase delays that was suggested by Roychoudhari is also one of the theories of thermoelasticity [
23] that is developed in this context.
In recent years, studies have given greater consideration to the Moore–Gibson–Thompson (MGT) equation, which is used in the thermal study of elastic materials. This model can be derived by applying a third-order differential equation that considers the importance of different aspects of fluid mechanics.
Although many generalized thermoelastic theories have been presented to solve the problem of heat wave propagation at infinite speeds, some generalized theories, such as the Green and Naghdi model of the third type (GN-III) [
15], in addition to Roychoudhari’s model [
23], show that the main component depends on infinity in some components of the energy equation and the revised Fourier law. As a result, relying on the solutions of these models often fails [
24]. Because of this problem, Quintanilla [
25,
26] developed a new and better heat transfer model by adding the effect of relaxation time to the third type of the GN-III model. Based on a system that includes the Moore–Gibson–Thompson (MGT) equation, Quintanilla [
25,
26] introduced the improved model for heat conduction. Many researchers have recently utilized different thermoelastic models by considering specific assumptions and hypotheses [
27,
28,
29,
30,
31,
32,
33,
34,
35,
36,
37,
38,
39].
As a result of the complexity of the model, the issue of microbeams resting on three-parameter foundations while being exposed to initial thermal stress has attracted less considerable attention. In addition, the thermoelastic behavior of microstructures that are constructed on visco-Pasternak foundations involving axial tension and thermal coupling has not been previously investigated, and the current research study is the first effort to concentrate on the MGT thermoelastic model. In reality, several different foundation models are explored, along with the use of numerical and analytical methods, to perform free vibration studies of most buildings. Therefore, adding thermal coupling to the equation is a new way to deal with the vibration problems of these kinds of buildings that are on elastic foundations.
This paper will be structured as follows. The fundamental equations and fundamental ideas of thermoelasticity are introduced in
Section 2. In
Section 3, the thermoelastic coupling and its influence on the transient behavior of microbeams that are on a three-parameter viscoelastic foundation that is illuminated by a very short laser pulse are considered. The strategy of the Laplace transform technique was used to solve the equilibrium system, as is shown in
Section 4, while the response of the studied field variables is analyzed and discussed in detail in
Section 5 through some sub-sections. In
Section 6 of the article, the most important conclusions that have been obtained are summarized.
2. Basic Equations of MGT Thermoelasticity
According to the theory of thermoelasticity, the basic system equations for homogeneous solids can be expressed as [
23]:
where
represents the stress tensor,
denotes the displacement components,
,
is the absolute temperature,
indicates the ambient temperature,
symbolizes the strain components,
stands for elastic parameters,
represents components of the thermal conductivity tensor,
denote the coupling coefficients,
indicates the specific heat,
,
denotes the heat source,
is the heat flow, and
represents the substance’s density.
A parabolic equation is produced when the energy conversion Equation (4) and Fourier’s law (1) are combined. This allows the waves to spread at nearly unlimited rates. This conclusion is physically unlikely, and many scientists have claimed that this model is incompatible with the idea of extremely fast thermal conductivity. For this reason, the MGT equation was employed in this investigation to propose a revised law of thermal conductivity to solve this physical discrepancy. In the proposed model, this problem is solved by producing a hyperbolic partial differential equation that represents the motion equation that is based on the displacement vector and another one that represents the thermal conductivity equation that is based on the temperature change.
The modified MGT thermoelastic equation can be written as [
25,
27]:
The Moore–Gibson–Thompson model of thermoelasticity, abbreviated as MGTTE, can be derived from four separate models of classical and generalized thermoelasticity, depending on the specific conditions of a particular case and the thermal constants. Here is a list of previous thermodynamic models that can be derived:
When and , it is possible to use the conventional thermoelastic theory (CTE).
The framework for Lord and Shulman (LS) may be produced if the value of is equal to zero.
The idea of Green-Naghdi (GN-II) may be derived if and
The model of Green-Naghdi (GN-II) of type III (GN-III) through the usage of .
3. Statement of the Problem
The illustration in
Figure 2 displays a very thin, flexible microbeam with a constant cross-section
that has been excited by an ultrashort laser. The
directions specify the beam axis, while the
and
variables specify the microbeam width and thickness, respectively. The microbeam is exposed to small-amplitude bending oscillations along the
-axis, resulting in a deflection that conforms to the linear Euler–Bernoulli theory (EBT). As shown below, the EBT with small deflections
is used to figure out the displacements and stresses.
The Pasternak type that has two parameters is the greatest straightforward extension of the Winkler model what has a single parameter. It functions by attaching the ends of the springs to a plate that employs a shear layer, which is thought to be incompressible and simulates the shear interaction that occurs between the various spring components. The microbeam that is under consideration is contained within a homogeneous, three-parameter, viscoelastic substrate. The Pasternak (shear) basis modulus , the damping modulus of the viscous medium and, finally, the linear Winkler modulus create the foundation model.
Taking into account the fact that there is infinite contact between the microbeam and the material, the reaction behaves according to a basic three-parameter visco-Pasternak model, according to the relationship.
This model may be reduced to the visco-Winkler model by setting
equal to zero. It is possible to eliminate the viscoelastic term by setting
. The is one constitutive equation for a one-dimensional thermoelastic medium, which is provided by:
where
,
denotes thermal expansion, E represents the modulus of elasticity,
denotes the Poisson’s ratio, and
indicates the thermal stress. The equation of motion for the microscale beams can be stated as follows when using Pasternak’s basis and Hamilton’s principle and taking into account the presence of an axial force:
where
represents the moment of bending, which may be calculated using the following equation:
When Equation (10) is put into Equation (12), the bending moment may be calculated as
in which the cross-section inertia moment is
, and
stands for the thermal moment. Equation (13), when it is inserted into Equation (11), yields
The following is a generalized form of the Moore–Gibson–Thompson heat transfer Equation (6), which may also be written as [
25,
27]:
The temperature increase for the present microbeam is supposed to vary in the following form:
Consequently, we obtain this by inserting Equation (16) into Equations (13) and (14):
After multiplying by
, the result of doing an integration concerning z along the microbeam thickness using Equation (15) is shown below:
We now offer the non-dimensional variables that are stated as follows.
When it comes to these non-dimensional values (after removing the primes), the equations that control the situation consist of the following:
where
,
, and
each stand for the amount of time it takes the elastic substrate to mechanically relax due to the viscoelasticity, the Winkler model, and the shear stiffness of the elastic substrate, respectively. The equation for calculating the energy absorption rate, denoted by
can be imposed as follows: [
40,
41]:
In which and represent the laser and surface effects, respectively, denotes the duration of a laser pulse, and is the laser penetration depth.
Moreover, initial conditions are expected to be:
We will proceed to the assumption that the microbiome meets the following boundary requirements in its entirety at its two ends:
There are several sinusoidal pulse actuations with different pulse widths. Pulsed thermal shock generators are used in several different applications to achieve higher machinery efficiency levels. At the end
, the microbeam is considered to be thermally stressed by a harmonically varying heat.
where
denotes the Heaviside step function,
represents the frequency of the system of the vibration that is caused by heat, and
is a constant. It should be noted that if the frequency of the thermal vibration
, the boundary is presumed to represent a thermal shock.
It is also supposed that at
, the function
meets the following requirements:
4. Solution in the Transformed Domain
The Laplace transform approach is a very effective and helpful tool that has been used in mathematical applications in various areas of science and engineering. It is utilized to solve differential and integral equations. However, identifying the reversal of a specific transformation is a difficult challenge to solve when one is using this strategy. To address the problem, we will utilize the Laplace Transform, which is described by the following formula:
When we apply the Laplace transform to Equations (16)–(19) together with (21), we find that:
where
By eliminating
from (31) and (32), we arrive at the following:
where
The governing Equation (36) can also be written as follows:
where
,
coefficients are equal to the roots of the equation:
The homogenous solution of Equation (38) is as follows:
where
. In Equation (40), the elements
,
will be computed by using the conditions of the system. Likewise, we may obtain the following DE by removing
between Equations (31) and (32) as:
where
and its solution may be written as:
where
It is possible to create a solution for the functions
,
,
, and
in the Laplace field with the assistance of Equations (40) and (42) by writing the solution as follows:
The energy that is put into a material in the form of work is referred to as the material stress energy. One may refer to the internal work that is done by the stress products as they travel through related distortions as stress energy. The material is expected to maintain its flexibility while working on it, ensuring that no energy is wasted, and no permanent deformation is caused by the material’s production. The total stress energy or work done by the microbeam can be represented as:
In the domain of Laplace transformation, the conditions (27)–(29) will have the form:
Substituting Equations (40) and (42) into the boundary conditions (49) yields the following results:
As a result of the governing equation being solved, we may deduce the indeterminate constants , .
The Riemann sum approximation technique is employed to get the numerical values of the studied field variables [
42]. Many algorithms are used for this purpose. In the current study, the following relationship will be utilized to convert the functions from the Laplace field to the time realm:
In the above relation, the abbreviation
stands for the imaginary number,
stands for the real component, and
stands for the number of truncated terms that were used in the numerical solution. It is important to point out that the accuracy of the Riemann sum approximation depends on the value of the term t. In contrast, the number of iterations governs the truncation error. Both of these facts are worth emphasizing. However, several numerical investigations have shown that the value
delivers the greatest results for a speedier convergence [
43]. These tests were carried out in a variety of ways.
5. Discussions and Numerical Results
In this part of the article, a comparative evaluation of the outcomes that have been obtained in this work and those of the previous literature will be included to validate the proposed model. For this reason, the influence of several parameters on the behavior of non-dimensional displacement , temperature change , axial heat stress , and deflection will be discussed.
The characteristics of the material are described as follows [
44]:
In addition to the above physical values, the microscopic beam is assumed to have a length-to-thickness ratio of
, where h is the beam thickness. Also, in the numerical calculations, the parameters of the microbial beam dimensions are taken into account, taking the values
,
and
. The results will be presented tabularly and graphically in
Figure 3,
Figure 4,
Figure 5,
Figure 6,
Figure 7,
Figure 8,
Figure 9 and
Figure 10 in the direction of the axis of the beam between the values of 0 and 1.0. The numerical results indicate that the distributions of different physical quantities are influenced not only by time and spatial coordinates but also by different effective parameters such as laser pulses and beam foundations. Both the graphs and numerical calculations will be organized into four different groups.
5.1. Results Validation
To verify the numerical results, a comparison was made between the formulation technique and the solutions that appeared in [
29,
30,
31,
32] with the thermal elastic behavior of the micro-beam using extended thermoelastic theories. The conduct of thermo-mechanical waves of various magnitudes is found to be very similar when comparing it to the results that were obtained in the present work with those that are discussed in the literature [
29,
30,
31,
32]. As a result, the model that is presented in this article is consistent with the results of previous models, but it is distinct from in that this one provides more accurate representations of the current world. It can be concluded that the inclusion of the relaxation time in the suggested heat transfer equation decreases the propagation of the thermal waves within the microbeams. Also, by making comparisons between the results of the proposed current framework with those in the case of Green and Naghdi theories [
16], the importance of the suggested version is shown in that the problem of infinite speeds of heat waves has been solved.
5.2. The Influence of visco-Pasternak Foundation Characteristics on the Results
As far as the author is aware, there has been no research that has been conducted examining the surface effects of microbeam vibration analysis when it is under the influence of a visco-Pasternak elastic basis and considering the effect of temperature change based on the generalized view of thermoelasticity. Therefore, in this section, the interest is in investigating how microbeams, which are grown under the influence of a visco-Pasternak elastic foundation, behave when they are exposed to laser pulses. These changes will be explained in
Table 1.
Table 1 shows the variation of the different dimensionless field variables concerning the dimensionless Winkler modulus
for different values of the shear stiffness of the elastic substrate (Pasternak’s modulus)
, as well as the effect of mechanical relaxation time due to viscosity (viscosity) or the so-called damping coefficient
to show the influence of the Winkler parameters and Pasternak on the behavior of a thermoplastic substrate. The Moore–Gibson–Thompson (MGTTE) generalized thermoelastic model was considered to investigate and obtain the numerical results of the physical distributions. For comparison purposes, we will also take the values of other parameters as constants, such as
. At single values on the phase of the beam axis at point
, numerical calculations are performed with different values of the foundation parameters. As a special case of the current work, in setting the foundation parameters
, the present results lead to the calculations of the problem of a thermoelastic beam that does not contain a flexible visco-Pasternak base.
According to the numerical values in
Table 1, the obvious influences of the Winkler and Pasternak parameters
and
on the values of deflection
are shown. Changes in the visco-Pasternak foundations
and
have a significant effect on both bending moment
and deflection
. Increasing the values of the fixed base factors reduces the absolute values of the deviation w for each analyzed case. This is due to the fact that the system becomes more robust as the microbeams become more rigid. This occurrence and behavior were previously discussed in the paper by Wattanasakulpong and Chaikittiratana [
45], and by comparison, it was found that there was great agreement in the results. Looking at
Table 1, it can be noticed that raising the values of the dimensionless Winkler foundation (
) and Pasternak foundation (
) moduli leads to an increase in the amount of both dimensionless heat stress
and the bending moment M in the microbeam. The reason for this occurring is because a rise in the values of
and
leads to a rise in structural rigidity and, as a consequence, the stability of the microbeam. On the other hand, the displacement
and the temperature change
are less influenced by the variations in the coefficients of the viscous basis of Pasternak. This is because the absolute value of the displacement
and the amount of temperature change are very sensitive to changes in the foundation as well as the viscous basis. The absolute value of the displacement
and the amount of temperature change decrease by raising the values of the basis factors
and
.
We also see from the table that the quantities of the various distributions of the physical areas of research decrease when the Winkler and Pasternak factors and remain unchanged, and the viscous damping coefficient undergoes an increase. Additionally, when there is a change in the viscous damping coefficient , there is a corresponding minor change in temperature as well as displacement. The absolute value of the deflection w increases in proportion to the decrease in the viscosity coefficient . The non-dimensional stress and the bending moment both grow as the viscous damping coefficient and other foundations are increased.
On the other hand, and by checking with the results that are in
Table 1, it becomes clear that the quantities of the different distributions of the studied physical fields change dramatically when the Winkler and Pasternak coefficients
and
remain unchanged, and the viscous damping coefficient
is subject to different values. In addition, when there is a rise in the value of the viscous damping coefficient
, there is an increase in the amount of both temperature change
, displacement
, and deflection
. When the viscosity coefficient
is lowered, the absolute values of the heat stress
and bending moment
rise to a higher level. The amount of temperature change grows when there are higher values of the viscous damping coefficient
and with the stability of other foundations. For this reason, this observation must be considered when manufacturing and designing small-sized devices.
It has been investigated that the reduction of the dimensionless thermophysical domains can be achieved by increasing the modulus of the viscoelastic structural damping of the beams. This is due to the fact that an increase in the structural damping factor causes the structure to be less rigid. The effect of the viscoelastic structural damping coefficient is amplified when the Winkler and Pasternak coefficients and are fixed.
5.3. The Influence of the Duration of the Laser Pulse
In recent years, investigations of laser applications have been necessary to gain knowledge of the internal structural properties of many materials. As a result, many contemporary applications of laser pulses have emerged in many fields, such as physical sciences, mechanical engineering, medicine, and others. When considering how extended thermoelastic concepts affect beams, the thermal influence of a non-Gaussian laser on a flexible material that is employed as a heat source is the most significant thing to consider.
This discussion subsection studies the phenomenon of small-scale heating beams using a pulsed laser. The microbeams in this study consist of aluminum material and are heated using a pulsed non-Gaussian laser for a few seconds . This allows the resonance to be damped by the thermal elasticity within the material. The basic energy dissipation process converts kinetic energy into stable thermal energy. This conversion occurs when temperature and pressure interact with each other within the microbeams.
Figure 3,
Figure 4,
Figure 5 and
Figure 6 show the results of an investigation into the influence that the duration of the laser pulse
has on the pattern of temperature change and deflection distribution as a function of the distance
. There is no variation in the values of the vital parameters
,
,
, or
. The main purpose of this part of the discussion is to apply the proposed thermal elastic MGT model to investigate the heat source-induced heat elastic wave behavior to a study where a non-Gaussian laser beam was used to heat the surfaces of the microscale beams. With the lengthening of the laser pulse, the intermolecular distances lengthen, and the temperature of the medium increases, thus reducing the intermolecular pressure. As the microbial beams sit on a flexible Winkler foundation, a higher laser pulse coefficient
results in lower temperature and deflection.
Figure 3.
Variation of deflection against laser pulse duration .
Figure 3.
Variation of deflection against laser pulse duration .
Figure 4.
Variation of temperature against laser pulse duration .
Figure 4.
Variation of temperature against laser pulse duration .
Figure 5.
Variation of displacement against laser pulse duration .
Figure 5.
Variation of displacement against laser pulse duration .
Figure 6.
Variation of thermal stress against laser pulse duration .
Figure 6.
Variation of thermal stress against laser pulse duration .
Figure 3 shows how a variety of laser pulse properties can affect the amount of microbeam deflection
. According to this figure, the fine beam endpoints do not show any deflection, which is in contrast to the deflection that is observed elsewhere along the axial axis between the two ends. This result is consistent with the boundary conditions that have been considered for the problem. It is also shown that the deflection
drops sharply to its lowest value at point
near the application of the heat source and the harmonic convection. After that, the deviation begins to rise gradually, making it close to zero at the point
, before fading out completely. A significant increase in the aberration contrast w is shown, referring to the increase in the laser pulse length parameter
, as shown in
Figure 3. This is due to the fact that the ultrashort laser pulse leads to very fast movements in the structural elements of the material, which in turn lead to significantly large inertial forces and, as a result, an exceptionally high vibration and, thus, deflection
.
The influence that the length of the laser pulse
has on the temperature distribution
is seen in
Figure 4. According to the results that have been shown, the temperature distributions become less extreme as the wave travels a given distance and continues to travel in the same direction as the wave travels. The thermoelastic beam is sent through the medium at a finite speed and does so in accordance with the temperature profiles. This discovery demonstrates that the experimental and theoretical predictions that are made by this novel model (MGTTE) correspond quite well with one another for metallic materials. The graphic demonstrates that an increase in the laser pulse duration
that was applied to the elastic medium results in a temperature rise. Because of the quickly varying contraction and expansion, those materials that are susceptible to heat transmission through conduction might experience temperature fluctuations [
46]. Because pulsed laser technologies are so prevalent in material processing, nondestructive testing, and characterization, this particular mechanism has garnered significant attention [
47].
Figure 5 is a graphic that illustrates the influence that the parameter
has on the displacement
sensitivity of an embedded microbeam that has been treated with an ultrafast laser. The information that is shown in this figure leads one to the conclusion that the impact of increasing the duration of the laser pulse results in an increase in the values of displacement
. The displacement begins with positive values and gradually decreases until it crosses the
-axis, at which point it begins to take positive values and eventually approaches zero in three different instances.
Figure 6 illustrates how the fluctuation in the thermal stress is affected by the parameter
. It is important to keep in mind that the amount of stress
will rise in proportion to the length of the laser pulse. According to [
48,
49], the laser pulse parameter
is an important variable to consider when analyzing the dispersed features of the microbeam. When the values of the
laser pulse are increased, this increases the values of the studied fields, which can be observed at the locations where there is a peak of the curve.
5.4. Comparative Analysis of Several Thermoelastic Theories
In these paragraphs, the outcomes of the mathematical modeling of the problem are presented, and some comparisons with previous research are made to verify the framework that is recommended in this article. Also, the most significant deductions that can be derived from the discussions are presented. If all other effective factors remain constant, then the studied fields along the -axis of the beam that is generated on a viscoelastic basis can be investigated if several different thermoelastic theories are used.
As indicated in
Section 2, some preceding theories of thermoelasticity can be obtained from the current model. First, if the thermal parameters vanish, i.e.,
, then the coupled theory of conventional thermoplastics (CTE) can be obtained, which predicts the infinite velocities of heat waves. Secondly, if the parameter
is omitted, Lord Shulman’s theory (LS) can be obtained, including the relaxation time. Third, the type III GN-III model of Green and Naghdi is accessible when the relaxation time parameter
disappears, while the type II GN-II model is formed when the first term on the right side of the generalized heat equation is omitted (6). Finally, in the case of the thermoelastic model (MGTE), which includes the Moore–Gibson–Thomson equation, the thermal parameters are
. In addition to
Figure 5,
Figure 6,
Figure 7 and
Figure 8, numerical results will also be shown in
Table 1,
Table 2,
Table 3 and
Table 4 so that researchers in this field can easily compare their results to them. It is shown how the laser pulse heating changes the thermomechanical behavior of the microbial beams with the change of the thermoelastic model that is used by taking various values of the thermal constants
and
at the same values for time
. It is clear from the numerical outcomes that the constant thermal factors
and
have a prominent effect on each field of physical journals.
After the beginning of the study of the thermoelastic microbeam, the numerical results in the case of the coupled model (CTE) and the generalized theories of thermoelasticity (LS, GNII, GNIII, and MGTE) are remarkably equivalent to one another. According to the CTE hypothesis, the reaction of the issues inside the elastic material may be very unalike. In the CTE model, thermal waves move at an unlimited speed, in contrast to the considered advanced models, which travel at a restricted rate. At the beginning points, such as for temperature and deformation , as well as at the peak places, such as for heat stress and thermal deflection , the differences between the models can be observed.
The numerical results in
Table 2,
Table 3,
Table 4 and
Table 5 and
Figure 7,
Figure 8,
Figure 9 and
Figure 10 show the differences between the predictions that were made by the new MGTE theory and those made by other thermoelastic models. According to the numerical analysis results, it is noted that the values of the physical fields that were predicted by the GN-III and CTE theories are much higher than the values that were predicted by the MGTE theory, as well as other generalized models of thermal elasticity. In addition, it is quite clear that the responses of the generalized thermoelastic model of type GN-III converge with the responses of the classical thermoelastic model (CTE), which does not decompose into heat as quickly as other generalized thermoelastic theories do. This shows the shortcomings of the GN-III model, which predicts the spread of thermal waves at unlimited speeds like the traditional model. When one is looking at the tables and figures, it can be seen that the effect of the relaxation coefficient in the MGTE and LS theories show a slight decrease in temperature change behavior. This displays the precision of the suggested MGTE model, which addresses the contradiction in the GN-III thermoelasticity model. The obtained results provide further evidence that the values of different physical domains will converge and behave similarly according to the LS, GN-II, and MGTE theories.
The positions and behaviors of viscoelastic materials in response to variations in temperature and other physical variables are illustrated in both the images and the table. These differences were examined in the MGTE and LS models of thermoelasticity. Although there are slight differences in magnitude, in most cases, both concepts predict roughly the same behavior patterns. These results are consistent with the results and concepts that are presented by the generalized thermoelastic theory. From the previous observations, it is clear that there is a great agreement with that which Abouelregal et al. [
27,
29,
31] said about this behavior in the case of different models.
Figure 7.
The thermal deflection for various thermoelastic models.
Figure 7.
The thermal deflection for various thermoelastic models.
Figure 8.
The temperature change for various thermoelastic models.
Figure 8.
The temperature change for various thermoelastic models.
Figure 9.
The displacement for various thermoelastic models.
Figure 9.
The displacement for various thermoelastic models.
Figure 10.
The axial stress for various thermoelastic models.
Figure 10.
The axial stress for various thermoelastic models.
6. Conclusions
This research has investigated the microbeams’ thermoelastic response while they were resting on a visco-Pasternak foundation using a novel model that was composed of three parameters. The Moore–Gibson–Thompson equation served as the foundation for the improved thermoelastic theory that was taken into consideration. Lord Shulman’s thermoelastic model (LS) was the starting point for the derivation of the heat equation. We established the mathematical model of a microbeam that was simply supported and sitting on a visco-Pasternak substrate with three parameters after being subjected to an ultrafast laser pulse. Thermal deflection, temperature change, axial displacement, bending moment, and axial thermal stress of the considered microbeam were all measured and analyzed. It was also found that the laser pulse duration parameter affects how the structure responds to changes in temperature.
According to the results of our research, it was demonstrated that the basis parameters, in particular the viscosity basis parameter, had a substantial influence on the thermal behavior of microbeams. Also, an increase in the parameters of the Winkler or shear institution lead to a decrease in both the axial displacement towards the first edge of the microbeams and the amount of microbeam deflection. On the other hand, the deflection and displacement values tended to move higher as the viscosity parameter became more difficult. Foundations are essential to any structural component, as they prevent the failures that are caused by mechanical stress.
In addition, a rise in the visco-Pasternak basis modulus reduced the microbeam’s dynamic response, making the beam balancing process occur significantly faster. In addition, both the amplitude of the physical domains and the reaction time decreased dramatically when the viscous damping coefficient increased. This is because of how viscous damping works. Also, since the laser pulse duration increased the dynamic reaction, this influence is important and must be considered when studying micro- and nanostructures. Due to the effect of the relaxation variable, both the LS theory and the MGTE theory predicted that the temperature and other mechanical fields would decay at a more gradual rate. The physical data that are presented in this study can be useful in designing and developing applications using the vibrations of micro-sensors and microresonator systems. In conclusion, it should be noted that the results that are reported in this study are intended to be utilized as criteria for the confirming grades that are obtained by other mathematical approaches and to assist in designing MEMS machines.