d’Alembertsches Prinzip

Das d’Alembertsche Prinzip (nach Jean-Baptiste le Rond d’Alembert) der klassischen Mechanik erlaubt die Aufstellung der Bewegungsgleichungen eines mechanischen Systems mit Zwangsbedingungen. Das Prinzip beruht auf dem Satz, dass die Zwangskräfte bzw. -momente in einem mechanischen System keine virtuelle Arbeit leisten.

Der Name „d’Alembertsches Prinzip“ wird von manchen Autoren für das Dynamische Gleichgewicht zwischen äußerer Kraft und d’Alembertscher Trägheitskraft verwendet, während andere Autoren dies mit heftigen Worten als eine unzulässige Verkürzung ablehnen.

Vorüberlegungen

Die Bewegungsgleichung für einen Massepunkt wird in einem Inertialsystem formuliert. Sie lautet nach dem zweiten newtonschen Gesetz:

$m{\vec a}={\vec F}$

Darin sind die Masse, ${\vec {a}}$ die Absolutbeschleunigung und ${\vec {F}}$ die äußere Kraft. Diese Grundgleichung der Mechanik kann auf die Form:

${\vec F}-m{\vec a}={\vec 0}$

gebracht werden. Der Term $-m{\vec a}$ wird als Kraft aufgefasst und als d'Alembertsche Trägheitskraft ${\vec {F}}_{T}$ bezeichnet.

${\vec F}+{\vec {F}}_{T}={\vec 0}$

Das dynamische Problem ist auf ein Gleichgewichtsproblem der Statik zurückgeführt. Man bezeichnet die Beziehung deshalb auch als dynamisches Gleichgewicht. Ein Problem der Dynamik kann somit auch mit Methoden der Statik behandelt werden, wenn Trägheitskräfte berücksichtigt werden. Beim d'Alembertschen Prinzip wird im Folgenden das Prinzip der virtuellen Arbeit ausgenutzt, das in der Statik zur Berechnung unbekannter Lagerkräfte eingesetzt werden kann.

Einführung

Bei einem System von N Massepunkten lautet die Bewegungsgleichung für die Masse i

$m_{i}{\ddot {{\vec r}}}_{i}={\vec F}_{i}$ .

Dabei ist ${\vec {F}}_{i}$ die resultierende äußere Kraft auf den Massepunkt i. Sie ist die Summe aus eingeprägter Kraft ${\vec {F_{i}^{e}}}$ und Zwangskraft ${\vec {F_{i}^{z}}}$ .

${\vec F}_{i}={\vec {F_{i}^{e}}}+{\vec {F_{i}^{z}}}\;.$

Eingesetzt in die Newtonsche Bewegungsgleichung:

$m_{i}{\ddot {\vec {r}}}_{i}={\vec {F_{i}^{e}}}+{\vec {F_{i}^{z}}}\;.$

Die Zwangskraft berechnet sich somit zu

${\vec {F}}_{i}^{z}=m_{i}{\ddot {\vec {r}}}_{i}-{\vec {F}}_{i}^{e}\;.$

Man bildet das Skalarprodukt der Zwangskräfte mit den virtuellen Verschiebungennote^{[Anmerkung 1]} $\delta {\vec {r}}_{i}$ . Wenn nach dem Prinzip der virtuellen Arbeit die Zwangskräfte insgesamt keine virtuelle Arbeit verrichten, verschwindet die Summe der Skalarprodukte von Zwangskräften und virtuellen Verschiebungen:

$\sum _{i=1}^{N}\left({\vec {F}}_{i}^{z}\cdot \delta {\vec {r}}_{i}\right)=0\;.$

Man erhält das d’Alembertsche Prinzip (in der Formulierung von Joseph-Louis Lagrange):

${\sum _{i=1}^{N}\left(m_{i}{\ddot {\vec {r}}}_{i}-{\vec {F}}_{i}^{e}\right)\cdot \delta {\vec {r}}_{i}=0}$

In der Gleichung treten die Zwangskräfte nicht mehr auf – nur die eingeprägten Kräfte. Die Zwangsbedingungen verstecken sich noch in den virtuellen Verschiebungen, denn es sind nur solche erlaubt, die mit den Zwangsbedingungen vereinbar sind.

Um daraus Bewegungsgleichungen zu gewinnen, geht man bei $\,k$ (holonomen) Zwangsbedingungen zu f=3N-k unabhängigen Koordinaten (Freiheitsgraden) $q=(q_{1}(t),...,q_{f}(t))$ über und drückt Lage, Geschwindigkeit, Beschleunigung, und virtuelle Verschiebungen der N Massen durch diese neuen Lagekoordinaten („generalisierte Koordinaten“) aus.

${\vec {r}}_{i}={\vec {r}}_{i}(q)\;,\quad {\dot {{\vec r}}}_{i}={\dot {{\vec r}}}_{i}(q,{\dot q})\;,\quad {\ddot {{\vec r}}}_{i}={\ddot {{\vec r}}}_{i}(q,{\dot q},{\ddot q})\;,\quad \delta {{\vec {r}}_{i}}=\sum _{{j=1}}^{f}{\frac {\partial {\vec {r}}_{i}}{\partial q_{j}}}\,\delta q_{j}\,.$

Da sich die neuen Koordinaten unabhängig variieren lassen, ergeben sich Differentialgleichungen zweiter Ordnung, die sich nach $\ddot{q}$ auflösen lassen. Die konkrete Vorgehensweise zur Aufstellung der Bewegungsgleichungen ist dem nächsten Abschnitt zu entnehmen.

Für holonome Zwangsbedingungen und konservative Kräfte (die sich aus einer Potentialfunktion ableiten lassen) ist das d’Alembert-Prinzip dann äquivalent zu den Lagrangegleichungen erster Art.

Gelegentlich wird schon die eingangs wiedergegebene einfache Umstellung der newtonschen Bewegungsgleichung als das d’Alembertsche Prinzip bezeichnet. Das übersieht aber wesentliche Folgerungen wie die Elimination von Zwangskräften, die keine virtuelle Arbeit leisten und kommt in den Worten von Georg Hamel fast einer Beleidigung von d’Alembert gleich. Es ist zudem zu beachten, dass das verwendete Prinzip der virtuellen Arbeit nicht aus den Newtonschen Axiomen folgt, sondern ein eigenes Grundpostulat darstellt.

Erweiterung auf Mehrkörpersysteme

Im allgemeinen Fall von Mehrkörpersystemen wird berücksichtigt, dass auch die virtuelle Arbeit der Zwangsmomente auf den virtuellen Verdrehungen verschwindet. Zur Berechnung der Zwangsmomente wird die Eulersche Gleichung verwendet.

${\sum _{{i=1}}^{N}\left(\left[m_{i}{\ddot {{\vec r}}}_{i}-{\vec {F}}_{i}^{e}\right]\delta {\vec {r}}_{i}^{{\,T}}+\left[I_{i}\,{\dot {{\vec \omega }}}_{i}+{\vec {\omega }}_{i}\times I_{i}\,{\vec {\omega }}_{i}-{\vec M}_{i}^{e}\right]\delta {\vec {\varphi }}_{i}^{{\,T}}\right)=0}.$

mit

$I_{i}$ Trägheitstensor des Körpers i

${\dot {{\vec \omega }}}_{i}$ Winkelbeschleunigung des Körpers i

${\vec \omega }_{i}$ Winkelgeschwindigkeit des Körpers i

${\vec M}_{i}^{e}$ eingeprägtes Moment auf den Körper i

$\delta {\vec {\varphi }}_{i}$ virtuelle Verdrehung des Körpers i.

Bei N Körpern und k Bindungen ergeben sich $f=6\,N-k$ Freiheitsgrade.

Die virtuellen Verschiebungen bzw. Verdrehungen erhält man aus den partiellen Ableitungen der translatorischen bzw. rotatorischen Lagekoordinaten nach den verallgemeinerten Koordinaten:

$\delta {{\vec {r}}_{i}}=\sum _{{j=1}}^{f}{\frac {\partial {\vec {r}}_{i}}{\partial q_{j}}}\,\delta q_{j}\,$

$\delta {{\vec {\varphi }}_{i}}=\sum _{{j=1}}^{f}{\frac {\partial {\vec {\varphi }}_{i}}{\partial q_{j}}}\,\delta q_{j}$

Die Beschleunigungen lassen sich in einen Teil, der nur von den zweiten Ableitungen der verallgemeinerten Koordinaten abhängt, und einen Restterm zerlegen:

${\ddot {{\vec r}}}_{i}=\sum _{{j=1}}^{f}{\frac {\partial {\vec {r}}_{i}}{\partial q_{j}}}\,{\ddot {q}}_{j}+{\vec {a}}_{i}^{{\,*}}$ und

${\ddot {{\vec {\alpha }}}}_{i}=\sum _{{j=1}}^{f}{\frac {\partial {\vec {\varphi }}_{i}}{\partial q_{j}}}\,{\ddot {q}}_{j}+{\vec {\alpha }}_{i}^{{\,*}}$ .

Damit lässt sich das Differentialgleichungssystem zweiter Ordnung in Matrixform darstellen.

${\mathbf {M}}\,{\ddot {{\vec {q}}}}={\vec {F}}^{*}+{\vec {M}}^{*}$

Dabei sind:

${\mathbf {M}}$ die f x f Massenmatrix

${\vec {F}}^{*}$ der Vektor der verallgemeinerten Kräfte

${\vec {M}}^{*}$ der Vektor der verallgemeinerten Momente

Die Elemente der Massenmatrix berechnen sich zu:

$M_{{m,n}}=\sum _{{i=1}}^{N}\left(m_{i}\,{\frac {\partial {\vec {r}}_{i}^{{\,T}}}{\partial q_{m}}}\cdot {\frac {\partial {\vec {r}}_{i}}{\partial q_{n}}}+{\frac {\partial {\vec {\varphi }}_{i}^{{\,T}}}{\partial q_{m}}}\cdot I_{i}\cdot {\frac {\partial {\vec {\varphi }}_{i}}{\partial q_{n}}}\right)$

Für die Komponenten verallgemeinerten Kräfte bzw. Momente ergibt sich:

$F_{m}^{*}=\sum _{{i=1}}^{N}\left({\frac {\partial {\vec {r}}_{i}^{{\,T}}}{\partial q_{m}}}\left[{\vec {F}}_{i}{^{e}}-m_{i}\,{\vec {a}}_{i}^{{\,*}}\right]\right)$

$M_{m}^{*}=\sum _{{i=1}}^{N}\left({\frac {\partial {\vec {\varphi }}_{i}^{{\,T}}}{\partial q_{m}}}\left[{\vec {M}}_{i}{^{e}}-I_{i}\,{\vec {\alpha }}_{i}^{{\,*}}-{\vec {\omega }}_{i}\times I_{i}\,{\vec {\omega }}_{i}\right]\right)$

Die Berechnung der Massenmatrix sowie der verallgemeinerten Kräfte und Momente kann numerisch im Rechner durchgeführt werden. Das Differentialgleichungssystem kann ebenfalls numerisch mit gängigen Programmen gelöst werden. Die Behandlung großer Mehrkörpersysteme mit kinematischen Bindungen wird so erst möglich.

Beispiel Fadenpendel

Fadenpendel:
$\varphi$ ist die Auslenkung aus der Gleichgewichtslage und generalisierte Koordinate

Beim ebenen Fadenpendel mit der Masse wird der Winkel $\varphi$ , mit dem der Faden aus der Ruheposition ausgelenkt ist, als Freiheitsgrad gewählt. Die konstante Fadenlänge stellt eine holonome Zwangsbedingung dar. Position, Geschwindigkeit und Beschleunigung der Masse können daher in Abhängigkeit dieses Winkels ausgedrückt werden:

${\vec {r}}={\begin{bmatrix}l\sin {\varphi }\\-l\cos {\varphi }\end{bmatrix}}$

${\dot {{\vec {r}}}}={\frac {\partial {\vec r}}{\partial \varphi }}\,{\dot {\varphi }}={\begin{bmatrix}l\cos {\varphi }\\l\sin {\varphi }\end{bmatrix}}{\dot {\varphi }}$

${\ddot {{\vec {r}}}}={\begin{bmatrix}l\cos {\varphi }\\l\sin {\varphi }\end{bmatrix}}{\ddot {\varphi }}+{\begin{bmatrix}-l\sin {\varphi }\\l\cos {\varphi }\end{bmatrix}}{\dot {\varphi }}^{2}$

Die virtuelle Verschiebung ergibt sich zu:

$\delta {\vec {r}}={\frac {\partial {\vec r}}{\partial \varphi }}\,\delta \varphi ={\begin{bmatrix}l\cos {\varphi }\\l\sin {\varphi }\end{bmatrix}}\delta \varphi$

Als eingeprägte Kraft wirkt die Gewichtskraft:

${\vec G}={\begin{bmatrix}0\\-m\,g\end{bmatrix}}$

Die Bewegungsgleichung ergibt sich aus der Bedingung, dass die virtuelle Arbeit der Zwangskräfte verschwindet.

$\left(m{\ddot {{\vec r}}}-{\vec {G}}\right)\delta {\vec {r}}=0\Rightarrow \left(m{\ddot {{\vec r}}}-{\vec {G}}\right){\frac {\partial {\vec r}}{\partial \varphi }}=0.$

Durch Auswertung der Skalarprodukte erhält man schließlich:

$m\,l^{2}\,{\ddot \varphi }=-m\,g\,l\,\sin \varphi$

Masse und Fadenlänge lassen sich kürzen, so dass man die bekannte Differentialgleichung:

${\ddot \varphi }=-{\frac {g}{l}}\,\sin \varphi$

erhält.

Die Vorgehensweise erscheint bei diesem einfachen Beispiel sehr umständlich. Da aber nur Skalarprodukte ausgewertet werden müssen, kann dies bei großen Systemen automatisiert werden und numerisch im Rechner durchgeführt werden. Dies erleichtert die Aufstellung von Bewegungsgleichungen wesentlich.

Anmerkungen

↑ Infinitesimale Verschiebungen heißen virtuell, wenn sie mit den Zwangsbedingungen verträglich sind. Außerdem sollen sie unmittelbar (oder instantan, zu einer festen Zeit) erfolgen.

Basierend auf einem Artikel in:

Wikipedia.de