Flac3D命令--完整经典版

实例分析命令：

1. X ，Y ，Z 旋转 Shift+ X ，Y ，Z 反向旋转

Gen zone ……；model ……；prop ……（材料参数）；set grav 0,0,-9.81（重力加速度） plot add block group red yellow 把在group 中的部分染成红色和黄色

plot add axes black 坐标轴线为黑色；print zone stress% K 单元应力结果输出 ini dens 2000 ran z a b （设置初始密度，有时不同层密度不同）；ini ……（设置初始条件）；fix ……（固定界面）

set plot jpg ；set plot quality 100 ；plot hard file 1.jpg 图像输出（格式、像素、名称） plot set magf 1.0视图的放大倍数为1.0；plo con szz z 方向应力云图

2. ini z add -1 range group one 群one 的所有单元，在z 方向上向下移动1m ；然后合并

命令 gen merge 1e-5 range z 0此命令是接触面单元合并成一个整体，1e-5是容差

3. (基坑开挖步骤)：Step 1: create initial model state （建立初始模型）Step 2: excavate

trench （开挖隧道）

4. group Top range group Base not 定义（群组Base 以外的为）群组Top

5. plot blo gro 使得各个群组不同颜色显示

6. （两个部分间设置界面；切割法）：gen separate Top 使两部分的接触网格分离

为两部分；interface 1 wrap Base Top 在（Base 和Top ）这两部分之间添加接触单元；plot create view_int 显示，并创建标题view_int ；plot add surface 显示表面；plot add interface red 界面颜色红色

7. （简单的定义函数及运行函数）new ；def setup 定义函数setup ；numy = 8定义常

量numy 为8；depth = 10.0 定义depth 为10；end 结束对函数的定义；setup 运行函数setup

8. （隧道生成）上部圆形放射性圆柱及下部块体单元体的建立，然后镜像。

9. 模拟模型的材料问题时为什么要去定义某个方向上的初始速度—

10. 渐变应力施加：apply nstress -1e6 gradient 0,0,1e5 range z 3.464,0 plane dip 60

dd 270 origin .1 0 0；施加法向应力：apply nstress -1e6 range plane dip 60 dd 270 origin .1 0 0

11. d ip dd 确定平面位置使用：（纠结）

12. p rint gp position range id=14647 输出节点坐标

13.

a pply sxx -10e6 gradient 0 , 0, 1e5 range z -100 , 0在这个求解方程中，z 为变量，所以xx σ为：65=-1010+10xx z σ⨯⨯ ；原点（0,0,0）

14. f ree x range x -.1 .1 z 6.9 10.1放松x=0 平面上，z=7，10 这一部分在x 方向的约

束（可以在此处产生破坏）

15. 体积模量K 和剪切模量G 与杨氏模量及泊松比v 之间的转换关系如下：

16. 一般而言，大多数问题可以采用FLAC 3D 默认的收敛标准（或称相对收敛标准），即当体

系最大不平衡力与典型内力的比率R 小于定值10-5；（也可由用户自定义该值，命令：

SET mech ratio ）

17. 所谓体系最大不平衡力，是指每一个计算循环(或称计算时步)中，外力通过网格节点

传递分配到体系各节点时，所有节点的外力与内力之差中的最大值；

所谓典型内力，则是指计算模型所有网格点力的平均值。

18.要遵循由简单到复杂、由少网格到多网格的思路。

19.混合离散化法的基本原理是通过适当调整四面体应变率张量中的第一不变量，来

给予单元更多体积变形方面的灵活性。

20.g en zone ref以z=0平面为对称面镜像生成网格; gen zone ref gen zone ref dip 90以y=0面

为对称面镜像生成网格; gen zone ref dip 90 dd 90以x=0平面为对称面镜像生成网格21.球体模型建立：make_sphere函数循环使用。计算球心至节点距离、获取节点坐

标、调整坐标值得到球面至球心距离、完成八分之一球模型。变换节点内存地址！

22.①变形云图：（位移、应力、速度等等）；plot add cont szz out on shade on打开szz

的等值线图（outline、effective网格、有效的…显示）；PLOT con szz out on ；

PLOT con zd ou on magf 20（块云图bcon不能跟magf放大倍数）；plo con zdisp outline on shade on 打开zdisp的等值线云图

②变形矢量图：（变形方向、大小PLOT sk dis scale 0.07其中scale用来改变箭头大小，

默认是0.05）；

③塑性区分布print zone state都可以用来显示模型的塑性区

命令显示那些应力符合屈服准则的区域（或称塑性区）

应力正位于屈服面，或说正处于破坏状态时以shear-n 或tension-n 标识；

曾进入过屈服状态，但现已经退出以shear-p 或tension-p 标识

根据塑性区标识判断破坏机制是否在起作用，可按下述步骤进行：

④变量监测（历史跟踪）：Hist gp zdisp 1；Hist gp szz 2输出2和1的hist关系图

用：plot his -2 vs 1 如果应力值都为负数，则加-2是为了应力坐标轴为正！默认id=1 2 3 4 …依次（hist id=? Gp zdis x y z）；hist write 7 v 8 file 6-3hist.txt使用该命令，程序会在默认文件夹中生成一个名为6-3hist.txt 的文本文件，读者可以打开这个文本文件观察所得的数据结果

⑤切片设置：剖面（plot set plane ori 0 1.5 0 norm 0 1 0）、云图命令后面加plane、

加网格（ske）加矢量(dis)、坐标系(axe)等修饰效果。

23.输出信息：print zone stress 在命令窗口中会输出模型全部单元的6个方向应力数值

print gp dis 该命令会显示模型中全部节点三个方向的变形大小

SET log on 在命令行中使用上述命令后，程序在默认目录自动建立一个名为flac3d.log的文件，并开始记录命令窗口中的所有信息，直到用户设置log状态关闭（set log off）为止rest 6-3.sav

set log on

set logfile 6-2.log

print zone stress

print gp dis range id a any id b any

set log off 之后重新打开6-2.log

24.最大不平衡力不代表没有了塑性流动，只是不平衡力趋于零，并不完全归零。

记录不平衡力随步数的关系。

hist unbal监测不平衡力，并保留历史记录；save t1.sav保存到文件t1；

hist gp zdisp 4,4,8监测网格坐标点（4，4，8）在z方向的位移，并保留历史纪录