[論文レビュー] Lattice Discrete Particle Model (LDPM): Comparison of Various Time Integration Solvers and Implementations

Erol Lale, Jan Eliáš+18 more|arXiv (Cornell University)|Mar 13, 2026

Numerical methods in engineering被引用数 0

ひとこと要約

本研究は、CPUとGPU実装の複数の暗黙的・明示的LDPMソルバを比較し、振動、圧縮、曲げ、張力、無拘束圧縮を含むベンチマーク試験での精度・安定性・効率性を評価します。公開可能なベンチマークデータと公開済み実装の参照も提供します。

ABSTRACT

This article presents a comparison of various implementations of the Lattice Discrete Particle Model (LDPM) for the numerical simulation of concrete and other heterogeneous quasibrittle materials. The comparison involves the use of transient implicit and explicit solvers and steady-state (static) solvers and implementations for Central Processing Unit (CPU) as well as Graphics Processing Unit (GPU). The various implementations are compared on the basis of a set of benchmarks tests describing behaviors of increasing computational complexity. They include elastic vibrations, confined strain-hardening compressive response, tensile fracture, and unconfined strain-softening compressive response. Metrics of interest extracted from the simulations include macroscopic stress versus strain responses, computational times, number of iterations, and energy balance error. Pairwise comparison of final crack patterns is provided through the correlation coefficient and normalized root mean square error of the crack opening vectors. Moreover, for the most numerically challenging case of unconfined compression with sliding boundary conditions, the stability of the strain-softening response is tested by perturbing the solutions as well as changing the convergence criteria and time step size. Attached to this paper is the complete input data of the benchmark tests; this will allow researchers to run the examples and compare them with their own implementations. In addition, most of the reported implementations are publicly available in open source packages.

研究の動機と目的

• 主要なソフトウェアプラットフォーム全体で、さまざまなLDPMソルバと実装の性能を動機づけ、定量化する。
• 増加する非線形複雑さの下で、LDPMの精度・収束性・エネルギーバランスを評価する。
• 今後のLDPM実装を検証・比較するためのベンチマークデータとガイダンスを提供する。

提案手法

• Facetベースのひずみとセル間相互作用を用いてLDPM支配方程式を定式化する。
• 複数のソフトウェアパッケージ（CPUとGPU）で明示・暗黙ソルバを用いたLDPMを実装する。
• 5つのベンチマーク試験（弾性振動、拘束圧縮、欠陥付き三点曲げ、直接張力、無拘束圧縮）を用いて応答・反復・エネルギーバランス・亀裂パターンを比較する。
• グローバル荷重-変位、応力-ひずみ応答、運動エネルギー、亀裂パターンの相関を用いて性能を評価する。
• 再現性のために完全な入力データを添付し、実装へのオープンソースアクセスを提供する。

実験結果

リサーチクエスチョン

• RQ1 暗黙的および明示的LDPMソルバは、複雑な準脆性挙動に対して精度と収束性でどのように比較されるか。
• RQ2 LDPMシミュレーションのさまざまな荷重シナリオに対して、どのソルバタイプと実装が最も効率的か。
• RQ3 CPUとGPUなど異なるソフトウェア環境がひずみ低減領域におけるLDPMの性能と安定性にどう影響するか。
• RQ4 提供されたベンチマークデータはツール間でLDPM実装の再現可能な検証を可能にするか。

主な発見

ID	ソフトウェア	ソルバタイプ	アルゴリズム	質量
AE	ABAQUS	Explicit	Cent. Difference	Lumped
JA	JAX-LDPM	Explicit	Cent. Difference	Lumped
CI	Chrono	Implicit	HHT	Consistent
OA	Open Academic Solver (OAS)	Implicit	Generalized-_alpha	Consistent
CA	CAST3M	Implicit	Static	N/A
MP	FEMultiPhys	Implicit	HHT	Lumped
JU	Julia LDPM	Implicit	Static, Arc-Length

• 暗黙ソルバは、多数の亀裂と複雑な非線形挙動を伴うLDPMにも有効に適用できる。
• 特定の非線形シナリオでは、収束性と時間積分に関連する注意点はあるものの、暗黙法は明示法より有利な場合がある。
• エネルギーバランス・反復回数・亀裂パターンの相関を実装間で詳細に比較している。
• 他研究者が再現・検証できるよう、完全なベンチマーク入力データセットを添付している。
• ほとんどの実装がオープンソースパッケージとして公開されており、再現性が高い。