料斗集沙过程与主要参数的离散元仿真研究_郝腾达.pdf

ISSN 1008-9446CN13-1265/TE承 德 石 油 高 等 专 科 学 校 学 报Journal of Chengde Petroleum College第 25 卷第 2 期,2023 年 4 月Vol.25,No.2,Apr.2023料斗集沙过程与主要参数的离散元仿真研究郝腾达a,高占凤b,郑明军a,吴文江c(石家庄铁道大学a.机械工程学院;b.工程训练中心;c.教务处,河北石家庄 050043)摘要:为了探究料斗与沙粒作用特性及转速、斗张角和除沙深度等主要参数对集沙量的影响,借助离散元法,在 EDEM 软件中创建料斗沙粒的离散元模型,对料斗集沙过程沙粒的运动特性进行仿真分析,阐述了料斗除沙过程的形成机理,为料斗强度提供了校核依据,并使用单一因素法对料斗的主要集沙参数进行探究。仿真结果表明料斗转速、斗张角和除沙深度对料斗集沙量有明显影响。当料斗转速为 25 r/min、斗张角 120和除沙深度 110 mm 时,料斗集沙量达到最大。研究结果对提高料斗的除沙效率具有实际意义,为后续集沙装置优化设计提供一定参考基础。关键词:料斗;离散元;除沙车;EDEM中图分类号:TH225文献标志码:A文章编号:1008-9446(2023)02-0057-06Discrete Element Simulation Study of Sand Collection Process andMain Parameters of HopperHAO Teng-daa,GAO Zhan-fengb,ZHENG Ming-juna,WU Wen-jiangc(a.School of Mechanical Engineering;b.Engineering Training Center;c.Office of Academic Affairs,Shijiazhuang Tiedao University,Shijiazhuang 050043,Hebei,China)Abstract:In order to explore the effect of interactive characteristics of hopper,sand particles andinfluence of main parameters such as speed,opening angle of the bucket and depth of sand removalon the amount of sand collection,using the discrete element method,a hopper-sand discrete elementmodel is created in EDEM software,movement characteristics of sand process for hopper set is ana-lyzed in the simulation,and the hopper load resistance figure is obtained.This paper expounds theformation mechanism of hopper bed cleaning process,and provides check basis for hopper strength,besides,the main sand collection parameters of the hopper is explored by using single factor method.Studies show that hopper speed bucket angle and depth of hopper bed cleaning set have significanteffects on the sediment.When the hopper speed reaches 25 r/min,the bucket angle is 120and bedcleaning depth is 110 mm,the sand hopper sets maximum amount.The study results have practicalsignificance for improving bed cleaning efficiency of hopper,and provide a certain reference basis foroptimization of subsequent sand device design.Key words:hopper;discrete element;sand truck;EDEM基金项目:河北省自然科学基金(基于 DEM 的机械结构与低粘结度沙土相互作用机理研究):E2017210166收稿日期:2021-12-07第一作者简介:郝腾达(1997-),男,河北石家庄人,硕士研究生,主要从事机械设计及理论研究,E-mail:243594432 。伴随着我国铁路事业的蓬勃发展,铁路线的分布范围也愈发广泛,其中不乏涉及到条件恶劣的铁路线,比如风沙频繁的西北地区等,常年的风沙使当地铁路线埋没沙土中,以致列车难以顺利通行。为了解决上述问题,石家庄铁道大学研制出料斗式除沙车,是用于清理轨道积沙的的一种铁路维护设备。为了使集沙料斗具备良好的除沙效果,探究料斗除沙机理效应,对作业过程中料斗沙粒作用的分析必不DOI:10.13377/ki.jcpc.2023.02.016承德石油高等专科学校学报2023 年第 25 卷第 2 期可少。由于沙粒的运动具有随机性,发散性等特点,传统的研究方法很难对沙粒的运动特性进行准确地分析,离散元法成为了探究料斗除沙运动的运动机理的主要手段,它是一种研究散体和机械相互接触的有效方法,EDEM 软件是应用离散元技术来模拟研究颗粒物料系统的 CAE 软件1-2。当前,我国已有学者使用离散元法分析了沙粒与机械之间的作用,如文献3等分析了不同集沙铲转速下沙粒能量的分布规律和不同集沙铲半径的受力情况;文献4等分析了沙粒在除沙车内的数值规律及最佳除沙性能,证明了除沙过程离散元仿真可行性。本文采用离散元法对料斗除沙运动的运动机理进行探究。本文以离散元法为基础,对料斗的集沙过程进行 DEM 仿真分析,探究沙粒的运动特性,了解沙粒与料斗的作用机理。分析比较料斗在不同转速、斗张角和除沙深度对集沙量的影响,选取最佳的工况参数以实现料斗集沙量的最大装载量。1离散元法理论颗粒离散元法是一种求解与分析离散颗粒系统动力学问题的数值计算方法5,在每步时长内,离散元法能根据各单元之间所选取接触模型的接触作用和牛顿第二定律,利用相关迭代方法计算出各个单元的受力和位移等状态,并对所有颗粒单位位置进行更新,从而得知整体的宏观运动规律。牛顿第二定律主要计算颗粒的位移、速度和加速度等,通过位移-力的关系能够从颗粒的位移得知颗粒单元受到的作用力,颗粒 i 的位移可由牛顿第二定律得出,然后计算出颗粒 i 受力情况。其求解计算原理如图 1 所示。颗粒接触模型能很好地反映颗粒的运动状态,选取合适的接触模型至关重要,是离散元法计算颗粒间位移和作用力的基础。由于铁路上的沙粒的黏连性不强,每个沙粒都具有独立性且沙粒之间的接触都是弹性碰撞,而 Hertz-Mindlin(no slip)接触模型具备高效运行效率、精准计算能力和弹性碰撞等特性,所以本文选用 Hertz-Mindlin(no slip)接触模型在离散元中进行运动分析,Hertz-Mindlin(no slip)接触模型示意图如图 2 所示。其中,R1 和 R2,Z1 和 Z2 分别为两颗粒的尺寸半径和接触半径,为接触圆半径,为接触变形量。其法向接触力为:Fn=kn3/2n+cnn(1)其切向接触力为:F=min-S-256SmFn|(2)式中,kn为法向刚度;n为法向重叠量;n为法向相对速度影响因子;cn是法向阻尼系数;S为切向刚度;为切向重叠量;为切向相对速度;是滑动摩擦系数。2集沙机构与仿真模型的建立2.1除沙车集沙原理与模型简化料斗式轨道除沙车主要由支撑装置、集沙装置、排沙装置、连接装置、驱动装置和行走装置组成,整85郝腾达,等:料斗集沙过程与主要参数的离散元仿真研究体结构如图 3 所示。集沙装置为主要部件,它由多个料斗和回转支承组成。集沙时,回转支承内圈固定在支撑装置上,外圈带动料斗进行转动,料斗通过旋转使自身与沙粒发生剪切作用进行铲沙运动,起到集沙目的。排沙装置是借助输送带把沙粒排出轨道外,支承装置用来固定集排沙装置,保证稳定工作,连接装置用来连接集排沙装置和驱动装置,驱动装置采用液压驱动方式驱使除沙车前进,行走装置是对底盘通过安装履带方式使除沙车能顺利移动。在 EDEM 中,需要创建沙粒与集沙部件模型进行集沙模拟仿真,为了缩减计算时间,对集沙装置仿真模型进行简化处理,只保留单个料斗和回转支承;单轨两侧均设有道钉,据测量道钉顶部距轨面的距离最小值为 75 mm,在轨面和道钉留有一定间隙的情况下,允许料斗清理的轨面以下深度最大可达70 mm,仿真模型如图 4 所示。2.2仿真运动定义2.2.1定义仿真参数风沙地区沙粒大小主要由 0.1 0.5 mm 中细沙组成,如果按照实际颗粒进行仿真会大大降低工作效率,基于现有工作站配置,使用粒径为 6 mm 的球形颗粒进行仿真。在进行模拟仿真之前,需要对沙粒与料斗的材料参数及接触参数进行正确标定,以确保其仿真时真实性。为了创建合适的风沙地区沙粒模型,采用实验安息角和模型安息角相结合的方法对标定参数加以验证,在 Fielke 建立的沙质土壤模型6基础上对沙土进行实验和模型标定,如图 5 所示。经测量,颗粒安息角为 34(2),实验安息角和模拟安息角吻合,可用于离散元仿真,通过实验和虚拟标定确定颗粒和几何体的物理参数如表 1 所示,颗粒与几何体之间的接触参数如表 2 所示。表 1沙粒与料斗材料性能属性材料泊松比剪切模量/(108Pa)密度/(kgm-3)沙粒0.31.02 650料斗0.37007 800表 2沙粒与料斗接触参数材料静摩擦因数滚动摩擦因数恢复系数沙粒沙粒0.530.080.6料斗沙粒0.450.050.695承德石油高等专科学校学报2023 年第 25 卷第 2 期2.2.2定义颗粒工厂由于生成的沙粒道床比较大,所以本文先生成缩小整数倍的小型沙粒床,小型沙粒床大小为原颗粒床的十分之一,在利用 EDEM2018 的快速填充效果生成仿真颗粒道床;设定产生沙粒 2 kg,生成速率4 kg/s;时间步长过大会导致运动失真,过小则会增加计算量,设置时间步长为 Rayleigh 时间的 25%;网格大小划分为最小颗粒半径 3 倍;仿真总时间设置为 1 s。2.2.3定义运动属性仿真前需对模型添加运动属性,由于料斗做旋转运动,赋予料斗线性转动(Linear Rotation)属性,其大小为现有料斗式除沙车料斗转速 20 r/min。3仿真结果与分析沙粒运动特性分析是计算集沙阻力的前提,沙粒运动特性分析能探究沙粒与料斗的作用机理,了解料斗集沙的过程,为研究料斗集沙量奠定一定基础。根据沙粒模型的运动状态可得知沙粒对料斗产生的集沙阻力。沙粒运动特性如图 6 所示,以现有料斗转速 20 r/min 为例,选取 0.1、0.25、0.5、0.8 s 时间节点的仿真图进行观察与分析,其中红色速度最大,蓝色速度最小。由图 6 可知,除沙料斗下部的沙粒受到的干扰比较小,沙粒颜色基本无变化呈速度最小的蓝色状态,扰动沙粒大多集中在料斗周围,并且随着旋转角度的增加,除沙料斗的受力逐渐增大。料斗在铲入沙粒堆时,沙粒堆积面受到一个剪切作用,使得与料斗前板直接接触的沙粒所受扰动较大,出现个别红色颗粒,速度最大为 2.166 7 m/s,由于沙粒之间的挤压运动使得与底部接触的沙粒也处于加速状态;随着料斗的继续转动,料斗