下载文档原格式
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、cha nge jcons=2，所以不连续结构面的缺省模型是Jcons=2。
10、可用以下命令检查材料号Plot block mat12、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（向下作用）。
通用离散元用户指导（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 ⼊门;new 是刷新udec窗⼝，从新调⽤⼀个程序;title 与heading代表标题，后⾯紧跟标题的名称。
如：titlehang dao mo ni;round 指块体与块体之间的圆⾓半径，默认值是0.5，其值要求⼩于模型中最⼩块体的最短那条边长的⼆分之⼀。
如：round 0.05set ovtol=0.5;此命令是指层与层之间的嵌⼊厚度block x1,y1 x2,y2 x3,y3 x4,y4;建⽴模型框架，crack x1,y1 x2,y2;两点划⼀线jregion id n x1,y1 x2,y2 x3,y3 x4,y4 deletejset 90,0 4,0 4,0 6,0 0,-50 range jreg 3;jset 倾⾓,0 线段长,0 线段与线段轴向间隔长,0 垂向间距,0 xm,ym range jregion n;其中xm,ym为起始点坐标,n为设置的区域标号gen quad 10 range xl xu yl yu;在指定的区域⽣成⼀定宽度的单元（xu为x⽅向的取值）zone model mo range xl xu yl yu;使指定的区域材料采⽤摩尔--库仑本构关系计算（即弹塑性）change jcons=2 range xl xu yl yu;使指定的区域节理遵循摩尔--库仑准则计算（即弹塑性）change mat=1 range xl xu yl yuchange mat=2 range xl xu yl yuchange mat=3 range xl xu yl yu;指定各岩层的材料标号change jmat=1 range xl xu yl yuchange jmat=2 range xl xu yl yuchange jmat=3 range xl xu yl yu;指定各岩层的节理标号prop mat=1 dens=2000prop mat=2 dens=2650prop mat=3 dens=2700;指定各材料的密度，⽐如1号材料dens=2000，即1⽴⽅⽶重2吨zone k=0.15e9,g=0.1e9,fric=10.00,coh=0.19e6,ten=0.09e6 range mat=1zone k=2.8e9,g=2.2e9,fric=30.00,coh=1.5e6,ten=0.4e6 range mat=2zone k=6.9e9,g=6.6e9,fric=38.62,coh=5.63e6,ten=3.20e6 range mat=3;k为材料的法向刚度，g为材料的切向刚度，friction为材料的内摩擦⾓，;cohesion为材料的内聚⼒，tension为材料的抗拉强度prop jmat=1 jkn=0.2e8,jks=0.1e7,jcoh=0,jfric=4,jten=0prop jmat=2 jkn=8e8,jks=5e7,jcoh=0.1e6,jfric=8,jten=0prop jmat=3 jkn=20e8,jks=16e7,jcoh=0.4e6,jfric=15,jten=0;jkn为节理的法向刚度，jks为节理的切向刚度，jfriction为节理的内摩擦⾓，;jcohesion为节理的内聚⼒，jtension为节理的抗拉强度set gravity 0,-9.81;设置重⼒加速度，x⽅向为0，y⽅向为-9.8bound xvel=0 range -0.1 3.00 -60.1 20.1bound xvel=0 range 97 100.1 -60.1 20.1bound yvel=0 range 0.1 100.1 -60.1 -58;采⽤位移法固定边界solve\step 5000\cycle 5000;执⾏计算save pingheng.sav;保存⽂件，⽂件的后缀为.sav，⽂件名可以⾃⼰命名。
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、cha nge jcons=2，所以不连续结构面的缺省模型是Jcons=2。
10、可用以下命令检查材料号Plot block mat12、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（向下作用）。
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、
SAVE
保存或恢复当前状态
HISTORY RESTORE 监测模型响应
PLOT
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、。
