1、侧向变形控制下的灰岩破坏及其峰后本构关系*王东1,2,刘长武1,2,王丁1,2,韩小刚1,2,周晓明1,2(1.四川大学 水力学与山区河流开发保护国家重点实验室,四川 成都 610065;2.四川大学水利水电学院,四川 成都 610065)摘要:利用煤矿现场钻取的灰岩岩芯,进行不同含水状态(干燥和饱和)和不同围压的三轴压缩全程试验、峰后加卸载循环试验,分析干燥和饱水状态岩石的变形破坏规律、围压对岩石峰值强度和残余强度及破裂面角度的影响,研究岩石峰前、峰后强度与变形之间关系。研究表明:不同围压水平引起不同破坏方式,泡水使岩石的脆性增加;残余强度相对峰值强度较小,且随峰后加卸载循环次数的增加继续衰
2、减;该类岩石在低围压下产生拉剪组合破坏,导致其破裂角度较Coulomb 准则预测值偏大,在高围压下发生完全剪切破坏,破裂角度与Coulomb 准则计算结果吻合;岩石轴向应力在峰值前是由其轴向应变决定的,峰后则由侧向变形控制,随侧向变形的增加逐渐衰减至残余强度。在此基础上建立了岩石本构模型,并推导了分段的岩石应力应变函数。利用该模型对岩石应力应变全过程进行计算,计算结果与试验结果吻合较好。关键词:岩石力学;围压;泡水;侧向变形;循环加载;峰后特性中图分类号:TU452 文献标识码:ALimestone failure and post-failure constitutive relation
3、in the control of lateral deformationWANG Dong1,2, LIU Changwu1,2, WANG Ding1,2, HAN Xiaogang1,2,ZHOU Xiaoming1,2 (1.State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, Sichuan 610065,China; 2.College of Water Resource and Hydropower, Sichuan University, C
4、hengdu, Sichuan 610065,ChinaAbstract: The mechanical tests were conducted under different moisture conditions (fully saturated and dry) and different confining pressures using limestone drilled from coal mine, and the cyclic loading tests were carried out after that, to analyze the law of deformatio
5、n and failure, the effect to peak and residual strength and angle of rupture brought by confining pressure, and make research on the law between strength and deformation of whole process. The study results show as follows: different confining pressure leads different form of failure, and pore water
6、will strengthen the brittleness of rock. The residual strength is much smaller than the peak strength, and decreases as the cyclic loading. The form of failure about this classification of rock is tension-shearing under low confining pressure which leads a higher angle of rupture, and fully shearing
7、 under high confining pressure which leads a angle as Coulomb criterion .The axial stress is caused by axial strain before failure, but controled by lateral strain after failure and decreases to the residual strength as the lateral strain increasing. Based on that, rock constitutive model and stress
8、-strain function were founded, which accord with the test results.Key words: rock mechanics; confining pressure; moisture; lateral deformation ;cyclic loading; post-failure behavior 基金项目:国家重点基础研究发展计划(973计划)(项目编号:2010CB226802); 国家自然科学基金资助项目(项目批准号:50879049;50574064)作者简介:王 东(1986-),男,汉族,四川大学水利水电学院博士研究生
9、,主要从事矿山岩石力学、地下工程方面的研究。E-mail:dongdongsky100深部煤炭资源是21世纪我国主体能源的后备力量。随着开采深度的增加,围岩由浅部开采中的弹性状态逐渐进入到峰后的塑性破坏状态,导致工程灾害日趋增多,对深部资源的安全高效开采造成了巨大威胁1。因此,深部资源开采过程中所产生的岩石力学问题已成为国内外研究的焦点2-8。深部开采力学问题就是研究开采围岩破坏后与支护系统相互作用达到二次稳定的复杂力学问题。而目前国内外对于岩石力学的研究,多集中于岩石应力应变的峰前段;损伤力学的引入较好的描述了岩石内部裂纹扩展对岩石强度的影响9-17,但欠缺对岩石宏观裂纹产生后的强度描述。本
10、文取华北煤田分布广泛的大青灰岩,进行不同含水状态(干燥和泡水饱和)的三轴压缩全程试验和峰后加卸载循环试验,结合弹性力学、塑性力学和损伤力学,对岩石的强度、变形和破坏过程进行了研究,以期为矿山深部开采和类似岩土工程提供一定的理论补充。1 试验设备、原理及方法试验所用岩样取自河北省邢台市葛泉煤矿大青灰岩顶板,由于灰岩取样困难,本次试验将取出的8块标准岩样分为两组,第一组做干燥状态下的三轴压缩试验,第二组做常温常压下饱水24h后的三轴压缩试验。利用三轴压缩试验中破坏的岩样进行峰后重复加卸载循环试验。试验采用四川大学水利水电学院的MTS815型电液伺服岩石力学试验系统,对两组岩样进行围压为5MPa、1
11、0MPa、15MPa和20MPa的三轴压缩试验。该套系统是高精度的岩石力学试验设备,可以进行单轴压缩、三轴压缩、孔隙水压及水渗透等试验,并且具有计算机控制、自动数据采集功能,试验机的自动数据采集系统可同步采集的数据包括时间、轴向荷载、试件的轴向变形、环向变形等。整个试验过程由计算机同时按照时间间隔3s、轴向位移0.001mm和环向位移0.005mm的方式自动采集记录应力-应变过程中的轴向压力、轴向位移、环向位移等数据。2 试验结果2.1 三轴压缩试验对灰岩试样(干燥)进行的三轴压缩试验结果如图1所示。在灰岩达到峰值强度前,应力应变曲线近似为一条直线,其弹性模量几乎不变。达到峰值强度后,低围压下
12、试样很快破坏,强度很快降至残余强度,岩石发生脆性破坏,并伴有轴向应变的弹性回弹;在围压较高的情况下,岩样在达到峰值强度后发生了一个塑性屈服过程,其强度逐渐衰减,最后稳定在岩石残余强度,岩石发生延性破坏。岩石的轴向峰值强度随着围压的增加而升高,但增加的幅度逐渐减小;残余强度随围压变化规律与峰值强度类似,但在低围压下增加的幅度更快。图1 岩石三轴压缩应力应变全程曲线Fig.1 Triaxial compression stress-strain curve of rock 2.2 破碎岩石循环加载试验将三轴压缩试验压坏的岩样进行峰后重复加卸载试验,得到残余强度随加卸载循环次数的关系如图2所示,可以
13、看出破坏后的岩石在干燥和饱水状态下,其强度随着加卸载循环的次数增加而逐渐降低,降低程度达15%左右。通过对试验结果的整理,同时得到了岩石抗剪强度参数粘聚力c和内摩擦角从破坏前到破坏后循环加载过程的变化曲线,如图3所示。无论干燥还是饱水状态,岩石破坏前后,抗剪强度指标均大幅下降:干燥状态下灰岩试样的粘聚力从破坏前的8.48Mpa经历破坏后下降到2.9Mpa,约为峰值的30%左右;内摩擦角则从峰值时的48.16下降为破坏后的28,约为峰值的58%。岩样在破坏后随着循环加载次数的增加,其抗剪强度有着进一步的降低,其中粘聚力在经历一次峰后加卸载循环后,下降了17%,从2.9Mpa降为2.4Mpa,二次
14、加卸载循环后降为2.1Mpa,第三次加卸载循环基本不再发生变化,保持在2.1Mpa左右;内摩擦角在岩石破坏后,基本不随加卸载循环而变化,维持在岩石破坏时的28左右。(a)干燥试样(b)饱水试样图2 岩石残余强度与循环加载次数关系曲线Fig.2 Relation between residual strength and cyclic loading 另外,从图3中还可以得出水对岩石峰前峰后抗剪强度的影响:水使灰岩在峰值时的粘聚力下降到干燥状态的50%,而在岩石破坏后与干燥状态下灰岩的峰后粘聚力基本相等;水对于内摩擦角的影响,无论是峰前还是峰后都不显著。 (a)粘聚力变化曲线(b)内摩擦角变化曲
15、线图3 岩石抗剪强度参数变化曲线Fig.3 Variation curve of shear strength parameters2.3 岩石承载能力和破坏形式岩石的极限承载能力与围压的关系如图4所示,可以发现灰岩试样的轴向承载能力随围压增加逐渐增大,且在围压达到一定临界值时,轴向承载能力增加幅度减小。饱水状态与干燥状态下的临界值有所区别:干燥状态下灰岩试样在围压达到10Mpa前,其峰值强度和残余强度随围压增加均快速增大,在围压超过10Mpa以后,由于岩石开始进入塑性变形区,其轴向承载能力随围压变化程度减缓;饱水状态下的灰岩试样在围压达到15Mpa前,其峰值强度和残余强度均随围压增加快速增大
16、,在围压超过15Mpa后开始放缓。可以看出,岩石中的水使岩石脆-延性破坏转化的临界围压强度增加。这个可以从孔隙水压力的角度来解释:令L表示岩石在干燥状态下的脆延性破坏的围压临界值,则在饱水状态下,岩石的脆延性转化临界值L*与孔隙水压力pw之间的关系为: (1)该试验过程不排水,由上式可知饱水灰岩试样在压缩过程中产生了5Mpa左右的孔隙水压力,这导致了泡水岩样进入延性破坏的临界围压强度的增加。孔隙水压力的产生使岩石的脆性增加。(a)干燥试样(b)泡水试样图4 岩石强度与围压关系曲线Fig.4 Relation between rock strength and confining pressur
17、e图5为干燥灰岩试样三轴压缩下的破裂断面示意图,可以看出围压增大,试样的破裂角度(破裂面与最小主应力夹角)呈现减小的趋势。在围压为5Mpa时,破裂角为76,比Coulomb准则预测的(45+/2)=69大了10%左右;而随着围压的增加,破裂角逐渐向69靠拢,与Coulomb准则对破裂面角度的预测相吻合。注:图上无单位数字为破裂尺寸,cm图5 岩石破裂面示意图Fig.5 Schematic diagram of rock fracture surface对于岩石破裂角有多种解释,文献18认为:由于岩石具有明显非均匀性,因而实际破裂角将受到岩样内部的层理、裂隙影响,具有一定随机性。这在解释“低围压
18、下岩石破裂角总是大于理论值”方面有不足之处。笔者认为,由于岩石颗粒应力分布的不均匀性,使得在低围压下岩石内部一部分微元的应力状况为:拉应变先达到极限值,产生拉破坏;而岩石内部其它微元的受力状况为:剪应力先达到极限值,发生了剪切破坏。因而在低围压下可视该类岩样的破坏为“拉剪组合破坏”。拉破坏下岩样近似为劈裂破坏,破裂角度大;剪切破坏下的破裂角可以认为是Coulomb强度准则计算出的(45+/2)。低围压下岩石的破坏是这两种破坏的组合,因此破裂角大于(45+/2)。随着围压的升高,较大的围压对拉破坏产生了较强的抑制作用,因而使岩石微元逐渐进入完全剪切破坏状态,破裂角逼近(45+/2)。3 岩石强度
19、、变形和破坏的讨论3.1岩石强度的决定因素和破坏控制条件图1反映的是岩石压缩过程中轴向应力与轴向应变之间的关系。可以看出在达到峰值强度前,轴向应力应变曲线近似一条直线,可以认为在破坏前满足线弹性广义胡克定律,如下式: (2)式中:1,2,3,1,2,3为岩石三个主应力和主应变,G是拉梅常数。令2=3,并代入泊松比u可以得到: (3) (4)对特定的岩石而言在压缩过程中的泊松比和弹性模量近似常数,可以认为(1-2u3)与轴向应变成线性关系。在围压一定的情况下,轴向应变控制着岩石的轴向应力;岩石的轴向承载能力,就是岩石的轴向变形产生的回弹力。对灰岩等脆性硬岩而言,岩样在破坏前的允许轴向变形越大,意
20、味着其轴向抗压强度越高。但是一旦岩石破坏,轴向应变的增加将不能产生新的回弹力,也就无法产生更大的轴向抗压强度。对于岩石破坏的控制条件,前人研究了很多,多数定义岩石的轴向应变临界值作为岩石破坏的控制条件13,14。从图1中可以看出,岩石在破坏时对应的轴向应变值随围压的变化而变化,取决于应力水平,不是定值。下图6是干燥灰岩试样压缩的轴向应力与环向应变的全程曲线。在不同围压下,试样在峰值强度时对应的环向应变几乎不变,保持在-0.001左右(压为正),说明在不同应力条件下该岩石破坏时的环向应变是基本固定的。图6 岩石轴向应力与环向应变关系曲线Fig.6 Axial Stress-lateral str
21、ain curve of rock同时,比较图1和图6可以发现:图1中的轴向应力-轴向应变曲线,在岩石达到峰值强度后,先是轴向应力产生了跌落,然后再产生的轴向应变的增加;而图6中的轴向应力与环向应变曲线,岩石达到峰值强度后,都是先产生的环向应变的增加(拉应变),然后发生的岩石轴向应力的跌落。也就是说,在岩石达到峰值强度后,先发生环向变形的增加,再产生轴向应力的降低,最后发生轴向应变的增加。这说明,引起岩石轴向应力跌落的原因是环向变形,而不是轴向变形;是环向应变的增加导致岩石的破坏,引起轴向应力的跌落以及轴向应变的增加。因此将岩石的环向应变作为岩石破坏的控制变量是合理而且可行的。由上述得出结论:
22、对于给定的岩石,轴向抗压强度是由岩石的轴向应变(峰前值)决定的;而岩石的破坏受其环向应变控制。如果岩石的环向变形得到控制,使其始终不超过容许值,岩石的轴向抗压强度随轴向应变的增加在理论上可以无限增大。3.2 本构模型的建立曹文贵,方祖烈等提出了“岩石损伤统计软化模型”13,该模型在岩土工程中得到广泛应用: (5) (6) (7)式中: D为岩石的损伤参量,无量纲;m,*分别为Weibull分布的分布标度和以应变形式表征的形态参数,它们反映了岩石的力学性质,岩石的损失就是由这些微元体的不断破坏引起的。为岩石微元体的轴向应变;()为轴向应变为时微元体破坏的破坏概率,是对岩石在加载过程中微元体损坏率
23、的一种量度,它从宏观上反映了试样的损伤程度;S和Sm分别为岩石损伤面积和无损时的全面积。其余符号意义同前述。利用该模型对岩石压缩过程的计算值与试验值比较如下图7所示,可以看出在岩石达到峰值前,该模型对应力应变的预测与实际试验结果较为吻合,但是对峰后段的描述存在不足:无法描述岩石轴向强度在达到峰值后的瞬间跌落;岩石在完全破坏后,损伤参量D将达到1,按此推算,岩石破坏后残余强度将变成2u3,这明显低于岩石实际残余强度值。本文认为岩石在达到峰值强度前,其轴向承载能力是由轴向应变决定的,且随着轴向应变的增加而逐渐增大;当岩石达到峰值强度以后,岩石虽然轴向应变继续增加,其承载能力是下降的,承载能力的降低
24、是由岩石的环向应变增加控制的;岩石在压缩过程中就有环向拉应变的产生,但在峰值前并未达到临界值或者并未达到影响岩石强度的程度,在达到临界值(即峰值强度)时,环向应变开始影响岩石强度,并且峰值前的累积环向应变也瞬间对岩石的强度产生衰减作用。岩石的应力应变本构模型示意图如图8所示。图7 灰岩三轴压缩应力应变过程Fig.7 Triaxial compression stress-strain curve of limestone0岩石没有损伤时的峰值强度;d岩石完全破坏后的残余强度;3c峰值时对应的的环向应变;3d残余强度对应的的环向应变图8 岩石峰后本构模型示意图Fig.8 Schematic di
25、agram of post-failure constitutive model 据此建立岩石的应力应变关系如下: (8) 其中,岩石在没有损伤情况下的峰值强度0用下式计算: (9) 式中:1c为岩石在三轴压缩达到峰值强度时对应的轴向应变值。3.3 模型的解释及验证该模型从岩石压缩破坏的实际情况出发,将应力应变曲线分段表达:峰前段岩石承载能力是岩石轴向应变的函数,由轴向应变控制;峰后段看做岩石强度在环向变形作用下的衰减,只是环向应变的函数,与轴向应变无关,直至最后稳定到残余强度阶段。利用本模型对试验数据进行了计算,理论值与实际值的吻合程度较高。(a)围压5Mpa(b)围压20Mpa图9 本构模
26、型对岩石应力应变全过程模拟结果Fig.9 The computing result of rock stress-strain by the constitutive model 4 结论本文通过对灰岩试样进行三轴压缩全程试验及破坏后的加卸载循环试验,结合弹性力学、塑性力学和损伤力学的相关知识,对岩石的应力应变过程进行了研究,得出以下结论:(1)岩石在三轴压缩过程中,围压水平的不同表现出不同的破坏方式:低围压下呈现出典型的脆性破坏;高围压下表现出延性破坏。岩石轴向承载能力随围压增加而增加,但当围压超过某一临界值后,轴向强度的增加幅度逐渐减小,该临界值可视为岩石脆延性转化的临界围压。泡水导致岩石
27、内部孔隙水压力的存在,加大了岩石脆延性转化的临界围压值,在某种程度上增加了岩石的脆性。(2)岩石的残余强度相对峰值强度下降较大,干燥状态下约为峰值的30%,泡水状态下约为峰值强度的40%;反映岩石抗剪强度指标的粘结力和内摩擦角在破坏时产生不同程度的降低。破碎岩石经历加卸载循环的过程中,粘聚力下降了17%左右;内摩擦角则变化不大。(3)岩石在低围压下的破坏属于“拉剪组合”破坏,这导致了其破裂角大于按照Coulomb准则预测的(45+/2);随着围压的升高,围压对“拉破坏”的抑制作用逐渐增加,使得岩石逐渐进入完全剪切破坏状态,破裂角趋近Coulomb准则的计算结果。(4)岩石的轴向抗压强度是由岩石
28、破坏前所能达到的轴向应变决定的;而岩石的破坏受其环向应变控制。基于此,建立了岩石的应力应变分段函数,峰前段其轴向应力是轴向应变的函数,峰后段则是受岩石的环向应变控制,是环向应变的函数。利用该模型对灰岩试样进行了计算,得出的计算值与试验值较为吻合,说明该模型能较好反映岩石应力应变的过程。参考文献1 何满潮,谢和平,彭苏萍,等.深部开采岩体力学研究J.岩石力学与工程学报,2005,24(16):2803-2813.(HE Manchao, XIE Heping, PENG Suping, JIANG Yaodong. Study on rock mechanics in deep mining e
29、ngineeringJ. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(16):2803-2813. (in Chinese)2 钱七虎.非线性岩石力学的新进展深部岩体力学的若干问题A.见:中国岩石力学与工程学会编.第八次全国岩石力学与工程学术大会论文集C.北京:科学出版社,2004.10-17.(QIAN Qihu. The current development of nonlinear rock mechanics: the mechanics problems of deep rock massA. In: Ch
30、inese Society of Rock Mechanics and Engineering ed. Proceedings of the 8th Rock Mechanics and Engineering ConferenceC.Beijing: Science Press, 2001.10-17.(in Chinese)3 KIDYBINSKI A, DUBINSKI J. Strata Control in Deep Mines M. Rotterdam: A.A.Balkema,1990.4 SUN J, WANG S J. Rock mechanics and rock engi
31、neering in China: developments and current state-of-artJ.International Journal of Rock of Rock Mechanics and Mining Science ,2000,(37):447-465.5 钱鸣高.20年来采场围岩控制理论和实践的回顾J.中国矿业大学学报,2000,19(1):-4.(QIAN Minggao. Review of the theory and practice of strata control around longwall face in recent 20 years J
32、.Journal of Chinese University of Mining and Technology,2000, 19(1):1-4.(in Chinese) 6 MALAN D F, SPOTTISWOODE S M. Time-dependent fracture zone behavior and seismicity surrounding deep level stopping operations A. In : Rock burst and Seismicity in MinestC.Rotterdam: A.A.Balkema, 1997.173-177.7 SWLL
33、ER E J, KLERCK P. Modeling of the effect of discontinuities on the extent of the fracture zonr surrounding deep tunnelsJ. Tunneling and Underground Space Technology,2000,15(4):463-469.8 谢和平.深部高应力下的资源开采现状、基础科学问题与展望A.见:香山科学会议编.科学前沿与未来(第六集)C.北京:中国环境科学出版社,2002.179-191.(XIE Heping. Resources development
34、under high ground stress: present state, base science problems and perspectiveA. In :Xiangshan Science Conference ed. Science Foreland and Future (Volume VI) C. Beijing: China Environment Science Press , 2002. 179-191.(in Chinese)9 谢和平,彭瑞东,鞠杨,等基于断裂力学与损伤力学的岩石强度理论研究进展J自然科学进展,2004, 14(10):10861092(XIE
35、Heping,PENG Ruidong,JU Yang,et alProgress in strength theory of rocks based on fracture mechanics and damage mechanicsJ.Progress in Nature,2004,14(10):1086-1092(in Chinese)10刘齐建,杨林德,曹文贵.岩石统计损伤本构模型及其参数反演J.岩石力学与工程学报,2005,24(4): 616-621. (LIU Qijian,YANG Linde,CAO Wengui. Statistical damage constitutiv
36、e model for rock and back analysis of its parametersJ.Chinese Journal of Rock Mechanics and Engineering,2005,24(4): 616-621. (in Chinese) )11 LEMAITRE J. How to Use Damage MechanicsM.S.1.:Nuclear Engineering and Design, 1984. 12 KRAJCINOVIC D,SILVA M A G. Statistical aspects of the continuous damage
37、 theoryJ.Int. J. Solids Structures,1982, 18(7): 551-56213曹文贵,方祖烈,唐学军.岩石损伤软化统计本构模型之研究J.岩石力学与工程学报,1998,17(6):628-633.(CAO Wengui, Fang Zulie, Tang Xuejun. A study of statistical constitutive model for soft and damage rockJ. Chinese Journal of Rock Mechanics and Engineering, 1998, 17(6):628-633.(in Chi
38、nese)14徐卫亚,韦立德.岩石损伤软化统计本构模型的研究J.岩石力学与工程学报,2002,21(6):787-791.(XU Weiya,WEI Lide. Study on statistical damage constitutive model of rockJ.Chinese Journal of Rock Mechanics and Engineering, 2002, 21(6):787-791.(in Chinese))15谢和平.岩石、混凝土损伤力学M.徐州: 中国矿业大学出版社, 1990.(XIE Heping. Damage Mechanics in Rock and
39、 ConcreteM. Xuzhou: China University of Mining and Technology Press,1990.(in Chiese))16唐春安.岩石破裂过程中的灾变M.北京:煤炭工业出版社,1993. (TANG Chunan. Catastrophe in Rock Unstable FailureM. Beijing:China Coal Industry Publishing House, 1993.(in Chinese)17余天庆,钱济成.损伤理论及其应用M.北京:国防工业出版社,1993.10(YU Tianqing, QIAN Jicheng
40、. Damage Theory and ApplicationM.Beijing: National Defense Industry Press,1993,10. (in Chinese)18尤明庆.岩石三轴压缩的破坏形式和Coulomb强度准则J. 地质力学学报,2002,8(2):179-185.(YOU Mingqing. Destroy character and coulomb criterion of rock specimen in pswudotriaxial compressionJ. Journal of Geomechanics. 2002,8(2):179-185.(in Chinese)作者简介:王东(1986.1-),男,汉族,四川大学水利水电学院博士研究生,主要研究方向为岩土工程、地下工程。手机:13402801920 办公室电话:02885463139邮箱:dongdongsky100通信地址:四川省成都市四川大学望江校区水利水电学院A座402邮编:610065