An omnidirectional multirotor has the maneuverability of decoupled translational and rotational motions, superseding the traditional multirotors' motion capability. Such maneuverability is achieved due to the ability of the omnidirectional multirotor to frequently alter the thrust amplitude and direction. In doing so, the rotors' settling time, which is induced by inherent rotor dynamics, significantly affects the omnidirectional multirotor's tracking performance, especially in aggressive flights. To resolve this issue, we propose a novel tracking controller that takes the rotor dynamics into account and does not require additional rotor state measurement. This is achieved by integrating a linear rotor dynamics model into the vehicle's equations of motion and designing a PD controller to compensate for the effects introduced by rotor dynamics. We prove that the proposed controller yields almost global exponential stability. The proposed controller is validated in experiments, where we demonstrate significantly improved tracking performance in multiple aggressive maneuvers compared with a baseline geometric PD controller.
論文ID : 2209.10024タイトル : Geometric Tracking Control of Omnidirectional Multirotors for Aggressive Maneuvers著者 : Hyungyu Lee, Sheng Cheng, Zhuohuan Wu, Jaeyoung Lim, Roland Siegwart, Naira Hovakimyan分類 : cs.RO(ロボティクス)、cs.SY(システムと制御)、eess.SY(電気工学システムと制御)発表時期 : 2022年9月初版投稿、2025年1月更新論文リンク : https://arxiv.org/abs/2209.10024 全方向マルチロータは、並進運動と回転運動を分離した機動能力を有し、従来のマルチロータの運動能力を超えている。この機動性は、推力の大きさと方向を頻繁に変更することで実現される。このプロセスにおいて、固有のプロペラダイナミクスに起因するプロペラ安定化時間は、特に積極的な飛行において全方向マルチロータの追従性能に大きな影響を与える。この問題に対処するため、本論文はプロペラダイナミクスを考慮しながら追加のプロペラ状態測定を必要としない新しい追従制御器を提案する。線形プロペラダイナミクスモデルを航空機の運動方程式に統合し、PD制御器を設計してプロペラダイナミクスがもたらす影響を補償する。著者らは、提案された制御器がほぼ全域指数安定性を有することを証明した。実験検証により、基準となる幾何学的PD制御器と比較して、本制御器は様々な積極的機動において追従性能を大幅に改善することが示された。
中核的問題 : 全方向マルチロータが積極的機動を実行する際、プロペラダイナミクスの安定化時間が追従性能に深刻な影響を与える技術的課題 : 従来の制御器はプロペラが理想的で高速応答を示すと仮定し、プロペラダイナミクスがシステム性能に与える影響を無視している実際の要求 : 動的環境では、全方向マルチロータは速度と方向を頻繁かつ正確に変更する必要があり、制御精度に対する要求は極めて高い全方向マルチロータは任意の姿勢で重力を相殺する推力を生成でき、応用範囲を大幅に拡張する 積極的機動能力は動的環境での操作に不可欠である 既存の方法ではプロペラダイナミクスの体系的考慮が不足しており、実際の応用性能を制限している プロペラダイナミクスの無視 : ほとんどの研究はプロペラが高速(したがって無視可能)なダイナミクスを有すると仮定している追加センサの必要性 : プロペラダイナミクスを考慮する一部の方法は、専用センサまたは特殊なESCを必要としてプロペラ状態を測定する安定性保証の欠如 : 既存の方法は厳密な安定性分析を提供していない革新的な制御アーキテクチャ : プロペラダイナミクスを考慮しながら追加のプロペラ状態測定を必要としない全方向マルチロータ制御器設計の初例を提案理論的保証 : プロペラダイナミクスを含む完全なシステムにおいて、提案された制御器がほぼ全域指数安定性を有することを証明実験検証 : 8ローター全方向マルチロータプラットフォームを使用して、積極的機動における制御器の優れた性能を検証汎用的設計 : 制御器は単方向および双方向プロペラシステムに適用可能入力 : 期望位置軌跡 p d ( t ) p_d(t) p d ( t ) 、期望姿勢軌跡 R d ( t ) R_d(t) R d ( t ) およびそれらの導関数
出力 : プロペラ推力指令 f c m d f_{cmd} f c m d 制約 : プロペラダイナミクスの影響を考慮した軌跡追従制御
推力ダイナミクス(TD)モデルを採用し、推力を1次システムとして簡略化:
f i ˙ = 1 α f , i ( f c m d , i − f i ) \dot{f_i} = \frac{1}{\alpha_{f,i}}(f_{cmd,i} - f_i) f i ˙ = α f , i 1 ( f c m d , i − f i )
ここで α f , i \alpha_{f,i} α f , i は第i番目のプロペラの推力時定数である。
プロペラダイナミクスを考慮したニュートン・オイラー方程式:
m v ˙ = − α m v ¨ − α R ( F × ω ) − m g z ^ I + R F c m d m\dot{v} = -\alpha m\ddot{v} - \alpha R(F \times \omega) - mg\hat{z}_I + RF_{cmd} m v ˙ = − α m v ¨ − α R ( F × ω ) − m g z ^ I + R F c m d J ω ˙ = − α ( ω ˙ × J ω + ω × J ω ˙ + J ω ¨ ) − ω × J ω + M c m d J\dot{\omega} = -\alpha(\dot{\omega} \times J\omega + \omega \times J\dot{\omega} + J\ddot{\omega}) - \omega \times J\omega + M_{cmd} J ω ˙ = − α ( ω ˙ × J ω + ω × J ω ˙ + J ω ¨ ) − ω × J ω + M c m d
誤差を定義:
位置誤差: e p = p − p d e_p = p - p_d e p = p − p d 速度誤差: e v = v − v d e_v = v - v_d e v = v − v d 力誤差: e F = F − F d e_F = F - F_d e F = F − F d 制御則の設計:
F c m d = F d + α F d ˙ F_{cmd} = F_d + \alpha \dot{F_d} F c m d = F d + α F d ˙
ここで期望力は:
F d = R T ( − k p e p − k v e v + m g z ^ I + m v d ˙ ) F_d = R^T(-k_p e_p - k_v e_v + mg\hat{z}_I + m\dot{v_d}) F d = R T ( − k p e p − k v e v + m g z ^ I + m v d ˙ )
誤差を定義:
姿勢誤差: e R = 1 2 [ R d T R − R T R d ] ∨ / 1 + t r ( R d T R ) e_R = \frac{1}{2}[R_d^T R - R^T R_d]^{\vee}/\sqrt{1 + tr(R_d^T R)} e R = 2 1 [ R d T R − R T R d ] ∨ / 1 + t r ( R d T R ) 角速度誤差: e ω = ω − R T R d ω d e_{\omega} = \omega - R^T R_d \omega_d e ω = ω − R T R d ω d トルク誤差: e M = M − M d e_M = M - M_d e M = M − M d 制御則の設計:
M c m d = M d + α M d ˙ M_{cmd} = M_d + \alpha \dot{M_d} M c m d = M d + α M d ˙
補償項の設計 : α F d ˙ \alpha\dot{F_d} α F d ˙ および α M d ˙ \alpha\dot{M_d} α M d ˙ 項を通じてプロペラダイナミクスに起因する遅延を能動的に補償追加センサ不要 : IMUを利用して力とトルク誤差を推定し、専用プロペラ状態センサを回避幾何学的制御フレームワーク : SO(3)多様体上で制御器を設計し、特異点問題を回避分離設計 : 並進制御器と回転制御器を独立に設計し、制御構造を簡略化構成 : 8ローター固定傾斜双方向プロペラ全方向マルチロータハードウェア : BrotherHobby LPD 2306.5 2000KVモーター、Gemfan 513D 3葉3Dプロペラパラメータ : 質量1.481kg、慣性モーメント J = diag(0.020, 0.021, 0.020) kg·m²プロペラ時定数 : αf = 0.07s3つの主要なテストシナリオ:
純並進軌跡 : 並進制御システムの検証単軸回転軌跡 : 回転制御システムの検証多軸回転軌跡 : 複合運動制御の検証位置二乗平均平方根誤差(RMSE) 姿勢二乗平均平方根誤差(RMSE) 追従精度の比較 軌跡 : x d ( t ) = − 0.4 cos ( 4 π 3 t ) x_d(t) = -0.4\cos(\frac{4\pi}{3}t) x d ( t ) = − 0.4 cos ( 3 4 π t ) 、y d ( t ) = 0.4 sin ( 4 π 3 t ) y_d(t) = 0.4\sin(\frac{4\pi}{3}t) y d ( t ) = 0.4 sin ( 3 4 π t ) 、z d ( t ) = 0.6 z_d(t) = 0.6 z d ( t ) = 0.6 m結果 : 提案手法の位置RMSEは基準手法比で31%低減分析 : x軸およびy軸方向で追従性能が大幅に改善され、z軸性能は同等軌跡 : ω z B = π 2 sin ( 2 π t ) \omega_z^B = \frac{\pi}{2}\sin(2\pi t) ω z B = 2 π sin ( 2 π t ) rad/s結果 : 姿勢RMSEは基準手法比で39%低減分析 : 高角加速度条件下で姿勢追従性能が大幅に改善軌跡 : ω z B = 2 π \omega_z^B = 2\pi ω z B = 2 π rad/s、ω y I ( t ) = π 2 sin ( π t ) \omega_y^I(t) = \frac{\pi}{2}\sin(\pi t) ω y I ( t ) = 2 π sin ( π t ) rad/s結果 : 姿勢RMSE改善11%、位置追従が大幅に改善分析 : 複雑な機動における総合的性能向上を検証プロペラダイナミクスの影響 : 積極的機動では、プロペラ安定化時間が性能に影響する主要因である補償効果 : 提案手法はプロペラダイナミクスに起因する位相遅延を効果的に補償汎用性 : 異なるタイプの積極的機動すべてで改善効果を示す並進システム :
V 1 = 1 2 k p ∥ e p ∥ 2 + 1 2 m ∥ e v ∥ 2 + 1 2 α ∥ e F ∥ 2 + c 1 e p ⋅ e v V_1 = \frac{1}{2}k_p\|e_p\|^2 + \frac{1}{2}m\|e_v\|^2 + \frac{1}{2}\alpha\|e_F\|^2 + c_1 e_p \cdot e_v V 1 = 2 1 k p ∥ e p ∥ 2 + 2 1 m ∥ e v ∥ 2 + 2 1 α ∥ e F ∥ 2 + c 1 e p ⋅ e v
回転システム :V 2 = 1 2 e ω ⋅ J e ω + k R Ψ ( R , R d ) + 1 2 α ∥ e M ∥ 2 + c 2 e R ⋅ e ω V_2 = \frac{1}{2}e_{\omega} \cdot Je_{\omega} + k_R\Psi(R,R_d) + \frac{1}{2}\alpha\|e_M\|^2 + c_2 e_R \cdot e_{\omega} V 2 = 2 1 e ω ⋅ J e ω + k R Ψ ( R , R d ) + 2 1 α ∥ e M ∥ 2 + c 2 e R ⋅ e ω
制御ゲインは以下を満たす必要がある:
k p > c 1 k v 2 + 2 c 1 k v − c 1 2 m ( 4 ( k v − c 1 ) − 1 ) 、 k v > c 1 + 1 4 k_p > \frac{c_1k_v^2 + 2c_1k_v - c_1^2}{m(4(k_v - c_1) - 1)}、\quad k_v > c_1 + \frac{1}{4} k p > m ( 4 ( k v − c 1 ) − 1 ) c 1 k v 2 + 2 c 1 k v − c 1 2 、 k v > c 1 + 4 1
k R > c 2 k ω 2 λ m ( 4 ( k ω − 1 2 c 2 ) − 1 ) 、 k ω > 1 2 c 2 + 1 4 k_R > \frac{c_2k_{\omega}^2}{\lambda_m(4(k_{\omega} - \frac{1}{2}c_2) - 1)}、\quad k_{\omega} > \frac{1}{2}c_2 + \frac{1}{4} k R > λ m ( 4 ( k ω − 2 1 c 2 ) − 1 ) c 2 k ω 2 、 k ω > 2 1 c 2 + 4 1
定理1 : 設計条件を満たす場合、完全なシステムの追従誤差ダイナミクスの零平衡点はほぼ全域指数安定である。
単方向プロペラシステム : 最低7個のプロペラまたは可変傾斜機構が必要双方向プロペラシステム : 機械構造を簡略化するが、電力効率が低い制御方法 : 幾何学的PIDから非線形MPCへの進化DCMDモデル : ブラシレスDCモーターダイナミクスに基づく精密モデルTDモデル : 精度と複雑さのバランスを取った簡略化推力ダイナミクスモデル応用 : 固定翼機と四ローターでの成功した応用プロペラダイナミクスを考慮した全方向マルチロータ制御器の設計に成功 ほぼ全域指数安定性を理論的に証明 積極的機動における性能改善を実験的に検証 追加センサ不要の実用的設計 プロペラ飽和 : 現在の設計ではプロペラ飽和制約を考慮していない反転遅延 : 双方向プロペラの反転遅延がモデル化されていないパラメータ適応 : プロペラ特性変化に対するロバスト性が考慮されていない制約付き最適化 : 制約付き最小二乗法によるプロペラ飽和の処理反転遅延モデリング : 特に双方向プロペラの詳細なモデリングロバスト制御 : プロペラ特性変化に適応するロバスト制御メカニズムの統合理論的厳密性 : 完全な安定性分析と証明を提供実用性 : 追加センサ不要で実装が容易実験の充実 : 3種類の異なる積極的機動で検証革新性 : 全方向マルチロータにおけるプロペラダイナミクスの体系的考慮の初例モデル簡略化 : TDモデルはDCMDモデルと比べて精度が低い実験規模 : 1つのプラットフォームでのみ検証、多様性に欠ける環境制限 : 理想的な室内環境での実験、外部干渉テストが不足計算複雑度 : リアルタイム計算負荷の詳細分析が不足学術的貢献 : 全方向マルチロータ制御理論に新しい視点を提供実用的価値 : 積極的機動応用に実行可能なソリューションを提供再現性 : 詳細なパラメータと実装詳細を提供拡張性 : 後続研究に良好な理論的基礎を提供捜索救助 : 複雑な環境での積極的機動実行が必要産業検査 : 狭い空間での正確な位置決めと検査研究プラットフォーム : 全方向飛行体制御研究の基礎プラットフォーム娯楽表演 : 複雑な飛行軌跡が必要な表演応用本論文は34篇の関連文献を引用しており、マルチロータモデリングと制御、全方向飛行体設計、プロペラダイナミクスモデリングなど主要分野の重要な研究をカバーし、本研究に堅実な理論的基礎を提供している。
総合評価 : これは全方向マルチロータ制御分野における重要な貢献を有する高品質論文であり、理論的に厳密で、実験が充実しており、実用的価値と学術的影響力が強い。