Program:
Język:

Morgenstern-Price

Die Spencer-Methode ist ein allgemeines Lamellenverfahren, das auf dem Grenzgleichgewicht beruht. Es benötigt ein befriedingendes Kräfte- und Momentengleichgewicht, das auf einzelne Blöcke wirkt. Die Blöcke (Lamellen) entstehen durch Aufteilen des Bodens über der Gleitfläche furch Trennebenen. Die Kräfte, die auf einen einzelnen Block wirken, sind in der folgenden Abbildung dargestellt:

Statisches Schema: Morgenstern-Price-Methode

Es wird angenommen, dass jede Lamelle die gleichen Kräfte besitzt, wie in der Spencer-Methode. Die folgenden Annahmen werden für die Morgenstern-Price-Methode eingeführt, um das Grenzgleichgewicht der Kräfte und Momente der einzelnen Blöcke zu berechnen:

  • Trennflächen zwischen den Lamellen sind immer vertikal
  • die Wirkungslinie des Blockgewichts Wi verläuft durch die Mitte des i-ten Abschnitts der Gleitfläche, den Punkt M
  • die Normalkraft Ni wirkt in der Mitte des i-ten Abschnitts der Gleitfläche, im Punkt M
  • die Neigung der Kräfte Ei wirkt zwischen den Blöcken und ist unterschiedlich für jeden Block (δi), nur am Ende der Gleitfläche ist δ = 0

Der einzige Unterschied zwischen Spencer und Morgenstern-Price ist in der obigen Liste der Annahmen aufgeführt. Die Wahl des Winkels δi der Kraft Ei, die zwischen den Blöcken wirkt, wird mithilfe einer Halb-Sinusfunktion getroffen - eine der Funktionen der nachfolgenden Grafik wird automatisch ausgewählt. Die Wahl der Form der Kurve hat kaum Einfluss auf das Endergebnis, allerdings kann eine sinnvolle Wahl die Konvergenz der Methode verbessern. Der Funktionswert der Halbsinusfunktion f(xi) am Schnittstellenpunkt xi multipliziert mit dem Parameter λ liefert den Wert des Neigungswinkels δi.

Halb-Sinusfunktion

Die Lösung herauskommt aus der Gleichungen (1) - (5) wie in der Spencer-Methode:

(1)

(2)

(3)

(4)

(5)

wo:

φi

-

Winkel der inneren Reibung auf dem Gleitflächenabschnitt

ci

-

Kohäsion des Bodens auf dem Gleitflächenabschnitt

αi

-

Neigung des Gleitflächenabschnitts

  • (1) Beziehung zwischen effektivem und Gesamtwert der Normalkraft, die auf die Gleitfläche wirkt
  • (2) Bedingungen nach Mohr-Coulomb, die die Beziehung zwischen Normal- und Scherkräften auf ein gegebenen Abschnitt der Gleitfläche darstellen (Ni und Ti)
  • (3) Kräftegleichgewicht senkrecht zum iten Abschnitt
  • (4) Gleichgewicht entlang des i-ten Abschnitts der Gleitfläche
  • (5) Momentengleichgewicht um den Punkt M

Das Umstellen der Gleichungen (3) und (4) ergibt rekursiver Formel (6):

(6)

Diese Formel erlaubt die Berechnung aller Kräfte Ei, die zwischen den Blöcken wirken, für gegebene Werte von δi und SF. Die Lösung geht davon aus, daß der Wert von E am Ursprung der Gleitfläche bekant und gleich E1=0 ist.

Aus den Momentengleichgewicht (5) ergibt sich eine weitere rekursive Formel (7):

(7)

Diese Formel erlaubt uns die Berechnung für gegebene Werte für δ aller Hebelarme z der Kräfte zwischen den Blöcken, bei Kenntnis des Werts am linken Rand der Gleitfläche, wo z1 = 0

Der Sicherheitsfaktor SF wird durch Anwendung des folgenden Iterationsprozesses bestimmt:

  1. Der Anfangswert des Winkels δi wird entsprechend der Halbsinusfunktion zu (δi = λ*f(xi)) bestimmt.
  2. Der Sicherheitsfaktor SF für einen Wert von δi folgt aus Gleichung (6), während man annimmt, dass En+1 = 0 am Ende der Gleitfläche ist.
  3. Den Wert δi liefert die Gleichung (7) indem man den Wert E, der im vorherigen Schritt ermittelt wurde mit der Anforderung, dass das Moment im letzten Block gleich Null ist. Die Funktionswerte f(xi) sind alle gleich während der Iteration, nur der Parameter λ wird iteriert. Gleichung (7) liefert nicht den Wert von zn+1 da dieser gleich Null ist. Für diesen Wert muss das Momentengleichgewicht (5) erfüllt sein.
  4. Schritt 2 und 3 werden solange wiederholt, bis sich δi (respektive λ) nicht mehr verändert.

Damit der Iterationsprozess stabil ist, ist es notwendig, dass instabile Lösungen vermieden werden. Solche Instabilitäten treten auf, wenn in (6) und (7) durch Null geteilt wird. In Gleichung (7) wird für δ = π/2 oder δ = -π/2 durch Null dividiert. Demnach muss der Wert des Winkels δ im Intervall (-π/2 ; π/2) liegen.

Division durch Null in (6) tritt auf, wenn:

Ein anderer Weg, numerische Instabilität zu vermeiden, ist die Überprüfung des Parameters mα  - folgende Bedingung muss erfüllt sein:

Demnach muss man, bevor man die Iteration durchführt, den höchsten kritischen Wert von SFmin finden, der obige Bedingung erfüllt. Werte kleiner dem kritischen Wert SFmin sind im Bereich der instabilen Lösungen, demnach beginnt die Iteration bei einem Wert von SF der "ein wenig" größer als SFmin ist und somit sind alle anderen Werte von SF ebenfalls größer SFmin.

Generell konvergieren rigorose Methoden schlechter als einfache Methoden (Bishop, Fellenius). Beispiele mit Konvergenzproblemen sind unter anderem zu steile Abschnitte der Gleitfläche, komplexe Geometrie, ein signifikanter Sprung in der Auflast, usw. Falls keine Lösung gefunden wird, empfehlen wir eine geringfügige Änderung der Eingangsdaten, z. B. eine weniger steile Gleitfläche, Eingabe von mehr Punkten für die Gleitfläche, etc, oder alternativ die Berechnung mit einfachen Methoden.

Literatur:

Morgenstern, N.R., and Price, V.E. 1965. The analysis of the stability of general slip surfaces. Géotechnique, 15(1): 79-93.

Morgenstern, N.R., and Price, V.E. 1967. A numerical method for solving the equations of stability of general slip surfaces. Computer Journal, 9: 388-393.

Zhu, D.Y., Lee, C.F., Qian, Q.H., and Chen, G.R. 2005. A concise algorithm for computing the factor of safety using the Morgenstern-Price method. Canadian Geotechnical Journal, 42(1): 272-278.

Wypróbuj GEO5. Bezpłatnie.