高压水射流冲击刚壁压力分布规律研究.pdf

高压水射流冲击刚壁压力分布规律研究 ① 刘智谋， 毛桂庭， 柯 熠， 欧阳邓培 （长沙矿冶研究院 深海矿产资源开发利用技术国家重点实验室，湖南 长沙 ４１００１２） 摘 要 为了探讨高压纯水射流和高压磨料射流冲击刚壁时射流冲击压力的分布规律，利用建立的专用测力工作台，研究了射流 方向与刚壁成 ９０角时冲击压力的分布规律及靶距、系统工作压力对冲击压力的影响，并对磨料射流的强化效果进行了分析。 结果 表明射流冲击力随径向距离的增加急速降低，是一非单调递减阶梯函数，连接各阶中心冲击力值，可视其为具非单调负指数下降 规律。 在靶距不变的情况下，纯水射流和磨料射流，其中心（最高压力）位置的冲击压力和系统压力呈正相关。 关键词 高压水射流； 冲击力； 分布规律； 表面清理 中图分类号 ＴＰ６０６文献标识码 Ａｄｏｉ１０．３９６９／ ｊ．ｉｓｓｎ．０２５３－６０９９．２０１６．０５．００４ 文章编号 ０２５３－６０９９（２０１６）０５－００１３－０４ Ｅｘｐｅｒｉｍｅｎｔａｌ Ｒｅｓｅａｒｃｈ ｏｆ Ｉｍｐａｃｔ Ｆｏｒｃｅ Ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｎ Ｒｉｇｉｄ Ｗａｌｌ ｂｙ Ｈｉｇｈ Ｐｒｅｓｓｕｒｅ Ｗａｔｅｒｊｅｔ ＬＩＵ Ｚｈｉ⁃ｍｏｕ， ＭＡＯ Ｇｕｉ⁃ｔｉｎｇ， ＫＥ Ｙｉ， ＯＵＹＡＮＧ Ｄｅｎｇ⁃ｐｅｉ （Ｓｔａｔｅ Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ ｏｆ Ｅｘｐｌｏｒａｔｉｏｎ ａｎｄ Ｕｔｉｌｉｚａｔｉｏｎ ｏｆ Ｄｅｅｐ⁃ｓｅａ Ｍｉｎｅｒａｌ Ｒｅｓｏｕｒｃｅｓ， Ｃｈａｎｇｓｈａ Ｒｅｓｅａｒｃｈ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｍｉｎｉｎｇ ＆ Ｍｅｔａｌｌｕｒｇｙ Ｃｏ Ｌｔｄ， Ｃｈａｎｇｓｈａ ４１００１２， Ｈｕｎａｎ， Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ Ｆｏｒ ｅｘｐｌｏｒｉｎｇ ｔｈｅ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｉｍｐａｃｔ ｆｏｒｃｅ ｏｎ ｒｉｇｉｄ ｗａｌｌ ｂｙ ｐｕｒｅ ｗａｔｅｒｊｅｔ ａｎｄ ａｂｒａｓｉｖｅ ｊｅｔ， ｔｈｅ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｉｍｐａｃｔ ｆｏｒｃｅ ｗｉｔｈ ｗａｔｅｒｊｅｔ ｃｕｔｔｉｎｇ ａｎｇｌｅ ａｔ ａ ９０ ａｎｇｌｅ ｔｏ ｔｈｅ ｒｉｇｉｄ ｗａｌｌ， ａｎｄ ｅｆｆｅｃｔｓ ｏｎ ｔｈｅ ｉｍｐａｃｔ ｆｏｒｃｅ ｂｙ ｔｈｅ ｄｉｓｔａｎｃｅ ｏｆ ｎｏｚｚｌｅ ｆｒｏｍ ｔｈｅ ｒｉｇｉｄ ｗａｌｌ ａｎｄ ｗａｔｅｒ ｐｒｅｓｓｕｒｅ ａｖａｉｌａｂｌｅ ｆｏｒ ｔｈｅ ｗａｔｅｒｊｅｔ ｓｙｓｔｅｍ ｗｅｒｅ ａｌｌ ｔｅｓｔｅｄ ｗｉｔｈ ａ ｓｐｅｃｉｆｉｃ ｆｏｒｃｅ⁃ ｍｅａｓｕｒｅｍｅｎｔ ｔａｂｌｅ， ａｎｄ ｔｈｅ ｉｎｔｅｎｓｉｆｙｉｎｇ ｅｆｆｅｃｔ ａｃｈｉｅｖｅｄ ｂｙ ａｂｒａｓｉｖｅ ｊｅｔｔｉｎｇ ｗａｓ ａｌｓｏ ａｎａｌｙｚｅｄ． Ｉｔ ｉｓ ｆｏｕｎｄ ｔｈａｔ ｉｍｐａｃｔ ｆｏｒｃｅ ｏｆ ｗａｔｅｒｊｅｔ ｄｅｃｒｅａｓｅｓ ｒａｐｉｄｌｙ ｗｉｔｈ ｔｈｅ ｉｎｃｒｅａｓｉｎｇ ｏｆ ｒａｄｉａｌ ｄｉｓｔａｎｃｅ， ａｓ ａ ｎｏｎ⁃ｍｏｎｏｔｏｎｉｃａｌｌｙ ｄｅｃｒｅａｓｉｎｇ ｓｔｅｐ ｆｕｎｃｔｉｏｎ， ｃｏｎｎｅｃｔｉｎｇ ｔｈｅ ｃｅｎｔｒａｌ ｐｏｓｉｔｉｏｎ ｏｆ ｗｏｒｋｂｅｎｃｈ ａｎｄ ｆｏｌｌｏｗｉｎｇ ａ ｎｏｎ⁃ｍｏｎｏｔｏｎｉｃ ｎｅｇａｔｉｖｅ ｅｘｐｏｎｅｎｔｉａｌ ｄｅｃｒｅａｓｉｎｇ ｒｕｌｅ． Ｗｉｔｈ ｔｈｅ ｓａｍｅ ｄｉｓｔａｎｃｅ ｂｅｔｗｅｅｎ ｔｈｅ ｎｏｚｚｌｅ ａｎｄ ｒｉｇｉｄ ｗａｌｌ ｆｏｒ ｂｏｔｈ ｐｕｒｅ ｗａｔｅｒｊｅｔ ａｎｄ ａｂｒａｓｉｖｅ ｊｅｔ， ｔｈｅ ｃｅｎｔｒａｌ ｐｏｓｉｔｉｏｎ （ｗｉｔｈ ｔｈｅ ｈｉｇｈｅｓｔ ｐｒｅｓｓｕｒｅ） ｃｏｒｒｅｌａｔｅｓ ｐｏｓｉｔｉｖｅｌｙ ｗｉｔｈ ｔｈｅ ｓｙｓｔｅｍ ｐｒｅｓｓｕｒｅ． Ｋｅｙ ｗｏｒｄｓ ｈｉｇｈ ｐｒｅｓｓｕｒｅ ｗａｔｅｒｊｅｔ； ｉｍｐａｃｔ ｆｏｒｃｅ； ｄｉｓｔｒｉｂｕｔｉｏｎ； ｓｕｒｆａｃｅ ｃｌｅａｎｉｎｇ 带钢表面氧化铁皮的酸洗去除给环境带来严重污 染，近年来国内外发展多种无酸除鳞方法来代替酸洗 除鳞。 在带钢表面无酸除鳞的众多方法中，利用高压 水混合磨料去除带钢表面氧化铁皮（又叫高压水喷丸 除鳞） ［１］，由于去除效率高，无废气废水排放，被广泛 推广应用。 从前对于射流特性的研究，多集中于水刀切割、岩 石破碎、光学材料抛光等小口径喷嘴射流［２－３］。 小口 径喷嘴射流要求射流能量更好地集中，从而有利于材 料的切割去除。 对于应用基本段［４］进行表面材料去 除的大口径喷嘴射流特性，目前尚未进行充分研究。 大口径喷嘴射流能量的高度集中并不能提高对于表面 材料的清理效率，反而可能影响被清理工件的表面质 量（出现条状凹坑），其射流的压力分布理想情况为降 低射流中心位置峰值，扩大在某一能量（压力值）以上 的分布区域。 在系统压力和其他条件一致的情况下， 射流分布中，位于某个能量临界值（刚好能够产生表 面材料去处并达到工业生产所需表面质量的冲击压 力）范围越广，喷头对于表面清理的效率越高，能量利 用率亦越高。 本文研究了特定喷头下射流冲击刚壁时 材料表面冲击压力的分布规律，系统工作压力、靶距对 冲击压力的影响及磨料射流的强化效果，为高压水射 流清理系统的参数优化、评价清理性能和后续改进设 计提供依据。 ①收稿日期 ２０１６－０２－１６ 作者简介 刘智谋（１９９０－ ），男，湖南长沙人，硕士研究生，主要研究方向为采矿工程和高压水射流理论。 第 ３６ 卷第 ５ 期 ２０１６ 年 １０ 月 矿矿 冶冶 工工 程程 ＭＩＮＩＮＧ ＡＮＤ ＭＥＴＡＬＬＵＲＧＩＣＡＬ ＥＮＧＩＮＥＥＲＩＮＧ Ｖｏｌ．３６ №５ Ｏｃｔｏｂｅｒ ２０１６ 万方数据 １ 高压射流冲击刚壁的流场分布 高压射流清理材料表面过程的流场按其流动特性 可分成 ３ 个区域自由射流区、冲击区和壁面射流区， 如图 １ 所示［５］。 自由射流区内，射流受壁面的影响非 常小，因此可认为该区域的流动特性与自由射流相同。 在冲击区，射流急速改变方向，由轴向流动变为径向流 动并且存在极大的压力梯度。 该区域是引起壁面发生 弹塑性变形，甚至塑性去除、脆性去除或断裂的集中区 域。 壁面射流区的流动速度减弱，而且相对平稳，因此 一般不会发生壁面材料变形或去除［６］。 r P 7*2 D A49 41 2 *09 649 τ 图 １ 高压射流冲击刚壁的流场分布 ２ 高压水射流试验系统 为研究高压水射流与高压磨料射流对刚壁冲击力 的分布规律，设计了高压水射流测力系统，见图 ２。 PC, 02 ; 5-; A42 4 14, 414, A4, 88 .5, ; 0 , , D9 7;, 7E 4;                              30 MPa 35 MPa 40 MPa 45 MPa 50 MPa 55 MPa 图 ５ 磨料射流强化系数对比 从图 ５ 可以看出，系统工作压力不大于 ５０ ＭＰａ、靶 距 １７０ ｍｍ 和系统压力大于 ５０ ＭＰａ、靶距 ２００ ｍｍ 的情 况，具有较高的强化系数。 其中系统压力 ４５ ＭＰａ、靶距 １７０ ｍｍ，系统压力 ５０ ＭＰａ、靶距 １７０ ｍｍ 和系统压力 ５５ ＭＰａ、靶距 ２００ ｍｍ 达到试验参数的最优值。 图 ６ 给出了在不同靶距情况下的平均强化系数随 系统工作压力的变化关系。 4MPa       0.4 0.3 0.2 0.1 0.0 3025354045505560 8/; 图 ６ 平均强化系数随系统工作压力的变化关系 ５１第 ５ 期刘智谋等 高压水射流冲击刚壁压力分布规律研究 万方数据 在绝大部分压力段（除４０～４５ ＭＰａ 范围），平均强 化系数随系统工作压力增大而增大，磨料强化效能与 系统工作压力呈正相关关系。 ６ 结 语 １） 建立了专用射流测力工作台，利用不同直径圆 盘承受冲击力的差值得出了纯水和磨料射流冲击刚壁 的压力分布。 磨料射流除冲击压力显著增加外，和纯 水射流具有相近的压力分布。 射流冲击力随径向距离 增加急速降低，是一非单调递减阶梯函数，可视其为具 非单调负指数下降规律，而非切割小口径喷嘴射流形 成的高斯分布规律。 ２） 对系统工作压力和靶距对冲击力的影响研究 表明在靶距不变的情况下，无论纯水射流还是磨料射 流，其中心（最高压力）位置的冲击压力和系统压力呈 正相关。 纯水射流在靶距 １４０ ｍｍ 时，冲击压力明显 大于其他 ３ 个靶距（１７０，２００ 和 ２３０ ｍｍ）的冲击压力； 而对于磨料射流，冲击压力随靶距增大，近似等速率降 低，呈负相关变化。 ３） 磨料的加入对射流冲击压力有显著的强化效 果。 这是由于磨料射流为两相流冲击，其密度远大于纯 水射流的缘故。 系统压力不大于５０ ＭＰａ、靶距１７０ ｍｍ 和系统压力大于 ５０ ＭＰａ、靶距 ２００ ｍｍ 的情况，具有较 高的强化系数。 其中系统压力 ４５ ＭＰａ、靶距 １７０ ｍｍ， 系统压力 ５０ ＭＰａ、靶距 １７０ ｍｍ 和系统压力 ５５ ＭＰａ、 靶距 ２００ ｍｍ 达到试验参数的最优值，可以认为，在这 些情况下，具有较佳的能量利用率和射流功效。 ４） 在绝大部分系统工作压力区段，平均强化系数 随工作压力增大而增大，其强化效能与系统工作压力 呈正相关关系。 参考文献 ［１］ Ｋｅｖｉｎ Ｖｏｇｅｓ，Ａｌａｎ Ｍｕｅｔｈ． Ｅｃｏ⁃ｐｉｃｋｌｅｄ ｓｕｒｆａｃｅ ａｎ ｅｎｖｉｒｏｎｍｅｎｔａｌｌｙ ａｄｖａｎｔａｇｅｏｕｓ ａｌｔｅｒｎａｔｉｖｅ ｔｏ ｃｏｎｖｅｎｔｉｏｎａｌ ａｃｉｄ ｐｉｃｋｌｉｎｇ［Ｊ］． Ｉｒｏｎ ＆ Ｓｔｅｅｌ Ｔｅｃｈｎｏｌｏｇｙ，２００８，５（８）８１－９６． ［２］ 侯亚康，毛桂庭，阳 宁． 淹没磨料水射流对花岗岩的冲蚀研究 ［Ｊ］． 矿冶工程， ２０１１，３１（３）１８－２１． ［３］ 薛胜雄． 高压水射流技术及应用［Ｍ］． 北京机械工业出版社， １９９５． ［４］ 刘新泰． 高压磨料射流切割的仿真研究［Ｄ］． 包头内蒙古科技大 学机械工程学院， ２０１２． ［５］ Ｂ Ｎｉｅ， Ｈ Ｗａｎｇ， ＬＩ Ｌｅｉ， ｅｔ ａｌ． Ｎｕｍｅｒｉｃａｌ ｉｎｖｅｓｔｉｇａｔｉｏｎ ｏｆ ｔｈｅ ｆｌｏｗ ｆｉｅｌｄ ｉｎｓｉｄｅ ａｎｄ ｏｕｔｓｉｄｅ ｈｉｇｈ⁃ｐｒｅｓｓｕｒｅ ａｂｒａｓｉｖｅ ｗａｔｅｒｊｅｔ ｎｏｚｚｌｅ［Ｊ］． Ｐｒｏｃｅｄｉａ Ｅｎｇｉｎｅｅｒｉｎｇ， ２０１１，２６（４）４８－５５． ［６］ 李兆泽． 磨料水射流抛光技术研究［Ｄ］． 长沙国防科学技术大学 研究生院， ２０１１． 􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎 （上接第 １２ 页） ２） ｗ／ ｒ ＝ １，可作为集矿车设计或运行时的参考 指标。 参考文献 ［１］ 简 曲． 我国已开始对大洋富钴结壳资源的开发研究［Ｊ］． 矿业 研究与开发， １９９９，１８（１）５０－５６． ［２］ 曹 颖． 加速国际海底区域资源开发产业化的战略思考［Ｊ］． 海 洋开发， ２００３（１）５７－６０． ［３］ 何清华，李爱强，邹湘伏． 大洋富钴结壳调查进展及开采技术［Ｊ］． 金属矿山， ２００５，４３（５）４－７． ［４］ Ｙａｍａｚａｋｉ Ｔ， Ｐａｒｋ Ｓ Ｈ， Ｓｈｉｍａｄａ Ｓ， ｅｔ ａｌ． Ｄｅｖｅｌｏｐｍｅｎｔ ｏｆ Ｔｅｃｈｎｉｃａｌ ａｎｄ Ｅｃｏｎｏｍｉｃａｌ Ｅｘａｍｉｎａｔｉｏｎ Ｍｅｔｈｏｄ ｆｏｒ Ｃｏｂａｌｔ⁃Ｒｉｃｈ Ｍａｎｇａｎｅｓｅ Ｃｒｕｓｔｓ［ Ｃ］ ∥ Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｏｆｆｓｈｏｒｅ ａｎｄ Ｐｏｌａｒ Ｅｎｇｉｎｅｅｒｉｎｇ Ｃｏｎｆｅｒｅｎｃｅ， ２００２４５４－４６１． ［５］ 倪建宇，周怀阳． 中国多金属结核开辟区的深海环境［Ｊ］． 海洋地 质与第四纪地质， ２００２（２２）１１２－１１６． ［６］ 梁 平，石海林，崔 波，等． 洋底富钴结壳的开采方法［Ｊ］． 金属 矿山， ２００２，１４（２）５３－６０． ［７］ 国家海洋局科技公司． ＤＹ９５⁃１１ 航次现场总结报告［Ｒ］． 中国大 洋矿产资源协会， １９９９． ［８］ 中国大洋矿产资源勘察． ＤＹ１０５⁃１１ 航次现场总结报告［Ｒ］． 中国 大洋矿产资源研究开发协会， ２０００． ［９］ Ｚｈｅｎｇ Ｘ Ｍ， Ｓｏｎａｌ Ｊ， Ｓｖｅｎ Ｋ， ｅｔ ａｌ． Ｍｕｌｔｉ⁃ｒｏｂｏｔ ｆｏｒｅｓｔ ｃｏｖｅｒａｇｅ［Ｃ］∥ Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ＩＥＥＥ Ｉｎｔ Ｃｏｎｆ ｏｎ Ｉｎｔｅｌｌｉｇｅｎｔ Ｒｏｂｏｔｓ ａｎｄ Ｓｙｓｔｅｍｓ． Ｅｄｍｏｎｔｏｎ， ２００５３４３－３５０． ［１０］ Ｅｎｒｉｑｕｅ Ｇ， Ｍａｕｒｉｃｉｏ Ａ， Ｐａｕｌａ Ａ， ｅｔ ａｌ． ＢＳＡ Ａ ｃｏｖｅｒａｇｅ ａｌｇｏｒｉｔｈｍ ［Ｃ］ ∥Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ２００３ ＩＥＥＥ／ ＲＳＪ Ｉｎｔ Ｃｏｎｆ ｏｎ Ｉｎｔｅｌｌｉｇｅｎｔ Ｒｏｂｏｔｓ ａｎｄ Ｓｙｓｔｅｍｓ． Ｌａｓ Ｖｅｇａｓ， ２００３１６７９－１６８４． ［１１］ Ｈｏｗｉｅ Ｃ． Ｃｏｖｅｒａｇｅ ｆｏｒ ｒｏｂｏｔｉｃｓＡ ｓｕｒｖｅｙ ｏｆ ｒｅｃｅｎｔ ｒｅｓｕｌｔｓ［Ｊ］． Ａｎｎａｌｓ ｏｆ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ａｒｔｉｆｉｃｉａｌ Ｉｎｔｅｌｌｉｇｅｎｃｅ， ２００１， ３１（１２４）１１３－１２６． ［１２］ Ｙｏａｖ Ｇ， Ｅｌｏｎ Ｒ． Ｓｐａｎｎｉｎｇ⁃ｔｒｅｅ ｂａｓｅｄ ｃｏｖｅｒａｇｅ ｏｆ ｃｏｎｔｉｎｕｏｕｓ ａｒｅａｓ ｂｙ ａ ｍｏｂｉｌｅ ｒｏｂｏｔ［Ｃ］∥Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ＩＥＥＥ Ｃｏｎｆ ｏｎ Ｒｏｂｏｔｉｃｓ ａｎｄ Ａｕｔｏｍａｔｉｏｎ． Ｓｅｏｕｌ， ２００１１９２７－１９３３． ［１３］ Ｄａｎａ Ｒ Ｙ， Ａｌｂｅｒｔ Ｍ Ｂ， Ｂａｒｒｉｅ Ｂ Ｗ， ｅｔ ａｌ． Ｆｉｎｅ⁃ｓｃａｌｅ ｓｅａｆｌｏｏｒ ｓｕｒｖｅｙ ｉｎ ｒｕｇｇｅｄ ｄｅｅｐ⁃ｏｃｅａｎ ｔｅｒｒａｉｎ ｗｉｔｈ ａｎ ａｕｔｏｎｏｍｏｕｓ ｒｏｂｏｔ［Ｃ］∥ ＩＥＥＥ Ｉｎｔ Ｃｏｎｆ ｏｎ Ｒｏｂｏｔｉｃｓ ａｎｄ Ａｕｔｏｍａｔｉｏｎ． Ｓａｎ Ｆｒａｎｃｉｓｃｏ， ２０００ １７８７－１７９２． ［１４］ Ｅｒｃａｎ Ｕ Ａ， Ｈｏｗｉｅ Ｃ， Ａｌｆｒｅｄ Ａ Ｒ， ｅｔ ａｌ． Ｍｏｒｓｅ ｄｅｃｏｍｐｏｓｉｔｉｏｎｓ ｆｏｒ ｃｏｖｅｒａｇｅ ｔａｓｋｓ［Ｊ］． Ｔｈｅ Ｉｎｔ Ｊｏｕｒｎａｌ ｏｆ Ｒｏｂｏｔｉｃｓ Ｒｅｓｅａｒｃｈ， ２００２，２１ （４）３３１－３４４． ［１５］ Ｇｒｏｓｓｂｅｒｇ Ｓ． Ｎｏｎｌｉｎｅａｒ ｎｅｕｒａｌ ｎｅｔｗｏｒｋｓ ｐｒｉｎｃｉｐｌｅｓ， ｍｅｃｈａｎｉｓｍ， ａｎｄ ａｒｃｈｉｔｅｃｔｕｒｅｓ［Ｊ］． Ｎｅｕｒａｌ Ｎｅｔｗｏｒｋｓ， １９８８（１）１７４９－１７６１． ［１６］ Ｈｏｄｇｋｉｎ Ａ Ｌ， Ｈｕｘｌｅｙ Ａ Ｆ． Ａ ｑｕａｎｔｉｔａｔｉｖｅ ｄｅｓｃｒｉｐｔｉｏｎ ｏｆ ｍｅｍｂｒａｎｅ ｃｕｒｒｅｎｔ ａｎｄ ｉｔｓ ａｐｐｌｉｃａｔｉｏｎ ｔｏ ｃｏｎｄｕｃｔｉｏｎ ａｎｄ ｅｘｃｉｔａｔｉｏｎ ｉｎ ｎｅｒｖｅｓ［Ｊ］． Ｐｈｙｓ Ｌｏｎｄ， １９５２，１１７（４）５００－５４４． ［１７］ Ｙａｎｇ Ｓ Ｘ， Ｌｕｏ Ｃ Ｍ． Ａ ｎｅｕｒａｌ ｎｅｔｗｏｒｋ ａｐｐｒｏａｃｈ ｔｏ ｃｏｍｐｌｅｔｅ ｃｏｖｅｒ⁃ ａｇｅ ｐａｔｈ ｐｌａｎｎｉｎｇ ［ Ｊ］． Ｓｙｓｔｅｍｓ， Ｍａｎ ａｎｄ Ｃｙｂｅｒｎｅｔｉｃｓ， Ｐａｒｔ Ｂ ＩＥＥＥ Ｔｒａｎｓａｃｔｉｏｎｓ， ２００４，３４（１）７１８－７２４． ６１矿 冶 工 程第 ３６ 卷 万方数据