Adaptive model reduction for nonsmooth discrete element simulation

A method for adaptive model order reduction for nonsmooth discrete element simulation is developed and analysed in numerical experiments. Regions of the granular media that collectively move as rigid bodies are substituted with rigid bodies of the corresponding shape and mass distribution. The method also support particles merging with articulated multibody systems. A model approximation error is defined and used to derive conditions for when and where to apply reduction and refinement back into particles and smaller rigid bodies. Three methods for refinement are proposed and tested: prediction from contact events, trial solutions computed in the background and using split sensors. The computational performance can be increased by 5–50 times for model reduction level between 70–95 %.

[1]  Christian Miehe,et al.  Homogenization and two‐scale simulations of granular materials for different microstructural constraints , 2010 .

[2]  T. Schwager,et al.  Computational Granular Dynamics: Models and Algorithms , 2005 .

[3]  Ming C. Lin,et al.  An efficient, error-bounded approximation algorithm for simulating quasi-statics of complex linkages , 2005, SPM '05.

[4]  C. Schütte,et al.  Homogenization of Hamiltonian systems with a strong constraining potential , 1997 .

[5]  Athanasios C. Antoulas,et al.  Approximation of Large-Scale Dynamical Systems , 2005, Advances in Design and Control.

[6]  P. Glösmann.  Reduction of discrete element models by Karhunen–Loève transform: a hybrid model approach , 2010 .

[7]  G. Kerschen,et al.  The Method of Proper Orthogonal Decomposition for Dynamical Characterization and Order Reduction of Mechanical Systems: An Overview , 2005 .

[8]  Peter Eberhard,et al.  Model Order Reduction in Elastic Multibody Systems using the Floating Frame of Reference Formulation , 2012 .

[9]  Ming C. Lin,et al.  Adaptive dynamics of articulated bodies , 2005, SIGGRAPH 2005.

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

[11]  C. Lacoursière,et al.  Examining the smooth and nonsmooth discrete element approaches to granular matter , 2014 .

[12]  A. Munjiza.  The Combined Finite-Discrete Element Method , 2004 .

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

[14]  C. Martin,et al.  A clumped particle model for rock , 2007 .

[15]  J. Moreau.  Numerical aspects of the sweeping process , 1999 .

[16]  P. Wriggers,et al.  A two-scale model of granular materials , 2012 .

[17]  H. Jaeger,et al.  Granular solids, liquids, and gases , 1996 .

[18]  Thorsten Pöschel,et al.  Computational Granular Dynamics , 2005 .

[19]  C. Lacoursière.  Regularized; Stabilized; Variational Methods for Multibodies , 2007 .

[20]  P. Alart,et al.  Numerical simulation of 2D steady granular ows in rotating drum: On surface ows rheology , 2005 .

[21]  Katta G. Murty,et al.  Linear complementarity, linear and nonlinear programming , 1988 .

[22]  Ning Guo,et al.  A coupled FEM/DEM approach for hierarchical multiscale modelling of granular media , 2014 .

[23]  A. Munjiza.  The Combined Finite-Discrete Element Method: Munjiza/Discrete Element Method , 2004 .

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