Der Grenzfall

 Obige Überlegungen gelten unter der Voraussetzung der „schwachen Dämpfung“ die wir durch unsere Rechnungen zur Bedingung

μ  <  2 m ω

präzisieren konnten. Die gedämpfte Kreisfrequenz ω* konvergiert gegen 0, wenn μ gegen 2 m ω strebt. Im Grenzfall μ = 2 m ω liegt keine unendliche Oszillation mehr vor, die Funktion ist, wenn man so sagen will, auf eine Periode gestreckt. Um Lösungen auch für größere Dämpfungen zu finden, modifizieren wir unseren Ansatz: Wir behalten die Funktion e−λt bei, ersetzen aber die Schwingung cos(ωt − φ) durch eine andere Funktion, die wir nun als zeitabhängige Skalierung der Funktion e−λt auffassen (während wir bisher e−λt als Skalierung einer Schwingung angesehen haben). Sei also q :    von der Form

q(t)  =  c(t) e−λt  für alle t  ∈  

mit einer zweimal differenzierbaren Funktion c :   . Dann gilt für alle t  ∈  :

q′(t) =  − λ q(t)  +  c′(t) e−λ t   (sodass c′(t) e−λt = λ q(t) + q′(t))
q″(t) =  − λ q′(t)  −  λ c′(t) e−λ t  +  c″(t) e−λt
=  − λ q′(t)  −  λ (λ q(t) + q′(t))  +  c″(t) e−λt
=  − λ2 q(t)  −  2 λ q′(t)  +  c″(t) e−λt

Für Schwingungen c(t) ist der dritte Summand ein skalares Vielfaches von q, sodass wir ihn in den ersten Summanden aufnehmen können (in der Form − ω*2 q). Nun sehen wir eine neue Möglichkeit: Der dritte Summand verschwindet, wenn c eine lineare Funktion der Form c(t) = a t + b ist. Es gilt dann

(+)  q″(t)  =  − λ2 q(t)  −  2 λ q′(t)  für alle t  ∈  

Für μ = 2 m ω wird die Differentialgleichung q″(t) = − ω2 q(t) − μ/m q′(t) zu

(#)  q″(t)  =  − ω2 q(t)  −  2 ω q′(t)

Mit λ = ω ist also q :    mit

q(t)  =  (a t + b) e−ωt  für alle t  ∈  

für alle a, b  ∈   eine Lösung von (#). Es gilt

q(0)  =  b,  q′(0)  =  − λ q(0) + c′(t)  =  − ω b + a

Das AWP mit den Startwerten x0 und v0 wird also gelöst durch

b  =  x0,  a  =  v0 + ω x0

Beispiel (aperiodischer Grenzfall)

Für m = 1, μ = 2, k = 1/2 gilt μ = 2 m ω. Für x0 = 1, v0 = 2 erhalten wir die Lösung q :    mit

q(t)  =  ((v0 + ω x0) t + x0) e−ωt  =  ((2 + 1/2) t + 1) e−ωt  für alle t  ∈  

ellipsen1-AbbIDoscillator_dampedvecfield_4

q ist maximal bei t = 4/(1 + 22) ∼ 1,04 und fällt dann exponentiell ab

ellipsen1-AbbIDoscillator_dampedvecfield_5

Der gleiche Ausschnitt der Lösung im Phasenraum (mit Farben wie für q). Die Lösung verbleibt im vierten Quadranten, sobald der Impuls negativ wird. Eine minimale Verkleinerung von μ führt zu einer unendlichen Spiralbewegung.