Dynamische Systeme höherer Ordnung: Unterschied zwischen den Versionen

Aus SystemPhysik
Inhalt hinzugefügt Inhalt gelöscht
 
(34 dazwischenliegende Versionen von 3 Benutzern werden nicht angezeigt)
Zeile 1: Zeile 1:
==Lernziele==
==Lernziele==
In dieser Vorlesung lernen Sie
In dieser Vorlesung lernen Sie
*wie man die Kreisfrequenz eines harmonischen Oszillators berechnet.
*wie man die Kreisfrequenz, Frequenz und Schwingungsdauer eines harmonischen Oszillators berechnet.
*wie die Potential-Zeit-, die Stromstärke-Zeit- und die Leistungs-Zeit-Funktion bei einem harmonischen Oszillator formuliert werden.
*wie die Potential-Zeit-, die Stromstärke-Zeit- und die Leistungs-Zeit-Funktion bei einem harmonischen Oszillator formuliert werden.
*wie man die Abklingzeit und die Kreisfrequenz bei einem gedämpften Oszillator berechnet.
*wie man die Abklingzeit und die Kreisfrequenz bei einem gedämpften Oszillator berechnet.
*unter welchen Bedingungen ein Serie-Schwingkreis bei harmonischer Anregung den stärksten Strom durchlässt
*unter welchen Bedingungen ein Serie-Schwingkreis bei harmonischer Anregung den stärksten Strom durchlässt
*ein paar Darstellungsformen für das Verhalten dynmischer Systeme kennen.
*wie das Verhalten dynamischer Systeme modelliert und analysiert wird.


==Problemstellung==
==Problemstellung==
Zeile 12: Zeile 12:
Das Verhalten des U-Rohrs kann durch zwei Speicher mit konstanter [[Kapazität]] und einem Leiter mit [[Widerstand]] und [[Induktivität]] modelliert werden. Analoge Systeme kennt man aus der Elektrodynamik (zwei über eine Spule verbunden Kondensatoren), aus der Translationsmechanik (zwei über eine Feder verbundene Luftkissenfahrzeuge) und aus der Rotationsmechanik (zwei über eine Drehfeder verbundene Schwungräder). Dominiert das Widerstandselement (Strömungswiderstand, elektrischer Widerstand der Spule, Dämpfer in Serie zur Feder), klingt der Potentialunterschied (Druck, Spannung, Geschwindigkeit, Winkelgeschwindigkeit) exponentiell ab, falls sich alle Elemente linear verhalten.
Das Verhalten des U-Rohrs kann durch zwei Speicher mit konstanter [[Kapazität]] und einem Leiter mit [[Widerstand]] und [[Induktivität]] modelliert werden. Analoge Systeme kennt man aus der Elektrodynamik (zwei über eine Spule verbunden Kondensatoren), aus der Translationsmechanik (zwei über eine Feder verbundene Luftkissenfahrzeuge) und aus der Rotationsmechanik (zwei über eine Drehfeder verbundene Schwungräder). Dominiert das Widerstandselement (Strömungswiderstand, elektrischer Widerstand der Spule, Dämpfer in Serie zur Feder), klingt der Potentialunterschied (Druck, Spannung, Geschwindigkeit, Winkelgeschwindigkeit) exponentiell ab, falls sich alle Elemente linear verhalten.


Dynamischen Systemen, die aus zwei Speichern bestehen und die über einen Leiter miteinander verbunden sind, können mathematisch zu einem System mit nur einem Speicher mit reduzierter Kapazität zusammenfasst werden. Systeme mit einem oder zwei Speicher und einer Induktivität, die also zwei unabhängig Energiespeicher besitzen, nennt man Systeme 2. Ordnung, weil man zur Berechnung des Verhaltens zweimal integrieren muss.
Dynamischen Systemen, die aus zwei Speichern bestehen und die über einen Leiter miteinander verbunden sind, können mathematisch zu einem System mit nur einem Speicher und reduzierter Kapazität zusammenfasst werden. Systeme mit Mengenspeicher und Induktivität nennt man Systeme 2. Ordnung, weil man zur Berechnung des Verhaltens zweimal integrieren muss. Systeme 2. Ordnung besitzen zwei unabhängige Energiespeicher.


==Speicher und Leiter==
==Speicher und Leiter==
Zeile 50: Zeile 50:
|}
|}


Die mechanischen Systemparameter sind oft reziprok definiert. So ist der Impulswiderstand gleich dem Reziprokwert der Dämpferkonstanten und die mechanische Induktivität gleich dem Reziprokwert der Federkonstanten.
Die mechanischen Systemparameter sind oft reziprok definiert: statt eines Impulswiderstandes führt man eine Dämpferkonstante ein; anstelle einer Impulsinduktivität definiert man eine Federkonstante


:<math>F_x=k\Delta v_x</math> statt <math>\Delta v_x= R_{px}F_x</math>
==harmonischer Oszillator==
:<math>F_x=D\Delta x</math> statt <math>\Delta v_x=L_{px}\dot F_x</math>
Ein harmonischer [[Oszillator]] ist ein schwingungsfähiges System, das einem oder zwei Speicher mit konstanter Kapazität und einem Leiter mit stromunabhängiger Induktivität besteht. Lenkt man ein solches System aus und lässt es los, es sinusförmig (=harmonisch) um seine Ruhelage, wobei die Schwingungsdauer unabhängig von der Grösse der Auslenkung ist. Beispiele für harmonische Oszillatoren sind Federpendel, elektrische Schwingkreise und U-Rohr mit sehr grossem Querschnitt.


==Harmonischer Oszillator==
Der harmonische Oszillator ist ein wichtiges Modellsystem der Physik. Er ist durch nur zwei Parameter vollständig beschrieben, die Gesamtkapazität und die Induktivität. Viele komplexere Systeme verhalten sich bei kleinen Auslenkungen näherungsweise wie harmonische Oszillatoren. Der harmonische Oszillator in der Quantenmechanik ist eines der wenigen quantenmechanischen Systeme, das sich ohne Näherungen berechnen lässt.
Ein harmonischer [[Oszillator]] ist ein schwingungsfähiges System, das aus einem oder zwei Speicher mit konstanter Kapazität und einem Leiter mit stromunabhängiger Induktivität besteht. Lenkt man ein solches System aus und lässt es los, schwingt es sinusförmig (harmonisch) um seine Ruhelage, wobei die Schwingungsdauer unabhängig von der Grösse der Auslenkung ist. Beispiele für harmonische Oszillatoren sind Federpendel, elektrische Schwingkreise und U-Rohr mit sehr grossem Querschnitt.


Der harmonische Oszillator ist ein wichtiges Modellsystem der Physik. Er ist durch nur zwei Parameter (Gesamtkapazität und die Induktivität) vollständig beschrieben. Viele komplexere Systeme verhalten sich bei kleinen Auslenkungen näherungsweise wie harmonische Oszillatoren. Der harmonische Oszillator der Quantenmechanik ist eines der wenigen quantenmechanischen Systeme, das sich ohne Näherungen berechnen lässt.
In der [[Elektrodynamik]] bildet ein über eine ideale Spule ([[Induktivität]] ''L'') kurz geschlossener [[Kondensator]] ([[Kapazität]] ''C'') einen harmonischen Oszillator. Lädt man den Kondensator auf die Spannung ''U<sub>0</sub>'' auf und verbindet ihn dann mit einer idealen Spule, ist die Spannung über beiden Elementen unabhängig vom Zustand gleich gross, oder im Sinne des Maschensatzes ist die Umlaufspannung gleich null

In der [[Elektrodynamik]] bildet ein über eine ideale Spule ([[Induktivität]] ''L'') kurz geschlossener [[Kondensator]] ([[Kapazität]] ''C'') einen harmonischen Oszillator. Lädt man den Kondensator auf die Spannung ''U<sub>0</sub>'' auf und verbindet ihn dann mit der idealen Spule, ist die Spannung über beiden Elementen zu jedem Zeitpunkt gleich gross (die Umlaufspannung muss immer gleich null sein)


:<math>U_C+U_L=0</math>
:<math>U_C+U_L=0</math>
Zeile 93: Zeile 96:
:<math>P=UI=U_0I_0cos(\omega_0t+\delta)sin(\omega_0t+\delta)</math>
:<math>P=UI=U_0I_0cos(\omega_0t+\delta)sin(\omega_0t+\delta)</math>


Setzt man die Phasenverschiebung ''&delta;'' gleich null (freie Wahl des Zeitnullpunktes), erhält man eine Leistung, die mit doppeter [[Frequenz]] (Kreisfrequenz durch 2''&pi;'') schwingt
Setzt man die Phasenverschiebung ''&delta;'' gleich null (freie Wahl des Zeitnullpunktes), erhält man eine Leistung, die mit doppelter [[Frequenz]] (Kreisfrequenz durch 2''&pi;'') schwingt


:<math>P=UI=\frac{U_0I_0}{2}sin(2\omega_0t)</math>
:<math>P=UI=\frac{U_0I_0}{2}sin(2\omega_0t)</math>


Während einer [[Periode]] oder Schwingungsdauer (2''&pi;'' durch Kreisfrequenz) ist der Kondensator zweimal geladen und der Strom erreicht zweimal seine maximale Stärke. Weil beide Bauteile zusammen mit der Ladung (Kondensator) oder mit dem Strom (Spule) Energie speichern, muss die Energie mit der doppelten Frequenz des Oszillators umgesetzt werden.
Während einer [[Periode]] oder Schwingungsdauer (2''&pi;'' durch Kreisfrequenz) ist der Kondensator zweimal geladen und der Strom erreicht zweimal seine maximale Stärke. Weil beide Bauteile zusammen mit der Ladung (Kondensator) oder mit dem Strom (Spule) Energie speichern, wird die Energie mit doppelter Oszillatorfrequenz hin und her verschoben.


'''Video:''' [https://www.youtube.com/watch?v=VkXZL-VAKxs harmonischer Oszillator]
'''Video:''' [https://www.youtube.com/watch?v=VkXZL-VAKxs harmonischer Oszillator]
Zeile 106: Zeile 109:
:<math>U_C+U_R+U_L=0</math>
:<math>U_C+U_R+U_L=0</math>


Setzt man wie beim harmonischen Oszillator die konstitutiven Gesetze für die [[Kapazität]], den [[Widerstand]] und die [[Induktivität]] ein und erstetzt die Stromstärke über die [[Bilanz]]gleichung durch die Änderungsrate der Kondensatorladung, erhält man eine etwas allgemeinere lineare Differentialgleichung 2. Ordnung mit konstanten Koeffizienten
Setzt man wie beim harmonischen Oszillator die konstitutiven Gesetze für die [[Kapazität]], den [[Widerstand]] und die [[Induktivität]] ein und ersetzt die Stromstärke über die [[Bilanz]]gleichung durch die Änderungsrate der Kondensatorladung, erhält man eine '''lineare Differentialgleichung 2. Ordnung mit konstanten Koeffizienten'''


:<math>\frac{Q}{C}+R\dot Q+L\ddot Q=0</math> oder <math>\frac{1}{LC}Q+\frac{R}{L}\dot Q+\ddot Q=0</math>
:<math>\frac{Q}{C}+R\dot Q+L\ddot Q=0</math> oder <math>\frac{1}{LC}Q+\frac{R}{L}\dot Q+\ddot Q=0</math>
Zeile 118: Zeile 121:
:<math>\left(\left(\frac{1}{LC}-\frac{R}{L\tau}+\frac{1}{\tau^2}-\omega^2\right)cos(\omega t+\delta)+\left(-\frac{R\omega}{L}+\frac{2\omega}{\tau}\right)sin(\omega t+\delta)\right)Q_0e^{-t/\tau}=0</math>
:<math>\left(\left(\frac{1}{LC}-\frac{R}{L\tau}+\frac{1}{\tau^2}-\omega^2\right)cos(\omega t+\delta)+\left(-\frac{R\omega}{L}+\frac{2\omega}{\tau}\right)sin(\omega t+\delta)\right)Q_0e^{-t/\tau}=0</math>


Die linke Seite dieser Gleichung kann nur zu allen Zeitpunkten gleich null sein, wenn beide Klammerausdrücke vor den Winkelfunktionen den Wert null annehmen. Damit gilt für die Abklingzeit
Die linke Seite dieser Gleichung kann nur zu allen Zeitpunkten gleich null sein, wenn beide Klammerausdrücke den Wert null annehmen. Damit gilt für die Abklingzeit


:<math>\tau =2\frac{L}{R}</math>
:<math>\tau =2\frac{L}{R}</math>
Zeile 126: Zeile 129:
:<math>\omega=\sqrt{\frac{1}{LC}-\frac{R^2}{4L^2}}=\sqrt{\omega_0^2-\frac{1}{\tau^2}}</math>
:<math>\omega=\sqrt{\frac{1}{LC}-\frac{R^2}{4L^2}}=\sqrt{\omega_0^2-\frac{1}{\tau^2}}</math>


Ohne Dämpfung ist die Abklingzeit unendlich gross und die Kreisfrequenz entspricht dem Wert des harmonischen Oszillators. Unterschreitet die Abklingzeit einen gewissen Wert, wird der Ausdruck unter der Wurzel negativ und das System schwingt gar nicht mehr, sondern kriecht - wie wenn die Induktivität nicht da wäre - gegen die Gleichgewichtslage. Die nachfolgende Tabelle zeigt noch die Analogie zwischen elektrischen und mechanischen Schwingkreisen
Ohne Dämpfung ist die Abklingzeit unendlich gross und die Kreisfrequenz entspricht dem Wert des harmonischen Oszillators. Unterschreitet die Abklingzeit einen gewissen Wert, wird der Ausdruck unter der Wurzel negativ und das System schwingt gar nicht mehr, sondern kriecht - wie wenn die Induktivität nicht da wäre - gegen die Gleichgewichtslage. Man unterscheidet drei Dämpfungsarten
*<math>\omega=\sqrt{\omega_0^2-\frac{1}{\tau^2}}>0</math>: '''unterkritische''' Dämpfung: System schwingt bezüglich der Gleichgewichtslage hin und her
*<math>\omega=\sqrt{\omega_0^2-\frac{1}{\tau^2}}=0</math>: '''kritische''' Dämpfung: System bewegt sich gegen die Gleichgewichtslage ohne zu überschwingen
*<math>\omega=\sqrt{\omega_0^2-\frac{1}{\tau^2}}<0</math>: '''überkritische''' Dämpfung: System kriecht gegen Gleichgewichtslage

Die nachfolgende Tabelle zeigt noch die Analogie zwischen elektrischen und mechanischen Schwingkreisen


{|class="wikitable"
{|class="wikitable"
Zeile 164: Zeile 172:
|-
|-
|Abklingzeit
|Abklingzeit
|<math>\tau=\frac{L}{R}</math>
|<math>\tau=\frac{2L}{R}</math>
|<math>\tau=\frac{L_V}{R_V}</math>
|<math>\tau=\frac{2L_V}{R_V}</math>
|<math>\tau=\frac{k}{D}</math>
|<math>\tau=\frac{2k}{D}</math>
|<math>\tau=\frac{k^*}{D^*}</math>
|<math>\tau=\frac{2k^*}{D^*}</math>
|-
|-
|}
|}
''k'' ist die Dämpferkonstante, ''D'' die Federkonstante oder Richtgrösse. Bei der entsprechenden Drehgrösse ist noch ein Stern zugefügt. In der Mechanik sind Feder und Dämpfer oft parallel statt in Reihe geschaltet. Der Strom verzweigt sich und man spricht man von einem Parallel-Schwingkreis. Das U-Ror verhält sich nur bedingt analog zu den andern Systemen, weil die Strömung eher turbulent als laminar ist und das System sich damit nichtlinear verhält.
''k'' ist die Dämpferkonstante, ''D'' die Federkonstante oder Richtgrösse. Bei der entsprechenden Drehgrösse ist noch ein Stern zugefügt. In der Mechanik sind Feder und Dämpfer oft parallel statt in Reihe geschaltet. Der Strom verzweigt sich und man spricht man von einem Parallel-Schwingkreis. Das U-Rohr verhält sich nur bedingt analog zu den andern Systemen, weil die Strömung eher turbulent als laminar ist und das System sich damit nichtlinear verhält.

'''Video:''' [https://www.youtube.com/watch?v=q1MuIE-TQDk Schwingkreis: Analogie und Dualität]


==Resonanz==
==Resonanz==
Ein elektrischer Schwingkreis kann mit einer Wechselspannungsquelle, die man in den Kreis einfügt, angeregt werden. Alle vier Elemente (Spannungsquelle, Kapazität, Widerstand und Induktivität) sind miteinander zu einem unverzweigten Kreis verbunden. Die beschreibende Gleichung lautet
Ein elektrischer Schwingkreis kann mit einer Wechselspannungsquelle, die man in den Kreis einfügt, angeregt werden. Alle vier Elemente (Spannungsquelle, Kapazität, Widerstand und Induktivität) sind miteinander zu einem nicht verzweigten Kreis verbunden. Die beschreibende Gleichung lautet


:<math>U_C+U_R+U_L=\hat Ucos(\hat{\omega}t)</math>
:<math>U_C+U_R+U_L=\hat Ucos(\hat{\omega}t)</math>
Zeile 191: Zeile 201:
:<math>U_L=L\dot I</math> also <math>U_L=-\hat\omega Lsin(\hat{\omega}t+\delta))</math>
:<math>U_L=L\dot I</math> also <math>U_L=-\hat\omega Lsin(\hat{\omega}t+\delta))</math>


Würde mann jedes der drei Elemente einzeln mit der Spannungsquelle verbinden, wäre die Lösung der Gleichung einfach zu finden: beim Widerstand schwingt der Strom in Phase mit der angelegten Spannung, bei der Kapazität läuft der Strom der Spannung eine Viertelperiode voraus (der Strom baut die Spannung auf) und bei der Induktivität eine Vierelperiode nach (Spannung ändert Stromstärke).
Würde man jedes der drei Elemente einzeln mit der Spannungsquelle verbinden, wäre die Lösung der Gleichung einfach zu finden: beim Widerstand schwingt der Strom in Phase mit der angelegten Spannung, bei der Kapazität läuft der Strom der Spannung eine Viertelperiode voraus (der Strom baut die Spannung auf) und bei der Induktivität eine Viertelperiode nach (Spannung ändert Stromstärke).


Strom- und Spannungsamplituden unterscheiden sich bei allen drei Elementen um einen Faktor, der entweder konstant ist, sich proportional oder rezirok-proportional mit der angelegten Kreisfrequenz ändert. Dieser Faktor, eigentlich eine Verallgemeinerung des Widerstandes, heisst [[Impedanz]] ''Z''. Die Verallgemeinerung des Leitwerts, des Kehrwerts der Impedanz, nennt man [[Admittanz]]. Betrag und Phase dieser Grössen lassen sich entweder als Zeiger oder als komplexe Zahl darstellen, was das Berechnen von Strom und Spannung stark erleichter. Der Betrag der Impedanz, also das Verhältnis zwischen den Amplituden der Spannung und der Stromstärke, Scheinwiderstand genannt, kann bei Serie- oder Reihenschaltung mit Hilfe des Pythagoras gerechnet werden
Strom- und Spannungsamplituden unterscheiden sich bei allen drei Elementen um einen Faktor, der entweder konstant ist, sich proportional oder reziprok-proportional mit der angelegten Kreisfrequenz ändert. Dieser Faktor, eigentlich eine Verallgemeinerung des Widerstandes, heisst [[Impedanz]] ''Z''. Die Verallgemeinerung des Leitwerts, des Kehrwerts der Impedanz, nennt man [[Admittanz]]. Betrag und Phase dieser Grössen lassen sich entweder als Zeiger oder als komplexe Zahl darstellen, was das Berechnen von Strom und Spannung stark erleichtert. Der Betrag der Impedanz, also das Verhältnis zwischen den Amplituden der Spannung und der Stromstärke, Scheinwiderstand genannt, kann bei Serie- oder Reihenschaltung mit Hilfe des Pythagoras gerechnet werden


:<math>Z=\sqrt{R^2+\left(\frac{1}{\hat\omega C}\right)^2-\left(\hat\omega L\right)^2}</math>
:<math>Z=\sqrt{R^2+\left(\frac{1}{\hat\omega C}\right)^2-\left(\hat\omega L\right)^2}</math>


Der stärkste Strom tritt dann auf, wenn sich die Blindwiderstände der Kapazität und der Indukditvität genau kompensieren. Dann begrenzt nur noch der Wirkwiderstand die Stromstärke.
Der stärkste Strom tritt dann auf, wenn sich die Blindwiderstände der Kapazität und der Induktivität genau kompensieren. Dann begrenzt nur noch der Wirkwiderstand die Stromstärke.


All diese Aussagen lassen sich auf mechanische Systeme übertragen, falls die zugehörigen Bauteile analog angeordnete sind. Oft sind Feder und Dämpfer aber parallel geschaltet (Impuls- bzw. Drehimpulsstrom verzweigen sich) und die Anregung erfolgt nicht mit konstanter Geschwindigkeits- oder Winkelgeschwindigkeitsamplitude. Dann muss man halt den dazu analogen elektrischen Kreis zur Untersuchung beiziehen. Die mechanisch-elektrische Analogie hat mehrere Vorteile. Erstens sind elektrische Systeme schneller aufgbeaut,reagieren schneller und sind billiger als mechanische. Zweitens stehen für elektrische Netzwerke leistungsstarke Simulatoren zur Verfügung. Drittens interessiert man sich in der Regelungstechnik nur für das Verhalten eines Systems und nicht für dessen physikalischen Eigenheiten.
All diese Aussagen lassen sich auf mechanische Systeme übertragen, falls die zugehörigen Bauteile analog angeordnet sind. Oft sind Feder und Dämpfer aber parallel geschaltet (Impuls- bzw. Drehimpulsstrom verzweigen sich) und die Anregung erfolgt nicht mit konstanter Geschwindigkeits- oder Winkelgeschwindigkeitsamplitude. Dann muss man halt den dazu analogen elektrischen Kreis zur Untersuchung beiziehen. Die mechanisch-elektrische Analogie hat mehrere Vorteile. Erstens sind elektrische Systeme schneller aufgebaut, reagieren schneller und sind billiger als mechanische. Zweitens stehen für elektrische Netzwerke leistungsstarke Simulatoren zur Verfügung. Drittens interessiert man sich in der Regelungstechnik nur für das Verhalten eines Systems und nicht für dessen physikalische Eigenheiten.


'''Video:''' [https://www.youtube.com/watch?v=BF1Wn40WlJk Resonanz]
==graphische Darstellungen==

Die [[System Dynamics|Systemdynamik]] ermöglicht eine intuitive Modellbildung. Ausgehend von der Mengenbilanz fügt man die konstitutiven Gesetze ein und verbindet das Ganze zu einem Wirkkreis, was dann mathematisch zu einer Differentialgleichung führt. Konkret berechnet man aus den Speichern die Potentiale und daraus dann wiederum die Stromstärken. Diese Methode, die auch oft Analogien benutz, hat auch ihre Grenzen. So muss die [[Kausalität]] immer festgelegt sein und es dürfen keine algebraischen Gleichungen formuliert werden (circular referende, algebraic loop). Auch sind die Analogien nicht immer offensichtlich. Dementsprechend kann man nicht unbesehen Resultate von einem Gebiet ins andere übertragen. Nehmen wir als Beispiel den elektrischen Schwingkreis mit Kondensator, Widerstand und Induktivität in Reihe. Die entsprechende mechanische Anordung sind Körper mit Feder und Dämper in Reihe. Damit das System geschlossen ist (keine Impulsaustausch mit der Erde), nehmen wir zwei frei bewegliche Gleiter einer Luftkissenbahn mit Feder und Dämpfer im selben Impulsfluss.
==Modellbildung==
Die [[System Dynamics|Systemdynamik]] ermöglicht eine intuitive Modellbildung. Ausgehend von der Mengenbilanz fügt man die konstitutiven Gesetze ein und verbindet das Ganze zu einem Wirkkreis, was dann mathematisch zu einer Differentialgleichung führt. Konkret berechnet man aus den Speichern die Potentiale und daraus dann wiederum die Stromstärken. Diese Methode, die auch oft Analogien benutz, hat auch ihre Grenzen. So muss die [[Kausalität]] immer festgelegt sein und es dürfen keine zirkuläre Abhängigkeiten(circular reference, algebraic loop) formuliert werden. Auch sind die Analogien nicht immer offensichtlich. Dementsprechend kann man nicht unbesehen Resultate von einem Gebiet ins andere übertragen. Nehmen wir als Beispiel den elektrischen Schwingkreis mit Kondensator, Widerstand und Induktivität in Reihe. Die entsprechende mechanische Anordnung sind Körper mit Feder und Dämpfer in Reihe. Damit das System geschlossen ist (keine Impulsaustausch mit der Erde), nehmen wir zwei frei bewegliche Gleiter einer Luftkissenbahn mit Feder und Dämpfer im selben Impulsfluss (siehe Video Schwingkreis: Analogie und Dualität).


===Modell===
===Modell===
Zeile 221: Zeile 233:
|<math>U=RI</math>
|<math>U=RI</math>
|<math>F=k\Delta v</math>
|<math>F=k\Delta v</math>
|''k'' entspricht einem Leitwert
|Dämpferkonstante ''k'' entspricht einem Leitwert
|-
|-
|Induktivität
|Induktivität
|<math>U=L\dot I</math>
|<math>U=L\dot I</math>
|<math>F=D\Delta x</math>
|<math>F=D\Delta x</math>
|Federkonstante entspricht einer reziproken Induktivität
|Federkonstante ''D'' entspricht einer reziproken Induktivität
|}
|}


Die Serischaltung von zwei Stromgliedern ist mit einem systemdynamischen Tool nicht ganz einfach zu modellieren (deshalb sind die Simulatoren für elektrische Netzwerke anders programmiert). Dazu muss die Kausalität des Widerstandselementes umgedreht werden, d.h. der Strom (Impulsstrom) bestimmt die Spannung (Geschwindigkeitsdifferenz) und diese muss vom kapazitiv berechneten Wert abgezogen werden (siehe Systemdiagramm rechts).
Die Serieschaltung von zwei Stromgliedern ist mit einem systemdynamischen Tool nicht ganz einfach zu modellieren (deshalb sind die Simulatoren für elektrische Netzwerke anders programmiert). Dazu muss die Kausalität des Widerstandselementes umgedreht werden, d.h. der Strom (Impulsstrom) bestimmt die Spannung (Geschwindigkeitsdifferenz) und diese muss vom kapazitiv berechneten Wert abgezogen werden.
'''Video:''' [https://www.youtube.com/watch?v=q1MuIE-TQDkSchwingkreis: Analogie und Dualität].


===Verhalten im Zeitbereich===
===Verhalten im Zeitbereich===
Zuerst übersetzen wir die oben aufgeführten elektrischen Grössen Kreisfrequen des ungedämpften Oszillars und Abklingzeit in mechanische Grössen
Nun übersetzen wir die Berechnungsformel für die Kreisfrequenz des nicht gedämpften Oszillators und die Abklingzeit von der Elektrizitätslehre in die Mechanik


:<math>\omega_0=\sqrt{\frac{1}{LC}}=\sqrt{\frac{D}{m}}</math>
:<math>\omega_0=\sqrt{\frac{1}{LC}}=\sqrt{\frac{D}{m}}</math>
Zeile 238: Zeile 251:
:<math>\tau=2\frac{L}{R}=2\frac{k}{D}</math>
:<math>\tau=2\frac{L}{R}=2\frac{k}{D}</math>


Dann gilt für die Kreisfrequenz des ungedämpften Oszillators
Dann gilt für die Kreisfrequenz des gedämpften Oszillators


:<math>\omega=\sqrt{\omega_0^2-\frac{4}{\tau^2}}=\sqrt{\frac{D}{m}-\frac{D^2}{k^2}}</math>
:<math>\omega=\sqrt{\omega_0^2-\frac{1}{\tau^2}}=\sqrt{\frac{D}{m}-\frac{D^2}{4k^2}}</math>


Die Kreisfrequenz ist gleich null, wenn die Dämpferkonstante gleich Wurzel aus Masse mal Federkonstante ist. Die ganze Analyse wäre mit komplexen Zahlen etwas übersichtlicher.
Die Kreisfrequenz ist gleich null, wenn die Dämpferkonstante gleich der Hälfte der Wurzel aus Masse mal Federkonstante ist (bitte nachrechnen; die ganze Analyse wäre mit komplexen Zahlen etwas übersichtlicher).

Zur Analyse der Resultate sollte man zuerst die Strom- und Potentialgrössen darstellen, also Kraft und Geschwindigkeit. Zur Kontrolle kann man dann noch das Zeitverhalten der gespeicherten Energie beiziehen. Bei mechanischen Schwingkreisen sind Dämpfer und Feder meist parallel statt wie hier in Reihe geschaltet, es liegt dann eine duale Analogie vor.


===Phasenraum===
===Phasenraum===
Statt alle Grössen als Funktion der Zeit darzustellen, kann man auch zwei Zustandgrössen (Topfgrössen in der Systemdynamik) gegeneinander auftragen. Man erhält dann eine [[Phasenraum]]darstellung. In der Translationsmechanik spannen [[Impuls]] und Ort den Phasenraum auf.

Die untenstehenden Bilder zeigen das Grundmodell des mechanischen Serie-Schwingkreises, die zugehörige Energieebene, das Geschwindigkeits- und das Impulsstromstärke-Zeit-Diagramm bei kritischer Dämpfung sowie die Phasenraumdarstellung für einen Körper bei drei verschiedenen Dämpfungen (unterkritisch bis kritisch).

<gallery>
MechSchwingkreisSerieBM.jpg|[[Systemdiagramm]] Grundmodell
MechSchwingkreisSerieEnergie.jpg|[[Systemdiagramm]] Energieebene
MechSchwingkreisSerie v F t.png|''v-t-'' und ''F-t-''Diagramm
MechSchwingkreisSerie p x.png|Phasenraum
</gallery>


==Systeme höherer Ordnung==
==Systeme höherer Ordnung==
Die Ordnung eines System ist durch die Zahl der unabhängigen Integrationen gegeben (zwei in Reihe geschaltete Kapazitäten oder Induktivitäten ergeben nur eine Ordnung). Weil sowohl Kapazität als auch Induktivität Energie speichern können, darf die Ordnung auch durch die Zahl der unabhängigen Energiespeicher ausgedrückt werden. Das Verhalten nichtlinearer Systeme (Blasenspeicher, Teller- oder Elastomerfedern, nichtlineare Dämpfer) kann in der Nähe eines Arbeitspunktes linearisiert werden. Dieses Verfahren wird oft in der Regeleungstechnik angewendet, weil man dem System nicht erlaubt, sich allzuweit von gewissen Zuständen zu entfernen. Federn sollte man aber nur dann als Induktivitäten behandeln, wenn deren Verhalten näherungsweise linear ist und die Analogie mit elektrischen oder hydraulischen Systemen gesucht wird. Ansonsten beschreibt man diese Elemente mit der Kraft-Verformungs-Funktion bzw. der Drehmoment-Verdrehungs-Funktion. Die folgende Tabelle führt nochmals analoge Grössen auf
Die Ordnung eines Systems ist durch die Zahl der unabhängigen Integrationen gegeben (zwei in Reihe geschaltete Kapazitäten oder Induktivitäten ergeben nur eine Ordnung). Weil sowohl Kapazität als auch Induktivität Energie speichern, darf die Ordnung auch durch die Zahl der unabhängigen Energiespeicher ausgedrückt werden. Das Verhalten nichtlinearer Systeme (Blasenspeicher, Teller- oder Elastomerfedern, nichtlineare Dämpfer) kann in der Nähe eines Arbeitspunktes linearisiert werden. Dieses Verfahren wird oft in der Regelungstechnik angewendet, weil man dem System nicht erlaubt, sich allzu weit von gewissen Zuständen zu entfernen. Federn sollte man aber nur dann als Induktivitäten behandeln, wenn deren Verhalten näherungsweise linear ist und die Analogie mit elektrischen oder hydraulischen Systemen gesucht wird. Ansonsten beschreibt man diese Elemente mit der Kraft-Verformungs-Funktion bzw. der Drehmoment-Verdrehungs-Funktion. Die folgende Tabelle führt nochmals analoge Grössen auf


{|class="wikitable"
{|class="wikitable"
Zeile 290: Zeile 315:
Ein Spannungsstoss (Integral der Spannung über die Zeit oder Fläche unter der Spannungs-Zeit-Kurve) ist oft die Folge einer magnetischen Flussänderung. Kraft oder Drehmoment sind Impuls- bzw. Drehimpulsstromstärken bezüglich eines Körpers. Wird eine Feder von einem Impulsstrom durchflossen, kann man die zugehörige Stromstärke am Ein- oder Ausgang messen. So erhält man die beiden Kräfte, welche die Feder im Gleichgewicht halten.
Ein Spannungsstoss (Integral der Spannung über die Zeit oder Fläche unter der Spannungs-Zeit-Kurve) ist oft die Folge einer magnetischen Flussänderung. Kraft oder Drehmoment sind Impuls- bzw. Drehimpulsstromstärken bezüglich eines Körpers. Wird eine Feder von einem Impulsstrom durchflossen, kann man die zugehörige Stromstärke am Ein- oder Ausgang messen. So erhält man die beiden Kräfte, welche die Feder im Gleichgewicht halten.


Systeme n-ter Ordnung werden mit einer Differentialgleichung der entsprechenden Ordnung beschrieben. Bei linearen, zeitinvarianten Systemen führt dies zu einer Differentialgleichung n-ter Ordnung mit konstanten Koeffizienten. Solche Systeme lassen sich auch und in gewissem Sinne anachaulicher mit n Differentialgleichungen erster Ordnung beschreiben. Die [[Systemphysik]] erlaubt nun auch in der Mechanik eine direkte Formulieren mittels Gleichungen erster Ordnung. Dabei liefert die Impulsbilanz bzw. die Drehimpulsbilanz für jedes System und jede Komponente eine Gleichung, die zugehörige Kinematik steuert dann den zweiten Partner bei. Die meisten Modelle, die Sie mit BerkeleyMadonna augebaut haben, erzeugen solche Gleichungen. So gesehen haben Sie in der Zustandsraumdarstellung modelliert.
Systeme n-ter Ordnung werden mit einer Differentialgleichung der entsprechenden Ordnung beschrieben. Bei linearen, zeitinvarianten Systemen führt dies zu einer Differentialgleichung n-ter Ordnung mit konstanten Koeffizienten. Solche Systeme lassen sich auch und in gewissem Sinne anschaulicher mit n Differentialgleichungen erster Ordnung beschreiben. Die [[Systemphysik]] erlaubt eine direkte Modellierung mechanischer Systeme mittels Gleichungen erster Ordnung. Die Impulsbilanz bzw. die Drehimpulsbilanz liefert so für jedes System und jede Komponente eine Gleichung, die zugehörige Kinematik steuert dann die zweite Gleichung bei. Die meisten Modelle, die Sie mit BerkeleyMadonna aufgebaut haben, erzeugen solche Gleichungen. So gesehen haben Sie ab erster Woche in der Zustandsraumdarstellung modelliert.


==Kontrollfragen==
==Kontrollfragen==
#Ein Körper schwingt an einer Feder auf und ab. Wie gross ist die Schwingungsdauer, wenn Masse und Federkonstante bekannt sind?
#Zwei Schwungräder, die über eine elastische Welle miteinander verbunden und reibungsfrei gelagert sind, schwingen gegeneinander. Wie berechnet man die Winkelrichtgrösse (Drehfederkonstante der Welle) aus Schwingungsdauer und den beiden Massenträgheitsmomenten?
#Wie gross muss der Widerstand in einem elektrischen Schwingkreis (Kapazität ''C'', Induktivität ''L'') gewählt werden, damit der Strom kritisch gedämpft abklingt?
#Ein gefedertes Fahrzeug belastet die vier Räder gleichmässig. Wie berechnet man die Dämpferkonstante für die kritische Dämpfung aus der Masse des Fahrzeugs und der Konstante für die Federung eines Rades (lineares Verhalten vorausgesetzt)?
#Um einen geraden Stahldraht um eine halbe Umdrehung zu verdrehen, muss man einen Drehimpulsstrom der Stärke 5 Nm durchfliessen lassen. Nun hängt man den Draht an der Decke auf und befestigt am unteren Ende einen Körper. Wie gross ist dessen Massenträgheitsmoment, wenn er für eine volle Drehschwingung 0.8 s benötigt?
#Zwei Schwungräder sind über eine elastische Welle miteinander verbunden. Zeichnen Sie ein elektrisch analoges Schaltbild unter der Annahme, dass die Reibung in beiden Lagern proportional zu Winkelgeschwindigkeit ist.
#Ein Körper liegt auf einem Feder-Dämpfer-System (parallel geschaltet) auf. Skizzieren Sie die elektrisch analoge Schaltung. Nun fügen sie in diesem Schaltbild eine Wechselspannungsquelle ein und überliegen sich, wie das mechanische Analogon aussehen müsste.


==Antworten zu den Kontrollfragen==
==Antworten zu den Kontrollfragen==
#Schwingungsdauer <math>T=\frac{2\pi}{\omega}=2\pi\sqrt{\frac{m}{D}}</math>
#Schwingungsdauer <math>T=2\pi\sqrt{\frac{J}{D^*}}</math> also gilt <math>D^*=J\frac{4\pi^2}{T^2}</math>
#Kritische Dämpfung tritt ein, wenn die Kreisfrequenz null ist, wenn gilt <math>\sqrt{\omega_0^2-\frac{1}{\tau^2}}=\sqrt{\frac{1}{LC}-\frac{R^2}{4L^2}}=0</math>. Somit gilt für den Widerstand bei kritischer Dämpfung <math>R=2\sqrt\frac{L}{C}</math>
#Übersetzen wir die Formel für die kritische Dämpfung vom elektrischen auf den mechanischen Schwingkreis <math>k=\frac{1}{2}\sqrt{mD}</math> (die Dämpferkonstante entspricht einem reziproken Widerstand, also einem Leitwert, die Federkonstante der reziproken Induktivität und die Masse der Kapazität). Nun verteilt sich der Impulsstrom zwischen Fahrzeug und Erde auf vier identische Stromglieder. Also ist nur ein Viertel der Masse zu nehmen <math>k=\frac{1}{4}\sqrt{m_{Fahrzeug}D}</math>
#<math>J=D^*\frac{T^2}{4\pi^2}=\frac{5 Nm}{\pi}\frac{(0.8 s)^2}{4\pi^2}=2.58 10^{-2} kgm^2</math>
#Zwei Kondensatoren einseitig geerdet mit Induktivität in erster Verbindung und zwei Widerständen in der zweiten Verbindung. Ein Punkt zwischen den beiden Widerständen ist geerdet.
#Kondensator und Wechselspannungsquelle einseitig geerdet. Induktivität und Widerstand parallel mit den andern beiden Enden von Kondensator und Wechselspannungsquelle verbunden. Mechanische Anregung erfolgt vom Boden her. Das ist ein einfaches Modell für ein Auto, das über Bodenwellen fährt. Die Gewichtskraft wäre als Konstantstromquelle zu modellieren, was aber keinen Einfluss auf die Dynamik hat.


'''[[Physik und Systemwissenschaft in Aviatik 2014]]'''
'''[[Physik und Systemwissenschaft in Aviatik 2014]]'''

Aktuelle Version vom 18. Mai 2017, 13:36 Uhr

Lernziele

In dieser Vorlesung lernen Sie

  • wie man die Kreisfrequenz, Frequenz und Schwingungsdauer eines harmonischen Oszillators berechnet.
  • wie die Potential-Zeit-, die Stromstärke-Zeit- und die Leistungs-Zeit-Funktion bei einem harmonischen Oszillator formuliert werden.
  • wie man die Abklingzeit und die Kreisfrequenz bei einem gedämpften Oszillator berechnet.
  • unter welchen Bedingungen ein Serie-Schwingkreis bei harmonischer Anregung den stärksten Strom durchlässt
  • wie das Verhalten dynamischer Systeme modelliert und analysiert wird.

Problemstellung

Zwei zylinderförmige Gefässe sind über ein bodennahes Röhrchen miteinander verbunden. Solange das eine Gefäss höher als das andere mit Wasser gefüllt ist, fliesst ein Ausgleichsstrom. Nimmt man Öl statt Wasser, bleibt die Strömung laminar und wir können den Ausgleichsvorgang mit einer Exponentialfunktion beschreiben. Nun denken wir uns das Röhrchen immer dicker bis ein mit Wasser gefülltes U-Rohr vor uns steht. Dann erfolgt der Ausgleichsvorgang nicht mehr über einen einmal abklingenden Volumenstrom, sondern über einen abklingend-oszillierenden Volumenstrom. Wäre die Flüssigkeit suprafluid, würde die Oszillation oder Schwingung überhaupt nicht mehr abklingen.

Das Verhalten des U-Rohrs kann durch zwei Speicher mit konstanter Kapazität und einem Leiter mit Widerstand und Induktivität modelliert werden. Analoge Systeme kennt man aus der Elektrodynamik (zwei über eine Spule verbunden Kondensatoren), aus der Translationsmechanik (zwei über eine Feder verbundene Luftkissenfahrzeuge) und aus der Rotationsmechanik (zwei über eine Drehfeder verbundene Schwungräder). Dominiert das Widerstandselement (Strömungswiderstand, elektrischer Widerstand der Spule, Dämpfer in Serie zur Feder), klingt der Potentialunterschied (Druck, Spannung, Geschwindigkeit, Winkelgeschwindigkeit) exponentiell ab, falls sich alle Elemente linear verhalten.

Dynamischen Systemen, die aus zwei Speichern bestehen und die über einen Leiter miteinander verbunden sind, können mathematisch zu einem System mit nur einem Speicher und reduzierter Kapazität zusammenfasst werden. Systeme mit Mengenspeicher und Induktivität nennt man Systeme 2. Ordnung, weil man zur Berechnung des Verhaltens zweimal integrieren muss. Systeme 2. Ordnung besitzen zwei unabhängige Energiespeicher.

Speicher und Leiter

Speicher mit konstanter Kapazität und Leiter mit linearem Widerstand und linearer Induktivität findet man in der Hydrodynamik, der Elektrodynamik, der Translationsmechanik und der Rotationsmechanik.

Gebiet Speicher Widerstand Induktivität Bemerkung
Hydrodynamik [math]\Delta p=\frac{\Delta V}{C_V}=\frac{\varrho g\Delta V}{A}[/math] [math]\Delta p=R_VI_V[/math] [math]\Delta p=L_V\dot I_V[/math] laminare Strömung
Elektrodynamik [math]U=\frac{Q}{C}[/math] oder [math]\dot U=\frac{I}{C}[/math] [math]U=RI[/math] [math]U=L\dot I[/math] Kondensator
Translationsmechanik [math]v_x=\frac{p_x}{m}[/math] [math]\Delta v_x=R_{px}F_x[/math] [math]\Delta v_x=L_{px}\dot F_x[/math] Federgesetz
Rotationsmechanik [math]\omega_x=\frac{L_x}{J_{xx}}[/math] [math]\Delta \omega_x=R_{Lx}M_x[/math] [math]\Delta \omega_x=L_{Lx}\dot M_x[/math] Drehfedergesetz

Die mechanischen Systemparameter sind oft reziprok definiert: statt eines Impulswiderstandes führt man eine Dämpferkonstante ein; anstelle einer Impulsinduktivität definiert man eine Federkonstante

[math]F_x=k\Delta v_x[/math] statt [math]\Delta v_x= R_{px}F_x[/math]
[math]F_x=D\Delta x[/math] statt [math]\Delta v_x=L_{px}\dot F_x[/math]

Harmonischer Oszillator

Ein harmonischer Oszillator ist ein schwingungsfähiges System, das aus einem oder zwei Speicher mit konstanter Kapazität und einem Leiter mit stromunabhängiger Induktivität besteht. Lenkt man ein solches System aus und lässt es los, schwingt es sinusförmig (harmonisch) um seine Ruhelage, wobei die Schwingungsdauer unabhängig von der Grösse der Auslenkung ist. Beispiele für harmonische Oszillatoren sind Federpendel, elektrische Schwingkreise und U-Rohr mit sehr grossem Querschnitt.

Der harmonische Oszillator ist ein wichtiges Modellsystem der Physik. Er ist durch nur zwei Parameter (Gesamtkapazität und die Induktivität) vollständig beschrieben. Viele komplexere Systeme verhalten sich bei kleinen Auslenkungen näherungsweise wie harmonische Oszillatoren. Der harmonische Oszillator der Quantenmechanik ist eines der wenigen quantenmechanischen Systeme, das sich ohne Näherungen berechnen lässt.

In der Elektrodynamik bildet ein über eine ideale Spule (Induktivität L) kurz geschlossener Kondensator (Kapazität C) einen harmonischen Oszillator. Lädt man den Kondensator auf die Spannung U0 auf und verbindet ihn dann mit der idealen Spule, ist die Spannung über beiden Elementen zu jedem Zeitpunkt gleich gross (die Umlaufspannung muss immer gleich null sein)

[math]U_C+U_L=0[/math]

Setzt man für die beiden Spannungen die zugehörigen konstitutiven Gesetze ein, erhält man

[math]\frac{Q}{C}+L\dot I=0[/math]

Nun ersetzt man noch die Stromstärke über die Bilanz durch die Änderungsrate der Ladung

[math]Q+LC\ddot Q=0[/math]

Die allgemeine Lösung dieser linearen Differentialgleichung 2. Ordnung mit konstanten Koeffizienten ist

[math]Q(t)=Q_0cos(\omega_0t+\delta)[/math]

wobei für die Kreisfrequenz gilt (Lösung in Differenzialgleichung einsetzen und vergleichen!)

[math]\omega_0=\sqrt{\frac{1}{CL}}[/math]

Die Ladungsamplitude Q0 und die Phasenverschiebung hängen von den Anfangsbedingungen ab. Mit Hilfe des kapazitiven Gesetzes kann auf die Spannungs-Zeit-Funktion umgerechnet werden

[math]U(t)=U_0cos(\omega_0t+\delta)[/math] mit [math]U_0=\frac{Q_0}{C}[/math]

Ersetzt man die Spannung mit Hilfe des induktiven Gesetzes erhält man eine Funktion für die Änderungsrate der Stromstärke

[math]\dot I(t)=\frac{U_0}{L}cos(\omega_0t+\delta)[/math].

Eine Integration über die Zeit liefert

[math]I(t)=I_0sin(\omega_0t+\delta)[/math] mit [math]I_0=\frac{U_0}{\omega_0L}=U_0\sqrt{\frac{C}{L}}=Q_0\sqrt{\frac{1}{CL}}=\omega_0 Q_0[/math]

Für die Prozessleistung über dem Kondensator oder der Spule gilt

[math]P=UI=U_0I_0cos(\omega_0t+\delta)sin(\omega_0t+\delta)[/math]

Setzt man die Phasenverschiebung δ gleich null (freie Wahl des Zeitnullpunktes), erhält man eine Leistung, die mit doppelter Frequenz (Kreisfrequenz durch 2π) schwingt

[math]P=UI=\frac{U_0I_0}{2}sin(2\omega_0t)[/math]

Während einer Periode oder Schwingungsdauer (2π durch Kreisfrequenz) ist der Kondensator zweimal geladen und der Strom erreicht zweimal seine maximale Stärke. Weil beide Bauteile zusammen mit der Ladung (Kondensator) oder mit dem Strom (Spule) Energie speichern, wird die Energie mit doppelter Oszillatorfrequenz hin und her verschoben.

Video: harmonischer Oszillator

gedämpfter Oszillator

Nun fügen wir beim elektrischen Schwingkreis noch einen Widerstand ein. Widerstand und Induktivität können als einfache Ersatzschaltung für eine reale Spule gesehen werden. Wieder gilt der Maschensatz, wonach Umlaufspannung gleich null sein muss

[math]U_C+U_R+U_L=0[/math]

Setzt man wie beim harmonischen Oszillator die konstitutiven Gesetze für die Kapazität, den Widerstand und die Induktivität ein und ersetzt die Stromstärke über die Bilanzgleichung durch die Änderungsrate der Kondensatorladung, erhält man eine lineare Differentialgleichung 2. Ordnung mit konstanten Koeffizienten

[math]\frac{Q}{C}+R\dot Q+L\ddot Q=0[/math] oder [math]\frac{1}{LC}Q+\frac{R}{L}\dot Q+\ddot Q=0[/math]

Diese Gleichung wird durch folgende allgemeine Funktion erfüllt

[math]Q(t)=Q_0e^{-t/\tau}cos(\omega t+\delta)[/math]

Setzt man diese Funktion in die Differentialgleichung ein, erhält man folgende Beziehung

[math]\left(\left(\frac{1}{LC}-\frac{R}{L\tau}+\frac{1}{\tau^2}-\omega^2\right)cos(\omega t+\delta)+\left(-\frac{R\omega}{L}+\frac{2\omega}{\tau}\right)sin(\omega t+\delta)\right)Q_0e^{-t/\tau}=0[/math]

Die linke Seite dieser Gleichung kann nur zu allen Zeitpunkten gleich null sein, wenn beide Klammerausdrücke den Wert null annehmen. Damit gilt für die Abklingzeit

[math]\tau =2\frac{L}{R}[/math]

und für die Kreisfrequenz

[math]\omega=\sqrt{\frac{1}{LC}-\frac{R^2}{4L^2}}=\sqrt{\omega_0^2-\frac{1}{\tau^2}}[/math]

Ohne Dämpfung ist die Abklingzeit unendlich gross und die Kreisfrequenz entspricht dem Wert des harmonischen Oszillators. Unterschreitet die Abklingzeit einen gewissen Wert, wird der Ausdruck unter der Wurzel negativ und das System schwingt gar nicht mehr, sondern kriecht - wie wenn die Induktivität nicht da wäre - gegen die Gleichgewichtslage. Man unterscheidet drei Dämpfungsarten

  • [math]\omega=\sqrt{\omega_0^2-\frac{1}{\tau^2}}\gt 0[/math]: unterkritische Dämpfung: System schwingt bezüglich der Gleichgewichtslage hin und her
  • [math]\omega=\sqrt{\omega_0^2-\frac{1}{\tau^2}}=0[/math]: kritische Dämpfung: System bewegt sich gegen die Gleichgewichtslage ohne zu überschwingen
  • [math]\omega=\sqrt{\omega_0^2-\frac{1}{\tau^2}}\lt 0[/math]: überkritische Dämpfung: System kriecht gegen Gleichgewichtslage

Die nachfolgende Tabelle zeigt noch die Analogie zwischen elektrischen und mechanischen Schwingkreisen

Thema Elektrodynamik Hydrodynamik Translation Rotation
Serie-Schwingkreis Kondensator, reale Spule U-Rohr Körper, Feder, Dämpfer, Körper Schwungrad, Feder, Dämpfer, Schwungrad
Systemparameter Kapazität C [math]C_V=\frac{A}{\varrho g}[/math] [math]C_{px}=m[/math] [math]C_{Lx}=J_{xx}[/math]
Widerstand R [math]R_V=\frac{8\pi\cdot Länge\cdot\eta}{A^2}[/math] [math]R_{px}=\frac{1}{k}[/math] [math]R_{Lx}=\frac{1}{k^*}[/math]
Induktivität L [math]L_V=\frac{\varrho A}{Länge}[/math] [math]L_{px}=\frac{1}{D}[/math] [math]L_{Lx}=\frac{1}{D^*}[/math]
Kreisfrequenz [math]\omega_0[/math] [math]\omega_0=\sqrt{\frac{1}{LC}}[/math] [math]\omega_0=\sqrt{\frac{1}{L_VC_V}}=\sqrt{\frac{g}{2\cdot Länge}}[/math] [math]\omega_0=\sqrt{\frac{D(m_1+m_2)}{m_1m_2}}[/math] [math]\omega_0=\sqrt{\frac{D^*(J_1+J_2)}{J_1J_2}}[/math]
Abklingzeit [math]\tau=\frac{2L}{R}[/math] [math]\tau=\frac{2L_V}{R_V}[/math] [math]\tau=\frac{2k}{D}[/math] [math]\tau=\frac{2k^*}{D^*}[/math]

k ist die Dämpferkonstante, D die Federkonstante oder Richtgrösse. Bei der entsprechenden Drehgrösse ist noch ein Stern zugefügt. In der Mechanik sind Feder und Dämpfer oft parallel statt in Reihe geschaltet. Der Strom verzweigt sich und man spricht man von einem Parallel-Schwingkreis. Das U-Rohr verhält sich nur bedingt analog zu den andern Systemen, weil die Strömung eher turbulent als laminar ist und das System sich damit nichtlinear verhält.

Video: Schwingkreis: Analogie und Dualität

Resonanz

Ein elektrischer Schwingkreis kann mit einer Wechselspannungsquelle, die man in den Kreis einfügt, angeregt werden. Alle vier Elemente (Spannungsquelle, Kapazität, Widerstand und Induktivität) sind miteinander zu einem nicht verzweigten Kreis verbunden. Die beschreibende Gleichung lautet

[math]U_C+U_R+U_L=\hat Ucos(\hat{\omega}t)[/math]

Das System durchläuft zuerst einen Einschwingvorgang und geht dann in ein stationäres Verhalten über. Der Einschwingvorgang hängt von den Anfangsbedingungen ab. Oft interessiert man sich nur für das stationäre Verhalten nach dem Einschwingvorgang. Dann übernimmt die Spannungsquelle das Zepter und das ganze System schwingt mit [math]\hat{\omega}=2\pi\hat f=\frac{2\pi}{\hat T}[/math].

Weil sich alle Elemente im selben Kreis befinden, ist der elektrische Strom überall gleich stark. Im stationären Zustand schwingt der Strom mit der gleichen Frequenz wie die Spannung. Einzig eine Phasenverschiebung ist noch möglich. Diese Überlegung bringt uns zu folgendem Ansatz

[math]I(t)=\hat Icos(\hat\omega t+\delta)[/math]

Nun können wir in der Differentialgleichung die Spannungen mit Hilfe der konstitutiven Gesetze durch die Stromstärken ersetzen

[math]\dot U_C=\frac{I}{C}[/math] also [math]U_C=\frac{\hat I}{\hat\omega C}sin(\hat{\omega}t+\delta)[/math]
[math]U_R=RI[/math] also [math]U_R=R\hat Icos(\hat{\omega}t+\delta))[/math]
[math]U_L=L\dot I[/math] also [math]U_L=-\hat\omega Lsin(\hat{\omega}t+\delta))[/math]

Würde man jedes der drei Elemente einzeln mit der Spannungsquelle verbinden, wäre die Lösung der Gleichung einfach zu finden: beim Widerstand schwingt der Strom in Phase mit der angelegten Spannung, bei der Kapazität läuft der Strom der Spannung eine Viertelperiode voraus (der Strom baut die Spannung auf) und bei der Induktivität eine Viertelperiode nach (Spannung ändert Stromstärke).

Strom- und Spannungsamplituden unterscheiden sich bei allen drei Elementen um einen Faktor, der entweder konstant ist, sich proportional oder reziprok-proportional mit der angelegten Kreisfrequenz ändert. Dieser Faktor, eigentlich eine Verallgemeinerung des Widerstandes, heisst Impedanz Z. Die Verallgemeinerung des Leitwerts, des Kehrwerts der Impedanz, nennt man Admittanz. Betrag und Phase dieser Grössen lassen sich entweder als Zeiger oder als komplexe Zahl darstellen, was das Berechnen von Strom und Spannung stark erleichtert. Der Betrag der Impedanz, also das Verhältnis zwischen den Amplituden der Spannung und der Stromstärke, Scheinwiderstand genannt, kann bei Serie- oder Reihenschaltung mit Hilfe des Pythagoras gerechnet werden

[math]Z=\sqrt{R^2+\left(\frac{1}{\hat\omega C}\right)^2-\left(\hat\omega L\right)^2}[/math]

Der stärkste Strom tritt dann auf, wenn sich die Blindwiderstände der Kapazität und der Induktivität genau kompensieren. Dann begrenzt nur noch der Wirkwiderstand die Stromstärke.

All diese Aussagen lassen sich auf mechanische Systeme übertragen, falls die zugehörigen Bauteile analog angeordnet sind. Oft sind Feder und Dämpfer aber parallel geschaltet (Impuls- bzw. Drehimpulsstrom verzweigen sich) und die Anregung erfolgt nicht mit konstanter Geschwindigkeits- oder Winkelgeschwindigkeitsamplitude. Dann muss man halt den dazu analogen elektrischen Kreis zur Untersuchung beiziehen. Die mechanisch-elektrische Analogie hat mehrere Vorteile. Erstens sind elektrische Systeme schneller aufgebaut, reagieren schneller und sind billiger als mechanische. Zweitens stehen für elektrische Netzwerke leistungsstarke Simulatoren zur Verfügung. Drittens interessiert man sich in der Regelungstechnik nur für das Verhalten eines Systems und nicht für dessen physikalische Eigenheiten.

Video: Resonanz

Modellbildung

Die Systemdynamik ermöglicht eine intuitive Modellbildung. Ausgehend von der Mengenbilanz fügt man die konstitutiven Gesetze ein und verbindet das Ganze zu einem Wirkkreis, was dann mathematisch zu einer Differentialgleichung führt. Konkret berechnet man aus den Speichern die Potentiale und daraus dann wiederum die Stromstärken. Diese Methode, die auch oft Analogien benutz, hat auch ihre Grenzen. So muss die Kausalität immer festgelegt sein und es dürfen keine zirkuläre Abhängigkeiten(circular reference, algebraic loop) formuliert werden. Auch sind die Analogien nicht immer offensichtlich. Dementsprechend kann man nicht unbesehen Resultate von einem Gebiet ins andere übertragen. Nehmen wir als Beispiel den elektrischen Schwingkreis mit Kondensator, Widerstand und Induktivität in Reihe. Die entsprechende mechanische Anordnung sind Körper mit Feder und Dämpfer in Reihe. Damit das System geschlossen ist (keine Impulsaustausch mit der Erde), nehmen wir zwei frei bewegliche Gleiter einer Luftkissenbahn mit Feder und Dämpfer im selben Impulsfluss (siehe Video Schwingkreis: Analogie und Dualität).

Modell

Zuerst zur Modellanalogie

Element elektrisch mechanisch Bemerkung
Kapazität [math]Q=CU[/math] [math]p=mv[/math] Gesamtladung des Kondensators ist null
Widerstand [math]U=RI[/math] [math]F=k\Delta v[/math] Dämpferkonstante k entspricht einem Leitwert
Induktivität [math]U=L\dot I[/math] [math]F=D\Delta x[/math] Federkonstante D entspricht einer reziproken Induktivität

Die Serieschaltung von zwei Stromgliedern ist mit einem systemdynamischen Tool nicht ganz einfach zu modellieren (deshalb sind die Simulatoren für elektrische Netzwerke anders programmiert). Dazu muss die Kausalität des Widerstandselementes umgedreht werden, d.h. der Strom (Impulsstrom) bestimmt die Spannung (Geschwindigkeitsdifferenz) und diese muss vom kapazitiv berechneten Wert abgezogen werden. Video: Analogie und Dualität.

Verhalten im Zeitbereich

Nun übersetzen wir die Berechnungsformel für die Kreisfrequenz des nicht gedämpften Oszillators und die Abklingzeit von der Elektrizitätslehre in die Mechanik

[math]\omega_0=\sqrt{\frac{1}{LC}}=\sqrt{\frac{D}{m}}[/math]
[math]\tau=2\frac{L}{R}=2\frac{k}{D}[/math]

Dann gilt für die Kreisfrequenz des gedämpften Oszillators

[math]\omega=\sqrt{\omega_0^2-\frac{1}{\tau^2}}=\sqrt{\frac{D}{m}-\frac{D^2}{4k^2}}[/math]

Die Kreisfrequenz ist gleich null, wenn die Dämpferkonstante gleich der Hälfte der Wurzel aus Masse mal Federkonstante ist (bitte nachrechnen; die ganze Analyse wäre mit komplexen Zahlen etwas übersichtlicher).

Zur Analyse der Resultate sollte man zuerst die Strom- und Potentialgrössen darstellen, also Kraft und Geschwindigkeit. Zur Kontrolle kann man dann noch das Zeitverhalten der gespeicherten Energie beiziehen. Bei mechanischen Schwingkreisen sind Dämpfer und Feder meist parallel statt wie hier in Reihe geschaltet, es liegt dann eine duale Analogie vor.

Phasenraum

Statt alle Grössen als Funktion der Zeit darzustellen, kann man auch zwei Zustandgrössen (Topfgrössen in der Systemdynamik) gegeneinander auftragen. Man erhält dann eine Phasenraumdarstellung. In der Translationsmechanik spannen Impuls und Ort den Phasenraum auf.

Die untenstehenden Bilder zeigen das Grundmodell des mechanischen Serie-Schwingkreises, die zugehörige Energieebene, das Geschwindigkeits- und das Impulsstromstärke-Zeit-Diagramm bei kritischer Dämpfung sowie die Phasenraumdarstellung für einen Körper bei drei verschiedenen Dämpfungen (unterkritisch bis kritisch).

Systeme höherer Ordnung

Die Ordnung eines Systems ist durch die Zahl der unabhängigen Integrationen gegeben (zwei in Reihe geschaltete Kapazitäten oder Induktivitäten ergeben nur eine Ordnung). Weil sowohl Kapazität als auch Induktivität Energie speichern, darf die Ordnung auch durch die Zahl der unabhängigen Energiespeicher ausgedrückt werden. Das Verhalten nichtlinearer Systeme (Blasenspeicher, Teller- oder Elastomerfedern, nichtlineare Dämpfer) kann in der Nähe eines Arbeitspunktes linearisiert werden. Dieses Verfahren wird oft in der Regelungstechnik angewendet, weil man dem System nicht erlaubt, sich allzu weit von gewissen Zuständen zu entfernen. Federn sollte man aber nur dann als Induktivitäten behandeln, wenn deren Verhalten näherungsweise linear ist und die Analogie mit elektrischen oder hydraulischen Systemen gesucht wird. Ansonsten beschreibt man diese Elemente mit der Kraft-Verformungs-Funktion bzw. der Drehmoment-Verdrehungs-Funktion. Die folgende Tabelle führt nochmals analoge Grössen auf

Grösse Elektrodynamik Hydrodynamik Translation Rotation
Menge elektrische Ladung Volumen Impuls Drehimpuls
Stromstärke elektrischer Strom Volumenstrom Kraft Drehmoment
Potential Spannung Druck Geschwindigkeit Winkelgeschwindigkeit
Extensum Spannungsstoss Druckstoss Ort Winkel
Prozessleistung [math]P=UI[/math] [math]P=\Delta pI_V[/math] [math]P=\Delta v_xF_x[/math] [math]P=\Delta\omega_xM_x[/math]

Ein Spannungsstoss (Integral der Spannung über die Zeit oder Fläche unter der Spannungs-Zeit-Kurve) ist oft die Folge einer magnetischen Flussänderung. Kraft oder Drehmoment sind Impuls- bzw. Drehimpulsstromstärken bezüglich eines Körpers. Wird eine Feder von einem Impulsstrom durchflossen, kann man die zugehörige Stromstärke am Ein- oder Ausgang messen. So erhält man die beiden Kräfte, welche die Feder im Gleichgewicht halten.

Systeme n-ter Ordnung werden mit einer Differentialgleichung der entsprechenden Ordnung beschrieben. Bei linearen, zeitinvarianten Systemen führt dies zu einer Differentialgleichung n-ter Ordnung mit konstanten Koeffizienten. Solche Systeme lassen sich auch und in gewissem Sinne anschaulicher mit n Differentialgleichungen erster Ordnung beschreiben. Die Systemphysik erlaubt eine direkte Modellierung mechanischer Systeme mittels Gleichungen erster Ordnung. Die Impulsbilanz bzw. die Drehimpulsbilanz liefert so für jedes System und jede Komponente eine Gleichung, die zugehörige Kinematik steuert dann die zweite Gleichung bei. Die meisten Modelle, die Sie mit BerkeleyMadonna aufgebaut haben, erzeugen solche Gleichungen. So gesehen haben Sie ab erster Woche in der Zustandsraumdarstellung modelliert.

Kontrollfragen

  1. Ein Körper schwingt an einer Feder auf und ab. Wie gross ist die Schwingungsdauer, wenn Masse und Federkonstante bekannt sind?
  2. Zwei Schwungräder, die über eine elastische Welle miteinander verbunden und reibungsfrei gelagert sind, schwingen gegeneinander. Wie berechnet man die Winkelrichtgrösse (Drehfederkonstante der Welle) aus Schwingungsdauer und den beiden Massenträgheitsmomenten?
  3. Wie gross muss der Widerstand in einem elektrischen Schwingkreis (Kapazität C, Induktivität L) gewählt werden, damit der Strom kritisch gedämpft abklingt?
  4. Ein gefedertes Fahrzeug belastet die vier Räder gleichmässig. Wie berechnet man die Dämpferkonstante für die kritische Dämpfung aus der Masse des Fahrzeugs und der Konstante für die Federung eines Rades (lineares Verhalten vorausgesetzt)?
  5. Um einen geraden Stahldraht um eine halbe Umdrehung zu verdrehen, muss man einen Drehimpulsstrom der Stärke 5 Nm durchfliessen lassen. Nun hängt man den Draht an der Decke auf und befestigt am unteren Ende einen Körper. Wie gross ist dessen Massenträgheitsmoment, wenn er für eine volle Drehschwingung 0.8 s benötigt?
  6. Zwei Schwungräder sind über eine elastische Welle miteinander verbunden. Zeichnen Sie ein elektrisch analoges Schaltbild unter der Annahme, dass die Reibung in beiden Lagern proportional zu Winkelgeschwindigkeit ist.
  7. Ein Körper liegt auf einem Feder-Dämpfer-System (parallel geschaltet) auf. Skizzieren Sie die elektrisch analoge Schaltung. Nun fügen sie in diesem Schaltbild eine Wechselspannungsquelle ein und überliegen sich, wie das mechanische Analogon aussehen müsste.

Antworten zu den Kontrollfragen

  1. Schwingungsdauer [math]T=\frac{2\pi}{\omega}=2\pi\sqrt{\frac{m}{D}}[/math]
  2. Schwingungsdauer [math]T=2\pi\sqrt{\frac{J}{D^*}}[/math] also gilt [math]D^*=J\frac{4\pi^2}{T^2}[/math]
  3. Kritische Dämpfung tritt ein, wenn die Kreisfrequenz null ist, wenn gilt [math]\sqrt{\omega_0^2-\frac{1}{\tau^2}}=\sqrt{\frac{1}{LC}-\frac{R^2}{4L^2}}=0[/math]. Somit gilt für den Widerstand bei kritischer Dämpfung [math]R=2\sqrt\frac{L}{C}[/math]
  4. Übersetzen wir die Formel für die kritische Dämpfung vom elektrischen auf den mechanischen Schwingkreis [math]k=\frac{1}{2}\sqrt{mD}[/math] (die Dämpferkonstante entspricht einem reziproken Widerstand, also einem Leitwert, die Federkonstante der reziproken Induktivität und die Masse der Kapazität). Nun verteilt sich der Impulsstrom zwischen Fahrzeug und Erde auf vier identische Stromglieder. Also ist nur ein Viertel der Masse zu nehmen [math]k=\frac{1}{4}\sqrt{m_{Fahrzeug}D}[/math]
  5. [math]J=D^*\frac{T^2}{4\pi^2}=\frac{5 Nm}{\pi}\frac{(0.8 s)^2}{4\pi^2}=2.58 10^{-2} kgm^2[/math]
  6. Zwei Kondensatoren einseitig geerdet mit Induktivität in erster Verbindung und zwei Widerständen in der zweiten Verbindung. Ein Punkt zwischen den beiden Widerständen ist geerdet.
  7. Kondensator und Wechselspannungsquelle einseitig geerdet. Induktivität und Widerstand parallel mit den andern beiden Enden von Kondensator und Wechselspannungsquelle verbunden. Mechanische Anregung erfolgt vom Boden her. Das ist ein einfaches Modell für ein Auto, das über Bodenwellen fährt. Die Gewichtskraft wäre als Konstantstromquelle zu modellieren, was aber keinen Einfluss auf die Dynamik hat.

Physik und Systemwissenschaft in Aviatik 2014