udec命令总结精华正宗.docx

1、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:111kaicai.sav, 注意其

2、中111文件必须提前建好，否则无法保存或者保存错误(4) 操作：Ctrl+Z选中图像可以放大，Ctrl+Z 双击复原，屏幕中会出现十字叉，按住鼠标左键不放，移动光标直到你满意的窗口为止；pause暂停，此时可以察看任何信息；continue 继续调用下面程序段。ESC可以随时进行停止，但不能继续；英文分号 ; 表示注释不运行命令。2. 基本命令2.1 基本设置Udec n 【new刷新窗口，从新调用一个程序，修改后的*.txt文件必须输入n，重新运行文本文件】Udec title 【或heading代表标题，后面紧跟标题的名称。如：hang dao mo ni 或济宁三号井围岩变形破坏规律研究

3、】Udec round d 【 “圆角”命令，Udec中所有的块体都有圆角，目的是为防止块体悬挂在有棱角的节点上，由于块体悬挂将产生应力集中。d指块体与块体之间的圆角半径，默认值是0.5，其值要求小于模型中最小块体的最短那条边长的二分之一，最大圆角长度不能超过块体平均棱长的1%。在block命令前指定圆角长度。如：round 0.05】图1 1010块体圆角Udecset ovtol 0.5 【此命令是指层与层之间的嵌入厚度，当提示为“overlap too large”时就需要修改此值更大一些，可以显示设计的块体，plot overlap！嵌入太大的原因可能为块体强度太小】Udecset l

4、og on print set log off 【命令用于导出数据，可以将数据导出至硬盘f:】命令：set log f:yuanyan.log print 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 【建立模型框架，四个坐标角点必须按“顺时针”排列，也可以为五个坐标点。】

5、图2 5个坐标点 pl bl显示块体；pl hold bl =pl bl +pause屏幕上直接显示划分的块体，关闭显示后将继续运行； del range x1 x2 y1 y2删除块体；fix range x1 x2 y1 y2 固定块体不动;画圆弧Udecarc xc,yc xb,yb theta ncracks【xc,yc圆弧中心，xb,yb指画圆弧的起始点，theta沿逆时针方向圆弧的角度;ncracks指构成指圆弧的折线段数，数目越多，圆弧越圆； 圆弧命令必须配合Crack命令，但tunnel 命令可以在输入后直接显示为圆形；先划分巷道，然后划分其它岩层节理】命令： ;开挖巷道节理，

6、形成巷道空间,巷道必须先开挖出来arc 52.17 85.59 53.97 85.59 180 8cra 50.37 83.59 50.37 85.59 cra 50.37 83.59 53.97 83.59 cra 53.97 83.59 53.97 85.59;建圆形隧道、圆形巷道Udectunnel 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方法

7、三：指定窗口 windows x1 x2 y1 y2 del range window，不一定能全部删除，可以继续删除编号方法四：定义区域，删除区域，jreg n x1,y1 x2,y2 x3,y3 x4,y4 delete jreg n】;节理划分、定义Udeccrack x1,y1 x2,y2【划分节理后，最好输入pause暂停，pl bl显示划分的块体是否正确或文本中输入pl hold bl】Udecjregion id n x1,y1 x2,y2 x3,y3 x4,y4 【节理区域定义，赋予节理区域地址号，可以删除任何区域delete jreg n 】命令：jset 90,0 4,0

8、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 ;定义节理区域4Udecchange mat 4 range 0 100 -65 -50 ;定义Udecvor edg 5 iter 20 rou 0.1 ran jreg 4【edge平均边长，最少20倍rounding长度；iterations 迭代次数，默认5；ro

9、und，产生块体的圆角长度，默认0.5；rtol 容许偏差，默认0.02*round；seed n】2.3 参数定义gen quad 10 range xl xu yl yu;在指定的区域生成一定宽度的单元zone model mo range xl xu yl yu;使指定的区域材料采用摩尔-库仑本构关系计算（即弹塑性）change jcons=2 range xl xu yl yu;使指定的区域节理遵循摩尔-库仑准则计算（即弹塑性）change mat=1 range xl xu yl yu change mat=2 range xl xu yl yu change mat=3 range

10、 xl xu yl yu ;指定各岩层的材料标号change jmat=1 range xl xu yl yu change jmat=2 range xl xu yl yu change jmat=3 range xl xu yl yu ;指定各岩层的节理标号prop mat=1 dens=2000 prop 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=1

11、zone 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,jte

12、n=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，定义范围，四点坐标，del

13、ete为常规语言 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.30.3826.87.0粉砂岩26.30.2215.610.8石灰岩28.50.2922.611.1页岩、泥岩11.10.298.84.3大理岩55.80.2537.222.3花岗岩73.80.2243.930.2岩石名称 内摩擦角

14、粘聚力抗拉强度砂岩27.827.21.17粉砂岩32.134.7页岩、泥岩14.438.4石英岩48.070.6石灰岩42.06.721.58花岗岩51.055.1玄武岩31.066.213.1体积模量K和剪切模量G之间的关系： gen 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

15、 0,0 0,15 300,15 300,0 ;定义物质1的范围，即赋予这个范围为物质1,只赋名，没有定义属性change mat=2 range reg 0,15 0,20 300,20 300,15change mat=3 range reg 0,20 0,23 300,23 300,20change 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 ;定义节理属性，使符合摩尔库伦准则普适。2代表摩尔库伦准则 ;参数设置prop

16、 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;定义节理的属性法向刚度、切向刚度、内摩擦角、粘聚力、抗拉强度1) D-质量密度 F-摩擦角 B-体积模量 C-内聚力(粘结力) S-剪切模量2) 应力-正号代表张力，负号代表压力； 应变-正的应变代表伸长，负的应变代表压缩 重力-正号的重力物质往下拉，负号的

17、重力将物质往上提密度必须设置density for block materical 1 is zero cannot cycle-块体密度未设2.4 初始应力set gravity 0,-9.81;设置重力加速度，x方向为0，y方向为-9.8bound xvel=0 range -0.1 3.00 -60.1 20.1 bound xvel=0 range 97 100.1 -60.1 20.1 bound yvel=0 range 0.1 100.1 -60.1 -58 ;采用位移法固定边界，让其速度=0。;边界条件boun stress 0 0 0 range 0,300 99.9 100

18、.1 ;定义边界条件，三向应力为0，上边界，范围，纵向波动范围0 0 0（x、xy、y）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应力梯度bound xvel=0 range -0.1 0.1 0 100 ;X位移边界,位移波动范围，范围bound xvel=0 range 2

19、99.9 300.1 0 100 ;X位移边界，同上bound yvel=0 range -0.1 300.1 -0.1,0.1 ;Y位移边界hist solve_ratio type 1 ;设置不平衡力的精度，普适hist xdisp 10,20 ;测该点的X方向位移hist sxx 10,20 ;测该点X方向的应力hist syy 10,20 ;测该点Y方向的应力reset dispreset hist ;位移清零，历史记录清零his unbalsolve ratio 1e-5 ;记录不平衡力，计算直到精度达到1e-5;执行计算solvestep 5000cycle 5000;进行开挖;删

20、除水平块体，开挖，不可再充填delete range x1 x2 y1 y2等价于delete range reg x1,y1 x1,y2 x2,y2 x2,y1;删除不规则四边形delete range region x1,y1 x2,y2 x3,y3 x4,y43.5 支护部分;加支柱，支架support x,y angle segment witdth mat;support是指加入支柱，x,y表明支柱所处的位置点。angle表明支柱的轴向与x轴正方向的夹角，以逆时针方向为正。segment指一组支柱的个数，默认值为5。witdth指一组支柱支撑的范围。;mat是指赋予支柱的材料性质属于

21、几号？下面的table指支柱的增阻特性。比如：sup 10,20 angle=90 seg=5 width=2 mat=6prop mat=6 sup_kn=-1 table 1 0 3e6 0.005 5e6 0.015 25e6 0.5 40e64.液压支架的命令supp xc yc wid l seg n mat j; (xc,yc)是指中心点坐标,l指支架的宽度,n指分段数,mat j指支柱材料性质为jprop mat j sup_kn -1 ;sup_kn指支架的刚度,-1与表1相对应del range x1 x2 y1 y2 ;挖掘范围(x1, y1)(x1 ,y2)(x2, y2

22、)(x2, y1)的块体,由支架支撑table 1 0 42.e6 0.05 5.0e6 0.1 6.0e6; 表1表示的是液压支架的(P-DS)特性曲线，后面为四个点的坐标。目前，液压支柱的最大支撑力为12*106N显示支架：pl bl support red ;加锚杆cable x1,y1 x2,y2 npoint mats asteel matg 等;x1,y1 x2,y2指锚杆的两个端点;npoint指锚杆的粘结点数;mats指锚杆材料性质,材料编号，按着前面给定的材料号继续往下编； asteel指锚杆横截面积；matg指锚固体性质；cable (hc1bx,hc1by) (hc1ex

23、,hc1ey) 25 9 314e-6 10 cable (hc2bx,hc2by) (hc2ex,hc2ey) 25 9 314e-6 10 ;锚杆参数prop m=9 cb_dens 7500 cb_ycomp 6.3e8 cb_yield 0.5e7 cb_ymod 0.5e11;锚固体的参数prop m=10 cb_kbond=1.6e9 cb_sbond=2e6 cable (73.4,2.91) (71.15,4.99) 10 12 380e-6 11 ;锚杆，起始点坐标，后面四个数值可固定prop mat 12 cb_dens 7800 cb_ycomp 430e6 cb_yie

24、ld 160e3 cb_ymod 1.3e11;锚固体的性质prop mat 11 cb_kbond 6.3e9 cb_sbond 6e5 ;定义锚杆的属性，密度、抗拉强度、屈服强度；巷道断面加梯子梁的命令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.3

25、5e9 if_ks=1.35e9 if_tens=0 if_fric=18 if_coh=0;interface-界面struct gen xc=20 yc=40.24 np=100 fang=-40 theta=260 mat=14 thick=0.2prop mat=14 st_ymod=13.5e9 st_prat=0.14 st_den=7800prop mat=14 st_yield=6e7 st_yresid=6e7 st_ycomp=2.5e7prop mat=14 if_kn=1.35e9 if_kf=1.35e9 if_tens=0 if_fric=18 if_coh=0图形

26、显示pl bl hist 1,2,. 观测点位移pl bl nu(mber) 显示块体标号pl bl cab red stru red cab-锚杆 stru-梁pl bl sxx/syy/szz 显示应力图pl bl dis 位移矢量图pl bound pl 显示塑性区pl bound dis 显示位移矢量图pl hold bl num 在屏幕上将显示建立的块体单元，并在按键后退出后，或直接继续执行命令 。相当于pau+pl bl num；删除块体 del bl +地址号plot block stress ;显示块体的应力plot block xdis ;显示块体的位移(有x和y方向)plo

27、t block pl ;显示块体的塑性区(plastic);屏幕观看下沉与应力等值线,int是指intervalpl bl ydisp int 0.1pl bl syy int 1e6udec 颜色关键词颜 色关键词颜 色black黑lgray浅灰blue蓝lblue浅蓝green绿lgreen浅绿cyan蓝绿lcyan浅蓝绿red红lred浅红magenta品红lmagenta浅品红orange橙lorange浅橙brown褐yellow黄dgray灰white白图形输出;抓取图像，在autocad内打开、编辑;输出应力分布pl bl pline 1 syy redset pl dxf 25

28、set out f:syy100.dxfpl pen bl pline 1 syy red ；=copy f:syy100.dxf1.把图形保存下来，能在AUTOCAD中 plot block cable red supp ye stru bl;显示块体、锚杆（红色）、支架（黄色）、梁（蓝色）的图形 set plot 256 dxf；设置图形为256色 set out c:ss.dxf copy c:ss.dxf；把这图形以ss.dxf文件保存在C磁盘下注：开始为dxf格式，移动后另存改为dwg格式，可以减少很大的空间！动画制作;打开图象抓拍功能，设置图片间隔时间以及窗口大小movie on

29、movie step 1000window x1 x2 y1 y2;movie off【注：只显示改区间内的内容,分别为x1，x2和y1，y2的坐标！window 6 27.5 -1 13pl bl syy movie file kaicai size 800 600 on step 1000del ran 40 70 12.0 15step 10000movie off；电影将被保存至C:Program FilesItascaudec310下面， 打开方式C:Program FilesItascasharedUtilitymovie，必须点击start才开始播放;【注：1)：电影的制作过程实

30、际上为抓频的过程，因此，在此期间最好将其它的程序先关掉！ 2): 在进行电影之前，输入pl bl syy,则整个过程将显示应力的变化过程，如果输入pl bl pline syy red则整个过程将显示测线应力变化的过程。输入先输入pl bl syy 后输入pl bl pline syy red，则将在电影中同时出现！输入pl bound，则将出现裂隙带的范围！4 实例分析: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