Evaluate the Fatigue Life of CFRC Subjected to Coupled Thermo–Mechanical Loading

Mechanical properties of composites manufactured by high-temperature polymer polyether ether ketone (PEEK) with continuous reinforced fibers are closely dependent on ambient temperature variations. In order to effectively study fatigue failure behaviors of composites under the coupled thermo–mechanical loading, a well-established microscopic model based on a representative volume element (RVE) is proposed in this paper. Stiffness degradation behaviors of the composite laminates at room and elevated temperatures are firstly investigated, and their failure strengths are compared with experimental data. To describe the fatigue behaviors of composites with respect to complex external loading and ambient temperature variations, a new fatigue equation is proposed. A good consistency between theoretical results and experimental data was found in the cases. On this basis, the temperature cycling effects on the service life of composites are also discussed. Microscopic stress distributions of the RVE are also discussed to reveal their fatigue failure mechanisms.


Introduction
Due to their excellent mechanical properties [1,2], such as fatigue resistance, good thermal stability, and fracture toughness [3,4], high-temperature polymer composites have been widely used in the aviation industry, and they are looking to replace metal materials in a variety of aircraft components, such as luggage racks and cable clips [5]. Airbus has discussed and intended to use high-temperature composite structures for wing structures and thicker fuselage structures [6].
In recent years, thousands of researchers have devoted themselves to deeply investigating the fatigue properties of composites by using experimental or theoretical methods [7][8][9][10]. Seyfullayev [7] and Li [8] established a model to investigate fatigue crack growth in composite structures. Carvalho et al. [8] studied the influence of reinforcement mechanisms of carbon nanotubes on wear, as well as fatigue tests on an aluminum-silicon composites. Huang et al. [10] employed the experimental method to study crack bridging in engineered cementitious composites under fatigue tensile loading. The fatigue properties have also been observed as the change in residual strength, cracks, external loading, stiffness, etc. [11][12][13][14][15]. Huang [11] and Zhao [12] established effective models to analyze the degradation properties and progressive fatigue damage in composite structures under fatigue tensile loading. At several temperatures, Eftekhari and Fatemi [13] investigated the cycling frequency on fatigue behaviors by conducting load-controlled fatigue tests. The results showed that fatigue behavior is also time-dependent for frequency-sensitive materials. Keulen et al. [14] presented a novel technique to predict the remaining life of composites under fatigue loading by employing the composite structures may be subjected to cycling temperature loading. Service life and microscopic failure investigations under coupled thermo-mechanical loading will lay a solid foundation for their usage and optimized structural design. Therefore, an effective microscopic model is urgently needed. To simplify the modelling procedure at the microscopic scale, herein, it was assumed that reinforced fibers are periodically distributed in matrix materials, and composites are composed of an infinite number of representative volume elements. On this basis, the RVE is further divided into N β * N γ sub-cells. According to the homogenization theory, the displacement vectors, u i (x,y), can also be expressed as the function of the average displacement term, u i (x, y), and their high-order terms, u (n) i (x, y)(n = 1, 2, . . .), that is: i (x, y) + δu (1) i (x, y) + O(δ 2 ), (1) where x = (x 2 , x 3 ) and y = (y 2 , y 3 ) are the macroscopic and microscopic coordinate systems, respectively. Due to its periodicity hypothesis of the RVE mentioned above, the boundary continuity condition of displacement terms, u (n) i (x,y), can be further expressed as: where k indicates the periodic dimension of the RVE at the microscopic scale.
It should be noted that the average displacement, u i (x, y), is unrelated to material parameters at the microscopic scale, and the component can be simplified as u i (x). Therefore, the displacement vectors, u i (x,y), in Equation (1) can be rewritten and expressed as follows: where u i = δ u (1) i . O(δ 2 ) indicates the higher-order terms of the fluctuating displacements. According to the elastic theory, average strain components, ε ij (x), and fluctuating strain components, ε ij (x, y), can easily be expressed as a function of displacement vectors, u i , and their high-order terms, u i , that is: Combining Equation (3) with Equations (4) and (5), the strain components, ε ij (x, y), in composites can be expressed as follows: Based on the strain relations shown in Equation (6), as well as the modelling scheme of the high-fidelity generalized method of cells (HFGMC), the relations between sub-cell displacements, , and average sub-cell displacements, W (βγ) , in a local coordinate system can be written as [34]: where h β and l γ are the sub-cell height and length in the RVE, respectively. W (βγ) (pq) indicates higher-order displacement terms. Similarly, temperature components in composites can be asymptotically expanded as the function of the small scaling parameter, δ, when the steady heat transfer problem is considered, that is: where T and T (0) indicate temperature and average temperature components, respectively. T (n) (n = 1, 2, · · · ) indicates the high-order terms of temperature components.
According to Fourier's law, the relations between sub-cell heat flux vectors, q (βγ) i , and thermal conductivity coefficients, k ij , can be written as follows: where H (βγ) i is related to the sub-cell temperature gradient. For investigating the thermal conductivity problem of the CFRC under steady-state conditions, the local equilibrium equation in the RVE can be expressed as: Similar to solving procedures of the sub-cell displacements, sub-cell temperature components, T (βγ) , will be expressed as a function of the heat flux vectors, q (βγ) i , and thermal conductivity coefficients, ij . To solve the coupled thermo-mechanical problems, sub-cell displacement components, W (βγ) , were further extended as a function of the sub-cell displacements and temperatures, that is: In the procedure of thermo-mechanical modelling, the constitutive relation of the CFRC is employed and the sub-cell average stress, σ . In order to simplify the calculation and improve efficiency, the surface average displacements, tractions, and heat flux in each sub-cell were employed. In addition, the continuity between adjacent sub-cell displacements and tractions was employed to solve the coupled thermo-mechanical equations. Detailed modelling procedures were reported by Aboudi and Ye [34,35]. Finally, the macroscopic average component, Y = [σ 11 σ 22 σ 33 σ 23 σ 13 σ 12 q 1 q 2 q 3 ] T , can be expressed as a function of the macroscopic average strains, ε ij (i, j = 1, 2, 3), and temperature gradients, H i (i = 1, 2, 3), according to the homogenization theory, that is:Ȳ where X = [ε 11 ε 22 ε 33 ε 23 ε 13 ε 12 − H 1 − H 2 − H 3 ]. A (βγ) indicates a 9 × 9 matrix.

Specimen Preparations
In general, material properties are closely dependent on the ambient temperature. It means that their parameters, such as elastic constants and effective strength, are temperature dependent. To analyze the effective behaviors and fatigue properties of AS4/PEEK composites (Xiechuang Composite Materials Co. LTD, Dongguan, China) subjected to coupled thermo-mechanical loading, constituent material parameters were considered to be a function of the ambient temperature. The specimens were prepared by stacking and cutting the composite laminates, and their dimensions were determined according to the standard of ASTM D3039/D3039M [36] and STP1330-98 [37]. The manufacturing procedures of composite specimens, which are shown in Figure 1, can be summarized as the following steps: (1) Cut and lay the carbon fiber prepreg according to the desired sequences of composite laminates.
(2) Clean mold and apply the release agent uniformly on it. In this procedure, we should try to keep the mold surface as smooth as possible. constituent material parameters were considered to be a function of the ambient temperature. The specimens were prepared by stacking and cutting the composite laminates, and their dimensions were determined according to the standard of ASTM D3039/D3039M [36] and STP1330-98 [37]. The manufacturing procedures of composite specimens, which are shown in Figure 1, can be summarized as the following steps: (1) Cut and lay the carbon fiber prepreg according to the desired sequences of composite laminates.
(2) Clean mold and apply the release agent uniformly on it. In this procedure, we should try to keep the mold surface as smooth as possible.   Table 1 indicates the matrix material properties of the Polyetheretherketone (PEEK, Victrex plc, Lancashire, UK) at room temperature (23 °C ) and elevated temperatures. It is clearly indicated that material properties are temperature dependent. A higher ambient temperature tends to further reduce the elastic modulus and effective strength. The glass transition temperature of the PEEK resin  Table 1 indicates the matrix material properties of the Polyetheretherketone (PEEK, Victrex plc, Lancashire, UK) at room temperature (23 • C) and elevated temperatures. It is clearly indicated that material properties are temperature dependent. A higher ambient temperature tends to further reduce the elastic modulus and effective strength. The glass transition temperature of the PEEK resin is 143 • C. Composites with continuous fibers AS4 (Hercules company, Wilmington, DE, USA), which afford excellent high-temperature properties, have good chemical resistance.

Tensile Test
In the tensile test procedure, the test system was composed of Material Test System (MTS, Instron 8802J5219, 250NK, Norwood, MA, USA) and the PC is shown in Figure 2a,b. The temperature chamber (Instron CP100358, maximum temperature 350 • C, Norwood, MA, USA) was adopted to control and keep the ambient temperature at specified levels, i.e., room temperature and 100 • C, during the test. Before tensile tests at the elevated temperature, it was kept for 3 min when the temperature increased to 100 • C in order to control the uniform distribution over the cross-section of the specimens. For each case, three specimens were employed. The stress-strain behaviors of the PEEK matrix under a strain rate of 2 × 10 −5 are seen in Figure 2c. The experiment results indicate that the PEEK matrix at room and high temperatures presents high ductility. However, a higher ambient temperature tends to sharply decrease the loading capacity. In detail, the longitudinal tensile strength of the PEEK matrix was 94.58 and 50.37 MPa when room temperature and 100 • C were considered, respectively.
A proper RVE plays an important role in studying the mechanical properties of composite structures. In order to verify the effectiveness of the microscopic model, two typical ply methods of [0/90] 4s and [0 composite laminates were considered in the cases. The material parameters of AS4 fibers, such as elastic moduli, Poisson's ratio, and failure strength, were provided by Soden et al. [38][39][40][41]. It can be seen that the PEEK matrix was considered to be isotropic. In this material system, the fiber volume fraction (VF), which was determined by using the micrograph of the fiber cross-section, is approximate to 0.60. It should be pointed out that the shear strength of PEEK at 100 • C was acquired by employing the linear interpolation method based on the failure properties of the PEEK at room temperature [41] and a high temperature of 90 • C derived from the experimental methods [42]. Table 2 indicates a comparison of the longitudinal moduli of composite laminates between the experimental data and theoretical results at both room temperature and high temperatures. In a word, a high prediction accuracy of the effective moduli can be obtained. The maximum relative error is approximate to 2.90% when composite laminates [0  To further confirm the rationality of this chosen RVE, the stiffness degeneration analysis and failure strength comparisons of composite laminates were also examined. To describe constituent material damages at room and elevated temperatures, the 3D Tsai-Hill and Maximum stress criteria [43,44] were developed and proposed to study microscopic damages in the matrix and fibers, respectively. Figure 3 indicates the stress-strain of the lamina with 0 • -reinforced fibers at room temperature. A linear stress-strain behavior can easily be found until the longitudinal strain approximates to 1.5%. Moreover, chamber (Instron CP100358, maximum temperature 350 °C , Norwood, MA, USA) was adopted to control and keep the ambient temperature at specified levels, i.e., room temperature and 100 °C , during the test. Before tensile tests at the elevated temperature, it was kept for 3 min when the temperature increased to 100℃ in order to control the uniform distribution over the cross-section of the specimens. For each case, three specimens were employed. The stress-strain behaviors of the PEEK matrix under a strain rate of 5 2 10   are seen in Figure 2c. The experiment results indicate that the PEEK matrix at room and high temperatures presents high ductility. However, a higher ambient temperature tends to sharply decrease the loading capacity. In detail, the longitudinal tensile strength of the PEEK matrix was 94.58 and 50.37 MPa when room temperature and 100 °C were considered, respectively.  [38][39][40][41]. It can be seen that the PEEK matrix was considered to be isotropic. In this material system, the fiber volume fraction (VF), which was determined by using the micrograph of the fiber cross-section, is approximate to 0.60. It should be pointed out that the shear strength of PEEK at 100 °C was acquired by employing the linear interpolation method based on the failure properties of the PEEK at room temperature [41] and a high temperature of 90 °C derived from the experimental methods [42]. Table 2 indicates a comparison of the longitudinal moduli of composite laminates between the experimental data and theoretical results at both room temperature and high temperatures. In a word, a high prediction accuracy of the effective moduli can be obtained. The maximum relative error is approximate to 2.90% when composite laminates   Based on validations of the micromechanical method in analyzing the lamina, herein, stress-strain properties of [0/90] 4s and [0/45/90/ − 45] 4s composite laminates were also considered. To this end, the implemented procedures were executed by employing the microscopic model at each integration point of the lamina along the thickness direction [40]. The longitudinal tensile strength of the fibers was assumed to be 2900 MPa at 100 • C. The stress-strain results of composite laminates at room and elevated temperatures are shown in Figure 4. It can be seen that the laminates [0/90] 4s always provide higher stiffness behaviors and final failure strengths than the laminates [0/45/90/ − 45] 4s whether the ambient temperature is considered or not. Compared with the unidirectional (UD) lamina as show in Figure 3, obvious differences in the stress-strain behaviors were found when the layer sequence was considered in the example. In detail, obvious turning points, F and M, in the stress-strain curves of composite laminates at room temperature can easily be found as shown Figure 4a. Axial tensile loading tends to linearly increase their stress-strain behaviors at the initial stage. However, a sharp downward phenomenon was found when longitudinal strains increased to 0.76% for the [

Proposed Fatigue Damage Model
The stress amplitude versus number of cycles (S-N) curve is one of the most popular methods to describe fatigue life in the engineering field. In order to effectively predict the service life of composite structures, a new set of fatigue equations based on the computational micromechanics' RVE is proposed. Here, AS4 fibers were considered to be linearly degenerated with respect to the cycling number. Matrix fatigue can be described by Chaboch's fatigue equation [46][47][48] as shown in Table 4.

Proposed Fatigue Damage Model
The stress amplitude versus number of cycles (S-N) curve is one of the most popular methods to describe fatigue life in the engineering field. In order to effectively predict the service life of composite structures, a new set of fatigue equations based on the computational micromechanics' RVE is proposed. Here, AS4 fibers were considered to be linearly degenerated with respect to the cycling number. Matrix fatigue can be described by Chaboch's fatigue equation [46][47][48] as shown in Table 4. It indicates that the fatigue life, N f , is a function of the maximum stress amplitude, σ max , and average stress amplitude, σ. σ u and σ f l indicate the maximum stress and fatigue limit of the polymer matrix. The fatigue parameters, a, β, and, will be confirmed according to the S-N curve of the PEEK matrix and composites. It can easily be found that this fatigue equation is restricted to the loading condition. If fatigue tests were executed with a stress ratio of σ min /σ max = 0.3, the denominator, σ max − σ − σ f l , may be a negative number. To extend its application range, a new fatigue pertinent equation is proposed as shown in Table 4. Here, the parameter, w, is defined as a fatigue control parameter, which can be acquired according to the fatigue experimental data of the matrix and composite materials.

Fatigue Model Validation
In the fatigue test, fatigue properties of constituent materials were closely related to the loading number and amplitude. Based on the tested results by Zhu [16], the fatigue limit of 1300 MPa of AS4 fibers was considered in the example. According to the fatigue test S-N curves of the PEEK matrix at room temperature, as well as AS4/PEEK composites at room and elevated temperatures executed by McKeen et al. [49], the fatigue limit of 61.21 and 21.47 MPa of the PEEK matrix at room temperature and 100 • C were employed in predicting the service life of composites, respectively. Here, the MTS was employed to execute T-T fatigue loading, and the environment chamber was employed to execute fatigue tests at 100 • C as a further comparison. A constant stress amplitude, by controlling the load mode with the stress ratio of σ min /σ max = 0.1, was employed, and the frequency was 2 Hz. According to the experimental curves of the PEEK and composites at room temperature [49] as shown in Figure 5a,b, the fatigue parameters in the proposed fatigue equation can be confirmed as follows: β = 6, a = 0.06, M = 145 MPa and w = 0.01. It can easily be seen that the theoretical prediction shows good consistency with the experimental data at room temperature, and a higher cycling loading will sharply decrease their fatigue endurance. It is interesting to mention that similar conclusions were found at 100 • C as shown in Figure 5c when the identical fatigue parameters mentioned above were employed in the example. In addition, the raised ambient temperature tended to reduce their service life to a great degree.
To further illuminate the VF influences, a different volume fraction of reinforced fibers of 0.5 was also considered in the example. From the figure, it is shown that the higher VF tended to increase their service life at both room and elevated temperatures. It means that AS4 fibers play an important role in increasing the fatigue life and failure strength for composite structures. In detail, the fatigue limits were 578.59 and 342.12 MPa when the VF of reinforced fibers was 0.6 and 0.5, respectively.
to reduce their service life to a great degree.
To further illuminate the VF influences, a different volume fraction of reinforced fibers of 0.5 was also considered in the example. From the figure, it is shown that the higher VF tended to increase their service life at both room and elevated temperatures. It means that AS4 fibers play an important role in increasing the fatigue life and failure strength for composite structures. In detail, the fatigue limits were 578.59 and 342.12 MPa when the VF of reinforced fibers was 0.6 and 0.5, respectively.

Evaluations of Fatigue Life under Coupled Thermo-Mechanical Loading
The material parameters are always temperature dependent. To reveal coupled thermo-mechanical effects on stiffness degeneration and fatigue failure properties, a small temperature variation between 13 and 33 • C was investigated. In the study, a temperature cycling frequency of 0.05 Hz was considered as shown in Figure 6a. For comparisons, two T-C mechanical loadings of 0.5 and 2.5 Hz were considered as shown in Figure 6b. Figure 7a indicates the stress amplitude of [0/45/90/ − 45] 4s composite laminates with respect to the cycling number in coupled thermo-mechanical loading. It can easily be seen that temperature cycling coupled with a mechanical loading of 2.5 Hz tends to sharply decrease the fatigue life of composite structures. However, a different conclusion was obtained when a mechanical frequency of 0.5 Hz was considered in the example. In detail, it is interesting to mention that coupled temperature cycling tends to firstly increase and then decrease the fatigue life with a decrease of the stress amplitude. To further clearly discern this phenomenon, a local enlargement graph was employed as shown in Figure 7b. In detail, the temperature cycling tends to extend the fatigue life of composite structures when the stress amplitude derived from mechanical loading is higher than 686.3 MPa. However, an opposite conclusion was acquired once a small mechanical amplitude was considered. That is, the temperature cycling tends to decrease the fatigue life when a small stress amplitude was considered. In addition, the cycling frequency of mechanical loading effects on the fatigue life of composite laminates under the thermo-mechanical loading can be ignored when the stress amplitude is less than 529.9 MPa.
shown in Figure 7b. In detail, the temperature cycling tends to extend the fatigue life of composite structures when the stress amplitude derived from mechanical loading is higher than 686.3 MPa. However, an opposite conclusion was acquired once a small mechanical amplitude was considered. That is, the temperature cycling tends to decrease the fatigue life when a small stress amplitude was considered. In addition, the cycling frequency of mechanical loading effects on the fatigue life of composite laminates under the thermo-mechanical loading can be ignored when the stress amplitude is less than 529.9 MPa.

Local Stress Distribution under Coupled Thermo-Mechanical Loading
To reveal the fatigue failure mechanisms under coupled thermo-mechanical loading, it is critical to accurately capture stress distributions at the microscopic scale. Based on the failure and fatigue validations mentioned above, herein, local stress fields in the RVE were employed. In the following cases, it should be noted that the frequency ratio of 0.1 of temperature cycling vs. mechanical cycling was considered. Meanwhile, the mechanical loading ratio of min max /  was −1. Figure 8 indicates the longitudinal stress distributions due to a rising temperature from 25 to 31 °C . It can easily be found that a higher temperature tends to increase the maximum stress, 11  , at the microscopic scale, and the thermal stress exhibits compression in the matrix and tension in fibers due to the mismatch of thermal expansion coefficients between reinforced fibers and the matrix. In addition, the maximum compressive stress was found in the matrix materials. Figure 9 indicates the microscopic stress distributions of 22  , 33  , and 23  subjected to a temperature loading at 33 °C . It was observed that the maximum tensile and shear stresses are both located at the interface between the fibers and matrix materials.

Local Stress Distribution under Coupled Thermo-Mechanical Loading
To reveal the fatigue failure mechanisms under coupled thermo-mechanical loading, it is critical to accurately capture stress distributions at the microscopic scale. Based on the failure and fatigue validations mentioned above, herein, local stress fields in the RVE were employed. In the following cases, it should be noted that the frequency ratio of 0.1 of temperature cycling vs. mechanical cycling was considered. Meanwhile, the mechanical loading ratio of σ min /σ max was −1. Figure 8 indicates the longitudinal stress distributions due to a rising temperature from 25 to 31 • C. It can easily be found that a higher temperature tends to increase the maximum stress, σ 11 , at the microscopic scale, and the thermal stress exhibits compression in the matrix and tension in fibers due to the mismatch of thermal expansion coefficients between reinforced fibers and the matrix. In addition, the maximum compressive stress was found in the matrix materials. Figure 9 indicates the microscopic stress distributions of σ 22 ,σ 33 , and σ 23 subjected to a temperature loading at 33 • C. It was observed that the maximum tensile and shear stresses are both located at the interface between the fibers and matrix materials. microscopic scale, and the thermal stress exhibits compression in the matrix and tension in fibers due to the mismatch of thermal expansion coefficients between reinforced fibers and the matrix. In addition, the maximum compressive stress was found in the matrix materials. Figure 9 indicates the microscopic stress distributions of 22  , 33  , and 23  subjected to a temperature loading at 33 °C . It was observed that the maximum tensile and shear stresses are both located at the interface between the fibers and matrix materials.   It can be seen from the figures that the maximum stresses are located at the interface between the matrix materials and fibers, and a distinct difference of stress distributions can easily be found. In detail, the reinforced fibers, which support the tensile or compressive loadings, are decided by the external loading status. To quantitatively explain coupled cycling temperature effects on their effective properties, cycling temperatures at 0.5 s and 1.5 s were also considered as shown in Figure 10c,d. Similarly, it can also be found that the maximum stress at the interface region is higher than that in the fiber or matrix region. In other words, this finite cycling temperature is insufficient to alter the stress distributions and maximum stress region. In addition, it was revealed that the maximum tensile stress of the composites under mechanical loading at 0.5 s was 74.849 MPa as shown in Figure 10a. While a higher stress of 75.129 MPa was found at 0.5 s when the cycling temperature was considered as shown in Figure 10c. This slight increase in maximum stress may be increase the probability of the initial damages at the interface with composites.  Figure 10a,b indicates the microscopic stress distributions, σ 22 , of composites under the pure mechanical loadings at 0.5 and 1.5 s, respectively. It can be seen from the figures that the maximum stresses are located at the interface between the matrix materials and fibers, and a distinct difference of stress distributions can easily be found. In detail, the reinforced fibers, which support the tensile or compressive loadings, are decided by the external loading status. To quantitatively explain coupled cycling temperature effects on their effective properties, cycling temperatures at 0.5 s and 1.5 s were also considered as shown in Figure 10c,d. Similarly, it can also be found that the maximum stress at the interface region is higher than that in the fiber or matrix region. In other words, this finite cycling temperature is insufficient to alter the stress distributions and maximum stress region. In addition, it was revealed that the maximum tensile stress of the composites under mechanical loading at 0.5 s was 74.849 MPa as shown in Figure 10a. While a higher stress of 75.129 MPa was found at 0.5 s when the cycling temperature was considered as shown in Figure 10c. This slight increase in maximum stress may be increase the probability of the initial damages at the interface with composites. Figure 11a,b indicates the microscopic stress distributions, σ 23 , under the mechanical loading at 0.5 and 1.5 s, respectively. Similar to the stress distributions, σ 22 , mentioned above, it can be seen that the maximum shear stress, σ 23 , is also located at the interface region, and the temperature cycling will not affect the microscopic stress distributions as shown in Figure 11c,d. In addition, it can be found that a higher temperature tends to increase the maximum shear stress in composites. In detail, when a loading time of 0.5 s was considered, the maximum shear stresses were 11.174 and 11.218 MPa under mechanical loading and coupled thermo-mechanical loading, respectively. coupled cycling temperature effects on their effective properties, cycling temperatures at 0.5 s and 1.5 s were also considered as shown in Figure 10c,d. Similarly, it can also be found that the maximum stress at the interface region is higher than that in the fiber or matrix region. In other words, this finite cycling temperature is insufficient to alter the stress distributions and maximum stress region. In addition, it was revealed that the maximum tensile stress of the composites under mechanical loading at 0.5 s was 74.849 MPa as shown in Figure 10a. While a higher stress of 75.129 MPa was found at 0.5 s when the cycling temperature was considered as shown in Figure 10c. This slight increase in maximum stress may be increase the probability of the initial damages at the interface with composites.   , mentioned above, it can be seen that the maximum shear stress, 23  , is also located at the interface region, and the temperature cycling will not affect the microscopic stress distributions as shown in Figure 11c,d. In addition, it can be found that a higher temperature tends to increase the maximum shear stress in composites. In detail, when a loading time of 0.5 s was considered, the maximum shear stresses were 11.174 and 11.218 MPa under mechanical loading and coupled thermo-mechanical loading, respectively.

Conclusions
This study focused on the fatigue life of composite structures under coupled thermo-mechanical loading. To validate the proposed microscopic model and the chosen RVE, the prediction results of effective moduli and failure properties of composites at room and elevated temperatures were both investigated and compared with experimental data. It was shown that the maximum relative error of the failure strength was equal to 5.51% when the composite laminate 4s [0 / 90] was considered at an ambient temperature of 100 °C . In order to extend the application range of the fatigue model, the

Conclusions
This study focused on the fatigue life of composite structures under coupled thermo-mechanical loading. To validate the proposed microscopic model and the chosen RVE, the prediction results of effective moduli and failure properties of composites at room and elevated temperatures were both investigated and compared with experimental data. It was shown that the maximum relative error of the failure strength was equal to 5.51% when the composite laminate [0/90] 4s was considered at an ambient temperature of 100 • C. In order to extend the application range of the fatigue model, the control parameter, u, was introduced. Compared with experimental data at room temperature and 100 • C, it was shown that the proposed fatigue equation can be effectively employed to describe the fatigue life. A higher VF tends to sharply increase their service life at both room and elevated temperatures.
It is interesting to mention that thermal cycling will deeply affect their service life. When a mechanical frequency of 0.5 Hz was considered, the coupled temperature cycling tended to first increase and then decrease the fatigue life of composites with the decreasing of stress amplitude. However, the finite cycling temperature studied in this paper was insufficient to alter the stress distributions, and the maximum stress was located at the interphase between reinforced fibers and matrix materials. The proposed method in this paper is also suitable for on-line life evaluation of composite structures of whether the loads can be measured exactly by the strain sensors. However, a higher accuracy of the microscopic stress field in the RVE can be acquired when parameter elements are considered as discrete to the RVE at the microscopic scale. Therefore, our research group is trying to construct parameter elements to achieve an effective life evaluation with a higher accuracy. In addition, the SEM will be used to further reveal failure mechanisms of composites at the elevated temperature.