Abstract
As aviation and space technology advances rapidly, fiber-reinforced polymer (FRP) composites are finding growing applications under extreme service conditions, particularly under high-thermal environments. Assessing the fatigue resistance of FRP composites across different temperature ranges has become a major concern. By applying the Force-Heat Equivalent Energy Density (FHEED) principle and damage failure analysis, a temperature and damage dependent fatigue strength model (TDDFS) was established for FRP composites, grounded in physical mechanisms. The TDDFS model considers the comprehensive influence of temperatures, evolution of constituent properties, load level, layup configurations, especially the damage evolution with temperature on the fatigue strength. The well comparisons between model predictions and experimental data demonstrate the reliability and reasonability of TDDFS model. The study also analyzed the fatigue performance and temperature-related damage evolution of FRP composites. These insights contribute to a reliable theoretical approach for TDDFS prediction and deepen the understanding of fatigue degradation and damage mechanisms in FRP materials at elevated temperatures.
Highlights
(1) The physics-based temperature dependent fatigue strength model of FRP composites is developed. (2) The model considers the effects of temperature, constituents’ properties and damage evolution. (3) The TDDFS of FRP composites can be well predicted by the model. (4) This study provides the understanding on fatigue analysis and strength prediction method.
Introduction
Fiber-reinforced polymer (FRP) composites have obtained significant attention in aerospace, defense, and transportation industries due to their exceptional properties, including high specific mechanical properties, strong designability, and superior processability.1–9 The mechanical properties and fatigue behavior of FRP composites are predominantly governed by the type of reinforcing fibers, while the toughness of the matrix resin significantly influences their fatigue crack propagation resistance and damage tolerance characteristics.10–14 With the escalating service temperature requirements in these advanced engineering applications, investigating the thermomechanical properties of FRP composites has evolved into a crucial scientific challenge.15–18 Notably, fatigue strength, a pivotal indicator of material durability, exhibits pronounced sensitivity to thermal conditions and progressive damage accumulation. The comprehensive investigations into the fatigue response at different temperatures and damage conditions are significant for evaluating the durability of FRP composites. Particularly, establishing a robust quantitative theoretical framework to characterize the synergistic effects of thermal degradation and microstructural damage on fatigue strength has become a fundamental scientific challenge, driving substantial research in materials science and mechanical engineering.19–21
Extensive research has been devoted to understanding the thermomechanical characteristics and damage evolution in FRP composites. Monte et al. 22 and Sonsino et al. 23 systematically studied the axial fatigue performance of short glass fiber-reinforced polyamide 6.6 across temperature gradients. Their experimental findings revealed a temperature-induced degradation of fatigue strength while demonstrating near-constant S-N curve slopes regardless of thermal conditions. Loos et al. 24 investigated the fatigue strength of carbon nanotube/polyurethane composites, which demonstrated significant enhancement in fatigue resistance across both high-cycle and low-cycle regimes upon CNT incorporation. Ren et al.25,26 performed the research regarding the fatigue damage modes of single-walled carbon nanotubes reinforced epoxy (CNT/epoxy) composites, and used S-N curves to describe the fatigue behavior of nanocomposites, which suggests that CNT bridging serves as the primary fatigue resistance mechanism. Guo et al.27,28 proposed a stress-normalized method for fatigue life prediction that accounts for the effect of temperature on fatigue damage, and enables life prediction under high-temperature and variable-amplitude loading. Furthermore, by considering the microstructure and ply thickness of multidirectional laminates, they developed a fatigue life prediction model suitable for high-temperature conditions where compression-dominated fatigue failure prevails. Ramakrishnan et al. 29 proposed a multi-scale stiffness-based damage model that quantifies composites degradation through constituents failure mechanisms, and establishing explicit correlations between damage accumulation and mechanical property deterioration. These investigations have significantly advanced the understanding of the fatigue properties evolution in FRP composites under combined thermal and mechanical damage, thereby establishing a critical foundation for physics-based modeling frameworks. Moreover, there is an urgent need to develop a temperature-dependent fatigue strength prediction model capable of estimating thermal fatigue performance based on fundamental material parameters. However, under coupled thermomechanical loading, the fatigue damage progression follows a sequential mechanism involving matrix cracking, crack bridging, interfacial failure between fiber and matrix, and ultimately fiber breakage. Currently, the existing models fail to account for the synergistic effects of real-time damage states and their temperature-dependent evolution on the fatigue resistance of FRP composites.
Consequently, the present investigation focuses on characterizing the fatigue performance and damage evolution of FRP composites under thermal variations, with the goal of developing a mechanistic model for fatigue strength prediction under multi-stage damage accumulation at elevated temperatures. To this end, we propose a TDDFS model by integrating the FHEED principle with established damage progression laws. In particular, this study innovatively introduces the FHEED principle into the field of fatigue failure of FRP composites at different temperatures. By establishing the equivalent relationship between heat energy and strain energy, a novel TDDFS theoretical model can be developed. The model systematically incorporates the synergistic effects of temperature, cyclic loading parameters, layup configurations, and crucially, the temperature-dependent damage evolution. The validation was conducted through experimental data and predictions from the TDDFS model, complemented by parametric analyses to quantify the sensitivity of fatigue properties to critical damage and temperature variables. This model can quantitatively characterize the entire process of fatigue damage evolution without relying on any empirical fitting parameters. Moreover, compared with the empirical model proposed by Jen et al.,30,31 the presented model has a stronger physical foundation and greater universality. This study establishes a robust theoretical framework for fatigue strength prediction and damage tolerance assessment of FRP composites in elevated temperature environments.
Theoretical derivations
The damage evolution in FRP composites under varying temperatures involves a progressive and intricate process, characterized by the cumulative development of micro-scale damage mechanisms, including interfacial debonding, matrix microcracking, and fiber breakage.21,32 Exploring the fatigue damage behavior and theoretical characterizing the temperature-dependent fatigue strength of FRP composites under thermo-mechanical fatigue load are the key issue for safety use and reliability evaluation.
Extending the FHEED principle established in our prior research,33–36 a novel theoretical framework was developed to quantitatively characterize the temperature-dependent fatigue behavior of FRP composites, representing the innovative theoretical model of its kind. According to this principle, the underlying assumption is that a critical energy storage value exists for a given FRP composites, i.e. the fatigue failure of composites is associated with a fixed maximum energy storage value. Meanwhile, since the failure of FRP composites is mainly a result of the combined effect of strain energy and heat energy, the critical energy storage value consists of the critical strain energy and heat energy of the composites. In order to clarify and deepen the modeling understanding on the fatigue strength, the schematic illustration of energy equivalence principle of FRP composites under fatigue load is depicted in Figure 1. Therefore, the governing equation can be derived as follows: Schematic illustration of energy equivalence principle and damage accumulation of FRP composites under fatigue load.
It is assumed that the density of the composite material is independent of temperature due to the weak temperature sensitivity. The heat energy density of composites at temperature T is given by:
The influence of temperature on the critical strain energy density at failure of FRP composites, denoted by
A primary objective is to derive the expressions for the critical strain energy density of the fiber and matrix components. When the service temperature is lower than the glass transition temperature, it can be assumed that the matrix approximately exhibits linear elastic behavior within the analyzed temperature range. Building on prior work, we derive the relationship between the thermally influenced internal stress distribution in constituent materials and the external load applied to the composite, yielding the following formulation
38
:
To account for the damage mechanisms including interfacial degradation, matrix cracking, and fiber breakage in fatigue strength prediction, the damage parameters were integrated into the average stress equations of constituents. Building on the mechanics of continuous damage, the intrinsic damage in composites model
39
was further generalized to accommodate temperature-dependent behavior:
Assuming a fiber direction of
Equation (5) and Equation (6) denote the mean stresses in the fiber and matrix of FRP composites under external load without considering damage effects. Combining the definition of damage with the introduction of corresponding damage parameters, the temperature-dependent mean stresses of the fiber (
The mean stress formulas for fiber and matrix varying with temperature and damage were established. The fiber and matrix can be considered to have reached the ultimate bearing state when the ultimate fatigue load (
The critical strain energy and heat energy versus temperature for FRP composites taking into account the damage effects are given, combining equations (1), (2), and (4), and
Based on equations (14) and (15), we obtain the ratio coefficients
Subsequently, combing equations (1), (4), (14), and (16), at any temperature
Substituting equations (8), (12), and (13) into equation (17) to derive
The TDDFS model (equation (18)) provides the theoretical framework for predicting the fatigue strength of unidirectional FRP composites at different temperatures and different damage levels. The model incorporates the synergistic influence of thermal and mechanical loading conditions, with particular emphasis on temperature dependent damage evolution and its impact on the fatigue performance of FRP composites. It should be noted that the present model is established under the premise of a constant ambient temperature T, and is applicable for predicting the fatigue strength of FRP composites under this condition. The potential influence of material self-heating during cyclic loading was not considered in the model. Moreover, the model can predict the temperature-dependent fatigue strength of composites under constant load but varying fiber content, requiring only a limited set of commonly available material parameters from standard handbooks or literature. The detailed analysis regarding the various factors affecting the fatigue properties of FRP composites can be conducted by employing the model, and it contributes to enhanced comprehension on the fatigue strength degradation and damage evolution of FRP composites at elevated temperatures.
Considering the relationship between damage and layup, a temperature dependent damage parameter is introduced
Specifically, the rule of mixtures can be applied to orthogonal woven composites, resulting in:
Note that
Substituting equations (12), (13), and (19) into equation (17) yields
The TDDFS model (equation (21)) was developed to predict the fatigue strength of orthogonal woven composites under varying temperatures and progressive damage conditions. A fundamental assumption of this model is that the composites are at a constant ambient temperature T, and the model does not incorporate self-heating effects during cyclic loading. Based on this premise, this model explicitly considers the combined effects of temperature, load level, especially the damage evolution with temperature on the fatigue strength. This study establishes a theoretical framework to elucidate the temperature dependent damage progression and its dominant role in composite fatigue degradation. The proposed model utilizes standard material properties readily available in technical databases, ensuring broad applicability. In this study, the TDDFS model supplies a theoretical prediction method for fatigue strength at different temperatures, and enables a systematic analysis of significant influencing factors about fatigue properties evolution in orthogonal woven composites, thereby advancing the understanding of their strength degradation and damage accumulation in different temperature environments.
Results and discussion
Fatigue strength prediction of polymer composites at different temperatures
The TDDFS model (equation (18)) of FRP composites, developed in the above section, aims to predict the fatigue strength at different temperatures and multi-factor conditions. To verify the scientific validity and precision of the developed model, this section compares the experimental fatigue strength with the predictions of TDDFS model across varying fiber contents, load levels, and temperatures for systematically assessing their agreement.
CNTs-GNP and epoxy composites with different fiber contents and load levels
Material parameters of 10% CNTs/epoxy, 10%GNP/epoxy and 9% GNP-1%CNTs/epoxy.

TDDFS model predictions and experimental data for 10% CNTs/epoxy at (a) 60% load level, (b) 62.5% load level, (c) 65% load level.

TDDFS model predictions and experimental data for 10% GNP/epoxy at (a) 57.5% load level, (b) 60% load level, (c) 65% load level.

TDDFS model predictions and experimental data for 9% GNP-1%CNTs/epoxy at (a) 60% load level, (b) 65% load level, (c) 70% load level.
CNTs and epoxy composites with different fiber contents and load levels
Material parameters of 0.5% CNTs/epoxy and 1% CNTs/epoxy.

TDDFS model predictions and experimental data for 0.5% CNTs/epoxy at (a) 40% load level, (b) 45% load level, (c) 50% load level, and (d) 55% load level.

TDDFS model predictions and experimental data for 1% CNTs/epoxy at (a) 35% load level, (b) 40% load level, (c) 45% load level, and (d) 50% load level.
AS-4 fiber and PEEK composites with different layups and loading levels
Material parameters of AS-4/PEEK composites (from ICI Fiberite Co. (USA)).

TDDFS model predictions and experimental data for AS-4/PEEK at (a) Cross-ply [0/90]4S, (b) Quasi-isotropic [0/+ 45/90/−45]2S.
Analysis of influencing factors of FRP composite materials
This study develops a TDDFS theoretical model based on the FHEED principle, which establishes the equivalence between strain energy and heat energy. By incorporating a damage parameter into the critical strain energy equations for both fiber and matrix, the model accounts for damage effects, including matrix microcracking and fiber breakage, on fatigue strength. Furthermore, this section analyzes key factors influencing the fatigue strength of FRP composites and their temperature-dependent variations.
The influence of temperature on the degree of damage and fatigue strength
Building upon the proposed TDDFS model, we systematically examined temperature effects on both damage parameters and fatigue strength in FRP composites, employing AS-4/PEEK cross-ply laminates as a representative case study. As shown in Figure 8, the fatigue strength exhibits a gradual decline with rising temperature, whereas the damage parameters increase.
29
This trend can be attributed to the temperature-induced degradation of material properties, which reduces the load-bearing capacity of composites. Meanwhile, the high-temperature properties degradation of fiber and matrix promotes premature failure in the composites, reflecting an elevated damage state. Consequently, temperature elevation directly diminishes the fatigue strength of the composites while simultaneously amplifying their damage degree, and then a dual mechanism that synergistically accelerates the fatigue strength reduction. Temperature-dependent evolution of damage accumulation and fatigue resistance.
The influence of interface bonding performance on damage at different temperatures
The cross-sectional analysis characterizes the relationship between fiber/matrix bonding properties and damage degree of composites at varying temperatures, as shown in Figure 9. The results demonstrate that at a fixed temperature, the damage parameters escalate as the decrease of interfacial parameter f *, which can be interpreted that the fiber/matrix bonding transitions from strong to weak caused the increasingly damage degree of FRP composites. Moreover, under the identical bonding condition, the damage parameters rise with increasing temperature. These analysis results suggest that weakened fiber/matrix bonding and elevated temperatures exacerbates microscopic damage mechanisms, including matrix cracking, fiber fracture, and interfacial defects. In particular, at elevated temperatures, FRP composites with inferior fiber/matrix bonding will achieve the significant damage degree, which further leads to the decreasing fatigue strength. Temperature dependent variation of damage level with fiber/matrix bonding properties.
The influence of young's modulus at different temperatures on fatigue strength
Figure 10 demonstrates that fatigue strength exhibits distinct temperature-dependent correlations: a positive relationship with matrix Young’s modulus versus a negative relationship with fiber Young’s modulus. This contrast is particularly pronounced in cross-ply FRP composites, where fiber modulus typically exceeds matrix modulus by an order of magnitude. Continuously increasing the fiber Young’s modulus causes its high-stiffness, and exacerbate the stress concentration. Further, It will reduce the interlaminar shear strength, and accelerate the crack propagation while lowering the fracture toughness, which leads to the decrease of temperature dependent fatigue strength of FRP composites under fatigue load. Meanwhile, the detrimental effect of increasing fiber Young’s modulus on fatigue strength becomes progressively more pronounced with rising temperature. In addition, the Young’s modulus level of the polymer matrix is relatively low, and increasing its modulus can significantly improve the mechanical properties, which is beneficial to the strength enhancement of FRP composites under fatigue load. Moreover, the enhancement of fatigue strength diminishes progressively with rising temperature. This can be interpreted that the lower matrix Young’s modulus at elevated temperatures has a negative effect on the fatigue strength. Consequently, this concludes that moderately increasing the matrix Young’s modulus enhances the fatigue resistance of FRP composites, whereas high level of the fiber Young’s modulus detrimentally impacts the fatigue performance. The influence of Young’s modulus of fibers and matrix on the fatigue strength of FRP composites at different temperatures.
Conclusions
In this study, we developed a TDDFS model for FRP composites based on the FHEED principle and damage failure analysis. The model enables accurate prediction of the fatigue strength of FRP composites under varying temperatures and damage levels. The combined effects of layup configurations, temperature, constituents’ properties evolution, and thermal induced damage variation on fatigue strength were systematically investigated and quantified in the model. The comparisons between model predictions and experimental data demonstrate that the TDDFS theoretical model reliably captures the temperature dependent fatigue behavior of FRP composites.
In addition, the analysis of factors affecting the fatigue properties of FRP composites was conducted, and providing the following insights: (1) As temperature increases, the damage degree escalates, and resulting in a decline in fatigue strength; (2) At a constant temperature, weakened bonding properties at the fiber/matrix interface promotes the damage accumulation, thereby reducing the fatigue strength; (3) While the increase of matrix Young’s modulus enhances the fatigue strength at lower temperatures, this beneficial effect diminishes with rising temperature. Conversely, excessively high Young’s modulus of the fiber induces the stress concentration phenomena, ultimately degrading the fatigue performance. This research develops a mechanistic theoretical framework for predicting temperature-dependent fatigue behavior in FRP composites, integrating damage progression physics with strength degradation mechanisms. The model provides both predictive capability and fundamental understanding of composite performance under thermo-mechanical fatigue conditions.
At the current stage, the validation and application scope of this model is explicitly confined to the fatigue strength prediction for unidirectional or orthogonally woven composites under constant ambient temperature. Accordingly, the model is not applicable to the following scenarios that have not been validated: (1) unconventional 3D woven or other complex braided composites; (2) structures with stress concentration effects such as holes, notches, or joints; (3) laminates subject to hygrothermal aging or moisture absorption effects; and (4) loading conditions with very high stress ratios.
Due to current experimental limitations, the test data of fatigue strength of FRP composites under high temperature conditions remain relatively scarce. Nevertheless, a systematic validation of the model was conducted within the temperature range covered by the available experimental data. The results demonstrate that the model predictions are in good agreement with the experimental measurements. It should be noted, however, that when the model is applied to predict the fatigue strength at elevated temperatures beyond the validation range, the predictions may tend to be higher than actual values. This overestimation is primarily attributed to the fact that the current model does not fully account for key factors such as the modulus reduction of CNTs/GNPs and the softening of interfacial strength at high temperatures. Therefore, the established model is explicitly applicable only to fatigue strength predictions within the temperature range from room temperature up to the glass transition temperature. The future work will focus on two main aspects: first, obtaining additional experimental data over a broader high-temperature range to further enhance the model’s applicability and generalizability. Afterwards, improving the physical foundation of the model under high-temperature conditions by incorporating mechanisms such as interfacial property degradation, thereby increasing its predictive accuracy and generalization ability.
Footnotes
Declaration of Conflicting Interests
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (Nos. 12202086, 12572238, 12302271); Natural Science Foundation of Chongqing, China (No. CSTB2025NSCQ-GPX0853); Science and Technology Research Program of Chongqing Municipal Education Commission (No. KJZD-K202500709); Opening Project of Key Laboratory of Testing Technology for Manufacturing Process, Southwest University of Science and Technology (No. 24kfzk03); Research and Innovation Program for Graduate Student in Chongqing (No. 2025S0097).
Data Availability Statement
The data that support the findings of this study are available from the corresponding author upon reasonable request.
