flac學(xué)習(xí)吐血經(jīng)驗

1、精選優(yōu)質(zhì)文檔-傾情為你奉上第一天 看版主的初學(xué)者建議！ 初學(xué)者將FLAC/FLAC3D版所有帖子都瀏覽一遍；學(xué)有所得后，再瀏覽一遍，會發(fā)現(xiàn)又有新的收獲。 Flac3D中為巖土工程問題的求解開發(fā)了特有的本構(gòu)模型，總共包含了10種材料模型： 1. 開挖模型 2. 3個彈性模型（各向同性，橫觀各向同性和正交各向同性彈性模型） 3. 6個塑性模型（Drucker-Prager模型、Morh-Coulomb模型、應(yīng)變硬化/軟化模型、遍布節(jié)理模型、雙線性應(yīng)變硬化/軟化遍布節(jié)理模型和修正的cam粘土模型）。 因為我是學(xué)采礦專業(yè)的，我就先學(xué)習(xí)一下開挖模型吧！ 第二天 學(xué)習(xí)建模 FLAC3D建模的一般步驟如下：

2、 一、建模 new ge zo b po x y z & p1 x y z &# Z/ E7 g/ L( F) B% A p2 x y z &7 F. m$ U- g5 h' S7 u p3 x y z &1 i* d1 K% C; c9 y6 R; a size x y z plotShow add surface color add axes color 這樣一個立體模型就建立出來了！ 2、 定義材料參數(shù)mo moh/null/. pro bu= sh= fric= c= ten= 其中bulk是指體積模量也就是K=E/3(1-2v),E是彈模 v是

3、泊松比，當(dāng)v接近0.5時，不能盲目使用；shear是指剪切模量也就是G=E/2(1+v）； cohesion是指內(nèi)聚力； tension是指抗拉強(qiáng)度。 3、 加載及邊界 set g 0 0 -10 in de=fix x/y/z ra x/y/z 其中要有一個面為自由面。四、求解7 R( : m9 s% n C. fh n=10h unb h gp x/y/zdisp x y z set mec fo 30so 此為顯式求解過程。 五、結(jié)果分析監(jiān)控變量：print h plot history 1history 2位移等值線圖： plot cr add con disp/. add axes

4、color show 應(yīng)力等值線圖： Clearadd bcon szz/. add axes剖面上的等值線圖： cr set p dip= dd= o=x,y,z add bo beadd bcon szz/. Padd axesshow 有這些命令可以查看你想分析的FLAC圖 七、開挖 pro c= ten= mo null ra x=, y=, z=,set large/small1in xdisp=0 ydisp=0 zdisp=0 系統(tǒng)中的所有節(jié)點位移全部清零。8、 求解與后處理step plot cr dispcont plot copy dispcont settingsplot

5、 add con disp p be shade on plot add axes plot show開挖模型終于學(xué)習(xí)完了，經(jīng)過學(xué)習(xí)我可以建立一個簡單的開挖模型并計算、分析了。 下一步的工作：我想建立一個曲面的立體模型來反映礦山地形地貌。研究中. 經(jīng)過學(xué)習(xí)，我可以建立一個分層的地質(zhì)模型，怎樣把地表設(shè)置成曲面呢？形成礦體上方的地形地貌! 一個分層的地質(zhì)模型： 命令流如下：new gen zone b size 40 10 1 p0 0 0 0 p1 400 0 0 p2 0 200 0 p3 0 0 1.21 group 泥巖1 gen zone b size 40 10 2 p0 0 0 1.

6、21 p1 400 0 1.21 p2 0 200 1.21 p3 0 0 7.07 group 灰?guī)r1 gen zone b size 40 10 3 p0 0 0 7.07 p1 400 0 7.07 p2 0 200 7.07 p3 0 0 13.47 group 粉砂巖1 gen zone b size 40 10 1 p0 0 0 13.47 p1 400 0 13.47 p2 0 200 13.47 p3 0 0 14.96 group 灰?guī)r2 gen zone b size 40 10 4 p0 0 0 14.96 p1 400 0 14.96 p2 0 200 14.96 p3

7、 0 0 31.4 group 粉砂巖2gen zone b size 40 10 1 p0 0 0 31.4 p1 400 0 31.4 p2 0 200 31.4 p3 0 0 32.46 group 泥巖2 gen zone b size 40 10 3 p0 0 0 32.46 p1 400 0 32.46 p2 0 200 32.46 p3 0 0 39.73 group 粉砂巖3 gen zone b size 40 10 6 p0 0 0 39.73 p1 400 0 39.73 p2 0 200 39.73 p3 0 0 51.93 group 細(xì)砂巖1 gen zone b

8、size 40 10 7 p0 0 0 51.93 p1 400 0 51.93 p2 0 200 51.93 p3 0 0 65 group 粉砂巖gen zone b size 40 10 3 p0 0 0 65 p1 400 0 65 p2 0 200 65 p3 0 0 67.93 group 10煤gen zone b size 40 10 4 p0 0 0 67.93 p1 400 0 67.93 p2 0 200 67.93 p3 0 0 71.97 group 泥巖gen zone b size 40 10 1 p0 0 0 71.97 p1 400 0 71.97 p2 0

9、200 71.97 p3 0 0 72.35 group 天然焦 gen zone b size 40 10 2 p0 0 0 72.35 p1 400 0 72.35 p2 0 200 72.35 p3 0 0 75.15 group 巖漿巖gen zone b size 40 10 2 p0 0 0 75.15 p1 400 0 75.15 p2 0 200 75.15 p3 0 0 78.85 group 細(xì)砂巖2 gen zone b size 40 10 25 p0 0 0 78.85 p1 400 0 78.85 p2 0 200 78.85 p3 0 0 132.6 group

10、泥巖4. gen zone b size 40 10 4 p0 0 0 132.6 p1 400 0 132.6 p2 0 200 132.6 p3 0 0 136.79 group 細(xì)砂巖3 gen zone b size 40 10 3 p0 0 0 136.79 p1 400 0 136.79 p2 0 200 136.79 p3 0 0 140.24 group 泥巖5, gen zone b size 40 10 3 p0 0 0 140.24 p1 400 0 140.24 p2 0 200 140.24 p3 0 0 146.29 group 粉砂巖5 gen zone b si

11、ze 40 10 4 p0 0 0 146.29 p1 400 0 146.29 p2 0 200 146.29 p3 0 0 150 group 泥巖6 plotShowadd surface greenadd axes red1 FLAC建模方法1.1 建模（1） 設(shè)計計算模型的尺寸（2） 規(guī)劃計算網(wǎng)格數(shù)目和分布（3） 安排工程對象（開挖、支護(hù)等） （4） 給出材料的力學(xué)參數(shù) （5） 確定邊界條件（6） 計算模擬 1.2 網(wǎng)格生成：Grid i,j 例如：grid 30,20 1.3 網(wǎng)格規(guī)劃： Gen x1,y1 x2,y2 x3,y3 x4,y4 例如：Gen 0,0 0,10 10,

12、20 20,0 1.4 分區(qū)規(guī)劃網(wǎng)格。例如：Gen xI1,yI1 xI2,yI2 xI3,yI3 xI4,yI4 i=1,10 j=1,21 （I區(qū)）Gen xII1,yII1 xII2,yII2 xII3,yII3 xII4,yII4 i=10,20 j=1,21 (II區(qū)） 1.5 特殊形狀的網(wǎng)格（1） 圓形 gen circle xc,yc rad（2） 弧線 gen arc xc,yc xb,yb theta（3） 直線 gen line x1,y1 x2,y2 （4） 任意形狀 tab 1 x1,y1, x2,y2, ¼,xn,yn, x1,y1 gen ta

13、b 11.6 賦給單元材料性質(zhì) mod e (彈性) prop d 1800e-6 bu 12.5 sh 5.77 i=1,20 j=1,10 prop d 2400e-6 bu 1250 sh 577 i=1,20 j=11,20mod m （彈塑性Mohr-Coulumb準(zhǔn)則）prop d 1800e-6 bu 12.5 sh 5.77 c 0 fri 20 ten 0.015 reg i,j 1.7 賦給模型邊界條件（1） 固定邊界 （結(jié)點） Fix x i=1, j=1,21 Fix y i=1,21 j=1（2） 施加邊界力 (結(jié)點) apply yf -10 i=1,21 j=21

14、 或 apply syy -10 i=1,21 j=21 或 apply xf -5 i=21, j=1,21 或 apply sxx -5 i=21, j=1,21 （3） 賦單元內(nèi)應(yīng)力 （單元） ini sxx -10 i=1,20 j=1,20ini syy -5 var 0 4 i=1,21 j=1,21. 1.8 計算 Set grav 9.81 Set large Step 1000 Save test.sav1.9 結(jié)果顯示 Plot grid 顯示網(wǎng)格Plot bo 顯示邊界 Plot plas 顯示塑性區(qū)Plot sig1 fi 顯示最大主應(yīng)力s1Plot sig2 fi 顯

15、示最小主應(yīng)力s2Plot sdif fi 顯示主應(yīng)力差(s1- s2) Plot str 顯示主應(yīng)力矢量場Plot xdis fi 顯示X方向位移Plot ydis fi 顯示Y方向位移 Plot disp 顯示位移矢量場1.10 保存與調(diào)用結(jié)果命令Call test.txt (或ca test.dat) 調(diào)用數(shù)據(jù) Save test.sav 保存結(jié)果 New 重新開始Rest test.sav 調(diào)用結(jié)果 Quit 退出程序 2. FLAC運(yùn)用技巧 2.1 模型尺寸 2.2 模擬開挖 Mod nu i=6,15 j=5,12 (或 region i,j) 2.3 模擬錨桿支護(hù)（端錨）struc

16、t cable begin grid i1,j1 end grid i2,j2 seg n prop 1 stru prop 1 e 2e5 yield 0.5 a 0.235e-3 sbond 0 kbond 0.0001 stru prop 1 sfri 30 peri 0.2723 den 7.5e-3 2.4 模擬錨桿支護(hù)（全長錨固）struct cable begin grid i,j end x,y seg n prop 2stru prop 1 e 2e5 yield 0.5 a 0.235e-3 sbond 0.42 kbond 5.37 stru prop 1 sfri 30

17、 perimeter 0.2723 den 7.5e-32.5 模擬錨桿支護(hù)（預(yù)應(yīng)力錨固）struct cable begin grid i1,j1 end i2,j2 seg n ten 10 prop 1struct cable begin node n1 end x,y seg n prop 2 stru prop 1 e 2e5 yield 0.5 a 0.235e-3 sbond 0. kbond 0.1 stru prop 1 sfri 0 perimeter 0.2723 den 7.5e-3stru prop 2 e 2e5 yield 0.5 a 0.235e-3 sb 0.

18、42 kb 5.37stru prop 2 sfri 30 perimeter 0.2723 den 7.5e-3 2.6 模擬斷層 Mod nu j=37 int 1 as from i1,j1 to ,i1,j2 bs from i2,j1 to i2,j2 int 1 kn 1000 ks 400 fri 25 c 0.01 t 0.0001 2.7 其它技巧（1） 網(wǎng)格優(yōu)化 （2） 先彈性、后塑性（3） 分步開挖（4） 動力學(xué)問題模擬. （5） 固流耦合問題模擬（6） 合理解釋結(jié)果（7） 多種形式輸出結(jié)果（8） 報告格式 第3、 四天學(xué)習(xí)彭文斌的FLAC 3D實用教程第5天 運(yùn)用FLA

19、C做一個小項目 第六天：分享一些FLAC3D命令的FAQ資料！ 1. 如何調(diào)整視圖的大小、角度？,答：綜合使用x, y, z, m, Shift鍵，配合使用Ctrl+R，Ctrl+Z等快捷鍵2. 最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始條件、邊界條件、后處理和求解的命令。3. 怎樣看模型的樣子？ 答：plo blo gro可以看到不同的group的顏色分布 4. 怎樣看模型的邊界情況？ 答：plo gpfix red sk5. 怎樣看模型的體力分布？答：plo fap red sk 6. 怎樣看模型的云圖？ 答：位移：plo con

20、dis (xdis, ydis, zdis) 應(yīng)力：plo con sz (sy, sx, sxy, syz, sxz) 7. 怎樣看模型的矢量圖？ 答：plo dis (xdis, ydis, zdis8. 怎樣看模型有多少單元、節(jié)點？ 答：plo info9. 怎樣輸出模型的后處理圖？答：File/Print type/Jpg file，然后選擇File/Print，將保存格式選擇為jpg文件10. 怎樣調(diào)用一個文件？ 答：File/call或者call命令 11. 如何施加面力？! 答：app nstress12. 如何進(jìn)行邊界約束？ 答：fix x ran （約束的是速度，在初始情況下

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

22、鍵跳過本次計算，自動進(jìn)入下一步20. Fish是什么東西？答：是FLAC3D的內(nèi)置語言，可以用來進(jìn)行參數(shù)化模型、完成命令本身不能進(jìn)行的功能21. initial 與 apply 有何區(qū)別？答：initial初始化命令，如初始化計算體的應(yīng)力狀態(tài)等 apply邊界條件限制命令，如施加邊界的力、位移等約束等。 initial的應(yīng)力狀態(tài)會隨計算過程的發(fā)生而發(fā)生改變，一般體力需要初始化，而apply施加的邊界條件不會發(fā)生變化。 22. 如何查看各個時段不平衡力的具體數(shù)值？ 答：采用his來記錄計算，包括位移應(yīng)力等 命令his unbal his gp(zone) zdis range (0 0 0) 或

23、者id=？ 導(dǎo)出數(shù)據(jù)命令his write n vs m begin 時步 end 時步 file filename.his.n表示紀(jì)錄的id m表示時步要導(dǎo)出不平衡力的具體數(shù)值 his unbal step or solve7 k: his write 1 vs step begin 1 end 1000 file 123.his使用上述命令就可以查看各個時步下的不平衡力的具體數(shù)值 23. 如何調(diào)用一些可選模塊？ 答：config dyn (fluid, creep, cppudm) 24 .如何在圓柱體四周如何施加約束條件？ 答：可以用fix . ran cylinder end1 end

24、2 radius r1 cylinder end1 end2 radius r2 not，其中r2<r1，其實就是選擇range的時候選兩個圓柱的差，即得到邊界。 25. 如何能把一個PLOT的圖像數(shù)據(jù)導(dǎo)出來以便用其他軟件繪圖？ 答：用set log on 命令，把數(shù)據(jù)導(dǎo)出來，轉(zhuǎn)到excel里處理一下，然后用surfer或者什么作圖軟件繪制就行了。 26. 關(guān)于dd、dip的理解： dip-direction angle, measured in the global xy-plane clockwise from the positive y-axis 這個是從y軸到x順時針方向所轉(zhuǎn)過

25、的角度；dip dip angle, measured in the negative z-direction from the global xy-plane 這個是xy平面與z軸負(fù)方向的夾角。 dip表示對稱參照面與xy平面的夾角，對稱參照面與xy平面的夾角在xy平面的投影是一條射線 ，dd表示y軸正向順時針到那條射線的夾角。首先應(yīng)該按照dd的方向大體確定這個面的朝向，dd 指的是從y軸正方向按順時針（clockwise）方向轉(zhuǎn)向所要確定面的法線方向在xy平面上的投影的夾角，然后再確定dip，dip指的是從xy平面轉(zhuǎn)向所要確定的平面的角27. 如何調(diào)用FISH文件？ 答：新手調(diào)用FISH文

26、件出錯的原因：是：文件名為 fishcall.fis.txt(因為在記事本里這個txt的后綴是默認(rèn)的）。修改方法是：將 fishcall.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表示pre

27、vious,即分別表示在現(xiàn)在和以前的循環(huán)當(dāng)中處于塑性狀態(tài)29. 繪制相關(guān)曲線命令的一些解釋1) 在plot hist m vs n的形式里,m代表y軸,n代表x軸(不管m,n的正負(fù))2) "-"表示對其值作"mirror",如果大家使用AUTOCAD的話就明白"mirror"的含義了,比如上面的 hist 1 gp ydisp的值是0100,那么vs -1就變成-1000.以此類推. 3) 為什么要這么做呢?主要是為了符合巖土工程的習(xí)慣需要,例如樁載荷試驗曲線,沉降曲線等都不是畫在第一象限內(nèi). 30. 如何查看剖面上任一點的位移值？p

28、rint gp disp range x() y( ) z( )或!print gp disp range id plot gp disp range x() y( ) z( )或plot gp disp range id 31. fix x y z range x -0.1 0.1 .a fix x range x -0.1 0.1.b 請問a式與b式有何區(qū)別？ a 表示 在x=0 這個平面 是固定的 b 表示在x=0 這個平面 只x方向不能發(fā)生運(yùn)動32. 如何顯示變形輪廓線的命令？ plo ske magf 10 其中10為放大系數(shù)2011年的第一天： 首先介紹一下FLAC3D4.0新版本的安