Abstract
This study predicted the fatigue life (N) of UHPC incorporated with different volume fractions ( = 0.0%, 0.5%, 1.0%, 1.5% and 2.0%) of steel fiber under flexural cyclic loading at various stress levels (S). The Weibull distribution, a two-parameter model, was utilized to estimate the distribution of fatigue life in UHPC. Subsequently, three methods were employed to calculate the parameters: the graphical method, the method of moments, and the method of maximum likelihood. The averaged values of these parameters were then obtained to enhance the accuracy of the estimation. The results are presented in the form of S-N diagrams, which depict the quantitative relationship between stress (S) and fatigue life (N). This relationship was determined using the Wohler equation, the modified Wohler equation, and the power equation. By employing these equations, the flexural fatigue strength of UHPC can be accurately predicted. Subsequently, the fatigue failure probability (Pf) was incorporated to enhance the reliability of the S-N quantitative relation. The fatigue testing results were presented in the form of S-N-Pf curves, which comprehensively reflect the relationship between stress, fatigue life, and failure probability. Furthermore, the mathematical relation of the S-N-Pf curves was derived to predict the fatigue life of UHPC with a given failure probability, providing a more comprehensive and accurate assessment of its fatigue behavior.
1 Introduction
UHPC has emerged as a novel cementitious composite material, renowned for its high strength and exceptional durability (; ; ). As a result, it has found widespread application in the construction of large bridges and high-rise buildings (). Notably, a significant number of these structures operate under cyclic loading conditions (; ; ). However, due to the inherent discreteness resulting from the multiphase and inhomogeneous nature of concrete materials at the meso-level, there is considerable variability in the fatigue test results for concrete. Therefore, it is crucial to accurately characterize the fatigue life distribution of UHPC to ensure its reliable performance in such structures. Ganesh () delved into the impact of stress levels (S) on the fatigue life (N) of UHPC and established that the scattered fatigue life (N) data can be statistically analyzed using the two-parameter Weibull distribution. This finding was further corroborated by Niu (), providing a robust framework for characterizing the fatigue behavior of UHPC. Therefore, due to its physically valid assumptions and robust experimental verification, the two-parameter Weibull distribution is chosen as an efficient method for quantitatively describing the distribution of fatigue life. Typically, the parameters of the Weibull distribution (the shape parameter α and the characteristic life u) are directly obtained using the graphical method. However, Singh () observed significant differences between the parameters derived from the graphical method and those obtained through the method of moments () and maximum likelihood estimation (). Consequently, it is imperative to utilize the graphical method, the method of moments, and the method of maximum likelihood to calculate the parameters α and u, enabling a more accurate assessment of their variations. This comprehensive approach will enhance our understanding of the fatigue life distribution of UHPC, providing a more reliable basis for structural design and evaluation.
The number of fatigue cycles that UHPC can endure under a specific stress level can be represented by S-N curves. employed the Wohler equation to plot the S-N diagram in order to predict the fatigue life of UHPC under constant amplitude tensile fatigue cycles. However, it is important to note that the applied minimum fatigue stress, , also has a significant impact on the fatigue life (N). Furthermore, in practical structures, the minimum value of the repeatedly applied stress is never zero. refined the Wohler equation by introducing a stress ratio R (R = Pmin/Pmax), which allowed for a more realistic simulation of loading conditions in actual structures. assessed the fatigue behavior of UHPC across seven different loading levels by utilizing the modified Wohler equation and derived the fatigue strength. However, a significant limitation of the modified Wohler equation is that it fails to meet the extreme boundary condition, where the stress level S approaches zero as the fatigue life N tends towards infinity. This limitation highlights the need for further improvements and considerations in accurately predicting the fatigue life of UHPC under various loading conditions (). Therefore, proposed a power formula to address the aforementioned shortcomings of the modified Wohler equation. applied the power equation to determine the flexural fatigue endurance limits of UHPC. The Wohler equation, modified Wohler equation, and power equation are commonly used methods for evaluating the relationship between S and N, and predicting the fatigue life of UHPC. These three equations mentioned above are employed to conduct a comprehensive quantitative analysis of the S-N relationship of UHPC under flexural fatigue loading.
Due to the heterogeneity of the UHPC composition, especially with the incorporation of steel fibers, the fatigue life exhibits significant fluctuations, which poses challenges in accurately predicting the fatigue behavior of UHPC. Therefore, it is imperative to incorporate the failure probability () into the fatigue testing data and present S-N-Pf curves to address the issue of discrete testing data. established the S-N-Pf curves using fatigue data and derived a mathematical equation to estimate the intricate relationship between stress (S), fatigue life (N), and failure probability (). Zhao () further explored the distribution of fatigue life under various stress levels in UHPC and predicted the fatigue life with different failure probabilities () by developing a numerical analysis model. conducted an analysis of the fatigue life of fiber reinforced concrete (FRC) with varying fiber content, considering different survival probabilities. Xu () took into account failure probabilities () and developed a single-logarithm fatigue equation specifically for FRC. The integration of probabilistic concepts for predicting flexural fatigue strength has become a widely adopted approach in the study of FRC. However, it is noteworthy that there are significant differences in both material composition and structure between FRC and UHPC. Therefore, the numerical prediction model developed for FRC cannot be directly applied to UHPC. As a result, it is crucial and meaningful to establish the S-N-Pf curves and develop a numerical prediction model specifically tailored for UHPC, based on its fatigue testing data. This approach will enable a more accurate understanding and prediction of the fatigue behavior of UHPC.
This study conducted an investigation into UHPC with different fiber volume fractions under varying applied stress levels. Subsequently, the parameters of fatigue equations representing the S-N curves were estimated, rendering these equations suitable for predicting the flexural fatigue strength of UHPC. The failure probability () was then incorporated into the fatigue testing data, leading to the development of a family of S-N- curves. Finally, a numerical prediction model tailored specifically for UHPC was established to quantitatively evaluate its flexural fatigue strength.
2 Experimental program
2.1 Materials and mixtures composition
The cementitious materials employed in this study were PII 52.5 portland cement and silica fume, the size distribution is shown in Figure 1 and the chemical composition is presented in Table 1. The silica fume contains 98% SiO2 with the averaged diameter of 0.2 μm. By using particle packing theory, the fine aggregate with two ranges (0.160–0.315 mm and 0.63–1.25 mm) was adopted to increase the stacking density. Superplasticizer, solid content 30%, was introduced to improve the workability for a low water-binder ratio fresh matrix. The properties of the steel fiber were presented in Table 2. The UHPC mixtures adopted to cast the tested specimens are listed in Figure 2. The water binder ratio (W/C) was 0.16 and the binder-sand ratio was 1:1, as shown in Table 3.
FIGURE 1
TABLE 1
| Cementitious materials | SiO2 | Al2O3 | Fe2O3 | CaO | MgO | K2O | Na2O | SO3 | Others | LOI |
|---|---|---|---|---|---|---|---|---|---|---|
| Cement | 21.60 | 4.35 | 2.95 | 63.81 | 1.76 | 0.51 | 0.16 | 2.06 | 1.61 | 1.19 |
| Silica fume | 98.07 | — | 0.12 | 0.51 | 0.31 | 0.53 | 0.14 | 0.12 | 0.19 | 0.01 |
Chemical compositions of the cementitious materials.
Note: LOI, loss on ignition.
TABLE 2
| Diameter, df (mm) | 0.2 |
| Length, lf (mm) | 6.0/13.0/20.0 |
| Aspect ratio, lf/df | 30.0/65.0/100.0 |
| Density, ρ (g/cm3) | 7.8 |
| Tensile strength, fst (MPa) | 2,500 |
| Elastic modulus, Ef (GPa) | 200 |
Properties of the steel fiber.
FIGURE 2
TABLE 3
| Cement | Silica fume | Sand | Volume fraction of steel fiber(%) | Superplasticizer | Water | |
|---|---|---|---|---|---|---|
| 0.16–0.315 mm | 0.63–1.25 mm | |||||
| 770 | 230 | 300 | 700 | 0.0 | 35 | 135.5 |
| 0.5 | ||||||
| 1.0 | ||||||
| 1.5 | ||||||
| 2.0 | ||||||
Proportions of components in the UHPC mixture ().
The mixing process was conducted in a double horizontal shaft mortar mixer with a volume fraction of 100 L. At the beginning of the mixing process, cement and silica fume were premixed for 2 min. Two types of silica sands with different particle diameters were added and stirred for another 5 min. Water and superplasticizer were then poured into the dry mixture and stirred for 10 min to provide the mortar matrix with considerable fluidity and viscosity. Finally, a small amount of steel fibers were added gradually into the fresh matrix and stirred for 5 min to improve the dispersion of steel fibers. The mixture was then poured along a single side of a mold to achieve a similar orientation distribution of the steel fibers, which can greatly improve the bending property of the specimen. The mold was vibrated for 90 s to enhance the compactness of the mixture. After the fresh concrete poured into the mould, the plastic sheets were covered on the top surface of the prepared specimens to reduce the water evaporation, then the hardened specimens were cured at room temperature for 24 h before demolding; Subsequently, the tested specimens were maintained into the water tank for 90 days to minimize the effect of the strength increase of the UHPC during the fatigue testing.
2.2 Fatigue testing
The dimensions of the tested specimens were 100 mm × 100 mm × 400 mm. Three beams were selected to measure the flexural strength and their values were averaged. A total of ten beams were utilized to assess the flexural fatigue strength. Based on the fatigue testing data, S-N diagrams were plotted for UHPC under various stress levels. To establish the appropriate stress levels, initial static flexural testing was conducted using a 100-kN MTS testing machine. The span for the flexural testing was set at 300 mm, and the loading frame operated in displacement-control mode with a constant rate of 0.02 mm/min, as shown in Figure 2. The static flexural stress of UHPC specimens containing fiber volume fractions of 0.5%, 1.0%, 1.5%, and 2.0% were 42.31, 54.47, 58.42, and 69.05 kN, respectively. For the fatigue testing, a loading frequency of 8 Hz was used with a constant-amplitude sinusoidal waveform. At present, there is no unified regulation for fatigue loading frequency. According to the performance of the testing machine, the cost and period of the test, the loading frequency of 5–15 Hz is usually adopted. When the loading frequency is between 100–900 times per minute, it has no obvious effect on the fatigue strength of concrete (). If the loading speed is too slow, the creep will increase, and the fatigue strength or life will be reduced. In this paper, the loading frequency of 480 times per minute is selected, which is within a reasonable range.
Initially, the load was linearly increased to the average value, , and then the specimens were subjected to cyclic loading with a wave range from to . The stress levels employed in the fatigue testing were 0.80, 0.75, 0.70, and 0.65, with a constant stress ratio () of 0.1. For each specimen under these given loading conditions, the number of cycles to failure was measured by the cycle counter of the testing machine. A maximum number of allowed cycles was specified, with a default value of two million cycles (). If fatigue failure did not occur within this limit, the testing was terminated.
3 Results and discussion
3.1 Fatigue test results
To establish the appropriate cyclic load stress level, it is crucial to first evaluate the flexural strength of the UHPC. The flexural strength of UHPC exhibits a positive correlation with fiber content. Specifically, when the fiber volume content () is 0.0%, the bending strength of the UHPC matrix stands at 8.8 MPa. However, as the fiber volume content increases to 0.5%, 1.0%, 1.5%, and 2.0%, the bending strength of UHPC rises significantly by 44.20%, 85.68%, 99.21%, and 135.00%, respectively.
The fatigue lives of UHPC specimens containing different volumes of steel fibers ( = 0.0%, 0.5%, 1.0%, 1.5%, and 2.0%) under various stress levels are summarized in Tables 4. To verify the fatigue life distribution of the UHPC using Weibull distribution function, the fatigue life (N) of the UHPC under a certain stress level (S) allowed to be sorted in ascending order. When comparing specimens under the same stress level, it was observed that as the fiber content increased, the fatigue life of the UHPC specimens gradually improved. Conversely, for specimens containing the same quantity of steel fibers, a decrease in the stress level led to a gradual increase in their fatigue life.
TABLE 4
| Fiber contents | Stress levels | S = 0.85 | S = 0.80 | S = 0.75 | S = 0.70 | S = 0.65 |
|---|---|---|---|---|---|---|
| 0.0% | Fatigue life | 54 | 445 | 6,467 | 21,718 | |
| 18 | 87 | 589 | 7,726 | 31,118 | ||
| 19 | 91 | 645 | 8,936 | 44,456 | ||
| 21 | 92 | 789 | 9,262 | 55,499 | ||
| 24 | 117 | 803 | 11,569 | 61,623 | ||
| 25 | 118 | 946 | 12,798 | 64,227 | ||
| 36 | 129 | 948 | 14,732 | 65,620 | ||
| 37 | 135 | 1,078 | 19,214 | 73,200 | ||
| 40 | 151 | 1,278 | 21,433 | 87,943 | ||
| 43 | 216 | 1743 | 22,903 | 93,666 | ||
| 0.5% | Fatigue life | 787 | 2,345 | 34,345 | ||
| 100 | 788 | 3,634 | 10,526 | 58,789 | ||
| 105 | 1,199 | 5,467 | 14,068 | 88,754 | ||
| 111 | 1,328 | 7,689 | 15,114 | 107,865 | ||
| 161 | 1,522 | 8,823 | 21,393 | 124,560 | ||
| 193 | 2,202 | 9,631 | 26,533 | 156,074 | ||
| 200 | 2,338 | 10,231 | 29,215 | 185,461 | ||
| 242 | 2,840 | 11,416 | 32,806 | 227,923 | ||
| 325 | 3,234 | 12,526 | 40,806 | 236,811 | ||
| 344 | 3,634 | 14,327 | 74,015 | 282,145 | ||
| 1.0% | Fatigue life | 1751 | 14,843 | 25,578 | 35,359 | |
| 143 | 1869 | 15,160 | 30,426 | 76,218 | ||
| 226 | 2,917 | 18,645 | 44,672 | 108,236 | ||
| 247 | 3,404 | 25,184 | 55,622 | 144,496 | ||
| 270 | 4,406 | 26,400 | 76,218 | 189,782 | ||
| 324 | 5,679 | 33,527 | 91,246 | 207,436 | ||
| 426 | 6,941 | 42,345 | 105,376 | 242,496 | ||
| 537 | 7,654 | 51,640 | 134,984 | 289,212 | ||
| 681 | 8,912 | 61,616 | 179,240 | 365,789 | ||
| 834 | 10,245 | 76,117 | 207,856 | 420,385 | ||
| 1.5 | Fatigue life | 616 | 4,445 | 19,916 | 25,476 | 105,530 |
| 772 | 6,549 | 25,528 | 45,649 | 184,434 | ||
| 837 | 7,439 | 32,983 | 59,727 | 205,010 | ||
| 1,131 | 8,862 | 46,033 | 63,186 | 388,448 | ||
| 1,385 | 9,038 | 53,317 | 105,000 | 487,886 | ||
| 1,578 | 15,411 | 67,054 | 135,000 | 590,789 | ||
| 1735 | 19,800 | 74,271 | 183,456 | 620,115 | ||
| 2,211 | 25,066 | 98,885 | 204,900 | 752,633 | ||
| 2,585 | 26,601 | 101,377 | 227,900 | 811,793 | ||
| 3,163 | 28,865 | 181,823 | 290,789 | 1,020,000 | ||
| 2.0% | Fatigue life | 10,269 | 26,983 | 65,111 | — | |
| 2072 | 12,533 | 44,651 | 76,334 | — | ||
| 2,245 | 13,560 | 69,074 | 89,074 | — | ||
| 2,456 | 16,795 | 82,586 | 122,586 | — | ||
| 3,786 | 20,140 | 123,402 | 189,320 | — | ||
| 4,356 | 24,632 | 142,398 | 226,616 | — | ||
| 5,924 | 42,126 | 177,635 | 306,375 | — | ||
| 7,730 | 44,651 | 187,447 | 475,452 | — | ||
| 7,854 | 63,181 | 212,538 | 588,342 | — | ||
| 9,012 | 65,708 | 299,859 | — |
Fatigue life data for UHPC with various of fiber contents.
3.2 Distribution of fatigue life
Due to the inherent uncertainties arising from internal defects in concrete structures and the inhomogeneous distribution of fibers, the fatigue life of UHPC exhibits significant variability. Consequently, numerous mathematical probability models have been proposed to statistically analyze the distribution of UHPC’s fatigue life. Presently, the Weibull distribution stands as a commonly employed mathematical probability statistical model for analyzing the fatigue life distribution of concrete, offering a robust framework for quantifying and understanding this complex phenomenon.
The Weibull probability distribution function f(n) and cumulative distribution function are expressed as Equations 1, 2 respectively ():n: specific value of random variable N; α and u: shape parameter of Weibull distribution and characteristic value of fatigue life, respectively; : minimum value of position parameter or fatigue life.
According to Equation 2, the survival function of fatigue life, that is, the reliability function , can be obtained as Equation 3: ()
For steel-fiber-reinforced concrete, assumed that the minimum fatigue life of the Weibull distribution function in practice is = 0, the equation can be simplified as:
After logarithm fetch on Equation 4 twice, the linear relation can be demonstrated as ():
For the Equation 5, the relationship between and is linear, which can be applied to assess the distribution of the fatigue life of UHPC whether is in accord with the Weibull function.
To verify the fatigue life distribution of the UHPC using Equation 5, the fatigue life (N) of the UHPC under a certain stress level (S) allowed to be sorted in ascending order. And the reliability function can be obtained by using Equation 6:i: sequence number of a fatigue data value in the sequence of fatigue life; k: number of the tested specimens under cyclic loading with an given stress level; if there existed a linear relationship between and , it is reasonably concluded that the fatigue life distribution of UHPC behave according to the Weibull function. For example, the relationships between the UHPC fatigue life and reliability function values for UHPC with 1.5% steel fiber are presented in Tables 5.
TABLE 5
| Stress level | i | |||||
|---|---|---|---|---|---|---|
| 0.85 | 1 | 616 | 0.9091 | −2.3506 | 6.4232 | 0.0909 |
| 2 | 772 | 0.8182 | −1.6061 | 6.6489 | 0.1818 | |
| 3 | 837 | 0.7273 | −1.1443 | 6.7298 | 0.2727 | |
| 4 | 1,131 | 0.6364 | −0.7941 | 7.0308 | 0.3636 | |
| 5 | 1,385 | 0.5455 | −0.5007 | 7.2335 | 0.4545 | |
| 6 | 1,578 | 0.4545 | −0.2377 | 7.3639 | 0.5455 | |
| 7 | 1735 | 0.3636 | 0.0115 | 7.4588 | 0.6364 | |
| 8 | 2,211 | 0.2727 | 0.2618 | 7.7012 | 0.7273 | |
| 9 | 2,585 | 0.1818 | 0.5334 | 7.8575 | 0.8182 | |
| 10 | 3,163 | 0.0909 | 0.8746 | 8.0593 | 0.9091 | |
| 0.80 | 1 | 4,445 | 0.9091 | −2.3506 | 8.3995 | 0.0909 |
| 2 | 6,549 | 0.8182 | −1.6061 | 8.7871 | 0.1818 | |
| 3 | 7,439 | 0.7273 | −1.1443 | 8.9145 | 0.2727 | |
| 4 | 8,862 | 0.6364 | −0.7941 | 9.0895 | 0.3636 | |
| 5 | 9,038 | 0.5455 | −0.5007 | 9.1092 | 0.4545 | |
| 6 | 15,411 | 0.4545 | −0.2377 | 9.6428 | 0.5455 | |
| 7 | 19,800 | 0.3636 | 0.0115 | 9.8934 | 0.6364 | |
| 8 | 25,066 | 0.2727 | 0.2618 | 10.1293 | 0.7273 | |
| 9 | 26,601 | 0.1818 | 0.5334 | 10.1887 | 0.8182 | |
| 10 | 28,865 | 0.0909 | 0.8746 | 10.2704 | 0.9091 | |
| 0.75 | 1 | 19,916 | 0.9091 | −2.3506 | 9.8993 | 0.0909 |
| 2 | 25,528 | 0.8182 | −1.6061 | 10.1475 | 0.1818 | |
| 3 | 32,983 | 0.7273 | −1.1443 | 10.4038 | 0.2727 | |
| 4 | 46,033 | 0.6364 | −0.7941 | 10.7371 | 0.3636 | |
| 5 | 53,317 | 0.5455 | −0.5007 | 10.8840 | 0.4545 | |
| 6 | 67,054 | 0.4545 | −0.2377 | 11.1133 | 0.5455 | |
| 7 | 74,271 | 0.3636 | 0.0115 | 11.2155 | 0.6364 | |
| 8 | 98,885 | 0.2727 | 0.2618 | 11.5017 | 0.7273 | |
| 9 | 101,377 | 0.1818 | 0.5334 | 11.5266 | 0.8182 | |
| 10 | 181,823 | 0.0909 | 0.8746 | 12.1108 | 0.9091 | |
| 0.70 | 1 | 25,476 | 0.9091 | −2.3506 | 10.1455 | 0.0909 |
| 2 | 45,649 | 0.8182 | −1.6061 | 10.7287 | 0.1818 | |
| 3 | 59,727 | 0.7273 | −1.1443 | 10.9975 | 0.2727 | |
| 4 | 63,186 | 0.6364 | −0.7941 | 11.0538 | 0.3636 | |
| 5 | 105,000 | 0.5455 | −0.5007 | 11.5617 | 0.4545 | |
| 6 | 135,000 | 0.4545 | −0.2377 | 11.8130 | 0.5455 | |
| 7 | 183,456 | 0.3636 | 0.0115 | 12.1197 | 0.6364 | |
| 8 | 204,900 | 0.2727 | 0.2618 | 12.2303 | 0.7273 | |
| 9 | 227,900 | 0.1818 | 0.5334 | 12.3367 | 0.8182 | |
| 10 | 290,789 | 0.0909 | 0.8746 | 12.5804 | 0.9091 | |
| 0.65 | 1 | 105,530 | 0.9091 | −2.3506 | 11.5668 | 0.0909 |
| 2 | 184,434 | 0.8182 | −1.6061 | 12.1251 | 0.1818 | |
| 3 | 205,010 | 0.7273 | −1.1443 | 12.2308 | 0.2727 | |
| 4 | 388,448 | 0.6364 | −0.7941 | 12.7021 | 0.3636 | |
| 5 | 487,886 | 0.5455 | −0.5007 | 13.0978 | 0.4545 | |
| 6 | 590,789 | 0.4545 | −0.2377 | 13.2892 | 0.5455 | |
| 7 | 620,115 | 0.3636 | 0.0115 | 13.3377 | 0.6364 | |
| 8 | 752,633 | 0.2727 | 0.2618 | 13.5313 | 0.7273 | |
| 9 | 811,793 | 0.1818 | 0.5334 | 13.6070 | 0.8182 | |
| 10 | 1,020,000 | 0.0909 | 0.8746 | 13.8353 | 0.9091 |
Relationships between and of UHPC with different stress levels.
Using the same method, the relationship between the fatigue life of the UHPC and the reliability function can be obtained under different stress levels with different steel fiber volume content ( = 0.0%, 0.5%, 1.0% and 2.0%). Figure 3 presents a linear regression analysis of the fatigue life of UHPC under different stress levels using the image method. As evident from the figure, the fatigue life exhibits a predominantly linear distribution, with a minimum correlation coefficient of 90.74%. This indicates that the fatigue life of UHPC aligns with the Weibull parameter distribution, validating the use of the Weibull model for analyzing its fatigue behavior.
FIGURE 3
3.3 Determination of the weibull distribution parameters
3.3.1 Graphical method
The graphical method not only determines whether the fatigue life of UHPC adheres to the Weibull distribution function but also enables the extraction of critical parameters such as the shape parameter α and the characteristic fatigue life u. Through linear regression analysis, the slope and intercept of the line correspond to α and u, respectively, as illustrated in Figure 3.
3.3.2 Method of moment
To evaluate the distribution parameters of the Weibull functions using the method of moments, it is required an appropriate sample moments (mean and variance of the samples), which can be calculated by the following Equations 7, 8 ():
Sample moments:
E(n): average expected value; T(x): gamma function.
Sample variance:
μ = E(n): the average value of the fatigue tested data under an given stress level.
According to Equations 7, 8, the coefficient of variation (CV) of the fatigue life of the UHPC can be calculated as follows:
, σ: standard deviation.
Equations 10, 11 can be used to evaluate the Weibull function distribution parameters (α and u), respectively.
3.3.3 Maximum likelihood equation
The maximum likelihood equation is expressed as follows (): and : maximum likelihood estimates of α and θ, respectively. The shape parameter α can be calculated iteratively using Equation 13. In the iterative process, the shape parameter α, which can be obtained using the image method or moment method, is used as the first estimated value, and the maximum likelihood estimate is then iterated step by step. By substituting the maximum likelihood estimate of the shape parameter into Equation 12, can be calculated. Finally, the value of the characteristic fatigue life u is calculated according to Equation 14. The Weibull distribution function parameters obtained by different mathematical calculation methods are shown in Tables 6.
TABLE 6
| Stress levels | 0.85 | 0.80 | 0.75 | 0.70 | 0.65 | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| Vf = 0.0% | α | u | α | u | α | u | α | u | α | u |
| Graphical method | 2.77 | 29 | 2.59 | 135 | 2.48 | 1,053 | 2.21 | 15,490 | 2.14 | 69,494 |
| Method of moments | 3.28 | 33 | 2.91 | 133 | 2.64 | 1,042 | 2.45 | 15,227 | 2.33 | 67,248 |
| Method of maximum likelihood | 3.58 | 33 | 2.98 | 133 | 2.80 | 925 | 2.64 | 15,271 | 2.41 | 67,059 |
| Average | 3.21 | 32 | 2.83 | 134 | 2.64 | 1,007 | 2.43 | 15,329 | 2.29 | 67,934 |
| Vf = 0.5% | α | u | α | u | α | u | α | u | α | u |
|---|---|---|---|---|---|---|---|---|---|---|
| Graphical method | 2.02 | 229 | 1.75 | 2,311 | 1.69 | 10,172 | 1.58 | 33,897 | 1.49 | 176,472 |
| Method of moments | 2.10 | 222 | 1.96 | 2,243 | 1.77 | 9,713 | 1.56 | 32,714 | 1.42 | 169,422 |
| Method of maximum likelihood | 2.38 | 223 | 2.02 | 2,254 | 1.88 | 9,698 | 1.75 | 33,282 | 1.54 | 169,810 |
| Average | 2.13 | 225 | 1.91 | 2,269 | 1.78 | 9,861 | 1.63 | 33,298 | 1.48 | 171,901 |
| Vf = 1.0% | α | u | α | u | α | u | α | u | α | u |
|---|---|---|---|---|---|---|---|---|---|---|
| Graphical method | 1.89 | 475 | 1.55 | 6,273 | 1.45 | 42,311 | 1.38 | 109,908 | 1.30 | 245,178 |
| Method of moments | 1.81 | 461 | 1.88 | 6,059 | 1.68 | 41,124 | 1.65 | 105,957 | 1.47 | 233,414 |
| Method of maximum likelihood | 2.09 | 465 | 1.92 | 6,092 | 1.87 | 41,491 | 1.67 | 101,091 | 1.49 | 233,846 |
| Average | 1.97 | 467 | 1.78 | 6,141 | 1.63 | 41,624 | 1.55 | 105,652 | 1.42 | 237,479 |
| Vf = 1.5% | α | u | α | u | α | u | α | u | α | u |
|---|---|---|---|---|---|---|---|---|---|---|
| Graphical method | 1.81 | 1,850 | 1.45 | 17,710 | 1.44 | 80,528 | 1.29 | 179,328 | 1.32 | 309,474 |
| Method of moments | 1.82 | 1,853 | 1.71 | 17,056 | 1.58 | 77,630 | 1.55 | 149,117 | 1.48 | 296,008 |
| Method of maximum likelihood | 1.93 | 1,876 | 1.74 | 17,209 | 1.64 | 78,946 | 1.61 | 150,083 | 1.52 | 297,645 |
| Average | 1.82 | 1,859 | 1.64 | 17,325 | 1.55 | 79,035 | 1.48 | 159,509 | 1.39 | 301,045 |
| Vf = 2.0% | α | u | α | u | α | u | α | u | α | u |
|---|---|---|---|---|---|---|---|---|---|---|
| Graphical method | 1.62 | 5,831 | 1.39 | 36,308 | 1.32 | 159,869 | 1.16 | 278,516 | — | — |
| Method of moments | 1.60 | 5,591 | 1.55 | 30,638 | 1.51 | 152,983 | 1.48 | 219,230 | — | — |
| Method of maximum likelihood | 1.83 | 5,619 | 1.67 | 35,366 | 1.61 | 153,641 | 1.58 | 269,403 | — | — |
| Average | 1.70 | 5,680 | 1.54 | 34,104 | 1.48 | 155,498 | 1.32 | 255,716 | — | — |
Weibull distribution parameters of UHPC.
With the same fiber volume content, the shape parameter α exhibited a decreasing trend as the stress level increased. Specifically, under higher stress levels, the shape parameter α decreased significantly, indicating a narrower scatter range in the fatigue life data. Conversely, at lower stress levels, the decrease in α was more gradual, suggesting a wider variation in fatigue life. This contrasting behavior underscores the complexity of fatigue life distribution in UHPC, particularly across different stress levels. Figure 4 presents the mean change of the shape parameter α in the Weibull function of UHPC with various steel fiber content. Notably, the average value of the shape parameter α for UHRC ( = 0.0%) is higher than that for UHPC containing steel fibers under different stress levels. This observation suggests that incorporating steel fibers enhances the dispersion of fatigue life data. Furthermore, under the same stress level, the shape parameter α decreases as the fiber volume content increases, indicating that the degree of discretization increases with a higher steel fiber volume fraction.
FIGURE 4
3.3.4 Distribution of fibers on fracture surface
It is evident that the distribution characteristics of fibers on the fracture surface of UHPC are closely related to fatigue life. To evaluate the fiber distribution, an image technique was used, as shown in Figure 5. The fiber number was quantified by the fiber dispersion coefficient, as shown in Equation 15 () and the distribution characteristics of fiber orientation and pullout length were evaluated by Equations 16, 17, respectively.
FIGURE 5
Where α is the fiber dispersion coefficient, x0 is the average number of steel fibers in each unit, xi is the measured number per unit, and n is the total number of units.
Where is the packing density; θ, and l are the fiber inclined angle, diameter, and fiber length, respectively.
Where p(θ) is the probability distribution function, p and q are shape parameters with values greater than 1/2, and the values of θ range from 0 to π/2.
The fiber dispersion coefficients (α) of the UHPC contained 0.5%, 1.0%, 1.5% and 2.0% steel fiber, were 0.86, 0.83, 0.80 and 0.74, respectively. The fiber dispersion coefficient was close to 1, indicating that the fibers were well distributed in UHPC. With the increase of fiber contents, the dispersion coefficient gradually decreased, which means the uniformity of fiber distribution on the fracture surface decreased, thus leading to an increase in the discreteness of UHPC fatigue life.
The probability density distribution of fiber orientation for UHPC is represented in Figure 6. The fiber orientation was more concentrated at 20°–60°. This, in turn, affected the fatigue life of the UHPC. A comparison with the results of Xia et al. () showed that their proposed two-parameter exponential could accurately predict the probability density distribution of fiber orientation. The probability distribution function of the pull-out lengths of UHPC was expressed in Figure 7. After UHPC cracking, fibers across the cracks provided a bridging action, improving the fatigue life of UHPC, especially the stable crack propagation stages, which accounts for about 75% of the entire fatigue life.
FIGURE 6
FIGURE 7
3.3.5 Test of goodness of fit
The shape parameter α and characteristic fatigue life u in the Weibull distribution function were quantitatively assessed using the image method, moment method, and maximum likelihood estimation method. To enhance the credibility of the α and u values obtained through these methods and to verify the efficacy of the Weibull probability distribution function in analyzing the fatigue life distribution of UHPC, the Kolmogorov-Smirnov (K-S) test was employed to conduct a goodness-of-fit analysis on the fatigue life data as shown in Equation 18. This approach, as outlined in (), allowed for a rigorous evaluation of how well the Weibull distribution matched the actual fatigue life data, providing further validation for its application in modeling the fatigue behavior of UHPC.
: cumulative histogram; : hypothesized cumulative distribution function.
Table 7 presents the results obtained by the K-S method for the case of S = 0.80 for the UHPC incorporated with the steel fiber volume = 1.5%.
TABLE 7
| i | ||||
|---|---|---|---|---|
| 1 | 4,445 | 0.1 | 0.1015 | 0.0015 |
| 2 | 6,549 | 0.2 | 0.1831 | 0.0169 |
| 3 | 7,439 | 0.3 | 0.2207 | 0.0793 |
| 4 | 8,862 | 0.4 | 0.2828 | 0.1172 |
| 5 | 9,038 | 0.5 | 0.2906 | 0.2094 |
| 6 | 15,411 | 0.6 | 0.5617 | 0.0382 |
| 7 | 19,800 | 0.7 | 0.7121 | 0.0121 |
| 8 | 25,066 | 0.8 | 0.8403 | 0.0403 |
| 9 | 26,601 | 0.9 | 0.8677 | 0.0323 |
| 10 | 28,865 | 1.0 | 0.9011 | 0.0989 |
K-S test for fatigue life.
When the stress level is S = 0.80, the maximum difference in the Kolmogorov-Smirnov (K-S) test table for the fatigue life of UHPC with a steel fiber content of = 1.5% is 0.2094. In this study, the sample size was fixed at 10, and the significance level was set at 5%. Under these conditions, the critical value is 0.4092. Since the maximum difference is less than the critical value (), it indicates that the Weibull probability distribution function is well-suited for assessing the fatigue life distribution of UHPC ().
3.4 S-N relationship for UHPC
Due to the discrete fatigue life of concrete, different mathematical models based the fatigue testing data was developed the fatigue life. Kaushik () and applied the Wholer equation, Equation 19, to build up the mathematical dependence of the stress level S and fatigue life N. Flexural strength of concrete; : Maximum fatigue stress; A and B are experimental coefficients.
Tepfers () introduced the concept of stress ratio R, which was based on Wohler’s equation. This stress ratio allowed for a redefined relationship between the stress level S and fatigue life N, as follows:β: material coefficient. , which is introduced to describe the real situation of cyclic loads () on engineering structures.
proposed using a power function to clarify the relationship between the S and N, as shown in Equation 21:, : test coefficients. The advantage of this method is that it includes an extreme case wherein the concrete fatigue life is infinite when S is close to zero. With a purpose to simplify calculation, taking Log on both sides of Equation 22:
By conducting a linear regression analysis on the S-N curves, the coefficients A and B in Equation 16 can be determined, as illustrated in Figure 8. Utilizing the fatigue testing data, a quantitative evaluation of the relationship between S, N, and R is performed, enabling the determination of the value of β. Subsequently, Equation 20 can be employed to predict the fatigue life of UHPC. Similarly, C1 and C2 can be calculated through linear regression analysis, as demonstrated in Figure 8.
FIGURE 8
3.5 Failure probability using S-N relationships for UHPC
The actual test results reveal that the fatigue life data of UHPC still demonstrated large dispersion when the fatigue test procedure is strictly controlled. Therefore, the S-N relationship becomes more convincing with the introducing the failure probability .
First, the fatigue life test data of UHPC with different fiber content under a certain stress level sorted in descending order, as presented in Tables 4. Then, the reliability function can be obtained by using Equation 6, the failure probability () can also be described as 1-; the quantitative relationship between and N was obtained. Next, the fatigue life of the UHPC under different stress levels (S-N) corresponding to the failure probability () was calculated according to Equations 23, 24. Finally, based on the S-N regression line, Equation 6 was used to calculate the value of failure probability () value matched with different stress levels S under the condition of the same fatigue life (S-).
The family of S–N– curves developed for UHPC with different fiber volume fractions ( = 0.0%, 0.5%, 1.0%, 1.5% and 2.0%) is presented in Figures 9–13, respectively.
FIGURE 9
FIGURE 10
FIGURE 11
FIGURE 12
FIGURE 13
Figures 9–13Figures 14 present the family of S-N- diagrams, offering a qualitative analysis of the relationship between stress (S), fatigue life (N), and failure probability (Pf). However, it is evident that these diagrams alone are insufficient for accurately predicting the fatigue life of UHPC. Therefore, utilizing mathematical analysis to estimate the fatigue life of UHPC under various stress levels, considering the failure probability (), is crucial for practical engineering applications. This approach is formally expressed in Equation 25 (; ).a, b and c: experimental coefficients.
FIGURE 14
The S-N-
relationship of the UHPC with different fiber volume contents can be expressed as shown in
Equations 26–
30:
UHPC ( = 0.0%):
UHPC ( = 0.5%):
UHPC ( = 1.0%):
UHPC ( = 1.5%):
UHPC ( = 2.0%):
Based on the experimental coefficients obtained in this study, the predicted curves are presented alongside the fatigue test curves in Figures 9–13. A slight discrepancy exists between the predicted and experimental curves, indicating that Equation 25 is appropriately used for evaluating the fatigue life of UHPC under various stress levels with a specified survival or failure probability.
4 Conclusion
This study evaluated the fatigue life (N) of UHPC reinforced with different volume fractions (
= 0.0%, 0.5%, 1.0%, 1.5% and 2.0%) of steel fiber under flexural cyclic loading at various stress levels (S). The Weibull distribution was utilized to assess the distribution of fatigue life for UHPC. A quantitative relationship between S and N was established to predict the flexural fatigue strength of UHPC. Furthermore, a mathematical model describing the S-N-P
fcurves was developed to estimate the fatigue life of UHPC for a given failure probability. Based on these findings, the following conclusions were draw:
(1) The fatigue life is essentially linearly distributed, and the minimum correlation coefficient is 90.74%, that is, the fatigue life of UHPC conforms to the Weibull parameter distribution.
(2) The discretization of the fatigue life of UHPC increased with the increase of fiber content.
(3) The quantitative relation between S and N were also describe quantitatively and the fatigue life of UHPC under different stress levels can be predicated.
(4) The mathematical models aimed to predicated the fatigue life of the UHPC with different fiber content at a given failure probability were established.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.
Author contributions
CL: Writing–original draft, Writing–review and editing. PY: Data curation, Methodology, Writing–review and editing. YN: Writing–original draft, Writing–review and editing. YZ: Project administration, Validation, Writing–review and editing. CC: Writing–review and editing, Project administration, funding acquisition
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The financial assistance received from National Youth Natural Science Foundation of China (Project No. 52108200), GuangDong Basic and Applied Basic Research Foundation (2023A1515011444).
Conflict of interest
Authors CL and PY were employed by Guangzhou Guangjian Construction Engineering Testing Center Co. Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AmranM.HuangS.-shanOnaiziA. M.MakulN.AbdelgaderH. S.OzbakkalogluT. (2022). Recent trends in ultra-high performance concrete (UHPC): current status, challenges, and future prospects. Constr. Build. Mater.352, 129029. 10.1016/j.conbuildmat.2022.129029
2
Baek-SikT. R. (1979). Fatigue strength of plain, ordinary, and light weight concrete. Am. Soc. Civ. Eng.8 (42), 159–167.
3
EFNARC (2005). The European guidelines for self-compacting concrete, specification, production and use, 68.
4
GaneshP.MurthyR. (2022). Flexural fatigue strains of constituent materials in strengthened RC beams with UHPC strips. Int. J. Fatigue21, 107351. 10.1016/j.ijfatigue.2022.107351
5
GaoP.YeG.HuangH.QianZ.SchlangenE.WeiJ.et al (2022). Incorporating elastic and creep deformations in modelling the three-dimensional autogenous shrinkage of cement paste. Cem. Concr. Res.160, 106907. 10.1016/j.cemconres.2022.106907
6
GoelS.SinghS. P. (2014). Fatigue performance of plain and steel fibre reinforced self compacting concrete using S–N relationship. Eng. Struct.74, 65–73. 10.1016/j.engstruct.2014.05.010
7
GoelS.SinghS. P.SinghP. (2012). Flexural fatigue strength and failure probability of self compacting fibre reinforced concrete beams. Eng. Struct.40, 131–140. 10.1016/j.engstruct.2012.02.035
8
HacèneH.GhofraneB.GérardD. (2014). Flexural fatigue performance of metal steel fibre reinforced mortar – influence of fibre aspect ratio and type. Constr. Build. Mater.58, 166–170. 10.1016/j.conbuildmat.2014.02.016
9
HuangJ.QiuS.RodrigueD. (2022). Parameters estimation and fatigue life prediction of sisal fibre reinforced foam concrete. J. Mater. Res. Technol.20, 381–396. 10.1016/j.jmrt.2022.07.096
10
JunZ.LiG.WangZ.ZhaoX. L. (2019). Fatigue behavior of concrete beams reinforced with glass- and carbon-fiber reinforced polymer (GFRP/CFRP) bars after exposure to elevated temperatures. Compos. Struct.229, 111427. 10.1016/j.compstruct.2019.111427
11
LiL.GuanJ.PengY.YinY.LiY. (2021). A Weibull distribution-based method for the analysis of concrete fracture. Eng. Fract. Mech.256, 107964. 10.1016/j.engfracmech.2021.107964
12
LiL.XuL.HuangLeXuF.HuangY.CuiK.et al (2022). Compressive fatigue behaviors of ultra-high performance concrete containing coarse aggregate. Cem. Concr. Compos.128, 104425. 10.1016/j.cemconcomp.2022.104425
13
MakitaT.BrühwilerE. (2014). Tensile fatigue behaviour of ultra-high performance fibre reinforced concrete (UHPFRC). Mater. Struct.47 (3), 475–491. 10.1617/s11527-013-0073-x
14
NiuY.HuangH.WeiJ.JiaoC.MiaoQ. (2022b). Investigation of fatigue crack propagation behavior in steel fiber-reinforced ultra-high-performance concrete (UHPC) under cyclic flexural loading. Compos. Struct.282, 115126. 10.1016/j.compstruct.2021.115126
15
NiuY.WeiJ.JiaoC. (2022a). Multi-scale fiber bridging constitutive law based on meso-mechanics of ultra-high performance concrete under cyclic loading. Constr. Build. Mater.354, 129065. 10.1016/j.conbuildmat.2022.129065
16
OhB. H. (1986). Fatigue analysis of plain concrete in flexure. Struct. Eng.112 (2), 273–288. 10.1061/(asce)0733-9445(1986)112:2(273)
17
SahooS.LakavathC.PrakashS. S. (2021). Experimental and analytical studies on fracture behavior of fiber-reinforced structural lightweight aggregate concrete. J. Mater. Civ. Eng.33 (5), 126281. 10.1061/(asce)mt.1943-5533.0003680
18
SahuM.GhoshA.KumarJ.SinghS.Palit SagarS. (2022). Weibull parameter-based probability distribution for predicting creep life of power plant materials: a non-destructive approach. Mater. today59, 1111–1118. 10.1016/j.matpr.2022.03.024
19
SavastanoS.SantosS.RadonjicM.SoboyejoW. (2009). Fracture and fatigue of natural fiber-reinforced cementitious composites. Cem. Concr. Compos.31 (4), 232–243. 10.1016/j.cemconcomp.2009.02.006
20
ShaoYiClaudiaP. (2022). Ostertag, Bond-slip behavior of steel reinforced UHPC under flexure: experiment and Prediction. Cem. Concr. Compos. (133), 104724. 10.1016/j.cemconcomp.2022.104724
21
SinghS. P.KaushikS. K. (2001). Flexural fatigue analysis of steel fiber-reinforced concrete. ACI Mater J.98, 306–312.
22
SinghS. P.KaushikS. K. (2003). Fatigue strength of steel fibre reinforced concrete in flexure. Cem. & Concr. Compos.25, 779–786. 10.1016/s0958-9465(02)00102-6
23
SinghS. P.MohammadiY.KaushikS. K., Flexural fatigue analysis of steel fibrous concrete containing mixed fibers. ACI Mater J., 2005b,102:438–444.
24
SinghS. P.SinghB.KaushikS. K. (2005a). Probability of fatigue failure of steel fibrous concrete. Mag. Concr. Res.57, 65–72. 10.1680/macr.57.2.65.58701
25
SinghS. P.SinghB.KaushikS. K. (2005c). Probability of fatigue failure of steel fibrous concrete. Mag. Concr. Res.57, 65–72. 10.1680/macr.57.2.65.58701
26
SongQ. L.YuR.ShuiZ. H.WangX. P.RaoS. D.LinZ. W. (2018). Optimization of fiber orientation and distribution for a sustainable Ultra-High-Performance Fiber Reinforced Concrete (UHPFRC): experiments and mechanism analysis. Constr. Build. Mater.169, 8–19.
27
TengL.MengW.KhayatK. H. (2020). Rheology control of ultra-high-performance concrete made with different fiber contents. Cem. Concre. Res.138, 106222. 10.1016/j.cemconres.2020.106222
28
ThaiD.-K.NguyenD.-L.PhamT.-H.DoanQ. H. (2021). Prediction of residual strength of FRC columns under blast loading using the FEM method and regression approach. Constr. Build. Mater.276, 122253. 10.1016/j.conbuildmat.2021.122253
29
VesicA. S. (1969). Analysis of structural behavior of road test rigid pavements. Highw. Res., 291–294.
30
WuZ.ZhaoG.HuangC. (1995). Concrete fatigue fracture characteristics. J. Civ. Eng.28 (3), 59–65.
31
YinT.YuR.LiuK.WangZ.FanD.WangS.et al (2022). Precise mix-design of Ultra-High-Performance Concrete (UHPC) based on physicochemical packing method: from the perspective of cement hydration. Constr. Build. Mater.352, 128944. 10.1016/j.conbuildmat.2022.128944
32
ZhangB.YuJ.ChenW.SunH.ChenS.WangH. (2022). Interfacial properties between ultra-high performance concrete (UHPC) and steel: from static performance to fatigue behavior. Eng. Struct.273, 115145. 10.1016/j.engstruct.2022.115145
33
ZhaoB.GengC.SongZ.PanJ.ChenJ.XiaoP.et al (2024). A fatigue fracture phase field model considering the effect of steel fibers in UHPC. Eng. Fract. Mech.300, 109981. 10.1016/j.engfracmech.2024.109981
34
ZhouBoUchidaY. (2017). Influence of flowability, casting time and formwork geometry on fiber orientation and mechanical properties of UHPFRC. Cem. Concr. Res.95, 164–177. 10.1016/j.cemconres.2017.02.017
Summary
Keywords
fatigue life, UHPC, weibull distribution, probability failure, fatigue failure probability
Citation
Luo C, Yang P, Niu Y, Zhang Y and Cheng C (2024) Analytical method of incorporating failure probability to predict the fatigue life of ultra-high-performance concrete (UHPC). Front. Built Environ. 10:1403245. doi: 10.3389/fbuil.2024.1403245
Received
19 March 2024
Accepted
28 June 2024
Published
09 September 2024
Volume
10 - 2024
Edited by
Yanping Zhu, University of Mount Union, United States
Reviewed by
Peng Zhang, Zhengzhou University, China
Md. Mashfiqul Islam, Florida State University, United States
Updates
Copyright
© 2024 Luo, Yang, Niu, Zhang and Cheng.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yanfei Niu, yanfei.n@gzhu.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.