考虑岩块接触面流变的离散元模型.pdfVIP

黪釜一 生薹耋兰虹： 差墨迸篙 606 考虑岩块接触面流变的离散元模型 王刚金峰张楚汉王光纶 清单太学水刺水电工程系 采京 100084 擅要本文提出了一种模担岩块接■面藏空特性的膏散元数值模型．在模蠹中，引入积分垂藏 变本构方翟和桔性等效接触力增量研究块体一接■面系统反应，并通过平衡遗代计算离靛块体系 统由于接触面瘫变变形所引起的相应位移．通过对一个简单算倒计算值和理论解的对比，验证了 本模型有良好的耩度，具有应用于岩质高边坡和青敢块体系统流变分析的重要前景。 关■词离散单元法，藏变，接触面 1．引 言 离散单元法(DEM)是一种适用于模拟离散介质的数值方法。自CundaU于70年代提出以 来，离散元方法在岩石介质数值分析中得到了迅速发展，它与有限元、DDA等方法一起，成为 岩体变形与稳定分析的有力工具。 2．离散单元法静动力计算的基本原理 离散单元法的突出特点是它不但能够考虑离散块体本身的变形，而且还能模拟接触面的 张开和滑移等力学行为，材料可以应用非线性本构关系。由于采用显式时步迭代，对应于应力 增量，本构方程采用增量形式，离散元计算方法对地震、爆破等非线性动力学问题具有内在的 适应性。 如图1(a)所示，岩体被裂隙切割成块体系统。裂隙接触面的力一位移关系通过法向和切向 弹簧来模拟。接触力在法向应满足零拉力条件，在切向服从coulomb准则，即： F。一^“ 若“。≤0 F，一0 若。0 … 、一 F，，=五凶 若I再乙l一，lF二l+cL F。一sign(也)(，ILI+正)若l凡I，lF盘l+止 式中，F。，凡分别为接触力向量B的法向和切向分量I^。，l，分别为接触面的法向和切向刚 度；“．，“，是接触面法向和切向的相对位移；，，c是接触面的摩擦系数和粘着力；工表示接触面 作用长度。 块体本身为可变形体，块体内部划分三角形常应变差分单元。如图l(b)所示，网格结点N 为差分单元B。，B：，…，B．的公共结点。三角形差分单元把惯性质量平均分配给其三个结点， 即把虚线r。所围起的多边形的质量分配给N结点。因为每个三角形差分单元的应力为常量， ·373· 结点弹性力积分路径可简化为将每个三角 形单元对应边中点依次连接而成的实线路 径r。由于块体弹性变形引起的结点Ⅳ的 弹性力向量只按下式计算： r r E=I叮·^ds—I口·ndj(2) ≠。 { 鞫阚⑩力 其中，口是积分边界所在差分单元的应力张 h，¨· 量，n是积分边界的法向单位矢量，rl和11 (a)离散系统 (b)局部放大 分别为图1(b)中虚线和实线所示的多边形 图l可变形体离散元系统 边界。