flac3d5.0常用命令集錦

精品文檔建 模1、調(diào)用文件：文件與工程在同一個文件夾，只寫文件名即可：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、創(chuàng)建 旋轉(zhuǎn) 縮放視圖3、建模命令model mech mohr ;莫爾庫倫模型model mech elastic ;彈性模型set grav 0, 0, -9.81 ; 重力加速度 negative z-direction. （垂直向下！常用的）下面這代碼，是沿著-y方向的重力加速度，注意區(qū)別！gen zone brick size 6,8,8 p0 -10, -10, -20 . ;省略號表示寫不下后面繼續(xù)p1 10, -10, -20 .p2 -10, 10, -20 .p3 -10, -10, 0 plot zonegen zone brick size 6,8,8 p0 -10, -10, -20 . ;不規(guī)則六面體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 0 1 z 0 1 ;顯示指定范圍內(nèi)各單元的主應(yīng)力，結(jié)果如下Hist命令： 命令編號按順序從1開始：each history is numbered sequentially from 1 as it is entered via the HISTORY command.查找顯示所有的his命令：Return to the Flac3D prompt and type list hist for 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 ;記錄節(jié)點的位移 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 ; 設(shè)置大變形large-strain logic is specifiedini xdis 0 ydis 0 zdis 0 ;位移歸零solvesave trench ; 保持計算結(jié)果的文件 “TRENCH.F3SAV” 不需要后綴 4、group命令(注意新版本的區(qū)別 ，group zone)gen zone brick size 6 6 6group zone Tunnel range x 1 6 y 0 6 z 1 5 ;設(shè)置groupmodel mech mohr ; 莫爾庫倫模型model mech null range group Tunnel ;group zone Tunnel range cylinder end1 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.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 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用導來導去法，已經(jīng)在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.1pl add interfacepauseini zpos add -1.0 range group Top4、hist命令(5.0版本需要加add！)hist nstep=5 /設(shè)置采用間隔為5 等效于set hist_rep=5hist add id =2 zone szz 0,0,-1 hist unbal /最大非平衡力hist add gp zdisp 4 4 8 print histhist add gp zdis 0,0,0hist add gp zvel 0,0,0hist add gp ydisp 0,0,0;監(jiān)測網(wǎng)格點（0，0，0）點的y 方向上的位移hist add zone syy 0,1,0;監(jiān)測單元（0，1，0）點的y 方向上的應(yīng)力hist add zone syy 1,1,0;監(jiān)測單元（1，1，0）點的y 方向上的應(yīng)力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 1HISTORY add fish var ; fish var變量的監(jiān)測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 whiteplot 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 illustrateplot add bound behind plot add bcont szz plane plot add axes plot show這樣我們就建立了一個叫“Grav V”的圖形，將它設(shè)為當前圖層，對此圖形設(shè)置一個傾角為90的面，（在x-y面上，約定z的負方向為下），傾向為0。（在x-y面上，從y的正方向開始沿順時針方向計量），過x-y面上一點（x=3,y=4,z=0）,在平面后加一個網(wǎng)格，垂直方向的應(yīng)力Qzz，最后給模型加軸線，整體圖形輪廓與添加圖形相反，顯示每個區(qū)域中心的應(yīng)力計算值，每個區(qū)域的顏色直接對應(yīng)基本區(qū)域的應(yīng)力，當鍵入show,我們建立的圖形被顯示 Plot create DispcontPl add cont disp plane behindPl add axesPl show plot sketch;繪出輪廓線plot add vel;添加矢量plot add his 1 2 vs 3;繪出歷史記錄n，actLoad ， p_solcyc 15000plot Create TrenchAdd surface yellowAdd axes blackShowplot create szz_contour;創(chuàng)建文件名為szz_contour（z 方向的應(yīng)力云圖）plot set cent 4 4 5;視圖的中心坐標為（4，4，5）plot set rot 20 0 30;視圖的旋轉(zhuǎn)角度為（20，0，30）plot set mag 1.0;視圖的放大倍數(shù)為1.0plot add cont szz out on shade on;打開szz 的等值線圖plot add axes;打開坐標軸線plot show;顯示窗口;不同的區(qū)域，材料賦值不同pl zone colorby modelpl zone colorby group求解前對模型進行檢查網(wǎng)格檢查(如前所述)邊界條件檢查速度約束條件 plot gpfix 荷載條件 plot fap 模型檢查模型賦值 plo block model參數(shù)賦值 plo block prop *初始應(yīng)力檢查(如后所述)經(jīng)常檢查模型的響應(yīng)plo zcon szz (syy, sxx)應(yīng)力場plo con zdis (ydis, xdis)位移場plo (add) zone colorby state屈服狀態(tài)plo gpfix 速度約束條件plo fap 體力plo hist (unbal) 1不平衡力 （必需的）pl set movieactive onpl set movieextension jpgpl set movieinterval 20pl set movieprefix sihekuangkaiwapl set moviesize ix iypl set name sihedacaigaopl bitmap view 1 file myplot.png pl bitmap view 1 size 1024 768pl create pl current pl cut add plane dd xxx dip xxx pl cut add plane name xxx normal xxx origin xxx pl dxf pl excel pl svg pl set background white ;設(shè)置背景色plot add his 3 4 5 vs 15、網(wǎng)格的連接、對稱不同尺寸的zone，網(wǎng)格用attach連接，粘起來！6、plane命令6、range命令list rangerange name trench x=-1,1 y=0,4 z=-2,2 ;定義rangemodel mech null range nrange trench ; ; 開挖excavate the trench7、gp命令list gp position ; 查看節(jié)點編號、位置坐標print gp disp range id 105hist add gp zdis 0,0,0hist add gp zvel 0,0,0message(1)=string(gp_xpos(pnt) message(1)=message(1)+,+string(gp_zpos(pnt) message(1)=message(1)+,+string(abs(gp_zdisp(pnt)GP PNT find gp(INT)GP PNT gp find(INT)GP PNT gp headii = ngp total number of gridpoints更多節(jié)點變量函數(shù)，詳見fish reference8、zone命令list zoneprint zone stress range id 2381 any id 2361 any 2341 any ;輸出網(wǎng)格單元的應(yīng)力hist zone szz 0,0,-1Z PNT find zone(INT)INT z find(INT)Z PNT zone headii = nzone total number of zones更多zone變量函數(shù)，詳見fish reference9、config 9、邊界條件 Boundary and Initial ConditionsBoundary conditions are normally specified with the APPLY and FIX/FREE commands, and initial conditions are assigned with the INITIAL commandExample： a constant, compressive xx-stress component of 10 MPa can be applied to a boundary located at x = 10 with the commandapply sxx -10e6 range x 10 ; x=10的面上 10MPaapply syy = -20e6 grad 0,0,20e5 range y -20 z 0,10在斜面上施加法向力10 雜項（材料力學性質(zhì)）11.M-C本構(gòu)模型的選取 在進行模擬時，采用摩爾-庫侖（Mohr-Coulomb）強度準則來判斷礦體與頂?shù)装宓钠茐臋C理。該準則基本的內(nèi)涵為：巖石的剪切破壞是由巖石破裂面產(chǎn)生的破壞剪應(yīng)力受到巖石材料的內(nèi)聚力和內(nèi)摩擦力的抵抗，當破裂面上的破壞剪應(yīng)力大于巖石的內(nèi)聚力和內(nèi)摩擦力時，即 Fs0 時，材料發(fā)生剪切破壞4：有限差分法(FDM).Fish函數(shù)1、 if條件從句2、 loop條件從句 newdef derives_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.25derivelist b_mod s_modgen zone brick size (2,2,2)model mech elasticprop bulk=b_mod shear=s_modlist zone prop bulklist zone prop shearnewdef xxxsum = 0prod = 1loop n (1,10)sum = sum + nprod = prod * nend_loopendxxxlist sum prodnewdef abcif xx 0 thenabc = 33elseabc = 11end_ifendset xx=1list abcset xx=-1list abcnewgen zone brick size 10 3 5def place_cablesloop n (1,5)z_d = float(n) - 0.5commandsel cable beg (0.0,1.5,z_d) end (7.0,1.5,z_d) nseg 7end_commandend_loopendplace_cablesplot add zone trans 75 plot add sel geomnewgen zone brick size 10 3 5model mech mohrprop bulk 1e8 shear 0.3e8 fric 35prop coh 1e3 tens 1e3ini dens 1000set grav (0,0,-10) ; ; 重力加速度 fix x y z range z 0.0fix y range y 0.0fix y range y 3.0fix x range x 0.0fix x range x 10.0set largehist add unbalsolvesave cab_strini xdis 0 ydis 0 zdis 0hist add gp xdisp (0,1,5)def place_cablesloop n (1,5)z_d = 5.5 - float(n)z_t = z_d + 0.5z_b = z_d - 0.5commandfree x range x 0.0 z z_b,z_tsolvesel cable beg (0.0,0.5,z_d) end (7.0,0.5,z_d) nseg 7sel cable beg (0.0,1.5,z_d) end (7.0,1.5,z_d) nseg 7sel cable beg (0.0,2.5,z_d) end (7.0,2.5,z_d) nseg 7sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 .gr_k 2e10 gr_coh 1e10 gr_per 1.0end_commandend_loopendplace_cablessave cab_endplot zone trans 75 sel cabblock force line width 4newgen zone brick size 1 2 1model mech mohrprop shear=1e8 bulk=2e8 .cohes=1e5 tens=1e10fix x y z range y 0.0apply yvel -1e-5 range y 2.0def get_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_addef loadload=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4)endhist add fish loadhist add gp ydis (0,6,0)step 1000plot his 1 vs -2newgen zone brick size 10 3 5model mech mohrprop bulk 1e8 shear 0.3e8 fric 35prop coh 1e3 tens 1e3ini dens 1000set grav (0,0,-10) ; ; 重力加速度fix x y z range z 0.0fix y range y 0.0fix y range y 3.0fix x range x 0.0fix x range x 10.0set large ; ;大變形hist add unbalsolvesave cab_str ; ;保持計算結(jié)果文件，不用加后綴名.Savini xdis 0 ydis 0 zdis 0 ; ; 位移歸零hist add gp xdisp (0,1,5)def place_cablesloop n (1,5)z_d = 5.5 - float(n)z_t = z_d + 0.5z_b = z_d - 0.5commandfree x range x 0.0 z z_b,z_tsolvesel cable beg (0.0,0.5,z_d) end (7.0,0.5,z_d) nseg 7sel cable beg (0.0,1.5,z_d) end (7.0,1.5,z_d) nseg 7sel cable beg (0.0,2.5,z_d) end (7.0,2.5,z_d) nseg 7sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 .gr_k 2e10 gr_coh 1e10 gr_per 1.0end_commandend_loopendplace_cablessave cab_endplot zone trans 75 sel cabblock force line width 4def parmglobal rad=4.0 ; radius of spherical cavityglobal len=10.0 ; length of outer box edgeglobal in_size=6 ; number of zone