锚杆作用下岩体结构面的力学和变形特征.docxVIP

锚杆作用下岩体结构面的力学和变形特征 结构面在岩石工程中广泛存在，破坏了岩体的连续性和完整性，导致了不同的异质性和各向异性，并控制了岩体的稳定性。这对于研究结构面的剪切强度和变形特征具有重要意义。许多学者通过人工制备不同起伏程度的结构面对其特性进行研究,如蒋宇静等研究了不同接触状态下的岩石断裂节理试件剪切过程中力学性质的变化情况;李海波等利用锯齿状的混凝土岩石节理试样,研究不同剪切速率下各种岩石节理起伏角度岩石节理的强度特征;杜守继等基于不规则的人工岩石节理经历不同剪切变形历史的剪切试验,分析了岩石节理剪切变形特性及与变形历史的依存关系。这些研究中的结构面主要处于原始状态,未进行加固处理;而实际工程中,设计者为了稳定岩体往往采用各种手段对其进行支护,其中,全长黏结式锚杆作为岩体的重要加固手段已得到广泛的运用;但是,锚杆加固情况下结构面的力学和变形效应的研究还较少,主要是由于结构面锚杆加固模型制备条件的限制;而近年来,随着计算机技术的不断发展,采用数值模拟方法分析岩土工程问题成为新的趋势,从而为结构面锚固特性研究开辟了新的道路。其中,快速拉格朗日元法(FLAC3D)适用于绝大多数的工程力学问题,尤其适用于材料弹塑性以及几何边界复杂情况的数值模拟。因此,本文拟采用FLAC3D分析三维结构面试样与锚杆的耦合效应,研究不同表面形态和锚杆倾角对结构面强度和变形特征的影响。 1 结构面起伏角模型 由于FLAC3D在建模以及单元网格划分等前处理问题上存在一定困难,因此,本文首先利用ANSYS建立数值模型,然后,通过自编的ANSYS-FLAC3D的转换程序,导入FLAC3D进行计算。结构面试样模型尺寸为4 m×4 m×4 m,将结构面视为软弱结构面,其厚度为0.1 m;为了探讨不同几何形态结构面的特性,设置结构面起伏角β=0、17°,如图1所示,相应模型的节点数为1 951、1 673,单元数分别为9 420、7 485。 为了分析锚杆对于结构面剪切效应的影响,在结构面中央布置不同倾角的锚杆单元,岩体与锚杆的相互作用采用弹簧单元进行模拟,可分别考虑锚杆的轴向拉伸作用,轴向摩擦作用以及弯曲作用,具体单元性质见文献。 模型的边界条件:底部固定约束,下部岩体的侧面分别约束水平方向的位移,模型上部为自由边界。数值试验过程中,模型顶部加载法向应力,上部岩体的施加位移荷载,位移大小为1×10-3mm/计算时步,使其发生移动。 计算采用Mohr-Coulomb准则,岩体参数:弹性模量2.0 GPa,泊松比0.2,黏结力0.8 MPa,容重23.0 kN/m3,膨胀角32.0°,内摩擦角37.0°,抗拉强度0.4 MPa;结构面参数:弹性模量200 MPa,泊松比0.3,黏结力0.050 MPa,容重19.0 kN/m3,膨胀角20.0°,内摩擦角24.0°,抗拉强度0.05 MPa。锚杆长度3 m,倾角θ分别为15、30、45、60、75、90°。 2 分析与讨论 2.1 剪切位移对残余力的影响 直剪试验过程中,在加锚和未加锚情况下,结构面的应力变形特征(图2)显示,二者剪应力-剪位移关系曲线在达到峰值前较为一致,为接近直线的试样弹性特征阶段;当曲线达到峰值后,则表现出较为不一致的特征:对于无锚状态的结构面,曲线峰值后表现为应变软化特征,随着剪切位移的增大剪切应力逐渐减小;当黏结键均受到破坏后,试样发生滑动,此时提供的抵抗力为结构面和上下岩体之间的相互摩擦,因此出现滑移曲线,即随着剪切位移的增加,试样的剪切应力基本保持不变,此时的剪切应力为试样的残余强度。而对于加锚状态下的结构面,曲线峰值后表现为应变强化特征,结构面的峰值剪切应力和峰值剪切强度得到一定提高。 2.2 moh-coolomb法 为了对比结构面加锚状态下的剪切强度的变化规律,记录不同结构面起伏角下(0、17°),结构面正应力和试样剪切强度之间的关系,如图3所示。 从图3中可以看出,结构面加锚后试样的剪切强度得到提高,并且试样剪切强度τs与结构面正应力σn之间呈现显著的线性关系,可用线性Mohr-Coulomb准则进行拟合,建立拟合方程,即 τs=c+σntanφ(1)τs=c+σntanφ(1) 式中,c和φ分别为结构面的等效黏结力和内摩擦角。 通过曲线拟合得到表1。从表1可以看出,当起伏角为0和17°时,结构面加锚后主要提高结构面的等效黏结力,提高程度分别为26.82%和48.42%,而等效内摩擦角较加锚前略有降低;可见当结构面起伏角较小时,锚杆主要发挥拉拔作用,从而等效提高了结构面的黏结力。 2.3 锚杆抗剪强度 记录结构面起伏角为0时,不同锚杆倾角下试样的剪切强度如图4所示。对于结构面起伏角为0的试样剪切试验,锚杆倾角为15°时对应的试样剪切强度最大;随着锚杆倾角的增大,试样的剪切强度