岩石结构面注浆前后抗剪特性数值分析

岩石结构面注浆前后抗剪特性数值分析
吴乐文;赵奎
【摘要】岩石结构面力学特性研究一直是岩石力学的热点问题，为研究注浆前后岩石结构面力学特性，利用ANSYS建立了结构面剪切计算模型并导入FLAC3D 进行计算，采用等效参数模拟方法对比分析了3种结构面类型在台阶高宽比为0.5和0.3下结构面的力学特性.结果表明：台阶宽高比不同台阶所表现的破坏规律不同，高宽比大时为剪断破坏，高宽比小时为压切破坏；注浆对提高结构面抗剪强度作用明显.最后，根据试验结果提出了非控制性爬坡角概念.%The researches on the mechanical property of rock structure plane has always been a hot topic .In order to study the mechanical property of rock mass at pre-and-post grouting, a FLAC3D calculation model is constructed based on ANSYS to perform calculation. The mechanical properties of three structural planes and two aspect ratio types are analyzed with aspect ratio are 0.5 and 0.3.The results show that the different step-shaped aspect ratios have the different failure laws, which exhibited shear failure with large aspect ratios and pressure-shear failure with small aspect ratios, and the shearing strength increased markedly after grouting . The concept of uncontrolled climbing angle is put forward according to the results .
【期刊名称】《有色金属科学与工程》
【年(卷),期】2014(000)004
【总页数】5页(P96-100)
【关键词】数值模拟;注浆;结构面;剪切强度
【作者】吴乐文;赵奎
【作者单位】广西华锡集团股份有限公司，广西柳州 545006;江西理工大学资源
与环境工程学院，江西赣州341000
【正文语种】中文
【中图分类】TD315
0 引言
许多国内外学者针对控制着工程岩体稳定性的岩体结构面力学特性这一热点问题进行了持续性研究.众所周知，结构面的存在对岩体力学性质有着决定性影响，因此，研究岩体结构面力学性质是非常有意义的.
孙广忠教授在文献[1]中归纳了结构面的4种起伏形态分别为平直结构面、台阶状
结构面、锯齿状结构面、波状结构面，4种结构面的力学效应是不同的.除了锯齿
状结构面，其他3种结构面力学特性的研究已较多，一些学者对岩体结构面锚固
特性采用室内试验或理论分析进行了相关研究[2-4].但是对结构面注浆加固前后力学特性的研究还比较少，尤其是采用数值分析方法的研究更少.因此，本文通过3
种类型和2种台阶高度结构面剪切模拟试验，对结构面注浆前后抗剪特性进行了
分析研究.该研究对揭示其加固机制具有重要的意义.
1 试验与计算模型
通过ANSYS有限元分析软件建立数值分析模型，再通过转换导入FLAC3D软件
中进行分析和计算[5-12].为尽可能得到岩石结构面剪切特性的规律，模型尺寸：
长×宽×高为100 mm×100 mm×100 mm.为了模拟结构面，在模型中间设置厚
度为2 mm的结构面模拟其注浆前后效果.设置的结构面类型包括1个台阶、2个
台阶以及组合（1个台阶与1个锯齿）3种类型，台阶宽度为20 mm，每种类型
包括2种台阶高度，分别为10 mm、6 mm，则3种类型结构面模型高宽比分别
为0.5和0.3.如图1～图3所示分别为3种结构面模型示意图.
图1 1个台阶结构面模型
图2 2个台阶结构面模型
图3 组合结构面模型
边界条件设置为底部固定约束，下半部分模型四周约束水平方向位移，上半部分顶部为自由边界.
计算采用Mohr-Coulomb准则，在模型顶部施加法向荷载，上半部分岩石施加水平位移荷载，荷载大小2×10－4 mm/时步.模拟试验中只考虑法向荷载，因为法
向荷载远远大于岩石自重，所以计算略去岩石自重[13-15].如表1、表2所示，岩石和结构面参数是根据RMT-150C型岩石力学试验系统对岩石试块进行室内单轴和剪切试验的结果选取.
表1 岩石参数选取参数名称内摩擦角弹性模量/GPa参数取值 45.7 0.278 10.0 20.0 15.5/（°）泊松比抗拉强度/MPa 黏结力/MPa
表2 结构面参数选取参数名称参数取值未注浆注浆黏结力/MPa 0 15.0弹性模量/MPa 0.02 15.5泊松比 0.000 2 0.027 8抗拉强度/MPa 0.001 1.0内摩擦角 /（°）10.0 40.0
2 破坏规律分析
通过数值计算分析讨论台阶状结构面剪切破坏形式以及规律，探讨台阶数量及结构面形态与抗剪强度关系、施加法向应力大小与抗剪强度关系.由于试验结果图太多，此处只列出图4～图6，分别为试验中高宽比0.5和0.3时，结构面最大、最小主
应力云图以及组合剖面塑性云图.本文中所有图示除非特别说明，剪切方向均为从
通过分析图4、图5、图6等1个台阶、2个台阶、组合结构面模型剖面的最大主应力、最小主应力以及塑性区云图得到下述研究结果.
2.1 台阶高宽比为0.5情况分析
1）拉应力区域皆形成于3种结构面的台阶及锯齿下部，拉应力最大处位于底部. 2）结构面上半部分岩石在剪切作用下有往上抬起效应，导致台阶右边下部岩石受到明显压应力作用，台阶右边底角转折处形成压应力集中区域.
3）组合类型岩石在剪切作用下锯齿下方和锯齿顶部均表现为拉应力.而压应力较集中区域为锯齿正下方与台阶底部，且在锯齿右面底部达到最大压应力，说明剪切过程上部岩石有往上抬起效应.
图4 1个台阶剖面最大主应力云图
图5 2个台阶剖面最小主应力云图
图6 组合剖面塑性区云图
4）由塑性云图6可知，因岩石抗拉强度远远小于岩石抗压强度，导致台阶底部在剪切应力和法向应力的作用下最先表现为拉伸破坏.组合类型模型中锯齿和台阶的底部都已开始出现拉伸破坏，且最大拉应力区域与破坏区域一致，台阶最终因张剪滑移而破坏.
2.2 台阶高宽比为0.3情况分析
为了对比分析台阶高宽比不同的情况下3种结构面剪切破坏规律，进行了台阶高宽比为0.3时结构面剪切试验.
通过分析台阶高宽比为0.3的1个台阶结构面模型剖面的最大主应力、最小主应力以及塑性区云图可知：
1）台阶高度降低后，台阶高宽比为0.3时结构面主应力云图与台阶高宽比为0.5时的应力云图规律是一样的.台阶底部受到拉应力作用，拉应力集中区域出现在台
2）沿顺剪方向，台阶只在左下部分出现斜向上的拉伸破坏痕迹，表明当台阶高宽比较小时，台阶一般因剪切作用产生压切破坏，且对于2个台阶结构面模型压切破坏同时发生于2个台阶.因此，当台阶高宽比较小时，台阶因剪切作用而产生压切破坏.破坏区域分布与文献[1]分析一致.
3）组合结构面模型中拉应力集中区域位于台阶和锯齿的底部，台阶底部受弯矩作用明显，由此可见：高宽比越大，台阶状结构面张性破坏效应越明显.锯齿底部的拉应力明显小于台阶底部的拉应力；锯齿高度降低后，锯齿下部的拉应力集中效应降低了，由此可见：起伏角α越大，结构面张性破坏效应越明显.
4）2个台阶结构面与1个台阶结构面模型类似，沿受剪方向，出现于台阶左下角的拉伸破坏区逐步向台阶顶端扩展，最终导致台阶形成压切破坏.锯齿下方拉伸破坏区在左下与右下部位明显.
3 法向应力与抗剪强度分析
通过对3种类型结构面注浆前后的试验过程监测得到了法向应力为 1 MPa、2 MPa、3 MPa、4 MPa、5 MPa下峰值剪切应力的大小，如表3所示.
表3 不同法向应力对应峰值剪切应力/MPa结构面类型注浆情况法向应力1 2 3 4 5 1台阶未注浆 1.821 3.642 5.50 7.291 9.080注浆 2.078 4.024 5.973 7.921 9.870 2台阶未注浆 2.250 4.500 6.280 9.000 11.300注浆 2.515 4.875 7.238 9.601 11.970组合未注浆 1.868 3.740 5.609 7.482 9.352注浆 2.318 4.443 6.572 8.694 10.820
通过数值模拟验证了法向应力与结构面抗剪强度的关系特性.通过数值计算结果可知：抗剪强度和法向应力拟合相关系数近似为1.0.结构面抗剪强度随法向应力增大而呈线性增长.
4 注浆与未注浆抗剪强度对比
依据文献[4]的研究可知，数值模拟中等效参数方法可以模拟结构面注浆效果，因此，通过改变中间层参数模拟注浆效果，结果如图7～图9所示.
图7所示为1个台阶结构面注浆与未注浆结构面抗剪强度对比条形图，注浆后结
构面抗剪强度在5个等级法向应力下提高幅度明显，分别达到14.11%、10.49%、8.60%、8.64%、8.70%.注浆后强度平均提高了10.11%.
图8所示为2个台阶结构面注浆与未注浆结构面抗剪强度对比条形图，注浆后结
构面抗剪强度在5个等级法向应力下提高幅度明显，分别为11.78%、8.33%、13.24%、6.68%、5.93%.注浆后强度平均提高了9.19%.
图9所示为组合结构面注浆与未注浆抗剪强度对比条形图，注浆后结构面抗剪强
度在5个等级法向应力下提高幅度明显，分别为24.09%、7.03%、18.80%、16.20%、15.70%.注浆后强度平均提高了16.35%.
图7 1个台阶结构面注浆与未注浆强度对比
图8 2个台阶结构面注浆与未注浆强度对比
图9 组合结构面注浆与未注浆强度对比
通过对图7～图9的分析，进一步确定岩石结构面注浆效果的等效参数模拟方法是可行的.
5 结构面强度与几何形态关系
表4为5个等级法向应力下1个台阶模型和组合结构面模型抗剪强度大小关系.
表4 不同法向应力2种类型结构面抗剪强度比较/MPa结构面类型法向应力1 2 3 4 5 1个台阶 1.821 3.642 5.500 7.291 9.080组合 1.868 3.740 5.609 7.482
9.352增长率/% 2.580 2.690 1.980 2.620 2.990
由表4可以看出，组合结构面模型相对于1个台阶结构面在增加一个锯齿的情况下，结构面抗剪强度并没有明显的提高，说明大爬坡角对结构面抗剪强度大小起控制作用，定义其为控制性爬坡角；小爬坡角对结构面抗剪强度大小不起控制性作用，
定义其为非控制性爬坡角.
由表4中结构面抗剪强度增长率可以看出，组合情况下台阶（即爬坡角为90°）结构面模型对结构面抗剪强度起决定性作用，锯齿对抗剪强度影响较小，平均增长率仅为2.57%.
由此，结合前文分析，可以提出并得到组合情况下岩石结构面抗剪强度修正公式：
式（1）中：ki为非控制性爬坡角因子，以本文中锯齿状结构面模型爬坡角为30°
情况下，其非控制性爬坡角因子ki=0.025 7.
6 结论
在开展室内单轴、剪切试验的基础上，开展了结构面数值模拟试验.通过数值计算，分析了3种类型结构面抗剪力学特性，得出结论如下：
1）在剪切试验中，高宽比较大时，结构面一般发生剪断破坏；高宽比较小时，结构面一般发生压切破坏.
2）结构面抗剪强度随着法向应力的增加近似呈线性增长.
3）通过设置中间层，采用等效参数模拟方法进行结构面注浆模拟试验是可行的.数值模拟结果表明，注浆后结构面抗剪强度得到显著提高，从而在微观上揭示了采用注浆加固岩石结构面的可行性.
4）以数值试验为依据提出了非控制性爬坡角因子ki以及组合结构面抗剪强度修正公式.计算出本文中组合结构面的非控制性爬坡角因子为0.025 7.因此，可以推断
出同一岩体中结构面抗剪强度大小主要由其中较大的爬坡角控制.
参考文献：
[1]杜时贵．岩体结构面的工程性质[M]．北京：地震出版社，1999．
[2]白红杰，陆文，肖正学，等．岩石结构面抗剪强度方向性模拟研究[J]．矿业研
究与开发，2007，10（5）：61－63．
[3]林杭，曹平，周正义．FLAC3D模拟全长注浆锚杆的作用效果[J]．岩土力学，2005，26（增刊 1）：167－170．
[4]程东幸，潘炜，刘大安，等．锚固节理岩体等效力学参数三维离散元模拟[J]．岩土力学，2006，27（12）：2128-2132．
[5]王中文，方建勤．岩体结构面锚固剪切特性的数值分析[J]．煤炭学报，2010，35（5）：729-733．
[6]韩立军，宗义江，韩贵雷，等．岩石结构面注浆加固抗剪特性试验研究[J]．岩
土力学，2011，32（9）：2570-2576．
[7]杨明．岩体结构面力学特性及其锚固效应的数值计算研究[D]．长沙：中南大学，2009．
[8]刘自由，江学良，林杭．软弱结构面加桩特性的数值分析[J]．中南大学学报（自然科学版）：2011，42（5）：1461－1466．
[9]李庆勇，远方，童中明．混凝土剪切破坏数值模拟[J]．低温建筑技术，2011（5）：51－53．
[10]周莲君，彭振斌，何忠明，等．结构面剪切特性的试验与数值模拟分析[J]．科技导报，2009，27（4）：31－35．
[11]李大勇，潘军刚．土的直接剪切试验三维数值模拟研究[J]．山东科技大学学报（自然科学版），2008，27（6）：16-20．
[12]胡黎明，马杰，张丙印．直剪试验中接触面渐进破坏的数值模拟[J]．清华大学学报（自然科学版），2008，48（6）：943-946．
[13]吴乐文．岩石结构面注浆前后力学特性研究[D].赣州：江西理工大学，2012．
[14]杨松林，徐卫亚，朱焕春．锚杆在节理中的加固作用[J]．岩土力学，2002，23（5）：604－607．
[15]金解放，钟海兵，吴越等.静载荷与循环冲击作用下岩石损伤变量定义方法的
选择[J].有色金属科学与工程，2013，4（4）：85-90．。