d’Alembertsches Prinzip

aus Wikipedia, der freien Enzyklopädie
Wechseln zu: Navigation, Suche

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.[1][2][3]

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

Vorüberlegungen[Bearbeiten]

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 m 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

Die Gleichung ist formal von einem statischen Gleichgewicht nicht zu unterscheiden. Man bezeichnet sie 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[Bearbeiten]

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 Verschiebungen[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 Lagrange):[6][7][1]

 {\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 = 3 N - 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 f 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.[8] 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[Bearbeiten]

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[Bearbeiten]

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

Beim ebenen Fadenpendel mit der Masse m wird der Winkel \varphi, mit dem der Faden aus der Ruheposition ausgelenkt ist, als Freiheitsgrad gewählt. Die konstante Fadenlänge l 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.

Literatur[Bearbeiten]

  • Herbert Goldstein, Charles P. Poole, John L. Safko: Klassische Mechanik. VCH.
  • Friedhelm Kuypers: Klassische Mechanik. VCH, 5. Auflage 1997, ISBN 3-527-29269-1.
  • Georg Hamel: Theoretische Mechanik. Springer 1967.
  • Werner Schiehlen: Technische Dynamik. Teubner Studienbücher, Stuttgart, 1986.

Anmerkungen[Bearbeiten]

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

Einzelnachweise[Bearbeiten]

  1. a b  Jürgen Dankert, Helga Dankert: Technische Mechanik: Statik, Festigkeitslehre, Kinematik/Kinetik. 7. Auflage. Springer Vieweg, 2013, ISBN 978-3-8348-2235-2 (eingeschränkte Vorschau in der Google-Buchsuche).
  2. Klaus-Peter Schnelle: Simulationsmodelle für die Fahrdynamik von Personenkraftwagen unter Berücksichtigung der nichtlinearen Fahrwerkskinematik. VDI-Verlag, Düsseldorf 1990, ISBN 3-18-144612-2. (Fortschrittsberichte VDI Nr. 146), S. 73
  3. Skript TU Berlin, PDF 120 KB
  4. Hans J. Paus: Physik in Experimenten und Beispielen. Hanser 2007, S. 34.
  5. Istvan Szabo: Geschichte der Mechanischen Prinzipien. Springer-Verlag, 1987, S. 40.
  6. István Szabó: Geschichte Der Mechanischen Prinzipien Und Ihrer Wichtigsten Anwendungen.. Springer DE, 1987, ISBN 978-3-7643-1735-5, S. 39– (Zugriff am 8. Februar 2013).
  7.  Kurt Magnus, H. H. Müller-Slany: Grundlagen der Technische Mechanik. 7. Auflage. Vieweg+Teubner, 2005, ISBN 978-3-8351-0007-7, S. 258 (eingeschränkte Vorschau in der Google-Buchsuche).
  8. Hamel Theoretische Mechanik, Springer 1967, S.220