FLAC3D数值模拟基础.ppt

1、FLAC3DFLAC3D数值模拟基础数值模拟基础FLAC3D2主要内容主要内容n FLAC3D软件简介软件简介n1、基坑开挖、基坑开挖n2、浅基础的稳定性分析、浅基础的稳定性分析n3、网格的合并联结、网格的合并联结n4、界面的生成、界面的生成n5、隧道的生成、隧道的生成n6、模型材料问题、模型材料问题n7、Interface 合并（联结）问题合并（联结）问题n8、初始条件问题、初始条件问题n9、破坏问题分析破坏问题分析n10、综合实例、综合实例-煤巷应力分析煤巷应力分析FLAC3D3主要内容主要内容n FLAC3D软件简介软件简介n1、基坑开挖、基坑开挖n2、浅基础的稳定性分析、浅基础的稳定性

2、分析n3、网格的合并联结、网格的合并联结n4、界面的生成、界面的生成n5、隧道的生成、隧道的生成n6、模型材料问题、模型材料问题n7、Interface 合并（联结）问题合并（联结）问题n8、初始条件问题、初始条件问题n9、破坏问题分析破坏问题分析n10、综合实例、综合实例-煤巷应力分析煤巷应力分析FLAC3D4FLAC3D简介简介nFast Lagrangian Analysis of Continuan美国美国Itasca（依泰斯卡）咨询公司开发（依泰斯卡）咨询公司开发2D程序程序(1986)n1990年代初引入中国年代初引入中国n有限差分法有限差分法(FDM)nDOS版版2.0 2.1

3、3.0nItasca其他软件其他软件FLAC3D5FLAC3D简介简介n应用：应用：岩土力学/岩石力学分析，例矿体滑坡、煤矿开采沉陷预测、水利枢纽岩体稳定性分析、采矿巷道稳定性研究等岩土工程、采矿工程、水利工程、地质工程 n特色：特色：大应变模拟完全动态运动方程使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍 显示求解具有较快的非线性求解速度FLAC3D61 承受荷载能力与变形分析：用于边坡稳定和基础设计 2 渐进破坏与坍塌反演：用于硬岩采矿和隧道设计 3 断层构造的影响研究：用于采矿设计 4 施加于地质体锚索支护所提供的支护力研究：岩锚和土钉的设计 5 排水和不排水加载条件下全饱和

4、流体流动和孔隙压力扩散研究：挡土墙结构的地下水流动和土体固结研究 6 粘性材料的蠕变特性：用于碳酸钾盐矿设计 7 陡滑面地质结构的动态加载：用于地震工程和矿山岩爆研究 8 爆炸荷载和振动的动态响应：用于隧道开挖和采矿活动 9 结构的地震感应：用于土坝设计 10 由于温度诱发荷载所导致的变形和结构的不稳定：高辐射废料地下埋藏的性能评价 12 大变形材料分析：用于研究粮仓谷物流动及井巷和矿洞中材料的总体流动 FLAC3D简介简介FLAC3D7基本原理基本原理开开开开 始始始始生成网格并调整网格的形状；生成网格并调整网格的形状；生成网格并调整网格的形状；生成网格并调整网格的形状；持续的运动和连续的物

5、质属性；持续的运动和连续的物质属性；持续的运动和连续的物质属性；持续的运动和连续的物质属性；特定的边界条件和初始条件。特定的边界条件和初始条件。特定的边界条件和初始条件。特定的边界条件和初始条件。到达平衡状态到达平衡状态到达平衡状态到达平衡状态效果是否符合要求效果是否符合要求效果是否符合要求效果是否符合要求作如下改变：作如下改变：作如下改变：作如下改变：开挖模型的物质属性开挖模型的物质属性开挖模型的物质属性开挖模型的物质属性改变边界条件改变边界条件改变边界条件改变边界条件实施求解实施求解实施求解实施求解效果是否符合要求效果是否符合要求效果是否符合要求效果是否符合要求参数是否调整参数是否调整参数

6、是否调整参数是否调整结束结束结束结束网格确定问题的几何尺寸；持续的运动和连续的物质属性决定了模型的扰动（如由于开挖引起的变形）形式；边界条件和初始条件确定了模型的初始状态（没有引起扰动或变形的状态）。建立FLAC计算模型，必须进行以下三个方面的工作：1.有限差分网格 2.本构特性与材料性质 3.边界条件与初始条件 完成上述工作后，可以获得模型的初始平衡状态，也就是模拟开挖前的原岩应力状态。然后，进行工程开挖或改变边界条件来进行工程的响应分析。FLAC3D8基本原理基本原理n有限差分法有限差分法nLagrangian网格网格n空间混合离散技术空间混合离散技术nLagrangian格式动量平衡方程

7、格式动量平衡方程nFLAC3D的求解过程的求解过程nFLAC3D的本构模型的本构模型FLAC3D9有限差分法有限差分法n古老的方法古老的方法(上世纪上世纪40年代年代)n用差分格式转化控制方程中的微商格式用差分格式转化控制方程中的微商格式n流体力学；土工渗流问题；固结流体力学；土工渗流问题；固结nFDM&FEM的混合求解的混合求解nFDM的新进展的新进展FLAC3D10Lagrangian网格网格n源自流体力学中的拉格朗日法源自流体力学中的拉格朗日法跟踪流体质点的运动状态跟踪固体力学中结点，按时步用Lagrangian法研究网格节点的运动n节点和单元随材料移动，边界和接节点和单元随材料移动，边

8、界和接触面与单元的边缘一致触面与单元的边缘一致n固体力学大变形理论固体力学大变形理论法国数学家、物理学家拉格朗日 FLAC3D11FLAC3D的求解过程的求解过程平衡方程(动量方程)应力应变关系(本构模型)Gauss定律单元积分应变率速度节点力新的应力对所有的网格节点对所有单元FLAC3D12FLAC3D中的本构模型中的本构模型n开挖模型开挖模型null n3个弹性模型个弹性模型各向同性弹性横观各向同性弹性正交各向同性弹性 n8个塑性模型（个塑性模型（Drucker-Prager模型、模型、Morh-Coulomb模型、应变硬化模型、应变硬化/软化模型、遍布节理模软化模型、遍布节理模型、双线性

9、应变硬化型、双线性应变硬化/软化遍布节理模型、修正剑软化遍布节理模型、修正剑桥模型和桥模型和胡克布朗模型胡克布朗模型）FLAC3D13FLAC3D中的本构模型中的本构模型FLAC3D14FLAC3D的的前后处理前后处理n术语术语水平边界 压力结构线模型边界格网点格网区域水 平边 界压力滚 动 底端边界内部开挖边界区域（区域（Zone）有限差分划分的带在几何上是最小的区域，在在这个区域里的每一个现象的变化，如应力应变都可以估计出，。各种形状的多面体（立方体、楔形、锥体、四面体等）可用来构造模型并可用plot显示出来。每一个多面体可能有一套或两套表层设置，这由5个四面体组成。默认的情况下，两个表层

10、设置用在对计算精度要求高的情况下，区域的另外一种叫法是要素。栅格点（栅格点（GridPoint）栅格点是有限差分单元的角点。一个多面体可能有5个、6个、7个或8个网格点，主要取决于多面体的形状。给定每个节点的x，y和z值这样就具体确定了有限差分单元，。其他叫法有：节点，交点。有限差分栅格（有限差分栅格（Finite Difference Grid）有限差分网格是研究区域中一个或多个通过物理边界连接的有限差分单元的集合。另一个叫法是网格，有限差分网格也可以标识出模型中每个状态的存储位置，FLAC3D所生成的矢量都保存在节点上（如：受力、速度、位移）。标量和张量保存在单元的中心（如应力、材料属性）

11、。FLAC3D15FLAC3D的的前后处理前后处理n术语术语水平边界 压力结构线模型边界格网点格网区域水 平边 界压力滚 动 底端边界内部开挖边界模型边界（模型边界（Model Boundary）即有限差分网格的外围，内部边界也同样是模型边界（如网格中的空洞）。边界条件（边界条件（Boundary Condition）即模型边界的约束条件或控制条件的给定（如：限制位移、渗透条件、绝热条件）。初始条件（初始条件（Initial Conditions）即在对模型加载或开挖等作用前的各种参数状态。基本模型（基本模型（Constitutive Model)基本模型（材料模型）即规定了FLAC3D模型中

12、某一区域的变形或强度效应，可用大量基本模型去近视地质材料，可以单独定义FLAC3D模型中的基本模型和材料模型。空单元（空单元（Null Zone）空单元表示此区域为空（就象没有材料一样）。次级网格（次级网格（SUB-GRID）有限差分网格可由次级网格组成，它可用来在模型中创建不同形状的区域，次级网格是分别生成，可进行合并和连接。FLAC3D16FLAC3D的的前后处理前后处理n术语术语水平边界 压力结构线模型边界格网点格网区域水 平边 界压力滚 动 底端边界内部开挖边界附属接触面附属接触面ATTACHED FACES附属栅格面是由被划分的次级栅格组成的网格接触或合并的面栅格面，接触面必须是共面

13、或接触，每个面的节点不一定一样，不同总密度的次级网格可以接触。接触面接触面INTERFACE-即次级网格在计算过程中可以分开（滑动，开裂）的两部分之间的面，可表示不连续的物理特征，如，断层，节理面或材料性质突变的临界面。范围范围INTERFACE-范围是对FLAC3D模型空间值的一个描述，可给定一个命令的作用范围，即使模型发生运动-不影响模型中的区域和节点位置，一个范围或范围确定的单元也不发生改变。范围或范围内的单元所包含的区域，也可以用一个单元的ID号来表示，它与区域，节点，或结构单元密切联系。集合集合GROUP-在FLAC3D模型中他们有共同的名称，由于限定具体命令的对象，如model命令

14、对某一集合设置为某种材料，任何命令加于集合名称也就相当于作用于这一集合的所有区域。FLAC3D17FLAC3D的的前后处理前后处理n术语术语水平边界 压力结构线模型边界格网点格网区域水 平边 界压力滚 动 底端边界内部开挖边界ID号码号码ID NUMBER-FLAC3D模型中的单元以ID号加以区分，下面的单元有ID号；内部面、节点、区域、体积、历史、表格、显示项和结构单元的全部内容。这帮助用户确定模型中的单元，可用porint命令获得ID号，用户可给内部面、结构单元、历史等赋ID号。实体结构单元同样也有CID号，系统给每一个网格，单元都创建了一个CID号，这与梁，柱等不一样。结构单元结构单元S

15、TRUCTURAL ELEMENT-在FLAC3D中有两种结构单元。二结点，线性单元表示梁，柱作用。三结点，三角平面单元表示面状，结构单元用来模拟土体或岩体中结构支护的相互作用。非线性材料作用可用单元表示。每一结构单元实体（梁，柱，面体）包括三个内容：结点、单个单元（也叫sels）和网格连线，这些内容的不同可区别出梁、桩、面体的作用。步步STEP-因为在FLAC3D是具体代码，问题的计算须分步进行，随步长的增加，现象的有关信息在研究区域传递。对于静态分析，需要给一个具体的步，让其达到平衡状态，典型的问题计算在2000-4000步之间，其他叫法有时间步，循环次。FLAC3D18FLAC3D的的前

16、后处理前后处理n术语术语水平边界 压力结构线模型边界格网点格网区域水 平边 界压力滚 动 底端边界内部开挖边界静态解答静态解答STATIC SOLUTION-如果模型中动量变化率小于了某一可忽略的值，就认为静态或类静态出现了，这通过限定运动方程实现，静态就是模型达到应力平衡，或流体材料受外力后从不稳定到稳定。这种分析方法在FLAC3D中是默认的分析方法，机械的静态分析也可与地下水渗透或热传递问题结合（通过特定设置后，动态问题可由带约束的静态分析代替）。非平衡力非平衡力UNBALANCED FORCE-非平衡力标征静态分析达到机械平衡（或塑性变形前），严格的说平衡时每个节点上的应力矢量都为0。最

17、大应力会自动被监测，当击活step或solve命令时，其值会显示在屏幕上。最大网格力也叫非平衡力或抗平衡力，非平衡力在数值上永远也不能达到0，当最大非平衡力相对加载的力很小时，我们就认为模型达到了平衡状态，如果非平衡力一直保持某一非0值，这就说明模型中可能发生了破坏或塑性变形。FLAC3D19FLAC3D的的前后处理前后处理n术语术语水平边界 压力结构线模型边界格网点格网区域水 平边 界压力滚 动 底端边界内部开挖边界动态解答动态解答DYNAMIC SOLUTION-在动态分析中，求解所有动态方程，动量的产生和消耗都将产生直接影响，在高频率发生或持续时间很短的过程中用到，如地震或爆炸。动态计算

18、是FLAC3D的一个可选模块。（见附录K）大应变大应变/小应变小应变LARGE STRAIN/SMALL STRAIN-FLAC3D默认的都是小应变，也就是说，即使计算出来的位移相对通常区域尺寸很大，节点也不发生相应位移。大应变中，节点根据每一步计算出的位移量发生位移。几何非线性只有用大应变才能实现。FLAC3D20FLAC3D的的前后处理前后处理n命令驱动命令驱动(推荐推荐)程序控制图形界面接口计算模型输出指定本构模型及参数指定初始条件及边界条件，指定结构单元指定接触面指定自定义变量及函数(FISH)求解过程的变量跟踪进行求解模型输出FLAC3D21菜单驱动菜单驱动(计算模式计算模式)命令栏

19、FLAC3D22菜单驱动菜单驱动(Plot)FLAC3D23一个最简单的例子一个最简单的例子gen zon bri size 3 3 3 ;建立网格建立网格model elas ;材料参数材料参数prop bulk 3e8 shear 1e8ini dens 2000 ;初始条件初始条件fix z ran z-.1.1 ;边界条件边界条件fix x ran x-.1.1fix x ran x 2.9 3.1fix y ran y-.1.1fix y ran y 2.9 3.1set grav 0 0-10solve ;求解求解app nstr-10e4 ran z 3 x 1 2 y 1 2s

20、olveRUN FLAC3DFLAC3D24前后处理功能的优点前后处理功能的优点n多种多种zone类型类型n后处理快捷、方便、丰富后处理快捷、方便、丰富n计算过程中的计算过程中的hist变量动态显示变量动态显示nFISH可进行参数化模型设计可进行参数化模型设计n单元状态的可编程单元状态的可编程n计算暂停时的后处理与可保存计算暂停时的后处理与可保存FLAC3D25前后处理功能的缺点前后处理功能的缺点n复杂模型的建模功能不强复杂模型的建模功能不强可以编程导入其他软件形成的网格(比如：Ansys、Adina、GeoCAD)n无等值线的后处理功能无等值线的后处理功能(3D)可编程将.sav文件写入Te

21、cPlot等其他后处理软件n全命令操作，学习困难全命令操作，学习困难n鼠标功能单一鼠标功能单一(双击取击点坐标双击取击点坐标)FLAC3D26用tecplot绘制曲线 1.第一主应力 2.xdisp、ydisp、zdisp、disp用excel做曲线 隧道 1做地表沉降槽（zdisp）2地表横向位移（xdisp）3隧道中线竖向沉降曲线（zdisp）4提取位移矢量图，5显示初期支护结构内力 6显示state（找塑性区）基坑 1做地表沉降槽（zdisp）2提取位移矢量图，3显示初期支护结构内力 4显示state（找塑性区）边坡 做安全系数和应变图后处理后处理FLAC3D27需要掌握需要掌握gen,

22、gen,iniini,app,app,ploplo,solve,solve等建模、初始条件、边界条件、后处理和求解的等建模、初始条件、边界条件、后处理和求解的命令。命令。常用命令常用命令FLAC 3D 程序的编写步骤：程序的编写步骤：1 Config _2 Grid _3 Model _4 求起始的应力平衡_(1)建立x,y坐标与网格的关系，建议使用Gen指示：Gen x1,y1 x2,y2,x3,y3 x4,y4 i=i0,i1 j=j0,j1(2)设定材料性质：prop(3)设定外力：Set Grav,Apply Pressure,ini sxx,Syy(4)设定边界条件：fix,free

23、(5)求起始的应力平衡：solve(6)储存：Save5 求工程的影响_求出区域内的应力分布情况后，再依工程的流程及步骤阶段执行各工程进行过程的影响，建议使用以下的步骤：(1)调出起初的应力平衡：re_.sav(2)设定新的材料性质：model,prop(3)设定新的支撑性质：struct(4)设定新的外力(5)设定边界条件(6)求工程时的应力平衡(7)储存FLAC3D28主要内容主要内容n FLAC3D软件简介软件简介n1、基坑开挖、基坑开挖n2、浅基础的稳定性分析、浅基础的稳定性分析n3、网格的合并联结、网格的合并联结n4、界面的生成、界面的生成n5、隧道的生成、隧道的生成n6、模型材料问

24、题、模型材料问题n7、Interface 合并（联结）问题合并（联结）问题n8、初始条件问题、初始条件问题n9、破坏问题分析破坏问题分析n10、综合实例、综合实例-煤巷应力分析煤巷应力分析FLAC3D29Case-1：土体中挖了一个长宽为土体中挖了一个长宽为2m4m，深为，深为4m的沟，挖土体的沟，挖土体的同时监测周围土体的变形情况的同时监测周围土体的变形情况（模型（模型6m8m8m）初始化不同限定的格网，可以使用GENERATE（生成）命令：gen zone brick size 6 8 8这个命令会建立以一个初始化的格网，这个格网在X方向上有6个分区，Y方向上有8个分区，Z方向有8个分区。

25、所建模型的Z轴在垂直方向上。第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立显示模型格网 命令plot Create TrenchAdd surface yellowAdd axes blackShow创建一个名为“Trench”的块视图并将块表面设为黄色，把视图块的轴设为黑色。show关键字是显示当前视图自动地创建一个视图并设置为当前视图。模型的正面透视图会出现在窗口视图中。按住键盘的X、Y、Z键，图形能X、Y、Z方向旋转。使用M键可以放大视图，使用箭头键可以移动图形块（shift键然后按住相应的键可以反转和缩小图形 显示的图形块中的标题中给出了“Ce

26、nter:”，“Rotation:”，“Dist:”和“Mag:”的数据 FLAC3D30Case-1：土体中挖了一个长宽为土体中挖了一个长宽为2m4m，深为，深为4m的沟，挖土体的沟，挖土体的同时监测周围土体的变形情况（模型的同时监测周围土体的变形情况（模型6m8m8m）创建另一个视图黑体部分显示了新增加的命令。（我们假定Plot命令提示符仍然是被激活的，如果没有，在键入命令之前键入PLOT然后按回车键。）Create Trench2Add surface yellowAdd axes blackSet rotation 17.495 9.847 28.481Set center 3 4 4

27、Set dist 26.12Set magnification 0.8Show第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立FLAC3D31Case-1：土体中挖了一个长宽为土体中挖了一个长宽为2m4m，深为，深为4m的沟，挖土体的沟，挖土体的同时监测周围土体的变形情况（模型的同时监测周围土体的变形情况（模型6m8m8m）使用莫尔-库仑(Mohr-Coulomb)准则弹塑性模型 （在MODEL命令中不指定区域的范围，FLAC3D假定所有的区域都是Mohr-Coulomb性质）Prop bulk=1e8 shear=0.3e8 fric=35 coh=1

28、e10 tens=1e10Prop coh=1e10 tens=1e10 注意：关键参数的值可能被空格键或“=”分开 体积、剪应力、内摩檫角、内聚力和抗压强度是一定的。我们所见的内聚力和抗压强度都给得很大，这仅仅是在重力作用阶段给材料的初始值。实际中，在初始应力阶段材料是弹性的。set grav 0,0,-9.81ini dens=1000/初始化（ini=initial）密度为1000Kg/m3在z 的负方向加以9.81m/sec2的重力加速度（在坐标轴的正方向把重力视为正），为了给模型一个重力，材料密度需预先给定。命令INI用于把模型的每个区域的平均密度设为1000Kg/m3 第一步：初始

29、模型的建立第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立-模型属性模型属性模型属性模型属性FLAC3D32Case-1：土体中挖了一个长宽为土体中挖了一个长宽为2m4m，深为，深为4m的沟，挖土体的沟，挖土体的同时监测周围土体的变形情况（模型的同时监测周围土体的变形情况（模型6m8m8m）边界条件也可在FLAC3D的提示符下输入：fix x range x 0.1 0.1fix x range x 5.9 6.1fix y range y 0.1 0.1fix y range y 7.9 8.1fix z range z 0.1 0.11、前两个命令规定，沿着在x=0和x=

30、6两个平面的栅格边界上的节点在x方向被“固定”，这两个边界面通过Fix命令在”range”限定的范围内下降；2、在y=0和y=8上的节点在y方向上被固定。第3、4条命令规定了他们的下降范围；3、沿着底面边界上的节点在z方向上被固定（z=0），这个平面被第5条命令固定了范围。第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立-边界条件边界条件边界条件边界条件FLAC3D33Case-1：土体中挖了一个长宽为土体中挖了一个长宽为2m4m，深为，深为4m的沟，挖土体的沟，挖土体的同时监测周围土体的变形情况（模型的同时监测周围土体的变形情况（模型6m8m8m）在计

31、算过程中我们想监测所选参数的变化，History 命令可以帮助我们判断是处于平衡状态还是发生了不稳定的破坏。hist n=5 /hist=history设置步长为5hist unbal /最大非平衡力hist gp zdisp 4 4 8/测试栅格点4 4 8的位移我们设定监测参数的变化步长为5（默认步长是10），则每隔5步被选定的参数值就会自动保存在历史列表中，被保存的2个参数是：最大非平衡力和节点z方向上的位移（x=4,y=4,z=8点）；知道一个模型中的最大非平衡力很有用，如果最大非平衡力的值很小或位移很稳定，这说明模型已达到平衡状态。为了考查模型的重力作用，先规定达到平衡的步数。Sol

32、ve 命令可自动判断平衡，*（计算时间，不是实际时间）当键入：set mech force=50solve计算开始，步数和最大非平衡力将在屏幕上显示；当最大非平衡力小于极限值（set命令给定的大小为50N的力）时，程序将停止运行。第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立-监测网格点监测网格点监测网格点监测网格点FLAC3D34Case-1：土体中挖了一个长宽为土体中挖了一个长宽为2m4m，深为，深为4m的沟，挖土体的沟，挖土体的同时监测周围土体的变形情况（模型的同时监测周围土体的变形情况（模型6m8m8m）可以根据最大主应力和z方向位移的历史记录

33、判断它是否达到平衡。当键入：plot hist 1 可得到最大非平衡力的历史记录 最大非平衡力的历史记录图最大非平衡力的历史记录图 z方向上位移历史记录图方向上位移历史记录图(在在x=4,y=4,z=8这点上这点上)hist 2 当最大非平衡力的历史值达到0，而且位移的历史记录稳定。这两个方面都可以反应出模型已达到了平衡状态。得到z方向的位移历史记录;注意：通过hist命令，每一个历史都被连续地从一开始进行数值化，回到FLAC提示符下，键入；print hist得到一个历史列表和相应的数值。第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立-监测网格点监测

34、网格点监测网格点监测网格点FLAC3D35Case-1：土体中挖了一个长宽为土体中挖了一个长宽为2m4m，深为，深为4m的沟，挖土体的沟，挖土体的同时监测周围土体的变形情况（模型的同时监测周围土体的变形情况（模型6m8m8m）同样可以在模型中考察重力效应。现在我们建立一个包含多个图形的复杂窗口，只需键入：plot create GravVplot set plane dip=90 dd=0 origin=3,4,0plot set rot 15 0 20 ;this would be achieved interactivelyplot set center 2.5 4.2 4.0 ;and

35、is show here simply to illustrateplot add bound behind plot add bcont szz plane plot add axes plot show这样我们就建立了一个叫“Grav V”的图形，将它设为当前图层，对此图形设置一个倾角为90的面，（在x-y面上，约定z的负方向为下），倾向为0。（在x-y面上，从y的正方向开始沿顺时针方向计量），过x-y面上一点（x=3,y=4,z=0）,在平面后加一个网格，垂直方向的应力Qzz，最后给模型加轴线，整体图形轮廓与添加图形相反，显示每个区域中心的应力计算值，每个区域的颜色直接对应基本区域的应力

36、，当键入show,我们建立的图形被显示 在初始平衡时的重力应力等值线图 第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立-监测网格点监测网格点监测网格点监测网格点FLAC3D36Case-1：土体中挖了一个长宽为土体中挖了一个长宽为2m4m，深为，深为4m的沟，挖土体的沟，挖土体的同时监测周围土体的变形情况（模型的同时监测周围土体的变形情况（模型6m8m8m）保存初始状态，你可以在进行参数分析的时候提取；保存此状态后，回到FLAC3D的提示符下键入；save trench.Sav于是在默认路径下将建立一个叫“trench.Sav”的文件。现在在模型中建立

37、两个面，键入命令：plot print就可以显示这些视图的列表，动态图将会用箭头标出，图形“base”为草稿视图。键入命令：plot print view将显示当前图形的详细信息。在以上命令中加入图形名称或数字将产生一个图形的详细描述，包含所有设置和项目在图上。当我们把当前的“Grav V”图形改变为“trench”图形时，只要通过输入命令：plot current Trench/Trench视图必须存在（通过create命令生成）这样就让“Trench”接受项目或设置的改变，任何现存的视图都可直接显示出来，而不用先把它设为当前。plot show Trench将在主菜单下显示所有可视的图形

38、第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立第一步：初始模型的建立-监测网格点监测网格点监测网格点监测网格点FLAC3D37Case-1：土体中挖了一个长宽为土体中挖了一个长宽为2m4m，深为，深为4m的沟，挖土体的沟，挖土体的同时监测周围土体的变形情况（模型的同时监测周围土体的变形情况（模型6m8m8m）第二步第二步第二步第二步 开挖巷道开挖巷道开挖巷道开挖巷道 在”plot”下，现在我们在土体内开挖巷道，先键入：prop fric=35 coh=1e3 tens=1e3摩擦角35 粘聚力1103 抗拉强度1103这就给整个区域设置了一个大小为1000Pa的内聚力和抗拉

39、强度，这个强度足以防止在初始阶段发生破坏（即开挖前），但是你要选择一个步数，使它计算完后有破坏的可能性。进行开挖，键入：model null range x=2 4 y=2 6 z=5 10开挖的巷道尺寸（即空材料），通过x-,y-,z-定义一个合适的范围。一个低内聚力和无垂直支护的巷道将发生跨塌，因为我们想仔细考察这一过程，应使用大应变模式。通过键入以下命令实现：set large为了得到更好的显示效果，我们只想看巷道开挖后的位移变化，而不想从重力加载就开始观察，所以我们去掉刚开始x,y,z方向上由于重力发生的位移；*ini xdis=0 ydis=0 zdis=0我们有意设置足够小的内聚力

40、让其破坏，我们现在不用“solve”命令计算，因为用他时要给一个失去平衡的应力极限（为了平衡而选择的），而我们的模拟永远都不会达到平衡；此外，我们可以给定一个时间步，当此时跨塌时显示其结果，这只是一个间接方法，模拟不必在每个计算过程中都达到平衡，因为我们用不着解决一系列数字式 FLAC3D38Case-1：土体中挖了一个长宽为土体中挖了一个长宽为2m4m，深为，深为4m的沟，挖土体的沟，挖土体的同时监测周围土体的变形情况（模型的同时监测周围土体的变形情况（模型6m8m8m）5 5、开挖巷道、开挖巷道、开挖巷道、开挖巷道 用“step”命令step 2000此时FLAC3D 将计算2000步。这

41、样考察计算主要通过显示一些结果。比如，为了产生出一个面的彩色等值线而重复“plot”。但要显示位移等值线图，键入plotcreate Dispcontcopy GravV Dispcont settingadd cont disp plane behindadd axesshow 位移等值线图可显示出由于开挖引起的下沉量位移等值线图可显示出由于开挖引起的下沉量 FLAC3D39Case-1：土体中挖了一个长宽为土体中挖了一个长宽为2m4m，深为，深为4m的沟，挖土体的沟，挖土体的同时监测周围土体的变形情况（模型的同时监测周围土体的变形情况（模型6m8m8m）FLAC3D40主要内容主要内容n

42、FLAC3D软件简介软件简介n1、基坑开挖、基坑开挖n2、浅基础的稳定性分析、浅基础的稳定性分析n3、网格的合并联结、网格的合并联结n4、界面的生成、界面的生成n5、隧道的生成、隧道的生成n6、模型材料问题、模型材料问题n7、Interface 合并（联结）问题合并（联结）问题n8、初始条件问题、初始条件问题n9、破坏问题分析破坏问题分析n10、综合实例、综合实例-煤巷应力分析煤巷应力分析FLAC3D41Case-2：浅基础的稳定性分析浅基础的稳定性分析条形基础是指基础长度远远大于宽度的一种基础形式。按上部结构分为条形基础是指基础长度远远大于宽度的一种基础形式。按上部结构分为墙下条形基础墙下条

43、形基础和柱下条形基础。和柱下条形基础。基础的长度大于或等于基础的长度大于或等于10倍基础的宽度。倍基础的宽度。条形基础的特点是，布置在一条轴线上且与两条以上轴线条形基础的特点是，布置在一条轴线上且与两条以上轴线相交，有时也和相交，有时也和独立基础独立基础相连相连 在此，以条形基础作为研究对象进行说明在此，以条形基础作为研究对象进行说明FLAC3D42Case-2：浅基础的稳定性分析浅基础的稳定性分析FLAC3D43Case-2：浅基础的稳定性分析浅基础的稳定性分析可以用平面应力分析方法解决条形基础问题可以用平面应力分析方法解决条形基础问题土体参数土体参数二维的条形基础的塑性流动二维的条形基础的

44、塑性流动(Tresca 材料材料)（模型（模型20m1m10m）FLAC3D44Case-2：二维的条形基础的塑性流动二维的条形基础的塑性流动(Tresca 材料材料)（模型（模型20m1m10m）FLAC3D45Case-2：二维的条形基础的塑性流动二维的条形基础的塑性流动(Tresca 材料材料)（模型（模型20m1m10m）gen zone brick size 20 1 10;建立矩形的网格区域，其大小为20110plot Create TrenchAdd surface yellowAdd axes blackShow;显示模型格网pausemodel mohr;采用摩尔-库仑模型p

45、rop bul 2.e8 shea 1.e8 cohesion 1.e5;模型的材料性质：体积模量2.0108;剪切模量 1.0108 粘聚力1.0105prop friction 0.dilation 0.tension 1.e10;内摩擦角0 剪胀角0 抗拉强度1.01010fix x range x-.1.1;在x 方向上固定边界 x=0（即只在x 方向上受约束）fix x y z range z-.1.1;在x y z 方向上固定边界 z=0（即为固定约束）fix x y z range x 19.9 20.1;在x y z 方向上固定边界 x=20fix y;约束所有y 方向上运动f

46、ix x y z range x-.1 3.1 z 9.9 10.1;在x y z 方向上固定平面 x=0，3 z=10ini zvel-0.5e-5 range x-.1 3.1 z 9.9 10.1;在平面 x=0，3 z=10上初始化z方向的速度为0.5105FLAC3D46Case-2：二维的条形基础的塑性流动二维的条形基础的塑性流动(Tresca 材料材料)（模型（模型20m1m10m）def p_cons;用fish 语言定义函数p_conspdis1=gp_near(0.,0.,10.);靠近坐标（0，0，10）的节点的地址赋予pdis1pdis2=gp_near(0.,1.,1

47、0.);靠近坐标（0，0，10）的节点的地址赋予pdis2p_sol=(2.+pi);2+赋予 p_solend;用于fish 语言的结尾p_cons;运行函数p_cons;-;p_load:average footing pressure/c;场地的平均压力;c_disp:magnitude of vertical displacement at footing center/a;场地中心的垂直沉降量;-def actLoad;定义函数 actLoad（活荷载）（fish 语言的起始句）FLAC3D47Case-2：二维的条形基础的塑性流二维的条形基础的塑性流动动(Tresca 材料材料)（

48、模型（模型20m1m10m）pnt=gp_head;把网格点指针的首个地址（gp_head）赋予变量 pntpload=0.0;定义变量 pload 为0n=0;定义变量 n 为0loop while pnt#null;loop是循环命令，其条件是pnt 0if gp_zpos(pnt)9.9 then;当变量 pnt 所对应的 z 坐标大于9.9if gp_xpos(pnt)3.1 then;且x 坐标小于3.1 时pload=pload+gp_zfunbal(pnt);变量 pload 为上次循环中的值与变量 pnt 所对应;的z 方向上的不平衡力的和n=n+1;n 自我加1endif;结

49、束里面的ifendif;结束外面的ifpnt=gp_next(pnt);把网格点指针的下一个地址（gp_next）赋予变量 pntendloop;结束循环actLoad=pload/(3.5*z_prop(zone_head,cohesion);这里函数actLoad 的值，为后面的式子所赋予，;其中pload 为循环结束后的最终;值z_prop(zone_head,cohesion)为首个单元;的粘聚力的数值c_disp=-(gp_zdisp(pdis1)+gp_zdisp(pdis2)/7.0;同理，这里函数c_disp 的值，;也为后面的式子所赋予，其中gp_zdisp(pdis1)为;

50、变量pdis1 所对应的网格点在z 方向的位移，;gp_zdisp(pdis2)变量pdis2 所对应;的网格点在z 方向的位移end;fish 语言的结束句FLAC3D48Case-2：二维的条形基础的塑性流二维的条形基础的塑性流动动(Tresca 材料材料)（模型（模型20m1m10m）def p_err;定义函数 p_errp_err=100.*(actLoad-p_sol)/p_sol;函数 p_err 的值为后面式子所定义end;fish 语言的结束句hist n 50;对n 做历史记录hist actLoad;对actLoad 做历史记录hist p_sol;对p_sol 做历史记