Physics-Informed Single Atom Matching Pursuit: Guided-Waves Wavenumbers and Propagation Distance Estimation for Damage Localization in Structural Health Monitoring

技術要約：誘導波損傷局在化のための物理情報に基づくシングルアトム・マッチングパース

問題提起
薄肉構造物の構造ヘルスモニタリング（SHM）は、多くの場合、圧電素子（PZT）によって生成される超音波誘導波（GW）に依存している。この領域における主要な課題は、単一のアクチュエータによって複数の分散モード（具体的には $A_0$ モードと $S_0$ モード）が同時に励起されることである。これらの波包は伝播する過程で境界や潜在的な損傷と相互作用し、重複および分散した寄与を含む複雑な信号を作り出す。従来の信号処理技術には、以下のような重大な限界がある：

• 辞書ベースの手法（例：標準的なマッチングパース）は、あらゆる可能なモードの寄与をあらかじめ計算した辞書を必要とするため、連続的な構造物に対しては非現実的であり、適応性を制限する。
• 非辞書ベースの手法（例：EMD、VMD、またはデータ駆動型オートエンコーダ）は、波の伝播に関する基礎的な物理学との強い結びつきを欠いていることが多く、抽出された特徴量は数学的には妥当であっても、物理的な解釈可能性に欠ける。

本研究で対処する核心的な困難は、構造固有の分散曲線に関する事前知識なしに、モード波数関数と伝播距離を直接抽出することを可能にする、物理的な波の伝播原理に基づいた、かつ計算効率の高い分解フレームワークを構築することである。

手法：物理情報に基づくシングルアトム・マッチングパース（PISAMP）
著者らは、測定された信号を外部励起の遅延および分散したインパルス応答としてモデル化する信号分解手法である、物理情報に基づくシングルアトム・マッチングパース（PISAMP）を提案している。この手法は以下の原理に基づいている：

1. 物理的定式化： 本手法は周波数領域で信号をモデル化する。測定された信号 $\hat{s}(\omega)$ は、$N$ 個のアトムの線形結合として近似される。ここで、各アトムは特定の距離 $d_i$ を伝播する特定の波数関数 $k_i(\omega)$ を持つ波包を表す：
   $$\hat{s}(\omega) \approx \sum_{i=1}^{N} \alpha_i(d_i) \hat{x}(\omega) e^{-j k_i(\omega) d_i}$$
   ここで、$\hat{x}(\omega)$ は励起、$\alpha_i$ は振幅（減衰をモデル化）、指数項は波数関数 $k_i(\omega)$ を通じて伝播距離と分散効果の両方を符号化する。

2. 反復最適化： この問題は、貪欲な増分アプローチを用いて解決される。各項 $i$ について、アルゴリズムは振幅 $\alpha_i$、伝播距離 $d_i$、および波数関数 $k_i(\omega)$ を、周波数領域における再構成誤差を最小化するように共同で最適化する。

   • 初期化： 波数には、非分散的な推定値（例：$S_0$ モード）を用いて初期化を行う。
   • 振幅と距離： 波数が前回のイテレーションから固定されているものとして、距離に関する誤差を最小化することで決定される。
   • 波数更新： 波数関数は、一次のテイラー線形化を用いて更新される。補正項 $\Delta k_i(\omega)$ は、滑らかさと数値的効率を確保するために、第二種チェビシェフ多項式を用いて近似される。

3. スケール校正と損傷局在化： 指数項におけるスケールの不定性（すなわち、$k_i$ と $d_i$ が逆比例して一定となる場合、$k_i d_i$ が一定となる現象）が重要な課題として特定されている。これを解決するために：

   • 本手法では、損傷のない構造におけるアクチュエータからセンサへの直接経路をリファレンスとして使用し、グローバルなスケーリング係数を定義することで、数値的な距離を物理単位に変換する。
   • 全信号と比較して振幅が極めて小さい場合が多いものの、局在化において極めて重要となる「初動到着」波包を特定するために、修正された誤差基準が導入される。
   • 物理的な距離が確立された後、損傷位置は楕円局在化を用いて決定される。損傷からアクチュエータ、および損傷からセンサへの総経路長は、2つのユークリッド距離の和としてモデル化される。複数のセンサからの推定経路長に最も適合する損傷座標を見つけるために、ガウス・ニュートン法またはレーベンバーグ・マルカート法を用いて非線形最小二乗問題が解かれる。

主な結果
論文では、数値的な調査を通じてPISAMP法の妥当性を検証している：

• 信号近似： 等方性板を用いた学術的なケースにおいて、PISAMPは6つのサブモードを用いて4%の再構成誤差を達成し、非分散型（$S_0$）と分散型（$A_0$）のモードを判別することに成功した。
• 複雑な構造物： エアバスA380ファンカウル構造（補強材を持つ複合材プレート）の実験データに適用したところ、本手法は40個のサブモードを用いて20%の再構成誤差を達成し、反射や材料減衰を効果的に処理した。
• 損傷局在化： 30通りの異なる損傷構成を持つ2次元正方形プレートのシミュレーションにおいて、提案されたアプローチ（PISAMPによる距離推定と楕円局在化の組み合わせ）は、x座標で平均相対誤差2.35%、y座標で**2.26%**を達成した。本手法は、散乱波成分のみを用いて損傷位置を特定することに成功した。

意義と主張
著者らは、本研究の主要な新規性は、コンパクトな信号分解フレームワークの中に波の伝播物理学を直接統合したことにあると主張している。純粋なデータ駆動型手法や計算負荷の高い手法とは異なり、PISAMPは強い物理的制約を組み込むことで、低次元かつ計算効率の高い表現を実現している。

本研究の意義は以下の二点に集約される：

1. 解釈可能性： すべてのモードを網羅する事前計算済みの辞書を必要とすることなく、物理的に意味のある特徴、具体的にはモード波数関数と伝播距離を直接特定することを可能にする。
2. 局在化能力： 分解された信号の特徴を幾何学的な経路長に直接結びつけることにより、薄板構造における堅牢で解釈可能な損傷局在化の手法を提供する。

論文では、現在の実装はユークリッド距離が適用される平面表面に限定されていることを謙虚に記している。今後の研究課題として、曲面ジオメトリ（測地線距離の計算が必要）への拡張、および実世界の実験データによる検証が挙げられている。