flac3d5.0常用命令集錦

...wd......wd......wd...建模1、調(diào)用文件：①文件與工程在同一個(gè)文件夾，只寫(xiě)文件名即可：IfthecalledfileislocatedinthesamefolderastheFLAC3Dprojectfile,thenonlythefilenameneed

beenteredwiththeCALLcommand.②不在同一個(gè)文件夾，全路徑：Otherwise,thefilemaybecalledbyspecifyingitscompletepath(e.g.,c:\myfolder\file.dat).Undo;撤銷(xiāo)上一條命令2、創(chuàng)立旋轉(zhuǎn)縮放視圖3、建模命令modelmechmohr;莫爾庫(kù)倫模型modelmechelastic;彈性模型setgrav0,0,-9.81;重力加速度negativez-direction.〔垂直向下！常用的〕下面這代碼，是沿著-y方向的重力加速度，注意區(qū)別?。。。enzonebricksize6,8,8p0-10,-10,-20...;省略號(hào)表示寫(xiě)不下后面繼續(xù)p110,-10,-20...

p2-10,10,-20...

p3-10,-10,0plotzonegenzonebricksize6,8,8p0-10,-10,-20...;不規(guī)則六面體p110,-10,-20p2-10,10,-20...p3-10,-10,0p410,10,-20...p5-10,10,10p610,-10,0...p710,10,10plotcurrentplotPlot01plotclearplotzoneUndo;撤銷(xiāo)命令setlogfile127X1001.txtsetlogontruncatesetlogofflistzoneprinrangex01y01z01;顯示指定范圍內(nèi)各單元的主應(yīng)力，結(jié)果如下Hist命令：①命令編號(hào)按順序從1開(kāi)場(chǎng)：eachhistoryisnumberedsequentiallyfrom1asitisenteredviatheHISTORYcommand.

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

listhist

foralistingofthehistoriesandtheircorrespondingnumbers.histnstep5;每5步記錄1次。默認(rèn)是10步記錄1次histaddratio;記錄maximumunbalancedforceratiohistaddgpzdisp448;記錄節(jié)點(diǎn)的位移modelnullrangex=2,4y=2,6z=5,10;Toexcavatethetrench,enter

Thedimensionsoftheexcavationregion(i.e.,thenullmaterial)aredefinedbythex-,y-andz-coordinaterange.setlarge;設(shè)置大變形large-strainlogicisspecifiedinixdis0ydis0zdis0;位移歸零solvesavetrench;保持計(jì)算結(jié)果的文件“TRENCH.F3SAV〞不需要后綴4、group命令(注意新版本的區(qū)別，groupzone)genzonebricksize666

groupzoneTunnelrangex16y06z15;設(shè)置group

modelmechmohr;莫爾庫(kù)倫模型

modelmechnullrangegroupTunnel;groupzoneTunnelrangecylinderend1100105100end210095100radius3

groupzoneMatrixrangegroupTunnelnot

;namingRANGESthatarefixedinspaceandunrelated

;tothemodelgeometry-representsstratigraphy...

rangenameLayer1planedip0dd0ori000above

rangenameLayer2planedip0dd0ori000belowpropbulk1e8shear0.5e8coh0tens0fric35rangenrangeLayer1

propbulk1e7shear0.3e7coh1e7tens0fric0rangenrangeLayer2plzonecolorbygroup4、Interface命令以下是切割模型法生成interface：寺河模擬用的plzonegenzonebrickp0000p120000p20400p3004size100202group1;genzonebrickp0004p120004p20404p3008size100202group2;genseparategroup1interface1wrapfirstgroup1secondgroup2interface1maxedge0.5顯示interface：plotcreateplotPlot01plotaddsurface;顯示外表plotaddinterfacered;顯示交互面為紅色plotshow用導(dǎo)來(lái)導(dǎo)去法，已經(jīng)在5.0上運(yùn)行過(guò)，很好，值得借鑒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//設(shè)置采用間隔為5等效于sethist_rep=5histaddid=2zoneszz0,0,-1histunbal//最大非平衡力histaddgpzdisp448printhisthistaddgpzdis0,0,0histaddgpzvel0,0,0histaddgpydisp0,0,0;監(jiān)測(cè)網(wǎng)格點(diǎn)〔0，0，0〕點(diǎn)的y方向上的位移histaddzonesyy0,1,0;監(jiān)測(cè)單元〔0，1，0〕點(diǎn)的y方向上的應(yīng)力histaddzonesyy1,1,0;監(jiān)測(cè)單元〔1，1，0〕點(diǎn)的y方向上的應(yīng)力histaddzoneszz0,0,-1Listhistorysetlarge

historynstep2

historyaddfishtt

historyaddfishtheta

historyaddzonesxx(1,1,1)

historyaddzonesyy(1,1,1)

historyaddzoneszz(1,1,1)plotaddhis345vs1HISTORYaddfishvar;fishvar變量的監(jiān)測(cè)5、solve命令?????什么意思5、plot命令通過(guò)圖形操作界面的方式：plzonecolorbygr;查看各個(gè)分組通過(guò)命令流的方式：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這樣我們就建設(shè)了一個(gè)叫“GravV〞的圖形，將它設(shè)為當(dāng)前圖層，對(duì)此圖形設(shè)置一個(gè)傾角為90°的面，〔在x-y面上，約定z的負(fù)方向?yàn)橄隆?，傾向?yàn)?°。〔在x-y面上，從y的正方向開(kāi)場(chǎng)沿順時(shí)針?lè)较蛴?jì)量〕，過(guò)x-y面上一點(diǎn)〔x=3,y=4,z=0〕,在平面后加一個(gè)網(wǎng)格，垂直方向的應(yīng)力Qzz，最后給模型加軸線，整體圖形輪廓與添加圖形相反，顯示每個(gè)區(qū)域中心的應(yīng)力計(jì)算值，每個(gè)區(qū)域的顏色直接對(duì)應(yīng)根本區(qū)域的應(yīng)力，當(dāng)鍵入show,我們建設(shè)的圖形被顯示PlotcreateDispcontPladdcontdispplanebehindPladdaxesPlshowplotsketch;繪出輪廓線plotaddvel;添加矢量plotaddhis12vs3;繪出歷史記錄n，actLoad，p_solcyc15000plotCreateTrenchAddsurfaceyellowAddaxesblackShowplotcreateszz_contour;創(chuàng)立文件名為szz_contour〔z方向的應(yīng)力云圖〕plotsetcent445;視圖的中心坐標(biāo)為〔4，4，5〕plotsetrot20030;視圖的旋轉(zhuǎn)角度為〔20，0，30〕plotsetmag1.0;視圖的放大倍數(shù)為1.0plotaddcontszzoutonshadeon;翻開(kāi)szz的等值線圖plotaddaxes;翻開(kāi)坐標(biāo)軸線plotshow;顯示窗口;不同的區(qū)域，材料賦值不同plzonecolorbymodelplzonecolorbygroup求解前對(duì)模型進(jìn)展檢查網(wǎng)格檢查(如前所述)邊界條件檢查速度約束條件 plotgpfix荷載條件 plotfap模型檢查模型賦值 ploblockmodel參數(shù)賦值 ploblockprop***初始應(yīng)力檢查(如后所述)經(jīng)常檢查模型的響應(yīng)plozconszz(syy,sxx) 應(yīng)力場(chǎng)ploconzdis(ydis,xdis) 位移場(chǎng)plo(add)zonecolorbystate 屈服狀態(tài)plogpfix 速度約束條件plofap 體力plohist(unbal)1 不平衡力〔必需的〕plsetmovieactiveonplsetmovieextensionjpgplsetmovieinterval20plsetmovieprefixsihekuangkaiwaplsetmoviesizeixiyplsetnamesihedacaigaoplbitmapview1filemyplot.pngplbitmapview1size1024768plcreate<viewid>plcurrent<viewid>plcutaddplaneddxxxdipxxxplcutaddplanenamexxxnormalxxxoriginxxxpldxf<viewid><filefilename.dxf>plexcel<filefilename.xlsx>plsvg<viewid><filefilename.svg>plsetbackgroundwhite;設(shè)置背風(fēng)光plotaddhis345vs15、網(wǎng)格的連接、對(duì)稱(chēng)……不同尺寸的zone，網(wǎng)格用attach連接，粘起來(lái)?。。。。。?、plane命令6、range命令listrangerangenametrenchx=-1,1y=0,4z=-2,2;定義rangemodelmechnullrangenrangetrench;;開(kāi)挖excavatethetrench7、gp命令listgpposition;查看節(jié)點(diǎn)編號(hào)、位置坐標(biāo)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更多節(jié)點(diǎn)變量函數(shù)，詳見(jiàn)fishreference8、zone命令listzoneprintzonestressrangeid2381anyid2361any2341any;輸出網(wǎng)格單元的應(yīng)力histzoneszz0,0,-1ZPNTfindzone(INT)INTzfind(INT)ZPNTzoneheadii=nzonetotalnumberofzones更多zone變量函數(shù)，詳見(jiàn)fishreference9、config9、邊界條件BoundaryandInitialConditionsBoundaryconditionsarenormallyspecifiedwiththeAPPLYandFIX/FREEcommands,andinitialconditionsareassignedwiththeINITIALcommandExample：aconstant,compressivexx-stresscomponentof10MPacanbeappliedtoaboundarylocatedatx=10withthecommand

applysxx-10e6rangex10;x=10的面上10MPaapplysyy=-20e6grad0,0,20e5rangey-20z0,10在斜面上施加法向力10雜項(xiàng)〔材料力學(xué)性質(zhì)〕11.M-C本構(gòu)模型的選取在進(jìn)展模擬時(shí)，采用摩爾-庫(kù)侖〔Mohr-Coulomb〕強(qiáng)度準(zhǔn)則來(lái)判斷礦體與頂?shù)装宓钠茐臋C(jī)理。該準(zhǔn)則根本的內(nèi)涵為：巖石的剪切破壞是由巖石破裂面產(chǎn)生的破壞剪應(yīng)力受到巖石材料的內(nèi)聚力和內(nèi)摩擦力的抵抗，當(dāng)破裂面上的破壞剪應(yīng)力大于巖石的內(nèi)聚力和內(nèi)摩擦力時(shí)，即Fs>0時(shí)，材料發(fā)生剪切破壞[4]：有限差分法(FDM)Fish函數(shù)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);;重力加速度f(wàn)ixxyzrangez0.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);;重力加速度f(wàn)ixxyzrangez0.0fixyrangey0.0fixyrangey3.0fixxrangex0.0fixxrangex10.0setlarge;;大變形histaddunbalsolvesavecab_str;;保持計(jì)算結(jié)果文件，不用加后綴名.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