Simulation analysis of factors and secure thickness of horizontal pillar for backfilling
-
摘要: 为分析东同V号矿体充填体下水平矿柱的稳定性,采用AHP分析方法得出影响充填体下水平矿柱稳定性的因素权重排序,得出水平矿柱的厚度为最重要因素.利用小变形薄板理论的公式计算水平矿柱最小安全厚度为5 m,在此基础上,建立数值模型进一步模拟分析,对水平矿柱厚度为5 m至9 m进行应力、垂直位移及塑性区分析,最终得出,水平矿柱厚度小于7 m时,有破坏的危险,为矿山开采提供指导.Abstract: The factor weights for horizontal pillar stability were obtained to analyze the stability of horizontal pillar of the Dongtong' V ore by applying AHP analysis method. The results show that the thickness of the horizontal pillar is the most important factor. The minimum safety thickness of the horizontal pillars is 5 m by using the formula of small deformation of thin plate theory. A numerical model is established to make a further analysis. The stress, vertical displacement and plastic zone of horizontal pillar thickness of 5 m to 9 m are analyzed. It is concluded that the horizontal pillar is likely to collapse when its thickness is less than 7 m.
-
Keywords:
- AHP analysis /
- backfilling /
- horizontal pillars /
- secure thickness
-
电磁搅拌利用线圈激发的交变磁场渗透到金属熔体中, 在其中感生起感应电流, 感应电流与当地磁场互相作用产生电磁力。该搅拌方式不与坯体直接接触, 却能够有效的改善铸坯表面、皮下以及内部质量。目前在金属冶炼行业应用相对广泛的电磁搅拌器, 基本可以按照搅拌区域和激发磁场形态进行分类, 其中各模式磁场的电磁搅拌工艺针对不同的铸坯形状和材料特性有不同的效果。近年来, 许多钢铁企业、有色金属生产企业以及各大科研院所着力于开发新型的复合型材料, 这些材料在具有高性能的同时也对冶金工艺提出了新的要求[1-2]。单一磁场形态的电磁搅拌已经难以达到高性能复合材料的冶金要求, 由此湖南理工学院和湖南科美达电气股份有限公司联合设计开发了一款多模式磁场电磁搅拌器, 用来产生旋转磁场、行波磁场以及螺旋磁场, 可根据不同材料特性或铸坯形状, 在金属冶炼过程中选择不同的磁场形态进行电磁搅拌, 以求达到最佳的冶金效果。
电磁搅拌的过程受到磁流体力学和传热学等诸多领域的耦合作用[3-5], 是一个非常复杂的过程, 虽然多方面因素影响着电磁搅拌效果, 但整体而言搅拌磁场模式、励磁电流强度、励磁电流频率以及搅拌的时间是制约电磁搅拌效用的主要参数。为了研究分析各模式磁场的分布规律以及各主要工艺参数对磁场的影响, 本文利用有限元分析软件Ansoft Maxwell对多模式磁场电磁搅拌器建立了有限元数值分析模型, 主要对磁场中的磁感应强度和电磁力的分布特征以及数值变化规律进行研究, 以验证本设计的有效性, 并为后续设计提供参考。
1 多模式磁场电磁搅拌器的结构与原理
多模式磁场电磁搅拌器主体是由非磁性外壳、铁芯、旋转线圈绕组以及行波线圈绕组组成, 在铁芯和线圈之间加入了压圈和连接片等固定结构的零件。其装配图见图 1。共有12个线圈绕组, 其中6个线圈绕组垂直缠绕在铁芯上, 为线圈Ⅰ, 负责激发旋转磁场; 另外6个线圈绕组水平嵌于铁芯中, 为线圈Ⅱ, 负责激发行波磁场。2组线圈同时输入励磁电流, 则激发出螺旋磁场, 螺旋磁场的螺旋角度由旋转磁场线圈电流和行波磁场线圈电流的大小和方向共同确定。
2 基于Ansoft Maxwell的有限元分析
2.1 有限元仿真模型的建立
在保证计算精度的前提下, 为了提高模型的收敛性以及简化计算过程, 多模式磁场电磁搅拌器的计算模型采用如下假设条件:
1) 铁芯附近的磁感应强度分布均匀, 从截面穿出的磁力线垂直于截面;
2) 只计算铁芯表面部分对于搅拌器内腔磁场产生影响;
3) 结晶器内的金属为熔融状态;
4) 不考虑结晶器内流场变化的影响;
5) 结晶器与熔体之间有绝缘间隙。
以多模式磁场电磁搅拌器实验平台的结构参数为参考, 忽略细小结构件以及细小间隙构建三维几何模型, 以达到减小数值仿真模型的计算量和增加数值仿真模型收敛性的目的[6]。多模式磁场电磁搅拌器的结构参数如表 1所列, 结构模型如图 2所示。
表 1 电磁搅拌器结构尺寸与相关参数Table 1. Structure dimension and relevant parameters of electromagnetic stirrer2.2 磁场计算的数学模型
麦克斯韦方程组(Maxwell Equations)为描述电场、磁场与电流密度之间关系的方程组, 通过这个方程组可以描述出电磁场的规律。
(1) (2) (3) (4) 式(1)法拉第电磁感应定律, 说明了磁感应强度随时间的变化产生电磁场强度的环流; 式(2)安培定律, 表明传导电流密度可以产生磁场的环流, 其方向按右螺旋法则; 式(3)高斯定律, 表明磁力线皆为首尾相接的闭合曲线, 没有初始点, 也没有终止点[7-9]; 式(4)状态方程, 表明磁感应强度与磁场强度之间的关系; 结合总体方程组能表示出变化的磁场和电场之间相互作用形成稳态电磁场的基本原理[10-11]。利用这组方程理论上可以解决大部分宏观的电磁场问题。
2.3 材料属性定义与激励源设置
此分析为电磁场分析, 导磁区域为线圈、铁芯、空气、结晶器和金属熔体, 导电区域为线圈、结晶器和金属熔体, 材料属性定义见表 2。
表 2 有限元模型中材料的物理性质参数Table 2. Physical property parameters of materials in the finite element modelWD315-50硅钢片的相对磁导率的曲线如图 3所示。
图 3为WD315-50硅钢片在磁化过程中磁感强度(B)与磁场强度(H)之间关系的曲线, 当H逐渐增大时, B急剧增加, 几乎呈直线上升, 当H进一步增大时, B的增加开始变得缓慢, 到达一定值后, B几乎不再增加, 即达到了饱和。
对于多模式磁场电磁搅拌器, 在对旋转磁场线圈绕组和行波磁场线圈绕组同时通入三相交流电激励后, 在搅拌器内腔中激发出呈正弦变换的螺旋磁场, 如果我们减去旋转磁场型励磁线圈绕组的激励, 则得到行波磁场型的磁场形态。相反, 如果减去行波磁场型励磁线圈绕组的激励, 则可得到旋转磁场型的磁场形态[12]。若三相交流电之间的相位差为120°, 则施加至线圈截面上的电流激励可表示为:
(5) (6) (7) 式中:A1、A2、A3分别代表为第一相、第二相和第三相电流强度; A0为电流强度幅值; w为电源角频率, w=2πf; f为频率; t为时间[13-15]; 电磁场采用瞬态磁场计算方法, 时间步长设置为一个周期的1/20, 计算时长为两个周期。
3 多模式电磁搅拌数值模拟
在搭建完多模式磁场电磁搅拌器的有限元分析模型后, 通过改变电流幅值、电流频率和磁场类型等计算参数, 完成了多组不同计算参数的仿真计算。选取搅拌器内结晶器中心位置为原点, 以此选取结晶器内轴向和径向的磁感应强度数值和电磁力数值, 进行对比分析。
3.1 实测结果验证
对应多模式磁场电磁搅拌器中的旋转磁场和行波磁场, 利用高斯计对其结晶器中心轴线上以中心轴线坐标原点为中心等距离的5个点进行磁感应强度的测量, 之后利用数值仿真模型对相同位置的磁感应强度进行数值计算, 两者的结果对比如图 4、图 5所示。
从图 4、图 5中可以看出, 高斯计测量结果与数值模拟计算结果偏差较小, 两者基本一致, 说明所建数值分析模型是准确和可靠的[16]。由于螺旋磁场是由旋转磁场和行波磁场在空间中的叠加形成, 所以所建模型同样能对螺旋磁场进行准确可靠的数值模拟仿真。
3.2 各模式磁场矢量分析
选取电磁搅拌工艺参数为200 A、2 Hz, 分别对旋转磁场、行波磁场以及螺旋磁场进行数值仿真模拟, 截取结晶器中心位置横截面以及纵截面磁感应强度矢量图, 得出各模式磁场矢量分布特征, 见图 6、图 7和图 8。
由图 6可知, 旋转磁场磁感应强度呈现出的变换规律为成对变换, 电流随时间变化, 磁场波则按三相交流电相序沿圆周旋转, 从而导致金属溶液做旋转运动。
由图 7可知, 行波磁场磁感应强度呈现出的变化规律为环形纵向正弦变换且左右对称, 这种磁场分布使金属溶液的流速在纵向呈正弦变换, 从而形成一个左右对称的大环流, 流动中心与结晶器轴线重合。
由图 8可知, 螺旋磁场可以看成旋转磁场与行波磁场在空间上的叠加, 其方向由旋转磁场方向和行波磁场方向共同确定, 可以通过改变旋转磁场或者行波磁场的线圈电流的大小和方向来改变螺旋磁场的螺旋角[17-18]。
3.3 电流强度对磁场分布的影响
选取结晶器中轴线和结晶器中心横截面径向的磁感应强度数据, 电流频率固定为4 Hz, 分析各磁场模式下, 电流强度为50, 150, 200, 250 A时对应的磁感应强度分布情况。
图 9所示为搅拌频率为4 Hz时, 不同磁场模式下励磁电流强度对应结晶器中轴线以及径向各处磁感应强度值的影响。结果表明:随着励磁电流的增大, 不同磁场模式下的磁感应强度均有所增加, 不过磁场分布结构基本不改变; 相同的励磁电流下, 轴向分布呈现中心位置附近达到最大值, 向两侧逐渐减弱, 径向分布呈现出由中心向结晶器边缘逐渐增强的趋势[19-20]。
3.4 搅拌频率对磁场分布的影响
选取结晶器中轴线和结晶器中心横截面径向的实验数据, 分析各磁场模式下, 电流强度固定为200 A, 电流频率为0.5, 1, 2, 4, 6, 10 Hz时对应的磁感应强度分布如图 10所示。
图 10所示为励磁电流强度为200 A时, 不同磁场模式下结晶器中轴线以及径向各处磁感应强度随励磁电流频率的变化趋势。由图 10可见, 无论在哪种搅拌方式下, 在频率在0~2 Hz区间, 频率越大, 结晶器内的磁感应强度越大。而在2~10 Hz区间, 频率越大, 结晶器内的磁感应强度反而越小, 在2 Hz左右, 磁感应强度达到最大值。
3.5 电磁力计算结果分析
以前面定义的多模式电磁搅拌数值分析模型为例, 稀土铝熔体作为金属熔体, 选取200 A、2 Hz工艺参数进行仿真模拟, 计算结果分析如图 11和12所示。
由图 11可知, 铸坯的中心轴向电磁力在轴向上呈"两端小、中间大"的分布特征; 而铸坯在中心横截面上切向电磁力的径向分布, 呈现出从铸坯表面向中心缩减的特征, 铸坯中心受到的电磁力较小。
铸坯横截面电磁力随电流强度的变化曲线如图 12所示。可以看出, 电磁力的大小与电流强度大小呈正比关系; 电磁力的大小与电流频率的关系在0~3.5 Hz区间, 由于结晶器内磁感应强度有一个快速上升的过程, 且电流频率的增大导致感生电流渗入深度缩减至结晶器内有效工作区域, 导致电磁力大小随着频率增大而增大。而在3.5~10 Hz区间, 由于金属熔体受到集肤效应的影响[21], 熔体内的感生电流对金属熔体的渗入深度继续减小, 逐渐退出结晶器内有效工作区域, 搅拌器内磁感应强度也逐渐减弱, 随之产生的电磁力也随着频率增加而减小, 即电磁搅拌效率降低。
4 结束语
采用数值模拟仿真和实测相结合的方法, 对多模式磁场电磁搅拌器的磁场分布情况进行比较。对比分析了不同磁场模式下, 结晶器内的磁感应强度分布以及金属熔体所受到的电磁力情况, 得到如下结论:
1) 磁场分布特征:3种磁场类型作用下, 不同时刻的磁场分布形式相似, 只是磁场的相位发生了变化; 结晶器中心轴线磁感应强度数值的变化皆为"两头小, 中间大", 沿半径方向磁感应强度数值的变化则为结晶器边缘处最大, 向中心方向不断衰减, 与实际测量结果一致。
2) 电磁力分布规律:3种磁场类型作用下, 径向电磁力边缘最大, 逐渐向中心衰减; 轴向电磁力中下部较大, 沿轴线向两端逐渐减小。
3) 结晶器内磁感应强度在常规工况参数区间内, 与励磁电流强度成正比, 励磁电流强度越大则磁感应强度越大。磁感应强度与励磁电流频率的关系则表现为, 在0~2 Hz区间成正比关系, 在2~10 Hz区间成反比关系, 整个系统在2 Hz左右得到最大磁感应强度值。
4) 结晶器内金属熔体受到的电磁力大小与励磁电流强度成正比。在0~3.5 Hz区间, 频率越大, 结晶器内的磁感应强度越大; 而在3.5~10 Hz区间, 频率越大, 结晶器内的电磁力反而越小, 整个系统在3.5 Hz左右获得电磁力最大值。
-
表 1 充填体下水平矿柱稳定性影响因素综合评价体系
Table 1 Comprehensive evaluation system of stability of level under pillar filling body
目标层A 准则层B 指标层C 充填体下水平矿柱稳定性影响因素综合评价 地质条件B1 岩体结构C1 岩石物理力学性质C2 地应力C3 回采参数B2 采场暴露面积C4 间柱尺寸C5 水平矿柱厚度C6 工程条件B3 爆破震动影响C7 支护条件C8 上部充填体高度C9 表 2 AHP评价尺度
Table 2 AHP evaluation
阶数 内容 1 两个要素具有同等的重要性 3 认为一个要素稍微重要 5 强烈倾向于某一要素 7 非常倾向于某一要素 9 有证据确认袁某一要素比另一要素重要 2、4、6、8 上述标准的折中值 上述数值的倒数 若甲要素与乙要素比较袁被赋予某个标度值袁则乙要素与甲要素比较袁则被赋予此标度值的倒数 表 3 判断矩阵A-B
Table 3 Judgment matrix A - B
A-B 地质条件B1 回采参数B2 工程条件B3 地质条件B1 1 1/4 1/3 回采参数B2 4 1 1 工程条件B3 3 1 1 表 4 判断矩阵B1-C
Table 4 Judgment matrix B1-C
B1-C 岩体结构C1 岩石物理力学性质C2 地应力C3 岩体结构C1 1 1 2 岩石物理力学性质C2 1 1 1 地应力C3 1/2 1 1 表 5 判断矩阵B2-C
Table 5 Judgment matrix B2-C
B2-C 采场暴露面积C4 间柱尺寸C5 水平矿柱厚度C6 采场暴露面积C4 1 2 1/2 间柱尺寸C5 1/2 1 1/2 水平矿柱厚度C6 2 1 1 表 6 判断矩阵B3-C
Table 6 Judgment matrix B3-C
B3-C 爆破震动影响C7 上部充填体高度C8 支护条件C9 爆破震动影响C7 1 1 2 上部充填体高度C8 1/3 1 3 支护条件C9 1/2 1/3 1 表 7 判断矩阵计算结果
Table 7 Calculation results of judgment matrix
判断矩阵 A-B B1-C B2-C B3-C B1 B2 B3 C1 C2 C3 C4 C5 C6 C7 C8 C9 权重 0.126 3 0.457 7 0.416 0.411 0.328 0.261 0.311 9 0.197 6 0.490 5 0.387 3 0.442 9 0.169 8 特征向量 (0.126 3,0.457 7,0.416)T (0.411,0.328,0.261)T (0.311 9,0.197 6,0.490 5)T (0.387 3,0.442 9,0.169 8)T 一致性检验 λmax 3.009 0 3.054 0 3.053 7 3.018 3 C.I. 0.004 6 0.026 8 0.026 8 0.009 2 R.I. 0.520 0 0.520 0 0.520 0 0.520 0 C.R. 0.008 9 <0.1 0.051 6<0.1 0.051 7<0.1 0.017 6<0.1 表 8 随机一致性指标
Table 8 Random consistency index
阶数 1 2 3 4 5 6 7 R.I. 0 0 0.52 0.89 1.12 1.26 1.36 表 9 层次分析总权重排序
Table 9 Analytic hierarchy total weight sorting
A-B 地质条件B1 回采参数B2 工程条件B3 总权重 排序 B-C 0.126 3 0.457 7 0.416 0 1 / 岩体结构C1 0.411 1 / / 0.051 9 7 岩石物理力学性质C2 0.327 8 / / 0.041 4 8 地应力C3 0.261 1 / / 0.033 0 9 采场暴露面积C4 / 0.311 9 / 0.142 8 4 间柱尺寸C5 / 0.197 6 / 0.090 4 5 水平矿柱厚度C6 / 0.490 5 / 0.224 5 1 爆破震动影响C7 / / 0.387 3 0.161 1 3 上部充填体高度C8 / / 0.442 9 0.184 2 2 支护条件C9 / / 0.169 8 0.070 7 6 注:“/”指无对应参数. 表 10 计算模型材料力学参数
Table 10 Calculation model for material mechanics parameters
材料 内摩擦角/(°) 内聚力/MPa 泊松比 容重/(N·m-3) 弹性模量/GPa 抗拉强度/MPa 上盘围岩 30.0 5.988 0.24 27 500 37.200 2.675 下盘围岩 42.2 3.067 0.23 25 900 36.500 2.333 矿体 41.0 8.487 0.25 35 600 69.200 4.700 充填体 12.4 0.290 0.18 17 400 0.025 0.122 -
[1] 李其昌. 金阳公路保安矿柱开采技术研究[D].长沙:中南大学,2004. http://cdmd.cnki.com.cn/article/cdmd-10533-2004093624.htm [2] 罗小新,马玄恒,黄晓毅. 隔离矿柱合理参数的研究[J]. 重型机械,2010(增2):109-113. http://www.cnki.com.cn/Article/CJFDTOTAL-ZXJX2010S2030.htm [3] 宋华. 建筑物下保安矿柱回采地压分布规律与控制研究[D].武汉:武汉理工大学,2013. http://cdmd.cnki.com.cn/article/cdmd-10497-1013296131.htm [4] ZHOU J, LI X B, SHI X Z. Predicting pillar stability for underground mine using Fisher discriminant analysis and SVM methods [J]. Transactions of Nonferrous Metals Society of China, 2011(12):2734-2743. http://cn.bing.com/academic/profile?id=1999470637&encoded=0&v=paper_preview&mkt=zh-cn
[5] LIU Z X, DANG W G, HE X Q. Cancelling ore pillars in large-scale coastal gold deposit: A case study in Sanshandao gold mine, China [J]. Transactions of Nonferrous Metals Society of China, 2013(10): 3046-3054. http://cn.bing.com/academic/profile?id=1988503938&encoded=0&v=paper_preview&mkt=zh-cn
[6] 董凯程,吴春平,冯盼学. 充填体下水平矿柱开采技术研究[J]. 中国矿业,2012(9):103-105. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGKA201209031.htm [7] HU Y X, LI X B.Bayes discriminant analysis method to identify risky of complicated goaf in mines and its application [J]. Transactions of Nonferrous Metals Society of China, 2012(2): 425-431. http://cn.bing.com/academic/profile?id=2032602668&encoded=0&v=paper_preview&mkt=zh-cn
[8] 张炳江.层次分析法及其应用案例[M].北京:电子工业出版社,2014. [9] 李伯德. AHP判断矩阵一致性改进方法研究[D].兰州:兰州大学,2007. http://cdmd.cnki.com.cn/article/cdmd-10730-2008019065.htm [10] 唐晓玲,叶明亮. 薄煤层坚硬顶板的薄板理论分析及来压预报[J]. 矿山压力与顶板管理,2003,20(2):89-92. http://www.cnki.com.cn/Article/CJFDTOTAL-KSYL200302036.htm [11] 刘波. 复杂条件下水平矿柱稳定性分析及其开采方案研究[D].赣州:江西理工大学,2013. http://cdmd.cnki.com.cn/article/cdmd-10407-1014102035.htm [12] 黄敏. 铜绿山矿充填体下残矿回采关键参数数值模拟优化研究[D].长沙:中南大学,2012. http://cdmd.cnki.com.cn/article/cdmd-10533-1012477520.htm [13] 谢晓斌. 龙门山矿区L-23矿体上行式开采技术研究[D].赣州:江西理工大学,2012. http://cdmd.cnki.com.cn/article/cdmd-10407-1012507226.htm [14] 陈育民. FLAC及FLAC3D基础与工程实例 [M].北京:水利水电出版社,2013. [15] 蔡美峰. 岩石力学与工程 [M].北京:科学出版社,2009.