FLAC3D的基本知识介绍

FLAC3D的基本知识介绍岩土工程结构的数值解是建立在满足基本方程（平衡方程、几何方程、本构方程）和边界条件下推导的。

由于基本方程和边界条件多以微分方程的形式出现，因此，将基本方程近假发改用差分方程（代数方程）表示，把求解微分方程的问题改换成求解代数方程的问题，这就是所谓的差分法。

差分法由来已久，但差分法需要求解高阶代数方程组，只有在计算机的出现，才使该法得以实施和发展。

FLAC3D（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。

目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存（64K），所以，程序求解的最大结点数仅限于2000个以内。

1995年，FLAC2D已升级为V5.0的版本，其程序能够使用护展内存。

因此，大大发护展了计算规模。

FLAC3D是一个三维有限差分程序，目前已发展到V3.0版本。

并且其推出的FLAC SLOPE有了WINDOWS界面。

FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。

因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。

FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。

调整三维网格中的多面体单元来拟合实际的结构。

单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发变形和移动（大变形模式）。

FLAC3D采用的显式拉格朗日算法和混合-离散分区技术能够非常准确发模拟材料的塑性破坏和流动。

由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。

FLAC3D采用ANSI C++语言编写的。

FLAC3D有以下几个优点：1 对模拟塑性破坏和塑性流动采用的是“混合离散法“。

这种方法比有限元法中通常采用的“离散集成法“更为准确、合理。

2 即使模拟的系统是静态的，仍采用了动态运动方程，这使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍。

3 采用了一个“显式解“方案。

因此，显式解方案对非线性的应力-应变关系的求解所花费的时间，几乎与线性本构关系相同，而隐式求解方案将会花费较长的时间求解非线性问题。

面且，它没有必要存储刚度矩阵，这就意味着；采用中等容量的内存可以求解多单元结构；模拟大变形问题几互并不比小变形问题多消耗更多的计算时间，因为没有任何刚度矩阵要被修改。

当然，它也存在以下几个不足之处：1 对于线性问题的求解，FLAC3D比有限元程序运行得要慢；因此，当进行大变形非线性问题或模拟实际可能出现不稳定问题时，FLAC3D是最有效的工具。

2 用FLAC3D求解时间取决于最长的自然周期和最短的自然周期之比。

但某些问题对模型是无效的。

Flac3D中为岩土工程问题的求解开发了特有的本构模型，总共包含了10种材料模型：1. 开挖模型null ；2. 3个弹性模型（各向同性，横观各向同性和正交各向同性弹性模型）3 6个塑性模型（Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模型）。

现在版本的已经有14种本构：1.Elastic-isotropic,othotropic and transversely isotropic2.Drucker-Prager3.Mohr Coulomb4.Ubiquition Joint5.Strain hardening/softening6.Bilinear strain hardening/softening&ubiquitous joint7.Double Yield8.Modified Cam Clay(critical state)9.Hoek-Brown10.Classical viscoelasticity11.2-component power law12.Reference creep formulation(WIPP)蠕变模型13.Burger viscoplasticity14.WIPP viscoplasticity15.Crushed salt creep modelFlac3D网格中的每个区域可以给以不同的材料模型，并且还允许指定材料参数的统计分布和变化梯度。

还包含了节理单元，也称为界面单元，能够模拟两种或多种材料界面不同材料性质的间断特性。

节理允许发生滑动或分离，因此可以用来模拟岩体中的断层、节理或摩擦边界。

FLAC3D中的网格生成器gen，通过匹配、连接由网格生成器生成局部网格，能够方便地生成所需要的三维结构网格。

还可以自动产生交叉结构网格（比如说相交的巷道），三维网格由整体坐标系x,y,z系统所确定，不同于FLAC程序是由行列方式确定。

这就提供了比较灵活的产生和定义三维空间参数。

边界条件和初始条件：定义方式与FLAC相同。

在边界区域可以指定速度（位移）边界条件或应力（力）边界条件。

也可以给出初始应力条件，包括重力荷载以及地下水位线。

所有的条件都充许指定变化梯度。

FLAC3D还包含了模拟区域地下水流动、孔隙水压力的扩散以及可形的多孔隙固体和在孔隙内粘性流动流体的相互耦合。

流体被认为是服从各向同性的达西定律。

流体和孔隙固体中的颗粒是可变形的，将稳态流处理为紊态流可以模拟非稳态流。

同时能够考虑固定的孔隙压力和常流的边界条件，也能模拟源和井。

流体模型可以与结构的力学分析独立进行。

FLAC3D做计算分析的一般步骤：与大多数程序采用数据输入方式不同，FLAC采用的是命令驱动方式。

命令字控制着程序的运行。

在必要时，尤其是绘图，还可以启动FLAc用户交互式图形界面。

为了建立FLAC计算模型，必须进行以下三个方面的工作：1. 有限差分网格2. 本构特性与材料性质3. 边界条件与初始条件完成上述工作后，可以获得模型的初始平衡状态，也就是模拟开挖前的原岩应力状态。

然后，进行工程开挖或改变边界条件来进行工程的响应分析，类似于FLAC的显式有限差分程序的问题求解。

与传统的隐式求解程序不同，FLAC采用一种显式的时间步来求解代数方程。

进行一系列计算步后达到问题的解。

在FLAC中，达到问题所需的计算步能够通过程序或用户加以控制，但是，用户必须确定计算步是否已经达到问题的最终的解。

传说中的FLAC2.1版本的新特性在这里：Summary of New Features in FLAC3D Version 2.1：1. New fluid-flow logic with phreatic surface calculation2. Addition of double-yield constitutive model for volumetric yielding3. User-defined constitutive models in C++, loaded at run-time4. Conversion to Visual C++ compiler5. Linking FLAC3D to other Itasca codes via TCP/IP connection6. New printing/output facilities (BMP, JPEG, DXF, PCX)7. New SOLVE FOS command automatically performs factor-of-safety calculations for slopes and similar problems8. Modification to Finn model for dynamic pore pressure generation9. Implementation of Von Neumann and Landshoff artificial viscosity terms10. Addition of liner and geogrid structural elements11. Improvement of calculation for stress resultants in Structural Elements upgrade to facilitate use of structural elements, including additional example applications in manual12. Extension to pile elements to simulate effect of rockbolt reinforcement13. Additional FISH functions to assist access to FLAC3D variables14. Revised and expanded user's manual with several new examples and verification problems15. New creep models:Burger viscoelastic and viscoplastic models and power-law viscoplastic modelFLAC 软件的主要特点(1) 对硬件配置较低由于FLAC 采用的是显示有限差分法, 在内存较小的低档机上亦可进行较大规模的计算, 计算机内存分别为2、4、8、16MB 时, FL A C 可计算的单元数分别为1500、10000、30000、60000(2) 强大的前后处理功能FLAC 具有很强的前后处理功能。

只要设置某些控制点的坐标, 软件就以自动生成计算网格, 界面友好且美观。

用户可以根据实际情况通过某些命令修改网格, 如对于圆形巷道可采用全放射性网格; 对于其它非规则峒室及复杂地下洞群, 可采用局部密集周边疏松的网格。

各阶段的计算结果均可以数据文件的形式存盘, 一旦需要就可用Restart 命令恢复全部现场, 使用起来非常方便。

输出的图形包括各个施工期的主应力向量图, 各应力分量, 以及位移的等值线图, 锚杆锚索) 以及结构的受力向量图, 塑性区范围等, 并且可用14 种颜色显示, 在V GA 以上显示器上可显出色彩鲜艳, 表现力极强的图形。

可通过彩色打印机、拍照和摄象等手段获得信息量极大的各种图形。