Modeling shear strength of medium- to ultra-high-strength concrete beams with stirrups using SVR and genetic algorithm

This paper presents a data-driven machine learning approach of support vector regression (SVR) with genetic algorithm (GA) optimization approach called SVR-GA for predicting the shear strength capacity of medium- to ultra-high strength concrete beams with longitudinal reinforcement and vertical stirrups. One hundred and forty eight experimental samples collected with different geometric, material and physical factors from literature were utilized for SVR-GA with fivefold cross validation. Shear influence factors such as the stirrup spacing, the beam width, the shear span-to-depth ratio, the effective depth of the beam, the concrete compressive and tensile strength, the longitudinal reinforcement ratio, the product of stirrup ratio and stirrup yield strength were served as input variables. The simulation results show that SVR-GA model can achieve highest accuracy in shear strength prediction based on testing set with a coefficient of determination (R2) of 0.9642, root mean squared error of 1.4685 and mean absolute error of 1.0216 superior to that for traditional SVR model with 0.9379, 2.0375 and 1.4917, which both perform better than multiple linear regression and ACI-318. Furthermore, the sensitivity analysis reveals the most important variables affecting the result of shear strength prediction are shear span-to-depth ratio, concrete compressive strength, reinforcement ratio and the product of stirrup ratio and stirrup yield strength. Three-dimensional input/output maps are employed to reflect the nonlinear variation of the shear strength with the two coupling variables. All in all, the proposed SVR-GA model can achieve excellent accuracy in prediction the shear strength of medium- to ultra-high strength concrete beams with stirrups in comparison with results obtained by traditional SVR, MLP and ACI-318.


Introduction
The shear failure of reinforced concrete beams with stirrups is a common concern of structural engineers (Collins et al. 2008;Sagaseta and Vollum 2011;Słowik 2014). However, it is difficult to predict the shear failure accurately due to the influence of a large number of parameters, such as stirrup spacing, width and effective depth of the beam, shear span-to-depth ratio, stirrup ratio, longitudinal reinforcement ratio, tensile compressive strength of concrete, and stirrup yield strength. This difficulty is particularly evident in ultra-high strength concrete (UHSC) and ultra-high performance concrete (UHPC) beams (Hossain et al. 2017).
In order to accurately estimate the shear capacity of UHPC beams, the shear capacity is artificially divided into concrete shear capacity and stirrup yield shear capacity. An additional shear contribution of steel fiber would be added while adding steel fiber into concrete. Further, an additional shear contribution of the pin would be superimposed when taking the pin action of longitudinal reinforcement into account (Yoo and Yoon 2016;Marì Bernat et al. 2020). However, these factors are not independent of each other, and there is a coupling effect between them. For example, the residual tensile stress between cracks will be increased by appropriately adding the steel fiber while both increasing the shear contribution of the fiber and concrete.
A variety of normative formulas or models are proposed to solve problems in engineering applications. However, there are great differences in selecting main variables affecting the shear strength, such as the code of China for the design of concrete structures (GB 50010-2010, Press 2010 uses the concrete tensile strength to calculate the shear strength, while the American concrete structure design code (ACI 318-14, ACI Committee 318 2014) adopts the concrete compressive strength for that. The tensile strength and compressive strength of concrete are both considered by the Chinese highway and bridge code (JTG 3362-2018, Ministry of Transport of China 2018). Although there is little difference between UHPC and normal concrete in the ratio of tensile strength to compressive strength. The addition of steel fibers has greatly affected the tensile strength of UHPC while almost no influence on the compressive strength (Hassan et al. 2012;Krassowska et al. 2019). It is noteworthy that most of the formulas in the codes have been tested on the limited data which is just an extension of the existing empirical formulas for the shear strength of medium and high strength concrete beams, without fully considering and utilizing the ultra-high mechanical properties of ultra-high strength concrete.
Traditional models/equations with low accuracy mainly rely on basic expressions and step-by-step refinement process, it is urgent to propose more accurate method to calculate the shear strength of UHPC beams. Recently, the data-driven machine learning (ML) methods have attracted extensive attention, because of their inherent ability of overcoming the shortcomings of traditional algorithms by relying on manual design features (Liu et al. 2017). Artificial neural network (Açikgenç et al. 2015;Golafshani et al. 2015;Hossain et al. 2017), adaptive fuzzy neural network (Mansouri et al. 2016;Nguyen et al. 2020), Gaussian process regression (Hoang et al. 2016;Guo and Hesthaven 2018) and support vector machine (Pal and Deswal 2011;Farfani et al. 2015) have been widely applied to model the mechanical properties of concrete structures.
ML model uses samples with input and output for training models that can be used to predict the output of new inputs. At present, ML models have been employed to predict the shear capacity of normal concrete beams or ultra-high strength concrete beams without stirrups (Solhmirzaei et al. 2020;Zhang et al. 2020). However, few studies were found for applying machine learning approach to predict the shear capacity of ultra-high strength concrete beams with stirrups. Support vector machine (SVM) is one of ML methods based on statistical learning theory of structural risk minimization. Also known as maximum interval classifier, it is very popular for solving multi-class problems by establishing a group of 2-class classifiers or establishing a multi-class classifier (Teng et al. 2018). It is called support vector regression (SVR) for solving regression problems. SVR has the advantages of fast learning, global optimization, avoiding local minimization and excellent generalization ability (Ç evik et al. 2015), while it is difficult to determine proper parameters for SVR model (Yu et al. 2006;Ccoicca 2013). GA can not only be used to solve the optimization problem directly (Deveci and Demirel 2016), but also optimize the model parameters (Yan and Lin 2016). Here, GA is introduced to simulate biological evolution process to search globally optimal hyper-parameter for SVR.
The remainder of this paper is organized in the following order. Section 2 introduces the background overview, Sect. 3 shows the data structure including the collected experimental results, Sect. 4 describes the proposed SVR-GA model, Sect. 5 presents the model implementation and performance metrics, Sect. 6 presents a performance comparison of SVR-GA, SVR, MLP and ACI 318, the sensitivity analysis for the input variables and the experimental results and analysis. Finally, Sect. 7 provides the conclusions of this study.

Literature review
With the wide application of high strength concrete, the shear capacity of high-strength and ultra-high strength reinforced concrete beam is one of the most important calculation problems in the design of reinforced concrete frames. In order to reveal and predict the performance of high strength beam and to find the most important parameters affecting the shear capacity, a number of studies have been carried out in the last 2 decades. Zhang (2005) used the modified compression field theory model to predict the shear capacity of high strength concrete beams. The influence of steel fiber and stirrups on the shear strength of steel fiber high-strength concrete beams was analyzed. The results showed that it is not reasonable for steel fiber to completely replace stirrups to improve the shear capacity of beams. Compared with ordinary concrete beam, high strength concrete beam has higher shear strength and better deformation capacity (Ji et al. 2011), and the shear strength decreases with the increase in the shear span-to-depth ratio while increases with the increase in longitudinal reinforcement ratio and stirrup ratio. For the ultra-high performance fiber reinforced concrete (UHPFRC) beam, the stirrup spacing exceeds the spacing limits recommended by the design code of ACI 318, but the test values of shear strength are significantly larger than the calculated values of ACI 318 (Lim and Hong 2016). The stirrup spacing required by ACI 318 for reinforced concrete beams is not suitable for UHPFRC deep beams (Yousef et al. 2018). The tensile strength of UHPC can be increased when adding a certain amount of steel fiber, but that has little effect on the compressive strength. In addition, steel fiber can enhance the ultimate shear strength of reinforced concrete beams (Krassowska et al. 2019). Based on the concept of piecemeal superposition, Qi et al. (2020) established the theoretical calculation formula for the shear capacity of UHPC beams considering the contribution of concrete, stirrups and fibers. The analysis results show that the formula is in good agreement with the experimental data only containing 11 samples.
Pal and Deswal (2011) established a SVR model to predict the shear strength of prestressed reinforced concrete deep beams. The results showed that SVR model is superior to empirical relationship and back propagation neural network in terms of prediction ability. And the concrete cylinder strength and shear span-to-depth ratio are of great significance to the strength prediction of deep beams. Some researchers used intelligent algorithms to optimize model parameters, such as genetic algorithm (GA) to optimize a portion of the existing method (Zhang et al. 2020), and the results obtained by this method are improved compared to other machine learning methods including classical SVR, ANN and gradient boosted decision trees (GBDTs). Solhmirzaei et al. (2020) presented a data-driven ML framework for the failure mode and shear strength capacity of UHPC beams. The results showed that different ML algorithms SVM, ANN, and K-nearest neighbor (K-NN) were effective to classify the shear, flexure-shear, and flexure failure modes. SVM and K-NN algorithms were better than ANN in identifying shear failure mode and flexure-shear failure modes. The explicit expression proposed by GP can predict the shear capacity of UHPC beams with R 2 of 0.92.
One concern of ML model is that the input data are sparse matrices of high dimensions, which are mostly missing or zero (Luo et al. 2018a). To solve this problem of the automatic web-service selection, making highly accurate predictions for missing quality of service (QoS) data via building an ensemble of non-negative latent factor models (Luo et al. 2016). In order to further generate highly accurate predictions for missing QoS data, the second-order solver into latent factor models are proposed (Luo et al. 2018b). Recently, the momentum-incorporated parallel stochastic gradient descent algorithm has also been applied to solve sparse matrix problems by implementing parallelization and accelerating convergence rate (Luo et al. 2021).
Another main problem of SVR is to choose proper parameters and kernel function, which needs to take time for trial and error (Yu et al. 2006;Ccoicca 2013). One method is to solve the problem of optimal hyper-parameter determination through grid search (Zhang et al. 2018), which is time-consuming while traversing all grids within the parameter value range. Another efficient way to do this is through genetic algorithm (GA), which has been used as a powerful optimization tool to solve a variety of academic and engineering problems. For example, to solve the airline crew pairing problem, an improved dynamic-based GA has been proposed Demirel 2016, 2018;Demirel and Deveci 2017). To predict the failure of the RC structure, GA are explored the optimal initial weights and biases for ANN with single and multiple objective to avoid falling into local minima (Yan and Lin 2016;Chatterjee et al. 2017;Umeonyiagu and Nwobi-Okoye 2019).
It is concluded ML models can be useful tool in predicting shear capacity of concrete beams, and GA can be an effective method to optimize parameters for ML models, which is introduced to select optimal parameters for SVR. Few studies were conducted to apply ML approach to predict the shear strength of ultra-high strength concrete beams with stirrups due to the complexity of the calculation. Therefore, we proposed a hybrid SVR-GA model to solve this problem.
In this study, a hybrid SVR-GA is proposed to predict the shear strength of medium-to ultra-high-strength strength concrete beams with stirrups. The major contributions are as follows: 1) Employment of a ML model combining SVR with GA optimization approach to predict the shear strength of medium-to ultra-high-strength strength concrete beams with stirrups. 2) Comparison results of ACI-318 and ML models including SVR-GA, SVR and multiple linear regression (MLR) to obtain the model with highest precision for shear strength prediction. 3) Utility of parametric sensitivity analysis with threedimensional visualization map to reveal the inputs that have significant influence on shear strength prediction.
Our proposed model can provide a guide with high accuracy for the performance analysis and optimal design of the shear strength of medium-to ultra-high strength concrete beams with stirrups.

Data construction
The experimental dataset in this study was collected from previously published works discussing about the shear strength of the simply supported beam (Zhang 2005;Magureanu et al. 2010;Ji et al. 2011;Baby et al. 2014;Xu et al. 2014;Kamal et al. 2014;Zhou and Chen 2015;Hou et al. 2015;Jin et al. 2015;Lim and Hong 2016;Pansuk et al. 2017;Smarzewski 2018;Yousef et al. 2018;Mészöly and Randl 2018;Hasgul et al. 2019;Zheng et al. 2019;Krassowska et al. 2019;Qi et al. 2020;Wang et al. 2020). The typical RC beam with stirrups and its geometric parameters are illustrated in Fig. 1. The paper review helped determine the factors affecting the beam's shear capacity (Russo et al. 2004;Olalusi and Viljoen 2020). Various related factors used to predict ultimate shear capacity (V u ) are shown in Table 1, including different geometric, material and physical factors, such as stirrup spacing (s), beam width (b), shear span-to-depth ratio (a/d), effective depth of the beam (d), concrete compressive strength (f c ), concrete tensile strength (f t ), longitudinal reinforcement ratio (q), as well as the product of the stirrup ratio and the stirrup yield strength (q sv f yv ). In this study, the normalized ultimate shear strength [i.e., v u = V u /(bd) (Zhang et al. 2020)] was used as a measure for evaluating the shear resistance of the beam.

Support vector regression (SVR)
Suppose a sample is represented by {(x i , y i ), i = 1, 2, …, n}, where x i and y i correspond to its input and output, respectively. In SVR, the fundamental idea is to map xi to a high-dimensional feature space F through a nonlinear mapping function U(x i ). The goal of SVR is to perform linear regression in this space and find the following linear equation (Vapnik et al. 1997) where w and b are the normal vector and scalar, respectively. They can be derived from the minimized regularization risk function (Vapnik et al. 1997) showed in Eq. 2 as follows: where C is the positive penalty coefficient, which is the tradeoff between model training error and model flatness. e is the insensitivity coefficient, which determines the error tolerance.a i and a Ã i are the Lagrange multipliers related to the constraint.
Þ is the kernel function. After determining a i , a Ã i and b the linear function in Eq. 1 can be expressed explicitly by Eq. 3.
One of the key problems in the application of SVR model is to select an appropriate kernel function. Four kinds of kernel functions are commonly used, namely linear, polynomial, sigmoid and Gaussian (radial basis) (Ç evik et al. 2015). Considering the efficiency and reliability, especially in the face of a variety of parameters, Gaussian kernel function is selected here as shown below: where c is the kernel parameter. The generalization ability of SVR depends on the proper setting of parameters C, c and e. Genetic algorithm is used to determine the optimal value of parameters.

Genetic algorithm (GA)
In the past few years, GA has been applied as an effective tool to solve different optimization problems in engineering and academia (Yan and Lin 2016;Chatterjee et al. 2017; Umeonyiagu and Nwobi-Okoye 2019). GA is not only a global optimization search algorithm by simulating the process of natural selection and biological evolution (Taylor 1994), but also a probabilistic parallel optimization method, which can adjust the search direction adaptively without the need of certain rules. The main characteristic is that it acts directly on the target object and uses the fitness function instead of the cost function without the requirement of function continuity and derivability. The complex problem can be solved by three genetic operations: selection, crossover and mutation.
In the GA, search parameters are binary encoded to produce many binary strings called chromosomes. Multiple chromosomes form an initial population. The aim is to get For this purpose, a fitness function associated with chromosomes is defined. The higher fitness, the higher probability of chromosome selection, and then the selected chromosomes are crossed and mutated to produce a new population. Finally, a chromosome with the best fitness is obtained after finite generations, which is decoded to obtain optimal parameters. From this process, it can be seen that GA is independent of the specific domain of the problem and very robust to various optimization problems. The searching efficiency of GA is affected by some factors, such as population size, crossover probability and mutation probability, but it is out the scope of our paper.

The hybridized SVR-GA model
The main problem of SVR model is to select appropriate kernel function and hyper parameters, which is time consuming work (Ç evik et al. 2015). GA can be employed to quickly select the optimal parameters for SVR model. The flow diagram of proposed hybrid SVR-GA model was shown in Fig. 2. In the SVR-GA model, the population size was 50, while the evolutionary iteration was 200 generations, and the probability of crossover and mutation was 0.8 and 0.09, respectively. The SVR parameters, C and c were binary encoded while e was always 0.01. For better evolution, the fitness function here adopts the correlation coefficient, which is calculated by the experimental value and predicted value of the SVR model. Higher fitness value means higher ranking for a chromosome, while a lower ranking chromosome is less likely to be selected. According to the principles of survival of the fittest, the latest generation chromosome with the best fitness can be decoded as the approximate optimal parameter after generation after generation of selection, crossover, and mutation. In the end, the optimal parameters can be obtained for the SVR model used for data training and model validation, and the well trained SVR-GA model for prediction with high accuracy can be obtained. The key parameters of SVR named C and c can be determined by the pseudo-code combing GA (Deveci and Demirel 2018;Gao et al. 2019) and SVR shown in Algorithm 1 with initial parameters shown in Table 2.

Model implementation and performance metrics
The dataset containing 148 experimental samples was split into training set and testing set for SVR-GA. 16% of the dataset (24 samples) were randomly selected as testing set, which were not involved in the model training. The remaining 84% of the dataset (124 samples) were used to  (Zheng et al. 2019) train and validate the model. To accurately evaluate the prediction ability of the model preventing from over-fitting and under-fitting conditions, a fivefold cross-validation method was utilized. It means the training set was divided into five equally sized subsets. Each fold was used to validate and obtain the optimal hyper-parameters (i.e., C, c) for SVR model by GA, while the other fourfolds were used to train SVR model and get parameters w and b. In this process, one SVR-GA sub-model (f i ) is generated for each fold. Over all folds, the average training error was calculated through five SVR-GA sub-model using all training set. The final prediction model of SVR-GA takes the average of the five sub-model predictions. More details for the fivefold cross-validation method in modeling process were shown in Fig. 3. The experimental dataset was normalized with Eq. 5 before model training and testing. where x n is the normalized value of experimental data (x real ), between -1 and 1. x max and x min are the maximum and minimum values of x real , respectively. After model predictions, the inverse normalization is needed for the predicted data shown in Eq. 6.
x n ¼ 2 x real À x min x max À x min À 1 ð5Þ Here, four statistical metrics are used to measure performance of proposed model, such as coefficient of determination (R 2 ), mean absolute error (MAE), root mean squared error (RMSE), mean square error (MSE). The closer the value of R 2 gets to 1, the better the prediction achieves. The formula of R 2 can be shown as below.
where y i and b y i are the ith observed and predicted values, respectively. À y and À b y are mean values of observed and predicted, respectively. n is the total number of observations. The low values of MAE and RMSE indicate good prediction accuracy of the model. MAE, MSE and RMSE are given by Eqs. 8-10.

Results and discussion
In this section, we employed three ML models including MLR, SVR and also SVR-GA to estimate the shear strength of medium-to ultra-high concrete strength beams with stirrups. The results of the three proposed models are compared with the results of the last edition of ACI 318. Scatter plots with plumb error lines for the results of the proposed models based on training set are shown in Figs. 4, 5, 6 and 7. It is clear from Fig. 4 that SVR-GA has more accuracy than the other two ML models, as well as ACI 318. The box diagrams of the ratio of experimental value to predicted value for both training set and testing set are shown in Fig. 8. Compared to MLR and ACI 318, the upper and lower quartiles and the median of both two SVR and SVR-GA are all close to 1 with small dispersion, indicating that both SVR model and SVR-GA model can well predict the shear strength, and the SVR-GA model performs better than the SVR model for both training set and testing set. Furthermore, we also made comparison between experimental shear strength and predicted value of SVR-GA, SVR, MLR and ACI 318, the scatter plots are shown in Figs. 9 and 10 for the training set and testing set, respectively. It shows that all ML proposed models has better results than that for ACI 318 equations. A summary result of the proposed models are also presented in Table 3, which reveals that SVR-GA achieves the least RMSE, the highest R 2 and minimum MAE in both training set and testing set. Therefore, the SVR-GA model is a useful model for predicting the considered target.
According to the description above, SVR-GA had better results than other models, which is proposed for the considered prediction in this study. All figures indicated that the SVR-GA reflects the best performance, and further it had been applied for the sensitivity analysis to explore the

Performance comparison of various models
Considering that MLR model is the most classical prediction model, it is also used as a comparison model. Performances of the MLR, SVR and SVR-GA are investigated and discussed based on the training set and testing set. Here, these models predict results are compared with ACI 318. The developed models were all applied to learn the relationship between the normalized ultimate shear strength (v u ) and eight different input variables (s, b, h 0 , f c , f t , q and q sv f yv ). In order to visualize the results of the models, the experimental and predicted shear strength values versus the experiment number of training set are presented graphically in scatter diagrams shown in Figs. 4, 5, 6 and 7 for the models of SVR-GA, SVR, MLR and ACI 318, respectively. Plumb lines are drawn to indicate prediction errors between the values of experimental shear strength and predicted shear strength. It can be seen that most of the predicted points are closest to the experimental points with SVR-GA model relative to the other models. It also reveals that the predicted errors in the SVR-GA model are lowest following by SVR model. And the prediction result of ACI 318 is too conservative, and its prediction accuracy is worse than all three ML models above. Meanwhile, the predicted errors of SVR and MLR mainly appear in the region of high shear strength, and they also increase with the increase in shear strength. One possible reason is that there are fewer experimental data points available in this region compared to other regions.
The prediction error analysis box diagram of the ratio of experimental value to predicted value is shown graphically in Fig. 8. As a statistical diagram, the box diagram can provide key information about the location and dispersion of the data, and it can also clearly show the maximum, minimum, median, and upper and lower quartiles of the data. When the ratio of the predicted value to the test value is used as the input data of the box diagram, the closer the data point is to 1 and the lower the dispersion degree gets, the higher the accuracy the model can achieve. As seen in Fig. 8, compared to MLR and ACI 318, the upper and lower quartiles and the median of both two SVR and SVR-GA are close to 1 and the degree of dispersion is very small, especially for the SVR-GA model. It suggests that both SVR model and SVR-GA model can well predict the shear strength, and the SVR-GA model performs better than the SVR model.
The comparison of the experimental shear strength with the predicted value of the SVR-GA, SVR, MLR and ACI 318 is presented graphically in scatter plot shown in Figs. 9 and 10 for the training set and testing set, respectively. The fitting lines of the experimental and predicted value are also shown in the figures. It can be seen that the closer the fitting line is to the perfect line (the included angle with the x-axis is to 45°), the better predicted performance the model works. The angle slope of the fitting line is close to the perfect line in both the training set and testing set for SVR-GA and SVR. This validates the consistency of the models of SVR-GA and SVR. We can also derive that the predicted values of the SVR-GA model have less Fig. 7 The comparison of experimental strength to predicted strength for ACI 318 with training set To further determine the accuracies of the four models mentioned above, different metrics such as R 2 , RMSE, MSE and MAE are measured by Eqs. 7-10. The reported statistical results are shown in Table 3 Fig. 11. These results show that the SVR-GA model can achieve best performance in ultimate shear strength prediction. More details can be seen in Table 3.
It can be derived from the above analysis that the highest precision prediction value can be obtained by using SVR-GA, while results obtained by ACI 318 are worse than all of ML models mentioned above both in training data and testing data. It also demonstrates that SVR-GA model can achieve best prediction effect on the shear strength of medium-to ultra-high-strength concrete beams, which is superior to the traditional SVR model, MLR model and ACI 318. Therefore, we applied SVR-GA model for sensitivity analysis in the following section. Fig. 9 The comparison of experimental strength to predicted strength using the training data

Sensitivity analysis
The sensitivity analysis is implemented for exploring the important degree of each of the input variables using the SVR-GA model. Twentieth quantiles (from minimum to maximum, increased by 5%) of each input variable collated from the experimental dataset are used as a new dataset to calculate the normalized ultimate shear strength. To be more precise, the value of one input varies from minimum to maximum, while all other inputs remain with their average values. According to the statistical probability distribution of input variables, the influence of input variable changes on output results is explored shown in Table 4. For example, the minimum value (0.79) and maximum value (5.0) of a/d, and the median value of other parameters are taken to predict the shear strength as v max (15.7606 MPa) and v min (4.9361 MPa), then the important degree of this variable is (v maxv min )/v max = 0.6868. The Fig. 10 The comparison of experimental strength to predicted strength using the testing data  , and is consistent with the results of the literature that the most important parameters are the shear span-to-depth ratio and concrete compressive strength. However, Fig. 12 reveals that the stirrup spacing has little effect on the shear strength of the slender beam, and its important degree is only 0.1658. It is noteworthy that although the range of stirrup spacing in the experimental dataset is very large ranging from 20 to 1000 mm, the median value is only 151 mm, and most of the data are small. Consequently, a larger database of larger stirrup spacing should be considered in future research to explore the degree of importance of the stirrup spacing.
The sensitivity analysis described above shows that machine learning and artificial intelligence technologies can be helpful during the design of beam shear resistance phase. In addition to accurately predicting the ultimate shear strength, the SVR-GA model can also help create more informative input/output maps of the ultimate shear strength. Especially, the most important variables (a/d, f c , q and q sv f yv ) will be used here to illustrate the input and output maps of the ultimate shear strength. The values of the other variables remain the same as the average value. Six ultimate shear strength maps with the same color range are presented in Fig. 13, which shows the relationship diagrams of f c and a/d, f c and q sv f yv , f c and q, a/d and q sv f yv , a/d and q, q and q sv f yv , respectively. Three-dimensional maps indicate that there is a nonlinear behavior in the input and output relationship, which makes it difficult to detect their relationship only from the experimental data. Except that the increase in a/d will reduce the shear strength to some extent, the increase in other input variables will increase the shear strength to varying degrees. Figure 13a and f shows an interesting phenomenon that the ultimate shear strength is not the maximum, when the two most important factors (a/d and f c ) are both used as inputs. On the contrary, the ultimate shear strength reaches the maximum, when taking the two lower important factors (q and q sv f yv ) as inputs. This reflects from the side that the shear strength calculation of the beam with stirrups is extremely complex, which is the result of the coupling of multiple variables, and the variables cannot be considered separately.
It can be derived from the above analysis that the highest precision prediction value can be obtained by using SVR-GA while results obtained by ACI 318 are worse than all of ML models mentioned above both in training data and testing data. It also demonstrates that SVR-GA model can achieve best prediction effect on the shear strength of medium-to ultra-high-strength concrete beams, which is superior to the traditional SVR model, MLR model and ACI 318.

Conclusions
In this study, machine learning models for prediction goals is investigated to propose an effective model with high accuracy to calculate the shear strength of medium-to ultra-high concrete strength beams with stirrups. For this purpose, three ML models including MLR, SVR and SVR-GA are considered and discussed. All models are trained and tested based on an experimental database with 148 samples gathered from literatures, which consists of eight input variables including stirrup spacing, beam width, shear span-to-depth ratio, effective depth of the beam, concrete compressive strength, concrete tensile strength, longitudinal reinforcement ratio and the product of the stirrup ratio and the stirrup yield strength. These inputs are applied in order to determine the shear strength of the considered medium-to ultra-high concrete beams with stirrups. The performances of the proposed ML models are compared with each other and also with ACI 318. The experimental results indicated that all ML models under fivefold cross-validation can achieve lower errors than ACI 318 equations, and the SVR model with optimal parameters obtained by GA had best performance than the other models, which can be applied to predict the shear strength.
It is also concluded that SVR-GA can be used as a suitable tool to predict the shear strength of ultra-high concrete beams with stirrups with a high level of precision. Furthermore, the sensitivity analysis reveals that all input variables can affect the prediction accuracy of the shear strength, and the most important variables are shear spanto-depth ratio, concrete compressive strength, longitudinal reinforcement ratio and the product of the stirrup ratio and the stirrup yield strength. The proposed SVR-GA model here can be a guide with high prediction accuracy for the shear strength design of ultra-high concrete beams with stirrups. It is expected that the shear strength of medium-to ultra-high strength concrete beams with any combination of design parameters can be predicted accurately providing guides for optimal design. At the same time, there are still some issues that need to be resolved. For example, there is poor interpretability for the proposed SVR-GA model, making it difficult to be well understood by designers for optimal design of medium-to ultra-high strength concrete beams with stirrups. The future works can be done with more experimental data with other more interpretable models.