2023年FLAC3D各种命令笔记.pdf

《2023年FLAC3D各种命令笔记.pdf》由会员分享，可在线阅读，更多相关《2023年FLAC3D各种命令笔记.pdf（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、FLAC3D 各种命令笔记 FLAC3D 各种命令笔记 编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（FLAC3D 各种命令笔记）的内容能够给您的工作和学习带来便利。同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快 业绩进步，以下为 FLAC3D 各种命令笔记的全部内容。FLAC3D 各种命令笔记 hist reset；清空里面所有的监测点 nstress=-40e3 his

2、t ramp的意思,很简单 nstress=40e3 ramp hist相当于乘法，为了方便控制动荷载的幅值。plotitems-Add-strucktural elementsgeometry cables 显示锚杆的几何形状 Pl-add-vector-apply body force 显示施加的力,不过是在求解之后显示 Creat 创建一新视图 plotitems-add-vetor-apply body force 显示荷载的施加位置 Array var(4,3）二维数组 array var(5)一维 5 个元素的数组 Var（m,n）=urand;随机 var数组的名字 Pnt=z_

3、next(pnt)查找下一个单元体表 Pl plas 显示塑性区 Array x_b(x，y）z_b（x,y)创建两个 x 行 y 列的表格，可以给其赋值，以便调用.Array就是数组的意思。Hist id=1 gp xd 7 0 0 监测 7 0 0 点处的 x 方向位移，ID=1 材料参数 dilation 剪胀角 FISH语言：Degrade;/180 Pi;Ngp；节点总数 Nzone;单元体总数 单元节点遍历：P_z=zone_head 内容能够给您的工作和学习带来便利同时也真诚的希望收到您的建议和反馈这将是我们进步的源泉前进的动力本文可清空里面所有的监测点各种命令笔记的意思很简单相

4、当于乘法为了方便控制动荷载的幅值显示锚杆的几何形状显示施下一个单元体表显示塑性区创建两个行列的表格可以给其赋值以便调用就是数组的意思监测点处的方向位移材料参数FLAC3D 各种命令笔记 Loop while p_z null P_z=z_next(p_z)End loop Read（ar，n）;将数组 ar 的前 n 个记录写到文件 print zone state 显示塑性区域 编写 fish函数时，尽量采用浮点型，例如 1。0 而不是 1 Z_sig1 z_sig2 z_sig3 fish中的三个主应力，由大到小 Print gp poison rang id n 输出节点 n 的坐标 X

5、table（n，s)对 ID 号为 n 的表的第 s 行，x 列进行赋值,table表只有 x,y 两列 Pl table 1 line/mark/both 关键词描述曲线的类型，线、标记、二者综合 File_name=add_step+string(n)+.sav Gp_near（x，y,z）得到靠近坐标 x,y,z的节点地址 Gp_zdis(p_gp)地址为 p_gp 的节点的 z 向变形 String(n)将变量 n 转换成字符串格式 Dim 看基本网格里面的标准 attach face 常用命令,建模型之后连接面 sxx=syy=（/1)szz szz=H Sel fix x y z（

6、局部坐标系位移速度)/xrot yrot zrot（局部坐标系转动位移）Ini xdisp ydisp zdisp（全局坐标系位移)/xpos(节点全局坐标）/xrdisp(全局坐标系旋转位移)Xv(全局坐标系位移速度）/xrvel(全局坐标系旋转速度）如果边界是 FIX控制，则自动平衡，不用 apply 控制应力.对地下深部开挖的工作,一般忽略顶部到底部的地心引力的差别，set grav被省略 全局坐标系 x y z 局部坐标系 d（倾向）s（走向）n（法向）plo fap ran group tunnel；显示施加在组上的体力 apply remove 施加新的边界条件前，用以移除旧条件

7、Table n x1 y1 x2 y2 x3 y3 内容能够给您的工作和学习带来便利同时也真诚的希望收到您的建议和反馈这将是我们进步的源泉前进的动力本文可清空里面所有的监测点各种命令笔记的意思很简单相当于乘法为了方便控制动荷载的幅值显示锚杆的几何形状显示施下一个单元体表显示塑性区创建两个行列的表格可以给其赋值以便调用就是数组的意思监测点处的方向位移材料参数FLAC3D 各种命令笔记 建成一个点的坐标的表，以供 FLAC 引用，n 为是表号。材料参数附加关键字:Gauss_dev s 以输入数值为平均数，s 为偏差数,高斯正态分布计算 Gradient gx gy gz 渐变：值=输入值+gx*

8、x+gt y+gz*z Uniform_dev s 以输入数值为平均数，s 为偏差数,随机正态分布计算 his unbal set mech force 50 检测最大不平衡力到 50 M n rang group x x1 x2 开挖 Plot grid 显示网格 Ini z add-6 rang group pile *pile组整体在 z 方向平移-6 个单位 建立网格：实体多少米 p0 p1 p2就设多少米,单元数,在百个左右，各方向与实体成比例 Brick 六面块体 redbrick 放射状 cylinder 柱体网格 radcylinder 渐变放射 Radtunnel 六面体隧道

9、渐变放射 cshell 柱形壳体网格 Gen zon redbrick p0p8 size n1 n2 n3 n4 ratio r1 r2 r3 r4 Model mohr 赋予本构模型 Prop bu n shear n coh n fric n 赋予材料性质 弹性只赋值弹性、剪切模量 Fix 固定边界 Ini dens n 初始化密度 App 施加力 Reflect 镜像命令 gen zon reflect norm x y z ori x y z Grad 渐变 ini szz 30e3 grad 0 0 20e3 ran z 0 1 应力 ZZ方向的分量按（-30e3+0*x+0*y+

10、20e3*z）渐变 plot blo gro 显示分组 plot block group range group dam 只显示 dam组单元 plotitems_block group_modify_勾选 grid point ids_ok 显示节点ID 编号 Hist gp 会自动监测距离输入坐标最近的节点 内容能够给您的工作和学习带来便利同时也真诚的希望收到您的建议和反馈这将是我们进步的源泉前进的动力本文可清空里面所有的监测点各种命令笔记的意思很简单相当于乘法为了方便控制动荷载的幅值显示锚杆的几何形状显示施下一个单元体表显示塑性区创建两个行列的表格可以给其赋值以便调用就是数组的意思监测点

11、处的方向位移材料参数FLAC3D 各种命令笔记 Hist zone会自动监测距离输入坐标最近的单元 Plot hist 7 vs 8 Set Hist_rap n 监测间隔，越小，曲线越圆滑 hist Write 7 vs 8 file 6-3hist.txt 输出对比信息文件 print 输出结果 print gp dis id rang 517 any 输出节点 517 的变形值 记录文件：log 文件是 FLAC输出文件的格式之一 set log on 建立一个名为 flac3d。log 的文件，记录窗口中所以信息，直到 set log off 继续下一个 log 文件，set logf

12、ile 61。log set log on print gp zdisp rang id 2179 a id 2189 set log off plot con zd/zv 显示 z 方向位移/速度 plot sxx 显示应力等值阴影图（没有 contour)plot con zd ou on magf 20 附带网格线，变形量放大 20 plot sk 显示网格线 plot sk dis scale 0.07 网格上附加矢量图，scale 默认 0.05 越大，箭头越长 后处理：plot set rot 20 0 30 ;输出前 设置视图角度 plot sk dis scale 矢量图显示

13、剖面显示：plot set plane ori 0 1。5 0 norm 0 1 0 设定定一个垂直于 Y轴的剖面，ori 为切面上一点坐标，norm 为切面的法向（即从 0 1。5 0 引到 0 1 0 内容能够给您的工作和学习带来便利同时也真诚的希望收到您的建议和反馈这将是我们进步的源泉前进的动力本文可清空里面所有的监测点各种命令笔记的意思很简单相当于乘法为了方便控制动荷载的幅值显示锚杆的几何形状显示施下一个单元体表显示塑性区创建两个行列的表格可以给其赋值以便调用就是数组的意思监测点处的方向位移材料参数FLAC3D 各种命令笔记 连线为法线）pl add con disp plane be

14、hind shade on 剖面后面的区域一样显示 plot con zd plane 在设定的剖面上显示 z 方向位移 plot add ske 增加网格线 add命令意思是 在绘图区域增加 plot add dis plane 增加位移矢量 plot add axe 增加坐标系 动画制作：Plot set rot 20 0 30 设置视角 Set movie avi step 1 file 65.avi 格式设置 Movie star 动画开关 Solve 求解完成后，自动生成文件.初始地应力场的生成:弹性求解法，分阶段的弹塑性求解法 存在了 null 模型,不能采用 solve elas

15、tic的求解方法获得初始应力，所以采用分阶段的弹塑性求解方法.先将 Mohr 模型的凝聚力。俏和伉拉强度护赋俏为无穷人进行求解,保证在重力作用卜单儿不至十发生屈服，然后冉将 Mohr 模型参数赋值为真实值，再进行求解.接触面:分界面应附在小面区域,与密度大的区域.胶合模型:不允许张开和滑移 库伦剪切模型：抗拉粘结模型：粘结型和滑移型.Bslipe off 不滑动（默认）on 滑动 接触面几何模型建立：移来移去法：首一先将需要建立接触而的两个网格分开建立,然后在一个网格的指定位置建立接触而，然后把另外一部分网格移到特定的位置，这种方法称为移来移去法。要点就是把需要建立接触而的网格分开建立，然后在

16、相应的网格上建立接触而，再把没有建立接触而的网格内容能够给您的工作和学习带来便利同时也真诚的希望收到您的建议和反馈这将是我们进步的源泉前进的动力本文可清空里面所有的监测点各种命令笔记的意思很简单相当于乘法为了方便控制动荷载的幅值显示锚杆的几何形状显示施下一个单元体表显示塑性区创建两个行列的表格可以给其赋值以便调用就是数组的意思监测点处的方向位移材料参数FLAC3D 各种命令笔记 移动过来，这样做的目的就是使同一坐标点有不同的节点号，从而可以模拟结构物的滑移、分离等情祝.Gen zon radcyl p0-p11 size a（x）b（y）c（弧度上的单元数）d（放射方向的网格）ratio 1

17、1 1 1.15 interface 1 face（在指定范围内所有网格的面上，均创建分界面单元）range cyl end1 0 0 0 end2 0 0 5.1 radius.31&cyl end1 0 0 0 end2 0 0 5。1 radius.29 not radius（半径），后一个半径要小于前一个，即两个圆柱的差，得到边界.命令指的整个圆柱再抠掉后面的 End 圆柱体断面中心点的坐标，两个点确定一条线，再根据半径确定柱体。Not 就是否认前面的意思 Pl sel geom.显示创建的 sel 单元 Ini z add 6 rang group pile pile 组整体在 z

18、方向平移6 个单位 导来导去法：1、保存完整模型保存为 1.Sav 2、保留需要建接触面的分组，删掉其他所有单元,用到命令：del rang group 2 not,在分组上建立接触面，保存为 2。sav;3、rest 1。sav&del rang group 2 *导出删除后的模型 expgrid 1.flac3d；4、rest 2.sav&impgrid 1.flac3d Interface 1 prop kn（法向刚度）ks（切向刚度)c f ks kn 由实验确定，一般取接触面两侧较小的 E的 10 倍 注意语法 切割模型法：此方法是在己建好的模型上，把 group 1与周围组分开，然

19、后在 group 1 上建立接触内容能够给您的工作和学习带来便利同时也真诚的希望收到您的建议和反馈这将是我们进步的源泉前进的动力本文可清空里面所有的监测点各种命令笔记的意思很简单相当于乘法为了方便控制动荷载的幅值显示锚杆的几何形状显示施下一个单元体表显示塑性区创建两个行列的表格可以给其赋值以便调用就是数组的意思监测点处的方向位移材料参数FLAC3D 各种命令笔记 面。Gen separate group 1*分离出 group 1建面 Int 1 wrap 1 2 *int=interface 在 group 1 group 2之间建立接触单元。结构单元：ID 结构单元编号 CID 结构构件编

20、号 node 节点 link 结构连接的编号(结构单元间、结构与网格间）Sel beam beg end nseg Sel beam prop Sel pile id=1 cid=i beg 0 0 0 end 0 0 10 nseg=4 *建立桩单元，设立结构单元 ID和结构构件ID，结构构件 4 个,三个节点 Sel node id=*设置节点 ID Pl sel geo id on nod on scale 0。04 Sel shellsel id=I rang cyl end1 end2 radius not *在半径处建一个无厚度的面 Sel node ini zpos add-25

21、 *壳单元的位移,节点初始位置减少 25 Delete 删除实体单元 sel delete 删除结构单元 plot add sel geom node=off shrinkfac=0 magfac=10 red red 显 示 结 构 单 元 shrinkfac是指结构单元的尺寸，范围是 01,默认值为 0.1 桩单元模拟锚杆支护特性：sel pile prop rockbolt on 可以激活锚杆特性，能够计算支护四周的约束应力、锚杆与周围单元之间的应变软化特性以及锚杆拉断、拉裂等现象。Print sel node?Plot sel node？后处理 其中 fob 表示不平衡力的大小,loc

22、al表示结构节点的局部坐标系 锚杆：三种非全长锚固模拟方法 内容能够给您的工作和学习带来便利同时也真诚的希望收到您的建议和反馈这将是我们进步的源泉前进的动力本文可清空里面所有的监测点各种命令笔记的意思很简单相当于乘法为了方便控制动荷载的幅值显示锚杆的几何形状显示施下一个单元体表显示塑性区创建两个行列的表格可以给其赋值以便调用就是数组的意思监测点处的方向位移材料参数FLAC3D 各种命令笔记 1、通过删除锚索端头，即 Cable 端部的 node 和 zone 之间自动建立的连接，然后在他们之间建立刚性链接来模拟托盘。锚杆（锚索)自由段和锚固段通过设置不同的锚固剂参数来模拟,预紧力加在锚杆自由段

23、。sel del link rang id 1 *删除自动链接 sel link id=100 1 target zone 在 cable 单元与 zone 之间设置新链接 target 目标 sel link attach xdir=rigid ydir=rigid zdir=rigid&六个自由度 rigid 刚性的 xrdir=rigid yrdir=rigid zrdir=rigid rang id=100 设置链接为刚性 sel node fix x y z xr yr zr ran id=1 *固定六个方向的自由度（此句为桩 pile）sel cable id=1 pretensi

24、on 60e3 rang cid 1 10 *自由段施加预紧力 2、将锚杆(锚索）的端头、自由段、锚固段赋不同的属性来模拟非全长锚固预应力锚杆锚索）,端头的锚固参数设为极大来模拟托盘，这样在锚杆锚索)受力时,端头将不会滑动，相当于托盘的作用。预紧力加在锚杆（锚索)自由段。3、借助别的结构单元模拟托盘 内容能够给您的工作和学习带来便利同时也真诚的希望收到您的建议和反馈这将是我们进步的源泉前进的动力本文可清空里面所有的监测点各种命令笔记的意思很简单相当于乘法为了方便控制动荷载的幅值显示锚杆的几何形状显示施下一个单元体表显示塑性区创建两个行列的表格可以给其赋值以便调用就是数组的意思监测点处的方向位移材料参数