FLAC3D命令流(挺實(shí)用)參考模板

1、1、怎樣查看模型？ 答：plot grid 可以查看網(wǎng)格，plot grid num 可以查看節(jié)點(diǎn)號。 2、請問在圓柱體四周如何施加約束條件？ 答：可以用fix . ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2<r1，其實(shí)就是選擇range的時(shí)候選兩個(gè)圓柱的差，即得到邊界。命令流如下：fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 13、怎么能把一個(gè)PLOT的圖像數(shù)據(jù)導(dǎo)出來以便用其他軟件繪圖？ 答：用set l

2、og on 命令，把數(shù)據(jù)導(dǎo)出來，轉(zhuǎn)到excel里處理一下，然后用surfer或者什么作圖軟件繪制就行了。4、用命令建立模型后，如何顯示點(diǎn)的坐標(biāo)？ 答：使用plo blo gro gpnum on命令 5、關(guān)于gauss_dev對性質(zhì)進(jìn)行高斯正態(tài)分布的問題？ 答：根據(jù)手冊上的說明：下面的命令設(shè)定一個(gè)平均摩擦角為40度，標(biāo)準(zhǔn)方差是±5。則命令如下：prop friction 40 gauss_dev 2 問題：請問gauss_dev 2中的2是如何計(jì)算的？如果把±5改為±10，則命令應(yīng)如何寫？40×526、 reflect問題 問：gen zone radb

3、rick &p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &size 3,5,5,7 &ratio 1,1,1,1.5 &dim 1 4 2 fill plot surf gen zone reflect dip 0 dd 90 （對xy面做鏡像）gen zone reflect dip 90 dd 90 （對yz面做鏡像）（1）dd表示y軸正向順時(shí)針到那條射線的夾角，dip表示對稱參照面與xy平面的夾角，對稱參照面與xy平面的夾角在xy平面的投影是一條射線。首先應(yīng)該按照dd的方向大體確定這個(gè)面的朝向，dd 指的是

4、從y軸正方向按順時(shí)針（clockwise）方向轉(zhuǎn)向所要確定面的法線方向在xy平面上的投影的夾角，然后再確定dip，dip指的是從xy平面轉(zhuǎn)向所要確定的平面的角（及z軸負(fù)方向轉(zhuǎn)過角度）（2）命令：reflect normal xyz origin xyz （根據(jù)法線和過一點(diǎn)建的面）最常用，至于實(shí)際使用過程中有人問道：normal -1 0 0 與normal 1 0 0的區(qū)別，我試過兩者的效果是一樣的，沒區(qū)別，雖然方向不同，但表達(dá)的是同一個(gè)面。1 / 25reflect dd a dip b origin xyz 這個(gè)就是用里建立斜面鏡像的。有一點(diǎn)注意了，鏡像的時(shí)候如果是模型鏡像的話，相當(dāng)于在建

5、立一個(gè)模型，故，在reflect的前面要加建模的關(guān)鍵字，gen zone reflect 7、 關(guān)于兩種fix是否有區(qū)別？fix z range z -0.1 0.1 fix x range x -0.1 0.1 any x 1.9 2.1 any fix y range y -0.1 0.1 any y 1.9 2.1 any 和 fix z range z -0.1 0.1 fix x range x -0.1 0.1 fix x range 1.9 2.1 fix y range y -0.1 0.1 fix y range y 1.9 2.1 答：是一樣的，any表示并列連接關(guān)系，及a

6、nd的意思。8、 initial apply set 的區(qū)別答：initial 初始化命令，初始化計(jì)算應(yīng)力狀態(tài)，定義材料參數(shù)、重力等；apply邊界條件限制命令，如施加邊界的力、位移等約束等；set 用于設(shè)置計(jì)算的模式，控制計(jì)算過程。Ini 施加的初始力在計(jì)算過程中是會(huì)變的，一般體力需要初始化，而apply施加的力、邊界條件不會(huì)隨著計(jì)算變化。set命令作用于整個(gè)網(wǎng)格(賦值全局變量),ini則初始化指定范圍內(nèi)網(wǎng)格點(diǎn)、域的偏移量。初始應(yīng)力不僅僅是重力引起的。9、Tension-p shear-p tension-n shear-n表示的含義 答：在計(jì)算的循環(huán)里面,每個(gè)循環(huán)中,每個(gè)zone 都依據(jù)f

7、ailure criterion 處于不同的狀態(tài),shear和tension 分別表示因受剪和受拉而處于塑性狀態(tài).n 表示now,p表示previous,即分別表示在現(xiàn)在和以前的循環(huán)當(dāng)中處于塑性狀態(tài)。10、history保存數(shù)據(jù)，print顯示在錯(cuò)誤的命令流中“原點(diǎn)位移.dat”文件中只有：“Step 1 Y-Displace” 而沒有保存下數(shù)值，怎么回事？我怎么才能獲得table中的數(shù)據(jù)？ 輸入'print tabel n'能得到table中的數(shù)據(jù)（正確）hist gp ydisp 0,0,0 hist zone syy 0,1,0 hist zone syy 1,1,0 s

8、tep 3000 hist write 1 file 原點(diǎn)位移.dat hist write 2 file 中部單元應(yīng)力.dat hist write 3 file 邊緣單元應(yīng)力.dat hist write 1 table 1 hist write 2 table 2 hist write 3 table 3hist write 2 3 vs 1 file 單壓.dat （錯(cuò)誤）hist gp ydisp 0,0,0 hist zone syy 0,1,0 hist zone syy 1,1,0 hist write 1 file 原點(diǎn)位移.dat hist write 2 file 中部單

9、元應(yīng)力.dat hist write 3 file 邊緣單元應(yīng)力.dat hist write 1 table 1 hist write 2 table 2 hist write 3 table 3 hist write 2 3 vs 1 file 單壓.dat step 300011、plot顯示x、y坐標(biāo)圖hist n 1 hist gp ydisp 0,0,0 hist ax_str hist gp xdisp 1,1,0 step 1500 plot hist 2 vs -1 ;axial stress vs axial disp. plot hist 3 vs -1 ; circum

10、ferential disp. vs axial disp答：那2是不是代表hist gp ydisp 0,0,0，而-1 指得是hist gp ydisp 0,0,0為什么不是1而是負(fù)值？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í)慣需要,

11、例如樁載荷試驗(yàn)曲線,沉降曲線等都不是畫在第一象限內(nèi)。12、 關(guān)于water.fis的命令注釋（設(shè)置地下水位）答：（1）water.fis的功能是設(shè)置地下水位面，在二維計(jì)算中，通過water table命令即可設(shè)置水位線，在三維計(jì)算中由于地下水位是面，因此需要作一些插值運(yùn)算，以便形成水位面。這個(gè)FISH程序包含三個(gè)函數(shù)parm,aux1和aux2。parm用來設(shè)定水位的邊值，即設(shè)置水位的區(qū)域，aux1插值計(jì)算等值線的點(diǎn)，aux2調(diào)用aux1的計(jì)算結(jié)果，使用FLAC3D命令形成水位面。（2）首先設(shè)置參數(shù)值，然后定義一個(gè)table. def parm nptab1=7 nprof=10 end pa

12、rm table 1 (12 0) (18.25 12.5) (25 17) (35 21) (50 24) (75 25) (150 25) 產(chǎn)生table1，包含7組數(shù)據(jù)點(diǎn)。 3) 在解釋的兩個(gè)FISH函數(shù)之前，對table,xtable,ytable作一些總的說明。water.fis程序反復(fù)使用了這三個(gè)變量，理解了這三個(gè)變量的含義，water.fis這段程序就解決了大半問題了。 Itasca 的所有軟件都使用了相同的數(shù)據(jù)結(jié)構(gòu)儲存數(shù)據(jù)，稱之為鏈表結(jié)構(gòu)，類似于關(guān)系型數(shù)據(jù)庫中的數(shù)據(jù)鏈接方法。如果有數(shù)據(jù)結(jié)構(gòu)的知識，或掌握關(guān)系型數(shù)據(jù)庫的SQL 語言，這三個(gè)變量的含義就不難理解了，只不過這個(gè)操作與真

13、正的數(shù)據(jù)庫操作相比，只有兩個(gè)數(shù)據(jù)變量x,y。我用數(shù)據(jù)庫的術(shù)語作解釋：table 就是產(chǎn)生一個(gè)“表”，xtable,ytable的作用就是產(chǎn)生一個(gè)新表(Create)以及對表中的數(shù)據(jù)進(jìn)行取出(Select)、插入(Insert)、替換(Update)操作. 4）FISH函數(shù)aux1的作用是在table 1的基礎(chǔ)上產(chǎn)生出另外個(gè)table，table的ID 從11到17，每個(gè)table有10個(gè)數(shù)據(jù)集。alfa的計(jì)算方法是為適合(fit)曲面而采用的坐標(biāo)點(diǎn)算法。 def aux1 loop n(1,nptab1) rr=xtable(1,n) zz=ytable(1,n) loop k(1,nprof

14、) case_of k 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 1 xtable(n+10,k)=rr ytable(n+10,k)=-40.0 case 2 xtable(n+10,k)=rr ytable(n+10,k)=0.0 end_case end_loop end_loop end aux1 當(dāng)執(zhí)行這個(gè)FISH函數(shù)后，生成(nptab1+10,nprof)組數(shù)據(jù)，這些數(shù)據(jù)為Aux2調(diào)用。注意：本段程序產(chǎn)生的數(shù)據(jù)是為那

15、個(gè)特定的曲面邊坡制作的，不是一個(gè)general purpose的LIB 5）FISH函數(shù)aux2的作用是使用WATER TABLE face命令產(chǎn)生水 面，把table 1的ytable值作為z坐標(biāo)，aux1產(chǎn)生的數(shù)據(jù)作為x,y坐標(biāo)。 def aux2 command set grav 0 0 -10 water dens 1000 end_command loop n(1,nptab1-1) loop k (1,nprof-1) xx1=xtable(n+10,k) yy1=ytable(n+10,k) zz1=ytable(1,n) xx2=xtable(n+10,k+1) yy2=yta

16、ble(n+10,k+1) zz2=zz1 xx3=xtable(n+11,k+1) yy3=ytable(n+11,k+1) zz3=ytable(1,n+1) xx4=xtable(n+11,k) yy4=ytable(n+11,k) zz4=zz3 command WATER TABLE face xx1,yy1,zz1 xx2,yy2,zz2 xx3,yy3,zz3 WATER TABLE face xx1,yy1,zz1 xx3,yy3,zz3 xx4,yy4,zz4 end_command end_loop end_loop end aux213、 FLAC結(jié)果數(shù)據(jù)的導(dǎo)出 pl s

17、et back white ;設(shè)置背景 set plot bitmap size (800,510) ;設(shè)置圖片大小 plot set caption size 37 ; 設(shè)置圖例數(shù)字大小 (default=35 范圍10-50) mainwin size 1 1 position 0 0.8 ; command window plot set window size 0.8 0.8 position 0.1 0.0 ; graphics window plot set center 100 0 40 根據(jù)情況更改14、一段fish命令的解釋 def aux1 定義函數(shù) loop n(1,np

18、tab1) 開始一個(gè)循環(huán)，從1循環(huán)到變量nptab1 rr=xtable(1,n) 從table 1 的x列中讀第n個(gè)數(shù)付給變量rr zz=ytable(1,n) 從table 1 的y列中讀第n個(gè)數(shù)付給變量zz loop k(1,nprof) 開始第二個(gè)循環(huán)，從1循環(huán)到變量nprof case_of k 下面是一個(gè)選擇結(jié)構(gòu)，利用變量k判斷 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 1 當(dāng)k1的時(shí)候執(zhí)行下面語句 xtable(n+

19、10,k)=rr table n+10的x列第k個(gè)值等于變量rr ytable(n+10,k)=-40.0 table n+10的y列第k個(gè)值等于變量-40 case 2 當(dāng)k2的時(shí)候執(zhí)行下面語句 xtable(n+10,k)=rr table n+10的x列第k個(gè)值等于變量rr ytable(n+10,k)=0.0 table n+10的y列第k個(gè)值等于變量0 end_case 結(jié)束選擇結(jié)構(gòu) end_loop 結(jié)束第二個(gè)循環(huán) end_loop 結(jié)束第一個(gè)循環(huán) end 函數(shù)定義結(jié)束 aux1 執(zhí)行什么定義的函數(shù) 15、有關(guān)命令的解釋 （1）apply nvel 0.1 plane dip 60

20、 dd 270 range left_boun apply 施加邊界條件，plane法向速度為0.1，plane 面的位置：傾角60（在xoy平面從-Z軸方向逆時(shí)針旋轉(zhuǎn)），dd 走向270度，從+y方向繞z軸順時(shí)針旋轉(zhuǎn)； （2） ini szz -5e4 grad 0 0 -1e4 ini 初始邊界條件，szz=-5e4+x*0+y*0+-z*-1e4,注：-5e4是range范圍起點(diǎn)坐標(biāo)的SZZ值，x,y,z所施加應(yīng)力網(wǎng)格點(diǎn)的坐標(biāo)。16、滲流示例 ;-參數(shù)部分 bulk modulus, K 390 MPa 體積模量（土質(zhì)比較硬） shear modulus, G 280 MPa 剪切模量

21、soil dry density, d 1200 kg/m3 土的密度 water density, w 1000 kg/m3 水的密度 wall density, wal 1500 kg/m3 檔土墻的密度 permeability, k 1012 m2/Pa-s 滲透系數(shù) porosity, n 0.3 孔隙率 fluid bulk modulus, Kf 2.0 GPa 流體的體積模量 ;-以下是命令流 ;- ; Excavation in a saturated soil ;- config fluid ;設(shè)置流體算法 ; - geometrical model - 建立模型 gen

22、zone brick p1 12 0 0 p2 0 12 0 p3 0 0 12 size 12 12 12 rat 1 1 1 group soil ;設(shè)置土組 group excavate range x 0 4 y 0 4 z 0 5 ;設(shè)置開挖范圍組 group wal1 range x 4 5 y 0 5 z 0 7 ; 設(shè)置檔土墻 1 組 group wal2 range x 0 4 y 4 5 z 0 7 ;設(shè)置檔土墻 2 組 group wall range group wal1 any group wal2 any ;檔土墻1 和2 合并為一個(gè)檔土墻組 ; - fluid f

23、low model - 流體模型 model fl_iso ;各項(xiàng)同性流體 prop perm 1e-12 poro 0.3 ;設(shè)置滲透系數(shù)和？比 or ？率 ini fdensity 1e3 ; 初始化水密度 ini fmod 2.0e9 ftens -1e-3 ; 初始化水的體積模量 及 ？ model fl_null range group excavate ; 開挖部分沒有水流 model fl_null range group wall ;檔土墻 不滲水 ini pp 0 grad 0 0 1e4 ;初始化孔隙水壓力 fix pp range z -0.1 0.1 ;約束 應(yīng)該是整個(gè)

24、模型的底部吧？ fix pp range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;開挖部分也給約束住了，而其他的地方不管 ; - mechanical model - 力學(xué)模型 model elas ; 彈塑性 prop bul 3.9e6 shea 2.8e6 ;設(shè)置體積模量和 剪切模量，這兩個(gè)必須，有一次我只設(shè)置了體積模量，就錯(cuò)的一塌糊涂 model null range group excavate ;力學(xué)模型的開挖 ini density 1.2e3 ;初始化密度 ，應(yīng)該是土的密度 ini density 1.5e3 range group wall ;初始化

25、墻（檔土墻）的密度， fix x range x -.1 .1 ;施加約束，據(jù)說差分法約束0 的時(shí)候，必須取一個(gè)-0.1<0<0.1,范圍 fix x range x 11.9 12.1 ;施加約束，或者叫做邊界條件 fix y range y -.1 .1 ;施加約束 fix y range y 11.9 12.1 fix z range z 11.9 12.1 ; initial total stresses 出世總應(yīng)力 ini szz 0 grad 0 0 -1.5e4 ;z方向的應(yīng)力大一些 ini sxx 0 grad 0 0 -1.2e4 ;水平方向的，小一些 ini s

26、yy 0 grad 0 0 -1.2e4 ; apply nstress 0 grad 0 0 -1.2e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0 apply nstress 0 grad 0 0 -1.2e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0 apply nstress -7.5e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1 ; - setting - 設(shè)置 set gravity 0 0 10 ;設(shè)置重力加速度 ; - initial state - ;初始化狀態(tài) solve force

27、1 ; check initial equilibrium 求解平衡 ; - histories - 記錄 set hist_rep 40 ;記錄步為40 hist fltime ;記錄時(shí)間 hist gp pp 0 0 6 ;記錄某點(diǎn)的孔隙水壓力 hist gp xdis 4 0 0 ;水平位移 hist gp xdis 4 0 2 ;記錄水平位移 hist gp xdis 4 2 0 hist gp zdis 0 0 5 ;記錄豎向位移 hist gp zdis 2 0 5 hist gp zdis 4 0 5 hist gp zdis 2 2 5 hist gp zdis 4 2 5 h

28、ist gp zdis 4 4 5 hist gp zdis 10 0 1 hist gp zdis 10 0 2 ; - excavation - 開挖 set fluid off ;關(guān)閉水 ; apply pore pressure at walls 在墻上施加孔隙水壓力 apply nstress 0 grad 0 0 -1.e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0 apply nstress 0 grad 0 0 -1.e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0 apply nstress -5.e4 range x

29、 0.0 4.0 y 0.0 4.0 z 4.9 5.1 solve ;force 1 ;求解 save exc1.sav ; 保存 ; ; - drainage - 排水 apply remove nstress ;撤掉剛才的壓力 def relaxsetup ;定義一個(gè)釋放 函數(shù)的參數(shù) step0 = step end relaxsetup ;調(diào)用 def relax ;定義一個(gè)釋放函數(shù) rstep = step - step0 if rstep < ncyc then relax=1.0-(float(rstep)/float(ncyc) else relax = 0.0 endi

30、f end ;結(jié)束 set ncyc = 1000 ;先賦值，隨后調(diào)用 apply nstress 0 grad 0 0 -1.e4 hist relax & range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0 apply nstress 0 grad 0 0 -1.e4 hist relax & range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0 apply nstress -5.e4 hist relax & range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1 cyc ncyc solve ;求解 save

31、 exc2.sav ;保存 ；; - percolation - 滲透 set fluid on ;設(shè)置水算法 fix pp 0 range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;施加孔隙水壓力 cyc 9000 ;循環(huán) save exc3.sav ;保存 plot create excav ;建立一個(gè)顯示的視圖 一下內(nèi)容都是顯示內(nèi)容 plot set rot 200 0 195 plot set cent 6 6 6 plot set dist 39.18 plot set magn 0.8 plot set plane ori 0 0 0 normal 0 -1

32、 0 plot add cont pp out on plot add block group lgra range group wall plot add cont pp int 10000 max 110000 range x 0 4 y 0 4 z 5 7 plot add flow plane plot show Ret17、 用什么命令可以位移圖中顯示最大位移點(diǎn)的坐標(biāo)？答：用fish編了一個(gè)，算法好像有點(diǎn)煩，不過目的可以實(shí)現(xiàn) 你可以參考一下。 def get_gp_maxdisp gp0_disp = gp_xdisp(gp_head)*gp_xdisp(gp_head) gp0_d

33、isp = gp0_disp + gp_ydisp(gp_head)*gp_ydisp(gp_head) gp0_disp = gp0_disp + gp_zdisp(gp_head)*gp_zdisp(gp_head) gp0_disp = sqrt(gp0_disp) p_gp=gp_head ;找最大值 loop while p_gp # null gp_disp = gp_xdisp(p_gp)*gp_xdisp(p_gp) gp_disp = gp_disp + gp_ydisp(p_gp)*gp_ydisp(p_gp) gp_disp = gp_disp + gp_zdisp(p_

34、gp)*gp_zdisp(p_gp) gp_disp = sqrt(gp_disp) if gp_disp>gp0_disp gp0_disp=gp_disp endif p_gp = gp_next(p_gp) endloop ;找最大值的坐標(biāo) p_gp=gp_head ss=0 loop while p_gp # null gp_disp = gp_xdisp(p_gp)*gp_xdisp(p_gp) gp_disp = gp_disp + gp_ydisp(p_gp)*gp_ydisp(p_gp) gp_disp = gp_disp + gp_zdisp(p_gp)*gp_zdis

35、p(p_gp) gp_disp = sqrt(gp_disp) if gp_disp=gp0_disp ss=ss+1 xtable(1,ss)=gp_xpos(p_gp) ytable(1,ss)=gp_ypos(p_gp) endif p_gp = gp_next(p_gp) endloop end get_gp_maxdisp 18、 如何顯示計(jì)算后模型的變形狀態(tài)？答：pl sk magfac *(一個(gè)放大的數(shù)字)19、如何查看剖面上任一點(diǎn)的位移值 答：print gp disp range x() y( ) z( )或 print gp disp range id plot gp di

36、sp range x() y( ) z( )或 plot gp disp range id20、group命令劃分斜面 答：group groupname colour range plane above dd * dip * origin * plane below dd * dip * origin * dd 和dip就是兩個(gè)斜面的傾向和傾角,origin是這個(gè)斜面通過的某個(gè)點(diǎn),這個(gè)命令就可以group兩個(gè)斜面間的區(qū)域 21、ansys導(dǎo)入flac的時(shí)候，如何劃分一下group 答：在ansys中對不同group用不同的材料號（mat）定義，再用ansys2flac之類的程序轉(zhuǎn)換成flac

37、的.dat 文件就可，group自動(dòng)生成。 22、 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 這個(gè)平面 是固定的 b 表示在x=0 這個(gè)平面 只x方向不能發(fā)生運(yùn)動(dòng)。 23、如何知道每個(gè)單元的ID？ 答：用鼠標(biāo)雙擊單元的表面，可以知道單元的ID和坐標(biāo)。 24、如何進(jìn)行切片？ 答：plo set plane ori (點(diǎn)坐標(biāo)) norm (法向矢量) plo con sz plane (顯示z方向應(yīng)力的切片) 這兩行命令，用來顯示指定面上的位移變化云圖（把第二行改改，還可以顯示應(yīng)力

38、云圖，等其他）；plot con syy y向應(yīng)力云圖plot con szz outline on z向應(yīng)力云圖加顯示網(wǎng)格plot con szz effective on z向有效應(yīng)力plot con szz inter 1e4 云圖顯示增量（后面的值越小，圖上顏色條帶就越細(xì)，）plot con szz min -10e4 顯示最小值，圖上小于這個(gè)數(shù)的值都用同一色表示plot con szz max 10e4 顯示最大值，圖上大于這個(gè)數(shù)的值都用同一色表示plot fip outside on 顯示力的箭頭 現(xiàn)在看邊界，如果是不動(dòng)面，用fix x range如果邊界帶速率，或位移用inise

39、t grav 0 0，-9.81ini dens 1000設(shè)置重力就一定要設(shè)置密度，不考慮重力就不要設(shè)計(jì)密度。兩者共存的，設(shè)計(jì)程序的時(shí)候一定要記的。set mech force 50solve用SET命令設(shè)定了最大不平衡力的上限值為50N，即：計(jì)算的最大不平衡力一旦小于此值，我們便人為模型達(dá)到了平衡狀態(tài)。求解過程終止。設(shè)置這個(gè)目的是在一定精度要求下，加快計(jì)算（值越大，算的量越少，花的時(shí)間也越少）。25、如何在程序中進(jìn)行暫停，并可恢復(fù)計(jì)算？ 答：在命令中加入pause命令，用continue進(jìn)行繼續(xù)。在我們分步求解中想得到某一個(gè)過程中的結(jié)果，不用等到全求完，還可以在分布求解錯(cuò)誤的時(shí)候就進(jìn)行改正，

40、而不是等到結(jié)果出來。26、如何跳過某個(gè)計(jì)算步？ 答：在計(jì)算中按空格鍵跳過本次計(jì)算，自動(dòng)進(jìn)入下一步 27、Fish是什么東西？Fish是否一定要學(xué)？答：是FLAC3D的內(nèi)置語言，可以用來進(jìn)行參數(shù)化模型、完成命令本身不能進(jìn)行的功能。Fish可以不用學(xué)，需要的時(shí)候查Mannual獲得需要的變量就可以了。 28、如何調(diào)用一些可選模塊？ 答：config dyn (fluid, creep, cppudm) 29、range的命令： 取圓球面的范圍，就是同一個(gè)中心，不同半徑之間夾的那一部分，用在給球面施加力，約束時(shí)用該命令，命令：range annulus center xyz（中心坐標(biāo)）radius

41、r1 r2（兩半徑），想要表達(dá)的是球面時(shí)，r1r2取值相差小點(diǎn)即可，如radius 2.9 3.1 這時(shí)范圍就是半徑為3的面。圓柱的環(huán)面施加力或約束時(shí)也用到該命令。range cylinder end1 xyz end2 xyz radius r(end是表示柱體斷面的中心點(diǎn)坐標(biāo))range plane above 表示的是z大于0的上部結(jié)構(gòu)range plane blow 表示的是z小于0的下部結(jié)構(gòu)range dd a dip b origin xyz 表示的是走向角度a，傾向角度b的面dd表示，當(dāng)面與xoy面相交時(shí)的垂線與y軸正向的夾角，順為正，說白話，就是面的指向方向面是朝那個(gè)方向的意思

42、。dip面與xoy面的夾角，與z軸負(fù)向?yàn)檎琽rigin表示面經(jīng)過xyz點(diǎn)；range normal xyz 法線方向?yàn)閤yz的面range not 沒試用過（范圍外的東西，想想用到圓上取范圍非常好）。 30、如何調(diào)用FISH文件？ 答：新手調(diào)用FISH文件出錯(cuò)的原因：是：文件名為 fishcall.fis.txt(因?yàn)樵谟浭卤纠镞@個(gè)txt的后綴是默認(rèn)的）。修改方法是：將 fishcall.fis.txt里面的內(nèi)容全部復(fù)制到新的記事本里，再按保存，保存時(shí)候?qū)⒈４骖愋透臑樗形募?*,*),這時(shí)你便可在文件名欄內(nèi)輸入 fishcall.fis，類型也就正確了。 31、問a式與b式有何區(qū)別？答：f

43、ix x y z range x -0.1 0.1 .a fix x range x -0.1 0.1.b a 表示 在x=0 這個(gè)平面 是固定的； b 表示在x=0 這個(gè)平面 只x方向不能發(fā)生運(yùn)動(dòng)。32、如何查看各個(gè)時(shí)段不平衡力的具體數(shù)值？ 答：采用his來記錄計(jì)算，包括位移應(yīng)力等 命令his unbal；his gp(zone) zdis range (0 0 0) 或者id=？； 導(dǎo)出數(shù)據(jù)命令his write n vs m begin 時(shí)步 end 時(shí)步+ file filename.his n表示紀(jì)錄的id m表示時(shí)步； 要導(dǎo)出不平衡力的具體數(shù)值 his unbal step 100

44、000 or solve his write 1 vs step begin 1 end 1000 file 123.his 使用上述命令就可以查看各個(gè)時(shí)步下的不平衡力的具體數(shù)值。33、如何顯示變形輪廓線的命令？ plo ske magf 10 其中10為放大系數(shù)34、怎樣查看模型的矢量圖？pl grid vel (disp) redplo blo gro可以看到不同的group的顏色分布（看模型的樣子）；plo gpfix red sk 看模型的邊界情況plo fap red sk 看模型的體力分布plo con dis (xdis, ydis, zdis) 看模型的云圖：位移 plo co

45、n sz (sy, sx, sxy, syz, sxz) 看模型的云圖：應(yīng)力plo dis (xdis, ydis, zdis) 看模型的矢量圖 plo info 看模型有多少單元、節(jié)點(diǎn)plot grid 可以查看網(wǎng)格，plot grid num 可以查看節(jié)點(diǎn)號。plot fix 查看邊界條件plot grid mark 查看預(yù)先留設(shè)的開挖邊界pl syy (sxx xdis ydis sig1.) fill；查看模型的云圖。注意，帶fill關(guān)鍵字就是查看云圖，不帶fill 關(guān)鍵字就是查看等直線圖。plo add ske 剖面的沉降云圖上增加網(wǎng)格plo add dis plane 剖面上的矢量

46、plo add axe 剖面上坐標(biāo)值35、怎樣查看INTERFACE？答：pl grid iface id redprint iface id ; id 是interface的編號36、FLAC怎么調(diào)整視圖？答：ctrl+z 或shift+z ,然后點(diǎn)擊要放大部位的左上角和右下角就可以放大某個(gè)部位，如果要縮小，雙擊屏幕。37、判斷模型中網(wǎng)格點(diǎn)是一個(gè)還是兩個(gè)的方法。 答：顯示模型網(wǎng)格點(diǎn)的ID，如果有兩個(gè)ID號，說明是兩個(gè)。節(jié)理面上的網(wǎng)格點(diǎn)應(yīng)該有兩個(gè)ID號。38、flac中影響運(yùn)行時(shí)間的因素答：（1）單元材料，結(jié)構(gòu)面之間的剛度對比很大；（2）劃分單元尺寸對比很大，盡量避免長細(xì)比大于5：1

47、的細(xì)長單元。39、簡單命令解釋：答：1.plot add cont disp out on 時(shí)時(shí)顯示云圖2.plot block group gpnum on 顯示節(jié)點(diǎn)編號3.plot inter red ske         看接觸面的效果.4.del range group 1 not 不刪除group信息  5.del range group 2 刪除group信息6.ini z add * range group * （name） 7.interface 1 face range x（

48、*，*）y（*，*）z（*，*）或者 interface 1 face range plane orig(x y z 坐標(biāo)）norm（法線方向） 建立接觸面8.water table clear 清除水位線9.pl set back white 設(shè)置背景10.set plot bitmap size (800,510) 設(shè)置圖片大小11.plot set caption size 37 設(shè)置圖例數(shù)字大小12.plot set center 100 0 40 設(shè)置圖片中心13.fix x y z range x -0.1 0.1 表示在x=0 這個(gè)平面是固定的14.fix x range x -

49、0.1 0.1 表示這個(gè)平面只x方向不能發(fā)生運(yùn)動(dòng)40、 在固定面上施加位移時(shí)，fix與apply不能同時(shí)使用，應(yīng)該用ini和fix,才能施加上。41、在進(jìn)行三軸壓縮時(shí)候，如何施加圍壓？答：使用apply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius 1，其中X和Z軸是在我們眼睛正對的面。42、在進(jìn)行三軸壓縮時(shí)候，如何對軸向力進(jìn)行位移控制？答：使用 ini yvel 1e-7 range y -0.1 0.1 ini yvel 1e-7 range y 1.9 2.1 即可，數(shù)值可以變化。43、在FLAC3D計(jì)算中，如何記錄程

50、序中的值，如果不想時(shí)時(shí)記錄該怎么弄？答：用 hist wtite n1 vs n2 file.txt， 想間隔記錄可用set hist_rep n44、想記錄接觸面上的相關(guān)力學(xué)參數(shù)該怎么辦？答：hist inter i nstress /ndis/sstress/sdisp45、接觸面設(shè)置參數(shù)包括哪幾個(gè)？答：coh dial fric kn ks ten.46、如何輸出單元格的數(shù)據(jù)？答；expgrid .flacsd47、建模的過程中，或同一個(gè)文件中有重復(fù)應(yīng)用的命令，就可以用宏命令了。答：宏命令： macro 名稱 命令的內(nèi)容例：macro shiping plot con zdip outl

51、ine on用shiping 代替后面的引號里的東西。在建模中只要有相同的都可以用shiping代替，縮減了書寫命令的時(shí)間，整體上也能讓人看的更明白。注意，你命名的符號最好跟falc自帶的命令不要重復(fù)，處理方法是，名稱后帶數(shù)字。如shiping1；如果你想修改宏對象，改變命令的內(nèi)容，在名稱上加單引號就是替換的意思。 例macro shiping1 新的命令。48、group的用法：group直接加到模型的后面相當(dāng)于給整個(gè)模型命名；group name range x xyy zz給一個(gè)范圍的體命名 。例：在gen zone brick - size - group name49、plot中ad

52、d與set的區(qū)別add原來有的加新的用addset原來沒有的，設(shè)置用setplot create GravV plot set plane dip=90 dd=0 origin=3,4,0 plot set rot 15 0 20 plot set center 2.5 4.2 4.0 plot add bound behind plot add bcont szz plane plot add axes plot show 。我們想顯示別的圖像，又不想把這個(gè)弄好的圖像丟掉，則用plot create name建立一個(gè)圖像，在建立新的圖像時(shí)，plot create name1（新的名）顯示以前

53、的圖像，那么輸入如下命令：plot current name。另一種方法，輸入：plot show name 50、不同的區(qū)域定義不同的本構(gòu)模型model mohr range z=0,5 prop ClayeyGravel range z=0,5 在定義參數(shù)后面加個(gè)range范圍，為了避免在輸入命令時(shí)總要通過range全局坐標(biāo)域的方法來限定該命令執(zhí)行的區(qū)域，我們可以先通過range命令來給不同的區(qū)域命名，然后就可以用range區(qū)域名的方法來限定范圍了。這樣就避免了每次輸入命令要限定在該區(qū)域執(zhí)行時(shí)，都要輸入一遍坐標(biāo)范圍的麻煩。具體如何給區(qū)域命名，以及日后如何用區(qū)域名來限定該區(qū)域，可以看下面這個(gè)例子： range name trench x=-1,1 y=0,4 z=-2,2 model null range trench51、右手法則：建立模型時(shí)，各點(diǎn)的排放時(shí)有規(guī)律的。方法如下：拿出右手伸直了4指合并指向就是p0到p1的方向，掌心指向就是p0到p2的方向，拇指與4指垂直的方向就是p0到p3的方向。很多人說網(wǎng)格后的云