Numerical analysis of wave propagation through assemblies of elliptical particles

Wave propagation in granular materials is numerically studied using discrete element simulation. Primary interest is concerned with linking material microstructure with wave propagational behaviors for materials composed of elliptical particles. The discrete element (DEM) scheme uses a nonlinear hysteretic contact law which accounts for differences related to the radius of curvature at the interparticle point of contact. Modeling results yield information on wave speed and amplitude attenuation on two-dimensional, meso-domain model systems of both regular and random assemblies. Particulate models were numerically generated using a biasing scheme whereby partial control of particular fabric measures could be achieved. Three specific fabric measures which were used to characterize the granular material models include branch, contact normal and orienation vectors. DEM simulation results indicated that wave speed and attenuation generally correlated with vector distributions of these fabric variables. A power law relation was proposed between wave speed/attenuation and three averaged projected fabric variables based on orientation, contact normal and branch vectors. Predictions from this specific relation correlated reasonably well with DEM results.     

Three dimensional steady-state locus for dry monodisperse granular materials: DEM numerical results and theoretical modelling

In this paper, steady-state conditions for ideal monodisperse dry granular materials are both theoretically and numerically analysed. A series of discrete element (DEM) numerical simulations have been performed on a periodic cell by imposing stress paths characterized by different Lode angles, pressures, and deviatoric strain rates. The dependence of the material response on both inertial number and loading path has been discussed in terms of void ratio, fabric, and granular temperature. DEM numerical results have been finally compared with the prediction of an already conceived model based on both kinetic and critical state theories, here suitably modified to account for three-dimensional conditions.     

颗粒材料的微观临界状态理论模型

存在临界状态是颗粒材料的一个重要特性。基于孔隙胞元的颗粒离散元方法对二维颗粒体进行双轴加载数值试验,在详细分析数值模拟结果的基础上,从微观几何组构的角度揭示了临界状态的存在机制。基于剪胀性原理,提出了以接触价键表征的微观临界状态理论模型,得到了接触价键与塑性剪切应变的关系表达式,理论模型的结果和二维离散元数值模拟得到的结果吻合较好。通过比较不同情况下数值结果和理论模型中的参数,得到以下结论：表征微观临界状态的参数（临界接触价键和达到临界状态所需要的塑性剪切应变）依赖于颗粒体的微观特性,如颗粒形状、表面摩擦性质、颗粒体的围压和初始孔隙比。     

An efficient flexible membrane boundary condition for DEM simulation of axisymmetric element tests

The mechanical response of an assembly of particles depends on the applied boundary conditions. Robust calibration of numerical discrete systems to laboratory results is also a primary step in many studies of granular materials. In this study, a new membrane model was developed for simulating axisymmetric element tests. This membrane model uses a simple algorithm of an array of independently controlled walls and is computationally efficient. The effect of boundary flexibility on the system response was investigated by simulating a series of triaxial tests on dense and loose specimens. At the specimen scale, differences in shear strength and volume change of specimens were observed. It was shown that localization pattern depends on the applied boundary conditions. At the particle scale, particle‐membrane contact forces, coordination number, local void ratio, and anisotropy of fabric were all affected by the boundary flexibility.     

Shear wave propagation in granular assemblies

Shear wave velocity is a fundamental property of a granular assembly. It is a measure of the true elastic stiffness of a bulk specimen of discrete grains. Shear wave velocity is typically measured in the laboratory (e.g., using bender elements) or in-situ (e.g., using a seismic cone penetrometer, sCPT). In the current work, shear wave propagation is modeled numerically using the discrete element method (DEM). First, an appropriate method for measuring wave velocity is identified. Then the effects of particle size and elastic properties are investigated. Specimen fabric is then quantified before and after wave excitation and the elasticity of the response at the scale of the particle contacts is investigated. The results show that shear wave velocity may be robustly measured for discrete numerical specimens. The ability to measure shear wave velocity using DEM simulations may provide another tool for researchers seeking to link results from physical and numerical experiments.     

基于双胞元的颗粒材料应力-应变关系研究

基于细观力学,建立颗粒材料的宏观应力-应变与接触力、接触位移、枝矢量等细观量之间的关系。用改进的Voronoi-Delaunay法对颗粒材料进行几何和物理上划分,得到改进Bagi双胞元体系;以固体胞元为基础,运用牛顿第二定律和Gauss定理提出含有旋转矢量和重力的颗粒材料平均等效应力,避免了颗粒材料的准静态假设;在孔隙胞元区域内利用变形协调条件推导出含有孔隙面矢量等几何变量的颗粒材料平均等效应变。结合文献的二维颗粒材料宏观试验结果验证了双胞元平均等效应力-应变的正确性;在三维情形下,对比双胞元等效应变和最优拟合应变结果,同样验证了基于双胞元的颗粒材料应力-应变关系,因此,该颗粒材料应力-应变关系可以为数值模拟颗粒材料力学行为提供依据。     

抗转动对颗粒材料组构特性的影响研究

颗粒材料大多由不规则形状的颗粒组成,如砂土、谷物等,抵抗转动是不规则形状颗粒的固有特性。已有研究表明,颗粒抗转动特性对其宏观力学特性有显著影响。因此,在颗粒材料的细观数值模拟中或采用非圆颗粒,或在圆颗粒离散元模拟中采用考虑抗转动的接触模型。采用有限元-离散元耦合方法（FDEM）和离散元方法（DEM）分别对椭球形状颗粒和具有抗转动能力的圆球颗粒进行三轴剪切数值模拟,指出了采用抗转动接触模型考虑颗粒形状影响的局限性,并基于颗粒的局部排布结构揭示了形状影响的细观来源。峰值内摩擦角和剪胀均随着转动摩擦系数和形状偏离圆球程度而单调增加,但颗粒形状对它们的影响呈现出明显的收敛趋势。细观组构分析也表明,虽然颗粒形状和转动摩擦都能显著增强组构各向异性,但是组构各向异性的演化模式有明显的区别。造成以上结果的差异在于其抵抗转动的影响机制不同。转动摩擦是通过限制颗粒转动,增强了颗粒间的稳定承载能力,而非圆颗粒是通过咬合作用形成稳定的局部排列结构。由于椭球颗粒腹部比端部能够传递更大的接触力,颗粒受剪切后发生转动,颗粒长轴倾向于正交大主应力方向,呈现交错排列,颗粒间相互锁定。     

控制应力路径散体材料真三轴试验强度特征的离散元模拟与细观机制分析

通过一系列真三轴离散元数值试验,模拟了不同应力路径下的等b试验中散体材料的强度特征。根据模拟结果详细地分析了三维应力条件下中主应力和应力路径对散体材料峰值强度的影响,研究了峰值摩擦角、峰值应力比的变化规律,并根据真应力的概念和组构张量的演化结果分析了散体材料的强度成因。研究表明,在不同类型的数值试验中峰值偏应力随b参数的变化规律不同,但采用初始围压归一化后的应力-应变曲线规律一致。峰值强度线的斜率只与b值有关而与应力路径无关,且随着b值的增加,峰值应力比qf /pf逐渐减小,数值模拟结果与室内试验结果吻合较好;随着应变的发展,数值试样的组构也随之发生变化,产生了明显的应力诱发各向异性;散体的强度为颗粒摩擦及材料各向异性共同作用的结果;理论上,组构比-应力比坐标系中破坏点位置仅取决于颗粒摩擦角 ,而数值模拟结果与理论值的差异源于颗粒间咬合和滚动摩擦的影响,其影响与颗粒表面摩擦系数有关,也受空间应力状态的影响。     

A micro–macro investigation of the capillary strengthening effect in wet granular materials

Wet granular materials are three-dimensionally simulated by the discrete element method with water bridges incorporated between particles. The water bridges are simplified as toroidal shapes, and the matric suction is constantly maintained in the material. A comparison with experimental tests in the literature indicates that the toroidal shape approximation may be one of the best choices with high practicability and decent accuracy. Mechanical behaviours of wet granular materials are studied by triaxial tests. Effects of particle size distributions and void ratios are investigated systematically in this study. The hydraulic limit of the pendular state is also discussed. It gives the capillary cohesion function which is not only determined by the degree of saturation but also positively correlated to relative density and particle size polydispersity and inversely proportional to mean particle size. Furthermore, the capillary strengthening effect is also analysed microscopically in aid of the Stress–Force–Fabric relationship, mainly in fabric anisotropy, coordination number and stress transmission pattern, which revealed the micro-mechanisms of the additional effective stress induced by capillary effect.     

川藏铁路线V形深切河谷地形地震放大效应数值模拟

V形河谷场地能产生地震动放大效应,并加剧地震的破坏作用和灾害效应。基于高性能离散元软件MatDEM,对规则V形河谷以及真实河谷地形下的放大效应进行数值模拟研究。结果表明,对于规则V形河谷,在地震波水平入射的情况下,背波坡面产生明显的地震动加速度放大效应,并在临近谷底处最为强烈;随着河谷坡角的增大,放大效应增大。数值模拟结果和解析解具有一致的趋势和规律,并符合现场观测现象。将模型进一步应用于川藏铁路线日扎深切河谷的分析,发现地形微小的转折将会影响到波的传播,临近谷底处的加速度放大效应与规则地形较为一致。离散元法能有效地模拟河谷中地震波的传播、反射、散射等现象,可用于川藏铁路沿线复杂条件下的动力作用灾害效应分析。     

颗粒材料剪胀性的微观力学分析

剪胀性是颗粒材料在加载过程中表现出来的重要变形特性。以孔隙胞元描述颗粒材料内部结构的最小单元,通过对单个孔隙胞元进行剪切受力分析,探讨了剪切过程中颗粒材料体积的改变对应力比和单个孔隙胞元形状的依赖关系,解释了排列密实的颗粒材料在剪切过程中先压缩后剪胀的微观机制。用离散元数值模拟得到了在双轴剪切过程中单个孔隙胞元形状以及孔隙胞元体积变形的演化过程。离散元数值结果表明,加载过程中孔隙胞元形状由初始各向同性到沿大主应力方向变大变长、体积变形先压缩后膨胀,并且体积变形在加载过程中存在局部化现象,体积变化大的孔隙胞元在较大变形时,排列成倾斜的窄带。综合孔隙胞元的受力分析和离散元数值结果表明,致密排列颗粒材料的剪胀性与微观尺度上孔隙胞元的几何结构及其内部的力链传递方式密切相关。     

A computational procedure to assess the distribution of constriction sizes for an assembly of spheres

The effectiveness of filters to counteract internal erosion in earth structures is particularly related to their ability to capture fine particles moving under seepage flow through the porous material. More precisely, fine particles are likely to be trapped by the narrowest paths between pores: the constrictions. This paper proposes a methodology to compute the constriction size distribution of model granular filters taking into account the relative density of the material. The approach is based upon probabilistic methods which adopt stated simple geometric packing arrangements to represent the solid structure in the extreme density states. Two new models are proposed for the design of the constriction size distribution according to the type of filter grading: continuously graded or gap-graded materials. The models require the usual material characteristics: the grading curve, and the minimum and maximum void ratios for this material. Calibrated on the basis of statistical analyses over numerical assemblies of spheres generated by a discrete element method, the proposed new models constitute a promising tool to significantly improve the modeling of filtration processes in granular materials.     

Strength anisotropy of granular material consisting of perfectly round particles

The strength anisotropy of granular materials deposited under gravity has mostly been attributed to elongated particles' tendency to align long axes along the bedding plane direction. However, recent experiments on near‐spherical glass beads, for which preferred particle alignment is inapplicable, have exhibited surprisingly strong strength anisotropy. This study tests the hypothesis that certain amount of fabric anisotropy caused by the anisotropic stress during deposition under gravity can be locked in a circular‐particle deposit. Such locked‐in fabric anisotropy can withstand isotropic consolidation and leads to significant strength anisotropy. 2D discrete element method simulations of direct shear tests on circular‐particle deposits are conducted in this study, allowing for the monitoring of both stress and fabric. Simulations on both monodispersed and polydispersed circular‐particle samples generated under downward gravitational acceleration exhibit clear anisotropy in shear strength, thereby proving the hypothesis. When using contact normal‐based and void‐based fabric tensors to quantify fabric anisotropy in the material, we find that the intensity of anisotropy is discernible but low prior to shearing and is dependent on the consolidation process and the dispersity of the sample. The fact that samples with very low anisotropy intensity measurements still exhibit fairly strong strength anisotropy suggests that current typical contact normal‐based and void‐based second‐order fabric tensor formulations may not be very effective in reflecting the anisotropic peak shear strength of granular materials. Copyright © 2017 John Wiley & Sons, Ltd.     

Homogenization of granular material modeled by a three-dimensional discrete element method

A homogenization strategy for granular materials is presented and applied to a three-dimensional discrete element method (DEM), that uses superellipsoids as particles. Macroscopic quantities are derived from the microscopic quantities resulting from a DEM simulation by averaging over representative volume elements (RVEs). The implementation of an RVE is described in detail regarding the definition and discretization of the RVE boundary. The homogenization strategy is validated by DEM simulations of compression and shear tests of cohesionless granular assemblies. Finally, an elasto-plastic material model is fitted to the resulting stress–strain curves.     

基于耦合方法的挡土墙地震响应的数值模拟

为有效描述动力灾变过程中挡土墙附近土体的细观特性,并大量减小离散元模拟的颗粒数目,节省计算时间,利用离散-连续耦合动力数值模型,模拟了已有离心机试验中重力式挡土墙的地震响应。靠近挡土墙的土体区域用相互作用的离散颗粒模拟,远离挡土墙的土体采用连续模型模拟,编写结构的动力有限元程序并嵌入离散元软件中来模拟挡土墙。通过离散元软件PFC2D和有限差分软件FLAC2D的交互运算实现耦合过程。耦合方法的核心在于:①在离散-连续土体的交界面保证连续性;②模拟中离散区域土体的宏观性质与连续土体模型一致。为进一步满足以上两点,分别提出了新的边界耦合力提取方法和自振柱模拟方法。研究表明耦合方法可以从细观尺度上有效描述土体与结构的相互作用和关心区域的土体特性。     

Analysis of wave propagation in dry granular soils using DEM simulations

In this paper, a three-dimensional particle-based technique utilizing the discrete element method (DEM) is proposed to study wave propagation in a dry granular soil column. Computational simulations were conducted to investigate the soil response to sinusoidal motions with different amplitudes and frequencies. Three types of soil deposits with different void ratios were employed in these simulations. Different boundary conditions at the base such as rigid bedrock, elastic bedrock, and infinite medium were also considered. Analysis is done in time domain while taking into account the effects of soil nonlinear behavior. The computational approach is able to capture a number of essential characteristics of wave propagation including motion amplification and resonance. Dynamic soil properties were then extracted from conducted simulations and used to predict the response of the soil using the widely used equivalent linear method program SHAKE and compare its predictions to DEM results. Generally, there was a good agreement between SHAKE and DEM results except when the exciting frequency was close to the resonance frequency of the deposit where significant discrepancy in computed shear strains between SHAKE predictions and DEM results was observed.     

Discrete element modeling of direct shear tests for a granular material

A succinct 3D discrete element model, with clumps to resemble the real shapes of granular materials, is developed. The quaternion method is introduced to transform the motion and force of a clump between local and global coordinates. The Hertz–Mindlin elastic contact force model, incorporated with the nonlinear normal viscous force and the Mohr–Coulomb friction law, is used to describe the interactions between particles. The proposed discrete element model is used to simulate direct shear tests of the irregular limestone rubbles. The simulation results of vertical displacements and shear stresses with a mixture of clumps are compared well with that of laboratory tests. The bulk friction coefficients are calculated and discussed under different contact friction coefficients and normal stresses. Copyright © 2009 John Wiley & Sons, Ltd.     

A critical state sand plasticity model accounting for fabric evolution

Fabric and its evolution need to be fully considered for effective modeling of the anisotropic behavior of cohesionless granular sand. In this study, a three‐dimensional anisotropic model for granular material is proposed based on the anisotropic critical state theory recently proposed by Li & Dafalias [2012], in which the role of fabric evolution is highlighted. An explicit expression for the yield function is proposed in terms of the invariants and joint invariants of the normalized deviatoric stress ratio tensor and the deviatoric fabric tensor. A void‐based fabric tensor that characterizes the average void size and its orientation of a granular assembly is employed in the model. Upon plastic loading, the material fabric is assumed to evolve continuously with its principal direction tending steadily towards the loading direction. A fabric evolution law is proposed to describe this behavior. With these considerations, a non‐coaxial flow rule is naturally obtained. The model is shown to be capable of characterizing the complex anisotropic behavior of granular materials under monotonic loading conditions and meanwhile retains a relatively simple formulation for numerical implementation. The model predictions of typical behavior of both Toyoura sand and Fraser River sand compare well with experimental data. Copyright © 2013 John Wiley & Sons, Ltd.     

A numerical investigation of the incremental behavior of crushable granular soils

The mechanical behavior of granular materials is characterized by strong nonlinearity and irreversibility. These properties have been differently described by a variety of constitutive models. To test any constitutive model, experimental data relative to the nature of the incremental stress–strain response of the material is desirable. However, this type of laboratory data is scarce because of being expensive and difficult to obtain. The discrete element method has been used several times as an alternative to obtain incremental responses of granular materials. Crushable grains add one extra source of irreversibility to granular materials. Crushability has been variously incorporated into different constitutive models. Again, it will be helpful to obtain incremental responses of crushable granular materials to test these models, but the experimental difficulties are increased. Making use of a recently introduced crushing model for discrete element simulation, this paper presents a new procedure to obtain incremental responses in discrete analogs of granular crushable materials. The parallel probe approach, previously used for uncrushable discrete analogs, is here extended to account for the presence of crushable grains. The contribution of grain crushing to the incremental irreversible strain is identified and separately measured. Robustness of the proposed method is examined in detail, paying particular attention to aspects such as dynamic instability or crushing localization. The proposed procedure is later applied to map incremental responses of a discrete analog of Fontainebleau sand on the triaxial plane. The effect of stress ratio and granular state on plastic flow characteristics is highlighted. Copyright © 2016 John Wiley & Sons, Ltd.     

DEM study of fabric features governing undrained post-liquefaction shear deformation of sand

In an effort to study undrained post-liquefaction shear deformation of sand, the discrete element method (DEM) is adopted to conduct undrained cyclic biaxial compression simulations on granular assemblies consisting of 2D circular particles. The simulations are able to successfully reproduce the generation and eventual saturation of shear strain through the series of liquefaction states that the material experiences during cyclic loading after the initial liquefaction. DEM simulations with different deviatoric stress amplitudes and initial mean effective stresses on samples with different void ratios and loading histories are carried out to investigate the relationship between various mechanics- or fabric-related variables and post-liquefaction shear strain development. It is found that well-known metrics such as deviatoric stress amplitude, initial mean effective stress, void ratio, contact normal fabric anisotropy intensity, and coordination number, are not adequately correlated to the observed shear strain development and, therefore, could not possibly be used for its prediction. A new fabric entity, namely the Mean Neighboring Particle Distance (MNPD), is introduced to reflect the space arrangement of particles. It is found that the MNPD has an extremely strong and definitive relationship with the post-liquefaction shear strain development, showing MNPD’s potential role as a parameter governing post-liquefaction behavior of sand.