flac常用命令

flac3d常用命令
1、最先需要掌握的命令有哪些？ 答：需要掌握 gen, ini, app, plo, solve 等建模、初始条件、边界条件、后处理 和求解的命令。 2、怎样输出模型的后处理图？ 答：File/Print type/Jpg file，然后选择 File/Print，将保存格式选择为 jpe 文件。 3、怎样调用一个文件？ 答：File/call 或者 call 命令 4、如何施加面力？ 答：app nstress 5、如何调整视图的大小、角度？ 答：综合使用 x, y, z, m, Shift 键，配合使用 Ctrl+R，Ctrl+Z 等快捷键。 6、如何进行边界约束？ 答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束） 。 7、如何知道每个单元的 ID？ 答：用鼠标双击单元的表面，可以知道单元的 ID 和坐标。 8、如何进行切片？ 答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示 z 方向应力的切片) 9、如何保存计算结果？ 答：save +文件名 10、如何调用已保存的结果？ 答：rest +文件名；或者 File / Restor 11、如何暂停计算？ 答：Esc 12、如何在程序中进行暂停，并可恢复计算？ 答：在命令中加入 pause 命令，用 continue 进行继续。 在我们分步求解中想得到某一个过程中的结果，不用等到全求完，还可以在 分布求解错误的时候就进行改正，而不是等到结果出来。 13、如何跳过某个计算步？ 答：在计算中按空格键跳过本次计算，自动进入下一步 14、Fish 是什么东西？Fish 是否一定要学？ 答：是 FLAC3D 的内置语言，可以用来进行参数化模型、完成命令本身不能 进行的功能。 Fish 可以不用学， 需要的时候查 Mannual 获得需要的变量就可以了。 15、FLAC3D 允许的命令文件格式有哪些？ 答：无所谓，只要是文本文件，什么后缀都可以。 16、如何调用一些可选模块？ 答：config dyn (fluid, creep, cppudm) 17、如何在圆柱体四周如何施加约束条件？ 可以用 fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not， 其中 r2
19、关于 dd、dip 的理解： dip-direction angle, measured in the global xy-plane clockwise from the positive y-axis 这个是从 y 轴到 x 顺时针方向所转过的角度； dip dip angle, measured in the negative z-direction from the global xy-plane 这个是 xy 平面与 z 轴负方向的夹角； dip 表示对称参照面与 xy 平面的夹角，对称参照面与 xy 平面的夹角在 xy 平面的投影是一条射线； dd 指的是从 y 轴正方向按顺时针（clockwise）方向转向所要确定面的法线 方向在 xy 平面上的投影的夹角， 然后再确定 dip， 指的是从 xy 平面转向所要 dip 确定的平面的角。 20、range 的命令： ①取圆球面的范围，就是同一个中心，不同半径之间夹的那一部分，用在给 球面施加力，约束时用该命令，命令： range annulus center xyz

（中心坐标）radius r1 r2（两半径） ，想要表达的是球 面时，r1r2 取值相差小点即可，如 radius 2.9 3.1 这时范围就是半径为 3 的面。 ②圆柱的环面施加力或约束时也用到该命令。 range cylinder end1 xyz end2 xyz radius r(end 是表示柱体断面的中心点坐标) range plane above 表示的是 z 大于 0 的上部结构 range plane blow 表示的是 z 小于 0 的下部结构 range dd a dip b origin xyz 表示的是走向角度 a，倾向角度 b 的面 dd 表示，当面与 xoy 面相交时的垂线与 y 轴正向的夹角，顺为正，说白话， 就是面的指向方向面是朝那个方向的意思。 dip 面与 xoy 面的夹角，与 z 轴负向为正，origin 表示面经过 xyz 点； range normal xyz 法线方向为 xyz 的面 range not 没试用过（范围外的东西，想想用到圆上取范围非常好） 。 21、如何调用 FISH 文件？ 新手调用 FISH 文件出错的原因：是：文件名为 fishcall.fis.txt(因为在记事本 里这个 txt 的后缀是默认的） 。修改方法是：将 fishcall.fis.txt 里面的内容全部复 制到新的记事本里，再按保存，保存时候将保存类型改为所有文件(*,*),这时你 便可在文件名栏内输入 fishcall.fis，类型也就正确了。 22、如何理解 Tension-p shear-p tension-n shear-n 表示的含义？ 在计算的循环里面,每个循环中， 每个 zone 都依据 failure criterion 处于不同的 状态，shear 和 tension 分别表示因受剪和受拉而处于塑性状态；n 表示 now， p 表示 previous，即分别表示在现在和以前的循环当中处于塑性状态。 23、绘制相关曲线命令的一些解释。 1) 在 plot hist m vs n 的形式里,m 代表 y 轴,n 代表 x 轴(不管 m,n 的正负)； 2) "-"表示对其值作"mirror"镜像。比如上面的 hist 1 gp ydisp 的值是 0~100,那 么 vs -1 就变成-100~0.以此类推； 3) 为什么要这么做呢?主要是为了符合岩土工程的习惯需要，例如桩载荷试验 曲线，沉降曲线等都不是画在第一象限内。 24、如何查看剖面上任一点的位移值？ print gp disp range x() y( ) z( )或 print gp disp range id plot gp disp range x() y( ) z( )或 plot gp disp range id
25、问 a 式与 b 式有何区别？ fix x y z range x -0.1 0.1 ..........................a fix x range x -0.1 0.1................................b a 表示 在 x=0 这个平面 是固定的； b 表示在 x=0 这个平面 只 x 方向不能发生运动。 26、initial 与 apply 有何区别？ initial 初始化命令，如初始化计算体的应力状态等； apply 边界条件限制命令，如施加边界的力、位移等约束等； initial 的应力状态会随计算过程的发生而发生改变， 一般体力需要初始化， 而 apply 施加的边界条件不会发生变化。 27、如何查看各个时段不平衡力的具体数值？ 采用 his 来记录计算，包括位移应力等 命令 his unbal； his

gp(zone) zdis range (0 0 0) 或者 id=？； 导出数据命令 his write n vs m begin 时步 end 时步+ file filename.his n 表示纪录的 id m 表示时步； 要导出不平衡力的具体数值 his unbal step 100000 or solve his write 1 vs step begin 1 end 1000 file 123.his 使用上述命令就可以查看各个时步下的不平衡力的具体数值。 28、如何显示变形轮廓线的命令？ plo ske magf 10 其中 10 为放大系数 29、怎样查看模型的矢量图？ pl grid vel (disp) red plo blo gro 可以看到不同的 group 的颜色分布（看模型的样子） ； plo gpfix red sk 看模型的边界情况 plo fap red sk 看模型的体力分布 plo con dis (xdis, ydis, zdis) 看模型的云图：位移 plo con sz (sy, sx, sxy, syz, sxz) 看模型的云图：应力 plo dis (xdis, ydis, zdis) 看模型的矢量图 plo info 看模型有多少单元、节点 plot grid 可以查看网格，plot grid num 可以查看节点号。 plot fix 查看边界条件 plot grid mark 查看预先留设的开挖边界 pl syy (sxx xdis ydis sig1..) fill；查看模型的云图。注意，带 fill 关键字就是查 看云图，不带 fill 关键字就是查看等直线图。 plo add ske 剖面的沉降云图上增加网格 plo add dis plane 剖面上的矢量 plo add axe 剖面上坐标值 30、怎样查看 INTERFACE？ 答：pl grid iface id red print iface id ; id 是 interface 的编号 31、FLAC 怎么调整视图？ 答： ctrl+z 或 shift+z ,然后点击要放大部位的左上角和右下角就可以放大某个
部位，如果要缩小，双击屏幕。 32、判断模型中网格点是一个还是两个的方法。 显示模型网格点的 ID，如果有两个 ID 号，说明是两个。节理面上的网格点 应该有两个 ID 号。 33、改变模型应力边界时 FLAC 与 3DEC 的不同之处。 在 3DEC 中，当改变应力边界时，模型表面的应力直接被更换，FLAC 中改 变应力边界时，模型表面的应力为两者应力相加。 34、flac 中影响运行时间的因素 （1）单元材料，结构面之间的刚度对比很大； （2）划分单元尺寸对比很大，尽量避免长细比大于 5：1 的细长单元。 35、FLAC 中有多少中材料模型。 （1）零模型 ：用于开挖回填的区域（2）各向同性弹性模型（3）正交各向异性 弹性模型（4）横观各向同性弹性模型（5）德鲁克－普拉格塑性模型（6）莫尔 －库仑塑性模型 ，一般岩土体的力学行为（7）节理化塑性模型 ：用于研究材料 破坏后行为（如边破稳定，屈服煤矿开采） （8）应变硬化/软化模型－库仑塑性 模型：薄板层状材料破坏后力学行为（9）双线性应变硬化/软化节理化塑性模型 临近层状土层开挖（10）双屈服塑性模型（11）修正的剑桥粘土模型。 36、对不平衡力的理解。 flac 中用最大不平衡力和典型内力大的比值 R 表示模型的不平衡力的相对大小， 该

值为百分数，且从不会减小到零，因此，模型不会达到绝对稳定的平横状态， 可以根据要求的精度设定 R 值，模型默认的 R 值是 1e-5。注意 R 值小仅表示所 有网格点的力都平衡， 但模型还可能发生塑性流动， 通过 plot vel 命令查看网格 点速度，一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。 37、简单命令解释： 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。 38、在固定面上施加位移时，fix 与 apply 不能同时使用，应该用 ini 和 fix, 才能施加上。 39、你可以用键盘上的上、下、左、右肩在窗口中移动模型。 40、在进行三轴压缩时候，如何施加围压？ 答：使用 apply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius 1，其
中 X 和 Z 轴是在我们眼睛正对的面。 42、在进行三轴压缩时候，如何对轴向力进行位移控制？ 答：使用 ini yvel 1e-7 range y -.1 .1 ini yvel 1e-7 range y 1.9 2.1 即可，数值可 以变化。 43、在进行边坡计算的时候，为什么用 sol, fos 命令运行？ 答：应该输入 sol fos slope3dfos.sav associated(不要也行） ；再用 restore 命令 调用即可。 44、在 FLAC3D 计算中，如何记录程序中的值，如果不想时时记录该怎么弄？ 答：用 hist wtite n1 vs n2 file.txt， 想间隔记录可用 set hist_rep n 45、想记录接触面上的相关力学参数该怎么办？ 答：hist inter i nstress /ndis/sstress/sdisp 46、接触面设置参数包括哪几个？ 答：coh dial fric kn ks ten. 47、如何输出单元格的数据？ 答；expgrid .flacsd 48、如何做切片？ 答：plot set plane ori x y z norm x y z plot con zd plane 这两行命令，用来显示指定面上的位移变化云图（把第二行改改，还可以显 示应力云图，等其他） ； plot con syy y 向应力云图 plot con szz outline on z 向应力云图加显示网格 plot con szz effective on z 向有效应力 plot con szz inter 1e4 云图显示增量 （后面的值越小， 图上颜色条带就越细， ） plot con szz min -10e4 显示最小值，图上小于这个数的值都用同一色表示 plot con szz ma

x 10e4 显示最大值，图上大于这个数的值都用同一色表示 plot fip outside on 显示力的箭头 现在看边界， 如果是不动面， fix x range…如果边界带速率， 用 或位移用 ini… set grav 0 0，-9.81 ini dens 1000 设置重力就一定要设置密度，不考虑重力就不要设计密度。两者共存的，设 计程序的时候一定要记的。 set mech force 50 solve 用 SET 命令设定了最大不平衡力的上限值为 50N，即：计算的最大不平衡 力一旦小于此值，我们便人为模型达到了平衡状态。求解过程终止。设置这个目 的是在一定精度要求下，加快计算（值越大，算的量越少，花的时间也越少） 。 49、FLAC 结果数据的导出 pl set back white 设置背景 set plot bitmap size (800,510) 设置图片大小 plot set caption size 37 设置图例数字大小 (default=35 范围 10-50) mainwin size 1 1 position 0 0.8 ; command window plot set window size 0.8 0.8 position 0.1 0.0 ; graphics window plot set center 100 0 40 根据情况更改。
50、一段 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 判断 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 执行什么定义的函数 51、有关命令的解释 1）apply nvel 0.1 plane dip 60 dd 270 range left_boun apply 施加边界条件，plane 法向速度为 0.1，plane 面的位置：倾角 60（在 xoy 平面从-Z 轴方向逆时针旋转） ，dd 走向 270 度，从+y 方向绕 z 轴顺时针旋转； 2）ini szz -5e4 grad 0 0 -1e4 ini 初始边界条件，szz=-5e4+x*0+y*0+-z*-1e4,注：-5e4 是 range 范围起点坐 标的 SZZ 值，x,y,z 所施加应力网格点的坐标。 52、FLAC 程式使用前准备步骤 步骤 1：依比例画出所欲分析的资料 在纸上画出地点的位置、地层资料、并简单标示距离及深度资料。 步骤 2：换算输入资料成同一单位 需谨慎检查输入资料的单位，如因单位不同而造成过大或过小的值，将会 造成 FLAC 无法计算，而产生 ERROR 信息。 步骤 3：应用公式简略计算 应用公式或依据经验，简略算出 FLAC 输出资料的范围，

以作为 Debug 及 输出资料分析时验证。 步骤 4：建立 x, y 坐标与 node i, j 之间的关系 在图上距离及深度的关系，建立 x, y 坐标系统，再由 x, y 坐标系统，转换与 网格关系， 为了便于以后输出资料的分析，故应确实掌握网格的位置及其相对应 的 x, y 坐标。 建议在敏感区域使用较密的网格，其它地方则使用较疏的网格，刚开始执行 程序时，不宜使用网格太大的网格数目，因尽量使网格总数少于 1000，以节省 时间。
53、FLAC 3D 程序的编写步骤： (1)建立 x, y 坐标与网格的关系，建议使用 Gen 指示： Gen x1,y1 x2,y2 ,x3,y3 x4,y4 i=i0,i1 j=j0,j1 (2)设定材料性质：prop (3)设定外力：Set Grav, Apply Pressure, ini sxx, Syy (4)设定边界条件：fix, free (5)求起始的应力平衡：solve (6)储存：Save 54、建模的过程中，或同一个文件中有重复应用的命令，就可以用宏命令了， 宏命令： macro 名称 ?命令的内容? 例：macro shiping ? plot con zdip outline on? 用 shiping 代替后面的引号里的东西。在建模中只要有相同的都可以用 shiping 代替，缩减了书写命令的时间，整体上也能让人看的更明白。注意，你 命名的符号最好跟 falc 自带的命令不要重复，处理方法是，名称后带数字。 如 shiping1；如果你想修改宏对象，改变命令的内容，在名称上加单引号就是替 换的意思。 例 macro ?shiping1? ?新的命令?。 55、很多时候我们要对单一的对象进行操作，例如进行模拟运算后，只想单独 的显示一部分的内容，单一的土层，锚杆，等就要用到 group。 group 给对象命名， 例： gen zone brick ------ size ---- group name name 随便取， 在 （ 不要跟 flac 的已有的命令重合就行） 这里的意思就是对这个体进行了命名，接下来想把这一 块模型挖掉加： model null range group name group 的用法：①group 直接加到模型的后面相当于给整个模型命名； ②group name range x xyy zz 给一个范围的体命名 。 56、上面说过建模一般只建一半然后用 reflect 命令（镜像） 。 命令：reflect normal xyz origin xyz （根据法线和过一点建的面）最常用至于 实际使用过程中有人问道：normal -1 0 0 与 normal 1 0 0 的区别，我试过两者的 效果是一样的，没区别，虽然方向不同，但表达的是同一个面。 reflect dd a dip b origin xyz 这个就是用里建立斜面镜像的。 有一点注意了，镜像的时候如果是模型镜像的话，相当于在建立一个模型， 故，在 reflect 的前面要加建模的关键字，gen zone reflect…… plot set rotation 20 0 30 ；设置观测角度的 plot set center auto ；设置观测中心点 plot set magnification 1.5 ；设置放大倍数的 plot add surface yellow plot add axes plot show ；显示 这一整部分就是我开头说的要自己整理出来

的句段了， 等以后用到的时候直 接插入， 作用就是显示模型的。 任何模型的显示都可以用他来表示。 很好很强大。 如果还要显示其他功能就把相应的命令在插入上面那部分就行了。 这个就是我常 用来查看模型的命令，大家也可以自己整理一份。 57、add 与 set 的区别 add 原来有的加新的用 add
set 原来没有的，设置用 set plot create GravV plot set plane dip=90 dd=0 origin=3,4,0 plot set rot 15 0 20 plot set center 2.5 4.2 4.0 plot add bound behind plot add bcont szz plane plot add axes plot show 。 我们想显示别的图像，又不想把这个弄好的图像丢掉，则用 plot create name 建立一个图像， 在建立新的图像时， plot create name1 （新的名） 显示以前的图像， 那么输入如下命令：plot current name。另一种方法，输入：plot show name 58、不同的区域定义不同的本构模型 model mohr range z=0,5 prop ClayeyGravel range z=0,5 在定义参数后面加个 range 范围，为了避免在输入命令时总要通过 range＋ 全局坐标域的方法来限定该命令执行的区域，我们可以先通过 range 命令来给不 同的区域命名，然后就可以用 range＋区域名的方法来限定范围了。这样就避免 了每次输入命令要限定在该区域执行时，都要输入一遍坐标范围的麻烦。具体如 何给区域命名，以及日后如何用区域名来限定该区域，可以看下面这个例子： range name trench x=-1,1 y=0,4 z=-2,2 model null range trench 59、右手法则：建立模型时，各点的排放时有规律的。 方法如下： 拿出右手伸直了 4 指合并指向就是 p0 到 p1 的方向，掌心指向就 是 p0 到 p2 的方向，拇指与 4 指垂直的方向就是 p0 到 p3 的方向。 很多人说网格后的云图有带齿什么的，改进的方法就是重复用下面的命令 Gen ADJust 用此命令来微调网格, 使之离散平滑。该命令可以连续使用以增进平滑效果的。 60、如果想看模型在计算过程中的变化情况 Set movie avi step 1 file 1.avi； 定义动画， 动画记录的频率为 1 时步一频率 （计 算量大的时候就没必要 1 时步的记录，可以 10,100 时步的记录，减少电脑的使 用） ： Movie start；动画开始 附件里这个命令“hist write 5 vs 2 file dibiao.txt” 能将结果按时步讲就算结果 输出 hist n= 5 ：N(Nstep)=5 指定了每迭代计算 5 次记录一次相关的值（如果不 指定，则系统默认值为 10） ； hist unbal 命令的第二行指定记录最大不平衡力（the maximum unbalanced force） ； hist gp zdisp 4,4,8 记录坐标值为（4，4，8）的节点 Z 方向的位移（这个用 的多，查看收敛量时就是这个命令，上面两个个人认为用处不大，也很少用过） 。