1、Flac 3D 常用命令总结new 不退出flac3d，新建一个新的分析计算任务generate 简写为ge 产生网格zone 简写为 zo 指示为3维网格体brick 简写为 b 指示为长方形的网格体size x y z 指示长方形网格体在x、y和z轴所划分的网格数plot 显示，输入之后回车会出现plot base/0 可以输入plot的任何子命令（add添加条目（surface 网格体颜色、axes 坐标系颜色）等），若回到flac3d命令状态，则按一次回车或输入qiut按x ，y， z 小写字母可以旋转当前视图中的X， Y ，Z 轴，相应的大写字母则方向相反。m 放大，M缩小。Ctrl
2、+G可以切换彩色图和灰色图。Model 简写为mo 定义材料模型Mohr 简写为 moh 指摩尔-库伦模型Property 简写为pro 定义模型的材料参数，材料参数包括：密度density 简写为de，体积模量bulk简写为bu，切变模量shear 简写为sh，内摩擦角friction 简写为fric，内聚力cohesion 简写为c 抗拉强度 tension 简写为ten等Set 用于设置模拟条件或控制条件，如重力加速度gravity，mechanical设置静态力学分析的有关参数 简写为mec，force 最大不平衡力简写为fo等。比如set mechanical force 50 In
3、itial 简写为in 初始化网格的相关值。Fix 保持网格节点指定参数（速度、压力、温度等）的值不变。Range指定范围History 简写为h 采样或记录迭代时模型中变量的值，检验采样记录时可以在plot子命令中输入history 1 、history 2等，忘记编号时可以输入print history就可以看到什么记录号对应什么变量或参数Nstep 简写为n 迭代次数，默认值等于10，ubalance 简写为unb 最大不平衡力 ，gp 采样网格节点的有关值，zdisplacement简写为zdis z轴方向的位移Solve 简写为so 控制相关过程的自动时间步一般等值线图（包括位移、速
4、度、应力、压力、温度等），可以在一个新的视图中显示，首先要在plot子命令中新创建一个视图，可以输入“create 视图名”，回车，继续输入add contour/bcontour disp/szz（contour为在当前视图中显示等值线图 简写为con ，bcontour 绘制指定区的等值线图 简写为bcon，disp 指定为位移量等值线图，szz指定为应力等值线图）任意剖面上的等值线图，可以在刚创建的视图中再创建一个视图然后输入剖面参数，比如“set plane dip 90 dd 0 origin 3 4 0”（plane 简写为p 设置一个剖平面，dip 为剖面倾角，dd 为剖面的倾向
5、，origin简写为o 为剖面中的一点）Boundary 简写为bo 在视图中增加面的边界线框Behind 简写为be 为当前视图后面要转换视图则在plot子命令中输入“current 视图名”、show开挖把网格体的材料模型设置成空（null）模型即可，命令如下：Model null range x y zSet large 把材料设置成大变形求解时计算不收敛则用step命令，step 步数节点的ID号和重心的全局坐标可以由下面的命令得到：Print gp position 单元体的ID号和重心的全局坐标可由如下的命令得到：Print zone所有命令大小写一样，都可以附带若干个关键字和有关
6、的值，格式如下：Command keyword valuekeyword，value命令、关键字、数字可以用数个空格隔开，也可以用（），=隔开， ； 为注释符， ； 开始到行末都是注释 一个简单的命令行只能至多80个字符（包括注释），如果超过80个字符则在行末输入&，每个命令最多含1024个字符。Macro 定义宏对象。若模型中有多次重复且含有又长又复杂字符的对象时，使用宏对象特别有用。使用方法如下：Macro pt0 p0 15 15 15Macro pt1 p1 10 0 0Macro pt2 p2 0 10 0Macro pt3 p3 0 0 10Macro model _size si
7、ze 4 5 6Macro big_brick zone brick pt0 pt1 pt2 pt3 model_sizeGen big_brick重新定义一个宏时，必须在宏名称前后放置单引号Macro pt0 p0 15 15 15Gen big_brick用不同的对象完成相同的工作：用一个range对象：Range name big_brick x=(-3,3) y=(-2,2) z=(-1,1)Model null range big_brick用一个group对象：Group big_brick range x=(-3,3) y=(-2,2) z=(-1,1)Model null ra
8、nge group big_brick用一个macro对象：Macro big_brick x=(-3,3) y=(-2,2) z=(-1,1)Model null range big_brick网格原始形状：见flac3d使用教程.pdf网格生成可以用四个点的坐标生成：ge zo b p0 p1 p2 p3 size x y z也可以用参数add（相对坐标）生成：Ge zo b p0 0 0 0 p1 add x 0 0 p2 0 y 0 p3 0 0 z size x y z FLAC/FLAC3D常规问题的整理1FLAC3D命令的FAQ/ ?9 Y0 y7 B1 x0 _( h+ E4
9、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等建模、初始条件、边界条件、后处理和
10、求解的命令。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*
11、 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. 怎样输出模型的后处理
12、图？答：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答：综合
13、使用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
14、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. 如何暂
15、停计算？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的内置语言，可以用来进行参数化模型、完成
16、命令本身不能进行的功能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
17、?( 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 lm
18、6 _ 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
19、 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文
20、件？, 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都依据failur
21、e 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 -
22、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
23、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.
24、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/ U
25、5 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 % 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为放大系数
