Calibration of Four Nonlinear Failure Envelopes from Triaxial Test Data and Influence of Nonlinearity on Geotechnical Computations

Show more

1. Introduction

Comprehensive literature reviews on nonlinear envelopes can be found in [1] [2] and a repetition here will be superfluous. A nonlinear or curved failure envelope is shown in Figure 1. The friction angle *ϕ* is a continuously varying quantity. The non-linearity of failure envelopes is most marked at lower stress levels, the failure envelope being of quasi-parabolic shape [3]. The most commonly used nonlinear

Figure 1. Curved failure envelope.

envelope is the power-type envelope $\tau ={\tau}_{o}{\left(1+\frac{{\sigma}^{\prime}}{{{\sigma}^{\prime}}_{t}}\right)}^{1/m}$ or $\tau ={\left(a+b{\sigma}^{\prime}\right)}^{n}$ where

*τ* = shear stress on failure arc;
${\sigma}^{\prime}$ = effective normal stress on failure arc = *σ *− *u*; *σ* = total normal stress; *u* = pore water pressure; *σ _{t}* = tensile strength of geomaterial;

$a>{\left[{b}^{2}n\left(1-2n\right)\right]}^{\frac{1}{2\left(1-n\right)}}$ (1)

prior to this discovery the power-type equation was affirmed to be a valid envelope for *n* > 0.5 only [2] [5].

Additionally, it was discovered that the quadratic equation

$\tau ={\left(a+b{\sigma}^{\prime}\right)}^{0.5}$ (2)

can only be a legitimate failure envelope if

$a\ge \frac{{b}^{2}}{4}$ (3)

*a* and *b* are parameters in the equation of quadratic failure envelope.

The aims of this manuscript are to

1) Calibrate the modified Maksimovic nonlinear failure envelope from triaxial test data;

2) Develop a simplified procedure for calibrating the polynomial type failure envelope from triaxial test data;

3) Develop a simple methodology for calibrating the power-type failure envelope from triaxial test data;

4) Produce a simpler procedure for calibrating the power-type failure envelope equivalent to the Hoek-Brown failure envelope using triaxial test data;

5) Compute the lateral earth pressure for a material with modified Maksimovic failure criterion;

6) Determine the factor of safety of a slope made of material with modified Maksimovic criterion;

7) Determine the factor of safety of a slope made of material with Hoek-Brown criterion.

This work has never been presented in the literature before. As had been mentioned real soil envelopes are nonlinear and the calibration of these nonlinear envelopes will help in determining the true response of soils. The accurate calibration of the Maksimovic failure envelope from triaxial test data has never been attempted before. The influences of nonlinearity on lateral earth pressure and factor of safety of slope with material made of modified Maksimovic law have never been attempted before and are quite difficult to implement.

2. Methodology

The calibration of a nonlinear failure envelope from triaxial test data is a problem of considerable difficulty [2] [7] [8]. To ease the calibration the formulas derived first by Balmer [9] come in handy.

They are:

${\sigma}^{\prime}={{\sigma}^{\prime}}_{3}+\frac{{{\sigma}^{\prime}}_{1}-{{\sigma}^{\prime}}_{3}}{1+\frac{\partial {{\sigma}^{\prime}}_{1}}{\partial {{\sigma}^{\prime}}_{3}}}$ (4)

$\tau =\frac{{{\sigma}^{\prime}}_{1}-{{\sigma}^{\prime}}_{3}}{1+\frac{\partial {{\sigma}^{\prime}}_{1}}{\partial {{\sigma}^{\prime}}_{3}}}\sqrt{\frac{\partial {{\sigma}^{\prime}}_{1}}{\partial {{\sigma}^{\prime}}_{3}}}$ (5)

where *σ'* and *τ* have their previous meanings,
${{\sigma}^{\prime}}_{1}$ = effective major principal stress at failure,
${{\sigma}^{\prime}}_{3}$ = effective minor principal stress at failure.

Given a set of experimental determined ( ${{\sigma}^{\prime}}_{3}$, ${{\sigma}^{\prime}}_{1}$ ) values then the normal stress and shear stress on the failure plane can be calculated from Equations (4) and (5) respectively.

It could be shown that

$\frac{\partial {{\sigma}^{\prime}}_{1}}{\partial {{\sigma}^{\prime}}_{3}}={N}_{\phi}=\frac{1+\mathrm{sin}{\phi}^{\prime}}{1-\mathrm{sin}{\phi}^{\prime}}$ (6)

where ${\phi}^{\prime}$ = the effective instantaneous friction angle.

These authors deduced that the effective instantaneous cohesion *c'* is given by

${c}^{\prime}=\frac{1}{2}\left[{{\sigma}^{\prime}}_{1}\sqrt{\frac{\partial {{\sigma}^{\prime}}_{3}}{\partial {{\sigma}^{\prime}}_{1}}}-{{\sigma}^{\prime}}_{3}\sqrt{\frac{\partial {{\sigma}^{\prime}}_{1}}{\partial {{\sigma}^{\prime}}_{3}}}\right]$ (7)

Equation (7) does not seem to exist in the literature. Fu and Liao (2010) seems to have derived the effective instantaneous cohesion for Hoek-Brown criterion that requires iteration to obtain.

2.1. Calibration of Modified Maksumovic Failure Envelope

Maksimovic’s [10] failure criterion has been determined to be excellent for rockfill.

According to Srbulov [11] this failure envelope provides the best fit to experimental data over a very wide stress range unlike the power-type envelopes that give low angle of friction at large normal stresses.

Maksimovic [10] proposed a hyperbolic failure law for rock-fills that can be expressed as

$\tau ={\sigma}^{\prime}\mathrm{tan}\left({a}_{1}+\frac{1}{{a}_{2}+{a}_{3}{\sigma}^{\prime}}\right)$ (8)

which has three parameters that needs to be determined. *a*_{1}, *a*_{2} and *a*_{3} are related to the Maksimovic’s parameters by *φ _{B}* =

Therefore a modified Makumovic hyperbolic law, that provided an improved match, was proposed namely

$\tau ={\sigma}^{\prime}\mathrm{tan}\left({a}_{1}+\frac{{\sigma}^{\prime}}{{a}_{2}+{a}_{3}{\sigma}^{\prime}}\right)$ (9)

For a drained test on a granular material
${\sigma}^{\prime}=\sigma $. By using Equations (4) and (5) values of *σ* and *τ* are computed in Table 1. Simple finite difference approximation has been used to compute values of
$\partial {\sigma}_{1}/\partial {\sigma}_{3}$. These seem to be crude approximations but in practice have been found to give excellent results.

By substituting for the largest values of *σ* and *τ* in Table 1, a relationship is derived between the constants, namely.

${a}_{1}+\frac{{\sigma}_{m}}{{a}_{2}+{a}_{3}{\sigma}_{m}}=\mathrm{arctan}\left(\frac{{\tau}_{m}}{{\sigma}_{m}}\right)$ (10)

with *σ _{m}* = 1105.69 and

Table 1. First two columns: Test data O from page 20 of Holtz and Gibbs (1952).

To obtain the best values of*a*_{2} and*a*_{3} the method of least squares is used. To implement the least squares method Equation (9) should be expressed as

$\frac{\sigma}{\mathrm{arctan}\left(\frac{\tau}{\sigma}\right)-{a}_{1}}=y={a}_{2}+{a}_{3}\sigma $ (11)

Applying the method of least squares it is obtained that

${a}_{2}N+{a}_{3}{\displaystyle \sum {\sigma}_{i}}={\displaystyle \sum {y}_{i}}$ (12a)

${a}_{2}{\displaystyle \sum {\sigma}_{i}}+{a}_{3}{\displaystyle \sum {\sigma}_{i}^{2}}={\displaystyle \sum {\sigma}_{i}{y}_{i}}$ (12b)

where subscript*i* will run from 1 to *N* = number of data points which is 6 in this case. Equation (12) can be solved to obtain

${a}_{3}=\frac{N{\displaystyle \sum \sigma y}-{\displaystyle \sum \sigma}{\displaystyle \sum y}}{N{\displaystyle \sum {\sigma}^{2}}-{\left({\displaystyle \sum \sigma}\right)}^{2}}$, (13a)

${a}_{2}=\frac{{\displaystyle \sum y}-{a}_{3}{\displaystyle \sum \sigma}}{N}$ (13b)

The best fit values of *a*_{2} and *a*_{3} should be obtained as follows: A certain value of *a*_{1} near one is chosen and then values of *a*_{2} and *a*_{3} are calculated as shown above. Equation (10) is used to calculate another value of *a*_{1}. if the chosen and the calculated *a*_{1} are equal then the correct solutions have been obtained. If they are different, then, another iteration should be done.

When the above routine is implemented it is obtained that the parameters for the modified Maksimovic failure law for the triaxial data given in Table 1 are

${a}_{1}=1.0$, ${a}_{2}=-64.35$, ${a}_{3}=-2.6837$.

The Mohr-Coulomb approximation to Table 1 is

$\tau =11.42+0.7269\sigma $ (14a)

with a SEE of 28.73. SEE = the standard error of estimate calculated from

$\text{SEE}=\sqrt{\frac{1}{N}{{\displaystyle \sum}}^{\text{}}{\left({\sigma}_{1pred}-{\sigma}_{1}\right)}^{2}}$

The modified Maksumovic approximation to Table 1 is

$\tau =\sigma \mathrm{tan}\left(1.0-\frac{\sigma}{64.35+2.6837\sigma}\right)$ (14b)

with a slightly lower SEE of 28.04 and the envelope passing through *τ* = 0 as required. The modified Maksimovic failure envelope would give a better factor of safety when the stability of geomaterial of shallow depth is considered. In Figure 2 the M-C and M-M failure envelopes are compared. In Table 2 is shown the data for test W of Holtz and Gibbs [12] for a sample containing 20% gravel at a relative density of 70%. The Mohr-Coulomb approximation to Table 2 is

$\tau =31.95+0.7542\sigma $ (15a)

with a SEE of 31.82.

The modified Maksimovic approximation to Table 2 is

Figure 2. The M-C and M-M failure envelopes are compared wrt Mohr circles of Test data O.

Table 2. Test data W from page 21 of Holtz and Gibbs (1952).

$\tau =\sigma \mathrm{tan}\left(1.0-\frac{\sigma}{254.55+2.7241\sigma}\right)$ (15b)

with a lower SEE of 23.78. This time around the difference is clear. The two envelopes are compared in Figure 3.

2.2. Calibration of Polynomial Failure Envelope

Yuanming *et al*. [13] derived a cubic polynomial *σ* − *τ* failure envelope for frozen sandy clay from *σ*_{3} − *σ*_{1} triaxial test data shown in Table 3. In this case effective normal stress = total normal stress. They calibrated the coefficients of the polynomial via a not really straightforward procedure. It is the purpose of this section to illustrate a much simpler procedure for doing the calibration.

Yuanming *et al*. [13] ’s procedure was proceeded by a derivation by regression of a *σ*_{1} − *σ*_{3} relation (Equation (16)) that enables an analytical determination of ∂*σ*_{1}/∂*σ*_{3} but these authors used a finite difference approximation for ∂*σ*_{1}/∂*σ*_{3}.

${\sigma}_{1}={\left({K}_{o}\right)}^{{\sigma}_{3}/{P}_{atm}}{\sigma}_{c}{\left(1+\frac{{\sigma}_{3}}{{\sigma}_{T}}\right)}^{{b}_{o}}$ (16)

Figure 3. The M-C and M-M failure envelopes are compared wrt Mohr circles of Test data W.

Table 3. Test data for frozen sany clay at −6˚C from Yuanming *et al*. (2010) page 51.

where *σ _{c}*and

*K _{o}*and

Using Equations (4) and (5) on Table 3 a set of values of *σ* and*τ* are obtained.

Let the cubic polynomial failure envelope be

$\tau ={b}_{1}+{b}_{2}\sigma +{b}_{3}{\sigma}^{2}+{b}_{4}{\sigma}^{3}$ (17)

Applying the method of least squares regression directly to this the following 4 equations are obtained

$\sum \tau}=N{b}_{1}+{b}_{2}{\displaystyle \sum \sigma}+{b}_{3}{\displaystyle \sum {\sigma}^{2}}+{b}_{4}{\displaystyle \sum {\sigma}^{3}$ (18a)

$\sum \sigma \tau}={b}_{1}{\displaystyle \sum \sigma}+{b}_{2}{\displaystyle \sum {\sigma}^{2}}+{b}_{3}{\displaystyle \sum {\sigma}^{3}}+{b}_{4}{\displaystyle \sum {\sigma}^{4}$ (18b)

$\sum {\sigma}^{2}\tau}={b}_{1}{\displaystyle \sum {\sigma}^{2}}+{b}_{2}{\displaystyle \sum {\sigma}^{3}}+{b}_{3}{\displaystyle \sum {\sigma}^{4}}+{b}_{4}{\displaystyle \sum {\sigma}^{5}$ (18c)

$\sum {\sigma}^{3}\tau}={b}_{1}{\displaystyle \sum {\sigma}^{3}}+{b}_{2}{\displaystyle \sum {\sigma}^{4}}+{b}_{3}{\displaystyle \sum {\sigma}^{5}}+{b}_{4}{\displaystyle \sum {\sigma}^{6}$ (18d)

where *N* is number of data points and the subscripts have been deleted without loss of meaning.

When Equations (18a) to (18d) are solved for the constants *b*_{1} to *b*_{4} it is obtained that

${b}_{1}=0.6549,{b}_{2}=0.6690,{b}_{3}=-0.04815,{b}_{4}=1.0773\times {10}^{-3}$

Yuanming *et al*. [13] obtained
${b}_{1}=0.6667$,
${b}_{2}=0.6290$,
${b}_{3}=-0.0423$,
${b}_{4}=0.90\times {10}^{-3}$

The two sets of constants can be seen to be approximately the same.

According to this paper the polynomial envelope is $\tau =0.6549+0.6690\sigma -0.04815{\sigma}^{2}+1.0773\times {10}^{-3}{\sigma}^{3}$

The standard error of estimate was determined to be SEE = 0.3727. The Mohr-Coulomb envelope is determined to be
$\tau =2.036+0.0969\sigma $ with an SEE = 1.1870. Thus, the Polynomial envelope provides a better fit to the experimental data. It is believed that the procedure used to do the calibration here is much more simpler and readily understood than the method of Yuanming *et al*.’s. A Q basic computer program [14] was used for doing all the relevant calculations and this program is available from the corresponding author on request. Figure 4 shows the Mohr circles of the experimental data of Table 3 and the associated M-C and polynomial failure envelopes.

2.3. Calibration of the Power-Type Failure Envelope

The power-type failure envelope takes the form $\tau ={\left(a+b\sigma \right)}^{n}$. A somehow complicated procedure for calibrating the power-type failure envelope has been given in Baker [1] and Anyaegbunam [2]. In Table 4 is shown experimental data lifted from Anyaegbunam [2].

Figure 4. Mohr circles of the experimental data of Table 3 and the associated M-C and polynomial failure envelopes.

Equations (4) and (5) are used to obtain the corresponding *σ* and *τ* values at failure. The power-type equation is then expressed as

$a+b\sigma ={\tau}^{1/n}$ (19)

By assuming several values of *n *Equation (19) is subjected to least squares regression for each value of *n*. The value of *n* that yields the minimum SEE is chosen as the correct failure envelope. This has been programmed in the afore mentioned QBASIC 4.5 program for the automatic determination of *n*, *a* and *b*. For the data of Table 4 it is obtained that
$\tau ={\left(0.389+2.61\sigma \right)}^{0.748}$ with an SEE of 0.25 that is almost exact match to the experimental data. This agrees almost exactly with Anyaegbunam [2] that gave
$\tau ={\left(0.439+2.612\sigma \right)}^{0.748}$ using a more complicated procedure. The Mohr-Coulomb envelope for the data of Table 4 is
$\tau =29.44+0.3812\sigma $ that corresponds to
$\phi =\mathrm{arctan}\left(0.3812\right)=20.9\u02da$. Figure 5 shows the Mohr circles of the experimental data of Table 4 and the associated linear M-C and power-type failure envelopes.

2.4. Calibration of the Power-Type Failure Envelope for the Hoek-Brown Criterion

Hoek and Brown [15] presented a useful and practical equation for the insitu

Table 4. Consolidated undrained test data for a sample of laterite from Table 2 of Anyaegbunam (2015).

Figure 5. Mohr circles of the experimental data of Table 4 and the associated linear M-C and power-type failure envelopes.

strength of rock mass. The 2002 version of this equation [16] (the modified Hoek-Brown criterion) is:

${\sigma}_{1}={\sigma}_{3}+{\sigma}_{ci}{\left(\frac{{m}_{b}{\sigma}_{3}}{{\sigma}_{ci}}+{s}_{b}\right)}^{\eta}$ (20)

where ${m}_{b}={m}_{i}\mathrm{exp}\left(\frac{\text{GSI}-100}{28-14D}\right)$, ${s}_{b}=\mathrm{exp}\left(\frac{\text{GSI}-100}{9-3D}\right)$

$\eta =\frac{1}{2}+\frac{1}{6}\left[\mathrm{exp}\left(-\text{GSI}/15\right)-\mathrm{exp}\left(-20/3\right)\right]$

where *GSI* is the Geological Strength Index, *D* is the disturbance factor, *σ** _{ci}* = the compressive strength of the intact rock,

The Hoek-Brown criterion can be expressed simply as

${\sigma}_{1}={\sigma}_{3}+{\left(\kappa {\sigma}_{3}+\beta \right)}^{\eta}$ (21)

where $\kappa ={m}_{b}{\sigma}_{ci}^{\frac{1}{\eta}-1}$ and $\beta ={s}_{b}{\sigma}_{ci}^{\frac{1}{\eta}}$

Letting $F=\kappa {\sigma}_{3}+\beta $ and $den=\kappa \eta +2{F}^{1-\eta}$ (22a, b)

then it can be shown that

$\sigma ={\sigma}_{3}+\frac{F}{den}$ and $\tau =\frac{{F}^{\frac{1+\eta}{2}}}{den}\sqrt{den-{F}^{1-\eta}}$ (23a, b)

Equations (23a, b) are the exact parametric (implicit) equations of the Mohr envelope of the H-B criterion. These equations with *σ*_{3} as the paramaeter are not very much useful in practice because the value of *τ* are not easily calculated from that of *σ*, vice versa. Hence, it is necessary to determine a single equation connecting *τ* and *σ*.

The H-B criterion utilized herein have the following constants *s** _{ci}*= 40 MN/m

The Hoek-Brown equation cannot be directly used in slope stability analysis because it is defined in terms of principal stresses. Considerable difficulty is encountered when Hoek-Brown equation is used directly in strength-reduction finite element type of slope stability analysis [17] [18] hence it will useful to obtain its Mohr envelope for use in limit equilibrium analysis. Baker [1] shows that the power-type failure envelope
$\tau ={\left(a+b\sigma \right)}^{n}$ provides an excellent fit to the Hoek-Brown criterion in the Mohr (*σ*-*τ*) plane. Therefore the Mohr envelope obtained in this section is of the power-type. Equation (21) has been used to derive Table 5 and Equations (23a and b) was used in the program to obtain values of *σ* and *τ*.

Figure 6 shows the Mohr circles of the H-B criterion and the associated Power-type envelope and M-C envelope. The derived power-type envelope has the equation $\tau ={\left(353.3+22.343\sigma \right)}^{0.684}$ with a SEE of 6.63. The Mohr-Coulomb envelope has the equation $\tau =159.86+0.8014\sigma $ or $\tau =159.86+\sigma \mathrm{tan}\left(38.71\u02da\right)$ with

Table 5. Data for Hoek-Brown criterion specified in the text.

Figure 6. Mohr circles of the H-B criterion and comparison with linear M-C and nonlinear power-type failure envelopes.

a SEE = 123.94. The envelopes derived in Anyaegbunam [2] namely
$\tau ={\left(357.302+22.337\sigma \right)}^{0.684}$ with a SEE = 5.77 and the M-C fit of *c* = 157.0 KN/m^{2} and *φ* = 38.7˚ with a SEE = 131.77 can be seen to be in excellent agreement. The method used in this paper can be seen to be much simpler and readily understood than Anyaegbunam [2]. Deng *et al*. [17] presented an approximate limit equilibrium technique for slope stability analysis using the Hoek-Brown criterion.

3. Influence of Nonlinearity of Failure Envelope on Geotechnical Computations

3.1. Passive Pressure on a Smooth Wall Due to Soil with Modified Maksimovic Law Derived from Table 2

This shall be illustrated for a 2.0 m smooth high wall that is a portion of a wall embedded in homogeneous soil of unit weight *γ* = 18 KN/m^{3}. The vertical stress, which is the minimum principal stress, is given by *σ*_{3} = *γz*(*i*) where *z*(*i*) = vertical distance from the ground surface. Assume that the wall height is divided into a number of divisions such that *σ*_{3} can be denoted by *σ*_{3}(*i*) and that the number of divisions is num. On the basis of a Mohr-Coulomb law the passive pressure at a given point is given by
${\sigma}_{1MC}\left(i\right)={\sigma}_{3}\left(i\right){N}_{\phi}+2C\sqrt{{N}_{\phi}}$. where *C* and *φ* are the

M-C shear strength parameters and
${N}_{\phi}={\mathrm{tan}}^{2}\left(\frac{\pi}{4}+\frac{\phi}{2}\right)$. *C* = 31.95 KN/m^{2}, tan*φ* = 0.7542, *a*_{1} = 1.0, *a*_{2} = −254.55, *a*_{3} = −2.7241.

The pseudo code for evaluating the passive pressure according to the modified Maksimovic law will be as follows:

Dimension the variables

set *i* = 1: eps = 0.001

450 compute *σ*_{3}(*i*): set *σ*_{3} = *σ*_{3}(*i*)’ 450 is line numbering in the code

compute *σ*_{1MC}(*i*): set
${\sigma}_{1MC}={\sigma}_{1MC}\left(i\right)$ (24)

Let *σ _{n}* = normal stress on the failure arc. Estimate

${\sigma}_{n1}=0.45\left[\left({\sigma}_{1MC}+{\sigma}_{3}\right)-\left({\sigma}_{1MC}-{\sigma}_{3}\right)\mathrm{sin}{\phi}_{MC}\right]$ _{ }

where ${\phi}_{MC}$ = M-C friction angle $\phi $

500 $\tau ={\sigma}_{n1}\mathrm{tan}\left({a}_{1}+\frac{{\sigma}_{n1}}{{a}_{2}+{a}_{3}{\sigma}_{n1}}\right)$

Let $der=\frac{\text{d}\tau}{\text{d}\sigma}$

$der=\frac{\tau}{{\sigma}_{n1}}+\frac{{a}_{2}{\sigma}_{n1}}{{\left({a}_{2}+{a}_{3}{\sigma}_{n1}\right)}^{2}}\left[1+\frac{{\tau}^{2}}{{\sigma}_{n1}^{2}}\right]$

$t=\tau \sqrt{1+de{r}^{2}}$

From Figure 1 it could be shown that

${\sigma}_{n2}={\sigma}_{n}={\sigma}_{3}+t-\sqrt{{t}^{2}-{\tau}^{2}}$

If i = num + 1 then END

If $\left|{\sigma}_{n2}-{\sigma}_{n1}\right|>eps$ then ${\sigma}_{n1}={\sigma}_{n2}$ : goto 500

If $\left|{\sigma}_{n2}-{\sigma}_{n1}\right|\le eps$ then ${\sigma}_{1}\left(i\right)={\sigma}_{3}+2t$ : $i=i+1$ : goto 450

In Figure 7 are shown the Passive pressures calculated for M-C and M. M. material having test data W.

The M-C soil has a passive force and moment of 401.4 KN/m and 353.1 KNm/m at the base of the wall. The M.M. soil has a passive force and moment of 284.3 KN and 203.5 KNm/m at the base of the wall which are 29.2% and 42.4% less than the M-C values respectively. This is usually the case for short walls that are less than 6.0 m in height.

3.2. Factor of Safety of a Slope with Material Governed by a Modified Maksimovic Law and Comparison with M-C Equivalent

It is proposed to determine the factor of safety of a 45˚ homogeneous slope of

Figure 7. Passive Pressures of soil with M-C and M.M. envelopes for test data W.

height 10 m via a Maksimovic material model and M-C material model using the Bishop simplified technique [19].

The material models were those derived from the triaxial data of Table 2 and are follows: For the M-C model: *C *= 31.95 KN/m^{2}, *ϕ* = 37.02˚ and for the M. M. model: *a*_{1} = 1.0, *a*_{2} = −254.55, *a*_{3} = −2.7241.

After the soil mass is divided into vertical slices the Bishop method for M-C material is given by

${F}_{s}=\frac{1}{{\displaystyle \sum \gamma {h}_{i}\mathrm{sin}{\alpha}_{i}}}{\displaystyle \sum \frac{{{c}^{\prime}}_{i}+\left(\gamma {h}_{i}-{u}_{i}\right)\mathrm{tan}{\phi}^{\prime}}{\mathrm{cos}{\alpha}_{i}\left(1+\frac{\mathrm{tan}{\alpha}_{i}\mathrm{tan}{\phi}^{\prime}}{{F}_{s}}\right)}}$ (25)

where *F _{s}* = factor of safety,

The Bishop method for modified Maksimovic (M.M.) material is given by

${F}_{s}=\frac{1}{{\displaystyle \sum \gamma {h}_{i}\mathrm{sin}{\alpha}_{i}}}{{\displaystyle \sum}}^{\text{}}\frac{{{\sigma}^{\prime}}_{i}}{\mathrm{cos}{\alpha}_{i}}\mathrm{tan}\left[{a}_{1}+\frac{{{\sigma}^{\prime}}_{i}}{{a}_{2}+{a}_{3}{{\sigma}^{\prime}}_{i}}\right]$ (26)

where ${{\sigma}^{\prime}}_{i}$ is obtained for each slice from the equation

${{\sigma}^{\prime}}_{i}+\mathrm{tan}{\alpha}_{i}\frac{{{\sigma}^{\prime}}_{i}}{{F}_{s}}\mathrm{tan}\left[{a}_{1}+\frac{{{\sigma}^{\prime}}_{i}}{{a}_{2}+{a}_{3}{{\sigma}^{\prime}}_{i}}\right]-\left(\gamma {h}_{i}-{u}_{i}\right)=0$ (27)

The factor of safety*F _{s}* is determined by assuming an initial value of

Using the M-C material model and 10 slices, the following results are obtained: *F _{s}* = 2.40 with the critical toe circle of radius 14.75 m centered at (−0.64, 14.74) with a central angle of 68.8˚. Using the M.M. material model the

Table 6. Comparison of F_{s} via Deng *et al*.’s H-B LEM and this paper’s H-B power-type LEM Bishop method.

with the critical toe circle of radius 16.87 m centered at (−4.46, 16.27) with a central angle of 52.8˚. The use of the M.M. material model is seen to have a strong influence on the calculated factor of safety which is seen to be smaller.

Unfortunately the M.M. law has not been used for stability calculations by other authors and comparison of the results herein cannot be compared with previous results.

3.3. Factor of Safety of Slope Using Deng *et al*.’s Approximate LEM with Hoek-Brown Parameters Compared with This Paper’s Rigorous Hoek-Brown Method

Deng *et al*. [17] developed an approximate limit equilibrium method (LEM) with Hoek–Brown parameters for determining the factor of safety of a rock slope.

The H-B parameters of the rock slope are *D* = 0, *GSI* = 100, *m _{i}* = 10,

Also, ${\sigma}_{cm}={\sigma}_{ci}\frac{\left[{m}_{b}+4{s}_{b}-\alpha \left({m}_{b}-8{s}_{b}\right)\right]{\left[\frac{{m}_{b}}{4}+{s}_{b}\right]}^{\alpha -1}}{2\left(1+\alpha \right)\left(2+\alpha \right)}$ (28)

${\sigma}_{3\mathrm{max}}=0.72{\sigma}_{cm}{\left(\frac{{\sigma}_{cm}}{\gamma H}\right)}^{-0.91}$ (29)

${\sigma}_{3}$ values are listed from values of 0 to
${\sigma}_{3\mathrm{max}}$ and the H-B relation is used to calculate corresponding values of
${\sigma}_{1}$. From these the power-type equivalent to this Hoek-Brown law is deduced by section (3.4) to be
$\tau ={\left(134.53+14.431\sigma \right)}^{0.643}$. The derived SEE = 1.26. In Table 6 is shown the results of factors of safety of slope stability obtained via the two methods. The F.S. for a slope using the power-type envelope was obtained using a modification of the method of Charles and Soares [20]. The F.S. of a slope using Deng *et al*.’s method was obtained using Equation (18) of their paper. The calculated factors of safety can be seen to be close except when *β* > 60˚ when they differ appreciably as Deng *et al.* discovered.

4. Conclusion

Four Mohr failure envelopes have been presented and the methods of for deriving them from experimental triaxial test data have been explained. The method for computing the passive resistance for granular soil with M.M. (modified Maksimovic) envelope is presented. In addition, the influence of a M.M. envelope on the factor of safety of slope is presented and it is shown that the use of a M.M. material model results in a significant reduction in the calculated factor of safety. Also, it is shown that the factors of safety of rock slope obtained by this paper’s power-type approximation and LEM are close to Deng *et al.*’s factors of safety using H-B approximation and LEM. A computer program written in QBASIC version 4.5 for doing all the calibrations has been developed by the authors and is available on request from the corresponding author. This program is a useful contribution to geotechnical engineering practice.

Acknowledgements

The authors wish to thank the University of Nigeria for the payment of the 2019 earned allowances to its staff that helped in funding this research.

References

[1] Baker, R. (2004) Non-Linear Strength Envelopes Based on Triaxial Test Data. Journal of Geotechnical & Geoenvironmental Engineering, 130, 498-506.

https://doi.org/10.1061/(ASCE)1090-0241(2004)130:5(498)

[2] Anyaegbunam, A.J. (2015) Nonlinear Power-Type Failure Laws for Geomaterials: Synthesis from Triaxial Data, Properties and Applications. International Journal of Geomechanics, 15, Article ID: 04014036.

https://doi.org/10.1061/(ASCE)GM.1943-5622.0000348

[3] Collins, I.F., Gunn, C.I.M., Pender, M.J. and Wang, Y. (1988) Slope Stability Analyses for Materials with a Nonlinear Failure Envelope. International Journal for Numerical and Analytical Methods in Geomechanics, 12, 533-550.

https://doi.org/10.1002/nag.1610120507

[4] Zhang, X.-J. and Chen, W.-F. (1987) Stability Analysis of Slopes with General Nonlinear Failure Criterion. International Journal for Numerical and Analytical Methods in Geomechanics, 11, 33-50.

https://doi.org/10.1002/nag.1610110104

[5] Jiang, J.-C., Baker, R. and Yamagami, T. (2003) The Effect of Strength Envelope Non-linearity on Slope Stability Computations. The Canadian Geotechnical Journal, 40, 308-325.

https://doi.org/10.1139/t02-111

[6] Li, X. (2007) Finite Element Analysis of Slope Stability Using a Nonlinear Failure Criterion. Computers and Geotechnics, 34, 127-136.

https://doi.org/10.1016/j.compgeo.2006.11.005

[7] Ucar, R. (1986) Determination of Shear Failure Envelope in Rock Masses. Journal of Geotechnical Engineering, 112, 303-315.

https://doi.org/10.1061/(ASCE)0733-9410(1986)112:3(303)

[8] Sheorey, P.R. (1997) Empirical Rock Failure Criteria. A.A. Balkema, Rotterdam.

[9] Balmer, G. (1952) A General Analytical Solution for Mohr’s Envelope. Proceedings of American Society of Test Materials, Vol. 52, 1260-1271.

[10] Maksimovic, M. (1989) Nonlinear Failure Envelope for Soils. Journal of Geotechnical & Geoenvironmental Engineering, 115, 581-586.

https://doi.org/10.1061/(ASCE)0733-9410(1989)115:4(581)

[11] Srbulov, M. (1997) On the Influence of Soil Strength Brittleness and Nonlinearity on Slope Stability. Computers and Geotechnics, 20, 95-104.

https://doi.org/10.1016/S0266-352X(96)00014-6

[12] Holtz, W.G. and Gibbs, H.J. (1956) Shear Strength of Pervious Gravelly Soils. Journal of the Soil Mechanics and Foundations Division, 82, 1-22.

https://doi.org/10.1061/JSFEAQ.0000004

[13] Lai, Y.M., et al. (2010) Stress-Strain Relationships and Nonlinear Mohr Strength Criteria of Frozen Sandy Clay. Soils and Foundations, 50, 45-53.

https://doi.org/10.3208/sandf.50.45

[14] QBASIC 4.5 [Computer Software]. Redmond, WA, Microsoft.

[15] Hoek, E. and Brown, E.T. (1980) Empirical Strength Criterion for Rock Masses. Journal of Geotechnical & Geoenvironmental Engineering, 106, 1013-1035.

https://doi.org/10.1061/AJGEB6.0001029

[16] Hammah, R.E., Curran, J.H., Yacoub, T.E. and Corkum, B. (2004) Stability Analysis of Rock Slopes Using the Finite Element Method. In: Schubert, W., Ed., Proceedings ISRM Regional Symposium EUROCK 2004 and 53rd Geomechanics Colloquy: Rock Engineering—Theory and Practice, International Society for Rock Mechanics, Lisbon, 783-788.

[17] Deng, D.P., Zhao, L.H. and Li, L. (2017) Limit Equilibrium Analysis for Rock Slope Stability Using Hoek-Brown Strength Criterion. Journal of Central South University, 24, 2154-2163.

https://doi.org/10.1007/s11771-017-3624-4

[18] Fu, W. and Liao, Y. (2010) Non-Linear Shear Strength Reduction Technique in Slope Stability Calculation. Computers and Geotechnics, 37, 288-298.

https://doi.org/10.1016/j.compgeo.2009.11.002

[19] Bishop, A.W. (1955) The Use of the Slip Circle in the Stability Analysis of Slopes. Géotechnique, 5, 7-17.

https://doi.org/10.1680/geot.1955.5.1.7

[20] Charles, J.A. and Soares, M.M. (1984) The Stability of Slopes in Soils with Nonlinear Failure Envelopes. Canadian Geotechnical Journal, 21, 397-406.

https://doi.org/10.1139/t84-044