運動学（キネマティクス）— 位置・速度・加速度とFEM時間積分

カテゴリ: 物理の基礎 | 2026-03-25 | サイトマップ
NovaSolver Contributors

関連トピック

基礎物理ニュートンの運動法則基礎物理仕事・エネルギー・動力基礎物理応力とひずみの基礎基礎物理単振動と固有振動数構造解析構造解析へ進むシミュレーターインタラクティブツール
CAE visualization for kinematics - technical simulation diagram
Kinematics

1. 運動学とは — 「なぜ動くか」より「どう動くか」

🧑‍🎓

速度と加速度の違いって、わかるようでわからないんですが… ざっくり説明してもらえますか？

🎓

速度は位置の変化率、加速度は速度の変化率だ。数学で言うと微分の関係。「どれだけ速く動いているか」が速度で、「速さが変わる早さ」が加速度。これが分かるとFEMの時間積分アルゴリズム — Newmark-β法や中心差分法 — の数式がスッと読めるようになる。

運動学（Kinematics）は、「なぜ物体が動くか（力の原因）」を問わず、「物体がどのように動くか（位置・速度・加速度の記述）」だけを扱う分野です。一方、力と運動を結びつけるのは動力学（Dynamics）です。

FEM動解析における時間積分は、まさに運動学の離散化です。連続的な運動方程式を「時間を細かく刻んで一歩ずつ計算する」操作が時間積分であり、その精度と安定性は運動学の理解に直結します。

2. 位置・速度・加速度の微積分関係

1次元の運動を例にとると：

$$v(t) = \frac{dx}{dt}, \quad a(t) = \frac{dv}{dt} = \frac{d^2x}{dt^2}$$

積分方向（逆向き）には：

$$v(t) = v_0 + \int_0^t a(\tau)\,d\tau, \quad x(t) = x_0 + \int_0^t v(\tau)\,d\tau$$

この微積分の連鎖が、FEMの時間積分アルゴリズムの本質です。時間積分とは「加速度 $a$ から速度 $v$ を、速度 $v$ から変位 $x$ を数値的に積分する」操作に他なりません。

加速度計センサーの原理

スマートフォンや自動車に搭載されているMEMS加速度計は $a(t)$ を直接測定します。測定した加速度を時間積分すると速度 $v(t)$ が、さらに積分すると変位 $x(t)$ が得られます。自動車衝突試験ではダミー（人形）の胸部加速度を計測し、加速度波形からHIC（頭部傷害基準）などの傷害指標を計算します。

3. 等加速度運動の3公式

加速度が一定（$a = \text{const}$）の場合、積分を解析的に実行できます：

$$v = v_0 + at \tag{1}$$ $$x = v_0 t + \frac{1}{2}at^2 \tag{2}$$ $$v^2 = v_0^2 + 2ax \tag{3}$$

これらはすべて上の微分関係式から導けます。（3）式は（1）と（2）から $t$ を消去したものです。

自動車の急制動距離計算

初速 $v_0 = 100\,\text{km/h} = 27.8\,\text{m/s}$、急制動時の減速度 $a = -8\,\text{m/s}^2$、最終速度 $v = 0$ として（3）式を使うと：

$$0 = v_0^2 + 2ax \implies x = -\frac{v_0^2}{2a} = -\frac{(27.8)^2}{2\times(-8)} = 48.3\,\text{m}$$

約48mの制動距離が必要です。これが高速道路での車間距離確保の根拠の一つです。実際の自動車の制動距離計算では路面μ（タイヤ-路面摩擦係数）を考慮し、乾燥路では $a \approx -8\,\text{m/s}^2$、濡れた路面では $a \approx -4\,\text{m/s}^2$ 程度を使います。

4. 2次元運動とベクトル分解（放物運動）

水平方向（$x$）と鉛直方向（$y$）を独立に扱うことで、2次元運動を解けます：

$$x(t) = v_0\cos\theta \cdot t$$ $$y(t) = v_0\sin\theta \cdot t - \frac{1}{2}gt^2$$

最大飛距離は $y = 0$ となる時刻 $t = \frac{2v_0\sin\theta}{g}$ を $x(t)$ に代入すると：

$$R = \frac{v_0^2 \sin 2\theta}{g}$$

$\sin 2\theta$ が最大になる $\theta = 45°$ で最大射程になることがわかります（空気抵抗無視）。

🧑‍🎓

野球のホームランってどれくらいの速さで打てば外野フェンスを越えますか？ 計算できます？

🎓

放物運動で計算できるよ。打ち出し速度45 m/s、角度35°で飛距離を計算してみると $R = (45)^2\sin(70°)/9.81 \approx 194\,\text{m}$（空気抵抗無視）。実際は空気抵抗で130m前後になる。最近は打球速度センサーで実測してFEM衝撃解析と組み合わせる研究もあるよ。バット-ボール衝突の応力分布を計算したりね。

5. 角運動学（回転運動の記述）

回転運動も同じ微積分の枠組みで記述できます：

$$\omega = \frac{d\theta}{dt}, \quad \alpha = \frac{d\omega}{dt} = \frac{d^2\theta}{dt^2}$$

線速度との関係：$v = r\omega$、接線加速度：$a_t = r\alpha$、向心加速度：$a_c = r\omega^2 = v^2/r$。

タービン翼の相対運動

ガスタービンエンジンでは、回転する動翼（ロータ）と静止した静翼（ステータ）の間の相対速度を正確に計算することが重要です。絶対座標系と回転座標系の間の速度の合成（ベクトル三角形）が翼の設計の基礎です。回転機械のFEM解析では、回転体の角速度と角加速度を正確に入力することが計算精度の前提となります。

6. 相対運動と速度の合成

物体Aの速度を、物体Bに対して表すには：

$$\vec{v}_{A/B} = \vec{v}_A - \vec{v}_B$$

例えば川を渡るボートの問題（ボートの対水速度 + 川の流れ速度 = 対地速度）は、この速度合成の典型例です。

マルチボディ動力学（MBD）への接続

複数の剛体がジョイントで連結されたシステム（自動車サスペンション・ロボットアーム）を解析するマルチボディ動力学（MBD）では、各リンク間の相対運動を追跡することが基本です。Adams（MSC Software）などのMBDソルバーは、この相対運動の拘束方程式をリアルタイムに解いています。

7. FEM時間積分法への架け橋

連続的な運動を離散的な時間ステップで解くための数値手法が時間積分法です。

陽解法 — 中心差分法（Explicit）

現在時刻の情報のみで次の時刻の変位を計算します：

$$\ddot{u}_n \approx \frac{u_{n+1} - 2u_n + u_{n-1}}{\Delta t^2}$$ $$u_{n+1} = 2u_n - u_{n-1} + \ddot{u}_n \Delta t^2$$

長所：各ステップの計算が軽い（大規模連立方程式を解かない）。衝突・爆発など短時間の過渡現象に強い。
短所：安定性条件（クーラント条件）を満たす必要があり、時間刻みを細かくしなければならない。

陰解法 — Newmark-β法（Implicit）

次の時刻の情報も使って現在時刻を解く（連立方程式を解く）手法：

$$u_{n+1} = u_n + \dot{u}_n\Delta t + \left(\frac{1}{2}-\beta\right)\ddot{u}_n \Delta t^2 + \beta\ddot{u}_{n+1}\Delta t^2$$ $$\dot{u}_{n+1} = \dot{u}_n + (1-\gamma)\ddot{u}_n\Delta t + \gamma\ddot{u}_{n+1}\Delta t^2$$

$\beta = 0.25$, $\gamma = 0.5$（平均加速度法）は無条件安定です。$\beta = 0$（$\gamma = 0.5$）は中心差分法に対応します。

🧑‍🎓

FEMで動解析するとき、時間刻みはどう決めるんですか？ 適当に選んじゃダメですか？

🎓

陽解法では絶対ダメだ。安定性条件が厳しくて、クーラント数 $CFL = c\Delta t/\Delta x \leq 1$ を満たさないと解が発散する。衝突解析で要素サイズが2 mmなら、鋼の音速 $c \approx 5170\,\text{m/s}$ から $\Delta t_{max} \approx 0.39\,\mu\text{s}$。Abaqus Explicitは最小要素サイズから自動計算してくれる。陰解法は無条件安定だけど、時間刻みが荒すぎると精度が落ちるから結果を確認しながら刻み幅を決めよう。

安定性条件（クーラント条件）

$$CFL = \frac{c\,\Delta t}{\Delta x} \leq 1, \quad c = \sqrt{\frac{E}{\rho}}$$

最小要素サイズ $\Delta x = 2\,\text{mm}$, 鋼（$E=210\,\text{GPa}$, $\rho=7850\,\text{kg/m}^3$）の場合：$c = 5170\,\text{m/s}$ → $\Delta t_{max} = 0.387\,\mu\text{s}$。衝突時間が100 ms なら 258,000 ステップの計算が必要です。

8. 計測とシミュレーションの接続（V&V）

CAEシミュレーションの信頼性を検証（V&V: Verification and Validation）するためには、実験計測との比較が不可欠です。

加速度センサー（MEMSと圧電型）

衝突試験データとのFEM比較

衝突試験では通常50〜100点の加速度センサーを車体とダミーに取り付け、加速度波形を記録します。FEM解析でも同じ位置の節点加速度を出力し、波形の一致を確認します。一致が悪い場合は材料モデルの改善・接触条件の見直し・メッシュの細分化などの対策を講じます。

9. 練習問題と解答

問題1：最高点の計算

初速 $v_0 = 20\,\text{m/s}$ で真上に打ち上げた物体の最高点の高さを求めよ。($g = 9.8\,\text{m/s}^2$)

解答を表示

最高点では $v=0$。公式 $v^2 = v_0^2 - 2gh$ より $h = \frac{v_0^2}{2g} = \frac{400}{19.6} \approx 20.4\,\text{m}$

問題2：制動距離計算

初速 $v_0 = 72\,\text{km/h}$、制動加速度 $a = -6\,\text{m/s}^2$ での制動距離を求めよ。

解答を表示

$v_0 = 72/3.6 = 20\,\text{m/s}$。公式 $v^2 = v_0^2 + 2ax$ より $x = \frac{-v_0^2}{2a} = \frac{400}{12} \approx 33.3\,\text{m}$

問題3：Newmark-β法の手計算（1ステップ）

$m=1\,\text{kg}$, $k=100\,\text{N/m}$, $F=10\,\text{N}$（定常）、$\Delta t = 0.01\,\text{s}$、初期条件 $u_0=0$, $\dot{u}_0=0$。Newmark-β法（$\beta=0.25$, $\gamma=0.5$）で1ステップ後の変位を求めよ。

解答を表示

$\ddot{u}_0 = (F - ku_0)/m = 10\,\text{m/s}^2$。予測値：$u_1^* = 0 + 0 + (0.5-0.25)\times10\times0.0001 = 0.00025$。剛性修正：$\hat{k} = k + m/(\beta\Delta t^2) = 100 + 40000 = 40100$。$\hat{F} = 10 + 40000\times0 + \cdots$（詳細な手順はNewmark-β記事参照）

関連ページ・ツール

分野横断タグ（Cross-topics）

#数値解析   #動解析   #ニュートン力学   #構造解析   #単振動