Higher-Order Thermo-Elastic Analysis of FG-CNTRC Cylindrical Vessels Surrounded by a Pasternak Foundation

This study analyses the two-dimensional thermo-elastic response of functionally graded carbon nanotube-reinforced composite (FG-CNTRC) cylindrical pressure vessels, by applying the third-order shear deformation theory (TSDT). The effective properties of FG-CNTRC cylindrical pressure vessels are computed for different patterns of reinforcement, according to the rule of mixture. The governing equations of the problem are derived from the principle of virtual works and are solved as a classical eigenproblem under the assumption of clamped supported boundary conditions. A large parametric investigation aims at showing the influence of some meaningful parameters on the thermo-elastic response, such as the type of pattern, the volume fraction of CNTs, and the Pasternak coefficients related to the elastic foundation.


Introduction
Pressure vessels are one of the most important and expensive components for mechanical and structural engineering applications. Due to their importance, some international companies have published useful guidelines for design purposes, e.g., the ASME Standards [1]. At the present stage, however, the available standards do not deliver any suggestion for the design of reinforced pressure vessels. To this end, the structural behavior of reinforced pressure vessels has gained increasing attention among the scientific community, as detailed in what follows.
A preliminary work was proposed by El Mir et al. [2] to analyze cylindrical sandwich structures made of a weak orthotropic core under a certain patch loading. The authors presented a novel numerical formulation based on a high-order shear deformation theory (HSDT) to predict the structural response in static conditions. Arefi and Rahimi [3] and Arefi et al. [4][5][6] studied the thermo-elastic behavior of non-homogeneous functionally graded (FG) cylinders with clamped supports under a double mechanical and thermal loading condition, while applying the first-order shear deformation theory (FSDT). Menasria et al. [7] employed an innovative displacement field with unknown integral terms to analyze the thermal buckling behavior of a FG sandwich plate under a uniform thermal loading. They derived the governing equations of the problem, based on a variational principle. A further work by Dong et al. [8] studied the local buckling behavior of composite plates resting on a Wrinkler foundation under a uniform in-plane shear loading, while focusing on the effect of the ply angle and stiffness foundation on the critical buckling coefficients. Ahmad Bidi et al. [9] studied a reinforced curved steel structure with nanopolyurea, under a low velocity impact loading. A double experimental and numerical study was performed by the authors to explain the effect of the in-plain curvature and the nano-particle effect on the structural response. In a recent work by Rahmani et al. [10], the Hamilton's principle and the Navier solution were applied to examine different non-local HSDTs for FG nanobeams, including the size-dependent effects. Golmakani et al. [11], in addition, studied the non-linear thermo-elastic bending behavior of functionally graded carbon nanotube-reinforced composite (FG-CNTRC) plates resting on an elastic foundation. Nasihatgozar et al. [12] used a HSDT to study the free vibration response of thick-layer doubly curved sandwich panels resting on simply supports. This investigation delivered an optimum range for the core to face sheet stiffness ratio by considering the effective stress components, for the first time.
Arefi et al. [13] investigated the thermo-piezo-magneto-elastic behavior of a FG piezomagnetic cylinder under thermal, magnetic and mechanical loading, which was, in turn, FG in the radial direction. Mohammadimehr et al. [14] considered a viscoelastic piezoelectric polymeric nano-composite plate reinforced by FG single-walled carbon nanotubes (SWCNT) based on a modified strain gradient theory. The authors, in their work, studied the effect of the thickness-to-width ratio, as well as of the magnetic field, the applied voltage, the static loading factor, the viscoelastic parameters of the foundation, and the surface density constant, on the dynamic stability region. In addition, Hao et al. [15] focused on the thermo-mechanical stresses within the graded interphase between the fiber and the matrix in fiber-reinforced composites. Shen et al. [16] presented a non-linear bending analysis for nanocomposite cylindrical panels subjected to a transverse uniform or sinusoidal load resting on an elastic foundation in a thermal environment. They studied the effect of the CNT volume fraction distribution, the foundation stiffness, the temperature rise and the in-plane boundary conditions on the mechanical response of the structure. Moreover, Arefi and Rahimi [17] applied the FSDT to study the electro-elastic response of a FG piezoelectric cylinder under internal pressure.
Alibeigloo [18] investigated the free vibration behavior of FG-CNTRC cylindrical panels embedded in piezoelectric layers with simply supported boundary conditions. A coupled thermo-mechanical problem was also treated by Shen and Xiang [19], who considered the nanocomposite cylindrical shells reinforced by SWCNTs in thermal environments and studied the sensitivity of the response in terms of natural frequencies as well as of non-linear-to-linear frequency ratios for varying temperatures and CNT volume fractions. In line with the previous work, Arefi et al. [6] presented a two-dimensional thermo-elastic analysis of a FG thick-walled cylinder under thermal and mechanical loading resting on a Pasternak foundation. They applied the FSDT to describe the displacement field and found that an increased non-homogeneous index yields to a decrease of the radial and axial displacement. A finite element method (FEM)-based numerical analysis was applied by Chavan and Lal [20] to analyze the static behavior of SWCNT-FG plates. A large parametric investigation was aimed at investigating the effect of the width-to-thickness ratio, the stress distribution profile along the thickness, the boundary conditions, as well as the volume fractions on the overall structural behavior.
A further coupled thermo-piezoelectric problem was proposed by Khoshgoftar et al. [21], for a thick-walled cylinder with FG material, subjected to an inner and outer pressure under a constant temperature gradient. An interesting post-buckling coupled problem can be found in Kiani [22], who studied the thermo-mechanical behavior of CNTRC plates under a uniform increase of the thermal loading, and determined the effect of the reinforcement pattern on the critical temperature and the maximum post-buckling deflection for different distributions. Zhang et al. [23] investigated the free vibration response and flexural strength of CNTRC cylindrical panels with four different reinforcement distributions. Moreover, Rahimi et al. [24] focused on the electro-elastic behavior of a FG piezoelectric cylinder adopted as a physical sensor, and subjected to an internal pressure. A further work by Asadi et al. [25] aimed at investigating the aero-thermo-elastic behavior of supersonic FG-CNTRC flat panels in a thermal environment. To this end, a dynamic model was developed by the authors, according to the FSDT, whereas the presence of an aerodynamic pressure was found to play a key role not only for the onset of aero-thermal buckling instability, but also for the mode shapes of the composite structure.
In recent years, increased attention has been devoted by the scientific community to a comparative evaluation of FSDTs and HSDTs in a combined form with the novel generalized differential quadrature (GDQ) method for the study of the statics and dynamics of composite plates and shells of arbitrary shapes reinforced by agglomerated nanoparticles (see [26][27][28][29][30][31][32][33][34][35] among others). Despite the available works in the literature on the reinforced cylindrical pressure vessels, there is a general lack of analytical formulations, based on the third-order shear deformation theory (TSDT), to study FG-CNTRC cylindrical pressure vessels, in thermal conditions. In this context, the main task of the present work regards the novel application of the TSDT to study the thermo-mechanical response of a two-dimensional FG-CNTRC cylindrical pressure vessel resting on a Pasternak foundation. This represents an extended version of the FSDT-based formulation proposed in Arefi et al. [6] for the same problem, whose governing equations were tackled as a classical eigenproblem, for a clamped structure at the two extremities. In this work, a parametric analysis aims at investigating the sensitivity of the thermo-mechanical response of the composite structure to some important parameters such as the reinforcement pattern, the volume fraction of CNTs, or the Pasternak parameters related to the foundation. This represents a useful aspect to account for practical industry applications and optimization design. The paper is organized as follows. First, we present the problem formulation in Section 2, whose solution procedure is briefly reviewed in Section 3. Next, the parametric investigation and the main results are presented and discussed in Section 4. Finally, conclusions are drawn in Section 5.

Geometry and Material Properties
Let us consider a FG cylindrical shell with geometrical parameters defined in Figure 1. It is possible to assume different CNT distributions along the thickness, whose geometries are depicted in Figure 1, with the analytical expressions reported in Table 1 as a function of the effective volume fraction V * CNT , in which R, r and h indicate the average radius, arbitrary radius and the thickness of the cylinder which are clarified in Figure 1e. This last one is computed as follows [36][37][38][39] where W CNT is the mass fraction of the CNTs, while ρ CNT and ρ m refer to the density of the CNTs and matrix, respectively. In Table 1, all of the patterns are defined in terms of the same total volume fraction V * CNT . In a FG-X distribution of CNTs, the top and bottom surfaces reach the maximum values, whereas the mid-plane is completely free of CNTs, with a linear increase of CNTs along the thickness. The contrary occurs for a FG-O distribution of CNTs, where the external surfaces are free of CNTs, the mid-surface of the plate is enriched with CNTs and a linear decrease of CNTs is established from the mid-plane to the external surfaces. In a FG-V distribution, the top surface is enriched of CNTs, the bottom one is completely free of CNTs, with a linear increase in the amount of CNTs from the bottom to the top surfaces. In a UD type, the volume fraction of CNTs maintains constant along the thickness. Figure 1 illustrates the whole distributions of V CNT as a function of the thickness coordinate. Table 1. Volume fraction of CNTs as a function of thickness coordinate for various cases of CNTs distribution [22]. The rule of mixture is here applied to define the effective material properties, namely, the Young's moduli , , the shear modulus , the density , and the Poisson's ratio , as follows [40,41]:

CNTs Distribution
where , and refer to the efficiency parameters explaining the scale-dependent material properties, and are the volume fractions of the CNT and the matrix, respectively, which are related to each other as: The other effective mechanical properties , , , and , are described here below:  The rule of mixture is here applied to define the effective material properties, namely, the Young's moduli E 11 , E 22 , the shear modulus G 12 , the density ρ, and the Poisson's ratio ν 12 , as follows [40,41]:

CNTs Distribution
where η 1 , η 2 and η 3 refer to the efficiency parameters explaining the scale-dependent material properties, V CNT and V m are the volume fractions of the CNT and the matrix, respectively, which are related to each other as: The other effective mechanical properties E 33 , G 13 , ν 31 , and ν 32 , are described here below:

Basic Equations
According to the TSDT, the displacement field (see Figure 2) of the cylinder reads [42]: where u x and u z are the axial and radial displacement components, respectively, and u i {i = 1, 2, 3} and w i {i = 1, 2, 3} are functions of x. Thus, the strain field can be obtained by derivation, together with the axial component ε x , radial component ε z , circumferential component ε t and shear component The stress-strain relations read as follows: where, In Equation (7), T represents the increment in temperature from the reference state, and α 11 , α 22 and α 33 are the coefficients of thermal expansion.
On the other hand, the following constitutive relations must be considered, relating forces and displacements, due to the presence of the Pasternak foundation: k 1 and k 2 being the Pasternak coefficients. In Equation (9), the two kinematic terms can be expressed in the following expanded form: Substitution of Equation (10) into Equation (9) yields the following expression: Therefore, we can define the variation of the strain energy as follows: By the substitution of Equation (6) into Equation (12), Equation (12) can be rewritten as: For simplification purposes, Equation (13) can be rewritten in terms of internal stress resultants as follows: where the internal stress resultants can be introduced as follows: In addition, we can define the variation of the external work as follows: where P i is the internal pressure, here set to P i = 70 MPa. By partial integration of Equation (14) and by substitution of Equations (6) and (11) into Equation (16), get the following variations of the strain energy and energy of the external work: The application of the Hamilton's principle gives the following governing equations of the problem: and the substitution of Equations (17) and (18) into Equation (19) gives the following differential equations in terms of internal stress resultants: By substitution of Equations (15) and (11) into Equation (20), we obtain the following set of equations in the matrix form: where the displacement vector X, the coefficient matrixes (G 1 , G 2 and G 3 ) and the force vector F are defined as follows:

Solution Procedure
The problem above is solved including both homogeneous and particular solutions. Here below the complete solving procedure, that embraces the following 4 steps Step 1. The particular solution of Equation (18) is obtained as [43][44][45]: Step 2. The eigenvalues { = 1, 2, 3} are obtained by solving the characteristic equation: Step 3. The associated eigenvectors { = 1, 2, 3}are determined as: Step 4. The displacement components are computed by substitution of the eigenvalues and eigenvectors into the following equation: where refers to the unknown coefficients determined by enforcing the appropriate boundary conditions, here assumed as clamped-clamped at each side, namely:

Validation
Khoshgoftar et al. [21,45] studied the thermoelastic behaviour of a FG piezoelectric cylindrical structure, whose results have been here selected for validation purposes. In detail, the accuracy of the solution procedure is evaluated for = 0, where is the non-homogeneous parameter. A comparative evaluation of the results is shown in Figure 3 in terms of distribution of the electrical potential along the thickness for a thick walled cylinder made of piezoelectric material. The perfect

Solution Procedure
The problem above is solved including both homogeneous and particular solutions. Here below the complete solving procedure, that embraces the following 4 steps Step 1. The particular solution of Equation (18) is obtained as [43][44][45]: Step 2. The eigenvalues m i {i = 1, 2, 3} are obtained by solving the characteristic equation: Step 3. The associated eigenvectors v i {i = 1, 2, 3} are determined as: Step 4. The displacement components are computed by substitution of the eigenvalues m i and eigenvectors v i into the following equation: where c i refers to the unknown coefficients determined by enforcing the appropriate boundary conditions, here assumed as clamped-clamped at each side, namely:

Validation
Khoshgoftar et al. [21,45] studied the thermoelastic behaviour of a FG piezoelectric cylindrical structure, whose results have been here selected for validation purposes. In detail, the accuracy of the solution procedure is evaluated for l = 0, where l is the non-homogeneous parameter. A comparative evaluation of the results is shown in Figure 3 in terms of distribution of the electrical potential along the thickness for a thick walled cylinder made of piezoelectric material. The perfect agreement between our results and predictions from the literature confirms the reliability and accuracy of the proposed formulation.

Numerical results
We illustrate the use of the proposed formulation, here applied for a cylinder with inner radius = 0.04 m , outer radius = 0.05 m and length = 0.5 m . The structure is made of a polyethylmethacrylate (PMMA) material with the same mechanical properties for the matrix as found in Kiani [22], i.e., = 2.5 , = 45 × 10 , = 0.34. A SWCNT is here considered as fiber reinforcement, with geometrical, mechanical and thermal properties, as detailed in Table 2. Table 2. Thermo-mechanical properties of (10, 10) armchair single-walled carbon nanotube (SWCNT) (tube length = 9.26 nm, tube mean radius = 0.68 nm, tube thickness = 0.067 nm) [22].  Additionally, we consider different CNT distributions and different Pasternak coefficients for comparative purposes, while applying both the FSDT and TSDT to investigate the kinematic and static response of the reinforced structure.

Kinematic Response
First, we analyze the effect of the reinforcement distribution on the kinematic response of the structure. Figure 4 shows the longitudinal distribution of the axial displacement ( Figure 4a) and radial displacement (Figure 4b) at the middle surface ( = 0), for a UD, FG-X and FG-V pattern, as well as for a fixed volume fraction * = 0.12 and a null value of the Pasternak coefficients = = 0.
Based on the numerical results, the boundary conditions are clearly satisfied at two ending supports of the cylindrical pressure vessel. In addition, it is worth noting that the maximum radial

Numerical Results
We illustrate the use of the proposed formulation, here applied for a cylinder with inner radius r i = 0.04 m, outer radius r o = 0.05 m and length L = 0.5 m. The structure is made of a polyethylmethacrylate (PMMA) material with the same mechanical properties for the matrix as found in Kiani [22], i.e., E m = 2.5 GPa, α m = 45 × 10 −6 1 K , ν m = 0.34. A SWCNT is here considered as fiber reinforcement, with geometrical, mechanical and thermal properties, as detailed in Table 2. Table 2. Thermo-mechanical properties of (10, 10) armchair single-walled carbon nanotube (SWCNT) (tube length = 9.26 nm, tube mean radius = 0.68 nm, tube thickness = 0.067 nm) [22]. As a comparison, we assume three different volume fractions of CNTs, as follows [22]: Additionally, we consider different CNT distributions and different Pasternak coefficients for comparative purposes, while applying both the FSDT and TSDT to investigate the kinematic and static response of the reinforced structure.

Kinematic Response
First, we analyze the effect of the reinforcement distribution on the kinematic response of the structure. Figure 4 shows the longitudinal distribution of the axial displacement u x (Figure 4a) and radial displacement u z (Figure 4b) at the middle surface (z = 0), for a UD, FG-X and FG-V pattern, as well as for a fixed volume fraction V * CNT = 0.12 and a null value of the Pasternak coefficients k 1 = k 2 = 0. radial displacement reduces slightly, and attains the maximum values at the mid length of the structure (see Figure 6b).
The efficiency of the proposed TSDT formulation is also verified against the FSDT, through a comparative evaluation of the kinematic results in the axial and radial direction. As expected, more accurate results can be obtained by applying a TSDT compared to the other based on lower order theories [33], as clearly shown in Figure 7. In more detail, the TSDT-based axial displacement in the longitudinal direction is lower than the FSDT estimate (see Figure 7a), whereas the TSDT-based radial displacement is always higher than the FSDT prediction along the whole specimen (see Figure  7b).   Based on the numerical results, the boundary conditions are clearly satisfied at two ending supports of the cylindrical pressure vessel. In addition, it is worth noting that the maximum radial displacement is reached under a FG-X distribution, whereas the maximum axial displacement is obtained for a FG-V distribution.
A further check is evaluated on the sensitivity of the structural response to the volume fraction of reinforcement. Figure 5 illustrates the variation of the axial and radial displacements at the middle surface (i.e., for z = 0) along the longitudinal direction, by applying different volume fractions of FG-CNTRC, and considering a UD reinforcement pattern. Also in this case, we keep the Pasternak coefficients k 1 , k 2 equal to zero. A significant decrease in both axial and radial displacements is noticed for increasing volume fractions of CNTs, in line with findings by Arefi et al. [6]. An increased volume fraction of CNTs, indeed, yields to an increased structural stiffness, with a consistent reduction of the kinematic quantities. The variations of the longitudinal and radial displacements along the radial direction are listed in Tables 3 and 4 at various longitudinal positions. It is worth noticing that the maximum radial displacements are reached at the middle surface of the cylindrical shell. In addition, the maximum longitudinal displacements occur near one quarter of the cylindrical shell.
Due to the application of the TSDT, the displacement components are expected to vary along the radial direction, as represented in Figure 6 for a volume fraction V * CNT = 0.12, a uniform distribution of CNTs, and null values of the Pasternak coefficients k 1 = k 2 = 0. As visible in Figure 6a, the axial displacement assumes an asymmetric behavior at each surface along the longitudinal direction, while reaching the minimum value at the mid-surface (i.e., for z = 0), and the maximum values at the external lateral surfaces (i.e., for z = −h 2 and z = h 2 ). In addition, moving from the inner surface (z = −h 2 ) to the outer one (z = h 2 ) of the cylinder, the radial displacement reduces slightly, and attains the maximum values at the mid length of the structure (see Figure 6b).
The efficiency of the proposed TSDT formulation is also verified against the FSDT, through a comparative evaluation of the kinematic results in the axial and radial direction. As expected, more accurate results can be obtained by applying a TSDT compared to the other based on lower order theories [33], as clearly shown in Figure 7. In more detail, the TSDT-based axial displacement in the longitudinal direction is lower than the FSDT estimate (see Figure 7a), whereas the TSDT-based radial displacement is always higher than the FSDT prediction along the whole specimen (see Figure 7b).

Tensional Response
A similar parametric study is also repeated from a tensional point of view, in terms of axial, circumferential, radial and shear stresses along the longitudinal direction, and for different reinforcement patterns. Based on the numerical results in Figure 8, a FG-O pattern of CNTs seems to

Tensional Response
A similar parametric study is also repeated from a tensional point of view, in terms of axial, circumferential, radial and shear stresses along the longitudinal direction, and for different reinforcement patterns. Based on the numerical results in Figure 8, a FG-O pattern of CNTs seems to yield to the maximum values for the axial and radial stresses and to the minimum values for the circumferential and shear components except for the boundary zones. The stress response given by a UD and a FG-V distribution, instead, is almost similar for each component along the whole structure, with the highest numerical differences nearby the two clamped sides. The stress response can be affected significantly by the volume fraction of the reinforcement, whose variation is plotted in Figure 9, assuming volume fractions V * CNT = 0.12, 0.17, 0.28. This parametric investigation is here tackled for a UD pattern, but could be similarly repeated for all the other reinforcement distributions. Based on the numerical results, a general increase in magnitude is observed for the axial, circumferential, and radial stress components, and for a volume fraction V * CNT higher than 0.12 (see Figure 9a-c).   Only the shear stress seems to be almost insensitive to the volume fraction (see Figure 9d), in agreement with findings by Arefi et al. [6].
In addition, Figure 10 represents the longitudinal distribution of the stress components for the middle and external surfaces of the structure (i.e., for z = −h 2 , z = 0, z = h 2 , respectively). It is worth observing that the TSDT allows all stress components to assume different magnitudes through the thickness direction, as can be seen in Figure 10a  The last parametric study focuses on the sensitivity of the stress response to the Pasternak parameters of the elastic foundation. Figure 11 plots the main curves for different combinations of k 1 and k 2 , while keeping V * CNT = 0.12, FG − UD, z = 0. The numerical results in Figure 11 clearly state that lower magnitudes can be obtained for each stress component by increasing both Pasternak parameters. This is strictly related to a general increase in stiffness of the foundation for increasing values of k 1 and k 2 , with an expected sensitive reduction of the kinematic quantities. Nanomaterials 2019, 9,

Radial Distribution of Results
In this section, the radial distribution of numerical results are provided. For this purpose, the radial and axial deformation of the cylindrical shell are plotted along the radial direction at a defined longitudinal position. Figure 12a,b show the radial distribution of the longitudinal deformation at x = 0.1 and 0.4, respectively. One can conclude that the maximum longitudinal deformations occur at the external surfaces while the minimum values occur at the inner surfaces. Figure 13a,b plot the radial distribution of the radial deformation at x = 0.25 and 0.5, respectively. One can conclude that maximum radial deformations occur at the inner surfaces while the minimum deformations involve the outer surfaces.

Radial Distribution of Results
In this section, the radial distribution of numerical results are provided. For this purpose, the radial and axial deformation of the cylindrical shell are plotted along the radial direction at a defined longitudinal position. Figure 12a,b show the radial distribution of the longitudinal deformation at x = 0.1 and 0.4, respectively. One can conclude that the maximum longitudinal deformations occur at the external surfaces while the minimum values occur at the inner surfaces. Figure 13a,b plot the radial distribution of the radial deformation at x = 0.25 and 0.5, respectively. One can conclude that maximum radial deformations occur at the inner surfaces while the minimum deformations involve the outer surfaces.

Conclusions
The third-order shear deformation theory (TSDT) is employed in this work to study the thermo-elastic response of FG-CNTRC cylindrical pressure vessels resting on a Pasternak foundation. The rule of mixture is adopted herein for computing the effective material mechanical properties, whereby the governing equations of the problem are derived from Hamilton's principle and solved as a classical eigenvalue problem. A parametric investigation aims at investigating the sensitivity of the static and kinematic response to some important parameters such as the reinforcement distribution, its volume fraction, or the Pasternak parameters related to the elastic foundation. Based on the numerical results, the main conclusions can be pointed out as: • The accuracy of higher order theories, such as the TSDT, and lower order theories, such as the FSDT, must be determined comparatively with respect to the experimental results, but it is expected that a TSDT provides a more accurate structural response due to its capability to

Conclusions
The third-order shear deformation theory (TSDT) is employed in this work to study the thermo-elastic response of FG-CNTRC cylindrical pressure vessels resting on a Pasternak foundation. The rule of mixture is adopted herein for computing the effective material mechanical properties, whereby the governing equations of the problem are derived from Hamilton's principle and solved as a classical eigenvalue problem. A parametric investigation aims at investigating the sensitivity of the static and kinematic response to some important parameters such as the reinforcement distribution, its volume fraction, or the Pasternak parameters related to the elastic foundation. Based on the numerical results, the main conclusions can be pointed out as: • The accuracy of higher order theories, such as the TSDT, and lower order theories, such as the FSDT, must be determined comparatively with respect to the experimental results, but it is expected that a TSDT provides a more accurate structural response due to its capability to

Conclusions
The third-order shear deformation theory (TSDT) is employed in this work to study the thermo-elastic response of FG-CNTRC cylindrical pressure vessels resting on a Pasternak foundation. The rule of mixture is adopted herein for computing the effective material mechanical properties, whereby the governing equations of the problem are derived from Hamilton's principle and solved as a classical eigenvalue problem. A parametric investigation aims at investigating the sensitivity of the static and kinematic response to some important parameters such as the reinforcement distribution, its volume fraction, or the Pasternak parameters related to the elastic foundation. Based on the numerical results, the main conclusions can be pointed out as:

•
The accuracy of higher order theories, such as the TSDT, and lower order theories, such as the FSDT, must be determined comparatively with respect to the experimental results, but it is expected that a TSDT provides a more accurate structural response due to its capability