Micromechanical Modeling of Yield in Isotropic Non-Cohesive Particulate Materials

We perform a micromechanical analysis of general isotropic non-cohesive particulate materials idealized as three-dimensional random assemblies of uniform spheres with a simple linear elastic inter-particle contact force law and inter-particle Coulomb friction law. We obtain analytical relationships between the inter-particle friction coefficient \(\mu\) (or inter-particle friction angle \(\phi _\mu = \tan ^{-1} \mu\)) on the microscale and the material friction angle \(\phi\) on the macroscale. Our micromechanical analysis directly employs force and moment equilibrium (together with compatibility and the contact constitutive assumptions noted) rather than energy methods, and thus can account for the effects of particle rotation, and in particular the effects of mechanisms or zero-energy modes due to particle rotation. To explore the effects of particle rotation, we perform analyses with particle rotation either allowed or prohibited. To validate the analytical results obtained here, we compare the \(\phi\) versus \(\phi _\mu\) curves determined theoretically to those obtained by the discrete element method (DEM) for six randomly packed specimens of 3430–29, 660 uniform spherical elements with uniform inter-element Coulomb friction in Fleischmann et al. in Geotech Geol Eng 32(4):1081–1100, (2014). The \(\phi\) versus \(\phi _\mu\) curves derived here show remarkable agreement with those obtained via DEM simulations in Fleischmann et al. in Geotech Geol Eng 32(4):1081–1100, (2014), especially for the case in which particle rotation is not artificially restrained.     

A non-linear numerical model for soil mechanics

A new computer program (CONBAL-2) is developed for 2D numerical simulations of granular soil by random arrays of spheres. CONBAL-2 uses the discrete-element method and is based on 3D program TRUBAL, previously presented by Cundall. As in TRUBAL, the new program models a random array of elastic spheres in a periodic space. The main modification of TRUBAL is the implementation by the authors of a rigorous solution for the force–displacement relation at the interparticle contacts. This force-displacement relation is a function of the elastic constants, friction coefficient and sizes of the spheres, with the properties of quartz used to simulate sand. Other specific features of CONBAL-2 include its 2D character, the lack of particle rotation and its capability to simulate shear loading on any plane. Simulated laboratory test results are presented using CONBAL-2 and several random arrays of 531 spheres having two particle sizes. These simulations include monotonic loading drained and undrained (constant volume) ‘triaxial’ experiments, as well as a cyclic-loading, constant-volume ‘torsional shear’ test. The stress–strain curves, effective stress paths, volume changes, as well as the ‘pore water pressure’ build-up behaviour obtained in the simulations compare favourably—qualitatively and in some aspects quantitatively—with similar laboratory results on sands. However, the simulated soil is somewhat stiffer and stronger due to the perfectly rounded particles, limited range of grain sizes, lack of particle rotation and 2D character of the model.     

Effect of particle friction and polydispersity on the macroscopic stress–strain relations of granular materials

The macroscopic mechanical behavior of granular materials inherently depends on the properties of particles that compose them. Using the discrete element method, the effect of particle contact friction and polydispersity on the macroscopic stress response of 3D sphere packings is studied. The analytical expressions for the pressure, coordination number and fraction of rattlers proposed for isotropically deformed frictionless systems also hold when the interparticle coefficient of friction is finite; however, the numerical values of the parameters such as the jamming volume fraction change with varying microscopic contact and particle properties. The macroscopic response under deviatoric loading is studied with triaxial test simulations. Concerning the shear strength, our results agree with previous studies showing that the deviatoric stress ratio increases with particle coefficient of friction μ starting from a nonzero value for μ = 0 and saturating for large μ. On the other hand, the volumetric strain does not have a monotonic dependence on the particle contact friction. Most notably, maximum compaction is reached at an intermediate value of the coefficient of friction μ ≈ 0.3. The effect of polydispersity on the macroscopic stress–strain relationship cannot be studied independent of initial packing conditions. The shear strength increases with polydispersity when the initial volume fraction is fixed, but the effect of polydispersity is much less pronounced when the initial pressure of the packings is fixed. Finally, a simple hypoplastic constitutive model is calibrated with numerical test results following an established procedure to ascertain the relation between particle properties and material coefficients of the macroscopic model. The calibrated model is in good qualitative agreement with simulation results.     

3D analysis of kinematic behavior of granular materials in triaxial testing using DEM with flexible membrane boundary

This paper describes the constitutive behavior and particle-scale kinematics of granular materials in three-dimensional (3D) axisymmetric triaxial testing using discrete element method (DEM). PFC3D code was used to run the DEM simulations using a flexible membrane boundary model consisting of spherical particles linked through flexible contact bonds. The overall deformation behavior of the specimen was then compared with the specimen with rigid boundary and experimental measurements. Computed tomography was used to track the evolution of particle translation and rotation within a laboratory triaxial specimen in 3D. The DEM model of the flexible membrane specimen successfully predicted the stress–strain behavior when compared with laboratory experiment results at different confining pressures. The DEM results showed that the rigid specimen applies a uniform deformation and leads to non-uniformities in the confining stress along the particle-boundary interface in the lateral direction. In contrast, the flexible specimen better replicates the uniformly applied confining stress of a laboratory triaxial experiment. The 3D DEM simulations of the specimen with flexible membrane overpredicted particle translation and rotation in all directions when compared to a laboratory triaxial specimen. The difference between the particle translation and rotation distributions of DEM specimens with rigid and flexible membrane is almost negligible. The DEM specimen with flexible membrane produces a better prediction of the macroscopic stress–strain behavior and deformation characteristics of granular materials in 3D DEM simulations when compared to a specimen with rigid membrane. Comparing macroscale response and particle-scale kinematics between triaxial simulation results of rigid versus flexible membrane demonstrated the significant influence of boundary effects on the constitutive behavior of granular materials.     

Three-dimensional random network model for thermal conductivity in particulate materials

This paper describes a three-dimensional random network model to evaluate the thermal conductivity of particulate materials. The model is applied to numerical assemblies of poly-dispersed spheres generated using the discrete element method (DEM). The grain size distribution of Ottawa 20–30 sand is modeled using a logistic function in the DEM assemblies to closely reproduce the gradation of physical specimens. The packing density and inter-particle contact areas controlled by confining stress are explored as variables to underscore the effects of micro- and macro-scales on the effective thermal conductivity in particulate materials. It is assumed that skeletal structure of 3D granular system consists of the web of particle bodies interconnected by thermal resistor at contacts. The inter-particle contact condition (e.g., the degree of particle separation or overlap) and the particle radii determine the thermal conductance between adjacent particles. The Gauss–Seidel method allows evaluation of the evolution of temperature variation in the linear system. Laboratory measurements of thermal conductivity of Ottawa 20–30 sand corroborate the calculated results using the proposed network model. The model is extended to explore the evolution of thermal conduction depending on the nucleation habits of secondary solid phase as an anomalous material in the pore space. The proposed network model highlights that the coordination number, packing density and the inter-particle contact condition are integrated together to dominate the heat transfer characteristics in particulate materials, and allows fundamental understanding of particle-scale mechanism in macro-scale manifestation.     

砂堆底部压力特性室内实验和PFC2D数值模拟分析

采用室内砂雨试验和颗粒流软件PFC2D,研究了颗粒粒径和摩擦系数对砂粒堆积结构压力特性的影响。室内试验是通过自制砂雨装置,选取不同级配的干砂、湿砂、湿砂+黄土,研究颗粒粒径和颗粒间摩擦系数对砂堆底部压力的影响;数值模拟利用PFC软件内置命令（ball generate）生成3种不同粒径的颗粒,通过ball property命令给颗粒间赋予6种不同的摩擦系数,从力链角度分析其对砂堆结构内部的影响。结果表明:颗粒粒径和颗粒间摩擦系数对砂堆结构底部压力具有显著的影响;相同摩擦系数条件下,颗粒粒径越大,堆积结构内部力链分布越稀疏,其底部压力越小;相同粒径条件下,砂堆底部压力随摩擦系数的增大呈先减小后趋于稳定的趋势,即摩擦系数对砂堆底部压力的影响存在上限值。     

Exploring the influence of interparticle friction on critical state behaviour using DEM

Understanding the extent to which discrete element method (DEM) simulations can capture the critical state characteristics of granular materials is important to legitimize the use of DEM in geomechanics. This paper documents a DEM study that considered the sensitivity of the critical state response characteristics to the coefficient of interparticle friction (μ) using samples with gradings that are representative of a real soil. Most of the features that are typically associated with sand behaviour at the critical state were seen to emerge from the DEM simulation data. An important deviation occurs when high μ values (μ ≥ 0.5) are used, as has been the case in a number of prior DEM studies. While there is a systematic variation in the critical state behaviour with μ for μ < 0.5, when μ ≥ 0.5, the behaviour at the critical state seems to be insensitive to further increases in μ. In contrast to observations of conventional soil response, when μ ≥ 0.5, the void ratio at the critical state initially increases with increasing mean effective stress (p′). Analysis of the DEM data and use of simple models of isolated force chains enabled some key observations. When ‘floating’ particles that do not transmit stress are eliminated from the void ratio calculation, the void ratio at the critical state decreases consistently with increasing p′. There is a transition from sliding to rolling behaviour at the contact points as μ increases. Beyond a limiting value of μ, further increases in μ do not increase the buckling resistance of individual strong force chains. Copyright © 2014 John Wiley & Sons, Ltd.     

Macroscale friction of granular soils under monotonic and cyclic loading based upon micromechanical determination of dissipated energy

Macroscopic frictional behavior of granular materials is of great importance for studying several complex problems associated with fault slip and landslides. The main objective of this study is to model the macroscale frictional behavior of granular soils under monotonic and cyclic loadings based upon micromechanical determination of dissipated energy at particle contacts. This study is built on the general observation that the externally computed energy dissipation should be equal to the total internal energy dissipation derived from inter-particle sliding and rolling, energy losses from inter-particle collisions, and damping. For this purpose, the discrete element method is used to model a granular soil and determine the stored, dissipated, and damping energies associated with shear loading for applied monotonic and cyclic velocities. These energies are then related to the friction by an application of the Taylor-critical state power balance relationship. Also, the contributions of the different modes of energy dissipation (normal, shear, and rolling) to the total frictional resistance were studied. By changing the inter-particle friction, the simulations showed that the macroscopic friction was nearly constant, the slip friction increased almost linearly with increasing inter-particle friction, and the difference between the two was attributed to the non-energy dissipating dilatancy component. By providing a clear relationship between energy dissipated by micro-scale mechanisms versus the traditional engineering definition based on macro-scale (continuum) parameters, this study provides a means to develop a better understanding for the frictional behavior of granular media.

     

DEM modeling of aging or creep in sand based on the effects of microfracturing of asperities and evolution of microstructural anisotropy during triaxial creep

Aging- or creep-related phenomena in sand have been widely studied, and the discrete element method (DEM) has been frequently used to model the associated soil behavior and then to explore the associated underlying mechanisms. However, several difficulties involved in modeling still remain unsolved. To resolve these difficulties, a new approach based on the effect of the microfracturing of asperities is proposed in this study for the DEM modeling of the sand aging or creep process through several aging cycles of associated reduction in the mobilized friction resistance at particle contacts and subsequent particle rearrangement to reach a new equilibrium state. This approach can be easily incorporated into different contact models and DEM simulations of the loading, unloading, and/or reloading processes, in either drained or undrained conditions, before and/or after aging. This new approach is proven effective because the DEM simulations incorporated with this new approach can satisfactorily reproduce the experimental observations in the triaxial creep process, drained and undrained recompression after aging, and 1D secondary compression and rebound. The simulation results also indicate that, based on the stress–force–fabric relationship, the contribution from the contact normal anisotropy to the deviatoric stress q gradually increases, whereas the contribution from the tangential force anisotropy becomes less during triaxial creep under a constant q. Moreover, the contacts between particles are gradually away from the state where the frictional resistance is fully mobilized, and then become more stable. During the subsequent triaxial recompression after creep, the aged samples exhibit enhanced soil stiffness, which is also found to be associated with the evolution of the invariants of the anisotropy tensors. It is worthwhile noting that the aging or creep effects on the microstructural changes, e.g., the invariants of the anisotropy tensors, can be gradually erased upon further recompression. This explains why the stress–strain responses of the aged samples during recompression gradually rejoin the original stress–strain response obtained from the sample without being subjected to aging or creep.     

Predictability of long runout landslide motion: implications from granular flow mechanics

Quality of landslide motion prediction is directly linked to the understanding of the basic flow mechanisms. Although it is known that landslides are granular mass flows and granular flow mechanics is an established area of research, hypotheses on landslide motion are still based on simple geometrical relations and heuristic assumptions. New insights into the development of flow properties of high-speed, high-concentration granular flows are given by results of discrete particle simulations: rapid granular flows are self-organizing dynamic systems that are forced to develop a plastic body rheology. This behaviour must be described by a coefficient of internal friction μ_CM that refers to the center of mass of a flow. Coefficients of rapid granular flows of inelastic and rough particles, which are typical for common rock materials, do not vary significantly around μ _CM ≈0.45 that is definitively smaller than the friction coefficient of soil creep (≈0.6). The motion of the center of mass is superimposed by the spreading of the granular mass that is controlled by the same plastic body rheology. This combined motion is a scale-invariant self-similar process that depends only on the drop height of a landslide and its volume. This allows specification of implications that must be given special attention in the development of future models for landslide prediction.     

Improving constant-volume simulations of undrained behaviour in DEM

In order to simulate undrained conditions using the discrete element method, a constant sample volume is often assumed. There are well-recognised problems with these constant-volume triaxial simulations, particularly of dense samples, which inhibit quantitative comparison with laboratory experiments. In this paper, four possible explanations for these problems with conventional constant-volume simulations of ideal spherical particles are explored, each of which has a physical basis: particle crushing, the presence of highly compressible air within the sample, or the reduction in stiffness due to particle surface asperities or non-spherical particle shapes. These options are explored independently and in combination through implementation in the open-source LAMMPS code. In situations where a significant amount of particle crushing occurs, it is important to incorporate this in the simulations so that stresses are not over-estimated. There is experimental evidence that irregular particles have lower Young’s moduli than the Hertzian spheres often used in DEM. In the absence of particle crushing, the most effective method to achieve more realistic stress–strain responses is to reduce the particle shear modulus substantially. This approach has the added computational benefit of enabling an increase in the simulation time-step.

     

Experimental investigation of shear strength of sands with inherent fabric anisotropy

Loading direction-dependent strength of sand has been traditionally characterized in the principal stress space as a direct extension of the Mohr–Coulomb criterion. A recent study found that it is more appropriate to define anisotropic strength of sand on failure/shear planes, but this proposition has only been demonstrated with discrete element method (DEM) simulations. The present study experimentally investigates anisotropic shear strength of sands in this new framework. Three granular materials with distinct grain characteristics ranging from smooth and rounded particles to flaky and angular particles are tested with the bedding plane inclination angle ψ _b varying over the full range of 0°–180°. The main objective is to study how the peak friction angle ? _p of sand is affected by the ψ _b angle and how the ψ _b–? _p relationship evolves with the change of characteristics of constituent sand particles. We find that the general trend of ψ _b–? _p curves for real sands resembles what was predicted by DEM in a previous study, whereas rich anisotropic strength behavior is revealed by the laboratory data. The effects of normal stress and initial density, as well as shear dilation behavior at different shear directions, are also studied.     

颗粒大小对颗粒材料力学行为影响初探

利用一种特殊颗粒材料-玻璃珠进行了一系列室内直剪试验,研究颗粒大小对颗粒材料力学行为的影响。试验一共考虑了3条近乎平行的级配曲线和4种颗粒摩擦情况：干燥状态、水浸润状态、水淹没状态和油浸润状态。试验结果表明,颗粒大小对颗粒材料的力学行为有显著影响,剪胀性随着粒径的增大而增强。为考虑颗粒大小对剪胀性的影响,提出了一种新的剪胀关系式。在该剪胀关系式中,剪胀系数为依赖于颗粒大小和颗粒摩擦等颗粒基本性质的变量。试验研究同时表明临界状态摩擦角随着颗粒大小的增加而增加。此外,从颗粒细观运动的角度提出了颗粒滑动的功能模型,推导出了功能方程,并以此揭示了颗粒大小对临界状态摩擦角影响的细观机制。     

Composite dust grains: Modeling of infrared absorption bands

We investigate the influence of small-scale asphericity of the surfaces of dust grains on the characteristics of the two deepest absorption bands observed in the spectra of protostellar objects and stars (the 3.1 μm water-ice and 9.7 μm silicate bands). The model used has composite particles in the form of radially inhomogeneous spheres with intermediate layers in which the index of refraction changes. The observed band widths and the ratios of the optical depths at the band centers can be explained if the grains are composed of small particles consisting of silicate cores with thin ice mantles and rough surfaces. The grain surface roughness considerably broadens the profile of the silicate band.     

颗粒形态对离散介质剪切强度的影响

自然条件下,颗粒介质大多以非规则单元形态存在。非规则几何形态对颗粒介质的宏观力学性能有很大影响。针对颗粒单元的不同几何形态,采用团颗粒单元对离散介质的直剪试验过程进行了离散元数值计算,详细地讨论了颗粒形态对离散介质剪切强度的影响。该非规则颗粒由不同形态、不同数目、镶嵌尺寸、组合方位和颗粒大小的球形颗粒进行随机构造,其在局部与整体坐标之间的转动、力矩和方位关系通过4元素方法进行确定,基本球体颗粒之间的作用力采用具有Mohr-Coulomb摩擦定侓的Hertz-Mindlin非线性接触模型,并考虑了非线性法向粘滞力的影响。通过构造7种具有相同的质量概率分布的不同形态的团颗粒,在不同法向应力下,对团颗粒的直剪试验进行了离散元模拟,分析了不同形态颗粒的剪切强度。通过对不同形态颗粒介质剪切强度的数值分析,进一步揭示了非规则颗粒间的咬合互锁效应,为分析非规则颗粒的宏观动力特性提供了依据。     

Evaluation of Shear Modulus and Damping Ratio of Granular Materials Using Discrete Element Approach

In this paper, an attempt has been made to highlight the influence of different parameters such as number of cycles, confining pressure, void ratio, gradation, initial anisotropy and stress path on the dynamic properties of granular materials using Discrete Element Method (DEM). A series of strain controlled cyclic triaxial numerical simulations using three dimensional DEM have been carried out on an assembly of spheres. Dynamic properties such shear modulus (G) and damping ratio (D) were determined from the typical hysteresis loop obtained during cyclic triaxial test simulation. It has been observed from the test results that the numerical simulation using DEM has captured the variation of dynamic properties over a wide range of shear strain values for different parameters considered for the current investigation. Maximum shear modulus (G _max) was found to be influenced by initial confining pressure, void ratio, gradation and initial anisotropy. Whereas, the damping ratio (D) was found to be influenced by number of cycles, initial confining pressure, gradation and stress path. Further it has been shown that the variation of shear modulus with shear strain can be divided into three distinct zones such as Isotropic Zone (IZ), Anisotropic Zone (AZ) and Stable Anisotropic Zone (SAZ). A drastic reduction of shear modulus with shear strain has been observed in the Anisotropic Zone (AZ). In addition, the results obtained using numerical simulations have been compared with the laboratory experimental values.     

基于离散元法的砂土破碎演化规律研究

颗粒破碎是影响砂土宏-微观力学性质的重要因素。采用改进型的可破碎颗粒生成方法,通过设置不同强度的平行黏结键模拟不同强度的可破碎颗粒,并借用基于离散元方法（DEM）的双轴压缩试验详细研究了可破碎性土在剪切过程中颗粒破碎率/平均破碎程度、微观尺度上的能量耗散分配机制、剪切破碎带形成以及断裂键各向异性的演化过程。结果表明,颗粒破碎强烈地影响砂土在宏观尺度上的力学响应、颗粒尺度上的能量分配机制以及剪切过程中的颗粒的组织结构演化。颗粒破碎主要影响小应变阶段各能量耗散元的分配机制,而在临界状态下剪切带内的颗粒摩擦以及破碎耗能是消耗外界功的主要因素。数值结果亦表明,颗粒的破碎伴随着整个剪切过程,但破碎率的增长速度却随着剪切应变的发展逐渐降低。另外,在剪切过程中,对于低破碎性土,在临界状态下剪切破碎带基本形成,带内的原有组织结构被打乱,断裂键的各向异性也随之弱化。     

Refined modeling and movement characteristics analyses of irregularly shaped particles

Irregularly shaped (IRS) particles widely exist in many engineering and industrial fields. The macro physical and mechanical properties of the particle system are governed by the interaction between the particles in the system. The interaction between IRS particles is more complicated because of their complex geometric shape with extremely irregular and co‐existed concave and convex surfaces. These particles may interlock each other, making the sliding and friction of IRS particles more complex than that of particles with regular shape. In order to study the interaction of IRS particles more efficiently, a refined method of constructing discrete element model based on computed tomography scanning of IRS particles is proposed. Three parameters were introduced to control the accuracy and the number of packing spheres. Subsequently, the inertia tensor of the IRS particle model was optimized. Finally, laboratory and numerical open bottom cylinder tests were carried out to verify the refined modeling method. The influence of particle shape, particle position, and mesoscopic friction coefficient on the interaction of particles was also simulated. It is noteworthy that with the increase of mesoscopic friction coefficient, the fluidity of IRS particle assembly decreases, and intermittent limit equilibrium state may appear. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical studies of shear banding in interface shear tests using a new strain calculation method

Strain localization is closely associated with the stress–strain behaviour of an interphase system subject to quasi‐static direct interface shear, especially after peak stress state is reached. This behaviour is important because it is closely related to deformations experienced by geotechnical composite structures. This paper presents a study using two‐dimensional discrete element method (DEM) simulations on the strain localization of an idealized interphase system composed of densely packed spherical particles in contact with rough manufactured surfaces. The manufactured surface is made up of regular or irregular triangular asperities with varying slopes. A new simple method of strain calculation is used in this study to generate strain field inside a simulated direct interface shear box. This method accounts for particle rotation and captures strain localization features at high resolution. Results show that strain localization begins with the onset of non‐linear stress–strain behaviour. A distinct but discontinuous shear band emerges above the rough surface just before the peak stress state, which becomes more expansive and coherent with post‐peak strain softening. It is found that the shear bands developed by surfaces with smaller roughness are much thinner than those developed by surfaces with greater roughness. The maximum thickness of the intense shear zone is observed to be about 8–10 median particle diameters. The shear band orientations, which are mainly dominated by the rough boundary surface, are parallel with the zero extension direction, which are horizontally oriented. Published in 2007 by John Wiley & Sons, Ltd.     

Failure criteria for cohesive‐frictional materials based on Mohr–Coulomb failure function

This study presents two three‐parameter failure criteria for cohesive‐frictional materials based on the Mohr–Coulomb failure function. One proposed failure criterion can be linked to Mogi's empirical formula and incorporates the well‐known Von‐Mises, Drucker–Prager, and Linear Mogi criteria as special cases. Another one with smooth and convex cross sections contains a general Lode dependence in the deviatoric plane and includes the Matsuoka–Nakai and Lade–Duncan Lode dependences as special cases. The effect of the intermediate principal stress on the strength of the material can be taken into account in both criteria. The proposed criteria are numerically calibrated against polyaxial data sets of many different types of rocks and concrete_. The comparison results show that the performance of the proposed criteria is excellent, and the failure criterion with a general Lode dependence performs better than the other one for concrete. Copyright © 2015 John Wiley & Sons, Ltd.