flac命令流
1、FLAC3D常见命令：
1. FLAC3D是有限元程序吗？
答：不是！是有限差分法。
2. 最先需要掌握的命令有哪些？
答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。
3. 怎样看模型的样子？
答：plo blo gro可以看到不同的group的颜色分布
4. 怎样看模型的边界情况？
答：plo gpfix red
5. 怎样看模型的体力分布？
答：plo fap red
6. 怎样看模型的云图？
答：位移：plo con dis (xdis, ydis, zdis)
应力：plo con sz (sy, sx, sxy, syz, sxz)
7. 怎样看模型的矢量图？
答：plo dis (xdis, ydis, zdis)
8. 怎样看模型有多少单元、节点？
答：pri info
9. 怎样输出模型的后处理图？
答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件
10. 怎样调用一个文件？
答：File/call或者call命令
10. 如何施加面力？
答：app nstress
11. 如何调整视图的大小、角度？
答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键
12. 如何进行边界约束？
答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束）
13. 如何知道每个单元的ID？
答：用鼠标双击单元的表面，可以知道单元的ID和坐标
14. 如何进行切片？
答：plo set plane ori (点坐标) norm (法向矢量)
plo con sz plane (显示z方向应力的切片)
15. 如何保存计算结果？
答：save +文件名.
16. 如何调用已保存的结果？
答：rest +文件名；或者File / Restore
17. 如何暂停计算？
答：Esc
18. 如何在程序中进行暂停，并可恢复计算？
答：在命令中加入pause命令，用continue进行继续
19. 如何跳过某个计算步？
答：在计算中按空格键跳过本次计算，自动进入下一步
20. Fish是什么东西？
答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学？
答：可以不用，需要的时候查Mannual获得需要的变量就可以了
22. FLAC3D允许的命令文件格式有哪些？
答：无所谓，只要是文本文件，什么后缀都可以
23. 如何调用一些可选模块？
答：config dyn (fluid, creep, cppudm)
24 .如何在圆柱体四周如何施加约束条件？
答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2
25.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？
答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘
制就行了。
26.关于dd、dip的理解：
答：首先你要确定你定义面的法线方向，dip 是倾角的意思（xoy面和定义面的夹角），法线向上dip是锐角，法线向下dip 时钝角。dd是y轴顺时针转到法线在xoy面上的投影线时，所转过的角度。
27.如何调用FISH文件？
答：新手调用FISH文件出错的原因：是：文件名为fishcall.fis.txt(因为在记事本里这个txt 的后
缀是默认的）。修改方法是：将fishcall.fis.txt里面的内容全部复制到新的记事本里，再按保存
，保存时候将保存类型改为所有文件(*,*),这时你便可在文件名栏内输入fishcall.fis，类型也就
正确了。
28.如何理解Tension-p shear-p tension-n shear-n表示的含义？
答：在计算的循环里面,每个循环中,每个zone都依据failure criterion处于不同的状态,shear 和
tension分别表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即分别表示在现在和以前
的循环当中处于塑性状态.
29.绘制相关曲线命令的一些解释
在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负);
"-"表示对其值作"mirror",如果大家使用AUTOCAD的话就明白"mirror"的含义了
30.如何查看剖面上任一点的位移值？
答：print gp disp range x() y( ) z( )或
print gp disp range id
plot gp disp range x() y( ) z( )或
plot gp disp range id
31. a: fix x y z range x -0.1 0.1
b: fix x range x -0.1 0.1 请问a式与b式有何区别？
答：a 表示在x=0 这个平面是固定的；b 表示在x=0 这个平面只x方向不能发生运动
32.initial 与apply 有何区别？
答：初始化命令，如初始化计算体的应力状态等；
apply边界条件限制命令，如施加边界的力、位移等约束等。
initial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件
不会发生变化。
33.如何查看各个时段不平衡力的具体数值？
答：采用his来记录计算，包括位移应力等命令his unbal
his gp(zone) zdis range (0 0 0) 或者id=？
导出数据命令his write n vs m begin 时步end 时步+file filename.his
要导出不平衡力的具体数值
his unbal
step 100000 or solve
his write 1 vs step begin 1 end 1000 file 123.his
使用上述命令就可以查看各个时步下的不平衡力的具体数值
34.如何显示变形轮廓线的命令？
plo ske magf 10 其中10为放大系数
35、怎样开始入门？
答：初步学习首先看说明书User's gurid的Section2 Getting started和Section 3 Problem sloving with FLAC
36、哪些命令需要掌握？
答：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大不相
37、FLAC和FLAC3D建模有什么区别？
答：FLAC和FLAC3D基本的建模思想大致是一样的，不同的是FLAC是首先用GRID命令把整体的网格画出
来，然后用MODELL NULL命令把网格分成许多部分，再用GEN命令把这些部分一块一块搭接起来。然后
再用GEN LINE，GEN ARC等命令画出要开挖的边界，在FLAC中产生网格的命令只有一个GRID，而且只
能用一次。而FLAC3D是直接用GEN 命令一块一块搭接，GEN可以用很多次。
38、怎样查看模型？
答：plot grid 可以查看网格，plot grid num 可以查看节点号。
39、怎样查看边界条件？
答：plot fix
40、怎样查看预先留设的开挖边界？
答：plot grid mark。
41、怎么查看模型的云图？
答：pl syy (sxx xdis ydis sig1..) fill ，注意，带fill关键字就是查看云图，不带fill 关键字
就是查看等直线图。
42、怎样查看模型的矢量图？
答：pl grid vel (disp) red
43、怎样查看INTERFACE？
答：pl grid iface id red
print iface id ; id 是interface的编号
44、FLAC怎么调整视图？
答：ctrl+z 或shift+z ,然后点击要放大部位的左上角和右下角就可以放大某个部位，如果要缩小，
双击屏幕。
45、判断模型中网格点是一个还是两个的方法
显示模型网格点的ID，如果有两个ID号，说明是两个。节理面上的网格点应该有两个ID 号。
46、改变模型应力边界时FLAC与3DEC的不同之处
在3DEC中，当改变应力边界时，模型表面的应力直接被更换，FLAC中改变应力边界时，模型
表面的应力为两者应力相加。
47、flac中影响运行时间的因素
（1）单元材料，结构面之间的刚度对比很大；
（2）划分单元尺寸对比很大，尽量避免长细比大于5：1的细长单元。
48 FLAC中有多少中材料模型
内置11中本构模型%_4[ W/N,n R5q
（1）零模型用于开挖回填的区域
（2）各向同性弹性模型
（3）正交各向异性弹性模型
（4）横观各向同性弹性模型
（5）德鲁克－普拉格塑性模型
（6）莫尔－库仑塑性模型一般岩土体的力学行为
（7）节理化塑性模型用于研究材料破坏后行为（如边破稳定，屈服煤矿开采）
（8）应变硬化/软化模型－库仑塑性模型薄板层状材料破坏后力学行为
（9）双线性应变硬化/软化节理化塑性模型临近层状土层开挖
（10）双屈服塑性模型
（11）修正的剑桥粘土模型
49、对不平衡力的理解
flac中用最大不平衡力和典型内力大的比值R表示模型的不平衡力的相对大小，该值为百分数，且从
不会减小到零，因此，模型不会达到绝对稳定的平横状态，可以根据要求的精度设定R值，模型默认
的R值是1e-5。
注意R值小仅表示所有网格点的力都平衡，但模型还可能发生塑性流动，通过plot vel 命令查看
网格点速度，一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。
50、简单命令解释：
plot add cont disp out on 时时显示云图
plot block group gpnum on 显示节点编号
plot inter red ske 看接触面的效果.
del range group 1 not 不删除group信息
del range group 2 删除group信息
ini z add 数值大小range group name 顺z方向移动
建立接触面interface 1 face range x(范围）y(范围）z(范围）
或者interface 1 face range plane orig(x y z 坐标）norm（法线方向）
water table clear 清除水位线
pl set back white ;设置背景
set plot bitmap size (800,510) ;设置图片大小
plot set caption size 37 ; 设置图例数字大小(default=35 范围10-50)
plot set center 100 0 40
fix x y z range x -0.1 0.1 表示在x=0 这个平面是固定的fixed boundary
fix x range x -0.1 0.1表示在x=0 这个平面只x方向不能发生运动roller boundary
51、fish命令的解释：
def aux1 定义函数
loop n(1,nptab1) 开始一个循环，从1循环到变量nptab1
rr=xtable(1,n) 从table 1 的x列中读第n个数付给变量rr
zz=ytable(1,n) 从table 1 的y列中读第n个数付给变量zz
loop k(1,nprof) 开始第二个循环，从1循环到变量nprof
case_of k 下面是一个选择结构，利用变量k判断
i F 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 当k＝1的时候执行下面语句
xtable(n+10,k)=rr table n+10的x列第k个值等于变量rr
ytable(n+10,k)=-40.0 table n+10的y列第k个值等于变量-40
case 2 当k＝2的时候执行下面语句
xtable(n+10,k)=rr table n+10的x列第k个值等于变量rr
ytable(n+10,k)=0.0 table n+10的y列第k个值等于变量0
end_case 结束选择结构
end_loop 结束第二个循环
end_loop 结束第一个循环
end 函数定义结束
aux1 执行什么定义的函数
52. 在固定面上施加位移时，fix与apply不能同时使用，应该用ini和fix,才能施加上。
53.菜单plot items/add/structure elements/geometry/结构单元的形式，把对话框里的mag fac的值改成不同的数，就可以看到位移放大后的结构单元。
54.你可以用键盘上的上、下、左、右肩在窗口中移动模型。
55.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？
答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。
56 在进行三轴压缩时候，如何施加围压？
答：使用apply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius 1，其中X和Z 轴是在我们眼睛正对的面。
57 在进行三轴压缩时候，如何对轴向力进行位移控制？
答：使用ini yvel 1e-7 range y -.1 .1 ini yvel 1e-7 range y 1.9 2.1 即可，数值可以变化58 在进行边坡计算的时候，为什么用sol, fos 命令运行不行？
答：应该输入sol fos slope3dfos.sav associated(不要也行）再用restore 命令调用即可。
59 在FLAC3D计算中，如何记录程序中的值，如果不想时时记录该怎么弄？
答：用hist wtite n1 vs n2 file.txt，想间隔记录可用set hist_rep n。
60 想记录接触面上的相关力学参数该怎么办？
答：hist inter i nstress /ndis/sstress/sdisp
61 接触面设置参数包括哪几个？
答：coh dial fric kn ks ten.
62 如何输出单元格的数据？
答；expgrid .flacsd
63 如何在剖面的沉降云图上增加网格、剖面上的矢量及坐标值？D Q Z9s1} e4Y U%q 答：plo add ske
plo add dis plane
plo add axe
64 如何做切片？
答：plot set plane ori x y z norm x y z
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 求工程的影响________
求出区域内的应力分布情况后，再依工程的流程及步骤阶段执行各工程进行过程的影响，建议使用以下的步骤：
(1)调出起初的应力平衡：re_____ .sav
(2)设定新的材料性质：model,prop
(3)设定新的支撑性质：struct
(4)设定新的外力
(5)设定边界条件
(6)求工程时的应力平衡
(7)储存
FLAC命令流
FLAC命令流 ;**************************************** group命令 group soil range z 1 2 ;定义group prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group soil ;给group赋值model null range group soil model elastic range group soil range命令 range name trench x 0 1 y 0 4 z 0 2 model null range trench model null range x=2,4 y=2,6 z=5,10 range name Big_Brick x -3 3 y -2 2 z -1 1 model elastic range Big_Brick prop bulk 1e8 shear 1e8 range Big_Brick range name Layer1 plane dip 0 dd 0 ori 0 0 0 above range name Layer2 plane dip 0 dd 0 ori 0 0 0 below macro命令 macro Sand 'bulk 1e8 shear 0.5e8 coh 0 tens 0 fric 35' macro Clay 'bulk 1e7 shear 0.3e7 coh 1e7 tens 0 fric 0' prop sand range Layer1 prop clay range layer2 macro Pt0 'p0 0 0 0' macro Pt1 'p1 add 10 0 0' macro Pt2 'p2 add 0 10 0' macro Pt3 'p3 add 0 0 10' macro Model_Size 'size 4 5 6' macro Big_Brick 'gen zone brick Pt0 Pt1 Pt2 Pt3 Model_Size' Big_Brick macro 'Pt0' 'p0 15 15 15' gen Big_Brick ; this will cause an error
flac 命令流 参考 实用
;**************************************** group命令 group soil range z 1 2 ;定义group prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group soil ;给group赋值 model null range group soil model elastic range group soil plot block group range group dam ;只显示dam range命令 range name trench x 0 1 y 0 4 z 0 2 model null range trench model null range x=2,4 y=2,6 z=5,10 range name Big_Brick x -3 3 y -2 2 z -1 1 model elastic range Big_Brick prop bulk 1e8 shear 1e8 range Big_Brick range name Layer1 plane dip 0 dd 0 ori 0 0 0 above range name Layer2 plane dip 0 dd 0 ori 0 0 0 below range cylinder end1 x1 y1 z1 end2 x2 y2 z2 radius r cylindrical range with one end of the cylinder axis (end1) at location (x1, y1, z1), the other end (end2) at location (x2, y2, z2), and with a cylinder radius of r '由(x1, y1, z1) 、(x2, y2, z2)两点确定旋转轴 macro命令 macro Sand 'bulk 1e8 shear 0.5e8 coh 0 tens 0 fric 35' macro Clay 'bulk 1e7 shear 0.3e7 coh 1e7 tens 0 fric 0' prop sand range Layer1 prop clay range layer2 macro Pt0 'p0 0 0 0' macro Pt1 'p1 add 10 0 0' macro Pt2 'p2 add 0 10 0' macro Pt3 'p3 add 0 0 10' macro Model_Size 'size 4 5 6'
FLAC3D 实例命令流1
第1部分命令流按照顺序进行 2-1定义一个FISH函数 new def abc abc = 25 * 3 + 5 End print abc 2-2使用一个变量 new def abc hh = 25 abc = hh * 3 + 5 End Print hh Print abc 2-3对变量和函数的理解 new def abc hh = 25 abc = hh * 3 + 5 End set abc=0 hh=0 print hh print abc print hh new def abc abc = hh * 3 + 5 end set hh=25 print abc set abc=0 hh=0 print hh print abc print hh 2-4获取变量的历史记录 new gen zone brick size 1 2 1 model mohr prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10 fix x y z range y -0.1 0.1 apply yvel -1e-5 range y 1.9 2.1 plot set rotation 0 0 45
plot block group def get_ad ad1 = gp_near(0,2,0) ad2 = gp_near(1,2,0) ad3 = gp_near(0,2,1) ad4 = gp_near(1,2,1) end get_ad def load load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) end hist load hist gp ydis 0,2,0 step 1000 plot his 1 vs -2 2-5用FISH函数计算体积模量和剪砌模量 new def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive print b_mod print s_mod 2-6 在FLAC输入中使用符号变量 New def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive gen zone brick size 2,2,2 model elastic prop bulk=b_mod shear=s_mod print zone prop bulk print zone prop shear 2-7 控制循环 New def xxx sum = 0 prod = 1 loop n (1,10)
FLAC3D基础命令流解释
;模型镜像 gen zone radcylinder size 25 1 25 25 gen zone reflect normal -1 0 0 origin x y z(面上一点）;沿X轴镜像,通过对称平面法线向量确定对称面 gen zone reflect normal 0 0 -1 ;沿z轴镜像 ;绘图控制 pl
FLAC3D命令流(整理版)
实用标准文档 1、怎样查看模型？ 答：plot grid 可以查看网格，plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件？ 答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2 F L A C D命令流挺实用 Revised final draft November 26, 2020 1、怎样查看模型 答：plot grid 可以查看网格，plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件 答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2 6、 reflect问题 问：gen zone radbrick &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轴正向顺时针到那条射线的夹角，dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线。首先应该按照dd的方向大体确定这个面的朝向，dd 指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向 在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角．（及z轴负方向转过角度） （2）命令：reflect normal xyz origin xyz （根据法线和过一点建的面）最常用，至于实际使用过程中有人问道：normal -1 0 0 与normal 1 0 0的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面。reflect dd a dip b origin xyz 这个就是用里建立斜面镜像的。 有一点注意了，镜像的时候如果是模型镜像的话，相当于在建立一个模型，故，在reflect的前面要加建模的关键字，gen zone reflect…… 1、怎样查看模型？ 答：plot grid 可以查看网格，plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件？ 答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2 f l a c d软件隧道支护与开 挖命令流 The latest revision on November 22, 2020 ;----------------------------------------------------- ; ---- Excavation and Support for a Shallow Tunnel --- ;----------------------------------------------------- new ;新建项目 set fish autocreate off title 'Excavation and Support for a Shallow Tunnel' ;定义题目 ; generate primitive components of grid ; concrete liner - upper tunnel gen zon cshell p0 0 0 0 p1 7 0 0 p2 0 51 0 p3 0 0 5.5 & dim 5 5 5 5 size 2 51 10 group zone 'concrete liner' ; ; upper tunnel gen zon cylinder p0 0 0 0 p1 5 0 0 p2 0 51 0 p3 0 0 5 & size 5 51 10 group zone tunnel range group 'concrete liner' not ; ; lower tunnel & liner gen zone brick p0 0 0 -4.5 p1 add 7 0 0 p2 add 0 51 0 p3 add 0 0 4.5 & size 7 51 3 ; ; surrounding rock (8 primitives) gen zon radcyl p0 0 0 0 p1 27 0 0 p2 0 51 0 p3 0 0 25 & dim 7 5.5 7 5.5 size 5 51 10 8 rat 1 1 1 1.3 ; gen zone brick p0 7 0 -4.5 p1 27 0 -15 p2 add 0 51 0 p3 7 0 0 & p4 27 51 -15 p5 7 51 0 p6 27 0 0 p7 27 51 0 & size 8 51 3 ratio 1.3 1 1 ; gen zone brick p0 0 0 -15 p1 add 27 0 0 p2 add 0 51 0 p3 0 0 -4.5 & p4 27 51 -15 p5 0 51 -4.5 p6 7 0 -4.5 p7 7 51 -4.5 & size 7 51 8 rat 1 1 0.7692307692307692 ; gen zon brick p0 0 0 25 p1 add 27 0 0 p2 add 0 51 0 p3 add 0 0 10 & size 5 51 2 ; gen zon bric p0 27 0 25 p1 add 17 0 0 p2 add 0 51 0 p3 add 0 0 10 & size 2 51 2 rat 2 1 1 ; 第1部分命令流按照顺序进行2-1定义一个FISH函数 new def abc abc = 25 * 3 + 5 End print abc 2-2使用一个变量 new def abc hh = 25 abc = hh * 3 + 5 End Print hh Print abc 2-3对变量和函数的理解 new def abc hh = 25 abc = hh * 3 + 5 End set abc=0 hh=0 print hh print abc print hh new def abc abc = hh * 3 + 5 end set hh=25 print abc set abc=0 hh=0 print hh print abc print hh 2-4获取变量的历史记录 new gen zone brick size 1 2 1 model mohr prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10 fix x y z range y -0.1 0.1 apply yvel -1e-5 range y 1.9 2.1 plot set rotation 0 0 45 plot block group def get_ad ad1 = gp_near(0,2,0) ad2 = gp_near(1,2,0) ad3 = gp_near(0,2,1) ad4 = gp_near(1,2,1) end get_ad def load load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) end hist load hist gp ydis 0,2,0 step 1000 plot his 1 vs -2 2-5用FISH函数计算体积模量和剪砌模量 new def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive print b_mod print s_mod 2-6 在FLAC输入中使用符号变量 New def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive gen zone brick size 2,2,2 model elastic prop bulk=b_mod shear=s_mod print zone prop bulk print zone prop shear 1、怎样查瞧模型？ 答:plot grid 可以查瞧网格,plot grid num 可以查瞧节点号。 2、请问在圆柱体四周如何施加约束条件？ 答:可以用fix 、、、ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2 建模 1、调用文件： ①文件与工程在同一个文件夹，只写文件名即可：If the called file is located in the same folder as the FLAC3D project file, then only the file name need be entered with the CALL command. ②不在同一个文件夹，全路径：Otherwise, the file may be called by specifying its complete path (e.g., c:\myfolder\file.dat). Undo ;撤销上一条命令 2、创建旋转缩放视图 3、建模命令 model mech mohr ;莫尔库伦模型 model mech elastic ;弹性模型 set grav 0, 0, -9.81 ; 重力加速度negative z-direction. （垂直向下！常用的） gen zone brick size 6,8,8 p0 -10, -10, -20 ... ;省略号表示写不下后面继续 p1 10, -10, -20 ... p2 -10, 10, -20 ... p3 -10, -10, 0 plot zone gen zone brick size 6,8,8 p0 -10, -10, -20 ... ;不规则六面体 p1 10, -10, -20 p2 -10, 10, -20 ... p3 -10, -10, 0 p4 10, 10, -20 ... p5 -10, 10, 10 p6 10, -10, 0 ... p7 10, 10, 10 plot current plot Plot01 plot clear plot zone Undo ;撤销命令 set logfile 127X1001.txt set log on truncate set log off list zone prin range x 0 1 y 0 1 z 0 1 ;显示指定范围内各单元的主应力，结果如下 命令流 1 gen zone bri p0 0 0 0 pl 10 0 0 p2 0 10 0 p3 0 0 10 & p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 & size 10 10 10 rat 1.0 0.9 1.1 group brick_1 gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 & size 10 10 10 rat 1.0 0.9 1.1 group brick_2 gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0.9 1.1 group brick_3 plot sur 2-1 new gen zon bri size 3 3 3 model elas prop bulk 3e8 shear 1e8 ini dens 2000 fix z ran z -.1 .1 fix x ran x -.1 .1 fix x ran x 2.9 3.1 fix y ran y -.1 .1 fix y ran y 2.9 3.1 set grav 0 0 -10 solve app n stress -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3 hist gp vel 0 3 3 plo hist 1 red plo add hist 2 blue solve 3-1 ; ---------------------- 工程信息 ;Project Record Tree export ;Title:Simple test ; -------------------- 计算第一步 ;... STATE: STATE1 .... config grid 10,10 model elastic group 'User:Soil' notn ull model elastic notnull group 'User:Soil' prop den sity=1500.0 bulk=3E6 shear=1E6 not null group 'User:Soil' fix x y j 1 fix x i 1 fix x i 11 set gravity=9.81 history 999 un bala need solve save state1.sav ; -------------------- 计算第二步 ;... STATE: STATE2 .... in itial xdisp 0 ydisp 0 in itial xvel 0 yvel 0 model null i 4 7 j 8 10 group 'n ull' i 4 7 j 8 10 group delete 'n ull' history 1 xdisp i=4, j=11 solve save state2.sav ; ------------------- 绘图命令 ;*** plot comma nds **** ;plot n ame: syy plot hold grid syy fill ;plot n ame: Un bala need force plot hold history 999 ;plot n ame: grid plot hold grid magnify 20.0 lred grid displaceme nt ;plot n ame: Xdis-A plot hold history 1 line 5- 1 new ;定义球体半径和半径方向上单元网格数 def parm rad=10.0 rad_size=5 end parm ;建立八分之一球体外接立方体网格 F L A C D命令流挺实用 集团标准化办公室：[VV986T-J682P28-JP266L8-68PNN] 1、怎样查看模型 答：plot grid 可以查看网格，plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件 答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2 6、 reflect问题 问：gen zone radbrick &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轴正向顺时针到那条射线的夹角，dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线。首先应该按照dd的方向大体确定这个面的朝向，dd 指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向 在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角．（及z轴负方向转过角度） （2）命令：reflect normal xyz origin xyz （根据法线和过一点建的面）最常用，至于实际使用过程中有人问道：normal -1 0 0 与normal 1 0 0的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面。reflect dd a dip b origin xyz 这个就是用里建立斜面镜像的。 有一点注意了，镜像的时候如果是模型镜像的话，相当于在建立一个模型，故，在reflect的前面要加建模的关键字，gen zone reflect…… 1、FLAC3D常见命令： 1. FLAC3D是有限元程序吗？ 答：不是！是有限差分法。 2. 最先需要掌握的命令有哪些？ 答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。 3. 怎样看模型的样子？ 答：plo blo gro可以看到不同的group的颜色分布 4. 怎样看模型的边界情况？ 答：plo gpfix red 5. 怎样看模型的体力分布？ 答：plo fap red 6. 怎样看模型的云图？ 答：位移：plo con dis (xdis, ydis, zdis) 应力：plo con sz (sy, sx, sxy, syz, sxz) 7. 怎样看模型的矢量图？ 答：plo dis (xdis, ydis, zdis) 8. 怎样看模型有多少单元、节点？ 答：pri info 9. 怎样输出模型的后处理图？ 答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件 10. 怎样调用一个文件？ 答：File/call或者call命令 10. 如何施加面力？ 答：app nstress 11. 如何调整视图的大小、角度？ 答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键 12. 如何进行边界约束？ 答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束） 13. 如何知道每个单元的ID？ 答：用鼠标双击单元的表面，可以知道单元的ID和坐标 14. 如何进行切片？ 答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 15. 如何保存计算结果？ 答：save +文件名. 16. 如何调用已保存的结果？ 答：rest +文件名；或者File / Restore 17. 如何暂停计算？ 答：Esc 18. 如何在程序中进行暂停，并可恢复计算？ 答：在命令中加入pause命令，用continue进行继续 19. 如何跳过某个计算步？ 答：在计算中按空格键跳过本次计算，自动进入下一步 FLAC3D 常用命令流 ;模型镜像 gen zone radcylinder size 25 1 25 25 gen zone reflect normal -1 0 0 origin x y z(面上一点）;沿X轴镜像,通过对称平面法线向量确定对称面 gen zone reflect normal 0 0 -1 ;沿z轴镜像 ;绘图控制 pl gen zone brick size 6 8 8 p0 -10 -10 -20 & p1 10 -10 -20 & p2 -10 10 -20 & p3 -10 -10 0 plot surf plot his -2 vs 1 为输出一二的关系图，如果应力值都为负数，则加-2是为了应力坐标轴为正plot add cont szz out on shade on 打开szz的等值线图 plo con zdisp outline on shade on 打开zdisp的等值线云图 plot set magf 1.0视图的放大倍数为1.0 plo con szz z方向应力云图 gen merge 1e-5 range z 0此命令是接触面单元合并成一个整体，1e-5是容差 渐变应力施加：apply nstress -1e6 gradient 0,0,1e5 range z 3.464,0 plane dip 60 dd 270 origin .1 0 0；施加法向应力：apply nstress -1e6 range plane dip 60 dd 270 origin .1 0 print gp position range id=14647 输出节点坐标 print zone state都可以用来显示模型的塑性区 Hist gp zdisp 1；Hist gp szz 2输出2和1的hist关系图用：plot his -2 vs 1 如果应力值都为负数，则加-2是为了应力坐标轴为正！或者hist id=？ Gp zdis x y z 是记录xyz坐标的位移值并顶为ID=？ print zone stress 在命令窗口中会输出模型全部单元的6个方向应力数值 print gp dis 该命令会显示模型中全部节点三个方向的变形大小 hist unbal监测不平衡力，并保留历史记录 hist gp zdisp 4,4,8监测网格坐标点（4，4，8）在z方向的位移，并保留历史纪录 HIST gpvelFLACD命令流挺实用
FLAC3D命令流(挺实用)
flacd软件隧道支护与开挖命令流
FLAC3D 实例命令流1
FLAC3D命令流挺实用
flac3d 5.0常用命令集锦
FLAC3D命令流
FLACD命令流挺实用完整版
flac命令流
Flac3D常用命令流Word版
flac命令流整理