FLAC数值模拟介绍

1、标准FLAC-3D(Three Dimensional Fast Lagrangian Analysis of Continua) 是美国 Itasca Consulting Goup lnc开发的三维快速拉格朗日分析程序，该程序能较好地模拟地质材料在达到强度极限或屈服极限时，发生的破坏或塑性流动的力学 行为，特别适用于分析渐进破坏和失稳以及模拟大变形.FLAC3D分析的使用领域根据手册总结如下：(1)承受荷载能力与变形分析：用于边坡稳定和基础设计(2)渐进破坏与坍塌反演：用于硬岩采矿和隧道设计(3)断层构造的影响研究：用于采矿设计(4)施加于地质体锚索支护所提供的支护力研究：岩锚和土钉的设计

2、(5)排水和不排水加载条件下全饱和流体流动和孔隙压力扩散研究：挡土墙结构的地下水流动，和土体固结研究(6)粘性材料的蠕变特性：用于碳酸钾盐矿设计(7)陡滑面地质结构的动态加载：用于地震工程和矿山岩爆研究(8)爆炸荷载和振动的动态响应：用于隧道开挖和采矿活动(9)结构的地震感应：用于土坝设计(10)由于温度诱发荷载所导致的变形和结构的不稳定(11)大变形材料分析：用于研究粮仓谷物流动和放矿的矿石流动10种材料本构模型文案标准Flac3D中为岩土工程问题的求解开发了特有的本构模型，总共包含了 10种材料 模型：开挖模型null3个弹性模型(各向同性，横观各向同性和正交各向同性弹性模型)6个塑性模型

3、(Drucker-Prager 模型、Morh-Coulomb 模型、应变硬化/软 化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土 模型).Flac3D网格中的每个区域可以给以不同的材料模型，并且还允许指定材料参数 的统计分布和变化梯度.还包含了节理单元，也称为界面单元，能够模拟两种或 多种材料界面不同材料性质的间断特性.节理允许发生滑动或分离，因此可以用 来模拟岩体中的断层、节理或摩擦边界.FLAC3D中的网格生成器gen,通过匹配、连接由网格生成器生成局部网格，能 够方便地生成所需要的三维结构网格.还可以自动产生交岔结构网格(比如说相 交的巷道),三维网格由整体坐

4、标系x, y, z系统所确定，这就提供了比较灵活的产 生和定义三维空间参数.五种计算模式(l)静力模式：这是FLAC-3D默认模式，通过动态松弛方法得静态解.(2)动力模式:用户可以直接输人加速度、速度或应力波作为系统的边界条件或初 始条件，边界可以固定边界和自由边界.动力计算可以与渗流问题相藕合.(3)蠕变模式:有五种蠕变本构模型可供选择以模拟材料的应力-应变-时间关系:Maxwell模型、双指数模型、参考蠕变模型、粘塑性模型、脆盐模型.文案标准(4)渗流模式:可以模拟地下水流、孔隙压力耗散以及可变形孔隙介质与其间的粘 性流体的耦合.渗流服从各向同性达西定律，流体和孔隙介质均被看作可变形体

5、考虑非稳定流，将稳定流看作是非稳定流的特例.边界条件可以是固定孔隙压力 或恒定流，可以模拟水源或深井.渗流计算可以与静力、动力或温度计算耦合， 也可以单独计算.(5)温度模式:可以模拟材料中的瞬态热传导以及温度应力.温度计算可以与静 力、动力或渗流计算藕合，也可单独计算.模拟多种结构形式(l)对于通常的岩体、土体或其他材料实体，用八节点六面体单元模拟.(2) FIAC-3D包含有四种结构单元：梁单元、锚单元、桩单元、壳单元.可用来模 拟岩土工程中的人工结构如支护、衬砌、锚索、岩栓、土工织物、摩擦桩、板桩等.(3) FLAC-3D的网格中可以有界面，这种界面将计算网格分割为若干部分，界面 两边的

6、网格可以分离，也可以发生滑动，因此，界面可以模拟节理、断层或虚拟 的物理边界.有多种边界条件边界方位可以任意变化，边界条件可以是速度边界、应力边界，单元部可以给定 初始应力，节点可以给定初始位移、速度等，还可以给定地下水位以计算有效应 力、所有给定量都可以具有空间梯度分布.FLAC-3D 嵌语言 FISH文案标准FLAC-3D具有强大嵌语言FISH,使得用户可以定义新的变量或函数,以适应用 户的特殊需要,例如，利用HSH做以下事情：(l)用户可以自定义材料的空间分布规律，如非线性分布等.(2)用户可以定义变量，追踪其变化规律并绘图表示或打印输出.(3)用户可以自己设计FLAC-3D部没有的单元

7、形态.(4)在数值试验中可以进行伺服控制.(5)用户可以指定特殊的边界条件.(6)自动进行参数分析(7)利用FLAC-3D部定义的Fish变量或函数，用户可以获得计算过程中节点、 单元参数，如坐标、位移、速度、材料参数、应力、应变、不平衡力等.FLAC-3D前后处理功能FLAC-3D具有强大的自动三维网格生成器，部定义了多种单元形态，用户还可 以利用FISH自定义单元形态，通过组合基本单元，可以生成非常复杂的三维网 格，比如交叉隧洞等.在计算过程中的任何时刻用户都可以用高分辨率的彩色或灰度图或数据文件输 出结果，以对结果进行实时分析，图形可以表示网格、结构以及有关变量的等值 线图、矢量图、曲线

8、图等，可以给出计算域的任意截面上的变量图或等直线图，计 算域可以旋转以从不同的角度观测计算结果.FLAC3D计算分析一般步骤文案标准与大多数程序采用数据输入方式不同，FLAC采用的是命令驱动方式.命令字控 制着程序的运行.在必要时，尤其是绘图,还可以启动FLAc用户交互式图形界面.为了建立FLAC计算 模型，必须进行以下三个方面的工作：(1)有限差分网格(2)本构特性与材料性质(3)边界条件与初始条件完成上述工作后，可以获得模型的初始平衡状态，也就是模拟开挖前的原岩应力状态.然后，进行工程开挖或改变边界条彳来进行工程的响应分析，类似于FLAC的显式有限差分程序的问题求解.与传统的隐式求解程序不

9、同，FLAC采用 一种显式的时间步来求解代数方程.进行一系列计算步后达到问题的解.在FLAC中，达到问题所需的计算步能够通过程序或用户加以控制，但是，用户必须确定计算步是否已经达到问题的最终的解.后处理(一)用tecplot绘制曲线(1)第一主应力(2) xdisp、ydisp、zdisp、disp(二)用excel做曲线隧道(1)做地表沉降槽(zdisp)地表横向位移(xdisp)文案标准(3)隧道中线竖向沉降曲线(zdisp)(4)提取位移矢量图,(5)显示初期支护结构力(6)显示state(找塑性区)基坑(1)做地表沉降槽(zdisp)(2)提取位移矢量图，(3)显示初期支护结构力(4)

10、显示state(找塑性区)边坡(1)做安全系数和应变图模型最优化用FLAC3D解决问题时，为了得到最有效的分析使模型最优化是很重要的.(1)检查模型运行时间：一个FLAC3D例子的运行时间是区域数的4/3倍.这个规则适用于平衡条件下的 弹性问题.对于塑性问题，运行时间会有点改变，但是不会很大，但是如果发生 塑性流动，这个时间将会大的多.对一个具体模型检查自己机子的计算速度很重 要.一个简单的方法就是运行基准测试.然后基于区域数的改变，用这个速度评 估具体模型的计算速度.(2)影响运行时间的因素：FLAC3D有时会需要较长时间才可以收敛主要发生在下列情况下：文案标准（a）材料本身刚度变异或材料与

11、结构及接触面之间的刚度差异很大.（b）划分的区域尺寸相差很大.这些尺寸差异越大编码就越无效.在做详细分析 前应该研究刚度差异的影响.例如，一个荷载作用下的刚性板，可以用一系列顶 点固定的网格代替，并施以等速度.（记住FIX命令确定速度，而不是位移.）地下 水的出现将使体积模量发生明显的增加（流体-固体相互作用）.（3）考虑网格划分的密度：FLAC3D使用常应变单元.如果应力/应变曲线倾斜度比较高，那么你将需要许 多区域来代表多变的分区.通过运行划分密度不同的同一个问题来检查影响. FLAC3D应用常应变区域，因为当用多的少节点单元与用比较少的多节点单元 模拟塑性流动时相比更准确.应尽可能保持网

12、格，尤其是重要区域网格的统一.避免长细比大于5:1的细长单 元，并避免单元尺寸跳跃式变化（即应使用平滑的网格）.应用GENERATE命令中 的比率关键词，使细划分区域平滑过渡到粗划分区域.（4）自动发现平衡状态：默认情况下，当执行SOLVE命令时，系统将自动发现力的平衡.当模型中所有 网格顶点中所有力的平均量级与其中最大的不平衡力的量级的比率小于1*10时,认为达到了平衡状态.注意一个网格顶点的力由力（例如，由于重力）和外力（例如, 由于所加的应力边界条件）共同引起.因为比率是没有尺寸的，所以对于有不同 的单元体系的模型，在大多数情况下，不平衡力和所加力比率的限制给静力平衡 提供了一个精确的限

13、制.同时还提供了其他的比率限制；可以用SET ratio命令施加.如果默认的比率限 制不能为静力平衡提供一个足够精确的限制，那么应考虑可供选择的比率限制.文案标准默认的比率限制同样可用于热分析和流体分析的稳定状态求解 .对于热分析，是 对不平衡热流量和所加的热流量量级进行评估 ，而不是力.对于流体分析，对不 平衡流度和所加流度量级进行评估.(5)考虑选择阻尼:对于静力分析，默认的阻尼是局部阻尼,对于消除大多数网格顶点的速度分量周 期性为零时的动能很有效.这是因为质量的调节过程依赖于速度的改变.局部阻 尼对于求解静力平衡是一个非常有效的计算法则且不会引入错误的阻尼力(见Cundall 1987)

14、.如果在求解最后状态，重要区域的网格海域的速度分量不为零，那么说明默认的 阻尼对于达到平衡状态是不够的.有另外一种形式的阻尼,叫组合阻尼,相比局 部阻尼可以使稳定状态达到更好的收敛，这时网格将发生明显的刚性移动.例如, 求解轴向荷载作用下桩的承载力或模拟蠕变时都可能发生.使用SETmechanical damp combined 命令来调用组合阻尼.组合阻尼对于减小动 能方面不如局部阻尼有效，所以应注意使系统白动力激发最小化.可以用SET mechanical damp local 命令转换到默认阻尼.(6)检查模型反应：FLAC3D显示了一个相试的物理系统是怎样变化的.做一个简单的试验证明你

15、 在做你认为你在做的事情.例如，如果荷载和实体在几何尺寸上都是对称的，当 然反应也是对称的.改变了模型以后，执行几个时步(假如，5或10步)，证明初始 反应是正确的，并且发生的位置是正确的.对应力或位移的期望值做一个估计， 与FLAC3D的输出结果作比较.如果你对模型施加了一个猛烈的冲击，你将会得到猛烈的反应.如果你对模型作文案标准了一些看起来不合理的事情，你一定要等待奇怪的结果.如果在分析的一个给定 阶段，得到了意外值,那么回顾到这个阶段所用的时步.在进行模拟前很关键的是检查输出结果.例如，除了一个角点速度很大外，一切 都很合理，那么在你理解原因前不要继续下去.这种情况下，你可能没有给定适

16、当的网格边界.(7)初始化变量：在模拟基坑开挖过程时，在达到目的前通常要初始化网格顶点位移.因为计算次 序法则不要求位移，所以可以初始化位移，这只是由网格顶点的速度决定，并有 益于用户初始化速度却是一件难事.如果设定网格顶点的速度为一常数，那么这 些点在设置否则前保持不变.所以，不要为了清除这些网格的速度而简单的初始 化它们为零.这将影响模拟结果.然而，有时设定速度为零是有用的(例如，消 除所有的动能).(8)最小化静力分析的瞬时效应：对于连续性静力分析，经过许多阶段逐步接近结果是很重要的.即，当问题条 件突然改变时，通过最小化瞬时波的影响，使结果更加“静力”.使FLAC3D解 决办法更加静态

17、的方法有两种.(a)当突然发生一个变化时(例如，通过使区域值为零模拟开挖),设定强度性能 为很高的值以得到静力平衡.然后为了确保不平衡力很低，设定性能为真实值, 再计算，这样，由瞬时现象引起的失败就不会发生了 .(b)当移动材料时，用FISH函数或表格记录来逐步减少荷载.(9)改变模型材料：FLAC3D对一个模拟中所用的材料数没有限制.这个准则已经尺寸化，允许用文案标准户在自己所用版本的FLAC3D中最大尺寸网格的每个区域（假如设定的）使用不 同的材料.（10）运行在现场原位应力和重力作用下的问题：有很多问题在建模时需要考虑现场原位应力和重力的作用.这种问题的一个例子是深层矿业开挖：回填.此时

18、大多数岩石受很高的原位应力区的影响（即，自重 应力由于网孔尺寸的限制可以忽略不计），但是回填桩的放置使自重应力发展导 致岩石在荷载作用下可能坍塌.在这些模拟中要注意的重点（因为任何一种模拟 都有重力的作用）是网格的至少三个点在空间上应固定.否则，整个网格在重 力作用下将转动.如果你曾经注意到整个网格在重力加速度矢量方向发生转动，那么你可能忘记在空间上固定网格了 .FLAC3D主要适明模拟计算地质材料和岩土上程的力学行为。特别是材料达到屈 服极限后产生的塑性流动、材料通过单元和区域表示，根据计算对象的形状构成 相应的网格。每个单元在外载和边界约束条件下，按照约定的线性或非线性应力 一应变关系产生

19、力学响应。由于 FLAC3D程序主要是为岩上工程应用而开发的 岩石力学计算程序，程序中包括了反映地质材料力学效应的特殊计算功能，可计算地质类材料的高度非线性（包括应变硬化/软化）、不可逆剪切破坏和压密、粘弹 （蠕变）、孔隙介质的应力渗流耦合、热力耦合以及动力学行为等；FLAC3D程序设有多种本构模型： 各向同性弹性材料模型； 横观各向同性弹性材料模型； 莫尔一库仑弹塑材料模型； （4）应变软化/硬化塑性材科模型；文案标准(5)双屈服塑性材料模型；(6)遍布节理材料模型；(7)空单元模型，可用来模拟地下碉室的开挖和煤层开采；另外，程序设有界面单元，可以模拟断层、节理和摩擦边界的滑动、开和闭合行

20、为。支护结构如砌衬、锚杆、可缩性支架或板壳等与围岩的相互作用也可以在 FLAC3D中进行模拟。同时、用户可根据需要在 FLAC3D中创建自己的本构模 型，进行各种特殊修正和补充。FLAC3D程序建立在拉格朗日算法基础上，特别适合模拟大变形和扭曲。FLAC3D采用显式算法来获得模型全部运动方程 (包括变量)的时间步长解，从 而可以追踪材料的渐进破坏和垮落， 这对研究采矿设计是非常重要的。 此外，程 序允许输入多种材料类型，亦可在计算过程中改变某个局部的材料参数，增强了程序使甲的灵活性。极大地方便了在计算上的处理。FLAC3D程序具有强大的后 处理功能，用户可以直接在屏幕上绘制或以文件形式创建和输

21、出打印多种形式的 图形。用户还可根据需要，将若干个变量合并在同一幅图形中进行研究分析。基于上述计算功能与特点。本研究应用 FLAC3D程序计算分析付村煤矿209工作 面山上开采对山下回采巷道围岩变形的影响。FLAC ( Fast Lagrangian Analysis of Continua )是一个利用显式有限差分方法 求解的岩土、采矿工程师进行分析和设计的二维连续介质程序， 主要用来模拟土、 岩、或其他材料的非线性力学行为，可以解决众多有限元程序难以模拟的复杂的 工程问题，例如大变形、大应变、非线性及非稳定系统(甚至大面积屈服/失稳或完全塌方)等问题。文案标准FLAC的基本功能和特征为：?

22、允许介质出现大应变和大变形；? Interface单元可以模拟连续介质中的界面，并允许界面发生滑动和开裂；?显式计算方法，能够为非稳定物理过程提供稳定解， 直观反映岩土体工程 中的破坏；?地下水流动与力学计算完全耦合（包括负孔隙水压，非饱和流及相界面计 算）；?采用结构加固单元模拟加固措施，例如衬砌、锚杆、桩基等；?材料模型库（例如：弹性模型、莫尔库仑塑性模型、任意各向异性模型、 双屈服模型、粘性及应变软化模型）；?预定义材料性质，用户也可增加用户自己的材料性质设定并储存到数据库 中；? 一系列可选择模块，包括：热力学模块、流变模块、动力学模块、二相流 模块等，用户还可用C+建立自己的模型；?

23、边坡稳定系数计算满足边坡设计的要求；?用户可用部语言（FISH）增加自己定义的各种特性（如：新的本构模型， 新变量或新命令）；FLAC软件的优势：?连续体大应变模拟文案标准?界面单元用已代表不连续接触界面可能出现的完全不连续性质的开和滑动，因此可以模拟断层、节理和摩擦边界等?显式求解模式可以获得不稳定物理过程的稳定解?材料模型：o ”空（null ） ” 模型；o三种弹性模型（各向同性、横观各向异性、和正交各向异性）；o 七种非线性模型（Drucker-Prager、Mohr-Coulomb 、应变硬 化及应变软化、 节理化、双线性应变硬化/软化节理化、双屈服、 修正的Cam-clay模型）?

24、任何参数指标的连续变化或统计分布的模拟?外接口编程语言 （FISH）允许用户添加用户自定义功能?方便的边界定义和初始条件定义方式?可定义水位线/面进行有效应力计算?地下水渗流计算以及完全的应力场渗流场偶合计算（含负孔隙压力、非饱和流、井）?结构单元如隧道衬砌、桩、壳、梁锚杆、锚索、土工织物及其组合，可以 模拟不同的加固手段及其与围岩（土体）的相互作用?自选模块包括：o热和热力学分析模块；o流变计算模块；o动力分析模块实现真时间历程的瞬时动力响应模拟；o用C+编写的用户自定义本构模块文案标准开挖直立坡的喷射混凝土墙加土钉加固的模拟加（下）和不加（上）土工织物土坡的潜在破坏特征抗潜越支护蚂头的土文

25、案FLAC3D （Three Dimensional Fast Lagrangian Analysis ofContinua ）是美国 Itasca Consulting Goup lnc开发的三维快速拉格朗日分析程序，是二维的有限差分程序FLAC2D的扩展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动标准分析。该程序能较好地模拟地质材料在达到强度极限或屈服极限时 发生的破坏或塑性流动的力学行为，特别适用于分析渐进破坏和失 稳以及模拟大变形。FLAC3D调整三维网格中的多面体单元来拟合 实际的结构。单元材料可采用线性或非线性本构模型，在外力作用 下，当材料发生屈服流动后，网格能

26、够相应发生变形和移动（大变 形模式）。FLAC3D采用的显式拉格朗日算法和混合-离散分区技 术，能够非常准确的模拟材料的塑性破坏和流动。由于无须形成刚 度矩阵，因此，基于较小存空间就能够求解大围的三维问题。FLAC3D采用ANSI C+语言编写的。它包含10种弹塑性材料本 构模型，有静力、动力、蠕变、渗流、温度五种计算模式，各种模 式间可以互相藕合，可以模拟多种结构形式，如岩体、土体或其他 材料实体，梁、锚元、桩、壳以及人工结构如支护、衬砌、锚索、 岩栓、土工织物、摩擦桩、板桩、界面单元等，可以模拟复杂的岩 土工程或力学问题。FLAC3D的优点：1 .对模拟塑性破坏和塑性流动采用的是“混合离散

27、法”。这种 方法比有限元法常采用的”离散集成法”更为准确、合理。2 .即使模拟的系统是静态的，仍采用了动态运动方程，这使得 FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍。3 .采用了一个“显式解”方案。因此，显式解方案对非线性的应力-应变关系的求解所花费的时间，几互与线性本构关系文案标准相同，而隐式求解方案将会花费较长的时间求解非线性问题。面且，它没有必要存储刚度矩阵，这就意味着， 采用 中等容量的存可以求解多单元结构； 模拟大变形问题几乎 并不比小变形问题多消耗更多的计算时间， 因为没有任何刚 度矩阵要被修改。固fl向3D z.io-ioa Alpha 151 - View soi

28、lpileFile Edit Settings Motitenis Window_ Ifi1! X槽形开挖FLAC3D的功能和特征:1、包含10材料本构模型Flac3D中为岩土工程问题的求解开发了特有的本构模型，总共包含了 10种材料模型：? 开挖模型null ;文案标准? 3个弹性模型（各向同性，横观各向同性和正交各向同性弹 性模型）? 6 个塑性模型（Drucker-Prager 模型、Morh-Coulomb 模型、应变硬化/软化模型、遍布节理模型、双线性应变硬 化/软化遍布节理模型和修正的cam粘土模型）。Flac3D网格中的每个区域可以给以不同的材料模型，并且还允许 指定材料参数的统

29、计分布和变化梯度。还包含了节理单元，也称为界面单元，能够模拟两种或多种材料界面不同材料性质的间断特 性。节理允许发生滑动或分离，因此可以用来模拟岩体中的断层、 节理或摩擦边界。2、有五种计算模式?静力模式。这是FLAC-3D默认模式，通过动态松弛方法得 到静态解。?动力模式。用户可以直接输人加速度、速度或应力波作为系 统的边界条件或初始条件，边界可以固定边界和自由边界。 动力计算可以与渗流问题相耦合。?蠕变模式。有五种蠕变本构模型可供选择以模拟材料的应力 -应变-时间关系：Maxwell模型、双指数模型、参考蠕变模 型、粘塑性模型、脆延模型。?渗流模式。可以模拟地下水流、孔隙水压力耗散以及可变形孔隙介质与其间的粘性流体的耦合。渗流服从各向同性达西文案标准定律，流体和孔隙介质均被看作可变形体。考虑非稳定流，将稳定流看作是非稳定流的特例。边界条件可以是固定孔隙 压力或恒定流，可以模拟水源或深井。渗流计算可以与静力、 动力或温度计算耦合，也可以单独计算。?温度模式。可以模拟材料中的瞬态热传导以及温度应力。温度计算可以与静力、动力或渗流计算耦合，也可单独计算。3、可以模拟多种结构形式?对于通常的岩体