Effective simulation of flexible lateral boundaries in two- and three-dimensional DEM simulations

Discrete element method （DEM） models to simulate laboratory element tests play an important role in advancing our understanding of the mechanics of granular material response, including bonded or cemented, particulate materials. Comparisons of the macro-scale response observed in a real physical test and a ＂virtual＂ DEM-simulated test can calibrate or validate DEM models. The detailed, particle scale information provided in the DEM simulation can then be used to develop our understanding of the material behaviour. It is important to accurately model the physical test boundary conditions in these DEM simulations. This paper specifically considers triaxial tests as these tests are commonly used in soil mechanics. In a triaxial test, the test specimen of granular material is enclosed within a flexible latex membrane that allows the material to deform freely during testing, while maintaining a specified stress condition. Triaxial tests can only be realistically simulated in 3D DEM codes, however analogue, 2D, biaxial DEM simulations are also often considered as it is easier to visualize particle interactions in two dimensions. This paper describes algorithms to simulate the lateral boundary conditions imposed by the latex membrane used in physical triaxial tests in both 2D and 3D DEM simulations. The importance of carefully considering the lateral boundary conditions in DEM simulations is illustrated by considering a 2D biaxial test on a specimen of frictional unbonded disks and a 3D triaxial test on a bonded （cemented） specimen of spheres. The comparisons indicate that the lateral boundary conditions have a more significant influence on the local, particle-scale response in comparison with the overall macro-scale observations.     

Quantification of large and localized deformation in granular materials

Quantifying large deformation in granular assemblies using concepts originating from continuum mechanics is a challenging task because of (1) the discontinuous nature of granular displacement, which does not allow the definition of a continuum measure of deformation, and (2) the almost inevitable shear band localization. These problems exist in both real-world granular materials and their numerical idealizations using particle-based simulations. In this work a new method is developed in order to address these issues. Instead of creating a meshed equivalent continuum for quantifying small engineering strains, the new method performs independent random queries on the velocity gradient characteristics of arbitrary sub-domains in the assembly through the novel concept of overlapping reference triangles, thus, enabling rigorous handling of large deformations which are usually associated with localization. The proposed method is illustrated and validated by discrete element method (DEM) simulation of a biaxial compressive test, in which apparent shear banding takes place. The homogenized deformation quantifications based on the new method match the estimations from the imposed boundary conditions. The numerical examples are also applied to (1) quantifying the heterogeneous distribution of deformation over the specimen, (2) visualizing the nucleation process of shear bands, and (3) characterizing shear flow patterns in shear bands. An investigation on the effects of the reference triangle sizes yields some inspiring and practically significant results.     

Stress and fabric in granular material

It has been well recognized that, due to anisotropic packing structure of granular material, the true stress in a specimen is different from the applied stress. However, very few research efforts have been focused on quantifying the relationship between the true stress and applied stress. In this paper, we derive an explicit relationship among applied stress tensor, material-fabric tensor, and force-fabric tensor; and we propose a relationship between the true stress tensor and the applied stress tensor. The validity of this derived relationship is examined by using the discrete element simulation results for granular material under biaxial and triaxial loading conditions.     

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.     

COMPARISON OF GRANULAR MATERIAL BEHAVIOUR UNDER DRAINED TRIAXIAL AND PLANE STRAIN CONDITIONS USING 3D DEM SIMULATIONS

Three dimensional (3D) DEM (discrete element method) simulations of drained triaxial compression and plane strain tests are presented for both dense and loose assemblies of polydisperse spheres using a...     

Implementation of rotational resistance models: A critical appraisal

Contact models that simulate rotational resistance at particle contacts have been proposed as a means to capture the effect of shape in DEM simulations. This contribution critically explores some key issues related to the implementation of rotational resistance models; these include the need for physically meaningful model parameters, the impact of the model on the overall numerical stability/critical time increment for the DEM model, model validation, and the assessment of model performance relative to real physical materials. The discussion is centred around a rotational resistance model that captures the resistance provided by interlocking asperities on the particle surface. An expression for the maximum permissible integration time step to ensure numerical stability is derived for DEM simulations when rotational resistance is incorporated. Analytical solutions for some single-contact scenarios are derived for model validation. The ability of this type of model to provide additional fundamental insight into granular material behaviour is demonstrated using particle-scale analysis of triaxial compression simulations to examine the roles that contact rolling and sliding have on the stability of strong force chains.     

混凝土动态双轴拉压破坏准则细观数值模拟研究

正常服役期内的混凝土结构往往处于复杂应力状态,并且不可避免地会受到偶发动力载荷作用.对于复杂载荷作用下的混凝土力学性能研究,破坏准则是基础.受试验设备等条件限制,现有的动态双轴拉压破坏准则形式复杂、缺乏更高应变率和侧应力比范围且尚未综合考虑应变率和侧应力比的耦合作用.为进一步提出适用范围更高且更准确的混凝土动态双轴拉压...     

Influence of grading on the undrained behavior of granular materials

This paper aims at investigating the influence of grading on the undrained behavior of granular materials. Series of undrained triaxial tests were carried out with two different materials, glass balls and Hostun sand. For each material, samples with different gradings and similar relative densities were prepared. Experimental results show that the undrained shear strength decreases when the coefficient of uniformity C_u=d₆₀/d₁₀$C_{\mathrm{u}}=d_{60}/d_{10}$ increases from 1.1 to 20. The conditions of instability for the selected granular materials were also analyzed, based on the sign of the second-order work during undrained triaxial loading. The results demonstrate a significant influence of the grading: increasing the coefficient of uniformity heightens the potential of static liquefaction and the materials become more unstable. Furthermore, numerical tests using the three-dimensional discrete element method (DEM) were conducted on assemblies of spheres. The DEM inter-particle parameters were derived from the experimental test results on glass balls. The DEM simulations showed similar behaviors compared to experimental results and confirmed the influence of the grain size distribution on the stress–strain relationship and instability phenomena.     

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.     

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.     

含抗转能力散粒体的宏微观力学特性数值分析

颗粒间抗转因素对散粒体宏微观特性的影响十分显著。本文采用离散单元法,将蒋明镜（2005年）等提出的抗转动接触模型植入离散元分析软件PFC2D中,模拟了大量刚性边界下的砂土双轴试验,研究了压缩过程中试样的宏观力学性质及部分微观参量的变化情况,用以定性分析颗粒间抗转作用对散粒体的宏微观力学特性的影响。研究表明：颗粒间的抗转...     

Assessing continuum postulates in simulations of granular flow

Continuum mechanics relies on the fundamental notion of a mesoscopic volume “element” in which properties averaged over discrete particles obey deterministic relationships. Recent work on granular materials suggests that a continuum law may be inapplicable, revealing inhomogeneities at the particle level, such as force chains and slow cage breaking. Here, we analyze large-scale three-dimensional discrete-element method (DEM) simulations of different granular flows and show that an approximate “granular element” defined at the scale of observed dynamical correlations (roughly three to five particle diameters) has a reasonable continuum interpretation. By viewing all the simulations as an ensemble of granular elements which deform and move with the flow, we can track material evolution at a local level. Our results confirm some of the hypotheses of classical plasticity theory while contradicting others and suggest a subtle physical picture of granular failure, combining liquid-like dependence on deformation rate and solid-like dependence on strain. Our computational methods and results can be used to guide the development of more realistic continuum models, based on observed local relationships between average variables.     

Triaxial tension–compression tests for multiaxial cyclic plasticity

The aim of this paper is to present triaxial tension–compression tests and a new triaxial specimen devoted to the study of cyclic plasticity under non-proportional loadings. Because the stress state in the central part of the specimen is not homogeneous, the analysis of the tests need a 3D finite element computation. The constitutive equations used in the simulations have been deduced from complex tension–torsion tests. The comparison between the structural analysis and the experimental results allows to determine the accuracy of the set of constitutive equations and if needed to optimize this set by modifying some hardening rules.     

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.     

MICROMECHANICS ANALYSIS FOR UNSATURATED GRANULAR SOILS

This paper aims at establishing an anisotropic stress expression for unsaturated pendular-state granular soils. Using the second-order fabric tensor, we formulate a micromechanics scheme of soils with statistically averaging method, and reveal that the macroscopic average stress of unsaturated granular soils in pendular-state is not isotropic. Not only is the stress from contact forces anisotropic due to the fabric, but also the capillary stress is directional dependent, which is different from the common point that the capillary stress is isotropic. The capillary stress of unsaturated pendular-state granular soils is determined by the orientation distribution of contact normals, so it is closely related to the initial and induced anisotropy of soils. Finally, DEM numerical simulations of triaxial compression tests of pendular-state soils at different degrees of saturation are used to verify the existence of above anisotropy of stresses.     

DEM investigation on the effect of sample preparation on the shear behavior of granular soil

The effect of initial fabric anisotropy produced by sample preparation on the shear behavior of granular soil is investigated by performing discrete element method (DEM) simulations of fourteen biaxial tests in drained conditions. Numerical test specimens are prepared by three means: gravitational deposition, multi-layer compression, and isotropic compression, such that different initial inherent soil fabrics are created. The DEM simulation results show that initial fabric anisotropy exerts a considerable effect on the shear behavior of granular soil, and that the peak stress ratio and peak dilatancy increase with an increase in the fabric index a_n that is estimated from the contact orientations. The stress–dilatancy relationship is found to be independent of the initial fabric anisotropy. The anisotropy related to the contact orientation and contact normal force accounts for the main contribution to the mobilized friction angle. Also, the occurrence of contractive shear response in an initial shearing stage is accompanied by the most intense particle rearrangement and microstructural reorganization, regardless of the sample preparation method. Furthermore, the uniqueness of the critical state line in e–log p′ and q–p′ plots is observed, suggesting that the influence of initial fabric anisotropy is erased at large shear strains.