1次元・2次元の弾性・非弾性衝突をリアルタイムアニメーション。反発係数を変えて運動量保存・エネルギー損失を体験。衝撃力-時間グラフと力積を可視化。
反発係数: $e = \frac{v_2'-v_1'}{v_1-v_2}$
衝突現象の核心は「運動量保存の法則」です。外力が働かない場合、衝突前後の運動量の合計は保存されます。
$$m_1 v_1 + m_2 v_2 = m_1 v_1' + m_2 v_2'$$$m_1, m_2$: 物体1,2の質量 [kg], $v_1, v_2$: 衝突前の速度 [m/s], $v_1', v_2'$: 衝突後の速度 [m/s]
衝突の「跳ね返り具合」を定量的に表すのが「反発係数」です。衝突前後の相対速度の比で定義されます。
$$e = -\frac{v_2' - v_1'}{v_2 - v_1}= \frac{\text{離散相対速度}}{\text{接近相対速度}}$$$e$: 反発係数 (0 ≤ e ≤ 1)。e=1で弾性衝突（運動エネルギー保存）、e=0で完全非弾性衝突（合体）。
衝突時に働く力と時間の関係を表すのが「力積-運動量の定理」です。力積は運動量の変化に等しく、力-時間グラフの面積として可視化されます。
$$J = \int F \, dt = \Delta p = m \Delta v$$$J$: 力積 [N·s], $F$: 力 [N], $\Delta t$: 作用時間 [s], $\Delta p$: 運動量変化 [kg·m/s]
自動車衝突安全解析：CAEソフト（LS-DYNA等）を用いて車体のクラッシュシミュレーションを行います。フロント部分のつぶれ（クラッシャブルゾーン）により衝突時間を意図的に延ばし、乗員室への衝撃力のピークを低減します。エアバッグの展開タイミングもこの衝撃力波形に基づいて設計されます。
スポーツ用具・保護具設計：野球のグローブやヘルメットの内部緩衝材は、反発係数を低くし衝突時間を長くすることで、手や頭部への衝撃力を分散させます。同じ原理で、競技用フロアやマットも衝撃吸収性能が評価・設計されています。
ロボット・精密機械の動作設計：産業用ロボットアームがワークを把持・設置する際の衝撃や、精密機器内の可動部の接触を制御するために衝突解析が用いられます。非弾性衝突モデルを用いて、振動や位置ずれを最小化する動作条件をシミュレーションで事前に検討します。
ボールスポーツの解析：ゴルフクラブとボール、テニスラケットとボールなど、スポーツにおける衝突の研究に応用されます。反発係数（COR）は用具の規格にも用いられ、シミュレーションで衝突後のボールの飛距離や方向を予測するのに役立ちます。
このシミュレーターを使い始めるときに、特に気をつけてほしいポイントがいくつかあるよ。まず一つ目は、「反発係数eは材質だけで決まらない」ということ。よく「鋼はe=0.9、粘土はe=0」みたいに覚えがちだけど、実は衝突速度や物体の形状でも大きく変わるんだ。例えば、同じ鋼球でも、超高速で衝突すると局部が塑性変形してeは低下する。シミュレーターでは単一のパラメータとして設定しているけど、実務のCAEではもっと複雑なモデルを使うことを頭に入れておこう。
二つ目は、「運動量保存」と「力のつり合い」を混同しないこと。衝突の瞬間、二つの物体には作用・反作用の法則で必ず同じ大きさの力が働く。でも、これは「力」が等しいのであって、「運動量の変化」は質量が違えば異なる。重いボール（m1=10kg）と軽いボール（m2=1kg）が衝突すると、軽い方がはるかに大きく速度を変えるよね。運動量変化（力積）は両者で等しいけど、速度変化は質量に反比例するんだ。
最後に、シミュレーターの設定で陥りがちなのが、初期条件の現実性。例えば、質量を極端に大きく（1000kg）、速度も高速（100m/s）にすると、計算上は成立しても、現実ではとんでもないエネルギーが発生して物体は粉々だ。あくまで原理の理解が目的だから、ボールの質量は1〜10kg、速度は数m/s程度の現実的な範囲で遊ぶのが、現象を正しく理解するコツだよ。
この1D/2D衝突の原理は、思っている以上に幅広い工学分野の基礎になっているんだ。まず挙げるのは「自動車の衝突安全（クラッシュワースィネス）」。シミュレーターで見た力積と力のピークの関係は、まさに車体設計の核心だ。フロント部分を段階的につぶすことで衝突時間を数十ミリ秒から数百ミリ秒に延ばし、乗員への衝撃を緩和する。CAEツールではLS-DYNAやRADIOSSなどがこの多次元衝突を詳細にシミュレートしている。
次に「ロケットの多段分離」も衝突・運動量の応用だ。ロケットの一段目と二段目を分離する時、小型の爆発ボルトやスプリングで互いに押し合う。これは内部的な力だから運動量は保存され、軽い上段が大きく加速し、重い下段はゆっくり離れていく。シミュレーターで質量比を極端に変えて衝突させると、この現象がイメージしやすいよ。
もう一つ、「粉体工学（パウダーテクノロジー）」も見逃せない。製薬工場で薬の粒子を混合したり、セメントの製造で微粒子が衝突したりする過程を解析するのに、粒子法（DEM: Discrete Element Method）というシミュレーション手法が使われる。ここでは無数の粒子一つ一つの衝突（反発係数、摩擦）を計算して、マクロな粉の流れや混合度を予測するんだ。我々のシミュレーターは、そのごくごく初歩的な「一衝突」の物理を体感させてくれる。
このツールで基本に慣れたら、次はステップアップしてみよう。まずおすすめは、「2次元衝突（斜衝突）に立ち向かう」こと。ツールを2Dモードに切り替えて、ボールを斜めに衝突させてみて。その時、運動量保存則はx成分とy成分でそれぞれ独立に成り立つ。反発係数の式も、衝突面に垂直な速度成分だけに適用されるんだ。例えば、なめらかな床にボールを斜めにぶつける場合、水平成分の速度は変わらず、垂直成分だけが反発係数で変化する。これを理解できれば、ビリヤードの球の動きも原理的に説明できるようになる。
数学的には、衝突前後の速度を求めるのは、運動量保存則と反発係数の式の連立方程式を解く作業だ。1次元なら簡単だけど、2次元ではベクトル計算が必要になる。ここで強力な武器になるのが「はね返り後の相対速度ベクトル」の考え方だ。衝突面の法線ベクトルを考え、それに沿った成分だけを反発係数で処理する。この辺りを学ぶには、剛体の平面運動の教科書が役に立つよ。
さらに先へ進みたければ、「弾性衝突における運動エネルギー保存の証明」に挑戦してみよう。反発係数e=1の時、衝突前後の運動エネルギーの和が保存されることを、運動量保存則$$m_1 v_1 + m_2 v_2 = m_1 v_1' + m_2 v_2'$$と反発係数の式$$v_2' - v_1' = - (v_2 - v_1)$$から導出できる。この計算を一度自分の手で追ってみると、エネルギーと運動量の関係が深く理解できる。それができたら、次のトピックは「角運動量を考慮した3次元剛体衝突」や「塑性変形を伴う衝突のCAE」へと広がっていくんだ。