運動量と力積 — 衝突解析・エアバッグ設計からExplicit FEMまで
関連トピック
1. 運動量とは — 「勢い」の物理的定義
エアバッグって膨らんで人を守るじゃないですか、あれって力積と関係あるんですか？
まさにそれが力積の教科書的な応用例だ。$J = F\Delta t = \Delta p$ という式があって、同じ運動量変化 $\Delta p$ でも衝突時間 $\Delta t$ を長くすれば力 $F$ を小さくできる。エアバッグは「人体がハンドルに直接当たる」のではなく「クッションに沈み込む時間を引き延ばす」デバイスなんだ。その分、人体に加わる最大加速度が下がって傷害が軽減される。
運動量（momentum）は、物体の「動きの勢い」を定量化した量です：
単位はkg·m/s（= N·s）。ニュートンの第2法則を運動量で表現すると：
これが実はニュートンが元々書いた形式に近く、$F=ma$ はこの特殊ケース（質量一定）です。ロケット推進のように質量が変化する場合は $F = \frac{d(mv)}{dt} = m\dot{v} + \dot{m}v$ を使います。
2. 力積と運動量の変化
時刻 $t_1$ から $t_2$ の間に物体に力 $\vec{F}(t)$ が作用したとき、力積（impulse）は：
力が一定の場合：$\vec{J} = \vec{F}\Delta t$。力積-時間グラフの面積が運動量変化に等しいという関係は、衝撃試験データの解析に直接使えます。
エアバッグ展開時の力積計算例
60kgの人が 50km/h（13.9m/s）から 0 に減速する場合、運動量変化は：
エアバッグなしでハンドルに直撃（$\Delta t = 0.005\,\text{s}$）：$F = 834/0.005 = 166{,}800\,\text{N}$（約17トン相当）
エアバッグあり（$\Delta t = 0.030\,\text{s}$）：$F = 834/0.030 = 27{,}800\,\text{N}$（約2.8トン相当）
時間を6倍に伸ばすことで力を1/6に減らせます。これがエアバッグの物理的原理です。
3. 運動量保存の法則
系全体に外力が働かない（または無視できる）場合、系の全運動量は保存されます：
2体の衝突では：
ロケット推進
質量 $M$ のロケットが速度 $U$ で燃料ガスを噴射すると、運動量保存からロケット機体が加速します。ツィオルコフスキーのロケット方程式：
ここで $v_e$ は排気速度、$M_0$ は初期質量、$M_f$ は最終質量（燃料を使い切った後）。これも運動量保存から導けます。
4. 弾性衝突（はね返り係数 e = 1）
弾性衝突では運動量保存に加えて運動エネルギーも保存されます。2つの連立方程式を解くと：
特別なケース：$m_1 = m_2$（等質量）のとき $v_1' = v_2$, $v_2' = v_1$（完全速度交換）。ビリヤードの玉が止まって的の玉が動き出すのはこのためです。
5. 非弾性衝突と完全非弾性衝突（e = 0）
はね返り係数（coefficient of restitution）：
$e=1$：完全弾性衝突（KE保存）、$e=0$：完全非弾性衝突（一体化）、$0 < e < 1$：非弾性衝突
完全非弾性衝突（$e=0$）では衝突後に一体になって動きます：
失われた運動エネルギー（変形・熱・音に変換）：
自動車事故の解析でははね返り係数e=0を使うんですか？
実際は0.1〜0.3くらいの値になる。完全非弾性（e=0）ではない。でも最近の高精度な衝突解析ではExplicit FEMを使うから、eは「入力値」ではなく「計算結果」として得られるんだ。材料の弾塑性特性・接触摩擦・構造剛性を正確にモデル化すれば、eに相当する挙動はFEMが自動的に再現してくれる。
6. 爆発と分裂（逆衝突）
静止していた物体が爆発によって分裂する場合も運動量保存が成立します。初期運動量がゼロなら分裂後の運動量の和もゼロ：
砲弾発射の後座力計算
質量 $m=10\,\text{kg}$ の砲弾が $v=800\,\text{m/s}$ で発射される場合、大砲（$M=2000\,\text{kg}$）の後座速度：
後座力制御メカニズムの設計には、この後座速度とブレーキ機構の力積計算が必要です。
7. 2次元衝突と斜め衝突
2次元では $x$ 方向と $y$ 方向それぞれで運動量が保存されます：
交差点でのT字衝突（自動車事故解析）
$m_A = 1500\,\text{kg}$（速度 20m/s、東方向）と $m_B = 1200\,\text{kg}$（速度 15m/s、北方向）が完全非弾性衝突した場合の合成速度方向：
事故調査では衝突後のスリップ痕の方向からこのベクトル計算を逆算して衝突前速度を推定します。
8. Explicit FEMへの架け橋（衝突シミュレーション）
FEMの陽解法では、運動量の時間更新が時間積分の基本操作です：
ここで $\vec{F}^n_{int}$ は内力（隣接要素の応力から計算）、$\vec{F}^n_{ext}$ は外力（重力・接触力など）です。
接触アルゴリズム
- ペナルティ法：接触面の侵入量に比例した接触力を加算。実装が簡単で最も広く使われる。ばね剛性の設定が結果に影響する。
- ラグランジュ乗数法：非侵入条件を厳密に満たす。精度が高いが方程式が大きくなる。
- Augmented Lagrangian法：両者の折衷。
LS-DYNAでの衝突解析ワークフロー
- メッシュ作成：最小要素サイズが時間刻みを決定
- 材料定義：弾塑性（MAT24）、座屈（MAT36）など
- 接触設定：
*CONTACT_AUTOMATIC_SINGLE_SURFACEが汎用 - 荷重・拘束：初速度
*INITIAL_VELOCITYで設定 - 実行・後処理：エネルギーバランス確認が最優先
9. 実践：歩行者保護基準（EURO NCAP）
EURO NCAPの歩行者保護試験では、歩行者の頭部が車のボンネット（フード）に衝突した際の傷害リスクを評価します。
頭部傷害基準 HIC（Head Injury Criterion）
ここで $a(t)$ は頭部重心の合成加速度（単位：g）。積分区間 $[t_1, t_2]$ は最大15ms（15mmHg以内で最大値になるよう選択）。
- $HIC_{15} < 700$：良（緑）
- $700 \leq HIC_{15} < 1000$：許容（橙）
- $HIC_{15} \geq 1000$：危険（赤）
FEM解析との相関
頭部インパクターをボンネット各点に衝突させるFEM解析を行い、HIC分布マップを作成します。この解析で弱点となる領域（HICが高い部位）を特定し、フードヒンジ形状の変更・アクティブフード（衝突時に跳ね上がる機構）の採用などの対策に繋げます。実際の試験との相関係数 $r > 0.9$ を目標に、モデル更新を繰り返します。