flac3d5.0常用命令集锦

...wd......wd......wd...建模1、调用文件：①文件与工程在同一个文件夹，只写文件名即可：IfthecalledfileislocatedinthesamefolderastheFLAC3Dprojectfile,thenonlythefilenameneed

beenteredwiththeCALLcommand.②不在同一个文件夹，全路径：Otherwise,thefilemaybecalledbyspecifyingitscompletepath(e.g.,c:\myfolder\file.dat).Undo;撤销上一条命令2、创立旋转缩放视图3、建模命令modelmechmohr;莫尔库伦模型modelmechelastic;弹性模型setgrav0,0,-9.81;重力加速度negativez-direction.〔垂直向下！常用的〕下面这代码，是沿着-y方向的重力加速度，注意区别！！！！genzonebricksize6,8,8p0-10,-10,-20...;省略号表示写不下后面继续p110,-10,-20...

p2-10,10,-20...

p3-10,-10,0plotzonegenzonebricksize6,8,8p0-10,-10,-20...;不规则六面体p110,-10,-20p2-10,10,-20...p3-10,-10,0p410,10,-20...p5-10,10,10p610,-10,0...p710,10,10plotcurrentplotPlot01plotclearplotzoneUndo;撤销命令setlogfile127X1001.txtsetlogontruncatesetlogofflistzoneprinrangex01y01z01;显示指定范围内各单元的主应力，结果如下Hist命令：①命令编号按顺序从1开场：eachhistoryisnumberedsequentiallyfrom1asitisenteredviatheHISTORYcommand.

②查找显示所有的his命令：ReturntotheFlac3D>promptandtype

listhist

foralistingofthehistoriesandtheircorrespondingnumbers.histnstep5;每5步记录1次。默认是10步记录1次histaddratio;记录maximumunbalancedforceratiohistaddgpzdisp448;记录节点的位移modelnullrangex=2,4y=2,6z=5,10;Toexcavatethetrench,enter

Thedimensionsoftheexcavationregion(i.e.,thenullmaterial)aredefinedbythex-,y-andz-coordinaterange.setlarge;设置大变形large-strainlogicisspecifiedinixdis0ydis0zdis0;位移归零solvesavetrench;保持计算结果的文件“TRENCH.F3SAV〞不需要后缀4、group命令(注意新版本的区别，groupzone)genzonebricksize666

groupzoneTunnelrangex16y06z15;设置group

modelmechmohr;莫尔库伦模型

modelmechnullrangegroupTunnel;groupzoneTunnelrangecylinderend1100105100end210095100radius3

groupzoneMatrixrangegroupTunnelnot

;namingRANGESthatarefixedinspaceandunrelated

;tothemodelgeometry-representsstratigraphy...

rangenameLayer1planedip0dd0ori000above

rangenameLayer2planedip0dd0ori000belowpropbulk1e8shear0.5e8coh0tens0fric35rangenrangeLayer1

propbulk1e7shear0.3e7coh1e7tens0fric0rangenrangeLayer2plzonecolorbygroup4、Interface命令以下是切割模型法生成interface：寺河模拟用的plzonegenzonebrickp0000p120000p20400p3004size100202group1;genzonebrickp0004p120004p20404p3008size100202group2;genseparategroup1interface1wrapfirstgroup1secondgroup2interface1maxedge0.5显示interface：plotcreateplotPlot01plotaddsurface;显示外表plotaddinterfacered;显示交互面为红色plotshow用导来导去法，已经在5.0上运行过，很好，值得借鉴plzonecolorbygroupgenzonebricksize333&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)groupzoneBasegenzonebricksize333&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)groupzoneToprangegroupBasenotinterface1facerangeplanenorm(-1,0,1)origin(1.5,1.5,3)dist0.1pladdinterfacepauseinizposadd-1.0rangegroupTop4、hist命令(5.0版本需要加add！！！)histnstep=5//设置采用间隔为5等效于sethist_rep=5histaddid=2zoneszz0,0,-1histunbal//最大非平衡力histaddgpzdisp448printhisthistaddgpzdis0,0,0histaddgpzvel0,0,0histaddgpydisp0,0,0;监测网格点〔0，0，0〕点的y方向上的位移histaddzonesyy0,1,0;监测单元〔0，1，0〕点的y方向上的应力histaddzonesyy1,1,0;监测单元〔1，1，0〕点的y方向上的应力histaddzoneszz0,0,-1Listhistorysetlarge

historynstep2

historyaddfishtt

historyaddfishtheta

historyaddzonesxx(1,1,1)

historyaddzonesyy(1,1,1)

historyaddzoneszz(1,1,1)plotaddhis345vs1HISTORYaddfishvar;fishvar变量的监测5、solve命令?????什么意思5、plot命令通过图形操作界面的方式：plzonecolorbygr;查看各个分组通过命令流的方式：plotcreateplotPlot01plotclearplotaddhistory1linestylecolorblackplotpostplotPlot01file127F1002.psmanualplotcreateplotPlot02plotclearplotaddhistory2linestylecolorblackplotpostplotPlot02file127F1003.psmanualplotcreateplotGravVplotcutaddplanedip90dd0origin3,4.1,0;截面plotsetcenter333plotseteye121812plotaddzoneplanebehindplaneononplaneoff;addlabel"Zone"whiteplotaddbcontszzrampgreyplotaddaxesplotpostplotGravVfile127F1004.psmanualreskw50plotcreateplotszz_contourplotcutaddplanedip90dd180origin0,18.9,0plotsetorient1202200plotsetcenter45450plotaddzoneplanebehindplaneononplaneoff;plotaddbcontszzramprainbowplsetmag1.0plotaddaxesplotshowplotcreateGravVplotsetplanedip=90dd=0origin=3,4,0plotsetrot15020;thiswouldbeachievedinteractivelyplotsetcenter2.54.24.0;andisshowheresimplytoillustrateplotaddboundbehindplotaddbcontszzplaneplotaddaxesplotshow这样我们就建设了一个叫“GravV〞的图形，将它设为当前图层，对此图形设置一个倾角为90°的面，〔在x-y面上，约定z的负方向为下〕，倾向为0°。〔在x-y面上，从y的正方向开场沿顺时针方向计量〕，过x-y面上一点〔x=3,y=4,z=0〕,在平面后加一个网格，垂直方向的应力Qzz，最后给模型加轴线，整体图形轮廓与添加图形相反，显示每个区域中心的应力计算值，每个区域的颜色直接对应根本区域的应力，当键入show,我们建设的图形被显示PlotcreateDispcontPladdcontdispplanebehindPladdaxesPlshowplotsketch;绘出轮廓线plotaddvel;添加矢量plotaddhis12vs3;绘出历史记录n，actLoad，p_solcyc15000plotCreateTrenchAddsurfaceyellowAddaxesblackShowplotcreateszz_contour;创立文件名为szz_contour〔z方向的应力云图〕plotsetcent445;视图的中心坐标为〔4，4，5〕plotsetrot20030;视图的旋转角度为〔20，0，30〕plotsetmag1.0;视图的放大倍数为1.0plotaddcontszzoutonshadeon;翻开szz的等值线图plotaddaxes;翻开坐标轴线plotshow;显示窗口;不同的区域，材料赋值不同plzonecolorbymodelplzonecolorbygroup求解前对模型进展检查网格检查(如前所述)边界条件检查速度约束条件 plotgpfix荷载条件 plotfap模型检查模型赋值 ploblockmodel参数赋值 ploblockprop***初始应力检查(如后所述)经常检查模型的响应plozconszz(syy,sxx) 应力场ploconzdis(ydis,xdis) 位移场plo(add)zonecolorbystate 屈服状态plogpfix 速度约束条件plofap 体力plohist(unbal)1 不平衡力〔必需的〕plsetmovieactiveonplsetmovieextensionjpgplsetmovieinterval20plsetmovieprefixsihekuangkaiwaplsetmoviesizeixiyplsetnamesihedacaigaoplbitmapview1filemyplot.pngplbitmapview1size1024768plcreate<viewid>plcurrent<viewid>plcutaddplaneddxxxdipxxxplcutaddplanenamexxxnormalxxxoriginxxxpldxf<viewid><filefilename.dxf>plexcel<filefilename.xlsx>plsvg<viewid><filefilename.svg>plsetbackgroundwhite;设置背风光plotaddhis345vs15、网格的连接、对称……不同尺寸的zone，网格用attach连接，粘起来！！！！！！6、plane命令6、range命令listrangerangenametrenchx=-1,1y=0,4z=-2,2;定义rangemodelmechnullrangenrangetrench;;开挖excavatethetrench7、gp命令listgpposition;查看节点编号、位置坐标printgpdisprangeid105histaddgpzdis0,0,0histaddgpzvel0,0,0message(1)=string(gp_xpos(pnt))message(1)=message(1)+','+string(gp_zpos(pnt))message(1)=message(1)+','+string(abs(gp_zdisp(pnt)))GPPNTfindgp(INT)GPPNTgpfind(INT)GPPNTgpheadii=ngptotalnumberofgridpoints更多节点变量函数，详见fishreference8、zone命令listzoneprintzonestressrangeid2381anyid2361any2341any;输出网格单元的应力histzoneszz0,0,-1ZPNTfindzone(INT)INTzfind(INT)ZPNTzoneheadii=nzonetotalnumberofzones更多zone变量函数，详见fishreference9、config9、边界条件BoundaryandInitialConditionsBoundaryconditionsarenormallyspecifiedwiththeAPPLYandFIX/FREEcommands,andinitialconditionsareassignedwiththeINITIALcommandExample：aconstant,compressivexx-stresscomponentof10MPacanbeappliedtoaboundarylocatedatx=10withthecommand

applysxx-10e6rangex10;x=10的面上10MPaapplysyy=-20e6grad0,0,20e5rangey-20z0,10在斜面上施加法向力10杂项〔材料力学性质〕11.M-C本构模型的选取在进展模拟时，采用摩尔-库仑〔Mohr-Coulomb〕强度准则来判断矿体与顶底板的破坏机理。该准则根本的内涵为：岩石的剪切破坏是由岩石破裂面产生的破坏剪应力受到岩石材料的内聚力和内摩擦力的抵抗，当破裂面上的破坏剪应力大于岩石的内聚力和内摩擦力时，即Fs>0时，材料发生剪切破坏[4]：有限差分法(FDM)Fish函数1、if条件从句2、loop条件从句newdefderives_mod=y_mod/(2.0*(1.0+p_ratio))b_mod=y_mod/(3.0*(1.0-2.0*p_ratio))endset@y_mod=5e8@p_ratio=0.25@derivelist@b_mod@s_modgenzonebricksize(2,2,2)

modelmechelastic

propbulk=@b_modshear=@s_mod

listzonepropbulk

listzonepropshearnew

defxxx

sum=0

prod=1

loopn(1,10)

sum=sum+n

prod=prod*n

end_loop

end@xxx

list@sum@prodnewdefabcifxx>0thenabc=33elseabc=11end_ifendset@xx=1list@abcset@xx=-1list@abcnewgenzonebricksize1035defplace_cablesloopn(1,5)z_d=float(n)-0.5commandselcablebeg(0.0,1.5,@z_d)end(7.0,1.5,@z_d)nseg7end_commandend_loopend@place_cablesplotaddzonetrans75plotaddselgeomnewgenzonebricksize1035modelmechmohrpropbulk1e8shear0.3e8fric35propcoh1e3tens1e3inidens1000setgrav(0,0,-10);;重力加速度fixxyzrangez0.0fixyrangey0.0fixyrangey3.0fixxrangex0.0fixxrangex10.0setlargehistaddunbalsolvesavecab_strinixdis0ydis0zdis0histaddgpxdisp(0,1,5)defplace_cablesloopn(1,5)z_d=5.5-float(n)z_t=z_d+0.5z_b=z_d-0.5commandfreexrangex0.0z@z_b,@z_tsolveselcablebeg(0.0,0.5,@z_d)end(7.0,0.5,@z_d)nseg7selcablebeg(0.0,1.5,@z_d)end(7.0,1.5,@z_d)nseg7selcablebeg(0.0,2.5,@z_d)end(7.0,2.5,@z_d)nseg7selcablepropemod2e10ytension1e8xcarea1.0...gr_k2e10gr_coh1e10gr_per1.0end_commandend_loopend@place_cablessavecab_endplotzonetrans75selcabblockforcelinewidth4newgenzonebricksize121modelmechmohrpropshear=1e8bulk=2e8...cohes=1e5tens=1e10fixxyzrangey0.0applyyvel-1e-5rangey2.0defget_adad1=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_addefloadload=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4)endhistaddfishloadhistaddgpydis(0,6,0)step1000plothis1vs-2newgenzonebricksize1035modelmechmohrpropbulk1e8shear0.3e8fric35propcoh1e3tens1e3inidens1000setgrav(0,0,-10);;重力加速度fixxyzrangez0.0fixyrangey0.0fixyrangey3.0fixxrangex0.0fixxrangex10.0setlarge;;大变形histaddunbalsolvesavecab_str;;保持计算结果文件，不用加后缀名.Savinixdis0ydis0zdis0;;位移归零histaddgpxdisp(0,1,5)defplace_cablesloopn(1,5)z_d=5.5-float(n)z_t=z_d+0.5z_b=z_d-0.5commandfreexrangex0.0z@z_b,@z_tsolveselcablebeg(0.0,0.5,@z_d)end(7.0,0.5,@z_d)nseg7selcablebeg(0.0,1.5,@z_d)end(7.0,1.5,@z_d)nseg7selcablebeg(0.0,2.5,@z_d)end(7.0,2.5,@z_d)nseg7selcablepropemod2e10ytension1e8xcarea1.0...gr_k2e10gr_coh1e10gr_per1.0end_commandend_loopend@place_cablessavecab_endplotzonetrans75selcabblockforcelinewidth4defparm

globalrad=4.0;radiusofsphericalcavity

globallen=10.0;lengthofouterboxedg