Wie behandle ich die Lagrange-Funktion bei einem starren Körper?

Hier ist Übung 1.11 aus Goldsteins Classical Mechanics 3rd Edition (die erste, nachdem die Lagrange-Funktion im Wesentlichen abgeleitet wurde):

Aufgabe 1.11: Betrachten Sie eine gleichförmige dünne Scheibe, die ohne Schlupf auf einer horizontalen Ebene rollt. Eine horizontale Kraft wird auf die Mitte der Scheibe und in einer Richtung parallel zur Mitte der Scheibe ausgeübt.

  • Leiten Sie die Lagrange-Gleichungen ab und finden Sie die verallgemeinerte Kraft;

  • Diskutieren Sie die Bewegung, wenn die Kraft nicht parallel zur Scheibenebene wirkt.

Jetzt ist mir klar, dass der Ansatz, den ich vorschlagen werde, ein bisschen "primitiv" in dem Sinne ist, dass er nicht viel voraussetzt, aber ich wollte die Lagrange-Funktion von "ersten Prinzipien" ableiten (da ich gerade mit Lagrange-Funktionen begonnen habe). .

Die Scheibe besteht aus N Teilchen, so hat es 3 N Freiheitsgrade. Wir haben offensichtlich mehrere Einschränkungen, die die Bewegung der Scheibe stark einschränken:

  • Es muss sich bewegen X j Ebene;
  • Es ist ein starrer Körper;
  • Es rollt ohne zu rutschen.

Ich möchte den Winkel verwenden θ mit dem X Achse als meine verallgemeinerte Koordinate (erste Frage: Kann ich das tun? Es sieht so aus, als ob ich aus der Nicht-Schlupf-Bedingung leicht ableiten kann X C M aus θ ˙ und umgekehrt, aber die Bewegungsgleichungen erscheinen irgendwie seltsam, wenn ich das mache). Die kinetische Energie des Systems ist T = 1 2 M ( X ˙ C M 2 + θ ˙ 2 R 2 ) = M θ ˙ 2 R 2 wegen des rutschfesten Zustands. Jetzt kann ich Euler-Lagrange-Gleichungen in dieser Form verwenden:

D D T ( T θ ˙ ) T θ F θ = 0

Wo F θ ist die verallgemeinerte Kraft. Hier ist meine Hauptfrage: Wie behandle ich die generalisierte Kraft, wenn ich einen starren Körper habe? Was gebe ich ein F θ ? Ich kann die Euler-Lagrange-Gleichungen für jedes Teilchen im Körper schreiben, aber das scheint irgendwie nicht sehr nützlich zu sein. Ich habe darüber nachgedacht, die Kraft auf ein einzelnes Partikel über den Körper zu integrieren, aber irgendetwas fühlt sich irgendwo nicht richtig an. Was vermisse ich?

Da Sie Italiener sind, könnten Sie sich meine Vorlesungsunterlagen zur analytischen Mechanik für Studenten ansehen (ich sollte sie ins Englische übersetzen, aber ich habe keine Freizeit) ... science.unitn.it/~moretti/dispense.html Sie könnten sie herunterladen die erste "Dispense", die Sie auf dieser Seite finden. Ich denke, Sie sollten sich für Kapitel 8 interessieren, insbesondere für Abschnitt 8.3.
Super Valter, danke! Ich schau es mir mal an! :) (und es ist auch eine großartige Ressource für den Rest des Kurses!)

Antworten (3)

Da Sie die Lagrangefunktion aus Grundprinzipien ableiten möchten, lohnt es sich, einige Worte über das Prinzip virtueller Werke zu verlieren und wie sich daraus die Lagrangefunktionsgleichungen ergeben.

Der springende Punkt dabei ist, dass der schlupffreie Abrollzustand garantiert, dass die Zwangskräfte keine Arbeit verrichten , da der Bodenkontaktpunkt immer in Ruhe ist. Das D'Alèmbert-Prinzip ist also erfüllt. Daher gelten die verallgemeinerten Gleichungen von Lagrange:

D D T T Q ˙ ich T Q ich = Q ich
halten, wo die Q ich s sind bestimmte Funktionen, die nur von äußeren Kräften abhängen. Ich denke, dass Sie sich dessen bewusst sind.

Finden Sie nun die verallgemeinerten Koordinaten Q ich und die damit verbundenen verallgemeinerten Kräfte Q ich ist Teil der Modellierung des Problems. In unserem Fall lassen wir nur Verschiebungen in einer Dimension zu und der rutschfeste Zustand garantiert, dass der "Winkel" < θ < + , oder die Position des Massenschwerpunkts < X < sind beides gute (globale) Koordinaten.

Nun, was ist mit Q = F θ ? Wenn Sie das System als zusammengesetzt sehen möchten N massive Teilchen, mit Positionen X ( N ) = X ( N ) ( Q 1 , Q 2 , , Q N ) Erfüllung der Newtonschen Gleichungen

M ( N ) X ¨ ( N ) = F ( N ) + R ( N ) ,
bei dem die R ( N ) keine virtuelle Arbeit machen, ich denke, dass Goldstein das in diesem Fall erklärt
Q ich = N = 1 3 N F ( N ) X ( N ) Q ich .
Die rechte Seite dieser Gleichung ist die sogenannte "virtuelle Arbeit" der wirkenden Kräfte, d.h. Q ich δ Q ich , ist die Arbeit, die geleistet wird, wenn das System eine durch parametrisierte Verschiebung erfährt Q ich Q ich + δ Q ich , wobei alle anderen Koordinaten (sowie die Form der zeitlichen Beschränkung, die im vorliegenden Fall konstant ist) festgehalten werden. (Auch wenn F ( N ) = U X ( N ) , das sieht man leicht Q ich = U Q ich )

Im vorliegenden Fall ist das virtuelle Werk F δ X = δ ( F X ) , So Q ich Ist U Q ich , Wo U = F X . Schreibe den Lagrangian als L = T U , mit der richtigen Form von T , können Sie den Lagrangian leicht ableiten.

Ich möchte betonen, dass dieses Verfahren, obwohl es physikalisch rigoros ist, nur eine Heuristik ist, da es keinen a priori Grund (zumindest in meiner Diskussion) gibt, warum die intuitive Identifizierung von Q ich mit der virtuellen Arbeit sollen die korrekten Bewegungsgleichungen wiedergegeben werden. Wenn Sie die Bewegungsgleichungen bereits in irgendeiner Form haben, ist der einfachste Weg, um die Lagrangian herauszufinden, wahrscheinlich:

  1. Schreiben Sie die LHS der verallgemeinerten Gleichungen von Lagrange.
  2. Setzen Sie die (gegebenen) Bewegungsgleichungen ein.
  3. Eliminieren Sie die Zwangskräfte.
  4. Erfahren Q ich und möglicherweise der Lagrange.

Normalerweise sollte Kraft in Ihrem Lagrange enthalten sein:

L = M θ ˙ 2 R 2 + F θ R

Aber man hätte es auch als einbinden können F θ = F R

Beide entsprechen:

2 M R 2 D θ ˙ D T F R = 0

Wenn Sie sich nicht sicher sind, wie Kräfte nach Koordinatentransformationen aussehen, schlage ich vor, dass Sie mit einem Lagrange in kartesischen Koordinaten beginnen und die Transformation in die verallgemeinerten Koordinaten durchführen, die Sie verwenden möchten.

In kartesischen Koordinaten hat man dann

L = M 2 ( X ˙ 2 + j ˙ 2 + R 2 ω 2 ) F X

Da der Zylinder die ganze Zeit auf dem Boden liegt j ˙ = 0 , und die Bedingung "kein Rutschen" ergibt X = R θ Und θ ˙ = ω . Ändern der Koordinaten zu θ gibt das dann X ˙ = R θ ˙ So

L = M 2 ( R 2 θ ˙ 2 + R 2 θ 2 ) F R θ = M R 2 θ ˙ 2 F R θ

Die generalisierte Kraft ist jetzt F θ = L θ = F R und die Bewegungsgleichungen lautet

2 M R 2 θ ¨ = F R θ ¨ = F 2 M R