udec命令讲解

1、; influence of the placement of backfill; in deep long-wall excavation;round = 0.01 ;倒角长度set ovtol = 0.2 ;接触搭接长度的容差，缺省值为倒角长度的一半block 0,0 0,50 80,50 80,0 ;生成块体;命令格式 block ; real joints 真实节理crack 0,2.5 80,2.5crack 0,7.5 80,7.5crack 0,12.5 80,12.5crack 0,17.5 80,17.5; 命令格式 crack （x1,y1）（x2,y2） ; vertic

2、al joints defining the future chambers; will be assigned high properties；开挖分段垂直“节理”定义将要开挖的空间，设置高的参数。crack 10,0 10,2.5crack 20,0 20,2.5crack 30,0 30,2.5crack 40,0 40,2.5crack 50,0 50,2.5crack 60,0 60,2.5crack 70,0 70,2.5; artificial joints for better zoning of the near roof;为了更好的划分顶板附近的区域而设置的“虚拟节理”cr

3、ack 0,3.5 80,3.5crack 0,2 80,2;gen quad 2 1 range 0 80 0 2；命令格式 gen quad xw range（xl，xu）（yl，yu） 块体的变形特性命令gen quad 2 0.5 range 0 80 2 2.5gen quad 2 1 range 0 80 2.5 3.5gen quad 4 2 range 0 80 3.5 7.5gen quad 8 5 range 0 80 7.5 17.5gen quad 16 6.5 range 0 80 17.5 50; assign material numbers to chamber

4、s；对开挖的空间设置材料号change mat=1 range 0 80 0 50 ；range（xl，x2）（yl，y2）change mat 2 range 0 10 0 2.5change mat 3 range 10 20 0 2.5change mat 4 range 20 30 0 2.5change mat 5 range 30 40 0 2.5change mat 6 range 40 50 0 2.5; assign material numbers to real and artificial joints;对实际的节理和虚拟的节理设置材料号change jmat=2 ra

5、nge 0 80 0 2.4 angle=-1, 1change jmat=1 range 0 80 2.4 2.6 angle=-1, 1change jmat=2 range 0 80 2.6 4.9 angle=-1, 1change jmat=1 range 0 80 4.9 50 angle=-1, 1change jmat=2 range 0 80 0 50 angle=89,91set jmatdf 2 ；为新的接触设置默认的材料号2（循环一步后产生）; define material properties 设置材料参数prop mat 1 dens=2700prop jmat

6、1 jkn=1e11 jks=1e10 jcoh=0.05e6 jfric=30 jten=0prop jmat 2 jkn=1e11 jks=1e10 jcoh=1e20 jfric=0 jten=1e20 jkn、 jks切向、法向normal刚度、切向刚度；cohesion粘聚力； internal friction angle 内摩擦角、jfric节理面摩擦角； tensile strength 抗拉强度；prop mat 2 den 1000prop mat 3 den 1000prop mat 4 den 1000prop mat 5 den 1000prop mat 6 den

7、1000; rock mass is strain-softening material岩体是应变软化材料zone model ss range 0 80 0 50zone bulk=38.9e9 shear=29.9e9 fric=45 dil=0 range 0 80 0 50zone coh=20e6 ctable=1 range 0 80 0 50table 1 0,20e6 0.03,0; initial stresses 初始应力insitu stress -33.8e6 0 -67.5e6 szz -50e6 命令insitu在x方向始化所有的应力为-33.8mpa和在y方向初始

8、化应力为-67.5mpa。在平面之外的z方向也给予初值zz=-50mpa。对于弹性块体分析，z方向的应力未初始化并不影响平面应变问题的解。然而，对于塑性分析，z方向的应力可能影响到破坏状态，因此，应当慎重选择zz的应力初始化; boundary conditions 边界条件bound stress -33.8e6 0 -67.5e6 range -0.01 80.01 49.99 50.01bound xvel 0 range -0.01 0.01 -0.01 50.01 ；x=0bound xvel 0 range 79.99 80.01 -0.01 50.01 ；x=80bound yv

9、el 0 range -0.01 80.10 -0.01 0.01 ；y=0;hist solve_rat type 1 ；最大的力学比率hist unbal ；记录最大不平衡力历hist syy 10.05 2.2hist syy 20.05 2.2hist syy 30.05 2.2hist syy 40.05 2.2hist syy 50.05 2.2;hist ydis 5,2.5 ； 记录坐标（x=5, y=2.5）附近结点y方向的位移hist ydis 15,2.5hist ydis 25,2.5hist ydis 35,2.5hist ydis 45,2.5;set pline

10、0 1.25 10 1.25 10 ；set pline set pline 0 1.25 20 1.25 20set pline 0 1.25 30 1.25 30set pline 0 1.25 40 1.25 40;solve rat 1e-5save fill0.sav; excavate chamber 1 开挖第一个空间zone model null range 0 10 0 2.5 ；开挖zone she 600e6 bul 450e6 range 0 10 0 2.5 ；回填材料参数solve rat 1e-5save fill1.sav; backfill chamber 1

11、zone model dy range mat 2zone she 600e6 bul 450e6 fric 40 dil 5 range mat 2zone cap 1e4 cptable 2 range mat 2tab 2 0,1e4 0.02,.2e6 .04,.8e6 .06,1.3e6 .1,5.5e6tab 2 .12,8.5e6 .14,11.5e6 .16,19e6 .18,34e6 .2,50e6; excavate chamber 2zone model null range mat 3zone she 600e6 bul 450e6 range mat 3solve r

12、at 1e-5save fill2.savplot block stress pline 1 sig1 fill; backfill chamber 2zone model dy range mat 3zone she 600e6 bul 450e6 fric 40 dil 5 range mat 3zone cap 1e4 cptable 2 range mat 3; excavate chamber 3zone model null range mat 4zone she 600e6 bul 450e6 range mat 4solve rat 1e-5save fill3.savplot

13、 block stress pline 2 sig1 fill; backfill chamber 3zone model dy range mat 4zone she 600e6 bul 450e6 fric 40 dil 5 range mat 4zone cap 1e4 cptable 2 range mat 4; excavate chamber 4zone model null range mat 5zone she 600e6 bul 450e6 range mat 5solve rat 1e-5save fill4.savplot block stress pline 3 sig

14、1 fill; backfill chamber 4zone model dy range mat 5zone she 600e6 bul 450e6 fric 40 dil 5 range mat 5zone cap 1e4 cptable 2 range mat 5; excavate chamber 5zone model null range mat 6zone she 600e6 bul 450e6 range mat 6solve rat 1e-5save fill5.savplot block stress pline 4 sig1 fillret命令：;屏幕观看下沉与应力等值线

15、,int是指intervalpl bl ydisp int 0.01pl bl syy int 1e6set pline x1 y1 x2 y2 nplot pline 线号 syy sxx 或 ydis1.把图形保存下来，能在autocad中打开、编辑plot block cable red supp ye stru bl;显示块体、锚杆（红色）、支架（黄色）、梁（蓝色）的图形set plot 256 dxf；设置图形为256色set out c:ss.dxfcopy c:ss.dxf；把这图形以ss.dxf文件保存在c磁盘下2.把数据导出set log on print pline 1 y

16、dis;把pline 1的y方向位移的数据导出到udec目录下udec.txt文件中set log off 3.plot block stress ;显示块体的应力plot block dis ;显示块体的位移(有x和y方向)plot block pl ;显示块体的塑性区(plastic)savexx.save ;保存计算结果res xx.save ;调用4.液压支架的命令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

17、相对应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.5e9st_prat=0.14st_den=7800prop mat=16 st_yield=6e7 st_yr

18、esid=6e7 st_ycomp=2.5e7prop mat=16 if_kn=1.35e9 if_ks=1.35e9if_tens=0 if_fric=18if_coh=0;interface-界面1.把图形保存下来，能在autocad中打开、编辑plot block cable red supp ye stru bl;显示块体、锚杆（红色）、支架（黄色）、梁（蓝色）的图形set plot 256 dxf；设置图形为256色set out c:ss.dxfcopy c:ss.dxf；把这图形以ss.dxf文件保存在c磁盘下2.把数据导出set log on print pline 1 ydis;把pline 1的y方向位移的数据导出到udec目录下udec.txt文件中set log off 3.plot block stress ;显示块体的应力plot bloc