Abstract

The objective of this work was to examine the compressive strength behavior of ground bottom ash (GBA) concrete by using an artificial neural network. Four input parameters, specifically, the water-to-binder ratio (WB), percentage replacement of GBA (PR), median particle size of GBA (PS), and age of concrete (AC), were considered for this prediction. The results indicated that all four considered parameters affect the strength development of concrete, and GBA with a high fineness can act as a good pozzolanic material. The optimal ANN model had an architecture with two hidden layers, with six neurons in the first hidden layer and one neuron in the second hidden layer. The proposed ANN-based explicit equation represented a highly accurate predictive model, for which the statistical values of R2 were higher than 0.996. Moreover, the compressive strength behavior determined using the optimal ANN model closely followed the trend lines and surface plots of the experimental results.

1. Introduction

Artificial neural networks (ANNs) have been widely applied owing to the excellent performance of the associated high-accuracy predictive model in learning and analyzing the effects of the input and output variables. In particular, ANN models are valuable tools for predicting experimental results. In the field of concrete engineering, for the experimental datasets of concrete properties, an ANN model was proposed to predict the compressive strength of concrete containing pozzolanic materials such as fly ash, silica fume, metakaolin, and ground granulated blast furnace slag [13]. In addition, the ANN model has been extended for predicting the properties of concrete, for instance, workability, corrosion currents, split tensile strength, water permeability, and chloride permeability [4]. Many studies have reported that ANN models lead to a more accurate and precise prediction than that obtained using linear and nonlinear regression techniques, and it can be easily expanded to new additional databases, enabling the retraining of the network [24].

Bottom ash (BA) is a byproduct of the combustion of powder coal in an electric power plant. In particular, fly ash (FA) is melted at high temperatures, and it agglomerates to form BA at low temperatures. BA is similar to FA; however, BA has a considerably larger particle size, and thus it cannot be used as a cementitious material in concrete. Ghafoori and Cai [5, 6] studied BA as a fine aggregate and reported the behavior of laboratory-made roller-compacted concrete (RCC) containing BA. The RCC containing BA exhibited excellent strength, stiffness, durability, and deformation properties. Furthermore, Ghafoori and Bucholc [7, 8] used dry BA as a fine aggregate in concrete and indicated that the amount of required mixing water increased rapidly with the increase in the BA content. Other researchers used BA with low density as an aggregate in lightweight concrete to decrease the concrete density.

Compressive strength is one of the most important properties of mortar or concrete. Many researchers [9, 10] studied the mechanism analysis of hydration and pozzolanic reaction which found that when cement and pozzolanic material were mixed with water, cement clinker minerals hydrate first, to release Ca(OH)2 crystal, producing the alkalinity solution. Then, the active compound of pozzolanic material reacted with Ca(OH)2 to produce the calcium silicate hydrate and calcium aluminate hydrate, resulting in reducing of Ca(OH)2 content in the solution; thus the hydration of cement clinker was accelerated. The later subsequent hydration reaction equations could be expressed as

In the formula, if x ≤ 3,

Equations (1)–(3) express the pozzolanic reaction, in which both C-S-H and C-A-H enhance the compressive strength of concrete. The fineness of pozzolanic materials is one of the most important physical properties that affect the pozzolanic reaction. Therefore, increasing the fineness of the BA can help enhance the quality of concrete. Consequently, in this research, three sizes of GBA were considered to investigate the corresponding influence on the compressive strength of concrete.

An investigation of the original BA [9] indicated that the median particle size (d50) and percentage of retained particles on sieve number 325 were approximately 300 μm and 95%, respectively. The strength activity indices of the original BA mortar at the ages of 7 and 28 days were 46% and 49% of the standard mortar, respectively, and its water requirement was 114%. In its original form, the BA cannot be used as a pozzolanic material to replace Portland cement in concrete mixtures owing to its large particle size (compared to that of cement) and porous structure, as these aspects may lead to an undesirable degradation of the properties of fresh and hardened concrete. However, several studies [4, 10] reported that BA could be ground, thereby obtaining the necessary properties for BA to be used as partial replacement of cement in concrete mixtures. In particular, the decrease in the particle size of BA increases its pozzolanic activity and helps improve the mechanical and physical properties of concrete.

To this end, the objective of this work was to examine the effect of the water-to-binder ratio (WB), percentage replacement of GBA (PR), median particle size of GBA (PS), and age of concrete (AC) on the compressive strength of GBA concrete. Furthermore, on the basis of this experimental data, a neural network was established to derive an explicit equation to predict the compressive strength of GBA concrete.

2. Experimental Database

2.1. Materials

Experimental ASTM type-I ordinary Portland cement and BA from Mae Moh power plant in northern Thailand were employed in this work. The chemical compositions of the cementitious materials are presented in Table 1, and their physical properties are listed in Table 2. The specific gravity of cement was 3.14, with a median particle size (d50) of 11.3 µm. The portion of the major oxide (SiO2 + Al2O3 + Fe2O3) of the BA was 78.97%, which is higher than the minimum requirement of 70% specified by ASTM C 618 [11] for class F (as FA). To modify the properties of BA, it was ground using a ball mill until the GBA exhibited three sizes: small (SB), medium (MB), and large (LB). The corresponding median particle sizes were 6.3, 13.6, and 22.2 µm, and the size of the particles retained on sieve number 325 were 4.9%, 18.3%, and 33.6% by weight, respectively. The fine and coarse aggregate (crushed limestone) had specific gravity values of 2.58 and 2.61 with fineness moduli of 3.1 and 7.2, respectively. The maximum size for the coarse aggregate was 20 mm.

2.2. Mixture Proportions

The control concretes were designed with water to obtain cement ratios of 0.43, 0.48, 0.55, 0.62, 0.70, and 0.80. Portland cement was replaced by different sized GBA in proportions of 0%, 10%, 20%, 30%, 40%, and 50% by weight of the binder. In all the cases, the concrete mix included 195 kg/m3 of water as a constant value, and the workability was maintained by using a slump of fresh concrete with values ranging from 80 to 120 mm. The mix proportions of all concrete mixes are summarized in Table 3.

2.3. Compressive Strength Test

The concrete specimens were cast in 100 × 200 mm cylindrical molds. Then, these specimens were molded for 24 h and cured in water at 23 ± 2°C until the testing age. In all the mixtures, the tests were performed on triplicate specimens, the average values were reported, and all the samples were tested for compressive strength at 3, 7, 14, 28, 60, 90, and 180 days. The compressive strength of the specimen was calculated by dividing the maximum load attained during the test by the cross-sectional area of the specimen, as specified in ASTM C39/C39-12a [12].

3. Artificial Neural Network

The ANN can be defined as a mathematical model of reasoning, and it was biologically inspired and copied based on the human brain, whose architecture is commonly a basis of the feed-forward neural network (also known as a multilayer perceptron) with an input layer, one or more hidden layers, and an output layer. The learning rules of an ANN (training algorithm) are most commonly established using the gradient descent training algorithm (also known as the backpropagation technique) by modifying the weights and biases according to the error from a network to adjust the actual outputs near the target outputs. In such a method, each node receives weighted inputs from other nodes and communicates its outputs to other nodes by using an activation function. During the learning process, the nodes in each layer are connected to those in the immediate next layer via a cyclic link. The connection links that provide the weights are usually designed to minimize the error in a training algorithm. Finally, the target values of this training algorithm produce the output values of the predicted model. For more detailed information regarding the ANN, the readers can refer to textbooks [13, 14] and the existing literature [15, 16]. Nevertheless, for convenience, the explicit mathematical equation with the results of the output, weight, and bias values of the ANN model used in this study is described in the subsequent section.

3.1. Evaluation of Model Performance

Statistics-based methods [1719] can be used to compare the predictive performance of several models and select a robust model. Herein, the model prediction accuracy was determined by calculating the absolute fraction of variance (R2), root mean square error (RMSE), and mean absolute percentage error (MAPE). R2 is a measure of the proportion of the total variation in the dataset; R2 lies between 0 and 1, where a value of 1 corresponds to a higher predictive power of the model. The RMSE is the average distance of a data point from the expected value. The MAPE is a measure of the prediction accuracy of the model. The RMSE and MAPE range from 0 to ∞, and lower values indicate a better mode fit. However, the MAPE cannot be used if there exist values with no error, involving division by zero occurs. The three statistics measures are expressed as follows:where n is the number of experimental dataset, M is the measured value, and T is the target value.

3.2. ANN Architecture Design

The neural network models are developed using MATLAB software. Therefore, there are several learning algorithms, activation functions, and network performances in this software. The performance of the ANN model primarily depends on the network architecture and parameter settings. In this study, the parameter settings in the ANN model are shown in Table 4, and Table 5 shows the range values and numerable parameters of the experiment and prediction. Four factor parameters constitute the input layer: the water-to-binder ratio (WB of 0.43, 0.48, 0.55, 0.62, 0.70, and 0.80), percentage replacement of ground bottom ash (PR of 0, 10, 20, 30, 40, and 50%), median particle sizes of ground bottom ash (PS of 6.3, 13.6, and 22.2 μm), and ages of concrete (AC of 3, 7, 14, 28, 60, 90, and 180 days); the compressive strength of ground bottom ash (GBA) concrete is the output layer. The data in the ANN models were obtained from laboratory experiments and trained and tested using 631 and 125 sets of experimental data for the model, respectively.

Bayesian regularisation (BR) learning algorithm was used to perform the training process for the ANN model in the MATLAB software. Meanwhile, the network performance function of the sum squared error (SSE) was measured to search the best technique based on the statistic value of R2 to minimize the error between the actual output and the target output [14]. Log-sigmoid, tan-sigmoid, and purelin linear activation functions were used to evaluate their performance against the model mean squared values [20]. There activation functions are defined in

The optimal ANN model was determined among the architectures that consisted of one hidden layer with 4–14 neurons or two hidden layers with 4–14 neurons in the first hidden layer and 1-2 neurons in the second hidden layer. The optimal ANN model was obtained with high performance in 3 steps. First, the high-accuracy predictive model is guaranteed by both training and testing datasets with statistical values of R2 higher than of 0.99. Second, an explicit model equation is reported in this paper. Finally, the reliable compressive strength behavior is confirmed using the sensitivity analysis illustrated by two-dimension (2D) plots and three-dimension (3D) plots to compare the experimental and predictive results. This result shows that the optimal ANN model predicts the concrete compressive strength that closely matches and follows the trend lines and surface plots of the experimental results.

4. Results and Discussion

The relationships of the four factor parameters (WB, PR, PS, and AC) that can potentially affect the compressive strength (CS) of GBA concrete were considered as the four-parameter input and one-parameter output variables in the ANN-modeled analysis. The amount and range values of these experimental parameters are presented in Table 5.

4.1. Experimental Results

All the 177 graphs of 3D plots were established using 756 sets of experimental data. The number of experimental 3D plots can be calculated as follows:where N3D is the number of 3D plots and the numbers of WB, PR, PS, and AC parameters are 6, 6, 3, and 7, respectively.

The effects of the four input parameters and the response output parameter on the surface plots relative to the behaviors of the concrete compressive strength are shown by the 3D plots, each of which illustrates the effects of two factor parameters on each surface plot of the concrete strength. For example, Figure 1 shows six samples of 3D plots from all 177 experimental graphs with WB = 0.43, PR = 0, PS = 6.3, and AC = 3, where two parameters are fixed, and two parameters are independent. The compressive strength increases when WB decreases at all values of AC, PR, and PS (Figures 1(a)1(c)). In Figures 1(b), 1(d), and 1(f), when AC increases, the compressive strength also increases for all values of WB, PR, and PS. In Figures 1(c), 1(e), and 1(f), decreasing the PS increases the strength at all levels of AC, PR, and WB. However, the graphs of the 3D plots are not sufficiently smooth owing to the use of insufficient data. The use of more information can further smoothen the graph. These findings indicate that the ANN model can be used to predict the compressive strength and better clarify the relationships among all the factor parameters.

Furthermore, it was observed that, at the same WB and replacement of GBA, the concrete mixed with the coarser GBA exhibited a lower compressive strength than that of the concrete with the finer GBA. As shown in Figure 1(d), the compressive strength of GBA concrete at early ages (before 28 days) depended on the percent replacement of the GBA. A lower replacement of GBA yielded a higher compressive strength of concrete. This result is similar to that reported in other studies [9, 10], in which original and ground BA were used to replace Portland cement. For high fineness GBA, at a later age, the concretes with a higher percent replacement of GBA produced a higher compressive strength than those with lower replacement and the control concretes. Moreover, the results suggested that a higher cementitious material content of GBA concrete corresponded to a higher development rate of compressive strength, when the cementitious material contents were in the range of 244 to 454 kg/m3, especially at the early ages (before 28 days). These results further indicated that the fineness and percent replacement of GBA are important factors that affect the compressive strength of concrete, and GBA with a high fineness can be used as a pozzolanic material for concrete.

4.2. ANN Model Performance and Equation

The ANN models were developed by using the command line code in the command window of the MATLAB software. A graphical user interface was applied to create 2D and 3D plots by using the experimental results and predictive model values in MATLAB. In this study, 83.5% and 16.5% of the experimental dataset were randomly selected and used to train (631 sets of experimental data) and test (125 sets of experimental data) the model, respectively. The formats of the variables were run by the ANN model and are presented in Tables 6 and 7 for the input and output datasets, respectively.

The optimal ANN model was determined by investigating the model performance. The results showed that the optimal ANN model had two hidden layers with six neurons in the first hidden layer and one neuron in the second hidden layer. The optimal ANN model was used in the following manner. First, the high-accuracy predictive model was developed by considering both training and testing datasets with statistical values of R2 higher than of 0.99. Second, an explicit model equation was established, as reported in (9)–(11). Finally, the reliable compressive strength behavior was validated by performing the sensitivity analysis using two-dimensional (2D) and three-dimensional (3D) plots to compare the experimental and predictive results. In addition, all the ANN model architectures exhibited R2 values higher than 0.94. The number of hidden layers is not related to R2. Although the ANN model provided a high accuracy, the plotted graphs were not consistent with the concrete behavior. Therefore, when choosing the ANN model architecture, the behavior of concrete must be considered. Figure 2 shows a comparison of the 2D plots of the experimental results and predictive results of the ANN model architectures with one hidden layer and R2 higher than of 0.99. Figure 2(b) (one hidden layer with six neurons) shows that the trend lines more closely match the experimental lines than those shown in Figure 2(a) (one hidden layer with five neurons) and Figure 2(c) (one hidden layer with 14 neurons). Therefore, the ANN model architecture with one hidden layer with six neurons is the best model for one hidden layer. Figure 3 shows a comparison of the 2D plots of the experimental results and predictive results of the ANN model architectures with two hidden layers and R2 higher than 0.99. In Figure 3(b) (two hidden layers with six neurons in the first hidden layers and one neuron in the second hidden layers), the trend lines more closely match the experimental lines than those in Figure 3(a) (two hidden layers with six neurons in the first hidden layer and zero neurons in the second hidden layer) and Figure 3(c) (two hidden layer with six neurons in the first hidden layer and two neurons in the second hidden layer). Therefore, the ANN model architecture with two hidden layers with six neurons in the first hidden layer and one neuron in the second hidden layer is the optimal ANN model, which is used in this study.

The parameter settings in the optimal ANN models are presented in Table 8, and the model architecture is shown in Figure 4. The optimal ANN model exhibits an acceptable performance and high accuracy because the absolute fraction of variance (R2) is higher than 0.99, and the RMSE and MAPE are lower than 3.4569 for both the training and testing datasets, as shown in Table 9. The weight matrices and bias vectors of the network were extracted to adjust the parameters, as presented in Table 10. Therefore, the explicit ANN-modeled formulation of the compressive strength of GBA concrete can be mathematically expressed as follows:where X1, X2, …, X6 denote the summations of the weighted inputs with the addition of the bias values of the first hidden layer. Next, the log-sigmoid activation function in the first hidden layer is as follows:

Finally, the tan-sigmoid activation function in the second hidden layer and the sum of the output layer is

An ANN predicted model was thus developed to establish the reliable behavior and accurate estimation of the compressive strength of GBA concrete. The results indicated that the optimal ANN model, which has two hidden layers with six neurons in the first hidden layer and one neuron in the second hidden layer, is represented by four factor parameters (WB, PR, PS, and AC). In addition, the log-sigmoid and tan-sigmoid activation functions were employed in the first and second hidden layers, respectively. The BR was used as the learning algorithm, and the sum square error was used to evaluate the network performance. A high-accuracy predictive model with acceptable error is obtained, with the statistical values of R2, MAPE, and RMSE being 0.9969, 2.4586, and 3.339, respectively, for the training dataset and 0.9968, 2.6417, and 3.4569, respectively, for the testing dataset.

4.3. Sensitivity Analysis

The response surface method is among the most widely used approaches for concrete mixture design optimization [2125] owing to its simplicity in illustrating the response-related equations to analyze the relationship between the required factors and the required response by using 3D plots. In addition, a few studies [15, 16] pertaining to the field of concrete property predictive models popularized the use of the RMS to represent the reliable behavior pertaining to the concrete property in input-output predictive relationships. Although [26] discussed that the concrete and microstructure-property relationships are at the heart of modern material science and the analytical methods of material science can model and predict the microstructurally stable behavior, homogeneous materials do not appear to follow this trend in the case of concrete structures. Such a case also pertains to the reliable compressive strength relationship of the GBA concrete behavior. Therefore, the 2D and 3D plots from the sensitivity analysis were used to illustrate the reliable relative behavior of the measured values and predicted values for comparison. The 3D plots of the ANN-based prediction results show the response-related equations of the surface plot of the ANN model. The 2D plots show the relationship of a factor variable and a response variable in the form of the trend lines of the measured and predicted values, which involve a side view of the 3D plots and can clearly illustrate the reliable compressive strength behavior of the predictive model.

Table 5 presents the value of each input parameter at various levels in the model and the output concrete compressive strength. The input datasets of the predictive model were as follows: WB = [0.43 : 0.01 : 0.80], PR = [0 : 2 : 50], PS = [6.3 : 0.3 : 22.2], and AC = [3 : 3 : 180] in the form of [m :i: n], where m, i, and n denote the initial, incremental, and final values of the input variable, respectively. In total, 3,201,120 sets of input variable data were calculated from 38, 26, 54, and 60 variables of WB, PR, PS, and AC, respectively. In addition, the sensitivity analysis findings to compare the measured and predicted values were represented in the form of 2D and 3D plots.

Figure 5 shows six samples of 3D plots from all 177 predictive graphs of the optimal ANN models with identical parameters as shown in Figure 1. The relationship of the parameters is clearer, and the graph has smoother characteristics. All the 3D predictive plots shown in Figure 5 correspond to the surface plots of the experimental results in Figure 1, which have identical parameters. Furthermore, the 2D plots guarantee the reliable predictive behavior concerning the side view of the 3D plots, as shown by the samples shown in Figure 6. The predictive values of the optimal ANN model are remarkably close to the experimental results, which indicate a strong correlation between the predictive and measured values. In addition, it is demonstrated that the optimal ANN model can provide reliable intelligence and accurate analysis because the predictive results are better correlated with the plotting surface of the 3D plots and trend lines of the 2D plots than the experimental results, as shown in Figures 1(b), 5(b), 6(c), and 6(d). However, the remaining 2D and 3D plots can be validated by using the explicit equation of the optimal ANN model by creating the 2D and 3D plots with this equation. Moreover, it can be observed that, at the same WB ratio and the same percent replacement, concrete mixed with coarser GBA produces a compressive strength lower than that of concrete with the finer GBA. The compressive strengths of the GBA concretes at the early ages (before 28 days) depend on the percent replacement. A lower percent replacement corresponds to a higher compressive strength of concrete than that for the case with the higher percent replacement. This result was similar to the findings of existing studies [27, 28].

5. Conclusions

The following conclusions were drawn from the abovementioned investigation:(1)The compressive strengths of GBA concretes at the early ages depend on the percent replacement. A lower percent replacement of GBA yields a higher compressive strength of concrete compared to that in the case with a higher percent replacement. At the latter ages, the compressive strength depends on the fineness of the GBA; in particular, concrete mixed with coarser ground bottom ash yields a compressive strength lower than that of the concrete with the finer GBA.(2)The WB, PR of GBA, PS of GBA, and AC are the major parameters affecting the strength development of concrete, and a GBA with a high fineness can be used as a good pozzolanic material.(3)The optimal ANN model has an architecture with two hidden layers with six neurons in the first hidden layer and one neuron in the second hidden layer. The log-sigmoid and tan-sigmoid activation functions are used in the first and second hidden layers, respectively, BR is used as the learning algorithm, and the SSE is used to evaluate the network performance. The high-accuracy predictive model has the statistical values of R2, MAPE, and RMSE as 0.9969, 2.4586, and 3.339, respectively, for the training dataset and 0.9968, 2.6417, and 3.4569, respectively, for the testing dataset.(4)The results conclusively prove that the proposed ANN-based-explicit equations can reasonably generalize the actual output and the target output. The 2D and 3D plots of the compressive strength behavior from the optimal ANN model closely matched and correspondingly followed the trend lines and surface plots of the experimental results.

Data Availability

The data used to support the findings of this study are included within the article. Any additional data related to the paper may be requested from the corresponding author.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research received financial support from Faculty of Engineering, Mahasarakham University. In addition, the authors thank Mr. Anon Chuachun from the Concrete Research Laboratory (CRL) for his valuable assistance in the exhaustive experimental work.