先日公開した半円筒面からの滑走・離脱・放物運動シミュレータでは、半径・質量・摩擦係数・初速をスライダーで変えながら、物体が半円筒のてっぺんから滑り出し、途中のある角度でふわっと面から離れて放物線を描く様子をその場で確認できます（関西大学2026年度入学試験「物理」の設問を着想源に、独自の物理設定・数式で実装した教材です）。
発射する前に「半径を2倍にしたら、面から離れる角度は変わると思う？」と予想させる機能があるのですが、これが地味に面白くて、直感だと「大きい円のほうが遅く離れそう」という気がしてしまいます。実際にやってみると、離れる角度はまったく変わりません。今日はこの「なぜその瞬間に物体は面から離れるのか」を、高校で習うニュートン力学だけを使って一歩ずつ導いてみようと思います。
1. 半円筒面の運動を「エネルギー保存則」から立てる
半径 \(r\) の半円筒（底面が床に接し、中心 \(C\) は床面上）を考えます。質量 \(m\) の小物体を、円筒の最高点 \(Q\)（床から高さ \(r\)）に置き、静かに（初速0で）滑らせます。円筒面はなめらか（摩擦なし）とします。
物体の位置は、中心 \(C\) から見た角度 \(\theta\)（鉛直上方向 \(CQ\) を \(\theta=0\) とする）で表すと見通しがよくなります。このとき、床からの高さは
\[ h(\theta) = r\cos\theta \]
です。\(\theta=0\) で \(h=r\)（最高点 \(Q\)）、\(\theta\) が増えるにつれて高さは下がっていきます。
円筒面はなめらかなので、物体に仕事をする力は重力だけです。力学的エネルギー保存則（運動エネルギー＋位置エネルギーが一定）を、最高点 \(Q\)（\(v=0\)）と角度 \(\theta\) の位置の間で立てると、
\[ \frac{1}{2}mv^2 = mg\bigl(r – r\cos\theta\bigr) \]
右辺は、\(Q\) から角度 \(\theta\) の位置まで下がった高さ \(r(1-\cos\theta)\) ぶんだけ、重力がした仕事です。これを \(v^2\) について解くと、
\[v^2 = 2gr(1-\cos\theta)\]
という、シミュレータの内部でも使われている式が出てきます。質量 \(m\) が両辺から消えるので、速さは物体の重さによらないこともこの時点でわかります。
2. 垂直抗力 N を求める — 向心力の運動方程式
物体は半径 \(r\) の円弧上を動いているので、中心 \(C\) に向かう向き（半径方向）に、大きさ \(v^2/r\) の向心加速度を持っています。半径方向にはたらく力は、重力の半径方向成分 \(mg\cos\theta\)（中心向き）と、面から受ける垂直抗力 \(N\)（中心から離れる向き、面が物体を押し返す力）の2つです。
ニュートンの運動方程式 \(F=ma\) を半径方向に立てると、
\[ mg\cos\theta – N = \frac{mv^2}{r} \]
すなわち
\[ N = mg\cos\theta – \frac{mv^2}{r} \]
となります。ここに第1節で求めた \(v^2=2gr(1-\cos\theta)\) を代入すると、
\[N = mg(3\cos\theta – 2)\]
という、意外にすっきりした式にまとまります。よく見ると、この式には半径 \(r\) が一切登場しません。垂直抗力は物体の重さと角度だけで決まり、円筒がどれだけ大きいかにはまったく関係しないのです。\(\theta\) が増えるにつれて \(\cos\theta\) は単調に減るので、\(N\) も単調に減っていきます。シミュレータの垂直抗力ゲージがすぅっと短くなっていくのは、この式をそのままリアルタイムに描画しているだけです。
3. 離脱角 θ_S を解く — なぜ半径にも質量にもよらないのか
面が物体を「押せる」のは \(N \ge 0\) の間だけです。面は引っ張る力（負の \(N\)）を出せないので、\(N\) がちょうど0になった瞬間、物体は面から離れます。第2節の式で \(N=0\) と置くと、
\[ mg(3\cos\theta_S – 2) = 0 \quad\Longrightarrow\quad \cos\theta_S = \frac{2}{3} \]
これを解くと \(\theta_S \approx 48.19°\) です。この式をもう一度眺めてみると、質量 \(m\) は両辺から消え、半径 \(r\) はそもそも式に現れていません。
離脱角は、物体の重さにも、円筒の大きさにもよらず、常に約48.19°
これが、冒頭の「半径を2倍にしても離脱角は変わらない」の正体です。感覚的には「大きい円のほうがなだらかで、もっと奥まで持ちこたえそう」と思いたくなりますが、\(v^2\) も \(N\) も \(r\) に比例する項と \(r\) を含まない項がちょうど同じ形で打ち消し合うため、\(N=0\) となる角度そのものには \(r\) が残らないのです。離脱時の速さ \(v_S\) と、離脱点の床からの高さ \(h_S\) も、ここまでの式に \(\cos\theta_S=2/3\) を代入するだけで求まります。
\[ v_S = \sqrt{\frac{2gr}{3}}, \qquad h_S = r\cos\theta_S = \frac{2}{3}r \]
こちらは \(r\) に依存します。「離脱する角度」と「離脱するときの速さ・高さ」とで、\(r\) への依存性が違う、という点も面白いところです。
4. 離脱の瞬間の速度と、放物運動へのバトンタッチ
面から離れた瞬間、物体の速度は円の接線方向（中心 \(C\) と物体を結ぶ半径 \(CS\) に垂直な向き）を向いています。この向きを求めるには、位置ベクトル
\[ (x,y) = (r\sin\theta,\ r\cos\theta) \]
を \(\theta\) で微分すればよく、速度の向きは \((\cos\theta, -\sin\theta)\) に比例します。したがって、離脱時の速さ \(v_S\) を水平・鉛直（下向き）に分解すると、
\[V_x = v_S\cos\theta_S, \qquad V_y = v_S\sin\theta_S\]
となります。この式は、\(\theta\to 90°\)（円筒のいちばん側面、床に接する点）で検算できます。そこでは円の接線はちょうど鉛直になるはずですが、代入すると \(V_x=v_S\cos90°=0\)、\(V_y=v_S\sin90°=v_S\) となり、確かに鉛直方向だけの速度になります。逆に \(\theta\to0°\)（てっぺん付近）では接線はほぼ水平になるはずで、こちらも \(V_x\approx v_S\)、\(V_y\approx0\) と整合します。
離脱後は物体を押すものが何もなくなるので、あとは重力だけを受ける斜方投射と同じです。離脱点を原点、離脱の瞬間を \(t=0\) とすると、
\[ x(t) = V_x t, \qquad y_{\text{drop}}(t) = V_y t + \frac{1}{2}gt^2 \]
着地は、離脱点からの落下距離 \(y_{\text{drop}}(t)\) がちょうど離脱点の高さ \(h_S\) に達した瞬間なので、\(y_{\text{drop}}(t)=h_S\) という2次方程式を解けば、着地時刻と水平到達距離が求まります。ここは空気抵抗も摩擦もない、純粋な斜方投射の計算です。シミュレータのキャンバス上で、円運動（青緑の弧）から放物運動（橙の曲線）へキンク（折れ目）なく滑らかに繋がって見えるのは、離脱点での速度の向きが、まさにこの接線方向の式と一致しているからです。
5. 時間の関数として滑走を見る — θ(t) と「動き出せない」パラドックス
ここまでは、角度 \(\theta\) を主役にして速さや垂直抗力を求めてきました。アニメーションとして滑走を再現するには、今度は時刻 \(t\) の関数として \(\theta(t)\) を求める必要があります。
物体が円弧に沿って動く速さは \(v = r\,d\theta/dt\) なので、第1節の \(v^2=2gr(1-\cos\theta)\) と組み合わせると、
\[ \frac{d\theta}{dt} = \sqrt{\frac{2g(1-\cos\theta)}{r}} \]
という微分方程式が得られます。半角公式 \(1-\cos\theta = 2\sin^2(\theta/2)\) を使って整理すると、
\[ \frac{d\theta}{dt} = 2\sqrt{\frac{g}{r}}\,\sin\frac{\theta}{2} \]
変数分離して積分すると（\(\int \csc u\,du = \ln|\tan(u/2)|\) という積分公式を使います）、
\[ \tan\frac{\theta}{4} = \tan\frac{\theta_0}{4}\,\exp\!\left(t\sqrt{\frac{g}{r}}\right) \]
という閉じた式（解析解）にたどり着きます。ここで \(\theta_0\) は積分定数を決めるための初期角です。
ここで、ちょっと気になることが起きます。もし物体が本当に厳密な \(\theta_0=0\)（てっぺん、静止状態）から出発するとしたらどうなるでしょうか。\(\tan(0/4)=0\) なので、右辺は常に0になってしまい、どれだけ時間が経っても \(\theta\) は0のまま。上の式を時刻について解くと、
\[ t = \sqrt{\frac{r}{g}}\,\ln\!\left(\frac{\tan(\theta/4)}{\tan(\theta_0/4)}\right) \]
θ₀→0の極限では、どんな正の角度に到達するにも無限の時間が必要になる
ということが、この式から読み取れます。これは実は、\(\theta=0\)（てっぺん）が力学的に見て「不安定なつり合いの点」であることの帰結です。重力の接線方向成分がちょうどゼロになる特別な点なので、寸分の狂いもなく置かれた物体は、理論上は永遠にそこに留まり続けます。振り子を真上に立てたまま静止させるのと同じ状況です。現実には、ごくわずかな振動や置き方のズレが必ずあるため、この理想化はあくまで数式の上での話です。シミュレータでは、この特異点を避けるために、ごく小さな初期角 \(\theta_0\)（無限小の摂動を受けて滑り出したとみなす角度）から積分を始めています。
6. 摩擦を入れると何が起きるか — 運動方程式の書き換えとRK4
ここまでは円筒面がなめらかな場合でした。摩擦係数 \(\mu\) を持つ面の場合、接線方向の運動方程式には、動摩擦力 \(\mu N\)（運動を妨げる向き）が追加されます。
\[ mr\ddot\theta = mg\sin\theta – \mu N \]
ここに、第2節と同じ考え方で立てた半径方向のつり合い式 \(N=mg\cos\theta – mr\dot\theta^2\) を代入し、両辺を \(mr\) で割ると、
\[\ddot\theta = \frac{g}{r}\sin\theta – \mu\frac{g}{r}\cos\theta + \mu\dot\theta^2\]
という、\(\theta\) についての非線形な2階微分方程式が出てきます。ここでも質量 \(m\) はきれいに消えるので、摩擦があっても \(\theta(t)\) の形は物体の重さにはよりません。ただし、この式には第5節のような \(\int\csc u\,du\) の技が使えるきれいな閉じた解が存在しません。そこでシミュレータは、RK4（4次のルンゲ・クッタ法）という数値積分に切り替えて、ごく短い時間刻みで \(\theta\) と \(\dot\theta\) を少しずつ前に進めながら追跡します。
摩擦があると、同じ角度 \(\theta\) に到達するまでに運動エネルギーの一部が失われるので、\(v^2\)（したがって \(mv^2/r\)）は摩擦なしのときより小さくなります。第2節の式 \(N=mg\cos\theta – mv^2/r\) を思い出すと、\(v^2\) が小さいぶん \(N\) は摩擦なしのときより長く正のまま残るので、離脱角 \(\theta_S\) は \(48.19°\) よりも大きくなります。摩擦が強いほど、物体はより「粘って」から面を離れるわけです。
半径・質量への依存性はどうなるでしょうか。時刻を \(s=t\sqrt{g/r}\) という無次元の時間に置き換えると、上の運動方程式は
\[ \frac{d^2\theta}{ds^2} = \sin\theta – \mu\cos\theta + \mu\left(\frac{d\theta}{ds}\right)^2 \]
と、\(r\) も \(m\) もいっさい登場しない形に書き換えられます。つまり、\(\theta\) が時間とともにどう変化していくかという「軌跡の形」そのものは \(\mu\) だけで決まり、\(r\) は単に時間の目盛り \(\sqrt{r/g}\) を伸び縮みさせるだけです。だから、離脱角 \(\theta_S\)（摩擦の分だけ大きくなった値）は、摩擦がある場合でも半径・質量によらず一定、という性質が保たれます。
7. 摩擦だけでは動き出せないことがある — 静止摩擦と初速のトレードオフ
第5節で見たように、\(\theta=0\) のすぐ近くでは重力の接線方向成分は \(mg\sin\theta \approx mg\theta\) とごく小さい一方、摩擦がある場合はそこに最大静止摩擦力 \(\mu mg\cos\theta \approx \mu mg\) が対抗してきます。物体が動き出せるのは、重力の接線成分が最大静止摩擦を上回る、つまり
\[ \tan\theta > \mu \]
を満たすときだけです。摩擦なしの場合は、\(\theta=0\)でこそ重力の接線成分がゼロになりますが、\(\theta>0\) のどんな小さな角度でも \(\tan\theta>0=\mu\)（\(\mu=0\)）を満たすため、無限小の摂動さえあれば必ず動き出せました。ところが摩擦がある場合、話はまったく違います。\(\theta=0\) ちょうどでは重力の接線成分は厳密にゼロですが、静止摩擦は0から最大\(\mu N\)まで好きな大きさを出せるため、\(\mu\)がどれだけ小さくてもこのゼロをぴったり支え続けられます。しかも \(\theta_0\) を小さくとるほど \(\sin\theta_0\approx\theta_0\) は際限なく小さくなる一方、最大静止摩擦 \(\mu\cos\theta_0\approx\mu\) はほぼ一定のままなので、どんなに小さい \(\mu>0\) でも、それより小さい \(\theta_0\) を選べば必ず \(\sin\theta_0<\mu\cos\theta_0\) となってしまいます。第5節のような無限小の摂動を経由するトリックは、摩擦が少しでもあると原理的に機能しません。
初速0で置いた場合、摩擦係数がどれだけ小さくても、物体はθ=0ちょうどで静止摩擦につかまり、永久に動き出せない
というのが、正確な結論です。高校物理の教科書ではあまり踏み込まれない領域ですが、「なめらかな面ならどんな微小な摂動でも動き出せるのに、摩擦が少しでもあるとそうはいかない」というのは、静止摩擦と動摩擦の違いを考えるとごく自然な帰結です。シミュレータでも、摩擦係数を上げた状態で初速0のまま発射すると、物体はてっぺんから一歩も動けないまま静止し、その旨を画面に表示します。
ここで初速 \(v_0\) を与えてやると、状況が変わります。動摩擦力はあくまで「すでに動いている物体を減速させる」力であって、静止している物体を無理やり押しとどめるための力ではありません。物体がある程度の運動エネルギーを持って \(\tan\theta=\mu\) となる角度を通過してしまえば、そこから先は重力の接線成分のほうが勝つので、あとは自力で加速しながら滑り続けます。いわば、坂の途中にある摩擦の「丘」を、勢いをつけて乗り越えるイメージです。ただし、与えた初速が丘を乗り越えるには不十分だった場合、v_0>0でも結局は（てっぺんではなく）坂の途中のある角度でstuckします。逆にいえば、v_0を少しずつ大きくしていくと、ある閾値を境に急に「丘を越えられる」側に切り替わるということでもあります。
ただし、初速を導入すると、第6節で確認した「半径によらず一定」という性質は崩れます。初期の角速度は \(\omega_0=v_0/r\) なので、同じ \(v_0\) でも半径 \(r\) が変わればスタート時点の角速度そのものが変わってしまうからです。第6節の無次元化した運動方程式に立ち返ると、初期条件 \(\omega_0\) を無次元の角速度に直したときの値 \(v_0/\sqrt{gr}\) が \(r\) に依存してしまうため、離脱角 \(\theta_S\) も半径ごとに違う値になります。「半径を変えても離脱角は変わらない」という、第3節で見た性質は、実は「初速0」という条件とセットで初めて成り立つ話だったわけです。
まとめ
| 力学の操作 | 対応する結果 |
| エネルギー保存則を \(Q\) と角度 \(\theta\) の間で立てる | 速さ \(v^2 = 2gr(1-\cos\theta)\) |
| 半径方向に \(F=ma\)（向心力の式）を立てる | 垂直抗力 \(N = mg(3\cos\theta-2)\)（\(r\) を含まない） |
| \(N=0\) を解く | 離脱角 \(\cos\theta_S=2/3\)（\(\theta_S\approx48.19°\)、\(m,r\)によらず一定） |
| 位置ベクトルを \(\theta\) で微分し接線方向を求める | 離脱速度の水平・鉛直成分 \(V_x=v_S\cos\theta_S,\ V_y=v_S\sin\theta_S\) |
| \(d\theta/dt=v/r\) を変数分離して積分 | \(\theta(t)\) の解析解。\(\theta_0\to0\) で立ち上がりに無限の時間が必要という不安定平衡点の帰結 |
| 接線方向の運動方程式に動摩擦力 \(\mu N\) を追加 | 非線形2階ODEへ。解析解が失われRK4数値積分に切り替え／離脱角が48.19°より大きくなる |
| 時刻を \(s=t\sqrt{g/r}\) で無次元化 | 運動方程式から \(r,m\) が消え、離脱角の半径・質量不変性は摩擦ありでも保たれることを確認 |
| \(\tan\theta>\mu\)（動き出せる条件）を初速0で検討 | どんなに小さい摩擦でも静止摩擦につかまり物体は常に動き出せないことを発見。初速 \(v_0\) を与えると解消するが、\(\omega_0=v_0/r\) が \(r\) 依存のため半径不変性は崩れる |
この半円筒面の問題のおもしろさは、たった1つの式 \(N=mg(3\cos\theta-2)\) に \(r\) が現れないという「消える」ことの気持ちよさと、摩擦や初速を1つ加えるだけで、解析解の世界から数値積分の世界へ、そして「動き出せない」という現実的な帰結までが芋づる式に出てくるところにあると思います。教科書の公式をただ覚えるのではなく、\(F=ma\) から自分の手で導いてみると、\(\cos\theta_S=2/3\) という一見天下り的な式の背後にある構造が見えてきます。
ぜひシミュレータで半径・摩擦係数・初速を動かしながら、この記事の式と見比べてみてください。
最後まで読んでいただきありがとうございます。ご質問やコメントはお問い合わせからよろしくおねがいします。