UDEC指令12.4.doc

;new title moni ;定义名称moni round=0.025 set ovtol=0.4 ;常规设置 bl 0,0 0,100 300,100 300,0 ;定义范围，四点坐标，顺时针方向 cr 0,15 300,15 ;画直线，两点坐标 jreg id 1 0,0 0,15 300,15 300,0 delete ;定义节理，命名1，定义范围，四点坐标，delete为常规语言 jset 90,0 5,0 5,0 5,0 0,0 rang jreg 1 ;定义节理角度90°；节理的长度；节理的距离，即纵向间隔；横向距离，即横向隔5一个；起始点坐标为0,0 jset 90,0 5,0 5,0 5,0 2.5,5 rang jreg 1 ;定义节理的另外一项 gen quad 10 range 0 300 0 100 ;当0 300 0 100范围内的块体的对角长度小于10时就是可变形体。 zone model mo range 0 300 0 100 ;定义摩尔库伦模型的范围，X的范围，Y的范围 change jcons=5 range 0 300 0 100 ;定义节理的某个属性 change mat=1 range reg 0,0 0,15 300,15 300,0 ;定义物质1的范围，即赋予这个范围为物质1,只赋名，没有定义属性 change mat=2 range reg 0,15 0,20 300,20 300,15 change mat=3 range reg 0,20 0,23 300,23 300,20 change mat=4 range reg 0,23 0,28 300,28 300,23 change jmat=1 range reg 0,0 0,15 300,15 300,0 ;定义节理1的范围 set jmatdf 2 ;定义节理属性，使符合摩尔库伦准则普适。2代表摩尔库伦准则 save ch.sav ;保存文件 cal prop-25.txt ;调出文件 res ch.sav ;调出前面保存的文件接着计算 ;参数设置 prop mat=1 dens=2720 ;物质的密度 zone k=12.12e9 g=10.26e9 fric=33 coh=5.77e6 ten=3.6e6 range mat=1 ;定义物质1的体积模量、剪切模量、内摩擦角、粘聚力、抗拉强度 prop jmat=1 jkn=11e8 jks=12e7 jfric=12 jcoh=0 jten=0 ;定义节理的属性法向刚度、切向刚度、内摩擦角、粘聚力、抗拉强度 ;边界条件 boun stress 0 0 0 range 0,300 99.9 100.1 ;定义边界条件，三向应力为0，上边界，范围，纵向波动范围0 0 0（x、xy、y） boun stress 0 0 -8.6095397e6 range ;定义边界应力条件，竖直方向的应力为，负号表示方向向下 insitu str -7.5e6 0 -15e6 szz -7.5e6 ygrad 1.1018e4 0 2.203607e4 ;定义应力属性，sxx方向即水平应力，sxy方向即剪切应力，syy方向即竖直应力，szz方向的应力，ygrad应力梯度 set gravity 0,-10 ;重力加速度 bound xvel=0 range -0.1 0.1 0 100 ;X位移边界,位移波动范围，范围 bound xvel=0 range 299.9 300.1 0 100 ;X位移边界，同上 bound yvel=0 range -0.1 300.1 -0.1,0.1 ;Y位移边界 hist solve_ratio type 1 ;设置不平衡力的精度，普适 s