1.1 FLAC常识.................................................................................................................. 2 1.2常用命令流 ............................................................................................................... 13 1.3建模过程 .................................................................................................................. 21 2-1定义一个FISH函数.................................................................................................. 24 2-2使用一个变量 ........................................................................................................... 24 2-3对变量和函数的理解................................................................................................. 24 2-4获取变量的历史记录................................................................................................. 24 2-5用FISH函数计算体积模量和剪砌模量...................................................................... 25 2-6 在FLAC输入中使用符号变量 ................................................................................. 25 2-7 控制循环 ................................................................................................................. 26 2-8 拆分命令行.............................................................................................................. 26 2-9 变量类型 ................................................................................................................. 27 2-10 IF条件语句 ............................................................................................................ 27 2-11 索单元自动生成 ..................................................................................................... 27 2-12圆形隧道开挖模拟计算 ........................................................................................... 28 4-1数组 ......................................................................................................................... 30 4-2函数操作 .................................................................................................................. 31 4-3函数删除与重定义 .................................................................................................... 32 4-4字符串...................................................................................................................... 32 4-5马蹄形隧道网格........................................................................................................ 33 4-6复杂形状网格生成 .................................................................................................... 33 4-7网格连接 .................................................................................................................. 34 4-8立方体洞穴网格生成................................................................................................. 35 4-9球体洞穴网格生成 .................................................................................................... 36 4-10应力边界 ................................................................................................................ 37 4-11改变应力边界条件 .................................................................................................. 37 4-12位移边界 ................................................................................................................ 39 4-13不考虑重力影响的均匀应力 .................................................................................... 41 4-14考虑应力梯度的均匀材料 ........................................................................................ 41 4-15考虑应力梯度的非均匀材料 .................................................................................... 42 4-16非均匀网格应力初始化 ........................................................................................... 42 4-17不规则自由面应力初始化 ........................................................................................ 43 4-18非均网格内部压实 .................................................................................................. 43 4-19模型改变后初始应力变化 ........................................................................................ 44 4-20应力与孔隙压力的初始化 ........................................................................................ 44 4-21加载顺序 ................................................................................................................ 45 6-1 V级围岩施工过程模拟 ............................................................................................. 47 6.2 IV级围岩施工过程模拟 ............................................................................................ 50 6.3 III级围岩施工过程模拟 ............................................................................................ 52 第七章 命令流按照顺序进行 ..........................................................................................
;右隧道中地层网格......................................................................................................... ;右隧道上侧土体网格 ..................................................................................................... ;右隧道下侧土体网格 ..................................................................................................... 第七章 命令流按照顺序进行 .......................................................................................... 57
1.1 FLAC常识
1. FLAC3D是有限元程序吗？ 答：不是！是有限差分法。
2. 最先需要掌握的命令有哪些？
答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。
3. 怎样看模型的样子？
答：plo blo gro可以看到不同的group的颜色分布
4. 怎样看模型的边界情况？ 答：plo gpfix red
5. 怎样看模型的体力分布？ 答：plo fap red
6. 怎样看模型的云图？
答：位移：plo con dis (xdis, ydis, zdis) 应力：plo con sz (sy, sx, sxy, syz, sxz)
7. 怎样看模型的矢量图？ 答：plo dis (xdis, ydis, zdis)
8. 怎样看模型有多少单元、节点？
答：pri info
9. 怎样输出模型的后处理图？
答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件
10. 怎样调用一个文件？ 答：File/call或者call命令
10. 如何施加面力？ 答：app nstress
11. 如何调整视图的大小、角度？
答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键
12. 如何进行边界约束？
答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束）
13. 如何知道每个单元的ID？
答：用鼠标双击单元的表面，可以知道单元的ID和坐标
14. 如何进行切片？
答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片)
15. 如何保存计算结果？ 答：save +文件名.
16. 如何调用已保存的结果？ 答：rest +文件名；或者File / Restore
17. 如何暂停计算？ 答：Esc
18. 如何在程序中进行暂停，并可恢复计算？
答：在命令中加入pause命令，用continue进行继续
19. 如何跳过某个计算步？
答：在计算中按空格键跳过本次计算，自动进入下一步
20. Fish是什么东西？
答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能
21. Fish是否一定要学？
答：可以不用，需要的时候查Mannual获得需要的变量就可以了
22. FLAC3D允许的命令文件格式有哪些？ 答：无所谓，只要是文本文件，什么后缀都可以
23. 如何调用一些可选模块？
答：config dyn (fluid, creep, cppudm)
24 .如何在圆柱体四周如何施加约束条件？
答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其
中r2 答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘 制就行了。 26.关于dd、dip的理解： 答： 首先你要确定你定义面的法线方向，dip 是倾角的意思（xoy面和定义面的夹角），法线向上dip是锐角，法线向下dip 时钝角。dd是y轴顺时针转到法线在xoy面上的投影线时，所转过的角度。 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.绘制相关曲线命令的一些解释 在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负); \"-\"表示对其值作\"mirror\如果大家使用AUTOCAD的话就明白\"mirror\"的含义了 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. a: fix x y z range x -0.1 0.1 b: fix x range x -0.1 0.1 请问a式与b式有何区别？ 答：a 表示 在x=0 这个平面 是固定的；b 表示在x=0 这个平面 只x方向不能发生运动 32.initial 与 apply 有何区别？ 答：初始化命令，如初始化计算体的应力状态等； apply边界条件命令，如施加边界的力、位移等约束等。 initial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件 不会发生变化。 33.如何查看各个时段不平衡力的具体数值？ 答：采用his来记录计算，包括位移应力等 命令his unbal his gp(zone) zdis range (0 0 0) 或者id=？ 导出数据命令his write n vs m begin 时步end 时步 +file filename.his 要导出不平衡力的具体数值 his unbal step 100000 or solve his write 1 vs step begin 1 end 1000 file 123.his 使用上述命令就可以查看各个时步下的不平衡力的具体数值 34.如何显示变形轮廓线的命令？ plo ske magf 10 其中10为放大系数 35、怎样开始入门？ 答：初步学习首先看说明书 User's gurid的Section2 Getting started和Section 3 Problem sloving with FLAC 36、哪些命令需要掌握？ 答：grid gen prop interface ini fix set hist step solve pause cont sav 等，在学习的过程 中遇到这些命令和不会的命令应该查看说明书的Command Reference，里面讲的很清楚，特别需要注 意的是一定要仔细看一下Command Reference开头的Optional Range Phrases，这是讲范围的，因为 很多命令的施加都需要加到模型的一部分，这里就告诉你怎样来表示这一部分。这一原则同样适合与 其他5个ITASCA软件。还有一点就是六个软件（FLAC FALC3D UDEC 3DEC PFC PFC3D）的range大不相 37、FLAC和FLAC3D建模有什么区别？ 答：FLAC和FLAC3D基本的建模思想大致是一样的，不同的是FLAC是首先用GRID命 令把整体的网格画出 来，然后用MODELL NULL命令把网格分成许多部分，再用GEN命令把这些部分一块一块搭接起来。然后 再用GEN LINE，GEN ARC等命令画出要开挖的边界，在FLAC中产生网格的命令只有一个GRID，而且只 能用一次。而FLAC3D是直接用GEN 命令一块一块搭接，GEN可以用很多次。 38、怎样查看模型？ 答：plot grid 可以查看网格，plot grid num 可以查看节点号。 39、怎样查看边界条件？ 答：plot fix 40、怎样查看预先留设的开挖边界？ 答：plot grid mark。 41、怎么查看模型的云图？ 答：pl syy (sxx xdis ydis sig1..) fill ，注意，带fill关键字就是查看云图，不带fill 关键字 就是查看等直线图。 42、怎样查看模型的矢量图？ 答：pl grid vel (disp) red 43、怎样查看INTERFACE？ 答：pl grid iface id red print iface id ; id 是interface的编号 44、FLAC怎么调整视图？ 答：ctrl+z 或shift+z ,然后点击要放大部位的左上角和右下角就可以放大某个部位，如果要缩小， 双击屏幕。 45、判断模型中网格点是一个还是两个的方法 显示模型网格点的ID，如果有两个ID号，说明是两个。节理面上的网格点应该有两个ID号。 46、改变模型应力边界时FLAC与3DEC的不同之处 在3DEC中，当改变应力边界时，模型表面的应力直接被更换，FLAC中改变应力边界时，模型 表面的应力为两者应力相加。 47、flac中影响运行时间的因素 （1）单元材料，结构面之间的刚度对比很大； （2）划分单元尺寸对比很大，尽量避免长细比大于5：1的细长单元。 48 FLAC中有多少中材料模型 内置11中本构模型%_4[ W/N,n R5q （1）零模型 用于开挖回填的区域 （2）各向同性弹性模型 （3）正交各向异性弹性模型 （4）横观各向同性弹性模型 （5）德鲁克－普拉格塑性模型 （6）莫尔－库仑塑性模型 一般岩土体的力学行为 （7）节理化塑性模型 用于研究材料破坏后行为（如边破稳定，屈服煤矿开采） （8）应变硬化/软化模型－库仑塑性模型 薄板层状材料破坏后力学行为 （9）双线性应变硬化/软化节理化塑性模型 临近层状土层开挖 （10）双屈服塑性模型 （11）修正的剑桥粘土模型 49、对不平衡力的理解 flac中用最大不平衡力和典型内力大的比值R表示模型的不平衡力的相对大小，该值为百分数，且从 不会减小到零，因此，模型不会达到绝对稳定的平横状态，可以根据要求的精度设定R值，模型默认 的R值是1e-5。 注意R值小仅表示所有网格点的力都平衡，但模型还可能发生塑性流动，通过 plot vel 命令查看 网格点速度，一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。 50、简单命令解释： plot add cont disp out on 时时显示云图 plot block group gpnum on 显示节点编号 plot inter red ske 看接触面的效果. del range group 1 not 不删除group信息 del range group 2 删除group信息 ini z add 数值大小 range group name 顺z方向移动 建立接触面 interface 1 face range x(范围）y(范围） z(范围） 或者 interface 1 face range plane orig(x y z 坐标） norm（法线方向） water table clear 清除水位线 pl set back white ;设置背景 set plot bitmap size (800,510) ;设置图片大小 plot set caption size 37 ; 设置图例数字大小 (default=35 范围10-50) plot set center 100 0 40 fix x y z range x -0.1 0.1 表示在x=0 这个平面是固定的 fixed boundary fix x range x -0.1 0.1表示在x=0 这个平面只x方向不能发生运动roller boundary 51、fish命令的解释： def aux1 定义函数 loop n(1,nptab1) 开始一个循环，从1循环到变量nptab1 rr=xtable(1,n) 从table 1 的x列中读第n个数付给变量rr zz=ytable(1,n) 从table 1 的y列中读第n个数付给变量zz loop k(1,nprof) 开始第二个循环，从1循环到变量nprof case_of k 下面是一个选择结构，利用变量k判断 i F alfa=0.5*(pi+.1)*float(k-2)/float(nprof-2) xtable(n+10,k)=rr*cos(alfa) ytable(n+10,k)=rr*sin(alfa) case 1 当k＝1的时候执行下面语句 xtable(n+10,k)=rr table n+10的x列第k个值等于变量rr ytable(n+10,k)=-40.0 table n+10的y列第k个值等于变量-40 case 2 当k＝2的时候执行下面语句 xtable(n+10,k)=rr table n+10的x列第k个值等于变量rr ytable(n+10,k)=0.0 table n+10的y列第k个值等于变量0 end_case 结束选择结构 end_loop 结束第二个循环 end_loop 结束第一个循环 end 函数定义结束 aux1 执行什么定义的函数 52. 在固定面上施加位移时，fix与apply不能同时使用，应该用ini和fix,才能施加上。 53.菜单 plot items/add/structure elements/geometry/结构单元的形式，把对话框里的mag fac的值改成不同的数，就可以看到位移放大后的结构单元。 .你可以用键盘上的上、下、左、右肩在窗口中移动模型。 55.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？ 答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。 56 在进行三轴压缩时候，如何施加围压？ 答：使用apply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius 1，其中X和Z轴是 在我们眼睛正对的面。 57 在进行三轴压缩时候，如何对轴向力进行位移控制？ 答：使用 ini yvel 1e-7 range y -.1 .1 ini yvel 1e-7 range y 1.9 2.1 即可，数值可以变化 58 在进行边坡计算的时候，为什么用sol, fos 命令运行不行？ 答：应该输入sol fos slope3dfos.sav associated(不要也行） 再用restore 命令调用即可。 59 在FLAC3D计算中，如何记录程序中的值，如果不想时时记录该怎么弄？ 答：用 hist wtite n1 vs n2 file.txt， 想间隔记录可用set hist_rep n。 60 想记录接触面上的相关力学参数该怎么办？ 答：hist inter i nstress /ndis/sstress/sdisp 61 接触面设置参数包括哪几个？ 答：coh dial fric kn ks ten. 62 如何输出单元格的数据？ 答；expgrid .flac3d 63 如何在剖面的沉降云图上增加网格、剖面上的矢量及坐标值？ D Q Z9s1} e4Y U%q 答：plo add ske plo add dis plane plo add axe 如何做切片？ 答：plot set plane ori x y z norm x y z 65 举一个隧道支护开挖的例子？ gen zone radbrick p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) & size 3 5 5 7 rat 1 1 1 1.5 dim 1 4 2 fill model mohr prop bulk 1e8 shear .3e8 fric 35 prop coh 1e3 tens 1e3 fix x range x -0.1 0.1 fix z range z -0.1 0.1 fix y range y 9.9 10.1 set grav 0 10 0 ini dens 1000 ini syy=0.0 grad 0,-10000,0 ini sxx=0.0 grad 0,-5000,0 ini szz=0.0 grad 0,-5000,0 apply sxx=0.0 grad 0,-5000,0 range x 9.9,10.1 apply szz=0.0 grad 0,-5000,0 range z 9.9,10.1 set large hist unbal solve save tr_eq.sav ; monitor trench excavation ini xdis=0 ydis=0 zdis=0 hist gp xdisp 1,0,0 hist gp zdisp 0,0,2 ; ; excavation steps range name exc1 x 0 1 y 0.0 0.8 z 0 2 range name exc2 x 0 1 y 0.8 1.6 z 0 2 range name exc3 x 0 1 y 1.6 2.4 z 0 2 range name exc4 x 0 1 y 2.4 3.2 z 0 2 range name exc5 x 0 1 y 3.2 4.0 z 0 2 ; excavation step 1 model null range exc1 step 10000 save exc1.sav ; excavation step 2 model null range exc2 step 1000 save exc2.sav ; excavation step 3 model null range exc3 step 1000 save exc3.sav ; excavation step 4 model null range exc4 step 1000 save exc4.sav ; excavation step 5 model null range exc5 step 1000 save exc5.sav ; install cable support rest exc2.sav sel cable beg 1.0,0.4,1.5 end 5.0,0.4,1.5 nseg 4 sel cable beg 1.0,0.4,0.5 end 5.0,0.4,0.5 nseg 4 sel cable beg 1.0,1.2,1.5 end 5.0,1.2,1.5 nseg 4 sel cable beg 1.0,1.2,0.5 end 5.0,1.2,0.5 nseg 4 sel cable prop emod 2e9 ytension 1e8 xcarea 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0 model null range exc3 step 1000 save cab3.sav sel cable beg 1.0,2.0,1.5 end 5.0,2.0,1.5 nseg 4 sel cable beg 1.0,2.0,0.5 end 5.0,2.0,0.5 nseg 4 sel cable prop emod 2e9 ytension 1e8 xcarea 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0 model null range exc4; step 1000 save cab4.sav sel cable beg 1.0,2.8,1.5 end 5.0,2.8,1.5 nseg 4 sel cable beg 1.0,2.8,0.5 end 5.0,2.8,0.5 nseg 4 sel cable prop emod 2e9 ytension 1e8 xcarea 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0 model null range exc5 step 1000 save cab5.sav plot create TrenchSupport plot set rot 63., 334.35, 3.4 plot set cent 2.3, 3.446, 3.9 plot set size 4.136 plot add cont disp out on plot add sel cable force plot add axes black 66 举一个三轴压缩的例子？ gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4 gen zone reflect norm 1,0,0 gen zone reflect norm 0,0,19 model mohr prop bulk 1.19e10 shear 1.1e10 prop coh 2.72e5 fric 44 ten 2e5 fix x y z range y -.1 .1 fix x y z range y 1.9 2.1 apply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius 1 ini yvel 0 grad 0 -1e-4 0 range y 0 2 sol 1.FLAC程式使用前准备步骤 步骤1：依比例画出所欲分析的资料 在纸上画出地点的位置、地层资料、并简单标示距离及深度资料。 步骤2：换算输入资料成同一单位 将现有地层资料，如 Density, Bulk modulus, Young`s modulus, tension, cohesion, friction Angle 等资料，换算成同一单 位。 附注 ：需谨慎检查输入资料的单位，如因单位不同而造成过大或过小的值，将会造成 FLAC无法计算，而产生ERROR信息。 步骤3：应用公式简略计算 应用公式或依据经验，简略算出FLAC输出资料的范围，以作为Debug及输出资料分析时验证。 步骤4：建立x, y坐标与node i, j之间的关系 在图上距离及深度的关系，建立x, y坐标系统，再由x, y坐标系统，转换与网格间系，为了便于以后输出资料的分析，故应确实 掌握网格的位置及其相对应的x, y 坐标。 建议在敏感区域使用较密的网格，其它地方则使用较疏的网格，刚开始执行程序时，不宜使用网格太大的网格数目，因尽量使网 格总数少于1000，以节省时间。 2.FLAC 3D 程序的编写步骤： 1 Config ________ 2 Grid ________ 3 Model ________ 4 求起始的应力平衡________ (1)建立x, y坐标与网格的关系，建议使用Gen指示： Gen x1,y1 x2,y2 ,x3,y3 x4,y4 i=i0,i1 j=j0,j1 详细指令参见使用手册，FLAC程式可自动产生x, y 坐标与网格的关系，但由 于产生的网格坐标不易控制，将对其它的工作产生负面影响，故依使用前步骤4所建立的关系，将网格依其疏密程度需要的不同，实 际控制网格的坐标。 (2)设定材料性质：prop (3)设定外力：Set Grav, Apply Pressure, ini sxx, Syy (4)设定边界条件：fix, free (5)求起始的应力平衡：solve (6)储存：Save 5 求工程的影响________ 求出区域内的应力分布情况后，再依工程的流程及步骤阶段执行各工程进行过程的影响，建议使用以下的步骤： (1)调出起初的应力平衡：re_____ .sav (2)设定新的材料性质：model,prop (3)设定新的支撑性质：struct (4)设定新的外力 (5)设定边界条件 (6)求工程时的应力平衡 (7)储存 1.2常用命令流 ;模型镜像 gen zone radcylinder size 25 1 25 25 gen zone reflect normal -1 0 0 origin x y z(面上一点）;沿X轴镜像,通过对称平面法线向量确定对称面 gen zone reflect normal 0 0 -1 ;沿z轴镜像 ;绘图控制 pl ;在模型中显示位移-应变曲线 hist gp ydisp 0,0,0 hist zone syy 0,1,0 hist zone syy 1,1,0 pl his -2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负); \"-\"表示对其值作\"mirror\" ;对模型进行压缩实验的方法 ;即在模型两侧施加相反方向的速度 ini yvel 1e-7 range y -.1 .1 ini yvel -1e-7 range y 1.9 2.1 ;修改模型的坐标值 ini x add -100 y add -100 z add -100 ;显示云图的同时也显示模型网格轮廓 plot add cont disp outline on ;输入角度、弧度方法 pi=π，90°为90.0*degrad def set_vals ptA = 25.0 * sin(pi/2);ptA=25.0 ptB = 25.0 * cos( 60.0*degrad );ptB=12.5 ptC=pi;ptC=3.1415926 end set_vals print ptA ptB ptC ;施加结构单元方法 sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) & end2=(0.0,25.0,0.0) radius=24.5 not plot add sel geom black black cid on scale=0.03 sel node init zpos add -25.0 ;如何显示某一平面 plot create name_plane plot set plane origin 3 4 0 normal 1 0 0 plot add cont disp plane behind shade on plot add sel geom black plot add axes red plot current name_plane plot show ;如何施加法向应力并显示所施加的力 apply nstress -3.0E7 range cylinder end1 0 -1 0 end2010 radius 39.9 not;柱体外表面施加法向力 apply nstress -5.0E6 range cylinder end1 0 -1 0 end2010 radius 2.1;柱体内表面施加法向力 plot grid green fap black;显示施加的力 ;获取某网格点的ID，坐标值 def aa get_id=gp_id(gp_near(3,0,0)) end aa print get_id print gp position range id 3 ;保存命令流文件 set log on set logfile hua.log;文件存放的路径与调用的文本路径相同 ...... set log off ;fish函数的说明与解释 def find_add ;定义fish函数find_add head = null ;给head赋值 p_gp = gp_head ;第一个网格结点的指针赋给p_gp loop_while p_gp # null ;当p_gp值不为null时作循环 x_pos = gp_xpos(p_gp) ;将指针为p_gp的结点的x坐标值赋给x_pos if x_pos = 10.0 then ;如果x_pos = 10.0 则(执行) new = get_mem(2) ;从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址 mem(new) = head ;将head类型和数值置于地址为new的fish变量 mem(new+1)=p_gp ;第二个变量存储符合条件的zone地址 head = new ;将new值赋给head endif p_gp = gp_next(p_gp) ;将结点指针为p_gp的下一个结点的指针赋给p_gp endloop;结束循环 end ;结束fish函数 ;实际上这个fish函数为满足条件（x坐标为10的）的zone的地址开辟一定的地址空间，各地址之间存在一定的联系 ;找到第一个符合条件的zone地址后，用new = get_mem(2) 从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址，并用;下面的mem(new) = head ，将第一次开辟的两个变量的第一个变量存储地址head（注意第一次head＝null），第二个变量存储第一;个符合条件的zone地址，并将第一个变量的地址赋予head（head＝new），第一次循环结束；下次循环，同样开辟两个变量对象，第;一个变量对象记录上次循环开辟的第一个变量的地址，第二个记录第二个符合条件的zone地址，其余循环依次类推，这样子就建立了;一个符合条件的zone地址链条，方便以后使用。 ；FLAC后处理 ；输出JPG图片 set plot jpg set plot quality 100 plot hard file XX.jpg ;输出TXT文件 hist write X file xx.txt ；(x是历史记录ID） print con zd ;输出单元应力和节点位移 set log on set logfile xx.log print zo stress print gp dis set log off ；生成影片文件 plot set rot 20 0 30 plot con szz ou on magf 10 plot add hist 1 set movie avi step 1 file xx.avi movie start solve ;计算间隔命令 在solve命令之前 set hist_rep 1 ；(步数） ；检查模型赋值情况 plot bcon prop bu ;变形云图加网格 放大20 倍 plot con zd ou on magf 20 ;变形网格 plot sk magf 20 ;变形矢量图（加网格） plot sk dis scale 0.07 ；(scale是箭头长度，默认0.05） ；如何观察塑性区 plot block state plot block state sh_n plot block state sh_p ;切片操作 plot set plane ori 1 2 3 norm 0 0 -1 plot con zd plane plot add dis plane plot add sk plot add ax ;切片后方部分显示 plot con zd plane behind shade on !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 显示某一个group plot block group range group 1 显示监测的数据 plot hist 112 ;112为hist编号 FLAC有用的命令流 显示命令流 位移：plo con dis (xdis, ydis, zdis) plot contour smin range group \"地层\" plot add contour smin range group \"初衬上\" \"上台阶\" \"核心土\" \"中台阶\" \"下台阶\" \"初衬仰拱\" \"仰拱\" any 同时选择多个group plot contour smax range group \"初衬上\" any group \"初衬中\" any group \"初衬下\" any group \"初衬仰拱\" any plot contour smin range group \"初衬上\" any group \"初衬中\" any group \"初衬下\" any group \"初衬仰拱\" any plot contour smin range group \"初衬上\" not group \"初衬中\" not group \"初衬下\" not group \"初衬仰拱\" not 排除某个group plot block state range group \"地层\" not 开挖某个group中的一部分 model null range group 1 y yy_1 yy_2 ************************************************************* 拱顶沉降 set logfile 拱顶沉降节点位置.log set log on print gp position range group \"地层\" x -0.1 0.1 y 0 50 z 12 13 set log off set logfile 拱顶沉降.log set log on print gp displacement range group \"地层\" x -0.1 0.1 y 0 50 z 12 13 set log off ************************************************************* ************************************************************* ;上台阶水平收敛 set logfile 上台阶水平收敛节点位置.log set log on print gp position range group \"地层\" x -7 7 y 0 50 z 7.3 7.8 set log off set logfile 上台阶水平收敛.log set log on print gp displacement range group \"地层\" x -7 7 y 0 50 z 7.3 7.8 set log off ;下台阶水平收敛 set logfile 下台阶水平收敛节点位置.log set log on print gp position range group \"地层\" x -7.2 7.2 y 0 50 z 4.4 4.8 set log off set logfile 下台阶水平收敛.log set log on print gp displacement range group \"地层\" x -7.2 7.2 y 0 50 z 4.4 4.8 set log off ************************************************************* ************************************************************* ;下台阶水平收敛 set logfile 上台阶水平收敛节点位置.log set log on print gp position range group \"地层\" x -7 7 y 0 50 z 7.3 7.8 set log off set logfile 上台阶水平收敛.log set log on print gp displacement range group \"地层\" x -7 7 y 0 50 z 7.3 7.8 set log off ;下台阶水平收敛 set logfile 下台阶水平收敛节点位置.log set log on print gp position range group \"地层\" x -7.2 7.2 y 0 50 z 4.4 4.8 set log off set logfile 下台阶水平收敛.log set log on print gp displacement range group \"地层\" x -7.2 7.2 y 0 50 z 4.4 4.8 set log off ************************************************************* 应力：plo con sz (sy, sx, sxy, syz, sxz) 怎样看模型的矢量图？ 答：plo dis (xdis, ydis, zdis)5 14. 如何进行切片？ 答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 监测位移应力 history id 1 gp displacement x y z hist gp zdis 0.034 0 60.044 ;地表 history id 2 gp displacement id=4442 hist id 31 gp zdisplacement id 3161;桩底中点位移 hist id 32 zone szz id 6818;桩底左侧，桩底共8个单元 导出监测点数据的命令： his write 11 vs step begin 1 end 4000 file 123.his;导出11号监测点从时步1至时步4000的数据， 加边界条件 边界力，压应力为负，拉应力为正 设置了固定边界条件时，就不能再在边界上施加应力。 apply szz 8e6 range z 0.9 1.1 初始化地应力,压应力为负 ini sxx -3.5e6 syy -3.5e6 szz -3.5e6 画等值线 1、导出数据 +++++++++++++++++++++++++++++++++++++++ res zizhong.sav ;data_process set logfile 区域信息.log set log on print zone information range group \"初衬上\" any group \"初衬中\" any group \"初衬下\" any group \"初衬仰拱\" any set log off set logfile 挖前主应力.log set log on print zone principal range set log off set logfile 挖前应力.log set log on print zone stress set log off set logfile 挖前节点位置.log set log on print gp position range x -0.1 0.1 z 7.5 7.7 set log off ;set logfile 挖前位移.log ;set log on ;print gp displacement range x -0.1 0.1 z 7.5 7.7 ;set log off ***************************************************** 定义interface interface 1 face range plane norm (-1,0,1) origin (1.5,1.5,3) dist 0.1 inter 2 face range z 2.9,3.1 inter 1 prop kn 300e9 ks 300e9 tens 1e10 SBRATIO=1 ;移动group ini z add -1.0 range group Top ;Top是上文定义的一个group ;定义group group Material group Bin range group Material not range name=Bin group Bin ;等号可不要 range name=Material group Material model mohr range Material 定义一个新的group 支护结构，这个group 支护结构是group 隧道的一部分，并且位于两个圆柱体之间的区域 group 支护结构 range group 隧道 cyl end1 0 0 5.22 end2 0 50 5.22 rad 7.16 cyl end1 0 0 5.22 end2 0 50 5.22 rad 7.00 not group 支护结构 range group 隧道 cyl end1 0 0 13.516 end2 0 50 13.516 rad 13.52 cyl end1 0 0 13.516 end2 0 50 13.516 rad 13.37 not group 支护结构 range group 隧道 cyl end1 3.3684 0 4.13 end2 3.3684 50 4.13 rad 3.55 cyl end1 3.3684 0 4.13 end2 3.3684 50 4.13 rad 3.3900 not x 5.49 6.53 z 1.285 2.495 group 支护结构 range group 隧道 cyl end1 -3.3684 0 4.13 end2 -3.3684 50 4.13 rad 3.55 cyl end1 -3.3684 0 4.13 end2 -3.3684 50 4.13 rad 3.3900 not x -6.53 -5.49 z 1.285 2.495 ;创建显示图 plot create view_inter plot add surface plot add interface red plot show ;清零位移 ini xdis 0 ini ydis 0 ini zdis 0 隧道埋深的影响 apply szz -5.9e6 range z 52 52.2 开挖一个group中的一段 model null range group 1 y 0 90 ; 导出指定位置的数据 set logfile 挖前节点位置.log set log on print gp position range x -0.1 0.1 z 7.5 7.7 set log off ;set logfile 挖前位移.log ;set log on ;print gp displacement range x -0.1 0.1 z 7.5 7.7 ;set log off 1.3建模过程 1、在CAD中建模 a，模型中的所有线条都须改为直线段，不能用多段线，线段与线段之间不能重叠。线段数量不宜太多； b，将模型的高程与地形图高程对应。将模型端点坐标放置恰当位置（端点坐标最好为整数）。确保模型在第一象限。 c，将模型保存为“DXF”格式； 2、将DXF格式的模型导入algor中 a，打开algor，点击“SD2”－file－import－autocad－F10－enclose b，此时，algor中将出现模型的轮廓 c，将模型导出，file－export－iges d，此时得到igs格式的模型 新方法： 1、在CAD中将点数据导出：若地形线为多段线，则选中多段线，输入list命令，即可将多段线各点导出 2、用excel编辑点，用ansys创建点 3、用spline线创建线 3、在ansys中处理模型 a，打开ansys软件－file－import－iges；将模型导入ansys中 删除line only，创建spline线，创建大面，用divide将大面切割成小面。面被线divide时，线消失。 创建面出现问题时可能是点出现重复，merge一下即可。 b，设置单元类型：element type－add－solid－quad 4node 42－ c，设置材料属性：material props－material models;模型中有多少个group有建多少个材料属性，材料属性可选择structural-density,赋值应赋不同的值，以便merger。 d，创建面单元：create-areas-arbitrary-by lines e，给每个面赋材料属性：meshtool-给每个面赋材料属性－element attributes-areas-set；面单元编号与材料属性编号应对应 d，网分：smart size自动网分-Tri三角形－Free自由映射－mesh（大小一致）设置line的长度-size controls-lines-set 选择等边长网分时，注意用clear命令清楚等边长的间隔，如未clear，这线为虚线。 e，合并节点和单元，重新编号：numbering ctrls-merge items-notes/elements numbering ctrls-compress numbers-notes/elements 网分后不能合并所有的对象，只能notes/elements，如果all则所有材料合并为一个 当材料属性所附的参数不相等时，可以在网分后merge all。 f，导出节点和单元数据:create-nodes-write node file-node.dat-ok create-elements-write elemet file-elment.dat-ok 4、利用ansys与flac的转换程序，将模型转入flac 平面模型转换-solid 42 最新转换程序： 平面模型转换-Grid文件-flac模型文件，导出的是一个模型，而非命令流，但速度很快。 打开flac，file-imput grid-已转换好的grid模型文件。 group 冲积物 range group 2 group 支护 range group 隧道 group 开挖体 range group \"已开挖\" 5、在flac中计算： 参数：剪切模量G=E/2(1+u),体积模量K=E/3(1-2u) 各岩石的参数： model mohr prop bulk 28.33e7 shear 21.5e6 coh 5000 fric 30.4 dens 1950 range group 1 ;粉细砂，弹模34Mpa，泊松比0.3 prop bulk 27.778e7 shear 15.8736e7 coh 6000 fric 33.2 dens 2050 range group 2 ;粗砂，弹模40Mpa，泊松比0.26 prop bulk 2.57e9 shear 1.77e9 coh 1.2e6 fric 40 dens 2320 range group 3 ;泥岩，弹模4320Mpa，泊松比0.22 prop bulk 2.57e9 shear 1.77e9 coh 1.2e6 fric 40 dens 2320 range group 4 ;泥岩，弹模4320Mpa，泊松比0.22 prop bulk 16.67e9 shear 12.5e9 dens 2400 range group pile;C30号混凝土，弹模30Gpa，泊松比0.2 prop bulk 16.96e9 shear 11.934e9 coh 1.78e6 fric 47 dens 2530 range group 1 ;低品位石灰岩覆盖层泊松比0.215 prop bulk 0.67e9 shear 0.22e9 coh 0.03e6 fric 21 dens 2250 range group 2 ;粘土夹层泊松比0.35 prop bulk 17e9 shear 12e9 coh 1.8e6 fric 47 dens 2580 range group 3 ;石灰岩矿石层泊松比0.22 prop bulk 6.67e9 shear 2.22e9 coh 0.7e6 fric 39 dens 2450 range group 2 ;四类围岩，变模6000Mpa，泊松比0.35，内摩擦角39，粘聚力0.7Mpa，密度2450kg/m3 property she=0.09e9 bulk=0.2e9 coh=0.04e6 fric =27 ten=0.02e6 den=2020 range gro m-1;河谷堆积物 property she=0.5e9 bulk=1.4e9 coh=0.2e6 fric =27 ten=0.2e6 den=2250 range gro m-3;强风化带，5类围岩 property she=1.4e9 bulk=3.5e9 coh=0.7e6 fric =33 ten=0.45e6 den=2325 range gro m-5;弱风化带，4类围岩 property she=5.1e9 bulk=9.6e9 coh=1.1e6 fric =44.5 ten=7e6 den=2550 range gro m-7;基岩，3类围岩 property she=0.09e9 bulk=0.2e9 coh=0.04e6 fric =40 ten=0.02e6 den=2020 range gro m-1;河床冲积物 prop bulk 16.67e9 shear 12.5e9 dens 2400 range group m-1;C30号混凝土，弹模30Gpa，泊松比0.2 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) sum = sum + n prod = prod * n end_loop end xxx print sum, prod new gen zone brick p0 (0,0,0) p1 (-10,0,0) p2 (0,10,0) p3 (0,0,-10) model elas plot set rotation 0 0 45 plot block group def install pnt = zone_head loop while pnt #null z_depth = -z_zcen(pnt) y_mod = y_zero + cc * sqrt(z_depth) z_prop(pnt, ’shear’) = y_mod / (2.0*(1.0+p_ratio)) z_prop(pnt, ’bulk’) = y_mod / (3.0*(1.0-2.0*p_ratio)) pnt = zone_next(pnt) end_loop end set p_ratio=0.25 y_zero=1e7 cc=1e8 install 2-8 拆分命令行 new ;example of a sum of many things def long_sum temp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10 long_sum = temp + v11 + v12 + v13 + v14 + v15 end 2-9 变量类型 new def haveone aa = 2 bb = 3.4 cc = ’Have a nice day’ dd = aa * bb ee = cc + ’, old chap’ end haveone print fish 2-10 IF条件语句 new def abc if xx > 0 then abc = 1000 else abc = -1000 end_if end set xx = 10 print abc set xx = 0 print abc 2-11 索单元自动生成 new gen zone brick size 10 3 5 plot set rotation 0 0 45 plot block group def place_cables loop n (1,5) z_d = float(n) - 0.5 command sel cable beg 0.0,1.5,z_d end end_command end_loop end place_cables plot grid sel geom red 7.0,1.5,z_d nseg 7 new gen zone brick size 10 3 5 plot set rotation 15 0 60 plot block group mod mohr prop bulk 1e8 shear .3e8 fric 35 prop coh 1e3 tens 1e3 ini dens 1000 set grav 0,0,-10 fix x y z range z -.1 .1 fix y range y -.1 .1 fix y range y 2.9 3.1 fix x range x -.1 .1 fix x range x 9.9 10.1 set large hist unbal solve save cab_str.sav ini xdis 0 ydis 0 zdis 0 hist gp xdisp 0,1,5 def place_cables loop n (1,5) z_d = 5.5 - float(n) z_t = z_d + 0.5 z_b = z_d - 0.5 command free x range x -.1,.1 z z_b z_t solve sel cable beg 0.0,0.5,z_d end 7.0,0.5,z_d nseg 7 sel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7 sel cable beg 0.0,2.5,z_d end 7.0,2.5,z_d nseg 7 sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 & gr_k 2e10 gr_coh 1e10 gr_per 1.0 end_command end_loop end place_cables save cab_end.sav plot sketch sel cable force red 2-12圆形隧道开挖模拟计算 ;建立模型 gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 & size 4 2 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outsiderock gen zone cshell p0 0 0 0 p1 3 0 0 p2 0 1 0 p3 0 0 3 & size 1 2 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group concretliner fill group insiderock gen zon reflect dip 90 dd 90 orig 0 0 0 gen zon reflect dip 0 dd 0 ori 0 0 0 gen zon brick p0 0 0 6 p1 6 0 6 p2 0 1 6 p3 0 0 13 size 4 2 6 group outsiderock1 gen zon brick p0 0 0 -12 p1 6 0 -12 p2 0 1 -12 p3 0 0 -6 size 4 2 5 group outsiderock2 gen zon brick p0 6 0 0 p1 21 0 0 p2 6 1 0 p3 6 0 6 size 10 2 4 group outsiderock3 gen zon reflect dip 0 dd 0 orig 0 0 0 range group outsiderock3 gen zon brick p0 6 0 6 p1 21 0 6 p2 6 1 6 p3 6 0 13 size 10 2 6 group outsiderock4 gen zon brick p0 6 0 -12 p1 21 0 -12 p2 6 1 -12 p3 6 0 -6 size 10 2 5 group outsiderock5 gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z 6.1 13.1 gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z -6.1 -12.1 gen zon reflect dip 90 dd 90 orig 0 0 0 range x 6.1 21.1 z -12.1 13.1 ;绘制模型图 plot block group plot add axes red ;plot set rotation 0 0 45 用于显示三维模型 ;设置重力 set gravity 0 0 -10 ;给定边界条件 fix z range z -12.01,-11.99 fix x range x -21.01,-20.99 fix x range x 20.99,21.01 fix y range y -0.01 0.01 fix y range y 0.99,1.01 ;求解自重应力场 model mohr ini density 1800 ;围岩的密度 prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4 ;体积、剪切、摩擦角、凝聚力、抗拉强度 set mech ratio=1e-4 solve save Gravsol.sav plot cont zdisp outl on plot cont szz ;毛洞开挖计算 initial xdisp=0 ydisp=0 zdisp=0 model null range group insiderock any group concretliner any plot block group plot add axes red set mech ratio=5e-4 solve save Kaiwsol.sav plot cont zdisp plot cont sdisp plot cont szz plot cont xzz ;模筑衬砌计算 model elas range group concretliner any plot block group plot add axes red ini density 2500 range group concretliner any ;衬砌混凝土的密度 prop bulk=16.67e9,shear=12.5e9 range group concretliner any ;衬砌混凝土的体积弹模、剪切弹模 set mech ratio=1e-4 solve save zhihusol.sav plot cont zdisp plot cont sdisp plot cont szz plot cont xzz ;完成计算分析 第四章 命令流按照顺序进行 4-1数组 new def afill ;fill matrix with random numbers array var(4,3) loop m (1,4) loop n (1,3) var(m,n) = urand endloop endloop end def ashow ;display contents of matrix loop m (1,4) hed = ’ ’ msg = ’ ’+string(m) loop n (1,3) hed = hed + ’ ’+string(n) msg = msg + ’ ’+string(var(m,n)) endloop if m = 1 dum = out(hed) endif dum = out(msg) endloop end afill Ashow 4-2函数操作 new def xxx aa = 2 * 3 xxx = aa + bb end new def stress_sum stress_sum = 0.0 pnt = zone head loop while pnt # null Stress_sum = stress_sum + z_sxx(pnt) pnt=z_next(pnt) end_loop end new def stress_sum sum = 0.0 pnt = zone head loop while pnt # null sum = sum + z_sxx(pnt) pnt=z next(pnt) end loop stress_sum = sum end new def h_var_1 ipz = z_near(1,2,3) H_var_1 = sxx(ipz) + fac H_var_2 = syy(ipz) + fac H_var_3 = szz(ipz) + fac H_var_4 = sxy(ipz) + fac H_var_5 = syz(ipz) + fac H_var_6 = sxz(ipz) + fac end 4-3函数删除与重定义 new def joe ii=out(‘This is A Function’) end def fred joe end fred def joe ii=(‘This is A New Function’) end fred 4-4字符串 new def in_def xx = in(msg+’(’+’default:’+string(default)+’):’) if type(xx) = 3 in_def = default else in_def = xx endif end ; def moduli_data default = 1.0e9 msg=’Input Young‘s modulus ’ Y_mod = in_def default = 0.25 msg=’Input Poisson‘s ratio ’ p_ratio = in_def if p_ratio = 0.5 then ii = out(’ Bulk mod is undefined at Poisson‘s ratio = 0.5’) ii = out(’ Select a different value --’) p_ratio = in_def endif s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end ; moduli_data gen zone brick size 2,2,2 model elastic prop bulk = b_mod shear = s_mod print p_ratio y_mod b_mod s_mod pause print zone prop bulk pause print zone prop shear 4-5马蹄形隧道网格 new gen zone radcyl size 5 10 6 12 rat 1 1 1 1.2 p0 0,0,0 p1 100,0,0 p2 0,200,0 p3 0,0,100 gen zone radtun size 5 10 5 12 rat 1 1 1 1.2 p0 0,0,0 p1 0,0,-100 p2 0,200,0 p3 100,0,0 gen zone reflect dip 90 dd 270 origin 0,0,0 plot set rotation 30 0 30 plot block group plo add axes red 4-6复杂形状网格生成 new ; 主隧道网格 gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 size 4 10 6 4 dim 4 4 4 4 rat 1 1 1 1 fill gen zon reflect dip 90 dd 90 orig 15 0 0 gen zon reflect dip 0 ori 0 0 0 ; 辅助隧道网格 gen point id 1 (2.969848,0.0,-0.575736) gen point id 2 (2.969848,50.0,-0.575736) gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 & p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 gen zon radcyl p0 0 0 -1 p1 0 0 -8 p2 0 50 -1 p3 7 0 0 p4 0 50 -8 & p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 ;衬砌网格 sel shell range cyl end1 0 0 -1 end2 0 50 -1 rad 3 ; 隧道外围边界网格 gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 & p5 15 50 50 p6 50 0 50 p7 50 50 50 & p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 & size 6 10 3 10 rat 1 1 1 1.1 gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 & p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 & size 3 10 10 rat 1 1 1.1 ;对称得出1/2模型 gen zon reflect dip 0 ori 0 0 0 range x 0 23 y 0 50 z 8 50 gen zon reflect dip 0 ori 0 0 0 range x 23 50 y 0 50 z 0 50 ;建立主隧道和辅助隧道块名称 group service range cyl end1 0 0 -1 end2 0 50 -1 rad 3 group main range cyl end1 15 0 0 end2 15 50 0 rad 4 ;对称得出完整模型 gen zon reflect dip 90 dd 270 ori 0 0 0 ;网格显示 plot set rotation 30 0 30 plot block group 4-7网格连接 new gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,2 gen zone brick size 8 8 4 p0 0,0,2 p1 4,0,2 p2 0,4,2 p3 0,0,4 attach face range z 1.9 2.1 model elas prop bulk 8e9 shear 5e9 fix z range z -.1 .1 fix x range x -.1 .1 fix x range x 3.9 4.1 fix y range y -.1 .1 fix y range y 3.9 4.1 apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2 hist unbal ;网格显示 plot set rotation 30 0 30 plot block group solve save att.sav plot cont zdisp outl on 单一网格（小）程序 new gen zone brick size 8 8 8 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elas prop bulk 8e9 shear 5e9 fix z range z -.1 .1 fix x range x -.1 .1 fix x range x 3.9 4.1 fix y range y -.1 .1 fix y range y 3.9 4.1 apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2 hist unbal ;网格显示 plot set rotation 30 0 30 plot block group solve save noattx.sav plot cont zdisp outl on 单一网格（大）程序 new gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elas prop bulk 8e9 shear 5e9 fix z range z -.1 .1 fix x range x -.1 .1 fix x range x 3.9 4.1 fix y range y -.1 .1 fix y range y 3.9 4.1 apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2 hist unbal ;网格显示 plot set rotation 30 0 30 plot block group solve save noattx.sav plot cont zdisp outl on 4-8立方体洞穴网格生成 new def parm rad=4.0 len=10.0 in_size=6 rad_size=10 end parm gen zone radbrick edge len size in_size in_size in_size rad_size & rat 1.0 1.0 1.0 1.2 dim rad rad rad ;对称得出完整模型 gen zon reflect dip 0 dd 0 ori 0 0 0 gen zon reflect dip 90 dd 270 ori 0 0 0 ;网格显示 plot set rotation 30 0 30 plot block group 4-9球体洞穴网格生成 new def parm rad=4.0 len=10.0 in_size=6 rad_size=10 end parm def make_sphere p_gp=gp_head loop while p_gp#null px=gp_xpos(p_gp) py=gp_ypos(p_gp) pz=gp_zpos(p_gp) dist=sqrt(px*px+py*py+pz*pz) if dist>0 then k=rad/dist ax=px*k ay=py*k az=pz*k maxp=max(px,max(py,pz)) k=len/maxp bx=px*k by=py*k bz=pz*k u=(maxp-rad)/(len-rad) gp_xpos(p_gp)=ax+u*(bx-ax) gp_ypos(p_gp)=ay+u*(by-ay) gp_zpos(p_gp)=az+u*(bz-az) end_if p_gp=gp_next(p_gp) end_loop end gen zone radbrick edge len size in_size in_size in_size rad_size & rat 1.0 1.0 1.0 1.2 dim rad rad rad make_sphere ;对称得出完整模型 gen zon reflect dip 0 dd 0 ori 0 0 0 gen zon reflect dip 90 dd 270 ori 0 0 0 ;网格显示 plot set rotation 15 0 30 plot block group 4-10应力边界 new gen zone brick size (4,4,4) p0 (0,0,0) p1 (4,0,0) & p2 (0,4,0) p3 (2,0,3.4) model elastic prop bulk 1e8 shear .3e8 apply nstress -1e6 range plane dip 60 dd 270 origin 0.1,0,0 above step 1 ;网格显示 plot set rotation 30 0 30 plot block group plot grid red fap green 4-11改变应力边界条件 new gen zone brick size 6 6 6 model elas prop bulk 1e8 shear 7e7 fix x range x -0.1 0.1 def superstep loop ns (1,n_steps) x_stress = stress_inc command apply sxx add x_stress range x 5.9,6.1 y 0,6 z 0,2 step 100 end_command end_loop end set n_steps=100 stress_inc=-1e3 hist zone sxx 6,0,0 plot create sxx_hist plot add hist 1 plot show Superstep ;网格显示 plot set rotation 30 0 30 plot block group plot grid red fap green new gen zone brick size 6 6 6 model elas prop bulk 1e8 shear 7e7 fix x range x -0.1 0.1 def x_stress x_stress = stress_inc * step end set stress_inc = -1e3 apply sxx 1.0 hist x_stress range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0 hist x_stress step 100 ;网格显示 plot set rotation 30 0 30 plot grid red fap green new gen zone brick size 6 6 6 model elas prop bulk 1e8 shear 7e7 fix x range x -0.1 0.1 table 1 0,0 100,-1e5 apply sxx 1.0 hist table 1 range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0 step 100;网格显示 plot set rotation 30 0 30 plot grid red fap green new gen zone brick size 6,6,6 p1 6,0,-1 model elas prop bulk 8e9 shear 5e9 apply sxx -2e6 range x -0.1 0.1 apply sxx -2e6 range x 5.9 6.1 step 500 plot bound fap plot grid disp new gen zone brick size 5,5,5 model elas prop bulk 8e9 shear 5e9 set grav 0 0 -10 ini dens 1000 fix x range x -.1 .1 fix x range x 4.9 5.1 fix y range y -.1 .1 fix y range y 4.9 5.1 ini szz -5e4 grad 0 0 -1e4 app szz -5e4 range z -.1 .1 solve model null range x 1,4 y 1 4 z 3 5 step 100 plot set plane dip 90 dd 180 origin 0,2.5,0 plot add surf plane behind yell plot add vel plane behind black plot set rotation 30 0 30 plot block group plot show 4-12位移边界 new gen zone brick size 4,4,4 p3 2,0,3.4 model elas prop bulk 1e8 shear .3e8 macro left_boun ’plane dip 60 dd 270 origin 0.1,0,0 above’ macro right_boun ’plane dip 60 dd 270 origin 3.9,0,0 below’ apply nvel 0.1 plane dip 60 dd 270 range left_boun apply nvel 0.1 plane dip 120 dd 90 range right_boun step 1 ;网格显示 plot set rotation 30 0 30 plot block group plot grid vel new new gen zone brick size 4 4 4 model elastic prop bulk 1e8 shear .3e8 apply nstress -1e6 plane dip 0 dd 0 range y 3.9 4.1 apply nvel 0.0 plane dip 0 dd 0 range y -.1 .1 apply nvel 0.0 plane dip 0 dd 0 range x -.1 .1 step 10 ;网格显示 plot set rotation 30 0 30 plot block group plot grid vel new gen zone brick size 10 5 5 mod el prop shear 1e8 bulk 2e8 fix x y z range x -.1 .1 y 0 5 z 0 5 fix x y z range x 0 10 y 0 5 z -.1 .1 fix x y z range x 9.9 10.1 y 0 5 z 0 5 table 1 def find_add head = null p_gp = gp_head loop while p_gp # null x_pos = gp_xpos(p_gp) if x_pos = width then new = get_mem(2) mem(new) = head mem(new+1) = p_gp head = new endif p_gp = gp_next(p_gp) endloop end set width=10.0 find_add def apply_vel while_stepping ad = head loop while ad # null p_gp = mem(ad+1) gp_xvel(p_gp) = vel_max * gp_zpos(p_gp) / height gp_zvel(p_gp) = -vel_max * (gp_xpos(p_gp) - width) / height ad = mem(ad) endloop end set large vel_max=1e-2 height=5.0 step 100 ;网格显示 plot set rotation 30 0 30 plot block group plot grid vel 4-13不考虑重力影响的均匀应力 new ;set grav 0 0 0 gen zone brick size 6 6 6 model elas ini sxx=-5e6 syy=-1e7 szz=-2e7 apply sxx=-5e6 range x -0.1 0.1 apply sxx=-5e6 range x 5.9 6.1 apply syy=-1e7 range y -0.1 0.1 apply syy=-1e7 range y 5.9 6.1 apply szz=-2e7 range z -0.1 0.1 apply szz=-2e7 range z 5.9 6.1 prop bulk 8e9 shear 5e9 ini dens 2000 fix x range z -.1 .1 step 10 plot set rotation 30 0 30 ;网格三维显示 plot cont zdisp outl on 4-14考虑应力梯度的均匀材料 new gen zone brick size 10 10 10 p1 20,0,0 p2 0,20,0 p3 0,0,20 model mohr prop bulk 5e9 shear 3e9 fric 35 ini density 2500 set gravity 0,0,-10 fix x y z range z -0.1 0.1 ini szz -5.0e6 grad 0,0,2.5e4 ini syy -2.5e6 grad 0,0,1.25e4 ini sxx -2.5e6 grad 0,0,1.25e4 apply szz -4.5e6 range z 19.9 20.1 apply szz -5.0e6 range z -0.1 0.1 apply sxx -2.5e6 grad 0,0,1.25e4 range x -0.1 0.1 apply sxx -2.5e6 grad 0,0,1.25e4 range x 19.9 20.1 apply syy -2.5e6 grad 0,0,1.25e4 range y -0.1 0.1 apply syy -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1 step 10 ;网格三维显示 plot set rotation 30 0 30 plot cont zdisp outl on 4-15考虑应力梯度的非均匀材料 new gen zone brick size 10 10 10 & p0 0,-25,0 p1 20,-25,0 p2 0,0,0 p3 0,-25,20 model elas prop bulk 5e9 shear 3e9 ini density 1600 range y -10,0 ini density 2000 range y -15,-10 ini density 2200 range y -25,-15 set gravity 0,-10,0 fix x range x -.1 .1 fix x range x 19.9 20.1 fix z range z -.1 .1 fix z range z 19.9 20.1 fix y range y -25.1 -24.9 ini syy 0.0 grad 0,1.6e4,0 range y -10,0 ini syy 4e4 grad 0,2.0e4,0 range y -15,-10 ini syy 7e4 grad 0,2.2e4,0 range y -15,-25 step 10 ;网格三维显示 plot set rotation 30 0 30 plot cont zdisp outl on 4-16非均匀网格应力初始化 new gen zone radcyl size 3 8 4 5 fill p1 10,0,0 p2 0,10,0 p3 0,0,10 mode elastic prop shear 3e8 bulk 5e8 fix x range x -.1 .1 fix x range x 9.9 10.1 fix y range y -.1 .1 fix y range y 9.9 10.1 fix z range z -.1 .1 ini szz = -2.5e5 grad 0,0,2.5e4 ini density 2500 set grav 0,0,-10 step 10 ;网格三维显示 plot set rotation 30 0 30 plot cont zdisp outl on 4-17不规则自由面应力初始化 gen zone brick size 15 15 10 p0 0,0,0 edge=100.0 model elastic prop shear 3e8 bulk 5e8 def mountain gp = gp_head loop while gp # null zz = sqrt(gp_xpos(gp)ˆ2 + gp_ypos(gp)ˆ2) dz = 0.06 * sin(0.2 * zz + 100.0) ; Sum Fourier terms for dz = dz + 0.06 * sin(0.22 * zz - 20.3) ; quasi-random surface dz = dz - 0.04 * sin(0.33 * zz + 33.3) ; topology. gp_zpos(gp) = 0.5 * gp_zpos(gp) * (1.0 + dz) gp = gp_next(gp) end_loop end mountain fix x range x -.1 .1 fix x range x 99.9 100.1 fix y range y -.1 .1 fix y range y 99.9 100.1 fix z range z -.1 .1 set grav 0,0,-10 ini density=2000 ini szz=-2.0e6 (grad 0,0,2.0e4) sxx=-4.0e6 (grad 0,0,4.0e4) syy=-4.0e6 (grad 0,0,4.0e4) step 100 plot set rotation 30 0 30 plot cont zdisp outl on;网格三维显示 4-18非均网格内部压实 new gen zone brick size 8 8 10 ratio 1.2 1 1 model mohr ini dens 2000 prop bulk 2e8 shear 1e8 prop fric 30 fix x range x -.1 .1 fix x range x 7.9 8.1 fix y range y -.1 .1 fix y range y 7.9 8.1 fix z range z -.1 .1 ;ini szz -2.0e5 grad 0,0,2e4 ;ini sxx -1.5e5 grad 0,0,1.5e4 ;ini syy -1.5e5 grad 0,0,1.5e4 set grav 10 step 1000 ;pause ;prop tens 1e10 coh 1e10 ;step 750 ;prop tens 0 coh 0 ;step 250 ;网格三维显示 plot set rotation 30 0 30 plot cont zdisp outl on 4-19模型改变后初始应力变化 new gen zone brick size 5 5 5 model elas prop sh 2e8 bu 3e8 fix x y z range z -.1 .1 set grav 0 0 -10 ini dens 2000 Solve ;网格三维显示 plot set rotation 30 0 30 plot cont zdisp outl on ;模型改变 model mohr range x 0 2 y 0 5 z 0 2 prop sh 2e8 bu 3e8 fric 35 range x 0 2 y 0 5 z 0 2 Step 2000 ;网格三维显示 plot set rotation 30 0 30 plot cont zdisp outl on 4-20应力与孔隙压力的初始化 new config fluid gen zone brick size 8 5 10 model elas model fl_iso ini dens 2000 prop bulk 1e9 shear 5e8 prop poros 0.5 perm 1e-10 ini fmod 2e9 ini fdensity 1e3 ini sat 0 ini sat 1 range z -.1 5.1 set grav 0 0 -10 fix x range x -.1 .1 fix x range x 7.9 8.1 fix y range y -.1 .1 fix y range y 4.9 5.1 fix z range z -.1 .1 ini pp 5.e4 grad 0,0,-1.e4 range z 0.0 5. ini szz -20e4 grad 0,0,20e3 ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5 solve ;网格三维显示 plot set rotation 30 0 30 plot cont zdisp outl on new gen zone brick size 1 1 10 model elas ini dens 2500 range z 0 5 ini dens 2250 range z 5 6 ini dens 2000 range z 6 10 prop bulk 1e9 shear 5e8 set grav 10 water dens 1000 water table ori 0 0 5 normal 0 0 1 fix x range x -.1 .1 fix x range x 7.9 8.1 fix y range y -.1 .1 fix y range y 4.9 5.1 fix z range z -.1 .1 ini szz -20e4 grad 0,0,20e3 ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5 solve ;网格三维显示 plot set rotation 30 0 30 plot cont zdisp outl on 4-21加载顺序 new ; 主隧道建模 gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 & size 4 10 6 4 dim 4 4 4 4 rat 1 1 1 1 fill gen zon reflect dip 90 dd 90 orig 15 0 0 gen zon reflect dip 0 ori 0 0 0 ;辅助隧道建模 gen point id 1 (2.969848,0.0,-0.575736) gen point id 2 (2.969848,50.0,-0.575736) gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 & p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fill gen zon radcyl p0 0 0 -1 p1 0 0 -8 p2 0 50 -1 p3 7 0 0 p4 0 50 -8 & p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fill ;隧道边界网格 gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 & p5 15 50 50 p6 50 0 50 p7 50 50 50 & p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 & size 6 10 3 10 rat 1 1 1 1.1 gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 & p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 & size 3 10 10 rat 1 1 1.1 gen zon reflect dip 0 ori 0 0 0 range x 0 23 y 0 50 z 8 50 gen zon reflect dip 0 ori 0 0 0 range x 23 50 y 0 50 z 0 50 group main1 range cyl end1 15 0 0 end2 15 25 0 rad 4 group main2 range cyl end1 15 25 0 end2 15 50 0 rad 4 group service range cyl end1 0.0,0.0,-0.575736 & end2 0.0,50.0,-0.575736 rad 3.0 save tun0.sav rest tun0.sav ; 初始应力场模量 model mohr prop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5 apply szz -1.4e6 range z 49.9 50.1 fix z range z -50.1 -49.1 fix x range x -.1 .1 fix x range x 49.9 50.1 fix y range y -.1 .1 fix y range y 49.9 50.1 ini sxx -1.4e6 syy -1.4e6 szz -1.4e6 hist unbal hist gp xdis 3,0,-1 hist gp zdis 0,0,2 hist gp xdis 3,25,-1 hist gp zdis 0,25,2 step 1000 save tun1.sav ; 辅助隧道开挖25 m计算 rest tun1.sav ini xdis 0.0 ydis 0.0 zdis 0.0 model null range group service y 0,25 step 1000 save tun2.sav ; 在辅助隧道上施加衬砌 restore tun2.sav sel shell id=1 range cyl end1 0 0 -1 end2 0 25 -1 rad 3 sel shell prop iso=(25.3e9, 0.266) thick = 0.5 sel node fix y xr zr range y -0.1 0.1 ; symmetry cond. sel node fix x yr zr range x -0.1 0.1 ; symmetry cond. model mohr range group service y 25 50 prop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5 ; 主隧道开挖25m model null range group main1 step 1000 save tun3.sav plot set rotation 0 0 60 plot cont zdisp outl on 6-1 V级围岩施工过程模拟 （1）建立隧道1/4圆周模型 gen zon radcyl p0 0 0 0 p1 10.55 0 0 p2 0 1 0 p3 0 0 10.55 & size 5 2 10 4 dim 5.55 5.55 5.55 5.55 rat 1 1 1 1.2 group outsiderocksy plot block group plot add axes red gen zone cshell p0 0 0 0 p1 5.55 0 0 p2 0 1 0 p3 0 0 5.55 & size 1 2 10 4 dim 5.05 5.05 5.05 5.05 rat 1 1 1 1 group concretlinersy fill group insiderocksy gen zon radcyl p0 0 0 0 p1 0 0 -9 p2 0 1 0 p3 10.55 0 0 & size 5 2 10 4 dim 4 5.55 4 5.55 rat 1 1 1 1.2 group outsiderockxy gen zone cshell p0 0 0 0 p1 0 0 -4 p2 0 1 0 p3 5.55 0 0 & size 1 2 10 4 dim 3.5 5.05 3.5 5.05 rat 1 1 1 1 group concretlinerxy fill group insiderockxy plot block group plot add axes red gen zon brick p0 0 0 10.55 p1 10.55 0 10.55 p2 0 1 10.55 p3 0 0 25.55 & size 5 2 8 group outsiderock1 gen zon brick p0 0 0 -34.55 p1 10.55 0 -34.55 p2 0 1 -34.55 p3 0 0 -9 & size 5 2 12 group outsiderock2 gen zon brick p0 10.55 0 0 p1 50 0 0 p2 10.55 1 0 p3 10.55 0 10.55 & size 20 2 5 group outsiderock3 gen zon brick p0 10.55 0 -9 p1 50 0 -9 p2 10.55 1 -9 p3 10.55 0 0 & size 20 2 5 group outsiderock3 gen zon brick p0 10.55 0 10.55 p1 50 0 10.55 p2 10.55 1 10.55 p3 10.55 0 25.55 & size 20 2 8 group outsiderock4 gen zon brick p0 10.55 0 -34.55 p1 50 0 -34.55 p2 10.55 1 -34.55 p3 10.55 0 -9 & size 20 2 12 group outsiderock5 gen zon reflect dip 90 dd 90 orig 0 0 0 plot set rotation 30 0 30 ;显示三维图，绕x轴负方向转30度 set gravity 0 0 -10 ;设置重力加速度为z方向-10 fix z range z -34.56,-34. ;设置底边界 fix x range x -50.01,-49.99 ;设置左边界 fix x range x 49.99,50.01 ;设置右边界 fix y range y -0.01 0.01 ;设置前边界 fix y range y 0.99,1.01 ;设置后边界 model mohr ;莫尔～库仑模型 ini density 2000 ;围岩的密度 prop bulk=7.14e8 shear=3.333e9 fric=25 coh=2e5 tension=1.0e5 Setp 6000 ;求解6000次 Solve ;求解计算 save Gsol.sav ;计算结果保存在Gsol.sav文件中 plot cont zdisp plot cont xdisp ;绘制竖向位移场，如图6-14所示 ;绘制竖向位移场，如图6-15所示 plot cont szz ;绘制竖向应力场，如图6-16所示 plot cont sxx ;绘制竖向应力场，如图6-17所示 ;1-左上半断面开挖 rest Gsol.sav plot block group plot add axes red ini xdis 0.0 ydis 0.0 zdis 0.0 ini density 2200 range group outsiderocksy any group outsiderocksz any prop bulk=1.923e9, shear=4.167e9 fric=35, coh=0.5e6 & tension=5e5, range group outsiderocksy any group outsiderocksz any model null range group insiderocksz any group concretlinersz any step 2000 save stepp1.sav plot cont zdisp rest stepp1.sav sel shell id=1 range x -0.01 0.01 z -0.01 5.56 cyl end1 0 0 0 end2 0 1 0 rad 5.55 sel shell id=1 range x -0.01 0.01 z -0.01 5.56 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25 step 2000 plot cont zdisp save stepp2.sav ;3-左下半断面开挖 rest stepp2.sav plot block group plot add axes red model null range group insiderockxz any group concretlinerxz any step 2000 save stepp3.sav plot cont zdisp rest stepp3.sav sel shell id=1 range x -0.01 0.01 z -4.01 0.01 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25 step 2000 plot cont zdisp save stepp4.sav ;5-右上半断面开挖 rest stepp4.sav plot block group plot add axes red model null range group insiderocksy any group concretlinersy any step 2000 save stepp5.sav plot cont zdisp rest stepp5.sav sel shell id=1 range x -0.01 5.56 z -0.01 5.56 cyl end1 0 0 0 end2 0 1 0 rad 5.55 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25 step 2000 plot cont zdisp save stepp6.sav ;7-右下半断面开挖 rest stepp6.sav plot block group plot add axes red model null range group insiderockxy any group concretlinerxy any step 2000 save stepp7.sav plot cont zdisp rest stepp7.sav sel delete shell id=1 range x -0.01 0.01 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25 step 2000 save stepp8.sav rest stepp8.sav model elas range group concretlinerxy model elas range group concretlinerxz model elas range group concretlinersy model elas range group concretlinersz ini density 2600 range group concretlinerxy ini density 2600 range group concretlinerxz ini density 2600 range group concretlinersy ini density 2600 range group concretlinersz prop bulk=12.5e9, shear=16.667e9 range group concretlinerxy prop bulk=12.5e9, shear=16.667e9 range group concretlinerxz prop bulk=12.5e9, shear=16.667e9 range group concretlinersy prop bulk=12.5e9, shear=16.667e9 range group concretlinersz step 2000 plot cont zdisp save stepp9.sav plot block group plot add axes red 6.2 IV级围岩施工过程模拟 new ;1/4 上半部分 gen zon radcyl p0 0 0 0 p1 10.55 0 0 p2 0 1 0 p3 0 0 10.55 size 5 2 10 4& dim 5.55 5.55 5.55 5.55 rat 1 1 1 1.2 group outsiderocksy gen zone cshell p0 0 0 0 p1 5.55 0 0 p2 0 1 0 p3 0 0 5.55 size 1 2 10 4& dim 5.15 5.15 5.15 5.15 rat 1 1 1 1 group concretlinersy fill group insiderocksy ;1/4 下半部分 gen zon radcyl p0 0 0 0 p1 0 0 -9 p2 0 1 0 p3 10.55 0 0 & size 5 2 10 4 dim 4 5.55 4 5.55 rat 1 1 1 1.2 group outsiderockxy gen zone cshell p0 0 0 0 p1 0 0 -4 p2 0 1 0 p3 5.55 0 0 size 1 2 10 4& dim 3.6 5.15 3.6 5.15 rat 1 1 1 1 group concretlinerxy fill group insiderockxy ;1/2绘制上下地层网格 gen zon brick p0 0 0 10.55 p1 10.55 0 10.55 p2 0 1 10.55 p3 0 0 25.55 size 5 2 8 group outsiderock1 gen zon brick p0 0 0 -34.55 p1 10.55 0 -34.55 p2 0 1 -34.55 p3 0 0 -9 size 5 2 12 group outsiderock2 ;1/2右侧中间部分土体网格 gen zon brick p0 10.55 0 0 p1 50 0 0 p2 10.55 1 0 p3 10.55 0 10.55 size 20 2 5 group outsiderock3 gen zon brick p0 10.55 0 -9 p1 50 0 -9 p2 10.55 1 -9 p3 10.55 0 0 size 20 2 5 group outsiderock3 ;1/2右上侧和下侧土体网格 gen zon brick p0 10.55 0 10.55 p1 50 0 10.55 p2 10.55 1 10.55 p3 10.55 0 25.55 size 20 2 8 group outsiderock4 gen zon brick p0 10.55 0 -34.55 p1 50 0 -34.55 p2 10.55 1 -34.55 p3 10.55 0 -9 size 20 2 12 group outsiderock5 ;对称得出完整模型 gen zon reflect dip 90 dd 90 orig 0 0 0 ;设置重力加速度 set gravity 0 0 -10 ;给定边界条件 fix z range z -34.56,-34. fix x range x -50.01,-49.99 fix x range x 49.99,50.01 fix y range y -0.01 0.01 fix y range y 0.99,1.01 ;0-求解自重应力场 model mohr ini density 2200 ;围岩的密度 prop bulk=1.923e9, shear=4.167e9 fric=35, coh=0.5e6 tension=5e5 ;体积、剪切、摩擦角、凝聚力、抗拉强度 step 6000 save Gsol.sav ;绘制计算结果 plot cont zdisp plot cont xdisp plot cont szz plot cont sxx ;new ;1-上半断面开挖 rest Gsol.sav ini xdis 0.0 ydis 0.0 zdis 0.0 model null range group insiderocksz any group concretlinersz any & group insiderocksy any group concretlinersy any step 2000 save stepp1.sav plot cont zdisp ;2-初期支护 rest stepp1.sav sel shell id=1 range x -5.56 5.56 z -0.01 5.56 cyl end1 0 0 0 end2 0 1 0 rad 5.55 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.20 step 20 save stepp2.sav ;3-下半断面开挖 rest stepp2.sav model null range group insiderockxz any group concretlinerxz any & group insiderockxy any group concretlinerxy any step 2000 save stepp3.sav plot cont zdisp ;4-二次衬砌 rest stepp3.sav model elas range group concretlinerxy model elas range group concretlinerxz model elas range group concretlinersy model elas range group concretlinersz ini density 2600 range group concretlinerxy ini density 2600 range group concretlinerxz ini density 2600 range group concretlinersy ini density 2600 range group concretlinersz prop bulk=12.5e9, shear=16.667e9 range group concretlinerxy prop bulk=12.5e9, shear=16.667e9 range group concretlinerxz prop bulk=12.5e9, shear=16.667e9 range group concretlinersy prop bulk=12.5e9, shear=16.667e9 range group concretlinersz step 2000 save stepp4.sav 6.3 III级围岩施工过程模拟 new ;1/4 上半部分 gen zon radcyl p0 0 0 0 p1 10.55 0 0 p2 0 1 0 p3 0 0 10.55 & size 5 2 10 4 dim 5.55 5.55 5.55 5.55 rat 1 1 1 1.2 group outsiderocksy gen zone cshell p0 0 0 0 p1 5.55 0 0 p2 0 1 0 p3 0 0 5.55 size 1 2 10 4 & dim 5.25 5.25 5.25 5.25 rat 1 1 1 1 group concretlinersy fill group insiderocksy ;1/4 下半部分 gen zon radcyl p0 0 0 0 p1 0 0 -9 p2 0 1 0 p3 10.55 0 0 & size 5 2 10 4 dim 4 5.55 4 5.55 rat 1 1 1 1.2 group outsiderockxy gen zone cshell p0 0 0 0 p1 0 0 -4 p2 0 1 0 p3 5.55 0 0 size 1 2 10 4 & dim 3.7 5.25 3.7 5.25 rat 1 1 1 1 group concretlinerxy fill group insiderockxy ;1/2绘制上下地层网格 gen zon brick p0 0 0 10.55 p1 10.55 0 10.55 p2 0 1 10.55 p3 0 0 25.55 size 5 2 8 group outsiderock1 gen zon brick p0 0 0 -34.55 p1 10.55 0 -34.55 p2 0 1 -34.55 p3 0 0 -9 size 5 2 12 group outsiderock2 ;1/2右侧中间部分土体网格 gen zon brick p0 10.55 0 0 p1 50 0 0 p2 10.55 1 0 p3 10.55 0 10.55 size 20 2 5 group outsiderock3 gen zon brick p0 10.55 0 -9 p1 50 0 -9 p2 10.55 1 -9 p3 10.55 0 0 size 20 2 5 group outsiderock3 ;1/2右上侧和下侧土体网格 gen zon brick p0 10.55 0 10.55 p1 50 0 10.55 p2 10.55 1 10.55 p3 10.55 0 25.55 size 20 2 8 group outsiderock4 gen zon brick p0 10.55 0 -34.55 p1 50 0 -34.55 p2 10.55 1 -34.55 p3 10.55 0 -9 size 20 2 12 group outsiderock5 ;对称得出完整模型 gen zon reflect dip 90 dd 90 orig 0 0 0 ;设置重力加速度 set gravity 0 0 -10 ;给定边界条件 fix z range z -34.56,-34. fix x range x -50.01,-49.99 fix x range x 49.99,50.01 fix y range y -0.01 0.01 fix y range y 0.99,1.01 ;0-求解自重应力场 model mohr ini density 2400 ;围岩的密度 prop bulk=6.0e9, shear=10.0e9 fric=45, coh=1.2e6 tension=2e6 ;体积、剪切、摩擦角、凝聚力、抗拉强度 step 6000 save Gsol.sav ;1-全断面开挖 rest Gsol.sav ini xdis 0.0 ydis 0.0 zdis 0.0 model null range group insiderocksz any group concretlinersz any & group insiderocksy any group concretlinersy any model null range group insiderockxz any group concretlinerxz any & group insiderockxy any group concretlinerxy any step 2000 save stepp1.sav ;2-初期支护 rest stepp1.sav sel shell id=1 range group x -5.56 5.56 z -0.01 5.56 cyl end1 0 0 0 end2 0 1 0 rad 5.55 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.15 step 2000 save stepp2.sav ;3-二次衬砌 rest stepp2.sav model elas range group concretlinerxy model elas range group concretlinerxz model elas range group concretlinersy model elas range group concretlinersz ini density 2600 range group concretlinerxy ini density 2600 range group concretlinerxz ini density 2600 range group concretlinersy ini density 2600 range group concretlinersz prop bulk=12.5e9, shear=16.667e9 range group concretlinerxy prop bulk=12.5e9, shear=16.667e9 range group concretlinerxz prop bulk=12.5e9, shear=16.667e9 range group concretlinersy prop bulk=12.5e9, shear=16.667e9 range group concretlinersz step 2000 save stepp3.sav plot block group new 第七章 命令流按照顺序进行 new ;1/4 上半部分 gen zon radcyl p0 15.7 0 0 p1 24.5 0 0 p2 15.7 1 0 p3 15.7 0 8.8 & size 5 2 10 4 dim 5.8 5.8 5.8 5.8 rat 1 1 1 1.2 group outrockys gen zone cshell p0 15.7 0 0 p1 21.5 0 0 p2 15.7 1 0 p3 15.7 0 5.8 & size 1 2 10 4 dim 5.4 5.4 5.4 5.4 rat 1 1 1 1 group liningys fill group inrockys gen zon reflect dip 90 dd 270 orig 15.7 0 0 ;1/4 下半部分 gen zon radcyl p0 15.7 0 0 p1 15.7 0 -6.9 p2 15.7 1 0 p3 24.5 0 0 & size 5 2 10 4 dim 3.9 5.8 3.9 5.8 rat 1 1 1 1.2 group outrockyx gen zone cshell p0 15.7 0 0 p1 15.7 0 -3.9 p2 15.7 1 0 p3 21.5 0 0 & size 1 2 10 4 dim 3.5 5.4 3.5 5.4 rat 1 1 1 1 group liningyx fill group inrockyx gen zon reflect dip 90 dd 270 orig 15.7 0 0 range group liningyx any & group inrockyx any group outrockyx any ;右隧道中地层网格 gen zon brick p0 0 0 -6.9 p1 6.9 0 -6.9 p2 0 1 -6.9 p3 0 0 0 size 6 2 5 group outrock1 gen zon brick p0 24.5 0 -6.9 p1 65 0 -6.9 p2 24.5 1 -6.9 p3 24.5 0 0 size 20 2 5 group outrock1 gen zon brick p0 0 0 0 p1 6.9 0 0 p2 0 1 0 p3 0 0 8.8 size 6 2 5 group outrock2 gen zon brick p0 24.5 0 0 p1 65 0 0 p2 24.5 1 0 p3 24.5 0 8.8 size 20 2 5 group outrock2 ;右隧道上侧土体网格 gen zon brick p0 0 0 8.8 p1 6.9 0 8.8 p2 0 1 8.8 p3 0 0 36 size 6 2 10 group outrock3 gen zon brick p0 24.5 0 8.8 p1 65 0 8.8 p2 24.5 1 8.8 p3 24.5 0 36 size 20 2 10 group outrock3 gen zon brick p0 6.9 0 8.8 p1 24.5 0 8.8 p2 6.9 1 8.8 p3 6.9 0 36 size 10 2 10 group outrock3 ;右隧道下侧土体网格 gen zon brick p0 0 0 -24 p1 6.9 0 -24 p2 0 1 -24 p3 0 0 -6.9 size 6 2 8 group outrock4 gen zon brick p0 24.5 0 -24 p1 65 0 -24 p2 24.5 1 -24 p3 24.5 0 -6.9 size 20 2 8 groupoutrock4 gen zon brick p0 6.9 0 -24 p1 24.5 0 -24 p2 6.9 1 -24 p3 6.9 0 -6.9 size 10 2 8 group outrock4 gen zon reflect dip 90 dd 270 orig 0 0 0 plot add axes red plot block group ;设置重力加速度 set gravity 0 0 -10 ;给定边界条件 fix z range z -24.01,-23.99 fix x range x -65.01,-.99 fix x range x .99,65.01 fix y range y -0.01 0.01 fix y range y 0.99,1.01 ;0-求解自重应力场 model mohr ini density 2200 ;围岩的密度 prop bulk=1.923e9, shear=4.167e9 fric=35, coh=0.5e6 tension=5e5 ;体积、剪切、摩擦角、凝聚力、抗拉强度 step 6000 save Gsol.sav plot cont zdisp plot cont szz ;1-左隧道上半断面开挖 rest Gsol.sav ini xdis 0.0 ydis 0.0 zdis 0.0 ini density 2400 range group outrockzs any prop bulk=6.0e9, shear=10.0e9 fric=45, coh=1.2e6 tension=2e6 & range group outrockzs any model null range group inrockzs any group liningzs any step 2000 save stepp1.sav ;2-初期支护 rest stepp1.sav sel shell id=1 range x -21.51 -9. cyl end1 -15.7 0 0 end2 -15.7 1 0 rad 5.8 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.20 step 2000 save stepp2.sav ;3-左隧道下半断面开挖 rest stepp2.sav model null range group inrockzx any group liningzx any step 2000 save stepp3.sav plot add axes red plot block group plot cont zdisp plot cont szz ;4-左隧道二次衬砌 rest stepp3.sav model elas range group liningzs model elas range group liningzx ini density 2600 range group liningzs ini density 2600 range group liningzx prop bulk=12.5e9, shear=16.667e9 range group liningzs prop bulk=12.5e9, shear=16.667e9 range group liningzx step 2000 save stepp4.sav plot add axes red plot block group plot cont zdisp plot cont szz ;5-右隧道上半断面开挖 rest stepp4.sav ini density 2400 range group outrockys any prop bulk=6.0e9, shear=10.0e9 fric=45, coh=1.2e6 tension=2e6 & range group outrockys any model null range group inrockys any group liningys any step 2000 save stepp5.sav plot cont zdisp ;6-初期支护 rest stepp5.sav sel shell id=1 range x 9. 21.51 cyl end1 15.7 0 0 end2 15.7 1 0 rad 5.8 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.20 step 2000 plot cont zdisp save stepp6.sav plot add axes red plot block group plot cont zdisp plot cont szz ;7-左隧道下半断面开挖 rest stepp6.sav model null range group inrockyx any group liningyx any step 2000 save stepp7.sav plot cont zdisp ;8-左隧道二次衬砌 rest stepp7.sav model elas range group liningys model elas range group liningyx ini density 2600 range group liningys ini density 2600 range group liningyx prop bulk=12.5e9, shear=16.667e9 range group liningys prop bulk=12.5e9, shear=16.667e9 range group liningyx step 2000 save stepp8.sav plot add axes red plot block group plot cont zdisp plot cont szz new 第七章 命令流按照顺序进行 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 model elas range group linz y 0 6 ini density 2500 range group linz y 0 6 prop bulk=16.67e9,shear=12.5e9 range group linz y 0 6 step 3000 save step1.sav ;第2步 rest step1.sav model null range group inrz y 6 12 model elas range group linz y 6 12 ini density 2500 range group linz y 6 12 prop bulk=16.67e9,shear=12.5e9 range group linz y 6 12 step 3000 save step2.sav ;第3步 rest step2.sav model null range group inrz y 12 18 model elas range group linz y 12 18 ini density 2500 range group linz y 12 18 prop bulk=16.67e9,shear=12.5e9 range group linz y 12 18 step 3000 save step3.sav ;第4步 rest step3.sav model null range group inrz y 18 24 model elas range group linz y 18 24 ini density 2500 range group linz y 18 24 prop bulk=16.67e9,shear=12.5e9 range group linz y 18 24 step 3000 save step4.sav ;第5步 rest step4.sav model null range group inrz y 24 30 model elas range group linz y 24 30 ini density 2500 range group linz y 24 30 prop bulk=16.67e9,shear=12.5e9 range group linz y 24 30 step 3000 save step5.sav ;第6步 rest step5.sav model null range group inrz y 30 36 model elas range group linz y 30 36 ini density 2500 range group linz y 30 36 prop bulk=16.67e9,shear=12.5e9 range group linz y 30 36 step 3000 save step6.sav ;第1步 rest step6.sav model null range group inry y 0 6 model elas range group liny y 0 6 ini density 2500 range group liny y 0 6 prop bulk=16.67e9,shear=12.5e9 range group liny y 0 6 step 3000 save step7.sav ;第2步 rest step7.sav model null range group inry y 6 12 model elas range group liny y 6 12 ini density 2500 range group liny y 6 12 prop bulk=16.67e9,shear=12.5e9 range group liny y 6 12 step 3000 save step8.sav ;第3步 rest step8.sav model null range group inry y 12 18 model elas range group liny y 12 18 ini density 2500 range group liny y 12 18 prop bulk=16.67e9,shear=12.5e9 range group liny y 12 18 step 3000 save step9.sav ;第4步 rest step9.sav model null range group inry y 18 24 model elas range group liny y 18 24 ini density 2500 range group liny y 18 24 prop bulk=16.67e9,shear=12.5e9 range group liny y 18 24 step 3000 save step10.sav ;第5步 rest step10.sav model null range group inry y 24 30 model elas range group liny y 24 30 ini density 2500 range group liny y 24 30 prop bulk=16.67e9,shear=12.5e9 range group liny y 24 30 step 3000 save step11.sav ;第6步 rest step11.sav model null range group inry y 30 36 model elas range group liny y 30 36 ini density 2500 range group liny y 30 36 prop bulk=16.67e9,shear=12.5e9 range group liny y 30 36 step 3000 save step12.sav new 因篇幅问题不能全部显示，请点此查看更多更全内容
Copyright © 2019- efsc.cn 版权所有 赣ICP备2024042792号-1
违法及侵权请联系：TEL:199 1889 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务