[論文レビュー] On the similarity of Peridynamics and Smooth-Particle Hydrodynamics

Georg Ganzenmüller, Stefan Hiermaier+1 more|arXiv (Cornell University)|Jan 31, 2014

Numerical methods in engineering参考文献 20被引用数 1

ひとこと要約

この論文は、歪み勾配を用いた古典的材料モデルにペルィダイナミクス（Peridynamics）を適用した場合、ノード内積分を用いると、その離散化された方程式がスムーズ・パーティクル・ハイドロダイナミクス（Smooth-Particle Hydrodynamics, SPH）のそれと数学的に同等であることを示している。この同等性は、ペルィダイナミクスがSPHの不安定性を引き継ぐことを示唆しており、特にノード内積分におけるランク不足に起因する誤ったゼロエネルギーモードが生じる。このため、安定かつ正確なシミュレーションを実現するには、強化された積分スキームの導入が不可欠である。

ABSTRACT

This paper discusses the similarity of meshless discretizations of Peridynamics and Smooth-Particle-Hydrodynamics (SPH), if Peridynamics is applied to classical material models based on the deformation gradient. We show that the discretized equations of both methods coincide if nodal integration is used. This equivalence implies that Peridynamics reduces to an old meshless method and all instability problems of collocation-type particle methods apply. These instabilities arise as a consequence of the nodal integration scheme, which causes rank-deficiency and leads to spurious zero-energy modes. As a result of the demonstrated equivalence to SPH, enhanced implementations of Peridynamics should employ more accurate integration schemes.

研究の動機と目的

• 古典的材料モデルに基づく歪み勾配を用いた場合のペルィダイナミクスとスムーズ・パーティクル・ハイドロダイナミクス（SPH）の間の数学的関係を調査すること。
• ペルィダイナミクスにおけるノード内積分が、SPHと同様の数値的不安定性を引き起こすかどうかを特定すること。
• ノード内積分を用いる場合、ペルィダイナミクスがメッシュレス法として簡略化され、既知の不安定性問題を有することを確立すること。
• これらの不安定性を克服するため、ペルィダイナミクスの定式化においてより正確な積分スキームを採用する動機を提示すること。

提案手法

• 研究では、歪み勾配に基づく古典的材料モデルを用いてペルィダイナミクスを定式化する。
• 弱形式のペルィダイナミクス方程式を離散化するためにノード内積分を適用する。
• 得られた離散方程式を、同一の積分スキームを用いたSPHのそれと直接比較する。
• 同一の仮定と積分規則の下で、離散系の同等性を解析的に証明する。
• ノード内積分スキームにおけるランク不足が、誤ったゼロエネルギーモードの原因であると特定する。
• 強化された積分技術——例えば安定化またはガウス積分——を採用すべきであると主張する。

実験結果

リサーチクエスチョン

• RQ1 ペルィダイナミクスとSPHがノード内積分と古典的材料モデルを用いる場合、数学的同等性が存在するか？
• RQ2 SPHで観察される不安定性、例えば誤ったゼロエネルギーモードは、ノード内積分を用いたペルィダイナミクスに対しても同様に生じるか？
• RQ3 ノード内積分を用いたペルィダイナミクスの定式化における不安定性の原因は何か？
• RQ4 SPHとの同等性は、ペルィダイナミクスの実装改善にどのように寄与できるか？
• RQ5 SPHから引き継ぐ不安定性を回避するには、ペルィダイナミクスでどの積分スキームを優先すべきか？

主な発見

• ペルィダイナミクスとSPHは、ノード内積分と歪み勾配に基づく古典的材料モデルを用いる場合、同一の離散方程式を生成する。
• ペルィダイナミクスの不安定性は、ノード内積分スキームにおけるランク不足に起因し、誤ったゼロエネルギーモードを引き起こす。
• 同等性は、ペルィダイナミクスがSPHのようなコロケーション型粒子法のすべての数値的不安定性を引き継ぐことを示唆する。
• これらの不安定性の根本的原因は、ノード内積分スキームにおける一貫性の欠如であり、適切な積分精度を保証できないことにある。
• 本論文は、安定性と正確性を確保するため、ペルィダイナミクスのシミュレーションにおいて、強化された積分スキーム——例えば安定化またはガウス求積——を用いる必要があると結論づける。
• 誤ったモードを回避し、ペルィダイナミクス解析の信頼性を確保するには、より正確な積分手法の採用が不可欠である。