[論文レビュー] Shape derivative of potential energy and energy release rate in fracture mechanics

Masato Kimura, Isao Wakano|ArXiv.org|Jun 20, 2008

Numerical methods in engineering参考文献 17被引用数 6

ひとこと要約

この論文は、バナッハ空間解析と陰関数定理を用いて、破壊力学におけるポテンシャルエネルギーの形状微分について、厳密な数学的枠組みを確立する。ドメイン変形におけるフレシェ微分として、エネルギー解放率の明示的積分公式を導出し、パラメータ化された変分原理とリプシッツ連続なドメイン摂動を用いて、曲がったクラックや滑らかでない領域への古典的結果の一般化を実現する。

ABSTRACT

We study general mathematical framework for variation of potential energy with respect to domain deformation. It enables rigorous derivation of the integral formulas for the energy release rate in crack problems. Applying a technique of the shape sensitivity analysis, we formulate the shape derivative of potential energy as a variational problem with a parameter. Key tools of our abstract theory are a new parameter variational principle and the classical implicit function theorem in Banach spaces.

研究の動機と目的

• ドメイン変形下でのポテンシャルエネルギーの形状微分の一般数学的枠組みの構築。
• 抽象的パラメータ変動を用いたバナッハ空間におけるエネルギー解放率の厳密な導出。
• エネルギー解放率の古典的結果を一般化された曲がったクラックおよび滑らかでない領域へと拡張すること。
• 変分的手法を通じて形状微分とJ積分の関係を確立すること。
• 2次元線形弾性を超えて、多次元、非線形、および半線形楕円型系へとエネルギー解放率の公式を一般化すること。

提案手法

• バナッハ空間におけるパラメータ化された変分問題として、ポテンシャルエネルギーの形状微分を定式化する。
• ドメイン摂動下での最小ポテンシャルエネルギーのフレシェ微分を導出するために、バナッハ空間における陰関数定理を適用する。
• 変分問題の解の存在および一意性を保証するために、ラックス＝ミルグラムの定理を用いる。
• 合成関数の微分法則と発散定理を用いて、ドメイン積分表現として形状微分を導出する。
• リプシッツクラスのドメイン摂動を取り扱うために、新しいパラメータ変分原理を導入する。
• 仮想クラック延長をパラメータ化されたベクトル場によるドメイン変形としてモデル化することで、クラック問題に理論を適用する。

実験結果

リサーチクエスチョン

• RQ1 破壊力学におけるエネルギー解放率を、ポテンシャルエネルギーの形状微分としてどのように厳密に導出できるか。
• RQ2 一般化された曲がったクラックを対象とした多次元線形および半線形系において、エネルギー解放率を導出するための抽象的数学的枠組みは何か。
• RQ3 ドメイン摂動下で、形状微分がどのようにバナッハ空間設定においてフレシェ微分として定式化できるか。
• RQ4 陰関数定理は、エネルギー解放率の公式を導出する際に果たす役割は何か。
• RQ5 導出された公式は、古典的J積分表現およびグリフィスの基準を滑らかでない領域および高次元設定へとどのように一般化するか。

主な発見

• 形状微分はバナッハ空間設定においてフレシェ微分として導出され、形状感度解析の厳密な基礎が提供される。
• エネルギー解放率の明示的積分公式が得られる：$ E_*'( ho_0)[\mu] = \int_{\partial\Omega} \left( \frac{1}{2}|\nabla u|^2 - fu \right) \mu \cdot \nu \, d\mathcal{H}^{n-1} - \int_{\partial\Omega} (\nabla u \cdot \nu)(\nabla u \cdot \mu) \, d\mathcal{H}^{n-1} $、これは特定の境界条件下で簡略化される。
• 境界条件 $ \Gamma_D = \partial\Omega $ かつ $ g \equiv 0 $ の場合、エネルギー解放率は $ E_*'\!(\rho_0)[\mu] = -\frac{1}{2} \int_{\partial\Omega} |\nabla u|^2 \mu \cdot \nu \, d\mathcal{H}^{n-1} $ に簡略化され、境界における変位勾配に依存することが示される。
• ベクトル場 $ \mu $ がディリクレ境界およびクラック表面から離れてコンパクトな台を持つ場合、公式は $ E_*'\!(\rho_0)[\mu] = -\frac{1}{2} \int_{\Gamma_D} |\nabla u|^2 \mu \cdot \nu \, d\mathcal{H}^{n-1} - \int_{\partial\Omega \setminus \Gamma_D} f u \, \mu \cdot \nu \, d\mathcal{H}^{n-1} $ に一般化される。
• エネルギー差の直接的な極限評価を避ける代わりに、抽象的パラメータ変動と陰関数定理を用いることで、より一般的かつ安定な枠組みで結果を導出する。
• リプシッツクラスのドメイン摂動を用いることで、一般化された曲がったクラックおよび滑らかでない領域へと、従来の2次元線形弾性における直線クラックに限られた結果を拡張する。