令把整体的网格画出
来，然后用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的值改成不同的数，就可以看到位移放大后的结构单元。
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 在进行三轴压缩时候，如何对轴向力进行位移控制？
答：使用 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
64 如何做切片？
答：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