A nonlocal multiscale discrete-continuum model for predicting mechanical behavior of granular materials

Abstract

Summary A three‐dimensional nonlocal multiscale discrete‐continuum model has been developed for modeling mechanical behavior of granular materials. In the proposed multiscale scheme, we establish an information‐passing coupling between the discrete element method, which explicitly replicates granular motion of individual particles, and a finite element continuum model, which captures nonlocal overall responses of the granular assemblies. The resulting multiscale discrete‐continuum coupling method retains the simplicity and efficiency of a continuum‐based finite element model, while circumventing mesh pathology in the post‐bifurcation regime by means of staggered nonlocal operator. We demonstrate that the multiscale coupling scheme is able to capture the plastic dilatancy and pressure‐sensitive frictional responses commonly observed inside dilatant shear bands, without employing a phenomenological plasticity model at a macroscopic level. In addition, internal variables, such as plastic dilatancy and plastic flow direction, are now inferred directly from granular physics, without introducing unnecessary empirical relations and phenomenology. The simple shear and the biaxial compression tests are used to analyze the onset and evolution of shear bands in granular materials and sensitivity to mesh density. The robustness and the accuracy of the proposed multiscale model are verified in comparisons with single‐scale benchmark discrete element method simulations. Copyright © 2015 John Wiley & Sons, Ltd.

Multiscale modeling of large deformation in geomechanics

“…Single element test has been widely used for numerical model validations . It is used here to benchmark the multiscale modeling approach.…”

A Semi-Implicit Microplar Discrete-to-Continuum Method for Granular Materials

“…Previous work on hierarchical DEM-FEM coupling models have been mainly focused on establishing coupling between discrete element simulations and Cauchy-type-continuum finite element analysis via first-order homogenization procedure [17,6,18,19]. This first order homogenization is, nevertheless, only valid if the separation of scales exists [20,21].…”

“…When the particles of the granular assembly are neither bonded nor subjected to confining pressure, granular flow may occur [3]. On the other hand, the collective macroscopic responses of these interacting particles may exhibit traits that are similar to those of a solid continuum [4,5,6,1,7].…”

Restoring Mesh Independency in FEM-DEM Multi-scale Modelling of Strain Localization Using Second Gradient Regularization

“…Indeed, FEMxDEM methods allow to couple the advantages of Discrete Elements and the efficiency of Finite Elements. Later works have enhanced and extended this approach to the study of anisotropy [7,18], granular cohesion [19], material heterogeneity [25], real scale engineering applications [18,8], more realistic constitutive behaviors using 3D DEM [12,26], macroscale hydro-mechanical coupling [26,10]. More recently, [12] have embedded non-local regularization at the macro-scale and [9] has developed a full micromacro 3D approach.…”

“…Regularization techniques have been developed to overcome these problem; both nonlocal [5] and local [17,6] approaches exist. An example of nonlocal regularization in a FEMxDEM model is presented in [12]. Local formulations use a local relationship between stress and strains in the same manner as classical constitutive relations are defined.Second gradient model, as a particular case of the Germain theory [6] has been developed [4,21,1,14].…”