1、建 模1、调用文件：文件与工程在同一个文件夹，只写文件名即可：If the called file is located in the same folder as the FLAC3D project file, then only the file name needbe entered with the CALL command. 不在同一个文件夹，全路径：Otherwise, the file may be called by specifying its complete path (e.g., c:myfolderfile.dat).Undo ;撤销上一条命令2、创建 旋转 缩放视
2、图3、建模命令model mech mohr ;莫尔库伦模型model mech elastic ;弹性模型set grav 0, 0, -9.81 ; 重力加速度 negative z-direction. （垂直向下！常用的）下面这代码，是沿着下面这代码，是沿着下面这代码，是沿着 -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
3、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 clear plot zoneUndo ;撤销命令set logfile 127X1001.txtset log on truncateset log offlist zone prin range x 0 1 y
4、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 步记录
5、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 logi
6、c is specifiedini xdis 0 ydis 0 zdis 0 ;位移归零solvesave trench ; 保持计算结果的文件 “TRENCH.F3SAV” 不需要后缀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
7、100 105 100 end2 100 95 100 radius 3group 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
8、.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 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
9、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 监测网
10、格点（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)plo
11、t add his 3 4 5 vs 1HISTORY 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 l
12、inestyle 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“ whiteplot add bcont szz ramp greyplot add axesplot post pl
13、ot 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 GravVplo
14、t 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 illustrateplot 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 DispcontPl add cont disp plane behindPl add axesPl show plot sketch;绘出轮廓线plot add vel;添加矢量plot add his 1 2 vs 3
