FLAC3D笔记
FLAC3D笔记
FLAC3D内置12种模型，静力、动力、蠕变、渗流、温度5种计算模式，各模式可以耦合作用。
初始菜单栏
File new 新建
call 读入dat、txt命令
model 读入自定义模型dll
restore 读入保存文件sav
save 保存
import grid 导入节点
print 打印
print type 设置输出图片格式
print setup 设置图片大小质量
exit 退出程序
display
options
plot creat 创建视图
1.base/0 show 显示基本视图
其他增加某种后处理结果，参考彭文斌附录
Window title 水平平铺窗口，上下切换图形窗口和命令窗口
cascade 层叠窗口
Help 详细的版本信息
当前菜单栏
Flie show 显示
print 打印
print type 设置输出图片格式
print setup 设置图片大小质量
information 屏幕信息
close 关闭
Edit camara 设置视图的属性如CENTER(中心);ROTATION（旋转角度）；DIST(距离)；MAG（放大倍数）等，快捷键ctrl+c
cutting plane 切图的属性，如过某点orgin、切面法向normal、倾角dip、倾向dd。快捷键ctrl+l
reset 视图还原为默认值，即为正平面视图。快捷键ctrl+r
zoom 用鼠标左键放大视图，快捷键ctrl+z
copy to clipboard 命令可以将显示的图形，存放到WINDOWS 剪贴板上（没有任何文件生成）。
该图形接着就可以以加强型图元文件格式被粘贴到其它兼容该格式的WINDOWS 应用程序中去。
点开画图，然后粘贴进去。
Setting Grayscale 设置黑白色，快捷键ctrl+G
caption 设置左侧说明栏的大小，默认是30
title 标题显示，比如写成title love，标题就是love
foreground 前景色，网格色
background 背景色
export 输出记录plot 命令为dat格式的文件，需要自命名为abc.dat
Plotitems add 增加某种后处理结果，参考彭文斌附录
clear 清除当前视图
1..... activate 激活视图
modify 修改视图参数，如outline、max、min、interval
remove 移除视图
Window title 水平平铺窗口，上下切换图形窗口和命令窗口
cascade 层叠窗口
1.云图coutour和块云图zone coutour输出：
变形云图：plot contour zdisp
plot contour zdisp magf 20 放大沉降为20倍，便于观察
应力云图：plot contour szz Z向应力云图
plot contour szz shade on 显示阴影
plot contour szz outline on 显示网络
plot contour szz effective on 显示有效应力
plot contour szz max -10e3 显示应力最大数值不
超过-10e3
plot contour szz min -40e3 显示应力最小数值不小于-40e3
2.变形矢量图
plot sketch disp scale 0.07， scale 默认0.05，数值越大矢量箭头越大
3.塑性区分布，判断破标准 陈育民书2.4节13页
plot block state 显示塑性区，两种类型剪切破坏shear failure和拉伸破坏tensile failure
正处于屈服面上或破坏状态shear-n或tension-n标识，now
曾进入屈服状态或破坏状态shear-p或tension-p标识，past
plot block state shear 包括shear-n，shear-p
tension-p 过去拉伸屈服单元
state now 现在
state past 过去
print zone state 在log里显示塑性区的区域
4.变量监测：hist放在加载应力状态之后、solve之前
history urban 监测不平衡力
history id 5 zdisp x y z 监测z沉降
plot hist5 显示
plot hist 1 vs 5 两个监测变量的关系
history write 1 vs 5 file ss.his，记录输出文件名ss的数值
hist gp vel x y z 网络节点速度，判断体系是否进入稳定状态，plot vel 绘制完整的速度矢量图，（判断收敛标准，陈书14页）
set hist_rep 1 监测记录间隔为1
5.结果输出print和日志log
print zone stress 输出应力数值
print gp dis 输出节点变形
set log on 建立log文件
set logfile loglove.log 建立名为loglove的log文件，print输出文件放入log日志中。
6.切片操作
切图的属性：过某点origin、切面法向normal、倾角dip、倾向dd。快捷键ctrl+l
plot set plane origin 0 1.5 0 normal 0 1 0
plot con zd plane
plot add sketch 增加网格线
plot add dis plane
plot add axes
7.动画制作
set movie avi step 1 file love.avi
movie start
solve
收敛标准 set mach force 100 （判断收敛标准，陈书14页）
generate zone add 以p0为原点建模型
group 网格组名
dimension 内部区域尺寸
p0-p16 网格控制点
edge 网格边长
ratio 相邻网格单元尺寸比率
fill 内部填充
size 坐标上网格单元数目
merge 检查导入模型的网格正确性
interface 1 prop kn=4e9 ks=4e9 tens=1e10 coh=15e3 fric=14.0 bslip=on sbratio=20
sbratio 抗剪强度设为抗拉强度的多少倍。
体集模量K、剪切模量G、杨氏模量E、泊松比v二元相关
施加外荷载 apply nstress -10e4 ran z3 x 1 2 y 1 2 施加z=3面向下法向应力、
“Plot clipboard”命令可以将显示的图形，存放到WINDOWS 剪贴板上（没有任何文件生成）。该图形接着就可以以加强型图元文件格式被粘贴到其它兼容该格式的WINDOWS 应用程序中去。点开画图，然后粘贴进去。
“Set plot metafile”命令可以将图形以加强型图元格式存盘，以便作
为计算的参考或日后插入到文档中去。
可以运用菜单命令：plot-print，然后保存为任意位置为格式PS的图片，用画图打开。
通过命令：Set plot +关键词（pcx, bitmap, bmp 或者jpg）可以存储为许多图像格式（pcx,bmp,jpeg 等）。输出的这些位图的分辨率由命名行：Set plot
无论是黑白的还是彩色的 postscript 打印机，都需要通过“Set plot postscript”命令来指定。打印图形将存储为文件，这样支持postscript 格式的图形处理程序就可以读入并进行修改了。
命令流文件可以定义为任何文件名，但是最好设定其扩展名为“.dat”，以防止和flac3d其它类型的文件相混淆。要读入命令流文件可以使用以下命令：call c:\我的文件夹\ file.dat
想知道计算这个问题需要占用多少的内存，可以输入以下命令：
print mem sav mem sys
输入命令后软件便会提供一张在现阶段内存的统计列表。我们这个例子在现阶段大概需要占用0.6MB的内存来存储数据，大约需要4.2MB的内存来进行模型的计算。
粘聚力和抗拉强度值取得很大，这样取值的目的是为了模拟得到模型在自重作用下的初始应力状态*。取值很大是防止了模型在初始加载分析时就达到塑性极限。一旦你拥有了很多的模拟经验，你就会知道这样做的原因是很明显的。
plot create GravV
plot set plane dip=90 dd=0 origin=3,4,0
plot set rot 15 0 20
plot set center 2.5 4.2 4.0
plot add bound behind
plot add bcont szz plane
plot add axes
plot show
注意以上命令均在“FLAC3D>” 提示符状态下输入，如果先输入plot命令，进入“plot”相关的提示符状态，那么命令中所有的plot均可删除。执行以上命令后，便会出现一个名叫GravV的视窗，并且系统指定它为当前显示窗口。我们在命令中已经定义了当前视图的剖面：倾角为90度（dip=90），倾向为0度（dd=0），通过点（3，4，0）。在命令中我们还定了视角（plot set rot 15 0 20 ；plot set center2.5 4.2 4.0）并且使用bound（boundary）命令来增加边界线框，线框只在平面的后面显示（由关键字behind定义，注可简写为BE）。plot add bcont szz plane命令的作用是在当前平面中绘出垂直应力分布图，图形如图2.5所示。
ID NUMBER——编号（id号），模型中的各个元素都有一个数字编号，同类元素之间就是以这个编号加以区分和识别的。下面这些元素都具有ID标识号：接触面（interface）、节点（gp）、单元（zone）、体(volume)、参考点（reference points）、历史记录项（history）、表（table）、绘图项（plot items）以及结构单元（包括：梁、索、桩、壳、衬砌、土
工栅格等）可以用PRINT命令来得到元素的ID号。用户可以自行定义当前输入的接触面、参考点、历史记录项、表、绘图项的id号。每个结构单元中还有CID编号，CID号用来标识该结构体中各个元素，结构体中的节点，单元，以及之间联系部件等都有各自的CID号。
LARGE STRAIN/SMALL STRAIN——大变形/小变形，在软件中，系统默认为小变形模式，即：系统中的节点坐标始终不发生变化，即使计算所得的应变比较大（与一般的单元尺寸相比）。在大变形模式中，可能发生几何非线性的变化。
节点和单元也可通过ID号加以识别。ID号对应的节点全局坐标，可以通过如下命令得到：print gp position
ID号对应的单元重心全局坐标可以通过如下命令得到：print zone
如果我们使用了solve命令，默认情况下，系统通过最大不平衡力比率来控制计算过程，当最大不平衡力与最初施加的节点力的平均值之比小于1×10-5，计算就会停止。我们也可以在输入solve命令之前通过输入以下命令人为的设置这个比率：set mech ratio = f 此处，f是一个用户给定的比率限定值。
在FLAC3D计算过程中的任何时候，我们可以按下键盘上的“ESC”键，来暂停计算。这经常会在step命令求解的情况下用到，尤其设定的求解步数非常多。我们可以按下“ESC”键，来暂停计算，然后查看历史记录项，再可以恢复计算，直到计算达到平衡状态。
日志文件
如果用户输入set logo on命令，那么在问题的分析过程中，软件会生成一个
“FLAC3D.log”文件。它是一个格式化的ASCII文件，默认的名为“FLAC3D.log”，
用户也可以使用set logfile +文件名命令来指定所生成的log文件的名字。当使用
了set logo on命令后，接下来在软件分析过程中在屏幕上出现的所有文本信息都
会被保存在日志文件中。
、历史记录文件
“FLAC3D.his”文件是用户输入history write n命令创建的一个历史记录文
件。命令中的n指的是历史纪录项的记录序号（参见命令手册第一章中有关history
命令的说明）。它也是一个格式化的ASCII 文件。系统默认的文件名为
“FLAC3D.his”。用户也可以通过在history write命令的后面加关上键词file 文
件名.his的方法来指定历史纪录文件的文件名。用户可以在文件中看到有关历史
记录项随计算步变化的值。
plot sel recover st zz cid 100 inter 0.56e6
显示初衬szz
plot sel recover st xx cid 100 inter 0.56e6
显示初衬sxx
print sel recover st zz cid 100 inter 0.56e6
print sel recover st xx cid 100 inter 0.56e6
plot sel recover pstress 1 ;最大主应力
显示应力数值
print sel recover stress zz
print sel recover stress zz
history sel cablesel force x y z/cid
history sel cablesel stress x y z/cid
围岩应力图
plot add contour smax outline on
plot contour zdisp inter 2.5e-4
hist sel cablesel force 0 4.425 22.4272
plot add contour smax range group 8 any group 9 &
any group 10 any group 15 any group 16 any
plot contour smax inter 0.5e5
plot add contour smax range group 5 any group 6 any &
group 7 any group 13 any group 14 any
plot add block group
m n r g 18 any g 14 any g 15 any g 17 any g 16 any
plot set center -2.833 15.09 15.05 rot 20 0 30 mag 3.81
plot show
抓图读数据命令
;施工步骤抓图
plot add block group
;m n r g 18 any g 14 any g 15 any g 17 any g 16 any
m n r g 13
plot set center 0 15 0.86 rot 10 0 20 mag 4.2
plot show
plot set background white
;岩体应力抓图
plot add contour smin inter 2e5
plot show
plot set background white
;初支应力图 实体
plot add contour smin range group 9 any group 10 any &
group 11 any group 12 any
plot set center 0 15 0.86 rot 10 0 20 mag 4.2
plot set background white
plot show
;沉降
plot add con zdis
plot set background white
plot show