Spiralbewegung eines starren Körpers

Ich möchte zeigen, dass ein starrer Körper mit zwei Komponenten seines Winkelgeschwindigkeitsvektors und einer Komponente seines linearen Geschwindigkeitsvektors in Abwesenheit äußerer Kräfte und Drehmomente spiralförmige Trajektorien hat. Dies wird in verschiedenen Artikeln, die ich gelesen habe, normalerweise als selbstverständlich angesehen (siehe zum Beispiel diese Seite 4, Abschnitt 4 , oder diese Seite 206, erster Absatz ).

Ich betrachte einen bewegten starren Körper, dessen Schwerpunkt in einem Punkt liegt R 0 ( T ) in dem durch die Basisvektoren definierten Laborbezugssystem e X , e j , e z . Ich füge einen beweglichen Referenzrahmen bei { R 0 ( T ) ; e X ' ( T ) , e j ' ( T ) , e z ' ( T ) } zum Körper; dieser Rahmen ist zentriert R 0 ( T ) und seine Achsen drehen sich mit dem Rahmen des Körpers. Der rotierende Dreiklang ( e X ' , e j ' , e z ' ) kann durch drei Euler-Winkel charakterisiert werden θ 1 , θ 2 , θ 3 und die folgenden Transformationsregeln:

e ich ' = L ( θ 1 , θ 2 , θ 3 ) e ich

Ich übernehme die Tait-Bryan-Winkelkonvention mit an j X ' ' z ' intrinsische Definition: θ 1 ist das Gieren (gegen den Uhrzeigersinn um e j ), θ 2 die Tonhöhe (gegen den Uhrzeigersinn ca e X ' ) Und θ 3 die Bank (Rolle, gegen den Uhrzeigersinn herum e z ' ). Der Kürze halber definiere ich den Vektor θ = ( θ 1 , θ 2 , θ 3 ) der drei unabhängigen Euler-Winkel. Die Matrix der Transformation ist

L ( θ ) = R Z ( θ 3 ) R X ( θ 1 ) R Y ( θ 2 )
bei dem die R ich sind die Rotationsmatrizen in R 3 . Ich gehe davon aus, dass der Körper lineare Geschwindigkeit hat U ( T ) = U e z ' ( T ) entlang seiner anterior-posterioren Achse orientiert e z ' ( T ) , und Rotationsgeschwindigkeitsvektor ω ( T ) . In der gewählten Darstellung sind die Komponenten der Winkelgeschwindigkeit im Bezugssystem des Körpers

ω X = θ ˙ 2 cos θ 1 θ ˙ 3 Sünde θ 1 cos θ 2 ω j = θ ˙ 1 + θ ˙ 3 Sünde θ 2 ω z = θ ˙ 2 Sünde θ 1 + θ ˙ 3 cos θ 1 cos θ 2 .

Diese können gefunden werden, indem man das bemerkt

L T L ˙ = ( 0 ω z ω j ω z 0 ω X ω j ω X 0 )

Die Flugbahn des Körpers ist daher durch die Kurve seines Massenschwerpunkts gegeben R 0 während es sich durch den Raum bewegt. Im Laborrahmen lautet dies

R 0 ( T ) = 0 T U ( τ ) D τ = U 0 T e z ' ( τ ) D τ

Deutlich, e z ' dreht sich mit dem Körper und ist daher eine Funktion der Zeit durch die Winkelgeschwindigkeitskomponenten ω ich :

e z ' ( T ) = L ( T ) e z = ( cos θ 3 Sünde θ 1 + cos θ 1 Sünde θ 2 Sünde θ 3 ) e X + ( cos θ 1 cos θ 3 Sünde θ 2 + Sünde θ 1 Sünde θ 3 ) e j + ( cos θ 1 cos θ 2 ) e z

Wo    θ ich = 0 T θ ˙ ich D τ .

Das Problem ist kompliziert, denn um zu integrieren e z ' Ich muss die Beziehungen zwischen umkehren ω ich Und θ J . Beachten Sie, dass ich hier alle drei Komponenten der Winkelgeschwindigkeit behalte, um zu sehen, an welchem ​​Punkt und in welchem ​​Umfang es notwendig ist, nur zwei Rotationsfreiheitsgrade zu haben.

Allerdings auch dann, wenn ich davon ausgehe, dass die Winkelgeschwindigkeiten konstant sind θ ich = θ ich T , bekomme ich keine Ausdrücke, die enthalten ω ich ausdrücklich. Gibt es einen anderen Ansatz, der Gebrauch macht L T L ˙ = ω × ?

Ich bin mir nicht sicher, ob ich deine Frage verstehe. Sie wollen zeigen, dass 2 Komponenten der Winkelgeschwindigkeit und eine Komponente der Lineargeschwindigkeit immer eine Schraubenbahn haben. Gehen wir davon aus, dass kein Nettodrehmoment oder keine Nettokräfte aufgebracht werden? Oder sind wir ganz allgemein und arbeiten mit willkürlichen Kräften?
Danke für den Hinweis. Die äußere Nettokraft und das Drehmoment sind Null. Und ja, meine Frage ist genau diese. Warum in Gegenwart von zwei Komponenten ω X , ω z der Winkelgeschwindigkeit und eins U z der linearen Geschwindigkeit, sollte ich immer eine Schraubenbahn erhalten.
Es gibt also keine ω j Komponente noch U X Und U j Komponenten?
Das ist richtig. Aber im Prinzip, wenn es auch gab ω j , die Literatur scheint zu implizieren, dass ich sowieso eine Helix erhalten würde.
In Abwesenheit von Kräften ist die Bahn des Massenschwerpunkts eine gerade Linie, keine Spirale. Die Referenzen, die Sie geben, beziehen sich auf Mikroben, die im Wasser schwimmen, also haben sie äußere Kräfte? Zum freien Körper siehe Goldstein, Klassische Mechanik, Kapitel 5.6.
Was Sie sagen, gilt beispielsweise nicht für R e 1 Regime in Flüssigkeit.
Damit eine Kurve eine zylindrische Spirale nachzeichnet, benötigen Sie nur einen Winkel θ X und eine konstante Geschwindigkeit in z-Richtung. Um die Überschrift zu ändern, müssen Sie einen anderen Blickwinkel wählen θ z .
Sind die Geschwindigkeitskomponenten vorgeschrieben, weil sie sich ändern werden, wenn das Objekt umherwirbelt? Wenn sie immer konstant sind, ist eine Einschränkung (und damit eine äußere Kraft) erforderlich, um dies zu erzwingen. Wenn sie frei sind, dann sind sie nur Anfangsbedingungen.
Nirgendwo in Ihrer Arbeit werden Bewegungsgleichungen und Winkelbeschleunigungen erwähnt. Alles, womit Sie zu arbeiten versuchen, ist die Kinematik (zwischen Euler-Winkeln und Rotationsvektor), während Sie die Dynamik ignorieren. Ist das beabsichtigt? Wollen Sie eine Antwort, die die Kinematik eines spiralförmig bewegten Körpers beschreibt?
Beachten Sie, dass L L ˙ = ω × lässt sich leicht als zeitliche Ableitung der Richtungsvektoren auf einem rotierenden Koordinatensystem erklären.
L ˙ = ω × L
(Siehe Folie 14 in cs.unc.edu/~lin/COMP768-F07/LEC/rbd1.pdf )
@ ja72 Ja, ich interessiere mich nur für die lokale (zeitliche) Kinematik, vorausgesetzt, das schwimmende Objekt kann über interne Mechanismen eine konstante Geschwindigkeit und Winkelgeschwindigkeit aufrechterhalten.
Wenn Sie also zwei Rotationskomponenten festgelegt haben, stellt dies immer noch eine einzelne Rotation um eine schiefe Achse dar. Sie können ein neues Koordinatensystem definieren, das an der kombinierten Drehung ausgerichtet ist, und das Problem erheblich vereinfachen.
Können Sie mir bitte zeigen, wie das geht? Ich denke intuitiv, dass ich eine Drehung um eine Achse habe, die sich selbst dreht (was lokal eine spiralförmige Bewegung verursacht), anstatt eine Drehung um eine schräge Achse ...
Drehen Sie einfach das Koordinatensystem so, dass es bei zwei Komponenten eine Drehung von Null und bei einer Komponente ungleich Null hat. Zum Beispiel ( 0 ω j ω z ) mit einer Drehung um x von bräunen 1 ( ω j ω z ) der Vektor wird ( 0 0 ω j 2 + ω z 2 ) .

Antworten (4)

Das Problem ist kompliziert, denn um zu integrieren e z ' wir müssen die Beziehungen zwischen umkehren ω ich Und θ J . Nehmen wir an, dass die Winkelgeschwindigkeiten konstant sind und nehmen wir die Kleinwinkelnäherung:

θ ˙ ich = konst.         θ ich ( T ) = θ ˙ ich T
Die Kleinwinkelnäherung impliziert, dass die Winkelgeschwindigkeitskomponenten im Körperrahmen werden:
Ω X θ ˙ 2 Ω j θ ˙ 1 Ω z θ ˙ 3

Mit diesen Annahmen erhalten wir:

e z ' ( T ) = cos θ 3 Sünde θ 1 e X + Sünde θ 1 Sünde θ 3 e j + cos θ 1 e z = cos ( Ω z T ) Sünde ( Ω j T ) e X + Sünde ( Ω j T ) Sünde ( Ω z T ) e j + cos ( Ω j T ) e z
die zur Ableitung integriert werden können R 0 ( T ) :
R 0 ( T ) = U 0 T e z ' ( τ ) D τ = U Ω j Ω z [ Ω j Ω j cos ( Ω j T ) cos ( Ω z T ) Ω z Sünde ( Ω j T ) Sünde ( Ω z T ) ] e X + U Ω j Ω z [ Ω j cos ( Ω j T ) Sünde ( Ω z T ) + Ω z Sünde ( Ω j T ) cos ( Ω z T ) ] e j + U Ω j Sünde ( Ω j T ) e z
Dies ist im Allgemeinen eine komplexe und sehr faszinierende Kurve, aber wir werden unsere Studie hier auf den Fall beschränken, wo Ω j Ω z . Eine Serienerweiterung von R 0 um Ω j Erträge
R 0 ( T ) U Ω j Ω z [ 1 + cos ( Ω z T ) + Ω z T Sünde ( Ω z T ) ] e X + U Ω j Ω z [ Ω z T cos ( Ω z T ) + Sünde ( Ω z T ) ] e j + T e z = ( R cos ( Ω z T ) R Sünde ( Ω z T ) T ) kreisförmige Spirale + ( T   Ω z R Sünde ( Ω z T ) T   Ω z R cos ( Ω z T ) 0 ) Spiral- + R ( 1 0 0 ) Achsenverschiebung

Die Wachstumsrate der Spirale ist sehr klein, ρ ˙ = U Ω j / Ω z 2 U , aufgrund der Annäherung Ω j Ω z . Daher ist die Trajektorie eine quasi-kreisförmige Wendel mit Radius R .

Der Zweck dieses Beitrags ist es zu demonstrieren, wie man eine Inverse vom Körperrahmen zur Referenz berechnet, nachdem Drehungen im Körperrahmen unter Verwendung der Tait-Bryan-Winkel durchgeführt wurden https://upload.wikimedia.org/wikipedia/commons/5/53/ Taitbrianzyx.svg

Ein Teil des Materials ist identisch mit dem Beitrag von OP, aber der Schwerpunkt liegt auf der Reihenfolge der Rotationen und der Verfolgung der durch die Körperrotationen induzierten multiplen Körperrahmen.

Letztendlich erleichtert es die Berechnung der Rücktransformation vom Körper- zum Bezugssystem und der Winkelgeschwindigkeiten.

Auf der anderen Seite kann es sein, dass das Problem übertrieben ist.

Referenzrahmen: R = ( X , Y , Z ) .

Körperrahmen: R = ( X , j , z ) .

Euler-Winkel: ( ψ , θ , ϕ ) .

Erste Drehung: Kurswinkel ψ , Wo 0 ψ 2 π - ist die positive Rotation ψ über die X Achse zum Rahmen R ' = ( X ' , j ' , z ' ) . Der Körperrahmen R stimmt mit dem Bezugssystem überein R .

Zweite Drehung: Haltungswinkel θ , Wo π / 2 θ π / 2 - positive Drehung um j ' Achse zum Rahmen R = ( X , j , z ) .

Dritte Rotation: Querneigungswinkel ϕ , Wo 0 θ 2 π - positive Drehung ungefähr ungefähr X Achse zum Rahmen R = ( X , j , z )

[ ψ ] = [ C Ö S ( ψ ) S ich N ( ψ ) 0 S ich N ( ψ ) C Ö S ( ψ ) 0 0 0 1 ]

[ θ ] = [ C Ö S ( θ ) 0 S ich N ( θ ) 0 1 0 S ich N ( θ ) 0 C Ö S ( θ ) ]

[ ϕ ] = [ 1 0 0 0 C Ö S ( ϕ ) S ich N ( ϕ ) 0 S ich N ( ϕ ) C Ö S ( ϕ ) ]

[ R ' ] = [ ψ ] [ R ]
[ R ] = [ θ ] [ R ' ]
[ R ] = [ ϕ ] [ R ]

oder

[ R ] = [ ϕ ] [ θ ] [ ψ ] [ R ]

Wo [ ϕ ] , [ θ ] , Und [ ψ ] sind Matrizen, und R , R , R ' , R sind Vektoren.

Seit [ ϕ ] , [ θ ] , Und [ ψ ] sind in dem Sinne orthogonal D e T [ ψ ] = 1 , usw. dann [ ψ ] T = [ ψ ] 1 , etc. und die inverse Transformation ist

[ R ] = [ ψ ] T [ θ ] T [ ϕ ] T [ R ] .
Für ψ ˙ entlang z ' , ϕ ˙ entlang X Und X , θ ˙ entlang j ' Und j
ω = ψ ˙ + θ ˙ + ϕ ˙

[ ω ] = [ ϕ ] [ θ ] [ ψ ˙ ] + [ ϕ ] [ θ ˙ ] + [ ϕ ˙ ]

Wo

[ ϕ ˙ ] = [ ϕ ˙ , 0 , 0 ] T
[ θ ˙ ] = [ 0 , θ ˙ , 0 ] T
[ ψ ˙ ] = [ 0 , 0 , ψ ˙ ] T .
Die Winkelgeschwindigkeiten sind
ω X = ϕ ˙ ψ ˙ S ich N ( θ )
ω j = θ ˙ C Ö S ( ϕ ) + ψ ˙ C Ö S ( θ ) S ich N ( ϕ )
ω z = ψ ˙ C Ö S ( θ ) C Ö S ( ϕ ) θ ˙ S ich N ( ϕ ) .

Ich habe eine allgemeine Antwort, um dies zu beweisen. Der erste Schritt besteht darin, den Augenblick festzustellen, in dem sich das Objekt mit einer parallelen Translation (schraubenförmige Bewegung) um eine Rotationsachse bewegt. Da dann die Geschwindigkeitskomponenten spezifiziert und erzwungen werden (dafür werden Kräfte benötigt), bleiben die Bedingungen für eine spiralförmige Bewegung über die Zeit bestehen.

Nehmen Sie also ein rotierendes Objekt mit Rotationsgeschwindigkeit ω , und mit linearer Geschwindigkeit v an einem bestimmten Punkt.

Bild1

Zerlegen Sie als Nächstes den Geschwindigkeitsvektor entlang der Rotationsachse und senkrecht dazu

Bild2

und beachten Sie, dass der parallele Vektor durch einen Skalarwert beschrieben werden kann H so dass

v = H ω
Auch der senkrechte Vektor wird erklärt, indem die Rotationsachse außermittig an einen Ort verschoben wird R so dass
v = R × ω

Bild3

Die Geschwindigkeit entlang paralleler und senkrechter Vektoren zerlegen

v = R × ω + H ω
Verwenden Sie die Position der Drehung
R = ω × v ω 2

Ebenso die Skalartonhöhe H wird gefunden von

H = ω v ω 2

Bild4

Das Ergebnis ist, dass die Bewegung des starren Körpers immer in eine Rotation um eine versetzte Achse zerlegt werden kann, die mit einer parallelen Geschwindigkeit zu dieser Achse gekoppelt ist.

Die Eigenschaften der Spirale sind wie folgt:

Eigentum Wert Größe ω = ω Richtung z ^ = ω ω Tonhöhe H = ω v ω 2 Standort R = ω × v ω 2

und die Bewegung von der Helix ist:

v = ω ( R × z ^ + H z ^ )

Der letzte Teil besteht darin, zu beachten, dass die Richtung, Steigung und Position der Spirale unverändert bleiben muss, damit sie sich zeitlich nicht ändert. Sie können sich ausrechnen, welche Geschwindigkeits- und Rotationskomponenten angegeben werden müssen, damit dies der Fall ist. Hinweis: Sie müssen 5 von 6 Komponenten angeben (3 lineare und 3 winklige), da die Richtung 2 Größen hat (ihre Größe spielt keine Rolle), die Tonhöhe skalar ist (1 Größe) und die Position der Rotationsachse 2 Größen hat (da die Position entlang der Achse nicht zählt).

Danke, ich bin mir immer noch nicht sicher, ob dies eine spiralförmige Bewegung ist. Können Sie das erweitern? Vielleicht zeigen, warum die parametrischen 3D-Gleichungen einer solchen Trajektorie zwangsläufig die einer Helix sein müssen?
Es ist das, was Sie erhalten, wenn Sie eine Rotation um eine Achse in Kombination mit einer parallelen Translation haben. Das Ergebnis ist eine Spirale. In der Mechanik nennt man das Schraubenbewegung .
@ ja72: Je schneller sich das Objekt dreht, desto langsamer ist der vertikale Anstieg. Kannst du auch die Krümmung der Spirale berechnen? Wie sieht es mit der Torsion aus?

Frenet-Rahmen

Siehe [Kinematik der helikalen Bewegung von Organismen mit bis zu sechs Freiheitsgraden] (link.springer.com/article/10.1007%2FBF02460302), auf die vom OP verwiesen wird, das den Frenet-Rahmen verwendet.

Der Frenet-Rahmen oder Trieder wird verwendet, um das Eindrehen und Eindrehen von Kurven zu beschreiben R 3 .

Es gibt 6 Freiheitsgrade, nämlich 3 Freiheitsgrade, um einen Punkt zu fixieren R 3 , und 3 Freiheitsgrade, um die zu fixieren T Tangentenvektor oder ein Geschwindigkeitsvektor an den Punkt.

Aus diesen 6 Freiheitsgraden folgt alles andere, nämlich der Frenet-Rahmen oder Trieder.

Der Frenet-Rahmen besteht aus 3 Vektoren, nämlich T ' , N ' , B ' aus denen gerechnet wird T .

Der Vektor T ' misst die Krümmung der Kurve und den Vektor B ' misst die Torsion der Kurve.

Der Parameter S erzeugt die Trajektorie.

Lassen β : ICH R 3 sei die Abbildung der zylindrischen Helixkurve in R 3 mit S im offenen Intervall 0 < S < 1 In ICH

β ( S ) = ( A C Ö S S C , A S ich N S C , B S C )

Wo C = ( A 2 + B 2 ) 1 / 2 Und A > 0 .

Der Einheits-Tangentenvektor an die Kurve T Ist

T ( S ) = β ' ( S ) = ( A C S ich N S C , A C C Ö S S C , B C )

Wo β ' ( S ) ist eine Kurzform für D β ( S ) D S .

Dann

T ' ( S ) = β ( S ) = ( A C 2 C Ö S S C , A C 2 S ich N S C , 0 )

und die Krümmung ist

κ ( S ) = | | T ' ( S ) | | = A C 2 = A A 2 + B 2 > 0

Seit T ' = κ N , Wo N ist der Normalenvektor

N ( S ) = ( C Ö S S C , S ich N S C , 0 )

die immer auf die Achse des Zylinders zeigt, auf dem β ( S ) liegt unabhängig von Werten von A Und B .

Mit dem Kreuzprodukt B = T × N gibt

B ( S ) = ( B C S ich N S C , B C C Ö S S C , A C )

Um die Torsion zu berechnen, gilt per Definition: B ' = τ N

B ' ( S ) = ( B C 2 C Ö S S C , B C 2 S ich N S C , 0 )

was impliziert

τ = B C 2 = B A 2 + B 2

Zusammenfassend besteht der Frenet-Rahmen aus den 3 orthonormalen Vektoren

T ' = κ N
N ' = κ T + τ B
B ' = τ N ,

Die erste Gleichung ist die Definition der Krümmung und die letzte Gleichung ist die Definition der Torsion - beide Gleichungen wurden direkt aus berechnet T Vektor. Die zweite Gleichung wurde als orthonormale Entwicklung der ersten und letzten Gleichung geschrieben.