Flexural Strength of Concrete Beams Made of Recycled Aggregates: An Experimental and Soft Computing-Based Study

: The implementation of recycled concrete aggregates (RCAs) in the construction industry has been highlighted in the literature recently. This study aimed to propose an intelligent model for predicting the ultimate ﬂexural strength of recycled reinforced concrete (RRC) beams. For this reason, a database comprising experimental tests on concrete beams was compiled from the literature. Additionally, two experimental tests were performed in the laboratory to enhance the aforementioned database. The ﬂexural test results showed a 10% reduction in ﬂexural strength when the RRC beam was tested instead of a conventional beam (constructed with natural aggregates). Nevertheless, an artiﬁcial neural network (ANN) improved by particle swarm optimization (PSO), as well as an imperialist competitive algorithm (ICA), were utilized for developing the predictive model. The inputs of the hybrid predictive models of ﬂexural strength were the beam geometrical properties, reinforcement ratio, RCA percentage, compressive strength of concrete, and the yield strength of steel. The overall ﬁndings (e.g., correlation coefﬁcient values of 0.997 and 0.994 for the testing data) showed the feasibility of the PSO-based ANN predictive model, as well as the ICA-based ANN predictive model in the ﬂexural assessment of RRC beams. Furthermore, comparing the prediction performances of PSO-based ANN with ICA-based ANN and the conventional ANN showed that the PSO-based ANN model outperformed the predictive model built with the conventional ANN and the ICA-ANN.


Introduction
The space limitation for construction in urban areas, waste disposal demolition, and the high renting costs of landfills is becoming a major problem. Many researchers underlined the beneficial effect of recycling concrete in minimizing the aforementioned problems to some extent [1]. The demand for new aggregates can be reduced considerably if recycled concrete aggregates (RCAs) are used. In other words, the utilization of RCA can lead to the preservation of resources, forested areas, and/or river beds. In fact, the beneficial effects of RCA utilization from the environmental point of view are considerable. Hossain et al. [2] mentioned that the implementation of coarse RCA in Hong Kong reduced the greenhouse gases footprints considerably. Many researchers investigated the feasibility of using RCA in the construction industry. Soutsos et al. [3] reported the workability of waste construction materials as aggregates for producing precast concrete products. Silva et al. [4] proposed a new method for using construction waste to produce concrete using a statistical analysis. In countries such as Japan, a considerable amount of concrete is recycled [5]. However, Kang et al. [6] mentioned that almost 14% of waste concrete is recycled.
Yehia et al. [7] conducted a comprehensive study on the effectiveness of using recycled aggregates for concrete production. Katerusha [8] reported that in 2018, the average recovery rate of construction in the European Union was about 90%. In Zurich, 90% of the utilized concrete between 2005-2018 comprised recycled aggregates. Although the use of RCA in a Kamanli et al. [32] proposed an ANN-based predictive model for assessing the flexural behavior of reinforced concrete beams. The coefficient of determination (R 2 ) equal to 0.96 for testing data suggested that their recommended predictive model was good enough. It is worth mentioning that the input parameters of their model comprised the reinforcement ratio, compressive strength of concrete, cement dosage, elastic modulus of concrete, yielding stress of steel, ultimate strength of steel, and load. The deflection was set as the model output in their ANN-based predictive model. They used 82 experimental sets of data to train their intelligent model. Perera et al. [33] utilized ANN to predict the ultimate shear strength of reinforced concrete beams.
Kaczmarek and Szymanska [34] also highlighted the workability of ANN in predicting the deflection of reinforced concrete beams. They used 293 sets of data for their model development. The input parameters of their model included Young's modulus of concrete, Young's modulus of steel, the surface area of tension reinforcement, and the bending moment. The geometry of the beam was omitted from their study. Bai et al. [35] promoted implementations of ANN, support vector machine (SVM), and adaptive neuro-fuzzy inference system (ANFIS) for assessing the deflection of reinforced concrete beams. They used 120 sets of data for their model development. The input parameters of their proposed models were the tensile strength, beam length, compressive strength, shear span-depth ratio, applied load, percentage of stirrups, and percentage of reinforcement.
Suguna et al. [36] proposed an ANN-based predictive model for the flexural assessment of high-strength concrete. They used fiber-reinforced polymer (FRP) laminates for concrete strengthening. The input parameters of their developed model were the beam geometrical properties (length, depth, width), reinforcement ratio, concrete compressive strength, tensile strength, and elasticity modulus of FRP laminates. According to their study, different output parameters, including the service load and ultimate load, were used. Overall, their study showed that ANN is a useful technique for assessing the aforementioned problem. Another ANN-based study was performed by Metwally [37]. The output of his proposed model was the ultimate load. He used 14 input parameters to construct his intelligent model. Shanmugavelu et al. [38] also highlighted the workability of ANN in assessing the performance characteristics of reinforced concrete beams strengthened with glass FRP laminates.
Saadoon and Malik [39] developed an ANN-based model for predicting the ultimate load of a concrete beam that was strengthened with FRP bars. They used 199 sets of data collected from the literature for their model construction. The proposed model consisted of eight input parameters. The input parameters of their model consisted of the cross-sectional width and depth of the beams, cross-sectional area, elasticity modulus, tensile strength of FRP bars, compressive strength of concrete, effective span length, and shear-span ratio. Al-Jurmma [40] also showed that ANN is a feasible tool for predicting the ultimate load capacity of reinforced concrete beams. The database of their ANN-based predictive model was constructed using finite-element-based numerical analyses.
Erdem [41] suggested an ANN-based predictive model of moment capacity for reinforced concrete slabs. They used 294 sets of data for their model construction. Cai et al. [42] suggested an ANN-based predictive model that was enhanced with a genetic algorithm to assess the postfire flexural capacity of reinforced concrete beams. They used 480 sets of data, which were produced with the aid of finite element analyses. The input parameters of their proposed model comprised the height and width of the beams, cross-sectional area and tensile strength of reinforcement, compressive strength and cover of concrete, and fire time. The coefficient of determination value of 0.99 showed that a GA-based ANN was a quick and feasible tool for flexural capacity prediction.
This study was aimed toward proposing intelligent models that were optimized using particle swarm optimization (PSO) and imperialist competitive algorithm (ICA). The PSO-based ANN and ICA-based ANN models were developed for predicting the flexural strength of concrete beams that were made from recycled aggregate. Although the implementation of conventional ANNs for concrete beams is not new, it is well established that the conventional ANNs suffer from some disadvantages, such as getting trapped in local minima and a slow rate of learning. Optimization algorithms, such as PSO and ICA, can enhance the performance of ANNs and overcome their shortcomings [26]. The main difference between the presented study and previous works is the point that in this study, ANNs were coupled with ICA and PSO for the flexural assessment of RRC beams. Apart from the artificial intelligence components, this paper presents the flexural performance of two experimental tests. The results of the two loading tests on concrete beams were added to the compiled sets of data that were obtained from literature and used for the predictive models' construction. Presenting the results of experimental tests is always of interest as these results can be utilized for enriching the related databases for further soft computing works. Nevertheless, after reviewing the related works in this section, in Section 2 of the presented paper, the utilized artificial intelligence methods are discussed. Section 3 deals with the implemented experimental program. In Section 4, the PSO-based ANN and the ICA-based ANN modeling procedures are discussed. Section 5 deals with the main results and discussion, and finally, the last section deals with the summary of the paper and concluding remarks.

Methods
In this section, the proposed methodology for predicting the flexural strength of concrete beams is discussed. At first, some explanations about ANN, PSO, and ICA are provided, and finally, the hybridization of the mentioned methods is proposed.

Artificial Neural Network
The use of artificial neural networks (ANNs) has become so pervasive among researchers. Artificial neural networks can be used for any system that requires analysis, decision-making, estimation, forecasting, design, and construction. In all these models, a mathematical structure is considered. This general structure is tuned and optimized using a training algorithm so that it can behave accurately. In civil engineering, ANNs are mostly considered good function approximation tools, especially for nonlinear problems. One of the most basic neural models available is the multi-layer perceptron (MLP), which was described by Dreyfus [43]. In this type of neural network, which is also called a feed-forward neural network, the behavior of the human brain is mimicked. Feed-forward ANNs are mainly made of three layers: the input layer, hidden layer, and output layer. Each layer comprises one or more neurons (also called nodes). The main engine of the ANNs is the hidden layer, where the hidden nodes receive information from the input layer (the influential parameters on the model output). It is worth mentioning that the information is transferred between layers through connection weights, as well as threshold values, which are called bias values. The net input of each hidden node is the summation of input weights plus the bias value of the specific hidden node. The outputs of the hidden nodes are determined after applying transfer functions (often sigmoid functions) on the net inputs of the hidden nodes. This process continues until a reasonable output is achieved. It is worth mentioning that the connection weights have to be adjusted using training algorithms until the point that the error between the network outputs and target values is minimized. The aforementioned error is often assessed using the mean square error (MSE) of the system or root-mean-square error (RMSE) of the system.

Particle Swarm Optimization Algorithm
The PSO algorithm is one of the most significant intelligent optimization algorithms, which is based on swarm intelligence. The PSO algorithm was first developed by Kennedy and Eberhart [44]. The particle swarm optimization algorithm has a memory so that the knowledge of possible solutions is reserved for all particles. Each member of the swarm changes their position according to personal experiences and the experiences of the whole society. The social sharing of information between members of society has some advantages, and this hypothesis is the basis of the particle swarm optimization algorithm.
In the PSO method, after the initialization of the population, the particles usually converge in successive iterations of the local search process in a feasible space. PSO works with two simple equations. The position and the velocity of each particle are defined through Equations (1) and (2): where V k+1 i and X k i are the ith particle's velocity and position at k and k + 1 iterations, respectively; c 1 and c 2 are two acceleration constants; P k i is the best position of the ith particle up to iteration k; and P k g is the best position among all particles in the swarm up to iteration k. In Figure 1, the velocity updating process is shown.
The PSO algorithm is one of the most significant intelligent optimization algorithms, which is based on swarm intelligence. The PSO algorithm was first developed by Kennedy and Eberhart [44]. The particle swarm optimization algorithm has a memory so that the knowledge of possible solutions is reserved for all particles. Each member of the swarm changes their position according to personal experiences and the experiences of the whole society. The social sharing of information between members of society has some advantages, and this hypothesis is the basis of the particle swarm optimization algorithm.
In the PSO method, after the initialization of the population, the particles usually converge in successive iterations of the local search process in a feasible space. PSO works with two simple equations. The position and the velocity of each particle are defined through Equations (1) and (2): where 1 k i V  and k i X are the ith particle's velocity and position at k and k + 1 iterations, respectively; c1 and c2 are two acceleration constants; k i P is the best position of the ith particle up to iteration k; and k g P is the best position among all particles in the swarm up to iteration k. In Figure 1, the velocity updating process is shown.

Imperialist-Competitive-Algorithm-Based ANN
In order to study the effect of different optimization algorithms on the ANN performance and to enrich the soft computing part of the study, the imperialist competitive algorithm (ICA) was also utilized to improve the ANN performance. Details on the ICA are beyond the scope of this study and can be found elsewhere [45,46]. Nevertheless, the ICA as a global search algorithm starts with generating a user-defined number of so-called countries that play the role of particles in PSO. Like PSO, each country (particle) can be a potential solution to the problem. However, unlike PSO, in the ICA, the countries are categorized into imperialists and colonies. In fact, countries are sorted based on their fitness. The fitness function (sometimes it is called the cost function) in most engineering problems is the computed error between the predicted and measured values (often the mean square error (MSE). Nevertheless, the proper number of imperialists, which are the fittest countries (i.e., countries with the lowest MSE values) can be determined after a trial-anderror procedure. The remaining countries are so-called colonies.
A certain number of colonies are assigned to each imperialist (empire) based on the empire's power (normalized cost of each empire). Consequently, more powerful empires have more colonies. For a better understanding, the most powerful empire is the country with the lowest MSE value. In fact, designers set an initial number of countries and imperialists; subsequently, they investigate the model performance. Like PSO, an iteration

Imperialist-Competitive-Algorithm-Based ANN
In order to study the effect of different optimization algorithms on the ANN performance and to enrich the soft computing part of the study, the imperialist competitive algorithm (ICA) was also utilized to improve the ANN performance. Details on the ICA are beyond the scope of this study and can be found elsewhere [45,46]. Nevertheless, the ICA as a global search algorithm starts with generating a user-defined number of so-called countries that play the role of particles in PSO. Like PSO, each country (particle) can be a potential solution to the problem. However, unlike PSO, in the ICA, the countries are categorized into imperialists and colonies. In fact, countries are sorted based on their fitness. The fitness function (sometimes it is called the cost function) in most engineering problems is the computed error between the predicted and measured values (often the mean square error (MSE). Nevertheless, the proper number of imperialists, which are the fittest countries (i.e., countries with the lowest MSE values) can be determined after a trial-and-error procedure. The remaining countries are so-called colonies.
A certain number of colonies are assigned to each imperialist (empire) based on the empire's power (normalized cost of each empire). Consequently, more powerful empires have more colonies. For a better understanding, the most powerful empire is the country with the lowest MSE value. In fact, designers set an initial number of countries and imperialists; subsequently, they investigate the model performance. Like PSO, an iteration process (in the ICA, it is called a decade) is required for optimization. The proper number of decades can also be obtained using a trial-and-error procedure.
In essence, the ICA works with its three operators: assimilation, revolution, and competition. Colonies will be attracted to the empires using assimilation. Revolutions can suddenly alter the positions of the countries. It is worth mentioning that during assimilation and revolution, a colony may reach a state that is better than its empire state. In competition, imperialists try to obtain more colonies, and all empires attempt to take possession of the colonies of other imperialists. All the imperialists have the chance to take control of at least one colony of the weakest empire (based on the imperialists' power). Hence, during competition, the weak empires start to collapse, and on the other hand, the more powerful imperialists gain more power and expand their empires. This process often continues until a designer-defined termination criterion such as the maximum number of decades is met. At last, the most powerful empire (the country with the lowest MSE value) represents the solution to the problem of interest.

Hybrid ANNs (PSO-Based ANN and ICA-Based ANN)
The artificial neural network is often utilized for prediction and cannot perform search operations. In fact, the process of training in the ANN requires weight optimization. In supervised ANNs, the network receives a training algorithm and calculates the output using the weights in the randomly assigned network. However, many researchers highlighted that the slow rate of learning is one of the major drawbacks of ANNs (e.g., Momeni et al. [47]). The neural network can help to some extent but does not guarantee that the optimal answer is fully shown and converged toward. In fact, it is well established that getting trapped in local minima is another drawback of ANNs [48].
The idea of improving the ANN performance with the aid of global search algorithms, such as PSO, was proposed in previous studies for other applications, mainly due to the point that the neural network results are not necessarily the global optimal solutions [49]. In fact, PSO can be implemented for the weight optimization of ANNs. The PSO algorithm's initial population comprises the possible ANN weights (each particle contains the whole ANN weights). The PSO algorithm, which has a high power of global exploration, can easily obtain the optimal global answer by sharing information between particles. In order to better understand the hybrid process of PSO-based ANN, the optimization process is shown in the following flowchart (see Figure 2). It is worth mentioning that the hybrid process of an ICA-based ANN is similar to the PSO-based ANN, except that in an ICA-based ANN, the ICA algorithm is implemented to optimize the ANN weights.

Experimental-Based Dataset
A comprehensive database plays a crucial role in designing artificial-intelligencebased predictive models. It is well established that the size of the database can enhance the efficiency of the predictive models. Nevertheless, in civil engineering problems where performing numerous experiments is a difficult task to be accomplished, compiling different sets of data from literature is common, especially when the scope of a study is on highlighting the feasibility of soft computing techniques. Therefore, this study used mostly used the previously published data in the literature. However, in this study, to enrich the quality of the work, apart from compiled data from literature [1,5,[9][10][11][13][14][15][16]50], the results of two experimental tests performed in the Structural Laboratory of Lorestan University were added to the aforementioned database. The experimental procedure for the tests is highlighted in the next section. Overall, 107 sets of data were used in this study. It is worth mentioning that in the context of AI-based predictive models, setting several input parameters, as well as an output parameter(s), is essential. Any parameter can be set as an input for the predictive model if there is a meaningful relationship between the input and output parameters [51]. In the case of necessity, a sensitivity anal-

Experimental-Based Dataset
A comprehensive database plays a crucial role in designing AI-based predictive models. It is well established that the size of database can enhance the efficiency of the predictive models. Nevertheless, in civil engineering problems where performing numerous experiments is a difficult task to be accomplished, compiling different sets of data from literature is common, especially when the scope of a study is on highlighting the feasibility of soft computing techniques. Therefore, this study mostly used the previously published data in the literature.
However, in this study, to enrich the quality of the work, apart from compiled data from literature [1,5,[9][10][11][13][14][15][16]50], the results of two experimental tests performed in the Structural Laboratory of Lorestan University were added to the aforementioned database. The experimental procedure for the tests is highlighted in the next section. Overall, 107 sets of data were used in this study. It is worth mentioning that in the context of AI-based predictive models, setting several input parameters, as well as an output parameter(s), is essential. Any parameter can be set as an input for the predictive model if there is a meaningful relationship between the input and output parameters [51]. In the case of necessity, a sensitivity analysis can be implemented to reduce the number of influential parameters.
The input parameters used in this study comprised the compressive strength of concrete, longitudinal reinforcement ratio, the width of the beam, effective depth of the section, amount of RCA in terms of the percentage, span's-length-to-effective-depth ratio, yield strength of steel, and ratio of the distance of the applied load from the beam edge to the effective depth of the beam. The output parameter was the flexural strength of the beam. Table 1 shows the summary of the implemented dataset in this study.

Experimental Procedure of the Performed Tests
In order to investigate the flexural performances of the reinforced beams, first, the raw materials were prepared (cement, natural and recycled aggregates). It should be mentioned that in this study, ordinary Portland cement (ASTM type II) was used according to ASTM C150 [52]. The chemical properties of the utilized cement are shown in Table 2. The natural aggregates with a maximum size of 10 mm were obtained from the Khorramabad River. On the other hand, the preparation of the recycled concrete aggregate had several stages. The first stage dealt with crushing the old concrete specimens by using a crusher machine. Subsequently, these materials were washed and granulated into coarse and fine aggregate. In this research, to achieve a high-quality recycled aggregate (RA), four hundred cubic old concrete specimens with dimensions of 150 × 150 × 150 mm were chosen. These cubic specimens (see Figure 3) were collected from Structure Laboratory at Lorestan University and were crushed in accordance with ASTM D8038 [53] and ACI 555 [54]. In this study, the coarse RA had a maximum size of 9 mm according to ASTM C125 [55], and it had a fractured percentage of more than 87%, conforming to ASTM D5821 [56]. it had a fractured percentage of more than 87%, conforming to ASTM D5821 [56].  In the next stage, to investigate the workability of the RCA, as well as provide dataset enrichment, two 1500 mm reinforced concrete (RC) beams with a width of 200 mm and a height of 300 mm were constructed and tested. These specimens were subjected to a concentrated load in the middle of their span, and the loading process was continued until the failure occurred. The geometrical properties of these concrete beams are presented in Figure 4. As can be seen from Figure 4, three T14 steel bars were used as a longitudinal tension bar, two T10 steel bars were utilized as a compression bar, and ten T10 steel bars In the next stage, to investigate the workability of the RCA, two 1500 mm reinforced concrete (RC) beams with a width of 200 mm and a height of 300 mm were constructed and tested. These specimens were subjected to a concentrated load in the middle of their span, and the loading process was continued until the failure occurred. The geometrical properties of these concrete beams are presented in Figure 4. As can be seen from Figure 4, three T14 steel bars were used as a longitudinal tension bar, two T10 steel bars were utilized as a compression bar, and ten T10 steel bars were used as a stirrup for each specimen. The tensile stress test was performed on the steel bars in accordance with ASTM C370 [57]. were used as a stirrup for each specimen. The tensile stress test was performed on the steel bars in accordance with ASTM C370 [57]. The longitudinal bars with diameters of 14 and 10 mm had yielding stresses of 293 and 369 MPa, respectively. Furthermore, these bars had failure stress of 520 and 562 MPa. Based on the tensile test results, the yielding and failure stresses of the used stirrups were 281 and 508 MPa, respectively. For the purpose of brevity, the aforementioned results are not discussed here. In this research, to investigate the effects of RAs on the flexural behavior of reinforced concrete beams, the two abovementioned beams were tested under pure bending. The specimens had the same rebars and dimensions, but in terms of the aggregate, they were made of natural aggregate and RAs. It is worth mentioning that the aggregates used in both beams had the same gradation. Figure 5 shows the distribution curves of the recycled aggregates.  For the concrete construction, first, fine and coarse aggregates were added to the The longitudinal bars with diameters of 14 and 10 mm had yielding stresses of 293 and 369 MPa, respectively. Furthermore, these bars had failure stress of 520 and 562 MPa. Based on the tensile test results, the yielding and failure stresses of the used stirrups were 281 and 508 MPa, respectively. For the purpose of brevity, the aforementioned results are not discussed here. In this research, to investigate the effects of RAs on the flexural behavior of reinforced concrete beams, the two abovementioned beams were tested under pure bending. The specimens had the same rebars and dimensions, but in terms of the aggregate, they were made of natural aggregate and RAs. It is worth mentioning that the aggregates used in both beams had the same gradation. Figure 5 shows the distribution curves of the recycled aggregates.  For the concrete construction, first, fine and coarse aggregates were added to the mixer and after three minutes, the cement and water were added to the aggregates and mixed for two more minutes. After making the concrete, it was poured into the beam mold, and after 28 days of curing the specimens, they were tested. As shown schematically in Figure 6, the beams were subjected to a net bending test in the middle of the span. Three LVDTs were used to measure the displacement of the specimens. Moreover, using a load cell with a capacity of 500 kN, the applied force was measured. The load applied to the beams was in the form of displacement control. The loading process continued until the specimens collapsed. Figure 7 shows the mid-span displacement versus the applied force. For the concrete construction, first, fine and coarse aggregates were added to the mixer and after three minutes, the cement and water were added to the aggregates and mixed for two more minutes. After making the concrete, it was poured into the beam mold, and after 28 days of curing the specimens, they were tested. As shown schematically in Figure 6, the beams were subjected to a net bending test in the middle of the span. Three LVDTs were used to measure the displacement of the specimens. Moreover, using a load cell with a capacity of 500 kN, the applied force was measured. The load applied to the beams was in the form of displacement control. The loading process continued until the specimens collapsed. Figure 7 shows the mid-span displacement versus the applied force. The maximum load capacity of the RRC beam was approximately 10% less than the reference beam (RC beam) (including natural aggregate). Furthermore, this behavior was acceptable in accordance with the RA features. As we know, flexural behavior in concrete beams depends on the compressive and tension zones of specimens. According to the results, having a rough surface had a profound effect on flexural behavior. Additionally, the natural aggregates had rougher surfaces in comparison with the RAs; thus, this feature had a large effect on the aggregate-cement bond strength. The concrete ultimate strain rose with an increase in the interlock between the cement paste and aggregate's surface. In fact, increasing the adhesion of cement affected the ultimate flexural strength of the beam specimens. The maximum load capacity of the RRC beam was approximately 10% less than the reference beam (RC beam) (including natural aggregate). Furthermore, this behavior was acceptable in accordance with the RA features. As we know, flexural behavior in concrete beams depends on the compressive and tension zones of specimens. According to the results, having a rough surface had a profound effect on flexural behavior. Additionally, the natural aggregates had rougher surfaces in comparison with the RAs; thus, this feature had a large effect on the aggregate-cement bond strength. The concrete ultimate strain rose with an increase in the interlock between the cement paste and aggregate's surface. In fact, increasing the adhesion of cement affected the ultimate flexural strength of the beam specimens.   Based on the experimental results, the initial slopes of the load-displacement curve for the two beams were approximately similar. When the applied load reached 40% of the ultimate load capacity, the linear behavior changed to nonlinear behavior. The flexural behaviors of the beams showed that with an increase in the load up to 40% of its ultimate capacity, the longitudinal rebars could reach their yielding strength. Moreover, the yielding process was accompanied by crack propagation on the beam's surface.
Concerning the results, tri-linear behavior could be assumed for beam specimens, and this behavior was classified into three steps. The occurrence of the first cracks in the mid-span was identified as the first stage. The yielding and collapse processes were recognized as the second and third stages in this study, respectively. The failure mode of the RRC beam was different compared with another specimen. The compressive zone failure dimension, crack width, shear-bending cracks, and crushing nature of the concrete were seen as the main differences. In this regard, Figure 8 shows the crack propagation of the RRC beam and the conventional reinforced concrete (RC) beam after failure. Based on the experimental results, the initial slope of the load-displacement curve for the two beams was approximately similar. When the applied load reached 40% of the ultimate load capacity, the linear behavior changed to nonlinear behavior. The flexural behaviors of the beams showed that with an increase in the load up to 40% of its ultimate capacity, the longitudinal rebars could reach their yielding strength. Moreover, the yielding process was accompanied by crack propagation on the beam's surface.
Concerning the results, tri-linear behavior could be assumed for beam specimens, and this behavior was classified into three steps. The occurrence of the first cracks in the mid-span was identified as the first stage. The yielding and collapse processes were recognized as the second and third stages in this study, respectively. The failure mode of the RRC beam was different compared with another specimen. The compressive zone failure dimension, crack width, shear-bending cracks, and crushing nature of the concrete were seen as the main differences. In this regard, Figure 8 shows the crack propagation of the RRC beam and the conventional reinforced concrete (RC) beam after failure.
Overall, a comparison of the load capacity curves for the two beams showed that the corresponding displacement of the failure point for the RRC beam was almost 10% less than that of the RC beam. It is worth mentioning that for the purpose of brevity and owing to the nature of this work, the experimental results are not discussed in detail. As mentioned earlier, in PSO-based ANN modeling, the ANN weights are imported into the PSO algorithm, and after optimization, the optimized weights are exported into Overall, a comparison of the load capacity curves for the two beams showed that the corresponding displacement of the failure point for the RRC beam was almost 10% less than that of the RC beam. It is worth mentioning that for the purpose of brevity and owing to the nature of this work, the experimental results are not discussed in detail.

PSO-Based ANN Modeling Procedure
As mentioned earlier, in PSO-based ANN modeling, the ANN weights are imported into the PSO algorithm, and after optimization, the optimized weights are exported into ANN for prediction purposes. Hence, before the ANN structure configuration, sensitivity analyses should be performed on the parameters of the PSO algorithm. The first usual sensitivity analysis is often implemented to determine the optimal number of iterations. For this reason, the number of iterations was set to 100 for a default PSO-based ANN model with nine hidden nodes. It is worth mentioning that the PSO coefficients, i.e., C 1 and C 2 , were set as 2 for the primary sensitivity analysis. The primary number of particles was set to 200. Figure 9 shows the influence of the number of iterations on the performance of the PSO-based ANN predictive model. This figure suggests that after 40 iterations, there was a negligible change in MSE; hence, the number of iterations used during the sensitivity analyses was set to 40. After determining the optimal number of iterations, the optimal number of particles, as well as the proper values for C1 and C2, should be identified. For this reason, the effect of the number of particles on the model performance was investigated. For this reason, several models with 100, 150, 200, 300, 400, and 500 particles were run and it was found that the number of particles had no remarkable effect on the model performance, which was in good agreement with the authors' previous works [47]. On the other hand, to investigate the importance of PSO coefficients, the effect of considering different values for C1 and C2 was taken into consideration. For the purpose of brevity, details on the sensitivity results are not presented here; however, it was found that setting the values of C1 and C2 equal to 2 can lead to an acceptable prediction performance; hence, the aforementioned values were set to 2, as suggested in the literature [48]. It is worth mentioning that when the results of the sensitivity analyses were close to each other, the performance of the testing data was considered. It should be underlined that in all the analyses, 80% of the data were considered for training the predictive model and the remaining data were used to test the developed model. After identifying the proper values for PSO parameters, the best ANN structure should be defined. In other words, the optimal number of hidden nodes should be identified. Researchers often implement two approaches to determine the optimal number of hidden nodes: using suggested equations and/or a trial-and-error proce- After determining the optimal number of iterations, the optimal number of particles, as well as the proper values for C 1 and C 2 , should be identified. For this reason, the effect of the number of particles on the model performance was investigated. For this reason, several models with 100, 150, 200, 300, 400, and 500 particles were run and it was found that the number of particles had no remarkable effect on the model performance, which was in good agreement with the authors' previous works [47]. On the other hand, to investigate the importance of PSO coefficients, the effect of considering different values for C 1 and C 2 was taken into consideration. For the purpose of brevity, details on the sensitivity results are not presented here; however, it was found that setting the values of C 1 and C 2 equal to 2 can lead to an acceptable prediction performance; hence, the aforementioned values were set to 2, as suggested in the literature [48]. It is worth mentioning that when the results of the sensitivity analyses were close to each other, the performance of the testing data was considered. It should be underlined that in all the analyses, 80% of the data were considered for training the predictive model and the remaining data were used to test the developed model. After identifying the proper values for PSO parameters, the best ANN structure should be defined. In other words, the optimal number of hidden nodes should be identified. Researchers often implement two approaches to determine the optimal number of hidden nodes: using suggested equations and/or a trial-and-error procedure. In this study, a trial-and-error procedure was implemented to determine the optimal number of hidden nodes. The effects of 7, 8, 9, 10, 11, and 12 hidden nodes on the model performance was investigated. It was found that when nine hidden nodes were used in the PSO-based ANN predictive model, the model performed the best. For the purpose of brevity, only the results of the PSO-based ANN predictive models with 8 and 9 nodes are presented in Figures 10 and 11, respectively. It is worth mentioning that to reduce the likelihood of accidental results, each model was run five times. Nevertheless, as suggested in Figure 11, the second predictive model coefficient correlation of 0.997 and MSE of 0.08% for testing data outperformed other models. The overall results of different runs with nine hidden nodes showed that the PSO-based ANN predictive model worked well enough, as shown in Figures 10 and 11. The nine-node predictive model outperformed the eight-node predictive model, although the results of the latter were also reliable to some extent. The best results are illustrated in a better way in the Main Results and Discussion section.

ICA-Based ANN Modeling Procedure
The first step in developing the ICA-based ANN model dealt with identifying the optimal number of decades. For this reason, a default model had to be considered. The considered default model comprised 200 countries with nine hidden nodes in one hidden layer. A total of 10% of countries were set to be imperialists. The initial number of decades was set to be 100. Similar to the PSO-based ANN model, 80% of the data was assigned to train the model and the remaining data were used to test the prediction performance of the model. Figure 12 shows the variation of MSE values versus the number of decades. This figure suggests that after 36 decades, the change in MSE was negligible; hence, the optimal number of decades was set to 36 for further analysis. Identifying the optimal number of countries and imperialists was the next performed sensitivity analysis for developing the ICA-based ANN model. For this reason, the number of countries varied from 125 to 400, as shown in Table 3. This table also showed the number of imperialists' variations. Similar to the previous section, the utilized performance indices were the R and MSE values. When the prediction performances of different models were close to each other, the prediction performance of the testing data was preferred.

ICA-Based ANN Modeling Procedure
The first step in developing the ICA-based ANN model dealt with identifying the optimal number of decades. For this reason, a default model had to be considered. The considered default model comprised 200 countries with nine hidden nodes in one hidden layer. A total of 10% of countries were set to be imperialists. The initial number of decades was set to be 100. Similar to the PSO-based ANN model, 80% of the data was assigned to train the model and the remaining data were used to test the prediction performance of the model. Figure 12 shows the variation of MSE values versus the number of decades. This figure suggests that after 36 decades, the change in MSE was negligible; hence, the optimal number of decades was set to 36 for further analysis. Identifying the optimal number of countries and imperialists was the next performed sensitivity analysis for developing the ICA-based ANN model. For this reason, the number of countries varied from 125 to 400, as shown in Table 3. This table also showed the number of imperialists' variations. Similar to the previous section, the utilized performance indices were the R and MSE values. When the prediction performances of different models were close to each other, the prediction performance of the testing data was preferred.  Table 3 also shows the results of the aforementioned sensitivity analysis. The results presented in this table suggested that the best prediction performance of the ICA-based ANN predictive model was good enough when the numbers of countries and imperialists were 200 and 20, respectively. Hence, during the model development process, these default values were considered as the optimal number of countries and imperialists, respectively. Table 3. The effect of ICA parameters on the prediction performances of the models.

Model
Parameter Training Data Testing Data  Table 3 also shows the results of the aforementioned sensitivity analysis. The results presented in this table suggested that the best prediction performance of the ICA-based ANN predictive model was good enough when the numbers of countries and imperialists were 200 and 20, respectively. Hence, during the model development process, these default values were considered as the optimal number of countries and imperialists, respectively. At last, to obtain a better comparison, the ICA-based ANN model with nine hidden nodes in one hidden layer was iterated five times (similar to the PSO-based ANN model to reduce the likelihood of accidental results). The prediction performances of the nine-node ICA-based ANN predictive models are shown in Figure 13. As shown in this figure, the R and MSE values (0.994 and 0.023) of the testing data suggested that the second model performed the best. At last, to obtain a better comparison, the ICA-based ANN model with nine hidden nodes in one hidden layer was iterated five times (similar to the PSO-based ANN model to reduce the likelihood of accidental results). The prediction performances of the ninenode ICA-based ANN predictive models are shown in Figure 13. As shown in this figure, the R and MSE values (0.994 and 0.023) of the testing data suggested that the second model performed the best. Figure 13. Prediction performances of the 9-node ICA-based ANN predictive model.  Figure 14 shows the prediction performance of the PSO-based ANN predictive model of the flexural strength of the RRC beams. The correlation coefficient of 0.994 for the training data showed that the predictive model worked well enough. On the other hand, Figure 15 displays the predicted flexural strengths of the RRC beams versus the measured values for the testing dataset. As shown in this figure, the correlation coefficient of 0.997 suggested that the PSO-based ANN model was a capable and feasible tool for assessing the flexural behavior of RRC beams. It is worth mentioning that MSE values for the training and testing data were 0.06% and 0.08%, respectively.  Similarly, Figures 16 and 17 display the prediction performance of the ICA-based ANN predictive model for the training and testing data, respectively. The R values of 0.967 and 0.994 for the training and testing data, respectively, showed the workability of the ICA-based ANN model for assessing the flexural strength of RRC beams. Despite the relatively good prediction performance of the ICA-based ANN model, the results showed that the PSO-based ANN model outperformed the ICA-based ANN model. To support the aforementioned conclusion, the prediction performances of the ICA-based ANN and PSO-based ANN models for testing data were assessed using another performance index known as variance accounted for (VAF). VAF can be computed using the following equa-  Similarly, Figures 16 and 17 display the prediction performance of the ICA-based ANN predictive model for the training and testing data, respectively. The R values of 0.967 and 0.994 for the training and testing data, respectively, showed the workability of the ICA-based ANN model for assessing the flexural strength of RRC beams. Despite the relatively good prediction performance of the ICA-based ANN model, the results showed that the PSO-based ANN model outperformed the ICA-based ANN model. To support the aforementioned conclusion, the prediction performances of the ICA-based ANN and PSO-based ANN models for testing data were assessed using another performance index known as variance accounted for (VAF). VAF can be computed using the following equation: Similarly, Figures 16 and 17 display the prediction performance of the ICA-based ANN predictive model for the training and testing data, respectively. The R values of 0.967 and 0.994 for the training and testing data, respectively, showed the workability of the ICA-based ANN model for assessing the flexural strength of RRC beams. Despite the relatively good prediction performance of the ICA-based ANN model, the results showed that the PSO-based ANN model outperformed the ICA-based ANN model. To support the aforementioned conclusion, the prediction performances of the ICA-based ANN and PSO-based ANN models for testing data were assessed using another performance index known as variance accounted for (VAF). VAF can be computed using the following equation:

Main Results and Discussion
where α and α' are measured and predicted values, respectively, and var indicates the variance.    At last, for comparison purposes, the prediction performance of the PSO-based A predictive model was checked against a conventional ANN model (famous backprop tion neural network) with the same structure (i.e., nine hidden nodes). The obtained sults for the conventional ANN (i.e., R-value of 0.884 for testing data) showed that PSO-based ANN model also outperformed the conventional ANN-based predic model in assessing the flexural behavior of the RRC beams.     At last, for comparison purposes, the prediction performance of the PSO-based AN predictive model was checked against a conventional ANN model (famous backpropa tion neural network) with the same structure (i.e., nine hidden nodes). The obtained sults for the conventional ANN (i.e., R-value of 0.884 for testing data) showed that PSO-based ANN model also outperformed the conventional ANN-based predict model in assessing the flexural behavior of the RRC beams.    Figure 18. Predicted values using the recommended predictive models versus the measured values (testing data).

Summary and Conclusions
In this study, two flexural tests were performed on concrete beams. The length, width, and depth of reinforced concrete beams were 1500 mm, 200 mm, and 300 mm, respectively. The results of the flexural tests showed an almost 10% reduction in ultimate flexural strength when recycled aggregates were used instead of natural aggregates. However, using the experimental results of this study, as well as the compiled experimental data from the literature, a relatively comprehensive database with 107 sets of data was created for the soft computing part of the study. The soft computing part of the study comprised the implementation of ANNs that were improved with PSO, as well as ICA for the flexural assessment of RRC beams. Additionally, for comparison purposes, the PSObased ANN model was compared with the ICA-based ANN model and a conventional ANN, and it was concluded that the PSO-based ANN model worked better compared with the aforementioned models.
It is worth mentioning that the inputs of the predictive models were RCA%, B, d, a/d, L/d, ρ, fc, and Fy. Overall, the R-value of 0.997, VAF value of 99.35, and MSE value of 0.08% for the testing data showed that the PSO-based ANN predictive model with nine hidden nodes in one hidden layer was a feasible tool for predicting the ultimate flexural strength of the RRC beams. Despite the promising results, a word of caution is required regarding generalizing the prediction performance of every predictive model. In fact, the reliability of soft-computing-based predictive models depends on the quality, quantity, and range of the training data. Hence, when the range of future data is beyond the range of the implemented dataset, the predictability of intelligent models is open to question. In this regard, further research on enhancing the implemented dataset in this study is recommended.   At last, for comparison purposes, the prediction performance of the PSO-based ANN predictive model was checked against a conventional ANN model (famous backpropagation neural network) with the same structure (i.e., nine hidden nodes). The obtained results for the conventional ANN (i.e., R-value of 0.884 for testing data) showed that the PSO-based ANN model also outperformed the conventional ANN-based predictive model in assessing the flexural behavior of the RRC beams.

Summary and Conclusions
In this study, two flexural tests were performed on concrete beams. The length, width, and depth of reinforced concrete beams were 1500 mm, 200 mm, and 300 mm, respectively. The results of the flexural tests showed an almost 10% reduction in ultimate flexural strength when recycled aggregates were used instead of natural aggregates. However, using the experimental results of this study, as well as the compiled experimental data from the literature, a relatively comprehensive database with 107 sets of data was created for the soft computing part of the study. The soft computing part of the study comprised the implementation of ANNs that were improved with PSO, as well as ICA for the flexural assessment of RRC beams. Additionally, for comparison purposes, the PSO-based ANN model was compared with the ICA-based ANN model and a conventional ANN, and it was concluded that the PSO-based ANN model worked better compared with the aforementioned models.
It is worth mentioning that the inputs of the predictive models were RCA%, B, d, a/d, L/d, ρ, f c , and F y . Overall, the R-value of 0.997, VAF value of 99.35, and MSE value of 0.08% for the testing data showed that the PSO-based ANN predictive model with nine hidden nodes in one hidden layer was a feasible tool for predicting the ultimate flexural strength of the RRC beams. Despite the promising results, a word of caution is required regarding generalizing the prediction performance of every predictive model. In fact, the reliability of soft computing-based predictive models depends on the quality, quantity, and range of the training data. Hence, when the range of future data is beyond the range of the implemented dataset, the predictability of intelligent models is open to question. In this regard, further research on enhancing the implemented dataset in this study is recommended.