FLAC数值模拟介绍.docVIP

FLAC-3D(Three Dimensional Fast Lagrangian Analysis of Continua)是美国Itasca Consulting Goup lnc开发的三维快速拉格朗日分析程序, 该程序能较好地模拟地质材料在达到强度极限或屈服极限时, 发生的破坏或塑性流动的力学行为, 特别适用于分析渐进破坏和失稳以及模拟大变形. FLAC3D分析的使用领域 根据手册总结如下: (1) 承受荷载能力与变形分析: 用于边坡稳定和基础设计 (2) 渐进破坏与坍塌反演: 用于硬岩采矿和隧道设计 (3) 断层构造的影响研究: 用于采矿设计 (4) 施加于地质体锚索支护所提供的支护力研究: 岩锚和土钉的设计 (5) 排水和不排水加载条件下全饱和流体流动和孔隙压力扩散研究: 挡土墙结构的地下水流动, 和土体固结研究 (6) 粘性材料的蠕变特性: 用于碳酸钾盐矿设计 (7) 陡滑面地质结构的动态加载: 用于地震工程和矿山岩爆研究 (8) 爆炸荷载和振动的动态响应: 用于隧道开挖和采矿活动 (9) 结构的地震感应: 用于土坝设计 (10) 由于温度诱发荷载所导致的变形和结构的不稳定 (11) 大变形材料分析: 用于研究粮仓谷物流动和放矿的矿石流动 10种材料本构模型 Flac3D中为岩土工程问题的求解开发了特有的本构模型, 总共包含了10种材料模型: (1) 开挖模型null (2) 3个弹性模型(各向同性, 横观各向同性和正交各向同性弹性模型) (3) 6个塑性模型(Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模型). Flac3D网格中的每个区域可以给以不同的材料模型, 并且还允许指定材料参数的统计分布和变化梯度. 还包含了节理单元, 也称为界面单元, 能够模拟两种或多种材料界面不同材料性质的间断特性. 节理允许发生滑动或分离, 因此可以用来模拟岩体中的断层、节理或摩擦边界. FLAC3D中的网格生成器gen, 通过匹配、连接由网格生成器生成局部网格, 能够方便地生成所需要的三维结构网格. 还可以自动产生交岔结构网格(比如说相交的巷道), 三维网格由整体坐标系x, y, z系统所确定, 这就提供了比较灵活的产生和定义三维空间参数. 五种计算模式 (l) 静力模式:这是FLAC-3D默认模式, 通过动态松弛方法得静态解. (2) 动力模式:用户可以直接输人加速度、速度或应力波作为系统的边界条件或初始条件, 边界可以固定边界和自由边界. 动力计算可以与渗流问题相藕合. (3) 蠕变模式:有五种蠕变本构模型可供选择以模拟材料的应力-应变-时间关系:Maxwell模型、双指数模型、参考蠕变模型、粘塑性模型、脆盐模型. (4) 渗流模式:可以模拟地下水流、孔隙压力耗散以及可变形孔隙介质与其间的粘性流体的耦合. 渗流服从各向同性达西定律, 流体和孔隙介质均被看作可变形体. 考虑非稳定流, 将稳定流看作是非稳定流的特例. 边界条件可以是固定孔隙压力或恒定流, 可以模拟水源或深井. 渗流计算可以与静力、动力或温度计算耦合, 也可以单独计算. (5) 温度模式:可以模拟材料中的瞬态热传导以及温度应力. 温度计算可以与静力、动力或渗流计算藕合, 也可单独计算. 模拟多种结构形式 (l) 对于通常的岩体、土体或其他材料实体, 用八节点六面体单元模拟. (2) FIAC-3D包含有四种结构单元:梁单元、锚单元、桩单元、壳单元. 可用来模拟岩土工程中的人工结构如支护、衬砌、锚索、岩栓、土工织物、摩擦桩、板桩等. (3) FLAC-3D的网格中可以有界面, 这种界面将计算网格分割为若干部分, 界面两边的网格可以分离, 也可以发生滑动, 因此, 界面可以模拟节理、断层或虚拟的物理边界. 有多种边界条件 边界方位可以任意变化, 边界条件可以是速度边界、应力边界, 单元内部可以给定初始应力, 节点可以给定初始位移、速度等, 还可以给定地下水位以计算有效应力、所有给定量都可以具有空间梯度分布. FLAC-3D内嵌语言FISH FLAC-3D具有强大内嵌语言FISH, 使得用户可以定义新的变量或函数, 以适应用户的特殊需要, 例如, 利用HSH做以下事情:(l) 用户可以自定义材料的空间分布规律, 如非线性分布等. (2) 用户可以定义变量, 追踪其变化规律并绘图表示或打印输出. (3) 用户可以自己设计FLAC-3D内部没有的单元形态. (4) 在数值试验中可以进行伺服控制. (5) 用户可以指定特殊的边界条件. (6) 自动进行参数分析(7) 利用FLAC-3D内部定义的Fish变量或函数