Mathematical Modeling and Analytical Solution of Thermoelastic Stability Problem of Functionally Graded Nanocomposite Cylinders within Different Theories

: Revolutionary advances in technology have led to the use of functionally graded nanocomposite structural elements that operate at high temperatures and whose properties depend on position, such as cylindrical shells designed as load-bearing elements. These advances in technology require new mathematical modeling and updated numerical calculations to be performed using improved theories at design time to reliably apply such elements. The main goal of this study is to model, mathematically and within an analytical solution, the thermoelastic stability problem of composite cylinders reinforced by carbon nanotubes (CNTs) under a uniform thermal loading within the shear deformation theory (ST). The inﬂuence of transverse shear deformations is considered when form-ing the fundamental relations of CNT-patterned cylindrical shells and the basic partial differential equations (PDEs) are derived within the modiﬁed Donnell-type shell theory. The PDEs are solved by the Galerkin method, and the formula is found for the eigenvalue (critical temperature) of the functionally graded nanocomposite cylindrical shells. The inﬂuences of CNT patterns, volume fraction, and geometric parameters on the critical temperature within the ST are estimated by comparing the results within classical theory (CT).


Introduction
For the development of almost all sectors of the economy, materials with improved performance properties, high strength and light weight, and that are resistant to abrasion, pressure, and various types of radiation are needed. For example, in shipbuilding and civil and mechanical engineering, lightweight and wear-resistant materials are in demand, which are necessary for the creation of structural cladding elements. Developments continue for the aerospace and automotive industries, aiming to reduce the weight and fuel consumption of parts while maintaining the strength of the body and its parts. The nuclear industry needs new materials for power reactors that can withstand high temperatures, pressures, and radiation. In modern technology, the above problems can be solved by using carbon nanotubes or nanocomposites in which they are used as reinforcements. Carbon nanotubes (CNTs) are carbon allotropes with a cylindrical nanostructure and are classified into single-walled and multi-walled carbon nanotubes. Carbon nanotubes, which have unique electrical, thermal, and mechanical properties with high specific surface areas and aspect ratios, were discovered by Sumio Ijima in 1993 [1]. Depending on their chirality, both their metallic and semiconductor properties made it possible to use CNTs in the field of microelectronics and power engineering, as well as a reinforcing element in the formation of nanocomposites [2][3][4].
Nanocomposites are generally divided into three categories depending on the type of main matrix: ceramic, metal, and polymer-based nanocomposites. Among these nanocomposite species, polymer nanocomposites are recently in demand. Polymer-based nanocomposites differ from conventional polymer composite materials in lighter weight as well as higher impact and abrasion resistance as well as good resistance to chemical attack, making them suitable for use in military and aerospace developments. The automotive industry is a leader in the growth of the use of nanocomposite materials. Polymer nanocomposites have become the main composite materials for fuel lines, replacing conventional steel to prevent static build-up [5][6][7][8][9][10][11][12]. Today, new results are revealed by performing molecular dynamics (MD) simulations and experimental tests on CNT-reinforced polymers. For example, in the Duan [13] study, based on the hypothesis that only the tube walls in direct contact with the polymer matrix carry the external load, the concept of a true volume ratio for CNT was proposed and a new micromechanical model on the modulus of nanocomposites was developed based on this concept. All these demonstrate that there are still efforts to develop theories about nanocomposites and that interest in this field continues to increase.
With superior performance in very high temperature environments, nanocomposites are among the best new generation composites that can be used as thermal barriers in many advanced technology areas, typically in nuclear reactors, spacecraft, and aerospace. As cylindrical shell structures are one of the main structural elements made up of nanocomposites, perhaps the first, there is a further need to investigate their buckling behavior in the thermal environment and under thermal loading. After the first attempt by Shen [14] on thermal buckling of nanocomposite cylindrical shells, some research has been done on this topic . Most of these studies are related to the buckling problem of FG-CNTreinforced structural elements in the thermal environment or thermal loading, and their solutions have been realized by using various theories and numerical methods.
The literature review revealed that the investigation of the thermal stability behavior of cylindrical shells based on CNTs is insufficient. In this study, the analytical solution of the thermal stability problem of composite cylinders reinforced by CNTs under the uniform thermal loading is presented within different theories. This study will be beneficial in terms of updating the buckling analysis of cylindrical shells made of carbon nanotubereinforced composites in the thermal environment and developing calculation methods during the design. As one of the main structural elements consisting of nanocomposites is the cylindrical shell, examining the buckling behavior in the thermal environment and under thermal loading will help to eliminate the problems that may arise in industries where advanced technologies are used during the design phase. The analytical solution of the buckling problem of cylindrical shells made of nanocomposites within the framework of various theories in the thermal environment and obtaining original formulas is also an important achievement from a mathematical point of view.

Theoretical Development
Consider the CNT-reinforced composite cylindrical shell of radius a, length l, and thickness h under the uniform thermal loading. The Oxyz coordinate system is located at the left end of the mid-surface of the cylinder, and the x, y, z axes are oriented, as in

Theoretical Development
Consider the CNT-reinforced composite cylindrical shell of radius a , length l , and thickness h under the uniform thermal loading. The Oxyz coordinate system is located at the left end of the mid-surface of the cylinder, and the , x y , z axes are oriented, as in Figure 1. Three common approaches are used for the modeling of carbon nanotube-reinforced composites: (1) atomistic modeling, e.g., molecular dynamics (MD) simulations; (2) continuum mechanics modeling; and (3) multiscale modeling. Besides atomic modeling, continuum mechanics (or micromechanical models) has been extended to modeling of the mechanical properties of nanocomposites. Among the micromechanical models, the most popular are the rule of mixture (RM), the Mori-Tanaka (MT) model, and the representative volume element (RVE) model. In this study, the effective material properties, such as Young's modulus and Poisson's ratio of nanocomposite cylindrical shells, are obtained using the RM method and are expressed as follows [14,17,35 where V denotes the volume fraction, superscript "cnt" and " m " represent the corresponding property for CNT and matrix, respectively, the efficiency parameters are indicated by where Figure 1. CNT-reinforced composite cylindrical shell and coordinate system. Three common approaches are used for the modeling of carbon nanotube-reinforced composites: (1) atomistic modeling, e.g., molecular dynamics (MD) simulations; (2) continuum mechanics modeling; and (3) multiscale modeling. Besides atomic modeling, continuum mechanics (or micromechanical models) has been extended to modeling of the mechanical properties of nanocomposites. Among the micromechanical models, the most popular are the rule of mixture (RM), the Mori-Tanaka (MT) model, and the representative volume element (RVE) model. In this study, the effective material properties, such as Young's modulus and Poisson's ratio of nanocomposite cylindrical shells, are obtained using the RM method and are expressed as follows [14,17,35]: where V denotes the volume fraction, superscript "cnt" and "m" represent the corresponding property for CNT and matrix, respectively, the efficiency parameters are indicated by η i (i = 1, 2, 3), and the Young and shear moduli of the polymer (or matrix) and CNTs are indicated by Y m , Y cnt kk (k = 1, 2) and G cnt 12 , respectively. Here ν 12 and ρ are Poisson's ratio and density, respectively.
in which m cnt is a mass fraction of CNTs. The use of these relations is based on the following hypotheses: CNTs are reinforced into the polymer matrix as linear functions and CNTs reinforced into the polymer matrix carry the charge transferred to the matrix.
While the first line of the Equation (2) shows that the CNT distribution is uniform (U), the other distributions are in the form of linear functions and are represented by the symbols V, O, and X, respectively. Based on Equation (2), the patterns formed by the CNT distribution throughout the thickness of the polymer cylindrical shell are illustrated in Figure 2.
into the polymer matrix as linear functions and CNTs reinforced into the polymer matrix carry the charge transferred to the matrix.
While the first line of the Equation (2) shows that the CNT distribution is uniform (U), the other distributions are in the form of linear functions and are represented by the symbols V, O, and X, respectively. Based on Equation (2), the patterns formed by the CNT distribution throughout the thickness of the polymer cylindrical shell are illustrated in Figure 2. In this study, it is assumed that the thermal expansion coefficients, as well as the material properties, vary with temperature and location [14,16,17,34 where 11 22 , cnt cnt α α , and m α are thermal expansion coefficients of CNT and matrix, respectively.

Governing Equations
If we generalize the ST to the heterogenous nanocomposite cylindrical shells, the displacement components can be represented as follows [34,37- where u and v are the displacements of axial and circumferential directions on the reference surface, respectively, w is the deflection, and 1 ( , ) x y ϕ and 2 ( , ) x y ϕ are the transverse normal rotations about the y and x axes, respectively.
Considering the assumptions of the ST [39], the following relations can be written: The basic relationships of CNT-originated cylindrical shells under the temperature field are constructed in the scope of ST as follows [14,17,34,35,40]: In this study, it is assumed that the thermal expansion coefficients, as well as the material properties, vary with temperature and location [14,16,17,34]: where α cnt 11 , α cnt 22 , and α m are thermal expansion coefficients of CNT and matrix, respectively.

Governing Equations
If we generalize the ST to the heterogenous nanocomposite cylindrical shells, the displacement components can be represented as follows [34,[37][38][39] where u and v are the displacements of axial and circumferential directions on the reference surface, respectively, w is the deflection, and ϕ 1 (x, y) and ϕ 2 (x, y) are the transverse normal rotations about the y and x axes, respectively. Considering the assumptions of the ST [39], the following relations can be written: where e 33 is the normal strain in the direction z, σ i3 (i = 1, 2) are the transverse shear stresses, and f i (z), (i = 1, 2) are the shape functions of transverse shear deformations (TSDs) of the CNT-originated cylindrical shells [37][38][39].
To find the expression for the critical temperature, expressions (18) are substituted into the basic Equations (16) and after applying the Galerkin method, the determinant of the matrix of coefficients of algebraic equations is equal to zero: where To find the eigenvalue or the critical uniform temperature (CUT) for the CNT patterned cylinders within ST, Equations (15) and (19) are solved together and the following closed form solution is obtained: When TSDs are ignored in the basic relationships, the expression (24) transforms into the expression for CUT of the nanocomposite cylindrical shells within CT and it is shown as T cr cst . The critical uniform temperature of nanocomposite cylindrical shells within the ST is found by minimizing Equation (22) versus the (m, n).

Results
In this section, to confirm the accuracy of the analytical solution obtained in this study, first a comparison is made with an existing and reliable result, and then original analyses are performed.

Comparative Study
In this subsection, the results for the buckling load of cylindrical shells modeled with different types of CNTs in a thermal environment and under axial load are compared with solutions of ref. [41] within ST and given in Table 1. As Shen [41] used N cr shen = 2πaN cr sdt , we use the formula (22) by multiplying 2πa at N 0 11 = −N, where N is the uniform axial load. The nanocomposite cylindrical shell characteristics are: a/h = 30, h = 0.002 m, V * cnt = 0.17, l = 10 √ 3ah [41]. It is observed from Table 1 that current results for N cr shen = 2πaN cr sdt (in kN) of the cylinder with U and X profiles are in good agreement with the results of ref. [41] in the thermal environment.

Analysis and Interpretations
The thermomechanical properties of PMMA and CNT at room temperature, that is, for T = 300 K, are numerically obtained using the following temperature-dependent functions. [14,17]: The efficiency parameters are described by [41]: η 1 = 0.137, η 2 = 1.022, η 3 = 0.7η 2 at V * cnt = 0.12; η 1 = 0.142, η 2 = 1.626, η 3 = 0.7η 2 at V * cnt = 0.17; η 1 = 0.141, η 2 = 1.585, η 3 = 0.7η 2 at V * cnt = 0.28 The CNT sizes are as follows: l cnt = 9.26 nm, a cnt = 0.68 nm, h cnt = 0.067 nm, ν cnt 12 = 0.175, and f 1 (z) = f 2 (z) = z − 4 3h 2 z 3 [39][40][41]. Table 2 presents the variation of the critical uniform temperature when the CNT cylindrical shell parameters are l/a = 1.25, h/a = 0.05, h = 0.002 m. It was determined that the critical temperature values decrease significantly, while those change irregularly, as the volume fraction increases. In addition, although the critical temperature decreases significantly, the number of waves corresponding to the minimum critical temperature does not change. The highest influence of TSDs on CUT is 22.51%, occurring in the X-pattern at V * cnt = 0.28, while the smallest effect is 6.01%, occurring in the O-model at V * cnt = 0.17. V-, O-, and X-patterns exhibit different effects compared to the U-model. In the framework of ST, the largest pattern effect on the critical temperature occurs in the X-model (+24.39%) for V * cnt = 0.28, whereas the least influence occurs in the V-model (−10.16%) for V * cnt = 0.17. It is observed that the pattern effect on the CUT value in the framework of CT is more pronounced than ST, between 11% and 13% in the X-pattern, between 3% and 4% in the O-pattern, and between 3% and 5% in the V-pattern. In addition, among these patterns, only the critical temperature values in the X-pattern are greater than the values in the U-pattern, and therefore the (+) sign is placed in front of them. Figure 3 shows the changes of the critical uniform temperature depending on the a/h ratio within ST and CT when the CNT cylinder parameters are l/a = 1.25, h = 0.002 m, V * cnt = 0.28. It is observed that the magnitudes of the critical temperature continuously reduce with the increase in the a/h ratio. The influences of TSDs on the critical temperature decrease with the increasing of the ratio, a/h. The influences of TSDs on the critical temperature for the U-, V-, O-, and X-patterns are 14.42%, 9.13%, 8.94%, and 22.51%, respectively, as a/h = 20, while those effects are 2.04%, 0.72%, 0.78%, and 2.72%, respectively, as a/h = 50 (Figure 3). For some models, when a/h > 50, CT can be used instead of ST. When compared with the uniform pattern, the effects of the V-and O-patterns on the CUT decrease from (−11.75%) to (−4.17%) and from (−16.41%) to (−11.81%), respectively, with the increase in a/h from 20 to 50, while this effect changes irregularly in the X-model, around 24.4-26.6%. In addition, it is observed that the consideration of transverse shear deformations significantly reduces the effect of patterns on critical temperatures.     Figure 4 shows the variation of the CUT depending on the l/a ratio within different theories, as the CNT cylindrical shell characteristics are h/a = 0.05, h = 0.002 m, and V * cnt = 0.28. As can be observed, CUT values decrease continuously depending on the l/a ratio, and this decrease is evident. The fastest decrease is observed in the U-and X-patterned cylindrical shell in the framework of ST. The use of the shear deformation theory significantly reduces the critical temperature values compared to the CT. The shear deformations effect on the CUT demonstrates that it decreases with the increase in l/a. For example, at l/a = 0.75, the effects of TSDs on the CUT are 36.96%, 26.51%, 26.44%, and 49.3% for the U-, V-, O-, and X-patterns, respectively, while those effects are 2.02%, 1.44%, 0.94%, and 3.56%, respectively, as l/a = 2.5. Although CT can be used instead of ST for l/a > 2 ratios in some patterns, its use at high temperatures may lead to incorrect results. When compared with the uniform pattern, the influences of the V-and O-patterns on the CUT reduce from (−15.09%) to (0.588%) and (−16.4%) and (−7.06%), respectively, with the increase in l/a from 0.75 to 2.25, while this effect changes irregularly in the X-model, at around 16.26-24.4%. As can be observed from the calculations, using CT does not cause an error when l ≥ 2.5a in V-patterned cylindrical shells.

Conclusions
In this study, thermal buckling behaviors of nanocomposite cylindrical shells are discussed in the scope of the ST. The basic relations of CNT-patterned cylindrical shells are first constructed, and then the governing equations are derived within the framework of the modified Donnell-type shell theory, considering the effect of TSDs. By solving the basic PDEs with the Galerkin method, the closed-form solution for the CUT of nanocomposite cylindrical shells is found. The effects of CNT patterns, volume fraction, and cylinder parameters on the critical uniform temperature within the ST are estimated in comparison with the results in the framework of the CT.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
In Equation (16), the following definitions apply:

Conclusions
In this study, thermal buckling behaviors of nanocomposite cylindrical shells are discussed in the scope of the ST. The basic relations of CNT-patterned cylindrical shells are first constructed, and then the governing equations are derived within the framework of the modified Donnell-type shell theory, considering the effect of TSDs. By solving the basic PDEs with the Galerkin method, the closed-form solution for the CUT of nanocomposite cylindrical shells is found. The effects of CNT patterns, volume fraction, and cylinder parameters on the critical uniform temperature within the ST are estimated in comparison with the results in the framework of the CT.

Conflicts of Interest:
The authors declare no conflict of interest.