Shear Strength Estimation of Reinforced Concrete Deep Beams Using a Novel Hybrid Metaheuristic Optimized SVR Models

: This study looks to propose a hybrid soft computing approach that can be used to accurately estimate the shear strength of reinforced concrete (RC) deep beams. Support vector regression (SVR) is integrated with three novel metaheuristic optimization algorithms: African Vultures optimization algorithm (AVOA), particle swarm optimization (PSO), and Harris Hawks optimization (HHO). The proposed models, SVR-AVOA, -PSO, and -HHO, are designed and compared to reference existing models. Multi variables are used and evaluated to model and evaluate the deep beam’s shear strength, and the sensitivity of the selected variables in modeling the shear strength is assessed. The results indicate that the SVR-AVOA outperforms other proposed and existing models for the shear strength prediction. The mean absolute error of SVR-AVOA, SVR-PSO, and SVR-HHO are 43.17 kN, 44.09 kN, and 106.95 kN, respectively. The SVR-AVOA can be used as a soft computing technique to estimate the shear strength of the RC deep beam with a maximum error of ± 3.39%. Furthermore, the sensitivity analysis shows that the deep beam’s key parameters (shear span to depth ratio, web reinforcement’s yield strength, concrete compressive strength, stirrups spacing, and the main longitudinal bars reinforcement ratio) are efﬁciently impacted in the shear strength detection of RC deep beam.


Introduction
In many high-rise reinforced concrete (RC) buildings, as the use of areas is changed from one story to another, some columns in the upper stories are not permitted to reach the foundation.To solve this conflict, transfer girders with a considerable thickness named deep beams are required [1,2].Furthermore, deep beams are used in many other critical structures and play a significant role in delivering heavy loads to the bearing elements [2,3].Deep beams have high flexural stiffness, and the shear diagonal failure is the predominant mechanism as the loads are mainly transferred from their action points to the supports locations through a direct diagonal strut [1][2][3][4].Concrete compressive strength, the provided top/bottom reinforcements, and web reinforcements in terms of amount and spacing all form the shear resistance of these deep beams [1,2,4,5].In literature, plenty of analytical and numerical studies have been focused on the ultimate shear strength assessment of such beams with large depths compared to their spans [2,3,[6][7][8][9].Unavoidable discrepancies were found with the implementation of both analytical/numerical methods compared to the experimental results [4].This study aims to propose a novel soft computing approach that can be used to predict an accurate shear strength of RC deep beams based on a wide range of test results collected from different experimental studies.
Many researchers have used different regression methods to estimate the shear strength of RC deep beams [8,10].Recently, soft computing techniques have been proposed and improved the prediction techniques of shear strength of beams [3,7,[11][12][13][14][15]. A conventional artificial neural network (ANN) was applied to estimate the shear strength of the RC deep beam, and the accuracy of the designed model was high [7].The ultimate shear strength of the RC deep beam was also estimated using ANN and compared to different building codes, and the proposed model provided an accurate prediction of shear capacity [6].ANN, adaptive network-based fuzzy inference system (ANFIS), and group method of data handling (GMDH) approaches were used in predicting the shear strength of RC beam-column joints, and the performance of these models was high at a different range of shear strength [15].A probabilistic model was applied to estimate the shear strength of beams, and the determination of shear strength was shown to be accurate [14].Integrated genetic programming and simulated annealing (GSA) outperformed American concrete institute (ACI) and Canadian standard association (CSA) codes in modeling the shear strength of RC deep beams [2].The shear strength of beams reinforced by fiber was calculated using hybrid support vector regression (SVR) and firefly optimization algorithm (FFA), and the designed model was shown to be robust in shear strength prediction [13].ANN was integrated with the adaptive harmony search optimization (AHS) technique for modeling the shear strength of RC walls, and the proposed model accuracy was high [12].Multivariate adaptive regression splines (MARS) and artificial bee colony (ABC) were also integrated to design a model for predicting the shear strength of RC deep beams, and the performance of MARS-ABC was higher than different building codes in shear strength estimation [16].Generally, parameters of machine learning (ML) models are tuned using metaheuristic algorithms to improve the prediction efficiency of ML models [17][18][19][20][21][22].
Meanwhile, novel optimization algorithms have recently been developed, such as the African Vultures optimization algorithm (AVOA), particle swarm optimization (PSO), and Harris Hawks optimization (HHO).Although these techniques are used in different engineering applications [17][18][19], the AVOA optimization technique is not applied yet in shear strength prediction based on our literature.PSO was integrated with an adaptive neuro-fuzzy inference system (ANFIS) to predict the shear strength of high strength concrete for a slender beam, and the ANFIS-PSO attained the best modeling accuracy over ANFIS -ant colony optimizer (ANFIS-ACO), -differential evolution (ANFIS-DE), and -genetic algorithm (ANFIS-GA) [20].Teaching-learning-based optimization (TLBO), PSO, and HHO were integrated with SVR, and the results of SVR-PSO, SVR-HHO, and SVR-TLBO were robust and can be used to estimate an accurate shear strength prediction of RC shear walls [21].
Based on the above literature, the hybrid SVR models are more robust for predicting the shear strength of RC deep beams [3,16,21,23].This study aims to evaluate a new hybrid technique (SVR-AVOA) in predicting the shear strength of RC deep beams.To benchmark the proposed SVR-AVOA model, the hybrid known models SVR-PSO and SVR-HHO are proposed and compared; in addition, the recent mathematical studies and building codes are compared to the proposed model to assess its accuracy of it in modeling shear strength of RC beam.SVR-AVOA, SVR-PSO, and SVR-HHO are developed using different scenarios of input variables.For this study, 202 datasets, including 19 variables of experimental studies, were collected from literature to design and evaluate the proposed models.The sensitivity analysis of optimum input variables is proposed and evaluated.

Background of Variables Impacts the Shear Strength of RC Deep Beams
Figure 1i presents a real case of deep beam function in load transfer of buildings and variables that impact the shear strength value.Figure 1ii demonstrates the main parameters of the deep beam.Figure 1iii illustrates the failure mode of deep beams.In the figures, V represents the deep beam shear capacity, a is the horizontal distance from the load to the support, and d denotes the deep beam depth.Here, V depends on (1) Concrete quality f c (for the diagonal strut), (2) Main steel yield strength f y (for the main tie), and (3) Web reinforcement (horizontal and vertical).As some variables have a big role in forming the deep beam's shear strength, the main variables considered in this study are the shear span to depth ratio, the main reinforcement, ratio and yield strength, concrete compressive strength, and web reinforcement characteristics.Many researchers investigated the deep beam's shear strength evaluation and predictions [1][2][3][4]24,25].They concluded that the compressive strength of concrete ( f c ), the shear span to the beam's depth ratio (a/d), bottom longitudinal reinforcement ratio (ρ), and web reinforcement ratio (both vertical ρv and horizontal ρh) are the main key parameters.The previous analytical/experimental studies focused on the role of f c in forming V.In the case of beams with smaller a/d ratios, the role of the truss mechanism in transferring loads to the support location diminishes, and the direct diagonal strut is the main transferring load mechanism to the support's location.The effectiveness of such struts has a significant impact on V values [3].Smith and Vantsiotis [26] and Ahmed [27] observed an increase in V of the deep beam with increasing fc of concrete, but the relationship was not linearly proportional.In addition, they observed no improvement in V if f c was above a certain limit.The non-proportional increase in shear strength compared to the increase in concrete compressive strength can be attributed to two reasons.First, the limited contribution of the aggregate interlock mechanism in members with high strength concrete compared to the one with normal strength concrete, as the cracks cross the aggregate particles in high strength concrete and do not go around them as in normal strength concrete.Second, the formed tensile strains perpendicular to the main diagonal strut work on reducing the benefits of using high strengths.Oh and Shin [28] noticed a brittle failure of deep beams with concrete of 74 MPa without any warning, which is different from the failure of other beams with 23 MPa.They also observed a decrease in the rate of increase in the ultimate strength of beams with high-strength concrete.
The inclination angle of the main diagonal strut plays a significant role in determining the concrete efficiency in the diagonal strut.This angle is directly dependent on the shear span to depth ratio a/d.As this angle increases, the forces can go directly inside the diagonal strut to the support.The previous studies [3,29] noticed that the increase in shear strength could be detected by decreasing the a/d ratio.Kim and Park [30] found a trivial impact of this ratio on the shear strength of beams with ratios greater than three, and the contrary was noticed for beams with ratios less than three.Oh and Shin [28] concluded that the ratio of a/d is the governing key parameter in determining the shear strength of a deep beam with high-strength concrete.In addition to resisting the induced tensile force of the main horizontal tie, the main reinforcement bars play an important role in controlling the width enlargement of the diagonal main cracks by dowel action mechanism and enable the aggregate interlock to work more effectively.Many researchers investigated the impact of the main reinforcement ratio on the deep beam's shear strength [3,31,32].They observed a significant increase in the shear strength with increasing the main reinforcement ratio but up to a certain limit.Above a ratio of 1.5%, Ashour et al. [33] noticed a local concrete crushing damage due to compressive stress concentration at the top strut far from the main diagonal strut without enabling the diagonal strut to reach its ultimate resistance.
Web reinforcement has an important role in confining the concrete and delivering the tensile stresses at the main shear diagonal crack to the intact zones around the crack, which consequently increases the deep beam's shear strength [2].Both vertical and horizontal web reinforcement has a key role in resisting shear stresses and limiting the enlargement of the width of the main diagonal crack [3].In deep beams with higher a/d ratios, the contribution of vertical reinforcement is more obvious than the horizontal reinforcements, and the contrary is noticed for beams with a/d less than 1.0.As the deep beam's shear strength is dependent on many parameters, the increase of horizontal and vertical reinforcement above a certain limit does not influence the ultimate shear strength as other parameters may govern the situation without reaching the maximum capacity of the provided web reinforcement.
Table 1 presents the existing models used in this study compared to the developed models.The performance of these models was used significantly in shear strength determination for the RC deep beams of structures.
V CLZ is the shear resisted at the critical loading zone, V ci represents the contribution of aggregate interlock, V s is the shear resisted by web reinforcement and V d is the dowel action in the main longitudinal bars.
where: f c is the compressive strength of concrete, b is the beam width, d is the beam effective depth, Av is the vertical web reinforcement, f yv and f yh are the yield strength of vertical and horizontal web reinforcement respectively, s is the spacing between the vertical web reinforcement, ρ h and ρ v are the ratio of horizontal and vertical web reinforcement respectively, n is the modular ratio, ρ is the ratio of the main longitudinal bars.

Material and Data Collection
The Supplementary Material (Table S1) presents the data collected from the literature.For this study, 202 datasets were collected from the literature.In the current study, 19 input variables were used and divided into two categories, the main (8 variables) and other (11 variables) variables, as presented in Table S2.Here, the main variables were considered based on our literature in the previous section, Section 2; the other variables were considered while the impact on shear strength calculation was high.The direct relationship between each variable and the ultimate shear strength (Vu) of the RC deep beam is presented in Figure S1.Exponential, linear, logarithmic, and power functions were used to estimate the best direct relationship equation between Vu and input variables.Table 2 presents the summary of these functions.From Table 2, it can be observed that the power function has the best correlation with Vu in the case of using a/d and f yv of the main variables.The relationship of ρ v with Vu is exponential.The linear correlation can be detected between Vu and (s and ρ h ).ρ, f y , and f c are correlated with Vu based on logarithmic functions.The best R 2 between Vu and the main variables is 0.39 for the f c variable.These results indicate that the relationship between the main variables and Vu cannot be estimated directly, and a complex relationship may be detected by using all main variables.Similarly, for the other variables, the relationship between Vu and variables varies.The best R 2 is estimated using a number of main bars, R 2 = 0.37.The variation in R 2 indicates the complex relationship between Vu and all variables.Table 2 and Figure S1 show the increase of the resistance with beam effective width and height, number of main bars, and concrete strength up to a certain limit.
The statistical evaluation, range (RA), average (M), standard deviation (SD), kurtosis (KU), and skewness (SK) of the used datasets is presented in Table 3. From the table, the range of datasets varies and will affect the models' performances, so the normalized datasets are used to overcome the range change of variables.The data is normalized between 0 and 1 in this study.In addition, in the proposed models, it is recommended to use the given ranges of input variables.The average and standard deviation values show that the distortion of datasets is high.The kurtosis and skewness values indicate that the distribution of datasets is not normal.Figure 2 presents the histogram and distribution of main variables and Vu.The figure shows positive skewness for whole variables is observed; the distribution is also supported by the presented values in Table 3.   Pal and Deswal [23] and Mozumder et al. [34] proposed the SVR formulas and theory in shear strength to predict RC beams.It was found to be a powerful computation technique for predicting the shear strength of deep beams [3,23].Here, a summary of SVR is presented.SVR is the regression category of support vector machine (SVM), aiming to find a function that represents the relationship between inputs features to forecast the corresponding value when a new input is used.In SVR, "a fixed mapping procedure to map its input to n-dimensional feature space; then nonlinear functions are used to fit the high-dimensional features [35]".Vapnik [35] proposed a loss function to allow the concept of SVM margin to be used for regression solutions.The following equation represents the mathematical formula for the SVRs' approximation function [23,34]: where w and ϕ represent the weight vector and transformation functions, respectively; x and d are the input and output vectors, and b denotes a scalar.In Equation ( 1), the w and b are used to determine the normal and scalar vector, respectively, for the high-dimensional space, which is determined through ϕ(x).Terms 1 2 w and are the standard error and the penalty terms, respectively.The ε-insensitive loss function introduced by Vapnik [34] is commonly used to estimate the Equation ( 1) parameters through the following minimization function [23,34]: where C > 0, which controls the trade-off between the model complexity and the amount up to which deviations larger than ε are tolerated.Equation ( 3) can be transformed to a dual space using Lagrange multipliers solution.This solution can be expressed as follows [23,34]: where L denotes the Lagrangian and α i and α * i represent the Lagrange multiplier.Once Equation ( 4) is used to estimate the parameters of Equation (2), Equation (1) can be rewritten as [23,34]: where nsv represents the number of support vectors (x r , x s ).Here, the solution of this equation depends on the training pattern of Lagrange multipliers, which are only applied to estimate the w and b.Therefore, the kernel function is commonly used to solve the nonlinear regression problems in SVR.The Kernel functions can transform the nonlinear problems into linear problems, as presented in Yaseen et al. [36], which allows the SVR to solve more complex problems.The nonlinear SVR can be expressed as follows: where K is the kernel function; K(x i x) = (ϕ(x i )ϕ(x)).In the current study, the radial basis kernel (RBF) is applied.
It should be mentioned that the SVR is built with statistical theory and based on the minimization of structural risk principle [37].It is a popular method for a small count of data, high dimensional, and non-linear problems.Therefore, SVR is used in many applications [13,38,39] for prediction tasks.Generally, SVR is a type of convex optimization technique to search a local solution within a problem domain [37].The tuning of learning parameters of SVR greatly impacts the evaluation quality.Therefore, finding optimal values of SVR parameters from global searched cost is a difficult task [40].Nature-inspired algorithms are proved to be successful in finding the local best solution from the global one.This work applied and showed PSO, HHO, and AVOA for optimizing SVR parameters with faster convergence capability to provide better prediction accuracy of the SVR model in the deep beam's shear strength prediction.

Optimization Methods
To optimize the best parameters (C, ε, and KernelScale (γ)) of SVR, the PSO, HHO, and AVOA algorithms are used separately in the current work.4.2.1.PSO Kennedy and Eberhart (1995) proposed a metaheuristic population-based optimization algorithm named particle swarm optimization (PSO).The mechanism of the PSO is inspired by the fish schooling and foraging of the flock of birds while exploring an unknown region.Further, PSO is defined as a swarm of particles moving nearby the problem space by the influence of its global (G best ) and local best (P best ) position [41].Therefore, a population search algorithm is used in PSO in a nature-inspired manner to analyze the input data features.The best position of the whole swarm is required to optimize the parameters.These can be performed through nature-inspired behaviors and learning experiences of population particles.PSO was found to be a robust integrated technique with SVR and ANN to model the shear strength of concrete [36,41,42].PSO algorithm may be summarized as follows: • First, it initializes the particle of the swarm, then defines the maximum number of iterations, and finally defines the cost function.

•
After defining the cost function, it evaluates the swarm in order to identify the global and local best.

•
Lastly, it calculates the velocity of each particle and then updates its position using the following equations: where y i denotes the i-th particle, k = the k-th dimension of the particle, coe f 1 and coe f 2 represent the acceleration coefficients, w refers to the inertia weight, rand 1 and rand 2 represent the random coefficients, which are randomly limited between zero and one.More details for PSO theory can be found in [43,44].

HHO
Heaidari and Mirjalili (2019) proposed a gradient-free, population-based optimization technique named Harris hawks optimization (HHO) [45].The main inspiration behind HHO is surprised pounce, i.e., the chasing style and cooperative behavior of Harris' hawks in nature.According to the HHO optimization technique, numerous hawks cooperate to surprise a prey by pouncing it from multiple directions.They display a variety of pursuit patterns based on different scenarios and escaping patterns of the prey.The mechanism of HHO is that it mimics the Harris' hawk behavior in that trace, encircle, flush out, and attack the prey.It has been integrated with SVR and ANN to model the concrete characteristics and other engineering applications [21,[46][47][48][49].The main phases in the attacking of hawks are exploration (phase 1), transferring (phase 2), and exploitation (phase 3).In phase 1, the hawk depends on its position from the prey based on his waiting, seeking, and discovering; this can be expressed as follows: where Y ranm and Y prey indicate the random position for the selected hawk and prey's position, respectively.UL and LL indicate the upper and lower range; r i indicates a random number, having a value between 0 and 1; and Y m = 1/N ∑ N 1 Y i (iter).In phase 2, the prey energy is modeled as E = 2E 0 1 − iter T , where T and E 0 ∈ (−1, 1), and they indicate that energy falls for the prey with their escapes.Thus, the hawk can decide the solution based on the E computation and starting in phase 3 when |E| ≥ 1, and exploiting the neighborhood when |E| < 1.Once starting phase 3, hawks decide to apply a soft or hard besiege.|E| ≥ 0.5 indicates the prey still has enough energy to escape, but maybe some misleading jumps occur in it to fail, so a soft besiege works.On the other hand, in the case of |E| < 0.5, the prey is too fatigued to escape, so hard besiege works.
Here, the HHO is used to optimize the SVR parameters.

AVOA
Abdollahzadeha and Gharehchopogh (2021) recently introduced a metaheuristic algorithm named African vultures' optimization algorithm (AVOA) [17].The inspiration behind the development of the AVOA algorithm is the competing and searching behavior of vultures to acquire a large amount of food.To acquire a large amount of food, these vultures, 'N' (N denotes the population of vulture), were divided into categories based on their fitness to find food and eat.The solution with the highest fitness value is treated as the first-best vulture and the second-best solution as the second-best vulture.The rest of the vultures were trying to approach the best vulture.This is formulated as follows.
Step 1: To determine the best vulture in the group.Fitness of all solutions is determined, and the best solution is selected as the best vulture of the group and other solutions will move towards the best solution using: where the value of K 1 and K 2 lies between 0 and 1 with their sum equal to 1.
Step 2: Starvation rate of vultures.The starvation rate is the rate at which the vultures are satiated or hungry.The satiated rate has a declining trend, and to model, behavior Equation ( 11) is used, t = h sin w 0.5π iter i iter max + cos 0.5π iter i iter max − 1 (12) where F indicates the vultures are satiated.If the value of |F| is greater than 1, vultures search for food in different areas, and the algorithm enters the exploration phase, whereas if the value is less than 1, AVOA enters the exploitation phase, and vultures search for food in the neighborhood.iter i Indicates the iteration number, iter max indicates the total number of iterations, ranm 1 has a random value between 0 and 1.Here, z and h indicate random numbers with values lying between −1 to 1 and −2 to 2, respectively.If the value of z is less than 0, it indicates the vulture is starved, and if it increases to 0, it indicates the vulture is satiated.Here, w indicates the optimization operation disrupts the exploration and operation phases.By increasing the value of w, the probability of entering the exploration phase in the final optimization stages increases, and vice versa for decreasing the value of w.
Step 3: Exploration phase.In this phase, different random areas can be examined using two different strategies.To select the strategies in the ranm 1 exploration phase, a random number between 0 and 1 is generated.This procedure is shown in Equation (13).
where P(i + 1) indicates the vulture position vector in the next iteration, F indicates the rate of vulture being satiated, and R(i) indicates one of the best vultures, which is selected at step 1 in the current iteration.X indicates a coefficient vector that increases the random motion by changing in each iteration and is obtained using the formula X = 2 × ranm, where ranm is a random number between 0 and 1.
Step 4: Exploitation phase.In this phase, the efficiency stage of the algorithm is investigated.If the value of |F| is less than 1, the algorithm enters the exploitation phase.Here, the exploitation phase is categorized based on the |F| value.If |F| is between 1 and 0.5, the rotating flight strategy of the vulture is processed based on parameter P 2 .This can be processed as follows: If |F| is less than 0.5, the two vultures' movements accumulate several types of vultures over the food sources, and the siege and aggressive strife to find food are implemented using parameter P 3 .This can be defined as follows: where, 16) (17) In which, vulture best1 (i) and vulture best2 (i) are the best vulture of the first and second groups, respectively, in the iteration i; d is the problem dimensions, µ and v denotes a random number between 0 and 1, and β = 1.5.
In this work, AVOA is used to tune the SVR parameter set to find the efficient performance of the SVR model.

Models' Development and Accuracy Assessment
This study proposed a new hybrid AVOA metaheuristic algorithm-based SVR model to find the optimal parameters (C, ε, and γ) of SVR and compare its performance with other two metaheuristics nature-inspired algorithms (called PSO and HHO)-based SVR model.Generally, the robustness of the SVR model depends on an appropriate selection of the parameters named as the penalty parameter/"BoxConstraint" (C), insensitive loss function/"epsilon" (ε), and the kernel parameter/"KernelScale" (γ/gamma).The range of these parameters is large, and it is difficult to search for the optimal set of values for these three parameters.Therefore, this optimization task may be solved using optimization algorithms.The authors of this article used three metaheuristic algorithms (AVOA, PSO, and HHO) to find the optimal parameter set of the SVR model.Figure 3 shows the process flow of the proposed technique.Initially, the missing value from the dataset is replaced using the k-nearest neighbor (KNN) imputer method.This study used a Euclidean distance measure to fill the missing value.The whole dataset is bifurcated into a train (80%) and test (20%) set.The three-metaheuristic algorithm is used to train SVR parameters for all/selected featured datasets separately.Root means squared error (RMSE) is selected as the fitness function of all algorithms.Metaheuristic algorithms are sensitive to their different parameter set.The Hit and trail approach is used to select the initial parameter set of metaheuristic algorithms.Table 4 shows the initial value of all parameters of three algorithms for SVR training.Since the number of epochs and population size affect the convergence rate of metaheuristic algorithms, this research work aims to find a faster convergence rate of metaheuristic algorithm with 15 epochs and five population sizes (Table 4).To evaluate the accuracy of the proposed models, eight statistical indices are used: coefficient of determination (R 2 ), variance account factor (VAF), variance inflation factor (VIF), mean absolute error (MAE), root mean square error (RMSE), performance index (PI), mean bias error (MBE), and percentage error (PE).The R 2 and VAF are used to measure the correlation between the measured and predicted values.VIF is used to evaluate the collinearity between the measured and predicted values; VIF > 10 indicates high collinearity.The models' errors are evaluated using MAE, RMSE, and MBE, and the PE is used to estimate the accuracy of the proposed model error in predicting the shear strength of RC deep beams.The mathematical expression of these indices can be expressed as follows: where V i and V pi represent the measured and predicted shear strength, V mean , V max , and V min are the average, maximum, and minimum, respectively, of measured values, adj R 2 is the adjustment R 2 , and N is the number of the data sample.

Sensitivity Analysis
Cosine Amplitude Method (CAM) is used to analyze the strength of the relation between input the parameter and power factor [50].It can also be used to determine the express similarity relation between correlated parameters.To apply CAM, all the data pairs were stated in common X-space.The data pairs used to construct a data array defined K as: Every elements i.e., K i , in the data array K is a vector of lengths j, i.e.,: Therefore, each of the data pairs is represented as a point in m-dimensional space, where each point requires j-coordinates for its complete description.

Results and Discussion
Two scenarios are presented in this section.The first is the evaluation of the proposed models based on all variables and the study of the effect of all variables on Vu estimation.Second, the main variables are considered in modeling and evaluating the sensitivity of input variables on the best selection model.The best solution is compared to existing models in the shear strength determination of RC deep beams.

All Variables Impact on Vu Estimation
Table 5 presents the statistical indices of the proposed models.The numerical investigation of the statistical indices shows that the performance of the AVOA-SVR model is better to estimate Vu with R 2 = 0.98 and RMSE = 32.20 kN in the training stage.The comparison between all models using performance indices of all statistical indices in the training and testing stages shows that the AVOA-SVR model has a high index.However, the performance of HHO-SVR is better in terms of PI, RMSE, and PE in the testing stage.Moreover, Figure 4 illustrates the linear correlation between experimental and predicted Vu values for the proposed models.From Figure 4, it is shown that the performance of the AVOA-SVR model is more accurate than other proposed models in the training and testing stages.The distortion of data points around best fitting is small in modeling Vu with the AVOA-SVR model, and the VIF is higher than in other models, as presented in Table 5.This means that when we used all variables, the AVOA-SVR model can be used to estimate Vu with a model error approach of 6.95%.The comparison between the AVOA-SVR model and previous studies is presented in Table 6.The statistical evaluation of relative predicted shear strength (Vu measured/predicted (Vm/p)) is presented in Table 6; COV is the coefficient of variation.From this table, it can be observed that the AVOA-SVR model performance is better and slightly better than ACI and Russo algorithms, respectively.Although the distortion around the mean for the proposed model is lower than for the Russo technique, the range of datasets for Russo is better than the proposed models.The Liu technique is better than previous and proposed models when considering the whole variables in modeling shear strength through the AVOA-SVR model.The selected variables are used and evaluated in the next section to estimate more accurate Vu values.

Selected Variables Impact on Vu Estimation
Table 7 and Figure 5 present the performance evaluation of the proposed models.The performance of AVOA-SVR is high in the training stage.A high correlation, R 2 = 0.97, and low model error, PE = 2.25%, are observed.In the testing stage, a low distortion around best fitting is observed with the AVOA-SVR.In addition, the statistical correlation factors are high, R 2 = 0.97 and VAF = 94.46, as presented in Table 7 and Figure 5.The VIF values of AVOA-SVR in the training and testing stages are higher than other models.This means the accuracy of AVOA-SVR is acceptable with low distortion around the observed values.Although the PI and RMSE of the HHO-SVR models are lower than for the AVOA-SVR model, the distortion of HHO-SVR datasets is high, as presented in Table 7 and Figure 5.Meanwhile, the performance of the proposed models is shown to be low to estimate the high shear strength (as presented in Figures 4 and 5).However, the performance of AVOA-SVR is seen as more robust.This indicates that AVOA-SVR can overcome the variation change in the data used.Figure 6 also shows a faster convergence rate of the AVOA-SVR model compared to the other two hybrid models.Therefore, the AVOA-SVR can be used to estimate the shear strength of RC deep beams with 3.4% model accuracy.The statistical comparison indices in Tables 5 and 7 show that the performance of proposed models with the selected variables is better than that for using all variables in modeling the proposed techniques.This means that the selected variables are more influential in the shear strength of RC deep beams.

Comparison with Previous Studies and Codes
Table 8 and Figure 7 show the performance of the proposed model compared to the previous studies' formulas.As seen in Table 8, the used selected variable improved the AVOA-SVR performance by 60% in terms of COV.This indicates that the selected variables are significantly affected by the Vu values of the RC deep beams.The comparison between previous formulas and AVOA-SVR shows that the proposed model accuracy is high in estimating the shear strength of RC deep beams.The small range is estimated with AVOA-SVR, and the range is 0.57 kN.The small SD and COV of the statistical indices are observed with AVOA-SVR.This means the accuracy of AVOA-SVR is high compared to other models.The boxplot in Figure 7a shows that the median, red horizontal line of Russo, is close to the true value "1", followed by the AVOA-SVR and Liu models, respectively.The low interquartile range (IQR) value, the height of the box, is observed to be small with the AVOA-SVR model, followed by Liu and Russo models, respectively.The maximum and minimum quartiles, the black horizontal solid lines, are small with the AVOA-SVR model, followed by Liu and Russo models, respectively.The outliers are observed near the median of the AVOA-SVR and far to the median of the ACI model.From the visualization of boxplot results, it can be concluded that the performance of the AVOA-SVR model is more accurate than the previous studies for modeling the shear strength of RC deep beams.In addition, the following model is the Liu model, as this model considers more shear resistance mechanisms and shows a higher normal distribution and lower error than Russo's and ACI's models.The quantile-quantile (Q-Q) plot is presented in Figure 7b for the Liu and AVOA-SVR models for further investigation.The relative shear strength is presented versus the standard normal distribution.From this figure, both models have approximately followed the normal distribution; this indicates that both models can be used to estimate the shear strength of the RC deep beam.The AVOA-SVR model has more correlation with the standard normal distribution, and the VAOA-SVR model is more accurate than the Liu technique in modeling the shear strength of RC deep beams.The scatter plot presented in Figure 7c,d shows that the worst model for estimating the shear strength is the ACI's model.The variation in the best solution "1" is shown as small for Russo, Liu, and AVOA-SVR models, respectively.The most of relative shear strength of the AVOA-SVR model falls within ±20%.The comparison of the AVOA-SVR model and previous models shows that the developed model can be used accurately to model the shear strength of the RC deep beams.Therefore, the AVOA-SVR model is a potential soft computing technique that can be used in predicting the shear strength of RC deep beams.

Sensitivity Analysis of Input Variables
Figure 8 presents the most influential input variables in modeling the shear strength of the RC deep beam.The impact of the input variables is presented for the three models.
From Figure 8, it can be noticed that the significant impact of the ratio of vertical and horizontal web reinforcements is low compared to other variables.The sensitivity of the shear span to depth ratio is high, followed by the yield strength of the main steel, the ratio of the main tensile bars, yield strength of vertical web reinforcement, stirrups spacing, and concrete compressive strength, respectively.The impact of the input variables on output for the other developed models is similar.These results imply that the shear strength of the RC deep beam is highly influenced by the beam geometry, concrete strength, and yield strength of the steel bars.The stirrups spacing also has a large effect on the shear strength of RC deep beams.

Conclusions
This study investigated the use of new metaheuristic optimization algorithms integrated with SVR to model the shear strength of RC deep beams and evaluate the sensitivity of input parameters.SVR-AVOA, -PSO, and -HHO were designed and compared to existing models in the current study.In this study, 202 datasets, including 19 variables of experimental studies, were collected from literature to design and evaluate the proposed models.The common eight parameters (shear span to depth ratio, the ratio of the main tensile bars, yield strength of main bars, concrete compressive strength, the ratio of vertical web reinforcement (stirrups), stirrups spacing, yield strength of vertical web reinforcement, and the ratio of horizontal web reinforcement) are also used to evaluate the performance of the proposed models' in predicting the shear strength of RC deep beams.The performance of SVR-AVOA is high in the cases of the used 19 and 8 parameters for modeling the shear strength.The accuracy of SVR-AVOA is improved by 60%, in COV terms, using the common input variables.Thus, other parameters were found less significant in modeling the shear strength of RC deep beams.The comparison of the SVR-AVOA and the previous studies shows that the accuracy of the proposed model is higher than Liu [4], Russo [8], and ACI [10] by 46%, 63%, and 90%, respectively, in terms of COV.This indicates that SVR-AVOA is the more robust model and can be accurately used in modeling the shear strength of RC deep beams.The sensitivity of the input variables in modeling the shear strength of RC beams with the SVR-AVOA was assessed.This investigation shows the impact of the shear span on the beam's depth ratio, yield strengths of vertical and horizontal web reinforcement, concrete compressive strength, stirrups spacing, and the ratio of the main longitudinal bars on the deep beams' shear strength.
Furthermore, the sensitivity of AVOA algorithm parameters can be tested to balance between exploitation and exploration side for enhancing the SVR performance.In the future, to check the efficiency of the proposed model should be tested on other datasets and other civil engineering application areas.The AVOA algorithm can be combined with other machine learning models like an extreme learning machine, random forest, etc., for prediction tasks.

Figure 1 .
Figure 1.(i) Real case of using deep beams (left) and deep beam terminology (right); V: Shear strength, a: Shear span, d: Effective depth of beam.(ii) Basic reinforcement details of simple RC deep beam.(iii) Failure pattern in deep beam (left), Different mechanism components (middle), and Flow of forces in deep beam (right); where, C = compression force in concrete, T = tensile force in the main steel, Hs = tensile force in horizontal web reinforcement, Vs = tensile force in vertical web reinforcement; 1: Main diagonal strut, 2: Splitting tension force, 3: Main tensile force.

Figure 2 .
Figure 2. Distribution histogram of input and output variables.

Figure 3 .
Figure 3. Process flowchart of the proposed method.

Figure 4 .
Figure 4. Scatter plot of model's performances in the training (upper row) and testing (lower row) stages.

Figure 5 .
Figure 5. Scatter plot of model's performances in the training (upper row) and testing (lower row) stages.

Figure 7 .
Figure 7.Comparison of model's performances, (a) boxplot, (b) Q-Q plot, (c) scatter plot of relative shear strength with measured shear strength, and (d) zoom in for upper plot with +20% limits for the best models.

Figure 8 .
Figure 8. Sensitivity of input variable on model's prediction.

Table 1 .
Summary of existing models.

Table 2 .
Direct relationship functions between Vu and input variables.

Table 3 .
Statistical analysis of input and output variables.

Table 5 .
Statistical evaluation of the proposed models.

Table 6 .
Vm/p statistical evaluation for AVOA-SVR and previous studies.

Table 7 .
Statistical evaluation of the proposed models.
Figure 6.Convergence rate of three models.

Table 8 .
Vm/p statistical evaluation for AVOA-SVR and previous studies.