FLAC3D基础命令流解释
flac3d基础相关命令
Flac3D基础知识1.Generate 命令产生网格，简写为 GE。
2. zone 参数代表了三维网格体（单元体），简写为 zo。
3. （1）size 参数指示单元体各边小单元体的个数；（2）dimension 参数指示单元体内部巷道各边边长，如果没给出就默认为两点间距的20%；（3） edge 参数指示网格体各边的尺寸 Ge zone brick p0(a,b,c) edge x y z 边长为x y z 顶点为p0；（边长也可以用确定点坐标来利用点与点的间距来确定p0 p1 p2……）。
4.各种网格体名称及关键字：矩形网格（Brick），退化矩形网格（dbrick），楔形网格（wedge），金字塔网格（pyramid），四面体形网格（tetrahedron），圆柱体网格（cylinder），矩形体外环绕放射状网格（radbrick），平行六面体外环绕放射状网格（radtunnel），圆柱体外环绕放射状网格（radcylinder），柱形壳体网格（Cshell），交叉圆柱体网格（cylint），交叉平行六面体网格（tunint）。
5.绘制网格体语句：①Generate zone brick size a，b，c （其中的brick可以被其他的网格体替换）。
②Generate zone brick p0=（，，） &P1=（，，） &P2=（，，） & （边长由点点间距确定）P3=（，，） & → （ &符号为继续行的符号）Size a b c → （各边小单元体个数）③如果在p1 p2 p3 前面加 add 则各点坐标是相对于p0点的坐标。
④如果只给出p0 则由 edge a & 给出边长。
（见第 3 条）；利用ratio a，b，c 给出小单元体长度比例依次变化。
⑤内部有空单元体或者巷道的网格体，如果内部巷道需要填充为单元体，使用 fill group groupname 语句。
FLAC3D常见命令与使用技巧
FLAC3D常见命令与使用技巧FLAC3D常见命令与使用技巧1、FLAC3D常见命令：1.FLAC3D是有限元程序吗？答：不是！是有限差分法。
2.最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。
3.怎样看模型的样子？答：plo blo gro可以看到不同的group的颜色分布4.怎样看模型的边界情况？答：plo gpfix red5.怎样看模型的体力分布？答：plo fap red6.怎样看模型的云图？答：位移：plo con dis (xdis, ydis, zdis)应力：plo con sz (sy, sx,sxy, syz, sxz)7.怎样看模型的矢量图？答：plo dis (xdis, ydis, zdis)8.怎样看模型有多少单元、节点？答：pri info9.怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件10.怎样调用一个文件？答：File/call或者call命令10.如何施加面力？答：app nstress11.如何调整视图的大小、角度？答：综合使用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 / Restore17.如何暂停计算？答：Esc18.如何在程序中进行暂停，并可恢复计算？答：在命令中加入pause命令，用continue进行继续19.如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步20. Fish是什么东西？21. Fish是否一定要学？答：可以不用，需要的时候查Mannual 获得需要的变量就可以了22.FLAC3D允许的命令文件格式有哪些？答：无所谓，只要是文本文件，什么后缀都可以23.如何调用一些可选模块？答：config dyn (fluid, creep, cppudm)24 .如何在圆柱体四周如何施加约束条件？答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2radius r2 not，其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用set log on命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。
FLAC3D常见命令与使用技巧
FLAC3D常见命令与使用技巧1、FLAC3D常见命令：是有限元程序吗答：不是！是有限差分法。
2.最先需要掌握的命令有哪些答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。
3.怎样看模型的样子答：plo blo gro可以看到不同的group的颜色分布4.怎样看模型的边界情况答：plo gpfix red5.怎样看模型的体力分布答：plo fap red6.怎样看模型的云图答：位移：plo con dis (xdis, ydis, zdis)应力：plo con sz (sy, sx,sxy, syz, sxz)7.怎样看模型的矢量图答：plo dis (xdis, ydis, zdis)8.怎样看模型有多少单元、节点答：pri info9.怎样输出模型的后处理图答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件10.怎样调用一个文件答：File/call或者call命令10.如何施加面力答：app nstress11.如何调整视图的大小、角度答：综合使用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 / Restore17.如何暂停计算答：Esc18.如何在程序中进行暂停，并可恢复计算答：在命令中加入pause命令，用continue进行继续19.如何跳过某个计算步答：在计算中按空格键跳过本次计算，自动进入下一步20. Fish是什么东西答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学答：可以不用，需要的时候查Mannual获得需要的变量就可以了允许的命令文件格式有哪些答：无所谓，只要是文本文件，什么后缀都可以23.如何调用一些可选模块答：config dyn (fluid, creep, cppudm)24 .如何在圆柱体四周如何施加约束条件答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2radius r2 not，其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图答：用set log on命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。
FLAC3D命令流(整理版)
1、怎样查看模型？答：plot grid 可以查看网格，plot grid num 可以查看节点号。
2、请问在圆柱体四周如何施加约束条件？答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2<r1，其实就是选择range的时候选两个圆柱的差，即得到边界。
命令流如下：fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 13、怎么能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。
4、用命令建立模型后，如何显示点的坐标？答：使用plo blo gro gpnum on命令5、关于gauss_dev对性质进行高斯正态分布的问题？答：根据手册上的说明：下面的命令设定一个平均摩擦角为40度，标准方差是±5％。
则命令如下：prop friction 40 gauss_dev 2问题：请问gauss_dev 2中的2是如何计算的？如果把±5％改为±10％，则命令应如何写？40×5％＝26、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 fillplot surfgen zone reflect dip 0 dd 90 （对xy面做镜像）gen zone reflect dip 90 dd 90 （对yz面做镜像）（1）dd表示y轴正向顺时针到那条射线的夹角，dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线。
陈育民《FLAC3D基础与工程实例》全部命令流
《FLAC/FLAC 3D基础与工程实例》全部命令流1gen zone bri p0 0 0 0 p1 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_1gen 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_2gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1。
0 0。
9 1.1 group brick_3plot sur2—1newgen zon bri size 3 3 3model elasprop bulk 3e8 shear 1e8ini dens 2000fix z ran z -.1 .1fix x ran x —.1 .1fix x ran x 2。
9 3。
1fix y ran y —.1 。
1fix y ran y 2.9 3.1set grav 0 0 -10solveapp nstress -10e4 ran z 3 x 1 2 y 1 2hist gp vel 0 0 3hist gp vel 0 3 3plo hist 1 redplo add hist 2 bluesolve3—1；———----——----—--—-—-——-—-—-——--------工程信息;Project Record Tree export；Title:Simple test;-—-—-—-——-—---———-—---—-——--——--—计算第一步;。
.. STATE： STATE1 ...。
FLAC-命令流及其解释大全
FLAC-命令流及其解释大全1. FLAC3D是有限元程序吗？答：不是！是有限差分法。
2. 最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。
3. 怎样看模型的样子？答：plo blo gro可以看到不同的group的颜色分布4. 怎样看模型的边界情况？答：plo gpfix red sk5. 怎样看模型的体力分布？答：plo fap red sk6. 怎样看模型的云图？答：位移：plo con dis (xdis, ydis, zdis)应力：plo con sz (sy, sx, sxy, syz, sxz)7. 怎样看模型的矢量图？答：plo dis (xdis, ydis, zdis)8. 怎样看模型有多少单元、节点？答：plo info9. 怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件10. 怎样调用一个文件？答：File/call或者call命令10. 如何施加面力？答：app nstress11. 如何调整视图的大小、角度？答：综合使用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 / Restor17. 如何暂停计算？答：Esc18. 如何在程序中进行暂停，并可恢复计算？答：在命令中加入pause命令，用continue进行继续19. 如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步20. Fish是什么东西？答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学？答：可以不用，需要的时候查Mannual获得需要的变量就可以了22. FLAC3D允许的命令文件格式有哪些？答：无所谓，只要是文本文件，什么后缀都可以23. 如何调用一些可选模块？答：config dyn (fluid, creep, cppudm)后注：这个工作很繁琐，需要的时间很多，希望广大网友能够将自己曾经遇到的常见问题在后续跟贴，也为了将这个FAQ进行很好的充实。
FLAC3D命令流(整理版)
1、怎样查看模型？答：plot grid 可以查看网格，plot grid num 可以查看节点号。
2、请问在圆柱体四周如何施加约束条件？答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2<r1，其实就是选择range的时候选两个圆柱的差，即得到边界。
命令流如下：fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 13、怎么能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。
4、用命令建立模型后，如何显示点的坐标？答：使用plo blo gro gpnum on命令5、关于gauss_dev对性质进行高斯正态分布的问题？答：根据手册上的说明：下面的命令设定一个平均摩擦角为40度，标准方差是±5％。
则命令如下：prop friction 40 gauss_dev 2 问题：请问gauss_dev 2中的2是如何计算的？如果把±5％改为±10％，则命令应如何写？40×5％＝26、 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 surfgen zone reflect dip 0 dd 90 （对xy面做镜像）gen zone reflect dip 90 dd 90 （对yz面做镜像）（1）dd表示y轴正向顺时针到那条射线的夹角，dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线。
FLAC3D命令流(整理版)
1、怎样查看模型答：plotgrid可以查看网格，plotgridnum可以查看节点号。
2、请问在圆柱体四周如何施加约束条件答：可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not，其中r2<r1，其实就是选择range的时候选两个圆柱的差，即得到边界。
命令流如下：fixxrangeend1100end2140rad1end1100end2140rad13、怎么能把一个PLOT的图像数据导出来以便用其他软件绘图答：用setlogon命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。
4、用命令建立模型后，如何显示点的坐标答：使用ploblogrogpnumon命令5、关于gauss_dev对性质进行高斯正态分布的问题答：根据手册上的说明：下面的命令设定一个平均摩擦角为40度，标准方差是±5％。
则命令如下：propfriction40gauss_dev2问题：请问gauss_dev2中的2是如何计算的如果把±5％改为±10％，则命令应如何写40×5％＝26、reflect问题问：genzoneradbrick&p0(0,0,0)p1(10,0,0)p2(0,10,0)p3(0,0,10)&size3,5,5,7&rati o1,1,1,&dim142fillplotsurfgenzonereflectdip0dd90（对xy面做镜像）genzonereflectdip90dd90（对yz面做镜像）（1）dd表示y轴正向顺时针到那条射线的夹角，dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线。
首先应该按照dd 的方向大体确定这个面的朝向，dd指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角．（及z轴负方向转过角度）（2）命令：reflectnormalxyzoriginxyz（根据法线和过一点建的面）最常用，至于实际使用过程中有人问道：normal-100与normal100的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
;模型镜像
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 <add> contour szz outline on
;在模型中显示位移-应变曲线
hist gp ydisp 0,0,0
hist zone syy 0,1,0
hist zone syy 1,1,0
pl his -2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负); "-"表示对其值作"mirror"
;对模型进行压缩实验的方法
;即在模型两侧施加相反方向的速度
ini yvel 1e-7 range y -.1 .1
ini yvel -1e-7 range y 1.9 2.1
;修改模型的坐标值
ini x add -100 y add -100 z add -100
;显示云图的同时也显示模型网格轮廓
plot add cont disp outline on <gradient/average>;gradient更精确
;输入角度、弧度方法 pi=π，90°为90.0*degrad
def set_vals
ptA = 25.0 * sin(pi/2);ptA=25.0
ptB = 25.0 * cos( 60.0*degrad );ptB=12.5
ptC=pi;ptC=3.1415926
end
set_vals
print ptA ptB ptC
;施加结构单元方法
sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) &
end2=(0.0,25.0,0.0) radius=24.5 not
plot add sel geom black black cid on scale=0.03
sel node init zpos add -25.0
;如何显示某一平面
plot create name_plane
plot set plane origin 3 4 0 normal 1 0 0
plot add cont disp plane behind shade on
plot add sel geom black
plot add axes red
plot current name_plane
plot show
;如何施加法向应力并显示所施加的力
apply nstress -3.0E7 range cylinder end1 0 -1 0 end2010 radius 39.9 not;柱体外表面施加法向力
apply nstress -5.0E6 range cylinder end1 0 -1 0 end2010 radius 2.1;柱体内表面施加法向力
plot grid green fap black;显示施加的力
;获取某网格点的ID，坐标值
def aa
get_id=gp_id(gp_near(3,0,0))
end
aa
print get_id
print gp position range id 3
;保存命令流文件
set log on
set logfile hua.log;文件存放的路径与调用的文本路径相同
......
set log off
;fish函数的说明与解释
def find_add ;定义fish函数find_add
head = null ;给head赋值
p_gp = gp_head ;第一个网格结点的指针赋给p_gp
loop_while p_gp # null ;当p_gp值不为null时作循环
x_pos = gp_xpos(p_gp) ;将指针为p_gp的结点的x坐标值赋给x_pos
if x_pos = 10.0 then ;如果x_pos = 10.0 则(执行)
new = get_mem(2) ;从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址
mem(new) = head ;将head类型和数值置于地址为new的fish变量
mem(new+1)=p_gp ;第二个变量存储符合条件的zone地址
head = new ;将new值赋给head
endif
p_gp = gp_next(p_gp) ;将结点指针为p_gp的下一个结点的指针赋给p_gp endloop;结束循环
end ;结束fish函数
;实际上这个fish函数为满足条件（x坐标为10的）的zone的地址开辟一定的地址空间，各地址之间存在一定的联系
;找到第一个符合条件的zone地址后，用new = get_mem(2) 从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址，并用;下面的mem(new) = head ，将第一次开辟的两个变量的第一个变量存储地址head（注意第一次head ＝null），第二个变量存储第一;个符合条件的zone地址，并将第一个变量的地址赋予head（head＝new），第一次循环结束；下次循环，同样开辟两个变量对
象，第;一个变量对象记录上次循环开辟的第一个变量的地址，第二个记录第二个符合条件的zone地址，其余循环依次类推，这样子就建立了;一个符合条件的zone地址链条，方便以后使用。
；FLAC后处理
；输出JPG图片
set plot jpg
set plot quality 100
plot hard file XX.jpg
;输出TXT文件
hist write X file xx.txt ；(x是历史记录ID）
print con zd
;输出单元应力和节点位移
set log on
set logfile xx.log
print zo stress
print gp dis
set log off
；生成影片文件
plot set rot 20 0 30
plot con szz ou on magf 10
plot add hist 1
set movie avi step 1 file xx.avi
movie start
solve
;计算间隔命令在solve命令之前
set hist_rep 1 ；(步数）
；检查模型赋值情况
plot bcon prop bu
;变形云图加网格放大20 倍
plot con zd ou on magf 20
;变形网格
plot sk magf 20
;变形矢量图（加网格）
plot sk dis scale 0.07 ；(scale是箭头长度，默认0.05）
；如何观察塑性区
plot block state
plot block state sh_n
plot block state sh_p
;切片操作
plot set plane ori 1 2 3 norm 0 0 -1
plot con zd plane
plot add dis plane
plot add sk
plot add ax
;切片后方部分显示
plot con zd plane behind shade on
交集就是range group a group b
并集就是range group a any group b any not是差集。