The mesoscopic structures of dense granular materials

密集颗粒物质由大量颗粒组成的多体相互作用体系,在一定条件下,颗粒互相连接,形成相对稳定的介观尺度结构,其几何和动力学性质较大程度上决定了颗粒体系的宏观物理和力学性质,因此开展颗粒的介观结构研究具有重要的理论价值,是科学的前沿之一.自然界的堆石坝、堰塞体和碎屑流,以及工程中的高温气冷堆堆芯颗粒流和先进核裂变能系统（ADS嬗变）的颗粒散裂靶等都是典型的颗粒体系,研究颗粒体系宏观力学性质是灾害预测和调控技术的关键.本文首先介绍颗粒接触力理论和简化模型的研究进展,接着介绍介观尺度结构分析方法与测量技术,颗粒体系Jamming转变、软点和颗粒微位移测量技术等,最后列举了几个关键的科学问题.颗粒介质中很多基本力学问题的解决需要借鉴物理和数学等学科的最新成果,建立新的概念和范式,从新的角度、思路、理念去认识颗粒介质的基本问题.同时,颗粒介质的基础研究还要紧密结合工程应用领域的大量相关的核心技术,与工程领域专家共同合作,使得颗粒介质的研究有的放矢,更具生命力.     

基于颗粒物质流体动力学的三轴力学分析

构建了颗粒体系热力学特征函数, 确定了宏观尺度不可逆过程的迁移系数, 将颗粒物质流体动力学理论应用于密砂的三轴加载力学性质分析, 得到了应力-应变关系和体积应变-应变关系, 还得到了颗粒温度演化过程.      

密集颗粒物质的介观结构

密集颗粒物质由大量颗粒组成的多体相互作用体系,在一定条件下,颗粒互相连接,形成相对稳定的介观尺度结构,其几何和动力学性质较大程度上决定了颗粒体系的宏观物理和力学性质,因此开展颗粒的介观结构研究具有重要的理论价值,是科学的前沿之一.自然界的堆石坝、堰塞体和碎屑流,以及工程中的高温气冷堆堆芯颗粒流和先进核裂变能系统(ADS嬗变)的颗粒散裂靶等都是典型的颗粒体系,研究颗粒体系宏观力学性质是灾害预测和调控技术的关键.本文首先介绍颗粒接触力理论和简化模型的研究进展,接着介绍介观尺度结构分析方法与测量技术,颗粒体系Jamming转变、软点和颗粒微位移测量技术等,最后列举了几个关键的科学问题.颗粒介质中很多基本力学问题的解决需要借鉴物理和数学等学科的最新成果,建立新的概念和范式,从新的角度、思路、理念去认识颗粒介质的基本问题.同时,颗粒介质的基础研究还要紧密结合工程应用领域的大量相关的核心技术,与工程领域专家共同合作,使得颗粒介质的研究有的放矢,更具生命力.     

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

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

颗粒材料三维应力路径下的接触组构特性

颗粒材料的宏观应力变形特征与其微观接触力、组构等紧密相关.一般而言,强接触系统属于颗粒内部体系的传力结构,其对应的组构张量是影响宏观应力性质的重要因素.细观数值方法(如离散单元法)能够反映物理试验的基本规律,并且可以方便地提取宏微观数据来研究颗粒体系的应力变形机制.采用离散单元法(discrete element method,DEM)进行一系列等$p$等$b$应力路径下颗粒材料的真三轴试验,在此基础上研究了三维应力路径下颗粒材料的宏微观力学参数的演化过程、三维组构张量与应力张量多重联系以及强接触体系反映的宏观应力特征.研究表明:颗粒体系偏应力峰值状态和临界状态均存在与加载路径无关的宏微观特征;三维应力路径下组构张量与应力张量存在非共轴性,但其联合不变量演化过程表现出加载路径无关的特征;与弱接触系统的组构张量相比,强接触系统的组构张量更能反映宏观应力张量的特征;强弱接触体系的组构张量对颗粒体系宏观响应的贡献不同,其分界点存在一定取值范围,但采用平均接触力较为简单合理.      

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

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

颗粒介质固-流态转变的理论分析及实验研究

颗粒介质由大量离散的颗粒聚集而成,因而与传统固体和流体不同,运动过程中的颗粒介质中可能同时存在多种流态及其相互间复杂的转换过程. 颗粒介质弹性失稳机理、不可恢复应变量化是研究颗粒介质固态和流态及固-流态转变的关键. 在前期建立的双颗粒温度热力学(two-granular-temperature, TGT) 理论基础上,确定了颗粒介质的弹性稳定性条件,建立了不可恢复应变流动法则,搭建了描述颗粒固态-液态及其相互转化的简单模型. 颗粒堆积体坍塌过程是典型的颗粒介质固态和流态及其转变过程,因此本文首先开展了25 167 个陶颗粒堆积体坍塌过程的实验研究,并使用基于TGT 理论的物质点方法和离散元方法对物理实验进行了模拟. 结果表明,模型数值结果与物理实验在颗粒堆坍塌过程中的形态、速度分布等细节上吻合很好,同时也发现了现阶段所使用的物质点方法和TGT 理论的不足. 初步说明TGT 理论可以实现颗粒介质固态和流态,以及状态转变的描述.      

ANALYSIS OF PARAMETERS IN GRANULAR SOLID HYDRODYNAMICS FOR TRIAXIAL COMPRESSION TESTS

Granular materials are omnipresent in industries and in nature. For small strains, elastic-plastic and hypoplastic constitutive relations are widely used in engineering practice, but they are not a significant reflection of the underlying physics. Under a unified thermodynamics framework explaining the physics of materials, granular solid hydrodynamics （GSH） was an ex- tension towards describing granular materials, not only solid-like, but also fluid-like behaviors. In this paper, the fundamentals of GSH are briefly treated and then simplified to analyze quasi- static deformations in triaxial compressions. The calculated stress-strain relations and volumetric strain are compared with experimental results. The influences of the major parameters in GSH, especially their cross coupling influences, are analyzed and their physical meanings are further clarified. After parameters were calibrated, the calculated stress values in the characteristic stress state are found to be within 22% of tested values. Meanwhile, the energy dissipation during triaxial compression is analyzed. The above results support and partially quantify GSH.     

Measurement techniques of grain motion and inter-grain structures in dense granular materials

颗粒材料由大量粗颗粒堆积形成, 是复杂的多体相互作用体系, 呈现出颗粒尺度的结构不均匀和动力学不均匀性的基本特征, 这决定了颗粒材料具有很多独特的宏观性质. 借鉴学科历史的发展途径, 基于统计力学, 从颗粒结构和动力学开始建立颗粒材料体系的宏观连续介质力学理论框架是必然途径.但是, 颗粒材料的基本特征决定了从基本理论到实验手段上, 表征与建立颗粒材料结构与性能的相关性都极其困难.这是由于现有测试分析手段所描述的颗粒系统组织结构过于简单化, 缺乏对颗粒结构和动力学的真正认识, 从而制约了颗粒物质研究的发展.因此, 开展颗粒体系结构和动力学性质的测量, 是理解和认识颗粒材料重要物理和力学问题的基础和依据.笔者来自不同的科研院所, 近十年来开展了颗粒体系结构和动力学性质的测量研究, 主要集中于以下两个方向: (1)数字图像测速法、散斑能见度光谱法和X射线- CT等非侵入式测量技术在颗粒运动方面的应用; (2)体积响应谱、力学谱(有效质量和内耗等)和声速测量技术等直接或间接测量颗粒接触力和颗粒结构技术.本文综述了这些实验手段的基本原理及其特点、取得的主要成果, 以及国际最新进展和困难. 最后是对全文的总结, 结合笔者开展测量的经验和教训, 提出了自己的看法, 并试图展望颗粒材料测量技术研究的前景.     

颗粒物质中的多尺度问题

颗粒物质是大量离散的固体颗粒相互作用而组成的复杂体系. 依据颗粒排布的稀疏程度, 体系可分为颗粒气体、颗粒流体和颗粒固体,它们有不同本质的动量传递和能量耗散机制. 后两者属于密集颗粒物质体系,内部形成了颗粒$\to $力链$ \to$体系的多尺度结 构,并涉及多个特征时间尺度,是典型的多尺度体系. 合理分割体系结构层 次、正确理解不同层次的物理过程、并确定它们之间的关联是密集颗粒物质研究的核心任务. 本文依次分析了密集颗粒物质的内在物理图像、多尺度结构层次和特征时间等,并介绍了多 尺度研究框架.     

散料在料仓内流动特性的实验研究

用振动筛测量多种散料的筛分尺寸，测量散料的容积密度、安息角，用应变直剪仪进行散料内摩擦角、壁面摩擦角的测量，分析了散料颗料的形状、尺寸对散料的安息角、内摩擦角、壁面摩擦角的影响。用示踪法进行了圆筒内散料重力场中颗粒流线和等时线的试验，以及出口直径与流量的关系、筒仓内散料高度对流量的影响的试验，流量与筒仓内散料高度关系的试验。分析了散料各物性参数对散料流动性的影响，为散料存贮和运输设备的结构设计提供实验依据。     

DEM simulation of particle mixing in a sheared granular flow

Mixing behaviors of particles are simulated in a sheared granular flow using differently colored but otherwise identical glass spheres, with five different bottom wall velocities. By DEM simulation, the solid fractions, velocities, velocity fluctuations and granular temperatures are measured.The mixing layer thicknesses are compared with the calculations from a simple diffusion equation using the data of apparent self-diffusion coefficients obtained from the current simulation measurements. The calculations and simulation results showed good agreements, demonstrating that the mixing process of granular materials occurred through the diffusion mechanism.     

DEM simulation of particle mixing in a sheared granular flow

Mixing behaviors of particles are simulated in a sheared granular flow using differently colored but otherwise identical glass spheres, with five different bottom wall velocities. By DEM simulation, the solid fractions, velocities, velocity fluctuations and granular temperatures are measured. The mixing layer thicknesses are compared with the calculations from a simple diffusion equation using the data of apparent self-diffusion coefficients obtained from the current simulation measurements. The calculations and simulation results showed good agreements, demonstrating that the mixing process of granular materials occurred through the diffusion mechanism.     

非均匀颗粒材料的类固-液相变行为及本构方程

以非均匀颗粒介质为研究对象,采用三维离散元方法对其在不同密集度和剪切速率下的动 力过程进行了数值模拟,分析了其在由瞬时接触的快速流动向持续接触的准静态流动的转变 过程及其行为特点. 通过对不同材料性质下相变过渡区内颗粒材料的宏观应力、接触时间数、 配位数、团聚颗粒数量、有效摩擦系数等参量的计算,更加全面地描述了非均匀颗粒材料在 类固-液相变过程中的基本特征. 基于以上数值计算结果,建立了一个适用于颗粒材料 类固态、类液态以及其相变过程的本构方程,并通过剪切室实验结果验证了它的合理性.     

煤仓内煤散料流动状态与力学行为影响因素

针对煤仓内煤散料流动问题及其力学行为,采用三维颗粒流模拟程序PFC3D建立了某型号煤仓与某种煤散料的离散元模型,简述了其力学模型与求解步骤,模拟分析了煤仓内煤散料卸料流动状态。通过分析水平向侧压力、颗粒速度场和接触力场,重点讨论了煤仓下部锥体内壁面摩擦系数、锥仓倾角和卸料口径等对煤散料颗粒流动状态和力学行为的影响。结果显示,深仓卸料流动为整体流动与中心流动混合状态,煤仓内壁摩擦系数、锥体倾角和卸料孔开口半径均对煤散料流动和水平侧压力有较大影响。     

A thermomechanical constitutive model for cemented granular materials with quantifiable internal variables. Part II – Validation and localization analysis

We study the mechanical failure of cemented granular materials (e.g., sandstones) using a constitutive model based on breakage mechanics for grain crushing and damage mechanics for cement fracture. The theoretical aspects of this model are presented in Part I: Tengattini et al. (2014), A thermomechanical constitutive model for cemented granular materials with quantifiable internal variables, Part I – Theory (Journal of the Mechanics and Physics of Solids, 10.1016/j.jmps.2014.05.021). In this Part II we investigate the constitutive and structural responses of cemented granular materials through analyses of Boundary Value Problems (BVPs).The multiple failure mechanisms captured by the proposed model enable the behavior of cemented granular rocks to be well reproduced for a wide range of confining pressures. Furthermore, through comparison of the model predictions and experimental data, the micromechanical basis of the model provides improved understanding of failure mechanisms of cemented granular materials. In particular, we show that grain crushing is the predominant inelastic deformation mechanism under high pressures while cement failure is the relevant mechanism at low pressures. Over an intermediate pressure regime a mixed mode of failure mechanisms is observed. Furthermore, the micromechanical roots of the model allow the effects on localized deformation modes of various initial microstructures to be studied. The results obtained from both the constitutive responses and BVP solutions indicate that the proposed approach and model provide a promising basis for future theoretical studies on cemented granular materials.     

A micro-investigation on water bridge effects for unsaturated granular materials with constant water content by discrete element method

The most common state of surface soil is unsaturated. Changes in water content will substantially impact its strength, leading to geological and engineering catastrophes. This paper used LIGGGHTS software to simulate the water bridge effect of unsaturated granular materials with constant water content and verify the rationality of the simplification of the stress-force-fabric (SFF) relationship. The results showed that the capillary force was not isotropic, which was different from the previous study, thus it cannot be overlooked in the simplification of the SFF relationship. Moreover, the influence of water content on the macroscopic mechanical behavior of unsaturated granular materials was interpreted through the evolutions of coordination number, interparticle force, fabric and force anisotropy, and other microscopic parameters. Compared to the literature, we found that different water bridge models would not change the characteristics of the solid skeleton.     

含液颗粒材料液固耦合分析的离散颗粒模型及特征线SPH法

对含液颗粒材料流固耦合分析建议了一个基于离散颗粒模型与特征线SPH法的显式拉格朗日-欧拉无网格方案。在已有的用以模拟固体颗粒集合体的离散颗粒模型[1]基础上,将颗粒间间隙内的流体模型化为连续介质,对其提出并推导了基于特征线的SPH法。数值例题显示了所建议方案在模拟颗粒材料与间隙流相互作用的能力和性能以及间隙流体对颗粒结构承载能力及变形的影响。     

A thermodynamic model of turbulent motions in a granular material

This paper is devoted to a thermodynamic theory of granular materials subjected to slow frictional as well as rapid flows with strong collisional interactions. The microstructure of the material is taken into account by considering the solid volume fraction as a basic field. This variable is of a kinematic nature and enters the formulation via the balance law of the configurational momentum, including corresponding contributions to the energy balance, as originally proposed by Goodman and Cowin [1], but modified here. Complemented by constitutive equations, the emerging field equations are postulated to be adequate for motions, be they laminar or turbulent, if the resolved length scales are sufficiently small. On large length scales the sub-grid motion may be interpreted as fluctuations, which manifest themselves in correspondingly filtered equations as correlation products, like in the turbulence theory. We apply an ergodic (Reynolds) filter to these equations and thus deduce averaged equations for the mean motions. The averaged equations comprise balances of mass, linear and configurational momenta, energy, and turbulent kinetic energy as well as turbulent configurational kinetic energy. They are complemented by balance laws for two internal fields, the dissipation rates of the turbulent kinetic energy and of the turbulent configurational kinetic energy. We formulate closure relations for the averages of the laminar constitutive quantities and for the correlation terms by using the rules of material and turbulent objectivity, including equipresence. Many versions of the second law of thermodynamics are known in the literature. We follow the Müller-Liu theory and extend Müllers entropy principle to allow the satisfaction of the second law of thermodynamics for both laminar and turbulent motions. Its exploitation, performed in the spirit of the Müller-Liu theory, delivers restrictions on the dependent constitutive quantities (through the Liu equations) and a residual inequality, from which thermodynamic equilibrium properties are deduced. Finally, linear relationships are proposed for the nonequilibrium closure relations.Received: 21 March 2003, Accepted: 1 September 2003, Published online: 11 February 2004PACS: 05.70.Ln, 61.25.Hq, 61.30.-vCorrespondence to: I. Luca     

基于微观结构非均匀性的非晶合金力学行为

非晶合金弛豫/晶化、玻璃转变、塑性变形等热力学和动力学行为都与其固有的结构非均匀性密切相关. 但是, 由于淹没在亚稳的长程无序结构中, 探究非晶合金的结构非均匀性十分困难. 尤其, 非晶合金微观结构非均匀性与其力学性能之间的本征关联是一个亟待解决的重要科学问题. 本文基于多尺度时空下的力学激励阐述非晶合金微观结构非均匀性特征与演化规律. 从实验、理论和数值模拟方面出发, 梳理了非晶合金微观结构非均匀性与弛豫机制和力学行为之间的关联. 最后, 针对非晶合金微观结构非均匀性与其物理/力学性能研究的方向提出了建议和展望.