离散元方法揭示了扇形筒仓中谷物储存的流动机制及侧向压力特性

时间:2026年1月19日
来源:Powder Technology

编辑推荐:

分类储粮扇形粮仓充放粮流动模式与动态侧向压力分布研究。通过SolidWorks构建1000mm高250mm半径三区隔扇形粮仓模型,EDEM离散元仿真分析大豆颗粒分层充填与多向放粮过程中流动形态演变及侧向压力分布规律。研究发现静态储粮时角约束区力链强度较自由区高1.2倍,底部静压力达顶部的1.8倍,动态压力峰值出现在充粮中期及过渡带,94%监测点超国标GB50077-2017值。提出修正的侧向压力计算公式,为多区隔粮仓设计提供理论依据。

广告
   X   

刘海燕|曲洪波|卢学飞|郭金勇|何贤梅|李颖|范晓轩|徐向南
河南工业大学土木工程学院,中国郑州450001

摘要

为了实现分类、基于质量的粮食储存,粮食行业提出了一种新型筒仓设计,在传统的圆形筒仓内增加了隔板。添加这些隔板后,筒仓的横截面从圆形变为扇形,与传统中心排放式筒仓相比,其排放流型和侧向压力特性有所不同。本研究使用SolidWorks建立了一个高1000毫米、半径250毫米的扇形筒仓模型,并利用EDEM软件进行了三次平行填充和排放模拟。通过分层染色和力链分析等方法研究了流动模式和侧向压力分布。结果表明,在静态储存过程中,筒仓下部的力链强度约为上部的1.8倍,而受角度限制的区域则约为不受限制区域的1.2倍。这导致筒仓中下部出现峰值侧向压力,靠近角度限制区域的压力更高。排放过程分为三个阶段:“质量流–混合流–漏斗流”,高度分布呈现“内部低、外部高”的特点,受角度限制的影响。在排放过程中,力链会周期性形成和崩塌。约70%的测量点在质量流阶段显示出峰值动态压力,最大值集中在流动过渡深度附近。94%的点的平均动态压力超过了GB50077-2017标准的数值,凸显了多仓筒仓在标准制定方面的局限性。所有点都出现了超压现象,中部和上部的超压更为严重。本文提供了改进的侧向压力公式,为扇形筒仓设计提供了理论参考。

引言

筒仓作为一种具有优异储存能力的基础设施,在农业、化工、建筑、采矿和航运等多个行业中发挥着关键作用,连接生产、加工和运输环节。其核心功能是储存和转运散装物料[1]、[2]、[3]。在农业领域,筒仓常用于储存小麦、玉米、水稻、大豆等粮食作物及其副产品[4]、[5]、[6]。近年来,对不同等级粮食作物进行分类储存的需求日益增长,对筒仓的储存质量提出了更高要求。早在1961年,Shlifer就提出了一种将仓库划分为多个仓室的方法来储存不同类型的物料[7]。这种类型的筒仓被称为多仓筒仓,不仅有效节约了土地资源,还显著提高了经济效益,实现了“一个筒仓多重储存”的目标[8]。在筒仓内安装隔板以实现分区储存,不仅大大提高了储存容量的利用率,还为粮食作物的分类储存提供了更大的灵活性。添加隔板后,筒仓的横截面几何形状从圆形变为扇形,使得各个仓室更加细长。这种变化带来了诸如不可预测的排放流型和波动的侧向压力等挑战[9]。此外,现有的筒仓设计规范[10]对这些问题的理论基础较为缺乏。
目前,全球范围内已进行了大量关于筒仓内储存物料流动行为及装卸过程中壁压响应的实验和数值研究。研究对象包括各种结构形式(圆形、矩形、筒仓群、带有内部组件的筒仓等)和不同的储存物料(玉米、小麦、水稻、水泥原料、标准砂、玻璃颗粒、煤球等)。Khalil等人[11]系统回顾了数十年的筒仓相关研究,涵盖了结构行为[2]、[12]、颗粒-材料-结构相互作用[13]、[14]、[15]、[16]以及失效模式和原因[17]、[18]。总体而言,关于储存物料流动状态与筒仓壁压响应之间关系的现有研究成果可分为三类:(1)基于填充/排放试验的压力分布模式及其影响因素的研究。全尺寸和缩比试验表明,筒仓结构几何形状、料斗角度、填充方法和排放口尺寸显著影响流动模式和压力分布。特别是在料斗顶部和过渡区,容易发生明显的压力波动和集中现象[9]、[20]、[21]、[22]、[23]、[24]、[25]。(2)关于排放动态效应和超压机制的研究。多项研究表明,排放过程中存在明显的动态侧向压力放大和波动特征,峰值远超理论预测。其背后的机制通常与流动模式转变、自由落体拱形效应、拱形-卸载循环以及力链重构过程密切相关[26]、[27]、[28]、[29]。(3)数值模拟和耦合方法研究:离散元方法(DEM)和有限元方法(FEM)被广泛用于揭示颗粒流动和结构响应:DEM能有效揭示质量流和漏斗流等流动模式以及筒仓壁压力的变化[27]、[28]、[29]、[30],而FEM可以研究柔性壁筒仓在地震载荷下的非均匀侧向压力分布和动态超压[3]、[32]、[33]。将DEM和FEM相结合进一步提供了评估结构-颗粒相互作用的方法[34]。
总之,现有的筒仓压力研究主要集中在具有均匀横截面的筒仓(如圆形或矩形横截面)上,实验或数值模拟研究中使用的储存物料多为颗粒状材料(如粮食、砂或轻质球)。关于多仓不规则筒仓的研究主要使用水泥作为储存物料。然而,这些研究仅限于实验阶段,并受测试方法的限制,难以全面反映储存物料之间的相互作用、分区条件下的流动模式变化以及侧向压力分布情况。对于用于储存大豆等散装颗粒物料的扇形筒仓,目前尚未发现相关文献使用离散元方法结合分层染色和力链分析来研究流动模式和筒仓壁压力分布。本研究以大豆作为储存物料,使用SolidWorks建模软件构建了一个高1000毫米、半径250毫米的扇形仓模型。选择半径为3.45毫米、球形度为95.47%、粒径服从标准正态分布(标准差0.05毫米)的大豆颗粒作为储存物料。利用EDEM离散元分析软件进行了三次平行填充/排放模拟。储存物料被分层并标记成100毫米厚的层。最后,通过整合力链网络分析等方法,系统研究了储存物料的流动模式和筒仓壁上的压力分布。

部分摘录

粮食筒仓和储存物料模型

由于EDEM预处理模块的建模能力有限,本研究基于我们团队自行设计的丙烯酸筒仓模型[9],使用SolidWorks构建了一个三仓粮食筒仓模型。该模型保存为STEP格式,并导入离散元分析软件EDEM中。导入后,定义了筒仓壁、隔板、筒仓底部和挡板的三维几何模型作为颗粒接触和碰撞的固定边界。

理论研究

关于筒仓壁的侧向压力,许多学者进行了理论推导,代表性的模型包括Janssen理论[39]、Airy理论[40]和Reimbert理论[41]。目前各国现有的筒仓设计规范[10]、[42]、[43]均基于Janssen理论,并根据各国筒仓设计和应用经验进行了修改。
根据Janssen理论,储存物料对单位面积施加的水平压力

h0(kPa)为...

流动模式

为了直观观察筒仓内储存物料的流动模式,物料从上到下以0.1米的间隔分层标记,共10层,如图7a所示。从上到下分别标记为F1至F10层。筒仓壁表面(弧面a-c)和隔板表面(倒三角形表面a-b-c]分别被指定为观察表面A和B(图7b)。为了观察排放口的流动状态变化...

结论

本研究使用SolidWorks构建了一个高1000毫米、半径250毫米的扇形筒仓模型。选择半径为3.45毫米、粒径服从标准正态分布(标准差0.05毫米)的大豆颗粒作为储存物料。利用EDEM离散元分析软件进行了三次平行填充和排放模拟。储存物料被分层并标记成100毫米厚的层。最后,通过...

作者贡献声明

刘海燕:撰写 – 审稿与编辑,撰写 – 初稿,监督,方法论,正式分析,概念化。曲洪波:撰写 – 审稿与编辑,撰写 – 初稿,方法论,正式分析,数据整理,概念化。卢学飞:监督,资源获取,调查,资金申请。郭金勇:软件应用,方法论,数据整理。何贤梅:软件应用,方法论,正式分析,数据整理。李颖:方法论,正式分析,数据整理。

利益冲突声明

作者声明没有可能影响本文研究的财务利益或个人关系。

致谢

本研究得到了河南省粮油储存设施与安全重点实验室(编号:2024KF04)、2025年河南省重点科学技术研究项目(编号:252102111170)、2025年河南省重点科学技术研究项目(编号:252102111181)和2024年河南省重点科学技术研究项目(编号:242102110352)的支持。

生物通微信公众号
微信
新浪微博


生物通 版权所有