Institute of Fundamental Technological Research

This work presents a novel scheme to couple the Boundary Element Method (BEM) and the Discrete Element Method (DEM) in the time domain. The DEM captures discontinuous material behaviour, such as fractured and granular media. However, applying the method to real-life applications embedded into infinite domains is challenging. The authors propose a solution to this challenge by coupling the DEM with the BEM. The capability of the BEM to model infinite domains accurately and efficiently, without the need for numerical artifices, makes it the perfect complement to the DEM. This study proposes a direct monolithic interface-based coupling method that resolves any incompatibilities between the two methods in two dimensions. The benchmark results show that the proposed methodology consistently produces results that align with analytical solutions. The final example in the paper showcases the full potential of this innovative methodology, where the DEM models a fracturing process, and the BEM evaluates its far-field effect.

Discrete Element Method (DEM), Boundary Element Method (BEM), Discontinuous materials, Wave propagation, Infinite domain, Monolithic coupling

The discrete element method (DEM) is the most dominant method for the numerical prediction of dynamic behaviour at grain or particle scale. Nevertheless, due to its discontinuous nature, the DEM is inherently unable to describe microscopic features of individual bodies which can be considered as continuous bodies. To incorporate microscopic features, efficient numerical coupling of the DEM with a continuous method is generally necessary. Thus, a generalised multi-scale PD–DEM framework is developed in this work. In the developed framework, meshfree discretised Peridynamics (PD) is used to describe intra-particle forces within bodies to capture microscopic features. The inter-particle forces of rigid bodies are defined by the DEM whereas a hybrid approach is applied at the PD–DEM interface. In addition, a staggered multi-scale time integration scheme is formulated to allow for an efficient numerical treatment of both methods. Validation examples are presented and the applicability of the developed framework to capture the characteristics mixtures with rigid and deformable bodies is shown.

Peridynamics (PD),Discrete element method (DEM),Contact coupling,Multi-scale modelling,Deformable particles