岩体失稳破坏过程数值模拟的研究.pdfVIP

岩体失稳破坏过程数值模拟研究。 乔河唐春安徐小穆 东北太学岩百艘裂与失稳研究中心?!甲j‘oc：一6 摘要本文基于岩石介质的宏观非线性主要是出非均质性和各搏异性造成的．应用新的数值计 算软件RFPA，对采动覆岩失稳破坏进行了数值槿拟研究。 关键词 岩体失稳非均质性数值模拟 1引言 近年来，随着计算机运算速度和容量的提高，数值模拟方法在是矿工程领域内应用越来 越广泛。通过数值模拟计算，可以减少实验室变验和现场的工作量，缩小实验和试验范围，同 时耗资少，周期短，还可以发现实验室实验和现场试验难蹦认识的『砭题。 在过去的研究中．由于岩体地质结构和岩体工程结构的复杂性，人们不得不将复杂的同 题抽象成各种简单的问题来进行理论分析(如将岩层简化啦弹眭蜓二i构)．针对各种问题已 经出现各种理论，而且针对各种问题，辽会出现各种新自j毽l℃ 如针对裂纹扩展问题发展了 断裂力学，针对损伤问题发展了损伤力学，针对岩层移动问题发展丁关键层理论等等。就各 种问题发展各种理论固然是正确的，但是连续性和均争睦缓设常常会限制这些理论的发展。 只有通过数值计算，才可能对研充自j对象寻求用解析方法准以得到的解答，以帮助，、蜘对所 研究对象的更深入了解．或对工程实际育更薄人的指寻意辽 如果说．材料的本构关系是第一层；xi二的露池，球幺，电斫裂一々辜，损伤力学，菱键层理论 等。则立是第二层次上的理论。作为数值模拟技术，她!．最：篮呆用第一层次理论。而尽量避免 第二层7炙理论。准确地说，真正的数值模拟应识下依赖于抓些第二量’置理论．咀它应能遥过 第一层次理论模拟出第二屡次理论所描述的现象‘假定进些理论是正确的话J。 RFPA数值模拟软件就是基于这种学术慰想发展起来的岩石破裂过程分析系统。 由东北大学岩石破裂与失稳研究中心基于岩石介质的非均质性和各相异性是导致其变 形非线性的主要固素，用细观线性破坏模拟宏观变形{j线性的RFPA数值模拟软件，能对 岩体从变形、破坏直至运动过程进行研究。本文正是应用这种新的数值计算软件对采动引起 的岩体失稳破坏进行了初步研究， 2基本原理 目前对于岩石变形及其破裂过程的研究，陈了理论分析、相似模拟方法以外，数值模拟 ① 国家杰出青年科学基金资助项目(59,52408 主’法碍弛r越来越普遍的应甩。其中主要包括有限元方法和离散元方法。有限元方法是迄 今为止在应列和变形计算方面最为成熟的7r法之一。岩体既；f：是连续体也不是均匀体。由 丁漫长的地E：年代的作用，岩石内部存在大量微孔隙和微裂纹等内部缺陷，并在一定程度上 受到非连续面的切割。这些都给岩石变形过程的数值模拟分析带来了极大困难。对于具有 明显层面和芊‘理面并且只受低应：灯作旰j的岩体，进!串来得到较快发展的离散单元法得到了 较多的应用。如离散单元法：t：模拟由采矿引起的上覆岩层层间脱落现象等方面取得了一定 研究成果。但遗憾的是，离散元虽然7：分考虑了岩体结构的不连续性，实际上只模拟了岩体 破裂成块以后的运动过程，且这种计算难以预料由于局部破坏而引起的|立力再分布过程，因 而不能模拟由采矿引起的覆岩从经受救荷的变形到破坏的全过程。而且离散元方法更，：的 困难还在于如何将岩体离散成符合买译的块体单元。如牙2说有限元适合于岩体破坏前的变 形分析，那么可以说离散元法则比较适合于研究岩体破坏看的运动分析：， 本文应用的这种新的RFPA软搿就是基于岩百宏观非线性行为可能是具有弹脆性特 征的细观微：元体不断破裂造成的，而遗成微元体不断破裂(并：非某一时列同时破裂)的原因 是微元体材料性质(微元体强差、弹模眚口泊松比等参数)的非均匀性，：黾一个数学上相对简单 但能充分考虑岩石介质复杂性的岩石；陵裂过程分析系统。 3 RFt’A(2D)程序简：介 RFPA(RockFailureProcess Anal,IsisCode)是由东北大学岩石破裂与失稳研究中心开 发的一种基二F有限元计算原理而开发的，但叉不同传统的有限元基本思路，可用于岩石破裂 基理及实际工程失稳破坏全过程研究的大型计算程宁。它的最大特点是仅仅依赖于第一层 次的理论(岩石基本的本构关j暴)。其：蓦本原理： (1)通过考虑岩石的细观非均质髅模拟岩石变形破坏的非线性； (2)通过单元破坏后的材舛性质弱化模拟岩石变形、破裂的非连续性。 该程序的