FLACD常見命令與使用技巧

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,

2、sxy, syz, 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,乙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ì)算步？

4、答：在計(jì)算中按空格鍵跳過本次計(jì)算，自動進(jìn)入下一步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 cylinderendl end

5、2radius 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文件？28. 如何理解 Tension-p shear-p tension-n shear-n 表示的含義？答：在計(jì)算的循環(huán)

6、里面，每個(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 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 i

7、dplot gp disp range x() y( ) z( )或 plot gp disp range id31. a: fix x y z range x-0.10.1b: fix x range x-0.1 0.1 請問 a 式與 b 式有何區(qū)別？答：a表示在x=0這個(gè)平面是固定的；b表示在x=0這個(gè)平面只x方向不能發(fā)生運(yùn)動32.initial 與 apply 有何區(qū)別？答：初始化命令，如初始化計(jì)算體的應(yīng)力狀態(tài)等；apply邊界條件限制命令，如施加邊界的力、位移等約束等。initial的應(yīng)力狀態(tài)會隨計(jì)算過程的發(fā)生而發(fā)生改變，一般體力需要初始化，而apply施加的邊界條件不會發(fā)生變化。3

8、3. 如何查看各個(gè)時(shí)段不平衡力的具體數(shù)值？答：采用his來記錄計(jì)算，包括位移應(yīng)力等命令 his unbalhis gp(zone) zdis range (0 0或者id=？導(dǎo)出數(shù)據(jù)命令 his write n vs m beg in 時(shí)步 end 時(shí)步 +file file name.his要導(dǎo)出不平衡力的具體數(shù)值 his unbal step 100000 or solvehis write 1 vs step begin 1 end 1000 file.123.his使用上述命令就可以查看各個(gè)時(shí)步下的不平衡力的具體數(shù)值34. 如何顯示變形輪廓線的命令？ plo ske magf 10 其

9、中 10 為放大系數(shù)35、怎樣開始入門？答：初步學(xué)習(xí)首先看說明書 Users gurid 的 Section2 Gettingstarted 和 Section 3 Problemsloving with FLAC36、哪些命令需要掌握？答：grid gen prop in terface ini fix set hist step solve pause cont sa等，在學(xué)習(xí)的過程中遇至 U這些命令和 不會的命令應(yīng)該查看說明書的 Comma nd Refere nee里面講的很清楚，特別需要注意的是一定要仔細(xì) 看一下Comma nd Refere nee開頭的Optio nal Rang

10、e Phrases這是講范圍的，因?yàn)楹芏嗝畹氖┘佣夹?要加到模型的一部分，這里就告訴你怎樣來表示這一部分。這一原則同樣適合與其他5個(gè)ITASCA軟件。還有一點(diǎn)就是六個(gè)軟件(FLAC FALC3D UDEC 3DECPFC P)C的Dan ge大不相同37、FLAC和FLAC3D建模有什么區(qū)別？答：FLAC和FLAC3D基本的建模思想大致是一樣的，不同的是 FLAC是首先用GRID命令把整體的網(wǎng) 格畫出來，然后用MODELL NULI命令把網(wǎng)格分成許多部分，再用 GEN命令把這些部分一塊一塊搭接 起來。然后再用GEN LINE GEN ARC等命令畫出要開挖的邊界，在 FLAC中產(chǎn)生網(wǎng)格的命令

11、只有一個(gè) GRID而且只能用一次。而 FLAC3D是直接用GEN命令一塊一塊搭接，GEN可以用很多次。38、怎樣查看模型？答： plot grid 可以查看網(wǎng)格， plot grid num 可以查看節(jié)點(diǎn)號。39、 怎樣查看邊界條件？答：plot fix40、 怎樣查看預(yù)先留設(shè)的開挖邊界？答：plot grid mark。41、怎么查看模型的云圖？答： pl syy (sxx xdis ydis sig1.) fi，ll 注意，帶 fill 關(guān)鍵字就是查看云圖，不帶 fill 關(guān)鍵字就是查看等直線圖。42、 怎樣查看模型的矢量圖？答：pl grid vel (disp) red43、怎樣查看

12、INTERFACE 答：pl grid iface id redpri nt iface id ; id 是 in terface 的編號44、FLAC怎么調(diào)整視圖？答：ctrl+z或shift+z然后點(diǎn)擊要放大部位的左上角和右下角就可以放大某個(gè)部位，如果要縮小，雙 擊屏幕。45、 判斷模型中網(wǎng)格點(diǎn)是一個(gè)還是兩個(gè)的方法顯示模型網(wǎng)格點(diǎn)的ID,如果有兩個(gè)ID號，說明是兩個(gè)。 節(jié)理面上的網(wǎng)格點(diǎn)應(yīng)該有兩個(gè)ID號。46、改變模型應(yīng)力邊界時(shí)FLAC與3DEC的不同之處在3DEC中，當(dāng)改變應(yīng)力邊界時(shí)，模型表面的應(yīng)力直接被更換，F(xiàn)LAC中改變應(yīng)力邊界時(shí)，模型表面的應(yīng)力為兩者應(yīng)力相加。47、flac中影響運(yùn)行時(shí)

13、間的因素（1）單元材料，結(jié)構(gòu)面之間的剛度對比很大；（2）劃分單元尺寸對 比很大，盡量避免長細(xì)比大于 5：1 的細(xì)長單元。48 FLAC中有多少種材料模型內(nèi)置11中本構(gòu)模型（ 1）零模型用于開挖回填的區(qū)域（ 2）各向同性彈性模型（ 3）正交各向異性彈性模型（ 4）橫觀各向 同性彈性模型（ 5）德魯克普拉格塑性模型（ 6）莫爾庫侖塑性模型一般巖土體的力學(xué)行為（7）節(jié)理化塑性模型用于研究材料破壞后行為（如邊破穩(wěn)定，屈服煤礦開采）（ 8）應(yīng)變硬化 /軟化模型 庫侖塑性模型薄板層狀材料破壞后力學(xué)行為（ 9）雙線性應(yīng)變硬化 / 軟化節(jié)理化塑性模型臨近層狀土層開挖（ 10）雙屈服塑性模型（ 11）修正的劍橋

14、粘土模型49、對不平衡力的理解flac中用最大不平衡力和典型內(nèi)力大的比值R表示模型的不平衡力的相對大小，該值為百分?jǐn)?shù)，且從不會減小到零，因此，模型不會達(dá)到絕對穩(wěn)定的平橫狀態(tài)，可以根據(jù)要求的精度設(shè)定R值，模型默認(rèn)的R值是1e-5。注意R值小僅表示所有網(wǎng)格點(diǎn)的力都平衡，但模型還可能發(fā)生塑性流動，通過plotvel 命令查看網(wǎng)格點(diǎn)速度，一個(gè)低振幅隨即速度場是沒有塑性流動的平衡狀態(tài)的可靠的指示器。50、簡單命令解釋：plot add cont disp out on 時(shí)時(shí)顯示云圖plot block group gp num on 顯示節(jié)點(diǎn)編號plot inter red ske看接觸面的效果 .de

15、l range group 1 not不刪除 group 信息del range group 2刪除 group 信息ini z add 數(shù)值大小range group name順 z方向移動建立接觸面in terface 1 face range x范圍）y（范圍）z范圍）或者 in terface 1 face range pla ne 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 siz

16、e 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è)平面是固定的fixed boundary fix x range x-0.1 0.1表示在x=0這個(gè)平面只x方向不能發(fā)生運(yùn)動roller boundary51. fish 命令的解釋：def aux1 定義函數(shù)loop n(1, nptabl)開始一個(gè)循環(huán)，從1循環(huán)到變量nptab1rr=xtable(1,n)從table 1的x列中讀第n個(gè)數(shù)付給變量rrzz=ytable(1,n)從table 1的y列

17、中讀第n個(gè)數(shù)付給變量zzloop k(1,nprof)開始第二個(gè)循環(huán)，從1循環(huán)到變量nprofcase_of k下面是一個(gè)選擇結(jié)構(gòu)，利用變量 k 判斷 i Falfa=0.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) k= 1 的時(shí)候執(zhí)行下面語句 xtable(n+10,k)=rr table n+10的 x列第 k個(gè)值等于變量 rrytable(n+10,k)=-40.0 table n+10的y列第k個(gè)值等于變量-40case 2當(dāng)k= 2

18、的時(shí)候執(zhí)行下面語句 xtable(n+10,k)=rr table n+10的x 列第k個(gè)值等于變量rrytable(n+10,k)=0.0 table n+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,才能施加上。53. 菜單plot items/add/structure elements/geometry/結(jié)構(gòu)單元的形式，把對話框里的 mag fac的值改成 不同的數(shù)，就可以看

19、到位移放大后的結(jié)構(gòu)單元。54. 你可以用鍵盤上的上、下、左、右肩在窗口中移動模型。55. 如何能把一個(gè)PLOT勺圖像數(shù)據(jù)導(dǎo)出來以便用其他軟件繪圖？答：用 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

20、-.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 ks ten.62 如何輸出