FLAC讲义.ppt

本课程主要讲授内容 1. 绪论 2. 数值模拟方法：FLAC 3. 神经元网络（自学） 4. 采矿专家系统（自学） 5.矿山信息技术（矿山生产过程中的监测与监控、虚拟现实技术） 什么是数值模拟？ 数值模拟方法是通过对结构系统构造数学模型，利用计算机求解大规模的代数联立方程组来模拟结构系统的反映过程的。 数值方法的突出优点就是能够较好地考虑诸如介质的各向异性、非均质特性及其随时间的变化、复杂边界条件和介质不连续等复杂地质条件。 在地下工程问题分析中，常用的数值方法有：有限元、差分法、加权残位法、边界元、离散元、拉格朗日元法等。 2.1.1 什么是FLAC？ FLAC其实是快速拉格朗日分析 (Fast Lagrangian Analysis of Continua)的简称，是一种用于工程力学计算的显式有限差分程序。该程序可模拟土、岩石等材料的力学行为。它利用拉格朗日有限差分原理编制，以牛顿第二定律为基础，可以对连续介质进行非连续大变形分析，特别适用于求解岩土工程中的大变形问题。 2.1 .2 FLAC的发展 (1).由美国ITASCA咨询集团公司开发，首先由Cundall在20世纪80年代提出并将其程序化、实用化。 (2).在20世纪90年代，得以广泛应用。我国20世纪90年代初才引进该软件。 (1)零空模型(Null)：代表网格中的孔洞(开挖单元)； (2)应变硬化/软化模型(SS)：代表非线性，不可逆剪切破碎与压缩； (3)粘弹性蠕变模型； (4)界面模型(界面为平面，沿界面允许滑动和分开)：模拟断层、节理和摩擦边界； (5)水利模型：模拟可变形空隙体与粘性流动的全藕合； (6)结构单元模型：模拟岩土体加固、衬砌、锚杆、混凝土喷层、可缩支柱及钢拱等。 (7)轴对称几何模型：模拟围堰、船闸及层状材料侧向载荷的影响； (8)动态分析模型：其代码能用于各种工程动态问题，诸如地震分析、坝的稳定、主结构间的作用与液化、爆破载荷的影响等； (9)热力模型：模拟材料中的瞬态热流，热应力的发生以及进行热与力的藕合计算等； (10)绘图功能：通过其重复占位程序，用户能绘制各种图形与表格，其中计算时步函数关系曲线的绘制特别有助于弄清楚何时到达平衡与破裂状态，并在瞬态计算(如地下水流计算)或动态计算(如地震运动计算)中进行变量化监控。 2.1.3 FLAC的应用 在国内主要用与岩土力学分析；例如矿体滑坡、煤矿开采沉陷预测、水利枢纽岩体稳定性分析、采矿巷道稳定性研究等。现在已经逐步发展到应用于地质力学问题的研究，主要测定工程地质、构造地质学、大陆动力学。 在国外该软件开始主要泛用于岩土力学。现在主要应用于以下几个方面的研究，如工程地质、构造地质学、大陆动力学、成矿学。 FLAC软件的基本原理即是拉格朗日差分法，是一种利用拖带坐标系分析大变形问题的数值方法。并利用差分格式按时步积分求解。 拉格朗日元法源于流体力学，在流体力学中有两种主要的研究方法。一种是定点观察法，亦称欧拉法；另一种是随机观察法，称为拉格朗日法。后者是研究每个流体质点随时间而变化的状态，即研究一物体质点在一段时间内的运动轨迹、速度、压力等特征。 把拉格朗日元法移植到固体力学中，把所研究的区域划分成网格，其结点就相当于流体质点。然后按时步用拉格朗日元法来研究网格结点的运动。 拉格朗日元法是一种利用拖带坐标系分析大变形问题的数值方法并利用差分格式按时步积分求解。随着构形的不断变化、不断更新坐标，允许介质有较大的变形。模型经过网格划分，物理网格影射成数字网格，数字网格上的某个结点就与物理网格上相应的结点坐标相对应。对于某个结点而言在每个时刻它受到来自其周围区域的合力的影响。如果合力不等于零，结点就具有了失稳力，就要产生运动。假定结点上集中有临接该结点的质量，于是在失稳力的作用下，根据牛顿第二定律，结点就要产生加速度，进而可在一个时步中求得速度和位移的增量。对于每个区域而言，可以根据其周围结点的运动速度求得它的应变率。 然后根据材料的本构关系求得应力增量。由应力增量求出t 和t+△t时刻各个结点的不平衡力和各个结点在 t+△t时的加速度。对加速度进行积分，即可得到结点的新的位移值。由此可以求得各结点的新的坐标值。同时由于物体的变形，单元要发生局部的平均整旋和整旋。只要计算相应的应力改正值，最后通过应力叠加就可以得到新的应力值，到此计算为一个循环，然后按时步进行下一轮的计算。如此一直进行到问题收敛。FLAC程序采用最大不平衡力来刻划FLAC计算的收敛过程。如果单元的最大不平衡力随着时步增加而逐渐趋于极小值，则计算是稳定的，否则，计算就是不稳定的。 上图给出了FLAC显示静态分析求解流程。现说明如下： (1)建立FLAC模型：实施FLAC算法，首先要建立FLAC平面应变模型，包括生成网