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

Amaechi J. Anyaegbunam, Fidelis O. Okafor

Department of Civil Engineering, University of Nigeria, Nsukka, Nigeria.

DOI: 10.4236/gm.2021.112003                

Abstract

It is now recognized that many geomaterials have nonlinear failure envelopes. This non-linearity is most marked at lower stress levels, the failure envelope being of quasi-parabolic shape. It is not easy to calibrate these nonlinear failure envelopes from triaxial test data. Currently only the power-type failure envelope has been studied with an established formal procedure for its determination from triaxial test data. In this paper, a simplified procedure is evolved for the development of four different types of nonlinear envelopes. These are of invaluable assistance in the evaluation of true factors of safety in problems of slope stability and correct computation of lateral earth pressure and bearing capacity. The use of the Mohr-Coulomb failure envelopes leads to an overestimation of the factors of safety and other geotechnical quantities.

Keywords

Calibration of Nonlinear Failure Envelope, Triaxial Test Data, Modified Maksimovic Envelope, Power-Type Envelope, Polynomial-Type Envelope, Hoek-Brown Envelope, Standard Error of Estimate

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; τ_o = cohesion of geomaterial; m is a constant. In the alternative representation a, b, and n are constants with the case of a = 0 being common and denotes a purely granular geomaterial. This had been used in several studies on slope stability [4] [5] [6]. Investigations carried out by Anyaegbunam [2] revealed that the power-type equation is a valid envelope for soil for all n > 0 (except n = 0.5) if

$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+\sin {\phi }^{\prime }}{1-\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 }\tan \left(a_1+\frac{1}{a_2+a_3{\sigma }^{\prime }}\right)$ (8)

which has three parameters that needs to be determined. a₁, a₂ and a₃ are related to the Maksimovic’s parameters by φ_B = a₁, Δφ = 1/a₂, P_N = a₂/a₃. The determination of a₁, a₂ and a₃ is not straight forward at all. The parameters were determined for the experimental data presented in test No. O of Holtz and Gibbs [12] for a sample containing 20% gravel at a relative density of 50% shown in Table 1. It was discovered that the original Maksimovic [10] law Equation (8) does not provide an excellent match to the data.

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

$\tau ={\sigma }^{\prime }\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}=\arctan \left(\frac{{\tau }_m}{{\sigma }_m}\right)$ (10)

with σ_m = 1105.69 and τ_m = 814.85.

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

To obtain the best values ofa₂ anda₃ the method of least squares is used. To implement the least squares method Equation (9) should be expressed as

$\frac{\sigma }{\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 subscripti 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₂ and a₃ should be obtained as follows: A certain value of a₁ near one is chosen and then values of a₂ and a₃ are calculated as shown above. Equation (10) is used to calculate another value of a₁. if the chosen and the calculated a₁ 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 \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 \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 σ₃ − σ₁ 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 σ₁ − σ₃ relation (Equation (16)) that enables an analytical determination of ∂σ₁/∂σ₃ but these authors used a finite difference approximation for ∂σ₁/∂σ₃.

${\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 σ_cand σ_T are the uni-axial compressive and tensile strengths of the specimen respectively.

K_oand b_o are experimental parameters. Equation (16) is not needed in this study.

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

${\displaystyle \sum \tau }=N{b}_1+b_2{\displaystyle \sum \sigma }+b_3{\displaystyle \sum {\sigma }^2}+b_4{\displaystyle \sum {\sigma }^3}$ (18a)

${\displaystyle \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)

${\displaystyle \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)

${\displaystyle \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₁ to b₄ 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 =\arctan \left(0.3812\right)=20.9˚$. 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\exp \left(\frac{\text{GSI}-100}{28-14D}\right)$, $s_b=\exp \left(\frac{\text{GSI}-100}{9-3D}\right)$

$\eta =\frac{1}{2}+\frac{1}{6}\left[\exp \left(-\text{GSI}/15\right)-\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, m_i is a material constant 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 +2F^{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 σ₃ 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², m_i = 10, GSI = 45, D = 0.9 with the rest calculated to be m_b = 0.281, η = 0.508, s_b = 1.616 × 10⁻⁴.

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 \tan \left(38.71˚\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² 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³. The vertical stress, which is the minimum principal stress, is given by σ₃ = γ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 σ₃ can be denoted by σ₃(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 }={\tan }^2\left(\frac{\pi }{4}+\frac{\phi }{2}\right)$. C = 31.95 KN/m², tanφ = 0.7542, a₁ = 1.0, a₂ = −254.55, a₃ = −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 σ₃(i): set σ₃ = σ₃(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 σ_n as σ_n1.

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

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

500 $\tau ={\sigma }_{n1}\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², ϕ = 37.02˚ and for the M. M. model: a₁ = 1.0, a₂ = −254.55, a₃ = −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\sin {\alpha }_i}}{\displaystyle \sum \frac{{c^{\prime }}_i+\left(\gamma h_i-u_i\right)\tan {\phi }^{\prime }}{\cos {\alpha }_i\left(1+\frac{\tan {\alpha }_i\tan {\phi }^{\prime }}{F_s}\right)}}$ (25)

where F_s = factor of safety, h_i = midheight of slicei, α_i = inclination of base of slice i, u = pore pressure on base of slicei.

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

$F_s=\frac{1}{{\displaystyle \sum \gamma h_i\sin {\alpha }_i}}{{\displaystyle \sum }}^{\text{}}\frac{{{\sigma }^{\prime }}_i}{\cos {\alpha }_i}\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+\tan {\alpha }_i\frac{{{\sigma }^{\prime }}_i}{F_s}\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 safetyF_s is determined by assuming an initial value of F_s and calculating ${{\sigma }^{\prime }}_i$ from Equation (27). Thereafter ${{\sigma }^{\prime }}_i$ are substituted into Equation (26) to re-calculate F_s. The calculation of F_s using the Bishop method on a M.M. material model does not exist in the literature.

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 F_s= 1.64

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, σ_ci = 140 KNm⁻², α = η = 0.5. m_b and s_bare calculated from Equation (20). The rock unit weight is γ = 23 KN/m³.

Also, ${\sigma }_{cm}={\sigma }_{ci}\frac{\left[m_b+4s_b-\alpha \left(m_b-8s_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\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\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.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

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