Varianten der klassischen Mechanik/ Erweiterung auf mehrdimensionale Systeme

Die Erkenntnisse bei der Verallgemeinerung auf eindimensionale Vielteilchensysteme können wir beinahe unmittelbar auf mehrdimensionale Systeme übertragen. Beim d'Alembert'schen Prinzip der virtuellen Arbeit finden wir z.B.

$\underset{i=1}{\sum ^{D\cdot N}}({\dot{p}}_i-F_i)\delta x_i=0$.

wieder. Hier bezeichne D die Anzahl der Komponenten, die Orts-, Geschwindigkeits- oder Kraftvektoren besitzen mögen (in der dreidimensionlen Welt der Mechanik ist in der Regel D=3). Jede dieser Komponenten trägt im Allg. zur gesamten virtuellen Arbeit bei. Betrachten wir z.B. dreidimensionale Vielteilchensysteme, dann haben wir bei N Teilchen $D\cdot N=3N$ Beiträge zur virtuellen Arbeit. D.h. auch hier werden die Beiträge zur gesamten virtuellen Arbeit einfach addiert. Bei der Einführung verallgemeinerter Koordinaten, $x_i=x_i(q_1,...,q_M,t)(i=1,...,D\cdot N)$, ist es möglich, dass deren Anzahl M kleiner ist als $D\cdot N$: $M\le D\cdot N$. Denn evtl. gibt es $\nu$ sog. »Nebenbedingungen«, die die Anzahl der sog. »Freiheitsgrade« auf $f=D\cdot N-\nu$ reduzieren. Beim quasi-zweidimensionalen Problem eines Körpers, der eine schiefe Ebene herunter rutscht (d.h. D=2 und N=1), benötigt man zwei Komponenten zur Beschreibung der vektoriellen Größen wie z.B. der Kräfte, d.h. $D\cdot N=2$. Der Körper kann sich jedoch nur entlang der schiefen Ebene bewegen (was die Nebenbedingung darstellt), also entlang einer Dimension, d.h. es gibt nur einen Freiheitsgrad: f=1. Die Anzahl der Nebenbedingungen $D\cdot N-f$ beträgt also gleichermaßen Eins. Die Anzahl M der gewählten verallgemeinerten Koordinaten kann aber immer noch größer sein als die wirkliche Anzahl f der Freiheitsgrade des betrachteten Systems: $M\ge f$, wie z.B. $M=D\cdot N$, im betrachteten Beispiel also 2. Bei M>f, d.h. wenn wir mehr (nämlich M) Komponenten der virtuelle Verrückung $\delta q_i$ als Systemfreiheitsgrade f haben, hat dies zur Folge, dass wir in

$0=\underset{i=1}{\sum ^M}[\underset{j=1}{\sum ^{D\cdot N}}{f}_j\frac{\partial x_j}{\partial q_i}-(\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_i}-\frac{\partial L}{\partial q_i})]\delta q_i$.

nicht mehr darauf schließen können, dass die Koeffizienten der virtuellen Verrückungen $\delta q_i$ alle einzeln verschwinden. Die virtuellen Verrückungen sind dann nämlich nicht mehr alle voneinander linear unabhängig. In der Regel können aber z.B. mit Hilfe der Methode der sog. »Lagrange'schen Multiplikatoren« auch diese Schwierigkeiten überwunden, d.h. die Nebenbedingungen berücksichtigt werden. Bei M=f, also einer optimalen Koordinatenwahl (was aber evtl. nicht immer gelingen wird), treten diese Probleme hingegen nicht auf, d.h. wir können dann ganz analog zum Kapitel über die (eindimensionalen) Vielteilchensystemen vorgehen.

D'Alembert'sches Prinzip, Nebenbedingungen und Zwangskräfte

Bleiben wir bei einem Körper, der im Schwerefeld der Erde reibungslos eine schiefen Ebene herab gleitet. Der Versuchsaufbau ist in der Fig. 1 skizziert.

Datei:Varianten der klass Mechanik Schiefe Ebene.jpg Fig. 1: Körper, der sich im Gravitationsfeld entlang einer schiefen Ebene bewegt.

Die Bewegung vollziehe sich dabei in der xy-Ebene, wobei $x_1=x$ und $x_2=y$ seien. Das kartesiche Koordinatensystem werde durch die Einheitsvektoren ${\hat{e}}_1={\hat{e}}_x=\left(\begin{array}{c}1\\ 0\end{array}\right)$ und ${\hat{e}}_2={\hat{e}}_y=\left(\begin{array}{c}0\\ 1\end{array}\right)$ aufgespannt, d.h. ${\hat{e}}_i^2=1,i=1,2$ und ${\hat{e}}_1\cdot {\hat{e}}_2=0$. Man erkennt in der Fig. 1 sofort, dass die Bewegung am einfachsten in einem Koordinatensystem beschrieben werden kann, dessen 1-Achse entlang der schiefen Ebene (in der xy-Ebene) verläuft. Ein solches Koordinatensystem geht aus dem ursprünglichen durch eine Drehung um den Winkel $\phi$ hervor und ist gleichermaßen kartesisch. Es werde von den Einheitsvektoren ${\hat{e}}_{q_1}$ und ${\hat{e}}_{q_2}$ aufgespannt: ${\hat{e}}_{q_i}^2=1,i=1,2$ und ${\hat{e}}_{q_1}\cdot {\hat{e}}_{q_2}=0$. Bei Drehungen ist der Zusammenhang zwischen den Basisvektoren ${\hat{e}}_{q_1}$, ${\hat{e}}_{q_2}$ und ${\hat{e}}_1$, ${\hat{e}}_2$ der folgende:

${\hat{e}}_{q_1}=({\hat{e}}_{q_1}\cdot {\hat{e}}_1){\hat{e}}_1+({\hat{e}}_{q_1}\cdot {\hat{e}}_2){\hat{e}}_2=\underset{i=1}{\sum ^2}({\hat{e}}_{q_1}\cdot {\hat{e}}_i){\hat{e}}_i=\cos \phi {\hat{e}}_1+\sin \phi {\hat{e}}_2$, ${\hat{e}}_{q_2}=({\hat{e}}_{q_2}\cdot {\hat{e}}_1){\hat{e}}_1+({\hat{e}}_{q_2}\cdot {\hat{e}}_2){\hat{e}}_2=\underset{i=1}{\sum ^2}({\hat{e}}_{q_2}\cdot {\hat{e}}_i){\hat{e}}_i=-\sin \phi {\hat{e}}_1+\cos \phi {\hat{e}}_2$,

was man wegen der Skalarprodukt-Eigenschaft ${\hat{e}}_{q_i}\cdot {\hat{e}}_j=\cos ({\hat{e}}_{q_i},{\hat{e}}_j)$ für Einheitsvektoren durch Ablesen der Winkel $({\hat{e}}_{q_i},{\hat{e}}_j)$ zwischen den jeweiligen Vektoren ${\hat{e}}_{q_i}$ und ${\hat{e}}_j$ aus der Fig. 1 (und z.B. durch Anwenden der Additionstheoreme für Sinus und Kosinus) erhält. Wenn ein Punkt $(x_1,x_2)$ im Koordinatensystem der ${\hat{e}}_i$ durch einen Punkt $(q_1,q_2)$ im Koordinatensystem der ${\hat{e}}_{q_i}$ beschrieben werde, dann gilt also folgender Zusammenhang zwischen den beiden Koordinaten-Paaren bzw. den zugehörigen Ortsvektoren $\overrightarrow{x}$ und $\overrightarrow{q}$:

$\overrightarrow{q}=\underset{i=1}{\sum ^2}{q}_i{\hat{e}}_{q_i}=(q_1\cos \phi -q_2\sin \phi ){\hat{e}}_1+(q_1\sin \phi +q_2\cos \phi ){\hat{e}}_2=$ $\left(\begin{array}{c}{q}_1\cos \phi -q_2\sin \phi \\ q_1\sin \phi +q_2\cos \phi \end{array}\right)=\left(\begin{array}{cc}\cos \phi & -\sin \phi \\ \sin \phi & \cos \phi \end{array}\right)\left(\begin{array}{c}{q}_1\\ q_2\end{array}\right)$,

wobei dieser Vektor gleich

$\overrightarrow{x}=\underset{i=1}{\sum ^2}{x}_i{\hat{e}}_i=\left(\begin{array}{c}{x}_1\\ x_2\end{array}\right)$,

sein muss, da beide Ortsvektoren ja denselben Punkt beschreiben (s. Fig. 1), d.h. es gilt

$\left(\begin{array}{c}{x}_1\\ x_2\end{array}\right)=\underset{\_}{D}\left(\phi \right)\left(\begin{array}{c}{q}_1\\ q_2\end{array}\right)$

mit der Matrix $\underset{\_}{D}\left(\phi \right)=\left(\begin{array}{cc}\cos \phi & -\sin \phi \\ \sin \phi & \cos \phi \end{array}\right)$ für eine Drehung um die z- bzw. 3-Achse ${\hat{e}}_3={\hat{e}}_z$ (die auf der xy-Ebene senkrecht stehe; der Vektor zeige in Fig. 1 dabei in Richtung des Betrachters) mit dem Winkel $\phi$ (im Gegenuhrzeigersinn). Die Koordinaten von $\overrightarrow{q}$ lassen sich auch umgekehrt durch jene von $\overrightarrow{x}$ ausdrücken:

$\left(\begin{array}{c}{q}_1\\ q_2\end{array}\right)={\underset{\_}{D}}^{-1}\left(\phi \right)\left(\begin{array}{c}{x}_1\\ x_2\end{array}\right)=\left(\begin{array}{cc}\cos \phi & \sin \phi \\ -\sin \phi & \cos \phi \end{array}\right)\left(\begin{array}{c}{x}_1\\ x_2\end{array}\right)=\left(\begin{array}{c}{x}_1\cos \phi +x_2\sin \phi \\ -x_1\sin \phi +x_2\cos \phi \end{array}\right)$,

wobei das Inverse ${\underset{\_}{D}}^{-1}\left(\phi \right)$ von $\underset{\_}{D}\left(\phi \right)$ aus letzterer Matrix durch Vorzeichenumkehr beim Drehwinkel hervorgeht, da ja einfach nur um den gleichen Winkel »zurückgedreht« werden muss (sodass also diesmal die Drehung im Uhrzeigersinn, d.h. mit $-\phi$ erfolgt): ${\underset{\_}{D}}^{-1}\left(\phi \right)=\underset{\_}{D}(-\phi )$. Man überzeuge sich leicht, dass dann tatsächlich ${\underset{\_}{D}}^{-1}\left(\phi \right)\underset{\_}{D}\left(\phi \right)=\underset{\_}{D}\left(\phi \right){\underset{\_}{D}}^{-1}\left(\phi \right)=\underset{\_}{1}=\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right)$ gilt.

Ob wir jetzt das D'Alembert'sche Prinzip mit Hilfe der alten Koordinaten $\overrightarrow{x}$ oder der neuen Koordinaten $\overrightarrow{q}$ beschreiben, darf aber am Resultat dieser Betrachtungen, d.h. den Bewegungsgleichungen für den Körper auf der schiefen Ebene, nichts ändern. Das eine Koordinatensystem (nämlich jenes von $\overrightarrow{q}$) ist aber der Problemstellung angepasster als das andere. Die Nebenbedingung für die Bewegung entlang der schiefen Ebene lässt sich nämlich im System von $\overrightarrow{q}$ leicht beschreiben: $q_2=0$ , während im System von $\overrightarrow{x}$ die gleiche Bedingung folgendermaßen lautet: $0=q_2=-x_1\sin \phi +x_2\cos \phi$. Die Zeitableitung des Impulses kann in beiden Koordinatensystemen dargestellt werden:

$\dot{\overrightarrow{p}}=m\underset{i=1}{\sum ^2}{\ddot{x}}_i{\hat{e}}_i=m\underset{i=1}{\sum ^2}{\ddot{q}}_i{\hat{e}}_{q_i}$.

Weil die Basisvektoren ${\hat{e}}_{q_i}$ wie die Vektoren ${\hat{e}}_i$ zeitlich konstant sind (denn der Winkel $\phi$, der die Schräge der schiefen Ebene angibt, verändert sich ja nicht mit der Zeit), haben wir keine zusätzlichen Terme durch Differenziation der Basisvektoren nach der Zeit erhalten.

Auf die Masse m wirke eine Gravitationskraft ${\overrightarrow{F}}_G=-mg{\hat{e}}_2$ mit der konstanten Beschleunigung g anti-parallel zur 2-Achse, d.h. Letzterer entgegen gerichtet.

Das D'Alembert'sche Prinzip ergibt im System der neuen Koordinaten $\overrightarrow{q}$ unter der Nebenbedingung $q_2=0$ und somit auch $\delta q_2=0$ bzw. $\delta \overrightarrow{q}=\delta q_1{\hat{e}}_{q_1}$ (statt $\delta \overrightarrow{q}=\delta q_1{\hat{e}}_{q_1}+\delta q_2{\hat{e}}_{q_2}$) folgende Gleichung:

$0=(\dot{\overrightarrow{p}}-{\overrightarrow{F}}_G)\cdot \delta \overrightarrow{q}=(m{\ddot{q}}_1{\hat{e}}_{q_1}+m{\ddot{q}}_2{\hat{e}}_{q_2}+mg{\hat{e}}_2)\cdot \delta q_1{\hat{e}}_{q_1}=(m{\ddot{q}}_1+mg\sin \phi )\delta q_1$.

Unter der sinnvollen Forderung, dass $\delta q_1$ im Allg. nicht gleich Null ist, schließen wir hieraus:

$0=m{\ddot{q}}_1+mg\sin \phi$.

Die Bewegungsgleichung sieht übrigens auf dem ersten Blick so aus, als hätten wir darin einen Vorzeichenfehler: Man beachte aber, dass im Gegensatz zu der in der übrigen Literatur oft gewählten Koordinate unsere q₁-Achse die schiefe Ebene hinauf, d.h. sozusagen entgegen der Fahrtrichtung, zeigt.

Das D'Alembert'sche Prinzip, formuliert im System der alten Koordinaten $\overrightarrow{x}$, führt unter der Nebenbedingung $0=-\delta x_1\sin \phi +\delta x_2\cos \phi$ auf:

$0=(\dot{\overrightarrow{p}}-{\overrightarrow{F}}_G)\cdot \delta \overrightarrow{x}=m{\ddot{x}}_1\delta x_1+(m{\ddot{x}}_2+mg)\delta x_2=[m{\ddot{x}}_1+(m{\ddot{x}}_2+mg)\frac{\sin \phi }{\cos \phi }]\delta x_1$,

woraus wir

$0=m{\ddot{x}}_1\cos \phi +(m{\ddot{x}}_2+mg)\sin \phi =m{\ddot{q}}_1+mg\sin \phi$

erhalten, wobei wir darin $q_1=x_1\cos \phi +x_2\sin \phi$ ausgenutzt haben.

Als Alternative hierzu hätten wir auch die Nebenbedingung mit Hilfe sog. »Lagrange'scher Multiplikatoren« berücksichtigen können. Hierzu wird jede Nebenbedingung mit einem beliebigen Multiplikator versehen und von der D'Alembert'schen Gleichung abgezogen. In diesem Falle heiße dieser $\lambda$:

$0=(\dot{\overrightarrow{p}}-{\overrightarrow{F}}_G)\cdot \delta \overrightarrow{q}-\lambda \delta q_2=(m{\ddot{q}}_1{\hat{e}}_{q_1}+m{\ddot{q}}_2{\hat{e}}_{q_2}+mg{\hat{e}}_2)\cdot (\delta q_1{\hat{e}}_{q_1}+\delta q_2{\hat{e}}_{q_2})-\lambda \delta q_2$. $=(m{\ddot{q}}_1+mg\sin \phi )\delta q_1+(m{\ddot{q}}_2+mg\cos \phi -\lambda )\delta q_2$.

An der Gültigkeit dieser Gleichung ändert der zusätzliche Summand nichts, da er ja gleich Null ist. Weil $\lambda$ zunächst als beliebig angenommen wurde, können wir ihn auch so wählen, dass der Koeffizient von $\delta q_2$ verschwindet. Somit muss auch der Summand mit $\delta q_1$ gleich Null sein. Zudem wissen wir ja bereits, dass aufgrund der Nebenbedingung zwar $\delta q_2$ verschwindet, aber auch, dass dies im Allg. für $\delta q_1$ nicht gilt. Daher muss der Koeffizient von $\delta q_1$ gleich Null sein: $0=m{\ddot{q}}_1+mg\sin \phi$. Nicht immer sind die Verhältnisse so einfach wie hier: Der frei wählbare Parameter $\lambda$ trat nur in einem der beiden Summanden der Gleichung auf, was die Schlussfolgerung enorm erleichtert hat. Aber selbst dann, wenn dies nicht der Fall gewesen wäre, hätten wir durch das Einführen des zusätzlichen Parameters $\lambda$ sowieso davon ausgehen dürfen, dass aus $0=(a_1-\lambda A_1)\delta q_1+(a_2-\lambda A_2)\delta q_2$ $a_1-\lambda =a_2-\lambda =0$ folgt, wenn allgemein formuliert $0=a_1\delta q_1+a_2\delta q_2$ die D'Alembert'sche Gleichung und $0=A_1\delta q_1+A_2\delta q_2$ die Nebenbedingung darstellen. Denn wir können die letzteren beiden Gleichungen als ein lineares Gleichungssystem in Matrixform darstellen:

$0=\underset{\_}{A}\cdot \delta \overrightarrow{x}=\left(\begin{array}{cc}{a}_1& a_2\\ A_1& A_2\end{array}\right)\left(\begin{array}{c}\delta q_1\\ \delta q_2\end{array}\right)=\left(\begin{array}{c}{a}_1\\ A_1\end{array}\right)\delta q_1+\left(\begin{array}{c}{a}_2\\ A_2\end{array}\right)\delta q_2={\overrightarrow{A}}_1\delta q_1+{\overrightarrow{A}}_2\delta q_2$.

Es gilt, dass der Vektor $\delta \tilde{\overrightarrow{q}}={(\delta q_1,\delta q_2)}^T\ne \overrightarrow{0}$ ist (d.h. dass wenigstens eine seiner beiden Komponenten nicht verschwindet, was ja eine sinnvolle und notwendige Voraussetzung ist, damit wir überhaupt so etwas wie eine Nebenbedingung haben), wenn die Spaltenvektoren ${\overrightarrow{A}}_1$ und ${\overrightarrow{A}}_2$ der Matrix $\underset{\_}{A}$ linear abhängig sind. Die Determinante der Matrix $\underset{\_}{A}$ muss dann aber verschwinden: $det\left(\underset{\_}{A}\right)=0$. Bei Determinanten sind zudem einige Operationen möglich, ohne ihren Wert zu verändern, z.B. gilt $det\left(\underset{\_}{A}\right)=det\left({\underset{\_}{A}}^T\right)$, wobei ${\underset{\_}{A}}^T$ die Transposition von $\underset{\_}{A}$ ist, d.h. in ihr alle Zeilen und Spalten vertauscht sind. Wir können auch einer Spalte das Vielfache, z.B. $-\lambda$-fache, einer weiteren Spalte in der Determinate hinzu addieren, ohne den Wert der ursprünglichen Determinante zu verändern. Dies haben wir genau in dieser Reihenfolge an unserer Determinante $det\left(\underset{\_}{A}\right)$ so durchgeführt:

$0=det\left(\underset{\_}{A}\right)=\left|\begin{array}{cc}{a}_1& a_2\\ A_1& A_2\end{array}\right|=\left|\begin{array}{cc}{a}_1& A_1\\ a_2& A_2\end{array}\right|=\left|\begin{array}{cc}(a_1-\lambda A_1)& A_1\\ (a_2-\lambda A_2)& A_2\end{array}\right|=a_1{A}_2-a_2{A}_1$.

Der Parameter $\lambda$ ist aber eine beliebige Zahl, die wir auch so wählen können, dass z.B. $a_2-\lambda A_2=0$ gilt (zumindest, wenn bei $a_2\ne 0$ auch $A_2\ne 0$ ist). Aus

$0=det\left(\underset{\_}{A}\right)=\left|\begin{array}{cc}{a}_1-\lambda A_1& A_1\\ 0& A_2\end{array}\right|=A_2(a_1-\lambda A_1)$

folgt somit $a_1-\lambda A_1=0$ (da wir ja zuvor $A_2\ne 0$ vorausgesetzt haben).

Wir hätten also auch folgendermaßen argumentieren können: Mit der Forderung $a_2-\lambda A_2=0$ (da ja $\lambda$ frei wählbar ist, d.h. eine beliebige Zahl sein darf) können wir in $0=(a_1-\lambda A_1)\delta q_1+(a_2-\lambda A_2)\delta q_2$ direkt auf $0=(a_1-\lambda A_1)\delta q_1$ schließen. Da $\delta \tilde{\overrightarrow{q}}={(\delta q_1,\delta q_2)}^T\ne \overrightarrow{0}$ sein soll, muss mindestens eine seiner Komponente ungleich Null sein, z.B. $\delta x_1\ne 0$. Wegen $\delta x_1\ne 0$ können wir somit $a_1-\lambda A_1=0$ folgern. Wenn hingegen $\delta x_1=0$ gewesen wäre, dann hätten wir lieber $0=(a_1-\lambda A_1)$ gefordert. Da bei $\delta x_1=0$ nach Voraussetzung aber $\delta x_2\ne 0$ gelten muss, können wir aus $0=(a_2-\lambda A_2)\delta q_2$ auf $0=(a_2-\lambda A_2)$ schließen. D.h. durch eine entsprechende Umnumerierung $1\to 2$ bzw. $2\to 1$ der Indizes $i=1,2$ von $a_i,A_i,\delta q_i$ hätten wir auch bei der ersten Variante unserer Forderung verbleiben können.

Das D'Alembert'sche Prinzip, formuliert im System der alten Koordinaten $\overrightarrow{x}$, führt auf genau diese Situation:

$0=(\dot{\overrightarrow{p}}-{\overrightarrow{F}}_G)\cdot \delta \overrightarrow{x}-\lambda (-\delta x_1\sin \phi +\delta x_2\cos \phi )=(m{\ddot{x}}_1+\lambda \sin \phi )\delta x_1+(m{\ddot{x}}_2+mg-\lambda \cos \phi )\delta x_2$.

Hier dürfen wir also dennoch getrost fordern:

$0=m{\ddot{x}}_1+\lambda \sin \phi$, $0=m{\ddot{x}}_2+mg-\lambda \cos \phi$.

Multiplizieren wir die erstere Gleichung mit $\cos \phi$ und die Letztere mit $\sin \phi$ und addieren anschließend beide Gleichungen, dann erhalten wir erneut $0=m{\ddot{x}}_1\cos \phi +(m{\ddot{x}}_2+mg)\sin \phi =m{\ddot{q}}_1+mg\sin \phi$.

Wir können das D'Alembert'sche Prinzip statt mit der geometrisch motivierten Nebenbedingung $\delta q_2=0$ natürlich auch mit Hilfe der Zwangskraft $-{\overrightarrow{F}}_N=mg\cos \phi {\hat{e}}_{q_2}$ formulieren, die senkrecht auf der schiefen Ebene (und somit parallel zu ${\hat{e}}_{q_2}$) steht und verhindert, dass der betrachtete Körper durch Letztere einfach hindurch fällt:

$0=(\dot{\overrightarrow{p}}-({\overrightarrow{F}}_G-{\overrightarrow{F}}_N))\cdot \delta \overrightarrow{q}=(m{\ddot{q}}_1{\hat{e}}_{q_1}+m{\ddot{q}}_2{\hat{e}}_{q_2}+mg{\hat{e}}_2-mg\cos \phi {\hat{e}}_{q_2})\cdot (\delta q_1{\hat{e}}_{q_1}+\delta q_2{\hat{e}}_{q_2})$.

Wegen $mg{\hat{e}}_2\cdot (\delta q_1{\hat{e}}_{q_1}+\delta q_2{\hat{e}}_{q_2})=mg\sin \phi \delta q_1+mg\cos \phi \delta q_2$ erhalten wir daraus

$0=m{\ddot{q}}_1+mg\sin \phi$ und $0=m{\ddot{q}}_2$.

Kehren wir noch einmal zurück zu den Bewegungsgleichungen mit den Lagrange'schen Multiplikatoren, formuliert im System von $\overrightarrow{x}$. Multiplizieren wir $0=m{\ddot{x}}_1+\lambda \sin \phi$ dieses Mal mit $-\sin \phi$ und $0=m{\ddot{x}}_2+mg-\lambda \cos \phi$ mit $\cos \phi$ und addieren anschließend jene beiden Gleichungen, dann erhalten wir durch Berücksichtigen von $0=q_2=-x_1\sin \phi +x_2\cos \phi$ eine Aussage darüber, was der Langrange-Multiplikator $\lambda$ ist: $\lambda =mg\cos \phi$. Diesen Wert für $\lambda$ in die Gleichung $0=m{\ddot{q}}_2+mg\cos \phi -\lambda$ eingesetzt, ergibt natürlich wieder $0=m{\ddot{q}}_2$ (was wir wegen $0=q_2$ natürlich auch sofort erraten hätten). Wir erkennen aber auch an den letzten beiden Gleichungen, dass $\lambda {\hat{e}}_{q_2}=mg\cos \phi {\hat{e}}_{q_2}=-{\overrightarrow{F}}_N$ gleich der Zwangskraft ist. D.h. die Formulierung des D'Alembert'schen Prinzips mit Hilfe von Zwangskräften oder Lagrange'schen Multiplikatoren ist völlig gleichwertig. Es ist oft jedoch einfacher, aus der Geometrie des Versuchsaufbaus Nebenbedingungen (wie in diesem Beispiel $0=q_2$) abzulesen als die auftretenden Zwangskräfte zu bestimmen.

Ein Beispiel, in dem mehr als nur eine Nebenbedingung berücksichtigt werden muss, stellt z.B. das sog. physikalische Pendel dar, das aus zwei nicht notwendigerweise gleichen Massen ${\displaystyle }$ bestehen soll, die an den Enden einer starren (und als masselos gedachten) Stange in der xy-Ebene befestigt sind. Dieses System kann sich nur in der xy-Ebene um eine Achse senkrecht zu dieser Ebene drehen (d.h. die Drehachse ist parallel zur z-Achse). Die Achse verlaufe dabei durch einen Punkt auf der Stange, dem Ursprung des gewählten Koordinatensystems, der im Allg. nicht der Schwerpunkt des Systems sei: Siehe folgende Fig. 2. Auf jede Masse $m_i$ wirke eine Gravitationskraft entlang der y-Achse:

${\overrightarrow{F}}_{G,i}=-m_{i}g{\hat{e}}_y$

mit i=1,2.

Datei:Phys pendel.jpg Fig. 2: Physikalisches Pendel aus zwei Massen bestehend.

Die Masse $m_1$ am Ort ${\overrightarrow{x}}_1$ befinde sich also immer im festen Abstand vom Koordinatenursprung:

${\overrightarrow{x}}_1^2=const.⟹0={\overrightarrow{x}}_1\cdot \delta {\overrightarrow{x}}_1$.

Letztere Gleichung lässt sich durch eine infinitesimale Drehung um den Winkel $\delta {\overrightarrow{\phi }}_1$ erfüllen:

$\delta {\overrightarrow{x}}_1=\delta {\overrightarrow{\phi }}_1\times {\overrightarrow{x}}_1$.

Entsprechendes gilt für die Masse $m_2$ am Ort ${\overrightarrow{x}}_2$:

${\overrightarrow{x}}_2^2=const.⟹0={\overrightarrow{x}}_2\cdot \delta {\overrightarrow{x}}_2$.

Auch diese Gleichung lässt sich durch eine infinitesimale Drehung um den Winkel $\delta {\overrightarrow{\phi }}_2$ erfüllen:

$\delta {\overrightarrow{x}}_2=\delta {\overrightarrow{\phi }}_2\times {\overrightarrow{x}}_2$.

Die Drehungen erfolgen um die z- bzw. 3-Achse, d.h.

$\delta \overrightarrow{{\phi }_i}=\delta {\phi }_i{\hat{e}}_z$.

Daher erhalten wir aus $\delta x_i{\hat{e}}_x+\delta y_i{\hat{e}}_y=\delta {\overrightarrow{x}}_i=\delta \overrightarrow{\phi }\times {\overrightarrow{x}}_i=x_i\delta {\phi }_i{\hat{e}}_y-y_i\delta {\phi }_i{\hat{e}}_x$ die beiden Nebenbedingungen $\delta x_i=-y_i\delta {\phi }_i$ und $\delta y_i=x_i\delta {\phi }_i$.

Der Tatsache, dass beide Massen während der Drehbewegung fest an den Enden einer starren Stange sitzen (und die Rotation um einen Achse senkrecht zu Letzterer erfolgt), können wir durch die Forderung

$0=\delta {\phi }_1-\delta {\phi }_2$

Rechnung tragen, da die Drehwinkel bei beiden Massen immer gleich sein müssen.

Für unser 2-Teilchen-System in einer Ebene, d.h. Teilchenzahl N=2 und die Anzahl der Dimensionen D=2, haben wir somit drei Nebenbedingungen ($\nu =3$), woraus sich für die Anzahl der Freiheitsgrade wieder $f=D\cdot N-\nu =1$ ergibt: Der einzige Parameter, den wir benötigen, um die Bewegung des Systems zu beschreiben, ist ja der Drehwinkel $\phi$.

Diese Nebenbedingungen lassen sich in der Form $\delta x_i=-y_i\delta \phi$ bzw. $\delta y_i=x_i\delta \phi$ beispielsweise wieder direkt in die D'Alembert'sche Gleichung einsetzen:

$0=\underset{i=1}{\sum ^2}(\dot{{\overrightarrow{p}}_i}-{\overrightarrow{F}}_{G,i})\cdot \delta {\overrightarrow{x}}_i=\underset{i=1}{\sum ^2}[m_i({\ddot{x}}_i{\hat{e}}_x+{\ddot{y}}_i{\hat{e}}_y)+m_{i}g{\hat{e}}_y]\cdot (\delta x_i{\hat{e}}_x+\delta y_i{\hat{e}}_y)=$ $\underset{i=1}{\sum ^2}[m_i{\ddot{x}}_i\delta x_i+(m_i{\ddot{y}}_i+m_{i}g)\delta y_i]=\underset{i=1}{\sum ^2}[m_i(x_i{\ddot{y}}_i-y_i{\ddot{x}}_i)+g{m}_i{x}_i]\delta \phi$

Hierin stellt die Größe $\underset{i=1}{\sum ^2}{m}_i(x_i{\ddot{y}}_i-y_i{\ddot{x}}_i)=\frac{d}{dt}\underset{i=1}{\sum ^2}{m}_i(x_i{\dot{y}}_i-y_i{\dot{x}}_i)=\frac{d}{dt}{L}_z=M_z$ die z-Komponente des gesamten Drehmomentes $\overrightarrow{M}=\underset{i=1}{\sum ^2}{\overrightarrow{x}}_i\times m_i{\ddot{\overrightarrow{x}}}_i$ dar, der wiederum gleich der Zeitableitung des Gesamtdrehimpulses $\overrightarrow{L}=\underset{i=1}{\sum ^2}{\overrightarrow{x}}_i\times m_i{\dot{\overrightarrow{x}}}_i$ ist. Hier wird die Vielseitigkeit des D'Alembert'schen Prinzips erneut deutlich: Kräfte und Drehmomente treten dort nämlich gleichberechtigt auf. Denn wir sind von einem Term des Typs $\dot{\overrightarrow{p}}\cdot \delta \overrightarrow{x}$ ausgegangen und schließlich zu einem Term $\overrightarrow{M}\cdot \delta \overrightarrow{\phi }$ gelangt:

$\dot{\overrightarrow{p}}\cdot \delta \overrightarrow{x}=\dot{\overrightarrow{p}}\cdot (\delta \overrightarrow{\phi }\times \overrightarrow{x})=(\overrightarrow{x}\times \dot{\overrightarrow{p}})\cdot \delta \overrightarrow{\phi }=\overrightarrow{M}\cdot \delta \overrightarrow{\phi }$ .

Beiden Termen ist aber gemeinsam, dass sie die Dimension einer Energie besitzen.

Die z-Komponente des Drehimpulses lässt sich über den Parameter $\phi$ z.B. mittels Polarkoordinaten ausdrücken: $x_i=r_i\cos \phi$ und $y_i=r_i\sin \phi$, worin $r_i=\left|{\overrightarrow{x}}_i\right|$ bedeuten. Die zugehörigen Geschwindigkeiten lauten dann ${\dot{x}}_i=-\dot{\phi }{r}_i\sin \phi$ und ${\dot{y}}_i=\dot{\phi }{r}_i\cos \phi$, woraus sich

$L_z=\underset{i=1}{\sum ^2}{m}_i(x_i{\dot{y}}_i-y_i{\dot{x}}_i)=\underset{i=1}{\sum ^2}{m}_i{r}_i^2\dot{\phi }$

ergibt. Der Term $\underset{i=1}{\sum ^2}{m}_i{r}_i^2$ ist hierbei das Trägheitsmoment ${\mathrm{\Theta }}_{33}$.

Den Term $\underset{i=1}{\sum ^2}{m}_i{x}_i=\underset{i=1}{\sum ^2}{m}_i{r}_i\cos \phi$ gibt im Wesentlichen die x-Komponente des Schwerpunktes an. In der Literatur erscheint an dieser Stelle meist ein Sinus anstelle der Kosinus-Funktion hier. Dies hängt davon ab, wie der Winkel gemessen wird. Drücken wir unseren Winkel $\phi$ durch einen Winkel $\alpha$ aus, der zwischen der negativen y-Achse und der Stange gemessen wird, dann erhalten wir folgende Zusammenhänge: $\phi =\frac{3}{4}\pi +\alpha$ bzw. $\cos \phi =\sin \alpha$. In dieser Koordinate erhalten wir also folgende Bewegungsgleichung für unser physikalisches Pendel:

${\mathrm{\Theta }}_{33}\ddot{\alpha }+g\underset{i=1}{\sum ^2}{m}_i{r}_i\sin \alpha =0$.

Das D'Alembert'sche Prinzip, mittels Lagrange-Multiplikatoren formuliert, lautet:

$0=\underset{i=1}{\sum ^2}[(\dot{{\overrightarrow{p}}_i}-{\overrightarrow{F}}_{G,i})\cdot \delta {\overrightarrow{x}}_i+{\lambda }_i(x_i\delta x_i+y_i\delta y_i)]+{\lambda }_3(\delta {\phi }_1-\delta {\phi }_2)=$ $\underset{i=1}{\sum ^2}[m_i(x_i{\ddot{y}}_i-y_i{\ddot{x}}_i)+g{m}_i{x}_i+{\lambda }_i(y_i{x}_i-x_i{y}_i)]\delta {\phi }_i+{\lambda }_3(\delta {\phi }_1-\delta {\phi }_2)$

Für jede der drei Nebenbedingungen haben wir einen Lagrange-Multiplikator ${\lambda }_i$ eingeführt, wobei in obiger Gleichung die Koeffizienten von ${\lambda }_1$ und ${\lambda }_2$ aber offensichtlich verschwinden. Übrig bleibt nur der Parameter ${\lambda }_3$, was zu zwei Gleichungen führt:

$m_1(x_1{\ddot{y}}_1-y_1{\ddot{x}}_1)+g{m}_1{x}_1+{\lambda }_3=0$, $m_2(x_2{\ddot{y}}_2-y_2{\ddot{x}}_2)+g{m}_2{x}_2-{\lambda }_3=0$.

Beim Bilden der Summe aus beiden Gleichungen wird ${\lambda }_3$ eliminiert, woraus wir wieder die Gleichung

$\underset{i=1}{\sum ^2}[m_i(x_i{\ddot{y}}_i-y_i{\ddot{x}}_i)+g{m}_i{x}_i]=0$

erhalten. Der Lagrange'sche Multiplikator ${\lambda }_3$ besitzt hier offensichtlich die Bedeutung einer Zwangskraft im verallgemeinerten Sinne und drückt das gesamte Drehmoment auf ein Teilchen aus, das durch das jeweils andere Teilchen (übertragen durch die Stange) auf ihn einwirkt. Auch hier ist die »Kraft« in Wirklichkeit ein Drehmoment.

Zusammenfassend können wir feststelen, dass das D'Alembert'sche Prinzip insbesondere dann nützlich wird, wenn wir Nebenbedingungen an das betrachtete System stellen. Werden die Nebenbedingungen mit Hilfe von Lagrange'schen Multiplikatoren berücksichtigt, ermöglicht dies ein Aufstellen von Kräfte- oder Drehmoment-Bilanz-Gleichungen wie in der Newton'schen Dynamik, da die Terme mit jenen Multiplikatoren letztlich auch nichts andres als Zwangskräfte (oder -drehmomente) darstellen, die durch die Nebenbedingungen verursacht werden. Im Folgenden wollen wir dies ganz allgemein an einem System mit z.B. $M=D\cdot N$ Koordinaten $q_i$ diskutieren, das $\nu <M$ Nebenbedingungen der Form

$0=\underset{i=1}{\sum ^M}{A}_i^{\left(\kappa \right)}\delta q_i,\kappa =1,...,\nu$

unterliegt. Dies sind $\nu$ Gleichungen mit M Unbekannten $\delta q_i$. Im Idealfall können $\nu$ der $\delta q_i$ mit Hilfe der restlichen $M-\nu$ $\delta q_i$ ausgedrückt werden: Das System besitzt dann $f=M-\nu$ Freiheitsgrade. Z.B. könnten die $\nu$ $\delta q_i,i=(f+1),...,M$ dank der $\nu$ Nebenbedingungen durch die $f=M-\nu$ Variablen $\delta q_j,j=1,...,f$ mit Hilfe von in den letzteren Variablen linearen Funktionen $F_i,i=(f+1),...,M$ ausgedrückt werden: $\delta q_i=F_i(\delta q_1,...,\delta q_f),i=(f+1),...,M$.

Die $\nu$ Größen $\delta q_i,i=(f+1),...,M$ sind daher linear abhängig, während die übrigen $f=M-\nu$ Variablen $\delta q_j,j=1,...,f$ voneinander linear unabhängig sind. Lineare Unabhängigkeit der $\delta q_j,j=1,...,f$ bedeutet ja, dass aus $0=\underset{i=1}{\sum ^f}{b}_i\delta q_i$ das Verschwinden der f Koeffizienten $b_i$ folgt: $b_1=...=b_f=0$.

Die $\nu$ Nebenbedingungen multiplizieren wir jetzt jeweils mit einem beliebigen Paramater ${\lambda }_{\kappa }$, dem Lagrange'schen Multiplikator, und addieren die daraus entstehenden Gleichungen:

$0=\underset{\kappa =1}{\sum ^{\nu }}{\lambda }_{\kappa }\underset{i=1}{\sum ^M}{A}_i^{\left(\kappa \right)}\delta q_i$.

Letztere Gleichung ziehen wir jetzt von der D'Alembert'sche Gleichung ab, die im Allgemeinen folgende Form besitzt:

$0=\underset{i=1}{\sum ^M}{a}_i\delta q_i$

mit $a_i={\dot{p}}_i-F_i$.

D.h. wir erhalten dadurch eine Gleichung

$0=\underset{i=1}{\sum ^M}(a_i-\underset{\kappa =1}{\sum ^{\nu }}{\lambda }_{\kappa }{A}_i^{\left(\kappa \right)})\delta q_i$,

die wir auch als zwei Summen aus Termen in den f linear unabhängigen $\delta q_j,j=1,...,f$ bzw. aus Termen in den $\nu =M-f$ linear abhängigen $\delta q_i,i=(f+1),...,M$ darstellen können:

$0=\underset{j=1}{\sum ^f}(a_j-\underset{\kappa =1}{\sum ^{M-f}}{\lambda }_{\kappa }{A}_j^{\left(\kappa \right)})\delta q_j+\underset{i=f+1}{\sum ^M}(a_i-\underset{\kappa =1}{\sum ^{M-f}}{\lambda }_{\kappa }{A}_i^{\left(\kappa \right)})\delta q_i$.

Der zweite Summand in letzterer Gleichung enthält dann somit ausschließlich die $\nu$ linear abhängigen $\delta q_i,i=(f+1),...,M$. In deren Koeffizienten befinden sich die gleichermaßen $\nu$ beliebigen, d.h. noch frei wählbaren Lagrange-Multiplikatoren ${\lambda }_{\kappa }$. Wenn wir jene $\nu$ Koeffizienten jeweils gleich Null setzen, erhalten wir aus diesen ein lineares Gleichungssystem bestehend aus $\nu$ Gleichungen mit gleichermaßen $\nu$ Unbekannten ${\lambda }_{\kappa }$:

$0=a_i-\underset{\kappa =1}{\sum ^{M-f}}{\lambda }_{\kappa }{A}_i^{\left(\kappa \right)}$.

für $i=(f+1),...,M$.

Da dieses lineare Gleichungssystem (aufgrund des oben bereits angenommenen Idealfalls) erfüllt werden kann, ist es also möglich, durch eine entsprechende Wahl der $\nu$ beliebigen Parameter ${\lambda }_{\kappa }$ jene ebenfalls $\nu$ Koeffizienten zum Verschwinden zubringen. Es bleiben dann nur noch die Summanden in den $f=M-\nu$ voneinander linear unabhängigen $\delta q_j,j=1,...,f$ übrig:

$0=\underset{j=1}{\sum ^f}(a_j-\underset{\kappa =1}{\sum ^{M-f}}{\lambda }_{\kappa }{A}_j^{\left(\kappa \right)})\delta q_j$.

Aber gearde wegen der linearen Unabhängigkeit dieser f $\delta q_j$ dürfen wir auf

$0=a_j-\underset{\kappa =1}{\sum ^{M-f}}{\lambda }_{\kappa }{A}_j^{\left(\kappa \right)},j=1,...,f$

schließen. Insgesamt erhalten wir also M Gleichungen (die sog. Langrange'schen Gleichungen 1. Art)

$0=({\dot{p}}_l-F_l)-\underset{\kappa =1}{\sum ^{M-f}}{\lambda }_{\kappa }{A}_l^{\left(\kappa \right)},l=1,...,M$

aus der D'Alembert'schen Gleichung

$0=\underset{l=1}{\sum ^M}[({\dot{p}}_l-F_l)-\underset{\kappa =1}{\sum ^{M-f}}{\lambda }_{\kappa }{A}_l^{\left(\kappa \right)}]\delta q_l$.

In den vorangegangenen Beispielen sind wir bei den Nebenbedingungen von Gleichungen vom Typ

$0={\mathrm{\Phi }}_{\kappa }(q_1,...,q_M,t),\kappa =1,...,(M-f)$

ausgegangen, aus denen wir durch Anwenden der Variation $\delta$

$0=\underset{l=1}{\sum ^M}\frac{\partial {\mathrm{\Phi }}_{\kappa }}{\partial q_l}\delta q_l$

gewonnen haben (wobei in unseren Beispielen ${\mathrm{\Phi }}_{\kappa }$ nicht explizit von der Zeit abhing): d.h. $A_l^{\left(\kappa \right)}=\frac{\partial {\mathrm{\Phi }}_{\kappa }}{\partial q_l}$. Solche Nebenbedingungen werden »holonom« genannt. Denkbar wären aber z.B. auch solche, die sich nicht auf ein totales Differenzial zurückführen lassen (die dann »anholonom« heißen), oder solche, die Ungleichungen oder Differenzialgleichungen darstellen.

Die Nebenbedingungen lassen sich selbst für ein N-Teilchensystem (in D Dimensionen) auch bequem in vektorieller Form schreiben:

$0={\mathrm{\Phi }}_{\kappa }({\overrightarrow{q}}_1,...,{\overrightarrow{q}}_N),\kappa =1,...,(DN-f)$,

d.h.

$0=\underset{l=1}{\sum ^N}{\overrightarrow{\mathrm{\nabla }}}_{{\overrightarrow{q}}_l}{\mathrm{\Phi }}_{\kappa }\cdot \delta \overrightarrow{q_l},\kappa =1,...,(DN-f)$

mit dem Gradienten ${\overrightarrow{\mathrm{\nabla }}}_{{\overrightarrow{q}}_j}=\frac{\partial }{\partial {\overrightarrow{q}}_j}$.

Die Lagrange-Gleichungen 1. Art lauten dann unter diesen Nebenbedingungen folgendermaßen:

$0={\dot{\overrightarrow{p}}}_j-{\overrightarrow{F}}_j-\underset{\kappa =1}{\sum ^{DN-f}}{\lambda }_{\kappa }{\overrightarrow{\mathrm{\nabla }}}_{{\overrightarrow{q}}_j}{\mathrm{\Phi }}_{\kappa }({\overrightarrow{q}}_1,...,{\overrightarrow{q}}_N),j=1,...,N$

mit dem Impuls ${\overrightarrow{p}}_j=m_j{\dot{\overrightarrow{q}}}_j$ des j-ten Teilchens und der Kraft ${\overrightarrow{F}}_j$ auf das j-te Teilchen.

Wir können uns jetzt auch dem Problem widmen, dass wir in

$0=\underset{i=1}{\sum ^M}[(\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_i}-\frac{\partial L}{\partial q_i})-\underset{j=1}{\sum ^{D\cdot N}}{f}_j\frac{\partial x_j}{\partial q_i}]\delta q_i$.

nicht mehr darauf schließen können, dass die Koeffizienten der virtuellen Verrückungen $\delta q_i$ alle einzeln verschwinden, indem wir vom D'Alembert'schen Prinzip ausgehen, in dem die an das betrachtete System gestellten Nebenbedingungen berücksichtigt sind:

$0=\underset{l=1}{\sum ^M}[({\dot{p}}_l-F_l)-\underset{\kappa =1}{\sum ^{M-f}}{\lambda }_{\kappa }{A}_l^{\left(\kappa \right)}]\delta q_l$.

Wir können hier alle Überlegungen übernehmen, die wir schon für den eindimensionalen Fall angestellt haben, und mit denen wir von der D'Alembert'schen Gleichung zu den Euler-Lagrange-Gleichungen gelangt sind: Der zusätzliche Term, mit dem den Nebenbedingungen Rechnung getragen wird, stört hierbei nicht. D.h. es ergibt sich unmittelbar:

$0=\underset{i=1}{\sum ^{D\cdot N}}[(\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_i}-\frac{\partial L}{\partial q_i})-\underset{j=1}{\sum ^{D\cdot N}}{f}_j\frac{\partial x_j}{\partial q_i}-\underset{\kappa =1}{\sum ^{D\cdot N-f}}{\lambda }_{\kappa }{A}_i^{\left(\kappa \right)}]\delta q_i$.

Weil wir hierin die Nebenbedingungen ja berücksichtigt haben, dürfen wir wieder davon ausgehen, dass alle Koeffizienten der $\delta q_i$ verschwinden. Hieraus ergeben sich die sog. Langrange'schen Gleichungen 2. Art:

$0=(\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_i}-\frac{\partial L}{\partial q_i})-\underset{j=1}{\sum ^{D\cdot N}}{f}_j\frac{\partial x_j}{\partial q_i}-\underset{\kappa =1}{\sum ^{M-f}}{\lambda }_{\kappa }{A}_i^{\left(\kappa \right)},i=1,...,D\cdot N$.

In unserem Beispiel der schiefen Ebene lässt sich die natürliche Lagrangefunktion als $L=T-V$ sehr leicht angeben: Im Koordinatensystem von $\overrightarrow{x}$ erhalten wir die kinetische Energie

$T=\frac{1}{2}m({\dot{x}}_1^2+{\dot{x}}_2^2)$,

die potentielle Energie

$V=mg{x}_2$

sowie für die Nebenbedingung $A_1^{\left(1\right)}=-\sin \phi$ und $A_2^{\left(1\right)}=\cos \phi$.

Im Koordinatensystem von $\overrightarrow{q}$ sind dies

$T=\frac{1}{2}m({\dot{q}}_1^2+{\dot{q}}_2^2)$, $V=mg{x}_2=mg(q_1\sin \phi +q_2\cos \phi )$

und für die Nebenbedingung $A_1^{\left(1\right)}=0$ und $A_2^{\left(1\right)}=1$. In beiden Fällen (als auch im zweiten Beispiel) traten ja keine Nicht-Potenzialkräfte auf, d.h. $f_j=0$. Durch Einsetzen in die Lagrange'sche Gleichung 2. Art überzeuge man sich selbst davon, dass wieder die bereits bekannten Bewegungsgleichungen (d.h. die Lagrange'schen Gleichungen 1. Art) hieraus folgen. Im letzteren Fall hätte man natürlich wegen der Nebenbedingung $q_2=0$ auch direkt auf diese Koordinate und somit auf die Methode mit den Lagrange'schen Multiplikatoren verzichten können.

Im Beispiel des betrachteten physikalischen Pendels stellen wir zunächst fest, dass wegen $r_i^2={\overrightarrow{x}}_i^2=x_i^2+y_i^2=const.$ nicht nur $0={\overrightarrow{x}}_i\cdot \delta {\overrightarrow{x}}_i=x_i\delta x_i+y_i\delta y_i$ sondern auch $0={\overrightarrow{x}}_i\cdot {\dot{\overrightarrow{x}}}_i=x_i{\dot{x}}_i+y_i{\dot{y}}_i$ gilt, was ja (wie wir oben bereits gesehen haben) auf $\delta x_i=-y_i\delta {\phi }_i$ und $\delta y_i=x_i\delta {\phi }_i$ bzw. analog hierzu auf ${\dot{x}}_i=-y_i{\dot{\phi }}_i$ und ${\dot{y}}_i=x_i{\dot{\phi }}_i$ führt. Hieraus folgt für die kinetische Energie:

$T=\frac{1}{2}\underset{i=1}{\sum ^2}{m}_i{\dot{x}}_i^2=\frac{1}{2}\underset{i=1}{\sum ^2}{m}_i{r}_i^2{\dot{\phi }}_i^2$.

Die potentielle Energie ist wieder

$V=\underset{i=1}{\sum ^2}{m}_{i}g{y}_i\left({\phi }_i\right)$.

Wegen ${\dot{y}}_i=x_i{\dot{\phi }}_i$ gilt übrigens auch $\frac{\partial y_i}{\partial {\phi }_i}=x_i$, was wir beim Bilden von $\frac{\partial V}{\partial {\phi }_i}$ in den Euler-Lagrange-Gleichungen bzgl. der Variablen ${\phi }_i,i=1,2$ benötigen. Die Nebenbedingungen hinsichtlich dieser Variablen führen auf $A_1^{\left(1\right)}=A_2^{\left(2\right)}=0$ und $A_1^{\left(3\right)}=1$ sowie $A_2^{\left(3\right)}=-1$.

Natürlich hätten wir auch direkt ${\phi }_1={\phi }_2=\phi$ setzen und Polarkoordinaten $x_i=r_i\cos \phi$ bzw. $y_i=r_i\sin \phi$ einführen können. D.h. durch eine geschickte Wahl der neuen Koordinaten ist es auch hier wieder möglich, die an das System gestellten Nebenbedingungen direkt zu berücksichtigen, so dass das System nur noch durch genau so viele Variablen beschrieben wird, wie es Freiheitsgrade besitzt: Langrange'schen Multiplikatoren werden dann somit nicht benötigt. In solch einem Fall gilt also bei f Freiheitsgraden

$0=(\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_i}-\frac{\partial L}{\partial q_i})-\underset{j=1}{\sum ^{D\cdot N}}{f}_j\frac{\partial x_j}{\partial q_i},i=1,...,f$.

Dies wird aber voraussichtlich nicht für alle zu betrachtenden physikalischen Systeme immer so gelingen.