Quasi-static loading of granular media as a linear complementarity problem
Abstract
The discrete quasi-static response of rate-independent dissipative granular media is addressed. Granular systems are conventionally simulated with methods that are intrinsically dynamic, such as the discrete element (DEM) and discontinuous deformation (DDA) methods, with the particles' accelerations and damping being essential aspects. In contrast, quasi-static methods derive from the static stiffness relationships among the particles. With frictional contacts, an assembly's stiffness is incrementally non-linear and dependent on the direction of loading. The paper resolves this difficulty by casting the response as a linear complementarity problem (LCP). The approach benefits from a foundation of past research on existence, uniqueness, and stability of LCP solutions, which are expounded in a concise set of rules. The LCP of a granular system is derived, accounting for geometric effects that arise from curvatures of particles at their contacts, frictional contact stiffnesses, and displacement constraints on the particles. The paper describes eight aberrant conditions (pathologies) of granular systems, in the forms of various bifurcations and instabilities. When placed in the context of an LCP, the eight conditions are unambiguously defined. These pathologies include three types of bifurcation: discrete, continuous-bounded (not yet reported in the literature), and continuous-unbounded. Not yet reported in the literature is an abrupt discontinuous change in movements upon a continuous change of the data. Methods and results are illustrated with examples of granular systems. Results show that instability and bifurcation are pervasive conditions near failure; that failure and softening are sensitive to the geometric contours of the particles at contacts; and that quasi-static systems can encounter states that require a dynamic transition to resolve a lack of paths consistent with loading.
Versions
Related
- Intruder dynamics in granular media under localized surface loading
- Intruder dynamics in granular media under localized surface loading
- Loading history and window geometry bound compact-state slip ranking during granular shear startup
- Shape matters: DEM investigation of geometry-controlled mechanical response in irregular rock fragments under static and dynamic loading
- Linear stability and viscoelastic response in weakly jammed frictional granular systems
BibTeX
@misc{kuhn2026quasistatic,
title = {{Quasi-static loading of granular media as a linear complementarity problem}},
author = {Matthew R. Kuhn},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.00152},
eprint = {2610.00152},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.00152},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Quasi-static loading of granular media as a linear complementarity problem
AU - Matthew R. Kuhn
PY - 2026
DA - 2026/10/01
PB - arXiv
DO - 10.48550/arxiv.2610.00152
UR - https://arxiv.org/abs/2610.00152
AB - The discrete quasi-static response of rate-independent dissipative granular media is addressed. Granular systems are conventionally simulated with methods that are intrinsically dynamic, such as the discrete element (DEM) and discontinuous deformation (DDA) methods, with the particles' accelerations and damping being essential aspects. In contrast, quasi-static methods derive from the static stiffness relationships among the particles. With frictional contacts, an assembly's stiffness is incrementally non-linear and dependent on the direction of loading. The paper resolves this difficulty by casting the response as a linear complementarity problem (LCP). The approach benefits from a foundation of past research on existence, uniqueness, and stability of LCP solutions, which are expounded in a concise set of rules. The LCP of a granular system is derived, accounting for geometric effects that arise from curvatures of particles at their contacts, frictional contact stiffnesses, and displacement constraints on the particles. The paper describes eight aberrant conditions (pathologies) of granular systems, in the forms of various bifurcations and instabilities. When placed in the context of an LCP, the eight conditions are unambiguously defined. These pathologies include three types of bifurcation: discrete, continuous-bounded (not yet reported in the literature), and continuous-unbounded. Not yet reported in the literature is an abrupt discontinuous change in movements upon a continuous change of the data. Methods and results are illustrated with examples of granular systems. Results show that instability and bifurcation are pervasive conditions near failure; that failure and softening are sensitive to the geometric contours of the particles at contacts; and that quasi-static systems can encounter states that require a dynamic transition to resolve a lack of paths consistent with loading.
ER -