See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/254542119

The Shear Strength Reduction Method for the Generalized Hoek-Brown

Criterion

Article · January 2005

CITATIONS READS

29 731

4 authors, including:

Reginald Hammah Thamer Yacoub

Golder Associates Rocscience
20 PUBLICATIONS   416 CITATIONS    35 PUBLICATIONS   252 CITATIONS   

SEE PROFILE SEE PROFILE

Brent Corkum
Rocscience Inc.
21 PUBLICATIONS   1,549 CITATIONS   

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Discrete Element Modeling of Particulate Material Behavior View project

Anisotropic Rock Masses and Directional Shear Strength Models View project

All content following this page was uploaded by Brent Corkum on 16 September 2015.

The user has requested enhancement of the downloaded file.

ARMA/USRMS 05-810

The Shear Strength Reduction Method for the

Generalized Hoek-Brown Criterion
Hammah, R.E., Yacoub, T.E. and Corkum, B.C.
Rocscience Inc., Toronto, ON, Canada
Curran, J.H.
Lassonde Institute, University of Toronto, Toronto, ON, Canada
Copyright 2005, ARMA, American Rock Mechanics Association

This paper was prepared for presentation at Alaska Rocks 2005, The 40th U.S. Symposium on Rock Mechanics (USRMS): Rock Mechanics for Energy, Mineral and Infrastructure
Development in the Northern Regions, held in Anchorage, Alaska, June 25-29, 2005.
This paper was selected for presentation by a USRMS Program Committee following review of information contained in an abstract submitted earlier by the author(s). Contents of the paper,
as presented, have not been reviewed by ARMA/USRMS and are subject to correction by the author(s). The material, as presented, does not necessarily reflect any position of USRMS,
ARMA, their officers, or members. Electronic reproduction, distribution, or storage of any part of this paper for commercial purposes without the written consent of ARMA is prohibited.
Permission to reproduce in print is restricted to an abstract of not more than 300 words; illustrations may not be copied. The abstract must contain conspicuous acknowledgement of where
and by whom the paper was presented.

ABSTRACT: This paper describes a method that allows direct use of the Generalized Hoek-Brown criterion in Finite Element
(FE) Shear Strength Reduction (SSR) analysis of rock slopes. It describes the generation of shear envelopes for the criterion,
lowering of the shear envelope by a factor, and means for determining an equivalent Generalized Hoek-Brown curve that best
approximates the lowered envelope. The paper provides two examples that demonstrate the performance of the new method. As
well it outlines some of the benefits of the SSR technique for slope stability analysis.

1. INTRODUCTION 2. OVERVIEW OF THE SHEAR STRENGTH

The Finite Element Method (FEM) is increasingly REDUCTION METHOD
being applied to slope stability analysis. One of the The SSR technique for slope stability analysis
most popular techniques for performing FEM slope involves systematic use of finite element analysis to
analysis is the shear strength reduction (SSR)
determine a stress reduction factor (SRF) or factor
approach [1]. The SSR is simple in concept:
of safety value that brings a slope to the verge of
systematically reduce the shear strength envelope of failure. The shear strengths of all the materials in a
material by a factor of safety, and compute FEM FE model of a slope are reduced by the SRF.
models of the slope until deformations are Conventional FE analysis of this model is then
unacceptably large or solutions do not converge. performed until a critical SRF value that induces
A factor limiting broader application of the SSR instability is attained. A slope is considered unstable
approach to slope stability analysis has been its in the SSR technique when its FE model does not
restriction to Mohr-Coulomb materials. Most converge to a solution (within a specified
discussions of the method found in literature deal tolerance).
with this criterion (the paper by Shukra and Baker As mentioned earlier, most existing descriptions
[2] is one of the few known to the authors that and discussions of the SSR technique are based on
examines application to non-linear strength use of the Mohr-Coulomb strength models for
envelopes). For rock masses, the Generalized Hoek- materials. The criterion is readily used in the SSR
Brown criterion is the most commonly applied technique for the following reasons:
strength model. As a result, the authors found it
expedient to develop an SSR framework for the (i) It can be expressed either in terms of
Hoek-Brown criterion. The aim of this paper is to principal stresses, or in terms of shear
outline an approach for applying the Generalized and normal stresses (this makes it
Hoek-Brown criterion in SSR analysis. The paper amenable for use in both FE and limit-
will also demonstrate the capability and accuracy of equilibrium analyses)
the proposed approach through two examples.
 (ii) Its linearity that allows reduced dσ 1
parameters to be readily calculated when dσ 3
τ = (σ 1 − σ 3 ) . (3)
an original shear strength model is dσ 1
+1
reduced by a factor F (Griffith and Lane dσ 3
[1] provide simple, closed-form
For the Generalized Hoek-Brown criterion, the
equations for calculating reduced
parameters), and following equations relate σ n and τ to σ 1 and σ 3 :

(iii) It is readily available in many existing ⎛ σ ⎞⎟a−1

finite element software. 1 + amb ⎜⎜⎜mb 3 + s⎟⎟
⎜⎝ σci ⎠⎟⎟
Determining the parameters of a Generalized Hoek- (
τ = σ1 − σ3 ) ⎛ σ ⎞⎟a−1
(4)
Brown model, which is equivalent to a shear 2 + amb ⎜⎜⎜mb 3 + s⎟⎟
⎜⎝ σci ⎠⎟⎟
envelope reduced by a factor of safety, is not as
straightforward. This paper will develop an
approach for estimating these equivalent
⎛ σ ⎞⎟a−1
parameters. ⎜
amb ⎜⎜mb 3
+ s⎟⎟
⎝⎜ σci ⎠⎟⎟
3. FE SLOPE STABILITY FOR GENERALIZED
σn =
1
2
( 1
)
σ1 + σ3 − σ1 − σ3
2
( ) ⎛ σ ⎞⎟a−1
(5)
⎜
2 + amb ⎜⎜mb 3
+ s⎟⎟
⎠⎟⎟
HOEK-BROWN MATERIALS ⎜⎝ σci
The Generalized Hoek-Brown criterion [3] for rock
The SSR form of analysis involves the following
masses is non-linear, and defines material strength
steps:
in terms of major and minor principal stresses
through the equation (i) Reduction of the shear strength envelope
by a factor F
⎛ ⎞⎟a
⎜⎜ σ3 ⎟
σ1 = σ3 + σci ⎜⎜mb + s⎟⎟⎟ (1) (ii) Determination of new strength model
⎜⎜ σci ⎟ parameters that conform to the lowered
⎝ ⎠⎟
envelope, and
where σci is the uniaxial compressive strength of (iii) Use of the new parameters in
the intact rock material, while conventional FE elasto-plastic analysis.
⎛ GSI −100 ⎞⎟ ⎛ GSI −100 ⎞⎟
mb = mi exp ⎜⎜ , s = exp ⎜⎜⎜ To lower the Generalized Hoek-Brown shear
⎜⎝ 28 −14 D ⎠⎟⎟ ⎟, and
⎝ 9 − 3D ⎠⎟ strength envelope by the factor F, we simply divide
1 1
(
a = + e−GSI /15 − e−20 / 3 .
2 6
) Equation (4) by F.
⎛ σ ⎞⎟a−1
mi is an intact rock material property, GSI is the 1 + amb ⎜⎜⎜mb 3 + s⎟⎟
τ orig ⎜⎝ σci ⎠⎟⎟
geological strength index, and D is the disturbance
factor [3].
τ red =
F
(
= σ1 − σ3 ) ⎛ σ ⎞⎟ a−1
⋅
1
F
2 + amb ⎜⎜⎜mb 3 + s⎟⎟
Using relationships developed by Balmer [4, 3], a ⎜⎝ σci ⎠⎟⎟
shear-normal stress envelope equivalent to the (6)
red
−1
⎛ σ3 ⎞⎟a
Generalized Hoek-Brown principal stress envelope 1 + ambred ⎜⎜⎜mbred red + s red ⎟⎟
can be determined. The shear and normal stress pair ⎜⎝ σci ⎠⎟⎟
corresponding to a point on a principal stress (
= σ1 − σ3 ) red
−1
⎛ ⎞a
envelope can be determined from the following 2 + ambred ⎜⎜m red σ3 + s red ⎟⎟
⎜⎜ b ⎟
equations: ⎝ σcired ⎠⎟⎟

dσ 1 The geometric interpretation of lowering the shear

−1
1 1 dσ 3
σ n = (σ 1 + σ 3 ) − ( σ 1 − σ 3 ) (2) envelope by a factor is illustrated on Figure 1. It
2 2 dσ 1 shows the original shear envelope, τ orig , which has
+1
dσ 3 been reduced by F to produce a reduced envelope
τ red .
In the above equation, σ cired , mbred , s red , and a red , are approximates the reduced envelope. For any given
the parameters of the lowered (reduced) strength σ n value, the square of the error between the
envelope. σ 1 can be eliminated from the equation reduced and approximated envelopes is defined by
by replacing it with Equation (1). the equation
ε (σ n ) = (τ appr − τ red ) .
2 2
It is important to note that given a set of (7)
Generalized Hoek-Brown parameters, and a
specified σ 3 value, σ n can be determined from
0 .18
Equation (5) through replacement of σ 1 with the
definition of the criterion (Equation (1)).
0 .16
Examination of Equation (6) reveals that the Original shear envelope
determination of the parameters of the reduced after reduction
envelope is not trivial. In a previous paper [5], the 0 .14
authors presented a simplified approach that
approximated the reduced Generalized Hoek-Brown
shear envelope with a linear Mohr-Coulomb 0 .12
equivalent. Although a useful first approximation,
this approach was quite sensitive to the range of
0.1
normal stresses over which the linear approximation Approximate GHB
was determined. A new approach, which actually shear envelope
determines a best-fit Generalized Hoek-Brown τ 0 .08
strength model to the reduced shear envelope, is
described next.
Differences (error)
0 .06
between the curves
0.45

0.4 0 .04
Shear envelope for
original GHB criterion
0.35 σt
0 .02
0.3

σnmax
0.25 0
τ

-0.0 5 0.0 5 0.1 5 0.2 5 0.35 0 .45

0.2
σn
0.15
Envelope after τ
reduction by F Figure 2. Reduction of a Generalized Hoek-Brown shear
envelope by a factor results in a lowered curve that generally
0.1
can no longer be described by a Hoek-Brown curve. The
figure shows an approximation of the reduced curve with an
0.05 τ/F
equivalent Hoek-Brown. Notice the regions of error or
differences between the two curves.
0
-0.05 0.05 0.15 0.25 0.35

σn The total error of the fit of τ appr to τ red can be

obtained through integration of the squared error
Figure 1. A Generalized Hoek-Brown criterion drawn in
shear-normal stress space, and the resulting curve when the function
envelope is reduced by a factor F. σ n max

Total error = ∫ ε (σ n ) d σ n
2
σt
(8)

3.1. Estimating the parameters of the reduced over the range σ t (the tensile strength) to a
shear strength envelope maximum normal stress value, σ n max . Because the
Figure 2 shows the reduced shear envelope, τ red , squared error function does not explicitly relate σ n
and a new Generalized Hoek-Brown, τ appr , that
to τ , the integration is best performed using a which is performed numerically) using a
numerical approach such as gaussian quadrature. technique such as the Simplex method.
The variables of the function are mi and
The parameters of the best-fit Generalized Hoek-
Brown envelope to the reduced shear strength GSI . σ cired and D have the fixed values
envelope can be obtained through minimization of described above.
the total squared error. This minimization is best
attained with techniques such as the Simplex 4. EXAMPLES
method, which do not require derivatives of the
function being minimized. We illustrate the capabilities of the above-outlined
FE SSR technique for the Generalized Hoek-Brown
criterion on two examples. The technique was
3.2. Algorithm for computing reduced Generalized implemented in the finite element program Phase2
Hoek-Brown parameters [6]. The particular implementation tested in this
Based on the discussions above, the authors paper assumes elastic-fully plastic material
developed an algorithm for determining the behaviour. (This condition can be easily relaxed
parameters of a curve that best fit a Generalized though, and in a future paper the authors will
Hoek-Brown shear strength envelope, which has discuss the impact of different elasto-plastic
been reduced by a factor, F. To reduce the number assumptions on results.)
of parameters to be determined, it is assumed that For Example 1, we compare the factor of safety
the uniaxial compressive strength, σ cired , of the value computed by the FE SSR to those obtained
reduced curve can be simply calculated as from the equivalent Mohr-Coulomb approach
σ ci described in [5], and from conventional limit-
σ cired = . (9) equilibrium analysis [7]. In the second example we
F
compare the FE SSR result to values computed
This assumption simplifies curve-fitting procedures from limit-equilibrium analysis.
considerably, but introduces practically no
additional error.
Next, instead of directly fitting for the parameters 4.1. Example 1
mbred , s red , and a red , the procedure assumes the
Example 1 involves analysis of a 10 m high
homogeneous rock slope with a 45o slope angle
disturbance parameter D = 0 , and estimates values
(Figure 1). The Generalized Hoek-Brown
for mi and GSI. (As in the case of σ cired , assuming parameters of the slope rock mass are provided in
D = 0 simplifies calculations substantially at very Table 1. Stresses in the slope are assumed
minimal penalty to accuracy.) These are then used gravitational, with a horizontal to vertical stress
to calculate values for mbred , s red , and a red . ratio of 1.
The steps for estimating the Generalized Hoek- Table 1. Properties of the rock mass in the Example 1slope
Brown parameters of the reduced shear envelope Property Value
are then as follow:
Young’s modulus, E (MPa) 5000
(i) Establish the range of minor principal
Poisson’s ratio, ν 0.3
stresses acting in a slope. Since the
minimum stress is taken to be the tensile Weight, γ (MN/m3) 0.025
strength, σ t , it is only necessary to Uniaxial compressive strength 30
σci (MPa)
determine the maximum σ 3 value in the
GSI 5
slope. Intact rock parameter mi 2
(ii) Determine the corresponding value of Disturbance factor, D 0
normal stress, σ n max , using Equation Parameter mb 0.067
(5). Parameter s 2.5 x 10-5
Parameter a 0.619
(iii) Minimize the squared error function
over the range [σ t , σ n max ] (integration of
 different strength type is also present in a slope. The
slope in Example 1, but this time with a horizontal
layer of Mohr-Coulomb material passing through
the toe (shown in Figure 5), was analyzed. The
Mohr-Coulomb layer had zero cohesion, a 25 o
friction angle, and a thickness of 1m.

Figure 3. Geometry of the slope in Example 1.

Table 2 below shows the factor of safety result
obtained from the FE SSR method, and compares it
to that calculated by the equivalent Mohr-Coulomb
approximation described in [5]. It also shows the
factor of safety values calculated by the Bishop
simplified and Spencer limit-equilibrium methods.
Figure 5. Geometry of slope in Example 2.
The SSR method for the Generalized Hoek-Brown
strength model gave near identical results to those
given by the other methods. Its predictions of the
failure mechanism (slip surface) are shown on
Figure 3.
Table 2. Comparison of factor of safety results for Example 1.
Method Factor of Safety
FE SSR technique:
Generalized Hoek-Brown 1.15
Equivalent Mohr-Coulomb 1.15
Figure 6. Contours of maximum shear strain in Example 2
Limit equilibrium: slope. Note the failure mechanism passing through the toe of
Bishop’s simplified 1.153 the slope.
Spencer’s method 1.152
The factor of safety value obtained from SSR
analysis of the slope is compared to those of non-
circular Bishop simplified and Spencer limit-
equilibrium analysis in Table 3. Again the SSR
method for the Generalized Hoek-Brown is very
close to those computed from limit-equilibrium
analysis. Also the failure mechanism indicated by
the SSR method, which is shown on Figure 6, is as
expected.
Table 3. Comparison of factor of safety results for Example 2.

Figure 4. Contours of maximum shear strain in the slope at

Method Factor of Safety
failure. The contours reveal the failure mechanism predicted FE SSR technique:
by the SSR method.
Generalized Hoek-Brown 0.95
Limit equilibrium:
4.2. Example 2 Bishop’s simplified 0.934
Example 2 was selected to test the performance of Spencer’s method 0.963
the proposed approach when a material of a
5. CONCLUSION REFERENCES
The development of an SSR framework that allows 1. Griffiths, D.V., and P.A. Lane. 1999. Slope stability
direct modeling of Generalized Hoek-Brown analysis by finite elements. Geotechnique, vol. 49, no.
3, pp. 387-403.
materials allows the rock engineering community to
more fully exploit the advantages and power of 2. Shukra, R. and R. Baker. 2003. Mesh geometry effects
on slope stability calculation by FLAC strength
finite element analysis.
reduction method – linear and non-linear failure
The FE SSR method is a robust alternative to limit- criteria. In Proceedings of the 3rd International FLAC
equilibrium slope stability methods. It is Symposium, Sudbury, Ontario, Canada, eds. R.
Brummer et al, pp. 109-116.
particularly beneficial in situations in which stress
has a dominant influence on stability. It is able to 3. Hoek E., C. Carranza-Torres, and B. Corkum. 2002.
Hoek-Brown criterion – 2002 edition. In Proceedings
calculate the deformations of reinforcement
of the 5th North American Rock Mechanics Symposium
elements such as rock bolts, anchors, and piles, as and the 17th Tunnelling Association of Canada:
well as important quantities such as bending NARMS-TAC 2002, Toronto, Canada, eds. R.E.
moments. Hammah et al, Vol. 1, pp. 267-273.

Like limit-equilibrium methods, the SSR technique 4. Balmer G. 1952. A general analytical solution for
Mohr’s envelope. American Society for Testing and
can accommodate multiple material layers, phreatic
Materials, vol. 52, pp. 1260-1271.
surfaces and seepage results. Unlike its limit-
equilibrium counterparts, it does not require a priori 5. Hammah, R.E., J.H. Curran, T.E. Yacoub, and B.
Corkum. 2004. Stability analysis of rock slopes using
assumptions on failure mechanisms (the shapes of the Finite Element Method. In Proceedings of the ISRM
failure surfaces). Regional Symposium EUROCK 2004 and the 53rd
Geomechanics Colloquy, Salzburg, Austria.
When contour plots of stresses and displacements,
such as those shown in this paper, are arranged in 6. Rocscience Inc. 2005. Phase2 v6.0 – a two-dimensional
sequence (for ordered factor of safety values finite element analysis program.
ranging from stable to unstable) they provide 7. Rocscience Inc. 2003. Slide v5.0 – a slope stability
insightful information on the development of failure program based on limit-equilibrium analysis.
mechanisms.
In addition to allowing rock engineers to harness the
above-listed advantages, the Generalized Hoek-
Brown formulation of the SSR technique makes it
possible to analyze limit-equilibrium slope models
involving non-linear material strength envelopes
such as the power curve criterion. Such strength
criteria, popular for soil materials, are difficult to
use in elasto-plastic finite element analysis due to
the absence of flow rules.
Using techniques similar to that used to fit a
Generalized Hoek-Brown curve to a reduced shear
envelope, equivalent Hoek-Brown parameters can
be determined for a non-linear strength envelope.
SSR analysis can then be performed on a slope
model using the equivalent Generalized Hoek-
Brown envelope.
Given the benefits of the SSR technique, the hope
authors hope that it will be applied more frequently
to rock slope problems.

View publication stats