SOIL-NAIL WALL STABILITY ANALYSIS USING ANFIS Primož jelušic and bojan žLeNDeR about the authors corresponding author Primož Jelušič University of Maribor, Faculty of Civil Engineering Smetanova 17, 2000 Maribor, Slovenia E-mail: primoz.jelusic@um.si Bojan Žlender University of Maribor, Faculty of Civil Engineering Smetanova 17, 2000 Maribor, Slovenia E-mail: bojan.zlender@um.si Abstract We present the safety-factor optimization for a soil-nail wall. The optimization is performed using the non-linear programming (NLP) approach. For this purpose, the NLP optimization model OPTINC was developed. The safety factor and the optimal inclination of the soil nails from the horizontal direction depend on the design of the soil-nail wall. Based on these results the ANFIS-INC model was developed for the prediction of the optimal inclination of the soil nail for any design of soil-nail wall. Additionally, an ANFIS-SF model was developed to predict the safety factor for different inclinations of the wall, the slope angle of the terrain, the length of the nails, and the hole diameter. It was found that the inclination of the soil nail should be adjusted to the inclination of wall, the length of nail, the slope angle of the terrain and the hole diameter. With increasing inclination of the wall, the length of the soil nail and the hole diameter, the safety factor is increasing. On the other hand, the safety factor is decreasing with the increasing slope angle of the terrain. The use of nonlinear programming and an Adaptive Network Based Fuzzy Inference System allows a comprehensive analysis of the geotechnical problems. Keywords soil-nail, wall stability, optimization, NLP, ANFIS introduction The idea of combining passive steel reinforcement and shotcrete was first used as a support system for underground excavations in rock. This system is known as the New Austrian Tunneling Method [1,2,3]. The soil-nailing support technique relies on the mobilization of the tensile strength of the steel reinforcement at relatively small deformations in the surrounding soil. The soil-nailing system was first used in 1972 to support an 18-m-high cut-slope near Versailles [4]. Since then, soil nailing was used in many locations in France [5,6], Germany [7,8,9] North America [10] and other countries. The use of soil-nail walls has substantially increased, because it has been demonstrated that soil-nail walls are technically feasible and, in many cases, a cost-effective alternative to the conventional retaining walls used in top-to-bottom excavations in temporary and permanent applications [11]. This paper presents the design analysis of a soil-nail wall. The internal stability of the soil-nail wall was calculated using a plane slip surface. Based on the analytical solutions for the nailed slopes an ANFIS-SF model was built. The depth of the excavation, the material properties of the soil, the number of nails, and the vertical and horizontal spacing of nails are constants, while the inclination of the wall, the slope angle of terrain, the length of the nails and the hole diameter are variables. The ANFIS-SF model predicts a safety factor for the initial design of the soil-nail wall. The safety factor is calculated for each angle (from 1° to 89°) of the slip surface. The critical angle of the slip surface is the one that has the smallest safety factor and depends on the inclination of the soil nails. Non-linear programming (NLP) was used to calculate the optimal inclination of the nails from the horizontal direction for different inclinations of the wall, the slope angle of the terrain, the length of the nails and the hole diameter. Based on the results obtained with NLP an ANFIS-INC model was built that allows a calculation of the optimal inclination of the soil nails for the initial design of the soil-nail wall. ACTA GEOTECHNICñ SLOVENICA, 2013/l 61. P. ÎELUSIS & b. slender: soil-nail wall stability analysis using anfis Adaptive network fuzzy inference systems (ANFIS) have been applied to many geotechnical engineering problems and have demonstrated some degree of success [12,13,14, 15,16,17]. For braced excavations in soft clay, Goh et al. [18] developed a neural network model to provide initial estimates of the maximum wall deflections. With an efficient mathematical technique two ANFIS models were developed to predict the safety factor (ANFIS-SF) and the optimal soil-nail inclination (ANFIS-INC). 2 soil nailing Soil nailing is a construction technique that has many advantages over conventional retaining walls used in top-to-bottom excavations. However, passive anchors are mainly used for the temporary protection of the excavation, and where the displacement of retaining structures are permitted. The soil-nailing technique is carried out in stages: 1. excavation, 2. installation of the nails, 3. installation of the steel reinforcing bars or mesh and drainage, 4. shotcreting, 5. installation of the bearing plates. These stages are repeated until the bottom of the excavation. The typical spacing between the soil nails is 1-2 m in the vertical and horizontal directions. The inclination of the soil nails is 5-30° from the horizontal direction. The holes are either cased or uncased depending on the type of soil and are, on average, 50 to 150 mm in diameter [11,19]. Soil nailing is a cost-effective alternative to pre-stressed anchors [20]. Soil nails are cheaper, mainly because they are made from cheaper materials. The use of soil nails is not suitable when the soil properties are inappropriate and displacements of supporting wall are not permitted [21]. Suitable soil properties for the use of passive anchors are: 1. natural cohesive materials, such as silts and low plasticity clays not prone to creep, 2. glacial till, 3. cemented sand with little gravel, 4. fine-to-medium sand with silt to acts as a binder. 1. materials without cohesion, 2. loose granular soils, 3. silt and clay of high plasticity, 4. soft, cohesive soils that will not provide a high pullout resistance, 5. expansive clays, 6. the presence of groundwater. The corrosion protection in aggressive soils has to be examined very closely in the soil-nailing system. The standard EN 14490:2010 (Execution of special geotech-nical works - Soil nailing) establishes general principles for the execution, testing, supervision and monitoring of soil nailing. In aggressive soils, fully encapsulated nails are recommended. The ultimate limit state (ULS) and serviceability limit state (SLS) are considered in the design and analysis of soil-nail walls. Three ULSs must be verified: 1. the external failure mode, 2. the internal failure mode, 3. the facing failure mode. The major SLS is excessive wall deformation. According to Byrne [10], the external failure modes are those were failure surface do not intersect the nails and the internal failure modes intersect the nails. The external failure modes are: 1. global stability failure, 2. sliding stability failure, 3. bearing failure. The internal failure modes are: 4. nail-soil pullout failure, 5. bar-grout pullout failure, 6. nail tensile failure, 7. nail bending failure. The facing failure modes are: 8. facing flexure failure, 9. facing punching shear failure, 10. headed-stud failure. Soil-nailed walls can withstand large deformation in all directions; therefore, they perform well during an earthquake. When large wall deformations are not permitted, soil nailing is not an appropriate retaining wall system. Passive anchors are not recommended for the following soil conditions: 62. ACTA GeOTeCHNICA SLOVENICA, 2013/1 P. ÎELUSIS & b. slender: soil-nail wall stability analysis using anfis 3. plane slip surface model An internal stability of soil nail wall could be calculated with different types of slip surface, such as plane slip surface, broken slip surface, parabolic slip surface, circular slip surface and logarithmic spiral slip surface. In this paper, the plane slip surface was used to calculate the safety factor. The geometry and geo-mechanics of the model are presented. The model is designed so that it is possible to change the geometric and geo-mechanical parameters. After determining the geometric design, the safety factor is calculated for each angle (from 1° to 89°) of the slip surface 9. The critical angle of the slip surface is the one that has the smallest safety factor. 3.1 GEOMETRY OF THE SOIL-NAIL WALL The geometry of the soil-nail wall is shown in Fig. 1. In this model the depth of the excavation wall (H) and the spacings of the soil nails in the vertical direction (sv) and the horizontal direction (s^) were determined. The typical spacing between the anchors is 1-2 m and this is the same in the horizontal direction and the vertical direction. Based on these data it is possible to determine the number of anchors (No). The geometry is designed so that it is possible to change the inclination of the wall (a), the slope angle of the terrain (3), the length of the nails (l), the inclination of the soil nails (n), the spacing of the soil nails in the vertical direction (sv) and the horizontal direction (sp). 3.2 GEO-MECHANICAL MODEL OF SOIL-NAIL WALL In the soil-nailing system for each nail, three bearing capacities are calculated: 1. tensile strength of the nail, 2. pull-out nail-bearing capacity, 3. nail-cap bearing capacity. The strength characteristics of a nail represent the basic parameters to calculate the actual force in a nail. The tensile strength [11] of the nail is calculated with Eq. (1): Rit = dl f y 4 SFt (1) where: R is the strength against breaking (kN), dn is the nail diameter (mm), fy is the strength of the nail material (MPa), SFt is the safety factor against breaking. The pull-out resistance [11] is calculated with Eq. (2): p ■d q (2) tip = SF„ Figure 1. Internal stability analysis of a soil-nail wall using a plane slip surface. 62. ACTA GeOTeCHNICA SLOVENICA, 2013/1 P. ÎELUSIS & b. slender: soil-nail wall stability analysis using anfis where: Tip is the pull-out nail bearing capacity (kN/m), d is the hole diameter (mm), qs is the ultimate bond strength (kPa), SFp is the safety factor against pull-out. The nail head strength [11] is calculated with Eq. (3): R f = Min(Ru ip ■ l).((0.6 + 0.2-(smax -i))) 0) For the design, the tensile force distribution along the nail can be simplified, as shown in Fig. 2. The tensile force in the nail increases at a constant slope Tip (equal to the pullout capacity per unit length), reaches a maximum value, Rit , and then decreases at the rate Tip to the value Rf at the nail head. ■^U Rn R" where: R1f is the nail-cap bearing capacity (kN), l is the length of the soil nail (m), smax is the spacing of the soil nails in the horizontal or vertical direction, whichever is greater (m). In plane slip a specific slip surface is examined for a variation of the angle 9. In an optimization analysis the calculation is carried out for different angles of the slip surface (From 1° to 89°). The safety factor is calculated for each angle of the slip surface. The ratio of resisting and shear (driving) forces acting on a slip surface should be greater than the minimum safety factor. The forces acting on a slip surface are: 1. gravitational force parallel to the slip surface, 2. active earth pressure acting on the vertical part of the structure and parallel to the slip surface. The resisting forces are: 3. the soil friction and cohesion along the slip surface, 4. the sum of the forces transmitted by the nails. The nail force is determined based on the location of its intersection with the slip surface (see Fig. 1). If a nail is completely in front of the slip surface, then it does not enter the calculation. If a nail crosses the slip surface, then its force is determined with Eq. (4): F = Min(Ti p-x, Rit, R f + Tx p ■ y (4) where: Rif Rit T ip is the nail length behind the slip surface in the direction of the soil body (m), is the nail length in front of the slip surface (m), is the nail cap bearing capacity (kN), is the strength against breaking (kN), is the pull-out nail-bearing capacity (kN/m). I- Nail head Nail length Figure 2. Distribution of the tensile force along the nail. The safety factor [ii] against global failure SF is expressed as the ratio of the resisting and driving forces, which act at a tangent to the potential failure plane: SF = Fh ■ cos(J + h) + Fcd (G + Sa,sv )sin(J) + Sa,vod ■ cos( (5) Fh = E Fh,n (6) Fed =Et( cos(j) + Fh ■ sin(j + h))-tan(^, ) + "E di,si up ■ ci where: G is the gravitational force (kN/m), Sa,sv is the vertical component of active pressure (kN/m), Sa,vod is the horizontal component of active pressure (kN/m), di,siip is the length of the ith section slip surface (m), (7) dslip Fh,n ft s is the length of the slip surface (m), is the bearing capacity of the nth nail behind slip surface (kN/m), is the cohesion of the ith soil layer (kPa), is the angle of the internal friction of the ith layer (°), is the inclination of the slip surface (°), is the inclination of the nails from the horizontal direction (°). x 62. ACTA GeOTeCHNICA SLOVENICA, 2013/1 P. ÎELUSIS & b. slender: soil-nail wall stability analysis using anfis The active earth pressure (the Coulomb theory) is given by Eq. (8) = Sz ■ Ka - 2cef ■ Kac (8) where: az is the vertical geostatic stress (kPa), Ka is the coefficient of active earth pressure (-), Kac is the coefficient of active earth pressure due to cohesion (-), The horizontal and vertical components of the active earth pressure are: Sa,vod = sa ' cos(a + S) ~ (9) Sa,sv = S a ■ sin (a + S) H (10) where: S is the angle of the friction structure - soil (°). 3.3 OPTIMIZATION MODEL OPTINC As the optimization problem of the soil-nail wall is nonlinear, e.g., the objective function and the (in)equality constraints are non-linear, the non-linear programming (NLP) optimization approach is used and described in the paper. The retaining wall optimization using the NLP approach was presented by several authors [22, 23]. The general NLP optimization problem can be formulated as follows: Max z=f(x) subjected to: h(x)=0 g(x)<0 x e X = {x|x e Rn,xLo < x < x1 where x is a vector of the continuous variables defined within the compact set X. The functionsf(x), h(x) and g(x) are nonlinear functions involved in the objective function z, equality and inequality constraints, respectively. All the functionsf(x), h(x) and g(x) must be continuous and differentiable. In the context of structural optimization, the variables include dimensions, cross-section characteristics, strains, materials, stresses, etc. The equality and inequality constraints and the bounds of the variables represent a rigorous system UP of the design, loading, stress and resistance functions taken from the structural analysis and dimensioning. The optimization of the structures may include various objectives worthy of consideration. In this paper, a geo-mechanical objective function is proposed to maximize the safety factor of the soil-nail wall. 3.3.1 input data Input data represent the design data (constants) for the optimization. The design data (constants) comprise the height of the soil-nail wall H (m), the vertical sv (m) and the horizontal Sh (m) spacing of the nails, the number of nails No (-), the inclination of the wall a (°), the slope angle of the terrain ft (°), the length of the nails l (m) and the hole diameter d (mm). In addition, the geo-mechanical data comprise the soil properties, such as the unit weight y (kN/m3), the angle of the internal friction fef (°), the cohesion of the soil (kPa), the ultimate bond strength of the soil nails in the soil qs (kPa) and the structure-soil angle of friction S (°). 3.3.2 variable The inclination of the soil nails q (°) is declared as a variable in this optimization model (see Fig. 1). 3.3.3 Geo-mechanical objective function The objective function is defined with equation Eq. (11): SF = Max (SFmt) (11) The objective function SF includes the gravitational force G (kN/m), the vertical component of the active pressure Sa,sv (kN/m), the horizontal component of the active pressure Sa,vod (kN/m), the bearing capacity of the nails behind the slip surface Fh (kN/m), the resisting forces acting on a slip surface Fcd (kN/m). 3.3.4 Geo-mechnical equality constraint The critical angle of the slip surface 9crit (°) is the one that has the smallest safety factor and depends on the inclination of the soil nails. SFcrit = Min Fh ■ C0S(J + h) + Fcd G + Sa,sv ) Sin(J) + Sa v0d ■ C0s(J) (12) 3.3.5 Design (In)eQuallty constraint The design (in)equality constraint determines the inclination of the soil nails from the horizontal direction q (°) to be calculated inside the defined limits. Lo ^ ^ h £ h £ h UP (°) (13) 62. ACTA GeOTeCHNICA SLOVENICA, 2013/1 P. ÎELUSIS & b. slender: soil-nail wall stability analysis using anfis 3.3.6 Result of the optimization model In order to interpret the proposed optimization approach, the paper presents a numerical example. The soil properties of silty sand are given in Table 1. Table 1. Soil properties of silty sand. Unit weight Y kN/m3 18 Angle of internal friction fef O 30 Cohesion of soil cef kPa 5 Angle of friction structure - soil S O 20 Ultimate bond strength qs kPa 100 The design data (constants), the initial values, and the lower and upper bounds of the design variable are given in Table 2. Table 2. Design data, initial values, lower and upper bounds of the design variable (see Fig. 1). Height of wall H m 8 Vertical spacing sv m 1,5 Horizontal spacing sh m 1,5 Inclination of wall a 0 0 Slope angle of terrain P 0 0 Length of nail l m 6 Hole diameter d mm 100 Lower bound Initial value Upper bound Inclination of nails q (°) 5 10 35 Silty sand a = 0° ß = 0° H = 8 m 1=0.75 H (6m) d = 100 mm Tlopt 5 15 25 35 Incilination of nails n (°) Figure 3. Optimal inclination of the soil nails from the horizontal direction qopt. The value of the safety factor is changing with the inclination of the nails q The optimal inclination of the soil nails from the horizontal direction n0pt is defined with the maximum safety factor (see Fig. 3). In this numerical example the optimal inclination of the soil nails from the horizontal direction is q0pt = 19.1°. 4 adaptive-network-based fuzzy inference system The basic structure of the fuzzy inference system (FIS) was introduced by Zadeh [24]. In this type of FIS it is essential to predetermine the rule structure and the membership functions. Human-determined membership functions are subjective and are different from person to person. Standard methods that transform human knowledge or experience into fuzzy rules and membership functions do not exist. Usually, there is a collection of input/output data, which we would like to use to construct the FIS model. The effective method for tuning the membership functions and minimizing the output error measure is the Adaptive-Network-Based Fuzzy Inference System (ANFIS). The ANFIS [25] uses given input/output data to construct a FIS, whose membership function parameters are tuned (adjusted) using either a back-propagation algorithm alone or in combination with a least-squares type of method. This adjustment allows fuzzy systems to learn from the data they are modelling. ANFIS only supports Sugeno-Takagi-Kang [26] identification models, which should Table 3. Input and output data for two ANFIS models. ANFIS-SF ANFIS-INC Constants Height of wall H (m) 8 8 Vertical spacing sv (m) 1.5 1.5 Horizontal spacing s^ (m) 1.5 1.5 Sv,up (m) 1 1 sv,bottom (m) 1 1 Number of nails No (-) 5 5 Inclination of soil nail q (°) 10 - Input data Inclination of wall a (°) 0, 10, 20 0, 10, 20 Slope angle of terrain j (°) 0, 15, 30 0, 15, 30 Length of nail l (m) 4, 6, 8 4, 6, 8 Hole diameter d (m) 50, 100, 150 50, 100, 150 Output data Safety factor SF (-) Inclination of soil nail n (°) 62. ACTA GeOTeCHNICA SLOVENICA, 2013/1 p. ÎelusiS & b. slender: soil-nail wall stability analysis using anfis have only one output parameter. The adaptive network is a superset of all kinds of feedforward neural networks with a supervised learning capability [27]. ANFIS is a fuzzy inference system implemented in the framework of adaptive networks and uses the advantages of neural networks and fuzzy logic. The safety factor SF and the optimal inclination of the soil nails from the horizontal direction n0pt depend on the design of the soil-nail wall. Therefore, the optimal inclination was calculated with the optimization model OPTINC for a different inclination of the wall a (°), the slope angle of terrain ¡3 (°), the length of the nails l (m) and the hole diameter d (mm). Based on these results the ANFIS-INC model was developed for the prediction of the optimal inclination of the soil nail for any design of soil-nail wall. Additionally, the ANFIS-SF model was developed to predict the safety factor for different inclinations of the wall a (°), the slope angle of terrain 3 (°), the length of the nails l (m) and the hole diameter d (mm). 4.1 RESULTS OF GEO-MECHANICAL AND OPTIMIZATION MODEL The soil properties of the geo-mechanical model are given for silty sand (see Table 1). The design data (constants) are given in Table 3. The series of calculations for the safety factor SF (-) and the optimal inclination soil nails q (°) are given in Table 4. Table 4 Table 4. Output data for different combinations of input data. Input Output Input Output Input Output laß d SF nopt laß d SF nopt l a ß d SF nopt (m) (°) (°) (mm) (-) (°) (m) (°) (°) (mm) (-) (°) (m) (°) (°) (mm) (-) (°) 4 0 0 50 0.69 5.00 4 0 15 50 0.62 5.00 4 0 30 50 0.53 5.00 6 0 0 50 0.89 5.53 6 0 15 50 0.78 5.00 6 0 30 50 0.66 5.00 8 0 0 50 1.09 5.00 8 0 15 50 0.94 5.00 8 0 30 50 0.77 5.00 4 0 0 100 0.85 9.54 4 0 15 100 0.76 8.02 4 0 30 100 0.64 5.00 6 0 0 100 1.17 14.98 6 0 15 100 1.00 13.99 6 0 30 100 0.82 11.20 8 0 0 100 1.51 13.85 8 0 15 100 1.24 13.85 8 0 30 100 0.97 12.57 4 0 0 150 0.98 14.03 4 0 15 150 0.86 14.03 4 0 30 150 0.72 10.05 6 0 0 150 1.39 20.57 6 0 15 150 1.16 19.47 6 0 30 150 0.92 16.75 8 0 0 150 1.80 18.94 8 0 15 150 1.43 20.19 8 0 30 150 1.07 20.19 4 10 0 50 0.85 11.69 4 10 15 50 0.77 9.99 4 10 30 50 0.66 6.55 6 10 0 50 1.07 10.17 6 10 15 50 0.95 8.96 6 10 30 50 0.80 5.54 8 10 0 50 1.29 9.55 8 10 15 50 1.11 9.55 8 10 30 50 0.91 6.80 4 10 0 100 1.02 14.94 4 10 15 100 0.92 12.99 4 10 30 100 0.78 9.06 6 10 0 100 1.36 19.06 6 10 15 100 1.17 18.68 6 10 30 100 0.96 15.82 8 10 0 100 1.71 18.04 8 10 15 100 1.42 18.85 8 10 30 100 1.09 17.60 4 10 0 150 1.16 19.88 4 10 15 150 1.02 18.36 4 10 30 150 0.86 15.14 6 10 0 150 1.58 24.31 6 10 15 150 1.33 24.12 6 10 30 150 1.05 21.51 8 10 0 150 2.02 22.83 8 10 15 150 1.61 24.12 8 10 30 150 1.16 25.41 4 20 0 50 1.03 17.52 4 20 15 50 0.95 15.76 4 20 30 50 0.82 12.02 6 20 0 50 1.28 15.63 6 20 15 50 1.14 14.81 6 20 30 50 0.96 11.73 8 20 0 50 1.51 14.70 8 20 15 50 1.31 14.34 8 20 30 50 1.06 11.86 4 20 0 100 1.22 19.97 4 20 15 100 1.10 18.36 4 20 30 100 0.94 19.27 6 20 0 100 1.57 23.80 6 20 15 100 1.36 22.32 6 20 30 100 1.10 20.84 8 20 0 100 1.94 21.79 8 20 15 100 1.61 23.01 8 20 30 100 1.19 23.01 4 20 0 150 1.37 24.72 4 20 15 150 1.21 23.15 4 20 30 150 1.02 19.97 6 20 0 150 1.80 28.17 6 20 15 150 1.52 28.17 6 20 30 150 1.17 26.72 8 20 0 150 2.26 24.36 8 20 15 150 1.79 28.39 8 20 30 150 1.21 24.36 62. ACTA GeOTeCHNICA SLOVENICA, 2013/1 p. ÎelusiS & b. slender: soil-nail wall stability analysis using anfis shows that the calculations were made for 81 different combinations of design data. 4.2 ANFIS STRUCTURE For a Sugeno fuzzy model ANFIS-SF, a rule set with n fuzzy "If-then" is as follows: Rule 1: If l is In1MF1 and a is In2MF1 and ft is In3MF1 and d is In4MF1 then SF1 = k1 ■ l + k1 ■ a + ■ b + k\ ■ d + k1 (14) Rule 2: If l is In1MF2 and a is In2MF2 and ft is In3MF2 and d is In4MF2 then SF2 = k2 ■ l + k2 ■ a + k2 ■ b + K ■ d + K (15) 2 j , 1,2 where k0,kj,k^,k^,k\,kg,k2,k^,k2,k2 are the consequent parameters and l, a, ft and d are the input variables. The sign In1MF1 stands for the membership function 1 in input 1. The output of each rule is equal to a constant and the final output is the weighted average of each rule's output. SF = £ w, ■ SFt (16) i=1 n _ hopt wi ■ hopt,i (17) i=1 The weights are obtained from the Gaussian membership function. mA(x)=exp 42 < (18) where c is the position of the centre of the curve's peak and a is the width of the curve. The parameters c and a are the premise parameters. The first membership grade of a fuzzy set (In1MFi , In2MFi , In3MFi , In4MF) is calculated with the following equations: minlMF, (l) = exp mIn2MFi (a) = exp l - c In1 MFi V2 sIn1MFi In2MF: V2 sIn2MF (19) (20) mIn3MF: (b) = exp mIn4MF (d) = exp b-c In3MF V2 aIn3MR d 1 In4MF V2 sIn4MF (21) (22) where l, a, ft and d are the inputs to the Gaussian membership function. Next, the product of the membership function for every rule is calculated: wi = mIn1MF (l) ■ mIn2MFi (a) ■ mIn3MFi (b) ■ mIn4MFS (d) (23) where wi represents the fire strength of the rule i. The ratio of the ith rule's firing strength to the sum of all the rules' firing strengths is defined with: w , for i = 1, 2,..., n. (24) The input data represents a node on the left, and the right node represents the output data (Fig. 4). In order to achieve the desired input-output mapping, the consequent and premise parameters are updated according to the given training data (Table 4) and the hybrid learning procedure. This hybrid learning procedure [25] is composed of a forward pass and a backward pass. In the forward pass, the algorithm uses the least-squares method to identify the consequent parameters. In the backward pass the errors are propagated backwards and the premise parameters are updated with the gradient descent. The premise (a i, c,) and the consequent (k) parameters of the ANFIS-SF model and ANFIS-INC model are given in Table 5. By setting the premise and consequent parameters, the space of four variables is described. 4.3 TESTING THE ANFIS MODELS To test and validate the ANFIS models, a data set was selected that was not used during the training of the network (Table 6). The coefficient of determination (R-square) and the root-mean-squared error (RMSE) between the predicted and the calculated values is taken as a measure of the performance. The calculated and predicted safety factors SFs for the model ANFIS-SF are shown in Fig. 5. 62. ACTA GeOTeCHNICA SLOVENICA, 2013/1 P. ÎELUSIS & b. slender: soil-nail wall stability analysis using anfis Input Input membership Rule Output membership Output function function Figure 4. Structure of ANFIS-SF and ANFIS-INC model. Table 5. The premise and consequent parameter for the ANFIS-SF and ANFIS-INC model. ANFIS-SF In1MFi In1MF2 In2MF1 In2MF2 In3MF1 In3MF2 In4MF1 In4MF2 c 3.229 10.554 11.334 14.691 16.367 5.660 61.041 109.767 a 2.523 11.309 11.792 17.988 18.558 9.861 48.925 136.300 k1 k2 k3 k4 kc Rule 1 0.022 0.015 -0.005 0.001 0.485 Rule 2 -0.090 -0.006 -0.082 -0.002 5.079 ANFIS-INC In1MF1 In1MF2 In2MF1 In2MF2 In3MF1 In3MF2 In4MF1 In4MF2 c 2.488 7.787 12.256 5.282 17.411 21.318 43.831 99.530 a 2.572 5.583 11.995 4.191 18.402 25.186 54.367 50.416 k1 k2 k3 k4 k> Rule 1 -6.552 0.175 0.199 0.004 86.382 Rule 2 -3.885 0.449 -0.058 -0.355 27.041 Analytical snlutiDns SF (-) Figure 5. Calculated versus predicted safety factor SF. 62. ACTA GeOTeCHNICA SLOVENICA, 2013/1 P. ÎELUSIS & b. slender: soil-nail wall stability analysis using anfis Table 6. Testing data compared to ANFIS models. l a ß d SF ANFIS-SF Error OPTINC ANFIS-INC Error (m) (°) (°) (mm) (-) (-) (%) (°) (°) (%) 5 5 10 75 0.91 0.89 2.1 11.1 10.2 8.1 7 5 10 75 1.15 1.14 1.2 12.0 12.3 2.8 5 5 10 125 1.08 1.08 0.3 16.2 17.8 10.0 7 5 10 125 1.38 1.41 2.3 19.1 20.6 7.9 5 15 10 75 1.09 1.07 1.9 17.3 15.2 12.1 7 15 10 75 1.34 1.33 1.0 17.0 16.8 1.1 5 15 10 125 1.27 1.27 0.3 21.7 22.7 4.6 7 15 10 125 1.58 1.61 1.8 23.7 24.6 3.7 5 5 20 75 0.84 0.81 3.0 10.1 8.9 11.5 7 5 20 75 1.04 1.02 1.9 12.0 11.5 4.0 5 5 20 125 0.98 0.97 0.7 15.5 16.6 6.9 7 5 20 125 1.22 1.24 1.4 18.7 20.3 8.2 5 15 20 75 1.01 0.98 2.7 14.5 14.0 3.4 7 15 20 75 1.21 1.19 1.7 16.9 16.1 4.6 5 15 20 125 1.15 1.15 0.4 23.4 21.6 8.0 7 15 20 125 1.38 1.40 1.3 23.9 24.4 1.8 The comparison between the optimal inclination of the soil nail calculated with the optimization model OPTINC and the ANFIS-INC model is shown in Fig. 6. The higher coefficient of determination means that the correlation between the monitored and predicted data is high; it does not mean that the monitored data are close to the predicted data. To estimate the error of each model we used RMSE. 5 discussion The presented ANFIS models are able to predict the safety factor SF (-) and the optimal inclination of the soil nails q (°) for different inclinations of the wall a (°), the slope angle of the terrain ft (°), the length of the nails l (m) and the hole diameter d (mm). Fig. 7 shows the safety factor depending on the length of the soil nails and the inclination of the wall. With increasing inclination of 10 1Î 14 16 16 30 32 îi Optimal inclination i] (") otHained witfi modei OPTINC Figure 6. Calculated versus predicted optimal inclination of soil nails n. 62. ACTA GeOTeCHNICA SLOVENICA, 2013/1 p. ÎelusiS & b. slender: soil-nail wall stability analysis using anfis the wall, the length of the soil nail and the hole diameter, the safety factor is increasing. On the other hand, the safety factor is decreasing with the increasing slope angle of the terrain. The larger hole diameter increases the pull-out resistance; therefore, the safety factor is larger. Fig. 8 shows the optimal inclination of the soil nails depending on the length of the soil nails and the inclination of the wall. With increasing inclination of the wall, the optimal inclination of the soil nails increases. The lengths of the soil nails have a small influence on the optimal inclination of the soil nail in comparison with the inclination of the wall. 6 conclusions The paper presents the safety-factor optimization for a soil-nail wall. The optimization is performed using the non-linear programming (NLP) approach. For this purpose, the NLP optimization model OPTINC was developed. The model comprises the safety-factor objective function, which is subjected to geo-mechanical and design constraints. As the model was developed in a general form, the optimization of the system can be performed for different heights of the wall, spacings of the nails as well as for different soil environments. The output of the OPTINC model is an optimal inclination of the soil nails. The safety factor and the optimal inclination of the soil nails from the horizontal direction depends on the design of the soil-nail wall. Therefore, the safety factor SF (-) and the optimal inclination of the nails q (°) were calculated for different inclinations of the wall a (°), the slope angle of the terrain ft (°), the length of the nails l (m) and the hole diameter d (mm). Based on these results two ANFIS models were developed. The ANFIS- Inclination of wall a p) u 4 Lwigth of soil nail I