FLAC3D命令的FAQ

1、FLAC3D命令的FAQ!1. FLAC3D是有限元程序嗎？' q' C8 u2 X% / 2 O$ J答：不是！是有限差分法。* ) V' ?- G( C( J8 C+ j1 a: b+ X) , _: i2 _2. 最先需要掌握的命令有哪些？; O% g4 t( q1 A8 x5 D答：需要掌握gen, ini, app, plo, solve等建模、初始條件、邊界條件、后處理和求解的命令。" k+ Y1 E* x: ( z0 o- ) S: o7 q+ f: 8 c- 3. 怎樣看模型的樣子？& ?& T' . m/ Z; |8

2、M* G答：plo blo gro可以看到不同的group的顏色分布0 |: B) P9 x. X. y' s0 k# L/ ! g! v+ m# c* _* v, S& g4. 怎樣看模型的邊界情況？. E( a2 3 k8 答：plo gpfix red sk* o& w4 E% s7 r- N( t% 8 4 N2 U* a  0 H7 , y1 S+ r! 5. 怎樣看模型的體力分布？2 T+ O# N6 B: n9 _- D0 d2 c答：plo fap red sk1 G& |& c5 % ?0 d& f8 D2

3、9 Z+ f, Y/ / w' F# O. H6. 怎樣看模型的云圖？- S. y- " r( 1 V$ U" L& B7 c答：位移：plo con dis (xdis, ydis, zdis)/ % 1 ! ' O1 B1 t* J% m  Q    應(yīng)力：plo con sz (sy, sx, sxy, syz, sxz)8 M, U( n( X  e5 Z9 v; 4 H+ J    3 ( l/ r" m% j7. 怎樣看模型的矢量圖？' B)

4、 1 J/ A$ W1 X答：plo dis (xdis, ydis, zdis)9 D) c" R4 m* P4 Z' G: q3 _7 K1 t- c' U8 G% f/ 6 D  v) r8. 怎樣看模型有多少單元、節(jié)點(diǎn)？! q5 Q5 J- U: Z: b! E2 s% r. G8 答：plo info+ r! " k3 i9 x$ b( a5 Z"   U3 j: H1 u7 V* c1 & # x9. 怎樣輸出模型的后處理圖？% f, y6 g- U7 c9 z# |) B! $ r答：F

5、ile/Print type/Jpg file，然后選擇File/Print，將保存格式選擇為jpe文件6 6 " Y/ M+ r+ S) F% K; k: I3 e  P+ C3 Y. ?/ N10. 怎樣調(diào)用一個(gè)文件？$ e: f9 I5 H% Y* i: O答：File/call或者call命令: t/ d9 G7 i7 |& j6 b- w. s3 x7 |2 / D9 D  m. d7 " y10. 如何施加面力？) y4 |1 J  N8 5 r答：app nstress. N/ t2 V, C

6、" m+ q9 m' 2 P1 u: l9 K( s* W+ 11. 如何調(diào)整視圖的大小、角度？# n. t2 L  - B# G9 |  W答：綜合使用x, y, z, m, Shift鍵，配合使用Ctrl+R，Ctrl+Z等快捷鍵4 R3 0 |( P0 a& Z: D; m9 U4 H/ K4 g9 L  i6 q# q: m12. 如何進(jìn)行邊界約束？. H3 p% F4 Z, a/ B答：fix x ran （約束的是速度，在初始情況下約束等效于位移約束）" a2 i; n6 0 q- r%

7、 I7 r8 g1 w( V0 _5 e  M! Z1 d13. 如何知道每個(gè)單元的ID？% L) I! |2 h" s! G' N! t* |答：用鼠標(biāo)雙擊單元的表面，可以知道單元的ID和坐標(biāo)/ x/ & k7 E0 z6 H  y5 : & Q7 g" ; b: Z' j3 r$ 14. 如何進(jìn)行切片？# " l# t8 T( ; L* n答：plo set plane ori (點(diǎn)坐標(biāo)) norm (法向矢量)- t$ E( K" / |- l+ v% X0 t  &

8、#160; plo con sz plane (顯示z方向應(yīng)力的切片)+ s- ! 2 : . Y0 C0 y*     3 H3 C  C5 M/ p15. 如何保存計(jì)算結(jié)果？+ s: L+ S6 S2 Q: W* 答：save +文件名$ ( ( G% v' O* G( p: I$ n( & _' I0 C; |1 X4 B$ l- h) q# 9 P16. 如何調(diào)用已保存的結(jié)果？3 d( p3 U7 f2 v& V! K5 p( U答：rest +文件名；或者File / Restor* s2 : ?/ t, 7 W

9、1 G' x- l) b" W& p& q17. 如何暫停計(jì)算？; V% 0 P7 a9 a6 o0 A- s! b" / q' p, D答：Esc3 C4 ) f1 X; v. o+ L/ D! r4 n! 0 x8 r18. 如何在程序中進(jìn)行暫停，并可恢復(fù)計(jì)算？" N% _( ! & j9 X& a$ O  O答：在命令中加入pause命令，用continue進(jìn)行繼續(xù)! v/ D2 1 J. e9 R8 Q/ G8 Z( C0 F* 0 _& H# $ I* 0 19. 如何跳過某個(gè)計(jì)

10、算步？4 K4 / ' W/ # s答：在計(jì)算中按空格鍵跳過本次計(jì)算，自動(dòng)進(jìn)入下一步3 n1 P8 L( F7 V7 j# 9 H0 7 8 m/ f# Z! Z20. Fish是什么東西？" H  C1 L* i0 u7 U, a% o' m; u答：是FLAC3D的內(nèi)置語言，可以用來進(jìn)行參數(shù)化模型、完成命令本身不能進(jìn)行的功能& O% T  k+ K/ M3 T$ a$ _' , C+ y# D! D  _7 ?. H3 N21. Fish是否一定要學(xué)？, . n. |  

11、;j8 x5 K2 ) ' R2 k答：可以不用，需要的時(shí)候查Mannual獲得需要的變量就可以了9 M( f4 v& |4 x4 n( C0 q% o2 O/ z# F2 1 M22. FLAC3D允許的命令文件格式有哪些？' f# X2 C! q( l4 V答：無所謂，只要是文本文件，什么后綴都可以!   o6 1 N3 o& h6 f& * + S+ c# s* P( G! y( T* a, r1 z23. 如何調(diào)用一些可選模塊？+ q& L  0 h. a2 g% o: e7 T# R答：confi

12、g dyn (fluid, creep, cppudm)- s0 ; F1 Y/ X# s(F24 .如何在圓柱體四周如何施加約束條件？, W0 h  _$ l, H可以用fix . ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2<r1，其實(shí)就是選擇range的時(shí)候選兩個(gè)圓柱的差，即得到邊界。. N9 g$ . |. E! r5 L3 M) e3 Q3 0 y6 K% G5 A8 ?25.   如何能把一個(gè)PLOT的圖像數(shù)據(jù)導(dǎo)出來以便用其他軟件繪圖？

13、60; x/ 5 v/ t2 x8 R: 1 E回答：用set log on 命令，把數(shù)據(jù)導(dǎo)出來，轉(zhuǎn)到excel里處理一下，然后用surfer或者什么作圖軟件繪制就行了。, q$ N0 D8 Y5 j7 o" M4 L% W/ % g# K% f# ) k) Y: V1 S26.關(guān)于dd、dip的理解：+ V0 u0 C# H0 Z, V. U3 U: m& Q& pdip-direction angle, measured in the global xy-plane clockwise from the positive y-axis 這個(gè)是從y軸到x順

14、時(shí)針方向所轉(zhuǎn)過的角度；dip dip angle, measured in the negative z-direction fromthe global xy-plane 這個(gè)是xy平面與z軸負(fù)方向的夾角。8 K* A7 S2 p+ m' a* / d( k# cdip表示對(duì)稱參照面與xy平面的夾角，對(duì)稱參照面與xy平面的夾角在xy平面的投影是一條射線 ，dd表示y軸正向順時(shí)針到那條射線的夾角。首先應(yīng)該按照dd的方向大體確定這個(gè)面的朝向，dd 指的是從y軸正方向按順時(shí)針（clockwise）方向轉(zhuǎn)向所要確定面的法線方向在xy平面上的投影的夾角，然后再確定dip，dip指的是從xy平面

15、轉(zhuǎn)向所要確定的平面的角8 " b: c) L5 p/ q5 z, - G* m+ E: x% t0 9 W# a" 0 I27.  如何調(diào)用FISH文件？( E6 h8 f% $ W# F9 ! 新手調(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，類型也就正確了。& |; o5 N/ L* O

16、# e$ G+ O  T5 T2 B6 Z28.  如何理解Tension-p shear-p tension-n shear-n表示的含義？1 C- . N6 P* z( N# C2 x5 在計(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)./ O9 g* p. P. 8 |- U* F4 h6 t2 M5 L. s# w4 L3 ) 7 O9 s29.

17、0;  繪制相關(guān)曲線命令的一些解釋9 l- B% N$ c$ l5 A$ E1) 在plot hist m vs n的形式里,m代表y軸,n代表x軸(不管m,n的正負(fù));2 M: 5 _2 G' c. I. V7 S2) "-"表示對(duì)其值作"mirror",如果大家使用AUTOCAD的話就明白"mirror"的含義了,比如上面的 hist 1 gp ydisp的值是0100,那么vs -1就變成-1000.以此類推.; A: h0 t* c" B& _& S9 & e1 g) I3)

18、 為什么要這么做呢?主要是為了符合巖土工程的習(xí)慣需要,例如樁載荷試驗(yàn)曲線,沉降曲線等都不是畫在第一象限內(nèi). 2 R$ C! a5 w( Q$ N- L8 R) s- c: Y  c: |- b9 L# f4 r% e30.   如何查看剖面上任一點(diǎn)的位移值？4 x: W, V8 Y4 D# Eprint gp disp range x() y( ) z( )或) y0 I$ v% K%   cprint gp disp range id + 5 B5 H$ w  t) s6 gplot gp disp range

19、x() y( ) z( )或% 3 Q3 k4 Q) d: V& M7 a2 ( Bplot gp disp range id! ) J& X( t; Q0 z4 i, o: u. 5 J: u1 8 E' C5 b) Y$ K$ c% A( a7 m8 31.    fix x y z range x -0.1 0.1 .a2 i; 1 R& . k4 P3 W$ L. N         fix x range x -0.1 0.1.b' Z. F* : P!  

20、; z        請(qǐng)問a式與b式有何區(qū)別？0 V6 V. E) H# k3 a9 ma 表示 在x=0 這個(gè)平面 是固定的# q  d  N4 z! c- Q. Q7 7 Yb 表示在x=0 這個(gè)平面 只x方向不能發(fā)生運(yùn)動(dòng)( A8 t/ J. g2 f0 % q" I) y# z& T! m' I6 32.    initial 與 apply 有何區(qū)別？# M+ a1 d( C( V* M* M) x+ 9 Vinitial初始化命令

21、，如初始化計(jì)算體的應(yīng)力狀態(tài)等；' i% B5 R: v2 e/ q& l' k* japply邊界條件限制命令，如施加邊界的力、位移等約束等。3 D( z' S7 b& H- A( l# ?) Ainitial的應(yīng)力狀態(tài)會(huì)隨計(jì)算過程的發(fā)生而發(fā)生改變，一般體力需要初始化，而apply施加的邊界條件不會(huì)發(fā)生變化。* 0 v  h7 F# g2 T; 3 d5 c5 D' U6 B5 S1 5 g. Y33.   如何查看各個(gè)時(shí)段不平衡力的具體數(shù)值？, N, K+ 5 f. M4 U% s采用his來記錄計(jì)算，包

22、括位移應(yīng)力等 命令his unbal ) V0 t* V% u$ V# N5 M2 # L' ehis gp(zone) zdis range (0 0 0) 或者id=？4 U" g2 r% v- U2 e  _2 I6 ' _0 W: r3 K導(dǎo)出數(shù)據(jù)命令his write n vs m begin 時(shí)步  end  時(shí)步 file filename.his1 a( V% h/ ?) # w3 B; B! & a- cn表示紀(jì)錄的id m表示時(shí)步9 x6 c7 o9 Z. A8 V% G0 , h

23、. N要導(dǎo)出不平衡力的具體數(shù)值' C; u' $ Y7 C2 E& ' Rhis unbal3 u1 p" ) t% C- a" Xstep 100000 or solve/ m4 f6 b: b+ ; his write 1 vs step begin 1 end 1000 file 123.his1 - V2 x1 k8 G; Q使用上述命令就可以查看各個(gè)時(shí)步下的不平衡力的具體數(shù)值 ; o# v0 f; h$ ! t6 t" t9 j) O1 z2 O34.    如何顯示變形輪廓線的命令？7 v3 s6 W

24、2 Z( U2 C" P3 s$ plo ske magf 10   其中10為放大系數(shù)  FLAC2d命令的FAQ：' C) ?+ k: / y1 n2 H/ 1、怎樣開始入門？( I3 y: O. r% 4 z+ O0 X答：初步學(xué)習(xí)首先看說明書 User's gurid的Section2 Getting started和Section 3 Problem sloving with FLAC2 6 F) k  T/ U2、哪些命令需要掌握？! X; x  o  a+

25、 x) U答：grid gen prop interface 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 3DEC PFC PFC3D）的ra

26、nge大不相同 ，要注意區(qū)分。2 p- h  j6 M/ |. w3、FLAC和FLAC3D建模有什么區(qū)別？! Z) R( s9 + k' _8 D' ?. 答：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可以用很多次。/ |*

27、J) & ?7 2 Y+ j' T' d" A4、怎樣查看模型？$ C8 c. d; o2 p& e$ $ : q: M答：plot grid 可以查看網(wǎng)格，plot grid num 可以查看節(jié)點(diǎn)號(hào)。" b- Q1 U4 6 H$ F: s8 m; t/ 5、怎樣查看邊界條件？. H8 I9 U' Z( C' W+ O答：plot fix 1 x9 n7 C' - D; G- ?' 8 w6、怎樣查看預(yù)先留設(shè)的開挖邊界？/ k0 ! $ U% ?0 o5 Y* B答：plot grid mark。 

28、 U5 q8 o$ p! P. C! G/ U% L7、怎么查看模型的云圖？4 3 U9 w1 x7 s. ) R* f答：pl syy (sxx xdis ydis sig1.) fill ，注意，帶fill關(guān)鍵字就是查看云圖，不帶fill 關(guān)鍵字就是查看等直線圖。2 1 b* v5 L: b5 8 ' p8、怎樣查看模型的矢量圖？0 j/ q' 9 U: a- t, b* K! 答：pl grid vel (disp) red' V8 Q6 Q0 k) S! s* p9、怎樣查看INTERFACE？" m% m* G  E&#

29、160; Y( _答：pl grid iface id red & Z/ g, 6 A, w8 T5 B, s% t+ J      print iface id ; id 是interface的編號(hào)/ |2 C9 z2 U; d8 N* a& M# U; R10、FLAC怎么調(diào)整視圖？5 R( M*   z, t+ u  n- q答：ctrl+z 或shift+z ,然后點(diǎn)擊要放大部位的左上角和右下角就可以放大某個(gè)部位，如果要縮小，雙擊屏幕。1、判斷模型中網(wǎng)格點(diǎn)是一個(gè)還是兩個(gè)的方法&qu

30、ot; k) W% o* Q- c: a: W$ ?: B   顯示模型網(wǎng)格點(diǎn)的ID，如果有兩個(gè)ID號(hào)，說明是兩個(gè)。節(jié)理面上的網(wǎng)格點(diǎn)應(yīng)該有兩個(gè)ID號(hào)。% G( e: x8 r% h! m2、改變模型應(yīng)力邊界時(shí)FLAC與3DEC的不同之處5 N: z# A7 h& K' p)    在3DEC中，當(dāng)改變應(yīng)力邊界時(shí)，模型表面的應(yīng)力直接被更換，F(xiàn)LAC中改變應(yīng)力邊界時(shí)，模型4 Y  g8 P1 A" k! D表面的應(yīng)力為兩者應(yīng)力相加。4 3 Z8 k) T3 u) |0 F3、對(duì)DD,DIP的理解# g+ W+ Q

31、- R& T, 1 q    DD表示節(jié)理面正方向在XY面上的投影與Y軸的夾角，也等于節(jié)理面與XY面的交線與X軸的夾角；. " - t6 o. H    DIP表示節(jié)理面與XY面的夾角。; + x" + v: p. _" p4 塑性狀態(tài)的解釋# o) T1 M7 N( i7 , Z% g  -p  表示某一區(qū)域的應(yīng)力在模型運(yùn)行過程中進(jìn)入過屈服狀態(tài)，但現(xiàn)在已經(jīng)推出了屈服狀態(tài)；' X4 a, + V& w3 z$ L  n  

32、表示某一區(qū)域正處于屈服狀態(tài)%   r3 d" x2 i* z+ t) s6 U5 對(duì)網(wǎng)格點(diǎn)速度的理解2 t4 9 H$ N; . B: X# b  flac中網(wǎng)格點(diǎn)的速度表示的是單位時(shí)步網(wǎng)格點(diǎn)發(fā)生的位移，查看網(wǎng)格點(diǎn)的速度有兩種方法：% D+ z! N2 M3 Z( n' （1）繪制整個(gè)區(qū)域的速度圖' H* 4 9 J6 Z（2）記錄某一網(wǎng)格點(diǎn)的速度  X+ n" I# H& p6 e7 N6 flac中影響運(yùn)行時(shí)間的因素2 i) , h' o+ q  I; （1）

33、單元材料，結(jié)構(gòu)面之間的剛度對(duì)比很大；1 o: n* y% u+ ?3 b. N （2）劃分單元尺寸對(duì)比很大，盡量避免長(zhǎng)細(xì)比大于5：1的細(xì)長(zhǎng)單元。1 F# o: h& e, G$ + m5 n( x; N9 U $ . k% C3 W; J2 w% P4 M7 阻尼的選擇:   Y( # A9 : a# 7 o  （1）對(duì)于靜力分析，局部阻尼比較有效，是默認(rèn)狀態(tài)的(SET MECHANIC DANM LOCAL)。9 |7 J. B. s+ m1 ?  （2）組合阻尼，比局部阻尼更容易在大量網(wǎng)格點(diǎn)發(fā)生剛體運(yùn)動(dòng)時(shí)更好的收斂到穩(wěn)定狀

34、態(tài)，使用SET MECH+ / E*   g" N DAMP COMBINE來調(diào)用組合阻尼，組合阻尼對(duì)于減小動(dòng)能方面不如局部阻尼有效。! y- E+ t$ j) s3 i4 Z7 K+ 8 FLAC中有多少中材料模型& L, l. v( % L8 h$   內(nèi)置11中本構(gòu)模型$ O& _  X; s. G$ b7 T! O8 f: D% J  （1）零模型    用于開挖回填的區(qū)域% U$ d1 k. 3 X7 k4 J  （2）各向同性彈性模型9 X2 y; + X/ S! h8 y  （3）正交各向異性彈性模型) 9 v. w: x, 7 d5 J  （4）橫觀各向同性彈性模型4 I0 ?" F4 W/ P9 B  （5）德魯克普拉格塑性模型4 m: % U; i: 2 E: 4 w0 o, K" w  （6）莫爾庫侖塑性模型        一般巖土體的力學(xué)行為6 R# ?,