Brachistochrone

Experiment: Welche Bahn ist die schnellste? (Ausstellung Elementa im Landesmuseum für Technik und Arbeit, Mannheim)

Brachistochrone

Tautochronie der Brachistochrone – von jedem Startpunkt auf der Kurve erreichen die Kugeln das Ziel gleichzeitig.

Die Brachistochrone (griechisch brachistos kürzeste, chronos Zeit) ist die Bahn zwischen einem Anfangs- und einem gleich hoch oder tiefer gelegenen Endpunkt, auf der ein sich reibungsfrei bewegender Massenpunkt, der mit Geschwindigkeit null startet, unter dem Einfluss der Gravitationskraft am schnellsten zum Endpunkt gleitet. Der tiefste Punkt der Bahn kann tiefer liegen als der Endpunkt.

Der Körper gleitet auf einer solchen Bahn schneller zum Ziel als auf jeder anderen Bahn, auch wenn diese kürzer ist, beispielsweise geradlinig.

Wenn die Kurve am Ende waagrecht verläuft, also am tiefsten Punkt endet, dann ist sie zugleich eine Tautochrone, das heißt von jedem anderen Startpunkt auf der Kurve benötigt der Massepunkt die gleiche Zeit, um zum Endpunkt zu gelangen. Dieser Sachverhalt wird beim sogenannten Zykloidenpendel ausgenutzt, bei dem die Pendelmasse auf einer Tautochrone schwingt.

Form

Die Brachistochrone ist Teil einer Zykloide.

Geschichte

Johann I Bernoulli hat sich mit dem Problem des schnellsten Falles beschäftigt. Im Jahre 1696 fand er schließlich die Lösung in der Brachistochrone. Heute sieht man dies oft als die Geburtsstunde der Variationsrechnung.

Christiaan Huygens veröffentlichte 1673 in seiner Abhandlung Horologium Oscillatorium eine ganggenaue Pendeluhr mit einem Zykloidenpendel, bei dem er sich die Tatsache zunutze machte, dass die Evolute der Zykloide selbst wieder eine Zykloide ist. Der Vorteil der Ganggenauigkeit wird jedoch durch die erhöhte Reibung wett- bzw. zunichtegemacht.

Funktion

Die Brachistochrone lässt sich in einer Parameterdarstellung beschreiben, das heißt, man kann ihre Punkte als Ortsvektor darstellen, der sich mit einem Parameter ändert. Als Funktion des Winkels $\varphi$ (im Bogenmaß), um den sich das Rad mit Radius $R$ beim Abrollen gedreht hat, sind die $x$- und $y$-Koordinaten:

$x=R\cdot (\varphi -\sin \varphi )\,,$

$y=R\cdot (-1+\cos \varphi )\,.$

Hilfreich für das Verstehen dieser Kurve ist: Der Radius mal dem Winkel „Berührungspunkt des Kreises-Kreismittelpunkt-Brachistochronenpunkt“ ist die bereits abgerollte Strecke.

Herleitung

Betrachten wir in der $x$-$y$-Ebene eine Kurve $y(x)$, längs welcher der Massepunkt vom Start $(x,y)=(0,0)$ mit fortlaufender Zeit $t$ zum Ziel $({\overline {x}},{\overline {y}})$ gleite.

Er hat die kinetische Energie

$E_{\text{kin}}={\frac {1}{2}}\,m\,v^{2}={\frac {1}{2}}\,m\,({v_{x}}^{2}+{v_{y}}^{2})={\frac {1}{2}}\,m\,\left(\left({\frac {\mathrm {d} x}{\mathrm {d} t}}\right)^{2}+\left({\frac {\mathrm {d} y}{\mathrm {d} x}}\,{\frac {\mathrm {d} x}{\mathrm {d} t}}\right)^{2}\right)={\frac {1}{2}}\,m\,\left({\frac {\mathrm {d} x}{\mathrm {d} t}}\right)^{2}\left(1+\left({\frac {\mathrm {d} y}{\mathrm {d} x}}\right)^{2}\right)$

und die potentielle Energie

$E_{\text{pot}}=m\cdot g\cdot y(x)$

Dabei ist $y$ die Höhe im Gravitationsfeld und $g$ die Schwerebeschleunigung.

Gleitet der anfänglich ruhende Massepunkt vom Ursprung los, so ist längs seiner Bahn die Gesamtenergie erhalten und hat den anfänglichen Wert Null,

$0={\frac {1}{2}}\,m\,\left({\frac {\mathrm {d} x}{\mathrm {d} t}}\right)^{2}\left(1+\left({\frac {\mathrm {d} y}{\mathrm {d} x}}\right)^{2}\right)+m\cdot g\cdot y(x)$

Dies kann nach ${\tfrac {\mathrm {d} x}{\mathrm {d} t}}$ aufgelöst werden. Die Ableitung der Umkehrfunktion, $t(x)$, die angibt, zu welchem Zeitpunkt das Teilchen den Ort $(x,y(x))$ durchläuft, ist hierzu invers

${\frac {\mathrm {d} t}{\mathrm {d} x}}={\sqrt {\frac {1+({\frac {\mathrm {d} y}{\mathrm {d} x}})^{2}}{-2\,g\,y}}}$

Integrieren wir über den $x$-Bereich von 0 bis ${\overline {x}}$, so ergibt sich die zu minimierende Laufzeit $T$ als Funktional der Bahnkurve $y(x)$

$T[y]={\frac {1}{\sqrt {2\,g}}}\int _{0}^{\overline {x}}\,{\sqrt {\frac {1+({\frac {\mathrm {d} y}{\mathrm {d} x}})^{2}}{-y}}}\mathrm {d} x$

Um an die bei physikalischen Variationsproblemen üblichen Bezeichnungen anzuschließen, nennen wir die Integrationsvariable die "Zeit" $\tau$ (also ${\overline {\tau }}={\overline {x}}$), bezeichnen $-y$ mit $q$ und minimieren einfachheitshalber das mit ${\sqrt {2\,g}}$ multiplizierte Funktional.

Wir minimieren also die "Wirkung"

$W[q]=\int _{0}^{\overline {\tau }}\,{\sqrt {\frac {1+({\frac {\mathrm {d} q}{\mathrm {d} \tau }})^{2}}{q}}}\mathrm {d} \tau \quad$,

wobei die "Zeit"-Grenzen fest sind.

Der Integrand hat also die Rolle einer Lagrangefunktion:

${\mathcal {L}}(\tau ,q,v)={\sqrt {\frac {1+v^{2}}{q}}}\quad$

Da die Lagrangefunktion nicht vom Integrationsparameter, der "Zeit" $\tau$ abhängt, ist die nach dem Noether-Theorem zugehörige "Energie" / Hamilton-Funktion

$H=v\partial _{v}{\mathcal {L}}-{\mathcal {L}}=-{\frac {1}{\sqrt {q(1+v^{2})}}}$

auf der Bahn $q(\tau )$ erhalten, für die $W[q]$ minimal wird.

Die Funktion $q(\tau )$ erfüllt also mit einer positiven Konstanten $R$ die Gleichung

$\left(1+\left({\frac {\mathrm {d} q}{\mathrm {d} \tau }}\right)^{2}\right)\,q=2\,R$ oder $\left({\frac {\mathrm {d} q}{\mathrm {d} \tau }}\right)^{2}-{\frac {2\,R}{q}}=-1$

wie ein Teilchen, das im Keplerpotential $\propto -1/q$ senkrecht aus der Gipfelhöhe $2\,R$ fällt.

Statt diese Gleichung mit getrennten Veränderlichen nach ${\tfrac {\mathrm {d} q}{\mathrm {d} \tau }}$ aufzulösen und zu integrieren, bestätigt man einfach, dass

$\tau (\varphi )=R\,(\varphi -\sin \varphi )\,,\ q(\varphi )=R\,(1-\cos \varphi )$

eine parametrische Lösung dieser Gleichung ist, wobei man

${\frac {\mathrm {d} q}{\mathrm {d} \tau }}={\frac {\frac {\mathrm {d} q}{\mathrm {d} \varphi }}{\frac {\mathrm {d} \tau }{\mathrm {d} \varphi }}}={\frac {\sin \varphi }{1-\cos \varphi }}$

ausnutzt.

Also ist die gesuchte Bahn $(x,y(x))$ parametrisch gegeben durch

${\begin{pmatrix}x(\varphi )\\y(\varphi )\end{pmatrix}}=R\,{\begin{pmatrix}\varphi -\sin \varphi \\\cos \varphi -1\end{pmatrix}}=R\,{\begin{pmatrix}\varphi \\-1\end{pmatrix}}+{\begin{pmatrix}\cos \varphi &-\sin \varphi \\\sin \varphi &\cos \varphi \end{pmatrix}}{\begin{pmatrix}0\\R\end{pmatrix}}$

Dabei wird an der letzten Zerlegung deutlich, dass die Bahn $y(x)$ sich aus den Ortsvektoren $R\,(\varphi ,-1)$ der Nabe eines Rades mit Radius $R$ zusammensetzt, das unter der $x$-Achse rollt, plus dem Speichenvektor, der anfänglich nach oben zeigt und mit dem Winkel $\varphi$ gedreht wird. Die Kurve ist die Bahn eines Randpunktes eines rollenden Rades.

Spezielle Eigenschaften der Bahn

• Die Bahn ist unabhängig von der Masse und der Gewichtskraft des Körpers, also unabhängig von der Größe der Erdbeschleunigung.
• Ebenso ändert eine rollende Kugel, die Rotationsenergie aufnimmt, nichts an der Idealkurve.
• Die Tangente im Anfangspunkt ist senkrecht.
• Haben zwei Brachistochronen dasselbe Gefälle zwischen Anfangs- und Endpunkt, sind sie ähnlich.
• Ist das Gefälle nicht kleiner als 2/π (63,66 %), so ist der Endpunkt der tiefste Punkt der Kurve, bei kleinerem Gefälle liegt der Tiefpunkt zwischen Anfangs- und Endpunkt.
• Ist das Gefälle 0, also liegen Anfangs- und Endpunkt auf derselben Höhe, ist die Kurve symmetrisch.

Bilder

• 
  1. Phase: Beide Kugeln sind gleichauf.
• 
  2. Phase: Vordere Kugel beschleunigt durch stärkeres Gefälle.
• 
  3. Phase: Vordere Kugel liegt trotz längeren Wegs vorne.

Weblinks

Commons: Brachistochrone – Sammlung von Bildern, Videos und Audiodateien

• Brachistochrone Construction
• GeoGebra-modell zum interaktiven Ausprobieren mit zehn Stationen zwischen Start- und Endpunkt
• Cyclobahn, eine U-Bahn mit Zykloiden-Trasse
• Der Irrtum von Galileo Galilei - Vortrag von Rudolf Taschner
• Brachistochrone - Ableitung, Eigenschaften und lineare Approximation