数值分析方法（flac）.doc

连续介质快速拉格朗日分析(FLAC)方法 本文档目录： §1 HYPERLINK 连续介质快速拉格朗日差分法法概述 1.1 HYPERLINK FLAC方法的提出 1.2 HYPERLINK FLAC方法的发展现状 1.3 HYPERLINK FLAC特点 §2 HYPERLINK FLAC基本原理 2.1 HYPERLINK FLAC-2D基本原理 2.1.1 HYPERLINK 差分方程 2.1.2 HYPERLINK 运动方程和位移计算 2.1.3 HYPERLINK 应力应变的求解 2.1.4 HYPERLINK 失衡力的求解 2.2 HYPERLINK FLAC-3D基本原理 2.2.1 HYPERLINK 空间导数的有限差分近似 2.2.2 HYPERLINK 运动方程 2.2.3 HYPERLINK 阻尼力 2.2.5 HYPERLINK 计算循环 §3 HYPERLINK FLAC数值模拟软件综述 3.1 HYPERLINK FLAC软件概述 3.2 HYPERLINK FLAC软件的主要特点 3.3 HYPERLINK FLAC软件的应用范围 3.4 HYPERLINK 基本术语 3.5 HYPERLINK 软件设计的基本算法与算法流程 3.6 HYPERLINK FLAC3D概述 §4 HYPERLINK FLAC方法优缺点 HYPERLINK 参考文献 §1 连续介质快速拉格朗日差分法法概述 1.1 FLAC方法的提出 近几年来,由于计算机技术的迅猛发展,使得各种数值模拟方法在岩土工程中获得了长足的进步。连续介质快速拉格朗日差分法(Fast Lagrangian Analysis of Continua，简写FLAC) 正是基于该前提从众多方法中脱颖而出的一种新型数值分析方法，该方法是以岩石力学理论为基础，以介质物理力学参数和地质构造特性为计算依据，建立在客观反映原型(地质体的几何形态与物理状态)和仿真其动态演化过程力学效应基础上的一种新型数值方法(谢和平等，1999)，近年来逐步得以成熟完善。 连续介质快速拉格朗日差分法与通常的有限元法和边界元法的不同之处在于：前者是显式的方法，而后者则是隐式的。显式差分法求解时未知数集中在方程式的一边，无须形成刚度矩阵，不用求解大型联立方程，因而占用内存少，便于微机求解较大的工程问题。另外，由于该方法采用的是随流观察法，即研究每个流体质点随时间变化的情况，所以适合于解决非线性大变形问题，这一点恰好可以满足岩土工程中力学分析的需要。 首先将显式拉格朗日差分法程序化的是美国的ITASCA咨询集团公司，该公司利用其原理于20世纪80年代开发了专门求解岩土力学问题的大型商用显式有限差分程序，并于20世纪90年代得以广泛利用。 该程序采用拉格朗日算法及混合离散划分单元技术，应用节点位移连续的条件，使得该程序能够精确地模拟材料的塑性流动和破坏。FLAC内部含有多个力学模型，如摩尔-库伦模型、应变硬化/软化模型、节理模型及双屈服模型等，用以模拟高度非线性、不可拟等地质材料的变形特征。除此之外，FLAC还有多种特殊功能：含有界面单元,可以模拟岩层中的不连续面，如断层、节理及层理等滑动和离层。它不但能处理一般的大变形问题，而且可以在模型中加入节理、弱面等地质构造，同时还可以模拟锚固工程中经常用到的锚杆(索) 设置、喷射混凝土等加固措施。可以模拟地应力场生成、边坡或地下工程的开挖、混凝土衬砌、地下渗流等。FLAC程序自90年代初引入国内后,已广泛应用于地下采矿、地基工程和边坡稳定性分析等多个领域。 1.2 FLAC方法的发展现状 FLAC方法自提出以后，在国外已被广泛应用于工程地质、岩土力学以及构造地质学和成矿学等研究领域。国内于20世纪90年代初才引进该软件，主要应用于工程地质和岩土力学分析，如矿体滑坡、煤矿开采沉陷预测、水利枢纽岩体稳定性分析、采矿巷道稳定性研究等。 随着计算机技术的迅速发展，数值模拟仿真分析方法已经广泛用于岩土及地下工程的研究和设计中。FLAC3D(Fast Lagrangian Analysis of Continuain 3 Dimensions)是美国Itasca咨询公司于20世纪90年代中期在原有的二维分析软件FLAC基础上开发发行的一种工程计算处理程序。 FLAC3D不仅包括了所有FLAC的功能，并在该程序基础上进一步开发，能够模拟计算三维岩土体及其他介质中工程结构的受力与变形性态，计算过程中允许材料发生屈服及流变，适合于解决岩土工程中经常遇到的大变形。与FLAC类似，基于拉格朗日元法的FLAC3D适用于绝大