Stabilitätstheorie

Die mathematische Stabilitätstheorie beschäftigt sich mit der Entwicklung von Störungen, die als Abweichung von bestimmten Zuständen dynamischer Systeme auftreten. Ein solcher Zustand kann etwa eine Ruhelage oder ein bestimmter Orbit sein, z. B. ein periodischer Orbit. Ein System ist instabil, wenn eine kleine Störung zu großen und aufklingenden Abweichungen führt.

Neben ihrer theoretischen Bedeutung wird die Stabilitätstheorie in der Physik und in der Theoretischen Biologie angewendet sowie in technischen Gebieten, z. B. in der Technischen Mechanik oder der Regelungstechnik.

Die Lösungsansätze für die Probleme der Stabilitätstheorie sind gewöhnliche und partielle Differentialgleichungen.

Mathematische Stabilitätsbegriffe

Für die Charakterisierung der Stabilität der Ruhelage eines dynamischen Systems ${\dot {\vec {x}}}=f({\vec {x}})$ existieren mehrere Stabilitätsbegriffe mit jeweils etwas unterschiedlicher Aussage:

• Eine Ruhelage ${\vec {x}}_{R}$ heißt Ljapunow-stabil, wenn eine hinreichend kleine Störung auch stets klein bleibt. Präziser formuliert: Für jedes $\varepsilon >0$ existiert ein $\delta (\varepsilon )>0$ derart, dass für alle Zeiten $t\geq 0$ und alle Trajektorien ${\vec {x}}(t)$ mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \|\vec{x}(0) - \vec{x}_R\| < \delta(\varepsilon)} gilt: Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \|\vec{x}(t) - \vec{x}_R\| < \varepsilon} .
• Eine Ruhelage Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{x}_R} heißt attraktiv, wenn es ein Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \eta> 0} derart gibt, dass jede Trajektorie Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{x}(t)} mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \|\vec{x}(0)-\vec{x}_R\| < \eta} für alle Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle t\ge0} existiert und die folgende Grenzwertbedingung erfüllt: Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \lim_{t\to\infty}\vec{x}(t)=\vec{x}_R.}
• Eine Ruhelage heißt asymptotisch stabil, wenn sie Ljapunow-stabil und attraktiv ist.
• Eine Ruhelage heißt neutral stabil oder marginal stabil, wenn sie stabil, aber nicht asymptotisch stabil ist.

Für den Fall diskreter Systeme, die durch Differenzengleichungen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{x}_{k+1} = f(\vec{x}_k)} beschrieben werden, ist die Ruhelage gleichzeitig Fixpunkt der Rekursionsgleichung Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{x}_{k+1} = f(\vec{x}_k)} und es sind ähnliche Stabilitätsdefinitionen üblich.

Stabilitätsanalyse linearer zeitinvarianter Systeme

Datei:Polstellen in der s-ebene.png

Bedeutung der Pole und der konjugiert komplexen Polpaare in der linken und rechten s-Halbebene

Bei zeitkontinuierlichen linearen zeitinvarianten Systemen kann die Stabilität an der Übertragungsfunktion durch die Lage der Pole in der s-Ebene (Nennerpolynom der Laplace Übertragungsfunktion) abgelesen werden:

• Asymptotische Stabilität: wenn sämtliche Pole in der linken s-Halbebene liegen,
• Grenzstabilität: wenn kein Pol in der rechten s-Halbebene liegt und mindestens ein einfacher Pol, aber kein mehrfacher Pol, auf der imaginären Achse der s-Halbebene liegt,
• Instabilität: sonst (wenn mindestens ein Pol in der rechten s-Halbebene liegt oder wenn mindestens ein mehrfacher Pol auf der imaginären Achse der s-Ebene liegt).

Bei zeitdiskreten linearen zeitinvarianten Systemen kann die Stabilität durch die Lage der Pole in der z-Ebene (Nennerpolynom der z-Übertragungsfunktion) abgelesen werden.

• Asymptotische Stabilität: wenn sämtliche Pole im Einheitskreis liegen,
• Grenzstabilität: wenn mindestens ein Pol auf dem Einheitskreis liegt und alle anderen innerhalb,
• Instabilität: sonst (wenn mindestens ein Pol außerhalb des Einheitskreises in der z-Ebene liegt).

Achtung: der Begriff grenzstabil führt leicht zu Missverständnissen, da die Polstellenlage tatsächlich nach den meisten Stabilitätsdefinitionen instabile Systeme kennzeichnet.

Direkte Methode von Ljapunow und Ljapunow-Funktion

Ljapunow entwickelte 1883 die sogenannte Direkte oder Zweite Methode (die Erste Methode war die Linearisierung, siehe unten), um die oben genannten Stabilitätseigenschaften an konkreten Systemen zu überprüfen. Hierzu definiert man zunächst zu einem dynamischen System der Form Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot{\vec{x}} = f(\vec{x})} und einer reellwertigen differenzierbaren Funktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle V(\vec{x})} die orbitale Ableitung

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot{V}(\vec{x}):=\left\langle\operatorname{grad}\,V(\vec{x}),\dot{\vec{x}}\right\rangle=\left\langle\operatorname{grad}\,V(\vec{x}),f(\vec{x})\right\rangle} .

Eine reellwertige differenzierbare Funktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle V} heißt Ljapunow-Funktion (für das Vektorfeld Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle f} ), wenn Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot V(\vec{x})\le0} für alle Punkte ${\vec {x}}$ aus dem Phasenraum gilt. Eine Ljapunow-Funktion ist ein ziemlich starkes Hilfsmittel für einen Stabilitätsbeweis, wie die folgenden beiden Kriterien zeigen:

Erstes Kriterium von Ljapunow: Gegeben sei ein dynamisches System Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot{\vec{x}} = f(\vec{x})} . Gelten die Bedingungen

1. Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{x}_R} ist eine Ruhelage des Systems,
2. Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle V(\vec{x})} ist eine Ljapunow-Funktion für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle f} ,
3. Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle V(\vec{x})} besitzt an der Stelle Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{x}_R} ein striktes lokales Minimum,

dann ist die Ruhelage Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{x}_R} stabil.

Zweites Kriterium von Ljapunow: Gilt zusätzlich zu den Voraussetzungen des ersten Kriteriums noch

4. für ${\vec {x}}\neq {\vec {x}}_{R}$ in einer Umgebung der Ruhelage Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{x}_R} gilt Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot V(\vec{x})<0} ,

dann ist die Ruhelage asymptotisch stabil.

Die Verwendung einer Ljapunow-Funktion nennt man Direkte Methode, weil sich damit direkt aus dem Vektorfeld Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle f} ohne Kenntnis der Trajektorien (also ohne, dass man die Differentialgleichung lösen müsste) Aussagen über die Stabilität einer Ruhelage gewinnen lassen.

Ljapunowgleichung

Für den Fall linearer Systeme Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot{\vec{x}}=A\vec{x}} kann zum Beispiel immer eine positiv definite quadratische Form Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle v(\vec{x})=\vec{x}^T R \vec{x}} als Ljapunow-Funktion Verwendung finden. Sie erfüllt offensichtlich die obigen Bedingungen (1) und (2). Bedingung (3) führt auf die Ljapunow-Gleichung

$A^{T}R+RA=-Q$,

welche eine spezielle Form der Sylvester-Gleichung ist. Falls Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle Q} positiv definit ist, so ist Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle v(\vec{x})=\vec{x}^T R \vec{x}} eine Ljapunow-Funktion. Für stabile lineare Systeme lässt sich eine solche Funktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle v(\vec{x})=\vec{x}^T R \vec{x}} immer finden.

Stabilitätsanalyse linearer und nichtlinearer Systeme

Ein dynamisches System sei gegeben durch die Differentialgleichung Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot{\vec{x}} = f(\vec{x})} .

Wir betrachten eine Störung Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \delta = \vec{x}(t) - \vec{x}_R} zum Zeitpunkt Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle t} als Abweichung von der Ruhelage Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{x}_R} :

• wenn das System linear ist, kann diese Störung vollständig durch die Jacobi-Matrix Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \mathbf{J}} der ersten Ableitungen nach Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{x}} ausgedrückt werden.
• ist das System nichtlinear und die Störung klein genug, so kann man es „linearisieren“, d. h. die Funktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle f} nach $\delta$ um Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{x}_R} Taylor-entwickeln.

In beiden Fällen ergibt sich für die Zeitentwicklung von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \delta} :

${\dot {\delta }}=\mathbf {J} ({\vec {x}}_{R})\,\delta$

Diese Entwicklung wird demnach maßgeblich von den Eigenwerten der Jacobi-Matrix bestimmt. Konkret ergeben sich die folgenden drei Fälle:

• Der Realteil aller Eigenwerte der Jacobi-Matrix ist negativ. Dann fällt Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \delta} exponentiell ab, und die Ruhelage ist asymptotisch stabil.
• Der Realteil eines Eigenwertes der Jacobi-Matrix ist positiv. Dann wächst Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \delta} exponentiell an, und die Ruhelage ist instabil.
• Der größte Realteil aller Eigenwerte der Jacobi-Matrix ist Null. Dies bedeutet für ein lineares System:
  • falls für alle Eigenwerte mit verschwindendem Realteil die algebraische Vielfachheit gleich der geometrischen Vielfachheit ist: marginale Stabilität der Ruhelage.
  • sonst, d. h. falls nicht für alle Eigenwerte mit verschwindendem Realteil die algebraische Vielfachheit gleich der geometrischen Vielfachheit ist: Instabilität der Ruhelage.

Bei nichtlinearen Systemen, die nur um die Ruhelage linearisiert wurden, kann die Stabilität in diesem dritten Fall auch noch von Termen höherer Ordnung in der Taylorentwicklung bestimmt werden. In diesem Fall vermag die lineare Stabilitätstheorie keine Aussage zu machen.

Siehe auch Autonome Differentialgleichung.

Stabilitätsgefährdung im Bauwesen

→ Hauptartikel: Stabilitätsgefährdung

Im Bauwesen müssen druckbeanspruchte Stäbe auf Stabilitätsgefährdung (i. d. R. Knicken) geprüft werden und gegebenenfalls nach Theorie II. Ordnung nachgewiesen werden. Man braucht die Theorie II. Ordnung um Stabilitätsgefährdung beschreiben zu können. Im Stahlbau, im (Stahl-)Betonbau als auch im Holzbau sind laut aktueller Normung stabilitätsgefährdete Stäbe auf Knicken nachzuweisen.

Beispiele

Ein untersuchter Verformungszustand der Festigkeitslehre oder ein Bewegungszustand der Dynamik können ab einer zu bestimmenden Stabilitätsgrenze in einen anderen Zustand wechseln. Damit verbunden sind in der Regel nichtlinear ansteigende Verformungen oder Bewegungen, die zur Zerstörung von Tragwerken führen können. Um diese zu vermeiden, ist die Kenntnis der Stabilitätsgrenze ein wichtiges Kriterium zur Bemessung von Bauteilen.

Weitere Beispiele:

• Stabilitätsprüfung von Regelkreisen,
• Dynamik von Insektenpopulationen,
• Eulerscher Knickstab,
• Kippen von schlanken Trägern,
• Beulen von Platten und Schalen,
• Wachstum kleiner Störungen in einer Grenzschicht, die zum laminar-turbulenten Umschlag führen.
• Raumsonden oder (kleine) Himmelskörper auf Bahnen um Lagrange-Punkte eines Planets, der um Zentralgestirn kreist

Siehe auch

• Eigenbewegung (Regelungstechnik)
• Lineare Stabilitätstheorie zur Stabilitätstheorie in der Strömungslehre

Literatur

• Herbert Amann: Gewöhnliche Differentialgleichungen. 2. Auflage. de Gruyter, Berlin 1995, ISBN 3-11-014582-0. 
• W. Hahn: Stability of Motion. Springer, 1967.
• N. Rouche, P. Habets und M. Laloy: Stability Theory by Liapunov's Direct Method. Springer, 1977.
• Gerald Teschl: Ordinary Differential Equations and Dynamical Systems. American Mathematical Society, Providence 2012, ISBN 978-0-8218-8328-0 (freie Onlineversion). 

Weblinks

Wikiversity

Wikiversity: Eine Einführung in Stabilitätsbegriffe für Matrizen – Kursmaterialien

Commons: Stability theory – Sammlung von Bildern, Videos und Audiodateien

• Stability. In: Scholarpedia. (englisch, inkl. Literaturangaben)