YAMAGATA UNIVERSITY
回転座標系
と
慣性力
YAMAGATA UNIVERSITY
ちょっとした思考実験
x
y
① 質点が先端につけられている単振り子
振り子をゆらした時、床面に投影された質点の影の運動は?
YAMAGATA UNIVERSITY
x
y
振り子運動の軌跡
① 質点が先端につけられている単振り子
振り子をゆらした時、床面に投影された質点の影の運動は?
YAMAGATA UNIVERSITY
振り子の質点の影: 床が回転していると?
振り子は同じように振動している。
床が振り子の支点を通る鉛直線を軸として回転している場合、 質点の影はどのような軌跡を描く?
x
y
② 振り子の周期と床の回転周期が同じ場合
② 振り子の周期と床の回転周期が同じ場合
YAMAGATA UNIVERSITY
考え方:
回転版
角速度 ω2 で回転
y x = sin sin
1t t ⋅cos ⋅sin
2t t
ω 2 t
r
x
y
( x y ) = ( r r (t)cos ω (t )sin ω
22t t )
r (t)=r sin ω 1 t
原点からの距離 角速度 ω
1 の単振動で変動
YAMAGATA UNIVERSITY
回転版上の振り子の軌跡
1=2 , 2=2 1=2 , 2= 1=2 , 2=2 /3
1=2 , 2=2 /4 1=2 , 2=2 /5 1=2 , 2=2 /6
y x = sin sin
11t⋅cos t⋅sin
22t t
②
回転周期と 振子周期が
同じ
③
回転周期が 振子周期の
半分
YAMAGATA UNIVERSITY
では、地球上で同じ事を考える
北極点にいる観測者から、質点の運動は どの様に観察されるか?
北極点にいる観測者から、質点の運動は どの様に観察されるか?
x
y
北極点に振り子を設置する。 空気抵抗が無視できるように、 振り子の糸は十分長く
質点は十分重い。
北極点に振り子を設置する。 空気抵抗が無視できるように、 振り子の糸は十分長く
質点は十分重い。
YAMAGATA UNIVERSITY
フーコーの振り子
霞城セントラルの振り子は
どんな運動をするのだろうか?
地表にいる観察者は地球の自転を感じない。
観察者にとって、この振り子の先端に結びつけられた 質点の力学はどのように理解されるのか?
地表にいる観察者は地球の自転を感じない。
観察者にとって、この振り子の先端に結びつけられた 質点の力学はどのように理解されるのか?
霞城セントラルの フーコーの振り子 霞城セントラルの フーコーの振り子
ジャン・ベルナールレオン・フーコー
Wikipediaより
YAMAGATA UNIVERSITY
二つの座標系: 『静止座標系』と『回転座標系』
振り子が固定されている座標系 x, y
x
x ' y
y '
回転板に固定された座標系 x', y'
静止座標系 座標系回転
YAMAGATA UNIVERSITY
回転座標系
x
y
x '
y '
二つの座標系の関係は
x ' y ' = −sin cos cos sin
x
y
i = 1
0
j= 0
1
単位ベクトルで確認
(x, y) 系 → (x', y') 系
⃗i = ( 1
0 ) → (
cos θ
−sin θ )
j= 0
1
sin
cos
④
⑤
⃗r' = ( x '
y ' )
⃗r= ( x
y )
YAMAGATA UNIVERSITY
質点の速度
x ' y ' = −sin cos cos sin
x
y
d
dt =
一定質点の速度を求める
→ 上の式の両辺を時間で微分
x
y
x '
y '
⃗r' = ( x '
y ' )
⃗r= ( x
y )
今回は、等速回転している場合を考える 一般的には
座標の回転角度も時間に依存
θ(t)
YAMAGATA UNIVERSITY
回転座標系での見かけの速度
˙x ' ˙y ' = − ˙ cos − ˙sin −˙ sin ˙cos
x
y
cos sin
−sin cos ˙x ˙y
x ' y ' = −sin cos cos sin
x
y
d f t g t
dt =
d f t
dt gt f t dgt dt d f g t
dt =
d f t dt
dg t dt
sin /2=cos cos /2=−sin
( ˙x ' ˙y ' ) =− ω ( cos(θ+π /2) sin(θ+π /2)
−sin (θ+π/2 ) cos( θ+π /2 ) ) (
x
y ) + (
cosθ sin θ
−sin θ cosθ ) ( ˙x ˙y )
θ 回転 (θ + π/2) 回転
⑥ ⑦
⑧ ⑨
両辺時間微分
YAMAGATA UNIVERSITY
回転変換行列
⃗r' = R (θ) ⃗r
x
y
x '
y ' ⃗r= ( ⃗r' = ( x ' y ' )
x
y )
x ' y ' = −sin cos cos sin
x
y
※ 角度θだけ座標を回転させる変換
R(θ) の性質(和・積、微分)を確認してみる
YAMAGATA UNIVERSITY
微分微分 足し算足し算
回転変換行列の演算
˙R =−R /2 = R − /2
二階微分二階微分
¨R =
2R =−
2R
˙R =−R R /2= R R − /2
x
y
x '
y '
ϕ
y ' ' x ' ' ⃗r ' '=R(ϕ)⃗r '
⃗r ' =R(θ)⃗r
⃗r ' '=R(θ+ϕ)⃗r
⃗r ' ' =R(ϕ)R (θ)⃗r
R (θ+ϕ)=R (θ)R(ϕ)
˙R(θ)=
(
−sin θ cosθ−cosθ −sin θ
)
˙θR (θ)= ( cos θ sin θ
−sin θ cos θ )
˙R(θ)=−
(
cos(θ+π/ 2) sin (θ+π /2)−sin (θ+π/2) cos(θ+π /2)
)
ωYAMAGATA UNIVERSITY
回転座標系からみた速度
⃗r' = R (θ) ⃗r
x
y
x '
y ' ⃗r= ( ⃗r' = ( x ' y ' )
x
y )
x ' y ' = −sin cos cos sin
x
y
速度の変換 速度の変換
˙⃗r' =ω R ( θ−π/2 )⃗r+R(θ) ˙⃗r
˙⃗r' =R(θ) [ ω R(−π/ 2)⃗r+ ˙⃗r ]
˙⃗r' = ˙R ⃗r+R (θ) ˙⃗r ˙R(θ)=ω R( θ−π/2 )
YAMAGATA UNIVERSITY
方程式の意味する事
˙r' =R [ R− /2r ˙r ]
r' =R [ r ]
x
y
x ' y '
r
(x,y) 座標系での位置 (x',y') 座標系への変換
大きさが ωr で、位置ベクトルを -90度 回転。 (x',y') 座標系
での位置
(x',y') 座標系での速度 (x,y) 座標系
(x',y') 座標系への変換 での速度
この項が余分: みかけの速度
YAMAGATA UNIVERSITY
みかけの速度: 回転座標系での質点の運動
(x,y) 系で静止している物体
˙r=0
⑩(x',y') 系では
˙r' =R R − /2r
⑪x
y
x '
y '
r
速度の大きさ
v ' =∣˙⃗r '∣=ω r
⑫※ 並進する座標系での話と同じ
˙r' =R [ R − /2r ˙r ]
逆回転に速度 v' で 等速回転するように 見える
YAMAGATA UNIVERSITY
回転座標系で静止している質点
˙r' =0
回転座標系で静止
R − /2r ˙r=0
˙⃗r=ω R (π /2)⃗r
位置ベクトルを90°回転した方向に大きさ rω の速度を持つ。
˙r' =R [ R − /2r ˙r ]
コマの上に書いた『点』は、
コマが角速度ωで等速回転している時どのように見えるか?
R (−π /2)⃗r =− R (π/2)⃗r
∣ ˙⃗r ∣ = r ω
おおきさ
x '
y '
⃗r'
θ=ω t
r ω
角速度ωで等速回転しているように見える
YAMAGATA UNIVERSITY
二階微分: 等速回転系
r' =R r
˙⃗r' =R(θ) [ ω R(−π/ 2)⃗r+ ˙⃗r ]
¨⃗r' = ˙R(θ) [ ω R(−π /2)⃗r+ ˙⃗r ] +R (θ) [ ω R(−π /2) ˙⃗r+ ¨⃗r ]
¨⃗r' =R(θ)ω R (−π/2) [ ω R(−π /2)⃗r+ ˙⃗r ] +R (θ) [ ω R (−π/2) ˙⃗r+ ¨⃗r ]
¨⃗r' =R(θ) [ −ω
2⃗r+ω R (−π/2) ˙⃗r ] +R(θ) [ ω R(−π /2) ˙⃗r+ ¨⃗r ]
¨⃗r' =R(θ) [ −ω
2⃗r+2 ω R (−π/2) ˙⃗r + ¨⃗r ]
˙R(θ)=ω R(θ−π/2)=ω R (θ)R (−π/2)
R(θ+π)=R(θ)R(ϕ) R(−π)⃗r=−⃗r
YAMAGATA UNIVERSITY
(等速)回転座標系での運動方程式
x
y
x '
y '
r
−m
2r r' =R r
¨r' =R [ −
2r2 R− /2 ˙r ¨r ]
˙r' =R [ R− /2r ˙r ]
m ¨ r ' =R [ −m
2r 2 m R − /2 ˙r m ¨r ]
“みかけの”向心力
物体が (x, y) 系で 速度を持つと現れる
F ' ⃗ ≠R(θ) ⃗ F
F ⃗ =m ¨⃗r
F ' ⃗ =m ¨⃗r' '
力の単純な回転ではない これは一体なんだろうか?
YAMAGATA UNIVERSITY
等速回転座標系からみた質点: 質点が等速運動している場合
x
y
x '
y '
r
−m
2r
v
2 mR − /2r
“みかけの”向心力 質点が運動している時の みかけの力
F ' =R [ −m
2r 2 m R− /2 ˙r F ]
例) 原点から遠ざかる質点 例) 原点から遠ざかる質点
v ∝r
速度に対して垂直方向(-90度)に 力が作用しているように見える
見かけの力によって、進行方向に対して
見かけの力によって、進行方向に対して
速度と垂直方向の慣性力
コリオリの力
YAMAGATA UNIVERSITY
振り子の質点の影: 『コリオリの力』
振り子は同じように振動している。
床が振り子の支点を通る鉛直線を軸として回転している場合、 質点の陰はどのような軌跡を描く?
x
y
② 振り子の周期と床の回転周期が同じ場合
② 振り子の周期と床の回転周期が同じ場合
YAMAGATA UNIVERSITY
『コリオリの力』
『コーヒーカップ』
に乗りながら、キャッチボールをすると・・・
共に回転している人からは、ボールが 進行方向右向きに曲がるように見える
YAMAGATA UNIVERSITY
角速度ベクトルを使ってあらわす
x
y
θ=ω t
− ×r
˙⃗r' =R(θ) [ ω R(−π /2)⃗r + ˙⃗r ]
位置ベクトルを 角速度ベクトルを軸に -90°回転させて、角速度をかける。
− ω ⃗ ×⃗r
角速度ベクトル ω をもつ質点の速度
静止している質点を角速度 ω の回転系から見たみかけの速度 角速度ベクトル ω をもつ質点の速度
静止している質点を角速度 ω の回転系から見たみかけの速度
×r
− ×r=r×
R − /2 f − × f
一般的に
×r
位置ベクトルと垂直 角速度ベクトルと垂直
r と ω の外積 =⃗r× ω ⃗
r
YAMAGATA UNIVERSITY
外積による表現との対応
˙r' =R [ R− /2r ˙r ] ˙r' =R [ − ×r ˙r ]
¨r' =R [ −
2r2 R− /2 ˙r ¨r ]
¨r' =R [ −
2r −2 × ˙r ¨r ]
3次元ベクトルに拡張しても基本的には同じ
r' =R r r ' =R r
R − /2 f − × f
F ' =R [ −m
2r2 m R− /2 ˙r F ]
F ' =R [ −m
2r −2 m × ˙r F ]
コリオリの力
YAMAGATA UNIVERSITY
地表面に固定された座標系
山形市の緯度 北緯 38 度
sin 38° = 0.62 cos 38° = 0.79
地球の平均半径 R = 6.4 × 106 m
→@山形市では R cos 38° = 5.0×106 m 自転の周期 T = 24 h = 1440 min = 86400 s 自転の角速度 ω = 2π/(86400) = 7.3×10-5 s-1
'
0コリオリの力を考える:
鉛直方向の動きを無視できる場合、
『緯度 θ0 の位置に固定された座標系』は
『角速度 ω'=ωsinθ0 の回転座標系』 として近似できる。
YAMAGATA UNIVERSITY
台風はなぜいつも進行方向に対して右に曲がるのか?
地球の外側の系で直進(北上)する「台風」は 地表に固定された座標系(回転座標系)
からは曲がって見える。
= 見かけの力により曲げられたように見える 地球の外側の系で直進(北上)する「台風」は 地表に固定された座標系(回転座標系)
からは曲がって見える。
= 見かけの力により曲げられたように見える 見かけの力
F ' =R [ 2 m R− /2sin
0˙r ]
'
YAMAGATA UNIVERSITY
単振り子
−mg
F ∼−mg sin θ
L
x
F ' =−mg sin θ
重力の弦に垂直な成分は
−mg
角度 θ が十分小さい時は、水平方向の 力の大きさに等しい
F ∼− mgx
L
x
L
振動係数 mg/L の単振動に等しい
振動係数 mg/L の単振動に等しい
m ¨x=− mg
L x
YAMAGATA UNIVERSITY
フーコーの振り子
フーコーの振り子(1851年・フーコー)
弦の長さ L = 67 m
おもりの質量 m = 28 kg
周期 T = 16.4 s
フーコーの振り子@霞城セントラル
弦の長さ L = 20 m
おもりの質量 m = 32 kg 周期 T = 9.0 s
単振り子の周期 T=2
Lg
重力加速度 g=9.8 m/s2
⑬
⑭
YAMAGATA UNIVERSITY
等速回転系から見た単振り子
−mg
F ∼− mg
L x
L
x
m ¨ r=−k r−2 m × ˙r
バネ定数 k の単振動 地球の自転によるコリオリ力
(角速度 ω の回転運動) 北緯 θ0 では、角速度 ωsinθ0 と近似
振り子は、
・ 周期 T = 2π√(L/g) の単振動
・ 原点周りに 角速度 ωsinθ0 で回転 地表に設置された振り子の運動
地球の自転の周期が 24時間 とすると、
北緯 θ0 での振り子のコリオリ力による周期は 24/sinθ0 時間
F ⃗ (⃗r)=−k ⃗r (k =mg/ L)
YAMAGATA UNIVERSITY
1時間あたりの単振り子の振動数
周期 T = 9.0 s
60×60/9.0 = 400回
一回転するのに必要な振動数
38.8 h × 400 ~ およそ15000回
フーコーの振り子@霞城セントラル
山形(北緯38°)での地球の自転による コリオリ力
角速度 ω sin 38° 1時間あたりの回転角度
(360°/24)×sin 38° = 9.3°
周期(時間) 24 h /sin38° = 38.8 h
⑯
⑰
⑱
⑮
⑮
⑮
⑰ ⑱ ⑲
YAMAGATA UNIVERSITY
フーコーの振り子@高速回転する星
振り子の軌跡を書くのは無謀 霞城セントラルの振り子が、 角速度 0.01 s-1 で回転する 星の上に置かれていた時
(重力加速度は同じ)の軌跡は
YAMAGATA UNIVERSITY
等速回転系からみた単振り子
m ( ⃗r× ¨⃗r ) =−k ⃗r×⃗r−2 m ( ⃗r× ( ω× ˙⃗r ⃗ ) )
d
dt
r× ˙r
=r× ¨r ddt
(
⃗r× ˙⃗r)
= ˙⃗r× ˙⃗r+⃗r× ¨⃗r ddt
(
⃗r× ˙⃗r)
m ¨ r=−k r−2 m × ˙r
d
dt
r×
×r
=2
r ×
× ˙r
ddt
(
⃗r×(
ω×⃗r⃗) )
= ˙⃗r×(
ω×⃗r⃗)
+⃗r×(
˙⃗ω×⃗r)
+⃗r×(
ω× ˙⃗r⃗)
˙=0
˙⃗r×
(
ω×⃗r⃗)
=−⃗ω×(
⃗r× ˙⃗r)
−⃗r×(
˙⃗r×⃗ω)
r× ˙r ∝
˙r×
×r
=r×
× ˙r
d
dt r× ˙r =−
d
dt r× ×r
r× ˙r =− r× ×r
角速度 ω での回転
(角速度成分のみに注目) 両辺の、位置ベクトルとの外積を考える
d
dt
(
⃗r×(
ω×⃗r⃗) )
˙⃗r×