Flac3D实例分析

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 之间的转换关系如下： =3(1-2v)E K G=2(1+v) E 16. 一般而言，大多数问题可以采用FLAC 3D 默认的收敛标准（或称相对收敛标准），即当体 系最大不平衡力与典型内力的比率R 小于定值10-5；（也可由用户自定义该值，命令：

FLAC3D 实例命令流1

第1部分命令流按照顺序进行 2-1定义一个FISH函数 new def abc abc = 25 * 3 + 5 End print abc 2-2使用一个变量 new def abc hh = 25 abc = hh * 3 + 5 End Print hh Print abc 2-3对变量和函数的理解 new def abc hh = 25 abc = hh * 3 + 5 End set abc=0 hh=0 print hh print abc print hh new def abc abc = hh * 3 + 5 end set hh=25 print abc set abc=0 hh=0 print hh print abc print hh 2-4获取变量的历史记录 new gen zone brick size 1 2 1 model mohr prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10 fix x y z range y -0.1 0.1 apply yvel -1e-5 range y 1.9 2.1 plot set rotation 0 0 45

plot block group def get_ad ad1 = gp_near(0,2,0) ad2 = gp_near(1,2,0) ad3 = gp_near(0,2,1) ad4 = gp_near(1,2,1) end get_ad def load load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) end hist load hist gp ydis 0,2,0 step 1000 plot his 1 vs -2 2-5用FISH函数计算体积模量和剪砌模量 new def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive print b_mod print s_mod 2-6 在FLAC输入中使用符号变量 New def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive gen zone brick size 2,2,2 model elastic prop bulk=b_mod shear=s_mod print zone prop bulk print zone prop shear 2-7 控制循环 New def xxx sum = 0 prod = 1 loop n (1,10)

FLAC3D 实例命令流1

第1部分命令流按照顺序进行2-1定义一个FISH函数 new def abc abc = 25 * 3 + 5 End print abc 2-2使用一个变量 new def abc hh = 25 abc = hh * 3 + 5 End Print hh Print abc 2-3对变量和函数的理解 new def abc hh = 25 abc = hh * 3 + 5 End set abc=0 hh=0 print hh print abc print hh new def abc abc = hh * 3 + 5 end set hh=25 print abc set abc=0 hh=0 print hh print abc print hh 2-4获取变量的历史记录 new gen zone brick size 1 2 1 model mohr prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10

fix x y z range y -0.1 0.1 apply yvel -1e-5 range y 1.9 2.1 plot set rotation 0 0 45 plot block group def get_ad ad1 = gp_near(0,2,0) ad2 = gp_near(1,2,0) ad3 = gp_near(0,2,1) ad4 = gp_near(1,2,1) end get_ad def load load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) end hist load hist gp ydis 0,2,0 step 1000 plot his 1 vs -2 2-5用FISH函数计算体积模量和剪砌模量 new def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive print b_mod print s_mod 2-6 在FLAC输入中使用符号变量 New def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive gen zone brick size 2,2,2 model elastic prop bulk=b_mod shear=s_mod print zone prop bulk print zone prop shear

flac3D基础与工程实例命令

2-1 new gen zon bri size 3 3 3 model elas prop bulk 3e8 shear 1e8 ini dens 2000 fix z ran z -.1 .1 fix x ran x -.1 .1 fix x ran x 2.9 3.1 fix y ran y -.1 .1 fix y ran y 2.9 3.1 set grav 0 0 -10 solve app nstress -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3 hist gp vel 0 3 3 plo hist 1 red plo add hist 2 blue solve 3-1 ;-------------------------------------工程信息 ;Project Record Tree export ;Title:Simple test ;---------------------------------计算第一步 ;... STATE: STATE1 .... config grid 10,10 model elastic group 'User:Soil' notnull model elastic notnull group 'User:Soil' prop density=1500.0 bulk=3E6 shear=1E6 notnull group 'User:Soil' fix x y j 1 fix x i 1 fix x i 11 set gravity=9.81

history 999 unbalanced solve save state1.sav ;----------------------------------计算第二步 ;... STATE: STATE2 .... initial xdisp 0 ydisp 0 initial xvel 0 yvel 0 model null i 4 7 j 8 10 group 'null' i 4 7 j 8 10 group delete 'null' history 1 xdisp i=4, j=11 solve save state2.sav ;--------------------------------绘图命令 ;*** plot commands **** ;plot name: syy plot hold grid syy fill ;plot name: Unbalanced force plot hold history 999 ;plot name: grid plot hold grid magnify 20.0 lred grid displacement ;plot name: Xdis-A plot hold history 1 line

FLAC3D命令流(整理版)

1、怎样查看模型？ 答：plot grid 可以查看网格，plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件？ 答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2

FLAC3D快速入门及简单实例

FLAC3D快速入门 及简单实例 李佳宇编 LJY指南针教程

前言 FLAC及FLAC3D是由国际著名学者、英国皇家工程院院士、离散元的 发明人Peter Cundall博士在70年代中期开始研究的，主要面对岩土工程的通 用软件系统，目前已经在全球70多个国家得到广泛应用，在岩土工程学术界 和工业界赢得了广泛的赞誉。前国际岩石力学会主席 C.Fairhurst(1994)对 FLAC程序的评价是：“现在它是国际上广泛应用的可靠程序。” 我从研二(2010年)开始接触FLAC3D，最初的原因是导师要求每一个人至 少学会一个数值计算软件，而他嘴里每天念叨最多的就是FLAC，自己当时对数 值计算一无所知，便答应老师要学会FLAC3D。第一次打开软件界面，我心里 就凉了大半截，面对着一个操作界面跟记事本无异的所谓“功能强大”的岩土工 程专业软件，半点兴趣也提不起来。年底，从项目工地回到学校准备论文开题， 老师对我的开题报告非常不满意，当着全教研室师生的面，劈头盖脸大批一顿， 第二天又找谈话。在巨大的压力和强烈的自尊心驱使下，我硬着头皮开始啃 FLAC3D，一个半月之后，终于有了初步的计算结果，对老师有个交代，我也能 回家过年了。 前面这一段过程可能是大多数FLAC3D初学者的必经阶段，或者是即将 开始软件学习的人惧怕的事情。毫无疑问，FLAC3D极其不友好的界面是阻碍 初学者前进的很大障碍，当然还包括它是一个全英文的软件。但是当你费尽周折 的走进FLAC3D的世界，你就会发现它独特的魅力，比如简洁的界面，快捷的 命令流操作，高效的计算方法，不易报错等等。另外一个拿不上台面的优点就是 它非常小巧，包括Manual在内一共才几十兆大小，而且已经被破解成绿色版， 只要把它和命令流装进U盘，你就可以随便找一个身边功能最强大的电脑开始 计算了，如果你有过ANSYS、ABAQUS等大型软件痛苦的安装经历，你便能 毕业之后，本以为不用再接触数值计算，但工作需要使得我又一次开始与理解“绿色版”的含义，当然还请大家尊重知识产权，支持正版。 FLAC3D进行亲密接触。我的领导给了我很多新思路和很大的支持，如今我的 水平比研究生时有了不少提高。于是，我想把我的经验总结成文，希望对初学 者起

陈育民《FLAC3D基础与工程实例》全部命令流

《FLAC/FLAC 3D基础与工程实例》全部命令流 1 gen zone bri p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 & p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 & size 10 10 10 rat 1.0 0.9 1.1 group brick_1 gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 & size 10 10 10 rat 1.0 0.9 1.1 group brick_2 gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0.9 1.1 group brick_3 plot sur 2-1 new gen zon bri size 3 3 3 model elas prop bulk 3e8 shear 1e8 ini dens 2000 fix z ran z -.1 .1 fix x ran x -.1 .1 fix x ran x 2.9 3.1 fix y ran y -.1 .1 fix y ran y 2.9 3.1 set grav 0 0 -10 solve app nstress -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3 hist gp vel 0 3 3 plo hist 1 red plo add hist 2 blue solve 3-1 ;-------------------------------------工程信息 ;Project Record Tree export ;Title:Simple test ;---------------------------------计算第一步 ;... STATE: STATE1 .... config grid 10,10 model elastic group 'User:Soil' notnull model elastic notnull group 'User:Soil' prop density=1500.0 bulk=3E6 shear=1E6 notnull group 'User:Soil' fix x y j 1 fix x i 1 fix x i 11 set gravity=9.81 history 999 unbalanced solve save state1.sav ;----------------------------------计算第二步 ;... STATE: STATE2 .... initial xdisp 0 ydisp 0 initial xvel 0 yvel 0 model null i 4 7 j 8 10 group 'null' i 4 7 j 8 10 group delete 'null' history 1 xdisp i=4, j=11 solve save state2.sav ;--------------------------------绘图命令 ;*** plot commands **** ;plot name: syy plot hold grid syy fill ;plot name: Unbalanced force plot hold history 999 ;plot name: grid plot hold grid magnify 20.0 lred grid displacement ;plot name: Xdis-A plot hold history 1 line 5-1 new ; =============================== ; 定义球体半径和半径方向上单元网格数 ; =============================== def parm rad=10.0 rad_size=5 end parm ; =============================== ; 建立八分之一球体外接立方体网格

FLAC原理实例与应用指南FLAC3D实例分析教程

《FLAC 原理实例与应用指南》FLAC3D 实例分析教程 1 FLAC3D 实例分析教程 刘波韩彦辉（美国）编著《FLAC 原理实例与应用指南》 北京：人民交通出版社，2005.9 Appendix(附录) 版权所有：人民交通出版社，Itasca Consulting Group, USA 说明： 1．本实例分析教程是为方便读者学习、应用FLAC 和 FLAC3D 而编写的，作为《FLAC 原理、实例与应用指南》一书的 附录。 2．计算算例参考了Itasca Consulting Group 的培训算例，命 令流的解析旨在方便读者理解FLAC 和FLAC3D 建模及求解问题的一般原则与步骤。 3．实例分析的算例中，FLAC 算例是基于FLAC 5.00 版本、FLAC 3D 算例是基于FLAC 3D 3.00 版本实现计算分析的。读者在学习和研究相关算例时，请务必采用Itasca 授权的合法版本进行分析计算。 4．本实例分析教程仅供读者参考，读者在参考本教程算例 进行工程分析时，编者对可能产生的任何问题概不负责。 编者 《FLAC 原理实例与应用指南》FLAC3D 实例分析教程

2 实例分析1：基坑开挖 图1：基坑开挖的位移等值线图 例1： ; simple tutorial – trench excavation ;简单的隧道开挖方法 ; Step 1: create initial model state ;第一步：初始模型的建立 ; create grid ; 建立网格 gen zone brick size 6 8 8 ;建立矩形的网格区域，其大小为6×8×8 ; assign material ;指定材料的性质 model mohr ;采用摩尔-库仑模型 prop bulk 1e8 shear 0.3e8 ;模型的材料性质：体积模量1×108 剪切模量0.3×108 prop fric 35 coh 1e3 tens 1e3 ;摩擦角35°粘聚力1×103 抗拉强度1×103 ; set global conditions ;设置全局参数 set grav 0,0,-9.81 《FLAC 原理实例与应用指南》FLAC3D 实例分析教程 3 ;设置重力加速度，z 坐标正方向为正，故为-9.81 ini dens 2000 ;初始密度为2000 ; set boundary conditions ;设置初始边界情况 fix x range x -0.1 0.1 ;在x 方向上固定x = 0 边界，为滚动支撑 fix x range x 5.9 6.1 ;在x 方向上固定x = 6 边界，为滚动支撑 fix y range y -0.1 0.1 ;在y 方向上固定y = 0 边界，为滚动支撑 fix y range y 7.9 8.1 ;在y 方向上固定y = 8 边界，为滚动支撑 fix z range z -0.1 0.1

FLAC3D模拟实例循环开挖与支护

FLAC3D模拟实例循环开挖与支护n res ini.sav set geometry=0.001 ini ydis0 ini xdis0 ini zdis0 ini yvel0 ini xvel0 ini zvel0 m m prop bulk 4.0e9shear 2.5e9fri32coh 2.0e6& range grou diban-shayan ; prop bulk 1.8e9shear 1.2e9fri25coh 1.0e6& range grou diban-niyan any grou hangdao any ; prop bulk 1.2e9shear0.8e9fri22coh0.8e6& range grou diban-gentuyan ; prop bulk 1.9e9shear 1.3e9fri24coh 1.0e6& range grou diban-tniyan

; prop bulk0.7e9shear0.8e9fri21coh0.7e6& range grou mc any grou gzm any ; prop bulk 3.0e9shear 2.5e9fri30coh 1.8e6& range grou dingban-fenshayan ; prop bulk 1.5e9shear 1.2e9fri25coh 1.1e6& range grou dingban-niyan ; prop bulk 3.5e9shear 2.5e9fri34coh 1.4e6& range grou dingban-shayan ;添加接触面 gen separate gzm interface1wrap mc gzm interface2wrap dingban-fenshayan gzm interface1prop kn20e9ks10e9tens1e9 interface2prop kn20e9ks10e9tens1e9 set mech ratio=5e-4 def excav_mc loop n(excav_p,excav_p_z+cut_liang)

《FLAC3D基础与工程实例》全部命令流

《FLAC/FLAC 3D基础与工程实例》全部命令流1 gen zone bri p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 & p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 & size 10 10 10 rat 1.0 0.9 1.1 group brick_1 gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 & size 10 10 10 rat 1.0 0.9 1.1 group brick_2 gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0.9 1.1 group brick_3 plot sur

new gen zon bri size 3 3 3 model elas prop bulk 3e8 shear 1e8 ini dens 2000 fix z ran z -.1 .1 fix x ran x -.1 .1 fix x ran x 2.9 3.1 fix y ran y -.1 .1 fix y ran y 2.9 3.1 set grav 0 0 -10 solve app nstress -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3 hist gp vel 0 3 3 plo hist 1 red plo add hist 2 blue solve

;-------------------------------------工程信息 ;Project Record Tree export ;Title:Simple test ;---------------------------------计算第一步 ;... STATE: STATE1 .... config grid 10,10 model elastic group 'User:Soil' notnull model elastic notnull group 'User:Soil' prop density=1500.0 bulk=3E6 shear=1E6 notnull group 'User:Soil' fix x y j 1 fix x i 1 fix x i 11 set gravity=9.81 history 999 unbalanced solve save state1.sav ;----------------------------------计算第二步 ;... STATE: STATE2 .... initial xdisp 0 ydisp 0 initial xvel 0 yvel 0 model null i 4 7 j 8 10 group 'null' i 4 7 j 8 10 group delete 'null' history 1 xdisp i=4, j=11 solve save state2.sav ;--------------------------------绘图命令 ;*** plot commands **** ;plot name: syy plot hold grid syy fill ;plot name: Unbalanced force plot hold history 999 ;plot name: grid plot hold grid magnify 20.0 lred grid displacement ;plot name: Xdis-A plot hold history 1 line

FLAC及FLAC3D基础与工程实例51-3

FLAC及FLAC3D基础与工程实例51 1.2.1FLAC/FLAC3D的使用特征；FLAC/FLAC3D的使用特征主要表现为：；?命令驱动模式；FLAC/FLAC3D有两种输入模式：①人机交互；FLAC/FLAC3D专为岩土工程力学分析开发，；FLAC/FLAC3D可以模拟多种结构形式，如岩；借助其强大的绘图功能，用户能绘制各种图形和表格；?开放性；FLAC/FLAC3D几乎是一个全开放的系统， 1.2.1 FLAC/FLAC3D的使用特征 FLAC/FLAC3D的使用特征主要表现为： ? 命令驱动模式 FLAC/FLAC3D有两种输入模式：① 人机交互模式，即从键盘输入各种命令控制软件的运行；② 命令驱动模式，即写成命令流文件，由文件来控制软件的运行。其中，命令驱动模式为FLAC/FLAC3D的主要输入模式，尽管这种驱动方式对于简单问题的分析过于繁杂，对软件初学者而言相对较困难，但对于那些从事大型复杂工程问题分析而言，因涉及多次参数修改、命令调试，这种方式无疑是最有效、最经济的（当然，由于二维建模相对简单，照顾不少用户的使用习惯，在FLAC中也可以采用界面操作模式即GIIC模式进行分析计算）。 ? 专一性 FLAC/FLAC3D专为岩土工程力学分析开发，内置丰富的弹、塑性材料本构模型（其中FLAC内置11个，FLAC3D内置12个），有静力、动力、蠕变、渗流、温度5种计算模式，各种模式间可以相互耦合，以模拟各种复杂的工程力学行为。 FLAC/FLAC3D可以模拟多种结构形式，如岩体、土体或其它材料实体：梁、锚元、桩、壳以及人工结构，如支护、衬砌、锚索、岩栓、土工织物、摩擦桩、板桩等。通过设置界面单元，可以模拟节理、断层或虚拟的物理边界等。 借助其强大的绘图功能，用户能绘制各种图形和表格。用户可以通过绘制计算时步函数关系曲线来分析、判断体系何时达到平衡与破坏状态，

FLac3D接触面简单实例

有接触面模型 gen zone brick p0 0 0 0 p1 5 0 0 p2 0 5 0 p3 0 0 2 size 10 10 4 group 1 gen zone brick p0 0 0 2 p1 5 0 2 p2 0 5 2 p3 0 0 4 size 10 10 4 group 2 gen zone brick p0 0 0 4 p1 5 0 4 p2 0 5 4 p3 0 0 6 size 10 10 4 group 3 gen separate 2 int 1 wrap 1 2 gen separate 3 int 2 wrap 2 3 plot block group int red fix z range z -0.1 0.1 fix x range x -0.1 0.1 fix x range x 4.9 5.1 fix y range y -0.1 0.1 fix y range y 4.9 5.1 model mohr prop bulk 2.4e9 shear 1.5e9 coh 2.5e6 fri 32 tension 1.2e6 ran group 1 prop bulk 2.4e9 shear 1.5e9 coh 2.5e6 fri 32 tension 1.2e6 ran group 2 prop bulk 2.4e9 shear 1.5e9 coh 2.5e6 fri 32 tension 1.2e6 ran group 3 ini dens 2660 ran group 1 ini dens 2660 ran group 2 ini dens 2660 ran group 3 inter 1 prop kn 0.88e11 ks 0.88e11 fri 15 coh 2e6 ten 1e6 inter 2 prop kn 0.88e11 ks 0.88e11 fri 15 coh 2e6 ten 1e6 set grav 0 0 -10 ini szz -7e6 range z 0 6 ini sxx -8.4e6 range z 0 6 ini syy -8.4e6 range z 0 6 app nstress -7e6 range z 5.9 6.1 x 0 5 y 0 5 plot con zd plot add int black solve 无接触面对比模型 gen zone brick p0 0 0 0 p1 5 0 0 p2 0 5 0 p3 0 0 2 size 10 10 4 group 1 gen zone brick p0 0 0 2 p1 5 0 2 p2 0 5 2 p3 0 0 4 size 10 10 4 group 2 gen zone brick p0 0 0 4 p1 5 0 4 p2 0 5 4 p3 0 0 6 size 10 10 4 group 3 plot block group int red fix z range z -0.1 0.1 fix x range x -0.1 0.1 fix x range x 4.9 5.1 fix y range y -0.1 0.1 fix y range y 4.9 5.1 model mohr

FLAC3D条带开采案例(代码)

n gen zone brick p0 0 0 0 p1 2500 0 0 p2 0 2000 0 p3 0 0 56 size 100 80 2 group 煤层底板; gen zone brick p0 0 0 56 p1 2500 0 56 p2 0 2000 56 p3 0 0 60 size 100 80 1 group 煤层; gen zone brick p0 0 0 60 p1 2500 0 60 p2 0 2000 60 p3 0 0 76 size 100 80 1 group 直接顶;顶板gen zone brick p0 0 0 76 p1 2500 0 76 p2 0 2000 76 p3 0 0 100 size 100 80 1 group 老顶;顶板gen zone brick p0 0 0 100 p1 2500 0 100 p2 0 2000 100 p3 0 0 280 size 50 40 4 group 中粒、粉砂岩互层;顶板 gen zone brick p0 0 0 280 p1 2500 0 280 p2 0 2000 280 p3 0 0 480 size 50 40 4 group 泥岩、砂泥岩互层;顶板 gen zone brick p0 0 0 480 p1 2500 0 480 p2 0 2000 480 p3 0 0 640 size 50 40 4 group 表土层;顶板 modelmohr prop bulk 6.94e9 shear 5.21e9 coh 24.6e7 fric 38 tens 3e7 range group 煤层底板 prop bulk 6.25e9 shear 0.35e9 coh 1e7 fric 25 tens 1e7 range group 煤层 prop bulk 1.33e9 shear 0.44e9 coh 1.5e7 fric 32 tens 1e7 range group 直接顶 prop bulk 7.44e9 shear 5.12e9 coh 8.56e7 fric 37 tens 3e7 range group 老顶 prop bulk 8.26e9 shear 4.72e9 coh 6.25e7 fric 33 tens 2.6e7 range group 中粒、粉砂岩互层prop bulk 8.83e9 shear 1.89e9 coh 1.5e7 fric 33 tens 1e7 range group 泥岩、砂泥岩互层prop bulk 0.017e9 shear 0.0036e9 coh 0.1e7 fric 20 tens 1e7 range group 表土层 initial density 2.6e3 range group 煤层底板 initial density 1.4e3 range group 煤层 initial density 2.4e3 range group 直接顶 initial density 2.5e3 range group 老顶 initial density 2.64e3 range group 中粒、粉砂岩互层 initial density 2.6e3 range group 泥岩、砂泥岩互层 initial density 1.8e3 range group 表土层 fix z range z -0.1 0.1 fix x range x -0.1 0.1 fix x range x 2499.9 2500.1 fix y range y -0.1 0.1 fix y range y 1999.9 2000.1 setgrav 0 0 -9.8 histunbal plotcontzdis

FLAC3D 实例命令流5

第5部分命令流按照顺序进行 new gen zon radcyl p0 9 0 0 p1 18 0 0 p2 9 36 0 p3 9 0 9 & size 4 24 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outr gen zone cshell p0 9 0 0 p1 12 0 0 p2 9 36 0 p3 9 0 3 & size 1 24 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group liny fill group inr gen zon reflect dip 90 dd 90 orig 9 0 0 gen zon reflect dip 0 dd 0 ori 9 0 0 ;1/2隧道地层部分 gen zon brick p0 18 0 -9 p1 36 0 -9 p2 18 36 -9 p3 18 0 9 size 10 12 8 group outr1 gen zon brick p0 18 0 9 p1 36 0 9 p2 18 36 9 p3 18 0 13 size 10 12 2 group outr2 gen zon brick p0 0 0 9 p1 18 0 9 p2 0 36 9 p3 0 0 13 size 8 12 2 group outr2 gen zon brick p0 18 0 -23 p1 36 0 -23 p2 18 36 -23 p3 18 0 -9 size 10 12 8 group outr3 gen zon brick p0 0 0 -23 p1 18 0 -23 p2 0 36 -23 p3 0 0 -9 size 8 12 8 group outr3 gen zon reflect dip 90 dd 270 orig 0 0 0 group linz range group liny x -12.1 -5.99 ;绘制模型图 plot add axes red plot set rotation 20 0 45 plot block group ;0-自重应力场求解 set gravity 0 0 -10 ;给定边界条件 fix z range z -23.01,-22.99 fix x range x -36.01,-35.99 fix x range x 35.99,36.01 fix y range y -0.01 0.01 fix y range y 35.99,36.01 model mohr ;围岩的密度 ini density 1800 ;体积、剪切、摩擦角、凝聚力、抗拉强度 prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4 step 3000 save Gsol.sav ;第1步 rest Gsol.sav initial xdisp=0 ydisp=0 zdisp=0 model null range group inrz y 0 6