udec命令总结精华正宗.docx

udec命令总结精华正宗.docx

《udec命令总结精华正宗.docx》由会员分享，可在线阅读，更多相关《udec命令总结精华正宗.docx（51页珍藏版）》请在冰点文库上搜索。

udec命令总结精华正宗

《Udec命令总结》精华版

1.安装、打开、保存

（1）安装：

①执行Udec3.1→将Crack文件中的内容替换；

②复制Udec.exe,粘贴为快捷方式→属性，目标，加入空格256【开始内存8M，将内存改为256】，可以复制快捷方式至桌面或硬盘。

（2）打开：

①直接在udec>命令行输入；

②写好程序，udec>命令输入call，然后将*.txt文件拖入命令行，执行

（3）保存：

输入saved:

\kaicai.sav,调用命令rest，将kaicai.sav拖入命令行

【文件名最好不用汉字，有时候不识别】

注释：

如果保存为saved:

\111\kaicai.sav,注意其中111文件必须提前建好，否则无法保存或者保存错误

（4）操作：

Ctrl+Z选中图像可以放大，Ctrl+Z双击复原，屏幕中会出现十字叉，按住鼠标左键不放，移动光标直到你满意的窗口为止；pause暂停，此时可以察看任何信息；continue继续调用下面程序段。

ESC可以随时进行停止，但不能继续；英文分号;表示注释不运行命令。

2.基本命令

2.1基本设置

Udec>n【new刷新窗口，从新调用一个程序，修改后的*.txt文件必须输入n，重新运行文本文件】

Udec>title【或heading代表标题，后面紧跟标题的名称。

如：

hangdaomoni或济宁三号井围岩变形破坏规律研究】

Udec>roundd【“圆角”命令，Udec中所有的块体都有圆角，目的是为防止块体悬挂在有棱角的节点上，由于块体悬挂将产生应力集中。

d指块体与块体之间的圆角半径，默认值是0.5，其值要求小于模型中最小块体的最短那条边长的二分之一，最大圆角长度不能超过块体平均棱长的1%。

在block命令前指定圆角长度。

如：

round0.05】

图110×10块体圆角

Udec>setovtol0.5【此命令是指层与层之间的嵌入厚度，当提示为“overlaptoolarge”时就需要修改此值更大一些，可以显示设计的块体，plotoverlap！

嵌入太大的原因可能为块体强度太小】

Udec>setlogon

print

setlogoff【命令用于导出数据，可以将数据导出至硬盘f:

\】

命令：

setlogf:

\yuanyan.log

printpline1syy【记录测线1应力】

printpline2ydis

setlogoff

;设置观测线

setpline1x1,y1x2,y2n（n--观测线分的段数）

setpline25,20300,2010;定义测线，起始点坐标，测线分段

2.2图形划分

;块体命令

Udec>blockx1,y1x2,y2x3,y3x4,y4

【建立模型框架，四个坐标角点必须按“顺时针”排列，也可以为五个坐标点。

】

图25个坐标点

plbl显示块体；plholdbl=plbl+pause屏幕上直接显示划分的块体，关闭显示后将继续运行；delrangex1x2y1y2删除块体；fixrangex1x2y1y2固定块体不动

;画圆弧

Udec>arcxc,ycxb,ybthetancracks

【xc,yc圆弧中心，xb,yb指画圆弧的起始点，theta沿逆时针方向圆弧的角度;ncracks指构成指圆弧的折线段数，数目越多，圆弧越圆；

①圆弧命令必须配合Crack命令，但tunnel命令可以在输入后直接显示为圆形；

②先划分巷道，然后划分其它岩层节理】

命令：

;开挖巷道节理，形成巷道空间,巷道必须先开挖出来

arc52.1785.5953.9785.591808

cra50.3783.5950.3785.59

cra50.3783.5953.9783.59

cra53.9783.5953.9785.59

;建圆形隧道、圆形巷道

Udec>tunnelxc,ycrn

【xc,yc圆的中心，r半径，n构成圆的线段数。

如：

tunnel4,5245

①先输tunnel，然后输crack，裂隙不贯穿隧道;

②隧道开挖：

方法一：

删除编号plblnum→deleteblock1920

方法二：

deleterange-11-11

方法三：

指定窗口windowsx1x2y1y2→delrangewindow，不一定能全部删除，可以继续删除编号

方法四：

定义区域，删除区域，jregnx1,y1x2,y2x3,y3x4,y4→deletejregn】

;节理划分、定义

Udec>crackx1,y1x2,y2

【划分节理后，最好输入pause暂停，plbl显示划分的块体是否正确或文本中输入plholdbl】

Udec>jregionidnx1,y1x2,y2x3,y3x4,y4

【节理区域定义，赋予节理区域地址号，可以删除任何区域deletejregn】

命令：

jset90,04,04,06,00,-50rangejregn

【在定义的节理区域内进行节理的再划分;

jsetα、t、g、s（xo,yo）n:

倾角，节理长度,岩桥长度,垂距，起始点坐标,n区域标号】

图3节理组定义

；不规则块体划分

Udec>jregid414.65073.016078.1616019.960;定义节理区域4

Udec>changemat4range0100-65-50;定义

Udec>voredg5iter20rou0.1ranjreg4

【edge平均边长，最少20倍rounding长度；iterations迭代次数，默认5；round，产生块体的圆角长度，默认0.5；rtol容许偏差，默认0.02*round；seedn】

2.3参数定义

genquad10rangexlxuylyu

;在指定的区域生成一定宽度的单元

zonemodelmorangexlxuylyu

;使指定的区域材料采用摩尔--库仑本构关系计算（即弹塑性）

changejcons=2rangexlxuylyu

;使指定的区域节理遵循摩尔--库仑准则计算（即弹塑性）

changemat=1rangexlxuylyu

changemat=2rangexlxuylyu

changemat=3rangexlxuylyu

;指定各岩层的材料标号

changejmat=1rangexlxuylyu

changejmat=2rangexlxuylyu

changejmat=3rangexlxuylyu

;指定各岩层的节理标号

propmat=1dens=2000

propmat=2dens=2650

propmat=3dens=2700

;指定各材料的密度，比如1号材料dens=2000（kg/m3），即1立方米重2吨

zonek=0.15e9,g=0.1e9,fric=10.00,coh=0.19e6,ten=0.09e6rangemat=1

zonek=2.8e9,g=2.2e9,fric=30.00,coh=1.5e6,ten=0.4e6rangemat=2

zonek=6.9e9,g=6.6e9,fric=38.62,coh=5.63e6,ten=3.20e6rangemat=3

;k为材料的法向刚度，g为材料的切向刚度，friction为材料的内摩擦角，cohesion为材料的内聚力，tension为材料的抗拉强度

propjmat=1jkn=0.2e8,jks=0.1e7,jcoh=0,jfric=4,jten=0

propjmat=2jkn=8e8,jks=5e7,jcoh=0.1e6,jfric=8,jten=0

propjmat=3jkn=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--抗拉强度

jregid10,00,15300,15300,0delete;定义节理，命名1，定义范围，四点坐标，delete为常规语言

jset90,05,05,05,00,0rangjreg1;定义节理角度90°；节理的长度；节理的距离，即纵向间隔；横向距离，即横向隔5一个；起始点坐标为0,0

jset90,05,05,05,02.5,5rangjreg1;定义节理的另外一项

表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.8

0.22

43.9

30.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.0

66.2

13.1

体积模量K和剪切模量G之间的关系：

genquad10range03000100;当03000100范围内的块体的对角长度小于10时就是可变形体。

zonemodelmorange03000100;定义摩尔库伦模型的范围，X的范围，Y的范围

changejcons=5range03000100;定义节理的某个属性

changemat=1rangereg0,00,15300,15300,0;定义物质1的范围，即赋予这个范围为物质1,只赋名，没有定义属性

changemat=2rangereg0,150,20300,20300,15

changemat=3rangereg0,200,23300,23300,20

changemat=4rangereg0,230,28300,28300,23

changejmat=1rangereg0,00,15300,15300,0;定义节理1的范围

setjmatdf2;定义节理属性，使符合摩尔库伦准则普适。

2代表摩尔库伦准则

;参数设置

propmat=1dens=2720;物质的密度

zonek=12.12e9g=10.26e9fric=33coh=5.77e6ten=3.6e6rangemat=1

;定义物质1的体积模量、剪切模量、内摩擦角、粘聚力、抗拉强度

propjmat=1jkn=11e8jks=12e7jfric=12jcoh=0jten=0

;定义节理的属性法向刚度、切向刚度、内摩擦角、粘聚力、抗拉强度

1）D--质量密度F--摩擦角B--体积模量C--内聚力（粘结力）S--剪切模量

2）应力--正号代表张力，负号代表压力；应变--正的应变代表伸长，负的应变代表压缩

重力--正号的重力物质往下拉，负号的重力将物质往上提

密度必须设置densityforblockmaterical1iszerocannotcycle----块体密度未设

2.4初始应力

setgravity0,-9.81

;设置重力加速度，x方向为0，y方向为-9.8

boundxvel=0range-0.13.00-60.120.1

boundxvel=0range97100.1-60.120.1

boundyvel=0range0.1100.1-60.1-58

;采用位移法固定边界，让其速度=0。

;边界条件

bounstress000range0,30099.9100.1;定义边界条件，三向应力为0，上边界，范围，纵向波动范围000（x、xy、y）

bounstress00-8.6095397e6range;定义边界应力条件，竖直方向的应力为，负号表示方向向下

insitustr-7.5e60-15e6szz-7.5e6ygrad1.1018e402.203607e4;定义应力属性，sxx方向即水平应力，sxy方向即剪切应力，syy方向即竖直应力，szz方向的应力，ygrad应力梯度

boundxvel=0range-0.10.10100;X位移边界,位移波动范围，范围

boundxvel=0range299.9300.10100;X位移边界，同上

boundyvel=0range-0.1300.1-0.1,0.1;Y位移边界

histsolve_ratiotype1;设置不平衡力的精度，普适

histxdisp10,20;测该点的X方向位移

histsxx10,20;测该点X方向的应力

histsyy10,20;测该点Y方向的应力

resetdisp

resethist;位移清零，历史记录清零

hisunbal

solveratio1e-5;记录不平衡力，计算直到精度达到1e-5

;执行计算

solve\step5000\cycle5000

;进行开挖

;删除水平块体，开挖，不可再充填

deleterangex1x2y1y2等价于deleterangeregx1,y1x1,y2x2,y2x2,y1

;删除不规则四边形

deleterangeregionx1,y1x2,y2x3,y3x4,y4

3.5支护部分

;加支柱，支架

supportx,yanglesegmentwitdthmat

;support是指加入支柱，x,y表明支柱所处的位置点。

angle表明支柱的轴向与x轴正方向的夹角，以逆时针方向为正。

segment指一组支柱的个数，默认值为5。

witdth指一组支柱支撑的范围。

;mat是指赋予支柱的材料性质属于几号？

下面的table指支柱的增阻特性。

比如：

sup10,20angle=90seg=5width=2mat=6

propmat=6sup_kn=-1

table103e60.0055e60.01525e60.540e6

4.液压支架的命令

suppxcycwidlsegnmatj;（xc,yc）是指中心点坐标,l指支架的宽度,n指分段数,matj指支柱材料性质为j

propmatjsup_kn-1;sup_kn指支架的刚度,-1与表1相对应

delrangex1x2y1y2;挖掘范围（x1,y1）（x1,y2）（x2,y2）（x2,y1）的块体,由支架支撑

table1042.e60.055.0e60.16.0e6;表1表示的是液压支架的（P-DS）特性曲线，后面为四个点的坐标。

目前，液压支柱的最大支撑力为12*106N

显示支架：

plblsupportred

;加锚杆

cablex1,y1x2,y2npointmatsasteelmatg等

;x1,y1x2,y2指锚杆的两个端点;npoint指锚杆的粘结点数;mats指锚杆材料性质,材料编号，按着前面给定的材料号继续往下编；asteel指锚杆横截面积；matg指锚固体性质；

cable（hc1bx,hc1by）（hc1ex,hc1ey）259314e-610

cable（hc2bx,hc2by）（hc2ex,hc2ey）259314e-610

;锚杆参数

propm=9cb_dens7500cb_ycomp6.3e8cb_yield0.5e7cb_ymod0.5e11

;锚固体的参数

propm=10cb_kbond=1.6e9cb_sbond=2e6

cable（73.4,2.91）（71.15,4.99）1012380e-611;锚杆，起始点坐标，后面四个数值可固定

propmat12cb_dens7800cb_ycomp430e6cb_yield160e3cb_ymod1.3e11

;锚固体的性质

propmat11cb_kbond6.3e9cb_sbond6e5;定义锚杆的属性，密度、抗拉强度、屈服强度

；巷道断面加梯子梁的命令

strugenxcycnp100faathetrabmat=16thick=0.2;100指分100段,ab指角度

propmat=16st_ymod=13.5e9st_prat=0.14st_den=7800

propmat=16st_yield=6e7st_yresid=6e7st_ycomp=2.5e7

propmat=16if_kn=1.35e9if_ks=1.35e9if_tens=0if_fric=18if_coh=0;interface-界面

structgenxc=20yc=40.24np=100fang=-40theta=260mat=14thick=0.2

propmat=14st_ymod=13.5e9st_prat=0.14st_den=7800

propmat=14st_yield=6e7st_yresid=6e7st_ycomp=2.5e7

propmat=14if_kn=1.35e9if_kf=1.35e9if_tens=0if_fric=18if_coh=0

图形显示

plblhist1,2,...观测点位移

plblnu（mber）显示块体标号

plblcabredstruredcab--锚杆stru--梁

plblsxx/syy/szz显示应力图

plbldis位移矢量图

plboundpl显示塑性区

plbounddis显示位移矢量图

plholdblnum在屏幕上将显示建立的块体单元，并在按键后退出后，或直接继续执行命令。

相当于pau+plblnum；删除块体delbl+地址号

plotblockstress;显示块体的应力

plotblockxdis;显示块体的位移（有x和y方向）

plotblockpl;显示块体的塑性区（plastic）

;屏幕观看下沉与应力等值线,int是指interval

plblydispint0.1

plblsyyint1e6

udec颜色

关键词

颜色

关键词

颜色

black

黑

lgray

浅灰

blue

蓝

lblue

浅蓝

green

绿

lgreen

浅绿

cyan

蓝绿

lcyan

浅蓝绿

red

红

lred

浅红

magenta

品红

lmagenta

浅品红

orange

橙

lorange

浅橙

brown

褐

yellow

黄

dgray

灰

white

白

图形输出

;抓取图像，在autocad内打开、编辑

;输出应力分布

plblpline1syyred

setpldxf25

setoutf:

\syy100.dxf

plpenblpline1syyred；=copyf:

\syy100.dxf

1.把图形保存下来，能在AUTOCAD中

plotblockcableredsuppyestrubl;显示块体、锚杆（红色）、支架（黄色）、梁（蓝色）的图形

setplot256dxf；设置图形为256色

setoutc:

\ss.dxf

copyc:

\ss.dxf；把这图形以ss.dxf文件保存在C磁盘下

注：

开始为dxf格式，移动后另存改为dwg格式，可以减少很大的空间！

动画制作

;打开图象抓拍功能，设置图片间隔时间以及窗口大小

movieon

moviestep1000

windowx1x2y1y2

;movieoff

【注：

只显示改区间内的内容,分别为x1，x2和y1，y2的坐标！

window627.5-113

plblsyy

moviefilekaicaisize800600onstep1000

delran407012.015

step10000

movieoff

；电影将被保存至C:

\ProgramFiles\Itasca\udec310下面，

打开方式C:

\ProgramFiles\Itasca\shared\Utility\movie，必须点击start才开始播放;

【注：

1）：

电影的制作过程实际上为抓频的过程，因此，在此期间最好将其它的程序先关掉！

2）:

在进行电影之前，输入plblsyy,则整个过程将显示应力的变化过程，如果输入plblplinesyyred

则整个过程将显示测线应力变化的过程。

输入先输入plblsyy后输入plblplinesyyred，则将在电影中同时出现！

输入plbound，则将出现裂隙带的范围！

4实例分析:

round=0.1;方块的圆角块半径为0.1米

setovtol=1.0;块体与块体之间相互嵌入量最大值为1米

bl0,00,2650,2650,0;在（x1,y1）（x1,y2）（x2,y2）（x2,y1）生成块体范围

crack0,2550,25;在（x1,y1）（x2,y2）两点间画直线

jregionid10,00,1050,1050,0delete;删除此区域的块体

jset90,02.5,02.5,03,00,0