⚡ Elektrotechnik

Limit-state geometry of trailing-edge tensile cracks in soil slopes with a tensile strength threshold

Diese Studie präsentiert eine neuartige Grenzzustandsformulierung, die durch die Behandlung der Zugfestigkeit und des Scherwiderstands als unabhängige Versagensmechanismen gleichzeitig die Tiefe und Neigung von Zugrissen an der Hinterkante von zementierten Böschungen bestimmt und somit einen umfassenderen analytischen Rahmen bietet, der durch theoretische Benchmarks, Felddaten und diskrete Element-Simulationen validiert wurde.

Ursprüngliche Autoren: Jijie Du, Liansheng Tang, Jianing Huang, Zihua Cheng

Veröffentlicht 2026-11-02

📖 4 Min. Lesezeit☕ Kaffeepausen-Lektüre

CC BY 4.0

Ursprüngliche Autoren: Jijie Du, Liansheng Tang, Jianing Huang, Zihua Cheng

Originalarbeit lizenziert unter CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Dies ist eine KI-generierte Erklärung des untenstehenden Papers. Sie wurde nicht von den Autoren verfasst oder gebilligt. Für technische Genauigkeit konsultieren Sie das Originalpaper. Vollständigen Haftungsausschluss lesen

Hänge sind keine statischen Monumente; sie sind lebendige Systeme, die langsam der Schwerkraft nachgeben. Wenn ein Hang zu versagen beginnt, geschieht dies selten auf einmal. Stattdessen bereitet sich der Boden oft durch Risse an der obersten Kante, nahe dem Kamm, auf eine Rutschung vor, bevor der Hauptkörper der Erde nach unten gleitet. Diese anfänglichen Risse, bekannt als rückseitige Zugrisse (trailing-edge tensile cracks), sind mehr als nur oberflächliche Narben. Sie fungieren als Warnsignale und, entscheidend, als Pfade, die es Regenwasser ermöglichen, tief in den Boden einzusickern und den Kollaps zu beschleunigen. Für Ingenieure und Geologen ist das Verständnis darüber, wie tief diese Risse gehen und in welchem Winkel sie geneigt sind, von entscheidender Bedeutung, um vorherzusagen, wann ein Hang nachgeben könnte. Traditionell beruhten Modelle zur Vorhersage dieser Risse auf einer vereinfachenden Annahme: dass die Risse perfekt vertikal verlaufen, wie eine gerade senkrecht geschnittene Wand. Dies erleichtert die Mathematik, ignoriert jedoch die Realität, dass Boden oft in einem Winkel bricht und dass die Fähigkeit des Bodens, dem Auseinanderziehen zu widerstehen, ebenso wichtig ist wie seine Fähigkeit, dem Zusammengedrücktwerden zu widerstehen.

Ein Team von Forschern hat nun eine neue Methode zur Berechnung der Form dieser Risse entwickelt, die die Annahme der Vertikalität aufhebt. Anstatt den Riss zu zwingen, senkrecht nach oben zu verlaufen, behandelt ihre Methode sowohl die Tiefe als auch den Winkel des Risses als gleichzeitig zu lösende Unbekannte. Sie entwickelten ein Modell für Hänge aus zementiertem Boden – Erde, die durch natürliche oder künstliche Bindemittel eine gewisse Festigkeit besitzt – unter Bedingungen, bei denen Wasser frei abfließen kann. Das Modell balanciert die Kräfte der Schwerkraft, die den Boden nach unten ziehen, gegen die interne Festigkeit des Bodens aus, welche aus zwei Quellen stammt: der Reibung und Klebrigkeit, die den Boden zusammenhalten, wenn er zusammengedrückt wird, und einem separaten Schwellenwert, der bestimmt, wie viel Zugkraft der Boden aushalten kann, bevor er reißt. Durch das Lösen des Punktes, an dem diese Kräfte ein kritisches Limit erreichen, können die Forscher die wahrscheinlichste Geometrie des Risses bestimmen, ohne den Winkel im Voraus erraten zu müssen.

Um ihren neuen Ansatz zu testen, verglich das Team seine Berechnungen mit etablierten Theorien, realen Messungen eines Loess-Hanges in Nordwestchina und fortschrittlichen Computersimulationen, die das Verhalten einzelner Bodenpartikel nachahmen. Die Ergebnisse zeigten, dass ihre Methode gut mit bestehenden Theorien darüber übereinstimmt, wie tief Risse unter verschiedenen Bedingungen sein sollten. Als sie ihr Modell auf den realen Loess-Hang anwandten, an dem ein Riss mit einer Tiefe von fast 13 Metern gemessen worden war, ergaben ihre Berechnungen einen Bereich möglicher Tiefen, der die tatsächliche Messung einschloss, sofern eine vernünftige Schätzung der Zugfestigkeit des Bodens verwendet wurde. Die Computersimulationen, die das Brechen der Bindungen zwischen virtuellen Bodenpartikeln verfolgten, bestätigten, dass das Modell korrekt identifizierte, wann ein Riss entsteht und wann nicht. Obwohl die exakte Tiefe, die das Modell vorhersagte, oft tiefer war als das, was die Computersimulationen zeigten, stimmte der vom neuen Verfahren vorhergesagte Winkel der Risse mit bemerkenswerter Präzision mit den simulierten Rissen überein und wich im Durchschnitt nur um wenige Grad ab.

Die Studie zeigte auch auf, wie verschiedene Bodeneigenschaften die Form des Risses verändern. Die Forscher fanden heraus, dass die Risse flacher werden und steiler geneigt sind, wenn die Fähigkeit des Bodens, dem Auseinanderziehen zu widerstehen, zunimmt; manchmal sind sie so widerstandsfähig, dass sich gar kein Riss bildet. Umgekehrt neigen die Risse dazu, tiefer zu gehen, wenn der Boden klebriger oder kohäsiver ist. Auch das Gewicht des Bodens selbst spielt eine Rolle: Schwererer Boden führt zu flacheren Rissen, hat jedoch kaum Einfluss auf den Winkel. Interessanterweise hat der interne Reibungswinkel des Bodens keinen einfachen, einseitigen Effekt auf die Form des Risses; sein Einfluss hängt stark davon ab, wie stark der Boden gegen das Auseinanderziehen ist. Durch eine detaillierte Analyse der Wechselwirkungen dieser Faktoren stellten das Team fest, dass die Klebrigkeit des Bodens der dominanteste Faktor bei der Entscheidung ist, ob ein Riss entsteht und wie tief er wird, dicht gefolgt von der Zugfestigkeit des Bodens.

Diese Arbeit bietet ein vollständigeres Bild davon, wie Hänge versagen, indem sie die Geometrie des Risses als Ergebnis physikalischer Kräfte statt als feste Regel behandelt. Indem sie anerkennt, dass Risse geneigt sein können und dass der Widerstand des Bodens gegen Zug eine andere Eigenschaft ist als der Widerstand gegen das Gleiten, bietet das neue Modell ein nuancierteres Werkzeug zur Bewertung der Hangstabilität. Es erhebt nicht den Anspruch, jeden zukünftigen Erdrutsch vorherzusagen, bietet aber eine klarere, mechanisch begründete Referenz für das Verständnis der Anfangsphasen eines Hangversagens. Die Ergebnisse legen nahe, dass das Ignorieren der Zugfestigkeit des Bodens oder die Annahme, dass Risse immer vertikal verlaufen, zu ungenauen Risikobewertungen führen kann, insbesondere bei zementierten Böden, in denen die Fähigkeit, Zugspannungen zu widerstehen, ein wesentlicher Teil des Materialcharakters ist.

Technische Zusammenfassung: Grenzzustandsgeometrie von Zugrissen an der Hinterkante in Bodenhängen mit einer Zugfestigkeitsschwelle

Problemstellung
Zugrisse an der Hinterkante sind kritische Vorläufer von Zug-Scher-Böschungsbrüchen, die die Lastübertragungswege verändern und die Wasserinfiltration erleichtern. Während bestehende analytische Modelle die Risstiefe effektiv abschätzen können, legen sie die Rissneigung typischerweise als vertikale Ebene fest und vernachlässigen dabei oft die Zugfestigkeit ( $\sigma_t$ ) oder behandeln diese indirekt durch modifizierte Scherparameter. Diese Studie adressiert die Einschränkung, die Rissgeometrie als fest vorgegebenen Input statt als gekoppeltes Ergebnis zu behandeln. Die Autoren argumentieren, dass die Neigung eines Zugrisses an der Hinterkante die Geometrie der potenziell abgelösten Bodenmasse und die anschließende Stabilitätsbewertung signifikant beeinflusst, was eine Formulierung erforderlich macht, bei der sowohl Tiefe als auch Neigung gleichzeitig unter mechanischen Randbedingungen bestimmt werden.

Methodik
Die Studie entwickelt eine zweidimensionale Grenzzustandsformulierung für einen ursprünglich intakten, homogenen zementierten Bodenhang unter drainierten Bedingungen. Das Modell repräsentiert ein lokales Zug-Scher-Keil mit einer durch einen Zugriss an der Hinterkante (BF) und eine untere Scherfläche (EF) begrenzten Volumen.

Grundprinzipien: Der Scherwiderstand folgt dem Mohr-Coulomb-Kriterium ( $c, \phi$ ), während das Zugversagen unabhängig durch eine makroskopische Zugfestigkeitsschwelle ( $\sigma_t$ ) gesteuert wird.
Gekoppelte Unbekannte: Im Gegensatz zu früheren Ansätzen werden die Risstiefe ( $L$ ) und die Neigung ( $\delta$ ) als gekoppelte Unbekannte behandelt. Die Lösung muss das Kräftegleichgewicht, das Momentengleichgewicht, die Zugversagensbedingungen und die geometrische Zulässigkeit ( $0 < \theta < \delta < \pi/2$ ) erfüllen.
Lösungsverfahren: Aufgrund des Fehlens einer geschlossenen Lösung für die gekoppelten Winkelvariablen verwenden die Autoren einen hierarchischen, einschränkenden Suchalgorithmus. Dies beinhaltet die Diskretisierung des Winkels der Scherfläche ( $\theta$ ), das Lösen für zulässige Rissneigungen ( $\delta$ ), welche das Momentengleichgewicht erfüllen, und die Auswahl der minimalen positiven Grenzzustands-Tiefe ( $L$ ) als die maßgebliche Lösung.
Validierung: Die Formulierung wird durch folgende Methoden evaluiert:
1. Theoretische Benchmarks: Vergleich mit etablierten Lösungen für vertikale Risstiefen (z. B. Michalowski, Utili, Zhu et al.).
2. Feldvergleich: Eine gebundene Konsistenzprüfung anhand eines Loess-Hangfalls im Nordwesten Chinas.
3. Numerische Simulation: Vergleich mit MatDEM-Simulationen (Discrete Element Method), bei denen die Rissgeometrie aus dem Bindungsbruch und der Kraftkettenentwicklung extrahiert wird.
4. Sensitivitätsanalyse: Lokale Sensitivitätsanalyse und globale Sobol-Sensitivitätsanalyse zur Quantifizierung des Paramendereinflusses.

Wesentliche Ergebnisse

Benchmark-Konsistenz: Das Modell reproduziert etablierte Trends der Risstiefe in Bezug auf Kohäsion, Einheitsgewicht und Reibungswinkel. Entscheidend ist, dass es eine explizite Abhängigkeit von der Zugfestigkeit einführt und zeigt, dass eine höhere Zugfestigkeit die zulässige Risstiefe reduziert.
Feldvergleich: Für einen Loess-Hang mit einer gemessenen Risstiefe von 12,78 m sagt das Modell eine Spanne von Tiefen (9,90 m bis 21,80 m) voraus, abhängig vom Verhältnis der Zugfestigkeit ( $T = \sigma_t/c$ ). Die Feldmessung liegt innerhalb dieses berechneten Bereichs und schneidet die Modellkurve bei $T \approx 0,7\text{--}0,8$ .
DEM-Validierung: Vergleiche mit MatDEM-Simulationen zeigen, dass das analytische Modell die Unterscheidung zwischen Riss- und Nicht-Riss-Zuständen bewahrt. Während die analytische Tiefe ( $L$ ) tendenziell größer ist als die aus DEM extrahierten Tiefen (mittlerer absoluter Fehler $\approx$ 1,66 m), zeigen die Neigungswinkel eine engere Übereinstimmung (mittlerer absoluter Fehler $\approx$ 4,22°). Die Autoren führen die Diskrepanzen in der Tiefe auf den Unterschied zwischen einem theoretischen Grenzzustandsreferenzpunkt und der spezifischen Phase der Rissentwicklung in der DEM zurück.
Parametrische Effekte:
- Zugfestigkeit ( $\sigma_t$ ): Eine Erhöhung von $\sigma_t$ reduziert sowohl $L$ als auch $\delta$ . Bei hohen Schwellenwerten kann sie die Bildung einer endlichen Lösung mit positiver Tiefe vollständig unterdrücken.
- Kohäsion ( $c$ ): Erhöht im Allgemeinen sowohl $L$ als auch $\delta$ und erweitert das Gebiet, in dem eine endliche Lösung existiert.
- Einheitsgewicht ( $\gamma$ ): Reduziert primär $L$ mit vernachlässigem Effekt auf $\delta$ .
- Reibungswinkel ( $\phi$ ): Sein Effekt auf $L$ ist nicht-monoton und von $\sigma_t$ abhängig; bei hohen Zugfestigkeiten kann eine Erhöhung von $\phi$ $L$ nicht erhöhen, da die Zugschwelle restriktiv wirkt.
Sensitivität: Eine lokale Sensitivität nahe der Baseline identifiziert die Kohäsion als den einflussreichsten Parameter für die Tiefe. Eine globale Sobol-Analyse ordnet die Parameter für die Tiefe (einschließlich des Null-Wertes) wie folgt: Kohäsion > Zugfestigkeit > Reibungswinkel > Einheitsgewicht. Die Zugfestigkeit ist besonders signifikant für ihre Rolle beim Übergang zwischen Zuständen mit endlicher Tiefe und Zuständen der Rissunterdrückung.

Bedeutung und Ansprüche
Die Autoren behaupten, dass diese Studie eine mechanisch gebundene Grundlage für die Interpretation der Geometrie von Zugrissen an der Hinterkante liefert, indem sie das Problem von einer vorgegebenen geometrischen Bedingung zu einem Satz gekoppelter Unbekannter umformuliert. Der primäre Beitrag ist die explizite Integration der Zugfestigkeit als unabhängiges Versagenkriterium neben der Scherfestigkeit, was die gleichzeitige Bestimmung von Risstiefe und Neigung ermöglicht.

Die Arbeit wahrt einen bescheidenen Umfang und stellt klar, dass die Formulierung eine idealisierte lokale Grenzzustandsgeometrie liefert und keine zeitabhängige Vorhersage der Rissausbreitung oder eines vollständigen Böschungsversagensprozesses darstellt. Sie ist als Referenz für die Interpretation von Versagen, Rückanalysen oder nachfolgende Stabilitätsbewertungen gedacht, insbesondere für zementierte Böden, in denen der Zugwiderstand nicht vernachlässigbar ist. Die Studie beansprucht nicht, fortgeschrittene numerische Simulationen zu ersetzen, sondern bietet vielmehr ein effizientes, physikalisch interpretierbares analytisches Werkzeug, das den gekoppelten Einfluss von Zug- und Scherversagensmechanismen erfasst.

Mehr davon