基于流固耦合的搅拌磨磨矿离散元仿真及试验研究.pdf
《基于流固耦合的搅拌磨磨矿离散元仿真及试验研究.pdf》由会员分享,可在线阅读,更多相关《基于流固耦合的搅拌磨磨矿离散元仿真及试验研究.pdf(7页珍藏版)》请在咨信网上搜索。
1、Series No.568October 2023 金 属 矿 山METAL MINE 总 第568 期2023 年第 10 期收稿日期 2023-02-28基金项目 郑洛新国家自创区创新引领型产业集群专项(编号:201200210500)。作者简介 刘 伟(1999),男,硕士研究生。通信作者 薛玉君(1971),男,教授,博士,博士研究生导师。基于流固耦合的搅拌磨磨矿离散元仿真及试验研究刘 伟1,2 刘 俊3,4 程 波3,4 薛玉君1,2,4 李济顺1,2,4(1.河南科技大学机电工程学院,河南 洛阳 471003;2.河南省机械设计及传动系统重点实验室,河南 洛阳 471039;3.洛
2、阳矿山机械工程设计研究院有限责任公司,河南 洛阳 471039;4.智能矿山重型装备全国重点实验室,河南 洛阳 471039)摘 要 立式螺旋搅拌磨是矿石超细磨作业的关键设备,为研究其在实际工况下磨矿效果的影响因素,分析其磨矿特征,以立式螺旋搅拌试验磨机为研究对象,采用 CFD-DEM 方法建立矿浆与研磨介质的流固耦合仿真模型,对不同尺寸和形状的研磨介质在磨矿过程中的运动速度、碰撞力、碰撞次数以及搅拌器的扭矩进行分析,并使用试验磨机进行磨矿试验,提取搅拌器扭矩及不同粒度的筛下累积量。结果表明,立式螺旋搅拌磨的磨矿效果与研磨介质的尺寸和形状有着密切的关系,使用钢球介质时比使用钢段介质时的磨矿效果
3、好,小钢球比大钢球的磨矿效果好,且使用小钢球产生的能耗较低。仿真所提取的扭矩与实验所测得的扭矩基本吻合,误差在 5%以内,该仿真模型具有一定的合理性。关键词 立式螺旋搅拌磨 流固耦合 离散单元法 研磨介质 磨矿特征 中图分类号TD453 文献标志码A 文章编号1001-1250(2023)-10-189-07DOI 10.19614/ki.jsks.202310027Discrete Element Simulation and Experimental Study of Grinding of Stirring Mill Based on Fluid-Structure Interactio
4、nLIU Wei1,2 LIU Jun3,4 CHENG Bo3,4 XUE Yujun1,2,4 LI Jishun1,2,4(1.School of Mechanical and Electrical Engineering,Henan University of Science and Technology,Luoyang 471003,China;2.Key Laboratory of Mechanical Design and Transmission System in Henan Province,Luoyang 471039,China;3.Luoyang Mining Mac
5、hinery Engineering Design Research Institute Co.,Ltd.,Luoyang 471039,China;4.National Key Laboratory of Intelligent Mining Heavy Equipment,Luoyang 471039,China)Abstract Vertical spiral stirring mill is the key equipment of ore ultrafine grinding.In order to study the influencing fac-tors of the grin
6、ding effect under actual working conditions and analyze its grinding characteristics,taking the vertical spiral stir-ring test mill as the research object,the CFD-DEM method was used to establish the fluid-structure interaction simulation mod-el of the pulp and the grinding media,and the motion spee
7、d,collision force,collision number and the agitators torque of the grinding media of different sizes and shapes in the grinding process were analyzed,and the grinding test was carried out by the test mill to extract the torque of the agitator and the accumulation under sieve of different particle si
8、zes.The results show that the grinding effect of vertical spiral stirring mill is closely related to the size and shape of the grinding media,the grinding effect of the steel ball media is better than that of the steel section media,the grinding effect of the small steel ball is better than that of
9、the large steel ball,and the energy consumption generated by the use of small steel ball is lower.The torque extracted by the simulation is basically consistent with the torque measured in the experiment,and the error is within 5%,and the simulation model has certain rationality.Keywordsvertical spi
10、ral stirring mill,fluid-structure interaction,discrete element method,grinding media,grinding characteristics 立式螺旋搅拌磨是一种超细磨设备,常应用于精矿再磨和尾矿再磨工艺,其能量利用率高,产品粒度细,最小可将矿石研磨至 10 m 以下,且结构简单,占地面积小,对建设绿色矿山具有重要意义1,因此,对立式螺旋搅拌磨进行磨矿过程仿真分析也愈发重要。而研磨介质作为磨矿设备的重要组成部分,如981何在降低成本的同时,又能达到较好的磨矿效果一直以来都是选矿工作者的难题,选择合适的研磨介质对磨矿本
11、身具有重大意义2。近年来,离散单元法逐渐兴起,该方法通过建立固体颗粒体系的参数化模型,进行颗粒的行为模拟和分析,被广泛应用于散状物料和粉体工程领域3。但是,立式螺旋搅拌磨在实际工况下采用湿法研磨的方法,筒体内部呈多相流混合状态,在螺旋搅拌器的作用下,矿物和研磨介质做自转和螺旋上升运动,通过研磨介质的冲击、挤压和摩擦对物料进行研磨来实现矿石的粉碎。单独采用离散单元法并不能够对立式螺旋搅拌磨的磨矿过程进行准确模拟,而采用流固耦合的方法进行仿真分析能够更准确地分析其在实际工况下的磨矿特征。国内外学者利用流固耦合方法对立式螺旋搅拌磨的磨矿机理及研磨行为展开了大量研究。STRO-BEL 等4以小型搅拌磨
12、机为研究对象,分析了介质球尺寸和流体黏度对物料应力状态的影响,结果表明,较大的磨粒可得到较高的应力能,而黏度的增加会降低应力能。PRZIWARA 等5研究了在搅拌磨机中助磨剂对磨粒稳定性的影响,结果表明,磨机内部矿浆的流动行为会对磨矿效果产生较大的影响。卢世杰等6建立了立式螺旋搅拌磨的 CFD 单相流仿真模型,分析了不同转速和不同螺旋直径对立磨机内流场速度的影响,结果表明,转速或螺旋直径增大时,流场速度呈增大趋势。孙新明等7-8采用流固耦合的方法建立了立式螺旋搅拌磨仿真模型,得出了螺旋搅拌器运行的阻力矩,同时,还进行了螺旋搅拌器结构强度分析,为磨机的设计和制造提供了依据。孙小旭等9-10建立了
13、超细磨用搅拌装置流体力学仿真模型,得到了搅拌装置的流场分布状态和输入功率变化情况,为实验和工业中超细磨技术和设备的研发奠定了基础。研磨介质的尺寸和形状会对立式螺旋搅拌磨的磨矿效果产生很大的影响,采用 CFD-DEM(Computa-tional Fluid Dynamics-Discrete Element Method)的方法建立矿浆与研磨介质的流固耦合模型,探究立式螺旋搅拌磨在采用不同尺寸和形状的研磨介质时的磨矿特征,并结合立式螺旋搅拌试验磨机的磨矿试验结果,分析研磨介质的尺寸和形状与磨矿效果的关系,为立式螺旋搅拌磨的研磨介质选型提供依据。1 流固耦合模型的建立1.1 流体连续相控制方程以
14、立式螺旋搅拌磨内部矿浆为对象建立流体模型,由于矿浆处于恒温恒压下,且考虑到其与研磨介质的相互影响,在守恒方程中加入一个额外的体积分数项11,由流体力学守恒理论可知,其守恒方程有:质量守恒方程:llt+?llul=0,(1)式中,l为空隙率,l为矿浆密度,kg/m3;t 为时间,s;ul为矿浆流速,m/s。动量守恒方程:llult+?llull=?(l?ul)+llgl-?l-S,(2)式中,gl为矿浆重力加速度,m/s2;l为矿浆动力黏度,Pa s;S 为动量汇,由以下公式算得:S=ni=1FDV,(3)式中,FD为研磨介质所受矿浆作用力的总和,N;V 为控制体的体积,m3;n 表示控制体内研
15、磨介质的数量。1.2 颗粒离散相控制方程以研磨介质为对象建立离散元模型,根据牛顿第二定律可推导其平衡方程12为mpdupdt=Fg+Fd+Fb+Fc,(4)式中,mp为研磨介质质量,kg;up为研磨介质速度,m/s;Fg为研磨介质自身的重力,N;Fd为研磨介质所受的曳力,N;Fb为浮力,N;Fc为研磨介质之间以及研磨介质与壁面之间的接触力,N。曳力模型选择 Wen Yu&Ergun 模型13,其计算公式可表示为Fd=Vpufuf1-l,(5)式中,Vp为研磨介质体积,m3;uf为研磨介质与矿浆的相对流速,m/s。由下式计算获得:=1501-l()2llL2+1.751-l()lufL,l 0.
16、834CDll-1.651-l()uf,l 0.8 ,(6)式中,L 为研磨介质粒径,m;CD为曳力系数,由以下公式算得:091总第 568 期 金 属 矿 山 2023 年第 10 期CD=24Re,Re 0.524Re1+0.15Re0.687(),0.5 1 000,(7)式中,Re 为雷诺数,由以下公式算得:Re=llufLl.(8)除了曳力,还需考虑矿浆对研磨介质产生的浮力,浮力计算公式如下:Fb=lgVp.(9)研磨介质之间、研磨介质与壁面之间的接触模型选择 Hertz-Mindlin(no-slip)模型11,其计算公式如下:法向接触力 Fn及法向阻尼力 Fdn:Fn=43ER3
17、2n,(10)Fdn=-256lneln2e+2Snmureln,(11)式中,E为等效弹性模量,Pa;R为等效半径,m;n为法向重叠量,m;m为等效质量,kg;ureln为相对速度的法向分量,m/s;e 为恢复系数,Sn为法向刚度,N/m;由以下公式算得:Sn=2ERn,(12)切向接触力 Ft及切向阻尼力 Fdt:Ft=-Stt,(13)Fdt=-256lneln2e+2Stmurelt,(14)式中,t为切向重叠量,m;urelt为相对速度的切向分量,m/s;St为切向刚度,N/m;由以下公式算得:St=8GRn,(15)式中,G为等效剪切模量,Pa。滚动摩擦力 Fr由以下公式算得:Fr
18、=-FnRii,(16)式中,为摩擦系数,Ri为研磨介质质心到接触点的距离,m;i为研磨介质接触点的角速度,rad/s。2 磨机模型的建立2.1 三维模型及网格划分立式螺旋搅拌试验磨机的结构非常复杂,在建模过程中需要对其进行一定的简化。试验所采用的物料为赤铁矿,密度为 3 900 kg/m3,入料粒度为-1 mm,研磨介质的直径一般为 810 mm,两者直径相差很大,仿真无法准确模拟出赤铁矿的破碎过程。因此,将赤铁矿与水混合的矿浆简化成一种单一的流体14,默认磨机内部只存在矿浆和研磨介质两种物质,将筒体及搅拌器视为无厚度的壁面。由于在FLUENT 中要通过滑移网格来实现搅拌器的旋转运动,因此,
19、需要对立式螺旋搅拌磨的流体域划分旋转域和静止域,旋转域为主动旋转运动区域,静止域由旋转域带动其进行旋转运动。图 1 和表 1 分别为立式螺旋搅拌试验磨机部分结构参数及三维模型。图 1 立式螺旋搅拌试验磨机三维模型Fig.1 3D Model of vertical spiral stirring test mill表 1 立式螺旋搅拌试验磨机结构参数Table 1 Structural parameters of vertical spiral stirring test millmm参数搅拌轴直径 D1搅拌器直径 D2旋转域直径 D3搅拌筒内径 D4螺旋桨高度 H1旋转域高度 H2搅拌筒高度
20、 H3数值40280296376800816970在 DESIGN MODEL 中进行流体域抽取,利用FLUENT MESHING 对模型进行网格划分,静止域网格最大尺寸控制为 20 mm,旋转域网格最大尺寸控制为 15 mm,对靠近搅拌器的网格再进行适当加密,网格类型选择多面体网格。经过质量提升、节点移动,最终得到如图 2 所示的网格,网格数量为 4.75 万个,最小正交质量为 0.36。将旋转域和静止域的交界面定义为 interface,用于数据交换。2.2 耦合模型设置在 FLUENT 中输出边界网格,将边界网格导入到EDEM 中。在 EDEM 中对研磨介质进行预填充,预填充时搅拌器不进
21、行旋转运动,以保证与 FLUENT 耦合时网 格 初 始 状 态 相 同,EDEM 时 间 步 长 选 为0.000 001 s,点击 Start Coupling Server,使 EDEM 处于耦合接听状态。在 FLUENT 中选择瞬态求解方式,导入耦合接口191 刘 伟等:基于流固耦合的搅拌磨磨矿离散元仿真及试验研究 2023 年第 10 期图 2 立式螺旋搅拌试验磨机网格剖视图Fig.2 Section view of mesh of vertical spiral stirring test millUDF。在 FLUENT 与 EDEM 的耦合中,能够实现固体与流体双向耦合的模型有
22、两种:Lagrangian 模型和Eulerian 模型。Lagrangian 模型只考虑固体和流体之间的动量交换,FLUENT 采用单相流计算,一般情况下,该模型只适用于稀相流,即颗粒相体积分数小于10%,且颗粒之间没有相互作用。Eulerian 模型除了考虑固体和流体之间的动量交换外,在 FLUENT 中还需开启欧拉多相流模型,利用颗粒位置、体积等计算颗粒在网格内体积分数。在立式螺旋搅拌磨中,研磨介质在磨机内部呈堆积形态,矿浆与介质、介质与介质之间都会发生相互作用,因此,采用 Eulerian 模型进行耦合。曳力模型选择 Wen Yu&Ergun 模型,连接成功后 FLUENT 会自动打开
23、欧拉多相流,设置两种欧拉相,主相为流体相,次相为离散相。湍流模型选择 k-epsilon-RNG 模型,标准壁面函数。FLUENT 转速与 EDEM 相同,出口设置为压力出口,时间步长为 EDEM 时间步长的 100 倍,即 0.000 1 s。由于 FLUENT 在仿真开始时引入颗粒会对内流场湍流引起较大变化,极易引起仿真发散,因此,须将体积分数项松弛因子调小,在 FLUENT 内流场稳定后再逐渐调大。在仿真中,通过改变密度和黏度来控制矿浆的浓度,矿浆密度 l由以下公式算得:C=l-wm-wml%,(17)式中,C 为矿浆质量浓度,%;w为水的密度,kg/m3;m为干矿密度,kg/m3。矿浆
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 耦合 搅拌 磨磨 离散 仿真 试验 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。