矿柱回采崩矿方式数值模拟优化及炮孔三维可视化设计
|
摘要
矿柱采场的安全高效回采对提高矿床资源回收率,实现矿山可持续发展具有重要意义,合理的崩矿方式和爆破参数是提高矿柱回采质量的有效途径。本文针对凡口铅锌矿矿柱采场的工程实际,综合运用理论分析、数值模拟、现场探测相结合的方法,采用Surpac、Flac3D及LS-DYNA等软件工具,结合“凡口铅锌矿采场空区三维探测及可视化计算方法研究”科研课题,构建了基于空区实测的矿柱采场三维可视化模型,开展了具有复杂边界矿柱的崩矿方式数值模拟优化研究,优化了矿柱采场回采崩矿方式,实现了矿柱采场回采爆破炮孔三维可视化设计,并对充填体破坏影响最大的边孔进行装药量及装药结构中的不耦合系数进行数值模拟优化,为矿柱安全有效开采提供技术支持。主要研究内容如下：
(1)运用CMS空区探测系统实现了对矿柱采场周边矿房采场的三维探测,构建了具有复杂边界矿柱采场的三维可视化模型,从而为矿柱采场三维炮孔设计及矿柱采场崩矿方式优化提供技术支持；
(2)运用数值分析技术研究了后退式侧向崩矿和水平上向分层崩矿两种崩矿方式位移场、应力场和塑性区的变化规律,优化了矿柱回采的崩矿方式；
(3)运用三维可视化软件Surpac创建了矿柱采场的中深孔爆破设计数据库,实现了矿柱采场回采爆破炮孔的三维可视化设计；(4)运用显示动力分析软件LS-DYNA对炮孔装药量和装药结构进行了数值模拟优化,获得给定边孔距条件下的边孔单段最大安全药量及合理的装药不耦合系数。
本文紧密结合工程实践,运用数值计算与可视化建模技术对矿柱采场崩矿方式,矿柱回采爆破炮孔设计以及边孔装药量与装药结构中的不耦合系数等相关技术进行了较为深入的研究,为深部矿柱采场的高效回采奠定了基础,具有较为重要的理论意义和工程应用价值。
Safe and efficient pillar mining is important for improving mine resources recovery and achieving sustainable development. Proper pillar caving method and blasting parameters are two effective ways to improve the quality of pillar mining. According to engineering of Fankou lead-zinc mine, the methods including theoretical analysis, numerical simulation and on-site investigation were used comprehensively in this paper. Combined with "3D cavity monitor and the calculation method of visualization in FanKou lead-zinc mine" research project, and by adopting software tools, such as Surpac, Flac3D, LS-DYNA and etc.,3D visual pillar model based on cavity monitor has been built, and the research has been carried out about numerical simulation optimization of mining method of pillar with complex boundary. In the end,3D blast hole design of pillar has been achieved and the charge weight and decouple coefficients of periphery-hole which had the greatest impact on filling were optimized, which provided technical support for safe and effective pillar mining. The main research contents are as follows:
(1) 3D visual model with complex boundary is built by using CMS to detect the stope beside pillar, which provides technical support for 3D blast hole design of pillar and pillar mining method.
(2) Displacement field, stress field and the variation of plastic zone of the retreating side pillar mining method and the level caving of pillar have been analyzed by using numerical simulation software, thus the pillar mining method was optimized.
(3) Deep hole blasting database is created by 3D design software and then 3D visual blast hole design of pillar mining has been achieved.
(4) Using Explicit dynamic analysis software LS-DYNA to optimize hole charge weight and charge structure obtained the maximum safe dose and appropriate decouple coefficients.
This paper closely combining with the engineering practice, studied pillar mining method,3D blast hole design, the maximum safe dose, decouple coefficients by using numerical simulation and visual modeling, which provides a basis for safe and efficient mining of deep pillar, and it has important theoretical and engineering application value.
(1)运用CMS空区探测系统实现了对矿柱采场周边矿房采场的三维探测,构建了具有复杂边界矿柱采场的三维可视化模型,从而为矿柱采场三维炮孔设计及矿柱采场崩矿方式优化提供技术支持；
(2)运用数值分析技术研究了后退式侧向崩矿和水平上向分层崩矿两种崩矿方式位移场、应力场和塑性区的变化规律,优化了矿柱回采的崩矿方式；
(3)运用三维可视化软件Surpac创建了矿柱采场的中深孔爆破设计数据库,实现了矿柱采场回采爆破炮孔的三维可视化设计；(4)运用显示动力分析软件LS-DYNA对炮孔装药量和装药结构进行了数值模拟优化,获得给定边孔距条件下的边孔单段最大安全药量及合理的装药不耦合系数。
本文紧密结合工程实践,运用数值计算与可视化建模技术对矿柱采场崩矿方式,矿柱回采爆破炮孔设计以及边孔装药量与装药结构中的不耦合系数等相关技术进行了较为深入的研究,为深部矿柱采场的高效回采奠定了基础,具有较为重要的理论意义和工程应用价值。
Safe and efficient pillar mining is important for improving mine resources recovery and achieving sustainable development. Proper pillar caving method and blasting parameters are two effective ways to improve the quality of pillar mining. According to engineering of Fankou lead-zinc mine, the methods including theoretical analysis, numerical simulation and on-site investigation were used comprehensively in this paper. Combined with "3D cavity monitor and the calculation method of visualization in FanKou lead-zinc mine" research project, and by adopting software tools, such as Surpac, Flac3D, LS-DYNA and etc.,3D visual pillar model based on cavity monitor has been built, and the research has been carried out about numerical simulation optimization of mining method of pillar with complex boundary. In the end,3D blast hole design of pillar has been achieved and the charge weight and decouple coefficients of periphery-hole which had the greatest impact on filling were optimized, which provided technical support for safe and effective pillar mining. The main research contents are as follows:
(1) 3D visual model with complex boundary is built by using CMS to detect the stope beside pillar, which provides technical support for 3D blast hole design of pillar and pillar mining method.
(2) Displacement field, stress field and the variation of plastic zone of the retreating side pillar mining method and the level caving of pillar have been analyzed by using numerical simulation software, thus the pillar mining method was optimized.
(3) Deep hole blasting database is created by 3D design software and then 3D visual blast hole design of pillar mining has been achieved.
(4) Using Explicit dynamic analysis software LS-DYNA to optimize hole charge weight and charge structure obtained the maximum safe dose and appropriate decouple coefficients.
This paper closely combining with the engineering practice, studied pillar mining method,3D blast hole design, the maximum safe dose, decouple coefficients by using numerical simulation and visual modeling, which provides a basis for safe and efficient mining of deep pillar, and it has important theoretical and engineering application value.
引文
[1]杨彪.基于CMS实测的采空区危险度分析及其处理[硕士].中南大学,2008.
[2]鹿浩.基于空区实测的采场回采质量评价与矿柱安全高效回采研究[硕士].中南大学,2009.
[3]陈昌民.康家湾铅锌矿残留富矿回采综合技术研究[硕士].中南大学,2002.
[4]秦杨.张家滩磷矿保安矿柱开采技术研究[硕士].江西理工大学,2010.
[5]郑学敏.狮子山铜矿西山矿段特大空区下矿柱回采的稳定性研究[硕士].中南大学,2002.
[6]万海涛.画眉坳钨矿矿柱回采及采空区处理数值模拟研究[硕士].江西理工大学,2010.
[7]吴亚斌.基于CMS实测的采空区群稳定性数值模拟研究[硕士].中南大学,2007.
[8]Luozhouquan,liuxiaoming,Deposit 3D modeling and application.中南大学学报(英文版).2007,(3).
[9]寇向宇,贾明涛,王李管,等.基于CMS及DIMINE-FLAC-(3D)耦合技术的采空区稳定性分析与评价[J].矿业工程研究,2010,(01).
[10]陈广仁.利用CMS提高采场回采质量评价水平[J].科技导报,2010,(16)：8.
[11]张木毅.凡口铅锌矿盘区上向中深孔配套机械化采矿新工艺研究[硕士].中南大学,2002.
[12]许毓海,许新启,桑守勤.凡口铅锌矿盘区机械化非胶结充填采矿试验研究[J].矿冶,2000,(1)：1-4.
[13]刘晓明,罗周全,张保,等.CMS空区探测的环境影响及误差分析与纠正[J].矿冶工程.2008,(04).
[14]张高青.大运公路介休至霍州段灵石隧道采空区处理技术研究[硕士].长安大学,2010.
[15]解海军.煤矿积水采空区瞬变电磁法探测技术研究[博士].中国地质大学(北京),2009.
[16]丁海鹏.本溪市东山煤矿采空区的瑞雷波勘察方法应用[硕士].东北大学,2008.
[17]靳辉.地震横波反射方法在煤矿采空区探测中的应用研究[硕士].东北大学,2008.
[18]宗志刚.地震勘探方法在探测煤矿采空区中的应用研究[硕士].中国地质大 学(北京),2006.
[19]钟韬.超高密度电法在探测采空区中的应用研究[硕士].成都理工大学,2008.
[20]张志伟.大深度三维高密度电阻率法接收装置的研制与实验[硕士].吉林大学,2010.
[21]刘海生.高密度电法在探测煤矿地下采空区中的应用研究[硕士].太原理工大学,2006.
[22]刘晓明.基于实测的采空区三维建模及其衍生技术的研究与应用[硕士].中南大学,2007.
[23]罗周全,杨彪,刘晓明,等.基于CMS实测及Midas-FLAC-(3D)耦合的复杂空区群稳定性分析[J].矿冶工程,2010,(6)：1～5.
[24]彭欣.复杂采空区稳定性及近区开采安全性研究[博士].中南大学,2008.
[25]李一帆.数字矿山信息系统的研究及应用[博士].中国科学院研究生院(武汉岩土力学研究所),2007.
[26]张力岩.数字矿山中三维地质模拟与体视化研究[硕士].中国科学院研究生院(遥感应用研究所),2004.
[27]庄玲.数字矿山可视化技术的研究与实现[硕士].山东科技大学,2003.
[28]吴斐.数字矿山中三维空间数据模型及其应用研究[硕士].中国科学院研究生院(遥感应用研究所),2006.
[29]杨晓坤,秦德先,冯美丽,等.基于Surpac的矿山三维地学模型及综合信息成矿预测研究[J].地质与资源,2008,(01).
[30]侯恩科,吴立新.面向地质建模的三维体元拓扑数据模型研究[J].武汉大学学报(信息科学版),2002,(05).
[31]池秀文,赵旭,黄解军.基于GIS的三维矿山爆破信息系统开发与应用[J].金属矿山,2010,(1)：138～141.
[32]冯超东,杨鹏,基于SURPAC的复杂矿床三维建模技术.In中国新疆乌鲁木齐,2007.
[33]冯红刚,王李管,毕林.基于DIMINE软件的三维可视化技术及其工程应用[J].矿业工程研究,2010,(3)：7～10.
[34]房智恒,王李管,何远富.基于DIMINE软件的采矿方法真三维设计研究与实现[J].金属矿山,2009,(5)：129～131.
[35]肖英才,王李管,朱明海,等.基于DIMINE软件的爆破设计[J].现代矿业,2011,(1)：26～29.
[36]郭子强,郑新厂.锚固巷道稳定性数值分析[J].中州煤炭,2008,(6)：12～13.
[37]姜谙男,赵德孝,鲁炳强,等.金山店铁矿采区围岩及回填体稳定性数值模 拟及分析[J].有色矿山,2002,(6)：6～8.
[38]李必红,双娜,刘庆成.有限差分算法求解砂岩型铀矿体氡气场分布初探[J].东华理工学院学报,2006,(4)：311～313.
[39]Eissa E A, Rofail N B, Ashmawy L S, et al. Neutron activation analysis of an Egyptian monazite ore sample[J]. Czechoslovak Journal of Physics,1999, (1).
[40]Liu J, Wu A, Yang B, et al. Dynamic experiment and numerical simulation of solute transmission in heap leaching processing[J]. Journal of Central South University of Technology,2007, (6).
[41]沈慧明,许振华,朱利平,等.残矿回采顺序优化与复杂采空区稳定性的有限元模拟研究[J].中国矿业,2011,(1)：78～81.
[42]张世雄,任高峰,彭涛,等.矿柱回采设计的三维有限元数值模拟研究[J].西部探矿工程,2004,(7)：79～81.
[43]许江,彭守建,陶云奇,等.天府三汇一矿地应力场三维有限元分析及其分布特征[J].重庆大学学报,2010,(5)：105～113.
[44]王清来,许振华,朱利平,等.复杂采空区条件下残矿回收与采区稳定性的有限元数值模拟研究[J].金属矿山,2010,(7)：37～40.
[45]陈赞成,侯克鹏,杨八九.某矿采空区稳定性三维有限元数值模拟[J].有色金属,2010,(3)：142～145.
[46]曹平,潘长良,李立明,薄矿脉开采条件下竖井变形三维不连续位移边界元法分析.In 中国北京,1994.
[47]徐长佑,郑泽岱,正交异性岩体中地下开挖的边界元计算模型.In中国上海,1990.
[48]庞建勇.不均匀介质样条边界元数值计算方法及应用[J].安徽理工大学学报(自然科学版),2004,(1)：25～28.
[49]王泳嘉.边界元法在岩石力学中的应用[J].岩石力学与工程学报,1986,(2)：205-222.
[50]陈合龙.有效空间中的损伤力学边界元理论与应用[硕士].浙江大学,2008.
[51]陈庆发,周科平,胡建华.高峰矿105号矿体碎裂矿段采空区稳定性离散元分析[J].矿冶工程,2009,(4)：14～17.
[52]方新秋,万德钧,王庆,等.离散元法在分析综放采场矿压中的应用[J].湘潭矿业学院学报,2003,(4)：11～14.
[53]徐国元,赵建平.基于刚体离散元模型的矿岩散体流动规律[J].辽宁工程技术大学学报(自然科学版),2009,(1)：24～27.
[54]李作良,李一帆.某钨矿采空区稳定性的离散元数值模拟[J].中国钨业, 2008,(6)：1～4.
[55]王泳嘉,刘国兴,邢纪波.离散元法在崩落法放矿中应用的研究[J].有色金属,1987,(2)：20-26.
[56]张飞,周连春,王琰,等FLAC3D数值模拟技术在瓦斯抽放中的应用[J].煤矿机电,2010,(6)：32～34.
[57]曾壮丽,江学良,刘自由.地下洞室顶板安全厚度的Flac～(3D)分析[J].四川建材,2010,(3)：51～53.
[58]于振辉,刘增东,郑昌晶,等.基坑支护工程稳定性的FLAC-3D数值模拟分析[J].西部探矿工程,2010,(11)：32-34.
[59]田树昆,曹兰柱.基于CAD与ANSYS的FLAC～(3D)边坡模拟分析[J].微计算机信息,2010,(33)：259-260.
[60]叶海旺,常剑,周磊.基于FLAC～(3D)的采场结构参数优化[J].金属矿山,2010,(12)：6～8.
[61]王新民,王长军,张钦礼,等.基于ANSYS程序下的采场稳定性分析[J].金属矿山,2008,(8)：17-20.
[62]王忠昶,张文泉.软岩矿层条带开采合理采留宽度的ANSYS三维流变模拟[J].矿山压力与顶板管理,2003,(2)：110～112.
[63]胡建华,周科平,古德生,等.基于RFPA～(2D)的顶板诱导崩落时变效应数值模拟[J].中国矿业,2007,(10)：86～88.
[64]李连崇,唐春安,邢军,等.节理岩质边坡变形破坏的RFPA模拟分析[J].东北大学学报,2006,(5)：559～562.
[65]贾兰.基于RFPA～(2D)的宝日希勒露天矿采场边坡稳定性分析[硕士].辽宁工程技术大学,2008.
[66]Cavity Monitoring System User Manual(Version 2.3)[M].Canada:Optech System Corporation,1996.
[67]Surpac User Manual(Version 4.0).Australia: Surpac SoftwareInternationa Corporation,2000.
[68]邓念东,侯恩科,霍正东,等.应用SURPAC实现矿床三维可视化地质模型动态更新[J].陕西煤炭,2006,(3)：8-9.
[69]王斌,刘保顺,王涛,等.基于Surpac的矿山三维可视化地质模型的研究与应用[J].中国矿业,2011,(2)：106～109.
[70]肖兴平,阮俊,于小芳,等.基于Surpac的地下水污染场地调查三维模型的建立[J].地理空间信息,2008,(6)：44-46.
[71]杨晓坤,秦德先,冯美丽,等.基于Surpac的矿山三维地学模型及综合信息 成矿预测研究[J].地质与资源,2008,(1)：61～64.
[72]Itasca Consulting Group Inc. FLAC3D(Fast Lagrangian Analysis of Continua in 3 Dimensions),Version 2.10,Users manual.USA:Itasca Consulting Group Inc,2002
[73]Carranza-Torres C, Fairhurst C. The elasto-plastic response of underground excavations in rock masses that satisfy the Hoek-Brown failure criterion[J]. International Journal of Rock Mechanics and Mining Sciences,1999,36(6):777-809.
[74]Haibo L, Xiang X, Jianchun L, et al. Rock damage control in bedrock blasting excavation for a nuclear power plant[J]. International Journal of Rock Mechanics and Mining Sciences,2011,48(2):210-218.
[75]SUN X, CAI F, YANG J, et al. Numerical simulation of the effect of coupling support of bolt-mesh-anchor in deep tunnel[J]. Mining Science and Technology (China),2009,19(3):352-357.
[76]陈新泽,唐辉明,杨有成,等.基于FLAC～(3D)强度折减法滑坡三维稳定性研究——以三峡库区白果树古滑坡群为例[J].水文地质工程地质,2008,(2)：24～29.
[77]房晓敏,孟祥瑞,徐铖辉.3DGIS构模与FLAC-3D建模网格数据融合技术研究[J].煤炭技术,2010,(4)：161～163.
[78]寇向宇,贾明涛,王李管,等.用于FLAC～(3D)力学分析的矿床块段模型优化[J].金属矿山,2011,(1)：15～18.
[79]苏国韶,符兴义,李书东.基于FLAC～(3D)的三维地应力场反演分析[J].人民黄河,2011,(2)：142～145.
[80]薛雷,李维朝,孙强,等.基于数字图像像素单元建立准三维FLAC～(3D)模型[J].岩土力学,2010,(6)：2001～2005.
[81]刘科伟,李夕兵,宫凤强,等.基于CALS及Surpac-FLAC～3D耦合技术的复杂空区稳定性分析[J].岩石力学与工程学报.2008(09).
[82]林杭,曹平,李江腾,等.基于SURPAC的FLAC-(3D)三维模型自动构建[J].中国矿业大学学报.2008(03).
[83]罗周全,吴亚斌,刘晓明,等.基于SURPAC的复杂地质体FLAC-(3D)模型生成技术[J].岩土力学.2008(05).
[84]Rimantas Barauskasa, Ausra Abraitiene. Computational analysis of impact of a bullet against the multilayer fabrics in LS-DYNA[J]. International Journal of Impact Engineering 34 (2007) 1286-1305
[85]Wilhelm Rust, Karl Schweizerhof. Finite element limit load analysis of thin-walled structures by ANSYS (implicit), LS-DYNA (explicit) and in combination[J]. Thin-Walled Structures 41 (2003) 227-244
[86]Norman F. Knight Jr, Navin Jaunky, Robin E. Lawson, et al. Penetration simulation for uncontained engine debris impact on fuselage-like panels using LS-DYNA[J]. Finite Elements in Analysis and Design 36 (2000) 99-133
[87]丁黄平.节理裂隙岩体隧道爆破成型效果研究[博士].吉林大学,2009.
[88]张冬茵，金星,丁海平,等.(?)NSYS/LS-DYNA在地震工程中的应用[J].地震工程与工程振动,2005,(04).
[89]贾宪振,胡毅亭,董明荣,等.基于ANSYS/LS-DYNA模拟水下爆炸冲击波的等效质量法[J].弹箭与制导学报,2008,(03).
[90]崔伟峰,杜耀志,周腊久,等.Ls-dyna在聚能射流模拟中的应用[J].采矿技术,2007,(03).
[2]鹿浩.基于空区实测的采场回采质量评价与矿柱安全高效回采研究[硕士].中南大学,2009.
[3]陈昌民.康家湾铅锌矿残留富矿回采综合技术研究[硕士].中南大学,2002.
[4]秦杨.张家滩磷矿保安矿柱开采技术研究[硕士].江西理工大学,2010.
[5]郑学敏.狮子山铜矿西山矿段特大空区下矿柱回采的稳定性研究[硕士].中南大学,2002.
[6]万海涛.画眉坳钨矿矿柱回采及采空区处理数值模拟研究[硕士].江西理工大学,2010.
[7]吴亚斌.基于CMS实测的采空区群稳定性数值模拟研究[硕士].中南大学,2007.
[8]Luozhouquan,liuxiaoming,Deposit 3D modeling and application.中南大学学报(英文版).2007,(3).
[9]寇向宇,贾明涛,王李管,等.基于CMS及DIMINE-FLAC-(3D)耦合技术的采空区稳定性分析与评价[J].矿业工程研究,2010,(01).
[10]陈广仁.利用CMS提高采场回采质量评价水平[J].科技导报,2010,(16)：8.
[11]张木毅.凡口铅锌矿盘区上向中深孔配套机械化采矿新工艺研究[硕士].中南大学,2002.
[12]许毓海,许新启,桑守勤.凡口铅锌矿盘区机械化非胶结充填采矿试验研究[J].矿冶,2000,(1)：1-4.
[13]刘晓明,罗周全,张保,等.CMS空区探测的环境影响及误差分析与纠正[J].矿冶工程.2008,(04).
[14]张高青.大运公路介休至霍州段灵石隧道采空区处理技术研究[硕士].长安大学,2010.
[15]解海军.煤矿积水采空区瞬变电磁法探测技术研究[博士].中国地质大学(北京),2009.
[16]丁海鹏.本溪市东山煤矿采空区的瑞雷波勘察方法应用[硕士].东北大学,2008.
[17]靳辉.地震横波反射方法在煤矿采空区探测中的应用研究[硕士].东北大学,2008.
[18]宗志刚.地震勘探方法在探测煤矿采空区中的应用研究[硕士].中国地质大 学(北京),2006.
[19]钟韬.超高密度电法在探测采空区中的应用研究[硕士].成都理工大学,2008.
[20]张志伟.大深度三维高密度电阻率法接收装置的研制与实验[硕士].吉林大学,2010.
[21]刘海生.高密度电法在探测煤矿地下采空区中的应用研究[硕士].太原理工大学,2006.
[22]刘晓明.基于实测的采空区三维建模及其衍生技术的研究与应用[硕士].中南大学,2007.
[23]罗周全,杨彪,刘晓明,等.基于CMS实测及Midas-FLAC-(3D)耦合的复杂空区群稳定性分析[J].矿冶工程,2010,(6)：1～5.
[24]彭欣.复杂采空区稳定性及近区开采安全性研究[博士].中南大学,2008.
[25]李一帆.数字矿山信息系统的研究及应用[博士].中国科学院研究生院(武汉岩土力学研究所),2007.
[26]张力岩.数字矿山中三维地质模拟与体视化研究[硕士].中国科学院研究生院(遥感应用研究所),2004.
[27]庄玲.数字矿山可视化技术的研究与实现[硕士].山东科技大学,2003.
[28]吴斐.数字矿山中三维空间数据模型及其应用研究[硕士].中国科学院研究生院(遥感应用研究所),2006.
[29]杨晓坤,秦德先,冯美丽,等.基于Surpac的矿山三维地学模型及综合信息成矿预测研究[J].地质与资源,2008,(01).
[30]侯恩科,吴立新.面向地质建模的三维体元拓扑数据模型研究[J].武汉大学学报(信息科学版),2002,(05).
[31]池秀文,赵旭,黄解军.基于GIS的三维矿山爆破信息系统开发与应用[J].金属矿山,2010,(1)：138～141.
[32]冯超东,杨鹏,基于SURPAC的复杂矿床三维建模技术.In中国新疆乌鲁木齐,2007.
[33]冯红刚,王李管,毕林.基于DIMINE软件的三维可视化技术及其工程应用[J].矿业工程研究,2010,(3)：7～10.
[34]房智恒,王李管,何远富.基于DIMINE软件的采矿方法真三维设计研究与实现[J].金属矿山,2009,(5)：129～131.
[35]肖英才,王李管,朱明海,等.基于DIMINE软件的爆破设计[J].现代矿业,2011,(1)：26～29.
[36]郭子强,郑新厂.锚固巷道稳定性数值分析[J].中州煤炭,2008,(6)：12～13.
[37]姜谙男,赵德孝,鲁炳强,等.金山店铁矿采区围岩及回填体稳定性数值模 拟及分析[J].有色矿山,2002,(6)：6～8.
[38]李必红,双娜,刘庆成.有限差分算法求解砂岩型铀矿体氡气场分布初探[J].东华理工学院学报,2006,(4)：311～313.
[39]Eissa E A, Rofail N B, Ashmawy L S, et al. Neutron activation analysis of an Egyptian monazite ore sample[J]. Czechoslovak Journal of Physics,1999, (1).
[40]Liu J, Wu A, Yang B, et al. Dynamic experiment and numerical simulation of solute transmission in heap leaching processing[J]. Journal of Central South University of Technology,2007, (6).
[41]沈慧明,许振华,朱利平,等.残矿回采顺序优化与复杂采空区稳定性的有限元模拟研究[J].中国矿业,2011,(1)：78～81.
[42]张世雄,任高峰,彭涛,等.矿柱回采设计的三维有限元数值模拟研究[J].西部探矿工程,2004,(7)：79～81.
[43]许江,彭守建,陶云奇,等.天府三汇一矿地应力场三维有限元分析及其分布特征[J].重庆大学学报,2010,(5)：105～113.
[44]王清来,许振华,朱利平,等.复杂采空区条件下残矿回收与采区稳定性的有限元数值模拟研究[J].金属矿山,2010,(7)：37～40.
[45]陈赞成,侯克鹏,杨八九.某矿采空区稳定性三维有限元数值模拟[J].有色金属,2010,(3)：142～145.
[46]曹平,潘长良,李立明,薄矿脉开采条件下竖井变形三维不连续位移边界元法分析.In 中国北京,1994.
[47]徐长佑,郑泽岱,正交异性岩体中地下开挖的边界元计算模型.In中国上海,1990.
[48]庞建勇.不均匀介质样条边界元数值计算方法及应用[J].安徽理工大学学报(自然科学版),2004,(1)：25～28.
[49]王泳嘉.边界元法在岩石力学中的应用[J].岩石力学与工程学报,1986,(2)：205-222.
[50]陈合龙.有效空间中的损伤力学边界元理论与应用[硕士].浙江大学,2008.
[51]陈庆发,周科平,胡建华.高峰矿105号矿体碎裂矿段采空区稳定性离散元分析[J].矿冶工程,2009,(4)：14～17.
[52]方新秋,万德钧,王庆,等.离散元法在分析综放采场矿压中的应用[J].湘潭矿业学院学报,2003,(4)：11～14.
[53]徐国元,赵建平.基于刚体离散元模型的矿岩散体流动规律[J].辽宁工程技术大学学报(自然科学版),2009,(1)：24～27.
[54]李作良,李一帆.某钨矿采空区稳定性的离散元数值模拟[J].中国钨业, 2008,(6)：1～4.
[55]王泳嘉,刘国兴,邢纪波.离散元法在崩落法放矿中应用的研究[J].有色金属,1987,(2)：20-26.
[56]张飞,周连春,王琰,等FLAC3D数值模拟技术在瓦斯抽放中的应用[J].煤矿机电,2010,(6)：32～34.
[57]曾壮丽,江学良,刘自由.地下洞室顶板安全厚度的Flac～(3D)分析[J].四川建材,2010,(3)：51～53.
[58]于振辉,刘增东,郑昌晶,等.基坑支护工程稳定性的FLAC-3D数值模拟分析[J].西部探矿工程,2010,(11)：32-34.
[59]田树昆,曹兰柱.基于CAD与ANSYS的FLAC～(3D)边坡模拟分析[J].微计算机信息,2010,(33)：259-260.
[60]叶海旺,常剑,周磊.基于FLAC～(3D)的采场结构参数优化[J].金属矿山,2010,(12)：6～8.
[61]王新民,王长军,张钦礼,等.基于ANSYS程序下的采场稳定性分析[J].金属矿山,2008,(8)：17-20.
[62]王忠昶,张文泉.软岩矿层条带开采合理采留宽度的ANSYS三维流变模拟[J].矿山压力与顶板管理,2003,(2)：110～112.
[63]胡建华,周科平,古德生,等.基于RFPA～(2D)的顶板诱导崩落时变效应数值模拟[J].中国矿业,2007,(10)：86～88.
[64]李连崇,唐春安,邢军,等.节理岩质边坡变形破坏的RFPA模拟分析[J].东北大学学报,2006,(5)：559～562.
[65]贾兰.基于RFPA～(2D)的宝日希勒露天矿采场边坡稳定性分析[硕士].辽宁工程技术大学,2008.
[66]Cavity Monitoring System User Manual(Version 2.3)[M].Canada:Optech System Corporation,1996.
[67]Surpac User Manual(Version 4.0).Australia: Surpac SoftwareInternationa Corporation,2000.
[68]邓念东,侯恩科,霍正东,等.应用SURPAC实现矿床三维可视化地质模型动态更新[J].陕西煤炭,2006,(3)：8-9.
[69]王斌,刘保顺,王涛,等.基于Surpac的矿山三维可视化地质模型的研究与应用[J].中国矿业,2011,(2)：106～109.
[70]肖兴平,阮俊,于小芳,等.基于Surpac的地下水污染场地调查三维模型的建立[J].地理空间信息,2008,(6)：44-46.
[71]杨晓坤,秦德先,冯美丽,等.基于Surpac的矿山三维地学模型及综合信息 成矿预测研究[J].地质与资源,2008,(1)：61～64.
[72]Itasca Consulting Group Inc. FLAC3D(Fast Lagrangian Analysis of Continua in 3 Dimensions),Version 2.10,Users manual.USA:Itasca Consulting Group Inc,2002
[73]Carranza-Torres C, Fairhurst C. The elasto-plastic response of underground excavations in rock masses that satisfy the Hoek-Brown failure criterion[J]. International Journal of Rock Mechanics and Mining Sciences,1999,36(6):777-809.
[74]Haibo L, Xiang X, Jianchun L, et al. Rock damage control in bedrock blasting excavation for a nuclear power plant[J]. International Journal of Rock Mechanics and Mining Sciences,2011,48(2):210-218.
[75]SUN X, CAI F, YANG J, et al. Numerical simulation of the effect of coupling support of bolt-mesh-anchor in deep tunnel[J]. Mining Science and Technology (China),2009,19(3):352-357.
[76]陈新泽,唐辉明,杨有成,等.基于FLAC～(3D)强度折减法滑坡三维稳定性研究——以三峡库区白果树古滑坡群为例[J].水文地质工程地质,2008,(2)：24～29.
[77]房晓敏,孟祥瑞,徐铖辉.3DGIS构模与FLAC-3D建模网格数据融合技术研究[J].煤炭技术,2010,(4)：161～163.
[78]寇向宇,贾明涛,王李管,等.用于FLAC～(3D)力学分析的矿床块段模型优化[J].金属矿山,2011,(1)：15～18.
[79]苏国韶,符兴义,李书东.基于FLAC～(3D)的三维地应力场反演分析[J].人民黄河,2011,(2)：142～145.
[80]薛雷,李维朝,孙强,等.基于数字图像像素单元建立准三维FLAC～(3D)模型[J].岩土力学,2010,(6)：2001～2005.
[81]刘科伟,李夕兵,宫凤强,等.基于CALS及Surpac-FLAC～3D耦合技术的复杂空区稳定性分析[J].岩石力学与工程学报.2008(09).
[82]林杭,曹平,李江腾,等.基于SURPAC的FLAC-(3D)三维模型自动构建[J].中国矿业大学学报.2008(03).
[83]罗周全,吴亚斌,刘晓明,等.基于SURPAC的复杂地质体FLAC-(3D)模型生成技术[J].岩土力学.2008(05).
[84]Rimantas Barauskasa, Ausra Abraitiene. Computational analysis of impact of a bullet against the multilayer fabrics in LS-DYNA[J]. International Journal of Impact Engineering 34 (2007) 1286-1305
[85]Wilhelm Rust, Karl Schweizerhof. Finite element limit load analysis of thin-walled structures by ANSYS (implicit), LS-DYNA (explicit) and in combination[J]. Thin-Walled Structures 41 (2003) 227-244
[86]Norman F. Knight Jr, Navin Jaunky, Robin E. Lawson, et al. Penetration simulation for uncontained engine debris impact on fuselage-like panels using LS-DYNA[J]. Finite Elements in Analysis and Design 36 (2000) 99-133
[87]丁黄平.节理裂隙岩体隧道爆破成型效果研究[博士].吉林大学,2009.
[88]张冬茵，金星,丁海平,等.(?)NSYS/LS-DYNA在地震工程中的应用[J].地震工程与工程振动,2005,(04).
[89]贾宪振,胡毅亭,董明荣,等.基于ANSYS/LS-DYNA模拟水下爆炸冲击波的等效质量法[J].弹箭与制导学报,2008,(03).
[90]崔伟峰,杜耀志,周腊久,等.Ls-dyna在聚能射流模拟中的应用[J].采矿技术,2007,(03).