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、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（向下作用）。

dos命令大全黑客必知的DOS命令集合_DOS/BAT net use ipipc$ "" /user:"”建立IPC空链接net use ipipc$ "密码”/user:"用户名"建立IPC非空链接net use h: ipc$ "密码" /user:"用户名"直接登陆后映射对方C :到本地为H:net use h: ipc$登陆后映射对方C :到本地为H:net use ipipc$ /del删除IPC链接net use h: /del删除映射对方到本地的为H:的映射net user用户名密码/add 建立用户net user guest /active:yes激活guest用户net user查看有哪些用户net user帐户名查看帐户的属性net locaLGroup administrators用户名/add把”用户"添加到管理员中使其具有管理员权限注意: administrator后加s用复数net start查看开启了哪些服务net start服务名开启服务 ; (如:net start telnet，net start schedule)net stop服务名停止某服务net time目标ip查看对方时间net time目标ip /set设置本地计算机时间与目标IP”主机的时间同步加上参数/yes 可取消确认信息net view查看本地局域网内开启了哪些共享net view ip查看对方局域网内开启了哪些共享net config显示系统网络设置net logoff断开连接的共享net pause服务名暂停某服务net send ip "文本信息"向对方发信息net ver局域网内正在使用的网络连接类型和信息net share查看本地开启的共享net share ipc$开启ipc$共享net share ipc$ /del删除ipc$共享net share c$ /del删除C :共享net share c$ /del删除C:共享net user guest 12345用guest用户登陆后用将密码改为12345net password密码更改系统登陆密码netstat -a查看开启了哪些端口，常用netstat -annetstat -n查看端口的网络连接情况,常用netstat -annetstat -V查看正在进行的工作netstat -p协议名例: netstat -p tcq/ip查看某协议使用情况(查看tcp/ip协议使用情况)netstat -s查看正在使用的所有协议使用情况nBTstat -A ip对方136到139其中一一个端口开了的话,就可查看对方最近登陆的用户名( 03前的为用户名) -注意:参数-A要大写trAcert -参数ip(或计算机名)跟踪路由(数据包) , 参数: "-w数字”用于设置超时间隔。

UDEC60常用命令集锦rest falll.savdelete range -2,2 -2,2reset disp reset hist hist un bal hist ydis 0,2 step 2000plot hold block stress disp ；位移归0edgeExample 3.2 EJfta rtrundit^ latglh on crack getteraiionro D .2；ro D . 1bi 0.0 o f io i0f in lo.o cr 0,3.09.7f10]f the rounding length ia reduced (o O. I. then the crack will be located where sfyecified.I lie SET edgeThisconimana, the u^er 匚日n set a small rnunding length for solution accuracy bur avoid blocks with small edg亡lengths and, consequently, adversel} high aspect ratios (see below). For exajnpk. ifcommand allows the user to Ueiine a minimum blxKvd*亡kngch nianually. WithLhe SET edge 0.4 and ROUND 0.1 arc $pe(;ifk也then block 、tnal吐山an G4 VrilL ilol be cruaicMl. and the rounding length for blocks ill be OJ. These coinnunds be gi\ en beLor^tlie BLOCK coinmaiid.Joints jregionjset angle 135 trace 1 gap 2 spacing 1.732 origin 0.259,0.966 range jreg 1 ;jregion 2 jointsjset angle jset angle 0 trace 10 trace 1gap 2 spacing 1?732 origin Q a0 range jreg 2gap 2 spacing 1.732 origin 1?5.0.866 range jreg 2jset angle 60 trace 1 gap 2 spacing 1?732 origin 1.0 range jreg 2jset angle 60 trace 1 gap 2 spacing 1?732 origin. 1.1.732 range jreg 2Joints that are joined are assigned high values for normal and shear stiffness. It is tempting to give VCQ high values for stiffnesses to prevent movement along a fictitious joint. 1 lowcvcr, the timestep calculation in UDEC is based upon stiffnesses; the response (and solution convergence) will be very slow if very high siifTncsses are spexified. Tlic lowest MifTiiess consisieni with small joint deformaiion is used when joints are joined. The rule-of-thumb is that joint stifTnesses t I n and k s. should be set to a factor times the equivalent stiffness of the stiffest neighboring zone. The equation to calculate fictitious joint normal stiffness related to the equivalent stiffness (expressed in stress-per-distance units) of a zone in the normal direction is of the form k n = f actor x max(3.1) Figure 3.14 Zone dimension used ui stiffness calculationSystems of Units2.S Systems of L'nilsUDEC accepts any con listen f of engineering units, Examples of consis：enr sets of units forh阳it parameters are shown in Table 2.5. The user should be very careful when converting from one stem of uniK to unoiher. An excellent reference on the subject or units and conversion between the Imperiikl and SI s\ stems can be iimnd in the Juunuti。

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 21reset 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. —————————————————————————————————————UDEC 实例翻译与命令解析中铁隧道集团科研所——珠穆浪玛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.—————————————————————————————————————UDEC 实例翻译与命令解析中铁隧道集团科研所——珠穆浪玛施工顺序是:(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. 假定衬砌为线弹性材料。

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与横向层理有关。

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所示的数值计算模型。

;newtitle moni ;定义名称moni块体都有“圆角”，其目的在于避免块体悬挂在有棱角的节点上。

由于块体悬挂引起应力集中。

然而，圆角值存在与模型有关的上限值。

对于变形块体，最大圆角长度应当不超过块体平均棱长的1％。

圆角长度可以如下命令加以改变：round dround=0.025；d 是圆角距离（缺省值是d=0.5）。

模型中的所有圆角长度都是相同的。

For example, if the commands SET edge 0.4and ROUND 0.1 are specified, then block edge lengths smaller than 0.4 will not be created, and the rounding length for blocks will be 0.1. These commands must be given before the BLOCK command.set ovtol=1;块体与块体之间相互嵌入量最大值为1米bl 0,0 0,100 300,100 300,0 ;定义范围，四点坐标，顺时针方向cr 0,15 300,15;Crack 命令用于产生块体中单一直线特征的裂缝，裂缝由端点坐标（x1，y1）和（x2，y2）所确定。

jreg id 1 0,0 0,15 300,15 300,0 delete ;定义节理，命名1，定义范围，四点坐标，delete为常规语言jset 90,0 5,0 5,0 5,0 0,0 rang jreg 1;Jset 命令则是自动节理组生成器。

根据所给定的特征参数（即倾角、迹长、岩桥长度、间距和空间位置）产生一组裂缝。

定义节理角度90°；节理的长度；节理的距离，即纵向间隔；横向距离，即横向隔5一个；起始点坐标为0,0jset 90,0 5,0 5,0 5,0 2.5,5 rang jreg 1 ;定义节理的另外一项gen quad 10 6 range 0 300 0 100 ;定义块体最大变形,若没有此语句,刚所有块体均为刚性块体；automatic generation of diagonally opposed triangular zones to improve plastic flow calculation.(对角三角形区域的自动生成改善塑性流动计算) Parameters xw and yw are zone widths in the x-and y-directions zone model mo range 0 300 0 100 ;定义摩尔库伦模型的范围，X的范围，Y的范围change jcons=5 range 0 300 0 100 ;定义节理的某个属性change mat=1 range reg 0,0 0,15 300,15 300,0 ;定义物质1的范围，即赋予这个范围为物质1,只赋名，没有定义属性change mat=2 range reg 0,15 0,20 300,20 300,15change jmat=1 range reg 0,0 0,15 300,15 300,0 ;定义节理1的范围set jmatdf 2 ;定义节理属性，使符合摩尔库伦准则save ch.sav ;保存文件cal prop-25.txt ;调出文件res ch.sav ;调出前面保存的文件接着计算;参数设置prop mat=1 dens=2720 ;物质的密度zone k=12.12e9 g=10.26e9 fric=33 coh=5.77e6 ten=3.6e6 range mat=1;定义物质1的体积模量、剪切模量、内摩擦角、粘聚力、抗拉强度prop jmat=1 jkn=11e8 jks=12e7 jfric=12 jcoh=0 jten=0;定义节理的属性法向刚度、切向刚度、内摩擦角、粘聚力、抗拉强度;边界条件boun stress 0 0 0 range 0,300 99.9 100.1 ;定义边界条件，三向应力为0，上边界，范围，纵向波动范围boun stress 0 0 -8.6095397e6 range ;定义边界应力条件，竖直方向的应力为，负号表示方向向下insitu str -7.5e6 0 -15e6 szz -7.5e6 ygrad 1.1018e4 0 2.203607e4 ;定义应力属性，sxx方向即水平应力，sxy方向即剪切应力，syy方向即竖直应力，szz方向的应力，ygrad应力梯度set gravity 0,-10 ;重力加速度bound xvel=0 range -0.1 0.1 0 100 ;X位移边界,位移波动范围，范围bound xvel=0 range 299.9 300.1 0 100 ;X位移边界，同上bound yvel=0 range -0.1 300.1 -0.1,0.1 ;Y位移边界hist solve_ratio type 1 ;设置不平衡力的精度，普适solve ;计算save 25-1.sav ;保存cal 25-2.txt ;调出文件res 25-1.sav ;调出前面保存的文件del 80,110 15,20;永久删除这个范围，即开挖，不可再充填，X的范围，Y的范围model null range ;变性删除这个范围，即之后可以充填物质，X的范围，Y的范围model m range ;充填这个范围，solve ;计算save 25-301.sav ;再保存，周而复始set pline 25,20 300,20 10 ;定义测线，起始点坐标，测线分段, 前两组数字表示坐标最后一个数表示布置了多少测点reset histset log on ;定义hist,保存在一个位置，在程序中输入hist1即可调出hist xdisp 10,20 ;测该点的X方向位移hist sxx 10,20 ;测该点X方向的应力hist syy 10,20 ;测该点Y方向的应力reset dispreset hist ;位移清零，历史记录清零his unbalsolve ratio 1e-5 ;记录不平衡力，计算直到精度达到1e-5ch cons=3 range 77.5,80 0,3ch mat=7 range 77.5,80 0,3 ;应变软化模型，这个是做充填用的cable (73.4,2.91) (71.15,4.99) 10 12 380e-6 11 ;锚杆，起始点坐标，后面四个数值可固定prop mat 12 cb_dens 7800 cb_ycomp 430e6 cb_yield 160e3 cb_ymod1.3e11prop mat 11 cb_kbond 6.3e9 cb_sbond 6e5;定义锚杆的属性，密度、抗拉强度、屈服强度density for block materical 1 is zero cannot cycle----块体密度未设range 范围;density 密度;volume 体积;regedit 调出注册表;pl bound pl 显示塑性区;pl bound dis 显示位移矢量图continue pa以后继续算number 可以不关udec,继续使用udecres xx.save (先调用文件)set pline x1,y1 x2,y2 n (n--观测线分的段数)set log on ；固定的set log 文件名print+pline+n+ syy sxx or xdis (观测线名称任意数字)set log off1) pl bl hist 1,2,... 观测点位移2) pl bl nu(mber) 显示块体标号3) pl bl cab red stru red cab--锚杆stru--梁4) pl bl sxx/syy/szz 显示应力图5) pl bl dis 位移矢量图6) set pl dxf 256 把UDEC 图引入CAD 换行copy 文件名自定.dxf块体力学参数k--体积模量;g--剪切模量;fric--摩擦角;coh--粘结力（内聚力）;ten--抗拉强度;d--体积力接触面力学参数jkn--法向刚度;jks--切向强度;jc--粘结力;jf--摩擦角;jt--抗拉强度1) D--质量密度F--摩擦角B--体积模量; C--内聚力(粘结力) S--剪切模量2) 应力--正号代表张力，负号代表压力应变--正的应变代表伸长，负的应变代表压缩重力--正号的重力物质往下拉，负号的重力将物质往上提；内存赋值udec mudec 14print memhist yvel (20,20) type 1; type 是在屏幕上以指定的间隔显示其值step 一般典型问题要2000～4000次循环。

;newtitle moni ;定义名称moni块体都有“圆角”，其目的在于避免块体悬挂在有棱角的节点上。

由于块体悬挂引起应力集中。

然而，圆角值存在与模型有关的上限值。

对于变形块体，最大圆角长度应当不超过块体平均棱长的1％。

圆角长度可以如下命令加以改变：round dround=0.025；d 是圆角距离（缺省值是d=0.5）。

模型中的所有圆角长度都是相同的。

For example, if the commands SET edge 0.4and ROUND 0.1 are specified, then block edge lengths smaller than 0.4 will not be created, and the rounding length for blocks will be 0.1. These commands must be given before the BLOCK command.set ovtol=1;块体与块体之间相互嵌入量最大值为1米bl 0,0 0,100 300,100 300,0 ;定义范围，四点坐标，顺时针方向cr 0,15 300,15;Crack 命令用于产生块体中单一直线特征的裂缝，裂缝由端点坐标（x1，y1）和（x2，y2）所确定。

jreg id 1 0,0 0,15 300,15 300,0 delete ;定义节理，命名1，定义范围，四点坐标，delete为常规语言jset 90,0 5,0 5,0 5,0 0,0 rang jreg 1;Jset 命令则是自动节理组生成器。

根据所给定的特征参数（即倾角、迹长、岩桥长度、间距和空间位置）产生一组裂缝。

定义节理角度90°；节理的长度；节理的距离，即纵向间隔；横向距离，即横向隔5一个；起始点坐标为0,0jset 90,0 5,0 5,0 5,0 2.5,5 rang jreg 1 ;定义节理的另外一项gen quad 10 6 range 0 300 0 100 ;定义块体最大变形,若没有此语句,刚所有块体均为刚性块体；automatic generation of diagonally opposed triangular zones to improve plastic flow calculation.(对角三角形区域的自动生成改善塑性流动计算) Parameters xw and yw are zone widths in the x-and y-directions zone model mo range 0 300 0 100 ;定义摩尔库伦模型的范围，X的范围，Y的范围change jcons=5 range 0 300 0 100 ;定义节理的某个属性change mat=1 range reg 0,0 0,15 300,15 300,0 ;定义物质1的范围，即赋予这个范围为物质1,只赋名，没有定义属性change mat=2 range reg 0,15 0,20 300,20 300,15change jmat=1 range reg 0,0 0,15 300,15 300,0 ;定义节理1的范围set jmatdf 2 ;定义节理属性，使符合摩尔库伦准则save ch.sav ;保存文件cal prop-25.txt ;调出文件res ch.sav ;调出前面保存的文件接着计算;参数设置prop mat=1 dens=2720 ;物质的密度zone k=12.12e9 g=10.26e9 fric=33 coh=5.77e6 ten=3.6e6 range mat=1;定义物质1的体积模量、剪切模量、内摩擦角、粘聚力、抗拉强度prop jmat=1 jkn=11e8 jks=12e7 jfric=12 jcoh=0 jten=0;定义节理的属性法向刚度、切向刚度、内摩擦角、粘聚力、抗拉强度;边界条件boun stress 0 0 0 range 0,300 99.9 100.1 ;定义边界条件，三向应力为0，上边界，范围，纵向波动范围boun stress 0 0 -8.6095397e6 range ;定义边界应力条件，竖直方向的应力为，负号表示方向向下insitu str -7.5e6 0 -15e6 szz -7.5e6 ygrad 1.1018e4 0 2.203607e4 ;定义应力属性，sxx方向即水平应力，sxy方向即剪切应力，syy方向即竖直应力，szz方向的应力，ygrad应力梯度set gravity 0,-10 ;重力加速度bound xvel=0 range -0.1 0.1 0 100 ;X位移边界,位移波动范围，范围bound xvel=0 range 299.9 300.1 0 100 ;X位移边界，同上bound yvel=0 range -0.1 300.1 -0.1,0.1 ;Y位移边界hist solve_ratio type 1 ;设置不平衡力的精度，普适solve ;计算save 25-1.sav ;保存cal 25-2.txt ;调出文件res 25-1.sav ;调出前面保存的文件del 80,110 15,20;永久删除这个范围，即开挖，不可再充填，X的范围，Y的范围model null range ;变性删除这个范围，即之后可以充填物质，X的范围，Y的范围model m range ;充填这个范围，solve ;计算save 25-301.sav ;再保存，周而复始set pline 25,20 300,20 10 ;定义测线，起始点坐标，测线分段, 前两组数字表示坐标最后一个数表示布置了多少测点reset histset log on ;定义hist,保存在一个位置，在程序中输入hist1即可调出hist xdisp 10,20 ;测该点的X方向位移hist sxx 10,20 ;测该点X方向的应力hist syy 10,20 ;测该点Y方向的应力reset dispreset hist ;位移清零，历史记录清零his unbalsolve ratio 1e-5 ;记录不平衡力，计算直到精度达到1e-5ch cons=3 range 77.5,80 0,3ch mat=7 range 77.5,80 0,3 ;应变软化模型，这个是做充填用的cable (73.4,2.91) (71.15,4.99) 10 12 380e-6 11 ;锚杆，起始点坐标，后面四个数值可固定prop mat 12 cb_dens 7800 cb_ycomp 430e6 cb_yield 160e3 cb_ymod1.3e11prop mat 11 cb_kbond 6.3e9 cb_sbond 6e5;定义锚杆的属性，密度、抗拉强度、屈服强度density for block materical 1 is zero cannot cycle----块体密度未设range 范围;density 密度;volume 体积;regedit 调出注册表;pl bound pl 显示塑性区;pl bound dis 显示位移矢量图continue pa以后继续算number 可以不关udec,继续使用udecres xx.save (先调用文件)set pline x1,y1 x2,y2 n (n--观测线分的段数)set log on ；固定的set log 文件名print+pline+n+ syy sxx or xdis (观测线名称任意数字)set log off1) pl bl hist 1,2,... 观测点位移2) pl bl nu(mber) 显示块体标号3) pl bl cab red stru red cab--锚杆stru--梁4) pl bl sxx/syy/szz 显示应力图5) pl bl dis 位移矢量图6) set pl dxf 256 把UDEC 图引入CAD 换行copy 文件名自定.dxf块体力学参数k--体积模量;g--剪切模量;fric--摩擦角;coh--粘结力（内聚力）;ten--抗拉强度;d--体积力接触面力学参数jkn--法向刚度;jks--切向强度;jc--粘结力;jf--摩擦角;jt--抗拉强度1) D--质量密度F--摩擦角B--体积模量; C--内聚力(粘结力) S--剪切模量2) 应力--正号代表张力，负号代表压力应变--正的应变代表伸长，负的应变代表压缩重力--正号的重力物质往下拉，负号的重力将物质往上提；内存赋值udec mudec 14print memhist yvel (20,20) type 1; type 是在屏幕上以指定的间隔显示其值step 一般典型问题要2000～4000次循环。

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 。

rest falll.savdelete range -2,2 -2,2reset disp reset hist hist un bal hist ydis 0,2 step 2000plot hold block stress disp ；位移归0edgeExample 3.2 EJfta rtrundit^ latglh on crack getteraiionro D .2；ro D . 1bi 0.0 o f io i0f in lo.o cr 0,3.09.7f10]f the rounding length ia reduced (o O. I. then the crack will be located where sfyecified.I lie SET edgeThisconimana, the u^er 匚日n set a small rnunding length for solution accuracy bur avoid blocks with small edg亡lengths and, consequently, adversel} high aspect ratios (see below). For exajnpk. ifcommand allows the user to Ueiine a minimum blxKvd*亡kngch nianually. WithLhe SET edge 0.4 and ROUND 0.1 arc $pe(;ifk也then block 、tnal吐山an G4 VrilL ilol be cruaicMl. and the rounding length for blocks ill be OJ. These coinnunds be gi\ en beLor^tlie BLOCK coinmaiid.Joints jregionjset angle 135 trace 1 gap 2 spacing 1.732 origin 0.259,0.966 range jreg 1 ;jregion 2 jointsjset angle jset angle 0 trace 10 trace 1gap 2 spacing 1・732 origin Q a0 range jreg 2gap 2 spacing 1.732 origin 1・5.0.866 range jreg 2jset angle 60 trace 1 gap 2 spacing 1・732 origin 1.0 range jreg 2jset angle 60 trace 1 gap 2 spacing 1・732 origin. 1.1.732 range jreg 2Joints that are joined are assigned high values for normal and shear stiffness. It is tempting to give VCQ high values for stiffnesses to prevent movement along a fictitious joint. 1 lowcvcr, the timestep calculation in UDEC is based upon stiffnesses; the response (and solution convergence) will be very slow if very high siifTncsses are spexified. Tlic lowest MifTiiess consisieni with small joint deformaiion is used when joints are joined. The rule-of-thumb is that joint stifTnesses t I n and k s. should be set to a factor times the equivalent stiffness of the stiffest neighboring zone. The equation to calculate fictitious joint normal stiffness related to the equivalent stiffness (expressed in stress-per-distance units) of a zone in the normal direction is of the formk n = f actor x max(3.1) Figure 3.14 Zone dimension used ui stiffness calculationSystems of Units2.S Systems of L'nilsUDEC accepts any con listen f of engineering units, Examples of consis：enr sets of units forh阳it parameters are shown in Table 2.5. The user should be very careful when converting from one stem of uniK to unoiher. An excellent reference on the subject or units and conversion between the Imperiikl and SI s\ stems can be iimnd in the Juunuti </J Petndeuni lec /iimlti^y (DecCEiiber 1977). No conveniions are perfomicd in L DEC except for firictiun jnd diLatiun angles, M hich arc entered in degrees.Table 2.5 Sysleuta of units - mechanical parametersSI Imperial【工哩ih m m m cm ft inDen HQ kg/Hi'111? kg/in1IO6 kg/m1tlugt/ft3di bJ in"Foree N kN MN Mdyfies 呵IbfStress內kPa ME bar P点(Iraviij 7m/scc m/sec cnJs*ft/sec^in/sec^m/secwhere 1 bar = 10^ dynes! cm- = 1()^ N/ni- = 1(P Pa;1 atm = 1,013 bars = 14.7 psi = 2U6 lb(/ft2= 1.01325 x 10s Pa；1 slug = 1 Lbj —s- / ft = 14.59 kg；I snail = 1 Ib|- —in: andI gravity = 9.81 ni/s: = 9S1 cm/s2 =32.17 ft/s3.Histhist solve type 1solve rat 1e-5save tun l.savhist ydis 42 -67 sxx 42 -67hist ydis 42 -73 sxx 42 -73 reset disp jdispsolve rat 1e-5save tun 2.savPlotpl syy fill projected pl syy fill bfillJOS 1TTLE :UDEG \VwsiQfiFi却球 X27 Stress ertnfou rs using PLOT syy fill projectedpl ydisp fillpl bl id blackplot pen stress block显示应力pen作用是udec4.0窗口显示pl se ;看裂隙plsig1Lontours of major pnncipal shvss in the 工> -plant; (delunnable blocks only)sig2contours of minor principal stress in the xy-pl tine (deformable b kicks only)slip joints that are at limilin耳fricLioncontours of niuxiTTiuin prin^ipuJ slress in tbe XY-plane (defomiuble blocks only).Because 字亠皿□尸y鰹严.this is themiQ^^rincipiil stress I sm^lc = gig2，.J contours of ininimum pnriLiptd stress in the A v-plLint (deformublt 丿bkeks only).Because compressive stresses nrc negative, (his is the m ajar principal slrciis I smln = $ig1 )・ssi cmil门山*。

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 命令用来初始化应力。

udevadm命令详解udevadm 后接⼀个命令和命令指定选项。

它控制了udev运⾏的⾏为，处理内核事件，控制事件队列，并且提供简单的调试机制。

选项：--debug 打印错误信息--version 打印版本信息--help 帮助⽂档udevadm info options 查询udev数据库中的设备信息。

也可以从sysfs⽂件系统中查询到设备的属性以辅助创建udev规则。

--query=type 从数据库中查询指定类型的设备。

需要--path和--name来指定设备。

合法的查询⽂件是：设备名，链接，路径，属性--path=devpath 设备的路径--name=file 设备节点或者链接--attribute-walk 打印指定设备的所有sysfs记录的属性，以⽤来udev规则匹配特殊的设备。

该选项打印链上的所有设备信息，最⼤可能到sys⽬录。

--device-id-of-file=file　打印主/从设备号--export-db 输出udev数据库中的内容udevadm trigger [options] 接收内核发送来的设备事件。

主要⽤于重放coldplug事件信息（译者补充：内核在启动时已经检测到了系统的硬件设备，并把硬件设备信息通过sysfs内核虚拟⽂件系统导出。

udev扫描sysfs⽂件系统，根据硬件设备信息⽣成热插拔（hotplug）事件，udev再读取这些事件，⽣成对应的硬件设备⽂件。

由于没有实际的硬件插拔动作，所以这⼀过程被称为coldplug。

）--verbose 输出将要被触发的设备列表。

--dry-run 不真的触发事件--type=type 触发⼀个特殊的设备。

合法的类型：devices,subsystem,failed.默认是devices--action=action 被触发的事件，默认是change--subsystem-match=subsystem 触发匹配⼦系统的设备事件。

从CAD中提取坐标建UDEC命令流方法第一步在cad中选中一条（或多条）线段或多段线第二步在cad命令栏里输入“list”命令按回车键执行，如果点过多一次输出不完可根据窗口提示按enter键直至结束第三步复制窗口中的坐标到ultraedit文档编辑器中使用列模式去掉多余信息即可;;;;;;;;;;;查看模型基本命令plot block ;查询块体pl bl zone ;查询计算网格pl bl boundry;检查边界plot bound xcond ;查看x方向边界应力plot bound ycond ;查看y方向边界应力pl bl model;查看材料模型pl bl mat number ; 查看参数材料号pl bl mat block ; 查看块体材料pl bl mat joint ; 查看节理材料pl bl histry ; 查看预设的记录位移的点pl bl densty ; 按不同密度显示材料赋值结果plot block number ;显示包括块体地址号的块体图形plot contact num ;显示接触面编号;;;;;;计算结果查看常用命令RESET disp jdisp ;重设位移pl bl jdis ;查看节理位移plot block dis(yel) ;显示块体的位移矢量(有x和y方向)后加颜色可以自定义输出图形的颜色命令中不要括号，其他依次类推pl hist unbal ;显示最大不平衡力pl hist 1;显示输出第一个记录点的数据pl hist 2;显示输出第二个记录点的数据plot block velcity ;显示块体的速度矢量plot block pl ;显示块体的塑性区(plastic)plot block slip lmag open yel ;识别已经发生滑动（粗线）或张开（黄色）的节理长度plot block contact;查看块体本身接触PLOT overlap;查看块体嵌入pl pp ;显示渗流压力pl xdisp fill ;显示x方向位移云图pl ydisp fill ;显示y方向位移云图pl syy fill ;显示竖向应力云图，去掉fill即为应力等值线图pl sxx fill ;显示x方向应力云图pl bl id black;显示块体id号，黑色plot block stress ;显示块体的应力plot block ccs1 ; 显示最大主应力图plot block ccs2 ; 显示最小主应力图plot pen stress block;显示应力pen作用是udec4.0窗口显示plot pen block displace;显示位移plot block shear yel ;显示节理面的错动变形量（黄色）pl bl slip 显示滑动plot bl sh yel ;显示剪切位移pl bl syy fill ;显示y方向正应力云图pl bl xd fill ;显示x方向位移云图plot block cable red supp ye stru bl;显示块体、锚杆（红色）、支架（黄色）、梁（蓝色）的图形pause;暂停计算;动画，先输出图再开始记录输出图的图，再设定文件名，再设定记录步长plot block; 输出图形movie on ;动画开始movie file=kala.dcx;设置动画文件名movie step 500;动画记录步长VORONOI生成器产生随机大小的多边形块体VORONOI命令具有下列形式：VORONOI edge l < iterations n> <round v> < range … >vor edge 1.0 iter 100 range mat 1VORONOI e l <i n><r v><ra >newround 0.01bl 0 0 0 30 10 30 10 0cr 0 10 10 10cr0 20 10 20ch mat 3 range 0 10 0 10ch mat 2 range 0 10 10 20vor edge 2.0 range mat 3; delete stored joints before adding JSET jointsjreg id=1 (0,0) (0,30) (10,30) (10,0) deletejset 45 0 50 0 0 0 1 0 range mat 2动画文件默认保存在安装目录下reset disp; 位移清零，重设位移reset hist; 历史清零，重设历史记录;应变plot ssi/ssr/vsi/vsr;应力plot ccs1 fill bl; plot sxx/syy/szz/sxy/sxz/syz fill; plot sig1/sig2/ fill主应力;列表显示plot ccs1/ccs2应力矢量; print hist; plot; gp;速度plot xvel/yvel/vel;变形plot xdisp/ydisp/disp fill检查plot model块体本构模型;plot mat block/joint;block: plot property bulk/fric...(keyword)用照片形式输出想要的图形（窗口方式）1、首先在窗口中显示出图形2、点击file 中的print plot setup，选中弹出窗口中的JPEG3、点击file中的print plot即可，输出的图形在udec的安装目录下左边begin，end要设置数据输出自己想要的cad图形（命令方式）Set plot dxf 256pl bl 你想要的图Copy name.dxf1.把图形保存下来，能在AUTOCAD中打开、编辑plot block cable red supp ye stru bl;显示块体、锚杆（红色）、支架（黄色）、梁（蓝色）的图形set plot 256 dxf；设置图形为256色set out c:\ss.dxfcopy c:\ss.dxf；把这图形以ss.dxf文件保存在C磁盘下2.把数据导出set log onprint pline 1 ydis;把pline 1的y方向位移的数据导出到UDEC目录下udec.txt文件中set log off请教，如何将udec 3.1中的每一层设置为不同的颜色？我试了一下，plot block prop dens可以运行，其效果是按密度将不同的岩层划分；相同的，你也可以输入plot block prop fric 按角度来进行不同层的划分PLOT plas命令显示应力满足屈服条件的单元set pline x1 y1 x2 y2 nplot pline n <keyword>其中：x1 y1为监测线的起点坐标，x2，y2为监测线的终点坐标，第一个n为监测线划分的段数，第二个n为监测线的编号。

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=-.03 0 -.03insitu stress=-24.0, 0.0, -24.0 ygrad=-.03 0 -.03bound 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 21reset 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. 假定衬砌为线弹性材料。