GEO3701 Unit 3
GEO3701 Unit 3
GEO3701 Unit 3
7.3 Stiffness and Shear Strength of Discontinuity Planes ............................................... 25
7.3.1 Stiffness of a Discontinuity ................................................................................ 26
7.3.2 Clarification of Planar Smooth, Planar Rough and Dilation ............................... 26
7.3.3 Friction Angles of Rock (Basic, Peak and Residual) ............................................ 28
7.3.4 Shear Strength of Planar Smooth Discontinuity ................................................ 29
7.3.5 Shear Strength of Planar Rough Discontinuity .................................................. 30
7.3.6 Barton‐Choubey Failure Criterion ...................................................................... 31
7.3.7 Summary of Mohr‐Coulomb (Planar Smooth and Rough) and Barton‐Choubey
Failure Criteria .................................................................................................................. 33
8 Rock – As Intact Material .................................................................................................. 34
8.1 Differentiating between Isotropic and Anisotropic Rock .......................................... 34
8.2 Uniaxial Compression: σ‐ε Plot ................................................................................. 35
8.3 Four Stages of σ‐ε (Load‐Deformation) ..................................................................... 35
8.4 Factors Affecting UCS of Rock ................................................................................... 37
8.5 Isotropic Rock Strength Criteria ................................................................................ 37
8.5.1 Mohr‐Coulomb Shear Strength Criterion .......................................................... 37
8.5.2 Hoek‐Brown Failure Criterion ............................................................................ 39
8.6 Anisotropic Rock ........................................................................................................ 42
8.7 Time Dependant Behaviour ...................................................................................... 44
9 Rock – As a Mass ............................................................................................................... 44
9.1 Generalised Hoek‐Brown Peak Strength Criterion for Jointed Rock Mass ............... 47
9.1.1 Rock Mass Strength ........................................................................................... 47
Shear Strength Expressed in Terms of Principal Stresses ........................... 48
9.1.2 Shear Strength Expressed in Terms of Mohr‐Coulomb Failure Criterion .......... 51
9.1.3 Rock Mass Deformation ..................................................................................... 53
10 Worked Examples ............................................................................................................. 54
10.1 WORKED EXAMPLE 1: Point Load Index Test for UCS Determination ...................... 54
10.2 WORKED EXAMPLE 2: UCS from Schmidt Hammer Rebound Value ........................ 58
10.3 WORKED EXAMPLE 3: Determination of τp, τr, φp and φr of a Planar Discontinuity
from Direct Shear Test Results ............................................................................................. 60
10.4 WORKED EXAMPLE 4: Shear Strength of Clay Filled Discontinuity ........................... 63
10.5 WORKED EXAMPLE 5: Barton‐Choubey Shear Strength ........................................... 64
10.6 WORKED EXAMPLE 6: UCS (σci), Mohr‐Coulomb φ and c and from Triaxial Test Data
Conducted on Isotropic Rock Core ....................................................................................... 70
10.7 WORKED EXAMPLE 7: UCS (σci), mi, Tensile Strength (σti) and Hoek‐Brown &
Theoretical Failure Envelope ................................................................................................ 74
10.8 WORKED EXAMPLE 8: Hoek‐Brown Parameters s, a, mm, UCS and Uniaxial Tensile
Strength of Rock Mass .......................................................................................................... 78
10.9 WORKED EXAMPLE 9: Rock Mass Hoek‐Brown Equivalent Mohr‐Coulomb Friction
and Cohesion Parameters .................................................................................................... 80
10.10 WORKED EXAMPLE 10: Young’s Modulus of the Rock Mass for Deformation
Evaluation Purposes ............................................................................................................. 84
11 Self‐Assessment Activities ................................................................................................ 87
12 Further Reading ................................................................................................................ 87
13 Informative Website Links ................................................................................................ 87
14 References ........................................................................................................................ 88
Figure 28: Typical GSI Ranges for Jointed Rock (Marinos & Hoek, 2000) ............................................. 50
Figure 29: Point Load Index Size Correction Factor “F” (ASTM Standard D 5731, 1995) ..................... 56
Figure 30: Generalized Value of “C” (ASTM Standard D 5731, 1995) ................................................... 56
Figure 31: Schmidt Hammer Rebound Value to UCS (Deere & Miller, 1966) ...................................... 59
Figure 32: Shear Stress against Shear Displacement Plot (Worked Example) ...................................... 61
Figure 33: Shear Stress Shear against Normal Stress Plot (Worked Example) ..................................... 62
Figure 34: Discontinuity Plane Shear Strength – Barton‐Choubey Scenarios ....................................... 64
Figure 35: Variation in Shear Strength with Variation in Normal Stress .............................................. 66
Figure 36: Variation in Shear Strength with Variation in JRC ................................................................ 66
Figure 37: Variation in Shear Strength with Variation in JCS ................................................................ 67
Figure 38: Variation in Shear Strength with Variation in Residual Friction of the Rock ....................... 67
Figure 39: Plot of (σ1‐σ3)2 against σ3 with calculation of intercept and UCS ...................................... 71
Figure 40: Sigma 1 (Major Principal Stress σ1) vs Sigma 3 (Minor Principal Stress σ3) with Linear Fit . 72
Figure 41: hear Failure Occurs Along a Plane at Angle 𝜽 𝟒𝟓° 𝝋′𝟐 to the Major Principal Plane . 72
Figure 42: Mohr Circle and Cohesion (Shear Stress vs Normal Stress) (Franki, 2008, p. 53) ................ 73
Figure 43: Current Focus on Intact Rock Properties vs our Goal .......................................................... 74
Figure 44: Determination of mi and σci from Triaxial Data ................................................................... 75
Figure 45: Hoek‐Brown Failure Envelope ‐ Theoretical vs Lab Results ................................................. 77
Figure 46: Effect of GSI on UCS of the Rock Mass ................................................................................. 79
Figure 47: Current Focus on Intact Rock Properties vs our Goal .......................................................... 80
Figure 48: Brittle to Ductile Formation – Mogi’s σ1=3.4σ3 Limit (Hoek & Brown, 2018) ...................... 81
Figure 49: RockLab Software Package Analysis Export for Comparison Purposes ............................... 86
Figure 50: RockLab Software Package σ1‐σ3 and τ‐σn Plots for Illustrative Purposes .......................... 87
Basic Geology
(Learning Unit 1)
Basic Fundamentals
Geological and Introduction to
Engineering Projects
Structural Geology
(Learning Unit 11)
(Learning Unit 2)
Civil
Engineering
Design
Slopes Construction Geomorphology
(Learning Unit 9) (Learning Unit 4)
The study of the engineering properties of rock is termed “rock mechanics”. Rock mechanics
includes the study of the theoretical and applied science of the properties and mechanical
behaviour or rock and rock masses, in response to the forces acting on them.
The aim of studying the engineering properties of rock is to understand the behaviour of rocks
in different environments, and to predict how they will react or behave in response to stress
changes resulting from engineering works on or in the rock.
The properties of rock can be divided into two categories:
1
2
3
A number of worked examples are provided as part of this learning unit. You are not expected
to know the formulae by heart, but you should understand the concepts and be able to do
basic calculations, if formulae are provided.
4
Figure 2: The Importance of Scale – Intact Rock, Discontinuity and Rock Mass (Hoek,
2019)
The scale will determine the behaviour and the criteria we use to predict the strength and
behaviour of the rock or rock mass.
Examples:
1) Focus on the bench area (on a scale where jointing is limited) in Figure 2. If we want
to analyse the stability of this bench, that will most probably be a kinematical problem
(if the joints are spaced in such a manner that movement can occur). Thus, do we have
blocks or wedges, that are formed by discontinuities, that may slide into the pit? If so,
then the stability will depend on, and be dictated by, the shear strength of these
discontinuities.
2) If we look at the same rock, however, with the focus on the global (overall) stability of
the slope, the heavily jointed nature on this scale will result in a mass movement rather
than movement along a single joint set or even three joint sets. The strength and
behaviour of the rock mass will be different from that of the intact rock, or the
discontinuity-controlled scenario (discussed above). The material will undergo soil-like
(mass) behaviour in this scenario.
The different behaviours of the same rock at different scales (bench vs. overall slope) are
depicted in Figure 3. In the scenario depicted in Figure 3 a (scale of the bench with 1–3 joint
sets), the behaviour is dictated by the shear strength of the discontinuity. In the scenario
depicted in Figure 3 b, the behaviour is dictated by the shear strength of the rock mass
(different from that of the intact rock).
Different strength criteria are used to evaluate the different classes. That will become
apparent in this unit.
When dealing with a rock-related problem, consider the scale of the problem, as this will
dictate what you need to obtain regarding parameters/properties during your geotechnical
investigation, and which criteria will be used to analyse the problem.
5
Figure 3: Scale-effect – Discontinuity or Rock Mass Strength
6
Figure 4: Common Conditions in Analysis of a Rock Slope (Mah & Wyllie, 2004, pp. 75-
77)
7
In preliminary assessments, the physical (or index) properties are generally used as a first
estimation of the expected mechanical properties (strength and deformation parameters). The
strength parameters for intact rock are the uniaxial compression and tensile strength. The
strength parameters for discontinuities or rock masses are cohesion and friction. The
deformation parameters are the stiffness modulus (Young’s modulus or E-modulus) and
Poisson’s ratio.
Depending on the scale of the problem (intact rock vs. mass and stress conditions), some
parameters have little effect on the prediction of the expected behaviour, and it may be
estimated by means of visual assessment and basic field correlations and tests. This may be,
for example, the UCS of the intact rock in a shallow slope assessment, where the failure is
dictated by the shear strength of a discontinuity (or infill of that discontinuity).
Here, the strength of the intact rock is not significant and a good estimation of the UCS will be
acceptable for this specific scenario (strong rock in a low-stress regime). On the other side of
the spectrum, you may be faced with an underground excavation project in soft rock under
relatively high-stress conditions. Here, laboratory and in-situ testing and a determination of
strength and deformation will be critical for accurately predicting the behaviour of the rock or
rock mass (rock burst in hard rock or squeeze in relatively soft rock).
8
Your understanding of the project, the level of input for a specific evaluation phase and the
design philosophy are critical in a properly planned and executed geotechnical investigation,
analysis and design. This topic will be covered in the site investigation unit.
Some of the important physical properties used in determining the strength and deformation
characteristics of rock are discussed in the sections which follow.
Where:
γ is the unit weight of the rock (kN/m3)
ρ is the density (kg/m3)
9
𝑄 𝑘 ∙ 𝑖 ∙ 𝐴 Equation 6.4
Where:
Q is the flow per unit area
k is the hydraulic conductivity (m/s, cm/s, m/day)
i is the gradient of the potential head
For the purposes of this unit, the flow in intact rock can be considered to follow Darcy’s Law:
∆ℎ Equation 6.5
𝑞 𝑘∙ ∙ 𝐴
∆𝑥
Where:
qx is the volume of flow in the x direction
k is the hydraulic conductivity (m/s, cm/s, m/day)
∆h is the difference in hydraulic head
∆x is the volume of flow over time (volume/time)
A is the section perpendicular to the x direction (flow direction)
𝑘 𝐾∙ Equation 6.6
µ
Where:
K is the intrinsic permeability (only dependent on the characteristics of the material)
γw is the unit weight of the fluid (or water)
µ is the kinematic viscosity of the fluid (or water)
10
Where:
A is area (m2)
The UCS test is conducted in the laboratory on a cylindrical specimen of rock, generally
obtained from a rotary core sample. The UCS is simply the stress as determined by:
11
The UCS can also be estimated or correlated from indirect tests. Referred to as index test
methods, these tests include:
Point load tests;
Schmidt hammer rebound test;
Ultrasonic test (wave propagation); and
Field indices by means of visual observation and professional judgement.
12
σ1
∆l
li lf
σ1
13
What will the stress-strain curve for a strong and brittle rock look like? Do you expect post-
peak stress strain?
14
If the core diameter tested is approximately 50 mm-equivalent size, say 40–60 mm, F can be
determined by:
.
𝐷 Equation 6.13
𝐹
50
If there is a large correction to be made, say for example the core tested is PQ-size, thus
around 85 mm, F can be calculated by:
𝐷 Equation 6.14
𝐹
50
Generally, NWD4 and TNW-size core barrels are most often used in smaller civil projects in
this country, where the core diameter is in the region of 52.3 to ~53.0 mm.
The uniaxial compressive strength is then estimated by:
15
Thus a “C” value of 42. This correlation factor is commonly used in South Africa. Note,
however, that this “C” value should be determined for your site-specific conditions by means
of UCS tests, and correlation of the PLI with these UCS tests.
Numerous other correlations are published. For a more detailed description on the PLI test
and testing of non-cylindrical samples, refer to the ASTM Standard D 5731 (1995) test method,
which is freely available from the internet.
Figure 7: Point Load Index Test (Is)
16
Figure 8: Point Load Index Correlation with UCS after Beniawski (1973) Source
(Franki, 2008, p. 57)
17
Figure 9: Schmidt Hammer Rebound UCS Correlation Chart (Hoek, 2019)
18
Note: Soil is not defined by UCS of less than 1 MPa or even 0.75 MPa, as stipulated in certain
codes and standards. Keep in mind the application of the description. Numerous engineers
may define a rock with a UCS of 300 kPa as soil. The failure mechanism may, however, be
that of a rock, and the Hoek-Brown or Barton-Choubey criterion may be the applicable
evaluation approach. More refined ISRM ranges are recommended. You can always simplify
and adapt a detailed description, but you cannot elaborate on a basic description with some
UCS limit trying to force a boundary between a rock and a soil.
19
As with the UCS test, it is the maximum stress the material can take, under tension. The stress
state and shear stress vs. normal stress plot is presented in Figure 10.
Figure 10: Stress Conditions and Shear Stress – Normal Stress Plot
20
Figure 11: a) Indirect Tensile Test (Brazil) and b) Point Load Test
21
7 Discontinuities
Depending on the scale and geological-structural conditions, the strength and deformation of
rock may be dictated by either the properties of the 1) intact rock itself, 2) shear properties
along discontinuities, or 3) shear strength and deformation characteristics of the rock mass.
This section focuses on discontinuities. Our scale of interest is depicted in Figure 12.
Figure 12: Intact Rock to Rock Mass – Sample Size – Scale Discontinuities (Hoek,
2019)
22
the weaker and more deformable the rock becomes. The wider the spacing, generally
the stronger the rock and the less deformable it becomes.
Persistence (continuity): Persistence is the fraction of area that is discontinuous. Is
the discontinuity on the micro-scale, only a couple of millimetres, metres or does it
extend across and beyond the area of interest?
Roughness and waviness: The rougher the discontinuity surfaces, the more
resistance against sliding and the stronger the discontinuity. The higher the amplitude
of the wave, the higher the resistance against sliding. It is obviously important to
determine whether movement will be along or across the waviness.
Aperture: This is the perpendicular distance between the walls of an open
discontinuity (if not filled).
Discontinuity fill: The open discontinuity may be filled with mineralisation, sand, clay
or other fill materials. The discontinuity fill will have a significant impact on the shear
strength of the discontinuity and overall rock mass behaviour. Clay fill may result in the
low shear strength of a joint surface, and mineralisation may result in the higher shear
strength of a joint surface.
Discontinuity water: Water reduces shear strength and thus has an impact on the
shear strength of a discontinuity or the overall behaviour of the rock mass. Water has
an impact on the effective stress conditions and physical properties of infill. The impact
of water is covered at a later stage.
Discontinuity wall compressive strength: A low wall strength indicates weathering
and alteration. The lower the discontinuity wall strength, the lower the shear strength
of the discontinuity and overall rock mass strength.
Typical parameters describing a rock mass are presented in Figure 13.
24
Figure 13: Parameters Describing the Rock Mass (Mah & Wyllie, 2004, p. 55)
25
Figure 14: Planar Discontinuity: a) Shear Stress and Displacement Plots and b) Peak
and Residual Strength (de Vallejo & Ferrer, 2011)
The peak shear strength is where maximum shear stress is obtained on the stress-
displacement plot. When maximum shear stress is exceeded, the discontinuity will still have
26
strength, but it requires less stress (force) to move. The shear stress will taper down and reach
a “constant”. This strength is the residual shear strength.
Exercise – Actual (rough discontinuity surface): Put your phone flat on the table where
you are sitting (hopefully your phone has a rubber, silicon or leather cover). Start pushing
against the side, very slowly, at a low angle. Do not move it at first, but note the increase in
pressure (or rather stress) on your fingertip (before the phone starts to slide). Now note what
happens once it starts to slide. Did you note the sudden pressure (stress) release on your
finger, and that you can push it with less effort along the table once it is on the move? If you
plot the stress vs. shear displacement you did on your phone exercise, you will obtain the
same figure as presented in Figure 14 a, with reference to the upper plot with the peak shear
strength. Its shear stress – normal stress plot, before the phone moved – will be presented by
the upper plot (peak strength) in Figure 14 b. Once mobilised, the shear stress – normal stress
plot – will be presented by the lower plot (residual strength) in Figure 14 b.
Exercise – Visualisation (smooth discontinuity surface): Now, visualise two pieces of flat
glass with some water separating the glass. Visualise doing the same exercise. You will
probably not note the sudden release in pressure (stress) on your fingertip. You will merely
note a slight increase in pressure (stress) on your fingertip, up to a point where it starts to
slide, and that pressure (stress) sensation will be maintained as long as you slide the glass
along. The τ-δ plot of this scenario will result in the same plot as presented in Figure 14 a,
with reference to the lower plot. The entire τ-σn plot will be presented by the lower plot (residual
strength) in Figure 14 b.
Now, is the peak strength plot of the sliding phone scenario representative or “true”? Refer to
the peak strength in Figure 14 b. Slowly lift your phone from the table. Does it stick to the
table? Do you note a release in effort once it starts to lift, or not? Hopefully not! Thus, there is
no cohesion between the phone and table surface. Figure 14 b, however, indicates a cohesion
intercept that is significant on this plot. So, what better criteria can we come up with for rough
surfaces? We introduce a bi-linear failure criterion, combining the peak and residual plots
in Figure 14 b into one plot, as presented in Figure 15 b.
27
Figure 15: Bi-linear Failure Criterion for Rough Discontinuity Surfaces (de Vallejo &
Ferrer, 2011)
It should be clear at this stage that the peak shear strength (cell-phone-on-table exercise) is
not due to cohesion, but rather a frictional component that needs to be overcome (override)
due to the roughness of the surfaces. This overriding is, in simple terms, dilation. You have
to overcome the roughness to get the phone sliding.
Now, place your phone of a coarse brick or concrete slab and repeat the exercise. What do
you reckon will happen with the angle of the plots of the dilation and shearing phases? What
will happen with the peak on the τ-µ plot? Answer: The peak will increase, and both the angles
on plot b) will increase. Note that for this phone exercise, we have considered low applied
normal stresses. Thus, we need to dilate (override) the roughness component. If we increase
the normal stress (put a number of bricks on your phone), will it dilate? Or will the roughness
of the brick simply shear through the leather cover? Do you think you will see the initial
envelope (dilation phase) followed by the less steep slope (shearing phase)? We may get to
this one at a later stage. The importance of the stress regime in a project should be clear by
now. The same material in environments with different stress regimes will behave differently.
An example is shallow excavation with low induced stresses vs. a deep tunnel with high
induced stresses.
28
The basic friction angle is the free sliding angle measured on two smooth rock surfaces in a
tilt test, such as two smooth saw-cut samples or smooth core samples (illustrated in Figure
16 b).
The residual friction angle is the friction angle after significant shearing. It is thus the friction
angle calculated at the near constant applied shear stress on the τ-ε plot. It can also be
determined using the direct box shear test (see worked example).
Figure 16: Peak Friction vs. Basic Friction Angle of Rock (Barton, 2017)
29
Where:
τ is shear strength on the failure plane
c is the cohesion
σ is the normal stress
φ is the friction angle
This equation has two components: cohesive and frictional.
Cohesive Frictional
component component
(c) (σ and φ)
𝜏 𝑐 𝜎 · 𝑡𝑎𝑛𝜑
Cohesion is “constant”, while friction is proportional to the normal stress (σ).
For a smooth discontinuity with no cohesion (thus no cementation or mineralisation with actual
cohesive bonds), the shear strength can be expressed as:
Where:
φr is the residual angle of friction (notation “r” for residual)
σn is the normal stress (notation “n” for normal)
With no cohesion, the cohesion component falls away. Cohesion c=0.
The superscripts, ′ , merely indicate effective cohesion, effective normal stress and effective
angle of friction.
30
In this initial phase, before reaching peak shear strength, we make use of the peak friction
angle φp. The peak friction is the basic friction angle plus the angle of the irregularity in
relation to the discontinuity plane.
The peak friction angle is defined by:
𝜑 𝜑 𝑖 Equation 7.4
Where:
φp is the peak friction angle
φb is the basic friction angle (refer to the definition of a basic friction angle)
i is the angle of the irregularity (or average for irregularities) in relation to the discontinuity
plane.
The shear strength for a rough discontinuity thus becomes:
Fundamentally, cohesion is broken at small strains (if you have any cohesion at all!), while
friction is mobilised and then remains. Due to dilation (remember the initial overriding of rough
discontinuity surfaces?), discontinuities have very high friction angles at low stress, from
where they decrease after a certain amount of movement (refer back to Figure 15). The
mechanism of shear failure and deformation is thus ignored by adopting this approach. Barton
advocates the application of a criterion which is fundamentally more acceptable where the
cohesion is first used, the cohesion is lost, and then the frictional component is utilised (as
opposed to the c + φ, he recommends c and then φ). This approach is compatible with the
actual mechanism of shear failure and deformation.
Barton performed numerous shear tests across a range of materials. He realised that the peak
shear strength of rough and clearly unweathered tension fractures could be described by a
relationship involving the uniaxial compression strength of the wall and that it is related to the
applied normal stress:
𝜎
𝜏 𝜎 ∙ tan 20𝑙𝑜𝑔 30° Equation 7.7
𝜎
The equation was derived from the link between and and the peak dilation angles, from
where the following relationship was developed:
31
32
Do not consider cohesion. Joint roughness and wall compression strength incorporated in the formulation. This
Barton not in favour of “cohesion + friction” but approach incorporates stress condition and higher initial friction angles at low normal
rather “cohesion then friction” that is fundamentally loads (Dilation effect). See what happens if you increase σn or increase and decrease
correct. Very valid argument! roughness and/or “strength” of the wall (thus ability of roughness to restrict
movement).
High stress levels relative to rock wall strength (JCS), this term becomes zero. Thus
simulating asperities shearing off.
Low stress levels relative to rock wall strength (JCS), this term allows for roughness to
add to residual friction. Thus simulating overriding of asperities. Dilation effect.
Cohesive Frictional
component component
33
Figure 17: Intact Rock to Rock Mass – Sample Size – Scale Intact Rock (Hoek, 2019)
34
The strength and deformation of an isotropic rock will be similar in all directions, as the
material is homogeneous. The strength and deformation characteristics of an anisotropic
rock will be different in different directions. This concept is illustrated in Figure 18.
The focus will be on isotropic rock and its strength criteria (Figure 18 a and b). Only the
concepts of variation in strength and deformation in anisotropic rock will be introduced (Figure
18 c and d).
35
Figure 19: Stages in Crack Formation for a Rock Specimen Under Uniaxial
Compression
36
37
𝜏 𝑐 𝜎 · 𝑡𝑎𝑛𝜑
σ1 σ 𝑓 σ
τ
σ1
θ c 2θ
θ
σ3,2 σ3,2
σ3 σn σ1
ψ
θ σc
σ1
σ3
Where:
τ is shear strength on an arbitrary or potential failure plane
c is the cohesion of the intact rock
σn is the normal stress
φ is the friction angle of the intact rock
You can also express this criterion as a function of the principal stresses σ1 and σ3:
2𝑐 𝜎 𝑠𝑖𝑛2𝜃 𝑡𝑎𝑛𝜑 1 𝑐𝑜𝑠2𝜃 Equation 8.2
𝜎
𝑠𝑖𝑛2𝜃 𝑡𝑎𝑛𝜑 1 𝑐𝑜𝑠2𝜃
Where:
σ1 is the major principal stress
σ3 is the minor principal stress
c is the cohesion of the intact rock
θ is the critical plane
𝜑 Equation 8.3
𝜃 45
2
σc can be expressed as:
38
This linear failure criterion generally overestimates the tensile strength of rock. A tension cut-
off is thus applied when using this criterion (illustrated in Figure 21).
Tension cut-off
applied in order not
to overestimate the
τ Actual failure envelope
σn
Linear Mohr-Coulomb
Failure Envelope
Figure 21: Tension Cut-off Concept for Linear Mohr-Coulomb Failure Criterion
39
Hoek and Brown proposed a non-linear empirical failure criterion, which was introduced in
1980 and is for isotropic intact rock strength under triaxial conditions (thus strength at different
stress states can be evaluated). Hoek and Brown mainly used the crack propagation criterion
proposed by Griffith in 1921 and 1924, which was derived from brittle materials subjected to
tensile stresses (i.e., concrete and glass). Hoek and Brown modified this criterion by fitting
parameters (thus observational manipulation – fitting of curves) and extending them across
the compressional stress range. These fitting parameters thus have no fundamental
relationship to the physical characteristics of the rock.
In geotechnical engineering, failure generally occurs in shear, so we present the failure
criterion in terms of shear and normal stresses on the failure plane (as per the Mohr-Coulomb
failure criterion presented in Equation 8.1). In rock mechanics it is common to present failure
criteria in terms of the principal stresses (thus σ1 and σ3). The Hoek-Brown criterion is
expressed in terms of the major principal stresses (σ1 and σ3). Later in the unit we will convert
the principal stresses to shear and normal stresses, and derive equivalent cohesion and
frictional parameters in order to evaluate rock mass behaviour with the Mohr-Coulomb
criterion.
The Hoek-Brown criterion for intact rock can be defined by:
.
𝜎 Equation 8.8
𝜎 𝜎 𝜎 𝑚 ∙ 𝑠
𝜎
Where:
σci is the uniaxial compressive strength (UCS) of the “i” – intact rock.
The constants mi, s and a define the shape of the curvilinear failure plot. These are curve-
fitting parameters derived from visual fitting. It is important to note, upon closer inspection of
the “mi”, “s” and “a” curve-fitting parameters, that:
The constant “mi” is analogous to the frictional strength of the rock (rock type, angle of
inter-particle or inter-block friction and degree of particle interlock with the rock mass).
Parameter “s” is associated with the degree of fracturing, thus analogous with the
cohesive strength of the rock mass. Isotropic rock s = 1 thus simulating cohesion;
highly fractured rock s = 0, thus simulating no cohesion.
We will not go into detail on these curve-fitting parameters. It is important to note that although
these curve-fitting parameters are based on pure observational fitment to force a workable
criterion (bending the failure envelope), the cohesion, friction and dilatancy effects are actually
hidden in these parameters.
40
The value of mi can be determined from triaxial test data. If such data are not available,
preliminary estimates can be obtained from Table 3.
Table 3: Estimates of the Constant mi for Intact Rock (If Triaxial Data not Available)
Note 1 (Conglomerates and Breccias): Conglomerates and breccias may present a wide range of mi values
depending on the nature of the cementing material and the degree of cementation, so they may range from
values similar to sandstone to values used for fine‐grained sediments.
41
Note 2 (Foliated): The values indicated are for intact rock specimens tested normal to bedding or foliation. The
value of mi will be significantly different if failure occurs along a weakness plane.
𝜎 𝜎 𝑚𝜎 𝜎 𝜎 Equation 8.9
𝜎 𝑚 4𝑠 𝑚 Equation 8.10
𝜎 2
𝑚 4𝑠 𝑚
𝜎 𝜎 ⎛ ⎞ Equation 8.11
2
⎝ ⎠
As with the Mohr-Coulomb criterion, a tension cut-off was introduced (see Figure 21), and the
relationship between the compressive and tensile strength ratio and the Hoek-Brown
parameter mi was proposed as:
𝜎
0.81𝑚 7 Equation 8.12
|𝜎 |
Where:
|𝜎 | is the absolute value of the uniaxial tensile strength (determined by direct tensile test).
42
a) c)
Strength
b) d)
43
a) c)
Less Stiff
b) d)
Note the anisotropic behaviour and the impact it will have on the strength and deformation
characteristics of the rock or rock mass. Strength anisotropy can be evaluated through
systematic laboratory testing of specimens drilled in different directions from orientated block
samples. There are behavioural criteria for anisotropic rock, but those will not be covered in
this learning unit. For the present purposes, only take note of anisotropy and the effect on
strength and stiffness.
9 Rock, As a Mass
Depending on the scale and geological-structural conditions, the strength and deformation of
rock may be dictated by 1) the properties of intact rock itself, 2) shear properties along
44
discontinuities or 3) the shear strength and deformation characteristics of the rock mass. This
section focuses on the “rock mass”. Our scale of interest is depicted in Figure 24.
The processes followed in the Hoek-Brown criterion are presented in Figure 25 and Figure
26.
Figure 24: Intact Rock to Rock Mass – Sample Size – Scale Rock Mass (Hoek, 2019)
45
Figure 25: Flow Chart – Application of Hoek-Brown Criterion and GSI to an Excavation
Design (Hoek & Brown, 2018)
46
Figure 26: Data Entry Stream for Using the Hoek-Brown System for Estimating Rock
Mass Parameter for Numerical Analysis (Hoek, et al., 2013)
47
48
1 1 Equation 9.4
𝑎 𝑒 𝑒
2 6
Where:
mm is a reduced value of the material constant mi (determined from triaxial test results,
covered in intact rock section), but can also be estimated from Table 3. The rock mass is
weaker than the intact rock.
D is a disturbance factor which depends on the degree of disturbance to which the rock mass
has been subjected by blast damage/excavation and stress relaxation (0 for undisturbed; to
1 for highly disturbed).
GSI is the Geological Strength Index (qualitative classification system). Typical ranges can
be obtained from Figure 28.
49
Figure 28: Typical GSI Ranges for Jointed Rock (Marinos & Hoek, 2000)
50
The uniaxial compressive strength of the rock mass “m” (or broken rock “b”) can be
determined in terms of σci and parameter “s” and “a” as:
𝜎 𝜎 ∙𝑠 Equation 9.5
Where:
σcm is the uniaxial compressive strength of the rock mass
σci is the uniaxial compressive strength of the intact rock
s and a are the material parameters.
Set s = 1 and a = 0.5, as for intact rock, and see whether the relationship holds: 1 exp0.5 = 1,
thus σc = σci, so yes.
Marinos and Hoek (2001) also propose an empirical equation for σcm in terms of mi, σci and
GSI as:
𝜎 𝜎 ∙ 0.0034 ∙ 𝑚 .
∙ 1.029 0.025𝑒 . Equation 9.6
𝜎 𝑚 4𝑠 𝑚 Equation 9.7
𝜎 2
Hoek (1983) notes that for brittle materials, the uniaxial tensile strength is equal to the biaxial
tensile strength. Thus, substituting 𝜎 𝜎 𝜎 in Equation 9.1 you obtain the tensile
strength of the rock mass by:
𝑠𝜎
𝜎 Equation 9.8
𝑚
51
⎛ ⎞
𝜎 , 𝜎 ,
⎜ 1 2∙𝑎 ∙𝑠 1 𝑎 ∙𝑚 ∙ 𝑠 𝑚 ∙ ⎟
𝜎 𝜎
𝑐 𝜎 ⎜ ⎟ Equation 9.10
⎜ 𝜎 ⎟
,
⎜ 6∙𝑎∙𝑚 𝑠 𝑚 ∙ ⎟
𝜎
1 𝑎 2 𝑎 ∙ 1
⎝ 𝑎 𝑎 2 𝑎 ⎠
⎡ 𝜎 , ⎤
⎢ 6∙𝑎∙𝑚 ∙ 𝑠 𝑚 ∙ ⎥
𝜎 Equation 9.11
𝜑 sin ⎢ ⎥
⎢2 ∙ 1 𝜎 , ⎥
𝑎 ∙ 2 𝑎 6∙𝑎∙𝑚 𝑠 𝑚 ∙
⎣ 𝜎 ⎦
52
Ei is the intact rock deformation modulus (Young’s modulus). This can be determined by
laboratory testing, or estimated (reduction values as proposed by Deere, 1968) depending
on the level of accuracy required.
D is a disturbance factor which depends on the degree of disturbance to which the rock mass
has been subjected by blast damage/excavation and stress relaxation (0 for undisturbed, 1
for highly disturbed).
GSI is the Geological Strength Index (qualitative classification system).
53
Hoek and Diederichs (2006) also propose an alternative equation for estimating rock mass
modulus when Ei is not available:
𝐷
1
𝐸 10 2 Equation 9.16
75 25𝐷 𝐺𝑆𝐼
1 exp
11
Thus, you can now estimate the stiffness of the rock mass by Equation 9.15 or Equation
9.16, for use in deformation calculations.
10 Worked Examples
These worked examples illustrate the practical application of some of the aspects/concepts
covered in this learning unit. They serve as general guidance. Before applying any criteria,
ensure that you have done extensive research on the topic and know the
limitations/restrictions and intended range of conditions for each criterion. The worked
examples provide a sound foundation on which to build. If you are unsure, seek independent
input.
54
Calculation steps:
1. Convert failure load, recorded as kN, to MN.
2. Convert core diameter, mm, to m.
55
F Size Correction Factor (50 mm = 1 ‐‐> No Correction)
1.8
1.6
F Size Correction Factor
1.4
1.2
1.0
0.8
y = 2E‐07x 3 ‐ 9E‐05x 2 + 0.0182x + 0.2905
0.6
R² = 0.9999
0.4
0.2
0 25 50 75 100 125 150 175
De (Equivalent) Core Diameter (mm)
Figure 29: Point Load Index Size Correction Factor “F” (ASTM Standard D 5731, 1995)
Value of "C" (Generalized)
25
24
Correlation Factor "C"
23
22
21
20
19
y = 0.1848x + 13.676
18
R² = 0.9934
17
15 20 25 30 35 40 45 50 55 60 65
Diameter of the Core (mm)
56
8 9 10 11 12
Applied
Is(50) Is (50)
Test No. Test No Correlation UCS (σci)
Corrected Corrected
Factor
2 7.605 2
3 7.911 3
8 7.945 8 7.945 23 182.74
11 8.074 11 8.074 23 185.71
4 8.397 4 8.397 23 193.14
7 8.453 7 8.453 23 194.43
12 8.466 12 8.466 23 194.72
9 8.608 9 8.608 23 197.99
1 8.746 1 8.746 23 201.16
10 8.746 10 8.746 23 201.16
6 9.276 6
5 9.948 5
Min Is(50): 7.95 182.74 MPa
Max Is(50): 8.75 201.16 MPa
Average Is(50): 8.43 193.88 MPa
Stdev Is (50): 0.29 6.71 MPa
CoV: 3.46% 3.46%
For guidance on UCS correlations on non-cylindrical samples and tests conducted on rock
with an anisotropic nature, refer to ASTM designation: D 5731-95 test method. It is freely
available from the internet.
57
58
Figure 31: Schmidt Hammer Rebound Value to UCS (Deere & Miller, 1966)
59
60
Table 4: Direct Shear Box Test Results Obtained from the Laboratory
Shear Stress (kN/m^2) Shear Displacement (mm)
0 0
159 0.05
200 1.19
241 3.61
228 4.5
214 8.51
207 9.4
200 11.61
193 12.6
179 17.09
179 19.81
Direct Shear Tests
260
Peak τp
240
Shear Stress (kN/m^2)
220
200
Residual τr
180
160
140
0 2 4 6 8 10 12 14 16 18 20
Shear Displacement (mm)
Figure 32: Shear Stress against Shear Displacement Plot (Worked Example)
61
Direct Shear Tests
300
Peak Shear Stress
Residual Shear Stress
250
200; 241
Shear Stress (kN/m^2)
200
200; 180
150
100
50 Peak φp
Residual φr
0 0; 0
0 50 100 150 200 250 300
Normal Stress (σn, kN/m^2)
Figure 33: Shear Stress Shear against Normal Stress Plot (Worked Example)
62
The shear strength along a smooth discontinuity can be defined by the Mohr-Coulomb
failure criterion:
Where:
τ is shear strength on the failure plane
c is the cohesion
𝜑 is the angle of internal friction of the soil
You now have the shear strength of the discontinuity, with no pore water pressure. You need
to determine the shear stress caused by the block of rock on this discontinuity, and ensure
that the shear strength exceeds the shear stress with some safety factor.
You will have the opportunity to go into more detail in the relevant learning units still to come.
63
𝜎 𝑁𝑜𝑟𝑚𝑎𝑙 𝑆𝑡𝑟𝑒𝑠𝑠
Rock
Planar Discontinuity
𝑆ℎ𝑒𝑎𝑟 𝑆𝑡𝑟𝑒𝑠𝑠
𝑆ℎ𝑒𝑎𝑟 𝑆𝑡𝑟𝑒𝑠𝑠
Rock
The purpose of this worked example is to illustrate the Barton-Choubey criterion in terms of
variation of shear strength, with a variation of:
1. Normal stress σn
2. Joint roughness coefficient (JRC)
3. Joint wall compression strength of the discontinuity JCS
4. Residual friction angle of the rock φr
The shear strength may later be used in a practical slope stability problem.
The shear strength criterion (Barton-Choubey):
𝐽𝐶𝑆
𝜏 𝜎 ∙ 𝑡𝑎𝑛 𝐽𝑅𝐶 ∙ 𝑙𝑜𝑔 𝜑
𝜎
Take note of the ranges of the different parameters.
1. JRC ranges between 0 (for smooth discontinuities) and 20 (for rough and stepped
discontinuities).
2. JCS is the UCS for unweathered walls (JCS = UCS). For altered walls the JCS may
be less than 25% of the UCS (JCS is determined via testing, but let’s just consider
percentages for the purposes of this worked example).
3. Residual friction for unaltered walls is generally 25–35° and as low as 15–12°, if
altered.
64
65
Variation in τ and Frictional Component i with Change in σn
90.00 35
Roughness (Frictional Component i) degrees
σn (kN/m^2)
80.00 Friction Component (i) 30
70.00
25
Shear Strength (τ)
60.00
50.00 20
40.00 15
30.00
10
20.00
5
10.00
0.00 0
0 20 40 60 80 100 120
Normal Stress (σn) (kN/m^2)
Note the increase in shear strength of the discontinuity with an increase in normal stress, as
expected. Importantly, note the decrease in the friction component “i” with an increase in
normal stress. Barton thus incorporates the dilation effect (sharing of asperities with an
increase in normal stress). Thus, if normal stress increases, the asperities start to shear off.
Variation in τ and Frictional Component i with Change in JRC
45.00 40
Roughness (Frictional Component i) degrees
JRC
40.00 Friction Component (i) 35
30
35.00
Shear Strength (τ)
25
30.00
20
25.00
15
20.00
10
15.00 5
10.00 0
0 2 4 6 8 10 12 14 16 18 20
JRC
Figure 36: Variation in Shear Strength with Variation in JRC
Note that the frictional component increases with an increase in JRC, as expected. Thus, also
an increase in the shear strength of the joint. This is at constant normal stress.
66
Variation in τ and Frictional Component i with Change in JCS
18.20 13
Roughness (Frictional Component i) degrees
JCS (kN/m^2
18.00
Friction Component (i) 12
17.80
17.60 11
Shear Strength (τ)
17.40 10
17.20
9
17.00
16.80 8
16.60 7
16.40
6
16.20
16.00 5
0 1000 2000 3000 4000 5000 6000
JCS (kN/m^2)
Note the increase in the frictional component “i” with an increase in the wall strength JCS
component. Thus, there is also an increase in shear strength of the joint/discontinuity.
Variation in τ and Total Friction with Change in φr
45.00
Friction Component (φr)
40.00 Frictional Total (φr + i)
35.00
Shear Strength (τ)
30.00
25.00
20.00
15.00
10.00
5.00
12 16 20 24 28 32 36
Friction Angle (degrees)
Figure 38: Variation in Shear Strength with Variation in Residual Friction of the Rock
Note the increase in shear strength and the overall frictional component, with an increase in
the residual friction angle of the rock φr.
67
Variation of Normal Stress (σn): SCENARIO 1
Calculation No 1 2 3 4 5 6 7 8 9 10
σn (kN/m^2) 1 2 4 6 10 20 40 60 80 100
JRC 10 10 10 10 10 10 10 10 10 10
JCS (kN/m^2 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
φr 30 30 30 30 30 30 30 30 30 30
JCS/σn 1000.00 500.00 250.00 166.67 100.00 50.00 25.00 16.67 12.50 10.00
Friction Component (i) 30.0 27.0 24.0 22.2 20.0 17.0 14.0 12.2 11.0 10.0
Radians (i) 0.52 0.47 0.42 0.39 0.35 0.30 0.24 0.21 0.19 0.17
Friction Component (φr) 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00
Radians (φ) 0.52 0.52 0.52 0.52 0.52 0.52 0.52 0.52 0.52 0.52
Frictional Total (φr + i) 60.00 56.99 53.98 52.22 50.00 46.99 43.98 42.22 40.97 40.00
Frictional Total (radians) 1.05 0.99 0.94 0.91 0.87 0.82 0.77 0.74 0.72 0.70
Shear Strengh of Discontinuity (τ) kN/m^2 1.73 3.08 5.50 7.74 11.92 21.44 38.60 54.44 69.47 83.91
Variation of JRC: SCENARIO 2
Calculation No 1 2 3 4 5 6 7 8 9 10
σn (kN/m^2) 20 20 20 20 20 20 20 20 20 20
JRC 0 2 4 6 8 10 12 14 17 20
JCS (kN/m^2 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
φr 30 30 30 30 30 30 30 30 30 30
JCS/σn 50.00 50.00 50.00 50.00 50.00 50.00 50.00 50.00 50.00 50.00
Friction Component (i) 0.0 3.4 6.8 10.2 13.6 17.0 20.4 23.8 28.9 34.0
Radians (i) 0.00 0.06 0.12 0.18 0.24 0.30 0.36 0.42 0.50 0.59
Friction Component (φr) 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00
Radians (φ) 0.52 0.52 0.52 0.52 0.52 0.52 0.52 0.52 0.52 0.52
Frictional Total (φr + i) 30.00 33.40 36.80 40.19 43.59 46.99 50.39 53.79 58.88 63.98
Frictional Total (radians) 0.52 0.58 0.64 0.70 0.76 0.82 0.88 0.94 1.03 1.12
Shear Strengh of Discontinuity (τ) kN/m^2 11.55 13.19 14.96 16.90 19.04 21.44 24.17 27.31 33.13 40.97
68
Variation of JCS (% of UCS): SCENARIO 3 100% 90% 80% 70% 60% 50% 40% 35% 30% 25%
Calculation No 1 2 3 4 5 6 7 8 9 10
σn (kN/m^2) 20 20 20 20 20 20 20 20 20 20
JRC 5 5 5 5 5 5 5 5 5 5
JCS (kN/m^2 5000 4500 4000 3500 3000 2500 2000 1750 1500 1250
φr 30 30 30 30 30 30 30 30 30 30
JCS/σn 250.00 225.00 200.00 175.00 150.00 125.00 100.00 87.50 75.00 62.50
Friction Component (i) 11.99 11.76 11.51 11.22 10.88 10.48 10.00 9.71 9.38 8.98
Radians (i) 0.21 0.21 0.20 0.20 0.19 0.18 0.17 0.17 0.16 0.16
Friction Component (φr) 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00
Radians (φ) 0.52 0.52 0.52 0.52 0.52 0.52 0.52 0.52 0.52 0.52
Frictional Total (φr + i) 41.99 41.76 41.51 41.22 40.88 40.48 40.00 39.71 39.38 38.98
Frictional Total (radians) 0.73 0.73 0.72 0.72 0.71 0.71 0.70 0.69 0.69 0.68
Shear Strengh of Discontinuity (τ) kN/m^2 18.00 17.86 17.70 17.52 17.31 17.07 16.78 16.61 16.41 16.18
Variation of φr: SCENARIO 4
Calculation No 1 2 3 4 5 6 7 8 9 10
σn (kN/m^2) 20 20 20 20 20 20 20 20 20 20
JRC 5 5 5 5 5 5 5 5 5 5
JCS (kN/m^2 100 100 100 100 100 100 100 100 100 100
φr 36 34 32 30 28 26 24 22 20 15
JCS/σn 5.00 5.00 5.00 5.00 5.00 5.00 5.00 5.00 5.00 5.00
Friction Component (i) 3.49 3.49 3.49 3.49 3.49 3.49 3.49 3.49 3.49 3.49
Radians (i) 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06
Friction Component (φr) 36.00 34.00 32.00 30.00 28.00 26.00 24.00 22.00 20.00 15.00
Radians (φ) 0.63 0.59 0.56 0.52 0.49 0.45 0.42 0.38 0.35 0.26
Frictional Total (φr + i) 39.49 37.49 35.49 33.49 31.49 29.49 27.49 25.49 23.49 18.49
Frictional Total (radians) 0.69 0.65 0.62 0.58 0.55 0.51 0.48 0.44 0.41 0.32
Shear Strengh of Discontinuity (τ) kN/m^2 16.48 15.34 14.26 13.24 12.25 11.31 10.41 9.54 8.69 6.69
69
σ1f (MPa) σ3 (MPa)
1 78.00 0.0
2 124.50 5.0
3 145.50 10.0
4 196.00 20.0
5 230.50 30.0
6 262.50 40.0
We want 1) the UCS and 2) to evaluate τf (shear strength at failure, thus ultimate shear
strength). From Equation 10.1 we know we require cohesion “c” and the friction angle “φ” to
calculate τf. The normal load σn is specific to the project.
70
Determination of UCS (σc) and Friction Angle (φ) from Triaxial Test Data
60000
𝑈𝐶𝑆 𝜎 𝐼𝑛𝑡𝑒𝑟𝑐𝑒𝑝𝑡
50000
(σ1‐σ3)^2 (Units in MPa)
40000
Slope Intercept
30000
20000
y = 1070.4x + 7835
R² = 0.9937
10000
Intercept where σ3=0, thus on y-axis
0
‐10.0 0.0 10.0 20.0 30.0 40.0 50.0
σ3 (Minor Principal Stress) ‐ Units in MPa
Figure 39: Plot of (σ1-σ3)2 Against σ3 With Calculation of Intercept and UCS
71
Determination of UCS (σc) and Friction Angle (φ) from Triaxial Test Data
300
𝜑
tan 45° 𝑆𝑙𝑜𝑝𝑒 𝑚
2
250 𝜑 2 ∙ 45° 𝑡𝑎𝑛 1 √𝑚
(σ1f (Units in MPa)
200
Slope Intercept
150
100
y = 4.4337x + 95.244
Intercept where σ3=0,
R² = 0.9735
thus on y-axis
50
0
‐25.0 ‐15.0 ‐5.0 5.0 15.0 25.0 35.0 45.0
σ3f (Minor Principal Stress) ‐ Units in MPa
Figure 40: Sigma 1 (Major Principal Stress σ1) vs Sigma 3 (Minor Principal Stress σ3)
with Linear Fit
𝝋
Figure 41: Shear Failure Occurs Along a Plane at Angle 𝜽 𝟒𝟓° to the Major
𝟐
Principal Plane
72
Figure 42: Mohr Circle and Cohesion (Shear Stress vs. Normal Stress) (Franki, 2008,
p. 53)
73
10.7 WORKED EXAMPLE 7: UCS (σci), mi, Tensile Strength (σti) and
Hoek-Brown & Theoretical Failure Envelope
You need to apply the Hoek-Brown criterion in a geotechnical evaluation. You have a set of
triaxial data and need to define the material constant mi, the UCS (σci) and the tensile strength
(σti) of the intact rock. The laboratory test results are presented in terms of principal stresses
σ1f and σ3f.
You need these parameters to evaluate a slope stability problem in a highly jointed rock mass
section. The software package you have in your office can only conduct evaluations making
use of Mohr-Coulomb input parameters, namely cohesion and friction. Your objective is to
determine the Hoek-Brown equivalent Mohr-Coulomb cohesion and friction parameters, in
order to ultimately evaluate the highly jointed rock mass.
It will become clear over the next working example or two that you will require the mi and σci
(UCS) of the intact rock. These values are then reduced to provide parameters for the rock
mass, with some changes to the shape/fitment parameters, s and a, from where we can obtain
our Hoek-Brown equivalent Mohr-Coulomb cohesion and friction angles. The concept is
clarified in Figure 43.
We have already calculated UCS (σci) in the previous worked example. In this working
example we will determine mi, σci and use this information to generate a theoretical Hoek-
Brown failure envelope. We will then calculate σti (the tensile strength of the intact rock).
Figure 43: Current Focus on Intact Rock Properties vs. Our Goal
Solution:
74
By knowing mi and σci (UCS) and setting s = 1 (for intact rock) and a = 0.5 (for intact rock) you
can determine σ1f for any value of σ3f or, alternatively, any value of σ3f for any value of σ1f.
The material constant mi and 𝜎 is determined from a plot of 𝜎 𝜎 (on the y-axis)
against 𝜎 (on the x-axis) (see Figure 44). You determine the linear relationship to obtain the
y-intercept and the slope of the line (m). To recap from the previous working example, σci (UCS
of the intact rock) is the square root of the intercept; mi is the slope (m) divided by the UCS
(σci). The plot and calculations are illustrated in Figure 44 and Worksheet 7.
The lab provided σ1f and σ3f from the six triaxial tests conducted under different confining
stresses (pressures in this case, as confinement was by means of oil pressure) (see Table 7).
You calculate (σ1f - σ3f)2 and plot this on the y-axis against σ1f on the x-axis (see Figure 44).
Table 7: Triaxial Test Results Conducted on Six Representative Isotropic Rock Samples
Determination of σci and mi from Triaxial Test Data
50000
45000
40000
(σ1f‐σ3f)^2 (Units in MPa)
35000
30000
Slope Intercept
25000
20000
y = 1070.4x + 7835
15000 R² = 0.9937
10000
5000 But do not forget me, you have done an actual UCS test in
the triaxial cell (where σ3=0 thus no confinement thus UCS).
0
‐10 0 10 20 30 40 50
σ3f (Minor Principal Stress) ‐ Units in MPa
75
You have now determined the UCS (σci) of 88.5 MPa and the mi of 12.1.
So why bother calculating the UCS, if you can just take the UCS of 78 MPa directly from the
triaxial test where σ3 = 0 (uniaxial compression, thus UCS test)?
It will become clear when you work with larger datasets that include hundreds or even
thousands of test results. We also need to calculate the slope to determine mi, and in excel it
requires the click of a button, so set yourself up to be able to interpret large datasets from the
start.
These tests were conducted on middle- to upper-bound strong rock crystalline limestone. Now,
compare the calculated results with the UCS ranges based on the rock strength description
(refer to Table 2 in the learning unit). Compare the calculated mi with the expected range of
mi values for a medium- to coarse-grained limestone (refer to Table 3 in the learning unit).
How do they compare?
For preliminary assessments, good observations by a suitably qualified and experienced
professional (generally your geologist, engineering geologist or geotechnical engineer) are of
great value.
By having mi and σci for the intact rock, and we have s = 1 and a = 0.5 for intact rock, we can
use the mi and σci (UCS) to derive the theoretical Hoek-Brown failure envelope for various
stress ranges. This enables us to simulate failure envelopes in ranges not generally practical
to test. The actual application comes in the rock mass scenario, where we use this same
approach to simulate rock mass failure envelopes.
With regard to the tensile strength of the intact rock: The tensile strength is expected to be the
intercept of the failure envelope on the x-axis. From the linear fit presented (see Figure 44) it
should be σti = 7.319 MPa. This seems a bit high, at 8.27% of the UCS. The linear relationship
overestimates the tensional strength. A tension cut-off should be applied by the relationship
σti = σci / (0.81 ⨯ mi + 7) (Hoek & Brown, 2018) (refer to Figure 21). This will give you σti of
5.27 MPa at ~5.95% of the UCS. We will consider a tensile strength of 5.27 MPa.
76
Theoretical Hoek‐Brown Failure Envelope (Actual Triaxial Data Included)
350.00
300.00
250.00
σ1f (MPa)
200.00
150.00
Theoretical Hoek‐Brown Failure Envelope from mi and σci
100.00
Actual Laboratory Test Data
Log. (Theoretical Hoek‐Brown Failure Envelope from mi and σci)
Poly. (Theoretical Hoek‐Brown Failure Envelope from mi and σci)
50.00
77
1 1 1 1
𝑎 𝑒 𝑒 𝑒 𝑒 ~0.5028405 0.503
2 6 2 6
You have s and a, so you can plug them into the formula to calculate UCSmass (σcm):
.
𝜎 𝜎 ∙𝑠 𝜎 ∙ 0.012 9.472 𝑀𝑃𝑎
The UCS of the rock mass is ~9 MPa vs. the UCS of the intact rock of ~88 MPa – a significant
difference!
Consider the calculation and impact of GSI. Consider intact rock, GSI = 100, calculate 𝜎 =
88.5 MPa. Can you recall this value? This was the σci (USC of the intact rock) that you have
78
already calculated. Refer to the GSI chart (see Figure 28). You will note a GSI of 100 is intact
rock. Before you apply, know the limits of this relationship!
Consider the influence of GSI on the σcm (UCS of the rock mass). You do the same calculation
and vary GSI from 0 to 100. The decrease in UCSmass with a decrease in rock mass quality
(decrease in GSI) is illustrated in Figure 46.
UCS of Rock Mass variation with GSI (UCS of Intact Rock = 88.5 MPa)
90
80
70
UCS of the Rock Mass MPa
60
50
40
30
20
10
0
0 10 20 30 40 50 60 70 80 90 100
Geological Strength Index (GSI)
Note: There are restrictions! Ensure that you understand the restrictions and recommended
ranges for applying any formulation, before doing so in practice!
To determine the uniaxial tensile strength of the rock mass, apply the following relationship:
𝑠𝜎
𝜎
𝑚
You have already determined the parameter “s” and the UCS of the intact rock (σci). You have
already calculated the mi for the intact rock (previous worked example) and need to decrease
the mi to represent the rock mass. We call this reduced mi, “mm” or “mb” (notations “m” and “b”
for mass or blocks). From the learning unit, you will know:
𝐺𝑆𝐼 100 60 100
𝑚 𝑚 𝑒𝑥𝑝 12.093 𝑒𝑥𝑝 2.898
28 14𝐷 28 14 0
𝑀𝑁
𝑠𝜎 0.0117 ∗ 88.5 𝑜𝑟 𝑀𝑃𝑎
𝜎 𝑚 0.359 𝑀𝑃𝑎 ~359 𝑘𝑁/𝑚
𝑚 2.898
Rock mass σcm = 7.479 MPa vs σtm = -0.359 MPa.
79
Figure 47: Current Focus on Intact Rock Properties vs. our Goal
We need to derive equivalent “c” and “φ” for the rock mass. From the learning unit you will
recall that the cohesion and friction of the rock mass can be calculated by:
⎛ ⎞
𝜎 , 𝜎 ,
⎜ 1 2∙𝑎 ∙𝑠 1 𝑎 ∙𝑚 ∙ 𝑠 𝑚 ∙ ⎟
𝜎 𝜎
𝑐 𝜎 ⎜ ⎟
⎜ 𝜎 ⎟
,
⎜ 6∙𝑎∙𝑚 𝑠 𝑚 ∙ ⎟
𝜎
1 𝑎 2 𝑎 ∙ 1
⎝ 𝑎 𝑎 2 𝑎 ⎠
80
⎡ 𝜎 , ⎤
⎢ 6∙𝑎∙𝑚 ∙ 𝑠 𝑚 ∙ ⎥
𝜎
𝜑 sin ⎢ ⎥
⎢2 ∙ 1 𝜎 , ⎥
𝑎 ∙ 2 𝑎 6∙𝑎∙𝑚 𝑠 𝑚 ∙
⎣ 𝜎 ⎦
The relationships may look intimidating at first. It is, however, very easy to formulate in excel.
All you need to calculate cohesion, are: a, s, mm, σ3 and σci (UCS).
These were all covered in the previous worked example. The only possible new concept may
be that there is a value of 𝜎 , , the upper limit of confining stress over which the relationship
between the Hoek-Brown and the Mohr-Coulomb criteria can be applied. This simply means
that there is a limit of applicability of the Hoek-Brown criterion. This concept may be familiar
to you. Mogi’s (1966) transition from brittle to ductile failure is defined as σ1=3.4*σ3. You can
thus consider any stress ratio inside the failure envelope, but only up to a maximum ration of
3.4. From here onwards the behaviour change from brittle to ductile and the Hoek-Brown
criterion no longer apply. Thus, if you formulate in excel, consider using the function
“min(σ3,σ3max)” to ensure that you only formulate up the maximum allowable ratio.
Figure 48: Brittle to Ductile Formation – Mogi’s σ1=3.4σ3 Limit (Hoek & Brown, 2018)
82
Illustration of Hoek‐Brown Calculation
GSI = 60 60 60 Geological Strength Index
D = 0 0 0 Damage factor applied
σci = 88.515 88.515 88.515 MPa
a = 0.503 0.503 0.503 HB Shape parameter
s = 0.0117 0.0117 0.0117 HB Shape parameter
mi = 12.0934 12.0934 12.0934
mm = 2.898 2.898 2.898 HB Parameter property of rock mass
σ3' = 0.000 1.350 1.350 MPa (y*H for scenario)
σcm = n/a 9.472 9.472
y = 0.027 0.027 0.027 Unit weight of rock (MN/m^3)
H = 50 50 50 Height or Depth (m)
Ductile transition = 4.00 σ3'max calc. σ3'max calc. Stress range applicable
σ3'max = 22.1288 0.7132 1.1583 Max minor principal stress
c' = 5.462 1.200 1.318 MPa <‐‐ Formulation error
φ' = 35.11 59.32 56.92 °
You now have cohesion and friction parameters which you can consider in the Mohr-Coulomb
shear strength criterion of 𝜏 𝑐 𝜎 ∙ 𝑡𝑎𝑛𝜑 for rock mass strength calculations. It is critical to
know the limitations of the application. You may apply the Hoek-Brown failure criterion to intact
rock and highly fractured rock mass, as indicated in Figure 47. Where the strength and
behaviour are controlled by individual or limited discontinuity planes, you should not apply the
Hoek-Brown criterion. In those cases, use shear strength of discontinuities (Mohr-Coulomb
and Barton-Choubey criteria) and kinematic analysis to determine whether
failure/sliding/movement is kinematically possible.
Note: You can now define cohesion and friction at different stress ratios (i.e., the non-linear
failure envelope).
83
𝐸 and 𝜎
If you are interested in deformation (thus strain “ε” as one of the components of deformation),
you need to have some idea of the applied load or stress (σ) and the stiffness “E” of the
material. The stress will depend on your problem, the depth of excavation, the height of slope
or the applied loading of a mega-structure on your rock mass. The rock mass property you
need is Erockmass or Em. Remember: Ei is for intact rock (established by means of triaxial testing).
We will consider two approaches (as discussed in the relevant section of the learning unit).
For the first equation we need Ei (the rest, D and GSI, we already have). You can establish
the Ei of the intact rock from the triaxial test results and apply the first formula, you can estimate
Ei from available information on the site (or from published data) or you can calculate Ei based
on the modulus ratio of the rock by the relationship Ei = MR*σci (not covered in the current
learning unit).
If you do not want to consider Ei in the calculation, consider formula 2. You should ideally
determine Ei from the triaxial results and consider more than one approach to increase your
level of confidence in the determined stiffness.
Let’s say we have determined or estimated Ei for the rock as 60 000 MPa. Now, for the highly
jointed rock mass, you should obtain an Em that is significantly less than 60 000 MPa.
The inputs into formulae 1 and 2 will thus be:
Ei = 60 000 MPa
D = 0 (no damage factor applied)
GSI = 60 (as provided by the experienced engineering geologist or geotechnical engineer)
𝐷
1
𝐸 𝐸 0.02 2
60 15𝐷 𝐺𝑆𝐼
1 exp
11
𝐷
1
𝐸 10 2
75 25𝐷 𝐺𝑆𝐼
1 exp
11
84
Worksheet 9: Illustration Calculation – Young’s Modulus for the Rock Mass (Rock
Mass Stiffness as Em) from Ei and GSI
85
If you follow the worked examples provided, you should be able to calculate the strength and
deformation parameters for intact rock and rock as a mass, to apply in a rock-engineering
problem. This should serve as a strong foundation on which to build and apply your
knowledge.
A snippet from RocLab (Version V1.031 - RocScience software) is provided in Figure 49.
Compare our calculated parameters with the parameters obtained by this software package.
Also see exports from the software package in Figure 50 that define the full failure envelope
in terms of principal stresses (σ1-σ3) and shear strength vs. normal stress (τ-σn). You should
now be equipped to construct them yourself. Just be aware of the ductile transition limit and
the tension cut-off! This was covered in the learning unit and in the calculation illustration in
this worked example.
Figure 49: RockLab Software Package Analysis Export for Comparison Purposes
86
Principal Stresses Normal Stress vs. Shear Stress
12000000 4000000
3500000
Major principal stress (MPa)
10000000
3000000
Shear stress (MPa)
8000000
2500000
6000000 2000000
1500000
4000000
1000000
2000000
500000
0 0
‐500000 0 500000 1000000 1500000 2000000 2500000 ‐1000000 0 1000000 2000000 3000000 4000000 5000000
Minor principal stress (MPa) Normal stress (MPa)
Figure 50: RockLab Software Package σ1-σ3 and τ-σn Plots for Illustrative Purposes
11 Self-Assessment Activities
Refer to SAPEM Chapter 3 and familiarise yourself with the following:
Page 52: Non-destructive tests on hardened concrete
Page 56: UCS apparatus
Page 43: Indirect tensile test (ITS)
Page 44: Triaxial test setup
Page 13: Durability tests
Page 16: ACV and 10% FACT test
Page 16: Ethylene Glycol soak test.
12 Further Reading
Further reading will be beneficial to you. Address the aspects introduced in this learning unit.
You can refer to the sources listed under the references. The listed sources are not
compulsory reading, and you may refer to alternative sources.
87
References
ASTM Standard D 5731, 1995. Standard Test Method for Determination of the Point Load
Strength Index of Rock, ASTM Designation: D 5731 ‐95. <<where can this be located,
author?>>
Barton, N., 2017. Barton‐Bandis Synopsis. [Online]
Available at: http://www.nickbarton.com/
de Vallejo, L., & Ferrer, M., 2011. Geological Engineering. Oxford: CRC Press; Taylor &
Francis. <<is my insert correct, author?>>
Deere, D. U., & Miller, R. P., 1966. Engineering Classification and Index Properties for Intact
Rock, s.l.: Air Force Weapons Laboratory Technical Report (AFWL‐TR). <<more details?>>
Franki, 2008. A Guide to Practical Geotechnical Engineering in Southern Africa (4th ed.). s.l.:
Franki. <<location?>>
Hoek, E., 2019. Practical Rock Engineering ‐ Chapter 11 ‐ Rock Mass Properties. [Online]
Available at: https://www.rocscience.com/learning/hoeks‐corner
Hoek, E., 2019. Practical Rock Engineering ‐ Chapter 4 ‐ Shear Strength of Discontinuities.
[Online]
Available at: https://www.rocscience.com/learning/hoeks‐corner
Hoek, E. & Brown, E. T., 2018. The Hoek‐Brown Failure Criterion and GSI 2018 Edition.
Journal of Rock Mechanics and Geotechnical Engineering, pp. 1‐19. <<volume and issue
number?>>
Hoek, E., Carranza‐Torres, C., & Corkum, B., 2002. Hoek‐Brown Failure Criterion ‐ 2002
Edition. << Place: Publisher?>>
Hoek, E., Carter, T. G., & Diederichs, M. S., 2013. Quantification of the Geological Strength
Index Chart. San Francisco: ARMA (Americal Rock Mechanics Association).
ISRM, 1979. Suggested Methods for Determining the Uniaxial Compressive Strength and
Deformability of Rock Materials. s.l.:International Society for Rock Mechanics. <<
location?>>
Mah, C. W., & Wyllie, D. C., 2004. Rock Slope Engineering ‐ Civil and Mining. (4th ed.).
London and New York: Spon Press; Taylor & Francis Group.
Marinos, P., & Hoek, E., 2000. GSI: Geologically Friendly Tool for Rock Mass Strength
Estimation. Melbourne, s.n., pp. 1422‐1442. <<publisher?>>
88
Sivakugan, N., Shukla, S. & Das, B., 2013. Rock Mechanics ‐ An Introduction. Oxford: CRC
Press; Taylor & Francis.
89