Abstract:We summarize and numerically compare two approaches for modeling and simulating the dynamics of dry granular matter. The first one, the discrete-element method via penalty (DEM-P), is commonly used in the soft matter physics and geomechanics communities; it can be traced back to the work of Cundall and Strack [P. Cundall, Proc. Symp. ISRM, Nancy, France 1, 129 (1971); P. Cundall and O. Strack, Geotechnique 29, 47 (1979)GTNQA80016-850510.1680/geot.1979.29.1.47]. The second approach, the discrete-element method via complementarity (DEM-C), considers the grains perfectly rigid and enforces nonpenetration via complementarity conditions; it is commonly used in robotics and computer graphics applications and had two strong promoters in Moreau and Jean [J. J. Moreau, in Nonsmooth Mechanics and Applications, edited by J. J. Moreau and P. D. Panagiotopoulos (Springer, Berlin, 1988), pp. 1-82; J. J. Moreau and M. Jean, Proceedings of the Third Biennial Joint Conference on Engineering Systems and Analysis, Montpellier, France, 1996, pp. 201-208]. The DEM-P and DEM-C are manifestly unlike each other: They use different (i) approaches to model the frictional contact problem, (ii) sets of model parameters to capture the physics of interest, and (iii) classes of numerical methods to solve the differential equations that govern the dynamics of the granular material. Herein, we report numerical results for five experiments: shock wave propagation, cone penetration, direct shear, triaxial loading, and hopper flow, which we use to compare the DEM-P and DEM-C solutions. This exercise helps us reach two conclusions. First, both the DEM-P and DEM-C are predictive, i.e., they predict well the macroscale emergent behavior by capturing the dynamics at the microscale. Second, there are classes of problems for which one of the methods has an advantage. Unlike the DEM-P, the DEM-C cannot capture shock-wave propagation through granular media. However, the DEM-C is proficient at handling arbitrary grain geometries and solves, at large integration step sizes, smaller problems, i.e., containing thousands of elements, very effectively. The DEM-P vs DEM-C comparison is carried out using a public-domain, open-source software package; the models used are available online.

Classical and Non-Classical Kinematic Fields of Two-Dimensional Penetration Tests on Granular Ground by Discrete Element Method Analyses

Discrete Element Simulations of Drained Granular Material Response under Multidirectional Rotational Shear

Discrete Element Modeling of Sand Behavior in a Biaxial Shear Test

Undrained Cyclic Behavior of Granular Materials Considering Initial Static Shear Effect: Insights from Discrete Element Modeling

Assess the Effects of Loading Rate and Interface Roughness on T-bar Penetration Resistance

Modelling T-bar Penetration in Soft Clay Using Large-Displacement Sequential Limit Analysis

Progressive Failure Mechanism in Granular Materials Subjected to an Alternant Active and Passive Trapdoor

Numerical simulation of cone penetration test by using CFD–DEM coupled analysis

Discrete element model (DEM) simulations of cone penetration test (CPT) measurements and soil classification

Unified discrete‐element approach applying arbitrary undrained loading paths in element testing for granular soils

Discrete element modelling of thermal penetration test with heating and cooling

Discrete element modelling of rotary CPT and its applications

A numerical investigation on the effect of rotation on the Cone Penetration Test

Investigating projectile penetration into immersed granular beds via CFD-DEM coupling

On the Influence of Inter-Particle Friction and Dilatancy in Granular Materials: a Numerical Analysis

DEM numerical test method for granular matter under complex 3D loading

Mechanics of the Penetration and Filtration of Cement-Based Grout in Porous Media: New Insights from CFD–DEM Simulations

Assessing the Bearing Capacity of Backfills by Stress Wave Velocity and Cone Penetration Resistance

Experimental and Numerical Study of Rate Effect in Cone Penetration Tests

Nonparametric Terrain Estimation Based on the Interaction Simulation Between Planetary Penetrator and Soil

Compliant contact versus rigid contact: A comparison in the context of granular dynamics