Effect of Temperature on the Mechanical Properties of Carbon Composites

The purpose of this research is to investigate and analyze the effect of temperature on the mechanical properties of carbon composites, aluminum foam, and carbon nanotube reinforced aluminum foam. For this analysis, fatigue behavior of aluminum foam has been discussed, and interlaminar shear stress between single-walled carbon nanotube reinforced aluminum foam and multiwalled carbon nanotube reinforced aluminum foam has been compared. The comparison has shown that the interlaminar shear stress within the interface of single-walled carbon nanotube and aluminum foam is higher than that within the interface of multiwalled carbon nanotube and aluminum foam. Hence, the probability of stress concentration and crack initiation within the interface of single-walled carbon nanotube and aluminum foam is higher than that within the interface of multiwalled carbon nanotube and aluminum foam. Furthermore, thermal fatigue lives of different single-walled carbon nanotube reinforced matrix nanocomposites have been evaluated. Additionally, interaction between carbon and molten aluminum has been analyzed. Finally, a new relation for the thermal interlaminar shear stress intensity factor to predict the crack initiation sites on fiber/matrix interface has been introduced.


Introduction
Carbon nanotubes (CNTs) have been introduced by Lijima in 1991 for the first time [1]. Due to the excellent properties of CNT, nowadays its application is broad in many industries such as aerospace and automotive industries. Carbon nanotubes have great mechanical properties such as high tensile strength, high Young's modulus, and high aspect ratio which makes CNTone of the best materials for different applications. Furthermore, the electrical conductivity of CNT is high [2]. Hence, recently, for enhancing the knowledge regarding the CNTs, a few studies to further analyze the mechanics and forest synthesis of CNT have been published [3][4][5][6]. Additionally, for evaluating the effect of CNTs on the mechanical and/or thermal properties of composite materials, many studies have been provided [7][8][9][10].
ere are many methods to produce the CNTs such as catalyst arrays, chemical vapor deposition, electric arc discharge, sonochemistry, laser ablation, and electrolysis. It is important to know that each of these production methods will result in different CNTs with different geometries because the diameter and the length of CNTs will be different with each of these production methods. e different geometries of CNT can have effects on its reactivity, failure mechanism, surface interaction, and mechanical properties. erefore, it seems that, with different production methods, different CNTs with different properties can be expected [2]. e application of CNT is broad in many industries due to its multifunctional properties such as thermal management, enhanced stiffness, and improved toughness [11].
ere are several methods to disperse the CNTs within the matrix. Among these methods, the most effective one is the chemical treatment with ultrasonication process [12].
CNTs are the crystalline shape of carbon. ey have high aspect ratio due to the size of diameter which is a few nanometers and the size of length of it which is typically about a few microns. Furthermore, they offer high flexibility and high thermal conductivity. CNTs can be in the forms of multiwalled carbon nanotubes (MWCNTs) and singlewalled carbon nanotubes (SWCNTs) [12].

Mechanical Properties, Geometries, and Failure
Mechanism. As it has been mentioned previously, the discovery of CNTs occurred in 1991. After the discovery of CNTs, their application started to increase due to the unique and important properties such as high strength, high modulus, high stiffness, electrical conductivity, etc. Compared to steel which is one of the most popular metals in different industries, CNTs offer higher Young's modulus and tensile strength. Young's modulus of steel is about 200 GPa, while Young's modulus of CNTs is around 1 TPa. e tensile strength of steel is about 505 MPa while the tensile strength of CNTs is about 1163 GPa. Moreover, CNTs can be used in other applications such as energy storage devices, hydrogen storage media, sensors, and many more [13].
CNTs are composed of graphene layers rolled to form cylinders with fullerene caps. SWCNT has only one graphene sheet while MWCNT has several graphene sheets rolled over each other with the same axis. e geometries of CNTs could be in three states: 1. armchair; 2. zig-zag; and 3. chiral. CNT thickness is 0.34 nm which is equal to the thickness of graphene sheet, but the length and dimeter of the CNT can vary. In Figure 1, the three different geometries of CNTs are illustrated [13].

Ranking of Unidirectional Fibers/Matrix Composites Based on Their Interlaminar
Shear Stress e aim of this section of the study is to rank the unidirectional fibers/matrix composites based on the interlaminar shear stress existing between fibers and matrix due to thermal stress. ermal stress induced in composites is one of the main issues in many applications such as space structures. In this section, by applying an analytical method, it is attempted to estimate the ranking of unidirectional fibers/matrix composites exposed to thermal stress based on their interlaminar shear stress between the fibers and matrix. e application of the results of this research is very broad. ese results could be very advantageous in any industry using unidirectional fibers/matrix composites exposed to thermal stress such as aerospace automotive, etc. ermal stress is one of the main issues in application of composites in many applications such as aerospace industry. Space structures are exposed to thermal cycles that induce thermal stress in composite structures in space. In order to investigate this effect, numbers of studies have been performed by Funk and Sykes [14], Shin et al. [15], and Unigovski et al. [16]. e objectives of the mentioned studies were to understand the damage induced due to thermal cyclic exposures.
In 2012 [17], Park et al. have conducted an experiment to simulate the thermal cycling in Low Earth Orbit (LEO) environment and estimate the effect of that on Unidirectional Carbon Fiber/Epoxy Composite (UD CF/EP). UD CF/ EP material is one of the composite materials that is applied in aerospace structures such as satellites [17]. e objective of this experiment was to estimate the deterioration of UD CF/EP mechanical properties in LEO due to thermal cycling exposure.
Additionally, in 2018 [18], the experimental results achieved by Park et al. [17] have been analyzed. e analysis of the experimental results had indicated that Interlaminar Shear Strength (ILSS), Flexural Strength (FS), and Flexural Modulus (FM) of UD CF/EP were decreasing as thermal cycles were increasing in LEO simulation experiment. Furthermore, based on the experimental results, it has been approximated that the deterioration of ILSS in UD CF/EP can be the main cause of fracture in UD CF/EP exposed to thermal stress.
Nevertheless, it appears that there is no study to rank the unidirectional fibers/matrix composites (UFMCs) based on their interlaminar shear stress (ILSs). In the presented study, by applying an analytical method, the amounts of ILSs in UFMCs are calculated. Based on the assumption which has been approved in 2018 [18], the higher the ILSs is, the higher the probability to cause the fracture within the UFMC becomes. us, the best UFMC in terms of minimum probability of fracture or more durability when exposed to thermal cycling is the one with the minimum ILSs. As a result, this calculation would contribute to ranking the UFMCs based on their minimum ILSs (higher durability or higher thermal fatigue life).

Interlaminar Shear Stress between the Fibers and Matrix
Interfaces. In 2018 [18], a relation has been proposed to estimate the ILSs within the UFMCs interfaces between the fibers and matrix. is equation is indicated as follows: (1) In this part of the manuscript, equation (1) [18] is employed to estimate the maximum ILSs. Based on the results obtained in 2018 [18], deterioration of ILSS in UD CF/EP seems to be the main cause of failure of this composite material which has been exposed to thermal cycling. us, in this section, it has been tried to establish a system to calculate the ILSs in common UFMC materials.
is interpretation would verify that the risk of failure in UFMC is increased while ILSs is increased. erefore, in this part of the article, ILSs in common UFMCs are calculated in order to rank them. It is important to notice that equation (1) is only valid in cases which fibers are impeded and constrained within the matrix.
In order to calculate the maximum ILSs (ILSs max ), it is enough to replace the shear modulus with maximum shear modulus in equation (1). us, equation (1) becomes In the equation (2), α c is the Coefficient of ermal Expansion (CTE) of matrix in axial direction, α F is the CTE of fiber in axial direction, and ΔT is the difference between the stress-free or crack-free temperature of the UFMC and the environment temperature. Furthermore, G max is the maximum shear modulus. It is important to notice that G max is determined by comparison between the shear modulus of matrix and fiber. e largest value of G along the fibers' direction between the matrix and fiber needs to be identified and be substituted in equation (2), in order to obtain ILSs max .
In Tables 1 and 2, there are lists of all common matrices and fibers used to develop UFMCs, respectively. In the Results and Discussion section, with applying equation (2), ILSs max of possible UFMCs exposed to thermal stress has been calculated. Table 3, ILSs max for 54 types of UFMCs are indicated as a function of ΔT. In evaluation of the ILSs max for ranking the composite materials, the amount of ΔT is not required because it is assumed to be of the same value for all the materials and is constant for all of them. In order to use the values of ILSs max in Table 3, it is enough to have this knowledge that, based on the results obtained in 2018 [18], ILSs could be the main cause of fracture and/or de-bonding within the UFMCs. Consequently, as the value of ILSs increases, the probability of de-bonding or breaking bonds between fibers and matrix increases. is phenomenon will lead to crack initiation, propagation, and ultimately fracture within the UFMC.

Results and Discussion. In
Based on equation (2), as the difference between the CTEs of fibers and matrix increases, ILSs max increases. Furthermore, as the shear modulus increases, ILSs max increases. us, ILSs max is proportional to the difference of CTEs between fibers and matrix, and maximum shear modulus between fibers and matrix in axial direction.
It seems that the mismatch of CTEs between fibers and matrix is a great cause of failure, de-bonding, or fracture between those. e reason behind is that the axial CTE of the matrix is a positive value and the axial CTE of fiber, in most cases, is a negative value. Consequently, while UFMC is heating, the matrix is expanding in axial direction. On the other hand, fiber is contracting in axial direction. Additionally, while UFMC is cooling, matrix is contracting. On the other hand, fiber is expanding in this state in axial direction. is reverse or mismatch behavior can cause inducing ILSs between fibers and matrix. If this process continues, it will create thermal cycles due to repeated cooling and heating. As a result of thermal cycling due to this    [21] 16.60 [21] reverse behavior between fibers and matrix, bonds between them will start to break gradually and will create crack initiation along the axial direction of fibers within the fibers/ matrix interfaces. ese cracks will grow as thermal cycle numbers increase and, ultimately, will cause fracture or failure within the UFMC. In Table 3, unidirectional diamond fibers/CE339 epoxy matrix has the highest ILSs max . It means that the probability of fracture due to thermal stress in this composite is the maximum. is fact is due to the high value of the axial CTE of CE339 epoxy matrix, and low axial CTE of diamond fibers. Additionally, the shear modulus of diamond is very high. us, ILSs max , is highest for this UFMC material. erefore, application of this UFMC material in low and high temperature environments with respect to crack-free temperature is not recommended. On the other hand, application of Sic fibers/borosilicate glass matrix composite is recommended because it has a minimum ILSs max compared to other UFMCs. Moreover, as it is indicated in Table 3, ILSs max for most of the unidirectional Sic and diamond fibers composites is very high. It means that the application of diamond and Sic fibers, in most cases, is not recommended in thermal stress condition. e reason is the high shear modulus of these fibers and the mismatch of CTEs between these fibers and the matrix.

Fatigue Behavior of Aluminum Foam
Previously, compression-compression and tension-tension fatigue experiments for Alulight Aluminum foam have been evaluated. Stress-fatigue cycles curve has been derived for closed cell foam samples with compositions Al-1Mg-10Si and Al-1Mg-0.6Si (wt.%). e results have shown that the fatigue life is very sensitive to the value of the stress which is applied to the sample [23]. Additionally, for the evaluation of the fatigue behavior of the aluminum foam, "Uniaxial Stress-Strain Behavior of Aluminum Alloy Foam" is submitted by McCullough et al. [24], "Compression Fatigue of a Cellular Al Alloy" is presented by Sugimura et al. [25], "On the Mechanical Performance of Closed Cell Al Alloy Foams" is provided by Sugimura et al. [26], "Fatigue Failure of an Open Cell and Closed Cell Aluminum Alloy Foam" is investigated by Harte et al. [27], and "Fatigue of Aluminum Alloy Foam" is submitted by Olurin et al. [28].
Another compression fatigue experiment has been done on aluminum-silicon foams under uniaxial cyclic stresses ranging from 50 to 80% of the total compression strength of the samples. For producing these samples, aluminum powder was combined with 7 wt.% silicon powders. Furthermore, 0.5% titanium hydride has been added to the sample. e samples had a cylindrical shape with height   [20] 190.7 [20] equal to 60 mm and diameter equal to 44 mm. e produced foams had densities 0.55, 0.63, and 0.94 g/cm 3 [29]. Furthermore, another fatigue experiment with compression-compression loading was performed on closed-cell aluminum foam in 2008. Observations have shown that noticeable hardening under loading may occur which is accompanied with plastic deformation within quasi-static loading [30].
It seems that the main reason for applying Al foam in different industries is due to its lightweight. Aluminum foam Sandwich is one of the panels that could be used in structural applications. Hence, shear fatigue behavior of this material has been investigated. e experiment has shown that mechanical behavior of this Al foam Sandwich seems to be highly dependent on cells distribution, morphology, and manufacturing process [31].
Metal foams are produced and developed due to their unique properties such as chemical reactions, damping capacity, energy absorption, and low stiffness to density ratio. Fatigue can occur in foams even at low applied stress due to the existence of maximum stress located on the foam microstructure nodes. Fatigue can affect the material parameters and macroscopic behavior of the specimen due to the influence imposed by the different number of stress cycles and stress amplitude [32].
In 2013, another compression-compression fatigue experiment has been performed. e objective of the cyclic loading experiment was to determine the load amount for failure cycles of 10 5 . e results have revealed that application of 50% of compression strength for the cyclic load is very low for obtaining 10 5 cycles for the fatigue life [33].
Metal matrix foams can be applied as load bearing, energy absorber, and vibration damping in structural parts and machine beds. Metal matrix foams show good stability under compression cyclic loads. Deformation of these foams under compression cyclic loads occurs in three stages. First stage is the linear increase of strain with increasing fatigue cycles. Second stage is the accumulation of strain in many numbers of cycles. ird stage is the rapid strain accumulation during few cycles until failure occurs. e experiment also has indicated that fatigue limit in case of tension is twice less than in the case of compression loading [34].

Interlaminar Shear Stress Comparison between SWCNT/Aluminum Foam and MWCNT/Aluminum Foam
In the presented part of the study, ILSs within the interfaces between the fibers and matrix, by using experimental data and applying analytical method, between SWCNT/aluminum foam and MWCNT/aluminum foam are compared. In order to compare the ILSs of these two nanocomposites, the following procedure has been followed. e mismatch of coefficients of thermal expansion between fibers and matrix in both nanocomposites is compared. According to the results obtained, it appears that MWCNT/aluminum foam can offer lower ILSs when compared to SWCNT/aluminum foam. e reason is, based on the evidence achieved by the mentioned analytical method, SWCNT/aluminum foam, in average, offers higher mismatch between coefficients of thermal expansion and interlaminar shear stress within the fibers and matrix interfaces. us, further stress concentration in these areas may result in higher crack propagation rate that can cause fracture. e application of aluminum foam is nowadays very broad in many industries such as transport vehicles. e reason behind is the mechanical properties of aluminum foam such as fire resistance, acoustic insulator, shock absorption, low thermal conductivity, suppress vibrations, and mitigation explosion effects. Furthermore, aluminum foam is now a popular material for advanced material for structural purpose due to its lightweight, crashworthy, thermally insulate, and cheap price. Additionally, aluminum foam offers novel physical, mechanical, and electrical properties along with low density, which makes it a suitable material for aeronautical industry. Moreover, this material can be used as a cellular material that may be applied to different applications such as cores for Sandwich panels, energy absorption elements for automotive industries, and many multifunctional structural elements in vehicles [35].
Increasing ILSs within the fibers and matrix interfaces can result in stress concentration in these areas. e stress concentration may cause crack initiation, propagation, and de-bonding between the fibers and matrix.
It appears that there is no research related to comparison of ILSs between SWCNT/aluminum foam (SWCNTAF) and MWCNT/aluminum foam (MWCNTAF). erefore, in the presented part of the study, by using the assumptions in Section 2, and applying experimental data [32,33], ILSs between SWCNTAF and MWCNTAF are compared. e results of this research may contribute to selecting the best nanocomposite between these two, to apply in thermal cycling environment.

Experimental Procedures
4.1.1. CNT Experimental Procedures. "Spinnable MWCNT arrays were obtained by chemical vapor deposition using C 2 H 2 and FeCl 2 as the base material and the catalyst, respectively. e diameter of the MWCNTs was measured using a transmission electron microscope (TEM, JEOL JEM-2100F, Japan). A partially cured epoxy resin (B-stage epoxy) with a release paper was used as the starting material, where the epoxy resin comprised bisphenol-A type epoxy, novolactype epoxy, and an aromatic diamine curing agent. e epoxy resin was then impregnated into the MWCNT monolithic sheet at 90°C for 3 min between the steel plates of a hot press (AS ONE AH-4015, Japan). After peeling off the release paper from the MWCNT sheet now impregnated with the epoxy resin (prepreg sheet), the prepreg sheet was cured at 130°C for 1.5 h at a pressure of 1 MPa using the hot press, forming a film specimen" [36].

Al Foam Experimental
Procedures. "Another way of foaming the melts is to directly add a foaming agent to the molten metal. e foaming agent releases gas by Journal of Engineering decomposition under the influence of heat which then carries out the foaming process. Shinko Wire Company, Japan has been manufacturing foamed aluminum under the registered name "Alporas." is process initially involves addition of 1.5% of calcium to the molten aluminum for adjusting the viscosity and secondly addition of molten aluminum at 670°C to 690°C, stirred for 6 to 10 minutes in an ambient atmosphere. It thickens the aluminum melt which is then poured into a casting mould and stirred with addition of powdered TiH2 (foaming agent) by using a rotating impeller. Addition of a sufficient amount of hydride (usually 1.6%) decomposes the melt under the influence of heat and it releases hydrogen gas. As a result, the foam expands and fills up the mould within 15 to 20 minutes. It is cooled by air (at room temperature) in the mould and solidifies as a block with porosity between 85% and 95%. e production rate of Alporas is reported to be 1000 kg/day by application of the batch casting process.
is process is capable of producing large blocks of good quality foams. is process is capable of producing large blocks of good quality foams such as Blocks of 450 mm, the Alporas process, and GASCAR process. e required width, 2050 mm length, and 650 mm height can be produced with this technique. ese foams have uniform pore structure and they are better than the foams which are produced with addition of ceramic particles. e latter process makes foams brittle. However, addition of hydride, which requires complex processing equipment, makes the process more expensive than the one in which foaming of melts is carried out by gas injection. e density range of these foams is 0.18-0.24 gcm −3 and the mean pore size is about 4.5 mm. Addition of calcium powder in aluminum melt enhances the viscosity of the melt. With increased viscosity, flow of aluminum melt slows down, resulting in its reduced drainage before solidification" [37]. Table 4, shear modulus and CTEs of SWCNT, MWCNT, and aluminum foam (AF) at temperature range of −5 to 85°C are indicated. For performing the thermal analysis to compare the ILSs between MWCNTAF and SWCNTAF, the data in Table 4 and equation (1) [18] are used.

Problem Formulation. In
As it has been mentioned in Section 2 of the manuscript, equation (1) can contribute to estimate the probability of crack initiation and/or propagation within the SWCNTAF and MWCNTAF. With this assumption, by employing the data in Table 4 and using equation (1), ILSs in both SWCNTAF and MWCNTAF can be approximated. Figures 2-5, most of the results are the evidence of lower ILSs within the MWCNTAF in comparison to that within the SWCNTAF. e reason is that the average mismatch of CTEs between SWCNT and AF is 0.88% more than that between the MWCNT and AF. Hence, maximum ILSs within the SWCNTAF is 399% more than that within the MWCNTAF. Additionally, average ILSs within the SWCNTAF is 349% more than that within the MWCNTAF. e only result that makes the opposite conclusion is the difference of CTEs between SWCNT and AF that is 19% less than that between MWCNT and AF. It means that, based on the results obtained in this study, it appears that crack initiation and/or propagation on interfaces between the nanofiber and matrix within the SWCNTAF may occur with higher probability when compared to that within the MWCNTAF. is result suggests using MWCNTAF instead of SWCNTAF in cases dealing with thermal cycling condition for having a longer thermal life.

Enhancing Thermal Fatigue Resistance Using
Unidirectional SWCNT/Aluminum Foam Composite Material e significant purpose of this section of this research is to use a method to compare the thermal fatigue resistance of unidirectional SWCNT/aluminum foam with unidirectional SWCNT/epoxy. e reason is that the application of SWCNT/epoxy is nowadays very common in many industries such as aerospace and automobile due its lightweight and high strength. In this part of the study, by using analytical method, the interlaminar shear stresses on interfaces of fibers/matrix in both nanocomposites are derived and compared. With the analysis of these results, it is determined that the risk of crack initiation and/or propagation due to stress concentration on interfaces of fibers/matrix is higher within the SWCNT/epoxy in comparison to that within the SWCNT/aluminum foam. us, thermal fatigue resistance of these two nanocomposites could be compared with this analysis. e results have shown that the thermal fatigue life could be enhanced by application of unidirectional SWCNT/aluminum foam because interlaminar shear stress within the SWCNT/aluminum foam is less than that within the SWCNT/epoxy. ese results are very advantageous because replacing epoxy with aluminum foam could decrease the manufacturing costs due to lower price of aluminum foam.
In order to have the best design, it seems enough to try to have a best match for CTEs between the nanocomposites components such as fibers and matrix. It means that, according to this theory, fibers and epoxy should have similar values of CTEs to offer less interlaminar shear stress (ILSs) and higher thermal fatigue life.
Many studies have been performed to estimate the thermal fatigue life of nanocomposites and composites such as UD CF/EP. Nevertheless, it seems that there is no study to evaluate the thermal fatigue life of unidirectional SWCNT/ aluminum foam (USWCNTAF).
In the presented part of the study, by applying analytical method, it is tried to compare the thermal fatigue life of USWCNTAF with unidirectional SWCNT/epoxy (USWCNTE).
e results have shown that USWCNTAF could offer higher thermal fatigue resistance when compared with USWCNTE.

Problem Formulation.
According to the results of the research in 2018 [18], it seems that, with increasing ILSs, thermal fatigue life decreases. In Table 5, CTEs and shear modulus of SWCNT, epoxy, and aluminum foam at temperature range −5 to 85°C are indicated. With the presented data in Table 5 and using equation (1), in the Results and Discussion section of this part of the manuscript, thermal fatigue resistance of USWCNTAF and that of USWCNTE are compared.       Figures 6-9, it can be observed that at temperature range of −5 to 85°C the maximum mismatches of CTEs and ILSs are 49% decreased in SWCNTAF in comparison with that within the SWCNTE. Additionally, mean mismatch of CTEs within the SWCNTAF is 40% decreased in comparison with that within the SWCNTE. Moreover, average ILSs in SWCNTAF is 44% decreased in comparison with that in SWCNTE.

Results and Discussion. By analyzing the results which are shown in
Finally, by summing up all the results, it can be concluded that, due to the decrease of CTEs mismatches and ILSs in SWCNTAF in comparison with that in SWCNTE, thermal fatigue life in SWCNTAF seems to be higher than that in SWCNTE. e reason is that, by decreasing the CTEs mismatches and ILSs, the stress concentration on interfaces between SWCNT and AF decreases. By decreasing the stress concentration on interfaces between SWCNT and AF, the probability of crack initiation and/or propagation decreases. us, the chance of de-bonding between SWCNTs and AF decreases. is will result in postponing the fatigue failure or fracture. Hence, the thermal fatigue life of SWCNTAF could be higher than SWCNTE.

Interaction between Carbon and Molten Aluminum
In 1990, Metal-Matrix Composites (MMCs) were introduced to use in structural applications. ey were also recommended to use for thermal management applications. eir characteristics were 1. high temperature operating, 2. improved strength, 3. higher stiffness, 4. higher thermal conductivity, 5. abrasion resistance, 6. creep resistance, and 7. dimensional stability [38].

Aluminum-Matrix Composites. Aluminum-Matrix
Composites were used due to their lightweight, environmental resistance, and enough melting point for many   Journal of Engineering applications. Furthermore, aluminum could accommodate with different fibers, including AL 2 O 3 , Sic, boron, graphite fibers, and many other particles [38]. In 2000, Pippel et al. performed a research to characterize the interlayer structure of carbon fiber reinforced aluminum wires. According to the results of their research, the structure of the interlayer between fiber and matrix has a thickness of 2-4 nm and consists of carbon enriched with aluminum. Within the layer, 1-2 nm grains are sometimes observed. ese grains have a lattice plane distance around 0.26 nm. Furthermore, the interlayers contain needles or platelets which are aluminum carbide. Based on the experiments, "the growth of the carbide needles or platelets starts immediately from the thin amorphous aluminumcontaining surface layer of the fiber" [39].
In the condition that carbon and aluminum are in contact with each other for several minutes, chemical reaction could occur between them and contact angle between them will decrease.
is results in producing aluminum carbide (Al 4 C 3 ). is production is called wetting. is can happen at temperatures above 550°C, but in temperatures above 900°C reaction can happen at a faster rate and it decomposes to Al 2 O vapor [40].
Aluminum (Al) carbide is a hygroscopic ceramic which is brittle and has a high hardness. us, it can act as a stress concentrator. As a result, it fails at relatively low strain. Aluminum carbide also is vulnerable to galvanic corrosion. In axial direction along the fibers, the CTE of carbon fiber (CF) is negative, while the CTE of Al is the same positive value in all directions. As a result, in manufacturing process of CF/Al composites when the temperature drops below 340 o C thermal stresses will appear. In order to prevent forming Al carbide, fiber coating, less reactive fibers, and limiting time in manufacturing process are recommended [40].
e results have shown that melting point of CNT is around 3500°C and MWCNTs are in a stable condition even in temperatures higher than 2800°C. Nevertheless, two failure mechanisms are most likely to occur. 1. If MWCNT is stronger than matrix and interface, the nanocomposite will fail pulling out of MWCNTs 2. If matrix and interface are stronger than MWCNT, the failure occurs through the stochastic breakage of MWCNT at the weakest sites along its length [41].
In solidification of MMC, reinforcement is combined with the liquid metal. e solidification process can be divided into two parts: 1. stir mixing; 2. melt infiltration [42].
"Defects in activated carbon and in CNTs provide sites for chemical reactions to occur. Reactions of CNT with metal matrix leading to carbide formation have been observed by many researchers. e small size and amount of Al 4 C 3 formed was due to the smaller availability of defective sites and amorphous carbon. e formation of carbide also depends on the processing techniques. Some researchers have reported no formation of Al 4 C 3 in the case of solid-state processes such as extrusion. On the other hand, there would be no existence of CNT in the matrix when heated up to 1073 K and XRD results show that CNT fully converts to Al 4 C 3 . e reaction at the triple point between liquid alloy and CNT leads to formation of carbide and subsequent spreading of metal. Minimal reaction of CNT is desirable in order that efficient stress transfer can occur without much damage to the CNT structure" [43].
" e application of a sputter-coated layer of aluminum on the CNTs resulted in the encapsulation of all exposed nanotubes for light deposition and the formation of a 1 μm-thick composite layer for longer deposition, a promising result which strongly suggests an affinity between the two materials.
is observation highlights the potential for non-destructively impregnating the CNT forests with aluminum, a technique demonstrated previously with a polymer matrix" [44].
" e results show a good distribution of the matrix between the CFs when using higher casting pressures of 38 and 53 MPa, lower pouring temperature of (700°C), and lower die temperatures of 100 and 200°C. is will result in higher ultimate tensile strength of CF/Al-Si because of the increased bonding, decreased shrinkage defects, and fibers degradation" [45].
Due to the high specific strength, high modulus, and low CTE in CF/Al composites, they have been used in electronic and space industries. However, the reaction between CF and Al may take place in CF/Al composites [45]: Interaction between Al and CF can produce Al 4 C 3 which is an interface drawback that is related to MMCs manufactured by liquid state procedure. Formation of Al 4 C 3 on the interfaces degrades the mechanical properties of the CF/ Al composite [45].
CF/Al composites are not easily recyclable due to the reaction of carbon with molten aluminum or the reaction of carbide with molten aluminum. e results have shown changes of enthalpy at temperatures over 650°C due to the phase change as aluminum enters the molten state. Furthermore, the experiments have shown that, at temperature 850°C, molten aluminum interacts effectively with Sic on their interfaces.
is will cause dissolving silicon that will produce aluminum and aluminum carbide [46]. e results of experiments have shown that, at 1.5 wt.% MWCNTs, elongation percentage and ultimate tensile strength are maximum values in MWCNTs/A356 aluminum alloy. MWCNT/A356 aluminum alloy can be applied in machinery, aircraft, and automobile industries due to its low-cost aluminum and mechanical properties. In automobile industry, this nanocomposite can be used to replace steel components and applications such as pistons, cylinder liners, and connecting rods [47].
While producing CF/Al alloy composite, Al 4 C 3 carbides may develop on the surface of CFs. Furthermore, at pressure infiltration process, CFs can burn and release gasses that could cause crack initiation due to thermal stress. Fiber clustering, poor wettability, chemical reaction on fiber/liquid aluminum interfaces, porosities, and non-homogeneity in CFs distribution are the restrictions of applying CF/Al alloy composites in many industries [48].
Low interfacial strength in matrix/fibers areas, intermetallic formation, poor infiltration, and low tensile strength are considered as main issues in CF/Al composite fabrication. Al 4 C 3 particles present at interfaces can cause cracking in CF coating under high pressure. is is due to the nucleation of carbide and heterogeneous growing from CFs into the Al due to interfacial reaction. e other challenge is the reaction of CF with molten Al that will cause air bubbles arrested in Al. ese bubbles can act as crack formation sites under loading. In some circumstances, the formation of carbide (Al 4 C 3 ) on interfaces can cause flat fractures areas without fiber to pull out due to brittleness of the CF/Al composite [49].
Low wettability between CNTs and Al, oxidation of molten metal, and consequent reduction in mechanical properties are the challenges to produce CNT/Al composite. Argon as an inert gas can be applied in fusion process to prevent the oxidation in aluminum [50].
Experiments have shown that in CNTs/Al composite self-ion irradiation can convert the CNTs to Al 4 C 3 which is metastable. Uniform dispersion of CNTs in Al can build internal interfaces which are prolific in Al that may cause venues in radiation defects situations to self-heal (recombine). Irradiation experiments with Al ion and helium have shown that if CNTs can be dispersed uniformly in Al, they can decrease the embrittlement and hardening due to radiation [51]. e best vol.% of CNTs in CNTs/Al composite is 1% because the number of clusters and pores are less, and the hardness is maximum. Furthermore, at this amount of CNT, CNT clusters have smaller size. Moreover, higher value of CNT may cause CNTs non-uniformity in CNT/Al composite, and softening due to larger clusters. e existence of carbides on interfaces of CNT/Al can cause higher strengthening [52].
It is very significant to have knowledge regarding the interaction between CNTs and Al due to the load transfer characteristics of this region. Weak bonding in this interphase area can be due to the wetting. Furthermore, CFs cannot usually have chemical interaction with Al.
us, this can result in low-strength bonding on interfaces. However, in cases exposed to high loads and thermal stresses, application of continuous fibers to induce higher strength in composite is recommended [53].
Formation of carbides on CFs surface can develop pits which can act as stress concentration and crack initiation sites which may cause fracture. Application of some carbide coatings on CFs can protect them from molten Al attack at high temperatures. e layer of carbide deposited on CFs  can be 50-500 nm thick which depends on the CFs characteristics [54].
Experiments have shown that increasing load can result in crack growth due to coalescence of cavities oriented from the surface into the center of specimen.
ere are two reasons for crack initiation: 1. rupture of large particles, 2. decohesion. Decohesion can result due to the difference of physical properties in different phases. Al has higher CTE and lower elastic modulus when compared to Al 4 C 3 and Al 2 O 3 which can cause high stress gradients on interfaces that could result in thermal residual stresses. Consequently, on interfaces, further dislocation density may appear. Rupture as a result of crack initiation in loading direction will occur on the surface which has higher dislocation density with the presence of higher amount of oxide (Al 2 O 3 ) and carbide (Al 4 C 3 ) particles [55]. Al 4 C 3 and Al 2 O 3 can form as shells during the CF/Al reactions. Al 4 C 3 particles can be either 30 nm or 0.6 μm, and Al 2 O 3 particles may be as large as 1 μm. Al 4 C 3 and Al 2 O 3 particles are in the shape of elongated and spherical, respectively. CF/Al can remain stable up to temperature 620°C [56].
e results have indicated that the formation of Al 4 C 3 in CF/Al composite deteriorates the mechanical properties and thermal conductivity. Due to the reaction between CF and Al at high temperatures (773 K), Al 4 C 3 will form and the formation will continue over time which results in CFs damaged surface and reduction of diameter. e growth of carbide (Al 4 C 3 ) can occur during the infiltration and cooling process especially when the exposure time increases [57].

Thermal Interlaminar Shear Stress Intensity
Factor to Predict Crack Initiation Sites e application of Metal-Matrix Composites (MMCs) is very broad in many industries, especially aerospace and automotive. One of the most applicable MMCs in industries is aluminum matrix composite which contains carbon fiber and/or carbon nanotube. However, in manufacturing carbon fiber/aluminum composites or carbon nanotube/aluminum nanocomposites, carbon can react with molten aluminum. e results of this reaction can lead to developing carbide (Al 4 C 3 ) and/or aluminum oxide (AL 2 O 3 ) in carbon/ aluminum interfaces. Hence, on the interface of carbon/ aluminum, carbon, carbide, aluminum, and aluminum oxide could be in contact with each other. e goal of this part of the study is to develop a method to identify the crack initiation site within the carbon/aluminum composites when exposed to thermal stress. In order to achieve appropriate results in this study, the new concept of thermal interlaminar shear stress intensity factor is used. e results have shown that the probability of crack initiation on interfaces of aluminum/aluminum oxide is the highest. e aluminum symbol is Al and is located in boron group. Al metal is soft, white, and ductile. Al is the most abundant metal located in crust of Earth. Al has a low density and is a corrosion-resistant metal. Due to its mechanical properties, it has been widely used in structural elements, transportation, and aerospace industries [58].
MMCs offer structural and thermal applications. ese materials are nonflammable and do not outgas in vacuum environment. Furthermore, their interaction with solvents and fuels is minimum. Moreover, they have higher mechanical properties such as strength, thermal conductivity, and creep resistance when they are compared to their metal counterparts [38].
Application of CNTs in reinforcing MMC occurred in 1990s. e reason for using CNT is its excellent thermal, mechanical, and electrical properties. Nowadays, CNT is used in many applications such as supercapacitors, interconnects, nanoscale transistors, molecular vehicles, artificial muscle fibers, etc. CNT strongly interacts with aluminum to create MMC. Furthermore, experiments have proved that CNT/aluminum nanocomposite has higher electrical, thermal, and mechanical properties when compared to pure aluminum. However, the wetting angle between CNT and aluminum is 130°-140° [44].
Carbon fiber/aluminum composites offer high modulus, specific strength, and low CTE. Hence, they have been frequently used in electronics and space applications. However, it seems that developing carbide on interfaces of carbon fiber/ aluminum is a drawback. e reason is carbide acts as a brittle material. us, it has a deteriorating effect on interfaces which degrade the mechanical properties of carbon fiber/aluminum. Nevertheless, this composite material has many applications in automotive and aerospace industries due to its high strength-to-weight ratio and low cost [45]. e application of aluminum to manufacture MMC is very broad. Aluminum is being used in many industries such as nautical, automotive, and aeronautical due to its low weight [50].
Experiments have shown that, in manufacturing process of carbon fiber/aluminum or CNT/aluminum at around 600°C, carbon interacts with molten aluminum. As a result of this interaction, carbide (Al 4 C 3 ) and/or aluminum oxide (Al 2 O 3 ) can be developed on interfaces. Aluminum has a lower modulus (70 GPa) than carbide (445 GPa) and aluminum oxide (393 GPa). On the other hand, the CTE of aluminum (26.5e − 6 1/K) is higher than carbide (5e − 6 1/K) and aluminum oxide (8.3e − 6 1/K). Consequently, when carbon/aluminum composites are exposed to thermal stress, due to difference of physical properties between these materials, de-cohesion might occur. Due to the existence of thermal residual stresses and mismatch of CTEs between these materials, stress gradients can develop over the interfaces which could result in crack initiation [55].
As it is mentioned in the previous paragraphs, in manufacturing carbon fiber/aluminum composite or CNT/ aluminum nanocomposite, carbide and aluminum oxide may develop at temperatures around 600°C. As a result, these materials could be in contact with each other on carbon/ aluminum interfaces.
It appears that up to this time there is no research to investigate the crack initiation sites over the carbon/aluminum interfaces while they are exposed to thermal stress. Hence, the goal in this section of the study is to develop a method to identify the crack initiation sites. In order to perform this evaluation, the new concept and relation of thermal interlaminar shear stress intensity factor is employed.

Problem Formulation.
In 2018 [18], an equation was proposed to calculate the interlaminar shear stress (ILSs) in composites and nanocomposites. As it is obvious, ILSs can be developed on fibers/matrix interfaces due to temperature variation. is equation is indicated as follows: In equation (4), Δα is the difference between the axial coefficients of thermal expansion of contacting materials, G max is the maximum shear modulus between the contacting materials, and ΔT is the difference between the stress-free temperature and ambient temperature.
In this part of the study, the aim is to develop a method to identify the crack initiation site within the carbon/ aluminum composites exposed to temperature variation. Previously, it was mentioned that, due to the interaction between carbon and molten aluminum, carbide and aluminum oxide could be developed on fiber/matrix interface. us, carbon, aluminum, carbide, and aluminum oxide could exist over the interface of fiber/matrix. Hence, they can be in contact with each other. According to equation (4), ILSs due to the temperature variation can be calculated with having, first, the mismatch of CTEs between the two materials, second, with having the temperature variation from stress-free temperature, and, third, with maximum shear modulus between the two materials.
With using equation (4), ILSs for these four materials, carbon, aluminum, carbide, and aluminum oxide that are one by another in contact with each other, can be calculated. It seems obvious that higher ILSs can cause higher stress concentration over the interfaces and consequently can cause higher probability of crack initiation.
Additionally, with closer investigation in equation (4), it can be inferred that ΔT is the same value for all the contacting materials in this research as they are part of the carbon/aluminum composite interfaces.
us, based on equation (4), it appears that the only factor that can cause stress concentration due to the temperature variation is the multiplication of Δα with G max (Δα·G max ). According to this assumption, as Δα·G max increases, the chance of crack initiation and/or propagation increases. Hence, it seems that Δα·G max represents thermal stress intensity factor on interfaces. As a result, this relation could be called "thermal interlaminar shear stress intensity factor" and can be shown as K TILSs . erefore, K TILSs is equal to Δα·G max as it is shown in e unit of K TILSs is N/k·m 2 and can be used to evaluate the thermal stress intensity on interfaces of two contacting materials. In this research, the objective is to identify the crack initiation and/or propagation sites over the carbon/ aluminum interfaces. erefore, the CTEs and shear modulus of aluminum, carbon, carbide, and aluminum oxide are required to calculate K TILSs . ese values are indicated in Table 6. ese data, with applying (2), can contribute to calculating the thermal stress intensity on interface of contacting materials. e highest value of K TILSs can represent the location with highest probability of crack initiation and/or propagation in carbon/aluminum composite exposed to thermal fatigue, thermal stress, or temperature variation. In Table 7, the values of K TILSs for all contacting materials on carbon/aluminum interface are indicated. Table 7, all the values of K TILSs for all the materials interacting at carbon/aluminum interfaces are indicated. e ranking of K TILSs from highest to lowest values is as follows:

Results and Discussion. In
(1) Aluminum interacting with Al oxide with K TILSs � 2,984,800 N/k·m 2 (2) Aluminum interacting with carbide with K TILSs � 2,773,500 N/k·m 2 (3) Al oxide interacting with CNT with K TILSs � 1,607,200 N/k·m 2 (4) Carbide interacting with CNT with K TILSs � 838,500 N/k·m 2 (5) Al interacting with CNT with K TILSs � 728,000 N/ k·m 2 (6) Carbide interacting with Al oxide with K TILSs � 541,200 N/k·m 2 e mentioned ranking of K TILSs for two materials in contact may also represent the ranking of stress concentration and crack initiation and/or propagation sites, since K TILSs represents the intensity of ILSs between contacting materials. According to these results, it appears that the location on interface between Al and Al oxide has the highest probability of crack initiation and/or propagation. On the other hand, it seems that crack initiation and/or propagation on interface between carbide and Al oxide is less likely to occur when compared to other contacting materials. Furthermore, it appears that the value of K TILSs for Al contacting with carbide is high. Hence, it seems that the second location of crack initiation and/or propagation is on the interface between Al and carbide.
Moreover, these results can show that if manufacturing process of carbon/Al composite can take place without developing aluminum oxide and carbide, the chance of crack initiation and/or propagation in carbon/Al composite exposed to thermal stress decreases. e reason is, based on the results obtained from this study, the first location with highest probability of stress concentration is Al/Al oxide interface. Furthermore, according to the results in Table 7, the second location with highest thermal stress intensity is on Al/carbide interface.
K TILSs of Al interacting with Al oxide and carbide are 2,984,800 N/k·m 2 and 2,773,500 N/k·m 2 , respectively. By eliminating of Al oxide and carbide in manufacturing process, the next interacting materials with highest K TILSs are CNT/Al with K TILSs equal to 728,000 N/k·m 2 . It means that, by eliminating the production of Al oxide and carbide in manufacturing process, the value of thermal interlaminar shear stress intensity factor decreases by 75.6% which causes high reduction in stress concentration and crack propagation on CNT/Al interface.

Summary and Concluding Remarks
(i) In the second part of the article, with the application of mathematical method, the values of ILSs max in common UFMCs due to thermal stress are obtained and indicated in Table 3. With determination of the values of ILSs max in UFMCs, these composites are ranked based on the probability of failure due to thermal cycling exposure in this environment. It means that the UFMC which has the minimum ILSs is the safest in case of failure due to the deterioration of ILSS. On the other hand, UFMC with the highest ILSs max has the highest chance of de-bonding, delamination, or fracture due to the exposure to thermal stress and/or thermal cycles. According to the mentioned method, the probability of fracture due to the exposure to thermal stress in unidirectional Sic fibers/borosilicate glass matrix composite is minimum because the ILSs max is the minimum value. Inversely, the probability of fracture in unidirectional diamond fibers/CE339 epoxy matrix composite is maximum. (ii) In Section 4 of this research, the differences between the CTEs of fibers and matrix within the SWCNTAF and MWCNTAF at temperature range of −5 to 85°C are calculated. Furthermore, the ILSs over the interfaces between the fibers and matrix within the SWCNTAF and MWCNTAF at the same temperature range are derived. According to the results obtained with this calculation, it appears that the chance of crack propagation rate in axial direction and interfaces between the fibers and matrix within the SWCNTAF is higher than that within the MWCNTAF, when they are exposed to thermal cycling. is conclusion is due to the higher stress concentration within SWCNTs and AF interfaces in comparison of that within the MWCNTs and AF interfaces.
(iii) In Section 5 of this study, with using analytical method and equation (4), the mismatch of CTEs and ILSs between fibers and matrix for both SWCNTAF and SWCNTE in temperature range of −5 to 85°C is derived. e results have shown that, for both SWCNTAF and SWCNTE, Δα max and ILSs max occur at 85 o C. Furthermore, the results have indicated that Δα max and ILSs max in SWCNTAF are 49% lower than those in SWCNTE. Additionally, Δα mean in SWCNTAF is 40% lower than that in SWCNTE. Moreover, ILSs mean in SWCNTAF is 44% smaller than that in SWCNTE. us, according to these results, it can be concluded that the probability of stress concentration, crack initiation and/or propagation, and de-bonding on interfaces between matrix and fibers in SWCNTAF is less than that in SWCNTE. Hence, thermal fatigue life of SWCNTAF could be higher than that in SWCNTE. According to this conclusion, application of SWCNTAF in many industries such as aerospace and automotive that are dealing with thermal cycling and/or thermal stress is recommended. (iv) e research topics related to carbon composites, Al, and CNT/Al are very broad, and in this study, only the most important topics have been covered. It seems that the topics related to the mechanical properties and damage mechanism are the most significant topics which this study also has tried to cover. (v) Based on the results which have been obtained by this research, ILSs between the SWCNT and Al foam is higher than that between MWCNT and Al foam. Consequently, stress concentration and thermal fatigue are higher within the SWCNTAF interface when compared to MWCNTAF interface. Furthermore, ILSs between the SWCNT and epoxy is higher than that between SWCNT and AF. erefore, thermal fatigue life of SWCNTAF could be higher than that in SWCNTE. Moreover, chemical reaction between carbon and molten aluminum can induce aluminum oxide and carbide on carbon/Al interface which can cause brittleness and stress concentration in this area. Furthermore, by new relation for thermal interlaminar shear stress intensity factor, it has been proven that interface of Al and Al oxide can induce the highest and interface of carbide and Al  Journal of Engineering 13 oxide can induce the lowest stress concentrations in the cases they are exposed to thermal stress or thermal cycling.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e author declares no conflicts of interest.