[論文レビュー] Normal-fault stress and displacement through finite-element analysis

Antonietta Megna, S. Barba+1 more|ArXiv.org|Oct 17, 2005

earthquake and tectonic studies参考文献 15被引用数 8

ひとこと要約

本稿では、繰り返し力補正法を用いて正断層の応力および変位場をモデル化する有限要素法を提示する。解析解との間に良好な一致を示し、変位差は5%未満、方位誤差は15°未塔、応力誤差は35%未塔であり、地球力学的解析における滑り分布の滑らかな再現に有効であることが示された。

ABSTRACT

We compute displacement and stress due to a normal fault by means of two-dimensional plane-strain finite-element analysis. To do so, we apply a system of forces to the fault nodes and develop an iterative algorithm serving to determine the force magnitudes for any slip distribution. As a sample case, we compute the force magnitudes assuming uniform slip on a 10-km two-dimensional normal fault. The numerical model generates displacement and stress fields that compare well to the analytical solution. In fact, we find little difference in displacements (<5%), displacement orientation (<15 DEG), and stress components (<35%, half of which due to slip tolerance). We analyze such misfit, and discuss how the error propagates from displacement to stress. Our scheme provides a convenient way to use the finite-elements direct method in a trial-and-error procedure to reproduce any smooth slip distribution.

研究の動機と目的

• 正断層のすべりに起因する応力および変位場を正確にシミュレートするための数値的手法を開発すること。
• 数値モデルにおいて所望のすべり分布を再現するために、断層面上に適切な引っ張り力（トレクション力）を適用する課題に対処すること。
• 2次元正断層における一様なすべり分布に対して、解析解と比較して有限要素モデルを検証すること。
• 数値解と解析解の乖離を定量的に評価し、特に変位から応力への誤差伝播のメカニズムに注目すること。
• 任意の滑らかなすべり分布を再現するために、直接的な有限要素法を実用的に用いるための反復的フレームワークを提供すること。

提案手法

• 10kmの長さを有する正断層を一様なすべりでモデル化するため、2次元平面応力状態の有限要素解析を用いる。
• 断層面上にノード力の体系を適用してすべりを模擬し、力の大きさを反復的に決定する。
• 反復アルゴリズムを用いて力の大きさを調整し、得られるすべり分布が目標（例：一様すべり）と一致するようにする。
• 計算された変位場および応力場を解析解と比較して検証する。
• 変位の大きさ、方向、応力成分の差異を分析することで、誤差要因を評価する。
• 特にすべり許容誤差を考慮した場合に、変位誤差が応力計算にどのように伝播するかを評価する。

実験結果

リサーチクエスチョン

• RQ1 正断層に一様なすべりが生じる場合、有限要素モデルは変位および応力場をどの程度正確に再現できるか？
• RQ2 数値解と解析解との間に生じる誤差の大きさは、変位、変位の向き、応力成分においてそれぞれどの程度か？
• RQ3 変位の誤差は、計算された応力場にどのように影響するか？
• RQ4 繰り返し力補正法は、断層上に所望の滑らかなすべり分布を効果的に再現できるか？
• RQ5 このモデル化フレームワークにおいて、数値解と解析解との間で顕著な乖離を生じる主な要因は何か？

主な発見

• 有限要素モデルは、解析解との間で変位場の差が大きさにおいて5%未塔である。
• 変位の向きに関する誤差は15°未塔であり、数値解と解析解との間で良好な角度的一致が得られている。
• 応力成分の差は35%未塔であり、その約半分は数値設定におけるすべり許容誤差に起因している。
• 繰り返し力補正スキームは滑らかなすべり分布を効果的に再現でき、直接的な有限要素法を断層モデル化に活用可能であることを示している。
• 誤差伝播解析から、変位の不正確さが応力場の乖離に顕著に寄与していることが判明し、応力計算が変位誤差に極めて敏感であることが強調された。
• 本モデルは解析解と強く一致しており、正断層系の地球力学的シミュレーションへの応用が妥当であることが検証された。