12. Kapitel

Analytische Lösung von Differentialgleichungen 2. Ordnung

Lösungsweg für homogene Differentialgleichungen 2. Ordnung

Mit dem Kräftegleichgewicht eines Feder-Masse-Reibungs-Systems

()

wird eine Differentialgleichung 2. Ordnung gebildet:

()

Sie ist vom Typ wieder eine inhomogene Differentialgleichung, deren Lösung die Lösung der homogene Differentialgleichung voraussetzt. Daher betrachten wir zunächst die Lösung der homogenen DGl.:

()

bzw.

Diese Differentialgleichung beinhaltet die Ableitungen der Ortskoordinate ausschließlich in erster Potenz mit konstanten Faktoren. Für einen solchen Typ gibt es ein einfaches Lösungsverfahren, das auf der Eigenschaft der Exponentialfunktion, bei der Differentiation unverändert zu bleiben, beruht. Also versuchen wir mit dem Ansatz

()

eine Lösung für die vorliegende Differentialgleichung zu finden.
Um mit diesem Lösungsansatz bedienen zu können, benötigen wir die erste und zweite Ableitung von , deren Herleitung unter Beachtung der Kettenregel erfolgt:

()

in eingesetzt, ergibt:

()

Jetzt wird der Vorteil dieser Lösungsmethode sichtbar; die e-Funktion des Ansatzes e^λ·t kann keinesfalls dazu beitragen, dass die verschwindet und kann deshalb gekürzt werden:

()

wird charakteristische Gleichung genannt. Die Gleichung ist deshalb "charakteristisch", weil ihre Lösung das charakteristische Zeitverhalten der Lösung der DGl beschreibt. Für den von uns untersuchten Fall eines Feder-Masse-Reibungs-Systems, das durch eine DGl 2. Ordnung beschrieben wird, ist die charakteristische Gleichung eine quadratische Gleichung. Zu ihrer Lösung verwenden wir den Wurzelsatz von VIETA.

()

Dabei stellen wir fest, dass es nicht nur eine, sondern zwei Lösungen gibt! Das ist aber auch ganz o.k., da wir eine Differentialgleichung mit einer zweiten Ableitung lösen wollen, brauchen wir auch zwei Integrationsschritte und zwei Konstanten, mit deren Hilfe wir die Anfangsbedingungen für Geschwindigkeit und Ort einarbeiten können.

Im strengen mathematischen Sinn ist die Lösung

s = e^{λ \cdot t}

eine spezielle oder partikuläre Lösung, die eine von vielen Lösungsmöglichkeiten darstellt. Siehe dazu auch Lineare Differentialgleichungen

Mit dem Wissen, dass die charakteristische Gleichung zwei Lösungen hat, können wir die Lösung der DGl. allgemein angeben, indem entsprechend erweitert wird. Das geschieht durch eine Überlagerung beider Lösungen:

()

Dabei spielen die beiden Konstanten A un B die gleiche Rolle wie die Konstanten C oder K in den früheren Lösungen.

on/off

Wie üblich werden in einem zweiten Schritt die Konstanten mittels der Anfangsbedingungen für Ort und Geschwindigkeit auf den konkreten Fall abgestimmt.
Die Einarbeitung der Anfangsbedingungen läuft nun etwas komplizierter ab, als wir das bisher kennen gelernt haben. Da infolge des gewählten Lösungsansatzes die Lösung der Differentialgleichung sofort auf die Zeitfunktion der Ortskoordinate () geführt hat, müssen wir, um die Anfangsgeschwindigkeit einzuarbeiten, die Zeitfunktion der Geschwindigkeit herleiten. Dies geschieht durch Differentiation der Ortskomponente:

()

Jetzt können die Anfangsbedingungen (t = 0!) eingesetzt werden. Aus folgt:

(a)

und aus

(b)

Dies ist ein lineares Gleichungssystem, bestehend aus zwei Gleichungen mit zwei Unbekannten A und B. Ein solches Gleichungssystem zu lösen ist nicht schwer und liefert für die beiden gesuchten Konstanten die folgenden Werte:

(a)

und

(b)

Da wir uns ja in erster Linie für die Zeitfunktion der Ortskoordinate s interessieren, genügt es, die Ergebnisse von a) und b) in einzusetzen:

()

Recht unübersichtlich? Ja, darum greife ich jetzt etwas vor und führe die folgenden Abkürzungen ein (deren physikalische Bedeutung folgt weiter unten!).

()

worin

Mit diesen Abkürzungen lautet

()

Nun ist das mit den Wurzeln ein eigen Ding! Nach unserem Schulwissen kann eine Wurzel nur aus positiven Zahlen gezogen werden. Was aber, wenn das Argument unter der Wurzel negativ ist?
Untersuchen wir mal, welche Werte das Argument annehmen kann und welche Wurzeln dazu gehören:

1. Fall:

Das Argument ist positiv und beide Lösungen sind stets negativ.

2. Fall:

Das Argument ist negativ. Die Wurzel wird imaginär und es liegen zwei komplexe Nullstellen vor. Deren Realteil stets negative Werte hat.

3. Fall:

Das Argument ist exakt Null. Es liegt eine doppelte Nullstelle vor, die stets negative Werte hat.

Zusammen gefasst:

()

Nehmen wir uns die drei Fälle der Reihe nach vor. Der Einfachheit halber wählen wir die Startgeschwindigkeit v₀ = 0. Das vereinfachte die Rechnung, ist aber nicht zwingend!

1. Fall: Einsetzen beider Lösungen in und Vereinfachen führt auf

()

Die hier verwendeten Hyperbolicus-Funktionen sinh(x) und cosh(x) stehen für die Zusammenfassungen der Exponentialfunktionen e^±x.
Das Zeitverhalten eines gespannten Feder-Masse-Reibungs-Systems mit diesen Charakteristika ist ein einfaches exponentielles Ausschleichen nach dem Loslassen. Deshalb wird dieser Fall auch Kriech-Fall genannt.

2. Fall: Dieser Fall ist der interessanteste! Einsetzen der komplexen Lösung in ergibt

()

Die EULERsche-Formel

stellt den Zusammenhang zwischen einer Exponential-Funktion mit imaginärem Argument und den trigonometrischen Funktionen her. So erklären sich auch die Zusammenhänge zwischen den hyperbolicus- und den trigonometrischen-Funktionen:

Werden jetzt die Wurzelargumente in aus dem 1. Fall durch die des 2. Falles ersetzt, dann wird zu

Für imaginäre Argumente werden nämlich die Hyperbolicus-Funktionen sinh(x) und cosh(x) "fast" zu den ganz gewöhnlichen Winkelfunktionen sin(x) und cos(x). Grundlage hierfür bildet die EULERsche Formel.

on/off

Und damit erschließt sich auch die physikalische Bedeutung der verwendeten Abkürzungen (). Stell Dir vor, die Größe δ hat einen sehr kleinen Wert. Dann können in einige Summanden wegfallen, da sie viel kleiner als die anderen Summanden sind. Übrig bleibt dann

()

die Gleichung für eine gedämpfte, sinus-förmige Schwingung. Daher ist die Wahl der Bezeichnungen naheliegend:

δ die Dämpfung gemessen in s^-1,
ω₀ die Kreisresonanzfrequenz gemessen in s^-1
f₀ = ω₀/2π die Eigen- oder Resonanzfrequenz gemessen in Hz bzw.

Hier schwingt das System nach dem Loslassen exponentiell gedämpft aus. Dieser Fall wird darum Schwing-Fall genannt.

3. Fall: Dieser Fall ist wegen der doppelten Nullstelle ein Sonderfall. Hierfür ist ein anderer Lösungsansatz als ihn vorschlägt, erforderlich. In der Literatur wird ein spezieller Ansatz gezeigt, der ebenfalls zwei linear unabhängige Gleichungen zur Einarbeitung der Anfangsbedingungen erzeugt werden:

()

Zu diesem Zweck wird wieder differenziert (Produkt- und Kettenregel!):

()

Mit den Anfangsbedingungen

(a)

und

(b)

folgt die Bewegungsgleichung für den 3. Fall mit v₀ = 0

()

Dieser Fall wird aperiodischer Grenzfall genannt. Das Zeitverhalten unterscheidet sich kaum von dem des 1. Falles. Bei ganz genauem Hinschauen ist ein leichtes Unterschleifen vor dem finalen Ausschleichen zu beobachten.

Im Programmbeispiel kann die Bewegung eines Feder-Masse-Reibungs-Systems als Zeitfunktion nachvollzogen werden. Der blaue Graph stellt die Bewegung des Verknüpfungspunktes zwischen Federende und Masse dar. Der graue Graph zeigt den Dämpfungsverlauf, also den Abbau der initialen Energie des Systems.
Die beiden Schieberegler gestatten die Eingabe der Resonanzfrequenz ω₀ und der Dämpfung δ. Unterhalb des START-Buttons wird angezeigt, zu welchem Fall die jeweilige Kombination von ω₀ und δ führt.

Bitte einen Augenblick Geduld
während das Programm geladen wird!

Abb. Lösungen der DGl. 2. Ordnung

Soweit die Lösung ohne Störfunktion, die die Lösung der homogenen Differentialgleichung mit eingearbeiteten Anfangsbedingungen darstellt. Wie gehen wir aber vor, wenn das System von außen, z.B. eine sinusförmige Anregung gestört wird. Diese Aufgabe führt auf die Aufstellung und Lösung einer inhomogenen Differentialgleichung.

Lösungsweg für inhomogene Differentialgleichungen 2. Ordnung

Für die Lösung der inhomogenen Aufgabe einer DGl. 2. Ordnung ist die Methode der Variation der Konstanten aus Aufwandsgründen nicht geeignet. Eine alternative Lösungsmethode bietet sich an, wenn eine spezielle Lösung (auch partikuläre Lösung genannt) der inhomogenen DGl bekannt ist. Gibt es eine solche, kann die Lösung der inhomogenen DGl durch die Überlagerung von homogener und spezieller Lösung gefunden werden ()

()

Bei der Lösung der inhomogenen Aufgabe wollen wir uns auf den Fall einer sinus-förmigen Anregung beschränken. gibt die inhomogenen DGl wider, wobei â Spitzenamplitude und ω die Frequenz der Anregungsschwingung, 2δ die Dämpfung und ω₀ die Eigenfrequenz des Feder-Masse-Systems darstellen:

()

Wie findet man nun eine spezielle Lösung? Wir gehen davon aus, dass ein System, das sinus-förmig angeregt wird, auch sinus-förmig antwortet. Daher versuchen wir den folgenden Ansatz:

()

Da wir nicht wissen, ob die Lösung richtiger durch eine sinus- oder cosinus-Funktion beschrieben werden kann, verwenden wir beide Funktionen im Ansatz. Welche der beiden Funktionen am Ende die richtige ist, oder zu welchem Bruchteil eingeht, überlassen wir den Konstanten C und D, die natürlich im Folgenden noch zu bestimmen sind.
Dies geschieht, indem wir den Lösungsansatz und seine Ableitungen

()

in die einsetzen:

()

Ein Koeffizientenvergleich nach sin(ω·t) bzw. cos(ω·t) liefert uns die erforderlichen zwei Gleichungen zur Bestimmung der zwei Unbekannten C und D:

()

Vereinfachen und Ordnen

()

Nach der Lösung des Gleichungssystems stehen die beiden Konstanten C und Dzur Verfügung und werden in eingesetzt:

()

Damit haben wir nun eine spezielle Lösung der inhomogenen Aufgabe gefunden. Nun fehlt nur noch die Überlagerung, d.h. Addition, mit der homogenen Lösung, die uns aber schon aus den Gleichungen bekannt ist:

()

Um sich das Verhalten eines Systems im Überblick zu verschaffen, ist die Betrachtung der sog. Übertragungsfunktion Ü interessant. Sie stellt das Verhältnis der Schwingungsamplitude des Systems zur Anregungsamplitude dar ():

()

Ü = \frac{\overset{‸}{s}}{\overset{‸}{a}} = \frac{1}{\sqrt{{(2 δ \frac{ω}{ω_{0}})}^{2} + {(1 - {(\frac{ω}{ω_{0}})}^{2})}^{2}}}

Sehr gut ist hier die Resonanzüberhöhung zu verfolgen. Kommt nämlich die Anregungsfrequenz ω in die Nähe der Eigenfrequenz des Systems ω₀, dann wird der Nenner sehr klein, was dann zu einem großen Quotienten führt. Steigt aber die Anregungsfrequenz deutlich über die Eigenfrequenz, dann wird der Nenner sehr groß, was dann zu einer sehr kleinen Amplitude führt. Das ist das typische Verhalten eines Tiefpasses.
Es zeigt sich, dass die Resonanzüberhöhung umso größer ist, je kleiner die Dämpfung δ ist.

Abb. Übertragungsverhalten

Symbol	Bedeutung
F_W	Wirkkraft [N]
F_R	Reibungskraft [N] ( bzw. )
F_T	Trägheitskraft [N] ()
F_F	Federkraft [N]
s	Strecke, Weg [m]
s_h	Strecke homogene Lösung [m]
s_p	Strecke partikuläre (speziell) Lösung [m]
s₀ / v₀	Anfangswerte Ort [m] / Geschwindigkeit [m/s]
$\overset{‸}{a}$	Spitzenwert (der Amplitude)
m	Masse [kg]
r	Reibzahl [kg/s], z.B. für STOKESsche Reibung ()
n	Federkonstante [N/m]
δ	Dämpfung [1/s] ()
ω₀	Kreisfrequenz der Eigenschwingung [rad/s] ()
τ	Zeitkonstante des Ausklingvorgangs [s] ()
λ / λ₁ / λ₂	Lösungsansatz charakteristische Gleichung, Eigenwerte [m] ()
t	Zeit [s]

Analytische Lösung von Differential­gleichungen 2. Ordnung

Lösungsweg für homogene Differentialgleichungen 2. Ordnung

Lösungsweg für inhomogene Differentialgleichungen 2. Ordnung

Analytische Lösung von Differentialgleichungen 2. Ordnung