Evolutionary Modeling to Evaluate the Shear Behavior of Circular Reinforced Concrete Columns

Despite their frequent occurrence in practice, only limited studies on the shear behavior of reinforced concrete (RC) circular members are available in the literature. Such studies are based on poor assumptions about the physical model, often resulting in being too conservative, as well as technical codes that essentially propose empirical conversion rules. On this topic in this paper, an evolutionary approach named EPR is used to create a structured polynomial model for predicting the shear strength of circular sections. The adopted technique is an evolutionary data mining methodology that generates a transparent and structured representation of the behavior of a system directly from experimental data. In this study experimental data of 61 RC circular columns, as reported in the technical literature, are used to develop the EPR models. As final result, physically consistent shear strength models for circular columns are obtained, to be used in different design situations. The proposed formulations are compared withmodels available from building codes and literature expressions, showing that EPR technique is capable of capturing and predicting the shear behavior of RC circular elements with very high accuracy. A parametric study is also carried out to evaluate the physical consistency of the proposed models.


Introduction
It is well known that columns are the most vulnerable elements in reinforced concrete (RC) structures.These elements are principally designed to bear axial loads, but as a result of lateral loads, for example, wind pressure or earthquake ground motion, they could deal with relevant shear loads and thus should inevitably be designed to avoid possible shear failures [1][2][3][4][5][6].Composite and steel columns [7] are particularly suitable at this aim.
Circular columns are very popular for bridge pier design, due to simplicity of construction and because their strength characteristics under wind and seismic loads are similar in any direction.Circular elements are also used extensively as columns in buildings, or as piles for foundations, or as secant piling to form diaphragm walls.Despite their frequent occurrence in practice, only limited researches on the shear behavior of RC circular members have been carried out.
Most of the researchers and codes state that the shear strength  Rd of a beam is the sum of the contributions of concrete   and shear reinforcement   , if present.The first takes into account the shear stress transferred by the compressed zone of the beam, the dowel action, the aggregate interlock, and the arch effect and depends on many factors, such as the tensile longitudinal steel ratio, concrete grade, aggregate size, or shear span, among others.The evaluation of this term is debatable and originates on empirical methods.Traditional codes, therefore, propose empirical formulas derived from various tests developed in rectangular sections and then extrapolating results for arbitrary sections, such as circular ones.Also the evaluation of the contribution of the hoops obtained from the truss analogy for rectangular sections was later extended to arbitrary sections, but no specific address was developed for circular ones.
A very limited number of shear models for circular section members exist in the literature.Ghee et al. [8] proposed a model for columns subjected to cyclic loading, consisting of a purely empirical concrete and a classical shear reinforcement contribution, where shear crack inclination 2 Advances in Civil Engineering was obtained by using a lower bound plasticity model.The concrete contribution for low flexural ductility was defined as the "initial concrete strength" and calculated as strength at maximum lateral load; the transverse reinforcement capacity was calculated from 45-degree truss mechanisms assuming that transverse reinforcement exposed by a presumed 45degree diagonal crack is at yield.At higher flexural ductility, a degradation of the concrete term was assumed and the "final concrete strength" was proposed.Different adjustments to this basic model have been added over the years by various authors.Priestley et al. [9] modified the model proposed by Ghee et al. [8] and separated the shear capacity term of the axial load arch mechanism   , from the concrete term   .Other modifications to Ghee et al. original model [8] were developed by Wong et al. [10] and by Kowalsky and Priestley [11].In this work, the authors revised the model proposed by Priestley et al. [9] and the effect of the aspect ratio and longitudinal reinforcement ratio was incorporated.
In addition, the truss mechanism component developed by Ghee et al. [8] has been slightly modified assuming that, in the member compression zone, the cracks are by definition closed.
A different approach was due to Collins et al. [12] by means of the application of the Modified Compression Field Theory (MCFT) [13] to deal with shear strength prediction of circular members.Moreover, this type of approach is more complex than those usually based on strut and tie models, which are still used in many applications, till their first implementation in early 1900, thanks to their simplicity.The practical application of the MCFT requires its implementation as a calculation algorithm within a computer routine, thus strongly increasing the complexity level, also for a standard code formulation.
Bentz et al. [14] described a simpler version of the MCFT, from which the shear design rules described in AASHTO LRFD [15] were derived and which also contains rules for shear design of circular members.Turmo et al. [16] proposed an analytical model for evaluating the contribution of the transverse reinforcement in concrete members of solid and hollow circular cross section, consisting in an accurate formula for evaluating the shear transferred by spiral reinforcement in solid members.The calculation of hollow core circular columns with both vertical and spiral reinforcement was also deduced.
Recently Merta [17] proposed a model based on the truss analogy by adding a concrete contribution term to the capacity of the shear reinforcement.An additional deviatoric shear resisting mechanism of hoops, present exclusively in members with curved transverse reinforcement, was identified and expressed analytically.It is explained by the fact that a curved reinforcing bar under tension induces compression in radial direction as well.The component of this compressive force in the direction of external shear could thus be considered as an additional shear enhancing mechanism of the hoops.Its magnitude is expressed through the friction force that is present between the concrete and steel after the section is cracked and the bond is partially destroyed.Under these assumptions, the concrete shear capacity was derived by a parametric study.

Experimental Approaches to Shear Capacity of Circular Members Modeling.
As for theoretical researches, also experimental studies on the shear capacity of circular members are quite limited.Assuming that the extrapolation of tests on rectangular beams may not be accurate for evaluating circular members shear strength, some researchers performed different test campaigns in such structural elements.Moreover, all the above-mentioned models are derived starting from a mechanical simplified model, whose parameters come from the extrapolation of experimental tests for evaluating circular member shear strength.In their paper, Clarke and Birjandi [18] reported also test results from different previous experimental studies, in particular: (i) Capon and Cossio [19] who tested 21 members, mainly of 250 mm in diameter, having only longitudinal steel, but four that had stirrups; eleven of the specimens are reported as failing in shear, two with stirrups and nine without; (ii) Khalifa and Collins [20] who tested five columns, all of 445 mm in diameter; these were subject to an axial load of 1000 kN, equivalent to a stress of 6.4 MPa, and all specimens failed in shear; (iii) Nagato [21] whose tests involved 16 columns, all with a diameter of 300 mm, with and without transverse reinforcement, and with and without axial stress; all specimens are reported as failing in shear.
Later on, Ghee et al. [8] tested 25 circular columns under axial load and cyclic lateral inelastic displacements, 24 of them failing in shear.The columns had a diameter of 400 mm.Variables in the test included axial load level, longitudinal reinforcement ratio, transverse reinforcement ratio, and aspect ratio.They addressed, for the first time, the efficiency of circular hoops and deduced the formula for calculating   presuming a 45-degree diagonal tension crack and vertical reinforcement.The work of Ghee et al. [8] was completed by Priestley et al. [9,22] which takes into account the potential development of steeper angles of cracking other than 45-degree one.Various experimental researches on circular RC members loaded monotonically in shear were aimed at verifying the possibility to extend to circular sections the use of design equations developed for rectangular sections [18][19][20].
Lee et al. [23] tested 4 full scale circular column specimens, which show shear related behavior under cyclic lateral load with constant axial force.The test variables are aspect ratio, transverse steel configuration (with or without crossties), and longitudinal steel ratio.The columns showed flexural failure or flexure-shear failure depending on the test variables.The test results were compared with shear strength equations adopted by the bridge seismic design specifications or guidelines and proposed by other researchers.The study showed that the accuracy of each method for shear strength prediction also depends on the aspect ratio and the relative steel amount of longitudinal reinforcement and transverse reinforcement.
More recently, Jensen and Hoang [24] presented the results of a test series on heavily shear reinforced circular concrete members.The specimens had shear reinforcement percentages up to more than three times the maximum percentage found in existing tests.The test results indicate that it is possible to obtain shear strengths which exceed the upper limit usually imposed on rectangular members.The test results are compared with a recently developed plasticity-based shear model for circular members.Satisfactory agreement was found.Comparisons were also made with calculations using the AASHTO LRFD [15] design code.It was found that the AASHTO LRFD design code gives reasonable results for members with small amounts of shear reinforcement while it underestimates the shear strength for heavily shear reinforced members.
This paper proposes an innovative approach, mainly based on experimental data, in order to develop an explicit compact expression to assess the shear strength of RC circular columns.Such approach is based on the evolutionary polynomial regression (EPR) [25], which is a hybrid data-driven technique that combines the effectiveness of evolutionary search with the advantage of classical numerical regression for developing data-driven models (i.e., input-output relationships) featured by explicit mathematical expression.The EPR technique has been successfully applied to modeling a wide range of complex engineering problems including constitutive modeling of soils; stability of slopes; settlement of foundations; liquefaction of soils due to earthquake and a number of other applications in civil engineering [26][27][28].
The EPR implements a multiobjective genetic algorithms (MOGA) optimization paradigm in order to search for model expressions as trade-off between fitness to data and complexity of mathematical formulations [29].The main advantage of the MOGA EPR approach is the possibility to select models from a set of optimal solutions according to the specific application and to the physical insight of the analyst.

Proposed Models for Shear Strength and Contemporary Code Provisions of Circular Reinforced Concrete Columns
Shear strength of reinforced concrete has received considerable attention in research during the past century.Several models for column shear strength have been proposed and used for design of new buildings and assessment of existing buildings.Shear reinforcements are used to ensure that the element fails in flexure; the addition shear reinforcement affects the mechanism by which shear is carried by a beam in several ways: shear reinforcement carries tensile actions across cracks, confines the compression-zone concrete increasing its shear capacity, encloses the flexural reinforcement, and can prevent dowel-splitting of the concrete.Moreover, for a given applied load, equilibrium of a cracked section with stirrup requires a shorter crack length but a larger crack width, that one without stirrup and the shape of the crack will also differ.Shear transfer mechanism in beam with stirrups has not been examined in so much detail as that in beams without stirrups.
Truss analogies are most commonly used in the design.The assumed internal equilibrium state comprises tensile shear reinforcement and inclined compressive struts of concrete.The original Mörsch truss analogy [30] uses a 45 ∘ strut angle and predicts failure when the shear reinforcements are at yield.The modified truss analogy [31] established an optimal lower bound for the shear capacity by varying the compressive strut angle to give reinforcement yield and web concrete failure simultaneously.
Even if shear transfer mechanisms are qualitatively well known, there is no agreement on the quantification of the shear strength of concrete members.As the value of   can be easily calculated with rational models, such as the truss analogy, research focuses on the elaboration of methods for an accurate evaluation of   .The concrete shear capacity is the capacity of the section without shear reinforcement: due to the complex stress redistribution after cracking, the shear transfer mechanisms of RC element have not been clearly understood yet.Its evaluation is very controversial and always relies, somehow, on empirical methods considering its dependence from various factors, such as the tensile longitudinal steel ratio, concrete grade, aggregate size, or shear span, among others.Present shear researches concern use of self-training algorithms to predict   or the development and application of equilibrium methods and equilibriumcompatibility methods, such as the Modified Compression Field Theory [13].The same interest has not been dedicated to improve the evaluation of the shear reinforcement contribution.
In this framework the shear strength of circular concrete columns has not received as much attention as rectangular ones and it is interesting to note that there are not guidelines for shear capacity in some international codes; likewise, many others prescribe that the shear capacity of a circular section equals the capacity of an equivalent rectangular section [32].In effect, this kind of approach has not been fully validated by test evidence and it obviously is questionable because in circular sections hoops contribute differently to the shear strength compared with rectangular stirrups.
In this section, a review of the shear strength provisions of various contemporary design code and proposed models is developed for circular columns.Most design codes include contributions from concrete   and transverse reinforcement   to analyze the shear strength of circular columns.The two components are then summed to estimate the total shear strength in One should point out that, generally, codified methods of design cannot be considered as predictive equations, since they are intended to provide only a conservative and safe lower bound to strength.Therefore, most shear design codes use empirical or semiempirical models for predicting shear strengths.The American code [15] and the New Zealand one [33] explicitly refer to circular concrete members, but with simplified methods.These recommend the calculation of the shear carried by the truss mechanism by representing the section as an equivalent rectangle having a width equal to the diameter  and an effective depth equal to 0.8.
Based on a large number of tests on circular cantilever columns under uniaxial load and multidirectional cyclic displacements, Ghee et al. [8] proposed a model for the shear strength of circular sections under cyclic load by adopting an additive approach based on the contribution from concrete and transverse reinforcement.The concrete contribution for low flexural ductility ( < 2) was defined as the "initial concrete strength" and calculated as a strength at maximum lateral load as follows: where  is the axial load,   is the cross-sectional area, and    is the concrete compressive strength.The effective area   is proposed as 0.8 times the cross-sectional area   , which approximately corresponds to the area of the confined concrete core.
For low aspect ratios (/ < 2), the shear capacity enhancement factor  was proposed as where  is the shear span and  is the section's diameter.The transverse reinforcement capacity was calculated from 45-degree truss mechanisms as follows: where  sh is the area of the shear reinforcement,  is the hoop's spacing,  yh is their yield strength, and   is the core diameter measured at the center line of the transverse hoop or spiral.It was assumed that all transverse reinforcements exposed by a presumed 45 ∘ diagonal crack are at yield.An integral averaging was imposed assuming that the spacing of the shear reinforcement is sufficiently small compared to the diameter.It was noted that for low   / ratios the equation can be up to 10% nonconservative.At higher flexural ductilities ( > 2), the degradation of the concrete produces a reduced concrete contribution and then a "final concrete strength" is proposed.One should point out that this formulation, such as other ones, implicitly assumes that all hoops exposed by a crack are at yield and also consider a "smeared" distribution of transverse reinforcement, whereas in reality transverse reinforcement consists of a discrete number of hoops that traverse a diagonal crack plane as shown in Figure 1.The later assumption is valid just if the pitch is sufficiently small compared to the diameter   .Clarke and Birjandi [18] proposed to use for circular sections the same shear design approach as given by the British codes of practice BS 5400 [34] for rectangular sections.
However, a modification was suggested; that is, that the section's effective depth should be considered as the distance from the extreme compression fiber to the centroid of the tension reinforcement.
Following this, the effective shear area represents the area corresponding to the effective depth.
The design approach for rectangular section in British codes of practice [34] consists in adding together the concrete and transverse reinforcement contributions to the shear capacity.The concrete term is defined as where  sl is the area of longitudinal steel,   is the section's width,  is the effective depth, and  cu is the concrete cube strength.For loads applied at a distance  V closer than 2 from the support, the concrete shear capacity is increased by The shear force carried by transverse reinforcement is calculated as where  sh is the cross-sectional area of the link's both legs at the section's neutral axis.In a member with axial compression load , the shear capacity should be multiplied by (1 + 0.05/  ).Dancygier [35] modified Ghee's model by adopting a discrete computation of tension-force components developed by reinforcing hoops in the direction of the shear force, given by the following expression: where  is the index of the hoop that is crossed by the inclined crack,  is the angle between its force ( sh  yh ) and the direction of the shear force, given by sin and  is the number of the hoops that are crossed by the crack: where  is the pitch,  is crack's inclination angle, and  0 is the distance from the end of the crack to the first hoop that is crossed by it.
Kowalsky and Priestley [11] discussed Ghee's model considering that in the derivation of the truss mechanism equation it was assumed that a diagonal crack is able to mobilize transverse reinforcement along a crack length extending the full width of the confined core of the concrete.In the compression zone of the column, however, all cracks are by definition closed.Therefore, if the crack is closed, shear cannot be transferred across it by tension strain in the transverse reinforcement.Therefore, it is evident that a reduced column width − −  − cov − (Figure 1) is appropriate for calculating the number of hoops mobilized by the cracks between the compression struts and then a revised truss component equation should be adopted.The authors obtained an approximated equation given by where   is the area of the hoop's one leg,   is the hoop's yield strength,  is the compression zone's depth, and cov is the concrete cover.
With reference to the concrete mechanism, the following equation was discussed in the cited work: where the parameter  accounts for the column aspect ratio, the factor  is a modified factor that accounts for the longitudinal steel ratio, and  represents the reduction in strength of the concrete shear resisting mechanism with increasing ductility.Also Kim and Mander [36] reported a deficiency of the model originally proposed by Ghee et al. [8] and verified that the applied integral averaging by the reinforcement term restricts its use only to members with diameter at least four times the spacing of the shear reinforcement.For all other ratios, the formula could be even more than 50% nonconservative.For this reason, they proposed a corrective coefficient to convert the amount of shear implied by a continuous truss model to the discrete model.
The current ACI 318-M-08 [37] considers a portion of the design shear force to be carried by the concrete shear resisting mechanism,   , with the remainder carried by truss mechanism,   , involving transverse reinforcement.The code proposes the following equation to calculate   for members subjected to combined shear, moment, and axial compression: The transverse reinforcement contribution is also calculated as where the meaning of various terms has been explained before.
Merta [17] developed a simple shear capacity model of RC circular sections.It was based on the truss analogy by adding an empirical concrete contribution term to the capacity of the shear reinforcement, obtaining finally the following relation: where   represents the number of hoops across the crack and is evaluated considering only the hoops outside the compression zone: while  is assumed to be equal to 0.53 and   is given by being  ≅ 0.3.The influence of the main variables on the shear capacity, such as the longitudinal reinforcement ratio   , the axial load level /  , and the shear span-to-depth ratio /, has been determined empirically, based on a total of 44 pieces of data of circular cross-section specimens without shear reinforcement under monotonic load.The term 0.7  represents the section's effective shear area, while  = 1 if / > 2.5 and 1.25 if / ≤ 2.5.

A Data-Driven Approach for Estimation of Shear Strength of Reinforced Concrete Circular Columns: The Evolutionary Polynomial Regression Paradigm
Numerical Regression is the most powerful and commonly applied form of regression that provides a solution to the problem of finding the best model to fit the observed data (e.g., fitting a line through a set of points).However, the functional form (linear, exponential, logarithmic, etc.) has to be selected before fitting commences.On the other hand, genetic pogramming uses simple but very powerful artificial intelligence tactics for computer learning inspired by natural evolution to find the appropriate mathematical model to fit a set of points.The computer produces and evolves a whole population of functional expressions based on how closely each of them fits the data.The evolutionary polynomial regression (EPR) is a recently developed hybrid regression method by Giustolisi and Savic [25] that integrates the best features of numerical regression with genetic programming [39,40].An example of the general model structures that EPR can manage is reported in ES(,1) ⋅ . . .⋅ (X  )

ES(𝑗,𝑘)
⋅  ((X 1 ) ES(,+1) ⋅ . . .⋅ (X  ) ES(,2) ) , where  is the number of additive terms,   are numerical parameters to be estimated, X  are candidate explanatory variables, ES(, ) (with  = 1, . . ., 2) is the exponent of the th input within the th term in (18), and  is a user-selected function among a set of possible alternatives (including no function selection).The exponents ES(, ) are selected from a user-defined set of candidate values.Note that it is recommended to include 0 among such exponents in order to allow EPR to exclude some of the candidate explanatory variables in X  from those present in the returned models, thus selecting the most influencing input variables on the target prediction.
Starting from a general structure, such as that in ( 18), the final model expression is obtained through a population based strategy that mimics the evolution of the individuals in nature; it is a combination of a genetic algorithm (GA) [41] paradigm for finding the best exponents to be assigned to each candidate input and a least squares (LS) regression [42], for the identification of the constant values   [25].
Recently, the EPR has been upgraded to MOGA EPR through the implementation of a search pattern of the best models based on a multiobjective genetic algorithm (MOGA) [43] optimization scheme, being able to find several model expressions maximizing fitness to data and parsimony of mathematical formulations simultaneously [29].MOGA EPR explores the space of -terms formulas using at least three objectives among the following conflicting ones: (i) maximization of model accuracy, (ii) minimization of the number of model coefficients, and (iii) minimization of the number of actually used model inputs (i.e., whose exponent is not null in the resulting model structure).Based on the Pareto dominance criterion [44], MOGA EPR finally obtains a set of optimal solutions (i.e., the Pareto front) which can be considered as trade-offs between structural complexity and accuracy.MOGA EPR displays the returned formulas, within the Pareto front, as ranked according to their structural complexity; this way, the comparisons among formulas according to different criteria (e.g., selected inputs, structure, etc.) are facilitated, as well as the selection of models suited for different modeling purposes [29].
In particular, this possibility has been exploited in this paper in order to find the best model for the estimation of shear strength of reinforced concrete circular columns.Thus, the first step was to define the MOGA EPR search options (candidate model attributes, candidate exponents for attributes, maximum number of parameters, objective functions, etc.); then, the set of optimal models, with respect to the objective functions, produced by MOGA EPR has been analyzed by considering (i) the model structure with respect to physical insight of the phenomenon at stake; (ii) the similarities and/or recurrent groups of variables within the mathematical structures of the Pareto models; (iii) the generalization performance of models in terms of statistical indicators and mathematical parsimony.Moreover, the returned models have been also evaluated looking at their consistency with existing codes/theories about the physical phenomenon, as described above.All these issues can lead to a quite robust model selection.Figure 2 is illustrative of the decision support framework allowed by MOGA EPR, as here implemented.
The present analyses have been performed using the EPR MOGA-XL v.1 which is a Microsoft Office Excel add-in function freely available at www.hydroinformatics.it, which employs a MOGA optimizer named optimized multiobjective genetic algorithm (OPTIMOGA) [45].

Results and Comparison with Literature and Technical Code Models
The proposed model has been verified by comparing it to existing ones and to experimental tests.The experimental data herein considered have been obtained from Merta [17], Capon and Cossio [19], Khalifa and Collins [20], Clarke and Birjandi [18], and Kim [38] and include the results of several tests performed by different authors.The test specimens refer to circular columns characterized by six different diameters and having both shear and longitudinal reinforcements.
Longitudinal reinforcement consists mainly of high yield bars placed evenly at eight locations round the pitch circle.Shear reinforcement is generally in the form of 6 mm or 8 mm mild steel links.The compressive strength of concrete ranges from 13.2 MPa to 49.3 MPa, the stirrup percentage ranges from 0.1% to 0.45%, the yielding stress of shear reinforcement ranges from 250 MPa to 1728 MPa, and the longitudinal steel percentage ranges from 2.2% to 5.6%.A total of 61 circular columns have been selected with different geometric and mechanical parameters such as dimension of the section diameter (), area of concrete (  ), total cross-section of links at the neutral axis of the section ( sh ), spacing of links along the member (), concrete compressive strength (  ), yield strength of closed links ( yh ), total area of longitudinal reinforcement ( sl ), yield strength of longitudinal reinforcement ( yl ), concrete cover (cov), effective depth (), and axial compression load ().The range of datasets is listed in Table 1.
The main objective of this study is to investigate the application of EPR MOGA technique to get possible new formulations of shear strength of circular RC columns able to reproduce the experimental data better than the existing code expressions.In fact the correlation coefficient  2 between experimental data and previsions of ACI 318 M-08 [37] and BS 8100 [33] are 90.15% and 76.92%, respectively; thus, there is some room for possible improvements.
Since most of experimental data refer to pure-shear tests, the dependence of shear strength from the axial force  is not investigated.It can be accounted for by adding a term representative of the axial force contribution, as the one proposed in Merta or in technical code models.

EPR MOGA Search Settings. Bearing in mind what
above mentioned about MOGA EPR search paradigm, some compact expressions will be sought for, in terms of both number of parameters and explanatory variables involved.This, in turn, would facilitate their physical interpretation and technical soundness for practical purposes.The base  model structure reported in ( 18) is adopted with no function  selected; accordingly, each additive monomial term is assumed to be a combination of the input variables raised to relevant exponents.The candidate explanatory variables (i.e., model inputs, see Table 1) are , ,   ,  sh , ,  sl ,   , and  yh ; they have been selected by both looking at the existing literature and performing some preliminary sensitivity analyses.The candidate exponents are [−1, −0.5, −2/3, 0, 0.5, 2/3, 1] in order to represent both linear and nonlinear relationships as well as direct/inverse dependence between candidate input(s) and output (i.e., model target).The maximum number of additive term in the final expressions is assumed to be  = 3, so that could be easily compared with the existing formulations.
Although the EPR MOGA search for models is based on the coefficient of determination (COD) as defined in [25], model performances are reported in terms of  2 in order to allow a straight comparison with code provisions [33,37] and alternative formulations [17,18].
Proposed models are compared to ACI 318-M-08 [37], BS 8100 [33] predictions, and the two existing models of Clarke and Birjandi [18] and Merta [17]; the corresponding correlation coefficients  2 with respect to the experimental data herein considered are, respectively, equal to 90.15%, 76.92% 96.23%, and 95.29%.Figure 3 depicts the shear strength as predicted by the EPR models ( 19)-( 23) versus the experimental results; similarly, Figure 4 shows the same graphs with reference to the code formulations and the literature models mentioned above.As expected from the values of the coefficients  2 , expression (20) provides the smallest scatter in the plot of the experimental/theoretical ultimate shear strength of specimens; however, the complexity of expression (20) discourages its use for practical purposes.Formulas (19), (22), and (23) provide predictions similar to the Clarke and Birjandi ones and are anyway slightly better than the Merta ones.
Formula (23) predictions are more conservative than Clarke and Birjandi ones, with a similar or often closer correlation with experimental data.It is also evident that the ACI 318-M-08 and BS 8100 predictions, which neglect the influence of the longitudinal reinforcement contribution, are highly conservative.In fact, BS 8100 predictions in Figure 4 are calculated by adopting cot  = 2.5; otherwise, by putting cot  = 1 more conservative values of shear strength are obtained.
In expressions (19) and (20), the physical consistency of the problem is not easily recognizable, while on the contrary models ( 21)-( 23) are similar to the building code formulations and thus preserve the physical insight into the problem.In particular, (21) differs from the ACI 318-M-08 [37] prediction only for the multiplying factor 0.25274 of the concrete contribution which is higher than the code one, allowing a better agreement with the experimental values of shear capacity.
Equation (21) implies that the shear capacity of the member is the sum of the shear carried by concrete and transverse reinforcement; nevertheless, according to this approach, which is consistent with code formulation, a maximum correlation coefficient equal to 91.53% can be obtained.On the contrary, in the other proposed models also the contribution of longitudinal reinforcement is accounted for, allowing gaining higher values of the correlation coefficient.
In particular in ( 22) and ( 23) the concrete and transverse reinforcement contributions preserve the same shapes as ACI 318-M-08, while the longitudinal reinforcement contribution is function of the concrete tensile strength, which in fact is proportional to √  according to ACI 318-M-08 [37] or Vecchio and Collins [13] and to  2/3  according to Italian code [46].Moreover, the longitudinal reinforcement contribution in ( 23) is similar to the Italian code expression valid in the case of rectangular sections in the absence of stirrups [46].These results suggest that the third contribution to the shear capacity of circular sections is due to the dowel effect; longitudinal reinforcement affects the strain state of the element and the crack widths, influencing the concrete contribution by aggregate interlock to the shear strength.
One of the most significant advantages of EPR is that this technique is capable of representing the relationship between the parameters in the form of a polynomial expression, so parametric studies can be conducted to investigate the contribution of different parameters involved and to assess whether the developed model is able to capture the physical relationships between various parameters of the system.In the context of the present study, this also gives an insight to the user about the degree of influence that different parameters can have on the shear of circular sections.To achieve this, a parametric study is carried out where all parameters are set to their mean values except one that is changed from its minimum to maximum values in the testing set.The results of the parametric study for the shear strength model of circular elements are shown in Figure 5 with reference to code provisions (ACI 318-M-08, BS 8100) and ( 20), (22).By adopting formula (23), graphs quite similar to formula (22) ones are obtained, which are omitted here for sake of brevity.
The results show that the shear strength increases as the concrete compressive strength, the amount of transverse reinforcement, the amount of longitudinal reinforcement, and the section diameter increase.In particular, the section diameter has a significant influence on the value of the shear strength.This is consistent with the expected shear behavior of circular sections.It emerges that the EPR models, directly developed using experimental data, are able to capture the various aspects of shear behavior of circular section correctly and with higher accuracy than exiting technical codes.Moreover, differently from ACI 318-M-08 and BS 8100 predictions, the proposed EPR models can also capture the influence of longitudinal reinforcement on the shear strength.In fact, according to the theoretical models, shear strength should be proportional to the longitudinal V Merta (kN) V ACI 318 (kN) V BS 8100 (kN) reinforcement ratio; this is due to the circumstance that as the longitudinal reinforcement ratio decreases, flexural stresses and strains in concrete increase, thus leading to larger crack widths; in addition, dowel action is weakened due to the low longitudinal reinforcement ratio.So the proposed models are coherent with this mechanical behavior.The above considerations clearly indicate that the proposed models (see (22) and ( 23)) improve the prediction of the shear capacity preserving the physical insight into the problem; moreover, they represent sufficiently conservative design equations and their simplicity makes them suitable for incorporation in technical codes.

Conclusions
In the present paper, new formulas for calculating the shear strength of RC circular columns have been proposed.Differently from standard approaches, they are obtained by EPR.In the EPR approach, no preprocessing of the data is required and there is no need for normalization or scaling of the data.The efficiency of such approach has been tested by using experimental data of 61 circular RC elements reported in technical literature.The considered inputs are dimension of the section diameter, area of concrete, total cross-section of links at the neutral axis of the section, spacing of links along the member, concrete compressive strength, yield strength of closed links, total area of longitudinal reinforcement, yield strength of longitudinal reinforcement, and effective depth.Models obtained from EPR MOGA runs show different number of terms and explanatory variables as well as different accuracy.Among all the proposed models, the selected ones provide good predictions, particularly with respect to the code ones.The increased accuracy descends from including the contribution of the dowel effect exercised by the longitudinal reinforcement in the expression of the shear strength.

Figure 1 :
Figure 1: Geometrical and mechanical characteristics affecting the shear strength of a circular column.

1 Figure 2 :
Figure 2: Decision support framework for model selection based on EPR MOGA paradigm.

Figure 4 :
Figure 4: Shear strength of circular section members with shear reinforcement: code and literature models.

Table 1 :
Details and test results of specimens with circular cross-section.
4.2.Selection of EPR MOGA Models.Two EPR MOGA runs were finally selected, characterized by different sets of exponents and by the same set of input variables (, ,   ,  sh , ,  sl ,   ,  yh ).Models obtained from EPR MOGA runs show different number of terms and explanatory variables as well as different accuracy.The model selection procedure takes into account some key aspects: agreement with experimental data, model parsimony and consistency with both physical insight about the shear behavior of circular RC columns, and previous building code formulations.From among all models obtained, the following formulas have been selected (the unit for  Rd is [N] and all the mechanical and geometrical quantities are expressed in[N]and [mm] resp.): Rd = 4.1525 ℎ  yh + 0.13552  √  √  ,  Rd = + 7.121√   yh + 0.043393 ℎ   √ yh   + 0.013495    √ , ( 2 = 97.84%),