Next Article in Journal
Effect of Air Inlet Spacing on Air Distribution and Ventilation Performance of Stratum Ventilation Using Multiple Parallel Jets
Previous Article in Journal
Evaluating the Impact of 222 nm Far-UVC Radiation on the Aesthetic and Mechanical Properties of Materials Used in Public Bus Interiors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Design and Optimisation of Self-Compacting High Early-Strength Cement-Based Mortars

by
Nara Cangussu
1,2,
Ana Mafalda Matos
1,2,*,
Paula Milheiro-Oliveira
2,3 and
Lino Maia
1,2,4
1
CONSTRUCT-LABEST, Faculty of Engineering, University of Porto, 4200-465 Porto, Portugal
2
Faculty of Engineering, University of Porto, 4200-465 Porto, Portugal
3
CMUP, Faculty of Sciences, University of Porto, 4200-465 Porto, Portugal
4
Faculty of Exact Sciences and Engineering, Campus da Penteada, University of Madeira, 9020-105 Funchal, Portugal
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(7), 4142; https://doi.org/10.3390/app13074142
Submission received: 17 February 2023 / Revised: 20 March 2023 / Accepted: 21 March 2023 / Published: 24 March 2023
(This article belongs to the Section Civil Engineering)

Abstract

:
The use of SCC in Europe began in the 1990s and was mainly promoted by the precast industry. Precast companies generally prefer high early-strength concrete mixtures to accelerate their production rate, reducing the demoulding time. From a materials science point of view, self-compacting and high early-strength concrete mixes may be challenging because they present contradicting mixture design requirements. For example, a low water/binder ratio (w/b) is key to achieving high early strength. However, it may impact the self-compacting ability, which is very sensitive to Vw/Vp. As such, the mixture design can be complex. The design of the experimental approach is a powerful tool for designing, predicting, and optimising advanced cement-based materials when several constituent materials are employed and multi-performance requirements are targeted. The current work aimed at fitting models to mathematically describe the flow ability, viscosity, and mechanical strength properties of high-performance self-compacting cement-based mortars based on a central composite design. The statistical fitted models revealed that Vs/Vm exhibited the strongest (negative) effect on the slump-flow diameter and T-funnel time. Vw/Vp showed the most significant effect on mechanical strength. Models were then used for mortar optimisation. The proposed optimal mixture represents the best compromise between self-compacting ability—a flow diameter of 250 mm and funnel time equal to 10 s—and compressive strength higher than 50 MPa at 24 h without any special curing treatment.

1. Introduction

The building and construction industries have a heavy environmental impact. They account for 40% of natural resources, consume 70% of electrical power and 12% of potable water, and produce 45–65% of landfilled waste and 48% of GHG emissions [1]. In Europe, construction is the most significant economic sector, employing 18 million people directly and creating close to 10% of the European Union Gross Domestic Product. Thus, construction is an essential pillar in the context of the 2030 Agenda since it impacts the economic, environmental, and social spheres.
Concrete is a massively used construction material and the second most-consumed material worldwide, just after water [2]. This is mainly because the basic constituents, cement, aggregates, and water, are widely available worldwide; thus, concrete can be produced locally and consequently, it is inexpensive. Moreover, it hardens quickly and in all habitable environments, including underwater.
Over recent decades, efforts have been made to improve the behaviour of concrete materials. It took more than 150 years from the patenting of Portland cement (in 1824) until the discovery and development of self-compactable concrete (SCC) in the 1980s in Japan [3,4]. The development of superplasticisers allowed the efficient dispersal of cementitious particles and the emergence of SCC. According to the European standard EN 206-1 (specification, performance, production, and conformity of concrete), the constituent raw materials of SCC are similar to ordinary vibrated concrete (OVC). The main differences between OVC and SCC are due to the constituents’ proportions. Figure 1 presents the indicative volumetric proportions of SCC and OVC [5]. As can be seen, the main difference between SCC and OVC is a higher dosage of fine materials, i.e., cement and supplementary cementitious materials (SCM), and a smaller amount of coarse aggregate and smaller particle size distribution [6].
As such, the design of more advanced cement-based materials such as SCC can be complex since several requirements must co-exist, such as (i) engineering properties, (ii) architecture, (iii) cost, and (iv) eco-efficiency, among others. Even though the selection of the material should try to improve robustness, making the mixture more tolerant to raw material variability, for economic reasons, depends significantly on local availability; hence, there is no fixed rule for the amount/type of aggregates, cement, additions, and admixtures [6]. Therefore, a decisive factor for designing advanced cementitious materials, such as self-compacting concrete (SCC), is a clear understanding of each constituent raw material’s effect and interactions on the final product properties [7]. Additionally, the design and optimisation of such mixtures require intensive laboratory testing, particularly if new and unconventional supplementary cementitious materials or aggregates are incorporated [8], such as agricultural or industrial by-products [9]. Thus, a more scientific and multi-scale approach to mix-design is needed in which key mixture design variables can explain the composite properties [10].
The Design of Experiments (DoE) is a systematic approach to understanding how process and product parameters affect response variables, such as processability, physical properties, or product performance. DoE uses statistical methodology to analyse data and predict product property performance under all possible conditions within limits selected for the experimental design [10]. From a literature survey, the DoE approach [11] has shown promise for aiding the understanding of the effect of mixture parameters on key fresh and hardened SCC properties. It can also facilitate the test protocol required to optimise SCC [12] or quantify SCC mixtures’ robustness [6] to develop a more eco-friendly SCC that incorporates different types of SCM and meets multiple performance requirements [3,13,14,15,16,17,18,19,20,21].

Research Significance and Objectives

Previous research shows that statistical models and numerical optimisation techniques allow for finding the best combination of self-compacting high-strength concrete (SCHPC) constituents which (i) reduce the risk of cracking during the first days, ensuring a more impermeable concrete in the final structure [3]; (ii) study the influence of mixture parameters and their coupled effects on deformability, viscosity, compressive strength, resistivity and resistance to carbonation [15]; (iii) optimise chloride permeability and the strength of concrete containing metakaolin [22]; (iv) optimise the durability and service life of self-consolidating concrete containing metakaolin [19]; (v) optimise the mixture proportions and analyse the fresh and hardened properties [20]; and (vi) minimise the cost of the most unsustainable concrete components [21], including cost and sustainability indicators, in a reduced time [3,14,21,23]. These models and counterplot charts are beneficial for predicting and selecting the optimum mixture proportions for a given application simply and accurately.
From a materials science point of view, self-compacting and high early-strength concrete mixes may be challenging because they contradict mixture design requirements and fresh-state properties. For example, a low water/binder ratio (w/b) is key to achieving high early strength [24]. However, it may also impact self-compacting ability, which is very sensitive to Vw/Vp.
Maia [25] carried out a full central composite design on self-compacting high early-strength cement-based mortars (SCHSCM) based on key mixture design parameters and described the target engineering properties, namely, self-compacting ability and high early-age strength. The key five SCHSM mixture parameters chosen were [25]: (i) Vw/Vc—water to cement volume ratio; (ii) Sp/p—superplasticiser to powder mass ratio; (iii) Vw/Vp—water to powder volume ratio; (iv) Vs/Vm—sand to mortar volume ratio; and (v) Vfs/Vs—fine sand to total sand volume ratio to produce a total of 50 SCHSM mixtures. The current work aims at statistical analysis, model fitting, validation, and optimisation of the data experiments carried out by Maia [25]. As such, combining the statistical design methods and regression analysis allowed for creating a model to mathematically describe the flow ability, viscosity, and flexural and compressive strength of SCHSCM. The models obtained were validated and used to optimise self-compacting mixtures with additional high early-age strength. The goal was to develop SCHSM with a minimum slump flow of 250 mm, flow time of 10 s, and compressive strength of at least 50 MPa at 24 h without any special curing treatment.
The current work starts with a brief explanation of the data obtained by Maia [25] in Section 2. The model’s fitting and adequacy are detailed in Section 3, and Section 4 presents the optimisation of SCHSCM, considering the performance requirements established by the authors. Finally, Section 5 presents the main conclusions.

2. Central Composite Design

In general, the DoE involves the following key steps: (i) a choice of factors (mixture parameters), levels, and ranges; (ii) a selection of response variable(s); (iii) a choice of experimental design; (iv) performing the experiments; (v) a statistical analysis of the data (fitting a model); and (vi) other computations with response models. The current work focused on steps (iv) and (v). A brief explanation of the design adopted is presented in Section 2.1, and the output properties are discussed in Section 2.2.

2.1. Experiment Definition and Planning

The present work sought to analyse 50 mortar mixes published by Maia [25] using DoE, the statistical methodology to analyse data and predict properties. It was necessary to integrate simple and conventional statistical methods into the experimental design methodology to draw statistically sound conclusions from the experiment. Subsequently, the commercial software Design-Expert was employed in the current study to develop a polynomial regression model.
Maia’s data set [25] corresponds to a full central 25 composite design augmented by 10 axial runs plus 8 central runs. It allows for the development of a polynomial regression model for each response variable as a function of the five design variables considered in that study (key factors): (i) Vw/Vc—water to cement volume ratio; (ii) Sp/p—superplasticiser to powder mass ratio; (iii) Vw/Vp—water to powder volume ratio; (iv) Vs/Vm—sand to mortar volume ratio; and (v) Vfs/Vs—fine sand to total sand volume ratio. In addition to the 50 mixes previously described, 14 extra mixes were produced to allow for comparing and validating the results of the response models. The effect of each key factor was evaluated at five levels −α, −1, 0, +1, and +α, as shown in Table 1. To make the design rotatable, the value of α was equal to nF1/4, where nF is the number of points in the factorial part of the design. In the current study, this corresponds to α equal to 2.378.
The self-compacting mortars were produced with the following commercial materials available (in Portugal): (i) cement CEM I 42.5 R with a specific gravity of 3100 g/m3, (ii) limestone filler with a specific gravity of 2680 kg/m3, (iii) siliceous natural sands with a maximum particle size of 4 mm, specific gravity of 2620 kg/m3, and water absorption of 0.40%, (iv) fine siliceous sand with a maximum particle size of 2 mm, specific gravity of 2630 kg/m3, and water absorption of 0.20%, and (v) a polycarboxylate based superplasticiser with a specific gravity of 1080 kg/m3 and solid content of 40% [25].
Four dependent variables (i.e., response variables) were considered: (i) slump-flow diameter (D-flow); (ii) the time in the V-funnel (T-funnel); (iii) flexural strength at 24 h, determined by the three-point loading method (F,24 h); and (iv) uniaxial compressive strength at 24 h (Rc,24 h). Table 2 summarises the response variables and is followed by a brief explanation of the measurement procedures to access the dependent variables.
Immediately after production, the (mini) slump flow was assessed following the procedures of EFNARC [26]. The procedure consists, in brief, in filling the mini cone with the fresh composite and then taking the cone off and observing the flow capacity on a plane surface. When the fresh composite stops moving, two perpendicular diameters are measured, and the average value is considered the flow diameter (D-flow) in mm.
Additionally, the V-funnel test was performed to assess the viscosity and passing ability using the mini funnel according to EFNARC [26]. In brief, the fresh composite was poured into the funnel after production. Then, the trap door was opened, allowing the fresh composite to flow out under gravity, and the time for the discharge to complete (the flow time), i.e., when light is seen from above through the funnel, was recorded. These results showed a typical trend of an increase in slump-flow diameter with a reduction in flow time. After fresh tests, three specimens 40 × 40 × 160 {mm3} were moulded. The specimens were demoulded at 24 h, and the mechanical strength was assessed following EN 196-1 [27].

2.2. Results

The experimental results obtained by Maia concerning the D-flow, T-funnel, F,24 h, and Rc,24 h response variables used for model fitting (Section 3.2) are presented in Table 3, as well as 14 validation points (Vi). The D-flow results presented in Table 3 correspond to the average of the two results obtained by Maia [25]. The F,24 h corresponds to an average of three results for each CCD mortar point, and the Rc,24 h is the average of six results, based on Maia [25].

3. Response Models

3.1. Test Results Analysis

Table 4 presents the statistical summaries of all the experimental results obtained (Table 3) and the eight central point’s results only. The perception of the coefficient of variation on the central points (Ci) and all CCD points (all Ci, Fi, and CCi) predicts a good fit for the D-flow, T-funnel, and Rc,24 h. On the other hand, good fits are not expected for F,24 h given that the coefficient of variation evaluated on the central points (4.16%) is close to the coefficient of variation assessed on the 50 points of the CCD (6.93%).

3.2. Model Fitting

Design-Expert software was used to assist in analysing the results for each response variable by examining summary plots of the data; fitting a model using regression analysis and ANOVA; validating the model by examining the residuals for trends and outliers, leverage points, autocorrelation, and violation of statistical assumptions, in general; and interpreting the model graphically. A significance level of 5% was used throughout. A detailed description of this procedure is presented in Section 3.2 to Section 3.5. The central composite design adopted allows for estimating a full quadratic model, as shown in Equation (1) [28,29].
Y = β 0 + i = 1 k β i X i + i = 1 k β i i X i 2 + i < j β i j X i X j + ε
where Y represents the response variable; Xi corresponds to the design variables considered; β is used for model parameters (β0 is the independent term, βi represents the linear effect of Xi, βii represents the quadratic effect of Xi, and βij represents the linear-by-linear interaction between Xi and Xj); and ε is the fitting error. The model parameters (β0, βi, and βij) can be estimated through multilinear regression analysis. In the course of the analysis, some of the terms in Equation (1) may not be significant. The results of the fitting procedure for the four response variables are summarised in Table 5, Table 6, Table 7, Table 8, Table 9, Table 10, Table 11 and Table 12. The model generated with all observed values for D-flow and T-funnel resulted in a significant lack of fit. This fact led to the withdrawal of the axial runs for the D-flow results and also the inclusion of one term (Sp/p) × (Vfs/Vs) in the model because the term was considered the most significant among those excluded, due to its low p-value (0.0574). In the case of the T-funnel model, in addition to applying the inverse transformation as a recommended Box–Cox transformation to the response, one observation was removed (F15 = 26.33, see Table 3). The F,24 h and Rc,24 h fitting considered all runs and resulted in a non-significant lack of fit.

3.3. D-Flow Model

For the D-flow model, the authors decided not to include the axial points (CC1, CC2, CC3, CC4, CC5, CC6, CC7, CC8, CC9, and CC10) because the D-flow values contained a wide range of spread flow diameters, as presented in Section 2.2, from 168 to 398 mm. The attempt to fit a quadratic model led to a significant lack of fit when axial runs were included, due essentially to non-linearities that a quadratic model is not able to explain. Moreover, axial points are associated with high Cook’s distances and leverage points. Further research should be put into finding a non-linear model for D-flow response.
Table 5 presents the ANOVA for D-flow data considering the aforementioned conditions. The reduced quadratic model is significant. Examining the magnitude of the effects, it can be perceived that Vs/Vm (factor X4) and Vw/Vp (factor X3) are dominant compared to Sp/p (factor X2), Vw/Vc (factor X1), and Vfs/Vs (factor X5). The analysis of variance confirms the previous interpretation. In this case, Vw/Vc, Sp/p, Vw/Vp, and Vs/Vm are significant terms of the model, i.e., the main effects Vw/Vc, Sp/p, Vw/Vp, and Vs/Vm, and the second order effects, (Vw/Vc) × (Vw/Vp), (Vw/Vp) × (Vfs/Vs), and (Vw/Vc)2, are those that most significantly influence the response Y1: D-flow. The predicted R2 of 0.9683 is in reasonable agreement with the adjusted R2 of 0.9802 (i.e., the difference is less than 0.2). The coefficients of the estimated regression model for D-flow (in coded values) are presented in the second column of Table 6. Since the second-order effect (Vw/Vp) × (Vfs/Vs) is considered a significant factor, Vfs/Vs was added to the model. Furthermore, (Sp/p) × (Vfs/Vs) was forced to enter the model to avoid a significant lack of fit. Table 6 presents the model’s coefficients’ standard error and confidence intervals.
Table 5. ANOVA for the D-flow response.
Table 5. ANOVA for the D-flow response.
Test forSourceSum of SquaresDegrees of FreedomMean of SquareF-Valuep-Value
Significance of RegressionModel59,008.65115364.42176.44<0.0001
Residual851.292830.4
Total59,859.9439
Lack of FitLack of Fit774.792136.893.380.0523
Pure Error76.5710.93
Partial significance of each predictor variableVw/Vc2547.212547.283.78<0.0001
Sp/p670.71670.722.06<0.0001
Vw/Vp12,070.7112,070.7397.02<0.0001
Vs/Vm39,025.2139,025.21283.59<0.0001
Vfs/Vs64.7164.72.130.1558
(Vw/Vc) × (Sp/p)267.381267.388.790.0061
(Vw/Vc) × (Vw/Vp)919.131919.1330.23<0.0001
(Vw/Vc) × (Vs/Vm)309.381309.3810.180.0035
(Sp/p) × (Vfs/Vs)114.381114.383.760.0626
(Vw/Vp) × (Vfs/Vs)285.011285.019.370.0048
(Vw/Vc)22734.8912734.8989.95<0.0001
The three most significant parameters are typed in bold, and the most significant term is also underlined.
Table 6. D-flow model coefficients in terms of coded values.
Table 6. D-flow model coefficients in terms of coded values.
FactorCoefficient EstimateStandard Error95% CI Low95% CI High
Intercept343.751.95339.76347.74
Vw/Vc8.920.97476.9310.92
Sp/p4.580.97472.586.57
Vw/Vp19.420.974717.4321.42
Vs/Vm−34.920.9747−36.92−32.93
Vfs/Vs−1.420.9747−3.420.5748
Vw/Vc × Sp/p−2.890.9747−4.89−0.894
Vw/Vc × Vw/Vp−5.360.9747−7.36−3.36
Vw/Vc × Vs/Vm3.110.97471.115.11
Sp/p × Vfs/Vs1.890.9747−0.1063.89
Vw/Vp × Vfs/Vs2.980.97470.98774.98
(Vw/Vc)2−20.672.18−25.14−16.21
The three most significant parameters are typed in bold, and the most significant term is also underlined.
Figure 2a presents the normal probability plot of the model’s residuals. It is reasonable to consider that there are no relevant normality issues in the data. The plot of externally studentised residuals is presented in Figure 2b. One can notice an observation falls outside the −2 to +2 band, corresponding to the spreading diameter measurement of F5 with a value of 375 mm. It was decided to keep this observation since it did not constitute a leverage point, and the Cook’s distance was well below 50%.

3.4. T-Funnel Model

Regarding the T-funnel data, attempts to fit a quadratic model to all points lead to a severe lack of fit, potentially due to strong non-linearities. F15 exhibits an abnormal residual which compromises the fit; the following modelling steps exclude this point (see Figure 3).
The inverse transformation was applied to the T-funnel response, as found in previous works [3,14]. Examining the magnitude of the effects, it can be perceived that Vs/Vm (factor X4) and Vw/Vp (factor X3) are dominant compared to Sp/p (factor X2), Vw/Vc (factor X1), and Vfs/Vs (factor X5), as obtained in the D-flow response (see Table 7). All five factors are significant, as well as three of the interactions and the other three s-order terms. The predicted R2 of 0.9827 is in reasonable agreement with the adjusted R2 of 0.9897, i.e., the difference is negligible.
Table 7. ANOVA results for the T-funnel response.
Table 7. ANOVA results for the T-funnel response.
Test forSourceSum of SquaresDegrees of FreedomMean of SquareF-Valuep-Value
Significance of RegressionModel0.0640110.0058394.23<0.0001
Residual0.0005340
Total0.064545
Lack of FitLack of Fit0.00052703.140.0607
Pure Error075.457 × 10−6
Partial significance of each predictor variableVw/Vc0.004110.0039305.22<0.0001
Sp/p0.001010.001182.2<0.0001
Vw/Vp0.025110.02431895.62<0.0001
Vs/Vm0.031710.02712113.06<0.0001
Vfs/Vs0.000110.00016.90.0131
(Vw/Vc) × (Sp/p)0.000110.000111.610.0018
(Vw/Vc) × (Vs/Vm)0.000410.000323.66<0.0001
(Vw/Vp) × (Vs/Vm)0.001110.000973.86<0.0001
(Vw/Vc)20.000110.00016.790.0138
(Vw/Vp)20.000110.00014.840.0351
(Vfs/Vs)20.000110.00018.390.0067
The three most significant parameters are typed bold and the most significant term is also underlined.
The estimated coefficients of the regression model for the T-funnel (in coded values) are presented in the second column of Table 8. Additionally, Table 8 presents standard errors and confidence intervals. To check the regression assumptions on the residuals, a normal probability plot was examined in Figure 4a and other plots such as the one in Figure 4b. No outliers were detected, although F15 has been considered a leverage point. In the author’s view, this point does not considerably compromise the fitted model’s quality.
Table 8. T-funnel model coefficients (coded values).
Table 8. T-funnel model coefficients (coded values).
FactorCoefficient EstimateStandard Error95% CI Low95% CI High
Intercept0.06880.00110.06670.071
Vw/Vc0.01150.00070.01010.0129
Sp/p0.0050.00060.00380.0062
Vw/Vp0.02480.00060.02360.026
Vs/Vm−0.03020.0007−0.0315−0.0289
Vfs/Vs−0.00160.0006−0.0028−0.0004
(Vw/Vc) × (Sp/p)−0.0020.0007−0.0035−0.0006
(Vw/Vc) × (Vs/Vm)−0.00370.0007−0.0051−0.0022
(Vw/Vp) × (Vs/Vm)−0.00610.0007−0.0076−0.0047
(Vw/Vc)2−0.00170.0007−0.0031−0.0003
(Vw/Vp)20.00110.00050.00010.0021
(Vfs/Vs)2−0.00140.0005−0.0024−0.0003
The three most significant parameters are typed bold and the most significant term is also underlined.

3.5. Flexural Strength: F,24 h Model

No leverage points were found for the F,24 h model and the Cook’s distance was below 60%. Table 9 presents the ANOVA results for the F,24 h model. Examining the magnitude of the effects, it can be perceived that Vw/Vc (factor X1), Vw/Vp (factor X3), and Vs/Vm (factor X4) are dominant compared to Sp/p (factor X2) and Vfs/Vs (factor X5). The analysis of variance confirms the previous interpretation. The quadratic model is significant. In this case, Vw/Vc, Vw/Vp, and Vs/Vm are significant terms of the model, i.e., the main effects Vw/Vc, Vw/Vp, and Vs/Vm, and the second order effects (Vw/Vp) × (Vs/Vm) and (Vs/Vm)2 are those that most significantly influence the F,24 h response. The predicted R2 of 0.5671 is in reasonable agreement with the adjusted R2 of 0.6216 (i.e., the difference is less than 10%). Although R2 did not reach a value close to 1, the model was considered useful for navigating the design space.
Table 9. ANOVA results for the F,24 h response.
Table 9. ANOVA results for the F,24 h response.
Test forSourceSum of SquaresDegrees of FreedomMean of SquareF-Valuep-Value
Significance of RegressionModel18.6553.7316.44<0.0008
Residual9.53420.2269
Total28.1847
Lack of FitLack of Fit7.88350.22520.95740.5812
Pure Error1.6570.2352
Partial significance of each predictor variableVw/Vc7.1117.1131.36<0.0001
Vw/Vp2.9512.9513.010.0008
Vs/Vm2.9512.9512.980.0008
(Vw/Vp) × (Vs/Vm) 1.0711.074.700.0358
(Vs/Vm)24.5714.5720.16<0.0001
The three most significant parameters are typed in bold, and the most significant term is also underlined.
The estimated values of the regression model coefficients in the coded values are presented in the second column of Table 10. Additionally, Table 10 presents standard errors and confidence intervals. Figure 5a presents the normal probability plot of the residuals of the observed values of flexural strength. The externally studentised residuals are presented in Figure 5b and show no strong evidence of violating the white noise assumptions of the remaining error in the regression model. It is reasonable to consider that there are no normality issues in the data.
Table 10. F,24 h model coefficients in terms of coded values.
Table 10. F,24 h model coefficients in terms of coded values.
FactorCoefficient EstimateStandard Error95% CI Low95% CI High
Intercept11.430.089311.2511.61
Vw/Vc−0.40530.0724−0.5513−0.2592
Vw/Vp−0.2610.0724−0.4071−0.115
Vs/Vm0.26080.07240.11470.4068
(Vw/Vp) × (Vs/Vm)0.18260.08420.01270.3525
(Vs/Vm)2−0.28340.0631−0.4108−0.156

3.6. Compressive Strength: Rc,24 h Model

Table 11 presents the ANOVA results for Rc,24 h. Examining the magnitude of the effects, it can be perceived that Vw/Vc (factor X1), Vw/Vp (factor X3), and Vfs/Vs (factor X5) are dominant compared to Sp/p (factor X2). The reduced quadratic model is significant (p < 0.0001) and there is no lack of fit. The analysis of variance confirms the previous interpretation: Vw/Vc, Vw/Vp, and Vfs/Vs are significant terms of the model, i.e., the main effects Vw/Vc, Vw/Vp, and Vfs/Vs, the second order effects (Vs/Vm) × (Vfs/Vs), (Vw/Vc)2, and (Vw/Vp)2, are those that most significantly influence the response Rc,24 h. In addition, the predicted R2 of 0.9326 is reasonably close to the adjusted R2.
Table 11. ANOVA results for the Rc,24 h response.
Table 11. ANOVA results for the Rc,24 h response.
Test forSourceSum of SquaresDegrees of FreedomMean of SquareF-Valuep-Value
Significance of RegressionModel1850.3210185.03119.74<0.0001
Residual60.26391.35
Cor Total1910.5849
Lack of FitLack of Fit51.87321.621.350.3599
Pure Error8.4071.20
Partial significance of each predictor variableVw/Vc1650.7511650.75879.84<0.0001
Sp/p5.3715.373.480.698
Vw/Vp61.35161.3539.70<0.0001
Vs/Vm3.0513.051.980.1677
Vfs/Vs76.06176.0649.22<0.0001
(Vw/Vc) × (Vw/Vp)6.0816.083.930.0545
(Vw/Vc) × (Vfs/Vs)7.0917.094.590.0385
(Vs/Vm) × (Vfs/Vs)7.7917.795.040.0305
(Vw/Vc)218.43118.4311.930.0013
(Vw/Vp)217.44117.4411.290.0018
The three most significant parameters are typed in bold, and the most significant term is also underlined.
The estimated values of the coefficients of the regression model that will be used to predict Rc,24 h in the coded values are presented in the second column of Table 12. Additionally, Table 12 shows the standard errors and confidence intervals. Figure 6a illustrates the normal probability plot of the residuals of the observed compressive strength values. It is reasonable to consider that there are no normality issues in the data. An example plot of the residuals against the predicted values is presented in Figure 6b. There is no evidence of assumption violation regarding the regression model fit.
The subsequent diagnostics include the identification of possible outliers in the regression or observations that could be leverage points. Points with a Cook’s distance higher than 1.0 or high leverage were not found (see Figure 7).
Table 12. Rc,24 h model coefficients in terms of coded values.
Table 12. Rc,24 h model coefficients in terms of coded values.
FactorCoefficient EstimateStandard Error95% CI Low95% CI High
Intercept60.650.273360.0961.2
Vw/Vc−6.170.1889−6.56−5.79
Sp/p−0.35220.1889−0.73420.0299
Vw/Vp−1.190.1889−1.57−0.8081
Vs/Vm0.26550.1889−0.11650.6476
Vfs/Vs−1.330.1889−1.71−0.9431
(Vw/Vc) × (Vw/Vp)0.43570.2197−0.00880.8802
(Vw/Vc) × (Vfs/Vs)0.47060.21970.02610.9151
(Vs/Vm) × (Vfs/Vs)−0.49340.2197−0.9379−0.049
(Vw/Vc)20.56390.16330.23360.8941
(Vw/Vp)20.54860.16330.21830.8788

3.7. Significant Individual and Interaction Effects

The final models (in terms of coded factors) are illustrated in Figure 8 regarding the relative impact of each factor by comparing the corresponding coefficient with the remaining (higher values of the estimated coefficient indicate the higher influence of the design variable in the response). The sign (positive or negative) of the estimated coefficient is marked on each bar’s label. We recall that a positive coefficient means that the response (or transformed response) variable will increase if the given mixture parameter increases and vice-versa.
The response surfaces corresponding to the two main effects for each response of the mixtures are illustrated in Figure 9. For the D-flow and the T-funnel variation, the most significant factor is the volume of sand concerning the mortar (Vs/Vm) and the content of water/fines factors (Vw/Vp). In addition, it was possible to confirm the great influence of the water/cement factor (Vw/Vc) on the flexural (F,24 h) and compressive strength (Rc,24 h) of the mortars. In the order of inversely proportional influence, there is the fine sand concerning medium sand (Vfs/Vs) and the content of water/fines factors (Vw/Vp). The other factors have a smaller influence on the responses, exhibiting smaller estimated coded coefficients.
The equations in terms of actual factors can be used to predict the response for the given levels of each factor (Table 13 and Figure 10). Uncertainty in the fitted models (confidence interval in Table 6, Table 8, Table 10 and Table 12) is provided for a 95% confidence level.

3.8. Models Validation

In order to validate the statistical models presented in Table 13, Maia [25] produced more than 14 SCHSCM mixtures and performed the experiments described in Section 2. The mixture proportion (in coded values) for these mixtures is presented in the last 14 lines of Table 3. The experimental results obtained and the estimated results given by the regression models are illustrated in Figure 11, including the x = y diagonals. The ratio between the predicted-to-measured values ranges from 0.84 to 1.20. The most accurate predictions are given by the Rc,24 h and D-flow models, where predicted-to-measured values range from 0.96 to 1.02 and 0.92 to 1.05, respectively. In addition, as expected, the flexural strength obtained results that exhibited higher dispersions when compared to the other results.

4. Optimisation

The empirical models obtained for modelling and predicting the relevant SCHSM properties were also used for mixture optimisation purposes. The optimisation module in Design-Expert studied in this paper searches for a combination of factor levels that simultaneously satisfy the goals established for each response and factor (see Table 14). The procedure uses the desirability function (D), and the solutions with D ≥ 0.995 were considered optimal. The optimisation was developed in two different scenarios, as follows in Section 4.1 and Section 4.2.

4.1. Self-Compacting Ability

The first optimisation scenario was to find mortar mixtures with self-compacting ability. The self-compacting ability was translated by targeting optimisation restrictions of a D-flow = 260 mm and a T-funnel = 10 s. This optimisation problem is defined in the second column of Table 14 and is entitled “self-compacting ability”. The importance level for all constraints, both mixture parameters and response variables, was kept the same. The optimisation aimed only to obtain solutions within the CCD region. Thus, the constraint conditions for mixture parameters were “In CCD range”.
In total, 68 mixture SCHSM solutions were found with D ≥ 0.995 in the region 1.069   Vw / Vc   1.108 ; 0.019   Sp / p 0.029 ; 0.446 Vs / Vm 0.525 ; 0.540 Vw / Vp 0.706 ; and 0.088 Vsf / Vs 0.731 .

4.2. SCC with High Early Strength

The second optimisation scenario was intended to discover self-compacting and high early strength mixtures simultaneously (see the third column of Table 14, entitled “self-compacting ability and high early strength”). For that, the optimisation conditions were: D-flow in the range from 259 to 261 mm, T-funnel in the range from 9.8 to 10.2 s, and compressive strength higher than 50 MPa at 24 h. The importance level for all constraints, both mixture parameters and response variables were kept the same. The optimisation aimed only to obtain solutions within the CCD region; thus, the constraint conditions for mixture parameters were “In CCD range”.
The best combination of SCHSM mixture parameters was: Vc/Vp = 1.098; sp/p = 0.0243; Vw/Vp = 0.607; Vs/Vm = 0.478; and Vfs/vs = 0.047. The corresponding estimated properties and their 95% confidence interval are presented in Table 15.

5. Conclusions

The current study investigates the effect of SCHSCM mixture design variables. It selects the best combination of powder materials (in a binary mixture of Portland cement + MTCK), water/binder ratio, superplasticiser dosage, and aggregates to reach optimum values of self-compactability and high early-age strength. As such, the DoE approach was followed to obtain regression models that could describe and predict the relevant engineering properties. In addition, regression models were used for optimisation purposes. The following main conclusions were drawn:
-
Regression models were found to be adequate to describe SCHSCM properties over the experimental region, namely, slump-flow diameter, T-funnel time, and flexural and compressive strength at 24 h;
-
The Vs/Vm factor exhibited the strongest (negative) effect on the slump-flow diameter and T-funnel time;
-
The Vw/Vp factor was found to have the most significant effect on mechanical strength, corresponding to response variables F,24 h and Rc,24 h;
-
Vw/Vc was the second most influencing factor, with a positive effect on the slump-flow diameter and T-funnel time and a negative effect on F,24 h and Rc,24 h;
-
The variation of Sp/p used in the current CCD was small compared to the remaining factors. As such, Sp/p exhibited the lowest influence on SCHSCM properties when compared to other mixture parameters;
-
The proposed optimal mixture represented the best compromise between self-compacting ability—a flow diameter of 250 mm and funnel time equal to 10 s—and a compressive strength higher than 50 MPa at 24 h without any special curing treatment and was found for Vc/Vp = 1.098; sp/p = 0.0243; Vw/Vp = 0.607; Vs/Vm = 0.478; and Vfs/vs = 0.047.
Even though transposing the mortar study into SCC is out of the scope of the manuscript, SCC mixtures can be obtained by fixing the mortar mix proportions given in Section 4.2 and adding a coarse aggregate. However, tests on the concrete level are needed to optimise the aggregate skeleton. Moreover, based on regression models, contour plots or interaction diagrams can be computed and simplify the test protocol required to optimise a given SCC mixture, namely, to select the combination of powder materials with admixtures. In particular, they can help to compare the efficiency of different admixtures and alternative SCM. Optimised mortar mixtures can serve as reference mixtures in a quality control plan to detect variations in different deliveries of constituent materials, such as cement, filler, or superplasticiser. The derived models can also be useful to predict mortar properties when the material properties change, such as the cement source.
The procedure presented in this work can easily be implemented in any concrete laboratory of a production centre since it involves mortar tests which are easy to carry out with simple tests and equipment.

Author Contributions

Conceptualisation, N.C., A.M.M., P.M.-O. and L.M.; methodology, N.C, A.M.M., P.M.-O. and L.M.; software, N.C., A.M.M., P.M.-O. and L.M.; validation, A.M.M., P.M.-O. and L.M.; formal analysis, N.C., A.M.M., P.M.-O. and L.M.; investigation, N.C. and L.M.; resources, L.M.; data curation, A.M.M. and L.M.; writing—original draft preparation, N.C., A.M.M. and P.M.-O.; writing—review and editing, N.C., A.M.M., P.M.-O. and L.M.; visualisation, N.C., A.M.M., P.M.-O. and L.M.; supervision, A.M.M. and L.M.; project administration, A.M.M. and L.M.; funding acquisition, A.M.M. and L.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work is funded by: Base Funding—UIDB/04708/2020 of the CONSTRUCT—Instituto de I&D em Estruturas e Construções—funded by national funds through the FCT/MCTES (PIDDAC); by the national funds through FCT, Fundação para a Ciência e a Tecnologia, I.P., under the Scientific Employment Stimulus (Institutional Call, CEECINST/00049/2018); and the national funds through FCT, Fundação para a Ciência e a Tecnologia, I.P., under individual Scientific Employment Stimulus 2021.01765.CEECIND. The author Paula Milheiro-Oliveira was partially supported by CMUP, a member of LASI, which is financed by national funds through FCT—Fundação para a Ciência e a Tecnologia, I.P., under the project with reference UIDB/00144/2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Acknowledgments

Collaboration and materials supply by Secil, Omya Comital and Sika is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ANOVAAnalysis Of Variance
CCDCentral Composite Design
CiCentral points of CCD
CCiAxial points of CCD
dDays
DDesirability function
D-flowFlow Diameter (mm)
DoEDesign of Experiments
F,24 hFlexure strength at 24 h (MPa)
FiFactorial points of CCD
GDPGross Domestic Product
GHGGreenhouse Gas
hHours
LFLimestone Filler
pPowder Mass
PCPortland Cement
PSDParticle Size Distribution
Rc,24 hCompressive strength at 24 h (MPa)
SCCSelf-Compacting Concrete
SCHSCMSelf-Compacting High early Strength Cement-based Mortars
SCMSupplementary Cementitious Materials
SDGSustainable Development Goals
SpSuperplasticizer
Sp/pSuperplasticiser to powder mass ratio
RHRelative Humidity (%)
T-funnelFunnel Time
VcCement Volume
VCVibrated Concrete
ViValidation points
VfsFine Sand Volume
VIFValue of the Inflammation Factor
VmMortar Volume
VpPowder Volume
VsSand Volume
VwWater Volume
Vw/VcWater to cement volume ratio
Vw/VpWater to powder volume ratio
Vs/VmSand to mortar volume ratio
Vfs/VsFine sand to total sand volume ratio

References

  1. Omer, M.A.B.; Noguchi, T. A conceptual framework for understanding the contribution of building materials in the achievement of Sustainable Development Goals (SDGs). Sustain. Cities Soc. 2020, 52, 101869. [Google Scholar] [CrossRef]
  2. Antiohos, S.K.; Tapali, J.G.; Zervaki, M.; Sousa-Coutinho, J.; Tsimas, S.; Papadakis, V.G. Low embodied energy cement containing untreated RHA: A strength development and durability study. Constr. Build. Mater. 2013, 49, 455–463. [Google Scholar] [CrossRef]
  3. Nunes, S.; Matos, A.M.; Duarte, T.; Figueiras, H.; Sousa-Coutinho, J. Mixture design of self-compacting glass mortar. Cem. Concr. Compos. 2013, 43, 1–11. [Google Scholar] [CrossRef]
  4. Matos, A.M. Design of Eco-Efficient Ultra-High Performance Fibre Reinforced Cement-Based Composite for Rehabilitation/Strengthening Applications; University of Porto: Porto, Portugal, 2020. [Google Scholar]
  5. Okamura, H.; Ouchi, M. Self-Compacting Concrete. J. Adv. Concr. Technol. 2003, 1, 5–15. [Google Scholar] [CrossRef]
  6. Nunes, S.; Figueiras, H.; Milheiro Oliveira, P.; Coutinho, J.S.; Figueiras, J. A methodology to assess robustness of SCC mixtures. Cem. Concr. Res. 2006, 36, 2115–2122. [Google Scholar] [CrossRef]
  7. Yogitha, B.; Karthikeyan, M.; Muni Reddy, M.G. Progress of sugarcane bagasse ash applications in production of Eco-Friendly concrete—Review. Mater. Today Proc. 2020, 33, 695–699. [Google Scholar] [CrossRef]
  8. Gettu, R.; Pillai, R.G.; Santhanam, M.; Dhanya, B.S. Ways of Improving the Sustainability of Concrete Technology Through the Effective Use of Admixtures. Indian Inst. Technol. Madras Chennai 2018, 21. [Google Scholar]
  9. Loganayagan, S.; Chandra Mohan, N.; Dhivyabharathi, S. Sugarcane bagasse ash as alternate supplementary cementitious material in concrete. Mater. Today Proc. 2021, 45, 1004–1007. [Google Scholar] [CrossRef]
  10. Wagner, J.R.; Mount, E.M.; Giles, H.F. Design of Experiments. In Extrusion; William Andrew Publishing: Norwich, NY, USA, 2014; pp. 291–308. ISBN 978-1-4377-3481-2. [Google Scholar]
  11. Shi, C.; Wu, Z.; Lv, K.; Wu, L. A review on mixture design methods for self-compacting concrete. Constr. Build. Mater. 2015, 84, 387–398. [Google Scholar] [CrossRef]
  12. Khayat, K.H.; Ghezal, A.; Hadriche, M.S. Utility of statistical models in proportioning self-consolidating concrete. Mater. Struct. 2000, 33, 338–344. [Google Scholar] [CrossRef]
  13. Nunes, S.; Costa, C. Self-compacting concrete also standing for sustainable circular concrete. In Waste and Byproducts in Cement-Based Materials; Woodhead Publishing: Sawston, UK, 2021; pp. 439–480. [Google Scholar]
  14. Matos, A.M.; Maia, L.; Nunes, S.; Milheiro-Oliveira, P. Design of self-compacting high-performance concrete: Study of mortar phase. Constr. Build. Mater. 2018, 167, 617–630. [Google Scholar] [CrossRef]
  15. Dinakar, P. Design of self-compacting concrete with fly ash. Mag. Concr. Res. 2012, 64, 401–409. [Google Scholar] [CrossRef]
  16. Rojo-López, G.; Nunes, S.; González-Fonteboa, B.; Martínez-Abella, F. Quaternary blends of portland cement, metakaolin, biomass ash and granite powder for production of self-compacting concrete. J. Clean. Prod. 2020, 266, 121666. [Google Scholar] [CrossRef]
  17. Ghosh, D.; Abd-Elssamd, A.; Ma, Z.J.; Hun, D. Development of high-early-strength fiber-reinforced self-compacting concrete. Constr. Build. Mater. 2021, 266, 121051. [Google Scholar] [CrossRef]
  18. Al-Alaily, H.S.; Hassan, A.A.A. Refined statistical modeling for chloride permeability and strength of concrete containing metakaolin. Constr. Build. Mater. 2016, 114, 564–579. [Google Scholar] [CrossRef]
  19. Abouhussien, A.A.; Hassan, A.A.A. Optimising the durability and service life of self-consolidating concrete containing metakaolin using statistical analysis. Constr. Build. Mater. 2015, 76, 297–306. [Google Scholar] [CrossRef]
  20. Abouhussien, A.A.; Hassan, A.A.A. Application of Statistical Analysis for Mixture Design of High-Strength Self-Consolidating Concrete Containing Metakaolin. J. Mater. Civ. Eng. 2014, 26, 04014016. [Google Scholar] [CrossRef]
  21. Dvorkin, L.; Bezusyak, A.; Lushnikova, N.; Ribakov, Y. Using mathematical modeling for design of self compacting high strength concrete with metakaolin admixture. Constr. Build. Mater. 2012, 37, 851–864. [Google Scholar] [CrossRef]
  22. Matos, A.M.; Nunes, S.; Costa, C.; Barroso-Aguiar, J.L. Spent equilibrium catalyst as internal curing agent in UHPFRC. Cem. Concr. Compos. 2019, 104, 103362. [Google Scholar] [CrossRef]
  23. Nunes, S. Performance-Based Design of Self-Compacting Concrete (SCC): A Contribution to Enchance SCC Mixtures Robustness; Faculty of Enginnering of University of Porto: Porto, Portugal, 2008. [Google Scholar]
  24. Galobardes, I.; Salvador, R.P.; Cavalaro, S.H.P.; Figueiredo, A.; Goodier, C.I. Adaptation of the standard EN 196-1 for mortar with accelerator. Constr. Build. Mater. 2016, 127, 125–136. [Google Scholar] [CrossRef] [Green Version]
  25. Maia, L. Experimental dataset from a central composite design to develop mortars with self-compacting properties and high early age strength. Data Brief 2021, 39, 107563. [Google Scholar] [CrossRef] [PubMed]
  26. EFNARC. Specification and Guidelines for Self-Compacting Concrete; EFNARC: Farnham, UK, 2002. [Google Scholar]
  27. NP EN 196-1—Methods of Testing Cement. Part 1: Determination of Strength; Instituto Português da Qualidade: Lisbon, Portugal, 2006; pp. 1–37. [Google Scholar]
  28. De Barros Neto, B.; Scarminio, I.S.; Bruns, R.E. Como Fazer Experimentos, 4th ed.; Bookman, Ed.; Artmed Editora SA: Porto Alegre, Brazil, 2010; ISBN 978-85-7780-652-2. [Google Scholar]
  29. Montgomery, D.C. Introduction to Statistical Quality Control, 7th ed.; John Wiley and Sons Inc.: Hoboken, NJ, USA, 2013; ISBN 978-1-118-14681-1. [Google Scholar]
Figure 1. Mixture proportions of OVC and SCC (in % volume).
Figure 1. Mixture proportions of OVC and SCC (in % volume).
Applsci 13 04142 g001
Figure 2. (a) Normal probability plot of the externally studentised residuals for the D-flow fitted model, and (b) plot of the externally studentised residuals versus D-flow predicted values.
Figure 2. (a) Normal probability plot of the externally studentised residuals for the D-flow fitted model, and (b) plot of the externally studentised residuals versus D-flow predicted values.
Applsci 13 04142 g002
Figure 3. T-funnel results: F15 excluded (arrow).
Figure 3. T-funnel results: F15 excluded (arrow).
Applsci 13 04142 g003
Figure 4. (a) Normal probability plot of the externally studentised residuals for the T-funnel fitted model, and (b) plot of the externally studentised residuals versus T-funnel fitted values.
Figure 4. (a) Normal probability plot of the externally studentised residuals for the T-funnel fitted model, and (b) plot of the externally studentised residuals versus T-funnel fitted values.
Applsci 13 04142 g004
Figure 5. (a) Normal probability plot of the externally studentised residuals for the F,24 h fitted model, and (b) plot of the externally studentised residuals versus F,24 h fitted values.
Figure 5. (a) Normal probability plot of the externally studentised residuals for the F,24 h fitted model, and (b) plot of the externally studentised residuals versus F,24 h fitted values.
Applsci 13 04142 g005
Figure 6. (a) Normal probability plot of the externally studentised residuals for the Rc,24 h fitted model, and (b) plot of the externally studentised residuals versus Rc,24 h fitted values.
Figure 6. (a) Normal probability plot of the externally studentised residuals for the Rc,24 h fitted model, and (b) plot of the externally studentised residuals versus Rc,24 h fitted values.
Applsci 13 04142 g006
Figure 7. Rc,24 h: (a) Cook’s distance and (b) leverage.
Figure 7. Rc,24 h: (a) Cook’s distance and (b) leverage.
Applsci 13 04142 g007
Figure 8. Influence of each factor on fresh state responses, in which (+) Positive effect in response models; (−) Negative effect in response model (see Table 13).
Figure 8. Influence of each factor on fresh state responses, in which (+) Positive effect in response models; (−) Negative effect in response model (see Table 13).
Applsci 13 04142 g008
Figure 9. Effect of the Vs/Vm and Vw/Vp ratios on (a) D-flow and (b) T-funnel; effect of Vw/Vc and Vw/Vp on (c) F,24 h and (d) Rc,24 h.
Figure 9. Effect of the Vs/Vm and Vw/Vp ratios on (a) D-flow and (b) T-funnel; effect of Vw/Vc and Vw/Vp on (c) F,24 h and (d) Rc,24 h.
Applsci 13 04142 g009aApplsci 13 04142 g009b
Figure 10. Predicted versus actual plan: (a) D-flow, (b) T-funnel, (c) F,24 h, and (d) Rc,24 h.
Figure 10. Predicted versus actual plan: (a) D-flow, (b) T-funnel, (c) F,24 h, and (d) Rc,24 h.
Applsci 13 04142 g010aApplsci 13 04142 g010b
Figure 11. Predicted versus measured: (a) D-flow; (b) T-funnel; (c) F,24 h, and (d) Rc,24 h.
Figure 11. Predicted versus measured: (a) D-flow; (b) T-funnel; (c) F,24 h, and (d) Rc,24 h.
Applsci 13 04142 g011aApplsci 13 04142 g011b
Table 1. Correspondence between the coded values and actual values of the design variables.
Table 1. Correspondence between the coded values and actual values of the design variables.
Design Variables−2.378−10+1+2.378
X1: Vw/Vc0.6820.8050.8950.8941.108
X2: Sp/p0.0190.0220.0240.0250.029
X3: Vw/Vp0.4340.5130.5700.6270.706
X4: Vs/Vm0.3660.4320.4800.5280.594
X5: Vfs/Vs0.0430.2500.4000.5500.757
Table 2. Response variable summary.
Table 2. Response variable summary.
Response VariablesUnitsMeasurement Method
Y1: D-flowmmEFNARC
Y2: T-funnelsEFNARC
Y3: F,24 hMPaEN 196-1
Y4: Rc,24 hMPaEN-196-1
Table 3. The mixture proportion (coded values) and experimental results obtained by Maia [25] and used for model fitting and validation points.
Table 3. The mixture proportion (coded values) and experimental results obtained by Maia [25] and used for model fitting and validation points.
CCD PointCoded ValuesResults
Vw/VcSp/pVw/VpVs/VmVfs/VsD-FlowT-FunnelF,24 hRc,24 h
(mm)(s)(MPa)(MPa)
C10.000.000.000.000.00346.5014.9811.7759.89
C20.000.000.000.000.00346.0013.7712.2962.10
C30.000.000.000.000.00339.5014.1511.2859.31
C40.000.000.000.000.00348.0014.7412.0860.90
C50.000.000.000.000.00339.5014.5311.9660.00
C60.000.000.000.000.00341.0014.3811.8962.17
C70.000.000.000.000.00345.0013.5611.0560.18
C80.000.000.000.000.00344.5014.2311.0359.63
F1−1.00−1.00−1.00−1.00−1.00330.0018.8911.6072.54
F21.00−1.00−1.00−1.00−1.00349.0012.2611.0557.80
F3−1.001.00−1.00−1.00−1.00338.5016.7012.6470.09
F41.001.00−1.00−1.00−1.00354.5011.3210.5056.06
F5−1.00−1.001.00−1.00−1.00375.009.2411.2667.05
F61.00−1.001.00−1.00−1.00369.507.089.4355.72
F7−1.001.001.00−1.00−1.00376.007.8010.6068.47
F81.001.001.00−1.00−1.00375.006.9110.1853.40
F9−1.00−1.00−1.001.00−1.00253.00114.0611.6873.13
F101.00−1.00−1.001.00−1.00289.5038.9111.6458.49
F11−1.001.00−1.001.00−1.00260.5074.3210.4871.81
F121.001.00−1.001.00−1.00289.5033.0310.9159.30
F13−1.00−1.001.001.00−1.00295.0025.4011.8367.80
F141.00−1.001.001.00−1.00313.5017.5311.4155.93
F15−1.001.001.001.00−1.00303.5026.3311.4168.43
F161.001.001.001.00−1.00320.0014.1510.9356.48
F17−1.00−1.00−1.00−1.001.00300.0028.4311.0666.54
F181.00−1.00−1.00−1.001.00351.5012.8910.6955.93
F19−1.001.00−1.00−1.001.00334.5017.6711.5769.59
F201.001.00−1.00−1.001.00349.5011.7310.7355.10
F21−1.00−1.001.00−1.001.00378.009.9011.2065.92
F221.00−1.001.00−1.001.00379.507.169.9354.52
F23−1.001.001.00−1.001.00385.508.5410.7164.39
F241.001.001.00−1.001.00382.006.9410.2653.64
F25−1.00−1.00−1.001.001.00233.00*12.1067.34
F261.00−1.00−1.001.001.00275.0052.6411.4255.30
F27−1.001.00−1.001.001.00265.5067.3412.7968.54
F281.001.00−1.001.001.00285.0033.3010.4856.02
F29−1.00−1.001.001.001.00289.5026.2011.9865.12
F301.00−1.001.001.001.00315.0016.5110.9853.39
F31−1.001.001.001.001.00309.0018.5711.8860.96
F321.001.001.001.001.00314.0017.0710.6653.41
CC1−2.380.000.000.000.00168.00*11.9178.42
CC22.380.000.000.000.00342.0011.5510.2548.92
CC30.00−2.380.000.000.00330.5017.6411.4562.44
CC40.002.380.000.000.00337.0012.7111.2258.89
CC50.000.00−2.380.000.00295.0079.6312.1164.13
CC60.000.002.380.000.00368.507.2710.1863.04
CC70.000.000.00−2.380.00398.007.499.4359.77
CC80.000.000.002.380.00169.50*10.3262.62
CC90.000.000.000.00−2.38338.0016.55Na62.21
CC100.000.000.000.002.38330.5016.14Na57.75
V10.633.050.210.52−216.17320.5011.7058.45
V20.631.50.210.52−216.58327.0011.0960.53
V30.631.50.210.52−2.6718.33310.5011.1761.54
V4−0.130.21−0.66−0.620.0012.96353.0012.0261.33
V51.08−0.42−1.13−0.620.0014.09337.5011.4756.70
V6−1.340.52−0.14−0.620.0014.15347.0010.3570.59
V7−0.820.58−0.39−0.620.0013.64338.509.6865.96
V80.57−0.15−0.92−0.620.0013.02335.0011.1659.11
V9−0.13−0.42−0.66−0.620.0014.76350.0013.2963.51
V10−0.13−0.42−0.66−0.620.0015.01344.5011.1465.16
V111.08−0.42−1.13−0.620.0014.32346.5011.7657.25
V12−1.340.52−0.14−0.620.0015.65344.5013.9772.63
V13−0.820.58−0.39−0.620.0015.26345.5013.1268.96
V140.57−0.15−0.92−0.620.0015.08347.50Na59.54
* impossible to measure; Na—non-available result.
Table 4. The main statistics for all CCD points and central points.
Table 4. The main statistics for all CCD points and central points.
D-Flow (mm)T-Funnel (s)F,24 h (MPa)Rc,24 h (MPa)
All 50 mixtures
Minimum168.006.919.4348.92
Maximum398.00114.0612.7978.42
Mean323.3122.3911.1761.61
Standard deviation48.9721.630.776.24
Coefficient of variation (%)15.1596.626.9310.14
Central points (Ci)
Minimum339.5013.5611.0359.31
Maximum348.0014.9812.2962.17
Mean343.7514.2911.6760.52
Standard deviation3.310.470.481.10
Coefficient of variation (%)0.963.314.161.81
Table 13. Final equations in terms of actual factors: estimated coefficients.
Table 13. Final equations in terms of actual factors: estimated coefficients.
Model TermsD-Flow (mm)1/(T-Funnel) (s)F,24 h (MPa)Rc,24 h (MPa)
Independent−2154.81−1.35825.400293.77
Vw/Vc5375.511.2092−4.529−257.75
Sp/p15,060.0513.4743-−186.47
Vw/Vp1141.351.1255−36.616−289.80
Vs/Vm−1375.321.411985.48732.94
Vfs/Vs−365.990.0374-−7.31
(Vw/Vc) × (Vw/Vp)−1050.80--85.43
(Vw/Vc) × (Sp/p)−17,105.24−12.0954--
(Vw/Vc) × (Vs/Vm)723.96−0.8550--
(Vw/Vp) × (Vs/Vm)-−2.239066.741-
(Vw/Vp) × (Vfs/Vs)349.05--35.06
(Sp/p) × (Vfs/Vs)6673.75---
(Vs/Vm) × (Vfs/Vs)---−68.53
(Vw/Vc)2−2581.92−0.2148-70.43
(Vw/Vp)2-0.3375-168.84
(Vs/Vm)2--−123.018-
(Vfs/Vs)2-−0.0602--
Table 14. Optimisation criteria.
Table 14. Optimisation criteria.
Self-Compacting AbilitySelf-Compacting Ability and High Early Strength
Constraints for optimisation
Mixture parameters:
Vw/VcIn CCD rangeIn CCD range
Sp/pIn CCD rangeIn CCD range
Vw/VpIn CCD rangeIn CCD range
Vs/VmIn CCD rangeIn CCD range
Vfs/VsIn CCD rangeIn CCD range
Response Variables:
D-flow (mm)Target 260In range {259, 261}
T-funnel (s)Target 10In range {9.8, 10.2}
F,24 h (MPa)NoneNone
Rc,24 h (MPa)None>50
Table 15. Estimated SCHSM properties of the optimised mixture.
Table 15. Estimated SCHSM properties of the optimised mixture.
ResponsePredicted Results95% CI Low95% CI High
D-flow259.15237.78280.51
T-funnel9.879.110.79
F,24 h10.319.9310.7
Rc,24 h50.0046.8653.14
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cangussu, N.; Matos, A.M.; Milheiro-Oliveira, P.; Maia, L. Numerical Design and Optimisation of Self-Compacting High Early-Strength Cement-Based Mortars. Appl. Sci. 2023, 13, 4142. https://doi.org/10.3390/app13074142

AMA Style

Cangussu N, Matos AM, Milheiro-Oliveira P, Maia L. Numerical Design and Optimisation of Self-Compacting High Early-Strength Cement-Based Mortars. Applied Sciences. 2023; 13(7):4142. https://doi.org/10.3390/app13074142

Chicago/Turabian Style

Cangussu, Nara, Ana Mafalda Matos, Paula Milheiro-Oliveira, and Lino Maia. 2023. "Numerical Design and Optimisation of Self-Compacting High Early-Strength Cement-Based Mortars" Applied Sciences 13, no. 7: 4142. https://doi.org/10.3390/app13074142

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop