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.     

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.     

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.     

粒料固有各向异性的离散元模拟与细观分析

料在摊铺后形成的颗粒定向排列将导致其力学性质的固有各向异性. 依据粒料的实际不规则形状, 构建了可模拟粒间咬合嵌挤作用的三维离散元复杂形状颗粒; 生成了5 种不同沉积方向的各向异性试件和1种各向同性试件, 对比了各试件在三轴压缩试验中的宏观力学特性差异; 引入组构张量以量化各向异性程度, 利用玫瑰图表达接触法向与接触力的分布特征, 探究了粒料各向异性的细观发展规律. 结果表明: 颗粒长轴愈趋向水平排布, 峰值应力比愈大, 剪缩与剪胀程度愈明显; 相较于各向同性试件, 沉积角$\theta$为料在摊铺后形成的颗粒定向排列将导致其力学性质的固有各向异性. 依据粒料的实际不规则形状, 构建了可模拟粒间咬合嵌挤作用的三维离散元复杂形状颗粒; 生成了5 种不同沉积方向的各向异性试件和1种各向同性试件, 对比了各试件在三轴压缩试验中的宏观力学特性差异; 引入组构张量以量化各向异性程度, 利用玫瑰图表达接触法向与接触力的分布特征, 探究了粒料各向异性的细观发展规律. 结果表明: 颗粒长轴愈趋向水平排布, 峰值应力比愈大, 剪缩与剪胀程度愈明显; 相较于各向同性试件, 沉积角$\theta$为$0^\circ$时, 峰值应力比和最大体积压缩应变分别提高了12.6\%和18.8\%, 其原因在于加载过程中颗粒旋转和滑动百分比更小, 内部调整时间更短、更易被剪密; 固有各向异性对颗粒法向接触力分布的影响不大, 但显著影响接触法向分布特征; 剪切过程中, $\theta$为$90^\circ$时的接触法向各向异性系数先快速减小后逐渐增大, 而$\theta$为$0^\circ$到$60^\circ$时则呈现出增大后稍有回落或趋于稳定的趋势, 且$\theta$ 愈小的试件各向异性系数增大愈快.     

Structural changes in granular materials: The case of irregular polygonal particles

This paper presents a series of numerical simulations of biaxial tests performed on assemblies of two-dimensional irregular polygonal particles. Each sample is prepared with a technique similar to dry pluviation. Different aspect ratios (1–3) are considered and the behavior of granular samples is analyzed from both a global and a local point of view. More precisely, the influence of the particle aspect ratio on both inherent (initial) and induced anisotropy is investigated. New internal variables which are related to the orientation of particles are proposed. They give new insight into the specific mechanisms that control the behavior of irregular polygonal materials. Associated to global variables, they demonstrate the existence of a critical state irrespective of the investigated aspect ratios. However, for materials with higher aspect ratios (2 and 3), their inherent anisotropy prevents any extensive reorganization, this means that, within the range of usual strains considered in civil engineering, the particle reorientation remains in progress and considerable deformations are required to reach the critical state.     

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.     

Multi-slip gradient formulation for modeling microstructure effects on shear bands in granular materials

This paper presents a higher order gradient multi-slip formulation to model the effect of inhomogeneous deformation in granular materials. The effects of heterogeneity and porosity anisotropy within the multi-slip formulation are taken into consideration through the modification of the mobilized friction. The mobilized friction is assumed to be a direct function of either the gradient of the porosity distribution or the fabric tensor. The formulation with two active slip planes was implemented into a finite element code and used to simulate biaxial shear tests on dry sand. The analysis quantifies most of the shear band characteristics observed by past experimentation. It is shown that the localization and shear band characteristics in granular materials are very much dependent on the initial fabric and slip system arrangement.     

Multiscale approaches to describe mechanical responses induced by particle removal in granular materials

A multi-scale approach is proposed as an attempt to assess internal erosion induced effects on the mechanical properties of a granular medium. When submitted to internal flow, some particles tend to be removed from the initial granular material, leading to drastic changes in the microstructure. The medium is changing during its lifetime and this cannot be ignored in attempts to model its mechanical behavior. In this first analysis, the degradation of a granular assembly is simulated through the progressive removal of its finest particles for both isotropic and anisotropic stress states. A discrete element model as well as an analytical micromechanical model are used to compare induced deformations and properties changes along the erosion process. The shear strength and flow properties appear to be strongly modified by this extraction phase: both models highlight a change from a dilatant to a contractant behavior with the degradation as well as strong dependency of the overall stability on the mobilized friction level. In particular, results show how failure is triggered when particle removal occurs for mobilized frictions greater than the one reached at the critical state.     

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.     

Modeling interface shear behavior of granular materials using micro-polar continuum approach

Recently, the authors have focused on the shear behavior of interface between granular soil body and very rough surface of moving bounding structure. For this purpose, they have used finite element method and a micro-polar elasto-plastic continuum model. They have shown that the boundary conditions assumed along the interface have strong influences on the soil behavior. While in the previous studies, only very rough bounding interfaces have been taken into account, the present investigation focuses on the rough, medium rough and relatively smooth interfaces. In this regard, plane monotonic shearing of an infinite extended narrow granular soil layer is simulated under constant vertical pressure and free dilatancy. The soil layer is located between two parallel rigid boundaries of different surface roughness values. Particular attention is paid to the effect of surface roughness of top and bottom boundaries on the shear behavior of granular soil layer. It is shown that the interaction between roughness of bounding structure surface and the rotation resistance of bounding grains can be modeled in a reasonable manner through considered Cosserat boundary conditions. The influence of surface roughness is investigated on the soil shear strength mobilized along the interface as well as on the location and evolution of shear localization formed within the layer. The obtained numerical results have been qualitatively compared with experimental observations as well as DEM simulations, and acceptable agreement is shown.     

A discrete element model for simulating saturated granular soil

A numerical model is developed to simulate saturated granular soil, based on the discrete element method. Soil particles are represented by Lagrangian discrete elements, and pore fluid, by appropriate discrete elements which represent alternately Lagrangian mass of water and Eulerian volume of space. Macro-scale behavior of the model is verified by simulating undrained biaxial compression tests. Micro-scale behavior is compared to previous literature through pore pressure pattern visualization during shear tests. It is demonstrated that dynamic pore pressure patterns are generated by superposed stress waves. These pore-pressure patterns travel much faster than average drainage rate of the pore fluid and may initiate soil fabric change, ultimately leading to liquefaction in loose sands. Thus, this work demonstrates a tool to roughly link dynamic stress wave patterns to initiation of liquefaction phenomena.     

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.     

椭球颗粒体系剪切过程中自由体积的分布与演化

颗粒材料是一种复杂的多体相互作用体系, 由大量离散的颗粒和其周围的自由体积组成. 虽然颗粒的自由体积与颗粒材料的力学性能和变形特征的相关性已得到证实, 但是由于表征上的困难, 目前对非球形颗粒体系的局部自由体积的认识还不够充分. 本文采用连续离散耦合分析方法进行了不同主轴长度的椭球颗粒试样的三轴剪切数值模拟, 基于Set Voronoi算法对剪切过程中的颗粒试样进行了Voronoi元胞分割, 分析了颗粒试验在剪切过程中自由体积的统计分布特性和演化规律, 研究了颗粒形态对自由体积的影响. 剪切过程中Voronoi元胞的各向异性逐渐增强, 且各项异性增强程度随颗粒非球度的增加而增大, 表明非球颗粒在剪切过程中经历更加强烈的重排列. 具有不同非球度的椭球颗粒体系的局部孔隙比均服从k?Γ分布, 且这个分布仅与颗粒体系的全局孔隙比相关, 不受颗粒形态和剪切状态的影响. 局部孔隙比的波动呈现非对称拉普拉斯分布, 非对称参数刻画了局部自由体积收缩和膨胀的博弈, 其与全局孔隙比呈线性关系.      

Effect of particle breakage on the behavior of simulated angular particle assemblies

Many attempts have been made to find various relationships for different parameters and some kinds of constitutive models for studying the behavior of particulate media. All these models are based on concepts of continuous media. Using a numerical method such as discrete element method, one can figure out what is happening through a discontinuous media where soil particles play the main role in introducing the shear strength and deformation characteristics. The behavior of the media with breakable particles is studied in this paper and compared with that of the assembly with non-breakable particles. In this paper, the hyperbolic elastic model is investigated for the assembly of polygon shaped particles in two different test series. In addition, evolution of different macro parameters of the assembly such as volume strain, angle of friction, angle of dilatancy and elastic modulus is studied during the simulation tests both for non-breakable and breakable soil particles. At the end, a parametric study is performed on the effect of strength of particle breakage on the assembly behavior.     

On the rheology of dilative granular media: Bridging solid- and fluid-like behavior

A new rate-dependent plasticity model for dilative granular media is presented, aiming to bridge the seemingly disparate solid- and fluid-like behavioral regimes. Up to date, solid-like behavior is typically tackled with rate-independent plasticity models emanating from Mohr–Coulomb and Critical State plasticity theory. On the other hand, the fluid-like behavior of granular media is typically treated using constitutive theories amenable to viscous flow, e.g., Bingham fluid. In our proposed model, the material strength is composed of a dilation part and a rate-dependent residual strength. The dilatancy strength plays a key role during solid-like behavior but vanishes in the fluid-like regime. The residual strength, which in a classical plasticity model is considered constant and rate-independent, is postulated to evolve with strain rate. The main appeal of the model is its simplicity and its ability to reconcile the classic plasticity and rheology camps. The applicability and capability of the model are demonstrated by numerical simulation of granular flow problems, as well as a classical shear banding problem, where the performance of the continuum model is compared to discrete particle simulations and physical experiment. These results shed much-needed light onto the mechanics and physics of granular media at various shear rates.     

The role of a realistic volume constraint in modelling a two dimensional granular assembly

This paper presents a deceptively simple mathematical model for the deformation of granular materials composed of rigid particles. The model captures many of the diverse features of the behaviour of such a material and emphasises the importance of volume constraints in situations where the deformation is mainly by particle rearrangement. It is constructed using a simple dissipation function and a rather more complicated dilatancy rule containing an updateable reference strain. This allows the solid-like and fluid-like properties of granular materials to be reconciled in a single model.The model has been used to simulate experiments that use an analogue of an ideal granular material [Joer, H.A., Lanier, J., Fahey, M., 1998. Deformation of granular materials due to rotation of principal axes. Geotechnique 48 (5), 605-619] consisting of a two dimensional assembly of thin PVC rods. These experiments clearly illustrate: partially reversible dilatancy in direct shear tests; cyclic shearing leading to liquefaction in constant volume shear tests; and non-coaxiality of the principal axes of stress and strain increment in circular loading tests. These radically different modes of deformation provide a challenging data set that allows the model's potential to be clearly demonstrated.The authors believe that the comparison of these experimental results and our simulations give strong support to the assertion that volume changes associated with shear deformation are responsible for the rotational kinematic hardening seen in granular materials, and hence, the non-coaxiality of the stress and strain-rate tensors.     

Three-dimensional constitutive relations for granular materials based on the dilatant double shearing mechanism and the concept of fabric

An extension of a three-dimensional model proposed by Anand and Gu (2000) for amorphous granular materials to include the effects of initial and induced anisotropy is presented in this paper. The proposed model can also be considered as a three-dimensional generalization of a model recently developed by Zhu et al. (2005) for the planar deformation of granular materials. The main ingredients of the model include the dilatant double shearing mechanism (Spencer, 1964, Mehrabadi and Cowin, 1978), the concept of fabric (Oda, 1972), and an extension of the Mohr–Coulomb yield criterion (Shield, 1955, Spencer, 1982) to three dimensions.The constitutive equations are implemented in the finite element program ABAQUS/Explicit (ABAQUS, 2001) by developing a user-material subroutine to conduct numerical triaxial compression tests for samples of granular materials with different initial anisotropy. The numerical results agree with the observed behavior and show that the extended constitutive model is capable of capturing the strength anisotropy of granular materials. Employing the anisotropic model developed here, we have also repeated the numerical simulation of the stress state in a static conical sand pile conducted earlier by Anand and Gu (2000). We find that fabric has little or no influence on the vertical stress distribution except at the base of the sand pile where the peak value of this stress is slightly higher than that predicted by the model of Anand and Gu (2000) which does not include the effects of fabric. We also find that the direction of the principal compressive stress changes from vertical at points away from the center of the pile to almost horizontal at points close to the center of the pile. This result provides a possible explanation for the observed dip in the vertical stress distribution in sand piles.     

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.     

Modeling shear localization along granular soil–structure interfaces using elasto-plastic Cosserat continuum

The current study presents finite element simulations of shear localization along the interface between cohesionless granular soil and bounding structure under large shearing movement. Micro-polar (Cosserat) continuum approach is applied in the framework of elasto-plasticity in order to overcome the numerical problems of localization modeling seen in the conventional continuum mechanics. The effects of different micro-polar kinematic boundary conditions, along the interface, on the evolution and location of shear band are shown by the numerical results. Furthermore, shear band thickness is also investigated for its dependence on the initial void ratio, vertical pressure and mean grain size. Here, the distribution and evolution of static and kinematic quantities are the main focuses regarding infinite layer of micro-polar material during plane shearing, especially with advanced large movement of bounding structure. The influence of such movement has not been investigated yet in the literature. Based on the results obtained from this study, shear localization appears parallel to the direction of shearing. It occurs either in the middle of granular layer or near boundaries, regarding the assumed micro-polar kinematic boundary conditions at the bottom and top surfaces of granular soil layer. Narrower shear band is observed in lower rotation resistance of soil particles along the interface. It is emphasized that the displacement magnitude of bounding structure has significant effect on the distribution and evolution of state variables and polar quantities in the granular soil layer. However, continuous displacement has no meaningful effect on the thickness of shear band. Here, smooth distributions of void ratio and shear stress components are obtained within the shear band, what the other previous numerical investigations did not receive. Despite indirect linking of Lade’s model to the critical state soil mechanics, state variables tend towards asymptotical stationary condition in large shear deformation.