Rotierende Referenzrahmen

Ich versuche, die Gleichungen zu verstehen, die die Geschwindigkeit in einem rotierenden Referenzrahmen regeln ...

v ich = ( D R D T ) R + Ω × R .

Ich möchte eine einfache Simulation einer Rakete erstellen, die mit einer konstanten Trägheitsgeschwindigkeit von der Erde abhebt, sagen wir: v ich = [ 1 , 0 , 0 ] T .

Ich vermute einige Ω Wert, um die Rotation der Erde um z darzustellen, sagen wir Ω = [ 0 , 0 , 1 ] .

Und dann nach der wahrgenommenen Geschwindigkeit im Rotationsrahmen auflösen:

( D R D T ) R = v ich Ω × R .

Was ich erwarte, nachdem ich die Geschwindigkeit in die Position integriert habe, wäre eine nach außen rotierende Spirale, die die relative Position der "Rakete" zu einem Beobachter im rotierenden Erdrahmen zeigt. Was ich von einer einfachen Simulink-Simulation sehe, ist ganz anders.

Mein Sim:

Geben Sie hier die Bildbeschreibung ein

Die Ausgabe:

Geben Sie hier die Bildbeschreibung ein

Gedanken?

Könntest du posten, was du siehst, ein Bild oder so?
Wäre Computational Science ein besseres Zuhause für diese Frage?
Es ist eine physikalische Frage im Herzen. Die Simulation soll mir nur helfen zu verstehen, was zum Teufel passiert.
@Qmechanic Wenn es um den Algorithmus ginge, wäre es so. Aber das ist mehr Physik als Numerik, würde ich sagen.
Verwenden Sie eine Rotationsgeschwindigkeit, die entweder 10 oder 0,1 beträgt, sodass sie unterschiedliche Größenordnungen haben. Ich werde in ein paar Minuten ausführlicher antworten.
Es ist nicht erstaunlich klar, wie die Notation auf dieser Ableitung lautet.
Zeitliche Ableitung der Position im Rotationsrahmen
Mit Ihrer ersten Formel die Trägheitsgeschwindigkeit v ich berechnet - richtig. Aber es drückt sich im rotierenden Rahmen aus!

Antworten (3)

Koordinatenvektor eines Punktes im statischen Rahmen: R S

Koordinatenvektor desselben Punktes im rotierenden Rahmen: R R

(Reine Rotation, beide Frames haben den gleichen Ursprung.)

Koordinatentransformation (Rotationsmatrix): R

Die Matrix ist orthogonal, dh R T R = R R T = 1 (die Einheitsmatrix)

Wichtige Eigenschaft: 0 = D D T 1 = D D T ( R T R ) = R ˙ T R + R T R ˙

Das bedeutet die Matrix Ω := R T R ˙ ist antisymmetrisch Ω = Ω T (mit nur drei relevanten Komponenten Ω 1 := Ω 32 , Ω 2 := Ω 13 , Ω 3 := Ω 21 ) und die Produkte Ω v kann mit dem Vektor ausgedrückt werden Ω = ( Ω 1 , Ω 2 , Ω 3 ) als Ω × v .

Koordinatenvektor im Drehrahmen:

R S = R R R

Geschwindigkeit, zeitliche Ableitung im statischen Rahmen:

v S S := D R S D T = R ˙ R R + R R ˙ R

Anwenden R T zu dieser Gleichung:

R T v S S = R T R ˙ R R + R ˙ R

Sie sehen, wie Sie die Geschwindigkeit transformieren v S S in den Drehrahmen (dasselbe wo auch R R Leben). Der richtige Name in unserer Nomenklatur für die im statischen Koordinatensystem berechnete und in das rotierende übertragene Zeitableitung wäre v S R = R T v S S .

Damit erhalten Sie Ihre Formel

v S R = Ω × R R + R ˙ R .

Lassen Sie mich Ihnen zuerst die Mathematik zeigen: Ihre Position im rotierenden Rahmen sein X R wo im statischen Rahmen es sein X S . Jetzt ist die Geschwindigkeit:

v S = D D T X S .

Das sollte in Ordnung sein. Um nun zum rotierenden Rahmen zu gelangen, müssen wir die Notation ein wenig erweitern. Wir müssen uns die Einheitsvektoren ansehen und den Ortsvektor in Komponenten zerlegen:

X S = ich = 1 3 X S ich e ich ,
Wo e ich sind die Einheitsvektoren des Bezugssystems. Im statischen Rahmen sind dies:
e 1 = ( 1 0 0 ) , e 2 = ( 0 1 0 ) , e 3 = ( 0 0 1 )

Sie hängen nicht von der Zeit ab, also müssen Sie das nicht berücksichtigen, wenn Sie eine Zeitableitung machen.

Der rotierende Rahmen hat Einheitsvektoren, die zeitabhängig sind, Sie haben die folgende Produktregel:

v R = D D T ich = 1 3 X R ich e ich = ich = 1 3 ( X ˙ R ich e ich + X R ich e ˙ ich )

Der letzte Begriff ist, wo Ihr Ω × kommt ins Spiel. Siehe Ableitung der Zentrifugal- und Corioliskraft für die verbleibende Ableitung.

Wenn Sie dieses Problem auf andere Weise betrachten möchten, können Sie zur Notation komplexer Zahlen wechseln (da dies ein planares Problem ist).

Wenn im nicht rotierenden System die Position der Rakete ist z = γ ( T ) , dann in einem mit konstanter Winkelgeschwindigkeit rotierenden System ω , du hast z = e ich ω T γ ( T ) . Wenn γ ( T ) = v T , Dann

z = e ich ω T v T = k θ e ich θ , mit  k = v ω ,
das ist, wie Sie sagten, eine Spirale.