高围压下砂岩循环加-卸载损伤本构及损伤阈值

刘之喜

引用本文:
Citation:

高围压下砂岩循环加-卸载损伤本构及损伤阈值

    作者简介: 刘之喜(1992—),男,硕士,主要从事岩石非线性力学及能量演化研究. E-mail:17355486579@163.com;
  • 中图分类号: TU457; O347.3

Damage Constitutive of Cyclic Loading and Damage Threshold of Rock under High Confining Pressure

  • CLC number: TU457; O347.3

  • 摘要: 在深部巷道掘进过程中,岩体受周期荷载和高地应力共同作用,一般的强度准则无法对岩石的应力-应变曲线进行描述和表征,而高地应力下岩石的循环加-卸载本构是预测深部巷道在周期荷载作用下长期稳定性的关键,因此关于高围压下岩石的本构研究具有一定的工程意义。岩石受外部荷载直至破坏的过程是旧裂纹扩展和新裂纹增长、扩展的过程,已有的研究表明,单位体积岩石裂纹数服从Weibull分布,Griffith准则认为岩石的破坏是由岩石的裂纹扩展所导致的,为此通过对Weibull建立的统计损伤本构扩展,基于裂纹数服从Weibull分布以及裂纹扩展服从Griffith准则,建立了高围压下岩石循环加-卸载统计损伤本构。通过本构的等式变换,对本构损伤阈值进行了研究,并进一步以砂岩的偶数次循环加-卸载试验数据对本构参数进行拟合,得到参数的演化规律,将本构理论值与砂岩奇数次循环加-卸载试验数据对比,验证了本构模型的正确性,为高围压下岩石的循环加-卸载本构研究提供了新的思路。
  • 图 1  不同围压下${{\textit{σ}} _1}$ε1曲线

    Figure 1.  Axial stress-strain curves under different confining pressures

    图 2  第1次循环加载曲线拟合及验证

    Figure 2.  Curve fitting and verification of the first cyclic loading

    图 3  循环加-卸载本构参数演化规律

    Figure 3.  Evolution of cyclic loading and unloading constitutive parameters

    图 4  围压与本构模型参数的关系

    Figure 4.  Relationship between confining pressure and constitutive model parameters

    图 5  奇数次循环的加-卸载曲线理论值与原始数据对比

    Figure 5.  Comparison of theoretical values of loading and unloading curves with original data for odd cycles

    表 1  围压为60 MPa时本构模型参数

    Table 1.  Parameters of constitutive model with confining pressure of 60 MPa

    MLoad parametersUnload parameters
    F0mF0m
    25.870 04–0.202 9821.330 95–0.141 14
    46.648 45–0.216 0323.502 82–0.137 86
    67.103 79–0.222 7724.513 79–0.136 14
    87.426 85–0.229 5125.231 08–0.134 94
    107.677 44–0.231 4026.171 96–0.133 82
    下载: 导出CSV
  • [1] 王伟, 田振元, 朱其志, 等. 考虑残余应力的砂岩损伤理论模型 [J]. 岩石力学与工程学报, 2015, 34(Suppl 2): 3676–3682.
    WANG W, TIAN Z Y, ZHU Q Z, et al. Study of statistical damage onstitutive model for rock considering waterpressure [J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(Suppl 2): 3676–3682.
    [2] 刘涛, 杨鹏, 吕文生, 等. 岩石在不同应力幅值下受低频循环扰动的力学特性试验 [J]. 煤炭学报, 2017, 40(9): 2280–2286.
    LIU T, YANG P, LÜ W S, et al. Rock mechanical properties experiments with low-frequency circulation disturbance under different stress amplitudes [J]. Journal of China Coal Society, 2017, 40(9): 2280–2286.
    [3] 何明明, 陈蕴生, 李宁, 等. 单轴循环荷载作用下砂岩变形特性与能量特征 [J]. 煤炭学报, 2015, 40(8): 1805–1812.
    HE M M, CHEN Y S, LI N, et al. Deformation and energy characteristics of sandstone subjected to uniaxial cyclic loading [J]. Journal of China Coal Society, 2015, 40(8): 1805–1812.
    [4] 姚强岭, 刘亚鹏, 陈田, 等. 地下水库人工坝体强度损伤演化特征试验研究 [J]. 煤炭学报, 2018, 43(4): 1111–1117.
    YAO Q L, LIU Y P, CHEN T, et al. Experimental study of damage evolution of aritifical dam strength of underground reservoir [J]. Journal of China Coal Society, 2018, 43(4): 1111–1117.
    [5] 谢和平, 高峰, 鞠杨, 等. 深部开采的定量界定与分析 [J]. 煤炭学报, 2015, 40(1): 1–10.
    XIE H P, GAO F, JU Y, et al. Quantitative definition and investigation of deep mining [J]. Journal of China Coal Society, 2015, 40(1): 1–10.
    [6] DEHLER W, LABUZ J F. Stress path testing of an anisotropic sandstone [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2007, 133(1): 116–119. doi: 10.1061/(ASCE)1090-0241(2007)133:1(116)
    [7] 王者超, 赵建纲, 李术才, 等. 循环荷载作用下花岗岩疲劳力学性质及其本构模型 [J]. 岩石力学与工程学报, 2012, 31(9): 1888–1900. doi: 10.3969/j.issn.1000-6915.2012.09.021
    WANG Z C, ZHAO J G, LI S C, et al. Fatigue mechanical behavior of granite subject to cyclic load and its constitutive model [J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(9): 1888–1900. doi: 10.3969/j.issn.1000-6915.2012.09.021
    [8] 易其康, 马林建, 刘新宇, 等. 考虑频率影响的盐岩变参数蠕变损伤模型 [J]. 煤炭学报, 2015, 40(Suppl 1): 93–99.
    YI Q K, MA L J, LIU X Y, et al. A variable parameters-based creep damage model considering the effect of frequency for rock salt [J]. Journal of China Coal Society, 2015, 40(Suppl 1): 93–99.
    [9] 霍润科, 李宁, 刘汉东. 受酸腐蚀砂岩的统计本构模型 [J]. 岩石力学与工程学报, 2005, 24(11): 1852–1856. doi: 10.3321/j.issn:1000-6915.2005.11.006
    HUO R K, LI N, LIU H D. Statistical constitutive model of sandsyone subjected to hydrochloric acid attack [J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(11): 1852–1856. doi: 10.3321/j.issn:1000-6915.2005.11.006
    [10] 李西蒙, 刘长友, 鲁岩, 等. 单轴分级循环加载条件下砂岩疲劳变形特性与损伤模型研究 [J]. 中国矿业大学学报, 2017, 46(1): 8–16.
    LI X M, LIU C Y, LU Y, et al. Fatigue deformation characteristics and damage model of sandston subjected to uniaxial step cyclic loading [J]. Journal of China University of Mining and Technology, 2017, 46(1): 8–16.
    [11] 曹文贵, 张升. 基于Mohr-Coulomb准则的岩石损伤统计分析方法研究 [J]. 湖南大学学报(自然科学版), 2005, 32(1): 43–47.
    CAO W G, ZHANG S. Study on the statistical analysis of rock damage based on Mohr-Coulomb criterion [J]. Journal of Hunan University (Natural Sciences), 2005, 32(1): 43–47.
    [12] 张平阳, 夏才初, 周舒威, 等. 循环加-卸载岩石本构模型研究 [J]. 岩土力学, 2015, 36(12): 3354–3359.
    ZHANG P Y, XIA C C, ZHOU S W, et al. Constitutive model for rock under cyclic loading and unloading [J]. Rock and Soil Mechanics, 2015, 36(12): 3354–3359.
    [13] 邓华锋, 胡玉, 李建林, 等. 循环加卸载过程中砂岩能量耗散演化规律 [J]. 岩石力学与工程学报, 2016, 35(Suppl 1): 2869–2875.
    DENG H F, HU Y, LI J L, et al. The evolution of sandstone energy dissipation under cyclic loading and unloading [J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(Suppl 1): 2869–2875.
    [14] 刘树新, 刘长武, 韩小刚, 等. 基于损伤多重分形特征的岩石强度Weibull参数研究 [J]. 岩土工程学报, 2011, 33(11): 1786–1791.
    LIU S X,LIU C W, HAN X G, et al. Weibull distribution parameters of rock strength based on multi-fractal characteristics of rock damage [J]. Chinese Journal of Geotechnical Engineering, 2011, 33(11): 1786–1791.
    [15] 张媛, 许江, 杨红伟, 等. 循环荷载作用下围压对砂岩滞回环演化规律的影响 [J]. 岩石力学与工程学报, 2011, 30(2): 320–326.
    ZHANG Y, XU J, YANG H W, et al. Effect of confining pressure on evolution law of hysteresis loop of sandstone under cyclic loading [J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(2): 320–326.
    [16] 赵永川, 杨天鸿, 肖福坤, 等. 西部弱胶结砂岩循环载荷作用下塑性应变能变化规律 [J]. 煤炭学报, 2015, 40(8): 1813–1819.
    ZHAO Y C, YANG T H, XIAO F K, et al. The variation law of plastic strain energy of western weak cemented sandstone during cyclic loading experiment [J]. Journal of China Coal Society, 2015, 40(8): 1813–1819.
    [17] 何俊, 潘结南, 王安虎. 三轴循环加卸载作用下煤样的声发射特征 [J]. 煤炭学报, 2014, 39(1): 84–90.
    HE J, PAN J N, WANG A H. Acoustic emission characteristics of coal specimen under triaxial cyclic loading and unloading [J]. Journal of China Coal Society, 2014, 39(1): 84–90.
    [18] 丁长栋, 张杨, 杨向同, 等. 致密砂岩高围压和高孔隙水压下渗透率演化规律及微观机理研究 [J]. 岩土力学, 2019, 40(9): 15–23.
    DING C D, ZHANG Y, YANG X T, et al. Permeability evolution of tight sandstone under high confining pressure and high pore pressure and its microscopic mechanism [J]. Rock and Soil Mechanics, 2019, 40(9): 15–23.
  • [1] 罗吉安刘丰茂刘之喜马雷鸣陈烨开李欣慰 . 基于八面体理论的岩石循环加-卸载本构模型及修正. 高压物理学报, 2020, 34(2): 024202-1-024202-9. doi: 10.11858/gywlxb.20190797
    [2] 肖大武李英雷胡时胜 . 高温高应变率下纯锆的本构模型研究. 高压物理学报, 2009, 23(1): 46-50 . doi: 10.11858/gywlxb.2009.01.008
    [3] 李磊张先锋吴雪高飞刘闯 . 不同硬度30CrMnSiNi2A钢的动态本构与损伤参数. 高压物理学报, 2017, 31(3): 239-248. doi: 10.11858/gywlxb.2017.03.005
    [4] 李阳龙王伟平骆永全王海峰张大勇 . 1 064 nm激光对氧化铟锡薄膜的损伤研究. 高压物理学报, 2012, 26(1): 107-112. doi: 10.11858/gywlxb.2012.01.016
    [5] 李春宏周信达韩伟冯斌朱启华郑万国 . 纳秒激光辐照下石英基单层石墨烯的损伤特性. 高压物理学报, 2016, 30(5): 392-398. doi: 10.11858/gywlxb.2016.05.008
    [6] 王政倪玉山曹菊珍张文 . 冲击载荷下混凝土本构模型构建研究. 高压物理学报, 2006, 20(4): 337-344 . doi: 10.11858/gywlxb.2006.04.001
    [7] 王立闻庞宝君陈勇张凯 . 高温处理后活性粉末混凝土动力学行为及本构模型研究. 高压物理学报, 2012, 26(4): 361-368. doi: 10.11858/gywlxb.2012.04.001
    [8] 李雪艳李志斌张舵 . 考虑温度效应的泡沫铝准静态压缩本构模型. 高压物理学报, 2018, 32(4): 044103-1-044103-6. doi: 10.11858/gywlxb.20170642
    [9] 陈青山苗应刚郭亚洲李玉龙 . 比较93钨合金材料的3种本构模型. 高压物理学报, 2017, 31(6): 753-760. doi: 10.11858/gywlxb.2017.06.010
    [10] 邹学韬张晓晴姚小虎 . 压剪载荷作用下TB6钛合金的动态力学性能. 高压物理学报, 2019, 33(2): 024206-1-024206-9. doi: 10.11858/gywlxb.20190713
    [11] 李定远朱志武卢也森 . 冲击加载下42CrMo钢的动态力学性能及其本构关系. 高压物理学报, 2017, 31(6): 761-768. doi: 10.11858/gywlxb.2017.06.011
    [12] 王涛汪兵林健宇柏劲松李平钟敏陶钢 . 金属锡Rayleigh-Taylor不稳定性对模型参数敏感性的数值分析. 高压物理学报, 2020, 34(2): 022301-1-022301-9. doi: 10.11858/gywlxb.20190813
    [13] 雷忠琦姚小虎龙舒畅常建虎王海波 . 聚苯乙烯泡沫材料的动态压缩特性. 高压物理学报, 2019, 33(2): 024202-1-024202-10. doi: 10.11858/gywlxb.20180655
    [14] 李俊玲卢芳云傅华赵玉刚谭多望 . 某PBX炸药的动态力学性能研究. 高压物理学报, 2011, 25(2): 159-164 . doi: 10.11858/gywlxb.2011.02.012
    [15] 姚小虎任会兰林荣张晓晴 . 聚合物泡沫材料动态力学性能及其能量吸收研究. 高压物理学报, 2012, 26(5): 531-536. doi: 10.11858/gywlxb.2012.05.008
    [16] 贾宇刘彦梁晓璐郑腾 . 高强度合金钢30CrMnMoRE/30CrMnSi的动态力学性能. 高压物理学报, 2018, 32(4): 044102-1-044102-7. doi: 10.11858/gywlxb.20170659
    [17] 王政楼建锋勇珩梁龙河 . 岩石、混凝土和土抗侵彻能力数值计算与分析. 高压物理学报, 2010, 24(3): 175-180 . doi: 10.11858/gywlxb.2010.03.003
    [18] 高宁朱志武 . 5083铝合金宽应变率实验与基于损伤的本构模型研究. 高压物理学报, 2017, 31(1): 51-60. doi: 10.11858/gywlxb.2017.01.008
    [19] 张延耿楼建锋周婷婷洪滔张树道 . PBX炸药含各向异性损伤的黏弹性统计微裂纹本构模型初步研究. 高压物理学报, 2016, 30(4): 301-310. doi: 10.11858/gywlxb.2016.04.006
    [20] 张若棋丁育青汤文辉冉宪文 . 混凝土HJC、RHT本构模型的失效强度参数. 高压物理学报, 2011, 25(1): 15-22 . doi: 10.11858/gywlxb.2011.01.003
  • 加载中
图(5)表(1)
计量
  • 文章访问数:  252
  • 阅读全文浏览量:  214
  • PDF下载量:  7
出版历程
  • 收稿日期:  2019-07-11
  • 录用日期:  2019-07-26
  • 刊出日期:  2020-08-05

高围压下砂岩循环加-卸载损伤本构及损伤阈值

    作者简介:刘之喜(1992—),男,硕士,主要从事岩石非线性力学及能量演化研究. E-mail:17355486579@163.com
  • 安徽理工大学土木建筑学院,安徽 淮南 232001

摘要: 在深部巷道掘进过程中,岩体受周期荷载和高地应力共同作用,一般的强度准则无法对岩石的应力-应变曲线进行描述和表征,而高地应力下岩石的循环加-卸载本构是预测深部巷道在周期荷载作用下长期稳定性的关键,因此关于高围压下岩石的本构研究具有一定的工程意义。岩石受外部荷载直至破坏的过程是旧裂纹扩展和新裂纹增长、扩展的过程,已有的研究表明,单位体积岩石裂纹数服从Weibull分布,Griffith准则认为岩石的破坏是由岩石的裂纹扩展所导致的,为此通过对Weibull建立的统计损伤本构扩展,基于裂纹数服从Weibull分布以及裂纹扩展服从Griffith准则,建立了高围压下岩石循环加-卸载统计损伤本构。通过本构的等式变换,对本构损伤阈值进行了研究,并进一步以砂岩的偶数次循环加-卸载试验数据对本构参数进行拟合,得到参数的演化规律,将本构理论值与砂岩奇数次循环加-卸载试验数据对比,验证了本构模型的正确性,为高围压下岩石的循环加-卸载本构研究提供了新的思路。

English Abstract

  • 地下岩石工程经常受到加-卸载作用,如交通、地震、巷道掘进等,为此诸多学者对岩石的加-卸载性质展开研究,然而这些研究极少涉及深部高地应力下岩石的加-卸载本构。深部岩体开挖作为高围压下岩石加-卸载作用的经典例子,在其掘进期间,围岩受到高围压及加-卸载的反复作用,其动力响应随之改变,对围岩长期稳定性的影响颇大,目前关于高围压下岩石循环加-卸载本构模型的研究较少。许多专家主要针对深部开采的定量界定与分析、低围压岩石本构模型建立、循环加-卸载试验方法、疲劳变形特性、能量耗散等方面进行研究,取得了丰硕的成果[1-6]

    目前岩石循环加载损伤本构模型的建立主要有两种方法:一是通过加-卸载试验,假设岩石服从某些数学模型,从而建立岩石的数学本构模型;二是在岩石中引入损伤变量,为可靠地利用统计物理的概念奠定基础,这方面的开创性工作为Weibull的统计尺寸效应理论,即假设岩石的微元破坏服从Weibull公式,微元强度服从某一强度准则,从而建立岩石的统计损伤本构模型[7-9]。李西蒙等[10]在单轴分级加-载作用的研究中建立了循环加-卸载条件下轴向应变与循环次数之间的理论模型,并且通过逆函数实现了对岩石疲劳寿命的预测。曹文贵等[11]在岩石本构模型研究中基于Weibull公式对岩石的损伤模型进行了研究,将概率分布公式应用到岩石损伤力学中,通过大量的试验验证了其合理性与正确性,并且在工程实际中得到了应用。张平阳等[12]通过在循环加-卸载岩石本构模型中引入内疲劳的变形模量公式,提出了一种新的能够描述岩石循环加-卸载的本构模型。邓华锋等[13]发现能量耗散率与残余应变的相关性,为循环加-卸载中岩石能量参数的准确计算提供了较好的思路。刘树新等[14]基于岩石裂纹三重分布的特征,引入不均匀系数,对Weibull分布下的岩石参数确定方法进行了研究,明确了参数的物理意义和确定方法。以上研究都是关于单轴循环加-卸载或低围压下循环加-卸载研究,高围压下的研究则主要针对岩石渗水特性、能量演化、破坏机制以及声发射等[15-18],较少涉及高围压下岩石循环加-卸载本构模型。对于高地应力和周期荷载作用下岩石表现的非线性力学特性,使用传统的强度准则已无法进行描述和表征,因此高围压下岩石循环加-卸载损伤本构模型研究具有一定的工程意义。

    目前关于岩石的Weibull分布的本构模型,都是对单轴压缩下或常规三轴压缩下的本构模型描述,并且均假设岩石微元破坏服从某一强度准则,但是以最大正应力理论、最大正应变理论、最大剪应力理论、八面体剪应力理论、D-P准则以及摩尔库仑准则建立的Weibull分布,均假设岩石为连续的均匀介质,这些强度准则的建立都与岩石的实际情况不相符。Griffith认为岩石内部存在着许多细微裂隙,在外部荷载作用下细微裂隙产生应力集中现象,致使裂纹扩展引起破坏,因此Griffith准则建立在普遍适用的能量概念和含有裂纹的物体之上,突破了传统的局部分析方法。本研究通过对Weibull建立的统计损伤本构进行扩展,基于裂纹数服从Weibull分布以及裂纹扩展服从Griffith强度准则,建立一种能够表示高围压下岩石循环加-卸载统计损伤本构模型;对本构的损伤阈值进行研究,基于砂岩在高围压下偶数次循环加-卸载曲线,对岩石的参数演化规律进行分析,进一步利用奇数次循环的应力-应变曲线与本构理论值进行对比,进而验证本构模型的正确性与精确性。

    • 由曹文贵等[11]关于岩石统计损伤理论的研究,岩石的损伤本构关系为

      ${\sigma _i} = E{\varepsilon _i}\left( {1 - D} \right) + \mu \left( {{\sigma _j} + {\sigma _k}} \right)$

      式中:$E$为材料弹性模量;${\sigma _i}$为名义应力,${\sigma _i}$${\sigma _j}$${\sigma _k}$为主应力(i, j, k = 1, 2, 3);εi为应变;$\mu $为泊松比;D为损伤变量,对于微裂纹各向同性分布的情况,损伤变量可采用标量形式;对于微裂纹有规律地平面分布的情况,可用与裂纹垂直的矢量表示损伤;对于微裂纹各向异性分布的情况,损伤变量可采用张量形式。以Griffith准则描述岩石裂纹扩展强度,说明裂纹服从各向异性分布,因此损伤变量可以是张量。目前在岩石损伤本构模型研究中主要考虑最大主应力方向的应力-应变关系,并且用张量表示损伤时,虽然能够更真实地反映微观裂纹的排列状态及其力学特性,但是其数学表达式比较复杂,在工程应用方面存在较大难度,所以本研究中的损伤变量使用标量形式。

      研究表明,单位体积岩石裂纹数服从随机分布[9]。根据Griffith准则,当裂纹拉应力超过该点抗拉强度时,裂纹就会扩展,导致试件破坏。岩石加载过程往往是新裂纹生成、扩展和旧裂纹扩展的过程。基于上述研究,令裂纹扩展服从Griffith准则,设加载过程中裂纹扩展数为n,破坏时裂纹扩展总数为N,则可以把损伤变量D定义为nN的比值,n=0时岩石无损伤产生,n=N时裂纹全部张开,岩石破坏。因此做出以下假设:(1)材料内部存在众多互不影响的裂纹,忽略材料局部特性的变化;(2)裂纹三维空间形状和裂缝平面内应力的影响忽略不计;(3)从初始加载到损伤阈值处都是弹性阶段。

      根据Weibull分布,损伤变量D的表达式[16]

      $D = \int_0^F {P\left( x \right){\rm d}x} $

      式中:P(x)表示岩石损伤变量的密度函数,$P\left( F \right) = \dfrac{m}{{{F_0}}}{\left( {\dfrac{{{F}}}{{{{{F}}_{\rm{0}}}}}} \right)^{{m - 1}}}{\rm{exp}}{\left[ {- \left(\dfrac{{{F}}}{{{{{F}}_{\rm{0}}}}}\right)} \right]^m}$,其中F0m为本构模型参数。

      在高围压下岩石的最大主应力与最小主应力均为正值,所以${\sigma _1} + 3{\sigma _3} > 0$,其中:${\sigma _1} $为最大主应力;${\sigma _3} $为最小主应力,与围压${\sigma _{\rm{m}}}$相等,其中围压是不变的。因此在高围压下Griffith破坏准则为

      ${\sigma \rm_t} = \frac{{{{({\sigma _1} - {\sigma _3})}^2}}}{{8\left( {{\sigma _1} + {\sigma _3}} \right)}}$

      式中:${\sigma _{\rm{t}}}$为应力集中处裂纹扩展的拉应力。故在Griffith准则下岩石的裂纹张开强度为

      $F = \frac{{{{\left( {{{\bar \sigma }_1} - {{\bar \sigma }_3}} \right)}^2}}}{{8\left( {{{\bar \sigma }_1} + {{\bar \sigma }_3}} \right)}}$

      式中:${{\bar \sigma }_1}$${{\bar \sigma }_3}$为损伤变量修正后的名义最大主应力和名义最小主应力。其中

      ${\bar \sigma _i} = \frac{{{\sigma _i}}}{{1 - D}}$

      将式(4)代入式(2),可得

      $D = 1 - {\rm{exp}}\left\{ {{\rm{ - }}{{\left[ {\frac{{{{\left( {{{\bar \sigma }_1}{\rm{ - }}{{\bar \sigma }_3}} \right)}^2}}}{{8{F_0}\left( {{{\bar \sigma }_1} + {{\bar \sigma }_3}} \right)}}} \right]}^{{m}}}} \right\}$

      将式(6)再代入式(1),可得其本构模型的表达式为

      ${\sigma _1} = E{\varepsilon _1}{\rm{exp}}\left\{ {{\rm{ - }}{{\left[ {\frac{{{{\left( {{{\bar \sigma }_1}{\rm{ - }}{{\bar \sigma }_3}} \right)}^2}}}{{8{F_0}\left( {{{\bar \sigma }_1} + {{\bar \sigma }_3}} \right)}}} \right]}^{m}}} \right\} + 2\mu {\sigma _3}$

      Weibull分布的概率公式的损伤变量表达式中存在损伤变量,使得本构模型无法使用数据处理软件进行拟合,本研究将基于Weibull分布的本构模型进行简单变换,消去损伤变量表达式中的损伤变量,并对本构进行等式变换。

      将式(7)变换可得

      $\frac{{{\sigma _1} - 2\mu {\sigma _3}}}{{E{\varepsilon _1}}} = {\rm{exp}}\left\{ { - {{\left[ {\frac{{{{\left( {{{\bar \sigma }_1}{\rm{ - }}{{\bar \sigma }_3}} \right)}^2}}}{{8{F_0}\left( {{{\bar \sigma }_1} + {{\bar \sigma }_3}} \right)}}} \right]}^m}} \right\} = 1{\rm{ - }}D$

      根据式(5),将式(6)由含损伤变量的名义主应力表达式转换成由主应力和损伤变量共同描述的岩石裂纹张开强度

      $F = \frac{{{{({\sigma _1} - {\sigma _3})}^2}}}{{8({\sigma _1} + {\sigma _3})\left( {1{\rm{ - }}D} \right)}}$

      将式(8)代入式(9),消去损伤变量,可得

      $F = \frac{{{{\left( {{\sigma _1} - {\sigma _3}} \right)}^2}E{\varepsilon _1}}}{{8\left( {{\sigma _1} + {\sigma _3}} \right)\left( {{\sigma _1} - 2\mu {\sigma _3}} \right)}}$

      将式(10)代入本构方程式(7),可得

      ${\sigma _1} = E{\varepsilon _1}{\rm{exp}}\left\{ {{\rm{ - }}{{\left[ {\frac{{{{({\sigma _1}{\rm{ - }}{\sigma _3})}^2}E{\varepsilon _1}}}{{8{F_0}\left( {{\sigma _1} + {\sigma _3}} \right)\left( {{\sigma _1} - 2\mu {\sigma _3}} \right)}}} \right]}^{m}}} \right\} + 2\mu {\sigma _3}$

    • 图1为张媛等[15]获得的高围压下岩石循环加-卸载试验曲线。张媛等[15]进行了围压分别为20、40、60 MPa时的三轴循环加-卸载试验,每次加载的围压都相同,轴向荷载采用应力控制方式,加载波形为三角波。每次循环的加载上限为岩石强度的70%,下限为岩石强度的25%。砂岩的弹性模量为45 GPa,泊松比为0.142。从其使用数学模型得到的拟合曲线中,可以获得每个循环中岩石的应力-应变曲线。本研究使用高围压循环加-卸载岩石本构模型,对张媛等[15]获得的试验数据进行拟合。

      图  1  不同围压下${{\textit{σ}} _1}$ε1曲线

      Figure 1.  Axial stress-strain curves under different confining pressures

    • 为了研究本构模型的适用范围,需要根据高围压下循环加-卸载应力-应变曲线和本构模型,对岩石的损伤变量进行探讨,以此确定岩石的损伤阈值,确定本构模型的使用范围。将式(1)进行等式变换,可得

      $D = \frac{{E{\varepsilon _i} + \mu \left( {{\sigma _j} + {\sigma _k}} \right) - {\sigma _i}}}{{E{\varepsilon _i}}}$

      在围压相等的高围压下,式(12)可变为

      $D = \frac{{E{\varepsilon _1} + 2\mu {\sigma _3} - {\sigma _1}}}{{E{\varepsilon _1}}}$

      D=0时,岩石不发生损伤;当D > 0时,岩石开始发生损伤,并且损伤是一个连续过程。将岩石参数代入式(13),可知砂岩在围压为20、40、60 MPa时发生损伤,应变值分别为3.182 2×10–4、6.364 4×10–4、7.111 1×10–4。由图1可知:在应力较小处,岩石已经产生损伤,并且不同围压下循环下限的应变值均大于产生初始损伤时的应变阈值;随着围压增大,岩石的损伤阈值也随之增加,呈非线性增长。高围压下岩石本构模型产生损伤时的轴向应变值较小,即岩石的弹性阶段较短,从侧面反映了岩石在高围压下的非线性力学特征。

    • 通过推导得到了高围压下砂岩的本构模型,如式(11)所示。式(11)中包含$\exp \left[ { - {{\left( {F/F_0} \right)}^m}} \right]$ 形式,其中F0m需要通过拟合得到,且两个参数并不是相互独立的。一般采用Levenberg Marquardt算法对式(11)进行迭代。在使用Levenberg Marquardt算法进行迭代的过程中,如果对参数的拟合初值选取不当,容易造成本构无法拟合,为此对式(11)进行等式变换,变换后的形式为

      $y = {F_0}{x^{\textstyle\frac{1}{m}}}$

      式中:$y = \dfrac{{{{\left( {{\sigma _1} - {\sigma _3}} \right)}^2}E{\varepsilon _1}}}{{8\left( {{\sigma _1} + {\sigma _3}} \right)\left( {{\sigma _1} - 2\mu {\sigma _3}} \right)}}$$x = \ln \dfrac{{E{\varepsilon _1}}}{{ {{\sigma _1} - 2\mu {\sigma _3}} }}$F0m均为本构模型参数。与式(11)相比,式(14)在使用Levenberg Marquardt算法拟合的过程中,可以直接将两个拟合参数初值赋1进行迭代拟合,解决了拟合参数初值不易选取的问题。以下将以第1次加载阶段的参数拟合为例说明具体拟合过程。

      图2为岩石在60 MPa围压下第1次加载时的应力-应变拟合曲线。从图2可以看出,第1次加载的拟合效果相对良好,并且将拟合得到的参数代入本构模型中时,能够较好地描述岩石的应力-应变曲线。

      图  2  第1次循环加载曲线拟合及验证

      Figure 2.  Curve fitting and verification of the first cyclic loading

      本研究以围压为60 MPa时砂岩第2、4、6、8、10次循环加-卸载原始数据为例进行拟合,具体参数如表1所示,其中M为循环次数。拟合的校正系数均在0.99以上,说明拟合效果良好。

      MLoad parametersUnload parameters
      F0mF0m
      25.870 04–0.202 9821.330 95–0.141 14
      46.648 45–0.216 0323.502 82–0.137 86
      67.103 79–0.222 7724.513 79–0.136 14
      87.426 85–0.229 5125.231 08–0.134 94
      107.677 44–0.231 4026.171 96–0.133 82

      表 1  围压为60 MPa时本构模型参数

      Table 1.  Parameters of constitutive model with confining pressure of 60 MPa

    • 通过表1可以看出:无论是加载阶段还是卸载阶段,参数F0均随着循环次数的提高呈非线性增长;但是加载参数m在加载阶段随着循环次数的增加而逐渐减小,在卸载阶段随着循环次数的增加而增大。图3显示了循环加-卸载本构参数的演化规律。

      图  3  循环加-卸载本构参数演化规律

      Figure 3.  Evolution of cyclic loading and unloading constitutive parameters

      在同一组数据中,只有循环次数是不断变化的,使用数据处理软件进行拟合后发现,参数的演化趋势服从

      $y = a\ln \, M + b$

      式中:ab表示拟合参数,M表示循环次数。

      图3(a)图3(c)拟合得到的参数演化规律可以看出,随着循环次数的增多,无论是加载还是卸载,参数F0都随着加载次数的增多呈逐渐上升趋势。根据第1节推导过程中的假设,参数F0与裂纹扩展的强度有关,参数m与岩石裂纹分布有关。随着加载次数增多,岩石逐渐硬化,岩石裂纹扩展时所需的强度逐渐增高,说明参数F0越大,岩石裂纹扩展需要的强度越大,岩石的硬化度越高。岩石并非通常假设的各向同性均质体,其内部存在许多细微裂纹,因此施加外部荷载时,岩石中的部分微裂纹会闭合,当荷载增加到一定程度时,部分闭合的裂纹开始扩展,新裂纹生成、扩展。由于裂纹间并非光滑,存在界面摩擦,随着循环次数增多,部分可张开的裂纹在荷载作用下闭合所需克服的摩擦力减小,所以随着循环次数增多,在荷载作用下旧裂纹扩展、新裂纹生成及扩展越多,并且部分原始裂纹闭合也越多。新裂纹生成及扩展需要的力远远大于原始裂纹闭合需要克服的摩擦力,所以在加载阶段新裂纹生成及扩展数相比原始裂纹闭合的裂纹数要少。在卸载阶段,由于荷载减小,部分闭合的裂纹逐渐张开,而扩展的裂纹和新增裂纹扩展部分并不会因为外部荷载的减小而闭合,所以参数m呈上升趋势;而在加载阶段,由于在荷载作用下闭合的原始裂纹较多,所以整体呈下降趋势。基于上述分析,可以认为,本构参数m与裂纹扩展数呈正相关,裂纹张开数越多,参数m越大。

    • 为了便于应用岩石循环加-卸载本构模型,对围压与岩石参数之间的关系进行了研究,见图4

      图  4  围压与本构模型参数的关系

      Figure 4.  Relationship between confining pressure and constitutive model parameters

      图4中参数${\sigma_{\rm{m}}}$表示围压大小。通过对不同围压下本构模型参数的演化规律进行分析可以发现,参数F0与围压的大小呈线性关系,参数m与围压的大小呈非线性关系。利用围压与参数的关系可以得出第1次加-卸载参数,进而得到对应围压下岩石参数的演化规律,对于深部岩石在长期周期荷载作用下的变形预测具有一定的借鉴意义。

    • 基于参数的演化规律,将式(15)代入式(11),可以得到高围压下岩石的统计损伤本构模型

      ${\sigma _1} = E{\varepsilon _1}{\rm{exp}}\left\{ {{\rm{ - }}{{\left[ {\frac{{{{({\sigma _1}{\rm{ - }}{\sigma _3})}^2}E{\varepsilon _1}}}{{8\left( {C\ln M + D} \right)\left( {{\sigma _1} + {\sigma _3}} \right)\left( {{\sigma _1} - 2\mu {\sigma _3}} \right)}}} \right]}^{{{A{\rm ln} \;M}} + {{B}}}}} \right\} + 2\mu {\sigma _3}$

      式中:ABCD均为本构模型参数,加载阶段的参数与卸载阶段的参数不相同,加载阶段ABCD的具体数值如图3(a)图3(b)所示,卸载阶段的参数数值如图3(c)图3(d)所示。如果将拟合得到的参数代入本构模型中对比理论值与原始数据只能说明拟合的效果,而不能表明本构模型的正确性。为此基于上述偶数次加载参数的拟合公式,求出奇数次循环的参数值,将奇数次循环的参数代入本构模型中得出本构的理论值。通过对比加-卸载曲线的原始数据与本构对应的理论值来验证本构模型的正确性,如图5所示。

      图  5  奇数次循环的加-卸载曲线理论值与原始数据对比

      Figure 5.  Comparison of theoretical values of loading and unloading curves with original data for odd cycles

      对比理论值与原始数据发现,基于Griffith准则的统计损伤本构模型能够较好地描述和表征高围压下岩石循环加-卸载应力-应变曲线,说明了本构模型的合理性及正确性。对于围压为20、40 MPa的循环加-卸载,采用同样的方法进行研究,除了F0m的大小不同,不同循环下的参数演化规律也服从式(15),故在此不再一一列举。以Griffith准则建立的本构模型理论值与高围压下岩石变形的试验值较吻合,从而为描述岩石在高围压下的应力-应变曲线提供了新的思路。

    • (1)根据Griffith准则,岩石内部存在许多细微裂纹,岩石破坏是由于裂纹的扩展。已有的研究也表明,岩石内部裂纹数服从Weibull分布。因此,基于Griffith准则建立的统计损伤本构模型与岩石的实际情况更相符。该模型参数少,参数演化规律明确,且便于拟合。该研究方法对高围压下受周期荷载作用的岩石变形预测具有一定借鉴意义。

      (2)通过本构模型的等式变换,确定了不同围压下岩石的损伤阈值,损伤阈值处的轴向应变值较小,并且不同围压下产生损伤时的应变值随着围压增大呈现非线性增长。

      (3)对本构模型参数的演化规律进行分析发现:参数F0与围压呈线性关系,并且F0与裂纹扩展强度呈正相关性;参数m与裂纹张开数目相关,裂纹张开数目越多,m越大;参数F0随着围压的增大呈线性增长,参数m随围压的增加呈非线性增长。

      (4)采用偶数次循环的加载和卸载曲线对本构参数进行分析,得到本构参数在不同循环次数和不同围压下的演化规律,由此计算出奇数次循环的加载和卸载曲线本构理论值,并与试验曲线对比,验证了本构模型的合理性和正确性。

参考文献 (18)

目录

    /

    返回文章
    返回