FLAC3D学习笔记(自己总结版)

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。

导入dat、txt文件的方法
导入.sav文件的方法
输入命令流：Restore 文件名.sav
导入.f3prj文件的方法
视图窗口关闭后，如何再次呈现输入命令plot zone
~
当命令流输错可从这里撤销
点亮相应命令流后才可撤销
透明度的调整
保存项目文件
保存文件的重命名做切片
]
倾向和倾角
法向和过一个点
删除切片
Ctrl+R为查看模型后还原
显示接触面，双击I nterface显示接触面
接触面属性项被激活后，可以更改属性
隐藏网格显示
输出网格
！
如果点击了透明，将不显示Information
只要不点击透明，则显示Information
等同于
点击显示节点坐标
fish简介
@之前要有（空格）
而在之间没有（空格）
~
按下列步骤取消fish语言前要加@
点击fish
软件的刷新
如何显示软件最终的安全系数
显示软件节点.
查看应力及最大剪应力的步骤
:
查看位移云图的步骤
输出云图
设置出图边宽
设置分辨率可设为1080
Reflect镜像restore调用之前的计算结果
Origin原点normal法向量（该点与原点的坐标差）
`
对于结构单元如桩，若点击Zone将不出现任何模型
按如下步骤才可显示
使用该命令后保存的文件与说建模的文件保存在同一文件夹
查看约束反力的步骤
查看分组的步骤
;
位移云图的显示颜色相反如何调整
连接两个不同尺寸的单元。

FLAC3D常用命令总结第一篇：FLAC3D常用命令总结1、模型建好如何检验所建模型是否有问题注意：第一步建模一定不要建太大的模型，可以建非常小的模型，确保模型的准确性，这样可以尽快开始后边的步骤，等后边计算分析步骤快速做完验证没有大问题，再将准确的模型建好换上，重新计算得到更加准确的结果。

2、如何检验应力边界、位移边界、加载速率等是否准确Plot add gpfix：看位移边界Plo add fap：看应力边界，看加载方向是否准确 Plot add vel：看速率，加载速率是否准确3、如何查看运算结果Plo con szz outline on：看z向应力云图，并显示边界网格 Plo add disp：看位移Plo con zdisp outline on magf 5：看z向位移云图，变形扩大5倍 Plo add bo ：添加边界 Plo blo state：看是否屈服Plo add stensor：查看大主应力方向Plo con smin：小主应力Plo con smax：大主应力Plo con ssi（shear strain increment）：查看剪切应变增量云图Plo bcont smin：查看单元小主应力4、如何做切面并查看切面结果在编辑“Edit“菜单中设置切面位置，在上述3命令后加上plane （pl）即可。

如：Plo con smin pl：查看切面小主应力5、后处理将模型导入到tecplot中做切面，加等值线等，对结果进一步处理，得到更漂亮的一些结果。

其他重要命令总结：1、如何在圆柱体四周如何施加约束条件？可以用fix...ran cylinder end1 end2 radius r1 cylinder end1end2 radius r2 not，其中r22、如何能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

FLAC3D5.0导入dat、txt文件的方法
FLAC3D5.0导入.sav文件的方法
输入命令流：Restore 文件名.sav
FLAC3D5.0导入.f3prj文件的方法
视图窗口关闭后，如何再次呈现？输入命令plot zone
当命令流输错可从这里撤销
点亮相应命令流后才可撤销
透明度的调整
保存项目文件
保存文件的重命名
做切片
倾向和倾角
法向和过一个点
删除切片
Ctrl+R为查看模型后还原
显示接触面，双击I nterface显示接触面
接触面属性项被激活后，可以更改属性隐藏网格显示
输出网格
如果点击了透明，将不显示Information
只要不点击透明，则显示Information
等同于
点击显示节点坐标
fish简介
之前要有（空格）之前加是Flac3D3.0和5.0区别而在之间没有（空格）
按下列步骤取消fish语言前要加
点击fish
软件的刷新如何显示软件最终的安全系数
显示软件节点.
查看应力及最大剪应力的步骤
查看位移云图的步骤
输出云图
设置出图边宽
设置分辨率可设为1080 Reflect镜像restore调用之前的计算结果
Origin原点normal法向量（该点与原点的坐标差）
对于结构单元如桩，若点击Zone将不出现任何模型
按如下步骤才可显示
使用该命令后保存的文件与说建模的文件保存在同一文件夹
查看约束反力的步骤
查看分组的步骤
位移云图的显示颜色相反如何调整
连接两个不同尺寸的单元。

FLAC3D学习1.FLAC3D的基本知识介绍岩土工程结构的数值解是建立在满足基本方程（平衡方程、几何方程、本构方程）和边界条件下推导的。

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

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

FLAC3D （FatLagrangianAnalyiofContinua）由美国Itaca公司开发的。

目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存（64K），所以，程序求解的最大结点数仅限于2000个以内。

1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。

因此，大大发护展了计算规模。

FLAC3D是一个三维有限差分程序，目前已发展到V2。

1版本。

FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。

因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。

FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。

调整三维网格中的多面体单元来拟合实际的结构。

单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发变形和移动（大变形模式）。

FLAC3D采用的显式拉格朗日算法和混合-离散分区技术能够非常准确发模拟材料的塑性破坏和流动。

由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。

FLAC3D采用ANSIC++语言编写的。

1对模拟塑性破坏和塑性流动采用的是“混合离散法”。

学习flac3d的心得:经过几天的学习,我学会了很多关于flac3d的基础知识,我先从基本的符号入手,知道了很多词语的意思,如xdis-x的位移,xvel-x的速度,fri-摩擦角,bu sh coh -代表不同的模量,还有很多基本的词语.然后从最简单的如何生成网格学习,学会了简单网格的生成,如怎么生成立方体(长,宽,高为95 85 75),步鄹如下:gen zone bri p0 0 0 0 p1 95 0 0 p2 0 85 0 p3 0 0 75 p4 95 85 0 p5 0 85 75 p6 95 0 75 p7 95 85 75 size 2 3 4 ratio 1.2 1.2 3 如何在flac3d的软件显示立方体只需要在命令行输入plot block group 这样的命令.还有可以在图形上显示射线的,步鄹如下:gen zone rad 下面的如上所述,但是我对球体的形成还是不懂,再者我学习了解决实际工程问题的步鄹,第一,根据题目中的条件在纸上按比例画出图形,第二,计算出初始的应力平衡,第三,根据画出的图形在flac3d生成,第四,对生成的图形进行赋值还有定义dens ,grav等已知量,第五,进行数据的处理,得出想要的结果.最后我学习对数据的后处理,可以用图形得到各层的受力的情况,还可以得到很多对实际工程有用的数据.如输出xx轴的受力输入命令:plot con sxx(不显示网格) plot con sxx ou on(显示网格).下面例题:NewGen zon bri size 2 2 4Model mohrProp bu 3e7 sh 1e7 coh 1e10 fri 15Fix z ran z 0Fix z ran z 2Fix x ran x 0Fix x ran x 1Fix y ran y 0Fix y ran y 4Ini dens 3000Set grav 9.8Ini sxx -3e4 grad 3.7777e3 0 0 ran z 0 1 Ini syy -4e5 grad 4.6666e4 0 0 ran z 0 1 Ini szz -6e7 grad 5.8888e6 0 0 ran z 0 1 App nstress -4e7 ran z 1Hist id=1 gp zdis 0 0 0Hist id=2 gp zdis 1 0 0Hist id=3 gp xdis 0.5 0 0Hist id=4 gp ydis 2.5 0 0Hist id=5 zone szz 0.5 2.5 1Hist id=6 zone syy 0.5 1 1Hist id=7 zone sxx 0.5 3 1Save 1-1.savSolve。

FLAC3D学习笔记(自己总结版)一、建立模型笔记1键盘上的键的功能:X,Y,Z分别控制所建立的模型围绕X,Y,Z三个轴做逆时针旋转、如果打开大小写锁定键,分边按XYZ时建立的模型会做顺时针旋转、+号的功能是加大模型旋转的角度，-号的功能是减小模型旋转的角度。

2键盘上的上下左右四个键可以实现图形的向上向下,向左和向右的移动、3当在一个区域建立多个模型的时候,有时候为了镜像单个模块儿,需要预先编组,编组的方法如下:gen zon tunint p0 0 0 0 p1 010 0 p255 0 p3 0 05 p6 0105 dim sq2 sq21 sq2111 group(编组)z1(组名为:z1)&fill group tunnel4镜像整个模型的方法:首先判断即将镜像的模型在原模型的哪个方向(即三个轴的方向),然后在轴所在的方向上去两个点、其中一个是原点(origin(xyz))、5镜像部分模型的方法: gen zone reflect normal (x,y,z)origin (x,y,z)range group z1就可以了、6特别要注意的一点:交叉巷道的巷道充填和巷道的网格是两个组,所以在镜像的时候不要忘了给充填部分镜像、7 group radcyl range group cylint not (编组的时候不包括cylint)7 快捷键总结Ctrl+Z 局部缩放的功能。

Ctrl+R的功能是恢复到初始状态。

8模型建立的方法:方法1:利用默认值生成网格,各个默认值:(1)p0的默认值是(0,0,0),网格的每边的单元格数size默认值是10,网格每边的长度默认值为size的值、方法2:利用4个点的坐标来生成矩形网格、p0~p3 size方法3:利用edge来确定边长生成矩形网格、方法4:利用参数ratio来确定单元体几何边画边生成矩形网格、方法5:利用参数add(相对坐标)来生成矩形网格、方法6:利用八个点的左边来生成矩形网格、9、用户自定义模型的方法我来贡献一点自己的成果：; U" j/ ?7 L+ a0 uFLAC3D的二次开发环境提供了开放的用户接口，在软件安装文件中包含了软件自带所有本构模型的源代码，且给出了Mohr-Coulomb模型和应变软化模型的编译示例，因此可以方便地进行本构模型的修改与开发。

FLAC3D 软件使用心得体会张强 20104323 10级水利3班一,初识FLAVC3D软件：1，作用：它是一款数值计算软件，有二维和三维计算程序两个版本，在岩土工程领域用处很大，其特点是简单易懂，能够进行土质.岩石和其他材料的三维结构受力特性模拟和塑性流动分析。

调整三维网格中的多面体单元来拟合实际的结构。

2，主要应用：(1) 承受荷载能力与变形分析: 用于边坡稳定和基础设计(2) 渐进破坏与坍塌反演: 用于硬岩采矿和隧道设计(3) 断层构造的影响研究: 用于采矿设计(4) 施加于地质体锚索支护所提供的支护力研究: 岩锚和土钉的设计(5) 排水和不排水加载条件下全饱和流体流动和孔隙压力扩散研究: 挡土墙结构的地下水流动, 和土体固结研究(6) 粘性材料的蠕变特性: 用于碳酸钾盐矿设计(7) 陡滑面地质结构的动态加载: 用于地震工程和矿山岩爆研究(8) 爆炸荷载和振动的动态响应: 用于隧道开挖和采矿活动(9) 结构的地震感应: 用于土坝设计(10) 由于温度诱发荷载所导致的变形和结构的不稳定(11) 大变形材料分析: 用于研究粮仓谷物流动和放矿的矿石流动二．操作对象：Flac3D中为岩土工程问题的求解开发了特有的本构模型, 总共包含了10种材料模型:(1) 开挖模型null(2) 3个弹性模型(各向同性, 横观各向同性和正交各向同性弹性模型)(3) 6个塑性模型(Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模型).Flac3D网格中的每个区域可以给以不同的材料模型, 并且还允许指定材料参数的统计分布和变化梯度. 还包含了节理单元, 也称为界面单元, 能够模拟两种或多种材料界面不同材料性质的间断特性. 节理允许发生滑动或分离, 因此可以用来模拟岩体中的断层、节理或摩擦边界.FLAC3D中的网格生成器gen, 通过匹配、连接由网格生成器生成局部网格, 能够方便地生成所需要的三维结构网格. 还可以自动产生交岔结构网格(比如说相交的巷道), 三维网格由整体坐标系x, y, z系统所确定, 这就提供了比较灵活的产生和定义三维空间参数.三．操作命令：结合一些边界条件加上各个实际情况就可以使用该软件对岩石的力学等性质进行计算了。

Flac3d学习体会Flac3d是有限差分数值分析软件，它可以模拟岩土和其它材料的三维力学特性。

因此，在岩土工程力学中得到了广泛的应用，我们可以用该程序较好的模拟地质材料在达到强度极限或屈服极限时发生的破坏或塑性流动的力学特性，特别适用于分析渐进破坏失稳以及模拟大变形。

该软件中提供的强大丰富的本构模型可以进行静力、动力、蠕变、渗流、温度等多种计算模式，各种模式间可以互相耦合，也可以模拟多种结构形式，如梁、锚杆、土工格栅、衬砌、桩等。

作为一位初学者，主要了解了下flac3d中的弹性模型和莫尔库仑模型，通过一段时间的练习主要有以下心得；在分析工程案列时，我们可以依照flac3d的一般求解流程进行剖析案列，将命令流归纳为三个部分，即分析模型、模拟及求解、计算结果输出。

分析模型部分包括生成网格单元、设置初始条件和边界条件以及初始应力平衡；在模拟求解部分包括加载及连续建模、后续工程施工、场方程的有限差分求解；输出计算结果部分主要为图表的绘制、相关数据的输出等；通过这三个部分，我们可以写出一个完整的命令流，方便后续在运行过程中修改命令流以及查出具体在哪个环节出错了。

这里我们应该注意在flac3d的建立分析模型部分，材料性质的定义、初始条件和边界条件的设置并无明显的先后顺序。

对于工程实例中，我们在分析模型中，首先看是否对称，如果对称我们可以简化模型，然后用gen zone reflect dip 关键字dd 关键字ori 关键字这个命令来镜像模型。

这里我们需要注意的就是dip和dd的具体含义才能把模型建好；dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线，dd表示y轴正向顺时针到那条射线的夹角。

首先应该按照dd的方向大体确定这个面的朝向，dd 指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角．gen zone reflect dip 90 dd 90是以YZ面.gen zone reflect dip 90 dd 0 才是以xz为对称面的当初始应力平衡后，我们可以求工程的影响，求出区域内的应力分布情况后，再依工程的流程及步骤阶段执行各工程进行过程的影响，建议使用以下的步骤：(1)调出起初的应力平衡：re_____ .sav(2)设定新的材料性质：model,prop(3)设定新的支撑性质：struct(4)设定新的外力 apply(5)设定边界条件 apply,fix(6)求工程时的应力平衡 solve, step, cycle(7)储存 save file.sav在初始应力平衡过程中，我们要注意是否收敛，对于不收敛一般要检查模型是不是没有建好，相互之间是不是没有连接好；仔细检查参数，是不是已经发生了破坏；对于计算时间过长，我们可以适当修改模型的网格数以减少计算时间，或者使用step命令，这里的关键问题是要运行多少步。

FLAC3D各种命令笔记编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（FLAC3D各种命令笔记）的内容能够给您的工作和学习带来便利。

同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。

本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为FLAC3D各种命令笔记的全部内容。

hist reset ；清空里面所有的监测点nstress=-40e3 hist ramp的意思,很简单 nstress = —40e3 ＊ ramp hist相当于乘法，为了方便控制动荷载的幅值。

plotitems--Add--strucktural elements——geometry ——cables 显示锚杆的几何形状Pl-add-vector-apply body force 显示施加的力,不过是在求解之后显示Creat 创建一新视图plotitems—--—--add——-—-vetor-—-—apply body force 显示荷载的施加位置Array var(4,3）二维数组 array var(5)一维5个元素的数组Var（m,n）=urand ;随机 var—数组的名字Pnt=z_next(pnt) 查找下一个单元体表Pl plas 显示塑性区Array x_b(x，y） z_b（x,y) 创建两个x行y列的表格，可以给其赋值，以便调用.Array就是数组的意思。

Hist id=1 gp xd 7 0 0 监测7 0 0 点处的x方向位移，ID=1材料参数dilation 剪胀角FISH语言：Degrade;∏/180Pi;∏Ngp；节点总数Nzone;单元体总数单元节点遍历：P_z = zone_headLoop while p_z ＃ null………P_z = z_next(p_z)End loopRead（ar，n）;将数组ar的前n个记录写到文件print zone state 显示塑性区域编写fish函数时，尽量采用浮点型，例如1。

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

1. FLAC3D用户界面的进化自从2009年FLAC3D发布4.0版本以来，FLAC3D在操作上发生了巨大的变化，由原来的单纯以CALL命令执行数据文件的方式转变到以项目文件(Project)为核心的图形操作(GUI)方式。

尽管目前还保留着DOS系统下的命令行执行文件，但似乎没有太大的用处，因为在DOS系统下不再支持Plot命令。

2. FLAC3D版本的进化如同Itasca的其它产品一样，FLAC3D有一个检查程序版本的命令: List version，输入这个命令 可以查看FLAC3D的版本号码。

关注程序的版本号对于不经常使用的用户来说没有太多的意义，但是对于经常使用者，特别是对于数据文件的持续维护和二次开发者来说就显得非常重要了。

List version命令可以在CMD命令行状态下执行，如下图所示：也可以在GUI内命令行状态下执行，如下图所示：在Tool‐>Check Modifications下查看, 可以通过点击菜单Tools下面的Check Modifications来检查版本的在线更新历史。

3. FLAC3D的扩展功能FLAC3D除了集成在主程序中的基本功能外，还有一些选择性的模块，如动力学模块(Dynamic),蠕变模块(Creep), 热力学模块(Thermal)以及Barton‐Bandis模块等。

使用List Security这个命令可以查看硬件锁中存储的信息 (summary of information stored in the hardware-lock key)。

为了验证FLAC3D的选择性模块，在命令行输入Config 加上后面的关键字即可查看。

FLAC3D 5.0 REVISION HISTORYVisit the FLAC3D Updates pageV ERSION D ATE D ESCRIPTION12007/09/14◦Improve numerical stability for cysoil model when shear and cap yielding occursimultaneously in a cohesive material.◦Fix bug that might cause CaveHoek restriction to activate without reason.◦Fix possible invalid floating point errors for asin function in cysoil, chsoil ormhoek models.11906/25/14◦Added SEL INITIALIZE COUPLING command. This attempts to initialize forcesin coupling springs based on zone stresses, for all 2D elements in the range. 11806/19/14◦Added ADD and MULTIPLY optional keywords after SEL NODE APPLY FORCEand SEL NODE APPLY MOMENT.◦Made the GEOMETRY DIST range element more flexible when encounteringdegenerate polygons.◦Detect degenerate quad polygons on DXF import and make them triangles. 11705/15/14◦Add a keyword NOMERGE to GEN ZONE SPLIT so that newly created gridpointswill not merge outside the range11604/16/14◦Fix bug on INTERFACE and DFN plotting exporting.◦Add an optional keyword EXCLUDEQUAD for the GENERATE ZONE SPLITcommand;◦Face groups now be kept for the GENERATE ZONE SPLIT command.11504/03/14◦Add GENERATE ZONE SPLIT command - which will divide tetrahedral andwedge shaped zones into hexahedra.◦The range element GEOMETRY COUNT will now work with non-planar polygons. 11403/20/14◦Fix bug FISH defined flux apply conditions.◦Fix bug in anisotropic model that nu12 was not correctly copied in z_copyintrinsic◦Fix global segment mapping in GENERATE ZONE DENSIFY for some specialcases◦Extended gp_massadd() FISH intrinsic to allow an optional second argumentindicating degree of freedom being modified. This can now be used toimplement a Westergaard approximation in an arbitrary direction.11303/13/14◦Fix bug - water table density could have random value, and not give error thata value must be specified.11203/05/14◦Fix data file export of label plot.◦Fix command processing bug in ISOZONE plots.11102/13/14◦Isosurface plots that take a really long time to calculate should not cause theengine to timeout.11002/10/14◦Fix bug where the color of job or view title cannot be modified.◦◦Fix APPLY STRESS (and related) on wedge zone types, not correctly matchingunderlying tet distribution.◦Fix slot index bug in Z_FACEGROUP FISH intrinsic.◦Fix pore pressure distribution on gridpoints where to separate interfaces meet,when CONFIG FLUID is on.◦Add FISH functions of "baseline" and "filter" into the FISH library.◦Fix segment position zero comparison bug in extruder.◦Fix bug of overflow error in calculation due to an uninitialized vector.◦Added tests for join faces that already exists, indicating you are generatingzones inside each other.◦Enabled looking at gp local ratio values (plotting and fish) during cycling.◦Added the ability to set the maximum number of labels a label plotitem willshow, using MAXLABEL keyword or a control in Global Display Settings.◦Properties available to COLORBY PROPERTY in a zone plot are nowalphabetically sorted.◦During a NEW, FISH reference count deletion of Mem and Array objects isskipped since they will all be deleted anway. This speeds up a model resetwhen many of these objects exist.◦New gridpoints created by GEN SEPARATE copy the group and extra variablevalues from the original.◦Fix bug exporting octant cutting planes to a data file.◦Plotting user-defined vectors as disks now property supports the scalingattribute.◦'geom remove ?' no longer deletes all information in the current geometry set.◦Reset the 64 color palette.◦Properties on side 2 of embedded liners weren't initialized to zero.◦The *.f3prj desktop file association should be created during installation.◦Updated system to prevent recursive cycling, hopefully preventing incorrecterror.9902/26/13◦Fix bug that causes the associated flag in FOS solve to not work.◦Fix bug when a CableSEL node being rigid linking to a zone.◦Fix problem with demo version availability on 64 bit systems.◦Fix contour or vector plots of Out-Of-Balance-Force when attach conditions arepresent.◦Added local force ratio and log of local force ratio to gridpoint data used incontours, profile lines, zone field data, etc.◦Added fish function gp_localratio, returns local force ratio.◦Added LIST KEYINFO, to allow seeing security key values even if it isn'tcurrently valid.◦The MODEL LOAD command now returns a warning if the model has alreadybeen loaded, instead of causing an error.◦Updated the maximum number of labels in a label plot to 1024.◦Fixed error handling during DFN import, and fix problem with import ofFRACMAN files.9801/28/13◦GROUP range elements are no longer always case-sensitive.◦The extruder calculates a better default view for small geometric backgrounddata.◦The extruder now limits how much you can grow or shrink a backgroundimage.◦In the extruder, a right-click always brings up the "Select Geometric Data Set"and "Select Background Image" options even if not on the backgroundata.◦The extruder now shows an error message if the background image fails toload.◦The grid import dialog now includes an "All Files" filter option.◦The CHSOIL model now has initial friction angle as an input property, andplastic shear modulus as an output property.◦Prevent rendering small fonts on plots, as a workaround for a memorycorruption bug.9712/26/12◦Prevent label plots from accumulating more than 256 distinct labels.◦Changed how axes are rendered when fixed to a location on the screen (thedefault).◦Fix slot value when using z_facegroup FISH function.◦z_faceingroup FISH function had same slot error as z_facegroup, and was alsospelled z_faceisgroup internally.◦Optimized VRML output algorithms, greatly speeding up file generation forlarge plots.◦Activate DEMO mode.◦Add interactive mouseover information to the Face Group plotitem whentransparency is turned off.◦Cause error message if both NMD and multistepping are active at the sametime.◦Improve performance of GP based APPLY commands.9612/21/12◦If the console row limit was set to a very small number, FLAC3D could hang upon output.◦Prevent rendering of plots during restore of that plot.◦Fix how cell spaces were interacting with the notice system, causing freezeupsin the GUI.◦Zone contour plots smoothed by volume weighted averaging no longer usedeformed volume for the calculation.◦During plot export to DXF, suppress output of polygons in wireframe mode.◦During VRML export, reverse winding rule used by polygons, some importersdon't check the setting.◦Improve RANGE GEOMETRY COUNT accuracy when lots of surfaces are exactlyparallel to the projection vector.◦Attempt to merge face metadata (groups, extra variables) when faces aremerged by GEN MERGE.◦Fix bug in RANGE GROUP, greatly speeding performance of filtering.◦Changed method of calculation of tolerance used for interface contactdetection.◦Fix bug recording edgestate when exporting a geometry view to a datafile.◦ZGroup plots not also plot surfaces of the range applied.◦Cleaned up command processing of the two GEOMETRY RANGE elements. 9512/12/12◦Fix bug when exporting a view to a datafile that was specified by orientation.◦Fix bug causing possible crash when deleting range elements in the rangedialog.◦Fix bug causing a crash when using the geometry remove command.◦Add warnings if the range filters out all faces/gridpoints/zones for the APPLYcommand.◦The console row storage limit was not being saved to project file/registry.◦Fix problem with the dialog allowing you to enter a network key server name. 9411/29/12◦Added gp_stiff(p_gp) read/write FISH function for gridpoint stiffness.◦FISH fragments are now highlighted in the FISH symbol color in the editor.◦The CONTINUE keyword is now recognized inside FISH code in the editor.◦Fix potential crash restoring save files while the extrude view of the extrude pane is showing.◦Fix bug in z_find FISH intrinsic.◦Fixed bug in GEOM DIST and DFN DIST calculations.◦Fix bug rendering screen fixed axes directions.◦Fix bug in IMPGRID zone group name processing when leading spaces were present.◦ATTACH FACE TOL t now applies the tolerance (approx) to determining if Gps fall on faces for a valid attach.。

学习记录1.2009.8.16念及FLAC 学习起来很困难，且仅自己一人在学，无人交流，时常感到困惑迷茫，故写此FLAC-learn 来明确自己的目标。

此篇有三大作用：一来鞭笞自我，二来总结规纳，三来收集资料FLAC学习是岩土人一生的事情，不指望一下子学精，要知道在FLAC上学无止境，坚持学习。

八月份先初步整理好建模、后处理、接触面、结构构件、fish、边界条件、初始条件的学习总结、问题及资料（三本书及网上资料）先搭个架子，养成习惯，故一不深究，二不纠缠于生偏。

经过这个过程后，FLAC学习一阶段完成，达成目标：熟练掌握全部基本操作、基本命令，能独立处理简单问题，养成好的操作习惯，形成规范的FLAC思维方式。

以后的事情，以后再说，但FLAC的学习不应三天打鱼，两天晒网，要持之以恒。

高屋建瓴1.FLAC原理2.FLAC的特点3.FLAC的使用禁忌及大局观4.FLAC全局细节常用单位：符号规定：文件系统：建模1.常用的建模方法常用命令一，Print print gp position（range id 59）/print zone(显示节点或单元体的ID号和重心的全局坐标)二，plot Plot set rotation 20 0 30Plot set center autoPlot set magnification 1.5Plot add surface yellowPlot add axes我有话说：set 与add经常混肴，注意set是改变已有的，而add是加入里面没有的。

命令后可加range:Plot add surface yellow range x 1 2 y 1 2 z 1 2三，指定范围并命名的三个命令Range range name xxx x x1 x2 y y1 y2 z z1 z2 /model null range xxxGroup group xxx range x x1 x2 y y1 y2 z z1 z2/model null range group xxxMacro macro xxx ‘x x1 x2 y y1 y2 z z1 z2’ /model null range group xxx书中建议：一，把对象作为名字而不是有意义的实体看待二，对象名xxx注意与FLAC保留字的冲突我有话说：对macro宏命令不是很了解，以后要留意。

网上FLAC3D学习笔记21 请问在圆柱体四周如何施加约束条件？可以用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 10 prop 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.25range 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.25range name hangdaoliangce1 x 12.749 12.751 y 0 6 z -16.50 -15.75range name hangdaoliangce2 x 17.249 17.251 y 0 6 z -16.50 -15.75range name hangdao x 12.75 17.25 y 0 6 z -16.501 -16.499 apply 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）仔细检查你的参数，是不是已经发生了破坏。

FLAC3D5.0导入dat、txt文件的方法
FLAC3D5.0导入.sav文件的方法
输入命令流：Restore 文件名.sav
FLAC3D5.0导入.f3prj文件的方法
视图窗口关闭后，如何再次呈现？输入命令plot zone
当命令流输错可从这里撤销
点亮相应命令流后才可撤销
透明度的调整
保存项目文件
保存文件的重命名
做切片
倾向和倾角
法向和过一个点
删除切片
Ctrl+R为查看模型后还原
显示接触面，双击I nterface显示接触面
接触面属性项被激活后，可以更改属性
隐藏网格显示
输出网格
如果点击了透明，将不显示Information
只要不点击透明，则显示Information
等同于
点击显示节点坐标
fish简介
@之前要有（空格）之前加@是Flac3D3.0和5.0区别
而在之间没有（空格）
按下列步骤取消fish语言前要加@
点击fish
软件的刷新
如何显示软件最终的安全系数
显示软件节点.
查看应力及最大剪应力的步骤
查看位移云图的步骤
输出云图
设置出图边宽
设置分辨率可设为1080 Reflect镜像restore调用之前的计算结果
Origin原点normal法向量（该点与原点的坐标差）
对于结构单元如桩，若点击Zone将不出现任何模型
按如下步骤才可显示
使用该命令后保存的文件与说建模的文件保存在同一文件夹
查看约束反力的步骤
查看分组的步骤
位移云图的显示颜色相反如何调整
连接两个不同尺寸的单元。

FLAC3D学习笔记(自己总结版)一、建立模型笔记1键盘上的键的功能:X,Y,Z分别控制所建立的模型围绕X,Y,Z三个轴做逆时针旋转、如果打开大小写锁定键,分边按XYZ时建立的模型会做顺时针旋转、+号的功能是加大模型旋转的角度，-号的功能是减小模型旋转的角度。

2键盘上的上下左右四个键可以实现图形的向上向下,向左和向右的移动、3当在一个区域建立多个模型的时候,有时候为了镜像单个模块儿,需要预先编组,编组的方法如下:gen zon tunint p0 0 0 0 p1 010 0 p255 0 p3 0 05 p6 0105 dim sq2 sq21 sq2111 group(编组)z1(组名为:z1)&fill group tunnel4镜像整个模型的方法:首先判断即将镜像的模型在原模型的哪个方向(即三个轴的方向),然后在轴所在的方向上去两个点、其中一个是原点(origin(xyz))、5镜像部分模型的方法: gen zone reflect normal (x,y,z)origin (x,y,z)range group z1就可以了、6特别要注意的一点:交叉巷道的巷道充填和巷道的网格是两个组,所以在镜像的时候不要忘了给充填部分镜像、7 group radcyl range group cylint not (编组的时候不包括cylint)7 快捷键总结Ctrl+Z 局部缩放的功能。

Ctrl+R的功能是恢复到初始状态。

8模型建立的方法:方法1:利用默认值生成网格,各个默认值:(1)p0的默认值是(0,0,0),网格的每边的单元格数size默认值是10,网格每边的长度默认值为size的值、方法2:利用4个点的坐标来生成矩形网格、p0~p3 size方法3:利用edge来确定边长生成矩形网格、方法4:利用参数ratio来确定单元体几何边画边生成矩形网格、方法5:利用参数add(相对坐标)来生成矩形网格、方法6:利用八个点的左边来生成矩形网格、9、用户自定义模型的方法我来贡献一点自己的成果：; U" j/ ?7 L+ a0 uFLAC3D的二次开发环境提供了开放的用户接口，在软件安装文件中包含了软件自带所有本构模型的源代码，且给出了Mohr-Coulomb模型和应变软化模型的编译示例，因此可以方便地进行本构模型的修改与开发。