YAMAGATA UNIVERSITY
剛体の力学
コマはなぜ倒れないのか?
YAMAGATA UNIVERSITY
色々なコマ
コマの軸は『枠』に接している
・ 歳差運動
・ 自立する(極端なものは 『眠りコマ』)
・ 起き上がるものもあるし、不思議な動きをするものもある
・ 歳差運動
・ 自立する(極端なものは 『眠りコマ』)
・ 起き上がるものもあるし、不思議な動きをするものもある
回転体の力学 → 『地球ゴマ』 コマは 『枠』 の中に設置。
コマの回転軸と、床などが直接接しない。 回転体の力学 → 『地球ゴマ』
コマは 『枠』 の中に設置。
コマの回転軸と、床などが直接接しない。
YAMAGATA UNIVERSITY
コマの運動: 歳差運動
z
x
y
固定点
I
xx=I
yy= I I
zzコマには 『安定な回転軸』 がある。
重力
重力によるモーメントが発生
コマを水平な床の上に置く 重力が働く
重力のモーメントによりコマはどのように運動するか?
重力のモーメントによりコマはどのように運動するか?
N ⃗
N ⃗
YAMAGATA UNIVERSITY
水平に倒した地球ゴマの運動
固定点
z
x y
地球ゴマの回転軸を水平に倒した場合 糸で吊す
棒で支える
回転軸の一端を空間に『固定』する
重力によるモーメント 重力
N ⃗
N ⃗
コマに張り付けた軸(慣性主軸)の動きを
静止系から考える
この条件の下で、コマの運動を
YAMAGATA UNIVERSITY
水平に固定したコマの運動
⃗er
l
N
M ⃗g ω⃗ =ω ⃗er
N ⃗ =l M ⃗e
r× ⃗g
コマが高速で回転
⃗L = I ω ⃗e
r固定点と重心との距離 l
回転の角速度
ω ⃗ = ω ⃗e
r 回転軸の方向の単位ベクトル⃗e
r
コマの質量 M
によるモーメント
コマの慣性モーメント I
角運動量 重力
YAMAGATA UNIVERSITY
上から見ると
ジャイロスコープ効果
d
⃗er
l
N
M ⃗g ω⃗ =ω ⃗er
N ⃗ =l M ⃗e
r× ⃗g
重力によるモーメント
⃗L = I ω ⃗e
r角運動量
⃗L
N ⃗
角運動量の変化 角運動量の変化
d ⃗L
dt =N⃗
d ⃗ L = N ⃗ dt
d ⃗ L
角運動量の回転の 『角速度』 角運動量の回転の 『角速度』
= ∣⃗ N ∣dt
∣⃗L∣ ω '=
d ϕ
dt =
∣⃗ N ∣
∣⃗L∣
モーメントと角運動量が直交 モーメントと角運動量が直交
『ジャイロスコープ効果』
→ 角運動量は大きさを変えず、円運動
d ϕ= ∣d ⃗L∣
∣⃗L∣
YAMAGATA UNIVERSITY
ジャイロスコープ効果による歳差運動の角速度
d
⃗L
N ⃗
d ⃗ L
N ⃗ =l M ⃗e
r× ⃗g
重力によるモーメント
⃗L = I ω ⃗e
r角運動量
角運動量の回転の 『角速度』 角運動量の回転の 『角速度』
= ∣⃗ N ∣dt
∣⃗L∣ ω '=
d ϕ
dt =
∣⃗ N ∣
d ϕ= ∣d ⃗L∣ ∣⃗L∣
∣⃗L∣
ω '= ∣⃗ N ∣
∣⃗L∣ =
l M g
I ω =
l M
I g
ω 1
コマの角速度コマの形状
歳差運動の角速度が決まる 歳差運動の角速度が決まる
YAMAGATA UNIVERSITY
地球ゴマを使った実験
歳差運動の周期(秒)
歳差運動の角速度(radian/s)
ω ' = T 2 π
T
ω '= l M g
I ω
運動方程式から
測定により決定
測定結果と比較
YAMAGATA UNIVERSITY
地球ゴマの質量・サイズ
コマ 直径 枠 合計
大 60.1 g 5.5 cm 18.8 g 78.9 g 小 26.6 g 4.0 cm 9.9 g 36.5 g
コマ 直径 枠 合計
大 60.1 g 5.5 cm 18.8 g 78.9 g 小 26.6 g 4.0 cm 9.9 g 36.5 g
ほとんどが円板のの質量だと仮定 コマの軸周りの慣性モーメント
大 60 g × (2.8 cm)2 / 2 = 235 g cm2 小 26.6 g × (2.0 cm)2 / 2 = 53.2 g cm2 ほとんどが円板のの質量だと仮定
コマの軸周りの慣性モーメント
大 60 g × (2.8 cm)2 / 2 = 235 g cm2 小 26.6 g × (2.0 cm)2 / 2 = 53.2 g cm2
半径 r0、質量 M の円板の慣性モーメント I は
I = 1
2 M r
02
YAMAGATA UNIVERSITY
歳差運動の角速度を求める
支点から重心までの距離 大 l ~ 4.5 cm
小 l ~ 3.2 cm
支点から重心までの距離 大 l ~ 4.5 cm
小 l ~ 3.2 cm
コマの回転の角速度 ω
こまの軸の直径 4.5 mm
→ 円周は π×4.5 mm ~ 14 mm こまを回すのにつかった紐の長さ 約 60cm
紐を引くと軸は 600 / 14 = 42.9 回転 紐を引くのにかかる時間 ~1s
軸の角速度 ω = 42.9 × 2π / (1 s) ~ 270 radian/s コマの回転の角速度 ω
こまの軸の直径 4.5 mm
→ 円周は π×4.5 mm ~ 14 mm こまを回すのにつかった紐の長さ 約 60cm
紐を引くと軸は 600 / 14 = 42.9 回転 紐を引くのにかかる時間 ~1s
軸の角速度 ω = 42.9 × 2π / (1 s) ~ 270 radian/s
コマ 直径 枠 合計 慣性モーメント
大 60.1 g 5.5 cm 18.8 g 78.9 g 235 g cm2 小 26.6 g 4.0 cm 9.9 g 36.5 g 53.2 g cm2
コマ 直径 枠 合計 慣性モーメント 大 60.1 g 5.5 cm 18.8 g 78.9 g 235 g cm2 小 26.6 g 4.0 cm 9.9 g 36.5 g 53.2 g cm2
ω '= l M g
I ω
それぞれのコマの 歳差運動の角速度 を求める
YAMAGATA UNIVERSITY
コマの歳差運動: 斜めの場合
N ⃗
⃗L = I ω ⃗
N ⃗ =⃗l× W ⃗
⃗L= I ω ⃗
dL = N ⃗ dt
W ⃗
コマの軸が角度 θ 傾いている場合 慣性モーメント I
角速度 ω
横から
上から
重力によるモーメント |N| = lWsinθ
角速度ベクトルの変化
dL = N ⃗ dt
歳差運動の角速度ω '= l W sin θ
I ω
ω '= d ϕ
dt =
∣⃗ N ∣
∣⃗L∣
YAMAGATA UNIVERSITY
剛体の運動エネルギー
並進運動による運動エネルギー
R
r
i'r
iK = 1
2 ∑
i=1N
m
i ˙r
i
2K = 1
2 ∑
i=1N
m
i ˙R ˙r
i'
2K = 1
2 M ˙⃗ R
2
+ ( ∑
i=1 N
m
i˙⃗r
i') ⋅ ˙⃗R+ ∑
i=1
N
1
2 m
i˙⃗r
i' 2
=0
K = 1
2 M ˙⃗ R
2
回転による運動エネルギー 回転による運動エネルギー
T = 1
2 ∑
i=1N
m
iv
i⋅v
iT = 1
2 ∑
i=1N
m
i ×r
i ⋅ ×r
i
T = 1
2 ⋅ ∑
i=1N
m
ir
i× ×r
i
=⃗L
iT = 1
2 ⋅L
回転していなければ 回転していなければ
˙⃗r '=0
YAMAGATA UNIVERSITY
剛体の運動: まとめ
1
2 M ˙⃗ R
2 1
2 ω⋅⃗L ⃗
重心にある質量Mの 質点の並進運動
重心にある質量Mの 質点の並進運動
運動エネルギー 運動エネルギー
重心まわりの回転運動 重心まわりの回転運動
運動方程式 運動方程式
= p
2
˙⃗p= ⃗F ˙⃗L= ⃗N
YAMAGATA UNIVERSITY
r
0コマの運動: 3つの角度と回転軸
⃗n
⃗e
ze
xl
r
0⃗e
x⃗e
yψ
ω= ⃗
ϕ
⃗j= ⃗e
z× ⃗n
∣ ⃗e
z× ⃗n ∣
⃗j
⃗e
z⃗n
上から見た図 上から見た図
横から見た図
横から見た図 回転軸に沿って見た図回転軸に沿って見た図
水平面内 の回転
+ ˙ϕ ⃗n
回転軸まわりの高速回転
˙ψ ⃗e z
上下方向への回転
+˙θ ⃗j
慣性主軸慣性主軸
⃗e
ze
xe
y⃗n
⃗e
z⃗j
YAMAGATA UNIVERSITY
慣性主軸に角速度を分解 ω= ˙ψ ⃗ ⃗e
z
+ ˙θ ⃗j + ˙ϕ ⃗n
⃗n
e
zl
r
0⃗j
⃗n=sin θ⃗ez+cos θ ⃗e'
e
xe
y⃗e'
ψ⃗j
⃗e '=cos ψ⃗ex−sin ψ⃗ey
ω= ⃗ ( ˙θsin ψ+ ˙ϕ cosθcos ψ ) ⃗e
x+ ( ˙θcos ψ− ˙ϕ cosθsin ψ ) ⃗e
y⃗j=sin ψ⃗ex+cos ψ⃗ey
⃗n=sin θ⃗ez+cos θ
(
cos ψ⃗ex−sin ψ⃗ey)
⃗e'
YAMAGATA UNIVERSITY
角速度ベクトル → 角運動量ベクトル
回転体の慣性主軸 回転軸を z 軸
z 軸に垂直に x, y 軸
慣性乗積 0 慣性テンソル
Izz ( x 軸) I = Ixx = Iyy( y, z 軸) 回転体の慣性主軸
回転軸を z 軸
z 軸に垂直に x, y 軸
慣性乗積 0 慣性テンソル
Izz ( x 軸) I = Ixx = Iyy( y, z 軸)
⃗L= I ω ⃗
= ( ˙θsin ψ+ ˙ϕcosθ cos ψ ) I ⃗e
x+ ( ˙θcos ψ− ˙ϕ cosθsin ψ ) I ⃗e
y+ ( ˙ψ+ ˙ϕsin θ ) I
zz⃗e
zω= ⃗ ( ˙θsin ψ+ ˙ϕ cosθcos ψ ) ⃗e
x+ ( ˙θcos ψ− ˙ϕ cosθsin ψ ) ⃗e
y+ ( ˙ψ+ ˙ϕsin θ ) ⃗e
zYAMAGATA UNIVERSITY
力学的エネルギーの保存
ω=⃗ (˙θsin ψ+ ˙ϕ cosθcos ψ)⃗ex ω=⃗ +(˙θcos ψ− ˙ϕcos θsin ψ)⃗ey ω=⃗ +( ˙ψ+ ˙ϕ sin θ)⃗ez
L=I
⃗L=(˙θsin ψ+ ˙ϕ cosθcos ψ)I ⃗ex
⃗L+(˙θcos ψ− ˙ϕcosθ sin ψ)I ⃗ey
%L+⃗ ( ˙ψ+ ˙ϕ sin θ)Izz⃗ez
E = 1
2 ( ( ˙ψ+ ˙ϕsin θ )
2I
zz+ ( ˙θ
2+ ˙ϕ
2cos
2θ ) I ) +Mgl sin θ
=
z2=
2x
2y力学的エネルギーの保存より 力学的エネルギーの保存より
⃗n
e
zl
r
0⃗j
力学的エネルギー 力学的エネルギー
E =
回転の運動エネルギー
1
2 ω⋅⃗L ⃗
円板の位置エネルギー
+Mgl sin θ
YAMAGATA UNIVERSITY
剛体の運動方程式
n
ez
ex
l r0
O
j
M g
N ⃗ =(l ⃗e
z)×(M ⃗g)
重力による支点まわりのモーメント
N⃗ ⊥⃗ez N⃗ ⊥ ⃗n
空間に固定された系 n では ˙Li=Ni N z=0 Nn=0
˙Ln=0 水平周りの回転は等角速度
=(l ⃗e
z)×(−M g ⃗n )
コマに固定された系(慣性主軸) コマに固定された系(慣性主軸)
L
n=L⋅n
一定
ω=⃗ (˙θsin ψ+ ˙ϕ cosθcos ψ)⃗ex ω=⃗ +(˙θcos ψ− ˙ϕcos θsin ψ)⃗ey ω=⃗ +( ˙ψ+ ˙ϕ sin θ)⃗ez
L=I
⃗L=(˙θsin ψ+ ˙ϕ cosθcos ψ)I ⃗ex
⃗L+(˙θcos ψ− ˙ϕcosθ sin ψ)I⃗ey
%L+⃗ ( ˙ψ+ ˙ϕ sin θ)Izz⃗ez
L
n= I
zz( ˙ψ+ ˙ϕ sin θ ) sin θ+I ˙ϕ cos
2θ
N ⃗
YAMAGATA UNIVERSITY
剛体の運動方程式
n
ez
ex
l r0
O
j
M g
剛体に固定された系 ez では 剛体に固定された系 e
z では
オイラーの方程式
N
z= I
zz ˙
z I
yy− I
xx
y
x=0
˙ω z =0
z 軸に等角速度で回転 I = Ixx = Iyy N⃗ ⊥⃗ez N z=0一定
ω=⃗
(
˙θsin ψ+ ˙ϕ cosθcos ψ)
⃗ex ω=⃗ +(
˙θcos ψ− ˙ϕcos θsin ψ)
⃗eyω=⃗ +
(
˙ψ+ ˙ϕ sin θ)
⃗ez z=
˙ ˙sin
YAMAGATA UNIVERSITY
三個の保存量
E=1
2
(
( ˙ψ+ ˙ϕ sin θ)2 Izz+(
˙θ2+ ˙ϕ2cos2θ)
I)
+ Mgl sin θLn= I zz
(
˙ψ+ ˙ϕ sin θ)
sinθ+I ˙ϕ cos2θωz=
(
˙ψ+ ˙ϕsin θ)
初期条件 ( t = 0 )初期条件 ( t = 0 ) θ=0
˙ϕ=0
回転軸が水平
水平方向に回転していない
n
ez
ex
l r0
O
j
L
n(t=0)= I
zzω
zsin (θ)+ I ˙ϕ cos
2(θ)
力学的エネルギー
水平回転の角運動量
コマの回転軸周りの角速度
L
n(t=0)=0 ˙ϕ= −I
zzω
zsin θ
I cos
2θ
YAMAGATA UNIVERSITY
水平方向への回転
ex
ey
ψ
n
ϕ
n
ez j
n
ez
ex
l r0
O
j
=0 ˙ϕ=0
n 軸まわりの角速度
回転軸の水平からの傾きが大 → 角速度 大 回転軸まわりの角速度大 → 角速度 大
˙ϕ= −I
zzω
zsin θ
I cos
2θ
角速度なし → 回転しない
水平の場合は角速度が生まれない
YAMAGATA UNIVERSITY
回転軸の水平に対する角度の時間的変化
E=1
2
(
( ˙ψ+ ˙ϕsin θ)2 Izz+( ˙θ+ ˙ϕcosθ)2I)
+Mgl sin θエネルギー・n軸周りの角運動量は保存
Ln= I zz
(
˙ψ+ ˙ϕ sin θ)
sinθ+I ˙ϕ cos2θI ¨ θ=− ( M g l −I
zzω
z˙ϕ+I ˙ϕ
2sin θ ) cos θ
ωz: 一定 ωz: 一定
→ 時間微分すると 0
I ˙ θ ¨θ+ I ˙ϕ ¨ϕcos
2θ− I ˙ϕ
2˙θcosθsin θ+ M g l ˙θ cosθ=0
I ¨ ϕcos
2θ−2 I ˙θ ˙ϕ cos θsin θ+I
zzω
z˙θcos θ=0
YAMAGATA UNIVERSITY
コマの歳差運動
ex
ey
n
n
ez j
n
ez
ex
l r0
O
j
I ¨ θ=− ( M g l −I
zzω
z˙ϕ+I ˙ϕ
2sin θ ) cos θ
˙ϕ= −I
zzω
zsin θ
I cos
2θ
数値計算の例
YAMAGATA UNIVERSITY
参考資料
●
比較的簡単な読み物
●
「コマの不思議」 黒須 茂 著
– 山文社
– ISBN-10: 4879260762
– ISBN-13: 978-4879260765
●
演習問題として、コマの力学を解いてみたい人向け
●
「こまはなぜ倒れないか」 大槻義彦・小牧研一郎 編
– 物理学 One Point、 共立出版株式会社
– ISBN-10: 4320033620
– ISBN-13: 978-4320033627