基于结构两相阻力分析的泥石流堆积形态研究

针对泥石流的结构两相阻力特征,采用结构两相阻力对泥石流的堆积形态进行了静力学分析,推导出在一般条件下泥石流堆积形态的公式,并讨论了在不同条件下泥石流堆积形态的特点。分析表明采用结构两相阻力对泥石流堆积形态分析可以较合理的模化阻力项,能够表现泥石流体复杂的成分结构组成,较真实的反映实际泥石流的堆积形态特征。在特定材料条件下,结构两相阻力项转化为单一材料阻力项,所得到的相应结论与已有对单一材料的分析结论一致。     

基于DEM方法的不同密度颗粒材料砂堆形成过程的数值研究

颗粒堆积是地质工程中较为普遍的现象,而其堆积稳定后的安息角可以在一定程度上反映颗粒堆积的基本特征.用DEM方法研究不同密度分布的颗粒堆积过程及堆积结果,有助于我们在细观层面研究颗粒密度对颗粒堆积结果的影响和验证DEM方法分析颗粒堆积的合理性.本研究用圆形颗粒简化模拟散体颗粒:第一步,用EDEM软件建立相关模型,在模型内生成简化的颗粒,模拟颗粒在重力作用下自由下落和堆积稳定的过程;第二步,模型中生成颗粒下落堆积稳定后,保留模型底部边界,数值模拟研究移除模型上部边界后颗粒散落堆积稳定的过程并测定堆积体的安息角;第三步,对应数值模型,做相关的实验进行对比验证.研究结果表明:在只考虑堆积颗粒密度而忽略其他参数影响的条件下,堆积颗粒材料的密度越小,其堆积稳定后的安息角越大.     

转动驱动圆角立方体颗粒有序堆积的离散元模拟

宏观颗粒体系的有序堆积研究可以为热系统中微观粒子的自组装提供研究模型,也有助于工业生产中提高颗粒材料的填充率.实验发现圆筒中的圆角立方体颗粒受容器往复旋转剪切作用会实现有序堆积.为探究旋转圆筒中立方体颗粒有序堆积过程的内部结构演化过程和动力学机理,采用超二次曲面方程构造了圆角立方体颗粒,基于离散元方法对旋转剪切作用下圆...     

钟摆状态下湿颗粒坍塌流动行为及其动力相似性

基于建立的湿颗粒离散动力学模型,本文系统研究了钟摆状态下湿颗粒柱在重力驱动下的坍塌流动过程,主要考虑了颗粒粒径、液体表面张力系数和液体含量等参数对系统坍塌流动模式和动力学行为的影响。研究发现,在湿颗粒系统中,颗粒粒径和液体表面张力系数会改变颗粒间的毛细力大小,引起系统发生不同的坍塌流动模式,而液体含量仅定量影响颗粒坍塌后的堆积形态。在此基础上,进一步探讨了不同模式下系统坍塌流动行为与模型参数的相关性,发现无量纲Bond数是决定钟摆状态下湿颗粒物质坍塌流动动力学行为的本质因素。     

密集颗粒体系的颗粒运动及结构测量技术

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

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

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

密集颗粒体系的颗粒运动及结构测量技术

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

激光辐照下多层圆柱体中三维瞬态温度场的解析解

考虑外表面的气流影响和层间温度与传热的协调关系,建立了激光辐照下,层合圆柱体中的三维瞬态热传导解析模型。利用特征值法和Bessel函数,导出了各层柱体中三维瞬态温度场的封闭解析解。以一维轴对称问题为例计算了柱体中的瞬态温度场,给出了柱体内部温度随时间的变化和柱体表面换热系数对温度场的影响规律。本文的理论解可进一步用于分析层合圆柱体中的三维瞬态热-力效应,并可作为相应问题的数值模拟中数值模型的修正依据。     

ON INITIATION OF GULLY DEBRIS FLOW BY LIMIT EQUILIBRIUM ANALYSIS1)

沟道泥石流启动的力学本质为沟道中的堆积体所受剪切力克服摩擦力的过程。本文通过堆积体临界状态的受力平衡分析确立颗粒堆积体的临界平衡条件,发现堆积结构稳定性受坡度、颗粒内摩擦系数、堆积高度和水流速度共同作用。在此基础上通过文献中的系列实验数据分析拟合得到了一个泥石流启动的临界流速条件,相关结果和文献中的两个系列实验结果吻合良好,有望为沟道泥石流启动提供一定的支撑。     

基于极限平衡条件的沟道泥石流启动条件研究

沟道泥石流启动的力学本质为沟道中的堆积体所受剪切力克服摩擦力的过程。本文通过堆积体临界状态的受力平衡分析确立颗粒堆积体的临界平衡条件,发现堆积结构稳定性受坡度、颗粒内摩擦系数、堆积高度和水流速度共同作用。在此基础上通过文献中的系列实验数据分析拟合得到了一个泥石流启动的临界流速条件,相关结果和文献中的两个系列实验结果吻合良好,有望为沟道泥石流启动提供一定的支撑。     

Numerical simulations of granular free-surface flows using smoothed particle hydrodynamics

This paper presents a two-dimensional SPH model designed to simulate free-surface flows of dense granular materials. Smoothed particle hydrodynamics (SPH) is a mesh-free numerical method based on a Lagrangian discretization of the continuum mass and momentum conservation equations. The rheology of dense granular materials is modelled using a new local constitutive law recently proposed by Jop et al. (Nature, 2006). Of the viscoplastic class, this law is characterized by an apparent viscosity depending both on the local strain rate and local pressure. Validation test cases of the model in steady and unsteady configurations are presented. For steady cases (vertical chute flow and uniform free-surface layers on inclines), excellent agreement with analytical predictions is obtained. In the unsteady case, the simulations satisfactorily capture the dynamics of gravity-driven surges observed in experiments, including behaviours that are very specific to granular materials. Among the various parameters involved in the computations, the influence of SPH particle configuration within the flow and of the threshold viscosity used in the regularization of the constitutive yield criterion are particularly discussed.     

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.     

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

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

Transient finite element analysis of generalized Newtonian coextrusion flows in complex geometries

A two-dimensional transient finite element model capable of simulating problems related to two-layer polymer flows has been developed. This technique represents an effective tool which can be used to study the possibility of the onset of interfacial instability in coextrusion flows, considering melt rheology as well as the fluid–geometry interaction. A code has been developed to solve the transient problem of the flow of bi-component systems of Newtonian and generalized Newtonian fluids through parallel plates and complex geometries, such as: 2:1 abrupt expansion, 2:1 (30°) expansion, 4:1 abrupt contraction and 4:1 tapered (30°) contraction. Solutions are compared with experimental data from the literature and results provided by linear stability analysis (LSA) for the case of parallel plate flows. Numerical results are in agreement with LSA results for the parallel plate geometry cases studied. The expansion geometries tend to stabilize flows in the parallel plate section downstream of the expansion. Contractions may give rise to break-up of the interface depending on the flow conditions. © 1998 John Wiley & Sons, Ltd.     

基于应力波动的修正非局部流变模型

基于Pouliquen 提出的非局部流变模型,考虑颗粒流中某个位置重新排列引起的自激发过程,详细分析颗粒介质中应力波动幅值的概率密度分布形式以及剪切速率与体积分数的耦合作用,提出一种修正的非局部流变模型. 采用此修正非局部流变模型对斜面剪切颗粒流的流动特性进行了预测,颗粒流动的临界厚度、平均流动速度及剪切速率廓线的预测结果与实验结果定量吻合. 此修正模型的提出为复杂的密集颗粒流的描述和表征提供了一种新的研究思路.      

Normal stress effects in the gravity driven flow of granular materials

In this paper, we study the fully developed gravity-driven flow of granular materials between two inclined plates. We assume that the granular materials can be represented by a modified form of the second grade fluid where the viscosity depends on the shear rate and volume fraction and the normal stress coefficients depend on the volume fraction. We also propose a new isotropic (spherical) part of the stress tensor which can be related to the compactness of the (rigid) particles. This new term ensures that the rigid solid particles cannot be compacted beyond a point, namely when the volume fraction has reached the critical/maximum packing value. The numerical results indicate that the newly proposed stress tensor has obvious and physically meaningful effects on both the velocity and the volume fraction fields.     

Micromechanical definition of an entropy for quasi-static deformation of granular materials

A micromechanical theory is formulated for quasi-static deformation of granular materials, which is based on information theory. A reasoning is presented that leads to the definition of an information entropy that is appropriate for quasi-static deformation of granular materials. This definition is based on the hypothesis that relative displacements at contacts with similar orientations are independent realisations of a random variable. This hypothesis is made plausible based on the results of Discrete Element simulations. The developed theory is then used to predict the elastic behaviour of granular materials in terms of micromechanical quantities. The case considered is that of two-dimensional assemblies consisting of non-rotating particles with an elastic contact constitutive relation. Applications of this case are the initial elastic (small-strain) deformation of granular materials. Theoretical results for the elastic moduli, relative displacements, energy distribution and probability density functions are compared with results obtained from the Discrete Element simulations for isotropic assemblies with various average numbers of contacts per particle and various ratios of tangential to normal contact stiffness. This comparison shows that the developed information theory is valid for loose systems, while a theory based on the uniform-strain assumption is appropriate for dense systems.     

What do dry granular flows tell us about dense non-Brownian suspension rheology?

Recent years have seen significant progress in our understanding of the rheology of dry granular materials. A scale invariance of equations of motion in the rigid grain limit has helped identifying dimensionless quantities which govern flow. After a review of recent results on dry granular materials, we show how this same dimensional analysis carries over to the case of dense, non-Brownian suspensions. Our review is based on compiled data from various sets of numerical simulations, using both molecular and contact dynamics. It covers the breakdown of kinetic theory, which arises when contact forces dominate collisional forces in inertial flows, and the approach to dense flows up to the quasi-static limit. We show that simple invariance arguments permit to clarify the conditions of occurrence of viscous and inertial scaling in suspensions, and show in particular that both may occur in dense flows even in the presence of significant contact forces, up to the jamming limit. Some implications of the properties of steady uniform quasistatic granular flow under constant pressure for very dense suspensions near the maximum concentration are also discussed.     

DISCRETE AND CONTINUUM MODELLING OF GRANULAR FLOW

This paper analyses three popular methods simulating granular flow at different time and length scales： discrete element method （DEM）, averaging method and viscous, elastic-plastic continuum model. The theoretical models of these methods and their applications to hopper flows are discussed. It is shown that DEM is an effective method to study the fundamentals of granular flow at a particle or microscopic scale. By use of the continuum approach, granular flow can also be described at a continuum or macroscopic scale. Macroscopic quantities such as velocity and stress can be obtained by use of such computational method as FEM. However, this approach depends on the constitutive relationship of materials and ignores the effect of microscopic structure of granular flow. The combined approach of DEM and averaging method can overcome this problem. The approach takes into account the discrete nature of granular materials and does not require any global assumption and thus allows a better understanding of the fundamental mechanisms of granular flow. However, it is difficult to adapt this approach to process modelling because of the limited number of particles which can be handled with the present computational capacity, and the difficulty in handling non-spherical particles. Further work is needed to develoo an aoorooriate aooroach to overcome these problems.     

Transient and steady state behaviors of rapid granular shear flows

Rapid granular shear flow is a classical example in granular materials which exhibits fluid-like behavior solely or fluid-like and solid-like behaviors simultaneously. We have performed experiments on annular granular shear flows using monodisperse steel spheres with 2-mm and 3-mm diameters. We discuss some transient (static failure of material due to shear) and steady-state (completely sheared and partially sheared flows) features of rapid granular flows. Our results map out the boundary between overall compaction and dilation that characterizes multiphase behavior of the system. Dimensionless stresses rise as more material is loaded into the system. On the other hand, increasing the compressive force or shear rate is equal to easing the process of shearing by decreasing the ratio of shear to normal stress.