颗粒阻塞有限元分析：Beverloo 卸料率与料斗稳定颗粒拱 DEM 快照

项目过程

本项目为某粉体处理设备公司的料斗卸料流动性评估项目，需求是评估玻璃珠颗粒在不同出料口尺寸下的卸料率 Q、堵塞概率，以及稳定颗粒拱的 DEM 模拟快照，给出料斗出料口设计建议。客户给定：玻璃珠 d=4.4 mm、密度 ρp=940 kg/m³（接近聚苯乙烯）、料斗倾角 30°、出料口 Do=20 mm、可选 Do=30~80 mm、流动状态为重力驱动重力流。

本项目采用 LIGGGHTS + EDEM 离散元方法（DEM）做颗粒流仿真 + Python 自研脚本做 Beverloo 公式拟合。DEM 部分：Hertz-Mindlin 无黏附接触模型，恢复系数 e=0.8、摩擦系数 μ=0.3、滚动摩擦 μr=0.05；料斗 2D 模型 0.3 m 高度 + 颗粒初始填充 5000 颗，仿真时长 5 s；后处理：每 0.1 s 记录出料颗粒数与堵塞状态。Beverloo 公式 Q=Cpb√(g(Do−kd)^(5/2))，C≈0.58、k≈1.4 用于玻璃珠。

案例分析

Beverloo 公式（蓝色实线）与 EDEM 仿真散点高度吻合：当出料口宽度 Do/d<5 时（阴影”无流动”区），质量流率 Q≈0 g/s，颗粒完全堵塞（概率 95%）；Do/d=58 时，Q 从 1 g/s 快速上升到 7 g/s；Do/d=812 时，Q 加速上升（1321 g/s）；Do/d=1215 时，Q 增速放缓（21~35 g/s）；Do/d=15 时，Q≈35 g/s 接近极限。堵塞概率（绿色虚线）从 Do/d=5 时的 90% 单调下降到 Do/d=15 时的 <2%。

DEM 颗粒云图清晰显示料斗稳定颗粒拱形成过程：彩色条带代表颗粒深度（顶部黄色 → 底部红色），出料口上方形成稳定的”5~6 颗粒宽”颗粒拱，支撑上方颗粒不塌落；拱脚位于出料口两侧壁面，拱形呈典型”圆拱”。颗粒拱一旦形成，系统进入”堵塞稳态”——本快照对应 t=2 s（出料后 2 s），拱结构稳定 1.5 s 后因外力微扰坍塌，颗粒继续流动形成”间歇流”模式（jamming-rattling flow）。

总结教训

本项目最容易踩的坑是颗粒数量。5000 颗玻璃珠在 30 mm 出料口下不足以准确统计堵塞概率。本项目通过 30 次重复试验（每次更换初始随机种子）累加 15 万颗样本，堵塞概率误差 <5%。

另一坑是接触模型参数。客户最初给的摩擦系数 μ=0.5 偏高，导致仿真中颗粒过度自锁堵塞概率高估 20%。修正为 μ=0.3 后与 Beverloo 公式一致。

项目经验

经验一：客户最初设计出料口 Do=20 mm，对应 Do/d≈4.5（恰好在堵塞边界内）。本项目证明必须 Do≥28 mm（Do/d=6.4）才能保证 90% 不堵塞。建议客户将出料口改为 35 mm（Do/d=8），流动率 7 g/s 完全满足工况要求。

经验二：颗粒拱的形成与出料口宽度+料斗壁面摩擦直接相关。如果壁面摩擦系数从 μ=0.3 降到 μ=0.1（抛光壁面），临界 Do/d 从 5 降到 3.5，出料口可缩小到 15 mm。

经验三：Beverloo 公式与 DEM 误差 <5% 的情况仅适用于球形颗粒；对不规则颗粒（长条、扁平），C 系数从 0.58 降到 0.40，Q 下降 30%。

有限元仿真