Trägheitstensor - LinkFang.de





Trägheitstensor


Der Trägheitstensor [math]I[/math] eines Körpers gibt seine Trägheitsmomente an, also die Trägheit des Körpers bezüglich der Drehungen. Er spielt damit für Drehungen dieselbe Rolle wie die träge Masse für lineare Bewegungen.

Da nicht-kugelförmige Körper für Drehungen um verschiedene Achsen im Allgemeinen verschiedene Trägheitsmomente aufweisen (beispielsweise lässt sich die Drehung eines homogenen zigarrenförmigen Körpers leichter um seine Längsachse als um seine Querachse ändern), reicht – anders als bei der trägen Masse – für die Beschreibung der Trägheitsmomente eine einzelne Zahl nicht aus, sondern es muss ein Tensor verwendet werden.

Berechnung

Ist eine Gesamtmasse gegeben durch einzelne Massenpunkte, so ist der Trägheitstensor [math]I[/math] gegeben durch:

[math]I = \sum_{i} m_{i} \begin{pmatrix} y_i^2+z_i^2 & -x_i y_i & -x_i z_i \\ -y_i x_i & x_i^2+z_i^2 & -y_i z_i \\ -z_i x_i & - z_i y_i & x_i^2+y_i^2 \\ \end{pmatrix} [/math]

bzw. in Komponentenschreibweise

[math]I_{\alpha\beta}=I_{\beta\alpha}= \sum_{i} m_{i} (r_{i}^2 \delta _{\alpha\beta} - x_{i\alpha} x_{i\beta})[/math]

Hierbei bezeichnen

  • [math]m_i[/math] die Masse des [math]i[/math]-ten Massenpunkts
  • [math]x_{i1} = x_i[/math], [math]x_{i2} = y_i[/math] und [math]x_{i3} = z_i[/math] die Koordinaten seines Ortsvektors
  • [math]r_i = \sqrt{x_i^2+y_i^2+z_i^2}[/math] seinen Abstand vom Ursprung.

Die Indizes [math]\alpha[/math] und [math]\beta[/math] nehmen die Werte 1 bis 3 an, entsprechend den drei Raumdimensionen.

[math]\delta_{\alpha\beta}[/math] ist das Kronecker-Delta:

[math]\delta_{\alpha\beta}=\left\{\begin{matrix}1&\mbox{ für }\alpha=\beta\\0&\mbox{ sonst}\end{matrix}\right.[/math]

Ist die Masse kontinuierlich im Körper verteilt und sind die Massendichten [math]\rho(\vec{r})[/math] bekannt, so geht man zur Integration über:

[math]I_{\alpha\beta}=\int \limits_{V} \rho(\vec{r})(r^2 \delta_{\alpha\beta} - x_{\alpha} x_{\beta}) \,\mathrm{d}V[/math]

Man nennt die Diagonalelemente des Trägheitstensors Trägheitsmomente und die restlichen Elemente Deviationsmomente. Diese bewirken bei Rotation durch Fliehkräfte ein Drehmoment normal zur Drehachse. Man stelle sich z. B. eine Stange vor, die um eine Achse rotiert, die mit der Symmetrieachse der Stange einen spitzen Winkel einschließt.

Beispiel: Der homogene Würfel mit Kantenlänge 2a

Trägheitstensoren werden im Allgemeinen in kartesischen Koordinaten berechnet. Aus der Summe in der oben angegebenen Definition

[math]I_{\alpha\beta}=I_{\beta\alpha}= \sum_{i} m_{i} (r_{i}^2 \delta _{\alpha\beta} - x_{i\alpha} x_{i\beta})[/math]

lässt sich bei konstanter Massenverteilung mit überall bekannter Dichte [math] \varrho(r) [/math], also [math] dm = dV \varrho (r) [/math] ein Integral machen. Das ist ein gängiger Trick in der Kontinuumsmechanik. Statt über die zahlreichen Einzelmassen (Atome, Moleküle) zu summieren und ihren jeweiligen Abstand [math] r_i [/math] zum Bezugspunkt zu berücksichtigen, geht man zu einem Integral der Form

[math]I_{\alpha\beta}=I_{\beta\alpha}= \int \limits_{V} (r^2 \delta _{\alpha\beta} - x_{\alpha} x_{\beta}) \varrho(r) dV [/math]

über. In dem Spezialfall, dass der Körper homogen ist, also über eine durchgehend konstante Dichte verfügt, kann man die Dichte als Konstante vor das Integral ziehen und erhält [math]I_{\alpha\beta}=I_{\beta\alpha}= \varrho \int \limits_{V} (r^2 \delta _{\alpha\beta} - x_{\alpha} x_{\beta}) dV.[/math] Nun lassen sich die sechs (drei Haupt- und nur drei Deviationsmomente, da die Matrix symmetrisch ist, das heißt, dass [math] I_{xy} = I_{yx} [/math] gilt) unabhängigen Tensorkomponenten bestimmen. Hier wie auch im allgemeinen Fall ist zu beachten, dass der einzige intellektuelle Anspruch in der Bestimmung der Integrationsgrenzen liegt. Hier ist klar, dass in allen drei Richtungen von -a bis a integriert werden muss. Im Allgemeinen muss man einen Körper in Teilkörper mit geometrisch einfach beschreibbaren Formen zerlegen, dann ergibt sich die Beschreibung der Form durch Abhängigkeit der Integrationsgrenzen von den zu integrierenden Parametern (z. B.

[math] A = \int_{0}^{b} \int_{0}^{(b-y)} dx dy = b^2/2[/math]

für ein in der Diagonale zerschnittenes Quadrat). Für den Würfel ergibt sich:

[math] I_{xx} = \varrho \int \limits_{V} (y^2 + z^2) dV = \varrho \frac{16}{3} a^5 [/math]
[math] I_{yy} = \varrho \int \limits_{V} (x^2 + z^2) dV = \varrho \frac{16}{3} a^5[/math]
[math] I_{zz} = \varrho \int \limits_{V} (y^2 + x^2) dV = \varrho \frac{16}{3} a^5[/math]
[math] I_{xy} = I_{yx} = -\varrho \int \limits_{V} yx dV = 0 [/math]
[math] I_{yz} = I_{zy} = -\varrho \int \limits_{V} zy dV = 0 [/math]
[math] I_{zx} = I_{xz} = -\varrho \int \limits_{V} xz dV = 0 [/math]

Dabei wird benutzt, dass

[math] \int \limits_{V} x^2 dV = \int \limits_{-a}^{+a} \int \limits_{-a}^{+a} \int \limits_{-a}^{+a} x^2 dx dy dz = (\int \limits_{-a}^{+a} x^2 dx ) (2a)^2 = ( \frac{2}{3}a^3 ) 4a^2 = \frac{8}{3}a^5 [/math]

und analog für [math]\int \limits_{V} y^2 dV, \int \limits_{V} z^2 dV [/math].

Weiterhin ist zu berücksichtigen, dass die Masse des Würfels

[math] M = V \varrho = (2a)^3 \varrho = 8a^3 \varrho [/math]

ist. Dann hat der Tensor die Form

[math] I = \frac{M}{8} \begin{pmatrix} \frac{16}{3}a^2 & 0 & 0 \\ 0 & \frac{16}{3}a^2 & 0 \\ 0 & 0 & \frac{16}{3}a^2 \\ \end{pmatrix} = \frac{2}{3} M a^2 \begin{pmatrix} 1 & 0 & 0\\ 0 & 1 & 0 \\ 0 & 0 & 1 \\ \end{pmatrix} = \frac{2}{3} M a^2 E, [/math]

wobei E die Einheitsmatrix ist.

Setzt man nun noch die Kantenlänge 2a = d, so ergibt sich

[math]I = \frac{1}{6} M d^2 E[/math]

Berechnungen aus dem Trägheitstensor

Trägheitsmomente

Das Trägheitsmoment bezüglich einer Achse durch den Ursprung (also anschaulich die Schwierigkeit, den Körper um diese Achse zu drehen) berechnet man mit Hilfe des Trägheitstensors folgendermaßen:

Ist [math]\vec n[/math] ein Einheitsvektor (Vektor der Länge 1) in Richtung der Achse, so ist das zugehörige Trägheitsmoment

[math]I_n = \sum_{\alpha=1}^3 \sum_{\beta=1}^3 I_{\alpha\beta} n_\alpha n_\beta=\vec n ^T \cdot \left( I\cdot \vec n \right),[/math]

wobei [math]n_\alpha[/math] die Komponenten des Einheitsvektors sind. Insbesondere geben die Diagonalelemente [math]I_{11}[/math], [math]I_{22}[/math] und [math]I_{33}[/math] die Trägheitsmomente um die Koordinatenachsen an.

Beispiel: Trägheitsmoment entlang der Diagonalen

Der Trägheitstensor sei

[math]I=\mathrm{\left(\begin{matrix} 4 & 0 & 0\\ 0 & 2 & 0\\ 0 & 0 & 6 \end{matrix}\right)}kg \cdot m^2,[/math]

und die Richtung der Achse sei gegeben durch den Normalenvektor

[math]\vec n = \frac{1}{\sqrt{3}}\left(\begin{matrix}1\\1\\1\end{matrix}\right)[/math]

(das ist gerade eine der Raumdiagonalen). Dann ist das Trägheitsmoment gegeben durch

[math]\mathrm{I_n = \frac{1}{3} (1 \cdot 4 kg \cdot m^2 \cdot 1 + 1 \cdot 2 kg \cdot m^2 \cdot 1 + 1 \cdot 6 kg \cdot m^2 \cdot 1) = 4 kg \cdot m^2}.[/math]

Drehimpuls

Für eine beliebige Winkelgeschwindigkeit [math]\vec \omega[/math] lässt sich der Drehimpuls durch Matrizenmultiplikation des Trägheitstensors mit dem Spaltenvektor der Winkelgeschwindigkeit berechnen: [math]\vec L = I \ \vec \omega[/math]

Eigenschaften

Der Trägheitstensor eines ausgedehnten Körpers hat den Rang 3 und die Dimension eines Trägheitsmoments, also 1 kg·m² als SI-Einheit.

Er ist ein symmetrischer Tensor zweiter Stufe, also

[math]I_{ij} = I_{ji}.[/math]

Hauptträgheitsachsen und Hauptträgheitsmomente

Wählt man als Koordinatenachsen die Achse des größten und die des kleinsten Trägheitsmoments sowie die auf diesen beiden senkrecht stehende Achse, so ist der Trägheitstensor in diesen Koordinaten diagonal. (Sind alle Trägheitsmomente gleich, so können drei beliebige aufeinander senkrecht stehende Achsen gewählt werden.) Diese drei Achsen nennt man Hauptträgheitsachsen, die zugehörigen reellen und nichtnegativen Eigenwerte [math]I_1, I_2, I_3[/math] des Trägheitstensors Hauptträgheitsmomente. Per Definition gilt:

[math]I_{max} = I_{1} \geq I_{2} \geq I_{3} = I_{min}.[/math]

Außerdem:

[math]I_{1} + I_{2} \geq I_{3},[/math]

auch bei beliebiger Vertauschung untereinander (d. h. [math]I_{1} + I_{3} \geq I_{2}[/math] und [math]I_{2} + I_{3} \geq I_{1}[/math]).

  • Sind alle drei Hauptträgheitsmomente untereinander gleich: [math]I_1 = I_2 = I_3[/math], so ist der Körper punktsymmetrisch wie eine Kugel oder ein Würfel.
  • Sind nur zwei der drei Hauptträgheitsmomente untereinander gleich, z. B. [math]I_1 = I_2 \not = I_3[/math], so handelt es sich um einen rotationssymmetrischen Körper, wie etwa einen runden Zylinder, oder z. B. um einen Stab mit quadratischer oder hexagonaler Grundfläche.
  • Sind [math]I_1[/math], [math]I_2[/math] und [math]I_3[/math] paarweise voneinander verschieden, so liegt keine Symmetrie vor, und/oder der Koordinatenursprung liegt nicht im Massenmittelpunkt.

Die Drehimpulse bei Rotation um die Hauptträgheitsachsen können durch einfache Multiplikation des jeweiligen (skalaren) Hauptträgheitsmoments mit dem Winkelgeschwindigkeitsvektor berechnet werden, z. B.:

[math]\vec {L_1} = I_1 \ \vec \omega[/math]

Drehung nicht um die Hauptträgheitsachsen

Häufig liegt der Trägheitstensor nicht diagonal vor. Da er symmetrisch ist, lässt er sich jedoch immer diagonaliseren, d. h. für jeden Punkt im starren Körper lassen sich die Hauptträgheitsachsen finden. Im Grunde ist dies ein Eigenwertproblem: zunächst werden die Hauptträgheitsmomente als Eigenwerte des Trägheitstensors berechnet. Anschließend berechnet man die drei zugehörigen Richtungen bzw. Hauptträgheitsachsen als Eigenvektoren aus den Eigenwerten und [math]I[/math]. Die Eigenvektoren bilden eine orthogonale Basis im Raum, d. h. sie stehen senkrecht aufeinander. Sie lassen sich zu einer Transformationsmatrix [math]T[/math] zusammenfassen, welche die Informationen über unser neu ausgerichtetes Koordinatensystem trägt.

Für jeden nicht-diagonalen Tensor [math]I[/math] gibt es also eine Transformationsmatrix [math]T \in \mathbb{R}^{3\times 3}[/math], sodass gilt:

[math]I' = T^{-1} \cdot I \cdot T[/math]

wobei [math]I'[/math] nun diagonalisiert ist, seine Einträge den Hauptträgheitsmomenten entsprechen, und [math]T^{-1}[/math] die zu [math]T[/math] inverse Matrix ist.

Die freie Rotation (ohne dass die Drehachse durch Lager festgehalten wird) ist für Körper im Allgemeinen instabil, d. h. es treten Winkelbeschleunigungen auf. Nur die freie Rotation um die Hauptträgheitsachse mit dem größten oder dem kleinsten Trägheitsmoment ist stabil; die Rotation um die Hauptträgheitsachse mit dem mittleren Trägheitsmoment ist dagegen labil.

Dreht sich der Körper nicht um eine Hauptträgheitsachse, so zeigt der Drehimpuls nicht in Richtung der Drehachse. Daher bleibt die momentane Drehachse nicht fest im Raum, sondern läuft auf der Oberfläche des Rastpolkegels um die Drehimpulsachse herum. Der Körper führt dann eine Taumelbewegung (Präzession) aus.

Wirkt außerdem ein äußeres Drehmoment, so wird die Bewegung noch komplizierter. Man kann dies z. B. an einem Spielzeugkreisel beobachten: zusätzlich zu einer – hier durch das äußere Drehmoment induzierten – Präzession um eine senkrechte Achse (die Rotationsachse wandert am Präzessionskegel entlang) kann man durch einen kleinen Stoß eine überlagerte Nutationsbewegung erzeugen.

Herleitung

Für den Drehimpuls [math]\vec{L}[/math] eines Teilchens gilt

[math]\begin{align} \vec{L} & = \vec{r} \times \vec{p}\\ & = \vec{r} \times (m \cdot \vec{v})\\ & = m \cdot \vec{r} \times (\vec{\omega} \times \vec{r}) \end{align}[/math]

mit

  • [math]\vec{r}[/math] dem Abstand vom Ursprung
  • [math]\vec{p}[/math] dem Impuls
  • [math]m[/math] der Masse
  • [math]\vec{v}[/math] der Geschwindigkeit.

Dann gilt für die [math]i[/math]-te Komponente des Drehimpulses:

[math]L_i = (m\cdot\vec{r}\times(\vec{\omega}\times\vec{r}))_i[/math]
[math]=m~\varepsilon_{ijk}~r_j~\varepsilon_{klm}~\omega_l~r_m[/math]
[math]=m~\varepsilon_{kij}~\varepsilon_{klm}~r_j~\omega_l~r_m[/math]
[math]=m~(\delta_{il}\delta_{jm} - \delta_{im}\delta_{jl})~r_j~\omega_l~r_m[/math]
[math]=m~(\delta_{il}\delta_{jm}~r_j~\omega_l~r_m - \delta_{im}\delta_{jl}~r_j~\omega_l~r_m)[/math]
[math]=m~(r_j~\omega_i~r_j - r_j~\omega_j~r_i)[/math]
[math]=m~(\vec{r}\,^2~\omega_i - r_i~r_j~\omega_j)[/math]
[math]=\underbrace{m~(\vec{r}\,^2~\delta_{ij} - r_i~r_j)}_{=I_{ij}} \omega_j[/math]

Dabei ist [math]\varepsilon[/math] das Levi-Civita-Symbol und [math]\delta[/math] das Kronecker-Delta. Außerdem wurde die Einsteinsche Summenkonvention benutzt.

Damit ist nun

[math]I_{ij} = m\cdot(\vec{r}\,^2~\delta_{ij} - r_i~r_j).[/math]

Bei einem Körper aus mehreren Massepunkten muss über diese summiert oder integriert werden.

Alternative Herleitung aus der Definition des Trägheitsmoments

Sei ein rotierender Körper aus einzelnen Massepunkten aufgebaut. Für die Ortsvektoren [math]\vec r_i[/math] der Massepunkte wird ein kartesisches Koordinatensystem verwendet, dessen Ursprung mit dem Schwerpunkt des Körpers zusammenfällt. Sei [math]\vec r_n[/math] der Ortsvektor des [math]n[/math]-ten Massepunktes und [math]s_n[/math] sein Abstand von der Drehachse, dann gilt

[math]s_n=\sqrt{(\vec r_n)^2-\left(\frac{\vec \omega}{\omega}\cdot\vec r_n\right)^2}[/math]

und somit für das Trägheitsmoment

[math]I=\sum_{n}m_ns_n^2[/math]
[math]=\sum_{n}m_n\left((\vec r_n)^2-\left(\frac{\vec \omega}{\omega}\cdot\vec r_n\right)^2\right)[/math]
[math]=\frac{1}{\omega^2}\sum_{n}m_n\left((\vec \omega)^2(\vec r_n)^2-\left(\vec \omega\cdot\vec r_n\right)^2\right)[/math]
[math]=\frac{1}{\omega^2}\sum_{i,j=1}^3\left\{\sum_{n}m_n\left((\vec r_n)^2\delta_{ij}-r_{n,i}r_{n,j}\right)\right\}\omega_i\omega_j[/math]
[math]=\frac{1}{\omega^2}\sum_{i,j=1}^3 I_{ij} \; \omega_i\omega_j[/math]

mit dem Trägheitstensor

[math]I_{ij}=\sum_{n}m_n\left((\vec r_n)^2\delta_{ij}-r_{n,i}r_{n,j}\right).[/math]

Literatur


Kategorien: Keine Kategorien vorhanden!

Quelle: Wikipedia - http://de.wikipedia.org/wiki/Trägheitstensor (Vollständige Liste der Autoren des Textes [Versionsgeschichte])    Lizenz: CC-by-sa-3.0

Änderungen: Alle Bilder mit den meisten Bildunterschriften wurden entfernt. Ebenso alle zu nicht-existierenden Artikeln/Kategorien gehenden internen Wikipedia-Links (Bsp. Portal-Links, Redlinks, Bearbeiten-Links). Entfernung von Navigationsframes, Geo & Normdaten, Mediadateien, gesprochene Versionen, z.T. ID&Class-Namen, Style von Div-Containern, Metadaten, Vorlagen, wie lesenwerte Artikel. Ansonsten sind keine Inhaltsänderungen vorgenommen worden. Weiterhin kann es durch die maschinelle Bearbeitung des Inhalts zu Fehlern gerade in der Darstellung kommen. Darum würden wir jeden Besucher unserer Seite darum bitten uns diese Fehler über den Support mittels einer Nachricht mit Link zu melden. Vielen Dank!

Stand der Informationen: August 201& - Wichtiger Hinweis: Da die Inhalte maschinell von Wikipedia übernommen wurden, ist eine manuelle Überprüfung nicht möglich. Somit garantiert LinkFang.de nicht die Richtigkeit und Aktualität der übernommenen Inhalte. Sollten die Informationen mittlerweile fehlerhaft sein, bitten wir Sie darum uns per Support oder E-Mail zu kontaktieren. Wir werden uns dann innerhalb von spätestens 10 Tagen um Ihr Anliegen kümmern. Auch ohne Anliegen erfolgt mindestens alle drei Monate ein Update der gesamten Inhalte.