UDEC实例详解1.doc

实例1 荷载作用下边坡稳定问题 朔准线黄河大桥，边坡高105m，层理间距取3m，岩层视倾角6°；主节理间距取2m，视倾角85°，计算模型如下： DEM计算模型 相关计算参数取值如下：桥基长15m，宽12m，荷载180000kN；灰岩弹性模量E=4×104MPa，泊松比μ=0.23，岩体密度γ=2.7×103kg/m3。 new round 0.05 set delc off block 0 0 0 50 50 50 62.6 115.2 68.9 123.9 85.5 132.5 105.5 140.1 125.4 156.4 180 156.4 180 0 jset 85,0 300,0 0 0 3,0 (50 50) jset 6 0 300 0 0 0 4 0 (50 50) gen auto 5 change jmat=1 range angle 84 86 change jmat=2 range angle 5 7 prop mat=1 den=2700 b=1.5e9 s=0.6e9 prop jmat=1 jkn=2e9 jks=2e9 jcoh=1e9 jfr=25 jten=1e5 节理参数不同 prop jmat=2 jkn=2e9 jks=2e9 jcoh=1e10 jfr=35 grav 0 -10 hist n=100 xdis 68.9 123.9 ydis 68.9 123.9 hist n=100 xdis 125.4 156.4 ydis 125.4 156.4 hist unbal bound xvel=0 range -1 0.2 -1 50.1 bound xvel=0 range 224.9 225.1 -1 157 bound yvel=0 range -1 226 -1 0.1 cyc 5000 save nature.sav 变坡点x位移 变坡点y位移 坡顶点x位移 坡顶点y位移 不平衡力  实例2 边坡渗流稳定问题 问题定义： ; ================================================ ; --- fluid flow test run --- ; --- slope : 10 m high --- ; --- 2 joint 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 23,10 23,-5 生成一个块体，每次分析该命令只能用一次。 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 定义材料属性：密度(density)、体积模量(bulk)，剪切模量(shear)（当cons=3、6时，还需要cohesion粘聚力、dilation剪胀角、friction内摩擦角、tension拉张强度） prop mat=1 jkn=10000 jks=10000 jfric=45.0 prop mat=1 jperm=1.0e8 azero=0.0005 ares=0.0002 ; 定义材料属性：法向刚度系数，切向刚度系数，摩擦角（度），节理渗透系数、法向应力为零时的裂隙开口大小、高应力作用下裂隙开口大小（长度） Jcons=2 (面面接触，库仑滑动准则)还有一些参数：如：jcohesion，jdilation，jtension Jcons=1（点接触）参数：ccohesion(内聚力)，cdilation(剪胀角)，c