• ベストアンサー

運動方程式の導出で困っているので教えてくれませんか

中村 拓男(@tknakamuri)の回答

回答No.1

まず, X, Y, Z 座標は中心星を原点に据えた普通のデカルト座標 x, y, z 座標は中心星を原点に、木星を x軸上に固定した 回転座標です。 木星の交点面は XY 平面と一致していて、角速度Ωで回っているようですね。 (1) ω=(0, 0, Ω) は X, Y, Z座標で表した木星の角速度です。 (2) i', j', k' は木星の方向(単位ベクトル)をデカルト座標で表しています。   これは回転座標の基準となる3軸で、i'は木星の動径方向、j' は木星の進行方向、   k' は公転面に垂直な方向です。 ベクトル a が原点に対してωで回転すると da/dt = ω x a ですから、di'/dt = ω x i' です。 (3) 外積を計算しているだけです。計算に座標系のひねりがありますが、 線形代数か物理数学の公式をあたってみてください。   求めているのはデカルト座標(X, Y, Z)での値です。 (4) X, Y, Z 座標における回転座標の3軸の方向の変化率です。

関連するQ&A

  • 運動座標系における運動方程式の導出

    http://hooktail.sub.jp/mechanics/acCoordinates/ こちらのホームページに運動座標系における運動方程式の導出について説明がされていますが、疑問があります。 (4)の式について疑問があるのですが、この結果を導くには(dω/dt)×rというベクトル積が0である必要がありますよね? 私にはなぜこのベクトル積が0になるのかがよくわかりません。 またこの説明にも書かれているようにdω/dt=δω/δtとなりますよね。 だとすれば(δω/dt)×rというベクトル積はなぜ0にならないんですか? ここの意味がよくわかりません。 もしわかる方がいらっしゃいましたら教えていただけないでしょうか。 よろしくお願いします。

  • 回転&運動座標系のベクトル

    運動座標系に対する運動方程式の計算の仕方がわからないので、お知恵をお借りしたく思います。ベクトルは大文字で示すことにさせてください。(ωもベクトル) 運動座標系 S' と慣性系 S があるとして、S'の原点の S系でのベクトルを R0 とし、S'系での質点の座標を R' とすると R=R0+R'でS系での座標があらわせ、 S系での速度は V= dR/dt = dR0/dt + dR'/dt となります。(ここまでは問題なし) S'系の単位ベクトル(I', J', K')を使うと、 dR'/dt = DR'/dt + ω x R' --- (1) となります。ただし、R'=x'I'+y'J'+z'K' としたとき DR'/dt= dx'/dt I' + dy'/dt J' + dz'/dt K' とします。 (ここまではたぶん理解できていると思います。もし、間違えていたらご指摘ください。) ここからがよくわからないのですが、 加速度ですが、 dV/dt = d^2R0/dt^2 + d^2R'/dt^2 より、後半の d^2 R'/d t^2 = D/dt(dR'/dt) + ω x dR'/dt ----(2) = D/dt(DR'/dt + ω x R') + ω x (DR'/dt + ω x R') ---- (3) = D^2 R'/dt^2 + 2ω x DR'/dt + ω x (ω x R') + Dω/dt x R' ----(4) (dω/dt = Dω/dt + ω x ω および ω x ω = 0 より) 以上より、 F = m dV / dt = m d^2R0/dt^2 + m d^2R'/dt^2 とすると m D^2 R'/dt^2 = m d^2 R'/d t^2 - 2 m ω x DR'/dt - m ω x (ω x R') - m Dω/dt x R' ----(5) = F - m d^2R0/dt^2 - 2 m (ω x DR'/dt) - m ω x (ω x R') - m dω/dt x R' とのことですが、(2) - (4) のあたりがよくわかりません。おしえていただけると助かります。 まず (2) は (1)にd/dt をかけたのだと思いますが、 d/dt (DR'/dt) -> D/dt(dR'/dt) になって、 d/dt (ω x R')-> ω x dR'/dt ということでよろしいのでしょうか? また、(4) は (3) を一つずつ、 D/dt(DR'/dt + ω x R') -> D^2R'/dt + Dω/dt x R' + ω x DR'/dt ω x (DR'/dt + ω x R') -> ω x DR'/dt + ω x (ω x R') としたのだと思いますが、 Dω/dt x R' は ω x R'とはならないでしょうか。 (d/dt (ω x R')-> ω x dR'/dtとは対照的に) よろしくお願いいたします。

  • 困っているので教えてください。

    この式でコリオリ力と遠心力の項はどうやって変形して出しているのかわかりません。 途中式などやり方や考え方を教えてください。 Ωなどはどこから出てきているのですか? よろしくお願いします。 拡大画像↓ http://kie.nu/uF8 http://kie.nu/uFc

  • コンデンサとコイルのみの回路の回路方程式

    図1について,q/C =L(dI/dt),-I=dq/dt この2式より q/C =-L(d^2q/dt^2) 図2について,q/C =L(dI/dt),I=dq/dt この2式より q/C =L(d^2q/dt^2) 質問1 私の立てた式はこれで合っているでしょうか。 質問2 ある本を読んだ時に,-/C =-L(dI/dt),-I=dq/dt この2式より q/C =-L(d^2q/dt^2)という記述を見かけました。  両辺に負の符号が付いているので,結局私の立てた式と同じだと思うのですが,疑問をもったのはなぜわざわざ両辺に負の符号を付けているのでしょうか。この意味を理解したいのですが,よく分かりません。 質問3 図2と同じ図を使って,q=CV=C×(-L(dI/dt))と立てた式を見かけました。(見間違いかもしれません。)両辺をCで割れば,q/C=-L(dI/dt),I=dq/dtなので,q/C =-L(d^2q/dt^2)となり,私の立てた式と正負が異なります。なぜ,これについて説明していただけると助かります。

  • 強制振動の運動方程式の解き方

    f0/m cosωt=-kx-2кd²x/dt²-ω²0x という運動方程式をAcosθ+Bsinθ=√(A²+B²)cos(θ-β) ただしtanβ=B/Aを使って解きたいのですが解き方がわかりません。 β=arctan{2кω/(ω²0-ω)}とするようです。 ご教授お願い致します。 f0やω0の0は下付き文字です。

  • この式の意味や導出など教えてください。

    座標系?か何かの変形みたいなことをしていたとき d/dt i'(t)=ω×i' という式がでてきたのですが、何を表しているのかよくわかりません。 導出とか式の意味・概念など解説していただけないでしょうか。 図などを用いてなるべくわかりやすく教えていただけると助かります。

  • 物理の運動方程式の時間積分について

    物理をしていて、運動方程式をとく際、画像の意味がわかりませんでした。4C2と書いてあるところの横です。問題は「m・d^2x/dt^2=qEsinωtの運動方程式をとけ」でした。この画像の4C2の2行目の流れがわかりません。よろしくお願いします

  • 剛体のオイラーの運動方程式において、観測者の視点は何処にあるんですか?

    剛体のオイラーの運動方程式において、観測者の視点は何処にあるんですか? 慣性主軸をとったときを考えます。 この慣性主軸のことを(つまり慣性系において動く物体に固定された座標系のことを)、 名前がわからないので以下ではとりあえず物体系と呼ぶことにします。 この時、 L = Iω となりますが、 L(ベクトル),I(テンソル),ω(ベクトル)は全て「物体系における」量ですよね? これを時間微分して、つまり、 dL/dt = d'L/dt + ω×L によって得られるオイラーの方程式について、 1、そもそも上でいう時間微分とは、どの座標系における微分なのでしょうか。  自分は、物体系における微分であり、d'L/dtは「物体系に乗っている回転座標系」における微分だと解釈してるのですが、あっているでしょうか。 2、オイラーの方程式から、物体系に対するωを求めることが出来ますが(多分。質問1における自分の解釈が間違っていたらこの質問は破綻。)、では慣性系におけるωはどうやって求めるのでしょうか。 言い換えると、慣性系に対する物体系の動きはどうやって求めるのでしょうか。 重心の軌跡を求める、ということではなく、例えば重心の移動が無いとき、 慣性系に対して物体系はどのように回転するのか、ということです。 宜しくお願いします。

  • 連立微分方程式が解けません!

    a(t)dt=iωexp(iΔt)b(t)/2 b(t)dt=iωexp(-iΔt)a(t)/2 上記2式において、a(0)=1,b(0)=0 が成り立つとき、a(t)及びb(t)を求めよ。 畑違いの分野なので、さっぱり分かりません…… でも、これが解けないとまずいことになるんです。 どうか、どなたか教えていただけないでしょうか? 心の底よりお願い申し上げます。

  • 物理の微分方程式についてです

    振り子の振動の問題を解いているのですが 微分方程式 (d^2θ/dt^2)=-gθ/L この式から角振動数ωを求めるにはどのようにすればよいでしょうか