flacd入门指南.pptVIP

第二部分： FLAC3D软件学习+实例操作 学习内容 1、基础知识 2、建模、划分网格 3、本构模型及参数 4、边界条件、初始条件及加载 5、初始地应力的生成方法及初 始平衡求解 6、求解及结果输出 1、基础知识 1.1 通用菜单(界面，菜单，术语) 1.2一个简单分析计算样例 建模、划分网格 定义材料模型和参数 初始条件及边界条件 求解 结果分析 2、建模、划分网格 在FLAC3D程序中建立计算网格主要采用gen命令，该命令可生成点(point)、面(surface)和单元(zone)。由于点和面在三维计算程序中应用相对较少，此处主要介绍单元(zone)的生成和组合方法。主要命令： 　　　　　　　　　　 brick 　 gen zone +关键词　cshell　------ 　　　　　　　　　　 radcyl 对于任何形状的单元体，其建立单元模型时关键点的描述需遵从一定的顺序，如建立一规则的六面体，其命令如下：　 当然，对于不规则的六面体我们需要指定全所有关键点的坐标，即从p0至p7所有点的坐标。 同时，关键词size还可配合ratio来进行运用，使得各单元间的长度按照一定的比率逐渐增大或减小。 gen zone radcyl size 5 10 6 12 ratio 1 1 1 1.2 建立比较复杂的计算模型，即通过生成这些基本的单元来进行“拼凑”。如建立一个马蹄形断面的隧道 对于对称的模型也可以采用镜像命令： gen zone reflect norm -1 0 0 origin 0,0,0 网格单元间的连接 采用FLAC3D进行计算，所建立的模型需是一个连续的整体，否则计算结果将出现较大的误差甚至无法进行计算。 对于在建立模型时，各关键点的坐标是准确无误输入且各公共面的网格数和大小均完全一致的模型，无需进行任何操作，模型即自动完成相互间的连接。此处所讨论的是公共面上网格的大小和划分的份数不一致或者两相邻面间存在间隙的问题。对于公共面上网格的大小和划分的份数不一致，主要采用attach命令来进行连接，而对两相邻面间存在间隙的问题，则采用gen merge 命令来进行连接的操作。 需要指出：只有网格数之间存在整数倍的关系，才能采用attach命令来进行连接，否则连接会失败或仅能连接部分节点。 Gen zone brick size 4 4 4Gen zone brick size 4 4 4 p0 4.1 0 0 p1 8 0 0 p2 4.1 4 0 p3 4.1 0 4Gen merge 0.1 注意： 为了避免边界效应的影响，模型的外边界尺寸最小应为应为地下开挖洞室直径的3～5倍。 为了保证计算结果的准确性，洞室周围的网格划分应密集一些，但也不宜过细，否则会影响计算速度。同时，单元体的边长比值要控制在一定的范围内，尽量避免比较狭长的单元体出现。 为了保证网格的连续性，应保证相邻边界节点的匹配，例如，相邻网格有相同的单元体数和一致的单元体几何变化率。 3、本构模型及参数 FLAC3D中为岩土工程问题的求解开发了特有的本构模型，总共包含了11种力学材料模型：  1. 开挖模型null;  2. 3个弹性模型（各向同性，横观各向同性和正 交各向异性弹性模型）;  3. 7个塑性模型（Drucker-Prager模型、Mohr- Coulomb模型、应变硬化/软化模型、霍克-布朗模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模型）. 对所计算的模型选择合理的本构关系，采用Model命令来进行指定，如定义计算模型为摩尔-库仑模型，则命令流为： Model mohr（可简写为：m m） 确定了合理的本构计算模型后，即赋予模型相应的材料属性，主要采用roperity命令执行，同时对于不同的本构模型，其材料参数各异，如对于摩尔-库仑本构模型： prop bulk * shear * fric * coh * ten * dilation* 对于体积模量和剪切模量，其和弹性模量之间存在一转化公式： 对于材料的密度(干密度、湿密度等)则采用initial命令来设