Lektion
Die Theorie — Bode/Nyquist-Regelungslabor
Das Nyquist-Stabilitätskriterium entscheidet über die Stabilität eines Regelkreises durch Zählen statt durch Lösen. Man trägt den Frequenzgang der offenen Kette G(jω) in die komplexe Ebene ein, zählt, wie oft die Kurve den Punkt −1 + j0 umschlingt, und addiert die Zahl der instabilen Pole, die die Kette schon vor dem Schließen hatte: Z = N + P. Jedes Z außer null bedeutet, dass der geschlossene Kreis wegläuft.
Was die einzelnen Symbole bedeuten
P- die Zahl der Pole der offenen Kette
G(s), die bereits in der rechten Halbebene liegen — gezählt, bevor der Kreis geschlossen wird. Für die Strecke auf dieser Seite ist sie immer 0. N- die Zahl der im Uhrzeigersinn durchlaufenen Umschlingungen von
−1 + j0durch die Nyquist-Kurve. Umläufe gegen den Uhrzeigersinn zählen negativ,Nkann also unter null liegen. Z- die Zahl der Pole des geschlossenen Kreises in der rechten Halbebene — die Antwort, um die es eigentlich ging.
Z = 0heißt stabil; alles andere nicht. −1 + j0- der kritische Punkt, auf der Nyquist-Tafel rot eingezeichnet. Dort beträgt die Kreisverstärkung genau 1 bei genau 180° Phasennacheilung: ein einmal um den Kreis geschicktes Signal kehrt gleich groß und umgekehrt zurück.
Woher die Formel kommt
- Die Frage richtig stellen. Bei negativer Einheitsrückführung lautet die Übertragungsfunktion des geschlossenen Kreises
G/(1 + G); ihre Pole sind also die Nullstellen von1 + G(s). Stabilität schrumpft damit auf eine einzige Frage: Hat1 + G(s)Nullstellen in der rechten Halbebene? Sie zu finden hieße, ein Polynom zu faktorisieren, das man womöglich gar nicht besitzt. - Zählen statt suchen. Das Argumentprinzip besagt: Läuft
seinmal im Uhrzeigersinn um eine geschlossene Kontur, so umläuft das Bild von1 + G(s)den Ursprung genau so oft, wie darin Nullstellen liegen, abzüglich der darin liegenden Pole. Ein Umlauf um den Ursprung von1 + Gist dieselbe Bewegung wie ein Umlauf um−1vonG— daher liegt der kritische Punkt dort, wo er liegt, und heißt−1und nicht anders. - Als Kontur wählt man die gesamte rechte Halbebene: die imaginäre Achse hinauf und über einen unendlich großen Halbkreis zurück. Jetzt bedeuten „innerhalb der Kontur“ und „instabil“ dasselbe. Die eingeschlossenen Nullstellen sind
Z, die eingeschlossenen PoleP, die Umlaufzahl istN, und das Argumentprinzip lautetN = Z − P— umgestelltZ = N + P. Auf der imaginären Achse ists = jω, undG(jω)ist der Frequenzgang: genau die Kurve, die diese Seite ohnehin zeichnet. - Dieser letzte Satz ist der eigentliche Gewinn. Für das Zählen wird kein Modell gebraucht.
G(jω)lässt sich messen — echte Hardware mit Sinusschwingungen anregen und aufzeichnen, was zurückkommt —, sodass man ein Stabilitätsurteil über ein System bekommt, für das nie jemand eine Gleichung aufgeschrieben hat. Amplituden- und Phasenreserve sind dann zwei Arten zu fragen, wie weit die gemessene Kurve an−1vorbeigeht: einmal entlang der negativen reellen Achse, einmal auf dem Einheitskreis.
So liest du, was du siehst
Alle drei Tafeln zeigen die offene Kette — G(jω) für sich, mit noch nicht angeschlossenem Rückführzweig —, während die beiden Reserven darunter beschreiben, was nach dem Schließen passiert. Die Amplitudenüberhöhung nahe ωn und ebenso ωr und Mr in der Tabelle gehören daher zu G und nicht zu dem geschlossenen Kreis, den du entwirfst. Auf der Nyquist-Tafel ist die durchgezogene Kurve ω von null aufwärts, die gestrichelte ihr Spiegelbild für negative ω. Beide werden gebraucht: Eine Umschlingung entsteht erst, wenn sich die Kontur schließt, und ein Frequenzdurchlauf kann immer nur ihre Hälfte messen. Die Pfeile geben die Durchlaufrichtung an — sie entscheidet, ob ein Umlauf als im Uhrzeigersinn zählt.
- Setzt voraus
- Negative Einheitsrückführung und ein bekanntes
P. Diese Seite fragtPnie ab — sie setzt es auf null, was hier nur deshalb stimmt, weil der Dämpfungsregler beiζ = 0.05endet und beide Pole bei jeder erreichbaren Einstellung in der linken Halbebene bleiben. Gibt man demselben Kriterium eine Strecke, die offen bereits instabil ist, so bedeutet eine Kurve ohne Umschlingung genau das Gegenteil dessen, was sie hier bedeutet. Vorausgesetzt wird außerdem, dass unterwegs nichts gekürzt wurde: Ein durch eine Nullstelle gekürzter Pol der rechten Halbebene verschwindet ausG(s), bleibt aber in der Hardware — und die Zählung bescheinigt dann einem Kreis Stabilität, der trotzdem wegläuft. - Versagt, wenn
- Nyquist antwortet mit Ja oder Nein über ein einziges exaktes Modell, und Hardware ist weder exakt noch unveränderlich.
Z = 0sagt, dass der nominelle Kreis stabil ist, und überhaupt nichts darüber, wie weit eine Verstärkung driften oder eine Temperatur wandern muss, bis er es nicht mehr ist — deshalb druckt diese Seite zwei Reserven statt des einen Wortes „stabil“. Die Reserven selbst sind zu optimistisch, denn jede wird gemessen, während die andere festgehalten wird: Ein Kreis kann bequeme Amplituden- und bequeme Phasenreserve zeigen und trotzdem bei einer mäßigen Änderung beider zugleich mit der Kurve mitten durch−1laufen. Der ehrliche Wert ist der kürzeste Abstand der Kurve zum kritischen Punkt, in welcher Richtung auch immer, und den zeichnet hier keine Tafel.
Aufgabe vollständig gelöst
-
Eine Resonanzspitze von 1,0002 und eine unendliche Amplitudenreserve 6 Schritte
Ein Regelkreis zweiter Ordnung mit ζ = 0,7. Das Bedienfeld zeigt eine Resonanzüberhöhung von 1,0002 und einen unendlichen Amplitudenrand. Erklären Sie beides.
-
Die Standardform zweiter Ordnung mit den drei Parametern, die über die Schieberegler eingestellt werden. Alles Folgende wird am Nenner abgelesen.
-
Die Pole sind die Nullstellen. Der Realteil ist die Abklingrate und der Imaginärteil die Schwingungsfrequenz, und beide sind in der Sprungantwort sichtbar.
-
Die Resonanzfrequenz ist nicht die Eigenfrequenz. Die darin enthaltene Quadratwurzel kann zu Null werden; darum geht es im gesamten nächsten Schritt.
-
Sie wird bei ζ = 1/√2 zu Null. Bei 0,7 befindet sich das System knapp innerhalb des Bereichs mit Resonanzüberhöhung, und die Überhöhung beträgt 0,002 dB — algebraisch vorhanden und auf jedem Diagramm unsichtbar.
-
Nun zur Phase. Jeder Pol trägt zu einer Phasenverzögerung bei, die sich 90° nähert, sodass zwei Pole sich von oben her 180° nähern und diesen Wert nie schneiden.
-
Kein Schneiden bedeutet, dass es keinen nennenswerten Amplitudenrand gibt — und genau das zeigt das Bedienfeld an.
Antwort
ζ = 0,7 liegt knapp unter 1/√2, und zwei Pole können −180° nicht erreichen. Die Resonanzfrequenz enthält √(1 − 2ζ²), was bei ζ = 0,7071 zu Null wird — somit liegt 0,7 genau 0,001 innerhalb des Bereichs, in dem überhaupt eine Überhöhung existiert, und die vorhandene Überhöhung beträgt 1,0002 bzw. 0,002 dB. Dieser Schwellenwert ist die Butterworth-Bedingung, und das ist der Grund, warum 0,707 die Dämpfung ist, zu der jedes Filterlehrbuch greift: Es ist die reaktionsschnellste Einstellung, die noch keine Überhöhung aufweist. Der unendliche Amplitudenrand ist eine stärkere und gefährlichere Aussage. Jeder Pol trägt höchstens 90° zur Phasenverzögerung bei, sodass zwei von ihnen sich 180° nähern und diese nie erreichen; keine noch so große Verstärkung bringt diesen Regelkreis zum Schwingen. Reale Regelkreise schwingen letztlich immer — der dritte Pol, die Totzeit, die unmodellierte Resonanz —, sodass ein unendlicher Amplitudenrand ein zuverlässiges Zeichen dafür ist, dass dem Modell etwas fehlt, und nicht dafür, dass die Hardware sicher ist.
-
Quellen (2)
- The lesson — the encirclement count Z = N + P, in the paper that introduced it: H. Nyquist, "Regeneration Theory." Bell System Technical Journal 11(1), 126–147, 1932.
- Insight block 3 — phase margin as a delay budget: K. Ogata, Modern Control Engineering, 5th edition, ch. 7. Prentice Hall, 2010. ISBN 978-0-13-615673-4.