運動方程式｜力と質量から運動の法則を導く基礎

2025.09.01

運動方程式

運動方程式とは、物体や系の運動状態(位置、速度、加速度)と、それに作用する力やエネルギーとの関係を数学的に記述する基本式である。古典力学ではニュートン形式、ラグランジュ形式、ハミルトン形式が広く用いられ、剛体系、粒子系、連続体(流体・固体)にまで拡張される。工学では設計・制御・振動解析・数値シミュレーションの出発点となり、常微分方程式(ODE)や偏微分方程式(PDE)として定式化される。

Table Of Contents

定義と枠組み

系の状態を一般化座標 q(t) で表すとき、運動方程式は q, その時間微分 q̇, q̈ と外力や拘束条件の関係を与える。初期値(初期位置・初期速度)が与えられれば、時間発展を一意に求める初期値問題となる。線形・非線形、定常・非定常、保存・非保存などの性質に応じて解法や安定性が異なる。

ニュートンの運動方程式

質点系の基本は F = m a である。並進運動では質量 m と加速度 a を結び、回転運動では τ = I α (トルク τ、慣性モーメント I、角加速度 α) に対応する。減衰とばねを含む 1 自由度系は m ẍ + c ẋ + k x = f(t) と表せる。非慣性系では慣性力(遠心力、コリオリ力など)を付加して釣り合いを取る。自由物体図を描き、作用力を成分に分解して代数的に組み立てるのが実務上の手順である。

ラグランジュの運動方程式

座標系に依らない系統的な定式化として、ラグランジアン L = T − V (T:運動エネルギー, V:ポテンシャル) を用い、各座標 q_i に対し d/dt(∂L/∂q̇_i) − ∂L/∂q_i = Q_i(非保存力) を得る。関節を多く持つ機械、ロボット、リンク機構など、拘束付き多自由度系の導出に有利で、座標選択の自由度が大きい。

拘束とラグランジュ未定乗数

ホロノミック拘束(幾何学的制約)は φ(q,t)=0 として導入し、ラグランジュ未定乗数 λ を用いて拡張ラグランジアン L* = L + λ φ を構成する。非ホロノミック拘束(速度依存)はアペルの方程式やダランベールの原理に基づき一般化し、仮想仕事の枠組みで扱う。これにより反力を陽に計算せずに運動を記述できる。

ハミルトンの運動方程式

正準運動量 p_i = ∂L/∂q̇_i を導入し、ハミルトニアン H(q,p,t) に対して q̇_i = ∂H/∂p_i, ṗ_i = −∂H/∂q_i が成り立つ。エネルギー保存則や対称性(ノーターの定理)を明確に扱え、位相空間上の構造(シンプレクティック幾何)を保つ数値積分法(シンプレクティック積分)との親和性が高い。

連続体の運動方程式

密度 ρ の連続体では、質量保存と運動量保存から ρ a = ∇·σ + ρ b を得る(σ:応力テンソル, b:物体力)。構成式を与えると、固体力学では弾性体の運動方程式、流体力学ではナビエ–ストークス方程式に帰着する。波動方程式、熱弾性連成、圧電・磁気機械連成などは偏微分方程式系として解かれる。

数値解法と安定性

  • 離散化: 空間は有限要素法(FEM)、有限差分法(FDM)、有限体積法(FVM)で、時間は陽解法(オイラー陽、ルンゲ–クッタ)または陰解法(オイラー陰、Newmark-β、Crank–Nicolson)を用いる。
  • 安定条件: 陽解法では CFL 条件や固有値半径に基づく時間刻み制約が支配的である。高減衰・剛性の強い系は陰解法が有利である。
  • 制約安定化: 乗数法、ペナルティ法、Baumgarte 安定化などで拘束ドリフトを抑える。
  • 制御系: 状態空間表現 ẋ = A x + B u, y = C x + D u として設計・同定に接続する。

工学的応用

機械設計では多自由度振動系を M ẍ + C ẋ + K x = f(t) の行列表現にまとめ、固有値解析で共振やモード形状を把握する。基礎励振や地震応答、回転体の不釣り合い、ギヤかみあい励振などは周波数応答や時刻歴解析で評価する。メカトロニクスではモータ・センサ・制御器を含む物理–情報の統合モデルに運動方程式を埋め込み、観測ノイズや非線形性を含む現実的な挙動を再現する。

非線形・大変位・接触

剛性の変位依存(幾何学的非線形)、材料非線形(塑性・粘弾性)、摩擦・バックラッシュ・衝突などの不連続性は多様な現象を生む。ダフフィング方程式に代表される硬化・軟化ばね、内部共振、カオスは設計余裕度や安全率の設定に影響する。接触問題では補題条件(Complementarity)やペナルティ/正規反力–摩擦則を数値的に解く。

代表的なモデル化の流れ

  1. 境界・支持条件、荷重条件、座標系の定義(自由度の選定)。
  2. 運動エネルギー・位置エネルギーの導出、非保存力(減衰・駆動)のモデル化。
  3. 運動方程式の導出(ニュートン/ラグランジュ/ハミルトン)。
  4. 無次元化とパラメータ同定、線形化の可否判断。
  5. 数値解法の選択、検証・妥当性確認(実験・モーダル解析)。

簡単な例

単振動系: m ẍ + c ẋ + k x = 0。減衰比 ζ = c/(2√(mk))、固有角振動数 ω_n = √(k/m) を定義すると応答の性質(過減衰・臨界・不足減衰)が決まる。剛体回転では I θ̈ + c_r θ̇ + k_r θ = τ(t) としてトルク入力に対する角変位応答を解析する。投射運動に空気抵抗を加えれば m v̇ = m g − D v − C v|v| の非線形 ODE となり、数値解が一般的となる。

コメント（β版）