技術的概要：多接触歩行における振り子構造の創発

問題定義

線形逆振り子モデル（LIPM）は、脚型歩行制御において広く用いられる仮定であり、しばしば制御器の目的から自然に創発するのではなく、モデリング制約として課されている。高さの変動や角運動量への拡張は存在するが、振り子形式は通常「手動で」付加される。本論文はより狭い問いを検証する：どのような特定の目的関数が、振り子力パターンを最適解として創発させ、この図式がどこで破綻するか？

著者らは、接触力および法線加速度に関する標準的な項に加えて、角運動量の変化率（$\|\dot{H}_G\|$）を罰する小型の重心最適制御問題（OCP）を分析する。目的は、角運動量変化率の最小化が、系を自然に振り子多様体（$\dot{H}_G = 0$）へと導くかどうかを判断し、異なる立ち姿勢構成（フルランク対 $N=2$）およびタスク要件の下でこの挙動の限界を特徴づけることである。

手法

本論文は、以下のコスト関数から構成されるネット・レンチ OCP（定義 1）を定式化し、これを最小化する：

1. 角運動量ペナルティ: $\alpha \|\dot{H}_G\|^2$
2. タスク追跡: $\lambda \|\dot{H}_G - \dot{H}_{task}\|^2$（バランス以外のタスクの場合）
3. 法線加速度ペナルティ: $\beta (n^\top \ddot{c})^2$
4. 力正則化: $\gamma \sum \|f_i\|^2$

分析は、点質量四足歩行ロボットおよび Unitree Go1 ロボットを用いた MuJoCo 上での数値実験（オープンループ QP とクローズドループトルクレベル MPC の両方）によって検証され、3 つの理論的段階で進行する：

1. フルランク姿勢分析: モメントヤコビアンの特異値分解（SVD）を用いて、角運動量重み $\alpha \to \infty$ における漸近挙動を導出する。
2. $N=2$ 姿勢分析: モメント演算子がランク 1 となるランク不足の場合（例：対角トロット）を調査し、避けられない残差と摩擦円錐制約を明らかにする。
3. タスク駆動分析: タスク追跡項（$\lambda > 0$）を追加することが、最適解を振り子多様体からどのようにずらすかを検討する。

主要な貢献と理論的結果

1. フルランク姿勢における振り子構造の創発

• 定理 1: フルランク姿勢において、最適解は角運動量重み $\alpha$ に対して $O(1/\alpha)$ の速度で振り子力パターン（$\dot{H}_G \to 0$）へと漸近的に収束する。
• 幾何定数: 比例定数はモメントヤコビアンの SVD によって決定され、これは直接足幅の幾何学に対応する。
• 古典的モデルとの関連: この振り子多様体上では、ゼロモーメント点（ZMP）は自然に仮想ピボット点と等しくなり、発散運動成分（DCM）の軌道方程式は仮定ではなく帰結として現れる。

2. $N=2$ 姿勢（対角トロット）における根本的限界

• 定理 5（幾何学的フロア）: $N=2$ 姿勢では、モメントヤコビアンはランク 1 である。2 つの足を結ぶ線方向（$\hat{D}$）には、力の再配分によっても残る角運動量変化率の残差が存在する。これにより、$\|\dot{H}_G\|$ に対して非ゼロの下限（フロア）が生じる。
• 定理 6（摩擦円錐のキンク）: 2 番目の障壁は摩擦円錐から生じる。垂直モーメント成分を打ち消すために必要な力が摩擦円錐の外部にある場合、臨界水平加速度（$a_x^*$）において「非滑らかなキンク」が発生する。この点を越えると、達成可能な $\|\dot{H}_G\|$ は幾何学的フロアを厳密に超える。
• 含意: 対角トロットでは、振り子方向がしばしば摩擦円錐の外側に傾く（例：$\mu=0.6$ の場合、$40.6^\circ$ 対 $31^\circ$）。これにより振り子解は実行不可能となり、$\alpha$ に依存しないフロアが生じる。

3. タスク駆動による逸脱

• 定理 7: タスク項（$\lambda > 0$）が追加されると、最適角運動量変化率は滑らかな補間となる：$\dot{H}_G^* = \frac{\lambda}{\alpha + \lambda} \dot{H}_{task}^*$。
• 相図: 系は「バランス優位」領域（振り子的、$\lambda/\alpha \ll 1$）から「タスク優位」領域（非振り子的、$\lambda/\alpha \gg 1$）へと滑らかに遷移する。これにより、同一ロボットが静的バランス時には LIPM 的に見えるが、機敏な機動時には非振り子的になる理由が、制御器の切り替えなしに説明される。

実験結果

著者らは Unitree Go1 および点質量モデルを用いてこれらの主張を検証した：

• 速度検証（点質量）: $\alpha$ が増加するにつれて、振り子力パターンからの逸脱は 227 倍減少し、対数 - 対数グラフの傾きは -1 となった。これは $O(1/\alpha)$ の速度（定理 1）を確認するものである。
• 幾何定数（Go1 N=4）: 4 足姿勢において、実測のスケーリング定数（$K_e = 8.4$）は、足幅から導出された解析的予測（$K_a = 9.8$）と約 16% の範囲で一致した。モメントヤコビアンの特異値は、物理的な足幅と小数点以下 4 桁まで一致した。
• フロア検証（Go1 N=2）: 対角トロットにおいて、$\|\dot{H}_G\|$ は予測された幾何学的フロア $0.2835 \, \text{m}^2/\text{s}^2$ で飽和し、解析的限界と 4 桁まで一致した。
• クローズドループの現実: クローズドループトルクレベル MPC において、フロアは持続したが、静的予測（$\approx 1.71 \, \text{m}^2/\text{s}^2$）よりも高かった。著者らはこの乖離を、準静的分析で捕捉されていない動的効果（スイング時の衝撃、接触遷移、サーボ誤差）に起因すると帰着させた。
• ZMP 同一性: 数値実験は、振り子多様体上において、力ベースの ZMP が仮想ピボットに収束することを確認し、LIPM と OCP 最適解の間の理論的リンクを検証した。

意義と主張

本論文は、LIPM と最適制御の間のリンクを明示的かつ検証可能なものにすると主張する。以下の点を論じている：

1. LIPM は単なる選択ではなく帰結である: 十分な大きさの重み $\alpha$ が与えられれば、フルランク姿勢において角運動量変化率を最小化することで、振り子パターンは自然に創発する。
2. 限界は予測可能である: $N=2$ 姿勢における振り子仮定の破綻は、制御器の失敗ではなく、根本的な幾何学的および摩擦的限界である。本論文は、下限および摩擦制約が効く臨界加速度に対する閉形式の式を提供する。
3. バランスとタスクの統合的視点: 比 $\lambda/\alpha$ は、バランス中心（LIPM 的）からタスク中心（機敏な）挙動への遷移を説明する連続パラメータを提供する。

著者らはクローズドループの結果については謙虚であり、動的シミュレーションにおいて定性的な挙動（フロアの存在）は成立するものの、静的理論とクローズドループの現実との間の定量的な乖離は有意であり、おそらくモデル化されていない動的要素（コンプライアンス、スリップ、有限の制御レート）に起因すると指摘している。彼らはこれらの動的乖離を解決するものではなく、むしろそれらの動的要素が作用する静的境界を定義するものであると主張している。