Direct stress imaging from shear wave propagation

技術要約：せん断波伝播からの直接応力イメージング

問題提起
不均質な固体における応力場の定量的イメージングは、力学およびバイオメカニクスにおいて依然として大きな課題である。変位とは異なり、応力はテンソル状の状態変数であり、直接観測することはできない。従来のアプローチでは、通常、構成モデルと材料特性および基準構成に関する仮定を用いて、測定された変形から応力を間接的に推定する。これらの手法は、広範な材料校正を必要とし、軟らかい材料や不均質な材料、あるいは生物学的に関連のある材料への適用に苦慮することが多い。さらに、既存の非破壊技術（既存の応力が波の伝播を変化させるアコーストエラスティック効果に基づくもの）は、一般に局所的な波速測定、校正された高次弾性定数、または幾何学的形状や均質性に関する簡略化された仮定に依存している。その結果、複雑な材料における内部応力場の高解像度マッピングは、長らく困難なままであった。

手法：アコーストエラスティック・イメージング (AEI)
著者らは、構成モデルの明示的な指定や材料パラメータの校正を必要とせずに、せん断波伝播から直接応力場を再構成するために設計された非破壊フレームワークである、アコーストエラスティック・イメージング（AEI）を導入する。

• 理論的基礎: この手法は、既存の応力が材料の増分弾性応答を変化させるというアコーストエラスティック効果を利用している。前応力状態における微小変位の支配的な波動方程式には、応力状態に直接結びついた弾性係数が含まれている。具体的には、主座標系において、主応力の差（$\sigma_1 - \sigma_2$）は、オイラー弾性テンソルの特定の成分の差（$A^0_{1'2'1'2'} - A^0_{2'1'2'1'}$）に符号化されている。
• 逆問題の定式化: AEIは、応力回復を逆問題として扱う。本フレームワークは、材料パラメータを推定する代わりに、全視野の波の測定値から、支配的な波動方程式の空間的に変化する係数を再構成することを目指す。
• 物理情報に基づくディープラーニング: この逆問題を効率的に解くために、著者らは物理情報に基づくニューラルネットワーク（PINN）によって可能となる全せん断波形反転（FSWI）を採用している。
  • アーキテクチャは、未知の場を空間の連続関数として表現するために、6つの完全結合ニューラルネットワークを利用する。1つのネットワークは（非圧縮性を満たすためにストリーム関数を介してパラメータ化された）変位場とラグランジュ乗数を予測し、残りのネットワークは主応力の方向と波動方程式の係数を推論する。
  • 損失関数は、データ駆動型の項（予測された変位場と測定された変位場の間のミスマッチを最小化する）と、物理情報に基づく項（支配的なアコーストエラスティック波動方程式を強制する）を組み合わせたものである。
  • 反転プロセスは、残差を最小化するように、応力差（$\sigma_1 - \sigma_2$）と主方向（$\theta$）を反復的に更新する。

主な貢献と結果
著者らは、数値シミュレーションと実験的な超音波せん断波エラストグラフィ（SWE）測定を通じてAEIを検証している。

1. 非一様および残留応力の再構成:

   • 拘束された熱膨張によって生成される、空間的に非一様な残留応力を持つ材料を含む数値シミュレーションにおいて、AEIは応力の大きさと主方向の両方を正常に再構成した。
   • 本手法は、応力と材料の剛性を分離する能力を示した。硬い包含物を含む軟らかいマトリックス内の不均質材料において、AEIは弾性的不均質性にもかかわらず応力場を正確に回復した。対照的に、従来のSWEは、応力による波の変化を剛性の変化として解釈し、滑らかな弾性率マップを生成した。
   • 個々の弾性係数（$\alpha$ および $\gamma$）の再構成は、一様なオフセットを伴って非一意的であることが判明したが、それらの差（応力を決定するもの）は一意に回復可能であった。

2. サブ波長分解能:

   • 特筆すべき発見は、サブ波長空間分解能（< 0.28$\lambda$）の達成である。全波場（近接場成分や波面の歪みを含む）を利用することで、AEIはせん断波長よりも著しく小さい応力特徴を解像する。
   • この分解能は、励起周波数（6 kHz と 12 kHz の比較）にほとんど依存しないことが示された。これは、回折によって制限される従来の遠方場イメージングとは対照的である。ただし、著者らは、部分的な波場データ（例：単一の変位成分）のみを使用すると分解能が低下することを指摘しているが、より高い周波数が部分的に補償できる場合もある。

3. 実験的検証:

   • 本フレームワークは、PVAハイドロゲル試料を用い、標準的な超音波SWEプラットフォームを使用して実装された。せん断波はPZTアクチュエータによって生成され、超高速超音波イメージングによって記録された。
   • 波場における顕著な歪みはわずかであったが（平均速度が 4.3 m/s から 4.5 m/s へ増加）、AEIは有限要素参照解と定量的に一致する応力場と主方向を再構成した。本手法は実験ノイズに対して頑健であり、信号対雑音比（SNR）が約 30 dB を超えても精度を維持した。

意義と主張
著者らは、AEIを、運動学的イメージング（変位の測定）からテンソル状の状態変数イメージング（応力の測定）へとパラダイムを転換する、高解像度応力イメージングのための一般的なフレームワークとして位置づけている。

• 直接測定: 本研究は、応力が特定の構成則に依存せず、波の伝播に直接符号化された測定可能な機械的状態変数として扱えることを確立している。
• 相補性: AEIと従来のエラストグラフィは相補的なモダリティとして提示されている。エラストグラフィが材料の剛性を特徴付けるのに対し、AEIは材料の機械的な状態（応力）を特徴付ける。
• 適用可能性: 本手法は、複雑な材料や生物組織における内部機械状態のマッピングに向けた道筋を提供する。
• 限界: 著者らは、現在の限界として、比較的滑らかな応力変化の仮定や、簡略化された波場表現の使用を控えめに認めている。また、実験的検証が現在は制御された環境に限定されていること、および圧縮性媒体やより複雑な非線形性へとフレームワークを拡張することが今後の必要なステップであることを述べている。

要約すると、本論文は、全波形反転と物理情報に基づく学習を組み合わせることで、事前の材料校正なしに高空間分解能で内部応力場を再構成することが可能であることを示しており、軟質材料や生物組織の機械的挙動を理解するための新しいツールを提供している。