SoFFT: Spatial Fourier Transform for Modeling Continuum Soft Robots

技術的概要：SoFFT：連続体ソフトロボットのモデリングのための空間フーリエ変換

問題定義
柔軟な材料から構成される連続体ソフトロボット（CSRs）は、理論的に無限の自由度（DoFs）を有し、非構造化環境において高い適応性を示す。しかし、この無限次元性は制御器の設計とシミュレーションを複雑化する。コッセラットロッド理論（CRT）は、偏微分方程式（PDEs）を介して CSRs を時間変化する曲線（背骨）としてモデル化する効果的な枠組みを提供するが、実用的な実装には自由度を削減するための空間離散化が必要である。既存の手法、例えば区間定数曲率（PCC）、区間定数ひずみ（PCS）、幾何学的可変ひずみ（GVS）は、セグメント数や関数基底の選択（多項式、三角関数など）に関するヒューリスティックな選択に依存している。これらのアプローチは、精度と計算効率のバランスを取るために最適な基底を選択し、必要な最小セグメント数を決定するための統一的な理論的根拠を欠くことが多い。さらに、固有直交分解（POD）のようなデータ駆動型の削減手法は、すべてのひずみモードを結合させる傾向があり、モード固有のダイナミクスやスケーリングの問題を曖昧にする可能性がある。

手法
著者らは、CSR の背骨を空間信号として扱い、その変形を解析するためにフーリエ変換（FT）を適用する「SoFFT」という枠組みを提案する。核心的な手法は以下の通りである：

1. 空間信号定式化: ひずみ場 $\xi(s, t)$ を空間 $s \in [0, L]$ 上で定義された信号として扱う。連続空間フーリエ変換（SFT）を定義し、ひずみ場を空間周波数領域（波数 $k$）へマッピングする。
2. 周波数領域解析: 本論文は、空間周波数領域内で CRT の静力学および力学方程式を再解釈する。これにより、剛性行列がひずみ場に対する畳み込みフィルタとして作用することが明らかになる。著者らは、既存の離散化手法（例：PCS、区間アフィンひずみ）が、周波数領域において特定の信号復元器（例：ゼロ次ホールド、一次ホールド）として数学的に解釈できることを示す。
3. データ駆動型スペクトル抽出: 物理パラメータの事前知識なしに実験データから最適なモデリング基底を直接抽出する新たな手法が提案される。これには以下が含まれる：
   • センサ化: ロボットに沿ってひずみサンプルを測定する。
   • モーター・バブリング: アクチュエータに標準的な励起信号（ステップ、チャープ、白色雑音）を印加して作業空間を探索する。
   • スペクトル解析: 支配的な空間高調波（$k_{max}$）と動的特性（共振/反共振）を特定するために、空間 - 時間フーリエ変換（STFT）を計算する。
   • 基底選択: 定義された辞書（多項式、三角関数）から、ひずみ場を最良に再構成しつつ自由度を最小化するスパースな最適基底関数の集合を選択するために、**基底追跡ノイズ除去（BPD）**アルゴリズムを使用する。
4. 切り捨て基準: 保持されるエネルギーが特定の精度閾値を満たすことを保証する最小のセグメント数または基底関数を決定するための、エネルギーに基づく切り捨て指標を定義する。

主要な貢献

• 統一的枠組み: 本論文は、既存の CSR モデリング戦略（PCC、PCS、GVS）を空間信号復元問題として解釈することで、単一の CRT 枠組みの下で統合する。
• 理論的根拠: 空間信号に適用されたナイキスト - シェノンの定理に基づき、特定のロボット長とひずみプロファイルに対して必要な最小セグメント数を決定するなど、ヒューリスティックな手法に対する理論的基盤を提供する。
• モード固有の解析: 全てのモードに対して結合された基底を導出する POD と異なり、提案された FT アプローチは、ひずみモード（曲げ、捩れ、せん断、伸縮）を独立して扱い、スケーリングの問題を回避し、モード固有の切り捨てを可能にする。
• データ駆動型モデリング: 実験データに基づいて最適な関数基底を特定し、ロボットの物理パラメータに依存せずにモデル次数を削減する、STFT と BPD を使用するデータ駆動型ワークフローを導入する。
• 検証: このアプローチは、数値シミュレーション（SoRoSim ツールボックスを使用）およびモーションキャプチャシステムを搭載した H サポートプロトタイプを用いた実世界の実験を通じて検証された。

結果

• 数値的検証: 円筒形および円錐形の H サポートロボットを用いたシミュレーションにより、STFT が共振および反共振ピークを効果的に捉えることが示された。解析により、関数基底の選択が空間高調波全体にわたる変形エネルギーの分布に大きく影響することが明らかになった。例えば、多項式基底はエネルギーを高調波次数に反比例して分配するのに対し、フーリエ基底は特定の高調波にエネルギーを集中させる。
• 実験的検証: 物理プロトタイプを用いた実験により、この手法が高精度で背骨形状を再構成する能力が確認された。多項式と三角関数の混合辞書を用いた BPD を使用した結果、モデルは完全次数モデルと同等の平均二乗誤差（MSE）を達成しつつ、自由度を削減した（例、1% の切り捨て閾値で 30 から 24 の自由度へ削減）。
• 摂動に対する頑健性: この手法は、外部の衝撃力によって誘発された新しい空間高調波（6–8 $m^{-1}$ の範囲）を正常に同定し、接触相互作用や不均一な変形をモデル化する能力を実証した。
• 誤差解析: 位置と姿勢の誤差は、前方運動学における誤差伝播によりロボットの長さ方向に増加することが判明したが、許容範囲内（例、5% の切り捨てによる位置誤差はロボット長さの 5% 未満）に留まった。

意義と主張
本論文は、SoFFT 枠組みがひずみ場のスペクトル特性を活用することで、CSR の変形をコンパクトかつ効率的に表現することを主張する。このアプローチは以下の点を主張する：

• 多様なモデリング戦略を共通の数学的言語の下で統合する。
• 精度と計算効率のバランスを取りながら、ヒューリスティックなモデル選択に対する厳密でデータ駆動型の代替手段を提供する。
• アクチュエータのダイナミクスや外力の効果を含む、ロボットの特定の動的および静的挙動を捉える最適な基底関数の同定を容易にする。
• CRT の PDEs を置き換えるものではなく、GVS 枠組み内で基底を最適化するためにスペクトル情報を活用する。

著者らは、現在の実装は主にオフラインであるが、高速フーリエ変換（FFT）の計算効率から、モデルベースまたは学習ベースの制御、および変化する環境条件へのリアルタイム適応などの将来のリアルタイム応用の可能性が示唆されると結論づけている。