采矿工程问题数值模拟研究与分析.docVIP

第一章 绪论 随着实测技术、非连续介质力学、大变形理论、物理模拟技术和计算机数值模拟的发展，使得人们可以从理论分析、现场实测、物理模拟和计算机模拟等不同侧面来研究采矿岩层控制中各种力学机理及围岩控制技术。然而现场实测需要大量人力、物力，而且耗时较长，物理模拟的结果虽然比较直观，但受模型尺寸、信息提取和处理技术以及成本和时间的制约，难于分析各个因素的影响规律，而且只能进行定性分析，现有的力学理论只能解决圆形或者椭圆形巷道或硐室等具有简单形体的问题，对于几何形状复杂的巷道或硐室则无能为力。这些问题在数值模拟中都能得到很好的解决。数值力学分析正是为了解决这些问题应运而生的，它不仅能模拟岩体的复杂力学和结构特性，也可很方便地分析各种边值问题和施工工艺过程对硐室或巷道围岩稳定性的影响，并对工程岩体稳定性进行预测和预报。如果能从宏观上把握岩体的力学特性，通过地应力测试把握地应力场，数值力学分析结果完全可以用于指导工程实践。正因为上述各个方面的原因，数值模拟技术得到了大力发展，已成为解决采矿工程和其他岩土工程问题的重要研究手段之一。 第一节 数值模拟发展及其现状 一、数值模拟方法的发展 近几十年来，随着计算机应用的发展，数值计算方法在岩土工程问题分析中迅速得到了广泛应用，大大推动了岩土力学的发展。在岩土力学中所用的数值方法主要有：有限差分法、有限元法、边界元法、加权余量法、半解析元法、刚体元法、非连续变形分析法、离散元法、无界元法和流形元方法等。 有限差分法比较古老，在计算机出现以前就有。它的基本思想是将待解决问题的基本方程和边界条件(一般为微分方程)近似地用差分代数方程来表示，即由有一定规则的空间离散点处的场变量(应力、位移)的代数表达式代替，从而把求解微分方程的问题转化为求解代数方程的问题。也就是说，它将实际的物理过程在时间和空间上离散化，分成有限数量的有限差分量，近似假设这些差分量足够小，以致在差分量的变化范围内物体的性能和物理过程都是均匀的，并且可以应用描述物理现象的定律，只不过在差分量之间发生阶跃式变化。其原理就是将实际连续的物理过程离散化，近似地置换成一连串的阶跃过程，用函数在一些特定点的有限差商代替微商，建立与原微分方程相应的差分方程，从而将微分方程转化为一组代数方程，且有限差分法通常采用“显式”时间步进方法解算代数方程。有限差分法原理简单，可以处理一些相当复杂的问题，应用范围很广。 20世纪50年代出现并得到广泛应用的有限元法，使经典力学解析方法难以解决的工程力学问题都可以用有限元方法求解。它是基于最小总势能变分原理，能方便地处理各种非线性问题，能灵活地模拟岩土工程中复杂的施工过程，因而是目前工程技术领域中实用性最强、应用最为广泛的数值模拟方法，它将连续的求解域离散为一组有限个单元的组合体，解析地模拟或逼近求解区域。由于单元能按各种不同的联结方式组合在一起，且单元本身又可有不同的几何形状，所以可以适应几何形状复杂的求解域。它的特点是利用每个单元内假设的近似函数来表示求解区域上待求的未知场函数，单元内的近似函数由未知场函数在各个单元节点上的数值以及插值函数表达。这就使未知场函数的节点值成为新未知量，把一个连续的无限自由度问题变成离散的有限自由度问题。只要解出节点未知量，便可以确定单元组合体上的场函数，随着单元数目的增加，近似解收敛于精确解。按所选未知量的类型，即是以节点位移作为基本未知量还是以节点力作为基本未知量、或二者皆有，有限元法可分为位移型、平衡型和混合型。由于位移型有限元法在计算机上更易实现，便于电算求解，且易推广到非线性和动力效应等其他方面，所以，位移型有限元法比其他类型的有限元法应用更为广泛。比较著名的通用有限元程序有SAP5、 ADINA等。有限元法需形成总体刚度矩阵，常常需要巨大的存储容量，甚至大得无法计算。由于相邻界面上只能位移协调，对于奇异性问题(应力出现间断)的处理比较麻烦，这是有限单元法的不足之处。 边界元法是20世纪60年代发展起来的求解边值问题的一种数值方法。它是以拜特(Betti)互等定理为基础，有直接法与间接法两种。直接边界元法是以互等定理为基础建立起来的，而间接边界元法是以叠加原理为基础建立起来的。它把边值问题归结为求解边界积分方程问题，在边界上划分单元，求边界积分方程的数值解，进而求出区域内任意点的场变量，故又称为边界积分方程法。边界元法只需对边界离散和积分，与有限元法相比，具有降低维数(将三维降为二维、将二维降为一维)、输人数据较简单、计算工作量少、精度较高等优点。这一点在无限域或半无限域问题的求解中表现得尤为明显，尤其是对均质成等效均质围岩地下工程问题的分析更为方便。边界元法的基本解本身就有奇异性，可以比较方便地处理所谓奇异性问题，因此，边界元法得到人们的重视。目前还有一种方法将边