《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 ;定义节理属性,使符合摩尔库伦准则普适。