FLAC 3D与FLAC常见问题的整理

Flac 3D 学习经验积累(by seapie)用FLAC3D 解决问题时，为了得到最有效的分析使模型最优化是很重要的。

通过资料收集和作者本人经验，对改进模型的运行提供了一些方法建议。

一、学习经验积累1、模型建立1.1 考虑网格划分的密度FLAC3D 使用常应变单元。

如果应力/应变曲线倾斜度比较高，那么你将需要许多区域来代表多变的分区。

通过运行划分密度不同的同一个问题来检查影响。

FLAC3D 应用常应变区域，因为当用多的少节点单元与用比较少的多节点单元模拟塑性流动时相比更准确。

（见理论卷第二章和实例卷中的塑性例子）。

应尽可能保持网格，尤其是重要区域网格的统一。

避免长细比大于5:1的细长单元，并避免单元尺寸跳跃式变化（即应使用平滑的网格）。

应用GENERATE 命令中的比率关键词，使细划分区域平滑过渡到粗划分区域。

2、本构模型与材料属性2.1 改变模型材料FLAC3D 对一个模拟中所用的材料数没有限制。

这个准则已经尺寸化，允许用户在自己所用版本的FLAC3D 中最大尺寸网格的每个区域（假如设定的）使用不同的材料。

3、边界条件3.1 初始应力场的加载岩土工程除常应力场外一般均应加重力场set grav 0 0 -10.0，并应加单元的密度ini density dens1(kg/m 3)。

若还考虑构造应力场时，两种情况：1 常构造应力，命令用apply ，如apply sxx -1.0e52 递度构造应力，命令用ini ，如ini sxx -80 grad 0 0 8.0，注意gz=8.0注意：在加载SZZ 或重力场时，程序自动按)1(μμ-的倍数关系加水平应力场，所以在通常的加载水平应力场时应考虑这部分的值。

3.2 边界约束初始化变量在模拟基坑开挖过程时，在达到目的前通常要初始化网格顶点位移。

因为计算次序法则不要求位移，所以可以初始化位移，这只是由网格顶点的速度决定，并有益于用户初始化速度却是一件难事。

FLAC若干问题的解答把结果中的信息输出到指定文件分为两种：1、把单元信息输出到指定文件set log onset log on finame.dat;可以输入文件路径，否则按当前目录处理print zoneset log off通过反复使用该命令，可以把不同信息输出到不同的文件。

2、把节点信息输出到指定文件set log onset log on finame.dat;可以输入文件路径，否则按当前目录处理print gpset log off请教一下print gp里的keyword该如何填。

比如想输出某点的zdis谢谢print gp dispprint gp positionlog 文件中有，然后处理，这些论坛里都有的关于flac3d—内置fish语言精讲FISH语言是FLAC3D程序的内置编程语言因为FLAC3D的最佳操作方式是命令流文件方式这一点与ANSYS很相似，而FISH就相当于ANSYS的APDL语言。

它包括循环、判断等结构。

如果你还用过其它高级语言，那么从形式上讲你也可以把它理解为子函数。

FISH语言的引入极大的方便了用户进行复杂的程序建模它不但可以嵌入命令流文件里工作而且还可以引用FLAC3D本身的任何命令所以说它实现了对FLAC3D的完全控制引入FISH语言的作法值得其它通用软件效仿。

以上仅是个人的使用体会请版主与管理员指教也供初学者入门之用。

谢谢！Fish函数增加了以下几种新的特性：1.增加了fish变量来获取结点、单元和界面变量2.fish提供了获取结构单元变量的途径3.休单元和面单元性质目前可以通过单元变量名z_prop(i_z string)和界面单元名i_prop(i_z string)分别加以识别4.fish函数可以获取单元应变和应力速率，还提供了全应变增量张量和应变速率张量5.提供了fish绘图子程序函数够生用户定义的图形内容6.fish函数已经增加了从文件读、写数据的能力FLAC3D是一个强大的软件，但是不得不承认的，它的界面没有ansys好用，商业运行不是很强的说也来凑个热闹，呵呵注意点：1.fish函数可以嵌套使用;2.以save命令保存模型时，fish函数和变量也同时保存；3.fish函数不支持缩写，这与flac3d命令不同，另外所有的fish 函数或变量不区分大小写，程序同意转化为大写进行编译，当然也可以通过执行set case-sensitivity on来区分大小写；4.变量或函数名不能以一个数字开头也不能是下列字符：. , * / ^ = > < # ( ) [ ] @ ; “ '5.如果用命令set safe on 指定了编译安全模式，则用户调用fish 函数时，函数名前必须加@；6.如果变量不曾赋值，则系统默认为零（整形），如果赋值，其类型由值的类型决定；7.fish函数的调用方法：.可以出项在其他fish函数的单独行中；.可以出现在其他fish函数的表达式中；.出现在flac3d的命令行中；.作为命令set,print,hist的参数。

岩土工程结构的数值解是建立在满足基本方程（平衡方程、几何方程、本构方程）和边界条件下推导的。

由于基本方程和边界条件多以微分方程的形式出现，因此，将基本方程近假发改用差分方程（代数方程）表示，把求解微分方程的问题改换成求解代数方程的问题，这就是所谓的差分法。

差分法由来已久，但差分法需要求解高阶代数方程组，只有在计算机的出现，才使该法得以实施和发展。

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 对模拟塑性破坏和塑性流动采用的是“混合离散法“。

岩土工程结构的数值解是建立在满足基本方程（平衡方程、几何方程、本构方程）和边界条件下推导的。

由于基本方程和边界条件多以微分方程的形式出现，因此，将基本方程近假发改用差分方程（代数方程）表示，把求解微分方程的问题改换成求解代数方程的问题，这就是所谓的差分法。

差分法由来已久，但差分法需要求解高阶代数方程组，只有在计算机的出现，才使该法得以实施和发展。

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 对模拟塑性破坏和塑性流动采用的是“混合离散法“。

FLAC3D常见命令与使用技巧1、FLAC3D常见命令：1.FLAC3D是有限元程序吗？答：不是！是有限差分法。

2.最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。

3.怎样看模型的样子？答：plo blo gro可以看到不同的group的颜色分布4.怎样看模型的边界情况？答：plo gpfix red5.怎样看模型的体力分布？答：plo fap red6.怎样看模型的云图？答：位移：plo con dis (xdis, ydis, zdis)应力：plo con sz (sy, sx,sxy, syz, sxz)7.怎样看模型的矢量图？答：plo dis (xdis, ydis, zdis)8.怎样看模型有多少单元、节点？答：pri info9.怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件10.怎样调用一个文件？答：File/call或者call命令10.如何施加面力？答：app nstress11.如何调整视图的大小、角度？答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键12.如何进行边界约束？答：fix x ran（约束的是速度，在初始情况下约束等效于位移约束）13.如何知道每个单元的ID？答：用鼠标双击单元的表面，可以知道单元的ID和坐标14.如何进行切片？答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片)15.如何保存计算结果？答：save +文件名.16.如何调用已保存的结果？答：rest +文件名；或者File / Restore17.如何暂停计算？答：Esc18.如何在程序中进行暂停，并可恢复计算？答：在命令中加入pause命令，用continue进行继续19.如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步20. Fish是什么东西？答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学？答：可以不用，需要的时候查Mannual获得需要的变量就可以了22.FLAC3D允许的命令文件格式有哪些？答：无所谓，只要是文本文件，什么后缀都可以23.如何调用一些可选模块？答：config dyn (fluid, creep, cppudm)24 .如何在圆柱体四周如何施加约束条件？答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2radius r2 not，其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用set log on命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

目录1．FLAC3D的固流耦合计算模式 1 2．FLAC3D固流耦合学习小结 53．关于流固耦合的问题 64．也谈采用FLAC3D对地下采矿的模拟 8 5．FLAC3D本构模型开发 86．FLAC3D自定义本构模型 117．数值计算中初始应力场的模拟 13 8．FLAC3D应变分析 139．FLAC3D的调参 1410．开采沉陷垂直剖面等值线的生成 15 11．FLAC3D的应变硬化软化模型 16 12．FLAC3D的塑性流动格式 17 13．FLAC3D的动画制作 1714．地下连续墙基坑开挖支护 1815．一个汇的小例子 2116．用3DEC生成岩体随机节理网络 2317．固结小算例 24FLAC3D的固流耦合计算模式英文原文 261．FLAC3D的固流耦合计算模式FLAC3D的计算模式中是否需要做孔压分析取决于是否采用config fluid命令。

1 无渗流模式（不使用config fluid）即使不使用命令config fluid，仍然可以在节点上施加孔压。

这种模式下，孔压将保持为常量。

如果采用塑性本构模型的话，材料的破坏将由有效应力状态来控制。

节点上的孔压分布可由initial pp命令或water table命令来设定。

如果采用water table命令，由程序自动计算水位线以下的静水孔压分布。

此时，必须施加流体密度（water density）和重力（set gravity）。

流体密度值和水位位置可以用命令print water显示。

如果水位线是由face关键字来定义的，则可用命令plot water命令显示水位。

这两种情况，单元的孔压都由节点孔压值平均求出，并在本构模型计算中用作有效应力。

这种计算模式下，体积力中不反映流体的出现：用户必须根据水位线以上或以下相应地指定干密度和湿密度。

使用命令print gp pp和print zone pp 可分别得到节点或单元孔压。

plot contour pp命令可绘出节点孔压云图。

有限差分法FLAC3D功能,优缺点分析FLAC/FLAC3D系列——岩⼟体⼯程⾼级连续介质⼒学分析软件通知：FLAC3D 4.0隆重推出，了解详细情况，点击此处FLAC（Fast Lagrangian Analysis of Continua）软件是由美国Itasca公司开发的。

⽬前，FLAC有⼆维和三维计算程序两个版本，⼆维计算程序V3.0以前的为DOS版本， 1995年，FLAC2D已升级为V3.3的版本，其程序能够使⽤护展内存，⾄今已发展到V5.0版本。

FLAC3D是⼀个三维有限差分程序，⽬前已发展到V4.0版本。

并且其推出的FLAC SLOPE有了WINDOWS界⾯。

FLAC（Fast Lagrangian Analysis of Continua）是⼀个利⽤显式有限差分⽅法求解的岩⼟、采矿⼯程师进⾏分析和设计的⼆维连续介质程序，主要⽤来模拟⼟、岩、或其他材料的⾮线性⼒学⾏为，可以解决众多有限元程序难以模拟的复杂的⼯程问题，例如⼤变形、⼤应变、⾮线性及⾮稳定系统（甚⾄⼤⾯积屈服/失稳或完全塌⽅）等问题。

FLAC的基本功能和特征为：●允许介质出现⼤应变和⼤变形；●Interface 单元可以模拟连续介质中的界⾯，并允许界⾯发⽣滑动和开裂；●显式计算⽅法，能够为⾮稳定物理过程提供稳定解，直观反映岩⼟体⼯程中的破坏；●地下⽔流动与⼒学计算完全耦合（包括负孔隙⽔压，⾮饱和流及相界⾯计算）；●采⽤结构加固单元模拟加固措施，例如衬砌、锚杆、桩基等；●材料模型库（例如：弹性模型、莫尔库仑塑性模型、任意各向异性模型、双屈服模型、粘性及应变软化模型）；●预定义材料性质，⽤户也可增加⽤户⾃⼰的材料性质设定并储存到数据库中；●⼀系列可选择模块，包括：热⼒学模块、流变模块、动⼒学模块、⼆相流模块等，⽤户还可⽤C++建⽴⾃⼰的模型；●边坡稳定系数计算满⾜边坡设计的要求；●⽤户可⽤内部语⾔（FISH）增加⾃⼰定义的各种特性（如：新的本构模型，新变量或新命令）；FLAC软件的优势：连续体⼤应变模拟界⾯单元⽤已代表不连续接触界⾯可能出现的完全不连续性质的张开和滑动，因此可以模拟断层、节理和摩擦边界等显式求解模式可以获得不稳定物理过程的稳定解材料模型:“空（null）”模型；三种弹性模型（各向同性、横观各向异性、和正交各向异性）；七种⾮线性模型（Drucker-Prager、Mohr-Coulomb、应变硬化及应变软化、节理化、双线性应变硬化/软化节理化、双屈服、修正的Cam-clay模型）任何参数指标的连续变化或统计分布的模拟外接⼝编程语⾔（FISH）允许⽤户添加⽤户⾃定义功能⽅便的边界定义和初始条件定义⽅式可定义⽔位线/⾯进⾏有效应⼒计算地下⽔渗流计算以及完全的应⼒场渗流场偶合计算（含负孔隙压⼒、⾮饱和流、井）结构单元如隧道衬砌、桩、壳、梁锚杆、锚索、⼟⼯织物及其组合，可以模拟不同的加固⼿段及其与围岩（⼟体）的相互作⽤⾃选模块包括:热和热⼒学分析模块；流变计算模块；动⼒分析模块实现真时间历程的瞬时动⼒响应模拟；⽤C++编写的⽤户⾃定义本构模块开挖直⽴坡的喷射混凝⼟墙加⼟钉加固的模拟加（下）和不加（上）⼟⼯织物⼟坡的潜在破坏特征FLAC-3D（Three Dimensional Fast Lagrangian Analysis of Continua）是美国Itasca Consulting Goup lnc开发的三维快速拉格朗⽇分析程序，是⼆维的有限差分程序FLAC2D的扩展，能够进⾏⼟质、岩⽯和其它材料的三维结构受⼒特性模拟和塑性流动分析。

FLAC 讲义一、什么是FLAC1.1 FLAC之字义F(Fast)L(Lagrangian)A(Analysis of)C(Continua). Lagrangian相对于Eulerian为每一时阶(timestep)之位移在Lagrangian之公式中，需对网格之座标予以更新，而Eulerian之公式则不予更新。

1. 2 FLAC之运算流程1.3 FLAC 基本单元1.4 分析模式大小与RAM之关系1.5 单位1.6 正负号方向(1)应力－正号代表张力，负号代表压力(2)剪应力－详见下图，图中所示剪应力为正号(3)应变－正的应变表示伸长，负的应变代表压缩(4)剪应变－剪应变的正负号与剪应力相同(5)孔隙压力－孔隙压力永远为正(6)重力－正号的重力物质往下拉，负号的重力将物质往上提。

二、FLAC内建之组合律FLAC内建之组合律有：1．空洞模式(null model)使用于土壤被移除或开挖2．弹性模式3．塑性模式，包括a. Drucker －Prager modelb. Mohr－Coulomb modelc. ubiquitous－joint modeld. strain－hardening/softening modele. bilinear strain－hardening/softening modelf. double－yield modelg modified cam－clay model此外，另有选购(option)模式，包括：1. 动力模式(Dynamic Option)2. 热力模式(Thermal Option)3. 潜变模式 (Creep Option)使用者另可使用FISH语言去建构独特的组合律以符合所需。

三、FLAC－以命令为输入语法请查阅相关手册四、FLAC程式之使用步骤4.1 FLAC程式使用前准备步骤步骤1：依比例画出所欲分析之资料于纸上画出地点之位置、地层资料、并简标示距离及深度资料。

Flac3D常用命令及常规问题总结Flac 3D 常用命令总结new 不退出flac3d，新建一个新的分析计算任务generate 简写为ge 产生网格zone 简写为zo 指示为3维网格体brick 简写为b 指示为长方形的网格体size x y z 指示长方形网格体在x、y和z轴所划分的网格数plot 显示，输入之后回车会出现plot base/0 可以输入plot的任何子命令（add添加条目（surface 网格体颜色、axes 坐标系颜色）等），若回到flac3d命令状态，则按一次回车或输入qiut 按x ，y，z 小写字母可以旋转当前视图中的X，Y ，Z 轴，相应的大写字母则方向相反。

m 放大，M缩小。

Ctrl+G可以切换彩色图和灰色图。

Model 简写为mo 定义材料模型Mohr 简写为moh 指摩尔-库伦模型Property 简写为pro 定义模型的材料参数，材料参数包括：密度density 简写为de，体积模量bulk简写为bu，切变模量shear 简写为sh，内摩擦角friction 简写为fric，内聚力cohesion 简写为c 抗拉强度tension 简写为ten等Set 用于设置模拟条件或控制条件，如重力加速度gravity，mechanical 设置静态力学分析的有关参数简写为mec，force 最大不平衡力简写为fo等。

比如set mechanical force 50Initial 简写为in 初始化网格的相关值。

Fix 保持网格节点指定参数（速度、压力、温度等）的值不变。

Range 指定范围History 简写为h 采样或记录迭代时模型中变量的值，检验采样记录时可以在plot子命令中输入history 1 、history 2等，忘记编号时可以输入print history就可以看到什么记录号对应什么变量或参数Nstep 简写为n 迭代次数，默认值等于10，ubalance 简写为unb 最大不平衡力，gp 采样网格节点的有关值，zdisplacement简写为zdis z轴方向的位移Solve 简写为so 控制相关过程的自动时间步一般等值线图（包括位移、速度、应力、压力、温度等），可以在一个新的视图中显示，首先要在plot子命令中新创建一个视图，可以输入“create 视图名”，回车，继续输入add contour/bcontour disp/szz…（contour为在当前视图中显示等值线图简写为con ，bcontour 绘制指定区的等值线图简写为bcon，disp 指定为位移量等值线图，szz 指定为应力等值线图）任意剖面上的等值线图，可以在刚创建的视图中再创建一个视图然后输入剖面参数，比如“set plane dip 90 dd 0 origin 3 4 0”（plane 简写为p 设置一个剖平面，dip 为剖面倾角，dd 为剖面的倾向，origin 简写为o 为剖面中的一点）Boundary 简写为bo 在视图中增加面的边界线框Behind 简写为be 为当前视图后面要转换视图则在plot子命令中输入“current 视图名”、show开挖把网格体的材料模型设置成空（null）模型即可，命令如下：Model null range x y zSet large 把材料设置成大变形求解时计算不收敛则用step命令，step 步数节点的ID号和重心的全局坐标可以由下面的命令得到：Print gp position单元体的ID号和重心的全局坐标可由如下的命令得到：Print zone所有命令大小写一样，都可以附带若干个关键字和有关的值，格式如下：Command 〈keyword value…〉[keyword，value…]…命令、关键字、数字可以用数个空格隔开，也可以用（），=隔开，；为注释符，；开始到行末都是注释一个简单的命令行只能至多80个字符（包括注释），如果超过80个字符则在行末输入&，每个命令最多含1024个字符。

给FLAC.FLAC3D软件初学者的一些意见一、认识和了解FLAC/FLAC3D有关这两个软件的介绍，在这里不再赘述，网上有大量而详细的材料。

需要提示大家的是，这两者都是基于连续介质及有限差分原理的2D（或3D）程序，一般不用于断裂、崩塌、侵彻等不连续变形问题。

这两个软件的区别，字面上，FLAC一般指的是2D程序，FLAC3D当然是三维程序。

但是从功能上来说，2D的FLAC更强大，可以进行两相流的模拟，用于非饱和土的研究很适用；同时，2D的FLAC有比较成熟的GUI（图形化处理界面）和更高的执行效率，当然也更简单。

所以，在进行学术性研究时，若问题可以简化成二维，则尽量采用2D的FLAC进行研究和分析。

当然了，现在高版本的FLAC3D（V4.0及以上）也有了比较好用的GUI界面，但是大部分的命令还是要通过文本才能实现。

认识这两个软件还需要注意的问题，就是它们的适用性。

哪些问题可以用，哪些问题无能为力，这个在大家决定学习之前一定要了解清楚。

因为这些商业程序都比较复杂，想通过“速成班”快速掌握（注意：不仅仅是上手）基本上不是一件很容易的事情。

根据我学习和使用的体会，下列问题请慎重使用FLAC/FLAC3D程序：（1）小尺寸试样的动力试验模拟问题，如动三轴试验的模拟等；（2）复杂水头的长时间固结问题，如考虑水头实时变化，需要完全耦合模式来求解的长时间沉降变形计算；（3）长期动力分析问题，如高铁长期振动、交通荷载长期振动分析等。

这类问题，为何不能适用，这里就不再展开。

总之，以上的各类问题本身也是具有相当的分析难度，大家在着手此类研究时，一定要慎重。

二、如何快速地“上手”？学习比较复杂的数值软件，都有一个“上手难、进展慢”的阶段，所以要合理安排学习进程，既要在短期内就培养自己的学习兴趣，同时也需要在前期学习中能够获得“局部成功”的喜悦。

我想，应该多做“小题目”应该是个好办法。

比如在开始上手FLAC3D程序时，多用很少单元的模型来“试验”程序的功能，积累学习经验，而不是一开始就动辄上万的单元、各种复杂的边界条件，那样只会困住自己，丧失学习的兴趣。

FLAC3D常见命令与使用技巧1、FLAC3D常见命令：1.FLAC3D是有限元程序吗？答：不是！是有限差分法。

2.最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。

3.怎样看模型的样子？答：plo blo gro可以看到不同的group的颜色分布4.怎样看模型的边界情况？答：plo gpfix red5.怎样看模型的体力分布？答：plo fap red6.怎样看模型的云图？答：位移：plo con dis (xdis, ydis, zdis)应力：plo con sz (sy, sx,sxy, syz, sxz)7.怎样看模型的矢量图？答：plo dis (xdis, ydis, zdis)8.怎样看模型有多少单元、节点？答：pri info9.怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件10.怎样调用一个文件？答：File/call或者call命令10.如何施加面力？答：app nstress11.如何调整视图的大小、角度？答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键12.如何进行边界约束？答：fix x ran（约束的是速度，在初始情况下约束等效于位移约束）13.如何知道每个单元的ID？答：用鼠标双击单元的表面，可以知道单元的ID和坐标14.如何进行切片？答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片)15.如何保存计算结果？答：save +文件名.16.如何调用已保存的结果？答：rest +文件名；或者File / Restore17.如何暂停计算？答：Esc18.如何在程序中进行暂停，并可恢复计算？答：在命令中加入pause命令，用continue进行继续19.如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步20. Fish是什么东西？答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学？答：可以不用，需要的时候查Mannual获得需要的变量就可以了22.FLAC3D允许的命令文件格式有哪些？答：无所谓，只要是文本文件，什么后缀都可以23.如何调用一些可选模块？答：config dyn (fluid, creep, cppudm)24 .如何在圆柱体四周如何施加约束条件？答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2radius r2 not，其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用set log on命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

陈老师好，请问flac能模拟地裂缝吗？还有就是断层的上下2盘该怎么融合呢？答地裂缝问题很难，因为FLAC本身是连续介质的理论。

回复收起回复2楼2013-06-21 00:25举报|个人企业举报垃圾信息举报本楼含有高级字体我也说一句若惜青吧主12问陈教授您好，FLAC作动力分析采用的吸收边界效果怎么样，不知您是否关注过？答效果很好啊，欢迎使用。

我就是做动力分析的回复收起回复3楼2013-06-21 00:26举报|个人企业举报垃圾信息举报本楼含有高级字体我也说一句若惜青吧主12问老师你好还是设置SHELL单元的问题主要由于表面不规则范围不知道怎么确定比如要在刷坡体表面设置shell 该怎么设定范围答不规则没有关系，你只要找到正确的range，程序会自动识别这个range范围内所有的“面”，就可以建立正确的shell了。

问最大最小主应力迹线用哪个命令显示呢？用箭线的方向和长度反应应力的大小和方向答plot stensor陈老师您好隧道开挖后的最大最小主应力云图有什么作用呢答有助于判断周围土体的应力状态，包括大主应力方向，了解应力集中的区域，以应对周围土体破坏等工程问题。

问现在的研究生论文貌似都多多少少有点数值模拟，很多是不用模拟都知道结果的，这...答“多多少少”、“很多”，概念太模糊。

既然你已经知道这样，所以你应该选择一些未知的、有重要意义的问题来做数值模拟，ok？问再问陈老师，在模拟深部构造应力时，除了您书上说的SB法外，我个人想了一个思路：即模型四周加构造应力边界条件，底面固定，顶面施加一个应力边界条件来反演埋深自重应力，这样计算可以吗？答是否可行，一试便知。

有新想法很好，但是要对该方法的正确性进行验证方可。

问陈老师，您好！在您的PPT中提及：接触面有三种工作模式-粘结界面、粘接滑移、库仑滑移。

请问，在命令流中，分别控制哪个参数；其工作原理又是怎样的，谢谢！答1、三种模式主要是根据接触面参数来确定的，建模的命令都一样，但是参数赋值不同。

3D的动力（爆破）分析问题我刚在做一个隧道爆破掘进的模拟，想这样做：q0HG tZ9u!tIl1.先模拟原岩应力状态-V0pvT\mo mohr v)p4FvX;岩体参数V%ws~'c r}T!]1r7F+wpro bulk 4.2e6 she 9.23e5 fric 12 coh 1e4 ten 9.6e3 range z 18.1 17.2ini density 1800 range z 18.1 17.2pro bulk 7.02e6 she 1.83e6 fric 7.2 coh 4.7e3 ten 4.7e3 range z 17.2 13.6 {+u(gHf T^k\ ini density 1650 range z 17.2 13.6 R4kg8ek'_3\bpro bulk 9.9e6 she 2.82e6 fric 29 coh 0 ten 0 range z 13.6 7.1 E l Ny U7T X"K K+Ptini density 1900 range z 13.6 7.1f/IL;I)s1\-g.Lpro bulk 4.37e7 she 1.46e7 fric 30 coh 3e4 ten 3e4 range z 7.1 6.4ini density 1990 range z 7.1 6.4 J-q^!e&c \!cprop bulk 4.95e7 she 1.65e7 fric 38 coh 5e4 ten 5e4 range z 6.4 -0.4ini density 2100 range z 6.4 -0.4HB$O-g}gA#tM&i\prop bulk 6.52e7 she 2.67e7 fric 40 coh 4e4 ten 4e4 range z -0.4 -9.1ini density 2480 range z -0.4 -9.1PS]r0Vu0v4N/{prop bulk 6.78e7 she 3.12e7 fric 43 coh 5e4 ten 5e4 range z -9.1 -17.6ini density 2510 range z -9.1 -17.6prop bulk 6.53e8 she 3.37e8 fric 44 coh 5e6 ten 5e6 range z -17.6 -26.9ini density 2520 range z -17.6 -26.9;设定边界条件(上边界z自由)Y$o?W:e!Lmfix z range z -26.8 -270?cN9c Hmfix x range x -39.9 -40.1fix x range x 39.9 40.1,d%l2mJP QX"YRfix y range y -0.1 0.1fix y range y 39.9 40.19s@u @-e&[N:D3zset grav 0 0 -9.8pj;w#^J_?3L;施加地面面状荷载#O+L t6v)~gapply szz -50e3 range z 18 18.2sol rat 1e-5save geom.sav2.模拟爆破掘进new&L c-C1ayu)W,EuMrest geom.sav ]]^+d!B2W4@config dynset dyn onini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0kTG$[0hNah-~ini state 0def setup#Jt S~R~ {omega = 2.0 * pi * freq!F!d0Q(P[j8Pendsetup|e4GW){3YC;y向速度时程曲线z \\-a8m'kdef wavewave = 1/2.718282^(1.5*pi*dytime)*ampl*sin(omega * dytime) ;大概震动1.3sendset freq = 40.625 ampl = 0.04765setup;p1为隧道掘进面apply yvel=0.0 hist=wave range gro p1 y -0.01 0.01N%y&l(ue;k1@)}5U/]O-a;设定粘滞边界8{VG E$G2S H;b Dapply xquiet yquiet ran z -26.8 -27]q:A9d*O5gapply xquiet yquiet ran x -39.9 -40.1apply xquiet yquiet ran x 39.9 40.1apply xquiet yquiet ran y -0.1 0.1apply xquiet yquiet ran y 39.9 40.1B4Y;L(t*c8I3a1p2QOhist dytimehist waveset dyn damp local 0.15a;DfR%]V1SH"Rsol age 1.0m null range group p1 y 0 1sol age 2.0:B;z$n*j'i-J%d3U{-k;之后就是循环的向前掘进，施加爆破波，掘进……请大家帮忙看看这些命令流有没有什么错误有这样几个问题：（用flac没多久，也是刚刚才接触动力模块，所以应该都是很浅显的问题，希望大家不要不屑于赐教哈）1.施加动力荷载时，apply yvel=0.0 hist=wave range gro p1 y -0.01 0.01中的yvel=？我看好多例子似乎都是=1.0，不知道这个赋的值代表啥？EN2g7@~O o5CA&]l2.同样是这一句，apply yvel=0.0 hist=wave range gro p1 y -0.01 0.01，最后那施加的位置好像多是一个面吧，能否施加在一定体积内呢？比如range gro p1 y 0 1？或者一条线上（线装炸药）？3.粘滞边界的设定跟之前静力计算中的约束条件有冲突么？比如fix z range z -26.8 -27与apply xquiet yquiet ran z -26.8 -27？r1C u"N9w4.如果同时施加x、y、z三个方向的速度时程，是否就是这样写：9|)~+oa$Fapply xvel=0.0 hist=wave1 range gro p1 y -0.01 0.01O"p*|-^2P Bbapply yvel=0.0 hist=wave2 range gro p1 y -0.01 0.01apply zvel=0.0 hist=wave3 range gro p1 y -0.01 0.013@:^bl(p BkU6[M(f;S?主要就是设置粘滞边界条件和apply动力荷载这两个地方就是搞不清楚，希望版主或者明白的朋友能给些指点，谢谢啦。

flac3‎d常用命令‎1、最先‎需要掌握的‎命令有哪些‎？答：需‎要掌握 g‎e n, i‎n i, a‎p p, p‎l o, s‎o lve ‎等建模、初‎始条件、边‎界条件、后‎处理和求‎解的命令。

‎2、怎样‎输出模型的‎后处理图？‎答：Fi‎l e/Pr‎i nt t‎y pe/J‎p g fi‎l e，然后‎选择 Fi‎l e/Pr‎i nt，将‎保存格式选‎择为 jp‎e文件。

‎3、怎样‎调用一个文‎件？答：‎F ile/‎c all ‎或者 ca‎l l 命令‎4、如何‎施加面力？‎答：ap‎p nst‎r ess ‎5、如何调‎整视图的大‎小、角度？‎答：综合‎使用 x,‎y, z‎, m, ‎S hift‎键，配合‎使用Ct‎r l+R，‎C trl+‎Z等快捷‎键。

6、‎如何进行边‎界约束？‎答：fix‎x ra‎n（约束‎的是速度，‎在初始情况‎下约束等效‎于位移约束‎）。

7‎、如何知道‎每个单元的‎ID？‎答：用鼠标‎双击单元的‎表面，可以‎知道单元的‎ID 和‎坐标。

8‎、如何进行‎切片？答‎：plo ‎s et p‎l ane ‎o ri (‎点坐标) ‎n orm ‎(法向矢量‎) plo‎con ‎s z pl‎a ne (‎显示 z ‎方向应力的‎切片) 9‎、如何保存‎计算结果？‎答：sa‎v e +文‎件名 10‎、如何调用‎已保存的结‎果？答：‎r est ‎+文件名；‎或者 Fi‎l e / ‎R esto‎r 11、‎如何暂停计‎算？答：‎E sc 1‎2、如何在‎程序中进行‎暂停，并可‎恢复计算？‎答：在命‎令中加入‎p ause‎命令，用‎cont‎i nue ‎进行继续。

‎在我们分‎步求解中想‎得到某一个‎过程中的结‎果，不用等‎到全求完，‎还可以在‎分布求解错‎误的时候就‎进行改正，‎而不是等到‎结果出来。

‎13、如‎何跳过某个‎计算步？‎答：在计算‎中按空格键‎跳过本次计‎算，自动进‎入下一步‎14、Fi‎s h 是什‎么东西？F‎i sh 是‎否一定要学‎？答：是‎FLAC‎3D 的内‎置语言，可‎以用来进行‎参数化模型‎、完成命令‎本身不能‎进行的功能‎。

1、FLAC3D常见命令：1.FLAC3D是有限元程序吗？答：不是！是有限差分法。

2.最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。

3.怎样看模型的样子？答：plo blo gro可以看到不同的group的颜色分布4.怎样看模型的边界情况？答：plo gpfix red5.怎样看模型的体力分布？答：plo fap red6.怎样看模型的云图？答：位移：plo con dis (xdis, ydis, zdis)应力：plo con sz (sy, sx,sxy, syz, sxz)7.怎样看模型的矢量图？答：plo dis (xdis, ydis, zdis)8.怎样看模型有多少单元、节点？答：pri info9.怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件10.怎样调用一个文件？答：File/call或者call命令10.如何施加面力？答：app nstress11.如何调整视图的大小、角度？答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z 等快捷键12.如何进行边界约束？答：fix x ran（约束的是速度，在初始情况下约束等效于位移约束）13.如何知道每个单元的ID？答：用鼠标双击单元的表面，可以知道单元的ID和坐标14.如何进行切片？答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z 方向应力的切片)15.如何保存计算结果？答：save +文件名.16.如何调用已保存的结果？答：rest +文件名；或者File / Restore17.如何暂停计算？答：Esc18.如何在程序中进行暂停，并可恢复计算？答：在命令中加入pause命令，用continue进行继续19.如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步20. Fish是什么东西？答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学？答：可以不用，需要的时候查Mannual获得需要的变量就可以了22.FLAC3D允许的命令文件格式有哪些？答：无所谓，只要是文本文件，什么后缀都可以23.如何调用一些可选模块？答：config dyn (fluid, creep, cppudm)24 .如何在圆柱体四周如何施加约束条件？答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2radius r2 not，其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用set log on命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

FLAC3D与FLAC常见问题的整理FLAC/FLAC3D常规问题的整理1．FLAC3D命令的FAQlakewater整理看到其它板块上都有这个FAQ，也就是常见问题问答，今天抽了时间进行了整理，想到了就写下来了，因为看到很多初学者费了很多的时间，但是还是没有将常用的命令掌握，所以这个也可以作为入门的初级教材，使大家能够快速的上手，而不用为了某个小命令到处求助。

1. FLAC3D是有限元程序吗？答：不是！是有限差分法。

2. 最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。

3. 怎样看模型的样子？答：plo blo gro可以看到不同的group的颜色分布4. 怎样看模型的边界情况？答：plo gpfix red sk5. 怎样看模型的体力分布？答：plo fap red sk6. 怎样看模型的云图？答：位移：plo con dis (xdis, ydis, zdis)应力：plo con sz (sy, sx, sxy, syz, sxz)7. 怎样看模型的矢量图？答：plo dis (xdis, ydis, zdis)8. 怎样看模型有多少单元、节点？答：plo info（错，应该为print info）9. 怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件10. 怎样调用一个文件？答：File/call或者call命令10. 如何施加面力？答：app nstress11. 如何调整视图的大小、角度？答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键12. 如何进行边界约束？答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束）13. 如何知道每个单元的ID？答：用鼠标双击单元的表面，可以知道单元的ID和坐标14. 如何进行切片？答：plo set plane ori (点坐标) norm (法向矢量)plo con sz plane (显示z方向应力的切片)15. 如何保存计算结果？答：save +文件名16. 如何调用已保存的结果？答：rest +文件名；或者File / Restor17. 如何暂停计算？答：Esc18. 如何在程序中进行暂停，并可恢复计算？答：在命令中加入pause命令，用continue进行继续19. 如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步20. Fish是什么东西？答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学？答：可以不用，需要的时候查Mannual获得需要的变量就可以了22. FLAC3D允许的命令文件格式有哪些？答：无所谓，只要是文本文件，什么后缀都可以23. 如何调用一些可选模块？答：config dyn (fluid, creep, cppudm)后注：这个工作很繁琐，需要的时间很多，希望广大网友能够将自己曾经遇到的常见问题在后续跟贴，也为了将这个FAQ进行很好的充实。

flac3d学习有一段时间，过一段时间没学究忘了很多，所以学习最好不要放下来。

对自己总结，也是学习的一种。

flac3d命令很多，用到的不多，就那么几个，基础的就不多说了进入正题：个人认为如果学习flac3d最好把flac的命令进行总结一下。

做模型就好比拼盘：要地应力，就把地应力的那段命令摘入，要出图，加入plot的相关命令，要出动画，加个动画命令，所以总结好命令后，以后做模型，只要把命令组一条条的对应放入，一个程序就出来了，后期调试一下，改相关参数就可以了。

接下来，我把自己弄好的整理一下，按照我上面的思路进行建模什么人建模都会有自己的喜好，下面说一下我自己建模第一步就是建立一个全新的txt文件（不论命令大小都是如此），这样的好处是，模型可以修改方便快捷。

而且还可以自我查错建立新的txt后，第一行命令就是new，加new的目的是为了防止其他命令的影响（相当于刷新）再执行其他命令前刷新一下。

不用想，接下来就是建模了，建立模型，不同人有不同的方法，有的人喜欢建完后在对模型开挖最终得到自己想要的模型。

model null 等；有的人喜欢用各种带缺角的组合成最终的模型。

具体建模看各自的喜好，我觉得，不同点是建模的网格划分有区别，用带缺角的组合模型网格划分更好最好把常用的集中模型记下来组成各种句子，以后用的时候直接插入例如：gen zone brick/radbrick size *,*,*,* ratio *,*,*,*建模的过程中，或同一个文件中有重复应用的命令，就可以用宏命令了，宏命令macro 名称‘命令的内容’例：macro shiping ’plot con zdip outline on’用shiping 代替后面的引号里的东西。

在建模中只要有相同的都可以用shiping代替，缩减了书写命令的时间，整体上也能让人看的更明白注意点是，你命名的符号最好跟falc自带的命令不要重复，我的处理方法是，名称后带数字如shiping1如果你想修改宏对象，改变命令的内容，在名称上加单引号就是替换的意思例macro ’shiping1‘‘新的命令’这些很多课本里都有讲，不在叙述很多时候我们要对单一的对象进行操作，例如进行模拟运算后，只想单独的显示一部分的内容，单一的土层，锚杆，等就要用到groupgroup给对象命名，例：在gen zone brick ------ size ---- group name（name随便取，不要跟flac的已有的命令重合就行）这里的意思就是对这个体进行了命名接下来想把这一块模型挖掉加model null range group namegroup的用法：1。

网上FLAC3D学习笔记1.请问在圆柱体四周如何施加约束条件？可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2<r1，其实就是选择range的时候选两个圆柱的差，即得到边界。

命令流如下：newgen zone cyl p1 1 0 0 p2 0 4 0 p3 0 0 1 size 12 30 12gen zone refle norm 1 0 0gen zone refle norm 0 0 1mo ssprop den 2500 bulk 2e8 she 1e8 co 2e6 fric 45 ten 1e6 dil 10prop ftab 1 ctab 2 dtab 3table 1 0 45 0.05 42 0.1 40 1 40table 2 0 2e6 0.05 1e6 0.1 5e5 1 5e5table 3 0 10 0.05 3 0.1 0fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 0.5 0 0 end2 0.5 4 0 rad 0.5 回答：加上的约束不对！fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 0.5 0 0 end2 0.5 4 0 rad 0.5 回答：上面的rad 是不是相差太大了？是不是在1左右来限制了2.关于计算不收敛问题;初始网格gen zone radcyl p0 15 0 -15.75 p1 add 15 0 0 p2 add 0 6 0 p3 add 0 0 15.75 &size 3 4 8 5 dim 2.25 2.25 2.25 2.25 rat 1 1 1 1.1 fillgen zone radtun size 3 4 3 5 dim 0.75 2.25 0.75 2.25 rat 1 1 1 1.1 &p0 15 0 -15.75 p1 add 0 0 -14.25 p2 add 0 6 0 p3 add 15 0 0 fillgen zone reflect dip 90 dd 270 ori 15 0 -15.75; 摩尔库仑模型model mohrmacro niyan 'bulk 3.2825e3 shear 1.5150e3 fric 20 coh 0.07'macro meiceng 'bulk 2.1735e3 shear 8.8917e2 fric 19 coh 0.05'prop niyan range z =-13.5 0prop meiceng range z =-15.5 -13.5prop niyan range z = -30 -15.5;为开挖定义整体部分range name in z -15.75 0 cyl end1 15 0 -15.75 end2 15 6 -15.75 rad 2.25 range name hangdao x 12.75 17.25 y 0 6 z -16.50 -15.75;进行加载fix x y z range z -30.001 -29.999app sxx=-20 range x -0.001 0.001app sxx=-20 range x 29.999 30.001app syy=-20 range y -0.001 0.001app syy=-20 range y 5.999 6.001app szz=-20 range z -0.001 0.001solve （1）model null range inmodel null range hangdaorange name in z -15.75 0 cyl end1 15 0 -15.75 end2 15 6 -15.75 rad 2.25 range name hangdaoliangce1 x 12.749 12.751 y 0 6 z -16.50 -15.75 range name hangdaoliangce2 x 17.249 17.251 y 0 6 z -16.50 -15.75 range name hangdao x 12.75 17.25 y 0 6 z -16.501 -16.499apply nstress -0.05 plane dip 90 dd 0 range in; 给巷道底部加压力app nstress -0.05 plane dip 90 dd 0 range hangdao; 给巷道直墙部分（左侧）加压力apply nstress -0.05 plane dip 90 dd 0 range hangdaoliangce2; 给巷道直墙部分（右侧）加压力apply nstress -0.05 plane dip 90 dd 0 range hangdaoliangce1;对点进行监测hist unbalhist gp xdisp 12.75 0 -15.75 :监测左帮hist gp zdisp 15 0 -13.5 ;监测顶板hist gp zdisp 15 0 -16.5 ;监测底板solve问：第一个solve平衡了,是第二个solve不能平衡，这是为什么？如果我把支护力加大，也不能平衡回答：我遇到的不收敛的原因，有两个：1）检查你的模型是不是没有建好，相互之间是不是没有连接好；2）仔细检查你的参数，是不是已经发生了破坏。