[論文レビュー] On Mixed Isogeometric Analysis of Poroelasticity

Yared W. Bekele, Eivind Fonn+3 more|arXiv (Cornell University)|Jun 5, 2017

Advanced Numerical Analysis Techniques被引用数 5

ひとこと要約

本稿では、微小時間刻みにおける不適切な圧力振動を軽減するために、変位と孔内圧力に異なる多項式次数を用いる混合アイソジオメトリック定式化を提示する。この手法は、等次数の有限要素法と比較して、特に高次形状関数および材料界面でのC0連続性を有する場合、振動を顕著に低減し、一次元のテルツァギの問題および低透水性土層を含む層状構造モデルの両方において、より優れた安定性と収束性を示す。

ABSTRACT

Pressure oscillations at small time steps have been known to be an issue in poroelasticity simulations. A review of proposed approaches to overcome this problem is presented. Critical time steps are specified to alleviate this in finite element analyses. We present a mixed isogeometric formulation here with a view to assessing the results at very small time steps. Numerical studies are performed on Terzaghi's problem and consolidation of a layered porous medium with a very low permeability layer for varying polynomial degrees, continuities across knot spans and spatial discretizations. Comparisons are made with equal order simulations.

研究の動機と目的

• 微小時間刻みにおける有限要素法における既知の問題である、多孔性弾性解析における圧力振動を是正すること。
• 変位場と圧力場に異なる多項式次数を用いる混合アイソジオメトリック定式化の性能を評価すること。
• 多項式次数、 knot スパン間の連続性、メッシュの細分化が数値的安定性に与える影響を調査すること。
• ベンチマーク問題としての多孔性弾性問題において、混合次数アイソジオメトリック法と等次数定式化を比較すること。
• 低透水性層に近接する領域での振動低減に、メッシュのグレーディングがどの程度効果を発揮するかを評価すること。

提案手法

• NURBSに基づく形状関数を用い、変位と孔内圧力に異なる多項式次数を適用する混合アイソジオメトリック定式化を採用する。
• ガラーキン法を用いて、ビオの定着方程式の弱形式を一貫性のある積分法で実装する。
• 時間方向離散化に後退オイラー法を適用し、臨界値未満の時間刻みを設定して振動を誘発する。
• 材料界面における連続性を、knot スパン間でCpp−1またはC0連続性として強制し、材料特性の不連続性をモデル化する。
• 均一メッシュおよびグレーディングメッシュを用い、低透水性層に近接する領域を高精度に細分化して急峻な勾配を捉える。
• テルツァギの一次元定着問題および低透水性中間層を有する2次元層状土壌モデルに対して数値的検証を実施する。

実験結果

リサーチクエスチョン

• RQ1 微小時間刻みにおける孔内圧力の振動抑制に関して、混合次数アイソジオメトリック定式化は等次数定式化と比較してどの程度優れているか？
• RQ2 多項式次数の増加が孔内圧力解の安定性および精度に与える影響は何か？
• RQ3 材料界面における連続性（Cpp−1 対 C0）が多孔性弾性解析における振動挙動に与える影響は何か？
• RQ4 材料界面付近でのメッシュグレーディングが、低透水性領域における振動をどの程度低減するか？
• RQ5 混合アイソジオメトリック手法は、複雑な層状多孔性弾性媒体において振動を効果的に低減できるか？

主な発見

• 微小時間刻みにおいて、混合次数アイソジオメトリックシミュレーションは等次数定式化と比較して、圧力振動が顕著に低減された。
• 圧力の多項式次数を高めると振動が減少し、孔内圧力誤差の相対的L2ノルムは高次近似により収束した。
• 材料界面におけるC0連続性はCpp−1連続性よりも優れた結果を示し、特に低透水性層で顕著であった。
• 界面近傍で細分化されたグレーディングメッシュは、等次数および混合次数の両ケースにおいて振動を著しく低減した。
• Haga問題において、pp = 4およびC0連続性を有する混合次数シミュレーションでは、ほぼ振動が観察されず、高い安定性を示した。
• 結果から、材料の不連続性に基づく適応的メッシュ細分化戦略が、多孔性弾性解析における数値的性能を顕著に向上させ得ると示唆された。