Warum echte Schwingungen als A=A0e−βt abklingen, wie sich die Dämpfung von schwach über kritisch, bc=2km, bis stark steigert und wie kleine Stöße im Takt eine Resonanz mit der Amplitude QF0/k aufbauen.
Drei gleiche Feder-Masse-Systeme mit schwacher, kritischer und starker Dämpfung werden nach unten gezogen und losgelassen, während eine vierte Feder sanft mit ihrer Eigenfrequenz geschüttelt wird und ihre Schwingungen anwachsen.
Alles, was schwingt oder hüpft, hat eine Eigenfrequenz, also seinen eigenen Rhythmus; ein kürzeres Pendel oder eine steifere Feder schwingt schneller.
Reibung und Luftwiderstand nehmen bei jeder Schwingung etwas Energie und wandeln sie in Wärme um, daher schrumpfen die Schwingungen: das ist Dämpfung.
Schwache Dämpfung ergibt viele kleiner werdende Schwingungen; starke Dämpfung ergibt gar kein Hüpfen, nur ein langsames Zurückkriechen.
Kritische Dämpfung stoppt das Hüpfen gerade eben, sodass das Objekt am schnellsten zur Ruhe kommt, wie eine Tür mit gutem Türschließer.
Sanfte Stöße im Takt der Eigenfrequenz addieren sich zu einer großen Bewegung: das ist Resonanz.
Alles, was schwingen oder hüpfen kann, hat eine Eigenfrequenz: seinen eigenen Rhythmus, also die Anzahl der Hin- und Her-Bewegungen pro Sekunde, wenn man es sich selbst überlässt. Ein kurzes Pendel schwingt schneller als ein langes, und ein Block an einer steifen Feder hüpft schneller als an einer weichen.
Schritt 1: Dämpfung. Stoße eine Schaukel an und geh weg: Jeder Hingang ist ein bisschen kleiner als der vorherige. Luftwiderstand und Reibung nehmen bei jedem Hingang etwas Energie und wandeln sie in Wärme um, und sie bremsen umso stärker, je schneller sich das Objekt bewegt. Dieses allmähliche Abklingen heißt Dämpfung.
Schritt 2: Wie stark ist die Dämpfung? In der Animation hängen drei gleiche Blöcke an gleichen Federn, und alle drei werden nach unten gezogen und gleichzeitig losgelassen. Der linke bewegt sich in Luft (schwache Dämpfung) und hüpft viele Male, jedes Hüpfen kleiner. Der rechte steckt in dickem Sirup (starke Dämpfung): Er schießt nie über die Ruhelage hinaus, kriecht aber sehr langsam zurück. Der mittlere, in einer dünneren Flüssigkeit, hat gerade so viel Dämpfung, dass das Hüpfen aufhört. Das ist die kritische Dämpfung, und ihr Etikett „in Ruhe“ erscheint zuerst. Stoßdämpfer am Auto und Türschließer sind nahe daran eingestellt.
Schritt 3: Resonanz. Kleine Stöße im Takt der Eigenfrequenz addieren im richtigen Moment jeweils Energie, sodass die Bewegung wächst. Rechts wird das obere Ende einer Feder sanft im Rhythmus des Blocks geschüttelt, und bald hüpft der Block viel weiter als das Schütteln. Die Dämpfung ist es, die das Wachsen stoppt.
Teste dich selbst: Eine Waschmaschine schüttelt sich beim Schleudern oft für einige Sekunden stark, während die Trommel schneller wird, und beruhigt sich dann bei voller Drehzahl. Warum? Antwort: Beim Hochdrehen durchläuft die Drehzahl der Trommel die Eigenfrequenz der Maschine, sodass sie für einen Moment im Takt geschüttelt wird: Resonanz.
Mittelstufe
Kernideen
Die Dämpfungskraft F=−bv wirkt der Bewegung entgegen; bei schwacher Dämpfung klingt die Amplitude gemäß A=A0e−βt ab, mit β=b/(2m).
Kritische Dämpfung bei bc=2km: darunter schwingt das System, darüber kriecht es zurück, ohne zu schwingen.
Q=ω0/(2β)=mω0/b: Nach etwa Q Schwingungen ist ein frei klingender Oszillator auf etwa 4 % seiner Amplitude abgeklungen.
Bei Antrieb mit f0 ist die stationäre Amplitude A=F0/(bω0)=QF0/k: Bei schwacher Dämpfung werden kleine Stöße zu großer Bewegung.
Zwei Kräfte: die Feder und der Dämpfer
Ein Block an einer Feder, die gedehnt und losgelassen wird, spürt zwei Kräfte: die rücktreibende Federkraft −kx, die zurück zur Gleichgewichtslage zeigt, und eine Dämpfungskraft, die der Bewegung entgegenwirkt, etwa den Widerstand einer Flüssigkeit oder des Öls im Stoßdämpfer. Bei langsamer Bewegung in einer Flüssigkeit ist die Dämpfungskraft nahezu proportional zur Geschwindigkeit, sodass die Nettokraft F=−kx−bv ist.
Die Dämpfungskonstante b wird in Kilogramm pro Sekunde angegeben (kg/s, dasselbe wie N·s/m). Die beiden Kräfte erreichen ihre Maxima zu verschiedenen Zeitpunkten: die Federkraft an den Umkehrpunkten, wo der Block für einen Augenblick stillsteht, und die Dämpfungskraft, wenn sich der Block am schnellsten bewegt, was bei schwacher Dämpfung nahe der Gleichgewichtslage ist.
In der Animation zeigen die goldenen Pfeile −bv; sie kehren sich jedes Mal um, wenn ein Block umkehrt. Der rechte Block, in dickem Sirup mit b=30kg/s, hat den längsten Pfeil, obwohl er sich am langsamsten bewegt: Bei seiner Höchstgeschwindigkeit von 0,077 m/s hebt die Dämpfungskraft 30⋅0,077≈2,3N genau die Zugkraft der Feder auf.
Schwache Dämpfung: ein exponentielles Abklingen
Bei schwacher Dämpfung behält der Block seinen Eigenrhythmus fast bei, doch jede Schwingung ist um denselben Bruchteil kleiner als die vorherige, weil eine mit der Geschwindigkeit wachsende Dämpfungskraft in jedem Zyklus einen festen Anteil der verbleibenden Energie entzieht. Das ist ein exponentielles Abklingen, A=A0e−βt, mit β=2mb,ω0=mk Die Abklingrate β wird in s⁻¹ angegeben: Nach der Zeit 1/β ist die Amplitude auf 1/e gesunken, etwa 37 %, und sie halbiert sich alle ln2/β. Dabei ist ω0 die Eigenkreisfrequenz. Die vollständige Bewegung ist ein Kosinus innerhalb der schrumpfenden Hüllkurve: x=A0e−βtcos(ωdt+φ),ωd=ω02−β2
Der Dämpfer hält den Block etwas zurück, deshalb liegt die gedämpfte Kreisfrequenz ωd knapp unter ω0. Der Gütefaktor misst, wie schwach die Dämpfung ist: Q=2βω0=bmω0 Bei jedem Zyklus wird die Amplitude mit etwa e−π/Q multipliziert, sodass nach Q Zyklen nur e−π, also etwa 4 %, übrig bleiben. Eine Stimmgabel mit einem Q in den Tausenden klingt nach einem einzigen Anschlag noch sekundenlang nach.
Schwache, kritische und starke Dämpfung
Mit wachsendem b wird das Abklingen schneller, bis bei einem kritischen Wert der Block nicht mehr über die Ruhelage hinausschießt: bc=2km=2mω0 Dort ist β=ω0, und ωd fällt auf null. Der Vergleich von b mit bc ergibt drei Arten von Bewegung:
b<bc, schwache Dämpfung: Der Block schießt immer wieder über die Ruhelage hinaus, mit exponentiell abklingenden Schwingungen.
b=bc, kritische Dämpfung: die schnellstmögliche Rückkehr in die Gleichgewichtslage ohne jedes Überschwingen.
b>bc, starke Dämpfung: ebenfalls kein Überschwingen, aber der Dämpfer wirkt so stark, dass der Block nur kriecht.
Zu viel Dämpfung ist aus demselben Grund langsam wie ein Löffel im Honig: Die Dämpfungskraft wird bald so groß wie die Federkraft, was die Geschwindigkeit auf etwa kx/b begrenzt. Verdreifacht man b, kriecht der Block dreimal so langsam zurück.
Beispiel: die drei Blöcke der Animation
Alle drei Blöcke haben m=1kg und k=25N/m, und jeder wird 10 cm nach unten gezogen und aus der Ruhe losgelassen: ω0=25/1=5rad/s,bc=225⋅1=10kg/s Damit ist f0≈0,80Hz, eine Schwingung alle 1,26 s. Der mittlere Block ist genau kritisch; der rechte hat das Dreifache von bc, der linke ein Zehntel davon. Für den linken Block gilt β=2⋅11=0,5s−1,Q=2⋅0,55=5
Nach 3 s beträgt seine Amplitude 10cm⋅e−0,5⋅3≈2,2cm, also der Wert auf der blauen A-Anzeige, wenn die t-Anzeige, die mit dem Loslassen beginnt, 3 s zeigt. Die Amplitude halbiert sich alle ln2/0,5≈1,4s, jeder Zyklus multipliziert sie mit e−π/5≈0,53, und nach 5 Zyklen, etwa 6,3 s, ist nur noch 4 % übrig. Jeder Block bleibt nach etwa 1,0 s (kritisch), 3,8 s (stark) bzw. 6,4 s (schwach) innerhalb von 4 mm um seine gestrichelte Gleichgewichtslinie.
Resonanz: Stöße im Takt der Bewegung
Treibt man den Oszillator nun mit einer periodischen Kraft F0cosωt an, so schwingt er nach dem Einschwingen stetig mit der Antriebsfrequenz, mit der Amplitude A(ω)=(ω02−ω2)2+(2βω)2F0/m Fern von ω0 ist die Antwort klein; in der Nähe steigt die Amplitude steil an: das ist Resonanz. Bei ω0 heben sich Federkraft und Trägheit auf, und der Antrieb muss nur die Dämpfung ausgleichen: A(ω0)=bω0F0=kQF0 Eine konstante Kraft würde die Feder um F0/k dehnen; im Takt angelegt, bewegt sie den Block Q-mal so weit. In der Resonanz hinkt der Block um eine Viertelperiode hinter dem Stoß her, sodass er am schnellsten ist, wenn der Stoß am stärksten ist, und jeder Stoß verrichtet positive Arbeit.
Rechts in der Animation wird das obere Ende einer gleichen, schwach gedämpften Feder mit Y=2cm genau bei f0=0,80Hz geschüttelt. Verschiebt man das obere Ende der Feder um y, ändert sich ihre Zugkraft um ky; das Schütteln wirkt daher wie eine Kraft F0=kY=25⋅0,02=0,5N, und A=bω0F0=1⋅50,5=0,10m Das ist QY=10cm. Die A-Anzeige des angetriebenen Blocks steigt dort langsam an, etwa wie QY(1−e−βt): etwa 6 cm nach 2 s Schütteln, 9 cm nach 4,6 s und 9,9 cm am Ende der Schleife.
Beispiel: ein Auto mit verschlissenem Stoßdämpfer
Eine Ecke eines Autos trägt etwa 400 kg auf einer Federung mit k=40000N/m: ω0=40000/400=10rad/s,bc=240000⋅400=8000kg/s Damit hüpft die Karosserie mit etwa 1,6 Hz. Ein intakter Stoßdämpfer, etwas unter bc eingestellt, lässt die Karosserie nach einer Bodenwelle einmal überschwingen. Ein verschlissener mit b=800kg/s ergibt β=800/(2⋅400)=1,0s−1 und Q=10/(2⋅1,0)=5, wie der linke Block der Animation. Ein Hüpfer von 4 cm halbiert sich dann alle 0,7 s und ist nach 2 s noch 4⋅e−2≈0,5cm groß, drei Hüpfer später.
Auch die Straße kann das Auto antreiben. Bodenwellen im Abstand d=10m lassen das Rad v/d-mal pro Sekunde wackeln, was bei v≈16m/s, etwa 57 km/h, zur Eigenfrequenz von 1,59 Hz passt. Bei dieser Geschwindigkeit hebt und senkt sich das verschlissene Auto etwa Q=5-mal so stark wie die Bodenwellen und scheint zu schweben.
Häufige Fehler
Mehr Dämpfung ist nicht immer schneller: Der rechte Block, mit dem Dreifachen von bc, braucht etwa 3,8 s, bis er zur Ruhe kommt, der kritische nur 1,0 s.
Schwache Dämpfung ändert die Amplitude, nicht den Rhythmus: Bei Q=5 wächst die Periode nur von 1,257 s auf 1,263 s.
Die Zeit 1/β ist keine Halbwertszeit: Mit β=0,5s−1 fällt die Amplitude in 2 s auf 37 % und halbiert sich in 1,4 s; die Energie, die proportional zu A2 ist, halbiert sich in 0,7 s.
Resonanz ist nicht unbegrenzt: Die Dämpfung begrenzt die Amplitude auf QF0/k, 10 cm bei der geschüttelten Feder, wo die Leistung des Antriebs von 0,125 W gleich den Verlusten des Dämpfers ist.
Der kritische Wert ist 2km, nicht km: Vergisst man die 2, erhält man für die Blöcke 5 kg/s statt 10 kg/s.
Wo man es antrifft
Nahezu kritische Dämpfung steckt in Türschließern, Messzeigern und Autofederungen; nahe der Spitze des Taipei-101-Turms schwingt ein Stahlpendel von etwa 660 Tonnen gegen das Schwanken des Gebäudes. Resonanz wird in Quarzuhren genutzt, deren Kristall mit 32 768 Hz schwingt, in Musikinstrumenten und beim Radioabstimmen. In Bauwerken muss man sie vermeiden: Die Millennium Bridge in London schwankte in ihren Eröffnungstagen im Jahr 2000 seitwärts, weil Fußgänger im Takt mit ihr zu gehen begannen. Dieselbe Gleichung beschreibt Schaltkreise aus den Bauteilen in Kondensatoren und Elektromagnetische Induktion und jede Schwingungsform in Stehende Wellen.
Fortgeschritten
Kernideen
Schreibe x=Re(x^eiωt): Jedes d/dt wird zum Faktor iω, einer Vierteldrehung, sodass die lineare Gleichung zur Algebra wird.
Freie Bewegung: s=−β±β2−ω02; komplexe Wurzeln (ζ<1) klingen, eine doppelte Wurzel (ζ=1) ist kritisch, reelle Wurzeln kriechen.
Q = gespeicherte Energie ÷ Energieverlust pro Radiant: dE/dt=−ω0E/Q, also E∝e−2βt, und die Amplitude klingt halb so schnell ab.
A=∣x^∣ erreicht bei ω0 den Wert QF0/k und hinkt der Kraft um 90° nach; die Breite bei halber Leistung ist Δω≈ω0/Q.
Dieselbe Gleichung gilt anderswo: ein RLC-Schwingkreis (L, R, 1/C für m, b, k; Q=ω0L/R), Gezeiten der Atmosphäre, Ionen, die Infrarot absorbieren.
Die Gleichung der freien Bewegung
Dämpfung ist jede Kraft, die die Energie einer Schwingung abzieht: Luft, Öl, ein Magnet, der eine Metallplatte bremst. Dieser Abschnitt ermittelt alle möglichen freien Bewegungen eines gedämpften Feder-Masse-Systems. Das ist aus zwei Gründen wichtig: Diese Bewegungen zeigen, wie sich ein gestörtes System beruhigt, und sie sind zugleich der Übergangsvorgang, den jeder angetriebene Oszillator abbauen muss, bevor er seinen stationären Zustand erreicht.
Betrachte einen Block der Masse m an einer Feder mit der Federkonstante k, mit einem Dämpfer, der mit der Kraft −bx˙ gegen die Geschwindigkeit wirkt. Dabei ist x die Auslenkung aus der Gleichgewichtslage und b die Dämpfungskonstante. Das zweite newtonsche Gesetz liefert mx¨=−kx−bx˙. Teilt man durch m und führt zwei Verhältnisse ein, die Abklingrate β=b/(2m) und die Eigenkreisfrequenz ω0=k/m, nimmt die Gleichung ihre Standardform an x¨+2βx˙+ω02x=0 Der Faktor 2 ist eine Konvention, die die Wurzeln unten frei von Brüchen hält.
Die Gleichung ist linear, und ihre Koeffizienten ändern sich nicht mit der Zeit. Eine Exponentialfunktion ist die einzige Funktion, deren Ableitung ein Vielfaches von ihr selbst ist; daher probiere x=est mit einer konstanten s, die komplex sein darf. Jede zeitliche Ableitung multipliziert dann nur mit s. Jeder Term trägt den gleichen Faktor est, und Herausdividieren ergibt die charakteristische Gleichung s2+2βs+ω02=0. Die quadratische Formel liefert ihre beiden Wurzeln: s1=−β+β2−ω02,s2=−β−β2−ω02 Da die Gleichung linear ist, ist jede Summe c1es1t+c2es2t ebenfalls eine Lösung, und die beiden Konstanten werden durch Anfangsposition und Anfangsgeschwindigkeit festgelegt. Alles an der freien Bewegung entscheidet sich also daran, wo die beiden Wurzeln liegen.
Die Animation verwendet m=1kg und k=25N/m, also ω0=5rad/s (0,80 Hz), mit drei Dämpfern, b=1, 10 und 30kg/s, das heißt β=0,5, 5 und 15s−1.
Drei Arten freier Bewegung
Ob ein losgelassener Block hin und her schwingt oder nur zurückkriecht, hängt allein vom Vorzeichen von β2−ω02 unter der Wurzel ab. Eine Karte der Wurzeln zeigt alle drei Fälle auf einmal und außerdem, welche Dämpfung den Block am schnellsten zurückbringt.
Abbildung 1. Die komplexe Ebene der Wurzeln s: Realteil waagerecht, Imaginärteil senkrecht. Mit wachsender Dämpfung wandert das Paar −β±iωd auf dem blauen Bogen mit Radius ω0, trifft sich bei −ω0 (kritisch) und teilt sich dann auf der reellen Achse.Legende
s1
Wurzel mit dem Pluszeichen: der obere blaue Punkt, dann auf der reellen Achse, zurück in Richtung 0 wandernd.
s−1
s2
Wurzel mit dem Minuszeichen: der untere blaue Punkt, dann auf der reellen Achse, sich von 0 entfernend.
s−1
β
Abklingrate b/(2m): minus der Realteil des Paares, also sein Abstand links von der senkrechten Achse.
s−1
ωd
Gedämpfte Kreisfrequenz ω02−β2: der Imaginärteil von s1, die Höhe des oberen Punkts.
rad/s
ω0
Eigenkreisfrequenz k/m: der Radius des Bogens; bei −ω0 treffen sich die beiden Wurzeln.
rad/s
Abbildung 1 ist die komplexe Ebene von s: Realteil waagerecht, Imaginärteil senkrecht. Eine einzige Zahl ordnet die Fälle, das Dämpfungsverhältnis ζ=β/ω0. Es ist gleich b/bc, wobei bc=2km die kritische Dämpfungskonstante ist, der Wert von b, der β=ω0 macht. In der Animation ist bc=10kg/s, daher haben die drei Blöcke ζ=0,1, 1 und 3.
Schwache Dämpfung, ζ<1. Die Zahl unter der Wurzel ist negativ, daher sind die Wurzeln das komplexe Paar −β±iωd, die beiden blauen Punkte, mit der gedämpften Kreisfrequenz ωd=ω02−β2. Da e(−β±iωd)t=e−βte±iωdt gilt, legt der Realteil −β das Abklingen fest und der Imaginärteil ±ωd die Schwingung. Fasst man die beiden Exponentialfunktionen so zusammen, dass x reell wird, ergibt sich x=Ce−βtcos(ωdt+φ) eine Schwingung innerhalb einer Hüllkurve, die als e−βt schrumpft; die Konstanten C und φ folgen aus Anfangsposition und Anfangsgeschwindigkeit. Betrachte nun das rechtwinklige Dreieck, gezeichnet für ζ=0,6. Seine waagerechte Kathete ist β, seine senkrechte Kathete ωd, und seine Hypotenuse β2+ωd2=ω0. Eine komplexe Wurzel liegt also immer im Abstand ω0 vom Ursprung. Wenn b wächst, wandert das Paar auf dem blauen Bogen mit dem Radius ω0, die obere Wurzel nach unten und die untere nach oben, auf die reelle Achse zu.
Kritische Dämpfung, ζ=1. Die beiden Wurzeln treffen sich bei −ω0. Eine doppelte Wurzel liefert nur eine Exponentialfunktion, e−ω0t; die zweite unabhängige Lösung ist te−ω0t, wie das Einsetzen in die Gleichung bestätigt.
Starke Dämpfung, ζ>1. Die Quadratwurzel ist nun reell und kleiner als β, daher sind beide Wurzeln reell und negativ. Wenn b wächst, wandert s1 zurück auf 0 zu, und s2 läuft nach links weg. Nichts schwingt: Der Block kriecht zurück.
Der Graph x(t) der Animation im Expertenmodus verwendet das Blau, Grün und Violett aus Abbildung 1 und zeigt die drei Blöcke, die aus der Ruhe bei −10cm losgelassen werden; diese Auslenkung beim Loslassen nennen wir A0=10cm. Blau (ζ=0,1, s=−0,5±4,97is−1) schießt eine halbe Periode nach dem Loslassen, bei 0,63 s, bis +7,3cm über die Ruhelage hinaus. Grün (ζ=1) folgt x=−A0(1+ω0t)e−ω0t und liegt nach etwa 1,0 s innerhalb von 4 mm um null. Violett (ζ=3) hat die Wurzeln −0,86 und −29,1s−1. Sein schneller Term, e−29,1t, ist innerhalb einer Zehntelsekunde verschwunden. Die langsame Wurzel liegt bei starker Dämpfung nahe −ω02/(2β) und hinterlässt ein Kriechen mit einer Zeitkonstante von 1,2 s.
Welche Dämpfung bringt den Block am schnellsten zurück? Eine Summe von Exponentialfunktionen klingt nur so schnell ab wie ihr langsamster Term, also derjenige, dessen Wurzel in Abbildung 1 der senkrechten Achse am nächsten liegt. Unterhalb der kritischen Dämpfung ist dieser Abstand β; er wächst mit b, wenn sich das Paar auf dem Bogen nach links bewegt. Oberhalb der kritischen Dämpfung ist er ∣s1∣; er schrumpft, wenn b wächst, weil s1 zurück gegen 0 driftet. Die beiden Tendenzen treffen sich bei ζ=1: Dort liegt die nächste Wurzel so weit links wie möglich, und die Abklingrate erreicht ihren größten Wert, ω0. Wenn ein kleines Überschwingen akzeptabel ist, geht es mit etwas weniger Dämpfung noch schneller. Bei ζ=0,75 bleibt der Block nach 0,64 s innerhalb von 4 mm um null, weshalb Messzeiger oft etwas unter der kritischen Dämpfung abgestimmt werden.
Energie und die Bedeutung von Q
Die Wurzeln sagen, wie schnell die Bewegung abklingt; die Energie sagt, warum. Sie liefert außerdem die eine Zahl, die man für jeden Oszillator angibt: den Gütefaktor Q.
Abbildung 2. Ein leicht gedämpfter Block, aus der Ruhe bei −A0 losgelassen: Position x oben, Energie E unten. E nimmt mit bx˙2 ab, daher ist es an jedem Umkehrpunkt flach, wo es die glatte Abnahme E0e−2βt berührt, und am steilsten, wo der Block am schnellsten ist.Legende
x
Auslenkung des Blocks aus der Gleichgewichtslage.
m
t
Zeit seit dem Loslassen.
s
E
Gespeicherte Energie, kinetisch plus Feder: 21mx˙2+21kx2.
J
E0
Energie beim Loslassen, ganz in der Feder: 21kA02.
J
e−2βt
Abklingfaktor der gestrichelten Kurve E0e−2βt, die die Energie an jedem Umkehrpunkt berührt.
dE/dt
Änderungsrate der Energie, −bx˙2: am größten dort, wo der Block sich am schnellsten bewegt, und null an einem Umkehrpunkt.
W
Td
Gedämpfte Periode 2π/ωd: ein voller Zyklus, die grüne Spanne.
s
Abbildung 2 verfolgt den leichten Block der Animation (ζ=0,1), losgelassen aus der Ruhe bei x=−A0: seine Position x(t) oben und seine gespeicherte Energie E unten, auf einer gemeinsamen Zeitachse. Um das untere Feld zu lesen, ermittle zunächst, wie schnell sich E ändert. Multipliziere die Bewegungsgleichung mx¨+bx˙+kx=0 mit x˙. Das Produkt mx¨x˙ ist die zeitliche Ableitung von 21mx˙2, und kxx˙ ist die von 21kx2. Die Trägheits- und Federterme zusammen ergeben also die Ableitung von E=21mx˙2+21kx2, und was übrig bleibt, lautet dtdE=−bx˙2 Der Dämpfer entzieht Energie mit der Rate bx˙2. Diese Rate ist nie negativ und null, sobald der Block für einen Augenblick stillsteht.
Sieh dir nun das untere Feld an. Energie wird nur verloren, solange sich der Block bewegt, daher fällt E stufenweise. Jeder Punkt markiert einen Umkehrpunkt, an dem x˙=0 gilt und die gesamte Energie in der Feder steckt; dort ist die Steigung von E null, und die Kurve ist flach. Die mit dE/dt markierte Tangente berührt den steilsten Teil des ersten Abfalls, dort wo sich der Block am schnellsten bewegt, kurz bevor er die Gleichgewichtslage passiert. An einem Umkehrpunkt gilt ∣x∣=A0e−βt, also ist die Energie dort 21kA02e−2βt. Die flachen Stellen liegen daher genau auf der gestrichelten Kurve E0e−2βt, mit E0=21kA02=21⋅25⋅0,102=0,125J. Zwischen den Umkehrpunkten schwankt E etwa 10 % über und unter dieser Kurve. Von einem Umkehrpunkt zum nächsten, eine halbe gedämpfte Periode Td/2, wird die Energie mit e−βTd≈0,53 multipliziert.
Glättet man die Stufen, ergibt sich ein einfaches Gesetz. Über einen Zyklus bei schwacher Dämpfung ist die Bewegung nahezu Acosω0t mit einer langsam schrumpfenden Amplitude A. Der Mittelwert von x˙2 über den Zyklus ist dann 21ω02A2, während E=21mω02A2 gilt. Der Mittelwert von bx˙2 ist also (b/m)E=2βE, und dtdE=−2βE Die Energie klingt wie e−2βt ab. Die Energie ist proportional zum Quadrat der Amplitude, daher klingt die Amplitude halb so schnell ab, nämlich wie e−βt, genau wie die Wurzeln es vorhergesagt haben.
Der Gütefaktor ist definiert als Q=ω0/(2β). Damit lautet das Gesetz dE/dt=−ω0E/Q. Der Block durchläuft einen Radiant Phase in der Zeit 1/ω0 und verliert in dieser Zeit den Bruchteil 1/Q seiner Energie: Q ist die gespeicherte Energie geteilt durch den Energieverlust pro Radiant. Ein Oszillator mit hohem Q klingt über viele Zyklen nach; der leichte Block hat nur Q=5. Eine bekannte Faustregel besagt, dass pro Zyklus der Bruchteil 2π/Q verloren geht. Hier ergäbe das 126 %, was unmöglich ist: Die Regel ist nur das erste Glied des exakten Ausdrucks 1−e−2π/Q und nur für großes Q brauchbar. Das exponentielle Gesetz selbst gilt von einem Zyklus zum nächsten exakt. Über eine gedämpfte Periode Td=2π/ωd=1,263s, die grüne Spanne in Abbildung 2, behält die Energie den Anteil e−2βTd≈0,28, also gehen 72 % im ersten Zyklus verloren.
Der getriebene Oszillator
Treibe den Block nun mit einer periodischen Kraft F0cosωt an. Das ist die Situation einer Lautsprechermembran, einer Brücke im böigen Wind oder eines Elektrons in einer Lichtwelle, und die Fragen lauten, wie groß die Antwort ist und wie weit sie hinterherhinkt. Die Bewegungsgleichung erhält eine rechte Seite: x¨+2βx˙+ω02x=(F0/m)cosωt Da die Gleichung linear ist, ist ihre allgemeine Lösung eine beliebige partikuläre Lösung plus die oben bestimmte freie Bewegung. Der freie Anteil klingt ab, bei schwacher Dämpfung wie e−βt; er heißt Übergangsvorgang. Was übrig bleibt, ist der stationäre Zustand, der mit der Antriebsfrequenz ω schwingt, nicht mit ω0.
Komplexe Zahlen finden den stationären Zustand schnell. Schreibe die Kraft als Realteil von F0eiωt und suche x als Realteil von x^eiωt, wobei die komplexe Amplitude x^ eine noch zu bestimmende Konstante ist. Die Bildung des Realteils vertauscht mit Addition, mit Ableiten und mit Multiplikation mit reellen Zahlen, daher kann man die komplexe Gleichung lösen und am Ende den Realteil nehmen. Jedes d/dt, das auf eiωt wirkt, wird zum Faktor iω. Multiplikation mit i ist eine Vierteldrehung in der komplexen Ebene, daher eilt die Geschwindigkeit der Auslenkung um 90° voraus. Die Differentialgleichung wird zur Algebra: (ω02−ω2+2iβω)x^=F0/m Der Trick benötigt Linearität: Ein Term in x2 würde Real- und Imaginärteil vermischen, und der Realteil der komplexen Lösung würde die reelle Gleichung nicht mehr erfüllen. Linearität liefert außerdem die Superposition, sodass eine Kraft aus vielen Frequenzen die Summe der einzelnen Antworten erzeugt.
Das Zeigerdreieck: Amplitude und Phasenverschiebung
Die komplexe Gleichung sagt, dass zwei Pfeile einen dritten ergeben. Zeichnet man sie, werden Amplitude und Phasenverschiebung zu einfacher Geometrie.
Abbildung 3. Die Gleichung (ω02−ω2+2iβω)x^=F0/m als Pfeile in der komplexen Ebene, unterhalb der Resonanz. Grüner und oranger Term stehen im rechten Winkel und ergeben zusammen den Antrieb F0/m als Hypotenuse; δ ist die Verzögerung von x^ hinter dem Antrieb.Legende
F0/m
Der Antrieb, Kraft pro Masse, auf der reellen Achse.
m/s2
x^
Komplexe Amplitude der Auslenkung, im Winkel −δ zum Antrieb; im eigenen Maßstab gezeichnet.
m
(ω02−ω2)x^
Feder- und Trägheitsterm zusammen: ein reelles Vielfaches von x^, zeigt also in Richtung von x^.
m/s2
(2βω)ix^
Dämpfungsterm: x^, um eine Vierteldrehung nach vorn gedreht durch den Faktor i und mit 2βω skaliert.
m/s2
δ
Verzögerung der Auslenkung hinter dem Antrieb: der Winkel zwischen x^ und F0/m.
rad
Abbildung 3 zeichnet die Gleichung in der komplexen Ebene für ω=0,75ω0 und β=0,25ω0. Der Antrieb F0/m liegt auf der reellen Achse. Der goldene Pfeil ist x^, im eigenen Maßstab gezeichnet. Der Federterm ω02x^ und der Trägheitsterm −ω2x^ sind beide reelle Vielfache von x^, daher ergeben sie zusammen den grünen Pfeil (ω02−ω2)x^, der entlang x^ zeigt. Der Dämpfungsterm 2βωix^ trägt den zusätzlichen Faktor i, daher steht der orangefarbene Pfeil im rechten Winkel zum grünen. Grün plus Orange muss den Antrieb ergeben. Der Antrieb ist also die Hypotenuse eines rechtwinkligen Dreiecks, und x^, das entlang der grünen Kathete liegt, hinkt dem Antrieb um den Winkel δ zwischen dieser Kathete und der Hypotenuse hinterher.
Die Gleichung nach x^ aufgelöst ergibt x^=ω02−ω2+2iβωF0/m Sein Betrag ist die Amplitude A, und sein Argument ist −δ, wobei δ die Phasenverschiebung der Auslenkung hinter der Kraft ist. Das Dreieck liefert beides ohne weitere Rechnung. Seine Katheten haben die Längen (ω02−ω2)A und 2βωA, und seine Hypotenuse ist F0/m. Pythagoras liefert die Amplitude und das Verhältnis der Katheten die Phasenverschiebung: A(ω)=(ω02−ω2)2+(2βω)2F0/m tanδ=ω02−ω22βω Wenn ω über ω0 hinaus ansteigt, steigt δ von 0 auf 180∘. Abbildung 3 ist unterhalb der Resonanz gezeichnet, wo die grüne Kathete positiv ist und δ unter 90∘ liegt. Oberhalb von ω0 wird der Faktor ω02−ω2 negativ, der grüne Pfeil zeigt entgegen x^, und δ überschreitet 90∘.
Die Resonanzkurve lesen
Die beiden Formeln oben enthalten die gesamte stationäre Antwort. Trägt man sie gegen die Antriebsfrequenz auf, sieht man, wo der Oszillator stark reagiert, wie scharf diese Reaktion ist und was die Phase auf dem Weg dorthin tut.
Abbildung 4. Die stationäre Antwort eines Blocks mit Q=3 auf einen Antrieb der Kreisfrequenz ω: Amplitude A oben, Phasenverschiebung δ unten. Bei ω0, knapp rechts vom Maximum, ist A=QF0/k und δ=π/2; die Breite zwischen den Halbleistungsfrequenzen ω1 und ω2 ist ω0/Q.Legende
A
Stationäre Amplitude der Auslenkung.
m
δ
Verzögerung der Auslenkung hinter dem Antrieb.
rad
ω
Kreisfrequenz des Antriebs.
rad/s
ω0
Eigenkreisfrequenz: knapp rechts vom Amplitudenmaximum, wo δ=π/2 gilt.
rad/s
ω1
Untere Halbleistungsfrequenz, wo δ=π/4 gilt.
rad/s
ω2
Obere Halbleistungsfrequenz, wo δ=3π/4 gilt.
rad/s
Δω
Volle Breite bei halber Leistung: ω2−ω1=2β=ω0/Q.
rad/s
F0/k
Statische Dehnung: die Amplitude bei einem sich sehr langsam ändernden Antrieb.
m
QF0/k
Amplitude bei ω0: Q-mal die statische Dehnung.
m
π/2
Eine Vierteldrehung: die Verzögerung bei ω0.
rad
π
Eine halbe Drehung: die Verzögerung weit oberhalb von ω0, wo der Block sich gegen den Antrieb bewegt.
rad
Abbildung 4 skizziert die stationäre Antwort eines Blocks mit Q=3: die Amplitude A(ω) oben und die Phasenverschiebung δ(ω) unten, auf derselben Frequenzachse. Das Resonanzfeld der Animation zeigt A(ω) für die geschüttelte Feder, mit Q=5 und Q=1. Dort wird das obere Ende der Feder um Y=2cm bewegt, was auf den Block wie eine Kraft mit F0/m=ω02Y wirkt, also F0/k=Y.
Lies die Kurven von links nach rechts. Weit unterhalb von ω0 dominiert der Federterm: A≈F0/k, die statische Dehnung (2 cm in der Animation), δ≈0, und der Block folgt dem Schüttler. Weit oberhalb von ω0 dominiert die Trägheit: A≈F0/(mω2) schrumpft unter das gestrichelte Niveau F0/k, δ nähert sich π, und der Block bewegt sich gegen den Antrieb. Bei ω0, knapp rechts vom Maximum, heben sich Feder und Trägheit auf, und nur die Dämpfung hält dem Antrieb das Gleichgewicht: A=QF0/k und δ=π/2. Die Geschwindigkeit ist dann in Phase mit der Kraft, daher verrichtet der Antrieb während des ganzen Zyklus positive Arbeit. In der Animation erreicht die hohe Q=5-Kurve 10 cm, während die niedrige graue Q=1-Kurve bei ω0 nur 2 cm ergibt.
Die Breite des Maximums misst man mit der Leistung. Im stationären Zustand ist die mittlere Leistung, die dem Antrieb entnommen wird, gleich der mittleren Leistung, die der Dämpfer umsetzt, P=21bω2A2, was für Q=5 bei Resonanz 0,125 W beträgt. Schreibt man ω2A2=(F0/m)2/[(ω02−ω2)2/ω2+4β2], sieht man, dass P genau bei ω0 sein Maximum hat und auf die Hälfte fällt, wo ∣ω02−ω2∣=2βω gilt. Das sind die Frequenzen ω1 und ω2 der Abbildung; für Q=5 betragen sie 4,52 und 5,52 rad/s (0,72 und 0,88 Hz). Dort sind die beiden Katheten des Dreiecks aus Abbildung 3 gleich lang, daher gilt tanδ=±1: Die Phasenkurve geht bei ω1 durch π/4 und bei ω2 durch 3π/4. Die volle Breite ist Δω=ω2−ω1=2β=Qω0 sodass Q auch die Schärfe des Maximums misst, ω0/Δω. Das Amplitudenmaximum selbst liegt etwas unter ω0, bei ω02−2β2≈4,95rad/s, wo es 10,05 cm statt 10 cm erreicht.
Der Punkt auf der Kurve der Animation zeigt den Einschwingvorgang. Der Schüttler beginnt in seiner Mitte und bewegt sich als Ysinω0t, während der Block aus der Ruhe startet. Die stationäre Auslenkung hinkt dem Antrieb um eine Vierteldrehung hinterher und ist damit −QYcosω0t; anfangs hebt der Übergangsvorgang sie teilweise auf, sodass x≈−QY(1−e−βt)cosω0t gilt. Die Amplitude wächst auf der Zeitskala 1/β=2s und erreicht erst nach ln10/β≈4,6s 90 % des Maximums. Ein schmales Maximum und ein langes Ausklingen sind ein und dasselbe, zweimal betrachtet: Die Breite 2β mal die Abklingzeit 1/β ist immer 2, daher reagiert ein scharf abgestimmter Oszillator träge.
Durchgerechnete Aufgabe: Schütteln im falschen Takt
Schüttle dieselbe Feder (F0/m=0,5m/s2, β=0,5s−1, ω02=25s−2, b=1kg/s) mit 0,50 Hz und mit 1,00 Hz, jeweils unterhalb und oberhalb ihrer Eigenfrequenz von 0,80 Hz. Bestimme die stationäre Amplitude, die Phasenverschiebung und die mittlere Leistung in jedem Fall.
Unterhalb der Resonanz, bei 0,50 Hz: ω=2πf=π≈3,14rad/s, also ω02−ω2=25−9,87=15,13s−2 und 2βω=3,14s−2. Dann A=15,132+3,1420,5=0,032m,tanδ=15,133,14,δ≈12∘ Oberhalb der Resonanz, bei 1,00 Hz: ω=2π≈6,28rad/s, also ω02−ω2=25−39,48=−14,48s−2 und 2βω=6,28s−2. Der negative Nenner legt δ zwischen 90∘ und 180∘: A=14,482+6,2820,5=0,032m,tanδ=−14,486,28,δ≈157∘
Beide Amplituden betragen etwa 3,2 cm, ein Drittel der Resonanzamplitude von 10 cm. Der Unterschied liegt in der Phase: Bei 0,50 Hz fährt der Block mit dem Schüttler mit, bei 1,00 Hz bewegt er sich fast gegen ihn. Auch die mittlere Leistung 21bω2A2 unterscheidet sich: etwa 5 mW bei 0,50 Hz und 20 mW bei 1,00 Hz, gegenüber 125 mW bei Resonanz. Die Amplituden sind gleich, aber die Frequenz hat sich verdoppelt, daher ist die Geschwindigkeitsamplitude ωA bei 1,00 Hz doppelt so groß und die Leistung viermal so groß.
Auch der Start zählt. Aus der Ruhe gestartet, überlagert sich das freie Schwingen mit ωd=4,97rad/s der stationären Bewegung und schwebt mit der Periode 2π/∣ω−ωd∣, hier 3,4 s und 4,8 s. Startet der Schüttler, wie in der Animation, aus seiner Mitte, so erreichen die frühen Schwingungen 4,1 bzw. 4,5 cm, mehr als die stationären 3,2 cm. Nach 3/β=6s ist der Übergangsvorgang auf e−3 gesunken, also etwa 5 %.
Über lineare Dämpfung hinaus
Bisher nahmen wir an, dass die Dämpfungskraft proportional zur Geschwindigkeit ist. Reibung ist oft nicht so, und die Form des Abklingens verrät, welches Gesetz wirkt.
Trockene Reibung ist eine konstante Kraft μN gegen die Bewegung, mit dem Reibungskoeffizienten μ und der Normalkraft N. Sie lässt die Periode 2π/ω0 unverändert, verringert aber die Amplitude bei jeder halben Schwingung um denselben Betrag 2μN/k. Der Grund: Über eine halbe Schwingung von A1 nach A2 ist die verlorene Federenergie, 21k(A12−A22), gleich der Arbeit μN(A1+A2) gegen die Reibung. Die Hüllkurve ist eine Gerade, und der Block bleibt nach endlicher Zeit stehen, sobald die Federkraft die Reibung nicht mehr überwinden kann. Bei 0,1 N Reibung würde der Block der Animation bei jeder halben Schwingung 8 mm verlieren und nach 12 davon stehen bleiben, also nach etwa 7,5 s. Ein mit v2 wachsender Widerstand dämpft große Schwingungen schnell und lässt kleine länger bestehen, sodass die Amplitude etwa wie 1/(1+t/t1) abnimmt, mit einer Zeit t1, die von der Anfangsamplitude abhängt. Ein praktischer Test: Trage den Logarithmus der Amplitude gegen die Zeit auf. Nur lineare Dämpfung ergibt eine Gerade.
Auch die Feder kann nichtlinear sein. Eine sich versteifende Feder schwingt bei großer Amplitude schneller, daher neigt sich ihre Resonanzkurve zu hohen Frequenzen, und beim Durchstimmen der Antriebsfrequenz kann die Amplitude zwischen zwei stationären Werten springen. Energie kann auch ohne periodische Kraft hineingelangen. Eine Schaukel anzutreiben, indem man pro Zyklus zweimal aufsteht und in die Hocke geht, verändert eine Eigenschaft des Oszillators, seine effektive Länge, statt ihn zu stoßen: Das ist parametrische Resonanz. Eine angestrichene Saite oder die Hemmung einer Uhr führt Energie im Takt der Bewegung selbst zu. Das wirkt wie negative Dämpfung, die die Schwingung anwachsen lässt, bis Nichtlinearität sie begrenzt. Die Tacoma-Narrows-Brücke, 1940 zerstört, versagte durch solches selbsterregtes Flattern bei stetigem Wind, nicht durch einfache Resonanz mit einer periodischen Kraft.
Grenzen der Gültigkeit
Der lineare Dämpfer ist exakt für einen elektrischen Widerstand, nahezu exakt für eine Wirbelstrombremse und gut für langsame Bewegung durch eine zähe Flüssigkeit (Strömungswiderstand und Endgeschwindigkeit). Eine echte Feder ist nur bei kleinen Dehnungen linear (Elastizität: hookesches Gesetz, Spannung und Dehnung). Auch die Resonanzamplitude braucht Zeit zum Aufbau: Der Antrieb muss viele Male 1/β=2Q/ω0 andauern, was bei einer Stimmgabel mit Q≈2000 bei 440 Hz 1,4 s sind.
Dämpfung kommt nie allein. Die Molekülstöße, die die Energie abziehen, stoßen den Block auch zufällig an, daher behält ein Oszillator im Gleichgewicht bei der Temperatur T eine mittlere Energie kBT, wobei kB die Boltzmann-Konstante ist. Die Hälfte davon steckt in der Feder, daher ist der Mittelwert von 21kx2 gleich 21kBT, und der Block zittert um kBT/k (quadratischer Mittelwert). Für die Feder mit 25 N/m bei 300 K sind das 1,3×10−11m, weit unterhalb alles Sichtbaren. Für den Cantilever eines Rasterkraftmikroskops mit k=0,1N/m sind es 0,2 nm, eine Untergrenze für das, was das Instrument messen kann. Dieser Zusammenhang zwischen Dissipation und zufälliger Kraft heißt Fluktuations-Dissipations-Theorem. Quanteneffekte treten in Stufen von ℏω0≈5×10−34J auf; gegenüber den 0,125 J des Blocks sind sie völlig vernachlässigbar, und sie zählen nur für Moleküle, Atome und sehr kalte Mikroresonatoren.
Wohin das führt
Wenn man q anstelle von x setzt, beschreibt dieselbe Gleichung einen Reihenschwingkreis, Lq¨+Rq˙+q/C=V0cosωt. Die Induktivität L, der Widerstand R und 1/C übernehmen die Rollen von m, b und k, daher ω0=1/LC und Q=ω0L/R. Ein Radio, das Mittelwellensender im Abstand von 10 kHz bei etwa 1 MHz trennt, braucht eine Breite ω0/Q, die nicht größer als der Abstand ist, also ein Q in der Größenordnung von 100; siehe Kondensatoren, Elektromagnetische Induktion und Elektromagnetische Wellen.
Licht treibt die in Atomen gebundenen Elektronen auf dieselbe Weise an. Unterhalb ihrer ultravioletten Resonanzen wächst die Antwort mit der Frequenz, weshalb Glas Violett stärker bricht als Rot (Prismen, Dispersion und Farbe); bei einer Resonanz wird Energie absorbiert, wie von den Ionen eines Salzkristalls im Infraroten. Dieselbe Kurve, jedes Mal mit eigenem ω0 und Q, passt auf Gezeiten der Atmosphäre, die Schwingungsformen in Stehende Wellen und die Kerne in einem MRT-Gerät. In der Quantenphysik kehrt die Breite mal Abklingzeit als Zusammenhang zwischen der Breite einer Spektrallinie und der Lebensdauer des Zustands zurück, der sie aussendet.
Wichtige Formeln
Nettokraft auf einen gedämpften Oszillator
F=−kx−bv
Die Feder zieht den Block mit −kx zurück zur Gleichgewichtslage, und der Dämpfer widersetzt sich seiner Bewegung mit −bv. Die beiden Kräfte erreichen ihre Maxima zu verschiedenen Zeitpunkten: die Federkraft an den Umkehrpunkten, die Dämpfungskraft, wenn sich der Block am schnellsten bewegt, was bei schwacher Dämpfung nahe der Gleichgewichtslage ist. Beide sind linear in der Bewegung, was die Gleichung exakt lösbar macht.
Symbole
F
Nettokraft auf den Block
N
k
Federkonstante
N/m
x
Auslenkung aus der Gleichgewichtslage
m
b
Dämpfungskonstante
kg/s
v
Geschwindigkeit des Blocks
m/s
Gilt, wenn: Setzt das hookesche Gesetz für die Feder und eine zur Geschwindigkeit proportionale Dämpfungskraft voraus, was für langsame Bewegung in einer Flüssigkeit, Ölstoßdämpfer und Wirbelstrombremsen gilt. Trockene Reibung (konstante Größe) und schneller Luftwiderstand (etwa v2) passen nicht zu dieser Form.
Beispiel: Bei seiner Höchstgeschwindigkeit von 0,077 m/s spürt der schwere Block (b=30kg/s) eine Dämpfungskraft von 30⋅0,077≈2,3N, die genau die Zugkraft der Feder aufhebt, während er noch 9,2 cm von der Gleichgewichtslage entfernt ist.
Abklingrate und Eigenkreisfrequenz
β=2mb,ω0=mk
Diese beiden Raten enthalten die gesamte Physik der freien Bewegung. ω0 sagt, wie schnell der ungedämpfte Oszillator schwingen würde: Steifigkeit gegen Trägheit; β sagt, wie schnell ihm der Dämpfer Energie entzieht: Dämpfung gegen Trägheit. Ihr Verhältnis ζ=β/ω0 entscheidet, ob das System überhaupt schwingt.
Symbole
β
Abklingrate
s−1
b
Dämpfungskonstante
kg/s
m
Masse des Blocks
kg
ω0
Eigenkreisfrequenz
rad/s
k
Federkonstante
N/m
Gilt, wenn: Für jede lineare Feder und jeden linearen Dämpfer. ω0 ist nur ohne Dämpfung die tatsächliche Schwingungsfrequenz; ein gedämpftes System schwingt mit ω02−β2.
Beispiel: Die Blöcke der Animation: ω0=25/1=5rad/s, und β beträgt 0,5, 5 und 15 s⁻¹ für Dämpfungskonstanten von 1, 10 und 30 kg/s.
Leicht gedämpfte freie Schwingung
x=A0e−βtcos(ωdt+φ),ωd=ω02−β2
Unterhalb der kritischen Dämpfung schwingt der Block trotzdem, aber innerhalb einer Hüllkurve A0e−βt, die bei jedem Zyklus um denselben Bruchteil schrumpft. Die Schwingung läuft mit ωd, etwas langsamer als ω0, weil der Dämpfer den Block zurückhält. Die Konstanten A0 und φ sind durch Anfangsposition und Anfangsgeschwindigkeit festgelegt.
Symbole
x
Auslenkung aus der Gleichgewichtslage zur Zeit t
m
A0
Amplitude der Hüllkurve bei t=0
m
β
Abklingrate
s−1
t
Zeit seit dem Loslassen
s
ωd
gedämpfte Kreisfrequenz
rad/s
φ
Phasenkonstante
rad
ω0
Eigenkreisfrequenz
rad/s
Gilt, wenn: Nur für schwache Dämpfung, β<ω0. Bei β=ω0 fällt die Frequenz ωd auf null, und die Bewegung wird zu einer Rückkehr ohne Schwingung.
Beispiel: Der leichte Block: ωd=25−0,25≈4,97rad/s, und nach 3 s beträgt die Hüllkurve 10cm⋅e−1,5≈2,2cm.
Konstante der kritischen Dämpfung
bc=2km=2mω0
Kritische Dämpfung ist erreicht, wenn die Abklingrate b/(2m) gleich der Eigenkreisfrequenz k/m ist, sodass die gedämpfte Frequenz auf null fällt. Nach b aufgelöst ergibt sich 2km. Bei diesem Wert kehrt ein ausgelenktes System so schnell wie möglich in die Gleichgewichtslage zurück, ohne darüber hinauszuschießen.
Symbole
bc
kritische Dämpfungskonstante
kg/s
k
Federkonstante
N/m
m
Masse
kg
ω0
Eigenkreisfrequenz
rad/s
Gilt, wenn: Für eine lineare Feder und einen linearen Dämpfer. Unterhalb von bc schwingt das System, darüber kriecht es zurück; wenn ein kleines Überschwingen akzeptabel ist, erreicht etwas weniger Dämpfung (ein Dämpfungsverhältnis nahe 0,75) sogar noch schneller einen Wert innerhalb weniger Prozent.
Beispiel:m=1kg und k=25N/m: bc=225⋅1=10kg/s, die Dämpfungskonstante des mittleren Blocks.
Gütefaktor
Q=2βω0=bmω0
Q vergleicht die natürliche Rate ω0 mit der Energieabklingrate 2β: Es ist die gespeicherte Energie geteilt durch den Energieverlust pro Radiant Bewegung. Ein hohes Q bedeutet schwache Dämpfung, ein langes Nachklingen nach einem Anschlag und einen hohen, schmalen Resonanzgipfel. Ein frei klingender Oszillator behält von einem Zyklus zum nächsten etwa e−π/Q seiner Amplitude.
Symbole
Q
Gütefaktor
ω0
Eigenkreisfrequenz
rad/s
β
Abklingrate
s−1
m
Masse
kg
b
Dämpfungskonstante
kg/s
Gilt, wenn: Nützlich bei schwacher Dämpfung. Faustregeln wie ein Anteil 2π/Q der pro Zyklus verlorenen Energie setzen ein großes Q voraus; kritische Dämpfung entspricht Q=1/2.
Beispiel: Der leichte Block: Q=5/(2⋅0,5)=5, sodass nach 5 Zyklen nur e−π, etwa 4 %, der Amplitude übrig ist.
Stationäre Amplitude eines getriebenen Oszillators
A(ω)=(ω02−ω2)2+(2βω)2F0/m
Sobald der Übergangsvorgang abgeklungen ist, schwingt ein gedämpfter, durch F0cosωt angetriebener Oszillator mit der Frequenz des Antriebs und dieser Amplitude. Der Term ω02−ω2 ist die Abweichung zwischen Feder und Trägheit, und 2βω ist die Dämpfung. Weit unterhalb von ω0 strebt die Amplitude gegen die statische Dehnung F0/k, weit oberhalb fällt sie wie 1/ω2, und nahe ω0 begrenzt nur der Dämpfungsterm sie.
Symbole
A(ω)
stationäre Amplitude bei der Antriebsfrequenz ω
m
ω
Antriebskreisfrequenz
rad/s
F0
Amplitude der Antriebskraft
N
m
Masse
kg
ω0
Eigenkreisfrequenz
rad/s
β
Abklingrate
s−1
Gilt, wenn: Nur im stationären Zustand, einige Vielfache von 1/β nach Beginn des Antriebs, und nur für eine lineare Feder und einen linearen Dämpfer mit sinusförmiger Kraft.
Beispiel: Die geschüttelte Feder (F0/m=0,5m/s2) bei 1,00 Hz angetrieben: A=0,5/14,482+6,282≈0,032m, ein Drittel der 10 cm bei Resonanz.
Amplitude bei Resonanz
A(ω0)=bω0F0=kQF0
Bei ω=ω0 heben sich Federkraft und Trägheit gegenseitig auf, sodass der Antrieb nur die Dämpfungskraft ausgleichen muss: F0=bω0A. Die Amplitude ist dann Q-mal die statische Dehnung F0/k. Schwache Dämpfung macht aus einer kleinen Kraft eine große Bewegung, und die Auslenkung hinkt der Kraft genau um eine Viertelperiode hinterher.
Symbole
A(ω0)
stationäre Amplitude bei Antrieb mit ω0
m
F0
Amplitude der Antriebskraft
N
b
Dämpfungskonstante
kg/s
ω0
Eigenkreisfrequenz
rad/s
Q
Gütefaktor
k
Federkonstante
N/m
Gilt, wenn: Für den stationären Zustand eines linearen Oszillators. Streng genommen liegt das Amplitudenmaximum etwas unter ω0, bei ω02−2β2, aber der Unterschied ist vernachlässigbar, außer wenn Q klein ist; sehr große Amplituden machen reale Federn nichtlinear.
Beispiel: Die geschüttelte Feder der Animation: F0=kY=0,5N, also A=0,5/(1⋅5)=0,10m, was QY=10cm entspricht.
Glossar
Eigenfrequenz
Die Frequenz, mit der ein System schwingt, wenn es ausgelenkt und sich selbst überlassen wird, ohne Dämpfung: ω0=k/m für einen Block an einer Feder, oder f0=ω0/(2π) in Hertz. Für die Blöcke der Animation beträgt sie 5 rad/s oder 0,80 Hz.
Dämpfung
Der allmähliche Verlust der Energie einer Schwingung, meist als Wärme, durch Reibung, Strömungswiderstand oder elektrischen Widerstand. Schwache Dämpfung verkleinert die Amplitude, lässt den Rhythmus aber fast unverändert.
Viskose Dämpfung
Dämpfung durch eine Kraft, die proportional und entgegengesetzt zur Geschwindigkeit ist, F=−bv, wie bei langsamer Bewegung durch eine Flüssigkeit oder im Stoßdämpfer eines Autos. Sie ergibt ein exponentielles Abklingen; trockene Reibung, eine Kraft konstanter Größe, ergibt dagegen ein lineares.
Schwache Dämpfung
Dämpfung unterhalb des kritischen Werts, b<bc, auch Unterdämpfung genannt. Das System schießt immer wieder über die Ruhelage hinaus, und seine Schwingungen schrumpfen bei jedem Zyklus um denselben Bruchteil.
Kritische Dämpfung
Die Dämpfung bc=2km, die gerade ausreicht, um jedes Überschwingen zu verhindern. Ein ausgelenktes System kehrt dann in kürzestmöglicher Zeit ohne Überqueren in die Gleichgewichtslage zurück, was Türschließer anstreben.
Starke Dämpfung
Dämpfung oberhalb des kritischen Werts, b>bc, auch Überdämpfung genannt. Es gibt keine Schwingung, nur ein langsames Zurückkriechen in die Gleichgewichtslage, und je größer b, desto langsamer das Kriechen.
Getriebene Schwingung
Die Bewegung eines Oszillators, der durch eine periodische äußere Kraft wie F0cosωt angestoßen wird, auch erzwungene Schwingung genannt. Nach dem Abklingen des Einschwingvorgangs schwingt er mit der Antriebsfrequenz ω, nicht mit seiner eigenen.
Übergangsvorgang
Der Teil der Bewegung eines getriebenen Oszillators, der von seinem eigenen freien Nachschwingen stammt. Er klingt wie e−βt ab und ist nach einigen Vielfachen von 1/β verschwunden, sodass der stationäre Zustand bleibt.
Stationärer Zustand
Die dauerhafte Bewegung eines getriebenen Oszillators, nachdem der Übergangsvorgang abgeklungen ist: eine Schwingung mit der Antriebsfrequenz, mit konstanter Amplitude A(ω) und konstanter Phasenverschiebung gegenüber der Kraft.
Resonanz
Das starke Ansprechen eines Oszillators, der nahe seiner Eigenfrequenz angetrieben wird, wo jeder Stoß Energie hinzufügt. Bei einer Dämpfungskraft −bv beträgt die stationäre Amplitude bei Resonanz QF0/k, also Q-mal die Dehnung, die dieselbe Kraft als konstante Kraft bewirken würde.
Resonanzkurve
Das Diagramm der stationären Amplitude oder der aufgenommenen Leistung gegen die Antriebsfrequenz. Sie erreicht ihr Maximum nahe ω0, mit einer Höhe proportional zu Q und einer Breite bei halber Leistung von ω0/Q.
Phasenverschiebung
Wie weit die Auslenkung der antreibenden Kraft nachhinkt, angegeben als Winkel δ, wobei ein ganzer Zyklus als 360∘ gezählt wird. Sie ist weit unterhalb der Resonanz nahe 0, bei ω0 genau 90∘ und weit oberhalb nahe 180∘.
Größen
b
Dämpfungskonstante
kg/s
Die Dämpfungskraft pro Geschwindigkeit; 1 kg/s entspricht 1 N·s/m. Die drei Blöcke der Animation haben 1, 10 und 30 kg/s.
β
Abklingrate
s−1
Die halbe Dämpfungskonstante pro Masseeinheit, b/(2m): In der Zeit 1/β sinkt die Amplitude um den Faktor e. Für den leichten Block beträgt sie 0,5s−1.
ω0
Eigenkreisfrequenz
rad/s
k/m, die Kreisfrequenz des ungedämpften Oszillators: 2π-mal die Eigenfrequenz f0 in Hertz.
ωd
Gedämpfte Kreisfrequenz
rad/s
Die Kreisfrequenz einer leicht gedämpften freien Schwingung, ω02−β2, etwas unter ω0: 4,97 statt 5 rad/s für den leichten Block.
bc
Kritische Dämpfungskonstante
kg/s
Die Dämpfungskonstante, bei der β=ω0 gilt, gleich 2km: die Grenze zwischen Schwingen und Zurückkriechen. Für die Blöcke der Animation beträgt sie 10 kg/s.
ζ
Dämpfungsverhältnis
β/ω0=b/bc: Unter 1 schwingt das System, bei 1 ist es kritisch gedämpft, über 1 kriecht es zurück. Die Blöcke der Animation haben 0,1, 1 und 3.
Q
Gütefaktor
ω0/(2β), die gespeicherte Energie geteilt durch den Energieverlust pro Radiant. Eine freie Schwingung klingt etwa Q Zyklen lang nach, und der Resonanzgipfel ist Q-mal so hoch wie die statische Antwort.
ω
Antriebskreisfrequenz
rad/s
Die Kreisfrequenz der periodischen Kraft, 2π-mal die Anzahl der Stöße pro Sekunde. Im stationären Zustand schwingt der Oszillator mit dieser Frequenz.
F0
Amplitude der Antriebskraft
N
Der größte Wert der periodischen Kraft F0cosωt. Das Schütteln des oberen Endes einer Feder um Y entspricht F0=kY, was in der Animation 0,5 N sind.
A
Amplitude
m
Die größte Auslenkung in einem Zyklus. Beim freien Abklingen schrumpft sie von ihrem Anfangswert A0 wie A0e−βt; unter einem Antrieb stellt sie sich auf den stationären Wert A(ω) ein.
Quiz
Einsteiger Drei gleiche Blöcke an gleichen Federn werden gemeinsam nach unten gezogen und losgelassen: einer schwach gedämpft, einer kritisch gedämpft, einer stark gedämpft. Welcher kommt zuerst zur Ruhe?
Der schwach gedämpfte
Der kritisch gedämpfte
Der stark gedämpfte
Alle drei im selben Moment
Antwort zeigen
Antwort: B. Der kritisch gedämpfte
Kritische Dämpfung ist genau so stark, dass das Überschwingen gerade aufhört. Bei weniger Dämpfung hüpft der Block weiter; bei mehr verlangsamt der Dämpfer die Rückkehr so sehr, dass der Block nur kriecht.
Einsteiger Eine Schaukel, die einmal angestoßen und dann sich selbst überlassen wird, kommt allmählich zum Stillstand. Warum?
Die Schwerkraft zieht irgendwann nicht mehr an ihr
Ihre Energie wird bei jeder Schwingung ein wenig vernichtet
Reibung und Luftwiderstand wandeln ihre Energie in Wärme um
Ihr natürlicher Rhythmus verlangsamt sich, bis sie stillsteht
Antwort zeigen
Antwort: C. Reibung und Luftwiderstand wandeln ihre Energie in Wärme um
Energie geht nie verloren: Reibung am Drehpunkt und Luftwiderstand wandeln die Energie der Schaukel in Wärme um, bei jeder Schwingung ein bisschen. Dieser Verlust ist Dämpfung. Der Rhythmus bleibt fast gleich, während die Schwingungen kleiner werden.
Einsteiger Du willst eine Schaukel nur mit kleinen Stößen höher bringen. Wie solltest du stoßen?
Einmal pro Schwingung, im Takt des Rhythmus, wenn sie sich von dir wegbewegt
So oft und so schnell du kannst
In zufälligen Momenten, damit sich die Stöße ausgleichen
Nur wenn sie den tiefsten Punkt passiert und sich auf dich zubewegt
Antwort zeigen
Antwort: A. Einmal pro Schwingung, im Takt des Rhythmus, wenn sie sich von dir wegbewegt
Ein Stoß pro Schwingung, jeweils wenn sie sich von dir wegbewegt, führt jedes Mal Energie zu, sodass kleine Stöße eine große Schwingung aufbauen: das ist Resonanz. Derselbe Stoß, während sie sich auf dich zubewegt, nimmt Energie weg.
Mittelstufe Ein 1 kg schwerer Block hängt an einer Feder mit k = 25 N/m. Welche Dämpfungskonstante b macht ihn kritisch gedämpft?
5 kg/s
25 kg/s
10 kg/s
50 kg/s
Antwort zeigen
Antwort: C. 10 kg/s
b_c = 2√(k·m) = 2 × √(25 × 1) = 10 kg/s. Vergisst man den Faktor 2, erhält man 5 kg/s; lässt man die Wurzel weg, ergibt sich 2 × 25 = 50 kg/s, und lässt man beides weg, 25 kg/s.
Mittelstufe Ein schwach gedämpfter Oszillator hat β = b/(2m) = 0,5 s⁻¹ und eine Amplitude von 10 cm bei t = 0. Wie groß ist seine Amplitude bei t = 2 s?
6,1 cm
1,4 cm
5,0 cm
3,7 cm
Antwort zeigen
Antwort: D. 3,7 cm
A = A₀·e^(−βt) = 10 cm × e^(−1) ≈ 3,7 cm. 1/β = 2 s ist die Zeit, um auf 1/e ≈ 37 % abzufallen, nicht auf die Hälfte; behandelt man sie als Halbwertszeit, ergibt sich 5 cm. Die Energie, proportional zu A², klingt wie e^(−2βt) ab: Mit diesem Exponenten erhält man 1,4 cm, und mit e^(−βt/2) ergeben sich 6,1 cm.
Mittelstufe Eine Schaukel mit etwas Reibung wird genau mit ihrer Eigenfrequenz angestoßen, jedes Mal mit dem gleichen kleinen Stoß. Was passiert mit ihrer Amplitude?
Sie wächst ohne Grenze
Sie wächst und bleibt dann konstant, wenn der Energieverlust pro Schwingung der zugeführten Energie entspricht
Sie bleibt so groß wie bei der ersten Schwingung
Sie wächst und schrumpft dann wieder auf null
Antwort zeigen
Antwort: B. Sie wächst und bleibt dann konstant, wenn der Energieverlust pro Schwingung der zugeführten Energie entspricht
Die Dämpfungsverluste wachsen mit der Amplitude, daher wächst die Schwingung nur, bis die Verluste die Zufuhr ausgleichen. Bei einer Dämpfungskraft −b·v ist die stationäre Amplitude A = F₀/(b·ω₀) = Q·F₀/k: groß, aber endlich.
Fortgeschritten Ein schwach gedämpfter Oszillator wird genau mit seiner Eigenkreisfrequenz ω₀ angetrieben. Wie verhält sich seine Auslenkung im stationären Zustand zur antreibenden Kraft?
In Phase mit der Kraft
Eine halbe Periode (180°) hinterher
Eine Viertelperiode voraus (90°)
Eine Viertelperiode hinterher (90°)
Antwort zeigen
Antwort: D. Eine Viertelperiode hinterher (90°)
tan δ = 2βω/(ω₀² − ω²) wird bei ω₀ unendlich, also ist δ = 90°: Die Auslenkung hinkt um eine Viertelperiode hinterher, und die Geschwindigkeit ist in Phase mit der Kraft, daher verrichtet der Antrieb während des ganzen Zyklus positive Arbeit.
Fortgeschritten Zwei Oszillatoren mit gleicher Masse und gleicher Feder unterscheiden sich nur in ihrer Dämpfung: Der eine hat Q = 50, der andere Q = 5. Wie verhält sich der Resonanzgipfel des Oszillators mit Q = 50, wenn beide von der gleichen Kraft angetrieben werden?
Etwa 10-mal höher und 10-mal schmaler
Etwa 10-mal höher und genauso breit
Genauso hoch, aber 10-mal schmaler
Etwa 100-mal höher und 10-mal schmaler
Antwort zeigen
Antwort: A. Etwa 10-mal höher und 10-mal schmaler
Bei Resonanz gilt A = Q·F₀/k, also wächst die Höhe proportional zu Q, während die Breite Δω ≈ ω₀/Q proportional zu 1/Q schrumpft.