《FLAC3D各种命令笔记》由会员分享，可在线阅读，更多相关《FLAC3D各种命令笔记（5页珍藏版）》请在金锄头文库上搜索。
1、hist reset ;清空里面所有的监测点nstress=-40e3 hist ramp 的意思，很简单 nstress = -40e3 * ramp hist 相当于乘法，为了方便控制动荷载的幅值。plotitemsAddstrucktural elementsgeometry cables 显示锚杆的几何形状Pl-add-vector-apply body force 显示施加的力，不过是在求解之后显示Creat 创建一新视图plotitems-add-vetor-apply body force 显示荷载的施加位置Array var(4,3) 二维数组 array var(5)一维 5
2、 个元素的数组 Var(m,n)=urand ;随机 var数组的名字 Pnt=z_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 Loop while p_z # nul
3、l 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 的坐标 Xtable(n,s)对 ID 号为 n 的表的第 s 行，x 列进行赋值，table 表只有 x,y 两列 Pl table 1 line/mark/both 关键词描述曲线的类型，线、标记、二者综合 File_n
4、ame=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= HSel fix x y z (局部坐标系位移速度) / xrot yrot zrot(局部坐标系转动位移) Ini xdisp ydisp zdisp(全局坐标系位移) / xpos(节点全局坐标) / xrdisp(全局坐标系旋转位移)
5、Xv(全局坐标系位移速度) / xrvel (全局坐标系旋转速度)如果边界是 FIX 控制，则自动平衡，不用 apply 控制应力。对地下深部开挖的工作，一般忽略顶部到底部的地心引力的差别，set grav 被省略全局坐标系 x y z 局部坐标系 d（倾向） s（走向） n（法向）plo fap ran group tunnel ;显示施加在组上的体力apply remove 施加新的边界条件前，用以移除旧条件Table n x1 y1 建成一个点的坐标的表, 以供 FLAC 引用，n 为是表号。材料参数附加关键字：Gauss_dev s 以输入数值为平均数，s 为偏差数，高斯正态分布计算G
6、radient gx gy gz 渐变：值=输入值+gx*x+gt*y+gz*zUniform_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 柱体网格 r
7、adcylinder 渐变放射 Radtunnel 六面体隧道渐变放射 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
8、z 0 1 应力 ZZ 方向的分量按（-30e3+0*x+0*y+20e3*z ）渐变plot blo gro 显示分组 plot block group range group dam 只显示 dam 组单元plotitems_block group_modify_勾选 grid point ids_ok 显示节点 ID 编号Hist gp 会自动监测距离输入坐标最近的节点 Hist zone 会自动监测距离输入坐标最近的单元 Plot hist 7 vs 8 Set Hist_rap n 监测间隔，越小，曲线越圆滑 hist Write 7 vs 8 file 6-3hist.txt 输出
9、对比信息文件 print 输出结果 print gp dis id rang 517 any 输出节点 517 的变形值 记录文件： log 文件是 FLAC 输出文件的格式之一 set log on 建立一个名为 flac3d.log 的文件，记录窗口中所以信息，直到 set log off 继续下一个 log 文件，set logfile 6-1.log set log on print gp zdisp rang id 2179 a id 2189 set log off plot con zd / zv 显示 z 方向位移/速度 plot sxx 显示应力等值阴影图（没有 contou
10、r） 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 矢量图显示剖面显示： plot set plane ori 0 1.5 0 norm 0 1 0 设定定一个垂直于 Y 轴的剖面，ori 为切面上一点坐标， norm 为切面的法向（即从 0 1.5 0 引到 0 1 0 连线为 法线） pl add con disp
11、 plane behind 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 6-5.avi 格式设置Movie star 动画开关Solve求解完成后，自动生成文件。 初始地应力场的生成：弹性求解法，分阶段的弹塑性求解法 存在了 null 模型，不能采用 solv
12、e elastic 的求解方法获得初始应力，所以采用分阶段的 弹塑性求解方法。先将 Mohr 模型的凝聚力。俏和伉拉强度护赋俏为无穷人进行求解，保 证在重力作用卜单儿不至十发生屈服，然后冉将 Mohr 模型参数赋值为真实值，再进行求 解。 接触面：分界面应附在小面区域，与密度大的区域。胶合模型：不允许张开和滑移 库伦剪切模型： 抗拉粘结模型：粘结型和滑移型。Bslipe off 不滑动（默认） on 滑动 接触面几何模型建立：移来移去法：首一先将需要建立接触而的两个网格分开建立，然后在一个网格的指定位 置建立接触而，然后把另外一部分网格移到特定的位置，这种方法称为移来移去法。要点 就是把需要建
13、立接触而的网格分开建立，然后在相应的网格上建立接触而，再把没有建立 接触而的网格移动过来，这样做的目的就是使同一坐标点有不同的节点号，从而可以模拟 结构物的滑移、分离等情祝。Gen zon radcyl p0-p11 size a（x） b（y） c（弧度上的单元数）d（放射方向的网格） ratio 1 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 radiu
14、s（半径） ，后一个半径要小于前一个，即两个圆柱的差，得到边界。命令指的整 个圆柱再抠掉后面的 End 圆柱体断面中心点的坐标,两个点确定一条线，再根据半径确定柱体。 Not 就是否认前面的意思 Pl sel geom. 显示创建的 sel 单元 Ini z add -6 rang group pile *pile 组整体在 z 方向平移-6 个单位导来导去法：1、保存完整模型保存为 1.Sav 2、保留需要建接触面的分组，删掉其他所有单元，用到命令：del rang group 2 not，在分组上建立接触面，保存为 2.sav； 3、rest 1.sav & del rang group
15、2 & *导出删除后的模型 expgrid 1.flac3d； 4、rest 2.sav & impgrid 1.flac3d Interface 1 prop kn(法向刚度) ks(切向刚度) c f ks kn 由实验确定，一般取接触面两侧 较小的 E 的 10 倍 *注意语法切割模型法：此方法是在己建好的模型上，把 group 1 与周围组分开，然后在 group 1 上建立接 触面。Gen separate group 1 *分离出 group 1 建面Int 1 wrap 1 2 *int=interface 在 group 1 group 2 之间建立接触单元。 结构单元：ID 结构单元编号 CID 结构构件编号 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 *在半径处建一个无厚度的面
