Parallel implicit contact algorithm for soft particle systems

Numerical simulations of the compaction of assemblies of rubberlike particles: A quantitative comparison with experiments

“…Although, in these studies the global properties were quite appropriately estimated, the particle-scale parameters were not well predicted [7,8]. Recently some numerical works using methodologies combining features of the DEM with appropriate methods for the simulation of particle deformations were carried out [9][10][11][12].…”

Compaction Model for Highly Deformable Particle Assemblies

“…The compaction of deformable matter has been addressed in experiments using ceramic, metallic, and pharmaceutical powders [1][2][3][4][5], gels [6][7][8], rubberlike particles [9][10][11][12], and even blood cells [13]. More recently, developments of numerical approaches based on meshless methods [14][15][16], the discrete element method [17,18], or the coupled finite element-discrete element methods [11,[19][20][21][22][23][24] have enabled the exploration of the physics of deformable granular media.…”

“…Kim and Carroll [3] arrived at a similar relationship while studying metallic powders, but also introducing three constants yet to be defined by experiments and data fitting (Y 0 , μ, and γ). Nezamabadi et al [16] also introduced a model of the form P ∝ −ϕ lnðϕÞ employing an incremental formulation of the mean pressure and the macroscopic effective P-wave modulus, but standing, nonetheless, upon fitting constants (c 1 and c 2 ). Finally, it is also notable that Zhang et al [25], using a multiparticle finite element method, nicely reproduced the compaction curve after fitting a double logarithmic equation on their data.…”

“…Note that the connectivity Z of an assembly of deformable grains is closely linked to the evolution of solid fraction as denser particle configurations lead to higher coordination numbers. Additionally, it has been systematically reported in the literature that the relation Z − ϕ obeys a power law of the form ðZ − Z 0 Þ ¼ kðϕ − ϕ 0 Þ α , with exponent α ≃ 0.5 [11,16,44], and k ≃ 5.1 a structural parameter fully defined as the P=E → ∞, and ϕ and Z reach maximal values. We found the same proportionality in our experiments independently of the coefficient of friction (see Fig.…”

Compaction of elastic granular materials: inter-particles friction effects and plastic events

“…15 Another approach allowing taking into account the realistic behavior of soft particles, is based on an implicit Material Point Method (MPM) coupled with the Contact Dynamics (CD) method. [15][16][17] In the MPM, each particle is discretized by a set of material points and the implicit MPM formulation permits to couple efficiently with implicit modeling of unilateral contacts and friction between particles as in the CD method. 18,19 An alternative technique consists in combining the Finite Element Method (FEM) with DEM.…”