Enhanced input stacking for non-square MIMO modal identification of aeronautical structures via Fast and Relaxed Vector Fitting

技術概要：航空機構造の非正方形 MIMO モダル同定のための高速かつ緩和されたベクトルフィット法による強化入力スタッキング

問題定義
構造モダルパラメータ（固有振動数、減衰比、モード形状）の正確な同定は、航空機の認証、有限要素モデルの更新、および空力弾性安定性の予測にとって不可欠である。地上振動試験（GVT）は必要な実験データを提供するが、現代の航空宇宙構造は、結合したおよび密接に配置された振動モードを捉えるために、多入力多出力（MIMO）フレームワークを必要とすることが多い。

高速かつ緩和されたベクトルフィット（FRVF）アルゴリズムの適用には重大な制限が存在する。FRVF は元々電気回路網向けに開発された計算効率の高い周波数領域システム同定手法であるが、機械システムへの適用には課題がある。標準的な FRVF 実装は、入力チャネル数と出力チャネル数が等しい正方形の周波数応答関数（FRF）行列を前提としている。しかし、航空機の実験的モダル解析では、測定された応答チャネル数（センサー）が通常、励振位置数（アクチュエータ）を大幅に上回るため、非正方形の FRF 行列が生じる。この不一致により、既存の FRVF 定式化を現実的な GVT シナリオに直接適用することが不可能となっている。

手法
著者らは、強化された入力スタッキング戦略を通じて非正方形 MIMO 構成に対応する FRVF の再定式化を提案する。手法は 3 つの段階で進行する。

1. 強化入力スタッキング：核心的な革新は、3 次元 FRF テンソル（$N_{out} \times N_{in} \times N_{freq}$）を FRVF と互換性のある 2 次元配列に再構成することである。正方形システムに用いられる従来の上三角スタッキングに代わり、提案アルゴリズムは入力チャネルの周波数ごとの重ね合わせを行う。各出力チャネルについて、すべての入力チャネルからの FRF が加算される。
   • 理論的基礎：線形時不変システムのモダル重ね合わせの原理に基づき、この加算は仮想の一貫した励振によって駆動される実効 FRF を生成する。これは残留重み付けを変更するが、システム極（固有振動数および減衰比）を保持し、モード形状の空間的比率を維持するため、入力と出力間の物理的相関を変化させることなく、弱く励起されたモードの信号対雑音比（SNR）を向上させる。
2. システム同定：スタックされた FRF は FRVF アルゴリズムを用いて処理される。これには、QR 法に基づく残留推定を用いた反復有理近似、ノイズ処理のための弱い逆数重み付け、および効率的な極安定化が含まれる。
3. パラメータ抽出：得られた有理モデルから安定化ダイアグラムを用いてモダルパラメータを抽出する。候補極は、振動数の安定性（1 Hz 未満の変動）、減衰の一貫性（0.05 未満の変動）、およびモード形状の相関（モーダル・アシュアランス・クリテリオン > 0.95）に基づいて選択される。

主要な貢献

• FRVF の一般化：本作業は、FRVF アルゴリズムを正方形システムから任意の非正方形 MIMO システムへ拡張し、電気回路網モデリング手法を構造力学に応用する際の重要なギャップを埋める。
• 入力スタッキングアルゴリズム：複数の励振寄与を集約し、相互性仮定や同置センサー/アクチュエータを必要とせずにモードの観測性を向上させる特定のスタッキング手順（アルゴリズム 1）の導入。
• 検証フレームワーク：本手法は、カンチレバー梁の数値シミュレーションと実機航空機の実験的試験の両方を通じて検証された。

結果

• 数値的検証（カンチレバー梁）：
  • 本手法は、ノイズのない条件下で、模擬された梁の 10 個の全モードを同定することに成功し、振動数誤差は無視できるレベル（<0.1%）であり、モード形状の相関は完全（MAC = 1）であった。
  • ノイズ耐性：1% までの加算白色ガウスノイズ（AWGN）下でも、本手法は高い精度を維持した。3〜5% のノイズレベルであっても、固有振動数は正確（誤差 <0.9%）であり、モード形状は高い相関（MAC ≥ 0.99）を示したが、高次モードの減衰比推定値は誤差が増大し（最大 40%）、見られた。
• 実験的検証（BAE システムズ ホーク T1A）：
  • 本手法は、ホーク T1A 練習機（5 入力、91 出力）からの GVT データに、5〜165 Hz の帯域幅で適用された。
  • iLF および LSCE との比較：FRVF は、改良されたロエウナーフレームワーク（iLF）が報告した 32 モードのうち 27 モードを同定し、振動数偏差は一般的に 1% 以内であった。一方、最小二乗複素指数法（LSCE）は、振動数および減衰の両方で顕著な偏差を伴う 4 つのモードのみを同定した。
  • モード形状：視覚的検査および MAC 値により、FRVF が物理的に意味のあるモード形状、特に低周波で分離の良いモードを回復したことが確認された。40 Hz 以上のモードでは、高モダル密度および観測性の低下により性能がわずかに低下したが、これは混雑した帯域における等間隔極初期化の制限と一致する。

意義と主張
本論文は、強化された入力スタッキングを伴う拡張 MIMO 定式化の FRVF が、航空機構造のモダル解析に対する信頼性が高く、正確で、ノイズに耐性のあるアプローチを提供すると主張している。著者らは、このアプローチが以下の点を実現すると述べている。

• 標準 FRVF の正方形システム制限を克服し、センサー数がアクチュエータ数を上回る現実的な航空機 GVT キャンペーンに適した手法とする。
• 確立された広帯域手法である iLF と同等の性能を示しつつ、複雑な MIMO 環境において古典的な時間領域手法である LSCE よりも優れた性能を発揮する。
• 測定ノイズに対する高い堅牢性が求められるアプリケーションに適した、構造力学の同定のための計算効率の高い代替手段を提供する。

著者らは、本手法が非常に効果的である一方で、将来の研究では、高度に混雑した周波数帯域における収束をさらに改善するために、自動モード選択基準および動的極配置に焦点を当てることができると結論付けている。