flac3d入门指南.ppt

《flac3d入门指南.ppt》由会员分享，可在线阅读，更多相关《flac3d入门指南.ppt（49页珍藏版）》请在优知文库上搜索。

1、学习内容 1、基础知识 2、建模、划分网格 3、本构模型及参数 4、边界条件、初始条件及加载 5、初始地应力的生成方法及初 始平衡求解 6、求解及结果输出1、基础知识1.1 通用菜单(界面，菜单，术语)1.2一个简单分析计算样例u 建模、划分网格u 定义材料模型和参数u 初始条件及边界条件u 求解u 结果分析 2、建模、划分网格 在FLAC3D程序中建立计算网格主要采用gen命令，该命令可生成点(point)、面(surface)和单元(zone)。由于点和面在三维计算程序中应用相对较少，此处主要介绍单元(zone)的生成和组合方法。主要命令： brick gen zone +关键词cshel

2、l- radcyl 对于任何形状的单元体，其建立单元模型时关键点的描述需遵从一定的顺序，如建立一规则的六面体，其命令如下：gen zone brick p0 * * * p1 * * * p2 * * * & p3 * * * size * * * rat * * * 当然，对于不规则的六面体我们需要指定全所有关键点的坐标，即从p0至p7所有点的坐标。同时，关键词size还可配合ratio来进行运用，使得各单元间的长度按照一定的比率逐渐增大或减小。 gen zone radcyl size 5 10 6 12 & ratio 1 1 1 1.2 建立比较复杂的计算模型，即通过生成这些基本的单元

3、来进行“拼凑”。如建立一个马蹄形断面的隧道 gen zone radcyl size 5 10 6 12 & ratio 1 1 1 1.2 p0 0,0,0 p1 100, 0,0 & p2 0,200,0 p3 0,0,100 gen zone radtun size 5 10 5 12 & ratio 1 1 1 1.2 p0 0,0,0 p1 0,0, -100 & p2 0,200,0 p3 100,0,0 这是通过radtun和 radcyl来组合生成所需要的模型。它们两者的生成关键点的描述存在较大的区别。 对于这两种基本的网格，其公共面上的关键点的对应关系更需校核好，否则将出现杂

4、乱错误的网格。 对此马蹄形隧道，其公共面处，p0 p0,p1p3, p2p2,p4p5 , p8p9,p10 p11对于对称的模型也可以采用镜像命令： gen zone reflect norm -1 0 0 & origin 0,0,0网格单元间的连接 采用FLAC3D进行计算，所建立的模型需是一个连续的整体，否则计算结果将出现较大的误差甚至无法进行计算。 对于在建立模型时，各关键点的坐标是准确无误输入且各公共面的网格数和大小均完全一致的模型，无需进行任何操作，模型即自动完成相互间的连接。此处所讨论的是公共面上网格的大小和划分的份数不一致或者两相邻面间存在间隙的问题。对于公共面上网格的大小和

5、划分的份数不一致，主要采用attach命令来进行连接，而对两相邻面间存在间隙的问题，则采用gen merge 命令来进行连接的操作。需要指出：只有网格数之间存在整数倍的关系，才能采用attach命令来进行连接，否则连接会失败或仅能连接部分节点。1.2634e61.2822e6Gen 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尺寸和份数要求为？注意：*为了避免边界效应的影响，模型的外边界尺寸最小应为应为地下开挖洞室直径的35倍。*为了保证

6、计算结果的准确性，洞室周围的网格划分应密集一些，但也不宜过细，否则会影响计算速度。同时，单元体的边长比值要控制在一定的范围内，尽量避免比较狭长的单元体出现。*为了保证网格的连续性，应保证相邻边界节点的匹配，例如，相邻网格有相同的单元体数和一致的单元体几何变化率。3、本构模型及参数 FLAC3D中为岩土工程问题的求解开发了特有的本构模型，总共包含了11种力学材料模型： 1. 开挖模型null; 2. 3个弹性模型（各向同性，横观各向同性和正 交各向异性弹性模型）; 3. 7个塑性模型（Drucker-Prager模型、Mohr- Coulomb模型、应变硬化/软化模型、霍克-布朗模型、遍布节理模

7、型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模型）. 对所计算的模型选择合理的本构关系，采用Model命令来进行指定，如定义计算模型为摩尔-库仑模型，则命令流为： Model mohr（可简写为：m m） 确定了合理的本构计算模型后，即赋予模型相应的材料属性，主要采用roperity命令执行，同时对于不同的本构模型，其材料参数各异，如对于摩尔-库仑本构模型：prop bulk * shear * fric * coh * ten * dilation* 对于体积模量和剪切模量，其和弹性模量之间存在一转化公式：或 对于材料的密度(干密度、湿密度等)则采用initial命令来设置，即：

8、ini density * (range -) 材料若考虑密度，则必须设置重力加速度，重力加速度的设置采用set命令，即： set gravity 0 0 10 若重力方向为沿Z轴正向，则为10，若沿Z轴负向，则为-10，若沿Y轴负向，则应设置为： set gravity 0 10 0注意:材料的本构模型必须先定义，严格按照本构模型所要求的材料参数关键字赋予材料参数，否则会导致命令流出错。若本构模型需要的材料参数没有指定值时，系统使用默认值，默认值为0。数值模拟的计算参数与现场勘测值并不完全一致，需要进行适当的调试，使其计算结果与真实情况更为符合。考虑材料密度时，施加重力加速度要注意坐标轴的方

9、向。4、边界条件及初始条件 在FLAC3D中，包含多种边界条件，边界方位可以任意变化，边界条件可以是速度边界、应力边界，单元内部可以给定初始应力，节点可以给定初始位移、速度等，还可以给定地下水位以计算有效应力等。这众多的边界条件主要通过apply或fix命令来进行设置。而初始条件则主要通过initial命令来执行，对所提的这两个命令必须严格区分并了解其差异。通常我们所计算的模型均采用力学边界，初始条件也基本是初始地应力的输入，对此两种不同的力，其设置存在差别，同时在计算过程中，该二者的变化情况也各不相同。 Apply 命令施加的边界力，其在计算过程中保持不变，且任意时刻均存在，除非再次执行ap

10、ply命令进行修改或网格移除，如地表超载等；而ini命令施加的初始地应力，在计算过程中，随着计算模型位移的产生，其内部的应力将不断进行调整和重分配。其基本格式为： apply szz * range z 60.1 59.9(施加在边界面上) ini szz * range z 0 60(施加于体上) 同时可配合gradient施加线形变化的荷载，如： apply sxx -10e6 gradient 0,0,1e5 apply sxx -10e6 gradient 0,0,1e5 range z -100,0 上述命令流所施加的荷载值为： 其具体求解方法为： 如对于z=-50处，sxx=-10

11、e6+1e5*(-50)=-15e6 对于存在多个不同变化梯度的初始应力的设置，必须保证在梯度变化处的应力值相等，若出现应力的跳跃将会产生较大的误差，其计算结果将不可取。 对于位移边界，通常情况下是设置固定边界，即约束各边界在法线方向不发生任何位移，采用fix命令进行设置，例如： fix x range x 49.9 50.1 (平面x=50处) 若在计算过程中，需对原先约束的平面接触约束，则采用free命令进行设置。 对于速度边界，则相应可理解为间接位移边界，设置的边界节点的速度，将随着计算时步产生相应的位移，对于日常生产中的伺服控制以及一些常规的试验，如单轴三轴实验等可采用速度边界来进行较

12、准确的位移控制和定位。 速度边界设置主要采用apply命令进行设置，相应的可设置的速度变量为xvel、yvel、zvel，例如： apply xvel 2e-7 range x 0.1 0.1 其表示将x=0面上所有的节点速度均设置为2e-7，若随着计算，需移除初始边界，则采用apply remove 命令进行删除操作。 在进行速度边界设置时，速度的取值是关键，若已知固定时间内的位移，速度的求解为： 速度=位移/计算时步 再者，进行速度和位移清零时，则可采用ini命令，即：ini xvel 0 yvel 0 zvel 0 xdisp 0 & ydisp 0 zdisp 0注意： 要特别注意fi

13、x x y z range z 与fix z range z 之间的差别，前者是x、y、z三方向的位移全被限制，后者只限制了z方向的位移。 范围要准确。 ini命令赋予梯度应力要注意初始值的计算。5、初始地应力的生成方法及初始平衡求解 求解分两种类型，一类是初始平衡求解，另一类是进行了相应操作后(如开挖)的计算求解。对于任何三维数值计算模型，其初始平衡计算必不可少，待模型建立完成后，设置了本构和材料参数以及相应的边界条件和初始应力后，即可进行初始平衡求解，其目的是为了建立和原始地层所处环境相接近的计算模型，以此可进行相应的后续操作。 初始平衡的求解，最关心的问题是初始地应力的获得，这里介绍两种

14、常用的方法去获得初始地应力(基于无现场地应力测试资料)。方法一：重力自动分配法 该法仅指定材料的密度和重力加速度，边界进行合理约束的情况下，由程序自身计算和分配应力，以得到连续的三向地应力场，采用命令solve elastic执行。其缺点是：求解速度慢，检错性弱，前期变形大。方法二：手动设置法 该法即手动进行初始地应力的设置，若有充足的勘察资料，即可按实际资料进行详细设置，若无现场资料，则垂直应力按土层自重设置，水平应力与垂直应力间考虑侧压力系数，可按v/1-v(v为泊松比)进行选择，采用solve命令。优点：求解速度快，检错性强，前期变形合理。初始地应力的生成为什么要单独列出？分析过程出现的

15、很多问题都与初始应力是否合理有关手册中的例子五花八门生成方法弹性求解更改强度参数的弹塑性求解设置初始应力的弹塑性求解存在水压力的初始应力生成水下建筑的初始应力生成弹性求解：gen zon bri size 1 1 2m elasprop bulk 3e7 shear 1e7fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solveStep = 162s sz = -40e3s sx = -21.54e3更改强度参数的弹塑性求解：gen zon bri si

16、ze 1 1 2model mohrprop bulk 3e7 shear 1e7 c 1e10 f 15 ten 1e10fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solveprop bulk 3e7 shear 1e7 c 10e3 f 15 ten 0solveStep = 163s sz = -40e3s sx = -21.54e3设置初始应力的弹塑性求解：gen zon bri size 1 1 2model mohrprop bulk 3e7 shear 1e7 c 10e3 f 15 ten 0fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000ini szz -40e3 grad 0 0 20e3 ran z 0 2ini syy -20e3 grad 0 0 10e3 ran z 0 2ini sxx -20e3 grad