Hybrid Dynamics of Rocking Blocks Beyond Overturning: Saltation Analysis, Bifurcations, and Stability Characterization

Technisches Resümee: Hybride Dynamik von kippenden Blöcken jenseits des Umkippens

Problemstellung
Kippende Blöcke sind fundamentale Modelle zum Verständnis des seismischen Antwortverhaltens starrer Strukturen wie Mauerwerkswände, Statuen und Hochbehälter. Während das klassische Housner-Modell seit langem als Eckpfeiler für die Analyse solcher Systeme dient, haben jüngere Studien Diskrepanzen zwischen dessen Vorhersagen und experimentell beobachteten Aufprallverhalten aufgezeigt. Insbesondere der klassische Restitutionskoeffizient (COR), der aus der Drehimpulserhaltung abgeleitet wird, versagt oft dabei, die Mechanismen der Energiedissipation präzise zu erfassen, insbesondere bei Blöcken mit geringen Schlankheitsgraden. Obwohl alternative Formulierungen, wie das von Mao et al. vorgeschlagene Modell der kinetischen Energieumverteilung, eingeführt wurden, um experimentelle Daten besser abzubilden, bleiben die Auswirkungen der Verwendung solcher alternativer Restitutionsgesetze auf die globale nichtlineare Dynamik von Kippsystemen weitgehend unerforscht. Diese Arbeit untersucht, wie die Wahl der Restitutionsmodellierung die Bifurkationsstruktur, die Stabilität und die Organisation der Einzugsbereiche (Basins of Attraction) von Kippsystemen beeinflusst, wobei sie über den traditionellen Fokus auf Umkippschwellen hinausgeht, um komplexe oszillatorische Regime zu charakterisieren.

Methodik
Die Studie formuliert den kippenden Block als ein hybrides, nicht-glattes dynamisches System, das einer harmonischen Basisanregung unterliegt. Das System ist definiert durch einen rechteckigen starren Körper mit Masse $M$, Breite $2b$ und Höhe $2h$, der ohne Gleiten kippt. Die Dynamik wird durch stückweise differenzierbare Gleichungen bestimmt, die basierend auf der Winkelverschiebung $\theta$ umschalten, wobei Aufprallereignisse bei $\theta = 0$ auftreten.

Es werden zwei unterschiedliche Restitutionsmodelle unter identischen geometrischen und Erregungsbedingungen verglichen:

1. Housner-Modell ($e_H$): Die klassische Formulierung, bei der der COR ausschließlich von der Geometrie des Blocks (kritischer Winkel $\theta_c$) abhängt.
2. Mao et al. Modell ($e_M$): Eine alternative Formulierung, die auf einer vektorbasierten Analyse des Geschwindigkeitsfeldes basiert und die elastische Rückfederung sowie die Energiedissipation an der Kontaktstelle berücksichtigt.

Die Analyse verwendet drei primäre numerische Werkzeuge, um die Dynamik über verschiedene Schlankheitsgrade ($h/b = 1$ und $h/b = 2$) zu charakterisieren:

• Bifurkationsdiagramme: Konstruiert mittels Poincaré-Schnitte, um Übergänge zwischen periodischen, multiperiodischen und chaotischen Regimen bei variierender Erregungsamplitude ($\alpha$) zu identifizieren.
• Lyapunov-Exponenten: Der maximale Lyapunov-Exponent ($\lambda_{max}$) wird unter Verwendung des Benettin–Wolf-Algorithmus mit Saltationsmatrizen berechnet, um die Diskontinuitäten beim Aufprall zu handhaben und zwischen stabiler periodischer Bewegung ($\lambda_{max} < 0$) und chaotischer Dynamik ($\lambda_{max} > 0$) zu unterscheiden.
• Einzugsbereiche (Basins of Attraction): Der Zustandsraum wird diskretisiert, um Anfangsbedingungen den asymptotischen Ergebnissen zuzuordnen (gedämpfte Bewegung, symmetrische/asymmetrische Oszillationen, multiperiodische Antworten oder Umkippen). Dies ermöglicht eine quantitative Bewertung der Robustheit begrenzter Bewegungen gegenüber der Umkippanfälligkeit.

Kernergebnisse
Der Vergleich zeigt, dass die Wahl des Restitutionsmodells die vorhergesagte dynamische Landschaft signifikant verändert, wobei diese Effekte stark vom Schlankheitsgrad abhängen:

• Divergenz bei Blöcken mit geringem Schlankheitsgrad ($h/b = 1$): Die Formulierung von Mao et al. sagt einen höheren Restitutionskoeffizienten (geringere Energiedissipation) voraus als das Housner-Modell. Dies führt zu einem früheren Einsetzen komplexer oszillatorischer Verhaltensweisen und einem breiteren Parameterbereich, der mit chaotischer Dynamik assoziiert ist. Entscheidend ist, dass die Analyse der Einzugsbereiche einen drastischen Unterschied zeigt: Das Mao-Modell reduziert das Gebiet stabiler symmetrischer Oszillationen drastisch (von 55,55 % auf 19,64 % der Anfangsbedingungen) und erhöht die Häufigkeit des Umkippens signifikant (von 44,45 % auf 80,36 %). Die Übereinstimmung zwischen den beiden Modellen für diese Geometrie beträgt nur 64,09 %, was bedeutet, dass über ein Drittel der Anfangsbedingungen allein aufgrund des Aufprallgesetzes qualitativ unterschiedliche Ergebnisse liefert.
• Konvergenz bei Blöcken mit hohem Schlankheitsgrad ($h/b = 2$): Mit zunehmendem Schlankheitsgrad verringert sich die Diskrepanz zwischen den beiden Restitutionskoeffizienten. Folglich konvergieren die durch beide Modelle erzeugten dynamischen Merkmale progressiv. Obwohl das Mao-Modell auch innerhalb des begrenzten Bereichs multiperiodische Attraktoren einführt, gleichen sich die allgemeine Organisation des Phasenraums und der Anteil der Umkippereignisse zwischen beiden Formulierungen immer mehr an.
• Art der Übergänge: Das Mao-Modell fördert einen abrupten Verlust der regulären periodischen Bewegung und verschiebt Übergänge zu komplexen Regimen hin zu niedrigeren Erregungsamplituden. Im Gegensatz dazu zeigt das Housner-Modell einen graduelleren Übergang mit alternierenden Fenstern von Stabilität und Chaos vor dem Umkippen.

Bedeutung und Behauptungen
Die Arbeit postuliert, dass der Restitutionskoeffizient nicht bloß als lokaler Aufprallparameter betrachtet werden sollte, sondern als eine Modellannahme, die in der Lage ist, die globale dynamische Struktur von Kippsystemen grundlegend zu verändern. Die Studie demonstriert, dass alternative Restitutionsformulierungen – obwohl sie potenziell genauer bei der Vorhersage lokaler Aufprallreaktionen sind – im Vergleich zum klassischen Housner-Modell grundlegend andere Langzeitverhalten erzeugen können.

Insbesondere behaupten die Autoren, dass für Blöcke mit geringem Schlankheitsgrad die Anwendung der Mao et al. Formulierung eine reduzierte Robustheit sicherer Kippreaktionen suggeriert, wobei kleinere Störungen der Anfangsbedingungen eher zum Umkippen führen. Umgekehrt sinkt die Sensitivität gegenüber dem spezifischen Restitutionsmodell bei höheren Blöcken. Diese Ergebnisse liefern eine dynamische Perspektive darauf, warum Diskrepanzen zwischen klassischen Vorhersagen und experimentellen Beobachtungen bestehen, und betonen, dass die Behandlung der Aufprall-Dissipation ein kritischer Faktor bei der Bestimmung der Stabilitätsgrenzen und der Zugänglichkeit von Attraktoren in hybriden Kippsystemen ist. Die Arbeit kommt zu dem Schluss, dass das Mao-Modell zwar eine engere Annäherung an experimentelle Aufpralldaten bietet, die inhärenten idealisierten Starrkörperannahmen beider Modelle jedoch ihre Fähigkeit einschränken, komplexe Kontaktmechaniken wie lokale Deformation und Reibgleiten, wie sie in realen Struktursystemen vorkommen, zu erfassen.