FLAC3D常见命令与使用技巧

合集下载

FLAC3D命令流（整理版）

FLAC3D命令流（整理版）1、怎样查看模型？答：plot grid 可以查看网格，plot grid num 可以查看节点号。

2、请问在圆柱体四周如何施加约束条件？答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2<r1,其实就是选择range的时候选两个圆柱的差，即得到边界。

命令流如下：< p="">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 bio gro gpnum on 命令5、关于gauss_dev对性质进行高斯正态分布的问题？答：根据手册上的说明：下面的命令设定一个平均摩擦角为40 度，标准方差是± 5％。

则命令如下：prop friction 40 gauss_dev 2 问题：请问gauss_dev 2 中的2 是如何计算的？如果把± 5％改为± 10%，则命令应如何写？40 X 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常见命令与使用技巧

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是什么东西？答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能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常用命令总结第一篇：FLAC3D常用命令总结1、模型建好如何检验所建模型是否有问题注意：第一步建模一定不要建太大的模型，可以建非常小的模型，确保模型的准确性，这样可以尽快开始后边的步骤，等后边计算分析步骤快速做完验证没有大问题，再将准确的模型建好换上，重新计算得到更加准确的结果。

2、如何检验应力边界、位移边界、加载速率等是否准确Plot add gpfix：看位移边界Plo add fap：看应力边界，看加载方向是否准确 Plot add vel：看速率，加载速率是否准确3、如何查看运算结果Plo con szz outline on：看z向应力云图，并显示边界网格 Plo add disp：看位移Plo con zdisp outline on magf 5：看z向位移云图，变形扩大5倍 Plo add bo ：添加边界 Plo blo state：看是否屈服Plo add stensor：查看大主应力方向Plo con smin：小主应力Plo con smax：大主应力Plo con ssi（shear strain increment）：查看剪切应变增量云图Plo bcont smin：查看单元小主应力4、如何做切面并查看切面结果在编辑“Edit“菜单中设置切面位置，在上述3命令后加上plane （pl）即可。

如：Plo con smin pl：查看切面小主应力5、后处理将模型导入到tecplot中做切面，加等值线等，对结果进一步处理，得到更漂亮的一些结果。

其他重要命令总结：1、如何在圆柱体四周如何施加约束条件？可以用fix...ran cylinder end1 end2 radius r1 cylinder end1end2 radius r2 not，其中r22、如何能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

(完整版)FLAC3D常用命令

1. apply（缩写：app）可用来定义边界条件及初始条件：1）添加应力格式1：apply szz -0.3395e6 range z -0.1 0.1格式1：apply szz -0.3395e6 range group pile格式3：apply szz -0.3395e6 range z -0.1 0.1 group pile格式4：apply nstress 数值range z 2.9 3.1（或3）x 1 2 y 1 22）以一定速度施加位移边界格式1：apply yvel -1e5 range y -1.9 2.1 ;施加y方向速度-1e5/step3）添加边界条件格式：apply szz <数值> grad <梯度> range <范围>示例1：apply szz -1e9 grad 0 0 8.3e5 range z 0 120示例2：apply szz -0.6e6 range z 0.05 0.15 group pile注：<数值>是梯度方向坐标0点的数值，可通过定义坐标范围的上下值与梯度计算得到。

2. range(缩写：ran)通过range功能，可以使命令作用在一定指定范围的目标上；如果一个命令没有使用range来确定范围，则命令对整个模型有效。

1)利用坐标指定一定的范围格式1：range z 0 1格式2：range z 2.9 3.1 x 1 2 y 1 22)利用分组来指定范围格式：range group 13)以上两种的复合格式：range z -0.1 0.1 group pile4）利用id号来指定一定的范围格式：range id 0 10该命令后跟起始id和结束id，这里的id可以是实体单元、网格、结构单元、接触面和节点的编号。

例：model elastic range id 1 10 ;指定id为1到10的单元为各向同性弹性本构。

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常用命令

1. apply（缩写：app)可用来定义边界条件及初始条件：1）添加应力格式1：apply szz -0。

3395e6 range z -0.1 0。

1格式1：apply szz -0.3395e6 range group pile格式3：apply szz —0.3395e6 range z —0。

1 0.1 group pile格式4：apply nstress 数值 range z 2.9 3.1（或3） x 1 2 y 1 22）以一定速度施加位移边界格式1:apply yvel —1e5 range y -1.9 2.1 ；施加y方向速度—1e5/step3）添加边界条件格式：apply szz <数值> grad <梯度〉 range 〈范围>示例1：apply szz -1e9 grad 0 0 8。

3e5 range z 0 120示例2:apply szz -0。

6e6 range z 0.05 0。

15 group pile注:〈数值〉是梯度方向坐标0点的数值，可通过定义坐标范围的上下值与梯度计算得到.2. range(缩写：ran)通过range功能，可以使命令作用在一定指定范围的目标上；如果一个命令没有使用range 来确定范围，则命令对整个模型有效。

1)利用坐标指定一定的范围格式1：range z 0 1格式2:range z 2。

9 3.1 x 1 2 y 1 22)利用分组来指定范围格式：range group 13）以上两种的复合格式：range z -0.1 0。

1 group pile4）利用id号来指定一定的范围格式：range id 0 10该命令后跟起始id和结束id，这里的id可以是实体单元、网格、结构单元、接触面和节点的编号。

例：model elastic range id 1 10 ;指定id为1到10的单元为各向同性弹性本构。

FLAC3D命令流(整理版)

1、怎样查看模型？答：plot grid 可以查看网格，plot grid num 可以查看节点号。

则命令如下：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的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面。

FLAC3D命令流(挺实用)

1、怎样查看模型？答：plot grid 可以查看网格，plot grid num 可以查看节点号。

Flac3D常用命令及常规问题总结

Flac3D常用命令及常规问题总结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+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 50Initial 简写为in 初始化网格的相关值。

Fix 保持网格节点指定参数（速度、压力、温度等）的值不变。

Range 指定范围History 简写为h 采样或记录迭代时模型中变量的值，检验采样记录时可以在plot子命令中输入history 1 、history 2等，忘记编号时可以输入print history就可以看到什么记录号对应什么变量或参数Nstep 简写为n 迭代次数，默认值等于10，ubalance 简写为unb 最大不平衡力，gp 采样网格节点的有关值，zdisplacement简写为zdis z轴方向的位移Solve 简写为so 控制相关过程的自动时间步一般等值线图（包括位移、速度、应力、压力、温度等），可以在一个新的视图中显示，首先要在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 为剖面的倾向，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所有命令大小写一样，都可以附带若干个关键字和有关的值，格式如下：Command 〈keyword value…〉[keyword，value…]…命令、关键字、数字可以用数个空格隔开，也可以用（），=隔开，；为注释符，；开始到行末都是注释一个简单的命令行只能至多80个字符（包括注释），如果超过80个字符则在行末输入&，每个命令最多含1024个字符。

FLAC3D常见命令与使用技巧

FLAC3D常见命令与使用技巧1、FLAC3D常见命令：1.FLAC3D是有限元程序吗？答：不是！是有限差分法。

2.最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。

(总结452类)FLAC3D命令流挺实用

、怎样查看模型？答：可以查看网格，可以查看节点号。

、请问在圆柱体四周如何施加约束条件？答：可以用... ，其中<，其实就是选择的时候选两个圆柱的差，即得到边界。

命令流如下：、怎么能把一个的图像数据导出来以便用其他软件绘图？答：用命令，把数据导出来，转到里处理一下，然后用或者什么作图软件绘制就行了。

、用命令建立模型后，如何显示点的坐标？答：使用命令、关于对性质进行高斯正态分布的问题？答：根据手册上的说明：下面的命令设定一个平均摩擦角为度，标准方差是±％。

则命令如下：问题：请问中的是如何计算的？如果把±％改为±％，则命令应如何写？×％＝、问题问：() () () ()（对面做镜像）（对面做镜像）（）表示轴正向顺时针到那条射线的夹角，表示对称参照面与平面的夹角，对称参照面与平面的夹角在平面的投影是一条射线。

首先应该按照的方向大体确定这个面的朝向，指的是从轴正方向按顺时针（）方向转向所要确定面的法线方向在平面上的投影的夹角，然后再确定，指的是从平面转向所要确定的平面的角．（及轴负方向转过角度）（）命令：（根据法线和过一点建的面）最常用，至于实际使用过程中有人问道：与的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面。

这个就是用里建立斜面镜像的。

有一点注意了，镜像的时候如果是模型镜像的话，相当于在建立一个模型，故，在的前面要加建模的关键字，……7、关于两种是否有区别？和答：是一样的，表示并列连接关系，及的意思。

8、的区别答：初始化命令，初始化计算应力状态，定义材料参数、重力等。

边界条件限制命令，如施加边界的力、位移等约束等。

用于设置计算的模式，控制计算过程。

施加的初始力在计算过程中是会变的，一般体力需要初始化，而施加的力、边界条件不会随着计算变化。

命令作用于整个网格(赋值全局变量)则初始化指定范围内网格点、域的偏移量。

初始应力不仅仅是重力引起的。

FLAC-3D小技巧

1、如果是外观上看模型，要切换视角的话除了上述的快捷键，命令是：plot set rotation 20 0 60 rotation后面的3个数字是绕x、y、z轴转动的角度。

使用快捷键时，每次操作转动的角度，或者move的大小，可以通过plot set rotincrement和plot set moveincrement来设置。

2、从surfer把数据倒入FLAC3D!不知道各位在做象边坡地形表面这种不规则面的模型时有没有什么好招。

我最近发现了一种方法。

在FLAC3D的fish函数库中有一个叫TOPO的函数，可以生成自由表面模型的网格，关键是它要用到的那个数据表做起来非常麻烦，具体就是在地形图上平行y轴画一系列的等距的直线，然后读出每条与每条直线相交的等高线的x值和高程，然后生成一张表，格式是table y x z。

用surfer可以很方便的生成这张表，在DATA菜单下的convert...可以将GRD文件生成以上要求的那种格式。

3、按找手册中的解释，print gp position 是用来输出gp的位置坐标的；print zone gp是用来输出与zone相关的gp的id值的。

这样我就有一个上面的疑问，两个命令中的gp都有id值，是不是不同命令中的gp的id值相同的话，他们的坐标就可以看作是一样的啊？答：如果要画等值线的话，先输出gp的坐标来，将其另存为一个文件，再输出zone的应力值，将与zone对应的gp的id与前面输出的gp的id值对应起来，找到gp的坐标。

用excel打开上面的文件进行编辑，最后用surfer画等值线。

这样的办法太苯，不知到有没有更好的办法啊？4、plot stensor plane local on out off能够显示的主应力矢量图，但仍有两个问题1）显示的矢量图中应力好像没有变化，如何通过设置符号的大小表示应力的大小变化？2）应力符号好像总在图形没有发生变形前的单元中心，如何使应力符号在发生变形了的单元中心？5、个人经验，无论有限元还是flac，如果支护措施仅仅只考虑其刚度的贡献，岩体的变形和塑性区变化往往都不大，可能需要适当提高岩体的参数，才能反映支护的作用6、＜锚杆支护的数值模拟＞．中国矿业大学岩石力学与土木工程学报前段时间有篇文章专门关于flac锚杆单元应用模型修改的文章，7、如何模拟预应力自由式锚索？首先我想说说我在模拟预应力自由式锚索的做法,然后在说出现的问题,希望大家帮忙解决一下!我采用先整根生成的方法来生成锚索,然后给自由段和锚固段分别赋不同的参数,自由段的灌浆参数统统赋为零,锚固段的按试验值确定.在模拟垫墩时,把自由段外端的交接删除,然后重新生成刚性铰接(之所以这样做是因为模型介质较软,不这样做就收敛不了,介质内发生了流变.而全长粘结式锚索不需这样处理也和试验值拟合的很好),使得该结点不会和周围介质发生相对位移;然后在垫墩作用的面积上再施加一个均布力,大小等于预应力除以垫墩面积.问题是位移值拟合锝还可以,但是锚索的预应力增加值比实测的大的多,有没有其他的方法来模拟自由式锚索???8、Fish函数增加了以下几种新的特性：1.增加了fish变量来获取结点、单元和界面变量2.fish提供了获取结构单元变量的途径3.休单元和面单元性质目前可以通过单元变量名z_prop(i_z string)和界面单元名i_prop(i_z string)分别加以识别4.fish函数可以获取单元应变和应力速率，还提供了全应变增量张量和应变速率张量5.提供了fish绘图子程序函数够生用户定义的图形内容6.fish函数已经增加了从文件读、写数据的能力FLAC3D是一个强大的软件，但是不得不承认的，它的界面没有ansys好用，商业运行不是很强的说注意点：1.fish函数可以嵌套使用;2.以save命令保存模型时，fish函数和变量也同时保存；3.fish函数不支持缩写，这与flac3d命令不同，另外所有的fish函数或变量不区分大小写，程序同意转化为大写进行编译，当然也可以通过执行set case-sensitivity on来区分大小写；4.变量或函数名不能以一个数字开头也不能是下列字符：. , * / ^ = > < # ( ) [ ] @ ; “ '5.如果用命令set safe on 指定了编译安全模式，则用户调用fish函数时，函数名前必须加@；6.如果变量不曾赋值，则系统默认为零（整形），如果赋值，其类型由值的类型决定；7.fish函数的调用方法：.可以出项在其他fish函数的单独行中；.可以出现在其他fish函数的表达式中；.出现在flac3d的命令行中；.作为命令set,print,hist的参数。

FLAC3D常用命令

2. range(缩写：ran)通过range功能，可以使命令作用在一定指定范围的目标上；如果一个命令没有使用range来确定范围，则命令对整个模型有效。

例：model elastic range id 1 10 ;指定id为1到10的单元为各向同性弹性本构。

FLAC3D常见命令与使用技巧

FLAC3D常见命令与使用技巧1.基本命令：- !new: 创建一个新的FLAC3D模型- !zone generate: 通过定义边界和离散网格参数生成岩石模型- zone struct properties: 设置岩石模型的物理性质，如弹性模量、泊松比等- zone gridpoint initialize: 初始化模型中的网格节点- zone cmodel assign: 在指定区域内分配一种材料模型- zone face apply stress: 在模型表面施加一定的应力条件- zone solve: 求解模型，得到模型的应力应变分布2.材料模型：- zone cmodel initialize: 初始化材料模型的参数- zone cmodel property assign: 为材料模型分配属性- zone cmodel model assign: 为材料模型指定模型类型3.自定义模型：- zone cmodel new: 创建一个自定义材料模型- zone cmodel property new: 创建一个自定义的材料模型属性- zone cmodel property set: 设置材料模型属性的值- zone cmodel initialize: 初始化材料模型的参数4.辅助命令：- zone gridpoint fix: 固定指定区域的网格节点- zone gridpoint value prescribe: 为指定区域的网格节点预设值- zone gridpoint stress initialize: 初始化模型中的网格节点应力- zone gridpoint velocity assign: 为指定区域的网格节点分配速度- zone gridpoint delete: 删除指定区域的网格节点5.使用技巧：-使用脚本：可通过编写脚本，将一系列命令自动执行，提高工作效率-模型可视化：FLAC3D提供了强大的可视化功能，可以实时观察模型行为-使用参数化模型：可以通过设定参数的方式进行模型设置，便于修改和调整-多个窗口：可以同时打开多个FLAC3D窗口，方便对不同模型进行比较和分析-参考案例：可以参考软件自带的案例，了解不同类型的建模和分析方法这些是FLAC3D的一些常见命令和使用技巧，希望对用户在使用该软件时有所帮助。

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仄里的投影是一条射线.最先该当依照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……7、关于二种fix是可有辨别？fix z range z -0.1 0.1fix x range x -0.1 0.1 any x 1.9 2.1 anyfix y range y -0.1 0.1 any y 1.9 2.1 any战fix z range z -0.1 0.1fix x range x -0.1 0.1fix x range 1.9 2.1fix y range y -0.1 0.1fix y range y 1.9 2.1问：是一般的，any表示并列对接关系，及and的意义.8、initial apply set 的辨别问：initial 初初化下令，初初化估计应力状态，定义资料参数、重力等；apply鸿沟条件节造下令，如施加鸿沟的力、位移等拘束等；set 用于树坐估计的模式，统造估计历程.Ini 施加的初初力正在估计历程中是会变的，普遍体力需要初初化，而apply施加的力、鸿沟条件不会随着估计变更.set下令效率于所有网格(赋值局部变量),ini则初初化指定范畴内网格面、域的偏偏移量.初初应力不但是仅是重力引起的.9、Tension-p shear-p tension-n shear-n表示的含意问：正在估计的循环内里,每个循环中,每个zone 皆依据failure criterion 处于分歧的状态,shear战tension 分别表示果受剪战受推而处于塑性状态.n 表示now,p表示previous,即分别表示正在当前战往日的循环核心处于塑性状态.10、history保存数据，print隐现正在过失的下令流中“”文献中惟有：“Step 1 Y-Displace”而不保存下数值，怎么回事？尔怎么才搞赢得table中的数据？输进'print tabel n'能得到table中的数据（粗确）hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0step 3000hist write 1 file 本面位移.dathist write 2 file 中部单元应力.dathist write 3 file 边沿单元应力.dathist write 1 table 1hist write 2 table 2hist write 3 table 3hist write 2 3 vs 1 file 单压.dat（过失）hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0hist write 1 file 本面位移.dathist write 2 file 中部单元应力.dathist write 3 file 边沿单元应力.dathist write 1 table 1hist write 2 table 2hist write 3 table 3hist write 2 3 vs 1 file 单压.datstep 300011、plot隐现x、y坐标图hist n 1hist gp ydisp 0,0,0hist ax_strhist gp xdisp 1,1,0step 1500plot hist 2 vs -1 ;axial stress vs axial disp.plot hist 3 vs -1 ; circumferential disp. vs axial disp问：那2是不是代表hist gp ydisp 0,0,0，而-1 指得是hist gp ydisp 0,0,0为什么不是1而是背值？1) 正在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正背);2) "-"表示对于其值做"mirror",如果大家使用AUTOCAD的话便明黑"mirror"的含意了,比圆上头的 hist 1 gp ydisp的值是0~100,那么vs -1便形成-100~0.以此类推.3) 要那样搞主假如为了切合岩土工程的习惯需要,比圆桩载荷考查直线,重落直线等皆不是画正在第一象限内.12、关于water.fis的下令注释（树坐天下火位）问：（1）water.fis的功能是树坐天下火位里，正在二维估计中，通过water table下令即可树坐火位线，正在三维估计中由于天下火位是里，果此需要做一些插值运算，以便产死火位里.那个FISH步调包罗三个函数parm,aux1战aux2.parm用去设定火位的边值，即树坐火位的天区，aux1插值估计等值线的面，aux2调用aux1的估计截止，使用FLAC3D 下令产死火位里.（2）最先树坐参数值，而后定义一个table.def parmnptab1=7nprof=10endparmtable 1 (12 0) (18.25 12.5) (25 17) (35 21) (50 24) (75 25) (150 25)爆收table1，包罗7组数据面.3) 正在阐明的二个FISH函数之前，对于table,xtable,ytable做一些总的证明.water.fis步调反复使用了那三个变量，明黑了那三个变量的含意，water.fis那段步调便办理了大半问题了.Itasca 的所有硬件皆使用了相共的数据结构储藏数据，称之为链表结构，类似于关系型数据库中的数据链接要领.如果罕见据结构的知识，大概掌握关系型数据库的SQL 谈话，那三个变量的含意便不易明黑了，只不过那个支配与真真的数据库支配相比，惟有二个数据变量x,y.尔用数据库的术语做阐明：table 便是爆收一个“表”，xtable,ytable的效率便是爆收一个新表(Create)以及对于表中的数据举止与出(Select)、拔出(Insert)、替换(Update)支配.4）FISH函数aux1的效率是正在table 1的前提上爆收出其余个table，table的ID 从11到17，每个table有10个数据集.alfa的估计要领是为切合(fit)直里而采与的坐标面算法.def aux1loop n(1,nptab1)rr=xtable(1,n)zz=ytable(1,n)loop k(1,nprof)case_of kalfa=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 1xtable(n+10,k)=rrytable(n+10,k)=-40.0case 2xtable(n+10,k)=rrytable(n+10,k)=0.0end_caseend_loopend_loopendaux1当真止那个FISH函数后，死成(nptab1+10,nprof)组数据，那些数据为Aux2调用.注意：本段步调爆收的数据是为那个特定的直里边坡创造的，不是一个general purpose的LIB5）FISH函数aux2的效率是使用WATER TABLE face下令爆收火里，把table 1的ytable值动做z坐标，aux1爆收的数据动做x,y坐标.def aux2commandset grav 0 0 -10water dens 1000end_commandloop n(1,nptab1-1)loop k (1,nprof-1)xx1=xtable(n+10,k)yy1=ytable(n+10,k)zz1=ytable(1,n)xx2=xtable(n+10,k+1)yy2=ytable(n+10,k+1)zz2=zz1xx3=xtable(n+11,k+1)yy3=ytable(n+11,k+1)zz3=ytable(1,n+1)xx4=xtable(n+11,k)yy4=ytable(n+11,k)zz4=zz3commandWATER TABLE face xx1,yy1,zz1 xx2,yy2,zz2 xx3,yy3,zz3WATER TABLE face xx1,yy1,zz1 xx3,yy3,zz3 xx4,yy4,zz4end_commandend_loopend_loopendaux213、 FLAC截止数据的导出pl set back white ;树坐背景set plot bitmap size (800,510) ;树坐图片大小plot set caption size 37 ; 树坐图例数字大小 (default=35 范畴10-50) mainwin size 1 1 position 0 0.8 ; command windowplot set window size 0.8 0.8 position 0.1 0.0 ; graphics windowplot set center 100 0 40根据情况变动14、一段fish下令的阐明def aux1 定义函数loop n(1,nptab1) 启初一个循环，从1循环到变量nptab1rr=xtable(1,n) 从table 1 的x列中读第n个数付给变量rrzz=ytable(1,n) 从table 1 的y列中读第n个数付给变量zzloop k(1,nprof) 启初第二个循环，从1循环到变量nprofcase_of k 底下是一个采用结构，利用变量k推断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个值等于变量rrytable(n+10,k)=-40.0 table n+10的y列第k个值等于变量-40case 2 当k＝2的时间真止底下语句xtable(n+10,k)=rr table n+10的x列第k个值等于变量rrytable(n+10,k)=0.0 table n+10的y列第k个值等于变量0end_case 中断采用结构end_loop 中断第二个循环end_loop 中断第一个循环end 函数定义中断aux1 真止什么定义的函数15、有关下令的阐明（1）apply nvel 0.1 plane dip 60 dd 270 range left_bounapply 施加鸿沟条件，plane法背速度为0.1，plane 里的位子：倾角60（正在xoy仄里从-Z轴目标顺时针转动），dd 走背270度，从+y目标绕z轴顺时针转动；（2） ini szz -5e4 grad 0 0 -1e4ini 初初鸿沟条件，szz=-5e4+x*0+y*0+-z*-1e4,注：-5e4是range范畴起面坐目标SZZ值，x,y,z所施加应力网格面的坐标.16、渗流示例;------------------------------参数部分bulk modulus, K 390 MPa 体积模量（土量比较硬）shear modulus, G 280 MPa 剪切模量soil dry density, ρd 1200 kg/m3 土的稀度water density, ρw 1000 kg/m3 火的稀度wall density, ρwal 1500 kg/m3 档土墙的稀度permeability, k 10−12 m2/Pa-s 渗透系数porosity, n 0.3 孔隙率fluid bulk modulus, Kf 2.0 GPa 流体的体积模量;-----------------------以下是下令流;------------------------------------------------------; Excavation in a saturated soil;------------------------------------------------------config fluid ;树坐流体算法; --- geometrical model --- 建坐模型gen zone brick p1 12 0 0 p2 0 12 0 p3 0 0 12 size 12 12 12 rat 1 1 1group soil ;树坐土组group excavate range x 0 4 y 0 4 z 0 5 ;树坐启掘范畴组group wal1 range x 4 5 y 0 5 z 0 7 ; 树坐档土墙 1 组group wal2 range x 0 4 y 4 5 z 0 7 ;树坐档土墙 2 组group wall range group wal1 any group wal2 any ;档土墙1 战2 合并为一个档土墙组; --- fluid flow model --- 流体模型model fl_iso ;各项共性流体prop perm 1e-12 poro 0.3 ;树坐渗透系数战？？比 or ？？率ini fdensity 1e3 ; 初初化火稀度ini fmod 2.0e9 ftens -1e-3 ; 初初化火的体积模量及？？model fl_null range group excavate ; 启掘部分不火流model fl_null range group wall ;档土墙不渗火ini pp 0 grad 0 0 1e4 ;初初化孔隙火压力fix pp range z -0.1 0.1 ;拘束该当是所有模型的底部吧？fix pp range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;启掘部分也给拘束住了，而其余的场合不管; --- mechanical model --- 力教模型model elas ; 弹塑性prop bul 3.9e6 shea 2.8e6 ;树坐体积模量战剪切模量，那二个必须，有一次尔只树坐了体积模量，便错的一塌糊涂model null range group excavate ;力教模型的启掘ini density 1.2e3 ;初初化稀度，该当是土的稀度ini density 1.5e3 range group wall ;初初化墙（档土墙）的稀度，fix x range x -.1 .1 ;施加拘束，传闻好分法拘束0 的时间，必须与一个-0.1<0<0.1,范畴fix x range x 11.9 12.1 ;施加拘束，大概者喊搞鸿沟条件fix y range y -.1 .1 ;施加拘束fix y range y 11.9 12.1fix z range z 11.9 12.1; initial total stresses 出世总应力ini szz 0 grad 0 0 -1.5e4 ;z目标的应力大一些ini sxx 0 grad 0 0 -1.2e4 ;火仄目标的，小一些ini syy 0 grad 0 0 -1.2e4 ;apply nstress 0 grad 0 0 -1.2e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0apply nstress 0 grad 0 0 -1.2e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0apply nstress -7.5e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1; --- setting --- 树坐set gravity 0 0 10 ;树坐重力加速度; --- initial state --- ;初初化状态solve force 1 ; check initial equilibrium 供解仄稳; --- histories --- 记录set hist_rep 40 ;记录步为40hist fltime ;记录时间hist gp pp 0 0 6 ;记录某面的孔隙火压力hist gp xdis 4 0 0 ;火仄位移hist gp xdis 4 0 2 ;记录火仄位移hist gp xdis 4 2 0hist gp zdis 0 0 5 ;记录横背位移hist gp zdis 2 0 5hist gp zdis 4 0 5hist gp zdis 2 2 5hist gp zdis 4 2 5hist gp zdis 4 4 5hist gp zdis 10 0 1hist gp zdis 10 0 2;; --- excavation --- 启掘set fluid off ;关关火; apply pore pressure at walls 正在墙上施加孔隙火压力apply nstress 0 grad 0 0 -1.e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0 apply nstress 0 grad 0 0 -1.e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0 apply nstress -5.e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1solve ;force 1 ;供解save exc1.sav ; 保存;; --- drainage --- 排火apply remove nstress ;撤掉刚刚才的压力def relaxsetup ;定义一个释搁函数的参数step0 = stependrelaxsetup ;调用def relax ;定义一个释搁函数rstep = step - step0if rstep < ncyc thenrelax=1.0-(float(rstep)/float(ncyc))elserelax = 0.0endifend ;中断set ncyc = 1000 ;先赋值，随后调用apply nstress 0 grad 0 0 -1.e4 hist relax &range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0apply nstress 0 grad 0 0 -1.e4 hist relax &range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0apply nstress -5.e4 hist relax &range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1cyc ncycsolve ;供解save exc2.sav ;保存；; --- percolation --- 渗透set fluid on ;树坐火算法fix pp 0 range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;施加孔隙火压力cyc 9000 ;循环save exc3.sav ;保存plot create excav ;建坐一个隐现的视图一下真量皆是隐现真量plot set rot 200 0 195plot set cent 6 6 6plot set dist 39.18plot set magn 0.8plot set plane ori 0 0 0 normal 0 -1 0plot add cont pp out onplot add block group lgra range group wallplot add cont pp int 10000 max 110000 range x 0 4 y 0 4 z 5 7 plot add flow planeplot showRet17、用什么下令不妨位移图中隐现最大位移面的坐标？问：用fish编了一个，算法佳像有面烦，不过脚段不妨真止您不妨参照一下.def get_gp_maxdispgp0_disp = gp_xdisp(gp_head)*gp_xdisp(gp_head)gp0_disp = gp0_disp + gp_ydisp(gp_head)*gp_ydisp(gp_head) gp0_disp = gp0_disp + gp_zdisp(gp_head)*gp_zdisp(gp_head) gp0_disp = sqrt(gp0_disp)p_gp=gp_head;找最大值loop while p_gp # nullgp_disp = gp_xdisp(p_gp)*gp_xdisp(p_gp)gp_disp = gp_disp + gp_ydisp(p_gp)*gp_ydisp(p_gp)gp_disp = gp_disp + gp_zdisp(p_gp)*gp_zdisp(p_gp)gp_disp = sqrt(gp_disp)if gp_disp>gp0_dispgp0_disp=gp_dispendifp_gp = gp_next(p_gp)endloop;找最大值的坐标p_gp=gp_headss=0loop while p_gp # nullgp_disp = gp_xdisp(p_gp)*gp_xdisp(p_gp)gp_disp = gp_disp + gp_ydisp(p_gp)*gp_ydisp(p_gp)gp_disp = gp_disp + gp_zdisp(p_gp)*gp_zdisp(p_gp)gp_disp = sqrt(gp_disp)if gp_disp=gp0_dispss=ss+1xtable(1,ss)=gp_xpos(p_gp)ytable(1,ss)=gp_ypos(p_gp)endifp_gp = gp_next(p_gp)endloopendget_gp_maxdisp18、怎么样隐现估计后模型的变形状态？问：pl sk magfac **(一个搁大的数字)19、怎么样查看剖里上任一面的位移值问：print gp disp range x() y( ) z( )大概 print gp disp range idplot gp disp range x() y( ) z( )大概 plot gp disp range id20、group下令区分斜里问：group groupname colour range plane above dd ** dip ** origin *** plane below dd ** dip ** origin ***dd 战dip便是二个斜里的倾背战倾角,origin是那个斜里通过的某个面,那个下令便不妨group二个斜里间的天区21、ansys导进flac的时间，怎么样区分一下group问：正在ansys中对于分歧group用分歧的资料号（mat）定义，再用ansys2flac之类的步调变换成flac的.dat 文献便可，group自动死成. 22、 fix x y z range x -0.1 0.1 ..........................afix x range x -0.1 0.1................................b请问a式与b式有何辨别？问：a 表示正在x=0 那个仄里是牢固的b 表示正在x=0 那个仄里只x目标不克不迭爆收疏通.23、怎么样了解每个单元的ID？问：用鼠标单打单元的表面，不妨了解单元的ID战坐标.24、怎么样举止切片？问：plo set plane ori (面坐标) norm (法背矢量)plo con sz plane (隐现z目标应力的切片)那二止下令，用去隐现指定里上的位移变更云图（把第二止改改，还不妨隐现应力云图，等其余）；plot con syy y背应力云图plot con szz outline on z背应力云图加隐现网格plot con szz effective on z背灵验应力plot con szz inter 1e4 云图隐现删量（后里的值越小，图上颜色条戴便越细，）plot con szz min -10e4 隐现最小值，图上小于那个数的值皆用共一色表示plot con szz max 10e4 隐现最大值，图上大于那个数的值皆用共一色表示plot fip outside on 隐现力的箭头当前瞅鸿沟，如果是不动里，用fix x range…如果鸿沟戴速率，大概位移用ini…ini dens 1000树坐重力便一定要树坐稀度，不思量重力便不要安排稀度.二者共存的，安排步调的时间一定要记的.set mech force 50solve用SET下令设定了最大不仄衡力的上限值为50N，即：估计的最大不仄衡力一往小于此值，咱们便人为模型达到了仄稳状态.供解历程末止.树坐那个脚段是正在一定粗度央供下，加快估计（值越大，算的量越少，花的时间也越少）.25、怎么样正在步调中举止久停，并可回复估计？问：正在下令中加进pause下令，用continue举止继承.正在咱们分步供解中念得到某一个历程中的截止，不必等到齐供完，还不妨正在分集供解过失的时间便举止改正，而不是等到截止出去. 26、怎么样跳过某个估计步？问：正在估计中按空格键跳过本次估计，自动加进下一步27、Fish是什么物品？Fish是可一定要教？问：是FLAC3D的内置谈话，不妨用去举止参数化模型、完毕下令自己不克不迭举止的功能.Fish不妨不必教，需要的时间查Mannual赢得需要的变量便不妨了.28、怎么样调用一些可选模块？问：config dyn (fluid, creep, cppudm)29、range的下令：①与圆球里的范畴，便是共一个核心，分歧半径之间夹的那一部分，用正在给球里施加力，拘束时用该下令，下令：range annulus center xyz（核心坐标）radius r1 r2（二半径），念要表黑的是球里时，r1r2与值出进小面即可，如radius 2.9 3.1 那时范畴便是半径为3的里.②圆柱的环里施加力大概拘束时也用到该下令.range cylinder end1 xyz end2 xyz radius r(end是表示柱体断里的核心面坐标)range plane above 表示的是z大于0的上部结构range plane blow 表示的是z小于0的下部结构range dd a dip b origin xyz 表示的是走背角度a，倾背角度b的里dd表示，劈里与xoy里相接时的垂线与y轴正背的夹角，顺为正，道心语，便是里的指背目标里是往那个目标的意义.dip里与xoy里的夹角，与z轴背背为正，origin表示里通过xyz面；range normal xyz 法线目标为xyz的里range not 出试用过（范畴中的物品，念念用到圆上与范畴非常佳）. 30、怎么样调用FISH文献？问：新脚调用FISH文献堕落的本果：是：文献名为fishcall.fis.txt(果为正在记事本里那个txt的后缀是默认的）.建改要领是：将fishcall.fis.txt 内里的真量局部复造到新的记事本里，再按保存，保存时间将保存典型改为所有文献(*,*),那时您即可正在文献名栏内输进fishcall.fis，典型也便粗确了.31、问a式与b式有何辨别？问：fix x y z range x -0.1 0.1 ..........................afix x range x -0.1 0.1................................ba 表示正在x=0 那个仄里是牢固的；b 表示正在x=0 那个仄里只x目标不克不迭爆收疏通.32、怎么样查看各个时段不仄衡力的简直数值？问：采与his去记录估计，包罗位移应力等下令his unbal；his gp(zone) zdis range (0 0 0) 大概者id=？；导出数据下令his write n vs m begin 时步 end 时步+ file filename.hisn表示记录的id m表示时步；要导出不仄衡力的简直数值his unbalstep 100000 or solvehis write 1 vs step begin 1 end 1000 file 123.his使用上述下令便不妨查看各个时步下的不仄衡力的简直数值.33、怎么样隐现变形表面线的下令？plo ske magf 10 其中10为搁大系数34、何如查看模型的矢量图？pl grid vel (disp) redplo blo gro不妨瞅到分歧的group的颜色分集（瞅模型的格式）；plo gpfix red sk 瞅模型的鸿沟情况plo fap red sk 瞅模型的体力分集plo con dis (xdis, ydis, zdis) 瞅模型的云图：位移plo con sz (sy, sx, sxy, syz, sxz) 瞅模型的云图：应力plo dis (xdis, ydis, zdis) 瞅模型的矢量图plo info 瞅模型有几单元、节面plot grid 不妨查看网格，plot grid num 不妨查看节面号.plot fix 查看鸿沟条件plot grid mark 查看预先留设的启掘鸿沟pl syy (sxx xdis ydis sig1..) fill；查看模型的云图.注意，戴fill关键字便是查看云图，不戴fill 关键字便是查看等直线图.plo add ske 剖里的重落云图上减少网格plo add dis plane 剖里上的矢量plo add axe 剖里上坐标值35、何如查看INTERFACE？问：pl grid iface id redprint iface id ; id 是interface的编号36、FLAC怎么安排视图？问：ctrl+z 大概shift+z ,而后面打要搁大部位的左上角战左下角便不妨搁大某个部位，如果要缩小，单打屏幕. 37、推断模型中网格面是一个仍旧二个的要领.问：隐现模型网格面的ID，如果有二个ID号，证明是二个.节理里上的网格面该当有二个ID号.38、flac中效率运止时间的果素问：（1）单元资料，结构里之间的刚刚度对于比很大；（2）区分单元尺寸对于比很大，尽管预防少细比大于5：1的细少单元.39、简朴下令阐明：问：1.plot add cont disp out on 常常隐现云图2.plot block group gpnum on 隐现节面编号3.plot inter red ske 瞅交战里的效验.4.del range group 1 not 不简略group疑息5.del range group 2 简略group疑息6.ini z add ** range group ** （name）7.interface 1 face range x（*，*）y（*，*）z（*，*）大概者 interface 1 face range plane orig(x y z 坐标）norm（法线目标）建坐交战里8.water table clear 扫除火位线9.pl set back white 树坐背景10.set plot bitmap size (800,510) 树坐图片大小11.plot set caption size 37 树坐图例数字大小12.plot set center 100 0 40 树坐图片核心13.fix x y z range x -0.1 0.1 表示正在x=0 那个仄里是牢固的14.fix x range x -0.1 0.1 表示那个仄里只x目标不克不迭爆收疏通40、正在牢固里上施加位移时，fix与apply不克不迭共时使用，该当用ini战fix,才搞施加上.41、正在举止三轴压缩时间，怎么样施加围压？问：使用apply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius 1，其中X战Z轴是正在咱们眼睛正对于的里.42、正在举止三轴压缩时间，怎么样对于轴背力举止位移统造？问：使用 ini yvel 1e-7 range y -0.1 0.1 ini yvel 1e-7 range y 1.9 2.1 即可，数值不妨变更.43、正在FLAC3D估计中，怎么样记录步调中的值，如果不念常常记录该怎么弄？问：用 hist wtite n1 vs n2 file.txt，念隔断记录可用set hist_rep n44、念记录交战里上的相关力教参数该怎么办？问：hist inter i nstress /ndis/sstress/sdisp45、交战里树坐参数包罗哪几个？问：coh dial fric kn ks ten.46、怎么样输出单元格的数据？问；expgrid .flacsd47、建模的历程中，大概共一个文献中有重复应用的下令，便不妨用宏下令了.问：宏下令： macro 称呼‘下令的真量’例：macro shiping ’ plot con zdip outline on’用shiping 代替后里的引号里的物品.正在建模中只消有相共的皆不妨用shiping代替，紧缩了书籍写下令的时间，完齐上也能让人瞅的更明黑.注意，您命名的标记最佳跟falc自戴的下令不要重复，处理要领是，称呼后戴数字.如shiping1；如果您念建改宏对于象，改变下令的真量，正在称呼上加单引号便是替换的意义. 例macro ’shiping1‘‘新的下令’.48、group的用法：①group间接加到模型的后里相称于给所有模型命名；②group name range x xyy zz给一个范畴的体命名 .例：正在gen zone brick ------ size ---- group name49、plot中add与set的辨别add本去有的加新的用addset本去不的，树坐用setplot create GravVplot set plane dip=90 dd=0 origin=3,4,0plot set rot 15 0 20plot set center 2.5 4.2 4.0plot add bound behindplot add bcont szz planeplot add axesplot show .咱们念隐现别的图像，又不念把那个弄佳的图像拾掉，则用plot create name建坐一个图像，正在建坐新的图像时，plot create name1（新的名）隐现往日的图像，那么输进如下下令：plot current name.另一种要领，输进：plot show name50、分歧的天区定义分歧的本构模型model mohr range z=0,5prop ClayeyGravel range z=0,5正在定义参数后里加个range范畴，为了预防正在输进下令时总要通过range＋局部坐标域的要领去规定该下令真止的天区，咱们不妨先通过range下令去给分歧的天区命名，而后便不妨用range＋天区名的要领去规定范畴了.那样便预防了屡屡输进下令要规定正在该天区真止时，皆要输进一遍坐标范畴的贫苦.简直怎么样给天区命名，以及日后怎么样用天区名去规定该天区，不妨瞅底下那个例子：range name trench x=-1,1 y=0,4 z=-2,2model null range trench51、左脚规则：建坐模型时，各面的排搁时有顺序的.要领如下：拿出左脚伸直了4指合并指背便是p0到p1的目标，掌心指背便是p0到p2的目标，拇指与4指笔直的目标便是p0到p3的目标.很多人道网格后的云图有戴齿什么的，矫正的要领便是重复用底下的下令Gen ADJust用此下令去微调网格, 使之失集仄滑.该下令不妨连绝使用以删进仄滑效验的.52、如果念瞅模型正在估计历程中的变更情况Set movie avi step 1 file 1.avi；定义动画，动画记录的频次为1时步一频次（估计量大的时间便出需要1时步的记录，不妨10,100时步的记录，缩小电脑的使用）：Movie start；动画启初附件里那个下令“”能将截止准时步道便算截止输出hist n= 5 ：N(Nstep)=5指定了每迭代估计5次记录一次相关的值（如果不指定，则系统默认值为10）；hist unbal 下令的第二止指定记录最大不仄衡力（the maximum unbalanced force）；hist gp zdisp 4,4,8 记录坐标值为（4，4，8）的节面Z目标的位移（那个用的多，查看支敛量时便是那个下令，上头二个部分认为用处不大，也很少用过）.53、应变硬化/硬化模型SS，怎么样采用参数表Pro ftab 1 ctab 2 dtab 3Table 1 0 45 .05 42 .1 40 1 40。

(完整word版)flac命令流

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。

Flac3D命令--完整经典版（2）

Flac3D命令--完整经典版（2）预览说明:预览图片所展示的格式为文档的源格式展示,下载源文件没有水印,内容可编辑和复制实例分析命令：1. X ，Y ，Z 旋转 Shift+ X ，Y ，Z 反向旋转Gen zone ……；model ……；prop ……（材料参数）；set grav 0,0,-9.81（重力加速度）plot add block group red yellow 把在group 中的部分染成红色和黄色plot add axes black 坐标轴线为黑色；print zone stress% K 单元应力结果输出 ini dens 2000 ran z a b （设置初始密度，有时不同层密度不同）；ini ……（设置初始条件）；fix ……（固定界面）set plot jpg ；set plot quality 100 ；plot hard file 1.jpg 图像输出（格式、像素、名称） plot set magf 1.0视图的放大倍数为1.0；plo con szz z 方向应力云图2. ini z add -1 range group one 群one 的所有单元，在z 方向上向下移动1m ；然后合并命令 gen merge 1e-5 range z 0此命令是接触面单元合并成一个整体，1e-5是容差3. (基坑开挖步骤)：Step 1: create initial model state （建立初始模型）Step 2: excavatetrench （开挖隧道）4. group Top range group Base not 定义（群组Base 以外的为）群组Top5. plot blo gro 使得各个群组不同颜色显示6. （两个部分间设置界面；切割法）：gen separate Top 使两部分的接触网格分离为两部分；interface 1 wrap Base Top 在（Base 和Top ）这两部分之间添加接触单元；plot create view_int 显示，并创建标题view_int ；plot add surface 显示表面；plot add interface red 界面颜色红色7. （简单的定义函数及运行函数）new ；def setup 定义函数setup ；numy = 8定义常量numy 为8；depth = 10.0 定义depth 为10；end 结束对函数的定义；setup 运行函数setup8. （隧道生成）上部圆形放射性圆柱及下部块体单元体的建立，然后镜像。

FLAC3D与FLAC常见命令的整理

FLAC/FLAC3D常规问题的整理1．FLAC3D命令的FAQlakewater整理看到其它板块上都有这个FAQ，也就是常见问题问答，今天抽了时间进行了整理，想到了就写下来了，因为看到很多初学者费了很多的时间，但是还是没有将常用的命令掌握，所以这个也可以作为入门的初级教材，使大家能够快速的上手，而不用为了某个小命令到处求助。

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进行很好的充实。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

FLAC3D常见命令与使用技巧1、FLAC3D常见命令：1.FLAC3D是有限元程序吗？答：不是！是有限差分法。

2.最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。

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 idplot gp disp range x() y( ) z( )或plot gp disp range id31.a: fix x y z range x -0.10.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 solvehis write 1 vs step begin 1 end 1000 file.123.his使用上述命令就可以查看各个时步下的不平衡力的具体数值34.如何显示变形轮廓线的命令？plo ske magf 10其中10为放大系数35、怎样开始入门？答：初步学习首先看说明书User's gurid的Section2 Gettingstarted和Section 3 Problemsloving with FLAC36、哪些命令需要掌握？答：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 3DECPFC 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 fix40、怎样查看预先留设的开挖边界？答：plot grid mark。

41、怎么查看模型的云图？答：pl syy (sxx xdis ydis sig1..) fill，注意，带fill关键字就是查看云图，不带fill关键字就是查看等直线图。

42、怎样查看模型的矢量图？答：pl grid vel (disp) red43、怎样查看INTERFACE？答：pl grid iface id redprint 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中本构模型（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 gp num 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 40fix x y z range x -0.1 0.1表示在x=0这个平面是固定的fixed boundaryfix x range x -0.1 0.1表示在x=0这个平面只x方向不能发生运动roller boundary51、fish命令的解释：def aux1定义函数loop n(1,nptab1)开始一个循环，从1循环到变量nptab1rr=xtable(1,n)从table 1的x列中读第n个数付给变量rrzz=ytable(1,n)从table 1的y列中读第n个数付给变量zzloop k(1,nprof)开始第二个循环，从1循环到变量nprofcase_of k下面是一个选择结构，利用变量k判断i Falfa=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个值等于变量rrytable(n+10,k)=-40.0 table n+10的y列第k个值等于变量-40case 2当k＝2的时候执行下面语句xtable(n+10,k)=rr table n+10的x 列第k个值等于变量rrytable(n+10,k)=0.0 table n+10的y列第k个值等于变量0end_case结束选择结构end_loop 结束第二个循环end_loop结束第一个循环end函数定义结束aux1执行什么定义的函数52.在固定面上施加位移时，fix与apply不能同时使用，应该用ini和fix,才能施加上。