1、8A版Flac3d50常用命令集锦建模1、调用文件：文件与工程在同一个文件夹，只写文件名即可：IfthecalledfileislocatedinthesamefolderastheFLAC3Dprojectfile,thenonlythefilenameneedbeenteredwiththeCALLcommand.不在同一个文件夹，全路径：Otherwise,thefilemaybecalledbyspecifyingitscompletepath(e.g.,c:myfolderfile.dat).Undo;撤销上一条命令2、创建旋转缩放视图3、建模命令modelmechmohr;莫尔库伦
2、模型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,-
3、10,0.p710,10,10plotcurrentplotPlot01plotclearplotzoneUndo;撤销命令setlogfile127G1001.tGtsetlogontruncatesetlogofflistzoneprinrangeG01y01z01;显示指定范围内各单元的主应力，结果如下Hist命令：命令编号按顺序从1开始：eachhistoryisnumberedsequentiallyfrom1asitisenteredviatheHISTORYcommand.查找显示所有的his命令：ReturntotheFlac3Dpromptandtypelisthistfor
4、alistingofthehistoriesandtheircorrespondingnumbers.histnstep5;每5步记录1次。默认是10步记录1次histaddratio;记录maGimumunbalancedforceratiohistaddgpzdisp448;记录节点的位移modelnullrangeG=2,4y=2,6z=5,10;ToeGcavatethetrench,enterThedimensionsoftheeGcavationregion(i.e.,thenullmaterial)aredefinedbytheG-,y-andz-coordinaterange.
5、setlarge;设置大变形large-strainlogicisspecifiediniGdis0ydis0zdis0;位移归零solvesavetrench;保持计算结果的文件“TRENCH.F3SAV”不需要后缀4、group命令(注意新版本的区别，groupzone)genzonebricksize666groupzoneTunnelrangeG16y06z15;设置groupmodelmechmohr;莫尔库伦模型modelmechnullrangegroupTunnel;groupzoneTunnelrangecylinderend1100105100end210095100rad
6、ius3groupzoneMatriGrangegroupTunnelnot;namingRANGESthatarefiGedinspaceandunrelated;tothemodelgeometry-representsstratigraphy.rangenameLayer1planedip0dd0ori000aboverangenameLayer2planedip0dd0ori000belowpropbulk1e8shear0.5e8coh0tens0fric35rangenrangeLayer1propbulk1e7shear0.3e7coh1e7tens0fric0rangenran
7、geLayer2plzonecolorbygroup4、Interface命令以下是切割模型法生成interface：寺河模拟用的plzonegenzonebrickp0000p120000p20400p3004size100202group1;genzonebrickp0004p120004p20404p3008size100202group2;genseparategroup1interface1wrapfirstgroup1secondgroup2interface1maGedge0.5显示interface：plotcreateplotPlot01plotaddsurface;显示表面
8、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)groupz
9、oneToprangegroupBasenotinterface1facerangeplanenorm(-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,0histaddgpyd
10、isp0,0,0;监测网格点（0，0，0）点的y方向上的位移histaddzonesyy0,1,0;监测单元（0，1，0）点的y方向上的应力histaddzonesyy1,1,0;监测单元（1，1，0）点的y方向上的应力histaddzoneszz0,0,-1Listhistorysetlargehistorynstep2historyaddfishtthistoryaddfishthetahistoryaddzonesGG(1,1,1)historyaddzonesyy(1,1,1)historyaddzoneszz(1,1,1)plotaddhis345vs1HISTORYaddfishv
11、ar;fishvar变量的监测5、solve命令?什么意思5、plot命令通过图形操作界面的方式：plzonecolorbygr;查看各个分组通过命令流的方式：plotcreateplotPlot01plotclearplotaddhistory1linestylecolorblackplotpostplotPlot01file127F1002.psmanualplotcreateplotPlot02plotclearplotaddhistory2linestylecolorblackplotpostplotPlot02file127F1003.psmanualplotcreateplotGr
12、avVplotcutaddplanedip90dd0origin3,4.1,0;截面plotsetcenter333plotseteye121812plotaddzoneplanebehindplaneononplaneoff;addlabelZonewhiteplotaddbcontszzrampgreyplotaddaGesplotpostplotGravVfile127F1004.psmanualreskw50plotcreateplotszz_contourplotcutaddplanedip90dd180origin0,18.9,0plotsetorient1202200plotse
13、tcenter45450plotaddzoneplanebehindplaneononplaneoff;plotaddbcontszzramprainbowplsetmag1.0plotaddaGesplotshowplotcreateGravVplotsetplanedip=90dd=0origin=3,4,0plotsetrot15020;thiswouldbeachievedinteractivelyplotsetcenter2.54.24.0;andisshowheresimplytoillustrateplotaddboundbehindplotaddbcontszzplaneplo
14、taddaGesplotshow这样我们就建立了一个叫“GravV”的图形，将它设为当前图层，对此图形设置一个倾角为90的面，（在G-y面上，约定z的负方向为下），倾向为0。（在G-y面上，从y的正方向开始沿顺时针方向计量），过G-y面上一点（G=3,y=4,z=0）,在平面后加一个网格，垂直方向的应力Qzz，最后给模型加轴线，整体图形轮廓与添加图形相反，显示每个区域中心的应力计算值，每个区域的颜色直接对应基本区域的应力，当键入show,我们建立的图形被显示PlotcreateDispcontPladdcontdispplanebehindPladdaGesPlshowplotsketch;绘
15、出轮廓线plotaddvel;添加矢量plotaddhis12vs3;绘出历史记录n，actLoad，p_solcyc15000plotCreateTrenchAddsurfaceyellowAddaGesblackShowplotcreateszz_contour;创建文件名为szz_contour（z方向的应力云图）plotsetcent445;视图的中心坐标为（4，4，5）plotsetrot20030;视图的旋转角度为（20，0，30）plotsetmag1.0;视图的放大倍数为1.0plotaddcontszzoutonshadeon;打开szz的等值线图plotaddaGes;打开
16、坐标轴线plotshow;显示窗口;不同的区域，材料赋值不同plzonecolorbymodelplzonecolorbygroup求解前对模型进行检查网格检查(如前所述)边界条件检查速度约束条件 plotgpfiG荷载条件 plotfap模型检查模型赋值 ploblockmodel参数赋值 ploblockpropGGG初始应力检查(如后所述)经常检查模型的响应plozconszz(syy,sGG) 应力场ploconzdis(ydis,Gdis) 位移场plo(add)zonecolorbystate 屈服状态plogpfiG 速度约束条件plofap 体力plohist(unbal)1
17、不平衡力（必需的）plsetmovieactiveonplsetmovieeGtensionjpgplsetmovieinterval20plsetmovieprefiGsihekuangkaiwaplsetmoviesizeiGiyplsetnamesihedacaigaoplbitmapview1filemyplot.pngplbitmapview1size1024768plcreateplcurrentplcutaddplaneddGGGdipGGGplcutaddplanenameGGGnormalGGGoriginGGGpldGfpleGcelplsvgplsetbackground
18、white;设置背景色plotaddhis345vs15、网格的连接、对称不同尺寸的zone，网格用attach连接，粘起来！6、plane命令6、range命令listrangerangenametrenchG=-1,1y=0,4z=-2,2;定义rangemodelmechnullrangenrangetrench;开挖eGcavatethetrench7、gp命令listgpposition;查看节点编号、位置坐标printgpdisprangeid105histaddgpzdis0,0,0histaddgpzvel0,0,0message(1)=string(gp_Gpos(pnt)m
19、essage(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(
20、INT)ZPNTzoneheadii=nzonetotalnumberofzones更多zone变量函数，详见fishreference9、config9、边界条件BoundaryandInitialConditionsBoundaryconditionsarenormallyspecifiedwiththeAPPLYandFIG/FREEcommands,andinitialconditionsareassignedwiththeINITIALcommandEGample：aconstant,compressiveGG-stresscomponentof10MPacanbeappliedto
21、aboundarylocatedatG=10withthecommandapplysGG-10e6rangeG10;G=10的面上10MPaapplysyy=-20e6grad0,0,20e5rangey-20z0,10在斜面上施加法向力10杂项（材料力学性质）11.M-C本构模型的选取在进行模拟时，采用摩尔-库仑（Mohr-Coulomb）强度准则来判断矿体与顶底板的破坏机理。该准则基本的内涵为：岩石的剪切破坏是由岩石破裂面产生的破坏剪应力受到岩石材料的内聚力和内摩擦力的抵抗，当破裂面上的破坏剪应力大于岩石的内聚力和内摩擦力时，即Fs0时，材料发生剪切破坏4：有限差分法(FDM)Fish函数
22、1、if条件从句2、loop条件从句newdefderives_mod=y_mod/(2.0G(1.0+p_ratio)b_mod=y_mod/(3.0G(1.0-2.0Gp_ratio)endsety_mod=5e8p_ratio=0.25derivelistb_mods_modgenzonebricksize(2,2,2)modelmechelasticpropbulk=b_modshear=s_modlistzonepropbulklistzonepropshearnewdefGGGsum=0prod=1loopn(1,10)sum=sum+nprod=prodGnend_loopend
23、GGGlistsumprodnewdefabcifGG0thenabc=33elseabc=11end_ifendsetGG=1listabcsetGG=-1listabcnewgenzonebricksize1035defplace_cablesloopn(1,5)z_d=float(n)-0.5commandselcablebeg(0.0,1.5,z_d)end(7.0,1.5,z_d)nseg7end_commandend_loopendplace_cablesplotaddzonetrans75plotaddselgeomnewgenzonebricksize1035modelmech
24、mohrpropbulk1e8shear0.3e8fric35propcoh1e3tens1e3inidens1000setgrav(0,0,-10);重力加速度fiGGyzrangez0.0fiGyrangey0.0fiGyrangey3.0fiGGrangeG0.0fiGGrangeG10.0setlargehistaddunbalsolvesavecab_striniGdis0ydis0zdis0histaddgpGdisp(0,1,5)defplace_cablesloopn(1,5)z_d=5.5-float(n)z_t=z_d+0.5z_b=z_d-0.5commandfreeGr
25、angeG0.0zz_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)nseg7selcablepropemod2e10ytension1e8Gcarea1.0.gr_k2e10gr_coh1e10gr_per1.0end_commandend_loopendplace_cablessavecab_endplotzonetrans75selcabblockforce
26、linewidth4newgenzonebricksize121modelmechmohrpropshear=1e8bulk=2e8.cohes=1e5tens=1e10fiGGyzrangey0.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)endget_addefloadload=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4)endhistaddfi
27、shloadhistaddgpydis(0,6,0)step1000plothis1vs-2newgenzonebricksize1035modelmechmohrpropbulk1e8shear0.3e8fric35propcoh1e3tens1e3inidens1000setgrav(0,0,-10);重力加速度fiGGyzrangez0.0fiGyrangey0.0fiGyrangey3.0fiGGrangeG0.0fiGGrangeG10.0setlarge;大变形histaddunbalsolvesavecab_str;保持计算结果文件，不用加后缀名.SaviniGdis0ydis0
28、zdis0;位移归零histaddgpGdisp(0,1,5)defplace_cablesloopn(1,5)z_d=5.5-float(n)z_t=z_d+0.5z_b=z_d-0.5commandfreeGrangeG0.0zz_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)nseg7selcablepropemod2e10ytension1e8Gcarea
29、1.0.gr_k2e10gr_coh1e10gr_per1.0end_commandend_loopendplace_cablessavecab_endplotzonetrans75selcabblockforcelinewidth4defparmglobalrad=4.0;radiusofsphericalcavitygloballen=10.0;lengthofouterboGedgeglobalin_size=6;numberofzonesalongoutercubeedgeglobalrad_size=10;numberofzonesinradialdirectionendparmgenzoneradbrickedgelensizein_sizein_sizein_sizerad_size&rat1.01.01.01.2dimradradraddefmountainlocalgp=gp_headloopwhilegp#nulllocalzz=sqrt(gp_Gpos(gp)2+gp_ypos(gp)2)
