OpenGeoRock
3.10

Support forces

Describes OGR Slip2D 0.1.235

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, F(d)F(d), with dd 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 SS.

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 α\alpha:

TS=F⋅t=s (Fxcos⁡α+Fysin⁡α),TN=F⋅n=Fxsin⁡α−Fycos⁡αT_S = \mathbf{F}\cdot\mathbf{t} = s\,(F_x\cos\alpha + F_y\sin\alpha), \qquad T_N = \mathbf{F}\cdot\mathbf{n} = F_x\sin\alpha - F_y\cos\alpha

where

Fx, FyF_x,\ F_y
components of the support force, xx to the right and yy upwards (kN/m)
ss
sense of sliding, ±1, chosen so that t=s (cos⁡α, sin⁡α)\mathbf{t} = s\,(\cos\alpha,\ \sin\alpha) is the unit vector along the base that opposes sliding
n=(sin⁡α, −cos⁡α)\mathbf{n} = (\sin\alpha,\ -\cos\alpha)
unit normal pointing from the slice into the ground below its base
TST_S
tangential component, positive when it resists sliding (kN/m)
TNT_N
normal component, positive when it presses the slice onto its base (kN/m)
A support crossing a slip surface: the force acts at the crossing P, on the base of the slice that contains it, and is split into a tangential component T_S that resists sliding and a normal component T_N that presses the slice onto its base slip surface slice that contains P head (slope face) tail (anchor) d L − d P F TS TN sliding
Figure 3.10.1 A support crossing a slip surface. The force FF available at the crossing PP, here parallel to the support (head to tail), acts on the base of the slice that contains PP and is resolved into TST_S, along the base against sliding, and TNT_N, normal to the base and into the ground. dd is the distance from the head to the crossing and L−dL - d the length beyond it.

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 RR and DD the resisting and driving terms of a method:

F=R+TNtan⁡ϕ′D−TS(active)F = \dfrac{R + T_N\tan\phi'}{D - T_S} \qquad \text{(active)}
F=R+TNtan⁡ϕ′+TSD(passive)F = \dfrac{R + T_N\tan\phi' + T_S}{D} \qquad \text{(passive)}

where

R, DR,\ D
resisting and driving terms of the method without the support (kN/m, or moments over the radius on a circle)
TS, TNT_S,\ T_N
tangential and normal components of Eq. (3.10.1), summed over the supports and slices crossed
ϕ′\phi'
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

Sm=c′l+(N−u l)tan⁡ϕ′F+TS act+TS pasFS_m = \dfrac{c'l + (N - u\,l)\tan\phi'}{F} + T_S^{\,act} + \dfrac{T_S^{\,pas}}{F}

with NN the base normal force, which already carries TNT_N through the equilibrium of the slice. The two readings differ only in that line: summed over a mass in force equilibrium, it turns F=R/(D−T)F = R/(D - T) into F=(R+T)/DF = (R + T)/D. 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 TST_S if active and TS/FT_S/F 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 mαm_\alpha (or nαn_\alpha) like the weight. Janbu's methods, which balance shear along the base, take TSsec⁡αT_S\sec\alpha on the driving side.
  • The ordinary method resolves forces normal to the base, so TNT_N 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 TST_S; a support tangent to the slip surface takes the arm RR, like a base shear.

On a planar slip surface the mass is a single free body and the two readings have closed forms,

Fact=c′L+(Wcos⁡β+TN)tan⁡ϕ′Wsin⁡β−TS,Fpas=c′L+(Wcos⁡β+TN)tan⁡ϕ′+TSWsin⁡βF_{act} = \dfrac{c'L + (W\cos\beta + T_N)\tan\phi'}{W\sin\beta - T_S}, \qquad F_{pas} = \dfrac{c'L + (W\cos\beta + T_N)\tan\phi' + T_S}{W\sin\beta}

with WW the weight of the wedge and β, L\beta,\ L 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.

In OGR Slip2D

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 TST_S approaches DD 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.py

3.10.3Orientation of the force#

The direction of FF 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 xx 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.

Table 3.10.1 The nine support types of OGR Slip2D, their defaults and the force at the crossing (per metre of slope). SS is the out-of-plane spacing.
TypeOrientationApplicationForce at the crossing
End anchoredparallel to supportactiveConstant: anchor capacity / SS
Grouted tiebackparallel to supportactiveSmallest of pullout, tensile and stripping; bond strength per metre over a bonded length at the tail
Grouted tieback with frictionparallel to supportactiveAs the grouted tieback, with the bond from an interface law integrated around the grout
Soil nailparallel to supportpassiveAs the grouted tieback, bonded over its whole length
Pile / micropiletangent to slip surfacepassiveShear: pile shear strength / SS; or the lateral thrust of Ito & Matsui (1975)
Geosyntheticparallel to supportpassiveSmallest of pullout, tensile and stripping, both faces of the sheet, times the coverage
User definedparallel to supportactiveTable of force against distance from the head, / SS
Retaining wall (EFP)horizontalactivePressure profile integrated from the crest down to the crossing; no spacing
Helical anchortangent to slip surfaceactiveSmallest 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 qbq_b per metre over a bonded length at the tail end:

F(d)=1S min⁡{ qb Lo,  T,  P+qb Li }F(d) = \dfrac{1}{S}\,\min\left\{\, q_b\,L_o,\ \ T,\ \ P + q_b\,L_i \,\right\}

where

LoL_o
bonded length beyond the crossing, towards the tail: the whole bond when the crossing is in the free length (m)
LiL_i
bonded length between the head and the crossing: zero when the crossing is in the free length (m)
T, PT,\ P
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, Lo=L−dL_o = L - d and Li=dL_i = d. In the grouted tieback with friction the bond is the integral of an interface strength τ\tau over the surface of a grout cylinder of diameter DD,

F(d)=1S min⁡{ πD ⁣∫Lo ⁣τ ds,  T,  P+πD ⁣∫Li ⁣τ ds }F(d) = \dfrac{1}{S}\,\min\left\{\, \pi D\!\int_{L_o}\!\tau\,\mathrm{d}s,\ \ T,\ \ P + \pi D\!\int_{L_i}\!\tau\,\mathrm{d}s \,\right\}

and in a geosynthetic, which mobilises both faces of the sheet (Berg, Christopher & Samtani, 2009), with AA the percentage of the width of slope its strips cover (100 % for a continuous sheet),

F(d)=min⁡{ 2A100 ⁣∫Lo ⁣τ ds,  A100 T,  A100[C+2 ⁣∫Li ⁣τ ds] }F(d) = \min\left\{\, \dfrac{2A}{100}\!\int_{L_o}\!\tau\,\mathrm{d}s,\ \ \dfrac{A}{100}\,T,\ \ \dfrac{A}{100}\Big[C + 2\!\int_{L_i}\!\tau\,\mathrm{d}s\Big] \,\right\}

with TT the tensile capacity per metre of sheet, CC the strength of the connection at the face, LiL_i from the face to the crossing and LoL_o from the crossing to the embedded end. A sheet whose embedded end is anchored cannot be pulled out, and its first term is dropped.

Available force along a soil nail and along a grouted tieback with the default parameters of OGR Slip2D: pullout, tensile and stripping capacities and the applied force, which is their minimum Soil nail, L = 10 m 0 50 100 150 200 0 2 4 6 8 10 pullout tensile stripping applied force distance from head (m) Grouted tieback, L = 12 m 0 100 200 300 400 0 2 4 6 8 10 12 pullout tensile stripping applied force distance from head (m) force (kN/m)
Figure 3.10.2 Force available at the crossing along a 10 m soil nail and a 12 m grouted tieback with the default parameters of OGR Slip2D (nail: T=100T = 100 kN, P=50P = 50 kN, qb=30q_b = 30 kN/m, S=1.5S = 1.5 m; tieback: T=600T = 600 kN, P=300P = 300 kN, qb=150q_b = 150 kN/m over the last 30 % of its length, S=2S = 2 m). The applied force (solid) is the smallest capacity at each point. In the nail, stripping governs near the head, the tendon in the middle and pullout near the tail; in the tieback, the plate governs over the free length and the start of the bond, and pullout beyond. Computed with OGR Slip2D 0.1.235.

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:

τ=a+σv′tan⁡δ(linear)\tau = a + \sigma'_v\tan\delta \qquad \text{(linear)}
τ=a∞ σv′tan⁡ϕ0a∞+σv′tan⁡ϕ0(hyperbolic)\tau = \dfrac{a_\infty\,\sigma'_v\tan\phi_0}{a_\infty + \sigma'_v\tan\phi_0} \qquad \text{(hyperbolic)}

where

a, δa,\ \delta
adhesion (kPa) and friction angle (°) of the interface, the Mohr–Coulomb bond of Jewell (1996)
a∞, ϕ0a_\infty,\ \phi_0
limiting strength at large normal stress (kPa) and tangent friction angle at zero normal stress (°), after Esterhuizen, Filz & Duncan (2001)
σv′\sigma'_v
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 τ\tau as a fraction CiC_i of the surrounding soil's own strength, evaluated with the material's strength model, or as τ=F∗σv′\tau = F^{*}\sigma'_v with a pullout friction factor F∗F^{*}, constant or varying linearly with depth (Berg, Christopher & Samtani, 2009); the scale-effect factor of that guidance is entered multiplied into F∗F^{*}.

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):

Qp=Ap(1.3 c Nc+σv′Nq)Q_p = A_p\left(1.3\,c\,N_c + \sigma'_v N_q\right)
Nq=eπtan⁡ϕtan⁡2 ⁣(45∘+ϕ2),Nc=(Nq−1)cot⁡ϕN_q = e^{\pi\tan\phi}\tan^2\!\left(45^\circ + \tfrac{\phi}{2}\right), \qquad N_c = (N_q - 1)\cot\phi

with ApA_p the plate area less the area of the shaft, NqN_q and NcN_c the factors of Prandtl (1921) and Reissner (1924) (Nc→2+πN_c \to 2 + \pi as ϕ→0\phi \to 0), and c, ϕc,\ \phi 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.

Note

At ϕ=0\phi = 0 the plate bearing becomes 1.3×5.14=6.681.3 \times 5.14 = 6.68 times the undrained strength, whereas Mooney, Adamczak & Clemence (1985) use Nc=9N_c = 9 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):

p=c A[E−2Nϕtan⁡ϕ−1Nϕtan⁡ϕ+G]−c[D1G−2D2Nϕ]+σv′Nϕ(A E−D2)p = c\,A\left[\dfrac{E - 2\sqrt{N_\phi}\tan\phi - 1}{N_\phi\tan\phi} + G\right] - c\left[D_1 G - \dfrac{2D_2}{\sqrt{N_\phi}}\right] + \dfrac{\sigma'_v}{N_\phi}\left(A\,E - D_2\right)

with

Nϕ=tan⁡2 ⁣(π4+ϕ2),A=D1(D1D2)Nϕtan⁡ϕ+Nϕ−1E=exp⁡ ⁣[D1−D2D2 Nϕtan⁡ϕ tan⁡ ⁣(π8+ϕ4)],G=2tan⁡ϕ+2Nϕ+Nϕ−1/2Nϕtan⁡ϕ+Nϕ−1\begin{gathered} N_\phi = \tan^2\!\left(\tfrac{\pi}{4} + \tfrac{\phi}{2}\right), \qquad A = D_1\left(\dfrac{D_1}{D_2}\right)^{\sqrt{N_\phi}\tan\phi + N_\phi - 1} \\[4pt] E = \exp\!\left[\dfrac{D_1 - D_2}{D_2}\,N_\phi\tan\phi\,\tan\!\left(\tfrac{\pi}{8} + \tfrac{\phi}{4}\right)\right], \qquad G = \dfrac{2\tan\phi + 2\sqrt{N_\phi} + N_\phi^{-1/2}}{\sqrt{N_\phi}\tan\phi + N_\phi - 1} \end{gathered}

where

D1D_1
centre-to-centre spacing of the piles in the row (m)
D2=D1−dpD_2 = D_1 - d_p
opening between neighbouring piles, with dpd_p the pile diameter (m)
c, ϕc,\ \phi
equivalent Mohr–Coulomb parameters of the soil at that depth, from the linearisation the methods use
σv′\sigma'_v
effective vertical stress, which takes the place of the γz\gamma z of the original paper (kPa)

For ϕ\phi below 10−810^{-8} rad the program evaluates the cohesive solution of the same paper, p=c[D1(3ln⁡D1D2+D1−D2D2tan⁡π8)−2(D1−D2)]+σv′(D1−D2)p = c\left[D_1\left(3\ln\frac{D_1}{D_2} + \frac{D_1 - D_2}{D_2}\tan\frac{\pi}{8}\right) - 2(D_1 - D_2)\right] + \sigma'_v (D_1 - D_2), 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 pp from the top of the pile down to the crossing, divided by the spacing, as Cai & Ugai (2000) write it:

F=1D1∫0zsp(z) dzF = \dfrac{1}{D_1}\int_0^{z_s} p(z)\,\mathrm{d}z

with zsz_s 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 pp 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 (D2≤0D_2 \le 0, 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).

Note

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 p(s)p(s) given over its length, from its upper end down to the crossing,

F(d)=∫0dp(s) ds,triangular profile: p=γe s ⇒ F(d)=12 γe d2F(d) = \int_0^{d} p(s)\,\mathrm{d}s, \qquad \text{triangular profile: } p = \gamma_{e}\,s \ \Rightarrow\ F(d) = \tfrac12\,\gamma_{e}\,d^2

with γe\gamma_e the equivalent fluid pressure (kN/m³) and dd 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 γeL\gamma_e L over (1−f)/2(1 - f)/2 of the length at each end and constant over the fraction ff 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.

In OGR Slip2D

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.py

3.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

Pw=12 γw hw2,horizontal, at hw3 above the base of the waterP_w = \tfrac12\,\gamma_w\,h_w^2, \quad\text{horizontal, at } \tfrac{h_w}{3} \text{ above the base of the water}

where

hwh_w
height of water in the crack, from the crack boundary at the wall up to the water level (m)
γw\gamma_w
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.

A slip surface truncated by a tension crack: the sliding mass is closed by a vertical wall from the crack boundary up to the ground, and water in the crack pushes on that wall with the resultant one half gamma_w h_w squared, at h_w/3 above the base of the water tension crack boundary crack zone surface inside the zone: removed wall water level Pw = ½γwhw2 hw hw/3 first slice: carries Pw sliding
Figure 3.10.3 A slip surface truncated by a tension crack. The part of the surface inside the crack zone is removed and the mass is closed by a vertical wall from the crack boundary up to the ground. The hydrostatic pressure of the water in the crack (here filled to 70 % of its height) has the resultant PwP_w of Eq. (3.10.17), applied to the first slice at hw/3h_w/3 above the base of the water.
In OGR Slip2D

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.py

3.10.6Back analysis of the required support force#

In the back analysis the factor of safety is fixed at a target FtF_t and the unknown is the horizontal force TT 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:

Tact=1a(D−RFt),Tpas=1a(Ft D−R)T_{act} = \dfrac{1}{a}\left(D - \dfrac{R}{F_t}\right), \qquad T_{pas} = \dfrac{1}{a}\left(F_t\,D - R\right)

where

R, DR,\ D
resisting sum evaluated at FtF_t, and driving sum, of the method
aa
arm factor: (yc−yT)/R(y_c - y_T)/R in Bishop's simplified method on a circle, the moment arm of a horizontal force at elevation yTy_T about the centre over the radius; 1 in Janbu's methods, whose force balance does not depend on where a horizontal force acts

Since Tpas=Ft TactT_{pas} = F_t\,T_{act}, the passive force is the larger of the two for any target above one, and the two coincide at Ft=1F_t = 1, 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 f0f_0, so its sums and force are taken at Ft/f0F_t/f_0.

In OGR Slip2D

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.py

3.10.7References#