Lagrange- und Hamilton-Formalismus für gedämpfte Pendel

Angenommen, ich habe ein vereinfachtes Pendel (masselose Schnur der Länge l , Kugel aus Masse M und etwas Nicht-Null-Volumen). Ich möchte die Bewegungsgleichungen herleiten, aber ich möchte den Luftwiderstand berücksichtigen. Ich weiß, wie es geht, indem ich Newtons 2. Gesetz anwende, aber gibt es eine Möglichkeit, dies mit Lagrange oder Hamilton zu tun? Was wäre in diesem Fall die Form von Lagrangian/Hamiltonian?

Antworten (3)

Für eine allgemeine Dämpfungskraft gibt es kein Prinzip der kleinsten Wirkung. Sie können jedoch die Bewegungsgleichungen aus dem d'Alembert-Prinzip erhalten. Das Verfahren wird hier erklärt und die Bewegungsgleichungen sind

D D T L Q ˙ ich L Q ich = Q ich P ,
Wo L = T v ist die Lagrange-Funktion und
Q ich P a F a R a Q ich ,
ist die sogenannte verallgemeinerte Kraft, die von der Dämpfungskraft abhängt F a Wirkung von Teilchen a .

Da wir für die allgemeine Dämpfung kein Variationsprinzip formulieren können, L allein reicht nicht aus, um die dynamische Entwicklung des Systems zu bestimmen. Daher gibt es keinen Hamiltonoperator H befriedigend

δ S = δ T 1 T 2 [ ich P ich Q ˙ ich H ( Q , P , T ) ] D T = 0.
Denken Sie daran, dass die kanonischen Gleichungen auch aus dem obigen Variationsprinzip erhalten werden können, sodass es keine Funktion gibt H die kanonischen Gleichungen und die Bewegungsgleichungen des Systems (die die Dämpfung berücksichtigen sollen) gleichzeitig erfüllen.

Im Fall der linearen Dämpfung können Sie der Bateman-Idee folgen, die in der Antwort von ZeroTheHero erwähnt wird, oder die Rayleigh-Dissipationsfunktion verwenden

Es hängt von der Form der Dämpfung ab, aber wenn Sie schreiben

L = e a T ( 1 2 M Q ˙ 2 v ( Q ) )
dann erhält man eine Bewegungsgleichung der Form
e a T ( M Q ¨ + M a Q ˙ + D v D Q ) = 0
die je nach Reibungsmodell für die Reibung verantwortlich sein können .

Der Hamilton-Operator würde auf die übliche Weise, aber eindeutig berechnet werden L hängt explizit davon ab T Die Energie wird also in Ihrem System nicht konserviert.


Bearbeiten: geändert a a nach einem korrekten Kommentar.

Warum e A T ? Ich sehe nicht, woher es kam. Ich habe auch vergessen zu erwähnen, dass ich möchte, dass die Dämpfung vom Quadrat der Geschwindigkeit abhängt. Nun, der Grund, warum ich eine solche Frage stelle, ist, dass ich letztendlich ein Doppelpendel mit Luftwiderstand simulieren möchte und das Ableiten von Gleichungen aus Lagrangian einfacher zu sein scheint als die Kraftmethode.
Ich denke, Sie wollen das entgegengesetzte Zeichen für a .
@JG In der Tat. Ich habe gerade die Korrektur vorgenommen.
@G.Fil Soweit ich weiß, das e a T Faktor ist phänomenologisch, dh es wird beobachtet, dass er den geeigneten Dämpfungsterm erzeugt. Ich kenne keinen Trick, der eine Dämpfung quadratisch einfängt Q ˙ : Möglicherweise müssen Sie die verallgemeinerten Bewegungsgleichungen verwenden, bei denen der Term, der nicht aus einem Potential stammt, "von Hand" hinzugefügt wird.
Es ist nicht phänomenologisch; Es gibt einen tiefen Grund für die Wahl. Das ist ein bisschen übertrieben, aber die e a T Faktor ist analog zu dem | G | Faktor in der Allgemeinen Relativitätstheorie, wenn wir die Zeit als einen betrachten 1 -dimensionale Riemannsche Mannigfaltigkeit. Dann Q ist ein Feld auf dieser Mannigfaltigkeit, mit der Bewegungsgleichung Q ˙ Term proportional zum Christoffel-Term Γ 00 0 = 1 2 G 00 0 G 00 = 0 ln G , woraus die G = e 2 a T Wahl folgt. (Eigentlich ist es ein de Sitter-Verteiler.)
@JGwow! Ich habe gerade etwas gelernt! Aber warum sollten Sie auf diese Weise über Zeit nachdenken, und bietet dies zusätzliche Erkenntnisse, die bei der Bearbeitung zusätzlicher (komplizierterer) Fälle helfen würden?

In gewisser Weise ja; Sie müssen die Kraft des Luftwiderstands so schreiben, als ob sie von einem geschwindigkeitsabhängigen Potential abgeleitet werden könnte.

Ein gängiges Beispiel hierfür ist ein Teilchen in einem Magnetfeld . Die Kraft des Magnetfelds auf das Teilchen hängt von seiner Geschwindigkeit ab, daher wird das Potential geschrieben als A v .