基于CFD超细磨用搅拌装置关键参数影响研究.pdf

7 6 有色金属 选矿部分2 0 2 0 年第1 期 d o i 1 0 ．3 9 6 9 ／j ．i s s n ．1 6 7 1 9 4 9 2 ．2 0 2 0 ．0 1 ．0 1 4 基于C F D 超细磨用搅拌装置关键参数影响研究 孙小旭1 ’2 ，祝启恒1 ’2 ，姚建超1 ’2 ，何建成1 ’2 1 ．北矿机电科技有限责任公司，北京1 0 0 1 6 0 ； 2 ．北京矿冶科技集团有限公司，北京1 0 0 1 6 0 摘要在分析超细磨技术要求的前提下，通过分析盘式、螺旋式、叶轮式、棒式四种常用搅拌装置型式，确定了盘式是 较为适合的超细磨用搅拌装置。利用计算机仿真流体力学建立了仿真分析模型，进行仿真分析。通过仿真得到了搅拌装置 的流场分布状态和输入功率变化情况，重点分析了转速、搅拌元件直径、间距、面积比与流场、功率间的影响关系及变化范围。 通过仿真明确了搅拌转速的范围为3 0 0 ～8 0 0r ／m i n ，搅拌元件直径的范围在1 5 0 ～1 7 0I T l r f l ，搅拌元件间距的范围为3 0m m 左 右，空隙面积比的范围为0 ．1 5 ～o ．3 5 。各参数的确定为实验和工业中超细磨技术和设备的研发奠定了基础。 关键词超细磨；盘式搅拌装置；C F D ；流场；输入功率 中图分类号T D4 5 3 ．9文献标志码A文章编号1 6 7 1 9 4 9 2 2 0 2 0 0 1 0 0 7 6 0 6 A n a l y s i so fK e yP a r a m e t e r sI n f l u e n c eo nU l t r a f i n eG r i n d i n gA g i t a t o rB a s e do nC F D S U NX i a o x “1 ～，Z H U Q i h e n 9 1 ～，Y A OJ i a n c h a 0 1 一，H EJ i a n c h e n 9 1 ’2 j ．B G R I M MM a c h i n e r ya n dA u t o m a t i o nT e c h n o l o g yC o ．，L t d ．，B e i j i n g1 0 0 1 6 0 ，C h i n a ； 2 ．B G R I M MT e c h n o l o g yG r o u p ，B e i j i n g1 0 0 1 6 0 ，C h i n a A b s t r a c t B a s e do nt h ea n a l y s i so ft h et e c h n i c a lr e q u i r e m e n t so fu l t r a f i n eg r i n d i n g ，f o u rt y p e so f c o m m o ns t i r r e da g i t a t o r si n c l u d i n gd i s c ，i m p e l l e r ，s c r e wa n db a ra r ea n a l y z e d ，a n di ti sc o n f i r m e dt h a tt h e d i s cs t i r r e da g i t a t o ri ss u i t a b l ef o ru l t r a f i n eg r i n d i n g ．T h es i m u l a t i o nm o d e lo ft h ea g i t a t o ri se s t a b l i s h e db y c o m p u t a t i o n a lf l u i dd y n a m i c s ，a n dt h em o d e li ss i m u l a t e d ．T h r o u g ht h es i m u l a t i o n ，t h ed i s t r i b u t i o n0 ft h e f l o wf i e l da n dt h ec h a n g eo ft h ep o w e rc o n s u m p t i o na r ea n a l y z e d ，a n dt h er e l a t i o n s h i pb e t w e e ns p e e d ， a g i t a t o rd i a m e t e r ，d i s t a n c e ，a r e ar a t i oa n dt h ef l o wf l i e d ，p o w e rc o n s u m p t i o na r es t u d i e de m p h a t i c a l l y ． A c c o r d i n gt ot h er e s u l t s ，t h er a n g eo ft h es p e e di s3 0 0 ～8 0 0r ／m i n ，t h er a n g eo ft h ed i a m e t e ri s 15 0 ～ 1 7 0m m ，t h er a n g eo ft h ed i s t a n c ei sa b o u t3 0m m ，a n dt h er a n g eo ft h eg a pa r e ar a t i oi s0 ．1 5 ～0 ．3 5 ．T h e c o n f i r m a t i o no ft h ep a r a m e t e r sm a k e saf o u n d a t i o nf o rt h ed e v e l o p m e n to fu l t r a f i n eg r i n d i n gt e c h n o l o g yi n e x p e r i m e n t a la n di n d u s t r i a lf i e l d ． K e yw o r d s u l t r a f i n eg r i n d i n g ；d i s cs t i r r e da g i t a t o r ；C F D ；f l o wf i e l d ；p o w e rc o n s u m p t i o n 伴随着国内外矿物的深度开采和加工，贫、细、 杂矿物逐渐进入矿物采选冶领域，有用矿物的嵌布 粒度越来越细，且普遍存在于国内外有色、黑色、非 金属矿物及化工，石油等领域。为有效提高有用矿 物品位，提升产品质量，提高产品适用性等，高效超 细磨技术及设备逐渐进入粉磨领域，且其作用也越 来越重要。以国内某种需要超细磨的浆料为例，为 满足工艺需求和产品性能，该浆料必须细磨至d 。。≤ 1 0 肚m ，且其固含量低 ≤3 5 ％ ，黏度高，而现存装备 具有效率低、能耗高、实用性差等缺点[ 1 ] 。因此，亟 需研究用于高效超细磨的技术及装备。 众所周知，用于超细磨技术及装备的关键参数 主要包括结构和工艺参数，本文主要研究了搅拌装 置关键结构参数与超细磨装备流场及功率间的相互 影响关系[ 2 ] 。 基金项目国家自然科学基金资助项目 5 1 6 0 4 0 2 7 ；5 1 6 7 4 0 3 3 收稿日期2 0 1 9 - 0 5 2 4 作者简介孙小旭 1 9 8 9 一 ，男，河北易县人，硕士。工程师，主要从事粉磨技术的研究及推广工作。 万方数据 2 0 2 0 年第1 期孙小旭等基于C F D 超细磨用搅拌装置关键参数影响研究 7 7 1 搅拌装置选择 搅拌装置是细磨、超细磨技术及装备的关键部 件，搅拌装置一般由搅拌轴和数量不等的搅拌元件 组成，搅拌元件带动介质从而产生细磨作用，其结构 形式对介质运动特性具有重要影响。常见搅拌装置 如图1 所示，主要有盘式、螺旋式、叶轮式和棒式四 类，四类搅拌装置对应的搅拌元件分别为搅拌盘、螺 旋、叶轮和搅拌棒。通过对常用搅拌装置分析可知， 盘式搅拌装置由搅拌轴和多层盘形搅拌元件组成。 搅拌盘工作面积较大，容易使介质产生较强的轴向 和径向运动。其搅拌作用很强，承受磨损的能力也 较强，而且在较高的输人转速下，内部可以获得较高 的能量输入，介质和矿浆的运动更加紊乱，比较符合 高效超细磨技术的要求口] 。同时，搅拌盘形状简单， 加工制造容易，T 作阻力小，适用性好，既可以立式 1 二作，也可以卧式工作，是目前较适合高效超细磨技 术的装置。 謇罪如豢 a 盘式 b 螺旋式 c 叶片式 d 棒式 图1常用搅拌装置类型 2 模型构建及分析 2 ．1 计算流体力学 计算流体力学 C F D C o m p u t a t i o n a lF l u i d D y n a m i c s 是2 0 世纪7 0 年代初逐渐形成的仿真计 算方法，将流体分析采用数值模拟的方法进行求解， 建立各种条件下的守恒方程‘⋯，通过计算机求解各 种条件下的偏微分方程，得到流体控制方程的近似 解，进而实现对流场及相关参数的仿真分析⋯。 2 ．2模型建立 以盘式搅拌磨机为模型进行建模，适当简化，以 水作为媒介，流体运动形式为湍流∽】，仿真模型结构 如图2 所示。 图2 盘式搅拌磨机 F i g ．2 S t r u c t u r eo fd i s cs t i r r e dm i l l 2 ．3 网格划分 四面体网格适合较为复杂的几何图形，简单方便； 六面体适合表征两相交界面，形貌清楚，质量较高。本 模型采用四面体网格进行划分，共划分出4 5 53 9 4 个网 格，如图3 所示。 a 搅拌装置网格 b 简体网格 图3 模型网格 F i g ．3M o d e lm e s h 万方数据 7 8 有色金属 选矿部分2 0 2 0 年第1 期 仿真分析及结果讨论 3 ．1 转速影响分析 选取多个速度值进行仿真分析，其中1 5 0 、3 0 0 、 5 0 0 、8 0 0r ／r a i n 下的仿真结果如图4 所示。 由图4 中知，不同转速下的速度分布情况较为 围及到筒壁之间区域，在其它区域，细磨效果较差。 在流线图中，可以看出随着转速增加，流场的流线形 状基本一致，在搅拌装置附近形成明显紊流运动，验 证了搅拌装置附近是细磨作用主要发生区域的 观点。 从图4 可知速度越大，细磨效果越佳，但是通过 类似，速度大小的分布很不均匀在搅拌装置附近以图5 分析可知，随着速度增加，输入功率呈几何倍数 及其外缘到筒壁间，速度梯度较大；其余部分，包括增加，因此，在确定搅拌装置速度时要综合考虑研磨 搅拌元件之间及搅拌装置上部，其速度梯度较小。强度、输入功率等多方面因素，实验室常规应用的速 因此可认为磨机的细磨作用主要发生在搅拌装置周 良 岜 { j { L 督 群篙飘m3 们即”。 } 12 9 5 e 0 0 0 ㈠．2 1 3 3 1 4 e e 0 0 0 0 0 0 篷9 1 ．7 0 5 协2 e 0 0 1 0 0 0 E 0 9 2 9 0 1 e e _ 0 - 0 0 0 ■7 ．2 8 3 e - 0 0 1 雕．4 6 洲7 4 e - 0 0 0 艮器；器{ l ；．4 2 2 3 ㈣7 e 。 器0 ■1 ．6 1 8 e - 0 0 1 ■．0 0 0 9 0 e 2 e - 0 0 0 0 2 0 【ms “- 1 1 v e c V e l o l o r c l P3 忏⋯。 9 ‘o j 酬o o 几13 3 5 e 0 0 0 U 瞄10 0 1 e 0 0 0 _ ■6 ．6 7 6 e ．0 0 1 _ ■3 ．3 3 8 0 ．0 0 1 _ 1 00 0 0 e 0 0 0 【ms “_ 1 】 a 1 门 1 5 0r ／r a i n 恻V e l o 。c r l 妒8 斛。⋯。‘9 ‘叫。。㈣ V c 。e №l o c ⋯i t y i n s l n F r a m e l } 4 4 ．。3 4 1 6 6 。e 。0 。0 鹾；, 5 7 0 7 7 6 e e 0 0 0 0 0 0 簟；．眦2 3 7 e 0 0 0 0 0 0 E 4 6 2 9 洲7 e 0 0 0 0 0 ■2 1 5 8 e 0 0 0 l { ．6 8 8 协8 e 0 0 0 0 0 0 l ．0 3 4 7 9 e 9 e 0 0 0 0 I i5 8 ．0 3 9 9 5 2 e e - 0 0 0 0 l ；．0 6 0 0 9 7 e e 0 0 - 0 0 1 0 f ms 1J e n 5 0 0r ／m i n 度值一般在3 0 0 ～8 0 0r ／r a i n 。 v e c V e l o l o r c l ”3 ”” 9 ‘0 j 删。“’黜蒿飘们叶“⋯8 厂T26 7 0 e 0 0 0 ■2 ．0 0 3 e 0 0 0 _ ■13 3 5 e 0 0 0 _ ■6 ．6 7 6 e O O l _ 1 00 0 0 e 0 0 0 f ms “1 j I 豳．5 2 8 8 e 9 e 十 0 0 0 0 0 0 I 基2 2 ．2 1 0 6 4 6 e e 0 0 0 0 0 0 I I ．7 9 4 8 0 2 e e 0 0 0 0 0 0 l { ．4 6 5 1 8 7 e e 0 0 0 0 一．2 1 3 3 9 5 e o e 0 0 0 0 0 蟹．纰7 1 0 e - 0 0 - 0 0 l 怒岩{ I I ；．。2 3 协7 e - 。0 。0 { 00 0 0 e 0 0 0 【r r ls ‘一1 】 b n 3 0 0r ／r a i n v V e c e l o I o c ，々m s t n F r a m 。 P ‘0 j 。硝帕“’2 1 盛％i “s t n F r a m e r r71 2 1 e 0 0 0 U 叫53 4 1 e 0 0 0 - ■3 ．5 6 0 e 0 0 0 _ ■17 8 0 e 0 0 0 _ 1 00 0 0 e 0 0 0 【ms “一1 】 图4横截面速度云图和流线图 l 嘞．9 0 7 4 5 e e 0 0 0 0 b } 60 4 2 0 0 0 0 譬i ．6 1 7 9 1 0 e e 0 0 O G o O k ．7 3 4 1 6 e 7 e 0 0 0 0 ■3 ．8 8 4 e 0 0 0 I ；．4 0 2 1 5 3 e e 伽 0 0 0 I ；．5 1 5 8 8 9 e 0 e 0 0 0 0 0 l ．7 。2 6 。。e 0 。。0 l 4 8 ．3 6 3 1 6 1 枷e - 0 1 0 0 0 0 e 0 0 0 【ms “- 1 】 d n 8 0 0r ／r a i n F i g ．4 C r o s ss e c t i o n a lv e l o c i t yn e p h o g r a m sa n df l o wc h a r t s 转速／ r m i n 。 图5 不同速度下的输入功率 F i g ．5 P o w e rc o n s u m p t i o no fd i f f e r e n ts p e e d 3 ．2 搅拌元件间距影响分析 为保证盘与盘之问有充足的空间使介质通过， 且达到较好的相互作用效果，在相同条件下，选取多 个不同的盘间距进行仿真分析，其中盘问距为3 0 、 4 0 、5 0m m 的仿真结果如图6 所示。 从图6 可知，随着盘问距增大，流场流速与梯度 基本保持不变，即盘间距对流速影响不大；但在流线 图中，随着盘问距变大，流体紊流运动增强，即细磨 作用更强烈，细磨效果更好。 如图7 所示，功率随着盘问距的变化而变化，但 稳定在一定的区问内，在间距为3 0m m 时达到最大 值。故在满足制造工艺与实际应用的前提下，适当 增大间距可以提升细磨效果。 3 万方数据 a 间距3 0m m c 间距5 0m m 图6 不同间距速度云图和流线图 f b l h 3 距4 0m m F i g ．6V e l o c i t yn e p h o g r a m sa n df l o wc h a r t so fd i f f e r e n td i s t a n c e 直径／r a m 图7 不同间距的输入功率 F i g ．7 P o w e rc o n s u m p t i o no fd i f f e r e n td i s t a n c e s 3 ．3 搅拌元件直径影响分析 为确定搅拌盘直径的影响，改变搅拌盘直径并 逐一进行仿真，得出筒体内速度流线、云图及输入功 率，如图8 和图9 所示。 由图8 分析可知，当直径为1 5 0 、1 7 0m m 时，在 搅拌盘之间及搅拌盘外缘区域，流体运动强烈且存 在较大速度差，有利于介质的碰撞摩擦，因此可以认 为此区域的细磨效果最佳。 在层流区输入功率与搅拌元件直径的d 3 成正 比，在湍流区搅拌功率与搅拌元件直径的d 5 成正 比，但由于直径变化会导致搅拌装置、筒体之间的几 何关系 d ／D 发生变化，当d ／D 发生变化时，搅拌功 率准数N ，也会发生变化，因此必须给予修正，对于 常用搅拌装置其修正系数如下 1 则搅拌装置的计算功率输入如下 P 7 一刷，p N 3 d 5 2 万方数据 ．8 0 ．有色金属 选矿部分 2 0 2 0 年第1 期 富 ≤ 褂 愿 a 直径1 5 0m m v e l 9 e 耐i nS t nF 阳m e P r o j e c t i o n V e c t o r l 53 4 1 e 0 0 0几． 4855e0005 7 9 e 0 0 0 l } ． 1 } ．1 f h l 直径1 7 0m m c 直径2 0 0m m 图8不同直径的速度流线及云图 d 直径2 2 0m m F i g ．8 V e l o c i t yc o n t o u r sa n dC R O S Ss e c t i o nf l o wc h a r to fd i f f e r e n td i a m e t e r s 直径／m m 图9 直径与输入功率的关系 F i g ．9R e l a t i o n s h i pb e t w e e nd i a m e t e ra n dp o w e r 图9 为直径与输入功率的关系曲线图，随着直 径d 的增加，输入功率与d 5 基本成正比，表明仿真 分析对功率预测具有指导意义。 3 ．4 搅拌元件空隙面积比影响分析 空隙面积比 空隙与搅拌元件的面积比 是影响 细磨效果的重要参数之一，在不同比值下进行仿真 分析，并讨论分析结果。 图1 0 为不同比值下的速度流线图及云图，分析 可知，随着比值增加，径向流动变剧烈，相邻两搅拌 元件间的湍流运动也随之增强，因此增大面积比可 以有效强化作用区域并增加径向流动，强化细磨 效果。 万方数据 2 0 2 0 年第1 期孙小旭等基于C F D 超细磨用搅拌装置关键参数影响研究 V e l o c i t yi nS t nF r a m e P r o j e c t i o n V e l o c i t yi nS t nF 阳m e V e c t o r1C o n t o u r1 4855e0005 7 9 e 0 0 0 l } ． 1 } ．1 a 比值0 ．0 6 v e l o c i t yi nS t nF r a m e P r o j e c t i o n V e c t o r l 53 4 1 e 0 0 0 几 n S t nF 阳m e c 比值0 ．2 7 V e l o c i t yi nS t nF r a m e P r o j e c t i o n V e C ∞r 1 4807e0005 3 4 1e 0 0 0 I Ⅲ、} ． 厂r ． a l C 值0 ．1 7 图1 0不同比值下的速度云图和流线图 V e I o c i l yi nS t nF r a m e P r o j e c t i o n V e c t o r l 1 75 ．3 4 1 e 0 0 0 U I I4 ．0 0 6 e 0 0 0 _ ■2 ．6 7 0 e 0 0 0 _ ■1 ，3 3 5 e 0 0 0 _ I O ．0 0 0 e 0 0 0 【ms “．1 】 d I C 值0 ．3 6 F i g ．1 0V e l o c i t yc o n t o u r sa n dc r o s ss e c t i o nf l o wc h a r t0 fd i f f e r e n tr a t i o 图1 1 为不同比值与输入功率之间的关系，从图 中可以看出随着面积比的增大，输入功率也随着增 加，且比值在0 ．1 5 t O ．3 5 左右时影响更为明显。 良 褂 髁 面积比 图1 1不同比值与输入功率之间的关系 F i g ．11R e l a t i o n s h i pb e t w e e nd i f f e r e n t r a t i oa n dp o w e r 4结论 针对嵌布粒度较细矿物对超细磨技术及设备的 要求，通过对现有搅拌装置类型进行分析，确定了适 用于超细磨技术及装备的盘式搅拌装置，利用C F D 仿真分析技术，明确了盘式搅拌装置的细磨机理，分 析了直径、转速、间距、面积比等多个关键参数对细 磨效果的影响规律，得到了各参数的取值范围，为实 验及工业技术及装备的研究与开发奠定了基础。 参考文献 E l i 卢世杰，孙小旭．大型立式螺旋搅拌磨机应用现状E J ] ．铜 业工程，2 0 1 4 2 3 8 4 2 ． L US h i ji e ，S U NX i a o x u ．A p p l i c a t i o no f l a r g ev e r t i c a l s p i r a ls t i r r i n gm i l l [ J ] ．C o p p e rE n g i n e e r i n g ，2 0 1 4 2 3 8 4 2 ． 下转第9 7 页 万方数据 2 0 2 0 年第1 期韩登峰等，气泡负载测量装置的创新及应用 9 7 E 5 ] 史帅星，韩登峰，张跃军，等．浮选柱内紊流强度对气泡与 颗粒碰撞概率的影1 1 i 向[ J ] ．有色金属 选矿部分 ，2 0 1 6 4 8 1 8 6 ． S H IS h u a i x i n g ，H A ND e n g f e n g ，Z H A N GY u e j u n ，e ta 1 ． E f f e c t so ft u r b u l e n c ei n t e n s i t yo nb u b b l e sa n dp a r t i c l e s c o l l i s i o np r o b a b i l i t yi nf l o t a t i o nc o l u m n [ J ] ．N o n f e r r o u s m e t a l s M i n e r a lP r o c e s s i n gS e c t i o n ，2 0 1 6 4 8 1 8 6 ． [ 6 ] L I UTY ，S C H W A R ZMP ．C F D - b a s e dm o d e l l i n go f b u b b l e p a r t i c l ec o l l i s i o ne f f i c i e n c yw i t h m o b i l e b u b b l e s u r f a c ei nat u r b u l e n te n v i r o n m e n t [ J ] ．I n t e r n a t i o n a l J o u r n a lo fM i n e r a lP r o c e s s i n g ，2 0 0 9 ，9 0 1 4 5 5 5 ． [ 7 ] Z H A N GW ．E v a l u a t i o no fe f f e c to fv i s c o s i t yc h a n g e so n b u b b l es i z e i nam e c h a n i c a lf l o t a t i o nc e l l [ J ] ．T r a n s a c t i o n s o fN o n f e r r o u sM e t a l sS o c i e t yo fC h i n a ，2 0 1 4 ，2 4 9 2 9 6 4 2 9 6 8 ． [ 8 ] Y I A N A T O SJB ，M O Y SMH ，C O N T R E R A SF ，e ta 1 ． F r o t hr e c o v e r yo fi n d u s t r i a lf l o t a t i o nc e l l s [ J ] ．M i n e r a l s E n g i n e e r i n g ，2 0 0 8 ，2 1 1 2 1 4 8 1 7 - 8 2 5 ． [ 9 ] N A K H A E IF ，M O S A V IMR ，S A MA ，e ta 1 ．R e c o v e r y a n dg r a d ea c c u r a t e p r e d i c t i o n o fp i l o t p l a n t f l o t a t i o n c o l u m nc o n c e n t r a t e N e u r a ln e t w o r ka n ds t a t i s t i c a l t e c h n i q u e s [ J ] ． I n t e r n a t i o n a lJ o u r n a lo fM i n e r a l P r o c e s s i n g ，2 0 1 2 1 1 0 1 1 1 1 4 0 1 5 4 ． E 1 0 ] 韩登峰，史帅星，吴峰，等．浮选过程差异性的气泡负载特 性对比分析[ J ] ．有色金属 选矿部分 ，2 0 1 7 增刊1 1 1 4 1 1 7 ． H A ND e n g f e n g ，S H IS h u i x i n g ，W UF e n g ，e ta 1 ． C o m p a r i s o na n dA n a l y s i so fB u b b l eL o a dC h a r a c t e r i s t i c s i nt h eP r o c e s so fF l o t a t i o n [ J ] ．N o n f e r r o u sm e t a l s M i n e r a lP r o c e s s i n gS e c t i o n ，2 0 1 7 S 1 1 1 4 1 1 7 ． [ 11 ] F A L U T S UM ，D O B B YGS ．F r o t hp e r f o r m a n c ei n c o m m e r c i a ls i z e d f l o t a t i o nc o l u m n s [ J ] ．M i n e r a l s E n g i n e e r i n g ，1 9 9 2 ，6 1 0 1 2 0 7 1 2 2 3 ． r 1 2 - ] S E A M A NDR ，F R A N Z I D I SJP ，M A N L A P I GEV ． B u b b l el o a dm e a s u r e m e n ti nt h ep u l pz o n eo fi n d u s t r i a l f l o t a t i o nm a c h i n e s an e wd e v i c ef o rd e t e r m i n i n gt h e f r o t hr e c o v e r yo fa t t a c h e dp a r t i c l e s [ J ] ．I n t e r n a t i o n a l J o u r n a lo fM i n e r a lP r o c e s s i n g ，2 0 0 4 ，7 4 1 1 - 1 3 ． [ 1 3 ] M O Y SMH ，Y I A N A T O SJ ，I ．A R E N A SJ ． M e a s u r e m e n to f p a r t i c l el o a d i n g o nb u b b l e si nt h e f l o t a t i o np r o c e s s [ J ] ．M i n e r a l sE n g i n e e r i n g 。2 0 1 0 ， 2 3 2 1 3 卜1 3 6 ． [ 1 4 ] Y I A N A T O SJ ，V I N N E T TL ，C A R R A S C OC ，e ta 1 ． E f f e c t o fe n t r a i n m e n ti nb u b b l el o a dm e a s u r e m e n to n f r o t h r e c o v e r ye s t i m a t i o n a ti n d u s t r i a ls c a l e [ J ] ． M i n e r a l sE n g i n e e r i n g ，2 0 1 5 ，7 2 3 1 3 5 ． 上接第8 1 页 [ 2 ] 孙小旭．G J M 型棒式搅拌磨机工业试验研究[ J ] ．有色金 属 选矿部分 ，2 0 17 3 6 6 6 9 ． S U NX i a o x u ．R e s e a r c ho nG J Mb a rs t i r r e dm i l lt h r o u g h c o m m e r c i a lt e s t [ J ] ．N o n f e r r o u sM e t a l s M i n e r a l P r o c e s s i n gS e c t i o n ，2 0 1 7 3 6 6 6 9 ． [ 3 ] R O I T T OI ，L E H T OH ，P A ZA ，e ta 1 ．S t i r r e dm i l l i n g t e c h n o l o g y - a n e w c o n c e p t i nf i n e g r i n d i n g [ J ] ． M e t a l l u r g i c a lP l a n tD e s i g na n dO p e r a t i n gS t r a t e g i e s M e t p l a n t2 0 1 3 ，2 0 1 3 1 1 9 0 2 0 1 ． [ 4 ] 何建成，卢世杰．基于A N S Y S 立磨机螺旋搅拌机构的分 析与优化[ J ] ．有色金属 选矿部分 ，2 0 1 4 3 7 2 7 5 ． H EJ i a n c h e n g ，I ．US h i j i e ．A n a y l s i sa n do p t i m i z a i t o no f v e r t i c a ls t i r r i n gm i l l ’Ss p i r a ls t i r r i n gm e c h a n i s mb a s e do n A N S Y S [ J ] ．N o n f e r r o u sM e t a l s M i n e r a lP r o c e s s i n g S e c t i o n ，2 0 1 4 3 7 2 - 7 5 ． [ 5 ] 张宏伟．介质搅拌磨单相流流场C F D 仿真模拟和实验验 证研究[ D ] ．广州华南理工大学，2 0 1 3 ． Z H A N GH o n g w e i ．C F DS i m u l a t i o na n de x p e r i m e n t a l v e r i f i c a t i o no fs i n g l ep h a s ef l o wf i e l d i nm e d i u ms t i r r e d m i l l [ D ] ．G u a n g z h o u S o u t hC h i n aU n i v e r s i t yo f T e c h n o l o g y ，2 0 1 3 ． [ 6 ] 张兆顺，崔桂香．流体力学[ M ] ．北京清华大学出版社， 2 0 0 6 3 1 3 3 1 4 ． Z H A N GZ h a o s h u n ，C U IG u i x i a n g ．F l u i dd y n a m i c s [ M ] ． B e i j i n g T s i n g h u aU n i v e r s i t yP r e s s ，2 0 0 6 3 1 3 3 1 4 ． 万方数据