圆柱管道内运动的小球受到的黏滞阻力的数值分析

为了更好地运用落球法测量研究流体的黏滞系数,研究小球在黏性流体中下落的受力情况,本文对小球在充满黏性流体的圆柱管道的下落过程进行分析。利用COMSOL4.4仿真模拟,建立了合理的仿真模型,并分析了小球受到的黏滞阻力与小球的大小、下落位置的关系。结果表明:选择速度项二阶近似、压强项一阶近似的离散化方法,可以得到和理论值非常相符的仿真结果;当下落过程中小球球心始终在圆柱轴线上时,小球受到的黏滞阻力相对于Stokes力的修正系数,是小球半径与圆柱管道半径的比例函数,本文得到了更大范围的符合理论解的修正系数;当下落过程中小球的球心偏离圆柱轴线时,对于同样大小的小球,黏滞阻力、压强力、黏性力均随着球心到轴线的距离先减小后增大,且具有不同的极小值点。     

考虑等效曲率的超二次曲面单元非线性接触模型

基于连续函数包络的超二次曲面单元可有效地描述自然界和工业生产中的非球体颗粒形态, 并通过非线性迭代方法精确计算单元间的接触力. 对于具有复杂几何形态的超二次曲面单元, 线性接触模型不能准确地计算不同接触模式下的作用力. 考虑超二次曲面单元相互作用时不同颗粒形状及表面曲率的影响, 本文发展了相应的非线性黏弹性接触模型. 该模型将不同接触模式下的法向刚度和黏滞力统一表述为单元间局部接触点处等效曲率半径的函数; 切向接触作用则借鉴基于Mohr-Coulomb摩擦定律的球体单元非线性接触模型的计算方法. 为检验超二次曲面单元接触模型的可靠性, 对球形颗粒间的法向碰撞、椭球体颗粒间的斜冲击过程、圆柱体的静态堆积和椭球体的动态卸料过程进行离散元模拟, 并与有限元数值结果及试验结果进行对比验证. 计算表明, 考虑接触点处等效曲率半径的超二次曲面非线性接触模型可准确地计算单元间的接触碰撞作用, 并合理地反映非球形颗粒体系的运动规律. 在此基础上进一步分析了不同长宽比和表面尖锐度对卸料过程中颗粒流动特性的影响, 为非球形颗粒材料的流动特性分析提供了一种有效的离散元方法.      

考虑等效曲率的超二次曲面单元非线性接触模型

基于连续函数包络的超二次曲面单元可有效地描述自然界和工业生产中的非球体颗粒形态,并通过非线性迭代方法精确计算单元间的接触力.对于具有复杂几何形态的超二次曲面单元,线性接触模型不能准确地计算不同接触模式下的作用力.考虑超二次曲面单元相互作用时不同颗粒形状及表面曲率的影响,本文发展了相应的非线性黏弹性接触模型.该模型将不同接触模式下的法向刚度和黏滞力统一表述为单元间局部接触点处等效曲率半径的函数;切向接触作用则借鉴基于Mohr-Coulomb摩擦定律的球体单元非线性接触模型的计算方法.为检验超二次曲面单元接触模型的可靠性,对球形颗粒间的法向碰撞、椭球体颗粒间的斜冲击过程、圆柱体的静态堆积和椭球体的动态卸料过程进行离散元模拟,并与有限元数值结果及试验结果进行对比验证.计算表明,考虑接触点处等效曲率半径的超二次曲面非线性接触模型可准确地计算单元间的接触碰撞作用,并合理地反映非球形颗粒体系的运动规律.在此基础上进一步分析了不同长宽比和表面尖锐度对卸料过程中颗粒流动特性的影响,为非球形颗粒材料的流动特性分析提供了一种有效的离散元方法.     

基于介观力学信息的颗粒材料损伤——愈合与塑性宏观表征

本文在二阶计算均匀化框架下提出了颗粒材料损伤--愈合与塑性的多尺度表征方法. 颗粒材料结构在宏观尺度模型化为梯度Cosserat连续体,在其有限元网格的每个积分点处定义具有离散颗粒介观结构的表征元. 建立了表征元离散颗粒系统的非线性增量本构关系. 表征元周边介质作用于表征元边界颗粒的增量力与增量力偶矩以表征元边界颗粒的增量线位移与增量转动角位移、当前变形状态下表征元离散介观结构弹性刚度、以及凝聚到表征元边界颗粒的增量耗散摩擦力表示. 基于平均场理论与Hill定理,导出了基于介观力学信息的梯度Cosserat连续体增量非线性本构关系. 在等温热动力学框架下定义了表征颗粒材料各向异性损伤--愈合和塑性的损伤、愈合张量因子与综合损伤、愈合效应的净损伤张量因子和塑性应变. 此外,定义了损伤和塑性耗散能密度与愈合能密度,以定量比较材料损伤、愈合、塑性对材料失效的效应. 应变局部化数值例题结果显示了所建议的颗粒材料损伤--愈合--塑性表征方法的有效性.      

基于介观力学信息的颗粒材料损伤-愈合与塑性宏观表征

本文在二阶计算均匀化框架下提出了颗粒材料损伤-愈合与塑性的多尺度表征方法.颗粒材料结构在宏观尺度模型化为梯度Cosserat连续体,在其有限元网格的每个积分点处定义具有离散颗粒介观结构的表征元.建立了表征元离散颗粒系统的非线性增量本构关系.表征元周边介质作用于表征元边界颗粒的增量力与增量力偶矩以表征元边界颗粒的增量线位移与增量转动角位移、当前变形状态下表征元离散介观结构弹性刚度、以及凝聚到表征元边界颗粒的增量耗散摩擦力表示.基于平均场理论与Hill定理,导出了基于介观力学信息的梯度Cosserat连续体增量非线性本构关系.在等温热动力学框架下定义了表征颗粒材料各向异性损伤-愈合和塑性的损伤、愈合张量因子与综合损伤、愈合效应的净损伤张量因子和塑性应变.此外,定义了损伤和塑性耗散能密度与愈合能密度,以定量比较材料损伤、愈合、塑性对材料失效的效应.应变局部化数值例题结果显示了所建议的颗粒材料损伤-愈合-塑性表征方法的有效性.     

弹塑性自黏结颗粒聚合体碰撞破损的离散元法模拟研究

利用颗粒离散元方法,并结合接触力学理论,对自黏连、弹塑性颗粒聚合体碰撞破损的细观力学机理进行了模拟研究。在颗粒间塑性变形存在的条件下,研究了颗粒物碰撞损伤。将弹性颗粒体碰撞与弹塑性颗粒体碰撞结果作了比较。结果表明:其他同等条件下,与弹性颗粒聚合体相比,弹塑性颗粒聚合体的损伤模式多为衰变或解体,而非断裂,其原由可以归结为塑性变形引起的额外能量损耗;其他同等条件下,与弹性颗粒聚合体比较而言,弹塑性颗粒聚合体碰撞加载阶段平板撞击力波动振幅较小,加载时间较长,撞击力峰值较大,且峰值出现的时间较晚,而损伤率却较大,其原因可以归结为颗粒间塑性变形对动能损耗、局部结构变化、内部力传播的影响。     

颗粒层横向推移过程离散元仿真与参数分析

针对颗粒层横向推移过程,展开了物理试验及离散元仿真分析。采用物理试验方法对颗粒的物性参数进行了测定,并实测了颗粒层推移过程的推移力大小,进行了颗粒层横向推移过程离散元仿真模型的调试,通过仿真-试验结果的对比分析验证了仿真模型的合理性。利用所建立的仿真模型分析出:推移力斜率与静摩擦系数呈近似正比关系,颗粒堆积坡角度也随静摩擦系数的增大而增大;发现了摩擦系数增大造成推移力增大的力链机理,推板倾角增大会降低推移力增大速率,使推移力转变点位置靠近终点墙及增大竖直方向附加力。     

材料黏滞系数与损耗因子的频率相关性研究

针对黏弹性材料KV阻尼模型的黏滞系数与复阻尼模型的损耗因子间的关系,由单自由度体系的结构动力学分析,并基于结构振动响应的一致性,推导建立了黏滞系数与损耗因子在结构线性稳态简谐振动和自由振动时的一般关系式;并利用该关系式,试验研究了纤维混凝土材料黏滞系数和损耗因子的频率相关性.结果表明,黏滞系数与损耗因子间的关系在稳态简谐振动和自由振动时的表达形式相同,只是频率取值不同;纤维混凝土的损耗因子和黏滞系数都随频率增加而降低,且在0.5～1.0Hz频段降幅显著,而后渐趋平缓;相比于素混凝土,纤维混凝土的黏滞系数和损耗因子与激振频率的相关性更强.试验所得纤维混凝土频率相关的黏滞系数、损耗因子及推导所建立的两参数关系式为构建物理意义明确且又便于结构振动反应分析的阻尼系数或阻尼矩阵奠定了基础.      

无胶结材料的离散元微观参数交互作用研究

为研究无胶结材料离散元参数之间的交互影响,采用抗转动模型进行三轴排水试验模拟。首先,采用完全试验设计方法和析因试验设计方法,分别对宏微观参数相依规律以及参数之间的交互作用进行定性分析;然后,采用多元非线性回归分析方法对其进行定量分析,并提出了考虑交互作用的宏微观参数换算理论公式。着重探讨了无胶结材料微观参数中的刚度比α=k_s/k_n(0.3~0.9)、滚动刚度系数β(0.1~3)、颗粒内摩擦角θ(10~40)、法向接触刚度k_n(1×10~7~7×10~7)以及塑性弯矩系数η(0.1~3)与材料宏观力学性质峰值应力σ_p、残余强度σ_r、弹性模量E、泊松比ν以及内摩擦角φ之间的相关关系。研究结果表明,(1)颗粒材料的宏观强度参数主要受颗粒内摩擦角以及抗转动参数影响,而其变形参数主要由颗粒接触刚度及接触刚度比控制;(2)各个微观参数之间就存在不同程度的交互作用,且交互作用对试样宏观性质将产生一定影响;(3)对峰值强度、残余强度和宏观内摩擦角而言,两抗转动参数之间以及抗转动参数与颗粒内摩擦角之间存在较强交互作用。对于宏观弹性模量和泊松比而言,颗粒接触刚度和颗粒接触刚度比与其他微观参数产生较强的交互作用。研究成果对合理科学地确定离散元模型的微观参数具有一定的指导意义。     

附设黏滞阻尼器的装配式悬臂挡土墙动力响应研究

为提高装配式悬臂挡土墙抗震性能,将黏滞阻尼器附设在装配式悬臂挡土墙上,形成一种新型装配式悬臂阻尼挡土墙,通过黏滞阻尼器耗散地震能量达到减震目的。为研究该新型挡土墙的动力响应,采用ABAQUS软件,建立了两组三维有限元对比模型,一组为无滑体的阻尼悬臂挡土墙及对比模型,另一组为包含滑坡体的新型阻尼挡土墙及对比模型。对模型施以Vrancea地震波,研究了附设黏滞阻尼器的装配悬臂挡土墙抗震性能,并进行参数分析,讨论了混凝土强度、装配尺寸、阻尼器关键参数对其力学性能的影响。结果表明:附设黏滞阻尼器的装配悬臂挡土墙相较于对比模型,立板最大荷载减小约50%,最大应力减小约60%,立板最大水平位移由6.69mm减小至1.24mm,抗震性能显著提高;随着立板厚度、阻尼系数的增大,附设黏滞阻尼器的装配悬臂式挡土墙抗震性能有所提高,而混凝土强度、阻尼指数需要合理选择。     

Revisiting rolling and sliding in two-dimensional discrete element models

It has long been recognized that the rotation of single particles plays a very important role in simulations of granular flow using the discrete element method (DEM). Many researchers have also pointed out that the effect of rolling resistance at the contact points should be taken into account in DEM simulations. However, even for the simplest case involving two-dimensional circular particles, there is no agreement on the best way to define rolling and sliding, and different definitions and calculations of rolling and sliding have been proposed. It has even been suggested that a unique rolling and sliding definition is not possible. In this paper we assess results from previous studies on rolling and sliding in discrete element models and find that some researchers have overlooked the effect of particles of different sizes. After considering the particle radius in the derivation of rolling velocity, all results reach the same outcome: a unique solution. We also present a clear and simple derivation and validate our result using cases of rolling. Such a decomposition of relative motion is objective, or independent of the reference frame in which the relative motion is measured.     

Effect of solid mass flux on anisotropic gas–solid flow in risers determined with an LES-SOM model

Within the framework of the two-fluid approach, gas was treated with a large-eddy simulation and a sub-grid-scale (SGS) turbulent kinetic energy model while particles were treated with a second-order-moment method to describe the anisotropy of the fluctuating velocity. A modified Simonin model was derived for the gas–solid interphase fluctuating energy transfer. The anisotropic gas–solid flow in a circulating fluidized bed was investigated. Predictions were in good agreement with experimental data. The distributions of the second-order moment of particles and SGS-turbulent kinetic energy of gas were simulated at different solid mass fluxes. The effects of the solid mass flux on the particle second-order moment, particle anisotropic behavior, gas SGS-turbulent kinetic energy and gas SGS energy dissipation were analyzed for the circulating fluidized bed.     

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

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

基于碰撞模型的斜坡滚石颗粒速度预测

滑坡滚石灾害是西部山区常见的地质灾害类型,具有突发性和随机性强的特点,是山区地质灾害预测和防治工作的重点和难点. 本文基于颗粒接触理论,考虑影响斜坡滚石碰撞过程中的随机因素,建立了用于预测斜坡滚石颗粒碰撞后速度的理论模型. 根据冲量及冲量矩定理建立滚石颗粒碰撞基本方程,得到斜坡滚石颗粒碰撞后反弹速度的解析解. 结果表明:斜坡滚石碰撞后反弹速度的解析解包含了坡角、坡体上被碰颗粒速度以及角度、入射速度和角度以及撞击角度等随机因素;当考虑入射滚石颗粒与坡体上被碰颗粒的撞击角度变化时,模型预测结果与试验结果吻合较好;本文进一步预测了滚石颗粒碰撞后颗粒反弹线速度和角度以及反弹旋转角速度的概率分布情况. 结果显示,反弹颗粒速度和角度以及反弹旋转角速度的概率分布均服从高斯分布;当坡体上被碰颗粒速度和坡角发生变化时,其对反弹颗粒速度和角度以及反弹旋转角速度概率分布定性上没有影响,但是对概率分布的中心参数有显著影响.      

Direct numerical simulation of a particle-laden flow in a flat plate boundary layer

In this paper, a direct numerical simulation of particle-laden flow in a flat plate boundary layer is performed, using the Eulerian–Lagrangian point-particle approach. This is, as far as we know, the first simulation of a particle-laden spatially-developing turbulent boundary layer with two-way coupling. A local minimum of the particle number density is observed in the close vicinity of the wall. The present simulation results indicate that the inertial particles displace the quasi-streamwise vortices towards the wall, which, in turn, enhance the mean streamwise fluid velocity. As a result, the skin-friction coefficient is increased whereas the boundary layer integral thicknesses are reduced. The presence of particles augments the streamwise fluctuating velocity in the near-wall region but attenuates it in the outer layer. Nevertheless, the wall-normal and spanwise velocity fluctuations are significantly damped, and so is the Reynolds stress. In addition, the combined effect of a reduced energy production and an increased viscous dissipation leads to the attenuation of the turbulent kinetic energy.     

稠密气固两相流各向异性颗粒相矩方法

基于气体分子动力学和颗粒动理学方法,考虑颗粒速度脉动各向异性,建立颗粒相二阶矩模型.应用初等输运理论,对三阶关联项进行模化和封闭.考虑颗粒与壁面之间的能量传递和交换,建立颗粒相边界条件模型.采用Koch等计算方法模拟气固脉动速度关联矩.考虑气体-颗粒间相互作用,建立稠密气体-颗粒流动模型.数值模拟提升管内气固两相流动特性,模拟结果表明提升管内颗粒相湍流脉动具有明显的各向异性.预测颗粒速度、浓度和颗粒脉动速度二阶矩与Tartan等实测结果相吻合.模拟结果表明轴向颗粒速度脉动强度约为平均颗粒相脉动强度的1.5倍,轴向颗粒脉动能大约是径向颗粒脉动能3.0倍.     

基于离散单元法的颗粒阻尼耗能减振特性研究

采用离散单元法并从能量耗散的角度研究颗粒阻尼对系统减振特性的影响。建立了颗粒介质细观下的法向、切向和滚动方向的粘弹性接触模型和能量耗散模型,通过冲击激励和简谐激励下系统振动响应的多参数能量耗散分析来研究颗粒阻尼的耗能机理和减振特性。数值试验表明,颗粒介质可以在一个较宽的振动幅值范围内有效的发挥其阻尼效应,其耗能具有阶梯状周期性的特点。填充率是影响颗粒阻尼耗能减振效果的主要工程可控参数并对系统共振频率产生重大影响,当填充率接近极值时,系统出现无阻尼共振及共振频率超出无颗粒系统固有频率的现象。系统在最优填充率下共振时,颗粒与箱体保持恒定相位差的超振幅稳态运动。较小粒径的颗粒可以提高能量耗散率并使振动系统更快趋向静平衡状态,而恢复系数和摩擦系数则对法向和切向耗能的比值有较大影响。     

Numerical simulation of polygonal particles moving in incompressible viscous fluids

A two-dimensional coupled lattice Boltzmann immersed boundary discrete element method is introduced for the simulation of polygonal particles moving in incompressible viscous fluids. A collision model of polygonal particles is used in the discrete element method. Instead of a collision model of circular particles, the collision model used in our method can deal with particles of more complex shape and efficiently simulate the effects of shape on particle–particle and particle–wall interactions. For two particles falling under gravity, because of the edges and corners, different collision patterns for circular and polygonal particles are found in our simulations. The complex vortexes generated near the corners of polygonal particles affect the flow field and lead to a difference in particle motions between circular and polygonal particles. For multiple particles falling under gravity, the polygonal particles easily become stuck owing to their corners and edges, while circular particles slip along contact areas. The present method provides an efficient approach for understanding the effects of particle shape on the dynamics of non-circular particles in fluids.     

Transient rolling friction model for discrete element simulations of sphere assemblies

The rolling resistance between a pair of contacting particles can be modeled with two mechanisms. The first mechanism, already widely addressed in the DEM literature, involves a contact moment between the particles. The second mechanism involves a reduction of the tangential contact force, but without a contact moment. This type of rotational resistance, termed creep-friction, is the subject of the paper. Within the creep-friction literature, the term “creep” does not mean a viscous mechanism, but rather connotes a slight slip that accompanies rolling. Two extremes of particle motions bound the range of creep-friction behaviors: a pure tangential translation is modeled as a Cattaneo–Mindlin interaction, whereas prolonged steady-state rolling corresponds to the traditional wheel–rail problem described by Carter, Poritsky, and others. DEM simulations, however, are dominated by the transient creep-friction rolling conditions that lie between these two extremes. A simplified model is proposed for the three-dimensional transient creep-friction rolling of two spheres. The model is an extension of the work of Dahlberg and Alfredsson, who studied the two-dimensional interactions of disks. The proposed model is applied to two different systems: a pair of spheres and a large dense assembly of spheres. Although creep-friction can reduce the tangential contact force that would otherwise be predicted with Cattaneo–Mindlin theory, a significant force reduction occurs only when the rate of rolling is much greater than the rate of translational sliding and only after a sustained period of rolling. When applied to the deviatoric loading of an assembly of spheres, the proposed creep-friction model has minimal effect on macroscopic strength or stiffness. At the micro-scale of individual contacts, creep-friction does have a modest influence on the incremental contact behavior, although the aggregate effect on the assembly's behavior is minimal.     

Impact comminution of solids due to local kinetic energy of high shear strain rate: I. Continuum theory and turbulence analogy

The modeling of high velocity impact into brittle or quasibrittle solids is hampered by the unavailability of a constitutive model capturing the effects of material comminution into very fine particles. The present objective is to develop such a model, usable in finite element programs. The comminution at very high strain rates can dissipate a large portion of the kinetic energy of an impacting missile. The spatial derivative of the energy dissipated by comminution gives a force resisting the penetration, which is superposed on the nodal forces obtained from the static constitutive model in a finite element program. The present theory is inspired partly by Grady's model for expansive comminution due to explosion inside a hollow sphere, and partly by analogy with turbulence. In high velocity turbulent flow, the energy dissipation rate gets enhanced by the formation of micro-vortices (eddies) which dissipate energy by viscous shear stress. Similarly, here it is assumed that the energy dissipation at fast deformation of a confined solid gets enhanced by the release of kinetic energy of the motion associated with a high-rate shear strain of forming particles. For simplicity, the shape of these particles in the plane of maximum shear rate is considered to be regular hexagons. The particle sizes are assumed to be distributed according to the Schuhmann power law. The condition that the rate of release of the local kinetic energy must be equal to the interface fracture energy yields a relation between the particle size, the shear strain rate, the fracture energy and the mass density. As one experimental justification, the present theory agrees with Grady's empirical observation that, in impact events, the average particle size is proportional to the (−2/3) power of the shear strain rate. The main characteristic of the comminution process is a dimensionless number B_a (Eq. (37)) representing the ratio of the local kinetic energy of shear strain rate to the maximum possible strain energy that can be stored in the same volume of material. It is shown that the kinetic energy release is proportional to the (2/3)-power of the shear strain rate, and that the dynamic comminution creates an apparent material viscosity inversely proportional to the (1/3)-power of that rate. After comminution, the interface fracture energy takes the role of interface friction, and it is pointed out that if the friction depends on the slip rate the aforementioned exponents would change. The effect of dynamic comminution can simply be taken into account by introducing the apparent viscosity into the material constitutive model, which is what is implemented in the paper that follows.