• ベストアンサー

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

d↑i'/dt=↑ω×↑i'=‥‥ などと下の3式が出てくるのはなぜか説明していただけませんか。 (1)ωとΩはそれぞれ何を表しているのでしょうか。 (2)また、d↑i'/dt=↑ω×↑i'となるのはなぜですか。 (3)なぜ行列式が出てきてこのような表現になるのですか? (4)d↑i'/dtというのはどういうことを表現しているのですか? 質問が多いですが教えてくださると幸いです。 拡大画像↓ http://kie.nu/uF8

質問者が選んだベストアンサー

  • ベストアンサー
回答No.3

Ωは定数と仮定します. (1)XYZ系で↑ω=(0,0,Ω)=Ω↑k'(つまりΩは↑ωのZ成分であり,大きさでもあります) (2)XYZ系はxyz系に対し角速度↑ωで回転運動する運動座標系です.任意のベクトル↑aに対して公式 d↑a/dt=d*↑a/dt+↑ω×↑a があります.ここで ↑a=a_x↑i+a_y↑j+a_z↑k=a_X↑i'+a_Y↑j'+a_Z↑k' とするとき, d*↑a/dt=(da_X/dt)↑i'+(da_Y/dt)↑i'+(da_Z/dt)↑k' は運動座標系XYZからみた↑aの時間変化です.↑a=↑i=↑i'+0↑j'+0↑k'とすると,これはXYZ系での定点X=1,Y=Z=0を表しますから当然d*↑a/dt=↑0です.すなわち d*↑i'/dt=(d1/dt)↑i'+(d0/dt)↑i'+(d0/dt)↑k'=↑0 したがって d↑i'/dt=↑ω×↑i' となります. (3)これは外積を直角座標系の成分で表すための便宜的な表現です.第1行に関して展開しなさいという意味です. (4)(2)でも示した通り,回転座標系XYZにおける点(1,0,0)を静止座標系xyzでみるとこの点が時間的に変化しているということを表します.↑i'を静止座標系xyzの成分で ↑i'=x↑i+y↑j+z↑k とかくと(1)で示したように↑ω=Ω↑k'=Ω↑k(↑k'=↑kに注意,今の場合Z軸はそのままz軸)ですから, d↑i'/dt=(dx/dt)↑i+(dy/dt)↑j+(dz/dt)↑k ↑ω×↑i'=Ω↑k×(x↑i+y↑j+z↑k)=Ω{x(↑k×↑i)+y(↑k×↑j)+z(↑k×↑k)}=Ω(x↑j-y↑i) ∴(dx/dt)↑i+(dy/dt)↑j+(dz/dt)↑k=Ω(x↑j-y↑i) ⇔dx/dt=-Ωy,dy/dt=Ωx,dz/dt=0 ⇔d(x+iy)=iΩ(x+iy),z=z(0)(定数) ここにiは虚数単位です.複素数x+iyは第1式からx+iy={x(0)+iy(0)}e^{iΩt}です. x(0)=1,y(0)=z(0)=0 ですから,x+iy=e^{iΩt},z=0すわなち ↑i'=cosΩt↑i+sinΩt↑j=(cosΩt,sinΩt)(xyz系から) つまり,静止座標xyz系から運動座標XYZ系のX方向の基本ベクトルをみると,角速度Ωでz=Z軸の周りを回転していることがわかります. ※http://www.sp.u-tokai.ac.jp/~yasue/ffn/riki-12.pdf

全文を見る
すると、全ての回答が全文表示されます。

その他の回答 (3)

回答No.4

>なぜこの粒子の運動方程式ではコリオリ力と遠心力の項が出てくるのでしょうか? この式の右辺は,木星と一緒に回転する座標系で記述した加速度の式だからです. >どうしてコリオリ力や遠心力の項と言えるのかわかりません。 詳細は省きますが,遠心力は-mω↑×(ω↑×r↑) (平面運動ではmω^2 r↑),コリオリの力は-2 m ω↑×v'↑なので,角速度と(回転系の)速度の積になっているのがコリオリ項,角速度の2乗になっているのが遠心力項です. >第1,2,3項は何を表現している項なのでしょうか? そのまま回転している座標系から見た加速度です. 質問されている式は,左辺のdV/dtが静止した座標系から見た加速度で, (静止系の加速度)=(回転系の加速度)+(コリオリ項)+(遠心力項) という式です.

全文を見る
すると、全ての回答が全文表示されます。
回答No.2

(1) ω↑は角速度ベクトル。回転軸の方向を向き、大きさが角速度のベクトルです。 質問の図ではk↑ベクトルの方向(z軸)を向いた大きさがΩのベクトルで、 z軸回りに角速度Ωで回転している事を意味します。 (2)(4) i'は図のx(小文字)方向の単位ベクトルで、木星の運動とともに中心星まわりに角速度Ωで回転しているということのようですね。i'が時間変化するのでi'の時間微分も必要になります。 (3) 外積の計算方法の一つです http://kagennotuki.sakura.ne.jp/la/node16.html

happy_lucky3368
質問者

補足

回答ありがとうございます。なんとなくわかった気がします。 http://oshiete.goo.ne.jp/qa/7760279.html の質問を締め切ってしまったのですが、なぜこの粒子の運動方程式ではコリオリ力と遠心力の項が出てくるのでしょうか?本当に物理が苦手なのでよくわかりません。どうしてコリオリ力や遠心力の項と言えるのかわかりません。 また、第4,5項がコリオリ力で第6,7項が遠心力なのだとしたら、第1,2,3項は何を表現している項なのでしょうか?

全文を見る
すると、全ての回答が全文表示されます。
回答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 この式から角振動数ωを求めるにはどのようにすればよいでしょうか

プリントできずに困っています
このQ&Aのポイント
  • PCからのプリントがブルーの印字で飛んで印刷され、コピーはモノクロが白紙状態で印刷できません。
  • 問題のある製品名はMFC-J877Nです。
  • Windows10のパソコンで有線LAN接続しています。関連するソフト・アプリは特にありません。電話回線は不明です。
回答を見る