www.wikidata.de-de.nina.az
Die Methode der harmonischen Balance angewandt in der Regelungstechnik besteht darin die Dauerschwingungen eines nichtlinearen zuruckgekoppelten dynamischen Systems durch eine harmonische Schwingung anzunahern Dabei wird das nichtlineare Teilsystem durch eine Beschreibungsfunktion dargestellt Mit der Methode werden die Parameter des Schwingungszustandes und die Stabilitatsgrenze berechenbar N M Krylow und N N Bogoljubow verwendeten bereits 1937 den Begriff Harmonische Balance fur dieses von ihnen entwickelte Verfahren Ein nichtlineares dynamisches System kann nach dem Hammerstein Modell in ein statisches nichtlineares System und ein dynamisches lineares System zerlegt werden Wird der Ausgang des linearen Systems negativ auf den Systemeingang des nichtlinearen Systems zuruckgefuhrt und damit zu einem Regelkreis geschaltet kann das Gesamtsystem schwingen Mit dem Verfahren der Harmonischen Linearisierung wird von einem schwingenden nichtlinearen Regelkreis ausgegangen dessen Ausgangssignal durch das Tiefpass Verhalten der Regelstrecke eine angenaherte harmonische Schwingung ausfuhrt die in den Systemeingang negativ zuruckgefuhrt wird Nur unter diesen Bedingungen darf das nichtlineare statische System als lineares Ubertragungsglied mit der Beschreibungsfunktion N A displaystyle N A definiert werden welche nur von der sinusformigen harmonischen Eingangsschwingung mit der Amplitude A displaystyle A und nicht von der komplexen Frequenz j w displaystyle mathrm j omega abhangt Mit der Gleichung der Harmonischen Balance werden die Beziehungen der Beschreibungsfunktion N A displaystyle N A des statischen nichtlinearen Systems und des dynamischen linearen Systems G 0 j w displaystyle G 0 mathrm j omega in ein Verhaltnis gesetzt Daraus lassen sich die zwei kritischen Systemgrossen des harmonisch schwingenden Regelkreises die Eingangsamplitude A displaystyle A und die kritische Frequenz w k r i t displaystyle omega mathrm krit an der Stabilitatsgrenze errechnen oder grafisch nach dem Zwei Ortskurven Verfahren bestimmen Die Anwendung der Harmonischen Balance zur Prufung von nichtlinearen Regelkreisen mit dem anschaulichen Zwei Ortskurven Verfahren darauf wann Dauerschwingungen auftreten und wie diese sich vermeiden lassen erfordert keine besonderen mathematischen Kenntnisse Die benotigte Beschreibungsfunktion des nichtlinearen statischen Systems zur Konstruktion der Ortskurve ist in vielen Varianten in der Fachliteratur der Regelungstechnik dargestellt Dies gilt insbesondere fur die verschiedenen Formen der Kennlinienregler Inhaltsverzeichnis 1 Grundlagen des linearen und nichtlinearen schwingenden Systems 1 1 Ruhelagen und Stabilitat 1 2 Stabilitatsgrenze des linearen dynamischen Systems 1 3 Konstruktion einer Ortskurve des Frequenzgangs 1 4 Verhalten des nichtlinearen dynamischen Systems 2 Harmonische Linearisierung 2 1 Festlegung der Signalgrossen Parameter und Eckwerte der Kennlinien Strukturen 2 2 Methode der harmonischen Linearisierung 2 2 1 Vereinfachung der nichtlinearen Beschreibungsfunktion im Zustand der Dauerschwingung 2 3 Beschreibungsfunktion des nichtlinearen Systems 2 4 Aufstellung bekannter Beschreibungsfunktionen von Kennlinien Elementen 3 Gleichung der Harmonischen Balance 3 1 Formelmassige Losung der Gleichung der Harmonischen Balance 3 2 Zwei Ortskurven Verfahren der Harmonischen Balance 3 3 Stabilitat der Grenzschwingungen 4 Vermeidung von Dauerschwingungen 4 1 Kompensation trager PT1 Glieder durch PD1 Glieder 4 2 Beispiel eines schwingungsfreien Regelkreises mit einem Dreipunktregler 5 Fazit zur Anwendung der Harmonischen Balance 6 Literatur 7 EinzelnachweiseGrundlagen des linearen und nichtlinearen schwingenden Systems BearbeitenStreng genommen sind fast alle Regelkreise nichtlineare Systeme Die haufigste statische Nichtlinearitat ist die Sattigungseigenschaft der Stellgrosse des Reglers So ist bei einem P Regler trotz erhohter Regelabweichung keine Signalerhohung der Stellgrosse moglich weil eine Signalbegrenzung vorliegt Fur diesen einfachen Fall der Signalbegrenzung wird die Regelgrosse bei einem Fuhrungsgrossensprung oder einer plotzlich auftretenden Storgrosse der Sollwert der Regelgrosse nicht so schnell erreicht wie ohne Stellgrossen Begrenzung Wird die P Verstarkung eines stabilen Regelkreises weiter erhoht dann wird vorausgesetzt die Regelstrecke ist ein lineares System hoherer als 2 Ordnung der Regelkreis an der Stabilitatsgrenze f P Verstarkung instabil Die Regelgrosse wird dann mit konstanter oder zunehmender Amplitude schwingen Dabei ist es nicht gleichgultig ob eine Stellgrossen Begrenzung vorliegt oder nicht Dauerschwingungen auch Grenzschwingungen und Grenzzyklen sind periodische sinusformige Zeitvorgange im eingeschwungenen Zustand mit konstanter Amplitude Sie unterscheiden sich von auf und abklingenden Schwingungen Ruhelagen und Stabilitat Bearbeiten nbsp Modell der System Ruhelagen Zustand 1 ist gegenuber kleinen Storungen stabil und geht bei grossen Storungen in Zustand 3 uber Zustand 2 ist labil Eine Ruhelage eines dynamischen Systems im schwingungsfreien Zustand ist asymptotisch stabil wenn nach einer Signalstorung die Ausgangsgrosse y t displaystyle y t nbsp wieder in die Ruhelage zuruckkehrt Ein stabiles System neigt dazu seinen momentanen Zustand beizubehalten auch wenn Storungen von aussen einwirken Zu jedem Stabilitatsbereich gehort eine Ruhelage und ein Einzugsbereich der Ruhelage Hat ein nichtlinearer Regelkreis nur eine Ruhelage so kann man annehmen dass seine Ausgangsgrosse y t displaystyle y t nbsp global asymptotisch stabil ist sofern keine Dauerschwingungen auftreten Nichtlineare dynamische Systeme sind im Gegensatz zu linearen Systemen von Anfangsbedingungen und des Eingangssignals der statischen Nichtlinearitat abhangig Je nach Funktion der statischen Nichtlinearitat kann stabiles und instabiles Verhalten auftreten Nichtlineare Systeme konnen bei Einsatz von sogenannten Kennlinien Reglern im Regelkreis gewollte Schwingungen hervorrufen Solche stationaren Dauerschwingungen werden auch als Grenzschwingungen bezeichnet Grenzschwingungen konnen stabil instabil und semistabil sein Nichtlineare Regelsysteme konnen im Gegensatz zu linearen Systemen mehrere Ruhelagen haben Ein Pendel wird mit seinen Massen der potentiellen und kinetischen Energien durch eine Differenzialgleichung 2 Ordnung beschrieben in der der Winkel f displaystyle varphi nbsp als 2 Ableitung und nichtlinear auftritt Bewertet man den Winkel des Befestigungspunktes mit 90 dann hat das Pendel eine Ruhelage bei 270 und einen Einzugsbereich der Ruhelage gt 90 und lt 90 und eine instabile Ruhelage von 90 im Koordinatensystem Ein nichtlineares dynamisches System mit der Ausgangsgrosse y t displaystyle y t nbsp hat folgendes Verhalten bei Auslenkung aus einer instabilen Ruhelage y t displaystyle y t nbsp geht in eine andere Ruhelage uber y t displaystyle y t nbsp strebt in Richtung displaystyle infty nbsp y t displaystyle y t nbsp macht eine Dauerschwingung Mit der Harmonischen Balance lasst sich klaren ob Grenzschwingungen auftreten konnen und ob sie stabil oder instabil sind und welche Frequenz und Eingangs Amplitude sie haben Stabilitatsgrenze des linearen dynamischen Systems Bearbeiten Ein dynamisches System verhalt sich linear wenn die Wirkungen zweier linear uberlagerter Eingangssignale sich am Ausgang des Systems in gleicher Weise linear uberlagern Das System verhalt sich linear wenn es das Superpositionsprinzip Siehe auch Superpositionsprinzip der Physik und das Verstarkungsprinzip erfullt Ein Ubertragungssystem ist intern stabil wenn alle Teil Ubertragungsfunktionen nur Pole in der linken s Halbebene haben Ein Ubertragungssystem gilt als extern stabil wenn jedes beliebige beschrankte Eingangssignal an dem System auch ein beschranktes Ausgangssignal hervorruft Siehe BIBO Stabilitat nbsp Blockschaltbild eines einfachen StandardregelkreisesDie Stabilitatsgrenze eines linearen Regelkreises ist nach dem vereinfachten Stabilitatskriterium von Nyquist erreicht wenn die Ortskurve des Frequenzgangs Siehe auch PT2 Glied des offenen Regelkreises genau den kritischen Punkt 1 der Abszisse des Realteils schneidet Fur Werte des Schnittpunktes Re lt 1 displaystyle operatorname Re lt 1 nbsp entstehen aufklingende Schwingungen des geschlossenen Regelkreises Damit gilt fur die Stabilitat der Grenzschwingung der kritischen Frequenz w k r i t displaystyle omega mathrm krit nbsp des offenen Regelkreises F 0 j w k r i t 1 displaystyle F 0 mathrm j omega mathrm krit 1 nbsp Zur Ermittlung der Ortskurve wird der komplexe Frequenzgang in Real und Imaginarteil zerlegt und fur verschiedene Werte der Frequenz w displaystyle omega nbsp in der Gaussschen Zahlenebene eingetragen Aus einer gegebenen Ubertragungsfunktion G s displaystyle G s nbsp oder einem Frequenzgang G j w displaystyle G mathrm j omega nbsp in Produktdarstellung wird der Nenner des Frequenzgangs als Polynom ausmultipliziert und in Realteil und Imaginarteil geordnet Damit der imaginare Anteil im Nenner des Frequenzgangs verschwindet muss der Real und Imaginarteil des Nenners mit sich selbst aber mit umgekehrten Vorzeichen mit der Frequenzganggleichung multipliziert werden Beispiel x j w x j w x 2 j 2 w 2 x 2 w 2 displaystyle x mathrm j omega cdot x mathrm j omega x 2 mathrm j 2 omega 2 x 2 omega 2 nbsp Anmerkung j w 2 w 2 j w 3 j w 3 displaystyle mathrm j omega 2 omega 2 mathrm j omega 3 mathrm j omega 3 nbsp Der Zahler wird mit dem Real und Imaginarteil des Nenners multipliziert damit die Gleichung mit der ursprunglichen Aussage unverandert ist Der Zahler wird nach Realteil und Imaginarteil geordnet nbsp Ortskurve des FrequenzgangsKonstruktion einer Ortskurve des Frequenzgangs Bearbeiten Gegeben Ubertragungsfunktion 3 Ordnung mit globalem I Verhalten und der Verstarkung K displaystyle K nbsp G s K s T 1 s 1 T 2 s 1 displaystyle G s frac K s cdot T 1 cdot s 1 T 2 cdot s 1 nbsp G j w K j w T 1 j w 1 T 2 j w 1 K T 1 T 2 w 2 j w T 1 T 2 w 3 displaystyle G mathrm j omega frac K mathrm j omega cdot T 1 cdot mathrm j omega 1 T 2 cdot mathrm j omega 1 frac K T 1 T 2 cdot omega 2 mathrm j cdot left omega T 1 cdot T 2 cdot omega 3 right nbsp Wird der Realteil des Nenners mit der Hilfsgrosse a T 1 T 2 w 2 displaystyle a T 1 T 2 omega 2 nbsp und der Imaginarteil des Nenners mit der Hilfsgrosse b w T 1 T 2 w 3 displaystyle b omega T 1 T 2 omega 3 nbsp bezeichnet kann der Frequenzgang G j w displaystyle G mathrm j omega nbsp fur den Realteil und den Imaginarteil wie folgt berechnet werden Zunachst werden die Zahlenwerte der Hilfsgrossen a displaystyle a nbsp und b displaystyle b nbsp als Funktion von w 0 displaystyle omega 0 nbsp bis w displaystyle omega infty nbsp in Stufen von w displaystyle omega nbsp Werten errechnet und in die unten stehende Gleichung des Frequenzgangs G j w displaystyle G mathrm j omega nbsp eingesetzt Die so errechneten realen und imaginaren Anteile des Frequenzgangs werden in das Diagramm der Gaussschen Zahlenebene als Ortskurve des Frequenzgangs eingetragen Ortskurve des Frequenzgangs G j w displaystyle G j omega nbsp G j w K a a 2 b 2 j K b a 2 b 2 displaystyle G mathrm j omega frac K cdot a a 2 b 2 mathrm j frac K cdot b a 2 b 2 nbsp nbsp Sprungantwort vor und an der Stabilitatsgrenze Zahlenbeispiel der Stabilitatsgrenze eines linearen Regelkreises Gegeben Offener Regelkreis mit einem P Regler K 8 displaystyle K 8 nbsp und einer Regelstrecke 3 Ordnung T 1 displaystyle T 1 nbsp G 0 s 8 T s 1 3 displaystyle G 0 s frac 8 T cdot s 1 3 nbsp Gesucht Verhalten der Sprungantwort fur die Fuhrungsgrosse w t 1 displaystyle w t 1 nbsp des geschlossenen Regelkreises Nach dem Schliessen des Regelkreises mit dieser ungunstigen Regelstrecke mit drei gleichen Zeitkonstanten wird bereits bei der Verstarkung K 8 displaystyle K 8 nbsp die Stabilitatsgrenze fur eine konstante Amplitude der Regelgrosse erreicht Die Ortskurve durchlauft als f w 2 p f displaystyle f omega 2 pi f nbsp beginnend bei f w 0 displaystyle f omega 0 nbsp mit dem Wert Re 8 w 0 Im 0 displaystyle operatorname Re 8 omega 0 operatorname Im 0 nbsp den 4 3 und 2 Quadranten und schneidet die Abszisse des Realteils bei genau Re 1 Im 0 displaystyle operatorname Re 1 operatorname Im 0 nbsp Die Form der Ortskurve andert sich nicht solange die drei Zeitkonstanten einen beliebigen aber untereinander den gleichen Wert T 1 T 2 T 3 T displaystyle T 1 T 2 T 3 T nbsp haben Die Grosse der Zeitkonstante T displaystyle T nbsp hat dann eine Bedeutung wenn das Zeitverhalten des Regelkreises z B als Sprungantwort der Fuhrungsgrosse dargestellt werden soll Drei gleiche Zeitkonstanten in einer Regelstrecke mit drei Verzogerungsgliedern stellen regelungstechnisch das ungunstigste Zeitverhalten bezuglich der Stabilitat des Regelkreises dar Wird die Verstarkung K gt 8 displaystyle K gt 8 nbsp erhoht dann treten progressiv zunehmend grossere Amplituden mit steigender Zeit auf Die Ortskurve des Frequenzgangs schneidet die Abszisse bei Werten lt 1 Wird die Verstarkung K lt 8 displaystyle K lt 8 nbsp reduziert dann nehmen die Amplituden der Regelgrosse einen aperiodischen gedampften Verlauf und die Regelgrosse nimmt nach der Einschwingzeit einen konstanten Wert an Die Ortskurve schneidet die Abszisse bei Werten gt 1 Besteht die Regelstrecke nur aus einem Verzogerungssystem 2 Ordnung mit einem P Regler beliebig hoher Verstarkung kann es theoretisch nicht zu Dauerschwingungen kommen weil die Ortskurve des Frequenzgangs nur den 4 und 3 Quadranten der gaussschen Zahlenebene durchlauft und den kritischen Punkt Re 1 displaystyle operatorname Re 1 nbsp nicht treffen kann Anmerkung Fur die Berechnung des Schwingungszustandes eines nichtlinearen dynamischen Systems nach dem Zwei Ortskurven Verfahren ist die Kenntnis der Konstruktion der Ortskurve des Frequenzgangs des linearen Systems erforderlich Verhalten des nichtlinearen dynamischen Systems Bearbeiten nbsp Allgemeine Darstellung des nichtlinearen RegelkreisesDynamische Systeme werden durch Differenzialgleichungen DGL beschrieben Enthalt die DGL oder die gesuchte Funktion oder deren Ableitungen eine Potenz n gt 1 displaystyle n gt 1 nbsp oder Produkte der gesuchten Funktion oder Winkelfunktionen Logarithmen usw in den Argumenten handelt es sich um eine nichtlineare DGL und das System verhalt sich nichtlinear Nichtlineare Differenzialgleichungen sind nur in sehr seltenen Ausnahmefallen analytisch losbar Sie konnen aber mittels der numerischen zeitdiskreten Methoden oder mit kommerziellen Rechenprogrammen Simulink relativ einfach berechnet werden Bei einem nichtlinearen statischen System kann die Ubertragungskennlinie als eine in weiten Grenzen sich andernde Verstarkung angesehen werden Das System antwortet auf ein sinusformiges Eingangssignal meist mit einer Verzerrung der Schwingung des Eingangssignals gleicher Frequenz und einer Phasenverschiebung Lediglich ein symmetrisches Zweipunkt Element Zweipunktregler ohne Hysterese ohne Totzone und Signalbegrenzung antwortet auf eine sinusformige Eingangsschwingung mit einer Rechteckschwingung der gewahlten Rechteck Amplitude u t U m a x d displaystyle u t pm U mathrm max pm d nbsp der gleichen Frequenz ohne eine Phasenverschiebung d displaystyle d nbsp ist der Kennlinien Hilfspunkt der nichtlinearen Kennlinie zu einer sinusformigen Eingangs Schwingung e t displaystyle e t nbsp Das einfache symmetrische Zweipunkt Element mit der Rechteck Amplitudenhohe d displaystyle pm d nbsp definiert sich als Signumfunktion mit u t displaystyle u t nbsp als Ausgangsgrosse u t f e t d sign e t d fur e t gt 0 stabil 0 fur e t 0 instabil d fur e t lt 0 stabil displaystyle u t f e t d cdot operatorname sign e t begin cases d amp text fur e t gt 0 text stabil 0 amp text fur e t 0 text instabil d amp text fur e t lt 0 text stabil end cases nbsp Das Zweipunkt Element hat zwei Ruhelagen Es erzeugt in Abhangigkeit von dem Eingangssignals e t displaystyle e t nbsp zwei Stellgrossen Der Zustand u t 0 displaystyle u t 0 nbsp fur e t 0 displaystyle e t 0 nbsp ist nur dann ein stationarer Wert wenn dem Zweipunkt Element eine Totzone T Z c displaystyle mathit TZ c nbsp c displaystyle c nbsp ist Kennlinien Hilfspunkt eingerichtet wird Damit wurde aus dem Zweipunkt Element ein Dreipunkt Element entstehen Sonderfalle sind statische nichtlineare Systeme mit quadratischer oder exponentieller Kennlinie Ein sinusformig erregtes System mit quadratischer Kennlinie antwortet mit einer Sinusschwingung der doppelten Frequenz und Gleichanteil verschobener Arbeitspunkt ein System mit exponentieller Kennlinie antwortet mit einem sinusahnlichen Impuls mit grossem Oberwellenanteil Wird ein nichtlineares dynamisches System 2 und hoherer Ordnung durch ein sinusformiges Eingangssignal e t displaystyle e t nbsp angeregt antwortet es meistens mit einer angenahert sinusformigen Schwingung gleicher Frequenz aber anderer Amplitude und einer Phasenverschiebung Ist das statische nichtlineare System ein schaltender Regler so wird in einem Regelkreis trotz steiler Impulsflanken des rechteckformigen Regler Ausgangssignals U m a x displaystyle pm U mathrm max nbsp bei genugend grosser Frequenz von e t displaystyle e t nbsp und u t displaystyle u t nbsp durch das Tiefpassverhalten des linearen dynamischen Systems ein angenahert sinusformiges Ausgangssignal y t displaystyle y t nbsp erreicht Durch dieses Verhalten ist die Harmonische Linearisierung begrundet die Rechteckschwingung mit ihrem grossen Oberwellen Anteil als Eingangsgrosse u t displaystyle u t nbsp wird nur in ihrer Grundschwingung entsprechender Amplitude berucksichtigt Damit kann das nichtlineare statische System N displaystyle N nbsp in Abhangigkeit von der Amplitude e m a x A displaystyle e mathrm max A nbsp der Sinusschwingung des Eingangssignals e t displaystyle e t nbsp beschrieben werden Nichtlineare Systeme sind haufig einzigartig jedoch konnen aus dem nichtlinearen statischen System Verhaltensweisen definiert werden die bei gegebenen Parametern eine System Beschreibungsfunktion erlauben Darunter fallen Elemente mit Begrenzung Totzone Hysterese progressiver Kennlinie und die sogenannten Kennlinienregler wie Zweipunkt und Mehrpunktregler die sich in zahlreichen Varianten wie folgt darstellen Zweipunkt Element Zweipunktregler symmetrisch und unsymmetrisch mit und ohne Hysterese Dreipunkt Element Dreipunktregler symmetrisch und unsymmetrisch mit und ohne Hysterese Mehrpunkt Element Element mit progressiver Kennlinie symmetrisch und unsymmetrisch Element mit degressiver Kennlinie symmetrisch und unsymmetrisch Elemente mit Begrenzung symmetrisch und unsymmetrisch Elemente mit Totzone symmetrisch und unsymmetrisch Elemente mit Totzone und Begrenzung Elemente mit Vorspannung Gleichrichter Element BetragselementLaut der Fachliteratur konnen bis zu 40 unterschiedliche Beschreibungsfunktionen nichtlinearer statischer Systeme entnommen werden die das nichtlineare Systemverhalten N A displaystyle N A nbsp nur fur ein periodisches sinusformiges Eingangssignal e t e max sin w t displaystyle e t e text max sin omega t nbsp mit der Amplitude A e m a x displaystyle A e mathrm max nbsp und das Ausgangssignal u t displaystyle u t nbsp nur mit der Grundwelle beschreiben Harmonische Linearisierung BearbeitenFestlegung der Signalgrossen Parameter und Eckwerte der Kennlinien Strukturen Bearbeiten In der bekannten Fachliteratur der Regelungstechnik sowie in den aktuellen Vorlesungsmanuskripten deutscher Hochschulen findet man fur die Harmonische Linearisierung kaum eine identische mathematische Ableitung aus der Fourieranalyse noch identische System und Signalbezeichnungen Ursache sind die verschiedenen Signalbezeichnungen z B mit Xe Xa E U X Y der Gleichungsaufbau der Fourieranalyse und die Bezeichnungen der geometrischen Eckpunkte der nichtlinearen Kennlinien Die Gleichung der Fourieranalyse und deren Vereinfachung unterscheiden sich durch die unterschiedliche Verwendung der harmonischen Sinus oder Kosinusschwingung mit und ohne Gleichanteile unterschiedliche Integrationsgrenzen und unterschiedliche Berechnung der Koeffizienten Auch sind die Ergebnisse der Beschreibungsfunktionen unter verschiedenen Namen ofter N A displaystyle N A nbsp N U 0 displaystyle N U 0 nbsp B X e displaystyle B hat X e nbsp N X e displaystyle N hat X e nbsp als Funktion der Amplitude der Eingangsgrosse des nichtlinearen Systems mit zahlreichen Kennlinien Regler Funktionen nicht ubereinstimmend weil die geometrischen Eckwerte der Kennlinien Beschreibung mit den Hilfsgrossen meist a b c d displaystyle a b c d nbsp in unterschiedlicher Weise benutzt werden Trotz der unterschiedlichen mathematischen Wege die zu der Gleichung der Harmonischen Balance fuhren ist das Ergebnis uberall identisch Das Verstandnis der Ortskurve des Frequenzgangs des linearen Systems ist erforderlich und ebenso das Verstandnis der Ortskurve der Beschreibungsfunktion des nichtlinearen statischen Systems Aus der Fachliteratur konnen die zahlreichen Beschreibungsfunktionen der nichtlinearen statischen Systeme mit den zugehorigen Skizzen der Kennlinien und der Ortskurven entnommen werden Die unterschiedlichen geometrischen Eckwerte mit der Kennlinienskizze a b c d displaystyle a b c d nbsp der Beschreibungsfunktion konnen leicht den verschiedenen Quellen der Fachliteratur zugeordnet werden Nachfolgend wird von den international verwendeten Systemgrossen des Regelkreises ausgegangen Fuhrungsgrosse w t 0 displaystyle w t 0 nbsp Reglereingangsgrosse e t displaystyle e t nbsp Reglerausgangsgrosse Nichtlinearitat u t displaystyle u t nbsp und linearer Regelstrecken Ausgang y t displaystyle y t nbsp Die folgenden dargestellten mathematischen Ableitungen sind der Fachliteratur 1 nachempfunden Methode der harmonischen Linearisierung Bearbeiten Die Methode der Harmonischen Linearisierung geht gemass dem Hammerstein Modell von der Annahme aus dass die von dem nichtlinearen System erzeugten Oberwellen durch das nachfolgende lineare System mit Tiefpassverhalten als unwirksam gefiltert werden Aus der Sicht der Regelungstechnik wird fur ein nichtlineares dynamisches System haufig das Hammerstein Modell im Gegensatz zum Wiener Modell eingesetzt weil in vielen Fallen die zum Regler gehorende Stellgrosse nichtlinear ist Das Verhalten des Reglers als nichtlineares statisches System N displaystyle N nbsp kann durch eine in weiten Grenzen sich andernde Verstarkung angesehen werden wodurch bei einem sinusformigen Eingangssignal e t e m a x sin w t displaystyle e t e mathrm max sin omega t nbsp eine Signalverzerrung und Phasenverschiebung f displaystyle varphi nbsp des Ausgangssignals u t f e t N e m a x displaystyle u t f e t N e mathrm max nbsp hervorgerufen wird Haufig ist die Abhangigkeit von N e m a x displaystyle N e mathrm max nbsp durch bekannte Parameter oder durch die Grosse der Amplitude e m a x displaystyle e mathrm max nbsp des sinusformigen Eingangssignals e t displaystyle e t nbsp bestimmt deshalb ist das Grosssignal Verhalten des Systems zu untersuchen Stabilitat im Grossen Manche Kennlinienregler sprechen erst an wenn das Eingangssignal eine Mindestgrosse erreicht Fur eindeutige nichtlineare Kennlinien wie das Zweipunkt Element d h fur jedes statische Eingangssignal e t displaystyle e t nbsp gibt es nur ein statisches Ausgangssignal u t displaystyle u t nbsp ist die Beschreibungsfunktion reell d h die Phasenverschiebung f 0 displaystyle varphi 0 nbsp nbsp Darstellung des nichtlinearen Regelkreises im Zustand der Harmonischen BalanceVereinfachung der nichtlinearen Beschreibungsfunktion im Zustand der Dauerschwingung Bearbeiten Die Eingangsgrosse e t displaystyle e t nbsp des statischen nichtlinearen Systems soll eine Sinusschwingung ohne Gleichanteil sein e t e m a x sin w t displaystyle e t e mathrm max cdot sin omega cdot t nbsp Die Amplitude e m a x displaystyle e mathrm max nbsp der Eingangs Sinusschwingung wird haufig in der Fachliteratur mit e m a x A displaystyle e mathrm max A nbsp bezeichnet Die verzerrte Ausgangsschwingung u t displaystyle u t nbsp lasst sich als Fourier Reihe schreiben als eine Summe aus einem Gleichwert a 0 displaystyle a 0 nbsp und aus den Sinussignalen der Grundfrequenz a 1 sin w 0 t ϕ 1 displaystyle a 1 sin omega 0 t phi 1 nbsp und den harmonischen Oberschwingungen als Vielfache der Grundfrequenz u t a 0 n 1 a n cos n 2 p t T b n sin n 2 p t T mit w 2 p T 2 p f displaystyle u t a 0 sum n 1 infty left a n cdot cos left frac n2 pi t T right b n cdot sin left frac n2 pi t T right right quad text mit quad omega frac 2 pi T 2 pi f nbsp Die Ausgangsgrosse u t displaystyle u t nbsp kann je nach Verhalten des nichtlinearen Systems einen Gleichanteil a 0 displaystyle a 0 nbsp enthalten Es interessiert aber nur das Schwingverhalten Wenn das nachgeschaltete lineare System G j w displaystyle G mathrm j omega nbsp die verzerrte sinusahnliche Eingangsschwingung u t displaystyle u t nbsp als Tiefpassfilter stark genug dampft werden die hohen Frequenzen im Vergleich zu den niedrigeren sehr stark unterdruckt Fur diese Betrachtung werden die Anteile der Oberschwingungen als vernachlassigbar angenommen u t a 0 a 1 cos w t b 1 sin w t displaystyle u t approx a 0 a 1 cos omega cdot t b 1 sin omega cdot t nbsp Die Beschreibungsfunktion des statischen nichtlinearen Systems ergibt sich durch das Verhaltnis der Grundwellen des Ausgangssignals u 1 t displaystyle u 1 t nbsp zum Eingangssignal e 1 t displaystyle e 1 t nbsp Sie berucksichtigt nur die Grundschwingung des Ausgangssignals ohne den Gleichanteil a 0 displaystyle a 0 nbsp Eingangssignal e 1 t e m a x sin w t displaystyle e 1 t e mathrm max cdot sin omega cdot t nbsp Ausgangssignal u 1 t a 1 cos w t b 1 sin w t displaystyle u 1 t a 1 cos omega cdot t b 1 sin omega cdot t nbsp Mit der komplexen Darstellung der beiden Signale ergibt sich die Beschreibungsfunktion nach der Euler Formel cos w t j sin w t e j w t displaystyle cos omega cdot t mathrm j cdot sin omega cdot t mathrm e mathrm j omega t nbsp cos w t j sin w t e j p 2 j e j w t displaystyle cos omega cdot t mathrm j cdot sin omega cdot t cdot mathrm e mathrm j pi 2 mathrm j cdot mathrm e mathrm j omega t nbsp e 1 t Re j e m a x e j w t e 1 t j e m a x e j w t displaystyle e 1 t operatorname Re mathrm j cdot e mathrm max cdot mathrm e mathrm j omega t Rightarrow tilde e 1 t mathrm j cdot e mathrm max cdot mathrm e mathrm j omega t nbsp u 1 t Re a 1 j b 1 e j w t u 1 t a 1 j b 1 e j w t displaystyle u 1 t operatorname Re a 1 mathrm j cdot b 1 cdot mathrm e mathrm j omega t Rightarrow tilde u 1 t a 1 mathrm j cdot b 1 cdot mathrm e mathrm j omega t nbsp Beschreibungsfunktion des nichtlinearen Systems Bearbeiten N e m a x u 1 j w t e 1 j w t a 1 j b 1 e j w t j e m a x e j w t displaystyle N e mathrm max frac tilde u 1 mathrm j omega t tilde e 1 mathrm j omega t frac a 1 mathrm j cdot b 1 cdot mathrm e mathrm j omega t mathrm j cdot e mathrm max cdot mathrm e mathrm j omega t nbsp Die Eingangsamplitude e m a x displaystyle e mathrm max nbsp der Sinusschwingung des nichtlinearen Systems wird haufig auch mit A displaystyle A nbsp bezeichnet N A b 1 j a 1 A b 1 A j a 1 A displaystyle N A frac b 1 mathrm j a 1 A frac b 1 A mathrm j frac a 1 A nbsp a 1 displaystyle a 1 nbsp und b 1 displaystyle b 1 nbsp sind die Koeffizienten der ersten Grundschwingungen des nichtlinearen Systems Fur die zahlreichen Formen der sogenannten Kennlinienregler lassen sich anhand der Geometrie der Kennlinien mit diesen Koeffizienten haufig mit a b c d displaystyle a b c d nbsp bezeichnet die Fourier Koeffizienten a 1 displaystyle a 1 nbsp und b 1 displaystyle b 1 nbsp der nichtlinearen Beschreibungsfunktionen N A displaystyle N A nbsp errechnen a 1 2 T T 2 T 2 u t cos w t d t mit w 2 p T displaystyle a 1 frac 2 T cdot int T 2 T 2 u t cdot cos omega t mathrm d t qquad text mit quad omega frac 2 pi T nbsp b 1 2 T T 2 T 2 u t sin w t d t mit w 2 p T displaystyle b 1 frac 2 T cdot int T 2 T 2 u t cdot sin omega t mathrm d t qquad text mit quad omega frac 2 pi T nbsp Die Beschreibungsfunktionen sind abhangig von der Funktion des nichtlinearen Kennlinien Elements Der Ausgang eines symmetrischen Zweipunkt Elements schwingt als Funktion der Eingangsschwingung synchron mit den Nulldurchgangen ohne Phasenverschiebung Unsymmetrische Mehrpunkt Elemente Amplitude d d displaystyle d neq left d right nbsp erzeugen einen Gleichanteil Bei Kennlinien Elementen mit Hysteresefunktionen sind die Beschreibungsfunktionen komplex Die Beschreibungsfunktionen kann man als einen Ersatzfrequenzgang eines nichtlinearen Systems bezeichnen Sie hangt von der Amplitude A displaystyle A nbsp der Eingangsschwingung e t displaystyle e t nbsp und der Frequenz der Dauerschwingung w 2 p f displaystyle omega 2 pi f nbsp ab und wird mit N A w displaystyle N A omega nbsp definiert Zum Unterschied des linearen Frequenzgangs mit der Abhangigkeit von w displaystyle omega nbsp ist die Beschreibungsfunktion im Wesentlichen von der Amplitude A displaystyle A nbsp als unabhangige Variable abhangig Bei rein statischer Nichtlinearitat entfallt die Abhangigkeit von der Frequenz Nachfolgend werden die Beschreibungsfunktionen als frequenzunabhangig fur statische Nichtlinearitaten betrachtet als N A displaystyle N A nbsp Die Aufstellung und umfangreichen Ableitungen der Beschreibungsfunktionen N A displaystyle N A nbsp der zahlreichen Kennlinien Elemente sind der Fachliteratur der Regelungstechnik zu entnehmen Die jeweiligen Beschreibungsfunktionen N A displaystyle N A nbsp der nichtlinearen statischen Systeme sind alle als algebraische Funktionen abgeleitet lediglich die relevanten geometrischen Eckpunkte mit den Buchstaben a b c d displaystyle a b c d nbsp unterscheiden sich konnen aber leicht anhand der beigefugten Kennlinien Skizzen zugeordnet werden Aufstellung bekannter Beschreibungsfunktionen von Kennlinien Elementen Bearbeiten Die nachfolgenden Beschreibungsfunktionen N A displaystyle N A nbsp gelten nur fur Anwendungen im Zustand der Harmonischen Balance 2 BenennungNichtlineares System BeschreibungsfunktionN A Re N A j Im N A displaystyle N A operatorname Re left N A right mathrm j operatorname Im left N A right nbsp Nichtlineare Kennlinie undinvertierte Ortskurve 1 N A displaystyle frac 1 N A nbsp Zweipunkt Regleru d fur e lt 0 0 fur e 0 d fur e gt 0 displaystyle u begin cases d amp text fur e lt 0 0 amp text fur e 0 d amp text fur e gt 0 end cases nbsp Re N A 4 d p A Im N A 0 displaystyle begin aligned operatorname Re left N A right amp frac 4 cdot d pi cdot A operatorname Im left N A right amp 0 end aligned nbsp nbsp Zweipunkt Regler mit Hystereseu d fur e c d d fur c lt e lt c d fur c e displaystyle u begin cases d amp text fur e leq c d d amp text fur c lt e lt c d amp text fur c leq e end cases nbsp Re N A 4 d p A 1 c A 2 Im N A 4 c d p A 2 mit A c displaystyle begin aligned operatorname Re left N A right amp frac 4 cdot d pi cdot A cdot sqrt 1 left frac c A right 2 operatorname Im left N A right amp frac 4 cdot c cdot d pi cdot A 2 text mit A geq c end aligned nbsp 3 nbsp Element mit Begrenzung Sattigung k tan a d c displaystyle k tan alpha d c nbsp u d fur e lt c k e fur c e c d fur c lt e displaystyle u begin cases d amp text fur e lt c k cdot e amp text fur c leq e leq c d amp text fur c lt e end cases nbsp Re N A 2 d p c arcsin c A 2 d p A 1 c A 2 Im N A 0 mit A c displaystyle begin aligned operatorname Re left N A right amp frac 2 cdot d pi cdot c cdot arcsin left frac c A right frac 2 cdot d pi cdot A cdot sqrt 1 left frac c A right 2 operatorname Im left N A right amp 0 text mit A geq c end aligned nbsp nbsp Element mit Totzonek tan a d c 0 c displaystyle k tan alpha d c 0 c nbsp u k e c fur e lt c 0 fur c e c k e c fur c lt e displaystyle u begin cases k cdot e c amp text fur e lt c 0 amp text fur c leq e leq c k cdot e c amp text fur c lt e end cases nbsp Re N A k 2 k p arcsin c e 2 k c p e 1 c A 2 Im N A 0 mit A c displaystyle begin aligned operatorname Re left N A right amp k frac 2 cdot k pi cdot arcsin left frac c e right frac 2 cdot k cdot c pi cdot e cdot sqrt 1 left frac c A right 2 operatorname Im left N A right amp 0 text mit A geq c end aligned nbsp nbsp Dreipunkt Regler Zweipunkt Regler mit Totzone u d fur e c 0 fur c lt e lt c d fur c e displaystyle u begin cases d amp text fur e leq c 0 amp text fur c lt e lt c d amp text fur c leq e end cases nbsp Re N A 4 d p A 1 c A 2 Im N A 0 mit A c displaystyle begin aligned operatorname Re left N A right amp frac 4 cdot d pi cdot A cdot sqrt 1 left frac c A right 2 operatorname Im left N A right amp 0 text mit A geq c end aligned nbsp nbsp Gleichung der Harmonischen Balance BearbeitenFur das Verfahren der Harmonischen Balance wird der geschlossene Regelkreis an der Stabilitatsgrenze fur die Fuhrungsgrosse w t 0 displaystyle w t 0 nbsp betrachtet Die Ausgangsgrosse des linearen Systems gilt Y j w G j w U j w displaystyle Y mathrm j omega G mathrm j omega cdot U mathrm j omega nbsp Die Nichtlinearitat verhalt sich im Schwingungsgleichgewicht wie ein lineares Ubertragungsglied U j w N A E j w N A Y j w displaystyle U mathrm j omega N A cdot E mathrm j omega N A cdot Y mathrm j omega nbsp Mit E j w Y j w displaystyle E mathrm j omega Y mathrm j omega nbsp ergibt sich aus den 3 Gleichungen die charakteristische Gleichung des nichtlinearen Regelkreises die mit Gleichung der Harmonischen Balance bezeichnet wird G j w N A 1 displaystyle G mathrm j omega cdot N A 1 nbsp Die haufigste Schreibweise der Gleichung der Harmonischen Balance lautet G j w 1 N A displaystyle G mathrm j omega frac 1 N A nbsp Wenn eine Dauerschwingung des nichtlinearen Regelkreises existiert dann stellt sich die Gleichung der Harmonischen Balance als die charakteristische Gleichung der nichtlinearen Regelung dar Die Losung der Gleichung der Harmonischen Balance kann formelmassig numerisch oder grafisch nach dem Zwei Ortskurven Verfahren erfolgen Formelmassige Losung der Gleichung der Harmonischen Balance Bearbeiten Bei komplizierten nichtlinearen Beschreibungsfunktionen oder linearen Systemen hoherer Ordnung kann die formelmassige Losung sehr schwierig werden In diesem Fall wahlt man die anschauliche grafische Losung des Zwei Ortskurven Verfahrens Die gesuchten unbekannten Grossen sind die Amplitude A displaystyle A nbsp der Nichtlinearitat und die Frequenz w displaystyle omega nbsp mit der das nichtlineare dynamische System schwingt 4 Die Gleichung der Harmonischen Balance wird umgeformt um zunachst A displaystyle A nbsp berechnen zu konnen N A 1 G j w displaystyle N A frac 1 G mathrm j omega nbsp Schreibt man die Gleichung in Realteil und Imaginarteil Re N A Re 1 G j w displaystyle operatorname Re N A operatorname Re frac 1 G mathrm j omega nbsp Im N A Im 1 G j w displaystyle operatorname Im N A operatorname Im frac 1 G mathrm j omega nbsp Damit ergeben sich zwei reelle Gleichungen Handelt es sich um eindeutige nichtlineare Kennlinien d h fur jede Eingangsgrosse existiert eine eindeutige Ausgangsgrosse wie bei den Zweipunkt und Dreipunkt Elementen so wird die Beschreibungsfunktion mit dem Imaginarteil zu Null Im 1 G j w 0 displaystyle operatorname Im frac 1 G mathrm j omega 0 nbsp Die Losung des Gleichungssystems vereinfacht sich auf die Losung zweier getrennter Gleichungen Zuerst werden die Pole von w displaystyle omega nbsp des linearen Systems aus Im 1 G j w 0 displaystyle operatorname Im 1 G j omega 0 nbsp bestimmt und in die rechte Seite der nachfolgenden Gleichung eingesetzt So kann A displaystyle A nbsp berechnet werden mit Re N A Re 1 G j w displaystyle operatorname Re N A operatorname Re frac 1 G mathrm j omega nbsp Zwei Ortskurven Verfahren der Harmonischen Balance Bearbeiten Die grafische Methode ist sehr anschaulich Man erhalt die Amplitude A displaystyle A nbsp und die Frequenz der Dauerschwingung w displaystyle omega nbsp des nichtlinearen Regelkreises wenn sich ein Schnittpunkt der beiden Ortskurven des Frequenzgangs G j w displaystyle G mathrm j omega nbsp und der negativen invertierten Beschreibungsfunktion 1 N A displaystyle 1 N A nbsp ergibt Es lasst sich feststellen ob es in einem Regelkreis zu Dauerschwingungen kommt und ob mehrere Dauerschwingungen moglich sind nbsp Darstellung des Zwei Ortskurven Verfahrens der Harmonischen BalanceDie in der Gleichung der Harmonischen Balance dargestellte negative invertierte Ortskurve der Beschreibungsfunktion N A displaystyle N A nbsp kann in das gleiche Diagramm der Ortskurve des Frequenzgangs G j w displaystyle G mathrm j omega nbsp als 1 N A displaystyle frac 1 N A nbsp eingetragen werden mit den Ortskurven Parametern von A displaystyle A nbsp Schwingt das System ergibt sich ein Schnittpunkt der beiden Ortskurven Daraus kann die Amplitude A displaystyle A nbsp und die Frequenz w 2 p f displaystyle omega 2 pi f nbsp abgelesen werden Je nach Art des nichtlinearen Systems ist die Lage der negativen inversen Ortskurve auf der realen Achse des Koordinatensystems eine Gerade und hat keine imaginaren Anteile Bei nichtlinearen Systemen mit einer Hysterese handelt es sich um eine im negativen imaginaren Bereich versetzte Gerade Bei nichtlinearen Systemen mit Dreipunktregler und Hysterese handelt es sich um ein Kennlinienfeld Je nach Art der Ortskurve des Frequenzgangs und der Ortskurve der Beschreibungsfunktion kann es zu mehreren Dauerschwingungen kommen weil es zu mehreren Schnittpunkten kommen kann Bei dem Dreipunktregler ist die Ortskurve der Beschreibungsfunktion auf der realen Achse gefaltet die beiden Aste der Ortskurve uberdecken sich d h fur eine gegebene Frequenz w displaystyle omega nbsp konnen im Schnittpunkt Schwingungen fur zwei verschieden grosse Amplituden A displaystyle A nbsp moglich sein Andererseits kann die Ortskurve des Frequenzgangs fur lineare Systeme grosserer Ordnung den Punkt Re 0 und Im 0 umschlingen Enthalt das lineare System eine Totzeit dann verlauft die Ortskurve des Frequenzgangs spiralformig um den Punkt Re 0 und Im 0 so dass es zu vielen Schnittpunkten und damit zu mehreren Dauerschwingungen unterschiedlicher Amplituden A displaystyle A nbsp und unterschiedlichen Frequenzen w displaystyle omega nbsp kommen kann Handelt es sich bei einem Regelkreis um beispielsweise einem Zweipunktregler und einem linearen System G 0 s displaystyle G 0 s nbsp erster oder zweiter Ordnung kann es laut dem Zwei Ortskurven Verfahren nicht zu Dauerschwingungen kommen weil die Ortskurve des Frequenzgangs nicht in den zweiten Quadranten kommen kann Leider weiss man durch numerische Berechnung oder aus Erfahrung mit Zweipunktreglern dass der Regelkreis mit einer Dauerschwingung antwortet Die Begrundung dieses Verhaltens wird damit erklart dass die Bedingung der harmonischen Schwingung der Ausgangsgrosse y t displaystyle y t nbsp des Regelkreises mit Tiefpassverhalten nicht erfullt wird Die Ausgangsgrosse y t displaystyle y t nbsp ist in diesem Fall keine harmonische Schwingung Dagegen konnen sich die Ortskurven eines Zweipunktreglers mit Hysterese durch die Lage der 1 N A Ortskurve im negativen imaginaren Bereich des 3 Quadranten sehr wohl mit der Ortskurve des linearen Systems G j w displaystyle G mathrm j omega nbsp 2 Ordnung treffen nbsp Sprungantwort des nichtlinearen Systems mit Dreipunktregler als DauerschwingungPraxisbeispiel der Harmonischen BalanceRegelstrecke mit globalem I Verhalten 5 G S s 2 s 2 s 1 0 5 s 1 displaystyle G S s frac 2 s cdot 2 cdot s 1 cdot 0 5 cdot s 1 nbsp Symmetrischer Dreipunktregler mit der Beschreibungsfunktion Umax 2 Totzone 0 5 Fuhrungsgrosse w t 1 Geometrische Eckpunkte Umax d Totzone c N A 4 d p A 1 c 2 A 2 displaystyle N A frac 4 cdot d pi cdot A cdot sqrt 1 c 2 A 2 nbsp Ergebnis Die Werte der Ortskurve der negativen invertierten Beschreibungsfunktion 1 N A displaystyle frac 1 N A nbsp werden als f A displaystyle f A nbsp aus der zugehorigen Gleichung N A displaystyle N A nbsp des Dreipunktreglers errechnet und auf der reellen Achse Re der Ortskurve des Frequenzgangs eingetragen Die zwei Aste der Ortskurve beginnen bei A c displaystyle A c nbsp enden bei A displaystyle A infty nbsp und kehren die Richtung um bei A 1 41 c displaystyle A 1 41 c nbsp Sie schneiden sich bei 1 N A 0 8 displaystyle frac 1 N A 0 8 nbsp mit A 1 97 displaystyle A 1 97 nbsp und A 0 517 displaystyle A 0 517 nbsp mit der Ortskurve des Frequenzgangs bei Realteil von 0 8 w 1 displaystyle omega 1 nbsp Die Dauerschwingung mit A 0 517 displaystyle A 0 517 nbsp ist instabil Die Amplitude A 0 517 displaystyle A 0 517 nbsp liegt im Grenzbereich der Totzone c 0 5 displaystyle c 0 5 nbsp in der noch keine Dauerschwingung moglich ist Die Beschreibungsfunktion mit der Amplitude A 1 97 displaystyle A 1 97 nbsp ergibt eine stabile Dauerschwingung Allgemein gilt die stabile Dauerschwingung stellt sich auf die grossere Amplitude ein Das mittels numerischer zeitdiskreter Methoden berechnete Diagramm mit dem Zeitverhalten der Systemgrossen stimmt bezuglich der Amplitude und der Frequenz nach dem Einschwingen nach ca 16 s vollig mit den Ergebnissen der Harmonischen Balance uberein Stabilitat der Grenzschwingungen Bearbeiten Die Ortskurve des Frequenzgangs eines linearen Systems als offener Regelkreis zeigt Dauerschwingungen an wenn nach dem Nyquist Verfahren in der komplexen Ebene die Ortskurve den Punkt R e 1 displaystyle mathit Re 1 nbsp und I m 0 displaystyle mathit Im 0 nbsp trifft Es handelt sich hierbei um keine stabilen Dauerschwingungen weil bei geringsten Anderungen der Pole von G j w displaystyle G mathrm j omega nbsp der kritische Punkt R e 1 D R e displaystyle mathit Re 1 pm Delta mathit Re nbsp nicht getroffen wird Die Folge ist bei einem Wert von R e 1 D R e displaystyle mathit Re 1 Delta mathit Re nbsp eine gedampfte Schwingung tritt auf oder bei einem Wert von R e 1 D R e displaystyle mathit Re 1 Delta mathit Re nbsp entsteht eine zunehmend aufklingende Schwingung Mit dem Zwei Ortskurven Verfahren lassen sich Dauerschwingungen eines nichtlinearen dynamischen zuruckgekoppelten Systems ermitteln wenn in der komplexen Ebene sich die Ortskurven der Beschreibungsfunktion des nichtlinearen Systems und die Ortskurve des Frequenzgangs des linearen Systems treffen Treffen sich die beiden Ortskurven nicht handelt es sich um einen schwingungsfreien Regelkreis Dabei ist zu berucksichtigen dass das lineare verzogernde System G j w displaystyle G mathrm j omega nbsp von mindestens zweiter Ordnung sein muss anderenfalls durchlauft die Ortskurve des Frequenzganges nicht den zweiten Quadranten und kann eine auf der realen Achse liegende Ortskurve N A displaystyle N A nbsp nicht treffen Bei einer Dauerschwingung treten Grenzbedingungen fur die Beschreibungsfunktion N A displaystyle N A nbsp auf Dabei werden die Grenzbedingungen der Schnittpunkte der Ortskurve der Beschreibungsfunktion fur drei Falle betrachtet bei der die Amplitude A displaystyle A nbsp durch transiente Storungen um einen kleinen Wert D A displaystyle Delta A nbsp erhoht oder vermindert wird und dann jeweils wieder in die Ausgangslage zuruckkehrt Es stellt sich die Frage was geschieht wenn eine Dauerschwingung besteht und durch eine kleine vorubergehende Amplitudenanderung positiv oder negativ auf den Gleichgewichtszustand der Dauerschwingung einwirkt Die Gleichung der harmonischen Balance wird fur diese Falle zu einer Ungleichung Folgende Regel gilt fur die meisten praktischen Anwendungen und insbesondere fur Schnittpunkte mit mehreren w Werten fur den kleinsten Wert von w 6 Der Schnittpunkt der beiden Ortskurven stellt eine stabile Grenzschwingung dar wenn mit wachsender Amplitude A displaystyle A nbsp der Betrag der Beschreibungsfunktion N A displaystyle left N A right nbsp abnimmt Eine instabile Grenzschwingung ergibt sich wenn der Betrag von N A displaystyle left N A right nbsp mit A displaystyle A nbsp der Beschreibungsfunktion zunimmt Beispiel fur einen Regelkreis mit einer Ortskurve eines Dreipunktreglers und einer PT3 Regelstrecke Die beiden Zweige der Ortskurve der reellen Beschreibungsfunktion N A displaystyle N A nbsp liegen negativ invertiert 1 N A displaystyle frac 1 N A nbsp ubereinander auf der reellen Achse nbsp Analyse der Grenzschwingungen bei transienten Abweichungen von der Amplitude A displaystyle A nbsp um D A displaystyle pm Delta A nbsp fur drei Falle getroffen beruhrt nicht getroffen Es werden folgende drei Falle unterschieden bei denen sich die beiden Ortskurven treffen beruhren oder nicht treffen 7 8 Fall 1 Die beiden Ortskurven treffen sichDie negativ invertierte Ortskurve der Beschreibungsfunktion N A displaystyle N A nbsp wird von der Ortskurve des Frequenzgangs G j w displaystyle G mathrm j omega nbsp getroffen Es ergeben sich theoretisch zwei Schnittpunkte und damit zwei Dauerschwingungen gleicher Frequenz mit unterschiedlicher Amplitude Dies gilt fur den Fall dass auch der kurze Arm der Beschreibungsfunktion A c displaystyle A c nbsp bis A 2 c displaystyle A sqrt 2 c nbsp Umkehrpunkt getroffen wird Nachfolgend wird der Zweig der Beschreibungsfunktion vom Umkehrpunkt A 2 c displaystyle A sqrt 2 c nbsp im Abstand zur Ordinate mit 1 N A p c 2 d displaystyle frac 1 N A frac p c 2 d nbsp mit zunehmender Grosse der Amplitude bis A displaystyle A infty nbsp betrachtet Wird die Amplitude infolge einer transienten Storung um einen kleinen Betrag D A displaystyle Delta A nbsp vergrossert liegt der neue Wert A displaystyle A nbsp links von der Ortskurve des Frequenzgangs G j w displaystyle G mathrm j omega nbsp und liegt ausserhalb der Ortskurve nicht umschlungen Nach dem Nyquist Kriterium ist der geschlossene Regelkreis stabil Die Schwingung klingt ab Nach Ende der Storung wird wieder der Zustand der stabilen Dauerschwingung erreicht Wird die Amplitude infolge einer transienten Storung um einen kleinen Betrag D A displaystyle Delta A nbsp verkleinert liegt der neue Wert A displaystyle A nbsp rechts von der Ortskurve des Frequenzgangs G j w displaystyle G mathrm j omega nbsp Der neue Wert A displaystyle A nbsp wird von der Ortskurve G j w displaystyle G mathrm j omega nbsp umschlungen und liegt innerhalb des Bereichs der der Ortskurve Nach dem Nyquist Kriterium ist der geschlossene Regelkreis instabil Die Schwingung klingt auf die Amplitude vergrossert sich auf ihren ursprunglichen Wert Nach Verschwinden der Storungen ist die Grenzschwingung stabil Fall 2 Die beiden Ortskurven beruhren sich am Umkehrpunkt A 2 c displaystyle A sqrt 2 c nbsp der beiden AsteSind die Schnittpunkte nur Beruhrungspunkte dann wird die Grenzschwingung als semistabil bezeichnet Wird die Amplitude infolge einer transienten Storung um einen kleinen Betrag vergrossert liegt der neue Wert A displaystyle A nbsp links von der Ortskurve des Frequenzgangs und ausserhalb der Ortskurve nicht umschlungen Nach dem Nyquist Kriterium ist der geschlossene Regelkreis stabil Die Schwingung klingt ab Nach Ende der Storung wird wieder der Zustand der stabilen Dauerschwingung erreicht Wird die Amplitude infolge einer transienten Storung um einen kleinen Betrag D A displaystyle Delta A nbsp verkleinert liegt der neue Wert A displaystyle A nbsp rechts von der Ortskurve des Frequenzgangs G j w displaystyle G mathrm j omega nbsp und wird von der Ortskurve G s displaystyle G s nbsp umschlungen Die semistabile Grenzschwingung klingt auf Das Regelsystem ist instabil Fall 3 Die beiden Ortskurven treffen sich nicht und beruhren sich nichtDer geschlossene Regelkreis ist stabil Es treten keine Dauerschwingungen auf Vermeidung von Dauerschwingungen BearbeitenNichtlineare dynamische Systeme konnen z B bei einem Einsatz eines einfachen Zweipunktreglers gewollt stabile Dauerschwingungen hervorrufen Die Lage der Ortskurve der Beschreibungsfunktion 1 N A displaystyle frac 1 N A nbsp liegt auf der realen Achse der Ortskurve des Frequenzgangs am Punkt R e 0 displaystyle mathit Re 0 nbsp und startet bei A 0 displaystyle A 0 nbsp und endet bei A displaystyle A infty nbsp Deshalb besteht keine Moglichkeit mit dem Zweipunktregler Dauerschwingungen zu vermeiden Anmerkung Ein Zweipunktregler schwingt bereits in Verbindung mit einem linearen Verzogerungsglied erster Ordnung obwohl die beiden Ortskurven sich nicht treffen konnen Die Harmonische Balance ist fur diesen Fall nicht gultig weil die Voraussetzung der harmonischen Schwingungen sagezahnahnliche Schwingung an y t displaystyle y t nbsp nicht gegeben sind Die gilt auch noch fur ein Verzogerungsglied zweiter Ordnung PT2 Glied Mit dem Einsatz anderer statischer nichtlinearer Systeme wie z B durch den Dreipunktregler lassen sich schwingungsfreie Regelkreise mit asymptotischem Verhalten realisieren Voraussetzung ist eine Regelstrecke mit globalem I Verhalten Es stellt sich die Frage welche Massnahmen mit Hilfe des Zwei Ortskurven Verfahrens realisierbar sind 9 Folgende Massnahmen sind gegeben damit die Gleichung der Harmonischen Balance zu einer Ungleichung wird bei der sich die beiden Ortskurven nicht treffen konnen Die Ungleichung der Harmonischen Balance konnte lauten 1 G j w gt N A oder 1 gt N A G j w displaystyle left frac 1 G mathrm j omega right gt left N A right quad text oder quad 1 gt left N A cdot G mathrm j omega right nbsp Reduzierung der Verstarkung des linearen Systems G j w displaystyle G mathrm j omega nbsp Als einfachste Methode konnte die Ortskurve des Frequenzgangs durch Andern des linearen Systems so gestaltet werden z B Reduzierung der Verstarkung dass sie die Ortskurve der Beschreibungsfunktion nicht schneidet Damit wurden die dynamischen Eigenschaften des Regelkreises verschlechtert werden Vergrosserung des Abstandes der Ortskurve der Beschreibungsfunktion von der Ordinate bei R e 0 displaystyle mathit Re 0 nbsp Die Ortskurve der Beschreibungsfunktion konnte geandert werden in dem der Abstand des Umkehrpunktes zur Ordinate der zwei Arme der Beschreibungsfunktion p c 2 d displaystyle frac pi cdot c 2 cdot d nbsp vergrossert wird Moglich ist die Verkleinerung der Amplitude d displaystyle d nbsp oder die Vergrosserung der Totzone c displaystyle c nbsp Aus regelungstechnischen Grunden in Bezug auf die Genauigkeit der Regelgrosse y t displaystyle y t nbsp muss die Totzone moglichst klein sein denn auf Regelabweichungen innerhalb der Totzone reagiert der Regler nicht Die Verkleinerung der Amplitude d displaystyle d nbsp Regler Ausgangsgrosse U m a x displaystyle pm U mathrm max nbsp fordert eine Verkleinerung der Fuhrungsgrosse w t displaystyle w t nbsp Deshalb fuhrt diese Betrachtungsweise zu keinem nennenswerten Ergebnis Teil Kompensation der Verzogerungsglieder des linearen Systems G j w displaystyle G mathrm j omega nbsp Anderung Teil Kompensation der linearen Regelstrecke durch Einsatz eines realen PD1 Gliedes Diese Methode der Verbesserung der Systemgeschwindigkeit durch Kompensation trager PT1 Glieder durch reale PD1 Glieder ist auch in der linearen Regelungstechnik ublich Sie fuhrt auch bei nichtlinearen Regelkreisen zu einem akzeptablen Ergebnis Kompensation trager PT1 Glieder durch PD1 Glieder Bearbeiten Durch Austausch trager Verzogerungsglieder des linearen dynamischen Systems erhalt das System ein schnelleres Zeitverhalten Beispiel lineare Regelstrecke G S s K s T 1 s 1 T 2 s 1 displaystyle G S s frac K s cdot T 1 cdot s 1 cdot T 2 cdot s 1 nbsp dd Lineares reales PD1 Glied Kompensierglied G KOMP s T 1 s 1 T N s 1 displaystyle G text KOMP s frac T 1 cdot s 1 T N cdot s 1 nbsp dd Die Zeitkonstante T 1 displaystyle T 1 nbsp der Strecke soll grosser als T 2 displaystyle T 2 nbsp der Strecke sein Das PD Glied mit Zeitkonstante T 1 displaystyle T 1 nbsp des Reglers soll das Verzogerungsglied mit der Zeitkonstante T 1 displaystyle T 1 nbsp kompensieren Die Zeitkonstante T N displaystyle T N nbsp muss kleiner T 1 displaystyle T 1 nbsp sein Damit ergibt sich die neue Ubertragungsfunktion des linearen Regelstrecke in Frequenzgang Darstellung zu G S 1 j w K j w T N j w 1 T 2 j w 1 displaystyle G S1 mathrm j omega frac K mathrm j omega cdot T N cdot mathrm j omega 1 cdot T 2 cdot mathrm j omega 1 nbsp dd nbsp Sprungantwort des nichtlinearen Regelkreises mit Dreipunktregler ohne Dauerschwingungen Darstellung des EinschwingvorgangsMit dieser Massnahme schmiegt sich der Verlauf des Frequenzgangs G j w displaystyle G mathrm j omega nbsp naher an die Ordinate beim Durchtreten vom 3 in den 2 Quadranten an Sind die Parameter der Beschreibungsfunktion so gewahlt dass die beiden Ortskurven sich nicht beruhren ist der Regelkreis asymptotisch stabil Die Minimierung der Totzone und damit die Optimierung des schwingungsfreien Regelkreises mit einem Dreipunktregler lasst sich nur mit der numerischen Simulation bestimmen Beispiel eines schwingungsfreien Regelkreises mit einem Dreipunktregler Bearbeiten Fur einen schwingungsfreien Regelkreis mit einem Dreipunktregler wird vorausgesetzt dass das lineare dynamische System die Regelstrecke ein globales I Verhalten zeigt Typisch fur solche Regelstrecken sind Stellantriebe die durch den Dreipunktregler auf die Bewegungen Vorlauf Rucklauf oder Ruhezustand gebracht werden Der Stillstand entspricht der Soll Positionierung Werden die Parameter des Berechnungsbeispiels des mit Dauerschwingungen behafteten Regelkreises mit dem Dreipunktregler wie folgt modifiziert entsteht wie in der Grafik dargestellt ein schwingungsfreier Regelkreis Die beiden Verzogerungen der PT1 Glieder werden durch PD1 Glieder um den Faktor 10 reduziert G 0 j w 2 j w 0 1 j w 1 0 05 j w 1 displaystyle G 0 mathrm j omega frac 2 mathrm j omega cdot 0 1 cdot mathrm j omega 1 cdot 0 05 cdot mathrm j omega 1 nbsp dd Die bezogene Grosse 1 100 der Rechteck Schwingamplitude U max displaystyle pm U text max nbsp wird auf 1 2 displaystyle pm 1 2 nbsp gesetzt Die bezogene Grosse der Totzone wird auf c 0 1 displaystyle pm c pm 0 1 nbsp gesetzt die maximale bezogene Fuhrungsgrosse w t 1 displaystyle w t 1 nbsp gesetzt Die Grosse der Totzone bestimmt bei welcher Regelabweichung der Regler ansprechen soll In diesem Fall richtet sich das Zeitverhalten Integrationskonstante T N displaystyle T N nbsp der auszulegenden Regelstrecke nach der Totzone Nur ein langsamer Stellantrieb erlaubt fur einen schwingungsfreien Betrieb eine kleine Totzone und damit eine genaue Positionierung Je kleiner die Totzone umso haufiger und kurzer sind die Rechteckschwingungen u t displaystyle u t nbsp wahrend des Einschwingvorgangs der Regelgrosse y t displaystyle y t nbsp Fazit zur Anwendung der Harmonischen Balance BearbeitenDie Harmonische Balance zeigt interessante Aspekte des Zustandekommens von Dauerschwingungen in nichtlinearen Regelkreisen und auch Moglichkeiten der Verhinderung von Dauerschwingungen Es gibt kein Berechnungsverfahren der Analyse des Systemverhaltens von nichtlinearen Regelkreisen mit Totzeit und Hysterese Elementen das annahernd so leistungsfahig und gleichzeitig so einfach ist wie die numerische rekursive Berechnung eines nichtlinearen Systems nach dem Euler Streckenzug Verfahren mit der diskreten Zeit D t displaystyle Delta t nbsp und der Berechnungsfolge k 0 1 2 3 k max displaystyle k 0 1 2 3 ldots k text max nbsp Ob die statische Nichtlinearitat mit logischen Programmbefehlen und die Linearfaktoren der linearen Teilsysteme mit Differenzengleichungen berechnet werden oder die kommerziellen numerischen Rechenprogramme wie Matlab oder Simulink zur Anwendung kommen die Erkenntnisse aus der Harmonischen Balance konnen im Zeitbereich unmittelbar simuliert und bestatigt werden Literatur BearbeitenHolger Lutz Wolfgang Wendt Taschenbuch der Regelungstechnik mit MATLAB und Simulink 12 Auflage Europa Lehrmittel Haan Gruiten 2021 ISBN 978 3 8085 5870 6 Gerd Schulz Regelungstechnik 1 Lineare und Nichtlineare Regelung Rechnergestutzter Reglerentwurf 3 Auflage Oldenbourg 2007 ISBN 3 486 58317 4 Otto Follinger Nichtlineare Regelungen II Mit Harmonische Balance Popow und Kreiskriterium Hyperstabilitat Synthese im Zustandsraum 8 Auflage Oldenbourg Munchen Wien 1993 ISBN 3 486 22503 0 Heinz Unbehauen Regelungstechnik II Mit Zustandsregelungen digitale und nichtlineare Regelsysteme 9 Auflage Vieweg Teubner Wiesbaden 2007 ISBN 978 3 528 83348 0 Einzelnachweise Bearbeiten Holger Lutz Wolfgang Wendt Taschenbuch der Regelungstechnik Kapitel Harmonische Linearisierung mit der Beschreibungsfunktion Holger Lutz Wolfgang Wendt Taschenbuch der Regelungstechnik Kapitel Nichtlineare Regelungen Tabelle Nichtlineare statische Elemente Jurgen Adamy Nichtlineare Systeme und Regelungen 2018 S 68 doi 10 1007 978 3 662 55685 6 springer com abgerufen am 21 Juni 2020 Otto Follinger Nichtlineare Regelungen II Harmonische Balance Kapitel Losung der Gleichung der Harmonischen Balance Gerd Schulz Regelungstechnik 1 Kapitel Nichtlineare Regler Praxisbeispiel ist identisch mit Fachbuch Heinz Unbehauen Regelungstechnik II Kapitel Nichtlineare Regelsysteme Unterkapitel Stabilitatsuntersuchung mittels der Beschreibungsfunktion Holger Lutz Wolfgang Wendt Taschenbuch der Regelungstechnik Kapitel Nichtlineare Regelungen Unterkapitel Stabilitat von Grenzschwingungen Otto Follinger Nichtlineare Regelungen II Kapitel Nichtlineare Regelungen Unterkapitel Stabilitatsverhalten von Dauerschwingungen Otto Follinger Nichtlineare Regelungen II Harmonische Balance Kapitel Stabilisierung nichtlinearer Regelungen Abgerufen von https de wikipedia org w index php title Methode der harmonischen Balance amp oldid 229344161