Beschleunigtes Bezugssystem

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

Beschleunigte Bezugssysteme sind alle Bezugssysteme, die kein Inertialsystem sind.

Obwohl in beschleunigten Bezugssystemen die physikalischen Gesetze im Allgemeinen komplizierter aussehen (in der Mechanik müssen z. B. bei der Aufstellung von Bewegungsgleichungen Trägheitskräfte berücksichtigt werden), können diese Bezugssysteme in manchen Fällen dennoch die Lösung bestimmter Probleme vereinfachen.

Das ist meist dann der Fall, wenn sie an die Bewegungen angepasst sind, die Bewegungen im beschleunigten Bezugssystem also einfach werden:

  • Rotierende Kreis- oder Spiralbewegungen um ein gemeinsames Zentrum lassen sich z. B. in um das Zentrum gleichförmig rotierenden Bezugssystemen oft gut beschreiben: Der kreiselnde bzw. spiralende Körper ruht dann oder bewegt sich entlang einer Geraden.
  • Das Foucaultsche Pendel wird meist in einem Bezugssystem berechnet, das die Erddrehung mitvollführt.

Kinematik[Bearbeiten]

Inertialsystem K und beschleunigtes Koordinatensystem K'.

Der Ortsvektor eines Punktes P sei im Inertialsystem K gleich \vec r. Der Koordinatenursprung eines bewegten Bezugssystems K' sei von K aus betrachtet bei \vec R. Der Ortsvektor des Punktes P in K' sei \vec r\,'.

\vec{r}=\vec{R}+\vec{r}\,'

Im Folgenden sollen die zeitlichen Ableitungen, also Geschwindigkeit und Beschleunigung berechnet werden. Der Ortsvektor von P komponentenweise bzgl. K. Schreibe die Komponenten, die in K' dargestellt sind gestrichen.

\vec{r}=\sum_{i}x_{i}\hat{e}_{i}
\vec{r}\,'=\sum_{i}x'_{i}\hat{e}'_{i}

Nun wird zuerst die Zeitableitung für Vektoren, die im bewegten Bezugssystem dargestellt sind, berechnet. Da K' rotieren darf, ändern sich auch dessen Basisvektoren und müssen somit auch differenziert werden (Produktregel):

\frac{\mathrm{d}}{\mathrm{d}t}\vec{r}\,'=\sum_{i}\frac{\mathrm{d}x'_{i}}{\mathrm{d}t}\hat{e}'_{i}+\sum_{i}x'_{i}\frac{\mathrm{d}\hat{e}'_{i}}{\mathrm{d}t}

Die Länge der normierten Basisvektoren ändert sich nicht, deren Richtung kann sich aber von K aus gesehen ändern. Dem entspricht eine Rotation um die durch \vec{\omega}/\omega=:\hat{n} vorgegebene Achse mit der Winkelgeschwindigkeit \omega. Der Zusammenhang zwischen Inertialbasis und Basis des bewegten Bezugssystems ist gegeben durch die lineare (längen- und orientierungstreue) Abbildung:

\hat{e}'_{i}=\sum_{j}\left[R_{\vec{\omega}}\right]_{ij}\hat{e}_{j}

Die Komponenten der dreidimensionalen Drehmatrix R_{\vec{\omega}} sind dabei (mit Kronecker-Delta \delta_{ij} und Epsilon-Tensor \varepsilon_{ikj})

[R_{\vec{\omega}}]_{ij}=\left[1-\cos(\omega t)\right]n_{i}n_{j}+\cos(\omega t)\,\delta_{ij}+\sin(\omega t)\,[J_{\hat{n}}]_{ij}=\left[\exp\left(\omega t J_{\hat{n}}\right)\right]_{ij}

Die schiefsymmetrische Matrix J_{\hat{n}} ist dabei die Erzeugende der passiven Drehung um die \hat n-Achse.

[J_{\hat{n}}]_{ij}=\sum_{k}-\varepsilon_{ikj}n_{k}=\sum_{k}\varepsilon_{kij}n_{k}

Da die Basisvektoren gedreht werden, handelt es sich um eine passive Drehung; deswegen das negative Vorzeichen in der J-Matrix. Die zeitliche Ableitung der Drehmatrix lässt sich als Matrizenmultiplikation darstellen:

\frac{\mathrm{d}}{\mathrm{d}t}R_{\vec{\omega}}=\omega J_{\hat{n}}R_{\vec{\omega}}

Dabei entspricht \omega J_{\hat{n}} der Matrixdarstellung des Kreuzprodukts -\vec{\omega} \times.

\sum_{i}x'_{i}\frac{\mathrm{d}\hat{e}'_{i}}{\mathrm{d}t}=\sum_{i,j,l}x'_{i}\omega\left[J_{\hat{n}}\right]_{ij}\left[R_{\vec{\omega}}\right]_{jl}\hat{e}_{l}=\sum_{i,j,k}\varepsilon_{kij}\underbrace{\omega n_{k}}_{\omega_{k}}x'_{i}\hat{e}'_{j}=\vec{\omega}\times\vec{r}\,'

Somit lässt sich die zeitliche Ableitung der gestrichenen Basisvektoren in drei Dimensionen als Kreuzprodukt darstellen:

\frac{\mathrm{d}\hat{e}'_{i}}{\mathrm{d}t}=\vec{\omega}\times\hat{e}'_{i}

Wirkt eine Ableitung nur auf die gestrichen Komponenten von \vec r\,', so schreibe dafür die gestrichene Ableitung. Dies ist die zeitliche Änderung des Vektors \vec{r}\,', also die Geschwindigkeit \vec{v}\,', die man im bewegten Bezugssystem K' misst:

\frac{\mathrm{d}'}{\mathrm{d}t}\vec{r}\,'=\sum_{i}\frac{\mathrm{d}x'_{i}}{\mathrm{d}t}\hat{e}'_{i}=\vec{v}\,'

Zusammen ergibt sich also:

\frac{\mathrm{d}}{\mathrm{d}t}\vec{r}\,'=\sum_{i}\frac{\mathrm{d}x'_{i}}{\mathrm{d}t}\hat{e}'_{i}+\vec{\omega}\times\sum_{i}x'_{i}\hat{e}'_{i} =\vec{v}\,'+\vec{\omega}\times\vec{r}\,'

Somit entspricht einer Zeitableitung in K die gestrichene Zeitableitung in K' plus der Rotation von K':

\frac{\mathrm{d}}{\mathrm{d}t}=\left(\frac{\mathrm{d}'}{\mathrm{d}t}+\vec{\omega}\times\right)

Speziell für den Vektor \omega gilt:

\frac{\mathrm{d}}{\mathrm{d}t}\vec{\omega}=\frac{\mathrm{d}'}{\mathrm{d}t}\vec{\omega}+\underbrace{\vec{\omega}\times\vec{\omega}}_{=0}=\frac{\mathrm{d}'}{\mathrm{d}t}\vec{\omega}

Um die zweite Ableitung zu berechnen, verwende zweimal obige Relation:

\frac{\mathrm{d}^{2}}{\mathrm{d}t^{2}}\vec{r}\,'=\left(\frac{\mathrm{d}'}{\mathrm{d}t}+\vec{\omega}\times\right)\left(\frac{\mathrm{d}'}{\mathrm{d}t}+\vec{\omega}\times\right)\vec{r}\,'=\underbrace{\frac{\mathrm{d}'^{2}\vec{r}\,'}{\mathrm{d}t^{2}}}_{\vec{a}\,'}+2\vec{\omega}\times\underbrace{\frac{\mathrm{d}'\vec{r}\,'}{\mathrm{d}t}}_{\vec{v}\,'}+\underbrace{\frac{\mathrm{d}'\vec{\omega}}{\mathrm{d}t}}_{\dot{\vec{\omega}}}\times\vec{r}\,'+\vec{\omega}\times\left(\vec{\omega}\times\vec{r}\,'\right)
\frac{\mathrm{d}^{2}}{\mathrm{d}t^{2}}\vec{r}\,'=\vec{a}\,'+2\vec{\omega}\times\vec{v}\,'+\dot{\vec{\omega}}\times\vec{r}\,'+\vec{\omega}\times\left(\vec{\omega}\times\vec{r}\,'\right)


Für den Punkt P mit \vec{r}=\vec{R}+\vec{r}\,' erhält man schließlich die Geschwindigkeit und Beschleunigung:

\vec{v}=\frac{\mathrm{d}}{\mathrm{d}t}\vec{r}=\frac{\mathrm{d}}{\mathrm{d}t}\vec{R}+\frac{\mathrm{d}}{\mathrm{d}t}\vec{r}\,'=\dot{\vec{R}}+\vec{v}\,'+\vec{\omega}\times\vec{r}\,'
\vec{a}=\frac{\mathrm{d}^{2}}{\mathrm{d}t^{2}}\vec{r}=\frac{\mathrm{d}^{2}}{\mathrm{d}t^{2}}\vec{R}+\frac{\mathrm{d}^{2}}{\mathrm{d}t^{2}}\vec{r}\,'=\ddot{\vec{R}}+\vec{a}\,'+2\vec{\omega}\times\vec{v}\,'+\dot{\vec{\omega}}\times\vec{r}\,'+\vec{\omega}\times\left(\vec{\omega}\times\vec{r}\,'\right)

Dynamik[Bearbeiten]

Die newtonsche Bewegungsgleichung gilt nur für Inertialsysteme und lautet:

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

In allen Bezugssystemen ist eine externe Kraft gleich, diese muss nur bzgl. der jeweiligen Basisvektoren ausgedrückt werden: \vec{F}=\vec{F}\,'. Einsetzen in obige Gleichung für die Beschleunigung und Umstellen nach \vec{a}\,' liefert:

m\vec{a}\,'=\vec{F}\,'-m\ddot{\vec{R}}\underbrace{-2m\vec{\omega}\times\vec{v}\,'}_{\text{Coriolis-}}-m\dot{\vec{\omega}}\times\vec{r}\,'\underbrace{-m\vec{\omega}\times\left(\vec{\omega}\times\vec{r}\,'\right)}_{\text{Zentrifugalkraft}}

Somit erhält man diese vier Trägheitskräfte als zusätzliche Terme in der Bewegungsgleichung bzgl. K'. Im Folgenden werden verschiedene Spezialfälle diskutiert.

Beschleunigte Translationsbewegung[Bearbeiten]

Hier gilt \vec{\omega}=0, so dass sich die Bewegungsgleichung vereinfacht zu:

m\vec{a}\,'=\vec{F}\,'-m\ddot{\vec{R}}

Dies ist z. B. der Fall eines mit einem geradlinig bewegten Fahrzeug verbundenen Bezugssystems. Sei F\,'=0, also es soll keine externe Kraft wirken. Bremst das Fahrzeug, so ist die Beschleunigung negativ (= Verzögerung) \ddot{\vec{R}} \dot {\vec{R}}<0 und somit \vec{a}\,' \dot {\vec{R}}>0. Ein im Fahrzeug sich befindlicher Körper wird also in Fahrtrichtung beschleunigt (z. B. Autofahren: „Kopfnicker“ beim kurzen starken Bremsen).

Rotierendes Bezugssystem[Bearbeiten]

Es soll \ddot{\vec{R}}=0 gelten, d. h. der Ursprung von K' bewegt sich gleichförmig gegenüber dem Ursprung von K:

m\vec{a}\,'=\vec{F}\,'-2m\vec{\omega}\times\vec{v}\,'-m\dot{\vec{\omega}}\times\vec{r}\,'-m\vec{\omega}\times\left(\vec{\omega}\times\vec{r}\,'\right)

Bezugssystem an der Erdoberfläche[Bearbeiten]

Die Winkelgeschwindigkeit der Erde ist konstant, d. h. \dot{\vec{\omega}}=0. Hier rotiert \vec R (= Vektor vom Erdmittelpunkt zum Ursprung von K' an der Erdoberfläche) mit derselben Winkelgeschwindigkeit wie K':

m\vec{a}\,'=\vec{F}\,'-m\ddot{\vec{R}}-2m\vec{\omega}\times\vec{v}\,'-m\vec{\omega}\times\left(\vec{\omega}\times\vec{r}\,'\right)

Stellt man \vec{R} bzgl. K' dar, so ergibt die zweite Zeitableitung (\vec{R} ist bzgl. K' konstant):

\frac{\mathrm{d}^{2}}{\mathrm{d}t^{2}}\vec{R}=\underbrace{\frac{\mathrm{d}'^{2}\vec{R}}{\mathrm{d}t^{2}}}_{=0}+2\vec{\omega}\times\underbrace{\frac{\mathrm{d}'\vec{R}}{\mathrm{d}t}}_{=0}+\underbrace{\frac{\mathrm{d}'\vec{\omega}}{\mathrm{d}t}}_{=0}\times\vec{R}+\vec{\omega}\times\left(\vec{\omega}\times\vec{R}\right)=\vec{\omega}\times\left(\vec{\omega}\times\vec{R}\right)

Somit ergibt sich die Bewegungsgleichung:

m\vec{a}\,'=\vec{F}\,'-m\vec{\omega}\times\left(\vec{\omega}\times\vec{R}\right)-2m\vec{\omega}\times\vec{v}\,'-m\vec{\omega}\times\left(\vec{\omega}\times\vec{r}\,'\right)

Für Bewegungen, die in der Nähe der Erdoberfläche verlaufen, kann man den letzten Term vernachlässigen, da hier |\vec{r}\,'| \ll |\vec{R}| gilt.

Setze als Kraft die Gewichtskraft \vec{F}=m\vec{g} ein:

\vec{a}\,'=\underbrace{\vec{g}-\vec{\omega}\times\left(\vec{\omega}\times\vec{R}\right)}_{\vec{g}_{\text{eff}}}-2\vec{\omega}\times\vec{v}\,'-\vec{\omega}\times\left(\vec{\omega}\times\vec{r}\,'\right)

Man fasst normalerweise die Gravitationsbeschleunigung (\vec{g} wirkt in radiale Richtung) und die Zentrifugalbeschleunigung (-\vec{\omega}\times(\vec{\omega}\times\vec{R}) wirkt senkrecht zur Erdachse) zusammen zu einer effektiven Schwerebeschleunigung (die Richtung folgt aus der Vektorsummenbildung). Da die Zentrifugalbeschleunigung von der geographischen Breite abhängt (an den Polen Null und am Äquator maximal), ist die effektive Schwerebeschleunigung von der geographischen Breite abhängig; die Erdoberfläche ist näherungsweise eine Äquipotentialfläche der effektiven Schwerebeschleunigung, nämliche ein Ellipsoid, das im Vergleich zur Kugel an den Polen abgeplattet ist. \vec{g}_{\text{eff}} bestimmt die Vertikale von der Erdoberfläche, die von der radialen Richtung etwas abweicht.

Man betrachte ein mitbewegtes Koordinatensystem K' auf der Erdoberfläche, das so ausgerichtet ist, dass \hat{e}'_{x} in Richtung Osten, \hat{e}'_{y} in Richtung Norden und \hat{e}'_{z} zum Zenit zeigt. Die Winkelgeschwindigkeit der Erde lautet in K', wobei \phi die geographische Breite ist,

\vec{\omega}=\Omega\sin\phi\,\hat{e}'_{z}+\Omega\cos\phi\,\hat{e}'_{y}=:\Omega_{v}\,\hat{e}'_{z}+\Omega_{h}\,\hat{e}'_{y}

Somit lautet die Coriolisbeschleunigung

\vec{a}^{\,\prime}_{c}=-2\vec{\omega}\times\vec{v}\,'=-2\begin{pmatrix}0\\ \Omega_{h}\\ \Omega_{v}\end{pmatrix}\times\begin{pmatrix}\dot{x}'\\ \dot{y}'\\ \dot{z}'\end{pmatrix} =2\begin{pmatrix}\Omega_{v}\dot{y}'-\Omega_{h}\dot{z}'\\ -\Omega_{v}\dot{x}'\\ \Omega_{h}\dot{x}'\end{pmatrix}

Beispiel: Foucaultsches Pendel[Bearbeiten]

Für weitere Erklärung siehe Foucaultsches Pendel, hier soll explizit die Berechnung durchgeführt werden.

Man betrachte ein mathematisches Pendel auf der Erdoberfläche. Für kleine Auslenkungen (Auslenkung viel kleiner als Pendellänge) gilt die Näherung z'\approx\text{const} bzw. \dot{z}'\approx 0 und die Pendelmasse schwingt in der x'-y'-Ebene; somit kann man die Bewegung zweidimensional betrachten. Es wirkt auf die Pendelmasse m die rücktreibende Kraft \vec{F}_r=-m\omega_{0}^{2}\vec{x}^{\,\prime} mit \omega_{0}=\sqrt{g/l} und die Corioliskraft \vec{F}_c=2m\begin{pmatrix}\Omega_{v}\dot{y}'\\ -\Omega_{v}\dot{x}'\end{pmatrix}. Die Zwangskraft der Pendelaufhängung führt dazu, dass aus der eigentlich wirkenden Gewichtskraft \vec{F}_G=-mg\hat{e}_{z}^{\prime} die rücktreibende Kraft \vec{F}_r resultiert und dass die z'-Komponente der Corioliskraft kompensiert wird. Die Bewegungsgleichung der Pendelmasse lautet:

\begin{pmatrix}\ddot{x}'\\ \ddot{y}'\end{pmatrix}=-\omega_{0}^{2}\begin{pmatrix}x'\\ y'\end{pmatrix}+2\begin{pmatrix}\Omega_{v}\dot{y}'\\ -\Omega_{v}\dot{x}'\end{pmatrix}

Die zwei gekoppelten gewöhnlichen Differentialgleichungen zweiter Ordnung lassen sich einfach in der komplexen Darstellung lösen, wobei u:=x'+\mathrm{i}y' definiert wird.

\ddot{u}=-\omega_{0}^{2}u-\mathrm{i}2\Omega_{v}\dot{u}

Diese Differentialgleichung hat die allgemeine Lösung

u(t)=e^{-i\Omega_{v}t}\left(u_{1}e^{i\omega t}+u_{2}e^{-i\omega t}\right)

Der Term in runden Klammern beschreibt die "normale" Schwingung des Pendels mit der leicht verschobenen Frequenz \omega:=\sqrt{\omega_{0}^{2}+\Omega_{v}^{2}}\approx\omega_{0}, wobei im Allgemeinen die Bahnkurven der normalen zweidimensionalen Schwingung Ellipsen sind (je nach Anfangsbedingung, ist auch die Bewegung entlang einer Gerade möglich). Diese wird von einer weiteren Schwingung e^{-i\Omega_{v}t} überlagert, nämlich der Drehung der Schwingungsebene. Aus der Überlagerung beider Schwingungen kommen die bekannten Rosettenbahnen zustande. Mit \Omega=360°/Tag=15°/Stunde dreht sich bei einer geographische Breite \phi die Schwingungsebene um

-\Omega_v=-\Omega\,\sin\phi=-2\pi\,\sin\phi/\text{Tag}

auf der Nordhalbkugel (0\leq \phi \leq \pi/2 und 0\leq \sin\phi \leq 1) also im Uhrzeigersinn (in der Mathematik wird ein positiver Winkel im Gegenuhrzeigersinn definiert, negativer Winkel bedeutet also im Uhrzeigersinn). Für Deutschland mit ca. \phi=50° nördlicher Breite dreht sich die Schwingungsebene um etwa -11,5° pro Stunde.

Die Konstanten u_{1},u_{2}\in\mathbb{C} sind aus den Anfangsbedingungen x'(0)=x'_{0},\ y'(0)=y'_{0},\ \dot{x}'(0)=\dot{x}'_{0},\ \dot{y}'(0)=\dot{y}'_{0} zu bestimmen:

\begin{align}
u_{1} & =\frac{1}{2\omega}\left[(\omega+\Omega_{v})x'_{0}+\dot{y}'_{0}\right]+\frac{\mathrm{i}}{2\omega}\left[(\omega+\Omega_{v})y'_{0}-\dot{x}'_{0}\right]\\
u_{2} & =\frac{1}{2\omega}\left[(\omega-\Omega_{v})x'_{0}-\dot{y}'_{0}\right]+\frac{\mathrm{i}}{2\omega}\left[(\omega-\Omega_{v})y'_{0}+\dot{x}'_{0}\right]\end{align}

Literatur[Bearbeiten]