Modeling Slump of Ready Mix Concrete Using Genetically Evolved Artificial Neural Networks

Artificial neural networks (ANNs) have been the preferred choice formodeling the complex and nonlinearmaterial behavior where conventional mathematical approaches do not yield the desired accuracy and predictability. Despite their popularity as a universal function approximator and wide range of applications, no specific rules for deciding the architecture of neural networks catering to a specific modeling task have been formulated.The research paper presents a methodology for automated design of neural network architecture, replacing the conventional trial and error technique of finding the optimal neural network. The genetic algorithms (GA) stochastic search has been harnessed for evolving the optimum number of hidden layer neurons, transfer function, learning rate, and momentum coefficient for backpropagation ANN. The methodology has been applied for modeling slump of ready mix concrete based on its design mix constituents, namely, cement, fly ash, sand, coarse aggregates, admixture, and water-binder ratio. Six different statistical performance measures have been used for evaluating the performance of the trained neural networks. The study showed that, in comparison to conventional trial and error technique of deciding the neural network architecture and training parameters, the neural network architecture evolved throughGAwas of reduced complexity and provided better prediction performance.


Introduction
Cement concrete is one of the most widely used construction materials in the world today.The material modeling of concrete is a difficult task owing to its composite nature.Although various empirical relationships in the form of regression equations have been derived from experimental results and are widely in use, these do not provide accuracy and predictability wherein the interactions among the number of variables are either unknown, complex, or nonlinear in nature.Artificial neural networks (ANNs), touted as the next generation of computing, have been the preferred choice since the last few decades for modeling unstructured problems pertaining to material behavior.The notable applications of ANN in modeling properties of concrete are implementations in predicting and modeling compressive strength of high performance concrete [1], self-compacting concrete [2], recycled aggregate concrete [3], rubberized concrete [4], fibre reinforced polymer (FRP) confined concrete [5], durability of high performance concrete [6], predicting drying shrinkage of concrete [7], concrete mix design [8], and prediction of elastic modulus of normal and high strength concrete [9].
Known for its design mix precision leading to better quality concrete and ease of transportation and laying at the construction site, ready mix concrete (RMC) has emerged as a preferred concrete construction product catering to the requirement of the end user.One of the important properties of concrete that play an important role in the success of RMC is its workability.Workability of concrete is determined by the effort to lay and compact the freshly prepared concrete at construction site with minimum loss of homogeneity.The workability of concrete quantitatively measured in terms of slump value is an important quality parameter in RMC industry.The slump of concrete depends on the concrete's design mix proportion.As noticed for other properties of concrete, the slump value also exhibits a highly nonlinear and complex functional relationship with concrete's constituents.Recent studies, like prediction of slump and strength of ready mix concrete containing retarders and high strength concrete containing silica fume and plasticizers [10], predicting slump of fly ash and slag concrete [11], modeling slump of high strength concrete [12], modeling slump of high performance concrete [13], and modeling and analysis of concrete slump using laboratory test results [14], have proved the effectiveness of ANN in modeling slump of concrete.
Besides ANN applications in modeling behavior of concrete discussed above, there are many multidisciplinary applications of ANN which are beyond the scope of this paper.The reason for this rapid growth in the field of neural networks is attributed to its "black box" nature which allows it to be applied to almost every available problem, without seeking the knowledge of underlying relationships among the input and output variables.In spite of its popularity as a universal function approximator, the neural networks are still being designed using trial and error approach.Generally, iterative techniques using different combinations of number of hidden layers and hidden layer neurons are employed in conjunction with different learning rate, momentum coefficient, and transfer function to arrive at optimal neural network design.This technique of designing a neural network is therefore time consuming and relies heavily on the experience of the designer.In order to reduce the effort and time in designing the optimal neural network architecture and its training parameters, various studies for automatic design of neural network have been successfully performed in the past by harnessing the stochastic search ability of genetic algorithms [15][16][17][18][19].The autodesign of neural network using genetic algorithms (GA) has not been so far employed for modeling the slump of concrete.The research paper presents a methodology for evolving an optimal architecture of neural network and its training parameters using GA for modeling the slump of concrete based on concrete's design mix proportions.
The study has been organized into sections.Section 2 deals with data collection and its description.Section 3 deals with the methodology, in which trial and error neural network modeling of concrete slump, evolving neural network architecture, and training parameters using genetic algorithms and statistical performance measures have been discussed.Results, discussions, and conclusions have been dealt with in Sections 4, 5, and 6, respectively.

Collection of Data and Its Description
The exemplar data for ANN was collected from the same RMC plant to mitigate any chance of change caused in the slump data due to change in physical and chemical properties of the concrete design mix constituents.The collected data comprised 560 concrete design mix proportions and their corresponding slump test values.The design mix proportions included weight per m 3 of cement, pulverized fly ash (PFA), sand (as fine aggregate), coarse aggregate 20 mm, coarse aggregate 10 mm, admixture, and water-binder ratio.The range (maximum and minimum values) of the RMC data used in the study is shown in Table 1.In order to scale down all the inputs and outputs in a particular bound range preferably −1 to +1 or 0 to +1, data normalization is performed.This type of normalization has the advantage of preserving exactly all relationships in the data and it does not introduce any bias [20].This facilitates the learning speed, as these values fall in the region of sigmoid transfer function where the output is most sensitive to the variations of the input values [21].Linear scaling in the range −1 to +1 has been used in the present study having function where  norm is the normalized value of the variable  and  max and  min are the minimum and maximum values of variable , respectively.

Processing done by a neuron
Transfer function Figure 1: Mathematical model of an artificial neuron.

Neural Network Architecture and Training Parameters.
An artificial neural network is an information processing paradigm which presents a computational analogy inspired by human brain.ANN consists of processing elements called the artificial neurons which are arranged in layers.The computational structure of an artificial neuron comprises several inputs  1 ,  2 , . . .,   and an output .Every input is assigned a weight factor ( 1 ,  2 , . . .,   ) signifying the importance of the input.The synaptic weights carry both positive and negative values.A positive weight value leads to the forward propagation of information, whereas a negative value inhibits the signal.The neuron has an additional input called bias  and is interpreted as an additional weight.Bias has a constant value and its incorporation within the neuron helps in learning of the neural network.The inputs multiplied by corresponding weights are summed up and form the argument of the neuron's activation function.Figure 1 represents the mathematical model of an artificial neuron.
The architecture of a neural network consists of three basic layers denoted by "input layer, " "output layer, " and a number of intermediate "hidden layer/s".In case of multilayer feedforward neural networks (MFNN), the neurons in each layer are connected in the forward direction only and no intralayer connections between the neurons are permitted.The input features form the neurons of the input layer and output features are represented by the neurons of the output layer.The number of hidden layers and hidden layer neurons depends on the number of training cases, the amount of noise, and the degree of complexity of the function or the classification desired to be learnt [22].Hornik et al. [23] concluded that a three-layer feedforward neural network with backpropagation algorithm can map any nonlinear relationship with desired degree of accuracy.Some "rules of thumb" acting as initial guidelines for choosing neural network architecture have been suggested by [24][25][26][27].Nevertheless, the selection of hidden layers and hidden layer neurons is a trial and error process and generally started by choosing a network with a minimum number of hidden layers and hidden neurons.
MFNNs trained using backpropagation (BP) algorithms are commonly used for tasks associated with function approximation and pattern recognition.Backpropagation algorithm, in essence, is a means of updating neural network synaptic weights by backpropagating a gradient vector in which each element is defined as the derivative of an error measure with respect to a parameter [28].BP algorithm is a supervised learning algorithm wherein exemplar patterns of associated input and output data values are presented to the neural network during training.The information from input layer to output layer proceeds in forward direction only and with the help of BP algorithm, the error computed is passed backwards and the weights and biases are iteratively adjusted to enable learning of neural network by reducing network error to a threshold value.
A suitable learning rate and momentum coefficient are employed for efficient learning of the network.A higher learning rate leads to faster training, but by doing so it produces large oscillations in the weight change which may force the ANN model to overshoot the optimal weight values.On the other hand, a lower learning rate makes convergence slower and increases the probability of ANN model to get trapped in local minima.The momentum term effectively filters out the high frequency variations of the error surface in the weight space, since it adds the effect of the past weight changes on the current direction of movement in the weight space [29].A combined use of these parameters helps the BP algorithm to overcome the effect of local minima.The utility of transfer functions in neural networks is to introduce nonlinearity into the network.A consequence of the nonlinearity of this transfer function in the operation Advances in Artificial Neural Systems of the network, when so introduced, is that the network is thereby enabled to deal robustly with complex, undefined relations between the inputs and the output [30].

Evolving Neural Network Architecture and Training
Parameters Using Trial and Error.In the present study, the input layer consists of seven neurons, namely, cement, fly ash, sand, coarse aggregate (20 mm), coarse aggregate (10 mm), admixture content, and water-binder ratio.The output layer comprises a single neuron representing the slump value corresponding to the seven input neurons defined above.In this study eleven single hidden layer neural network architectures of different complexities with hidden layer neurons varying in the range 5 to 20 have been used for evolving the optimal neural network architecture.The neural network architecture with five hidden layer neurons for the present study is shown in Figure 2. The learning rate and momentum coefficient have been varied in the range 0 to +1.Hyperbolic tangent and log-sigmoid transfer functions have been used in the hidden layers along with linear transfer function in the output layer.Different combinations of learning rate and momentum coefficient with hidden layer transfer function have been tried for effective training of neural networks.The neural networks were trained using the training data set.The information to the neural network was presented through input layer neurons.The information propagated in the forward direction through hidden layers and was processed by the hidden layer neurons.The network's response at the output layer was evaluated and compared with actual output.The error between the actual and the predicted response of the neural network was computed and propagated in the backward direction to adjust the weights and biases of the neural network.Using the BP algorithm the weights and biases were adjusted in a manner to render error to a minimum value.In the present study, Levenberg-Marquardt backpropagation algorithm which is the fastest converging algorithm preferred for supervised learning has been used as training algorithm.During training process, the neural network has a tendency to overfit the training data.This leads to poor generalization when the trained neural network is presented with unseen data.The validation data set is used to test the generalization ability of trained neural at each iteration cycle.Early stopping of neural network training is generally undertaken to avoid overfitting or overtraining of the neural network.In this technique the validation error is also monitored at each iteration cycle along with training error and the training is stopped once the validation error begins to increase.The neural network architecture having the least validation error is selected as the optimal one.

Evolving Neural Network Architecture and Training
Parameters Using Genetic Algorithms.Genetic algorithm (GA) inspired by Darwin's theory "survival of the fittest" is a global search and optimization algorithm which involves the use of genetic and evolution operators.GA is a population based search technique that simultaneously works on a number of probable solutions to a problem at a time and uses probabilistic operators to narrow down the search to the region where there is maximum possibility of finding an optimal solution.GA presents a perfect blend of exploration and exploitation of the solution search space by harnessing computational models of evolutionary processes like selection, crossover, and mutation.GAs outperform  the efficiency of conventional optimization techniques in searching nonlinear and noncontinuous spaces, which are characterized by abstract or poorly understood expert knowledge [31].
Automatic design of neural network architecture and its training parameters is accomplished by amalgamating GA with ANN during its training process.The methodology uses GA to evolve neural network's hidden layer neurons and transfer function along with its learning rate and momentum coefficient.The BP algorithm then uses these ANN design variables to compute the training error.The training process is monitored at each iteration by computing the validation error.The training of neural network is stopped once validation error starts to increase.This process is repeated number of times till optimum neural network architecture and its training parameters are evolved.The steps of this methodology are presented as flow chart in Figure 3 and are summarized as below.

Advances in Artificial Neural Systems
(1) Initialization of Genetic Algorithm with Population of Chromosomes.In GA the chromosomes contain the information regarding the solution to a problem.The chromosomes contain number of genes which represent potential solutions to a problem.Being a population based heuristic, the GA starts with a number of chromosomes representing initial guesses to the possible solutions to a problem.In the present study, the chromosomes represent the information regarding number of hidden layer neurons, transfer function, learning rate, and momentum coefficient.The numbers of hidden layer neurons are initialized in the range 5 to 20 neurons and are represented as discrete integer number.The GA is allowed to select the transfer function between tangent hyperbolic and log-sigmoid transfer functions.The range of learning rate and momentum coefficient is bounded in the range 0 to +1 and is represented by real numbers.
The size of the population is chosen in such a way to promote evolving of optimal set of solutions to a particular problem.A large initial population of chromosomes tends to increase the computational time, whereas a small population size leads to poor quality solution.Therefore population size must be chosen to derive a balance between the computational effort and the quality of solution.In the present study an initial population size of 50 chromosomes is used.
(2) Creating the Neural Network.The neural network is created using randomly generated hidden layer neurons.As discussed above, the number of input layer neurons is 7 and output layer neurons is 1; hence the neural network architecture is initialized as 7-N-1, where N is the number of hidden layer neurons determined using GA.The transfer function for hidden layers is also initialized using GA.
(3) Training of Neural Network and Evaluating Fitness of Chromosomes.The neural network is trained using Levenberg-Marquardt backpropagation training algorithm.The learning parameters, namely, learning rate and momentum coefficient, initialized using GA are used during the training process for systematic updating of neural network weights and biases.BP algorithm updates the weight and biases through backpropagation of errors.Early stopping technique is applied to improve the generalization of the neural network.The fitness function acts as measure of distinguishing optimal solution from numerous suboptimal solutions by evaluating the ability of the possible solutions to survive or biologically speaking it tests the reproductive efficiency of chromosomes.The fitness of each chromosome is computed by evaluating the root mean square error (RMSE) using where   and   denote the target or observed values and ANN predicted concrete slump values, respectively.
(4) Selecting the Fitter Chromosomes.The evolution operator of GA called as selection operator helps in selecting the fitter chromosomes based on the value of the fitness function.
A selection operator performs the function synonymous to a filtering membrane, allowing the fitter chromosomes to survive to create new offspring.As one moves from one generation of population to the next, the selection operator gradually increases the proportion of the fitter chromosomes in the population.The present study uses roulette wheel selection strategy which allows probability of selection proportional to the fitness of the chromosome.The basic advantage of roulette wheel selection is that it discards none of the individuals in the population and gives a chance to all of them to be selected [32].
(5) Creating Next Generation of Population.The new generation of population is created through two genetic operators of GA called crossover and mutation.The crossover operator is a recombination operator and is applied to a pair of parent chromosomes with the hope to create a better offspring.This is done by randomly choosing a crossover point and copying the information before this point from the first parent and then copying the information from the second parent beyond the crossover point.The present study utilized the scattered crossover with probability 0.9 for recombining the two parent chromosomes for producing a fitter child.The mutation operator modifies the existing building blocks of the chromosomes maintaining genetic diversity in the population.It therefore prevents GA from getting trapped at a local minimum.In contrast to crossover which exploits the current solution, the mutation aids the exploration of the search space.Too high mutation rate increases the search space to a level that convergence or finding global optima becomes a difficult issue whereas a lower mutation rate drastically reduces the search space and eventually leads genetic algorithm to get stuck in a local optima.The present study uses uniform mutation with mutation rate 0.02.The procedure for creating new population of chromosomes is continued till maximum generation limit is achieved or the fitness function reaches a saturation level.Maximum number of generations used for present study is 150.

Evaluating Performance of the Trained Models.
The study uses six different statistical performance metrics for evaluating the performance of the trained models.The statistical parameters are mean absolute error (MAE), root mean square error (RMSE), mean absolute percentage error (MAPE), coefficient of correlation (), Nash-Sutcliffe efficiency (), and root mean square to standard deviation ratio (RSR).The above performance statistics were evaluated using ) , where

Results
The neural network architecture was evolved through trial and error process by analyzing 30 different combinations of hidden layer neurons, transfer function, learning rate, and momentum coefficient.The optimal neural network architecture (BPNN) was evolved as 7-11-1 having eleven hidden layer neurons with learning rate 0.45, momentum coefficient 0.85, and tangent hyperbolic hidden layer transfer function.The same operation was performed by incorporating GA during the training of ANN.The GA was able to evolve the optimal neural network architecture and training parameters in 92 generations (Figure 4).The time taken by GA to reach the saturation level of fitness function 2.2357 mm was evaluated as 305.3412 seconds.The GA evolved neural network architecture (ANN-GA) comprised 9 hidden layer neurons and tangent hyperbolic transfer function.The optimal learning rate and momentum coefficient for backpropagation neural network were computed as 0.3975 and 0.9385, respectively.Both ANN-GA and BPNN models subsequent to training were validated and tested.The results in terms of the performance statistics are presented in Table 2.
The entire RMC data was also used for evaluating the prediction ability of the trained models, namely, BPNN and ANN-GA.The regression plots showing the prediction of trained BPNN and ANN-GA models are exhibited in Figures 5(a) and 5(b), respectively.The statistical performance for the entire data set is tabulated at Table 3.

Discussions
The results of the study show that amalgamation of GA with ANN during its training phase leads to evolving of optimal neural network architecture and training parameters.In comparison to trial and error BPNN neural network having architecture 7-11-1, the hybrid ANN-GA automatically evolved a less complex architecture 7-9-1.Moreover the   to trained BPNN model.Overall, the performance metrics show that the ANN-GA model has consistently outperformed the BPNN model.

Conclusions
The study presented a methodology of designing the neural networks using genetic algorithms.Genetic algorithms population based stochastic search was harnessed during the training phase of the neural networks to evolve the number of hidden layer neurons, type of transfer function, and the values of learning parameters, namely, learning rate and momentum coefficient for backpropagation based ANN.The performance metrics show that ANN-GA model outperformed the prediction accuracy of BPNN model.Moreover the GA was able to automatically determine the number of hidden layer neurons which were found to be less than those evolved using trial and error methodology.The hybrid ANN-GA provided a good alternative over time consuming conventional trial and error technique for evolving optimal neural network architecture and its training parameters.
The proposed model based on past experimental data can be very handy for predicting the complex material behavior of concrete in quick time.It can be used as a decision support tool, aiding the technical staff to easily predict the slump value for a particular concrete design mix.This technique will considerably reduce the effort and time to design a concrete mix for a customized slump without undertaking multiple trials.Despite the effectiveness and advantages of this methodology, it is also subjected to some limitations.Since the mathematical modeling of concrete slump is dependent on the physical and chemical properties of the design mix constituents, hence the same trained model may or may not be applicable for accurate modeling of slump on the basis of design mix data obtained from other RMC plants deriving its raw material from a different source.

Figure 2 :
Figure 2: Single hidden layer neural network with five hidden layer neurons.

Figure 3 :
Figure 3: Evolving neural network architecture and training parameters using GA.

Figure 5 :
Figure 5: Regression plot of BPNN and ANN-GA predicted slump versus observed slump.

Table 1 :
Range of RMC data used for neural network modeling.
and   denote the target or observed values and ANN predicted values and  and  represent the mean observed and mean ANN predicted values, respectively. represents the total number of data.A lower value of MAE, RMSE, MAPE, and RSR indicates good performance of the model.A higher value of  and  statistics above 0.90 indicates good prediction of the model.

Table 2 :
Statistical performance of ANN models for training, validation, and testing data sets.

Table 3 :
Statistical performance of the trained ANN models for the entire data set. and  during training, validation, and testing of trained models.The performance statistics computed for the entire data set using the trained ANN-GA model shows a lower MAE, RMSE, MAPE, and RSR value of 1.5754 mm, 2.3345 mm, 1.0841%, and 0.1902, respectively, and higher  and  values of 0.9819 and 0.9638, respectively, in comparison Advances in Artificial Neural Systems