岩石力学的数值模拟(讲义).doc

第10章 岩石力学的数值模拟
随着计算机软硬件技术的迅速发展，使岩石力学有了长足的进步，特别在岩石力学的数值计算和模拟方面发展尤为迅速，使得许多岩石力学解析方法难于解决的问题得以重新认识。正如钱学森在给中国力学学会“力学——迎接21世纪新的挑战”的一封信中对力学发展趋势总结的那样“今日力学是一门用计算机计算去回答一切宏观的实际科学技术问题，计算方法非常重要”。岩石力学和其他力学学科一样，需要数值计算方法并推动岩石力学的发展。
岩石介质不同于金属材料，在数值计算方面具有其独特的特点[205]：
（1）岩石介质是赋存于地壳中的各向异性天然介质。
（2）岩石介质被众多的节理、裂缝等弱面所切割而呈现高度的非均质性，而其物理、化学及力学性质具有随机性特点。
（3）岩石介质赋存时以受压为主，而且抗压强度远大于抗拉强度。
（4）岩石力学与工程问题在时空分布上较广，从本质上讲都是三维问题。
（5）岩石工程一般无法进行原型试验，而实验室测得的数据不能直接应用于工程设计和计算。
（6）岩石力学与工程具有数据有限问题。
数值计算方法经过几十年的发展，目前已形成许多种岩石力学计算方法，主要有有限元法、边界元法、有限差分法、离散元法、流形元法、拉格朗日元法、不连续变形法及无单元法等。它们各有优缺点，有限元的理论基础和应用比较成熟，在金属材料和构件的计算中应用十分成功，但它是以连续介质为基础，似乎与岩体的非连续性有一定差距，流形元等数值方法虽然考虑了岩体中节理效应，但其理论基础还不完全成熟。相信在不久的将来，肯定会出现完全适合于岩体材料和工程的数值计算方法[206~208]。
岩石力学的有限元分析[209~213]
有限元法（finite element method，FEM）是岩石力学数值计算方法中最为广泛应用的一种。自20世纪50年代发展至今，有限元已成功地求解了许多复杂的岩石力学与工程问题。被广大岩石力学研究与工程技术人员喻为解决岩石工程问题的有效工具。有限元法是根据变分原理求解数学物理方程的一种数值方法。有限元法把连续体离散成有限个单元，每个单元的场函数只包含有限个节点参量的简单场函数，这些有限个单元的场函数集合构成整个结构连续体场函数。根据能量方程和加权函数方程可建立有限个求解参数的方程组，求解这些离散方程组，就是有限元法的精髓所在。虽然求解时把连续函数转化为求解有限个离散点处的函数值，但只要单元划分得充分小时，足可以满足计算要求。
有限元法求解问题时一般遵循以下步骤：
（1）有限元计算模型的建立，包括模型单元的划分、确定边界条件。
（2）对单元体进行力学分析，包括求解节点位移、单元应变和单元应力。
（3）对计算模型进行分析。
（4）进行计算分析。
线弹性有限元法的基本方程
线弹性有限元是弹塑性有限元、损伤有限元、流变有限元等非线性有限元的基础。线弹性有限元假定岩石介质连续、均质、小变形和完全弹性。
有限元法求解弹性力学问题时通常以位移作为基本未知量，单元位移是以单元节点位移为基本未知量，选择合理的位移插值函数，将单元位移表达为节点坐标的连续函数，插值函数也可称为形函数。不同形状的单元具有不同的形函数。
图10-1为三种最常见单元形式，即三角形、四边形及四面体单元。它们的形函数分别为：
图10-1 有限元的三种基本单元形式
（a）三角形单元（b）四边形单元（c）四面体单元
三角形的形函数
式中，S为三角形面积；；；。
四边形的形函数
式中，位移量为；；；。
四面体的形函数


式中，V为四面体的体积。
单元在直角坐标轴中位移分量分别为u，v，，因此单元的位移矩阵为
（）
式中，为单元位移矩阵；[N]为形函数矩阵；为单元节点位移列阵。
根据几何方程，对位移矩阵求偏导数，可以得到应变矩阵
（）
式中，[B]为连续单元节点位移和单元应变的矩阵，也称为应变矩阵。对于三角形单元，[B]为常数矩阵，元素值取决于单元节点坐标差。
根据本构方程，可以得到单元节点位移与单元应力矩阵之间的关系
（）
式中，[D]为弹性矩阵。应用虚功原理和最小势能原理可以推导出单元刚度矩阵的表达式
（）
各单元的体积力和面力按照静力的等效原则移置到各单元的节点上，其等效节点力为
（）
式中，{Pe}为作用于单元体积力{P}的等