ANN-Based Shear Capacity of Steel Fiber-Reinforced Concrete Beams without Stirrups

Abstract

Comparing experimental results of the shear capacity of steel fiber-reinforced concrete (SFRC) beams without stirrups to the capacity predicted using current design equations and other available formulations shows that predicting the shear capacity of SFRC beams without mild steel shear reinforcement is still difficult. The reason for this difficulty is the complex mechanics of the problem, where the steel fibers affect the different shear-carrying mechanisms. Since this problem is still not fully understood, we propose the use of artificial intelligence (AI) to derive an expression based on the available experimental data. We used a database of 430 datapoints obtained from the literature. The outcome is an artificial neural network-based expression to predict the shear capacity of SFRC beams without shear reinforcement. For this purpose, many thousands of artificial neural network (ANN) models were generated, based on 475 distinct combinations of 15 typical ANN features. The proposed “optimal” model results in maximum and mean relative errors of 0.0% for the 430 datapoints. The proposed model results in a better prediction (mean V_test/V_ANN = 1.00 with a coefficient of variation 1 × 10⁻¹⁵) as compared to the existing code expressions and other available empirical expressions, with the model by Kwak et al. giving a mean value of V_test/V_pred = 1.01 and a coefficient of variation of 27%. Until mechanics-based models are available for predicting the shear capacity of SFRC beams without shear reinforcement, the proposed model thus offers an attractive solution for estimating the shear capacity of SFRC beams without shear reinforcement. With this approach, designers who may be reluctant to use SFRC because of the large uncertainties and poor predictions of experiments, may feel more confident using the material for structural design.

1. Introduction

Because concrete is strong in compression but weak in tension, adding fibers to the material can be a solution for the limited strength in tension. In structural applications, fiber-reinforced concrete is combined with regular reinforcement steel. The type of fibers that are used, are most often steel fibers. These fibers help to distribute cracks and keep the crack widths small [1].

One failure mode where crack shape and width is essential, is shear failure. When steel fibers are included in the concrete mix, and the reinforced concrete element built with this concrete mix is tested in shear, then the addition of the steel fibers influences all mechanisms that contribute to the shear-carrying capacity of the member [2]. Since the mechanics of the problem are still not fully understood, it is interesting to study the behavior of steel fiber-reinforced concrete (SFRC) with longitudinal steel reinforcement and without shear reinforcement. As such, we can study the contribution of steel fibers to the shear capacity of structural concrete, without the influence of the stirrups as shear reinforcement. Once this problem is understood, we can find an optimal combination of steel fiber reinforcement and regular stirrups to act as shear reinforcement. Such a combination is particularly interesting in joints where rebar congestion can make concreting difficult [3], and for bridge girders, where the addition of steel fibers can lead to a more durable structure as cracks will be more distributed and crack widths will be smaller.

Even though the mechanics of the shear resistance of SFRC is not fully understood, several design expressions are available in the literature and current codes. These expressions are mostly semi-empirical expressions, with the exception of extensions of the Modified Compression Field Theory [4,5,6,7,8,9,10,11,12], the Dual Potential Capacity Model [13,14], and plasticity-based models [15,16,17,18,19].

Table 1. Shear prediction equations from literature and available codes, adapted from [20].

Authors	Ref	Expression
Sarveghadi et al.	[22]	$V_u=\left[\rho +\frac{\rho }{v_b}+\frac{1}{\frac{a}{d}}\left(\frac{\rho f_t\prime \left(\rho +2\right)\left(f_t\prime \frac{a}{d}-\frac{3}{v_b}\right)}{\frac{a}{d}}+f_t\prime \right)+v_b\right]b_{w}d$	(2)
		$f_t\prime =0.79\sqrt{f_c\prime }$	(3)
		$v_b=0.41\tau F$ with τ = 4.15 MPa	(4)
Kwak et al.	[23]	$V_u=\left[3.7e{f}_{spfc}^{2/3}{\left(\rho \frac{d}{a}\right)}^{1/3}+0.8v_b\right]b_{w}d$	(5)
		$f_{spfc}=\frac{f_{cuf}}{\left(20-\sqrt{F}\right)}+0.7+1.0\sqrt{F}$ in MPa	(6)
		$e=\{\begin{array}{l}1\mathrm{for}\frac{a}{d}>3.4\\ 3.4\frac{d}{a}\mathrm{for}\frac{a}{d}\le 3.4\end{array}$	(7)
Greenough and Nehdi	[24]	$V_u=\left[0.35\left(1+\sqrt{\frac{400}{d}}\right){\left(f_c\prime \right)}^{0.18}{\left(\left(1+F\right)\rho \times 100\times \frac{d}{a}\right)}^{0.4}+0.9{\eta }_o\tau F\right]b_{w}d$	(8)
Kuntia et al.	[25]	$V_u=\left[\left(0.167+0.25F\right)\sqrt{f_c\prime }\right]b_{w}d$	(9)
Sharma	[26]	$V_u=\left(\frac{2}{3}\times 0.8\sqrt{f_c\prime }{\left(\frac{d}{a}\right)}^{0.25}\right)b_{w}d$	(10)
Mansur et al.	[27]	$V_u=V_c+{\sigma }_{tu}{b}_{w}d$	(11)
		$V_c=\left(0.16\sqrt{f_c\prime }+17.2\frac{\rho Vd}{M}\right)b_{w}d\le 0.29\sqrt{f_c\prime }{b}_{w}d$	(12)
		${\sigma }_{tu}=3.2{\eta }_o{\eta }_{l}F\tau \mathrm{with}\tau =2.58\mathrm{MPa}$	(13)
		${\eta }_l=1-\frac{\tanh \left(\beta \frac{l_f}{2}\right)}{\beta \frac{l_f}{2}}$	(14)
		$\beta =\sqrt{\frac{2\pi G_m}{E_f{A}_f\ln \left(\frac{S}{r_f}\right)}}$	(15)
		$S=25\sqrt{\frac{d_f}{V_f{l}_f}}$	(16)
Ashour et al.	[28]	$V_u=\left[\left(0.7\sqrt{f_c\prime }+7F\right)\frac{d}{a}+17.2\rho \frac{d}{a}\right]b_{w}d$	(17)
		$V_u=\left[\left(2.11\sqrt[3]{f_c\prime }+7F\right){\left(\rho \frac{d}{a}\right)}^{0.333}\right]b_{w}d\mathrm{for}\frac{a}{d}\ge 2.5$	(18)
		$V_u=\left[\left(\left(2.11\sqrt[3]{f_c\prime }+7F\right){\left(\rho \frac{d}{a}\right)}^{0.333}\right)\frac{2.5}{\frac{a}{d}}+v_b\left(2.5-\frac{a}{d}\right)\right]b_{w}d\mathrm{for}\frac{a}{d}<2.5$	(19)
Arslan et al.	[29]	$V_u=\left[\left(0.2{\left(f_c\prime \right)}^{2/3}\frac{c}{d}+\sqrt{\rho \left(1+4F\right)f_c\prime }\right)\sqrt[3]{\frac{3}{\frac{a}{d}}}\right]b_{w}d$	(20)
		${\left(\frac{c}{d}\right)}^2+\left(\frac{600\rho }{f_c\prime }\right)\left(\frac{c}{d}\right)-\frac{600\rho }{f_c\prime }=0$	(21)
Imam et al.	[30]	$V_u=\left[0.6\psi \sqrt[3]{\omega }\left({\left(f_c\prime \right)}^{0.44}+275\sqrt{\frac{\omega }{{\left(\frac{a}{d}\right)}^5}}\right)\right]b_{w}d$	(22)
		$\psi =\frac{1+\sqrt{\frac{5.08}{d_a}}}{\sqrt{1+\frac{d}{25d_a}}}$	(23)
		$\omega =\rho \left(1+4F\right)$	(24)
Yakoub	[31]	$V_u=\left[0.83\xi \sqrt[3]{\rho }\left(\sqrt{f_c\prime }+249.28\sqrt{\frac{\rho }{{\left(\frac{a}{d}\right)}^5}}+0.405\frac{l_f}{d_f}{V}_f{R}_g\frac{d}{a}\sqrt{f_c\prime }\right)\right]b_{w}d\mathrm{for}\frac{a}{d}\le 2.5$	(25)
		$V_u=\left[0.83\xi \sqrt[3]{\rho }\left(\sqrt{f_c\prime }+249.28\sqrt{\frac{\rho }{{\left(\frac{a}{d}\right)}^5}}+0.162\frac{l_f}{d_f}{V}_f{R}_g\sqrt{f_c\prime }\right)\right]b_{w}d\mathrm{for}\frac{a}{d}\ge 2.5$	(26)
		$\xi =\frac{1}{\sqrt{1+\frac{d}{25d_a}}}$	(27)
		$V_u=2.5\left(\frac{0.40}{1+1500{\epsilon }_x}\times \frac{1300}{1000+s_{xe}}\right)\sqrt{f_c\prime }\left(1+0.7\frac{l_f}{d_f}{V}_f{R}_g\right)\frac{d}{a}{b}_w{d}_v\mathrm{for}\frac{a}{d}\le 2.5$	(28)
		$V_u=\left(\frac{0.40}{1+1500{\epsilon }_x}\times \frac{1300}{1000+s_{xe}}\right)\sqrt{f_c\prime }\left(1+0.7\frac{l_f}{d_f}{V}_f{R}_g\right)b_w{d}_v\mathrm{for}\frac{a}{d}\ge 2.5$	(29)
		$d_v=\max \left(0.9d,0.72h\right)$	(30)
		${\epsilon }_x=\frac{\frac{M}{d_v}+V}{2E_s{A}_s}$	(31)
		$s_{xe}=\frac{35s_x}{16+d_a}\ge 0.85s_x\mathrm{and}{s}_x\approx d_v$	(32)
Association Française de Génie Civil	[32]	$V_{Rd}=V_{Rd,c}+V_{Rd,f}$	(33)
		$V_{Rd,c}=\frac{0.21}{{\gamma }_{cf}{\gamma }_E}{f}_{ck}^{1/2}{b}_{w}d$ with γcfγE = 1.5	(34)
		$V_{Rd,f}=\frac{A_{vf}{\sigma }_{Rd,f}}{\tan \theta }$ with θ ≥ 30o	(35)
		${\sigma }_{Rd,f}=\{\begin{array}{l}\frac{1}{K{\gamma }_{cf}}\frac{1}{w_{\lim }}{\displaystyle {\int }_0^{w_{\lim }}{\sigma }_f\left(w\right)dw}\mathrm{for}\mathrm{low}\mathrm{strain}\mathrm{hardening}\mathrm{or}\mathrm{softening}\\ \frac{1}{K{\gamma }_{cf}}\frac{1}{{\epsilon }_{\lim }-{\epsilon }_{el}}{\displaystyle {\int }_{{\epsilon }_{el}}^{{\epsilon }_{\lim }}{\sigma }_f\left(\epsilon \right)d\epsilon }\mathrm{for}\mathrm{high}\mathrm{strain}\mathrm{hardening}\end{array}$ with K = 1.25 or based on tension tests on the SFRC mix	(36)
		$w_{\lim }=\max \left(w_u,w_{\max }\right)$	(37)
		${\epsilon }_{\lim }=\max \left({\epsilon }_u,{\epsilon }_{\max }\right)$	(38)
		$A_{vf}=b_{w}z$	(39)
DAfStB	[33]	$V_{Rd,c}^f=V_{Rd,c}+V_{Rd,cf}$	(40)
		$V_{Rd,c}=\frac{C_{Rd,c}}{{\gamma }_c}k{\left(100\rho f_{ck}\right)}^{1/3}{b}_{w}d>V_{Rd,c,\min }$ with CRd,c = 0.15 and γc = 1.5, ρ ≤ 2%	(41)
		$V_{Rd,cf}=\frac{{\alpha }_c^f{f}_{ctR,u}^f{b}_{w}h}{{\gamma }_{ct}^f}$ with ${\gamma }_{ct}^f=1.25$ and ${\alpha }_c^f=0.85$	(42)
		$f_{ctR,u}^f=k_F^f{k}_G^{f}0.37f_{cfIk,L2}^f$ with $k_F^f=0.5$	(43)
		$k_G^f=1.0+0.5A_{ct}^f\le 1.7$	(44)
		$A_{ct}^f=b_w\times \min \left(d,1.5m\right)$	(45)
		$k=1+\sqrt{\frac{200mm}{d}}$	(46)
RILEM	[34]	$V_{Rd}=V_{cd}+V_{fd}$	(47)
		$V_{cd}=0.12k{\left(100\rho f_{ck}\right)}^{\frac{1}{3}}{b}_{w}d$ with ρ ≤ 2%	(48)
		$V_{fd}=0.7k_{f}k{\tau }_{fd}{b}_{w}d$	(49)
		$k_f=1+n\left(\frac{h_f}{b_w}\right)\left(\frac{h_f}{d}\right)\le 1.5$	(50)
		$n=\frac{b_f-b_w}{h_f}\le 3\mathrm{and}n\le \frac{3b_w}{h_f}$	(51)
		${\tau }_{fd}=0.12f_{Rk,4}$	(52)
fib	[35]	$V_{Rd}=V_{Rd,f}=\frac{C_{Rd,c}}{{\gamma }_c}k{\left(100\rho \left(1+7.5\frac{f_{Ftuk}}{f_{ctk}}\right)f_{ck}\right)}^{1/3}{b}_{w}d$ with CRd,c = 0.18, γc = 1.5, and ρ ≤ 2%	(53)
		$f_{ctk}=\{\begin{array}{l}0.3{\left(f_{ck}\right)}^{2/3}\mathrm{for}\mathrm{concrete}\mathrm{grades}\le \mathrm{C}50\\ 2.12\ln \left(1+0.1\left(f_{ck}+8MPa\right)\right)\mathrm{for}\mathrm{concrete}\mathrm{grades}>\mathrm{C}50\end{array}$	(54)
CNR-DT	[36]	$V_{Rd}=V_{Rd,f}\ge V_{min}$ with VRd,f from Equation (53)	(55)
		$V_{min}=0.035k^{3/2}{f}_{ck}^{1/2}{b}_{w}d$	(56)

gives an overview of the expressions for determining the shear capacity of SFRC beams without shear reinforcement from the literature that were considered in this study for comparison [20]. In these expressions, the properties of the fibers are often described by the fiber factor F [19]:

$$F=V_f\frac{l_f}{d_f}{\rho }_f$$

(1)

with V_f the fiber volume fraction, l_f the length of the fiber, d_f the diameter of the fiber, and ρ_f the fiber bond factor, which depends on the fiber type. An overview of the notations used in Table 1 is given in the List of Notations at the end of this article.

For small values of a/d (typically a/d < 2.5) the shear behavior is different than for larger shear span to depth ratios. The reason for this difference in behavior is that for short shear spans a compressive strut can develop between the point of application of the load and the support [21]. This additional load-carrying mechanism enhances the shear capacity. In Table 1, this behavior is reflected by the inclusion of a factor, such as e in Equation (7) to enhance the shear capacity for short shear spans. Not all expressions from the literature include this effect. Where this effect is not included, the predicted shear capacities tend to be conservative for short shear spans.

Comparing experimental results on the shear capacity of SFRC beams without stirrups to the capacity predicted using the expressions from Table 1 shows that predicting the shear capacity of SFRC beams without mild steel shear reinforcement is still difficult [20]. The reason for this difficulty is the complex mechanics of the problem, since, as previously mentioned, the steel fibers affect each of the shear-carrying mechanisms [2].

Since the problem of the shear capacity of SFRC elements without shear reinforcement is still not fully understood, we propose the use of machine learning, and in particular the use of artificial neural networks (ANNs or neural nets) to derive an expression based on the available experimental data. Machine learning is a subdiscipline of artificial intelligence (AI) in which computers learn how to carry out tasks based on examples of how they should be done [37]. When we have a lot of data but no theoretical foundations to explain these, machine learning can be a suitable tool. An illustration of the several possible scenarios is presented in [38], see Figure 1, where the shadowed areas represent regions where any of the contiguous tools might be used. The ANN is the oldest [39] and most powerful [40] machine learning technique. Generally speaking, an ANN is an analytical model for a singular task, functioning similarly to the human brain by using neurons. ANN is more powerful than traditional approaches (e.g., multi-variate nonlinear regression) and does not require prior knowledge of the shape of the function that will be approximated [41].

Several the expressions in Table 1 have been developed with methods of AI. The equation by Sarveghadi et al. [22] is derived from a more general expression developed with the use of artificial neural networks. The expression by Greenough and Nehdi [24] results from genetic programming. Other studies that used AI to evaluate the shear capacity of SFRC beams are the study by Hossain et al. [42] based on 173 experiments, resulting in an ANN model with 5 hidden neurons. The mean square error of the comparison between the experiments and the proposed model was 3.0665 and the root mean square error 1.7512. The study by Kara [43] used gene expression programming. In total, 101 experiments from the literature were used. Comparison between the proposed model and the experiments gave an absolute average error of 11.39% and a coefficient of variation of 15.42%. Other recent research works that have considered databases of experiments have used smaller databases than this study: 122 experiments [11] or 171 experiments (of which 93 with Fiber Reinforced Concrete) [44].

The goal of this study is to derive an ANN-based model to predict the shear capacity of SFRC elements without stirrups within the ranges of the available experimental data. While similar studies have been carried out in the past, as discussed in the previous paragraph, our study is an improvement of the state-of-the-art for the following reasons: (1) the dataset used for this study contains 430 experiments, which is significantly larger than the datasets used in the previous studies; (2) a large number of ANN features were varied, resulting in 475 combinations of features, with which the model that best predicts the experimental results could be selected; and (3) the error of the comparison between the proposed model and the experimental results is 0%, showing that the model proposed with this study is a significant improvement as compared to the previously cited models.

2. Materials and Methods

2.1. Data Gathering

We used a database of test results reported in the literature as input for the model. The database is earlier reported in [20]. To create unique datapoints, the outcome of repeat tests is averaged. Therefore, the original dataset of 488 test results is reduced to 430 unique datapoints considered for this study. As a result, some of the inherent scatter observed in repeat tests is removed from the database. The reader should keep this preprocessing in mind when evaluating the performance of the proposed model. The experimental results are taken from the literature [3,17,23,24,26,27,28,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102].

Table 2. Overview of input and output variables considered in the dataset, including ranges of values.

Input Parameters			Input Number	Min.	Max.
Geometry	b (mm)	width	1	50	610
	d (mm)	effective depth	2	85.3	1118
	av/d (-)	clear shear span to depth ratio	3	0.2	6.0
Properties of reinforcement	ρ (-)	reinforcement ratio	4	0.004	0.057
	fy (MPa)	yield strength of steel	5	257.9	900
Concrete properties	da (mm)	maximum aggregate size	6	0.4	22
	fc,cyl (MPa)	average concrete compressive strength	7	9.8	215
Fiber properties	F (-)	fiber factor	8	0.1	2.9
	ftenf (MPa)	tensile strength of fiber	9	260	4913
Output	Vutot (kN)	sectional shear capacity	1	12.9	1480.9

gives an overview of the input and output values considered in this study. The geometry is described based on the effective depth d, the width b and the clear shear span to effective depth ratio a_v/d. Figure 2 shows the geometry of a typical beam specimen. The properties of the steel reinforcement are determined based on the reinforcement ratio ρ, which is defined as follows:

$$\rho =\frac{A_S}{bd}$$

(57)

with A_s the area of longitudinal reinforcement. The other parameter determining the reinforcement is the yield strength of the steel f_y. The properties of the concrete mix are described by the maximum aggregate size d_a and the measured average concrete cylinder compressive strength f_c. The properties of the fibers are given through the fiber factor F according to Equation (1) as well as the tensile strength of the fibers f_tenf. The output is the sectional shear capacity V_utot, determined as shown in Figrue 2b. V_utot includes the effect of the self-weight of the beam. In total, nine input variables and one output variable are selected.

Table 2 also gives an overview of the ranges of the input variables in the dataset used for the development of the model. The currently available test results come from relatively small specimens, as can be seen in Table 2. With a maximum effective depth of 1118 mm, the available experimental data may not sufficiently address the size effect in shear for SFRC [103,104,105,106]. The data moreover show that most experimental data come from small-scale specimens. A wide range values for a_v/d is covered, so that the model resulting from the database can be used for deep members as well as for slender members.

The reinforcement ratios used in the experiments cover a wide range. The majority of SFRC shear experiments are carried out on heavily reinforced beams, to avoid a flexural failure. These large reinforcement ratios are not commonly used in practice. Considering the range of aggregate sizes available in the database, we can observe that both mortars and concretes are used for the experiments in the literature. The concrete compressive strength range shows that mixes from low strength to ultra-high-strength concrete are used.

The fiber types used in the experiments reflect all commercially available fiber types (hooked-ended, corrugated, crimped, straight smooth, round), as well as some fiber types that were explored for research purposes (flat-ended, flat, chopped with butt ends, straight mild steel, mill-cut, recycled fibers, brass-coated high-strength steel), and mixes of different types (hooked-ended and straight). The bond factor ρ_f in the expression for the fiber factor F takes into account the effect of the fiber type. For the less commonly used fiber types, there may be discussion about which value to use for the bond factor. Most of the experiments (63% of all experiments) analyzed used hooked-end fibers. Most specimens use a fiber factor of 0.5–1; higher values result in concrete mixes with low workability.

Since the proposed matrix-based model using artificial neural networks is only as good as the input data used for the model, the reader should keep the aforementioned limitations regarding the parameters used in the experiments in mind when applying the resulting expression for the design of members with SFRC. Extrapolation of the proposed model outside of the ranges of data points considered in this study does not guarantee a good approximation. The dataset used for this research, as well as the calculated values with our proposed model, can be found online [107].

2.2. Artificial Neural Networks

2.2.1. Introduction

The neural net, see Figure 3, consists of L layers of several nodes, with the first layer the input layer, layers 2 to L-1 the hidden layers, and layer L the output layer. The ANNs in this work are feed-forward: neurons connect to nodes in the layers further down the net, see Figure 3. Each node, except those in the input layer, has the following unknowns associated with it: the (non)linear transfer function, the synaptic weights W, and the bias b. The transfer function is determined by trying out different possible functions. W and b are determined through learning: finding a (local) minimum solution so that pre-determined requirements for performance of the neural net are met. Learning has three stages: training, validation, and testing. The input dataset is subdivided into a training dataset (for determining the unknowns of the neural net), validation dataset (for checking the generalization performance—loss of this characteristic is called overfitting, which can be caused by the net learning properties of the noise on the available data [108]), and testing dataset (independent verification of performance of neural net).

2.2.2. Implemented ANN Features

This work considers 15 ANN features, including data pre/post-processing features. We used parametric analysis (using nine sub-analyses) to determine the features that best describe the problem at hand.

Table 3. ANN features 1 through 5.

Feature Method	F1	F2	F3	F4	F5
	Qualitative Var Represent	Dimensional Analysis	Input Dimensionality Reduction	% Train-Valid-Test	Input Normalization
1	Boolean Vectors	Yes	Linear Correlation	80-10-10	Linear Max Abs
2	Eq Spaced in [0, 1]	No	Auto-Encoder	70-15-15	Linear [0, 1]
3	-	-	-	60-20-20	Linear [−1, 1]
4	-	-	Ortho Rand Proj	50-25-25	Nonlinear
5	-	-	Sparse Rand Proj	-	Lin Mean Std
6	-	-	No	-	No

Abbreviations: MLPN = multi-layer perceptron net, RBFN = radial basis function net.

,

Table 4. ANN features 6 through 10.

Feature Method	F6	F7	F8	F9	F10
	Output Transfer	Output Normalization	Net Architecture	Hidden Layers	Connectivity
1	Logistic	Lin [a, b] = 0.7[φmin, φmax]	MLPN	1 HL	Adjacent Layers
2	-	Lin [a, b] = 0.6[φmin, φmax]	RBFN	2 HL	Adj Layers + In-Out
3	Hyperbolic Tang	Lin [a, b] = 0.5[φmin, φmax]	-	3 HL	Fully Connected
4	-	Linear Mean Std	-	-	-
5	Bilinear	No	-	-	-
6	Compet	-	-	-	-
7	Identity	-	-	-	-

Abbreviations: MLPN = multi-layer perceptron net, RBFN = radial basis function net.

and

Table 5. ANN features 11 through 15.

Feature Method	F11	F12	F13	F14	F15
	Hidden Transfer	Parameter Initialization	Learning Algorithm	Performance Improvement	Training Mode
1	Logistic	Midpoint (W) + Rands (b)	BP	-	Batch
2	Identity-Logistic	Rands	BPA	-	Mini-Batch
3	Hyperbolic Tang	Randnc (W) + Rands (b)	LM	-	Online
4	Bipolar	Randnr (W) + Rands (b)	ELM	-	-
5	Bilinear	Randsmall	mb ELM	-	-
6	Positive Sat Linear	Rand [−Δ, Δ]	I ELM	-	-
7	Sinusoid	SVD	CI ELM	-	-
8	Thin-Plate Spline	MB SVD	-	-	-
9	Gaussian	-	-	-	-
10	Multiquadratic	-	-	-	-
11	Radbas	-	-	-	-

Abbreviations: SVD = singular value decomposition, MB SVD = mini-batch SVD, BP = back propagation, BPA = back propagation with adaptive learning rate, LM = Levenberg-Marquardt, ELM = extreme learning machine, mb ELM = mini-batch ELM, I ELM = incremental ELM, CI ELM = convex incremental ELM, NNC = neural network composite.

give an overview of the different formats used for the 15 features, based on information in the literature. Previous work (e.g., [109]) contains full descriptions and references to the literature used to select the different methods for the feature. MATLAB [110] is used for programming these routines. The neural network toolbox of MATLAB is used for the most commonly used algorithms for learning (1–3 in Table 5). In each sub-analysis (SA) the software runs all possible combinations of neural nets for preselected approaches for the features. The output is then the performance of each trial net. The optimal net is then the net with the best performance: the net with the smallest average relative error. In addition to average relative error, we also evaluate maximum error, and percentage of errors larger than 3%. The definitions used to assess the performance of the resulting neural net are defined in [109]. The developed software has been validated with several benchmark datasets and functions—a full validation report is available in the public domain [111].

With respect to the ANN formulation used in [109], two fewer changes were carried out for this work: (1) the elimination of performance improvements (Feature 14), and (2) the algorithm used in Feature 4. For the current study, four distributions of data p_t-p_v-p_tt (percentage for training, validation, and testing) were implemented (Feature 4). The following algorithm was implemented to divide the dataset into the training, validation, and testing subsets:

• Reduce p_t-p_v-p_tt values by 10 units each.
• Compute minimum and maximum values for each variable q (row) of the full input dataset.
• Define patterns where each variable takes its minimum or maximum value from the full input dataset. These patterns ought to be included in the training dataset. If the number of patterns is lower than p_t * P/100 (rounded off), more patterns should be added to the training set in the following way:

  (a)

  Compute the number of patterns (Lp_t) that need to be added to the initially selected training patterns to equal round (p_t * P/100).

  (b)

  Randomly select 10.000 combinations of Lp_t patterns from all those not included in the training set defined prior to (a).

  (c)

  For each combination/scenario in (b), add those Lp_t patterns to the set of training patterns defined prior to (a), and label all remaining learning patterns as “validation + testing”.

  (d)

  For each scenario in (c), and for each pattern labeled as “validation + testing”, check if that pattern has at least one input variable that equals a value not included in any pattern in the training set. If it hasn’t, then that pattern should be moved to the training set.

  (e)

  Among all 10,000 scenarios of training and “validation + testing” subsets addressed in (b) till (d), the selected scenario should be the one guaranteeing the amount of training data (P_t*) closest to round (p_t * P/100).

• If the training set selected in (e) guarantees |P_t*/P − p_t| ≤ 0.2, then that becomes the training data to be taken for simulation. Otherwise, the training data should be selected according to [112].
• Increase p_t-p_v-p_tt values by 10 units each (to re-obtain the original input values—See step 1).
• Randomly select p_v/(p_v + p_tt) of those patterns not belonging to the training dataset for the validation patterns. The remaining data then forms the testing dataset.

The distribution p_t-p_v-p_tt in the simulation can differ from the one chosen a priori (before step 1).

2.2.3. Parametric Analysis Results

The software runs nine SAs, of preselected ANN features. Feature 7 takes a single value only in each of the SAs. The SAs serve the purpose of determining the best method for the features in a consecutive way. Further information on how the SAs are defined can be found in [109]. 475 combinations of features were explored through the SAs.

Table 6. For best combo of each SA: F methods used.

SA	F1	F2	F3	F4	F5	F6	F7	F8	F9	F10	F11	F12	F13	F15
1	1	2	6	2	5	7	1	1	1	1	3	2	3	3
2	1	2	6	2	3	7	1	1	1	1	3	2	5	3
3	1	2	1	1	5	3	1	1	1	1	3	2	3	3
4	1	2	6	2	5	1	2	1	1	1	3	2	3	3
5	1	2	6	3	5	1	3	1	1	1	3	2	3	3
6	1	2	6	3	5	7	4	1	1	1	3	2	3	3
7	1	2	6	4	5	7	5	1	1	1	3	2	3	3
8	1	2	6	4	5	7	5	1	1	1	1	5	3	3
9	1	2	6	4	5	7	5	1	3	3	1	5	3	3

shows the best feature methods (that led to the combination with the best performance) from the different SAs. The numbers refer to those given in Table 3 through Table 5.

Table 7. For each SA: performance results of best net.

SA	ANN
	Max Error (%)	Performance all Data (%)	Errors > 3% (%)	Total Hidden Nodes	Running Time/Data Point (s)
1	24.2	0.7	5.6	36	1.88 × 10−4
2	1375.2	21.6	83.7	120	9.96 × 10−5
3	15.4	0.5	4.0	36	1.31 × 10−4
4	11.7	0.5	4.0	36	1.14 × 10−4
5	15.9	0.7	7.0	36	1.06 × 10−4
6	12.7	0.5	3.0	36	9.58 × 10−5
7	67.0	5.3	40.0	36	1.07 × 10−4
8	90.0	4.0	24.0	36	1.10 × 10−4
9	0.0	0.0	0.0	36	9.72 × 10−5

shows the performance results of the best combos of each SA. The results are obtained in the original format, compared to the actual values of the dataset. The microprocessor used in this work has the following features: RAM: 48 GB, OS: Win10Home 64bits, CPU: Intel^® Core™ i7 8700K @ 3.70-4.70 GHz, Local Disk Memory: 1 TB.

3. Results

3.1. Proposed ANN-Based Model

The proposed model is the one with the best performance among the best nets of each SA. As can be seen in Table 7, the best performance is obtained for the best net of SA9. For any user to apply our proposed model, we will provide all relevant expressions in the following subparagraphs and have provided the W and b arrays for download from the public domain. The proposed model uses five layers, and has the following distribution of nodes per layer: 9 in the input layer, 12 in each hidden layer, 1 in the output layer, see Figure 4, which also shows the connectivity of the network. The performance results of the proposed model are detailed in §3.2. The calculated solution with our proposed model is also available in the public domain for easy comparison [113].

3.1.1. Preprocessing of Input Data

Features 2, 3, and 5 deal with preprocessing of the data. We found that dimensional analysis nor reduction were necessary. The expression for input normalization is given as follows, with $Y_{1,sim}$ the input data:

$$\begin{array}{l}{\left\{Y_{1,sim}\right\}}_n^{after}=\left(\left\{Y_{1,sim}\right\}-\mathrm{INP}(:,1)\right)./\mathrm{INP}(:,2)\\ \mathrm{INP}=\left[\begin{array}{cc}149.674651162791& 66.7396705241561\\ 262.705930232558& 153.595523149765\\ 2.53778427906977& 0.958947876821730\\ 0.0244412479069768& 0.0104087109797821\\ 480.789007209302& 90.8692510105096\\ 11.0945581395349& 4.95663746228321\\ 49.7411932558140& 26.2694346953084\\ 0.555063267441860& 0.364880290571624\\ 1261.49069767442& 476.799170124293\end{array}\right]\end{array}$$

(58)

3.1.2. ANN-Based Analytical Model

The next part is the actual description of the ANN-based expression of our proposed model. They are presented here completely, with the W and b arrays in the public domain [113], so that all our results are reproducible. The following equations transfer the preprocessed input ${\left\{Y_{1,sim}\right\}}_n^{after}$ to the outputs of the hidden layers 2 through 4 and then to the preprocessed output ${\left\{Y_{5,sim}\right\}}_n^{after}$

$$\begin{array}{l}{Y}_2={\phi }_2\left(W_{1-2}^T{\left\{Y_{1,sim}\right\}}_n^{after}+b_2\right)\\ Y_3={\phi }_3\left(W_{1-3}^T{\left\{Y_{1,sim}\right\}}_n^{after}+W_{2-3}^T{Y}_2+b_3\right)\\ Y_4={\phi }_4\left(W_{1-4}^T{\left\{Y_{1,sim}\right\}}_n^{after}+W_{2-4}^T{Y}_2+W_{3-4}^T{Y}_3+b_4\right)\\ {\left\{Y_{5,sim}\right\}}_n^{after}={\phi }_5\left(W_{1-5}^T{\left\{Y_{1,sim}\right\}}_n^{after}+W_{2-5}^T{Y}_2+W_{3-5}^T{Y}_3+W_{4-5}^T{Y}_4+b_5\right)\end{array}$$

(59)

The following transfer functions are used:

$$\begin{array}{l}{\phi }_2(s)={\phi }_3(s)={\phi }_4(s)=\frac{1}{1+e^{-s}}\\ {\phi }_5(s)=s\end{array}$$

(60)

3.1.3. Output Data Post-Processing

Since in our proposed model, no output normalization or dimensional analysis is used, the output from Equation (59) {Y_5,sim}_n^after is the final result $Y_{5, sim}$.

3.2. Performance Indicators of Results

The performance of the proposed model is presented in this subsection. Figure 5 shows the regression plot of the relation between output of the dataset and network targets together with the Pearson Correlation Coefficient R. Figure 6 shows the average error values of the training, validation, and testing datasets, as well as of all data. Figure 7 shows the maximum error of all data as well as the percentage of data with an error above 3%. The reader should recall here that the outcome of repeat tests was averaged, so that some of the inherent material heterogeneity is removed to obtain unique datapoints.

3.3. Comparison between ANN-Based and Existing Methods

In this section, we compare the ANN-based model and the existing methods introduced in Table 1. The results of the ANN-based model are reported in [107]. The comparison between the experimental results and the existing methods are repeated here from [114]. Figure 8 shows the comparison between the existing methods proposed in the literature and our proposed model and the experimental results from the database. Figure 9 shows the comparison based on the existing code models.

Table 8. Statistical properties of V_utot/V_pred for all datapoints, with AVG = average of V_utot/V_pred, STD = standard deviation on V_utot/V_pred, and COV = coefficient of variation of V_utot/V_pred. This table repeats results from [114] for comparison to our proposed model.

Model	AVG	STD	COV	Min	Max
Proposed model	1.00	1.08 × 10−15	1.08 × 10−15	1.00	1.00
Sarveghadi et al. [22]	1.03	0.29	28%	0.23	2.49
Kwak et al. [23]	1.01	0.28	27%	0.27	2.39
Greenough and Nehdi [24]	1.34	0.48	36%	0.31	3.11
Khuntia et al. [25]	1.81	0.85	47%	0.18	6.53
Imam et al. [30]	0.97	0.36	37%	0.06	2.51
Sharma [26]	1.24	0.49	39%	0.18	3.59
Mansur et al. [27]	1.30	0.60	46%	0.15	3.85
Ashour et al. [28] 1	1.08	0.38	35%	0.24	3.14
Ashour et al. [28] 2	1.29	0.37	29%	0.31	3.22
Arslan et al. [29]	1.17	0.37	31%	0.43	3.24
Yakoub [31] 1	1.90	0.76	40%	0.28	7.50
Yakoub [31] 2	2.97	1.37	46%	0.51	17.48
French code [32]	1.85	0.88	48%	0.22	5.95
German code [33]	1.12	0.31	27%	0.21	2.13
fib [35]	1.24	0.36	29%	0.30	2.33
RILEM [34]	1.16	0.33	29%	0.23	2.28

gives the statistical results of V_utot/V_pred for all methods from Table 1 as well as our proposed model. The statistical properties of V_utot/V_pred result from all experimental results, and thus cover experiments on beams with a short shear span and slender beams. In [20], the analysis is further subdivided to evaluate the existing expressions for only slender beams, since expressions that did not include the enhancement factor for short shear spans will result in overly conservative predictions for beams with short shear spans. However, the overall conclusions regarding scatter on the results remains as discussed here.

From the presented results, we can conclude that our proposed model is a significant improvement as compared to the existing methods for determining the shear capacity of SFRC concrete members without shear reinforcement, for the dataset used in this study. The reason for this improvement is that our proposed model uses the available information from the literature in an optimal way. Moreover, comparing the predictions with our model to the model by Greenough and Nehdi [24] and Sarveghadi et al. [22], which were also based on soft computing methods, shows that using a larger database and evaluating a large number of ANN features results in a significant better fit of the experimental data.

4. Discussion

An important step is developing the proposed ANN-based expression was the selection of variables. While the initial data collection [20] focused on gathering as much information from the experiments as possible, it is not desirable to use all possible input values to develop the ANN-based expressions. Doing so increases the computational time of the algorithm that evaluates all 15 ANN features to find the optimal neural net. To select input variables, we included a number of dimensionless values, such as the clear shear span to depth ratio a_v/d, the reinforcement ratio ρ, and the fiber factor F. These inputs can take into account the combined effect of different parameters, and are more widely applicable since they are dimensionless. The factor a_v/d is the clear span to depth ratio, which takes into account the enhancement of the shear capacity due to direct load transfer for loads close to the support [115,116]. In addition, the clear shear span a_v, which is taken from the face of the loading plate to the face of the support, takes into account the dimensions of the support and the loading plate [117,118]. Moreover, from earlier shear research it was concluded that the distance that should be considered in the shear span to depth ratio is the clear shear span a_v and not the shear span a [119].

The reinforcement ratio ρ is the dimensionless equivalent of the area of longitudinal steel A_s. As a dimensionless parameter it allows to include a wider range of possible reinforcement layouts. The influence of the longitudinal steel needs to be considered as an input for a model for the shear capacity of SFRC, because the longitudinal steel resists shear through dowel action [120,121,122,123,124].

The fiber factor F [125,126,127,128] takes into account the geometry of the fibers (length and diameter) [53,129], the amount of fibers (fiber volume fraction) [69,70,130], and the bond properties of the fibers, which depends on the fiber type [4,131,132,133]. A challenge here was to ascribe bond properties to the less common fiber types that were encountered in the literature. Especially from older references [80,81,86], in which researchers were experimenting with many different fiber types, the fiber bond had to be estimated.

As mentioned before, the proposed model is limited to the ranges of the variables used for developing the model. These ranges are given in Table 2. While we can see that the experiments considered in the database cover a wide range for all material properties of the concrete [96,134], mild steel reinforcement [74,75], and fibers [135], the range of heights of specimens tested in the literature (expressed in Table 2 based on the effective depth d) is limited. We can see that the maximum effective depth in the experiments is 1118 mm. The number of specimens with an effective depth of d ≥ 600 mm for the considered dataset is 20 experiments on a total of 430 datapoints. As such, only 4.7% of the data from the literature considers realistic and large sizes of beams. All other specimens are laboratory-sized specimens. This finding underlines the need for further experiments on large SFRC beams failing in shear. The need for large-sized beams is especially important given the fact that the size effect in shear affects the shear capacity of large members. The size effect in shear [104,105,106,136,137] is the effect observed in experiments where the shear stress at failure decreases as the member height increases. The cause of the size effect is still under discussion. For SFRC, there is a discussion on whether or not the size effect in shear occurs as well [2,138]. Since the number of experiments on large beams is limited, and the size effect in shear in SFRC may be related to the amount of fibers (or fiber factor), further research is needed on this topic. Given the lack of data in this regard, we may need to reevaluate our proposed model when further experimental results on larger-sized beams become available.

Since the proposed ANN-based model simulates the experimental results with a very good fit, parameter studies with the ANN-based model give the same result as parameter studies with the original data. Such parameter studies were carried out earlier [20]. The outcome of these parameter studies can be summarized as follows:

• the shear strength strongly depends on the longitudinal reinforcement ratio [139], as a result of the larger contribution of dowel action [64,120,123,124,140,141] for larger amounts of reinforcement,
• for loads close to the support, the capacity increases as a result of direct load transfer [115,116,142],
• the shear strength strongly depends on the fiber factor, as a result of the additional tensile capacity across the crack for an increasing value of the fiber factor [5,25,126,128,143,144],
• the influence of the size of the aggregates on the shear capacity is relatively small, but should not negligible. Larger aggregates reduce the shear capacity of SFRC elements, as the mix becomes less uniform and the bond between the matrix and fibers becomes less [14,145].

With our proposed model, which can simulate well the available experimental results, the following practical applications are possible: the model can be used to prepare further laboratory studies, the model can be used for comparison to mechanical models, and the model can be used to predict the shear capacity of structural elements for design purposes, provided that the input parameters are within the specified ranges from Table 2. One of the main advantages of the proposed model as compared to other numerical approaches is the short amount of time required to obtain the prediction of the shear capacity. As can be seen in Table 7, the computation time is less than 0.1 millisecond. Therefore, compared to nonlinear finite element models, the proposed method is very fast.

At this moment, a theory that predicts the shear capacity of SFRC elements without mild steel shear reinforcement based on the different shear-carrying mechanisms is not available yet. Given the lack of theoretical understanding of the problem under study, our ANN-based proposed model can bridge the gap and give an optimal prediction based on the available data. However, research on the shear-carrying mechanisms in SFRC is still necessary to understand the actual mechanics of the problem.

5. Conclusions

This paper proposes a neural networks-based method to predict the shear capacity of SFRC elements without mild steel shear reinforcement based on the available data in the literature. To derive the ANN-based expression, we did the following:

• We used a database with 430 datapoints from the literature.
• For the analysis, we selected nine input parameters related to the geometry, properties of the concrete, the flexural steel reinforcement, and the fibers, and one output parameter, the maximum sectional shear force caused by the applied load in the experiment and self-weight of the beam.
• To find the optimal ANN-based model, different combinations of 15 features of ANN models were analyzed.
• The optimal model resulted in a maximum error of 0% and a mean relative error of 0.0% for the 430 datapoints, respectively.

The main advantages of our proposed model and the main outcomes of this study are:

• Our proposed model outperforms the available models and expressions for the shear capacity of SFRC.
• Our model can be used to prepare experiments, for design (within the input parameter ranges), and to support further development of mechanical models through robust parameter studies.
• The computational time of a datapoint with our model is less than 0.1 millisecond.

The limitations of the proposed model are as follows:

• The proposed model can only be used for the ranges of the variables available in the dataset.
• The model does not cover large-sized beams as a result of a lack of data on such specimens. As such, we recommend further experiments on large SFRC beams failing in shear and further studies on the size effect in SFRC.
• This study does not answer the question about the mechanics underlying the problem of shear in SFRC, but we can explore various influences with parametric studies using our proposed ANN-based model. Our model also facilitates the evaluation and improvement of existing and future mechanical models, based on the currently available experimental results.