岩石缺陷位置不同时剪切带图案启动、演变及试样的宏观力学行为.docxVIP

岩石缺陷位置不同时剪切带图案启动、演变及试样的宏观力学行为 1 数值模拟方法在基本粒子局部化中的应用进展 近年来，局部适应性、间隙、适应性软化和稳定性一直是国际固体力学和岩浆岩力学领域的研究重点之一。这是描述岩石介质实际破坏过程的理论。虽然这项研究起步不久，但对判断岩土工程的失稳与破坏起着重大的作用，因而必将成为岩土塑性理论中的重要组成部分。国内的许多学者在这方面也做了许多有益的工作[2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17]，文献对应变局部化的多种理论及数值模拟方法研究进展进行了比较详细的论述。最近，有人基于梯度塑性理论，在非均质工程材料（岩石、金属及混凝土）的变形、破坏及稳定性等方面进行了一些探索，但所采用的单轴压缩、单轴拉伸及直接剪切试验的力学模型相对都比较简单，离真实岩土材料受力状态还有不少差距。 为了深入理解平面应变状态下岩样在压缩载荷作用下的破坏机制及过程，文献[20,21,22,23,24,25,26]采用拉格朗日元法（FLAC）研究了岩样尺寸、围压、应变率、扩容及端面约束对局部化图案的影响。 本文通过在岩样左边界及右边界的不同位置预制材料缺陷，激发了不均匀的变形场。从剪切应变增量的二维等值线图出发，研究了右缺陷位置不同时剪切带图案的启动、演化规律及全程应力-变形曲线的性状，发现了剪切带的转折及反射现象，对剪切带图案与应力-变形曲线软化段性状的关系进行了探讨。 2 本构模型及材料缺陷的强度面的确定 计算模型的几何尺寸、单元划分及边界条件见图1。试样下端面被约束，在试样上端面施加常速度，v=2×10-9m/时步，计算是在小变形模式及平面应变状态下进行。试样上、下端面的岩石质点仅允许向下运动，其它方向的运动都被约束，这对应于上、下端面与试验机压头之间存在较大摩擦力的情形。 在弹性阶段，岩石的本构关系为线弹性，体积模量取15 GPa，剪切模量取11 GPa。峰值强度后岩石的本构模型取莫尔-库仑剪破坏与拉破坏复合的应变软化模型，抗拉强度取200 kPa。凝聚力、内摩擦角与塑性应变的关系见图2，初始凝聚力取275 kPa，初始内摩擦角取44°。 为了得到不对称的剪切带图案，在试样左边界靠近下端面的位置及右边界的不同位置预制材料缺陷。具有材料缺陷的单元是一个空单元。计算采用10个方案。在方案1中，第21及780号单元被置空，在方案2中，第21及700号单元被置空，依此类推，直到在方案10中，第21及60号单元被置空。 3 对结果的分析与讨论 3.1 应力-变形本构模型 图3～10分别是方案1、方案3～6和方案8～10不同时间步的剪切应变增量二维等值线图的数值结果，图3～10的时间步（个）分别为1 500,2 000,2 200,2 400,2 600,2 800,3 000及3 200。限于篇幅，方案2及方案的数值结果不再给出。 图3～10中某些位置的等值线越密，表示位于这些位置的单元具有的剪切应变增量越大。时间步数t与试样上端面的压缩位移s的关系，s=vt。因此，很容易将时间步数换算成试样上端面的位移或轴向应变。图3～10中各图的横、纵坐标轴上的数字（1～20,1～40）分别表示在水平及垂直方向上各单元。 图11给出了方案1～6的应力-变形曲线的一部分。图12给出了方案6～10的应力-变形曲线的一部分。在这两个图中，岩样上端面的压缩位移小于4×10-6m的曲线均被省略，是由于缺陷的具体位置对应力-变形曲线的这一部分影响甚小，图中数字表示各个方案。 3.2 折射波发生的折 对于方案1（图3），在1 500个时间步以下（包含1 500个时间步），剪切应变增量场是比较均匀的。当时间步达到2 000个时（图3(d)），一条剪切应变集中带（称之为局部化带或剪切带）出现在试样右边界缺陷附近。位于试样左边界缺陷附近的单元也存在一定的剪切应变集中现象，但集中程度较弱。随着时步的增加，见图3(a)～图3(f)，使得试样右边界缺陷附近的剪切带长度不断增加，沿两缺陷连线方向传播，直到贯穿试样的左边界。 对于方案3，启动于试样右边界缺陷附近的剪切带向试样内部传播过程中，见图4(d)～图4(e)，先保持直线状态，剪切带的倾角保持为一个值，这一点与方案1相同。但是，当剪切带的前端接近试样的纵向对称线时，剪切带的方向明显发生了转折，剪切带的切向与水平轴之间所夹的锐角（剪切带倾角）明显变小，见图4(f)～图4(h)），这一点与方案1有所不同。笔者认为，出现这种现象的原因是强烈的端面约束的作用，使接近岩样下端面的三角形区域（“Δ”）内部的变形的发展受到阻碍，因而，剪切带在这一位置不能形成，从而迫使剪切带的方向发生了转折。本文称上述现象为剪切带的转折现象（或称之为折射现象）。 对于方案4（图5(c)），在试样右边界缺陷