Three-dimensional granular contact dynamics with rolling resistance

A three-dimensional granular contact dynamics scheme is presented. The scheme is variational in structure, thus making it possible to solve the governing equations by means of mathematical programming methods. To facilitate the modeling of natural grains using spherical geometries, a rolling resistance model is developed. A number of static and dynamic benchmark examples are considered including the granular column collapse problem where the agreement between simulation and previously published experimental results is found to be very good.

[1]  Mihai Anitescu,et al.  A computational study of the use of an optimization-based method for simulating large multibody systems , 2009, Optim. Methods Softw..

[2]  Jean-Pierre Vilotte,et al.  Spreading of a granular mass on a horizontal plane , 2004 .

[3]  Tami C. Bond,et al.  Export efficiency of black carbon aerosol in continental outflow: Global implications , 2005 .

[4]  W. L. Wood Practical Time-Stepping Schemes , 1990 .

[5]  Rich R. Kerswell,et al.  Planar collapse of a granular column: Experiments and discrete element simulations , 2008 .

[6]  Jørgen Nielsen,et al.  Pressures from flowing granular solids in silos , 1998, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[7]  M. Gutierrez,et al.  Comprehensive study of the effects of rolling resistance on the stress–strain and strain localization behavior of granular materials , 2010 .

[8]  Andrei V. Lyamin,et al.  Computational Cam clay plasticity using second-order cone programming , 2012 .

[9]  Erling D. Andersen,et al.  On implementing a primal-dual interior-point method for conic quadratic optimization , 2003, Math. Program..

[10]  M. Anitescu,et al.  A matrix-free cone complementarity approach for solving large-scale, nonsmooth, rigid body dynamics , 2011 .

[11]  G. Bartels,et al.  The effect of contact torques on porosity of cohesive powders , 2004, cond-mat/0403110.

[12]  R. Kerswell,et al.  Dam break with Coulomb friction: a model for granular slumping , 2005 .

[13]  Sean McNamara,et al.  Measurement of indeterminacy in packings of perfectly rigid disks. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  Alfredo Edmundo Huespe,et al.  A nonlinear optimization procedure for limit analysis , 1996 .

[15]  R. Luciano,et al.  Stress-penalty method for unilateral contact problems: mathematical formulation and computational aspects , 1994 .

[16]  Pierre-Etienne Gautier,et al.  Modelling ballast behaviour under dynamic loading. Part 1: A 2D polygonal discrete element method approach , 2006 .

[17]  Jeffrey C. Trinkle,et al.  Complementarity formulations and existence of solutions of dynamic multi-rigid-body contact problems with coulomb friction , 1996, Math. Program..

[18]  D. Stewart,et al.  AN IMPLICIT TIME-STEPPING SCHEME FOR RIGID BODY DYNAMICS WITH INELASTIC COLLISIONS AND COULOMB FRICTION , 1996 .

[19]  Pierre-Yves Lagrée,et al.  The granular column collapse as a continuum: validity of a two-dimensional Navier–Stokes model with a μ(I)-rheology , 2011, Journal of Fluid Mechanics.

[20]  M. Anitescu,et al.  A Convex Complementarity Approach for Simulating Large Granular Flows , 2010 .

[21]  Michel Saint Jean,et al.  The non-smooth contact dynamics method , 1999 .

[22]  Armando N. Antão,et al.  Upper bound limit analysis with a parallel mixed finite element formulation , 2008 .

[23]  Jeffrey C. Trinkle,et al.  An implicit time-stepping scheme for rigid body dynamics with Coulomb friction , 2000, Proceedings 2000 ICRA. Millennium Conference. IEEE International Conference on Robotics and Automation. Symposia Proceedings (Cat. No.00CH37065).

[24]  Scott W. Sloan,et al.  Three-dimensional Mohr-Coulomb limit analysis using semidefinite programming , 2007 .

[25]  Roberto Zenit,et al.  Computer simulations of the collapse of a granular column , 2005 .

[26]  M. Anitescu,et al.  A Time-stepping Method for Stii Multibody Dynamics with Contact and Friction ‡ , 2022 .

[27]  Teruo Nakai,et al.  STRESS-DEFORMATION AND STRENGTH CHARACTERISTICS OF SOIL UNDER THREE DIFFERENT PRINCIPAL STRESSES , 1974 .

[28]  Armando N. Antão,et al.  A non‐linear programming method approach for upper bound limit analysis , 2007 .

[29]  János Török,et al.  Minimal dissipation theory and shear bands in biaxial tests , 2011 .

[30]  Andrei V. Lyamin,et al.  Granular contact dynamics using mathematical programming methods , 2012 .

[31]  E. J. Hinch,et al.  The spreading of a granular mass: role of grain properties and initial conditions , 2006 .

[32]  Roux,et al.  Force Distributions in Dense Two-Dimensional Granular Systems. , 1996, Physical review letters.

[33]  Gert Lube,et al.  Collapses of two-dimensional granular columns. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[34]  S. Sloan,et al.  Formulation and solution of some plasticity problems as conic programs , 2007 .

[35]  C. Thornton NUMERICAL SIMULATIONS OF DEVIATORIC SHEAR DEFORMATION OF GRANULAR MEDIA , 2000 .

[36]  José Herskovits,et al.  An iterative algorithm for limit analysis with nonlinear yield functions , 1993 .

[37]  F. Stillinger,et al.  Improving the Density of Jammed Disordered Packings Using Ellipsoids , 2004, Science.

[38]  Eric Vincens,et al.  Influence of particle shape and angularity on the behaviour of granular materials: a numerical analysis , 2003 .

[39]  Dietrich E. Wolf,et al.  Influence of particle elasticity in shear testers , 2006 .

[40]  Hans J. Herrmann,et al.  A micromechanical model of collapsing quicksand , 2010, 1012.2696.

[41]  Scott W. Sloan,et al.  Lower bound limit analysis with adaptive remeshing , 2005 .

[42]  J. M. Duncan,et al.  Elastoplastic Stress-Strain Theory for Cohesionless Soil , 1975 .

[43]  R. Kerswell,et al.  Axisymmetric granular collapse: a transient 3D flow test of viscoplasticity. , 2008, Physical review letters.

[44]  Farhang Radjai,et al.  Rheology, force transmission, and shear instabilities in frictional granular media from biaxial numerical tests using the contact dynamics method , 2005 .

[45]  Armando N. Antão,et al.  A novel augmented Lagrangian‐based formulation for upper‐bound limit analysis , 2012 .

[46]  H. Huppert,et al.  Axisymmetric collapses of granular columns , 2004, Journal of Fluid Mechanics.

[47]  Vincent Richefeu,et al.  Contact dynamics as a nonsmooth discrete element method , 2009 .

[48]  Jos F. Sturm,et al.  A Matlab toolbox for optimization over symmetric cones , 1999 .

[49]  F. Radjai,et al.  Identification of rolling resistance as a shape parameter in sheared granular media. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[50]  Andrei V. Lyamin,et al.  Granular contact dynamics with particle elasticity , 2012, Granular Matter.

[51]  E. J. Hinch,et al.  Study of the collapse of granular columns using two-dimensional discrete-grain simulation , 2005, Journal of Fluid Mechanics.

[52]  J. Moreau,et al.  Unilateral Contact and Dry Friction in Finite Freedom Dynamics , 1988 .

[53]  P. Cundall,et al.  A discrete numerical model for granular assemblies , 1979 .

[54]  Arkadi Nemirovski,et al.  Lectures on modern convex optimization - analysis, algorithms, and engineering applications , 2001, MPS-SIAM series on optimization.

[55]  Alexander Ries,et al.  Shear zones in granular media: three-dimensional contact dynamics simulation. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[56]  Armando N. Antão,et al.  An upper bound-based solution for the shape factors of bearing capacity of footings under drained conditions using a parallelized mixed f.e. formulation with quadratic velocity fields , 2012 .

[57]  Kristian Krabbenhoft,et al.  A general non‐linear optimization algorithm for lower bound limit analysis , 2003 .

[58]  Pierre Alart,et al.  A parallel version of the non smooth contact dynamics algorithm applied to the simulation of granular media , 2004 .

[59]  Jean-Pierre Vilotte,et al.  On the use of Saint Venant equations to simulate the spreading of a granular mass , 2005 .

[60]  Antonio Huerta,et al.  Upper and lower bounds in limit analysis: Adaptive meshing strategies and discontinuous loading , 2009 .

[61]  I Zuriguel,et al.  Effect of particle shape on the stress dip under a sandpile. , 2007, Physical review letters.