离散单元法高炉无钟炉顶布料模拟研究.docVIP

离散单元法高炉无钟炉顶布料模拟研究 　　　　林成城　杜鹤桂 (宝山钢铁(集团)公司)　(东北大学) 摘　要　从动力学角度，采用离散单元法研究炉料颗粒的受力和运动，建立高炉无钟炉顶布料数值计算模型，模拟无钟炉顶布料过程和料面形状。模拟与实测对比表明，两者基本吻合，证明离散单元法在高炉无钟布料应用中可行，并为高炉无钟炉顶布料的机理研究开辟了一条新的途径。关键词　离散单元法　无钟炉顶　布料　料面形状　数值模拟 ON SIMULATION OF CHARGING IN BELL-LESS TOP BLAST FURNACE WITH DISTINCT ELEMENT METHOD　　　　LIN Chengcheng 　　　　　　　　　DU Hegui　　　　 (Baoshan Iron and Steel Corp.)　　(Northeastern University) ABSTRACT　In this paper,the distinct element method is adopted to study the force and the motion of granular burden from a kinematic view point.The numerical model of bell-less top charging is established.The charging process and burden surface profile have been imitated.The comparison made between simulation results and measured results,where the two results basically coincide with each other.It proves that the application of the distinct element method to blast furnace bell-less top charging is feasible which has also provided a new approach in theoretical study of bell-less top charging.KEY WORDS　distinct element method,bell-less top,charging,burden surface profile,numerical simulation 1　前言 　　高炉炉料都是散料体颗粒介质，在布料过程中表现出复杂特性，很难求解其运动方程的精确解，只能求其数值解，目前常用的数值计算方法对非连续介质有很大局限性。　　离散单元法由Cundall于1971年提出［1］，是研究散料体介质颗粒动力学的有力工具［2］，已在边坡工程，采矿工程和岩土力学等方面得到广泛应用。本文首次应用离散单元法建立炉料颗粒运动的数值模型，结合计算机实施技术，分别对高炉无钟炉顶单环布料和多环布料进行模拟 ，并对比模拟结果与实测结果，确立离散单元法在高炉无钟炉顶布料中应用的研究地位，对高炉无钟炉顶布料机理研究进行初步探索［3］。 2　离散单元法基本模型［4，5］ 2.1　基本原理　　离散单元法基本原理就是对松散介质的每一个单元，以牛顿第二运动定律为基础，在充分考虑单元的几何形状及其相邻单元相互作用关系前提下，结合不同的本构关系，运用动态松弛法进行循环迭代计算，运用各种形式的阻尼吸收单元的动能，按时步进行计算，并遍及整个介质单元，直至每一个单元不再出现不平衡力和力矩为止。2.2　基本假设　　(1) 把高炉炉料按大小分为若干等级，并简化成等当量直径的圆形准刚性颗粒。　　(2) 适用于二维平面问题研究。　　(3) 炉料颗粒之间只有摩擦力而无拉力。　　(4) 炉料颗粒接触模型是一弹性力学模型。如图1，用弹性元件Kn、Ks提供接触点的法向力和切向力；与刚度有关的粘性阻尼元件Dn,Ds吸收颗粒相对运动的动能；与质量和速度有关的阻尼元件Dm吸收颗粒的绝对动能。 图 1　颗粒接触力学模型　　Fig.1　Mechanical model for the interaction between particles 2.3　本构方程　　颗粒之间的法向力Fn与其之间的“叠合”位移un成正比，即 Fn＝Kn.un (1) 　　由于切向力与其运动和加载的途径有关，用其增量表示，切向力增量ΔFs与切向位移Δus增量成正比，即 ΔFs＝Ks.Δus (2) 式中，Kn、Ks分别为法向和切向刚度系数。 　　规定法向力Fn以使颗粒受压为正，切向力Fs以使颗粒产生逆时针旋转为正。如图2，