from pybimstab.slope import AnthropicSlope
from pybimstab.slipsurface import CircularSurface
from pybimstab.slices import MaterialParameters, Slices
slope = AnthropicSlope(slopeHeight=7.5,
slopeDip=[1, 1.5],
crownDist=5,
toeDist=5)
surface = CircularSurface(slopeCoords=slope.coords,
dist1=2,
dist2=10,
radius=9)
material = MaterialParameters(cohesion=15,
frictAngle=23,
unitWeight=17,
blocksUnitWeight=21,
wtUnitWeight=9.8)
slices = Slices(material=material,
slipSurfCoords=surface.coords,
slopeCoords=slope.coords,
numSlices=10,
watertabCoords=None,
bim=None)
fig = slices.plot()
from pybimstab.slope import AnthropicSlope
slope = AnthropicSlope(slopeHeight=12,
slopeDip=[1, 1.5],
crownDist=10,
toeDist=10)
fig = slope.plot()
# Example Case 5 (Fredlund & Krahn, 1977)
from numpy import array
from pybimstab.slope import AnthropicSlope
from pybimstab.slipsurface import CircularSurface
from pybimstab.watertable import WaterTable
from pybimstab.slices import MaterialParameters, Slices
from pybimstab.slopestabl import SlopeStabl
slope = AnthropicSlope(slopeHeight=40,
slopeDip=[2, 1],
crownDist=60,
toeDist=30,
depth=20)
surface = CircularSurface(slopeCoords=slope.coords,
dist1=45.838,
dist2=158.726,
radius=80)
material = MaterialParameters(cohesion=600,
frictAngle=20,
unitWeight=120,
wtUnitWeight=62.4)
watertable = WaterTable(slopeCoords=slope.coords,
watertabDepths=array([[0, 140], [20, 0]]))
slices = Slices(material=material,
slipSurfCoords=surface.coords,
slopeCoords=slope.coords,
numSlices=50,
watertabCoords=watertable.coords,
bim=None)
stabAnalysis = SlopeStabl(slices, seedFS=1, Kh=0)
fig = stabAnalysis.plot()