corex3000竖炉布料的离散元模型-炼铁技术炼钢技术.docVIP

COREX3000竖炉布料的离散元模型 李强1，冯明霞2，邹宗树1 (1．东北大学材料与冶金学院，辽宁沈阳110004；2．辽宁科技学院冶金工程学院，辽宁本溪117004) 摘 要：竖炉内煤气流的分布主宰着其内温度分布、煤气利用率和金属化率的高低，其上部调节方式仅有布料模式。通过引入离散颗粒动力学理论，基于经典牛顿力学和颗粒碰撞软球模型建立了针对COREX3000竖炉布料过程的离散元数值模拟模型，模型直观可视化再现装料过程，可定量获得料流轨迹、炉料落点及形成料堆的过程。应用模型进一步分析混装2种不同粒径颗粒，获得了炉料运动及形成的料堆过程，发现混装布料过程粒度偏析严重。建立的模型可为寻找和优化合理的布料模式提供重要研究基础。 关 键 词：COREX3000；竖炉；布料；离散元模型 COREX工艺是现阶段唯一工业化实践的熔融还原炼铁技术。对于COREX工业实践而言，控制竖炉内煤气流分布是COREX工艺生产中极为重要的操作技术之一，也是竖炉上部重要调剂手段。文献对竖炉布料模式的研究鲜有报道[1-3]，但关于高炉布料一直是研究热点[4-15]，可以说竖炉的布料与高炉布料既有类似又有不同。与高炉常用的溜槽布料相比，COREX竖炉布料器采用的是西门子公司研发的万向节布料器(Gimbal)，从功能上两者具有极大类似性。从布料过程的特点来看，高炉是矿焦分层布料，而竖炉是混合装入。 根据对布料过程研究方法的不同，主要可以分为4大类：物理模型模拟研究[4-7]、基于经典力学的料流质点轨迹的计算[8-12]、离散单元法的数值模拟[13-14]和现场布料实测[15]。其中离散单元法在处理颗粒物质工程应用领域发挥了其他数值算法不可替代的作用，离散单元方法最早于20世纪70年代由Cundall首先提出[16]，最初，它的研究对象主要是岩石等非连续介质的力学行为。采用离散单元法进行模拟分析，可以直接获得离散物质大量复杂行为信息和不易测量的颗粒尺度行为信息，能给出颗粒流宏观和微观力学的信息，对颗粒流的运动、受力、热量和质量传递提供有效的解决途径。离散单元方法已经成为最近几年的研究高炉布料和颗粒流动热点手段之一[13-14]。本文正是在这样的背景下提出基于离散颗粒动力学原理，建立COREX3000预还原竖炉布料过程的离散单元数值模型，研究竖炉顶部装料装置的特点及竖炉布料特性，为优化适宜于宝钢COREX3000预还原竖炉的布料模式提供模型基础。 1 离散元数学模型 在应用离散单元法进行数值模拟过程中，把炉料中的每个颗粒单独作为一个粒子单元建立数学模型，每个粒子单元自身都具有质量、尺寸、形状、刚度和阻尼等物理性质。各个粒子间存在接触和分离两种关系，每个粒子的本构关系体现于力与位移的关系，运动方程为经典的牛顿第二定律。在模拟过程中，通过相邻颗粒与颗粒间或者颗粒与壁面间的碰撞产生的接触力、力矩，计算每个颗粒在特定时刻的速度、加速度和位置等运动特征。通过每个颗粒的特征反映整个系统的运动特征。按处理接触方式的不同，在离散元中通常有“硬球”模型和“软球”模型之分。本文采用可同时考虑多个颗粒碰撞的软球模型离散单元法。 采用Tsuji提出的线性“弹簧─阻尼器─摩擦板”模型(见图1)来分析接触作用力。当颗粒i和颗粒j在发生碰撞时，首先在接触点处发生弹性变形。颗粒在前进方向上受到阻力，阻力的大小与法向变形位移、颗粒刚度成正比。达到最大位移变形时，颗粒停止运动，沿原来运动的方向反弹。碰撞时颗粒的动能会产生一定的损失，损失的大小取决于颗粒的弹性阻尼系数和颗粒间的相对速度。当2个粒子发生偏心碰撞时，碰撞点的接触力可分解为法向分力和切向分力，其中切向力的极值受到颗粒表面摩擦系数与法向力乘积的限制。当切向力大于该极值时，两颗粒在接触表面将发生滑动。同时颗粒间还存在滚动摩擦力矩。该力矩与滑动摩擦力类似，受到滚动摩擦系数的限制。 在从料仓经布料器进入竖炉过程中，颗粒主要受到2种力的作用，即自身重力mig及颗粒间和／或颗粒与壁面间的接触力。此外，颗粒i还受到2种力矩的作用，即切向力造成的力矩和滚动摩擦力矩。则根据牛顿第二定律，每个颗粒的平动和转动方程可描述为： 式(3)等号右侧第1个方括号内表示的是2个颗粒之间的法向力，第2个方括号内表示的是2个颗粒间的切向力。其中法向力包括2项，分别为接触弹力和法向阻尼力；切向力也包括2项，分别为切应力和切向阻尼力。发生滑动摩擦时(即︱Ft，ij︱μf︱F，ij︱)，切向力可采用如下公式来计算：Ft，ijμf︱F，ij︱ij。本文计算中忽略了颗粒与流体之间的作用力F。考虑到颗粒受到的唯一外力是重力Fmig。其中，法向和切向刚性系数取 2 结果与讨论 2．1 初始条件与边界条件 图2给出了竖炉DEM的计算区域几何和放大的布料器示意图。实验室具备IBMP55