The impact of cables on local scouring of bridge piers using experimental study and ANN, ANFIS algorithms

The main purpose of this study was to analyze and predict scour depth and hydraulic performance of piers using soft computing methods and to estimate scour depths using artificial neural networks and ANFIS methods. In the present study, three situations were studied: a rectangular pier without a cable (type I), a rectangular pier with a cable that has a diameter equal to 10% of the pier diameter, and a cable twist angle of 15 degrees (type II), and a rectangular pier with a cable 15% of the pier diameter and an angle of twist of 12 degrees (type III). Tests were carried out with different flow-approach angles: zero, 5, 10, and 15 degrees. Dimensional analysis based on the π Buckingham method was performed. Then, the effect of different parameters, and their importance for estimating scour depth was investigated. Piers with an angle of 15 degrees with respect to the direction of flow had the greatest depth of scouring. Cables can reduce scour depths at this angle; for the second and third classes of piers, the scouring is 10 and 22% compared to the first pier classification. For the second type of pier, angles of 5, 10, and 15 degrees led to increases in scouring depths of 3, 21, and 37% compared to the zero-angle situation.


INTRODUCTION
A major cause of bridge failure is the lack of attention paid to hydraulic standards in their design, implementation, and maintenance. To improve our understanding of bridge failure, a hydraulic study of piers reinforced with cables was performed with a focus on the scouring depths that can occur. There are various ways scouring around piers can be reduced. The use of slots, riprap, collars, submerged plates, and modern methods such as cables and nano clays have been studied for their ability to reduce scouring of sediments (Rezaie et al. 2018).
These issues have been studied in prior research. For example, Kumar et al. (1999), Osroush et al. (2019), and Saleh et al. (2020) studied the performance of a slot with changes to its shape and size, the arrangement and placement height of the slot, and piers with different shapes. Chiew (1995), Lauchlan & Melville (2001), Zarrati et al. (2010), and Froehlich (2013) investigated the stability of riprap around pier foundations and the temporal changes in geometric and hydraulic parameters. They protected the area around the pier from scouring and reduced scouring. maximum scour depth and the GA method showed more accurate results than MLR. El-Hady Rady et al. (2020) predicted the scour depth around a bridge pier using artificial intelligence that was based on conventional regression modeling. Their results showed that the gene programming model is superior to the neuro-fuzzy inference system model. Norouzi et al. (2019) used multilayer perceptron (MLP) networks, radial basis function networks, and support vector machines with different kernel functions to estimate the discharge coefficient labyrinth weirs. The results showed that the MLP model had a much higher accuracy than other models. Daneshfaraz et al. (2021aDaneshfaraz et al. ( , 2021b investigated the hydraulic performance of vertical drop using a support vector machine, therefore they analyzed hydraulic parameters from 120 different experimental data and 12 support vector machine models.
In previous studies, the phenomenon of scouring near bridge piers has been studied, however the effects of cable supports have not been considered. In addition, while the prediction of scour depth has been investigated previously, the present investigation goes further. Here we wrap the cable around the pier and change the collision angle of the flow to alter the intensity of the downward flow. Other objectives of the research include the following: 1. Investigating the effect of changing the diameter and torsion angle of the cable. 2. The effect of changing the collision angle of the flow to protect the pier and to reduce the scouring depth and also to increase the time required for maximum scouring depth. 3. The use of ANN algorithms and ANFIS to estimate scour depths and then the subsequent comparison of the performance of an artificial neural network with ANFIS. 4. A comparison of the results of the ANNs and ANFIS with experimental results.

Dimensional analysis
The effective parameters governing scour depth around a bridge pier outfitted with a cable and at a different angle of impact are as follows: F 1 (B, L, P, V, V c , D 50 , d s , t, t e , g, y, r, r s , y, a, u) ¼ 0 The symbol B is the pier width, L refers to pier length, P denotes the cable diameter, V is the flow velocity, V c represents velocity threshold conditions, D 50 is the media diameter size, d s is the final scour depth, t refers to the scour start time, t e is the scour equilibrium time, g is gravitational acceleration, y refers to the flow depth, ρ is water density, ρ s is the sediment particle density, υ is kinematic viscosity, α is the flow collision angle, and θ is the angle of twist of the cable. Using Buckingham's theory of dimensional analysis, Equation (1) can be written in dimensionless format as: Due to the constant values of the parameters gy/V 2 , ρ s /ρ, V/V c , L/B, d 50 /B, and y/B, and the turbulence of the flow from the above parameters and the Reynolds number υ/ρVD, certain simplifications are possible. By simplifying and combining dimensionless parameters, Equation (2) can be rewritten as: In the above relation, d s /B is the final dimensionless scour depth, t/t e dimensionless scouring time, P/B is the dimensionless cable diameter, α refers to the collision angle, and θ is the cable torsion angle.

Experimental methodologies
The experiments were performed in a rectangular experimental flume with a total length of 14 m, a useful length of 10 m, a width of 0.4 m, and a depth of 0.6 m with a fixed horizontal floor without a slope. The floor of the flume was made of PVC and the walls were made of Plexiglas with a thickness of 10 mm. There is an adjustable gate downstream of the test site that is able to control the water depth in the channel. A constant flow rate of 23 liters per second was regulated by a calibrated valve next to the pump. The range of the experimental portion of the flume is 2 meters long and 0.2 meters wide. To measure the scouring depth, a depth gauge was used. The gauge was movable in all three directions. Table 1 shows the criteria used in the present study.
Sand with a relative density of 2.65 was used to create the sediment. To ensure the uniformity of sediments, these materials have been passed through two consecutive sieves with 0.59 mm and 0.85 mm pore sizes. Residual sediments between the two sieves were selected for this experiment. These sediments were larger than 0.59 mm and smaller than 0.85 mm in diameter. In the present study, the sediments have d 50 ¼ 0.7 mm and σ ¼ 1.23. Raudkivi & Ettema (1983) found timewise changes in scour depth over three consecutive hours to be not more than one millimeter. In the present study, more than 90% of scouring depth occurs in the first 8 hours, so each experiment was limited to an 8-hour duration.

Physical models and experiments
Rectangular piers made of Teflon were used in the experimental fabrication. The width and length of the rectangular pier were 30 and 90 mm, respectively. The experiments were performed according to Figure 1, and the piers were divided into three groups according to Table 2.

Artificial neural network
An ANN is a very simple approach that was developed to model the biological structures found in the human brain (Kim et al. 2001;Pourtaghi & Lotfollahi-Yaghin 2013).  Raudkivi & Ettema (1983) The ratio of channel width to pier width (W/B) 6.25 7.65 Due to the lack of effect of the channel walls Melville & Sutherland (1988) The ratio of pier width to average grains diameter (b/d 50 ) Larger than 25 to 25 42.85 Due to the elimination of the effect of sediment grains size on scouring depth Melville & Sutherland (1988) The geometric standard deviation of sediment grains (σ ¼ d 84 . 1 /d 15.9 ) Less than 1.3 1.23 Preventing the non-uniformity of sediment grains reduces the depth of scouring equilibrium Chiew (1995), Ettema (1980) and Shen et al. (1969) Flow depth (y) More than 4 times the pier diameter 21 To prevent the increase of scouring depth in shallow flows by increase the flow depth Raudkivi & Ettema (1983) The average diameter of grains (d 50 ) Greater than or equal to 0.7 mm 0.7 mm To prevent the formation of ripple phenomenon in non-sticky sediments Raudkivi & Ettema (1983) Flow velocity (Vc) 0.9Vc 0.9Vc To start local scouring in clearwater ANN has a multilayered structure that includes many of the key elements of interconnected processing to mimic brain neurons. The neural network contains computational elements called neurons that combine inputs and determine their weight. Because ANN stores data as a template in a set of processing elements adapted from connection weights, complex maps can be identified through distributed representation features (Kim et al. 2001). Figure 2 is a representation of the ANN architecture.
In the forward step, the weighted value of the input components is calculated as follows: where net j is the weighted sum of the j th neuron for the input received from the preceding layer with n neurons, w ij is the weight between the j th neuron and the i th neuron in the preceding layer; x i is the output of the i th neuron in the preceding layer. The output of the j th neuron (out j ) is calculated using a sigmoid function as follows: This type of layered ANN has been utilized to solve many complex problems (Kartam et al. 1997). Details of the learning algorithm can be found in Karayiannis & Venetsanopoulos (2013) and interested readers are directed there.
The adaptive neuro-fuzzy inference system The compatible ANFIS system consists of 5-layer networks and node connectors. Reasonably structured neuro-fuzzy systems are proportional to input data, membership, rules, and membership functions output. Figure 3 is a representation of neuralfuzzy network architecture with two inputs, one output, and two laws (Aghdam et al. 2017).

Determine the structure of the proposed models
For the Levenberh-Margaret neural network algorithm, the number of hidden layers, and the number of neurons in each layer is determined by trial and error. Typically, 8 to 12 neurons are selected. To validate the created neural network, 70% Corrected Proof of the data were selected for training, and the remaining 30% were maintained for validation. The output of the sigmoid excitation function used in the reverse diffusion algorithm is in the range between zero and one, Therefore, the output should be in the same range. Mapping was done using the following general rule: The values x max , and x min are the maximum and minimum values in the x i range, respectively. Network training and network repetition rates of 0.1 and 400 are used, respectively. As the mean square error (MSE) increases, the network training process automatically stops. An ANFIS model to study two membership functions with 100 reps for a variety of membership function is intended. Considering the effect of the membership function and the use of trial and error, the psigmf membership function provided the best results.

Evaluation criteria
To evaluate the performance accuracy of the proposed intelligent models, statistical parameters measuring the correlation between observed and computed values were identified. The selected metrics are: root mean square error (RMSE), mean absolute error (MAE), correlation coefficient (CC), and demonstration coefficient (DC). The relationships of these indices are given in Equations (7)-(10).
t k and y k are the observed results (experiments) and the output of the neural network, respectively. First, the data were normalized using a Gaussian normalization technique. In the present study, the data were normalized using Equation (6), where x min , x max , respectively, refer to the minimum and maximum data and x n is the normalized data. To investigate the ability of cables to reduce the final scour depth, scour depths in the presence and absence of cables were determined and compared. The final depth was obtained at the end of each experiment. In general, for the second class of piers, the final scour depth is 9.09% less than the first type pier. For the third type of pier, the final scour depth is 10.59 and 18.72% less than the second and first pier types.
Temporal development of scour depth for the first pier and the influence of angle The temporal development of scour depth for the first pier type and angles of zero, 5, 10, and 15 degrees are provided in Figure 4. Over time, the appearance of horseshoe vortices has also dropped dramatically. These vortices, which enhanced vertical water motions, are primarily responsible for moving and transporting particles. For the first type of pier, a comparison of the 10-and 15-degree cases reveals very small differences, particularly for the first 90 minutes. After 90 minutes, the 15-degree case has a greater slope than the other curves. Also, as the collision angle increases, the scouring depth will increase.
Increasing the collision angle of the flow makes the pier wider and ultimately increases the strength of the horseshoe vortex and increases scouring. For experiments on the first type of pier and at angles of 5, 10, and 15 degrees increases the scouring depth by 4, 20, and 35%, respectively compared to the zero-degree case.
Temporal development of scour depth for the second pier type and the influence of flow angle In Figure 5, as the angle varies from zero to 5°, there is not much effect on the scour depth. The two scour depth curves coincide from the initial moment until equilibrium. Conversely, for the 10-and 15-degree angles, the scouring experiences equal slopes for 90 minutes. Thereafter, 15 degree pier has a steeper slope. Approximately 3 hours after the start of scouring, they become parallel again. The second type of pier, with an angle of 5 degrees experiences a 12% depth reduction, compared to the first type. The second type pier at an angle of 10 degrees leads to a 10% reduction in final scour depth; a collision angle flow of 15 degrees has the greatest scouring depth. The effect of cables for reducing the scour depth for this angle and for the second type pier is 10%.
For angles of 5, 10, and 15 degrees and for the second class of pier, there is an increase in horseshoe vortex activity and a consequent depth increase of 3, 21, and 37%, respectively (compared to the 0-degree case).
Temporal development of scour depth for the third type of pier and the influence of flow angle Figure 6 shows photographs of the experimental facility and Figure 7 shows the temporal development of scouring. As seen in the figure, for the first and second pier types, the scouring depths at angles of 10 and 15 degrees and for d s /b ¼ 0.2 differ. The piers with angles of zero and 5 degrees are also almost identical.
For the third pier type, it is observed that at an angle of 5 degrees, the scouring depth is reduced 18% compared to the first pier type. At an angle of 10 degrees, the third class of pier has a 17% reduction in the final scour depth compared to the first pier class. All experiments at the angle of 15 degrees have the most scouring and the effect of cables for reducing scour depth Corrected Proof (15-degree angle and the third pier type) is 22%. Table 3 shows the reduction of the scour depth for the second and third pier types, compared to the first type of the pier. By comparing Figures 4, 5 and 7, it is concluded that increasing the cable the diameter will reduce the initial scouring depth and the maximum depth.
For the third type of pier, a cable with a diameter 15% that of the pier diameter is used and the cable has a smaller torsion angle, with angles of 5, 10, and 15 degrees compared to the pier, there is a depth increase of 3, 20, and 27%, compared to the   Corrected Proof zero-angle case. Consequently, the third pier type with a collision angle of 15 degrees has a lesser scour depth than the first and second pier types. It can be seen that the third pier leads to less total scouring. As evident in Figures 4, 5 and 7 the scour depth for the third pier with a zero-degree angle is reduced by ds/ b ¼ 1.3, 1.4, and 1.6, respectively. Such decreases can be seen from other angles as well. The final scour depth for the third pier is 17% less than for the first pier.
Temporal development of scour depth according to the type of pier Temporal development of scour depth for first, second, and third pier types at zero degrees Typically, increasing the collision angle of the flow makes the pier wider and increases the power of the horseshoe vortices. Consequently, as the collision angle increases, the depth of scouring also increases. For the second pier type, a thin cable was used. For the first 90 minutes of scouring, the results are similar to the first pier without a cable. The initial scour depth for the third pier is less than the scour for the second pier. Furthermore, the second pier experiences less scouring than the first pier. This finding demonstrates the influence of a cable in reducing scour depth.

Temporal development of scour depth for all three piers at a 5 degree angle
For the three pier types, it was observed that with a flow collision angle of 5 degrees, the scour depth is increased marginally, by 4, 3, and 3%, respectively (compared to the zero-angle situation). As can be seen from Figure 9, when scouring initiates at the third pier, the scour depth is much less than that for the first pier. Therefore, it can be concluded that increasing the diameter of the cable has a significant effect on reducing the depth of scour in the early moments. A comparison of these curves in Figure 9 shows the depth of scouring at this angle for the first type of pier. It is much more than for piers equipped with a cable. The third pier initially prevented the scoured hole from deepening. The final scour depth is closer to the second and third types of piers and is most different compared to the first pier. For all three piers, and for at t/te ¼ 0.2, about 90% of the scour depth has developed. Thereafter, the scouring process continues slowly. The results were obtained with a collision angle of the flow of zero and 5 degrees (Figures 8 and 9). The dimensionless depth (ds/B) for the first to third type piers has increased from 1.34, 1.45, and 1.6 to 1.47, 1.5, and 1.67, respectively.

Temporal development of scour depths for all pier types, at a 10-degree angle
For piers with a 10-degree collision angle, the first, second, and third piers experience a scouring depth increase of 20, 21, and 20% respectively, compared to the zero-degree case. This issue can be seen in Figure 10, when compared to Figures 8 and 9.

Corrected Proof
The second and third pier curves are similar until t/t e ¼ 0.1. Thereafter, the first pier experiences more scouring than other piers. Increasing the collision angle of the flow increases the width of the pier and intensifies the downward flow and horseshoe vortices, as already discussed. This leads to an enhancement of vertical flow caused by horseshoe vortices.

Temporal development of scour depths for all three piers, at a 15 degree angle
For all three piers, a cable with a diameter that is 15% of the pier diameter was employed. The cable has a low torsion angle. There are depth increases of 35%, 37%, and 27.5%, respectively, for the three piers, compared to the zero-flow collision angle    Figure 11 for the second and third pier types shows that the scour depth for piers one and two are very similar, until t/t e ¼ 0.2. A slight difference is observed for the second pier and is probably due to the insufficiency of the thin cable (with a diameter of 10% compared to the pier). Table 4 shows the increase of the final scour depth with collision angles relative to the zero-degree angle case.
To compare the present study with the experimental studies of previous researches, similar models in the same (laboratory) flow conditions should be used. Therefore, similar studies using cables, collars and gaps were used to compare with the present study. For this purpose, the best results of the present study and previous researches have been used in Table 5.

Prediction of scour depth by ANN and ANFIS methods
To predict the scour depth using ANN methods and ANFIS, the 11 models with input parameters listed in Tables 6 and 8 were used. From the 11 models. The optimal structure was determined based on performance criteria for both the learning and testing processes. Next, a membership function was adjusted using a recursive diffusion algorithm in combination with the recursive least squares method. To reduce the MAE, the number of membership functions is gradually increased by reducing the penetration amplitude of the cluster centers using a trial-and-error method. To compare ANN and ANFIS methods, the RMSE, MAE, correlation coefficient (CC), and coefficient of determination (DC) are utilized (see Tables 6 and 8).

Estimations of the final depth of scouring using the ANN method
The 11 models encompass different combinations of input variables to be able to provide the best combination for scouring estimation. The goal is to adjust the models so that the effect of single variables can be analyzed and compared with models that have other input variables. By reviewing the results and considering the evaluation criteria for different models, it was observed that model number 1 is superior compared to the others. This model has the lowest RMSE and MAE (close to zero) and the largest CC and DC (close to one). Table 6 has been prepared to help compare results. It is found that by considering the ratio of cable diameter to pier width (P/B) and comparing scour start time to scour equilibrium time (t/t e ), the effect of the flow collision angle parameter (α) is much more significant than the cable torsion angle parameter (θ). It should be noted that for predictions made with an ANN, models 2 and 4 differ from model 1 by less than 3% in the CC criterion. In other words, the cable diameter parameter to pier width (P/B) has very little effect on scouring.
The diagram in Figure 15 is obtained using the values in Table 6. Values of the performance criteria CC and RMSE are shown as a percentage; values of CC ∼ 1 and RMSE ∼ 0 signify better performance. According to Table 6, and as mentioned in Figure 12, it can be seen that model number 1 (the superior model) is the best, and model number 2 is next in performance  Corrected Proof quality. Therefore, these two models have the best accuracy. In addition, model 10 has minimal accuracy according to Table 6 and Figure 12, These observations show that the P/B and θ are unimportant for estimating the final scour depth. According to the performance criteria studied by the researchers in Table 7, RMSE and R 2 criteria are the most used compared to other criteria. By comparing RMSE and R 2 criteria in the study of Bateni et al. (2007) with 98% accuracy, Ismail et al. (2013) estimated the scour depth around the bridge pier with 83% accuracy. The present study with RMSE and MAE of 0.075 and 0.035, respectively, has good accuracy. Many studies have been done to estimate scour depth, which is of particular importance to designers because of the importance of scour depth around the bridge pier.

Results of final scour depth estimation by ANFIS method
In this study, 11 models were performed (see Table 8) to estimate the scouring depth with ANFIS. The results and evaluation criteria for different models showed that model 2 is superior. This model has the lowest RMSE and MAE (close to zero) and the highest CC and DC (close to one) compared to other models. According to the results of Table 8, and by comparing models 7 and 9 and also models 8 and 10 with constant values of parameters (t/t e ) and (P/B), the effect of the flow angle parameter (α) is found to be much more significant than the cable torsion angle parameter (θ).
The diagram in Figure 13, corresponds to the values in Table 8. The values of CC and RMSE are shown in the figure as a percentage. The greater the distance between the two lines indicates an improved performance. According to Table 8 and seen in Figure 13, model 2 (a superior model) has the lines with the greatest distance and the maximum accuracy. Model 10 has the lowest distance and minimum accuracy according to Table 8 and Figure 13.
According to Tables 6 and 8, the input parameters are the same in the methods used and in the proposed models. Therefore, the comparison of two methods with the same parameters in the input data has been made. Figure 14, shows a comparison of the top models of the neural network and ANFIS method, As can be seen, the ANN method (see Table 6 and Figure 14) (first   Laursen and Toch (1956) 0.8160 -1.0138 -- Shen (1971) 0.7119 -0.5994 -- Hancu (1971) 0.7303 -1.1360 -- Coleman (1971) 0.6786 -0.4765 -- Ettema et al. (1999) 0.7523 -1.5825 --BPN (Lee et al. 2007) 0  If it is necessary to use the proposed methods in practice, it is enough to present the maximum amount of error obtained from the method to the proposed value. The Table 9 is the results of validation data for researchers who used ANFIS to predict scour depth. Due to the large effect of the number of data and input variables on the inference and according to the fuzzy system computational algorithm, the results of this method, like the neural network method, are in an acceptable range, but it still has less adaptation than the neural network. Due to the fact that in this study only nine different experimental models have been investigated in different conditions, there is not suitable scope for adaptation to the ANFIS computational algorithm. For this reason, it has shown poorer performance in estimation, which can be remedied by performing more and various experiments.

Sensitivity analysis
To evaluate the effect of input variables on the output predictions, the neural network sensitivity analysis (superior method) was used. With this method, the inputs to a statistical model are changed in a systematic way so that the effects of these changes can be identified: By removing each parameter from the series of input parameters, the impact of the parameters could be assessed. The effect of parameter removal on the accuracy of the overall model was quantified using the statistical parameters RMSE, MAE, CC, and DC.

Corrected Proof
As can be seen from the results of the sensitivity analysis in Table 10, by removing the α parameter, the accuracy of the results was drastically reduced. The α parameter has the largest effect among the input parameters. the effect of this parameter compared to the total input parameters is a 40.31% reduction in accuracy. Data corresponding to Table 10 are shown in Figures 15 and 16. Also, the L/B parameter according to Table 10 has less sensitivity in estimating the final scour depth.

CONCLUSION
In this study, the performance and efficiency of cable in controlling local scour around the pier of a rectangular bridge by changing the flow collision angle were investigated and presented in experimental and intelligent algorithms. Based on the results, by increasing the cable diameter, the initial and final scouring depth can be reduced. The reason for this is to reduce the downward flow and weaken the horseshoe vortices. Using different experimental models, the scouring method and the final scouring depth reduction by cable are described. The general results are as follows:

1.
Wrapping the cable around the bridge pier (with diameters of 10 and 15% of the pier diameter) and angles of 15°and 12°w ill reduce the depth of the scour pit by up to 22%. 2. For the first, second, and third foundation types, scouring starts upstream of the bridge pier and after about 3 hours, it extends downstream of the pier at the nose. 3. At the beginning of scouring, the presence of horseshoe vortices on both sides of the pier deposit sediments downstream of the pier. After three hours, it reaches the initial depth and is washed. 4. For the second and third pier types, the rate of scouring is somewhat slower compared to the first pier type. 5. Reducing the angle of twist of the cable around the pier reduces the depth of scouring. If the cable twists are close together, the number of turns of available cables increases and this reduces the scouring depth. 6. By increasing the flow collision angle to the pier of the bridge, the scour depth increases. Changing the angle by more than 10 degrees is not recommended. 7. Analysis using ANN and ANFIS methods showed a high capability for estimating scour depth. It was also observed that the model with input parameters α, θ, t/t e , P/B and using a neural network method has the highest accuracy and efficiency with the highest correlation coefficient and coefficient of explanation (CC and DC) and the lowest error (RMSE). 8. The results of the sensitivity analysis showed that by removing the input parameter α, the model error rate increases to a large extent (40.31%). Therefore, this parameter is the most important variable for estimating scour depth.

DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.