2017, Vol. 39
随着我国大力发展寒区经济, 寒区的土木工程建设也越来越多.寒区的岩石由于经历低温冻融循环, 岩石的稳定性与耐久性将会遭受很大影响, 从而严重影响了寒区的工程建设.因此, 研究岩石经历不同温度冻融循环后的力学特性具有重大的工程意义.目前, 国内外学者对此问题的研究已经获得了一定的成果.Tan等[1]分析了花岗岩在不同冻融循环次数下的单轴和三轴力学性质, 发现抗压强度、围压和冻融循环次数间的关系可利用莫尔-库仑强度准则进行描述.文献[2]对西班牙Laspra白云岩细观孔隙结构在不同冻融循环次数下的冻融破坏特征和细观缺陷演化机制进行了研究.Mutlutürk等[3]通过研究发现, 岩石在经过循环冻融后, 完整性大受影响, 不同的岩石类型, 岩石的受损速率不同, 提出了岩石经过冻融循环后的损伤数学模型.Pellet等[4]研究了时间效应对岩石蠕变损伤的影响, 通过对P波波速的测量来监测长期蠕变试验过程中岩石的损伤.Bayram[5]对寒区自然岩石经历冻融循环后的强度损失进行了研究, 建立了适合特定石灰石冻融循环后单轴抗压强度的预测方程式.Huang等[6]通过卸荷试验指出, 卸荷过程中初始围压和卸荷速率会影响岩样的破坏和应变能转换模式.周科平等[7]研究岩石在冻融循环作用下微观结构的变化特征, 对冻融循环后的花岗岩进行核磁共振测量, 得到了不同冻融循环次数后岩样的横向弛豫时间T2分布及核磁共振成像图像.袁小清等[8]通过对节理岩体的研究, 提出了冻融细观损伤、受荷细观损伤和节理宏观损伤的概念, 推导出了复合损伤变量, 并且用试验进行验证.吴安杰等[9]对饱水状态下的泥质白云岩进行了不同冻融循环条件下的单轴压缩试验, 建立了以冻融循环次数和应变为控制变量的本构模型.吴刚等[10]将大理岩分成饱水、干燥和对比组, 对其进行冻融循环处理后再进行单轴压缩试验, 归纳出冻融循环后大理岩的主要物理力学特性.刘红岩等[11]通过循环冻融和相似材料试验, 研究了节理岩体在冻融条件下的损伤破坏机制及其相应的力学特性.徐光苗等[12]对红砂岩和页岩进行了单轴和三轴压缩试验, 研究岩石力学特性与温度、岩石含水率的关系.
虽然国内外对于岩石冻融循环后的力学特性的研究已经取得了一定的成果, 但是, 总结前人的成果发现, 大部分学者研究时都将冻融循环的次数作为变量研究, 而忽视了-40 ℃以下温度对岩石力学特性的影响.本文将冻融循环的温度作为变量, 通过试验来研究随着温度降低花岗岩的力学特性的变化.这一问题的研究结果对解决寒区岩石类材料损伤问题具有一定的指导意义.
1 试验概况 1.1 试样制作花岗岩试块取自于福建宁德, 从现场取得一块完整的花岗岩大岩块, 然后将其加工成若干直径50 mm、高度100 mm的圆柱形试块, 加工完以后先将表面有缺陷的试块剔除, 然后通过动弹模仪测出试块的波速, 选择波速相近的试块, 从而减小误差.花岗岩的主要矿物质成分有石英、长石和云母, 测得其密度为2.7 g/cm3.
1.2 试验方案将花岗岩试块分为5组, 对应冻融循环最低温度为-20, -30, -40, -50 ℃和对照组各5块, 共25块.试验前先测出弹性模量, 然后将除对照组以外的4组试块放入水中浸泡48 h以上, 浸泡完毕取出后开始进行冻融循环试验.在高低温试验箱进行冻融循环试验, 设置程序, 最低温度分别为-20, -30, -40, -50 ℃, 冻结4 h, 融化4 h, 循环30次, 待试验结束后取出试块, 测出试块的弹性模量, 然后再对所有试块进行单轴压缩试验.
1.3 试验仪器试验采用的仪器分别为混凝土超声波测试仪、建恒仪器的高低温试验箱WGD501以及西安力创的电液伺服混凝土三轴流变试验机, 如图 1~3所示.
|
图 1 混凝土超声波测试仪 Fig. 1 Ultrasonic concrete tester |
|
图 2 高低温试验箱 Fig. 2 High-low temperature test chamber |
|
图 3 电液伺服混凝土三轴流变试验机 Fig. 3 Electro-hydraulic servo concrete triaxial rheological testing machine |
冻融循环后的花岗岩经过单轴压缩试验后, 主要的破坏模式有两种, 分别为柱状劈裂模式和锥形破坏模式[13], 如图 4和图 5所示.经多个试块破坏模式的分析发现, 不管花岗岩冻融循环的最低温度为多少, 花岗岩单轴压缩破坏时柱状劈裂模式和锥形破坏模式都有可能出现, 因此, 本文认为花岗岩冻融循环后的单轴压缩破坏模式与温度无关.
|
图 4 柱状劈裂模式 Fig. 4 Columnar fracturing mode |
|
图 5 锥形破坏模式 Fig. 5 Tapered damage mode |
通过对混凝土超声波测试仪测出的弹性模量数据的分析可得[13], 当冻结温度为-20 ℃时, 弹性模量降低了0.95 GPa, 相对弹性模量降低了2.4%;当冻结温度为-30 ℃时, 弹性模量降低了5.04 GPa, 相对弹性模量降低了12.53%;当冻结温度为-40 ℃时, 弹性模量降低了7.22 GPa, 相对弹性模量降低了17.95%;当冻结温度为-50 ℃时, 弹性模量降低了7.68 GPa, 相对弹性模量降低了19.09%.由此分析可知, 弹性模量随着冻结温度的降低而呈不同程度的降低, 并且花岗岩弹性模量的劣化速度随着冻结温度的降低先快速降低而后逐渐缓慢降低.弹性模量降低值与温度的关系曲线如图 6所示.ΔE表示弹性模量降低值, T表示冻结温度.
|
图 6 弹性模量降低值与温度的关系曲线图 Fig. 6 Curve of the reduced value of elastic modulus versus the temperature |
冻融循环后花岗岩的单轴压缩应力应变曲线大致可以分为3个阶段[14]:a.孔隙裂隙压密阶段, 岩石中的微裂隙逐渐被压密, 形成非线性变形; b.弹性变形阶段, 应力应变曲线大致呈直线型, 为弹性变形, 并且该阶段曲线的末端为峰值应力, 即该岩石的单轴抗压强度; c.破坏阶段, 岩石达到破坏荷载后, 内部结构完全破坏, 应变也迅速增长.
如图 7所示, σ表示应力, ε表示应变, 由此应力应变曲线图分析可知:a.冻融循环后花岗岩的单轴压缩应力应变曲线形状比较类似, 都为脆性材料应力应变曲线, 而且随着冻融循环温度的降低, 花岗岩的单轴压缩应力应变曲线趋于平缓, 这是由于在相同的循环次数下, 冻结温度越低, 花岗岩的内部裂隙越多, 所以, 压密阶段越长.b.峰值应力和峰值应变分别随着冻结温度的降低而呈线性下降趋势.
|
图 7 冻融循环后花岗岩的单轴应力应变曲线图 Fig. 7 Uniaxial stress-strain curves of granite after freeze-thaw cycles |
在表 1中, σc表示单轴抗压强度, ε表示峰值应变.由表 1可知, 当冻结温度为-20 ℃时, 花岗岩单轴抗压强度降低了46.2 MPa, 相对峰值应力降低了30.8%;当冻结温度为-30，-40，-50 ℃时花岗岩单轴抗压强度分别降低了68.1, 86.8, 102 MPa, 相对峰值应力分别降低了45.4%, 57.7%, 68%.
|
表 1 冻融循环后花岗岩的单轴抗压强度及峰值应变 Table 1 Uniaxial compressive strength and peak strain of granite after freeze-thaw cycles |
经分析可知, 冻结温度越低, 花岗岩的单轴抗压强度降低得越多, 并且大致呈现线性关系, 与文献[12]中饱和红砂岩单轴抗压强度和温度的关系类似.花岗岩的单轴抗压强度与冻结温度的拟合曲线关系如图 8所示, 对拟合曲线作趋势预测, 大致可以预测, 当冻结温度为-60 ℃时的单轴抗压强度约为25 MPa.
|
图 8 花岗岩单轴抗压强度与冻结温度的关系曲线 Fig. 8 Curve of the uniaxial compressive strength versus the freezing temperature of granite |
经分析发现, 花岗岩的冻结温度越低, 单轴压缩峰值应变越大, 并且也是大致呈线性关系, 花岗岩的单轴压缩峰值应变与冻结温度的拟合曲线关系如图 9所示, 对冻结温度为-60 ℃时的峰值应变作了预测, 大约为0.017.
|
图 9 花岗岩峰值应变与冻结温度的关系曲线 Fig. 9 Curve of the peak strain versus the freezing temperature of granite |
为了研究花岗岩的抗冻性, 本文根据文献[15]中给出的岩石冻融系数的公式进行计算.
|
(1) |
式中:Kf为岩石的冻融系数; 
由式(1)计算出花岗岩的冻融系数, 如表 2所示.
|
|
表 2 不同冻结温度的花岗岩的冻融系数 Table 2 Freezing and thawing coefficient of granite at different freezing temperatures |
由表 2的数据分析可知, 花岗岩的抗冻性随着冻结温度的降低而降低, 且花岗岩的冻融系数受冻融循环冻结温度的影响, 呈线性逐渐降低.
2.7 冻融受荷后的损伤扩展力学特性分析张慧梅等[16]提出了冻融损伤、受荷损伤以及总损伤的概念, 并且推导出岩石冻融受荷后的损伤演化方程为
|
(2) |
式中:Dm为岩石冻融受荷总损伤变量; E0为基准损伤状态; En为以冻融循环次数为变量的冻融损伤状态的弹性模量; εf, σf分别为应力峰值和所对应的应变; m为表征材料损伤演化特征的材料参数.
在文献[16]中, 是以冻融循环的次数为变量, 而本文是以冻结温度为变量, 因此, 可将式(2)推广为
|
(3) |
式中, Et为以冻结温度为变量的冻融损伤状态的弹性模量.
将本文试验所得到的数据代入推广的式(3)中, 再将结果绘制成曲线图, 如图 10所示.
|
图 10 损伤变量与应变的关系曲线 Fig. 10 Curves of the damage variables versus the strain |
由图 10可知, 在相同的应变情况下, 花岗岩的总损伤变量随着冻结温度的降低而增大, 当冻结温度为-50 ℃时, 峰值应变值大约是0.015, 总损伤变量约为0.95, 接近于1, 表明花岗岩受冻融循环损伤的影响比较大, 特别是冻结温度的影响; 花岗岩的总损伤变量随着应变的增大而增大, 并且呈现先快速增大而后趋于平缓的趋势.可以预测, 即使温度再低, 峰值应变也不会超过0.017, 而且花岗岩将直接破坏.
3 结论a. 冻融循环后的花岗岩经过单轴压缩试验后, 主要的破坏模式有两种, 分别为柱状劈裂模式和锥形破坏模式, 随着花岗岩冻结温度的降低, 花岗岩单轴压缩的破坏模式不会改变.
b. 冻融循环后花岗岩的单轴压缩应力应变曲线形状比较类似, 都为脆性材料应力应变曲线, 而且随着冻融循环冻结温度的降低, 花岗岩的单轴压缩应力应变曲线也趋于平缓.
c. 冻融循环冻结温度越低, 花岗岩的单轴抗压强度降低得越多, 单轴压缩峰值应变越大.
d. 花岗岩的抗冻性随着冻融循环冻结温度的降低而降低, 并且花岗岩的冻融系数受冻融循环冻结温度的影响, 呈线性逐渐降低.
e. 根据总损伤变量公式得出花岗岩的冻融受荷损伤演化模型, 即损伤变量与应变关系曲线.分析得出, 当应变相同时, 花岗岩总损伤变量随着冻结温度的降低而不断增大, 可以预测, 即使冻结温度不断下降, 应变都存在一个临界值, 一旦超过临界值, 试块将破坏.
| [1] | TAN X J, CHEN W Z, YANG J P, et al. Laboratory investigations on the mechanical properties degradation of granite under freeze-thaw cycles[J]. Cold Regions Science and Technology, 2011, 68(3): 130–138. DOI:10.1016/j.coldregions.2011.05.007 |
| [2] | DE ARGANDOÑ A V G R, REY A R, CELORIO C, et al. Characterization by computed X-ray tomography of the evolution of the pore structure of a dolomite rock during freeze-thaw cyclic tests[J]. Physics and Chemistry of the Earth, Part A:Solid Earth and Geodesy, 1999, 24(7): 633–637. DOI:10.1016/S1464-1895(99)00092-7 |
| [3] | MUTLUTVRK M, ALTINDAG R, TVRK G. A decay function model for the integrity loss of rock when subjected to recurrent cycles of freezing-thawing and heating-cooling[J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(2): 237–244. DOI:10.1016/S1365-1609(03)00095-9 |
| [4] | PELLET F L, FABRE G. Damage evaluation with P-wave velocity measurements during uniaxial compression tests on argillaceous rocks[J]. International Journal of Geomechanics, 2007, 7(6): 431–436. DOI:10.1061/(ASCE)1532-3641(2007)7:6(431) |
| [5] | BAYRAM F. Predicting mechanical strength loss of natural stones after freeze-thaw in cold regions[J]. Cold Regions Science and Technology, 2012, 83-84: 98–102. DOI:10.1016/j.coldregions.2012.07.003 |
| [6] | HUANG D, LI Y R. Conversion of strain energy in Triaxial Unloading Tests on marble[J]. International Journal of Rock Mechanics and Mining Sciences, 2014, 66: 160–168. DOI:10.1016/j.ijrmms.2013.12.001 |
| [7] | 周科平, 李杰林, 许玉娟, 等. 冻融循环条件下岩石核磁共振特性的试验研究[J]. 岩石力学与工程学报, 2012, 31(4): 731–737. |
| [8] | 袁小清, 刘红岩, 刘京平. 冻融荷载耦合作用下节理岩体损伤本构模型[J]. 岩石力学与工程学报, 2015, 34(8): 1602–1611. |
| [9] | 吴安杰, 邓建华, 顾乡, 等. 冻融循环作用下泥质白云岩力学特性及损伤演化规律研究[J]. 岩土力学, 2014, 35(11): 3065–3072. |
| [10] | 吴刚, 何国梁, 张磊, 等. 大理岩循环冻融试验研究[J]. 岩石力学与工程学报, 2006, 25(S1): 2930–2938. |
| [11] | 刘红岩, 刘冶, 邢闯锋, 等. 循环冻融条件下节理岩体损伤破坏试验研究[J]. 岩土力学, 2014, 35(6): 1547–1554. |
| [12] | 徐光苗, 刘泉声, 彭万巍, 等. 低温作用下岩石基本力学性质试验研究[J]. 岩石力学与工程学报, 2006, 25(12): 2502–2508. DOI:10.3321/j.issn:1000-6915.2006.12.016 |
| [13] | 蒋立浩, 陈有亮, 刘明亮. 高低温冻融循环条件下花岗岩力学性能试验研究[J]. 岩土力学, 2011, 32(S2): 319–323. |
| [14] | 刘佑荣, 唐辉明. 岩体力学[M]. 北京: 中国地质大学出版社, 1999: 43-44. |
| [15] | 长江水利委员会. SL 264-2001水利水电工程岩石试验规程[S]. 北京: 中国水利水电出版社, 2001. |
| [16] | 张慧梅, 杨更社. 冻融与荷载耦合作用下岩石损伤模型的研究[J]. 岩石力学与工程学报, 2010, 29(3): 471–476. |