Differenzialgleichungen zur Beschreibung von Federschwingungen

Ein Körper, der an einer Feder befestigt ist, führt nach einer Auslenkung eine Schwingung durch. Der Ort des Körpers wird durch die zeitabhängige Ortskoordinate y(t) beschrieben, deren Gleichung gefunden werden soll.
Im Folgenden werden mit einer derartigen Anordnung gedämpfte und ungedämpfte Schwingungen untersucht.

Federschwinger

Auf den Körper wirken folgende Kräfte:

die Trägheitskraft: $m a = m y ″ (t)$ mit $m$ - Masse, $a$ - Beschl.
die geschwindigkeitsproportionale Reibungskraft:
$β v = β y' (t)$ mit $β$ - Reibungskoeffizient, $v$ - Geschwindigkeit
die rücktreibende Federkraft (hookesches Gesetz):
$D y (t)$ mit $D$ - Federkonstante, $y$ - Auslenkung

Die Summe der betrachteten Kräfte ist gleich der von außen wirkenden Kraft, hier also 0. Es gilt demzufolge:
$m y ″ (t) + β y' (t) + D y (t) = 0$

Dividiert man diese Differenzialgleichung 2. Ordnung durch m, so ergibt sich
$y ″ (t) + \frac{β}{m} y' (t) + \frac{D}{m} y (t) = 0 (*)$ .

Diese Gleichung entspricht mit den Zuordnungen $q = \frac{β}{m} u n d r = \frac{D}{m}$ der Gleichung $f ″ (x) + q f' (x) + r f (x) = 0$ .

Gleichung $(*)$ wird mit dem Lösungsansatz $y (t) = e^{k ‌ t}$ gelöst.

Es gilt $y' (t) = k e^{k ‌ t} u n d y ″ (t) = k^{2} e^{k ‌ t}$ .

Durch Einsetzen in die Differenzialgleichung $(*)$ erhält man
$k^{2} e^{k t} + \frac{β}{m} k e^{k t} + \frac{D}{m} e^{k t} = e^{k t} (k^{2} + \frac{β}{m} k + \frac{D}{m}) = 0$
und daraus nach Division durch $e^{k t}$ die charakteristische Gleichung
$k^{2} + \frac{β}{m} k + \frac{D}{m} = 0$
mit den Lösungen
$k_{1 / 2} = - \frac{β}{2 m} \pm \sqrt{{(\frac{β}{2 m})}^{2} - \frac{D}{m}} .$

Ist der Reibungskoeffizient $β$ zwar nicht vernachlässigbar, aber dennoch hinreichend klein, so überwiegt im Radikand der Wurzel der Anteil $\frac{D}{m}$ , der Radikand wird negativ.

Wegen $q = \frac{β}{m}$ (s. o.) ergibt sich als allgemeine Lösung
$y (t) = e^{- \frac{β}{2 m} t} (c_{1} \cos ω t + c_{2} \sin ω t) m i t c_{1}, c_{2} \in ℝ u n d ω = \sqrt{\frac{D}{m} - {(\frac{β}{2 m})}^{2}}$ .

Da $\frac{β}{2 m}$ als Dämpfungskonstante $δ$ bezeichnet wird, heißt die Lösung schließlich $y (t) = e^{- δ t} (c_{1} \cos ω t + c_{2} \sin ω t) m i t c_{1}, c_{2} \in ℝ$ .

Diese Gleichung besagt, dass die Federanordnung in eine gedämpfte Schwingung gebracht werden kann.

Beispiel: Es soll nun eine Federanordnung mit $m = 1 k g, D = 1,04 k g \cdot s^{- 2} u n d β = 0,4 k g \cdot s^{- 1}$ betrachtet werden.

Hier gilt $\begin{array}{l} \frac{D}{m} = 1,04 \cdot s^{- 2} > > {(\frac{β}{2 m})}^{2} = {(0,2 s^{- 1})}^{2} = 0,04 s^{- 2} \\ δ = \frac{β}{2 m} = 0,2 s^{- 1} \\ ω = \sqrt{1,04 s^{- 2} - 0,04 s^{- 2}} = 1 s^{- 1} . \end{array}$

Die Schwingungsgleichung für diese Federanordnung lautet:
$y (t) = e^{- 0,2 s^{- 2} \cdot t} (c_{1} \cos (1 s^{- 1} \cdot t) + c_{2} \sin (1 s^{- 1} \cdot t))$

Wird die Masse aus der Ruhelage heraus $(y (0) = 0)$ mit einer Geschwindigkeit von $10 m s^{- 1}$ in die positive y-Richtung weggestoßen $(y' (0) = 10 m s^{- 1}),$ so lassen sich die Parameter $c_{1} u n d c_{2}$ der allgemeinen Lösung bestimmen:

$y (0) = 0$ in die allgemeine Lösung einsetzen:
$y (0) = 0 = e^{- δ \cdot 0} (c_{1} \cos (ω \cdot 0) + c_{2} \sin (ω \cdot 0)) = 1 \cdot (c_{1} \cdot 1 + c_{2} \cdot 0) = c_{1} = 0.$

Erste Ableitung der Lösung:
$\begin{array}{l} y' (t) = [e^{- δ t} (c_{1} \cos ω t + c_{2} \sin ω t)]' \\ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ ‌ = - δ e^{- δ t} (c_{1} \cos ω t + c_{2} \sin ω t) + e^{- δ t} (- c_{1} ω \sin ω t + c_{2} ω \cos ω t) \\ y' (0) = 10 m s^{- 1} = - δ e^{- δ \cdot 0} (c_{1} \cos (ω \cdot 0) + c_{2} \sin (ω \cdot 0)) \\ + e^{- δ \cdot 0} (- c_{1} ω \sin (ω \cdot 0) + c_{2} ω \cos (ω \cdot 0)) \\ = - δ (c_{1} + 0) + ω (0 + c_{2}) \\ y' (0) = 10 m s^{- 1} = - δ \cdot c_{1} + ω c_{2} . \\ D a c_{1} = 0 i s t, g i l t ω c_{2} = 10 m s^{- 1} b z w . c_{2} = \frac{10 m s^{- 1}}{ω} = \frac{10 m s^{- 1}}{1 s^{- 1}} = 10 m . \end{array}$

Die Lösung des Anfangswertproblems lautet:
$y (t) = 10 m \cdot e^{- 0,2 s^{- 1} \cdot t} \cdot \sin (1 s^{- 1} \cdot t)$

Lösung des Anfangswertproblems einer gedämpften Schwingung

Nachdem eine gedämpfte Schwingung untersucht wurde, soll nun für die gleiche Federanordnung geklärt werden,

welche Schwingung ausgeführt wird, wenn keine Reibung vorhanden wäre und
wie groß der Reibungskoeffizient bzw. die Dämpfungskonstante sein müsste, damit die Anordnung überhaupt nicht schwingen kann.

Im ersten Fall gilt $β = 0$ und damit auch $δ = 0$ Aus diesen Werten errechnet sich die Kreisfrequenz $ω = 1,020 s^{- 1}$ und damit ergibt sich für die allgemeine Lösung der Schwingungsgleichung
$y (t) = e^{- 0 \cdot t} (c_{1} \cos ω t + c_{2} \sin ω t) = c_{1} \cos 1,020 s^{- 1} \cdot t + c_{2} \sin 1,020 s^{- 1} \cdot t .$

Unter den bereits betrachteten Anfangsbedingungen ergibt sich die ungedämpfte Schwingung $y (t) = 9,806 m \cdot \sin 1,020 s^{- 1} \cdot t$ (siehe Graph 1 in der folgenden Abbildung).

Die Frequenz der Schwingung ist hier etwas höher als bei der gedämpften Schwingung (siehe Graph 2 in der folgenden Abbildung).

Die Federanordnung kann nicht mehr schwingen, wenn die Reibung und damit der Wert für $β$ mindestens so groß geworden ist, dass die charakteristische Gleichung $k^{2} + \frac{β}{m} k + \frac{D}{m} = 0$ eine reelle Doppellösung besitzt.

In diesem Falle gilt $\frac{D}{m} = {(\frac{β}{2 m})}^{2} o d e r β = 2 \sqrt{D \cdot m}$ .

Für die betrachtete Federanordnung liegt dieser Grenzwert bei $β = 2,040 k g s^{- 1} b z w . b e i δ = 1,020 s^{- 1}$ .

Als allgemeine Lösung erhält man dann $y (t) = c_{1} e^{- 1,020 s^{- 1} \cdot t} + c_{2} t \cdot e^{- 1,020 s^{- 1} \cdot t}$ .

Werden auch hier die bekannten Anfangsbedingungen benutzt, so erhält man $y (t) = 10 m \cdot t \cdot e^{- 1,020 s^{- 1} \cdot t}$ (siehe Graph 3 in der folgenden Abbildung).

Bild

Stand: 2010
Dieser Text befindet sich in redaktioneller Bearbeitung.

Mathe Note verbessern?

Kostenlos bei Duden Learnattack registrieren und ALLES 48 Stunden testen.

Kein Vertrag. Keine Kosten.

40.000 Lern-Inhalte in Mathe, Deutsch und 7 weiteren Fächern
Hausaufgabenhilfe per WhatsApp
Original Klassenarbeiten mit Lösungen
Deine eigene Lern-Statistik
Kostenfreie Basismitgliedschaft

Mathe kostenlos lernen

Verwandte Artikel

Differenzialgleichungen zur Beschreibung von Federschwingungen

Rechtliches

Einloggen