Discrete element method for modelling solid and particulate materials

The discrete element method (DEM) is developed in this study as a general and robust technique for unified two-dimensional modelling of the mechanical behaviour of solid and particulate materials, including the transition from solid phase to particulate phase. Inter-element parameters (contact stiffnesses and failure criteria) are theoretically established as functions of element size and commonly accepted material parameters including Young’s modulus, Poisson’s ratio, ultimate tensile strength, and fracture toughness. A main feature of such an approach is that it promises to provide convergence with refinement of a DEM discretization. Regarding contact failure, an energy criterion based on the material’s ultimate tensile strength and fracture toughness is developed to limit the maximum contact forces and inter-element relative displacement. This paper also addresses the issue of numerical stability in DEM computations and provides a theoretical method for the determination of a stable time-step. The method developed herein is validated by modelling several test problems having analytic solutions and results show that indeed convergence is obtained. Moreover, a very good agreement with the theoretical results is obtained in both elastic behaviour and fracture. An example application of the method to high-speed penetration of a concrete beam is also given. Copyright q 2006 John Wiley & Sons, Ltd