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>setlogonprintsetlogoff【命令用于导出数据，可以将数据导出至硬盘f:\】命令：setlogf:\yuanyan.logprintpline1syy【记录测线1应力】printpline2ydissetlogoff;设置观测线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.591808cra50.3783.5950.3785.59cra50.3783.5953.9783.59cra53.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区域标号】ααgts图3节理组定义；不规那么块体划分Udec>jregid414.65073.016078.1616019.960;定义节理区域4Udec>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=1rangexlxuylyuchangemat=2rangexlxuylyuchangemat=3rangexlxuylyu;指定各岩层的材料标号changejmat=1rangexlxuylyuchangejmat=2rangexlxuylyuchangejmat=3rangexlxuylyu;指定各岩层的节理标号propmat=1dens=2000propmat=2dens=2650propmat=3dens=2700;指定各材料的密度，比方1号材料dens=2000(kg/m3)，即1立方米重2吨zonek=0.15e9,g=0.1e9,fric=10.00,coh=0.19e6,ten=0.09e6rangemat=1zonek=2.8e9,g=2.2e9,fric=30.00,coh=1.5e6,ten=0.4e6rangemat=2zonek=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=0propjmat=2jkn=8e8,jks=5e7,jcoh=0.1e6,jfric=8,jten=0propjmat=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,0jset90,05,05,05,02.5,5rangjreg1;定义节理的另外一项表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岩石名称内摩擦角粘聚力抗拉强度砂岩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之间的关系：genquad10range03000100;当03000100范围内的块体的对角长度小于10时就是可变形体。zonemodelmorange03000100;定义摩尔库伦模型的范围，X的范围，Y的范围changejcons=5range03000100;定义节理的某个属性changemat=1rangereg0,00,15300,15300,0;定义物质1的范围，即赋予这个范围为物质1,只赋名，没有定义属性changemat=2rangereg0,150,20300,20300,15changemat=3rangereg0,200,23300,23300,20changemat=4rangereg0,230,28300,28300,23changejmat=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.8boundxvel=0range-0.13.00-60.120.1boundxvel=0range97100.1-60.120.1boundyvel=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方向的应力resetdispresethist;位移清零，历史记录清零hisunbalsolveratio1e-5;记录不平衡力，计算直到精度到达1e-5;执行计算solve\step5000\cycle5000;进行开挖;删除水平块体，开挖，不可再充填deleterangex1x2y1y2等价于deleterangeregx1,y1x1,y2x2,y2x2,y1;删除不规那么四边形deleterangeregionx1,y1x2,y2x3,y3x4,y43.5支护局部;加支柱，支架supportx,yanglesegmentwitdthmat;support是指参加支柱，x,y说明支柱所处的位置点。angle说明支柱的轴向与x轴正方向的夹角，以逆时针方向为正。segment指一组支柱的个数，默认值为5。witdth指一组支柱支撑的范围。;mat是指赋予支柱的材料性质属于几号？下面的table指支柱的增阻特性。比方：sup10,20angle=90seg=5width=2mat=6propmat=6sup_kn=-1table103e60.0055e60.01525e60.540e64.液压支架的命令suppxcycwidlsegnmatj;(xc,yc)是指中心点坐标,l指支架的宽度,n指分段数,matj指支柱材料性质为jpropmatjsup_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-610cable(hc2bx,hc2by)(hc2ex,hc2ey)259314e-610;锚杆参数propm=9cb_dens7500cb_ycomp6.3e8cb_yield0.5e7cb_ymod0.5e11;锚固体的参数propm=10cb_kbond=1.6e9cb_sbond=2e6cable(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=7800propmat=16st_yield=6e7st_yresid=6e7st_ycomp=2.5e7propmat=16if_kn=1.35e9if_ks=1.35e9if_tens=0if_fric=18if_coh=0;interface-界面structgenxc=20yc=40.24np=100fang=-40theta=260mat=14thick=0.2propmat=14st_ymod=13.5e9st_prat=0.14st_den=7800propmat=14st_yield=6e7st_yresid=6e7st_ycomp=2.5e7propmat=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在屏幕上将显示建立的块体单元，并在按<Esc>键后退出后，或直接<Enter>继续执行命令。相当于pau+plblnum；删除块体delbl+地址号plotblockstress;显示块体的应力plotblockxdis;显示块体的位移(有x和y方向)plotblockpl;显示块体的塑性区(plastic);屏幕观看下沉与应力等值线,int是指intervalplblydispint0.1plblsyyint1e6udec颜色关键词颜色关键词颜色black黑lgray浅灰blue蓝lblue浅蓝green绿lgreen浅绿cyan蓝绿lcyan浅蓝绿red红lred浅红magenta品红lmagenta浅品红orange橙lorange浅橙brown褐yellow黄dgray灰white白图形输出;抓取图像，在autocad内翻开、编辑;输出应力分布plblpline1syyredsetpldxf25setoutf:\syy100.dxfplpenblpline1syyred；=copyf:\syy100.dxf1.把图形保存下来，能在AUTOCAD中plotblockcableredsuppyestrubl;显示块体、锚杆〔红色〕、支架〔黄色〕、梁〔蓝色〕的图形setplot256dxf；设置图形为256色setoutc:\ss.dxfcopyc:\ss.dxf；把这图形以ss.dxf文件保存在C磁盘下注！！！：开始为dxf格式，移动后另存改为dwg格式，可以减少很大的空间！动画制作;翻开图象抓拍功能，设置图片间隔时间以及窗口大小movieonmoviestep1000windowx1x2y1y2;movieoff【注：只显示改区间内的内容,分别为x1，x2和y1，y2的坐标！window627.5-113plblsyymoviefilekaicaisize800600onstep1000delran407012.015step10000movieoff；电影将被保存至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,0rangejregion1;jsetA,0a,0b,0c,0x0,y0rangejregion1jset90,02.5,02.5,03,01.5,2.5rangejregion1可画成列的线段(见图2)jset0,026,00,02.5,00,0rangejregion1pa;当程度运行至此时暂停,可以看看你所建的局部模型,用continue继续运行下面局部savet.save;建模保存在t.save中岩层赋属性及原岩力平衡计算rest.save;调用已建好的模型genquad12;定义块体最大变形,假设没有此语句,刚所有块体均为刚性块体zonemodelmorange0,1800,45.2;在范围(x1,y1)(x1,y2)(x2,y2)(x2,y1)的块体符合库仑准那么changejcons=2range0,1800,45.2;节理面间接触-coulmb滑移;下面主要是讲岩层赋属性changemat=1rangereg0,00,3180,3180,0;岩层1的范围changemat=2rangereg0,30,4.5180,4.5180,3changemat=3rangereg0,4.50,10.2180,10.2180,4.5changemat=4rangereg0,10.20,11.2180,11.218010.2propmat=1dens=2500;岩层1的密度propmat=2dens=2500propmat=3dens=1300zonek=2e10,g=1.1e10,fric=30,coh=2e6,ten=4e6rangemat=1;岩层1的块体力学参数(见表1)zonek=8.7e9,g=4.2e9,fric=25,coh=1.5e6,ten=1.5e6rangemat=2zonek=13.05e9,g=6.3e9,fric=31,coh=4e6,ten=2e6rangemat=3changejmat=1rangereg0,00,3180,3180,0;在这范围1内的块体间相互接触的接触面力学参数changejmat=2rangereg0,30,4.5180,4.5180,3changejmat=3rangereg0,4.50,10.2180,10.2180,4.5propjmat=1jkn=7e9,jks=2e9,jcoh=0.1e6,jfric=0,jten=0.1e6;在这范围1内的块体间相互接触的接触面力学参数propjmat=2jkn=6e9,jks=2e9,jcoh=1e6,jfric=20,jten=1e6propjmat=3jkn=6e9,jks=2e9,jcoh=1e6,jfric=20,jten=1e6;工作面埋深550米setgravity0,-10;地下岩层主要受重力,还有构造应力(水平应力=垂直应力*侧压系数)insitustr-1.5144e70-1.262e7szz=-1.5144e7ygrad3e402.5e4zgrad03e4;initiasituatestress初始设置应力,侧压系数为1.2;(550-48.2)*2500*(-10)=-1.262e7-1.262e7*1.2=-1.5144e7,1*2500*10=2.5e4,2.5e4*1.2=3.0e4boundstress0,0,-1.375e7range018045.145.3;boundarystresssxx0,sxy0,syy0rangex1x2y1y2;550*2500*(-10)=-1.375e7boundxvel=0range-0.10.1045.2;固定左边界(xvel-也就是x方向的速度),见图3boundxvel=0range179.9180.1045.2;固定右边界boundyvel=0range0180-0.10.1;固定下边界solve;计算savefyuanyan1.saveμ—侧压系数grad表示梯度,即在此方向的变化量(τ=c+σ.tgφ,c值,反映岩石剪切时的粘结阻力,故称岩石的内聚力(或粘结力),φ,值反映岩石剪切时摩擦阻力的大小,故称岩石的内摩擦角,tgφ,相当于摩擦系数f。可得用c、φ两个数在应力圆中判断某种应力状态下的岩石是否产生破坏，通常，岩石愈坚硬，c、φ两个值超大，反之亦然。1G=103M，表1块体力学参数岩层体积模量

K/Gpa剪切模量

G/Gpa体积力

d/N.m-3摩擦角

f/(o)粘结力

C/Mpa抗拉强度

t/Mpa上覆岩层2011250003024根本顶40202700040410直接顶2011250003024煤层13.056.3139003142直接底8.710250003771.5根本底2011250003024接触面力学参数岩层法向刚度

jkn/Gpa切向刚度

jks/Gpa粘结力

jc/Mpa摩擦角

jf/(o)抗拉强度

jt/Mpa上覆岩层721301根本顶20150400直接顶72000煤层621201直接底621201根本底720.100.1;打入锚杆后的锚固体，其摩擦角和粘结力以用弹性模量都会增大，所以锚固体处的属性要改变，否那么锚杆的效果显不著。1全软建模;建模setlogonnewtitle济宁三号特厚煤层放顶煤条件下上行开采机理研究数值模拟round0.02setedg0.1【与genedg有什么区别】setedgn是设置模型中最小块体的边界，要求大于或等于round设置的参数的两倍，具体参见程序自带的commandreference中的setedg命令〔P1—205〕。Genedgedmax是对模型中任意块体产生域,edmax是产生的域中的三角块体的最大边界长度。具体参考commandreference中的genedge命令〔P1-73〕，其中还有一个类似的genquad命令，这个命令只能对模型中长方形块体产生于。setovtol0.55;模型250×110bl00065.525065.52500;按岩性自下而上;第10层,8.0m粉细砂岩互层cr065.525065.5;第9层,2.0m泥岩cr057.525057.5;第8层,1.5m3上煤层cr055.525055.5;第5层,24.0m中细砂岩互层-老顶cr054.025054.0;第4层,12.0m粉砂岩-下直接顶cr030.025030.0;第34层,6.0m3下煤层【为什么加这一层】我当时模拟的就是这层煤的开挖，其中的34并不是真正的第34层，而是我随意设置的，cr018.025018.0;第3层,6.0m3下煤层cr017.025017.0;第2层,2.0m泥岩,直接底cr012.025012.0;第1层,10.0m粉砂岩,老底cr010.025010.0;切眼cr5012.05017.0;右边界cr20012.020017.0;埋设测线setpline027.025027.0133setpline049.7525049.75133【参数133的意义】133〔n〕是一个具体的参数，是设置的测线上的测点数目。具体参见commandreference中的setpline命令。;节理区定义;第8层,8.0m粉细砂岩互层jregid8057.5065.525065.525057.5;第7层,2.0m泥岩jregid7055.5057.525057.525055.5;第6层,1.5m3上煤层jregid6054.0055.525055.525054.0;第5层,24m中细砂岩互层-老顶jregid5030.0054.025054.025030.0;第4层,12.0m粉砂岩-下直接顶jregid4018.0030.025030.025018.0;第34层,6.0m3下煤层jregid34017.0018.025018.025017.0;第3层,6.0m3下煤层jregid3012.0017.025017.025012.0jregid31012.0017.05017.05012.0jregid3220012.020017.025017.025012.0jregid335012.05017.020017.020012.0;第2层,2.0m泥岩,直接底jregid2010.0012.025012.025010.0;第1层,10.0m粉砂岩,老底jregid100010.025010.02500;划节理;第10层,8.0m粉细砂岩互层jset0025000040057.5ranjreg8jset900404050057.5ranjreg8jset9004040502.561.5ranjreg8;第9层,2.0m泥岩jset002500003.00055.5ranjreg7jset9003.003.003.00055.5ranjreg7;第8层,1.5m3上煤层jset002500001.50054.0ranjreg6jset9001.501.502.00054.0ranjreg6;第5层,24m中细砂岩互层-老顶jset0025000020030.0ranjreg5jset900200020030.0ranjreg5;第4层,12.0m粉砂岩-下直接顶jset002500002.00018.0ranjreg4jset9002.00002.00018.0ranjreg4;第34层,6.0m3下煤层遗留jset002500001.00017.0ranjreg34jset9001.001.001.00017.0ranjreg34;第3层下,6.0m3下煤层-采jset00500005.00012.0ranjreg31jset00500005.0020012.0ranjreg32jset001500005.005012.0ranjreg33jset9005.0000505012.0ranjreg33;jset9001.001.001005513.0ranjreg33;第2层,2.0m泥岩,直接底jset002500002.00010.55ranjreg2;第1层,10.0m粉砂岩,老底jset002500005.0000ranjreg1genedg5;保存建模文件sav1.savcalfb.txtqnewres1.savchangecons=1changejcons=2setjcondf=2【命令的含义】这个命令是对新产生的节理设置参数，一般可以不用。具体可以参见commandreference中的setjcondfn这个命令，一般可以不用。;定义岩石材料标号;煤changemat=1ranjreg3changemat=2ranjreg6;粉砂岩changemat=3ranjreg1changemat=4ranjreg4;泥岩changemat=5ranjreg2changemat=6ranjreg7;粉细砂岩互层changemat=7ranjreg8;中细砂岩互层changemat=8ranjreg5;遗留煤层changemat=9ranjreg34;定义节理材料标号;煤changejmat=1ranjreg3changejmat=2ranjreg6;粉砂岩changejmat=3ranjreg1changejmat=4ranjreg4;泥岩changejmat=5ranjreg2changejmat=6ranjreg7;粉细砂互层changejmat=7ranjreg8;中细砂岩互层changejmat=8ranjreg5;遗留煤层changejmat=9ranjreg34;定义岩石材料力学特性;煤propmat=1d=1.4e3k=4e9g=3.5e9propmat=2d=1.4e3k=4e9g=3.5e9;粉砂岩promat=3d=2.60e3k=15e9g=15e9promat=4d=2.60e3k=15e9g=15e9;泥岩promat=5d=2.50e3k=6e9g=9e9promat=6d=2.50e3k=6e9g=9e9;粉细砂岩互层promat=7d=2.65e3k=12e9g=12e9;中细砂岩互层promat=8d=2.75e3k=15e9g=15e9;遗留煤层promat=9d=1.50e3k=5e9g=4e9;定义节理材料力学特性;煤propjmat=1jkn=2e9,jks=0.5e9,jcoh=0.2e6,jfric=3,jtens=0.1e6propjmat=2jkn=3e9,jks=0.5e9,jcoh=0.2e6,jfric=6,jtens=0.1e6;粉砂岩projmat=3jkn=6e9,jks=1.5e9,jcoh=0.2e6,jfric=15,jtens=0.3e6projmat=4jkn=6e9,jks=1.5e9,jcoh=0.2e6,jfric=15,jtens=0.3e6;泥岩projmat=5jkn=7e9,jks=0.8e9,jcoh=0.2e6,jfric=8,jtens=0.2e6projmat=6jkn=7e9,jks=0.8e9,jcoh=0.2e6,jfric=8,jtens=0.2e6;粉细砂岩互层projmat=7jkn=9e9,jks=1.2e9,jcoh=0.3e6,jfric=12,jtens=0.4e6;中细砂岩互层projmat=8jkn=9e9,jks=1.5e9,jcoh=0.3e6,jfric=15,jtens=0.4e6;遗留煤层projmat=9jkn=3e9,jks=0.5e9,jcoh=0,jfric=3,jtens=0setgrav0-10;边界条件boundstress0,0,-12e6range025065.465.6boundxvel=0range-.1,.10,65.5boundxvel=0range249.9,250.10,65.5boundyvel=0range0,250-.1,.1应力应该施加多大？;原岩应力insitustress-6e60.0-12e6szz=-6e6【为什么施加Z向应力】首先你要知道你所用的二维UDEC，所建立的模型的力学模型是平面应变模型而不是平面应力模型，平面应变模型在z轴方向是没有位移的，因此在z轴方向要施加应力。;初次平衡histsolve_rattype1【什么样的命令】这个命令是记录UDEC的运行模式，每次当你运行UDEC时后，光标处不是有一行不断变化的参数吗？这个命令就是指的这个模式，solverat是指在模型运行的时候模型到达平衡的时候的不平衡力的极限比率。具体参见commandreference中的solveratio命令〔P1-219〕。UDEC运算默认的都是这个模式类型，你用这个就行了，大家都习惯用这个的。histvmax;第一测线histxdis40.549.75histxdis80.549.75histxdis100.549.75histxdis120.549.75histxdis160.549.75histxdis180.549.75histydis40.549.75histydis80.549.75histydis100.549.75histydis120.549.75histydis160.549.75histydis180.549.75histxvel40.549.75histxvel80.549.75histxvel100.549.75histxvel120.549.75histxvel160.549.75histxvel180.549.75histyvel40.549.75histyvel80.549.75histyvel100.549.75histyvel120.549.75histyvel160.549.75histyvel180.549.75histsxx40.549.75histsxx80.549.75histsxx100.549.75histsxx120.549.75histsxx160.549.75histsxx180.549.75histsyy40.549.75histsyy80.549.75histsyy100.549.75histsyy120.549.75histsyy160.549.75histsyy180.549.75histnstr40.549.75【是法向应力吗？意义】在这个点40.549.75〔xy〕附近接触的法向应力。具体可以参见commandreference中的histnstrxy命令〔P1-81〕histnstr80.549.75histnstr100.549.75histnstr120.549.75histnstr160.549.75histnstr180.549.75histsstr40.549.75histsstr80.549.75histsstr100.549.75histsstr120.549.75histsstr160.549.75histsstr180.549.75histsdis40.549.75histsdis80.549.75histsdis100.549.75histsdis120.549.75histsdis160.549.75histsdis180.549.75histsvel40.549.75histsvel80.549.75histsvel100.549.75histsvel120.549.75histsvel160.549.75histsvel180.549.75histsxy40.549.75histsxy80.549.75histsxy100.549.75histsxy120.549.75histsxy160.549.75histsxy180.549.75histndis40.549.75histndis80.549.75histndis100.549.75histndis120.549.75histndis160.549.75histnvel40.549.75histnvel80.549.75histnvel100.549.75histnvel120.549.75histnvel160.549.75histnvel180.549.75;第二测线histxdis40.527histxdis80.527histxdis100.527histxdis120.527histxdis160.527histxdis180.527histydis40.527histydis80.527histydis100.527histydis120.527histydis160.527histydis180.527histxvel40.527histxvel80.527histxvel100.527histxvel120.527histxvel160.527histxvel180.527histyvel40.527histyvel80.527histyvel100.527histyvel120.527histyvel160.527histyvel180.527histsxx40.527histsxx80.527histsxx100.527histsxx120.527histsxx160.527histsxx180.527histsyy40.527histsyy80.527histsyy100.527histsyy120.527histsyy160.527histsyy180.527histnstr40.527histnstr80.527histnstr100.527histnstr120.527histnstr160.527histnstr180.527histsstr40.527histsstr80.527histsstr100.527histsstr120.527histsstr160.527histsstr180.527histsdis40.527histsdis80.527histsdis100.527histsdis120.527histsdis160.527histsdis180.527histsvel40.527histsvel80.527histsvel100.527histsvel120.527histsvel160.527histsvel180.527histsxy40.527histsxy80.527histsxy100.527histsxy120.527histsxy160.527histsxy180.527histndis40.527histndis80.527histndis100.527histndis120.527histndis160.527histndis180.527histnvel40.527histnvel80.527histnvel100.527histnvel120.527histnvel160.527histnvel180.527;应力函数轨迹histxhishistyhissolv;保存属性文件sav21.savcalff0.txtqnewres21.savdelran507512.017.0step20000save01.savcalff1.txtqnewres01.savdelran7510012.017.0step30000sav02.savcalff2.txtqnewres02.savdelran10012512.017.0step50000sav03.savcalff3.txtqnewres03.savdelran12515012.017.0step50000sav04.savdelran15017512.017.0step50000sav05.savdelran17520012.017.0solvesav06.savq红涛：这是你模拟的全软的几个命令，局部地方没有看懂，帮我解释一下！另外，我想模拟沿推进方向的可以吗?像这样的效果，来自论文：浅埋煤层采动裂缝损害机理及控制方法研究可以吗？现在黄师兄让我做兴隆庄煤矿的模拟！你这个说的不太明白，由于二维UDEC模拟的模型限制，在模型选择方面，一般只能模拟沿煤层走向方向的，其实我们模拟的就象在工作面中部沿煤层走向〔长壁开采时就是沿推进方向〕进行平剖图，就像你下面的两个图所反映的一样。二维UDEC也可以模拟沿倾斜方向〔也就是工作面方向〕的模型，但是开挖是一次完成的。这种研究一般是研究回采巷道两侧的煤柱问题，或者端头放煤问题。你下面两个图所表示的意思就是沿走向，即推进方向模拟的，一般模拟都是这样的。评：授之以鱼，不如授之以渔！2全硬建模;建模setlogonnewtitle济宁三号特厚煤层放顶煤条件下上行开采机理研究数值模拟round0.02setedg0.1setovtol0.35;模型250×110bl00065.525065.52500;按岩性自下而上;第7层,8.0m粉细砂岩互层cr065.525065.5;第6层,2.0m泥岩cr057.525057.5;第5层,1.5m3上煤层cr055.525055.5;第4层,36.0m中细砂岩互层-老顶cr054.025054.0;第34层,6.0m3下煤层cr018.025018.0;第3层,6.0m3下煤层cr017.025017.0;第2层,2.0m泥岩,直接底cr012.025012.0;第1层,10.0m粉砂岩,老底cr010.025010.0;切眼cr5012.05017.0;右边界cr20012.020017.0;埋设测线setpline032.025032.0133setpline054.7525054.75133;节理区定义;第7层,8.0m粉细砂岩互层jregid7057.5065.525065.525057.5;第6层,2.0m泥岩jregid6055.5057.525057.525055.5;第5层,1.5m3上煤层jregid5054.0055.525055.525054.0;第4层,36m中细砂岩互层-老顶cr70187054cr94189454cr1181811854cr1421814254cr1661816654cr1901819054jregid4018.0054.025054.025018.0jregid41018.0054.07054.07018.0jregid4219018.019054.025054.025018.0jregid437018.07054.019054.019018.0jregid4317018.07054.09454.09418.0jregid4329418.09454.011854.011818.0jregid43311818.011854.014254.014218.0jregid43414218.014254.016654.016618.0jregid43516618.016654.019054.019018.0;第34层,6.0m3下煤层jregid34017.0018.025018.025017.0;第3层,6.0m3下煤层jregid3012.0017.025017.025012.0jregid31012.0017.05017.05012.0jregid3220012.020017.025017.025012.0jregid335012.05017.020017.020012.0;第2层,2.0m泥岩,直接底jregid2010.0012.025012.025010.0;第1层,10.0m粉砂岩,老底jregid100010.025010.02500;划节理;第7层,8.0m粉细砂岩互层jset0025000040057.5ranjreg7jset900404050057.5ranjreg7jset9004040502.561.5ranjreg7;第6层,3.0m泥岩jset002500003.00055.5ranjreg6jset9003.003.003.00055.5ranjreg6;第5层,1.5m3上煤层jset002500001.50054.0ranjreg5jset9001.501.502.00054.0ranjreg5;第4层,36m中细砂岩互层-老顶jset007000040018.0ranjreg41jset9004000140018.0ranjreg41jset00700004019018.0ranjreg42jset900400012019018.0ranjreg42jset0024000407018.0ranjreg431jset90040001207018.0ranjreg431jset0024000409418.0ranjreg432jset90040001209418.0ranjreg432jset00240004011818.0ranjreg433jset900400012011818.0ranjreg433jset00240004014218.0ranjreg434jset900400012014221.0ranjreg434jset00240004016618.0ranjreg435jset900400012016618.0ranjreg435;第34层,6.0m3下煤层遗留jset002500001.00017.0ranjreg34jset9001.001.001.00017.0ranjreg34;第3层下,6.0m3下煤层-采jset00500005.00012.0ranjreg31jset00500005.0020012.0ranjreg32jset001500005.005012.0ranjreg33jset9005.0000505012.0ranjreg33;jset9001.001.001005513.0ranjreg33;第2层,2.0m泥岩,直接底jset002500002.00010.55ranjreg2;第1层,10.0m粉砂岩,老底jset002500005.0000ranjreg1genedg5;保存建模文件sav1.savcalfb.txtqnewres1.savchangecons=1changejcons=2setjcondf=2;定义岩石材料标号;煤changemat=1ranjreg3changemat=2ranjreg5;粉砂岩changemat=3ranjreg1;泥岩changemat=4ranjreg2changemat=5ranjreg6;粉细砂岩互层changemat=6ranjreg7;中细砂岩互层changemat=7ranjreg41changemat=7ranjreg42changemat=31ranjreg431changemat=32ranjreg432changemat=33ranjreg433changemat=34ranjreg434changemat=35ranjreg435;遗留煤层changemat=8ranjreg34;定义节理材料标号;煤changejmat=1ranjreg3changejmat=2ranjreg5;粉砂岩changejmat=3ranjreg1;泥岩changejmat=4ranjreg2changejmat=5ranjreg6;粉细砂互层changejmat=6ranjreg7;中细砂岩互层changejmat=7ranjreg41changejmat=7ranjreg42changejmat=31ranjreg431changejmat=32ranjreg432changejmat=33ranjreg433changejmat=34ranjreg434changejmat=35ranjreg435;遗留煤层changejmat=8ranjreg34;定义岩石材料力学特性;煤propmat=1d=1.4e3k=4e9g=3.5e9propmat=2d=1.4e3k=4e9g=3.5e9;粉砂岩promat=3d=2.60e3k=10e9g=10e9;泥岩promat=4d=2.50e3k=6e9g=9e9promat=5d=2.50e3k=6e9g=9e9;粉细砂岩互层promat=6d=2.65e3k=12e9g=12e9;中细砂岩互层promat=7d=2.75e3k=18e9g=18e9promat=31d=2.75e3k=18e9g=18e9promat=32d=2.75e3k=18e9g=18e9promat=33d=2.75e3k=18e9g=18e9promat=34d=2.75e3k=18e9g=18e9promat=35d=2.75e3k=18e9g=18e9;遗留煤层promat=8d=1.50e3k=5e9g=4e9;定义节理材料力学特性;煤propjmat=1jkn=2e9,jks=0.5e9,jcoh=0.2e6,jfric=3,jtens=0.1e6propjmat=2jkn=3e9,jks=0.5e9,jcoh=0.2e6,jfric=6,jtens=0.1e6;粉砂岩projmat=3jkn=6e9,jks=1.0e9,jcoh=0.2e6,jfric=10,jtens=0.3e6;泥岩projmat=4jkn=7e9,jks=0.8e9,jcoh=0.2e6,jfric=8,jtens=0.2e6projmat=5jkn=7e9,jks=0.8e9,jcoh=0.2e6,jfric=8,jtens=0.2e6;粉细砂岩互层projmat=6jkn=9e9,jks=1.2e9,jcoh=0.3e6,jfric=12,jtens=0.4e6;中细砂岩互层projmat=7jkn=9e9,jks=5e9,jcoh=0.3e6,jfric=20,jtens=0.4e6projmat=31jkn=9e9,jks=3e9,jcoh=0.3e6,jfric=20,jtens=0.4e6projmat=32jkn=9e9,jks=3e9,jcoh=0.3e6,jfric=20,jtens=0.4e6projmat=33jkn=9e9,jks=3e9,jcoh=0.3e6,jfric=20,jtens=0.4e6projmat=34jkn=9e9,jks=3e9,jcoh=0.3e6,jfric=20,jtens=0.4e6projmat=35jkn=9e9,jks=3e9,jcoh=0.3e6,jfric=20,jtens=0.4e6;遗留煤层projmat=8jkn=3e9,jks=0.5e9,jcoh=0,jfric=3,jtens=0setgrav0-10;边界条件boundstress0,0,-12e6range025065.465.6boundxvel=0range-.1,.10,65.5boundxvel=0range249.9,250.10,65.5boundyvel=0range0,250-.1,.1;原岩应力insitustress-6e60.0-12e6szz=-6e6;初次平衡histsolve_rattype1histvmax;第一测线histxdis40.549.75histxdis80.549.75histxdis100.549.75histxdis120.549.75histxdis160.549.75histxdis180.549.75histydis40.549.75histydis80.549.75histydis100.549.75histydis120.549.75histydis160.549.75histydis180.549.75histxvel40.549.75histxvel80.549.75histxvel100.549.75histxvel120.549.75histxvel160.549.75histxvel180.549.75histyvel40.549.75histyvel80.549.75histyvel100.549.75histyvel120