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

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

Structural stability and jamming of self-organized cluster conformations in dense granular materials

We examine emergent, self-organized particle cluster conformations in quasistatically deforming dense granular materials from the perspective of structural stability. A structural mechanics approach is employed, first, to devise a new stability measure for such conformations in equilibrium and, second, to use this measure to explore the evolving stability of jammed states of specific cluster conformations, i.e. particles forming force chains and minimal contact cycles. Knowledge gained on (a) the spatial and temporal evolution of stability of individual jammed conformations and (b) their relative stability levels, offer valuable clues on the rheology and, in particular, self-assembly of granular materials. This study is undertaken using data from assemblies of nonuniformly sized circular particles undergoing 2D deformation in two biaxial compression tests: a discrete element simulation of monotonic loading under constant confining pressure, and cyclic loading of a photoelastic disk assembly under constant volume. Our results suggest that the process of self-assembly in these systems is realized at multiple length scales, and that jammed force chains and minimal cycles form the basic building blocks of this process. In particular, 3-cycles are stabilizing agents that act as granular trusses to the load-bearing force chain columns. This co-evolutionary synergy between force chains and 3-cycles proved common to the different materials under different loading conditions. Indeed, the remarkable similarities in the evolution of stability, prevalence and persistence of minimal cycles and force chains in these systems suggest that these structures and their co-evolution together form a generic feature of dense granular systems under quasistatic loading.     

颗粒物质中的多尺度问题

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

颗粒材料接触特性之接触点的确定

颗粒物质在自然界中普遍存在,如砂砾、谷物、堆石体和碎屑流等,许多自然现象和工业过程都与颗粒物质的运动状态有关。弄清颗粒之间的接触特征是理解颗粒变形与运动规律的关键,与之相关联的接触点确定是获取全场接触力(包括力链)信息的难点。近年来,光弹法逐渐成为这类问题研究的有效手段之一,原因在于该方法具有原理简单、操作方便、图像直观等优点。然而,已有的接触点判定方法主要依赖于圆盘圆心距离的计算,对于小变形圆盘接触特征难以有效识别。为此,本文发展了一种基于去除零级条纹区域及图像矩阵转换的算法,可以有效地确定接触点的位置。主要步骤如下:首先运用主成分分析方法对图像进行分析,提取主成分,并重构原图像;然后利用模板匹配方法提取各个光弹圆盘,接着对提取出来的各个圆盘进行图像矩阵变换;最后对变换后的图像局部进行垂直方向灰度梯度平方计算,描绘梯度变化曲线并提取曲线峰值,从而识别出接触点(位置)。     

离散元法铁粉末压制中粒径分布对力链演化机制的影响

为阐明粒径分布对铁粉压制中体系内部细观力学行为的影响, 基于离散元理论, 通过改变铁粉颗粒粒径分布建立压制模型, 结合力链提取方法, 通过对力链空间分布、力链数目、力链长度和力链方向性的分析, 探究粒径分布对力链演化的影响机理. 研究结果表明: 不同粒径分布的粉体压制时形成的力链空间分布具有差异, 粒径分布范围越小, 形成的力链分布越集中, 反之, 粒径分布范围越大, 形成的力链分布越松散且均匀; 在粉末压制时, 粒径分布对力链数目也有影响, 具体表现为随着粉体的粒径分布范围变大, 力链总数逐渐减少; 粉体的粒径分布对颗粒形成短力链的数目起着显著影响, 而对力链长度的影响较为有限; 随着粒径分布范围的增大, 力链的方向由均匀分布逐渐集中在特定角度方向, 表现出一定各向异性, 形成的交叉力链网络结构有利于提高粉体致密化程度. 本文为从粉体粒径分布影响层面拓展粉末压制细观力学理论提供基础, 亦为进一步结合粉体粒径分布及体系内力链演变过程改善粉末致密化行为提供指导.      

High-amplitude elastic solitary wave propagation in 1-D granular chains with preconditioned beads: Experiments and theoretical analysis

Elastic solitary waves resulting from Hertzian contact in one-dimensional (1-D) granular chains have demonstrated promising properties for wave tailoring such as amplitude-dependent wave speed and acoustic band gap zones. However, as load increases, plasticity or other material nonlinearities significantly affect the contact behavior between particles and hence alter the elastic solitary wave formation. This restricts the possible exploitation of solitary wave properties to relatively low load levels (up to a few hundred Newtons). In this work, a method, which we term preconditioning, based on contact pre-yielding is implemented to increase the contact force elastic limit of metallic beads in contact and consequently enhance the ability of 1-D granular chains to sustain high-amplitude elastic solitary waves. Theoretical analyses of single particle deformation and of wave propagation in a 1-D chain under different preconditioning levels are presented, while a complementary experimental setup was developed to demonstrate such behavior in practice. The experimental results show that 1-D granular chains with preconditioned beads can sustain high amplitude (up to several kN peak force) solitary waves. The solitary wave speed is affected by both the wave amplitude and the preconditioning level, while the wave spatial wavelength is still close to 5 times the preconditioned bead size. Comparison between the theoretical and experimental results shows that the current theory can capture the effect of preconditioning level on the solitary wave speed.     

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

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

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

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

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

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

Dynamic stress chain formation in a two-dimensional particle bed

Cure cast plastic bonded explosives (PBXs) consist of relatively hard particles in a soft binder. Under compressive loading, the explosive cyrstals come into contact that causes high stress concentrations. The lines along which the crystals are loaded are called stress chains. Damage done to these particle beds during compressive loading can lead to reaction. The photoelastic effect of PMMA is exploited to examine the stress state within a two-dimensional particle bed. Stress chain development within the bed is recorded and is shown to increase the stress state within some particles while leaving others unloaded. These concentrations form early in the loading process, leading to fracture along the stress bridges and generating likely reaction initiation sites. Through material point method simulations, contact friction is shown to have a large effect on the stress distribution within the particle bed.     

Superquadric DEM-SPH coupling method for interaction between non-spherical granular materials and fluids

The research on the coupling method of non-spherical granular materials and fluids aims to predict the particle–fluid interaction in this study. A coupling method based on superquadric elements is developed to describe the interaction between non-spherical solid particles and fluids. The discrete element method (DEM) and the smoothed particle hydrodynamics (SPH) are adopted to simulate granular materials and fluids. The repulsive force model is adopted to calculate the coupling force and then a contact detection method is established for the interaction between the superquadric element and the fluid particle. The contact detection method captures the shape of superquadric element and calculates the distance from the fluid particle to the surface of superquadric element. Simulation cases focusing on the coupling force model, energy transfer, and large-scale calculations have been implemented to verify the validity of the proposed coupling method. The coupling force model accurately represents the water entry process of a spherical solid particle, and reasonably reflects the difference of solid particles with different shapes. In the water entry process of multiple solid particles, the total energy of the water entry process of multiple solid particles tends to be stable. The collapse process of the partially submerged granular column is simulated and analyzed under different parameters. Therefore, this coupling method is suitable to simulate fluid–particle systems containing solid particles with multiple shapes.     

EFFECT OF INTERGRANULAR FRICTION ON THE BOTTOM FORCE CHAIN OF TWO-DIMENSIONAL PARTICLE PILES WITH POINT DEFECTS1)

力链对颗粒物质的宏观与微观力学性质起着决定性的作用。在离散元法的基础上,建立二维规则排列的颗粒物质系统,分别研究无缺陷颗粒系统在集中载荷变化与有缺陷颗粒系统在缺陷区域改变时,粒间摩擦系数对颗粒系统底部接触力分布规律的影响。结果表明：在无缺陷颗粒系统中,颗粒系统底部接触力的分布形式受摩擦系数和集中载荷的大小影响。随摩擦系数的增大,底部接触力由双峰形式经平台过渡,逐渐向单峰形式转变。在有缺陷颗粒系统中,摩擦系数和缺陷尺寸对底部接触力分布均有影响。同种载荷作用下,随缺陷尺寸的增大,底部接触力峰值显著增大;底部平均力被明显削弱,力向边界的传递增强。系统中轴线上缺陷的存在使底部中间区域受力削弱,当缺陷尺寸超过2 层以上时,底部中间力随摩擦系数的变化特征由递增曲线演变为线形衰减曲线。     

基于持续同调的边坡土颗粒应力链失稳分析

采用颗粒离散元方法和持续同调理论,研究了内排土场堆叠至不同高度时的边坡稳定性。为便于研究,现采用一水平金属板向下施加压力,代替不同厚度土层的重力荷载,对边坡在竖向荷载作用下的失稳破坏过程进行了颗粒离散元模拟。研究了二维边坡土颗粒速度总矢量、边坡失稳破坏时滑移开裂面的角度以及边坡坡顶y方向的平均速度等宏观响应过程,并构建了自然堆积下边坡堆积体颗粒的法向力链无向网络模型。最后,用持续同调方法对边坡坡顶颗粒接触力链网络的拓扑特征进行分析,获得条码图,建立了岩体结构持续同调特征与失稳演化的关系。本文为研究边坡失稳拓扑识别提供了一种新方法,从而可以有效预测边坡的失稳破坏。     

Dynamic Inter-Particle Force Inference in Granular Materials: Method and Application

Inter-particle force transmission in granular media plays an important role in the macroscopic static and dynamic behavior of these materials. This paper presents a method for inferring inter-particle forces in opaque granular materials during dynamic experiments. By linking experimental measurements of particle kinematics and volume-averaged strains to forces, the method provides a new tool for quantitatively studying force transmission and its relation to macroscopic behavior. We provide an experimental validation of the method, comparing inter-particle forces measured in a simple impact test on two-dimensional rubber grains to a finite-element simulation. We also provide an application of the method, using it to study inter-particle forces during impact of an intruder on a granular bed. We discuss the current challenges for applying the method to both model materials and real geologic materials.     

Evolution of mesoscopic granular clusters in comminution systems: a structural mechanics model of grain breakage and force chain buckling

A major scientific challenge in establishing a micromechanics theory for complex materials is the characterisation and modelling of emergent mesoscopic phenomena. This study demonstrates the key elements of a structural mechanics approach to the modelling of mesoscopic dissipative phenomena in comminution systems where grain breakage and force chain buckling coexist. Given the many degrees of freedom in these systems, there are multitude of possible configurations and configurational transitions accessible even for a small particle cluster (e.g. a particle and its immediate neighbours). Here, we develop a model of the evolution of a 6-particle cluster undergoing breakage and force chain buckling, in sequence. The analysis lays bare the intricate connections between the contact topology, the relative kinematics arising from the interactions of particles at the bonded versus non-bonded contacts, and the collective dynamics of these interactions as the cluster is monotonically compressed under confinement. The stress-displacement response profiles at the cluster scale exhibit qualitatively similar properties to those seen in macroscopic assemblies under confined compression. A parametric analysis is undertaken to explore the effects of grain-scale resistances to breakage and buckling with respect to the overall force-displacement behaviour of the granular cluster. The study casts light on open problems for future research into the micromechanics of emergent cluster behaviour germane to comminution systems.     

Nonlinear dynamics of excited plate immersed in granular matter

The interaction between granular matter and the elastic body is a complex issue due to the complex properties of granular matter. An experiment involving a sinusoidally excited plate buried in glass bead particles contained in a box is conducted. The motion behavior of the plate is observed and recorded by the strain gauge. The amplitude–frequency and phase–frequency curves are recorded to study the natural property of the plate in granular matter. In this experiment, jump phenomena are found in both the amplitude–frequency and phase–frequency planes in circumstances with smaller particle sizes, lower buried depths, and larger amplitudes of the excitation force. Otherwise, the period-doubling bifurcation, especially 3T, is found with the increase in the excitation force. These bifurcations usually occur in specific buried depth and excitation frequency band and require smaller particle sizes. The experiments with random-shaped particles exhibit no-jump phenomenon, but period-doubling bifurcation and chaos. These phenomena are sensitive to parameters and closely related to the varying process of the excitation frequency and force. Reasonable mechanisms are summarized qualitatively through some of our recent researches in this paper.     

冲击载荷作用下颗粒材料动态力学响应的近场动力学模拟

颗粒材料在冲击载荷作用下的动态力学行为是学术界关注的热点问题. 新近问世的近场动力学(peridynamics)理论将材料视为由大量有限体积和有限质量的物质点组成,基于非连续性和非局部作用假定建模,建立空间积分形式的运动方程,自然适应于颗粒材料动态力学行为的描述与分析. 发展了描述颗粒间接触作用的物质点尺度的排斥力模型,考虑近场动力学方法中非局部长程力特征,改进了近场动力学中的初始微观弹脆性(prototype microelastic brittle, PMB) 模型的本构力函数,并消除了原PMB 模型中存在的“边界效应” 问题. 计算分析了冲击载荷作用下碳化钨陶瓷颗粒体系的动态力学响应,得到了不同冲击速度下颗粒体系的冲击波速,PD计算结果与试验结果高度一致;通过颗粒物质点尺度作用描述单颗粒尺度的接触作用,很好地再现了颗粒的转动与平动、颗粒挤压变形以及颗粒破碎等现象;刚性冲击板附近同时存在严重的颗粒破碎与轻微的颗粒损伤,远离冲击板的部分颗粒出现破损,且颗粒破碎主要是由颗粒间挤压、碰撞以及相对滑动剪切作用造成的. 研究结果表明,所发展的计算模型和分析方法能很好地反映颗粒材料动态力学行为,具有广泛的应用价值.      

Influence of aspect ratio and arrangement direction on the shear behavior of ellipsoids

This work studies the macroscopic and microscopic behaviors of ellipsoids under triaxial tests using 3D discrete element method (DEM) simulation. To avoid the boundary effect, a novel stress servo-controlled periodic boundary condition is proposed to maintain the confining pressure of samples during testing. The shape features of ellipsoids are investigated, including the aspect ratio of elongated/oblate ellipsoids and the initial arrangement directions of ellipsoids. The macroscopic properties of ellipsoidal particle samples, such as the deviatoric stress, volumetric strain, internal friction angle, as well as dilatancy angles are explored. Elongated and oblate ellipsoids with varying aspect ratios are investigated for the occurrence of stick-slips. In addition, it is demonstrated that the initial arrangement direction has a significant impact on the coordination number and contact force chains. The corresponding anisotropy coefficients of the entire contact network are analyzed to probe the microscopic roots of macroscopic behavior.     

Force Chain Characteristics and Effects of a Dense Granular Flow System in a Third Body Interface During the Shear Dilatancy Process

In order to investigate the characteristics of force chains in a granular flow system, a parallel plate shear cell is constructed to simulate the shear movement of an infinite parallel plate and observe variations in relevant parameters. The shear dilatancy process is divided into three stages, namely, plastic strain, macroscopic failure, and granular recombination. The stickslip phenomenon is highly connected with the evolution of force chains during the shear dilatancy process. The load–distribution rate curves and patterns of the force chains are utilized to describe the load-carrying behaviors and morphologic changes of force chains separately. Force chains, namely, “diagonal gridding,” “tadpole-shaped,” and “pinnate” are defined according to the form of the force chains in the corresponding three stages.