Acessibilidade / Reportar erro

Application of Support Vector Machine and Finite Element Method to predict the mechanical properties of concrete

Abstract

Compressive strength and Young’s modulus are the main properties used in the design of concrete structures. They are responsible for the cost, safety and dimensioning of the structure, and are generally measured in expensive and time-demanding tests. This fact encourages researches for fast and cost-effective methods to investigate the concrete’s properties. Among the concrete types, structural Lightweight Aggregate Concrete (LWAC) is one of the most employed worldwide, but it presents limited studies and mix design techniques. Thus, this work evaluates and compares the performances of two methods to predict the compressive strength of LWAC samples: Support Vector Machine and Finite Element Method. To this end, both strategies use the LWAC’s mix proportions and the Young’s modulus, and the compressive strength of mortars and aggregates obtained from an experimental program from the literature. The results encourage further researches towards the development of a numerical tool that may assist engineers for practical purposes, since both methods show good agreement with the validation data.

Keywords:
Lightweight Aggregate Concrete; Finite Element Method; Computational Intelligence; Support Vector Regression

1. Introduction

Concrete is one of the most widely employed materials in the world. This composite comprises a matrix (the mortar - cement, fine aggregates and water) and a particulate solid phase (the coarse aggregates). When the concrete reaches a compressive strength (f c ) higher than 20 MPa, it is considered structural concrete (ABNT, 2015ABNT, 2015. NBR 8953 - Concrete for structural use - Density, strength and consistence classification, Rio de Janeiro: Associação Brasileira de Normas Técnicas.). They can be classified according to their dry unit weight (ρ) in lightweight (ρ < 2000 kg/m3), normal weight (2000 kg/m3 < ρ < 2800 kg/m3) and heavyweight (ρ > 2800 kg/m3) (ABNT, 2015ABNT, 2015. NBR 8953 - Concrete for structural use - Density, strength and consistence classification, Rio de Janeiro: Associação Brasileira de Normas Técnicas.). Lightweight concrete can be achieved through foaming-inducing admixtures, air-entraining admixtures or using lightweight aggregates (LWA). In turn, lightweight aggregates (or low-density aggregates) are those with bulk density less than 1120 kg/m3; e.g. pumice, volcanic cinders, tuff, expanded or sintered clay, shale, slate and vermiculite (ASTM, 2018ASTM, 2018. ASTM C125 - Standard Terminology Relating to Concrete and Concrete Aggregates, West Conshohocken: ASTM International.).

Nowadays structural lightweight aggregate concrete (LWAC) is used on a large scale to reduce the self-weight of structures. In addition, it also presents transport advantages, allows longer spams, has thermoacoustic insulation properties, presents reduced autogenous shrinkage and a negligible interfacial transition zone due to internal curing (Mehta & Monteiro, 2014Mehta, P. & Monteiro, P. J. M., 2014. Concrete: Microstructure, Properties, and Materials. 2ª ed. São Paulo: IBRACON.) (Bentur et al., 2001Bentur, A., Igarashi, S. I. & Kovler, K., 2001. Prevention of autogenous shrinkage in high-strength concrete by internal curing using wet lightweight aggregates. Cement and concrete research, 31(11), pp. 1587-1591.). Differently from conventional concretes, in which the aggregate has a higher compressive strength than the mortar and the interfacial transition zone is the weak area of the composite; in LWAC, the aggregate strength is the main factor controlling the fracture of concrete (Zhang & Gjvorv, 1991Zhang, M. H. & Gjvorv, O. E., 1991. Mechanical properties of high-strength lightweight concrete. ACI Materials Journal, 88(3), pp. 240-247.) (Yang & Huang, 1998Yang, C. C. & Huang, R., 1998. Approximate strength of lightweight aggregate using micromechanics method. Advanced Cement Based Materials, 7(3-4), pp. 133-138.).

The compressive strength and static Young’s Modulus of a concrete are amongst its most important mechanical properties, affecting the design, cost and safety of structures. These parameters are a direct function of cement type, water/cement ratio, aggregate type and the proportion of these ingredients, also called mix design. The mechanical properties are generally obtained by cost-demanding testing of concrete specimens after a standard curing period of 28 days, which means that 1) a significant amount of material must be spent for preparing the specimens and 2) only 28 days after poured into the structure they can be evaluated. This scenario is remarkably more sensitive for LWAC, from which there is not as many research projects and mix design techniques.

In this sense, several authors have been studying solutions to evaluate the concrete properties prior to its pouring, in order to avoid structural issues and material waste. Among these solutions, artificial intelligence algorithms have been used since the 2000s, especially the ones based on artificial neural networks (ANN) or support vector machines (SVM). For instance, Tenza-Abril et al. (2018Tenza-Abril, A. J., Villacampa, Y., Solak, A. M. & Baeza-Brotons, F., 2018. Prediction and sensitivity analysis of compressive strength in segregated lightweight concrete based on artificial neural network using ultrasonic pulse velocity. Construction and Building Materials, Volume 189, pp. 1173-1183.), Sadrmomtazi et al. (2013Sadrmomtazi, A., Sobhani, J. & Mirgozar, M. A., 2013. Modeling compressive strength of EPS lightweight concrete using regression, neural network and ANFIS. Construction and Building Materials, Volume 42, pp. 205-216.), Gupta, Kewalramani e Goel (2006Gupta, R., Kewalramani, M. A. & Goel, A., 2006. Prediction of concrete strength using neural-expert system. Journal of materials in civil engineering, 18(3), pp. 462-466.), Kim et al. (2005Kim, D. K., Lee, J. J., Lee, J. H. & Chang, S. K., 2005. Application of probabilistic neural networks for prediction of concrete strength. Journal of Materials in Civil Engineering, 17(3), pp. 353-362.), Lee (2003Lee, S. C., 2003. Prediction of concrete strength using artificial neural networks. Engineering Structures, 25(7), pp. 849-857.), and Ni e Wang (2000Ni, H. G. & Wang, J. Z., 2000. Prediction of compressive strength of concrete by neural networks. Cement and Concrete Research, 30(8), pp. 1245-1250.) employed artificial neural networks with satisfactory performance; while Abd & Abd (2017Abd, A. M. & Abd, S. M., 2017. Modelling the strength of lightweight foamed concrete using support vector machine (SVM). Case studies in construction materials, Volume 6, pp. 8-15.), Akande et al. (2014Akande, K. O., Owolabi, T. O., Twaha, S. & Olatunji, S. O., 2014. Performance comparison of SVM and ANN in predicting compressive strength of concrete. IOSR Journal of Computer Engineering, 16(5), pp. 88-94.), Yan e Shi (2010Yan, K. & Shi, C., 2010. Prediction of elastic modulus of normal and high strength concrete by support vector machine. Construction and building materials, 24(8), pp. 1479-1485.), among others, studied support vector machine (SVM), reaching promising results. The results obtained from ANN or SVM have, in general, presented similar performances; and their reliability is proportional to the quantity and quality of training examples available.

Non-conventional concretes such as high-strength concrete (Cheng et al., 2012Cheng, M. Y., Chou, J. S., Roy, A. F. & Wu, Y. W., 2012. High-performance concrete compressive strength prediction using time-weighted evolutionary fuzzy support vector machines inference model. Automation in Construction, Volume 28, pp. 106-115.) (Gupta, 2007Gupta, S. M., 2007. Support vector machines based modelling of concrete strength. World Academy of Science, Engineering and Technology, Volume 36, pp. 305-311.), lightweight foaming concrete (Abd & Abd, 2017Abd, A. M. & Abd, S. M., 2017. Modelling the strength of lightweight foamed concrete using support vector machine (SVM). Case studies in construction materials, Volume 6, pp. 8-15.) and even geopolymer concrete (Nazari & Sanjayan, 2015Nazari, A. & Sanjayan, J. G., 2015. Modelling of compressive strength of geopolymer paste, mortar and concrete by optimized support vector machine. Ceramics International, 41(9), pp. 12164-12177.) have also been evaluated via SVM, but papers concerning LWAC are rather scarce.

Also, seeking to simulate the concrete behavior in complex civil engineering structures, the Finite Element Method (FEM) is extensively adopted. Fracture and damage are some of the main phenomena studied in the literature (Yu et al., 2010Yu, T. T. J. G., Teng, J. G., Wong, Y. L. & Dong, S. L., 2010. Finite element modeling of confined concrete-I: Drucker-Prager type plasticity model. Engineering Structures, 32(3), pp. 665-679.) (Al-Rub & Kim, 2010Al-Rub, R. K. A. & Kim, S. M., 2010. Computational applications of a coupled plasticity-damage constitutive model for simulating plain concrete fracture. Engineering Fracture Mechanics, 77(10), pp. 1577-1603.) (Schlangen & Van Mier, 1992Schlangen, E. & Van Mier, J. G. M., 1992. Simple lattice model for numerical simulation of fracture of concrete materials and structures. Materials and Structures, 25(9), pp. 534-542.), along with the application of these models for the investigation of reinforced concrete (Genikomsou & Polak, 2015Genikomsou, A. S. & Polak, M. A., 2015. Finite element analysis of punching shear of concrete slabs using damaged plasticity model in ABAQUS. Engineering Structures, Volume 98, pp. 38-48.) (Hassan et al., 2015Hassan, N. Z., Sherif, A. G. & Zamarawy, A. H., 2015. Finite element analysis of reinforced concrete beams with opening strengthened using FRP. Ain Shams Engineering Journal.), concrete under fire conditions (Gao et al., 2013Gao, W. Y., Dai, J. G. T. J. G. & Chen, G. M., 2013. Finite element modeling of reinforced concrete beams exposed to fire. Engineering structures, Volume 52, pp. 488-501.) (Ellobody & Bailey, 2009Ellobody, E. & Bailey, C. G., 2009. Modelling of unbonded post-tensioned concrete slabs under fire conditions. Fire Safety Journal, 44(2), pp. 159-167.), among others. Limited research was found, though, investigating the use of FEM models to predict the mechanical properties of lightweight concretes out of their raw materials. A relevant research involving this subject was carried out by Chung et al (2018Chung, S. Y., Elrahman, M. A., Stephan, D. & Kamm, P. H., 2018. The influence of different concrete additions on the properties of lightweight concrete evaluated using experimental and numerical approaches. Construction & Building Materials, Volume 189, pp. 314-322.), where a commercial FE package was used to evaluate mechanical and thermal properties of LWACs.

In this context, an import question emerges - which kind of methodology allows better prediction of the mechanical properties of LWACs: machine learning techniques or classical mechanical-based models? Therefore, seeking to improve the safety and contribute to the design of LWAC structures, the present work evaluates two methods to predict the compressive strength (f c ) and Young’s Modulus (E c ) of this material from the characteristics of its mortar and aggregates. They are: Support Vector Machine with Regression (SVR) and Finite Element Method (FEM).

2. Methodology

2.1 Finite Element Method

In order to apply the Finite Element Method, LWAC was modeled as a biphasic material: the mortar (m) and the lightweight aggregates (a). The experimental data used as input and for validation were obtained from an experimental study with cylindrical samples of LWAC. The tests were performed by Ke et al. (2009Ke, Y. et al., 2009. Influence of volume fraction and characteristics of lightweight aggregates on the mechanical properties of concrete. Construction and Building Materials, 23(8), pp. 2821-2828.) (2014), who analyzed compressive strength and modulus of elasticity of LWAC samples at 28 days. In this experimental program, three types of mortars were used: normal (40 MPa), high performance (64 MPa), and very high performance (86 MPa). Regarding the aggregates, five types were used: 0/4-650-A, 4/10-550-A, 4/10-430-A, 4/10-520-S and 4/8-750-S; in which each name matches the characteristic diameters (d/D, in mm), the bulk density (kg/m3), and the type of aggregate (A - expanded clay, S - expanded shale), respectively. The combination of these materials in proportions of 12.5%, 25.0%, 37.5% and 45.0% of LWA in volume generated 45 different formulations and over 180 specimens were tested. Each LWAC sample is also characterized by a water/cement factor, a quantity of cement, aggregate volume, aggregate density, Young’s modulus and compressive strength.

From the aggregate and mortar properties we simulated the behavior of LWAC samples at 28 days when subjected to compressive load. The free processing software CAST3M was used, which includes built-in pre-processing and post-processing tools. The parameters used in the FEM simulation were:

  • Compressive strength (f m ) and Young modulus (E m ) of the mortar matrix (experimentally obtained by Ke et al. (2014Ke, Y., Ortola, S., Beaucour, A. L. & Dumontet, H., 2014. Micro-stress analysis and identification of lightweight aggregate’s failure strength by micromechanical modeling. Mechanics of Materials, Volume 68, pp. 176-192.));

  • Tensile strength (t m ) of the mortar, based on empirical correlations with its compressive strength proposed by Chust e Figueiredo Filho (2004Chust, C. R. & Figueiredo Filho, J. R., 2004. Cálculo e Detalhamento de Estruturas Usuais de Concreto armado. São Carlos: Editora da UFSCar.);

  • Volume fraction of lightweight aggregate adopted in the concrete’s mixture (12.5%, 25.0%, 37.5% and 45.0%);

  • Particle size distribution of the aggregates (Ke, 2008Ke, Y., 2008. Caractérisation du comportement mécanique des bétons de granulats légers: expérience, Cergy: University of Cergy-Pontoise.);

  • Young modulus of the LWA (E a ), obtained empirically, based on its dry density (Ke et al., 2014Ke, Y., Ortola, S., Beaucour, A. L. & Dumontet, H., 2014. Micro-stress analysis and identification of lightweight aggregate’s failure strength by micromechanical modeling. Mechanics of Materials, Volume 68, pp. 176-192.);

  • Compressive strength of the LWA (f a ), obtained in an analytical inverse method by Ke et al. (2014Ke, Y., Ortola, S., Beaucour, A. L. & Dumontet, H., 2014. Micro-stress analysis and identification of lightweight aggregate’s failure strength by micromechanical modeling. Mechanics of Materials, Volume 68, pp. 176-192.).

  • Tensile strength (t a ) of the LWA, from (Ke et al., 2009Ke, Y. et al., 2009. Influence of volume fraction and characteristics of lightweight aggregates on the mechanical properties of concrete. Construction and Building Materials, 23(8), pp. 2821-2828.);

  • For comparison purposes, compressive strength (f c,exp ) and Young modulus (E c,exp ) measured on hardened LWAC samples (Ke et al., 2014Ke, Y., Ortola, S., Beaucour, A. L. & Dumontet, H., 2014. Micro-stress analysis and identification of lightweight aggregate’s failure strength by micromechanical modeling. Mechanics of Materials, Volume 68, pp. 176-192.).

The FEM model considers a longitudinal cross-section of a standard cylindrical sample with 16 cm of diameter and 32 cm high (Fig. 1a). For computational efficiency, 1/4 of the section was used in this work, considering the boundary conditions of geometric symmetry of the sample. In this specimen, the LWA were modelled as spheres randomly distributed within the homogeneous mortar matrix (Fig. 1b), with the same particle size distribution as the experimental program. In order to account for the dispersion of results, three models with random aggregate scattering (within the particle size distribution) were generated from each synthetic sample.

An elastoplastic material behavior was adopted for the FEM simulations, with Drucker-Prager yield criterion (Commissariat a l’Energie Atomic, 2003Commissariat a l’Energie Atomic, 2003. Notices | cast3m - mode. [Online] Available at: http://www-cast3m.cea.fr/index.php?page=notices&notice=MODE [Accessed 15 07 2018].
http://www-cast3m.cea.fr/index.php?page=...
). Limits for compressive and tensile strength for each material were assumed according to the input parameters, and a perfect plastic yield was considered. The materials properties adopted for simulation and validation of these models are shown in Table 1.

To simulate the compressive strength test, a downward linear displacement of 0.5 mm/s was imposed to the upper surface. The rate of 100 computational steps per 0.1 second was adopted. The compressive stress values were obtained from the support conditions at the base of the sample. The concrete is considered to reach failure at a 3.0‰ strain (Sussekind, 1985Sussekind, J., 1985. Curso de concreto, Volume 1. Porto Alegre: Editora Globo.).

Fig. 1
Numerical representation of one LWAC sample

Table 1
Summary of materials properties used in the FEM simulation

2.2 Support Vector Machine with Regression

The models based on SVR method were applied to the set of experimental results for 12.5%, 25.0%, 37.5% and 45.0% volume of LWA (total of 180 mixtures of LWAC). From the total data, 2/3 were used to train the SVR model, and prediction results were obtained from the remaining 1/3. Four attributes were used as input variables: water/cement ratio, quantity of cement, volume of aggregate and density of aggregate; and two as output variables: LWAC’s compressive strength (f c ) and Young’s modulus (E c ). They are shown in Table 2.

Table 2
Summary of attributes used in the SVR evaluation

To minimize bias associated with random data samples for training and testing, cross validation was applied by the k-folding technique. This method divides the total data set into k mutually exclusive subsets of the same size and, subsequently, a subset is used for testing and the remaining k-1 are employed to estimate the parameters, what allows for evaluating the accuracy of the model (Kohavi, 1995aKohavi, R., 1995a. The power of decision tables. Berlin, Springer, pp. 174-189.). This process is performed k times alternating the test subsets. The k-fold technique is widely used in modeling problems concerning value prediction, such as the concrete properties (Chou & Tsai, 2012Chou, J. S. & Tsai, C. F., 2012. Concrete compressive strength analysis using a combined classification and regression technique. Automation in Construction, Volume 24, pp. 52-60.) (Fushiki, 2011Fushiki, T., 2011. Estimation of prediction error by using K-fold cross-validation. Statistics and Computing, 21(2), pp. 137-146.) (Chou et al., 2010Chou, J. S., Chiu, C. K., Farfoura, M. & Al-Taharwa, I., 2010. Optimizing the prediction accuracy of concrete compressive strength based on a comparison of data-mining techniques. Journal of Computing in Civil Engineering, 25(3), pp. 242-253.).

In the present work, prediction methods were applied in a 10-fold cross-validation, which are favored from a computational standpoint (Kohavi, 1995bKohavi, R., 1995b. A study of cross-validation and bootstrap for accuracy estimation and model selection. s.l., Morgan Kaufmann Publishers Inc., p. 1137-1143.) (Fushiki, 2011Fushiki, T., 2011. Estimation of prediction error by using K-fold cross-validation. Statistics and Computing, 21(2), pp. 137-146.). This procedure was repeated 100 times with a different test set each round, aiming to obtain dispersion results.

The parameters used to execute the SVR for output variable f c were: kernel is RBF (Radial Basis Function), C = 100 and γ = 0.03162278; and for output variable E c were: kernel is RBF, C = 10 and γ = 0.06309573; where C is the cost and γ is the gamma. These parameters were optimized via grid search technique, which performs an exhaustive search to find the best setting to switch all input values (Chen et al., 2000Chen, D. et al., 2000. Grid search as applied to the determination of Mark-Houwink parameters. Journal of applied polymer science, 76(4), pp. 481-487.).

For comparison purposes, besides SRV, three other supervised learning methods based on adjustment of parameters in a set experimental results were used. These are: Multiple Linear Regression (MLR), Quadratic Multiple Linear Regression (MR2) and Cubic Multiple Linear Regression (MR3). All the methods in this study were developed using the R software.

In addition, to evaluate the performance results, two different statistical parameters were employed: Root Mean Square Error (RMSE) and Mean Absolute Percentage Error (MAPE). They are mathematically defined by Equations 1 and 2, respectively, where n is the number of samples, y i is the observed value and p i is the predicted value.

R M S E = 1 n i = 1 n ( y i p i ) 2 (1)

M A P E = 1 n i = 1 n ( | y i p i y i | ) (2)

The RMSE compares observed values with their predicted counterparts and calculates the average of the squared residual errors. However, the RMSE weights large errors more significantly. Alternatively, MAPE statistical metric is dimensionless and provides an effective means of residual error compared to each observed value with their respective predicted value. The lower the RMSE and MAPE values, the better the performance of the prediction model.

4. Results

2.1 Finite Element Method

In the FEM simulations, the compressive strength and Young's modulus were evaluated from 180 models (45 for each aggregate volume). In every step of computing time, the imposed displacement generates a stress increase in the elements of each phase (mortar and aggregates). Therefore, the stress developed in the elements can be normalized, according to each yield strength, in levels from 0 to 1. This evolution can be followed in Fig. 2, which represents the tension levels related to deformations of 0.065‰, 0.34‰, 0.81‰, 1.10‰, 1.50‰ and 3.00‰ of a specimen with 12.5% of LWA type 4/10-550-A and normal mortar.

Fig. 2
Evolution of the stress levels in FEM simulations related to imposed displacement. Each material is standardized in accordance with their respective tensile and compressive strength.

In turn, the stress-strain curve of a LWAC sample with 12.5% of LWA type 4/10-550-A is shown in Fig. 3, along with the isolated curves from its mortar and aggregates. It can be seen from Fig. 3 that for strains up to approximately 1.3‰, the LWAC behavior is elastic, and Fig. 2a, Fig. 2b, Fig. 2c and Fig. 2d show the evolution of the stress levels in this phase. At strains from 1.3‰ to 2.5‰, there is a transitional phase (Fig. 2e) and after 2.5‰, the specimen has reached the yield point that remains constant until the end of imposed displacement. Fig. 2f shows the final configuration assumed for the imposed displacement, when most elements have reached their yield strength.

Fig. 3
FEM results of stress-strain curve of a sample with 12.5% of LWA type 4/10-550-A with normal mortar, along with the isolated curves from its mortar and aggregates

The simulation results for all LWAC models are summarized in Table 3, which includes minimum and maximum values, mean values, and the mean of relative variations from reference parameters (represented by the MAPE value). The results obtained are close to the experimental reference: the larger mean deviation found for f c was 7.25%. On the other hand, the E c simulations presented results more similar to the experimental values, with mean deviation lower than 3.32%. No bias was detected regarding the increasing volume of aggregate.

Table 3
Results of compressive strength (f c ) and Young’s modulus (E c ) obtained from FEM simulations and relative variation to the experimental reference

The fact that LWA are already approximately rounded, thus similar to the model adopted for the aggregate phase, possibly contributed to the positive results. Although within a satisfactory range, the deviations observed are probably a function of the simplicity of model regarding the behavior of the materials, the reliability and significance of the experimental methods used to obtain the input parameters from the materials, and the heterogeneous nature of real concrete specimens (Lamond & Pielert, 2006Lamond, J. F. & Pielert, J. H., 2006. Significance of tests and properties of concrete and concrete-making materials. West Conshohocken: ASTM International.). Since the fracture mechanism close to the material yield stress is a more complex and cascading phenomenon, it is understandable that the numerical simulation obtained better results in the predictions regarding the initial elastic phase of the LWAC model (E c ).

In addition, the experimental procedure for E c requires a significant larger effort in terms of time, cost and specimens than testing for f c . In this sense, the FEM simulation can contribute to speeding up laboratory results in the development of new mix designs. Finally, the results obtained from FEM also allow extrapolation of these concrete models to more complex structures such as beams, slabs, sealing elements, among others, in order to simulate real world phenomena.

Table 4 shows the residual errors for the FEM simulations. It is observed that the results’ dispersion is comparable to that obtained in laboratory experiments, which denotes the adequacy of the numerical procedures for practical purposes. Furthermore, Fig. 4 display residual errors histograms with a high degree of symmetry for each output variable, indicating that they are correspond to a noise.

Table 4
Residual errors of prediction for 180 simulations for FEM.

Fig. 4
Residual errors histogram for 180 simulations for FEM.

2.2 Support Vector Machine with Regression

The result of SVR predictions obtained from 100 runs of 10-fold technique with all experimental data can be seen in Table 5. Table 6 shows the relative variations (predicted versus experimental) from two statistic metrics of results obtained with SVR, MLR, MR2 and MR3.

Table 5
Results of compressive strength (f c ) and Young’s modulus (E c ) obtained from SVR simulations and relative variation to the experimental reference

Table 6
Comparison of the relative variations of the prediction of compressive strength (f c ) and Young’s modulus (E c ) obtained with SVR, MLR, MR2 and MR3.

It can be seen from Table 5 that the SVR results are also close to the experimental references: the mean deviation of f c was below 7.4% and for E c , lower than 3.63%. Again, the increasing volume of aggregate did not explicitly influence the results. The small discrepancy is possibly due to minor adjustments issues regarding the prediction model (kernel, C e γ parameters) and the limitations of data used as training. In addition, care must be taken in seeking excessive error minimization, as overfitting could bring limitations to the prediction model, which could lose its ability to predict data different from those used in this study.

It can be observed that the values obtained by MR2 and MR3 methods were significantly better than the MLR method, since the lower the RMSE and MAPE, the better the result. Therefore, it can be concluded that there is nonlinearity between each input variable to the output variables, what reinforces the choice of the SRV technique.

Finally, it can be seen from Table 7 that the minimum and maximum values of residual errors of prediction have no discrepancy and no bias for both output variables. Furthermore, Fig. 5 shows a high degree of symmetry of the residual errors for each output variable, indicating that they are similar to a noise.

Table 7
Residual errors of prediction for 100 runs of 10-fold for SVR.

Fig. 5
Residual errors histogram for 100 runs of 10-fold for SVR.

5. Discussions and Conclusions

In this study, two techniques for predicting the mechanical properties of Lightweight Aggregate Concrete (LWAC) have been developed and compared. The mechanical properties of interest were the compressive strength (f c ) and static Young’s Modulus (E c ). The Finite Element Method (FEM) technique sought to reproduce computationally, through a 2-phase model, laboratory tests for these mechanical properties. In turn, the Support Vector Machine with Regression (SRV) technique is based on parameter settings to a set of experimental results, i.e. a supervised learning technique.

The advantages of the SVM to predict the concrete properties are its speed and the possibility for continuous re-training with new information. On the other hand, this technique requires a comprehensive amount of input data to be reliable. Regarding FEM, although slower, the method does not require previous final results. It provides satisfactory predictions with fewer inputs, given that the computational models are suitably constructed.

In order to better compare both techniques, Fig. 6 and Fig. 7 present a summary of all achieved results, by applying FEM and SVM, respectively. These figures exhibit experimental data on the horizontal axis and mean numerical results on the vertical axis. It is possible to observe that both techniques allow predictions for f c and E c that may be considered as identical for design purposes. These results allow us to answer the question highlighted in the Introduction of this paper (retyped): “which kind of methodology allows better prediction of the mechanical properties of LWAC: machine learning techniques or classical mechanical-based models?” The answer is: both.

Both techniques of computer modeling and computational intelligence showed good agreement with the experimental reference. SVR presented slightly better performance in predicting the mechanical properties of LWAC, with average deviations from the experimental values of 5.46% for f c and 2.48% for E c . These results are lower than the 6.55% and 2.91% obtained by FEM for f c and E c , respectively. The deviations of up to 7.25% for f c in the FEM model can be justified by the simplifications assumed in the modeling of the sample and the reliability of the input parameters. Alternatively, the deviations up to 7.37% on the SRV analysis are possibly due to minor adjustments issues regarding the prediction model and the data used. Overall, the SVR technique showed a lower average error compared to the FEM technique and both presented a better prediction of the E c than the f c . This is expected, since the E c is obtained in the initial elastic stage of the specimens and is thus less subjected to intricate failure mechanisms and the consequences of random heterogeneities.

The two different statistical parameters used - RMSE and MAPE - yielded a comprehensive comparison of the prediction methods employed. The comparison among SVR, MLR, MR2 and MR3 confirm that the LWAC prediction problem is nonlinear and thus cannot be solved by traditional linear regression techniques (e.g. MLR). According to the methodology adopted, the best method of computational intelligence to predict the mechanical properties of LWAC is SVR, with an average MAPE of 4.5% and 2.4% for f c and E c , respectively; and average RMSE of 2.3% and 0.7% for f c and E c , respectively.

In summary, both FEM and SVR were technically feasible to estimate the compressive strength and modulus of elasticity of LWAC with minor deviations. Considering that experimental procedures for obtaining these properties require a significant amount of time, cost and materials, these methods can contribute to a more safe, economic and sustainable design of concrete mixtures.

Fig. 6
Comparison between experimental and predicted FEM results

Fig. 7
Comparison between experimental and predicted SVR results

Acknowledgments

The authors gratefully acknowledge CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior), FAPEMIG (Fundação de Amparo à Pesquisa do Estado de Minas Gerais), CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico), UFJF (Universidade Federal de Juiz de Fora), and UFOP (Universidade Federal de Ouro Preto) for providing infrastructure and financial support.

References

  • Abd, A. M. & Abd, S. M., 2017. Modelling the strength of lightweight foamed concrete using support vector machine (SVM). Case studies in construction materials, Volume 6, pp. 8-15.
  • ABNT, 2015. NBR 8953 - Concrete for structural use - Density, strength and consistence classification, Rio de Janeiro: Associação Brasileira de Normas Técnicas.
  • Akande, K. O., Owolabi, T. O., Twaha, S. & Olatunji, S. O., 2014. Performance comparison of SVM and ANN in predicting compressive strength of concrete. IOSR Journal of Computer Engineering, 16(5), pp. 88-94.
  • Al-Rub, R. K. A. & Kim, S. M., 2010. Computational applications of a coupled plasticity-damage constitutive model for simulating plain concrete fracture. Engineering Fracture Mechanics, 77(10), pp. 1577-1603.
  • ASTM, 2018. ASTM C125 - Standard Terminology Relating to Concrete and Concrete Aggregates, West Conshohocken: ASTM International.
  • Bentur, A., Igarashi, S. I. & Kovler, K., 2001. Prevention of autogenous shrinkage in high-strength concrete by internal curing using wet lightweight aggregates. Cement and concrete research, 31(11), pp. 1587-1591.
  • Chen, D. et al., 2000. Grid search as applied to the determination of Mark-Houwink parameters. Journal of applied polymer science, 76(4), pp. 481-487.
  • Cheng, M. Y., Chou, J. S., Roy, A. F. & Wu, Y. W., 2012. High-performance concrete compressive strength prediction using time-weighted evolutionary fuzzy support vector machines inference model. Automation in Construction, Volume 28, pp. 106-115.
  • Chou, J. S., Chiu, C. K., Farfoura, M. & Al-Taharwa, I., 2010. Optimizing the prediction accuracy of concrete compressive strength based on a comparison of data-mining techniques. Journal of Computing in Civil Engineering, 25(3), pp. 242-253.
  • Chou, J. S. & Tsai, C. F., 2012. Concrete compressive strength analysis using a combined classification and regression technique. Automation in Construction, Volume 24, pp. 52-60.
  • Chung, S. Y., Elrahman, M. A., Stephan, D. & Kamm, P. H., 2018. The influence of different concrete additions on the properties of lightweight concrete evaluated using experimental and numerical approaches. Construction & Building Materials, Volume 189, pp. 314-322.
  • Chust, C. R. & Figueiredo Filho, J. R., 2004. Cálculo e Detalhamento de Estruturas Usuais de Concreto armado. São Carlos: Editora da UFSCar.
  • Commissariat a l’Energie Atomic, 2003. Notices | cast3m - mode. [Online] Available at: http://www-cast3m.cea.fr/index.php?page=notices&notice=MODE [Accessed 15 07 2018].
    » http://www-cast3m.cea.fr/index.php?page=notices&notice=MODE
  • Ellobody, E. & Bailey, C. G., 2009. Modelling of unbonded post-tensioned concrete slabs under fire conditions. Fire Safety Journal, 44(2), pp. 159-167.
  • Fushiki, T., 2011. Estimation of prediction error by using K-fold cross-validation. Statistics and Computing, 21(2), pp. 137-146.
  • Gao, W. Y., Dai, J. G. T. J. G. & Chen, G. M., 2013. Finite element modeling of reinforced concrete beams exposed to fire. Engineering structures, Volume 52, pp. 488-501.
  • Genikomsou, A. S. & Polak, M. A., 2015. Finite element analysis of punching shear of concrete slabs using damaged plasticity model in ABAQUS. Engineering Structures, Volume 98, pp. 38-48.
  • Gupta, R., Kewalramani, M. A. & Goel, A., 2006. Prediction of concrete strength using neural-expert system. Journal of materials in civil engineering, 18(3), pp. 462-466.
  • Gupta, S. M., 2007. Support vector machines based modelling of concrete strength. World Academy of Science, Engineering and Technology, Volume 36, pp. 305-311.
  • Hassan, N. Z., Sherif, A. G. & Zamarawy, A. H., 2015. Finite element analysis of reinforced concrete beams with opening strengthened using FRP. Ain Shams Engineering Journal.
  • Ke, Y., 2008. Caractérisation du comportement mécanique des bétons de granulats légers: expérience, Cergy: University of Cergy-Pontoise.
  • Ke, Y. et al., 2009. Influence of volume fraction and characteristics of lightweight aggregates on the mechanical properties of concrete. Construction and Building Materials, 23(8), pp. 2821-2828.
  • Ke, Y., Ortola, S., Beaucour, A. L. & Dumontet, H., 2014. Micro-stress analysis and identification of lightweight aggregate’s failure strength by micromechanical modeling. Mechanics of Materials, Volume 68, pp. 176-192.
  • Kim, D. K., Lee, J. J., Lee, J. H. & Chang, S. K., 2005. Application of probabilistic neural networks for prediction of concrete strength. Journal of Materials in Civil Engineering, 17(3), pp. 353-362.
  • Kohavi, R., 1995a. The power of decision tables. Berlin, Springer, pp. 174-189.
  • Kohavi, R., 1995b. A study of cross-validation and bootstrap for accuracy estimation and model selection. s.l., Morgan Kaufmann Publishers Inc., p. 1137-1143.
  • Lamond, J. F. & Pielert, J. H., 2006. Significance of tests and properties of concrete and concrete-making materials. West Conshohocken: ASTM International.
  • Lee, S. C., 2003. Prediction of concrete strength using artificial neural networks. Engineering Structures, 25(7), pp. 849-857.
  • Mehta, P. & Monteiro, P. J. M., 2014. Concrete: Microstructure, Properties, and Materials. 2ª ed. São Paulo: IBRACON.
  • Nazari, A. & Sanjayan, J. G., 2015. Modelling of compressive strength of geopolymer paste, mortar and concrete by optimized support vector machine. Ceramics International, 41(9), pp. 12164-12177.
  • Ni, H. G. & Wang, J. Z., 2000. Prediction of compressive strength of concrete by neural networks. Cement and Concrete Research, 30(8), pp. 1245-1250.
  • Sadrmomtazi, A., Sobhani, J. & Mirgozar, M. A., 2013. Modeling compressive strength of EPS lightweight concrete using regression, neural network and ANFIS. Construction and Building Materials, Volume 42, pp. 205-216.
  • Schlangen, E. & Van Mier, J. G. M., 1992. Simple lattice model for numerical simulation of fracture of concrete materials and structures. Materials and Structures, 25(9), pp. 534-542.
  • Sussekind, J., 1985. Curso de concreto, Volume 1. Porto Alegre: Editora Globo.
  • Tenza-Abril, A. J., Villacampa, Y., Solak, A. M. & Baeza-Brotons, F., 2018. Prediction and sensitivity analysis of compressive strength in segregated lightweight concrete based on artificial neural network using ultrasonic pulse velocity. Construction and Building Materials, Volume 189, pp. 1173-1183.
  • Yang, C. C. & Huang, R., 1998. Approximate strength of lightweight aggregate using micromechanics method. Advanced Cement Based Materials, 7(3-4), pp. 133-138.
  • Yan, K. & Shi, C., 2010. Prediction of elastic modulus of normal and high strength concrete by support vector machine. Construction and building materials, 24(8), pp. 1479-1485.
  • Yu, T. T. J. G., Teng, J. G., Wong, Y. L. & Dong, S. L., 2010. Finite element modeling of confined concrete-I: Drucker-Prager type plasticity model. Engineering Structures, 32(3), pp. 665-679.
  • Zhang, M. H. & Gjvorv, O. E., 1991. Mechanical properties of high-strength lightweight concrete. ACI Materials Journal, 88(3), pp. 240-247.

Edited by

Editor:

Adnan Ibrahimbegovic; Paulo de Mattos Pimenta.

Publication Dates

  • Publication in this collection
    12 Aug 2019
  • Date of issue
    2019

History

  • Received
    20 Sept 2018
  • Reviewed
    28 Feb 2019
  • Accepted
    04 June 2019
  • Published
    13 June 2019
Individual owner www.lajss.org - São Paulo - SP - Brazil
E-mail: lajsssecretary@gmsie.usp.br