浅析拉格朗日有限元数值计算法

1、成绩中国矿业大学2011级硕士研究生课程考试试卷考试科目岩土工程数值计算法考试时间2011.11。27学生姓名夏明学 号 ZS11020068所在院系矿'业工程学院任课教师徐志伟中国矿业大学研究生院培养管理处印制对论文基于拉格朗日差分法的露天边坡稳定性研究中数值计算法浅析基于拉格朗日差分法的露天边坡稳定性研究 一文以易门铜矿露天开采境 界优化方案下的边坡为工程背景，通过岩体构造调查 ,质量分类，室内力学试验， 力学参数的工程处理。建立铜厂露天开采边坡的三维地质模型，采用 DIMINE 数字矿山软件的耦合集成技术四面体网格化将地质模型转换为力学模型， 应用 基于拉格朗日法的有限差分（ F

2、LAC3D ）大变形方法对铜厂露天开采边坡的稳定 性进行了数值模拟， 分析了基于强度折减理论计算出的边坡安全系数以及基于莫 尔库伦屈服准则的边坡开采后的位移、应力等的变化状况 ,得出了易门铜厂露天 矿露天境界优化方案下的边坡的稳定性状况。岩质边坡稳定性评价的方法分主要有 :极限平衡法、数值模拟计算、地质力 学物理模拟试验和其它新方法。 随着计算机技术和计算方法的发展， 复杂的工程 问题可以采用离散化的数值计算方法并借助计算机得到满足工程要求的数值解， 数值模拟技术是现代工程学形成和发展的重要动力之一 .通过计算模拟 ,可以模拟 并得到模拟体内部的应力 -应变关系 ,再现其变形甚至破坏过程及其机

3、制。在岩土 工程数值分析中最常用的数值方法有有限元法、离散元法、边界元法等。拉格朗日差分法（ （FLAC 法）源于流体力学。它首先是 Cundail 在 80 年代提 出 来的，其基本原理类似于离散元法， 但它却能像有限元那样适用于多种材料模式 与边界条件的非规则区域的连续问题求解。在求解的过程中， FLAC 又采用了离 散元的动态松弛法，不需要求解大型联立方程，便于在微机上实现。另一方面， 同以往的差分分析相比， FLAC 在以下几个方面做了较大的改进和发展 :它不但能 处理一般的大变形问题， 而且能模拟岩体沿某一弱面产生的滑移变形。 一般有限 单元法可以用来解决材料非线性问题 ,但对于大变

4、形的几何非线性问题，有限单 元法和边界元法都无能为力 .拉格朗日法是分析非线性大变形问题的数值方法 , 它依然遵循连续介质的假设， 基于拖带坐标系的基本原理。 用差分法或按时步显 式迭代求解，不但可以解决几何非线性，也能解决材料非线性问题。下面是原文的部分描述： 一、露天矿边坡三维地质模型的建立 1。1计算几何模型范围及模拟方案设计a 模拟计算方案设计决定边坡稳定性的因素很多 ,其中包括岩体强度 ,结构特征，水理性质和边 坡几何参数等。根据现场工程地质调查 ,力学试验参数，对露天境界边坡进行修改，修改方案同第五章简化实体模型的建立。b模型范围模型范围需考虑两方面的内容：首先，必须能够覆盖铜厂露

5、天矿边坡岩体的 主要特征，即能容纳前面建立的实体模型；另外，为配合数值计算的需要，要考 虑边坡岩体稳定影响范围以及数值模拟软件以南北向为丫方向，东西为x向, 高程方向为z向。确定的数值计算范围参数见表 71表71模型尺寸参数表X方向尺寸(m)Y方向尺寸52方向尺寸5)490470324c单元划分及计算几何模型单元划分需要考虑三方面的问题。首先，能体现铜厂的地质特征，在岩性，断层和地表高度变化的位置需要细分单元块， 以确保计算模型能够真实的反映岩 体的几何形态；其次，保证数值计算中应力、位移等计算结果的连续性；第三， 由于计算机计算能力和内存大小的限制，单元大小必须能保证计算程序能正常的 运行。

6、根据以上原则，对四面体单元单元块边长以及体积进行限制，四面体单元最小半径一边比设为1.414,落雪组四面体单元块最大体积设为 500,火成胶 结角砾岩组四面体单元块最大体积设为 10000已经建立了边坡地质模型，根据地学实体模型，按照模型范围和单元划分规则，使用DIMINE软件四面体化功能，完成了从地质模型到数值计 算模型的转换。网格的划分由先前地学软件 DIMINE确定，网格的大小为245762格单元，42840个节点阳卜3慮常岐件凹囲体岡对化6的边城權陛圈7-4就1计算惟扯申迪城模型从以上两图中可以看出，通过转换程序的转换，得到了能完全反映地质结构 以及岩性在空间上分布的精确三维力学模型，

7、实现了三维地学模拟与数值模拟的 耦合从图7-4中可以看出，建立的计算模型能较准确地反映地形的微地貌形态以 及能充分体现岩体、断层的尺寸和在空间中的位置 ，露天坑的坡度也在数值计算 模型中也得以显现。1。2计算参数与约束条件a。屈服准则的选择屈服准则是判断材料进入塑性阶段的准则.岩土材料的屈服准则经过几十年 的研究，提出的表达式不下几十种。以莫尔一库伦定律为基础的摩擦屈服准则在 岩石力学与工程的实践中经受了考验，至今仍被广泛的应用。FLAC3D所采用的屈服准则之一为莫尔库伦与拉破坏准则结合的复合准则.b. 计算力学参数根据现场地质调查和岩石力学试验结果可确定岩体力学参数，由于FLAC3D中采用体

8、积模量和剪切模量描述弹性模量和泊松比，所以根据式7-3 5式计算体积模量和剪切模量。(E, v)与(K, G)的转换关系如下E -3(1-2巧 | E |*均2(3)其中：K-体积模量 G剪切模量 E-弹性模量 V-泊松比考虑到计算机计算能力和内存大小的限制，在建立块段模型时将落雪组地层中的 三种岩性(泥质白云岩，灰白色白云岩，青灰色白云岩)合并成一种岩性，这样可 以减少各种岩性间边界上的单元块次分，从而做到减少单元块数目.同时由第四章力学试验可得，落雪组地层的三种岩性力学参数很相似，所以在力学强度上将 其合并也是可行的.计算中所需的力学参数如表72所示。表72各类岩石物理力学参数组名密度Ms

9、/3)体积模童 K (Gpa)剪切模量 G (Gpa)抗拉强度粘聚力 C (Mpa)内摩擦角 申2.80920.19810,2640.9421,297断星纽2,7600.0210.00010.3K3火成胶蜡角砾岩銅17945,6475JI70,069042933.36c。边界约束模型边界约束采用位移(在FLAC3D中实质上是速度约束)约束的方式。底 部所有节点取X、丫、Z三个方向的固定约束；南北边界为 丫方向约束，东西 边界进行X方向约束(即各面的法向方向的约束).坡面为自由面。二、计算结果与分析按前述约束条件，在只考虑重力作用的情况下，进行本构模型为Mohr-Coulomb模型的弹塑性求解，

10、直至系统达到平衡。图77为数值计算过程 中弹塑性求解阶段的系统不平衡力演化全过程曲线.设定体系最大不平衡力与典型内力比值下限为Ixlo 5,迭代计算3683步后，系统达到近似平衡（见图7 一 7）.为便于分析，根据边坡的几何形态变化沿 丫 轴方分别在 丫=2734650m、丫=2734750m、丫=2734850m 和 丫=2735000m 切出 4 个典型剖面。下面就数值分析计算结果从安全系数、变形（位移 ）情况、应力情 况等三个方面来分析边坡岩土体的力学响应特性，以及可能的内在变形破坏机 理。浅析论文中的数值计算方法作者利用FLAC3D数值模拟软件对铜厂露天境界优化方案下的边坡稳定 问题进行了数值模拟分析，并对开挖后的边坡安全系数，位移、应力的变化状况 进行了分析。这为寻找露天矿坑高边坡的稳定性与潜在破坏规律提供了捷径.数值模拟结果可作为生产现场必要的参考.但是，作者未考虑露天边坡的开挖以及台阶局部破坏同时顾及到计算规模和 网格数目的庞大，没有模拟整个露天开采边坡 ，研究结果具有一定的局限性，在 后期的研究中，可考虑更全面可行的模拟技术与方法；对建模耦合技术的理论研究不够深入，对其他方式的耦合方法未进行更深入的讨论和研究，其他耦合方法可能是后期研究的方向；对于四面体网格化还需进一步研究改进， 使得网格更加 符合数值模拟计算网格的需求。在这个学期的学期中，我深刻感受到