udec命令总结精华-正宗
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命令大全黑客必知的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常用命令集锦
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。
(完整word版)UDEC模拟实例与解析
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. 假定衬砌为线弹性材料。
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模拟[技巧]](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所示的数值计算模型。
UDEC指令
set gravity 0,-10 ;重力加速度
cr 0,15 300,15 ;画直线，两点坐标
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 ;定义节理角度90°；节理的长度；节理的距离，即纵向间隔；横向距离，即横向隔5一个；起始点坐标为0,0
change mat=3 range reg 0,20 0,23 300,23 300,20
change mat=4 range reg 0,23 0,28 300,28 300,23
change jmat=1 range reg 0,0 0,15 300,15 300,0 ;定义节理1的范围
model null range ;变性删除这个范围，即之后可以充填物质，X的范围，Y的范围
model m range ;充填这个范围，
solve ;计算
;new
title
moni ;定义名称moni
round=0.025
set ovtol=0.4 ;常规设置
bl 0,0 0,100 300,100 300,0 ;定义范围，四点坐标，顺时针方向
set jmatdf 2 ;定义节理属性，使符合摩尔库伦准则普适
save ch.sav ;保存文件
《udec版本整理版》PPT课件
UDEC常见问题
USTB
1
运算速度缓慢
2
自动退出运算
3
节理生成失败
University of Science & Technology Beijing
UDEC软件应用－－边界条件命令总结
USTB
命 令
BOUNDARY Stress Xload
效
果
施加总应力到刚体或变形体块体的边界 上 施加刚体或变形体边界的x方向的荷载 施加刚体或变形体边界的y方向的荷载 施加变形体边界的x方向的速度（位移） 施加变形体边界的y方向的速度（位移） 固定刚体边界的速度（位移）
◆ 美国Itasca Consulting Group Inc. 开发的离散单元法的计算软件。
◆ 离散单元法理论由Cundall（1971）首次提出。 ◆ 岩体被看作一种不连续的离散介质。其内部可存在大位移、旋转和滑动乃至块体 的分离，从而可以较真实地模拟节理岩体中的非线性大变形特征。
UDEC简介－应用范畴
Yload
Xvel Yvel FIX
FREE
LOAD Xload Yload
释放刚体的速度（位移）
施加x方向的荷载到刚体的边界 施加y方向的荷载到刚体的边界
UDEC软件应用－－建模
USTB 划分范围： Jr id n x1,y1 x2,y2 x3,y3 x4,y4
模型内应力： ini sxx range x1,x2 y1,y2
注意：量纲统一
UDEC软件应用－－建模
USTB
边界条件：Boundary xvel=0 range x x1,x2 Boundary yvel=0 range y y1,y2
UDEC软件应用－－建模
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《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【建立模型框架，四个坐标角点必须按“顺时针”排列，也可以为五个坐标点。
】图2 5个坐标点✈pl bl显示块体；pl hold bl =pl bl +pause屏幕上直接显示划分的块体，关闭显示后将继续运行；del range x1 x2 y1 y2删除块体；fix range x1 x2 y1 y2固定块体不动;画圆弧Udec>arc xc,yc xb,yb theta ncracks【xc,yc圆弧中心，xb,yb指画圆弧的起始点，theta沿逆时针方向圆弧的角度;ncracks指构成指圆弧的折线段数，数目越多，圆弧越圆；✈①圆弧命令必须配合Crack命令，但tunnel 命令可以在输入后直接显示为圆形；②先划分巷道，然后划分其它岩层节理】命令：;开挖巷道节理，形成巷道空间,巷道必须先开挖出来arc 52.17 85.59 53.97 85.59 180 8cra 50.37 83.59 50.37 85.59cra 50.37 83.59 53.97 83.59cra 53.97 83.59 53.97 85.59;建圆形隧道、圆形巷道Udec>tunnel xc,yc r n【xc,yc 圆的中心，r 半径，n 构成圆的线段数。
如：tunnel 4,5 2 45①先输tunnel ，然后输crack ，裂隙不贯穿隧道;②隧道开挖：方法一：删除编号pl bl num → delete block 1920 方法二：delete range -1 1 -1 1方法三：指定窗口 windows x1 x2 y1 y2 →del range window ，不一定能全部删除，可以继续删除编号方法四：定义区域，删除区域，jreg n x1,y1 x2,y2 x3,y3 x4,y4→ delete jreg n 】;节理划分、定义 Udec>crack x1,y1 x2,y2【划分节理后，最好输入pause 暂停，pl bl 显示划分的块体是否正确或文本中输入pl hold bl 】 Udec>jregion id n x1,y1 x2,y2 x3,y3 x4,y4【节理区域定义，赋予节理区域地址号，可以删除任何区域delete jreg n 】 命令：jset 90,0 4,0 4,0 6,0 0,-50 range jreg n【在定义的节理区域内进行节理的再划分;jset α、t 、g 、s (xo,yo) n :倾角，节理长度, 岩桥长度,垂距，起始点坐标,n 区域标号】图3 节理组定义；不规则块体划分Udec> jreg id 4 14.65 0 73.0 160 78.16 160 19.96 0 ;定义节理区域4 Udec>change mat 4 range 0 100 -65 -50 ;定义 Udec>vor edg 5 iter 20 rou 0.1 ran jreg 4【edge 平均边长，最少20倍rounding 长度；iterations 迭代次数，默认5；round ，产生块体的圆角长度，默认0.5；rtol 容许偏差，默认0.02*round ；seed n 】 2.3 参数定义gen quad 10 range xl xu yl yu ;在指定的区域生成一定宽度的单元αg tszone 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(kg/m3)，即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为节理的抗拉强度块体力学参数k--体积模量；g--剪切模量；fric--摩擦角；coh--粘结力（内聚力）；ten--抗拉强度；d--体积力接触面力学参数jkn--法向刚度；jks--切向强度；jc--粘结力；jf--摩擦角；jt--抗拉强度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 ;定义节理角度90°；节理的长度；节理的距离，即纵向间隔；横向距离，即横向隔5一个；起始点坐标为0,0jset 90,0 5,0 5,0 5,0 2.5,5 rang jreg 1 ;定义节理的另外一项表1 岩石块体参数岩石名称 弹性模量E 泊松比υ 体积模量K 剪切模量G 砂岩 19.3 0.38 26.8 7.0 粉砂岩 26.3 0.22 15.6 10.8 石灰岩 28.5 0.29 22.6 11.1 页岩、泥岩 11.1 0.29 8.8 4.3 大理岩 55.8 0.25 37.2 22.3 花岗岩73.80.2243.930.2岩石名称 内摩擦角 粘聚力 抗拉强度 砂岩 27.8 27.2 1.17 粉砂岩 32.1 34.7 \ 页岩、泥岩 14.4 38.4 \ 石英岩 48.0 70.6 \ 石灰岩 42.0 6.72 1.58 花岗岩 51.0 55.1 \ 玄武岩31.066.213.1体积模量K 和剪切模量G 之间的关系：)21(3υ-=E K )1(2υ+=E G υυ+-=1)21(5.1K Ggen quad 10 range 0 300 0 100 ;当0 300 0 100范围内的块体的对角长度小于10时就是可变形体。
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,15 change mat=3 range reg 0,20 0,23 300,23 300,20 change mat=4 range reg 0,23 0,28 300,28 300,23change jmat=1 range reg 0,0 0,15 300,15 300,0 ;定义节理1的范围set jmatdf 2 ;定义节理属性，使符合摩尔库伦准则普适。