建 模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 方方向向的的重重力力加加速速度度,,注注意意区区别别!!!!!!!!gen zone brick size 6,8,8 p0 -10, -10, -20 ... ;省略号表示写不下后面继续p1 10, -10, -20 ...p2 -10, 10, -20 ...p3 -10, -10, 0plot zonegen 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, 10plot current plot Plot01plot clearplot zoneUndo ;撤销命令set logfile 127X1001.txtset log on truncateset log offlist 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 typelist histfor a listing of the histories and their corresponding numbers.hist nstep 5 ;每5步记录1次。
默认是10步记录1次hist add ratio ; 记录maximum unbalanced force ratiohist add gp zdisp 4 4 8 ;记录节点的位移model null range x=2,4 y=2,6 z=5,10 ;To excavate the trench,enterThe 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 specifiedini xdis 0 ydis 0 zdis 0 ;位移归零solvesave trench ; 保持计算结果的文件“TRENCH.F3SA V”不需要后缀4、group命令(注意新版本的区别,group zone)gen zone brick size 6 6 6group zone Tunnel range x 1 6 y 0 6 z 1 5 ;设置groupmodel 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 aboverange name Layer2 plane dip 0 dd 0 ori 0 0 0 belowprop bulk 1e8 shear 0.5e8 coh 0 tens 0 fric 35 range nrange Layer1prop bulk 1e7 shear 0.3e7 coh 1e7 tens 0 fric 0 range nrange Layer2pl zone colorby group4、Interface命令以下是切割模型法生成interface:寺河模拟用的pl zonegen 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 1interface 1 wrap first group 1 second group 2interface 1 maxedge 0.5显示interface:plot create plot Plot01plot add surface;显示表面plot add interface red;显示交互面为红色plot show用导来导去法,已经在5.0上运行过,很好,值得借鉴pl zone colorby groupgen zone brick size 3 3 3 &p0 (0,0,0) p1 (3,0,0) p2 (0,3,0) p3 (0,0,1.5) &p4 (3,3,0) p5 (0,3,1.5) p6 (3,0,4.5) p7 (3,3,4.5)group zone Basegen zone brick size 3 3 3 &p0 (0,0,2.5) p1 (3,0,5.5) p2 (0,3,2.5) p3 (0,0,7) &p4 (3,3,5.5) p5 (0,3,7) p6 (3,0,7) p7 (3,3,7)group zone Top range group Base notinterface 1 face range plane norm (-1,0,1) origin (1.5,1.5,3) dist 0.1 pl add interfacepauseini zpos add -1.0 range group Top4、hist命令(5.0版本需要加add)hist nstep=5 //设置采用间隔为5 等效于set hist_rep=5 hist add id =2 zone szz 0,0,-1hist unbal //最大非平衡力hist add gp zdisp 4 4 8print histhist add gp zdis 0,0,0hist add gp zvel 0,0,0hist add gp ydisp 0,0,0;监测网格点(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,-1List historyset largehistory nstep 2history add fish tthistory add fish thetahistory 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 1HI STORY add fish var ;fish var变量的监测5、solve命令什么意思5、plot命令通过图形操作界面的方式:pl zone colorby gr ;查看各个分组通过命令流的方式:plot create plot Plot01plot clearplot add history 1 linestyle color blackplot post plot Plot01 file 127F1002.ps manualplot create plot Plot02plot clearplot add history 2 linestyle color blackplot post plot Plot02 file 127F1003.ps manualplot create plot GravVplot cut add plane dip 90 dd 0 origin 3,4.1,0 ;截面plot set center 3 3 3plot set eye 12 18 12plot add zone plane behindplane on onplane off ;addlabel "Zone" white plot add bcont szz ramp greyplot add axesplot post plot GravV file 127F1004.ps manualres kw50plot create plot szz_contourplot cut add plane dip 90 dd 180 origin 0,18.9,0plot set orient 120 220 0plot set center 45 45 0plot add zone plane behindplane on onplane off ;plot add bcont szz ramp rainbowpl set mag 1.0plot add axesplot showplot create GravVplot set plane dip=90 dd=0 origin=3,4,0plot set rot 15 0 20 ;this would be achieved interactivelyplot set center 2.5 4.2 4.0 ;and is show here simply to illustrate plot add bound behindplot add bcont szz planeplot add axesplot show这样我们就建立了一个叫“Grav V”的图形,将它设为当前图层,对此图形设置一个倾角为90°的面,(在x-y面上,约定z的负方向为下),倾向为0°。