flac3d5.0常用命令集锦汇编

建 模 1 调用文件 文件与工程在同一个文件夹 只写文件名即可 If the called file is located in the same folder as the FLAC3D project file then only the file name need be entered with the CALL command 不在同一个文件夹 全路径 Otherwise the file may be called by specifying its complete path e g c myfolder file dat Undo 撤销上一条命令 2 创建 旋转 缩放视图 3 建模命令 model mech mohr 莫尔库伦模型 model mech elastic 弹性模型 set grav 0 0 9 81 重力加速度 negative z direction 垂直向下 常用的 下面这代码 是沿着 下面这代码 是沿着下面这代码 是沿着 y y y 方向的重力加速度 注意区别 方向的重力加速度 注意区别 方向的重力加速度 注意区别 gen zone brick size 6 8 8 p0 10 10 20 省略号表示写不下后面继续 p1 10 10 20 p2 10 10 20 p3 10 10 0 plot zone gen zone brick size 6 8 8 p0 10 10 20 不规则六面体 p1 10 10 20 p2 10 10 20 p3 10 10 0 p4 10 10 20 p5 10 10 10 p6 10 10 0 p7 10 10 10 plot current plot Plot01 plot clear plot zone Undo 撤销命令 set logfile 127X1001 txt set log on truncate set log off list zone prin range x 0 1 y 0 1 z 0 1 显示指定范围内各单元的主应力 结果如下 Hist 命令 命令编号按顺序从 1 开始 each history is numbered sequentially from 1 as it is entered via the HISTORY command 查找显示所有的 his 命令 Return to the Flac3D prompt and type list hist for a listing of the histories and their corresponding numbers hist nstep 5 每 5 步记录 1 次 默认是 10 步记录 1 次 hist add ratio 记录 maximum unbalanced force ratio hist add gp zdisp 4 4 8 记录节点的位移 model null range x 2 4 y 2 6 z 5 10 To excavate the trench enter The dimensions of the excavation region i e the null material are defined by the x y and z coordinate range set large 设置大变形 large strain logic is specified ini xdis 0 ydis 0 zdis 0 位移归零 solve save trench 保持计算结果的文件 TRENCH F3SAV 不需要后缀 4 group 命令 注意新版本的区别 group zone gen zone brick size 6 6 6 group zone Tunnel range x 1 6 y 0 6 z 1 5 设置 group model mech mohr 莫尔库伦模型 model mech null range group Tunnel group zone Tunnel range cylinder end1 100 105 100 end2 100 95 100 radius 3 group zone Matrix range group Tunnel not naming RANGES that are fixed in space and unrelated to the model geometry represents stratigraphy range name Layer1 plane dip 0 dd 0 ori 0 0 0 above range name Layer2 plane dip 0 dd 0 ori 0 0 0 below prop bulk 1e8 shear 0 5e8 coh 0 tens 0 fric 35 range nrange Layer1 prop bulk 1e7 shear 0 3e7 coh 1e7 tens 0 fric 0 range nrange Layer2 pl zone colorby group 4 Interface 命令 以下是切割模型法生成 interface 寺河模拟用的 pl zone gen zone brick p0 0 0 0 p1 200 0 0 p2 0 40 0 p3 0 0 4 size 100 20 2 group 1 gen zone brick p0 0 0 4 p1 200 0 4 p2 0 40 4 p3 0 0 8 size 100 20 2 group 2 gen separate group 1 interface 1 wrap first group 1 second group 2 interface 1 maxedge 0 5 显示 interface plot create plot Plot01 plot add surface 显示表面 plot add interface red 显示交互面为红色 plot show 用导来导去法 已经在用导来导去法 已经在 5 0 上运行过 很好 值得借鉴上运行过 很好 值得借鉴 pl zone colorby group gen zone brick size 3 3 3 监测网格点 0 0 0 点的 y 方向上的位移 hist add zone syy 0 1 0 监测单元 0 1 0 点的 y 方向上的应力 hist add zone syy 1 1 0 监测单元 1 1 0 点的 y 方向上的应力 hist add zone szz 0 0 1 List history set large history nstep 2 history add fish tt history add fish theta history add zone sxx 1 1 1 history add zone syy 1 1 1 history add zone szz 1 1 1 plot add his 3 4 5 vs 1 HISTORY add fish var fish var 变量的监测变量的监测 5 solve 命令 什么意思 5 plot 命令 通过图形操作界面的方式 pl zone colorby gr 查看各个分组 通过命令流的方式 plot create plot Plot01 plot clear plot add history 1 linestyle color black plot post plot Plot01 file 127F1002 ps manual plot create plot Plot02 plot clear plot add history 2 linestyle color black plot post plot Plot02 file 127F1003 ps manual plot create plot GravV plot cut add plane dip 90 dd 0 origin 3 4 1 0 截面 plot set center 3 3 3 plot set eye 12 18 12 plot add zone plane behindplane on onplane off addlabel Zone white plot add bcont szz ramp grey plot add axes plot post plot GravV file 127F1004 ps manual res kw50 plot create plot szz contour plot cut add plane dip 90 dd 180 origin 0 18 9 0 plot set orient 120 220 0 plot set center 45 45 0 plot add zone plane behindplane on onplane off plot add bcont szz ramp rainbow pl set mag 1 0 plot add axes plot show plot create GravV plot set plane dip 90 dd 0 origin 3 4 0 plot set rot 15 0 20 this would be achieved interactively plot set center 2 5 4 2 4 0 and is show here simply to illustrate plot add bound behind plot add bcont szz plane plot add axes plot show 这样我们就建立了一个叫 Grav V 的图形 将它设为当前图层 对此图形设 置一个倾角为 90 的面 在 x y 面上 约定 z 的负方向为下 倾向为 0 在 x y 面上 从 y 的正方向开始沿顺时针方向计量 过 x y 面上一点 x 3 y 4 z 0 在平面后加一个网格 垂直方向的应力 Qzz 最后给模型加轴 线 整体图形轮廓与添加图形相反 显示每个区域中心的应力计算值 每个区 域的颜色直接对应基本区域的应力 当键入 show 我们建立的图形被显示 Plot create Dispcont Pl add cont disp plane behind Pl add axes Pl show plot sketch 绘出轮廓线 plot add vel 添加矢量 plot add his 1 2 vs 3 绘出历史记录 n actLoad p sol cyc 15000 plot Create Trench Add surface yellow Add axes black Show plot create szz contour 创建文件名为 szz contour z 方向的应力云图 plot set cent 4 4 5 视图的中心坐标为 4 4 5 plot set rot 20 0 30 视图的旋转角度为 20 0 30 plot set mag 1 0 视图的放大倍数为 1 0 plot add cont szz out on shade on 打开 szz 的等值线图 plot add axes 打开坐标轴线 plot show 显示窗口 不同的区域 材料赋值不同 pl zone colorby model pl zone colorby group 求解前对模型进行检查 网格检查 如前所述 边界条件检查 速度约束条件 plot gpfix 荷载条件 plot fap 模型检查 模型赋值 plo block model 参数赋值 plo block prop 初始应力检查 如后所述 经常检查模型的响应 plo zcon szz syy sxx 应力场 plo con zdis ydis xdis 位移场 plo add zone colorby state屈服状态 plo gpfix 速度约束条件 plo fap 体力 plo hist unbal 1不平衡力 必需的 pl set movieactive on pl set movieextension jpg pl set movieinterval 20 pl set movieprefix sihekuangkaiwa pl set moviesize ix iy pl set name sihedacaigao pl bitmap view 1 file myplot png pl bitmap view 1 size 1024 768 pl create pl current pl cut add plane dd xxx dip xxx pl cut add plane name xxx normal xxx origin xxx pl dxf pl excel pl svg pl set background white 设置背景色 plot add his 3 4 5 vs 1 5 网格的连接 对称 不同尺寸的 zone 网格用 attach 连接 粘起来 6 plane 命令 6 range 命令 list range range name trench x 1 1 y 0 4 z 2 2 定义 range model mech null range nrange trench 开挖 excavate the trench 7 gp 命令 list gp position 查看节点编号 位置坐标 print gp disp range id 105 hist add gp zdis 0 0 0 hist add gp zvel 0 0 0 message 1 string gp xpos pnt message 1 message 1 string gp zpos pnt message 1 message 1 string abs gp zdisp pnt GP PNT find gp INT GP PNT gp find INT GP PNT gp head ii ngp total number of gridpoints 更多节点变量函数 详见 fish reference 8 zone 命令 list zone print zone stress range id 2381 any id 2361 any 2341 any 输出网格单元的应力 hist zone szz 0 0 1 Z PNT find zone INT INT z find INT Z PNT zone head ii nzone total number of zones 更多 zone 变量函数 详见 fish reference 9 config 9 边界条件 Boundary and Initial Conditions Boundary conditions are normally specified with the APPLY and FIX FREE commands and initial conditions are assigned with the INITIAL command Example a constant compressive xx stress component of 10 MPa can be applied to a boundary located at x 10 with the command apply sxx 10e6 range x 10 x 10 的面上 10MPa apply syy 20e6 grad 0 0 20e5 range y 20 z 0 10 在斜面上施加法向力 10 杂项 材料力学性质 11 M C 本构模型的选取 在进行模拟时 采用摩尔 库仑 Mohr Coulomb 强度准则来判断矿体与 顶底板的破坏机理 该准则基本的内涵为 岩石的剪切破坏是由岩石破裂面产 生的破坏剪应力受到岩石材料的内聚力和内摩擦力的抵抗 当破裂面上的破坏 剪应力大于岩石的内聚力和内摩擦力时 即 Fs 0 时 材料发生剪切破坏 4 有限差分法 FDM Fish 函数 1 if 条件从句 2 loop 条件从句 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 list b mod s mod gen zone brick size 2 2 2 model mech elastic prop bulk b mod shear s mod list zone prop bulk list zone prop shear new def xxx sum 0 prod 1 loop n 1 10 sum sum n prod prod n end loop end xxx list sum prod new def abc if xx 0 then abc 33 else abc 11 end if end set xx 1 list abc set xx 1 list abc new gen zone brick size 10 3 5 def place cables loop n 1 5 z d float n 0 5 command sel cable beg 0 0 1 5 z d end 7 0 1 5 z d nseg 7 end command end loop end place cables plot add zone trans 75 plot add sel geom new gen zone brick size 10 3 5 model mech mohr prop bulk 1e8 shear 0 3e8 fric 35 prop coh 1e3 tens 1e3 ini dens 1000 set grav 0 0 10 重力加速度 fix x y z range z 0 0 fix y range y 0 0 fix y range y 3 0 fix x range x 0 0 fix x range x 10 0 set large hist add unbal solve save cab str ini xdis 0 ydis 0 zdis 0 hist add 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 0 0 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 plot zone trans 75 sel cabblock force line width 4 new gen zone brick size 1 2 1 model mech mohr prop shear 1e8 bulk 2e8 cohes 1e5 tens 1e10 fix x y z range y 0 0 apply yvel 1e 5 range y 2 0 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 add fish load hist add gp ydis 0 6 0 step 1000 plot his 1 vs 2 new gen zone brick size 10 3 5 model mech mohr prop bulk 1e8 shear 0 3e8 fric 35 prop coh 1e3 tens 1e3 ini dens 1000 set grav 0 0 10 重力加速度 fix x y z range z 0 0 fix y range y 0 0 fix y range y 3 0 fix x range x 0 0 fix x range x 10 0 set large 大变形 hist add unbal solve save cab str 保持计算结果文件 不用加后缀名 Sav ini xdis 0 ydis 0 zdis 0 位移归零 hist add 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 0 0 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 plot zone trans 75 sel cabblock