#### Abstract

The shear strength prediction of fiber-reinforced polymer- (FRP-) reinforced concrete beams is one of the most complicated issues in structural engineering applications. Developing accurate and reliable prediction models is necessary and cost saving. This paper proposes three new prediction models, utilizing artificial neural networks (ANNs) and gene expression programming (GEP), as a recently developed artificial intelligent techniques, and nonlinear regression analysis (NLR) as a conventional technique. For this purpose, a large database including 269 shear test results of FRP-reinforced concrete members was collected from the literature. The performance of the proposed models is compared with a large number of available codes and previously proposed equations. The comparative statistical analysis confirmed that the ANNs, GEP, and NLR models, in sequence, showed excellent performance, great efficiency, and high level of accuracy over all other existing models. The ANNs model, and to a lower level the GEP model, showed the superiority in accuracy and efficiency, while the NLR model showed that it is simple, rational, and yet accurate. Additionally, the parametric study indicated that the ANNs model defines accurately the interaction of all parameters on shear capacity prediction and have a great ability to predict the actual response of each parameter in spite of its complexity and fluctuation nature.

#### 1. Introduction

Although the behavior of steel**-**reinforced concrete (RC) members in shear has been extensively studied for more than one hundred years, and after numerous research challenges on shear behavior and on identification of the complex shear resistance mechanism of RC members, still there is an absence of comprehension and agreement on the mechanisms and behavior of shear resistance. While the behavior of steel**-**RC members in shear is an area of concern, the behavior of fiber-reinforced polymer**-**reinforced concrete (FRP**-**RC) members has additional complications due to their different mechanical properties [1]. Many important issues related to shear problems of FRP**-**RC members still remain without in**-**depth investigations and become open to discussion. Moreover, past investigations [2–4] inferred that the present shear design equations are exceptionally conservative in predicting the shear capacity of FRP**-**RC beams. Consequently, the extra amount of FRP used to resist shear could be both costly and likely to create reinforcement congestion problems, and therefore, this study was motivated**.**

The shear strength of RC beams without web reinforcement was observed to be influenced by several variables, such as concrete strength (), shear span-to-depth ratio (), beam depth (), longitudinal reinforcement ratio (), and beam width () [5, 6]. Additionally, for FRP**-**RC members, the type of FRP rebar and their variable mechanical properties, like low elastic modulus (), high tensile strength, and low transverse shear strength, are other factors that should be taken into account.

Albeit many theoretical and empirical studies have been conducted to predict the shear strength of FRP-RC members, and to investigate the interaction between influencing parameters on shear strength mechanisms, up to date there is no general agreement on a specific model. The vast majorities of the existing shear design equations have different forms and do not give a reliable factor of safety against shear failure; additionally, there is no general agreement between the existing codes and equations on the parameters affecting the shear capacity. This controversy in predicting models may be attributed to two main reasons; firstly, there is no comprehension theory on shear failure due to its complicated behavior and failure type which is sudden and catastrophic, and secondly, the existence wide variability in mechanical properties of FRP bars. Therefore, some existing predicting models were developed based on modifying the different theories which is originally proposed for steel-RC members and others are empirical or semiempirical models which were proposed based on a limited data; for instance, the ACI 440 algorithm was calibrated based on the test results from 326 steel-RC specimens and 44 FRP-RC specimens and the latter having a maximum effective depth of 360 mm. An alternative to classical and conventional methods for the prediction of complicated problems in various disciplines, such as the shear strength of FRP-RC members, is the utilization of artificial intelligence techniques like artificial neural networks (ANNs), fuzzy inference system (FIS), and genetic programming (GP). Throughout the most recent two decades, these methodologies have been prominent and effectively utilized as a part in an extensive variety of scientific applications particularly in civil engineering disciplines.

Recently, several studies have been performed on using different artificial intelligence techniques for predicting the shear capacity of FRP-RC members without stirrups. Kara [7] utilized gene expression programming (GEP) to obtain a simple model based on a set of 104 databases. Bashir and Ashour [8] suggested a model based on ANN using a set of 128 databases. Nasrollahzadeh and Basiri [9] developed an FIS to predict a model depending on a 128 database. Lee and Lee [10] used a database of 110 samples to propose a theoretical model utilizing ANN. Golafshani and Ashour [11] introduced a model using a new technique of biogeography-based programming based on an experimental database of 138 test specimens. Most of these studies proposed good models with a high level of accuracy; however, they used relatively a small size of a database. The size of the database used for training and testing the models plays a great role in the success of proposed models.

The present research aims to propose three new models to predict the shear strength, , of FRP**-**RC slender beams without stirrups based on utilizing ANNs and GEP as artificial intelligence techniques and nonlinear regression analysis (NLR) as a conventional technique. These techniques were chosen due to the fact that the ANNs and GEP are powerful, accurate, and highly efficient tools in the model prediction of complicated problems, while the NLR is a classical tool which can be used to derive a simple, rational, and yet accurate model. The main affected parameters on shear strength are considered in modeling process which are , , , , , and . The three models are derived using the largest collected database from the literature of (269) test specimens. The results of obtained models are compared to a large number of codes and proposed equations from the literature to examine their accuracy, validity, and efficiency. Additionally, a parametric study was conducted to indicate and compare the effect of all parameters on shear capacity predicted by the proposed models and different well-known codes.

#### 2. Reviews of Available Shear Design Equations

The last two decades witnessed development and rapid increase of using FRP-reinforcement for concrete structures; therefore, there are universal endeavors to create prediction models. These endeavors have brought about distributing a few codes and outline rules. However, there is no agreement in regards to shear configuration models among all available codes and guidelines. Table 1 provides a summary of the shear design equations for FRP**-**RC members without stirrups.

Most of the available shear design equations of FRP**-**RC members have been developed by modifying the existing equations adopted by design codes for steel**-**RC members. The methods of modification tried to cover the major changes in material properties by taking into account the effect of the lower elastic modulus, higher tensile strength, lower transverse shear strength, and no yielding criteria in the FRP rebar. Nevertheless, the relatively large variations in mechanical properties such as tensile strength, elastic modulus, bond properties, and transverse shear strength between the different FRP types have made it hard to precisely predict the shear capacity. Furthermore, the procedure that adopted by various codes, manuals, and guidelines to predict the shear strength are significantly different, and most of them are empirical and based on a limited set of test results. Consequently, some of the design models are quite conservative, whereas others sometimes yield unconservative results [3, 22]; as well as these equations are differing from each other in considered parameters. On the contrary, the design model that does not consider the effect of known parameters would lack generality, and its applicability to general design situations would be dubious [23]. Therefore, there is a need to develop shear design equations which will consider the effect of all parameters; as well as it will be accurate, consistent, and simple to use for the general application.

#### 3. Collected Databases

For the purpose of this study, an extensive survey of the open literature was conducted and a large database containing the test results of a large number of 269 FRP**-**RC beams and one**-**way slabs were assembled. The database of the specimens were collected from 42 different references, reinforced with FRP rebars without transverse reinforcement, and only those that failed in the shear mode was compiled [1–3, 20, 22, 24–60]. Only slender beams and one way slabs with ≥ 2.5 were considered in this study. The specimens included 230 beams and 39 one way slabs; all were simply supported and were tested either under three points or four bending points. The details of the specimens that were used as a database are shown in Table 2.

#### 4. Artificial Neural Networks (ANNs)

##### 4.1. General Background

In its most broad frame, a neural network is a machine that is intended to demonstrate the manner by which the brain plays out a specific errand or function; the system is generally executed by utilizing electronic parts or is mimicked in programming on an advanced PC. The ANNs perform helpful calculations and simulate complex modeling through a procedure of learning similar to that of the human brain. ANNs do not represent the sophistication, complication, and multifaceted nature of the brain, as the simulated neurons are substantially more straightforward than natural neurons in the area of recognition, efficiency, control, and learning. Nonetheless, ANNs are similar to the brain in two main points; first, the building squares of the two systems are straightforward computational tools which are exceedingly interconnected; second, the connections between neurons decide the task of the networks [61, 62].

ANNs are made out of numerous interconnected processing units operating in parallel. Each processing unit keeps some data locally, can play out some basic calculations, can have numerous input data and, however, can send just a single output. The ANNs have the ability to react and create the relating response and to adjust to the changing condition by learning from experience [63].

Generally, ANNs are trained, with the goal that specific input data prompt a particular target output. The network is balanced, in view of an examination of the output and the target, until the networks yield the objective. Gradually training the network modifies the weights and biases required after the introduction of every individual input vector. ANNs have been prepared to perform complex functions in different fields of application including biological, engineering, business, financial, manufacturing, medical, and military. Today ANNs can be prepared to tackle issues that are troublesome for regular PCs or human beings [64]. The drawback of the ANNs modeling is that it is not as simple as ordinary design provision, where it is consist of regularly numerous small units; therefore, simple calculations of the ANNs models are excessively monotonous, and usually, computers are needed for calculations.

##### 4.2. Back Propagation ANNs

The multilayer perceptron, trained through the back-propagation algorithm, is nowadays the most broadly utilized neural networks. The term back propagation neural networks (BPNNs) alludes to the way in which the gradient is processed for nonlinear multilayer networks [64]. The BPNNs algorithm follows a technique to minimize the errors through the calibration of the weights in each cycle by a small amount for a specific training scheme.

BPNNs is one of the straightforward and most appropriate networks being utilized in the prediction of civil engineering problems [63, 65, 66], principally because of that it can modify the weights of each layer in view of the errors introduced at the output. A common structure of BPNNs display comprises of an input layer, at least one hidden layer and an output layer, and each layer comprises of many neurons. Every neuron is a processing unit that gets at least one source of input and produces an output response through the transfer function. Related to every connection is a weighting that communicates the impact on the present procedure component of an information set or another procedure component in the past layer, as shown in Figure 1. The connection weighting and bias values are at first picked as arbitrary numbers and afterward settled by the consequences of the training procedure [67].

The calculations that carried out inside the BPNNs include summation process of weighted input values with the bias of the input layer, as in Equation (1), and then the results passed into the hidden layer neurons. Inside neurons, the process of sigmoid transformed function was performed by using one of the transformed functions such as that in Equation (2). Another linear process, using the purelin function, was carried out by summing the weighted output from the hidden layer neurons with the bias of the hidden layer as in Equation (3):where is the weighted sum generated at the *j*th hidden neuron; is the input value from the *i*th input neuron; and are the weights added to the hidden layer and the output layer neurons, respectively; and are the biases added to the hidden layer and the output layer neurons, respectively; is the processed output from the *j*th hidden neuron; is the processed output from the *k*th output neuron; and is the number of input neurons, and is the number of neurons in the hidden layer.

##### 4.3. BPNNs Model Used for the Prediction of Shear Capacity

In this study, the BPNNs were utilized to produce a predictive model for shear strength of FRP-reinforced concrete members without stirrups. The key parameters affecting the shear strength which considered in this study are , , , , , and . The most important issue in BPNNs modeling is selecting the number of hidden layers and the number of neurons in each hidden layer, which depend on several factors like type of problem of interest, the association of input parameters with the output pattern, size of instances, number of input parameters, and number of outputs. However, there is no basic acknowledged criteria to pick the number of hidden layers and the number of neurons; therefore, a trial and error method was achieved based on training and testing up to obtain a reasonable fit. Consequently, a BPNNs training algorithm was utilized in a feed-forward network trained using nonlinear hyperbolic log-sigmoid transfer function in the hidden layer and a linear purelin transfer function in the output layer, as in Equations (1)–(3).

The trial and error procedure resulted in that the best fitting model was that obtained using thirteen neurons in a single hidden layer (6-13-1), as demonstrated in Figure 2. The predictive model was trained through numerous iterations to find the best epoch value, momentum rate, and learning rate which give the optimum model for shear strength prediction. The data were divided arbitrarily into “training” 70%, “validating” 15%, and “testing” 15%, and the partitioning of data was chosen based on the results of model performance after numerous trials.

It is recommended by many authors [11, 63, 68, 69] to normalize the input and output data before using them to predict the model. In this study, the data of input and target were scaled by the linear normalization function in the range of (0 to 1), as shown in the following equation:where is the data sample, is the normalized data sample, and is the minimum and maximum values of the data for the interested parameter.

The input layer weights (ILWs), the input layer biases (ILBs), the hidden layer weights (HLWs), and the hidden layer bias (HLB) of the (6-13-1) predictive BPNNs model are shown below:

#### 5. Genetic Programming

##### 5.1. General Background

Genetic programming (GP) is a naturally inspired machine learning (evolutionary) technique utilized for randomly reproducing a population of computer programs in light of Darwin’s advancement hypothesis. GP, which was developed by Koza [70], can be regarded as a developed form of the genetic algorithm (GA), which is an evolutionary optimization tool that depends on the principles of genetics and normal determination. The systems of GA incorporate altering and changing the extent of chromosomes by genetic operators to solve function optimization problems. The fundamental contrast amongst GA and GP is that the solution of GA is a binary string of fixed length utilized for parameters optimization of a specific arrangement of model parameters, while the solution of GP is a development program represented as a subset of parse trees which use input variables and corresponding outputs for generating optimized models [71–73].

The GP techniques for modeling incorporate production of an initial population of chromosomes (models) which include an arrangement of functions and an arrangement of terminals that selected and arranged randomly in the form of the parse tree (computer model). The models composed of a root node, functional nodes, and terminals, as shown in Figure 3.

At that point, the produced population models are assessed and the fitness of the models is evaluated using the provided data of the particular problem. After that, other generations of models are produced by primary genetic operators: reproduction, crossover, and mutation, as shown in Figure 4. The evolutionary procedure of assessing the fitness of the present population and producing new population is continued until the point when an end termination criterion is approached, which can be either a particular acceptable error or a specified maximum number of generations. The model with the best solution which developed after each generation assigns the outcome of the GP program [76].

##### 5.2. Gene Expression Programming

The most widely used GP in civil engineering problems is gene expression programming (GEP), which was developed by Ferreira [71]. GEP utilizes development of numerical conditions that is encoded linearly in chromosomes of a consistent length and considered nonlinearly as expression trees (ETs) with various shapes and sizes. The chromosomes are composed of several genes; each gene contains a smaller subprogram or subexpression trees (Sub-ETs). Every gene has a consistent length and comprises a head and a tail. The schematic of the GEP flowchart is shown in Figure 5. The implementation of GEP predictive modeling incorporates five main components including function set, terminal set, fitness function, control parameters, and termination condition [7, 75].

The preferred quality of GEP is that the genetic operators work at the chromosome level and the production of genetic diversity is to a great degree simplified. Another favorable position is credited to the remarkable multigenetic nature of GEP, which permits the advancement of all the more capable models made out of a few subprograms [71].

##### 5.3. GEP Model for Predicting Shear Capacity

For predicting the shear strength of FRP-reinforced concrete beams without stirrups by GEP-based formulation, six influencing parameters were considered: , , , , , and , through training and testing the experimental results. The data are evaluated on the GeneXproTools 5 [78] to develop the empirical model which is the most widely used software for GEP model prediction. Similar to the ANNs model, the data are partitioned randomly to 70% for training and remaining 30% for testing such that both portions of data have balanced statistical parameters like mean, standard deviation, and max and min values.

To obtain the best GEP prediction model, many models with a different number of genes were evolved throughout a set of genetic operators (mutation, transposition, and crossover). Initially, a model composed of two genes with multiplication linking functions and head sizes of three was selected and evolved many times. Then, the model parameters were changed, step by step, through increasing the number of genes, head size, number of chromosomes, and weights of function sets. The program was evolved numerous times for various models, and the outcome models were monitored and compared to evaluate their performance. Other parameters like mutation rate, inversion, and points of recombination were chosen based on previous works [7, 79–82] and then evaluated to get their optimum effects. After many trials, the final mathematical model was chosen, for which the utilized parameters, fitness function, and function sets are given in Table 3. The final model was selected based on criteria of the best fitness and less complicity of mathematical formulation; the expression trees are shown in Figure 6 and the mathematical formulation is given in the following equation:where: refer to , , , , , and respectively, = 6.56723946184092, = −241413.053630399, = −12.1479701959083, = 3.07014791103389, and = 22.0241396443835.

**(a)**

**(b)**

**(c)**

**(d)**

#### 6. Nonlinear Regression Model (NLR)

Regression analysis is considered as an important tool that can be used in the modeling process for solving complicated problems of various engineering disciplines. However, the proposed models based on nonlinear regression (NLR) are less accurate than the ANN and GEP; however, it is characterized by that it is simple and easy for application without the need of computers for calculations. A nonlinear regression model based on the experimental database collected from the literature was suggested. The proposed equation was derived based on developing a model similar to that proposed by Zsutty [83] which is the most accurate and simple empirical equation that was proposed for shear prediction of steel-RC members. In order to predict the shear strength of FRP-RC beams more accurately, considering the effect of all influencing parameters, the experimental shear strength was analyzed with most common influencing parameters: , , , , , and . Based on data fitting analysis, and after some simplification, the following equation was proposed:

It can be observed that the equation is simple, rational, and more generalized than most codes and proposed equations from the literature as it considers the effect of all the influencing parameters on the shear strength of FRP-RC members without stirrups.

#### 7. Results and Discussion

The shear strength results of FRP-reinforced concrete members for the three proposed models, the ANNs, GEP, and NLR, are evaluated and compared with the common design provisions (ACI 440-15, JSCE-97, CSA 806-12) and seven proposed equations [1, 4, 7, 10, 19–21]. The comparison was examined based on statistical terms such as mean, standard deviation (SD), and coefficient of variation (COV) in addition to error values, such as mean absolute error (MAE), root mean square error (RMSE), and , as shown below:where and are the experimental and predicted shear force in kN, respectively, is the average of the experimental values, and is the total number of the considered samples. The model that minimized the error values (MAE and RMSE) and maximized is selected as the optimum model.

Figure 7 presents the scatter of the experimental shear force, *V*_{exp}, versus the predicted shear strength, *V*_{pred}, for the proposed models, codes, and proposed equations from the literature. It is clear that the best models for predicting the shear strength are those suggested in this study, which in a sequence are ANNs, GEP, and NLR, as these models predict the shear strength very accurately with the optimum statistical performance of , RMSE, and MAE. High correlation coefficients with low error rates indicate that the proposed three models are excellent models and exhibit generalization performance in predicting shear strength of FRP-reinforced concrete beams without stirrups. In contrast, the code design provisions and other proposed equations from the literature yielded dispersion results with lower and larger error rates as given in the scatter relationships.

Table 4 presents the statistical parameters (mean, standard deviation (SD), coefficient of variation (COV), minimum, maximum, range, RMSE, MAE, and ) of the ratio of experimental shear strength to the predicted shear strength, *V*_{exp}/*V*_{pred}, for the proposed models in this study, seven well-known codes, and seven proposed equations from literature based on the (269) collected database. It is clear that the proposed equations yielded superior accuracy between all the codes and proposed equation from the literature because their optimum statistical parameters such as the mean, SD, range, RMSE, and are 0.998, 0.113, 0.72, 5.3, and 0.9877, respectively, for ANNs; 1.007, 0.162, 0.88, 8.2, and 0.970, respectively, for GEP; and 1.003, 0.174, 0.99, 10.7, and 0.95, respectively, for NLR.

The ANNs and GEP models are artificial intelligent models derived based on soft computing; as verified above, these models are accurate, superior, and highly efficient in predicting the complicated problem such as shear strength. However, the drawback of artificial intelligent models is that it cannot be calculated manually and it needs computers for calculation. Besides, it is considered as a powerful and efficient tool which can be used in soft programs specialized for analysis and design of reinforced concrete structures. On the contrary, the NLR model is simply accurate and rational which can be used manually with sufficient accuracy more than the design codes and previously proposed equations. Nevertheless, the regression equations cannot properly define, with the same level of efficiency of artificial intelligent techniques, the interactions between the input and output parameters due to that the regression models usually are produced based only on a few predefined equations [79]. Finally, it can be said that all the three proposed equations can be considered as accurate, efficient, and useful models, which were trained on a largest collected database, for predicting the shear strength of FRP-reinforced concrete members.

#### 8. Parametric Analyses

For further verification of the proposed ANNs, GEP, and NLR prediction models, a parametric analysis was performed in this study. The main goal was to find the response of the three proposed models to predict the shear strength for a set of input parameters which are selected taking into account the average value of each parameter in the collected database. For further showing the performance of the proposed models and for comparison purpose, the response of three well-known codes (ACI.1R 440-15, JSCE-97, and CSA S806-12) was also considered. The approach depends on keeping all input parameters at their selected values, except that, in each time, a single parameter value was changed incrementally within the expected practical values. Consequently, an arrangement of manufactured data for the single considered parameter is produced; this procedure was repeated for all parameters so that all models are examined for the entire input parameters. The selected input parameters are = 46 MPa, = 200 mm, = 400 mm, = 4, = 0.011, and = 74 GPa. It is worthy to note that the ANNs model cannot properly predict the shear strength for input parameters outside the range of training data; therefore, it is only valid for prediction within the range of input parameters of the training data [8].

Figure 8 indicates the tendency of the predicted shear strength by different models to the varieties of the main design parameters (, , , , , and ). Generally, it can be observed that the proposed models, particularly the ANNs model, yielded larger shear capacities for all parameters, while ACI 440-15 yielded the lowest shear capacities which indicate that the ACI 440-15 is very conservative.

The relationship of beam depth with the predicted shear strength showed that all the proposed models and codes increased with the increase in the beam depth throughout the range of 100–1000 mm. The ACI 440-15 considers the effect of beam depth linearly, while all other models represent the effect of depth nonlinearly; the ANNs, GEP, and NLR models consider the depth raised to a power of 0.83, 0.67, and 0.66 respectively. Also, it can be observed that the GEP and NLR yield almost the same results. All the proposed models and codes showed an increase in shear strength with the increase in the concrete compressive strength but in different rates. It was found that is raised to a power of 0.66, 0.30, and 0.2 for ANNs, GEP, and NLR, respectively.

The effect of the flexure reinforcement ratio on shear strength indicates a similar trend of the NLR model and the design codes, but the ANNs and GEP models showed a slightly decreased trend for reinforcement ratios greater than 0.03; this is not in agreement with codes and previous studies because of that only few data were available with larger than 0.03, and these data actually showed decreasing shear capacity. Additionally, these models were derived based only on the artificial technology of data training without manually choosing the general mathematical shape of the model, in contrary to the NLR model and design codes which was derived based on the theoretical or semiempirical criteria.

The shear strength also increased with the increase in the elastic modulus of the flexure reinforcement ratio. It can be observed that all models showed a similar trend line for the entire range of data, while the ANNs model showed an increased rate up to about 110 GPa, and then the increase rate reduced significantly; similar results also obtained by Bashir and Ashour [8]. On the contrary, the trend lines of the shear span to depth ratio are decreased for GEP and NLR models for the entire range of data; meanwhile for ANNs model, it decreases up to of 5, and then it started to increase slightly. This behavior is in agreement with the few available experimental results [32, 52, 56] which were studied as greater than 5. However, other codes do not consider the effect of as the trend lines for ACI 440-15 and JSCE-97 are constant, while the CSA S806 trend line remains constant beyond the value of 5. Finally, the effect of beam width on shear strength is linear as indicated by trend lines for all models except the ANNs model which showed a lower increase rate beyond 700 mm width.

To investigate how the proposed models consider the size effect on shear capacity, the same parametric data used for studying the beam depth are used. Figure 9 presents the relationship between normalized shear strength as () and beam depth. It is clear that the normalized shear strength is decreased with the increase in the beam depth for ANNs, GEP, and NLR showing that the size effect is efficiently presented in these models. The trends of GEP and NLR are very close and compatible and decreased from a depth of 100 mm to 1000 mm, while the ANNs model seems to be more accurate as it shows an increase in the trend up to the depth of 200 mm, and then it decreased similarly to that of GEP and NLR; after that, the trend starts to increase slightly.

It can be summed up to that the artificial intelligent models, due to their efficiency and powerful ability in recognizing the parameters interaction, demonstrate various rates of change and fluctuation for the entire scope of input parameters. The ANNs and GEP models define accurately the interaction of each parameter on shear strength prediction and have a great ability to represent the actual response of each parameter in spite of its complexity and fluctuation nature. Furthermore, these models are evolved by training numerous preliminary linear and nonlinear models without any prior assumption regarding the shape and structure of the mathematical model.

#### 9. Conclusions

An experimental database including 269 shear test results of FRP-reinforced concrete members without stirrups available in the literature was established to propose three prediction models to estimate the shear capacity of FRP-RC members using three powerful tools which are ANNs, GEP, and NLR. The following conclusions can be drawn from this study:(i)The proposed models for predicting the shear capacity of FRP-RC members based on ANNs and GEP techniques showed excellent performance, great efficiency, and high level of accuracy and yielded consistent results throughout a wide range of influencing parameters; moreover, the ANNs model is more superior than the GEP, as the SD and for ANNs are 0.113 and 0.9877 and for GEP are 0.162 and 0.97, respectively.(ii)A simple, rational, and yet accurate empirical equation was proposed for predicting shear strength of FRP-RC members by NLR. The equation considers all the important influencing parameters and simplified to be easy for use. In addition, it gives consistent and accurate results with SD and of 0.173 and 0.95, respectively.(iii)Through a comparative statistical analysis for the proposed models with seven shear design codes and six proposed equations from the literature, it has been found that the most accurate and efficient models are ANNs and GEP and to a lower dgree NLR. Among the codes and available proposals, the CSA S806-12 and Ashour and Kara [20] equations were appearing to be more accurate than the others.(iv)From a parametric study considering all the influencing parameters on shear capacity, it was found that the ANNs model, and to a lower level GEP model, defines accurately the interaction of each parameter on shear capacity prediction and has a great ability to represent the actual response of each parameter in spite of its complexity and fluctuation nature.(v)The ANNs and GEP models demonstrate that the effect of flexure reinforcement ratio on shear capacity decreased beyond 0.03, as a real response to the available few data in this range; therefore, more studies are needed to confirm this case.(vi)The ANNs model demonstrated that the effect of does not absolutely reduce the shear capacity, while beyond the value of 5, its effect could become reverse and contribute to increase in the shear capacity as confirmed by few previous studies.(vii)The size effect is considered properly by the proposed models, particularly the ANNs model, which yields more accurate shear capacity prediction to the experimental results.

#### Data Availability

The (database of shear strength of beams without stirrups) data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.