FLAC3d如何操作

1、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, syz, sxz)7. 怎樣看模型的矢量圖？答：pl

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

3、元的ID？答：用鼠標(biāo)雙擊單元的表面，可以知道單元的ID和坐標(biāo)14. 如何進行切片？答：plo set plane ori (點坐標(biāo)) norm (法向矢量) plo con sz plane (顯示z方向應(yīng)力的切片) 15. 如何保存計算結(jié)果？答：save +文件名.16. 如何調(diào)用已保存的結(jié)果？答：rest +文件名；或者File / Restore17. 如何暫停計算？答：Esc18. 如何在程序中進行暫停，并可恢復(fù)計算？答：在命令中加入pause命令，用continue進行繼續(xù)19. 如何跳過某個計算步？答：在計算中按空格鍵跳過本次計算，自動進入下一步20. Fish是什么東西？答：是F

4、LAC3D的內(nèi)置語言，可以用來進行參數(shù)化模型、完成命令本身不能進行的功能21. Fish是否一定要學(xué)？答：可以不用，需要的時候查Mannual獲得需要的變量就可以了22. FLAC3D允許的命令文件格式有哪些？答：無所謂，只要是文本文件，什么后綴都可以23. 如何調(diào)用一些可選模塊？答：config dyn (fluid, creep, cppudm)24 .如何在圓柱體四周如何施加約束條件？答：可以用fix . ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2r1，其實就是選擇range的時候選兩個圓柱

5、的差，即得到邊界。25.如何能把一個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 時鈍角。dd是y軸順時針轉(zhuǎn)到法線在xoy面上的投影線時，所轉(zhuǎn)過的角度。27.如何調(diào)用FISH文件？答：新手調(diào)用FISH文件出錯的原因：是：文件名為 fishcall.fis.txt(因為在記事本里這個txt的后綴是默認(rèn)的）。修改方法是：將 fis

6、hcall.fis.txt里面的內(nèi)容全部復(fù)制到新的記事本里，再按保存，保存時候?qū)⒈４骖愋透臑樗形募?*,*),這時你便可在文件名欄內(nèi)輸入 fishcall.fis，類型也就正確了。28.如何理解Tension-p shear-p tension-n shear-n表示的含義？答：在計算的循環(huán)里面,每個循環(huán)中,每個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的形式

7、里,m代表y軸,n代表x軸(不管m,n的正負(fù));-表示對其值作mirror,如果大家使用AUTOCAD的話就明白mirror的含義了30.如何查看剖面上任一點的位移值？答：print gp disp range x() y( ) z( )或print gp disp range idplot gp disp range x() y( ) z( )或plot gp disp range id31. a: fix x y z range x -0.1 0.1 b: fix x range x -0.1 0.1 請問a式與b式有何區(qū)別？答：a 表示 在x=0 這個平面 是固定的；b 表示在x=0 這

8、個平面 只x方向不能發(fā)生運動32.initial 與 apply 有何區(qū)別？答：初始化命令，如初始化計算體的應(yīng)力狀態(tài)等；apply邊界條件限制命令，如施加邊界的力、位移等約束等。initial的應(yīng)力狀態(tài)會隨計算過程的發(fā)生而發(fā)生改變，一般體力需要初始化，而apply施加的邊界條件不會發(fā)生變化。33.如何查看各個時段不平衡力的具體數(shù)值？答：采用his來記錄計算，包括位移應(yīng)力等 命令his unbalhis gp(zone) zdis range (0 0 0) 或者id=？導(dǎo)出數(shù)據(jù)命令his write n vs m begin 時步end 時步 +file filename.his要導(dǎo)出不平衡力

9、的具體數(shù)值his unbalstep 100000 or solvehis write 1 vs step begin 1 end 1000 file 123.his使用上述命令就可以查看各個時步下的不平衡力的具體數(shù)值34.如何顯示變形輪廓線的命令？plo ske magf 10 其中10為放大系數(shù)35、怎樣開始入門？答：初步學(xué)習(xí)首先看說明書 Users gurid的Section2 Getting started和Section 3 Problemsloving with FLAC36、哪些命令需要掌握？答：grid gen prop interface ini fix set hist s

10、tep solve pause cont sav 等，在學(xué)習(xí)的過程中遇到這些命令和不會的命令應(yīng)該查看說明書的Command Reference，里面講的很清楚，特別需要注意的是一定要仔細(xì)看一下Command Reference開頭的Optional Range Phrases，這是講范圍的，因為很多命令的施加都需要加到模型的一部分，這里就告訴你怎樣來表示這一部分。這一原則同樣適合與其他5個ITASCA軟件。還有一點就是六個軟件（FLAC FALC3D UDEC 3DEC PFC PFC3D）的range大不相37、FLAC和FLAC3D建模有什么區(qū)別？答：FLAC和FLAC3D基本的建模思想大

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

12、l 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的編號44、FLAC怎么調(diào)整視圖？答：ctrl+z 或shift+z ,然后點擊要放大部位的左上角和右下角就可以放大某個部位，如果要縮小，雙擊屏幕。45、判斷模型中網(wǎng)格點是一個還是兩個的方法 顯示模型網(wǎng)格點的ID，如果有兩個ID號，說明是兩個。

13、節(jié)理面上的網(wǎng)格點應(yīng)該有兩個ID號。46、改變模型應(yīng)力邊界時FLAC與3DEC的不同之處在3DEC中，當(dāng)改變應(yīng)力邊界時，模型表面的應(yīng)力直接被更換，F(xiàn)LAC中改變應(yīng)力邊界時，模型表面的應(yīng)力為兩者應(yīng)力相加。47、flac中影響運行時間的因素（1）單元材料，結(jié)構(gòu)面之間的剛度對比很大；（2）劃分單元尺寸對比很大，盡量避免長細(xì)比大于5：1的細(xì)長單元。48 FLAC中有多少中材料模型內(nèi)置11中本構(gòu)模型%_4 W/N,n R5q（1）零模型 用于開挖回填的區(qū)域（2）各向同性彈性模型（3）正交各向異性彈性模型（4）橫觀各向同性彈性模型（5）德魯克普拉格塑性模型（6）莫爾庫侖塑性模型 一般巖土體的力學(xué)行為（7）節(jié)

14、理化塑性模型 用于研究材料破壞后行為（如邊破穩(wěn)定，屈服煤礦開采）（8）應(yīng)變硬化/軟化模型庫侖塑性模型 薄板層狀材料破壞后力學(xué)行為（9）雙線性應(yīng)變硬化/軟化節(jié)理化塑性模型 臨近層狀土層開挖（10）雙屈服塑性模型（11）修正的劍橋粘土模型49、對不平衡力的理解flac中用最大不平衡力和典型內(nèi)力大的比值R表示模型的不平衡力的相對大小，該值為百分?jǐn)?shù)，且從不會減小到零，因此，模型不會達到絕對穩(wěn)定的平橫狀態(tài)，可以根據(jù)要求的精度設(shè)定R值，模型默認(rèn)的R值是1e-5。注意R值小僅表示所有網(wǎng)格點的力都平衡，但模型還可能發(fā)生塑性流動，通過 plot vel 命令查看網(wǎng)格點速度，一個低振幅隨即速度場是沒有塑性流動的平

15、衡狀態(tài)的可靠的指示器。50、簡單命令解釋：plot add cont disp out on 時時顯示云圖plot block group gpnum on 顯示節(jié)點編號plot inter red ske 看接觸面的效果.del range group 1 not 不刪除group信息 del range group 2 刪除group信息ini z add 數(shù)值大小 range group name 順z方向移動建立接觸面 interface 1 face range x(范圍）y(范圍） z(范圍） 或者 interface 1 face range plane orig(x y z 坐

16、標(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 這個平面是固定的 fixed boundaryfix x range x -0.1 0.1表示在x=0 這個平面只x方向不能發(fā)生運動roller boundary51、fish命

17、令的解釋：def aux1 定義函數(shù)loop n(1,nptab1) 開始一個循環(huán)，從1循環(huán)到變量nptab1rr=xtable(1,n) 從table 1 的x列中讀第n個數(shù)付給變量rrzz=ytable(1,n) 從table 1 的y列中讀第n個數(shù)付給變量zzloop k(1,nprof) 開始第二個循環(huán)，從1循環(huán)到變量nprofcase_of k 下面是一個選擇結(jié)構(gòu)，利用變量k判斷i F alfa=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

18、 1 當(dāng)k1的時候執(zhí)行下面語句xtable(n+10,k)=rr table n+10的x列第k個值等于變量rrytable(n+10,k)=-40.0 table n+10的y列第k個值等于變量-40case 2 當(dāng)k2的時候執(zhí)行下面語句xtable(n+10,k)=rr table n+10的x列第k個值等于變量rrytable(n+10,k)=0.0 table n+10的y列第k個值等于變量0end_case 結(jié)束選擇結(jié)構(gòu)end_loop 結(jié)束第二個循環(huán)end_loop 結(jié)束第一個循環(huán)end 函數(shù)定義結(jié)束aux1 執(zhí)行什么定義的函數(shù)52. 在固定面上施加位移時，fix與apply不能同

19、時使用，應(yīng)該用ini和fix,才能施加上。53.菜單 plot items/add/structure elements/geometry/結(jié)構(gòu)單元的形式，把對話框里的mag fac的值改成不同的數(shù)，就可以看到位移放大后的結(jié)構(gòu)單元。54.你可以用鍵盤上的上、下、左、右肩在窗口中移動模型。55.如何能把一個PLOT的圖像數(shù)據(jù)導(dǎo)出來以便用其他軟件繪圖？答：用set log on 命令，把數(shù)據(jù)導(dǎo)出來，轉(zhuǎn)到excel里處理一下，然后用surfer或者什么作圖軟件繪制就行了。56 在進行三軸壓縮時候，如何施加圍壓？答：使用apply sxx =-1e5 szz=-1e5 range cyl end1 0

20、,0,0 end2 0,2,0 radius 1，其中X和Z軸是在我們眼睛正對的面。57 在進行三軸壓縮時候，如何對軸向力進行位移控制？答：使用 ini yvel 1e-7 range y -.1 .1 ini yvel 1e-7 range y 1.9 2.1 即可，數(shù)值可以變化58 在進行邊坡計算的時候，為什么用sol, fos 命令運行不行？答：應(yīng)該輸入sol fos slope3dfos.sav associated(不要也行） 再用restore 命令調(diào)用即可。59 在FLAC3D計算中，如何記錄程序中的值，如果不想時時記錄該怎么弄？答：用 hist wtite n1 vs n2 f

21、ile.txt， 想間隔記錄可用set hist_rep n。60 想記錄接觸面上的相關(guān)力學(xué)參數(shù)該怎么辦？答：hist inter i nstress /ndis/sstress/sdisp61 接觸面設(shè)置參數(shù)包括哪幾個？答：coh dial fric kn ks ten.62 如何輸出單元格的數(shù)據(jù)？答；expgrid .flacsd63 如何在剖面的沉降云圖上增加網(wǎng)格、剖面上的矢量及坐標(biāo)值？ D Q Z9s1 e4Y U%q 答：plo add skeplo add dis planeplo add axe64 如何做切片？ 答：plot set plane ori x y z norm x

22、 y z65 舉一個隧道支護開挖的例子？ gen zone radbrick p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &size 3 5 5 7 rat 1 1 1 1.5 dim 1 4 2 fillmodel mohrprop bulk 1e8 shear .3e8 fric 35prop coh 1e3 tens 1e3fix x range x -0.1 0.1fix z range z -0.1 0.1fix y range y 9.9 10.1set grav 0 10 0ini dens 1000ini syy=0.0 gra

23、d 0,-10000,0ini sxx=0.0 grad 0,-5000,0ini szz=0.0 grad 0,-5000,0apply sxx=0.0 grad 0,-5000,0 range x 9.9,10.1apply szz=0.0 grad 0,-5000,0 range z 9.9,10.1set largehist unbalsolvesave tr_eq.sav; monitor trench excavationini xdis=0 ydis=0 zdis=0hist gp xdisp 1,0,0hist gp zdisp 0,0,2 ; excavation steps

24、range name exc1 x 0 1 y 0.0 0.8 z 0 2range name exc2 x 0 1 y 0.8 1.6 z 0 2range name exc3 x 0 1 y 1.6 2.4 z 0 2range name exc4 x 0 1 y 2.4 3.2 z 0 2range name exc5 x 0 1 y 3.2 4.0 z 0 2; excavation step 1model null range exc1step 10000save exc1.sav; excavation step 2model null range exc2step 1000sav

25、e exc2.sav; excavation step 3model null range exc3step 1000save exc3.sav; excavation step 4model null range exc4step 1000save exc4.sav; excavation step 5model null range exc5step 1000save exc5.sav; install cable supportrest exc2.savsel cable beg 1.0,0.4,1.5 end 5.0,0.4,1.5 nseg 4sel cable beg 1.0,0.

26、4,0.5 end 5.0,0.4,0.5 nseg 4sel cable beg 1.0,1.2,1.5 end 5.0,1.2,1.5 nseg 4sel cable beg 1.0,1.2,0.5 end 5.0,1.2,0.5 nseg 4sel cable prop emod 2e9 ytension 1e8 xcarea 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0model null range exc3step 1000save cab3.savsel cable beg 1.0,2.0,1.5 end 5.0,2.0,1.5 nseg 4sel ca

27、ble beg 1.0,2.0,0.5 end 5.0,2.0,0.5 nseg 4sel cable prop emod 2e9 ytension 1e8 xcarea 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0model null range exc4;step 1000save cab4.savsel cable beg 1.0,2.8,1.5 end 5.0,2.8,1.5 nseg 4sel cable beg 1.0,2.8,0.5 end 5.0,2.8,0.5 nseg 4sel cable prop emod 2e9 ytension 1e8 xc

28、area 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0model null range exc5step 1000save cab5.savplot create TrenchSupportplot set rot 63.89, 334.35, 3.4plot set cent 2.893, 3.446, 3.899plot set size 4.136plot add cont disp out onplot add sel cable forceplot add axes black66 舉一個三軸壓縮的例子？gen zone cyl p0 0 0 0 p1 1

29、0 0 p2 0 2 0 p3 0 0 1 size 4 5 4gen zone reflect norm 1,0,0gen zone reflect norm 0,0,19model mohrprop bulk 1.19e10 shear 1.1e10prop coh 2.72e5 fric 44 ten 2e5fix x y z range y -.1 .1fix x y z range y 1.9 2.1apply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius 1ini yvel 0 grad 0 -1e-4 0 range y 0 2sol1.FLAC程式使用