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种材料 模型:(1)

3、开挖模型 null(2) 3 个弹性模型 (各向同性 , 横观各向同性和正交各向同性弹性模型 )6 个塑性模型(Drucker-Prager 模型、Morh-Coulomb模型、应变硬化/软化 模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模 型).Flac3D 网格中的每个区域可以给以不同的材料模型 , 并且还允许指定材料参数 的统计分布和变化梯度 . 还包含了节理单元 , 也称为界面单元 , 能够模拟两种 或多种材料界面不同材料性质的间断特性 . 节理允许发生滑动或分离 , 因此可 以用来模拟岩体中的断层、节理或摩擦边界 .FLAC3D中的网格生成器gen,通过匹配

4、、连接由网格生成器生成局部网格,能够 方便地生成所需要的三维结构网格 . 还可以自动产生交岔结构网格 (比如说相交 的巷道), 三维网格由整体坐标系 x, y, z 系统所确定, 这就提供了比较灵活的 产生和定义三维空间参数 .五种计算模式(1) 静力模式：这是FLAC-3D默认模式，通过动态松弛方法得静态解(2) 动力模式 : 用户可以直接输人加速度、速度或应力波作为系统的边界条件或 初始条件 , 边界可以固定边界和自由边界 . 动力计算可以与渗流问题相藕合 .(3) 蠕变模式 ： 有五种蠕变本构模型可供选择以模拟材料的应力-应变-时间关系 :Maxwell 模型、双指数模型、参考蠕变模型、

5、粘塑性模型、脆盐模型.(4) 渗流模式: 可以模拟地下水流、孔隙压力耗散以及可变形孔隙介质与其间的 粘性流体的耦合 . 渗流服从各向同性达西定律 , 流体和孔隙介质均被看作可变 形体 . 考虑非稳定流 , 将稳定流看作是非稳定流的特例 . 边界条件可以是固定 孔隙压力或恒定流 , 可以模拟水源或深井 . 渗流计算可以与静力、 动力或温度计 算耦合 , 也可以单独计算 .(5) 温度模式: 可以模拟材料中的瞬态热传导以及温度应力 . 温度计算可以与静 力、动力或渗流计算藕合 , 也可单独计算 .模拟多种结构形式(l) 对于通常的岩体、土体或其他材料实体 , 用八节点六面体单元模拟 .(2) FI

6、AC-3D 包含有四种结构单元 : 梁单元、锚单元、桩单元、壳单元 . 可用来 模拟岩土工程中的人工结构如支护、衬砌、锚索、岩栓、土工织物、摩擦桩、板 桩等 .(3) FLAC-3D的网格中可以有界面,这种界面将计算网格分割为若干部分,界面 两边的网格可以分离 , 也可以发生滑动 , 因此, 界面可以模拟节理、 断层或虚拟 的物理边界 .有多种边界条件边界方位可以任意变化 , 边界条件可以是速度边界、 应力边界, 单元部可以给定 初始应力 , 节点可以给定初始位移、 速度等 , 还可以给定地下水位以计算有效应 力、所有给定量都可以具有空间梯度分布 .FLAC-3D嵌语言 FISHFLAC-3D

7、具有强大嵌语言FISH,使得用户可以定义新的变量或函数,以适应用 户的特殊需要，例如，利用HSHf故以下事情：(1) 用户可以自定义材料的空间分布规律 , 如非线性分布等 .(2) 用户可以定义变量 , 追踪其变化规律并绘图表示或打印输出 . 用户可以自己设计FLAC-3D部没有的单元形态.(4) 在数值试验中可以进行伺服控制 .(5) 用户可以指定特殊的边界条件 .(6) 自动进行参数分析利用FLAC-3D部定义的Fish变量或函数，用户可以获得计算过程中节点、 单元参数 , 如坐标、位移、速度、材料参数、应力、应变、不平衡力等 .FLAC-3D前后处理功能FLAC-3D具有强大的自动三维网

8、格生成器，部定义了多种单元形态，用户还可 以利用FISH自定义单元形态，通过组合基本单元，可以生成非常复杂的三维网 格, 比如交叉隧洞等 .在计算过程中的任何时刻用户都可以用高分辨率的彩色或灰度图或数据文件输 出结果 , 以对结果进行实时分析 , 图形可以表示网格、 结构以及有关变量的等值 线图、矢量图、曲线图等 , 可以给出计算域的任意截面上的变量图或等直线图 计算域可以旋转以从不同的角度观测计算结果 .FLAC3D计算分析一般步骤与大多数程序采用数据输入方式不同 , FLAC 采用的是命令驱动方式 . 命令字控 制着程序的运行 . 在必要时,尤其是绘图,还可以启动FLAc用户交互式图形界面

9、为了建立FLAC计算模型, 必须进行以下三个方面的工作 :(1) 有限差分网格(2) 本构特性与材料性质(3) 边界条件与初始条件 完成上述工作后 , 可以获得模型的初始平衡状态 , 也就是模拟开挖前的原岩应 力状态 . 然后, 进行工程开挖或改变边界条件来进行工程的响应分析 , 类似于FLAC的显式有限差分程序的问题求解.与传统的隐式求解程序不同,FLAC采用 一种显式的时间步来求解代数方程 . 进行一系列计算步后达到问题的解 .在FLAC中，达到问题所需的计算步能够通过程序或用户加以控制，但是，用户 必须确定计算步是否已经达到问题的最终的解 .后处理( 一) 用 tecplot 绘制曲线(

10、1) 第一主应力(2) xdisp 、 ydisp 、zdisp 、 disp ( 二 ) 用 excel 做曲线隧道(1)做地表沉降槽 (zdisp)(2)地表横向位移 (xdisp)(3)隧道中线竖向沉降曲线 (zdisp)(4)提取位移矢量图 ,(5)显示初期支护结构力(6)显示 state( 找塑性区 )基坑(1)做地表沉降槽 (zdisp)(2)提取位移矢量图 ,(3)显示初期支护结构力(4)显示 state( 找塑性区 )边坡(1)做安全系数和应变图模型最优化 用FLAC3D军决问题时,为了得到最有效的分析使模型最优化是很重要的(1) 检查模型运行时间 :一个FLAC3D例子的运行

11、时间是区域数的4/3倍.这个规则适用于平衡条件下的 弹性问题 . 对于塑性问题 , 运行时间会有点改变 , 但是不会很大 , 但是如果发 生塑性流动 , 这个时间将会大的多 . 对一个具体模型检查自己机子的计算速度 很重要 . 一个简单的方法就是运行基准测试 . 然后基于区域数的改变 , 用这个 速度评估具体模型的计算速度 .(2) 影响运行时间的因素 :FLAC3D有时会需要较长时间才可以收敛主要发生在下列情况下：(a) 材料本身刚度变异或材料与结构及接触面之间的刚度差异很大 .(b) 划分的区域尺寸相差很大 . 这些尺寸差异越大编码就越无效 . 在做详细分析 前应该研究刚度差异的影响 .

12、例如, 一个荷载作用下的刚性板 , 可以用一系列 顶点固定的网格代替，并施以等速度(记住FIX命令确定速度，而不是位移.) 地下水的出现将使体积模量发生明显的增加 (流体-固体相互作用 ).(3) 考虑网格划分的密度 ：FLAC3D使用常应变单元.如果应力/应变曲线倾斜度比较高,那么你将需要许 多区域来代表多变的分区 . 通过运行划分密度不同的同一个问题来检查影响 .FLAC3D应用常应变区域,因为当用多的少节点单元与用比较少的多节点单元模 拟塑性流动时相比更准确 .应尽可能保持网格 , 尤其是重要区域网格的统一 . 避免长细比大于 5：1 的细长 单元,并避免单元尺寸跳跃式变化(即应使用平滑

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

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

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

16、物理系统是怎样变化的.做一个简单的试验证明你在 做你认为你在做的事情 . 例如, 如果荷载和实体在几何尺寸上都是对称的 , 当 然反应也是对称的 . 改变了模型以后 , 执行几个时步 (假如, 5 或 10步), 证明 初始反应是正确的 , 并且发生的位置是正确的 . 对应力或位移的期望值做一个 估计,与FLAC3D的输出结果作比较.如果你对模型施加了一个猛烈的冲击 , 你将会得到猛烈的反应 . 如果你对模型 作了一些看起来不合理的事情 , 你一定要等待奇怪的结果 . 如果在分析的一个 给定阶段 , 得到了意外值 , 那么回顾到这个阶段所用的时步 . 在进行模拟前很关键的是检查输出结果 . 例

17、如, 除了一个角点速度很大外 , 一 切都很合理 , 那么在你理解原因前不要继续下去 . 这种情况下 , 你可能没有给 定适当的网格边界 .(7) 初始化变量 :在模拟基坑开挖过程时 , 在达到目的前通常要初始化网格顶点位移 . 因为计算 次序法则不要求位移 , 所以可以初始化位移 , 这只是由网格顶点的速度决定 , 并有益于用户初始化速度却是一件难事 . 如果设定网格顶点的速度为一常数 , 那么这些点在设置否则前保持不变 . 所以, 不要为了清除这些网格的速度而简 单的初始化它们为零 . . . 这将影响模拟结果 . 然而, 有时设定速度为零是有用 的(例如, 消除所有的动能 ).(8) 最

18、小化静力分析的瞬时效应 :对于连续性静力分析 , 经过许多阶段逐步接近结果是很重要的 . . . 即, 当问题 条件突然改变时 , 通过最小化瞬时波的影响 , 使结果更加“静力” . 使 FLAC3D 解决办法更加静态的方法有两种 .(a) 当突然发生一个变化时 ( 例如, 通过使区域值为零模拟开挖 ), 设定强度性 能为很高的值以得到静力平衡 . 然后为了确保不平衡力很低 , 设定性能为真实 值, 再计算, 这样, 由瞬时现象引起的失败就不会发生了 .(b) 当移动材料时 , 用 FISH 函数或表格记录来逐步减少荷载 .(9) 改变模型材料 :FLAC3D对一个模拟中所用的材料数没有限制.

19、这个准则已经尺寸化,允许用户 在自己所用版本的FLAC3D中最大尺寸网格的每个区域(假如设定的)使用不同的 材料.(10) 运行在现场原位应力和重力作用下的问题 : 有很多问题在建模时需要考虑现场原位应力和重力的作用 . 这种问题的一个例 子是深层矿业开挖 :回填. 此时大多数岩石受很高的原位应力区的影响 (即, 自 重应力由于网孔尺寸的限制可以忽略不计 ), 但是回填桩的放置使自重应力发展 导致岩石在荷载作用下可能坍塌 . 在这些模拟中要注意的重点 (因为任何一种模 拟都有重力的作用 )是网格的至少三个点在空间上应固定 . . . 否则, 整个网格 在重力作用下将转动 . 如果你曾经注意到整

20、个网格在重力加速度矢量方向发生 转动, 那么你可能忘记在空间上固定网格了 .FLAC3吐要适明模拟计算地质材料和岩土上程的力学行为。特别是材枓达到屈 服极限后产生的塑性流动、 材料通过单元和区域表示， 根据计算对象的形状构成 相应的网格。 每个单元在外载和边界约束条件下， 按照约定的线性或非线性应力应变关系产生力学响应。由于 FLAC3D?序主要是为岩上工程应用而开发的岩 石力学计算程序， 程序中包括了反映地质材料力学效应的特殊计算功能， 可计算 地质类材料的高度非线性 (包括应变硬化 / 软化) 、不可逆剪切破坏和压密、粘弹 ( 蠕变)、孔隙介质的应力渗流耦合、热力耦合以及动力学行为等；FL

21、AC3D程序设有多种本构模型：(1) 各向同性弹性材料模型；(2) 横观各向同性弹性材枓模型；(3) 莫尔库仑弹塑材料模型；(4) 应变软化 / 硬化塑性材科模型；(5) 双屈服塑性材料模型；(6) 遍布节理材料模型；(7) 空单元模型，可用来模拟地下硐室的开挖和煤层开采； 另外，程序设有界面单元，可以模拟断层、节理和摩擦边界的滑动、开和闭合行 为。支护结构如砌衬、锚杆、可缩性支架或板壳等与围岩的相互作用也可以在 FLAC3D中进行模拟。同时、用户可根据需要在 FLAC3冲创建自己的本构模型， 进行各种特殊修正和补充。FLAC3D1序建立在拉格朗日算法基础上，特别适合模拟大变形和扭曲。FLAC

22、3D采用显式算法来获得模型全部运动方程 (包括变量) 的时间步长解， 从而可以追 踪材料的渐进破坏和垮落， 这对研究采矿设计是非常重要的。 此外，程序允许输 入多种材枓类型， 亦可在计算过程中改变某个局部的材料参数， 增强了程序使甲 的灵活性。极大地方便了在计算上的处理。FLAC3D?序具有强大的后处理功能， 用户可以直接在屏幕上绘制或以文件形式创建和输出打印多种形式的图形。 用户 还可根据需要， 将若干个变量合并在同一幅图形中进行研究分析。 基于上述计算 功能与特点。本研究应用FLAC3D?序计算分析付村煤矿209工作面山上开采对 山下回采巷道围岩变形的影响。FLAC(Fast Lagran

23、gian Analysis of Continua)是一个利用显式有限差分方法求解的岩土、采矿工程师进行分析和设计的二维连续介质程序， 主要用来模拟土、 岩、或其他材料的非线性力学行为， 可以解决众多有限元程序难以模拟的复杂的 工程问题，例如大变形、大应变、非线性及非稳定系统(甚至大面积屈服/ 失稳或完全塌方)等问题。FLAC的基本功能和特征为：? 允许介质出现大应变和大变形；? Interface 单元可以模拟连续介质中的界面，并允许界面发生滑动和开 裂；? 显式计算方法， 能够为非稳定物理过程提供稳定解， 直观反映岩土体工程 中的破坏；? 地下水流动与力学计算完全耦合 (包括负孔隙水压，

24、非饱和流及相界面计 算)；? 采用结构加固单元模拟加固措施，例如衬砌、锚杆、桩基等；? 材料模型库(例如：弹性模型、莫尔库仑塑性模型、任意各向异性模型、 双屈服模型、粘性及应变软化模型)；? 预定义材料性质， 用户也可增加用户自己的材料性质设定并储存到数据库 中；? 一系列可选择模块，包括：热力学模块、流变模块、动力学模块、二相流 模块等，用户还可用C+建立自己的模型；? 边坡稳定系数计算满足边坡设计的要求；?用户可用部语言（FISH）增加自己定义的各种特性（如：新的本构模型， 新变量或新命令）；FLAC软件的优势：? 连续体大应变模拟? 界面单元用已代表不连续接触界面可能出现的完全不连续性质

25、的开和滑 动，因此可以模拟断层、节理和摩擦边界等? 显式求解模式可以获得不稳定物理过程的稳定解? 材料模型 :o “空（ null ）” 模型；o 三种弹性模型（各向同性、横观各向异性、和正交各向异性）；o 七种非线性模型 （ Drucker-Prager 、 Mohr-Coulomb、 应变硬化及 应变软化、 节理化、双线性应变硬化 / 软化节理化、双屈服、 修 正的 Cam-clay 模型）? 任何参数指标的连续变化或统计分布的模拟?外接口编程语言 （FISH）允许用户添加用户自定义功能? 方便的边界定义和初始条件定义方式? 可定义水位线 / 面进行有效应力计算? 地下水渗流计算以及完全的

26、应力场渗流场偶合计算 （含负孔隙压力、 非饱 和流、井）? 结构单元如隧道衬砌、桩、壳、梁锚杆、锚索、土工织物及其组合，可以 模拟不同的加固手段及其与围岩（土体）的相互作用? 自选模块包括 :o 热和热力学分析模块；o 流变计算模块；o 动力分析模块实现真时间历程的瞬时动力响应模拟；o用C+编写的用户自定义本构模块o开挖直立坡的喷射混凝土墙加土钉加固的模拟加（下）和不加（上）土工织物土坡的潜在破坏特征抗脅桩吏护码头的土FLACd（ Three Dimensional Fast Lagrangian Analysis ofContinua ）是美国 Itasca Consulting Goup

27、Inc开发的三维快速拉格朗日分析程序，是二维的有限差分程序FLAC2D勺扩展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分 析。该程序能较好地模拟地质材料在达到强度极限或屈服极限时发 生的破坏或塑性流动的力学行为，特别适用于分析渐进破坏和失稳 以及模拟大变形。FLAC3D调整三维网格中的多面体单元来拟合实 际的结构。单元材料可采用线性或非线性本构模型，在外力作用下， 当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模 式）。FLAC3D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。由于无须形成刚度矩阵，因此，基于较小存空间就能够求

28、解大围的三维问题。FLACD采用ANSI C+语言编写的。它包含10种弹塑性材料本构模 型，有静力、动力、蠕变、渗流、温度五种计算模式，各种模式间 可以互相藕合，可以模拟多种结构形式，如岩体、土体或其他材料 实体，梁、锚元、桩、壳以及人工结构如支护、衬砌、锚索、岩栓、 土工织物、摩擦桩、板桩、界面单元等，可以模拟复杂的岩土工程 或力学问题。flaCd的优点：1. 对模拟塑性破坏和塑性流动采用的是”混合离散法”。这种 方法比有限元法常采用的”离散集成法”更为准确、合理。2. 即使模拟的系统是静态的，仍采用了动态运动方程，这使得 FLAC3匪模拟物理上的不稳定过程不存在数值上的障碍。3. 采用了一

29、个”显式解”方案。因此，显式解方案对非线性的 应力-应变关系的求解所花费的时间，几互与线性本构关系 相同，而隐式求解方案将会花费较长的时间求解非线性问题。面且，它没有必要存储刚度矩阵，这就意味着，采用中等容量的存可以求解多单元结构； 模拟大变形问题几乎 并不比小变形问题多消耗更多的计算时间， 因为没有任何刚 度矩阵要被修改。槽形开挖PjFi A 匸吕 n ?.io-ioa Alphd 151 - View oilpile挣 File Ediit Settings Plottern WindowFLAC3D勺功能和特征: 1、包含10材料本构模型Flac3D 中为岩土工程问题的求解开发了特有的本

30、构模型，总共包 含了 10 种材料模型：? 开挖模型 null ；? 3 个弹性模型（各向同性，横观各向同性和正交各向同性弹 性模型）? 6 个塑性模型（ Drucker-Prager 模型、 Morh-Coulomb 模型、 应变硬化 / 软化模型、 遍布节理模型、 双线性应变硬化 / 软化 遍布节理模型和修正的cam粘土模型）。Flac3D 网格中的每个区域可以给以不同的材料模型，并且还允许 指定材料参数的统计分布和变化梯度。 还包含了节理单元，也称 为界面单元，能够模拟两种或多种材料界面不同材料性质的间断特 性。节理允许发生滑动或分离，因此可以用来模拟岩体中的断层、 节理或摩擦边界。2、

31、有五种计算模式?静力模式。这是FLAC-3D默认模式，通过动态松弛方法得到 静态解。? 动力模式。用户可以直接输人加速度、 速度或应力波作为系 统的边界条件或初始条件，边界可以固定边界和自由边界。 动力计算可以与渗流问题相耦合。? 蠕变模式。有五种蠕变本构模型可供选择以模拟材料的应力- 应变 - 时间关系 :Maxwell 模型、双指数模型、参考蠕变模 型、粘塑性模型、脆延模型。? 渗流模式。 可以模拟地下水流、 孔隙水压力耗散以及可变形 孔隙介质与其间的粘性流体的耦合。 渗流服从各向同性达西 定律，流体和孔隙介质均被看作可变形体。考虑非稳定流， 将稳定流看作是非稳定流的特例。 边界条件可以是固定孔隙 压力或恒定流，可以模拟水源或深井。 渗流计算可以与静力、 动力或温度计算耦合，也可以单独计算。? 温度模式。可以模拟材料中的瞬态热传导以及温度应力。 温