《flac3d命令--完整经典版》由会员分享，可在线阅读，更多相关《flac3d命令--完整经典版（5页珍藏版）》请在金锄头文库上搜索。
1、实例分析命令：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 ；
2、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: excavate trench（开挖隧道）4. group Top range grou
3、p 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；num
4、y = 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 60 dd 270 origin .1 0 0；施加法向应力：apply nstress -1e6 range plane dip 60 dd 270 origin .1 0
5、 011.dip dd 确定平面位置使用：（纠结）12.print gp position range id=14647 输出节点坐标 13.apply sxx -10e6 gradient 0 , 0, 1e5 range z -100 , 0 在这个求解方程中，z 为变量，所以 为： ；原点（0,0,0）x65=-1+x14.free x range x -.1 .1 z 6.9 10.1放松x=0 平面上，z=7 ，10 这一部分在x 方向的约束（可以在此处产生破坏）15.体积模量K和剪切模量G与杨氏模量及泊松比v之间的转换关系如下：=3(1-2)E(+v)16.一般而言，大多数问题可以
6、采用 FLAC3D 默认的收敛标准（或称相对收敛标准） ，即当体系最大不平衡力与典型内力的比率 R 小于定值 10-5；（也可由用户自定义该值，命令：SET mech ratio ）17.所谓体系最大不平衡力，是指每一个计算循环(或称计算时步) 中，外力通过网格节点传递分配到体系各节点时，所有节点的外力与内力之差中的最大值；所谓典型内力，则是指计算模型所有网格点力的平均值。18.要遵循由简单到复杂、由少网格到多网格的思路。19.混合离散化法的基本原理是通过适当调整四面体应变率张量中的第一不变量，来给予单元更多体积变形方面的灵活性。20.gen zone ref 以 z=0 平面为对称面镜像生成
7、网格 ; 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 ；PLO
8、T 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 标识根
9、据塑性区标识判断破坏机制是否在起作用，可按下述步骤进行：变量监测（历史跟踪）：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
10、 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.sav set log on set logfile 6-2.log print zone
11、stress print 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
12、种方式观察网格节点速度：HIST gpvel记录某个节点的速度；PLOT vel 命令绘制完整的速度矢量场图；力平衡时，由于网格节点力不为零就存在网格节点速度，仍有可能产生较大位移。26.初始地应力场：通常用的是以下三种方法，即弹性求解法、改变参数的弹塑性求解法以及分阶段弹塑性求解法。 由于为弹性求解，在体系达到平衡时，岩、土体中并未有产生屈服的区域 model elas 更改强度参数的弹塑性求解法生成是指求解过程中始终采用塑性模型 model mohr 此法与前述弹性求解方法的不同之处在于，计算达到最终平衡时，岩、土体中可能有产生屈服的区域，较合理。 （先把基点处总的力算出，然后渐变，一般是
13、递减） ，高度 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 渐变的过程饱和密度、干密度、孔隙率以及饱和度间关系式为： （存在问题）？=+sdwn？27.ID号查询通过 标题栏查看 或者执行【Plotitems 】/【1 Block group】/【Modify
14、 】命令中grid_point IDs 28.建立接触面单元的各类方法（桩的接触面桩端和桩侧最好采用不同id号的接触面单元）接触面参数的选取见文件155页9.51。接触面有关的常用命令9.7 162页。 移来移去法：两个模型分开，一个建立接触面单元，然后移动使两个模型合并。ini z add -1 range group one 移动命令。接触面单元 interface () face range cylinder end1 (0,0,-4.9) end2 (0,0,-5.1) radius .31 查看施加结构单元那部分命令（施加接触面单元）， 柱体两端面圆心及半径。 导来导去法：总模型；删
15、掉其他群组，留下一个群组建立接触面单元（文件1）；总模型；删除留下的群组 del ran group 2 ， expgrid 导出（文件 2）；restore 接触单元文件 1，导入文件 2，impgrid 文件 2。. interface 1 face range x 1 y 1 2 z 1 2 interface 1 face range x 2 y 1 2 z 1 2 interface 1 face range x 1 2 y 1 z 1 2 interface 1 face range x 1 2 y 1 2 z 1 interface 1 face range x 1 2 y 1
16、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 group1 group2 ， plo int red 红色接触面。29.range name=Bin group Bin range name=Material group Material range name= 组 group 组30.建立各类结构单元： id on node on scale （）改变节点、id大小。两单元共用一节点时，需要连接时用 node-
