基于DEM-MBD的高压辊磨机仿真模型及应用研究.pdf

Series No.567September 2023 金 属 矿 山METAL MINE 总 第567 期2023 年第 9 期收稿日期 2023-03-28基金项目 东北大学航空动力装备振动及控制教育部重点实验室研究基金资助项目(编号:VCAME202102)。作者简介 任云鹏(1974),男,教授,博士,硕士研究生导师。基于 DEM-MBD 的高压辊磨机仿真模型及应用研究任云鹏1 李安帅1 李旭东2 宋 芳1(1.沈阳建筑大学机械工程学院,辽宁 沈阳 110168;2.山东黄金矿业(莱州)有限公司三山岛金矿,山东 莱州 261442)摘 要 针对常用数学模型和目前离散元仿真模型在预测高压辊磨机处理量、单位能耗和破碎比等方面误差较大的问题,提出采用离散单元法(DEM)和多体动力学(MBD)联合仿真技术对高压辊磨机进行仿真建模分析。以山东黄金矿业(莱州)有限公司三山岛金矿正在使用的 KHD3751400-1000 型高压辊磨机为研究对象,根据相关理论和实际设备运行参数进行关键仿真参数设定,利用 EDEM 和 MotionView 进行联合仿真,并将仿真结果与数学模型、单一EDEM 仿真模型和矿山实际运行数据进行对比。结果表明:DEM-MBD 联合仿真模型能够更真实、更准确地模拟高压辊磨机的运行状态,在高压辊磨机预测性能指标上有更好的效果,DEM-MBD 联合仿真模型比 DEM 仿真模型在处理量、单位能耗、破碎比上的误差值分别减小了 13.41、38.10、6.74 个百分点,有效降低了现有仿真模型的误差值。该研究也为此类复杂破碎设备的仿真提供了合理参考,具有一定的工程应用价值。关键词 高压辊磨机 离散单元法 多体动力学 联合仿真 模型对比 中图分类号TD453 文献标志码A 文章编号1001-1250(2023)-09-164-09DOI 10.19614/ki.jsks.202309022Simulation Model and Application Research of High-Pressure Grinding Rolls Based on DEM-MBDREN Yunpeng1 LI Anshuai1 LI Xudong2 SONG Fang1(1.School of Mechanical Engineering,Shenyang Jianzhu University,Shenyang 110168,China;2.Sanshandao Gold Mine,Shandong Gold Mining(Laizhou)Co.,Ltd.,Laizhou 261442,China)Abstract Aiming at the problems of the commonly used mathematical models and the current discrete element simula-tion models in predicting the processing capacity,unit energy consumption and crushing ratio of high-pressure grinding rolls(HGPR)with large errors,this paper proposes the simulation modeling analysis of HGPR using the co-simulation technology of discrete element method(DEM)and multi-body dynamics(MBD).Taking the KHD3751400-1000 HGPR used in Sans-handao Gold Mine,Shandong Gold Mining(Laizhou)Co.,Ltd.,as the research object,the key simulation parameters are set according to the relevant theory and actual equipment operation parameters,and the co-simulation is carried out by EDEM and MotionView.The simulation results are compared with the mathematical model,single EDEM simulation model and the actual operation data of the mine.The results showed that the DEM-MBD co-simulation model can simulate the operation of the HGPR more realistically and accurately.Compared with the DEM simulation model,the DEM-MBD co-simulation model reduces the error values of processing capacity,unit energy consumption and crushing ratio by 13.41,38.10 and 6.74 percentage points re-spectively,which effectively reduces the error values of the existing simulation model.The study also provides a reasonable ref-erence for the simulation of such complex crushing equipment,which has some engineering application value.Keywords high-pressure grinding rolls,discrete element method,multi-body dynamics,co-simulation,model comparison 高压辊磨机是一种新型破磨粉碎设备,具有处理能力强、单位破碎能耗低、设备工作效率高等特征1。该设备最先应用于水泥领域,随着技术的不断发展,已经广泛应用于金属矿山领域2-3。高压辊磨机的相关研究涉及机械运动学、动力学、离散元等诸多方面,当前研发设计中仍存在现有模型不准确、破碎过程中重要参数无法定量表示等问题。虚拟仿真技术是改进粉碎设备设计和操作的一项重要方法,它能够比较不同参数设置而不需要物理原型或生产损失的相关成本,其仿真结果的可用性取决于模型的保真度。BILGILI 等4用严格的数学模型研究了多个接触颗粒的粉碎,得出了粉碎过程矩461阵,但无法考虑到设备参数、颗粒材料性质对粉碎效果的影响。TORRES 等5以矿石性质、高压辊磨机工作参数、粉碎工艺为基础,研究了一系列高压辊磨机粉碎的数学模型,包括处理能力模型、能耗模型,但在预测准确度上有所不足,无法应用到工业领域。DANIEL 等6通过大量实验开展了关于高压辊磨机数学模型及应用的专项研究,验证和扩展了高压辊磨机的产品粒度分布、处理量以及能量利用率的数学模型。KUMAR7使用 EDEM 模拟了实验室规格的高压辊磨机,研究了颗粒破碎过程,以及工艺参数和操作参数对处理量和能耗的影响。祖大磊等8和何剑伟等9根据高压辊磨机破碎过程中的工艺参数、矿石单轴抗压强度、给料粒度及粉碎产品粒度分布对比能耗的影响,提出了利用高压辊磨机比能耗进行计算评估的 DUCS 数学模型。鲍诺10利用 EDEM 软件研究了高压辊磨机工作过程中颗粒运动状态与工作参数之间的关系,得到了不同辊速和粒径下颗粒所受压力随时间变化的规律。张乐11利用 EDEM 软件模拟了高压辊磨机破碎矿石过程,通过正交试验方法和单因素试验方法研究压辊转速、压辊间隙以及压辊与矿石间的静摩擦因数对矿石破碎效果的影响,但单一的DEM 仿真模型忽视了设备的动态响应,无法真实模拟实际工作状态。在高压辊磨机中,压辊间隙不是一个直接控制的变量,是所施加的压力和颗粒床之间相互作用的结果。以上基于理论的数学模型和基于离散元法的仿真模型由于忽略了颗粒对高压辊磨机运行的影响,无法描述辊隙的动态变化对仿真结果的影响,导致在预测处理量、单位能耗和破碎比时均存在较大误差。针对高压辊磨机现有模型存在误差和无法描述设备的动态辊隙随颗粒载荷变化规律的问题,为了减小误差,降低企业研发设计成本,在离散元法基础上,提出采用离散元(DEM)和多体动力学(MBD)联合仿真技术对高压辊磨机进行仿真分析,通过对比数学模型、单一离散元仿真模型以及实际工作数据,证实联合仿真模型有着较为理想的结果。1 高压辊磨机数学模型1.1 动态辊隙模型动态辊隙模型是在动辊动力学模型基础上推导而来的,高压辊磨机主要由两个压辊和液压系统组成,液压系统部分蓄能器中的压缩氮气充当液压系统的弹簧。在破碎过程中,动辊和定辊的旋转速度是恒定的,但动辊会因受到物料的挤压力和液压系统提供的力往复平移。动辊的横向移动可以由弹簧阻尼系统来描述,如图 1 所示。图 1 高压辊磨机动态辊隙模型Fig.1 Dynamic roll gap model of high-pressure grinding rolls高压辊磨机的辊隙动态变化可由下式表示:xt()=F t()+cxt()+kx t()m,(1)xt+t()=xt()+xt()t,x t+t()=x t()+xt+t()t,(2)S t()=x t()+s,(3)式中,x t(),xt(),xt(),x t+t(),xt+t(),xt+t()分别是当前时间和下一时间动辊的位移,m;速度,m/s;加速度,m/s2;F(t)是当前时间物料作用在动辊上的总力,N;t 是时间间隔,s;m 为动辊质量,kg;c 为阻尼系数,N/(m/s);k 为弹簧刚度,N/m;S(t)为当前时间辊隙大小,m;s 为静态辊隙大小,m。1.2 处理量模型高压辊磨机的处理量数学模型可以基于设备的设计参数和物料的物理性质进行推导。最常见的处理量数学模型是基于高压辊磨机的辊隙宽度、压辊直径、压辊旋转速度、压辊面长度、料饼密度等参数来进行建模的。图 2 是高压辊磨机的粉碎示意图。图 2 高压辊磨机粉碎示意Fig.2 High-pressure grinding rolls crushing diagram根据质量守恒定律,物料在挤压区域内,质量不发生改变,并且假设辊面与物料之间不存在相对滑动,则高压辊磨机的单位处理量 G 可以表示为:G()=3 600()S()Lvcos,(4)式中,G()为任意圆心角 下高压辊磨机的单位处理量,t/s;()为任意圆心角 下物料密度,t/m3;S()为任意圆心角 下的压辊间隙,m;L 为辊宽,m;v 为压辊线速度,m/s。在实际应用中,高压辊磨机的处理量一般指出料口位置处,即 =0 时,高压辊磨机的理论处理量 G为12:G=3 600sLv.(5)561 任云鹏等:基于 DEM-MBD 的高压辊磨机仿真模型及应用研究 2023 年第 9 期式中,为料饼密度,一般为物料密度的 0.8513,t/m3;s 为出料口间隙,m。1.3 单位能耗模型图 3 是高压辊磨机压辊受力示意图。图 3 高压辊磨机压辊受力示意Fig.3 Force diagram of pressure roller of high-pressure grinding rolls如图所示,阴影区域是工作压力的影响区域,在该区域内压辊对物料的挤压力 F 可以通过压辊的工作压力求得。由于高压辊磨机的工作特点,物料只在压缩区受到压辊的挤压力,故相应的受力区域为压辊的上半部分,因而压辊的挤压力 F 为:F=1 000PD2L,(6)式中,F 为压辊对物料的挤压力,kN;P 为压辊的工作压力,MPa;D 为辊径,m。则垂直方向上的挤压力 F 引起的转矩 T 为:T=Fsin()D2,(7)式中,T 为转矩,kN;为发生层压粉碎角度的一半。高压辊磨机具有 2 个压辊,那么总功率是压辊垂直方向上挤压力 F 和压辊线速度 v 乘积的 2 倍,故功率 P 为:P=2Fsin()v,(8)单位能耗为总功率 P 与处理量 G 的比值,故高压辊磨机的理论单位能耗 W 为14:W=PG,(9)式中,W 为理论单位能耗,kW h/t;G 为处理量,t/h。1.4 破碎比模型平均破碎比是破碎机的一个重要性能参数,用来衡量高压辊磨机对物料的破碎效果。平均破碎比 i通常定义为物料经过高压辊磨机破碎后,进料的平均粒径与出料的平均粒径之比,该种方法求得的破碎比能较为真实地反映破碎程度。将颗粒划分为体积或质量相等的两部分的粒径 d50亦被称为中位径或中值粒径,常用来表征颗粒样品的平均粒径15。高压辊磨机出料的平均粒径可以用 d50表示,d50可通过产品粒径分布函数罗逊拉姆勒(Rosin-Rammler)求得。该函数是一种经验公式,可以用来描述颗粒大小的分布情况。该函数的形式如下16:(d)=1-exp-c(ddmax-d)n,0 d dmax,(10)式中,(d)为累计粒度产率分布函数;d 为颗粒粒径,mm;c,n,dmax是根据经验数据确定的参数,其中dmax为破碎产品中最大粒径,mm。则高压辊磨机的理论平均破碎比 i 为:i=D50d50,(11)式中,D50和 d50分别表示破碎前后累计粒度分布百分数达到 50%时所对应的粒径。2 DEM-MBD 联合仿真模型建立考虑到高压辊磨机是一个复杂的多体动力学系统,其运行状态会受到多个因素的影响,当物料与压辊之间的挤压作用力导致物料材料性质发生变化时,如粉碎过程中的物料磨损和破碎,则必须考虑设备的动态响应,如辊隙的动态变化,因为它会影响使用DEM 进行模拟仿真的效果。对于单一 DEM 仿真模型,只能考虑压辊固定间隙,无法反映辊隙的动态变化对系统的影响。而辊隙的变化又是影响高压辊磨机处理量、单位能耗、破碎比的关键因素,所以现有的DEM 仿真模型在预测高压辊磨机的性能上具有较大误差。高压辊磨机的 MBD 仿真需要对压辊施加瞬态载荷,以便描述辊隙随载荷变化而变化的特性。DEM-MBD 联合仿真利用两软件之间的数据交互功能,将物料破碎过程中引起的设备实时动态响应又反馈给破碎过程。DEM-MBD 联合仿真模型引入了辊隙的动态变化,可以更真实、更准确地模拟高压辊磨机的运行状态和性能,从而提高仿真结果的可靠性和有效性,故本文提出采用 DEM-MBD 联合仿真方法预测高压辊磨机性能指标。2.1 联合仿真方法高压辊磨机联合仿真的流程图见图 4。当仿真开始后,MBD 部分将定义完材料参数和运动参数的动力学模型传递给 DEM 部分,DEM 部分颗粒工厂产生颗粒,然后由 EDEM 软件计算确定每个颗粒的位置和速度,以及颗粒之间和颗粒与几何体之间的作用力。DEM 部分结束后,EDEM 将几何体上的等效粒子力和力矩再反馈给 MBD 部分,MBD 部分由 Motion-view 软件实现求解 MBD 的控制方程,并计算几何体的位置和速度。2.2 多体动力学模型高压辊磨机几何机构的精确构建是联合仿真的关键,以 KHD3751400-1000 型高压辊磨机为研究对661总第 567 期 金 属 矿 山 2023 年第 9 期图 4 DEM-MBD 联合仿真流程图Fig.4 DEM-MBD co-simulation flow chart象,利用 SolidWorks 软件构建几何机构模型,包括定辊、动辊、机架、液压缸、氮气蓄能器、喂料装置六大部分,如图 5 所示。图 5 高压辊磨机几何模型Fig.5 Geometry model of high-pressure grinding rolls2.2.1 约束关系构建将构建好的 KHD3751400-1000 型高压辊磨机几何模型导入 Motionview,此时模型之间只存在空间位置关系而不存在运动学上的约束关系。因此为保证虚拟高压辊磨机的准确运动,首先应根据各零部件的实际运动关系构建模型之间的约束关系。其约束关系如表 1 所示。表 1 约束关系Table 1 Binding relationships名称约束固定副机架-地固定副机架-喂料装置固定副定辊座-机架旋转副定辊-定辊座旋转副动辊-动辊座移动副动辊座-机架 高压辊磨机零部件的几何关系和约束关系构建完后,需要进一步施加作用在压辊上的旋转运动,从而使高压辊磨机能够实现破碎作业。2.2.2 液压系统参数理论计算高压辊磨机中的液压系统由 3 个主要部件组成:液压缸、氮气蓄能器和液压油。蓄能器中的压缩氮气充当液压系统的弹簧,其刚度可通过初始氮气压力进行调整,初始油压和氮气压力的组合有效地决定了弹簧的非线性刚度响应。图 6 为液压系统力学模型图,当力从 f1增加到 f2时,它会引起压辊位移 s,该位移被传递到横截面积为 ap的活塞上,使氮气蓄能器中的压力从 p1增加到 p2。图 6 高压辊磨机液压系统力学模型Fig.6 Mechanical model of hydraulic system of high-pressure grinding rolls液压系统部件参数与弹簧刚度 k 和阻尼系数 c有关,可通过以下公式估算17:k=nf2p0V0,(12)c=kgf,(13)f=npapp1,(14)式中,n 是多变指数;p0是初始氮气压力,MPa;V0是初始氮气体积,m3;p1是初始液压,MPa;ap是液压缸的横截面积,m2;np是气缸数,个;g 是重力加速度,m/s2。该模型考虑了液压系统的主要部件,但未考虑其他部件,如阀门、油管、油黏度等。根据实际高压辊磨机运行参数,依据上述公式计算出液压系统的弹簧刚度 k 和阻尼系数 c,通过在Motionview 中添加弹簧系统部件并设置弹簧刚度 k和阻尼系数 c,实现液压系统的代替。2.3 离散元模型2.3.1 颗粒破碎模型的选用目前 EDEM 中有 2 种用于模拟颗粒破碎的模型:Bonding 键模型和 Tavares 破碎能模型。Bonding键模型是将小颗粒通过 Bonding 键粘结成大颗粒,以键的断裂来模拟颗粒的破碎。Bonding 键模型在颗粒形状和大小相对均匀的情况下表现较好,但在颗粒形状和大小差异较大时,可能会出现较大误差,并且计算量比较大,破碎后无法对颗粒粒径进行统计。761 任云鹏等:基于 DEM-MBD 的高压辊磨机仿真模型及应用研究 2023 年第 9 期Tavares 破碎能模型是一种基于经验的破碎模型,当冲击能大于破碎能时,颗粒发生破碎,没有破碎的颗粒也会受到损伤,产生更小的破碎能。该模型基于经典的 Kick、Rittinger 和 Bond 破碎理论,考虑了影响颗粒内部的应力分布和破碎的各种因素,如颗粒的强度、形状、大小、载荷等18。Tavares 破碎模型适用于各种形状和大小的颗粒,且在颗粒破碎前后的形状不会发生太大变化的情况下表现较好,并且计算速度快,可对破碎后颗粒进行粒径分布、获取破碎比等后处理,能直观展示破碎效果。高压辊磨机仿真破碎过程颗粒数量巨大,且需要分析颗粒破碎后的破碎比、处理量和单位能耗等。相比 Bonding 键破碎模型,Tavares 破碎模型能够方便对颗粒粒径进行统计,并且能同时满足其他的仿真需求。Tavares 破碎模型原理如图 7 所示。图 7 Tavares 破碎原理Fig.7 Tavares crushing principle在 Tavares 破碎模型中,每个颗粒都有特定的破碎能量,根据其大小、平均值和标准偏差进行分配。这个能量根据公式19-20所描述的分布而变化。P(E)=121-erflnE-lnE5022E(),(15)E=EmaxEEmax-E,(16)式中,E 为破碎能量分布,J/kg;Emax为能量分布的上限值,J/kg;E50为能量分布的中值,J/kg;E为标准差。破碎能量的中值由以下公式给出:E50=E1+KP/Kst1+d0dp(),(17)式中,E为极限破碎能,J/kg;是根据试验数据调整的参数;KP和 Kst分别为颗粒和钢的硬度,GPa;dp为颗粒的大小,mm;d0为矿石中特征颗粒的大小,mm。当颗粒受到挤压而未发生破碎时,它将受到损伤,之后颗粒会产生新的破碎能 Ef,新的破碎能会低于其原破碎能 Ef:Ef=Ef(1-D),(18)D=22-5D+5eEkEf25,(19)式中,Ef为颗粒的破碎能,J;eEk为有效碰撞能量,J;D 为损伤系数;为损伤累积系数。Tavares 破碎模型主要是基于参数 t10,t10表示产品中小于母粒径 1/10 的碎片比例,由公式描述为:t10=A 1-exp-beEkEf(),(20)式中,A 和 b 是根据试验数据调整的参数。t10值越大,新产生的颗粒就越细。根据 Tavares 团队研究成果,参考 RODRIGUEZ等21校核的 Tavares 破碎模型参数来定义物料的破碎参数,具体破碎参数如表 2 所示。表 2 Tavares 破碎模型参数Table 2 Tavares crushing model parameters参数数值特性损伤累积系数 5计算物料损伤的常数极限破碎能 E/(J/kg)500E50参数,极限破碎能特征颗粒尺寸 d0/mm75E50参数,裂纹萌生到裂纹扩展的过渡物料粒径E50拟合参数 0.35E50参数拟合参数,用于计算物料破碎能标准差 0.8破碎能分布的标准差,计算物料破碎的概率t10参数 A/%68单次冲击破碎中 t10最大值,百分比t10参数 b0.03t10拟合参数Dmin/mm1.2用户自定义的破碎中允许生成的最小粒径Emin/J110-4纳入统计的最小冲击能量,推荐110-4Ct0切向冲击能百分比Truncation_ratio/%100破碎能对数正态分布的上截断2.3.2 仿真时间步长理论计算在 EDEM 仿真中首先需要合理地设定时间步长,时间步长的设定需要考虑多个因素,如模拟的物理现象、计算机的性能、仿真时间等等。时间步长的大小对于仿真结果的准确性和计算速度都有很大的影响,如果设置时间步长过大,则会令仿真结果不收敛,并且对于硬球颗粒,过大的时间步长可能会使某些颗粒穿过设备,使得仿真结果不准确。如果设置的时间步长太小,则会要求计算机较高的配置并且耗费时间较长。因此,针对不同的仿真工况需要选取适当的时间步长。颗粒在发生碰撞时,百分之七十的能量是通过瑞利波消耗的,因此临界的仿真时间步长需要根据瑞利波的传播速度决定。瑞利波是指当两个单元发生碰撞时,因单元表面受到变应力而产生的沿单元表面传播的偏振波。瑞利波传播时,单元表面的能量最强,861总第 567 期 金 属 矿 山 2023 年第 9 期并且沿着球心方向上显著减弱。各向同性材料中,瑞利波呈指数形式衰减;各向异性材料中,则随着深度呈振荡衰减,振荡幅度的包络线呈指数形式。通常情况下,如果传播深度超过两倍波长时,振幅就已非常小22。弹性固体颗粒表层的瑞利波速度计算公式为:vR=1G,(21)式中,G 是颗粒的剪切模量,Pa;是颗粒的密度,kg/m3。G/()是此弹性颗粒横波的速度,m/s;1为瑞利波方程的根。则有:(2-12)4=16(1-21)1-1-22(1-)21,(22)式中,是颗粒的泊松比,则由上式求出 1的近似解为:1=0.163+0.877,(23)瑞利波速度进一步写成:vR=(0.163+0.877)G.(24)由于两接触颗粒的接触作用力仅限于发生在碰撞的两颗粒上,而不应该通过瑞利波传递到其他颗粒上,因此时间步长必须小于瑞利波传递半球面的时间,故可求得时间步长的理论计算公式为:t=RvR=R0.163+0.877G.(25)当由不同粒径的颗粒组成的系统时,时间步长应为:t=R0.163+0.877Gmin.(26)上述公式是在两颗粒为静态或者相对运动速度较小的情况下得出的结论,由于高压辊磨机的颗粒运行并不剧烈,故通过上式取得步长可使仿真模型稳定。但是考虑到高压辊磨机仿真破碎颗粒数量巨大,为使计算稳定可取上式的 1%10%。2.4 仿真参数设定三山岛金矿 KHD3751400-1000 高压辊磨机设备参数和操作参数如表 3 所示。表 4 为该仿真中使用的材料参数和接触参数,这些参数参考了 BARRIOS 等23和 RODRIGUEZ 等21的研究,并根据 EDEM 软件中材料数据库(通用材料模型数据库-GEMM)的安息角结果进行验证。3 联合仿真模型验证与分析为了验证 DEM 和 MBD 联合仿真在高压辊磨机预测处理量、破碎比和单位能耗上的准确性,以及证明 该 联 合 仿 真 方 法 的 优 势,以 三 山 岛 金 矿表 3 高压辊磨机设备参数和操作参数Table 3 Equipment parameters and operating parameters of high-pressure grinding rolls参数类型参数数值设备参数辊径/mm1 400辊宽/mm1 100静态辊隙/mm25操作参数辊速/(r/min)16.8最大喂料粒度/mm35液压系统参数静态氮气压力/MPa5.5工作液压/MPa8.3静态氮气体积/L13油缸数量/个2表 4 离散元仿真参数Table 4 Discrete element simulation parameters类别密度/(kg/m3)剪切模量/GPa泊松比接触参数恢复系数静摩擦系数滚动摩擦系数矿石2 70011080.25 矿石-矿石 0.200.550.51设备7 800700.3矿石-设备 0.150.490.47KHD3751400-1000 型高压辊磨机实际工作数据作为对比目标,对辊隙的动态特性进行验证与分析,以及对前文所述的数学模型、单一离散元仿真模型和联合仿真模型进行对比分析。试验设备、控制系统和整体仿真模型如图 8 所示,在 EDEM 软件中出料口位置设置 GeometryBin(几何选择区域)和 Mass Flow sensor(质量流传感器),用于监测破碎出料颗粒粒径大小和排料流量。高压辊磨机进料的一个重要的特性是,需要一定高度的料柱给位于压辊之间的颗粒床施加压力24。为了避免进料的堵塞,根据 CLEARY 等25的相关结论,通过 DEM 模拟证明,料斗中的物料高度至少应等于其宽度,超过这个高度,物料额外的压力就会逐步传递到料斗壁上,而不影响高压辊磨机的性能。因此,仿真中采用的喂料高度为压辊的宽度,只有当DEM 颗粒达到此高度时,仿真的结果才视为有效结果。并且在仿真过程中,保持相同的颗粒产生速度,以保持破碎过程的稳定。3.1 辊隙的动态变化通过在 Motionview 中添加动辊和定辊之间的位移输出,从而得到辊隙在破碎过程中的动态变化,仿真结果具体如图 9 所示。从图 9 可以看出,物料在 0.5 s 时接触压辊,1.5 s 后喂料高度到达理论要求,在稳定后辊隙保持在35.5 mm 左右。根据矿山提供的数据,破碎过程中KHD3751400-1000 型高压辊磨机的辊隙在 30 40 mm之间,该仿真结果符合工程实际。通过DEM-961 任云鹏等:基于 DEM-MBD 的高压辊磨机仿真模型及应用研究 2023 年第 9 期图 8 试验设备和仿真模型Fig.8 Test equipment and simulation model图 9 高压辊磨机辊隙变化Fig.9 Roll gap change of high-pressure grinding rolls MBD 联合仿真可以模拟高压辊磨机辊隙的动态响应,研究不同工艺参数下辊隙的变化规律,能够帮助设计人员更好地了解系统的性能和特性,优化设计方案。3.2 处理量通过 EDEM 导出的质量流量数据和根据理论计算数据以及高压辊磨机工作数据,绘制各模型预测处理量,结果如图 10 所示。图 10 各模型处理量结果曲线Fig.10 Results curve of treatment volume of each model 比较了 DEM-MBD 联合仿真、DEM 仿真以及处理量数学模型各自的预测值与矿山实际处理量数据。由图 10 可以明显看出,通过数学模型预测的处理量数据和实际工作数据偏差较大,可能是数学模型无法考虑到工作压力、辊隙变化、喂料粒径大小、物料物理性质等其他参数对设备的影响。DEM-MBD 仿真结果与 DEM 仿真结果和实际工作数据比较接近,但DEM 仿真结果比实际数据偏低,原因可能是设备的处理量和辊隙有直接关系,在破碎过程中,辊隙因受到物料的作用力而发生改变,DEM 模型无法仿真辊隙的动态变化,造成误差。为减小误差和避免偶然性,每个模型仿真 2 次,各模型预测处理量数据误差如表 5 所示。表 5 各模型处理量结果对比Table 5 Comparison of treatment volume results by each model方案试验序号平均值/(t/h)误差/%DEM-MBD1597.27-8.112587.76-9.58DEM1498.68-23.282511.91-21.24数学模型1355.68-45.282355.68-45.28实际数据1650026500 分析表 5 可知,DEM-MBD 联合仿真模型预测处理量准确度最高,其次是 DEM 模型,数学模型准确度最低。DEM-MBD 联合仿真模型预测值与实际测定值相比,虽然误差更小,但仍然偏低,可能是因为在实际生产中高压辊磨机辊面是柱钉面,而仿真中采取的是光滑面,降低了辊面对矿石的一部作用。在前期仿真试验过程中发现,数量大且尺寸较小的柱钉会极大增加仿真中细小单元格数量,影响仿真进度,考虑到计算机性能,在不对仿真结果造成较大误差的情况下,将柱钉辊面简化成光滑辊面。3.3 单位能耗在 EDEM 软件中导出压辊的转矩,再利用公式(8)和(9)求得设备的单位能耗,将此结果与单位能耗数学模型求得的单位能耗,以及工业生产中的实际单位能耗数据进行对比,结果如图 11 所示。由图 11 可以看出,DEM-MBD 联合仿真模型预测的单位能耗和实际数据拟合度非常高,DEM 仿真模型和数学模型预测值比较接近,但都比实际数据大。可能是因为在其他条件一致的情况下,单位能耗和辊隙大小成反比,辊隙越大,单位能耗越低。DEM模型和数学模型无法模拟辊隙的变化,所以造成预测值偏大。仿真 2 次,各模型预测单位能耗数据误差见表 6。071总第 567 期 金 属 矿 山 2023 年第 9 期图 11 各模型单位能耗曲线Fig.11 Unit energy consumption curve for each model表 6 各模型单位能耗结果对比Table 6 Comparison of energy consumption results per unit for each model方案试验序号平均值/(kW h/t)误差/%DEM-MBD10.60+5.2620.59+3.51DEM10.82+43.8620.81+42.11数学模型10.93+63.1620.93+63.16实际数据10.57020.570 分析表 6 可知,DEM-MBD 联合仿真模型在预测单位能耗时,准确度最高,最小误差仅有 3.51%,具有比较满意的结果,为后续预测单位能耗提供了一种更为准确的方式。3.4 破碎比表 7 为各模型预测平均破碎比与实际设备测定平均破碎比的相对误差。比较表中数据,发现 DEM-MBD 联合仿真模型预测的平均破碎比与实际数据拟合度最好,其次是 DEM 模型,但 DEM 模型预测数据比实际数据要高,这是因为在高压辊磨机中,压辊通过高压将物料挤压成片,使其在辊隙中受到高压力的挤压和剪切作用。辊隙越小,物料在挤压和剪切过程中受到的力度越大,破碎效果越好,破碎产品的细度也就越高。DEM 模型没有考虑到辊隙的动态变化,一直以初始辊隙运行,所以造成破碎比偏大。表 8 分析总结了高压辊磨机各性能指标的不同模型仿真结果与实际数据的平均误差值,在初步建模阶段,应该认为仿真结果误差在一定范围内并显示正确的趋势,则 DEM-MBD 联合仿真模型是成功的。4 结 论(1)以 KHD3751400-1000 型高压辊磨机为研究对象,为提高高压辊磨机处理量、单位能耗、破碎比预测的准确度,本文首次将 DEM-MBD 联合仿真技术应用在高压辊磨机虚拟仿真上,并详细介绍联合仿真表 7 各模型平均破碎比对比Table 7 Comparison of the average crushing ratio of each model方案序号破碎比误差/%DEM-MBD12.38-10.8622.33-12.73DEM13.18+19.1023.15+17.98数学模型14.80+79.7824.80+79.78实际数据12.67022.670表 8 各模型的性能指标相对误差Table 8 Relative error of performance indicators for each model%方案处理量单位能耗破碎比DEM-MBD8.854.3911.80DEM22.2642.4918.54数学模型45.2863.1679.78方法和相关参数设定,研究表明 DEM-MBD 仿真能够更真实地模拟设备的运行状态。(2)为评价该联合仿真模型对高压辊磨机性能指标的预测准确度,将 DEM-MBD 模型分别与 DEM模型以及前文所述的处理量、单位能耗、破碎比数学模型和实际工作数据进行比较。通过对比处理量、单位能耗、平均破碎比三大性能指标,综合考虑各模型的精确度和稳定性,结果表明 DEM-MBD 模型预测准确度最为理想。DEM-MBD 联合仿真模型比现有的单一 DEM 仿真模型在预测处理量、单位能耗、破碎比上的误差值分别减小了 13.41、38.10、6.74 个百分点。(3)高压辊磨机 DEM-MBD 联合仿真相比于DEM 仿真能够更真实、更准确地模拟设备的运行状态和性能,从而使高压辊磨机仿真模型的准确度得到显著改善,为此类仿真提供合理参考,降低了高压辊磨机的研发成本,也为高压辊磨机的设计和优化提供了更科学的依据。然而,对高压辊磨机产品粒径分布的仿真还有未尽之处,后续将针对此开展进一步研究。参 考 文 献1 魏波,张宪伟,李丽匣,等.高压辊磨机粉碎工艺国外应用进展与发展趋势J.金属矿山,2022(2):10-18.WAN Bo,ZHANG Xianwei,LI Lixia,et al.Progress and develop-ment trend of high pressure grinding rollers comminution circuits a-broadJ.Metal Mine,2022(2):10-18.2 邵泽富,曹品安,李亚娥.试析高压辊磨机在矿物加工工程中的171 任云鹏等:基于 DEM-MBD 的高压辊磨机仿真模型及应用研究 2023 年第 9 期应用J.世界有色金属,2021(9):215-216.SHAO Zefu,CAO Pinan,LI Yae.Application of high pressure roller mill in mineral processing engineeringJ.World Nonferrous Met-als,2021(9):215-216.3 袁满,卢冀伟,张俊飞,等.高压辊磨机在有色金属矿山的应用进展与展望J.中国矿业,2023,32(2):90-96.YUAN Man,LU Jiwei,ZHANG Junfei,et al.Research progress and prospect of high pressure grinding rolls used in non-ferrous meta