[論文レビュー] A Mixed Finite Element Method for Multi-Cavity Computation in Incompressible Nonlinear Elasticity

Weijie Huang, Zhiping Li|arXiv (Cornell University)|Jan 8, 2018

Advanced Numerical Methods in Computational Mathematics被引用数 3

ひとこと要約

本稿では、非圧縮性非線形弾性における多重キャビテーション成長を計算するため、$Q_2$-$P_1$ および $P_2^+$-$P_1$ の等パrametric要素を組み合わせた混合有限要素法を提示する。この手法は、解析的にロックフリーかつ収束性が保証されており、数値実験により、負荷が減少する際に非対称から対称な配置に遷移するという、新たな分岐現象が、臨界伸び閾値 $ olambda_c olinebreak circ olinebreak circ (1.175, 1.180)$ で支配されることを明らかにした。

ABSTRACT

A mixed finite element method combining an iso-parametric $Q_2$-$P_1$ element and an iso-parametric $P_2^+$-$P_1$ element is developed for the computation of multiple cavities in incompressible nonlinear elasticity. The method is analytically proved to be locking-free and convergent, and it is also shown to be numerically accurate and efficient by numerical experiments. Furthermore, the newly developed accurate method enables us to find an interesting new bifurcation phenomenon in multi-cavity growth.

研究の動機と目的

• 非圧縮性非線形弾性材料における大変形と体積不変制約が大きな課題をもたらす中で、多重キャビテーション発生を安定にシミュレートするための数値的手法を開発すること。
• 特にキャビテーション付近で顕著に現れるロック効果を回避するため、標準的な有限要素法が非圧縮性弾性に適用される際に生じる問題を克服すること。
• 多重キャビテーション問題に対して、数値的安定性と収束性を両立する混合有限要素式を確立すること。
• 負荷条件の変化に伴う多重キャビテーション解の分岐挙動、特に対称的・非対称的キャビティ配置間の遷移を数値的に調査すること。
• 提案手法が、エネルギー駆動によるキャビティ形状選択を含む複雑な変形パターンを正確に捉えられることを、解析的・数値的根拠で示すこと。

提案手法

• 本手法は、変位と圧力に $Q_2$-$P_1$ を、変位場の高精度近似に $P_2^+$-$P_1$ を用いる2種類の等パrametric要素を組み合わせた混合有限要素式を採用する。
• 混合式では、体積不変制約 $\det \nabla \mathbf{u} = 1$ を弱く扱うために、圧力に類似したラグランジュ乗数 $p$ を導入し、強制的強制を回避する。
• 離散化された変分問題は、鞍点汎関数から導かれる非線形オイラー＝ラグランジュ方程式に、減衰ニュートン法を適用して解く。
• 先行研究における二パramトリック混合要素の手法を拡張することで、解析的にロックフリーであることを証明し、適切な正則性仮定の下で最適収束率を保証する。
• 保存エネルギー密度 $W(F) = \mu|F|^s + d(\det F)$ を用い、非圧縮性下で定数となる $d(\det F)$ が収束証明において重要な役割を果たす。
• 変形勾配およびその行列式の $L^2$ および $L^1$ ノルムの数値的収束性を検証し、最適収束率が得られていることを確認した。

実験結果

リサーチクエスチョン

• RQ1 多重キャビテーション解が非対称から対称な配置に分岐する臨界負荷パラメータ $\lambda_c$ は何か？
• RQ2 変形が増大するに従い、軸対称的および非軸対称的キャビティ解間のエネルギー差 $E_s - E_r$ はどのように変化するか？
• RQ3 非圧縮性非線形弾性下の多重キャビティ系において、あるキャビティが他のキャビティよりも大きく成長する要因は何か？
• RQ4 高応答異方性を示す多重キャビティ問題に対し、混合有限要素法が同時にロックフリーかつ収束性を有できるか？
• RQ5 初期条件が多重キャビティ成長における最終的なキャビティ形状に与える影響は何か？また、エネルギー最小化は安定配置の選択にどのような役割を果たすか？

主な発見

• 提案された混合有限要素法は、解析的にロックフリーかつ収束性が保証されており、$L^2$ および $L^1$ ノルムの両方において、数値的に最適収束率が確認された。
• 数値実験により、$\lambda > 1.18$ の場合、軸対称的解と右側優位解のエネルギー差 $E_s - E_r$ が急激に増加し、非対称配置が好まれることを示した。
• 非軸対称解における2つのキャビティの体積比 $v_2/v_1$ は、$\lambda > 1.18$ の範囲で急激に増大し、キャビティ成長における強い非対称性を示した。
• 臨界遷移点 $\lambda_c$ は区間 $(1.175, 1.180)$ 内に位置し、負荷が減少するに従い、非対称解が好まれる状態から対称的で分離された配置に移行する。
• 本手法は、理論的予測（Henao & Serfatyによる多重キャビティ挙動）に整合するエネルギー的に好ましい配置の自発的選択を的確に捉えた。
• 本研究では、新たな分岐現象を発見した：負荷パラメータ $\lambda$ が $\lambda_c$ を下回る際、複数の安定な形態（非対称または融合したキャビティ）から、対称的かつ明確に分離された配置へと系が遷移する。