[论文解读] 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

一句话总结

本文通过巴拿赫空间分析与隐函数定理，为断裂力学中势能的形状导数建立了严格的数学框架。它在域变形下以弗雷chet导数形式推导出能量释放率的显式积分公式，通过参数化变分原理与利普希茨类域摄动，将经典结果推广至曲裂纹和非光滑域的情形。

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-积分之间的联系。
• 将能量释放率公式推广至多维、非线性及半线性椭圆系统，超越二维线性弹性理论。

提出的方法

• 将势能的形状导数表述为巴拿赫空间中的参数化变分问题。
• 在巴拿赫空间中应用隐函数定理，推导最小势能的弗雷chet导数。
• 利用拉克斯-米爾格拉姆定理确保变分问题解的存在性与唯一性。
• 通过链式法则与散度定理，推导形状导数的域积分表达式。
• 引入新的参数变分原理，以处理利普希茨类的域摄动。
• 通过将虚拟裂纹扩展建模为参数化向量场的域变形，将理论应用于裂纹问题。

实验结果

研究问题

• RQ1 如何将断裂力学中的能量释放率严格推导为势能的形状导数？
• RQ2 何种抽象数学框架可使一般曲裂纹在多维线性与半线性系统中实现能量释放率的推导？
• RQ3 在域摄动下，形状导数如何在巴拿赫空间设定中表述为弗雷chet导数？
• RQ4 隐函数定理在推导能量释放率公式中起到何种作用？
• RQ5 所推导的公式如何将经典J-积分表达式与格里菲斯判据推广至非光滑与高维情形？

主要发现

• 在巴拿赫空间设定下，势能的形状导数被推导为弗雷chet导数，为形状敏感性分析提供了严格的理论基础。
• 获得能量释放率的显式积分公式：$ 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_*'\left(\rho_0\right)[\mu] = -\frac{1}{2} \int_{\partial\Omega} |\nabla u|^2 \mu \cdot \nu \, d\mathcal{H}^{n-1} $，表明其对边界上位移梯度的依赖性。
• 当向量场 $ \mu $ 的紧支集与狄利克雷边界及裂纹面不相交时，公式推广为 $ E_*'\left(\rho_0\right)[\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} $。
• 该方法避免了对能量差值的直接极限计算，转而通过抽象参数变分与隐函数定理，在更一般且更稳定的框架下推导结果。
• 通过利普希茨类域摄动，该方法可适用于一般曲裂纹与非光滑域，扩展了以往仅限于二维线性弹性中直裂纹的研究结果。