FLAC3D与FLAC常见命令的整理.doc

FLAC/FLAC3D常规问题的整理1FLAC3D命令的FAQ/ ?9 Y0 y7 B1 x0 _( h+ E4 lakewater整理) , a0 k( t2 I Q: |0 s看到其它板块上都有这个FAQ，也就是常见问题问答，今天抽了时间进行了整理，想到了就写下来了，因为看到很多初学者费了很多的时间，但是还是没有将常用的命令掌握，所以这个也可以作为入门的初级教材，使大家能够快速的上手，而不用为了某个小命令到处求助。1. FLAC3D是有限元程序吗？1 u# v* u k& 4 b. K% s答：不是！是有限差分法。2. 最先需要掌握的命令有哪些？/ B& g Y6 V. U# 8 q Y8 w2 r9 b答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。3 5 A5 y$ Z2 N5 R6 E$ i* L4 _- g8 X- O, _3. 怎样看模型的样子？0 D% s! R9 ! E2 QY6 e) E; Z& A答：plo blo gro可以看到不同的group的颜色分布0 g) b7 M. T6 H l; C4. 怎样看模型的边界情况？# o- d( j8 , v/ , G答：plo gpfix red sk. x! q8 e7 c, E3 m* u8 P& u+ C0 . g7 j8 ?( i5. 怎样看模型的体力分布？. g# W! 3 F; 答：plo fap red sk- A& p/ U% 7 ) b6 l4 m* h- W3 b- Y) q8 ) f6. 怎样看模型的云图？2 R) G# Z$ K8 o* S答：位移：plo con dis (xdis, ydis, zdis)4 c2 j- l( u- B; n( 应力：plo con sz (sy, sx, sxy, syz, sxz)8 a4 ! F3 Z8 d* F: i4 ?! + eD7. 怎样看模型的矢量图？. g$ b) G% ON2 S8 答：plo dis (xdis, ydis, zdis)5 P/ E* _6 J% Z6 z/ l# u4 y8. 怎样看模型有多少单元、节点？答：plo info: m19. 怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件1 4 s& 9 c* O1 e6 j z; / k9 4 J, U3 e0 e D Y2 x10. 怎样调用一个文件？2 H, 0 ?9 Q3 k) Z) m2 Y答：File/call或者call命令) ZMJd% N N 8 s- V- N6 R, t9 N/ 10. 如何施加面力？. r2 H4 b6 ?4 R0 m答：app nstress& T4 Y+ 4 q/ s# V! S5 j11. 如何调整视图的大小、角度？0 I$ T$ T1 l! D) 3 z* A答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键5 t) % E/ 6 U8 J/ N! U7 b12. 如何进行边界约束？, z: T/ N. B& 1 a答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束）( # # . X) % d# J/ t, D; b3 Y0 ! N7 Q1 B13. 如何知道每个单元的ID？+ o/ j& T! R d% d: h答：用鼠标双击单元的表面，可以知道单元的ID和坐标! B* Y1 6 y7 U2 . k K2 G4 z14. 如何进行切片？ 5 6 T1 d6 B/ J7 p答：plo set plane ori (点坐标) norm (法向矢量)7 7 V! j6 + y . L plo con sz plane (显示z方向应力的切片)5 g0 b ?4 E% k8 * F15. 如何保存计算结果？$ r$ e( H+ k! ( o; D% W7 S答：save +文件名& C, w+ K$ l/ C; _! R+ i z8 N, ; i& Y16. 如何调用已保存的结果？9 Q9 0 3 E u( o! m2 v- _% B答：rest +文件名；或者File / Restor/ B$ 1 I: Z5 , i: O5 h7 A5 I5 H; F& _4 g4 IN17. 如何暂停计算？2 U) F( B- D答：Esc0 ( a9 y% ! v$ d7 R h18. 如何在程序中进行暂停，并可恢复计算？1 H/ h) d: A1 M5 + g; H7 答：在命令中加入pause命令，用continue进行继续1 |: p: K6 Y8 / n2 s/ |19. 如何跳过某个计算步？) _# Z4 E5 i4 ) j; f% u3 v答：在计算中按空格键跳过本次计算，自动进入下一步; t, t, Q2 V5 D! s- b* m a# A0 R20. Fish是什么东西？6 X+ w1 ?% _0 A7 y e答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能8 u5 . d2 NV, R. i6 G4 G4 b W2 I2 B: q6 W1 x! H6 JG% A: L21. Fish是否一定要学？$ l1 R|. K4 B答：可以不用，需要的时候查Mannual获得需要的变量就可以了 u& q! Y: p) n0 D/ A w0 X3 n( R# j5 e7 t5 X- P22. FLAC3D允许的命令文件格式有哪些？& c. P$ Z1 ?Z8 y- A0 _答：无所谓，只要是文本文件，什么后缀都可以 T: Z# K8 k, 7 V0 g* g/ D* P1 C- G( e F$ 23. 如何调用一些可选模块？; N5 ?( o5 B+ f8 p答：config dyn (fluid, creep, cppudm) ( l0 Y1 q, c: O# K$ f) U3 c. z: u5 u后注：这个工作很繁琐，需要的时间很多，希望广大网友能够将自己曾经遇到的常见问题在后续跟贴，也为了将这个FAQ进行很好的充实。24 .如何在圆柱体四周如何施加约束条件？G( $ x( k% i可以用fix . ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2r1，其实就是选择range的时候选两个圆柱的差，即得到边界。2 |2 T8 lm6 _ G- r) j2 p25. 如何能把一个PLOT的图像数据导出来以便用其他软件绘图？. r m6 N: R+ T% # R回答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。 n, E 0 m& B j26.关于dd、dip的理解：( G: ! 3 NTdip-direction angle, measured in the global xy-plane clockwise from the positive y-axis 这个是从y轴到x顺时针方向所转过的角度；dip dip angle, measured in the negative z-direction fromthe global xy-plane 这个是xy平面与z轴负方向的夹角。3 w7 U# C* , _2 h+ vq9 Pdip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线 ，dd表示y轴正向顺时针到那条射线的夹角。首先应该按照dd的方向大体确定这个面的朝向，dd 指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角+ B( M( r# c! e. t: D |27.如何调用FISH文件？, t. I9 g7 L& o新手调用FISH文件出错的原因：是：文件名为 fishcall.fis.txt(因为在记事本里这个txt的后缀是默认的）。修改方法是：将 fishcall.fis.txt里面的内容全部复制到新的记事本里，再按保存，保存时候将保存类型改为所有文件(*,*),这时你便可在文件名栏内输入 fishcall.fis，类型也就正确了。2 G# x( N2 Re) T/ b% A28.如何理解Tension-p shear-p tension-n shear-n表示的含义？. w7 l q% K, C3 l4 m在计算的循环里面,每个循环中,每个zone都依据failure criterion处于不同的状态,shear 和tension分别表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即分别表示在现在和以前的循环当中处于塑性状态.8 r, d1 q# N4 q5 u/ h$ r E0 - 29. 绘制相关曲线命令的一些解释4 e# z6 u1 Q9 a1) 在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负);5 L0 D# _1 1 C- L2) -表示对其值作mirror,如果大家使用AUTOCAD的话就明白mirror的含义了,比如上面的 hist 1 gp ydisp的值是0100,那么vs -1就变成-1000.以此类推. E; k& y6 r 3) 为什么要这么做呢?主要是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限内. 5 L! z- iK: L9 f2 j) r, ?0 ; I* V* N: L7 + 1 u8 Y2 M2 w30. 如何查看剖面上任一点的位移值？7 p3 C# b) y0 ; print gp disp range x() y( ) z( )或! + P/ b0 r0 H! z Aprint gp disp range id - b q* d( H1 3 Tplot gp disp range x() y( ) z( )或1 T: f- F6 S( hg& E5 Z. ! fplot gp disp range id! i* ?6 P7 z5 * x& C, w4 V $ _8 C7 P1 ( I31. fix x y z range x -0.1 0.1 .a6 S7 t- q* q$ c! R fix x range x -0.1 0.1.b7 L2 D- B! |% A$ 请问a式与b式有何区别？( f R- C# Y7 W0 E. |9 w( a 表示 在x=0 这个平面 是固定的& x. J2 q2 |1 2 L+ ib 表示在x=0 这个平面 只x方向不能发生运动9 Z! q K/ S3 & v32. initial 与 apply 有何区别？$ R% G$ L( M5 g, M$ Q. Ninitial初始化命令，如初始化计算体的应力状态等； W. Z% C( - + r r7 1 U; 3 Papply边界条件限制命令，如施加边界的力、位移等约束等。 _5 a) G3 g) B+ y. Dinitial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件不会发生变化。6 I/ I4 ; 5 g5 Z33. 如何查看各个时段不平衡力的具体数值？: H z5 4 B2 R采用his来记录计算，包括位移应力等 命令his unbal / P) ( : M/ U5 C5 bhis gp(zone) zdis range (0 0 0) 或者id=？0 + z9 Q. 2 nY8 p. u导出数据命令his write n vs m begin 时步end时步 file filename.his0 B: o4 u+ Q( a# N: d9 o5 Z8 Kn表示纪录的id m表示时步7 C ( i, b/ b6 s: E0 I要导出不平衡力的具体数值5 L( V; I7 E3 Vy* B) ahis unbal, E. v8 pH# DT% E; Cstep 100000 or solve! N; D, y8 8 ehis write 1 vs step begin 1 end 1000 file 123.his. e4 I# t2 S1 R% gP使用上述命令就可以查看各个时步下的不平衡力的具体数值 # |) P) P7 a7 C: t) c! R& N, m$ R$ E% O, | 34. 如何显示变形轮廓线的命令？/ E: p, Y: ( d* G6 v! H+ O9 Z$ d- plo ske magf 10 其中10为放大系数2FLAC2d命令的FAQ：+ m. y4 b+ Y$ u1、怎样开始入门？ a# t9 P( B0 V, Q9 / X. |4 i4 s答：初步学习首先看说明书 Users gurid的Section2 Getting started和Section 3 Problem sloving with FLAC2 s# s+ ?$ g2 q0 s0 W4 2、哪些命令需要掌握？( S) ! P/ b& r答：grid gen prop interface ini fix set hist step solve pause cont sav 等，在学习的过程中遇到这些命令和不会的命令应该查看说明书的Command Reference，里面讲的很清楚，特别需要注意的是一定要仔细看一下Command Reference开头的Optional Range Phrases，这是讲范围的，因为很多命令的施加都需要加到模型的一部分，这里就告诉你怎样来表示这一部分。这一原则同样适合与其他5个ITASCA软件。还有一点就是六个软件（FLAC FALC3D UDEC 3DEC PFC PFC3D）的range大不相同 ，要注意区分。l$ E& u1 y2 9 F7 M8 l- I3、FLAC和FLAC3D建模有什么区别？: u u; ) f- : X2 L* e答：FLAC和FLAC3D基本的建模思想大致是一样的，不同的是FLAC是首先用GRID命令把整体的网格画出来，然后用MODELL NULL命令把网格分成许多部分，再用GEN命令把这些部分一块一块搭接起来。然后再用GEN LINE，GEN ARC等命令画出要开挖的边界，在FLAC中产生网格的命令只有一个GRID，而且只能用一次。而FLAC3D是直接用GEN 命令一块一块搭接，GEN可以用很多次。5 u9 Z3 M8 & I* ny O% D u4、怎样查看模型？& n; P 6 p) A答：plot grid 可以查看网格，plot grid num 可以查看节点号。& q3 % Cw$ J/ ) 5、怎样查看边界条件？2 D* o) a1 K% 1 Q答：plot fix 6 B) & o5 * D6 m! F0 6、怎样查看预先留设的开挖边界？% Y# q. d. o# R; E; |* S$ f答：plot grid mark。 F, n, ! k0 W: q6 n3 u* G b7、怎么查看模型的云图？% B9 U$ w6 S5 ?; s) 8 P- X答：pl syy (sxx xdis ydis sig1.) fill ，注意，带fill关键字就是查看云图，不带fill 关键字就是查看等直线图。# e n7 X, F* s+ G8、怎样查看模型的矢量图？5 v3 o0 q8 W% A2 j. j6 X b答：pl grid vel (disp) red0 D; e- , N4 v7 3 o9、怎样查看INTERFACE？8 Q r; . t, D( U 0 g答：pl grid iface id red + z% _( t5 n3 P print iface id ; id 是interface的编号, _1 & d( X4 e6 , Ui6 L10、FLAC怎么调整视图？7 h0 _0 z6 V( r# l# Z, Y* y答：ctrl+z 或shift+z ,然后点击要放大部位的左上角和右下角就可以放大某个部位，如果要缩小，双击屏幕。3 FLAC 使用步骤FLAC 使用步骤1.FLAC程式使用前准备步骤 步骤1：依比例画出所欲分析的资料在纸上画出地点的位置、地层资料、并简单标示距离及深度资料。 步骤2：换算输入资料成同一单位将现有地层资料，如 Density, Bulk modulus, Youngs modulus, tension, cohesion, friction Angle 等资料，换算成同一单位。附注 ：需谨慎检查输入资料的单位，如因单位不同而造成过大或过小的值，将会造成 FLAC无法计算，而产生ERROR信息。 步骤3：应用公式简略计算应用公式或依据经验，简略算出FLAC输出资料的范围，以作为Debug及输出资料分析时验证。步骤4：建立x, y坐标与node i, j之间的关系 在图上距离及深度的关系，建立x, y坐标系统，再由x, y坐标系统，转换与网格间系，为了便于以后输出资料的分析，故应确实掌握网格的位置及其相对应的x, y 坐标。建议在敏感区域使用较密的网格，其它地方则使用较疏的网格，刚开始执行程序时，不宜使用网格太大的网格数目，因尽量使网格总数少于1000，以节省时间。2.FLAC 3D 程序的编写步骤：1 Config_ 2 Grid_ 3 Model _4求起始的应力平衡 (1)建立x, y坐标与网格的关系，建议使用Gen指示：Gen x1,y1 x2,y2 ,x3,y3 x4,y4 i=i0,i1 j=j0,j1 详细指令参见使用手册，FLAC程式可自动产生x, y 坐标与网格的关系，但由于产生的网格坐标不易控制，将对其它的工作产生负面影响，故依使用前步骤4所建立的关系，将网格依其疏密程度需要的不同，实际控制网格的坐标。 (2)设定材料性质：prop (3)设定外力：Set Grav, Apply Pressure, ini sxx, Syy(4)设定边界条件：fix, free(5)求起始的应力平衡：solve (6)储存：Save 5求工程的影响 求出区域内的应力分布情况后，再依工程的流程及步骤阶段执