FLAC3D数值模拟基础

发布于:2021-06-20 14:27:53

FLAC3D

3D数值模拟基础 FLAC

刘升贵

中国矿业大学力学系
liushg2002@163.com

FLAC3D

主要内容
?
? ? ?

?
? ? ? ? ? ?

FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并(联结)问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析

2

FLAC3D

主要内容
?
? ? ?

?
? ? ? ? ? ?

FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并(联结)问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析

3

FLAC3D

FLAC3D简介
?
? ? ? ? ?

Fast Lagrangian Analysis of Continua 美国Itasca(依泰斯卡)咨询公司开发2D程序(1986) 1990年代初引入中国 有限差分法(FDM) DOS版→2.0 →2.1 →3.0 Itasca其他软件

4

FLAC3D

FLAC3D简介
?

应用:
? 岩土力学/岩石力学分析,例矿体滑坡、煤矿开采沉陷预

测、水利枢纽岩体稳定性分析、采矿巷道稳定性研究等 ? 岩土工程、采矿工程、水利工程、地质工程
?

特色:
? 大应变模拟

? 完全动态运动方程使得FLAC3D在模拟物理上的不稳定过

程不存在数值上的* ? 显示求解具有较快的非线性求解速度

5

FLAC3D

FLAC3D简介
1 承受荷载能力与变形分析:用于边坡稳定和基础设计 2 渐进破坏与坍塌反演:用于硬岩采矿和隧道设计 3 断层构造的影响研究:用于采矿设计 4 施加于地质体锚索支护所提供的支护力研究:岩锚和土钉的设计 5 排水和不排水加载条件下全饱和流体流动和孔隙压力扩散研究:挡土 墙结构的地下水流动和土体固结研究 6 粘性材料的蠕变特性:用于碳酸钾盐矿设计 7 陡滑面地质结构的动态加载:用于地震工程和矿山岩爆研究 8 爆炸荷载和振动的动态响应:用于隧道开挖和采矿活动 9 结构的地震感应:用于土坝设计 10 由于温度诱发荷载所导致的变形和结构的不稳定:高辐射废料地下埋 藏的性能评价 12 大变形材料分析:用于研究粮仓谷物流动及井巷和矿洞中材料的总体 流动
6

FLAC3D

基本原理
建立FLAC计算模型,必须进行以下三个方 面的工作: 1. 有限差分网格 2. 本构特性与材料性质 3. 边界条件与初始条件 完成上述工作后,可以获得模型的初始* 衡状态,也就是模拟开挖前的原岩应力状态 。然后,进行工程开挖或改变边界条件来进 行工程的响应分析。





生成网格并调整网格的形状; 持续的运动和连续的物质属性; 特定的边界条件和初始条件。 到达*衡状态

效果是否符合要求

作如下改变: 开挖模型的物质属性 改变边界条件

实施求解

网格确定问题的几何尺寸;持续的运动和 连续的物质属性决定了模型的扰动(如由 于开挖引起的变形)形式;边界条件和初 始条件确定了模型的初始状态(没有引起 扰动或变形的状态)。
参数是否调整 效果是否符合要求





7

FLAC3D

基本原理
?
? ? ? ? ?

有限差分法 Lagrangian网格 空间混合离散技术 Lagrangian格式动量*衡方程 FLAC3D的求解过程 FLAC3D的本构模型

8

FLAC3D

有限差分法
?
? ? ? ?

古老的方法(上世纪40年代) 用差分格式转化控制方程中的微商格式 流体力学;土工渗流问题;固结 FDM & FEM的混合求解 FDM的新进展

9

FLAC3D

Lagrangian网格
?

源自流体力学中的拉格朗日法
? 跟踪流体质点的运动状态 ? 跟踪固体力学中结点,按时步用

Lagrangian法研究网格节点的运动
?

?

节点和单元随材料移动,边界和接 触面与单元的边缘一致 固体力学大变形理论

法国数学家、物理学家拉格朗日
10

FLAC3D

FLAC3D的求解过程
对所有的网格节点
速度 *衡方程 (动量方程) 节点力

Gauss定律

单元积分

对所有单元
应力—应变关系 (本构模型)

应变率

新的应力

11

FLAC3D

FLAC3D中的本构模型
?
?

开挖模型null 3个弹性模型
? 各向同性弹性

? 横观各向同性弹性
? 正交各向同性弹性

?

8个塑性模型(Drucker-Prager模型、MorhCoulomb模型、应变硬化/软化模型、遍布节理 模型、双线性应变硬化/软化遍布节理模型、修正 剑桥模型和胡克布朗模型)
12

FLAC3D

FLAC3D中的本构模型
模型 空模型 线弹性模型 正交各向同性 弹性 横观各向同性 弹性 德鲁克-普拉格 模型 摩尔-库仑 模型 应变硬化/软化摩尔-库仑 模型 遍布解理模型 双线性应变硬化/软化遍 布解理模型 双屈服面塑性模型 修正剑桥模型 胡克-布朗模型 材料特性 空 均匀各向同性的线形本构关系 正交各向同性材料 横观各向同性弹性(即板岩) 极限分析,底摩擦角的软粘土 松散或胶结的粒状材料:土,岩石,混 凝土 存在非线性硬化或软化的粒状材料 具有强度各向异性的层状材料(即板岩) 具有非线性材料硬化或软化的层状材料 轻胶结的粒状材料,在压力作用下导致 永久体积减小 变形和抗剪强度是体变的函数 各向同性的岩石材料 粘土 岩石 实际应用 孔洞,开挖,后续施工材料(如回填) 低于强度极限的人工材料(如钢 铁);安全系数计算 不超过强度极限的柱状玄武岩 不超过强度极限的层压材料 与隐式有限元程序相比的常用模型 岩土力学通用模型(边坡稳定性分 析,地下开挖) 破坏后研究(失稳过程,立柱屈服, 顶板崩落) 松散沉积地层中的开挖 层状材料破坏后研究

13

FLAC3D

FLAC3D的前后处理
?

术语
格网点 结构线

区域(Zone)……有限差分划分的带在几何上是最小的区 域,在在这个区域里的每一个现象的变化,如应力应变都可以 估计出,。各种形状的多面体(立方体、楔形、锥体、四面体 格网区域 等)可用来构造模型并可用plot显示出来。每一个多面体可能 有一套或两套表层设置,这由5个四面体组成。默认的情况下, 水 * 边 界 压力 两个表层设置用在对计算精度要求高的情况下,区域的另外一 种叫法是要素。 栅格点(GridPoint)……栅格点是有限差分单元的角点。 一个多面体可能有5个、6个、7个或8个网格点,主要取决于多 面体的形状。给定每个节点的x,y和z值这样就具体确定了有限 差分单元,。其他叫法有:节点,交点。 有限差分栅格(Finite Difference Grid)……有限差分 网格是研究区域中一个或多个通过物理边界连接的有限差分单 元的集合。另一个叫法是网格,有限差分网格也可以标识出模 型中每个状态的存储位置,FLAC3D所生成的矢量都保存在节 点上(如:受力、速度、位移)。标量和张量保存在单元的中 心(如应力、材料属性)。

模型边界

水*边 界 压力

内部开挖边界 滚动底 端边界

14

FLAC3D

FLAC3D的前后处理
?

术语
格网点 结构线

模型边界(Model Boundary)……即有限差分网格的 外围,内部边界也同样是模型边界(如网格中的空洞)。 边界条件(Boundary Condition)……即模型边界的格网区域 约束条件或控制条件的给定(如:限制位移、渗透条件、绝 水 * 热条件)。 边 界 初始条件(Initial Conditions)……即在对模型加载 压力 或开挖等作用前的各种参数状态。 基本模型(Constitutive Model)……基本模型(材料 模型)即规定了FLAC3D模型中某一区域的变形或强度效应, 可用大量基本模型去*视地质材料,可以单独定义FLAC3D 模型中的基本模型和材料模型。 空单元(Null Zone)……空单元表示此区域为空(就 象没有材料一样)。 次级网格(SUB-GRID)……有限差分网格可由次级网 格组成,它可用来在模型中创建不同形状的区域,次级网格 是分别生成,可进行合并和连接。

模型边界

水*边 界 压力

内部开挖边界 滚动底 端边界

15

FLAC3D

FLAC3D的前后处理
?

术语
格网点 结构线 格网区域 模型边界 水 * 边 界 压力

附属接触面ATTACHED FACES……附属栅格面是由被 划分的次级栅格组成的网格接触或合并的面栅格面,接触面 必须是共面或接触,每个面的节点不一定一样,不同总密度 的次级网格可以接触。
接触面INTERFACE------即次级网格在计算过程中可以 分开(滑动,开裂)的两部分之间的面,可表示不连续的物 理特征,如,断层,节理面或材料性质突变的临界面。 范围INTERFACE------范围是对FLAC3D模型空间值的 一个描述,可给定一个命令的作用范围,即使模型发生运动--不影响模型中的区域和节点位置,一个范围或范围确定的 单元也不发生改变。范围或范围内的单元所包含的区域,也 可以用一个单元的ID号来表示,它与区域,节点,或结构单 元密切联系。 集合GROUP-----在FLAC3D模型中他们有共同的名称, 由于限定具体命令的对象,如model命令对某一集合设置为 某种材料,任何命令加于集合名称也就相当于作用于这一集 合的所有区域。

水*边 界 压力

内部开挖边界 滚动底 端边界

16

FLAC3D

FLAC3D的前后处理
?

术语
格网点 结构线

ID号码ID NUMBER-----FLAC3D模型中的单元以ID号 加以区分,下面的单元有ID号;内部面、节点、区域、体积、 历史、表格、显示项和结构单元的全部内容。这帮助用户确 定模型中的单元,可用porint命令获得ID号,用户可给内部格网区域 面、结构单元、历史等赋ID号。 实体结构单元同样也有CID号,系统给每一个网格, 水 * 边 界 单元都创建了一个CID号,这与梁,柱等不一样。 压力 结构单元STRUCTURAL ELEMENT------在FLAC3D 中有两种结构单元。二结点,线性单元表示梁,柱作用。三 结点,三角*面单元表示面状,结构单元用来模拟土体或岩 体中结构支护的相互作用。非线性材料作用可用单元表示。 每一结构单元实体(梁,柱,面体)包括三个内容: 结点、单个单元(也叫sels)和网格连线,这些内容的不同 可区别出梁、桩、面体的作用。 步STEP------因为在FLAC3D是具体代码,问题的计算 须分步进行,随步长的增加,现象的有关信息在研究区域传 递。对于静态分析,需要给一个具体的步,让其达到*衡状 态,典型的问题计算在2000-4000步之间,其他叫法有时 间步,循环次。

模型边界

水*边 界 压力

内部开挖边界 滚动底 端边界

17

FLAC3D

FLAC3D的前后处理
?

术语
格网点 结构线

静态解答STATIC SOLUTION-----如果模型中动量变化 率小于了某一可忽略的值,就认为静态或类静态出现了,这格网区域 通过限定运动方程实现,静态就是模型达到应力*衡,或流 体材料受外力后从不稳定到稳定。这种分析方法在FLAC3D 水 * 边 界 中是默认的分析方法,机械的静态分析也可与地下水渗透或 压力 热传递问题结合(通过特定设置后,动态问题可由带约束的 静态分析代替)。 非*衡力UNBALANCED FORCE-------非*衡力标征 静态分析达到机械*衡(或塑性变形前),严格的说*衡时 每个节点上的应力矢量都为0。最大应力会自动被监测,当 击活step或solve命令时,其值会显示在屏幕上。最大网格 力也叫非*衡力或抗*衡力,非*衡力在数值上永远也不能 达到0,当最大非*衡力相对加载的力很小时, 我们就认为 模型达到了*衡状态,如果非*衡力一直保持某一非0值, 这就说明模型中可能发生了破坏或塑性变形。

模型边界

水*边 界 压力

内部开挖边界 滚动底 端边界

18

FLAC3D

FLAC3D的前后处理
?

术语
格网点 结构线 格网区域

动态解答DYNAMIC SOLUTION------在动态分析中, 水 * 求解所有动态方程,动量的产生和消耗都将产生直接影响, 边 界 压力 在高频率发生或持续时间很短的过程中用到,如地震或爆炸。 动态计算是FLAC3D的一个可选模块。(见附录K) 大应变/小应变LARGE STRAIN/SMALL STRAIN -------FLAC3D默认的都是小应变,也就是说,即使计算出来 的位移相对通常区域尺寸很大,节点也不发生相应位移。大 应变中,节点根据每一步计算出的位移量发生位移。几何非 线性只有用大应变才能实现。
滚动底 端边界

模型边界

水*边 界 压力

内部开挖边界

19

FLAC3D

FLAC3D的前后处理
?

命令驱动(推荐)
? 程序控制 ? 图形界面接口 ? 计算模型输出 ? 指定本构模型及参数 ? 指定初始条件及边界条件,指定结构单元 ? 指定接触面 ? 指定自定义变量及函数(FISH)

? 求解过程的变量跟踪
? 进行求解 ? 模型输出
20

FLAC3D

菜单驱动(计算模式)

命令栏

21

FLAC3D

菜单驱动(Plot)

22

FLAC3D

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

RUN FLAC3D
23

FLAC3D

前后处理功能的优点
?
? ? ? ? ?

多种zone类型 后处理快捷、方便、丰富 计算过程中的hist变量动态显示 FISH可进行参数化模型设计 单元状态的可编程 计算暂停时的后处理与可保存

24

FLAC3D

前后处理功能的缺点
?

复杂模型的建模功能不强
? 可以编程导入其他软件形成的网格(比如:Ansys、

Adina、GeoCAD)
?

无等值线的后处理功能(3D)
? 可编程将.sav文件写入TecPlot等其他后处理软件

?

?

全命令操作,学*困难 鼠标功能单一(双击取击点坐标)

25

FLAC3D

后处理
用tecplot绘制曲线 1.第一主应力 2.xdisp、ydisp、zdisp、disp 用excel做曲线 隧道 1做地表沉降槽(zdisp) 2地表横向位移(xdisp) 3隧道中线竖向沉降曲线(zdisp) 4提取位移矢量图, 5显示初期支护结构内力 6显示state(找塑性区) 基坑 1做地表沉降槽(zdisp) 2提取位移矢量图, 3显示初期支护结构内力 4显示state(找塑性区) 边坡 做安全系数和应变图
26

FLAC3D

常用命令
需要掌握gen, ini, app, plo, 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 (5)求起始的应力*衡:solve (6)储存:Save 5 求工程的影响________ 求出区域内的应力分布情况后,再依工程的流程及步骤阶段执行各工程进行过 程的影响,建议使用以下的步骤: (1)调出起初的应力*衡:re_____ .sav (2)设定新的材料性质:model,prop (3)设定新的支撑性质:struct (4)设定新的外力 (5)设定边界条件 (6)求工程时的应力*衡 (7)储存

27

FLAC3D

主要内容
?
? ? ?

?
? ? ? ? ? ?

FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并(联结)问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析

28

FLAC3D

Case-1: 土体中挖了一个长宽为2m×4m,深为4m的沟,挖土体 的同时监测周围土体的变形情况 (模型6m×8m×8m)
第一步:初始模型的建立
初始化不同限定的格网,可以使用 GENERATE(生成)命令: gen zone brick size 6 8 8 这个命令会建立以一个初始化的格 网,这个格网在X方向上有6个分区, Y方向上有8个分区,Z方向有8个分区。 所建模型的Z轴在垂直方向上。

显示模型格网 命令plot
Create Trench Add surface yellow Add axes black Show

显示的图形块中的标题中给出了“Center:” , “Rotation:”,“Dist:” 和 “Mag:”的 数据 创建一个名为“Trench”的块视图并将块表面设为黄 色,把视图块的轴设为黑色。show关键字是显示当 前视图——自动地创建一个视图并设置为当前视图。 模型的正面透视图会出现在窗口视图中。按住键盘的 X、Y、Z键,图形能X、Y、Z方向旋转。使用M键可 以放大视图,使用箭头键可以移动图形块(shift键然 后按住相应的键可以反转和缩小图形

29

FLAC3D

Case-1: 土体中挖了一个长宽为2m×4m,深为4m的沟,挖土体 的同时监测周围土体的变形情况(模型6m×8m×8m)
第一步:初始模型的建立
创建另一个视图——黑体部分显示了新增加的命令。(我们假定Plot>命令提示符仍 然是被激活的,如果没有,在键入命令之前键入PLOT然后按回车键。) Create Trench2 Add surface yellow Add axes black Set rotation 17.495 9.847 28.481 Set center 3 4 4 Set dist 26.12 Set magnification 0.8 Show

30

FLAC3D

Case-1: 土体中挖了一个长宽为2m×4m,深为4m的沟,挖土体 的同时监测周围土体的变形情况(模型6m×8m×8m)
第一步:初始模型的建立-模型属性
使用莫尔-库仑(Mohr-Coulomb)准则弹塑性模型
(在MODEL命令中不指定区域的范围,FLAC3D假定所有的区域都是Mohr-Coulomb性质) Prop bulk=1e8 shear =0.3e8 fric=35 coh=1e10 tens=1e10 Prop coh=1e10 tens=1e10

注意:关键参数的值可能被空格键或“= ”分开 体积、剪应力、内摩檫角、内聚力和抗压强度是一定的。我们所见的内聚力和抗 压强度都给得很大,这仅仅是在重力作用阶段给材料的初始值。实际中,在初始 应力阶段材料是弹性的。
set grav 0, 0, -9.81 ini dens = 1000//初始化(ini=initial)密度为1000Kg/m3

在z 的负方向加以9.81m/sec2的重力加速度(在坐标轴的正方向把重力视为正) ,为了给模型一个重力,材料密度需预先给定。命令INI用于把模型的每个区域的 *均密度设为1000Kg/m3
31

FLAC3D

Case-1: 土体中挖了一个长宽为2m×4m,深为4m的沟,挖土体 的同时监测周围土体的变形情况(模型6m×8m×8m)
第一步:初始模型的建立-边界条件
边界条件也可在FLAC3D的提示符下输入: fix fix fix fix fix x x y y z range range range range range x x y y z –0 .1 5 .9 –0 .1 7 .9 –0 .1 0. 1 6. 1 0. 1 8. 1 0. 1

1、前两个命令规定,沿着在x=0和x=6两个*面的栅格边界上的节点在x方向被“固 定”,这两个边界面通过Fix命令在”range”限定的范围内下降;
2、在y=0和y=8上的节点在y方向上被固定。第3、4条命令规定了他们的下降范围; 3、沿着底面边界上的节点在z方向上被固定(z=0),这个*面被第5条命令固定了范 围。

32

FLAC3D

Case-1: 土体中挖了一个长宽为2m×4m,深为4m的沟,挖土体 的同时监测周围土体的变形情况(模型6m×8m×8m)
第一步:初始模型的建立-监测网格点
在计算过程中我们想监测所选参数的变化,History 命令可以帮助我们判断是 处于*衡状态还是发生了不稳定的破坏。
hist n=5 //hist=history设置步长为5 hist unbal //最大非*衡力 hist gp zdisp 4 4 8 //测试栅格点4 4 8的位移

我们设定监测参数的变化步长为5(默认步长是10),则每隔5步被选定的参数值就会自动 保存在历史列表中,被保存的2个参数是:最大非*衡力和节点z方向上的位移 (x=4,y=4,z=8点);知道一个模型中的最大非*衡力很有用,如果最大非*衡力的值很 小或位移很稳定,这说明模型已达到*衡状态。
为了考查模型的重力作用,先规定达到*衡的步数。 Solve 命令可自动判断*衡,*(计算 时间,不是实际时间)当键入: set mech force=50 solve 计算开始,步数和最大非*衡力将在屏幕上显示;当最大 非*衡力小于极限值(set命令给定的大小为50N的力) 时,程序将停止运行。
33

FLAC3D

Case-1: 土体中挖了一个长宽为2m×4m,深为4m的沟,挖土体 的同时监测周围土体的变形情况(模型6m×8m×8m)
第一步:初始模型的建立-监测网格点
可以根据最大主应力和z方向位移的历史记录判断它是否达到*衡。当键入: plot hist 1 可得到最大非*衡力的历史记录 最大非*衡力的历史记录图 hist 2 得到z方向的位移历史记录; 当最大非*衡力的历史值达到0,而且位移的历史记录稳定。 这两个方面都可以反应出模型已达到了*衡状态。 注意:通过hist命令,每一个历史都被连续地从一开 始进行数值化,回到FLAC提示符下,键入; print hist 得到一个历史列表和相应的数值。
z方向上位移历史记录图(在 x=4,y=4,z=8这点上)

34

FLAC3D

Case-1: 土体中挖了一个长宽为2m×4m,深为4m的沟,挖土体 的同时监测周围土体的变形情况(模型6m×8m×8m)
第一步:初始模型的建立-监测网格点
同样可以在模型中考察重力效应。现在我们建立一个包含多个图形的复杂窗口,只需键入: plot create GravV plot set plane dip=90 dd=0 origin=3,4,0 plot set rot 15 0 20 ;this would be achieved interactively plot set center 2.5 4.2 4.0 ;and is show here simply to illustrate plot 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,最后给模型加 轴线,整体图形轮廓与添加图形相反,显示每个区 域中心的应力计算值,每个区域的颜色直接对应基 本区域的应力,当键入show,我们建立的图形被显 示

在初始*衡时的重力应力等值线图
35

FLAC3D

Case-1: 土体中挖了一个长宽为2m×4m,深为4m的沟,挖土体 的同时监测周围土体的变形情况(模型6m×8m×8m)
第一步:初始模型的建立-监测网格点

保存初始状态,你可以在进行参数分析的时候提取;保存此状态后,回到FLAC3D的提示符下键入; save trench.Sav 于是在默认路径下将建立一个叫“trench.Sav”的文件。 现在在模型中建立两个面,键入命令: plot print 就可以显示这些视图的列表,动态图将会用箭头标出,图形“ base”为草稿视图。 键入命令: plot print view 将显示当前图形的详细信息。在以上命令中加入图形名称或数字将产生一个图形的详细描述, 包含所有设置和项目在图上。当我们把当前的“ Grav V”图形改变为“trench”图形时,只要通 过输入命令: plot current Trench//Trench视图必须存在(通过create命令生成) 这样就让“Trench”接受项目或设置的改变,任何现存的视图都可直接显示出来,而不用先把 它设为当前。 plot show Trench 将在主菜单下显示所有可视的图形
36

FLAC3D

Case-1: 土体中挖了一个长宽为2m×4m,深为4m的沟,挖土体 的同时监测周围土体的变形情况(模型6m×8m×8m)
第二步 开挖巷道
在”plot”下,现在我们在土体内开挖巷道,先键入: prop fric=35 coh=1e3 tens=1e3

摩擦角35° 粘聚力1×103 抗拉强度1×103
这就给整个区域设置了一个大小为1000Pa的内聚力和抗拉强度,这个强度足以防止在初始阶段 发生破坏(即开挖前),但是你要选择一个步数,使它计算完后有破坏的可能性。进行开挖,键入: 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
我们有意设置足够小的内聚力让其破坏,我们现在不用“solve”命令计算,因为用他时要给一个失去*衡的应力极 限(为了*衡而选择的),而我们的模拟永远都不会达到*衡;此外,我们可以给定一个时间步,当此时跨塌时显 示其结果,这只是一个间接方法,模拟不必在每个计算过程中都达到*衡,因为我们用不着解决一系列数字式 37

FLAC3D

Case-1: 土体中挖了一个长宽为2m×4m,深为4m的沟,挖土体 的同时监测周围土体的变形情况(模型6m×8m×8m)
5、开挖巷道
用“step”命令 step 2000 此时FLAC3D 将计算2000步。这样考察计算主要通过显示一些结果。比如,为了产 生出一个面的彩色等值线而重复“plot”。但要显示位移等值线图,键入 FLAC3D 2.10 plot Step 1687 Model Perspective 16:17:54 Sat May 07 2011 create Dispcont Rotation: copy GravV Dispcont setting Center: X: 2.500e+000 X: 15.000 Y: 4.200e+000 Y: 0.000 add cont disp plane behind Z: 4.000e+000 Z: 20.000 Dist: 3.387e+001 Mag.: 1 Ang.: 22.500 add axes Plane Origin: Plane Normal: show X: 3.000e+000 X: 0.000e+000
Y: 4.000e+000 Z: 0.000e+000 Y: 1.000e+000 Z: 6.123e-017

Contour of Displacement Mag.

位移等值线图可显示出由于开挖引起的下沉量

Plane: on behind Magfac = 1.000e+000 0.0000e+000 to 0.0000e+000 0.0000e+000 to 2.5000e-002 2.5000e-002 to 5.0000e-002 5.0000e-002 to 7.5000e-002 7.5000e-002 to 1.0000e-001 1.0000e-001 to 1.2500e-001 1.2500e-001 to 1.5000e-001 1.5000e-001 to 1.7500e-001 1.7500e-001 to 2.0000e-001 2.0000e-001 to 2.2500e-001 2.2500e-001 to 2.5000e-001 2.5000e-001 to 2.5498e-001 Interval = 2.5e-002 Itasca Consulting Group, Inc. Minneapolis, MN USA

Z Y X

38

FLAC3D

Case-1: 土体中挖了一个长宽为2m×4m,深为4m的沟,挖土体 的同时监测周围土体的变形情况(模型6m×8m×8m)

39

FLAC3D

主要内容
?
? ? ?

?
? ? ? ? ? ?

FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并(联结)问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析

40

FLAC3D

Case-2:浅基础的稳定性分析

在此,以条形基础作为研究对象进行说明 条形基础是指基础长度远远大于宽度的一种基础形式。按上部结构分为墙下条形基础和柱下条形基础。 基础的长度大于或等于10倍基础的宽度。 条形基础的特点是,布置在一条轴线上且与两条以上轴线 相交,有时也和独立基础相连

41

FLAC3D

Case-2:浅基础的稳定性分析

42

FLAC3D

Case-2:浅基础的稳定性分析
二维的条形基础的塑性流动(Tresca 材料) (模型20m×1m×10m)

可以用*面应力分析方法解决条形基础问题

土体参数

43

FLAC3D

Case-2: 二维的条形基础的塑性流动(Tresca 材料) (模型20m×1m×10m)
FLAC3D 2.10
Settings: Model Perspective 17:44:23 Sat May 07 2011 Center: X: 1.000e+001 Y: 5.000e-001 Z: 5.000e+000 Dist: 5.580e+001 Rotation: X: 40.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500

Surface
Magfac = 0.000e+000

Axes
Linestyle

Y Z X

Itasca Consulting Group, Inc. Minneapolis, MN USA

44

FLAC3D

Case-2: 二维的条形基础的塑性流动(Tresca 材料) (模型20m×1m×10m)
gen zone brick size 20 1 10 ;建立矩形的网格区域,其大小为20×1×10 plot Create Trench Add surface yellow Add axes black Show ;显示模型格网 pause
FLAC3D 2.10
Settings: Model Perspective 17:50:59 Sat May 07 2011 Center: X: 1.000e+001 Y: 5.000e-001 Z: 5.000e+000 Dist: 5.580e+001 Rotation: X: 0.000 Y: 0.000 Z: 0.000 Mag.: 1 Ang.: 22.500

Surface
Magfac = 0.000e+000

Axes
Linestyle

Z

Y

X

model mohr ;采用摩尔-库仑模型 prop bul 2.e8 shea 1.e8 cohesion 1.e5 ;模型的材料性质:体积模量2.0×108 ;剪切模量 1.0×108 粘聚力1.0×105 prop friction 0. dilation 0. tension 1.e10 ;内摩擦角0° 剪胀角0° 抗拉强度1.0×1010 fix 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 =20 fix y ;约束所有y 方向上运动 fix x y z range x -.1 3.1 z 9.9 10.1 ;在x y z 方向上固定*面 x =0,3 z= 10 ini zvel -0.5e-5 range x -.1 3.1 z 9.9 10.1 ;在*面 x =0,3 z= 10上初始化z方向的速度为?0.5×10?5

Itasca Consulting Group, Inc. Minneapolis, MN USA

45

FLAC3D

Case-2: 二维的条形基础的塑性流动(Tresca 材料) (模型20m×1m×10m)
def p_cons ;用fish 语言定义函数p_cons pdis1 = gp_near(0.,0.,10.) ;靠*坐标(0,0,10)的节点的地址赋予pdis1 pdis2 = gp_near(0.,1.,10.) ;靠*坐标(0,0,10)的节点的地址赋予pdis2 p_sol = (2. + pi) ; 2 +π 赋予 p_sol end ;用于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 语言的起始句)
46

FLAC3D

pnt = gp_head ;把网格点指针的首个地址(gp_head)赋予变量 pnt pload = 0.0 ;定义变量 pload 为0 n=0 ;定义变量 n 为0 loop while pnt # null ;loop是循环命令,其条件是pnt ≠ 0 if gp_zpos(pnt) > 9.9 then ;当变量 pnt 所对应的 z 坐标大于9.9 if gp_xpos(pnt) < 3.1 then ;且x 坐标小于3.1 时 pload = pload + gp_zfunbal(pnt) ;变量 pload 为上次循环中的值与变量 pnt 所对应 ;的z 方向上的不*衡力的和 n = n+ 1 ;n 自我加1 endif ;结束里面的if endif ; 结束外面的if

Case-2: 二维的条形基础的塑性流 动(Tresca 材料) (模型20m×1m×10m)

pnt = gp_next(pnt) ; 把网格点指针的下一个地址(gp_next)赋予变量 pnt endloop ;结束循环 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)为 ;变量pdis1 所对应的网格点在z 方向的位移, ;gp_zdisp(pdis2) 变量pdis2 所对应 ;的网格点在z 方向的位移 end ; fish 语言的结束句
47

FLAC3D

Case-2: 二维的条形基础的塑性 流动(Tresca 材料) (模型20m×1m×10m)

def p_err ;定义函数 p_err p_err = 100. * (actLoad - p_sol) / p_sol ;函数 p_err 的值为后面式子所定义 end ;fish 语言的结束句 hist n 50 ;对n 做历史记录 hist actLoad ;对actLoad 做历史记录 hist p_sol ;对p_sol 做历史记录 hist c_disp ;对c_disp 做历史记录 hist unbal ;对unbal(不*衡力)做历史记录 plot sketch ;绘出轮廓线 plot add vel ;添加矢量 plot add his 1 2 vs 3 ;绘出历史记录n,actLoad , p_sol cyc 15000 ;运行15000 时步 print p_err ;输出函数 p_err save pran.sav ;形成sav 文件,并保存为pran.sav 48

FLAC3D

Case-2: 二维的条形基础的塑性流动(Tresca 材料) (模型20m×1m×10m)
FLAC3D 2.10
Step 15000 Model Perspective 18:09:41 Sat May 07 2011 Center: X: 1.000e+001 Y: 5.000e-001 Z: 5.000e+000 Dist: 5.580e+001 Rotation: X: 0.000 Y: 0.000 Z: 0.000 Mag.: 1 Ang.: 22.500
5.0

4.5

4.0

Surface
Magfac = 0.000e+000

3.5

Axes
Linestyle

3.0

Velocity
Maximum = 5.002e-006 Linestyle

2.5

History
1 actLoad (FISH function) Linestyle 2.318e-001 <-> 5.275e+000 2 p_sol (FISH symbol) Linestyle 5.142e+000 <-> 5.142e+000

Z
2.0

1.5

Y

X

1.0

Vs.
3 c_disp (FISH symbol) 7.143e-005 <-> 2.143e-002 Itasca Consulting Group, Inc. Minneapolis, MN USA
0.5

0.2

0.4

0.6

0.8

1.0 x10^-2

1.2

1.4

1.6

1.8

2.0

49

FLAC3D

Case-2: 二维的条形基础的塑性流动(Tresca 材料) (模型20m×1m×10m)

50

FLAC3D

主要内容
?
? ? ?

FLAC3D 2.10
Settings: Model Perspective 20:06:49 Sat May 07 2011 Center: X: 0.000e+000 Y: 0.000e+000 Z: -8.500e+000 Dist: 7.360e+001 Rotation: X: 10.000 Y: 0.000 Z: 20.000 Mag.: 0.8 Ang.: 22.500
Z

Y

X

?
? ? ? ? ? ?

FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并(联结)问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析
Block Group Axes
lower_block upper_block Linestyle

Itasca Consulting Group, Inc. Minneapolis, MN USA

51

FLAC3D

Case-3:网格的合并联结

FLAC3D 2.10
Settings: Model Perspective 19:59:41 Sat May 07 2011 Center: X: 0.000e+000 Y: 0.000e+000 Z: -8.500e+000 Dist: 7.360e+001 Rotation: X: 10.000 Y: 0.000 Z: 20.000 Mag.: 0.8 Ang.: 22.500
Z

plot create grid4 plot add block group red yellow ;把在group 中的部分染成红色和黄色 Block Group plot add axes black ;添加坐标轴线为黑色 Axes pause gen zone brick size 8 8 8 p0 -10,-10,-20 p1 10,-10,-20 & p2 -10,10,-20 p3 -10,-10,0 & group lower_block ;建立块体网格,大小为8×8×8 ,其尺寸是p0 、p1、 p2、p3来决定,其中p0为 ;起始点的坐标,p1 为x 方向,p2 为y 方向,p3 位z 方向的坐标,并把这些区域 FLAC3D 2.10 ;建立成一个群,名为lower_block。 Settings: Model Perspective 20:11:32 Sat May 07 2011 gen zone brick size 2 ,2,1 p0 0,0,1 p1 5,0,1 & Center: Rotation: p2 0,5,1 p3 0,0,3 & X: 0.000e+000 X: 40.000 group upper_block Y: 0.000e+000 Y: 0.000 Z: -8.500e+000 Z: 20.000 ;建立另外一个块体网格,其大小为2× 2×1,其尺寸及位置是p0 、p1、 0.8 p2、p3 Dist: 7.360e+001 Mag.: Ang.: 22.500 ;来决定,具体如上所述,并把这个区域归为一个群,名为upper_block。 plot show Block Group ;在视图窗中显示 lower_block upper_block pause Axes ;暂停 Linestyle ini z add -1 range group upper_block ;群upper_block 的所有单元,在z 方向上向下移动1m ;gen zone brick size 2 ,2,1 p0 0,0,0 p1 5,0,0 & ; p2 0,5,0 p3 0,0,2 & ; group upper_block ; 重新建立另外一个块体网格,其大小为2×2×1,其尺寸及位置是p0 、p1、 ;p2、p3 来决定,具体如上所述,并把这个区域归为一个群,名为upper_block。 gen merge 1e-5 range z 0 ;此命令是接触面单元合并成一个整体。
Linestyle

Y

X

lower_block upper_block

Itasca Consulting Group, Inc. Minneapolis, MN USA

Z Y X

Itasca Consulting Group, Inc. Minneapolis, MN USA

52

FLAC3D

主要内容
?
? ? ?

?
? ? ? ? ? ?

FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并(联结)问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析

53

FLAC3D

Case-4:界面的生成
; Create Base gen zone brick size 3 3 3 p0 (0,0,0) p1 (3,0,0) p2 (0,3,0) p3 (0,0,1.5) p4 (3,3,0) p5 (0,3,1.5) p6 (3,0,4.5) p7 (3,3,4.5) group Base ;建立另外一个块体网格,其大小为3×3×3,其尺寸及位置是p0 、p1、 p2、p3、 ;p4、p5、p6 和p7 等楔型体的8 个角点坐标来决定,并把这个区域归为一个群, ;名为Base(即下面的底座)。 ; Create Top - 1 unit high for initial spacing gen zone brick size 3 3 3 p0 (0,0,1.5) p1 (3,0,4.5) p2 (0,3,1.5) p3 (0,0,6) p4 (3,3,4.5) p5 (0,3,6) p6 (3,0,6) p7 (3,3,6) group Top range group Base not ;建立另外一个块体网格,其大小为3×3×3,其尺寸及位置是p0 、p1、 p2、p3、 ;p4、p5、p6 和p7 等楔型体的8 个角点坐标来决定,并把这个区域归为一个群, ;名为Top(即上面的部分)。 FLAC3D 3.00 gen separate Top Settings: Model Perspective 10:46:08 Sun May 08 2011 ;使两部分的接触网格分离为两部分 Center: ; Create interface elements on the top surface of the base1.500e+000 Rotation: X: X: 20.000 Y: 1.500e+000 Y: 0.000 Z: 3.000e+000 Z: 20.000 ;在两部分之间添加交界面单元 Dist: 1.809e+001 Mag.: 1 Ang.: 22.500 interface 1 wrap Base Top Surface ;在这两部分之间添加接触单元 Magfac = 0.000e+000 ; Interface Locations plot create view_int ;显示,并创建标题view_int plot add surface ;显示表面 plot add interface red ;显示交互面为红色 plot show ;打开图形 Itasca Consulting Group, Inc. Minneapolis, MN USA save int.sav ; 形成sav 文件,并保存为int.sav

54

FLAC3D

主要内容
FLAC3D 2.10

?
? ? ?

?
? ? ? ? ? ?

FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并(联结)问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析
Surface

Settings: Model Perspective 10:31:50 Sun May 08 2011 Center: X: 0.000e+000 Y: 5.000e+000 Z: 0.000e+000 Dist: 6.030e+001 Rotation: X: 10.000 Y: 0.000 Z: 10.000 Mag.: 1 Ang.: 22.500

Magfac = 0.000e+000

Itasca Consulting Group, Inc. Minneapolis, MN USA

55

FLAC3D

Case-5:隧道的生成
new ;新建 def setup ;定义函数setup numy = 8 ;定义常量numy 为8 depth = 10.0 ;depth 为10 end ;结束对函数的定义 setup ;运行函数setup gen zone radcyl size 5 numy 8 10 & p0 0 0 0 p1 10 0 2 p2 0 depth 0 p3 0 0 10 & p4 10 depth 2 p5 0 depth 10 p6 10 0 10 p7 10 depth 10 ;建立模型的单元为放射性圆柱,其内部圆柱大小尺寸为5×unmy×8, 其总体的 ;具体的空间位置由p0,p1,p2,p3,p4,p5,p6,p7(它们的具 ;体意思参看generate 命令)决定。 plo crea qqq plo add surf red plo sho pause ; ;暂停 gen zone radcyl size 5 numy 8 10 & p0 0 0 0 p1 10 0 -10 p2 0 depth 0 p3 10 0 2 & p4 10 depth -10 p5 10 depth 2 p6 10 0 -4 p7 10 depth -4 ;建立模型的单元为放射性圆柱,其内部圆柱大小尺寸为5×unmy×8, 其中数值 ;10 的意义是扩展网格的个数。其总体的具体的空间位置由 ;p0,p1,p2,p3,p4,p5,p6,p7(它们的具体意思参看generate 命令) 决定。 plo crea qqq plo add surf red plo sho pause ; gen zone brick size 5 numy 10 & p0 0 0 -10 p1 10 0 -10 p2 0 depth -10 p3 0 0 2 p4 10 depth -10 & p5 0 depth -2 p6 2 0 -2 p7 2 depth -2 ;建立块体单元网格,其大小为5×unmy×10, ;空间位置由p0,p1,p2,p3,p4,p5,p6,p7 ;(它们的具体意思参看generate 命令)决定。 gen zone reflect orig 0 0 0 norm 1 0 0 ;对单元进行镜像,镜像面为以(0,0,0)为原点, ;以(1,0,0 )为法向向量。 plo crea qqq ;创建名为qqq plo add surf red ;表面颜色设置为红色 FLAC3D 2.10 plo show ;显示图形
Rotation: X: 20.000 Y: 0.000 Z: 10.000 Mag.: 1 Ang.: 22.500

Settings: Model Perspective 10:39:21 Sun May 08 2011 Center: X: 0.000e+000 Y: 5.000e+000 Z: 0.000e+000 Dist: 6.030e+001

Surface
Magfac = 0.000e+000

Itasca Consulting Group, Inc. Minneapolis, MN USA

56

FLAC3D

主要内容
?
? ? ?

?
? ? ? ? ? ?

FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并(联结)问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析

FLAC3D 3.00
Settings: Model Perspective 11:47:19 Sun May 08 2011 Center: X: 0.000e+000 Y: 1.000e+000 Z: 0.000e+000 Dist: 6.530e+000 Rotation: X: 300.000 Y: 0.000 Z: 0.000 Mag.: 0.8 Ang.: 22.500

Z X

Y

Surface
Magfac = 0.000e+000

Axes
Linestyle

Itasca Consulting Group, Inc. Minneapolis, MN USA

57

FLAC3D

Case-6:模型材料问题
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4 ;建立圆柱的网格单元,大小为4×5× 4 gen zone reflect norm 1,0,0 ;对单元进行镜像,法线为(1,0,0),原点为(0,0,0) gen zone reflect norm 0,0,1 ;对以上所有单元再进行镜像,法线为(0,0,1),原点为(0,0,0) plot Create Trench Add surface yellow Add axes black Show FLAC3D 3.00 Settings: Model Perspective ;显示模型格网 11:47:19 Sun May 08 2011 pause Center: Rotation: X: 0.000e+000 X: 300.000 model mohr Y: 1.000e+000 Y: 0.000 Z: 0.000e+000 Z: 0.000 ;采用摩尔-库仑模型 Dist: 6.530e+000 Mag.: 0.8 Ang.: 22.500 prop bulk 1.19e10 shear 1.1e10 ;为模型赋予参数,体积模量为1.19×1010 ,剪切模量为1.1×1010。 Surface Magfac = 0.000e+000 prop coh 2.72e5 fric 44 ten 2e5 Axes ; 内聚力为2.72×105,内摩擦角为44°,抗拉强的为2×105。 Linestyle fix x y z range y -.1 .1 ;边界y=0 为固定支座 fix x y z range y 1.9 2.1 ;边界y=2 为固定支座 ini yvel 1e-7 range y -.1 .1 ;在y=0的边界上,初始y方向上的速度为1×10?7 ini yvel -1e-7 range y 1.9 2.1 ;在y=2的边界上,初始y方向上的速度为-1×10?7 ; ini pp 1e5 hist gp ydisp 0,0,0 ;监测网格点(0,0,0)点的y 方向上的位移 hist zone syy 0,1,0 ;监测单元(0,1,0)点的y 方向上的应力 hist zone syy 1,1,0 ;监测单元(1,1,0)点的y 方向上的应力 step 3000 ;运行3000 步

Z X

Y

58

FLAC3D

Case-6:模型材料问题
FLAC3D 3.00
Step 3000 11:42:37 Sun May 08 2011
x10^-4 3.0

History
1 Y-Displacement Gp 1 Linestyle 1.000e-006 <-> 3.000e-004

2.5

Vs.
Step 1.000e+001 <-> 3.000e+003
2.0

1.5

1.0

0.5

监测网格点(0,0,0)点的y 方向上的位移
Itasca Consulting Group, Inc. Minneapolis, MN USA
0.5 1.0 1.5 x10^3 2.0 2.5 3.0

59

FLAC3D

Case-6:模型材料问题
监测单元(0,1,0)点的y 方向上的应力
FLAC3D 3.00
Step 3000 11:43:19 Sun May 08 2011
-0.5 x10^6

History
2 SYY Stress Zone 29 Linestyle -5.248e+006 <-> -4.916e+001
-1.0

Vs.
Step 1.000e+001 <-> 3.000e+003

-1.5

-2.0

-2.5

-3.0

-3.5

-4.0

-4.5

-5.0

Itasca Consulting Group, Inc. Minneapolis, MN USA

0.5

1.0

1.5 x10^3

2.0

2.5

3.0

60

FLAC3D

Case-6:模型材料问题
FLAC3D 3.00
Step 3000 11:43:50 Sun May 08 2011
x10^6

监测单元(1,1,0)点的y 方向上的应力
-0.2

History
3 SYY Stress Zone 172 Linestyle -1.374e+006 <-> -4.046e+002

Vs.
Step 1.000e+001 <-> 3.000e+003

-0.4

-0.6

-0.8

-1.0

-1.2

Itasca Consulting Group, Inc. Minneapolis, MN USA

0.5

1.0

1.5 x10^3

2.0

2.5

3.0

61

FLAC3D

主要内容
?
? ? ?

FLAC3D 3.00
Step 1626 Model Perspective 17:06:34 Thu May 26 2011 Center: X: 2.000e+000 Y: 2.000e+000 Z: 2.000e+000 Dist: 1.306e+001 Rotation: X: 20.000 Y: 0.000 Z: 70.000 Mag.: 0.8 Ang.: 22.500

?
? ? ? ? ? ?

FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并(联结)问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析
Surface Axes
Linestyle

Magfac = 0.000e+000

Z X Y

Itasca Consulting Group, Inc. Minneapolis, MN USA

62

FLAC3D

Case-7:Interface 合并(联结)问题
new ;新建 ; glued interface grid ;粘合交互网格 gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,2 ;建立块体网格单元,大小为4×4× 4,其具体位置由 p0,p1,p2,p3 的坐标决定。 gen zone brick size 8 8 4 p0 0,0,3 p1 4,0,3 p2 0,4,3 p3 0,0,5 ;建立块体网格单元,大小为8×8×4,其具体位置由 p0,p1,p2,p3 的坐标决定。 call plot.dat ;调用plot.dat 文件 pause ;暂停 inter 1 face range z 2.9,3.1 ;在z= 3 的面上建立交互单元 inter 1 prop kn 300e9 ks 300e9 tens 1e10 SBRATIO =1 ;赋予交界面参数,法向刚度为3.0×1011,剪切刚度为 3.0×1011, ;抗拉强度为1.0×1010,网格的比率为1 ini z add -1.0 range z 2.9,5.1 ;z 方向上在2.9 到5.1 的范围内,所有z 的坐标都减1 model elas ;材料模型为弹性模型 prop bulk 8e9 shear 5e9 ;材料的参数为,体积模量为8×109,剪切模量为5×109 fix z range z -.1 .1 ;固定边界z=0 ,为滚动支座 fix x range x -.1 .1 ;固定边界x=0 ,为滚动支座 fix x range x 3.9 4.1 ;固定边界x=4 ,为滚动支座 fix y range y -.1 .1 ;固定边界y=0 ,为滚动支座 fix y range y 3.9 4.1 ;固定边界y=4 ,为滚动支座 apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2 ;在z=4 的*面中,x(0,2),y(0,2)的范围内施加 z方向的应力为?1×106 hist unbal ;监测不*衡力 solve ;计算,一般计算到ratio=1×10?5为止 save inter.sav ;将文件保存为inter.sav FLAC3D 3.00 call plot.dat 1626 Model Perspective Step 08:04:56 Fri May ;调用plot.dat 文件27 2011 Center: Rotation: X: 2.000e+000 X: 20.000 pause Y: 2.000e+000 Y: 0.000 FLAC3D 3.00 Z: 2.000e+000 Z: 30.000 ;暂停 Dist: 1.306e+001 Mag.: 0.64
Ang.: 22.500
Rotation: X: 0.000 Y: 0.000 Z: 0.000 Mag.: 1 Ang.: 22.500

Settings: Model Perspective 08:15:45 Fri May 27 2011 Center: X: 2.000e+000 Y: 2.000e+000 Z: 2.500e+000 Dist: 1.583e+001

Surface
Magfac = 0.000e+000

Axes
Linestyle

Surface
Magfac = 0.000e+000

Axes
Linestyle

Z

Z Y X
X

Y

Itasca Consulting Group, Inc. Minneapolis, MN USA

Itasca Consulting Group, Inc.

63

FLAC3D

Case-7:Interface 合并(联结)问题
FLAC3D 3.00
Step 1626 Model Perspective 08:44:59 Fri May 27 2011 Center: X: 2.000e+000 Y: 2.000e+000 Z: 2.500e+000 Dist: 1.583e+001 Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 0.8 Ang.: 22.500

Surface
Magfac = 0.000e+000

Axes
Linestyle

Contour of Z-Displacement
Magfac = 0.000e+000 -2.2073e-004 to -2.0000e-004 -2.0000e-004 to -1.7500e-004 -1.7500e-004 to -1.5000e-004 -1.5000e-004 to -1.2500e-004 -1.2500e-004 to -1.0000e-004 -1.0000e-004 to -7.5000e-005 -7.5000e-005 to -5.0000e-005 -5.0000e-005 to -2.5000e-005 -2.5000e-005 to 0.0000e+000 0.0000e+000 to 0.0000e+000 Interval = 2.5e-005

Z Y X

History Location
Itasca Consulting Group, Inc. Minneapolis, MN USA

64

FLAC3D

Case-7:Interface 合并(联结)问题
FLAC3D 3.00 new ;新建文档 History ; attached grid Vs. ;连接网格 gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 Surface p3 0,0,2 Axes ;建立块体的网格单元,其大小为4×4×4,其空间位置由 p0,p1,p2,p3决定 gen zone brick size 8 8 4 p0 0,0,2 p1 4,0,2 p2 0,4,2 p3 0,0,4 ;建立块体的网格单元,其大小为8×8×4,其空间位置由 p0,p1,p2,p3决定 attach face range z 1.9 2.1 ;在z=2 的*面上连接交界面 model elas apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2 ;材料模型为弹性模型 ;在z=4 的*面中,x(0,2),y(0,2)的范围内施加z方向 prop bulk 8e9 shear 5e9 的应力为?1×106 ;体积模量为8×109,剪切模量为5×109 hist unbal fix z range z -.1 .1 ;监测不*衡力 ;固定边界z=0 ,为滚动支座 Solve fix x range x -.1 .1 ;计算,一般计算到ratio=1×10?5为止 ;固定边界x=0 ,为滚动支座 save att.sav fix x range x 3.9 4.1 ;将文件保存为att.sav ;固定边界x=4 ,为滚动支座 call plot.dat fix y range y -.1 .1 ;调用plot.dat 文件 ;固定边界y=0 ,为滚动支座 pause fix y range y 3.9 4.1 ;;暂停 ;固定边界y=4 ,为滚动支座
Step 1220 07:58:09 Fri May 27 2011
2.4

FLAC3D 3.00

x10^5

Settings: Model Perspective 08:21:59 Fri May 27 2011 Center: X: 2.000e+000 Y: 2.000e+000 Z: 2.000e+000 Dist: 1.306e+001

2.2

Rotation: X: 0.000 Y: 0.000 Z: 0.000 Mag.: 1 Ang.: 22.500

1 Max. Unbalanced Force Linestyle 1.268e+001 <-> 2.499e+005 Step 1.000e+001 <-> 1.220e+003

2.0

1.8

1.6

Magfac = 0.000e+000 Linestyle

1.4

1.2

1.0

0.8

0.6

0.4

Z

0.2

Y

X

Itasca Consulting Group, Inc. Minneapolis, MN USA

Itasca Consulting Group, Inc. Minneapolis, MN USA

0.2

0.4

0.6

0.8

1.0

1.2

x10^3

65

FLAC3D

Case-7:Interface 合并(联结)问题
new ;新建文本 FLAC3D 3.00 ; one grid Step 1649 Model Perspective 08:02:51 Fri May 27 2011 ;一个网格 Center: 0,4,0 p3 gen zone brick size 8 8 8 p0 0,0,0 p1 4,0,0 2.000e+000 Rotation: 0,0,4 p2 X: X: 30.000 Y: 2.000e+000 Y: 0.000 ;建立块体的网格单元,其大小为8×8×8,其空间位置由p0,p1,p2,p3决定 Z: 2.000e+000 Z: 30.000 model elas Dist: 1.306e+001 Mag.: 0.64 Ang.: 22.500 ;材料模型为弹性模型 prop bulk 8e9 shear 5e9 Surface ;材料的体积模量为8×109,剪切模量为5×109 Magfac = 0.000e+000 Axes fix z range z -.1 .1 Linestyle ;固定边界z=0 ,为滚动支座 fix x range x -.1 .1 ;固定边界x=0 ,为滚动支座 fix x range x 3.9 4.1 ;固定边界x=4 ,为滚动支座 fix y range y -.1 .1 ;固定边界y=0 ,为滚动支座 fix y range y 3.9 4.1 ;固定边界y=4 ,为滚动支座 apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2 Itasca Consulting Group, Inc. ;在z=4 的*面中,x(0,2),y(0,2)的范围内施加z方向的应力为?1×106 Minneapolis, MN USA hist unbal ;监测不*衡力 solve ;计算,一般计算到ratio=1×10?6为止 save noatt.sav ;将文件保存为noatt.sav call plot.dat ;调用plot.dat 文件

Y

Z X

66

FLAC3D

主要内容
?
? ? ?

?
? ? ? ? ? ?

FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并(联结)问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析

FLAC3D 3.00
Step 1114 Model Perspective 09:21:57 Fri May 27 2011 Center: X: 4.000e+000 Y: 4.000e+000 Z: 5.000e+000 Dist: 4.316e+001 Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500

Contour of SZZ
Magfac = 0.000e+000 Gradient Calculation -1.9998e+005 to -1.8000e+005 -1.8000e+005 to -1.6000e+005 -1.6000e+005 to -1.4000e+005 -1.4000e+005 to -1.2000e+005 -1.2000e+005 to -1.0000e+005 -1.0000e+005 to -8.0000e+004 -8.0000e+004 to -6.0000e+004 -6.0000e+004 to -4.0000e+004 -4.0000e+004 to -2.0000e+004 -2.0000e+004 to 0.0000e+000 0.0000e+000 to 3.0348e+000 Interval = 2.0e+004

Axes
Linestyle
Y

Z X

Itasca Consulting Group, Inc. Minneapolis, MN USA

67

FLAC3D

Case-8:初始条件问题
APPLY 命令可以在其后面选用关键词gradient,这时,我们就可以在一个指定的范围内线 性地改变应力或应变。在关键词gradient 后面的参数为gx , gy , gz ,它们可以用来说 明应力或者应变在x , y , z 上的变化趋势。应力或者应变线性改变的距离,是以全局坐标 系下的原点(0,0,0)为准的。例如:

对于这个公式的理解,我们最好用下面的这得例子加以解释: apply sxx -10e6 gradient 0 , 0, 1e5 range z -100 , 0

特别指出,我们可以应用此命令来模拟由重力引起的,随深度而有所增加的应力变化。 我们在说明INITIAL 命令时,确定应力场是随坐标而产生坡度变化的,并且重力加速 度的值已有SET gravity 命令给出。
68

FLAC3D

Case-8:初始条件问题
gen zone brick size 8 8 10 ratio 1.2 1 1 ;产生块体的网格单元,其到小为8×8×10,在x,y, z方向上的放大比率依 ;次为1.2 , 1 ,1 model mohr ;材料模型为摩尔-库仑模型 ini dens 2000 ;所有材料的初始密度为2000 prop bulk 2e8 shear 1e8 ;材料的体积模量为2×108,剪切模量为1×108。 prop fric 30 ;材料的内摩擦角为30° 。 fix x range x -.1 .1 ;固定边界x=0 ,为滚动支座 fix x range x 7.9 8.1 ;固定边界x=8 ,为滚动支座 fix y range y -.1 .1 ;固定边界y=0 ,为滚动支座 fix y range y 7.9 8.1 ;固定边界y=8 ,为滚动支座 fix z range z -.1 .1 ;固定边界z=0 ,为滚动支座 set grav 10 ;设置重力加速度为10 solve ;运算

save ini1.sav ;形成sav 文件,保存为ini1.sav plot create szz_contour ;创建文件名为szz_contour(z 方向的应力云图) plot set cent 4 4 5 ;视图的中心坐标为(4,4,5) plot set rot 20 0 30 ;视图的旋转角度为(20,0,30) plot set mag 1.0 ;视图的放大倍数为1.0 plot add cont szz out on shade on ;打开szz 的等值线图 plot add axes ;打开坐标轴线 plot show ;显示窗口 pause ;暂停
FLAC3D 3.00
Step 910 Model Perspective 09:06:59 Fri May 27 2011 Center: X: 4.000e+000 Y: 4.000e+000 Z: 5.000e+000 Dist: 4.316e+001 Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500

Contour of SZZ

Magfac = 0.000e+000 Gradient Calculation -2.1011e+005 to -2.0000e+005 -2.0000e+005 to -1.7500e+005 -1.7500e+005 to -1.5000e+005 -1.5000e+005 to -1.2500e+005 -1.2500e+005 to -1.0000e+005 -1.0000e+005 to -7.5000e+004 -7.5000e+004 to -5.0000e+004 -5.0000e+004 to -2.5000e+004 -2.5000e+004 to 0.0000e+000 0.0000e+000 to 1.5899e+003 Interval = 2.5e+004 Linestyle

Axes
Y

Z X

Itasca Consulting Group, Inc. Minneapolis, MN USA

69

FLAC3D

Case-8:初始条件问题
FLAC3D 3.00
Step 910 Model Perspective 09:06:59 Fri May 27 2011 Center: X: 4.000e+000 Y: 4.000e+000 Z: 5.000e+000 Dist: 4.316e+001 Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500

Contour of SZZ
Magfac = 0.000e+000 Gradient Calculation -2.1011e+005 to -2.0000e+005 -2.0000e+005 to -1.7500e+005 -1.7500e+005 to -1.5000e+005 -1.5000e+005 to -1.2500e+005 -1.2500e+005 to -1.0000e+005 -1.0000e+005 to -7.5000e+004 -7.5000e+004 to -5.0000e+004 -5.0000e+004 to -2.5000e+004 -2.5000e+004 to 0.0000e+000 0.0000e+000 to 1.5899e+003 Interval = 2.5e+004

Axes
Linestyle
Y

Z X

Itasca Consulting Group, Inc. Minneapolis, MN USA

70

FLAC3D

Case-8:初始条件问题
new ;新建文本 gen zone brick size 8 8 10 ratio 1.2 1 1 ;产生块体的网格单元,其到小为8×8×10,在x,y, ;z方向上的放大比率依次为1.2 , 1 ,1 model mohr ;材料模型为摩尔-库仑模型 ini dens 2000 ;所有材料的初始密度为2000 prop bulk 2e8 shear 1e8 ;材料的体积模量为2×108,剪切模量为1×108。 prop fric 30 ;材料的内摩擦角为30° 。 fix x range x -.1 .1 ;固定边界x=0 ,为滚动支座 fix x range x 7.9 8.1 ;固定边界x=8,为滚动支座 fix y range y -.1 .1 ;固定边界y=0 ,为滚动支座 fix y range y 7.9 8.1 ;固定边界y=8,为滚动支座 fix z range z -.1 .1 ;固定边界z=0 ,为滚动支座 ini szz -2.0e5 grad 0,0,2e4 ;施加初始应力,因为0 ;zz zz σ =σ + gx× x + gy× y + gz × z,所以, ;在此命令下2.0 105 2 104 zz σ = ? × + × × z ini sxx -1.5e5 grad 0,0,1.5e4 ;同理, 1.5 105 1.5 104 xx σ = ? × + × × z ini syy -1.5e5 grad 0,0,1.5e4 ;同理, 1.5 105 1.5 104 yy σ = ? × + × × z set grav 10 ;设置重力加速度为10 solve ;运算 save ini2.sav ;形成sav 文件,保存为ini2.sav plot create szz_contour ;创建文件名为szz_contour(z 方向的应力云图) plot set cent 4 4 5 ;视图的中心坐标为(4,4,5) plot set rot 20 0 30 ;视图的旋转角度为(20,0,30) plot set mag 1.0 ;视图的放大倍数为1.0 FLAC3D 3.00 plot add cont szz out on shade on ;打开szz 的云图 plot add axes ;打开坐标轴线 plot Contour of SZZ show ;显示窗口 pause ;暂停
Step 1 Model Perspective 09:09:39 Fri May 27 2011 Center: X: 4.000e+000 Y: 4.000e+000 Z: 5.000e+000 Dist: 4.316e+001 Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500 Magfac = 0.000e+000 Gradient Calculation -2.0000e+005 to -2.0000e+005 -2.0000e+005 to -1.7500e+005 -1.7500e+005 to -1.5000e+005 -1.5000e+005 to -1.2500e+005 -1.2500e+005 to -1.0000e+005 -1.0000e+005 to -7.5000e+004 -7.5000e+004 to -5.0000e+004 -5.0000e+004 to -2.5000e+004 -2.5000e+004 to 0.0000e+000 0.0000e+000 to 6.0594e-008 Interval = 2.5e+004 Linestyle
Y

Axes

Z X

71

FLAC3D

Case-8:初始条件问题
FLAC3D 3.00
Step 1 Model Perspective 09:09:39 Fri May 27 2011 Center: X: 4.000e+000 Y: 4.000e+000 Z: 5.000e+000 Dist: 4.316e+001 Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500

Contour of SZZ
Magfac = 0.000e+000 Gradient Calculation -2.0000e+005 to -2.0000e+005 -2.0000e+005 to -1.7500e+005 -1.7500e+005 to -1.5000e+005 -1.5000e+005 to -1.2500e+005 -1.2500e+005 to -1.0000e+005 -1.0000e+005 to -7.5000e+004 -7.5000e+004 to -5.0000e+004 -5.0000e+004 to -2.5000e+004 -2.5000e+004 to 0.0000e+000 0.0000e+000 to 6.0594e-008 Interval = 2.5e+004

Axes
Linestyle
Y

Z X

Itasca Consulting Group, Inc. Minneapolis, MN USA

72

FLAC3D

FLAC3D 3.00
Step 1114 Model Perspective 09:21:57 Fri May 27 2011

Case-8:初始条件问题
new ;新建文本 gen zone brick size 8 8 10 ratio 1.2 1 1 ;产生块体的网格单元,其到小为8×8×10,在x,y, ;z方向上的放大比率依 ;次为1.2 , 1 ,1 model mohr ;材料模型为摩尔-库仑模型 ini dens 2000 ;所有材料的初始密度为2000 prop bulk 2e8 shear 1e8 ;材料的体积模量为2×108,剪切模量为1×108。 prop fric 30 ;材料的内摩擦角为30° 。 fix x range x -.1 .1 ;固定边界x=0 ,为滚动支座 fix x range x 7.9 8.1 ;固定边界x=8,为滚动支座 fix y range y -.1 .1 ;固定边界y=0 ,为滚动支座 fix y range y 7.9 8.1 ;固定边界y=8,为滚动支座 fix z range z -.1 .1 ;固定边界z=0 ,为滚动支座 set grav 10 ;设置重力加速度为10

Center: X: 4.000e+000 Y: 4.000e+000 Z: 5.000e+000 Dist: 4.316e+001

Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500

Contour of SZZ
Magfac = 0.000e+000 Gradient Calculation -1.9998e+005 to -1.8000e+005 -1.8000e+005 to -1.6000e+005 -1.6000e+005 to -1.4000e+005 -1.4000e+005 to -1.2000e+005 -1.2000e+005 to -1.0000e+005 -1.0000e+005 to -8.0000e+004 -8.0000e+004 to -6.0000e+004 -6.0000e+004 to -4.0000e+004 -4.0000e+004 to -2.0000e+004 -2.0000e+004 to 0.0000e+000 0.0000e+000 to 3.0348e+000 Interval = 2.0e+004

Axes
Linestyle
Y

Z X

solve elastic Itasca Consulting Group, Inc. ;进行弹性计算 Minneapolis, MN USA save ini3.sav ;形成sav 文件,保存为ini3.sav plot create szz_contour ;创建文件名为szz_contour(z 方向的应力云图) plot set cent 4 4 5 ;视图的中心坐标为(4,4,5) plot set rot 20 0 30 ;视图的旋转角度为(20,0,30) plot set mag 1.0 ;视图的放大倍数为1.0 plot add cont szz out on shade on ;打开szz 的云图 plot add axes ;打开坐标轴线 plot show ;显示窗口 73

FLAC3D

Case-8:初始条件问题
FLAC3D 3.00
Step 1114 Model Perspective 09:21:57 Fri May 27 2011 Center: X: 4.000e+000 Y: 4.000e+000 Z: 5.000e+000 Dist: 4.316e+001 Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500

Contour of SZZ
Magfac = 0.000e+000 Gradient Calculation -1.9998e+005 to -1.8000e+005 -1.8000e+005 to -1.6000e+005 -1.6000e+005 to -1.4000e+005 -1.4000e+005 to -1.2000e+005 -1.2000e+005 to -1.0000e+005 -1.0000e+005 to -8.0000e+004 -8.0000e+004 to -6.0000e+004 -6.0000e+004 to -4.0000e+004 -4.0000e+004 to -2.0000e+004 -2.0000e+004 to 0.0000e+000 0.0000e+000 to 3.0348e+000 Interval = 2.0e+004

Axes
Linestyle
Y

Z X

Itasca Consulting Group, Inc. Minneapolis, MN USA

74

FLAC3D

主要内容
?
? ? ?

FLAC3D 3.00
Step 2836 Model Perspective 11:31:06 Fri May 27 2011 Center: X: 1.000e+001 Y: 5.000e-001 Z: 5.000e+000 Dist: 5.580e+001 Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500

?
? ? ? ? ? ?

FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并(联结)问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析

Surface
Magfac = 0.000e+000 Linestyle

Axes

Z Y

X

Itasca Consulting Group, Inc. Minneapolis, MN USA

75

FLAC3D

Case-9:破坏问题分析
gen zone brick size 20 1 10 ;建立块体网格,其大小为20×1×10。 model mohr ;材料模型为摩尔-库仑模型 prop dens 2000 bul 2.e8 shea 1.e8 cohesion 0.0 ;材料的性质: 密度为 2000, 体积模量为2.0×108,剪切 模量为1.0×108, ;粘聚力为0 prop friction 30. dilation 0. tension 0 ;内摩擦角为30° ,剪胀角为0,抗拉强度为0 fix x range x -.1 .1 ;固定边界x=0 ,为滚动支座 save fail1.sav fix x range x 19.9 20.1 ;保存为fail1.sav ;固定边界x=20 ,为滚动支座 free x range x -.1 .1 z 6.9 10.1 fix x y z range z -.1 .1 ;放松x=0 *面上,z=7,10 这一部分在x 方向的约束 ;固定边界z=0 ,为固端支座 step 2000 fix y ;运算2000 个时步 ;固定所有单元y 方向的运动 save fail2.sav set gravity 0 0 -10 ;保存为fail2.sav ;设置重力加速度为(0,0,-10) plot hist unbal Create Trench ;监测不*衡力 Add surface yellow solve Add axes black ;运算

Show

76

FLAC3D

Case-9:破坏问题分析
FLAC3D 3.00
Step 2836 Model Perspective 11:53:10 Fri May 27 2011 Center: X: 1.000e+001 Y: 5.000e-001 Z: 5.000e+000 Dist: 5.580e+001 Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500

土体的速度矢量

Surface
Magfac = 0.000e+000

Axes
Linestyle

Velocity
Maximum = 8.161e-004 Linestyle

Z Y

X

Itasca Consulting Group, Inc. Minneapolis, MN USA

77

FLAC3D

Case-9:破坏问题分析
FLAC3D 3.00
Step 2836 Model Perspective 11:57:29 Fri May 27 2011 Center: X: 1.000e+001 Y: 5.000e-001 Z: 5.000e+000 Dist: 5.580e+001 Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500

土体的剪切应力等值线图

Surface
Magfac = 0.000e+000

Axes
Linestyle

Contour of Shear Strain Increment
Magfac = 0.000e+000 Gradient Calculation 9.5205e-008 to 2.5000e-002 2.5000e-002 to 5.0000e-002 5.0000e-002 to 7.5000e-002 7.5000e-002 to 1.0000e-001 1.0000e-001 to 1.2500e-001 1.2500e-001 to 1.5000e-001 1.5000e-001 to 1.7500e-001 1.7500e-001 to 2.0000e-001 2.0000e-001 to 2.2500e-001 2.2500e-001 to 2.4193e-001 Interval = 2.5e-002 Itasca Consulting Group, Inc. Minneapolis, MN USA

Z Y

X

78

FLAC3D

Case-9:破坏问题分析
FLAC3D 3.00
Step 2836 Model Perspective 12:01:54 Fri May 27 2011 Center: X: 1.000e+001 Y: 5.000e-001 Z: 5.000e+000 Dist: 5.580e+001 Rotation: X: 20.000 Y: 0.000 Z: 30.000 Mag.: 1 Ang.: 22.500

土体的塑性状态

Block State
None shear-n shear-p shear-n shear-p tension-p shear-p shear-p tension-p

Itasca Consulting Group, Inc. Minneapolis, MN USA

79


相关推荐

最新更新

猜你喜欢