Ist die üblicherweise gelehrte Lösung für erzwungene harmonische Bewegung nur eine spezielle Lösung?

Nehmen wir an, wir haben eine Masse an einer Feder, die von einer Zwangsfunktion angetrieben wird. Angesichts des Hook-Gesetzes F = k X , und eine Zwangsfunktion von

F ( T ) = F 0 Sünde ( ω T ) .
Wir können schreiben:

M D 2 X D T 2 = k X + F 0 Sünde ( ω T )

Alle physikalischen Ressourcen, auf die ich gestoßen bin, gehen davon aus, dass die Bewegung der Feder der ausgeübten Kraft folgt, und präsentieren die Lösung als eine Form von:

X = C Sünde ( ω T )

Sie gehen dann typischerweise auf den Ersatz x in der Differentialgleichung ein und erhalten:

C = F 0 M ( ω 0 2 ω 2 )

Das ist eine ziemlich coole Formel. Ich wollte wirklich verstehen, warum, wenn ω > ω 0 , C wird negativ und unsere Bewegung ist genau außer Phase mit unserer Kraft. Dies ist für mich nicht intuitiv, und um es besser zu verstehen, habe ich beschlossen, eine numerische Analyse durchzuführen.

Ich habe meine Masse zunächst im Ruhezustand bei Position Null gestartet und zu meinem Entsetzen ergab meine numerische Analyse folgende Ergebnisse:Erzwungene harmonische Oszillationslösung

Offensichtlich kann die Bewegung der Masse nicht durch einen einzigen Sinus beschrieben werden! Was ist denn hier los? Nachdem ich mich ein wenig an den Haaren gezogen hatte, stellte ich fest, dass meine numerische Analyse tatsächlich richtig war und es die analytische Lösung war, die fehlte. Die vollständige Lösung unserer Bewegungsgleichung lautet:

X ( T ) = A Sünde ( ω 0 T ) + B cos ( ω 0 T ) + F 0 Sünde ( ω T ) M ( ω 0 2 ω 2 )

Und wenn wir die Anfangsbedingungen richtig aufstellen, stimmt diese analytische Lösung mit der numerischen Lösung überein! Unsere frühere Lösung ist ein Sonderfall dieser Lösung - aber die Anfangsbedingungen müssen auf ganz bestimmte Werte gesetzt werden, damit dies geschieht.

Meine Frage ist also - was zum Teufel ist hier los? Ist diese Lösung einfach nicht wichtig oder nicht relevant? Es scheint mir, dass Sie nicht einmal das Verhalten sehen werden, das in den meisten Physikressourcen gezeigt wird, wenn die Anfangsbedingungen nicht genau richtig sind. Kommt die gelehrte Lösung in realen Anwendungen normalerweise einfach nicht zustande, oder übersehe ich etwas?

Gute Arbeit, dies selbst numerisch zu erkunden. Ich muss sagen, dass ich überrascht bin, dass die homogene Komponente in Ihrem Kurs/Lehrbuch nicht deutlich gemacht wurde.

Antworten (2)

Normalerweise ist der (sinusförmig) angesteuerte harmonische Oszillator gedämpft , und die ersten beiden Teile Ihrer Lösung (die von den Anfangsbedingungen abhängen, während der dritte Term dies nicht tut) sind transient, dh nach kurzer Zeit nicht relevant. Das ist die Lösung

X ( T ) = F 0 Sünde ( ω T ) M ( ω 0 2 ω 2 )
nicht die "vollständige" Lösung der Bewegungsgleichung sein kann, wird schon dadurch deutlich, dass sie nicht von den Anfangsbedingungen abhängt X ( 0 ) Und X ˙ ( 0 ) , wie es die vollständige Lösung einer Differentialgleichung zweiter Ordnung sollte.

In Ihrem Fall eines ungedämpften angetriebenen harmonischen Oszillators sind die vom Anfangszustand abhängigen Teile nicht transient, da keine Dämpfung vorhanden ist, um sie zu entfernen.

Da das Fehlen von Dämpfung in fast allen Fällen eine Idealisierung ist, wird fast kein erzwungener Oszillator in der realen Welt das Verhalten zeigen, das Sie simuliert haben.

Offensichtlich kann die Bewegung der Masse nicht durch einen einzigen Sinus beschrieben werden! Was ist denn hier los?

Die allgemeine Lösung des einfachen harmonischen Oszillators ist die Summe der ungezwungenen Antwort (homogene Lösung) und der erzwungenen Antwort.

Die homogene Lösung ist

X H ( T ) = X H ( 0 ) cos ( ω 0 T ) + X ˙ H ( 0 ) ω 0 Sünde ( ω 0 T )

Somit ist für Anfangsbedingungen ungleich Null die allgemeine Lösung für eine sinusförmige Zwangsfunktion mit Kreisfrequenz ω ω 0 ist die Summe zweier Sinuskurven, eine bei der Eigenfrequenz und eine bei der Antriebsfrequenz.

Bei Dämpfung zerfällt die homogene Lösung mit der Zeit als e T / τ es bleibt im Wesentlichen nur die Antriebsfrequenzkomponente für T > 5 τ