Software Requirements Specification for Slope Stability analysis Program

Henry Frankis and Brooks MacLachlan

Table of Contents

An outline of all sections included in this SRS is recorded here for easy reference.

Reference Material

This section records information for easy reference.

Table of Units

The unit system used throughout is SI (Système International d'Unités). In addition to the basic units, several derived units are also used. For each unit, the Table of Units lists the symbol, a description, and the SI name.

Symbol Description SI Name
° angle degree
kg mass kilogram
m length metre
N force newton
Pa pressure pascal
s time second

Table of Units

Table of Symbols

The symbols used in this document are summarized in the Table of Symbols along with their units. Throughout the document, a subscript i indicates that the value will be taken at, and analyzed at, a slice or slice interface composing the total slip mass. For vector quantities, the units shown are for each component of the vector.

Symbol Description Units
A Area: A part of an object or surface. m2
a(t) Acceleration: The rate of change of a body's velocity. \(\frac{\text{m}}{\text{s}^{2}}\)
b Base Width of Slices: The width of each slice in the x-direction. m
Cden Proportionality Constant Denominator: Values for each slice that sum together to form the denominator of the interslice normal to shear force proportionality constant. N
Cnum Proportionality Constant Numerator: Values for each slice that sum together to form the numerator of the interslice normal to shear force proportionality constant. N
c′ Effective Cohesion: The internal pressure that sticks particles of soil together. Pa
const_f Decision on F: A Boolean decision on which form of f the user desires: constant if true, or half-sine if false. --
Fn Total Normal Force: Component of a force in the normal direction. N
Frot Force Causing Rotation: A force in the direction of rotation. N
FS Factor of Safety: The global stability metric of a slip surface of a slope, defined as the ratio of resistive shear force to mobilized shear force. --
Ft Tangential Force: Component of a force in the tangential direction. N
Fx x-coordinate of the Force: The force acting in the x-direction. N
Fy y-coordinate of the Force: The force acting in the y-direction. N
F Force: An interaction that tends to produce change in the motion of an object. N
FxG Sums of the Interslice Normal Forces: The sums of the normal forces acting on each pair of adjacent interslice boundaries. N
FxH Sums of the Interslice Normal Water Forces: The sums of the normal water forces acting on each pair of adjacent interslice boundaries. N
f Interslice Normal to Shear Force Ratio Variation Function: A function of distance in the x-direction that describes the variation of the interslice normal to shear ratio. --
G Interslice Normal Forces: The forces per meter in the z-direction exerted between each pair of adjacent slices. \(\frac{\text{N}}{\text{m}}\)
g Gravitational Acceleration: The approximate acceleration due to gravity on Earth at sea level. \(\frac{\text{m}}{\text{s}^{2}}\)
H Interslice Normal Water Forces: The normal water forces per meter in the z-direction exerted in the x-direction between each pair of adjacent slices. \(\frac{\text{N}}{\text{m}}\)
h Height: The distance above a reference point for a point of interest. m
h y-direction Heights of Slices: The heights in the y-direction from the base of each slice to the slope surface, at the x-direction midpoint of the slice. m
hL Heights of the Left Side of Slices: The heights of the left side of each slice, assuming slice surfaces have negative slope. m
hR Heights of the Right Side of Slices: The heights of the right side of each slice, assuming slice surfaces have negative slope. m
hz Heights of Interslice Normal Forces: The heights in the y-direction of the interslice normal forces on each slice. m
hz,w Heights of the Water Table: The heights in the y-direction from the base of each slice to the water table. m
i Index: A number representing a single slice. --
Unit Vector: A vector that has a magnitude of one. --
Kc Seismic Coefficient: The proportionality factor of force that weight pushes outwards; caused by seismic earth movements. --
Lb Total Base Lengths of Slices: The lengths of each slice in the direction parallel to the slope of the base. m
Ls Surface Lengths of Slices: The lengths of each slice in the direction parallel to the slope of the surface. m
M Moment: A measure of the tendency of a body to rotate about a specific point or axis. Nm
m Mass: The quantity of matter in a body. kg
N Normal Forces: The total reactive forces per meter in the z-direction for each slice of a soil surface subject to a body resting on it. \(\frac{\text{N}}{\text{m}}\)
N′ Effective Normal Forces: The forces per meter in the z-direction for each slice of a soil surface, subtracting pore water reactive force from total reactive force. \(\frac{\text{N}}{\text{m}}\)
n Number of Slices: The number of slices into which the slip surface is divided. --
P Resistive Shear Force: The Mohr Coulomb frictional force that describes the limit of mobilized shear force that can be withstood before failure. N
P Resistive Shear Forces: The Mohr Coulomb frictional forces per meter in the z-direction for each slice that describe the limit of mobilized shear force the slice can withstand before failure. \(\frac{\text{N}}{\text{m}}\)
p Pressure: A force exerted over an area. Pa
Q External Forces: The forces per meter in the z-direction acting into the surface from the midpoint of each slice. \(\frac{\text{N}}{\text{m}}\)
R Resistive Shear Forces Without the Influence of Interslice Forces: The resistive shear forces per meter without the influence of interslice forces in the z-direction for each slice. \(\frac{\text{N}}{\text{m}}\)
r Length of the Moment Arm: The distance between a force causing rotation and the axis of rotation. m
r Position Vector: A vector from the origin of the Cartesian coordinate system defined to the point where the force is applied. m
S Mobilized Shear Force: The shear force in the direction of potential motion. N
S Mobilized Shear Force: The mobilized shear force per meter in the z-direction for each slice. \(\frac{\text{N}}{\text{m}}\)
T Mobilized Shear Forces Without the Influence of Interslice Forces: The mobilized shear forces per meter without the influence of interslice forces in the z-direction for each slice. \(\frac{\text{N}}{\text{m}}\)
Ub Base Hydrostatic Forces: The forces per meter in the z-direction from water pressure within each slice. \(\frac{\text{N}}{\text{m}}\)
Ug Surface Hydrostatic Forces: The forces per meter in the z-direction from water pressure acting into each slice from standing water on the slope surface. \(\frac{\text{N}}{\text{m}}\)
u Pore Pressure: The pressure that comes from water within the soil. Pa
u Displacement: The change in an object's location relative to a reference point. m
V Volume: The amount of space that a substance or object occupies. m3
Vdry Volumes of Dry Soil: The amount of space occupied by dry soil for each slice. m3
Vsat Volumes of Saturated Soil: The amount of space occupied by saturated soil for each slice. m3
v Local Index: Used as a bound variable index in calculations. --
W Weight: The gravitational force acting on an object. N
W Weights: The downward force per meter in the z-direction on each slice caused by gravity. \(\frac{\text{N}}{\text{m}}\)
X Interslice Shear Forces: The shear forces per meter in the z-direction exerted between adjacent slices. \(\frac{\text{N}}{\text{m}}\)
x x-coordinate: The x-coordinate in the Cartesian coordinate system. m
xslipmaxEtr Maximum Entry x-coordinate: The maximum potential x-coordinate for the entry point of a slip surface. m
xslipmaxExt Maximum Exit x-coordinate: The maximum potential x-coordinate for the exit point of a slip surface. m
xslipminEtr Minimum Entry x-coordinate: The minimum potential x-coordinate for the entry point of a slip surface. m
xslipminExt Minimum Exit x-coordinate: The minimum potential x-coordinate for the exit point of a slip surface. m
xcs,ycs Critical Slip Surface Coordinates: The set of x-coordinates and y-coordinates that describe the vertices of the critical slip surface. m
xslip x-coordinates of the Slip Surface: x-coordinates of points on the slip surface. m
xslope x-coordinates of the Slope: x-coordinates of points on the soil slope. m
xwt x-coordinates of the Water Table: X-positions of the water table. m
y y-coordinate: The y-coordinate in the Cartesian coordinate system. m
yslipmax Maximum y-coordinate: The maximum potential y-coordinate of a point on a slip surface. m
yslipmin Minimum y-coordinate: The minimum potential y-coordinate of a point on a slip surface. m
yslip y-coordinates of the Slip Surface: Heights of the slip surface. m
yslope y-coordinates of the Slope: y-coordinates of points on the soil slope. m
ywt y-coordinates of the Water Table: Heights of the water table. m
z z-coordinate: The z-coordinate in the Cartesian coordinate system. m
α Base Angles: The angles between the base of each slice and the horizontal. °
β Surface Angles: The angles between the surface of each slice and the horizontal. °
γ Specific Weight: The weight per unit volume. \(\frac{\text{N}}{\text{m}^{3}}\)
γdry Soil Dry Unit Weight: The weight of a dry soil/ground layer divided by the volume of the layer. \(\frac{\text{N}}{\text{m}^{3}}\)
γsat Soil Saturated Unit Weight: The weight of saturated soil/ground layer divided by the volume of the layer. \(\frac{\text{N}}{\text{m}^{3}}\)
γw Unit Weight of Water: The weight of one cubic meter of water. \(\frac{\text{N}}{\text{m}^{3}}\)
λ Proportionality Constant: The ratio of the interslice normal to the interslice shear force. --
π Ratio of Circumference to Diameter for Any Circle: The ratio of a circle's circumference to its diameter. --
ρ Density: The mass per unit volume. \(\frac{\text{kg}}{\text{m}^{3}}\)
σ Total Normal Stress: The total force per area acting on the soil mass. Pa
σ′ Effective Stress: The stress in a soil mass that is effective in causing volume changes and mobilizes the shear strength arising from friction; represents the average stress carried by the soil skeleton. Pa
σN Effective Normal Stress: The normal stress in a soil mass that is effective in causing volume changes; represents the average normal stress carried by the soil skeleton. Pa
τ Tangential Stress: The shear force per unit area. Pa
τf Shear Strength: The strength of a material against shear failure. Pa
τ Torque: A twisting force that tends to cause rotation. Nm
Φ First Function for Incorporating Interslice Forces Into Shear Force: The function for converting resistive shear without the influence of interslice forces, to a calculation considering the interslice forces. --
φ′ Effective Angle of Friction: The angle of inclination with respect to the horizontal axis of the Mohr-Coulomb shear resistance line. °
Ψ Second Function for Incorporating Interslice Forces Into Shear Force: The function for converting mobile shear without the influence of interslice forces, to a calculation considering the interslice forces. --
ω Imposed Load Angles: The angles between the external force acting into the surface of each slice and the vertical. °

Table of Symbols

Abbreviations and Acronyms

Abbreviation Full Form
2D Two-Dimensional
3D Three-Dimensional
A Assumption
DD Data Definition
GD General Definition
GS Goal Statement
IM Instance Model
LC Likely Change
PS Physical System Description
R Requirement
RefBy Referenced by
Refname Reference Name
SRS Software Requirements Specification
SSP Slope Stability analysis Program
TM Theoretical Model
UC Unlikely Change
Uncert. Typical Uncertainty

Abbreviations and Acronyms

Introduction

A slope of geological mass, composed of soil and rock and sometimes water, is subject to the influence of gravity on the mass. This can cause instability in the form of soil or rock movement. The effects of soil or rock movement can range from inconvenient to seriously hazardous, resulting in significant life and economic losses. Slope stability is of interest both when analysing natural slopes, and when designing an excavated slope. Slope stability analysis is the assessment of the safety of a slope, identifying the surface most likely to experience slip and an index of its relative stability known as the factor of safety.

The following section provides an overview of the Software Requirements Specification (SRS) for a slope stability analysis problem. The developed program will be referred to as the Slope Stability analysis Program (SSP). This section explains the purpose of this document, the scope of the requirements, the characteristics of the intended reader, and the organization of the document.

Purpose of Document

The primary purpose of this document is to record the requirements of SSP. Goals, assumptions, theoretical models, definitions, and other model derivation information are specified, allowing the reader to fully understand and verify the purpose and scientific basis of SSP. With the exception of system constraints, this SRS will remain abstract, describing what problem is being solved, but not how to solve it.

This document will be used as a starting point for subsequent development phases, including writing the design specification and the software verification and validation plan. The design document will show how the requirements are to be realized, including decisions on the numerical algorithms and programming environment. The verification and validation plan will show the steps that will be used to increase confidence in the software documentation and the implementation. Although the SRS fits in a series of documents that follow the so-called waterfall model, the actual development process is not constrained in any way. Even when the waterfall model is not followed, as Parnas and Clements point out parnasClements1986, the most logical way to present the documentation is still to "fake" a rational design process.

Scope of Requirements

The scope of the requirements includes stability analysis of a two-dimensional (2D) soil mass, composed of a single homogeneous layer with constant material properties. The soil mass is assumed to extend infinitely in the third dimension. The analysis will be at an instant in time; factors that may change the soil properties over time will not be considered.

Characteristics of Intended Reader

Reviewers of this documentation should have an understanding of undergraduate level 4 physics and undergraduate level 2 or higher solid mechanics. It would be an asset to understand soil mechanics. The users of SSP can have a lower level of expertise, as explained in Sec:User Characteristics.

Organization of Document

The organization of this document follows the template for an SRS for scientific computing software proposed by koothoor2013, smithLai2005, smithEtAl2007, and smithKoothoor2016. The presentation follows the standard pattern of presenting goals, theories, definitions, and assumptions. For readers that would like a more bottom up approach, they can start reading the instance models and trace back to find any additional information they require.

The goal statements are refined to the theoretical models and the theoretical models to the instance models. The instance models provide the set of algebraic equations that must be solved.

General System Description

This section provides general information about the system. It identifies the interfaces between the system and its environment, describes the user characteristics, and lists the system constraints.

System Context

Fig:sysCtxDiag shows the system context. A circle represents an external entity outside the software. A rectangle represents the software system itself (SSP). Arrows are used to show the data flow between the system and its environment.

System Context
System Context

The responsibilities of the user and the system are as follows:

  • User Responsibilities
    • Provide the input data related to the soil layer(s) and water table (if applicable), ensuring conformation to input data format required by SSP
    • Ensure that consistent units are used for input variables
    • Ensure required software assumptions are appropriate for the problem to which the user is applying the software
  • SSP Responsibilities
    • Detect data type mismatch, such as a string of characters input instead of a floating point number
    • Verify that the inputs satisfy the required physical and other data constraints
    • Identify the critical slip surface within the possible input range
    • Find the factor of safety for the slope
    • Find the interslice normal force and shear force along the critical slip surface

User Characteristics

The end user of SSP should have an understanding of undergraduate Level 1 Calculus and Physics, and be familiar with soil and material properties, specifically effective cohesion, effective angle of friction, and unit weight.

System Constraints

The Morgenstern-Price method morgenstern1965, which involves dividing the slope into vertical slices, will be used to derive the equations for analysing the slope.

Specific System Description

This section first presents the problem description, which gives a high-level view of the problem to be solved. This is followed by the solution characteristics specification, which presents the assumptions, theories, and definitions that are used.

Problem Description

A system is needed to evaluate the factor of safety of a slope's slip surface and identify the critical slip surface of the slope, as well as the interslice normal force and shear force along the critical slip surface.

Terminology and Definitions

This subsection provides a list of terms that are used in the subsequent sections and their meaning, with the purpose of reducing ambiguity and making it easier to correctly understand the requirements.

  • Factor of safety: The global stability metric of a slip surface of a slope, defined as the ratio of resistive shear force to mobilized shear force.
  • Slip surface: A surface within a slope that has the potential to fail or displace due to load or other forces.
  • Critical slip surface: Slip surface of the slope that has the lowest factor of safety, and is therefore most likely to experience failure.
  • Water table: The upper boundary of a saturated zone in the ground.
  • Stress: The ratio of an applied force to a cross-sectional area.
  • Strain: A measure of deformation representing the displacement between particles in the body relative to a reference length.
  • Normal force: A force applied perpendicular to the plane of the material.
  • Shear force: A force applied parallel to the plane of the material.
  • Mobilized shear force: The shear force in the direction of potential motion, thus encouraging motion along the plane.
  • Resistive shear force: The shear force in the direction opposite to the direction of potential motion, thus hindering motion along the plane.
  • Effective forces and stresses: The normal force or normal stress carried by the soil skeleton, composed of the effective force or stress and the force or stress exerted by water.
  • Cohesion: An attractive force between adjacent particles that holds the matter together.
  • Isotropy: A condition where the value of a property is independent of the direction in which it is measured.
  • Plane strain: A condition where the resultant stresses in one of the directions of a three-dimensional material can be approximated as zero. This condition results when a body is constrained to not deform in one direction, or when the length of one dimension of the body dominates the others, to the point where it can be assumed as infinite. Stresses in the direction of the dominant dimension can be approximated as zero.

Physical System Description

The physical system of SSP, as shown in Fig:PhysicalSystem, includes the following elements:

PS1: A slope comprised of one soil layer.

PS2: A water table, which may or may not exist.

An example slope for analysis by SSP, where the dashed line represents the water table
An example slope for analysis by SSP, where the dashed line represents the water table

Morgenstern-Price analysis morgenstern1965 of the slope involves representing the slope as a series of vertical slices. As shown in Fig:IndexConvention, the index i is used to denote a value for a single slice, and an interslice value at a given index i refers to the value between slice i and adjacent slice i+1.

Index convention for slice and interslice values
Index convention for slice and interslice values

A free body diagram of the forces acting on a slice is displayed in Fig:ForceDiagram. The specific forces and symbols will be discussed in detail in Sec:General Definitions and Sec:Data Definitions.

Free body diagram of forces acting on a slice
Free body diagram of forces acting on a slice

Goal Statements

Given the shape of the soil mass, the location of the water table, and the material properties of the soil, the goal statements are:

Identify-Crit-and-FS: Identify the critical slip surface and the corresponding factor of safety.

Determine-Normal-Forces: Determine the interslice normal forces between each pair of vertical slices of the slope.

Determine-Shear-Forces: Determine the interslice shear forces between each pair of vertical slices of the slope.

Solution Characteristics Specification

The instance models that govern SSP are presented in the Instance Model Section. The information to understand the meaning of the instance models and their derivation is also presented, so that the instance models can be verified.

Assumptions

This section simplifies the original problem and helps in developing the theoretical models by filling in the missing information for the physical system. The assumptions refine the scope by providing more detail.

Slip-Surface-Concave: The slip surface is concave with respect to the slope surface. The (xslip, yslip) coordinates of a slip surface follow a concave up function. (RefBy: IM:crtSlpId.)

Factor-of-Safety: The factor of safety is assumed to be constant across the entire slip surface. (RefBy: GD:mobShr.)

Soil-Layer-Homogeneous: The soil mass is homogeneous, with consistent soil properties throughout. (RefBy: GD:sliceWght, GD:resShr, and LC:Calculate-Inhomogeneous-Soil-Layers.)

Soil-Properties: The soil properties are independent of dry or saturated conditions, with the exception of unit weight. (RefBy: GD:resShr.)

Soil-Layers-Isotropic: The soil mass is treated as if the effective cohesion and effective angle of friction are isotropic properties. (RefBy: GD:resShr.)

Interslice-Norm-Shear-Forces-Linear: Following the assumption of Morgenstern and Price (morgenstern1965), interslice normal forces and interslice shear forces have a proportional relationship, depending on a proportionality constant (λ) and a function (f) describing variation depending on x position. (RefBy: GD:normShrR, IM:fctSfty, IM:nrmShrFor, and UC:Normal-And-Shear-Linear-Only.)

Plane-Strain-Conditions: The slope and slip surface extends far into and out of the geometry (z coordinate). This implies plane strain conditions, making 2D analysis appropriate. (RefBy: GD:srfWtrF, GD:sliceWght, GD:resShr, GD:effNormF, and GD:baseWtrF.)

Effective-Norm-Stress-Large: The effective normal stress is large enough that the shear strength to effective normal stress relationship can be approximated as a linear relationship. (RefBy: TM:equilibrium and UC:2D-Analysis-Only.)

Surface-Base-Slice-between-Interslice-Straight-Lines: The surface and base of a slice are approximated as straight lines. (RefBy: GD:srfWtrF, GD:sliceWght, GD:baseWtrF, TM:mcShrStrgth, DD:slcHeight, DD:angleB, and DD:angleA.)

Edge-Slices: The interslice forces at the 0th and nth interslice interfaces are zero. (RefBy: IM:intsliceFs, IM:fctSfty, and IM:nrmShrFor.)

Seismic-Force: There is no seismic force acting on the slope. (RefBy: IM:fctSfty, IM:nrmShrFor, and LC:Calculate-Seismic-Force.)

Surface-Load: There is no imposed surface load, and therefore no external forces, acting on the slope. (RefBy: IM:fctSfty, IM:nrmShrFor, and LC:Calculate-External-Force.)

Water-Intersects-Base-Edge: The water table only intersects the base of a slice at an edge of the slice. (RefBy: GD:sliceWght and GD:baseWtrF.)

Water-Intersects-Surface-Edge: The water table only intersects the slope surface at the edge of a slice. (RefBy: GD:srfWtrF and GD:sliceWght.)

Negligible-Effect-Surface-Slope-Seismic: The effect of the slope of the surface of the soil on the seismic force is assumed to be negligible. (RefBy: GD:momentEql.)

Hydrostatic-Force-Slice-Midpoint: The resultant surface hydrostatic forces act into the midpoint of each slice surface and the resultant base hydrostatic forces act into the midpoint of each slice base. (RefBy: GD:srfWtrF, GD:momentEql, and GD:baseWtrF.)

Theoretical Models

This section focuses on the general equations and laws that SSP is based on.

Refname TM:factOfSafety
Label

Factor of safety

Equation \[{F_{\text{S}}}=\frac{P}{S}\]
Description
  • FS is the factor of safety (Unitless)
  • P is the resistive shear force (N)
  • S is the mobilized shear force (N)
Source

fredlund1977

RefBy

GD:mobShr

Refname TM:equilibrium
Label

Equilibrium

Equation \[\displaystyle\sum{{F_{\text{x}}}}=0\] \[\displaystyle\sum{{F_{\text{y}}}}=0\] \[\displaystyle\sum{M}=0\]
Description
  • Fx is the x-coordinate of the force (N)
  • Fy is the y-coordinate of the force (N)
  • M is the moment (Nm)
Notes

For a body in static equilibrium, the net forces and moments acting on the body will cancel out. Assuming a 2D problem (A:Effective-Norm-Stress-Large), the x-coordinate of the force Fx and y-coordinate of the force Fy will be equal to 0. All forces and their distance from the chosen point of rotation will create a net moment equal to 0.

Source

fredlund1977

RefBy

GD:normForcEq, GD:momentEql, and GD:bsShrFEq

Refname TM:mcShrStrgth
Label

Mohr-Coulumb shear strength

Equation \[{τ^{\text{f}}}={σ_{N}}' \tan\left(φ'\right)+c'\]
Description
  • τf is the shear strength (Pa)
  • σN is the effective normal stress (Pa)
  • φ′ is the effective angle of friction (°)
  • c′ is the effective cohesion (Pa)
Notes

In this model the shear strength τf is proportional to the product of the effective normal stress σN on the plane with its static friction in the angular form tan(φ′). The τf versus σN relationship is not truly linear, but assuming the effective normal forces is strong enough, it can be approximated with a linear fit (A:Surface-Base-Slice-between-Interslice-Straight-Lines) where the effective cohesion c′ represents the τf intercept of the fitted line.

Source

fredlund1977

RefBy

GD:resShr

Refname TM:effStress
Label

Effective stress

Equation \[σ'=σ-u\]
Description
  • σ′ is the effective stress (Pa)
  • σ is the total normal stress (Pa)
  • u is the pore pressure (Pa)
Notes

σ is defined in DD:normStress.

Source

fredlund1977

RefBy

GD:effNormF

Refname TM:NewtonSecLawMot
Label

Newton's second law of motion

Equation \[\symbf{F}=m \symbf{a}\text{(}t\text{)}\]
Description
  • F is the force (N)
  • m is the mass (kg)
  • a(t) is the acceleration (\(\frac{\text{m}}{\text{s}^{2}}\))
Notes

The net force F on a body is proportional to the acceleration a(t) of the body, where m denotes the mass of the body as the constant of proportionality.

Source

--

RefBy

GD:weight

General Definitions

This section collects the laws and equations that will be used to build the instance models.

Refname GD:normForcEq
Label

Normal force equilibrium

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[{\symbf{N}}_{i}=\left({\symbf{W}}_{i}-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \cos\left({\symbf{ω}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)+\left(-{K_{\text{c}}} {\symbf{W}}_{i}-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)\]
Description
  • N is the normal forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • W is the weights (\(\frac{\text{N}}{\text{m}}\))
  • X is the interslice shear forces (\(\frac{\text{N}}{\text{m}}\))
  • Ug is the surface hydrostatic forces (\(\frac{\text{N}}{\text{m}}\))
  • β is the surface angles (°)
  • Q is the external forces (\(\frac{\text{N}}{\text{m}}\))
  • ω is the imposed load angles (°)
  • α is the base angles (°)
  • Kc is the seismic coefficient (Unitless)
  • G is the interslice normal forces (\(\frac{\text{N}}{\text{m}}\))
  • H is the interslice normal water forces (\(\frac{\text{N}}{\text{m}}\))
Notes

This equation satisfies TM:equilibrium in the normal direction. W is defined in GD:sliceWght, Ug is defined in GD:srfWtrF, β is defined in DD:angleB, and α is defined in DD:angleA.

Source

chen2005

RefBy

IM:fctSfty

Normal force equilibrium is derived from the free body diagram of Fig:ForceDiagram in Sec:Physical System Description.

Refname GD:bsShrFEq
Label

Base shear force equilibrium

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[{\symbf{S}}_{i}=\left({\symbf{W}}_{i}-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \cos\left({\symbf{ω}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)-\left(-{K_{\text{c}}} {\symbf{W}}_{i}-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)\]
Description
  • S is the mobilized shear force (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • W is the weights (\(\frac{\text{N}}{\text{m}}\))
  • X is the interslice shear forces (\(\frac{\text{N}}{\text{m}}\))
  • Ug is the surface hydrostatic forces (\(\frac{\text{N}}{\text{m}}\))
  • β is the surface angles (°)
  • Q is the external forces (\(\frac{\text{N}}{\text{m}}\))
  • ω is the imposed load angles (°)
  • α is the base angles (°)
  • Kc is the seismic coefficient (Unitless)
  • G is the interslice normal forces (\(\frac{\text{N}}{\text{m}}\))
  • H is the interslice normal water forces (\(\frac{\text{N}}{\text{m}}\))
Notes

This equation satisfies TM:equilibrium in the shear direction. W is defined in GD:sliceWght, Ug is defined in GD:srfWtrF, β is defined in DD:angleB, and α is defined in DD:angleA.

Source

chen2005

RefBy

IM:fctSfty

Base shear force equilibrium is derived from the free body diagram of Fig:ForceDiagram in Sec:Physical System Description.

Refname GD:resShr
Label

Resistive shear force

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[{\symbf{P}}_{i}={\symbf{N'}}_{i} \tan\left({φ'}_{i}\right)+{c'}_{i} {\symbf{L}_{b,i}}\]
Description
  • P is the resistive shear forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • N′ is the effective normal forces (\(\frac{\text{N}}{\text{m}}\))
  • φ′ is the effective angle of friction (°)
  • c′ is the effective cohesion (Pa)
  • Lb is the total base lengths of slices (m)
Notes

Lb is defined in DD:lengthLb.

Source

chen2005

RefBy

GD:mobShr

Derived by substituting DD:normStress and DD:tangStress into the Mohr-Coulomb shear strength, TM:mcShrStrgth, and multiplying both sides of the equation by the area of the slice in the shear-z plane. Since the slope is assumed to extend infinitely in the z-direction (A:Plane-Strain-Conditions), the resulting forces are expressed per metre in the z-direction. The effective angle of friction φ′ and the effective cohesion c′ are not indexed by i because they are assumed to be isotropic (A:Soil-Layers-Isotropic) and the soil is assumed to be homogeneous, with constant soil properties throughout (A:Soil-Layer-Homogeneous, A:Soil-Properties).

Refname GD:mobShr
Label

Mobilized shear force

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[{\symbf{S}}_{i}=\frac{{\symbf{P}}_{i}}{{F_{\text{S}}}}=\frac{{\symbf{N'}}_{i} \tan\left({φ'}_{i}\right)+{c'}_{i} {\symbf{L}_{b,i}}}{{F_{\text{S}}}}\]
Description
  • S is the mobilized shear force (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • P is the resistive shear forces (\(\frac{\text{N}}{\text{m}}\))
  • FS is the factor of safety (Unitless)
  • N′ is the effective normal forces (\(\frac{\text{N}}{\text{m}}\))
  • φ′ is the effective angle of friction (°)
  • c′ is the effective cohesion (Pa)
  • Lb is the total base lengths of slices (m)
Notes

Lb is defined in DD:lengthLb.

Source

chen2005

RefBy

IM:fctSfty

Mobilized shear forces is derived by dividing the definition of the P from GD:resShr by the definition of the factor of safety from TM:factOfSafety. The factor of safety FS is not indexed by i because it is assumed to be constant for the entire slip surface (A:Factor-of-Safety).

Refname GD:effNormF
Label

Effective normal force

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[{\symbf{N'}}_{i}={\symbf{N}}_{i}-{\symbf{U}_{\text{b},i}}\]
Description
  • N′ is the effective normal forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • N is the normal forces (\(\frac{\text{N}}{\text{m}}\))
  • Ub is the base hydrostatic forces (\(\frac{\text{N}}{\text{m}}\))
Notes

Ub is defined in GD:baseWtrF.

Source

chen2005

RefBy

Derived by substituting DD:normStress into TM:effStress and multiplying both sides of the equation by the area of the slice in the shear-z plane. Since the slope is assumed to extend infinitely in the z-direction (A:Plane-Strain-Conditions), the resulting forces are expressed per metre in the z-direction.

Refname GD:resShearWO
Label

Resistive shear force, without interslice normal and shear forces

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[{\symbf{R}}_{i}=\left(\left({\symbf{W}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)+\left(-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)-{\symbf{U}_{\text{b},i}}\right) \tan\left({φ'}_{i}\right)+{c'}_{i} {\symbf{L}_{b,i}}\]
Description
  • R is the resistive shear forces without the influence of interslice forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • W is the weights (\(\frac{\text{N}}{\text{m}}\))
  • Ug is the surface hydrostatic forces (\(\frac{\text{N}}{\text{m}}\))
  • β is the surface angles (°)
  • α is the base angles (°)
  • H is the interslice normal water forces (\(\frac{\text{N}}{\text{m}}\))
  • Ub is the base hydrostatic forces (\(\frac{\text{N}}{\text{m}}\))
  • φ′ is the effective angle of friction (°)
  • c′ is the effective cohesion (Pa)
  • Lb is the total base lengths of slices (m)
Notes

W is defined in GD:sliceWght, Ug is defined in GD:srfWtrF, β is defined in DD:angleB, α is defined in DD:angleA, H is defined in DD:intersliceWtrF, Ub is defined in GD:baseWtrF, and Lb is defined in DD:lengthLb.

Source

chen2005 and karchewski2012

RefBy

IM:intsliceFs and IM:fctSfty

Refname GD:mobShearWO
Label

Mobilized shear force, without interslice normal and shear forces

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[{\symbf{T}}_{i}=\left({\symbf{W}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)-\left(-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)\]
Description
  • T is the mobilized shear forces without the influence of interslice forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • W is the weights (\(\frac{\text{N}}{\text{m}}\))
  • Ug is the surface hydrostatic forces (\(\frac{\text{N}}{\text{m}}\))
  • β is the surface angles (°)
  • α is the base angles (°)
  • H is the interslice normal water forces (\(\frac{\text{N}}{\text{m}}\))
Notes

W is defined in GD:sliceWght, Ug is defined in GD:srfWtrF, β is defined in DD:angleB, α is defined in DD:angleA, and H is defined in DD:intersliceWtrF.

Source

chen2005 and karchewski2012

RefBy

IM:intsliceFs and IM:fctSfty

Refname GD:normShrR
Label

Interslice shear forces

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[\symbf{X}=λ \symbf{f} \symbf{G}\]
Description
  • X is the interslice shear forces (\(\frac{\text{N}}{\text{m}}\))
  • λ is the proportionality constant (Unitless)
  • f is the interslice normal to shear force ratio variation function (Unitless)
  • G is the interslice normal forces (\(\frac{\text{N}}{\text{m}}\))
Notes

Mathematical representation of the primary assumption for the Morgenstern-Price method (A:Interslice-Norm-Shear-Forces-Linear). f is defined in DD:ratioVariation.

Source

chen2005

RefBy

IM:fctSfty and IM:nrmShrFor

Refname GD:momentEql
Label

Moment equilibrium

Units

N

Equation \[0=-{\symbf{G}}_{i} \left({\symbf{h}_{\text{z},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+{\symbf{G}}_{i-1} \left({\symbf{h}_{\text{z},i-1}}-\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)-{\symbf{H}}_{i} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+{\symbf{H}}_{i-1} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i-1}}-\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+\frac{{\symbf{b}}_{i}}{2} \left({\symbf{X}}_{i}+{\symbf{X}}_{i-1}\right)+\frac{-{K_{\text{c}}} {\symbf{W}}_{i} {\symbf{h}}_{i}}{2}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right) {\symbf{h}}_{i}+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right) {\symbf{h}}_{i}\]
Description
  • G is the interslice normal forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • hz is the heights of interslice normal forces (m)
  • b is the base width of slices (m)
  • α is the base angles (°)
  • H is the interslice normal water forces (\(\frac{\text{N}}{\text{m}}\))
  • hz,w is the heights of the water table (m)
  • X is the interslice shear forces (\(\frac{\text{N}}{\text{m}}\))
  • Kc is the seismic coefficient (Unitless)
  • W is the weights (\(\frac{\text{N}}{\text{m}}\))
  • h is the y-direction heights of slices (m)
  • Ug is the surface hydrostatic forces (\(\frac{\text{N}}{\text{m}}\))
  • β is the surface angles (°)
  • Q is the external forces (\(\frac{\text{N}}{\text{m}}\))
  • ω is the imposed load angles (°)
Notes

This equation satisfies TM:equilibrium for the net moment. b is defined in DD:lengthB, α is defined in DD:angleA, W is defined in GD:sliceWght, h is defined in DD:slcHeight, Ug is defined in GD:srfWtrF, and β is defined in DD:angleB.

Source

chen2005

RefBy

IM:nrmShrFor

Moment is equal to torque, so the equation from DD:torque will be used to calculate moments:

\[\symbf{τ}=\symbf{u}\times\symbf{F}\]

Considering one dimension, with moments in the clockwise direction as positive and moments in the counterclockwise direction as negative, and replacing the torque symbol with the moment symbol, the equation simplifies to:

\[M={F_{\text{rot}}} r\]

where Frot is the force causing rotation and r is the length of the moment arm, or the distance between the force and the axis about which the rotation acts. To represent the moment equilibrium, the moments from each force acting on a slice must be considered and added together. The forces acting on a slice are all shown in Fig:ForceDiagram. The midpoint of the base of a slice is considered as the axis of rotation, from which the length of the moment arm is measured. Considering first the interslice normal force acting on slice interface i, the moment is negative because the force tends to rotate the slice in a counterclockwise direction, and the length of the moment arm is the height of the force plus the difference in height between the base at slice interface i and the base at the midpoint of slice i. Thus, the moment is expressed as:

\[-{\symbf{G}}_{i} \left({\symbf{h}_{\text{z},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)\]

For the i−1th slice interface, the moment is similar but in the opposite direction:

\[{\symbf{G}}_{i-1} \left({\symbf{h}_{\text{z},i-1}}-\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)\]

Next, the interslice normal water force is considered. This force is zero at the height of the water table, then increases linearly towards the base of the slice due to the increasing water pressure. For such a triangular distribution, the resultant force acts at one-third of the height. Thus, for the interslice normal water force acting on slice interface i, the moment is:

\[-{\symbf{H}}_{i} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)\]

The moment for the interslice normal water force acting on slice interface i−1 is:

\[{\symbf{H}}_{i-1} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i-1}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)\]

The interslice shear force at slice interface i tends to rotate in the clockwise direction, and the length of the moment arm is the length from the slice edge to the slice midpoint, equivalent to half of the width of the slice, so the moment is:

\[{\symbf{X}}_{i} \frac{{\symbf{b}}_{i}}{2}\]

The interslice shear force at slice interface i−1 also tends to rotate in the clockwise direction, and has the same length of the moment arm, so the moment is:

\[{\symbf{X}}_{i-1} \frac{{\symbf{b}}_{i}}{2}\]

Seismic forces act over the entire height of the slice. For each horizontal segment of the slice, the seismic force is KcWi where Wi can be expressed as γ bi y using GD:weight where y is the height of the segment under consideration. The corresponding length of the moment arm is y, the height from the base of the slice to the segment under consideration. In reality, the forces near the surface of the soil mass are slightly different due to the slope of the surface, but this difference is assumed to be negligible (A:Negligible-Effect-Surface-Slope-Seismic). The resultant moment from the forces on all of the segments with an equivalent resultant length of the moment arm is determined by taking the integral over the slice height. The forces tend to rotate in the counterclockwise direction, so the moment is negative:

\[-\int_{0}^{{\symbf{h}}_{i}}{{K_{\text{c}}} γ {\symbf{b}}_{i} y}\,dy\]

Solving the definite integral yields:

\[-{K_{\text{c}}} γ {\symbf{b}}_{i} \frac{{\symbf{h}}_{i}^{2}}{2}\]

Using GD:weight again to express γ bihi as Wi, the moment is:

\[-{K_{\text{c}}} {\symbf{W}}_{i} \frac{{\symbf{h}}_{i}}{2}\]

The surface hydrostatic force acts into the midpoint of the surface of the slice (A:Hydrostatic-Force-Slice-Midpoint). Thus, the vertical component of the force acts directly towards the point of rotation, and has a moment of zero. The horizontal component of the force tends to rotate in a clockwise direction and the length of the moment arm is the entire height of the slice. Thus, the moment is:

\[{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right) {\symbf{h}}_{i}\]

The external force again acts into the midpoint of the slice surface, so the vertical component does not contribute to the moment, and the length of the moment arm is again the entire height of the slice. The moment is:

\[{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right) {\symbf{h}}_{i}\]

The base hydrostatic force and slice weight both act in the direction of the point of rotation (A:Hydrostatic-Force-Slice-Midpoint), therefore both have moments of zero. Thus, all of the moments have been determined. The moment equilibrium is then represented by the sum of all moments:

\[0=-{\symbf{G}}_{i} \left({\symbf{h}_{\text{z},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+{\symbf{G}}_{i-1} \left({\symbf{h}_{\text{z},i-1}}-\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)-{\symbf{H}}_{i} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+{\symbf{H}}_{i-1} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i-1}}-\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+\frac{{\symbf{b}}_{i}}{2} \left({\symbf{X}}_{i}+{\symbf{X}}_{i-1}\right)+\frac{-{K_{\text{c}}} {\symbf{W}}_{i} {\symbf{h}}_{i}}{2}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right) {\symbf{h}}_{i}+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right) {\symbf{h}}_{i}\]
Refname GD:weight
Label

Weight

Units

N

Equation \[W=V γ\]
Description
  • W is the weight (N)
  • V is the volume (m3)
  • γ is the specific weight (\(\frac{\text{N}}{\text{m}^{3}}\))
Source

Definition of Weight

RefBy

GD:sliceWght and GD:momentEql

Detailed derivation of weight:

Under the influence of gravity, and assuming a 2D Cartesian coordinate system with down as positive, an object has an acceleration vector of:

\[\symbf{a}\text{(}t\text{)}=\begin{bmatrix} 0\\ \symbf{g} \symbf{\hat{j}} \end{bmatrix}\]

Since there is only one non-zero vector component, the scalar value W will be used for the weight. In this scenario, Newton's second law of motion from TM:NewtonSecLawMot can be expressed as:

\[W=m \symbf{g}\]

Mass can be expressed as density multiplied by volume, resulting in:

\[W=ρ V \symbf{g}\]

Substituting specific weight as the product of density and gravitational acceleration yields:

\[W=V γ\]
Refname GD:sliceWght
Label

Slice weight

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[{\symbf{W}}_{i}={\symbf{b}}_{i} \frac{1}{2} \begin{cases} \left({\symbf{y}_{\text{slope},i}}-{\symbf{y}_{\text{slip},i}}+{\symbf{y}_{\text{slope},i-1}}-{\symbf{y}_{\text{slip},i-1}}\right) {γ_{\text{sat}}}, & {\symbf{y}_{\text{wt},i}}\gt{}{\symbf{y}_{\text{slope},i}}\lor{}{\symbf{y}_{\text{wt},i-1}}\gt{}{\symbf{y}_{\text{slope},i-1}}\\ \left({\symbf{y}_{\text{slope},i}}-{\symbf{y}_{\text{wt},i}}+{\symbf{y}_{\text{slope},i-1}}-{\symbf{y}_{\text{wt},i-1}}\right) {γ_{\text{dry}}}+\left({\symbf{y}_{\text{wt},i}}-{\symbf{y}_{\text{slip},i}}+{\symbf{y}_{\text{wt},i-1}}-{\symbf{y}_{\text{slip},i-1}}\right) {γ_{\text{sat}}}, & {\symbf{y}_{\text{slope},i}}\geq{}{\symbf{y}_{\text{wt},i}}\geq{}{\symbf{y}_{\text{slip},i}}\land{}{\symbf{y}_{\text{slope},i-1}}\geq{}{\symbf{y}_{\text{wt},i-1}}\geq{}{\symbf{y}_{\text{slip},i-1}}\\ \left({\symbf{y}_{\text{slope},i}}-{\symbf{y}_{\text{slip},i}}+{\symbf{y}_{\text{slope},i-1}}-{\symbf{y}_{\text{slip},i-1}}\right) {γ_{\text{dry}}}, & {\symbf{y}_{\text{wt},i}}\lt{}{\symbf{y}_{\text{slip},i}}\lor{}{\symbf{y}_{\text{wt},i-1}}\lt{}{\symbf{y}_{\text{slip},i-1}} \end{cases}\]
Description
  • W is the weights (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • b is the base width of slices (m)
  • yslope is the y-coordinates of the slope (m)
  • yslip is the y-coordinates of the slip surface (m)
  • γsat is the soil saturated unit weight (\(\frac{\text{N}}{\text{m}^{3}}\))
  • ywt is the y-coordinates of the water table (m)
  • γdry is the soil dry unit weight (\(\frac{\text{N}}{\text{m}^{3}}\))
Notes

This equation is based on the assumption that the surface and the base of a slice are straight lines (A:Surface-Base-Slice-between-Interslice-Straight-Lines). The soil dry unit weight γdry and the soil saturated unit weight γsat are not indexed by i because the soil is assumed to be homogeneous, with constant soil properties throughout (A:Soil-Layer-Homogeneous). b is defined in DD:lengthB.

Source

fredlund1977

RefBy

GD:resShearWO, GD:normForcEq, GD:momentEql, GD:mobShearWO, and GD:bsShrFEq

For the case where the water table is above the slope surface, the weights come from the weight of the saturated soil. Substituting values for saturated soil into the equation for weight from GD:weight yields:

\[{\symbf{W}}_{i}={\symbf{V}_{\text{sat},i}} {γ_{\text{sat}}}\]

Due to A:Plane-Strain-Conditions, only two dimensions are considered, so the areas of saturated soil are considered instead of the volumes of saturated soil. Any given slice has a trapezoidal shape. The area of a trapezoid is the average of the lengths of the parallel sides multiplied by the length between the parallel sides. The parallel sides in this case are the interslice edges and the length between them is the width of the slice. Thus, the weights are defined as:

\[{\symbf{W}}_{i}={\symbf{b}}_{i} \frac{1}{2} \left({\symbf{y}_{\text{slope},i}}-{\symbf{y}_{\text{slip},i}}+{\symbf{y}_{\text{slope},i-1}}-{\symbf{y}_{\text{slip},i-1}}\right) {γ_{\text{sat}}}\]

For the case where the water table is below the slip surface, the weights come from the weight of the dry soil. Substituting values for dry soil into the equation for weight from GD:weight yields:

\[{\symbf{W}}_{i}={\symbf{V}_{\text{dry},i}} {γ_{\text{dry}}}\]

A:Plane-Strain-Conditions again allows for two-dimensional analysis so the areas of dry soil are considered instead of the volumes of dry soil. The trapezoidal slice shape is the same as in the previous case, so the weights are defined as:

\[{\symbf{W}}_{i}={\symbf{b}}_{i} \frac{1}{2} \left({\symbf{y}_{\text{slope},i}}-{\symbf{y}_{\text{slip},i}}+{\symbf{y}_{\text{slope},i-1}}-{\symbf{y}_{\text{slip},i-1}}\right) {γ_{\text{dry}}}\]

For the case where the water table is between the slope surface and slip surface, the weights are the sums of the weights of the dry portions and weights of the saturated portions of the soil. Substituting values for dry and saturated soil into the equation for weight from GD:weight and adding them together yields:

\[{\symbf{W}}_{i}={\symbf{V}_{\text{dry},i}} {γ_{\text{dry}}}+{\symbf{V}_{\text{sat},i}} {γ_{\text{sat}}}\]

A:Plane-Strain-Conditions again allows for two-dimensional analysis so the areas of dry soil and areas of saturated soil are considered instead of the volumes of dry soil and volumes of saturated soil. The water table is assumed to only intersect a slice surface or base at a slice edge (A:Water-Intersects-Surface-Edge, A:Water-Intersects-Base-Edge), so the dry and saturated portions each have trapezoidal shape. For the dry portion, the parallel sides of the trapezoid are the lengths between the slope surface and water table at the slice edges. For the saturated portion, the parallel sides of the trapezoid are the lengths between the water table and slip surface at the slice edges. Thus, the weights are defined as:

\[{\symbf{W}}_{i}={\symbf{b}}_{i} \frac{1}{2} \left(\left({\symbf{y}_{\text{slope},i}}-{\symbf{y}_{\text{wt},i}}+{\symbf{y}_{\text{slope},i-1}}-{\symbf{y}_{\text{wt},i-1}}\right) {γ_{\text{dry}}}+\left({\symbf{y}_{\text{wt},i}}-{\symbf{y}_{\text{slip},i}}+{\symbf{y}_{\text{wt},i-1}}-{\symbf{y}_{\text{slip},i-1}}\right) {γ_{\text{sat}}}\right)\]
Refname GD:hsPressure
Label

Hydrostatic pressure

Units

Pa

Equation \[p=γ h\]
Description
  • p is the pressure (Pa)
  • γ is the specific weight (\(\frac{\text{N}}{\text{m}^{3}}\))
  • h is the height (m)
Notes

This equation is derived from Bernoulli's equation for a slow moving fluid through a porous material.

Source

Definition of Pressure

RefBy

GD:srfWtrF and GD:baseWtrF

Refname GD:baseWtrF
Label

Base hydrostatic force

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[{\symbf{U}_{\text{b},i}}={\symbf{L}_{b,i}} {γ_{w}} \frac{1}{2} \begin{cases} {\symbf{y}_{\text{wt},i}}-{\symbf{y}_{\text{slip},i}}+{\symbf{y}_{\text{wt},i-1}}-{\symbf{y}_{\text{slip},i-1}}, & {\symbf{y}_{\text{wt},i}}\gt{}{\symbf{y}_{\text{slip},i}}\lor{}{\symbf{y}_{\text{wt},i-1}}\gt{}{\symbf{y}_{\text{slip},i-1}}\\ 0, & {\symbf{y}_{\text{wt},i}}\leq{}{\symbf{y}_{\text{slip},i}}\land{}{\symbf{y}_{\text{wt},i-1}}\leq{}{\symbf{y}_{\text{slip},i-1}} \end{cases}\]
Description
  • Ub is the base hydrostatic forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • Lb is the total base lengths of slices (m)
  • γw is the unit weight of water (\(\frac{\text{N}}{\text{m}^{3}}\))
  • ywt is the y-coordinates of the water table (m)
  • yslip is the y-coordinates of the slip surface (m)
Notes

This equation is based on the assumption that the base of a slice is a straight line (A:Surface-Base-Slice-between-Interslice-Straight-Lines). Lb is defined in DD:lengthLb.

Source

fredlund1977

RefBy

GD:resShearWO, GD:effNormF, and GD:baseWtrF

The base hydrostatic forces come from the hydrostatic pressure exerted by the water above the base of each slice. The equation for hydrostatic pressure from GD:hsPressure is:

\[p=γ h\]

The specific weight in this case is the unit weight of water γw. The height in this case is the height from the slice base to the water table. This height is measured from the midpoint of the slice because the resultant hydrostatic force is assumed to act at the slice midpoint (A:Hydrostatic-Force-Slice-Midpoint). The height at the midpoint is the average of the height at slice interface i and the height at slice interface i−1:

\[\frac{1}{2} \left({\symbf{y}_{\text{wt},i}}-{\symbf{y}_{\text{slip},i}}+{\symbf{y}_{\text{wt},i-1}}-{\symbf{y}_{\text{slip},i-1}}\right)\]

Due to A:Plane-Strain-Conditions, only two dimensions are considered, so the base hydrostatic forces are expressed as forces per meter. The pressures acting on the slices can thus be converted to base hydrostatic forces by multiplying by the corresponding length of the slice base Lb,i, assuming the water table does not intersect a slice base except at a slice edge (A:Water-Intersects-Base-Edge). Thus, in the case where the height of the water table is above the height of the slip surface, the base hydrostatic forces are defined as:

\[{\symbf{U}_{\text{b},i}}={\symbf{L}_{b,i}} {γ_{w}} \frac{1}{2} \left({\symbf{y}_{\text{wt},i}}-{\symbf{y}_{\text{slip},i}}+{\symbf{y}_{\text{wt},i-1}}-{\symbf{y}_{\text{slip},i-1}}\right)\]

This equation is the non-zero case of GD:baseWtrF. The zero case is when the height of the water table is below the height of the slip surface, so there is no hydrostatic force.

Refname GD:srfWtrF
Label

Surface hydrostatic force

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[{\symbf{U}_{\text{g},i}}={\symbf{L}_{s,i}} {γ_{w}} \frac{1}{2} \begin{cases} {\symbf{y}_{\text{wt},i}}-{\symbf{y}_{\text{slope},i}}+{\symbf{y}_{\text{wt},i-1}}-{\symbf{y}_{\text{slope},i-1}}, & {\symbf{y}_{\text{wt},i}}\gt{}{\symbf{y}_{\text{slope},i}}\lor{}{\symbf{y}_{\text{wt},i-1}}\gt{}{\symbf{y}_{\text{slope},i-1}}\\ 0, & {\symbf{y}_{\text{wt},i}}\leq{}{\symbf{y}_{\text{slope},i}}\land{}{\symbf{y}_{\text{wt},i-1}}\leq{}{\symbf{y}_{\text{slope},i-1}} \end{cases}\]
Description
  • Ug is the surface hydrostatic forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • Ls is the surface lengths of slices (m)
  • γw is the unit weight of water (\(\frac{\text{N}}{\text{m}^{3}}\))
  • ywt is the y-coordinates of the water table (m)
  • yslope is the y-coordinates of the slope (m)
Notes

This equation is based on the assumption that the surface of a slice is a straight line (A:Surface-Base-Slice-between-Interslice-Straight-Lines). Ls is defined in DD:lengthLs.

Source

fredlund1977

RefBy

GD:srfWtrF, GD:resShearWO, IM:nrmShrForNum, GD:normForcEq, GD:momentEql, GD:mobShearWO, and GD:bsShrFEq

The surface hydrostatic forces come from the hydrostatic pressure exerted by the water above the surface of each slice. The equation for hydrostatic pressure from GD:hsPressure is:

\[p=γ h\]

The specific weight in this case is the unit weight of water γw. The height in this case is the height from the slice surface to the water table. This height is measured from the midpoint of the slice because the resultant hydrostatic force is assumed to act at the slice midpoint (A:Hydrostatic-Force-Slice-Midpoint). The height at the midpoint is the average of the height at slice interface i and the height at slice interface i−1:

\[\frac{1}{2} \left({\symbf{y}_{\text{wt},i}}-{\symbf{y}_{\text{slope},i}}+{\symbf{y}_{\text{wt},i-1}}-{\symbf{y}_{\text{slope},i-1}}\right)\]

Due to A:Plane-Strain-Conditions, only two dimensions are considered, so the surface hydrostatic forces are expressed as forces per meter. The pressures acting on the slices can thus be converted to surface hydrostatic forces by multiplying by the corresponding length of the slice surface Ls,i, assuming the water table does not intersect a slice surface except at a slice edge (A:Water-Intersects-Surface-Edge). Thus, in the case where the height of the water table is above the height of the slope surface, the surface hydrostatic forces are defined as:

\[{\symbf{U}_{\text{g},i}}={\symbf{L}_{s,i}} {γ_{w}} \frac{1}{2} \left({\symbf{y}_{\text{wt},i}}-{\symbf{y}_{\text{slope},i}}+{\symbf{y}_{\text{wt},i-1}}-{\symbf{y}_{\text{slope},i-1}}\right)\]

This equation is the non-zero case of GD:srfWtrF. The zero case is when the height of the water table is below the height of the slope surface, so there is no hydrostatic force.

Data Definitions

This section collects and defines all the data needed to build the instance models.

Refname DD:intersliceWtrF
Label

Interslice normal water forces

Symbol

H

Units

\(\frac{\text{N}}{\text{m}}\)

Equation \[\symbf{H}=\begin{cases} \frac{\left({\symbf{y}_{\text{slope},i}}-{\symbf{y}_{\text{slip},i}}\right)^{2}}{2} {γ_{w}}+\left({\symbf{y}_{\text{wt},i}}-{\symbf{y}_{\text{slope},i}}\right)^{2} {γ_{w}}, & {\symbf{y}_{\text{wt},i}}\geq{}{\symbf{y}_{\text{slope},i}}\\ \frac{\left({\symbf{y}_{\text{wt},i}}-{\symbf{y}_{\text{slip},i}}\right)^{2}}{2} {γ_{w}}, & {\symbf{y}_{\text{slope},i}}\gt{}{\symbf{y}_{\text{wt},i}}\land{}{\symbf{y}_{\text{wt},i}}\gt{}{\symbf{y}_{\text{slip},i}}\\ 0, & {\symbf{y}_{\text{wt},i}}\leq{}{\symbf{y}_{\text{slip},i}} \end{cases}\]
Description
  • H is the interslice normal water forces (\(\frac{\text{N}}{\text{m}}\))
  • yslope is the y-coordinates of the slope (m)
  • i is the index (Unitless)
  • yslip is the y-coordinates of the slip surface (m)
  • γw is the unit weight of water (\(\frac{\text{N}}{\text{m}^{3}}\))
  • ywt is the y-coordinates of the water table (m)
Source

fredlund1977

RefBy

GD:resShearWO, IM:nrmShrForNum, and GD:mobShearWO

Refname DD:angleA
Label

Base angles

Symbol

α

Units

°

Equation \[\symbf{α}=\arctan\left(\frac{{\symbf{y}_{\text{slip},i}}-{\symbf{y}_{\text{slip},i-1}}}{{\symbf{x}_{\text{slip},i}}-{\symbf{x}_{\text{slip},i-1}}}\right)\]
Description
  • α is the base angles (°)
  • yslip is the y-coordinates of the slip surface (m)
  • i is the index (Unitless)
  • xslip is the x-coordinates of the slip surface (m)
Notes

This equation is based on the assumption that the base of a slice is a straight line (A:Surface-Base-Slice-between-Interslice-Straight-Lines).

Source

fredlund1977

RefBy

GD:resShearWO, IM:nrmShrForNum, GD:normForcEq, GD:momentEql, GD:mobShearWO, GD:bsShrFEq, DD:lengthLb, DD:convertFunc2, and DD:convertFunc1

Refname DD:angleB
Label

Surface angles

Symbol

β

Units

°

Equation \[\symbf{β}=\arctan\left(\frac{{\symbf{y}_{\text{slope},i}}-{\symbf{y}_{\text{slope},i-1}}}{{\symbf{x}_{\text{slope},i}}-{\symbf{x}_{\text{slope},i-1}}}\right)\]
Description
  • β is the surface angles (°)
  • yslope is the y-coordinates of the slope (m)
  • i is the index (Unitless)
  • xslope is the x-coordinates of the slope (m)
Notes

This equation is based on the assumption that the surface of a slice is a straight line (A:Surface-Base-Slice-between-Interslice-Straight-Lines).

Source

fredlund1977

RefBy

GD:resShearWO, IM:nrmShrForNum, GD:normForcEq, GD:momentEql, GD:mobShearWO, GD:bsShrFEq, and DD:lengthLs

Refname DD:lengthB
Label

Base width of slices

Symbol

b

Units

m

Equation \[\symbf{b}={\symbf{x}_{\text{slip},i}}-{\symbf{x}_{\text{slip},i-1}}\]
Description
  • b is the base width of slices (m)
  • xslip is the x-coordinates of the slip surface (m)
  • i is the index (Unitless)
Source

fredlund1977

RefBy

GD:sliceWght, IM:nrmShrForNum, IM:nrmShrForDen, GD:momentEql, DD:lengthLs, and DD:lengthLb

Refname DD:lengthLb
Label

Total base lengths of slices

Symbol

Lb

Units

m

Equation \[{\symbf{L}_{b}}={\symbf{b}}_{i} \sec\left({\symbf{α}}_{i}\right)\]
Description
  • Lb is the total base lengths of slices (m)
  • b is the base width of slices (m)
  • i is the index (Unitless)
  • α is the base angles (°)
Notes

b is defined in DD:lengthB and α is defined in DD:angleA.

Source

fredlund1977

RefBy

GD:resShr, GD:resShearWO, GD:mobShr, and GD:baseWtrF

Refname DD:lengthLs
Label

Surface lengths of slices

Symbol

Ls

Units

m

Equation \[{\symbf{L}_{s}}={\symbf{b}}_{i} \sec\left({\symbf{β}}_{i}\right)\]
Description
  • Ls is the surface lengths of slices (m)
  • b is the base width of slices (m)
  • i is the index (Unitless)
  • β is the surface angles (°)
Notes

b is defined in DD:lengthB and β is defined in DD:angleB.

Source

fredlund1977

RefBy

GD:srfWtrF

Refname DD:slcHeight
Label

y-direction heights of slices

Symbol

h

Units

m

Equation \[\symbf{h}=\frac{1}{2} \left({{\symbf{h}^{\text{R}}}}_{i}+{{\symbf{h}^{\text{L}}}}_{i}\right)\]
Description
  • h is the y-direction heights of slices (m)
  • hR is the heights of the right side of slices (m)
  • i is the index (Unitless)
  • hL is the heights of the left side of slices (m)
Notes

This equation is based on the assumption that the surface and base of a slice are straight lines (A:Surface-Base-Slice-between-Interslice-Straight-Lines).

hR and hL are defined in DD:sliceHghtRightDD and DD:sliceHghtLeftDD, respectively.

Source

fredlund1977

RefBy

IM:nrmShrForNum and GD:momentEql

Refname DD:normStress
Label

Total normal stress

Symbol

σ

Units

Pa

Equation \[σ=\frac{{F_{\text{n}}}}{A}\]
Description
  • σ is the total normal stress (Pa)
  • Fn is the total normal force (N)
  • A is the area (m2)
Source

huston2008

RefBy

GD:resShr, GD:effNormF, and TM:effStress

Refname DD:tangStress
Label

Tangential stress

Symbol

τ

Units

Pa

Equation \[τ=\frac{{F_{\text{t}}}}{A}\]
Description
  • τ is the tangential stress (Pa)
  • Ft is the tangential force (N)
  • A is the area (m2)
Source

huston2008

RefBy

GD:resShr

Refname DD:torque
Label

Torque

Symbol

τ

Units

Nm

Equation \[\symbf{τ}=\symbf{r}\times\symbf{F}\]
Description
  • τ is the torque (Nm)
  • r is the position vector (m)
  • F is the force (N)
Notes

The torque on a body measures the tendency of a force to rotate the body around an axis or pivot.

Source

--

RefBy

GD:momentEql

Refname DD:ratioVariation
Label

Interslice normal to shear force ratio variation function

Symbol

f

Units

Unitless

Equation \[\symbf{f}=\begin{cases} 1, & \mathit{const\_f}\\ \sin\left(π \frac{{\symbf{x}_{\text{slip},i}}-{\symbf{x}_{\text{slip},0}}}{{\symbf{x}_{\text{slip},n}}-{\symbf{x}_{\text{slip},0}}}\right), & \neg{}\mathit{const\_f} \end{cases}\]
Description
  • f is the interslice normal to shear force ratio variation function (Unitless)
  • π is the ratio of circumference to diameter for any circle (Unitless)
  • xslip is the x-coordinates of the slip surface (m)
  • i is the index (Unitless)
  • n is the number of slices (Unitless)
  • const_f is the decision on f (Unitless)
Source

fredlund1977

RefBy

IM:nrmShrForDen, GD:normShrR, DD:convertFunc2, and DD:convertFunc1

Refname DD:convertFunc1
Label

First function for incorporating interslice forces into shear force

Symbol

Φ

Units

Unitless

Equation \[\symbf{Φ}=\left(λ {\symbf{f}}_{i} \cos\left({\symbf{α}}_{i}\right)-\sin\left({\symbf{α}}_{i}\right)\right) \tan\left(φ'\right)-\left(λ {\symbf{f}}_{i} \sin\left({\symbf{α}}_{i}\right)+\cos\left({\symbf{α}}_{i}\right)\right) {F_{\text{S}}}\]
Description
  • Φ is the first function for incorporating interslice forces into shear force (Unitless)
  • λ is the proportionality constant (Unitless)
  • f is the interslice normal to shear force ratio variation function (Unitless)
  • i is the index (Unitless)
  • α is the base angles (°)
  • φ′ is the effective angle of friction (°)
  • FS is the factor of safety (Unitless)
Notes

f is defined in DD:ratioVariation and α is defined in DD:angleA.

Source

chen2005 and karchewski2012

RefBy

IM:intsliceFs, IM:fctSfty, and DD:convertFunc2

Refname DD:convertFunc2
Label

Second function for incorporating interslice forces into shear force

Symbol

Ψ

Units

Unitless

Equation \[\symbf{Ψ}=\frac{\left(λ {\symbf{f}}_{i} \cos\left({\symbf{α}}_{i}\right)-\sin\left({\symbf{α}}_{i}\right)\right) \tan\left(φ'\right)-\left(λ {\symbf{f}}_{i} \sin\left({\symbf{α}}_{i}\right)+\cos\left({\symbf{α}}_{i}\right)\right) {F_{\text{S}}}}{{\symbf{Φ}}_{i-1}}\]
Description
  • Ψ is the second function for incorporating interslice forces into shear force (Unitless)
  • λ is the proportionality constant (Unitless)
  • f is the interslice normal to shear force ratio variation function (Unitless)
  • i is the index (Unitless)
  • α is the base angles (°)
  • φ′ is the effective angle of friction (°)
  • FS is the factor of safety (Unitless)
  • Φ is the first function for incorporating interslice forces into shear force (Unitless)
Notes

f is defined in DD:ratioVariation, α is defined in DD:angleA, and Φ is defined in DD:convertFunc1.

Source

chen2005 and karchewski2012

RefBy

IM:intsliceFs and IM:fctSfty

Refname DD:nrmForceSumDD
Label

Sums of the interslice normal forces

Symbol

FxG

Units

N

Equation \[{{\symbf{F}_{\text{x}}}^{\text{G}}}={\symbf{G}}_{i}+{\symbf{G}}_{i-1}\]
Description
  • FxG is the sums of the interslice normal forces (N)
  • G is the interslice normal forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
Source

fredlund1977

RefBy

Refname DD:watForceSumDD
Label

Sums of the interslice normal water forces

Symbol

FxH

Units

N

Equation \[{{\symbf{F}_{\text{x}}}^{\text{H}}}={\symbf{H}}_{i}+{\symbf{H}}_{i-1}\]
Description
  • FxH is the sums of the interslice normal water forces (N)
  • H is the interslice normal water forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
Source

fredlund1977

RefBy

Refname DD:sliceHghtRightDD
Label

Heights of the right side of slices

Symbol

hR

Units

m

Equation \[{\symbf{h}^{\text{R}}}={\symbf{y}_{\text{slope},i}}-{\symbf{y}_{\text{slip},i}}\]
Description
  • hR is the heights of the right side of slices (m)
  • yslope is the y-coordinates of the slope (m)
  • i is the index (Unitless)
  • yslip is the y-coordinates of the slip surface (m)
Source

fredlund1977

RefBy

DD:slcHeight

Refname DD:sliceHghtLeftDD
Label

Heights of the left side of slices

Symbol

hL

Units

m

Equation \[{\symbf{h}^{\text{L}}}={\symbf{y}_{\text{slope},i-1}}-{\symbf{y}_{\text{slip},i-1}}\]
Description
  • hL is the heights of the left side of slices (m)
  • yslope is the y-coordinates of the slope (m)
  • i is the index (Unitless)
  • yslip is the y-coordinates of the slip surface (m)
Source

fredlund1977

RefBy

DD:slcHeight

Instance Models

This section transforms the problem defined in the problem description into one which is expressed in mathematical terms. It uses concrete symbols defined in the data definitions to replace the abstract symbols in the models identified in theoretical models and general definitions.

The goals GS:Identify-Crit-and-FS, GS:Determine-Normal-Forces, and GS:Determine-Shear-Forces are met by the simultaneous solution of IM:fctSfty, IM:nrmShrFor, and IM:intsliceFs. The goal GS:Identify-Crit-and-FS is also contributed to by IM:crtSlpId.

The Morgenstern-Price method is a vertical slice, limit equilibrium slope stability analysis method. Analysis is performed by breaking the assumed slip surface into a series of vertical slices of mass. Static equilibrium analysis is performed, using two force equations and one moment equation as in TM:equilibrium. The problem is statically indeterminate with only these 3 equations and one constitutive equation (the Mohr Coulomb shear strength of TM:mcShrStrgth) so the assumption GD:normShrR and corresponding equation GD:normShrR are used. The force equilibrium equations can be modified to be expressed only in terms of known physical values, as done in GD:resShearWO and GD:mobShearWO.

Refname IM:fctSfty
Label

Factor of safety

Input

xslope, yslope, ywt, c′, φ′, γdry, γsat, γw, xslip, yslip, const_f

Output

FS

Input Constraints
Output Constraints
Equation \[{F_{\text{S}}}=\frac{\displaystyle\sum_{i=1}^{n-1}{{\symbf{R}}_{i} \displaystyle\prod_{v=i}^{n-1}{{\symbf{Ψ}}_{v}}}+{\symbf{R}}_{n}}{\displaystyle\sum_{i=1}^{n-1}{{\symbf{T}}_{i} \displaystyle\prod_{v=i}^{n-1}{{\symbf{Ψ}}_{v}}}+{\symbf{T}}_{n}}\]
Description
  • FS is the factor of safety (Unitless)
  • R is the resistive shear forces without the influence of interslice forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • Ψ is the second function for incorporating interslice forces into shear force (Unitless)
  • v is the local index (Unitless)
  • n is the number of slices (Unitless)
  • T is the mobilized shear forces without the influence of interslice forces (\(\frac{\text{N}}{\text{m}}\))
Notes

R is defined in GD:resShearWO, Ψ is defined in DD:convertFunc2, and T is defined in GD:mobShearWO

Source

chen2005 and karchewski2012

RefBy

IM:intsliceFs, IM:fctSfty, IM:nrmShrFor, FR:Display-Interslice-Shear-Forces, FR:Display-Interslice-Normal-Forces, FR:Display-Factor-of-Safety, and FR:Determine-Critical-Slip-Surface

The mobilized shear force defined in GD:bsShrFEq can be substituted into the definition of mobilized shear force based on the factor of safety, from GD:mobShr yielding Equation (1) below:

\[\left({\symbf{W}}_{i}-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \cos\left({\symbf{ω}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)-\left(-{K_{\text{c}}} {\symbf{W}}_{i}-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)=\frac{{\symbf{N'}}_{i} \tan\left(φ'\right)+c' {\symbf{L}_{b,i}}}{{F_{\text{S}}}}\]

An expression for the effective normal forces, N′, can be derived by substituting the normal forces equilibrium from GD:normForcEq into the definition for effective normal forces from GD:resShearWO. This results in Equation (2):

\[{\symbf{N'}}_{i}=\left({\symbf{W}}_{i}-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \cos\left({\symbf{ω}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)+\left(-{K_{\text{c}}} {\symbf{W}}_{i}-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)-{\symbf{U}_{\text{b},i}}\]

Substituting Equation (2) into Equation (1) gives:

\[\left({\symbf{W}}_{i}-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \cos\left({\symbf{ω}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)-\left(-{K_{\text{c}}} {\symbf{W}}_{i}-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)=\frac{\left(\left({\symbf{W}}_{i}-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \cos\left({\symbf{ω}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)+\left(-{K_{\text{c}}} {\symbf{W}}_{i}-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)-{\symbf{U}_{\text{b},i}}\right) \tan\left(φ'\right)+c' {\symbf{L}_{b,i}}}{{F_{\text{S}}}}\]

Since the interslice shear forces X and interslice normal forces G are unknown, they are separated from the other terms as follows:

\[\left({\symbf{W}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \cos\left({\symbf{ω}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)-\left(-{K_{\text{c}}} {\symbf{W}}_{i}-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)-\left(-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}\right) \cos\left({\symbf{α}}_{i}\right)+\left(-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}\right) \sin\left({\symbf{α}}_{i}\right)=\frac{\left(\left({\symbf{W}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \cos\left({\symbf{ω}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)+\left(-{K_{\text{c}}} {\symbf{W}}_{i}-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)+\left(-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}\right) \sin\left({\symbf{α}}_{i}\right)+\left(-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}\right) \cos\left({\symbf{α}}_{i}\right)-{\symbf{U}_{\text{b},i}}\right) \tan\left(φ'\right)+c' {\symbf{L}_{b,i}}}{{F_{\text{S}}}}\]

Applying assumptions A:Seismic-Force and A:Surface-Load, which state that the seismic coefficient and the external forces, respectively, are zero, allows for further simplification as shown below:

\[\left({\symbf{W}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)-\left(-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)-\left(-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}\right) \cos\left({\symbf{α}}_{i}\right)+\left(-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}\right) \sin\left({\symbf{α}}_{i}\right)=\frac{\left(\left({\symbf{W}}_{i}+{\symbf{U}_{\text{g},i}} \cos\left({\symbf{β}}_{i}\right)\right) \cos\left({\symbf{α}}_{i}\right)+\left(-{\symbf{H}}_{i}+{\symbf{H}}_{i-1}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)\right) \sin\left({\symbf{α}}_{i}\right)+\left(-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}\right) \sin\left({\symbf{α}}_{i}\right)+\left(-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}\right) \cos\left({\symbf{α}}_{i}\right)-{\symbf{U}_{\text{b},i}}\right) \tan\left(φ'\right)+c' {\symbf{L}_{b,i}}}{{F_{\text{S}}}}\]

The definitions of GD:resShearWO and GD:mobShearWO are present in this equation, and thus can be replaced by Ri and Ti, respectively:

\[{\symbf{T}}_{i}+\left(-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}\right) \sin\left({\symbf{α}}_{i}\right)-\left(-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}\right) \cos\left({\symbf{α}}_{i}\right)=\frac{{\symbf{R}}_{i}+\left(\left(-{\symbf{X}}_{i-1}+{\symbf{X}}_{i}\right) \cos\left({\symbf{α}}_{i}\right)+\left(-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}\right) \sin\left({\symbf{α}}_{i}\right)\right) \tan\left(φ'\right)}{{F_{\text{S}}}}\]

The interslice shear forces X can be expressed in terms of the interslice normal forces G using A:Interslice-Norm-Shear-Forces-Linear and GD:normShrR, resulting in:

\[{\symbf{T}}_{i}+\left(-λ {\symbf{f}}_{i-1} {\symbf{G}}_{i-1}+λ {\symbf{f}}_{i} {\symbf{G}}_{i}\right) \sin\left({\symbf{α}}_{i}\right)-\left(-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}\right) \cos\left({\symbf{α}}_{i}\right)=\frac{{\symbf{R}}_{i}+\left(\left(-λ {\symbf{f}}_{i-1} {\symbf{G}}_{i-1}+λ {\symbf{f}}_{i} {\symbf{G}}_{i}\right) \cos\left({\symbf{α}}_{i}\right)+\left(-{\symbf{G}}_{i}+{\symbf{G}}_{i-1}\right) \sin\left({\symbf{α}}_{i}\right)\right) \tan\left(φ'\right)}{{F_{\text{S}}}}\]

Rearranging yields the following:

\[{\symbf{G}}_{i} \left(\left(λ {\symbf{f}}_{i} \cos\left({\symbf{α}}_{i}\right)-\sin\left({\symbf{α}}_{i}\right)\right) \tan\left(φ'\right)-\left(λ {\symbf{f}}_{i} \sin\left({\symbf{α}}_{i}\right)+\cos\left({\symbf{α}}_{i}\right)\right) {F_{\text{S}}}\right)={\symbf{G}}_{i-1} \left(\left(λ {\symbf{f}}_{i-1} \cos\left({\symbf{α}}_{i}\right)-\sin\left({\symbf{α}}_{i}\right)\right) \tan\left(φ'\right)-\left(λ {\symbf{f}}_{i-1} \sin\left({\symbf{α}}_{i}\right)+\cos\left({\symbf{α}}_{i}\right)\right) {F_{\text{S}}}\right)+{F_{\text{S}}} {\symbf{T}}_{i}-{\symbf{R}}_{i}\]

The definitions for Φ and Ψ from DD:convertFunc1 and DD:convertFunc2 simplify the above to Equation (3):

\[{\symbf{G}}_{i} {\symbf{Φ}}_{i}={\symbf{Ψ}}_{i-1} {\symbf{G}}_{i-1} {\symbf{Φ}}_{i-1}+{F_{\text{S}}} {\symbf{T}}_{i}-{\symbf{R}}_{i}\]

Versions of Equation (3) instantiated for slices 1 to n are shown below:

\[{\symbf{G}}_{1} {\symbf{Φ}}_{1}={\symbf{Ψ}}_{0} {\symbf{G}}_{0} {\symbf{Φ}}_{0}+{F_{\text{S}}} {\symbf{T}}_{1}-{\symbf{R}}_{1}\] \[{\symbf{G}}_{2} {\symbf{Φ}}_{2}={\symbf{Ψ}}_{1} {\symbf{G}}_{1} {\symbf{Φ}}_{1}+{F_{\text{S}}} {\symbf{T}}_{2}-{\symbf{R}}_{2}\] \[{\symbf{G}}_{3} {\symbf{Φ}}_{3}={\symbf{Ψ}}_{2} {\symbf{G}}_{2} {\symbf{Φ}}_{2}+{F_{\text{S}}} {\symbf{T}}_{3}-{\symbf{R}}_{3}\]

...

\[{\symbf{G}}_{n-2} {\symbf{Φ}}_{n-2}={\symbf{Ψ}}_{n-3} {\symbf{G}}_{n-3} {\symbf{Φ}}_{n-3}+{F_{\text{S}}} {\symbf{T}}_{n-2}-{\symbf{R}}_{n-2}\] \[{\symbf{G}}_{n-1} {\symbf{Φ}}_{n-1}={\symbf{Ψ}}_{n-2} {\symbf{G}}_{n-2} {\symbf{Φ}}_{n-2}+{F_{\text{S}}} {\symbf{T}}_{n-1}-{\symbf{R}}_{n-1}\] \[{\symbf{G}}_{n} {\symbf{Φ}}_{n}={\symbf{Ψ}}_{n-1} {\symbf{G}}_{n-1} {\symbf{Φ}}_{n-1}+{F_{\text{S}}} {\symbf{T}}_{n}-{\symbf{R}}_{n}\]

Applying A:Edge-Slices, which says that G0 and Gn are zero, results in the following special cases: Equation (8) for the first slice:

\[{\symbf{G}}_{1} {\symbf{Φ}}_{1}={F_{\text{S}}} {\symbf{T}}_{1}-{\symbf{R}}_{1}\]

and Equation (9) for the nth slice:

\[-\left(\frac{{F_{\text{S}}} {\symbf{T}}_{n}-{\symbf{R}}_{n}}{{\symbf{Ψ}}_{n-1}}\right)={\symbf{G}}_{n-1} {\symbf{Φ}}_{n-1}\]

Substituting Equation (8) into Equation (4) yields Equation (10):

\[{\symbf{G}}_{2} {\symbf{Φ}}_{2}={\symbf{Ψ}}_{1} \left({F_{\text{S}}} {\symbf{T}}_{1}-{\symbf{R}}_{1}\right)+{F_{\text{S}}} {\symbf{T}}_{2}-{\symbf{R}}_{2}\]

which can be substituted into Equation (5) to get Equation (11):

\[{\symbf{G}}_{3} {\symbf{Φ}}_{3}={\symbf{Ψ}}_{2} \left({\symbf{Ψ}}_{1} \left({F_{\text{S}}} {\symbf{T}}_{1}-{\symbf{R}}_{1}\right)+{F_{\text{S}}} {\symbf{T}}_{2}-{\symbf{R}}_{2}\right)+{F_{\text{S}}} {\symbf{T}}_{3}-{\symbf{R}}_{3}\]

and so on until Equation (12) is obtained from Equation (7):

\[{\symbf{G}}_{n-1} {\symbf{Φ}}_{n-1}={\symbf{Ψ}}_{n-2} \left({\symbf{Ψ}}_{n-3} \left({\symbf{Ψ}}_{1} \left({F_{\text{S}}} {\symbf{T}}_{1}-{\symbf{R}}_{1}\right)+{F_{\text{S}}} {\symbf{T}}_{2}-{\symbf{R}}_{2}\right)+{F_{\text{S}}} {\symbf{T}}_{n-2}-{\symbf{R}}_{n-2}\right)+{F_{\text{S}}} {\symbf{T}}_{n-1}-{\symbf{R}}_{n-1}\]

Equation (9) can then be substituted into the left-hand side of Equation (12), resulting in:

\[-\left(\frac{{F_{\text{S}}} {\symbf{T}}_{n}-{\symbf{R}}_{n}}{{\symbf{Ψ}}_{n-1}}\right)={\symbf{Ψ}}_{n-2} \left({\symbf{Ψ}}_{n-3} \left({\symbf{Ψ}}_{1} \left({F_{\text{S}}} {\symbf{T}}_{1}-{\symbf{R}}_{1}\right)+{F_{\text{S}}} {\symbf{T}}_{2}-{\symbf{R}}_{2}\right)+{F_{\text{S}}} {\symbf{T}}_{n-2}-{\symbf{R}}_{n-2}\right)+{F_{\text{S}}} {\symbf{T}}_{n-1}-{\symbf{R}}_{n-1}\]

This can be rearranged by multiplying both sides by Ψn−1 and then distributing the multiplication of each Ψ over addition to obtain:

\[-\left({F_{\text{S}}} {\symbf{T}}_{n}-{\symbf{R}}_{n}\right)={\symbf{Ψ}}_{n-1} {\symbf{Ψ}}_{n-2} {\symbf{Ψ}}_{1} \left({F_{\text{S}}} {\symbf{T}}_{1}-{\symbf{R}}_{1}\right)+{\symbf{Ψ}}_{n-1} {\symbf{Ψ}}_{n-2} {\symbf{Ψ}}_{2} \left({F_{\text{S}}} {\symbf{T}}_{2}-{\symbf{R}}_{2}\right)+{\symbf{Ψ}}_{n-1} \left({F_{\text{S}}} {\symbf{T}}_{n-1}-{\symbf{R}}_{n-1}\right)\]

The multiplication of the Ψ terms can be further distributed over the subtractions, resulting in the equation having terms that each either contain an R or a T. The equation can then be rearranged so terms containing an R are on one side of the equality, and terms containing a T are on the other. The multiplication by the factor of safety is common to all of the T terms, and thus can be factored out, resulting in:

\[{F_{\text{S}}} \left({\symbf{Ψ}}_{n-1} {\symbf{Ψ}}_{n-2} {\symbf{Ψ}}_{1} {\symbf{T}}_{1}+{\symbf{Ψ}}_{n-1} {\symbf{Ψ}}_{n-2} {\symbf{Ψ}}_{2} {\symbf{T}}_{2}+{\symbf{Ψ}}_{n-1} {\symbf{T}}_{n-1}+{\symbf{T}}_{n}\right)={\symbf{Ψ}}_{n-1} {\symbf{Ψ}}_{n-2} {\symbf{Ψ}}_{1} {\symbf{R}}_{1}+{\symbf{Ψ}}_{n-1} {\symbf{Ψ}}_{n-2} {\symbf{Ψ}}_{2} {\symbf{R}}_{2}+{\symbf{Ψ}}_{n-1} {\symbf{R}}_{n-1}+{\symbf{R}}_{n}\]

Isolating the factor of safety on the left-hand side and using compact notation for the products and sums yields Equation (13), which can also be seen in IM:fctSfty:

\[{F_{\text{S}}}=\frac{\displaystyle\sum_{i=1}^{n-1}{{\symbf{R}}_{i} \displaystyle\prod_{v=i}^{n-1}{{\symbf{Ψ}}_{v}}}+{\symbf{R}}_{n}}{\displaystyle\sum_{i=1}^{n-1}{{\symbf{T}}_{i} \displaystyle\prod_{v=i}^{n-1}{{\symbf{Ψ}}_{v}}}+{\symbf{T}}_{n}}\]

FS depends on the unknowns λ (IM:nrmShrFor) and G (IM:intsliceFs).

Refname IM:nrmShrFor
Label

Normal and shear force proportionality constant

Input

xslope, yslope, ywt, γw, xslip, yslip, const_f

Output

λ

Input Constraints
Output Constraints
Equation \[λ=\frac{\displaystyle\sum_{i=1}^{n}{{\symbf{C}_{\text{num},i}}}}{\displaystyle\sum_{i=1}^{n}{{\symbf{C}_{\text{den},i}}}}\]
Description
  • λ is the proportionality constant (Unitless)
  • Cnum is the proportionality constant numerator (N)
  • i is the index (Unitless)
  • Cden is the proportionality constant denominator (N)
Notes

Cnum is defined in IM:nrmShrForNum and Cden is defined in IM:nrmShrForDen.

Source

chen2005 and karchewski2012

RefBy

IM:nrmShrForNum, IM:nrmShrForDen, IM:intsliceFs, IM:fctSfty, IM:nrmShrFor, FR:Display-Interslice-Shear-Forces, FR:Display-Interslice-Normal-Forces, FR:Display-Factor-of-Safety, and FR:Determine-Critical-Slip-Surface

From the moment equilibrium of GD:momentEql with the primary assumption for the Morgenstern-Price method of A:Interslice-Norm-Shear-Forces-Linear and associated definition GD:normShrR, Equation (14) can be derived:

\[0=-{\symbf{G}}_{i} \left({\symbf{h}_{\text{z},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+{\symbf{G}}_{i-1} \left({\symbf{h}_{\text{z},i-1}}-\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)-{\symbf{H}}_{i} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+{\symbf{H}}_{i-1} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i-1}}-\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+λ \frac{{\symbf{b}}_{i}}{2} \left({\symbf{G}}_{i} {\symbf{f}}_{i}+{\symbf{G}}_{i-1} {\symbf{f}}_{i-1}\right)+\frac{-{K_{\text{c}}} {\symbf{W}}_{i} {\symbf{h}}_{i}}{2}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right) {\symbf{h}}_{i}+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right) {\symbf{h}}_{i}\]

Rearranging the equation in terms of λ leads to Equation (15):

\[λ=\frac{-{\symbf{G}}_{i} \left({\symbf{h}_{\text{z},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+{\symbf{G}}_{i-1} \left({\symbf{h}_{\text{z},i-1}}-\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)-{\symbf{H}}_{i} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+{\symbf{H}}_{i-1} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i-1}}-\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+\frac{-{K_{\text{c}}} {\symbf{W}}_{i} {\symbf{h}}_{i}}{2}+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right) {\symbf{h}}_{i}+{\symbf{Q}}_{i} \sin\left({\symbf{ω}}_{i}\right) {\symbf{h}}_{i}}{-\frac{{\symbf{b}}_{i}}{2} \left({\symbf{G}}_{i} {\symbf{f}}_{i}+{\symbf{G}}_{i-1} {\symbf{f}}_{i-1}\right)}\]

This equation can be simplified by applying assumptions A:Seismic-Force and A:Surface-Load, which state that the seismic and external forces, respectively, are zero:

\[λ=\frac{-{\symbf{G}}_{i} \left({\symbf{h}_{\text{z},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+{\symbf{G}}_{i-1} \left({\symbf{h}_{\text{z},i-1}}-\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)-{\symbf{H}}_{i} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i}}+\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+{\symbf{H}}_{i-1} \left(\frac{1}{3} {\symbf{h}_{\text{z,w},i-1}}-\frac{{\symbf{b}}_{i}}{2} \tan\left({\symbf{α}}_{i}\right)\right)+{\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right) {\symbf{h}}_{i}}{-\frac{{\symbf{b}}_{i}}{2} \left({\symbf{G}}_{i} {\symbf{f}}_{i}+{\symbf{G}}_{i-1} {\symbf{f}}_{i-1}\right)}\]

Taking the summation of all slices, and applying A:Edge-Slices to set G0, Gn, H0, and Hn equal to zero, a general equation for the proportionality constant λ is developed in Equation (16), which combines IM:nrmShrFor, IM:nrmShrForNum, and IM:nrmShrForDen:

\[λ=\frac{\displaystyle\sum_{i=1}^{n}{{\symbf{b}}_{i} \left({{\symbf{F}_{\text{x}}}^{\text{G}}}+{{\symbf{F}_{\text{x}}}^{\text{H}}}\right) \tan\left({\symbf{α}}_{i}\right)+{\symbf{h}}_{i} -2 {\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right)}}{\displaystyle\sum_{i=1}^{n}{{\symbf{b}}_{i} \left({\symbf{G}}_{i} {\symbf{f}}_{i}+{\symbf{G}}_{i-1} {\symbf{f}}_{i-1}\right)}}\]

Equation (16) for λ is a function of the unknown interslice normal forces G (IM:intsliceFs) which itself depends on the unknown factor of safety FS (IM:fctSfty).

Refname IM:nrmShrForNum
Label

Normal and shear force proportionality constant numerator

Input

xslope, yslope, ywt, γw, xslip, yslip

Output

Cnum

Input Constraints
Output Constraints
Equation \[{\symbf{C}_{\text{num},i}}=\begin{cases} {\symbf{b}}_{1} \left({\symbf{G}}_{1}+{\symbf{H}}_{1}\right) \tan\left({\symbf{α}}_{1}\right), & i=1\\ {\symbf{b}}_{i} \left({{\symbf{F}_{\text{x}}}^{\text{G}}}+{{\symbf{F}_{\text{x}}}^{\text{H}}}\right) \tan\left({\symbf{α}}_{i}\right)+\symbf{h} -2 {\symbf{U}_{\text{g},i}} \sin\left({\symbf{β}}_{i}\right), & 2\leq{}i\leq{}n-1\\ {\symbf{b}}_{n} \left({\symbf{G}}_{n-1}+{\symbf{H}}_{n-1}\right) \tan\left({\symbf{α}}_{n-1}\right), & i=n \end{cases}\]
Description
  • Cnum is the proportionality constant numerator (N)
  • i is the index (Unitless)
  • b is the base width of slices (m)
  • G is the interslice normal forces (\(\frac{\text{N}}{\text{m}}\))
  • H is the interslice normal water forces (\(\frac{\text{N}}{\text{m}}\))
  • α is the base angles (°)
  • FxG is the sums of the interslice normal forces (N)
  • FxH is the sums of the interslice normal water forces (N)
  • h is the y-direction heights of slices (m)
  • Ug is the surface hydrostatic forces (\(\frac{\text{N}}{\text{m}}\))
  • β is the surface angles (°)
  • n is the number of slices (Unitless)
Notes

b is defined in DD:lengthB, H is defined in DD:intersliceWtrF, α is defined in DD:angleA, h is defined in DD:slcHeight, Ug is defined in GD:srfWtrF, and β is defined in DD:angleB.

Source

chen2005 and karchewski2012

RefBy

IM:nrmShrFor

See IM:nrmShrFor for the derivation of Cnum.

Refname IM:nrmShrForDen
Label

Normal and shear force proportionality constant denominator

Input

xslip, const_f

Output

Cden

Input Constraints
Output Constraints
Equation \[{\symbf{C}_{\text{den},i}}=\begin{cases} {\symbf{b}}_{1} {\symbf{f}}_{1} {\symbf{G}}_{1}, & i=1\\ {\symbf{b}}_{i} \left({\symbf{f}}_{i} {\symbf{G}}_{i}+{\symbf{f}}_{i-1} {\symbf{G}}_{i-1}\right), & 2\leq{}i\leq{}n-1\\ {\symbf{b}}_{n} {\symbf{G}}_{n-1} {\symbf{f}}_{n-1}, & i=n \end{cases}\]
Description
  • Cden is the proportionality constant denominator (N)
  • i is the index (Unitless)
  • b is the base width of slices (m)
  • f is the interslice normal to shear force ratio variation function (Unitless)
  • G is the interslice normal forces (\(\frac{\text{N}}{\text{m}}\))
  • n is the number of slices (Unitless)
Notes

b is defined in DD:lengthB and f is defined in DD:ratioVariation.

Source

chen2005 and karchewski2012

RefBy

IM:nrmShrFor

See IM:nrmShrFor for the derivation of Cden.

Refname IM:intsliceFs
Label

Interslice normal forces

Input

xslope, yslope, ywt, c′, φ′, γdry, γsat, γw, xslip, yslip, const_f

Output

G

Input Constraints
Output Constraints
Equation \[{\symbf{G}}_{i}=\begin{cases} \frac{{F_{\text{S}}} {\symbf{T}}_{1}-{\symbf{R}}_{1}}{{\symbf{Φ}}_{1}}, & i=1\\ \frac{{\symbf{Ψ}}_{i-1} {\symbf{G}}_{i-1}+{F_{\text{S}}} {\symbf{T}}_{i}-{\symbf{R}}_{i}}{{\symbf{Φ}}_{i}}, & 2\leq{}i\leq{}n-1\\ 0, & i=0\lor{}i=n \end{cases}\]
Description
  • G is the interslice normal forces (\(\frac{\text{N}}{\text{m}}\))
  • i is the index (Unitless)
  • FS is the factor of safety (Unitless)
  • T is the mobilized shear forces without the influence of interslice forces (\(\frac{\text{N}}{\text{m}}\))
  • R is the resistive shear forces without the influence of interslice forces (\(\frac{\text{N}}{\text{m}}\))
  • Φ is the first function for incorporating interslice forces into shear force (Unitless)
  • Ψ is the second function for incorporating interslice forces into shear force (Unitless)
  • n is the number of slices (Unitless)
Notes

T is defined in GD:mobShearWO, R is defined in GD:resShearWO, Φ is defined in DD:convertFunc1, and Ψ is defined in DD:convertFunc2.

Source

chen2005

RefBy

IM:intsliceFs, IM:fctSfty, IM:nrmShrFor, FR:Display-Interslice-Shear-Forces, FR:Display-Interslice-Normal-Forces, FR:Display-Factor-of-Safety, and FR:Determine-Critical-Slip-Surface

This derivation is identical to the derivation for IM:fctSfty up until Equation (3) shown again below:

\[{\symbf{G}}_{i} {\symbf{Φ}}_{i}={\symbf{Ψ}}_{i-1} {\symbf{G}}_{i-1} {\symbf{Φ}}_{i-1}+{F_{\text{S}}} {\symbf{T}}_{i}-{\symbf{R}}_{i}\]

A simple rearrangement of Equation (3) leads to Equation (17), also seen in IM:intsliceFs:

\[{\symbf{G}}_{i}=\frac{{\symbf{Ψ}}_{i-1} {\symbf{G}}_{i-1}+{F_{\text{S}}} {\symbf{T}}_{i}-{\symbf{R}}_{i}}{{\symbf{Φ}}_{i}}\]

The cases shown in IM:intsliceFs for when i = 0, i = 1, or i = n are derived by applying A:Edge-Slices, which says that G0 and Gn are zero, to Equation (17). G depends on the unknowns FS (IM:fctSfty) and λ (IM:nrmShrFor).

Refname IM:crtSlpId
Label

Critical slip surface identification

Input

xslope, yslope, xwt, ywt, c′, φ′, γdry, γsat, γw, const_f

Output

FSmin

Input Constraints
Output Constraints
Equation \[{\symbf{x}_{\text{slope}}}={\symbf{y}_{\text{slope}}}\]
Description
  • xslope is the x-coordinates of the slope (m)
  • yslope is the y-coordinates of the slope (m)
Notes

The minimization function must enforce the constraints on the critical slip surface expressed in A:Slip-Surface-Concave and Sec:Properties of a Correct Solution. The sizes of xwt and ywt must be equal and not 1. The sizes of xslope and yslope must be equal and at least 2. The first and last xwt values must be equal to the first and last xslope values. xwt and xslope values must be monotonically increasing. xslipmaxExt, xslipmaxEtr, xslipminExt, and xslipminEtr must be between or equal to the minimum and maximum xslope values. yslipmax cannot be below the minimum yslope value. yslipmin cannot be above the maximum yslope value. All x values of xcs,ycs must be between xslipminEtr and xslipmaxExt. All y values of xcs,ycs must not be below yslipmin. For any given vertex in xcs,ycs the y value must not exceed the yslope value corresponding to the same x value. The first and last vertices in xcs,ycs must each be equal to one of the vertices formed by xslope and yslope. The slope between consecutive vertices must be always increasing as x increases. The internal angle between consecutive vertices in xcs,ycs must not be below 110 degrees.

Source

li2010

RefBy

FR:Display-Graph and FR:Determine-Critical-Slip-Surface

Data Constraints

The Data Constraints Table shows the data constraints on the input variables. The column for physical constraints gives the physical limitations on the range of values that can be taken by the variable. The uncertainty column provides an estimate of the confidence with which the physical quantities can be measured. This information would be part of the input if one were performing an uncertainty quantification exercise. The constraints are conservative to give the user of the model the flexibility to experiment with unusual situations. The column of typical values is intended to provide a feel for a common scenario.

Var Physical Constraints Typical Value Uncert.
c′ c′ > 0 10000 Pa 10%
xslipmaxEtr -- 20 m 10%
xslipmaxExt -- 100 m 10%
xslipminEtr -- 0 m 10%
xslipminExt -- 50 m 10%
xslope -- 0 m 10%
xwt -- 0 m 10%
yslipmax -- 30 m 10%
yslipmin -- 0 m 10%
yslope -- 0 m 10%
ywt -- 0 m 10%
γdry γdry > 0 20000 \(\frac{\text{N}}{\text{m}^{3}}\) 10%
γsat γsat > 0 20000 \(\frac{\text{N}}{\text{m}^{3}}\) 10%
γw γw > 0 9800 \(\frac{\text{N}}{\text{m}^{3}}\) 10%
φ′ 0 < φ′ < 90 25 ° 10%

Input Data Constraints

Properties of a Correct Solution

The Data Constraints Table shows the data constraints on the output variables. The column for physical constraints gives the physical limitations on the range of values that can be taken by the variable.

Var Physical Constraints
FS FS > 0

Output Data Constraints

Requirements

This section provides the functional requirements, the tasks and behaviours that the software is expected to complete, and the non-functional requirements, the qualities that the software is expected to exhibit.

Functional Requirements

This section provides the functional requirements, the tasks and behaviours that the software is expected to complete.

Read-and-Store: Read the inputs, shown in the table Required Inputs, and store the data.

Verify-Input: Verify that the input data lie within the physical constraints.

Determine-Critical-Slip-Surface: Determine the critical slip surface for the input slope, corresponding to the minimum factor of safety, by using IM:fctSfty, IM:nrmShrFor, and IM:intsliceFs to calculate the factor of safety for a slip surface and using IM:crtSlpId to find the slip surface that minimizes it.

Verify-Output: Verify that the minimum factor of safety and critical slip surface satisfy the physical constraints shown in Properties of a Correct Solution.

Display-Input: Display as output the user-supplied inputs listed in Tab:inputsToOutputTable.

Display-Graph: Display the critical slip surface of the 2D slope, as determined from IM:crtSlpId, graphically.

Display-Factor-of-Safety: Display the value of the factor of safety for the critical slip surface, as determined from IM:fctSfty, IM:nrmShrFor, and IM:intsliceFs.

Display-Interslice-Normal-Forces: Using IM:fctSfty, IM:nrmShrFor, and IM:intsliceFs, calculate and graphically display the interslice normal forces.

Display-Interslice-Shear-Forces: Using IM:fctSfty, IM:nrmShrFor, and IM:intsliceFs, calculate and graphically display the interslice shear forces.

Write-Results-To-File: Provide the option of writing the output result data, as given in FR:Display-Input, FR:Display-Graph, FR:Display-Factor-of-Safety, FR:Display-Interslice-Normal-Forces, and FR:Display-Interslice-Shear-Forces, to a file.

Symbol Description Units
(x,y) Cartesian position coordinates m
c′ Effective cohesion Pa
const_f Decision on f --
xslipmaxEtr Maximum entry x-coordinate m
xslipmaxExt Maximum exit x-coordinate m
xslipminEtr Minimum entry x-coordinate m
xslipminExt Minimum exit x-coordinate m
xslope x-coordinates of the slope m
xwt x-coordinates of the water table m
yslipmax Maximum y-coordinate m
yslipmin Minimum y-coordinate m
yslope y-coordinates of the slope m
ywt y-coordinates of the water table m
γdry Soil dry unit weight \(\frac{\text{N}}{\text{m}^{3}}\)
γsat Soil saturated unit weight \(\frac{\text{N}}{\text{m}^{3}}\)
γw Unit weight of water \(\frac{\text{N}}{\text{m}^{3}}\)
φ′ Effective angle of friction °

Required Inputs following FR:Read-and-Store

Symbol Name
const_f decision on f
xslipmaxExt maximum exit x-coordinate
xslipmaxEtr maximum entry x-coordinate
xslipminExt minimum exit x-coordinate
xslipminEtr minimum entry x-coordinate
yslipmax maximum y-coordinate
yslipmin minimum y-coordinate

Inputs to be Returned as Output following FR:Display-Input

Non-Functional Requirements

This section provides the non-functional requirements, the qualities that the software is expected to exhibit.

Correct: The outputs of the code have the properties described in Sec:Properties of a Correct Solution.

Understandable: The code is modularized with complete module guide and module interface specification.

Reusable: The code is modularized.

Maintainable: If a likely change is made to the finished software, it will take at most 10% of the original development time, assuming the same development resources are available.

Likely Changes

This section lists the likely changes to be made to the software.

Calculate-Inhomogeneous-Soil-Layers: A:Soil-Layer-Homogeneous - The system currently assumes the soil mass is homogeneous. In the future, calculations can be added for inconsistent soil properties throughout.

Calculate-Seismic-Force: A:Seismic-Force - The system currently assumes no seismic force. In the future, calculations can be added for the presence of seismic force.

Calculate-External-Force: A:Surface-Load - The system currently assumes no external forces. In the future, calculations can be added for an imposed surface load on the slope.

Unlikely Changes

This section lists the unlikely changes to be made to the software.

Normal-And-Shear-Linear-Only: Changes related to A:Interslice-Norm-Shear-Forces-Linear are not possible due to the dependency of the calculations on the linear relationship between interslice normal forces and interslice shear forces.

2D-Analysis-Only: A:Effective-Norm-Stress-Large allows for 2D analysis with these models only because stress along the z-direction is zero. These models do not take into account stress in the z-direction, and therefore cannot be used without manipulation to attempt three-dimensional analysis.

Traceability Matrices and Graphs

The purpose of the traceability matrices is to provide easy references on what has to be additionally modified if a certain component is changed. Every time a component is changed, the items in the column of that component that are marked with an "X" should be modified as well. Tab:TraceMatAvsA shows the dependencies of the assumptions on each other. Tab:TraceMatAvsAll shows the dependencies of the data definitions, theoretical models, general definitions, instance models, requirements, likely changes, and unlikely changes on the assumptions. Tab:TraceMatRefvsRef shows the dependencies of the data definitions, theoretical models, general definitions, and instance models on each other. Tab:TraceMatAllvsR shows the dependencies of the requirements and goal statements on the data definitions, theoretical models, general definitions, and instance models.

A:Slip-Surface-Concave A:Factor-of-Safety A:Soil-Layer-Homogeneous A:Soil-Properties A:Soil-Layers-Isotropic A:Interslice-Norm-Shear-Forces-Linear A:Plane-Strain-Conditions A:Effective-Norm-Stress-Large A:Surface-Base-Slice-between-Interslice-Straight-Lines A:Edge-Slices A:Seismic-Force A:Surface-Load A:Water-Intersects-Base-Edge A:Water-Intersects-Surface-Edge A:Negligible-Effect-Surface-Slope-Seismic A:Hydrostatic-Force-Slice-Midpoint
DD:intersliceWtrF
DD:angleA X
DD:angleB X
DD:lengthB
DD:lengthLb
DD:lengthLs
DD:slcHeight X
DD:normStress
DD:tangStress
DD:torque
DD:ratioVariation
DD:convertFunc1
DD:convertFunc2
DD:nrmForceSumDD
DD:watForceSumDD
DD:sliceHghtRightDD
DD:sliceHghtLeftDD
TM:factOfSafety
TM:equilibrium X
TM:mcShrStrgth X
TM:effStress
TM:NewtonSecLawMot
GD:normForcEq
GD:bsShrFEq
GD:resShr X X X X
GD:mobShr X
GD:effNormF X
GD:resShearWO
GD:mobShearWO
GD:normShrR X
GD:momentEql X X
GD:weight
GD:sliceWght X X X X X
GD:hsPressure
GD:baseWtrF X X X X
GD:srfWtrF X X X X
IM:fctSfty X X X X
IM:nrmShrFor X X X X
IM:nrmShrForNum
IM:nrmShrForDen
IM:intsliceFs X
IM:crtSlpId X
FR:Read-and-Store
FR:Verify-Input
FR:Determine-Critical-Slip-Surface
FR:Verify-Output
FR:Display-Input
FR:Display-Graph
FR:Display-Factor-of-Safety
FR:Display-Interslice-Normal-Forces
FR:Display-Interslice-Shear-Forces
FR:Write-Results-To-File
NFR:Correct
NFR:Understandable
NFR:Reusable
NFR:Maintainable
LC:Calculate-Inhomogeneous-Soil-Layers X
LC:Calculate-Seismic-Force X
LC:Calculate-External-Force X
UC:Normal-And-Shear-Linear-Only X
UC:2D-Analysis-Only X

Traceability Matrix Showing the Connections Between Assumptions and Other Items

The purpose of the traceability graphs is also to provide easy references on what has to be additionally modified if a certain component is changed. The arrows in the graphs represent dependencies. The component at the tail of an arrow is depended on by the component at the head of that arrow. Therefore, if a component is changed, the components that it points to should also be changed. Fig:TraceGraphAvsA shows the dependencies of assumptions on each other. Fig:TraceGraphAvsAll shows the dependencies of data definitions, theoretical models, general definitions, instance models, requirements, likely changes, and unlikely changes on the assumptions. Fig:TraceGraphRefvsRef shows the dependencies of data definitions, theoretical models, general definitions, and instance models on each other. Fig:TraceGraphAllvsR shows the dependencies of requirements and goal statements on the data definitions, theoretical models, general definitions, and instance models. Fig:TraceGraphAllvsAll shows the dependencies of dependencies of assumptions, models, definitions, requirements, goals, and changes with each other.

TraceGraphAvsA
TraceGraphAvsA
TraceGraphAvsAll
TraceGraphAvsAll
TraceGraphRefvsRef
TraceGraphRefvsRef
TraceGraphAllvsR
TraceGraphAllvsR
TraceGraphAllvsAll
TraceGraphAllvsAll

For convenience, the following graphs can be found at the links below:

Values of Auxiliary Constants

There are no auxiliary constants.

References