[論文レビュー] A Simple Multi-Directional Absorbing Layer Method to Simulate Elastic Wave Propagation in Unbounded Domains

Jean‐François Semblat, Luca Lenti+1 more|arXiv (Cornell University)|Sep 6, 2010

Geotechnical Engineering and Underground Structures参考文献 41被引用数 4

ひとこと要約

本稿では、既存の有限要素ソフトウェアに容易に統合可能な、レイノルズ/コーガー減衰を用いたシンプルで効率的な多方向吸収層手法を提案する。この手法は、さまざまな波型および入射角において不自然な反射を効果的に低減し、2次元の不均一なモデルにおいてPMLと同等の性能を示す。

ABSTRACT

The numerical analysis of elastic wave propagation in unbounded media may be difficult due to spurious waves reflected at the model artificial boundaries. This point is critical for the analysis of wave propagation in heterogeneous or layered solids. Various techniques such as Absorbing Boundary Conditions, infinite elements or Absorbing Boundary Layers (e.g. Perfectly Matched Layers) lead to an important reduction of such spurious reflections. In this paper, a simple absorbing layer method is proposed: it is based on a Rayleigh/Caughey damping formulation which is often already available in existing Finite Element softwares. The principle of the Caughey Absorbing Layer Method is first presented (including a rheological interpretation). The efficiency of the method is then shown through 1D Finite Element simulations considering homogeneous and heterogeneous damping in the absorbing layer. 2D models are considered afterwards to assess the efficiency of the absorbing layer method for various wave types and incidences. A comparison with the PML method is first performed for pure P-waves and the method is shown to be reliable in a more complex 2D case involving various wave types and incidences. It may thus be used for various types of problems involving elastic waves (e.g. machine vibrations, seismic waves, etc).

研究の動機と目的

• 弾性波伝播の数値シミュレーションにおける人工境界での不自然な波の反射を解決すること。
• 複雑で不均一な、または層状の媒体において、従来の吸収境界条件や完全整合PML（PML）の限界を克服すること。
• 既存のレイノルズ／コーガー減衰式を活用できる、標準的な有限要素ソフトウェアと互換性のある手法を開発すること。
• 2次元構成において、P波・S波を含む多様な波型および入射角に対して、安定した性能を発揮すること。
• 純粋なP波と複雑なマルチ波状況の両方において、PMLと比較して手法の妥当性を検証すること。

提案手法

• コーガー減衰式を多方向吸収層として適応し、出射する弾性波を減衰させる。
• 計算領域の周囲に減衰ゾーンとして吸収層を配置し、空間的に変化する減衰パラメータを適用する。
• 粘弾性材料としての rheological 解釈を用いて、周波数依存減衰を持つ吸収層をモデル化する。
• 既存のレイノルズ減衰機能を活用することで、標準的な有限要素フレームワーク内に手法を実装する。
• 複数の方向および波型にわたる波の反射を最小限に抑えるために、層内の減衰パラメータを最適化する。
• 異なる材料の不均一性と波の入射角を想定した1次元および2次元の有限要素シミュレーションを通じて、手法を検証する。

実験結果

リサーチクエスチョン

• RQ1 Caugheyに基づく吸収層は、不均一領域における弾性波シミュレーションで不自然な反射を効果的に抑制できるか？
• RQ2 均一および不均一な減衰分布を有する1次元および2次元モデルにおいて、この手法はどのように性能を示すか？
• RQ3 純粋なP波に対して、広く用いられるPMLと比較して、この手法の反射抑制性能はいかがなものか？
• RQ4 複数の波型（P波およびS波）と斜め入射を含む複雑な波干渉を処理できるか？
• RQ5 地震波モデリングや機械振動解析などの実世界の応用において、この手法は実用的か？

主な発見

• コーガー吸収層手法は、均一および不均一な減衰を有する1次元シミュレーションにおいて、不自然な反射を顕著に低減した。
• 2次元モデルにおいて、P波およびS波を含む多様な波型および入射角（斜め入射を含む）に対して、反射を効果的に抑制した。
• 純粋なP波ケースにおいて、この手法はPMLと同等の反射レベルを達成し、その正確性を裏付けた。
• 複数の波モードと不均一な材料を含む複雑な2次元構成においても、安定した性能を維持した。
• 計算的に効率的であり、標準的なレイノルズ減衰機能を活用することで、既存の有限要素ソフトウェアへの直接的な実装が可能である。
• 吸収層の rheological 解釈により、物理的解釈が明確になり、パラメータチューニングが容易になった。