FLAC3D常見命令與使用技巧

1、.FLAC3D常見命令與使用技巧1、FLAC3D常見命令：1.FLAC3D是有限元程序嗎？答：不是！是有限差分法。2.最先需要掌握的命令有哪些？答：需要掌握 gen, ini, app, plo, solve等建模、初始條件、邊界條件、后處理和求解的命令。3.怎樣看模型的樣子？答：plo blo gro 可以看到不同的group 的顏色分布4.怎樣看模型的邊界情況？答：plo gpfix red5.怎樣看模型的體力分布？答：plo fap red6.怎樣看模型的云圖？答：位移：plo con dis (xdis, ydis, zdis)應(yīng)力： plo con sz (sy, sx,sxy, s

2、yz, sxz)7.怎樣看模型的矢量圖？答：plo dis (xdis, ydis, zdis)8.怎樣看模型有多少單元、節(jié)點(diǎn)？答：pri info9.怎樣輸出模型的后處理圖？答： File/Print type/Jpg file ，然后選擇 File/Print ，將保存格式選擇為jpe 文件10.怎樣調(diào)用一個(gè)文件？答：File/call 或者 call 命令10.如何施加面力？答： app nstress11.如何調(diào)整視圖的大小、角度？答：綜合使用x, y, z, m, Shift鍵，配合使用 Ctrl+R，Ctrl+Z等快捷鍵12.如何進(jìn)行邊界約束？答：fix x ran（約束的是速度，

3、在初始情況下約束等效于位移約束）13.如何知道每個(gè)單元的ID？答：用鼠標(biāo)雙擊單元的表面，可以知道單元的ID 和坐標(biāo)14.如何進(jìn)行切片？答： plo set plane ori (點(diǎn)坐標(biāo) ) norm (法向矢量 )plo con sz plane (顯示 z 方向應(yīng)力的切片 )15.如何保存計(jì)算結(jié)果？答：save +文件名 .16.如何調(diào)用已保存的結(jié)果？答：rest +文件名；或者 File / Restore17.如何暫停計(jì)算？答： Esc18.如何在程序中進(jìn)行暫停，并可恢復(fù)計(jì)算？答：在命令中加入pause 命令，用 continue 進(jìn)行繼續(xù)19.如何跳過某個(gè)計(jì)算步？答：在計(jì)算中按空格鍵跳

4、過本次計(jì)算，自動(dòng)進(jìn)入下一步1 / 7.20. Fish 是什么東西？答：是 FLAC3D的內(nèi)置語言，可以用來進(jìn)行參數(shù)化模型、完成命令本身不能進(jìn)行的功能21. Fish 是否一定要學(xué)？答：可以不用，需要的時(shí)候查Mannual 獲得需要的變量就可以了22.FLAC3D允許的命令文件格式有哪些？答：無所謂，只要是文本文件，什么后綴都可以23.如何調(diào)用一些可選模塊？答：config dyn (fluid, creep, cppudm)24 .如何在圓柱體四周如何施加約束條件？答：可以用 fix . ran cylinder end1 end2 radius r1 cylinder end1 end2r

5、adius r2 not，其中 r225.如何能把一個(gè) PLOT的圖像數(shù)據(jù)導(dǎo)出來以便用其他軟件繪圖？答：用 set log on 命令，把數(shù)據(jù)導(dǎo)出來，轉(zhuǎn)到 excel 里處理一下，然后用 surfer 或者什么作圖軟件繪制就行了。26.關(guān)于 dd、 dip 的理解：答：首先你要確定你定義面的法線方向， dip 是傾角的意思（ xoy 面和定義面的夾角），法線向上 dip 是銳角，法線向下 dip 時(shí)鈍角。 dd 是 y 軸順時(shí)針轉(zhuǎn)到法線在 xoy 面上的投影線時(shí)，所轉(zhuǎn)過的角度。27.如何調(diào)用 FISH文件？答：新手調(diào)用 FISH文件出錯(cuò)的原因是：文件名為因?yàn)樵谟浭卤纠镞@個(gè) txt 的后綴是默認(rèn)

6、的）。修改方法是：將里面的內(nèi)容全部復(fù)制到新的記事本里，再按保存，保存時(shí)候?qū)⒈４骖愋透臑樗形募?(*,*), 這時(shí)你便可在文件名欄內(nèi)輸入 fishcall.fis，類型也就正確了。28.如何理解 Tension-p shear-p tension-n shear-n表示的含義？答：在計(jì)算的循環(huán)里面 ,每個(gè)循環(huán)中 ,每個(gè) zone 都依據(jù) failure criterion 處于不同的狀態(tài) ,shear 和 tension 分別表示因受剪和受拉而處于塑性狀態(tài) .n 表示 now,p 表示 previous,即分別表示在現(xiàn)在和以前的循環(huán)當(dāng)中處于塑性狀態(tài) .29.繪制相關(guān)曲線命令的一些解釋在plot

7、 hist m vs n 的形式里 ,m 代表 y 軸,n 代表 x 軸 (不管 m,n 的正負(fù) );"-" 表示對其值作 "mirror", 如果大家使用 AUTOCAD的話就明白 "mirror" 的含義了30.如何查看剖面上任一點(diǎn)的位移值？答： print gp disp range x() y( ) z( )或 print gp disp range idplot gp disp range x() y( ) z( ) 或 plot gp disp range id b: fix x range x -0.1 0.1 請問 a

8、 式與 b 式有何區(qū)別？答： a 表示在 x=0 這個(gè)平面是固定的； b 表示在 x=0 這個(gè)平面只 x 方向不能發(fā)生運(yùn)動(dòng)2 / 7.32. initial 與 apply 有何區(qū)別？答：初始化命令，如初始化計(jì)算體的應(yīng)力狀態(tài)等；apply 邊界條件限制命令，如施加邊界的力、位移等約束等。 initial 的應(yīng)力狀態(tài)會(huì)隨計(jì)算過程的發(fā)生而發(fā)生改變，一般體力需要初始化，而 apply 施加的邊界條件不會(huì)發(fā)生變化。33.如何查看各個(gè)時(shí)段不平衡力的具體數(shù)值？答：采用 his 來記錄計(jì)算，包括位移應(yīng)力等命令his unbalhis gp(zone) zdis range (0 0 0)或者 id=？導(dǎo)出數(shù)

9、據(jù)命令 his write n vs m begin 時(shí)步 end 時(shí)步 +file filename.his要導(dǎo)出不平衡力的具體數(shù)值his unbal step 100000 or solve使用上述命令就可以查看各個(gè)時(shí)步下的不平衡力的具體數(shù)值34.如何顯示變形輪廓線的命令？plo ske magf 10其中 10 為放大系數(shù)35、怎樣開始入門？答：初步學(xué)習(xí)首先看說明書User's gurid 的 Section2 Gettingstarted 和 Section 3 Problemsloving with FLAC36、哪些命令需要掌握？答：grid gen prop interf

10、ace ini fix set hist step solve pause cont sav等，在學(xué)習(xí)的過程中遇到這些命令和不會(huì)的命令應(yīng)該查看說明書的Command Reference，里面講的很清楚，特別需要注意的是一定要仔細(xì)看一下 Command Reference開頭的 Optional Range Phrases，這是講范圍的，因?yàn)楹芏嗝畹氖┘佣夹枰拥侥Ｐ偷囊徊糠?，這里就告訴你怎樣來表示這一部分。這一原則同樣適合與其他 5 個(gè) ITASCA軟件。還有一點(diǎn)就是六個(gè)軟件（ FLAC FALC3D UDEC 3DECPFC PFC3D）的range 大不相同37、FLAC和 FLAC3D

11、建模有什么區(qū)別？答： FLAC和 FLAC3D基本的建模思想大致是一樣的，不同的是FLAC是首先用 GRID命令把整體的網(wǎng)格畫出來，然后用MODELL NULL命令把網(wǎng)格分成許多部分，再用GEN命令把這些部分一塊一塊搭接起來。然后再用GEN LINE，GEN ARC等命令畫出要開挖的邊界，在FLAC中產(chǎn)生網(wǎng)格的命令只有一個(gè)GRID，而且只能用一次。而FLAC3D是直接用 GEN命令一塊一塊搭接， GEN可以用很多次。38、怎樣查看模型？答： plot grid 可以查看網(wǎng)格， plot grid num 可以查看節(jié)點(diǎn)號(hào)。39、怎樣查看邊界條件？答：plot fix40、怎樣查看預(yù)先留設(shè)的開挖邊

12、界？答：plot grid mark 。3 / 7.41、怎么查看模型的云圖？答： pl syy (sxx xdis ydis sig1.) fill，注意，帶 fill 關(guān)鍵字就是查看云圖，不帶 fill 關(guān)鍵字就是查看等直線圖。42、怎樣查看模型的矢量圖？答：pl grid vel (disp) red43、怎樣查看 INTERFACE？答： pl grid iface id redprint iface id ; id 是 interface 的編號(hào)44、FLAC怎么調(diào)整視圖？答： ctrl+z 或 shift+z ,然后點(diǎn)擊要放大部位的左上角和右下角就可以放大某個(gè)部位，如果要縮小，雙擊

13、屏幕。45、判斷模型中網(wǎng)格點(diǎn)是一個(gè)還是兩個(gè)的方法顯示模型網(wǎng)格點(diǎn)的ID，如果有兩個(gè) ID 號(hào)，說明是兩個(gè)。節(jié)理面上的網(wǎng)格點(diǎn)應(yīng)該有兩個(gè)ID 號(hào)。46、改變模型應(yīng)力邊界時(shí) FLAC與 3DEC的不同之處在 3DEC中，當(dāng)改變應(yīng)力邊界時(shí)，模型表面的應(yīng)力直接被更換， FLAC中改變應(yīng)力邊界時(shí)，模型表面的應(yīng)力為兩者應(yīng)力相加。47、flac 中影響運(yùn)行時(shí)間的因素（ 1）單元材料，結(jié)構(gòu)面之間的剛度對比很大；（ 2）劃分單元尺寸對比很大，盡量避免長細(xì)比大于 5： 1 的細(xì)長單元。48 FLAC中有多少種材料模型內(nèi)置11 中本構(gòu)模型（ 1）零模型用于開挖回填的區(qū)域（ 2）各向同性彈性模型（ 3）正交各向異性彈性模

14、型（ 4）橫觀各向同性彈性模型（ 5）德魯克普拉格塑性模型（ 6）莫爾庫侖塑性模型一般巖土體的力學(xué)行為（ 7）節(jié)理化塑性模型用于研究材料破壞后行為（如邊破穩(wěn)定，屈服煤礦開采）（8）應(yīng)變硬化 / 軟化模型庫侖塑性模型薄板層狀材料破壞后力學(xué)行為 （9）雙線性應(yīng)變硬化 / 軟化節(jié)理化塑性模型臨近層狀土層開挖（ 10）雙屈服塑性模型（ 11）修正的劍橋粘土模型49、對不平衡力的理解flac 中用最大不平衡力和典型內(nèi)力大的比值R 表示模型的不平衡力的相對大小，該值為百分?jǐn)?shù)，且從不會(huì)減小到零，因此，模型不會(huì)達(dá)到絕對穩(wěn)定的平橫狀態(tài)，可以根據(jù)要求的精度設(shè)定 R 值，模型默認(rèn)的 R 值是 1e-5。注意 R 值

15、小僅表示所有網(wǎng)格點(diǎn)的力都平衡， 但模型還可能發(fā)生塑性流動(dòng)， 通過 plot vel 命令查看網(wǎng)格點(diǎn)速度，一個(gè)低振幅隨即速度場是沒有塑性流動(dòng)的平衡狀態(tài)的可靠的指示器。50、簡單命令解釋：plot add cont disp out on 時(shí)時(shí)顯示云圖plot block group gp num on 顯示節(jié)點(diǎn)編號(hào)plot inter red ske看接觸面的效果 .del range group 1 not 不刪除 group 信息4 / 7.del range group 2刪除 group 信息ini z add 數(shù)值大小range group name順 z 方向移動(dòng)建立接觸面 inte

16、rface 1 face range x(范圍） y(范圍） z(范圍）或者 interface 1 face range plane orig(x y z坐標(biāo)） norm （法線方向）water table clear 清除水位線pl set back white ;設(shè)置背景set plot bitmap size (800,510) ;設(shè)置圖片大小plot set caption size 37 設(shè);置圖例數(shù)字大小 (default=35 范圍 10-50)plot set center 100 0 40fix x y z range x -0.1 0.1表示在 x=0 這個(gè)平面是固定的

17、fixed boundaryfix x range x -0.1 0.1表示在 x=0 這個(gè)平面只 x 方向不能發(fā)生運(yùn)動(dòng) roller boundary51、fish 命令的解釋：def aux1 定義函數(shù)loop n(1,nptab1)開始一個(gè)循環(huán)，從1 循環(huán)到變量nptab1rr=xtable(1,n)從 table 1 的 x 列中讀第 n 個(gè)數(shù)付給變量rrzz=ytable(1,n)從 table 1 的 y 列中讀第 n 個(gè)數(shù)付給變量zzloop k(1,nprof)開始第二個(gè)循環(huán)，從1 循環(huán)到變量 nprofcase_of k下面是一個(gè)選擇結(jié)構(gòu)，利用變量k 判斷 i Falfa=0

18、.5*(pi+.1)*float(k-2)/float(nprof-2)xtable(n+10,k)=rr*cos(alfa)ytable(n+10,k)=rr*sin(alfa)case 1 當(dāng) k1 的時(shí)候執(zhí)行下面語句xtable(n+10,k)=rr table n+10的 x 列第 k 個(gè)值等于變量 rrytable(n+10,k)=-40.0 tablen+10 的 y 列第 k 個(gè)值等于變量 -40case 2當(dāng) k 2 的時(shí)候執(zhí)行下面語句xtable(n+10,k)=rr table n+10 的 x列第 k 個(gè)值等于變量 rrytable(n+10,k)=0.0 table n

19、+10的 y 列第 k 個(gè)值等于變量 0end_case結(jié)束選擇結(jié)構(gòu)end_loop結(jié)束第二個(gè)循環(huán)end_loop 結(jié)束第一個(gè)循環(huán) end 函數(shù)定義結(jié)束 aux1 執(zhí)行什么定義的函數(shù)52.在固定面上施加位移時(shí)，fix 與 apply 不能同時(shí)使用，應(yīng)該用ini 和 fix,才能施加上。5 / 7.53.菜單 plot items/add/structure elements/geometry/ 結(jié)構(gòu)單元的形式， 把對話框里的 mag fac的值改成不同的數(shù)，就可以看到位移放大后的結(jié)構(gòu)單元。54.你可以用鍵盤上的上、下、左、右肩在窗口中移動(dòng)模型。55.如何能把一個(gè) PLOT的圖像數(shù)據(jù)導(dǎo)出來以便用

20、其他軟件繪圖？答：用 set log on 命令，把數(shù)據(jù)導(dǎo)出來，轉(zhuǎn)到 excel 里處理一下，然后用 surfer 或者什么作圖軟件繪制就行了。56 在進(jìn)行三軸壓縮時(shí)候，如何施加圍壓？答：使用 apply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius，1其中 X 和 Z 軸是在我們眼睛正對的面。57 在進(jìn)行三軸壓縮時(shí)候，如何對軸向力進(jìn)行位移控制？答：使用 ini yvel 1e-7 range y -.1 .1 ini yvel 1e-7 range y 1.9 2.即1可，數(shù)值可以變化58 在進(jìn)行邊坡計(jì)算的時(shí)候，為什么用sol, fos 命令運(yùn)行不行？答：應(yīng)該輸入 sol fos slope3dfos.sav associated(不要也行）再用restore 命令調(diào)用即可。59 在 FLAC3D計(jì)算中，如何記錄程序中的值，如果不想時(shí)時(shí)記錄該怎么弄？答：用 hist wtite n1 vs n2 file.txt ，想間隔記錄可用set hist_rep n。60想記錄接觸面上的相關(guān)力學(xué)參數(shù)該怎么辦？答： hist inter i nstress /ndis/sstress/sdisp61接觸面設(shè)置參數(shù)包括哪幾個(gè)？答： coh dial fric kn k