Jointed rock slope stability

Jointed rock slope stability

Product: Abaqus/Standard  

This example illustrates the use of the jointed material model in the context of geotechnical applications. We examine the stability of the excavation of part of a jointed rock mass, leaving a sloped embankment. This problem is chosen mainly as a verification case because it has been studied previously by Barton (1971) and Hoek (1970), who used limit equilibrium methods, and by Zienkiewicz and Pande (1977), who used a finite element model. This example also has been extended to study the slope stability of excavated soil medium with the same geometry, by using the Mohr-Coulomb plasticity model with and without the tension cutoff feature.

Geometry and model

The plane strain model analyzed is shown in Figure 1.1.6�1 together with the excavation geometry and material properties. The rock mass contains two sets of planes of weakness: one vertical set of joints and one set of inclined joints. We begin from a nonzero state of stress. In this problem this consists of a vertical stress that increases linearly with depth to equilibrate the weight of the rock and horizontal stresses caused by tectonic effects: such stress is quite commonly encountered in geotechnical engineering. The active “loading” consists of removal of material to represent the excavation. It is clear that, with a different initial stress state, the response of the system would be different. This illustrates the need of nonlinear analysis in geotechnical applications—the response of a system to external “loading” depends on the state of the system when that loading sequence begins (and, by extension, to the sequence of loading). We can no longer think of superposing load cases, as is done in a linear analysis.

Practical geotechnical excavations involve a sequence of steps, in each of which some part of the material mass is removed. Liners or retaining walls can be inserted during this process. Thus, geotechnical problems require generality in creating and using a finite element model: the model itself, and not just its response, changes with time—parts of the original model disappear, while other components that were not originally present are added. This example is somewhat academic, in that we do not encounter this level of complexity. Instead, following the previous authors' use of the example, we assume that the entire excavation occurs simultaneously.

Solution controls

The jointed material model includes a joint opening/closing capability. When a joint opens, the material is assumed to have no elastic stiffness with respect to direct strain across the joint system. Because of this, and also as a result of the fact that different combinations of joints may be yielding at any one time, the overall convergence of the solution is expected to be nonmonotonic. In such cases the use of *CONTROLS, ANALYSIS=DISCONTINUOUS is generally recommended to prevent premature termination of the equilibrium iteration process because the solution may appear to be diverging.

As the end of the excavation process is approached, the automatic incrementation algorithm reduces the load increment significantly, indicating the onset of failure of the slope. In such analyses it is useful to specify a minimum time step to avoid unproductive iteration.

For the nonassociated flow case UNSYMM=YES is used on the *STEP option. This is essential for obtaining an acceptable rate of convergence since nonassociated flow plasticity has a nonsymmetric stiffness matrix.

Results and discussion

In this problem we examine the effect of joint cohesion on slope collapse through a sequence of solutions with different values of joint cohesion, with all other parameters kept fixed. Figure 1.1.6�2 shows the variation of horizontal displacements as cohesion is reduced at the crest of the slope (point A in Figure 1.1.6�1) and at a point one-third of the way up the slope (point B inFigure 1.1.6�1). This plot suggests that the slope collapses if the cohesion is less than 24 kPa for the case of associated flow or less than 26 kPa for the case of nonassociated flow. These compare well with the value calculated by Barton (26 kPa) using a planar failure assumption in his limit equilibrium calculations. Barton's calculations also include “tension cracking” (akin to joint opening with no tension strength) as we do in our calculation. Hoek calculates a cohesion value of 24 kPa for collapse of the slope. Although he also makes the planar failure assumption, he does not include tension cracking. This is, presumably, the reason why his calculated value is lower than Barton's. Zienkiewicz and Pande assume the joints have a tension strength of one-tenth of the cohesion and calculate the cohesion value necessary for collapse as 23 kPa for associated flow and 25 kPa for nondilatant flow.

Figure 1.1.6�3 shows the deformed configuration after excavation for the nonassociated flow case and clearly illustrates the manner in which the collapse is expected to occur. Figure 1.1.6�4shows the magnitude of the frictional slip on each joint system for the nonassociated flow case. A few joints open near the crest of the slope.

The study of soil slope stability using the Mohr-Coulomb plasticity model is performed for two cases: one without tension cutoff and one including the tension cutoff feature. The tension cutoff feature limits the stress carrying capacity of soil in tension. It can be seen that the maximum principal stress without tension cutoff (see contour plot in Figure 1.1.6�5) is higher than the limiting maximum principal stress (see Figure 1.1.6�6) with tension cutoff as expected. With tension cutoff, one also observes the appearance of the equivalent plastic strain in tension, PEEQT in the region of maximum principal stress. In this case it is also seen that the equivalent plastic strain, PEEQ on the cohesion failure surface is higher compared to the case without tension cutoff. The contour plots for the equivalent plastic strains are not shown.

Input files

jointrockstabil_nonassoc_30pka.inp

Nonassociated flow case problem; cohesion = 30 kPa.

jointrockstabil_assoc_25kpa.inp

Associated flow case; cohesion = 25 kPa.

mc_slopestabil.inp

Slope stability analysis, Mohr-Coulomb plasticity without tension cutoff

mctc_slopestabil.inp

Slope stability analysis, Mohr-Coulomb plasticity with tension cutoff

References

Barton, N., “Progressive Failure of Excavated Rock Slopes,” Stability of Rock Slopes, Proceedings of the 13th Symposium on Rock Mechanics, Illinois, pp. 139–170, 1971.

Hoek,  E., “Estimating the Stability of Excavated Slopes in Open Cast Mines,” Trans. Inst. Min. and Metal., vol. 79, pp. 109�132, 1970.

Zienkiewicz,  O. C., and G. N. Pande, “Time-Dependent Multilaminate Model of Rocks – A Numerical Study of Deformation and Failure of Rock Masses,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 1, pp. 219�247, 1977.

Figures

Figure 1.1.6�1 Jointed rock slope problem.

Figure 1.1.6�2 Horizontal displacements with varying cohesion.

Figure 1.1.6�3 Deformed configuration (nonassociated flow).

Figure 1.1.6�4 Contours of frictional slip magnitudes (nonassociated flow).

Figure 1.1.6�5 Maximum principal stress without tension cutoff.

Figure 1.1.6�6 Maximum principal stress with tension cutoff.