05第五章 岩石力学数值分析方法.doc

第五章 岩石力学数值分析方法 5.1 概 述 岩体不仅为一般材料，更重要的是一种地质结构体，它具有非均质、非连续、非线性以及复杂的加卸载条件和边界条件。这使得岩石力学问题通常无法用解析方法简单地求解。相比之下，数值法具有较广泛的适用性。它不仅能模拟岩体的复杂力学与结构特性，也可很方便地分析各种边值问题和施工过程，并对工程进行预测和预报。因此，岩石力学数值分析方法是解决岩土工程问题的有效工具之—。 岩石力学数值分析方法主要用于研究岩土工程活动和自然环境变化过程中岩体及其加固结构的力学行为和工程活动对周围环境的影响 目前较为常用的方法有：有限元法、边界元法、有限差分法、加权余量法、离散元法、刚体元法、不连续变形分析法、流形方法等其中前四种方法是基于连续介质力学的方法，随后的三种方法则是基于非连续介质力学的方法，而最后一种方法具有这两大类方法的共性。 有限元法基于最小总势能变分原理[1，6~9]，以其能方便地处理各种非线性问题，能灵活地模拟岩土工程中复杂的施工过程，因而成为岩石力学领域中应用最广泛的数值分析方法。 边界元法以表述拜特（Betti）互等定理的积分方程为基础[2]，建立了直接法的基本方程，而基于叠加原理建立了间接法的总体方程；因其前处理工作量少、能有效模拟运场效应而普遍应用于无界域或半无界域问题的求解。 有限差分法是将问题的基本方程和边界条件以简单、直观的差分形式来表述，使得其更易于在工程实际中应用。尤其是近年来FIAC程序在国内外的广泛应用，使得有限差分法在解决岩石力学问题时又获得了新生。 离散单元法是康德尔（Cundall）（1971）[3]以刚性离散单元为基本单元，根据牛顿第二定律，提出的一种动态分析方法。随后又将其发展为变形离散单元（简单变形离散单元和充分变形离散单元），使之既能模拟块体受力后的运动；又能模拟块体本身受力变形状态。自20世纪80年代中期引入我国后，在边坡工程、隧道工程、采矿工程及基础工程等方面都有重要应用。 流形方法是由石根华等人近期（1992，1996）发展的一种新的数值分析方法。这种方法以拓扑学中的拓扑流形和微分流形为基础，在分析域内建立可相互重叠、相交的数学覆盖和覆盖材料全域的物理覆盖，在每一物理覆盖上建立独立的位移函数（覆盖函数），将所有覆盖上的独立覆盖函数加权求和，即可得到总体位移函数。然后，根据总势能最小原理,建立可以用于处理包括非连续和连续介质的耦合问题、小变形、大位移、大变形等多种问题的求解格式[4，5]。它是一种具有一般形式的通用数值分析方法，有限元法和不连续变形分析法（DDA）都可以看作是它的特例。 还有不少其他的数值分析方法［6］和各种数值方法的相互耦合方法。同时，由于实际岩土具有模糊性、随机性和各种不确定性，因此，非确定性数值分析方法也有长足的发展。 所有这些数值分析方法的共同特点是：将带有边值条件的常微分方程或偏微分方程离散为线形代数方程组，采用适当的求解方法解方程组，获得基本未知量，进而根据几何方程和物理方程，求出研究范围内的其他未知量。 对于岩石力学与工程问题，还应重视随工程活动和时间变化岩土体及加固结构中的节理、裂隙等各种不连续面及其扩展过程和自然环境变化过程的数值仿真模拟。与此同时，应根据实际岩体的赋存环境和工程因素影响情况，确定用于数值分析的计算力学参数。因此，利用工程现场的实测位移，进行唯一反分析法来获取数值分析的力学参数[7，8]，将仍是数值分析获取较符合实际结果的有效方法之一。 鉴于篇幅所限，本章仅简要介绍有限元法、边界元法、有限差分法、离散元法和位移反分析法的基本内容。 5.2 有限元法 有限元法是目前已广为应用的岩土工程与结构分析的有力工具。该法是把一个实际的结构物或连续作用一种由多个彼此相联系的单元体所组成的近似等价物理模型来代替。通过结构及连续体力学的基本原理及单元的物理特性建立起表征力和位移关系的方程组。解方程组求其基本未知物理量，并由此求得各单元的应力、应变以及其他辅助量值。有限元法按其所选未知量的类型，即以节点位移作为基本未知量，还是以节点力作为基本未知量，或二者皆有，可分为位移型的、平衡型的和混合型的有限元法。由于位移型有限元法在计算机上更易实现复杂问题的系统化，且便于电算求解，更易推广到非线性和动力效应等其他方面。所以，位移型有限元法比其他类型的有限元法应用更为广泛。因此，本节以平面三节点三角形单元为例，简单介绍位移型有限元法的基本方程。 5.2.1 有限元法的