地下工程的岩石力学数值模拟综述

《地下工程的岩石力学数值模拟综述》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、地下工程的岩石力学数值模拟综述自20世纪70年代以来，数值模拟方法作为一种科学有效而又快速简便的分析方法，被逐步引人到交通、水工、采矿、建筑等地下工程的力学分析、稳定性评价、方案比较等工作中。目前，这种方法已作为一种常用的手段在实际工程中大量采用。1．数值模拟方法选择现今可用于对岩体工程结构进行力学分析的数值方法多种多样，每一种方法有其针对性和特点，对一个具体的问题用数值模拟方法进行分析时，应选择一种最适合该问题的方法进行研究。数值模拟方法的选择，取决于研究对象即岩体工程结构的岩石力学性质和数值模拟的目的。严格地讲，岩石除具有弹性性质外，还具有塑性性质和粘性性质，

2、只不过在特定情况下，某种性质占主导地位而已。在岩体工程实践中，硬岩及应力水平不甚高的中硬岩，其力学性质主要呈现为弹性或弹塑性；高应力环境下的软岩，其力学性质主要呈现为塑性或粘塑性；对于服务时间较长的地下工程，岩石极软或软且应力水平很高，则在计算分析中不能忽视岩石的流变性质。为了达到了解整个岩体工程系统的应力及变形规律的目的，各种数值方法均可采用，但以弹、塑性有限单元法或拉格朗日元法最为适宜。这两种方法的单元划分灵活，计算所需参数较少且易获得，软件也易于得到，成本较低。局部工程结构的应力及变形分析，若岩石中硬以上，则各种方法均可采用；若岩石软弱，则宜采用能进行大变形

3、分析的拉格朗日元法；若岩体可能发生非连续破坏，则宜采用离散单元法、非连续变形法或流形元法或大变形分析的拉格朗日元法；若岩体可能发生非连续破坏，则宜采用离散单元法、非连续变形法或流形元法。要实现了解岩体结构的破坏特征及动态破坏过程的目的，则只能采用离散单元法、非连续变形法或流形元法，因为这些方法就是针对岩体介质的非连续性而提出的。此外，对于一个具体问题，是进行平面分析或是进行三维分析，也需做出恰当选择。严格地讲，所有的问题都是三维的，但如果采用平面分析既能达到目的，计算结果误差也不大，为了降低费用和快速方便起见，则以采用平面分析为宜。反之，则应采用三维程序代码进行计

4、算分析。表1给出了几种数值模拟方法所依据的基本原理、求解方式、离散化方法及其适用条件，可供选择模拟方法时参考。表1几种数值模拟方法的对比数值模拟方法基本原理求解方法离散方式适用条件有限单元法最小势能原理解方程组全区域划分单元岩石中硬以上，小变形，岩体不会发生非连续性破坏如滑动、转动、分离等边界单元法Betti互等定理解方程组边界上划分单元同上离散单元法牛顿运动定律显式差分按结构弱面分布特征划分单元岩石中硬以上，低应力水平，大变形，岩体沿弱面发生非连续性破坏非连续变形法最小势能原理解方程组按主要结构弱面实际情况划分单元大变形，岩体发生非连续性破坏数值流形法最小势能原

5、理解方程组全区域划分单元中硬以上岩体的连续或非连续变形拉各朗日元法牛顿运动定律显式差分全区域划分单元岩石软弱，大变形，岩体的破坏以变形为主2计算模型的建立数值模拟方法确定以后，建立一个能充分符合原型几何特征和物理性质的计算模型至关重要。建立的模型与实际原型越接近，则计算结果越能代表原型的实际情况。但是，实际工程的几何形状及岩体的物理力学性质是极其复杂的，要使建立的计算模型与工程原型在几何和物理性质上完全相同几乎是不可能的，因此，应根据工程的实际情况和研究目的对原型进行适当的工程处理，以便进行计算分析。2．1几何形状若需对较大范围内的地下工程进行力学分析，可按照工程

6、施工设计建立数值模拟分析所需的几何模型。如要对具体部位岩体的应力及变形进行分析，则应按照开挖形成的实际几何形状或统计平均情况建立模型。2．2计算范围计算范围原则上按照开挖(采)活动所引起的围岩应力重新分布的显著影响范围确定。对于隧道型或洞室型开挖工程，理论分析表明，开挖活动的显著影响范围为最小开挖尺寸的3—5倍。对较大范围的地下工程进行力学分析时，其显著影响范围取决于工程的赋存深度、开挖规模与赋存深度之比等。一般情况下，计算范围向上取至地表，向下取开挖空间较小尺寸的O．5～1．0倍即可，水平方向的影响范围自开挖边界向外推最大开挖深度的距离。2．3边界条件在地下工程

7、岩石力学数值模拟计算分析中，计算范围通常取至开挖影响范围之外，除地表为应力边界(一般为零应力边界)外，下边界可为全约束或铅垂方向约束边界。前后左右边界视情况不同而进行相应处理：若仅考虑岩体的自重作用，则前后左右边界可为垂直于边界平面方向上约束或全约束；若有水平构造应力存在，且计算分析所用代码不能对地应力初始化，则前后左右边界应为应力边界。2．4初始状态初始状态包括原始地应力状态、大型地质构造面(岩层层面、岩性界面、断层破碎带、软弱夹层、裂隙等)、地下水赋存及其径流状态等。对一个具体的地下工程而言，其所处区域的原始地应力可以认为其大小和方向是恒定不变的，三个主应力的

8、性质为压应