flac学习吐血体会

第一天看版主的初学者建议！
初学者将FLAC/FLAC3D版所有帖子都阅读一遍；学有所得后，再阅读一遍，会发觉又有新的收成。

Flac3D中为岩土工程问题的求解开发了特有的本构模型，总共包括了10种材料模型：
1.开挖模型]
2.3个弹性模型（各向同性，横观各向同性和正交各向同性弹性模型）
3.6个塑性模型（Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍及节理模型、双线性应变硬化/软化遍及节理模型和修正的cam粘土模型）。

因为我是学采矿专业的，我就先学习一下开挖模型吧！
第二天学习建模
FLAC3D建模的一样步骤如下：
一、建模
new
ge zo b po x y z &
p1 x y z &# Z/ E7 g/ L( F) B% A
p2 x y z &7 F. m$ ^ U- g5 h' S7 u
p3 x y z &1 i* d1 K% C; c9 y6 R; a
size x y z
plot
Show
add surface color
add axes color
如此一个立体模型就成立出来了！
二、概念材料参数
mo moh/null/...
pro bu=~ sh=~ fric=~ c=~ ten=~
其中bulk是指体积模量也确实是K=E/3(1-2v),E是弹模v是泊松比，当v接近0.5时，不能盲目利用；
shear是指剪切模量也确实是G=E/2(1+v）；
cohesion是指内聚力；
tension是指抗拉强度。

三、加载及边界
set g 0 0 -10
in de=~
fix x/y/z ra x/y/z ~ ~
其中要有一个面为自由面。

四、求解7 R( ^: m9 s% n C. f
h n=10
h unb
h gp x/y/zdisp x y z
set mec fo 30
so
此为显式求解进程。

五、结果分析
监控变量：
print h
plot
history 1
history 2
位移等值线图：
plot
cr ~
add con disp/...
add axes color
show
应力等值线图：
Clear
add bcon szz/...
add axes
剖面上的等值线图：
cr ~
set p dip=~ dd=~ o=x,y,z
add bo be
add bcon szz/... P
add axes
show
有这些命令能够查看你想分析的FLAC图
七、开挖
pro c=~ ten=~
mo null ra x=~,~ y=~,~ z=~,~
set large/small1
in xdisp=0 ydisp=0 zdisp=0
系统中的所有节点位移全数清零。

八、求解与后处置
step~
plot cr dispcont
plot copy ~ dispcont settings
plot add con disp p be shade on
plot add axes
plot show
开挖模型终于学习完了，通过学习我能够成立一个简单的开挖模型并计算、分析了。

下一步的工作：我想成立一个曲面的立体模型来反映矿山地形地貌。

研究中......
通过学习，我能够成立一个分层的地质模型，如何把地表设置成曲面呢？形成矿体上方的地形地貌!
一个分层的地质模型：
命令流如下：
new
gen zone b size 40 10 1 p0 0 0 0 p1 400 0 0 p2 0 200 0 p3 0 0 1.21 group 泥岩1
gen zone b size 40 10 2 p0 0 0 1.21 p1 400 0 1.21 p2 0 200 1.21 p3 0 0 7.07 group 灰岩1
gen zone b size 40 10 3 p0 0 0 7.07 p1 400 0 7.07 p2 0 200 7.07 p3 0 0 13.47 group 粉砂岩1
gen zone b size 40 10 1 p0 0 0 13.47 p1 400 0 13.47 p2 0 200 13.47 p3 0 0 14.96 group 灰岩2
gen zone b size 40 10 4 p0 0 0 14.96 p1 400 0 14.96 p2 0 200 14.96 p3 0 0 31.4
group
粉砂岩2
gen zone b size 40 10 1 p0 0 0 31.4 p1 400 0 31.4 p2 0 200 31.4 p3 0 0 32.46 group 泥岩2
gen zone b size 40 10 3 p0 0 0 32.46 p1 400 0 32.46 p2 0 200 32.46 p3 0 0 39.73 group 粉砂岩3
gen zone b size 40 10 6 p0 0 0 39.73 p1 400 0 39.73 p2 0 200 39.73 p3 0 0 51.93 group 细砂岩1
gen zone b size 40 10 7 p0 0 0 51.93 p1 400 0 51.93 p2 0 200 51.93 p3 0 0 65 group
粉砂岩
gen zone b size 40 10 3 p0 0 0 65 p1 400 0 65 p2 0 200 65 p3 0 0 67.93 group 10煤
gen zone b size 40 10 4 p0 0 0 67.93 p1 400 0 67.93 p2 0 200 67.93 p3 0 0 71.97 group 泥岩
gen zone b size 40 10 1 p0 0 0 71.97 p1 400 0 71.97 p2 0 200 71.97 p3 0 0 72.35 group 天然焦
gen zone b size 40 10 2 p0 0 0 72.35 p1 400 0 72.35 p2 0 200 72.35 p3 0 0 75.15 group 岩浆岩
gen zone b size 40 10 2 p0 0 0 75.15 p1 400 0 75.15 p2 0 200 75.15 p3 0 0 78.85 group 细砂岩2
gen zone b size 40 10 25 p0 0 0 78.85 p1 400 0 78.85 p2 0 200 78.85 p3 0 0 132.6 group
泥岩4.
gen zone b size 40 10 4 p0 0 0 132.6 p1 400 0 132.6 p2 0 200 132.6 p3 0 0 136.79 group 细砂岩3
gen zone b size 40 10 3 p0 0 0 136.79 p1 400 0 136.79 p2 0 200 136.79 p3 0 0 140.24 group 泥岩5,
gen zone b size 40 10 3 p0 0 0 140.24 p1 400 0 140.24 p2 0 200 140.24 p3 0 0 146.29 group 粉砂岩5
gen zone b size 40 10 4 p0 0 0 146.29 p1 400 0 146.29 p2 0 200 146.29 p3 0 0 150 group 泥岩6
plot
Show
add surface green
add axes red
1 FLAC建模方式
1.1 建模
（1）设计计算模型的尺寸
（2）计划计算网格数量和散布
（3）安排工程对象（开挖、支护等）
（4）给出材料的力学参数
（5）确信边界条件
（6）计算模拟
1.2 网格生成：Grid i,j 例如：grid 30,20
1.3 网格计划：Gen x1,y1 x2,y2 x3,y3 x4,y4 例如：Gen 0,0 0,10 10,20 20,0 1.4 分区计划网格。

例如：Gen xI1,yI1 xI2,yI2 xI3,yI3 xI4,yI4 i=1,10 j=1,21 （I区）
Gen xII1,yII1 xII2,yII2 xII3,yII3 xII4,yII4 i=10,20 j=1,21 (II区）
1.5 特殊形状的网格
（1）圆形gen circle xc,yc rad
（2）弧线gen arc xc,yc xb,yb theta
（3）直线gen line x1,y1 x2,y2
（4）任意形状tab 1 x1,y1, x2,y2, ¼,xn,yn, x1,y1 gen tab 1
1.6 赋给单元材料性质
mod e (弹性)
prop d 1800e-6 bu 12.5 sh 5.77 i=1,20 j=1,10
prop d 2400e-6 bu 1250 sh 577 i=1,20 j=11,20
mod m （弹塑性Mohr-Coulumb准那么）
prop d 1800e-6 bu 12.5 sh 5.77 c 0 fri 20 ten 0.015 reg i,j
1.7 赋给模型边界条件
（1）固定边界（结点）Fix x i=1, j=1,21 Fix y i=1,21 j=1
（2）施加边界力(结点) apply yf -10 i=1,21 j=21
或apply syy -10 i=1,21 j=21
或apply xf -5 i=21, j=1,21
或apply sxx -5 i=21, j=1,21
（3）赋单元内应力（单元）
ini sxx -10 i=1,20 j=1,20
ini syy -5 var 0 4 i=1,21 j=1,21.
1.8 计算
Set grav 9.81
Set large
Step 1000
1.9 结果显示
Plot grid 显示网格
Plot bo 显示边界
Plot plas 显示塑性区
Plot sig1 fi 显示最大主应力s1
Plot sig2 fi 显示最小主应力s2
Plot sdif fi 显示主应力差(s1- s2))
Plot str 显示主应力矢量场
Plot xdis fi 显示X方向位移
Plot ydis fi 显示Y方向位移
Plot disp 显示位移矢量场
1.10 保留与挪用结果命令
Call test.txt (或ca test.dat) 挪用数据
Save test.sav 保留结果
New 从头开始
Rest test.sav 挪用结果Quit 退出程序
2.FLAC运用技术
2.1 模型尺寸
2.2 模拟开挖Mod nu i=6,15 j=5,12 (或region i,j) 2.3 模拟锚杆支护（端锚）
struct cable begin grid i1,j1 end grid i2,j2 seg n prop 1 stru prop 1 e 2e5 yield 0.5 a 0.235e-3 sbond 0 kbond 0.0001 stru prop 1 sfri 30 peri 0.2723 den 7.5e-3
2.4 模拟锚杆支护（全长锚固）
struct cable begin grid i,j end x,y seg n prop 2
stru prop 1 e 2e5 yield 0.5 a 0.235e-3 sbond 0.42 kbond 5.37
2.5 模拟锚杆支护（预应力锚固）
struct cable begin grid i1,j1 end i2,j2 seg n ten 10 prop 1 struct cable begin node n1 end x,y seg n prop 2
stru prop 1 e 2e5 yield 0.5 a 0.235e-3 sbond 0. kbond 0.1
s
stru prop 2 sfri 30 perimeter 0.2723 den 7.5e-3
2.6 模拟断层
Mod nu j=37
int 1 as from i1,j1 to ,i1,j2 bs from i2,j1 to i2,j2
int 1 kn 1000 ks 400 fri 25 c 0.01 t 0.0001
2.7 其它技术
（1）网格优化
（2）先弹性、后塑性
（3）分步开挖
（4）动力学问题模拟.
（5）固流耦合问题模拟
（6）合理说明结果
（7）多种形式输出结果
（8）报告格式
第三、四天学习彭文斌的《FLAC 3D有效教程》
第5天运用FLAC做一个小项目
第六天：分享一些FLAC3D命令的FAQ资料！
1.如何调整视图的大小、角度？,
答：综合利用x, y, z, m, Shift键，配合利用Ctrl+R，Ctrl+Z等快捷键
2.最先需要把握的命令有哪些？
答：需要把握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处置和求解的命令。

3. 如何看模型的样子？
答：plo blo gro能够看到不同的group的颜色散布
4.如何看模型的边界情形？
答：plo gpfix red sk
5.如何看模型的体力散布？
答：plo fap red sk
6.如何看模型的云图？
答：位移：plo con dis (xdis, ydis, zdis) 应力：plo con sz (sy, sx, sxy, syz, sxz)
7.如何看模型的矢量图？
答：plo dis (xdis, ydis, zdis
8.如何看模型有多少单元、节点？
答：plo info
9.如何输出模型的后处置图？
答：File/Print type/Jpg file，然后选择File/Print，将保留格式选择为jpg文件
10.如何挪用一个文件？
答：File/call或call命令
11.如何施加面力？!
答：app nstress
12.如何进行边界约束？
答：fix x ran （约束的是速度，在初始情形下约束等效于位移约束）
13.如何明白每一个单元的ID？
答：用鼠标双击单元的表面，能够明白单元的ID和坐标
14.如何进行切片？
答：plo set plane ori (点坐标) norm (法向矢量)
plo con sz plane (显示z方向应力的切片）
15. 如何保留计算结果？
答：save +文件名
16.如何挪用已保留的结果？
答：rest +文件名；或File/Restor
17.如何暂停计算？
答：Esc
18.如安在程序中进行暂停，并可恢复计算？"
答：在命令中加入pause命令，用continue进行继续
19.如何跳过某个计算步？
答：在计算中按空格键跳过本次计算，自动进入下一步
20.Fish是什么东西？
答：是FLAC3D的内置语言，能够用来进行参数化模型、完成命令本身不能进行的功能21.initial 与apply 有何区别？
答：initial初始化命令，如初始化计算体的应力状态等
apply边界条件限制命令，如施加边界的力、位移等约束等。

initial的应力状态会随计算进程的发生而发生改变，一样体力需要初始化，而apply 施加的边界条件可不能发生转变。

22. 如何查看各个时段不平稳力的具体数值？
答：采纳his来记录计算，包括位移应力等命令his unbal
his gp(zone) zdis range (0 0 0)
或id=？导出数据命令his write n vs m begin 时步end 时步file filename.his.
n表示纪录的id m表示时步
要导出不平稳力的具体数值
his unbal
step 100000 or solve7 k:
利用上述命令就能够够查看各个时步下的不平稳力的具体数值
23.如何挪用一些可选模块？
答：config dyn (fluid, creep, cppudm)
24 .如安在圆柱体周围如何施加约束条件？
答：能够用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2<r1，其实确实是选择range的时候选两个圆柱的差，即取得边界。

25.如何能把一个PLOT的图像数据导出来以便用其他软件画图？
答：用set log on 命令，把数据导出来，转到excel里处置一下，然后用surfer或什么作图软件绘制就好了。

26.关于dd、dip的明白得：dip-direction angle, measured in the global xy-plane clockwise from the positive y-axis 那个是从y轴到x顺时针方向所转过的角度；dip dip angle, measured in the negative z-direction from the global xy-plane 那个是xy平面与z轴负方向的夹角。

dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线，dd表示y轴正向顺时针到那条射线的夹角。

第一应该依照dd的方向大体确信那个面的朝向，dd 指的是从y轴正方向按顺时针（clockwise）方向转向所要确信面的法线方向在xy平面上的投影的夹角，然后再确信dip，dip指的是从xy平面转向所要确信的平面的角．
27.如何挪用FISH文件？
答：新手挪用FISH文件犯错的缘故：是：文件名为fishcall.fis.txt(因为在记事本里那个txt 的后缀是默许的）。

修改方式是：将fishcall.fis.txt里面的内容全数复制到新的记事本里，再按保留，保留时候将保留类型改成所有文件(*,*),这时你即可在文件名栏内输入fishcall.fis，类型也就正确了。

28.如何明白得Tension-p shear-p tension-n shear-n表示的含义？
答：在计算的循环里面,每一个循环中,每一个zone都依据failure criterion处于不同的状态,shear 和tension别离表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即别离表示在此刻和以前的循环当中处于塑性状态
29.绘制相关曲线命令的一些说明
1)在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负)
2)"-"表示对其值作"mirror",若是大伙儿利用AUTOCAD的话就明白"mirror"的含义了,比如上面的hist 1 gp ydisp的值是0~100,那么vs -1就变成-100~0.以此类推.
3)什么缘故要这么做呢?主若是为了符合岩土工程的适应需要,例如桩载荷实验曲线,沉降曲
线等都不是画在第一象限内.
30. 如何查看剖面上任一点的位移值？
print gp disp range x() y( ) z( )或!
print gp disp range id
plot gp disp range x() y( ) z( )或
plot gp disp range id
31.fix x y z range x -0.1 0.1 ..........................a
fix x range x -0.1 0.1................................b 请问a式与b式有何区别？
a 表示在x=0 那个平面是固定的
b 表示在x=0 那个平面只x方向不能发生运动
32.如何显示变形轮廓线的命令？plo ske magf 10 其中10为放大系数
2020年的第一天：
第一介绍一下FLAC3D4.0新版本的安装吧！
第一步：放DVD安装光盘到光驱；
第二步：打开光盘，安装点START安装程序，按提示一步一步的进行，建议利用32位的！第三步：插上密码狗；
第四步：就能够够为所欲为的利用了！
第一讲！断层的加入
断层用INTERFACE来建
思路：1：别离成立上、下盘
2：成立INTERFACE 并设置相关参数
3：移动到一块
4：重点是法向量（能够看看高中的数学讲义）
5：如此只能成立平面的断层
6：关于曲面的断层应该进行更深的研究
第二讲：露天高边坡的建模
能够用CYLINDER和CSHELL成立
思路：1：划分好台阶
2：台阶的成立很重要必然要矩形形式的CSHELL 不能用楔形形式的CSHELL成立台阶面很关键
3：一层一层的摞起来就形成了具有台阶的高边坡
4：边坡脚位置到左侧界为1.5坡高，边坡顶位置到右边界的距离为2.5倍坡高，坡顶到底面为2倍坡高
5：运用自身的功能成立的边坡比较规那么，而且工作量专门大
6：下一步的工作确实是想通过其他软件进建模来导入l
第三讲：地下采空区的稳固性分析
能够通过成立剖面图来进行分析研究
思路：1：选取剖面具体位置
2：用BRICK进行建模
3：按实际开采进程开挖吧
4：开挖完了确实是空区了
5：运用PL的那些相关命令进行分析研究
6：重要的不是模拟的进程，而是分析的重要性
第四讲：采煤工作面的推动进程中围岩压力分析研究
用BRICK建模
思路：1：成立模型，尤其要学会用GROUP那个命令很壮大
2：模拟开挖进程NULL命令的利用很关键
3：分析顶板和围岩的应力和位移转变
4：能分析出第一次来压和周期来压的大致范围
5：能够进行离层的分析模拟
6：关于工作面通过断层的数值模拟我会在下一专题介绍
第五讲:简化模型与符合现实
1.这是一个辩证的问题,鱼和熊掌不可兼得;
2.思路要明确,你要明白自己想分析什么东西;
3.预期一下可能会是个什么样的结果;
4.此刻随着 4.0的发布网格密一些是能够的,但也要尽可能的减少没必要要的网格设置;
5.FLAC3D只是一个工具，不要太依托模拟的结果，要和其他的方式结合利用；
6.就到这吧！这是自己最近做的几个课题以后的一些方式体会！希望能对开始接触项目研究的人员提供一些建议和参考！。