Abstract

Compressive strength of alkali-activated slag (AAS) concrete is influenced by multi-factors in a nonlinear way. Both artificial neural network (ANN) and alternating conditional expectation (ACE) models of 3-day (3 d) and 28-day (28 d) compressive strength of AAS were established in this study by using the data reported in related literature, where alkali concentration of activator (Na2O%), modulus of activator (Ms), water/binder ratio (W/B), surface area of slag (SA), and basicity index of slag (Kb) were taken as input parameters. The models were employed later to predict 3 d and 28 d compressive strength of AAS concretes, respectively, and the results were validated by experimental work. The results show that both the ANN and the ACE models had adequate accuracy, no matter 3 d or 28 d compressive strength was considered. Compared to the 3 d compressive strength, due to data scattering that increased with the increase of data size, both the models did not yield a higher accuracy in the case of 28 d strength. However, also due to the increase in data size, both the models were more feasible to implement 28 d strength prediction as a result of sufficient learning and training during modeling. In addition, based on ACE analysis, the weight-influencing compressive strength of AAS decreased in a sequence of Na2O% > Ms > W/B > Kb > SA. If data size was sufficiently large, it was more suitable to establish an ANN model for compressive strength prediction of AAS concretes. Otherwise, ACE could be considered as an alternative to yield an acceptable result.

1. Introduction

Alkali-activated slag (AAS), in which alkali activator (such as sodium silicate solution, known as water glass, WG) blends with ground granulated blast-furnace slag, has received widespread attention due to its extremely low CO2 emission and superior performance [1]. Concrete manufactured by using AAS instead of Portland cement (PC) is potentially employed in structural construction. On the other hand, performance testing of concrete is usually time-consuming via experimental work. In this context, prediction based on mathematic models has shown its practical significance. As one of the most basic and the most important performances of concrete, compressive strength is largely dependent on material components and their characteristics. Different from PC, whose compressive strength mainly depends on water-binder ratio, compressive strength dependence of AAS comes from both alkali activator and slag, as well as water-binder ratio. Such multi-factor conditions will certainly bring much more challenges for compressive strength prediction of AAS.

Bagheri et al. [2] established an artificial neural network (ANN) by using the data obtained from the authors’ experimental work, followed by using genetic programming (GP) to conduct a prediction for compressive strength of boron-based alkali-activated fly ash-slag system. To develop the ANN model, percentages of fly ash and slag, as well as ratios of B, Si, and Na in the alkali activator were taken as input variables. Totally 114 series of data were collected and 70%, 15%, and 15% of the data were used for training, testing, and validation, respectively. A very low error value of less than 0.1 confirmed the high accuracy of the model. After a GP based on 20 series of experiments, R2 of 0.95 and root mean square error (RMSE) of 0.07 were obtained, respectively, implying a good feasibility of the model for such a prediction. Also, an ANN model was developed for compressive strength of alkali-activated fly ash-slag concretes in the study carried out by Nagajothi and Elavenil [3]. Twenty groups of data obtained from their own experiment were used for establishment of the model, where slag-fly ash ratio and river sand–manufactured sand ratio were used as input parameters. Finally, mean absolute error (MAE), RMSE, and mean absolute percentage error (MAPE) of 0.042, 0.094, and 0.001 were obtained, respectively, exhibiting that the model had high accuracy. Faridmehr et al. [4] also developed an ANN model for compressive strength of alkali-activated fly ash-slag self-compacting concretes, where contents of fly ash and slag, as well as curing age were taken as input parameters. Within the total 6 groups of data based on their own experiment, 70% and 30% of the data were used for training and testing, respectively. Reliability of the model was evaluated by comparing its error to the ones obtained from multiple linear regression (MLR) model and genetic algorithm combined with ANN (GA-ANN) model. Finally, relatively lower MAE of 1.78, mean squared error (MSE) of 5.98, RMSE of 2.45, and average absolute error (AAE) of 0.04 were obtained, respectively. Shariatmadari et al. [5] used alkali-activated volcanic ash-slag-PC to stabilize sandy soil, and developed ANN and evolutionary polynomial regression (EPR) models based on their own experimental results. Replacements of volcanic ash and slag in PC, concentration of alkali activator, alkali to binder ratio, Na/Al, Si/Al, curing time, and curing temperature were employed as input parameters. Training, testing, and validation used 75%, 10%, and 15% of data, respectively. After several trials, based on RMSE and MAE of 0.0439 and 0.0336, respectively, the ANN model with architecture of 8-5-10-1 was proved as the most accurate one. Zhang et al. [6] proposed a chemical-informed machine learning model for compressive strength of alkali-activated fly ash-slag system. To build the model, relatively comprehensive features were considered as input parameters, including contents of fly ash, slag, NaOH, WG, water, fine aggregates and coarse aggregates, reactivity modulus, hydraulic modulus, silica modulus, alumina modulus, lime modulus, Na2O in WG, SiO2 in WG, water in WG, relative humidity, as well as age. Training and testing used 70% and 30% of data (676 groups in total from literature), respectively. Eventually, high accuracy was achieved with a low MAE of 3.228 MPa.

From the above, it can be seen that ANN has been successfully used to estimate compressive strength of alkali-activated materials. Most studies have been carried out on alkali-activated fly ash-slag and other alkali-activated systems rather than AAS, which, however, has more important significance in engineering practice currently. Furthermore, the data used in most studies were from their own experimental work. This could have resulted in a relatively small data size and then followed by a decrease in accuracy and feasibility when using the local optimum model to solve a global issue.

In this context, this study will focus on compressive strength of AAS concretes by using the ANN approach. In the development of the ANN model, data will be collected as fully as possible from literature. After the model with sufficient accuracy has been established, its feasibility will be evaluated by using own experimental work. For comparison, alternating conditional expectation (ACE) analysis will also be conducted and the weight of the influencing factors on the compressive strength of AAS concretes will be determined elementarily by using such method.

2. Data Collection and Modeling

2.1. Data Collection

ANN and ACE models were established in this study to predict compressive strength of AAS concretes with WG as alkali activator. Data used for the model establishment was collected from related literature published publicly. The influence on strength of AAS should be from both alkali activator and slag, as well as their proportion. Considering the importance of the influencing parameters, alkali concentration of WG (Na2O%), modulus of WG (Ms), water/binder ratio (W/B), surface area of slag (SA), and basicity index of slag (Kb) were taken as the input parameters of the models [713]. Other parameters, such as binder content, which are less important and then not typically considered in the literature, are not adopted in this study for model establishment. Na2O%, Ms, W/B, and Kb should be calculated according to formulae (1)–(4) [1417]. If NaOH was used in the alkali activator as well together with WG, NaOH should be treated as a part of WG. Therefore, “WG” in the formulae means the WG after NaOH adjustment. Compressive strengths of AAS at curing ages of 3 d and 28 d were considered as output parameters. Finally, 171 and 309 sets of data were successfully selected for 3 d and 28 d compressive strength prediction, respectively. Literature used and input/output parameter values listed in Tables 1 and 2 gives the number of tests, range of variation, and average and standard deviation values for each of the references presented in Table 1. Parts of the data in the Table 1 have same inputs but different output values. This is because other influencing factors rather than the five mentioned above were also considered in the literature, but here only the five factors were concerned.

2.2. ANN Modeling

ANN is a kind of general mathematical model to solve nonlinear system problems. It processes information by imitating human neural network [67]. Due to its strong generalization ability, it has been applied successfully in the field of civil engineering in recent years [6872]. ANN models consist of three layers: input layer, hidden layer(s), and output layer. The number of neurons in both the input and the output layers is just the number of parameters in the corresponding layers. Determination of the number of layers and number of neurons in each layer is relatively complex for the hidden layer(s), which is generally done by trial-and-error. Connection between neurons should be implemented by using a linear/nonlinear transfer function. Input parameters in the input layer are fed to the output layer through the hidden layer(s). If outputs were not desired, they would bring error signals to back-propagate into the input layer to minimize the error by modifying the weight of each neuron there. In this context, training is very important in ANN modeling, and sufficient data size is needed to guarantee its accuracy.

In this study, the five parameters, Na2O%, Ms, W/B, SA, and Kb, were treated as the neutrons in the input layer, and compressive strength of AAS was the neutron in the output layer. The determination of the hidden layer was made by trial and error. Both single-layer and double-layer were tried in this study. For both cases, number of the neutrons in the hidden layer was equal to root of the sum of input and output node numbers, plus a constant value less than 10 in general [73]. After several trials (a sample of performance comparison is given in Table 3), one hidden layer was determined to be used, and the number of the neutrons in the hidden layer was 10 and 13 for 3 d and 28 d compressive strength, respectively, according to the MAE and the RMSE error evaluations (to be described in section 2.4). Therefore, ANN model structures of 5-10-1 and 5-13-1 were used eventually in this study for compressive strength prediction of AAS at curing ages of 3 d and 28 d, respectively, as shown in Figure 1. For both the networks, 70%, 15%, and 15% of the data collected from literature were randomly selected for training, testing, and validation, respectively. Number of training period, learning rate, and minimum error of the training were set as 1000, 0.01, and 0.00001, respectively. After a comparison between Gradient Descent method and Levenberg-Marquardt algorithm, the latter was used for training due to its higher accuracy. The results to be reported were the optimum solution after 50 times of training.

2.3. ACE Modeling

ACE nonparametric regression technique, which has strong feasibility and high accuracy, is essentially a kind of response surface function. It is usually used to establish a mapping relationship between inputs and outputs which are already known by the approach of data fitting, as given in Ref. [74].where , are input parameters; are response surface functions of , ; is output parameter; is response surface function of ; is inverse function of ; and are accuracy errors.

As the accuracy error is difficult to be determined, correction coefficient is usually introduced to modify formula (6), and finally, formula (7) is obtained as given below.

Furthermore, influencing weight of each input parameter on the output parameter could be analyzed elementarily by using ACE approach as well during the calculation.

In this study, , were nominated as Na2O%, Ms, W/B, SA, and Kb, respectively, as input parameters. Their values were collected from literature as detailed in Section 2.1. Input parameters were imported in ACE calculation software S-PLUS to calculate their corresponding response surface function , . After that, a mapping relationship was established between the response surface functions and output parameters, namely, compressive strength of AAS at curing ages of 3 d and 28 d, respectively, by data fitting regression. The values of the strength were also collected from literature as detailed in Section 2.1.

2.4. Error Evaluation

Both mean absolute error (MAE) and root mean square error (RMSE) were used in this study to evaluate accuracy of both the ANN and the ACE models. MAE and RMSE values are to be zero when the values predicted and tested are exactly the same. Error increases with the increase of MAE and RMSE values. Calculations of MAE and RMSE are given in formulae (8) and (9) according to the literature [26]. Furthermore, agreement between predicted and measured values together with coefficient of determination was also applied to better assess the applicability of the proposed networks.where (, , ......) are the values tested; (, , ......) are the values predicted; n is the number of data.

2.5. Experimental Validation

Eleven AAS concrete mixes were prepared in this study, their compressive strength was tested at curing ages of 3 d and 28 d, respectively, to validate the results predicted by using both the ANN and the ACE models.

The slag used in the experimental work was from Qujing, Yunnan Province, China, whose Kb and SA were 0.96 and 405 m2/kg, respectively. Its chemical compositions are given in Table 4.

Industrial water glass (with Na2O% and SiO2% of 11.49% and 31.25%, respectively) supplied by Kunming, Yunnan Province was used as activator. NaOH (with purity above 96%) from Fuchen (Tianjin) Chemical Reagent Co., Ltd. was used to adjust the modulus of the WG to the values required. Crushed limestone with size of 5–25 mm (continuous grading) and manufactured sand with fineness modulus of 2.82 were used as coarse and fine aggregates, respectively. Tap water was used for mixing.

By using the mix proportion given in Table 5, concrete specimens were manufactured according to Chinese standard JGJ/T 439–2018 [75]. After mixing, the concrete mixture was cast into molds with size of 100 mm × 100 mm × 100 mm in two layers. Vibration was carried out after each layer was cast. When the concrete surface was stiff, the specimens with molds were moved into curing room with temperature of 20 ± 2°C and RH ≥ 95%. The specimens were de-molded one day after mixing, and then were moved into the curing room again for continuous curing to the ages required.

When the specimens reached curing ages of 3 d and 28 d, respectively, compressive strength test was carried out according to Chinese national standard GB/T 50081- 2019 [76]. The strength results to be reported are an average value of three duplicated specimens.

3. Results and Discussion

3.1. ANN

Regression of ANN modeling for prediction of 3 d and 28 d compressive strength of AAS is given in Figures 2 and 3, respectively. Correlation coefficient R is usually used to evaluate the correlation between the data samples. The closer the R-value is to 1, the stronger the correlation. From the figures, it can be seen that both the predictions have a good correlation between the data samples in general no matter training, testing, validation, or all of them are considered as the R-values are higher than 0.8. However, the correlation in the case of 28 d compressive strength is relatively lower in comparison with the 3 d one. It is known that 28 d compressive strength is the most elementary property of concrete, and therefore has been extensively reported in the literature. Consequently, large amounts of data samples on the 28 d compressive strength of AAS have been collected in this study, namely, 309 sets. This could have increased the scattering of the data as a result of different conditions of the specimen, such as different types and sizes of the specimen, different sources of slag, different aggregates, different curing conditions, used in the different studies.

Also, from MAE and RMSE error evaluation results, as given in Table 6, it can be seen that the error of the ANN models is at a low level, with MAE and RMSE values less than 20 MPa, in both the cases of 3 d and 28 d compressive strength, which means the models have sufficient accuracy.

Later, 3 d and 28 d compressive strengths of the 11 AAS concrete mixes were experimentally measured. Simultaneously, the material parameters used in the experimental work were put in the ANN model to run a prediction modeling. Comparison of the results measured and predicted is reported in Table 7, together with the error between them (100% × |predicted value-measured value|/measured value). From the table, it can be seen that in comparison to 28 d compressive strength, although the ANN model of 3 d compressive strength had an equivalent accuracy as discussed previously, it is less feasible to be employed to run the prediction. As stated previously, compared to the 309 sets of data that have been collected for the establishment of the 28 d compressive strength model, only 171 sets were available in the case of 3 d strength. Furthermore, the material parameters reported in the 171 data sets may not have covered the ones used in the experimental work, resulting in a lack of corresponding training, and finally in a high error between measurement and prediction. As can be seen from Table 1, few literatures consider the effects of 5 parameters on the 3D compressive strength of AAS at the same time. Therefore, as a result of insufficient training, the error of the mix 11 at 3 days is as high as 77.1%, while the larger database of 28 d compressive strength has supported effective training. Consequently, the error decreases dramatically to a low value of 4.6%. In this context, in order to improve feasibility of the ANN model for prediction, it is necessary to enlarge the scale of database enabling which could cover circumstances as full as possible, to ensure that a sufficient training could be implemented.

3.2. ACE

Relationship between input parameters , and their corresponding response surface function , is shown in Figures 4 and 5. ACE fitting curves of compressive strength of AAS at curing ages of 3 d and 28 d were plotted in Figure 6, which have been adjusted by introducing correction coefficient as stated in Section 2.3. After trial-and-error, the values of were determined to be 0.90 and 0.85, respectively, for 3 d and 28 d compressive strength. Eventually, fitting goodness R2 of 0.986 and 0.950 were obtained, respectively, in the case of 3 d and 28 d compressive strength, indicating that the ACE model has very high accuracy for both the cases. While probably due to the much bigger data size and the resulting increase in scattering, the ACE model for 28 d compressive strength is less accurate compared to the 3 d case, which could be seen from MAE and RMSE error evaluations as given in Table 6. However, the value is still acceptable for practical application.

In addition, from Figures 4 and 5, it can be seen that the relationship between (Na2O%) and is much clearer compared to the others, and with the change of (Na2O%), varies to a considerable scale. These indicate that among the five influencing factors, Na2O%, Ms, W/B, SA, and Kb, Na2O% took the largest weight-influencing compressive strength of AAS, followed by Ms > W/B > Kb > SA.

Similar to the ANN prediction, the ACE model was also run to predict the 3 d and 28 d compressive strength of the 11 AAS concrete mixes by inputting the material parameters used in the experimental work into the model. Comparison of the results measured and predicted is reported in Table 6, together with the error between them. Similarly, it is also found that no matter the mean error, max. error, or min. error, the value for 28 d compressive strength is lower than that for 3 d compressive strength, which indicates that compared to 3 d compressive strength, the ACE modeling was much more feasible to predict 28 d compressive strength, although the model for the latter was less accurate than the model for the former.

3.3. Comparison between ANN and ACE

From the discussion above it can be seen that when data size is relatively small, such as in the case of 3 d compressive strength (171 sets of data), although ANN model with a high precision could be established based on the given data size, its feasibility to predict nontrained data would be limited as ANN is a kind of statistics of known information in essence. While ACE is based on fitting technique, its function could be close to the true values by adjusting the coefficient of the input parameters. As a result, from Table 6 it can be seen that in the case of 3 d compressive strength, nearly half of the error values of ACE prediction is at a very low level and the average one is also acceptable although very few predictions have yielded a high error value.

When data size is sufficiently large, such as in the case of 28 d compressive strength (309 sets of data), the error value between the results predicted and measured dramatically decreases as shown in Table 7, for both ANN and ACE, indicating a great improvement of their feasibility for the prediction. When applying ANN, there are more samples to train the model more efficiently. When a prediction is being processed, the related memory could be aroused to give a precise prediction. On the other hand, an increase in data size is also helpful to improve the precision of fitting, as a result of which the ACE error value decreases as well. However, it seems that there is an optimum data size for the fitting of ACE. When the size is beyond the upper bound, the increase in fitting precision would be limited. Consequently, in this study, ACE performed not as well as ANN for the prediction of 28 d compressive strength, but it is worthy to note that the ACE error values are still in an acceptable range for practical application. Therefore, from the discussion above it can be seen that when massive data are available, it is more suitable to establish an ANN model for compressive strength prediction of AAS concretes. If data size was insufficient to support a precise ANN modeling, ACE could be considered as an alternative to yield an acceptable result.

4. Conclusions

Based on the data collected and modeling used in this study, the following conclusions could be drawn:(1)Both the ANN and the ACE models had adequate accuracy, no matter 3 d or 28 d compressive strength was considered.(2)Compared to the 3 d compressive strength, due to increased data scattering and data size, both the ANN and the ACE models did not yield a higher accuracy in the case of 28 d strength. However, also due to the increase in data size, both the models were more feasible to implement a 28 d strength prediction as a result of sufficient learning and training during modeling.(3)Based on the ACE analysis, among the five influencing factors, Na2O%, Ms, W/B, SA, and Kb,, Na2O% took the largest weight-influencing compressive strength of AAS, followed by Ms > W/B > Kb > SA.(4)When massive data are available, it is more suitable to establish an ANN model for compressive strength prediction of AAS concretes. If data size was insufficient to support a precise ANN modeling, ACE could be considered as an alternative to yield an acceptable result.

Data Availability

All the data used to support the finding of this study have been included in the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by National Natural Science Foundation of China (grant no. 52068038) and Yunnan Provincial Department of Science and Technology (grant no. 202101AT070089). The authors are grateful for the support of both institutions.