Numerical and experimental direct shear tests for coarse-grained soils

The presence of particles larger than the permissible dimensions of conventional laboratory specimens causes difficulty in the determination of shear strength of coarse-grained soils. In this research, the influence of particle size on shear strength of coarse-grained soils was investigated by resorting to experimental tests in different scale and numerical simulations based on discrete element method (DEM). Experimental tests on such soil specimens were based on using the techniques designated as "parallel" and "scalping" to prepare gradation of samples in view of the limitation of laboratory specimen size. As a second approach, the direct shear test was numerically simulated on assemblies of elliptical particles. The behaviors of samples under experimental and numerical tests are presented and compared, indicating that the modification of sample gradation has a significant influence on the mechanical properties of coarse-grained soils. It is noted that the shear strengths of samples produced by the scalping method are higher than samples by the parallel method. The scalping method for preparing specimens for direct shear test is therefore recommended. The micromechanical behavior of assemblies under direct shear test is also discussed and the effects of stress level on sample behavior are investigated.     

Discrete element modeling of inherently anisotropic granular assemblies with polygonal particles

In the present article, we study the effect of inherent anisotropy, i.e., initial bedding angle of particles and associated voids on macroscopic mechanical behavior of granular materials, by numerical simulation of several biaxial compression tests using the discrete element method (DEM). Particle shape is considered to be irregular convex-polygonal. The effect of inherent anisotropy is investigated by following the evolution of mobilized shear strength and volume change during loading. As experimental tests have already shown, numerical simulations also indicate that initial anisotropic condition has a great influence on the strength and deformational behavior of granular assemblies. Comparison of simulations with tests using oval particles, shows that angularity influences both the mobilized shear strength and the volume change regime, which originates from the interlocking resistance between particles.     

DISCRETE AND CONTINUUM MODELLING OF GRANULAR FLOW

This paper analyses three popular methods simulating granular flow at different time and length scales： discrete element method （DEM）, averaging method and viscous, elastic-plastic continuum model. The theoretical models of these methods and their applications to hopper flows are discussed. It is shown that DEM is an effective method to study the fundamentals of granular flow at a particle or microscopic scale. By use of the continuum approach, granular flow can also be described at a continuum or macroscopic scale. Macroscopic quantities such as velocity and stress can be obtained by use of such computational method as FEM. However, this approach depends on the constitutive relationship of materials and ignores the effect of microscopic structure of granular flow. The combined approach of DEM and averaging method can overcome this problem. The approach takes into account the discrete nature of granular materials and does not require any global assumption and thus allows a better understanding of the fundamental mechanisms of granular flow. However, it is difficult to adapt this approach to process modelling because of the limited number of particles which can be handled with the present computational capacity, and the difficulty in handling non-spherical particles. Further work is needed to develoo an aoorooriate aooroach to overcome these problems.     

ROTATIONAL RESISTANCE AND SHEAR-INDUCED ANISOTROPY IN GRANULAR MEDIA

This paper presents a micromechanical study on the behavior of granular materials under confined shear using a three-dimensional Discrete Element Method （DEM）. We consider rotational resistance among spherical particles in the DEM code as an approximate way to account for the effect of particle shape. Under undrained shear, it is found rotational resistance may help to increase the shear strength of a granular system and to enhance its resistance to liquefaction. The evolution of internal structure and anisotropy in granular systems with different initial conditions depict a clear bimodal character which distinguishes two contact subnetworks. In the presence of rotational resistance, a good correlation is found between an analytical stress-force-fabric relation and the DEM results, in which the normal force anisotropy plays a dominant role. The unique properties of critical state and liquefaction state in relation to granular anisotropy are also explored and discussed.     

A numerical analysis of the shear behavior of granular soil with fines

Shear behavior of granular soil with fines is investigated using the discrete element method(DEM) and particle arrangements and inter-particle contacts during shear are examined.The DEM simulation reveals that fine particles play a vital role in the overall response of granular soil to shearing.The occurrence of liquefaction and temporary reduction of strength is ascribed mainly to the loss of support from the fine particle contacts(S-S) and fine particle-to-large particle contacts(S-L) as a consequence of the removal of fine particles from the load-carrying skeleton.The dilative strain-hardening response following the strain-softening response is associated with the migration of fine particles back into the load-carrying skeleton,which is thought to enhance the stiffness of the soil skeleton.During shear,the unit normal vector of the large particle-to-large particle(L-L) contact has the strongest fabric anisotropy,and the S-S contact unit normal vector possesses the weakest anisotropy,suggesting that the large particles play a dominant role in carrying the shear load.It is also found that,during shear,fine particles are prone to rolling at contacts while the large particles are prone to sliding,mainly at the S-L and L-L contacts.     

Segregation of particulate solids: Experiments and DEM simulations

Segregation of granular materials is a complex phenomenon, difficult to measure quantitatively and to predict. Discrete element method (DEM) can be a useful tool to predict segregation effects and to support the industrial design. In this context, a very challenging idea is the characterization of the granular solids to provide the key parameters needed for a successful DEM simulation of segregation processes. Rolling friction, sliding friction and the coefficient of restitution are the critical parameters to be studied. These microscopic simulation parameters are calibrated by comparing the macroscopic behavior of granular matter in standard bulk experiments, which have the advantage of being highly repeatable and reliable.An experimental method is presented to characterize free surface segregation. The effects of different particle properties, particularly, shape and size, on segregation of cohesionless materials were investigated. From the experiments, particle size demonstrated a stronger effect on segregation than particle shape. Finally, the corresponding DEM simulations of the segregation experiments were presented. The parameters obtained by calibration were validated by the comparison of the modeled segregation behavior with the experimental results. Thus, calibrated DEM simulations are capable of predicting segregation effects.     

Segregation of particulate solids＂ Experiments and DEM simulations

Segregation of granular materials is a complex phenomenon, difficult to measure quantitatively and to predict. Discrete element method （DEM） can be a useful tool to predict segregation effects and to support the industrial design. In this context, a very challenging idea is the characterization of the granular solids to provide the key parameters needed for a successful DEM simulation of segregation processes. Rolling friction, sliding friction and the coefficient of restitution are the critical parameters to be studied. These microscopic simulation parameters are calibrated by comparing the macroscopic behavior of granular matter in standard bulk experiments, which have the advantage of being highly repeatable and reliable. An experimental method is presented to characterize free surface segregation. The effects of different particle properties, particularly, shape and size, on segregation of cohesionless materials were investi- gated. From the experiments, particle size demonstrated a stronger effect on segregation than particle shape. Finally, the corresponding DEM simulations of the segregation experiments were presented. The parameters obtained by calibration were validated by the comparison of the modeled segregation behav- ior with the experimental results. Thus, calibrated DEM simulations are capable of predicting segregation effects.     

Assessing mixing characteristics of particle-mixing and granulation devices

The mixing of particulates such as powders is an important process in many industries including pharmaceuticals, plastics, household products (such as detergents) and food processing. The quality of products depends on the degree of mixing of their constituent materials which in turn depends on both geometric design and operating conditions. Unfortunately, due to lack of understanding of the interaction between mixer geometry and the granular material, limited progress has been made in optimizing mixer design. The discrete element method (DEM) is a computational technique that allows particle systems to be simulated and mixing to be predicted. Simulation is an effective way of acquiring information on the performance of different mixers that is difficult and/or expensive to obtain using traditional experimental approaches. Here we demonstrate how DEM can be used to unravel flow dynamics and assess mixing in several different types of devices. These devices used for mixing and/or granulation of particulates, are classified broadly as gravity controlled, bladed and high shear. We also explore the role of particle shape in mixing performance and use DEM to test whether Froude number scaling is suitable for predicting scale performance of rotating mixers.     

Calibration of granular material parameters for DEM modelling and numerical verification by blade-granular material interaction

The discrete element method (DEM) is a promising approach to model blade-granular material interactions. The accuracy of DEM models depends on the model parameters. In this study, a calibration process was developed to determine the parameter values. The particle size was the same as the real material and the particle shape was modelled using two spherical particles rigidly clumped together to form a single grain. Laboratory shear tests and compressions tests were used to determine the material internal friction angle and stiffness, respectively. These tests were replicated numerically using DEM models with different sets of particle friction coefficients and particle stiffness values. The shear test results are found to be dependent on both the particle friction coefficient and the particle stiffness. The compression test results show that it is only dependent on the particle stiffness. The combination of shear test and compression test results can be used to determine a unique set of particle friction and particle stiffness values. The calibration process was validated experimentally and numerically by modelling a blade moving through granular material. Results show that the forces acting on the blade can be accurately modelled with DEM and the maximum error is found to be 26%. The relative particle-blade displacements were used to predict the position and shape of the shear lines in front of the blade. A good qualitative correlation was achieved between the experiments and the DEM simulations.     

Assessing mixing characteristics of particle-mixing and granulation devices

The mixing of particulates such as powders is an important process in many industries including pharmaceuticals, plastics, household products （such as detergents） and food processing. The quality of products depends on the degree of mixing of their constituent materials which in turn depends on both geometric design and operating conditions. Unfortunately, due to lack of understanding of the interaction between mixer geometry and the granular material, limited progress has been made in optimizing mixer design. The discrete element method （DEM） is a computational technique that allows particle systems to be simulated and mixing to be predicted. Simulation is an effective way of acquiring information on the performance of different mixers that is difficult and/or expensive to obtain using traditional experimental approaches. Here we demonstrate how DEM can be used to unravel flow dynamics and assess mixing in several different types of devices. These devices used for mixing and/or granulation of particulates, are classified broadly as gravity controlled, bladed and high shear. We also explore the role of particle shape in mixing performance and use DEM to test whether Froude number scaling is suitable for predicting scale performance of rotating mixers.     

DEM investigation on conveying of non-spherical particles in a screw conveyor

Screw conveyors are extensively used in modern industry such as metallurgy, architecture and pharmaceutical due to their high-efficiency in the transportation of granular materials. And substantial efforts have been devoted to the study of the screw conveyors. Numerical method is an effective way to study screw conveyor. However, previous studies have mainly focused in the regime of spherical particles while the in-depth investigations for non-spherical particles that should be the most encountered in practical applications are still limited. In view of the above situations, discrete element method (DEM), which has been widely accepted in simulating the discrete systems, is utilized to investigate the conveying process of non-spherical particles in a horizontal screw conveyor, with particles being modeled by super-ellipsoids. In addition, a wear model called SIEM (Shear Impact Energy Model) is incorporated into DEM to predict the wear of screw conveyor. The DEM simulation results demonstrate that the particle shape is influential for the flow behaviors of particles and the wear of conveyor. The conveying performance evaluated quantitatively of both mass flow rate and power consumption is subsequently obtained to investigate the effect of sphericity of particle with different operation parameters. Moreover, particle collision frequency and collision energy consumption are acquired to investigate the possible particle breakage between particles and screw blade. The comparisons between particle–particle collision and particle–wall collision reveal that particles with large shape index have more possibility to be damaged in particle–wall impingement.     

Discrete simulation and micromechanical analysis of two-dimensional saturated granular media

In this study, a novel approach to incorporate the pore water pressure in the discrete element method (DEM) to comprehensively model saturated granular media was developed. A numerical model was constructed based on the DEM by implanting additional routines in the basic DEM code; pore water pressure calculations were used with a two-dimensional (2D) model to simulate the undrained behavior of saturated granular media. This model coupled the interaction of solid particles and the pore fluid in saturated granular media. Finally, several 2D undrained shear tests were simulated. The test results showed that the model could predict the response of the saturated granular soil to shear loading. The effect of initial compaction was investigated. Biaxial tests on dense and loose specimens were conducted, and the effect of the initial density on the change in shear strength and the volume change of the system was investigated. The overall behavior of loose and dense specimens was phenomenologically similar to the real granular material. Constant volume tests were simulated, and the results were compared to those from the coupled model. Induced anisotropy was micromechanically investigated by studying the contact force orientation. The change in anisotropy depended on the modeling scheme. However, the overall responses of the media obtained using the coupled and constant volume methods were similar.     

Discrete simulation and micromechanical analysis of two-dimensional saturated granular media

In this study, a novel approach to incorporate the pore water pressure in the discrete element method （DEM） to comprehensively model saturated granular media was developed. A numerical model was constructed based on the DEM by implanting additional routines in the basic DEM code; pore water pressure calculations were used with a two-dimensional （2D） model to simulate the undrained behavior of satu- rated granular media. This model coupled the interaction of solid particles and the pore fluid in saturated granular media. Finally, several 2D undrained shear tests were simulated. The test results showed that the model could predict the response of the saturated granular soil to shear loading. The effect of initial compaction was investigated. Biaxial tests on dense and loose specimens were conducted, and the effect of the initial density on the change in shear strength and the volume change of the system was inves- tigated. The overall behavior of loose and dense specimens was phenomenologically similar to the real granular material. Constant volume tests were simulated, and the results were compared to those from the coupled model. Induced anisotropy was micromechanically investigated by studying the contact force orientation. The change in anisotropy depended on the modeling scheme. However, the overall responses of the media obtained usinz the couoled and constant volume methods were similar.     

Assessment of assembly shape and geometric effects on amplification of shear waves using discrete element method

The discrete element method (DEM) is a capable tool used to simulate shear wave propagation in granular assemblies for many years. Researchers have studied assembly shapes such as rectangles (in 2D simulations) or cylinders and cubes (in 3D simulations). This paper aimed to qualify the effect of assembly shape on the shear wave propagation and maximum amplification in the vertical plane (horizontal and vertical directions) caused by this propagation. To this end, shear wave propagations in different assembly shapes such as rectangle, trapezium, and triangle with rigid boundary conditions were simulated. A sine wave pulse was applied with a point source by moving a particle as the transmitter particle. To evaluate the shear wave velocity of the assemblies, the transmitter and receiver particles were simulated. All the simulations were performed with 2D DEM which is a useful tool to determine the amount and location of the maximum amplification factor of the assembly in both horizontal and vertical directions. An advantage of this study was assessing the effect of parameters such as input wave frequency, assembly height, shape, and aspect ratios on the amplification of the input waves.     

Numerical modeling of fluid-particle interaction in granular media

Fluid-particle interaction underpins important behavior of granular media. Particle-scale simulation may help to provide key microscopic information governing the interaction and offer better understanding of granular media as a whole. This paper presents a coupled computational fluid dynamics and discrete element method (CFD-DEM) approach for this purpose. The granular particle system is modeled by DEM, while the fluid flow is simulated by solving the locally averaged Navier–Stokes equation with CFD. The coupling is considered by exchanging such interaction forces as drag force and buoyancy force between the DEM and CFD. The approach is benchmarked by two classic geomechanics problems for which analytical solutions are available, and is further applied to the prediction of sand heap formation in water through hopper flow. It is demonstrated that the key characteristic of granular materials interacting with pore water can be successfully captured by the proposed method.     

An extended finite element method for the simulation of particulate viscoelastic flows

We present an extended finite element method (XFEM) for the direct numerical simulation of the flow of viscoelastic fluids with suspended particles. For moving particle problems, we devise a temporary arbitrary Lagrangian–Eulerian (ALE) scheme which defines the mapping of field variables at previous time levels onto the computational mesh at the current time level. In this method, a regular mesh is used for the whole computational domain including both fluid and particles. A temporary ALE mesh is constructed separately and the computational mesh is kept unchanged throughout the whole computations. Particles are moving on a fixed Eulerian mesh without any need of re-meshing. For mesh refinements around the interface, we combine XFEM with the grid deformation method, in which nodal points are redistributed close to the interface while preserving the mesh topology. Our method is verified by comparing with the results of boundary fitted mesh problems combined with the conventional ALE scheme. The proposed method shows similar accuracy compared with boundary fitted mesh problems and superior accuracy compared with the fictitious domain method. If the grid deformation method is combined with XFEM, the required computational time is reduced significantly compared to uniform mesh refinements, while providing mesh convergent solutions. We apply the proposed method to the particle migration in rotating Couette flow of a Giesekus fluid. We investigate the effect of initial particle positions, the Weissenberg number, the mobility parameter of the Giesekus model and the particle size on the particle migration. We also show two-particle interactions in confined shear flow of a viscoelastic fluid. We find three different regimes of particle motions according to initial separations of particles.     

Determining a representative element volume for DEM simulations of samples with non-circular particles

Numerical studies on the number of particles or system size required to attain a representative element volume (REV) for discrete element method (DEM) simulations of granular materials have almost always considered samples with spherical or circular particles. This study considers how many particles are needed to attain a REV for 2D samples of 2-disc cluster particles where the particle aspect ratio (AR) was systematically varied. Dense and loose assemblies of particles were simulated. The minimum REV was assessed both by considering the repeatability of static packing characteristics and the shearing behaviour in biaxial compression tests, and by investigating the effect of sample size on the measured characteristics and observed shearing behaviour. The repeatability of the data considered generally improved with increasing sample size. The packing characteristics of the dense samples were more repeatable suggesting that the minimum REV reduces with increasing packing density. The minimum REV was observed to be sensitive to the characteristic measured. Although the overall responses of the samples during shear deformation were similar irrespective of the sample sizes, the smaller the sample size, the higher the fluctuations observed in the responses. Analysis of the coefficient of variation of the fluctuations around the critical state stress ratio can provide insight as to whether a REV is attained. The particle AR influences the effect of sample size on shearing characteristics and thus the minimum number of particles required to attain a REV; this can be explained by the influence of AR on the number of contacts within the samples.     

Prediction of draft forces in cohesionless soil with the Discrete Element Method

The Discrete Element Method (DEM) is applied to predict draft forces of a simple implement in cohesionless granular material. Results are compared with small-scale laboratory tests in which the horizontal force is measured at a straight blade. This study is focused on the case of cohesionless material under quasi-static conditions.The DEM requires the calibration of the local contact parameters between particles to adjust the bulk material properties. The most important bulk property is the angle of internal friction ?. In the DEM, the shear resistance is limited in the case of spherical particles due to excessive particle rotations. This is cured by retaining rotations of the particles. Although this is known to prevent the material from developing shear bands, the model still turns out to be capable of predicting the reaction force on the blade.In contrast to empirical formulas for this kind of application, the DEM model can easily be extended to more complex tool geometries and trajectories. This study helps to find a simple and numerically efficient setup for the numerical model, capable of predicting draft forces correctly and so allowing for large-scale industrial simulations.     

Modeling particle sedimentation in viscous fluids with a coupled immersed boundary method and discrete element method

Numerical techniques have increasingly been used to model fluid–particle two-phase flows. Coupling the immersed boundary method (IBM) and discrete element method (DEM) is one promising approach for modeling particulate flows. In this study, IBM was coupled with DEM to improve the reliability and accuracy of IBM for determining the positions of particles during the sedimentation process within viscous fluids. The required ratio of the particle diameter to the grid size (D/dx) was determined by comparing the simulation results with the analytical solution and experimental data. A dynamic mesh refinement model was utilised in the IBM model to refine the computational fluid dynamics grid near the particles. In addition, an optimum coupling interval between the IBM and DEM models was determined based on the experimental results of a single particle sedimentation within silicon oil at a Reynolds number of 1.5. The experimental results and the analytical solution were then utilised to validate the IBM–DEM model at Reynolds numbers of 4.1, 11.6, and 31.9. Finally, the validated model was utilised to investigate the sedimentation process for more than one particle by modeling the drafting-kissing-tumbling process and the Boycott phenomenon. Benchmark tests showed that the IBM–DEM technique preserves the advantages of DEM for tracking a group of particles, while the IBM provides a reliable and accurate approach for modeling the particle–fluid interaction.     

利用连接尺度方法的颗粒材料破碎数值分析

基于已有的颗粒材料连接尺度方法(BSM)[1-2],发展了在细尺度上采用离散颗粒集合体模型与离散单元法(DEM)并引入了颗粒破碎模型,而在粗尺度上采用Cosserat连续体模型与有限单元法(FEM)的BSM。仅在有限局部区域内采用DEM从细观层次关注颗粒材料破碎现象,而在全域采用储存空间和花费时间较少的FEM,同时在粗细两个尺度采用不同的时间步长。讨论了颗粒材料发生破碎时,颗粒材料结构的承载能力与微结构的演变。数值算例结果说明了所提出可模拟破碎的BSM的可用性和优越性,以及颗粒破碎对颗粒材料微观力学行为的影响。