Differentialgleichungen · Dämpfung

Dämpfung

Dämpfung entzieht einer Schwingung Energie. Verändere ζ und vergleiche schwingende Rückkehr mit einer Rückkehr ohne Überqueren der Ruhelage. Die Zeit wird in T₀ gemessen, der ungedämpften Periodendauer.

01

Dämpfung in dieser Visualisierung

Dämpfung ist der Energieverlust, durch den eine Schwingung abklingt. Für eine Masse an Feder und Dämpfer gilt mx′′+cx′+kx=0m x'' + c x' + k x = 0, geschrieben als x′′+2ζω0x′+ω02x=0x'' + 2\zeta\omega_0 x' + \omega_0^2 x = 0 mit der Eigenkreisfrequenz ω0=k/m\omega_0 = \sqrt{k/m} und dem Dämpfungsgrad ζ=c/(2mk)\zeta = c/(2\sqrt{mk}). Unter 1 ist die Bewegung unterkritisch gedämpft (sie schwingt), bei 1 kritisch gedämpft, über 1 überkritisch gedämpft (sie kriecht zurück). Bei einer senkrecht hängenden Masse verschiebt die Schwerkraft lediglich die Ruhelage um mg/k nach unten; x wird von dort gemessen. Die Seite lässt die Masse bei x = 1 aus der Ruhe los und fragt, ob mehr Dämpfung sie früher zur Ruhe bringt.

Die Wurzeln bestimmen die Bewegung; eine doppelte Wurzel liefert auch einen Term t·e^(λt)
λ2+2ζω0λ+ω02=0,λ=ω0(−ζ±ζ2−1),σ=−max⁡Re⁡λ\lambda^2 + 2\zeta\omega_0\lambda + \omega_0^2 = 0, \qquad \lambda = \omega_0\left(-\zeta \pm \sqrt{\zeta^2 - 1}\right), \qquad \sigma = -\max \operatorname{Re}\lambda

Rechts im Bild liegt die Wurzelebene. Das Produkt der beiden Wurzeln ist ω02\omega_0^2. Solange die Bewegung schwingt, liegen sie deshalb auf dem Kreis mit Radius ω₀: Der Realteil −ζω0-\zeta\omega_0 bestimmt das Abklingen, der Imaginärteil ωd=ω01−ζ2\omega_d = \omega_0\sqrt{1-\zeta^2} die Frequenz. Dämpfung verlängert also auch jede Schwingung. Bei ζ = 1 treffen sich die Wurzeln bei −ω₀. Danach läuft eine nach links, während die andere langsam gegen 0 zurückkehrt, für große ζ ungefähr zu −ω0/(2ζ)=−k/c-\omega_0/(2\zeta) = -k/c. Der langsamste Anteil der Bewegung klingt wie e−σte^{-\sigma t} ab. Diese Abklingrate σ, der Abstand der am weitesten rechts liegenden Wurzel von der imaginären Achse, ist genau bei ζ = 1 maximal. Im Diagramm von x(t) zeigt die gestrichelte blaue Kurve diesen Anteil: bei unterkritischer Dämpfung die Einhüllende ±e−ζω0t/1−ζ2\pm e^{-\zeta\omega_0 t}/\sqrt{1-\zeta^2} (die Kurve berührt sie zwischen ihren Extrema), bei überkritischer Dämpfung den Term der langsamen Wurzel.

„Zur Ruhe kommen“ braucht eine Definition. Die Seite verwendet T1%T_{1\%}, den frühesten Zeitpunkt, ab dem |x| dauerhaft innerhalb von 1 % der Anfangsauslenkung bleibt. Der Streifen unter x(t) vergrößert dieses Band. Die Extrema liegen bei tn=nπ/ωdt_n = n\pi/\omega_d mit ∣x(tn)∣=e−ζω0tn|x(t_n)| = e^{-\zeta\omega_0 t_n}. Daraus folgt exakt, welche Extrema außerhalb des Bandes liegen. Die letzte Rückkehr wird anschließend durch Bisektion der geschlossenen Lösung bestimmt; die Anzeige ist daher numerisch. Bei kleinerem ζ verlängert sich die Einschwingzeit sprunghaft, sobald ein weiteres Extremum das Band verlässt. Daher hat die Kurve von T1%T_{1\%} gegen ζ Sprünge. Ihr Minimum liegt nicht bei kritischer Dämpfung, sondern bei ζ∗=ln⁡100/π2+ln⁡2100≈0.826\zeta^* = \ln 100/\sqrt{\pi^2 + \ln^2 100} \approx 0.826, wo das einzige Überschwingen genau −1 % berührt: 0.667 T₀ statt 1.057 T₀. Kritische Dämpfung ist nur unter den Bewegungen am schnellsten, die x = 0 nie überqueren, und nur beim Start aus der Ruhe. Ein anderes Band verschiebt ζ* (bei 2 % auf ungefähr 0.780).

02

Was auffällt

  • Zweimal „am schnellsten“: Die Abklingrate σ ist genau bei ζ = 1 maximal; die Einschwingzeit ist bei ζ* etwas darunter am kürzesten. Der eingegebene Wert 0.826 ist nicht ζ*: Er führt zu 1.0015 % Überschwingen, und die Einschwingzeit springt auf 0.896 T₀. Die Seite bewahrt ζ* exakt.
  • Verdoppelt man einen großen Dämpfungsgrad, verdoppelt sich ungefähr die Rückkehrzeit: Die Zeitkonstante der langsamen Wurzel beträgt ungefähr c/k.
  • x(t) ist in allen drei Bereichen die exakte Lösung. Der Bereich eines eingegebenen Dämpfungsgrads wird anhand der exakt eingegebenen Zahl bestimmt; 0.99999999999999999 ist also unterkritisch. Die Flüssigkeit im Zylinder ist nur schematisch dargestellt.
03

Bedienung

  • Verschiebe den ζ-Regler (er rastet bei 0.1, 1 und 10 ein), gib einen Bruch oder eine Dezimalzahl ein oder drücke ζ = 1 beziehungsweise ζ*. Eine Wurzel lässt sich entlang ihrer Bahn ziehen. Bei −ω₀ rastet sie ein und wechselt nur dort den Ast. Ziehe den roten Zeitzeiger, um die Masse zu bewegen.
  • Der Präsentationsmodus öffnet eine Vollbildansicht: ← → ändern ζ um einen Reglerschritt; H blendet Abklingrate und Einschwingzeiten aus.
04

Verwandte Themen

Im Spur-Determinanten-Diagramm der Phasenporträt-Seite entspricht der gedämpfte Oszillator der waagerechten Linie Δ = ω₀².

Weiterführende Literatur: Wikipedia: Damping (Englisch); OpenStax, University Physics Volume 1, 15.5 Damped Oscillations (Englisch).