udec命令總結(jié)精華-正宗

《Udec命令總結(jié)》精華版1.安裝、翻開、保存(1)安裝：①執(zhí)行Udec3.1→將Crack文件中的內(nèi)容替換；②復(fù)制Udec.exe,粘貼為快捷方式→屬性，目標(biāo)，參加空格256【開始內(nèi)存8M，將內(nèi)存改為256】，可以復(fù)制快捷方式至桌面或硬盤。(2)翻開：①直接在udec>命令行輸入；②寫好程序，udec>命令輸入call，然后將*.txt文件拖入命令行，執(zhí)行(3)保存：輸入saved:\kaicai.sav,調(diào)用命令rest，將kaicai.sav拖入命令行【！??！文件名最好不用漢字，有時(shí)候不識(shí)別】注釋：如果保存為saved:\111\kaicai.sav,注意其中111文件必須提前建好，否那么無法保存或者保存錯(cuò)誤(4)操作：Ctrl+Z選中圖像可以放大，Ctrl+Z雙擊復(fù)原，屏幕中會(huì)出現(xiàn)十字叉，按住鼠標(biāo)左鍵不放，移動(dòng)光標(biāo)直到你滿意的窗口為止；pause暫停，此時(shí)可以觀察任何信息；continue繼續(xù)調(diào)用下面程序段。ESC可以隨時(shí)進(jìn)行停止，但不能繼續(xù)；英文分號(hào);表示注釋不運(yùn)行命令。2.根本命令2.1根本設(shè)置Udec>n【new刷新窗口，從新調(diào)用一個(gè)程序，修改后的*.txt文件必須輸入n，重新運(yùn)行文本文件】Udec>title【或heading代表標(biāo)題，后面緊跟標(biāo)題的名稱。如：hangdaomoni或濟(jì)寧三號(hào)井圍巖變形破壞規(guī)律研究】Udec>roundd【“圓角〞命令，Udec中所有的塊體都有圓角，目的是為防止塊體懸掛在有棱角的節(jié)點(diǎn)上，由于塊體懸掛將產(chǎn)生應(yīng)力集中。d指塊體與塊體之間的圓角半徑，默認(rèn)值是0.5，其值要求小于模型中最小塊體的最短那條邊長(zhǎng)的二分之一，最大圓角長(zhǎng)度不能超過塊體平均棱長(zhǎng)的1%。在block命令前指定圓角長(zhǎng)度。如：round0.05】圓角圓角圖110×10塊體圓角Udec>setovtol0.5【此命令是指層與層之間的嵌入厚度，當(dāng)提示為“overlaptoolarge〞時(shí)就需要修改此值更大一些，可以顯示設(shè)計(jì)的塊體，plotoverlap！嵌入太大的原因可能為塊體強(qiáng)度太小】Udec>setlogonprintsetlogoff【命令用于導(dǎo)出數(shù)據(jù)，可以將數(shù)據(jù)導(dǎo)出至硬盤f:\】命令：setlogf:\yuanyan.logprintpline1syy【記錄測(cè)線1應(yīng)力】printpline2ydissetlogoff;設(shè)置觀測(cè)線setpline1x1,y1x2,y2n(n--觀測(cè)線分的段數(shù))setpline25,20300,2010;定義測(cè)線，起始點(diǎn)坐標(biāo)，測(cè)線分段2.2圖形劃分;塊體命令Udec>blockx1,y1x2,y2x3,y3x4,y4【建立模型框架，四個(gè)坐標(biāo)角點(diǎn)必須按“順時(shí)針〞排列，也可以為五個(gè)坐標(biāo)點(diǎn)?！繄D25個(gè)坐標(biāo)點(diǎn)plbl顯示塊體；plholdbl=plbl+pause屏幕上直接顯示劃分的塊體，關(guān)閉顯示后將繼續(xù)運(yùn)行；delrangex1x2y1y2刪除塊體；fixrangex1x2y1y2固定塊體不動(dòng);畫圓弧Udec>arcxc,ycxb,ybthetancracks【xc,yc圓弧中心，xb,yb指畫圓弧的起始點(diǎn)，theta沿逆時(shí)針方向圓弧的角度;ncracks指構(gòu)成指圓弧的折線段數(shù)，數(shù)目越多，圓弧越圓；①圓弧命令必須配合Crack命令，但tunnel命令可以在輸入后直接顯示為圓形；②先劃分巷道，然后劃分其它巖層節(jié)理】命令：;開挖巷道節(jié)理，形成巷道空間,巷道必須先開挖出來arc52.1785.5953.9785.591808cra50.3783.5950.3785.59cra50.3783.5953.9783.59cra53.9783.5953.9785.59;建圓形隧道、圓形巷道Udec>tunnelxc,ycrn【xc,yc圓的中心，r半徑，n構(gòu)成圓的線段數(shù)。如：tunnel4,5245①先輸tunnel，然后輸crack，裂隙不貫穿隧道;②隧道開挖：方法一：刪除編號(hào)plblnum→deleteblock1920方法二：deleterange-11-11方法三：指定窗口windowsx1x2y1y2→delrangewindow，不一定能全部刪除，可以繼續(xù)刪除編號(hào)方法四：定義區(qū)域，刪除區(qū)域，jregnx1,y1x2,y2x3,y3x4,y4→deletejregn】;節(jié)理劃分、定義Udec>crackx1,y1x2,y2【劃分節(jié)理后，最好輸入pause暫停，plbl顯示劃分的塊體是否正確或文本中輸入plholdbl】Udec>jregionidnx1,y1x2,y2x3,y3x4,y4【節(jié)理區(qū)域定義，賦予節(jié)理區(qū)域地址號(hào)，可以刪除任何區(qū)域deletejregn】命令：jset90,04,04,06,00,-50rangejregn【在定義的節(jié)理區(qū)域內(nèi)進(jìn)行節(jié)理的再劃分;jsetα、t、g、s(xo,yo)n:傾角，節(jié)理長(zhǎng)度,巖橋長(zhǎng)度,垂距，起始點(diǎn)坐標(biāo),n區(qū)域標(biāo)號(hào)】ααgts圖3節(jié)理組定義；不規(guī)那么塊體劃分Udec>jregid414.65073.016078.1616019.960;定義節(jié)理區(qū)域4Udec>changemat4range0100-65-50;定義Udec>voredg5iter20rou0.1ranjreg4【edge平均邊長(zhǎng)，最少20倍rounding長(zhǎng)度；iterations迭代次數(shù)，默認(rèn)5；round，產(chǎn)生塊體的圓角長(zhǎng)度，默認(rèn)0.5；rtol容許偏差，默認(rèn)0.02*round；seedn】2.3參數(shù)定義genquad10rangexlxuylyu;在指定的區(qū)域生成一定寬度的單元zonemodelmorangexlxuylyu;使指定的區(qū)域材料采用摩爾--庫侖本構(gòu)關(guān)系計(jì)算〔即彈塑性〕changejcons=2rangexlxuylyu;使指定的區(qū)域節(jié)理遵循摩爾--庫侖準(zhǔn)那么計(jì)算〔即彈塑性〕changemat=1rangexlxuylyuchangemat=2rangexlxuylyuchangemat=3rangexlxuylyu;指定各巖層的材料標(biāo)號(hào)changejmat=1rangexlxuylyuchangejmat=2rangexlxuylyuchangejmat=3rangexlxuylyu;指定各巖層的節(jié)理標(biāo)號(hào)propmat=1dens=2000propmat=2dens=2650propmat=3dens=2700;指定各材料的密度，比方1號(hào)材料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為材料的內(nèi)摩擦角，cohesion為材料的內(nèi)聚力，tension為材料的抗拉強(qiáng)度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為節(jié)理的法向剛度，jks為節(jié)理的切向剛度，jfriction為節(jié)理的內(nèi)摩擦角，jcohesion為節(jié)理的內(nèi)聚力，jtension為節(jié)理的抗拉強(qiáng)度塊體力學(xué)參數(shù)k--體積模量；g--剪切模量；fric--摩擦角；coh--粘結(jié)力〔內(nèi)聚力〕；ten--抗拉強(qiáng)度；d--體積力接觸面力學(xué)參數(shù)jkn--法向剛度；jks--切向強(qiáng)度；jc--粘結(jié)力；jf--摩擦角；jt--抗拉強(qiáng)度jregid10,00,15300,15300,0delete;定義節(jié)理，命名1，定義范圍，四點(diǎn)坐標(biāo)，delete為常規(guī)語言jset90,05,05,05,00,0rangjreg1;定義節(jié)理角度90°；節(jié)理的長(zhǎng)度；節(jié)理的距離，即縱向間隔；橫向距離，即橫向隔5一個(gè)；起始點(diǎn)坐標(biāo)為0,0jset90,05,05,05,02.5,5rangjreg1;定義節(jié)理的另外一項(xiàng)表1巖石塊體參數(shù)巖石名稱彈性模量E泊松比體積模量K剪切模量G砂巖19.30.3826.87.0粉砂巖26.30.2215.610.8石灰?guī)r28.50.2922.611.1頁巖、泥巖11.10.298.84.3大理巖55.80.2537.222.3花崗巖73.80.2243.930.2巖石名稱內(nèi)摩擦角粘聚力抗拉強(qiáng)度砂巖27.827.21.17粉砂巖32.134.7\頁巖、泥巖14.438.4\石英巖48.070.6\石灰?guī)r42.06.721.58花崗巖51.055.1\玄武巖31.066.213.1體積模量K和剪切模量G之間的關(guān)系：genquad10range03000100;當(dāng)03000100范圍內(nèi)的塊體的對(duì)角長(zhǎng)度小于10時(shí)就是可變形體。zonemodelmorange03000100;定義摩爾庫倫模型的范圍，X的范圍，Y的范圍changejcons=5range03000100;定義節(jié)理的某個(gè)屬性changemat=1rangereg0,00,15300,15300,0;定義物質(zhì)1的范圍，即賦予這個(gè)范圍為物質(zhì)1,只賦名，沒有定義屬性changemat=2rangereg0,150,20300,20300,15changemat=3rangereg0,200,23300,23300,20changemat=4rangereg0,230,28300,28300,23changejmat=1rangereg0,00,15300,15300,0;定義節(jié)理1的范圍setjmatdf2;定義節(jié)理屬性，使符合摩爾庫倫準(zhǔn)那么普適。2代表摩爾庫倫準(zhǔn)那么;參數(shù)設(shè)置propmat=1dens=2720;物質(zhì)的密度zonek=12.12e9g=10.26e9fric=33coh=5.77e6ten=3.6e6rangemat=1;定義物質(zhì)1的體積模量、剪切模量、內(nèi)摩擦角、粘聚力、抗拉強(qiáng)度propjmat=1jkn=11e8jks=12e7jfric=12jcoh=0jten=0;定義節(jié)理的屬性法向剛度、切向剛度、內(nèi)摩擦角、粘聚力、抗拉強(qiáng)度1)D--質(zhì)量密度F--摩擦角B--體積模量C--內(nèi)聚力(粘結(jié)力)S--剪切模量2)應(yīng)力--正號(hào)代表張力，負(fù)號(hào)代表壓力；應(yīng)變--正的應(yīng)變代表伸長(zhǎng)，負(fù)的應(yīng)變代表壓縮重力--正號(hào)的重力物質(zhì)往下拉，負(fù)號(hào)的重力將物質(zhì)往上提?。。∶芏缺仨氃O(shè)置densityforblockmaterical1iszerocannotcycle----塊體密度未設(shè)2.4初始應(yīng)力setgravity0,-9.81;設(shè)置重力加速度，x方向?yàn)?，y方向?yàn)?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;定義邊界條件，三向應(yīng)力為0，上邊界，范圍，縱向波動(dòng)范圍000〔x、xy、y〕bounstress00-8.6095397e6range;定義邊界應(yīng)力條件，豎直方向的應(yīng)力為，負(fù)號(hào)表示方向向下insitustr-7.5e60-15e6szz-7.5e6ygrad1.1018e402.203607e4;定義應(yīng)力屬性，sxx方向即水平應(yīng)力，sxy方向即剪切應(yīng)力，syy方向即豎直應(yīng)力，szz方向的應(yīng)力，ygrad應(yīng)力梯度boundxvel=0range-0.10.10100;X位移邊界,位移波動(dòng)范圍，范圍boundxvel=0range299.9300.10100;X位移邊界，同上boundyvel=0range-0.1300.1-0.1,0.1;Y位移邊界histsolve_ratiotype1;設(shè)置不平衡力的精度，普適histxdisp10,20;測(cè)該點(diǎn)的X方向位移histsxx10,20;測(cè)該點(diǎn)X方向的應(yīng)力histsyy10,20;測(cè)該點(diǎn)Y方向的應(yīng)力resetdispresethist;位移清零，歷史記錄清零hisunbalsolveratio1e-5;記錄不平衡力，計(jì)算直到精度到達(dá)1e-5;執(zhí)行計(jì)算solve\step5000\cycle5000;進(jìn)行開挖;刪除水平塊體，開挖，不可再充填deleterangex1x2y1y2等價(jià)于deleterangeregx1,y1x1,y2x2,y2x2,y1;刪除不規(guī)那么四邊形deleterangeregionx1,y1x2,y2x3,y3x4,y43.5支護(hù)局部;加支柱，支架supportx,yanglesegmentwitdthmat;support是指參加支柱，x,y說明支柱所處的位置點(diǎn)。angle說明支柱的軸向與x軸正方向的夾角，以逆時(shí)針方向?yàn)檎egment指一組支柱的個(gè)數(shù)，默認(rèn)值為5。witdth指一組支柱支撐的范圍。;mat是指賦予支柱的材料性質(zhì)屬于幾號(hào)？下面的table指支柱的增阻特性。比方：sup10,20angle=90seg=5width=2mat=6propmat=6sup_kn=-1table103e60.0055e60.01525e60.540e64.液壓支架的命令suppxcycwidlsegnmatj;(xc,yc)是指中心點(diǎn)坐標(biāo),l指支架的寬度,n指分段數(shù),matj指支柱材料性質(zhì)為jpropmatjsup_kn-1;sup_kn指支架的剛度,-1與表1相對(duì)應(yīng)delrangex1x2y1y2;挖掘范圍(x1,y1)(x1,y2)(x2,y2)(x2,y1)的塊體,由支架支撐table1042.e60.055.0e60.16.0e6;表1表示的是液壓支架的(P-DS)特性曲線，后面為四個(gè)點(diǎn)的坐標(biāo)。目前，液壓支柱的最大支撐力為12*106N顯示支架：plblsupportred;加錨桿cablex1,y1x2,y2npointmatsasteelmatg等;x1,y1x2,y2指錨桿的兩個(gè)端點(diǎn);npoint指錨桿的粘結(jié)點(diǎn)數(shù);mats指錨桿材料性質(zhì),材料編號(hào)，按著前面給定的材料號(hào)繼續(xù)往下編；asteel指錨桿橫截面積；matg指錨固體性質(zhì)；cable(hc1bx,hc1by)(hc1ex,hc1ey)259314e-610cable(hc2bx,hc2by)(hc2ex,hc2ey)259314e-610;錨桿參數(shù)propm=9cb_dens7500cb_ycomp6.3e8cb_yield0.5e7cb_ymod0.5e11;錨固體的參數(shù)propm=10cb_kbond=1.6e9cb_sbond=2e6cable(73.4,2.91)(71.15,4.99)1012380e-611;錨桿，起始點(diǎn)坐標(biāo)，后面四個(gè)數(shù)值可固定propmat12cb_dens7800cb_ycomp430e6cb_yield160e3cb_ymod1.3e11;錨固體的性質(zhì)propmat11cb_kbond6.3e9cb_sbond6e5;定義錨桿的屬性，密度、抗拉強(qiáng)度、屈服強(qiáng)度；巷道斷面加梯子梁的命令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,...觀測(cè)點(diǎn)位移plblnu(mber)顯示塊體標(biāo)號(hào)plblcabredstruredcab--錨桿stru--梁plblsxx/syy/szz顯示應(yīng)力圖plbldis位移矢量圖plboundpl顯示塑性區(qū)plbounddis顯示位移矢量圖plholdblnum在屏幕上將顯示建立的塊體單元，并在按<Esc>鍵后退出后，或直接<Enter>繼續(xù)執(zhí)行命令。相當(dāng)于pau+plblnum；刪除塊體delbl+地址號(hào)plotblockstress;顯示塊體的應(yīng)力plotblockxdis;顯示塊體的位移(有x和y方向)plotblockpl;顯示塊體的塑性區(qū)(plastic);屏幕觀看下沉與應(yīng)力等值線,int是指intervalplblydispint0.1plblsyyint1e6udec顏色關(guān)鍵詞顏色關(guān)鍵詞顏色black黑lgray淺灰blue藍(lán)lblue淺藍(lán)green綠lgreen淺綠cyan藍(lán)綠lcyan淺藍(lán)綠red紅lred淺紅magenta品紅lmagenta淺品紅orange橙lorange淺橙brown褐yellow黃dgray灰white白圖形輸出;抓取圖像，在autocad內(nèi)翻開、編輯;輸出應(yīng)力分布plblpline1syyredsetpldxf25setoutf:\syy100.dxfplpenblpline1syyred；=copyf:\syy100.dxf1.把圖形保存下來，能在AUTOCAD中plotblockcableredsuppyestrubl;顯示塊體、錨桿〔紅色〕、支架〔黃色〕、梁〔藍(lán)色〕的圖形setplot256dxf；設(shè)置圖形為256色setoutc:\ss.dxfcopyc:\ss.dxf；把這圖形以ss.dxf文件保存在C磁盤下注?。。。洪_始為dxf格式，移動(dòng)后另存改為dwg格式，可以減少很大的空間！動(dòng)畫制作;翻開圖象抓拍功能，設(shè)置圖片間隔時(shí)間以及窗口大小movieonmoviestep1000windowx1x2y1y2;movieoff【注：只顯示改區(qū)間內(nèi)的內(nèi)容,分別為x1，x2和y1，y2的坐標(biāo)！window627.5-113plblsyymoviefilekaicaisize800600onstep1000delran407012.015step10000movieoff；電影將被保存至C:\ProgramFiles\Itasca\udec310下面，翻開方式C:\ProgramFiles\Itasca\shared\Utility\movie，必須點(diǎn)擊start才開始播放;【注：1)：電影的制作過程實(shí)際上為抓頻的過程，因此，在此期間最好將其它的程序先關(guān)掉！2):在進(jìn)行電影之前，輸入plblsyy,那么整個(gè)過程將顯示應(yīng)力的變化過程，如果輸入plblplinesyyred那么整個(gè)過程將顯示測(cè)線應(yīng)力變化的過程。輸入先輸入plblsyy后輸入plblplinesyyred，那么將在電影中同時(shí)出現(xiàn)！輸入plbound，那么將出現(xiàn)裂隙帶的范圍！4實(shí)例分析: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)兩點(diǎn)間畫直線jregionid10,00,1050,1050,0delete;刪除此區(qū)域的塊體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;當(dāng)程度運(yùn)行至此時(shí)暫停,可以看看你所建的局部模型,用continue繼續(xù)運(yùn)行下面局部savet.save;建模保存在t.save中巖層賦屬性及原巖力平衡計(jì)算rest.save;調(diào)用已建好的模型genquad12;定義塊體最大變形,假設(shè)沒有此語句,剛所有塊體均為剛性塊體zonemodelmorange0,1800,45.2;在范圍(x1,y1)(x1,y2)(x2,y2)(x2,y1)的塊體符合庫侖準(zhǔn)那么changejcons=2range0,1800,45.2;節(jié)理面間接觸-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的塊體力學(xué)參數(shù)(見表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內(nèi)的塊體間相互接觸的接觸面力學(xué)參數(shù)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內(nèi)的塊體間相互接觸的接觸面力學(xué)參數(shù)propjmat=2jkn=6e9,jks=2e9,jcoh=1e6,jfric=20,jten=1e6propjmat=3jkn=6e9,jks=2e9,jcoh=1e6,jfric=20,jten=1e6;工作面埋深550米setgravity0,-10;地下巖層主要受重力,還有構(gòu)造應(yīng)力(水平應(yīng)力=垂直應(yīng)力*側(cè)壓系數(shù))insitustr-1.5144e70-1.262e7szz=-1.5144e7ygrad3e402.5e4zgrad03e4;initiasituatestress初始設(shè)置應(yīng)力,側(cè)壓系數(shù)為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;計(jì)算savefyuanyan1.saveμ—側(cè)壓系數(shù)grad表示梯度,即在此方向的變化量(τ=c+σ.tgφ,c值,反映巖石剪切時(shí)的粘結(jié)阻力,故稱巖石的內(nèi)聚力(或粘結(jié)力),φ,值反映巖石剪切時(shí)摩擦阻力的大小,故稱巖石的內(nèi)摩擦角,tgφ,相當(dāng)于摩擦系數(shù)f。可得用c、φ兩個(gè)數(shù)在應(yīng)力圓中判斷某種應(yīng)力狀態(tài)下的巖石是否產(chǎn)生破壞，通常，巖石愈堅(jiān)硬，c、φ兩個(gè)值超大，反之亦然。1G=103M，表1塊體力學(xué)參數(shù)巖層體積模量

K/Gpa剪切模量

G/Gpa體積力

d/N.m-3摩擦角

f/(o)粘結(jié)力

C/Mpa抗拉強(qiáng)度

t/Mpa上覆巖層2011250003024根本頂40202700040410直接頂2011250003024煤層13.056.3139003142直接底8.710250003771.5根本底2011250003024接觸面力學(xué)參數(shù)巖層法向剛度

jkn/Gpa切向剛度

jks/Gpa粘結(jié)力

jc/Mpa摩擦角

jf/(o)抗拉強(qiáng)度

jt/Mpa上覆巖層721301根本頂20150400直接頂72000煤層621201直接底621201根本底720.100.1;打入錨桿后的錨固體，其摩擦角和粘結(jié)力以用彈性模量都會(huì)增大，所以錨固體處的屬性要改變，否那么錨桿的效果顯不著。1全軟建模;建模setlogonnewtitle濟(jì)寧三號(hào)特厚煤層放頂煤條件下上行開采機(jī)理研究數(shù)值模擬round0.02setedg0.1【與genedg有什么區(qū)別】setedgn是設(shè)置模型中最小塊體的邊界，要求大于或等于round設(shè)置的參數(shù)的兩倍，具體參見程序自帶的commandreference中的setedg命令〔P1—205〕。Genedgedmax是對(duì)模型中任意塊體產(chǎn)生域,edmax是產(chǎn)生的域中的三角塊體的最大邊界長(zhǎng)度。具體參考commandreference中的genedge命令〔P1-73〕，其中還有一個(gè)類似的genquad命令，這個(gè)命令只能對(duì)模型中長(zhǎng)方形塊體產(chǎn)生于。setovtol0.55;模型250×110bl00065.525065.52500;按巖性自下而上;第10層,8.0m粉細(xì)砂巖互層cr065.525065.5;第9層,2.0m泥巖cr057.525057.5;第8層,1.5m3上煤層cr055.525055.5;第5層,24.0m中細(xì)砂巖互層-老頂cr054.025054.0;第4層,12.0m粉砂巖-下直接頂cr030.025030.0;第34層,6.0m3下煤層【為什么加這一層】我當(dāng)時(shí)模擬的就是這層煤的開挖，其中的34并不是真正的第34層，而是我隨意設(shè)置的，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;埋設(shè)測(cè)線setpline027.025027.0133setpline049.7525049.75133【參數(shù)133的意義】133〔n〕是一個(gè)具體的參數(shù)，是設(shè)置的測(cè)線上的測(cè)點(diǎn)數(shù)目。具體參見commandreference中的setpline命令。;節(jié)理區(qū)定義;第8層,8.0m粉細(xì)砂巖互層jregid8057.5065.525065.525057.5;第7層,2.0m泥巖jregid7055.5057.525057.525055.5;第6層,1.5m3上煤層jregid6054.0055.525055.525054.0;第5層,24m中細(xì)砂巖互層-老頂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;劃節(jié)理;第10層,8.0m粉細(xì)砂巖互層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中細(xì)砂巖互層-老頂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【命令的含義】這個(gè)命令是對(duì)新產(chǎn)生的節(jié)理設(shè)置參數(shù)，一般可以不用。具體可以參見commandreference中的setjcondfn這個(gè)命令，一般可以不用。;定義巖石材料標(biāo)號(hào);煤changemat=1ranjreg3changemat=2ranjreg6;粉砂巖changemat=3ranjreg1changemat=4ranjreg4;泥巖changemat=5ranjreg2changemat=6ranjreg7;粉細(xì)砂巖互層changemat=7ranjreg8;中細(xì)砂巖互層changemat=8ranjreg5;遺留煤層changemat=9ranjreg34;定義節(jié)理材料標(biāo)號(hào);煤changejmat=1ranjreg3changejmat=2ranjreg6;粉砂巖changejmat=3ranjreg1changejmat=4ranjreg4;泥巖changejmat=5ranjreg2changejmat=6ranjreg7;粉細(xì)砂互層changejmat=7ranjreg8;中細(xì)砂巖互層changejmat=8ranjreg5;遺留煤層changejmat=9ranjreg34;定義巖石材料力學(xué)特性;煤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;粉細(xì)砂巖互層promat=7d=2.65e3k=12e9g=12e9;中細(xì)砂巖互層promat=8d=2.75e3k=15e9g=15e9;遺留煤層promat=9d=1.50e3k=5e9g=4e9;定義節(jié)理材料力學(xué)特性;煤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;粉細(xì)砂巖互層projmat=7jkn=9e9,jks=1.2e9,jcoh=0.3e6,jfric=12,jtens=0.4e6;中細(xì)砂巖互層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應(yīng)力應(yīng)該施加多大？;原巖應(yīng)力insitustress-6e60.0-12e6szz=-6e6【為什么施加Z向應(yīng)力】首先你要知道你所用的二維UDEC，所建立的模型的力學(xué)模型是平面應(yīng)變模型而不是平面應(yīng)力模型，平面應(yīng)變模型在z軸方向是沒有位移的，因此在z軸方向要施加應(yīng)力。;初次平衡histsolve_rattype1【什么樣的命令】這個(gè)命令是記錄UDEC的運(yùn)行模式，每次當(dāng)你運(yùn)行UDEC時(shí)后，光標(biāo)處不是有一行不斷變化的參數(shù)嗎？這個(gè)命令就是指的這個(gè)模式，solverat是指在模型運(yùn)行的時(shí)候模型到達(dá)平衡的時(shí)候的不平衡力的極限比率。具體參見commandreference中的solveratio命令〔P1-219〕。UDEC運(yùn)算默認(rèn)的都是這個(gè)模式類型，你用這個(gè)就行了，大家都習(xí)慣用這個(gè)的。histvmax;第一測(cè)線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【是法向應(yīng)力嗎？意義】在這個(gè)點(diǎn)40.549.75〔xy〕附近接觸的法向應(yīng)力。具體可以參見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;第二測(cè)線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;應(yīng)力函數(shù)軌跡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紅濤：這是你模擬的全軟的幾個(gè)命令，局部地方?jīng)]有看懂，幫我解釋一下！另外，我想模擬沿推進(jìn)方向的可以嗎?像這樣的效果，來自論文：淺埋煤層采動(dòng)裂縫損害機(jī)理及控制方法研究可以嗎？現(xiàn)在黃師兄讓我做興隆莊煤礦的模擬！你這個(gè)說的不太明白，由于二維UDEC模擬的模型限制，在模型選擇方面，一般只能模擬沿煤層走向方向的，其實(shí)我們模擬的就象在工作面中部沿煤層走向〔長(zhǎng)壁開采時(shí)就是沿推進(jìn)方向〕進(jìn)行平剖圖，就像你下面的兩個(gè)圖所反映的一樣。二維UDEC也可以模擬沿傾斜方向〔也就是工作面方向〕的模型，但是開挖是一次完成的。這種研究一般是研究回采巷道兩側(cè)的煤柱問題，或者端頭放煤?jiǎn)栴}。你下面兩個(gè)圖所表示的意思就是沿走向，即推進(jìn)方向模擬的，一般模擬都是這樣的。評(píng)：授之以魚，不如授之以漁！2全硬建模;建模setlogonnewtitle濟(jì)寧三號(hào)特厚煤層放頂煤條件下上行開采機(jī)理研究數(shù)值模擬round0.02setedg0.1setovtol0.35;模型250×110bl00065.525065.52500;按巖性自下而上;第7層,8.0m粉細(xì)砂巖互層cr065.525065.5;第6層,2.0m泥巖cr057.525057.5;第5層,1.5m3上煤層cr055.525055.5;第4層,36.0m中細(xì)砂巖互層-老頂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;埋設(shè)測(cè)線setpline032.025032.0133setpline054.7525054.75133;節(jié)理區(qū)定義;第7層,8.0m粉細(xì)砂巖互層jregid7057.5065.525065.525057.5;第6層,2.0m泥巖jregid6055.5057.525057.525055.5;第5層,1.5m3上煤層jregid5054.0055.525055.525054.0;第4層,36m中細(xì)砂巖互層-老頂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;劃節(jié)理;第7層,8.0m粉細(xì)砂巖互層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中細(xì)砂巖互層-老頂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;定義巖石材料標(biāo)號(hào);煤changemat=1ranjreg3changemat=2ranjreg5;粉砂巖changemat=3ranjreg1;泥巖changemat=4ranjreg2changemat=5ranjreg6;粉細(xì)砂巖互層changemat=6ranjreg7;中細(xì)砂巖互層changemat=7ranjreg41changemat=7ranjreg42changemat=31ranjreg431changemat=32ranjreg432changemat=33ranjreg433changemat=34ranjreg434changemat=35ranjreg435;遺留煤層changemat=8ranjreg34;定義節(jié)理材料標(biāo)號(hào);煤changejmat=1ranjreg3changejmat=2ranjreg5;粉砂巖changejmat=3ranjreg1;泥巖changejmat=4ranjreg2changejmat=5ranjreg6;粉細(xì)砂互層changejmat=6ranjreg7;中細(xì)砂巖互層changejmat=7ranjreg41changejmat=7ranjreg42changejmat=31ranjreg431changejmat=32ranjreg432changejmat=33ranjreg433changejmat=34ranjreg434changejmat=35ranjreg435;遺留煤層changejmat=8ranjreg34;定義巖石材料力學(xué)特性;煤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;粉細(xì)砂巖互層promat=6d=2.65e3k=12e9g=12e9;中細(xì)砂巖互層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;定義節(jié)理材料力學(xué)特性;煤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;粉細(xì)砂巖互層projmat=6jkn=9e9,jks=1.2e9,jcoh=0.3e6,jfric=12,jtens=0.4e6;中細(xì)砂巖互層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;原巖應(yīng)力insitustress-6e60.0-12e6szz=-6e6;初次平衡histsolve_rattype1histvmax;第一測(cè)線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