Estimating the Slope Safety Factor Using Simple Kinematically Admissible Solutions

Determining soil and water conditions is essential for designing the optimal foundation and safely transferring loads, including the self-weight of structures, to the ground. Excessive or uneven settlement of the subsoil may ultimately lead to the formation of structural cracks in buildings or the loss of slope stability. In extreme cases, the damage results in structural failure. This paper presents the application of simple solutions from plasticity theory—an evaluation of the upper and lower bounds of the exact solution—to estimate the slope safety factor. It is demonstrated that simple kinematically admissible mechanisms for the non-associated flow rule provide solutions are close to those obtained from the traditional Fellenius method.


Introduction
Slope stability is one of the issues related to geotechnical design.It is particularly relevant when the building is founded on a site with a significant slope or one is in its immediate vicinity.In such cases, there is a danger that the slope may slip together with the building or trigger a landslide onto the building, which may lead to its damage or even structural failure.Slope failures can also be caused by the disruption of the slope's natural stability (e.g., earthquake-induced landslides), which can result in property damage, injury, loss of human life, as well as contamination and other serious environmental problems [1][2][3].In the context of preventing these accidents, it is of great importance to assess the safety factor.
The current standard, Eurocode 7 (EN 1997-1) [4], does not include a single selected computational approach according to which the stability of slopes should be checked.However, it is recommended to check for the total moment and vertical stability of the wedge of detachment.Instead, it is stated what to pay attention to when adopting a particular calculation method and what equilibrium equations it verifies.In addition, in the case of not verifying the equilibrium/horizontal stability, it is recommended to assess the horizontal direction of forces among different slices of soil.The characteristic features in terms of verified equilibrium states of the models summarized in a paper [5] or analyzed in [6][7][8][9] are helpful in this case.The authors of [5] recommend using Spencer's method and Bishop's simplified method for slope analysis.Analyses using these two methods to determine slope stability can be found in many studies, including [10][11][12][13].
Numerous other computational methods can be found in the literature and engineering practice.In addition to the traditional approaches discussed in [5], such as the Fellenius method, Janbu, or Sarma, more advanced computational approaches can be considered, for example, those based on probabilistic methods [14][15][16][17][18][19]; semi-probabilistic methods [20], including the stochastic finite difference method [21]; and those based on the finite element method [10,11,22].
At the same time, it should be noted that, currently, with the development of digital methods, there are more and more computer programs and digital computational models that allow the analysis of even very complex cases, as evidenced by, e.g., [1][2][3].
Despite this development, relatively simple methods are still used and recommended for less complex issues [23].
In the present paper, we discuss applying simple solutions from plasticity theorywith an evaluation of the upper and lower bounds of the exact solution-to estimate the slope safety factor.We demonstrate that simple kinematically admissible solutions for the non-associated flow rule give solutions close to those obtained using the Fellenius method [4].
Critical loads and failure mechanisms in the case of passive earth pressure (earthmoving machinery tools) were analyzed using plasticity theory.Due to the adoption of a rigid-perfectly-plastic model, it is difficult or even impossible to obtain exact solutions that satisfy all the conditions of statics and kinematics.Hence, to estimate admissible loads, statically or kinematically admissible solutions are built as the lower or upper bounds of the exact solution.Since the magnitudes for the exact solution are between the lower and upper bounds, kinematically admissible solutions give the lowest force estimate (upper bound), and statically admissible solutions give the highest force estimate (lower bound).
Plasticity theory was applied to solve boundary issues for the plane deformation state [24], assuming associated or non-associated flow rules.Experimental results for the pressure problem of walls with different shapes (smooth wall, loader bucket models) in a loose medium (sand) show good agreement with calculations for the non-associated flow rule [25].
Ongoing studies have shown that for a cohesive medium, the kinematically admissible mechanisms [26] describe the kinematics of the failure process quite well for both active and passive pressures.Even complex processes related to the working problem of loader bucket models can be described in the same manner [25,[27][28][29].Associated and non-associated flow rules have been applied to describe such processes.For example, the non-associated flow rule has been used to describe a quite advanced wall pushing process [27,29].
In the present study, plasticity theory and simple kinematically admissible mechanisms were used to evaluate slope stability.The obtained results, easy to calculate, were compared with those obtained using the Fellenius method.The kinematically admissible mechanisms were applied within the landslide area, as determined by the method mentioned above.

Materials and Methods
The issue of slope stability can be analyzed using plasticity theory and a stress function G σ ij called the plastic potential, where the relationship between strain and stress rates is defined by .
where λ is a certain scalar multiplier that has a constant value at a given state and point; it has different values at different points.
If the plastic potential is identified with the boundary condition (e.g., Mohr-Coulomb), Equation ( 1) is called an associated flow rule.If it is a different condition (Tresca condition was adopted in this paper), we deal with a non-associated flow rule.
In the general spatial issue of quasi-static motion, we have 16 unknowns to be determined, and we have 16 equations.
This system is quite challenging to solve, but effective solutions can be obtained for special states, such as the planar flow state or axisymmetric flows.
In the case of planar surface flow, we have: This system is quite challenging to solve, but effective solutions can be obtained for special states, such as the planar flow state or axisymmetric flows.
In the case of planar surface flow, we have: where γ is a volumetric weight; This system is quite challenging to solve, but effective solutions can be obtained for special states, such as the planar flow state or axisymmetric flows.
In the case of planar surface flow, we have: where γ is a volumetric weight;  a limit state condition, e.g., Mohr-Coulomb;  The relationships between the components of the stress vector and tensor, which take the form: A practical method of solving plane flow problems involves first determining the state of stress and then the velocity field.Figure 1 shows an example of its application to the analysis of smooth wall pressure (passive pressure).Despite the simplifications, these solutions are labor-intensive and have limited practical application.In addition, adopting a rigid-perfectly-plastic model makes it difficult to obtain exact solutions, especially specific solutions that satisfy all the conditions of statics and kinematics.Hence, statically or kinematically admissible solutions are built as approximations, which, for the associated flow rule, represent the lower bound (statically admissible solutions) and upper bound (kinematically admissible solutions) to the exact solution.The solution is located between these quantities.
The upper bound means that, for example, the allowable loads are no greater than those specified, and the lower bound is not less than that specified.The smaller the difference between these bounds, the closer to the exact solution we are.This system is quite challenging to solve, but effective solutions can be obtained for special states, such as the planar flow state or axisymmetric flows.
In the case of planar surface flow, we have: where γ is a volumetric weight;  a limit state condition, e.g., Mohr-Coulomb;  The relationships between the components of the stress vector and tensor, which take the form: A practical method of solving plane flow problems involves first determining the state of stress and then the velocity field.Figure 1 shows an example of its application to the analysis of smooth wall pressure (passive pressure).Despite the simplifications, these solutions are labor-intensive and have limited practical application.In addition, adopting a rigid-perfectly-plastic model makes it difficult to obtain exact solutions, especially specific solutions that satisfy all the conditions of statics and kinematics.Hence, statically or kinematically admissible solutions are built as approximations, which, for the associated flow rule, represent the lower bound (statically admissible solutions) and upper bound (kinematically admissible solutions) to the exact solution.The solution is located between these quantities.
The upper bound means that, for example, the allowable loads are no greater than those specified, and the lower bound is not less than that specified.The smaller the difference between these bounds, the closer to the exact solution we are.
The relationships between the components of the stress vector and tensor, which take the form:  A practical method of solving plane flow problems involves first determining the state of stress and then the velocity field.Figure 1 shows an example of its application to the analysis of smooth wall pressure (passive pressure).
This system is quite challenging to solve, but effective solutions can be obtained for special states, such as the planar flow state or axisymmetric flows.
In the case of planar surface flow, we have: where γ is a volumetric weight;  a limit state condition, e.g., Mohr-Coulomb;  The relationships between the components of the stress vector and tensor, which take the form: A practical method of solving plane flow problems involves first determining the state of stress and then the velocity field.Figure 1 shows an example of its application to the analysis of smooth wall pressure (passive pressure).Despite the simplifications, these solutions are labor-intensive and have limited practical application.In addition, adopting a rigid-perfectly-plastic model makes it difficult to obtain exact solutions, especially specific solutions that satisfy all the conditions of statics and kinematics.Hence, statically or kinematically admissible solutions are built as approximations, which, for the associated flow rule, represent the lower bound (statically admissible solutions) and upper bound (kinematically admissible solutions) to the exact solution.The solution is located between these quantities.
The upper bound means that, for example, the allowable loads are no greater than those specified, and the lower bound is not less than that specified.The smaller the difference between these bounds, the closer to the exact solution we are.Despite the simplifications, these solutions are labor-intensive and have limited practical application.In addition, adopting a rigid-perfectly-plastic model makes it difficult to obtain exact solutions, especially specific solutions that satisfy all the conditions of statics and kinematics.Hence, statically or kinematically admissible solutions are built as approximations, which, for the associated flow rule, represent the lower bound (statically admissible solutions) and upper bound (kinematically admissible solutions) to the exact solution.The solution is located between these quantities.
The upper bound means that, for example, the allowable loads are no greater than those specified, and the lower bound is not less than that specified.The smaller the difference between these bounds, the closer to the exact solution we are.
In the case of the non-associated flow rule, there is no evidence for the limit analysis theorems, but we proceed via analogy.
In the statically admissible solution, static conditions must be satisfied, but this requirement does not apply to kinematic conditions [30], which allows determining lower stress values and searching for the solution giving the highest values [30].The solution is based on the analysis of separated slices of, for example, a slope and the determination of their interaction by checking the state of stress.In each slice, the values of forces (loads) are determined, as a result of which plasticization of the material can occur based on, for example, the Mohr-Coulomb condition.The conditions that the surfaces (lines) of discontinuous stresses must meet are derived from the equilibrium conditions (Figure 2).
theorems, but we proceed via analogy.
In the statically admissible solution, static conditions must be satisfied, but this re quirement does not apply to kinematic conditions [30], which allows determining lowe stress values and searching for the solution giving the highest values [30].The solution based on the analysis of separated slices of, for example, a slope and the determination o their interaction by checking the state of stress.In each slice, the values of forces (load are determined, as a result of which plasticization of the material can occur based on, fo example, the Mohr-Coulomb condition.The conditions that the surfaces (lines) of di continuous stresses must meet are derived from the equilibrium conditions (Figure 2).From this, it follows that the stresses on both sides of the discontinuity surface mu be the same in such a structure.The stress parallel to this line may experience a spike.
On the other hand, the upper limit of load values can be determined on the basis o kinematically admissible mechanisms arbitrarily adopted for the system.To apply kinematically admissible solution, the mechanism adopted for consideration must mee the following criteria:

•
The continuity of the medium must be achieved.

•
The rule of flow must be assumed.
Static conditions need not be satisfied in this case.
In addition, since this is not an exact solution and depends on the kinematic mech anism adopted by the researcher, the solution with the smallest value of estimated force (loads) that is "closer" to the exact solution should be searched for.The kinematicall admissible solution (the number of mechanisms is practically infinite) assumes that in th medium under consideration, deformation occurs along the lines of discontinuity sepa rating areas of different kinematics defined by the researcher.
Consider an element undergoing plane deformation, which consists of the move ment with two rigid blocks moving along the discontinuity line L.
In the case of a Mohr-Coulomb medium and the associated flow rule, the velocit increase vector along the discontinuity line deviates from it by an angle φ (internal fri tion) (Figure 3a).From this, it follows that the stresses on both sides of the discontinuity surface must be the same in such a structure.The stress parallel to this line may experience a spike.
On the other hand, the upper limit of load values can be determined on the basis of kinematically admissible mechanisms arbitrarily adopted for the system.To apply a kinematically admissible solution, the mechanism adopted for consideration must meet the following criteria:

•
The continuity of the medium must be achieved.

•
The rule of flow must be assumed.
Static conditions need not be satisfied in this case.
In addition, since this is not an exact solution and depends on the kinematic mechanism adopted by the researcher, the solution with the smallest value of estimated forces (loads) that is "closer" to the exact solution should be searched for.The kinematically admissible solution (the number of mechanisms is practically infinite) assumes that in the medium under consideration, deformation occurs along the lines of discontinuity separating areas of different kinematics defined by the researcher.
Consider an element undergoing plane deformation, which consists of the movement with two rigid blocks moving along the discontinuity line L.
In the case of a Mohr-Coulomb medium and the associated flow rule, the velocity increase vector along the discontinuity line deviates from it by an angle ϕ (internal friction) (Figure 3a).In the case of a Mohr-Coulomb medium and the non-associated flow rule sume the Tresca potential), the vector of the velocity increase on the discontinui tangent to it (Figure 3).
It can be demonstrated that the unit energy dispersion  ̇ on the discontinui determined by the correlation: In the case of a Mohr-Coulomb medium and the non-associated flow rule (we assume the Tresca potential), the vector of the velocity increase on the discontinuity line is tangent to it (Figure 3).
It can be demonstrated that the unit energy dispersion .D on the discontinuity line is determined by the correlation: where c is cohesion for the associated flow rule, and: .
for the non-associated flow rule.
In the general case, it is possible to distinguish some typical mechanisms relating to the adopted division of the embankment or slope, along with the adopted behavior of the separate parts and the planes of their mutual displacement, as described in more detail in [24,30].
Regardless of the mechanism adopted, a hodograph-a velocity plan incorporating velocity vectors related to a common pole-should be made in each case, using the properties of the velocity field [30].
Figure 4 shows three simple kinematically admissible mechanisms along with velocity hodographs for the associated flow rule.In the case of a Mohr-Coulomb medium and the non-associated flow rule sume the Tresca potential), the vector of the velocity increase on the discontinuit tangent to it (Figure 3).
It can be demonstrated that the unit energy dispersion  ̇ on the discontinuit determined by the correlation: where c is cohesion for the associated flow rule, and: for the non-associated flow rule.
In the general case, it is possible to distinguish some typical mechanisms rel the adopted division of the embankment or slope, along with the adopted behavi separate parts and the planes of their mutual displacement, as described in more [24,30].
Regardless of the mechanism adopted, a hodograph-a velocity plan incorp velocity vectors related to a common pole-should be made in each case, u properties of the velocity field [30].
Figure 4 shows three simple kinematically admissible mechanisms along w locity hodographs for the associated flow rule.The solution for determining the force required to implement the adopted kinematic mechanism can be found by applying the block equilibrium method or the energy method.
The first method involves determining the forces acting on the various parts of the considered system, and the graphical or mathematical solution is based on the estimation of the unknown force values.The energy method involves determining the maximum force based on the energy balance of the system under consideration.
An example of the application of kinematically admissible methods can be found in the work of [1], which analyzed the passive pressure acting on the retaining wall, which results from the operation of machinery during earthworks, as well as in the work of [31], in which the stability of slopes was estimated based on the analysis of the equilibrium of curved slices.

Statically and Kinematically Admissible Solution-Simple Example
As mentioned in Section 2, to calculate the upper and lower bounds of the exact solution, both statically and kinematically admissible solutions can be constructed.The concept of such simple solutions for a vertical slope is presented below.

Statically Admissible Solution-Lower Bound
It is assumed that in the areas shown in Figure 5, the state of stress is determined by Equations ( 7)- (11).These relationships must be statically admissible and meet the conditions of discontinuity of stress-dashed lines indicate the lines of discontinuity of stress (Figure 2).

method.
The first method involves determining the forces acting on the various parts o considered system, and the graphical or mathematical solution is based on the estim of the unknown force values.The energy method involves determining the maxi force based on the energy balance of the system under consideration.
An example of the application of kinematically admissible methods can be fou the work of [1], which analyzed the passive pressure acting on the retaining wall, w results from the operation of machinery during earthworks, as well as in the work of in which the stability of slopes was estimated based on the analysis of the equilibriu curved slices.

Statically and Kinematically Admissible Solution-Simple Example
As mentioned in Section 2, to calculate the upper and lower bounds of the exa lution, both statically and kinematically admissible solutions can be constructed concept of such simple solutions for a vertical slope is presented below.

Statically Admissible Solution-Lower Bound
It is assumed that in the areas shown in Figure 5, the state of stress is determine Equations ( 7)- (11).These relationships must be statically admissible and meet the c tions of discontinuity of stress-dashed lines indicate the lines of discontinuity of s (Figure 2).There is no load along line AB and in Areas I and IV, it was assumed that   = 0. Areas III, II, and V, horizontal stress   has to be equal and it was assumed that   ( − ).In all areas, stresses are statically admissible, and conditions of discontinu stress between areas are fulfilled.
Area I: Area IV: Area III: There is no load along line AB and in Areas I and IV, it was assumed that σ x = 0.In Areas III, II, and V, horizontal stress σ x has to be equal and it was assumed that σ x = γ(y − L).In all areas, stresses are statically admissible, and conditions of discontinuity of stress between areas are fulfilled.
Area I: Area IV: Area III: Area II: Area V: It can be noticed that Area I gives the most unsafe results.
The results of the use of the Mohr-Coulomb failure criterion are shown in Figure 6. (σ where c is cohesion, and ϕ is the internal friction angle.The safety factor for unloaded slope q = 0.
where  is the allowable stress, and φ is the internal friction angle.
It can be demonstrated that the stress state in Areas III, II, and V for q < qdop, w qdop is an allowable load, does not reach an arbitrarily deep boundary condition (at in ity).The safety factor for unloaded slope q = 0.
where σ dop is the allowable stress, and ϕ is the internal friction angle.
It can be demonstrated that the stress state in Areas III, II, and V for q < q dop , where q dop is an allowable load, does not reach an arbitrarily deep boundary condition (at infinity).

Kinematically Admissible Solution-Associated Flow Rule-Upper Bound
Energy dissipation for the associated flow rule, where P = q•BC, and the energy equilibrium condition is defined by the relation: Hence, for the data adopted: q = 41.43 kN m The safety factor (for the unloaded slope) is defined as the ratio of the energy dispersed at the slip line to the energy of the sliding block.
To calculate the forces and the safety factor, the block balance method (Figure 7), equivalent to the energy method, or a simpler method, i.e., energy dissipation, can be used.
and the energy equilibrium condition is defined by the relation: Hence, for the data adopted: The safety factor (for the unloaded slope) is defined as the ratio of the energy persed at the slip line to the energy of the sliding block.
To calculate the forces and the safety factor, the block balance method (Figur equivalent to the energy method, or a simpler method, i.e., energy dissipation, ca used.

Kinematically Admissible Solution-Non-Associated Flow Rule-Lower Boun
From the energy equilibrium condition,  0 ( + ) =  1 •  Hence, for the data adopted: The safety factor for unloaded slope is: Table 1 summarizes the load values and safety factors estimated using the t calculation methods.

Kinematically Admissible Solution-Non-Associated Flow Rule-Lower Bound
From the energy equilibrium condition, Hence, for the data adopted: q = 18.52 kN m The safety factor for unloaded slope is: Table 1 summarizes the load values and safety factors estimated using the three calculation methods.The distinction between the non-associated (Figure 8) and associated flow rule (Figure 7) exists in the intrinsic property of a material and has to be established/assumed.Not knowing this material property, to describe the slope sliding process, both rules were applied, and safety factors were calculated for comparison.In the case of material described using a non-associated flow rule, a solution closer to a statically admissible solution is obtained, which is thus closer to a precise solution.

Safety factor
1.33 2.93 1.86 The distinction between the non-associated (Figure 8) and associated flow (Figure 7) exists in the intrinsic property of a material and has to be established/assu Not knowing this material property, to describe the slope sliding process, both rules applied, and safety factors were calculated for comparison.In the case of materia scribed using a non-associated flow rule, a solution closer to a statically admissible tion is obtained, which is thus closer to a precise solution.

Slope Stability-The Fellenius Method
As the first computational example, the bearing capacity of a slope was analy inclined at an angle of 60° and built of a single type of soil (clayed sand), with ang internal friction φ = 20°, cohesion c = 20 kPa, and volumetric density  = 2.15 g cm 3 .Figure 9 shows a diagram of the analyzed slope, where, according to the Felle method, the point of theoretical rotation O is marked, with the following coordinate 0.25 m and y = 7.75 m, and the cylindrical slip plane is plotted.Figure 9 shows a diagram of the analyzed slope, where, according to the Fellenius method, the point of theoretical rotation O is marked, with the following coordinates x = 0.25 m and y = 7.75 m, and the cylindrical slip plane is plotted.
The distinction between the non-associated (Figure 8) and associated flow (Figure 7) exists in the intrinsic property of a material and has to be established/assu Not knowing this material property, to describe the slope sliding process, both rules applied, and safety factors were calculated for comparison.In the case of materia scribed using a non-associated flow rule, a solution closer to a statically admissible tion is obtained, which is thus closer to a precise solution.

Slope Stability-The Fellenius Method
As the first computational example, the bearing capacity of a slope was analy inclined at an angle of 60° and built of a single type of soil (clayed sand), with ang internal friction φ = 20°, cohesion c = 20 kPa, and volumetric density  = 2.15 g cm 3 .Figure 9 shows a diagram of the analyzed slope, where, according to the Felle method, the point of theoretical rotation O is marked, with the following coordinate 0.25 m and y = 7.75 m, and the cylindrical slip plane is plotted.The area thus defined was divided into 11 (Figure 10) slices with a width B of no more than 0.73 m according to Equation (1): where B is the width of the block (slice); R is the radius of rotation of the block.
Hence, the slope safety factor F (7) is: Thus, the slope can be considered stable.The values of the forces acting on each block (slice) are determined using Equations ( 2)-( 5).
where G i is the weight of a single block, P i is the area of a separate block, and γ is the volumetric weight of the soil (21.9 kN/m 3 ).
N i is the normal component.cos α i is the cosine of angle alpha.
c is cohesion; l i is the length the section of contact between the block and the base; tgφ is the tangent of the internal friction angle.
The holding moment M u and the moment of causing rotation M o are given by the following correlation: Hence, the slope safety factor F (7) is: Thus, the slope can be considered stable.

Kinematically Admissible Solution-One Slip Line Associated Flow Rule
An analysis was also conducted for a simpler case, where the slope landslide was modeled as a single block (for which the parameters are summarized in Table 2).In this case, too, both associated and non-associated flow rules were used.For the associated flow rule, the distribution of forces and velocities and hodograph are shown in Figures 11 and 12.
Materials 2023, 16, x FOR PEER REVIEW 3.2.2.Kinematically Admissible Solution-One Slip Line Associated Flow Rule An analysis was also conducted for a simpler case, where the slope landsl modeled as a single block (for which the parameters are summarized in Table 2 case, too, both associated and non-associated flow rules were used.3.2.2.Kinematically Admissible Solution-One Slip Line Associated Flow R An analysis was also conducted for a simpler case, where the slope modeled as a single block (for which the parameters are summarized in T case, too, both associated and non-associated flow rules were used.Hence: 3.2.3.Kinematically Admissible Solution-One Slip Line (Non-Associated Flow Rule) For the non-associated flow rule, the velocity distribution and the hodograph are shown in Figures 13 and 14.
Materials 2023, 16, x FOR PEER REVIEW 3.2.3.Kinematically Admissible Solution-One Slip Line (Non-Associated Flow R For the non-associated flow rule, the velocity distribution and the hodogr shown in Figures 13 and 14.

Kinematically Admissible Solution for an Associated Flow Rule
The stability analysis of the selected slope was also carried out by assum rigid blocks as the failure mechanism.The area at risk of the landslide was assum as in the Fellenius method, entering the rigid area-the triangle ∆ABC and ∆ACD Figures 15 and 16 show the slip lines and velocities of each block and th acting on the slip lines.The stability analysis of the selected slope was also carried out by assum rigid blocks as the failure mechanism.The area at risk of the landslide was assum as in the Fellenius method, entering the rigid area-the triangle ∆ABC and ∆AC Figures 15 and 16 show the slip lines and velocities of each block and th acting on the slip lines.Hence:

Kinematically Admissible Solution for an Associated Flow Rule
The stability analysis of the selected slope was also carried out by assuming two rigid blocks as the failure mechanism.The area at risk of the landslide was assumed to be as in the Fellenius method, entering the rigid area-the triangle ∆ABC and ∆ACD.
Figures 15 and 16 show the slip lines and velocities of each block and the forces acting on the slip lines.The geometric parameters of the adopted fields are given in Table 3. Table 3. Geometric parameters of the adopted fields.The geometric parameters of the adopted fields are given in Table 3.

Block
Table 3. Geometric parameters of the adopted fields.The geometric parameters of the adopted fields are given in Table 3.

Block
Table 3. Geometric parameters of the adopted fields.The geometric parameters of the adopted fields are given in Table 3.

Block
Hence: 3.2.6.Kinematically Admissible Solution-Associated Flow Rule The stability analysis of the selected slope was also carried out using a system of four blocks in the rigid area-the triangles ∆ABC, ∆ACD, ∆ADE, and ∆AEF.
Figure 20 shows the slip lines and velocities of each block, as well as the forces acting on the slip lines, whose geometric parameters for the associated flow rule are given in Hence:

Kinematically Admissible Solution-Associated Flow Rule
The stability analysis of the selected slope was also carried out using a system of four blocks in the rigid area-the triangles ∆ABC, ∆ACD, ∆ADE, and ∆AEF.
Figure 20 shows the slip lines and velocities of each block, as well as the forces acting on the slip lines, whose geometric parameters for the associated flow rule are given in Table 4.  Figure 21 shows the velocity hodograph.The energy method was used to determine the safety of the slope.Figure 21 shows the velocity hodograph.The energy method was used to determine the safety of the slope.The energy method was used to determine the safety of the slope.

𝑛𝑛
3.2.7.The Kinematically Admissible Solution-The Associate Flow Rule An analysis was also performed for the non-associated flow rule.Figure 22 shows the slip and velocity lines of each block and the forces acting on the slip lines.An analysis was also performed for the non-associated flow rule.Figure 22 show the slip and velocity lines of each block and the forces acting on the slip lines.The energy method was used to determine the safety of the slope.The energy method was used to determine the safety of the slope.The energy method was used to determine the safety of the slope.
A summary of the results is presented in the Table 5.As mentioned before, the distinction between the non-associated and associated flow rule exists in the intrinsic property of a material and has to be established/assumed.Not knowing such material property, to describe the slope sliding process, both rules were applied, and safety factors were calculated for comparison.
In this paper, where the Mohr-Coulomb medium is discussed, we assume the Tresca potential for the non-associated flow rule.It is worth mentioning that potentials other than the Tresca potential can be taken into account to describe material behavior.
In the case of the Fellenius method, the safety factor was built in terms of the moment ratio, since the "block rotation" mechanism was assumed, and the comparison of the holding moment and moment-rotation is most appropriate.
In the case of kinematically admissible mechanisms, the mechanism of sliding blocks is assumed, and the comparison of work rate was used to calculate the safety factor.

Conclusions
This paper presents the application of simple plasticity theory solutions that give a basic assessment of slope stability and bearing capacity for simple cases.Using very simple kinematically admissible solutions, it is easy to assess the stability or bearing capacity of a slope.The calculations presented in Section 3.1 show that for a material described using the non-associated flow rule, the solution is close to that obtained with the Fellenius method.This way one can have a simple and quick rough estimation of slope stability.
Even a very simple kinematically admissible mechanism (one sliding block), inscribed in the landslide area determined with the Fellenius method, can be taken into consideration (Table 5).
The authors plan to carry on further research on this subject, including experimental studies.

Figure 1 .
Figure 1.Velocity grids and the corresponding hodograph for the problem of wall pressure in a Coulomb medium for (a) associated flow rule and (b) non-associated flow rule.

Figure 1 .
Figure 1.Velocity grids and the corresponding hodograph for the problem of wall pressure in a Coulomb medium for (a) associated flow rule and (b) non-associated flow rule.

Figure 1 .
Figure 1.Velocity grids and the corresponding hodograph for the problem of wall pressure in a Coulomb medium for (a) associated flow rule and (b) non-associated flow rule.

Figure 1 .
Figure 1.Velocity grids and the corresponding hodograph for the problem of wall pressure in a Coulomb medium for (a) associated flow rule and (b) non-associated flow rule.

Figure 3 .
Figure 3. Line of discontinuity of velocity: (a) associated flow rule and (b) non-associated

Figure 3 .
Figure 3. Line of discontinuity of velocity: (a) associated flow rule and (b) non-associated flow rule.

Figure 3 .
Figure 3. Line of discontinuity of velocity: (a) associated flow rule and (b) non-associated f

Figure 7 .
Figure 7. Diagram of force distribution according to the associated flow rule.

Figure 7 .
Figure 7. Diagram of force distribution according to the associated flow rule.

Figure 8 .
Figure 8. Diagram of the distribution of forces according to the non-associated flow rule.

Figure 9 .
Figure 9. Diagram of the analyzed slope.

Figure 8 .
Figure 8. Diagram of the distribution of forces according to the non-associated flow rule.

3. 2 .
Simple Kinematically Admissible Solutions for the Slope Stability Problem-Comparison with the Classical Fellenius Method 3.2.1.Slope Stability-The Fellenius Method As the first computational example, the bearing capacity of a slope was analyzed, inclined at an angle of 60 • and built of a single type of soil (clayed sand), with angle of internal friction ϕ = 20 • , cohesion c = 20 kPa, and volumetric density ρ = 2.15 g cm 3 .

Figure 8 .
Figure 8. Diagram of the distribution of forces according to the non-associated flow rule.

Figure 9 .
Figure 9. Diagram of the analyzed slope.

Figure 9 .
Figure 9. Diagram of the analyzed slope.

Figure 10 .
Figure 10.(a) Division of the slope into slices with a diagram of the forces acting on a slice; (b) distribution of forces in a single slice.

Figure 10 .
Figure 10.(a) Division of the slope into slices with a diagram of the forces acting on a slice; (b) distribution of forces in a single slice.

= 1 For
the associated flow rule, the distribution of forces and velocities and hod are shown in Figures 11 and 12.

Figure 11 .
Figure 11.Distribution of forces and velocities acting on the slip line (associated flow rule

Figure 11 .
Figure 11.Distribution of forces and velocities acting on the slip line (associated flow rule).
For the associated flow rule, the distribution of forces and velocities an are shown in Figures11 and 12.

Figure 11 .
Figure 11.Distribution of forces and velocities acting on the slip line (associated flo

Figure 13 .
Figure 13.Distribution of velocities acting in the slope according to the adopted kine admissible solution (non-associated flow rule).

Figure 13 .
Figure 13.Distribution of velocities acting in the slope according to the adopted kinematically admissible solution (non-associated flow rule). .

Figure 13 .
Figure 13.Distribution of velocities acting in the slope according to the adopted kine admissible solution (non-associated flow rule).

Figure 15 .
Figure 15.Dividing the slope into blocks for a kinematically admissible solution.

Figure 16 .
Figure 16.Distribution of forces acting in the slope according to the adopted kinematically admissible solution.

Figure 17
Figure17shows the hodograph where for the associated flow rule α = 33.4°.

Figure 17 .
Figure 17.Velocity hodograph for associated flow rule with two rigid blocks.

Figure 15 .Figure 15 .
Figure 15.Dividing the slope into blocks for a kinematically admissible solution.

Figure 16 .
Figure 16.Distribution of forces acting in the slope according to the adopted kinematicall sible solution.

Figure 17 shows
Figure 17 shows the hodograph where for the associated flow rule α = 33.4°.

Figure 17 .
Figure 17.Velocity hodograph for associated flow rule with two rigid blocks.

Area Pi 2 Side
Lengths Li Weight G = Pi•γ Tangential Compo

Figure 16 .
Figure 16.Distribution of forces acting in the slope according to the adopted kinematically admissible solution.

Figure 17
Figure 17 shows the hodograph where for the associated flow rule α = 33.4• .

Figure 15 .
Figure 15.Dividing the slope into blocks for a kinematically admissible solution.

Figure 16 .
Figure 16.Distribution of forces acting in the slope according to the adopted kinem sible solution.

Figure 17 Figure 17 .
Figure17shows the hodograph where for the associated flow rule α =

Figure 17 .
Figure 17.Velocity hodograph for associated flow rule with two rigid blocks.

5 .
Figure 18 shows the velocity distribution on the BC, AC, and DC slip lines, and Figure 19 shows the hodograph.

Figure 18 .
Figure 18.Distribution of velocities acting in the slope.Figure 18. Distribution of velocities acting in the slope.

Figure 18 . 20 Figure 19 .
Figure 18.Distribution of velocities acting in the slope.Figure 18. Distribution of velocities acting in the slope.Materials 2023, 16, x FOR PEER REVIEW 15 of 20

Figure 19 .
Figure 19.Velocity hodograph for non associated flow rule with two rigid blocks.

Figure 20 .
Figure 20.Distribution of forces acting in the slope according to the adopted kinematically admissible solution.

Figure 20 .
Figure 20.Distribution of forces acting in the slope according to the adopted kinematically admissible solution.

Figure 21
Figure 21 shows the velocity hodograph.

Figure 20 .
Figure 20.Distribution of forces acting in the slope according to the adopted kinematically sible solution.

Figure 21 .
Figure 21.Velocity hodograph for associated flow rule with tree rigid blocks.

7 .
The Kinematically Admissible Solution-The Associate Flow Rule

Figure 22 .
Figure 22.Distribution of forces acting in the slope according to the adopted kinematically adm sible solution.

Figure 23
Figure 23 shows the hodograph of velocity for the non-associated flow rule.

Figure 23 .
Figure 23.Velocity hodograph for non associated flow rule with tree rigid blocks.

Figure 22 .
Figure 22.Distribution of forces acting in the slope according to the adopted kinematically admissible solution.

Figure 23 7 .
Figure 23 shows the hodograph of velocity for the non-associated flow rule.

Figure 22 .
Figure 22.Distribution of forces acting in the slope according to the adopted kinematically a sible solution.

Figure 23
Figure 23 shows the hodograph of velocity for the non-associated flow rule.

Figure 23 .
Figure 23.Velocity hodograph for non associated flow rule with tree rigid blocks.

Figure 23 .
Figure 23.Velocity hodograph for non associated flow rule with tree rigid blocks.

Table 1 .
Load values and safety factors obtained using three calculation methods.

Table 2 .
Parameters of the separated block.

Table 2 .
Parameters of the separated block.

Table 2 .
Parameters of the separated block.

Table 4 .
Geometric parameters of adopted fields.

Table 5 .
Summary of the test results.