FLAC3D简介(word文档良心出品).docxVIP

FLAC3D 知识基本介绍 SimWe 岩土工程结构的数值解是建立在满足基本方程 （平衡方程、几何方程、 本构方程） 和边 界条件下推导的。 由于基本方程和边界条件多以微分方程的形式出现， 因此， 将基本方程近 假发改用差分方程（代数方程）表示，把求解微分方程的问题改换成求解代数方程的问题， 这就是所谓的差分法。 差分法由来已久， 但差分法需要求解高阶代数方程组， 只有在计算机 的出现，才使该法得以实施和发展。 FLAC3D （Fast Lagrangian Analysis of Continua ）由美国 Itasca 公司开发的。 目前， FLAC 有二维和三维计算程序两个版本，二维计算程序 V3.0 以前的为 DOS 版本， V2.5 版本仅仅 能够使用计算机的基本内存 （ 64K ），所以，程序求解的最大结点数仅限于 2000 个以内。1995 年， FLAC2D 已升级为 V3.3 的版本，其程序能够使用护展内存。因此，大大发护展了计算 规模。 FLAC3D 是一个三维有限差分程序，目前已发展到 V2.1 版本。 FLAC3D 的输入和一般的数值分析程序不同， 它可以用交互的方式， 从键盘输入各种命 令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。因此，采用 FLAC 程序 进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排 列，形成可以完成一定计算任务的命令文件。 FLAC3D 是二维的有限差分程序 FLAC2D 的护展， 能够进行土 质、 岩石和其它材料的 三维结构受力特性模拟和塑性流动分析。调整三维网格中的多面体单元来拟合实际的结构。 单元材料可采用线性或非线性本构模型， 在外力作用下， 当材料发生屈服流动后， 网格能够 相应发变形和移动（大变形模式）。FLAC3D采用的显式拉格朗日算法和混合 -离散分区技术 能够非常准确发模拟材料的塑性破坏和流动。 由于无须形成刚度矩阵， 因此， 基于较小内存 空间就能够求解大范围的三维问题。 FLAC3D 采用 ANSI C++ 语言编写的。 FLAC3D 有以下几个优点： 对模拟塑性破坏和塑性流动采用的是 “混合离散法 “。这种方法比有限元法中通常采 用的 “离散集成法 “更为准确、合理。 即使模拟的系统是静态的，仍采用了动态运动方程，这使得 FLAC3D 在模拟物理上 的不稳定过程不存在数值上的障碍。 采用了一个 “显式解 “方案。因此，显式解方案对非线性的应力 -应变关系的求解所花 费的时间，几互与线性本构关系相同， 而隐式求解方案将会花费较长的时间求解非线性问题。 面且，它没有必要存储刚度矩阵，这就意味着；采用中等容量的内存可以求解多单元结构； 模拟大变形问题几互并不比小变形问题多消耗更多的计算时间， 因为没有任何刚度矩阵要被 修改。 当然，它也存在以下几个不足之处： 对于线性问题的求解， FLAC3D 比有限元程序运行得要慢；因此，当进行大变形非 线性问题或模拟实际可能出现不稳定问题时， FLAC3D 是最有效的工具。 用 FLAC3D 求解时间取决于最长的自然周期和最短的自然周期之比。 但某些问题对 模型是无效的。 Flac3D 中为岩土工程问题的求解开发了特有的本构模型，总共包含了 10 种材料模型： 开挖模型 null 3 个弹性模型（各向同性，横观各向同性和正交各向同性弹性模型） 6 个塑性模型（ Drucker-Prager 模型、 Morh-Coulomb 模型、应变硬化 /软化模型、遍 布节理模型、双线性应变硬化 /软化遍布节理模型和修正的 cam 粘土模型）。 Flac3D 网格中的每个区域可以给以不同的材料模型，并且还允许指定材料参数的统计 分布和变化梯度。 还包含了节理单元， 也称为界面单元， 能够模拟两种或多种材料界面不同 材料性质的间断特性。 节理允许发生滑动或分离， 因此可以用来模拟岩体中的断层、 节理或 摩擦边界。 FLAC3D中的网格生成器gen,通过匹配、连接由网格生成器生成局部网格，能够方便地生成所需要的三维结构网格。还可以自动产生交叉结构网格（比如说相交的巷道） ，三维 网格由整体坐标系 x,y,z 系统所确定,不同于 FLAC 程序是由行列方式确定。这就提供了比 较灵活的产生和定义三维空间参数。 边界条件和初始条件： 定义方式与 FLAC 相同。在边界区域可以指定速度（位移）边界条件或应力（力）边 界条件。 也可以给出初始应力条件, 包括重力荷载以及地下水位线。 所有的条件都充许指定 变化梯度。 FLAC3D 还包含了模拟区域地下水流动、 孔隙水压力的扩散以及可形的多孔隙固体和在 孔隙内粘性流动流体的相互耦合。 流体被认为是服从各向同性的达西定律。 流体和孔隙固体 中的颗粒是可变形的