文档介绍：第1０章 岩石力学的数值模拟
随着计算机软硬件技术的迅速发展，使岩石力学有了长足的进步，特别在岩石力学的数值计算和模拟方面发展尤为迅速，“力学——迎接2１世纪新的挑战"的一封信中对力学发展趋势总结的那样“今日力学是一门用计算机计算去回答一切宏观的实际科学技术问题，计算方法非常重要".岩石力学和其他力学学科一样，需要数值计算方法并推动岩石力学的发展.
岩石介质不同于金属材料,在数值计算方面具有其独特的特点［２0５]：
（１）岩石介质是赋存于地壳中的各向异性天然介质。
（2)岩石介质被众多的节理、裂缝等弱面所切割而呈现高度的非均质性，而其物理、化学及力学性质具有随机性特点。
(３）岩石介质赋存时以受压为主，而且抗压强度远大于抗拉强度。
(4）岩石力学与工程问题在时空分布上较广，从本质上讲都是三维问题。
(5）岩石工程一般无法进行原型试验，而实验室测得的数据不能直接应用于工程设计和计算。
(6)岩石力学与工程具有数据有限问题。
数值计算方法经过几十年的发展,目前已形成许多种岩石力学计算方法，主要有有限元法、边界元法、有限差分法、离散元法、流形元法、拉格朗日元法、不连续变形法及无单元法等。它们各有优缺点,有限元的理论基础和应用比较成熟，在金属材料和构件的计算中应用十分成功,但它是以连续介质为基础，似乎与岩体的非连续性有一定差距，流形元等数值方法虽然考虑了岩体中节理效应，但其理论基础还不完全成熟。相信在不久的将来，肯定会出现完全适合于岩体材料和工程的数值计算方法
[２0６～20８]。
１0。１ 岩石力学的有限元分析［２０9~２13］
有限元法（ｆinite element metｈｏd,FＥM），有限元已成功地求解了许多复杂的岩石力学与工程问题。被广大岩石力学研究与工程技术人员喻为解决岩石工程问题的有效工具。有限元法是根据变分原理求解数学物理方程的一种数值方法。有限元法把连续体离散成有限个单元，每个单元的场函数只包含有限个节点参量的简单场函数，，求解这些离散方程组,就是有限元法的精髓所在。虽然求解时把连续函数转化为求解有限个离散点处的函数值，但只要单元划分得充分小时，足可以满足计算要求。
有限元法求解问题时一般遵循以下步骤:
(1)有限元计算模型的建立，包括模型单元的划分、确定边界条件。
（2）对单元体进行力学分析,包括求解节点位移、单元应变和单元应力。
(３）对计算模型进行分析。
(4）进行计算分析。
线弹性有限元法的基本方程
线弹性有限元是弹塑性有限元、损伤有限元、、均质、小变形和完全弹性。
有限元法求解弹性力学问题时通常以位移作为基本未知量,单元位移是以单元节点位移为基本未知量，选择合理的位移插值函数，将单元位移表达为节点坐标的连续函数,。
图10－１为三种最常见单元形式，即三角形、四边形及四面体单元。它们的形函数分别为:
图10-1 有限元的三种基本单元形式
（a）三角形单元（b）四边形单元（ｃ）四面体单元
三角形的形函数
式中，S为三角形面积；；;。
四边形的形函数
式中,位移量为；；；。
四面体的形函数
式中，Ｖ为四面体的体积。
单元在直角坐标轴中位移分量分别为ｕ,v,，因此单元的位移矩阵为
(10。１)
式中，为单元位移矩阵;[Ｎ］为形函数矩阵;为单元节点位移列阵。
根据几何方程，对位移矩阵求偏导数，可以得到应变矩阵
（10．2）
式中,[B］为连续单元节点位移和单元应变的矩阵，，[B]为常数矩阵,元素值取决于单元节点坐标差.
根据本构方程，可以得到单元节点位移与单元应力矩阵之间的关系
（10。3）
式中,［Ｄ］为弹性矩阵。应用虚功原理和最小势能原理可以推导出单元刚度矩阵的表达式
(10．4）
各单元的体积力和面力按照静力的等效原则移置到各单元的节点上，其等效节点力为
（1０。5）
式中,{Pe}为作用于单元体积