UDEC经典学习总结 张科学修改 完成

1.把图形保存下来，能在AUTOCAD(图片)中打开、编辑plot block cable red supp ye stru bl;显示块体、锚杆（红色）、支架（黄色）、梁（蓝色）的图形set plot dxf(jpg) 256；设置图形为256色set out c:\ss.dxf(jpg) (set out 1.dxf 这个文件保存在当前带数值模拟的文件夹里)copy c:\ss.dxf(jpg) ( copy 1.dxf 这个文件保存在当前带数值模拟的文件夹里)；把这图形以ss.dxf文件保存在C磁盘下2.把数据导出set log onprint pline 1 ydis;把pline 1的y方向位移的数据导出到UDEC目录下udec.txt文件中set log off3.plot block stress ;显示块体的应力plot block dis ;显示块体的位移(有x和y方向)plot block pl ;显示块体的塑性区(plastic)save xx.save ;保存计算结果res xx.save ;调用4.液压支架的命令supp xc yc wid l seg n mat j;(xc,yc)是指中心点坐标,l指支架的宽度,n指分段数,mat j指支柱材料性质为j prop mat j sup_kn -1 ;sup_kn指支架的刚度,-1与表1相对应del range x1 x2 y1 y2 ;挖掘范围(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)的块体,由支架支撑table 1 0 42.e6 0.05 5.0e6 0.1 6.0e6;表1表示的是液压支架的(P-DS)特性曲线5 巷道施工中断面加梯子梁的命令(见图1)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.35e9 if_ks=1.35e9 if_tens=0 if_fric=18 if_coh=0;interface-界面例子讲解:建模(以米为单位,;后为解释部分)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 range jregion 1 ;jset A,0 a,0 b,0 c,0 x0,y0 range jregion 1jset 90,0 2.5,0 2.5,0 3,0 1.5,2.5 range jregion 1 可画成列的线段(见图2)jset 0,0 26,0 0,0 2.5,0 0,0 range jregion 1pa ;当程度运行至此时暂停,可以看看你所建的部分模型,用continue继续运行下面部分save t.save ;建模保存在t.save中岩层赋属性及原岩力平衡计算res t.save ;调用已建好的模型gen quad 12 ;定义块体最大变形,若没有此语句,刚所有块体均为刚性块体zone model mo range 0,180 0,45.2 ;在范围(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)的块体符合库仑准则change jcons=2 range 0,180 0,45.2 ;节理面间接触-coulmb滑移;下面主要是讲岩层赋属性change （范围） mat=1 range reg 0,0 0,3 180,3 180,0 ;岩层1的范围 change mat=2 range reg 0,3 0,4.5 180,4.5 180,3change mat=3 range reg 0,4.5 0,10.2 180,10.2 180,4.5 change mat=4 range reg 0,10.2 0,11.2 180,11.2 180 10.2 prop （赋值的意思） mat=1 dens=2500 ;岩层1的密度prop mat=2 dens=2500 prop mat=3 dens=1300zone （块体） k=2e10, g=1.1e10, fric=30,coh=2e6,ten=4e6 range mat=1 ;岩层1的块体力学参数(见表1) zone k=8.7e9, g=4.2e9, fric=25,coh=1.5e6,ten=1.5e6 range mat=2 zone k=13.05e9, g=6.3e9, fric=31,coh=4e6,ten=2e6 range mat=3change jmat(节理)=1 range reg 0,0 0,3 180,3 180,0 ;在这范围1内的块体间相互接触的接触面力学参数 change jmat=2 range reg 0,3 0,4.5 180,4.5 180,3change jmat=3 range reg 0,4.5 0,10.2 180,10.2 180,4.5prop jmat=1 jkn=7e9, jks=2e9, jcoh=0.1e6 , jfric=0,jten=0.1e6 ;在这范围1内的块体间相互接触的接触面力学参数 prop jmat=2 jkn=6e9, jks=2e9, jcoh=1e6 , jfric=20,jten=1e6 prop jmat=3 jkn=6e9, jks=2e9, jcoh=1e6 , jfric=20,jten=1e6;工作面埋深550米set gravity 0,-10 ;地下岩层主要受重力,还有构造应力(水平应力=垂直应力*侧压系数)bound stress 0,0,-1.375e7 range 0 180 45.1 45.3 ; boundary stress sxx0 ,sxy0, syy0 range x1 x2 y1 y2;(550—煤层到上边界的距离)*2500*(-10)= -1.375e7 insitu str -1.5144e7 0 -1.262e7 szz=-1.5144e7 ygrad 3e4 0 2.5e4 zgrad 0 3e4;initi asituate stress 初始设置应力 ,侧压系数为1.2;syy=(550+煤层到下边界的距离)*2500*(-10)= -1.262e7（在下边界上施加的力）-1.262e7*1.2=-1.5144e7 ,y 方向的梯度为=1*2500*10=2.5e4, x 、z 方向的梯度为=2.5e4*1.2=3.0e4bound xvel=0 range -0.1 0.1 0 45.2 ;固定左边界(xvel-也就是x 方向的速度),见图3 bound xvel=0 range 179.9 180.1 0 45.2 ;固定右边界bound yvel=0 range 0 180 -0.1 0.1 ;固定下边界（因为上边界有上覆岩层故不需固定了） solve ;计算save fyuanyan1.save注: insitu str sxx sxy syy szz=-1.6638e8 ygrad sxxy sxyy syyy zgrad szzx szzy (xgrad sxxx sxyx syyy,上面没有写,表示在x 方向没有变化 )(σx=sxx ,σy=syy, σz=szz,τxy=sxy ,τxz=sxz, τyz=syz 弹性力学;grad 表示梯度,即在此方向的变化量) (xgrad sxxx sxyx syyy ygrad sxxy sxyy syyy zgrad szzx szzy(sxx=sxx0+sxxx.x+sxxy .y,syy=syy0+syyx.x+syyy.y ,szz=szz0+szzx.x+szzy.y) (μ—侧压系数,sxx=syy . μ,szz=syy . μ,sxxy=syyy . μ,szzy=syyy. μ)(τ=c+σ.tg φ,c 值,反映岩石剪切时的粘结阻力,故称岩石的内聚力(或粘结力) ,φ,值反映岩石剪切时摩擦阻力, 的大小,故称岩石的内摩擦角,tg φ,相当于摩擦系数f 。

可得用c 、φ两个数在应力圆中判断某种应力状态下的岩石是否产生破坏，通常，岩石愈坚硬，c 、φ两个值超大，反之亦然。

1G=103M=109()21E G υ=+，()312EK υ=-】锚杆的效果显不著。

round=0.01set ovtol=3bl 0,0 0,45.2 180,45.2 180,0;基本底3m(1.5*1)crack 0,3 180,3crack 57,3 57,0crack 73.8,3 73.8,0jregion id 11 0,0 0,3 57,3 57,0 deletejset 0,0 57,0 0,0 1,0 0,0 range jregion 11jset 90,0 1,0 1,0 1.5,0 0,0 range jregion 11jset 90,0 1,0 1,0 1.5,0 0.75,1 range jregion 11jregion id 12 57,0 57,3 73.8,3 73.8,0 deletejset 0,0 16.8,0 0,0 0.25,0 57,0 range jregion 12jset 90,0 0.25,0 0.25,0 0.5,0 57,0 range jregion 12jset 90,0 0.25,0 0.25,0 0.5,0 57.25,0.25 range jregion 12 jregion id 13 73.8,0 73.8,3 180,3 180,0 deletejset 0,0 106.2,0 0,0 1,0 73.8,0 range jregion 13jset 90,0 1,0 1,0 1.5,0 73.8,0 range jregion 13jset 90,0 1,0 1,0 1.5,0 74.55,1 range jregion 13;直接底1.5m(1*0.5/0.5*0.25)crack 0,4.5 180,4.5crack 57,3 57,4.5crack 73.8,3 73.8,4.5;前57米(1*0.5)jregion id 21 0,3 0,4.5 57,4.5 57,3 deletejset 0,0 57,0 0,0 0.5,0 0,3 range jregion 21jset 90,0 0.5,0 0.5,0 1,0 0.5,3 range jregion 21jset 90,0 0.5,0 0.5,0 1,0 0,3.5 range jregion 21;57到73.8米(0.5*0.25)jregion id 22 57,3 57,4.5 73.8,4.5 73.8,3 deletejset 0,0 16.8,0 0,0 0.25,0 57,3 range jregion 22jset 90,0 0.25,0 0.25,0 0.5,0 57,3 range jregion 22jset 90,0 0.25,0 0.25,0 0.5,0 57.25,3.25 range jregion 22 ;后73.8到180米(1*0.5)jregion id 23 73.8,3 73.8,4.5 180,4.5 180,3 deletejset 0,0 106.2,0 0,0 0.5,0 73.8,3 range jregion 23jset 90,0 0.5,0 0.5,0 1,0 73.8,3 range jregion 23jset 90,0 0.5,0 0.5,0 1,0 74.3,3.5 range jregion 23;5.7米的煤层，采3.5米crack 0,10.2 180,10.2crack 0,9.5 57,9.5crack 73.8,9.5 180,9.5crack 57,4.5 57,10.2crack 73.8,4.5 73.8,10.2crack 63,4.5 63,7.75crack 67.8,4.5 67.8,10.2crack 57,7.75 67.8,7.75crack 67.8,8 180,8;前57米(1*0.5)jregion id 31 0,4.5 0,9.5 57,9.5 57,4.5 deletejset 0,0 57,0 0,0 0.5,0 0,4.5 range jregion 31jset 90,0 0.5,0 0.5,0 1,0 0,4.5 range jregion 31jset 90,0 0.5,0 0.5,0 1,0 0.5,5 range jregion 31jregion id 30 0,9.5 0,10.2 57,10.2 57,9.5 deletejset 90,0 0.7,0 0.7,0 1,0 0,9.5 range jregion 30;57到73.8米(0.5*0.25)jregion id 321 57,4.5 57,7.75 63,7.75 63,4.5 deletejset 0,0 6,0 0,0 0.25,0 57,4.5 range jregion 321jset 90,0 0.25,0 0.25,0 0.5,0 57,4.5 range jregion 321 jset 90,0 0.25,0 0.25,0 0.5,0 57.25,4.75 range jregion 321 jregion id 322 63,4.5 63,7.75 67.8,7.75 67.8,4.5 delete jset 0,0 3.2,0 0,0 0.25,0 63,4.5 range jregion 322jset 90,0 0.25,0 0.25,0 0.5,0 63,4.5 range jregion 322 jset 90,0 0.25,0 0.25,0 0.5,0 63.25,4.75 range jregion 322 jregion id 323 57,7.75 57,10.2 67.8,10.2 67.8,7.75 delete jset 0,0 10.8,0 0,0 0.25,0 57,7.75 range jregion 323jset 90,0 0.25,0 0.25,0 0.5,0 57.25,7.75 range jregion 323 jset 90,0 0.25,0 0.25,0 0.5,0 57,8 range jregion 323 jregion id 324 67.8,4.5 67.8,8 73.8,8 73.8,4.5 deletejset 0,0 6,0 0,0 0.25,0 67.8,4.5 range jregion 324jset 90,0 0.25,0 0.25,0 0.5,0 68,4.5 range jregion 324jset 90,0 0.25,0 0.25,0 0.5,0 68.25,4.75 range jregion 324 jregion id 325 67.8,8 67.8,10.2 73.8,10.2 73.8,8 delete jset 0,0 6,0 0,0 0.25,0 67.8,4.5 range jregion 325jset 90,0 0.25,0 0.25,0 0.5,0 68,8 range jregion 325jset 90,0 0.25,0 0.25,0 0.5,0 68.25,8.25 range jregion 325 ;后73.8到180米(1*0.5)jregion id 331 73.8,4.5 73.8,9.5 180,9.5 180,4.5 delete jset 0,0 106.2,0 0,0 0.5,0 73.8,4.5 range jregion 331jset 90,0 0.5,0 0.5,0 1,0 74.3,4.5 range jregion 331jset 90,0 0.5,0 0.5,0 1,0 73.8,5 range jregion 331jregion id 300 73.8,9.5 73.8,10.2 180,10.2 180,9.5 delete jset 90,0 0.7,0 0.7,0 1,0 74.3,9.5 range jregion 300;直接顶1m(1.5*1)crack 0,11.2 180,11.2crack 57,10.2 57,11.2crack 73.8,10.2 73.8,11.2;前57米(1*0.5)jregion id 41 0,10.2 0,11.2 57,11.2 57,10.2 delete;jset 0,0 57,0 0,0 0.5,0 0,10.2 range jregion 41jset 90,0 1,0 1,0 1.5,0 0.75,10.2 range jregion 41;jset 90,0 0.5,0 0.5,0 1,0 0,10.7 range jregion 41;57到73.8米(0.5*0.25)jregion id 42 57,10.2 57,11.2 73.8,11.2 73.8,10.2 delete jset 0,0 16.8,0 0,0 0.25,0 57,10.2 range jregion 42jset 90,0 0.25,0 0.25,0 0.5,0 57.25,10.2 range jregion 42 jset 90,0 0.25,0 0.25,0 0.5,0 57,10.45 range jregion 42;后73.8到180米(1.5*1)jregion id 43 73.8,10.2 73.8,11.2 180,11.2 180,10.2 delete ;jset 0,0 106.2,0 0,0 0.5,0 73.8,10.2 range jregion 43;jset 90,0 0.5,0 0.5,0 1,0 73.8,10.2 range jregion 43jset 90,0 1,0 1,0 1.5,0 73.8,10.2 range jregion 43;基本顶10m(12*5)crack 0,21.2 180,21.2jregion id 5 0,11.2 0,21.2 180,21.2 180,11.2 deletejset 0,0 180,0 0,0 5,0 0,11.2 range jregion 5jset 90,0 5,0 5,0 12,0 0,11.2 range jregion 5jset 90,0 5,0 5,0 12,0 6,16.2 range jregion 5;上覆岩层24m(6*4)crack 0,45.2 180,45.2jregion id 6 0,21.2 0,45.2 180,45.2 180,21.2 deletejset 0,0 180,0 0,0 4,0 0,21.2 range jregion 6jset 90,0 4,0 4,0 5,0 3,21.2 range jregion 6jset 90,0 4,0 4,0 5,0 0,25.2 range jregion 6save inir.saveres inir.savegen quad 1.6 1.2 range 0,180 0,11.2gen quad 13 5.5 range 0,180 11.2,45.2zone model mo range 0,180 0,45.2 ;( change cons=3 range 0,180 0,45.2)change jcons=2 range 0,180 0,45.2;设定不同岩层、煤层change mat=1 range reg 0,0 0,3 180,3 180,0change mat=2 range reg 0,3 0,4.5 180,4.5 180,3change mat=3 range reg 0,4.5 0,10.2 180,10.2 180,4.5change mat=4 range reg 0,10.2 0,11.2 180,11.2 180 10.2change mat=5 range reg 0,11.2 0,21.2 180,21.2 180,11.2change mat=6 range reg 0,21.2 0,45.2 180,45.2 180,21.2change jmat=1 range reg 0,0 0,3 180,3 180,0change jmat=2 range reg 0,3 0,4.5 180,4.5 180,3change jmat=3 range reg 0,4.5 0,10.2 180,10.2 180,4.5change jmat=4 range reg 0,10.2 0,11.2 180,11.2 180 10.2change jmat=5 range reg 0,11.2 0,21.2 180,21.2 180,11.2change jmat=6 range reg 0,21.2 0,45.2 180,45.2 180,21.2prop mat=1 dens=2500prop mat=2 dens=2500prop mat=3 dens=1390prop mat=4 dens=2500prop mat=5 dens=2700prop mat=6 dens=2500zone k=2e10, g=1.1e10, fric=40, coh=8.5e6, ten=5.6e6 range mat=1 (prop mat=1 k=2e10, g=1.1e10, fric=40, coh=8.5e6, ten=5.6e6 )zone k=17e9, g=8.9e9, fric=34, coh=6e6, ten=2.4e6 range mat=2zone k=13.05e9, g=6.3e9, fric=28, coh=4e6, ten=2.1e6 range mat=3 zone k=2e10, g=1.1e10, fric=33, coh=6.2e6, ten=2.3e6 range mat=4 zone k=2.4e10, g=1.3e10, fric=40, coh=9.5e6, ten=5.4e6 range mat=5 zone k=2e10, g=1.1e10, fric=35, coh=8.4e6, ten=4e6 range mat=6prop jmat=1 jkn=14e9, jks=8.7e9, jcoh=6.5e6 , jfric=35,jten=0prop jmat=2 jkn=9.5e9, jks=6.4e9, jcoh=4.8e6 , jfric=28,jten=0prop jmat=3 jkn=6e9, jks=4.7e9, jcoh=2.8e6 , jfric=20,jten=0prop jmat=4 jkn=8.9e9, jks=5.9e9, jcoh=4.9e6, jfric=27, jten=0prop jmat=5 jkn=16e9, jks=9.5e9, jcoh=7.8e6 , jfric=36,jten=0prop jmat=6 jkn=14e9, jks=8.7e9, jcoh=6.5e6 , jfric=29,jten=0;工作面埋深550米set gravity 0,-10bound stress 0,0,-1.375e7 range 0 180 45.1 45.3insitu str -1.722e7 0 -1.435e7 szz=-1.722e7 range 0 180 21.2 45.2insitu ygrad 3e4 0 2.5e4 zgrad 0 3e4 range 0 180 21.2 45.2insitu str -1.7544e7 0 -1.462e7 szz=-1.7544e7 range 0 180 11.2 21.2insitu ygrad 3.24e4 0 2.7e4 zgrad 0 3.24e4 range 0 180 11.2 21.2insitu str -1.7574e7 0 -1.4645e7 szz=-1.7574e7 range 0 180 10.2 11.2insitu ygrad 3e4 0 2.5e4 zgrad 0 3e4 range 0 180 10.2 11.2insitu str -1.7669076e7 0 -1.472423e7 szz=-1.7669076e7 range 0 180 4.5 10.2 insitu ygrad 1.668e4 0 1.39e4 zgrad 0 1.668e4 range 0 180 4.5 10.2insitu str -1.7804076e7 0 -1.483673e7 szz=-1.7804076e7 range 0 180 0 4.5 insitu str 3e4 0 2.5e4 zgrad 0 3e4 range 0 180 0 4.5bound xvel=0 range -0.1 0.1 0 45.2bound xvel=0 range 179.9 180.1 0 45.2bound yvel=0 range 0 180 -0.1 0.1solvesave yuanyan.savecal kh.txtres yuanyan.savereset dispreset vel;开掘巷道并支护巷道del range reg 63,4.5 63,7.75 67.8,7.75 67.8,4.5;加顶梯子梁(顶梁4.6米)stru gen xc 65.4 yc 4.95 np 100 fa 0 theta 180 mat=16 thick 0.2 connect prop mat=16 st_ymod=13.5e9 st_prat=0.14 st_den=7800prop mat=16 st_yield=60e6 st_yresid=60e6 st_ycomp=35e6prop mat=16 if_kn=1.35e9 if_ks=1.35e9 if_tens=0 if_fric=18 if_coh=0;两帮梯子梁(2.75米);stru gen xc 65.4 yc 6.15 np 100 fa 150.64 theta 59.62 mat=17 thick 0.2;prop mat=17 st_ymod=21e9 st_prat=0.15 st_den=7800;prop mat=17 st_yield=40e6 st_yresid=40e6 st_ycomp=40e6;prop mat=17 if_kn=1e9 if_ks=1e9 if_tens=0 if_fric=50 if_coh=0;stru gen xc 65.4 yc 6.15 np 100 fa 329.75 theta 59.62 mat=18 thick 0.2;prop mat=18 st_ymod=21e9 st_prat=0.15 st_den=7800;prop mat=18 st_yield=40e6 st_yresid=40e6 st_ycomp=40e6;prop mat=18 if_kn=1e9 if_ks=1e9 if_tens=0 if_fric=50 if_coh=0;顶锚杆(L2.4米,6根)cable (63.2,7.75) (62.379,10.005) 5 10 314e-6 11 4e4 connectcable (64.08,7.75) (64.08,10.15) 5 10 314e-6 11 4e4 connectcable (64.96,7.75) (64.96,10.15) 5 10 314e-6 11 4e4 connectcable (65.84,7.75) (65.84,10.15) 5 10 314e-6 11 4e4 connectcable (66.72,7.75) (66.72,10.15) 5 10 314e-6 11 4e4 connectcable (67.6,7.75) (68.421,10.005) 5 10 314e-6 11 4e4 connect;左帮锚杆(L2米,4根)cable (63,4.95) (61,4.95) 5 10 314e-6 11 4e4 connectcable (63,5.8) (61,5.8) 5 10 314e-6 11 4e4 connectcable (63,6.65) (61,6.65) 5 10 314e-6 11 4e4 connectcable (63,7.5) (61.121,8.184) 5 10 314e-6 11 4e4 connect;右帮锚杆(L2米,4根)cable (67.8,4.95) (69.8,4.95) 5 10 314e-6 11 4e4 connectcable (67.8,5.8) (69.8,5.8) 5 10 314e-6 11 4e4 connectcable (67.8,6.65) (69.8,6.65) 5 10 314e-6 11 4e4 connectcable (67.8,7.5) (69.679,8.184) 5 10 314e-6 11 4e4 connect;锚杆参数20*2000prop m=10 cb_dens 7500 cb_ycomp 430e6 cb_yield 260e3 cb_ymod 1.2e11 ;锚固体的参数prop m=11 cb_kbond=6.3e9 cb_sbond=6e6cable (64.5,7.75) (64.5,14.95) 5 12 314e-6 13 4e4cable (65.8,7.75) (65.8,14.95) 5 12 314e-6 13 4e4;锚索参数20*7200prop m=12 cb_dens 7500 cb_ycomp 700e6 cb_yield 500e3 cb_ymod 1.2e11 ;锚固体的参数prop m=13 cb_kbond=6.3e9 cb_sbond=6e6chang mat=7 range region 60,4.5 60,10.2 71,10.2 71,4.5prop mat=7 dens=1390zone k=14.05e9, g=6.9e9,fric=37, coh=4.4e6, ten=2.1e6 range mat=7solvesave kh.savecall kgx.txtres kh.savereset dispreset vel;采3.5米煤层del range reg 67.8,4.5 67.8,8 167,8 167,4.5stru del range 67.7,68 4.5,7.75cable del range 67.8,70 4.5 7.8def aaxx=68.58supp xx 4.6 wid 1.5 seg 3 mat 1 endcommandxx=xx+1.55end_loopendaasupp 65.4 4.8 angle 90 wid 4.2 seg 3 mat 2 prop m 1 sup_kn -1 sup_tmax=1000prop m 2 sup_kn -2table 1 0 4.2e6 0.05 5.0e6 0.1 6.0e6table 2 0 2.5e6 0.2 3.2e6set pline 60,11.2 72,11.2 12set pline 63,7.75 67.8,7.75 12set pline 60,11.2 72,11.2 12set pline 63,4.5 63,7.75 8set pline 67.8,4.5 67.8,7.75 8;顶板下沉量hist ydis 63.3,7.75hist ydis 67.5,7.75hist ydis 65.4,7.75;底板鼓起量hist ydis 63.3,4.5hist ydis 67.5,4.5hist ydis 65.4,4.5;左帮水平移近量hist xdis 63,5hist xdis 63,6.1hist xdis 63,7.2;右帮水平移近量hist xdis 67.8,5hist xdis 67.8,6.1hist xdis 67.8,7.2step 60000save kgx.savecall kgs.txtres kgx.save;放顶煤2.2米def aaaxx=73.23supp xx 4.6 delendcommandxx=xx+1.55end_loopendaaadel range reg 157,4.5 157,8 180,8 180,4.5del range reg 72.5,8 72.5,10.2 180,10.2 180,8set pline 60,11.2 72,11.2 12set pline 63,7.75 67.8,7.75 12set pline 60,11.2 72,11.2 12set pline 63,4.5 63,7.75 8set pline 67.8,4.5 67.8,7.75 8;顶板下沉量hist ydis 63.3,7.75hist ydis 67.5,7.75hist ydis 65.4,7.75;底板鼓起量hist ydis 63.3,4.5hist ydis 67.5,4.5hist ydis 65.4,4.5;左帮水平移近量hist xdis 63,5hist xdis 63,6.1hist xdis 63,7.2;右帮水平移近量hist xdis 67.8,5hist xdis 67.8,6.1hist xdis 67.8,7.2set sup_delsolvesave kgs.savecall kgss.txt用的较多的是:1.changeCH ANGE keyword <keyword><ra nge. . . >Block, joint and cable element material characteristics are prescribed and changed with the CHANGE command. All blocks with centroids lying within the optional range (see Section 1.1.3) have block material characteristics changed. Likewise, all joints with contact coordinates lying within the optional range, or cables with nodes lying within the optional range, have material characteristics changed. If no range is specified, all blocks and joints will have characteristics changed according to the keywords given below. Extended zone models anduser-defined models are assigned using the extended ZONE command. User-defined joint models are defined using the JMODEL command. Local storage joint models are assigned using the JOINT command.The following keywords are used to change characteristics.1. Block Characteristicscons nConstitutive number n is assigned to designated deformable blocks(see Table 1.2).m at nMaterial property number n is assigned to designated rigid or deformableblocks. (All blocks initially default to mat = 1. The maximum value for n is 50.)Table 1.2 Constitutive models for deformable blockscons Model Description0 null material (The null model is used to model excavated material.The stresses within the null block are automatically set to zero.)1 linearly elastic, isotropic (default)3 elastic/plastic, Mohr-Coulomb failure (This model should be used withcaution since accurate solution to plasticity problems requires thatthe triangular zoning have a gridpoint at the centroid of each block.GENERATE quad zoning should be used whenever possible to improveplasticity analyses. However, no significant errors have been noted inproblems for which the above criterion has not been met.)6 elastic/plastic, Drucker-Prager failure (The same caution discussed abovefor the Mohr-Coulomb model (cons = 3) also applies for the Drucker-Prager model.)CH ANGENOTES:1. Constitutive models may also be assigned to regions of zones within blocks withthe ZONE command.2. All block constitutive models are described in Section 2 in Theory and Background.2. Joint Characteristicsjc ons nConstitutive number n is assigned to designated contacts (see Table1.3).jm at nMaterial property number n is assigned to designated contacts. (All contacts initially default to jmat = 1. The maximum value for n is 50.)Table 1.3 Joint constitutive modelsjcons Model Description1 point contact elastic/plastic with Coulomb slip failure (units are [force/displacement] for contact stiffnesses, and [force] for cohesion and tension)2 joint area* contact elastic/plastic with Coulomb slip failure (units are [stress] for cohesion and tension, and[stress/displacement] for joint stiffnesses) (default)3 continuously yielding joint model (see Section 3 in Theory and Background for a detailed explanation)5 same as jcons = 2, except that the internal fracture flag is set for each joint segment when joint shear or tensile strength is exceeded. If the fracture flag is set, residual values for friction, cohesion and tension are used in all subsequent calculations.7 optional Barton-Bandis (BB) joint model. See Section 3 in Special Features for details.* NOTE: The minimum joint area is limited to twice the rounding length, so thatit is not necessary to specify point contact properties if both point contacts andarea contacts occur between blocks.3. Cable Characteristicscab le matg <mats>Material property number matg is assigned to designated cable nodes.Material property number mats is assigned to designated cable elements.Midpoint of element must lie within the range to be changed.(All cable elements and nodes default to mat=1. The maximum valuefor matg or mats is 50.)4. Domain Characteristicsdm at nMaterial property number n is assigned to designated domains. (Alldomains initially default to dmat=1. The maximum value for n is 50.)CAB LE x1 y1 x2 y2 npoint mats matg <preten><keyword>delete <range>Execution of this command creates reinforcing elements that explicitly model theshear behavior of a grout annulus. These elements are generated between endpoints(x1,y1) and (x2,y2).NOTE: For excavation problems, if point (x1,y1) is inside the excavation periphery,the first nodal point will be on the excavation periphery. The point (x2,y2) shouldalways be located in the rock mass. This reinforcing logic can only be used withdeformable blocks, and the CABLE command must be invoked after the GENERATEcommand. Reinforcement properties and grout properties are specified and storedusing the PROPERTY command. The following parameters are also required.matg material property number for grout. (The PROPERTY commandshould be used to specify cb kbond and cb sbond for grout.)NOTE: The units for cb kbond and cb sbond are [force / cable length/ displacement] and [force / cable length], respectively.mats material property number for reinforcing (i.e., steel). (The PROPERTYcommand should be used to specify cb ymod, cb yield, cb ycomp,cb fstrain, and cb density for the material.)npoint number of lumped mass nodal points (npoint ≥2)Pre-tensioning of the reinforcing can be specified by providing a value (in unitsof force) for the optional parameter preten. Always use the PRINT property cablecommand to check property assignment.The following keywords may be used.connect The optional connect keyword will shift the cable node end closest to a structural (beam) element node to coincide with the beam node.The structural element nodes must be created before the cable nodes.The connection is not allowed to fail. The common, or shared, nodewill not appear in the list of cable nodes—it does appear in the list of structural element nodes. The cable element list will show the commonnode ID number at one end of the cable element and the normalcable node at the other. Beam and cable nodes are entirely separateentities. This is because being that the beam nodes are in contactwith the surface of a block, whereas the cable nodes are connected tothe model in the interior of a zone. All common nodes are thereforestructural element nodes, rather than cable nodes. As a consequence,it is possible for the cable element information display to indicateidentical node IDs at either end when the connect keyword is used.UDEC Version 4.0delete The optional delete keyword will remove all cable elements from the specified range. If no range is given, all cable elements will bedeleted.extend The extend keyword causes the end of the current cable to connect to the end node of an existing cable element.The cable reinforcing model is described in Section 1 in Special Features.The material properties for segments of an existing cable element can be changed with the CHANGE cable command.Plot ca ble <keyword><n>plots location of cable elements. For most of the keywords (shown below),the identification number n for the cable (defined by the CABLEcommand) may be used to plot only the cable elements associatedwith the number. The corresponding value of n for the cable maybe found by using the command PLOT cable number. The following keywords may be used.af ail cable element axial failureax ial <n> axial forcePL OT ca ble ele mentele ment cable element numbersfa il shows failure modes in color.gf ail cable node grout failureno de cable node numbersnu mber the identification number of the cable group used, for example, for line plotssd isp nodal point displacement vectorssh ear <n> relative shear force between cable node and host material st rain <n> axial strain in cable elementssvel nodal point velocity vectorsxd isp <n> x-displacement of cable nodesxvel <n> x-velocity of cable nodesyd isp <n> y-displacement of cable nodesyvel <n> y-velocity of cable nodesNOTE:(1) The line plot switch can be used to plot a cable variable as a line plot (see plot category 4).(2) The sense of the cable variable plots can be reversed by giving the switch yrev after the cable variable keyword.Plot stru ct <keyword><n>plots location of structural (beam) elements and associated variables. (See the STRUCT command and Section 1 in Special Features.) For certain keywords (shown below), the identification number, n, for the beam (defined by the STRUCT command) may be used to plot only the structural elements associated with that number. The corresponding value of n for the beam may be found by using the command PLOT struct number. The following keywords may be used.af ail structural element axial failureav el <n> angular velocity of elementax ial <n> axial forceif ail interface failurein ormal <n> normal force at structural interfaceint erface structural interface locationis hear <n> shear force at structural interfacemo ment <n> momentnu mber ID number of the element group used, for example, forline plotssd isp nodal point displacement vectorssh ear <n> shear forcesv el velocity vectors of structural nodesth ick structural element is plotted with thickness rather thana single linexd isp <n> x-displacement of structural nodesxv el <n> x-velocity of structural nodesyd isp <n> y-displacement of structural nodesyv el <n> y-velocity of structural nodesNOTE:(1) The line plot switch can be used to plot a structural element variableas a line plot (see plot category 4).(2) The sense of structural element variable plots can be reversed bygiving the switch yrev after the structural element keyword.Plot su pport structural support elementsPrint ca ble information on cable-reinforcing elements. The output is divided into cable-element information and cable-node information. (SeeSection 1 in Special Features for further descriptions.) The columnheadings for cable elements are:(1) cable element address(2) cable element segment identification number(3) cable node identification number — first node(4) cable node identification number — second node(5) material number for cable element(6) cross-sectional area of cable element(7) axial force in cable element (tension is negative)(8) length of cable element(9) axial strain in cable element (tension is negative)(10) failure indicator for element:0 elastic1 at yield in tension2 elastic, tension yield in past3 at yield in compression4 elastic, compression yield in past5 failure in extensional strainThe column headings for cable nodes are:(1) cable node address(2) cable node identification number(3) x-coordinate of cable node(4) y-coordinate of cable node(5) x-velocity of cable node(6) y-velocity of cable node(7) x-displacement of cable node(8) y-displacement of cable node(9) material number for cable node(10) x-force applied to node(11) y-force applied to node(12) shear force in grout(13) node fixed in x if = 1, else 0(14) node fixed in y if = 1, else 0(15) failure indicator for node:0 elastic1 yield2 elastic, yield in past(16) address of zone in which node is currently locatedPlot stru ct keywordinformation on structural elements. The following keywords are available.el ement <keyword>information on structural element segments. The following optional keywords are available to print information separately.d isp displacements. The column headings are:(1) structural element address(2) x-coordinate of midpoint(3) y-coordinate of midpoint(4) x-displacement of element(5) y-displacement of elementf orce forces. The column headings are:(1) structural element ID number(2) axial force (compression positive)。