# I. Introduction p to the present, the limit equilibrium method (LEM) is still the most popular method as used by engineers and researchers for slope stability analysis. In general, LEM can be classified under two major groups: "simplified" methods and "rigorous" methods. Traditionally, the LEM is taken to be a statically indeterminate problem, and assumptions on the distributions of internal forces are required for the solution of the factor of safety (Cheng and Lau 2014). Various methods of analysis are adopted for various engineering applications, and the Spencer method appears to be the most popular at present. Cheng et al. (2010), Cheng et al. (2011) and Cheng et al. (2013) have pointed out that if the concept of extrema (extremum principle or equivalently numerical variational principle) or the ultimate state is considered, then there will be sufficient condition to solve a slope stability problem without the use of internal force distribution function or any other arbitrary assumption, and the LEM will become a statistically determinate problem. Furthermore, Cheng et al. (2010) have also found that the convergence problem using the Spencer method may affect the determination of the lowest factor of safety and the location of critical failure surface. Cheng et al. (2010) Cheng et al. (2011) and Cheng et al. (2013) have however also pointed out that for normal problems, the extrema will be close to the classical solutions by Spencer so that the determination of the extrema is necessary only for complicated problems. The power of the LEM is finally illustrated by the equivalence between the bearing capacity, lateral earth pressure and slope stability problems. The extremum from the LEM is equivalent to the results from plasticity solutions, and the results from the LEM can be a good approximation of the solution of a general geotechnical problem. The strength reduction method (SRM), implemented through the finite element method, was applied for slope stability analysis as early as 1975 by Zienkiewicz et al. Later the SRM was applied by Naylor (1982), Donald and Giam (1988), Matsui and San (1992), Ugai and Leshchinsky (1995), Dawson et al. (1999), Griffiths and Lane (1999), Zheng et al. (2005), Cheng et al. (2007a), Wei et al.(2009), Cheng (2009a, 2009b), Wei and Cheng (2010) and Nian et al. (2012). More recently, the SRM was implemented by other numerical procedures such as the mesh-free method (MFM) and the spectral-element method (SEM) (Tiwari, 2015). SRM technique has also been implemented into several commercial geotechnical finite element programs for engineering applications. A detailed discussion and study about the use of SRM in slope stability analysis has been given by Cheng et al. (2007a). Various problems including sensitivity to mesh design, size of solution domain, dilation angle, numerical instability with different SRM computer programs have been identified by Cheng et al. (2007a), and some program developers have updated their programs in accordance with the identified problems. Griffiths and Lane (1999 2015) compared LEM, SRM and limit analysis methods, and showed that for steep slopes with low factors of safety, the flow rule may have a significant influence on the comparisons, while numerical instabilities may occur in the case of nonassociated plasticity with large differences between the friction angle and dilation angle. Tshuchnigg et al (2015b) further investigated this phenomenon and proposed various approaches, based on the work on plasticity by Davis (1968), to overcome such obstacle. ) Shen and Karakus (2014) and Zhao et al. (2015) implemented the SRM with nonlinear failure criteria to study rock and soil slope stability, respectively, but they adopted different 'strength reduction strategies', and Zhao et al. (2015) concluded that the factors of safety obtained by SRM will be substantially influenced by these strategies, i.e., whether the 'cohesive' and 'stressdependent' components of shear strength are factored separately or simultaneously in the SRM analyses. Limit analysis does not require the interslice force function and is free of convergence problem which are unavoidable for the classical limit equilibrium method (except for the extremum principle by ). It has the advantages similar to the LEM in that no constitutive model and initial conditions are required, a flow rule is however required to specified which is usually not critical towards the factor of safety (similar to the interslice force function). For limit analysis, the upper bound approach is the more popular approach, and recently some commercial programs are available for the limit analysis of the stability geotechnical problems. The equivalency between limit analysis and LEM has been demonstrated by Leshchinsky et al. (1985). It is usually considered that the LEM methods cannot satisfy all of the equilibrium requirements. This understanding is true for the classical LEM, but has been demonstrated to be not true with the extremum principle by Cheng et al. (2010), Cheng et al. (2011) and Cheng et al. (2013). The uses of limit analysis for simple geotechnical stability problems have been discussed by Chen (1975), but such analytical approach is not practical for real problems with complex geometry and soil/geologic conditions. For limit analysis, a new approach called the discontinuity layout optimization method (DLO) has attracted the attention of some engineers and researchers. DLO procedure expresses the limit analysis problem entirely in terms of lines of discontinuity instead of elements as in the classical continuum problem (Smith and Gilbert 2007). Using DLO, a large number of potential discontinuities are set up at different orientations; while the continuum based element formulations, discontinuities are typically restricted to lie only at the edges of elements. With the use of modern optimization algorithms, an optimized solution can be achieved easily. After the initial success by Smith and Gilbert (2007), there are different works in DLO by Clarke et al. (2013), Smith and Gilbert (2013), Bauer and Lackner (2015), Al-Defae and Knappett (2015), Leshchinsky (2015), Vahedifard et al. (2014), Leshchinsky and Ambauen (2015). The original DLO formulation suffers from the limitation that only the translation mechanism can be considered. In view of such limitation, Gilbert et al. (2010) and later Smithy and have extended the DLO formulation to cover the rotational formulation. Since DLO is actually a numerical form of limit analysis, the basic limitation of limit analysis is similar to that for DLO. a) Some Case Studies with LEM and DLO Yu et al. (1998) have given a very detailed comparison between the use of limit analysis and LEM, and it is found that the results from the two methods are similar and comparable in most cases for relatively simple problems. Recently, DLO has been adopted for slope stability analysis by Leshchinsky and Ambauen (2015), and it is found that the results by DLO and LEM are comparable in general. Leshchinsky and Ambauen (2015) have however found some cases for which there are noticeable differences between the DLO and LEM, and they have concluded that DLO requires less assumption on the location of collapse, and therefore may be more preferable than LEM, especially for complex, yet realistic geotechnical problems. After reviewing the examples by Leshchinsky and Ambauen (2015), the authors tend to disagree with the results and comments by Leshchinsky and Ambauen (2015). There are some limitations in the works by Leshchinsky and Ambauen (2015) which include: 1) use of classical LEM method which are greatly affected by convergence problem (Cheng et al. 2008); 2) critical failure surface has not been determined (Fig. 12 from Leshchinsky and Ambauen 2015 has only considered 151 surfaces); 3) interslice force function can be critical in complex problems. As discussed by Cheng (2003) 2013) have overcome these problems and can provide solutions similar to some classical plasticity problems which are not possible with the classical LEM. A fair comparison and commentary on these methods must be based on reliable and robust analyses that identify the differences between DLO and LEM. Some problems with the DLO have been previously identified by Cheng (2018), and more studies will be carried out in this paper. With reference to Fig. 1 which is Fig. 5a by Leshchinsky and Ambauen (2015), the soil parameters are unit weight=19 kN/m3, c'=28 kPa and ?'=20°. The critical result by DLO pass below the toe of the slope at the right hand side of Fig. 1 by Leshchinsky and Ambauen (2015). On the other hand, the critical result by the authors using the heuristic optimization method and Spencer method developed by Cheng et al. (2007a) and Cheng and Lau (2014) pass through the toe of the slope. The critical result by Baker (1980) also pass through the toe of the slope while the critical result by Krahn and Fredlund (1997) (not shown for clarity) is similar to that by DLO but extends further to the right of the toe. The result by Krahn and Fredlund (1997) is not determined by the use of advanced optimization algorithm, and the adequacy of the result has not been confirmed. The authors have tried several updated commercial programs and have obtained results similar to that by the authors as shown in Fig. 1. Since the friction angle of the soil is 20° which is not a small value, the critical result by limit analysis will pass through the toe of slope as demonstrated by Chen (1975) using limit analysis. In views of the above discussion, the authors will suggest that the results by DLO cannot give the critical solution for such a simple case which is surprising to the authors. With reference to Fig. 2 where there is a 0.5m thickness of soft material for soil layer 2, the soil parameters are unit weight=19 kN/m 3 , c'=28 kPa and ?'=20° for layer 1 and unit weight=19 kN/m 3 , c'=0 kPa and ?'=10° for soil layer 2. The results by DLO and the (2015), Baker (1980), and Krahn and Fredlund (1997). When the authors increase the thickness of the soft layer to 1.5m, the critical result will still lie at the bottom of the soft layer with a factor of safety 1.14. Since the shear strength parameters at soil layer 2 are low, the weight of the soil tends to push the soft material to the right so that the critical slip surface should lie within the soft band, and the results by the authors are more reasonable as compared with other results. As discussed by Cheng (2007) and Cheng et al. (2012), the presence of a soft band is mathematically equivalent to a Dirac function, for which many optimization algorithms fail to work. The domain transformation technique by Cheng (2007) and the coupled optimization algorithm by Cheng et al. (2012) have effectively overcome this problem without any special precaution required by the engineers in the analysis. In Fig. 3 which is same as that for Fig. 2 with a pore pressure ratio 0.25 (Fig. 5d by Leshchinsky and Ambauen, 2015), the critical result by Leshchinsky and Ambauen (2015) lies at the top of the soft band while the critical results by the authors, Baker (1980), and Krahn and Fredlund (1997) (mistaken to be at the top of the soft band by Leshchinsky and Ambauen 2015) lie at the bottom of the soft band, and the inability to locate the critical result for a soft band by DLO is clearly illustrated. In Fig. 4 which is same as that for Fig. 2 with a prescribed water table (Fig. 5f In Fig. 5 (Fig. 12 by Leshchinsky and Ambauen, 2015), there are great differences between the critical result by the authors and Leshchinsky and Ambauen (2015). The soil parameters are unit weight=20 kN/m 3 , c'=0 kPa and ?'=30° for soil layer 1, unit weight=19 kN/m 3 , c'=0 kPa and ?'=45° for soil layer 2 and unit weight=19 kN/m 3 , c'=10 kPa and ?'=0° for soil layer 3. On the left hand side of the critical slip surface by Leshchinsky and Ambauen (2015), there is a very sudden change in the slope of the critical failure surface which seems unlikely to happen. At the right hand side of the critical slip surface by Leshchinsky and Ambauen (2015), the critical slip surface is nearly vertical, which is also highly unlikely, as the friction angle of soil layer 1 and 2 are 30° and 40° respectively with zero cohesive strength. When this same slip surface by Leshchinsky and Ambauen (2015) is considered with the M-P method using f(x)=sin(x), the authors actually get a factor of safety of 1.05, which is significantly greater than the result of 0.95 by Leshchinsky and Ambauen (2015). This problem is then reanalyzed by the authors using LEM to locate the critical slip surface. For this problem, the use of f(x)=1 is poor in convergence, and the authors get a slightly different critical slip surface and a factor of safety of 0.97 by using f(x)=1.0. As mentioned by Cheng et al. (2008Cheng et al. ( , 2010)), f(x) can be critical in some cases which will affect the optimized solution. In this respect, the authors have also adopted the extremum principle (Cheng et The authors have adopted an accuracy of 0.001 in all the global optimization search in the present study, and the global minima of each example has been tested with different optimization algorithms for confirmation. Based on the above case studies, it can be concluded that some of the past reported results in literature which are not optimized with the modern optimization algorithms may not be reliable enough for comparisons. In particular, for the presence of a soft band which is a difficult problem, the present study and the works by Cheng (2007), Cheng et al. (2012) have demonstrated that great care must be taken in order to obtain a good result. Furthermore, as a relatively new computational method, DLO has been demonstrated to be affected by the soft band or local minima problem. Overall, the authors view that the problems presented in this section are not fundamental deficiencies of DLO. Instead, they highlight the limitations of the numerical technique in implementing the DLO up to the present moment. With refined and improved numerical technique coupled with DLO, the authors expect that better results will be produced by DLO in the future. On the other hand, it is dangerous to compare the advantages and limitations of different stability methods based on old results or computer programs with limitations. Some of the comments in previous literature are possibly distorted by the limitations of the computational technique in computer programs instead of being the actual comparisons of different stability analysis methods. # b) Further Study on DLO, SRM and LEM Cheng et al. (2007a) and many others have conducted comparisons between LEM and SRM, and it is generally found that the factor of safety and the critical failure surface are not sensitive to the dilation angle, and the results from SRM are comparable to LEM in most cases. Cheng et al. (2007a) have however found many minor problems in several commercial SRM programs in the previous study, and many of these commercial programs have updated the programs with reference to the case studies by Cheng et al. (2007a). With reference to the 45°slope as shown in Fig. 6 which has been studied by Cheng et al. (2006), the authors have also found that results from DLO are comparable to LEM and SRM in many cases, but there are some cases where greater differences (more than 5% in Table 1) are observed between DLO and other methods, which are worth consideration. In Table 1, the critical factors of safety by LEM are obtained by the Spencer method, and the results are close to that by the extremum principle except for the three values for the case of zero friction angle which are marked with * in Table1. In the full comparisons between the three methods, it is found that the factors of safety from DLO are always greater than those by the other methods, and the differences become greater with smaller friction angle, but the differences between the critical failure surfaces from the three methods are however minor. For SRM, the authors have found some surprising results from another program (new version) for which the results are given by cases 10 to 13 in Table 1. The SRM2 analysis is very sensitive to the dilation angle when the friction angle approaches 45°, and the factors of safety (as shown in bracket in Table 1) from this program are particularly low for SRM2 analysis. In view of the surprising SRM2 results for ?'=45° for that particular SRM program, the dilation angle is varied and the results are shown in Table 2. It is noticed that the computer program is very sensitive to the dilation angle case, and a small change in the dilation angle will give a significant change in the factor of safety which is obviously not correct. Furthermore, it is also noticed that a smaller dilation angle sometimes result in a higher factor of safety, which is again obviously wrong. For the critical failure surface, there are great differences between the case for c'=2kPa, ?'=40° and c'=2kPa, ?'=35° as shown in Fig. 7. In fact, the result in Fig7b is similar to that by LEM, DLO and SRM1, and the factor of safety from it is also close to the other three methods. It appears that SRM program determine a wrong critical failure surface and factor of safety, but the reason behind such problem is unknown and surprising. Cheng et al. (2007a) have found many limitations in the commercial SRM programs, and it appears that the updated version of some SRM programs may still face numerical problems under some cases which should be addressed. Being a new numerical method for DLO, the authors have considered another interesting case for this method. For a slope with very low to zero cohesive strength, the critical failure surface will be a shallow face failure. If the friction angle is equal to the slope angle, then the critical factor of safety of the slope should be equal to 1.0. From Table 3, it is however found that if c' is 0.03 kPa to 0, the critical factor of safety is much greater 1.0 while the critical failure surface is not a near surface failure. As long as c' is not too small, the results from DLO will then be normal. The authors view that the surprising results from DLO as shown in Tables 1 and 3 are the problems of the numerical implementation instead of the problem of DLO itself. It is possible that these kinds of problems may be overcome in the future, and the reason for the numerical problems behind DLO must be investigated. It is also interesting to note that the authors have never found such problem for LEM and SRM programs so far. For the problem with a soft band at soil layer 2 as discussed by Cheng et al. (2007a), surprising results are again obtained by DLO. The unit weight of the soils are 19 kN/m 3 , and c'=20 kPa and ?'=35° for soil layer 1, c'=0 kPa and ?'=25° for soil layer 2 and c'=10 kPa and ?'=35° for soil layer 3. As discussed by Cheng et al. (2007a), it appears that some SRM programs are affected by the size of the solution domain. The factor of safety for LEM is obtained as 0.927 by the Spencer method by Cheng et al. (2007a), and this value lie within Year 2022 ( ) J the SRM1 and SRM2 results by Plaxis and the new version of Phase (8.0). On the other hand, the factor of safety appears to be highly dependent on the nodal number adopted in the analysis. Even if 2000 nodal number is adopted, the factor of safety from DLO still appears to be unsatisfactory which is given in Table 4. The results by DLO are higher than those by LEM or SRM under all cases in Table 4, and the differences are not minor. Surprisingly, the critical failure surface from DLO as shown in Fig. 9 is similar to that by LEM or SRM (Cheng et al. 2007a). From Table 4, it can be concluded that the most influential factor in a proper DLO analysis is the nodal number. x y 0,0 0,5 If the third layer of soil instead of the second layer of soil is a soft material, the factor of safety has been established to be 1.29 from Spencer method, 1.27 from Extremum principle and 1.33 for SRM2 for all the SRM programs as discussed by Cheng et al. (2007a), and f(x) is relatively important for the present case (in general f(x) is not negligible if the friction angle is low). On the other hand, the result by DLO will approach the above factor of safety when the nodal number is large enough (Table 5). However, while the critical failure surfaces from LEM and SRM agree quite well as shown in Fig. 10a and 10b, the critical failure surface from DLO extends further to the right in Fig. 10c. To further examine these results, the authors have found another local minimum 1.29 with the Spencer method for the failure surface as shown in Fig. 11, which is very similar to that one by DLO as shown in Fig. 10. The authors view that a local minimum has been obtained from the DLO analysis. It should also be notes that the failure surfaces in Fig. 10a and Fig. 11 bear virtually the same factors of safety, and the differences between the two values are small so that it can be viewed that there are two global minimum for this problem. LEM can analyzed such problem easily while it takes more effort for SRM to detect the result in Fig. 11. Actually, without the previous knowledge about the existence of another global minimum, engineers will miss the result in Fig. 11 easily. For DLO, the failure surface as given by Fig. 10a or 10b cannot be obtained even increasing nodal number as given in Table 5. In this respect, there are some inherent limitation in the present development of DLO. It should be noted that c' is zero and ?' is smaller for the soft band layer. In the parametric study, different shear strength are used and DLO, LEM, SRM 1 and SRM 2 are carried out. SRM 1 is a non-associated flow rule analysis with a dilation angle = 0, while SRM 2 is an associated flow rule analysis with a dilation angle = friction angle. The thickness of soft band is set to be 500mm, 5mm and 2mm with three different slopes, while the soil properties are kept to be the same for three cases. In this analysis, number of nodes are set to be 500, 1000, 2000 in DLO and number of mesh are set to be 2000, 5000, 10000 in SRM. For the LEM, number of slices is set to be 50. Additionally, the number of mesh in SRM is increased (maximum 10000) to resolve the cases which the difference between SRM 1 and SRM 2 is significantly large (> 8%). From the analysis, the FOS by SRM have great differences from those by DLO, particularly for the soft band problem with a thickness of 2mm. In model 1, the FOS are determined to be 1.163 and 1.396 from DLO (c' = 5 and 10 kPa for Soil A respectively) with 1000 nodes, while the FOS are 1.23 (c' = 5 kPa for Soil A) and 1.02 (c' = 10 kPa for Soil A) in SRM 1 with 5000 and 10000 meshes respectively. For SRM 2, the FOS are found to be 1.28 and 0.74 in the same cases, and FOS are 1.129 and 1.3362 from LEM. For the slope analysis by DLO, one thing as found is that the number of nodes is a factor affecting the result. There can be no solutions for some cases in the three models. The number of nodes is required to adjust to fix the problem. For example in model 2 with soft band of 2mm thickness, the result cannot be determined if the number of nodes is set to be 1000. If it is 900, the solution can be determined. Another thing found in DLO is that the FOS changes slightly with the large change in the number of nodes (from 500 increases to 2000) in all cases. In the SRM, the FOS obtained from SRM 2 should be larger than those from SRM 1 theoretically. Interestingly, the result from a SRM program shows that the FOS obtained from SRM 1 is larger than that from SRM 2 when the soft band becomes thinner or steeper. Apart from that, the difference between FOS from SRM 1 and SRM 2 is supposed to be small (around 2%), however, the differences are more than 2% for all the cases. In model 2 (c' = 10kPa for Soil A and the soft band has a thickness of 2mm), the FOS is found to be 1.05 from SRM 1 and 0.42 from SRM 2 with 10000 elements. The difference of FOS between them is 60.0%. If we look careful into the results in Tables in 6 to 8, we may be disappointed to find that commercial programs cannot perform well under some cases, and the errors can be significant or even ridiculous! For a new problem without any known solution, how should an engineer assess and accept the results from the computer programs is a difficulty issue. The reasons for these numerical problems should also be rectified in the updated versions of the commercial programs. In fact, the authors have tested several versions of the commercial programs and find that each version may give rise to different problems. So far, the authors cannot find a commercial SRM program which can be correct/reasonable under all cases! Actually, the authors ![and Cheng et al. (2007a) have demonstrated that there are no major differences between the factors of safety from SRM and LEM for normal problems. Similar conclusions have been reached by Lu et al. (2014) and Liu et al. (2015) for 2D and 3D slope geometries. On the other hand,](image-2.png "") ![, Cheng et al. (2010) and Cheng et al. (2013), these three problems can lead to relatively poor solution in some cases by the classical LEM, and the extremum principle by Cheng et al. (2010), Cheng et al. (2011) and Cheng et al. (](image-3.png "") 1![Fig. 1: Comparisons of DLO and LEM for Fig.5a by Leshchinsky and Ambauen (2015), FOS=2.03 by DLO and 1.94 by Spencer method](image-4.png "Fig. 1 :") 21434![Fig. 2: Comparisons of DLO and LEM for Fig.5b by Leshchinsky and Ambauen (2015), FOS=1.28 by DLO and 1.21 by Spencer method. The thickness of the soft band has been increased to 1.25m for the result in Fig.2b with FOS=1.14](image-5.png "Fig. 2 : 14 Fig. 3 :Fig. 4 :") 5![Fig.5: Comparisons of DLO and LEM for Fig.12by Leshchinsky and Ambauen (2015), FOS=0.95 by DLO and 1.0 by Spencer method by Leshchinsky and Ambauen (2015) using 151 trial surface, 0.92 by M-P method using f(x)=sin(x) and 0.97 by f(x)=1.0 by the authors using about 20000 trials with simulated annealing optimization method](image-6.png "Fig. 5 :") 6![Fig. 6: Comparisons between LEM, SRM and DLO (based on Cheng et al. 2007a)](image-7.png "Fig. 6 :") 7![Fig. 7: Critical failure surface for SRM2 analysis (c'=2 kPa), ?'=40° for (a) [cursor in (a)?] and 35° for (b)](image-8.png "Fig. 7 :") 8![Fig.8: A slope problem with a soft band as discussed byCheng et al. (2007a) ](image-9.png "Fig. 8 :") 9![Fig. 9: Critical failure surface from DLO for the problem in Fig. 8](image-10.png "Fig. 9 :") 2422Leshchinsky-Ambauen (2015)2018Spencer Method16141210A B8Bedrock642005101520253035404550Leshchinsky-Ambauen (2015)Spencer MethodLeshchinsky-Ambauen (2015)Spencer Methodand have obtained a critical solution 0.915 which is Year 2022 © 2022 Global Journals ( ) J close to that by using f(x)=sin(x). 1Casec'(kPa)?'FOS(LEM)FOS(DLO)FOS(SRM1)FOS(SRM2)1550.410.440.420.4221050.650.710.670.67310150.981.040.981.0142051.061.221.131.14520151.481.611.511.53620251.851.961.871.887500.2 (0.21*)0.240.210.2181000.4 (0.42*)0.470.440.4492000.8 (0.85*)0.950.890.89102451.351.401.421.44 (fail)115451.651.701.681.74 (0.98)1210452.042.092.052.15 (1.1)1320452.692.792.672.83 (1.59)between DLO, LEM and SRM for Fig.6(SRM1 means zero dilation angle, SRM2 means dilation angle=friction angle) 2Casec' (kPa)?' ( °)?' ( °)fos by SRM2fos by LEM124545no solution2 32 245 4544.9 40no solution 1.221.354245351.45545450.986 75 545 4544.9 401.19 1.471.658545351.6991045451.110 1110 1045 4544.9 401.49 2.052.04121045352.06132045451.5914 1520 2045 4544.9 402.03 2.622.69162045352.67 3Casec' (kPa)?' (°)FOS10.03301.4220.1301.0130.03351.6240.1351.0250.03401.8460.1401.02 4Year 2022Fig. 828m domain size, solution tolerance 0.01, different nodal densityCaseNodal No.FOS by DLOFOS by LEMFOS difference with LEM (DLO %)12501.3560.927-46.2825001.0690.927-15.32310001.0820.927-16.72420001.0550.927-13.8128m domain size, nodal density 500, different solution toleranceCaseSolution toleranceFOS by DLOFOS by LEMFOS difference with LEM (DLO %)10.011.0690.927-15.3220.0011.0690.927-15.3230.0041.0690.927-15.3240.0051.0690.927-15.32Solution tolerance 0.01, nodal density 500, different domain sizeCaseDomain Size (m)FOS by DLOFOS by LEMFOS difference with LEM (DLO %)1281.0690.927-15.322201.0930.927-17.913121.0250.927-10.57© 2022 Global Journals 5Year 2022( ) JDLO analysis of soft soil layer 3 withdifferent nodal density (28 domain, solution tolerance 0.01)Case Nodal No.FOS by DLOFOS by SRM2FOS by LEMFOS difference with LEM (DLO%)FOS difference with LEM (SRM2%)12501.4051.331.27-8.91-3.1025001.3581.331.27-5.27-3.10310001.351.331.27-4.65-3.10 6Year 2022Factor of safety FOS by DLO, LEM and SRM for case 1Number of Nodes (DLO)Number of Mesh (SRM)Thickness of soft band (mm)Soil Type c' (kPa)?' (°)FOS (DLO)FOS (LEM)FOS (SRM 1)FOS (SRM 2)500Soil A Soil B5 035 251.0961.10331.021.09500Soil A Soil B10 035 251.2751.18421.171.2450020005 5Soil A Soil B Soil A Soil B5 0 10 035 25 25 351.165 1.3961.1279 1.33471.2 0.981.28 0.9( ) J2 2Soil A Soil B Soil A Soil B5 0 10 035 25 35 251.166 1.3971.129 1.33621.24 0.961.28 0.75Soil A Soil B5 035 251.1621.12791.221.28100050005 2Soil A Soil B Soil A Soil B10 0 5 035 25 35 251.395 1.1631.3347 1.1291.22 1.230.95 1.282Soil A Soil B10 035 251.3961.33620.930.852000100005 2Soil A Soil B Soil A Soil B10 0 10 035 25 35 25No Solution 1.395 (1000 nodes) 1.396 (1000 nodes) No Solution1.3347 1.33621.21 1.020.96 0.74© 2022 Global Journals 7 8 of conditions.Number of Nodes (LEM)Number of Mesh (SRM)Thickness of soft band (mm)Soil Type c' (kPa)?' (°)FOS (DLO)FOS (LEM)FOS (SRM 1)FOS (SRM 2)500Soil A Soil B5 035 251.061.05220.981.04500Soil A Soil B10 035 251.2231.19371.111.1750020005 5Soil A Soil B Soil A Soil B5 0 10 035 25 25 351.128 1.3381.0924 1.28291.2 1.161.27 0.942Soil A Soil B5 035 25No Solution 1.127 (600 nodes)1.09351.211.272Soil A Soil B10 035 251.341.28331.050.45Soil A Soil B5 035 251.1231.09241.221.28100050002 5Soil A Soil B Soil A Soil B10 0 5 035 35 25 251.126 (900 nodes) No Solution 1.3371.0935 1.28291.2 1.371.27 0.792Soil A Soil B10 035 25No Solution 1.339 (900 nodes)1.28331.480.352000100005 2Soil A Soil B Soil A Soil B10 0 10 035 25 35 25No Solution 1.337 (1000 nodes) 1.339 (900 nodes) No Solution1.2829 1.28331.37 1.050.74 0.42Factor of safety FOS by DLO, LEM and SRM for case 3Number of Nodes (LEM)Number of Mesh (SRM)Thickness of soft band (mm)Soil Type c' (kPa)?' (°)FOS (DLO)FOS (LEM)FOS (SRM 1)FOS (SRM 2)500Soil A Soil B5 035 251.0211.02880.931500Soil A Soil B10 035 251.1691.1381.061.1350020005 5Soil A Soil B Soil A Soil B5 0 10 035 25 35 25No Solution 1.083 (600 nodes) No Solution 1.276 (600 nodes)1.0528 1.2171.1 1.060.99 0.742Soil A Soil B5 035 251.0841.05481.161.272Soil A Soil B10 035 25No Solution 1.278 (600 nodes)1.21061.420.295Soil A Soil B5 035 25No Solution 1.082 (700 nodes)1.05281.10.98100050005 2Soil A Soil B Soil A Soil B10 0 5 035 25 25 35No Solution 1.276 (950 nodes) 1.0821.217 1.05481.1 0.9950.93 1.262Soil A Soil B10 035 251.2771.21061.440.85Soil A Soil B5 035 25No Solution 1.082 (700 nodes)1.05281.080.772000100005 2Soil A Soil B Soil A Soil B10 0 5 035 25 35 25No Solution 1.276 (950 nodes) 1.081 (1900 nodes) No Solution1.217 1.05480.995 1.170.74 1.262Soil A Soil B10 035 25No Solution 1.276 (1800 nodes)1.21061.51.1 © 2022 Global Journals ( ) J © 2022 Global Journals authors are very similar except that the critical result by the authors lies at the bottom of the soft layer while the results by Leshchinsky and Ambauen (2015) lie at the top of the soft layer. The critical results by Baker (1980) and Krahn and Fredlund (1997) are also at the bottom of the soft band but are mistaken to be at the top of the soft band by Leshchinsky and Ambauen (2015). The authors reduce the thickness of the soft layer to 1mm, and the factor of safety as well as the critical result will then be equal to that by Leshchinsky and Ambauen ## Acknowledgement The present project is funded by the National Natural Science Foundation of China (Grant No. 51778313) and Cooperative Innovation Center of Engineering Construction and Safety in Shangdong Blue Economic Zone. have constructed more problem cases, but the examples as shown here are sufficient to illustrate the difficulty to develop a numerical algorithm/computer program which can pass through all types of material parameters and geometry. ## II. Discussion In this study, the three major methods for slope stability analysis have been studied and compared. In general, all the three methods will give similar results under normal cases. Through this study, several problems are however identified which must be considered with care. It is found that DLO always gives a higher factor of safety as compared with LEM or SRM. Even though the differences are small in most cases, there are also cases where the differences are appreciable (> 10%), for which the results by DLO have to be used with some caution. DLO has been demonstrated to be sensitive to the nodal number, and an adequate nodal number should always be adopted in practice. The authors have also found some surprising results from SRM analysis. The authors have not tested all the SRM programs in market, and there are only limited case studies (including those not shown in this paper) to conclude the performance of the SRM programs. The sensitivity of a SRM analysis when the friction angle and dilation angle are high should be noted, and the results should be carefully assessed. Cheng et al. (2007a) have mentioned the difficulty in the nonlinear solution scheme and the assessment of the critical condition from SRM, and it appears that some current SRM programs may still fail to work properly under some special cases. The authors have also noted that other programs occasionally give SRM1 factor of safety slightly higher than that for SRM2, but the differences are usually small and not critical. There are some cases where the results from SRM are highly unacceptable, but without the knowledge on the acceptable results, how can an engineer judge and accept the results from any SRM program? An interesting issue to note is the occurrence of multiple local minimum in a slope analysis. For the results in Fig. 10 and Fig. 11, there are actually two critical failure surfaces with the same factor of safety 1.29 using the Spencer method which can be considered as multiple global minimum problem. Using the Extremum principle, the authors have however obtained only one global minimum for this problem, but it is possible (though rare) that there are multiple global minimum (with the same minimum value) even with the extremum principle. As mentioned by Cheng et al. (2007a), the use of LEM for such case is simple. The authors simply choose to view all the trial failure surfaces with a factor of safety within 2% (or other value) from the global minimum, and the problem of local minimum or even multiple global minimum can be detected immediately. On the other hand, such application in SRM is still not automatic in general, and usually the users need to exercise some kind of tricks in order to detect the other local minimum or global minimum. It is very easy for the users to miss the other global or local minimum from SRM, and very few users carry out this check in routine design work as SRM is much more time consuming as compared with LEM. For DLO, it appears to be trapped by the local minimum in the present study. The authors are also not aware of any simple method to detect all local/global minima by DLO, at least up to the present development of the method. In this respect, the authors view that DLO is still green at present, and there are still plenty of works ahead for enhancing DLO. ## III. Conclusion The purpose of the present study is not an assessment of the slope stability programs, but an assessment of the slope stability methods. In general, the authors view that DLO, LEM and SRM can be effective for normal cases. There are, however, many problems identified in SRM and DLO which engineers and researchers should consider. LEM has been developed for many years, and the problem of convergence, location of critical failure surface and local minimum have been studied in depth by Cheng (2003), Cheng (2007) (2013) and many others. LEM can be considered to be a very mature and robust tool to the engineers for various difficult problems. On the other hand, there are still some minor problems for SRM and some major problems for DLO. The authors tend to view that the problems as found are the results of the numerical implementation instead of the nature of SRM and DLO. With the continuous development in SRM and DLO, the authors believe that these two methods can become robust tools to the engineers in the future. The authors would also like to point out that some previous research works compare the results by LEM with other methods, and such comparisons may be misleading since many previous LEM results are not accurate enough due to the convergence problem or issues with finding the global minimum. This problem is clearly illustrated for the cases in Figs. 1 to 4 in this study as well as some other research works. The previous LEM results by other researchers (without the use of modern optimization method) are not the critical results for these cases, and the comparisons based on such results can be misleading. Although the authors view that DLO, LEM and SRM can all perform well in general, the authors tend to prefer LEM at present for normal routine analysis and design. LEM is simple to operate and robust for a variety * Newmark sliding block model for pile-reinforced slopes under earthquake loading AHAl-Defae JAKnappett Soil Dynamics and Earthquake Engineering 75 2015 * Gradient-based adaptive discontinuity layout optimization for the prediction of strength properties in matrix-inclusion materials SBauer RLackner International Journal of Solids and Structures 63 2015 * Nonlinear finite element analysis of safety factors RB JBrinkgreve HLBakker Proceedings of the Seventh International Conference on Computer Methods and Advances in Geomechanics the Seventh International Conference on Computer Methods and Advances in GeomechanicsCairns, Australia 1991. 6-10 May 1991 * WFChen Limit Analysis and Soil Plasticity USA Elsevier 1975 * Locations of Critical Failure Surface and some Further Studies on Slope Stability Analysis YMCheng Computers and Geotechnics 30 2003 * Global optimization analysis of slope stability by simulated annealing with dynamic bounds and Dirac function, Engineering Optimization YMCheng 2007 39 * Twodimensional Slope Stability Analysis by Limit Equilibrium and Strength Reduction Methods Y MCheng TLansivaara WWei Computers and Geotechnics 34 2007a. 2007 * Studies on six heuristic global optimization methods in the location of critical slip surface for soil slopes YMCheng LLi SCChi Computers and Geotechnics 34 2007b * Impact of Convergence on Slope Stability Analysis and Design YMCheng TLansivaara JSiu Computers and Geotechnics 35 1 2008 * Evaluation of interslice force function and discussion on convergence in slope stability analysis by the lower bound method YMCheng ZHZhao YJSun Journal of Geotechnical and Geoenvironmental Engineering 136 8 2010 ASCE * Limit equilibrium method based on approximate lower bound method with variable factor of safety that can consider residual strength YMCheng DZLi LLi YJSun RBaker YYang Computers and Geotechnics 38 2011 * A coupled particle swarm and harmony search optimization algorithm for difficult geotechnical problems, Structural and Multidisciplinary Optimization YMCheng LLi YJSun SKAu 2012 45 * The use of internal and external variables and extremum principle in limit equilibrium formulations with application to bearing capacity and slope stability problems, Soils and Foundation YMCheng TLansivaara RBaker NLi 2013 53 * YMCheng CKLau Soil Slope stability analysis and stabilization -new methods and insights Spon Press 2014 2nd edition * Equivalence between Bearing Capacity, Lateral Earth Pressure and Slope Stability Problems by Slip-line and Extremum Limit Equilibrium Methods YMCheng NLi International Journal of Geomechanics 17 12 4017113 2017 ASCE * YMCheng Slope stability and reliability analysis Nova Publishing 2018 * Modelling discrete soil reinforcement in numerical limit analysis SDClarke CCSmith MGilbert Canadian Geotechnical Journal 50 2013 * Theories of plasticity and failure of soil masses EHDavis Soil mechanics: selected topics 1968 * Slope stability 2000: proceedings of sessions of Geo-Denver EDawson FMotamed SNesarajah WRoth 2000. 2000. August 5-8, 2000 Denver, Colorado Geotechnical stability analysis by strength reduction * Slope stability analysis by strength reduction EMDawson WHRoth ADrescher Geotechnique 49 6 1999 * Application of the nodal displacement method to slope stability analysis IBDonald SKGiam Proceedings of the 5th Australia-New Zealand conference on geomechanics the 5th Australia-New Zealand conference on geomechanicsSydney, Australia 1988 * Masonry arch analysis using discontinuity layout optimisation MGilbert CSmith TPritchard Proc. Instn Civ. Engrs -Engng Comput. Mech 163 3 2010 * Slope stability analysis by finite elements DVGriffiths PALane Geotechnique 49 3 1999 * JKrahn The 2001 R.M. Hardy Lecture; the limits of limit equilibrium analysis 2003 40 * Evaluation of the University of Saskatchewan slope stability program JKrahn DGFredlund R.T.A. 1997 Quebec City, Canada C Forum, Roads and Transportation Association of Canada * Bearing Capacity of Footings Placed Adjacent to c'-?' Slopes BLeshchinsky Journal of Geotechnical and Geoenvironmental Engineering 6 4015022 2015 ASCE * Limit Equilibrium and Limit Analysis: Comparison of Benchmark Slope Stability Problems BLeshchinsky SAmbauen Journal of Geotechnical and Geoenvironmental Engineering 10 4015043 2015 ASCE * Slope stability analysis using the limit equilibrium method and two finite element methods SYLiu LTShao HJLi Computers and Geotechnics 63 2015 * Comparison between 3D limit equilibrium and shear strength reduction methodologies HHLu LMXu MDFredlund DGFredlund ASCE GSP 234 2014 * Finite element slope stability analysis by shear strength reduction technique TMatsui KCSan Soils and Foundations 32 1 1992 * The Evaluation of Slope Stability -A 25 Year Perspective, Stability and Performance of Slopes and Embankments-II NRMorgenstern Geotechnical Special Publication 31 1992 ASCE * Finite elements and slope stability DJNaylor Proc. Of the NATO Advanced Study Institute Of the NATO Advanced Study InstituteLisbon, Portugal 1982. 1981 * Three-dimensional strength-reduction finite element analysis of slopes: geometric effects TKNian RQHuang SSWan GQChen Canadian Geotechnical Journal 49 2012 * Slope stability evaluations using the finite element method. NGM MGSaeterbo Glamen SNordal XIV Nordic Geotechnical Meeting 1 2004. 2004 * Three-dimensional numerical analysis for rock slope stability using shear strength reduction method JShen MKarakus Canadian Geotechnical Journal 51 2014 * Identification of rotational failure mechanisms in cohesive media using discontinuity layout optimization CCSmith MGilbert Geotechnique 63 14 2013 * Application of discontinuity layout optimization to plane plasticity problems CSmith MGilbert Proc. R. Soc. A 463 2007. 2086 * Simplified Numerical Implementation in Slope Stability Modeling RCTiwari International Journal of Geomechanics 15 3 4014051 2015 * Comparison of finite-element limit analysis and strength reduction techniques FTschuchnigg HFSchweiger SWSloan AVLyamin IRaissakis Geotechnique 65 4 2015a * Slope stability analysis by means of finite element limit analysis and finite element strength reduction techniques. Part I: Numerical studies considering non-associated plasticity FTschuchnigg HFSchweiger SWSloan Computers and Geotechnics 70 2015b * Threedimensional limit equilibrium and finite element analysis: a comparison of results KUgai DLeshchinsky Soils and Foundations 35 4 1995 * Impact of Cohesion on Seismic Design of Geosynthetic-Reinforced Earth Structures FVahedifard BLeshchinsky SSehat DLeshchinsky Journal of Geotechnical and Geoenvironmental Engineering 140 6 4014016 2014 ASCE * Threedimensional slope failure by strength reduction and limit equilibrium methods WBWei YMCheng LLi Computers and Geotechnics 36 2009 * Soil nailed slope by strength reduction and limit equilibrium methods WBWei YMCheng Computers and Geotechnics 37 2009a * Strength reduction analysis for slope reinforced with one row of piles WBWei YMCheng Computers and Geotechnics 36 2009b * Stability analysis of slope with water flow by strength reduction method, Soil and Foundations WBWei YMCheng 2010 50 * Limit analysis versus limit equilibrium for slope stability HSYu RSalgado SWSloan KimJ M Journal of Geotechnical and Geoenvironmental Engineering 125 10 1998 ASCE * Effects of shear strength reduction strategies on safety factor of homogeneous slope based on a general nonlinear failure criterion LZhao FYang YZhang HDan WLiu Computers and Geotechnics 63 2015 * Geotechnical engineering limit analysis using finite element method. Rock and Soil Mechanics YRZheng SYZhao WXKong CJDeng 2005 26 * Associated and non-associated viscoplasticity and plasticity in soil mechanics OCZienkiewicz CHumpheson RWLewis Geotechnique 25 4 1975