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快速入门及简单实例
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命令流(挺实用)
1、怎样查看模型？ 答：plot grid 可以查看网格，plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件？ 答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2 《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模拟实例循环开挖与支护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) 有接触面模型 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 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 第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 作一个盾构隧道开挖，地下水位设为z=12.6m，（隧道中心线设为z=0m）。分了7层土，给定初始条件的代码如下： config fluid （建模过程省） model fl_isotropic prop perm 1e-12 poro 0.48 range z 9.1 12.6 prop perm 1e-12 poro 0.55 range z -0.9 9.1 prop perm 1e-12 poro 0.5 range z -7.4 -0.9 prop perm 1e-12 poro 0.5 range z -11.4 -7.4 prop perm 1e-12 poro 0.45 range z -18.4 -11.4 prop perm 1e-12 poro 0.45 range z -30.4 -18.4 ini fdens 1e3 ini fmod 8.5e7 ini sat 0 range z 12.6 15.1 ini sat 1 range z -30.4 12.6 ini pp 0 grad 0 0 -1e4 range z -30.4 12.6 fix pp range x -.1 .1 fix pp range x 39.9 40.1 fix pp range y -.1 .1 fix pp range y 119.9 120.1 fix pp range z 12.5 12.7 fix pp range z -30.5 -30.3 ;material mechanic parameters model mohr def derive s_mod1=E_mod1/(2.0*(1.0+p_ratio1)) b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1)) s_mod2=E_mod2/(2.0*(1.0+p_ratio2))FLAC原理实例与应用指南FLAC3D实例分析教程
FLAC3D模拟实例循环开挖与支护
FLac3D接触面简单实例
FLAC3D条带开采案例(代码)
FLAC3D 实例命令流5
FLAC3D隧道耦合实例