UDEC学习整理资料
UDEC3.1资料
通用离散元用户指导（U D E C 3.1）北京科技大学2004.9目录1 引言 (1)1.1 总论 (1)1.2 与其他方法的比较 (2)1.3 一般特性 (2)1.4 应用领域 (3)2 开始启动 (4)2.1 安装和启动程序 (4)2.1.7 内存赋值 (4)2.1.9 运行UDEC (5)2.1.10 安装测试程序 (5)2.2 简单演示－通用命令的应用 (5)2.3 概念与术语 (6)2.4 UDEC模型：初始块体的划分 (8)2.5 命令语法 (9)2.6 UDEC应用基础 (10)2.6.1 块体划分 (10)2.6.2 指定材料模型 (16)2.6.2.1 块体模型 (16)2.6.2.2 节理模型 (17)2.6.3 施加边界条件和初始条件 (19)2.6.4 迭代为初始平衡 (21)2.6.5 进行改变和分析 (24)2.6.6 保存或恢复计算状态 (25)2.6.7 简单分析的总结 (25)2.8 系统单位 (26)3 用UDEC求解问题 (27)3.1 一般性研究 (27)3.1.1 第1步：定义分析模型的对象 (28)3.1.2 第2步：产生物理系统的概念图形 (28)3.1.3 第3步：建造和运行简单的理想模型 (28)3.1.4 第4步：综合特定问题的数据 (29)3.1.5 第5步：准备一系列详细的运行模型 (29)3.1.6 第6步：进行模型计算 (29)3.1.7 第7步：提供结果和解释 (30)3.2 产生模型 (30)3.2.1 确定UDEC模型合适的计算范围 (30)3.2.2 产生节理 (32)3.2.2.1 统计节理组生成器 (32)3.2.2.2 VORONOI多边形生成器 (34)3.2.2.3 例子 (34)3.2.3 产生内部边界形状 (35)3.3 变形块体和刚体的选择 (38)3.4 边界条件 (42)3.4.1 应力边界 (42)3.4.1.1 施加应力梯度 (43)3.4.1.2 改变边界应力 (44)3.4.1.3 打印和绘图 (44)3.4.1.4 提示和建议 (45)3.4.2 位移边界 (46)3.4.3 真实边界－选择合理类型 (46)3.4.4 人工边界 (46)3.4.4.1 对称轴 (46)3.4.4.2 截取边界 (46)3.4.4.3 边界元边界 (49)3.5 初始条件 (50)3.5.1 在均匀介质中的均匀应力：无重力 (50)3.5.2 无节理介质中具有梯度变化的应力：均匀材料 (51)3.5.3 无节理介质中具有梯度变化的应力：非均匀材料 (51)3.5.4 具有非均匀单元的密实模型 (52)3.5.5 随模型变化的初始应力 (53)3.5.6 节理化介质的应力 (54)3.5.7 绘制应力等值线图 (55)3.6 加载与施工模拟 (57)3.7 选择本构模型 (62)3.7.1 变形块体材料模型 (63)3.7.2 节理材料模型 (64)3.7.3 合理模型的选择 (65)3.8 材料性质 (71)3.8.1 岩块性质 (71)3.8.1.1 质量密度 (71)3.8.1.2 基本变形性质 (71)3.8.1.3 基本强度性质 (72)3.8.1.4 峰后效应 (73)3.8.1.5 现场性质参数的外延 (77)3.8.2 节理性质 (80)3.9 提示和建议 (81)3.9.1 节理几何形状的选择 (81)3.9.2 设计模型 (81)3.9.3 检查模型运行时间 (82)3.9.4 对允许时间的影响 (82)3.9.5 单元密度的考虑 (83)3.9.6 检查模型响应 (83)3.9.7 检查块体接触 (83)3.9.8 应用体积模量和剪切模量 (83)3.9.9 选择阻尼 (84)3.9.10 给块体和节理模型指定模型和赋值 (84)3.9.11 避免圆角误差 (85)3.9.12 接触嵌入 (85)3.9.13 非联结块体 (86)3.9.14 初始化变量 (86)3.9.15 确定坍塌荷载 (86)3.9.16 确定安全系数 (86)3.10 解释 (88)3.10.1 不平衡力 (88)3.10.2 块体/网格结点的速度 (88)3.10.3 块体破坏的塑性指标 (89)3.11 模拟方法 (89)3.11.1 有限数据系统模拟 (89)3.11.2 混沌系统的模拟 (90)3.11.3 局部化、物理的不稳定性和应力路径 (91)1 引言1.1 总论通用离散元程序（UDEC，Universal Distinct Element Code）是一个处理不连续介质的二维离散元程序。
《udec版本版》课件
软件下载与安装
从udec官网下载最新版本，按照安装向导逐步完成安装过程。
配置参数
根据实际需求，配置udec版本的相关参数，如端口号、数据库连 接等。
udec版本的团队协同工作流程
团队协作平台搭建
01
建立团队协作平台，包括项目管理、任务分配、进度跟踪等功
面对激烈的市场竞争，udec版本需要不断提升产品品质和服务 水平，以赢得用户信任和支持。
新兴市场拓展
拓展新兴市场，发掘潜在用户需求，为udec版本的发展提供更 多机遇。
udec版本的未来发展方向
拓展应用领域
将udec版本的应用领域不断拓展，满足更多行业和用户的需求。
提升用户体验
持续优化udec版本的用户界面和功能，提高用户满意度和忠诚度 。
开发效率。
持续集成
udec支持持续集成功能，能够将 项目的构建、测试和部署等环节 集成在一起，实现自动化和持续
化的软件开发流程。
自动化测试
udec支持自动化测试功能，能够 自动执行测试用例并生成测试报 告，帮助开发者及时发现和修复
问题。
03
udec版本的实施 与部署
udec版本的安装与配置
安装环境要求
代码导航
代码重构
udec支持代码重构功能，能够帮助开 发者优化和改进代码结构，提高代码 的可维护性和可读性。
udec提供了代码导航功能，方便用户 快速找到和跳转到代码的特定部分， 提高代码阅读和编辑的便捷性。
自动化构建功能
自动化构建
udec支持自动化构建功能，能够 根据项目的需求自动完成项目的 编译、打包和部署等任务，提高
(完整word版)Udec计算算例
1.3 UDEC算例1.3.1工程概况某隧道位于一包含高角度连续节理岩体内，节理倾角为50度，平均间距为7m，隧道为一半径为9m的圆形隧道。
贯穿于开挖面内的一垂直断层，在隧道拱顶形成了一个三角楔形体。
本算例使用UDEC的结构单元逻辑来模拟喷射混凝土和锚杆联合支护的圆形隧道开挖问题。
1.3.2构建模型隧道埋深451m，为半径9m圆形隧道，本次计算模型左右边界取41m，隧道至上下边界也取41m。
总的来说，模型长100m，宽100m。
计算模型如图1-1所示。
图1-1 UDEC计算模型1.3.3计算参数在包含高倾角节理和垂直断层的岩体内进行圆形开挖的UDEC模型岩体、节理和断层参数如下所示：表1-1 完整岩石物理力学参数密度Dens （kg/m3）体积模量K（Gpa）剪切模量G（Gpa）2500 1.5 0.6表1-2 节理、断层物理力学参数表1-3 喷射混凝土物理力学参数表1-4 岩体和喷射混凝土接触面物理力学参数表1-5 锚杆物理力学参数作为演示的目的，隧道开挖和支护是瞬时发生的。
本算例共两种支护分析被计算：第一，只施加喷射混凝土衬砌；第二，喷射混凝土和锚杆联合提供支护。
为了在第二种支护情况分析中可以更清晰的看到锚杆提供的支护，算例采用喷射混凝土的抗压强度被设置成一个很低的值，且厚度仅取为10cm。
1.3.4模拟步骤1.建立模型在UDEC中输入以下命令可建立隧道结构模型及边界。
如图1-2所示。
newround 0.1block -50,-50 -50,50 50,50 50,-50 jset -50,0 100,0 0,0 7,0 ；刷新UDEC窗口，重新调用一个新程序；块与块之间的圆角半径，必须小于块体最小边的1/2 ；建立模型框架；设置节理crack -6 -50 -6 50 tunnel (0,0) 9,16 del range area 0.08 gen edge 10 ；设置断层；模拟开挖隧道边界；删除面积小于0.08的块体；自动划分单元，单元最大边长不超过10图1-2 初步模型图2.设置单元属性和材料特征在UDEC中输入如下命令设置单元属性和材料特征。
学习资料整理.doc
学习资料整理—：TD-SCDMA系统概述（重点）第三代公众移动通信系统的工作频段：（一）主要工作频段：频分双工（FDD）方式：1920-1980MHZ / 2110-2170MHZo 时分双工（TDD）方式：1880-1920MHz> 2010-2025MHZ （二）补充工作频率：频分双工（FDD）方式：1755-1785MHz/1850-1880MHz时分双工（TDD）方式：2300-2400MHZo （三）卫星移动通信系统工作频段：1980-2010MHZ / 2170-2200MHZTD-SCDMA优势：频谱利用率高不需成对的频谱，能够满足未来扩展需求，为频谱分眠带来极大的灵活性相对于FDD运营商，TDD运营商频谱获取成本低，同时在业务方面，提高语音和非对称数据应用的频谱效率TD系统分配非对称上下行传输，经济高效地支持互联网接入业务结合智能天线技术，可以提供快速精确定位业务（LCS）TD-SCDMA基本原理：时分多址——在时间轴上，上行和下行分开，实现了TDD模式。
频分多址——TDD模式反映在频率上，是上行下行共用一个频点，节省了带宽。
在频率轴上，不同频点的载波可以共存。
码分多址——在能最轴上，每个频点的每个时隙可以容纳16个码道。
空分多址——通过使用智能天线技术，针对不同的用户使用不同的赋形波束覆盖。
智能天线由于采用了波束赋形技术，可以有效的降低干扰，提高系统的容量。
系统网络接曰：无线接口从协议结构上可以划分为三层：物理层（L1）数据链路层（L2）网络层（L3）L2和L3划分为控制平面（C-平面）和用户平面（U-平面）。
RLC和MAC之间的业务接入点（SAP）提供逻辑信道，物理层和MAC之间的SAP提供传输信道。
RRC与下层的PDCP、BMC、RLC和物理层之间都有连接，用以对这些实体的内部控制和参数配置UE只监听PICH信道和接收广播信道信息。
空闲模式UE由非接入层标识，如IMSI、TMSI和P-TMSI,此时在UTRAN中没有单独的空闲模式的UE信息。
udec学习总结
1.改变水压，可以通过将渗透系数设为0来实现。
即prop jmat=2 azero=1e-4 ares=0.5e-4 jperm=0
Udec学习与总结
1 round的设置，
Round的值为模型中最小块体最小边长的1%。
太大会引起模型不正确反映。
2 要删除一部分区域，最好在建模前将这一区域划分出来。
3 change mat=1等设置材料号时，最好是一层一个号，相同材料同一个材料号也不是不可以，但是计算时运算步会变大，也不利于修改材料属性、
4，bulk和shear控制模型材料的总体软硬程度，bulk调大模型总体变硬，bulk调小，模型总体变软；内聚力控制模型的抗剪程度，udec中的塑性区变形大部分跟它有关系。
Ten控制抗拉强度，拉变形与它有关。
5 jkn对模型节理的竖向变形有关，jks与横向层理有关。
udec命令总结精华-正宗
《Udec 命令总结》精华版1. 安装、打开、保存(1) 安装：①执行Udec 3.1→将Crack文件中的内容替换；②复制Udec.exe,粘贴为快捷方式→属性，目标，加入空格 256【开始内存8M，将内存改为256】，可以复制快捷方式至桌面或硬盘。
(2) 打开：① 直接在udec> 命令行输入；② 写好程序，udec>命令输入call，然后将*.txt文件拖入命令行，执行(3) 保存：输入 save d:\kaicai.sav, 调用命令rest，将kaicai.sav 拖入命令行【！！！文件名最好不用汉字，有时候不识别】注释：如果保存为save d:\111\kaicai.sav, 注意其中111文件必须提前建好，否则无法保存或者保存错误(4) 操作：Ctrl+Z选中图像可以放大，Ctrl+Z 双击复原，屏幕中会出现十字叉，按住鼠标左键不放，移动光标直到你满意的窗口为止；pause暂停，此时可以察看任何信息；continue 继续调用下面程序段。
ESC可以随时进行停止，但不能继续；英文分号 ; 表示注释不运行命令。
2. 基本命令2.1 基本设置Udec> n 【new刷新窗口，从新调用一个程序，修改后的*.txt文件必须输入n，重新运行文本文件】Udec> title 【或heading代表标题，后面紧跟标题的名称。
如：hang dao mo ni 或济宁三号井围岩变形破坏规律研究】Udec> round d 【“圆角”命令，Udec中所有的块体都有圆角，目的是为防止块体悬挂在有棱角的节点上，由于块体悬挂将产生应力集中。
d指块体与块体之间的圆角半径，默认值是0.5，其值要求小于模型中最小块体的最短那条边长的二分之一，最大圆角长度不能超过块体平均棱长的1%。
在block命令前指定圆角长度。
如：round 0.05】图1 10×10块体圆角Udec>set ovtol 0.5 【此命令是指层与层之间的嵌入厚度，当提示为“overlap too large”时就需要修改此值更大一些，可以显示设计的块体，plot overlap！嵌入太大的原因可能为块体强度太小】Udec>set log onprintset log off 【命令用于导出数据，可以将数据导出至硬盘f:\】命令：set log f:\yuanyan.logprint pline 1 syy 【记录测线1应力】print pline 2 ydisset log off;设置观测线set pline 1 x1,y1 x2,y2 n (n--观测线分的段数)set pline 25,20 300,20 10 ;定义测线，起始点坐标，测线分段2.2 图形划分;块体命令Udec> block x1,y1 x2,y2 x3,y3 x4,y4【建立模型框架，四个坐标角点必须按“顺时针”排列，也可以为五个坐标点。
UDEC模拟[技巧]
![UDEC模拟[技巧]](https://img.taocdn.com/s3/m/aca4360ca9114431b90d6c85ec3a87c240288aee.png)
1 模型的建立建立数学模型是数值模拟工作的首要任务, 模型建立正确与否, 是能否获得符合实际计算结果的前提, 模型的设计, 必须遵循下列原则:采动覆岩移动的影响因素很多, 模型的设计,必须突出影响采动覆岩移动的主要因素, 并尽可能多地考虑其它重要因素。
模型是由实体简化的, 但应不失一般性。
模型的设计, 必须能很好地反映材料的物理力学特性,如材料的均匀性, 弱面影响及非线性等。
地下工程实际上是半无限域问题, 但数值模拟只能是在有限的范围内进行。
因此, 模型的设计,必须考虑其边界效应, 选择适当的边界条件。
任何地下工程, 也都是一个时空问题, 采动围岩移动也是如此。
因此, 模型的设计,必须能体现工作面的推进与接续, 能体现出覆岩冒落、底板膨胀鼓起及变形移动的时间过程。
模型的设计, 应尽可能便于数值模拟计算, 在模型范围及受力分析方面, 既要满足弹塑性理论对应力分析的基本要求, 又要顾及现有计算机的容量。
2 模型的基本参数各岩层物理力学参数按表2.1选取，表中抗拉强度、泊松比参考附近矿区岩层实际参数，由于该矿并没有各岩层粘聚力和摩擦角等参数，粘聚力、摩擦角和弹性模量按该岩性岩体平均参数选取，体积模量和剪切模量由泊松比和弹性模量按公式计算得出。
νE K=3（1-2） νEG=2（1+）式中：K 为岩体体积模量；G 为岩体剪切模量；E 为岩体的弹性模量；ν为岩体的泊松比。
表2.1 模型中采用的岩体物理力学参数岩层名称岩层厚度/m体积模量/GPa 剪切模量 /GPa 抗拉强度 /MPa 粘聚力 /MPa 内摩擦 角/° 砂质泥岩或粉砂岩15.63 3.125 2.542 0.8 2.5 35 11煤10.94 2.381 1.163 0.65 1.3 32.9 砂质泥岩或粉砂岩2 43.67 3.571 2.459 0.74 2.5 35 砂质泥岩或粉砂岩37.92 6.667 2.222 0.76 2.5 35 9煤6.4 2.381 1.163 0.65 1.3 32.9 砂质泥岩或粉砂岩421.83 2.857 2.609 1 2.5 35 8煤3.52 2.381 1.163 0.65 1.3 32.9 砂质泥岩或粉砂岩527.77 10 2.143 0.72 2.5 35 6煤6.35 2.381 1.163 0.65 2.3 32.9 砂质泥岩或粉砂岩611.42 10 2.1430.88 2.5 35 砾岩502.1351.6680.82.234依据工作面的地质条件, 建立图1所示的数值计算模型。
UDEC3.0中文手册33-63
－33－
特殊 的几何 参数 对力学 特性 产生的影响 可以进行 定量 描述 。 同时，在获 得 现 场观测的节 理模 式情况下，采 用人工模 式与观 测模式 相匹配方法决定 生成的参数。 一个节理组可以通过 8 个生成参数表述： 4 个几何 参数的 均值和 4 个随机参数的均 方差 。由 JSET 命令给出 的参数如下： JSET a m、a d 、t m、t d 、g m、g d 、s m、s d < x0、y 0 >< ad 0 > 在此， a －节理与 x 轴的 夹角；
(K + 4 3 G max ∆z min
在式 中， K 、 G － 分别 是体积模量和剪 切模量；
－34－
理论与背景 第 1.1.3 节中的 range 关键词的描述。在大 部分情况下， 通过 jregion n 关键 词定义各种 区域，在此 n 是 指 定 JREGION 命令的参考序号 id。 JREGION 定义 了一个 凸多边形区域来 限制节理组的生成范围 。该命令参数如 下： JREGION id n x1 y1 x2 y2 x3 y3 x4 y4 <delete> 每 一个节理 区域是通过 id 序号识 别。区域的坐标 按顺时 针方向定义了节理产生的边 界。如果给出选择 的关键词 delete ，在此之前 由 JSET、VORONOI 或 CRACK 命令生成 的节理将 全部 被 删除。 这 可避免当指 定多 个节理 区域情况下 ， 被相邻区域 产生的节理所 切割。 其他相对于 JSET 命令（即 mat n ）的关键词 range 也能 用来限制 块体被 切割。如果 对 JSET 命令没 有指定范围，将产生 遍及整个区域 的节理。 JSET 生 成器也可能通过设定很大的 t m 和 s m 值用来产生贯穿 UDEC 块体的一条节 理。请记住， JSET 能够产生不连续节理。不连续节理完全处在块体里边不能看见。当 一个变形块体开始生成单元或当一个刚体模型开始运行时，所有的节理段将被删除 。 3.2.2.2 VORONOI 多边形生成器 VORONOI 生成器产生 随机大小的 多边形 块体。 UDEC 模型中的一个或 多个块体能 够 在 细 分 成 任意 大小 的 Voronoi 子 块 。 该 节理生 成器 对 模 拟 裂缝扩展 是 有 用的。当 Voronoi 块体间的节理强度 被超过 时将发生 断裂。 VORONOI 命令具有 下列形式： VORONOI edge l < iterations n> <round v> < range … > 对 于 Voronoi 多边形 指 定平 均棱 长。 多边形 具 有 随机尺寸， 但具 有 平 均棱 长 l 。 Voronoi 块体 的 大小尺寸 能够通过增加 迭代次 数 变得 均匀缺省 值 n=5 。也 能 指 定 圆 角长 度。 圆角长 度， v 必须 至少小 于块体棱长 l 的 20 倍。 当应用 VORONOI 命令时 ，多边形区域应当略大 于被细 分的块体区域。 这将抑制边 界影响。Range 关键词在产生 Voronoi 块体区域的用 法，与 JSET 命令具有 相同的方 式。 Voronoi 算法 根据 随机 分布 点 从多边形 区域 开始。然后允 许内 部点移 动。迭代 过程 运 动到这 些 点。 迭代 步 越 高 ， 点 间距 越均匀 。接 下来，所有 点 产生 三角 形。 最 后， 通过 做具 有公共 边所有 三角形的 垂直平分线，生成 Voronoi 多边形。多边形在所 指定区域被 其边 界所截 取。 3.2.2.3 例子 为 了说明节理的生成过程 ，下面给出了产生节理的 几个例子。 例 3.3 四组 规则节理组 new ro 0.01 blo 0 0 0 20 20 20 20 0
udec讲义_外文
ST041 - Ecole des Mines de Nancy - Thierry Verdel
plot his 0
plot his 1
plot his 2
ST041 - Ecole des Mines de Nancy - Thierry Verdel
Fichier de commandes
udec讲义外文udec命令疯狂android讲义疯狂java讲义上海外文书店纯爱讲义外文文献数据库外文书店费曼物理学讲义外文出版社外文局
Méthodes des éléments distincts
Thierry Verdel
ST041 - Ecole des Mines de Nancy - Thierry Verdel
ST041 - Ecole des Mines de Nancy - Thierry Verdel
ST041 - Ecole des Mines de Nancy - Thierry Verdel
Exercices
ST041 - Ecole des Mines de Nancy - Thierry Verdel
Poutre en flexion sous une charge uniforme
Comparaison avec une solution analytique
ST041 - Ecole des Mines de Nancy - Thierry Verdel
Fichier de commandes
ST041 - Ecole des Mines de Nancy - Thierry Verdel
Equilibre initial
ST041 - Ecole des Mines de Nancy - Thierry Verdel
UDEC3.0中文手册64-101
图 3.34 直剪试验模型 首先， 施加具有代表性的 10MPa 法向应 力。然后 ，在块体端部施 加水平 应力产生剪
－67－
切位 移。此值也具有代表性的。为了说明这 个目的， 仅施 加小于 1mm 的剪切 位移。 采用 FISH 函数（ av_atr ） 记录了节理的 平均法向和 剪切 应力、法 向和剪切位移。基 于此 记录结果，我们能够确 定不同模型的 峰值和残余剪切强度以及剪胀特性。例 3.22 包 含了应用库仑滑移模型进行 该试验 模拟的数 据文件。 图 3.35 和 图 3.36 分别给出了节理的平 均剪应力 与 位移 的关心曲 线和平 均法向 位移 对 剪切位 移 的 图 形。 这 两个 图形 表明 ：对于给定的模型性质和 边界条件 ，节理 已经 发 生 滑移 。在峰 值剪切强 度接近 6MPa 前，图 3.35 中的 加载 斜率为 线性。正 如在图 3.36 所显 示的，当节理剪切位 移约为 0.15mm 时 ， 不能抗 剪 ，节理开始 剪胀。其剪 胀一致 延 续到 其极限剪切 位移（ zdilation=0.4mm）。最大 的平均剪 胀接近 于 0.027mm。
面接触 位移弱化的面接触 连续屈服
UDEC经典学习总结
UDEC经典学习总结1.把图形保存下来，能在AUTOCAD(图⽚)中打开、编辑plot block cable red supp ye stru bl;显⽰块体、锚杆（红⾊）、⽀架（黄⾊）、梁（蓝⾊）的图形set plot dxf(jpg) 256；设置图形为256⾊set out c:\ss.dxf(jpg) (set out 1.dxf 这个⽂件保存在当前带数值模拟的⽂件夹⾥)copy c:\ss.dxf(jpg) ( copy 1.dxf 这个⽂件保存在当前带数值模拟的⽂件夹⾥)；把这图形以ss.dxf⽂件保存在C磁盘下2.把数据导出set log onprint pline 1 ydis;把pline 1的y⽅向位移的数据导出到UDEC⽬录下udec.txt⽂件中set log off3.plot block stress ;显⽰块体的应⼒plot block dis ;显⽰块体的位移(有x和y⽅向)plot block pl ;显⽰块体的塑性区(plastic)save xx.save ;保存计算结果res xx.save ;调⽤4.液压⽀架的命令supp xc yc wid l seg n mat j;(xc,yc)是指中⼼点坐标,l指⽀架的宽度,n指分段数,mat j指⽀柱材料性质为j prop mat j sup_kn -1 ;sup_kn指⽀架的刚度,-1与表1相对应del range x1 x2 y1 y2 ;挖掘范围(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)的块体,由⽀架⽀撑table 1 0 42.e6 0.05 5.0e6 0.1 6.0e6;表1表⽰的是液压⽀架的(P-DS)特性曲线5 巷道施⼯中断⾯加梯⼦梁的命令(见图1)stru gen xc yc np 100 fa a thetra b mat=16 thick=0.2 ;100 指分100段,a b 指⾓度prop mat=16 st_ymod=13.5e9 st_prat=0.14 st_den=7800prop mat=16 st_yield=6e7 st_yresid=6e7 st_ycomp=2.5e7prop mat=16 if_kn=1.35e9 if_ks=1.35e9 if_tens=0 if_fric=18 if_coh=0;interface-界⾯例⼦讲解:建模(以⽶为单位,;后为解释部分)round=0.1 ;⽅块的圆⾓块半径为0.1⽶set ovtol=1.0 ;块体与块体之间相互嵌⼊量最⼤值为1⽶bl 0,0 0,26 50,26 50,0 ;在(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)⽣成块体范围crack 0,25 50,25 ;在(x1, y1)(x2 ,y2)两点间画直线jregion id 1 0,0 0,10 50,10 50,0 delete ;删除此区域的块体jset 90,0 2.5,0 2.5,0 3,0 0,0 range jregion 1 ;jset A,0 a,0 b,0 c,0 x0,y0 range jregion 1jset 90,0 2.5,0 2.5,0 3,0 1.5,2.5 range jregion 1 可画成列的线段(见图2)jset 0,0 26,0 0,0 2.5,0 0,0 range jregion 1pa ;当程度运⾏⾄此时暂停,可以看看你所建的部分模型,⽤continue继续运⾏下⾯部分save t.save ;建模保存在t.save中岩层赋属性及原岩⼒平衡计算res t.save ;调⽤已建好的模型gen quad 12 ;定义块体最⼤变形,若没有此语句,刚所有块体均为刚性块体zone model mo range 0,180 0,45.2 ;在范围(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)的块体符合库仑准则change jcons=2 range 0,180 0,45.2 ;节理⾯间接触-coulmb滑移;下⾯主要是讲岩层赋属性change （范围） mat=1 range reg 0,0 0,3 180,3 180,0 ;岩层1的范围change mat=2 range reg 0,3 0,4.5 180,4.5 180,3change mat=3 range reg 0,4.5 0,10.2 180,10.2 180,4.5change mat=4 range reg 0,10.2 0,11.2 180,11.2 180 10.2prop （赋值的意思） mat=1 dens=2500 ;岩层1的密度prop mat=2 dens=2500prop mat=3 dens=1300zone （块体） k=2e10, g=1.1e10, fric=30,coh=2e6,ten=4e6 range mat=1 ;岩层1的块体⼒学参数(见表1)zone k=8.7e9, g=4.2e9, fric=25,coh=1.5e6,ten=1.5e6 range mat=2zone k=13.05e9, g=6.3e9, fric=31,coh=4e6,ten=2e6 range mat=3change jmat(节理)=1 range reg 0,0 0,3 180,3 180,0 ;在这范围1内的块体间相互接触的接触⾯⼒学参数change jmat=2 range reg 0,3 0,4.5 180,4.5 180,3change jmat=3 range reg 0,4.5 0,10.2 180,10.2 180,4.5prop jmat=1 jkn=7e9, jks=2e9, jcoh=0.1e6 , jfric=0,jten=0.1e6 ;在这范围1内的块体间相互接触的接触⾯⼒学参数 prop jmat=2 jkn=6e9, jks=2e9, jcoh=1e6 , jfric=20,jten=1e6prop jmat=3 jkn=6e9, jks=2e9, jcoh=1e6 , jfric=20,jten=1e6;⼯作⾯埋深550⽶set gravity 0,-10 ;地下岩层主要受重⼒,还有构造应⼒(⽔平应⼒=垂直应⼒*侧压系数)bound stress 0,0,-1.375e7 range 0 180 45.1 45.3 ; boundary stress sxx0 ,sxy0, syy0 range x1 x2 y1 y2;(550—煤层到上边界的距离)*2500*(-10)= -1.375e7insitu str -1.5144e7 0 -1.262e7 szz=-1.5144e7 ygrad 3e4 0 2.5e4 zgrad 0 3e4;initi asituate stress 初始设置应⼒ ,侧压系数为1.2;syy=(550+煤层到下边界的距离)*2500*(-10)= -1.262e7（在下边界上施加的⼒）-1.262e7*1.2=-1.5144e7 ,y ⽅向的梯度为=1*2500*10=2.5e4, x 、z ⽅向的梯度为=2.5e4*1.2=3.0e4bound xvel=0 range -0.1 0.1 0 45.2 ;固定左边界(xvel-也就是x ⽅向的速度),见图3bound xvel=0 range 179.9 180.1 0 45.2 ;固定右边界bound yvel=0 range 0 180 -0.1 0.1 ;固定下边界（因为上边界有上覆岩层故不需固定了）solve ;计算save fyuanyan1.save注: insitu str sxx sxy syy szz=-1.6638e8 ygrad sxxy sxyy syyy zgrad szzx szzy(xgrad sxxx sxyx syyy,上⾯没有写,表⽰在x ⽅向没有变化 )(σx=sxx ,σy=syy, σz=szz,τxy=sxy,τxz=sxz, τyz=syz 弹性⼒学;grad 表⽰梯度,即在此⽅向的变化量)(xgrad sxxx sxyx syyy ygrad sxxy sxyy syyy zgrad szzx szzy(sxx=sxx0+sxxx.x+sxxy.y,syy=syy0+syyx.x+syyy.y,szz=szz0+szzx.x+szzy.y)(µ—侧压系数,sxx=syy. µ,szz=syy. µ,sxxy=syyy. µ,szzy=syyy. µ)(τ=c+σ.tg υ,c 值,反映岩⽯剪切时的粘结阻⼒,故称岩⽯的内聚⼒(或粘结⼒) ,υ,值反映岩⽯剪切时摩擦阻⼒,的⼤⼩,故称岩⽯的内摩擦⾓,tg υ,相当于摩擦系数f 。
udec注意事项
3.9.5 单元密度的考虑
UDEC 模型的变形块体采用的是常应变单元。如果具有很高梯度的应力或应变， 描述变量分布则需要很多单元。可采用不同单元密度运行同一问题来检查其效果。在 UDEC 中采用常应变单元是因为当模拟塑性流动时，较多低阶单元比较少高阶单元具 有更高的精度（见理论与背景的第 1.2.5 节） 。 期望使单元尽可能的均匀，尤其在兴趣的区域，尽量避免单元的边长比大于 5:1 的长瘦单元。
连续介质分析是完善的。当存在连续分析是否能代表不连续效应时的问题时，应采用 连续和非连续介质计算。 第二个考虑是包含详细地质结构的计算模型的范围。最关键节理结构的详细 描述（见上述第 1 点）通常仅需要包含感兴趣的有限区域，即在几倍隧道半径的围岩 范围内。一般地，考虑详细节理的范围从感兴趣的区域扩展到足够的距离，包围可能 产生破坏的区域。详细的地质结构应在延伸到节理滑移和张开的范围之外。
3K (1 2 ) E 2G (1 ) E
(3.38)
涉及两组常数。然而，只有当我们所研究的参数接近（但没有达到）这种极限情况下， 方程才成立。对于物理试验，我们并不需要涉及可能或不可能行不通的问题。该方程 是两种定义比例系数方式的简化结果。假设我们有一种逐渐较小抵抗变形的力，但仍 保持常体积变化材料。在此情况下， 接近于 0.5。方程 3K (1 2 ) E 仍必须被满足。 有两种可能性（关于代数基础，而不是物理的争论） ：或者 E 保持为有限值（非零）和 K 趋向于任意大值，或 K 保持为有限值和 E 趋向于零。我们排除第一种可能，因为， 对于所有的材料，都存在压缩的极限值（即对于水是 2GPa ，在此，泊松比为 =0.5） 。 这只可能是第二种情况，即使我们假设材料的主要弹性抗力模式是卸载，E 值也戏剧性 变化。我们推论，用（ E ， ）表达材料特性是不合适的。
UDEC学习
New Features in 3DEC 4.1
Released in January 2008
1. 2. Accelerated interactive graphics (OpenGL based plotting) with new graphical structure New Mixed Discretization scheme for tetrahedral elements “Nodal Mixed Discretization” provides more accurate solution of plasticity problems using tetrahedral grids 64-bit version Factor of Safety calculation mode Improvements to user-defined constitutive models and addition of user-defined joint constitutive models Improvements to structural element logic including installation at tunnel intersections Help File containing Command Reference, FISH Reference and Example Applications New PGEN interface Improved FISH language
full and sectioned views of a tower modeled for stability using the dynamic opt ion
Model of aห้องสมุดไป่ตู้ch dam using finite-element blocks
UDEC模拟实例与解析
UDEC 实例翻译与命令解析翻译：珠穆朗玛1 地震诱发地层坍塌 Seismic-Induced Groundfall1.1 问题描述本例展示使用 UDEC 模拟分析地震诱发地层坍塌的一类的问题,模型见图 1.1,该模型基 于加拿大安大略省萨德伯里市鹰桥公司弗雷则矿 34-1-554 切割断面的一个剖面图的结构和 尺寸. 用二维平面应变模型代表垂直于超采轴向方向的平面效应,超采面高 5m,宽 10m.假定两个连续节理交叉平面分析:一个角度为 45 度,另一个为-9 度,两者节理间距均为 5m,为了演示的目的,一个近似垂直的“虚拟节理”也被添加到块体内开挖面顶部以增强不稳 定性。
围岩参数来自试验室平均测试数值，假定岩石块体参数如下：假定块体仅具有弹性行为，节理假定符合库伦滑动准则，选择典型的教课书数值作为节 理参数，如下：初始应力状态按各向同性估计为24Mpa（假定垂直荷载由覆盖深度大约800m 的岩层产生）。
1.2 UDEC 分析UDEC 模拟顺序分三个阶段,首先,模型在初始应力状态下进行无超采固结.其次,进行开挖并且模型循环至平衡状态.本阶段超采面周围的应力分布见图1.2.超采正上方和下方的块体滑动后稳定.在第三阶段.估计了两个不同的峰值速度的地震事件.对所有地震模拟,在问题域的外周边界引入粘滞边界用以消除波的反射.从而模拟有限的岩体,地震事件用施加到模型顶部y 方向的正弦应力波表现.应力波被叠加到已存在的初始地应力上.在第一个模拟中,施加1.25Mpa 的峰值应力,应当注意的是,由于粘滞边界条件实际是在模型顶部, 施加的有效影响应力应该是1.25 MPa/2, or 0.625 MPa.0.02 秒后的开挖面拱顶的应力分布见图1.3,两点的位移被监测,1 点位于开挖面的左角,点2 位于拱顶块体的右角, 图1.4 的位移时间曲线显示两点本质上是弹性反应.本例关心的问题是在模型顶部施加的速度和计算速度的对比,下面的公式可以用以估计施加的波速.使用这个方程,施加的最大波速大概是0.04m/sec,图1.5 显示的峰值波速小于0.06m/sec. 估计的波速和监测波速的不同在于使用的围岩模量.而是没有考虑节理变形的相等变形模量.在第二个案例中,施加应力波峰值12.5 Mpa(有效应力6.25Mpa).0.02 秒后的开挖拱顶应力分布见图1.6.该图显示出拱顶岩体不受力,表面该块体已经松散并正在下落.对于关心的问题,后来三个时间的几何体和应力分布见图1.8 至图1.10.在问题的顶部预测的波速(从上面的方程)是0.4m/sec.从模型中计算的波速见图1.11,再次,由于使用的是原岩弹性模量而不是岩体的变形模量导致预测和监测的波速之间的差异.1.3 节包含了该模型的数据列表,该列表包含了一个FISH 函数(show)被用来创建坍塌的动画文件,每隔0.02 秒俘获一个显示的图片.通过改变FISH 参数time_int 可以改变动画帧的间隔.视图的总数也可以通过改变snap_shot 的数值进行改变.为了显示80 帧的显示图片而创建的该电影文件需要大概13MB 的硬盘空间.1.3 数据文件列表Example 1.1 SEISMIC.DATtitleSEISMIC INDUCED ROOF COLLAPSE 地震诱发拱顶坍塌;round 0.01; define original boundary of modeled region 定义模型区域的原始边界block -25,-20 -25,20 25,20 25,-20; generate joint pattern over entire original region 在整个原始区域生成节理形态jregion id 1 -25,-25 -25,25 25,25 25,-25jset 45,0 200,0 0,0 5.0,0 (0,0) range jreg 1jset -9,0 200,0 0,0 5.0,0 (0,0) range jreg 1; put in joints needed for the later excavation 为了后面开挖而设置的节理crack -5.01,-2.51 5.01,-2.51crack -5.01, 2.51 5.01, 2.51crack -5,-2.5 -5,2.5crack 5,-2.5 5,2.5crack 2.25,2.5 1.93,5.0; generate fdef zones and assign joint properties (mat=1 & jmat=1;default) 生成单元和设置节理参数generate edge 9.0 range -30,30 -30,30prop mat=1 d=0.00300 k=39060 g=31780prop jmat=1 jkn=20000 jks=20000prop jmat=1 jf=30.0; apply boundary conditions and initial conditions to 在地应力下施加边界条件和初始条件; consolidate model under field stressesbound stress=-24.0, 0.0, -24.0 ygrad=-0.3 0 -0.3insitu stress=-24.0, 0.0, -24.0 ygrad=-0.3 0 -0.3bound yvel 0.0 range -26,26 -21,-19grav 0.0 -10.0; track the x-displacement, and y-displacement over time 追踪位移hist solvehist xdis=0,7 ydis=0,7 type 1solve rat 1e-5; save consolidated statesave seismic1.sav; make excavationdelete range -5,5 -2.5,2.5solve rat 1e-5; save excavated statesave seismic2.sav;rest seismic2.sav; apply seismic load from top (peak velocity=0.04 m/sec);; set up nonreflecting boundarybound mat=1bound xvisc range -26 -23 -21 21bound xvisc range 23 26 -21 21bound xvisc yvisc range -26 26 -21 -19bound xvisc yvisc range -26 26 19 21; apply sinusoidal stress wavebound stress 0 0 -1.25 yhist=cos(100.0,0.0195) range -26 26 19 21 ;reset time hist disp rothist ydis (-4.48,2.57)hist ydis (0,2.57) yvel (0,2.57) yvel (4,2.57) yvel(-4.48,2.57)hist yvel (0,20) yvel (25,10) yvel (25,-10) yvel (0,-20)hist yvel (-25,-10) yvel (-25,10)hist sxx (25,10) sxx (25,-10) sxx (-25,-10) sxx (-25,10)hist syy (0,20);damp 0.1 1.0 mass; 0.02 sec.cyc time 0.02save seismic3.sav;rest seismic2.sav; apply seismic load from top (peak velocity=0.4 m/sec); set up nonreflecting boundarybound mat=1bound xvisc range -26 -23 -21 21bound xvisc range 23 26 -21 21bound xvisc yvisc range -26 26 -21 -19bound xvisc yvisc range -26 26 19 21; apply sinusoidal stress wavebound stress 0 0 -12.5 yhist=cos(100.0,0.0195) range -26 26 19 21 reset time hist disphist ydis (-4.48,2.57)hist ydis (0,2.57) yvel (0,2.57) yvel (4,2.57) yvel(-4.48,2.57)hist yvel (0,20) yvel (25,10) yvel (25,-10) yvel (0,-20)hist yvel (-25,-10) yvel (-25,10)hist sxx (25,10) sxx (25,-10) sxx (-25,-10) sxx (-25,10)hist syy (0,20);damp 0.1 1.0 masssave seismov.sav;; 0.02 sec.cyc time 0.02save seismic4.sav; 0.25 sec.cyc time 0.23save seismic5.sav; 0.50 sec.cyc time 0.25save seismic6.sav; 0.75 seccyc time 0.25save seismic7.sav;rest seismov.sav; make a movie of the groundfall;wind -12 12 -12 12set ovtol 0.05plot block vel max 2.0 blue stress max 50movie onmovie file = seismic.dcxmovie step 1000step 400003 隧道支护荷载Tunnel Support Loading3.1 问题陈述本例模拟展示了UDEC 在检查衬砌隧道方面的应用,着重强调了荷载在混凝土衬砌中的发展,本例也解释了模拟连续建造操作中独立阶段的模拟程序.隧道系统的理想几何体见图3.1.系统包含在海床下大约70m(中线)深度,中线间距12m 的两个隧道, 初始水位在隧道中线上方110m 处.服务隧道直径5.24m,衬砌厚度37cm.主隧道直径8.22m,衬砌厚度46cm.服务隧道先于主隧道开挖和衬砌.随后设置主隧道衬砌,水位上升增加到100m.施工顺序是:(1)开挖服务隧道excavation of the service tunnel;(2)衬砌服务隧道lining of the service tunnel; (3)开挖主隧道excavation of the main tunnel; (4)衬砌主隧道lining of the main tunnel; and (5)升高水位raising of the water level.分析的目的是评价每个施工阶段服务隧道和主隧道支护状况.本例的材料参数见下:岩体——开挖隧道的围岩参数为：弹性模量elastic modulus 0.89 GPa泊松比Poisson’s ratio 0.35单轴抗压强度uniaxial compressive strength 3.5 MPa粘聚力cohesion 1 MPa密度density 1340 kg/m3混凝土衬砌——弹性模量为24 GPa ，泊松比为0.19. 假定衬砌为线弹性材料。
UDEC学习整理资料
1、角点必须按顺时针方向排列；2、Crack命令用于产生块体中单一直线特征的裂缝。
裂缝由端点坐标（x1，y1）和（ x2， y2）所确定。
3、Jset 命令则是自动节理组生成器。
根据所给定的特征参数（即倾角、迹长、岩桥长度、间距和空间位置）产生一组裂缝。
4、round d---d 是圆角距离，建议在block 命令前指定圆角长度。
5、DELETE命令，能从模型中删除一个块体。
例如，为了删除槽口块体，delete range 4.5，5.5 8，10。
6、GEN 命令激活三角形网格有限单元自动生成器。
命令 GEN edge v 将作用于任意形状的块体。
其v 值定义三角形单元的最大边长，即v 值越小，块体中的单元越小。
应当注意的是：具有高的边长比值的块体并不能产生单元，其极限的比重近似为 1:10。
7、采用命令 GEN quad v，指定模型为塑性材料模型的单元。
该类型的单元提供了对于塑性问题的精确解。
然而， GEN quad 命令可能对某些形状的块体不起作用。
在此情况下，应当采用 GEN edge8、Change 命令改变块体为指定的变形块体。
Cons=0 意味着模型块体材料被移出或开挖。
Cons=1 改变块体为各向同性弹性特性；而 Cons=3 则改变块体为摩尔－库仑模型，考虑塑性特性。
缺省值为所有变形体则自动改变为Cons=1。
P219、change jcons=2，所以不连续结构面的缺省模型是Jcons=2。
10、可用以下命令检查材料号Plot block mat11、边界条件命令总结命令效果Stress施加总应力到刚体或变形体块体的边界上Xload施加刚体或变形体边界的x方向的荷载BOUNDARY Yload施加刚体或变形体边界的y方向的荷载Xvel施加变形体边界的x方向的速度（位移）Yvel施加变形体边界的y方向的速度（位移）FIX固定刚体边界的速度（位移）FREE释放刚体的速度（位移）Xload施加 x方向的荷载到刚体的边界LOADYload施加 y方向的荷载到刚体的边界12、INSITU 命令用来初始化应力。
udec版本[整理版]
![udec版本[整理版]](https://img.taocdn.com/s3/m/36d40ea831b765ce04081493.png)
USTB
UDEC软件应用
建模
变化分析
迭代计算
结果输出
UDEC软件应用－－建模
USTB
设计模型尺寸
划分块体和节理
定义本构模型
建模
给出材料力学参数
确定模型边界条件
迭代使之平衡
USTB
UDEC软件应用－－建模
块体生成： Block x1,y1 x2,y2 x3,y3 x4,y4
UDEC软件应用－－建模
USTB
边界条件：Boundary stress sxx sxy syy range x1,x2 y1,y2 Boundary ygrad sxxy sxyy syyy
UDEC软件应用－－边界条件命令总结
USTB
命令
BOUNDARY
Stress
Xload Yload
Xvel Yvel
FIX
FREE
LOAD
Xload Yload
效果
施加总应力到刚体或变形体块体的边界 上 施加刚体或变形体边界的x方向的荷载
施加刚体或变形体边界的y方向的荷载 施加变形体边界的x方向的速度（位移）
施加变形体边界的y方向的速度（位移）
固定刚体边界的速度（位移）
释放刚体的速度（位移） 施加x方向的荷载到刚体的边界 施加y方向的荷载到刚体的边界
udec版本[整理版]
UDEC简介－应用范畴
USTB ◆ 便于处理所有非线性变形和破坏都集中在节理面上为特征的岩体破坏问题。 ◆ 求解连续介质及连续介质向非连续介质转化的力学问题。 ◆ 广泛应用于模拟边坡，滑坡和节理岩体地下水渗流。 ◆ 土壤开裂，分离等非连续现象。 ◆ 围岩洞室破坏诱发断裂，滑移。 ◆ 特别适用于模拟计算节理岩体破坏与大变形，地下结构和深部高辐射废料的储存 研究。
【VIP专享】UDEC经典学习总结
set log off
3.plot block stress ;显示块体的应力
plot block dis ;显示块体的位移(有 x 和 y 方向)
plot block pl ;显示块体的塑性区(plastic)
save xx.save ;保存计算结果
res xx.save ;调用
4.液压支架的命令
table 1 0 42.e6 0.05 5.0e6 0.1 6.0e6;表 1 表示的是液压支架的(P-DS)特性曲线
5 巷道施工中断面加梯子梁的命令(见图 1)
stru gen xc yc np 100 fa a thetra b mat=16 thick=0.2 ;100 指分 100 段,a b 指角度
supp xc yc wid l seg n mat j;(xc,yc)是指中心点坐标,l 指支架的宽度,n 指分段数,mat j 指支柱材料性质为 j
prop mat j sup_kn -1
;sup_kn 指支架的刚度,-1 与表 1 相对应
del range x1 x2 y1 y2
;挖掘范围(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)的块体,由支架支撑
prop mat=16 st_ymod=13.5e9 st_prat=0.14 st_den=7800
6.培养学生观察、思考、对比及分析综合的能力。过程与方法1.通过观察蚯蚓教的学实难验点，线培形养动观物察和能环力节和动实物验的能主力要；特2.征通。过教对学观方察法到与的教现学象手分段析观与察讨法论、，实对验线法形、动分物组和讨环论节法动教特学征准的备概多括媒，体继课续件培、养活分蚯析蚓、、归硬纳纸、板综、合平的面思玻维璃能、力镊。子情、感烧态杯度、价水值教观1.和通过学理解的蛔1.虫过1观适、察于3程.观六 阅蛔寄内列察、 读虫生出蚯让 材标容生3根.常蚓学 料本教活 2据 了见身生 ：，师的、问 解的体巩总 鸟看活形用 作题 蛔线的固结 类雌动态手 业自 虫形练本 与雄学、三4摸 ：学 对动状习节 人蛔生结、、一 收后 人物和同课 类虫活构请一蚯摸 集回 体并颜步重 关的动、学、蚓蚯 鸟答 的归色学点 系形教生生让在蚓 类问 危纳。习并 从状学理列学平的题 害线蚯四线归 人、意特出四面体 生以形蚓、形纳 类大图点常、五观玻表 存及动的鸟请动本 文小引以见3引、察璃， 现预物身类学物节 明有言及的、导巩蚯上是 状防的体之生和课 历什根蚯环怎学固蚓和干 ，感主是所列环学 史么据蚓节二样生练引牛燥 鸟染要否以举节到 揭不上适动、区回习导皮还 类的特分分蚯动的 晓同节于物让分答。学纸是 减方征节布蚓物一 起，课穴并?学蚯课生上湿 少法。广的教些 ，体所居归4在生蚓前回运润 的；.泛益鸟色学生纳4靠物完的问答动的 原.蛔，处目类和习活环近在成前题蚯的？ 因了虫以。标生 就体的节身其实端并蚓快触 及解寄上知物 同表内特动体结验和总利的慢摸 我蚯生适识学 人有容点物前构并后结用生一蚯 国蚓在于与的 类什，的端中思端线问活样蚓的人? 飞技基 有么引进主的的考形题环吗十生体能本 着特出行要几变以动1境？节 大活，的.方 密征本“特节化下物让并为近 珍习会形理法 切。课生征有以问的2小学什腹 稀性引态解。 的.课物。什游题主结生出么面 鸟和起结蛔关观题体么戏：要利明蚯？处 类适哪构虫系察：的特的特用确蚓， 等于些特适。蛔章形殊形征板，这是 资穴疾点于可虫我态结式。书生种光 料居病寄是的们结构，五小物典滑 ，生？重生鸟内学构5，学、结的型还 以活.要生类部习与其习巩鸟结是 爱的如原活生结了功颜消固类构线粗 鸟形何因的存构腔能色化练适特形糙 护态预之结的，肠相是系习于点动？ 鸟、防一构现你动适否统。飞都物为结蛔。和状认物应与的行是。主构虫课生却为和”其结的与题、病本理不蛔扁的他构8特环以生？8特乐虫形观部特征境小理三页观的动点位点梳相组等、这；，哪物教相，2理适为方引.鸟些，育同师知应单面导鸟掌类结了；？生2识位的学你握日构解.互。办特生认线益特了通动手征观识形减点它过，抄；察吗动少是们理5生报.蛔？物，与的解参一了虫它和有寄主蛔与份解结们环些生要虫其。蚯构都节已生特对中爱蚓。会动经活征人培鸟与飞物灭相。类养护人吗的绝适这造兴鸟类？主或应节成趣的为要濒课的情关什特临？就危感系么征灭来害教；？；绝学，育，习使。我比学们它生可们理以更解做高养些等成什的良么两好。类卫动生物习。惯根的据重学要生意回义答；3的.情通况过，了给解出蚯课蚓课与题人。类回的答关：系线，形进动行物生和命环科节学动价环值节观动的物教一育、。根教据学蛔重虫点1病.引蛔出虫蛔适虫于这寄种生典生型活的线结形构动和物生。理二特、点设；2置.问蚯题蚓让的学生生活思习考性预和习适。于穴居生活的形态、结构、生理等方面的特征；3.线形动物和环节动物的主要特征。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1、角点必须按顺时针方向排列；
2、Crack 命令用于产生块体中单一直线特征的裂缝。
裂缝由端点坐标（x1，y1）和（x2，y2）所确定。
3、Jset 命令则是自动节理组生成器。
根据所给定的特征参数（即倾角、迹长、岩桥长度、间距和空间位置）产生一组裂缝。
4、round d---d是圆角距离，建议在block命令前指定圆角长度。
5、DELETE 命令，能从模型中删除一个块体。
例如，为了删除槽口块体，delete range 4.5，5.5 8，10。
6、GEN命令激活三角形网格有限单元自动生成器。
命令GEN edge v 将作用于任意形状的块体。
其v值定义三角形单元的最大边长，即v值越小，块体中的单元越小。
应当注意的是：具有高的边长比值的块体并不能产生单元，其极限的比重近似为1:10。
7、采用命令GEN quad v，指定模型为塑性材料模型的单元。
该类型的单元提供了对于塑性问题的精确解。
然而，GEN quad 命令可能对某些形状的块体不起作用。
在此情况下，应当采用GEN edge
8、Change 命令改变块体为指定的变形块体。
Cons=0意味着模型块体材料被移出或开挖。
Cons=1 改变块体为各向同性弹性特性；而Cons=3则改变块体为摩尔－库仑模型，考虑塑性特性。
缺省值为所有变形体则自动改变为Cons=1。
P21
9、cha nge jcons=2，所以不连续结构面的缺省模型是Jcons=2。
10、可用以下命令检查材料号Plot block mat
12、INSITU命令用来初始化应力。
采用该命令，可以赋值初始应力。
13、hist xvel 5, 5 hist ydisp 0, 11 第一个是记录位移坐标（x=5，y=5）附近结点x方向的速度，而第二个是记录接近坐标（x=0，y=11）位置处y方向的位移。
14、set grav 0.0 , -9.81第一个是x方向的加速度，第二个值为y方向的加速度为9.81m/sec2（向下作用）。
15、UDEC计算可在任何时间通过按<Esc>被中断。
更方便的是使用STEP命令进行高次数的计算和周期的中断和再次分析，以确保达到平衡状态。
16、用DELETE命令或CHANGE cons=0命令模拟材料开挖。
用BOUNDARY xload，yload 或stress 命令施加荷载和应力。
通过采用BOUNDARY xvel 或yvel 命令固定边界角点。
通过BOUNDARY xfree 和yfree 命令移去边界约束。
用CHANGE命令改变变形块体和不连续面的材料模型。
而用PROPERTY命令可改变材料性质参数。
17、Z放大，ctrl+z恢复原样。
Pl ot blo ck stress disp lacement（块体应力位移）；pl ot hist 1（最大不平衡力）；
Pl+con+szz（z方向应力云图），Pl+ x disp（位移云图）。
18、在一个阶段的结尾（即初始平衡），采用Save file.sav命令，可以保存模型状态。
这个文件可以采用rest file.sav命令进行恢复。
式中，file.sav 是一个用户定义的文件名。
扩展名.sav定义这个文件是一个保存文件。
19、
20、SET edge命令运行用户人为定义的最小块体棱长。
基于这个命令，用户为精确求解而避免块体具有小的棱长，或反之具有很大的边长比，能够设定一个较小的圆角长度。
例如，如果设置命令SET edge 0.4 和ROUND 0.1，则不会产生小于0.4的块体棱长和块体的圆角长度将为0.1。
这些命令必须在BLOCK命令之前给出。
21、模型边界必须具有足够的远，以致模型对边界不产生影响。
一般地，对于单一地下开挖工程，边界离开挖边界的距离应当大于开挖跨度的5倍左右。
然而，合适的距离取决于分析的目的。
如果分析目的主要用于考虑破坏，然而，模型边界可以靠近一些。
如果关注的是位移（变形），则距离边界的距离需要增加。
22、一个节理组可以通过8个生成参数表述：4个几何参数的均值和4个随机参数的均方差。
由JSET命令给出的参数。
P38
23、JREGION定义了一个凸多边形区域来限制节理组的生成范围。
该命令参数如下：JREGION id n x1 y1 x2 y2 x3 y3 x4 y4 <delete>
24、JSET生成器也可能通过设定很大的t m和s m 用来产生贯穿UDEC块体的一条节理。
请记住，JSET能够产生不连续节理。
不连续节理完全处在块体里边不能看见。
当一个变形块体开始生成单元或当一个刚体模型开始运行时，所有的节理段将被删除。
25、用下列命令实施边界形状的产生：CRACK ，JSET ，TUNNEL ，ARC。
每一个命令切割块体成一条或多条裂缝，集合组成所要求的形状。
CRACK和JSET命令产生节理的直线段。
TUNNEL和ARC命令分别形成圆形线段和弧段。
26、隧道节理性质由下列命令进行改变：change jmat=2 range ann 0，0 1.99，2.01 在圆心为（0,0）、半径r1=1.99和r2=2.01之间的环状节理将具有节理性质号为2。
该变化将由以下命令进行查看：plot mat joint
27、泊松效应将定义为：当在垂直方向施加荷载的条件，不允许在水平方向产生应变（位移）时水平应力与垂直应力之比。
28、UDEC模型的缺省边界是无约束的自由边界。
力或应力可以通过BOUNDARY命令施加到任意整个边界或部分边界上。
用stress关键词可以指定平面应力张量（σyy、σxy、σxx）的每一个单独应力分量。
例如，命令
Boundary stress 0，-1e6，-2e6 range 0，10 -1，1
将施加σyy=0、σxy =-106、σxx =-2×106到位于坐标窗口0< x <10，-1<y<1范围内模型边界上。
用户应检查窗体周围的所有边界角点所指定的边界条件，这可用下面的命令实现：Print boundary
29、UDEC中按照一般符合惯例，压应力为负号。
而且，UDEC实际上施加应力分量作为力或者是产生作用到给定的边界平面上的应力张量的摩擦力。
摩擦力被分成两个分量：永久的和瞬时的。
30、BOUNDARY命令可以增加关键词xgrad和ygrad，此关键词允许应力或力在指定的边界上按线性变化。
31、瞬时荷载可用history关键词，为动力分析施加边界条件。
对于静力分析，在UDEC模拟过程中，改变边界的应力值也是必要的。
例如，作用到基础的荷载可能发生改变。
为了考虑应力或荷载突然变化的影响，需要给出新的BOUNDARY命令，对原来边界角点施加改变后的应力或荷载。
32、在这种情况下，新的边界应力或荷载值将被叠加到原有的值上。
如果应力被移动，当前的值应当以相反的符号给出。
如果一瞬时荷载被改变（即用history 关键词指定荷载），用相同的history 所给出新的荷载被叠加到存在的荷载上。
然而，用不同的history所给出的一个新的瞬时荷载，将取代老的荷载。
33、边界应力和荷载可以通过命令PRINT bound 和PLOT bound xcond或PLOT bound ycond进行验证。
PRINT bound 命令列出每一角点所指定的边界角点地址和对应的当前的值，PLOT bound xcond 或ycond 命令通过符号表示施加在x-或y-方向边界条件的类型。
需要施加边界条件的边界角点地址能够通过以下命令输出显示。
Print bound state
34、应力边界影响所有的自由度。
所以速度边界条件必须在应力边界条件之后给出。
如果应力边界在速度边界之后施加，给出速度的响应将会消失。
35、通过在指定的速度关键词之后增加关键词gvel，速度也允许在指定的边界范围按线性变化。
Gvel关键词后有六个参数，描述了速度分量在x或y方向的变化：gvel vx0 vy0 vxx vxy vyx vyy P50
36、。