Research

Paper

TESTING March 13, 2026

Lattice Discrete Particle Model (LDPM): Comparison of Various Time Integration Solvers and Implementations

Authors

Erol Lale, Jan Eliáš, Ke Yu, Matthew Troemner, Monika Středulová, Julien Khoury, Tianju Xue, Ioannis Koutromanos, Alessandro Fascetti, Bahar Ayhan, Baixi Chen, Giovanni Di Luzio, Yuhui Lyu, Madura Pathirage, Gilles Pijaudier-Cabot, Lei Shen, Alessandro Tasora, Lifu Yang, Jiawei Zhong, Gianluca Cusatis

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.

Metadata

arXiv ID: 2603.13190
Provider: ARXIV
Primary Category: cs.CE
Published: 2026-03-13
Fetched: 2026-03-16 06:01

Related papers

Raw Data (Debug)
{
  "raw_xml": "<entry>\n    <id>http://arxiv.org/abs/2603.13190v1</id>\n    <title>Lattice Discrete Particle Model (LDPM): Comparison of Various Time Integration Solvers and Implementations</title>\n    <updated>2026-03-13T17:24:51Z</updated>\n    <link href='https://arxiv.org/abs/2603.13190v1' rel='alternate' type='text/html'/>\n    <link href='https://arxiv.org/pdf/2603.13190v1' rel='related' title='pdf' type='application/pdf'/>\n    <summary>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.</summary>\n    <category scheme='http://arxiv.org/schemas/atom' term='cs.CE'/>\n    <category scheme='http://arxiv.org/schemas/atom' term='cond-mat.mtrl-sci'/>\n    <published>2026-03-13T17:24:51Z</published>\n    <arxiv:comment>28 pages, 15 tables, 8 figures</arxiv:comment>\n    <arxiv:primary_category term='cs.CE'/>\n    <arxiv:journal_ref>International Journal for Numerical and Analytical Methods in Geomechanics, 2026</arxiv:journal_ref>\n    <author>\n      <name>Erol Lale</name>\n    </author>\n    <author>\n      <name>Jan Eliáš</name>\n    </author>\n    <author>\n      <name>Ke Yu</name>\n    </author>\n    <author>\n      <name>Matthew Troemner</name>\n    </author>\n    <author>\n      <name>Monika Středulová</name>\n    </author>\n    <author>\n      <name>Julien Khoury</name>\n    </author>\n    <author>\n      <name>Tianju Xue</name>\n    </author>\n    <author>\n      <name>Ioannis Koutromanos</name>\n    </author>\n    <author>\n      <name>Alessandro Fascetti</name>\n    </author>\n    <author>\n      <name>Bahar Ayhan</name>\n    </author>\n    <author>\n      <name>Baixi Chen</name>\n    </author>\n    <author>\n      <name>Giovanni Di Luzio</name>\n    </author>\n    <author>\n      <name>Yuhui Lyu</name>\n    </author>\n    <author>\n      <name>Madura Pathirage</name>\n    </author>\n    <author>\n      <name>Gilles Pijaudier-Cabot</name>\n    </author>\n    <author>\n      <name>Lei Shen</name>\n    </author>\n    <author>\n      <name>Alessandro Tasora</name>\n    </author>\n    <author>\n      <name>Lifu Yang</name>\n    </author>\n    <author>\n      <name>Jiawei Zhong</name>\n    </author>\n    <author>\n      <name>Gianluca Cusatis</name>\n    </author>\n    <arxiv:doi>10.1002/nag.70286</arxiv:doi>\n    <link href='https://doi.org/10.1002/nag.70286' rel='related' title='doi'/>\n  </entry>"
}