1、实例 1 荷载作用下边坡稳定问题朔准线黄河大桥，边坡高 105m，层理间距取 3m，岩层视倾角 6；主节理间距取 2m，视倾角 85，计算模型如下：DEM 计算模型相关计算参数取值如下：桥基长 15m，宽 12m，荷载 180000kN；灰岩弹性模量 E=4104MPa，泊松比 =0.23 ，岩体密度 =2.7103kg/m3。newround 0.05set delc offblock 0 0 0 50 50 50 62.6 115.2 68.9 123.9 85.5 132.5 =; - fluid flow test run -; - slope : 10 m high -; - 2 j
2、oint sets : 20 and 80 deg.; friction = 30 deg.; - r.h.s. water level : 6 m - no failure -; 8 m; 9 m; 10 m - failure -; =;round 0.05定义刚体角边圆角半径，缺省值0.5，建议为块体边长的1%，设定后不可更改。set delc off如果接触块体张开距离大于 0.55 倍 round 定义的长度时，接触将被自动删除（缺省值 on） ，off 则表示不自动删除。block 0,-5 0,0 5,0 11,10 22,10 22,-5生成一个块体，每次分析该命令只能用一次。
3、jset 20,0 100,0 0,0 2,0 (5,1)设置节理性质，20 为角度，100 为长度，2 为间距， （5，1 表示起点位置）jset 80,0 100,0 0,0 3,0 (5,0);del area 0.1;删除面积小于 0.1 的块体gen auto 10;prop mat=1 dens=0.0025 k=16667 g=10000定义材料属性：密度(de nsity)、体积模量( bulk)，剪切模量(shear)（当cons=3、6 时，还需要 cohesion 粘聚力、dilation 剪胀角、friction 内摩擦角、tension 拉张强度）prop mat=1
4、 jkn=10000 jks=10000 jfric=45.0prop mat=1 jperm=1.0e8 azero=0.0005 ares=0.0002; 定义材料属性：法向刚度系数，切向刚度系数，摩擦角（度） ，节理渗透系数、法向应力为零时的裂隙开口大小、高应力作用下裂隙开口大小（长度）Jcons=2 (面面接触，库仑滑动准则 )还有一些参数：如：jcohesion，jdilation，jtensionJcons=1（点接触）参数： ccohesion(内聚力)，cdilation(剪胀角)，cfriction(摩擦系数)， cperm(渗透系数)， ctension(拉张强度)，kn(
5、法向刚度)，ks( 切向刚度)Jcons=5(面面接触，库仑滑动准则，残余强度) ，与 jcons=2 相比，多了以下参数：jrescoh(残余内聚力 )，jrfric(残余摩擦角)，jrtens(残余拉张强度)对于 jcons=2 和 5，可以指定非线性应力应变关系，用 property nstable n 命令，Table n (x,y) (x1,y2).change jmat=1 jcons=2;设置材料对应关系，如节理采用编号为 1 的材料，节理性质为2（1，2，3，5，7）insitu str -0.125 0 -0.25 ygrad 0.0125 0 0.025;初始化应力，str
6、 表示(0,0)处应力，sxx0,sxy0,syy0, 分别表示 x 方向应力，剪应力，y 方向应力Ygrad 表示变化梯度，sxxy,sxyy,syyySxx=sxx0+sxxx*x+sxxy*ySxy=sxy0+sxyx*x+sxyy*ySyy=syy0+syyx*x+syyy*y还可以给出 Xgrad，同上。bound xvel=0 range -1,1 -6 1bound yvel=0 range -1,31 -6 -4bound xvel=0 range 21,23 -6 11;定义边界条件外部边界条件：1）荷载边界：bound xload fx range x,ybound ylo
7、ad fy range x,y,2) 自由边界：bound xfreebound yfree3) 应力边界：bound stress sxx0,sxy0,syy0 xgrad sxxx,sxyx,syyx ygrad sxxy,sxyy,xyyy4) 位移速度边界 bound gvel vx0,vy0,vxx,vxy,vyx,vyyVx=vx0+vxx*x+vxy*yVy=vy0+vyx*x+vyy*yBound nvel vn 变形块体法向位移速度Bound svel vs 变形块体剪切位移速度Bound xvel vx 变形体 x 方向位移速度Bound yvel vy 变形体 y 方向位
8、移速度内部边界条件：1）荷载边界：bound interior xload fx range x,ybound interior yload fy range x,y,2) 自由边界：bound interior xfreebound interior yfree3) 应力边界：bound interior stress sxx0,sxy0,syy0 xgrad sxxx,sxyx,syyx ygrad sxxy,sxyy,xyyy4) 位移速度边界 bound interior gvel vx0,vy0,vxx,vxy,vyx,vyyVx=vx0+vxx*x+vxy*yVy=vy0+vyx*
9、x+vyy*yBound interior nvel vn 变形块体法向位移速度Bound interior svel vs 变形块体剪切位移速度Bound interior xvel vx 变形体 x 方向位移速度Bound interior yvel vy 变形体 y 方向位移速度grav 0 -10;x、y 方向的重力加速度，同 set gravity x, yhist n=100 xdis 11,10 ydis 11,10保存变量历史，n=100 表示第循环 100 时间步，存储一次，Xdis 11,10 表示存储距 11，10 最近的点的 x 方向位移；Ydis 11,10 表示存储
10、距 11，10 最近的点的 y 方向位移；hist unbal; 记录最大不平衡力Solve求解save slfl1.sav;保存文件; -; flow - r.h.s. water at y=6m; - no failure -; -fluid dens 0.001定义流体密度，注意单位要统一bound imperm range -1,31 -6 -4 设置无渗流的边界范围，x 最小最大范围，y 最小最大范围bound pygrad -0.01 range -1,1 -6 1 设置流体压力梯度范围bound pp=0.06 pygrad -0.01 range 21,23 -6 6; 设置流
11、体压力梯度范围，pp 表示压力初值，pygrad 表示梯度 range 表示范围。set flow steady;设置稳定流模式hist xdis 5.92 1.54 ydis 5.92 1.54记录位移历史reset displ jdis设置块体和节理位移为零prop mat=1 jfric=25设置材料 1 的节理摩擦角为 25 度。cyc 500执行时间 500 步。save slfl2.sav; -; flow - r.h.s. water raised to y=8m; - some slip, but no block failure -; -bound pp=0.08 pygrad -0.01 range 21,23 -6 8cyc 500save slfl3.sav; -; flow - r.h.s. water raised to y=9m; - no failure -; -bound pp=0.09 pygrad -0.01 range 21,23 -6 9 cyc 1000save slfl4.sav; -; flow - r.h.s. water raised to y=10m; - failure -; -bound pp=0.10 pygrad -0.01 range 21,23 -6 11cyc 15000save slfl5.savreturn
