A reinforcement (an anchor, a nail, a pile, a geosynthetic sheet or a wall) acts on the sliding mass where it crosses the slip surface. The force depends on the support type (OGR Slip2D has nine) and enters each method as an active or a passive force on the base of the slice crossed. Two related calculations follow: the water thrust in a tension crack and the back analysis of the force a slope needs to reach a target factor of safety.
3.10.1A support on a slip surface#
A support is a straight segment from its head, the end at the slope face where its plate or connection is, to its tail, the anchored end. OGR Slip2D finds the first point where the segment crosses the slip surface, counting from the head, and evaluates there the force the support can deliver, , with the distance from the head to the crossing. The force acts at that point, on the base of the slice that contains it (Figure 3.10.1). Everything is per metre of slope: the capacity of a discrete element is divided by its out-of-plane spacing .
The force, in the direction its orientation gives it (Section 3.10.3), is resolved on the base of that slice, whose chord is inclined at :
where
- components of the support force, to the right and upwards (kN/m)
- sense of sliding, ±1, chosen so that is the unit vector along the base that opposes sliding
- unit normal pointing from the slice into the ground below its base
- tangential component, positive when it resists sliding (kN/m)
- normal component, positive when it presses the slice onto its base (kN/m)
A support that never reaches the surface, or that develops no capacity where it crosses it, puts no force on that surface. The factor of safety is then the one of the slope without it, which looks like any other result, so the analysis notes say how many of the supports placed contributed nothing and why.
3.10.2Active and passive supports#
Whether a reinforcement force exists before the mass moves, or develops only as it moves, decides how it enters a factor of safety defined on strength. Following Duncan & Wright (2005), OGR Slip2D offers both readings. An active support, such as a tensioned anchor, is a force already acting on the mass, and it reduces the driving force. A passive support, such as an untensioned nail, a pile or a geosynthetic, is mobilised with the soil, so its tangential force is added to the resistance and divided by the factor of safety like the shear strength. With and the resisting and driving terms of a method:
where
- resisting and driving terms of the method without the support (kN/m, or moments over the radius on a circle)
- tangential and normal components of Eq. (3.10.1), summed over the supports and slices crossed
- friction angle at the base of the slice crossed
In both expressions the normal component raises the frictional resistance of the base it presses; only the tangential component changes side. Written for one base, the shear mobilised on a base crossed by a support becomes
with the base normal force, which already carries through the equilibrium of the slice. The two readings differ only in that line: summed over a mass in force equilibrium, it turns into . Each family of methods applies it in its own equations (Section 3.2):
- Spencer, GLE / Morgenstern–Price, the Corps of Engineers methods and Lowe–Karafiath take the normal part of the force as a load on the slice and add the tangential part to the mobilised shear of the base, at if active and if passive, as in Eq. (3.10.4).
- Bishop's and Janbu's simplified methods obtain the base normal force from the vertical equilibrium of the slice, so the vertical component of the support force enters that equilibrium and its frictional share is divided by (or ) like the weight. Janbu's methods, which balance shear along the base, take on the driving side.
- The ordinary method resolves forces normal to the base, so adds directly to the normal force.
- On a circle, the methods with a moment equation take the support's moment about the centre, divided by the radius, in place of ; a support tangent to the slip surface takes the arm , like a base shear.
On a planar slip surface the mass is a single free body and the two readings have closed forms,
with the weight of the wedge and the inclination and length of the plane. Every method that closes global force equilibrium must give these values; the program's source records that the ordinary method, Janbu's simplified method, the two Corps of Engineers methods and Lowe–Karafiath reproduce them.
Force Application (Active or Passive) and Force Orientation are set for each support type in Properties → Define Support; a support placed without them takes the defaults of its type (Table 3.10.1). In the ordinary, Bishop and Janbu methods a surface whose active support force reaches or exceeds the driving term is marked inadmissible and is not reported as critical: as approaches the active expression grows without bound and ceases to have physical meaning.
Source: ogr_slip2d/support_integration.py · ogr_slip2d/interslice.py · ogr_slip2d/methods/bishop.py3.10.3Orientation of the force#
The direction of at the crossing is one of six options:
- Tangent to slip surface: along the base of the slice crossed, in the sense that resists sliding.
- Parallel to support: along the support axis from head to tail, the direction in which a bolt in tension pulls the mass. This option depends only on the geometry: a support whose anchor lies downhill of the slip surface pushes the mass, and the program applies it as drawn.
- Bisector: the bisector of the two directions above.
- Horizontal: horizontal, against the direction of failure.
- Perpendicular to pile: perpendicular to the support axis, on the side that resists sliding.
- User-defined angle: at a given angle from the positive axis.
The resisting senses follow the failure direction set for the model (Section 3.1.5). Four types (the grouted tieback, the grouted tieback with friction, the soil nail and the helical anchor) also declare a shear capacity. When it is not zero, their shear capacity divided by the spacing is added as a second vector, perpendicular to the support axis on the side that opposes sliding, and the resultant replaces the axial force in Eq. (3.10.1).
Head and tail are not interchangeable: the stripping length is measured from the head, the pullout length from the tail, and a parallel force points from head to tail. A support drawn with its head in the stable ground and its tail inside the sliding mass is reported in the analysis notes, because drawn that way round all three are reversed and the factor of safety does not show it.
3.10.4Support types#
OGR Slip2D has nine support types. Table 3.10.1 lists them with their default orientation and application and the force each one delivers at the crossing.
| Type | Orientation | Application | Force at the crossing |
|---|---|---|---|
| End anchored | parallel to support | active | Constant: anchor capacity / |
| Grouted tieback | parallel to support | active | Smallest of pullout, tensile and stripping; bond strength per metre over a bonded length at the tail |
| Grouted tieback with friction | parallel to support | active | As the grouted tieback, with the bond from an interface law integrated around the grout |
| Soil nail | parallel to support | passive | As the grouted tieback, bonded over its whole length |
| Pile / micropile | tangent to slip surface | passive | Shear: pile shear strength / ; or the lateral thrust of Ito & Matsui (1975) |
| Geosynthetic | parallel to support | passive | Smallest of pullout, tensile and stripping, both faces of the sheet, times the coverage |
| User defined | parallel to support | active | Table of force against distance from the head, / |
| Retaining wall (EFP) | horizontal | active | Pressure profile integrated from the crest down to the crossing; no spacing |
| Helical anchor | tangent to slip surface | active | Smallest of seven capacities: three in pullout, tensile, three in stripping |
Capacity along a support#
For the types with failure modes, the force available at the crossing is the smallest of the capacities that apply there, so it varies with the position of the crossing along the support. For a grouted tieback, with a bond strength per metre over a bonded length at the tail end:
where
- bonded length beyond the crossing, towards the tail: the whole bond when the crossing is in the free length (m)
- bonded length between the head and the crossing: zero when the crossing is in the free length (m)
- tensile capacity of the tendon and capacity of the plate at the head (kN)
The three terms are the pullout of the anchored length, the failure of the tendon and the stripping of the moving mass off the plate together with the bond on the head side. A soil nail is the same with the bond over its whole length, and . In the grouted tieback with friction the bond is the integral of an interface strength over the surface of a grout cylinder of diameter ,
and in a geosynthetic, which mobilises both faces of the sheet (Berg, Christopher & Samtani, 2009), with the percentage of the width of slope its strips cover (100 % for a continuous sheet),
with the tensile capacity per metre of sheet, the strength of the connection at the face, from the face to the crossing and from the crossing to the embedded end. A sheet whose embedded end is anchored cannot be pulled out, and its first term is dropped.
Interface strength#
Where the bond depends on the stress in the ground, OGR Slip2D samples the support once per analysis in 50 equal segments, from head to tail, and evaluates the interface strength at the middle of each; the integrals above are taken over that profile. The normal stress on the interface is the effective vertical stress at the point: the weight of the soil column above it, plus any ponded water and the vertical pressure of the distributed loads over it, minus the pore pressure. It is the stress of the pullout equation of Berg, Christopher & Samtani (2009). Two envelopes are available:
where
- adhesion (kPa) and friction angle (°) of the interface, the Mohr–Coulomb bond of Jewell (1996)
- limiting strength at large normal stress (kPa) and tangent friction angle at zero normal stress (°), after Esterhuizen, Filz & Duncan (2001)
- effective vertical stress at the point of the support, not less than zero (kPa)
The same pair of numbers describes two different envelopes in the two forms: in the hyperbolic one the adhesion is the limit at large stress, not the strength at zero stress. A geosynthetic can also take as a fraction of the surrounding soil's own strength, evaluated with the material's strength model, or as with a pullout friction factor , constant or varying linearly with depth (Berg, Christopher & Samtani, 2009); the scale-effect factor of that guidance is entered multiplied into .
Helical anchors#
A helical anchor carries steel plates along its shaft, and its capacity steps as the crossing passes each plate. The bearing capacity of one plate is that of a circular footing without the self-weight term (Terzaghi, 1943), as used for helices by Perko (2009):
with the plate area less the area of the shaft, and the factors of Prandtl (1921) and Reissner (1924) ( as ), and the equivalent Mohr–Coulomb pair of the soil at the plate, linearised at the effective vertical stress. Pullout is the smallest of three mechanisms, evaluated for the plates beyond the crossing: shallow failure, in which a cylinder of soil is sheared from the crossing to the farthest plate; cylindrical shear between the plates plus the bearing of the nearest one (Mitsch & Clemence, 1985; Mooney, Adamczak & Clemence, 1985); and the sum of the individual bearing capacities. Stripping is the same three mechanisms for the plates between the head and the crossing, plus the capacity of the head assembly. The shear on the cylinders is the soil's own strength at the effective vertical stress, with no separate earth-pressure or adhesion coefficient.
At the plate bearing becomes times the undrained strength, whereas Mooney, Adamczak & Clemence (1985) use for deep plates in clay. The program keeps the form above because it is the one its published verification values were produced with.
Piles#
A pile resists across its axis. In shear mode it applies a constant force, the shear strength of its section divided by the spacing, wherever the slip surface crosses it. In Ito & Matsui mode the force comes from the soil squeezing between the piles of a row, which loads each pile with a lateral force per unit thickness of layer (Ito & Matsui, 1975):
with
where
- centre-to-centre spacing of the piles in the row (m)
- opening between neighbouring piles, with the pile diameter (m)
- equivalent Mohr–Coulomb parameters of the soil at that depth, from the linearisation the methods use
- effective vertical stress, which takes the place of the of the original paper (kPa)
For below rad the program evaluates the cohesive solution of the same paper, , which is the limit of Eq. (3.10.13) and avoids its cancellation in floating point. The force on the sliding mass, per metre of slope, is the integral of from the top of the pile down to the crossing, divided by the spacing, as Cai & Ugai (2000) write it:
with the depth of the crossing measured along the pile from its top; the integral, not the samples, is clipped at zero. Its resultant acts at the crossing or at the centroid of the diagram above it. The theory is written for rigid vertical piles in a dry layer, and the analysis notes say when a model steps outside it: piles that touch (, where the row applies no force), piles more than 2° off vertical (integrated along the shaft), pore pressure along the pile, and negative pressures near the surface (integrated as published).
A published discussion of the solution (Ukritchon & Keawsawasvong, 2017) argues that it underestimates the pile force, increasingly with the friction angle. OGR Slip2D keeps the original equations, which are the ones the published factors of safety it is checked against were computed with.
Retaining wall (equivalent fluid pressure)#
The wall type does not derive an earth pressure: it integrates a pressure profile given over its length, from its upper end down to the crossing,
with the equivalent fluid pressure (kN/m³) and measured from the crest of the wall, whichever end was drawn first (Terzaghi & Peck, 1967). Four profiles are offered: uniform; triangular, zero at the crest; trapezoidal, zero at both ends, rising to over of the length at each end and constant over the fraction between (60 % by default); and a custom table of pressure against relative distance from the crest, linearly interpolated. A wall is continuous, so its force is already per metre and it has no spacing; the force is horizontal and active by default and can act at the crossing or at the centroid of the diagram above it. A wall drawn with both ends at the same elevation has no crest to measure from and is left out of the analysis, with a note. The couple left by moving the resultant to the centroid only enters the four methods with a moment equation: ordinary, Bishop, Spencer and GLE.
Support types and their parameters are defined in Properties → Define Support and placed with Support → Add Support or Support → Add Support Pattern, head first. The pullout laws that depend on the stress are sampled once per analysis, not once per trial surface, because the stress along a support depends only on what lies above it. A model that cannot be priced at a support (for instance, a soil whose strength model cannot be read there) leaves that support out and says so, instead of silently analysing the slope without it.
Source: ogr_core/support/support.py · ogr_core/support/bond.py · ogr_core/support/ito_matsui.py · ogr_core/support/helical_anchor.py · ogr_core/support/retaining_wall.py3.10.5Tension cracks#
A tension crack is drawn as a boundary, and the zone above it is soil that cannot carry tension, so no shear strength may be counted along a slip surface inside it. When the crest end of a slip surface lies in that zone, OGR Slip2D cuts the surface at its first crossing with the crack boundary, counting from the crest, and closes the sliding mass with a vertical wall from that point up to the ground surface; the soil between the wall and the slip surface keeps its whole weight (Duncan & Wright, 2005). Only the crest end is cut, because the crack forms at the head of the slide, where the interslice forces go into tension, while the toe is in compression. A surface whose crest lies in the zone but never reaches the boundary lies entirely in soil that cannot resist, and it is discarded.
Water in the crack pushes on the wall with the hydrostatic resultant
where
- height of water in the crack, from the crack boundary at the wall up to the water level (m)
- unit weight of water, the pore-fluid unit weight of the project (9.81 kN/m³ by default)
The thrust (Terzaghi, 1943) is applied to the outermost slice, next to the wall, and pushes the mass away from the intact ground. The water level is set by one of seven options: dry; filled, which is the default; filled to a percentage of the crack height; filled up to a given elevation; filled to a level a given depth below the top of the crack; at the water table; or at a piezometric line. The level is always clipped to the height of the crack, and a dry crack still truncates the surface.
The crack is drawn with Boundaries → Add Tension Crack and its water is set in Properties → Define Tension Crack; one crack boundary per model is used. A separate mechanism, for circles only: a circle whose ground entry point lies above its own centre would overhang, and with the search option that creates tension cracks on such reverse-curvature circles (on by default) it is ended by a vertical crack where it turns vertical. That crack is a consequence of the search geometry and is always dry (Section 3.4).
Source: ogr_slip2d/slicer.py · ogr_core/geometry/tension_crack.py · ogr_slip2d/surface.py3.10.6Back analysis of the required support force#
In the back analysis the factor of safety is fixed at a target and the unknown is the horizontal force that would have to be added, at a chosen elevation, to reach it. With the factor fixed, the resisting and driving sums of the method can be evaluated directly, and Eq. (3.10.2) and Eq. (3.10.3) give the force in closed form, with no iteration:
where
- resisting sum evaluated at , and driving sum, of the method
- arm factor: in Bishop's simplified method on a circle, the moment arm of a horizontal force at elevation about the centre over the radius; 1 in Janbu's methods, whose force balance does not depend on where a horizontal force acts
Since , the passive force is the larger of the two for any target above one, and the two coincide at , an identity the program's tests check. Negative values, for a surface that already stands above the target, are reported as zero. Every valid surface of the search is back-analysed and the one that needs the largest force is reported, with both values. Janbu's corrected method solves its uncorrected factor and multiplies it by , so its sums and force are taken at .
Support → Back Analysis of Support Force takes the target factor (1.3 by default), the elevation of the force and the method: Bishop's simplified method, on circular surfaces only, or Janbu's simplified or corrected method. The sums are the method's own, with the supports already in the model applied as the method applies them, so the force found is the one to add to the existing reinforcement. Composite surfaces are left out and counted in the notes. The calculation neither uses nor changes the results of the stability analysis.
Source: ogr_slip2d/back_analysis.py · ogr_slip2d/analysis_runner.py3.10.7References#
- Berg, R.R., Christopher, B.R. & Samtani, N.C. (2009). Design and Construction of Mechanically Stabilized Earth Walls and Reinforced Soil Slopes, Vol. I. Report FHWA-NHI-10-024. Federal Highway Administration, Washington, DC.
- Cai, F. & Ugai, K. (2000). Numerical analysis of the stability of a slope reinforced with piles. Soils and Foundations, 40(1), 73–84.
- Duncan, J.M. & Wright, S.G. (2005). Soil Strength and Slope Stability. John Wiley & Sons, Hoboken, NJ.
- Esterhuizen, J.J.B., Filz, G.M. & Duncan, J.M. (2001). Constitutive behavior of geosynthetic interfaces. Journal of Geotechnical and Geoenvironmental Engineering, 127(10), 834–840.
- Ito, T. & Matsui, T. (1975). Methods to estimate lateral force acting on stabilizing piles. Soils and Foundations, 15(4), 43–59.
- Jewell, R.A. (1996). Soil Reinforcement with Geotextiles. CIRIA Special Publication 123. CIRIA and Thomas Telford, London.
- Mitsch, M.P. & Clemence, S.P. (1985). The uplift capacity of helix anchors in sand. In Clemence, S.P. (ed.), Uplift Behavior of Anchor Foundations in Soil. ASCE, New York, 26–47.
- Mooney, J.S., Adamczak, S. & Clemence, S.P. (1985). Uplift capacity of helical anchors in clay and silt. In Clemence, S.P. (ed.), Uplift Behavior of Anchor Foundations in Soil. ASCE, New York, 48–72.
- Perko, H.A. (2009). Helical Piles: A Practical Guide to Design and Installation. John Wiley & Sons, Hoboken, NJ.
- Prandtl, L. (1921). Über die Eindringungsfestigkeit (Härte) plastischer Baustoffe und die Festigkeit von Schneiden. Zeitschrift für Angewandte Mathematik und Mechanik, 1(1), 15–20.
- Reissner, H. (1924). Zum Erddruckproblem. Proceedings of the First International Congress for Applied Mechanics, Delft, 295–311.
- Terzaghi, K. (1943). Theoretical Soil Mechanics. John Wiley & Sons, New York.
- Terzaghi, K. & Peck, R.B. (1967). Soil Mechanics in Engineering Practice, 2nd ed. John Wiley & Sons, New York.
- Ukritchon, B. & Keawsawasvong, S. (2017). Error in Ito and Matsui's limit-equilibrium solution of lateral force on a row of stabilizing piles. Journal of Geotechnical and Geoenvironmental Engineering, 143(9), 02817004.