Simulation of railway ballast using crushable polyhedral particles

Abstract The need to understand and predict the long- and short term behavior of railway ballast calls for extensive experimental programs as well as the development of reliable numerical models. From the modeling point of view, the granular nature of ballast in connection with the high onset speed of dynamic loading makes it ideal for the application of the Discrete Element Method (DEM). The paper describes the employment of the DEM to simulate ballast behavior in large-scale oedometric testing. Ballast grains are represented by convex polyhedral particles whose shapes are randomly generated via Voronoi tessellation. The effect of the aspect ratio of particles is investigated. A novel algorithm to compute repulsive contact force based on the intersecting volume of polyhedrons is presented. Crushing of grains is included via splitting the particles into smaller polyhedrons when a certain stress-based criterion is fulfilled. The results obtained from the model are compared with the published experimental investigations.

[1]  Dawei Zhao,et al.  A fast contact detection algorithm for 3-D discrete element method , 2004 .

[2]  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.

[3]  Harald Kruggel-Emden,et al.  Comparison of the multi-sphere and polyhedral approach to simulate non-spherical particles within the discrete element method: Influence on temporal force evolution for multiple contacts , 2011 .

[4]  Magnus Evertsson,et al.  The contribution of DEM to the science of comminution , 2013 .

[5]  Fulvio Tonon,et al.  Explicit Exact Formulas for the 3-D Tetrahedron Inertia Tensor in Terms of its Vertex Coordinates , 2005 .

[6]  Katalin Bagi,et al.  Stress and strain in granular assemblies , 1996 .

[7]  Glenn R. McDowell,et al.  A method to model realistic particle shape and inertia in DEM , 2010 .

[8]  Siegmar Wirtz,et al.  A study on the influence of particle shape and shape approximation on particle mechanics in a rotating drum using the discrete element method , 2014 .

[9]  Mohammad Hossein Abbaspour-Fard,et al.  Modeling nonspherical particles using multisphere discrete elements , 2001 .

[10]  Luis E. Vallejo,et al.  Crushing a weak granular material: experimental numerical analyses , 2005 .

[11]  Qicheng Sun,et al.  Probability-based contact algorithm for non-spherical particles in DEM , 2011 .

[12]  Jun Yang,et al.  A discrete element analysis of elastic properties of granular materials , 2013, Granular Matter.

[13]  E. Ramm,et al.  On the application of a discrete model to the fracture process of cohesive granular materials , 2002 .

[14]  J. Roux,et al.  3D numerical simulation study of quasistatic grinding process on a model granular material , 2013 .

[15]  Soheil Mohammadi,et al.  Micromechanics of breakage in sharp-edge particles using combined DEM and FEM , 2008 .

[16]  I. N. Sneddon.  The relation between load and penetration in the axisymmetric boussinesq problem for a punch of arbitrary profile , 1965 .

[17]  W. L. Lim,et al.  Mechanics of railway ballast behaviour , 2004 .

[18]  J. Bullard,et al.  Defining shape measures for 3D star-shaped particles: Sphericity, roundness, and dimensions , 2013 .

[19]  J. Goddard.  Nonlinear elasticity and pressure-dependent wave speeds in granular media , 1990, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences.

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

[21]  H. Hertz.  Ueber die Berührung fester elastischer Körper. , 1882 .

[22]  Buddhima Indraratna,et al.  DEM analysis of angular ballast breakage under cyclic loading , 2007 .

[23]  Erol Tutumluer,et al.  Discrete element modelling of ballasted track deformation behaviour , 2013 .

[24]  Richard Turton,et al.  Contact detection algorithms for DEM simulations of tablet-shaped particles , 2006 .

[25]  Richard J. Bathurst,et al.  Micromechanical features of granular assemblies with planar elliptical particles , 1992 .

[26]  G. McDowell,et al.  Discrete element modelling of railway ballast , 2005 .

[27]  Pierre Alliez,et al.  Computational geometry algorithms library , 2008, SIGGRAPH '08.

[28]  L. Vallejo,et al.  Discrete Element Method Analysis of Railtrack Ballast Degradation during Cyclic Loading , 2006 .

[29]  Lutz Kettner,et al.  Using generic programming for designing a data structure for polyhedral surfaces , 1999, Comput. Geom..

[30]  Jennifer S. Curtis,et al.  Discrete element method modeling of bi-convex pharmaceutical tablets: Contact detection algorithms and validation , 2012 .

[31]  Anthony Wachs,et al.  Grains3D, a flexible DEM approach for particles of arbitrary convex shape — Part I: Numerical model and validations , 2012 .

[32]  D. Asahina,et al.  Geometry of irregular particles: Direct surface measurements by 3-D laser scanner , 2011 .

[33]  Feras Y. Fraige,et al.  Vibration induced flow in hoppers： DEM 2D polygon model , 2008 .

[34]  David E. Muller,et al.  Finding the Intersection of two Convex Polyhedra , 1978, Theor. Comput. Sci..

[35]  John E. Bolander,et al.  Voronoi-based discretizations for fracture analysis of particulate materials , 2011 .

[36]  Fabrice Emeriault,et al.  Statistical homogenization for assemblies of elliptical grains: effect of the aspect ratio and particle orientation , 2004 .

[37]  Algis Dziugys,et al.  A new approach to detect the contact of two‐dimensional elliptical particles , 2001 .

[38]  J. Kozicki,et al.  A new open-source software developed for numerical simulations using discrete modeling methods , 2008 .

[39]  P. Eberhard,et al.  A discrete element model predicting the strength of ballast stones , 2012 .

[40]  Guy T. Houlsby,et al.  A new algorithm for contact detection between convex polygonal and polyhedral particles in the discrete element method , 2012 .

[41]  Hai-Sui Yu,et al.  Particle shape effects in discrete element modelling of cohesive angular particles , 2011 .

[42]  P. Cundall,et al.  FORMULATION OF A THREE-DIMENSIONAL DISTINCT ELEMENT MODEL - PART II. MECHANICAL CALCULATIONS FOR MOTION AND INTERACTION OF A SYSTEM COMPOSED OF MANY POLYHEDRAL BLOCKS , 1988 .

[43]  P. A. Cundall,et al.  FORMULATION OF A THREE-DIMENSIONAL DISTINCT ELEMENT MODEL - PART I. A SCHEME TO DETECT AND REPRESENT CONTACTS IN A SYSTEM COMPOSED OF MANY POLYHEDRAL BLOCKS , 1988 .

[44]  Youssef M A Hashash,et al.  Shortest link method for contact detection in discrete element method , 2006 .

[45]  Siegmar Wirtz,et al.  Experimental and numerical investigation on the influence of particle shape and shape approximation on hopper discharge using the discrete element method , 2013 .

[46]  Hai Huang,et al.  Discrete element modeling of railroad ballast using imaging based aggregate morphology characterization , 2010 .

[47]  G. Combe,et al.  Polygons vs. clumps of discs: a numerical study of the influence of grain shape on the mechanical behaviour of granular materials , 2011, 1208.0565.