Numerical simulation study on particle breakage behavior of granular materials in confined compression tests

In order to study the fragmentation law, the confined compression experiment of granular assemblies has been conducted to explore the particle breakage characteristic by DEM approach in this work. It is shown that contact and contact force during the loading process gradually transform from anisotropy to isotropy. Meanwhile, two particle failure modes caused by different contact force states are analyzed, which are single-through-crack failure and multi-short-crack failure. Considering the vertical distribution of the number of cracks and the four characteristic stress distributions (the stress related to the maximum contact force, the major principal stress, the deviatoric stress and the mean stress), it is pointed out that the stress based on the maximum contact force and the major principal stress can reflect the distribution of cracks accurately. In addition, the size effect of particle crushing indicates that small size particles are prone to break. The lateral pressure coefficient of four size particles during the loading process is analyzed to explain the reason for the size effect of particle breakage.     

Dilation and breakage dissipation of granular soils subjected to monotonic loading

Dilation and breakage energy dissipation of four different granular soils are investigated by using an energy balance equation. Due to particle breakage, the dilation curve does not necessarily pass through the origin of coordinates. Breakage energy dissipation is found to increase significantly at the initial loading stage and then gradually become sta-bilised. The incremental dissipation ratio between breakage energy and plastic work exhibits almost independence of the confining pressure. Accordingly, a plastic flow rule consid-ering the effect of particle breakage is suggested. The critical state friction angle is found to be a combination of the basic friction between particles and the friction contributed by par-ticle breakage.     

可破碎颗粒体在动力载荷下的耗能特性

采用离散元的数值方法, 通过连接键将若干小颗粒绑定为一个具有不规则外形的大颗粒体, 设置不同连接键强度模拟了颗粒体在外加动力载荷下破碎过程, 并探讨其中系统能量耗散特性. 计算结果表明, 颗粒体的破碎程度决定了系统能量耗散率, 即内部耗能占外界输入能量的比例. 破碎率越高, 颗粒间相互摩擦和碰撞越剧烈,系统能量耗散率越高. 同时, 在循环载荷下系统内颗粒体破碎绝大部分发生在加载初期, 随着颗粒体的分解破碎速率逐渐减小, 系统耗能能力也随之降低.      

A coupled sliding and rolling friction model for DEM calibration

The accuracy of dense Discrete Element Method (DEM) simulations is sensitive to initial density, contact orientation, particle size and shape, and interparticle interaction parameters including contact stiffness, cohesion, coefficients of friction, and coefficients of restitution. Although studies have characterized the effects of individual particle interaction parameters on mechanical responses of loaded granular material, research combining DEM parameters for calibration is scarce. Robust DEM calibration methodology combining sliding and rolling friction coefficients was developed and validated to predict bulk residual soil strength of initially dense DEM particle assemblies.     

基于水平集-离散元方法的球谐函数颗粒材料缓冲性能分析

在自然环境与工业领域中,颗粒材料是一种常见的缓冲材料,其中大量形态各异的非球形颗粒表现出复杂的力学特性并应用于不同工程领域。本文采用球谐函数构造不同球面度和表面凹凸特性的非规则颗粒,通过水平集方法计算球谐函数颗粒间的接触点和碰撞力,并对冲击过程中球形和凹形颗粒的缓冲性能进行离散元分析。数值结果表明,颗粒床厚度、冲击速度和颗粒形状显著影响球谐函数颗粒材料的缓冲性能。颗粒床底部的冲击力峰值随着颗粒床厚度和表面凹凸性的增加而降低,同时冲击力峰值随着冲击速度和颗粒球面度的增加而增加。与球形颗粒相比,球谐函数颗粒具有凹凸表面和多接触点特性,这有利于冲击荷载向四周扩展并提高凹形颗粒的缓冲效果。     

A nonlocal contact formulation for confined granular systems

We present a nonlocal formulation of contact mechanics that accounts for the interplay of deformations due to multiple contact forces acting on a single particle. The analytical formulation considers the effects of nonlocal mesoscopic deformations characteristic of confined granular systems and, therefore, removes the classical restriction of independent contacts. This is in sharp contrast to traditional contact mechanics theories, which are strictly local and assume that contacts are independent regardless the confinement of the particles. For definiteness, we restrict attention to elastic spheres in the absence of gravitational forces, adhesion or friction. Hence, a notable feature of the nonlocal formulation is that, when nonlocal effects are neglected, it reduces to Hertz theory. Furthermore, we show that, under the preceding assumptions and up to moderate macroscopic deformations, the predictions of the nonlocal contact formulation are in remarkable agreement with detailed finite-element simulations and experimental observations, and in large disagreement with Hertz theory predictions—supporting that the assumption of independent contacts only holds for small deformations. The discrepancy between the extended theory presented in this work and Hertz theory is borne out by studying periodic homogeneous systems and disordered heterogeneous systems.     

Effect of particle size distribution on micro- and macromechanical response of granular packings under compression

The role of particle size heterogeneity on micro- and macromechanical properties of assemblies of spherical particles was studied using DEM simulations. The response to an imposed load of a granular material composed of non-uniformly sized spheres subjected to uniaxial confined compression was investigated. A range of geometrical and micro-mechanical properties of granular packings (e.g., void fraction, contact force distribution, average coordination number and degree of mobilisation of friction at contacts between particles) were examined, and provided a more accurate interpretation of the macroscopic behaviour of mixtures than has previously been available. The macromechanical study included stress transmission, stiffness and angle of internal friction of the granular assemblies.The degree of polydispersity showed slight effect on both, the void fraction and the elastic properties of the system. The tendency for increase in the lateral-to-vertical pressure ratios was observed with an increasing degree of particle size heterogeneity; however, the different pressure ratios calculated for samples with various degrees of polydispersity lay within the range of data scatter.     

基于颗粒物质力学的铁粉末压制中摩擦特性对力链演化影响

结合颗粒物质力学理论,通过离散元法实现铁粉末压制过程模拟并通过压制方程进行验证,针对粉末体系中的力链演化问题,提出力链特征定量分析方式,进一步通过分析不同颗粒间摩擦系数、侧壁摩擦系数与颗粒运动状态转变的方式,探讨摩擦特性对力链量化特征的影响,从而建立摩擦行为与力链演化间的联系. 研究结果表明:随颗粒间摩擦系数增大,整体力链数目变少,力链方向系数、承载不均匀度及单位屈曲度均变大,而随侧壁摩擦系数增大,力链特征差异较小,则颗粒间摩擦系数较侧壁摩擦系数对力链特征演化具有更显著影响. 同时发现,颗粒接触状态的改变与力链特征演化间具有对应性. 研究成果将进一步拓展粉末压制中考虑摩擦行为及力链演化过程在内的粉体致密化行为理论.      

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.     

Dynamic coupling of fluid and structural mechanics for simulating particle motion and interaction in high speed compressible gas particle laden flow

In this work, structural finite element analyses of particles moving and interacting within high speed compressible flow are directly coupled to computational fluid dynamics and heat transfer analyses to provide more detailed and improved simulations of particle laden flow under these operating conditions. For a given solid material model, stresses and displacements throughout the solid body are determined with the particle–particle contact following an element to element local spring force model and local fluid induced forces directly calculated from the finite volume flow solution. Plasticity and particle deformation common in such a flow regime can be incorporated in a more rigorous manner than typical discrete element models where structural conditions are not directly modeled. Using the developed techniques, simulations of normal collisions between two 1 mm radius particles with initial particle velocities of 50–150 m/s are conducted with different levels of pressure driven gas flow moving normal to the initial particle motion for elastic and elastic–plastic with strain hardening based solid material models. In this manner, the relationships between the collision velocity, the material behavior models, and the fluid flow and the particle motion and deformation can be investigated. The elastic–plastic material behavior results in post collision velocities 16–50% of their pre-collision values while the elastic-based particle collisions nearly regained their initial velocity upon rebound. The elastic–plastic material models produce contact forces less than half of those for elastic collisions, longer contact times, and greater particle deformation. Fluid flow forces affect the particle motion even at high collision speeds regardless of the solid material behavior model. With the elastic models, the collision force varied little with the strength of the gas flow driver. For the elastic–plastic models, the larger particle deformation and the resulting increasingly asymmetric loading lead to growing differences in the collision force magnitudes and directions as the gas flow strength increased. The coupled finite volume flow and finite element structural analyses provide a capability to capture the interdependencies between the interaction of the particles, the particle deformation, the fluid flow and the particle motion.     

Macro–micro relations in granular mechanics

In granular mechanics, macroscopic approaches treat a granular material as an equivalent continuum at macro-scale, and study its constitutive relationship between macro-quantities, such as stresses and strains. On the other hand, microscopic approaches consider a granular material as an assembly of individual particles interacting with each other at micro-scale (i.e., particle-scale), and the physical quantities under study are forces and displacements. This paper aims at linking up the findings from these two scales and to establish the macro–micro relations in granular mechanics.Three aspects of the macro–micro relations are investigated. They are about the internal structure, the stress tensor and the strain tensor. The internal structure is described with geometrical systems at the particle scale. Micro-structural definitions of the stress and strain tensors are derived, which link the macro-stress tensor with the contact forces, and the macro-strain tensor with the relative displacements at contact. In addition to a brief review of the past research work on these topics, further generalizations are made in this paper. In particular, the two cell systems proposed by Li and Li (2009), namely the solid cell system and the void cell system, are introduced and used for the derivation of the macro-structural expressions. The stress expression is derived based on Newton’s second law of motion. The result is valid for both static and dynamic cases. The strain expression is derived based on the compatibility requirement. And the expression is valid for any tessellation subdividing the granular assembly into polyhedral elements.The homogenization for deriving a macroscopic constitutive relationship from microscopic behaviour is discussed. Attention is placed on the macroscopic quantification of the internal structure in terms of a second rank tensor, known as the fabric tensor. Existing definitions of the fabric tensors have been reviewed. The correlations among different fabric tensors and their relations with the stress–strain behaviours have been investigated.     

岩土类颗粒物质宏-细观力学研究进展

岩土类颗粒物质在自然界、工程建设以及日常生活中普遍存在,其运动特性的研究在力学界已经开展了几十年.在近20年开展的一系列小尺寸物理实验中,颗粒物质表现出许多新奇现象,人们从物理角度开展了系统研究,在统计力学中,颗粒固体的流体动力学等理论研究以及实验检测技术等方面都取得突破性进展,深刻地揭示了颗粒材料的物理机制,促使力学...     

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.     

颗粒材料破碎演化路径细观热力学机制

颗粒材料在高应力环境下会发生颗粒破碎现象,颗粒破碎不仅影响颗粒材料的力学特性,同时与大量工程问题密切相关.目前的相关研究主要集中在唯象地描述颗粒破碎的演化以及破碎对力学特性的影响层面,对颗粒破碎演化路径的物理机制研究较少.本文基于热力学框架,采用细观力学中细观-宏观的均匀化方法推导了颗粒体系弹性能和破碎能量耗散,并在最大能量耗散的假设下,在热力学框架内,建立了理想化的无摩擦球体颗粒等向压缩过程的弹性-破碎模型,阐述了颗粒材料破碎演化路径细观热力学机制.由于模型的推导不依赖任何唯象的经验公式,因此模型中包含的参数均有明确的物理意义.模型预测与前人试验结果对比表明,材料的初始级配对弹性压缩模量和破碎应力的影响并不相同:不同分形维数级配对应的弹性体变模量存在极大值,而破碎应力却随着分形维数的增大单调递增;颗粒破碎的演化符合最大能量耗散原理,且颗粒材料的压缩曲线可以分为弹性-破碎-拟弹性3个机制不同的阶段.      

颗粒物质在竖直管中下行蠕动流动的实验测量

颗粒流蠕动行为是颗粒物质在竖直管中流动时常见的一种流动现象,其产生机理较复杂。为此,本文在在内径为150mm、高为5000mm的竖直管实验装置上,以FCC催化剂为固体颗粒物料,采用PV6型颗粒速度测量仪,测量不同颗粒流率下竖直管中的颗粒下行蠕动流动速度以及颗粒固含率,系统地考察了颗粒物质在竖直管中下行流动时的蠕动流动特性及产生机理。实验结果表明,颗粒物质在竖直管中下行流动时的流动行为可划分为两种形式。在颗粒流率较小时,颗粒物质下行速度呈现脉冲式变化,有速度停滞,可称之为蠕动I型流动。随着颗粒流率的增加,颗粒下行速度停顿消失,但仍是起伏变化,为蠕动II型流动。当颗粒流率增大到一定值后,颗粒物质下行蠕动行为消失,转变为流化流动。颗粒物质下行的蠕动行为是出口区颗粒成拱与崩塌、颗粒与器壁滑动摩擦和颗粒力链作用的综合反映。     

颗粒流动力学及其离散模型评述

颗粒流是由众多颗粒组成的具有内在相互作用的非经典介质流动. 自然界常见颗粒流都是密集流, 颗粒间接触形成力链, 诸多力链相互交接构成支撑整个颗粒流重量和外载荷的网络, 其局部构型及强度在外载荷下演化, 是颗粒流摩擦特性和接触应力的来源.本文介绍球形颗粒间无粘连作用时的Hertz法向接触理论和Mindlin-Deresiewicz切向接触理论. Campbell依据是否生成较为稳定的力链把颗粒流分为弹性流和惯性流两大类, 其中弹性-准静态流和惯性-碰撞流分别对应准静态流和快速流, 作为两种极端流动情况通常处理成连续体, 分别采用摩擦塑性模型和动理论予以描述, 但是表征接触力链的颗粒弹性参数并不出现这两个模型和理论框架中, 如何进一步考虑颗粒弹性参数将非常困难. 目前离散动力学方法逐渐成为复现其复杂颗粒流动现象、提取实验不可能获得的内部流动信息进而综合起来探索颗粒流问题的一种有效工具, 其真实性强于连续介质理论的描述. 软球模型对颗粒间接触力简化处理, 忽略了切向接触力对法向接触力及其加载历史的依赖, 带来了法向和切向刚度系数如何标度等更艰难的物理问题, 但由于计算强度小而广泛应用于工程问题中. 硬球模型不考虑颗粒接触变形, 因而不能描述颗粒流内在接触应变等物理机理, 仅适用于快速颗粒流, 这不仅仅是由于两体碰撞的限制. 因此基于颗粒接触力学的离散颗粒动力学模型是崭新的模型,适用于准静态流到快速流整个颗粒流态的模拟, 可以细致考虑接触形变及接触力的细节,建立更为合理的颗粒流本构关系, 进而有力的促进颗粒流这一非经典介质流动的研究.      

TENSILE STRENGTH FOR SPLITTING FAILURE OF BRITTLE PARTICLES WITH CONSIDERATION OF POISSON''''S RATIO

1. Introduction The core mechanism of comminution could be reduced to the breakage of individual particles that occurs through contact with other particles or with the grinding media, or with the solid walls of the mill (Potapov & Campbell, 1994). When a particle is subjected to a load, the nature of the stress field around and within the brittle particle, its material properties, and the size and distribution of micro-flaws within the particle govern the size and shape distribution of the fra…     

Crushing of particles in idealised granular assemblies

Four idealised assemblies of equally sized spherical particles are subjected to a range of macroscopic compressive principal stresses and the contact forces on individual particles are determined. For each set of contact forces the stress fields within individual particles are studied. A failure criterion for brittle materials is imposed and indicates that crushing (or rupture) occurs when the maximum contact force reaches a threshold particle strength value, irrespective of the presence and magnitude of other lesser contact forces acting on the particle and the material properties of the particle. Combining the crushing mechanism with an assembly instability mechanism enables failure surfaces to be drawn in the three-dimensional stress space. A simple spatial averaging technique has been applied to the failure surfaces to remove the effects of assembly anisotropies. Sections of the failure surfaces on π planes have similarities to those commonly used in sand modelling.     

基于离散单元法和人工神经网络的近壁颗粒动力学特征研究

颗粒与壁面的相互作用往往对颗粒流动具有显著影响. 为研究颗粒与壁面作用机理, 对滚筒内颗粒流动过程进行离散单元法(DEM)数值模拟. 基于模拟结果统计分析靠近壁面处颗粒的运动特征, 结果表明, 小摩擦系数时颗粒平动和旋转速度均近似满足正态分布, 但由于壁面影响, 摩擦系数增大时颗粒沿滚筒轴向的旋转速度偏离正态分布, 颗粒动力学理论推导壁面边界条件时应考虑速度正态分布的修正及速度脉动的各向异性. 采用人工神经网络(ANN)构建了颗粒无因次旋转温度、滑移速度和平动温度之间的函数模型, 进而可以在常规双流模型壁面边界条件中考虑颗粒旋转的影响. 基于DEM模拟及结果分析可以为壁面边界条件的理论构造和半经验修正提供基础数据和封闭模型.      

On flows of granular materials

This article reviews the behavior of materials made up of a large assemblage of solid particles under rapid and quasi static deformations. The focus is on flows at relatively high concentrations and for conditions when the interstitial fluid plays an insignificant role. The momentum and energy exchange processes are then primarily governed by interparticle collisions and Coulomb-type frictional contact. We first discuss some physical behavior —dilatancy, internal friction, fluidization and particle segregation — that are typical to the understanding of granular flows. Bagnold's seminal Couette flow experiments and his simple stress analysis are then used to motivate the first constitutive theories that use a microstructural variable — the fluctuation energy or granular temperature — governing the subscale fluctuating motion. The kinetic theories formalize the derivation of the field equations of bulk mass, momentum and energy, and permit derivation of constitutive relations for stress, flux of fluctuation energy and its dissipation rate for simple particle assemblages and when frictional rubbing contact can be ignored. These statistical considerations also show that formulation of boundary conditions needs special attention. The frictional-collisional constitutive behavior in which both Coulomb-type rubbing contact and collisional encounters are significant are discussed. There is as yet no rigorous formulation. We finally present a phenomenological approach that describes rapid flows of granular materials under simultaneous transport of heat and close with a summary of stability analyses of the basic flow down an inclined plane.Dedicated to Professor Dr.-Ing. Franz Gustav Kollmann on the occasion of his sixtieth brithday