Optimization of Ultra-High Performance Concrete Based on Response Surface Methodology and NSGA-II
- 1. School of Materials Science and Engineering, University of Jinan, China
- 2. Tiezheng Testing Technology Co., Ltd, China
ABSTRACT
This study systematically investigated three influential factors—water-to-binder ratio, cement-sand ratio, and steel fiber content—that significantly impact the performance of ultra-high-performance concrete (UHPC). Utilizing the Response Surface Methodology (RSM) with a Central Composite Design (CCD), 20 carefully designed mix proportions underwent comprehensive experimental testing. Through rigorous statistical analysis, models were established to elucidate the complex relationships between the specified factors and the overall properties of UHPC. Variance analysis reveals significant effects of the three factors on UHPC performance, with workability and compressive strength increasing with higher cement-sand ratios while flexural strength decreases. Moreover, increased water-to-binder ratios exhibit substantial negative impacts on both 28-day compressive and flexural strengths. Despite adversely affecting workability, higher steel fiber dosages contribute positively to mechanical performance. Furthermore, Monte Carlo sampling and the multi-objective non-dominated sorting genetic algorithm-II (NSGA-II) were employed to validate the reliability of the statistical model and to conduct multi-objective optimization. The final UHPC mix design obtained consists of a cement-sand ratio of 1.12, a water-binder ratio of 0.16, and a steel fiber content of 2.94%. Experimental results yielded a slump flow of 802 mm, compressive strength of 122.7 MPa, and flexural strength of 24.3 MPa.
KEYWORDS
- Response Surface Methodology
- UHPC
- Mechanical performance
- Workability
- NSGA-II
CITATION
Ye Z, Wu J, Wang Z, Su L, Gao Z (2024) Optimization of Ultra-High Performance Concrete Based on Response Surface Methodology and NSGA-II. J Materials Applied Sci 5(1): 1009.
INTRODUCTION
Since its inception, Ultra-High-Performance Concrete (UHPC) [1], has witnessed a growing significance in the realm of construction engineering, evolving over time. Notably, the realm of construction has articulated a demand for the high-performance characteristics of building materials to effectively address challenges posed by adverse conditions [2,3]. The extensive utilization of UHPC in practical engineering is attributable to its remarkable mechanical properties and exemplary durability, including its capacity to withstand the deleterious effects of substance erosion [4,5].
In comparison to ordinary concrete, UHPC utilizes high cement content (800-1000 kg/m³) [6], with the dosage of silica fume accounting for 25-35% of the cement content [7]. Fine aggregates ranging from 120-800 μm are chosen, and UHPC typically employs a low water-to-binder ratio. The use of efficient superplasticizers ensures workability [8], and 13mm long, 0.2mm diameter steel fibers are added. The mix design is conducted based on the closely packed theory to achieve the densest state and consequently, an ultra-low porosity, leading to enhanced performance [9-14]. The water-to-binder ratio, cement-sand ratio, and steel fiber content are critical factors exerting significant influence on the mechanical properties of UHPC [15-17].
The Response Surface Method (RSM) is an analytical approach that combines mathematical and statistical methods for optimizing multi-factor parameters. In contrast to single- factor and orthogonal experimental methods, RSM considers interactions among multiple factors to predict the optimal response values. The optimal results obtained from RSM are not limited to a specific set of experiments but rather represent the optimal points on the entire three-dimensional surface. RSM analyzes the interactions between major influencing factors by fitting models and establishes specific functional relationships between influencing factors and response values, thereby enabling discussions on the extent of the impact of multiple factors on response values. Consequently, it finds extensive application in the fields of multi-parameter optimization and material design [18-21]. In recent years, with the rapid and continuous development of the concrete industry, RSM has found widespread use in the optimization of concrete configurations for multiple performance requirements. It achieves this by adjusting certain nonlinear factors to define an optimal domain [22]. Many scholars have utilized RSM to study the optimization of UHPC.
Ferdosian et al., [23] employing the Response Surface Methodology (RSM), investigated the impact of sand, silica fume, and ultrafine fly ash on the primary properties of UHPC. They utilized multi-objective optimization to obtain the optimal mix proportions for enhanced performance. Wille et al. [24], studied the influence of parameters such as cement, silica fume, and water-to-binder ratio on UHPC performance, ultimately determining the optimal range for these parameters. Mohammed et al. [25], explored the effects of different water-to-binder ratios on UHPC performance, revealing that lower ratios result in higher mechanical performance. Steel fibers exert a significant influence on the toughness and flexural performance of UHPC, with current research indicating that the highest acceptable fiber content is 6% [26]. The type of steel fiber also imparts varying effects on UHPC, for a given content, shaped steel fibers exhibit a greater impact on tensile and flexural strength than straight steel fibers [27]. Additionally, different lengths of steel fibers contribute differently to UHPC performance [28]. Mu et al. [29], applied a magnetic field to orient the steel fibers in UHPC, enhancing its mechanical properties. Wang et al. [30], adjusted the water-to-binder ratio or the dosage of superplasticizers to achieve different rheological properties, studying the influence of rheological properties on fiber distribution and yield stress.
Due to the complexity of UHPC compositions, minor perturbations can result in significant performance deviations. Currently, multiple experiments are typically conducted to reduce experimental errors, leading to a substantial increase in the number of experiments required. NSGA-II is a commonly used multi-objective optimization algorithm. It simulates the process of natural evolution to find the optimal solution set while maintaining a balance between different objectives. In recent years, the powerful optimization capability of NSGA-II has led to its widespread application in the research of cement-based materials [31-33].
Many researchers have extensively studied the influencing factors of UHPC using Response Surface Methodology (RSM). However, there is a lack of systematic research utilizing both RSM and NSGA-II algorithms for optimizing UHPC mix proportions. This study primarily focuses on exploring the collective effects of the water-to-binder ratio, cement-sand ratio, and steel fiber content. The research objectives involve optimizing these factors. using the Central Composite Design (CCD) method within the framework of Response Surface Methodology, with workability, compressive strength, and flexural strength selected as response variables. Regression analysis, residual analysis, and model fitting tests were employed for a thorough analysis of each factor, combining the NSGA-II method for multi-objective optimization, the ultimate aim was to derive a well-optimized mix proportion for UHPC.
EXPERIMENTAL STUDY
Materials
In this study, the raw materials used include PO 52.5 Portland cement, microspheres, silica fume (SF), expanding agent, steel fiber, quartz sand, and a polycarboxylate acid water reducer. The chosen fibers are straight steel fibers with a diameter of 0.2mm and a length of 1.3mm. Two grades of quartz sand, namely 20-40 mesh and 40-80 mesh were used. Table 1 provides the chemical compositions of cement, microspheres, and silica fume.
Table 1. Chemical composition of cementitious materials (%).
Code |
CaO |
SiO2 |
Al2O3 |
Fe2O3 |
SO3 |
K2O |
Na2O |
MgO |
L.OL |
Cement |
63.06 |
21.01 |
3.76 |
3.34 |
3.63 |
0.89 |
0.19 |
3.03 |
2.1 |
Silica fume |
0.35 |
97.43 |
0.42 |
0.13 |
0.33 |
0.16 |
0.33 |
0.22 |
1.2 |
microspheres |
7.9 |
45.28 |
39.76 |
2.9 |
0.47 |
0.91 |
0.42 |
0.18 |
1.8 |
Specimens Preparation
The preparation of test specimens involves the following steps: initially, mix cement and aggregates, such as quartz sand, for 5 minutes. Subsequently, add the high-efficiency water- reducing agent and water, and stir for 10 minutes. Once the mixture reaches an appropriate consistency, uniformly disperse the steel fibers by adding them and stirring for an additional 5 minutes. After the mixing process is completed, conduct slump flow tests and mold specimens according to the T/CECS 864 standard.
Experimental Design and RSM
The RSM statistical method can achieve optimized results for target performance under given influencing conditions and can significantly minimize the number of experiments required. In RSM, CCD is one of the most commonly used optimization methods in the field of concrete. CCD includes a 2k factorial design, a central point K0 (3 ≤ K0 ≤ 5), and a star design method. The central point design is considered a preferable alternative in full factorial design, with the distance from test points to the central point denoted as α. The parameters α and K0 are designed to satisfy rotatability. In this study, the CCD method in RSM was adopted for design and analysis, and the first-order regression model used is [34]:
In this context, where Y represents the anticipated response value, β0 is the intercept term, βi is the linear term, Xi is an independent variable, K is the number of variables, and ε is the associated experimental error.
The second-order polynomial model is given by:
In this context, where Y represents the anticipated response value, β? is the intercept term, βi is the linear term, βii is the quadratic term, βij is the quadratic interaction term, and Xi and Xj are independent variables.
As shown in Table 2, in this study, the variables selected are cement-sand ratio (coded as A), water-to-binder ratio (coded as B), and steel fiber content (coded as C). The response values are slump flow, compressive strength, and flexural strength, as indicated in Table 3.
Table 2. Range and Codes of Variables in the CCD Method.
Variables |
Coded |
Actual |
||
-1 |
0 |
1 |
||
cement-sand ratio |
A |
1 |
1.1 |
1.2 |
Water-to-binder ratio |
B |
0.14 |
0.16 |
0.18 |
Steel Fiber Content (%) |
C |
2 |
2.5 |
3 |
Table 3. Mix Design in CCD Method.
No |
Actual |
C |
SF |
FA |
Sand |
SP |
Steel Fiber |
w |
|||
|
A |
B |
C |
20-40 |
40-80 |
||||||
1 |
1.2 |
0.18 |
3 |
865 |
110 |
110 |
275 |
629 |
28 |
234 |
195.3 |
2 |
1.1 |
0.16 |
1.66 |
865 |
110 |
110 |
300 |
687 |
28 |
129.5 |
173.6 |
3 |
1.2 |
0.18 |
2 |
865 |
110 |
110 |
275 |
629 |
28 |
156 |
195.3 |
4 |
1.1 |
0.16 |
2.5 |
865 |
110 |
110 |
300 |
687 |
28 |
195 |
173.6 |
5 |
1 |
0.14 |
3 |
865 |
110 |
110 |
330 |
755 |
28 |
234 |
151.9 |
6 |
1.1 |
0.13 |
2.5 |
865 |
110 |
110 |
300 |
687 |
28 |
195 |
141.1 |
7 |
1 |
0.14 |
2 |
865 |
110 |
110 |
330 |
755 |
28 |
156 |
151.9 |
8 |
1.1 |
0.19 |
2.5 |
865 |
110 |
110 |
300 |
687 |
28 |
195 |
206.2 |
9 |
1.1 |
0.16 |
2.5 |
865 |
110 |
110 |
300 |
687 |
28 |
195 |
173.6 |
10 |
1.2 |
0.14 |
3 |
865 |
110 |
110 |
275 |
629 |
28 |
234 |
151.9 |
11 |
1.1 |
0.16 |
2.5 |
865 |
110 |
110 |
300 |
687 |
28 |
195 |
173.6 |
12 |
1.27 |
0.16 |
2.5 |
865 |
110 |
110 |
260 |
595 |
28 |
195 |
173.6 |
13 |
0.93 |
0.16 |
2.5 |
865 |
110 |
110 |
355 |
812 |
28 |
195 |
173.6 |
14 |
1 |
0.18 |
2 |
865 |
110 |
110 |
330 |
755 |
28 |
156 |
195.3 |
15 |
1.1 |
0.16 |
2.5 |
865 |
110 |
110 |
300 |
687 |
28 |
195 |
173.6 |
16 |
1.1 |
0.16 |
3.34 |
865 |
110 |
110 |
300 |
687 |
28 |
160.5 |
173.6 |
17 |
1.2 |
0.14 |
2 |
865 |
110 |
110 |
275 |
629 |
28 |
156 |
151.9 |
18 |
1.1 |
0.16 |
2.5 |
865 |
110 |
110 |
300 |
687 |
28 |
195 |
173.6 |
19 |
1 |
0.18 |
3 |
865 |
110 |
110 |
330 |
755 |
28 |
234 |
195.3 |
20 |
1.1 |
0.16 |
2.5 |
865 |
110 |
110 |
300 |
687 |
28 |
195 |
173.6 |
Experimental Methods
Workability: Following the Chinese T/CECS 864 Standard for test method of ultra-high performance concrete, the workability performance of freshly mixed UHPC is evaluated using the slump flow test. A slump cone is filled in one continuous operation, and excess material is struck off level with the top of the cone. The cone is then smoothly and vertically lifted, and the diameter of the spread of UHPC is measured in two perpendicular directions using a steel ruler along the edge of the cone.
Compressive Strength: The compressive strength test, following the Chinese T/CECS 864 Standard for test method of ultra-high performance concrete, involves casting the specimens into molds with dimensions of 100×100×100 mm³. Immediately after casting, the specimens are covered with plastic film and cured for 24 hours at 20±2°C. After demolding, the specimens are further cured at a temperature of 20±2°C and a relative humidity exceeding 95%. The compressive strength test for UHPC cubic specimens is conducted at a curing age of 28 days, with six specimens per batch.
Flexural Strength: Following the Chinese T/CECS 864 Standard for test method of ultra-high performance concrete, form 100×100×400 mm³ prismatic specimens. After demolding, place the specimens under curing conditions at a temperature of 20±2°C and a relative humidity exceeding 95%. When the curing period reaches 28 days, conduct the flexural strength test, with three specimens per batch.
RESULTS AND ANALYSIS
In this study, the CCD model consists of three factors with three levels each, and three response values: slump flow, compressive strength at 28 days, and flexural strength at 28 days. The summarized test results are presented in Table 4. The compressive strength data represent the arithmetic mean of six specimen measurements, meeting the standard requirement that the individual test values deviate no more than 10% from the mean. The flexural strength data are the arithmetic mean of three measurements, with deviations not exceeding the 15% standard requirement. Additionally, the relative deviations between the experimental data and the predicted values are also listed in Table 4.
Table 4. Experimental Results and Predicted Values of UHPC.
Sample |
Workability(mm) |
|
28d-Compressive strength(MPa) |
|
28d-Flexural strength(MP) |
|
|||
Actual |
Predicted |
RD(%) |
Actual |
Predicted |
RD(%) |
Actual |
Predicted |
RD(%) |
|
1 |
900 |
875.79 |
2.76 |
112.6 |
111.35 |
1.12 |
21.59 |
20.66 |
4.5 |
2 |
860 |
868.06 |
-0.93 |
116.6 |
116.15 |
0.39 |
23.17 |
20.78 |
11.5 |
3 |
930 |
950.81 |
-2.19 |
110.5 |
111.6 |
-0.99 |
18.53 |
17.88 |
3.64 |
4 |
720 |
752.26 |
-4.29 |
114.3 |
113.29 |
0.89 |
24.68 |
23.11 |
6.79 |
5 |
670 |
645.94 |
3.72 |
124.6 |
123.07 |
1.24 |
27.51 |
28.34 |
-2.93 |
6 |
700 |
720.38 |
-2.83 |
122 |
122.07 |
-0.06 |
30.55 |
27.54 |
10.93 |
7 |
720 |
720.96 |
-0.13 |
116.1 |
115.52 |
0.5 |
24.02 |
25.57 |
-6.06 |
8 |
920 |
939.58 |
-2.08 |
101.1 |
104.5 |
-3.25 |
16.33 |
18.69 |
-12.63 |
9 |
740 |
752.26 |
-1.63 |
111.2 |
113.29 |
-1.84 |
23.54 |
23.11 |
1.86 |
10 |
740 |
745.45 |
-0.73 |
118.6 |
121.8 |
-2.63 |
23.09 |
25.92 |
-10.92 |
11 |
760 |
752.26 |
1.03 |
114 |
113.29 |
0.63 |
24.85 |
23.11 |
7.53 |
12 |
840 |
835.95 |
0.48 |
120.1 |
115.5 |
3.98 |
21.5 |
21.07 |
2.04 |
13 |
650 |
668.58 |
-2.78 |
110 |
111.07 |
-0.96 |
24.93 |
25.15 |
-0.87 |
14 |
870 |
851.3 |
2.2 |
106.5 |
105.07 |
1.36 |
19.68 |
20.31 |
-3.1 |
15 |
750 |
752.26 |
-0.3 |
112.2 |
113.29 |
-0.96 |
21.54 |
23.11 |
-6.79 |
16 |
710 |
741.89 |
-4.3 |
121.8 |
122.29 |
-0.4 |
26.84 |
25.45 |
5.46 |
17 |
860 |
820.47 |
4.82 |
120.4 |
122.05 |
-1.35 |
20.54 |
23.15 |
-11.27 |
18 |
780 |
752.26 |
3.69 |
116 |
113.29 |
2.39 |
23.99 |
23.11 |
3.81 |
19 |
810 |
776.28 |
4.34 |
113.9 |
112.62 |
1.14 |
23.34 |
23.08 |
1.13 |
20 |
745 |
752.26 |
-0.97 |
111.9 |
113.29 |
-1.23 |
22.06 |
23.11 |
-4.54 |
The relative deviation of the slump flow and 28- day compressive strength is less than 5%. Although the relative deviation for flexural strength is larger due to its smaller values, it remains below 15%. In Table 5, the results of the regression model and the variance analysis are presented. Parameters associated with significance, such as the F-value, P-value, R-squared (R2), and Adjusted R-squared (Adj-R2), are outlined in the ANOVA table (factors with excessively high P-values have been removed in Table 5. A P-value less than 0.05 indicates that the model term is significant; conversely, a P-value greater than 0.1 indicates that the model term is not significant, with all model P-values being <0.0001. Utilizing regression analysis, the final model is derived.
Table 5. Results of the Full Regression Model and Variance Analysis from Numerical Simulation.
ANOVA |
Responses |
|||||
Slump flow |
P-value |
28d-Compressive strength(MPa) |
P-value |
28d-Flexural strength(MPa) |
P-value |
|
intercept |
752.26 |
|
113.29 |
|
23.11 |
|
Linear terms |
|
|
|
|
||
A |
49.76 |
<0.0001 |
1.32 |
0.051 |
-1.21 |
0.0245 |
B |
65.17 |
<0.0001 |
-5.22 |
<0.0001 |
-2.63 |
<0.0001 |
C |
-37.51 |
<0.0001 |
1.83 |
0.0103 |
1.39 |
0.0119 |
Quadratic terms |
|
|
|
|
||
A2 |
|
|
|
|
|
|
B2 |
27.48 |
0.001 |
|
|
|
|
C2 |
18.64 |
0.0141 |
2.1 |
0.0034 |
|
|
Interaction terms |
|
|
|
|
||
AB |
|
|
|
|
|
|
AC |
|
|
-1.95 |
0.0298 |
|
|
BC |
|
|
|
|
|
|
Other terms |
|
|
|
|
|
|
R2 |
0.9333 |
|
0.8806 |
|
0.7298 |
|
Adj-R2 |
0.9094 |
|
0.838 |
|
0.6792 |
|
Pred-R2 |
0.8245 |
|
0.7531 |
|
0.5405 |
|
F-Value |
20.66 |
|
20.66 |
|
14.41 |
|
Lack of Fit |
1.92 |
|
1.92 |
|
2.08 |
|
Model p-Value |
<0.0001 |
|
<0.0001 |
|
<0.0001 |
|
Analysis of Working Performance
In Table 4, the actual and predicted values of UHPC workability analyzed through RSM are presented. The non-significant data for Lack of Fit (F=4.71) suggests that the model is suitable for predictive analysis (R2=0.9333; Adj. R2=0.9094). Figure 1(a), illustrates the normal distribution plot of standardized residuals for slump flow. The residuals closely align with the straight line, indicating a satisfactory fit of the model. According to the ANOVA results in Table 5, the relationship between workability and the three factors has been determined. The ANOVA model used is as follows: Slump flow of UHPC mm = 752.26+49.76*A+65.17*B- 37.51*C+27.48*B2+18.64*C2.
(Where A represents the cement-sand ratio, B represents the water-to-binder ratio, and C represents the steel fiber content.)
As depicted in Figure 1(b), the component effects curve illustrates the variation trend of model-predicted response values with changes in variables. The impact of factor A on slump flow is linear, while the effects of factors B and C are nonlinear. With an increment in the cement-sand ratio and water-to-binder ratio, the slump flow of freshly mixed UHPC increases. The rise in the water-to-binder ratio leads to increased moisture surrounding the cement particles, resulting in an augmented water film thickness and reduced inter-particle friction, consequently enhancing the slump flow [35]. From Figure 2, it can be inferred that, with a steel fiber content of 3.0% and a cement-sand ratio of 1.2, the maximum increase in slump flow is 160mm with an increasing water-to-binder ratio.
Due to the mechanical fragmentation of quartz sand, its particles exhibit numerous sharp edges and irregular shapes, impeding the flow of UHPC. With an increase in the cement-sand ratio, the volume of quartz sand per unit of UHPC decreases, mitigating this hindrance and resulting in an augmented slump flow. From Figure 1(b) and Figure 3, the impact of factor C, steel fiber content, on slump flow is detrimental. As the fiber content increases, slump flow decreases, with a maximum reduction of up to 120mm. This phenomenon arises from the crisscross arrangement of steel fibers within the UHPC, forming a distinct network that impedes the flow. Greater steel fiber content leads to a denser network with smaller mesh sizes, thereby intensifying the hindering effect.
Figure 1 Normal plot of residuals value (a) and perturbation curve of Slump Flow (b).
Figure 2 Slump flow response surface with water-to-binder ratio: (a) Low water-to-binder ratio(0.14), (b) High water-to-binder ratio(0.18).
Figure 3 Slump flow response surface with steel fiber content: (a) Low steel fiber content(2%), (b) High steel fiber content(3%).
Analysis of Compressive Strength Performance
In Table 4, the actual and predicted values of UHPC’s 28- day compressive strength are presented through RSM analysis. The non-significant data for Lack of Fit (F=1.92) suggests the applicability of the model for predictive analysis (R2=0.8806; Adj. R2=0.8380). Figure 4(a) depicts the normal distribution plot of standardized residuals for the 28-day compressive strength. The proximity of residual values to the straight line indicates favorable results. Based on the ANOVA results in Table 5, the relationship between workability and the three factors is established. The ANOVA model employed is as follows:
The 28-day compressive strength of UHPC MPa = 113.29+1.32*A-5.22*B+1.83*C -1.95*AC+2.10*C2.
(The P value of factor A in the model is 0.0510, which is less than the predetermined significance level, but considering the small difference between P value and 0.05, it still has certain statistical significance, so it is not ignored)
As shown in Figure 4(b), the component effect curves illustrate the trend of the predicted response values with varying variables. The impact of factors A and B on compressive strength is linear, while the influence of C on compressive strength is nonlinear. With an increase in the cement-sand ratio, the 28-day compressive strength exhibits a linear improvement. This result may be related to the Interfacial Transition Zone (ITZ) around the aggregate, which is considered a weak area in UHPC. An increase in the cement-sand ratio leads to a decrease in the volume fraction of aggregates. As the amount of aggregates in UHPC decreases, the ITZ formed in the cured UHPC reduces. Although the reduction in aggregates weakens the compressive strength of UHPC, the decrease in the ITZ enhances the compressive strength. These two effects counteract and combine, increasing compressive strength [36].
Figure 4 Normal plot of residual value(a) and perturbation curve of 28-day compressive strength(b).
With an increase in the water-to-binder ratio, the compressive strength of all UHPC specimens shows a significant decrease trend. This is consistent with the research results of Appa Rao [37], Mohamad, et al., [38], and Fernandes et al., [39]. An increase in the water-to-binder ratio leads to a thicker water film on the particle surfaces, resulting in more pores in the cured UHPC, increasing the porosity, and weakening the ITZ. As shown in Figure 5, with a constant steel fiber content and cement-sand ratio, an increase in the water-to-binder ratio causes the 28-day compressive strength to decrease from 124.6 MPa to 113.9 MPa, a reduction of 10.7 MPa.
Figure 5 Response Surface of 28-day Compressive Strength as a Function of water-to-binder Ratio: (a) Low water-to-binder Ratio(0.14), (b) High water-to-binder Ratio(0.18).
The influence of steel fibers on the 28-day compressive strength shows a decreasing-then-increasing trend with an increase in the content. The lowest point occurs at a content of 2.50%, as shown in (Figure 6). This result may be attributed to the hindrance effect of steel fibers on the flow of UHPC, leading to a decrease in compactness, thereby reducing strength when the fiber content increases from 2.0% to 2.5% [40]. When the content increases from 2.5% to 3.0%, the increase in strength due to steel fibers outweighs the impact of voids, increasing compressive strength.
Figure 6 Response Surface of 28-day Compressive Strength as a Function of Steel Fiber Content: (a) Low Steel Fiber Content(2%), (b) Medium Steel Fiber Content(2.5%), (c) High Steel Fiber Content(3%).
The Flexural Strength Performance Analysis
In Table 4, the actual values and predicted values of UHPC flexural strength through Response Surface Methodology (RSM) analysis are presented. The non-significant Lack of Fit data (F=2.08) indicates the suitability of the model for predictive analysis (R2=0.7298; Adj.R2=0.6792). Figure 7(a), illustrates the normal distribution of standardized residuals for flexural strength. The residuals closely align with the straight line, indicating satisfactory results. Based on the ANOVA results in Table 5, the relationship between flexural strength and the three factors is established. A linear relationship is observed between flexural strength and the three factors, and the ANOVA model used is as follows:
Flexural Strength of UHPC MPa =23.11-1.21*A-2.63*B+1.39*C. As shown in Figure 7(b), the component effect curves illustrate the trend of predicted response values with changes in components. The influence of the three factors on flexural strength is linear. Flexural strength exhibits a negative correlation with factors A and B, where the cement-sand ratio and water- to-binder ratio gradually increase from low levels (coded as -1) to high levels (coded as +1), resulting in a decrease in flexural strength. An increase in the water-to-binder ratio leads to more internal pore structures and, a decrease in the strength of the ITZ. The increase in the cement-sand ratio leads to a decrease in flexural strength, which can be explained by the crack bridging mechanism [41-42]. The interlocking of sand particles inside the specimen significantly consumes the energy of the crack surface. When the crack encounters materials similar to sand during the expansion process, it cannot pass directly through the sand particles but has to transfer energy along the transition zone on the surface of the sand particles. This increases the energy transfer path, requiring the consumption of more energy. Therefore, an increase in the cement-sand ratio leads to a decrease in sand content, resulting in a reduction in flexural strength.
Figure 7 Normal plot of residual value(a) and perturbation curve(b) of Flexural Strength.
The impact of factor C, in contrast, exhibits an opposing trend to A and B, showing a positive correlation between flexural strength and steel fiber content. The increasing content of steel fiber results in smaller inter-fiber spacing and denser microstructures near the fibers. Consequently, a better bond strength forms between the fibers and the matrix [43]. When cracks occur, more steel fibers intersect each crack, increasing the bonded area between the matrix and the fibers. This impedes the extension and propagation of cracks, leading to the formation of more microcracks and delaying the growth of macroscopic cracks [44,45]. As a result, a higher initial cracking strength is achieved, leading to higher flexural strength of UHPC. Consequently, achieving higher flexural strength in UHPC can be realized by introducing a greater concentration of steel fibers. This observation is consistent with the pattern illustrated in Figure 8, where flexural strength rises with an augmentation in steel fiber content while keeping other variables constant. For instance, in comparison between UHPC samples with 2% and 3% steel fiber content (cement-sand ratio of 1.0, water-to-binder ratio of 0.14), the flexural strength increases from 24.02 MPa to 27.51 MPa, representing a notable increase of 3.49 MPa.
Figure 8 Response Surface of 28-day Flexural Strength with Steel Fiber Content: (a) Low Steel Fiber Content, (b) High Steel Fiber Content.
Optimization Validation
Constraints in Multi-Objective Optimization: The main purpose of multi-objective optimization is to minimize the mean square deviation while satisfying constraints and maximizing the reliability of the optimization results. Optimization methods effectively transform deterministic optimization into reliability- based optimization to find the optimal design space between the target performance and reliability probability of UHPC. Therefore, the design variables, optimization objectives, and constraints for the multi-objective optimization design based on robustness for UHPC are determined as follows:
Design Variables: Water-to-binder ratio, cement-sand ratio, and steel fiber content. Optimization Objectives: (1) Slump flow, (2) Maximum 28-day compressive strength, (3) Maximum 28-day flexural strength. Constraint Conditions: Deterministic constraints for slump flow range from 750mm to 850mm, for 28-day compressive strength range from 115 MPa to 125 MPa, and for 28-day flexural strength range from 21 MPa to 31 MPa. Therefore, based on the design variables, optimization objectives, and constraint ranges, the multi-objective optimization model is determined as follows:
Reliability Analysis: The reliability analysis was conducted using the Monte Carlo sampling method. Random factors such as noise are selected, and 10,000 random samples are drawn to determine their distribution characteristics, including mean, standard deviation, and Sigma levels. As shown in Figure 9, the Sigma levels for slump flow, 28-day compressive strength, and 28-day flexural strength were determined to be 0.99, 1.125, and 0.99, respectively, with reliabilities of 67.93%, 86.94%, and 83.27%, all of which are below 90%. This indicates that the reliability of results within the deterministic constraint range based on the response surface model is relatively low.
Figure 9 The Monte Carlo sampling analysis results.
Pareto Frontier of Response Variables and Reliability Analysis: Based on the mathematical models in Sections 3.1, 3.2, and 3.3, this study employed NSGA-II for multi-objective optimization. During optimization, the initial population size of NSGA-II was set to 200, with 200 iterations. The crossover distribution index was set to 0.8, mutation percentage to 0.7, and mutation probability to 0.4. NSGA-II operates by searching for the optimal Pareto solutions in the solution space and then performs Pareto front ranking to select the best solutions. Concurrently, to validate the reliability of the NSGA-II multi- objective optimization, reliability analysis was conducted on the sample data to obtain the mean, variance, and 6-sigma level, as shown in Figure 10. It can be observed that the sigma level values of the selected sample points for UHPC multi-objective optimization are all not less than 6. This indicates that the reliability of multi-objective optimization is not less than 99.99%, satisfying the reliability requirements of 6-sigma robust multi- objective optimization principles.
Figure 10 Reliability requirements of NSGA-II.
Multi-Objective Optimization Result Verification: As shown in Table 6, based on the NSGA-II multi-objective optimization results, the most ideal point was selected for validation experiments, with a cement-sand ratio of 1.12, a water- to-binder ratio of 0.16, and a steel fiber content of 2.94. This scheme was subjected to three repeated validation experiments. From the validation experiment results, it can be observed that there is a slight deviation between the experimental results and the NSGA-II optimization results. Although there is some deviation in the fitted model, the validation experiment results partially demonstrate the rationality of the optimization method.
Table 6. The validation experiment results of NSGA-II.
Response quantity |
Slump flow/mm |
Relative error |
28d-Compressive strength/ MPa |
Relative error |
28d-Flexural strength/MPa |
Relative error |
Predicted value |
815 |
— |
123.6 |
— |
24.4 |
— |
|
|
|
|
|
|
|
Test 1 |
800 |
-1.84% |
122.3 |
-1.05% |
25.2 |
3.28% |
Test 2 |
795 |
-2.45% |
124.1 |
0.40% |
24.1 |
-1.23% |
Test 3 |
810 |
-0.60% |
121.7 |
-1.54% |
23.6 |
-4.51% |
average |
802 |
-1.60% |
122.7 |
-0.73% |
24.3 |
-0.41% |
SUMMARY AND DISCUSSION
This study employed the response surface methodology to analyze the influence of water-to-binder ratio, cement- sand ratio, and steel fiber content on the workability, 28- day compressive strength, and 28-day flexural strength of UHPC. The study also conducted multi-objective optimization. Polynomial models for workability, compressive strength, and flexural strength were established. The fitting effect of the models was found to be satisfactory through various correlation coefficients and predicted values. Simultaneously, NSGA-II was employed for multi-objective optimization, and experimental validation confirmed that the optimized results closely matched the validation results. The water-to-binder ratio exhibited a negative impact on 28-day compressive strength and flexural strength, with a more significant effect as the water-to-binder ratio increased, leading to a gradual reduction in mechanical properties. The cement-sand ratio demonstrated a good linear correlation with the three response variables. Specifically, workability and compressive strength increased with an increase in the cement-sand ratio, while flexural strength decreased. The study revealed that using a high dosage of steel fibers, while negatively affecting workability, proved beneficial for mechanical properties, especially achieving higher flexural strength. The addition of 3% steel fibers was considered the most suitable dosage in UHPC, providing higher mechanical performance while maintaining adequate workability.
ACKNOWLEDGMENTS
This work was supported by the Shandong Province Key R&D Program (Majoy Technological Innovation) (2021CXGC010301) and the National Natural Science Foundation of China (No. 52172018, 52302021).
REFERENCES
- De Larrard F, Sedran T. Optimization of ultra-high-performance concrete by the use of a packing model. J of CCR. 1994; 24: 997-1009.
- Jiang Du, Weina Meng, Kamal HK, Yi Bao, Pengwei G, Zhenghua Lyu, et al. New development of ultra-high-performance concrete (UHPC). Compos Part B Eng. 2021; 224: 109220.
- Mugahed A, Shan-Shan H, Ali MO, Natt Makul, Hakim S. Abdelgader, et al. Recent trends in ultra-high performance concrete (UHPC): Current status, challenges, and future prospects. Const and Buil Materials. 2022; 352: 129029.
- Russell HG, Graybeal BA. Ultra-high performance concrete: A state- of-the-art report for the bridge community. United States. Federal Highway Administration. Office of Infrastructure Research and Development. 2013.
- Habel K, Gauvreau P. Response of ultra-high performance fiber reinforced concrete (UHPFRC) to impact and static loading. Cem Concr Compos. 2008; 30: 938-946.
- Amin M, Zeyad AM, Tayeh BA, Ibrahim SA. Effects of nano cotton stalk and palm leaf ashes on ultrahigh-performance concrete properties incorporating recycled concrete aggregates. Const and Buil Materials. 2021; 302: 124196.
- Graybeal BA. Material property characterization of ultra- high performance concrete. United States. Federal Highway Administration. Office of Infrastructure Research and Development. 2006.
- Hirschi T, Wombacher F. Influence of different superplasticizers on UHPC[C]. Proceedings of the Second International Symposium on Ultra High Performance Concrete, Kassel University Press, Kassel. 2008; 77-84.
- Yang Du, Lifeng Z, Shaoqin R, Xiaoqian Q, Kuangliang Q, Pengyun C, et al. Powder gradation effect on the fresh, mechanical, and early-age shrinkage properties of UHPC. Build Eng. 2023; 67: 105958.
- Yin T, Liu K, Fan D, Rui Yu. Derivation and verification of multilevel particle packing model for Ultra-High Performance Concrete (UHPC): Modelling and experiments. Cem Concr Compos. 2023; 136: 104889.
- Xinpeng Wang , Di Wu, Jinrui Zhang, Rui Yu, Dongshuai Hou, Zhonghe Shui. Design of sustainable ultra-high performance concrete: A review. Const and Buil Materials. 2021: 307: 124643.
- Bahmani H, Mostofinejad D. Microstructure of ultra-high- performance concrete (UHPC)–A review study. J Build Eng. 2022; 50: 104118.
- Huang W, Kazemi-Kamyab H, Sun W, Karen Scrivener. Effect of cement substitution by limestone on the hydration and microstructural development of ultra-high performance concrete (UHPC). Cem Concr Compos. 2017; 77: 86-101.
- Fan Dingqiang, Rui Yu, Liu Kangning, Tan Junhui, Shui Zhonghe, Wu Chunfeg, et al. Optimized design of steel fibres reinforced ultra- high performance concrete (UHPC) composites: Towards to dense structure and efficient fibre application. Constru and Build Mat, 2021; 273: 121698.
- Shamsad A, Amin Al-F, Ashraf AB, Moruf OY. Fracture toughness of UHPC mixtures: Effects of w/b ratio, cement and silica fume contents. Constr Buil Mat. 2024; 417: 135327.
- Prem PR, Ramachandra MA, Bharatkumar BH. Influence of curing regime and steel fibres on the mechanical properties of UHPC. Mag Concr Res. 2015; 67: 988-1002.
- Jiang Du, Weina Meng, KamalH. Khayat, Yi Bao, Pengwei Guo, et al. New development of ultra-high-performance concrete (UHPC). Comp Part B Eng. 2021; 224: 109220.
- Yuhong Zhang, Jiaqi Xu, Xianbo Dong, Jiabao Wang, Changfeng Liu, Jiaju Liu. Optimization of nitrogen removal conditions based on response surface methodology and nitrogen removal pathway of Paracoccus sp. QD-19. Sci Total Environ. 2024; 908: 168348.
- Xia Zhang, Panjie Su, Weichao Wang, Wencai Yang, Yuanyuan Ge, Kuaile Jian, et al. Optimized carbonization of coffee shell via response surface methodology: A circular economy approach for environmental remediation. Environmental Pollution. 2023; 123018.
- Marcos AB, Ricardo ES, Eliane PO, Leonardo SV, Luciane AE. Response surface methodology (RSM) as a tool for optimization in analytical chemistry. Talanta. 2008; 76: 965-977.
- CM Ho, SI Doh, X Li, SC Chin, T Ashraf. RSM-based modelling of cement mortar with various water to cement ratio and steel slag content. Phy Che of Earth A/B/C. 2022; 128: 103256.
- Li Z, Lu D, Gao X. Multi-objective optimization of gap-graded cement paste blended with supplementary cementitious materials using response surface methodology. Const Bui Mate. 2020; 248: 118552.
- Ferdosian I, Camões A. Eco-efficient ultra-high performance concrete development by means of response surface methodology. Cem Concr Compos. 2017; 84: 146-156.
- Wille K, Naaman AE, Parra-Montesinos GJ. Ultra-High Performance Concrete with Compressive Strength Exceeding 150 MPa (22 ksi): A Simpler Way. ACI Materials Journal. 2011; 108.
- Barham HM, Aryan Far HS, Rabar HF, Hakar HQ, Khaleel HY. Mechanical properties and ductility behavior of ultra-high performance fiber reinforced concretes: Effect of low water-to-binder ratios and micro glass fibers. Ain Shams Eng J. 2021; 12: 1557-1567.
- Prem PR, Ramachandra MA, Bharatkumar BH. Influence of curing regime and steel fibres on the mechanical properties of UHPC. Magaz of Conc Res. 2015; 67: 988-1002.
- Doo-Yeol Y, Soonho K, Gi-Joon P, Jung-Jun P, Sung-Wook K. Effects of fiber shape, aspect ratio, and volume fraction on flexural behavior of ultra-high-performance fiber-reinforced cement composites. Composite Structures. 2017, 174: 375-388.
- Kinda H, Hui B, William P-A, Balaji R. Effect of different types of fibers on the microstructure and the mechanical behavior of ultra-high performance fiber-reinforced concretes. Compos B Eng. 2016; 86: 214-220.
- Ru Mu, Jiao C, Xiangshang C, Chengran D, Xiaowei W, Longbang Q. Effect of the orientation of steel fiber on the strength of ultra-high- performance concrete (UHPC). Constr Buil Mat. 2023; 406: 133431.
- Wang R, Gao X, Huang H, Guangshuai H. Influence of rheological properties of cement mortar on steel fiber distribution in UHPC. Constr Buil Mat. 2017; 144: 65-73.
- Fan M, Li Y, Shen J, Kaikai Jin, Junjie Shi. Multi-objective optimization design of recycled aggregate concrete mixture proportions based on machine learning and NSGA-II algorithm. Adv Eng Softw. 2024; 192: 103631.
- Zhou D, Mei Y, Ke X, Ziyang Liu. Preparation and use of new green underwater coagulation grouting materials based on NSGA-II algorithm-A case study. Const Buil Mate. 2024; 136407.
- Tanhadoust A, Madhkhan M, Nehdi ML. Two-stage multi-objective optimization of reinforced concrete buildings based on non- dominated sorting genetic algorithm (NSGA-III). J Buil Eng. 2023, 75: 107022.
- Sankar B, Ramadoss P. Modeling the compressive strength of high- performance concrete containing metakaolin using distinctive statistical techniques. Res Con Opti. 2023; 100241.
- Li Q, Fan Y. Rheological evaluation of nano-metakaolin cement pastes based on the water film thickness. Constr Build Mater. 2022; 324: 126517.
- Amparano FE, Yunping Xi, Young-Sook R. Experimental study on the effect of aggregate content on fracture behavior of concrete. Eng Fract Mech. 2000: 67: 65-84.
- Rao G A. Generalization of Abrams’ law for cement mortars. Cem Concr Res. 2001; 31: 495-502.
- Mohamad G, Lourenço PB, Camões A, Humberto R. Mechanical characterization study of laying mortars for structural masonry. 2007.
- Cordon WA, Gillespie HA. Variables in concrete aggregates and Portland cement paste which influence the strength of concrete. Engi Mate Sci. 1963; 60: 1029-1052.
- Pourbaba M, Asefi E, Sadaghian H, Amir M. Effect of age on the compressive strength of ultra-high-performance fiber-reinforced concrete. Const Bui Mat. 2018; 175: 402-410.
- Giaccio G, Zerbino R. Failure mechanism of concrete: combined effects of coarse aggregates and strength level. Adv Cem Based Mater. 1998; 7: 41-48.
- Fernandes V, Silva L, Ferreira VM, JA Labrincha. Evaluation of mixing and application process parameters of single-coat mortars. Ceme conc res. 2005; 35: 836-841.
- Wu Z, Khayat KH, Shi C. How do fiber shape and matrix composition affect fiber pullout behavior and flexural properties of UHPC?. Cem Conc Comp. 2018; 90: 193-201.
- Abed M, Fo?t J, Rashid K. Multicriterial life cycle assessment of eco- efficient self-compacting concrete modified by waste perlite powder and/or recycled concrete aggregate. Cons Buil Mate. 2022; 348: 128696.
- Doo-Yeol Y, Soonho K, Gi-Joon P, Jung-Jun P, Sung-Wook K. Effects of fiber shape, aspect ratio, and volume fraction on flexural behavior of ultra-high-performance fiber-reinforced cement composites. Comp Struc. 2017; 174: 375-388.ces