颗粒物质的分形特征与其物理力学性质的关系探讨

从颗粒流的一般本构方程出发,利用室内大型相对密度仪和直接剪切仪,针对粒径分形分布对孔隙率和内摩擦角的影响进行了实验研究。结果表明,在相同实体密度、弹性恢复系数和几何平均中值粒径条件下,随粒径分维数的增大,孔隙率增大,而内摩擦角、颗粒浓度、剪切应力系数与正应力系数均减小。进而分析了粒径分形特性对颗粒流动行为的影响,得出粒径分维数小的颗粒系统剪切应力大,系统的流动性能差,鲁棒性强;反之,分维数大的颗粒系统剪切应力小,系统的流动性能好,敏感性强。对系统组构的分形特征影响系统颗粒物质崩塌流动行为的机制进行了初步探索,为颗粒流系统的研究提供了新思路。     

颗粒振动分离的浮力计算折减修正法及数值模拟

颗粒材料在振动或流动情况下的分离行为如巴西果效应已经广为人知,为了解释这种分离行为已提出了渗透、对流、凝聚等一系列解释机制,但是尚未达成共识。本文研究颗粒上浮时所受的浮力,发现已有的大颗粒浮力计算方法无法解释大颗粒粒径改变时上浮行为不同这一典型现象。为此引入粒径比对已有的大颗粒浮力计算方法进行了修正,并证明了其合理性。此外,还研究了局部阻尼系数、振幅、颗粒密度对压强分布的影响。结果表明:所建议的修正办法能正确地计算不同粒径比下大颗粒受到的浮力;随着局部阻尼系数的增大,压强最大值增大,压强曲线下降趋势由直线变为曲线;随着颗粒密度的增大,压强最大值增大,但压强曲线下降趋势变化并不明显;随着加载振幅的增大,压强最大值与下降趋势变化均不明显。     

颗粒物对变压器油表面张力的影响研究

根据ISO4406标准,采用均匀设计法,配制了24组含Cu、Fe、Si O2的不同粒径、不同污染度的颗粒污染物油样并对其进行了表面张力测试,得到随着颗粒物含量的增加,油液表面张力值呈现增加的趋势。采用Takagi-Sugeno模糊原理对油液中不同粒径和含量的Cu、Fe、Si O2混合颗粒物与油样的表面张力性能的关系进行了建模。结果表明:根据所建立的不同粒径和含量的Cu、Fe、Si O2颗粒物与表面张力性能关系的T_S模糊模型,获取了验证集的决定系数和均方根误差分别是0.9124、0.016;油中Cu、Fe、Si O2混合颗粒物在粒径为5μm、15μm时对变压器油表面张力的影响较为显著;由于油液分子作用力与颗粒的吸附能出现饱和-动态平衡作用,随着颗粒数的增加油液表面张力的变化较大,油液表面张力的变化量分别为0.067MN/m、0.042MN/m。     

可侵蚀地表上方含尘大气运动的数值模型

在稀相气固两相流的双流体模型框架下，导出含尘大气运动的基本方程，其中的源项考虑相间双向耦合作用；通过对基本方程进行无量纲化，求得控制气相和弥散相动力学行为的相似准则，包括弗洛德数、罗斯贝数、颗粒雷诺数、颗粒惯性参数、颗粒剪切横移参数、颗粒旋转横移参数和沙尘质量载荷率等. 作为一个模型问题，研究可侵蚀地表上方充分发展湍流大气边界层流动对土壤颗粒的气动卷扬作用，其中计及气动阻力、萨弗曼升力和重力等因素的影响. 为了克服由于沙尘跃移运动轨迹交叉导致流动参数多值性所造成的困难，引进拉格朗日坐标下弥散相连续方程. 对两种不同风速和4种不同粒径条件下的沙尘运动进行数值模拟，给出相应情况中沙尘运动特性和密度分布剖面并讨论风速和粒径等参数的影响，还细致分析了含尘大气边界层中沙尘与气流之间的能量传递过程. 所得的方程、准则和方法可以为土壤风蚀和沙尘暴等自然灾害的预报提供理论基础.     

球形和方形湿颗粒团的碰撞模拟研究

本研究利用颗粒离散元方法,并结合接触力学理论,对微米级颗粒组成的湿颗粒聚合体碰撞破损的细观力学机理进行了模拟研究。模拟的颗粒聚合体分成球形和方形。将两种形状的颗粒聚合体碰撞结果做比较,颗粒之间的接触采用干、湿接触两种模型并存的方式来模拟湿颗粒团碰撞中的不可恢复变形过程。在碰撞速度保持不变的情况,考察了不同液体黏度对颗粒聚合体碰撞的影响;在液体黏度保持不变的情况,考察了不同碰撞速度对碰撞结果的影响。研究结果表明,湿颗粒团表现出了与干颗粒团完全不同的损伤模式,没有产生类似于弹性颗粒团的裂解损伤。湿颗粒团在碰撞中,它的内部液桥数量会在与平板碰撞接触的初期快速减少,然后慢慢累计,或会超过湿颗粒团在碰撞前的液桥数量。这一现象在此仿真研究中首次被发现。比如,液体黏度为50mPa·s的方形湿颗粒团以2.0m/s撞击平板时,其内部液桥数可由约9400降低至约8700后,逐渐增加至11300左右;由低黏度液体组成的湿颗粒团在碰撞后会发生反弹现象。而随着液体黏度的增加,恢复系数会趋近于0。     

两空间耦合下齿轮传动系统多稳态特性研究

通过将系统参数定义为参数变量, 构成参数空间,研究齿轮传动系统在参数空间和状态空间耦合下的非线性全局动力学特性,以及多参数、多初值和多稳态行为之间的关联特性.首先设计了一个两空间耦合下非线性系统多稳态行为的计算和辨识方法.其次,基于该方法并结合相图、Poincaré映射图、分岔图、最大Lyapunov指数、吸引域等,研究齿轮传动系统在不同参数平面上多稳态行为的存在区域和分布特性,以及多稳态行为在状态平面上的分布特性,揭示了参数平面和状态平面上系统可能隐藏的多稳态行为和分岔,并分析了多稳态行为的形成机理. 结果发现,两空间耦合下系统在参数平面上存在大量多稳态行为并呈"带状"分布, 状态平面上多稳态行为出现两种不同的侵蚀现象, 即内部侵蚀和边界侵蚀.分岔点或分岔曲线对初值的敏感性导致多稳态行为的出现.当齿侧间隙和误差波动在较小的范围内变化时,系统全局动力学特性受间隙和误差扰动的影响较小,受啮合频率的影响较大.两空间耦合下系统全局动力学特性变得丰富和复杂.      

尼龙粉末在SLS预热温度下的离散元模型参数确定及其流动特性分析

尼龙粉末是增材制造中常用的粉体材料,温度对其流动性有重要影响. 探索尼龙粉末增材制造预热温度下的流动性是研究选择性激光烧结(selective laser sintering, SLS)工艺中粉体铺展成形的基础. 选取SLS技术中的尼龙粉末为原材料,采用离散元数值方法,研究尼龙粉末的流动行为,是增材制造工艺数值模拟和铺粉工艺优化的研究热点. 以Hertz-Mindlin模型为基础,基于Hamaker理论模型和库伦定律,在尼龙粉末的接触动力学模型中引入范德华力和静电力,建立预热温度下尼龙粉末流动的离散元模型(discrete element method, DEM),通过对比相应实验结果,标定了该模型的参数. 对加热旋转圆筒中尼龙粉末流动过程进行了DEM数值模拟,校核了所建模型的正确性,并研究了粉体粒径分布对尼龙粉末流动特性的影响规律. 研究表明,尼龙粉末黏附力是静电力与范德华力的共同作用结果;随着粉体粒径的增大,尼龙粉末崩塌角增大,流动性增强;相对于高斯粒径分布,粒径均匀分布的尼龙粉末颗粒流动性更强. 研究结果可指导SLS中铺粉工艺的优化.      

两空间耦合下齿轮传动系统多稳态特性研究

通过将系统参数定义为参数变量,构成参数空间,研究齿轮传动系统在参数空间和状态空间耦合下的非线性全局动力学特性,以及多参数、多初值和多稳态行为之间的关联特性.首先设计了一个两空间耦合下非线性系统多稳态行为的计算和辨识方法.其次,基于该方法并结合相图、Poincaré映射图、分岔图、最大Lyapunov指数、吸引域等,研究齿轮传动系统在不同参数平面上多稳态行为的存在区域和分布特性,以及多稳态行为在状态平面上的分布特性,揭示了参数平面和状态平面上系统可能隐藏的多稳态行为和分岔,并分析了多稳态行为的形成机理.结果发现,两空间耦合下系统在参数平面上存在大量多稳态行为并呈"带状"分布,状态平面上多稳态行为出现两种不同的侵蚀现象,即内部侵蚀和边界侵蚀.分岔点或分岔曲线对初值的敏感性导致多稳态行为的出现.当齿侧间隙和误差波动在较小的范围内变化时,系统全局动力学特性受间隙和误差扰动的影响较小,受啮合频率的影响较大.两空间耦合下系统全局动力学特性变得丰富和复杂.     

尼龙粉末在SLS预热温度下的离散元模型参数确定及其流动特性分析

尼龙粉末是增材制造中常用的粉体材料,温度对其流动性有重要影响.探索尼龙粉末增材制造预热温度下的流动性是研究选择性激光烧结(selective laser sintering, SLS)工艺中粉体铺展成形的基础.选取SLS技术中的尼龙粉末为原材料,采用离散元数值方法,研究尼龙粉末的流动行为,是增材制造工艺数值模拟和铺粉工艺优化的研究热点.以Hertz-Mindlin模型为基础,基于Hamaker理论模型和库伦定律,在尼龙粉末的接触动力学模型中引入范德华力和静电力,建立预热温度下尼龙粉末流动的离散元模型(discrete element method,DEM),通过对比相应实验结果,标定了该模型的参数.对加热旋转圆筒中尼龙粉末流动过程进行了DEM数值模拟,校核了所建模型的正确性,并研究了粉体粒径分布对尼龙粉末流动特性的影响规律.研究表明,尼龙粉末黏附力是静电力与范德华力的共同作用结果;随着粉体粒径的增大,尼龙粉末崩塌角增大,流动性增强;相对于高斯粒径分布,粒径均匀分布的尼龙粉末颗粒流动性更强.研究结果可指导SLS中铺粉工艺的优化.     

时滞非线性悬架系统的混沌动力学特性研究

针对时滞减振控制的非线性悬架系统,建立其二自由度系统的动力学方程。首先,对动力系统进行了数值模拟,通过不同控制参数下系统的动力学行为的分岔图、相轨迹、庞加莱截面、功率谱图来研究时滞非线性悬架系统的混沌动力学行为。研究表明,基于系统参数和外在激励,选择适当的时滞控制参数,可避免系统在运行过程中出现混沌现象,改善系统的运行品质。然后,以主系统幅值均方根为目标函数,对系统进行优化得出减振效果最优时的时滞和反馈增益系数,并与无时滞时非线性悬架系统的主振幅响应进行比较。结果表明,时滞对非线性悬架系统减振和系统品质的改善是能够同时实现的。最后,研究了时滞控制参数变化对系统动力学行为的影响,研究发现,同一系统在不同时滞参数下其分岔形式以及通往混沌的形式具有着多样性,会出现倍周期分岔、Hopf分岔、阵发性分岔以及它们各自通往混沌的不同演化模式,这为实现悬架参数的优化控制提供了理论依据。     

三维圆柱型颗粒堆坍塌问题的全相态数值模拟

静态颗粒堆在重力作用下的坍塌问题, 是认识和理解许多人为过程和自然现象的基础. 采用传统方法进行模拟存在单颗粒追踪数量大、宏观模拟流变特性明显和相态演变复杂等计算难点. 本文从颗粒介质表现出不同相态的物理机理出发, 对全相态概念进行了定义并进行了区域划分. 根据颗粒介质的应力-应变关系及体积分数的不同, 通过确定不同相态之间的耦合关系和转化准则, 将描述各相态的现有理论有效结合起来, 建立了描述颗粒介质经历全部相态的耦合模型理论. 采用光滑离散颗粒流体动力学方法和离散单元法相耦合的策略, 对颗粒介质物理模型求解, 实现了对不同长径比下的三维圆柱型颗粒堆坍塌过程的数值模拟. 计算结果与实验结果吻合较好, 同时与离散单元法相比, 计算量得到了控制. 不仅捕捉到了不同参数影响下颗粒堆坍塌后沉积的不同现象, 同时获得了不同条件参数对颗粒堆坍塌后铺展特性的影响规律, 为揭示工业和自然界中广泛存在的颗粒介质复杂运动机理提供有效的支撑.      

基于介观结构的饱和与非饱和多孔介质有效应力

基于描述含液颗粒材料介观结构的Voronoi 胞元模型和离散颗粒集合体与多孔连续体间的介-宏观均匀化过程, 定义饱和与非饱和多孔介质有效应力. 导出了计及孔隙液压引起之颗粒体积变形的饱和多孔介质广义有效应力. 用以定义广义有效应力的Biot 系数不仅依赖于颗粒材料的多孔连续体固体骨架及单个固体颗粒的体积模量(材料参数),同时与固体骨架当前平均广义有效应力及单个固体颗粒的体积应变(状态量) 有关. 提出了描述非饱和多孔介质中非混和固体颗粒、孔隙液体和气体等三相相互作用的具介观结构的Voronoi 胞元模型.具体考虑在低饱和度下双联(binary bond) 模式的摆动(pendular) 液桥系统介观结构. 导出了基于介观水力-力学模型的非饱和多孔介质的各向异性有效应力张量与有效压力张量. 考虑非饱和多孔介质Voronoi 胞元模型介观结构的各向同性情况,得到了与非饱和多孔连续体理论中唯象地假定的标量有效压力相同的有效压力形式.但本文定义的与确定非饱和多孔介质有效应力和有效压力相关联的Bishop 参数由基于三相介观水力-力学模型, 作为饱和度、孔隙度和介观结构参数的函数导出,而非唯象假定.      

Separation efficiency of liquid–solid undergoing vibration based on breakage of liquid bridge

Vibrating separation is a significant method for liquid–solid separation. A typical example is the vibrating screen to dewater wet granular matter. The properties of granular matter and the vibrating parameters significantly affect the separation efficiency. This study investigates the effect of vibration parameters in separation based on the breakage of large-scale liquid bridge numerically by using a calibrated simulation model. Through analysing the simulation results, the liquid bridge shape and the volume between two sphere particles for various particle sizes and particle distances were studied in the static condition under the effect of gravity. The results show a general reducing trend of liquid bridge volume when the radius ratio of two particles increases, particularly when the ratio increases to 5. Additionally, a set of vibrating motion was applied to the liquid bridge in the simulation model. A group of experiments were also performed to validate the simulation model with vibration. Then, the effect of vibrating peak acceleration, distance between spheres and radius on the separation efficiency which was reflected by the residual water were investigated. It is found that separation efficiency increased obviously with the peak acceleration and the increase slowed down after the peak acceleration over 1 m/s².     

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.     

Two-fluid smoothed particle hydrodynamics simulation of submerged granular column collapse

In this paper, a two-fluid smoothed particle hydrodynamics (SPH) model, based on the mixture theory, is employed to investigate the complex interactions between the solid particles and the ambient water during the process of submerged granular column collapse. From the simulation, two regimes of the collapse, one being quick and the other being slow, are identified and the reasons of formation are analyzed. It is found that, a large internal friction angle of the granular phase, representing large drag force between solid particles, helps form the slow regime. Small hydraulic conductivity, representing large inter-phase drag force, also retards the collapse dramatically. Good agreements between our numerical results and other researchers’ numerical and experimental results are observed, which demonstrates the capability of the proposed two-fluid SPH approach in dealing with saturated water–soil mixture flows.     

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

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

A novel hybrid solid-like fluid-like (SLFL) method for the simulation of dry granular flows

This paper proposes a novel hybrid method to simulate the dry granular flow of materials over a wide range of inertial numbers that simultaneously covers the quasi-static and dense granular flow regimes. To overcome the lack of incremental objectivity whenever large deformations occur in solid-like regimes and to remove computational singularities in fluid-like regimes close to rest, the elastic–perfectly plastic theory based on the Drucker–Prager yield criterion is combined with the theory of dense granular flows. By implementing some new modifications at the boundaries and removing all ghost particles, smoothed particle hydrodynamics (SPH) is used as the framework for the method. A number of benchmark problems have been solved to show the capabilities of the new modified SPH method. Precise prediction of both location and pressure makes the modifications comparable with the previous works on SPH. Finally, the method is used to solve the classic 2D dry granular cliff collapse problem and to model dry granular material flow inside a rotary drum. The outcomes of the numerical simulation show good agreement with tabletop experiments and published results.     

Discrete element parameter calibration and the modelling of dragline bucket filling

The Discrete Element Method (DEM) is useful for modelling granular flow. The accuracy of DEM modelling is dependent upon the model parameter values used. Determining these values remains one of the main challenges. In this study a method for determining the parameters of cohesionless granular material is presented. The particle size and density were directly measured and modelled. The particle shapes were modelled using two to four spheres clumped together. The remaining unknown parameter values were determined using confined compression tests and angle of repose tests. This was done by conducting laboratory experiments followed by equivalent numerical experiments and iteratively changing the parameters until the laboratory results were replicated. The modelling results of the confined compression tests were mainly influenced by the particle stiffness. The modelling results of the angle of repose tests were dependent on both the particle stiffness and the particle friction coefficient. From these observations, the confined compression test could be used to determine the particle stiffness and with the stiffness known, the angle of repose test could be used to determine the particle friction coefficient. Usually DEM codes do not solve the equations of motion for so-called walls (non-granular structural elements). However, in this study a dynamic model of a dragline bucket is developed and implemented in a commercial DEM code which allows the dynamics of the walls to be modelled. The DEM modelling of large systems of particles is still a challenge and procedures to simplify and speed up the modelling of dragline bucket filling are presented. Using the calibrated parameters, numerical results of bucket filling are compared to experimental results. The model accurately predicted the orientation of the bucket. The model also accurately predicted the drag force over the first third of the drag, but predicted drag forces too high for the subsequent part of the drag.     

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

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

Investigation on effect of drag models on flow behavior of power-law fluid–solid two-phase flow in fluidized bed

In this study, a Eulerian-Eulerian two-fluid model combined with the kinetic theory of granular flow is adopted to simulate power-law fluid–solid two-phase flow in the fluidized bed. Two new power-law liquid–solid drag models are proposed based on the rheological equation of power-law fluid and pressure drop. One called model A is a modified drag model considering tortuosity of flow channel and ratio of the throat to pore, and the other called model B is a blending drag model combining drag coefficients of high and low particle concentrations. Predictions are compared with experimental data measured by Lali et al., where the computed porosities from model B are closer to the measured data than other models. Furthermore, the predicted pressure drop rises as liquid velocity increases, while it decreases with the increase of particle size. Simulation results indicate that the increases of consistency coefficient and flow behavior index lead to the decrease of drag coefficient, and particle concentration, granular temperature, granular pressure, and granular viscosity go down accordingly.