剛体力学の本質的理解: 基礎概念から人体歩行モデルまでの段階的構築

Yoshinori Nomura

Published: Dec. 2, 2025, 5:38 a.m. UTC / Updated: Dec. 2, 2025, 6:49 a.m. UTC

👍 0

🔖 0 Bookmarks

日本語

Executive Summary

目的と到達目標

本レポートは、剛体力学の初心者が人体歩行のような複雑な多体システムの運動方程式を記述できるようになるまでの、体系的な学習カリキュラムを提示する。制御理論を除外し、運動の数理的記述能力の獲得に焦点を絞った教科書設計である。

カリキュラムの構造と教育戦略

本カリキュラムは9つの段階的ステップで構成され、「概念の深さ・正確性最優先」の方針に基づく。各ステップでは、問題意識の明確化、理論的解決方法の提示、鍵概念の懇切丁寧な解説を行う。

学習の進行構造

基礎段階（ステップ1-3） : 座標系・参照枠の概念から始まり、質点の運動記述、ニュートン力学の基本原理を確立する。観察者依存性や座標系選択の本質的意味を重視。
回転の導入（ステップ4-5） : 並進運動では捉えきれない回転運動の複雑性を認識し、角速度ベクトル、慣性テンソル、オイラー方程式を導入。変形可能な物体から剛体への抽象化を、物理的必然性と数学的帰結の両面から解説。
多体系への拡張（ステップ6-8） : 内力と外力の区別、接触力のモデル化、拘束条件の分類（ホロノミック/非ホロノミック）を経て、ラグランジュ形式による多体システムの一般的記述法を習得。
実践的検証（ステップ9） : 人体歩行モデルを具体例として、理論の統合的適用を実証。モデル化の意思決定プロセスと計算実行可能性を検証。

主要な教育的特徴

概念理解の徹底

剛体の定義 : 「任意の二点間距離が時間不変」という数学的定義から、自由度削減の帰結までを厳密に導出
座標系の本質 : 単なる計算道具ではなく、物理法則の記述基盤としての座標系の役割を明確化
拘束条件の意味 : 自由度削減の数学的意味と、拘束力が未知変数として現れるメカニズムを解明

段階的複雑性の管理

各ステップは前段階の知識のみを前提とし、論理的飛躍を排除。質点→単一剛体→二体系→多体系という自然な流れで、学習者の概念的負荷を制御。

実用性の保証

最終ステップで人体歩行の3剛体モデル（胴体-大腿-下腿）の運動方程式導出を完遂し、理論の適用可能性を具体的に実証。

到達される能力

本カリキュラムを完了した学習者は以下を獲得する：

モデル化能力 : 複雑な実在物体を適切な剛体系として抽象化し、座標系・自由度・拘束条件を定義できる
方程式導出能力 : Newton-Euler法またはLagrange法により、多体システムの完全な運動方程式を記述できる
物理的洞察 : 数式の背後にある物理的意味を理解し、結果の妥当性を判断できる

強みと限界

主要な強み

数学的厳密性 : 初心者向けでありながら、概念定義の曖昧さを排除
体系的構成 : 各概念の必然性が前段階から自然に導かれる論理構造
実証的検証 : 抽象理論が人体モデルという具体例で機能することを確認

認識すべき限界

学習負荷 : 概念の厳密な理解を要求するため、全体のボリュームは相当量となる
計算複雑性 : 高自由度システムでは、方程式の複雑さが実用的障壁となりうる
制御との接続 : 運動方程式の記述で終わり、制御設計への橋渡しは含まない

推奨される使用法

自習者向け

各ステップの「問題意識」セクションを熟読し、なぜその概念が必要かを理解してから数式に進むこと。MathJax数式と文章説明を相互補完的に活用すること。

教育者向け

ステップ4（回転）とステップ7（拘束）で学習者の負荷が増大するため、演習問題や具体例の追加を推奨。計算ツール（MATLAB, Python）との併用が効果的。

カリキュラム設計者向け

本構成は1学期（15週）での実施を想定。時間制約がある場合は、ステップ8の一般化理論を簡略化し、ステップ9の具体例に時間を配分する選択肢もある。

結論

本カリキュラムは、概念の正確な理解を基盤とした長期的応用可能性を重視する設計である。人体歩行を含む複雑な多体システムの運動記述という目標に対し、必要十分な理論的基盤と実践的検証の両方を提供する。初学者が段階的に理解を積み上げ、最終的に独力で複雑構造の運動方程式を導出できるようになる、実現可能な学習パスを構成している。

基礸記述能力の前提条件

この章の問題意識

物体の運動を記述するという行為は、私たちが日常的に観察する現象を数学的な言語に翻訳する作業である。しかし、この翻訳作業を始める前に、私たちは「何を記述するのか」「どのように記述するのか」という根本的な問いに答える必要がある。

多くの初学者が剛体力学の学習で挫折する原因は、 記述の対象と記述の方法についての共通理解が曖昧なまま、複雑な数式の操作に入ってしまう ことにある。例えば、「物体の位置」という一見単純な概念でさえ、実は多くの暗黙の前提を含んでいる。位置とは何に対しての位置なのか? 物体のどの部分の位置なのか? 時間とともに変化する位置をどう捉えるのか?

本章では、剛体力学という建物を建てる前に、その土台となる概念的な基礎を確立する。ここで扱う内容は一見当たり前に思えるかもしれないが、これらの「当たり前」を明確に言語化し、理解することが、後の複雑な理論を正確に理解するための不可欠な準備となる。

観察者と座標系：記述の出発点

なぜ座標系が必要なのか

物体の運動を記述する際、最初に直面する根本的な問いは「 誰が見ているのか 」である。この問いは哲学的に聞こえるかもしれないが、物理学における最も実践的な問題の一つである。

日常生活で「ボールが動いている」と言うとき、私たちは暗黙のうちに「地面に対して」動いていることを想定している。しかし、電車の中でボールを投げた場合、電車内の観察者と駅のホームにいる観察者では、ボールの運動の記述が全く異なる。両者とも正しい観察をしているが、異なる 参照枠 を使っている。

この「参照枠」を数学的に具体化したものが 座標系 である。座標系とは、空間内の任意の点に数の組(通常は3つの実数)を対応させる体系的な方法である。

座標系の構成要素

座標系を完全に定義するには、以下の3つの要素が必要である:

原点 (Origin) : 位置を測る基準となる点
基底ベクトル (Basis Vectors) : 方向を定義する独立なベクトルの組
スケール (Scale) : 距離を測る単位

最も標準的な座標系は デカルト座標系 (直交座標系)である。原点を $$O$$ とし、互いに直交する単位ベクトル $\mathbf{e}_x, \mathbf{e}_y, \mathbf{e}_z$ を選ぶと、空間内の任意の点 $$P$$ の位置は 位置ベクトル として表現できる:

\mathbf{r} = x\mathbf{e}_x + y\mathbf{e}_y + z\mathbf{e}_z

ここで $$(x, y, z)$$ は点 $$P$$ の座標と呼ばれる実数の組である。

座標系の選択の自由と責任

重要な認識は、 座標系の選択は任意である という点だ。同じ物理現象を記述するのに、無限に多くの座標系を選ぶことができる。しかし、この自由には責任が伴う:

一貫性 : 一つの問題を解く間、座標系を勝手に変えてはならない(または変える際は明示的に変換を行う必要がある)
適切性 : 問題の対称性や境界条件に合った座標系を選ぶことで、計算が大幅に簡略化される
明示性 : どの座標系を使っているかを常に明確にしなければならない

例えば、円運動を記述する際、デカルト座標系よりも 極座標系 や 円筒座標系 を用いる方が自然である。人体の歩行を記述する際には、地面に固定された座標系と、身体の特定部位(例えば骨盤)に固定された座標系を使い分けることが有効である。

時間パラメータ：変化の記述

時間という特別なパラメータ

運動とは「 時間とともに位置が変化すること 」である。この定義は単純に見えるが、時間を数学的にどう扱うかは慎重に考える必要がある。

物理学における時間は、実数直線 $\mathbb{R}$ 上のパラメータ $$t$$ として扱われる。つまり、時間の各瞬間は実数に対応し、その順序関係と連続性が保たれる。この数学的抽象化により、微分積分学の道具を運動の記述に適用できるようになる。

時刻と時間間隔の区別

初学者がしばしば混同するのが「時刻」と「 時間間隔 」である:

時刻 (Time instant) : 時間軸上の特定の点。記号 $$t$$ で表されることが多い。「午後3時」や「実験開始から5秒後」などが該当する。
時間間隔 (Time interval) : 2つの時刻の間の長さ。記号 $\Delta t$ や $\tau$ で表される。「3秒間」などが該当する。

数学的には、時刻 $$t_1$$ と $$t_2$$ に対して、時間間隔は $\Delta t = t_2 - t_1$ である。

時間の原点の任意性

空間座標と同様に、 時間の原点(いつを $$t=0$$ とするか)も任意に選べる。これは実用的に非常に重要である。例えば:

ボールを投げる問題では、投げた瞬間を $$t=0$$ とするのが自然
周期運動では、特定の基準位置を通過する瞬間を $$t=0$$ とすることが多い
人体歩行の解析では、片足が地面に接地した瞬間を $$t=0$$ とすることがある

時間原点の選択は、初期条件の設定と密接に関連する。

位置・速度・加速度：運動状態の階層構造

位置：運動記述の第一段階

物体の位置は、選択した座標系において、物体がどこにあるかを示す。数学的には、位置は時間の関数として表現される:

\mathbf{r}(t) = x(t)\mathbf{e}_x + y(t)\mathbf{e}_y + z(t)\mathbf{e}_z

ここで重要なのは、位置ベクトル $\mathbf{r}(t)$ が時間 $$t$$ の関数であるという点だ。これは 軌道 (Trajectory) を定義する。

位置の記述には、暗黙の前提がある: 我々は物体を「点」として扱っている。これは厳密には近似であり、後に剛体を扱う際には、物体が空間的な広がりを持つことを考慮する必要がある。しかし、まずは点として扱える物体(または物体の特定の点)を考えることから始める。

速度：変化の割合

物体が運動しているとき、私たちは「どれくらい速く動いているか」を知りたい。この問いに答えるのが速度の概念である。

直感的には、速度は「単位時間あたりの位置の変化」である。微小時間間隔 $\Delta t$ の間に位置が $\mathbf{r}(t)$ から $\mathbf{r}(t + \Delta t)$ に変化したとき、平均的な速度は:

\mathbf{v}_{\text{avg}} = \frac{\mathbf{r}(t + \Delta t) - \mathbf{r}(t)}{\Delta t}

しかし、運動は時々刻々と変化するため、ある瞬間 $$t$$ における速度を知る必要がある。これは、 $\Delta t \to 0$ の極限をとることで得られる 瞬間速度 (Instantaneous velocity) である:

\mathbf{v}(t) = \lim_{\Delta t \to 0} \frac{\mathbf{r}(t + \Delta t) - \mathbf{r}(t)}{\Delta t} = \frac{d\mathbf{r}}{dt}

これは位置ベクトルの 時間微分 である。成分で書けば:

\mathbf{v}(t) = \frac{dx}{dt}\mathbf{e}_x + \frac{dy}{dt}\mathbf{e}_y + \frac{dz}{dt}\mathbf{e}_z = \dot{x}\mathbf{e}_x + \dot{y}\mathbf{e}_y + \dot{z}\mathbf{e}_z

ここで $\dot{x} = \frac{dx}{dt}$ という記法を導入した。これは ニュートンのドット記法 と呼ばれ、時間微分を簡潔に表すのに便利である。

重要な概念的理解 : 速度は ベクトル量 であり、大きさ(速さ)と方向の両方の情報を持つ。速度ベクトルの方向は、その瞬間の運動の方向を示す。

加速度：変化の変化

速度が一定であれば、物体は等速直線運動をする。しかし、多くの興味深い運動では速度自体が変化する。この「 速度の変化の割合 」を表すのが 加速度 (Acceleration) である。

速度の定義と同様に、加速度は速度ベクトルの時間微分として定義される:

\mathbf{a}(t) = \frac{d\mathbf{v}}{dt} = \frac{d^2 \mathbf{r}}{dt^2}

加速度は位置ベクトルの 第2次時間微分 でもある。ニュートンのドット記法では:

\mathbf{a}(t) = \ddot{\mathbf{r}}(t) = \ddot{x}\mathbf{e}_x + \ddot{y}\mathbf{e}_y + \ddot{z}\mathbf{e}_z

加速度の物理的意味 : 加速度がゼロでない場合、何らかの力が物体に作用している。実際、ニュートンの第2法則は、力と加速度を直接結びつける。これは後の章で詳しく扱う。

運動学的記述の階層性

ここまでで、運動を記述する3つの階層が明らかになった:

位置 $\mathbf{r}(t)$ : 物体がどこにあるか
速度 $\mathbf{v}(t) = \dot{\mathbf{r}}(t)$ : どれくらい速く、どの方向に動いているか
加速度 $\mathbf{a}(t) = \ddot{\mathbf{r}}(t)$ : 速度がどのように変化しているか

これらは微分の関係で結ばれている:

\text{位置} \xrightarrow{\text{微分}} \text{速度} \xrightarrow{\text{微分}} \text{加速度}

逆に、積分の関係もある:

\text{加速度} \xrightarrow{\text{積分}} \text{速度} \xrightarrow{\text{積分}} \text{位置}

この階層構造の理解は、運動方程式を解く際の基本となる。

ベクトルとスカラー：量の性質による分類

スカラー量とベクトル量の本質的違い

物理量は、その変換性質によって スカラー量 と ベクトル量 に分類される。この区別は単なる便宜上のものではなく、物理法則の記述において本質的である。

スカラー量 (Scalar quantity) は、大きさのみで完全に記述される量である。座標系の選び方によらず、その値は変わらない。例:

質量 $$m$$
時間 $$t$$
エネルギー $$E$$
温度 $$T$$

ベクトル量 (Vector quantity) は、大きさと方向の両方を持つ量である。座標系を変えると、その成分表示は変わるが、ベクトルそのものの幾何学的性質は不変である。例:

位置 $\mathbf{r}$
速度 $\mathbf{v}$
加速度 $\mathbf{a}$
力 $\mathbf{F}$

ベクトルの座標系非依存性

ベクトル量を理解する上で最も重要な概念は、 ベクトルそのものは座標系の選択に依存しない幾何学的対象である という点だ。

例えば、速度ベクトル $\mathbf{v}$ を考える。あるデカルト座標系で:

\mathbf{v} = v_x \mathbf{e}_x + v_y \mathbf{e}_y + v_z \mathbf{e}_z

別の座標系(例えば、原点は同じだが軸を回転させた座標系)では:

\mathbf{v} = v'_x \mathbf{e}'_x + v'_y \mathbf{e}'_y + v'_z \mathbf{e}'_z

成分 $$(v_x, v_y, v_z)$$ と $$(v'_x, v'_y, v'_z)$$ は異なるが、これらは同じベクトル $\mathbf{v}$ を表している。物体の実際の運動(速度の大きさと方向)は座標系の選択によらない。

この性質により、物理法則をベクトル形式で書くと、 座標系によらない普遍的な形 になる。これは理論の美しさと実用性の両方に寄与する。

ベクトルの演算

ベクトルに対しては以下の基本演算が定義される:

1. ベクトルの加法 :

\mathbf{a} + \mathbf{b} = (a_x + b_x)\mathbf{e}_x + (a_y + b_y)\mathbf{e}_y + (a_z + b_z)\mathbf{e}_z

幾何学的には、平行四辺形の法則に従う。

2. スカラー倍 :

c\mathbf{a} = (ca_x)\mathbf{e}_x + (ca_y)\mathbf{e}_y + (ca_z)\mathbf{e}_z

ベクトルの大きさを $$|c|$$ 倍し、 $$c < 0$$ なら向きを反転する。

3. 内積(スカラー積) :

\mathbf{a} \cdot \mathbf{b} = a_x b_x + a_y b_y + a_z b_z = |\mathbf{a}||\mathbf{b}|\cos\theta

ここで $\theta$ は2つのベクトルのなす角である。内積の結果はスカラー量である。

4. 外積(ベクトル積) :

\mathbf{a} \times \mathbf{b} = \begin{vmatrix} \mathbf{e}_x & \mathbf{e}_y & \mathbf{e}_z \\ a_x & a_y & a_z \\ b_x & b_y & b_z \end{vmatrix}

外積の結果は、 $\mathbf{a}$ と $\mathbf{b}$ の両方に垂直なベクトルである。その大きさは $|\mathbf{a}||\mathbf{b}|\sin\theta$ であり、これは $\mathbf{a}$ と $\mathbf{b}$ が張る平行四辺形の面積に等しい。

外積は特に回転運動の記述において中心的な役割を果たす。これは後の章で詳しく扱う。

質量と慣性：運動の「重さ」

質量の概念

質量 (Mass) は物体の基本的な属性であり、2つの重要な側面を持つ:

慣性質量 : 物体の運動状態を変えることの困難さを表す
重力質量 : 重力相互作用の強さを決める

一般相対性理論が示すように、これら2つは等価であるが、剛体力学では主に慣性質量の側面が重要である。

質量の単位はSI単位系ではキログラム(kg)である。質量はスカラー量であり、常に正の値を持つ。

慣性の意味

ニュートンの第1法則(慣性の法則)は、「外力が働かない限り、物体は静止または等速直線運動を続ける」と述べる。この「運動状態を維持しようとする性質」が慣性である。

質量が大きいほど、同じ力を加えても加速度は小さくなる。これはニュートンの第2法則:

\mathbf{F} = m\mathbf{a}

から明らかである。同じ力 $\mathbf{F}$ に対して、質量 $$m$$ が大きいほど加速度 $\mathbf{a} = \mathbf{F}/m$ は小さい。

質点から剛体へ：質量分布の重要性

質点 (Point mass) とは、質量を持つが空間的な広がりを持たない理想化された物体である。多くの基礎的な問題では、物体を質点として扱うことが有効である。

しかし、人体の歩行のような複雑な運動を記述するには、 物体の空間的な広がりと質量の分布 を考慮する必要がある。例えば:

人体の各部位(頭、胴体、腕、脚)はそれぞれ異なる質量を持つ
これらの部位の空間的配置が時間とともに変化する
回転運動において、質量がどのように分布しているかが決定的に重要になる

この質量分布の記述が、後に導入する重心や 慣性モーメント の概念につながる。

質量保存の原則

古典力学においては、 物体の質量は時間によらず一定 であると仮定する。これは質量保存則と呼ばれる。数学的には:

\frac{dm}{dt} = 0

この仮定は、相対論的効果が無視できる速度範囲や、質量の流入・流出がない閉じた系において有効である。剛体力学では常にこの仮定が成り立つ。

力の概念：運動変化の原因

力とは何か

力 (Force) は、物体の運動状態を変化させる原因である。より厳密には、力は物体に加速度を生じさせる物理的作用である。

力の定義には2つのアプローチがある:

運動学的定義 : ニュートンの第2法則 $\mathbf{F} = m\mathbf{a}$ により、力を加速度から定義する
相互作用的定義 : 物体間の相互作用として力を定義する(ニュートンの第3法則が関連)

どちらのアプローチも有効だが、実用的には両方の視点を持つことが重要である。

力のベクトル性

力は ベクトル量 である。これは力が以下の性質を持つことを意味する:

大きさ : どれだけ強く押す/引くか(単位: ニュートン N)
方向 : どの向きに押す/引くか
作用点 : どこに力が作用するか(後に重要になる)

同じ点に複数の力が作用する場合、それらは ベクトル和 として合成される:

\mathbf{F_{\text{net}} }= \mathbf{F_1} + \mathbf{F_2} + \cdots + \mathbf{F_n} = \sum_{i=1}^{n} \mathbf{F_i}

この合力 $\mathbf{F}_{\text{net}}$ が物体の運動を決定する。

力の分類

物理学における力は、その起源によって分類できる:

基本的な力 :

重力
電磁気力
強い核力
弱い核力

見かけの力(接触力など) :

垂直抗力
摩擦力
張力
弾性力

剛体力学では、主に後者の「見かけの力」を扱う。これらは実際には電磁気力の巨視的な現れだが、その微視的起源を考えずに現象論的に扱うことができる。

力と運動の因果関係

初学者がしばしば混同するのは、「力」と「運動」の関係である。重要な点:

力がなければ加速度はゼロ だが、速度はゼロとは限らない(等速運動は力なしで持続する)
静止している物体 には、合力がゼロであるか、全く力が作用していないかのどちらかである
一定速度で動いている物体 にも、合力はゼロである

つまり、力は「運動」そのものではなく「運動の変化」を引き起こす。

微分方程式としての運動方程式

なぜ微分方程式なのか

ニュートンの第2法則 $\mathbf{F} = m\mathbf{a}$ を位置の2階微分で書き直すと:

\mathbf{F} = m\frac{d^2\mathbf{r}}{dt^2}

これは 2階常微分方程式 である。なぜ運動方程式は微分方程式の形をとるのか?

答えは、 物理法則は瞬間的な因果関係を記述する からである。力 $\mathbf{F}$ がある瞬間 $$t$$ において物体に作用すると、その瞬間の加速度 $\mathbf{a}(t)$ が決まる。未来の運動は、この瞬間的な関係を時間的に積分することで予測される。

初期値問題

微分方程式を解くには、 初期条件 が必要である。2階微分方程式の場合、以下の2つの情報が必要:

初期位置 : $\mathbf{r}(t_0) = \mathbf{r}_0$
初期速度 : $\mathbf{v}(t_0) = \mathbf{v}_0$

これらが与えられると、運動方程式を解いて、すべての時刻 $$t > t_0$$ における位置 $\mathbf{r}(t)$ を求めることができる(原理的には)。

物理的意味 : 初期条件は「過去の歴史」を要約する。物体が時刻 $$t_0$$ にどこにあり、どのように動いていたかが分かれば、それ以降の運動が(力が既知なら)完全に決定される。これが古典力学の 決定論的 性質である。

解の一意性と存在

適切な条件下では、初期値問題の解は 一意に存在する。これは数学的には常微分方程式の存在一意性定理によって保証される。

物理的には、これは「同じ初期状態と同じ力の下では、常に同じ運動が生じる」ことを意味する。この再現性が、実験物理学と理論物理学を結びつける。

解析解と数値解

微分方程式を解く方法には2種類ある:

解析解 (Analytical solution) : 閉じた数学的形式で解を表現する。例えば、一定の力 $\mathbf{F}$ の下での運動:

\mathbf{r}(t) = \mathbf{r}_0 + \mathbf{v}_0 t + \frac{1}{2}\frac{\mathbf{F}}{m}t^2

数値解 (Numerical solution) : コンピュータを用いて、離散的な時刻での解を計算する。複雑な問題(例えば人体の歩行)では、解析解を得ることはほぼ不可能であり、数値解に頼る。

本教科書の最終目標である複雑な多体系の記述では、運動方程式を正しく定式化することが主眼であり、その解は数値的に求めることを前提とする。

次章への橋渡し

本章では、運動を記述するための概念的・数学的基礎を確立した:

座標系と時間 : 観察の枠組み
位置・速度・加速度 : 運動状態の階層的記述
ベクトルとスカラー : 物理量の分類
質量 : 慣性の尺度
力 : 運動変化の原因
微分方程式 : 物理法則の数学的表現

これらの概念は、一見基礎的に見えるが、複雑な多体系を扱う際にも常に立ち返る基盤となる。

次章では、これらの道具を使って 単一物体の運動 を具体的に記述する。まずは質点として扱える物体の並進運動から始め、運動方程式の立て方と解き方の基本的なパターンを学ぶ。その後、回転運動という新しい複雑性を導入し、最終的には空間的広がりを持つ剛体の記述へと進む。

各ステップで、「なぜこの概念が必要なのか」「どのような問題を解決するのか」を明確にしながら、複雑な多体系の記述能力を段階的に構築していく。

単一物体の運動記述における基本要素

3.1 なぜ単一物体から始めるのか：問題意識の明確化

複雑な多体システム（例えば人体の歩行）を理解するための最初のステップは、 単一の物体がどのように運動するか を正確に記述することです。これは単なる教育的便宜ではなく、物理学における本質的な方法論です。

問題の階層性

人体のような複雑なシステムも、究極的には「複数の物体が相互作用しながら運動している集合体」です。この複雑さに直面する前に、まず以下の基本的な問いに答える必要があります:

一つの物体の「位置」を数学的にどう表現するのか?
物体の「運動」とは、数学的に何を意味するのか?
物体に働く「力」と運動との関係はどう記述されるのか?

これらの問いに正確に答えられなければ、複数の物体が絡み合う系を記述することは不可能です。

3.2 質点近似：最も単純化された物体の表現

3.2.1 質点とは何か：概念の本質

質点（point mass） とは、物体が持つすべての質量が一点に集中していると仮定した理想化モデルです。この仮定には深い意味があります:

空間的広がりの無視 : 物体の形状、サイズ、内部構造をすべて無視する
回転の不在 : 大きさがないため、回転という概念が存在しない
位置の一意性 : 物体全体の位置が一つの点の座標だけで完全に記述できる

この極端な単純化が有効なのは、以下のような状況です:

物体のサイズが運動のスケールに比べて十分小さい場合
- 例: 太陽系における惑星の運動（惑星の直径 << 惑星間距離）
回転運動が問題の本質でない場合
- 例: 放物運動する物体の軌道計算
並進運動のみに着目したい場合
- 例: 物体の重心の軌道追跡

3.2.2 質点の位置と運動の記述

質点の位置は、選択した座標系における 位置ベクトル $\mathbf{r}(t)$ によって記述されます。3次元空間では:

\mathbf{r}(t) = \begin{pmatrix} x(t) \\ y(t) \\ z(t) \end{pmatrix}

ここで重要なのは、 時刻 $$t$$ の関数として記述されている 点です。運動とは、時間の経過に伴う位置の変化そのものです。

速度は位置ベクトルの時間微分として定義されます:

\mathbf{v}(t) = \frac{d\mathbf{r}(t)}{dt} = \dot{\mathbf{r}}(t)

加速度 はさらにその時間微分です:

\mathbf{a}(t) = \frac{d\mathbf{v}(t)}{dt} = \frac{d^2\mathbf{r}(t)}{dt^2} = \ddot{\mathbf{r}}(t)

3.2.3 ニュートンの運動方程式：力と運動の関係

質点に働く力 $\mathbf{F}$ と加速度 $\mathbf{a}$ の関係は、ニュートンの第二法則によって記述されます:

\mathbf{F} = m\mathbf{a} = m\ddot{\mathbf{r}}

この式は 2階の常微分方程式 です。これを解くことは、与えられた力の下で質点がどのように運動するかを決定することを意味します。

方程式の意味を正確に理解する:

左辺 $\mathbf{F}$ : 質点に働くすべての力の合力
右辺 $m\ddot{\mathbf{r}}$ : 質量と加速度の積（慣性力）
等号 : 力が加わると、それに応じた加速度が生じるという因果関係

3.2.4 質点モデルの限界：なぜ拡張が必要か

質点モデルは以下の状況を記述できません:

物体が回転する場合
- 例: 体操選手の宙返り、こまの回転
物体の異なる部分に異なる力が働く場合
- 例: 長い棒の両端に異なる力が作用
物体の形状や向きが運動に影響する場合
- 例: 歩行における足の向き、飛行機の姿勢制御

人体の歩行は明らかに回転運動を含みます（足の振り、腕の振り、胴体の回転など）。したがって、質点モデルの拡張が不可欠です。

3.3 質点から有限サイズ物体への拡張：何が変わるのか

3.3.1 物体の広がりがもたらす新しい自由度

物体に有限の大きさを認めると、その位置を記述するには 位置だけでは不十分 になります。なぜなら、物体は空間内で 向き（姿勢、orientation） を変えることができるからです。

例で理解する:

質点としての本: 本の重心の3次元座標 $$(x, y, z)$$ のみ → 3自由度
有限サイズの本: 重心の位置 $$(x, y, z)$$ + 本の向き（3つの角度）→ 6自由度

この6自由度は、一般的な3次元空間における剛体の運動を完全に記述するために必要な最小限の情報です。

3.3.2 並進と回転の分離：重心の特別な役割

有限サイズの物体の運動を理解する鍵は、 並進運動 と 回転運動 の分離です。

重心（center of mass）の定義

物体を構成する微小要素の質量を $$dm$$ 、その位置を $\mathbf{r}'$ とすると、重心の位置 $\mathbf{r}_{\text{cm}}$ は:

\mathbf{r}_{\text{cm}} = \frac{1}{M}\int \mathbf{r}' \, dm

ここで $$M$$ は物体の全質量です。離散的な質点系の場合:

\mathbf{r}_{\text{cm}} = \frac{1}{M}\sum_{i} m_i \mathbf{r}_i, \quad M = \sum_{i} m_i

なぜ重心が特別なのか

重心には以下の重要な性質があります:

性質1: 重心の運動方程式

物体に働くすべての外力の合力を $\mathbf{F}_{\text{ext}}$ とすると:

\mathbf{F}_{\text{ext}} = M\ddot{\mathbf{r}}_{\text{cm}}

これは、 物体全体の並進運動は、全質量が重心に集中した質点の運動と同じ方程式に従う ことを意味します。

性質2: 重心周りの運動の独立性

物体の回転運動は、重心を基準として記述するのが最も自然です。なぜなら、外力による回転（トルク）と並進が数学的に分離できるからです。

3.3.3 並進運動と回転運動の記述的独立性

物体の任意の点の位置 $\mathbf{r}$ は、以下のように分解できます:

\mathbf{r} = \mathbf{r}_{\text{cm}} + \mathbf{r}'

ここで:

$\mathbf{r}_{\text{cm}}$ : 重心の位置（並進成分）
$\mathbf{r}'$ : 重心から測った相対位置（回転成分）

この分解によって:

並進運動 : $\mathbf{r}_{\text{cm}}$ の時間発展として記述
回転運動 : $\mathbf{r}'$ の時間発展として記述

という独立した扱いが可能になります。

3.4 回転運動の記述：新しい数学的道具の必要性

3.4.1 回転を記述する困難さ

回転運動を記述するには、質点の運動記述では不要だった新しい概念が必要です:

角度の表現 : 物体の向きをどう数値化するか?
角速度 : 回転の速さをどう定義するか?
角運動量 : 回転の慣性をどう表現するか?
トルク : 回転を引き起こす力の効果をどう測るか?

3.4.2 角速度ベクトル：回転の瞬間的記述

物体が回転しているとき、その 角速度ベクトル $\boldsymbol{\omega}$ は以下の性質を持ちます:

方向 : 回転軸の方向（右手の法則による）
大きさ : 単位時間あたりの回転角度 $|\boldsymbol{\omega}| = \omega$

重心からの相対位置 $\mathbf{r}'$ にある点の速度は:

\mathbf{v}' = \boldsymbol{\omega} \times \mathbf{r}'

この外積（クロス積）の形は、回転運動の幾何学的性質を反映しています:

速度は回転軸に垂直
速度の大きさは回転軸からの距離に比例

3.4.3 角運動量：回転運動の慣性を表現する量

並進運動における運動量 $\mathbf{p} = m\mathbf{v}$ に対応する概念が、回転運動における 角運動量 $\mathbf{L}$ です。

重心周りの角運動量は:

\mathbf{L} = \int \mathbf{r}' \times \mathbf{v}' \, dm

角速度 $\boldsymbol{\omega}$ を用いると:

\mathbf{L} = \int \mathbf{r}' \times (\boldsymbol{\omega} \times \mathbf{r}') \, dm

この積分計算を実行すると、 慣性テンソル $\mathbf{I}$ という新しい量が現れます:

\mathbf{L} = \mathbf{I} \boldsymbol{\omega}

3.4.4 慣性テンソル：回転における「質量」の一般化

並進運動では質量 $$m$$ が一つのスカラー量でしたが、回転運動ではそれに対応する 慣性テンソル $\mathbf{I}$ は $3 \times 3$ の行列です:

\mathbf{I} = \begin{pmatrix} I_{xx} & I_{xy} & I_{xz} \\ I_{yx} & I_{yy} & I_{yz} \\ I_{zx} & I_{zy} & I_{zz} \end{pmatrix}

各成分は以下のように定義されます:

\begin{aligned} I_{xx} &= \int (y'^2 + z'^2) \, dm \\ I_{yy} &= \int (x'^2 + z'^2) \, dm \\ I_{zz} &= \int (x'^2 + y'^2) \, dm \\ I_{xy} &= I_{yx} = -\int x'y' \, dm \\ I_{xz} &= I_{zx} = -\int x'z' \, dm \\ I_{yz} &= I_{zy} = -\int y'z' \, dm \end{aligned}

なぜテンソルが必要なのか：概念的理解

質量がスカラーである理由は、並進運動の方向によらず同じ慣性を示すからです。しかし回転運動では、 回転軸の方向によって回転のしやすさが異なる ため、方向依存性を表現する行列（テンソル）が必要になります。

例で理解する:

細長い棒を長軸周りに回転させる: 回しやすい（小さな慣性モーメント）
同じ棒を短軸周りに回転させる: 回しにくい（大きな慣性モーメント）

この方向依存性が慣性テンソルの非対角成分（ $I_{xy}$ など）にも反映されます。

3.4.5 回転運動の方程式：オイラー方程式への道

角運動量の時間変化は、物体に働く トルク （力のモーメント） $\boldsymbol{\tau}$ によって決まります:

\frac{d\mathbf{L}}{dt} = \boldsymbol{\tau}

これは並進運動における $\frac{d\mathbf{p}}{dt} = \mathbf{F}$ に対応します。

$\mathbf{L} = \mathbf{I}\boldsymbol{\omega}$ を代入すると:

\frac{d(\mathbf{I}\boldsymbol{\omega})}{dt} = \boldsymbol{\tau}

ここで重要な複雑さが現れます: 慣性テンソル $\mathbf{I}$ 自体が時間とともに変化する可能性があります（物体が回転すると、空間固定座標系から見た慣性テンソルの成分が変わる）。

この複雑さを扱うために、後の章で 物体固定座標系 という概念を導入します。

3.5 トルク：回転を引き起こす力の効果

3.5.1 トルクの定義と物理的意味

力 $\mathbf{F}$ が位置 $\mathbf{r}$ に作用するとき、基準点周りの トルク （モーメント）は:

\boldsymbol{\tau} = \mathbf{r} \times \mathbf{F}

トルクのベクトルは:

方向 : 回転軸の方向を示す（右手の法則）
大きさ : $|\boldsymbol{\tau}| = |\mathbf{r}||\mathbf{F}|\sin\theta$ （ $\theta$ は $\mathbf{r}$ と $\mathbf{F}$ の間の角度）

3.5.2 トルクの加法性と合力

複数の力 $\mathbf{F}_i$ が位置 $\mathbf{r}_i$ に作用する場合、全トルクは:

\boldsymbol{\tau}_{\text{total}} = \sum_{i} \mathbf{r}_i \times \mathbf{F}_i

重要な注意点: トルクは 基準点の選択に依存 します。ただし、重心周りのトルクを計算する場合、以下の重要な定理が成り立ちます:

定理 : 外力による重心周りのトルクは、重心の加速度に依存しない（内力のトルクは互いに打ち消し合う）。

3.6 単一物体記述のまとめ：何を獲得したか

この章で構築した枠組みによって、以下が可能になりました:

3.6.1 完全な運動記述の要素

並進運動 :
- 記述変数: 重心位置 $\mathbf{r}_{\text{cm}}$
- 運動方程式: $\mathbf{F}_{\text{ext}} = M\ddot{\mathbf{r}}_{\text{cm}}$
- パラメータ: 全質量 $$M$$
回転運動 :
- 記述変数: 角速度 $\boldsymbol{\omega}$ （または角運動量 $\mathbf{L}$ ）
- 運動方程式: $\frac{d\mathbf{L}}{dt} = \boldsymbol{\tau}$
- パラメータ: 慣性テンソル $\mathbf{I}$

3.6.2 未解決の課題：次章への橋渡し

しかし、以下の重要な問題がまだ残っています:

物体の向きの記述 : 角速度だけでは物体の現在の姿勢が分からない
座標系の選択問題 : 慣性テンソルが時間変化する複雑さ
剛体の厳密な定義 : 「形が変わらない」とはどういう条件か?

これらの問題に対処するため、次章では 回転運動のより深い数学的構造 と 剛体の厳密な定義 に踏み込みます。特に、物体の姿勢を記述するための数学的道具（回転行列、オイラー角、四元数など）と、物体固定座標系の概念を導入します。

3.6.3 学習の確認：この章で理解すべき核心

以下の問いに自分の言葉で答えられるか確認してください:

質点モデルと有限サイズ物体の決定的な違いは何か?
なぜ重心を基準にすると並進と回転が分離できるのか?
慣性テンソルがスカラーではなく行列である物理的理由は?
トルクが「力×距離」以上の意味を持つのはなぜか?

これらの理解が、複雑な多体システムを記述するための確固たる基盤となります。

回転運動の導入と複雑性の認識

問題意識：なぜ回転運動が特別なのか

前章までで、質点および質点系の並進運動を記述する枠組みを構築してきた。重心運動の方程式 $M\ddot{\boldsymbol{r}}_G = \boldsymbol{F}_{\text{ext}}$ は、系全体がまるで一点に質量が集中しているかのように運動する様子を記述する。

しかし、 この記述だけでは物体の運動を完全には捉えきれない。考えてみよう：

回転する独楽（コマ）は、重心が静止していても複雑な運動を続ける
人体の歩行では、各部位が回転しながら全体として前進する
投げられた棒は、重心が放物線を描きながら、同時に回転する

重心運動の方程式は「どこに移動するか」を教えてくれるが、「どの向きを向いているか」「どのように回転しているか」については何も語らない。 物体の姿勢と回転の記述が欠けている のである。

本章では、この欠落を埋めるために回転運動の記述に踏み込む。ただし、ここで重要なのは、 回転運動が並進運動よりも本質的に複雑である理由を理解すること である。この複雑性を認識することが、剛体という概念の必然性を理解する鍵となる。

回転とは何か：角度と角速度の幾何学

回転の記述における根本的課題

並進運動では、位置ベクトル $\boldsymbol{r}$ とその時間微分である速度 $\boldsymbol{v} = \dot{\boldsymbol{r}}$ で運動を記述できた。では、回転運動はどうだろうか？

問題は、「回転の量」をどう定義するかである。平面内の回転なら角度 $\theta$ で記述できそうだが、三次元空間での回転は単純ではない。物体がどの軸の周りに、どれだけ回転したかを指定する必要がある。

角速度ベクトルの導入

回転運動を記述する鍵となる概念が 角速度ベクトル $\boldsymbol{\omega}$ である。これは単なる角度の時間微分ではなく、以下の情報を同時に含む：

大きさ $|\boldsymbol{\omega}|$ ：単位時間あたりの回転角度（rad/s）
向き：回転軸の方向（右ねじの進む方向）

この定義により、三次元空間での任意の回転運動を一つのベクトルで表現できる。

重要な性質 ：角速度ベクトル $\boldsymbol{\omega}$ を持つ物体上の点 $\boldsymbol{r}$ （回転軸からの位置ベクトル）の速度は：

\boldsymbol{v} = \boldsymbol{\omega} \times \boldsymbol{r}

この外積の関係は、回転運動の幾何学を端的に表現している。ベクトル $\boldsymbol{v}$ は $\boldsymbol{\omega}$ と $\boldsymbol{r}$ の両方に垂直で、大きさは $|\boldsymbol{v}| = |\boldsymbol{\omega}||\boldsymbol{r}|\sin\phi$ （ $\phi$ は $\boldsymbol{\omega}$ と $\boldsymbol{r}$ のなす角）となる。

角速度の非可換性

注意すべき重要な性質 ：三次元空間での回転は 順序依存的 である。つまり、異なる軸周りの二つの回転を施す順序を入れ替えると、結果が変わる。

例えば、本を手に取って以下を試してみよう：

x軸周りに90度回転してから、y軸周りに90度回転する
y軸周りに90度回転してから、x軸周りに90度回転する

最終的な本の向きが異なることが分かる。数学的には、回転変換が非可換であることを意味する。

この非可換性は、角速度ベクトルの足し算には影響しない（瞬間的には可換）が、 有限角度の回転を扱う際には本質的な複雑さ となる。この点は、後に姿勢の記述で再び問題となる。

回転運動エネルギーと慣性モーメントの必然性

なぜ新しい概念が必要か

並進運動の運動エネルギーは $T_{\text{trans}} = \frac{1}{2}M\boldsymbol{v}_G^2$ と、質量という一つのスカラー量で記述できた。では、回転運動のエネルギーはどうだろうか？

物体が角速度 $\boldsymbol{\omega}$ で回転しているとき、物体内の各質点 $$i$$ （質量 $$m_i$$ 、位置 $\boldsymbol{r}_i$ ）は速度 $\boldsymbol{v}_i = \boldsymbol{\omega} \times \boldsymbol{r}_i$ を持つ。したがって、回転運動エネルギーは：

T_{\text{rot}} = \sum_i \frac{1}{2}m_i |\boldsymbol{v}_i|^2 = \sum_i \frac{1}{2}m_i |\boldsymbol{\omega} \times \boldsymbol{r}_i|^2

この式を展開すると、単純に「質量 × 角速度の2乗」の形にはならない。 質点の位置 $\boldsymbol{r}_i$ がどのように分布しているかが本質的に関わってくる のである。

単純な例：固定軸周りの回転

まず、z軸周りに回転する場合を考えよう。このとき $\boldsymbol{\omega} = \omega\boldsymbol{e}_z$ （ $\boldsymbol{e}_z$ はz軸方向の単位ベクトル）とすると、質点 $$i$$ の回転軸からの距離を $r_{i\perp} = \sqrt{x_i^2 + y_i^2}$ として：

|\boldsymbol{v}_i|^2 = \omega^2 r_{i\perp}^2

したがって：

T_{\text{rot}} = \frac{1}{2}\left(\sum_i m_i r_{i\perp}^2\right)\omega^2 = \frac{1}{2}I_z\omega^2

ここで定義された量：

I_z = \sum_i m_i r_{i\perp}^2

が、z軸周りの 慣性モーメント （moment of inertia）である。

慣性モーメントの物理的意味 ：これは「質量が回転軸からどれだけ離れて分布しているか」を表す量である。同じ質量でも、軸から遠くに分布していれば慣性モーメントは大きくなる。つまり、 回転のしにくさ を表す量と解釈できる。

一般の回転：慣性テンソルの必要性

しかし、回転軸が座標軸と一致しない一般的な場合、状況はより複雑になる。任意の角速度ベクトル $\boldsymbol{\omega} = (\omega_x, \omega_y, \omega_z)$ に対して、外積 $\boldsymbol{\omega} \times \boldsymbol{r}_i$ を計算し、その2乗を取ると：

|\boldsymbol{\omega} \times \boldsymbol{r}_i|^2 = \omega_x^2(y_i^2 + z_i^2) + \omega_y^2(z_i^2 + x_i^2) + \omega_z^2(x_i^2 + y_i^2) - 2\omega_x\omega_y x_i y_i - 2\omega_y\omega_z y_i z_i - 2\omega_z\omega_x z_i x_i

この式を全質点について和をとると、 $\omega_x^2$ , $\omega_y^2$ , $\omega_z^2$ だけでなく、 $\omega_x\omega_y$ , $\omega_y\omega_z$ , $\omega_z\omega_x$ といった交差項が現れる。

これを行列形式で表現すると：

T_{\text{rot}} = \frac{1}{2}\boldsymbol{\omega}^T \mathbf{I} \boldsymbol{\omega} = \frac{1}{2}\sum_{j,k}\omega_j I_{jk} \omega_k

ここで $\mathbf{I}$ は 慣性テンソル （inertia tensor）と呼ばれる $3 \times 3$ 行列で、その成分は：

I_{jk} = \sum_i m_i \left(\delta_{jk}|\boldsymbol{r}_i|^2 - r_{ij}r_{ik}\right)

または、より明示的に：

\mathbf{I} = \begin{pmatrix} I_{xx} & -I_{xy} & -I_{xz} \\ -I_{yx} & I_{yy} & -I_{yz} \\ -I_{zx} & -I_{zy} & I_{zz} \end{pmatrix}

ここで：

\begin{align} I_{xx} &= \sum_i m_i(y_i^2 + z_i^2), \quad I_{yy} = \sum_i m_i(z_i^2 + x_i^2), \quad I_{zz} = \sum_i m_i(x_i^2 + y_i^2) \\ I_{xy} &= I_{yx} = \sum_i m_i x_i y_i, \quad I_{yz} = I_{zy} = \sum_i m_i y_i z_i, \quad I_{zx} = I_{xz} = \sum_i m_i z_i x_i \end{align}

対角成分 $I_{xx}, I_{yy}, I_{zz}$ は各軸周りの慣性モーメント、非対角成分 $I_{xy}, I_{yz}, I_{zx}$ は 慣性乗積 （products of inertia）と呼ばれる。

重要な認識 ：回転運動のエネルギーを記述するには、もはや一つのスカラー量（質量）では不足で、 9個の成分を持つテンソル量 が必要になる。これが、回転運動が並進運動よりも本質的に複雑である第一の理由である。

角運動量：回転運動の「勢い」

角運動量の定義と物理的意味

並進運動における運動量 $\boldsymbol{p} = M\boldsymbol{v}$ に対応する、回転運動の「勢い」を表す量が 角運動量 （angular momentum）である。

基準点 $$O$$ に対する質点 $$i$$ の角運動量は：

\boldsymbol{L}_i = \boldsymbol{r}_i \times m_i\boldsymbol{v}_i

系全体の角運動量は：

\boldsymbol{L} = \sum_i \boldsymbol{r}_i \times m_i\boldsymbol{v}_i

物理的意味 ：角運動量の大きさは、「どれだけ遠くで（ $|\boldsymbol{r}_i|$ ）、どれだけ速く（ $|\boldsymbol{v}_i|$ ）、どれだけの質量（ $$m_i$$ ）が運動しているか」を表す。外積の性質から、動径ベクトルと速度が平行な成分は寄与せず、 回転的な運動成分のみが角運動量に寄与する。

角運動量と角速度の関係

物体が回転運動している場合、 $\boldsymbol{v}_i = \boldsymbol{\omega} \times \boldsymbol{r}_i$ を代入すると：

\boldsymbol{L} = \sum_i m_i \boldsymbol{r}_i \times (\boldsymbol{\omega} \times \boldsymbol{r}_i)

ベクトル三重積の公式 $\boldsymbol{a} \times (\boldsymbol{b} \times \boldsymbol{c}) = \boldsymbol{b}(\boldsymbol{a} \cdot \boldsymbol{c}) - \boldsymbol{c}(\boldsymbol{a} \cdot \boldsymbol{b})$ を用いると：

\boldsymbol{r}_i \times (\boldsymbol{\omega} \times \boldsymbol{r}_i) = \boldsymbol{\omega}|\boldsymbol{r}_i|^2 - \boldsymbol{r}_i(\boldsymbol{r}_i \cdot \boldsymbol{\omega})

これを成分表示で整理すると：

\boldsymbol{L} = \mathbf{I}\boldsymbol{\omega}

つまり、 角運動量は角速度に比例するが、その比例係数が慣性テンソルというテンソル量になる。

重要な帰結 ：一般に $\boldsymbol{L}$ と $\boldsymbol{\omega}$ は 異なる方向を向く。並進運動では運動量と速度は常に同じ方向だったが、回転運動では、物体の質量分布に応じて角運動量の方向が角速度の方向からずれる。この非直感的な性質が、回転運動の複雑さの核心である。

主軸と主慣性モーメント

慣性テンソル $\mathbf{I}$ は対称行列であるため、 適切な座標系を選べば対角化できる。この特別な座標系の軸を主軸（principal axes）と呼び、対角化された慣性テンソルの固有値を 主慣性モーメント $$I_1, I_2, I_3$$ と呼ぶ：

\mathbf{I}_{\text{principal}} = \begin{pmatrix} I_1 & 0 & 0 \\ 0 & I_2 & 0 \\ 0 & 0 & I_3 \end{pmatrix}

主軸座標系では、慣性乗積がすべてゼロになり、角運動量と角速度の関係が単純化される：

\boldsymbol{L} = (I_1\omega_1, I_2\omega_2, I_3\omega_3)

主軸の物理的意味 ：主軸の一つの周りに回転するとき、そのときに限り、角運動量と角速度が同じ方向を向く。物体の対称性がある場合（例：球、円柱、直方体）、主軸は対称軸と一致する。

主軸の概念は、回転運動の解析を大幅に簡略化する。しかし、一般的な物体では主軸は空間に固定されておらず、物体とともに回転する。この事実が、後述するオイラーの運動方程式の複雑さにつながる。

回転運動方程式：トルクの役割

トルクの定義と並進運動との対応

並進運動では $\boldsymbol{F} = \dot{\boldsymbol{p}}$ が基本方程式だった。回転運動における対応する方程式を導出しよう。

基準点 $$O$$ に関する トルク （torque）は：

\boldsymbol{N} = \sum_i \boldsymbol{r}_i \times \boldsymbol{F}_i

角運動量の時間微分を計算すると：

\frac{d\boldsymbol{L}}{dt} = \sum_i \frac{d}{dt}(\boldsymbol{r}_i \times m_i\boldsymbol{v}_i) = \sum_i \left(\dot{\boldsymbol{r}}_i \times m_i\boldsymbol{v}_i + \boldsymbol{r}_i \times m_i\dot{\boldsymbol{v}}_i\right)

第一項は $\dot{\boldsymbol{r}}_i = \boldsymbol{v}_i$ なので $\boldsymbol{v}_i \times \boldsymbol{v}_i = \boldsymbol{0}$ となり消える。第二項で $m_i\dot{\boldsymbol{v}}_i = \boldsymbol{F}_i$ を用いると：

\frac{d\boldsymbol{L}}{dt} = \sum_i \boldsymbol{r}_i \times \boldsymbol{F}_i = \boldsymbol{N}

これが 角運動量の変化の法則 （回転運動の第二法則）である：

\boxed{\boldsymbol{N} = \frac{d\boldsymbol{L}}{dt}}

この式は $\boldsymbol{F} = \frac{d\boldsymbol{p}}{dt}$ と完全に対応している。トルクが角運動量を変化させる原因であり、トルクがゼロなら角運動量は保存される。

オイラーの運動方程式

角速度で表現された回転運動方程式を導出しよう。 $\boldsymbol{L} = \mathbf{I}\boldsymbol{\omega}$ を時間微分すると：

\boldsymbol{N} = \frac{d\boldsymbol{L}}{dt} = \frac{d(\mathbf{I}\boldsymbol{\omega})}{dt} = \mathbf{I}\frac{d\boldsymbol{\omega}}{dt} + \frac{d\mathbf{I}}{dt}\boldsymbol{\omega}

ここに大きな問題がある ：慣性テンソル $\mathbf{I}$ は座標系に依存する量であり、物体が回転すると、空間固定座標系で見た $\mathbf{I}$ の成分は時間変化する（ $\frac{d\mathbf{I}}{dt} \neq \boldsymbol{0}$ ）。

この困難を回避するため、 物体に固定された座標系 （物体系）で方程式を記述する。物体系では $\mathbf{I}$ の成分は時間不変だが、代わりにベクトルの時間微分が複雑になる。空間固定系での微分と物体系での微分の関係は：

\left(\frac{d\boldsymbol{A}}{dt}\right)_{\text{space}} = \left(\frac{d\boldsymbol{A}}{dt}\right)_{\text{body}} + \boldsymbol{\omega} \times \boldsymbol{A}

これを $\boldsymbol{L}$ に適用すると：

\boldsymbol{N} = \dot{\boldsymbol{L}}_{\text{body}} + \boldsymbol{\omega} \times \boldsymbol{L}

主軸座標系を物体系に選ぶと、 $\boldsymbol{L} = (I_1\omega_1, I_2\omega_2, I_3\omega_3)$ なので：

\dot{\boldsymbol{L}}_{\text{body}} = (I_1\dot{\omega}_1, I_2\dot{\omega}_2, I_3\dot{\omega}_3)

そして $\boldsymbol{\omega} \times \boldsymbol{L}$ を計算すると、 オイラーの運動方程式 （Euler's equations of motion）が得られる：

\boxed{ \begin{align} N_1 &= I_1\dot{\omega}_1 - (I_2 - I_3)\omega_2\omega_3 \\ N_2 &= I_2\dot{\omega}_2 - (I_3 - I_1)\omega_3\omega_1 \\ N_3 &= I_3\dot{\omega}_3 - (I_1 - I_2)\omega_1\omega_2 \end{align} }

この方程式の重要な特徴 ：

非線形性 ：右辺に $\omega_i\omega_j$ という積の項が現れる。これは、ある軸周りの回転が他の軸周りの角加速度に影響を与えることを示す。
連成性 ：三つの方程式が互いに結合している。一つの軸だけを独立に扱うことはできない。
主慣性モーメントの差が本質的 ： $$(I_2 - I_3)$$ などの項が、異なる軸間の相互作用を引き起こす。

この非線形連成方程式こそが、回転運動が並進運動よりも解析的に困難である本質的理由である。

回転運動の複雑性の本質

ここまでの考察から、 回転運動が並進運動と比べて本質的に複雑である理由 を整理しよう。

複雑性の源泉

記述量の増大
- 並進運動：質量（スカラー1個）と速度ベクトル（3成分）
- 回転運動：慣性テンソル（独立な成分6個）と角速度ベクトル（3成分）
方向依存性
- 並進運動：運動量 $\boldsymbol{p}$ と速度 $\boldsymbol{v}$ は常に同じ方向
- 回転運動：角運動量 $\boldsymbol{L}$ と角速度 $\boldsymbol{\omega}$ は一般に異なる方向
座標系依存性
- 並進運動：質量は座標変換で不変
- 回転運動：慣性テンソルは座標回転で変換される
運動方程式の非線形性
- 並進運動： $\boldsymbol{F} = M\dot{\boldsymbol{v}}$ （線形）
- 回転運動：オイラー方程式（角速度に関して非線形）
有限回転の非可換性
- 並進運動：変位の順序は結果に影響しない
- 回転運動：異なる軸周りの回転の順序が結果を変える

なぜこれらの複雑性が生じるのか

根本的な理由は、 回転が幾何学的に複雑な変換だから である。

並進運動は、空間のすべての点を同じ方向に同じ距離だけ移動させる。これは最も単純な変換（平行移動）である。

一方、回転運動は、ある軸上の点は固定し、軸から離れるほど大きく移動させる。さらに、三次元空間での回転は三つの独立なパラメータ（どの軸の周りに、どれだけ回転するか）が必要で、これらのパラメータの組み合わせ方が結果に影響する。

質量分布の影響 も決定的である。並進運動では、質量がどこに分布していようと、総質量さえ分かれば運動が決まる。しかし回転運動では、 同じ総質量でも、その空間分布によって運動が全く異なる。軸から遠くに質量が分布していれば回転しにくく（大きな慣性モーメント）、軸の近くなら回転しやすい。

剛体概念への道筋

変形の問題

ここまでの議論で、実は暗黙のうちに重要な仮定をしてきた。それは、 物体の形状が時間変化しない という仮定である。

もし物体が変形すれば：

各質点の位置 $\boldsymbol{r}_i$ が時間変化する
慣性テンソル $\mathbf{I}$ が時間変化する
主軸の方向が時間変化する

このとき、オイラーの運動方程式は使えなくなる。なぜなら、 $\mathbf{I}$ が時間不変であることを前提にして導出されたからである。

剛体という抽象化の必然性

複雑な物体の回転運動を実用的に記述するためには、 物体内の質点間の距離が変化しない という理想化が必要になる。これが剛体（rigid body）の定義である：

剛体とは、任意の二点間の距離が時間的に不変な物体である。

この制約のもとでは：

物体の形状を記述する内部座標は不要になる
物体の配置は、重心位置（3自由度）と姿勢（3自由度）の計6自由度で完全に記述できる
慣性テンソルは物体系で時間不変となり、オイラー方程式が適用可能になる

剛体の仮定は、回転運動の複雑性を管理可能にするための抽象化 なのである。実在の物体は多少変形するが、その変形が運動に比べて十分小さければ、剛体近似は極めて有効である。

次章への展望

本章で明らかにしたのは、回転運動を記述するには並進運動とは質的に異なる概念（慣性テンソル、角運動量、トルク）が必要であり、その運動方程式は非線形で複雑だということである。

しかし、この複雑性は剛体という抽象化によって管理可能になる。次章では、剛体の定義を厳密に与え、剛体の配置を記述する数学的枠組み（姿勢の表現）を構築する。それによって、一つの剛体の運動を完全に記述する能力を獲得する。

その後、複数の剛体が連結されたシステム（人体など）へと進むが、その際には本章で学んだ回転運動の複雑性が、システム全体のレベルで再び現れることになる。しかし、一つ一つの剛体については本章の理論で記述できるため、問題は「どのように連結されているか」（拘束条件）と「それをどう定式化するか」に絞られる。これが、後の章での主題となる。

まとめ：本章で獲得した概念

本章では以下の重要な概念を導入し、回転運動の本質的な複雑性を明らかにした：

角速度ベクトル $\boldsymbol{\omega}$ ：回転の速さと軸を同時に表現する量
慣性テンソル $\mathbf{I}$ ：質量分布が回転運動に及ぼす影響をテンソル量として表現
主軸と主慣性モーメント ：慣性テンソルを対角化する特別な座標系
角運動量 $\boldsymbol{L} = \mathbf{I}\boldsymbol{\omega}$ ：回転の「勢い」を表す

変形可能から剛体への抽象化

本章の目的と問題意識

これまで私たちは、質点の並進運動と回転運動を個別に扱ってきました。しかし現実の物体は、質点のような「大きさのない点」ではありません。木製の椅子、金属のロボットアーム、そして人体の骨格など、実際の物体は 空間的な広がり を持ち、無数の原子や分子から構成されています。

ここで根本的な問いが生じます：

なぜ私たちは、無数の粒子から構成される複雑な物体を、たった数個の変数で記述できるのか？

椅子が倒れる運動を記述する際、私たちは椅子を構成する $10^{23}$ 個もの原子それぞれの運動を追跡する必要はありません。わずか6つの変数（位置3つ、姿勢3つ）で椅子全体の運動が記述できます。この驚くべき簡略化を可能にする概念が 剛体（rigid body） です。

本章では、この抽象化のプロセスを丁寧に追跡します。具体的には：

なぜ変形を無視できるのか という物理的根拠
剛体とは何か という数学的・物理的定義
剛体の仮定が何を可能にするか という帰結

この理解なくして、複雑な多体システムの運動方程式を正しく導出することはできません。

現実の物体が持つ複雑性：変形の普遍性

すべての物体は変形する

物理的真実を直視しましょう： 完全に変形しない物体は自然界に存在しません。

ゴム製のボールを握れば、大きく変形します
鋼鉄の梁も、重量物を載せれば微小にたわみます
ダイヤモンドでさえ、十分な力を加えれば変形します
人体の骨は歩行時に微小に曲がります

変形の度合いは材料の 弾性率 （stiffness）に依存しますが、原理的にはすべての物体が力に応じて形状を変えます。

変形を考慮した場合の記述の困難性

もし物体の変形を完全に記述しようとすれば、私たちは 連続体力学 （continuum mechanics）の領域に入ります。この枠組みでは：

物体内の各点 $\mathbf{r}$ における変位 $\mathbf{u}(\mathbf{r}, t)$ を関数として記述
応力テンソル $\sigma_{ij}(\mathbf{r}, t)$ とひずみテンソル $\epsilon_{ij}(\mathbf{r}, t)$ の関係を記述
偏微分方程式系を解く必要がある

例えば、弾性体の運動方程式は次のような形になります：

\rho \frac{\partial^2 \mathbf{u}}{\partial t^2} = \nabla \cdot \boldsymbol{\sigma} + \mathbf{f}

ここで $\rho$ は密度、 $\boldsymbol{\sigma}$ は応力テンソル、 $\mathbf{f}$ は体積力です。

この記述は数学的には美しいですが、 人体の歩行のような複雑な運動を解析する という私たちの目的には、あまりにも複雑すぎます。変形の詳細まで追跡すれば、問題は解析不可能なほど複雑になります。

剛体という抽象化：本質的な運動のみを抽出する

何を捨て、何を残すか

ここで戦略的な問いを立てます：

変形を無視しても、物体の「大域的な運動」は正しく記述できるだろうか？

この問いへの答えが、 剛体近似 （rigid body approximation）です。

剛体近似の核心は：

捨てるもの ：物体内部の相対的な位置変化（変形）
残すもの ：物体全体の並進運動と回転運動

この近似が有効なのは、 変形のタイムスケールやエネルギースケールが、全体運動のそれと比較して十分小さい 場合です。

具体例：歩行する人体の場合

人間が歩行する場合を考えましょう。

変形に関わるスケール：

大腿骨の微小なたわみ： $10^{-6}$ m オーダー
関節軟骨の変形： $10^{-4}$ m オーダー
これらの変形による弾性エネルギー：全運動エネルギーの $< 0.1\%$

全体運動に関わるスケール：

重心の移動距離：1 m オーダー（1歩あたり）
脚の回転角度：数十度オーダー
運動エネルギーの主要部分：並進と回転

このスケール分離があるため、「大腿骨は剛体である」「脛骨は剛体である」と近似しても、歩行運動の本質的な特徴（歩幅、周期、関節トルクなど）は正確に捉えられます。

剛体の厳密な定義

定義1：距離不変性による定義

剛体を最も明快に定義する方法は、 物体内の任意の2点間の距離が時間によらず一定 という条件です。

定義： 物体 $$B$$ が剛体であるとは、物体内の任意の2点 $$P, Q$$ に対して、それらの間の距離 $|\mathbf{r}_P - \mathbf{r}_Q|$ が時刻 $$t$$ によらず一定であることをいう。

|\mathbf{r}_P(t) - \mathbf{r}_Q(t)| = \text{const.} \quad \forall P, Q \in B, \, \forall t

この定義は直感的でありながら、数学的に厳密です。

定義の帰結：自由度の劇的な削減

この定義から、極めて重要な帰結が導かれます。

定理： 3次元空間における剛体の配置は、 6個の独立なパラメータ で完全に記述される。

証明の概略：

物体内の任意の1点 $$O$$ の位置を指定するのに3つのパラメータが必要（並進の自由度）
点 $$O$$ を固定した状態で、物体の向きを指定するのに3つのパラメータが必要（回転の自由度）
距離不変性の条件により、これ以上の自由度は存在しない

したがって、 $$N$$ 個の質点から構成される物体は本来 $$3N$$ 個の自由度を持ちますが、剛体では $$(3N - 6)$$ 個の拘束条件が課され、実質的な自由度は6個に削減されます。

具体例： アボガドロ数 $N_A \approx 10^{23}$ 個の原子からなる物体でも、剛体近似すれば6変数で記述可能です。これが剛体力学の驚異的な単純化の源泉です。

剛体配置の数学的記述

剛体の配置を記述する標準的な方法を確立しましょう。

参照配置と現在配置

剛体の状態を記述するため、2つの配置を導入します：

参照配置（reference configuration） ：ある基準時刻 $$t_0$$ における物体の配置
現在配置（current configuration） ：時刻 $$t$$ における物体の配置

記述に必要な要素

剛体の現在配置を指定するには：

（1）基準点の位置ベクトル

物体内に固定された基準点 $$O'$$ （通常は重心を選ぶ）の、空間固定座標系（慣性系）における位置ベクトル：

\mathbf{R}(t) = (X(t), Y(t), Z(t))

これが 並進自由度 （3個）を表します。

（2）姿勢を表す回転

参照配置から現在配置への回転を表す 回転行列 （rotation matrix） $\mathbf{R}(t)$ ：

\mathbf{R}(t) \in SO(3)

ここで $$SO(3)$$ は3次元特殊直交群（special orthogonal group）です。これが 回転自由度 （3個）を表します。

物体内の任意点の位置

物体内の任意の点 $$P$$ の位置は、次のように記述されます：

\mathbf{r}_P(t) = \mathbf{R}(t) + \mathbf{R}(t) \mathbf{r}'_P

ここで：

$\mathbf{r}_P(t)$ ：慣性系における点 $$P$$ の位置
$\mathbf{R}(t)$ ：基準点 $$O'$$ の位置（並進）
$\mathbf{R}(t)$ ：回転行列
$\mathbf{r}'_P$ ：物体固定座標系における点 $$P$$ の位置（時間不変）

重要な注意： $\mathbf{r}'_P$ は物体固定座標系で測った座標であり、 時間によって変化しません。これが剛体の定義の数学的表現です。

剛体近似の物理的正当化

いつ剛体近似は有効か

剛体近似が正当化される条件を明確にしましょう。

条件1：変形の空間スケールが小さい

物体の代表的なサイズを $$L$$ とし、最大変形量を $\delta$ としたとき：

\frac{\delta}{L} \ll 1

例：大腿骨（ $L \sim 0.5$ m）のたわみ（ $\delta \sim 10^{-6}$ m）では $\delta/L \sim 10^{-6}$

条件2：変形のエネルギーが支配的でない

弾性変形エネルギー $E_{\text{elastic}}$ と運動エネルギー $E_{\text{kinetic}}$ の比：

\frac{E_{\text{elastic}}}{E_{\text{kinetic}}} \ll 1

条件3：変形の時間スケールが速い

物体の固有振動周期 $T_{\text{vib}}$ と、注目する運動の周期 $T_{\text{motion}}$ の比：

\frac{T_{\text{vib}}}{T_{\text{motion}}} \ll 1

これは、変形の応答が瞬時に平衡状態に達し、運動に追従することを意味します。

剛体近似が破綻する例

逆に、剛体近似が 使えない 状況も理解することが重要です：

高速衝突 ：衝突時間が材料の応力波伝播時間と同程度の場合
共振現象 ：外力の周波数が物体の固有振動数に近い場合
座屈・大変形 ：圧縮荷重下での細長い棒の座屈など
柔軟構造 ：ロープ、布、柔軟ロボットアームなど

剛体としての質量分布の記述

剛体近似により配置の記述は簡略化されましたが、物体の 質量分布 は依然として重要です。

連続的質量分布

剛体内の質量分布は、密度関数 $\rho(\mathbf{r}')$ によって記述されます。ここで $\mathbf{r}'$ は物体固定座標系での位置です。

重要な点： 剛体では $\rho(\mathbf{r}')$ は時間によらず一定です（質量は再配分されない）。

全質量

剛体の全質量 $$M$$ は：

M = \int_V \rho(\mathbf{r}') \, dV'

ここで積分は物体の体積 $$V$$ 全体にわたります。

重心の定義

重心（center of mass）の位置 $\mathbf{r}'_G$ は、物体固定座標系で：

\mathbf{r}'_G = \frac{1}{M} \int_V \mathbf{r}' \rho(\mathbf{r}') \, dV'

慣習： 通常、物体固定座標系の原点を重心に選びます。この場合 $\mathbf{r}'_G = \mathbf{0}$ となります。

慣性モーメントテンソルの導入

回転運動を記述するには、質量分布の 回転に対する抵抗 を定量化する必要があります。これが 慣性モーメントテンソル （moment of inertia tensor）です。

物体固定座標系における慣性モーメントテンソル $\mathbf{I}'$ の成分は：

I'_{ij} = \int_V \rho(\mathbf{r}') \left( |\mathbf{r}'|^2 \delta_{ij} - r'_i r'_j \right) dV'

ここで $\delta_{ij}$ はクロネッカーのデルタです。

重要な性質：

$\mathbf{I}'$ は対称行列： $I'_{ij} = I'_{ji}$
$\mathbf{I}'$ は正定値（固有値はすべて正）
剛体では $\mathbf{I}'$ は 時間不変 （物体固定座標系で測る限り）

この最後の性質が、剛体の回転運動を扱う上で決定的に重要です。変形する物体では、慣性モーメントテンソルが時々刻々と変化し、運動方程式が著しく複雑になります。

剛体近似がもたらす数学的恩恵

座標系の選択の自由

剛体では、物体内に 固定された座標系 （body-fixed frame）を自由に設定できます。この座標系では：

すべての点の座標 $\mathbf{r}'$ が時間不変
質量分布 $\rho(\mathbf{r}')$ が時間不変
慣性モーメントテンソル $\mathbf{I}'$ が時間不変

これにより、複雑な回転運動を 物体固定座標系 と 空間固定座標系 を使い分けて記述する強力な手法が使えるようになります。

運動方程式の分離

剛体の運動方程式は、次のように 並進と回転に分離 できます：

並進運動（重心の運動）：

M \ddot{\mathbf{R}}_G = \mathbf{F}_{\text{ext}}

回転運動（重心まわりの回転）：

\frac{d\mathbf{L}_G}{dt} = \boldsymbol{\tau}_{\text{ext}}

ここで $\mathbf{L}_G$ は重心まわりの角運動量、 $\boldsymbol{\tau}_{\text{ext}}$ は重心まわりの外力モーメントです。

この分離は、 ケーニッヒの定理 （König's theorem）と呼ばれ、剛体力学の基礎となります。変形する物体では、このような明快な分離は一般には成立しません。

離散質点モデルと連続体モデルの統一

2つの視点の等価性

剛体を記述する際、2つのアプローチがあります：

（1）離散質点モデル

物体を $$N$$ 個の質点の集合とみなし、それらの間の距離を固定：

|\mathbf{r}_i - \mathbf{r}_j| = c_{ij} \quad (\text{const.})

（2）連続体モデル

物体を連続的な質量分布とみなし、密度関数 $\rho(\mathbf{r}')$ で記述。

これらは本質的に同等であり、連続体モデルは離散モデルの極限 $N \to \infty$ と理解できます。

実用上の使い分け

理論的導出 ：連続体モデルを用いることが多い（積分による表現）
数値計算 ：離散モデルを用いることが多い（有限要素法など）
多体システム ：各剛体を1つの単位として扱い、接続する

人体モデルでは、各骨を1つの剛体（離散単位）として扱い、それらを関節で接続するアプローチが標準的です。

剛体の内部応力の役割

内部応力の存在と不可視性

剛体内部には、依然として 内部応力 （internal stress）が存在します。剛体近似は「変形が小さい」ことを仮定するのであって、「内部応力が存在しない」ことを意味するわけではありません。

例えば、静止している剛体棒を両端から引っ張る場合：

棒内部には引張応力が発生
しかし伸びは無視できるほど小さい
剛体近似では、この微小な伸びを無視

運動方程式における内部力の相殺

剛体全体の運動方程式を導出する際、 内部力は対になって相殺 されます。これはニュートンの第3法則（作用・反作用の法則）の帰結です。

物体内の質点 $$i, j$$ 間の内部力を $\mathbf{f}_{ij}, \mathbf{f}_{ji}$ とすると：

\mathbf{f}_{ij} = -\mathbf{f}_{ji}

全質点について和をとれば：

\sum_{i,j} \mathbf{f}_{ij} = 0

したがって、 剛体全体の運動は外力のみで決定 され、内部応力の詳細を知る必要はありません。これが剛体力学のもう1つの大きな簡略化です。

抽象化のレベルの階層性

ここまでの議論を整理すると、物体の記述には複数の抽象化レベルがあることがわかります：

レベル1：量子力学的記述

個々の原子・電子の波動関数
自由度：無限大（場の理論）

レベル2：古典粒子系

$$N$$ 個の粒子の位置・速度
自由度： $$6N$$

レベル3：連続体力学

変位場 $\mathbf{u}(\mathbf{r}, t)$
自由度：無限大（連続関数）

レベル4：剛体力学

位置 $\mathbf{R}(t)$ と姿勢 $\mathbf{R}(t)$
自由度：6

レベル5：質点力学

重心位置 $\mathbf{R}_G(t)$ のみ
自由度：3

私たちの目標である「人体歩行の記述」には、 レベル4（剛体力学） が適切です。各骨を剛体として扱い、関節で接続することで、歩行運動の本質を捉えつつ、計算可能な複雑さに抑えることができます。

本章のまとめ

本章では、現実の変形可能な物体から剛体への抽象化について、その物理的根拠と数学的定式化を詳述しました。

主要な概念

すべての物体は変形するが、多くの場合その変形は全体運動に比べて無視できる
剛体は「物体内の任意の2点間距離が一定」という拘束で定義される
この拘束により、自由度が $$3N$$ から6に劇的に削減される
剛体の配置は、並進（3自由度）と回転（3自由度）で完全に記述される
物体固定座標系では質量分布が時間不変となり、計算が大幅に簡略化される

次章への展望

剛体という概念の確立により、私たちは以下を扱う準備が整いました：

単一剛体の運動方程式の導出（並進と回転の結合）
複数の剛体が接続されたシステム（多体システム）
関節による拘束条件の定式化
最終的な目標：人体のような複雑な構造の運動記述

次章では、複数の物体（剛体）が互いに力を及ぼし合う状況を、どのように体系的に記述するかを学びます。これが多体システムへの第一歩となります。

複数物体間の相互作用の記述枠組み

6.1 問題意識:単一物体から複数物体へ

これまで私たちは単一の剛体の運動を記述してきました。しかし、人体の歩行を記述するという目標を考えると、現実は遥かに複雑です。人体は大腿骨、脛骨、足部など、多数の剛体が関節で接続された 多体系(multi-body system) です。

ここで重要な問いが生じます: 複数の物体が存在するとき、それぞれの運動方程式をどのように記述し、相互にどう関連づけるのか?

単純に考えれば、物体Aと物体Bがあるとき、それぞれに対してニュートン・オイラー方程式を書けばよいように思えます:

\begin{align} m_A \ddot{\boldsymbol{r}}_A &= \boldsymbol{F}_A \\ \boldsymbol{I}_A \dot{\boldsymbol{\omega}}_A + \boldsymbol{\omega}_A \times (\boldsymbol{I}_A \boldsymbol{\omega}_A) &= \boldsymbol{\tau}_A \\ m_B \ddot{\boldsymbol{r}}_B &= \boldsymbol{F}_B \\ \boldsymbol{I}_B \dot{\boldsymbol{\omega}}_B + \boldsymbol{\omega}_B \times (\boldsymbol{I}_B \boldsymbol{\omega}_B) &= \boldsymbol{\tau}_B \end{align}

しかし、この素朴な定式化には重大な問題があります。 物体間に力が作用する場合、 $\boldsymbol{F}_A$ と $\boldsymbol{F}_B$ は独立ではありません。AがBを押せば、BもAを押し返します(作用・反作用の法則)。この相互作用をどう扱うかが、複数物体系の記述における中心的課題です。

本章の目標

作用・反作用の法則を数学的に正確に表現する
内力と外力の概念的区別を明確にする
接触力のモデル化手法を理解する
複数物体系の運動方程式の系統的な構築法を習得する

6.2 内力と外力の厳密な区別

6.2.1 概念的定義

物体系を扱う際、力を 外力(external force) と 内力(internal force) に区別することが本質的に重要です。この区別は単なる便宜的分類ではなく、運動方程式の構造そのものに深く関わります。

外力とは、 考察している系の外部から作用する力 です。例えば:

重力(地球という外部物体からの力)
地面からの反力
風による空気抵抗

内力とは、 系内の物体間で相互に作用する力 です。例えば:

関節で接続された二つの骨の間に働く力
連結された振り子の各部分間の張力
接触している物体間の接触力

6.2.2 内力の根本的性質:作用・反作用の法則

ニュートンの第三法則(作用・反作用の法則)は、内力に対して極めて特殊な構造を与えます。物体Aが物体Bに力 $\boldsymbol{F}_{AB}$ を及ぼすとき、物体Bは物体Aに力 $\boldsymbol{F}_{BA}$ を及ぼし、これらは次の関係を満たします:

\boldsymbol{F}_{BA} = -\boldsymbol{F}_{AB}

この関係は 力の大きさが等しく方向が逆 という以上の深い意味を持ちます。数学的には、系全体の運動量保存則の根拠となっています。

物体Aの質量中心の運動方程式が:

m_A \ddot{\boldsymbol{r}}_A = \boldsymbol{F}_A^{\text{ext}} + \boldsymbol{F}_{BA}

物体Bの質量中心の運動方程式が:

m_B \ddot{\boldsymbol{r}}_B = \boldsymbol{F}_B^{\text{ext}} + \boldsymbol{F}_{AB}

であるとき、これらを足し合わせると:

m_A \ddot{\boldsymbol{r}}_A + m_B \ddot{\boldsymbol{r}}_B = \boldsymbol{F}_A^{\text{ext}} + \boldsymbol{F}_B^{\text{ext}} + (\boldsymbol{F}_{BA} + \boldsymbol{F}_{AB})

作用・反作用の法則により $\boldsymbol{F}_{BA} + \boldsymbol{F}_{AB} = \boldsymbol{0}$ ですから:

\frac{d}{dt}(m_A \dot{\boldsymbol{r}}_A + m_B \dot{\boldsymbol{r}}_B) = \boldsymbol{F}_A^{\text{ext}} + \boldsymbol{F}_B^{\text{ext}}

系全体の運動量変化は外力の総和のみで決まり、内力は寄与しません。この事実は、多体系を扱う際の基本原理です。

6.2.3 トルクに関する作用・反作用

力だけでなく、トルク(モーメント)についても作用・反作用を考える必要があります。しかしここで注意すべき微妙な点があります。

物体Aが物体Bに及ぼす力 $\boldsymbol{F}_{AB}$ が作用点 $\boldsymbol{r}_P$ で働くとき、物体Bの質量中心周りのトルクは:

\boldsymbol{\tau}_B = (\boldsymbol{r}_P - \boldsymbol{r}_B) \times \boldsymbol{F}_{AB}

一方、物体Aが受けるトルクは:

\boldsymbol{\tau}_A = (\boldsymbol{r}_P - \boldsymbol{r}_A) \times \boldsymbol{F}_{BA} = -(\boldsymbol{r}_P - \boldsymbol{r}_A) \times \boldsymbol{F}_{AB}

ここで重要なのは、 $\boldsymbol{\tau}_A + \boldsymbol{\tau}_B \neq \boldsymbol{0}$ となりうる ということです:

\boldsymbol{\tau}_A + \boldsymbol{\tau}_B = \{(\boldsymbol{r}_A - \boldsymbol{r}_B) \times \boldsymbol{F}_{AB}\}

この量がゼロになるのは、 $\boldsymbol{F}_{AB}$ が二つの質量中心を結ぶ線上にある場合のみです。一般には系全体の角運動量は内力によっても変化しうるのです。

ただし、 力とトルクの対(force-couple) として内力を考える場合、物体Aが物体Bに力 $\boldsymbol{F}_{AB}$ とトルク $\boldsymbol{M}_{AB}$ を及ぼすなら:

\boldsymbol{F}_{BA} = -\boldsymbol{F}_{AB}, \quad \boldsymbol{M}_{BA} = -\boldsymbol{M}_{AB}

という形で作用・反作用が成立します。この形式は関節力のモデル化で頻繁に使われます。

6.3 接触力のモデル化

6.3.1 接触の物理的本質と抽象化の必要性

物体同士が「接触する」とは、微視的には複雑な電磁気的相互作用です。しかし剛体力学では、この複雑さを抽象化し、 接触点での力とトルク として記述します。

この抽象化には二つのレベルがあります:

レベル1: 分布力から集中力へ

実際の接触では、力は有限の面積にわたって分布しています。しかし剛体の仮定の下では、この分布を 接触点での単一の力ベクトル と トルクベクトル に集約できます。これは第5章で学んだ「剛体では力の作用点を移動できる」という性質の応用です。

レベル2: 力の成分分解

接触力 $\boldsymbol{F}_{\text{contact}}$ は、接触面の法線方向成分(垂直抗力)と接線方向成分(摩擦力)に分解されます:

\boldsymbol{F}_{\text{contact}} = \boldsymbol{F}_N + \boldsymbol{F}_f

ここで:

$\boldsymbol{F}_N = F_N \boldsymbol{n}$ (法線ベクトル $\boldsymbol{n}$ 方向の垂直抗力)
$\boldsymbol{F}_f$ は接触面内の摩擦力

6.3.2 垂直抗力の特性

垂直抗力には基本的な 不等式拘束 があります:

F_N \geq 0

これは「引っ張る接触はできない」という物理的事実を表します。物体は押し合うことはできても、接着剤がなければ引っ張り合うことはできません。

この不等式性は、接触状態の 離散的な変化 を引き起こします:

$$F_N > 0$$ のとき:接触継続中
$$F_N = 0$$ になった瞬間:離脱が起こる

人体の歩行モデルでは、足が地面から離れる瞬間がこの状態に対応します。運動方程式の形式そのものが変わるため、この遷移の扱いは後の章で詳しく学びます。

6.3.3 摩擦力のモデル化

摩擦力のモデル化には複数のレベルがあり、目的に応じて選択します。

クーロン摩擦モデル

最も標準的なモデルで、次の性質を持ちます:

静止摩擦 : 物体が相対静止しているとき( $\boldsymbol{v}_{\text{rel}} = \boldsymbol{0}$ )、摩擦力は運動を引き起こそうとする力を打ち消す方向に働き、その大きさには上限があります:

\|\boldsymbol{F}_f\| \leq \mu_s F_N

ここで $\mu_s$ は静止摩擦係数です。

動摩擦 : 物体が滑っているとき( $\boldsymbol{v}_{\text{rel}} \neq \boldsymbol{0}$ )、摩擦力は相対速度と逆向きで、大きさは垂直抗力に比例します:

\boldsymbol{F}_f = -\mu_k F_N \frac{\boldsymbol{v}_{\text{rel}}}{\|\boldsymbol{v}_{\text{rel}}\|}

ここで $\mu_k$ は動摩擦係数で、一般に $\mu_k < \mu_s$ です。

モデルの数学的課題

このモデルには重要な数学的特性があります:

静止と滑りの境界( $\|\boldsymbol{F}_f\| = \mu_s F_N$ )で摩擦力が不連続になる
相対速度ゼロでの摩擦力が一意に定まらない(集合値になる)

これらは数値計算上の困難を引き起こしますが、物理的には正しい記述です。実用的には、 正則化(regularization) と呼ばれる手法で滑らかな関数に近似することもあります:

\boldsymbol{F}_f = -\mu F_N \frac{\boldsymbol{v}_{\text{rel}}}{\sqrt{\|\boldsymbol{v}_{\text{rel}}\|^2 + \epsilon^2}}

ここで $\epsilon$ は小さな正定数です。 $\|\boldsymbol{v}_{\text{rel}}\| \gg \epsilon$ のとき動摩擦に近づき、 $\|\boldsymbol{v}_{\text{rel}}\| \ll \epsilon$ のとき粘性的な性質を示します。

6.3.4 ヒンジ関節(回転関節)の力学モデル

人体の関節の多くは ヒンジ関節(hinge joint) として近似できます。これは一つの軸周りの回転のみを許し、他の自由度を拘束する接続です。

ヒンジ関節の力学的記述には二つのアプローチがあります:

アプローチ1: 拘束力による記述

関節点で物体AとBが接続されているとき、関節が及ぼす力 $\boldsymbol{F}_{\text{joint}}$ とトルク $\boldsymbol{M}_{\text{joint}}$ を未知数として扱います。

物体Aの運動方程式:

\begin{align} m_A \ddot{\boldsymbol{r}}_A &= \boldsymbol{F}_A^{\text{ext}} + \boldsymbol{F}_{\text{joint},A} \\ \boldsymbol{I}_A \dot{\boldsymbol{\omega}}_A + \boldsymbol{\omega}_A \times (\boldsymbol{I}_A \boldsymbol{\omega}_A) &= \boldsymbol{\tau}_A^{\text{ext}} + \boldsymbol{M}_{\text{joint},A} + (\boldsymbol{r}_{\text{joint}} - \boldsymbol{r}_A) \times \boldsymbol{F}_{\text{joint},A} \end{align}

物体Bの運動方程式:

\begin{align} m_B \ddot{\boldsymbol{r}}_B &= \boldsymbol{F}_B^{\text{ext}} + \boldsymbol{F}_{\text{joint},B} \\ \boldsymbol{I}_B \dot{\boldsymbol{\omega}}_B + \boldsymbol{\omega}_B \times (\boldsymbol{I}_B \boldsymbol{\omega}_B) &= \boldsymbol{\tau}_B^{\text{ext}} + \boldsymbol{M}_{\text{joint},B} + (\boldsymbol{r}_{\text{joint}} - \boldsymbol{r}_B) \times \boldsymbol{F}_{\text{joint},B} \end{align}

ここで作用・反作用の法則により:

\boldsymbol{F}_{\text{joint},B} = -\boldsymbol{F}_{\text{joint},A}, \quad \boldsymbol{M}_{\text{joint},B} = -\boldsymbol{M}_{\text{joint},A}

さらに、ヒンジ関節の幾何学的拘束条件として、回転軸を $\boldsymbol{e}_{\text{axis}}$ とすると:

\begin{align} \boldsymbol{r}_{\text{joint},A} &= \boldsymbol{r}_{\text{joint},B} \quad \text{(位置の一致)} \\ (\boldsymbol{\omega}_A - \boldsymbol{\omega}_B) \cdot \boldsymbol{e}_{\perp} &= 0 \quad \text{(回転軸以外の回転の禁止)} \end{align}

この拘束条件と運動方程式を連立させることで、拘束力と運動を同時に決定します。

アプローチ2: 縮約座標系による記述

拘束条件を満たす座標系(後述の一般化座標)を最初から選ぶことで、拘束力を陽に扱わない定式化も可能です。これは次章以降で詳しく学びます。

6.4 複数物体系の運動方程式の系統的構築

6.4.1 記述の基本フレームワーク

$$N$$ 個の剛体からなる系を考えます。各物体 $$i$$ ( $i = 1, 2, \ldots, N$ )に対して、次の量を定義します:

質量中心位置: $\boldsymbol{r}_i$
姿勢(回転行列): $\boldsymbol{R}_i$ または四元数 $\boldsymbol{q}_i$
角速度: $\boldsymbol{\omega}_i$
質量: $$m_i$$
慣性テンソル: $\boldsymbol{I}_i$

各物体の運動は、ニュートン・オイラー方程式で記述されます:

\begin{align} m_i \ddot{\boldsymbol{r}}_i &= \boldsymbol{F}_i^{\text{total}} \\ \boldsymbol{I}_i \dot{\boldsymbol{\omega}}_i + \boldsymbol{\omega}_i \times (\boldsymbol{I}_i \boldsymbol{\omega}_i) &= \boldsymbol{\tau}_i^{\text{total}} \end{align}

ここで $\boldsymbol{F}_i^{\text{total}}$ と $\boldsymbol{\tau}_i^{\text{total}}$ は物体 $$i$$ に働く すべての力とトルクの合計 です。

6.4.2 力の系統的分類と記法

$\boldsymbol{F}_i^{\text{total}}$ を次のように分解します:

\boldsymbol{F}_i^{\text{total}} = \boldsymbol{F}_i^{\text{ext}} + \sum_{j \neq i} \boldsymbol{F}_{ji}

ここで:

$\boldsymbol{F}_i^{\text{ext}}$ : 物体 $$i$$ に働く外力(重力など)
$\boldsymbol{F}_{ji}$ : 物体 $$j$$ が物体 $$i$$ に及ぼす内力

記法の注意: $\boldsymbol{F}_{ji}$ は「 $$j$$ から $$i$$ へ」という順序で、作用・反作用の法則は:

\boldsymbol{F}_{ji} = -\boldsymbol{F}_{ij}

となります。

トルクについても同様に:

\boldsymbol{\tau}_i^{\text{total}} = \boldsymbol{\tau}_i^{\text{ext}} + \sum_{j \neq i} \left[ (\boldsymbol{r}_{ji} - \boldsymbol{r}_i) \times \boldsymbol{F}_{ji} + \boldsymbol{M}_{ji} \right]

ここで:

$\boldsymbol{r}_{ji}$ : 物体 $$j$$ が物体 $$i$$ に力を及ぼす作用点
$\boldsymbol{M}_{ji}$ : 物体 $$j$$ が物体 $$i$$ に及ぼす純トルク(force-couple)

6.4.3 接続関係の表現:グラフ構造

複数物体系の接続関係は、数学的には グラフ(graph) として表現できます:

ノード(node) : 各剛体
エッジ(edge) : 物体間の相互作用(関節、接触など)

例えば、人体の下肢を考えると:

ノード: 大腿(thigh)、下腿(shank)、足部(foot)
エッジ: 膝関節(大腿-下腿)、足首関節(下腿-足部)、地面接触(足部-地面)

このグラフ構造は、運動方程式を効率的に構築するための 隣接行列(adjacency matrix) や 接続リスト として実装されます。

物体 $$i$$ と物体 $$j$$ が接続されているとき、接続行列 $$C$$ の要素を:

C_{ij} = \begin{cases} 1 & \text{物体 $i$ と $j$ が接続されている} \\ 0 & \text{それ以外} \end{cases}

と定義すれば、内力の和は:

\sum_{j \neq i} \boldsymbol{F}_{ji} = \sum_{j=1}^{N} C_{ij} \boldsymbol{F}_{ji}

と書けます。この形式は計算機実装で有用です。

6.4.4 構築の具体的手順

複雑な多体系の運動方程式を構築する際の系統的手順を示します:

ステップ1: 系の自由度の特定

各剛体は本来6自由度(位置3 + 姿勢3)を持ちます。 $$N$$ 個の剛体があれば、総自由度は $$6N$$ です。

ステップ2: 接続関係の列挙

すべての相互作用をリスト化します:

関節の種類と接続する物体のペア
接触の可能性がある物体のペア

ステップ3: 各物体の自由物体図(Free Body Diagram)の作成

各物体に働く力を以下のように分類して図示します:

重力
各接続点からの力とトルク
外部からの作用力

ステップ4: 各物体の運動方程式を記述

物体 $$i$$ に対して:

\begin{align} m_i \ddot{\boldsymbol{r}}_i &= m_i \boldsymbol{g} + \sum_{j \in \mathcal{N}(i)} \boldsymbol{F}_{ji} + \boldsymbol{F}_i^{\text{external}} \\ \boldsymbol{I}_i \dot{\boldsymbol{\omega}}_i + \boldsymbol{\omega}_i \times (\boldsymbol{I}_i \boldsymbol{\omega}_i) &= \sum_{j \in \mathcal{N}(i)} \left[ (\boldsymbol{r}_{ji} - \boldsymbol{r}_i) \times \boldsymbol{F}_{ji} + \boldsymbol{M}_{ji} \right] + \boldsymbol{\tau}_i^{\text{external}} \end{align}

ここで $\mathcal{N}(i)$ は物体 $$i$$ と接続されている物体の集合です。

ステップ5: 作用・反作用の法則を適用

すべての内力のペアに対して:

\boldsymbol{F}_{ji} = -\boldsymbol{F}_{ij}, \quad \boldsymbol{M}_{ji} = -\boldsymbol{M}_{ij}

ステップ6: 接続の物理モデルを導入

各接続に対して適切なモデル(6.3節)を選択し、未知の内力を既知の状態変数で表現するか、拘束条件として定式化します。

ステップ7: 方程式の整理

最終的に、状態変数を $\boldsymbol{x} = (\boldsymbol{r}_1, \boldsymbol{R}_1, \boldsymbol{\omega}_1, \ldots, \boldsymbol{r}_N, \boldsymbol{R}_N, \boldsymbol{\omega}_N)^T$ とすると、次の形式の方程式系が得られます:

\boldsymbol{M}(\boldsymbol{x}) \ddot{\boldsymbol{x}} = \boldsymbol{f}(\boldsymbol{x}, \dot{\boldsymbol{x}}, t)

ここで $\boldsymbol{M}$ は質量行列、 $\boldsymbol{f}$ は力の項です。拘束がある場合は次章で学ぶより精緻な扱いが必要になります。

6.5 具体例:二重振り子系

理論を具体化するため、 二重振り子(double pendulum) の運動方程式を詳細に構築します。この系は単純ながら、多体系の本質的な特徴を持っています。

6.5.1 系の設定

物理的構成 :

物体1(長さ $$l_1$$ 、質量 $$m_1$$ ):固定点Oにヒンジで接続
物体2(長さ $$l_2$$ 、質量 $$m_2$$ ):物体1の端点にヒンジで接続
重力: $\boldsymbol{g} = -g \boldsymbol{e}_z$
回転軸: すべて $$z$$ 軸(平面運動)

座標の選択 :

物体1の質量中心: $\boldsymbol{r}_1 = (x_1, y_1, 0)^T$
物体2の質量中心: $\boldsymbol{r}_2 = (x_2, y_2, 0)^T$
物体1の角度: $\theta_1$ (鉛直下向きからの角度)
物体2の角度: $\theta_2$ (鉛直下向きからの角度)

ただし、この系は拘束されているため、実際の自由度は2(角度 $\theta_1, \theta_2$ )です。まず拘束力を陽に含む形式で記述します。

6.5.2 幾何学的関係

固定点を原点とすると:

\begin{align} \boldsymbol{r}_1 &= \frac{l_1}{2}(\sin\theta_1, -\cos\theta_1, 0)^T \\ \boldsymbol{r}_{\text{joint}} &= l_1(\sin\theta_1, -\cos\theta_1, 0)^T \\ \boldsymbol{r}_2 &= l_1(\sin\theta_1, -\cos\theta_1, 0)^T + \frac{l_2}{2}(\sin\theta_2, -\cos\theta_2, 0)^T \end{align}

角速度(平面運動なので $$z$$ 成分のみ):

\boldsymbol{\omega}_1 = \dot{\theta}_1 \boldsymbol{e}_z, \quad \boldsymbol{\omega}_2 = \dot{\theta}_2 \boldsymbol{e}_z

6.5.3 各物体の自由物体図

物体1に働く力 :

重力: $m_1 \boldsymbol{g}$
固定点からの拘束力: $\boldsymbol{F}_O$ (未知)
関節点から物体2が及ぼす力: $\boldsymbol{F}_{21}$ (未知)

物体2に働く力 :

重力: $m_2 \boldsymbol{g}$
関節点から物体1が及ぼす力: $\boldsymbol{F}_{12} = -\boldsymbol{F}_{21}$ (作用・反作用)

6.5.4 運動方程式の記述

物体1の並進運動 :

m_1 \ddot{\boldsymbol{r}}_1 = m_1 \boldsymbol{g} + \boldsymbol{F}_O + \boldsymbol{F}_{21}

**物

拘束条件下での運動記述

本章の問題意識と背景課題

これまで私たちは、自由に空間を運動する剛体、あるいは複数の剛体が相互作用する系を扱ってきました。しかし現実の機械システムや生物の運動では、物体は完全に自由に動けるわけではありません。

典型的な例を考えてみましょう：

扉はヒンジ（蝶番）によって、ある軸周りの回転しかできません
車輪は地面に接触しながら転がり、滑らないという条件下で運動します
人体の膝関節は、大腿骨と脛骨が離れないという条件下で曲げ伸ばしされます
ロボットアームの各関節は、特定の自由度のみを持つように設計されています

これらの運動は「 拘束条件 （constraint）」の下での運動です。拘束条件とは、物体の位置や速度が満たさなければならない数学的な制約のことです。

なぜ拘束条件の扱いが重要なのか

人体の歩行を記述する際、各関節は骨同士が離れないという拘束を持ちます。もし拘束を無視して運動方程式を立てると、私たちは以下の困難に直面します：

自由度の過剰設定 ：実際には動けない方向にも運動の可能性を考慮してしまう
物理的に不可能な解 ：骨が離れる、物体が貫通するなどの非現実的な運動が数学的に許されてしまう
拘束力の未知性 ：ヒンジが及ぼす力、関節が及ぼす力など、拘束によって生じる力が未知のままになる

本章の目標は、 拘束条件を数学的に定式化し、拘束下での運動方程式を導出する体系的な方法 を習得することです。これにより、人体や機械などの複雑な多体系の運動を正確に記述できるようになります。

拘束条件の分類と数学的表現

拘束条件とは何か：本質的理解

拘束条件を理解する上で、まず明確にすべき核心的な問いがあります： 「自由度が減る」とは数学的に何を意味するのか？

3次元空間における単一剛体は、並進3自由度と回転3自由度の計6自由度を持ちます。この「自由度」とは、 系の配置を指定するために必要な独立なパラメータの数 を意味します。

拘束条件が加わると、これらのパラメータ間に関係式が生じ、独立なパラメータの数が減少します。この関係式を数学的に表現したものが 拘束方程式 です。

ホロノミック拘束と非ホロノミック拘束

拘束条件は、その数学的性質によって2つの根本的に異なるタイプに分類されます：

ホロノミック拘束（holonomic constraint）

定義：座標のみの関数として表現できる拘束

ホロノミック拘束は次の形で表されます：

\Phi(\mathbf{q}, t) = 0

ここで、 $\mathbf{q}$ は一般化座標のベクトル、 $$t$$ は時間です。

重要な特徴 ：この拘束は座標の値そのものに制約を課すもので、速度には直接的な制約を課しません（ただし座標の制約から間接的に速度にも制約が生じます）。

具体例1：固定長の棒

質点Aと質点Bが長さ $$l$$ の剛体棒で結ばれている場合：

\Phi = |\mathbf{r}_B - \mathbf{r}_A|^2 - l^2 = (x_B - x_A)^2 + (y_B - y_A)^2 + (z_B - z_A)^2 - l^2 = 0

この拘束は2つの質点の位置座標のみで表現されており、典型的なホロノミック拘束です。

具体例2：曲面上を運動する質点

質点が半径 $$R$$ の球面上を運動する場合：

\Phi = x^2 + y^2 + z^2 - R^2 = 0

この拘束により、3次元空間の3自由度が2自由度（球面上の2次元）に減少します。

非ホロノミック拘束（non-holonomic constraint）

定義：速度に関する制約として表現され、座標のみの関数に積分できない拘束

非ホロノミック拘束は次の形で表されます：

\sum_i a_i(\mathbf{q}, t) \dot{q}_i + a_t(\mathbf{q}, t) = 0

この関係式が $\Phi(\mathbf{q}, t) = 0$ の形に積分できないとき、これを非ホロノミック拘束と呼びます。

具体例：滑らない車輪の転がり

地面を滑らずに転がる車輪を考えます。車輪の中心位置を $$(x, y)$$ 、車輪の向きを $\theta$ 、回転角を $\phi$ とします。

「滑らない」という条件は：

\dot{x} = R \dot{\phi} \cos\theta, \quad \dot{y} = R \dot{\phi} \sin\theta

これらは速度間の関係式です。この関係式を積分しても、 $x, y, \theta, \phi$ の間の代数的関係式（座標のみの関数）は得られません。なぜなら、車輪がどのような経路を通ったかによって、最終的な $$(x, y)$$ と $(\theta, \phi)$ の関係が異なるからです。

なぜこの分類が重要なのか

この2つの拘束の区別は、単なる数学的な分類ではありません。 運動方程式の導出方法と、系の性質が根本的に異なります ：

ホロノミック拘束 ：拘束を満たす座標系（一般化座標）を選ぶことで、拘束を陽に扱わずに済ませることができる
非ホロノミック拘束 ：拘束を満たす座標系への変換ができず、常に拘束を陽に扱う必要がある

人体の関節の多くはホロノミック拘束ですが、足と地面の接触（滑らない条件）は非ホロノミック拘束です。したがって歩行の記述には両方の扱いが必要です。

拘束力の概念と物理的意味

拘束力とは何か

拘束条件が存在するとき、その拘束を維持するために何らかの力が作用しています。この力を 拘束力 （constraint force）と呼びます。

具体的イメージ ：

棒で結ばれた2つの質点：棒が及ぼす張力または圧縮力が拘束力
球面上を運動する質点：球面からの垂直抗力が拘束力
ヒンジで固定された扉：ヒンジが及ぼす反力が拘束力
人体の関節：靭帯や関節面が及ぼす力が拘束力

拘束力の本質的性質

拘束力について、以下の重要な性質を理解する必要があります：

性質1：拘束力は未知である

通常の力（重力、バネの力など）は、物体の位置や速度の既知の関数として表現できます。しかし 拘束力の大きさは事前には分からず、運動の結果として決まります。

なぜでしょうか？拘束力は「拘束を維持するために必要なだけの大きさ」を持つからです。例えば、振り子の糸の張力は、振り子がどの位置にあるか、どのくらいの速さで運動しているかによって変化し、その値は運動方程式を解いて初めて決まります。

性質2：理想的拘束と仮想仕事の原理

多くの実用的な問題では、拘束力は 理想的拘束 （ideal constraint）であると仮定します。理想的拘束とは：

定義：拘束力が許される微小変位（仮想変位）に対して仕事をしない拘束

数学的には：

\sum_i \mathbf{F}_i^{(c)} \cdot \delta\mathbf{r}_i = 0

ここで、 $\mathbf{F}_i^{(c)}$ は質点 $$i$$ に作用する拘束力、 $\delta\mathbf{r}_i$ は拘束条件を満たす仮想変位です。

物理的意味 ：

滑らかな面上の運動：摩擦がないので、垂直抗力は面に垂直、変位は面に沿うため、仕事はゼロ
固定長の棒：棒の張力は棒の方向、許される変位は棒の方向に垂直なので、仕事はゼロ
ヒンジ：ヒンジの反力は並進を拘束する力で、許される変位は回転のみなので、仕事はゼロ

非理想的拘束の例 ：摩擦力は拘束力の一種ですが、変位に対して仕事をするため、理想的拘束ではありません。

仮想変位と実変位の違い：微妙だが決定的な概念

ここで、解析力学の最も微妙で重要な概念の一つを明確にします： 仮想変位 （virtual displacement）と 実変位 （actual displacement）の区別です。

実変位

実変位 $d\mathbf{r}$ は、実際の運動における無限小時間 $$dt$$ の間の位置の変化です：

d\mathbf{r} = \mathbf{v} \, dt

実変位は運動方程式に従い、時間の経過を伴います。

仮想変位

仮想変位 $\delta\mathbf{r}$ は、 時間を止めた瞬間に、拘束条件を満たしながら系の配置を微小に変化させる想像上の変位 です。

重要な性質 ：

時間を固定する： $\delta t = 0$
拘束条件を満たす： $\Phi(\mathbf{q} + \delta\mathbf{q}, t) = 0$ （一次の微小量まで）
複数の可能性がある：仮想変位は一般に多数存在する（拘束がない限り）

なぜこの区別が重要か ：

実変位は唯一つに決まりますが（運動方程式の解として）、仮想変位は拘束を満たす範囲で任意です。理想的拘束の定義「拘束力が仮想変位に対して仕事をしない」とは、「拘束を満たすあらゆる可能な変位に対して仕事をしない」という強い条件を意味します。

具体例 ：球面上の質点

質点が球面 $$x^2 + y^2 + z^2 = R^2$$ 上にあるとき、拘束条件は：

x \, \delta x + y \, \delta y + z \, \delta z = 0

これは、仮想変位 $(\delta x, \delta y, \delta z)$ が位置ベクトルに垂直、つまり球面の接平面内にあることを意味します。

垂直抗力 $\mathbf{N} = N(x, y, z)/R$ は半径方向なので：

\mathbf{N} \cdot \delta\mathbf{r} = \frac{N}{R}(x \, \delta x + y \, \delta y + z \, \delta z) = 0

これが「理想的拘束」の意味です。

ホロノミック拘束系の運動方程式：ラグランジュの未定乗数法

問題設定と基本戦略

$$N$$ 個の質点からなる系を考えます。拘束がない場合、系は $$3N$$ 個の自由度を持ちます。ここに $$m$$ 個のホロノミック拘束：

\Phi_\alpha(\mathbf{r}_1, \ldots, \mathbf{r}_N, t) = 0, \quad \alpha = 1, 2, \ldots, m

が課されているとします。

各質点 $$i$$ には、外力 $\mathbf{F}_i$ と拘束力 $\mathbf{F}_i^{(c)}$ が作用し、運動方程式は：

m_i \ddot{\mathbf{r}}_i = \mathbf{F}_i + \mathbf{F}_i^{(c)}

困難点 ：拘束力 $\mathbf{F}_i^{(c)}$ が未知です。

解決戦略の2つのアプローチ ：

拘束力を含めて解く ：拘束力を未知変数として扱い、運動方程式と拘束条件を連立して解く
拘束力を消去する ：拘束を満たす座標系（一般化座標）を選び、拘束力が現れない形の運動方程式を導く

まず第1のアプローチ、 ラグランジュの未定乗数法 を学びます。

ラグランジュの未定乗数法の導出

ステップ1：拘束条件の微分形式

ホロノミック拘束 $\Phi_\alpha = 0$ を時間で微分すると：

\frac{d\Phi_\alpha}{dt} = \sum_i \frac{\partial \Phi_\alpha}{\partial \mathbf{r}_i} \cdot \dot{\mathbf{r}}_i + \frac{\partial \Phi_\alpha}{\partial t} = 0

さらに微分して加速度の関係式を得ます（詳細は省略しますが、 $\ddot{\mathbf{r}}_i$ に関する線形関係式になります）。

ステップ2：仮想変位と拘束条件

時刻 $$t$$ を固定して拘束条件を考えると、仮想変位は：

\sum_i \frac{\partial \Phi_\alpha}{\partial \mathbf{r}_i} \cdot \delta\mathbf{r}_i = 0, \quad \alpha = 1, 2, \ldots, m

を満たします。これは、 $\delta\mathbf{r}_i$ が $$m$$ 個の線形制約を満たすことを意味します。

ステップ3：ダランベールの原理

運動方程式 $m_i \ddot{\mathbf{r}}_i = \mathbf{F}_i + \mathbf{F}_i^{(c)}$ を変形すると：

(m_i \ddot{\mathbf{r}}_i - \mathbf{F}_i) - \mathbf{F}_i^{(c)} = 0

両辺と仮想変位の内積をとり、全質点について和をとります：

\sum_i (m_i \ddot{\mathbf{r}}_i - \mathbf{F}_i) \cdot \delta\mathbf{r}_i - \sum_i \mathbf{F}_i^{(c)} \cdot \delta\mathbf{r}_i = 0

理想的拘束では第2項がゼロなので：

\sum_i (m_i \ddot{\mathbf{r}}_i - \mathbf{F}_i) \cdot \delta\mathbf{r}_i = 0

これを ダランベールの原理 （d'Alembert's principle）と呼びます。

ステップ4：ラグランジュの未定乗数の導入

ここで核心的なアイデアが登場します。仮想変位 $\delta\mathbf{r}_i$ は独立ではなく、 $$m$$ 個の拘束条件：

\sum_i \frac{\partial \Phi_\alpha}{\partial \mathbf{r}_i} \cdot \delta\mathbf{r}_i = 0

を満たします。

数学的な補題：ベクトル $\mathbf{A}_i$ ( $i = 1, \ldots, N$ ) について、 $$m$$ 個の線形制約を満たす任意の $\delta\mathbf{r}_i$ に対して $\sum_i \mathbf{A}_i \cdot \delta\mathbf{r}_i = 0$ が成り立つならば、適当な係数 $\lambda_\alpha$ が存在して：

\mathbf{A}_i = \sum_{\alpha=1}^{m} \lambda_\alpha \frac{\partial \Phi_\alpha}{\partial \mathbf{r}_i}

と表される。

この $\lambda_\alpha$ を ラグランジュの未定乗数 （Lagrange multiplier）と呼びます。

ダランベールの原理に適用すると：

m_i \ddot{\mathbf{r}}_i - \mathbf{F}_i = \sum_{\alpha=1}^{m} \lambda_\alpha \frac{\partial \Phi_\alpha}{\partial \mathbf{r}_i}

ステップ5：拘束力の表現

上式を運動方程式と比較すると：

\mathbf{F}_i^{(c)} = \sum_{\alpha=1}^{m} \lambda_\alpha \frac{\partial \Phi_\alpha}{\partial \mathbf{r}_i}

これは 拘束力が拘束条件の勾配の線形結合として表される ことを示しています。

幾何学的解釈 ： $\nabla_{\mathbf{r}_i} \Phi_\alpha$ は拘束面に垂直なベクトルです。したがって拘束力は拘束面に垂直な方向にのみ作用します（理想的拘束の帰結）。

ラグランジュの運動方程式（拘束系）

結局、拘束系の運動方程式は：

m_i \ddot{\mathbf{r}}_i = \mathbf{F}_i + \sum_{\alpha=1}^{m} \lambda_\alpha \frac{\partial \Phi_\alpha}{\partial \mathbf{r}_i}, \quad i = 1, \ldots, N

と拘束条件：

\Phi_\alpha(\mathbf{r}_1, \ldots, \mathbf{r}_N, t) = 0, \quad \alpha = 1, \ldots, m

の連立系として表されます。

未知変数 ： $$3N$$ 個の座標成分 $\mathbf{r}_i$ と $$m$$ 個の未定乗数 $\lambda_\alpha$ 、計 $$3N + m$$ 個

方程式 ： $$3N$$ 個の運動方程式と $$m$$ 個の拘束条件、計 $$3N + m$$ 個

したがって、系は数学的に閉じています（未知数と方程式の数が一致）。

具体例：2質点を結ぶ剛体棒

系の設定 ：

質量 $$m_1, m_2$$ の2つの質点が、長さ $$l$$ の質量のない剛体棒で結ばれています。重力下での運動を考えます。

拘束条件 ：

\Phi = |\mathbf{r}_2 - \mathbf{r}_1|^2 - l^2 = 0

勾配：

\frac{\partial \Phi}{\partial \mathbf{r}_1} = -2(\mathbf{r}_2 - \mathbf{r}_1), \quad \frac{\partial \Phi}{\partial \mathbf{r}_2} = 2(\mathbf{r}_2 - \mathbf{r}_1)

運動方程式 ：

m_1 \ddot{\mathbf{r}}_1 = m_1 \mathbf{g} - 2\lambda(\mathbf{r}_2 - \mathbf{r}_1)

m_2 \ddot{\mathbf{r}}_2 = m_2 \mathbf{g} + 2\lambda(\mathbf{r}_2 - \mathbf{r}_1)

物理的解釈 ：

$2\lambda(\mathbf{r}_2 - \mathbf{r}_1)$ の項は、棒に沿った方向の力です。質点1には引っ張られる方向（または押される方向）、質点2には逆向きの力が作用します。これが棒による拘束力（張力または圧縮力）です。

$\lambda$ の符号によって：

$\lambda > 0$ ：棒は張力を及ぼす（引っ張り）
$\lambda < 0$ ：棒は圧縮力を及ぼす（押し）

運動の解析 ：

重心 $\mathbf{R} = (m_1\mathbf{r}_1 + m_2\mathbf{r}_2)/(m_1 + m_2)$ と相対位置 $\mathbf{r} = \mathbf{r}_2 - \mathbf{r}_1$ を導入すると、運動方程式は分離されます。重心は自由落下し、相対運動は拘束 $|\mathbf{r}| = l$ の下での回転運動になります。

一般化座標による拘束の消去

一般化座標の導入：なぜ拘束力を扱わなくて済むのか

ラグランジュの未定乗数法は、拘束力を求めるには有効ですが、運動そのものを求めるには未定乗数という余分な変数を扱う必要があります。

より洗練されたアプローチ ：拘束を自動的に満たす座標系を選べば、拘束力は運動方程式に現れなくなります。

基本的アイデア

$$3N$$ 自由度の系に $$m$$ 個の独立なホロノミック拘束があれば、実質的な自由度は $$n = 3N - m$$ です。この $$n$$ 個の独立なパラメータ $q_1, q_2, \ldots, q_n$ を選び、すべての座標をこれらで表現します：

\mathbf{r}_i = \mathbf{r}_i(q_1, \ldots, q_n, t), \quad i = 1, \ldots, N

この $$q_j$$ を 一般化座標 （generalized coordinates）と呼びます。

重要な性質 ：この表現は拘束条件を自動的に満たします。なぜなら、拘束を考慮して座標変換を構成したからです。

ラグランジュ方程式の導出

一般化座標 $$q_j$$ を用いると、運動方程式は美しい統一的な形式になります：

\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}_j}\right) - \frac{\partial L}{\partial q_j} = 0, \quad j = 1, \ldots, n

ここで $$L = T - V$$ は ラグランジアン （Lagrangian）、 $$T$$ は運動エネルギー、 $$V$$ はポテンシャルエネルギーです。

導出の概略（詳細は複雑なので要点のみ）

速度を一般化座標で表現：

\dot{\mathbf{r}}_i = \sum_j \frac{\partial \mathbf{r}_i}{\partial q_j}\dot{q}_j + \frac{\partial \mathbf{r}_i}{\partial t}

運動エネルギーは $\dot{q}_j$ の2次形式になります
ダランベールの原理を一般化座標で表現し、 $\delta q_j$ が独立であることを用いると、上記の形が導かれます

この方程式の驚くべき点 ：

拘束力が現れない（理想的拘束の場合）
座標系の選び方によらず同じ形式（共変性）
運動エネルギーとポテンシャルエネルギーのみから導ける

具体例：単振り子

系の設定 ：

質量 $$m$$ の質点が、長さ $$l$$ の質量のない糸で吊るされています。糸は常に張った状態を保ちます（ホロノミック拘束）。

デカルト座標での拘束 ：

\Phi = x^2 + (y - l)^2 - l^2 = 0

（原点を吊り点とし、 $$y$$ 軸を鉛直下向きにとる）

一般化座標の選択 ：

鉛直下向きからの角度 $\theta$ を一般化座標として選びます：

x = l\sin\theta, \quad y = l(1 - \cos\theta)

この変換は拘束条件を自動的に満たします（ $\theta$ の値に関わらず質点は糸の先端にあります）。

速度：

\dot{x} = l\cos\theta \cdot \dot{\theta}, \quad \dot{y} = l\sin\theta \cdot \dot{\theta}

運動エネルギー ：

T = \frac{1}{2}m(\dot{x}^2 + \dot{y}^2) = \frac{1}{2}m l^2(\cos^2\theta + \sin^2\theta)\dot{\theta}^2 = \frac{1}{2}m l^2 \dot{\theta}^2

ポテンシャルエネルギー （吊り点を基準）：

V = -mgy = -mgl(1 - \cos\theta) = -mgl + mgl\cos\theta

定数項は無視できるので：

V = mgl\cos\theta

ラグランジアン ：

L = T - V = \frac{1}{2}m l^2 \dot{\theta}^2 - mgl\cos\theta

ラグランジュ方程式 ：

\frac{\partial L}{\partial \dot{\theta}} = ml^2\dot{\theta}

\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{\theta}}\right) = ml^2\ddot{\theta}

\frac{\partial L}{\partial \theta} = mgl\sin\theta

したがって：

ml^2\ddot{\theta} - mgl\sin\theta = 0

\ddot{\theta} = -\frac{g}{l}\sin\theta

これが単振り子の運動方程式です。拘束力（糸の張力）は現れ

多体システムへの一般化

8.1 多体システム記述の問題意識

8.1.1 なぜ一般化が必要か

これまでの章で、私たちは単一剛体の運動、二つの物体間の拘束、そして小規模な連結システムを扱ってきました。しかし、 人体の歩行 のような現実の複雑なシステムを記述するには、決定的な飛躍が必要です。

人体は約200個の骨から構成され、主要な運動に関与する関節だけでも数十に及びます。このような系を「個別に考える」アプローチでは限界があります。なぜなら:

スケーラビリティの欠如 : 物体が増えるたびに、記述方法を一から考え直すのは非現実的
体系性の不在 : 個別アプローチでは、異なるシステム間で知見を転用できない
計算実装の困難 : コンピュータでシミュレーションする際、統一的な枠組みがなければプログラム化できない

したがって、私たちは 任意の数の剛体から構成されるシステム を記述できる一般的な枠組みを必要としています。

8.1.2 多体システムの本質的特徴

多体システムには、単一剛体や2体系にはない本質的な特徴があります:

累積効果 : 一つの物体の運動が連鎖的に他の物体に影響を及ぼす。例えば、歩行中の腕の振りは、肩、胴体、骨盤、脚へと影響が伝播します。

配位の爆発 : $$n$$ 個の剛体が存在すると、理論上は $$6n$$ 個の自由度が存在します。しかし、拘束により実効的な自由度はこれより遥かに小さくなります。この「見かけの複雑さ」と「本質的な複雑さ」の乖離を正しく扱う必要があります。

階層構造 : 人体のような系には、「胴体から腕が生える」「腕から前腕が生える」という階層的な構造があります。この構造を活用できれば、記述が大幅に簡略化されます。

8.2 多体系の構成要素と記法

8.2.1 索引付けによる体系化

$$n$$ 個の剛体から成る系を考えます。各剛体に番号 $i = 1, 2, \ldots, n$ を割り当てます。

剛体 $$i$$ について:

質量 : $$m_i$$
慣性テンソル : $\mathbf{I}_i$ （物体固定座標系で記述）
位置 : 質量中心の位置ベクトル $\mathbf{r}_i$
姿勢 : 回転行列 $\mathbf{R}_i$ または四元数 $\mathbf{q}_i$
並進速度 : $\mathbf{v}_i = \dot{\mathbf{r}}_i$
角速度 : $\boldsymbol{\omega}_i$ （空間座標系で記述）

この記法により、任意の多体系を 統一的に表現 できます。

8.2.2 状態ベクトルの導入

システム全体の状態を一つのベクトルにまとめると便利です。これを 状態ベクトル と呼びます。

位置・姿勢を表す 配位ベクトル $\mathbf{q}$ :

\mathbf{q} = \begin{bmatrix} \mathbf{r}_1 \\ \boldsymbol{\theta}_1 \\ \mathbf{r}_2 \\ \boldsymbol{\theta}_2 \\ \vdots \\ \mathbf{r}_n \\ \boldsymbol{\theta}_n \end{bmatrix}

ここで $\boldsymbol{\theta}_i$ は剛体 $$i$$ の姿勢を表すパラメータ（オイラー角など）です。

速度を表す 速度ベクトル $\dot{\mathbf{q}}$ または $\mathbf{u}$ :

\mathbf{u} = \begin{bmatrix} \mathbf{v}_1 \\ \boldsymbol{\omega}_1 \\ \mathbf{v}_2 \\ \boldsymbol{\omega}_2 \\ \vdots \\ \mathbf{v}_n \\ \boldsymbol{\omega}_n \end{bmatrix}

重要な注意 : 姿勢パラメータ $\boldsymbol{\theta}_i$ の時間微分 $\dot{\boldsymbol{\theta}}_i$ は、一般に角速度 $\boldsymbol{\omega}_i$ とは 異なります （オイラー角の場合など）。したがって、 $\mathbf{u}$ と $\dot{\mathbf{q}}$ は単純に等しくない場合があります。この関係は 運動学的関係式 $\mathbf{u} = \mathbf{B}(\mathbf{q})\dot{\mathbf{q}}$ で表されます。

8.3 多体系の運動エネルギー

8.3.1 総運動エネルギーの構成

多体系の運動エネルギーは、各剛体の運動エネルギーの和です:

T = \sum_{i=1}^{n} T_i = \sum_{i=1}^{n} \left[ \frac{1}{2}m_i \mathbf{v}_i \cdot \mathbf{v}_i + \frac{1}{2}\boldsymbol{\omega}_i \cdot (\mathbf{I}_i^{\text{world}} \boldsymbol{\omega}_i) \right]

ここで $\mathbf{I}_i^{\text{world}}$ は空間座標系で表現された慣性テンソルです:

\mathbf{I}_i^{\text{world}} = \mathbf{R}_i \mathbf{I}_i \mathbf{R}_i^T

8.3.2 質量行列による表現

運動エネルギーを速度ベクトル $\mathbf{u}$ で表現すると:

T = \frac{1}{2}\mathbf{u}^T \mathbf{M}(\mathbf{q}) \mathbf{u}

ここで $\mathbf{M}(\mathbf{q})$ は 質量行列 （または慣性行列）と呼ばれる $6n \times 6n$ の対称正定値行列です。

質量行列のブロック構造:

\mathbf{M}(\mathbf{q}) = \begin{bmatrix} \mathbf{M}_{11} & \mathbf{M}_{12} & \cdots & \mathbf{M}_{1n} \\ \mathbf{M}_{21} & \mathbf{M}_{22} & \cdots & \mathbf{M}_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ \mathbf{M}_{n1} & \mathbf{M}_{n2} & \cdots & \mathbf{M}_{nn} \end{bmatrix}

対角ブロック $\mathbf{M}_{ii}$ は:

\mathbf{M}_{ii} = \begin{bmatrix} m_i \mathbf{I}_{3} & \mathbf{0} \\ \mathbf{0} & \mathbf{I}_i^{\text{world}} \end{bmatrix}

拘束がない場合 （すべての剛体が独立）、非対角ブロック $\mathbf{M}_{ij}$ ( $i \neq j$ ) はすべてゼロです。

拘束がある場合 （剛体が連結している場合）、非対角ブロックが 非ゼロ になることがあります。これは物体間の 動力学的結合 を表します。

8.3.3 配位依存性の理解

質量行列 $\mathbf{M}(\mathbf{q})$ が配位 $\mathbf{q}$ に依存する理由は、 慣性テンソルの姿勢依存性 です:

\mathbf{I}_i^{\text{world}} = \mathbf{R}_i(\mathbf{q}) \mathbf{I}_i \mathbf{R}_i(\mathbf{q})^T

姿勢 $\mathbf{R}_i$ が変われば、空間座標系での慣性テンソルも変わります。これが質量行列の配位依存性の本質です。

物理的意味 : 同じ角速度を与えても、物体の向きが違えば必要なトルクが異なる、という現象を数学的に表現しています。

8.4 拘束を持つ多体系の運動方程式

8.4.1 拘束なし系の運動方程式

拘束がない $$n$$ 個の剛体系の運動方程式は:

\mathbf{M}(\mathbf{q})\dot{\mathbf{u}} + \mathbf{C}(\mathbf{q}, \mathbf{u})\mathbf{u} = \mathbf{F}^{\text{ext}}

ここで:

$\mathbf{M}(\mathbf{q})$ : 質量行列（ $6n \times 6n$ ）
$\mathbf{C}(\mathbf{q}, \mathbf{u})$ : コリオリ・遠心力項を表す行列
$\mathbf{F}^{\text{ext}}$ : 外力ベクトル（重力、接触力など）

8.4.2 コリオリ・遠心力項の起源

$\mathbf{C}(\mathbf{q}, \mathbf{u})\mathbf{u}$ は、 配位の変化 に伴う運動エネルギーの非線形効果を表します。

その成分を具体的に書くと:

[\mathbf{C}(\mathbf{q}, \mathbf{u})\mathbf{u}]_k = \sum_{i,j} C_{kij}(\mathbf{q}) u_i u_j

ここで:

C_{kij} = \frac{1}{2}\left( \frac{\partial M_{ki}}{\partial q_j} + \frac{\partial M_{kj}}{\partial q_i} - \frac{\partial M_{ij}}{\partial q_k} \right)

これは クリストッフェル記号 と呼ばれる微分幾何学の概念と同等です。

物理的解釈 :

回転する系での遠心力・コリオリ力
質量分布が変化する系での見かけの力
連結された物体の「振り回し効果」

8.4.3 拘束力の導入

拘束 $\boldsymbol{\Phi}(\mathbf{q}) = \mathbf{0}$ を持つ系では、拘束力 $\mathbf{F}^{\text{const}}$ が追加されます:

\mathbf{M}(\mathbf{q})\dot{\mathbf{u}} + \mathbf{C}(\mathbf{q}, \mathbf{u})\mathbf{u} = \mathbf{F}^{\text{ext}} + \mathbf{F}^{\text{const}}

拘束力は ラグランジュの未定乗数 $\boldsymbol{\lambda}$ を用いて:

\mathbf{F}^{\text{const}} = \mathbf{J}_{\boldsymbol{\Phi}}^T(\mathbf{q}) \boldsymbol{\lambda}

ここで $\mathbf{J}_{\boldsymbol{\Phi}}$ は拘束のヤコビ行列です:

\mathbf{J}_{\boldsymbol{\Phi}} = \frac{\partial \boldsymbol{\Phi}}{\partial \mathbf{q}}

8.4.4 拘束系の完全な方程式系

拘束を持つ多体系は、以下の 微分代数方程式系 （DAE: Differential-Algebraic Equations）で記述されます:

\begin{cases} \mathbf{M}(\mathbf{q})\dot{\mathbf{u}} + \mathbf{C}(\mathbf{q}, \mathbf{u})\mathbf{u} = \mathbf{F}^{\text{ext}} + \mathbf{J}_{\boldsymbol{\Phi}}^T(\mathbf{q}) \boldsymbol{\lambda} \\ \boldsymbol{\Phi}(\mathbf{q}) = \mathbf{0} \\ \dot{\mathbf{q}} = \mathbf{B}^{-1}(\mathbf{q})\mathbf{u} \end{cases}

これは $$6n + m + 6n$$ 個の方程式（ $$m$$ は拘束数）から成り、未知数は:

$\mathbf{u}$ : $$6n$$ 個
$\boldsymbol{\lambda}$ : $$m$$ 個
$\mathbf{q}$ : $$6n$$ 個

この方程式系が多体系記述の標準形です。

8.5 階層的多体系の効率的記述

8.5.1 木構造とループ構造

多体系の構造は、グラフ理論の観点から分類できます:

木構造（Tree Structure） :

閉じたループがない連結構造
例: 人体の骨格（開脚時）、ロボットアーム
特徴: ルート（根）から各物体への経路が一意

ループ構造（Loop Structure） :

閉じたループを含む構造
例: 四足歩行動物の両足接地状態、パラレルリンク機構
特徴: 冗長な拘束を持つ

木構造は 数値的に扱いやすく、効率的なアルゴリズムが存在します。

8.5.2 相対座標と絶対座標

多体系を記述する座標系には二つの主要なアプローチがあります:

絶対座標（Absolute Coordinates） :

各剛体の位置・姿勢を基準座標系から直接記述
これまで使ってきた $\mathbf{q} = [\mathbf{r}_1, \boldsymbol{\theta}_1, \ldots, \mathbf{r}_n, \boldsymbol{\theta}_n]^T$
利点: 概念が単純、任意の構造に適用可能
欠点: 拘束方程式が必要、自由度が見かけ上多い

相対座標（Relative Coordinates / Joint Coordinates） :

隣接する物体間の相対的な位置・姿勢を記述
例: 関節角度ベクトル $\boldsymbol{\theta} = [\theta_1, \theta_2, \ldots, \theta_m]^T$
利点: 拘束が暗黙的に満たされる、自由度が最小
欠点: 運動方程式の導出が複雑

どちらを使うべきか :

木構造系で効率重視 → 相対座標
一般的な系、数値安定性重視 → 絶対座標
人体モデルの初期段階 → 相対座標が自然

8.5.3 再帰的運動方程式の構造

木構造系では、 再帰的アルゴリズム により効率的に運動方程式を構築できます。

基本的な考え方:

外向き伝播（Forward Kinematics） : ルートから末端へ、速度・加速度を伝播
内向き伝播（Backward Dynamics） : 末端からルートへ、力・トルクを伝播

この二段階プロセスにより、 $$O(n)$$ の計算量で運動方程式を解けます（一般的な方法は $$O(n^3)$$ ）。

代表的なアルゴリズム:

Featherstone's Articulated Body Algorithm : 最も効率的な $$O(n)$$ アルゴリズム
Composite Rigid Body Algorithm : 質量行列を効率的に構築

8.5.4 関節の数学的記述

相対座標アプローチでは、関節が中心的な役割を果たします。

関節 $$j$$ が剛体 $$i$$ と剛体 $$k$$ を連結する場合:

関節の自由度 : 関節が許す相対運動の次元数

回転関節（hinge）: 1自由度
球関節（ball joint）: 3自由度
平面関節: 3自由度

関節の運動学 : 親リンク $$i$$ から子リンク $$k$$ への変換

\mathbf{T}_k = \mathbf{T}_i \cdot \mathbf{T}_j(\theta_j)

ここで $\mathbf{T}_j(\theta_j)$ は関節 $$j$$ の変換行列、 $\theta_j$ は関節変数です。

関節速度の伝播 :

\mathbf{v}_k = \mathbf{v}_i + \boldsymbol{\omega}_i \times \mathbf{r}_{ik} + \mathbf{v}_j^{\text{rel}}

\boldsymbol{\omega}_k = \boldsymbol{\omega}_i + \boldsymbol{\omega}_j^{\text{rel}}

ここで $\mathbf{r}_{ik}$ は物体 $$i$$ から物体 $$k$$ の質量中心への位置ベクトル、 $\mathbf{v}_j^{\text{rel}}$ と $\boldsymbol{\omega}_j^{\text{rel}}$ は関節が生成する相対速度です。

8.6 多体系運動方程式の導出法

8.6.1 ラグランジュ法による導出

最も体系的な方法は ラグランジュの運動方程式 です。

ラグランジアン $$L = T - V$$ を定義し:

\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}_j}\right) - \frac{\partial L}{\partial q_j} = Q_j

ここで $$Q_j$$ は一般化力、 $$q_j$$ は一般化座標（配位ベクトル $\mathbf{q}$ の成分）です。

多体系への適用 :

運動エネルギー $T(\mathbf{q}, \dot{\mathbf{q}})$ を計算
位置エネルギー $V(\mathbf{q})$ を計算
各座標成分に対してラグランジュ方程式を適用

この方法の利点:

座標系の選択に依存しない
拘束を考慮しやすい （拘束付きラグランジュ方程式）
エネルギー保存則との対応が明確

8.6.2 ニュートン・オイラー法による導出

各剛体に対して個別にニュートン・オイラー方程式を適用する方法:

剛体 $$i$$ に対して:

\begin{cases} m_i \dot{\mathbf{v}}_i = \mathbf{F}_i^{\text{ext}} + \sum_{j} \mathbf{F}_{ji} \\ \mathbf{I}_i^{\text{world}} \dot{\boldsymbol{\omega}}_i + \boldsymbol{\omega}_i \times (\mathbf{I}_i^{\text{world}} \boldsymbol{\omega}_i) = \boldsymbol{\tau}_i^{\text{ext}} + \sum_{j} \boldsymbol{\tau}_{ji} \end{cases}

ここで $\mathbf{F}_{ji}$ , $\boldsymbol{\tau}_{ji}$ は剛体 $$j$$ から剛体 $$i$$ への拘束力・トルクです。

作用反作用の法則 :

\mathbf{F}_{ij} = -\mathbf{F}_{ji}, \quad \boldsymbol{\tau}_{ij} = -\boldsymbol{\tau}_{ji}

この方法の利点:

物理的意味が直感的
力の伝達が明示的
再帰的アルゴリズムと相性が良い

8.6.3 仮想仕事の原理による導出

仮想変位 $\delta \mathbf{q}$ に対する仮想仕事:

\delta W = \sum_{i=1}^{n} (\mathbf{F}_i \cdot \delta \mathbf{r}_i + \boldsymbol{\tau}_i \cdot \delta \boldsymbol{\theta}_i)

D'Alembertの原理により、慣性力を含めて:

\delta W - \sum_{i=1}^{n} (m_i \ddot{\mathbf{r}}_i \cdot \delta \mathbf{r}_i + [\mathbf{I}_i \dot{\boldsymbol{\omega}}_i + \boldsymbol{\omega}_i \times \mathbf{I}_i \boldsymbol{\omega}_i] \cdot \delta \boldsymbol{\theta}_i) = 0

拘束条件下では、仮想変位は拘束と矛盾しない方向のみに制限されます:

\mathbf{J}_{\boldsymbol{\Phi}} \delta \mathbf{q} = \mathbf{0}

この原理から拘束付き運動方程式が自然に導かれます。

8.6.4 三つの方法の比較と使い分け

方法	適用場面	利点	欠点
ラグランジュ法	一般的な理論導出、エネルギー保存系	体系的、座標系不変	計算量が多い
ニュートン・オイラー法	木構造系、数値シミュレーション	効率的、直感的	拘束の扱いが煩雑
仮想仕事の原理	複雑な拘束を持つ系	拘束力を陽に扱わない	抽象的

実用的な選択基準 :

理論的理解や手計算 → ラグランジュ法
効率的な数値計算（木構造） → ニュートン・オイラー法（再帰的）
複雑な拘束を持つ一般的な系 → 仮想仕事 + 数値的拘束解法

8.7 多体系の数値シミュレーション

8.7.1 時間積分の課題

多体系の運動方程式を数値的に解く際、いくつかの特有の課題があります:

剛性（Stiffness） : 大きく異なる時間スケールが共存する

例: 剛体の高速振動と全体のゆっくりした運動
解決策: 陰的積分法（implicit methods）

拘束のドリフト : 数値誤差により拘束条件 $\boldsymbol{\Phi}(\mathbf{q}) = \mathbf{0}$ が徐々に破れる

放置すると系が崩壊する
解決策: 拘束安定化（constraint stabilization）

エネルギー保存 : 長時間シミュレーションでエネルギーが人為的に増減する

解決策: シンプレクティック積分法（symplectic integrators）

8.7.2 拘束安定化手法

拘束 $\boldsymbol{\Phi}(\mathbf{q}) = \mathbf{0}$ とその速度レベル $\dot{\boldsymbol{\Phi}} = \mathbf{J}_{\boldsymbol{\Phi}} \dot{\mathbf{q}} = \mathbf{0}$ が数値誤差で破れるのを防ぐ手法:

Baumgarte安定化 :
加速度レベルの拘束方程式に安定化項を追加:

\ddot{\boldsymbol{\Phi}} + 2\alpha \dot{\boldsymbol{\Phi}} + \beta^2 \boldsymbol{\Phi} = \mathbf{0}

ここで $\alpha, \beta > 0$ は安定化パラメータ。これは拘束違反を フィードバック制御 で抑制する考え方です。

射影法（Projection Method） :
各時間ステップ後、配位を拘束多様体上に射影:

\mathbf{q}^{\text{corrected}} = \mathbf{q} - \mathbf{J}_{\boldsymbol{\Phi}}^+ \boldsymbol{\Phi}(\mathbf{q})

ここで $\mathbf{J}_{\boldsymbol{\Phi}}^+$ は擬似逆行列です。

8.7.3 代表的な時間積分スキーム

陽的オイラー法（Explicit Euler） :

\begin{cases} \mathbf{q}^{n+1} = \mathbf{q}^n + h \dot{\mathbf{q}}^n \\ \mathbf{u}^{n+1} = \mathbf{u}^n + h \mathbf{M}^{-1}(\mathbf{F} - \mathbf{C}\mathbf{u})^n \end{cases}

簡単だが数値的に不安定、時間刻み $$h$$ を小さくする必要

半陰的オイラー法（Semi-implicit Euler） :

\begin{cases} \mathbf{u}^{n+1} = \mathbf{u}^n + h \mathbf{M}^{-1}(\mathbf{F} - \mathbf{C}\mathbf{u})^n \\ \mathbf{q}^{n+1} = \mathbf{q}^n + h \dot{\mathbf{q}}^{n+1} \end{cases}

速度を先に更新し、その結果で位置を更新
エネルギー散逸が自然に起こり安定

Runge-Kutta法（4次） :

\mathbf{q}^{n+1} = \mathbf{q}^n + \frac{h}{6}(\mathbf{k}_1 + 2\mathbf{k}_2 + 2\mathbf{k}_3 + \mathbf{k}_4)

高精度だが計算コストが高い
剛体の姿勢更新には注意が必要（回転の非線形性）

8.7.4 実装上の注意点

姿勢表現の選択 :

オイラー角: 特異点（gimbal lock）がある
回転行列: 正規直交性を維持する必要
四元数: 特異点なし、ノルム制約 $|\mathbf{q}| = 1$ を維持

質量行列の逆行列計算 :

複雑構造物体への適用可能性の検証

9.1 この章の目的と位置づけ

9.1.1 なぜこの章が必要か

これまでの章では、剛体力学の基本的な概念から出発し、単一剛体の運動、複数剛体の連結系、拘束条件の扱い方、そして多体システムへの一般化された定式化まで、段階的に理論を構築してきた。しかし、理論の理解と「実際に複雑な物体の運動を記述できる」ことの間には、依然として大きなギャップが存在する。

このギャップは以下の問いに集約される：

理論の射程範囲 ：これまで学んだ枠組みは、人体のような複雑な構造にも本当に適用可能なのか？
モデル化の選択 ：実在物体を剛体系としてモデル化する際、何を保持し何を捨てるべきか？
計算実行可能性 ：理論的に正しい方程式が導出できても、それは実際に解ける形になっているか？

この章の目的は、これらの問いに答えることで、学習者が自信を持って複雑な物体の運動記述に取り組めるようにすることである。

9.1.2 検証の戦略

本章では、以下の三段階で適用可能性を検証する：

具体例による実証 ：人体歩行モデルという具体的な複雑システムを題材に、これまでの理論がどのように適用されるかを示す
モデル化の意思決定プロセスの明示 ：実在物体から数理モデルへの変換において、どのような判断が必要かを体系的に整理する
一般的な複雑系への拡張可能性の議論 ：人体以外の複雑構造（多関節ロボット、動物の運動など）への適用における共通原理と個別配慮を明確化する

9.2 人体歩行モデルの構築：理論の統合的適用

9.2.1 問題設定：人体という複雑構造をどう捉えるか

人体は、約200の骨、600以上の筋肉、無数の軟組織からなる極めて複雑なシステムである。この複雑性に対して、剛体力学のアプローチは以下の根本的な問いを投げかける：

「このような複雑な対象を、有限個の剛体の連結系として近似することは妥当なのか？」

この問いに答えるには、 モデルの目的 を明確にする必要がある。歩行運動を記述したい場合：

保持すべき特性 ：
- 重心の軌道（前進運動の本質）
- 主要関節（股関節、膝関節、足関節など）の角度変化
- 地面との接触・離脱のタイミング
- 全体の力学的エネルギーの流れ
捨象可能な特性 ：
- 個々の筋繊維の収縮ダイナミクス
- 軟組織の弾性変形
- 血流や呼吸による質量分布の微小変化
- 皮膚や筋肉表面の局所的な振動

この取捨選択の基準は、 時間スケールと空間スケール の分離にある。歩行の一周期（約1秒）において支配的な運動は、骨格系の大域的な配置変化であり、これは剛体近似で十分に捉えられる。

9.2.2 剛体分割の決定：セグメンテーション戦略

人体を剛体系としてモデル化する最初のステップは、 どこで身体を分割するか を決定することである。

分割の原理

実在の人体において、相対運動が顕著に生じる場所は関節である。したがって、基本的な分割戦略は：

「各関節で身体を分割し、関節間の骨格部分を一つの剛体とみなす」

例えば、下肢の二次元モデルでは：

剛体1（胴体） ：腰部から頭部までを一つの剛体
剛体2（大腿） ：股関節から膝関節まで
剛体3（下腿） ：膝関節から足関節まで
剛体4（足部） ：足関節から足先まで

左右を区別すれば、下肢だけで7剛体（胴体1 + 左右各3）となる。

分割の妥当性の検証

この分割が妥当であるためには：

剛性の仮定 ：各セグメント内部での相対変形が、セグメント間の相対運動に比べて十分小さいこと
- 大腿部の骨は歩行中ほとんど変形しない → 妥当
- 足部は実際には曲がるが、第一近似では剛体とみなせる
質量分布の一定性 ：各剛体内での質量分布が時間的に変化しないこと
- 筋肉の収縮による局所的な質量移動は、セグメント全体の慣性特性に与える影響は小さい
連結点の明確性 ：剛体間の連結点（関節）が明確に定義できること
- 解剖学的な関節中心は、X線や運動学的データから推定可能

9.2.3 座標系の設定

全体座標系（慣性系）

地面に固定された座標系を設定する：

原点：観察開始時の適当な地点（例：最初の着地点）
$$x$$ 軸 ：進行方向（前向きを正）
$$y$$ 軸 ：鉛直上向き
$$z$$ 軸 ：横方向（右手系を構成）

二次元近似（矢状面内運動）では、 $$z$$ 方向の運動を無視し、 $$(x, y)$$ 平面内の運動のみを考える。

各剛体の局所座標系

各剛体 $$i$$ に固定された座標系 $\{B_i\}$ を設定する。人体セグメントでは、以下の慣例が一般的：

原点：近位端の関節中心（例：大腿剛体なら股関節中心）
$$x_i$$ 軸 ：セグメントの長軸方向（近位から遠位へ）
$$y_i$$ 軸 ：横方向
$$z_i$$ 軸 ：矢状面に垂直な方向

この設定により、各剛体の慣性テンソルが局所座標系で対角化されやすくなる。

9.2.4 一般化座標の選択

自由度の特定

二次元平面内の人体モデル（下肢のみ、7剛体）の場合：

各剛体は平面内で3自由度（ $x, y, \theta$ ）
7剛体で最大 $3 \times 7 = 21$ 自由度

しかし、関節による拘束がある：

股関節（ヒンジ）：左右各1拘束 → 2拘束
膝関節（ヒンジ）：左右各1拘束 → 2拘束
足関節（ヒンジ）：左右各1拘束 → 2拘束

さらに、地面との接触による拘束：

接地脚がある場合、接地点での滑りなし条件など

一般化座標の設計

独立な一般化座標として、以下を選択する：

絶対座標方式 （胴体基準）：

\mathbf{q} = \begin{bmatrix} x_{\text{pelvis}} \\ y_{\text{pelvis}} \\ \theta_{\text{trunk}} \\ \theta_{\text{left hip}} \\ \theta_{\text{left knee}} \\ \theta_{\text{left ankle}} \\ \theta_{\text{right hip}} \\ \theta_{\text{right knee}} \\ \theta_{\text{right ankle}} \end{bmatrix}

ここで：

$(x_{\text{pelvis}}, y_{\text{pelvis}})$ ：骨盤（胴体）の位置
$\theta_{\text{trunk}}$ ：胴体の傾き角
$\theta_{\text{joint}}$ ：各関節角度（解剖学的参照位置からの偏差）

この選択により、自由度は9となる（二次元、非接地時）。

重要な概念的注意 ：

関節角度 $\theta_{\text{hip}}$ は、大腿部の絶対角度ではなく、 胴体に対する相対角度 として定義することが多い。これにより、座標間の依存関係が明示的になる：

\theta_{\text{thigh, absolute}} = \theta_{\text{trunk}} + \theta_{\text{hip}}

9.2.5 運動学的記述：位置・速度・加速度

各剛体の配置は、一般化座標 $\mathbf{q}$ から一意に決定される。第8章で学んだ一般理論を適用する。

剛体 $$i$$ の重心位置

例として、左大腿部の重心位置を求める：

\mathbf{r}_{\text{thigh}}(q) = \mathbf{r}_{\text{pelvis}} + R(\theta_{\text{trunk}} + \theta_{\text{hip}}) \mathbf{d}_{\text{thigh, local}}

ここで：

$\mathbf{r}_{\text{pelvis}} = (x_{\text{pelvis}}, y_{\text{pelvis}})^T$
$R(\theta)$ は二次元回転行列：

R(\theta) = \begin{bmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{bmatrix}

$\mathbf{d}_{\text{thigh, local}}$ は局所座標系での股関節から大腿重心への変位ベクトル

ヤコビ行列の構築

速度は：

\dot{\mathbf{r}}_{\text{thigh}} = \frac{\partial \mathbf{r}_{\text{thigh}}}{\partial \mathbf{q}} \dot{\mathbf{q}} = J_{\text{thigh}}(\mathbf{q}) \dot{\mathbf{q}}

ヤコビ行列 $J_{\text{thigh}}$ の各列は、対応する一般化座標の変化が大腿重心位置に与える影響を表す：

J_{\text{thigh}} = \begin{bmatrix} \frac{\partial \mathbf{r}_{\text{thigh}}}{\partial x_{\text{pelvis}}} & \frac{\partial \mathbf{r}_{\text{thigh}}}{\partial y_{\text{pelvis}}} & \frac{\partial \mathbf{r}_{\text{thigh}}}{\partial \theta_{\text{trunk}}} & \cdots \end{bmatrix}

具体的計算例 ：

\frac{\partial \mathbf{r}_{\text{thigh}}}{\partial x_{\text{pelvis}}} = \begin{bmatrix} 1 \\ 0 \end{bmatrix}, \quad \frac{\partial \mathbf{r}_{\text{thigh}}}{\partial \theta_{\text{hip}}} = \frac{\partial R(\theta_{\text{trunk}} + \theta_{\text{hip}})}{\partial \theta_{\text{hip}}} \mathbf{d}_{\text{thigh, local}}

回転行列の微分は：

\frac{\partial R(\theta)}{\partial \theta} = \begin{bmatrix} -\sin\theta & -\cos\theta \\ \cos\theta & -\sin\theta \end{bmatrix} = R(\theta) \begin{bmatrix} 0 & -1 \\ 1 & 0 \end{bmatrix}

9.2.6 運動エネルギーの構成

系全体の運動エネルギーは、各剛体の並進運動エネルギーと回転運動エネルギーの和：

T = \sum_{i=1}^{7} \left( \frac{1}{2} m_i \|\dot{\mathbf{r}}_i\|^2 + \frac{1}{2} I_i \dot{\theta}_i^2 \right)

ここで、 $$I_i$$ は剛体 $$i$$ の重心周りの慣性モーメント（二次元の場合はスカラー）。

質量行列形式 ：

第8章で見たように、これは一般化座標で表現すれば：

T = \frac{1}{2} \dot{\mathbf{q}}^T M(\mathbf{q}) \dot{\mathbf{q}}

質量行列 $M(\mathbf{q})$ の各要素は：

M_{jk}(\mathbf{q}) = \sum_{i=1}^{7} \left( m_i J_{i}^T J_i + I_i \frac{\partial \theta_i}{\partial q_j} \frac{\partial \theta_i}{\partial q_k} \right)

ここで $$J_i$$ は剛体 $$i$$ の重心位置に対するヤコビ行列。

重要な洞察 ：

質量行列は配置 $\mathbf{q}$ に依存する。これは、同じ一般化速度 $\dot{\mathbf{q}}$ でも、身体の姿勢によって実際の運動エネルギーが異なることを意味する。例えば、脚を伸ばした状態で回転するのと曲げた状態で回転するのでは、必要な運動エネルギーが異なる。

9.2.7 ポテンシャルエネルギー

重力場における位置エネルギー：

V = \sum_{i=1}^{7} m_i g y_i(\mathbf{q})

ここで $y_i(\mathbf{q})$ は剛体 $$i$$ の重心の鉛直座標、 $$g$$ は重力加速度。

9.2.8 拘束条件のモデル化

関節の運動学的拘束

前述の一般化座標の選択により、関節の拘束条件は既に座標定義に組み込まれている（ 自然座標 の利点）。例えば、膝関節をヒンジとして扱うことで、膝関節での横方向の相対移動や軸周り以外の回転は自動的に排除される。

地面接触拘束

歩行では、足部と地面の接触が重要な役割を果たす。接地相では：

位置拘束 （滑りなし条件）：

\phi_{\text{contact}}(\mathbf{q}) = \mathbf{r}_{\text{foot contact point}}(\mathbf{q}) - \mathbf{r}_{\text{ground}} = \mathbf{0}

これは ホロノミック拘束 であり、第7章で学んだラグランジュ未定乗数法で扱える。

片側拘束 （地面は押すが引かない）：

\lambda_{\text{normal}} \geq 0, \quad y_{\text{contact}} \geq 0, \quad \lambda_{\text{normal}} \cdot y_{\text{contact}} = 0

これは 不等式拘束 であり、より高度な扱いが必要だが、基本原理は同じである。

9.2.9 外力のモデル化

重力

既にポテンシャルとして組み込まれているが、一般化力の形式では：

Q_j^{\text{gravity}} = -\frac{\partial V}{\partial q_j} = -\sum_{i=1}^{7} m_i g \frac{\partial y_i}{\partial q_j}

地面反力

接地時、地面からの反力 $\mathbf{F}_{\text{ground}}$ が足部接地点に作用する。これを一般化力に変換：

Q_j^{\text{ground}} = \sum_{\text{contact points}} \mathbf{F}_{\text{ground}}^T \frac{\partial \mathbf{r}_{\text{contact}}}{\partial q_j}

これは第8章で学んだ 仮想仕事の原理 の直接的な応用である。

関節トルク（能動力）

人体では、筋肉が関節にトルクを生成する。関節 $$k$$ に作用するトルク $\tau_k$ は、対応する一般化力として直接現れる：

Q_k^{\text{muscle}} = \tau_k

注意：この章の目的は運動方程式の記述であり、 $\tau_k$ 自体がどう決定されるか（制御の問題）は扱わない。 $\tau_k$ は既知の入力または未知の変数として扱う。

9.2.10 運動方程式の導出

ラグランジュ方程式（拘束なしの場合）：

\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{\mathbf{q}}}\right) - \frac{\partial L}{\partial \mathbf{q}} = \mathbf{Q}

ここで $$L = T - V$$ はラグランジアン、 $\mathbf{Q}$ は非保存力による一般化力。

展開すると：

M(\mathbf{q})\ddot{\mathbf{q}} + C(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} + \mathbf{g}(\mathbf{q}) = \boldsymbol{\tau} + \mathbf{Q}^{\text{ext}}

ここで：

$M(\mathbf{q})$ ：質量行列
$C(\mathbf{q}, \dot{\mathbf{q}})$ ：コリオリ・遠心力項
$\mathbf{g}(\mathbf{q})$ ：重力項
$\boldsymbol{\tau}$ ：関節トルクベクトル
$\mathbf{Q}^{\text{ext}}$ ：地面反力などの外力

拘束がある場合 （接地相）：

M(\mathbf{q})\ddot{\mathbf{q}} + C(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} + \mathbf{g}(\mathbf{q}) = \boldsymbol{\tau} + J_{\text{contact}}^T(\mathbf{q}) \boldsymbol{\lambda}

\phi_{\text{contact}}(\mathbf{q}) = \mathbf{0}, \quad \dot{\phi}_{\text{contact}}(\mathbf{q}, \dot{\mathbf{q}}) = \mathbf{0}

ここで $\boldsymbol{\lambda}$ は拘束力（地面反力）、 $J_{\text{contact}}$ は接触点のヤコビ行列。

9.2.11 この記述の意味：何が達成されたか

上記の導出により、以下が達成された：

完全な数理モデル ：人体歩行の力学を記述する微分代数方程式系が得られた
計算可能性 ：各項（ $$M$$ , $$C$$ , $\mathbf{g}$ など）は、与えられた $\mathbf{q}$ に対して数値計算可能
解析可能性 ：この方程式系を数値積分することで、歩行運動のシミュレーションが可能
逆問題への適用 ：実測された歩行データから逆に $\boldsymbol{\tau}$ （関節トルク）を推定することも可能

9.3 モデル化における意思決定プロセス

9.3.1 抽象化レベルの選択

前節の人体モデルは、多くの意思決定の結果である。ここでは、そのプロセスを体系化する。

レベル0：完全な物理的記述

原理的には、人体を構成する全ての原子の運動を量子力学で記述することも可能である。しかし、これは：

計算不可能 ：原子数 $\sim 10^{28}$
情報過剰 ：歩行の理解に原子レベルの情報は不要
時間スケール不整合 ：原子振動 $\sim 10^{-15}$ 秒、歩行 $\sim 1$ 秒

レベル1：連続体力学

人体を弾性体・流体として扱う。筋肉の変形、血流などを記述可能。しかし：

依然として複雑 ：偏微分方程式系、境界条件の複雑さ
主要な運動様式との不一致 ：歩行の主要な運動は骨格の大域的配置変化

レベル2：剛体系近似

骨格を剛体とし、関節で連結する。筋肉は関節トルク源として抽象化。

計算可能 ：常微分方程式系、有限自由度
本質的情報の保存 ：重心軌道、関節角度など主要な運動学的量を再現
時間スケールの整合 ：骨格運動の時間スケールと一致

レベル3：さらなる簡約化

特定の問いに応じて、さらに簡略化することも可能：

倒立振子モデル ：胴体と脚全体を一つの剛体棒とみなす（歩行の基本的な重心ダイナミクスの理解）
質点モデル ：重心の運動のみを追跡（エネルギーコスト推定）

9.3.2 適切な抽象化レベルの選択基準

どのレベルを選ぶべきかは、 問いの性質 に依存する：

問い	推奨レベル	理由
歩行時の関節トルクは？	レベル2（剛体系）	関節ごとの情報が必要
歩行のエネルギーコストは？	レベル2または3	重心の運動で大部分が決まる
筋肉の損傷リスクは？	レベル1（連続体）	筋肉内の応力分布が必要
歩行の安定性原理は？	レベル3（簡約）	本質的メカニズムの理解優先

一般原則 ：

可能な限り単純なモデルを用いよ。ただし、答えるべき問いに必要な情報が失われない範囲で。（"Everything should be made as simple as possible, but not simpler." — Einstein）

9.3.3 パラメータ同定の問題

剛体系モデルが構築されても、数値的に解くには各パラメータの値が必要：

各セグメントの質量 $$m_i$$
慣性モーメント $$I_i$$
セグメント長、重心位置
関節の運動学的パラメータ

人体測定学的データの利用

これらのパラメータは、以下の方法で推定される：

直接測定 ：身長、体重、セグメント長などを測定
回帰式の利用 ：身長・体重から各セグメント質量を推定する経験式
画像データ ：CTスキャンやMRIから質量分布を推定
逆同定 ：実測された運動データから運動方程式を用いて逆算

Winter回帰式：実用的パラメータ推定手法

バイオメカニクス分野で広く採用されている実用的手法として、 Winter回帰式 がある。1970年代から1980年代にかけてDavid A. Winterらが開発したこの手法は、死体解剖データや放射線画像データを基に構築された経験式であり、非侵襲的に得られる基本的な人体計測値（身長、体重、セグメント長）のみから各セグメントの慣性特性を推定できる[^1]。

Winter回帰式の一般的な形式：

セグメント質量 ：全体重に対する比率として表現

m_i = \alpha_i \cdot m_{\text{total}}

ここで $\alpha_i$ は各セグメント固有の比率係数（例：大腿部は約0.1、下腿部は約0.0465）

重心位置 ：セグメント長に対する比率として表現

l_{\text{COM}, i} = \beta_i \cdot L_i

ここで $\beta_i$ は近位端からの距離比率、 $$L_i$$ はセグメント長

慣性モーメント ：質量とセグメント長の積に対する回転半径の二乗として表現

I_i = m_i \cdot (k_i \cdot L_i)^2

ここで $$k_i$$ は回転半径比

実務上の利点 ：

非侵襲性 ：高価な画像診断装置や複雑な測定機器が不要
実用的精度 ：歩行解析における逆動力学計算に十分な精度を提供
標準化 ：多くの研究で採用されており、結果の比較可能性が高い

ただし、Winter回帰式は主に白人成人男性のデータから構築されているため、異なる人種、年齢層、体型に適用する際には、モデル選択の影響を考慮する必要がある[^2]。個人ベースの人体計測データを用いた予測回帰モデルの研究では、集団平均データに依存する従来手法から、個人の特性をより反映した予測モデルへの移行が検討されている。

パラメータの不確実性と妥当性評価

重要な認識 ：

パラメータの不確実性は避けられない。したがって、モデル予測にも不確実性が伴う。モデルの妥当性は、「完全に正確な予測ができるか」ではなく、「問いに答えるのに十分な精度があるか」で判断される。

妥当性検証の実践的アプローチ ：

感度解析 ：パラメータの変動がモデル出力に与える影響を定量化
実測データとの比較 ：モーションキャプチャシステムで得られた運動学的データとの照合
力学的整合性チェック ：計算された地面反力が物理的に妥当な範囲にあるかの確認

9.4 計算実行可能性の検証

9.4.1 方程式の数値的性質

導出された運動方程式：

結論

教科書設計の到達点

本カリキュラムは、剛体力学の初心者が人体歩行のような複雑な多体システムの運動方程式を記述できるようになるまでの、体系的な学習パスを提供する。9つのステップを通じて、以下の能力獲得を実現する設計となっている。

達成される主要能力

複雑構造の数理モデル化能力 : 人体のような関節で連結された多剛体系を、座標系・拘束条件・運動方程式の組み合わせで完全に記述できる
物理的直観と数学的厳密性の統合 : 「剛体とは何か」「力のモーメントとは何か」といった基礎概念を正確に理解した上で、Lagrange方程式やNewton-Euler方程式を適用できる
問題分解と段階的構築 : 複雑なシステムを単純な要素に分解し、拘束条件で再構成するアプローチを習得できる

設計の特徴と妥当性

本教科書設計は「概念の深さ・正確性最優先」の方針に基づき、以下の特徴を持つ:

厳密性と理解可能性の両立 : 剛体の定義から始まり、各概念が「なぜ必要か」「何を可能にするか」を明確にしながら進行する
段階的複雑性の管理 : 単一剛体→回転導入→多体系→拘束条件という順序で、各段階の問題意識を明示しながら理論を積み上げる
制御理論の明確な除外 : 筋肉制御などには触れず、運動方程式の記述と解析に焦点を限定している

実現可能性の評価

強み :

人体歩行モデル（股関節・膝関節・足首を持つ3剛体系）を具体例として、実際に運動方程式を導出する最終ステップが含まれている
各ステップが前段階の知識のみを前提とし、論理的飛躍がない
MathJax数式と懇切丁寧な文章説明の組み合わせにより、初心者が独学可能

課題 :

概念説明の厳密性を優先した結果、全体のボリュームが増大する傾向がある
ステップ4（回転運動）からステップ6（多体相互作用）への移行で、学習負荷が急増する可能性がある
実際の人体のような高自由度システムでは、計算の複雑さが理解の障壁となりうる

教育的価値の総合評価

本カリキュラムは、以下の点で初心者向け剛体力学教科書として高い教育的価値を持つ:

長期的応用可能性 : 表面的な公式暗記ではなく、概念の本質理解を重視することで、様々な複雑システムへの応用が可能になる
誤解の事前予防 : 「剛体は完全に硬い物体」といった初学者の典型的誤解を、定義段階で正確に修正する
問題意識の明確化 : 各理論が「何を解決するために導入されるか」を常に提示し、学習動機を維持する

ただし、この設計は「完全性と正確性を重視する学習者」に最適化されており、「最短経路で実用的計算手法を習得したい」学習者には過剰に詳細な可能性がある。

推奨事項

カリキュラム実装に関する推奨

1. 段階的ボリューム調整の実施

推奨アクション :

ステップ4（回転運動）を2つのサブステップに分割: (4a) 平面内回転の導入、(4b) 3次元回転と慣性テンソル
ステップ7（拘束条件）で、「ホロノミック拘束」と「非ホロノミック拘束」を明確に区別し、初学者には前者に焦点を当てる
各ステップ末尾に「習得確認のための演習問題」を3-5問配置し、理解度を自己評価できるようにする

理由 : 概念の深い説明を維持しながら、学習負荷の急増を防ぎ、段階性を改善する。

2. 具体例と抽象理論のバランス強化

推奨アクション :

各ステップで「身近な具体例」→「数学的定式化」→「具体例への適用」のサイクルを徹底する
人体歩行以外にも、「振り子」「転がる円盤」「体操選手の宙返り」など、複数の応用例を各ステップで提示する
ステップ9（複雑構造物体への適用）で、人体歩行だけでなく「ロボットアーム」「多関節構造物」などの例も追加する

理由 : 厳密な説明を維持しつつ、抽象概念を複数の文脈で理解することで、応用範囲を拡大し学習動機を高める。

3. 計算複雑性への対処戦略

推奨アクション :

ステップ8（多体システムへの一般化）で、「手計算による理解」と「数値計算ツールの活用」を明確に区別する
3剛体以上のシステムについては、「運動方程式の導出方法」に焦点を当て、「実際の数値解法」は付録または発展的内容とする
Python/MATLABなどのコード例を付録として提供し、理論理解後の実装をサポートする（ただし本文では扱わない）

理由 : 複雑なシステムで計算が学習の障壁にならないよう、理解すべき本質（方程式の導出）と実用的手段（数値計算）を分離する。

教科書全体の構成に関する推奨

4. モジュール式構成の採用

推奨アクション :

「コア教科書」(ステップ1-6) と「応用編」(ステップ7-9) に分冊化を検討する
各ステップを独立した章として、必要に応じて順序を調整可能にする（例: 回転運動を後回しにする学習パスも許容）
「最短学習パス」（重要概念のみ）と「完全学習パス」（全詳細説明）を明示し、学習者が選択できるようにする

理由 : ボリューム増大のリスクを管理しつつ、異なる学習目標を持つ読者に対応する。

5. 概念マップと学習進捗の可視化

推奨アクション :

教科書冒頭に「概念依存関係図」を配置し、各ステップがどの概念に依存し、何を可能にするかを視覚化する
各章の冒頭に「このステップの位置付け」セクションを設け、全体像との関係を常に明示する
「人体歩行モデルの完成」までの進捗を示すマイルストーンを各ステップで提示する

理由 : 詳細な説明による「森を見失う」リスクを軽減し、学習者が常に全体像を意識できるようにする。

概念説明の質に関する推奨

6. 鍵概念の「3層説明法」の導入

推奨アクション :
各重要概念（剛体、慣性モーメント、一般化座標など）について:

第1層 : 日常的な言葉での直観的説明（2-3文）
第2層 : 数学的定義と厳密な説明（本文）
第3層 : よくある誤解とその訂正（コラム形式）

理由 : 初心者向けの理解可能性と数学的厳密性を、階層的説明により両立させる。

7. 「なぜこの理論が必要か」の明確化

推奨アクション :

各ステップの冒頭に「問題意識」セクションを設け、前段階の限界や新たに記述したい現象を具体例で示す
Lagrange形式主義など高度な理論の導入時には、「Newton方程式との比較表」を作成し、利点を明示する
「この理論を使わないと何が困るか」を常に示す

理由 : 理論の必要性を理解することで、学習動機が維持され、抽象的な数学も意味を持つ。

リスク管理に関する推奨

8. 冗長性の定期的監査

推奨アクション :

各章完成後に「説明の最小化チェック」を実施: 各段落が「新しい情報」を提供しているか確認
「懇切丁寧な説明」と「冗長な繰り返し」を区別する基準を設定（例: 異なる文脈での説明は許容、同じ内容の繰り返しは削除）
ベータリーダー（初学者）からのフィードバックを収集し、「理解に必要だった説明」と「不要だった説明」を識別

理由 : 厳密性追求によるボリューム増大を管理し、効率性を確保する。

9. 段階的リリースとフィードバック収集

推奨アクション :

ステップ1-3（基礎部分）を先行公開し、初学者からの理解度データを収集
「つまずきポイント」を特定し、後続ステップの説明方法に反映
最終ステップ（人体歩行モデル）を実際に学習者が実装できるか、ケーススタディで検証

理由 : 理論的完全性と実際の学習効果のギャップを早期に発見し、修正する。

最終推奨: 「剛体力学の階段」としての価値最大化

本教科書が目指すべき最終的な姿は、単なる知識の羅列ではなく、 初学者が確実に登れる階段 である。そのために:

各ステップの「習得証明」を明示 : 「このステップを終えたら、この問題が解ける」という具体的なベンチマークを設定
失敗からの学習経路を組み込む : 「よくある間違い」とその修正方法を積極的に提示
最終到達点の明確化 : 「人体歩行の運動方程式が書ける」という目標を、数式で具体的に示し、学習のゴールを可視化

これらの推奨事項を実装することで、「概念の深さ・正確性最優先」という方針を維持しながら、実用的で学習可能な教科書が完成すると期待される。