Gekoppelte Oszillatoren mit imaginärer Frequenz?

Ein typisches Problem in der Lagrange-Mechanik ist ein System aus zwei gekoppelten Körpern unter kleinen Schwingungen. Der typische Weg, diese Probleme zu lösen, ist:

  • Schreiben Sie die Lagrange-Funktion zB L ( θ 1 , θ 2 , θ ˙ 1 , θ ˙ 2 ) .
  • Erweitern Sie in kleinen Mengen alle Terme nach zweiter Ordnung.
  • Verwenden Sie die Euler-Lagrange-Gleichungen, um die Bewegungsgleichung zu finden.
  • Machen Sie den Ansatz θ 1 = A e ich ω T Und θ 2 = B e ich ω T .
  • Lösen für ω indem verlangt wird, dass die Determinante der relevanten Matrix auf Null geht.

Dieses Verfahren gibt Ihnen normalerweise einen Ausdruck der Form:

ω ± 2 = a ± β
für a > β (beides als positiv gewertet) das bekommen wir dann ω + Kann beides sein ± | ω + | aber diese führen nicht zu linear unabhängigen Lösungen (da Sünde ( | ω + | ) ist keine linear unabhängige Form Sünde ( | ω + | ) ) und die Lösung für θ 1 kann dann geschrieben werden als:
θ 1 = A + e ich | ω + | T + A e ich | ω | T
Allerdings wann a < β die Lösung für ω + funktioniert genauso wie vorher, aber z ω du hast ω = ± ich | ω | was Ihnen zwei Lösungen gibt exp ( | ω | T ) Und exp ( | ω | T ) die linear unabhängig sind. Meine Frage ist daher, wie schreiben wir in diesem Fall die Lösung für auf θ 1 vorausgesetzt, wir können nur zwei Integrationskonstanten haben, aber drei unabhängige Lösungen?

@DanielSank Danke, ich habe die Änderung vorgenommen. Es passiert mit Dämpfung, es passiert auch, wenn in rotierenden Koordinatensystemen gearbeitet wird. Aber ich denke, meine Methode sollte trotzdem so oder so gelten.
Ich denke, Sie sollten die Fälle diskutieren, in denen Sie in der Post eine imaginäre Häufigkeit erhalten, da diese Fälle die Randbedingungen des Problems informieren, die Ihnen wiederum sagen, welche Lösungen Sie beibehalten müssen.
@DanielSank Bist du sicher, dass wir Randbedingungen berücksichtigen müssen? Randbedingungen werden verwendet, um die Integrationskonstanten zu bestimmen - wenn wir sie verwenden, um zu bestimmen, welche Lösung beibehalten wird und welche nicht, dann haben wir effektiv eine dritte Integrationskonstante.
Ich vermute, dass entweder 1) eine der Lösungen die Randbedingungen nicht erfüllt oder 2) die beiden Exponentiale nicht linear unabhängig sind.
@DanielSank Ich denke dann ist es letzteres, dh die Exponentiale sind nicht linear unabhängig. Das heißt, die Exponenten exp ( ich a + β T ) , exp ( β a T ) , Und exp ( β a T ) sind nicht linear unabhängig. Das erscheint mir nicht offensichtlich - aber ich bin immer noch nicht davon überzeugt, dass wir Randbedingungen verwenden können, um zu bestimmen, welche beibehalten werden sollen.
In einem gedämpften System ist die exponentiell abnehmende Lösung erlaubt, aber die wachsende nicht.
Im Allgemeinen erhalten Sie in einem System aus zwei gedämpften Oszillatoren 4 verschiedene komplexe Frequenzen. Ihr Fall (zwei Paare einander entgegengesetzter Frequenzen) sieht entartet aus. Vielleicht hängt es irgendwie mit außergewöhnlichen Punkten nicht-hermitescher Hamiltonianer zusammen?
@DanielSank Ein instabiles Gleichgewicht hat imaginäre Frequenzen und erlaubt sowohl den wachsenden Modus als auch den abnehmenden Modus, die jeweils dem Partikel entsprechen, das aus dem Gleichgewicht rutscht, und einem Partikel, das sich ihm asymptotisch nähert. Bei 2 Freiheitsgraden kann es noch interessanter werden. Das hat nichts mit Verschwendung zu tun. Siehe meine Antwort für mehr.

Antworten (3)

Fixpunkte werden am elegantesten in einem Hamilton-Formalismus behandelt, also werde ich dies hier verwenden (ich werde kurz kommentieren, wie sich dies am Ende in den Lagrange-Formalismus übersetzen lässt).


Betrachten Sie einen Hamilton-Operator H ( θ 1 , θ 2 , P 1 , P 2 ) und einem Fixpunkt (Gleichgewicht) bei θ 10 , θ 20 , P 10 , P 20 , die definiert ist durch θ ˙ 1 = θ ˙ 2 = P ˙ 1 = P ˙ 2 = 0 . Unter Berücksichtigung der Hamilton-Gleichungen kann dies auch so angegeben werden D H = 0 am Gleichgewichtspunkt im Phasenraum. Nehmen wir dann einen Anfangszustand mit einer kleinen Störung aus dem Gleichgewicht θ 10 + δ θ 1 , θ 20 + δ θ 2 , P 10 + δ P 1 , P 20 + δ P 2 wir erhalten einen Satz von Gleichungen für diese Störung der Form

D D T ( δ θ 1 δ P 1 δ θ 2 δ P 2 ) = ( 2 H P 1 θ 1 2 H P 1 2 2 H P 1 θ 2 2 H P 1 P 2 2 H θ 1 2 2 H θ 1 P 1 2 H θ 1 θ 2 2 H θ 1 P 2 2 H P 2 θ 1 2 H P 2 P 1 2 H P 2 θ 2 2 H P 2 2 2 H θ 2 θ 1 2 H θ 2 P 1 2 H θ 2 2 2 H θ 2 P 2 ) ( δ θ 1 δ P 1 δ θ 2 δ P 2 )
Das sieht ein wenig überwältigend aus, aber Sie können dort Ihren bevorzugten Hamilton-Operator einsetzen und Sie werden sehen, dass normalerweise nur wenige der Komponenten dieser Störungsmatrix nicht Null sind. Die obige Differentialgleichung ist eine lineare Differentialgleichung der Form v ˙ = A v deren Lösung man in fast jedem Buch über lineare Algebra findet.

Der vollständige Lösungssatz wird durch einen Ansatz des Formulars gefunden v ( T ) = v 0 e λ T . Daraus ergibt sich, dass die Lösungen haben müssen v 0 ein Eigenvektor der Matrix A mit möglicherweise komplexem Eigenwert λ . Es gibt im Allgemeinen vier Eigenwerte und vier entsprechende linear unabhängige Eigenvektoren von A .

Nun sagt uns die mathematische Theorie der klassischen Mechanik, dass aufgrund der „symplektischen Struktur“ der klassischen Mechanik die Eigenwerte in einer sogenannten loxodromen Menge gruppiert werden. Ein loxodromer Satz von Zahlen in der komplexen Ebene ist ein Satz von Zahlen, der unter Reflexionen sowohl in Bezug auf die reelle als auch auf die imaginäre Achse symmetrisch ist. Dies ist beispielsweise ein loxodromer Satz von vier Zahlen

{ a + ich β , a + ich β , a ich β , a ich β }
Wo a , β sind reelle Zahlen. Ein Fixpunkt mit einem solchen Eigenwertsatz entspricht einem System mit einer außer Kontrolle geratenen Schwingung e ( a + ich β ) T . Dieser Fall kann nur in zwei Freiheitsgraden und höher auftreten, Sie werden dafür keine Intuition von einem Freiheitsgrad bekommen! Ein besonderer Fall, in dem dies passieren kann, sind rotierende Systeme, bei denen das Teilchen durch ein anziehendes Potential zur Mitrotation und Oszillation gezwungen wird, sich aber aufgrund der Zentrifugalkraft auch spiralförmig ausdehnt.

Eine der anderen Möglichkeiten für die Eigenwerte ist eine "entartete" loxodromische Menge der Form

{ a , a , ich β , ich β }
mit a , β wieder echt. Diese entsprechen einem System, bei dem eine Störung in einer Richtung eine stabile Schwingung verursacht, während sie in der anderen Richtung instabil ist und die Bahn zu driften beginnt. Die loxodrome Symmetrie der Eigenwerte entspricht der Reversibilität des klassisch-mechanischen Systems. Das bedeutet, dass man für jede Bewegung in eine allgemeine Richtung in gewisser Weise auch die Umkehrung finden kann; für jeden e a T Ausreißer, es gibt einen e a T langsame Annäherung an das Gleichgewicht (dies entspricht einem Teilchen, das genau die richtige Energie hat, um ein potentielles Maximum zu erreichen, und langsam darauf zukriecht).


Nun zu deiner Frage: dein Ansatz θ 1 = A e ich ω T , θ 2 = B e ich ω T entspricht dem allgemeineren Ansatz v ( T ) = v 0 e λ T im Hamiltonschen Formalismus und sollte Ihnen auch einen loxodromen Satz von Frequenzen geben.

Aber Sie sollten jetzt sehen, dass, wenn wir eine Anfangsbedingung an einer allgemeinen Position und Geschwindigkeit nahe dem Gleichgewicht haben, es eine Linearkombination der Basis aller vier linear unabhängigen Eigenrichtungen sein wird - unabhängig davon, ob die loxodromische Menge entartet ist oder nicht -degenerieren. (Denken Sie daran, dass Sie vier Zahlen angeben müssen θ 1 , θ 2 , θ ˙ 1 , θ ˙ 2 als Ausgangszustand!!)

Wenn Ihre Anfangsbedingung (und Ihr dynamisches System) real ist, stellt die Loxodromie der Menge sicher, dass sich Ihre Anfangsbedingungen tatsächlich in einer realen Basis von Evolutionen kombinieren, die für eine nicht entartete Menge wie folgt geschrieben wurden

1 2 ich ( e ( a + ich β ) T e ( a ich β ) T ) , 1 2 ( e ( a + ich β ) T + e ( a ich β ) T ) , 1 2 ich ( e ( a + ich β ) T e ( a ich β ) T ) , 1 2 ( e ( a + ich β ) T + e ( a ich β ) T )
und ähnlich für die entarteten Mengen. Obwohl also die Frequenzen für ein System mit zwei oder mehr Freiheitsgraden gemischt-komplex werden können (beachten Sie rein imaginäre oder reale), retten die Symmetrien des mechanischen Systems den Tag und Sie erhalten immer noch eine ordnungsgemäße, nicht komplexe Entwicklung.

Freie Schwingungen

Bei freier kleiner Schwingung um ein stabiles Gleichgewicht sind die Eigenfrequenzen des Systems immer positiv, ω 2 = a ± β > 0 .

Die Lagrangefunktion eines Systems mit N Freiheitsgrade und mit einem stabilen Gleichgewicht bei Q 0 = ( Q 0 , 1 , , Q 0 , N ) kann geschrieben werden als

L = 1 2 η ˙ T M ( Q 0 ) η ˙ + 1 2 η T v ( Q 0 ) η ,
Wo M ich J = A = 1 N M A R A Q ich R A Q J Und v ich J = 2 v Q ich Q J Sind N × N Matrizen u η , η ˙ Sind N -Komponenten-Spaltenmatrizen, die die kleinen Verschiebungen aus dem Gleichgewicht bezeichnen. Der erste Term im Lagrangian entspricht der kinetischen Energie, also ist er für jeden immer positiv η ˙ 0 , somit M ist positiv definit . Die Bedingung für ein stabiles Gleichgewicht impliziert, dass die Matrix v ( Q 0 ) ist auch positiv definit, also ist der zweite Term im Lagrangian positiv für η 0 .

Bewegungsgleichungen liest

M ( Q 0 ) η ¨ + v ( Q 0 ) η = 0.
Wir suchen nach Lösungen des Typs η = ρ cos ( ω T + ϕ ) und dann erhalten wir, dass die charakteristischen Vektoren ρ erfüllen
( v ( Q 0 ) ω 2 M ( Q 0 ) ) ρ = 0.
Links multiplizieren mit ρ T und rechts vorbei ρ wir erhalten
ω 2 = ρ T v ( Q 0 ) ρ ρ T M ( Q 0 ) ρ .
Da beide v ( Q 0 ) Und M ( Q 0 ) sind also positiv definite Matrizen ω 2 > 0 für jede nicht triviale Lösung ρ .

Wenn das Gleichgewicht neutral ist, dh es gibt wiederherstellende verallgemeinerte Kräfte in allen Richtungen des Konfigurationsraums, außer in einigen Richtungen, in denen die Kraft null ist, dann ist die Matrix v ( Q 0 ) ist nur positiv semidefinit . In diesem Fall, ω 2 0 .

Gedämpfte Schwingungen

Der Lagrange-Formalismus ist für dissipative Systeme nicht geeignet, da er davon ausgeht, dass das System anhand seiner kinetischen und potentiellen Energien beschrieben wird. Für einige Systeme kann man jedoch den Formalismus passend ändern, um dissipative Kräfte einzubeziehen. Insbesondere können dissipative Kräfte proportional zu den Teilchengeschwindigkeiten aus der sogenannten Dissipationsfunktion abgeleitet werden F . Euler-Lagrange-Gleichungen werden dann gelesen

D D T L η ˙ ich L η ich + F η ich = 0.

Die Dissipationsfunktion ist in den Geschwindigkeiten quadratisch und hängt mit der Rate der dissipierten Energie zusammen, was bedeutet, dass sie nicht negativ ist. Für kleine Schwingungen um das Gleichgewicht kann geschrieben werden als

F = 1 2 η ˙ T F ( Q 0 ) η ,
und die Bewegungsgleichungen sind
M ( Q 0 ) η ¨ + F ( Q 0 ) η ˙ + v ( Q 0 ) η = 0.
Unwahrscheinlich der freie Fall, dieses System ist nicht einfach zu lösen, da im Allgemeinen die Matrizen M , F Und v können nicht gleichzeitig diagonalisiert werden. Eine Ausnahme tritt auf, wenn die Dämpfung auch von den Teilchenmassen abhängt, die Matrizen gleichzeitig diagonalisiert werden können und die Bewegungsgleichungen in den Normalkoordinaten entkoppeln,
ζ ¨ ich + F ich ζ ˙ ich + ω ich 2 ζ ich = 0 ,
Wo F ich 0 , ich = 1 , , N sind die Eigenwerte von F ( Q 0 ) . Diese Gleichungen können leicht mit komplexen Funktionen gelöst werden z = z 0 exp ( ich Ω ich T ) und am Ende nur seinen wahren Teil nimmt. Stecken z in die Bewegungsgleichung gibt
Ω ich = ± ω ich 2 F ich 2 4 ich F ich 2 .
Bei kleiner Dämpfung können wir den quadratischen Term in vernachlässigen F ich damit die Lösung die schöne Form hat
ζ ich = C ich e ich F ich T / 2 cos ( ω ich T + φ ich ) .
Wenn F ich 2 / 4 ω ich 2 dann schwingt das System überhaupt nicht und diese Bereiche werden als kritisch gedämpft und überdämpft bezeichnet.

Siehe die Kommentare oben, wir werden negativ ω 2 B. bei vorhandener Dämpfung oder beim Arbeiten in einem rotierenden Koordinatensystem. Ich bin mir nicht sicher, ob Ersteres mit der Lagrange-Mechanik einfach zu handhaben ist, aber der Fall von rotierenden Koordinatensystemen ist es.
@Quantumspaghettification Alle diese Verfahren, die Sie beschreiben, gelten für konservative Systeme. In diesem Fall sind die charakteristischen Frequenzen aus dem von mir gezeigten Grund immer positiv. Bei Dämpfung sind die Dinge nicht so einfach, da für nicht-dissipative Systeme Euler-Lagrange-Gleichungen angenommen werden. Der Lagrangian selbst geht von einem monogenen System (durch Potentiale beschriebenes System) aus. Bei rotierenden Koordinaten können Sie beliebige verallgemeinerte Koordinaten verwenden, aber Sie müssen die Lagrangian in Bezug auf ein Inertialsystem schreiben.
@Diracology Der Punkt ist offensichtlich, dass das Gleichgewicht nicht stabil ist, Stabilität nirgendwo im OP diskutiert wird und dies daher nicht erforderlich ist ω 2 > 0 .
@Diracology Sie können Ihren Lagrangian beliebig schreiben und ihn und die entsprechenden Gleichungen rein mathematisch behandeln. Monogenität ist nur für die Äquivalenz mit dem d'Alembert-Prinzip erforderlich. Die Transformation in ein korotierendes System ermöglicht es Ihnen, die Zeitabhängigkeit zB eines rotierenden Potentials im Lagrange auf eine zeitunabhängige Beschreibung zu reduzieren, Sie müssen nur das richtige verallgemeinerte Potential finden, das alle Rotationseffekte erfasst (es gibt es).
@Void Wäre Stabilität nicht eine notwendige Bedingung, um Schwingungen zu untersuchen? Mindestens neutrale Stabilität (Fluchtkraft in eine Richtung des Konfigurationsraums, während die anderen Rückstellkraft haben) wäre erforderlich. Bei neutraler Stabilität gilt ω 2 0 .

Wenn es nur zwei unabhängige Lösungen gibt (die anderen beiden Lösungen sind eine wachsende Funktion und haben keine physikalische Bedeutung) mit den Anfangsbedingungen jeweils zwei unabhängige Entscheidungen und der Koeffizient vor der wachsenden Lösung ist Null. Wenn der Eigenwert einen negativen Realteil hat, dann ist die Lösung gedämpft. Wenn der Realteil des Eigenwerts Null ist, dann oszillierend.