隧道及地下工程FLAC解析方法常识、建模、常用命令流及其解释

1、隧道及地下工程FLAC解析方法 _常识、建模、常用命令流及其解释 TOC o 1-1 h z u HYPERLINK l _Toc320528144 1.1 FLAC常识 PAGEREF _Toc320528144 h 2 HYPERLINK l _Toc320528145 1.2常用命令流 PAGEREF _Toc320528145 h 13 HYPERLINK l _Toc320528146 1.3建模过程 PAGEREF _Toc320528146 h 21 HYPERLINK l _Toc320528147 2-1定义一个FISH函数 PAGEREF _Toc320528147 h 2

2、4 HYPERLINK l _Toc320528148 2-2使用一个变量 PAGEREF _Toc320528148 h 24 HYPERLINK l _Toc320528149 2-3对变量和函数的理解 PAGEREF _Toc320528149 h 24 HYPERLINK l _Toc320528150 2-4获取变量的历史记录 PAGEREF _Toc320528150 h 24 HYPERLINK l _Toc320528151 2-5用FISH函数计算体积模量和剪砌模量 PAGEREF _Toc320528151 h 25 HYPERLINK l _Toc320528152 2-

3、6 在FLAC输入中使用符号变量 PAGEREF _Toc320528152 h 25 HYPERLINK l _Toc320528153 2-7 控制循环 PAGEREF _Toc320528153 h 26 HYPERLINK l _Toc320528154 2-8 拆分命令行 PAGEREF _Toc320528154 h 26 HYPERLINK l _Toc320528155 2-9 变量类型 PAGEREF _Toc320528155 h 27 HYPERLINK l _Toc320528156 2-10 IF条件语句 PAGEREF _Toc320528156 h 27 HYPE

4、RLINK l _Toc320528157 2-11 索单元自动生成 PAGEREF _Toc320528157 h 27 HYPERLINK l _Toc320528158 2-12圆形隧道开挖模拟计算 PAGEREF _Toc320528158 h 28 HYPERLINK l _Toc320528159 4-1数组 PAGEREF _Toc320528159 h 30 HYPERLINK l _Toc320528160 4-2函数操作 PAGEREF _Toc320528160 h 31 HYPERLINK l _Toc320528161 4-3函数删除与重定义 PAGEREF _Toc

5、320528161 h 32 HYPERLINK l _Toc320528162 4-4字符串 PAGEREF _Toc320528162 h 32 HYPERLINK l _Toc320528163 4-5马蹄形隧道网格 PAGEREF _Toc320528163 h 33 HYPERLINK l _Toc320528164 4-6复杂形状网格生成 PAGEREF _Toc320528164 h 33 HYPERLINK l _Toc320528165 4-7网格连接 PAGEREF _Toc320528165 h 34 HYPERLINK l _Toc320528166 4-8立方体洞穴网

6、格生成 PAGEREF _Toc320528166 h 35 HYPERLINK l _Toc320528167 4-9球体洞穴网格生成 PAGEREF _Toc320528167 h 36 HYPERLINK l _Toc320528168 4-10应力边界 PAGEREF _Toc320528168 h 37 HYPERLINK l _Toc320528169 4-11改变应力边界条件 PAGEREF _Toc320528169 h 37 HYPERLINK l _Toc320528170 4-12位移边界 PAGEREF _Toc320528170 h 39 HYPERLINK l _T

7、oc320528171 4-13不考虑重力影响的均匀应力 PAGEREF _Toc320528171 h 41 HYPERLINK l _Toc320528172 4-14考虑应力梯度的均匀材料 PAGEREF _Toc320528172 h 41 HYPERLINK l _Toc320528173 4-15考虑应力梯度的非均匀材料 PAGEREF _Toc320528173 h 42 HYPERLINK l _Toc320528174 4-16非均匀网格应力初始化 PAGEREF _Toc320528174 h 42 HYPERLINK l _Toc320528175 4-17不规则自由面应

8、力初始化 PAGEREF _Toc320528175 h 43 HYPERLINK l _Toc320528176 4-18非均网格内部压实 PAGEREF _Toc320528176 h 43 HYPERLINK l _Toc320528177 4-19模型改变后初始应力变化 PAGEREF _Toc320528177 h 44 HYPERLINK l _Toc320528178 4-20应力与孔隙压力的初始化 PAGEREF _Toc320528178 h 44 HYPERLINK l _Toc320528179 4-21加载顺序 PAGEREF _Toc320528179 h 45 HY

9、PERLINK l _Toc320528180 6-1 V级围岩施工过程模拟 PAGEREF _Toc320528180 h 47 HYPERLINK l _Toc320528181 6.2 IV级围岩施工过程模拟 PAGEREF _Toc320528181 h 50 HYPERLINK l _Toc320528182 6.3 III级围岩施工过程模拟 PAGEREF _Toc320528182 h 52 HYPERLINK l _Toc320528183 第七章 命令流按照顺序进行 PAGEREF _Toc320528183 h 54 HYPERLINK l _Toc320528184 ;右

10、隧道中地层网格 PAGEREF _Toc320528184 h 54 HYPERLINK l _Toc320528185 ;右隧道上侧土体网格 PAGEREF _Toc320528185 h 54 HYPERLINK l _Toc320528186 ;右隧道下侧土体网格 PAGEREF _Toc320528186 h 54 HYPERLINK l _Toc320528187 第七章 命令流按照顺序进行 PAGEREF _Toc320528187 h 571.1 FLAC常识1. FLAC3D是有限元程序吗？答：不是！是有限差分法。2. 最先需要掌握的命令有哪些？答：需要掌握gen, ini,

11、app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。3. 怎样看模型的样子？答：plo blo gro可以看到不同的group的颜色分布4. 怎样看模型的边界情况？答：plo gpfix red5. 怎样看模型的体力分布？答：plo fap red6. 怎样看模型的云图？答：位移：plo con dis (xdis, ydis, zdis) 应力：plo con sz (sy, sx, sxy, syz, sxz)7. 怎样看模型的矢量图？答：plo dis (xdis, ydis, zdis)8. 怎样看模型有多少单元、节点？答：pri info9. 怎样输出模型的

12、后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件10. 怎样调用一个文件？答：File/call或者call命令10. 如何施加面力？答：app nstress11. 如何调整视图的大小、角度？答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键12. 如何进行边界约束？答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束）13. 如何知道每个单元的ID？答：用鼠标双击单元的表面，可以知道单元的ID和坐标14. 如何进行切片？答：plo set plane ori

13、(点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 15. 如何保存计算结果？答：save +文件名.16. 如何调用已保存的结果？答：rest +文件名；或者File / Restore17. 如何暂停计算？答：Esc18. 如何在程序中进行暂停，并可恢复计算？答：在命令中加入pause命令，用continue进行继续19. 如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步20. Fish是什么东西？答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学？答：可以不用，需要的

14、时候查Mannual获得需要的变量就可以了22. FLAC3D允许的命令文件格式有哪些？答：无所谓，只要是文本文件，什么后缀都可以23. 如何调用一些可选模块？答：config dyn (fluid, creep, cppudm)24 .如何在圆柱体四周如何施加约束条件？答：可以用fix . ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2r1，其实就是选择range的时候选两个圆柱的差，即得到边界。25.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用set log on 命令，把数据导

15、出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。26.关于dd、dip的理解：答： 首先你要确定你定义面的法线方向，dip 是倾角的意思（xoy面和定义面的夹角），法线向上dip是锐角，法线向下dip 时钝角。dd是y轴顺时针转到法线在xoy面上的投影线时，所转过的角度。27.如何调用FISH文件？答：新手调用FISH文件出错的原因：是：文件名为 fishcall.fis.txt(因为在记事本里这个txt的后缀是默认的）。修改方法是：将 fishcall.fis.txt里面的内容全部复制到新的记事本里，再按保存，保存时候将保存类型改为所有文件(*,*),这时你便可

16、在文件名栏内输入 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的话就明白mirro

17、r的含义了30.如何查看剖面上任一点的位移值？答：print gp disp range x() y( ) z( )或print gp disp range id plot gp disp range x() y( ) z( )或plot gp disp range id31. 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 有何区别？答：初始化命令，如初始化计算体的应力状态等；

18、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 unbalstep 100000 or solvehis write 1 vs step begin

19、 1 end 1000 file 123.his使用上述命令就可以查看各个时步下的不平衡力的具体数值34.如何显示变形轮廓线的命令？ plo ske magf 10 其中10为放大系数35、怎样开始入门？答：初步学习首先看说明书 Users gurid的Section2 Getting started和Section 3 Problem sloving with FLAC36、哪些命令需要掌握？答：grid gen prop interface ini fix set hist step solve pause cont sav 等，在学习的过程中遇到这些命令和不会的命令应该查看说明书的Com

20、mand Reference，里面讲的很清楚，特别需要注意的是一定要仔细看一下Command Reference开头的Optional Range Phrases，这是讲范围的，因为很多命令的施加都需要加到模型的一部分，这里就告诉你怎样来表示这一部分。这一原则同样适合与其他5个ITASCA软件。还有一点就是六个软件（FLAC FALC3D UDEC 3DEC PFC PFC3D）的range大不相同。37、FLAC和FLAC3D建模有什么区别？答：FLAC和FLAC3D基本的建模思想大致是一样的，不同的是FLAC是首先用GRID命令把整体的网格画出来，然后用MODELL NULL命令把网格分成

21、许多部分，再用GEN命令把这些部分一块一块搭接起来。然后再用GEN LINE，GEN ARC等命令画出要开挖的边界，在FLAC中产生网格的命令只有一个GRID，而且只能用一次。而FLAC3D是直接用GEN 命令一块一块搭接，GEN可以用很多次。38、怎样查看模型？答：plot grid 可以查看网格，plot grid num 可以查看节点号。39、怎样查看边界条件？答：plot fix40、怎样查看预先留设的开挖边界？答：plot grid mark。41、怎么查看模型的云图？答：pl syy (sxx xdis ydis sig1.) fill ，注意，带fill关键字就是查看云图，不带f

22、ill 关键字就是查看等直线图。42、怎样查看模型的矢量图？答：pl grid vel (disp) red43、怎样查看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中，当改变

23、应力边界时，模型表面的应力直接被更换，FLAC中改变应力边界时，模型表面的应力为两者应力相加。47、flac中影响运行时间的因素（1）单元材料，结构面之间的刚度对比很大；（2）划分单元尺寸对比很大，尽量避免长细比大于5：1的细长单元。48 FLAC中有多少中材料模型内置11中本构模型（1）零模型 用于开挖回填的区域（2）各向同性弹性模型（3）正交各向异性弹性模型（4）横观各向同性弹性模型（5）德鲁克普拉格塑性模型（6）莫尔库仑塑性模型 一般岩土体的力学行为（7）节理化塑性模型 用于研究材料破坏后行为（如边破稳定，屈服煤矿开采）（8）应变硬化/软化模型库仑塑性模型 薄板层状材料破坏后力学行为（9

24、）双线性应变硬化/软化节理化塑性模型 临近层状土层开挖（10）双屈服塑性模型（11）修正的剑桥粘土模型49、对不平衡力的理解flac中用最大不平衡力和典型内力大的比值R表示模型的不平衡力的相对大小，该值为百分数，且从不会减小到零，因此，模型不会达到绝对稳定的平横状态，可以根据要求的精度设定R值，模型默认的R值是1e-5。注意R值小仅表示所有网格点的力都平衡，但模型还可能发生塑性流动，通过 plot vel 命令查看网格点速度，一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。50、简单命令解释：plot add cont disp out on 时时显示云图plot block gr

25、oup 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 plo

26、t 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 boundaryfix x range x -0.1 0.1表示在x=0 这个平面只x方向不能发生运动roller boundary51、fish命令的解释：def aux1 定义函数loop n(1,nptab1) 开始一个循环，从1循环到变量nptab1rr=xtable

27、(1,n) 从table 1 的x列中读第n个数付给变量rrzz=ytable(1,n) 从table 1 的y列中读第n个数付给变量zzloop k(1,nprof) 开始第二个循环，从1循环到变量nprofcase_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 当k1的时候执行下面语句xtable(n+10,k)=rr table n+10的x列第k个值等于变量rrytable(n

28、+10,k)=-40.0 table n+10的y列第k个值等于变量-40case 2 当k2的时候执行下面语句xtable(n+10,k)=rr table n+10的x列第k个值等于变量rrytable(n+10,k)=0.0 table n+10的y列第k个值等于变量0end_case 结束选择结构end_loop 结束第二个循环end_loop 结束第一个循环end 函数定义结束aux1 执行什么定义的函数52. 在固定面上施加位移时，fix与apply不能同时使用，应该用ini和fix,才能施加上。53.菜单 plot items/add/structure elements/geo

29、metry/结构单元的形式，把对话框里的mag fac的值改成不同的数，就可以看到位移放大后的结构单元。54.你可以用键盘上的上、下、左、右肩在窗口中移动模型。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 在进行三轴压缩时候，如何对轴向力进行位

30、移控制？答：使用 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 int

31、er i nstress /ndis/sstress/sdisp61 接触面设置参数包括哪几个？答：coh dial fric kn ks ten.62 如何输出单元格的数据？答；expgrid .flac3d63 如何在剖面的沉降云图上增加网格、剖面上的矢量及坐标值？ D Q Z9s1 e4Y U%q 答：plo add ske plo add dis planeplo add axe64 如何做切片？ 答：plot set plane ori x y z norm x y z65 举一个隧道支护开挖的例子？gen zone radbrick p0 (0,0,0) p1 (10,0,0) p

32、2 (0,10,0) p3 (0,0,10) &size 3 5 5 7 rat 1 1 1 1.5 dim 1 4 2 fillmodel mohrprop bulk 1e8 shear .3e8 fric 35prop coh 1e3 tens 1e3fix x range x -0.1 0.1fix z range z -0.1 0.1fix y range y 9.9 10.1set grav 0 10 0ini dens 1000 ini syy=0.0 grad 0,-10000,0ini sxx=0.0 grad 0,-5000,0ini szz=0.0 grad 0,-5000

33、,0apply sxx=0.0 grad 0,-5000,0 range x 9.9,10.1apply szz=0.0 grad 0,-5000,0 range z 9.9,10.1set largehist unbalsolvesave tr_eq.sav; monitor trench excavationini xdis=0 ydis=0 zdis=0hist gp xdisp 1,0,0hist 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

34、 0.8 1.6 z 0 2range name exc3 x 0 1 y 1.6 2.4 z 0 2range name exc4 x 0 1 y 2.4 3.2 z 0 2range name exc5 x 0 1 y 3.2 4.0 z 0 2; excavation step 1model null range exc1step 10000save exc1.sav; excavation step 2model null range exc2step 1000save exc2.sav; excavation step 3model null range exc3step 1000s

35、ave exc3.sav ; excavation step 4model null range exc4step 1000save exc4.sav; excavation step 5model null range exc5step 1000 save exc5.sav; install cable supportrest exc2.savsel cable beg 1.0,0.4,1.5 end 5.0,0.4,1.5 nseg 4sel cable beg 1.0,0.4,0.5 end 5.0,0.4,0.5 nseg 4sel cable beg 1.0,1.2,1.5 end

36、5.0,1.2,1.5 nseg 4sel cable beg 1.0,1.2,0.5 end 5.0,1.2,0.5 nseg 4sel cable prop emod 2e9 ytension 1e8 xcarea 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0model null range exc3step 1000save cab3.savsel cable beg 1.0,2.0,1.5 end 5.0,2.0,1.5 nseg 4sel cable beg 1.0,2.0,0.5 end 5.0,2.0,0.5 nseg 4sel cable prop e

37、mod 2e9 ytension 1e8 xcarea 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0model null range exc4;step 1000save cab4.savsel cable beg 1.0,2.8,1.5 end 5.0,2.8,1.5 nseg 4sel cable beg 1.0,2.8,0.5 end 5.0,2.8,0.5 nseg 4sel cable prop emod 2e9 ytension 1e8 xcarea 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0model null range e

38、xc5step 1000save cab5.savplot create TrenchSupport plot set rot 63.89, 334.35, 3.4plot set cent 2.893, 3.446, 3.899plot set size 4.136plot add cont disp out onplot add sel cable forceplot add axes black66 举一个三轴压缩的例子？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,

39、0,0gen zone reflect norm 0,0,19model mohr prop bulk 1.19e10 shear 1.1e10prop coh 2.72e5 fric 44 ten 2e5fix x y z range y -.1 .1fix x y z range y 1.9 2.1apply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius 1ini yvel 0 grad 0 -1e-4 0 range y 0 2sol1.FLAC程式使用前准备步骤步骤1：依比例画出所欲分析的资料在纸上画出地点的位置、地

40、层资料、并简单标示距离及深度资料。步骤2：换算输入资料成同一单位将现有地层资料，如 Density, Bulk modulus, Youngs modulus, tension, cohesion, friction Angle 等资料，换算成同一单 位。附注 ：需谨慎检查输入资料的单位，如因单位不同而造成过大或过小的值，将会造成 FLAC无法计算，而产生ERROR信息。 步骤3：应用公式简略计算应用公式或依据经验，简略算出FLAC输出资料的范围，以作为Debug及输出资料分析时验证。步骤4：建立x, y坐标与node i, j之间的关系在图上距离及深度的关系，建立x, y坐标系统，再由x,

41、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 坐标与网格的关系，但由 于产生的

42、网格坐标不易控制，将对其它的工作产生负面影响，故依使用前步骤4所建立的关系，将网格依其疏密程度需要的不同，实 际控制网格的坐标。 (2)设定材料性质：prop (3)设定外力：Set Grav, Apply Pressure, ini sxx, Syy (4)设定边界条件：fix, free (5)求起始的应力平衡：solve (6)储存：Save5 求工程的影响_求出区域内的应力分布情况后，再依工程的流程及步骤阶段执行各工程进行过程的影响，建议使用以下的步骤： (1)调出起初的应力平衡：re_ .sav (2)设定新的材料性质：model,prop (3)设定新的支撑性质：struct (4

43、)设定新的外力 (5)设定边界条件 (6)求工程时的应力平衡(7)储存1.2常用命令流;模型镜像gen zone radcylinder size 25 1 25 25gen zone reflect normal -1 0 0 origin x y z(面上一点）;沿X轴镜像,通过对称平面法线向量确定对称面gen zone reflect normal 0 0 -1 ;沿z轴镜像;绘图控制pl contour szz outline on;在模型中显示位移-应变曲线hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0pl his

44、-2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负); -表示对其值作mirror;对模型进行压缩实验的方法;即在模型两侧施加相反方向的速度ini yvel 1e-7 range y -.1 .1ini 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 ;gradient更精确;输入角度、弧度方法 pi=，90为90.0*degraddef set_va

45、lsptA = 25.0 * sin(pi/2);ptA=25.0ptB = 25.0 * cos( 60.0*degrad );ptB=12.5ptC=pi;ptC=3.1415926endset_valsprint 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 notplot add sel geom black black cid on scale=0.03sel node init zpos add -25.0;如何显示某一平面

46、plot create name_planeplot set plane origin 3 4 0 normal 1 0 0plot add cont disp plane behind shade onplot add sel geom black plot add axes redplot current name_planeplot show;如何施加法向应力并显示所施加的力apply nstress -3.0E7 range cylinder end1 0 -1 0 end2010 radius 39.9 not;柱体外表面施加法向力apply nstress -5.0E6 range

47、 cylinder end1 0 -1 0 end2010 radius 2.1;柱体内表面施加法向力plot grid green fap black;显示施加的力;获取某网格点的ID，坐标值def aaget_id=gp_id(gp_near(3,0,0)endaaprint get_idprint gp position range id 3;保存命令流文件set log onset logfile hua.log;文件存放的路径与调用的文本路径相同set log off;fish函数的说明与解释def find_add ;定义fish函数find_add head = null ;给h

48、ead赋值p_gp = gp_head ;第一个网格结点的指针赋给p_gploop_while p_gp # null ;当p_gp值不为null时作循环x_pos = gp_xpos(p_gp) ;将指针为p_gp的结点的x坐标值赋给x_posif 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 = n

49、ew ;将new值赋给headendifp_gp = gp_next(p_gp) ;将结点指针为p_gp的下一个结点的指针赋给p_gpendloop;结束循环end ;结束fish函数;实际上这个fish函数为满足条件（x坐标为10的）的zone的地址开辟一定的地址空间，各地址之间存在一定的联系;找到第一个符合条件的zone地址后，用new = get_mem(2) 从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址，并用;下面的mem(new) = head ，将第一次开辟的两个变量的第一个变量存储地址head（注意第一次headnull），第二个变量存储第一;个符合条件的zo

50、ne地址，并将第一个变量的地址赋予head（headnew），第一次循环结束；下次循环，同样开辟两个变量对象，第;一个变量对象记录上次循环开辟的第一个变量的地址，第二个记录第二个符合条件的zone地址，其余循环依次类推，这样子就建立了;一个符合条件的zone地址链条，方便以后使用。；FLAC后处理；输出JPG图片set plot jpgset plot quality 100plot hard file XX.jpg;输出TXT文件hist write X file xx.txt ；(x是历史记录ID） print con zd;输出单元应力和节点位移set log onset logfile

51、 xx.logprint zo stressprint gp disset log off；生成影片文件plot set rot 20 0 30plot con szz ou on magf 10plot add hist 1 set movie avi step 1 file xx.avimovie start solve;计算间隔命令 在solve命令之前set hist_rep 1 ；(步数）；检查模型赋值情况plot bcon prop bu;变形云图加网格 放大20 倍plot con zd ou on magf 20;变形网格plot sk magf 20;变形矢量图（加网格）pl

52、ot sk dis scale 0.07 ；(scale是箭头长度，默认0.05）；如何观察塑性区plot block stateplot block state sh_nplot block state sh_p;切片操作plot set plane ori 1 2 3 norm 0 0 -1plot con zd planeplot add dis plane plot add sk plot add ax;切片后方部分显示plot con zd plane behind shade on!显示某一个groupplot block group range group 1显示监测的数据plo

53、t hist 112 ;112为hist编号FLAC有用的命令流显示命令流位移：plo con dis (xdis, ydis, zdis)plot contour smin range group 地层plot add contour smin range group 初衬上 上台阶 核心土 中台阶 下台阶 初衬仰拱 仰拱 any同时选择多个groupplot contour smax range group 初衬上 any group 初衬中 any group 初衬下 any group 初衬仰拱 anyplot contour smin range group 初衬上 any grou

54、p 初衬中 any group 初衬下 any group 初衬仰拱 anyplot contour smin range group 初衬上 not group 初衬中 not group 初衬下 not group 初衬仰拱 not 排除某个groupplot block state range group 地层 not开挖某个group中的一部分model null range group 1 y yy_1 yy_2*拱顶沉降set logfile 拱顶沉降节点位置.logset log onprint gp position range group 地层 x -0.1 0.1 y 0

55、50 z 12 13set log offset logfile 拱顶沉降.logset log onprint gp displacement range group 地层 x -0.1 0.1 y 0 50 z 12 13set log off*;上台阶水平收敛set logfile 上台阶水平收敛节点位置.logset log onprint gp position range group 地层 x -7 7 y 0 50 z 7.3 7.8set log offset logfile 上台阶水平收敛.logset log onprint gp displacement range gr

56、oup 地层 x -7 7 y 0 50 z 7.3 7.8set log off;下台阶水平收敛set logfile 下台阶水平收敛节点位置.logset log onprint gp position range group 地层 x -7.2 7.2 y 0 50 z 4.4 4.8set log offset logfile 下台阶水平收敛.logset log onprint gp displacement range group 地层 x -7.2 7.2 y 0 50 z 4.4 4.8set log off*;下台阶水平收敛set logfile 上台阶水平收敛节点位置.lo

57、gset log onprint gp position range group 地层 x -7 7 y 0 50 z 7.3 7.8set log offset logfile 上台阶水平收敛.logset log onprint gp displacement range group 地层 x -7 7 y 0 50 z 7.3 7.8set log off;下台阶水平收敛set logfile 下台阶水平收敛节点位置.logset log onprint gp position range group 地层 x -7.2 7.2 y 0 50 z 4.4 4.8set log offse

58、t logfile 下台阶水平收敛.logset log onprint gp displacement range group 地层 x -7.2 7.2 y 0 50 z 4.4 4.8set 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

59、 zhist gp zdis 0.034 0 60.044 ;地表history id 2 gp displacement id=4442hist 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

60、z 0.9 1.1初始化地应力,压应力为负ini sxx -3.5e6 syy -3.5e6 szz -3.5e6画等值线1、导出数据+res zizhong.sav;data_processset logfile 区域信息.logset log onprint zone information range group 初衬上 any group 初衬中 any group 初衬下 any group 初衬仰拱 anyset log offset logfile 挖前主应力.logset log onprint zone principal range set log offset logfi