Search results

Granular materials possess multiscale structures, i.e. micro-scales involving atoms and molecules in a solid particle, meso-scales involving individual particles and their correlated structure, and macroscopic assembly. Strong and abundant dissipations are exhibited due to mesoscopic unsteady motion of individual grains, and evolution of underlying structures (e.g. force chains, vortex, etc.), which defines the key differences between granular materials and ordinary objects. The purpose of this paper is to introduce the major studies have been conducted in recent two decades.

The main properties at individual scale are introduced, including the coordination number, pair-correlation function, force and mean stress distribution functions, and the dynamic correlation function. The relationship between meso- and macro-scales is analyzed, such as between contact force and stress, the elastic modulus, and bulk friction in granular flows. At macroscales, conventional engineering models (i.e. elasto-plastic and hypo-plastic ones) are introduced. In particular, the so-called granular hydrodynamics theory, derived from thermodynamics principles, is explained.

On the basis of recent study the authors conducted, the multiscales (both spatial and temporal) in granular materials are first explained, and a multiscale framework is presented for the mechanics of granular materials.

It would provide a paramount view on the multiscale studies of granular materials.

The purpose of this paper is to develop a discrete element (DE) and multiscale modeling methodology to represent granular media at their particle scale as they interface solid deformable bodies, such as soil‐tool, tire, penetrometer, pile, etc., interfaces.

A three‐dimensional ellipsoidal discrete element method (DEM) is developed to more physically represent particle shape in granular media while retaining the efficiency of smooth contact interface conditions for computation. DE coupling to finite element (FE) facets is presented to demonstrate initially the development of overlapping bridging scale methods for concurrent multiscale modeling of granular media.

A closed‐form solution of ellipsoidal particle contact resolution and stiffness is presented and demonstrated for two particle, and many particle contact simulations, during gravity deposition, and quasi‐static oedometer, triaxial compression, and pile penetration. The DE‐FE facet coupling demonstrates the potential to alleviate artificial boundary effects in the shear deformation region between DEM granular media and deformable solid bodies.

The research is being extended to couple more robustly the ellipsoidal DEM code and a higher order continuum FE code via overlapping bridging scale methods, in order to remove dependence of penetration/shear resistance on the boundary placement for DE simulation.

When concurrent multiscale computational modeling of interface conditions between deformable solid bodies and granular materials reaches maturity, modelers will be able to simulate the mechanical behavior accounting for physical particle sizes and flow in the interface region, and thus design their tool, tire, penetrometer, or pile accordingly.

A closed‐form solution for ellipsoidal particle contact is demonstrated in this paper, and the ability to couple DE to FE facets.

A microstructural finite element model (MFEM) for granular material considering the microstructure of material is presentented. In the MFEM method, a volume of large number of particles is represented by an element consisting of a few nodal points. The stiffness matrix of the element is then formulated based on the contact stiffness and the packing arrangement of the particles. The method can be easily applied in the framework of finite element technique to solve boundary value problems in practical situations. Applicability of the model is evaluated by comparing the results of MFEM model with that from DEM model.

The distinct element method has been used to perform a number of numerical experiments, analogous to physical tests routinely carried out for real geologic materials. The ability to prescribe exactly the material properties and the boundary conditions allows us to study in detail the influence of microscopic structure on macroscopic response. These macroscopic quantities are usually the only data obtained in physical experiments. We are now in a position to formulate a general constitutive model, based on physical mechanisms and appropriate for use in a continuum code.

The purpose of this paper is to extend the bridge scale method (BSM) developed for granular materials with only the solid phase to that taking into account the effects of wetting process in porous continuum. The granular material is modeled as partially saturated porous Cosserat continuum and discrete particle assembly in the coarse and fine scales, respectively.

Based on the mass and momentum conservation laws for the three phases, i.e. the solid skeleton, the pore water and the pore air, the governing equations for the unsaturated porous Biot-Cosserat continuum model in the coarse scale are derived. In light of the passive air pressure assumption, a reduced finite element model for the model is proposed. According to the decoupling of the fine and coarse scale calculations in the BSM, the unsaturated porous Cosserat continuum model using the finite element method and the discrete element model using the discrete element method for granular media are combined.

The numerical results for a 2D example problem of slope stability subjected to increasing rainfall along with mechanical loading demonstrate the applicability and performance of the present BSM. The microscopic mechanisms of macroscopic shear band developed in the slope are demonstrated.

Do not account for yet the effects of unsaturated pore water in the fine scale.

The novel BSM that couples the Biot-Cosserat porous continuum modeling and the discrete particle assembly modeling in both coarse and fine scales, respectively, is proposed to provide a micro-macro discrete-continuum two-scale modeling approach for numerical simulations of the hydro-mechanical coupling problems in unsaturated granular materials.

The purpose of this paper is to examine several calibration techniques that have been developed to determine the discrete element method (DEM) parameters for slow and rapid unconfined flow of granular conical pile formation. This paper also aims to discuss some of the methods currently employed to scale particle properties to reduce computational resources and time to solve large DEM models.

DEM models have been calibrated against simple bench‐scale experimental results to examine the validity of selected parameters for the contact, material and mechanical models to simulate the dynamic and static behaviour of cohesionless polyethylene pellets. Methods to determine quantifiable single particle parameters such as static friction and the coefficient of restitution have been highlighted. Numerical and experimental granular pile formation has been investigated using different slumping and pouring techniques to examine the dependency of the type of flow mechanism on the DEM parameters.

The proposed methods can provide cost effective and simple techniques to determine suitable input parameters for DEM models. Rolling friction and particle shape representation has shown to have a significant influence on the bulk flow characteristics via a sensitivity analysis and needs to be accessed based on the environmental conditions.

This paper describes several effective known and novel methodologies to characterise granular materials that are needed to accurately model granular flow using the DEM to provide valuable quantitative data. For the DEM to be a viable predictive tool in industrial applications which often contain huge quantities of particles with random particle shapes and irregular properties, quick and validated techniques to “tune” DEM models are necessary.

We have developed a general constitutive theory that estimates the effective elastic moduli of a cemented granular material by applying statistical mechanical averaging to a purely micromechanical model. We have also constructed a distinct element model of a cemented granular material, based on the same micromechanical model, which accounts for the elastic forces due to bonding between pairs of spherical particles, and which allows for the possibility of anisotropic damage to the bonds. In this paper, we use a model based on the distinct element method (DEM) to validate the predictions of the theory for various prescribed patterns of damage. In particular we impose several anisotropic patterns of damage on the bonds of a randomly generated assembly of particles. We then undertake numerical experiments, sending both p‐waves and s‐waves through the samples and measuring the wave velocities. The predictions of the theory for these velocities agree well with the results of the numerical model for a variety of damage patterns. We discuss the implications of our theory, as well as potential applications.

A consistent implementation of the general computational framework of unified second-order time accurate integrators via the well-known GSSSS framework in conjunction with the traditional Finite Difference Method is presented to improve the numerical simulations of reactive two-phase flows.

In the present paper, the phase interaction evaluation in the present implementation of the reactive two-phase flows has been derived and implemented to preserve the consistency of the correct time level evaluation during the time integration process for solving the two phase flow dynamics with reactions.

Numerical examples, including the classical Sod shock tube problem and a reactive two-phase flow problem, are exploited to validate the proposed time integration framework and families of algorithms consistently to second order in time accuracy; this is in contrast to the traditional practices which only seem to obtain first-order time accuracy because of the inconsistent time level implementation with respect to the interaction of two phases. The comparisons with the traditional implementation and the advantages of the proposed implementation are given in terms of the improved numerical accuracy in time. The proposed approaches provide a correct numerical simulation implementation to the reactive two-phase flows and can obtain better numerical stability and computational features.

The new algorithmic framework and the consistent time level evaluation extended with the GS4 family encompasses a multitude of past and new schemes and offers a general purpose and unified implementation for fluid dynamics.

The study aims to investigate the effects of pre-loading histories (pre-shearing and pre-consolidation) on the liquefaction behaviour of saturated loose sand via discrete element method (DEM) simulations.

The pre-shearing history is mimicked under drained conditions (triaxial compression) with different pre-shearing strain levels ranging from 0% to 2%. The pre-consolidation history is mimicked by increasing the isotropic compression to different levels ranging from 100 kPa to 300 kPa. The macroscopic and microscopic behaviours are analysed and compared.

Temporary liquefaction, or quasi-steady state (QSS), is observed in most samples. A higher pre-shearing or pre-consolidation level can provide higher liquefaction resistance. The ultimate state line is found to be unique and independent of the pre-loading histories in stress space. The Lade instability line prematurely predicts the onset of liquefaction for all samples, both with and without pre-loading histories. The redundancy index is an effective microscopic indicator to monitor liquefaction, and the onset of the liquefaction corresponds to the phase transition state where the value of redundancy index is one, which is true for all cases irrespective of the proportions of sliding contacts.

The liquefaction behaviour of granular materials still remains elusive, especially concerning the effects of pre-loading histories on soils. Furthermore, the investigation of the effects of pre-consolidation histories on undrained behaviour and its comparison to pre-sheared samples is rarely reported in the DEM literature.

This paper aims to determine how the viscosity and curing agent content affect the flowability of moist silica sand granules. In addition, a coating device was designed according to the flow properties of silica sand granules.

The flowability of silica sand granules premixed with two curing agents of different viscosities is studied using a Jenike shear apparatus. An open-ended device was used in discharge testing of sand granules with a design based on the variable dip angle of the two plates and variable outlet size.

The test results show that increasing the curing agent content would significantly decrease the flowability of silica sand granules, and a curing agent of higher viscosity has a greater effect on the flowability of silica sand. The presence of a curing agent strengthens the cohesion among sand granules, lubricates them and restrains their deformation. The shape function of the coating device was obtained by theoretical derivation.

The flow properties provide a valuable theoretical guidance for the design of coating device for sand mold printing.

This paper deals with experimental work on flow properties of silica sand granules with different viscosities and curing agent content. The shape function of a wedge-shaped coating device is obtained based on experimental data.