Flac3D命令--完整经典版 (2)
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&ratio1,1,1,1.5&di m142fillplotsurfgenzonereflectdip0dd90（对xy面做镜像）genzonereflectdip90dd90（对yz面做镜像）（1）dd表示y轴正向顺时针到那条射线的夹角，dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线。
首先应该按照dd的方向大体确定这个面的朝向，dd指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角．（及z轴负方向转过角度）（2）命令：reflectnormalxyzoriginxyz（根据法线和过一点建的面）最常用，至于实际使用过程中有人问道：normal-100与normal100的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面。
(完整word版)FLAC3D命令流
(完整word版)FLAC3D命令流命令流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 。
.configgrid 10，10model elastic(完整word版)FLAC3D命令流group 'User：Soil’ notnullmodel elastic notnull group 'User:Soil'prop density=1500.0 bulk=3E6 shear=1E6 notnull group ’User：Soil’fix x y j 1fix x i 1fix x i 11set gravity=9.81history 999 unbalancedsolvesave state1.sav;—-------—---——----—----------—---—计算第二步;... STATE：STATE2 ..。
(完整版)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 可以查看节点号。
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的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面。
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命令--完整经典版
27. ID号查询通过标题栏查看或者执行【Plotitems】/【1 Block group】/【Modify】命令中
grid_point ID’s 28. 建立接触面单元的各类方法（桩的接触面桩端和桩侧最好采用不同id号的接触面单元） 接触面参数的选取见文件155页9.51。接触面有关的常用命令9.7 162页。 ① 移来移去法：两个模型分开，一个建立接触面单元，然后移动使两个模型合并。
PLOT block group 查看计算模型 Ctrl+G 变成灰色图像
25. 采用 2 种方式观察网格节点速度： ： HIST gpvel<x y z>记录某个节点的速度； PLOT
vel 命令绘制完整的速度矢量场图；力平衡时，由于网格节点力不为零就存在网
格节点速度，仍有可能产生较大位移。 26. 初始地应力场：通常用的是以下三种方法，即弹性求解法、改变参数的弹塑性求解法
定义群组方法！
Gen zone brick size 3 3 3 ； group 2 range x 1 2 y 1 2 z 1 2 ； group 1 range group 2 not
删除其它群组 只保留群组2
del ran group 2 not
Plot blo gro range group dam只显示dam组的单元
6 5 所以 xx 为： xx =-10 10 +10 z ；原点（0,0,0）
14. free x range x -.1 .1 z 6.9 10.1放松x=0 平面上，z=7，10 这一部分在x 方向的约 束（可以在此处产生破坏） 15. 体积模量K和剪切模量G与杨氏模量及泊松比v之间的转换关系如下：
之后重新打开6-2.log
Flac 3D 常用命令总结
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.基本命令：- !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后处理命令集
;0.)实体建模的显示plot set rotation 0 0 45plot block group ;显示块体plot add axes redplot add surface yellowplot show;(1.形成切面的操作);以0为中心取切面plot set plane ori 0 0 0 norm 0 1 0 ;输出szz方向的云图plot con szz plane;添加网格线plot add ske;在切面上添加位移矢量plot add dis plane;在切面上增加坐标系plot add axe;在切面上绘制竖向位移云图plot con zd plane;(2.结果输出命令:print(输出关键词)，gp（节点），zone(单元)，history （监测历史）。
);单元应力结果输出print zone stress;节点变形的信息print gp dis;存储节点信息(可以用记事本打开.log)set log onset 文件名.logprint zone stressprint gp disset log off;(3.变形云图）;生产变形矢量图plot sk dis;输出变形网格plot sk magf 20;(变形倍数);变形后的网格与应力云图结合起来plot con szz ou on magef 20;（4.形成一半的剖面图plot set plane ori 0 0 0 norm 0 1 0plot add surf plane behind yell;在剖面后显示位移云图plot add cont disp plane behind shade on ;;显示单元的id的网格模型plot block group id on;显示最大、最小、中间主应力云图plot con smaxplot con sminplot con smid。
flac3d命令总结
6.加载前的位移清零
ini xdisp=0 ydisp=0 zdisp=0
7.开挖命令
m null range x 20 100 y 20 30 z 20 24.50
8.看切面及记录数据
pl set plane ori 0 0 76.85 norm 0 0 1
12.显示塑性区
plot block state
13.显示部分网格情况（位移上的范围）
plot block state/zadisp/szz range z 3.15 7.65
14.显示部分网格情况（id上的范围）
pl block state/zadisp/szz range id 501 1100
21.出图命令
set plot jpg
set plot quality 40
plot hard file name
22.看已知单元id号岩性参数方法
菜单栏-display-set range-"add"-ID-"在low里输入ID号"-确定-再点display-zones-property-bulk
15.修改间隔
pl con szz interval 2.5e06
16.显示网格
pl grid
17.查看节点号
pl grid num
18.查看模型有多少个单元、节点。
pl info
19.查看边界条件
plot fix
20.如何保存/名
1.建模型命令
new
plot block group
gen zone brick size 60 100 2 p0 0 0 0 p1 120 0 0 p2 0 200 0 p3 0 0 8.40 group 1
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. （隧道生成）上部圆形放射性圆柱及下部块体单元体的建立，然后镜像。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实例分析命令：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. （隧道生成）上部圆形放射性圆柱及下部块体单元体的建立，然后镜像。
9. 模拟模型的材料问题时为什么要去定义某个方向上的初始速度？—10. 渐变应力施加：apply nstress -1e6 gradient 0,0,1e5 range z 3.464,0 plane dip 60dd 270 origin .1 0 0；施加法向应力：apply nstress -1e6 range plane dip 60 dd 270 origin .1 0 011. d ip dd 确定平面位置使用：（纠结）12. p rint gp position range id=14647 输出节点坐标13.a pply sxx -10e6 gradient 0 , 0, 1e5 range z -100 , 0在这个求解方程中，z 为变量，所以xx σ为：65=-1010+10xx z σ⨯⨯ ；原点（0,0,0）14. f ree x range x -.1 .1 z 6.9 10.1放松x=0 平面上，z=7，10 这一部分在x 方向的约束（可以在此处产生破坏）15. 体积模量K 和剪切模量G 与杨氏模量及泊松比v 之间的转换关系如下：16. 一般而言，大多数问题可以采用FLAC 3D 默认的收敛标准（或称相对收敛标准），即当体系最大不平衡力与典型内力的比率R 小于定值10-5；（也可由用户自定义该值，命令：SET mech ratio <value>）17. 所谓体系最大不平衡力，是指每一个计算循环(或称计算时步)中，外力通过网格节点传递分配到体系各节点时，所有节点的外力与内力之差中的最大值；所谓典型内力，则是指计算模型所有网格点力的平均值。
18.要遵循由简单到复杂、由少网格到多网格的思路。
19.混合离散化法的基本原理是通过适当调整四面体应变率张量中的第一不变量，来给予单元更多体积变形方面的灵活性。
20.g en zone ref以z=0平面为对称面镜像生成网格; gen zone ref gen zone ref dip 90以y=0面为对称面镜像生成网格; gen zone ref dip 90 dd 90以x=0平面为对称面镜像生成网格21.球体模型建立：make_sphere函数循环使用。
计算球心至节点距离、获取节点坐标、调整坐标值得到球面至球心距离、完成八分之一球模型。
变换节点内存地址！22.①变形云图：（位移、应力、速度等等）；plot add cont szz out on shade on打开szz的等值线图（outline、effective网格、有效的…显示）；PLOT con szz out on ；PLOT con zd ou on magf 20（块云图bcon不能跟magf放大倍数）；plo con zdisp outline on shade on 打开zdisp的等值线云图②变形矢量图：（变形方向、大小PLOT sk dis scale 0.07其中scale用来改变箭头大小，默认是0.05）；③塑性区分布print zone state都可以用来显示模型的塑性区命令显示那些应力符合屈服准则的区域（或称塑性区）应力正位于屈服面，或说正处于破坏状态时以shear-n 或tension-n 标识；曾进入过屈服状态，但现已经退出以shear-p 或tension-p 标识根据塑性区标识判断破坏机制是否在起作用，可按下述步骤进行：④变量监测（历史跟踪）：Hist gp zdisp 1；Hist gp szz 2输出2和1的hist关系图用：plot his -2 vs 1 如果应力值都为负数，则加-2是为了应力坐标轴为正！默认id=1 2 3 4 …依次（hist id= Gp zdis x y z）；hist write 7 v 8 file 6-3hist.txt使用该命令，程序会在默认文件夹中生成一个名为6-3hist.txt 的文本文件，读者可以打开这个文本文件观察所得的数据结果⑤切片设置：剖面（plot set plane ori 0 1.5 0 norm 0 1 0）、云图命令后面加plane、加网格（ske）加矢量(dis)、坐标系(axe)等修饰效果。
23.输出信息：print zone stress 在命令窗口中会输出模型全部单元的6个方向应力数值print gp dis 该命令会显示模型中全部节点三个方向的变形大小SET log on 在命令行中使用上述命令后，程序在默认目录自动建立一个名为flac3d.log的文件，并开始记录命令窗口中的所有信息，直到用户设置log状态关闭（set log off）为止rest 6-3.savset log onset logfile 6-2.logprint zone stressprint gp dis range id a any id b anyset log off 之后重新打开6-2.log24.最大不平衡力不代表没有了塑性流动，只是不平衡力趋于零，并不完全归零。
记录不平衡力随步数的关系。
hist unbal监测不平衡力，并保留历史记录；save t1.sav保存到文件t1；hist gp zdisp 4,4,8监测网格坐标点（4，4，8）在z方向的位移，并保留历史纪录plot set rot 20 0 30视图的旋转角度为（20，0，30）PLOT block group 查看计算模型 Ctrl+G 变成灰色图像25. 采用2种方式观察网格节点速度：：HIST gpvel<x y z>记录某个节点的速度；PLOTvel 命令绘制完整的速度矢量场图；力平衡时，由于网格节点力不为零就存在网格节点速度，仍有可能产生较大位移。
26. 初始地应力场：通常用的是以下三种方法，即弹性求解法、改变参数的弹塑性求解法以及分阶段弹塑性求解法。
① 由于为弹性求解，在体系达到平衡时，岩、土体中并未有产生屈服的区域model elas ② 更改强度参数的弹塑性求解法生成是指求解过程中始终采用塑性模型model mohr此法与前述弹性求解方法的不同之处在于，计算达到最终平衡时，岩、土体中可能有产生屈服的区域，较合理。
③ （先把基点处总的力算出，然后渐变，一般是递减），高度z=3 密度2000；设置初始应力：ini szz 60e3 grad 0 0 20e3 ran z 0 3;ini sxx 30e3 grad 0 0 10e3 ran z 0 3; ini syy 30e3 grad 0 0 10e3 ran z 0 3;水平为竖向的一半。
如果都在水下，则孔隙水压力为ini pp 30e3 grad 0 0 -10e3 ran z 0 3渐变的过程 饱和密度、干密度、孔隙率以及饱和度间关系式为：=+s d w ns ρρρ（存在问题）？？ 27. I D 号查询通过标题栏查看或者执行【Plotitems 】/【1 Block group 】/【Modify 】命令中grid_point ID ’s28. 建立接触面单元的各类方法（桩的接触面桩端和桩侧最好采用不同id 号的接触面单元）接触面参数的选取见文件155页9.51。
接触面有关的常用命令9.7 162页。
① 移来移去法：两个模型分开，一个建立接触面单元，然后移动使两个模型合并。
ini z add -1 range group one 移动命令。
接触面单元interface () face r ange cylinder end1 (0,0,-4.9) end2 (0,0,-5.1) radius .31 查看施加结构单元那部分命令（施加接触面单元），—柱体两端面圆心及半径。
② 导来导去法：总模型；删掉其他群组，留下一个群组建立接触面单元（文件1）；总模型；删除留下的群组del ran group 2 ，expgrid 导出（文件2）；restore 接触单元文件1，导入文件2，impgrid 文件2。
.interface 1 face range x 1 y 1 2 z 1 2interface 1 face range x 2 y 1 2 z 1 2interface 1 face range x 1 2 y 1 z 1 2interface 1 face range x 1 2 y 1 2 z 1interface 1 face range x 1 2 y 1 2 z 2定义群组方法！Gen zone brick size 3 3 3 ；group 2 range x 1 2 y 1 2 z 1 2 ；group 1 range group 2 not删除其它群组 只保留群组2 del ran group 2 notPlot blo gro range group dam 只显示dam 组的单元③ 切割法：gen separate group1（表示在1组上建立面） ；int 1 wrap group1group2 ， plo int red 红色接触面。