ニュートンの運動法則 — 第1・第2・第3法則からFEM運動方程式まで

1. ニュートンの法則とCAE — なぜ今も使われるのか

アイザック・ニュートンが1687年に発表した「プリンキピア（自然哲学の数学的原理）」に記された運動3法則は、今から約340年前の成果でありながら、現代のCAE（コンピュータ支援工学）シミュレーションの数学的基盤そのものです。

古典力学の適用範囲は以下のとおりです：

• 速度が光速に対して無視できる場合（$v \ll c$）→ 特殊相対論は不要
• スケールが原子・電子より十分大きい場合（$L \gg 10^{-9}\,\text{m}$）→ 量子力学は不要
• 自動車・航空機・橋梁・建物などの構造・流体・熱解析はすべてこの範囲内

すなわち、機械・土木・航空宇宙・電機など工学分野で扱うほぼすべての問題に、ニュートン力学が直接適用できます。

2. 第1法則（慣性の法則）

「外力の合力がゼロであれば、静止物体は静止し続け、運動物体は等速直線運動を続ける。」

これが意味するのは、「物体には現在の運動状態を保とうとする性質（慣性）がある」ということです。質量 $m$ が大きいほど慣性は強く、同じ力を加えても加速しにくい。

急ブレーキ時に体が前に倒れる理由

電車が急停車するとき、乗客の体は前方に倒れます。これは体に前向きの「力が加わった」のではなく、電車が減速しても体は慣性で元の速度を保とうとするために相対的に前に出るのです。第1法則の直感的な例として最もわかりやすい現象です。

CAEにおける慣性の重要性

FEM静解析では慣性項 $[M]\{\ddot{u}\}$ をゼロとして扱います。つまり「物体は変形しても加速しない」という前提です。しかし衝突・振動・衝撃解析では慣性を無視できず、動的解析が必要になります。慣性（質量）の扱いを誤ると、解析結果が大きく外れることがあります。

3. 第2法則（運動の法則）F = ma

「物体の加速度は、加えられた合力に比例し、質量に反比例する。」

SI単位系では、$F$ の単位はニュートン（N）、$m$ はキログラム（kg）、$a$ は m/s² です。$1\,\text{N} = 1\,\text{kg}\cdot\text{m/s}^2$ という関係を常に意識しましょう。

自動車衝突の力の計算例

EURO NCAP前面フルラップ衝突試験（56 km/h バリア衝突）の場合：

• 車両質量：$m = 1{,}500\,\text{kg}$
• 衝突速度変化：$\Delta v = 50\,\text{km/h} = 13.9\,\text{m/s}$
• 衝突時間（変形による）：$\Delta t = 0.10\,\text{s}$

208 kNという力は、約21トンの重量物を持ち上げる力に相当します。この巨大な力を車体が安全に受け止め、乗員空間を守るためにクラッシュボックスや骨格構造を最適化するのが衝突CAE解析の役割です。

重力加速度と質量・重量の区別

日常会話で「60kgの人」と言うとき、これは質量です（kg）。地球上での重力による力（重力）は $W = mg = 60 \times 9.81 = 588.6\,\text{N}$ です。CAEでは重力加速度 $g = 9.80665\,\text{m/s}^2$ を体積力として入力することが多いため、単位の混同に特に注意が必要です。

4. 第3法則（作用反作用の法則）

「物体Aが物体Bに力を及ぼすとき、物体BはAに対して大きさが等しく向きが逆の力を及ぼす。」

ロケット推力の仕組み

ロケットエンジンが後方にガスを噴射すると（ガスへの作用）、その反力としてロケット機体は前方に推力を受けます（反作用）。宇宙空間では空気がないため推力が得られないと誤解する人もいますが、第3法則は2つの物体の間の力の関係であり、空気は不要です。

橋梁の支点反力

橋梁に車両が乗ると、車両は橋に下向きの力（荷重）を与えます。これに対して橋は車両に上向きの等しい力（接地反力）を返します。FEMでは境界条件（支点）で反力を計算し、支承の設計に用います。

FEM接触解析との関係

2つの構造物が接触するシミュレーション（例：歯車の噛み合い、衝突）では、第3法則を厳密に満たすことが数値安定性の鍵です。接触力の計算に使うペナルティ法やラグランジュ乗数法は、いずれも作用反作用を保証するように設計されています。

5. 力のベクトル合成と分解

力はベクトル量であり、大きさだけでなく方向も持ちます。3次元空間でのFEM解析では、各節点（ノード）に作用する力を $x, y, z$ 成分に分解して扱います。

斜面（角度 $\theta$）上に置かれた質量 $m$ の物体に働く力は：

• 斜面に沿う方向（下向き）：$F_{\parallel} = mg\sin\theta$
• 斜面に垂直な方向：$F_{\perp} = mg\cos\theta$

6. 慣性力と非慣性系（遠心力・コリオリ力）

ニュートンの法則は慣性系（加速していない座標系）で成立します。加速する座標系（非慣性系）では、見かけの力（擬似力）が現れます。

遠心力

角速度 $\omega$ で回転する座標系にいる質量 $m$ の物体には、外向きに以下の遠心力が見かけ上作用します：

ここで $r$ は回転軸からの距離です。ガスタービン翼や遠心ポンプの翼車など、高速回転する構造のCAE解析では遠心力が支配的な荷重になります。例えば航空機エンジンのタービン翼（$r=0.3\,\text{m}$, $\omega=1000\,\text{rad/s}$, $m=0.5\,\text{kg}$）では $F_c = 0.5 \times 1000^2 \times 0.3 = 150{,}000\,\text{N}$ もの遠心力が発生します。

コリオリ力

回転する座標系内で移動する物体には、コリオリ力 $\vec{F}_{Cor} = -2m\vec{\omega}\times\vec{v}$ も働きます。気象シミュレーション（台風の渦の回転方向）や、高精度な弾道計算に現れます。

7. FEM運動方程式への架け橋

ニュートンの第2法則を連続体に適用し、有限要素でアセンブルすると、以下のFEM運動方程式が得られます：

各項の物理的意味：

• $[M]\{\ddot{u}\}$：慣性力（$ma$ の行列版）— 質量マトリクス × 加速度ベクトル
• $[C]\{\dot{u}\}$：減衰力（構造減衰、材料内部摩擦）— 減衰マトリクス × 速度ベクトル
• $[K]\{u\}$：弾性復元力（$F=kx$ の行列版）— 剛性マトリクス × 変位ベクトル
• $\{F(t)\}$：外力ベクトル（荷重の時刻歴）

静解析との関係

静解析では時間変化がないため加速度も速度もゼロです：

これは運動方程式から $[M]\{\ddot{u}\}=0$, $[C]\{\dot{u}\}=0$ を除いた特殊ケースです。

固有値解析（モーダル解析）との関係

外力ゼロ・減衰ゼロの自由振動を考えると：

これが固有値問題（固有角振動数 $\omega^2$ と固有モード $\{\phi\}$ を求める問題）です。構造物の共振周波数を求めるモーダル解析はまさにこの式を解いています。

8. 実践：自動車バリア衝突シミュレーション

自動車の前面衝突解析は、ニュートンの法則を体現したCAEの代表的な応用例です。

使用するFEMソルバー

衝突解析では時間刻みが非常に小さい（～1 μs）ため、陽解法（Explicit FEM）を使います。代表的なソルバーは LS-DYNA（衝突解析のデファクト）と Abaqus/Explicit です。

解析の流れ

1. 入力パラメータ：初速（例：56 km/h）、車体質量・ジオメトリ、材料の弾塑性特性（応力-ひずみ曲線）
2. 時間積分：$\{u\}^{n+1} = \{u\}^n + \Delta t\{\dot{u}\}^n + \frac{\Delta t^2}{2}\{\ddot{u}\}^n$（中心差分法）
3. 出力：変形量・侵入量・加速度時刻歴・エネルギー吸収量

時間刻みの安定性条件

陽解法では安定性のために時間刻みに上限があります（クーラント条件）：

最小要素サイズ $L_{min} = 2\,\text{mm}$、鋼材（$E=210\,\text{GPa}$, $\rho=7850\,\text{kg/m}^3$）の場合：$c_d = 5170\,\text{m/s}$、$\Delta t_{crit} \approx 0.39\,\mu\text{s}$。衝突解析では数十万ステップを計算することになります。

9. よくある間違いと注意点