|
|
| Line 1: |
Line 1: |
| {{Scripthinweis|Mechanik|1}} | | {{Scripthinweis|Mechanik|1}} |
|
| |
|
| |
| = Das d'Alembertsche Prinzip =
| |
|
| |
|
| |
|
| |
| ==Zwangsbedingungen und Zwangskräfte==
| |
|
| |
| Ein System von N Massepunkten hat 3N Freiheitsgrade, wenn keine Zwangsbedingungen vorliegen. Die Zahl der Freiheitsgrade wird verringert durch
| |
|
| |
| ===Holonome (integrable) Zwangsbedingungen===
| |
|
| |
| Die Aufstellung der Zwangsbedingungen erfolgt derart, dass für eine Zwangsbedingung Lambda gilt:
| |
|
| |
|
| |
| <math>{{f}_{\lambda }}({{\vec{r}}_{1}},{{\vec{r}}_{2}},{{\vec{r}}_{3}},...{{\vec{r}}_{N}})=0\quad \quad \lambda =1,...\nu </math>
| |
|
| |
|
| |
| Betrachten wir als Beispiel einen Starren Körper aus 3 Teilchen, die jeweils den festen Abstand
| |
| [[Datei:Dreieck_graph.svg|miniatur|Punkte 1,2,3 mit Abständen <math>l_ij</math>]]
| |
|
| |
| <math>{{l}_{ij}}</math>
| |
| einhalten, so erhalten wir 3 Zwangsbedingungen:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{f}_{1}}=|{{{\vec{r}}}_{1}}-{{{\vec{r}}}_{2}}|-{{l}_{12}}=0 \\
| |
| & {{f}_{2}}=|{{{\vec{r}}}_{2}}-{{{\vec{r}}}_{3}}|-{{l}_{23}}=0 \\
| |
| & {{f}_{3}}=|{{{\vec{r}}}_{1}}-{{{\vec{r}}}_{3}}|-{{l}_{13}}=0 \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Die Zahl der Freiheitsgrade beträgt
| |
| <math>f=3N-\nu =9-3=6</math>
| |
|
| |
|
| |
| {| class="float-right wikitable"
| |
| |+ Starrer Körper aus N Teilchen
| |
| |- style="background: #DDFFDD;"
| |
| ! N
| |
| ! Hinzukommende
| |
| Einschränkungen
| |
| ! Zwangsbedingungen
| |
| (<math>\nu </math>)
| |
| ! Freiheitsgrade
| |
| <math>f=3N-\nu </math>
| |
| |-
| |
| ! style="background: #FFDDDD;"| 1
| |
| | 0
| |
| | 0
| |
| | 3
| |
| |-
| |
| ! style="background: #FFDDDD;"| 2
| |
| | 1
| |
| | 1
| |
| | 5
| |
| |-
| |
| ! style="background: #FFDDDD;"| 3
| |
| | 2
| |
| | 3
| |
| | 6
| |
| |-
| |
| ! style="background: #FFDDDD;"| 4
| |
| | 3
| |
| | 6
| |
| | 6
| |
| |-
| |
| ! style="background: #FFDDDD;"| 5
| |
| | 3
| |
| | 9
| |
| | 6
| |
| |-
| |
| ! style="background: #FFDDDD;"|...
| |
| | ..
| |
| | ..
| |
| | ..
| |
| |-
| |
| ! style="background: #FFDDDD;"|<math>N \ge 4</math>
| |
| | 3
| |
| | 3N-6
| |
| | 6
| |
| |}
| |
|
| |
| Allgemein könnte man nun für einen beliebigen starren Körper aus N Teilchen annehmen:
| |
|
| |
| <math>{{f}_{\lambda }}({{\vec{r}}_{1}},{{\vec{r}}_{2}},{{\vec{r}}_{3}},...{{\vec{r}}_{N}})=|{{\vec{r}}_{i}}-{{\vec{r}}_{j}}|-{{l}_{ij}}=0\quad i,j=1,2,...,N</math>
| |
| Jedoch sind diese Bedingungen nicht unabhängig. So gibt es für N=4 noch wie zu erwarten drei zusätzliche neue Einschränkungen. i,j kann von 1-4 laufen, 2 aus 4 sind gerade 6 Möglichkeiten und es gibt für N=4 auch genau 6 Zwangsbedingungen.
| |
|
| |
| Für N=5 kommen jedoch nicht 4 neue Zwangsbedingungen hinzu, sondern lediglich drei. Hier greift die Abhängigkeit einer Zwangsbedingung mit den anderen und man kann eine Zwangsbedingung der 2 aus 5 Kombinationen weglassen. Wäre dies nicht der Fall, so würde sich die Zahl der Freiheitsgrade des starren Körpers ja auch gerade wieder reduzieren, was unsinnig scheint.
| |
| Es zeigt sich, dass für jeden Massepunkt ab N=4 genau drei neue Einschränkungen hinzukommen. Die Zahl der Freiheitsgrade bleibt ab N=3 konstant, nämlich 6.
| |
| Unabhängigkeit bedeutet, dass für alle
| |
| <math>\lambda =1,...,\nu </math>
| |
| die Zwangsbedingungen ein linear unabhängiges Gleichungssystem bilden, also
| |
| <math>\operatorname{Rang}\left( \frac{\partial {{f}_{\lambda }}}{\partial {{r}_{i}}} \right)=\nu </math>
| |
| Somit gibt es genau 3N-6 Freiheitsgrade für N größer/ gleich drei:
| |
| Nun sucht man eine Lösung für die Bewegungsgleichung. Ohne Zwangsbedingung findet man für das i-te Teilchen eine Bahnkurve
| |
| <math>{{\vec{r}}_{i}}(t)</math>
| |
| . Alle Bahnen
| |
| <math>{{\vec{r}}_{i}}(t)</math>
| |
| müssen nun jedoch die
| |
| <math>\nu </math>
| |
| unabhängigen Zwangsbedingungen erfüllen:
| |
| <math>{{f}_{\lambda }}({{\vec{r}}_{1}}(t),{{\vec{r}}_{2}}(t),{{\vec{r}}_{3}}(t),...{{\vec{r}}_{N}}(t),t)=0\quad \quad \lambda =1,...\nu \quad f\ddot{u}r\ alle\ t</math>
| |
|
| |
|
| |
| Das totale Differenzial ( längs der Bahn
| |
| <math>{{\vec{r}}_{i}}(t)</math>
| |
| ) läßt sich schreiben:
| |
|
| |
|
| |
| <math>\frac{d{{f}_{\lambda }}}{dt}=\sum\limits_{i=1}^{N}{{{\nabla }_{ri}}{{f}_{\lambda }}\cdot {{{\vec{v}}}_{i}}+\frac{\partial {{f}_{\lambda }}}{\partial t}=}0\quad \quad \lambda =1,...\nu </math>
| |
|
| |
|
| |
| In differenzieller Schreibweise gewinnen wir das vollständige Differential:
| |
|
| |
|
| |
| <math>d{{f}_{\lambda }}=\sum\limits_{i=1}^{N}{{{\nabla }_{ri}}{{f}_{\lambda }}\cdot d{{{\vec{r}}}_{i}}+\frac{\partial {{f}_{\lambda }}}{\partial t}dt=}0\quad \quad \lambda =1,...\nu </math>
| |
|
| |
| ===Nichtholonome Zwangsbedingungen===
| |
| Nun sind jedoch Nichtholonome Zwangsbedingungen der Art:
| |
|
| |
|
| |
| <math>\sum\limits_{i=1}^{N}{{{{\vec{a}}}_{\lambda i}}({{{\vec{r}}}_{1}},{{{\vec{r}}}_{2}},...,{{{\vec{r}}}_{N}},t)\cdot d{{{\vec{r}}}_{i}}+{{{\vec{a}}}_{\lambda 0}}({{{\vec{r}}}_{1}},{{{\vec{r}}}_{2}},...,{{{\vec{r}}}_{N}},t)dt=}0\quad \quad \lambda =1,...\nu </math>
| |
|
| |
|
| |
| Dies ist eine Pfaffsche Differenzialform. Diese ist nicht integrabel, was gleichbedeutend ist damit, dass kein integrierender Faktor
| |
| <math>{{g}_{\lambda }}</math>
| |
| existiert, so dass
| |
|
| |
|
| |
| <math>\sum\limits_{i=1}^{N}{{{g}_{\lambda }}{{{\vec{a}}}_{\lambda i}}({{{\vec{r}}}_{1}},{{{\vec{r}}}_{2}},...,{{{\vec{r}}}_{N}},t)\cdot d{{{\vec{r}}}_{i}}+{{g}_{\lambda }}{{{\vec{a}}}_{\lambda 0}}({{{\vec{r}}}_{1}},{{{\vec{r}}}_{2}},...,{{{\vec{r}}}_{N}},t)dt=}d{{f}_{\lambda }}\quad \quad \lambda =1,...\nu </math>
| |
|
| |
|
| |
| Gleichbedeutend mit
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{g}_{\lambda }}{{{\vec{a}}}_{\lambda i}}={{\nabla }_{ri}}{{f}_{\lambda }} \\
| |
| & {{g}_{\lambda }}{{{\vec{a}}}_{\lambda 0}}=\frac{\partial {{f}_{\lambda }}}{\partial t} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Dies kann man wieder so interpretieren, dass beliebige Positionen der Teilchen, also
| |
| <math>{{\vec{r}}_{i}}(t)\quad i=1...N</math>
| |
| möglich sind, also
| |
| <math>{{\vec{r}}_{1}},{{\vec{r}}_{2}},...{{\vec{r}}_{N}}</math>
| |
| beliebig, jedoch ist die momentane Bewegungsrichtung eingeschränkt. Man sagt auch, die lokalen Bewegungen sind eingeschränkt ( längs der Bahn
| |
| <math>{{\vec{r}}_{i}}(t)\quad i=1...N</math>
| |
| )
| |
|
| |
|
| |
| <math>\sum\limits_{i=1}^{N}{{{{\vec{a}}}_{\lambda i}}({{{\vec{r}}}_{1}},{{{\vec{r}}}_{2}},...,{{{\vec{r}}}_{N}},t)\cdot {{{\vec{v}}}_{i}}+{{{\vec{a}}}_{\lambda 0}}({{{\vec{r}}}_{1}},{{{\vec{r}}}_{2}},...,{{{\vec{r}}}_{N}},t)=}0\quad \quad \lambda =1,...\nu </math>
| |
|
| |
|
| |
| Beispiel ist das Rangieren eines Autos auf einer freien Fläche. Jeder Punkt ist erreichbar, jedoch ist
| |
| <math>d{{\vec{r}}_{i}}(t)</math>
| |
| durch die momentane Radrichtung bestimmt
| |
|
| |
| Es ist weiter zu unterscheiden
| |
|
| |
| ====Zeitabhängigkeit====
| |
|
| |
| * zeitabhängige Zwangsbedingungen heißen '''rheonom'''
| |
| * zeitunabhängige ( nicht explizit zeitabhängige) , starre, ZB heißen '''skleronom'''
| |
|
| |
| ====Zwangsbedingungen als Ungleichungen====
| |
|
| |
| z.B. bei einem Gas in einem Behälter mit Wänden
| |
|
| |
| <u>'''Bewegungsgleichungen'''</u>
| |
|
| |
|
| |
| <math>{{m}_{i}}{{\ddot{\vec{r}}}_{i}}(t)={{\vec{F}}_{i}}+\sum\limits_{j}{{{{\vec{F}}}_{ij}}=:{{{\vec{X}}}_{i}}}\quad i=1...N</math>
| |
|
| |
|
| |
| diese Art ist bekannt. Auf der rechten Seite findet sich die Summe der '''Äußeren Kräfte''', eine äußere Kraft auf das i-te Teilchen und die Summe über die '''inneren Kräfte''' durch Wechselwirkung mit den weiteren j Teilchen, die anwesend sind. Die Summe aller Kräfte nennt man '''eingeprägte Kräfte'''.
| |
|
| |
| Diese Bewegungsgleichungen sind nun jedoch unter den '''Nebenbedingungen'''
| |
|
| |
|
| |
| <math>{{f}_{\lambda }}({{\vec{r}}_{1}},{{\vec{r}}_{2}},{{\vec{r}}_{3}},...{{\vec{r}}_{N}})=0\quad \quad \lambda =1,...\nu </math> (holonom)
| |
|
| |
| oder
| |
|
| |
| <math>\sum\limits_{i=1}^{N}{{{{\vec{a}}}_{\lambda i}}({{{\vec{r}}}_{1}},{{{\vec{r}}}_{2}},...,{{{\vec{r}}}_{N}},t)\cdot {{{\vec{v}}}_{i}}+{{{\vec{a}}}_{\lambda 0}}({{{\vec{r}}}_{1}},{{{\vec{r}}}_{2}},...,{{{\vec{r}}}_{N}},t)=}0\quad \quad \lambda =1,...\nu </math> (anholonom)
| |
|
| |
| zu lösen.
| |
|
| |
| Dazu soll die '''Beschreibung gewechselt''' werden.
| |
|
| |
| Wir nehmen an, dass die Nebenbedingungen (Zwangsbedingungen) durch Zwangskräfte
| |
| <math>{{\vec{Z}}_{i}}</math>
| |
| erzwungen werden.
| |
|
| |
| Damit folgt für unsere Bewegungsgleichung:
| |
|
| |
|
| |
| <math>{{m}_{i}}{{\ddot{\vec{r}}}_{i}}(t)={{\vec{F}}_{i}}+\sum\limits_{j}{{{{\vec{F}}}_{ij}}+{{{\vec{Z}}}_{i}}=:{{{\vec{X}}}_{i}}+{{{\vec{Z}}}_{i}}}\quad i=1...N</math>
| |
|
| |
|
| |
| Beim Beispiel der schiefen Ebene wirkt die Zwangskraft gerade der Normalkraft entgegen und verhindert somit das Fallen des Körpers durch die schiefe Ebene.
| |
|
| |
| [[Datei:Rownia.svg|miniatur|schiefe Ebene mit der Zwangskraft Z im Bild mit N Bezeichnet und Schwerkraft G=mg]]
| |
|
| |
| Es gilt:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \vec{Z}=mg\cos \vartheta \cdot \left( \begin{matrix}
| |
| \sin \vartheta \\
| |
| \cos \vartheta \\
| |
| \end{matrix} \right) \\
| |
| & \vec{F}=-mg\cdot \left( \begin{matrix}
| |
| 0 \\
| |
| 1 \\
| |
| \end{matrix} \right) \\
| |
| & \vec{Z}+\vec{F}=mg\sin \vartheta \cdot \left( \begin{matrix}
| |
| \cos \vartheta \\
| |
| -\sin \vartheta \\
| |
| \end{matrix} \right) \\
| |
| \end{align}</math>
| |
|
| |
| ==Virtuelle Verrückungen==
| |
| {{Def|
| |
| Unter einer virtuellen Verrückung
| |
| <math>\left\{ \delta {{{\vec{r}}}_{i}} \right\}</math>
| |
| versteht man die infinitesimale Änderung der Koordinaten, di zu fester Zeit
| |
| <math>\left\{ \delta t=0 \right\}</math>
| |
| die holonomen, bzw. nicht holonomen Zwangsbedingungen erfüllen.|virtuelle Verrückung
| |
| }}
| |
| Damit ist der Unterschied zu einer reellen Verrückung klar, die als
| |
| <math>d{{\vec{r}}_{i}}</math>
| |
| im Zeitintervall
| |
| <math>dt</math>
| |
| längs der Bahn geschieht.
| |
|
| |
| Die Zwangsbedingungen lassen sich jedoch nicht virtuell verrücken.
| |
|
| |
| Es gilt folglich
| |
|
| |
|
| |
| <math>\delta {{f}_{\lambda }}=\sum\limits_{i=1}^{N}{{{\nabla }_{ri}}{{f}_{\lambda }}\cdot \delta {{{\vec{r}}}_{i}}=}0\quad \quad \lambda =1,...\nu </math>
| |
|
| |
|
| |
|
| |
| <math>\sum\limits_{i=1}^{N}{{{{\vec{a}}}_{\lambda i}}({{{\vec{r}}}_{1}},{{{\vec{r}}}_{2}},...,{{{\vec{r}}}_{N}},t)\cdot \delta {{{\vec{r}}}_{i}}=}0\quad \quad \lambda =1,...\nu </math>
| |
|
| |
|
| |
| Die zeitabhängigen Anteile fallen raus, da ja nach Definition
| |
| <math>\left\{ \delta t=0 \right\}</math>
| |
| .
| |
|
| |
| Als Beispiel betrachten wir die Bewegung eines Massepunktes in einer Ebene:
| |
|
| |
|
| |
| <math>\vec{a}\cdot (\vec{r}-{{\vec{r}}_{o}}(t))=0</math>
| |
|
| |
|
| |
| Dabei ist
| |
| <math>{{\vec{r}}_{o}}(t)</math>
| |
| der Startpunkt des Teilchens, also ein fester Punkt in der Ebene und nicht notwendigerweise zeitunabhängig. a charakterisiert den Normalenvektor auf der Ebene Schließlich kann sich die Ebene bewegen, beispielsweise hoch und runter.
| |
|
| |
| Formuliert man nun holonome Zwangsbedingungen für N Massepunkte, so gilt:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & f({{{\vec{r}}}_{i}})=\vec{a}\cdot ({{{\vec{r}}}_{i}}-{{{\vec{r}}}_{o}}(t))=0\quad i=1,2,...,N \\
| |
| & df({{{\vec{r}}}_{i}})=\vec{a}\cdot (d{{{\vec{r}}}_{i}}-{{{\vec{v}}}_{o}}(t)dt)=0\quad i=1,2,...,N\quad {{{\vec{v}}}_{o}}(t)=\frac{d{{{\vec{r}}}_{o}}(t)}{dt} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| also gilt im Allgemeinen:
| |
|
| |
|
| |
| <math>\vec{a}\cdot d{{\vec{r}}_{i}}=\vec{a}\cdot {{\vec{v}}_{o}}(t)dt\ne 0\quad i=1,2,...,N\quad {{\vec{v}}_{o}}(t)=\frac{d{{{\vec{r}}}_{o}}(t)}{dt}</math>
| |
|
| |
|
| |
| aber:
| |
|
| |
|
| |
| <math>\delta f=\vec{a}\cdot \delta {{\vec{r}}_{i}}=0\quad i=1,2,...,N\quad {{\vec{v}}_{o}}(t)=\frac{d{{{\vec{r}}}_{o}}(t)}{dt}</math>
| |
|
| |
|
| |
| Das heißt, die virtuellen Verrückungen geschehen alle bei festgehaltenem
| |
| <math>{{\vec{r}}_{o}}(t)</math>
| |
| . Es gilt:
| |
| <math>\delta {{\vec{r}}_{i}}\bot \vec{a}</math>
| |
|
| |
| == D´Alembertsches Prinzip der virtuellen Arbeit==
| |
|
| |
| Gegeben sei ein System von N Massepunkten mit beliebigen ( holonomen oder nicht holonomen) Zwangsbed.
| |
|
| |
| Schreiben wir die Bewegungsgleichungen mit den Zwangskräften Zi als:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{m}_{i}}{{{\ddot{\vec{r}}}}_{i}}(t)-{{{\vec{X}}}_{i}}={{{\vec{Z}}}_{i}}\quad i=1...N \\
| |
| & \to \sum\limits_{i}{\left( {{m}_{i}}{{{\ddot{\vec{r}}}}_{i}}(t)-{{{\vec{X}}}_{i}} \right)\delta {{{\vec{r}}}_{i}}=}\sum\limits_{i}{{{{\vec{Z}}}_{i}}\delta {{{\vec{r}}}_{i}}} \\
| |
| \end{align}</math>
| |
|
| |
| {{Def|;Dabei versteht man
| |
| :<math>\sum\limits_{i}{{{{\vec{X}}}_{i}}\delta {{{\vec{r}}}_{i}}}</math> als virtuelle Arbeit der eingeprägten Kräfte und
| |
| :<math>\sum\limits_{i}{{{{\vec{Z}}}_{i}}\delta {{{\vec{r}}}_{i}}}</math> als virtuelle Arbeit der Zwangskräfte.|Virtuelle Arbeit}}
| |
|
| |
| {{Beispiel|<u>'''Beispiel: Bewegung auf einer Fläche'''</u>
| |
|
| |
|
| |
| <math>f({{\vec{r}}_{i}},t)=0</math>
| |
|
| |
|
| |
| das ist auf der Ebene gerade durch die Normale auszudrücken:
| |
|
| |
|
| |
| <math>\vec{a}\cdot (\vec{r}-{{\vec{r}}_{o}}(t))=0</math>
| |
|
| |
|
| |
| Annahme: Alle Zwangskräfte stehen senkrecht auf die Fläche:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{{\vec{Z}}}_{i}}={{\lambda }_{i}}({{{\vec{r}}}_{1}},{{{\vec{r}}}_{2}},...,{{{\vec{r}}}_{N}}){{\nabla }_{ri}}f \\
| |
| & {{\nabla }_{ri}}f\quad z.B.\vec{a}\ f\ddot{u}r\ Ebene \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Die Virtuelle Arbeit der Zwangskräfte verschwindet nun:
| |
|
| |
|
| |
| <math>{{\vec{Z}}_{i}}\delta {{\vec{r}}_{i}}=0={{\lambda }_{i}}({{\vec{r}}_{1}},{{\vec{r}}_{2}},...,{{\vec{r}}_{N}}){{\nabla }_{ri}}f\delta {{\vec{r}}_{i}}={{\lambda }_{i}}\delta f</math>
| |
|
| |
|
| |
| Begründung:
| |
|
| |
|
| |
| <math>{{\nabla }_{ri}}f\delta {{\vec{r}}_{i}}</math>
| |
| ist als Variation der Zwangsbedingung zu verstehen:
| |
|
| |
|
| |
| <math>{{\nabla }_{ri}}f</math>
| |
| ist ein Differenzial senkrecht auf die Fläche
| |
|
| |
|
| |
| <math>f\delta {{\vec{r}}_{i}}</math>
| |
| ein Differenzial parallel zur Fläche
| |
|
| |
| Also folgt:
| |
|
| |
|
| |
| <math>\sum\limits_{i}{{{{\vec{Z}}}_{i}}\delta {{{\vec{r}}}_{i}}}=0</math>
| |
|
| |
|
| |
| Die reale Arbeit der Zwangskräfte verschwindet dagegen im Allgemeinen nicht:
| |
|
| |
|
| |
| <math>\sum\limits_{i}{{{{\vec{Z}}}_{i}}d{{{\vec{r}}}_{i}}}\ne 0</math>
| |
| }}
| |
| {{Beispiel|
| |
| '''Beispiel: Starrer Körper'''
| |
|
| |
|
| |
| <math>{{f}_{\lambda }}=\left| {{{\vec{r}}}_{i}}-{{{\vec{r}}}_{j}} \right|-{{l}_{ij}}:={{r}_{ij}}-{{l}_{ij}}=0</math>
| |
|
| |
|
| |
| Annahme: Die Zwangskräfte wirken in Richtung
| |
| <math>{{\vec{r}}_{i}}-{{\vec{r}}_{j}}</math>
| |
|
| |
|
| |
|
| |
| <math>{{\vec{Z}}_{ij}}={{\lambda }_{ij}}\frac{{{{\vec{r}}}_{i}}-{{{\vec{r}}}_{j}}}{{{r}_{ij}}}</math>
| |
|
| |
|
| |
| Das Vorgehen läßt sich also folgendermaßen schematisieren:
| |
|
| |
| Bestimme die Richtung der Zwangskraft und multipliziere einen beliebigen skalaren Faktor mit dieser Richtung.
| |
|
| |
| Falls die Richtungen für verschiedene Zwangskräfte verschieden sind, so muss man diese indizieren ( mit einem Index kenntlich machen). Die Zwangskräfte erhalten dann ebenso indizierte skalare Faktoren.
| |
|
| |
| Mit Hilfe des 3. Newtonschen Axioms können wir weiter einschränken:
| |
|
| |
|
| |
| <math>{{\vec{Z}}_{ij}}=-{{\vec{Z}}_{ji}}\Rightarrow {{\lambda }_{i{{j}_{{}}}}}={{\lambda }_{ji}}</math>
| |
|
| |
|
| |
| Auf das Teilchen i wirkt also insgesamt die Zwangskraft:
| |
|
| |
|
| |
| <math>{{\vec{Z}}_{i}}=\sum\limits_{j\ne i}{{{Z}_{ij}}}=\sum\limits_{j}{{{\lambda }_{ij}}\frac{{{{\vec{r}}}_{i}}-{{{\vec{r}}}_{j}}}{{{r}_{ij}}}}</math>
| |
|
| |
|
| |
|
| |
| <math>{{\vec{Z}}_{i}}\delta {{\vec{r}}_{i}}\ne 0</math>
| |
| im Allgemeinen. Es verschwindet also nicht die virtuelle Arbeit für jede Masse einzeln.
| |
|
| |
| Jedoch gilt:
| |
|
| |
|
| |
| <math>\sum\limits_{i}{{{{\vec{Z}}}_{i}}{{\delta }_{{}}}{{{\vec{r}}}_{i}}}=\sum\limits_{i,j}{{{\lambda }_{ij}}\frac{{{{\vec{r}}}_{i}}-{{{\vec{r}}}_{j}}}{{{r}_{ij}}}{{\delta }_{{}}}{{{\vec{r}}}_{i}}}=\frac{1}{2}\sum\limits_{i,j}{{{\lambda }_{ij}}\frac{{{{\vec{r}}}_{i}}-{{{\vec{r}}}_{j}}}{{{r}_{ij}}}}{{\delta }_{{}}}{{({{\vec{r}}_{i}}-{{\vec{r}}_{j}})}_{{}}}=\frac{1}{2}\sum\limits_{i,j}{{{\lambda }_{ij}}}{{\delta }_{{}}}{{r}_{ij}}=0</math>
| |
|
| |
|
| |
| Beweis:
| |
|
| |
|
| |
| <math>\delta |r|=\delta {{(\vec{r}\cdot \vec{r})}^{\frac{1}{2}}}=\frac{1}{2}{{(\vec{r}\cdot \vec{r})}^{-\frac{1}{2}}}2\vec{r}\delta \vec{r}=\frac{\vec{r}\delta \vec{r}}{r}</math>
| |
|
| |
|
| |
| und
| |
|
| |
|
| |
| <math>{{\delta }_{{}}}{{r}_{ij}}=0</math>
| |
|
| |
| }}
| |
| Allgemein kann man fordern:
| |
|
| |
|
| |
| <math>\sum\limits_{i}{{{{\vec{Z}}}_{i}}{{\delta }_{{}}}{{{\vec{r}}}_{i}}}=0</math>
| |
| für alle betrachteten Zwangskräfte.
| |
|
| |
| Das bedeutet: Gleitreibungskräfte längs einer Fläche sind als Zwangskräfte ausgeschlossen.
| |
|
| |
| Somit folgt als d'Alembertsches Prinzip:
| |
|
| |
|
| |
| <math>\sum\limits_{i}{{{{\vec{Z}}}_{i}}{{\delta }_{{}}}{{{\vec{r}}}_{i}}}=\sum\limits_{i}{\left( {{m}_{i}}{{{\ddot{\vec{r}}}}_{i}}-{{{\vec{X}}}_{i}} \right)}\delta {{\vec{r}}_{i}}=0</math>
| |
|
| |
|
| |
| Das d´Alembertsche Prinzip gilt gleichermaßen für holonome und anholonome Zwangsbedingungen
| |
|
| |
| <u>'''Beispiel für ein Variationsprinzip:'''</u>
| |
|
| |
| '''Differentialprinzip: ( für infinitesimal kleine Variationen):'''
| |
|
| |
| Der wirklich angenommene Zustand eines Systems ist in Extremalzustand in dem Sinn, dass die gesamte virtuelle Arbeit Null ist. Dieser Zustand ist stabil gegen kleine Verrückungen der Bahn
| |
| <math>\left\{ \delta {{{\vec{r}}}_{i}} \right\}</math>
| |
| .
| |
|
| |
| ====Variationsprinzip mit Nebenbedingungen====
| |
|
| |
| Wir numerieren nun die Vektorkoordinaten um:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \vec{r}\to {{r}_{j}}(j=1...3) \\
| |
| & \vec{X}\to {{X}_{j}} \\
| |
| & \vec{a}\to {{b}_{j}}^{n} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Aus dem d´Alembertschen Prinzip gewinnen wir:
| |
|
| |
|
| |
| <math>\sum\limits_{i=1}^{3N}{{{Z}_{i}}{{\delta }_{{}}}{{r}_{i}}}=\sum\limits_{i=1}^{3N}{\left( {{m}_{i}}{{{\ddot{r}}}_{i}}-{{X}_{i}} \right)\delta {{r}_{i}}=0}</math>
| |
|
| |
|
| |
| Nebenbedingung:
| |
|
| |
|
| |
| <math>\sum\limits_{i=1}^{3N}{{{b}_{i}}^{n}\delta {{r}_{i}}=0\quad n=1,...,\nu }</math>
| |
|
| |
|
| |
| Nü charakterisiert auch hier die Zahl der Nebenbedingungen, der Index n steht für die n-te Nebenbedingung
| |
|
| |
| Dies ist lösbar mit der Methode der Lagrange-Multiplikatoren.
| |
|
| |
| Denn: Wenn die Vektorkomponenten
| |
| <math>{{r}_{i}}</math>
| |
| frei variierbar wären, also
| |
| <math>\delta {{r}_{i}}</math>
| |
| beliebig, so müsste gelten:
| |
|
| |
|
| |
| <math>{{m}_{i}}{{\ddot{r}}_{i}}-{{X}_{i}}=0</math>
| |
|
| |
|
| |
| Also wäre es sinnvoll, das lineare Gleichungssystem so umzuschreiben, dass ein Satz von Faktoren frei variierbar ist:
| |
|
| |
| Zuerst addieren wir die Nebenbedingungen mit noch beliebigen Lagrangemultiplikatoren
| |
| <math>{{\lambda }_{n}}</math>
| |
| :
| |
|
| |
| Wir erhalten:
| |
|
| |
|
| |
| <math>\sum\limits_{j=1}^{3N}{\left( {{m}_{j}}{{{\ddot{r}}}_{j}}-{{X}_{j}}-\sum\limits_{n=1}^{\nu }{{{\lambda }_{n}}{{b}_{j}}^{n}} \right)\delta {{r}_{j}}=0}</math>
| |
|
| |
|
| |
| Nun sind
| |
| <math>\delta {{r}_{1}},\delta {{r}_{2}},...,\delta {{r}_{\nu }}</math>
| |
| aus den Nebenbedingungen zu eliminieren.
| |
|
| |
| Die verbleibenden
| |
| <math>\delta {{r}_{\nu +1}},...,\delta {{r}_{3N}}</math>
| |
| sind nun frei variierbar.
| |
|
| |
| Nun kann das Summenzeichen weggelassen werden, da die verbleibenden Vektorkomponenten frei variiert werden können und dementsprechend jeder Summand für sich Null sein muss:
| |
|
| |
| Es lassen sich
| |
| <math>{{\lambda }_{1,}}...,{{\lambda }_{\nu }}</math>
| |
| derart bestimmen, dass
| |
|
| |
|
| |
| <math>{{m}_{j}}{{\ddot{r}}_{j}}-{{X}_{j}}-\sum\limits_{n=1}^{\nu }{{{\lambda }_{n}}{{b}_{j}}^{n}}=0\quad j=1,...,\nu </math>
| |
|
| |
|
| |
| Das heißt, wir suchen die
| |
| <math>{{\lambda }_{1,}}...,{{\lambda }_{\nu }}</math>
| |
| aus diesem gegebenen linearen Gleichungssystem für die
| |
| <math>{{\lambda }_{n}}(t)</math>
| |
| als Funktion der
| |
| <math>{{\ddot{r}}_{j}}(t)</math>
| |
| . Im stationären Fall ist dies direkt auflösbar.
| |
|
| |
|
| |
| <math>\sum\limits_{j=\nu +1}^{3N}{\left( {{m}_{j}}{{{\ddot{r}}}_{j}}-{{X}_{j}}-\sum\limits_{n=1}^{\nu }{{{\lambda }_{n}}{{b}_{j}}^{n}} \right)\delta {{r}_{j}}=0}</math>
| |
|
| |
|
| |
| Da hier jedoch die
| |
| <math>\delta {{r}_{j}}</math>
| |
| frei variierbar sind, gilt:
| |
|
| |
|
| |
| <math>{{m}_{j}}{{\ddot{r}}_{j}}-{{X}_{j}}-\sum\limits_{n=1}^{\nu }{{{\lambda }_{n}}{{b}_{j}}^{n}}=0</math>
| |
|
| |
|
| |
| Die Lagrange- Gleichung der 1. Art
| |
|
| |
|
| |
| <math>\sum\limits_{n=1}^{\nu }{{{\lambda }_{n}}{{b}_{j}}^{n}}</math>
| |
| kann als Zwangskraft interpretiert werden und taucht in der Statik als Lagrange- Parameter auf.
| |
|
| |
| ====Beispiel Atwoodsche Fallmaschine====
| |
|
| |
| [[Datei:Atwoods_machine_functionally.svg|miniatur|Atwoods Fallmaschine]]
| |
|
| |
| Aus der Schule bekannt ist die Kraft, die an m1 angreift, nämlich -m1g und die Kraft , die an m2 angreift, nämlich -m2g.
| |
|
| |
| Beginnen wir mit dem d´Alembertschen Prinzip:
| |
|
| |
|
| |
| <math>\sum\limits_{i}{{{{\vec{Z}}}_{i}}{{\delta }_{{}}}{{{\vec{r}}}_{i}}}=\sum\limits_{i}{\left( {{m}_{i}}{{{\ddot{\vec{r}}}}_{i}}-{{{\vec{X}}}_{i}} \right)}\delta {{\vec{r}}_{i}}=0</math>
| |
|
| |
|
| |
| so folgt:
| |
|
| |
|
| |
| <math>({{m}_{1}}{{\ddot{h}}_{1}}-{{X}_{1}})\delta {{h}_{1}}+({{m}_{2}}{{\ddot{h}}_{2}}-{{X}_{2}})\delta {{h}_{2}}=0</math>
| |
|
| |
|
| |
| Da der Aufbau nur ein Rädchen besitzt gilt ganz einfach:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{h}_{1}}+{{h}_{2}}=const. \\
| |
| & \delta {{h}_{1}}=-\delta {{h}_{2}} \\
| |
| & {{{\ddot{h}}}_{1}}=-{{{\ddot{h}}}_{2}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Also folgt:
| |
|
| |
|
| |
| <math>({{m}_{1}}{{\ddot{h}}_{1}}+{{m}_{1}}g)\delta {{h}_{1}}-(-{{m}_{2}}{{\ddot{h}}_{1}}+{{m}_{2}}g)\delta {{h}_{1}}=0</math>
| |
|
| |
|
| |
|
| |
| <math>{{m}_{1}}{{\ddot{h}}_{1}}+{{m}_{1}}g+{{m}_{2}}{{\ddot{h}}_{1}}-{{m}_{2}}g=0</math>
| |
|
| |
|
| |
|
| |
| <math>{{\ddot{h}}_{1}}=\frac{({{m}_{2}}-{{m}_{1}})}{{{m}_{1}}+{{m}_{2}}}g</math>
| |
|
| |
|
| |
| Also: Am bedeutendsten ist das d´Alembertsche Prinzip, welches sagt, dass die Summe über alle virtuellen Arbeiten der Zwangskräfte Null ist:
| |
|
| |
|
| |
| <math>\sum\limits_{i}{{{{\vec{Z}}}_{i}}{{\delta }_{{}}}{{{\vec{r}}}_{i}}}=\sum\limits_{i}{\left( {{m}_{i}}{{{\ddot{\vec{r}}}}_{i}}-{{{\vec{X}}}_{i}} \right)}\delta {{\vec{r}}_{i}}=0</math>
| |
|
| |
| ==Generalisierte Koordinaten==
| |
|
| |
| Problematischerweise liegen bei holonomen Zwangsbedingungen
| |
|
| |
|
| |
| <math>{{f}_{\lambda }}({{\vec{r}}_{1}}(t),{{\vec{r}}_{2}}(t),{{\vec{r}}_{3}}(t),...{{\vec{r}}_{N}}(t),t)=0\quad \quad \lambda =1,...\nu \quad f\ddot{u}r\ alle\ t</math>
| |
|
| |
|
| |
| gekoppelte Koordinaten vor ( die Koordinaten sind in den Zwangsbedingungen gekoppelt).
| |
|
| |
| Somit können die Punktkoordinaten
| |
|
| |
|
| |
| <math>\left\{ {{{\vec{r}}}_{1}}(t),{{{\vec{r}}}_{2}}(t),{{{\vec{r}}}_{3}}(t),...{{{\vec{r}}}_{N}}(t) \right\}</math>
| |
| nicht unabhängig voneinander variiert werden.
| |
|
| |
| <u>'''Ziel:'''</u>
| |
|
| |
| Suche einen Satz von f unabhängigen generalisierten Koordinaten. Diese sind optimal angepasst, wenn so viele unabhängige Koordinaten wie Freiheitsgrade existieren:
| |
|
| |
|
| |
| <math>\left\{ {{q}_{1}}(t),{{q}_{2}}(t),...{{q}_{f}}(t) \right\}\quad f=1,2,...3N-\nu </math>
| |
|
| |
|
| |
| Anschließend können Bewegungsgleichungen für die
| |
| <math>\left\{ {{q}_{1}}(t),{{q}_{2}}(t),...{{q}_{f}}(t) \right\}\quad f=1,2,...3N-\nu </math>
| |
| aus einfachen Extremalprinzipien ermittelt werden.
| |
|
| |
| Wesentlich: Die
| |
| <math>\left\{ {{q}_{1}}(t),{{q}_{2}}(t),...{{q}_{f}}(t) \right\}\quad f=1,2,...3N-\nu </math>
| |
| sind FREI variierbar ! Wegen
| |
|
| |
|
| |
| <math>{{\vec{r}}_{i}}={{\vec{r}}_{i}}\left( {{q}_{1}}(t),{{q}_{2}}(t),...{{q}_{f}}(t) \right)\quad f=1,2,...3N-\nu </math>
| |
| sind die Zwangsbedingungen identisch erfüllt.
| |
|
| |
| <u>'''Beispiel: Der Massenpunkt auf der bewegten Ebene:'''</u>
| |
|
| |
|
| |
| <math>\vec{a}\cdot (\vec{r}-{{\vec{r}}_{o}}(t))=0</math>
| |
|
| |
|
| |
|
| |
| Betrachten wir ein mitbewegtes Koordinatensystem
| |
| <math>\bar{e}{{\acute{\ }}_{1}},\bar{e}{{\acute{\ }}_{2}}</math>
| |
|
| |
|
| |
| Für den Radiusvektor existiert dann eine Verallgemeinerung:
| |
|
| |
|
| |
| <math>\bar{r}={{\bar{r}}_{o}}(t)+{{q}_{1}}\bar{e}{{\acute{\ }}_{1}}+{{q}_{2}}\bar{e}{{\acute{\ }}_{2}}</math>
| |
|
| |
|
| |
| Somit existiert eine injektive Abbildung der Koordinaten und wir können als generalisierte Koordinaten bestimmen:
| |
|
| |
|
| |
| <math>\left\{ {{q}_{1}},{{q}_{2}} \right\}</math>
| |
| , f=2
| |
|
| |
| <u>'''Beispiel: Massepunkt auf Kreis mit Radius R:'''</u>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \bar{r}=R(\cos \phi {{{\bar{e}}}_{1}}+\sin \phi {{{\bar{e}}}_{2}}) \\
| |
| & q=\phi \\
| |
| & f=1 \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| '''Virtuelle Verrückungen'''
| |
|
| |
| müssen nun auch in den generalisierten Koordinaten ausgedrückt werden, also:
| |
|
| |
|
| |
| <math>\delta {{\bar{r}}_{i}}</math>
| |
| wird ausgedrückt durch
| |
| <math>\delta {{q}_{1}},...,\delta qf</math>
| |
|
| |
|
| |
|
| |
|
| |
| Betrachten wir eine reale Verrückung ( in der Zeit), so gilt:
| |
|
| |
|
| |
| <math>{{\vec{v}}_{i}}=\frac{d}{dt}{{\bar{r}}_{i}}=\sum\limits_{j=1}^{f}{\left( \frac{\partial {{{\bar{r}}}_{i}}}{\partial {{q}_{j}}}{{{\dot{q}}}_{j}} \right)}+\frac{\partial }{\partial t}{{\bar{r}}_{i}}</math>
| |
|
| |
|
| |
| Daraus ergibt sich jedoch die Gleichung:
| |
|
| |
|
| |
| <math>\frac{\partial }{\partial {{{\dot{q}}}_{j}}}{{\vec{v}}_{i}}=\frac{\partial }{\partial {{{\dot{q}}}_{j}}}{{\left[ \sum\limits_{j=1}^{f}{\left( \frac{\partial {{{\bar{r}}}_{i}}}{\partial {{q}_{j}}}{{{\dot{q}}}_{j}} \right)}+\frac{\partial }{\partial t}{{{\bar{r}}}_{i}} \right]}_{{}}}=\frac{\partial }{\partial {{q}_{j}}}{{\bar{r}}_{i}}({{q}_{1}},...,{{q}_{f}},t)</math>
| |
|
| |
|
| |
| Mit diesen Gleichung kann die Virtuelle Arbeit der eingeprägten Kräfte gewonnen werden:
| |
|
| |
|
| |
| <math>\sum\limits_{i}^{{}}{{{{\vec{X}}}_{i}}\delta {{{\vec{r}}}_{i}}}=\sum\limits_{j}{\left\{ \sum\limits_{i}^{{}}{{{{\vec{X}}}_{i}}\frac{\partial {{{\bar{r}}}_{i}}}{\partial {{q}_{j}}}} \right\}\delta q_{j}^{{}}}=\sum\limits_{j=1}^{f}{{{Q}_{j}}\delta }q_{j}^{{}}</math>
| |
|
| |
|
| |
| Somit kann man als Ausdruck für die verallgemeinerte Kraft angeben:
| |
|
| |
|
| |
| <math>{{Q}_{j}}=\sum\limits_{i}^{{}}{{{{\vec{X}}}_{i}}\frac{\partial {{{\bar{r}}}_{i}}}{\partial {{q}_{j}}}}</math>
| |
|
| |
|
| |
| Sind die eingeprägten Kräfte konservativ:
| |
|
| |
|
| |
| <math>{{\vec{X}}_{i}}=-{{\nabla }_{\vec{r}i}}V({{\bar{r}}_{1}},{{\bar{r}}_{2}},...,{{\bar{r}}_{N}})</math>
| |
|
| |
|
| |
| So folgt:
| |
|
| |
|
| |
| <math>-\frac{\partial V}{\partial {{q}_{j}}}=-\sum\limits_{i}^{{}}{{{\nabla }_{\vec{r}i}}V({{{\bar{r}}}_{1}},{{{\bar{r}}}_{2}},...,{{{\bar{r}}}_{N}})\frac{\partial {{{\bar{r}}}_{i}}}{\partial {{q}_{j}}}}=\sum\limits_{i}^{{}}{{{{\vec{X}}}_{i}}\frac{\partial {{{\bar{r}}}_{i}}}{\partial {{q}_{j}}}}={{Q}_{j}}</math>
| |
|
| |
|
| |
| Somit besitzen auch die verallgemeinerten Kräfte ein Potenzial, natürlich das physikalisch gleiche wie die eingeprägten Kräfte !
| |
|
| |
| == Lagrangegleichungen 2. Art==
| |
|
| |
| Betrachten wir wieder das d Álembertsche Prinzip:
| |
|
| |
|
| |
| <math>\sum\limits_{i}{{{m}_{i}}{{{\ddot{\bar{r}}}}_{i}}\delta {{{\bar{r}}}_{i}}}=\sum\limits_{i}^{{}}{{{{\vec{X}}}_{i}}\delta {{{\bar{r}}}_{i}}=\sum\limits_{j}{{{Q}_{j}}}}\delta {{q}_{j}}</math>
| |
|
| |
|
| |
| Linke Seite:
| |
|
| |
|
| |
| <math>\sum\limits_{i}{{{m}_{i}}{{{\ddot{\bar{r}}}}_{i}}\delta {{{\bar{r}}}_{i}}}=\sum\limits_{j}{{}}\left( \sum\limits_{i}{{{m}_{i}}{{{\ddot{\bar{r}}}}_{i}}\frac{\partial }{\partial {{q}_{j}}}{{{\bar{r}}}_{i}}} \right)\delta {{q}_{j}}=\sum\limits_{j}{{}}\sum\limits_{i}^{{}}{\left\{ \frac{d}{dt}\left( {{m}_{i}}{{{\dot{\vec{r}}}}_{i}}\frac{\partial }{\partial {{q}_{j}}}{{{\bar{r}}}_{i}} \right)-{{m}_{i}}{{{\dot{\vec{r}}}}_{i}}\frac{d}{dt}\left( \frac{\partial }{\partial {{q}_{j}}}{{{\bar{r}}}_{i}} \right) \right\}\delta {{q}_{j}}_{{}}}</math>
| |
|
| |
|
| |
| Mit
| |
|
| |
|
| |
| <math>\frac{\partial }{\partial {{{\dot{q}}}_{j}}}{{\vec{v}}_{i}}=\frac{\partial }{\partial {{{\dot{q}}}_{j}}}{{\left[ \sum\limits_{j=1}^{f}{\left( \frac{\partial {{{\bar{r}}}_{i}}}{\partial {{q}_{j}}}{{{\dot{q}}}_{j}} \right)}+\frac{\partial }{\partial t}{{{\bar{r}}}_{i}} \right]}_{{}}}=\frac{\partial }{\partial {{q}_{j}}}{{\bar{r}}_{i}}({{q}_{1}},...,{{q}_{f}},t)</math>
| |
|
| |
|
| |
| und
| |
|
| |
|
| |
| <math>{{\dot{\bar{r}}}_{i}}=\sum\limits_{j=1}^{f}{\left( \frac{\partial {{{\bar{r}}}_{i}}}{\partial {{q}_{j}}} \right)}{{\dot{q}}_{j}}+\frac{\partial }{\partial t}{{\bar{r}}_{i}}=\sum\limits_{j=1}^{f}{\left( \frac{\partial }{\partial {{{\dot{q}}}_{j}}}{{{\vec{v}}}_{i}} \right)}{{\dot{q}}_{j}}+\frac{\partial }{\partial t}{{\bar{r}}_{i}}\Rightarrow \frac{d}{dt}\left( \frac{\partial {{{\bar{r}}}_{i}}}{\partial {{q}_{j}}} \right)=\frac{\partial }{\partial {{q}_{j}}}{{\vec{v}}_{i}}</math>
| |
|
| |
|
| |
| Beweis für die letzte Deduktion:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \frac{d}{dt}\left( \frac{\partial {{{\bar{r}}}_{i}}}{\partial {{q}_{j}}} \right)=\sum\limits_{k=1}^{{}}{\left( \frac{{{\partial }^{2}}{{{\bar{r}}}_{i}}}{\partial {{q}_{k}}\partial {{q}_{j}}} \right)}{{{\dot{q}}}_{k}}+\frac{{{\partial }^{2}}}{\partial {{q}_{j}}\partial t}{{{\bar{r}}}_{i}} \\
| |
| & \frac{\partial }{\partial {{q}_{j}}}{{{\vec{v}}}_{i}}=\frac{\partial }{\partial {{q}_{j}}}\left\{ \sum\limits_{k=1}^{{}}{\left( \frac{\partial {{{\bar{r}}}_{i}}}{\partial {{q}_{k}}} \right)}{{{\dot{q}}}_{k}}+\frac{\partial }{\partial t}{{{\bar{r}}}_{i}} \right\}=\sum\limits_{k=1}^{{}}{\left( \frac{{{\partial }^{2}}{{{\bar{r}}}_{i}}}{\partial {{q}_{k}}\partial {{q}_{j}}} \right)}{{{\dot{q}}}_{k}}+\frac{{{\partial }^{2}}}{\partial {{q}_{j}}\partial t}{{{\bar{r}}}_{i}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Somit ergibt sich für die linke Seite
| |
|
| |
|
| |
| <math>\sum\limits_{i}{{{m}_{i}}{{{\ddot{\bar{r}}}}_{i}}\delta {{{\bar{r}}}_{i}}}=\sum\limits_{i,j}{\left\{ \frac{d}{dt}\left( {{m}_{i}}{{{\vec{v}}}_{i}}\frac{\partial }{\partial {{{\dot{q}}}_{j}}}{{{\vec{v}}}_{i}} \right)-{{m}_{i}}{{{\vec{v}}}_{i}}\left( \frac{\partial }{\partial {{q}_{j}}}{{{\vec{v}}}_{i}} \right) \right\}\delta {{q}_{j}}}</math>
| |
|
| |
|
| |
|
| |
| Ziel ist es, diese Seite durch die gesamte KINETISCHE ENERGIE auszudrücken:
| |
|
| |
|
| |
| <math>T=\sum\limits_{i}{\frac{1}{2}{{m}_{i}}{{{\vec{v}}}_{i}}^{2}}</math>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{m}_{i}}{{{\vec{v}}}_{i}}\frac{\partial }{\partial {{{\dot{q}}}_{j}}}{{{\vec{v}}}_{i}}=\frac{\partial }{\partial {{{\dot{q}}}_{j}}}\left( \frac{1}{2}{{m}_{i}}{{{\vec{v}}}_{i}}^{2} \right) \\
| |
| & {{m}_{i}}{{{\vec{v}}}_{i}}\frac{\partial }{\partial {{q}_{j}}}{{{\vec{v}}}_{i}}=\frac{\partial }{\partial {{q}_{j}}}\left( \frac{1}{2}{{m}_{i}}{{{\vec{v}}}_{i}}^{2} \right) \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Somit folgt:
| |
|
| |
|
| |
| <math>\left\{ \frac{d}{dt}\left( {{m}_{i}}{{{\vec{v}}}_{i}}\frac{\partial }{\partial {{{\dot{q}}}_{j}}}{{{\vec{v}}}_{i}} \right)-{{m}_{i}}{{{\vec{v}}}_{i}}\left( \frac{\partial }{\partial {{q}_{j}}}{{{\vec{v}}}_{i}} \right) \right\}=\left\{ \frac{d}{dt}\left( \frac{\partial }{\partial {{{\dot{q}}}_{j}}}T \right)-\left( \frac{\partial }{\partial {{q}_{j}}}T \right) \right\}</math>
| |
|
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \sum\limits_{i}{{{m}_{i}}{{{\ddot{\bar{r}}}}_{i}}\delta {{{\bar{r}}}_{i}}}=\sum\limits_{i}^{{}}{{{{\vec{X}}}_{i}}\delta {{{\bar{r}}}_{i}}=\sum\limits_{j}{{{Q}_{j}}}}\delta {{q}_{j}} \\
| |
| & \Rightarrow \sum\limits_{j}{\left\{ \frac{d}{dt}\left( \frac{\partial }{\partial {{{\dot{q}}}_{j}}}T \right){{-}_{{}}}\left( \frac{\partial }{\partial {{q}_{j}}}T \right)-{{Q}_{j}} \right\}\delta {{q}_{j}}=0} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Der T-abhängige Ausdruck ist jedoch in qj völlig frei variierbar. Somit ist keine lineare Abhängigkeit der Variationen über verschiedene j gegeben.
| |
|
| |
| Jedes qj ist für sich frei variierbar, so dass der Ausdruck auf der linken Seite für sich Null wird:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \frac{d}{dt}\left( \frac{\partial }{\partial {{{\dot{q}}}_{j}}}T \right){{-}_{{}}}\left( \frac{\partial }{\partial {{q}_{j}}}T \right)-{{Q}_{j}}=0 \\
| |
| & \Rightarrow \frac{d}{dt}\left( \frac{\partial }{\partial {{{\dot{q}}}_{k}}}T \right){{-}_{{}}}\left( \frac{\partial }{\partial {{q}_{k}}}T \right)={{Q}_{k\quad \quad k=1,....,f}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| '''Lagrange- Gleichungen 2. Art:'''
| |
|
| |
| Die Lagrangegleichungen der zweiten Art können aus dem d ´Alembertschen Prinzip nur für HOLONOME Zwangsbedingungen gewonnen werden ( im Gegensatz zur Lagrangegleichung erster Art).
| |
|
| |
| Dies liegt daran, dass nur für HOLONOME Zwangsbedingungen generalisierte Koordinaten definiert werden können:
| |
|
| |
| <u>'''Spezialfall konservative Kräfte:'''</u>
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & -\frac{\partial V}{\partial {{q}_{j}}}={{Q}_{j}} \\
| |
| & V({{q}_{1}},...,{{q}_{f}},t)=V({{{\vec{r}}}_{1}}({{q}_{1,}}...,{{q}_{f}},t),...,{{{\vec{r}}}_{N}}({{q}_{1,}}...,{{q}_{f}},t)) \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Dies bedingt jedoch:
| |
|
| |
|
| |
| <math>\frac{\partial V({{q}_{1}},...,{{q}_{f}},t)}{\partial {{{\dot{q}}}_{k}}}=0</math>
| |
|
| |
|
| |
| Wir können uns die Lagrangefunktion derart definieren, dass:
| |
|
| |
|
| |
| <math>L({{q}_{1}},...,{{q}_{f}},{{\dot{q}}_{1}},...,{{\dot{q}}_{f}},t)=L({{q}_{k}},{{\dot{q}}_{k}},t)=T-V</math>
| |
|
| |
|
| |
| Es folgt:
| |
|
| |
|
| |
| <math>\frac{d}{dt}\left( \frac{\partial L}{\partial {{{\dot{q}}}_{k}}} \right)-\frac{\partial L}{\partial {{q}_{k}}}=0</math>
| |
|
| |
|
| |
| Die sagenumwobene Lagrangegleichung 2. Art für konservative Kräfte !
| |
|
| |
| Anmerkung:
| |
|
| |
| * die genannte Lagrangegleichung L ist nicht eindeutig festgelegt
| |
| * L=T-V ist nur EINE mögliche Form
| |
| *
| |
| <math>\begin{align}
| |
| & T({{q}_{k}},{{{\dot{q}}}_{k}},t)=\frac{1}{2}\sum\limits_{i=1}^{N}{{{m}_{i}}{{\left( \sum\limits_{k=1}^{f}{\frac{\partial {{{\vec{r}}}_{i}}}{\partial {{q}_{k}}}{{{\dot{q}}}_{k}}+\frac{\partial {{{\vec{r}}}_{i}}}{\partial t}} \right)}^{2}}} \\
| |
| & T({{q}_{k}},{{{\dot{q}}}_{k}},t)=a+\sum\limits_{k=1}^{f}{{{b}_{k}}{{{\dot{q}}}_{k}}}+\sum\limits_{k,l=1}^{f}{{{c}_{kl}}{{{\dot{q}}}_{k}}{{{\dot{q}}}_{l}}} \\
| |
| \end{align}</math>
| |
|
| |
| * Dabei ist die kinetische Energie nur für skleronome Zwangsbedingungen eine HOMOGENE Bilinearform in
| |
| <math>{{\dot{q}}_{k\quad }}(a={{b}_{k}}=0)</math>
| |
|
| |
|
| |
| Anwendungsschema für Lagrangegleichungen zweiter Art:
| |
|
| |
| Die Atwoodsche Fallmaschine
| |
|
| |
| Generalisierte Koordinate: q
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & T({{q}_{{}}},{{{\dot{q}}}_{{}}},t)=\frac{1}{2}({{m}_{{{1}_{{}}}}}+{{m}_{2}}){{{\dot{q}}}^{2}} \\
| |
| & V(q,\dot{q},t)={{m}_{1}}gq+{{m}_{2}}(l-q)g \\
| |
| & \frac{\partial L}{\partial q}={{m}_{1}}g-{{m}_{2}}g \\
| |
| & \frac{\partial L}{\partial \dot{q}}=({{m}_{1}}+{{m}_{2}})\dot{q} \\
| |
| & ({{m}_{1}}+{{m}_{2}})\ddot{q}+{{m}_{1}}g-{{m}_{2}}g=0 \\
| |
| & \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| <u>'''Beispiel 2:'''</u>
| |
|
| |
| Eine Masse m rotiert mit Winkelgeschwindigkeit w an einem Faden der Länge Ro, welcher mit Geschwindigkeit c durch ein Loch gezogen wird (rheonome Zwangsbedingung).
| |
|
| |
| Generalisierte Koordinate q ist der Winkel
| |
| <math>\phi </math>
| |
| :
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & T({{q}_{{}}},{{{\dot{q}}}_{{}}},t)=\frac{1}{2}m{{c}^{2}}+\frac{1}{2}m{{{\dot{q}}}^{2}}{{({{R}_{o}}^{{}}-ct)}^{2}} \\
| |
| & V(q,\dot{q},t)=0 \\
| |
| & L=\frac{1}{2}m{{c}^{2}}+\frac{1}{2}m{{{\dot{q}}}^{2}}{{({{R}_{o}}^{{}}-ct)}^{2}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Dahin kommt man im Übrigen aus:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & T=\frac{1}{2}m({{{\dot{x}}}^{2}}+{{{\dot{y}}}^{2}}) \\
| |
| & x=({{R}_{o}}-ct)\cos \phi \\
| |
| & \dot{x}=-c\cos \phi -({{R}_{o}}-ct)\dot{\phi }\sin \phi =-c\cos q-({{R}_{o}}-ct)\dot{q}\sin q \\
| |
| & y=({{R}_{o}}-ct)\sin \phi \\
| |
| \end{align}</math>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \frac{\partial L}{\partial \dot{q}}=m\dot{q}{{({{R}_{o}}^{{}}-ct)}^{2}} \\
| |
| & \frac{d}{dt}\frac{\partial L}{\partial \dot{q}}=m\ddot{q}{{({{R}_{o}}^{{}}-ct)}^{2}}-2cm\dot{q}({{R}_{o}}^{{}}-ct) \\
| |
| & \Rightarrow \ddot{q}({{R}_{o}}-ct)=2c{{{\dot{q}}}^{{}}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Somit haben wir eine Bewegungsgleichung für die Winkelgeschwindigkeit gefunden:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \frac{{\dot{\omega }}}{\omega }=\frac{2c}{{{R}_{o}}-ct} \\
| |
| & \int{\frac{d\omega }{\omega }=2c\int{\frac{dt}{{{R}_{o}}-ct}}} \\
| |
| & \ln \omega =-2\ln ({{R}_{o}}-ct)+const \\
| |
| & \ln \omega =\ln \frac{con\tilde{s}}{{{({{R}_{o}}-ct)}^{2}}} \\
| |
| & \omega =\frac{con\tilde{s}}{{{({{R}_{o}}-ct)}^{2}}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Bestimmung der Konstanten aus den Anfangsbedingungen liefert:
| |
|
| |
| Drehimpuls:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \vec{L}=m\vec{v}\times \vec{r} \\
| |
| & {{{\vec{L}}}_{o}}=m{{\omega }_{o}}^{{}}{{R}_{o}}^{2}\quad {{v}_{o}}={{\omega }_{o}}{{R}_{o}}\quad {{r}_{o}}={{R}_{o}} \\
| |
| & andererseits: \\
| |
| & {{\omega }_{o}}=\frac{con\tilde{s}}{{{({{R}_{o}})}^{2}}} \\
| |
| & \Rightarrow \omega =\frac{con\tilde{s}}{{{({{R}_{o}}-ct)}^{2}}}\Rightarrow con\tilde{s}=\frac{{{{\vec{L}}}_{o}}}{m} \\
| |
| & \omega =\frac{{{{\vec{L}}}_{o}}}{m{{({{R}_{o}}-ct)}^{2}}}=\dot{q} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Durch Integration gewinnt man:
| |
|
| |
|
| |
| <math>q={{q}_{o}}+\frac{{{{\vec{L}}}_{o}}}{cm({{R}_{o}}-ct)}</math>
| |
|
| |
|
| |
| Das heißt, wie zu erwarten war, die Masse dreht sich immer schneller, je kürzer der Faden wird ( Drehimpulserhaltung !)
| |
|
| |
| ==Normalschwingungen==
| |
|
| |
| Anwendung: Kleine Schwingungen eines Systems von Massepunkten
| |
| <math>{{m}_{i}}</math>
| |
|
| |
|
| |
| Die Zwangsbedingungen seien holonom und skleronom.
| |
|
| |
| Außerdem sei das Potenzial beliebig
| |
|
| |
|
| |
| <math>V({{\bar{r}}_{1}},{{\bar{r}}_{2}},...,{{\bar{r}}_{N}})</math>
| |
|
| |
|
| |
| es existiere lediglich eine stabile Ruhelage.
| |
|
| |
| Dazu wähle man generalisierte Koordinaten ( f Stück) mit der Ruhelage 0
| |
|
| |
| Man kann an dieses Problem herangehen, indem die potenzielle Energie um die Ruhelage entwickelt wird:
| |
|
| |
|
| |
| <math>V({{q}_{1}},...,{{q}_{f}})=V(0,....,0)+\sum\limits_{j}{{{\left( \frac{\partial V}{\partial {{q}_{j}}} \right)}_{0}}{{q}_{j}}+\frac{1}{2}\sum\limits_{j,k}{{{\left( \frac{{{\partial }^{2}}V}{\partial {{q}_{j}}\partial {{q}_{k}}} \right)}_{0}}{{q}_{j}}{{q}_{k}}+...}}</math>
| |
|
| |
|
| |
| Der erste Term kann gleich Null gesetzt werden ( Skalenverschiebung bei Potenzialen). Dies entspricht einer Skalenverschiebung der Energie.
| |
|
| |
| Im Zweiten Term tauchen jedoch die verallgemeinerten Kräfte ( von außen) auf. Wenn diese nicht existieren, so ist dieser Term ebenfalls Null:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & V(0,....,0)=0 \\
| |
| & \sum\limits_{j}{{{\left( \frac{\partial V}{\partial {{q}_{j}}} \right)}_{0}}{{q}_{j}}}=0\quad \left( \frac{\partial V}{\partial {{q}_{j}}} \right)=-{{Q}_{j}}=0 \\
| |
| & \frac{1}{2}\sum\limits_{j,k}{{{\left( \frac{{{\partial }^{2}}V}{\partial {{q}_{j}}\partial {{q}_{k}}} \right)}_{0}}{{q}_{j}}{{q}_{k}}}=\frac{1}{2}\sum\limits_{j,k}{{{V}_{jk}}{{q}_{j}}{{q}_{k}}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Für kleine Schwingungen hinreichend genau erhalten wir also in niedrigster Näherung grundsätzlich harmonische Schwingungen in einem q²- Potenzial :
| |
|
| |
| Das Potenzial ergibt eine positiv definite quadratische Form ( positiv definit, da Ruhelage stabil !)
| |
|
| |
|
| |
| <math>V({{q}_{1}},...,{{q}_{f}})\approx \frac{1}{2}\sum\limits_{j,k}{{{V}_{jk}}{{q}_{j}}{{q}_{k}}\ge 0}\quad \quad {{V}_{jk}}={{V}_{kj}}</math>
| |
|
| |
|
| |
| Ansatz für die kinetische Energie:
| |
|
| |
|
| |
| <math>T=\frac{1}{2}\sum\limits_{i}{{{m}_{i}}{{{\vec{v}}}_{i}}^{2}}\ge 0</math>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{{\vec{v}}}_{i}}=\sum\limits_{j}{{}}\left( \frac{\partial {{{\vec{r}}}_{i}}}{\partial {{q}_{j}}} \right){{{\dot{q}}}_{j}} \\
| |
| & T=\frac{1}{2}\sum\limits_{i}{{{m}_{i}}}\left( \sum\limits_{j,k}{\left( \frac{\partial {{{\vec{r}}}_{i}}}{\partial {{q}_{j}}} \right)\left( \frac{\partial {{{\vec{r}}}_{i}}}{\partial {{q}_{j}}} \right)}{{{\dot{q}}}_{j}}{{{\dot{q}}}_{k}} \right)\ge 0 \\
| |
| & T=\frac{1}{2}\sum\limits_{j,k}{{{T}_{jk}}}{{{\dot{q}}}_{j}}{{{\dot{q}}}_{k}} \\
| |
| & {{T}_{jk}}={{T}_{kj}}\approx \sum\limits_{i}{{{m}_{i}}{{\left( \frac{\partial {{{\vec{r}}}_{i}}}{\partial {{q}_{j}}} \right)}_{0}}{{\left( \frac{\partial {{{\vec{r}}}_{i}}}{\partial {{q}_{j}}} \right)}_{0}}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Die Auswertung der Ableitungen des Radiusvektor an der Ruhelage (0) gilt dann als niedrigste ( quadratische) Näherung für kleine Schwingungen.
| |
|
| |
| Auch die kinetische Energie ist in unserem Fall nun eine positiv definite quadratische Form.
| |
|
| |
| Die Lagrangegleichung 2. Art ist somit vollständig bestimmt:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & L=T-V=\frac{1}{2}\left( \sum\limits_{j,k}{{{T}_{jk}}}{{{\dot{q}}}_{j}}{{{\dot{q}}}_{k}}-\sum\limits_{j,k}{{{V}_{jk}}}{{q}_{j}}{{q}_{k}} \right) \\
| |
| & \frac{\partial L}{\partial {{{\dot{q}}}_{l}}}=\frac{1}{2}\sum\limits_{j,k}{{{T}_{jk}}}\frac{\partial }{\partial {{{\dot{q}}}_{l}}}\left( {{{\dot{q}}}_{j}}{{{\dot{q}}}_{k}} \right)=\frac{1}{2}\sum\limits_{j,k}{{{T}_{jk}}}\left( {{\delta }_{jl}}{{{\dot{q}}}_{k}}+{{\delta }_{kl}}{{{\dot{q}}}_{j}} \right)=\frac{1}{2}\sum\limits_{j,k}{{{T}_{lk}}}{{{\dot{q}}}_{k}}+{{T}_{lj}}{{{\dot{q}}}_{j}}=\sum\limits_{k}{{{T}_{lk}}}{{{\dot{q}}}_{k}}\quad mit\ {{T}_{jl}}={{T}_{lj}} \\
| |
| & \frac{d}{dt}\left( \frac{\partial L}{\partial {{{\dot{q}}}_{l}}} \right)=\sum\limits_{k}{{{T}_{lk}}}{{{\ddot{q}}}_{k}} \\
| |
| & \frac{\partial L}{\partial {{q}_{l}}}=-\sum\limits_{k}{{{V}_{lk}}{{q}_{k}}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| <u>'''Einschub: Transformation auf Kugelkoordinaten:'''</u>
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \left( r,\vartheta ,\phi \right)=\left( {{q}_{1}},{{q}_{2}},{{q}_{3}} \right) \\
| |
| & x=r\cos \phi \sin \vartheta \\
| |
| & y=r\sin \phi \sin \vartheta \\
| |
| & z=r\cos \vartheta \\
| |
| \end{align}</math>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{{\vec{v}}}_{{}}}=\sum\limits_{j}{{}}\left( \frac{\partial {{{\vec{r}}}_{{}}}}{\partial {{q}_{j}}} \right){{{\dot{q}}}_{j}} \\
| |
| & \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| In Komponenten ergibt sich somit:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{v}_{x}}=\frac{dx}{dt}=\frac{\partial x}{\partial r}\dot{r}+\frac{\partial x}{\partial \vartheta }\dot{\vartheta }+\frac{\partial x}{\partial \phi }\dot{\phi }=\sin \vartheta \cos \phi \dot{r}+r\cos \vartheta \cos \phi \dot{\vartheta }-r\sin \vartheta \sin \phi \dot{\phi } \\
| |
| & {{v}_{y}}=\frac{dy}{dt}=\frac{\partial y}{\partial r}\dot{r}+\frac{\partial y}{\partial \vartheta }\dot{\vartheta }+\frac{\partial y}{\partial \phi }\dot{\phi }=\sin \vartheta \sin \phi \dot{r}+r\cos \vartheta \sin \phi \dot{\vartheta }+r\sin \vartheta \cos \phi \dot{\phi } \\
| |
| & {{v}_{z}}=\frac{dz}{dt}=\frac{\partial z}{\partial r}\dot{r}+\frac{\partial z}{\partial \vartheta }\dot{\vartheta }+\frac{\partial z}{\partial \phi }\dot{\phi }=\cos \vartheta \dot{r}-r\sin \vartheta \dot{\vartheta } \\
| |
| & \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Es läßt sich eine Funktionalmatrix zusammenstellen:
| |
|
| |
|
| |
| <math>\left( \begin{matrix}
| |
| \frac{\partial x}{\partial r} & \frac{\partial x}{\partial \vartheta } & \frac{\partial x}{\partial \phi } \\
| |
| \frac{\partial y}{\partial r} & \frac{\partial y}{\partial \vartheta } & \frac{\partial y}{\partial \phi } \\
| |
| \frac{\partial z}{\partial r} & \frac{\partial z}{\partial \vartheta } & \frac{\partial z}{\partial \phi } \\
| |
| \end{matrix} \right)=\left( \begin{matrix}
| |
| \sin \vartheta \cos \phi & r\cos \vartheta \cos \phi & -r\sin \vartheta \sin \phi \\
| |
| in\vartheta \sin \phi & r\cos \vartheta \sin \phi & r\sin \vartheta \cos \phi \\
| |
| \cos \vartheta & -r\sin \vartheta & 0 \\
| |
| \end{matrix} \right)</math>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & T=\frac{1}{2}\sum\limits_{j,k}{{{T}_{jk}}}{{{\dot{q}}}_{j}}{{{\dot{q}}}_{k}} \\
| |
| & {{T}_{jk}}={{T}_{kj}}\approx \sum\limits_{i}{{{m}_{i}}{{\left( \frac{\partial {{{\vec{r}}}_{i}}}{\partial {{q}_{j}}} \right)}_{0}}{{\left( \frac{\partial {{{\vec{r}}}_{i}}}{\partial {{q}_{j}}} \right)}_{0}}} \\
| |
| & {{T}_{jk}}=m\left[ \left( \frac{\partial x}{\partial {{q}_{j}}} \right)\left( \frac{\partial x}{\partial {{q}_{k}}} \right)+\left( \frac{\partial y}{\partial {{q}_{j}}} \right)\left( \frac{\partial y}{\partial {{q}_{k}}} \right)+\left( \frac{\partial z}{\partial {{q}_{j}}} \right)\left( \frac{\partial z}{\partial {{q}_{k}}} \right) \right] \\
| |
| \end{align}</math>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{T}_{11}}=m\left( {{\sin }^{2}}\vartheta {{\cos }^{2}}\phi +{{\sin }^{2}}\vartheta {{\sin }^{2}}\phi +{{\cos }^{2}}\vartheta \right)=m \\
| |
| & {{T}_{22}}=m{{r}^{2}}\left( {{\cos }^{2}}\vartheta {{\cos }^{2}}\phi +{{\cos }^{2}}\vartheta {{\sin }^{2}}\phi +{{\sin }^{2}}\vartheta \right)=m{{r}^{2}} \\
| |
| & {{T}_{33}}=m{{r}^{2}}({{\sin }^{2}}\vartheta {{\sin }^{2}}\phi +{{\sin }^{2}}\vartheta {{\cos }^{2}}\phi )=m{{r}^{2}}{{\sin }^{2}}\vartheta \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Diese Wert hängen dabei von den gewählten Koordinaten, also den qj ab.
| |
|
| |
| Aus diesem Grund ( um dies zu erreichen) wurden ja gerade die qj so eingeführt.
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{T}_{12}}={{T}_{21}}=mr\left( \sin \vartheta \cos \phi \cos \vartheta \cos \phi +\sin \vartheta \sin \phi \cos \vartheta \sin \phi -\sin \vartheta \cos \vartheta \right)=0 \\
| |
| & {{T}_{13}}={{T}_{31}}=0 \\
| |
| & {{T}_{23}}={{T}_{32}}=0 \\
| |
| \end{align}</math>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{T}_{jk}}=\left( \begin{matrix}
| |
| m & 0 & 0 \\
| |
| 0 & m{{r}^{2}} & 0 \\
| |
| 0 & 0 & m{{r}^{2}}{{\sin }^{2}}\vartheta \\
| |
| \end{matrix} \right) \\
| |
| & T=\frac{1}{2}m\left( {{{\dot{r}}}^{2}}+{{r}^{2}}{{{\dot{\vartheta }}}^{2}}+{{r}^{2}}{{\sin }^{2}}\vartheta {{{\dot{\phi }}}^{2}} \right) \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Zurück:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & L=T-V=\frac{1}{2}\left( \sum\limits_{j,k}{{{T}_{jk}}}{{{\dot{q}}}_{j}}{{{\dot{q}}}_{k}}-\sum\limits_{j,k}{{{V}_{jk}}}{{q}_{j}}{{q}_{k}} \right) \\
| |
| & \frac{\partial L}{\partial {{{\dot{q}}}_{l}}}=\frac{1}{2}\sum\limits_{j,k}{{{T}_{jk}}}\frac{\partial }{\partial {{{\dot{q}}}_{l}}}\left( {{{\dot{q}}}_{j}}{{{\dot{q}}}_{k}} \right)=\frac{1}{2}\sum\limits_{j,k}{{{T}_{jk}}}\left( {{\delta }_{jl}}{{{\dot{q}}}_{k}}+{{\delta }_{kl}}{{{\dot{q}}}_{j}} \right)=\frac{1}{2}\sum\limits_{j,k}{{{T}_{lk}}}{{{\dot{q}}}_{k}}+{{T}_{lj}}{{{\dot{q}}}_{j}}=\sum\limits_{k}{{{T}_{lk}}}{{{\dot{q}}}_{k}}\quad mit\ {{T}_{jl}}={{T}_{lj}} \\
| |
| & \frac{d}{dt}\left( \frac{\partial L}{\partial {{{\dot{q}}}_{l}}} \right)=\sum\limits_{k}{{{T}_{lk}}}{{{\ddot{q}}}_{k}} \\
| |
| & \frac{\partial L}{\partial {{q}_{l}}}=-\sum\limits_{k}{{{V}_{lk}}{{q}_{k}}} \\
| |
| & \Rightarrow \sum\limits_{k}{{{T}_{lk}}}{{{\ddot{q}}}_{k}}+{{V}_{lk}}{{q}_{k}}=0\quad \quad l=1,...,f \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Somit haben wir ein System von f linearen Differenzialgleichungen gegeben.
| |
|
| |
| Bekanntlich eignet sich als Ansatz für die Lösung:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{q}_{k}}(t)={{A}_{k}}{{e}^{iwt}}\quad {{A}_{k}}\in C \\
| |
| & \sum\limits_{k}{({{V}_{lk}}-{{\omega }^{2}}{{T}_{lk}}){{A}_{k}}=0} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Dies ist eine Eigenwertgleichung für w²
| |
|
| |
| Bei gegebenen w² liegt ein lineares Gleichungssystem für die Ak vor:
| |
|
| |
| Eine nichttriviale Lösung existiert aber genau dann, wenn
| |
|
| |
|
| |
| <math>\det \left( {{V}_{lk}}-{{\omega }^{2}}{{T}_{lk}} \right)=0</math>
| |
|
| |
|
| |
| Dies ist die charakteristische Gleichung für w², die sogenannte Säkulargleichung, ein Polynom f-ten Grades.
| |
|
| |
|
| |
| <math>{{V}_{lk}},{{T}_{lk}}positiv\ definit\Rightarrow {{\omega }^{2}}>0</math>
| |
| für alle Nullstellen des charakteristischen Polynoms.
| |
|
| |
| <u>Beweis:</u>
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \sum\limits_{k}{({{V}_{lk}}-{{\omega }^{2}}{{T}_{lk}}){{A}_{k}}=0}\left| \cdot \sum\limits_{l}{{{A}_{l}}^{*}} \right. \\
| |
| & \sum\limits_{l,k}{{{V}_{lk}}{{A}_{l}}^{*}{{A}_{k}}-}{{\omega }^{2}}\sum\limits_{l,k}{{{T}_{lk}}{{A}_{l}}^{*}{{A}_{k}}}=0 \\
| |
| & {{\omega }^{2}}=\frac{\sum\limits_{l,k}{{{V}_{lk}}{{A}_{l}}^{*}{{A}_{k}}}}{\sum\limits_{l,k}{{{T}_{lk}}{{A}_{l}}^{*}{{A}_{k}}}} \\
| |
| & \sum\limits_{l,k}{{{V}_{lk}}{{A}_{l}}^{*}{{A}_{k}}}=\frac{1}{2}\sum\limits_{l,k}{{{V}_{lk}}{{A}_{l}}^{*}{{A}_{k}}}+\frac{1}{2}\sum\limits_{l,k}{{{V}_{kl}}{{A}_{k}}^{*}{{A}_{l}}=}\frac{1}{2}\sum\limits_{l,k}{{{V}_{lk}}\left( {{A}_{l}}^{*}{{A}_{k}}+{{A}_{k}}^{*}{{A}_{l}} \right)=}\frac{1}{2}\sum\limits_{l,k}{{{V}_{lk}}2\cdot \operatorname{Re}\left( {{A}_{l}}^{*}{{A}_{k}} \right)} \\
| |
| \end{align}</math>
| |
| Also handelt es sich hierbei um eine reelle quadratische Form. Nun sind Vlk und Tlk positiv definite Matrizen.
| |
|
| |
| Zähler und Nenner sind aber reelle quadratische Formen.
| |
|
| |
| Was zur Folge hat, dass w²>0
| |
|
| |
| Die Lösungen des Gleichungssystems
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{q}_{k}}(t)={{A}_{k}}{{e}^{iwt}}\quad {{A}_{k}}\in C \\
| |
| & \sum\limits_{k}{({{V}_{lk}}-{{\omega }^{2}}{{T}_{lk}}){{A}_{k}}=0} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| sind die Eigenfrequenzen
| |
| <math>{{\omega }^{2}}_{a}\quad a=1,...,f</math>
| |
|
| |
|
| |
| und die Eigenvektoren
| |
| <math>{{A}_{k}}^{(a)}\quad a=1,...,f</math>
| |
|
| |
|
| |
| Wobei die Eigenvektoren nur bis auf einen Normierungsfaktor bestimmt sind und reell gewählt werden können.
| |
|
| |
| Die allgemeine Lösung für die verallgemeinerten Kooridnaten lautet:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{q}_{k}}(t)=\operatorname{Re}\left\{ \sum\limits_{a=1}^{f}{{{C}_{a}}}{{A}_{k}}^{(a)}{{e}^{i{{w}_{a}}t}} \right\} \\
| |
| & \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Die
| |
| <math>{{C}_{a}}</math>
| |
| werden durch die Anfangsbedingungen
| |
| <math>{{q}_{k}}(0),{{\dot{q}}_{k}}(0)</math>
| |
| bestimmt
| |
|
| |
| <u>'''Normalkoordinaten'''</u>
| |
|
| |
| Ziel:
| |
|
| |
| Transformiere auf neue generalisierte Koordinaten, so dass die Bewegungsgleichungen für die Koordinaten entkoppeln.
| |
|
| |
| Seien diese neuen Koordinaten
| |
| <math>{{Q}_{j}}</math>
| |
| so soll gelten:
| |
|
| |
|
| |
| <math>{{\ddot{Q}}_{j}}+{{\omega }_{j}}^{2}{{Q}_{j}}=0\quad j=1,...,f</math>
| |
|
| |
|
| |
| Dies wird bekanntlich erreicht durch eine Hauptachsentransformation der symmetrischen Matrizen Vlk und Tlk
| |
|
| |
| Die Transformation ist das Diagonalisierungsverfahren. Dazu werden reell gewählte Eigenvektoren
| |
|
| |
|
| |
| <math>{{A}_{k}}^{(a)}</math>
| |
| eingesetzt. In diesen müssen sich dann die generalisierten Koordinaten mit den Normalkoordinaten als Entwicklungskoeffizienten darstellen lassen:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{q}_{k}}(t)=\sum\limits_{a=1}^{f}{{}}{{A}_{k}}^{(a)}{{Q}_{a}} \\
| |
| & \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Die diagonalisierte Matrix kann die Koordinatentransformation als Abbildung vollständig darstellen:
| |
|
| |
|
| |
|
| |
| <math>\vec{q}=A\vec{Q}\quad mit\ \vec{q},\vec{Q}\in {{R}^{f}}</math>
| |
|
| |
|
| |
|
| |
|
| |
| Bleibt zu zeigen, dass Vlk und Tlk durch das gleiche System von Eigenvektoren diagonalisiert werden:
| |
|
| |
| Es gelten die Eigenwertgleichungen:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \sum\limits_{k}{({{V}_{lk}}-{{\omega }_{a}}^{2}{{T}_{lk}}){{A}_{k}}^{a}=0}\left| \cdot \sum\limits_{l}{{{A}_{l}}^{b}} \right. \\
| |
| & \sum\limits_{l}{({{V}_{kl}}-{{\omega }_{b}}^{2}{{T}_{kl}}){{A}_{l}}^{b}=0\left| \cdot \sum\limits_{k}{{{A}_{k}}^{b}} \right.} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \sum\limits_{k,l}{{{A}_{l}}^{b}({{V}_{lk}}-{{V}_{kl}}){{A}_{k}}^{a}-{{A}_{l}}^{b}({{\omega }_{a}}^{2}{{T}_{lk}}-{{\omega }_{b}}^{2}{{T}_{kl}}){{A}_{k}}^{a}}=0 \\
| |
| & {{V}_{lk}}={{V}_{kl}} \\
| |
| & \sum\limits_{k,l}{({{\omega }_{a}}^{2}-{{\omega }_{b}}^{2}){{A}_{l}}^{b}{{T}_{kl}}{{A}_{k}}^{a}}=0 \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Die Annahme lautet nun noch:
| |
|
| |
|
| |
| <math>{{\omega }_{a}}^{2}-{{\omega }_{b}}^{2}\ne 0</math>
| |
| Die Eigenwerte sind nicht entartet, natürlich für verschiedene a/b
| |
|
| |
| Somit folgt jedoch
| |
|
| |
|
| |
| <math>\sum\limits_{k,l}{{}}{{A}_{l}}^{b}{{T}_{kl}}{{A}_{k}}^{a}={{\delta }_{ab}}</math>
| |
|
| |
|
| |
| Im wesentlichen ist dieser Ausruck ( die transformierte kinetische Energie)Null für verschiedene a und b. Bei geeigneter Normierung kann er für a=b gleich 1 gesetzt werden.
| |
|
| |
| Die Trafo ist also eine verallgemeinerte orthogonale Trafo.
| |
|
| |
| Es folgt wegen
| |
|
| |
|
| |
| <math>\sum\limits_{k}{({{V}_{lk}}-{{\omega }_{a}}^{2}{{T}_{lk}}){{A}_{k}}^{a}=0}\left| \cdot \sum\limits_{l}{{{A}_{l}}^{b}} \right.</math>
| |
|
| |
|
| |
| dass
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \sum\limits_{k,l}{({{A}_{l}}^{b}{{V}_{lk}}-{{\omega }_{a}}^{2}{{A}_{l}}^{b}{{T}_{lk}}){{A}_{k}}^{a}=0} \\
| |
| & \sum\limits_{k,l}{({{A}_{l}}^{b}{{V}_{lk}}{{A}_{k}}^{a})=\sum\limits_{k,l}{{{\omega }_{a}}^{2}{{A}_{l}}^{b}{{T}_{lk}}{{A}_{k}}^{a}}}={{\omega }_{a}}^{2}{{\delta }_{ab}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Also werden Tlk und Vlk durch das gleiche System von Eigenvektoren diagonalisiert.
| |
|
| |
| Lagrangefunktion:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & L=T-V=\frac{1}{2}\left( \sum\limits_{j,k}{{{T}_{jk}}}{{{\dot{q}}}_{j}}{{{\dot{q}}}_{k}}-\sum\limits_{j,k}{{{V}_{jk}}}{{q}_{j}}{{q}_{k}} \right) \\
| |
| & L=\frac{1}{2}\left( \sum\limits_{a,b}{\left( \sum\limits_{j,k}{{{A}_{j}}^{b}{{T}_{jk}}{{A}_{k}}^{a}{{{\dot{Q}}}_{a}}{{{\dot{Q}}}_{b}}-\sum\limits_{j,k}{{{A}_{j}}^{b}{{V}_{jk}}{{A}_{k}}^{a}{{Q}_{a}}{{Q}_{b}}}} \right)} \right) \\
| |
| & \sum\limits_{j,k}{{{A}_{j}}^{b}{{T}_{jk}}{{A}_{k}}^{a}={{\delta }_{ab}}} \\
| |
| & \sum\limits_{j,k}{{{A}_{j}}^{b}{{V}_{jk}}{{A}_{k}}^{a}={{\omega }_{a}}^{2}{{\delta }_{ab}}} \\
| |
| & L=\frac{1}{2}\left( \sum\limits_{a}{\left( {{{\dot{Q}}}_{a}}^{2}-{{\omega }_{a}}^{2}{{Q}_{a}}^{2} \right)} \right) \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| In der tat entkoppeln nun die Bewegungsgleichungen:
| |
|
| |
|
| |
| <math>{{\ddot{Q}}_{a}}+{{\omega }_{j}}^{2}{{Q}_{a}}=0\quad a=1,...,f</math>
| |
|
| |
|
| |
| <u>'''Beispiel: Pendel'''</u>
| |
|
| |
| Leicht kann man sich an einer Skizze klar machen:
| |
|
| |
|
| |
| <math>z=l(1-\cos \phi )</math>
| |
|
| |
|
| |
| Als verallgemeinerte Koordinate kann man die Bogenlänge wählen:
| |
|
| |
|
| |
| <math>q=s=\phi l</math>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & T=\frac{1}{2}m{{{\dot{q}}}^{2}} \\
| |
| & V=mgz=mgl(1-\cos \phi )\approx \frac{1}{2}mgl{{\phi }^{2}}=\frac{1}{2}\frac{g}{l}m{{q}^{2}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Die Entwicklung des Potenzials kann auführlich gezeigt werden.
| |
|
| |
| Nun seien zwei Pendel über eine Feder der Federkonstante k gekoppelt:
| |
|
| |
| <u>'''Zwei gekoppelte Pendel'''</u>
| |
|
| |
| Hier nehmen wir für beide Pendel generalisierte Koordinaten:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{q}_{1}}={{s}_{1}}={{\phi }_{1}}l \\
| |
| & {{q}_{2}}={{s}_{2}}={{\phi }_{1}}l \\
| |
| \end{align}</math>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & T=\frac{1}{2}m({{{\dot{q}}}_{1}}^{2}+{{{\dot{q}}}_{2}}^{2}) \\
| |
| & V=mg{{z}_{1}}+mg{{z}_{2}}+\frac{1}{2}k{{({{q}_{1}}-{{q}_{2}})}^{2}}=mgl(1-\cos \frac{{{q}_{1}}}{l})+\frac{1}{2}k{{({{q}_{1}}-{{q}_{2}})}^{2}}+mgl(1-\cos \frac{{{q}_{2}}}{l}) \\
| |
| & V\approx \frac{1}{2}mgl{{\phi }_{1}}^{2}+\frac{1}{2}mgl{{\phi }_{2}}^{2}+\frac{1}{2}k{{({{q}_{1}}-{{q}_{2}})}^{2}}=\frac{1}{2}\frac{g}{l}m{{q}_{1}}^{2}+\frac{1}{2}\frac{g}{l}m{{q}_{2}}^{2}+\frac{1}{2}k{{({{q}_{1}}-{{q}_{2}})}^{2}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Nun kann gefordert werden:
| |
|
| |
|
| |
| <math>V\approx \frac{1}{2}\frac{g}{l}m{{q}_{1}}^{2}+\frac{1}{2}\frac{g}{l}m{{q}_{2}}^{2}+\frac{1}{2}k{{({{q}_{1}}-{{q}_{2}})}^{2}}=\sum\limits_{j,k=1}^{2}{{{V}_{jk}}{{q}_{j}}{{q}_{k}}\quad Forderung!}</math>
| |
|
| |
|
| |
| Dies läßt sich direkt über die mehrdimensionale Taylorreihe zeigen, Mit Hilfe der Multiindizes:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{\left( \frac{{{\partial }^{2}}V}{\partial {{q}_{1}}^{2}} \right)}_{0}}={{\left( \frac{{{\partial }^{2}}V}{\partial {{q}_{2}}^{2}} \right)}_{0}}=m\frac{g}{l}+k \\
| |
| & \left( \frac{{{\partial }^{2}}V}{\partial {{q}_{1}}\partial {{q}_{2}}} \right)=mg\frac{\partial }{\partial {{q}_{1}}}(\sin \frac{{{q}_{2}}}{l})-k\frac{\partial }{\partial {{q}_{1}}}({{q}_{1}}-{{q}_{2}})=-k \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Somit läßt sich die kinetische Energie angeben:
| |
|
| |
| Somit lassen sich kinetische Energie und Potenzial als Matrizen angeben:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{T}_{lk}}=\left( \begin{matrix}
| |
| m & 0 \\
| |
| 0 & m \\
| |
| \end{matrix} \right) \\
| |
| & {{V}_{lk}}=\left( \begin{matrix}
| |
| m\frac{g}{l}+k & -k \\
| |
| -k & m\frac{g}{l}+k \\
| |
| \end{matrix} \right) \\
| |
| \end{align}</math>
| |
|
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & T=\frac{1}{2}m({{{\dot{q}}}_{1}}^{2}+{{{\dot{q}}}_{2}}^{2}) \\
| |
| & V\approx \frac{1}{2}mgl{{\phi }_{1}}^{2}+\frac{1}{2}mgl{{\phi }_{2}}^{2}+\frac{1}{2}k{{({{q}_{1}}-{{q}_{2}})}^{2}}=\frac{1}{2}\frac{g}{l}m{{q}_{1}}^{2}+\frac{1}{2}\frac{g}{l}m{{q}_{2}}^{2}+\frac{1}{2}k{{({{q}_{1}}-{{q}_{2}})}^{2}} \\
| |
| & L=T-V=\frac{1}{2}m({{{\dot{q}}}_{1}}^{2}+{{{\dot{q}}}_{2}}^{2})-\frac{1}{2}\frac{g}{l}m{{q}_{1}}^{2}-\frac{1}{2}\frac{g}{l}m{{q}_{2}}^{2}-\frac{1}{2}k{{({{q}_{1}}-{{q}_{2}})}^{2}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Die Bewegungsgleichungen ergeben sich als:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & m{{{\ddot{q}}}_{1}}+\frac{g}{l}m{{q}_{1}}+k({{q}_{1}}-{{q}_{2}})=0 \\
| |
| & m{{{\ddot{q}}}_{2}}+\frac{g}{l}m{{q}_{2}}-k({{q}_{1}}-{{q}_{2}})=0 \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Auch hier haben wir ein System gekoppelter Differenzialgleichungen.
| |
|
| |
| Als Loesungsansatz wählen wir:
| |
|
| |
|
| |
| <math>{{q}_{k}}={{A}_{k}}{{e}^{iwt}}</math>
| |
|
| |
|
| |
| Die resultierende Eigenwertgleichung lautet:
| |
|
| |
|
| |
| <math>\left( \begin{matrix}
| |
| m\frac{g}{l}+k-m{{\omega }^{2}} & -k \\
| |
| -k & m\frac{g}{l}+k-m{{\omega }^{2}} \\
| |
| \end{matrix} \right)\left( \begin{matrix}
| |
| {{A}_{1}} \\
| |
| {{A}_{2}} \\
| |
| \end{matrix} \right)=0</math>
| |
|
| |
|
| |
| Aus der charakteristischen Gleichung gewinnen wir das charakteristische Polynom
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & 0=\det ({{V}_{lk}}-{{\omega }^{2}}{{T}_{lk}})={{m}^{2}}\left| \begin{matrix}
| |
| \frac{g}{l}+\frac{k}{m}-{{\omega }^{2}} & -\frac{k}{m} \\
| |
| \frac{-k}{m} & \frac{g}{l}+\frac{k}{m}-{{\omega }^{2}} \\
| |
| \end{matrix} \right|=0 \\
| |
| & 0={{\omega }^{4}}-2\left( \frac{k}{m}+\frac{g}{l} \right){{\omega }^{2}}+\frac{{{g}^{2}}}{{{l}^{2}}}+2\frac{gk}{lm}={{\omega }^{4}}-2\left( \frac{k}{m}+\frac{g}{l} \right){{\omega }^{2}}+{{\left( \frac{g}{l}+\frac{k}{m} \right)}^{2}}-{{\left( \frac{k}{m} \right)}^{2}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
|
| |
| <math>{{\omega }_{1,2}}^{2}=\left( \frac{k}{m}+\frac{g}{l} \right)\pm {{\left( \frac{k}{m} \right)}^{{}}}=\left\{ \begin{matrix}
| |
| \frac{g}{l} \\
| |
| \frac{g}{l}+2\left( \frac{k}{m} \right) \\
| |
| \end{matrix} \right.</math>
| |
|
| |
|
| |
| Somit kennt das System die folgenden Eigenfrequenzen:
| |
|
| |
|
| |
| <math>{{\omega }_{1}}=\sqrt{\frac{g}{l}}:={{\omega }_{0}}</math>
| |
|
| |
|
| |
| ungestörte Pendelfrequenz
| |
|
| |
|
| |
| <math>{{\omega }_{2}}=\sqrt{\frac{g}{l}+2\frac{k}{m}}:={{\sqrt{{{\omega }_{0}}^{2}+2{{{\tilde{\omega }}}^{2}}}}_{{}}}</math>
| |
|
| |
|
| |
| Die zugehörigen Eigenvektoren lauten:
| |
|
| |
|
| |
| <math>\left( m\frac{g}{l}+k-m{{\omega }_{a}}^{2} \right){{A}_{1}}^{a}-k{{A}_{2}}^{a}=0</math>
| |
|
| |
|
| |
| Somit ergibt sich mit der ungestörten Pendelfrequenz w1:
| |
|
| |
|
| |
| <math>k{{A}_{1}}^{1}-k{{A}_{2}}^{1}=0\Rightarrow \left( \begin{matrix}
| |
| {{A}_{1}}^{1} \\
| |
| {{A}_{2}}^{1} \\
| |
| \end{matrix} \right)\propto \left( \begin{matrix}
| |
| 1 \\
| |
| 1 \\
| |
| \end{matrix} \right)</math>
| |
|
| |
|
| |
| Aus der Eigenfrequenz w2 ergibt sich:
| |
|
| |
|
| |
|
| |
| In Normalkoordinaten gilt für die Lösung des Ortes:
| |
|
| |
|
| |
| <math>{{q}_{k}}(t)={{A}_{k}}^{1}{{Q}_{1}}+{{A}_{k}}^{2}{{Q}_{2}}</math>
| |
|
| |
|
| |
| Bis auf einen konstanten Faktor.
| |
|
| |
| Die Umkehrung lautet:
| |
|
| |
|
| |
| <math>\left( \begin{matrix}
| |
| {{Q}_{1}} \\
| |
| {{Q}_{2}} \\
| |
| \end{matrix} \right)=\left( \begin{matrix}
| |
| {{A}_{1}}^{1} & {{A}_{2}}^{1} \\
| |
| {{A}_{1}}^{2} & {{A}_{2}}^{2} \\
| |
| \end{matrix} \right)\left( \begin{matrix}
| |
| {{q}_{1}} \\
| |
| {{q}_{2}} \\
| |
| \end{matrix} \right)</math>
| |
|
| |
|
| |
| Mit der zu oben transponierten Matrix ( Umkehrung)
| |
|
| |
| Die Eigenvektoren sind so zu normieren, dass:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & \sum\limits_{k,l}{{{A}_{l}}^{a}{{T}_{lk}}{{A}_{k}}^{a}=m\sum\limits_{k}{{{\left| {{A}_{k}}^{a} \right|}^{2}}}=1} \\
| |
| & \Rightarrow \left( \begin{matrix}
| |
| {{A}_{1}}^{1} \\
| |
| {{A}_{2}}^{1} \\
| |
| \end{matrix} \right)=\frac{1}{\sqrt{2m}}\left( \begin{matrix}
| |
| 1 \\
| |
| 1 \\
| |
| \end{matrix} \right) \\
| |
| & \left( \begin{matrix}
| |
| {{A}_{1}}^{2} \\
| |
| {{A}_{2}}^{2} \\
| |
| \end{matrix} \right)=\frac{1}{\sqrt{2m}}\left( \begin{matrix}
| |
| 1 \\
| |
| -1 \\
| |
| \end{matrix} \right) \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Es folgt für die Normalkoordinaten:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{Q}_{1}}=\frac{1}{\sqrt{2m}}({{q}_{1}}+{{q}_{2}})\quad SChwerpunktskoordinaten \\
| |
| & {{Q}_{2}}=\frac{1}{\sqrt{2m}}({{q}_{1}}-{{q}_{2}})\quad \operatorname{Re}lativkoordinaten \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| An Normalschwingungen existiert somit:
| |
|
| |
|
| |
| <math>\begin{align}
| |
| & {{\omega }_{1}}=\sqrt{\frac{g}{l}} \\
| |
| & {{\omega }_{2}}=\sqrt{\frac{g}{l}+2\frac{k}{m}} \\
| |
| \end{align}</math>
| |
|
| |
|
| |
| Dabei stellt ersteres die gleichphasige Schwerpunktsschwingung dar, letzteres repräsentiert die gegenphasige Relativschwingung.
| |
|
| |
| In Realität haben wir es mit einer beliebigen Überlagerung von Schwerpunktsschwingung und Relativschwingung zu tun.
| |
|
| |
| Dabei treten Überlagerungszustände als Schwebung auf.
| |
|
| |
| In Realität erhält man eine reine Schwerpunktschwingung, wenn die Anfangsbedingungen reine Lösung der Schwerpunktsskoordinaten sind.
| |
|
| |
| Eine Relativschwingung ergibt sich, wenn die Anfangsbedingung exakt eine Lösung der Relativkoordinaten repräsentieren.
| |
|
| |
| [[Kategorie:Mechanik]]
| |
