下载文档原格式
岩土工程软件FLAC3D的基本知识介绍岩土工程结构的数值解是建立在满足基本方程（平衡方程、几何方程、本构方程）和边界条件下推导的。
由于基本方程和边界条件多以微分方程的形式出现，因此，将基本方程近假发改用差分方程（代数方程）表示，把求解微分方程的问题改换成求解代数方程的问题，这就是所谓的差分法。
差分法由来已久，但差分法需要求解高阶代数方程组，只有在计算机的出现，才使该法得以实施和发展。
一、FLAC3D简介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++语言编写的。
1.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有以下几个优点：1 对模拟塑性破坏和塑性流动采用的是“混合离散法“。
1.基本介绍1.1.概述FLAC 3D是一个三维显式有限差分程序，主要应用于工程力学计算。
程序基于二维FLAC程序中已经建好的数值方程式。
FLAC 3D将FLAC的分析能力拓展到三维，用于模拟三维土体、岩体或其他材料的力学特性，尤其是达到屈服极限的塑性流变特性。
用户通过调整多面体单元的三维网格结构，来拟合要被建模的物体的实际形状。
每个单元体根据既定的线形/非线性的应力/应变规律对相应施加的力和边界约束条件作出响应。
并且当材料发生屈服流动后，网格也能够适应变形和移动（大变形模式）。
FLAC 3D采用的显式拉格朗日算法和混合——离散分区技术能够确保材料塑性坍塌破坏和流动过程的精确模拟。
由于无须形成刚度矩阵，因此采用较小的计算资源，就能够求解大范围的三维（岩土工程）计算问题。
通过自动惯性缩放及自动阻尼，显式方程式的缺点（小时间步限制和阻尼问题）已经被克服，且不会影响到物体的原有破坏行为。
FLAC 3D为三维岩土工程问题的解决提供了一个理想的分析工具。
FLAC 3D被设计，专门为了在装有Windows98及更高的版本的操作系统的IBM兼容的微型计算机上操作。
在岩土工程方面，实际的三维模型计算可以在合理的时间内被完成。
例如，创建一个包含大约140000个单元体的模型需要128M的内存。
对于一个有10000个单元体的摩尔——库伦模型，在2.4GHz的奔腾IV微型计算机上，完成5000个计算步需要大概18分钟。
对于显式计算求解，到达平衡状态的所需求解计算步数不定，但无论什么类型的模型，这个值都大概会在3000-5000步之内。
随着浮点数计算速度的提高，以及以低代价安装附加内存的能力，用FLAC 3D解决更大的三维问题成了可能。
FLAC 3D既可以通过命令行驱动，也可以通过图案菜单驱动。
默认的命令驱动模式和Itasca其他的软件产品是一样的。
你会发现其中大部分命令都是一样的。
在FLAC 3D中，菜单驱动的图形用户界面可用于绘图，显示工作。
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嵌语言FISHFLAC-3D具有强大嵌语言FISH, 使得用户可以定义新的变量或函数, 以适应用户的特殊需要, 例如, 利用HSH做以下事情:(l) 用户可以自定义材料的空间分布规律, 如非线性分布等.(2) 用户可以定义变量, 追踪其变化规律并绘图表示或打印输出.(3) 用户可以自己设计FLAC-3D部没有的单元形态.(4) 在数值试验中可以进行伺服控制.(5) 用户可以指定特殊的边界条件.(6) 自动进行参数分析(7) 利用FLAC-3D部定义的Fish变量或函数, 用户可以获得计算过程中节点、单元参数, 如坐标、位移、速度、材料参数、应力、应变、不平衡力等.FLAC-3D前后处理功能FLAC-3D具有强大的自动三维网格生成器, 部定义了多种单元形态, 用户还可以利用FISH自定义单元形态, 通过组合基本单元, 可以生成非常复杂的三维网格, 比如交叉隧洞等.在计算过程中的任何时刻用户都可以用高分辨率的彩色或灰度图或数据文件输出结果, 以对结果进行实时分析, 图形可以表示网格、结构以及有关变量的等值线图、矢量图、曲线图等, 可以给出计算域的任意截面上的变量图或等直线图, 计算域可以旋转以从不同的角度观测计算结果.FLAC3D计算分析一般步骤与大多数程序采用数据输入方式不同, FLAC采用的是命令驱动方式. 命令字控制着程序的运行. 在必要时, 尤其是绘图, 还可以启动FLAc用户交互式图形界面. 为了建立FLAC计算模型, 必须进行以下三个方面的工作:(1) 有限差分网格(2) 本构特性与材料性质(3) 边界条件与初始条件完成上述工作后, 可以获得模型的初始平衡状态, 也就是模拟开挖前的原岩应力状态. 然后, 进行工程开挖或改变边界条件来进行工程的响应分析, 类似于FLAC的显式有限差分程序的问题求解. 与传统的隐式求解程序不同, FLAC采用一种显式的时间步来求解代数方程. 进行一系列计算步后达到问题的解.在FLAC中, 达到问题所需的计算步能够通过程序或用户加以控制, 但是, 用户必须确定计算步是否已经达到问题的最终的解.后处理(一) 用tecplot绘制曲线(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例子的运行时间是区域数的4/3倍. 这个规则适用于平衡条件下的弹性问题. 对于塑性问题, 运行时间会有点改变, 但是不会很大, 但是如果发生塑性流动, 这个时间将会大的多. 对一个具体模型检查自己机子的计算速度很重要. 一个简单的方法就是运行基准测试. 然后基于区域数的改变, 用这个速度评估具体模型的计算速度.(2) 影响运行时间的因素:FLAC3D有时会需要较长时间才可以收敛主要发生在下列情况下:(a)材料本身刚度变异或材料与结构及接触面之间的刚度差异很大.(b)划分的区域尺寸相差很大. 这些尺寸差异越大编码就越无效. 在做详细分析前应该研究刚度差异的影响. 例如, 一个荷载作用下的刚性板, 可以用一系列顶点固定的网格代替, 并施以等速度. (记住FIX命令确定速度, 而不是位移. )地下水的出现将使体积模量发生明显的增加(流体-固体相互作用).(3) 考虑网格划分的密度:FLAC3D使用常应变单元. 如果应力/应变曲线倾斜度比较高, 那么你将需要许多区域来代表多变的分区. 通过运行划分密度不同的同一个问题来检查影响. FLAC3D应用常应变区域, 因为当用多的少节点单元与用比较少的多节点单元模拟塑性流动时相比更准确.应尽可能保持网格, 尤其是重要区域网格的统一. 避免长细比大于5:1的细长单元, 并避免单元尺寸跳跃式变化(即应使用平滑的网格). 应用GENERATE命令中的比率关键词, 使细划分区域平滑过渡到粗划分区域.(4) 自动发现平衡状态:默认情况下, 当执行SOLVE 命令时, 系统将自动发现力的平衡. 当模型中所有网格顶点中所有力的平均量级与其中最大的不平衡力的量级的比率小于1*10时, 认为达到了平衡状态. 注意一个网格顶点的力由力(例如, 由于重力)和外力(例如, 由于所加的应力边界条件)共同引起. 因为比率是没有尺寸的, 所以对于有不同的单元体系的模型, 在大多数情况下, 不平衡力和所加力比率的限制给静力平衡提供了一个精确的限制.同时还提供了其他的比率限制;可以用SET ratio 命令施加. 如果默认的比率限制不能为静力平衡提供一个足够精确的限制, 那么应考虑可供选择的比率限制. 默认的比率限制同样可用于热分析和流体分析的稳定状态求解. 对于热分析,是对不平衡热流量和所加的热流量量级进行评估, 而不是力. 对于流体分析,对不平衡流度和所加流度量级进行评估.(5) 考虑选择阻尼:对于静力分析, 默认的阻尼是局部阻尼, 对于消除大多数网格顶点的速度分量周期性为零时的动能很有效. 这是因为质量的调节过程依赖于速度的改变. 局部阻尼对于求解静力平衡是一个非常有效的计算法则且不会引入错误的阻尼力(见Cundall 1987).如果在求解最后状态, 重要区域的网格海域的速度分量不为零, 那么说明默认的阻尼对于达到平衡状态是不够的. 有另外一种形式的阻尼, 叫组合阻尼, 相比局部阻尼可以使稳定状态达到更好的收敛, 这时网格将发生明显的刚性移动. 例如, 求解轴向荷载作用下桩的承载力或模拟蠕变时都可能发生. 使用SETmechanical damp combined命令来调用组合阻尼. 组合阻尼对于减小动能方面不如局部阻尼有效, 所以应注意使系统的动力激发最小化. 可以用SETmechanical damp local命令转换到默认阻尼.(6) 检查模型反应:FLAC3D 显示了一个相试的物理系统是怎样变化的. 做一个简单的试验证明你在做你认为你在做的事情. 例如, 如果荷载和实体在几何尺寸上都是对称的, 当然反应也是对称的. 改变了模型以后, 执行几个时步(假如, 5或10步), 证明初始反应是正确的, 并且发生的位置是正确的. 对应力或位移的期望值做一个估计, 与FLAC3D 的输出结果作比较.如果你对模型施加了一个猛烈的冲击, 你将会得到猛烈的反应. 如果你对模型作了一些看起来不合理的事情, 你一定要等待奇怪的结果. 如果在分析的一个给定阶段, 得到了意外值, 那么回顾到这个阶段所用的时步.在进行模拟前很关键的是检查输出结果. 例如, 除了一个角点速度很大外, 一切都很合理, 那么在你理解原因前不要继续下去. 这种情况下, 你可能没有给定适当的网格边界.(7) 初始化变量:在模拟基坑开挖过程时, 在达到目的前通常要初始化网格顶点位移. 因为计算次序法则不要求位移, 所以可以初始化位移, 这只是由网格顶点的速度决定, 并有益于用户初始化速度却是一件难事. 如果设定网格顶点的速度为一常数, 那么这些点在设置否则前保持不变. 所以, 不要为了清除这些网格的速度而简单的初始化它们为零. . . 这将影响模拟结果. 然而, 有时设定速度为零是有用的(例如, 消除所有的动能).(8) 最小化静力分析的瞬时效应:对于连续性静力分析, 经过许多阶段逐步接近结果是很重要的. . . 即, 当问题条件突然改变时, 通过最小化瞬时波的影响, 使结果更加“静力”. 使FLAC3D 解决办法更加静态的方法有两种.(a) 当突然发生一个变化时(例如, 通过使区域值为零模拟开挖), 设定强度性能为很高的值以得到静力平衡. 然后为了确保不平衡力很低, 设定性能为真实值, 再计算, 这样, 由瞬时现象引起的失败就不会发生了.(b) 当移动材料时, 用FISH 函数或表格记录来逐步减少荷载.(9) 改变模型材料:FLAC3D 对一个模拟中所用的材料数没有限制. 这个准则已经尺寸化, 允许用户在自己所用版本的FLAC3D中最大尺寸网格的每个区域(假如设定的)使用不同的材料.(10) 运行在现场原位应力和重力作用下的问题:有很多问题在建模时需要考虑现场原位应力和重力的作用. 这种问题的一个例子是深层矿业开挖:回填. 此时大多数岩石受很高的原位应力区的影响(即, 自重应力由于网孔尺寸的限制可以忽略不计), 但是回填桩的放置使自重应力发展导致岩石在荷载作用下可能坍塌. 在这些模拟中要注意的重点(因为任何一种模拟都有重力的作用)是网格的至少三个点在空间上应固定. . . 否则, 整个网格在重力作用下将转动. 如果你曾经注意到整个网格在重力加速度矢量方向发生转动, 那么你可能忘记在空间上固定网格了.FLAC3D主要适明模拟计算地质材料和岩土上程的力学行为。
FLAC 3D基础知识介绍一、概述FLAC（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。
目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。
1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。
因此，大大发护展了计算规模。
FLAC3D是一个三维有限差分程序，目前已发展到V3.0版本。
FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。
因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。
FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。
FLAC3D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。
由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。
三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。
三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。
pfc基础知识FLAC 3D基础知识介绍一、概述FLAC（Fast Lagrangian Analysis of Continua）由美国Itasca 公司开发的。
目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。
1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。
因此，大大发护展了计算规模。
FLAC3D是一个三维有限差分程序，目前已发展到V3.0版本。
FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。
因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。
FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。
FLAC3D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。
由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。
三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。
三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。
《FLACFLAC3D基础与工程实例》阅读札记目录一、FLACFLAC3D软件概述 (2)1. 软件背景与简介 (3)1.1 FLACFLAC3D的发展历程 (4)1.2 软件的应用领域及特点 (5)2. 软件安装与运行环境 (6)2.1 系统要求 (7)2.2 安装步骤 (8)2.3 运行环境配置 (10)二、FLACFLAC3D基础知识 (11)1. 基本概念与术语 (13)1.1 有限元分析原理 (14)1.2 离散元法简介 (14)1.3 FLACFLAC3D中的相关术语解释 (15)2. 软件操作界面及功能模块 (17)2.1 操作界面介绍 (18)2.2 主要功能模块说明 (20)2.3 菜单功能详解 (20)三、工程实例分析 (22)1. 地质工程实例 (23)1.1 工程背景及问题定义 (25)1.2 模型建立与参数设置 (26)1.3 结果分析与讨论 (27)2. 土木工程实例 (29)2.1 工程概况与建模目的 (30)2.2 建模过程及计算步骤 (31)2.3 结果展示与工程应用 (32)四、FLACFLAC3D应用技巧与注意事项 (33)1. 建模技巧与优化方法 (34)1.1 建模策略及优化思路 (35)1.2 网格划分与模型简化技巧 (36)1.3 参数设置与模型验证方法 (38)2. 数据分析与处理方法 (40)2.1 数据采集与整理方法 (41)2.2 结果分析与图表展示技巧 (42)一、FLACFLAC3D软件概述3D是一种广泛使用的岩土力学与有限元分析软件。
它是一套专门用来分析连续介质中的物理力学现象的强大工具，主要应用于土木、矿山、隧道等领域，能针对各种复杂的工程问题进行数值建模和模拟分析。
3D以其高效、灵活的数值分析能力，为工程师提供了强大的技术支持。
其主要特点包括：多功能：3D能够模拟多种物理过程，包括应力分析、稳定性分析、流体流动分析等，适用于多种工程场景。
FLAC 3D基础知识介绍一、概述FLAC（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。
目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。
1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。
因此，大大发护展了计算规模。
FLAC3D是一个三维有限差分程序，目前已发展到V3.0版本。
FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。
因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。
FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。
FLAC3D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。
由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。
三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。
三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。
F L A C D基础知识介绍 Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】FLAC 3D基础知识介绍一、概述FLAC（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。
目前，FLAC有二维和三维计算程序两个版本，二维计算程序以前的为DOS版本，版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。
1995年，FLAC2D已升级为的版本，其程序能够使用护展内存。
因此，大大发护展了计算规模。
FLAC3D是一个三维有限差分程序，目前已发展到版本。
FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。
因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。
FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。
FLAC3D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。
由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。
三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。
FLAC3D软件原理及特点FLAC3D[1]是三维岩土力学有限差分计算机程序，是国际通用的岩土工程专业分析软件。
FLAC 代表连续介质快速拉格朗日分析，是由国际著名学者、英国皇家工程院院士、离散元法的发明人彼得库德尔Peter Cundall 博士在70年代中期开始研究开发的面向土木建筑、采矿、交通、水利、地质、核废料处理、石油及环境工程的通用软件系统，是美国艾塔斯Itasca国际咨询集团公司的软件核心产品最知名的软件系统之一。
自20世纪90年代中页，中国开始引进FLAC及FLAC3D等Itasca系列软件，许多工业部门都在应用FLAC系统进行工程设计、计算及科学研究；今天，FLAC已经成为我国岩土力学与工程界发展最快、最具有影响的数值分析软件系统。
FLAC程序将单元之间的不平衡力重新分配各节点之上，再进行下一步的迭代运算，直到不平衡力足够小或者各节点的唯一区域平衡为止，如下图1所示：图2 迭代求解过程图3 FLAC3D的特点3.1应用范围广泛FLAC3D是帮助土木、交通、采矿、水利工程师进行分析、测试及设计的连续介质程序。
由于其分析能力并不局限于某一类特殊问题或分析类型，FLAC3D得到了广泛的应用。
FLAC3D的设计思想是针对任何需要连续介质力学分析的岩土工程项目。
在分析岩土工程问题时，边界条件是最重要的考虑因素，而FLAC3D有多种边界条件，并且边界方位可以任意变化，边界条件可以是速度边界、应力边界，单元内部可以给定初始应力，节点可以给定初始位移、速度等，还可以给定地下水位以计算有效应力、所有给定量都可以具有空间梯度分布。
功能强大：FLAC3D是一个利用显示有限差分法为岩土工程提供精确有效分析的工具，可以解决诸多有限元程序难以模拟的复杂的工程问题。
另外，FLAC3D具有强大的内嵌语言FISH，使得用户可以定义新的变量或函数，以适应用户的特殊需要。
FISH可以做如下事情：（1）用户可以自定义材料的分布变化情况；（2）用户可以定义变量，追踪其变化规律并绘图表示；（3）用户可以自己设计FLAC3D内部所没有的单元形态；（4）在数值试验中可以检测控制；（5）用户可以指定特殊的边界条件；（6）自动进行参数分析；（7）利用FLAC3D内部定义的FISH变量或者函数，用户可以获得计算过程中的节点、单元参数，如坐标、位移、速度、材料参数、应力、应变和不平衡力等。
FLAC 3D学习笔记之网格生成器概述（一）FLAC生成网格一、网格生成器概述1. Generate命令可以调出网格生成器。
1) Generate zone：生成基本形状的网格。
gen zone可以调用的基本形状的网格关键词keyword 简写含义方向brick b 块状型(Brick-shaped mesh)； 3dbrick db 变形块状型(Degenerate Brick mesh) 3wedge w 楔型(Wedge-shaped mesh 3pyramid p 角锥型(Pyramid-shaped mesh) 3tetra tet 四面体(Tetrahedron-shaped mesh) 3cylinder c 圆柱型(Cylinder-shaped mesh) 3redbrick redb 放射状块状型(Radial graded mesh around Brick) 4radtunnel radt 放射状隧道型(Radial graded mesh around parallelepiped-shaped Tunnel) 4 radcylinder radc 放射状圆柱型(Radial graded mesh around Cylindrical-shaped tunnel) 4 cshell cs 圆柱型壳状物(Cylindrical Shell mesh) 4cylint cylint 圆柱交叉型(Cylinder Intersection) 5tunint t 隧道交叉型(Tunnel Intersection) 5Gen zon 生成基本形状网格的常用关键词关键词keyword 含义关键词keyword 含义dimension(dim) 内部区域的尺寸 p0-p16 各种形状网格的参考点edge 网格的边长 ratio(rat) 单元尺寸大小的比率fill 用网格填充的内部区域 size 网格在每个方向上的单元数目gen zon命令如下：generate zone <keyword> &p0 ( x0,y0,z0 ) p1 ( x1,y1,z1 ) p2 ( x2,y2,z2 ) p3 ( x3,y3,z3 ) &p4 ( x4,y4,z4 ) p5 ( x5,y5,z5 ) p6 ( x6,y6,z6 ) p7 ( x7,y7,z7 ) &p8 ( x8,y8,z8 ) p9 ( x9,y9,z9 ) p10 ( x10,y10,z10 ) p11 ( x11,y11,z11 ) &p12 ( x12,y12,z12 ) p13 ( x13,y13,z13 ) p14 ( x14,y14,z14 ) &p15 ( x15,y15,z15 ) p16 ( x16,y16,z16 )size n1<n2><n3><n4><n5> &ratio r1<r2><r3><r4><r5> &<keyword>其中表示为网格的型式;p0~p16为参考点,依不同型式之网格有不同数量之参考点,其中( xi,yi,zi )为各个参考点之座标; n1 <n2><n3><n4><n5>表对网格切分几个元素并针对不同区域切分,切分区域主要依不同型式之网格而定; r1<r2><r3><r4><r5>表对网格切分之元素依比例缩放。
有的行可省略，计算机按照默认的设定。
例1：gen zone radcyl size 5 10 6 12 dim 1 1 1 1 fillgroup tunnel range cyl end1 0 0 0 end2 0 10 0 rad 1plot add axe redplo blo gro见下图1图12) 连接简单形状得到网格：a) 通过设置不同的参考点参数，翻转基本网格，将不同基本网格连接例2：联合使用radcylinder和radtunnel建立马蹄形隧道模型见图2、3；radcyl建立圆柱体隧道周围的放射状网格gen zone radcyl size 5 10 6 12 rat 1 1 1 1.2 p0 0,0,0 p1 100,0,0 p2 0,200,0 p3 0,0,100 ；radtu 建立平行六面体隧道周边的放射状网格；radtu建立的网格通过制定不同的p1,p2,p3翻转了90°，从而和用radcyl建立的网格相连接gen zone radtun size 5 10 5 12 rat 1 1 1 1.2 p0 0,0,0 p1 0,0,-100 p2 0,200,0 p3 100,0,0图2图3b) gen zone copy 和gen zone reflect复制和镜像见图4、图5gen zone radcyl size 5 10 6 12 rat 1 1 1 1.2 p0 0,0,0 p1 100,0,0 p2 0,200,0 p3 0,0,100 gen zone radtun size 5 10 5 12 rat 1 1 1 1.2 p0 0,0,0 p1 0,0,-100 p2 0,200,0 p3 100,0,0 gen zone reflect dip 90 dd 270 origin 0,0,0图4dip是对称参考面与xy平面的夹角dd 是对称参考面的外法线在xy平面上投影的射线与y轴正方向的夹角或者，从y轴正向顺时针转向称参考面的外法线在xy平面上投影的射线的夹角。
2. Attach 命令attach 命令可以连接单元大小不同的基本网格，但是不同的单元大小应该是成整数倍的（2，3，4倍）。
例3，见图5、6。
gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,2gen zone brick size 8 8 4 p0 0,0,2 p1 4,0,2 p2 0,4,2 p3 0,0,4;attach命令attach face range z 1.9 2.1model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1；施加荷载的范围是在Z=4的平面上，0<x<2,0<y<2的小范围之内apply szz -1e6 range z 3.9,4.1 x 0,2 y 0,2；绘坐标轴plot add axe red；绘z方向位移云图plot con zdis；监测不平衡力hist unbalsolvesave attach.sav图5图6二、网格的组合来确定和实情形际类似的模型考虑：1，研究区域的单元密度；2，边界位置对模型计算影响gen point 生成的参考点用于建模期间确定原有网格连接范围是否合适。
确保连接尺寸的单元数量相同，单元尺寸变化相同比率一致。
例4：仰拱高程相同的两个隧道见图7;生成仰拱高程相同的两个隧道的网格;主隧道gen zone radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 &size 4 10 6 4 dim 4 4 4 4 rat 1 1 1 1 fillgen zone reflect dip 90 dd 90 ori 15 0 0gen zone reflect dip 0 ori 0 0 0;辅助隧道;radcyl生成辅助隧道的网格时辅助隧道半径的最低点point1必须进行调整;注意radcyl生成的网格形状只是四分之一的辅助隧道，其半径最低点，对应的是p8 p10;旋转网格90°后得到模型中辅助隧道右下四分之一部分，其p9 p11就是上面部分的p8 p10;x=2.969848和z=-0.575736是怎末得出的？答：而x/3=7/（50的平方根）;point1半径=3，辅助隧道和主隧道圆心距=7，辅助隧道圆心(0,0,-1),主隧道圆心(0，0，0)，辅助隧道的圆心距主隧道和辅助隧道分界线的中点：50的平方根gen point id 1 (2.969848,0.0,-0.575736)gen point id 2 (2.969848,50.0,-0.575736)gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 &p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 &size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1gen zon radcyl p0 0 0 -1 p1 0 0 -8 p2 0 50 -1 p3 7 0 0 p4 0 50 -8 &p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 &size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1;衬砌结构壳单元如何plo这种结构单元呢？;这里的衬砌和网格紧密相连，在网格上建立可以滑动或者分离的衬砌则须使用SEL liner 命令sel shell range cyl end1 0 0 -1 end2 0 50 -1 rad 3图7建立边界网格，见图8-13。
;outer boundary；用radtun建立边界网格如图8gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 &p5 15 50 50 p6 50 0 50 p7 50 50 50 &p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 &size 6 10 3 10 rat 1 1 1 1.1；brick建立边界网格如图9gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 &p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 &size 3 10 10 rat 1 1 1.1；做镜像见图10gen zon reflect dip 0 ori 0 0 0 range x 0 23 y 0 50 z 8 50；做镜像见图11gen zon reflect dip 0 ori 0 0 0 range x 23 50 y 0 50 z 0 50；最终成图见图12，组的显示见图13group service range cyl end1 0 0 -1 end2 0 50 -1 rad 3group main range cyl end1 15 0 0 end2 15 50 0 rad 4左图8 右图9图10 图11图12 图13群组建立方式:group <name> range cylinder end1 ( x1,y1,z1 ) end2 ( x2,y2,z2 ) r <rad>group <name> range y y1 y2或者group <name> range y y1,y2 x x1,x2 z z1, z2(中间的,要加就都加)群组建立方式中, <name>为设定群组之群组名称,cylinder表示以圆柱方式群组,end1 ( x1,y1,z1 )为起点座标 end2 ( x2,y2,z2 ) 为终点座标, <rad>为圆柱半径;下方为针对某一轴向群组,y1为起点座标 y2为终点座标.。
