充填管道布置方案优化中的数值模拟分析.pdf 4页VIP

2014年第8期／第35卷 黄 金 GoLD 充填管道布置方案优化中的数值模拟分析 羊羽翔，付建新，宋卫东 (北京科技大学土木与环境工程学院) 摘要：矿山在充填管道布置方案中参数的选取往往通过经验公式估算，而 目前尚无精确解的计 算方法。但参数的选取恰当与否直接影响到矿山充填系统的效率以及使用寿命，因此参数选取的 精确化便显得尤为重要。以石人沟铁矿的充填系统为工程背景，应用ANSYS分析软件中的FLOT— RANCFD模块建立管道模型，分析了管道 系统采用不同充填倍线时流速与压力的关系，并据此对 充填倍线的理想值进行了分析计算，最终得 出充填倍线为3．5，符合经验值的取值范围，对矿山充 填管道布置方案有一定借鉴意义。 关键词：FLOTRANCFD模块；充填倍线；管道模型；FLUD141单元；方案优化 中图分类号：TD853．34 文章编号：1001—1277(2014)08—0037—04 文献标志码 ：A doi：10．11792／hi201408O8 由于矿山管道输送工艺及其力学结构存在复杂 Y-／ 性，因此尚无对管道输送工程实际问题的精确解。目 前，通过环管试验或经验公式对两相流输送进行的有 a)FLUIDI41单元 b)FLUID142单元 关参数的计算已成为实际工程中最常采用的方法。 图1 FLOTRAN单元结构示意图 虽然环管试验得到的实测数据非常可靠，但是其试验 的开销过于昂贵且需要大量的矿样作为支撑，因此这 1)FLUID141单元特征。该单元为二维平面单元 种方法并没有普及开来…。近年来，随着计算机工 模型；呈四节点四边形或三节点三角形2种分布特 程的发展，通过对管道两相流的仿真分析和数值模 征；包括速度、温度、压力、素流动能、紊流能量耗散等 拟 J，对两相流参数的确定提供了行之有效的方 6种流体不同自由度，它们在综合质量中按照各 自不 法。采用大型有限元 ANSYS软件对充填管道进行相 同的质量所 占的份额也不同。 2)FLUID142单元特征。该单元为三维立体空间 关参数的数值仿真模拟便是其中之一，通过计算可以 得出流速、黏度及阻力损失等参数的变化与充填倍线 单元模型；呈四节点四面体或八节点六面体2种分布 之间的关系，从而对充填倍线数值的选择进行优化， 特征；包括速度、温度、压力、紊流动能、紊流能量耗散 为矿山充填系统选择相关参数提供依据。 等6种流体不同自由度，它们在综合质量中按照各 自 不同的质量所 占的份额也不同。 1 FLOTRANCFD简介 由于充填管道具有对称性，为了简化计算，充填 FLOTRANCFD是ANSYS分析软件中的一种模 管道一般采用FLUD141的平面单元进行模拟即可。 块，是用来分析流体流动性的工具，且不论对二维还 2 工程背景 是三维问题的解决都卓有成效，其中也包括一些弯管 中流体 的复杂 的三维 问题 J。它包括 FLUD141、 河北钢铁集团矿业有限公司石人沟铁矿 (下称 “石人沟铁矿”)位于河北省遵化市西北 10km，东南 FLUD1422种单元形式，这2种单元形式可以计算黏 性流体的各种参数分布情况，包括二维和三维流动、 距唐山市90km。石人沟铁矿 2001年开始转人地下 压力和温度分布等参数。 开采，根据三期工程设计，矿山设计规模 200万t／a， FLUID141、FLUID1422种单元形式在计算黏性 矿山最终产品为单一铁精矿，每年产出尾矿 130万t。 流体的各种参数分布情况有不同的形式和使用情况。 石人沟铁矿已建成 1座充填站，由中钢集团马鞍 山矿院工程勘察设计有限公司设计施工总承包。充 2种单元形式的特征如图1所示。 填站内设置4套充填系统，每套充填系统可以满足 收稿 日期 ：2014—05—30 基金项 目：国家 自然科学基金项 目；博士点基金项 目(20120006110022) 作者简介：羊羽翔(199l一)，男 ，江苏靖江人，硕士研究生，从事金属矿11I果矿、充填工艺研究；北京市海淀区学院路3O号，北京科技大学土木与环 境工程学院373信箱，100083 黄 金 50万t／a开采规模的要求，同时每套充填系统有2个 式中： 为雷诺数；p为流体密度(kg／m )；f为物体 砂仓、1个水泥仓。充填体的料浆浓度为72％。密 特征长度(m)；为特征速度(m／s)，如对圆管流动取 度为1570kg／m ，动力黏度为0．118Pa·s。矿山使 管内平均速度 ；r／为动力黏度(Pa·s)。 用的管材直径为 121mm。根据管径和充填能力，料 模型相似简化处理后，其管路长度见表 2。管径 浆的流速确定为3m／s。 不变D=121mm，垂直管道长度为25m。 表2 充填倍线与管长对应关系 3 模型建立与求解 充填倍线 Ⅳ 3．0 3．5 4．0 4．5 5．0 5．5 6．0 大型通用有限元分析软件ANSYS中采用FLOT· 管道总长／m 75 87．5 100 112．5 125 137．5 150 RAN模块分析步骤 以及各个步骤之间的关系如图2 所示 。 3．3 载荷的施加 自流输送充填料浆的载荷施加及分析，主要考虑 4个方面：给定浆体在管道中的进 口速度；对管道壁 施加无滑动位移边界条件；管道出口气压与管道进 口 的大气压相等，即相对气压为零；重力载荷，即在垂直 负方向上的重力加速度。在计算中施加的载荷将不 随单元方向变化而改变，始终保持它们最初的方向， 表面载荷作用在单元表面的法向。 3．4 模型求解结果 建立管路 自流输送模型，并施加载荷运算，各 自 由度收敛监测值如图3所示。 图2 FLOTRAN模块 7个主要步骤之间的相互关系 3．1 模拟流体的条件假设 在建立优化模型时，对流体以下几个方面赋予较 为理想化的处理 (对模拟结果影响不大)：黏性浆体 具有恒黏性，不随温度、时间的变化而变化；浆体为宾 汉体；不考虑热交换；不考虑震动、地压波等对管道输 送的影响。 3．2 建立简化模型 本次优化模拟采用充填管道的简化模型。假设 垂直管道长为 250m，且始终保持不变，管路总长随 图3 各 自由度收敛监测 图 着充填倍线 (Ⅳ)的改变而改变。充填倍线为 3．0～ 6．0，管路总长见表 1。 最终得到各充填倍线满管流状态下的流速变化 表 1 充填倍线与管长对应关系 图、压力变化图，以N=3．0为例，如图4所示。 充填倍线 N 3．0 3．5 4．0 4．5 5．0 5．5 6．0 4 充填管路方案的数据统计与分析 管道 总长／m 750 875 1000 1125 1250 1375 1500 由于网格划分的尺度有细微偏差，尤其在弯管处 的划分会有不同，故软件计算的精确度会有细小变 考虑到按照实际管长来模拟网格划分难度大、计 化，但结果的整体性和趋势还是比较明显的，对分析 算容易不稳定、结果说服力差等原因，本研究拟采用 不会造成很大的影响。 雷诺数相似准则公式 对模型相似处理，对于本文 根据模拟的过程，出口压力设置为0，所以进 口 分析的管路系统，JD、f、叩几个参数是不变的，给定特 压力即为整个管道输送时的阻力损失值 。进 口压力 征初速度相同。对于管路垂直高度和水平距离均缩 随着充填倍线的增大而增大，说明管路阻力的大小也 小 1O倍，模拟效果不变。 ’ 是与充填倍线成正比(如图5和表3所示)。损失的 ： f1 压力增大会导致流体受到管壁的摩阻加大，即管壁的 2014年第 8期／第35卷 磨损增加，这不仅会减少充填系统的使用时长，也会 弯管处的压力损失是选择充填倍线时需重点考察的 增加充填系统的使用成本。但当倍线太小时，依据前 因素。由图5以及表 3中可以看出，当充填倍线N= 文分析浆体对管道冲击压力增大，使得弯管局部压力 3．5时，进 口压力最小，也即阻力最小，作为优化值考 增大，也即增加了管道的磨损。由于竖直管道与水平 虑可行。 管道连接处是整个充填管道系统最脆弱的地方，所 以 速度值，佃 ·s) 。 3 216．~ 3 蝴 川 33l9抛 a)管道简化模型 b)弯道处流速变化 ^ ．Ⅲ)、缎 口哥 3 2 2 2 l l l O O O l 6 3 O 7 4 l 8 5 2 2 6 6 7 7 8 8 8 9 9 4 4 8 2 6 O 4 8 2 6 O 速度值 ／(m ·) 0 0．669576 133915 2-加873 2 83 0．334788 1．00436 1,67394 2．M351 3．3192 出口处宽度 ／cm c)出口处流速变化 d)出口处流速分布 图4 N：3．0时管道模型及流速变化图 表 3 充填倍线与进 口压力对应关系 进 口初始速度 0／(m ·S ) 3．O 3．0 3．O 3．O 3．O 3．0 3．O 出口速度 ／(m-S ) 3．124 2．965 4．065 3．097 5．625 4．668 5．173 沿程速度最大值 vma~／(m ·s ) 3．319 3．418 4．268 4．243 6．512 6．584 5．941 速度变化量 Av／(m ·s ) 0．319 O．418 1．268 1．243 3．512 3．584 2．94l 1 根据模拟所得数据，在初始速度相同的条件下， 要i 不同充填倍线下出口速度、沿程速度最大值及速度变 管 化量如表4和图6所示。 1)充填倍线与出口速度。随着充填倍线的增 加，出口处速度的最大值呈现出增大的趋势。而出口 充填倍线 最大速度增大，会明显引起管 口的振荡，使料浆输送 不稳定，直接影响充填体的质量。 2)充填倍线与沿程最大速度。随着充填倍线的 田 ● 增加，沿程最大速度也呈现增大的趋势。而且变化量 旦 、 较大，这样使得系统在局部的损失增大，磨损相应加 强，容易破管，从各种充填倍线下的速度变化图可以 明显看出，在弯管处的变化最大，也最容易破坏。沿 程速度变化大，同时给系统运行造成不稳定因素。在 图6 充填倍线与速度示意图 充填倍线 Ⅳ为3或3．5时，变化量较小，可以作为优 黄 金 化值考虑。 充填 自流输送管路系统理想充填倍线为3．5。对石 通过模拟和分析管路，锁定充填倍线 (Ⅳ)为优化 人沟铁矿在充填管路布置方案的选择上有一定的借 目标。一般工程经验认为，充填倍线 Ⅳ在 3．0—6．0， 鉴意义。 管道能够进行 自流输送。因此，对 J7v连续取值 ，并建 4)ANSYS分析软件在充填管道布置方案优化中 立模型模拟和分析。根据 Ⅳ对压力的影响，取得优 的重要作用十分明显。FLOTRANCFD模块的应用 化值 N=3．5。根据Ⅳ对速度的影响，取得优化值 较好地解决了管道输送工程实际问题的精确解难求， N=3和 N=3．5。综合考虑，石人沟铁矿充填 自流输 只能粗略给出范围的问题。通过对变量的连续取值， 送管路系统理想充填倍线为3．5。 并建立模型模拟和分析可以求得相对精确的解答。 但为了简化计算，模型中做了许多理想性假设。在工 5 结 论 程真实管道中的实际情况要复杂得多，这些假设条件 1)由数值模拟得出的结论可知，进 口压力的值 对最终解的影响还需进一步考量。 随着充填倍线的增加先减小后增大。充填倍线太小 [参 考 文 献] 时，浆体对弯管处的冲击压力较大，导致进 口压力较 [1] 李鹏程，韩文亮，田龙．高浓度管道输送参数计算模型的研 究[J]．金属矿山，2005(4)：60—62． 大。充填倍线增加时，流体受到管壁的摩阻加大，导 [2] 王新民，丁德强，吴亚斌，等．膏体充填管道输送数值模拟与分 致进 口压力增大。 析[J]．中国矿业 ，2006，15(7)：57—59． 2)随着充填倍线的增加，出口处速度的最大值 [3] 王海瑞．全尾砂胶结充填 自流输送管路改造及优化 [D]，长沙： 呈现出增大的趋势，沿程最大速度也呈现增大的趋 中南大学 ，2010． 势。而出口速度增加会导致管 口震荡，沿程最大速度 [4] 丁德强，王新民．膏体管道输送阻力损失研究[J]．矿业快报， 2007，23(1)：32—34． 增加会导致管道磨损加剧。 [5] 鲍鹏，岳建伟．流体力学[M]．郑州：黄河水利出版社，2006． 3)通过后期数据的处理分析，得出石人沟铁矿 Numericalsimulationinanalysisoftheoptimizationoffillingpipearrangements YangYuxiang，FuJianxin，SongWeidong (SchoolofCivilandEnvironmentalEngineering，BengUniversityofScienceandTechnology) Abstract：Selectionofparametersinthelayoutoffillingpipelinesisoftentheestimationofempiricalformulain- steadofexactalgorithm thathasnotyetbeensetup．However．itisimportanttoensurepreciseparametersSOthatfill— ingsystemscanbeusedefficientlyandcontinuously．Takethefillingsystem inShirengouIronMinefortheengineering background，thepipelinemodelhasbeenestablishedbytheapplicationofFLOTRAN CFD moduleinANSYS．There— lationbetweenvelocityandpressurehasbeenanalyzedindifferentfillingfoldlinesandtheresultleadstothecalcula- tionoftheidealvalueofthefillingfoldline，whichis3．5．TheresultgoeswiththevMuerangeinferredbyexperience， andalsohascertainref