Drehimpuls und Bewegungsgleichungen: Unterschied zwischen den Versionen

Aus PhysikWiki
Zur Navigation springen Zur Suche springen
Die Seite wurde neu angelegt: „<noinclude>{{Scripthinweis|Mechanik|6|3}}</noinclude> Betrachten wir eine infinitesimale Verrückung <math>\begin{align} & d\bar{r}=d{{{\bar{r}}}_{S}}+d\bar{x}…“
 
*>SchuBot
K Interpunktion, replaced: ! → ! (2), ( → ( (2)
 
(3 dazwischenliegende Versionen von 2 Benutzern werden nicht angezeigt)
Zeile 1: Zeile 1:
<noinclude>{{Scripthinweis|Mechanik|6|3}}</noinclude>
<noinclude>{{Scripthinweis|Mechanik|6|3}}</noinclude>


Betrachten wir eine infinitesimale Verrückung
====Drehimpuls====
<math>\begin{align}
# diskret:
  & d\bar{r}=d{{{\bar{r}}}_{S}}+d\bar{x}=d{{{\bar{r}}}_{S}}+d\bar{\phi }\times \bar{x} \\
& d\bar{\phi }:=\bar{n}d\phi  \\
\end{align}</math>


:<math>\begin{align}
  & \bar{l}=\sum\limits_{i=1}^{n}{{}}{{m}_{i}}{{{\bar{r}}}_{i}}\times {{{\dot{\bar{r}}}}_{i}}=\sum\limits_{i=1}^{n}{{}}{{m}_{i}}\left( {{{\bar{r}}}_{S}}+{{{\bar{x}}}^{(i)}} \right)\times \left( \bar{V}+\bar{\omega }\times {{{\bar{x}}}^{(i)}} \right) \\
& \bar{l}=\sum\limits_{i=1}^{n}{{}}{{m}_{i}}\left( {{{\bar{r}}}_{S}}\times \bar{V} \right)+\sum\limits_{i=1}^{n}{{}}{{m}_{i}}{{{\bar{x}}}^{(i)}}\times \bar{V}+{{{\bar{r}}}_{S}}\times \left( \bar{\omega }\times \sum\limits_{i}{{{m}_{i}}}{{{\bar{x}}}^{(i)}} \right)+\sum\limits_{i}{{{m}_{i}}}{{{\bar{x}}}^{(i)}}\times \left( \bar{\omega }\times {{{\bar{x}}}^{(i)}} \right) \\
& \sum\limits_{i=1}^{n}{{}}{{m}_{i}}{{{\bar{x}}}^{(i)}}\times \bar{V}=\sum\limits_{i=1}^{n}{{}}{{m}_{i}}\left( {{{\bar{x}}}^{(i)}} \right)=0 \\
& \bar{l}=\sum\limits_{i=1}^{n}{{}}{{m}_{i}}\left( {{{\bar{r}}}_{S}}\times \bar{V} \right)+\sum\limits_{i}{{{m}_{i}}}{{{\bar{x}}}^{(i)}}\times \left( \bar{\omega }\times {{{\bar{x}}}^{(i)}} \right)=M\left( {{{\bar{r}}}_{S}}\times \bar{V} \right)+\sum\limits_{i}{{{m}_{i}}}{{{\bar{x}}}^{(i)}}\times \left( \bar{\omega }\times {{{\bar{x}}}^{(i)}} \right) \\
\end{align}</math> Mit <math>M\left( {{{\bar{r}}}_{S}}\times \bar{V} \right)</math>
als Schwerpunktsdrehimpuls


In Kapitel 3.3 haben wir bereits mit infinitesimalen Drehungen gearbeitet.  Dort handelte es sich um passive Drehungen. Hier haben wir es nun mit aktiven Drehungen zu tun -> anderes Vorzeichen.


:<math>\sum\limits_{i}{{{m}_{i}}}{{\bar{x}}^{(i)}}\times \left( \bar{\omega }\times {{{\bar{x}}}^{(i)}} \right)</math>
als Relativdrehimpuls


<math>\bar{V}:=\frac{d{{{\bar{r}}}_{s}}}{dt}</math> Schwerpunktsgeschwindigkeit
# kontinuierliche Situation




<math>\bar{\omega }:=\frac{d\bar{\phi }}{dt}</math> Winkelgeschwindigkeit
:<math>\begin{align}
  & \bar{l}={{{\bar{r}}}_{s}}\times M\bar{V}+\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x})\bar{x}}\times \left( \bar{\omega }\times \bar{x} \right) \\
& \bar{L}=\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x})\bar{x}}\times \left( \bar{\omega }\times \bar{x} \right)=\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x})\left[ {{x}^{2}}\bar{\omega }-\left( \bar{x}\cdot \bar{\omega } \right)\bar{x} \right]} \\
\end{align}</math>


Damit ergibt sich die Geschwindigkeit eines beliebigen Aufpunktes des starren Körpers:


Also:


<math>\bar{v}=\bar{V}+\bar{\omega }\times \bar{x}</math>


:<math>\bar{L}=\bar{\bar{J}}\bar{\omega }</math>


Nebenbemerkungen:


Dies sieht man an der Komponentenschreibweise:


<math>\bar{\omega }</math> hängt von der Wahl von S ab.


Falls S der Schwerpunkt ist, so gilt:
:<math>{{L}_{m}}=\sum\limits_{n=1}^{3}{{}}\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x})}\left[ {{x}^{2}}{{\delta }_{mn}}-{{x}_{m}}{{x}_{n}} \right]{{\omega }_{n}}=\sum\limits_{n=1}^{3}{{}}{{J}_{mn}}{{\omega }_{n}}</math>




<math>\sum\limits_{i=1}^{n}{{{m}_{i}}{{{\bar{x}}}^{(i)}}}=0</math>
Nebenbemerkung:
nach Def. A) des starren Körpers


Im Allgemeinen ist
:<math>\bar{L}</math>
nicht parallel zu
:<math>\bar{\omega }</math>,
nur falls
:<math>\bar{\omega }</math>
in Richtung der Hauptträgheitsachse liegt!


<math>\int_{{}}^{{}}{{{d}^{3}}x\bar{x}\rho (\bar{x})}=0</math>
====Allgemeine Bewegungsgleichung für den Gesamtdrehimpuls====
Definition B) -> Schwerpunktsvektor im körperfesten System
<math>\bar{K}</math>
:
 
====Kinetische Energie:====


#<math>\frac{1}{2}\sum\limits_{i=1}^{n}{{{m}_{i}}{{v}^{(i)}}^{2}}=\frac{1}{2}\sum\limits_{i=1}^{n}{{{m}_{i}}}{{\left( V+\omega \times {{x}^{(i)}} \right)}^{2}}=\frac{1}{2}\sum\limits_{i=1}^{n}{{{m}_{i}}}{{V}^{2}}+V\cdot \sum\limits_{i=1}^{n}{{{m}_{i}}}\left( \omega \times {{x}^{(i)}} \right)+\frac{1}{2}\sum\limits_{i=1}^{n}{{{m}_{i}}}{{\left( \omega \times {{x}^{(i)}} \right)}^{2}}</math>


Mit den Beziehungen
:<math>\frac{d}{dt}\bar{l}=\sum\limits_{i}{{{{\bar{r}}}_{i}}\times {{{\bar{F}}}_{i}}}</math>.
<math>\begin{align}
  Dabei sind
  & \sum\limits_{i=1}^{n}{{{m}_{i}}}=M \\
:<math>{{\bar{F}}_{i}}</math>
& \bar{V}\cdot \sum\limits_{i=1}^{n}{{{m}_{i}}}\left( \bar{\omega }\times {{{\bar{x}}}^{(i)}} \right)=\left( \bar{V}\times \bar{\omega } \right)\sum\limits_{i=1}^{n}{{{m}_{i}}}{{{\bar{x}}}^{(i)}}=0,da\sum\limits_{i=1}^{n}{{{m}_{i}}}{{{\bar{x}}}^{(i)}}=0 \\
äußere, eingeprägte Kräfte. Die resultierende Kraft
  & {{\left( \bar{\omega }\times {{{\bar{x}}}^{(i)}} \right)}^{2}}={{\omega }^{2}}{{x}^{2}}{{\sin }^{2}}\alpha ={{\omega }^{2}}{{x}^{2}}(1-{{\cos }^{2}}\alpha )={{\omega }^{2}}{{x}^{2}}-{{\left( \bar{\omega }\cdot \bar{x} \right)}^{2}}=\sum\limits_{m=1}^{3}{\sum\limits_{n=1}^{3}{{}}{{\omega }^{m}}\left[ {{x}^{2}}{{\delta }_{mn}}-{{x}_{m}}{{x}_{n}} \right]}{{\omega }^{n}} \\
:<math>\bar{F}({{\bar{r}}_{S}})=\sum\limits_{i}{{}}{{\bar{F}}_{i}}</math>
\end{align}</math>
soll auf den Schwerpunkt wirken, so dass gilt:
Somit folgt:
<math>T=\frac{1}{2}\sum\limits_{i=1}^{n}{{{m}_{i}}{{{v}}^{(i)}}^{2}}=\frac{1}{2}M{{V}^{2}}+\frac{1}{2}\sum\limits_{m=1}^{3}{\sum\limits_{n=1}^{3}{{}}{{\omega }^{m}}{{J}_{mn}}}{{\omega }^{n}}</math>
mit dem Trägheitstensor




<math>{{J}_{mn}}:=\sum\limits_{i=1}^{n}{{}}{{m}_{i}}\left[ {{{{x}}}^{(i)}}^{2}{{\delta }_{mn}}-{{x}_{m}}^{(i)}{{x}_{n}}^{(i)} \right]</math>
:<math>\sum\limits_{i}{{}}\frac{{{{\bar{F}}}_{i}}}{{{m}_{i}}}=\frac{\bar{F}({{{\bar{r}}}_{S}})}{M}</math>




Der Trägheitstensor ist also durch die Massenverteilung bestimmt
Somit:


Im Sinne der Definition B) dagegen gilt:


:<math>\frac{d}{dt}\bar{l}=\sum\limits_{i}{{{m}_{i}}{{{\bar{r}}}_{i}}\times \frac{{{{\bar{F}}}_{i}}}{{{m}_{i}}}={{{\bar{r}}}_{S}}\times \bar{F}({{{\bar{r}}}_{S}})}</math>


<math>\begin{align}
  & T=\frac{1}{2}\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x}){{\left( \bar{V}+\bar{\omega }\times \bar{x} \right)}^{2}}}=\frac{1}{2}\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x}){{V}^{2}}+\left( \bar{V}\times \bar{\omega } \right)\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x})\bar{x}+\frac{1}{2}\sum\limits_{m=1}^{3}{\sum\limits_{n=1}^{3}{{}}{{\omega }^{m}}{{J}_{mn}}}{{\omega }^{n}}}} \\
& mit\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x})\bar{x}=0} \\
\end{align}</math>


Bekanntlich gilt für die Schwerpunktsbewegung:


und dem Trägheitstensor


:<math>M{{\ddot{\bar{r}}}_{S}}=\bar{F}({{\bar{r}}_{S}})</math>
(Newton)


<math>{{J}_{mn}}=\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x})\left[ {{x}^{2}}{{\delta }_{mn}}-{{x}_{m}}{{x}_{n}} \right]}</math>


:<math>\frac{d}{dt}\bar{l}=\frac{d}{dt}\left( {{{\bar{r}}}_{s}}\times M\bar{V}+\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x})\bar{x}}\times \left( \bar{\omega }\times \bar{x} \right) \right)=\frac{d}{dt}\left( {{{\bar{r}}}_{s}}\times M{{{\dot{\bar{r}}}}_{s}} \right)+\frac{d}{dt}\bar{L}=M{{\dot{\bar{r}}}_{s}}\times {{\dot{\bar{r}}}_{s}}+{{\bar{r}}_{S}}\times M{{\ddot{\bar{r}}}_{S}}+\frac{d}{dt}\bar{L}</math>


Also gilt die Zerlegung der kinetischen Energie:




<math>\begin{align}
:<math>\begin{align}
   & T=\frac{1}{2}\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x}){{\left( \bar{V}+\bar{\omega }\times \bar{x} \right)}^{2}}}=\frac{1}{2}\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x}){{V}^{2}}+\frac{1}{2}\sum\limits_{m=1}^{3}{\sum\limits_{n=1}^{3}{{}}{{\omega }^{m}}{{J}_{mn}}}{{\omega }^{n}}} \\
   & M{{{\dot{\bar{r}}}}_{s}}\times {{{\dot{\bar{r}}}}_{s}}=0 \\
\\
& {{{\bar{r}}}_{S}}\times M{{{\ddot{\bar{r}}}}_{S}}={{{\bar{r}}}_{S}}\times \bar{F}({{{\bar{r}}}_{S}}) \\
\end{align}</math>
\end{align}</math>






<math>\begin{align}
:<math>\frac{d}{dt}\bar{l}={{\bar{r}}_{S}}\times \bar{F}({{\bar{r}}_{S}})+\frac{d}{dt}\bar{L}</math>
  & T=\frac{1}{2}M{{V}^{2}}+\frac{1}{2}\bar{\omega }\bar{\bar{J}}\bar{\omega } \\
&  \\
\end{align}</math>




Dabei ist
Gleichzeitig gilt:
:<math>\frac{d}{dt}\bar{l}=\sum\limits_{i}{{{m}_{i}}{{{\bar{r}}}_{i}}\times \frac{{{{\bar{F}}}_{i}}}{{{m}_{i}}}={{{\bar{r}}}_{S}}\times \bar{F}({{{\bar{r}}}_{S}})}</math>




<math>{{T}_{trans}}=\frac{1}{2}M{{V}^{2}}</math>
Somit:
kinetische Energie der translatorischen Bewegung




<math>{{T}_{rot}}=\frac{1}{2}\bar{\omega }\bar{\bar{J}}\bar{\omega }</math>
:<math>\frac{d}{dt}\bar{L}=0</math>
kinetische Energie der Rotationsbewegung
Der Relativdrehimpuls ist im Schwerpunktsystem mit raumfesten Achsen konstant.


====Eigenschaften des Trägheitstensors====
Also: Es verschwindet die Zeitableitung des relativdrehimpulses im Schwerpunktsystem mit RAUMFESTEN Achsen.


<math>\bar{\bar{J}}</math>
Die Transformation muss nun noch ins körperfeste, rotatorisch mitbewegte System
ist ein Tensor zweiter Stufe. Das heißt unter Drehungen
:<math>\bar{K}</math>
<math>R\in SO(3)</math>
erfolgen:
transformiert er sich wie folgt:


R kennzeichnet dabei die Drehmatrizen im
Dabei sieht der Beobachter im RAUMFESTEN System neben der zeitlichen Änderung
<math>{{R}^{3}}</math>
:<math>\left( \frac{d}{dt} \right)\acute{\ }</math>,
mit Orthogonalitätseigenschaft:
die im mitbewegten System ebenfalls stattfindet noch die Rotation des mitbewegten Systems überlagert.
<math>R{{R}^{T}}=1,\det R=1</math>


Also:


Nun , er transformiert sich unter Drehungen wie folgt:


Wenn
:<math>\left( \frac{d}{dt} \right)={{\left( \frac{d}{dt} \right)}^{\acute{\ }}}+\bar{\omega }\times </math>
<math>{{x}_{m}}\to {{x}_{m}}\acute{\ }=\sum\limits_{n=1}^{3}{{{R}_{mn}}{{x}_{n}}}</math>




Dann:
Somit gilt für das körperfeste System
<math>{{J}_{mn}}\to {{J}_{mn}}\acute{\ }=\sum\limits_{l=1}^{3}{\sum\limits_{s=1}^{3}{{}}{{R}_{ml}}{{R}_{ns}}{{J}_{ls}}}</math>
:<math>\bar{K}</math>
:




Kompakt:
:<math>\begin{align}
 
   & \dot{\bar{L}}+\bar{\omega }\times \bar{L}=0 \\
 
  & {{\left( \frac{d}{dt} \right)}^{\acute{\ }}}\bar{L}+\bar{\omega }\times \bar{L}=0 \\
<math>\begin{align}
   & \bar{x}\to \bar{x}\acute{\ }=R\bar{x} \\
  & \bar{\bar{J}}\to \bar{\bar{J}}\acute{\ }=R\bar{\bar{J}}{{R}^{T}} \\
\end{align}</math>
\end{align}</math>




Dabei bemerken wir: Matrizen sind einfach Zahlenschemata mit Zeilen und Spalten. Aber erst das Transformationsverhalten definiert einen Tenor ( Im Gegensatz zu einer Matrix).
'''Mit '''
:<math>\bar{L}=\bar{\bar{J}}\bar{\omega }</math>
'''folgt '''im körperfesten System,wo gilt:
:<math>\dot{\bar{\bar{J}}}</math>
=0


Tensor 1. Stufe:
<math>{{x}_{m}}\acute{\ }=\sum\limits_{n=1}^{3}{{{R}_{mn}}{{x}_{n}}}</math>
= Vektor


Tensor 2. Stufe
:<math>\bar{\bar{J}}\dot{\bar{\omega }}+\bar{\omega }\times \bar{\bar{J}}\bar{\omega }=0</math>
<math>{{J}_{mn}}\acute{\ }=\sum\limits_{l=1}^{3}{\sum\limits_{s=1}^{3}{{}}{{R}_{ml}}{{R}_{ns}}{{J}_{ls}}}</math>




Tensor n-ter STufe:
Dies ist eine nichtlineare DGL in
<math>{{A}_{mn....x}}\acute{\ }=\sum\limits_{l,s,...,t=1}^{3}{{{R}_{ml}}{{R}_{ns}}...{{R}_{xt}}{{A}_{ls...t}}}</math>
:<math>\bar{\omega }</math>
wobei links n Indices stehen und rechts n mal die Drehmatrix angewendet wird ( und jeweils von 1-3 summiert !)
:
 
====Beweis des Transformationsverhaltens für====


<math>{{J}_{mn}}:=\sum\limits_{i=1}^{n}{{}}{{m}_{i}}\left[ \left( \sum\limits_{t}{{{x}_{t}}{{^{(i)}}^{2}}} \right){{\delta }_{mn}}-{{x}_{m}}^{(i)}{{x}_{n}}^{(i)} \right]</math>
Im Schwerpunktsystem ergeben sich die EULERSCHEN Gleichungen für den kräftefreien Kreisel, falls
:<math>\bar{\bar{J}}</math>
diagonal (Hauptträgheitsachsensystem):




Zunächst zum Skalarprodukt:
:<math>\begin{align}
 
   & {{J}_{1}}{{{\dot{\omega }}}_{1}}=\left( {{J}_{2}}-{{J}_{3}} \right){{\omega }_{2}}{{\omega }_{3}} \\
 
& {{J}_{2}}{{{\dot{\omega }}}_{2}}=\left( {{J}_{3}}-{{J}_{1}} \right){{\omega }_{3}}{{\omega }_{1}} \\
<math>\begin{align}
& {{J}_{3}}{{{\dot{\omega }}}_{3}}=\left( {{J}_{1}}-{{J}_{2}} \right){{\omega }_{1}}{{\omega }_{2}} \\
   & {{x}_{m}}\acute{\ }=\sum\limits_{n=1}^{3}{{{R}_{mn}}{{x}_{n}}} \\
& \Rightarrow \sum\limits_{t}{{{x}_{t}}{{\acute{\ }}^{2}}}=\sum\limits_{t}{\sum\limits_{l}{\sum\limits_{s}{{{R}_{tl}}{{R}_{ts}}{{x}_{l}}{{x}_{s}}}=}\sum\limits_{l}{\sum\limits_{s}{\left( \sum\limits_{t}{{}}{{R}_{lt}}^{T}{{R}_{ts}} \right){{x}_{l}}{{x}_{s}}}=\sum\limits_{l}{\sum\limits_{s}{{{\delta }_{ls}}{{x}_{l}}{{x}_{s}}}=\sum\limits_{l}{{{x}_{l}}^{2}}}}} \\
\end{align}</math>
\end{align}</math>




das Skalarprodukt ist also invariant
'''Beispiel:  Symmetrischer Kreisel: '''
 
:<math>{{J}_{1}}={{J}_{2}}\equiv J\ne {{J}_{3}}</math>
Aber auch das Delta- Element ist invariant:
 
 
<math>\sum\limits_{l}{\sum\limits_{s}{{{R}_{ml}}{{R}_{ns}}{{\delta  }_{ls}}}=}\sum\limits_{l}{{{R}_{ml}}{{R}_{nl}}=}\sum\limits_{l}{{{R}_{ml}}{{R}_{\ln }}^{T}={{\delta }_{mn}}}</math>




Kompakt:


 
:<math>{{\dot{\omega }}_{3}}=0</math>,
<math>R1{{R}^{T}}=R{{R}^{T}}=1</math>
also
 
:<math>{{\omega }_{3}}=const</math>
 
im mitrotierenden System
Also:




<math>\begin{align}
:<math>\begin{align}
   & {{J}_{mn}}\acute{\ }=\sum\limits_{l=1}^{3}{\sum\limits_{s=1}^{3}{{}}{{R}_{ml}}{{R}_{ns}}{{J}_{ls}}}=\sum\limits_{i=1}^{n}{{{m}_{i}}}\sum\limits_{l=1}^{3}{\sum\limits_{s=1}^{3}{{}}{{R}_{ml}}{{R}_{ns}}\left[ \left( \sum\limits_{t}{{{x}_{t}}{{^{(i)}}^{2}}} \right){{\delta }_{ls}}-{{x}_{l}}^{(i)}{{x}_{s}}^{(i)} \right]} \\
   & J{{{\dot{\omega }}}_{1}}=\left( {{J}_{{}}}-{{J}_{3}} \right){{\omega }_{2}}{{\omega }_{3}} \\
& {{J}_{mn}}\acute{\ }=\sum\limits_{i=1}^{n}{{{m}_{i}}}\left[ \left( \sum\limits_{t}{{{x}_{t}}^{(i)}{{\acute{\ }}^{2}}} \right){{\delta }_{mn}}-{{x}_{m}}^{(i)}\acute{\ }{{x}_{n}}^{(i)}\acute{\ } \right] \\
& {{{\ddot{\omega }}}_{1}}=\frac{\left( {{J}_{{}}}-{{J}_{3}} \right)}{J}{{{\dot{\omega }}}_{2}}{{\omega }_{3}}=-{{\left[ \frac{\left( {{J}_{{}}}-{{J}_{3}} \right)}{J}{{\omega }_{3}} \right]}^{2}}{{\omega }_{1}} \\
\end{align}</math>
\end{align}</math>




Der Trägheitstensor J´ in den neuen Koordinaten ist also gleich dem alten, was Transformationsverhalten eines Tensors zweiter Stufe belegt:
Diese Gleichung kann zweimal integriert werden. Mit den Integrationskonstanten


Dabei gilt:


:<math>{{\omega }_{\bot }},{{\phi }_{0}}</math>
und der Zusammenfassung
:<math>{{\omega }_{0}}:=\frac{\left( J-{{J}_{3}} \right)}{J}{{\omega }_{3}}</math>


<math>\left( \sum\limits_{t}{{{x}_{t}}^{(i)}{{\acute{\ }}^{2}}} \right){{\delta }_{mn}}</math>
ist der invariante Anteil


 
folgt:
<math>{{x}_{m}}^{(i)}\acute{\ }{{x}_{n}}^{(i)}\acute{\ }</math>
hängt von der Wahl des körperfesten koordinatensystems ab.
 
<u>'''Weitere Eigenschaften'''</u>
 
#
<math>{{J}_{mn}}</math>
enthält einen kugelsymmetrischen, also rotationsinvarianten Anteil
<math>\left( \sum\limits_{t}{{{x}_{t}}{{^{(i)}}^{2}}} \right){{\delta }_{mn}}</math>
 
#
<math>{{J}_{mn}}</math>
ist linear in der Massendichte. Der Trägheitstensor ist also additiv beim Zusammenfügen zweier starrer Körper
#
<math>{{J}_{mn}}</math>
ist ein reeller, symmetrischer Tensor, dargestellt durch die reelle, symmetrisch Matrix
 
 
<math>\bar{\bar{J}}=\int_{{}}^{{}}{{{d}^{3}}x\rho (\bar{x})\left( \begin{matrix}
  {{x}_{2}}^{2}+{{x}_{3}}^{2} & -{{x}_{1}}{{x}_{2}} & -{{x}_{1}}{{x}_{3}}  \\
  -{{x}_{1}}{{x}_{2}} & {{x}_{3}}^{2}+{{x}_{1}}^{2} & -{{x}_{2}}{{x}_{3}}  \\
  -{{x}_{1}}{{x}_{3}} & -{{x}_{2}}{{x}_{3}} & {{x}_{1}}^{2}+{{x}_{2}}^{2}  \\
\end{matrix} \right)}</math>
 
 
Der Tensor ist diagonalisierbar durch die orthogonale Transformation
<math>{{R}_{0}}\in SO(3):</math>
 
 
 
<math>\bar{\bar{J}}\acute{\ }={{R}_{0}}\bar{\bar{J}}{{R}_{0}}^{T}=\left( \begin{matrix}
  {{J}_{1}} & 0 & 0  \\
  0 & {{J}_{2}} & 0  \\
  0 & 0 & {{J}_{3}}  \\
\end{matrix} \right)</math>
 
 
Das heißt: Es existiert ein gedrehtes, körperfestes Koordinatensystem (y1,y2,y3)  in Richtung der '''Hauptträgheitsachsen:'''
 
 
<math>\bar{\bar{J}}\acute{\ }=\int_{{}}^{{}}{{{d}^{3}}y\rho (\bar{y})\left( \begin{matrix}
  {{y}_{2}}^{2}+{{y}_{3}}^{2} & 0 & 0  \\
  0 & {{y}_{3}}^{2}+{{y}_{1}}^{2} & 0  \\
  0 & 0 & {{y}_{1}}^{2}+{{y}_{2}}^{2}  \\
\end{matrix} \right)}</math>
 
 
Also:
<math>{{J}_{i}}\ge 0</math>
i=1,..,3, Matrix positiv semidefinit.
 
Die Diagonalisierung führt auf das Eigenwertproblem:
 
 
<math>\bar{\bar{J}}{{\hat{\bar{w}}}^{(i)}}={{J}_{i}}{{\hat{\bar{w}}}^{(i)}}</math>
mit Eigenvektoren
<math>{{\hat{\bar{w}}}^{(i)}}</math>
und Eigenwerten Ji. Ein homogenes, lineares Gleichungssystem
 
Ziel ist es nun, die Hauptachsenrichtung
<math>{{\hat{\bar{w}}}^{(i)}}</math>
so zu suchen, dass
<math>\bar{\bar{J}}</math>
diagonal wird:
 
 
<math>\Leftrightarrow \det \left( \bar{\bar{J}}-{{J}_{i}}1 \right)=0</math>




Somit ergeben sich 3 reelle, positiv semidefinite Eigenwerte Ji
:<math>{{\omega }_{1}}={{\omega }_{\bot }}\cos \left( {{\omega }_{0}}t+{{\phi }_{0}} \right)</math>


====Das Trägheitsmoment====


<u>'''Trägheitsmoment bezüglich Achse '''</u>
Dies kann in
<math>\bar{n}:J(\bar{n}):=\bar{n}\bar{\bar{J}}\bar{n}</math>
Diese quadratische Form ist positiv semidefinit.




<math>\bar{\omega }=\bar{n}\omega \Rightarrow T=\frac{1}{2}{{\omega }^{2}}J(\bar{n})</math>
:<math>J{{\dot{\omega }}_{1}}=\left( {{J}_{{}}}-{{J}_{3}} \right){{\omega }_{2}}{{\omega }_{3}}</math>
eingesetzt werden und es ergibt sich:




<u>'''Trägheitsellipsoid'''</u>
:<math>{{\omega }_{2}}=-{{\omega }_{\bot }}\sin \left( {{\omega }_{0}}t+{{\phi }_{0}} \right)</math>


Die Normierung des Trägheitsmomentes liefert eine Ellipsoidgleichung:
<math>\bar{n}\bar{\bar{J}}\bar{n}=1</math>
.


Die Lage des Ellipsoids sind ist durch die Eigenvektoren
Die Definitionen sind an folgender Figur ersichtlich:
<math>{{\hat{\bar{w}}}^{(i)}}</math>
, die Maße folgen aus den Ji derart, dass die zu
<math>{{\hat{\bar{w}}}^{(i)}}</math>
gehörige Achse die Länge
<math>\frac{1}{\sqrt{{{J}_{i}}}}</math>
trägt:


# Die Ji heißen Hauptträgheitsmomente ( Trägheitsmomente entlang der Eigenvektoren= Hauptachsen)
Dabei ist x3 die Figurenachse (Achse durch die Drehachse von J3)


Es gilt:
Es gilt:




<math>{{J}_{1}}\ne {{J}_{2}}\ne {{J}_{3}}</math>
:<math>{{\omega }_{1}}{{(t)}^{2}}+{{\omega }_{2}}{{(t)}^{2}}={{\omega }_{\bot }}^{2}=const</math>
unsymmetrischer Kreisel
 
 
<math>{{J}_{1}}={{J}_{2}}\ne {{J}_{3}}</math>
symmetrischer Kreisel ( axialsymmetrisch)
 
 
<math>{{J}_{1}}={{J}_{2}}={{J}_{3}}</math>
kugelsymmetrischer Kreisel ( nicht notwendigerweise Kugelform)
 
====Satz von Steiner====
 
Sei''' '''
<math>{{J}_{mn}}</math>
der Trägheitstensor in einem körperfesten System
<math>\bar{K}</math>
, welches im Schwerpunkt S zentriert ist. Sei nun
<math>\bar{K}\acute{\ }</math>
ein zu
<math>\bar{K}</math>
achsparalleles, um den Vektor
<math>\bar{a}</math>
verschobenes System. Dann ist
<math>{{J}_{mn}}\acute{\ }</math>
in
<math>\bar{K}\acute{\ }</math>
gegeben durch
 
 
<math>{{J}_{mn}}\acute{\ }={{J}_{mn}}+M\left[ {{a}^{2}}{{\delta }_{mn}}-{{a}_{m}}{{a}_{n}} \right]</math>
 
 
Die beiden Koordinatensystem dürfen dabei nur durch die Translation um
<math>\bar{a}</math>
unterschiedlich sein. Wesentlich ist vor allem, dass bei roationsvarianten Systemen keine Verdrehung der Achsen erfolgt !
 
Beweis:
 
 
<math>{{J}_{mn}}\acute{\ }=\int_{{}}^{{}}{{{d}^{3}}x\acute{\ }}\rho \acute{\ }(\bar{x}\acute{\ })\left[ \left( \sum\limits_{t}{{{x}_{t}}{{\acute{\ }}^{2}}} \right){{\delta }_{mn}}-{{x}_{m}}\acute{\ }{{x}_{n}}\acute{\ } \right]</math>
 
 
Bei uns:
<math>\bar{x}\acute{\ }=\bar{x}+\bar{a}</math>
 
 
 
<math>\begin{align}
  & {{J}_{mn}}\acute{\ }=\int_{{}}^{{}}{{{d}^{3}}x}\rho (\bar{x})\left[ \left( \sum\limits_{t}{{{\left( {{x}_{t}}+{{a}_{t}} \right)}^{2}}} \right){{\delta }_{mn}}-\left( {{x}_{m}}+{{a}_{m}} \right)\left( {{x}_{n}}+{{a}_{n}} \right) \right] \\
& {{J}_{mn}}\acute{\ }=\int_{{}}^{{}}{{{d}^{3}}x}\rho (\bar{x})\left[ \left( \sum\limits_{t}{\left[ {{\left( {{x}_{t}} \right)}^{2}}+2\left( {{a}_{t}}{{x}_{t}} \right)+{{a}_{t}}^{2} \right]} \right){{\delta }_{mn}}-{{x}_{m}}{{x}_{n}}-{{x}_{m}}{{a}_{n}}-{{x}_{n}}{{a}_{m}}-{{a}_{m}}{{a}_{n}} \right] \\
& \int_{{}}^{{}}{{{d}^{3}}x}\rho (\bar{x})\sum\limits_{t}{\left( {{a}_{t}}{{x}_{t}} \right)}=\int_{{}}^{{}}{{{d}^{3}}x}\rho (\bar{x})\left( {{x}_{m}}{{a}_{n}}+{{x}_{n}}{{a}_{m}} \right)=0\quad wegen\ \int_{{}}^{{}}{{{d}^{3}}x}\rho (\bar{x})\bar{x}=0 \\
\end{align}</math>
 


Somit:




<math>\begin{align}
:<math>{{\omega }_{1}}{{(t)}^{2}}+{{\omega }_{2}}{{(t)}^{2}}+{{\omega }_{3}}{{(t)}^{2}}={{\omega }_{\bot }}^{2}+{{\omega }_{3}}{{(t)}^{2}}=const</math>
  & {{J}_{mn}}\acute{\ }=\int_{{}}^{{}}{{{d}^{3}}x}\rho (\bar{x})\left[ \left( \sum\limits_{t}{\left( {{x}_{t}}^{2}+{{a}_{t}}^{2} \right)} \right){{\delta }_{mn}}-{{x}_{m}}{{x}_{n}}-{{a}_{m}}{{a}_{n}} \right] \\
& {{J}_{mn}}\acute{\ }={{J}_{mn}}+\int_{{}}^{{}}{{{d}^{3}}x}\rho (\bar{x})\left[ \left( \sum\limits_{t}{\left( {{a}_{t}}^{2} \right)} \right){{\delta }_{mn}}-{{a}_{m}}{{a}_{n}} \right]={{J}_{mn}}+M\left[ {{a}^{2}}{{\delta }_{mn}}-{{a}_{m}}{{a}_{n}} \right] \\
\end{align}</math>




'''Speziell im Hauptachsensystem:'''
Das heißt


keine Außerdiagonalelemente: m=n:=i


:<math>\bar{\omega }</math>
und  damit auch
:<math>\bar{L}</math> mit <math>{{L}_{i}}={{J}_{i}}{{\omega }_{i}}</math>
rotieren um die Figurenachse
:<math>\bar{f}||{{x}_{3}}</math>


<math>{{J}_{i}}\acute{\ }={{J}_{i}}+M({{a}^{2}}-{{a}_{i}}^{2})\quad i=1,..,3</math>


Veranschaulicht man diese Situation im Schwerpunktsystem mit '''raumfesten '''Achsen, so gilt mit
:<math>\frac{d}{dt}\bar{L}=0\Rightarrow \bar{L}\ fest</math>


mit
<math>({{a}^{2}}-{{a}_{i}}^{2})</math>
als Quadrat des Abstandes der beiden Drehachsen.


Dabei wird bei einer Verschiebung um
<math>\bar{a}</math>
nur der Abstand der Drehachsen berücksichtigt. das heißt, die Komponente der Verschiebung in Richtung der Drehachse wird wieder quadratisch subtrahiert:


:<math>\bar{\omega }</math> und <math>\bar{f}</math>
präzedieren um die raumfeste Achse
:<math>\bar{L}</math>.
Dabei müssen
:<math>\bar{\omega }</math>,


:<math>\bar{f}</math> und <math>\bar{L}</math>
stets in einer Ebene liegen.


<u>'''Beispiele'''</u>
'''Anwendung:'''


1. Kugelsymmetrische Massendichte:
Erde als abgeplattetes Rotationsellipsoid:




<math>\begin{align}
:<math>\begin{align}
   & \rho (\bar{x})=\rho (r) \\
   & \frac{\left( {{J}_{{}}}-{{J}_{3}} \right)}{J}\approx \frac{1}{300} \\
& {{J}_{1}}={{J}_{2}}={{J}_{3}}=:J \\
  & \frac{2\pi }{{{\omega }_{3}}}=1Tag \\
& 3J={{J}_{1}}+{{J}_{2}}+{{J}_{3}}=\int_{{}}^{{}}{{{d}^{3}}x}\rho (r)\left[ \left( {{x}_{2}}^{2}+{{x}_{3}}^{2} \right)+\left( {{x}_{1}}^{2}+{{x}_{3}}^{2} \right)+\left( {{x}_{1}}^{2}+{{x}_{2}}^{2} \right) \right]=\int_{{}}^{{}}{{{d}^{3}}x}\rho (r)2{{r}^{2}} \\
  & 3J=2\cdot 4\pi \int_{0}^{R}{dr{{r}^{4}}\rho (r)} \\
\end{align}</math>
\end{align}</math>




Bei homogener Massenverteilung:
Damit kann die Präzessionsperiode leicht berechnet werden:
 
 
<math>M=\frac{4\pi }{3}{{R}^{3}}</math>
bezüglich Schwerpunkt S
 
folgt:
 
 
<math>\begin{align}
  & J=\frac{8}{3}\pi \int_{0}^{R}{dr{{r}^{4}}\rho (r)}=\frac{2M}{{{R}^{3}}}\int_{0}^{R}{dr{{r}^{4}}} \\
& J=\frac{2}{5}M{{R}^{2}} \\
\end{align}</math>
 


2. Abrollende Kugel: Momentaner Auflagepunkt ist A


Das Trägheitsmoment bezüglich der momentanen Drehachse durch den Auflagepunkt A:
:<math>T=\frac{2\pi }{{{\omega }_{0}}}=\frac{2\pi J}{{{\omega }_{3}}(J-{{J}_{3}})}=\frac{J}{(J-{{J}_{3}})}\cdot 1Tag=300Tage</math>




<math>{{J}_{A}}=J+M{{R}^{2}}=\frac{2}{5}M{{R}^{2}}+M{{R}^{2}}=\frac{7}{5}M{{R}^{2}}</math>
Die Erde präzediert also einmal in 300 Tagen um ihre eigene Achse!
rac{2}{5}M{{R}^{2}}+M{{R}^{2}}=\frac{7}{5}M{{R}^{2}}</math>

Aktuelle Version vom 12. September 2010, 23:27 Uhr




Drehimpuls

  1. diskret:
l¯=i=1nmir¯i×r¯˙i=i=1nmi(r¯S+x¯(i))×(V¯+ω¯×x¯(i))l¯=i=1nmi(r¯S×V¯)+i=1nmix¯(i)×V¯+r¯S×(ω¯×imix¯(i))+imix¯(i)×(ω¯×x¯(i))i=1nmix¯(i)×V¯=i=1nmi(x¯(i))=0l¯=i=1nmi(r¯S×V¯)+imix¯(i)×(ω¯×x¯(i))=M(r¯S×V¯)+imix¯(i)×(ω¯×x¯(i)) Mit M(r¯S×V¯)
als Schwerpunktsdrehimpuls


imix¯(i)×(ω¯×x¯(i))
als Relativdrehimpuls
  1. kontinuierliche Situation


l¯=r¯s×MV¯+d3xρ(x¯)x¯×(ω¯×x¯)L¯=d3xρ(x¯)x¯×(ω¯×x¯)=d3xρ(x¯)[x2ω¯(x¯ω¯)x¯]


Also:


L¯=J¯¯ω¯


Dies sieht man an der Komponentenschreibweise:


Lm=n=13d3xρ(x¯)[x2δmnxmxn]ωn=n=13Jmnωn


Nebenbemerkung:

Im Allgemeinen ist

L¯

nicht parallel zu

ω¯,
nur falls
ω¯

in Richtung der Hauptträgheitsachse liegt!

Allgemeine Bewegungsgleichung für den Gesamtdrehimpuls

ddtl¯=ir¯i×F¯i.
Dabei sind
F¯i

äußere, eingeprägte Kräfte. Die resultierende Kraft

F¯(r¯S)=iF¯i

soll auf den Schwerpunkt wirken, so dass gilt:


iF¯imi=F¯(r¯S)M


Somit:


ddtl¯=imir¯i×F¯imi=r¯S×F¯(r¯S)


Bekanntlich gilt für die Schwerpunktsbewegung:


Mr¯¨S=F¯(r¯S)

(Newton)


ddtl¯=ddt(r¯s×MV¯+d3xρ(x¯)x¯×(ω¯×x¯))=ddt(r¯s×Mr¯˙s)+ddtL¯=Mr¯˙s×r¯˙s+r¯S×Mr¯¨S+ddtL¯


Mr¯˙s×r¯˙s=0r¯S×Mr¯¨S=r¯S×F¯(r¯S)


ddtl¯=r¯S×F¯(r¯S)+ddtL¯


Gleichzeitig gilt:

ddtl¯=imir¯i×F¯imi=r¯S×F¯(r¯S)


Somit:


ddtL¯=0
Der Relativdrehimpuls ist im Schwerpunktsystem mit raumfesten Achsen konstant.

Also: Es verschwindet die Zeitableitung des relativdrehimpulses im Schwerpunktsystem mit RAUMFESTEN Achsen.

Die Transformation muss nun noch ins körperfeste, rotatorisch mitbewegte System

K¯

erfolgen:

Dabei sieht der Beobachter im RAUMFESTEN System neben der zeitlichen Änderung

(ddt)´,
die im mitbewegten System ebenfalls stattfindet noch die Rotation des mitbewegten Systems überlagert.

Also:


(ddt)=(ddt)´+ω¯×


Somit gilt für das körperfeste System

K¯


L¯˙+ω¯×L¯=0(ddt)´L¯+ω¯×L¯=0


Mit

L¯=J¯¯ω¯

folgt im körperfesten System,wo gilt:

J¯¯˙

=0


J¯¯ω¯˙+ω¯×J¯¯ω¯=0


Dies ist eine nichtlineare DGL in

ω¯

Im Schwerpunktsystem ergeben sich die EULERSCHEN Gleichungen für den kräftefreien Kreisel, falls

J¯¯

diagonal (Hauptträgheitsachsensystem):


J1ω˙1=(J2J3)ω2ω3J2ω˙2=(J3J1)ω3ω1J3ω˙3=(J1J2)ω1ω2


Beispiel: Symmetrischer Kreisel:

J1=J2JJ3


ω˙3=0,
also
ω3=const

im mitrotierenden System


Jω˙1=(JJ3)ω2ω3ω¨1=(JJ3)Jω˙2ω3=[(JJ3)Jω3]2ω1


Diese Gleichung kann zweimal integriert werden. Mit den Integrationskonstanten


ω,ϕ0

und der Zusammenfassung

ω0:=(JJ3)Jω3


folgt:


ω1=ωcos(ω0t+ϕ0)


Dies kann in


Jω˙1=(JJ3)ω2ω3

eingesetzt werden und es ergibt sich:


ω2=ωsin(ω0t+ϕ0)


Die Definitionen sind an folgender Figur ersichtlich:

Dabei ist x3 die Figurenachse (Achse durch die Drehachse von J3)

Es gilt:


ω1(t)2+ω2(t)2=ω2=const


ω1(t)2+ω2(t)2+ω3(t)2=ω2+ω3(t)2=const


Das heißt


ω¯

und damit auch

L¯ mit Li=Jiωi
rotieren um die Figurenachse
f¯||x3


Veranschaulicht man diese Situation im Schwerpunktsystem mit raumfesten Achsen, so gilt mit

ddtL¯=0L¯fest


ω¯ und f¯

präzedieren um die raumfeste Achse

L¯.
Dabei müssen
ω¯,
f¯ und L¯

stets in einer Ebene liegen.

Anwendung:

Erde als abgeplattetes Rotationsellipsoid:


(JJ3)J13002πω3=1Tag


Damit kann die Präzessionsperiode leicht berechnet werden:


T=2πω0=2πJω3(JJ3)=J(JJ3)1Tag=300Tage


Die Erde präzediert also einmal in 300 Tagen um ihre eigene Achse! rac{2}{5}M{{R}^{2}}+M{{R}^{2}}=\frac{7}{5}M{{R}^{2}}</math>