## 1 Introduction

In recent decades, the use of nanocomposite structures has been growing rapidly and has attracted much attention from scientific research communities [1, 2]. Among the utilized nanofillers, CNT and graphene have become very popular due to their unique thermo-mechanical properties. Specifically, values of up to 3000 W/m·K [3, 4] and 5300 W/m·K [5] for the thermal conductivities of CNT and graphene have been reported, respectively. Because of such extraordinary thermal behavior, these two nanofiller are excellent candidates to improve thermal conductivity of polymeric nanocomposite materials. On the other hand, it is well established that the FG dispersion of nanofiller in nanocomposite materials results in more manageable thermal-mechanical responses of nanocomposite structures [6, 7, 8, 9, 10].

Since thermal analysis is very critical in engineering structures, many researchers have investigated the thermal conductivity of nanocomposite materials. These materials and structures generate one of the hot topics of research in several emerging technologies and, therefore, they are being extensively investigated by researchers because of considerably improved their thermo-mechanical behaviors [11, 12]. Several parameters including nanofiller size, volume fraction, aspect ratio, thermal conductivity and intersection thermal resistance as well as temperature and nanofiller-matrix interface stickiness have influence on the effective thermal conductivity of nanocomposites. Highly conductive fillers, like metallic or carbon-based materials, can be added into nanocomposite cylinders to improve their thermal conductivity [13, 14]. Dispersion/ aggregation state and filler shape are the critical factors influencing the thermal conductivity of nanofluids and solid composites. Regarding nanofillers, the resistance of thermal interfacial is another critical factor affecting thermal transport. Researchers believe that the main factor limiting the thermal performance of suspensions and CNT composites is interfacial thermal resistance [15, 16, 17]. Gu *et al*. [18] performed some tests and proved that the use of 21.4% graphene platelets in a type of polyethylene (PE) called ultra-high molecular weight polyethylene (UHMWPE) increased the thermal conductivity of host material by up to around nine times the original UHMWPE matrix. Zabihi *et al*. [19] studied the effects of different types of defects on the effective thermal conductivity of nanocomposite materials reinforced with defected CNT, graphene and hybrid CNT/graphene fillers. Moheimani *et al*. [20] also conducted extensive research on the development of a closed-form unit cell micromechanical model for determining the effective thermal conductivities of unidirectional CNTs reinforced polymer nanocomposites.

Nanocomposite materials reinforced with CNT or graphene have been extensively utilized in different shell and plate structures in order to analyze their structural thermo-mechanical responses by many researchers in recent years. Qin *et al*. [21] proposed unified solution for laminated FG shallow shells reinforced by graphene with general boundary conditions; upon which free vibration of FG cylindrical shells and plates reinforced by CNTs and graphene with both classic and non-classic boundary conditions were investigated [22, 23]. Shokri-Oojghaz [24] studied the effect of reinforcing with aggregated CNT on the stress distributions of nanocomposite sandwich cylinders based on finite element method. Temperature effect on the static and dynamic responses of nanocomposite cylinders reinforced with wavy CNT were studied by Moradi-Dastjerdi *et al*. [25, 26, 27]. Safaei *et al*. [28] analyzed the free vibration of nanocomposite plates with CNTs embedded in amorphous polyethylene. They used the rules of mixture according to different plate model theories to find the fundamental frequencies of nanocomposite plates. The macro-mechanical characteristics of nanocomposites are influenced by CNT and graphene volume fraction and microstructure. Safaei *et al*. [29, 30] examined the thermal and mechanical behaviors of lightweight polymeric nanocomposite porous sandwich plates resting on two-parameter elastic foundation. Damadam *et al*. [31] found that the increase of non-dimensional temperature led to ratcheting in plastic strain in subsequent cycles. Furthermore, they introduced Bree’s Diagram to more effectively plot the material behavior of thick-walled cylinders under different loading conditions and were able to identify safe loading boundary for different temperatures and pressures. Jalali *et al*. [32] studied mechanical analysis of micro FG beams considering thermal conditions using a numerical method to see the effect of temperature change in natural frequencies of FG micro beams. The effects of thermal conductivity, temperature change, heat wave speed, heat flux, CNTs volume fraction and end support conditions on the nonlinear vibration of beam were discussed in detail by Pourasghar and Chen [33]. They also combined the applications of differential quadrature and Newton Raphson methods and used the developed method to solve the non-Fourier heat conduction equations to obtain temperature, displacement and stress in CNT reinforced, nanocomposite cylindrical panels [34]. Safaei *et al*. [35] investigated transient heat transfer and stress wave propagation in polymeric sandwich plates with two nanocomposite face sheets were made of clusters of CNTs embedded in a polymeric matrix. Zargar *et al*. [36] investigated thermal response of porous circular fins using numerical/analytical methods to see the effect of different geometry and nondimensional parameters in this problem. Moradi-Dastjerdi and Behdinan [37] studied the effect of graphene on the static and free vibration responses of nanocomposite cylinders subjected to thermal gradient load. Moreover, steady state and transient heat transfer responses of CNT-reinforced cylinders were performed in [38].

As mentioned in the above literature review, the effect of graphene or CNT on the thermal conductivity of polymeric nanocomposite materials has widely been investigated in literature.However, the impacts of these outstanding nanofillers on the heat transfer responses of axisymmetric cylinders have not fully been considered yet. Therefore, in this work, we have tried to address this shortcoming by the investigation and comparison of the heat transfer responses of FG axisymmetric cylinders reinforced with randomly oriented graphene or CNT. To do so, nanocomposite cylinders subjected to different thermal boundary conditions were considered. It should be noted that the thermal conductivities of nanocomposite materials were estimated through an accurate micromechanical model. Moreover, governing thermal equation was derived using a developed mesh-free method which included MLS shape function and transformation approach. Using the developed methods, the effects of nanofiller content and distribution as well as thermal boundary conditions on the heat transfer responses of FG nanocomposite cylinders were presented.

## 2 Modeling of nanocomposite cylinders

As mentioned earlier, randomly orientated graphene and CNT were selected as nanoscale reinforcements of nanocomposite cylinders with length *L*, inner radius *r _{i}*, and outer radius

*r*, as illustrated in Figure 1. The boundaries of nanocomposite cylinders were assumed to be subjected to a constant temperature

_{o}*T*, heat flux

*Q*or thermal convection with ambient temperature

*T*and convective heat transfer coefficient

_{a}*h*. All the thermal boundary conditions are also assumed to be uniform on their surfaces. It should be noted that boundary conditions and structure are axisymmetric, therefore, three-dimensional analysis of such structures can be reduced to a two-dimensional analysis without losing the accuracy of in results.

### 2.1 Dispersion of nanofillers

With the aim of optimizing the amount of utilized nanofillers to provide the best design, the concept of FGMs was employed in this work. Therefore, four linear FG profiles were selected for the dispersions of nanofillers along the radial direction of cylinders. In addition, a nanocomposite cylinder with uniform dispersion (UD) profile was considered to compare its steady state and transient thermal responses with FG nanocomposite cylinders.

According to the dispersions of nanofillers, the four FG types of dispersion were labeled as V, *∧*, X and O. The schematic FG distributions of nanofillers are shown in Figure 2. The variations of the volume fraction of nanofillers *f _{r}*(

*r*) was also described as the following equations [37]:

where *r _{m}* is the mid radius of cylinder. Moreover, the difference between the maximum

*ᐃf*.

_{r}### 2.2 Thermal properties

The thermal conductivity of polymeric nanocomposite materials reinforced with graphene or CNT can be estimated by employing a micromechanical model developed in [39] and modified in [19, 40]. This method was based on the volume fractions *f _{r}*, thermal conductivity

*k*and aspect ratio

^{r}*p*of nanofillers. According to this method, the thermal conductivities of nanocomposite materials reinforced with graphene or CNT was estimated as follow [19, 40]:

where

where *k ^{m}* is the thermal conductivity of polymeric matrix and

*l*and

*R*are nanofiller length and thermal resistance between the interfaces of nanofiller-nanofiller and nanofiller-polymeric, respectively. It should be noted that the estimation of

_{k}*k*through this approach was already verified by experimental results of graphene [18] and CNT-reinforced [41] polymers.

Furthermore, the mass density *ρ* and specific heat conduction *c _{p}* of the resulted nanocomposite materials were roughly calculated as follow:

where the superscripts of *r* and *m* indicate nanofiller and matrix, respectively, and matrix volume fraction is shown by *f _{m}* = 1 −

*f*.

_{r}## 3 Meshless formulations

The governing equation for transient heat conduction in an FG axisymmetric cylinder was stated as [38, 42]:

where *t* is time and *T*(*r*, *z*) is temperature at point **X** = **X*** _{i}*(

*r*,

*z*). The initial and boundary conditions were also defined as:

where *Г* shows boundaries including inner, outer, upper and lower surfaces of axisymmetric cylinders.

To analyze the thermal conduction of the considered nanocomposite cylinders, the weak form of Eq. (7) was treated through the standard form of Galerkin method. Developing a mesh-free method based on moving least square (MLS) shape functions **Φ**, the vector of generalized values of temperature field

where

where *n* is node number located inside the support domain of the defined shape function. Moreover, MLS shape function at **X** = **X*** _{i}*(

*r*,

*z*) was described as [43]:

where **P**, *w* and **M** are base vector, cubic Spline weight function and moment matrix, respectively. In axisymmetric condition, **P** and **M** were determined as:

The submission of approximated temperature field Eq. (9) in the weak form of governing thermal equation Eq. (7) was resulted in following discrete mesh-free form [38]:

where *M _{T}* = 0 and

*S*=

*Q*for boundaries subjected to thermal flux, and

*M*=

_{T}*h*and

*S*=

*hT*for boundaries under convection environment. By applying such boundary conditions, Eq. (13) was rearranged as:

_{a}where

in which

Before the solution of Eq. (14), essential boundary conditions had to be applied to this equation. However, such condition could not directly be imposed because Eq. (14) was based on the generalized values of temperature field. This extra process was due to the shortcoming of MLS shape functions in the satisfaction of Kronecker delta property. To overcome this problem, a transformation matrix **R** was employed to modify

where

The substitution of Eq. (17) in Eq. (14) gave:

where,

The essential boundary conditions be applied in Eq. (19) as easy as FEM.

## 4 Results and discussions

Steady state and transient heat transfer responses of the considered FG nanocomposite cylinders are presented in this section. Nanocomposite materials were made of a mixture of graphene or CNT as reinforcement part embedded in polyethylene as polymeric matrix. The properties of the materials employed in the simulations were [18, 19, 41]:

According to Eq. (2), Figure 3 shows the variations of the thermal conductivities of CNT and graphene-reinforced nanocomposite materials as a function of nanofiller volume fraction. This figure shows that by the increase of nanofiller content, the thermal conductivities of both nanocomposites are significantly improved. The comparison between CNT and graphene shows that, graphene/PE nanocomposite enjoys higher thermal conductivity than CNT/PE nanocomposite because of graphene shape and PE-graphene thermal resistance *R _{k}*.

### 4.1 Validation study

Here, we examine the accuracy of the proposed solution framework for steady state and transient heat transfer responses of axisymmetric cylinders. First, the steady state responses of an isotropic axisymmetric long cylinder with inner temperature *T _{i}* = 300 K and outer temperature

*T*= 350 Khas been investigated. Hetnarski and Eslami [42] presented the exact solution of such cylinders as the following equation:

_{o}Figure 4 compares the temperature profiles along cylinder thickness obtained from our solution with exact results. An excellent agreement between the results can be observed in this figure which shows the accuracy of the proposed solution method.

In addition to steady state response, the accuracy of the proposed method has been examined for transient response of graphene and CNT reinforced nanocomposite cylinders. In this case, we tracked the transient solutions of temperature profiles until their steady state conditions were achieved and compared the final profile of each cylinder with their corresponding temperature profiles obtained from steady state solution. Hence, we considered cylinders at initial temperature *T*^{0} = 300 K subjected to inner heat flux *Q _{i}* = 8 kW/m

^{2}and outer heat convection environment with

*h*= 200 W/m

_{o}^{2}K and

*T*

_{a}_{−o}= 250 K when upper and lower faces of cylinders are thermally insulated i.e.

*Q*= 0 and

_{u}*Q*= 0. In these cylinders, the volume fractions of nanofillers are varied from

_{d}*L*= 0.04 m,

*r*= 0.1 m and

_{o}*r*= 0.08 m. Figure 5 shows that the transient responses of both cylinders are gradually converged to steady state temperature profiles. After 2500 seconds, perfect consistency between transient response and steady state solutions are observed for both cylinders.

_{i}### 4.2 Heat transfer response

Once the developed mesh-free method was validated, the transient heat transfer of FG nanocomposite cylinders reinforced with graphene or CNT were investigated under a variety of boundary conditions. The effects of nanofiller distribution on temperature profile and the time required to reach steady state condition have been studied. Therefore, nanocomposite cylinders with *L* = 0.04 m, *r _{i}* = 0.08 m,

*T*

^{0}= 300 K were considered. The inner and upper surfaces of cylinders are assumed under the heat fluxes of

*Q*= 8 kW/

_{i}*m*

^{2}and

*Q*= 4 kW/

_{u}*m*

^{2}. In addition, the bottom and outer surfaces are in heat convection environments with

*h*= 100 W/(m

_{o}^{2}K),

*h*= 200 W/(m

_{d}^{2}K) and

*T*

_{a}_{−o}=

*T*

_{a}_{−d}= 250 K. Figure 6a shows temperature time histories of FG and UD-CNT reinforced nanocomposite cylinders at their mid points. Figure 6b also shows temperature distributions of cylinders at their steady state conditions at

*z*=

*L*/2. For FG and UD-graphene reinforced cylinders under the same conditions, Figures 7a and 7b also illustrate temperature time histories and steady state temperature profiles, respectively. These figures disclose that nanofiller distribution not only affects stationary times but also governs the obtained temperature profiles such that X cylinder has the shortest stationary time and O cylinder has the longest one. Furthermore, it shows that temperature at locations with higher concentration of nanofiller are closer to that of the nearest boundary. Comparison between the effects of graphene and CNT on thermal responses shows that graphene reinforced cylinders have lower temperature variation than those reinforced with CNT because they have higher thermal conductivity as shown in Figure 3. It should be mentioned that cylinders with higher thermal conductivity have a better thermal flow which results in cylinders with lower temperature variation.

Moreover, Figure 8 illustrates the two-dimensional temperature profiles of the considered graphene reinforced cylinders. The significant effect of graphene dispersion on temperature profile are clearly visible in this figure such that the highest temperature is observed in V-type cylinder and it has a smallest high spread at its outer lower part, while *∧*-type cylinder has the biggest high temperature spread at its inner upper part.

In order to investigate the effect of nanofiller volume fraction on thermal responses, X-CNT/PE and X-graphene/PE nanocomposite cylinders with the same dimensions, *T _{i}* = 300 K,

*T*= 350 K and upper and lower faces were thermally insulated. Figures 9 and 10 show the thermal responses of CNT and graphene reinforced cylinders, respectively, in terms of temperature time histories and steady state temperature profiles. The results disclose that the increase of nanofiller contents significantly reduces the stationary time of temperature profiles. It also has a remarkable impact on steady state temperature profiles such that the use of more nanofiller increases the temperature gradient of cylinders.Moreover, it was observed that the use of graphene has stronger effects on stationary time than the use of CNT such that graphene/PE, with higher thermal conductivity than CNT/PE (see Figure 3), has shorter approaching time to steady state condition.

_{o}Here, the effects of different thermal boundary conditions on the steady state temperature profiles of graphene/PE with the same dimensions, *Q _{i}* = 5 kW/

*m*

^{2},

*Q*= −

_{o}*Q*(

_{i}*r*/

_{i}*r*) and

_{o}*T*=

_{u}*T*= 300 K. Due to symmetric thermal conditions, a symmetric temperature profile was observed for UD and X cylinders although temperature gradient in X cylinder was much higher than UD cylinder. However, V cylinder experiences an asymmetric temperature profile such that the inner half thickness of cylinder has wide temperature gradient and temperature at inner radius goes above 344 K.

_{d}Considering an outer thermal convection environment for previous cylinders, Figure 12 illustrates temperature profiles when *h _{o}* = 150 W/(m

^{2}K) and

*T*

_{a}_{−o}= 250 K. The obtained temperature profiles show that the effects of such heat convection environment are higher than employed internal heat flux because coldness (

*T*<300 K) spread is more visible in all figures.

Figure 13 also shows temperature profiles for such cylinders when thermal boundary conditions are *T _{i}* = 330 K,

*T*= 300 K,

_{o}*Q*= 5 kW/

_{u}*m*

^{2}and

*Q*= −5 kW/

_{d}*m*

^{2}. Although temperature profiles of UD and X cylinders are almost the same, temperature variation in X cylinder is higher than the other one due to its sharper variation in nanofiller dispersion. In V cylinder, coldness spread is dominant because of the allocation of nanofillers in the outer half thickness of cylinder.

## 5 Conclusions

In this paper, the impacts of graphene and CNT to steady state and transient heat transfer behaviors of FG nanocomposite axisymmetric cylinders were investigated using a mesh-free method under different types of thermal boundary conditions. By employing a micromechanical model, the thermal properties of nanocomposite materials were calculated. The effects of nanofiller content and distribution as well as thermal boundary conditions on the heat transfer behaviors of nanocomposite cylinders were studied which resulted in the following main conclusions:

- The use of nanofiller and the increase of its amount result in shorter stationary times such that the increase
of from 0.1 to 0.3 leads to the reduction of stationary time from 400 s to 100 s.${f}_{r}^{max}$ - The allocation of nanofiller close to the boundaries leads to approaching the temperatures of the nearest boundary.
- The FG allocation of nanofiller intensifies temperature gradient and thermal management of FG nanocomposite cylinders. In one case shown in Figure 6, it was observed that
*T*in V-CNT/PE cylinder was around 356 K while in UD and_{i}*∧*-CNT/PE cylinder*T*=344 K._{i} - The use of graphene in nanocomposite has stronger impact on thermal response than the use of CNT. Comparing Figure 6b and 7b shows that the use of graphene leads to better conduction such that
*T*in graphene/PE cylinders is 10 K less than CNT/ PE cylinders._{i}

The work described in this paper was supported by National Natural Science Foundation of China (Grant no. 11972204) and Natural Sciences and Engineering Research Council of Canada (NSERC under grant RGPIN-217525). The authors are grateful for their supports.

## References

- [1]↑
Power A.C., Gorey B., Chandra S., Chapman J., Carbon nanomaterials and their application to electrochemical sensors: A review, Nanotechnol. Rev., 2018, 7, 19-41.

- [3]↑
Fujii M., Zhang X., Xie H., Ago H., Takahashi K., Ikuta T., Abe H., Shimizu T., Measuring the thermal conductivity of a single carbon nanotube, Phys. Rev. Lett., 2005, 95.

- [4]↑
Samani M.K., Khosravian N., Chen G.C.K., Shakerzadeh M., Bail-largeat D., Tay B.K., Thermal conductivity of individual multi-walled carbon nanotubes, Int. J. Therm. Sci., 2012, 62, 40-43.

- [5]↑
Balandin A.A., Ghosh S., Bao W., Calizo I., Teweldebrhan D., Miao F., Lau C.N., Superior Thermal Conductivity of Single-Layer Graphene, Nano. Lett., 2008, 8, 902-907.

- [6]↑
Fattahi A.M., Safaei B., Ahmed N.A., A comparison for the non-classical plate model based on axial buckling of single-layered graphene sheets, Eur. Phys. J. Plus., 2019; 134, 555.

- [7]↑
Safaei B., Fattahi A.M., Chu F., Finite element study on elastic transition in platelet reinforced composites, Microsyst. Technol., 2018, 24, 2663-2671.

- [8]↑
Sahmani S., Safaei B., Nonlocal strain gradient nonlinear resonance of bi-directional functionally graded composite micro/nano-beams under periodic soft excitation, Thin-Walled Struct., 2019, 143, 106226.

- [9]↑
Ghanati P., Safaei B., Elastic buckling analysis of polygonal thin sheets under compression, Indian J. Phys., 2019, 93, 47-52.

- [10]↑
Moradi-Dastjerdi R., Behdinan K., Stability analysis of multifunctional smart sandwich plates with graphene nanocomposite and porous layers, Int. J. Mech. Sci., 2020, 167, 105283.

- [11]↑
Huang C., Qian X., Yang R., Thermal conductivity of polymers and polymer nanocomposites, Mater. Sci. Eng. R. Reports, 2018, 132, 1-22.

- [12]↑
Gori F., Corasaniti S., Effective thermal conductivity of composites, Heat Mass Transf., 2014, 77, 653-661.

- [13]↑
Balandin A.A., Thermal properties of graphene and nanostructured carbon materials, Nat. Mater., 2011, 10, 569-581.

- [15]↑
Biercuk M.J., Llaguno M.C., Radosavljevic M., Hyun J.K., Johnson A.T., Fischer J.E., Carbon nanotube composites for thermal management, Appl. Phys. Lett., 2002, 80, 2767-2769.

- [16]↑
Huxtable S.T., Cahill D.G., Shenogin S., Xue L., Ozisik R., Barone P., et al., Interfacial heat flow in carbon nanotube suspensions, Nat. Mater. 2003, 2, 731-734.

- [17]↑
Choi S.U.S., Zhang Z.G., Yu W., Lockwood F.E., Grulke E.A., Anomalous thermal conductivity enhancement in nanotube suspensions, Appl. Phys. Lett. 2001, 79, 2252-2254.

- [18]↑
Gu J., Li N., Tian L., Lv Z., Zhang Q., High thermal conductivity graphite nanoplatelet/UHMWPE nanocomposites, R.S.C. Adv., 2015,5, 36334-36339.

- [19]↑
Zabihi Z., Araghi H., Effective thermal conductivity of carbon nanostructure based polyethylene nanocomposite: Influence of defected, doped, and hybrid filler, Int. J. Therm. Sci., 2017, 120, 185-189.

- [20]↑
Moheimani R., Hasansade M., A closed-form model for estimating the effective thermal conductivities of carbon nanotube-polymer nanocomposites, Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci., 2019, 233, 2909-2919.

- [21]↑
Qin Z., Zhao S., Pang X., Safaei B., Chu F., A unified solution for vibration analysis of laminated functionally graded shallow shells reinforced by graphene with general boundary conditions, Int. J. Mech. Sci., 2020, 105341.

- [22]↑
Qin Z., Pang X., Safaei B., Chu F., Free vibration analysis of rotating functionally graded CNT reinforced composite cylindrical shells with arbitrary boundary conditions, Compos. Struct., 2019, 220, 847-860.

- [23]↑
Qin Z., Safaei B., Pang X., Chu F., Traveling wave analysis of rotating functionally graded graphene platelet reinforced nanocomposite cylindrical shells with general boundary conditions, Results Phys., 2019, 15, 102752.

- [24]↑
Shokri-Oojghaz R., Moradi-Dastjerdi R., Mohammadi H., Behdinan K., Stress distributions in nanocomposite sandwich cylinders reinforced by aggregated carbon nanotube, Polym. Compos., 2019, 40, E1918-E1927.

- [25]↑
Moradi-Dastjerdi R., Payganeh G., Tajdari M., Thermoelastic Analysis of Functionally Graded Cylinders Reinforced by Wavy CNT Using a Mesh-Free Method, Polym. Compos., 2018, 39, 2190-2201.

- [26]↑
Moradi-Dastjerdi R., Payganeh G., Thermoelastic dynamic analysis of wavy carbon nanotube reinforced cylinders under thermal loads, Steel Compos. Struct., 2017, 25, 315-326.

- [27]↑
Moradi-Dastjerdi R., Payganeh G., Thermoelastic Vibration Analysis of Functionally Graded Wavy Carbon Nanotube-Reinforced Cylinders, Polym. Compos. 2018, 39, E826-E834.

- [28]↑
Safaei B., Ahmed N.A., Fattahi A.M., Free vibration analysis of polyethylene/CNT plates, Eur. Phys. J. Plus., 2019, 134, 271.

- [29]↑
Safaei B., Moradi-Dastjerdi R., Behdinan K., Chu F., Critical buckling temperature and force in porous sandwich plates with CNT-reinforced nanocomposite layers, Aerosp. Sci. Technol., 2019, 91, 175-185.

- [30]↑
Safaei B., Moradi-Dastjerdi R., Behdinan K., Qin Z., Chu F., Thermoelastic behavior of sandwich plates with porous polymeric core and CNT clusters/polymer nanocomposite layers, Compos. Struct., 2019, 226, 111209.

- [31]↑
Damadam M., Moheimani R., Dalir H., Bree’s diagram of a functionally graded thick-walled cylinder under thermo-mechanical loading considering nonlinear kinematic hardening, Case Stud. Therm. Eng., 2018, 12, 644-654.

- [32]↑
Jalali M.H., Zargar O., Baghani M., Size-Dependent Vibration Analysis of FG Microbeams in Thermal Environment Based on Modified Couple Stress Theory, Iran J. Sci. Technol.- Trans. Mech. Eng., 2019, 43, 761-771.

- [33]↑
Pourasghar A., Chen Z., Effect of hyperbolic heat conduction on the linear and nonlinear vibration of CNT reinforced size-dependent functionally graded microbeams, Int. J. Eng. Sci., 2019, 137, 57-72.

- [34]↑
Pourasghar A., Zhao C., Hyperbolic heat conduction and thermoelastic solution of functionally graded CNT reinforced cylindrical panel subjected to heat pulse, Int. J. Solids Struct., 2019, 163, 117-129.

- [35]↑
Safaei B., Moradi-Dastjerdi R., Qin Z., Behdinan K., Chu F., Determination of thermoelastic stress wave propagation in nanocomposite sandwich plates reinforced by clusters of carbon nanotubes, J. Sandw. Struct. Mater., 2019, DOI: 10.1177/109963621984828.

- [36]↑
Zargar O., Mollaghaee-Roozbahani M., Bashirpour M., Baghani M., The application of Homotopy Analysis Method to determine the thermal response of convective-radiative porous fins with temperature-dependent properties, Int. J. Appl. Mech., 2019, https://doi.org/10.1142/S1758825119500893

- [37]↑
Moradi-Dastjerdi R., Behdinan K., Thermoelastic static and vibrational behaviors of nanocomposite thick cylinders reinforced with graphene, Steel Compos. Struct., 2019, 31, 529-539.

- [38]↑
Moradi-Dastjerdi R., Payganeh G., Transient heat transfer analysis of functionally graded CNT reinforced cylinders with various boundary conditions, Steel Compos. Struct., 2017, 24, 359-367.

- [39]↑
Deng F., Zheng Q.S., Wang L.F., Nan C.W., Effects of anisotropy, aspect ratio, and nonstraightness of carbon nanotubes on thermal conductivity of carbon nanotube composites, Appl. Phys. Lett. 2007, 90.

- [40]↑
Evans W., Prasher R., Fish J., Meakin P., Phelan P., Keblinski P., Effect of aggregation and interfacial thermal resistance on thermal conductivity of nanocomposites and colloidal nanofluids, Int. J. Heat Mass Transf., 2008, 51, 1431-1438.

- [41]↑
Haggenmueller R., Guthy C., Lukes J.R., Fischer J.E., Winey K.I., Single wall carbon nanotube/polyethylene nanocomposites: Thermal and electrical conductivity, Macromolecules, 2007, 40, 2417-2421.

- [42]↑
Hetnarski R.B., Eslami M.R., Thermal stresses-advanced theory and applications, The Netherlands: Springer; 2009.

- [43]↑
Lancaster P., Salkauskas K., Surface Generated by Moving Least Squares Methods, Math. Comput., 1981, 37, 141-158.