颗粒毛细效应影响因素的离散元分析

颗粒毛细效应是指将一根细管插入填充有颗粒物质的容器中并对管施加竖直振动时颗粒在管内上升并最终达到一个稳定的高度的现象,该现象为颗粒物料的逆重力输运提供了一种潜在的技术途径.为探究颗粒毛细效应的影响因素,采用离散元方法,模拟再现了颗粒毛细效应过程,展示了不同管径下颗粒竖直方向速度演变特性,考察了不同容器宽度和振动条件下颗粒最终毛细上升高度随管径的演变规律.结果表明,在容器宽度与粒径比为40、管振幅与粒径比为14.33、管振动频率为12 Hz情况下,管径与粒径比D/d=3.33时,管内颗粒堵塞严重,使得颗粒上升缓慢,并造成颗粒柱中断; D/d=8.33时,起初毛细上升高度增加迅速,随后毛细上升高度的增大逐渐减缓,管内颗粒在管径方向几乎不存在速度梯度; D/d=15时,随着颗粒毛细上升高度的增大,管内颗粒柱分离为速度截然不同的两层,上层颗粒在管径方向几乎不存在速度梯度,而下层颗粒存在明显的速度梯度.研究还发现,在毛细效应能够发生的管径范围内,存在一个对应于颗粒最终毛细上升高度最大值的临界管径,当管径小于临界管径时,颗粒最终毛细上升高度随管径的增大而增大,当管径大于临界管径时,颗粒最终毛细上升...     

基于水平集-离散元方法的球谐函数颗粒材料缓冲性能分析

在自然环境与工业领域中,颗粒材料是一种常见的缓冲材料,其中大量形态各异的非球形颗粒表现出复杂的力学特性并应用于不同工程领域。本文采用球谐函数构造不同球面度和表面凹凸特性的非规则颗粒,通过水平集方法计算球谐函数颗粒间的接触点和碰撞力,并对冲击过程中球形和凹形颗粒的缓冲性能进行离散元分析。数值结果表明,颗粒床厚度、冲击速度和颗粒形状显著影响球谐函数颗粒材料的缓冲性能。颗粒床底部的冲击力峰值随着颗粒床厚度和表面凹凸性的增加而降低,同时冲击力峰值随着冲击速度和颗粒球面度的增加而增加。与球形颗粒相比,球谐函数颗粒具有凹凸表面和多接触点特性,这有利于冲击荷载向四周扩展并提高凹形颗粒的缓冲效果。     

颗粒毛细效应影响因素的离散元分析

颗粒毛细效应是指将一根细管插入填充有颗粒物质的容器中并对管施加竖直振动时颗粒在管内上升并最终达到一个稳定的高度的现象,该现象为颗粒物料的逆重力输运提供了一种潜在的技术途径.为探究颗粒毛细效应的影响因素,采用离散元方法,模拟再现了颗粒毛细效应过程,展示了不同管径下颗粒竖直方向速度演变特性,考察了不同容器宽度和振动条件下颗粒最终毛细上升高度随管径的演变规律.结果表明,在容器宽度与粒径比为40、管振幅与粒径比为14.33、管振动频率为12 Hz情况下,管径与粒径比D/d=3.33时,管内颗粒堵塞严重,使得颗粒上升缓慢,并造成颗粒柱中断; D/d=8.33时,起初毛细上升高度增加迅速,随后毛细上升高度的增大逐渐减缓,管内颗粒在管径方向几乎不存在速度梯度; D/d=15时,随着颗粒毛细上升高度的增大,管内颗粒柱分离为速度截然不同的两层,上层颗粒在管径方向几乎不存在速度梯度,而下层颗粒存在明显的速度梯度.研究还发现,在毛细效应能够发生的管径范围内,存在一个对应于颗粒最终毛细上升高度最大值的临界管径,当管径小于临界管径时,颗粒最终毛细上升高度随管径的增大而增大,当管径大于临界管径时,颗粒最终毛细上升高度随管径的增大而趋于减小;增大容器宽度,临界管径有所增大;增大振幅、适当提高频率能够有效促进临界管径的增大.     

基于超二次曲面的非球形离散单元模型研究

通过对颗粒体系接触过程的运动学和动力学分析,建立了一种基于超二次曲面的非球形离散单元模型,该模型避免了球形接触模型描述颗粒形状的局限性,使离散单元法更接近物理事实,并在此基础上提出了计算求解模型的数值方法,实现了对复杂形状的颗粒体系的模拟计算。将所建立的数值计算方法进行了编程实现,并对模型和算法进行了算例测试,证实了本文所建立的非球形离散单元模型的可行性和正确性。测试结果表明,本文的模型能够比较准确地模拟复杂颗粒体系的真实运动,可为复杂颗粒体系的模拟研究提供一种新的数值计算方法。     

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

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

基于颗粒离散元的沙纹演变大涡模拟研究

应用大涡模拟、"点球" 浸入边界法和基于"事件驱动" 模型和"弹簧-阻尼" 模型的颗粒离散元法, 数值模拟了明渠湍流中沙纹的演变. 通过对不同谢尔兹数下无沙纹床面的推移质输沙率进行计算, 并与经典输沙率公式进行对比, 验证了模型的精度和可靠性. 随后, 对明渠湍流中沙纹床面的演变过程进行了数值模拟, 计算了有沙纹床面的推移质输沙率、沙纹长度和高度、等效床面高度的最大值、最小值和平均值、沙纹形状阻力、体积流速随时间的变化曲线. 研究发现:初始平整的床面在较短的时间内(tU_b=h≈100) 发展出数条沙纹, 随后沙纹逐渐发展, 在tU_b=h为1 600~2 000 时, 沙纹发生合并. 在沙纹数量不变的条件下, 沙纹高度随时间近似呈线性增长, 而沙纹的长度的平均值却保持恒定. 随着沙纹高度的增大, 输沙率和体积流速逐渐降低, 沙纹形状阻力则逐渐增大;当沙纹发生合并时, 沙纹高度快速增加, 输沙率、体积流速和沙纹形状阻力也出现了大幅跳跃. 在同等的水流强度条件下, 有沙纹床面的输沙率小于平整床面的输沙率.      

竖直振动管中颗粒毛细效应的离散元模拟

将一根细管插入填充有颗粒的静止容器中并对管施加竖直振动,颗粒将在管内发生上升运动,并最终稳定在一定高度,这一现象与液体毛细效应类似,被称为颗粒毛细效应.为探究颗粒毛细效应过程中伴随的颗粒尺度动力学行为及机理,基于离散元方法建立颗粒运动模型,对颗粒毛细效应动力学过程和特性开展数值模拟研究.模拟再现了文献中实验得到的颗粒毛细效应全过程,给出了管内颗粒柱高度随时间的演变规律,结果表明,受到颗粒系统参数的影响,本模拟条件下颗粒毛细效应过程呈现单周期上升、倍周期上升和倍周期稳定三个阶段,在倍周期上升阶段颗粒柱上升速度逐渐减小,平缓过渡到稳定阶段.在此基础上,分析了管内颗粒速度场和填充率分布随时间的演变特性,揭示了颗粒毛细效应过程中由容器传输到管内的颗粒的占比分布.研究发现,管内不同高度位置颗粒的运动并不同步,随着管的振动,管内出现速度波,速度波的传播引起管内颗粒出现膨胀和压缩交替的情况,从而管内颗粒填充率随时间发生周期性波动;在上升阶段,越接近管壁由容器传输到管内的颗粒占比越大,在稳定阶段,管内上层颗粒的对流引起容器传输到管内的颗粒占比发生反转.      

不同重力场下颗粒冲击过程的离散元分析

开展不同重力场下颗粒材料冲击动力学研究有助于加深对颗粒运动机制的理解和深空探测任务的实施。本文采用离散元模拟对颗粒介质受球形冲击物的冲击过程进行了数值模拟,并与地球重力下冲击的试验结果进行对比验证。在此基础上,进一步研究了重力加速度对冲击物动力学的影响规律。计算结果表明,在所有重力加速度下,冲击物的穿透深度d与冲击速度v₀的关系可以用Poncelet模型表达;d与冲击物下落的总高度H表现为d～Hⁿ的幂律关系,当H<10 m时,d与H的幂率标度为0.322,而H>10 m时,d与H的幂率标度下降到0.211。此外,穿透深度小于冲击物半径时,重力加速度对冲击物减速过程无影响。在所有的重力加速度下,当冲击速度大于5 m/s时,冲击物的持续碰撞时间t_c是恒定的,且与重力的-1/2次方呈正比。     

大规模颗粒系统的精确缩尺和粗粒化离散元方法

计算效率是制约工程尺度大规模颗粒系统离散元计算发展的重要因素,现有的粗粒化处理方法局限于特定应用并且缺少一般的理论依据。本文采用量纲分析方法,描述了在精确缩尺系统中各物理量应当满足的缩放定律;通过在粗粒化系统和原始系统的代表性体积单元之间建立质量、动量和能量的近似守恒关系,采用多尺度的描述方法得到了粗粒化系统与原始系统之间宏观和细观两种不同尺度的缩放关系,即双尺度粗粒化模型;精确缩尺系统中得到的缩放定律及离散元接触模型处理方法,完全适用于粗粒化系统中细观颗粒层面相关物理量的缩放,通过筒仓侧壁压力和休止角两个算例对精确缩尺模型在粗粒化系统中的有效性进行了验证。     

颗粒材料内部滑动变形局部化的离散元模拟

颗粒材料在加载过程中表现出非常复杂的变形性质.在微观尺度上,单个颗粒的运动、单个孔隙胞元内的变形都是非均匀的,但也不是完全随机、没有规律可循的,而是呈现一定结构性的分布,如微带、剪切带等.本文用基于孔隙胞元的离散元方法对颗粒体进行双轴加载数值试验,模拟了以滑动变形表征的变形局部化现象.数值结果发现,对应加载过程中的不同...     

回转窑多粒径颗粒流动与偏析的DEM模拟与实验研究

为明晰回转窑内颗粒的运动行为及偏析机理,以绿豆、黄豆和黑豆为颗粒介质,依次对3种装填顺序下的颗粒流动过程进行离散元模拟与实验研究,以颗粒质量分数和平均粒度为判据,对颗粒偏析进行评价。结果表明,回转窑内颗粒流动区可分为自由滚落区、渗流呆滞区以及窑壁携带区,自由滚落区颗粒流速最大,而渗流呆滞区流速最小。窑内颗粒沿轴向输运过程发生径向偏析,形成夹层结构,小颗粒受渗流作用在渗流呆滞区中心形成内核,大粒径和中等粒径颗粒集中在自由滚落区和窑壁携带区。窑内颗粒力链分布不均匀,强力链分布于近窑壁区,弱力链分布于自由滚落区和渗流呆滞区,且渗流呆滞区力链细而密集。当窑头附近不同粒径颗粒存在轴向速度差时,颗粒在轴向发生掺混,并产生径向偏析。     

基于Voronoi切割算法的碎冰区构造及离散元分析

离散元方法广泛应用于海冰,特别是碎冰区的动力过程及其对海洋结构作用过程的数值模拟。为构造碎冰区中的冰块几何特性,基于二维Voronoi图方法对计算域进行随机切割以生成碎冰区中冰块的几何形态,并采用球体单元对每个碎冰块单元进行填充,从而确定碎冰区的初始分布场。在采用Voronoi图进行碎冰区构造时,可对冰块尺寸、几何形态和密集度等海冰参数进行设定。为确定冰块的不同几何规则度,综合采用排斥法和扰动法以定量地控制碎冰块几何形态从完全随机分布到规则分布的连续变换。为分析不同几何规则度下碎冰块的几何特性概率分布规律,对计算域内冰块的面积和边数等参数进行统计分析,从而可更合理地参数化控制初始冰场中碎冰块的几何特性。在此基础上,本文基于粘接-破碎的球体离散元方法对不同冰况下锥体结构的冰荷载进行了数值计算,讨论分析了碎冰区的海冰密集度、冰块面积和几何规则度对冰载荷的影响。     

Formulation of a semi-analytical approach based on Gurtin variational principle for dynamic response of general thin plates

4 semi-analytical approach for the dynamic response of general thin plates which employes finite element discretization in space domain and a series of representation in time domain is developed on the basis of Gurtin variational principles. The formulation of time series is also investigated so that the dynamic response of plates with arbitrary shape and boundary constraints can be achieved with adequate accuracy.Project supported by the National Natural Science Foundation of China