DOI: 10.11817/j.issn.1672-7207.2019.06.020

基于Normal分布的岩石统计热损伤本构模型研究

朱振南¹，蒋国盛¹，田红¹，吴文兵^1,2，梁荣柱^1,2，窦斌¹

(1. 中国地质大学(武汉) 工程学院，湖北 武汉，430074；

2. 浙江大学 滨海和城市岩土工程研究中心，浙江 杭州，310058)

摘 要：

tre应变等价性理论的岩石损伤模型，假定受热损伤的岩石微元强度服从Normal分布，考虑温度对岩石力学参数的影响，引入热损伤变量，在微元破坏符合Mohr-Coulomb准则条件下，建立高温作用后岩石统计热损伤本构模型。依据岩石屈服的概念，利用极值法确定模型参数，并与不同温度下花岗岩的压缩试验结果进行对比分析。研究结果表明：本模型所得曲线与三轴压缩试验曲线较吻合，并能很好地反映高温作用后岩石软化特征；本模型不包含非常规岩石力学参数，便于工程应用；研究成果能够为高温岩体工程的开发和防护提供参考。

关键词：

高温；岩石；热损伤；Normal分布；本构模型；

中图分类号：TU45             文献标志码：A         文章编号：1672-7207(2019)06-1411-08

Study on statistical thermal damage constitutive model of rock based on normal distribution

ZHU Zhennan¹, JIANG Guosheng¹, TIAN Hong¹, WU Wenbing^{1, 2}, LIANG Rongzhu^{1, 2}, DOU Bin¹

(1. Faculty of Engineering, China University of Geosciences (Wuhan), Wuhan 430074, China;

2. Research Center of Coastal and Urban Geotechnical Engineering, Zhejiang University, Hangzhou 310058, China)

Abstract: Based on rock damage model founded by Lemaitre’s strain equivalent theory and hypothesis that rock particles intensity obeys normal random distribution after the thermal treatment, a statistical thermal damage constitutive model for rocks was established by adopting Mohr-Coulomb criterion with the thermal damage variable being introduced. The parameters of the model parameters were determined by the extreme value method according to the conception of yielding. The rationality of the model was verified by the comparison of the calculated results with those of the uniaxial experiments and conventional triaxial experiments of granite after different high temperature treatments. The results show that the constitutive model can reflect the characteristics of rock strain softening rock at different temperatures. The model does not contain any unconventional rock mechanics parameters, and can be used easily in engineering. The researches can provide reference for the development and protection of high-temperature rock projects.

Key words: high temperature; rock; thermal damage; Normal distribution; constitutive model

随着高放核废料地下深埋处置、干热岩地热资源开发和矿产资源深部开采等高温岩体工程的进一步发展，人们对高温作用下或高温作用后岩石物理力学特征的相关研究越来越多。与常温状态下相比，高温下或高温后岩石的力学性质通常会发生劣化，进而影响地下高温岩体工程的稳定与安全。研究者从不同角度和层次，通过宏观及微观试验研究了温度对岩石物理力学性质的影响。通过宏观力学试验研究发现花岗岩^[1-3]、大理岩^[4-6]、砂岩^[7-9]、黏土岩^[10]等的力学性能如弹性模量、抗压强度随温度升高而呈现出弱化的趋势。PENG等^[11-13]利用扫描电子显微镜(SEM)观测得出，不同岩石内部的微裂纹会随温度的升高而不断生成、扩展、增多，并且微裂纹类型由晶间裂纹逐渐转变为穿晶裂纹。这些试验研究揭示了温度对岩石具有热损伤作用，因此，在高温岩体工程设计与计算中，为了更精确地描述岩石的力学行为，采用考虑温度效应的岩石本构模型十分必要。目前，考虑温度作用的本构模型主要分3大类：1) 基于损伤理论和热弹性理论建立的热弹塑性损伤模型^[14-16]；2) 基于蠕变试验和热黏塑性理论建立的流变模型^[17-18]；3) 基于统计分布建立的岩石热损伤本构模型^[19-20]。统计损伤本构模型可以直接而又准确地描述岩石损伤演化过程的缺陷，从而更好地刻画岩石损伤的力学机制，因而该类模型较前2类模型具有明显的优势。在岩石统计损伤本构模型^[21-26]的基础上，徐银花等^[19]通过考虑温度对统计参数的影响建立了高温下岩石损伤本构模型，但仅考虑了单向压力，对三向应力状态的适用性没有进一步验证；李天斌等^[20]通过考虑温度对统计参数的影响建立了考虑温度启裂应力的热-力-损伤本构模型，但未验证其高温的适用性。为此，本文作者在传统岩石统计本构模型基础上，假设受热损伤的岩石微元强度服从Normal分布，引入热损伤变量，建立满足Mohr- Coulomb准则的岩石统计热损伤本构模型。通过与高温作用后花岗岩单轴和三轴压缩试验结果进行对比，验证所提模型的合理性和适用性。

1  模型建立

1.1  本构关系和损伤变量

岩石内部含有微裂隙和孔洞等天然缺陷，高温作用后会引起这些缺陷不断破坏、连接和贯通。依据损伤理论，高温对岩石的损伤为统计损伤，受热损伤的岩石微元体强度服从统计分布。本文在常温损伤变量D定义^[22]的基础上，提出热损伤变量D_T，为岩石经历高温作用后在某一应力状态下已破坏微元体数目N_Tf与总微元体数目N_T的比值，即

                (1)

假设经高温作用后微元体强度概率密度函数为P(F_T)，当应力水平F_T超过一定强度时，单元相继破裂。在区间[F_T，F_T+dF_T]内，若新破坏的微元体数目为N(F_T)dF_T，则D_T为

  (2)

高温作用后，岩石在荷载作用下，应变等效性假设示意图如图1所示(其中，A为岩石初始横截面总面积，A^*为岩石受损后的有效承载面积，A-A^*为在加载过程中出现损伤后的截面面积)。由式(1)可得

           (3)

图1  应变等效性假设示意图

Fig. 1  Illustration of strain equivalence hypothesis

根据Lemaitre应变等价理论和有效应力的概念^[27]，此假说认为名义应力下含损伤岩石产生的应变等价于有效应力下含损伤岩石产生的有效应变(其中名义应力即为试验测的应力)，则考虑温度作用后含损伤的名义应力(i=1，2和3)和有效应力(i=1，2和3)的关系为

i=1, 2, 3           (4)

经历高温作用后岩石应力-应变曲线依然存在明显的弹性阶段^[10]，且力学参数(弹性模量E、泊松比v)受温度影响^[1-10]。设依赖于损伤温度的弹性模量和泊松比分别为E_T和v_T，并假定高温作用后未受损伤的岩石部分处于完全弹性状态，则根据广义胡克定律可知

          (5)

式中：(i，j，k)为(1，2，3)，(2，3，1)或(3，1，2)；和分别为不同温度作用后岩石有效弹性模量和有效泊松比；为与有效应力所对应的有效应变。

由式(4)和(5)可得

     (6)

根据Lemaitre应变等价理论可得

i=1, 2, 3              (7)

其中：ε_i为与有效应力σ_i所对应的名义应变。岩石弹性模量定义为岩石处在弹性阶段轴向应力和轴向应变之比，因此，可认为

                 (8)

在不同温度下，岩石泊松比可以表示为v_T=,同时，岩石有效岩石泊松比表示为，因而，

                 (9)

将式(4)，(7)~(9)代入式(6)可得出

    (10)

对式(10)进行简化可得到岩石材料经历热损伤作用的各向同性弹性损伤本构方程：

      (11)

式中：(i，j，k)为(1，2，3)，(2，3，1)或(3，1，2)；E_T和分别为岩石经历不同温度作用后的弹性模量和泊松比。对于含初始损伤的岩石，D_T的取值范围为0≤D_T≤1，其中，D_T=0表示岩石材料未经高温作用处于初始无损状态，D_T=1表示岩石材料经高温作用后处于完全损伤状态。

应用此热损伤本构方程时，需要着重考虑以下2个问题：1) 高温作用后岩石微元体应力水平F_T的确定；2) 高温作用后微元强度概率密度函数P(F_T)的选择。

1.2  应力水平

经高温作用后岩石微观破坏准则可用有效应力和岩石材料参数表示：

             (12)

式中：K_T为与经高温作用后岩石材料参数有关的常数；为岩石内部微元体的有效应力。当≥K_T时，岩石微元单元产生破坏，因此，可以作为微元体应力，用来反映岩石微观破坏程度，即

              (13)

曹文贵等^[22]提出基于岩石破坏准则的新型岩石微元强度度量方法，该方法能直接反映应力状态对岩石微元强度的影响。试验结果表明，岩石的内摩擦角随温度的升高而增大，黏聚力随温度的升高而降低^[2]。Mohr-Coulomb准则具有参数形式简单且适用于岩石介质等特点，利用有效应力和岩石力学参数表示的考虑温度效应的Mohr-Coulomb准则^[28]可表示为

   (14)

式中：和c_T分别为依赖温度的岩石内摩擦角和黏聚力；(i=1，2和3)为有效应力下3个主应力分量。因此，经高温作用后岩石微元应力水平F_T可表示为

   (15)

式(15)包含有效应力，因此，需要得到用名义应力表示的应力。根据式(4)可得：

              (16)

              (17)

由式(11)可得

        (18)

将式(18)分别代入式(16)和式(17)可得：

          (19)

          (20)

将式(19)和式(20)代入式(15)则可得到基于Mohr- Coulomb准则的高温作用后微元体应力：

    (21)

1.3  统计热损伤本构模型

岩石内包含的裂隙和许多其他缺陷分布方向是随机的，因此，可以在计算方程中引入反映强度统计特点的函数^[29]。KRAJCINOVIC等^[30]将连续损伤理论和统计强度有机地结合起来，提出1个统计损伤本构模型。唐春安^[21]提出常温状态下岩石微元体强度符合Normal分布，本文假设高温作用后岩石微元体强度服从Normal分布，类比常温状态下概率密度函数^[24]，在高温下，概率密度函数表示为

   (22)

式中：F_T为微元强度随机分布的应力；S₀及F_n0为不同温度作用后Normal分布函数参数，它们反映岩石材料经历不同温度作用后的力学性质。

将式(22)代入式(2)可得

   (23)

将式(23)代入式(11)，根据广义胡克定律，基于Normal分布下高温作用后岩石统计热损伤的本构方程则可表示为

                          (24)

2  参数确定

本文所提出的岩石统计热损伤模型涉及到的力学参数包括E_T，v_T，c_T和，均可通过高温作用后岩石常规三轴试验测得。在不同温度下，Normal分布参数S₀和F_n0可以依据LI等^[25]提出的极值法确定，具体求解过程如下。

岩石于低围压条件下表现出屈服特征，应力达到峰值后随应变的增加而降低。因此，对应峰值点的导数等于0，即

             (25)

式中：和分别表示峰值应力及相对应的峰值应变。将峰值点坐标(，)代入式(21)和式(24)，可得

 (26)

      (27)

式中：F_sc为且时的岩石微元体强度。由式(19)，(20)和(24)可得出

         (28)

结合式(23)和(26)，D_T可以表示为

  (29)

式中：代表Normal分布标准函数。将峰值点坐标参数和代入式(29)，可得

    (30)

式中：D_sc为且时的热损伤变量；F_sc如式(26)所示。

结合式(26)，(27)和(30)，式(25)可以表示成

    (31)

假设式(30)中Z=(F_sc-F_n0)/S₀，Z可查看标准正态分布表获得，因此，结合式(30)和式(31)，不同高温作用后S₀可表示为

   (32)

同时，不同温度作用后F_n0可表示为

             (33)

3  模型验证与讨论

本文引用秦严^[31]通过高温后花岗岩单轴和三轴压缩试验的相关研究成果，以验证本文所提出的岩石统计热损伤本构模型的合理性及适用性。

在此项研究中，岩样常温下纵波波速为4 800 m/s，完整性和均匀性较好，平均密度为2.72 g/cm³，平均单轴抗压强度为83.33 MPa，泊松比为0.141，内摩擦角为48.38°，黏聚力为19.81 MPa。X线衍射分析其矿物成分的质量分数为：石英42%，斜长石39%，微斜长石14%，云母5%。在试验过程中，高温处置的升温速率为3 ℃/min，目标温度为100，200，300，400和500 ℃，岩样温度达到目标温度后恒温2 h，再在高温炉内自然冷却，三轴试验的围压为5，15和25 MPa。各个温度和围压作用下试验样本数均为3个。

根据秦严^[31]的试验结果，本文所建模型参数取值如表1所示。表1中，各温度和围压作用下Normal分布参数S₀和F_n0通过式(32)和式(33)得出。由式(24)可得不同温度和不同围压作用下岩石应力-应变关系的理论曲线，见图2。从图2可以看出：在同一围压下，花岗岩抗压强度和弹性模量随温度升高逐渐降低，说明高温作用后花岗岩力学性质发生劣化，高温对岩石具有损伤作用，岩石的峰值应变随温度增大而增大；应力-应变曲线随温度增大而向右移动，表明随温度的增大岩石塑性变形增大；不同温度作用后，花岗岩抗压强度和弹性模量随围压升高而增大，围压具有提高岩石强度的作用。通过对比试验数据与理论曲线可知：试验结果与理论结果拟合程度较好，理论结果可以反映花岗岩强度和变形随温度和围压变化的趋势，表明本文提出的考虑温度效应的统计热损伤本构模型能够较充分地反映岩石峰后软化特征和岩石强度依赖于温压状态的特征，同时能更好地表征不同温度和不同围压条件下岩石应力-应变曲线全过程，且围压越高，吻合程度越高。

在求解方法上，本文依据岩石屈服概念，利用极值法确定模型参数。此确定方法适用于不同温度和围压作用下加载试验数据的求解，不包含非常规岩石力学参数，工程应用更加方便。在岩石统计热损伤模型建立过程中，模型所用力学参数和统计分布参数与岩性无关，故本模型同样适用于其他不同岩性的岩石。

值得注意的是：在进行单轴试验时，岩石应力-应变曲线通常存在原生裂隙和孔隙的压密阶段，表现出明显的非线性特征；本文所提出的岩石统计热损伤本构模型是基于广义胡克定律建立的，认为岩石受压后即进入弹性阶段，忽略了岩石原生裂隙、孔隙及受热损伤所产生微裂隙闭合的影响，因此，在单轴压缩条件下，理论曲线和试验数据之间存在一定偏差，偏差在合理的范围内，这与徐银花等^[19]由所建立的单轴压缩条下模型所得结果及试验结果一致。进行三轴压缩试验时，在围压作用下，岩石中的微裂隙和孔隙发生闭合，应力-应变曲线压密阶段不明显，且随着围压升高，应力-应变曲线起始阶段越呈现线性特征，所以，理论曲线与三轴压缩试验曲线较吻合。

表1  不同温度下损伤岩石的模型计算参数

Table 1  Model parameters of thermally treated rocks in different temperatures

图2  不同温度下试验和理论曲线的对比

Fig. 2  Comparison between calculated and experimental curves after thermal treatment at different temperatures

4  结论

1) 本文所建立的本构模型能够充分反映了岩石的峰后软化和残余强度特征，很好地描述了不同温度和压力下全应力-应变曲线，围压越高，吻合程度越高。

2) 本文模型所需参数均可根据高温岩石单轴和三轴压缩试验结果确定，不包含非常规岩石力学参数，工程应用方便，并且与岩性无关，具有适用性广的特点。

3) 模型忽略了岩石原生裂隙、孔隙及受热损伤所产生微裂隙的影响。在单轴压缩条件下，理论结果和单轴压缩试验结果存在一定偏差；在三轴压缩条件 下，随着围压升高，岩石原生裂隙和孔隙及热损伤所产生的微裂隙逐渐被压实，理论结果与常规试验结果较吻合。

参考文献：

[1] 郤保平, 赵阳升. 600 ℃内高温状态花岗岩遇水冷却后力学特性试验研究[J]. 岩石力学与工程学报, 2010, 29(5): 892-898.

XI Baoping, ZHAO Yangsheng. Experimental research on mechanical properties of water-cooled granite under high temperatures within 600 ℃[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(5): 892-898.

[2] 徐小丽, 高峰, 张志镇. 高温作用后花岗岩三轴压缩试验研究[J]. 岩土力学, 2014, 35(11): 3177-3183.

XU Xiaoli, GAO Feng, ZHANG Zhizhen. Research on triaxial compression test of granite after high temperatures[J]. Rock and Soil Mechanics, 2014, 35(11): 3177-3183.

[3] 李二兵, 王永超, 陈亮, 等. 北山花岗岩热损伤力学特性试验研究[J]. 中国矿业大学学报, 2018, 47(4): 735-741.

LI Erbing, WANG Yongchao, CHEN Liang, et al. Experimental study of mechanical properties of Beishan granite’s thermal damage[J]. Journal of China University of Mining & Technology, 2018, 47(4): 735-741.

[4] 郭清露, 荣冠, 姚孟迪, 等. 大理岩热损伤声发射力学特性试验研究[J]. 岩石力学与工程学报, 2015, 34(12): 2388-2400.

GUO Qinglu, RONG Guan, YAO Mengdi, et al. Experimental study on acoustic emission behaviors and mechanical properties of thermal damaged marbles[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(12): 2388-2400.

[5] 陈国飞, 杨圣奇. 高温后大理岩力学性质及其破坏规律研究[J]. 工程力学, 2014, 31(8): 189-196.

CHEN Guofei, YANG Shengqi. Study on failure mechanical behavior of marble after high temperature[J]. Engineering Mechanics, 2014, 31(8): 189-196.

[6] SU Naijing, JING Hongwen, YIN Qian, et al. Strength and deformation behaviors of veined marble specimens after vacuum heat treatment under conventional triaxial compression[J]. Acta Mechanica Sinica, 2017, 33(5): 886-898.

[7] 吴刚, 王德咏, 翟松韬. 单轴压缩下高温后砂岩的声发射特征[J]. 岩土力学, 2012, 33(11): 3237-3242.

WU Gang, WANG Deyong, ZHAI Songtao. Acoustic emission characteristics of sandstone after high temperature under uniaxial compression[J]. Rock and Soil Mechanics, 2012, 33(11): 3237-3242.

[8] 苏承东, 韦四江, 秦本东, 等. 高温对细砂岩力学性质影响机制的试验研究[J]. 岩土力学, 2017, 38(3): 623-630.

SU Chengdong, WEI Sijiang, QIN Bendong, et al. Experimental study of influence mechanism of high temperature on mechanical properties of sandstone[J]. Rock and Soil Mechanics, 2017, 38(3): 623-630.

[9] 尹土兵, 李夕兵, 殷志强, 等. 高温后砂岩静、动态力学特性研究与比较[J]. 岩石力学与工程学报, 2012, 31(2): 273-279.

YIN Tubing, LI Xibing, YIN Zhiqiang, et al. Study and comparison of mechanical properties of sandstone under static and dynamic loadings after high temperature[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(2): 273-279.

[10] TIAN Hong, ZIEGLER M, KEMPKA T. Physical and mechanical behavior of claystone exposed to temperatures up to 1 000 °C[J]. International Journal of Rock Mechanics & Mining Sciences, 2014(70): 144-153.

[11] PENG Jun, RONG Guan, CAI Ming, et al. Physical and mechanical behaviors of a thermal-damaged coarse marble under uniaxial compression[J]. Engineering Geology, 2016, 200(12): 88-93.

[12] YANG Shengqi, XU Peng, LI Yanbao, et al. Experimental investigation on triaxial mechanical and permeability behavior of sandstone after exposure to different high temperature treatments [J]. Geothermics, 2017(69): 93-109.

[13] 吴刚, 翟松韬, 李玉寿, 等. 高温下大理岩受压破坏的细观结构分析[J]. 岩石力学与工程学报, 2012, 31(增2): 3579-3585.

WU Gang, ZHAI Songtao, LI Yushou, et al. Mesostructure investigation on marbles under uniaxial compression at high temperatures[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(Suppl 2): 3579-3585.

[14] HUECKEL T, PEANO A, PELLEGRINI R. A thermo-plastic constitutive law for brittle-plastic behavior of rocks at high temperatures[J]. Pure & Applied Geophysics, 1994, 143(1/2/3): 483-510.

[15] 朱元广, 刘泉声, 张程远, 等. 岩石非线性黏弹性蠕变特性的时温等效效应[J]. 岩土力学, 2012, 33(8): 2303-2309.

ZHU Yuanguang, LIU Quansheng, ZHANG Chengyuan, et al. Nonlinear viscoelastic creep property of rock with time- temperature equivalence effect[J]. Rock and Soil Mechanics, 2012, 33(8): 2303-2309.

[16] 张岩, 李宁, 于海鸣, 等. 温度应力对裂隙岩体强度的影响研究[J]. 岩石力学与工程学报, 2013, 32(S1): 2660-2668.

ZHANG Yan, LI Ning, YU Haiming, et al. Research on influence of thermal stress on fractured rock mass strength[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(S1): 2660-2668.

[17] 郤保平, 赵阳升, 万志军, 等. 热力耦合作用下花岗岩流变模型的本构关系研究[J]. 岩石力学与工程学报, 2009, 28(5): 956-967.

XI Baoping, ZHAO Yangsheng, WAN Zhijun, et al. Study of constitutive equation of granite rheological model with thermo-mechanical coupling effects[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(5): 956-967.

[18] 赵延林, 唐劲舟, 付成成, 等. 岩石黏弹塑性应变分离的流变试验与蠕变损伤模型[J]. 岩石力学与工程学报, 2016, 35(7): 1297-1308.

ZHAO Yanlin, TANG Jinzhou, FU Chengcheng, et al. Rheological test of separation between viscoelastic-plastic strains and creep damage model[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(7): 1297-1308.

[19] 徐银花, 徐小丽. 高温下岩石损伤本构模型研究[J]. 广西大学学报(自然科学版), 2017, 42(1): 226-235.

XU Yinhua, XU Xiaoli. Study on damage constitutive model of rock under high temperature[J]. Journal of Guangxi University(Natural Science Edition), 2017, 42(1): 226-235.

[20] 李天斌, 高美奔, 陈国庆, 等. 应脆性岩石热-力-损伤本构模型及其初步运用[J]. 岩土工程学报, 2017, 39(8): 1477-1484.

LI Tianbin, GAO Meiben, CHEN Guoqing, et al. A thermal-mechanical-damage constitutive model for hard brittle rocks and its preliminary application[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(8): 1477-1484.

[21] 唐春安. 岩石破裂过程中的灾变[M]. 北京: 煤炭工业出版社, 1993: 13-31.

TANG Chunan. Catastrophe in rock unstable failure[M]. Beijing: China Coal Industry Publishing House, 1993: 13-31.

[22] 曹文贵, 方祖烈, 唐学军. 岩石损伤软化统计本构模型之研究[J]. 岩石力学与工程学报, 1998, 17(6): 628-633.

CAO Wengui, FANG Zulie, TANG Xuejun. A study of statistical constitutive model for soft and damage rocks[J]. Chinese Journal of Rock Mechanics and Engineering, 1998, 17(6): 628-633.

[23] 曹文贵, 张超, 贺敏, 等. 考虑空隙压密阶段特征的岩石应变软化统计损伤模拟方法[J]. 岩土工程学报, 2016, 38(10): 1754–1761.

CAO Wengui, ZHANG Chao, HE Min, et al. Statistical damage simulation method of strain softening deformation process for rocks considering characteristics of void compaction stage[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(10): 1754–1761.

[24] CAO Wengui, LI Xiang, ZHAO Heng. Damage constitutive model for strain-softening rock based on normal distribution and its parameter determination[J]. Journal of Central South University of Technology, 2007, 14(5): 719-724.

[25] LI Xiang, CAO Wengui, SU Yonghua. A statistical damage constitutive model for softening behavior of rocks[J]. Engineering Geology, 2012(143/144): 1-17.

[26] ZHANG Hongwei, WAN Zhijun, WANG Chaoyi, et al. A Statistical damage constitutive model for geomaterials under plane-strain biaxial stress state[J]. Advances in Mathematical Physics, 2017(8): 1-10.

[27] LEMAITRE J. How to use damage mechanics[J]. Nuclear Engineering and Design, 1984, 80(2): 233-245.

[28] TIAN Hong. Development of a thermo-mechanical model for rocks exposed to high temperatures during underground coal gasification[D]. Aachen, Germany: RWTH Aachen University. Department of Geotechnical Engineering, 2013: 61-81.

[29] 皮萨林科F C, 列别捷夫 A A. 复杂应力状态下的材料变形与强度[M]. 江明行, 译. 北京: 科学出版社, 1983: 129-140.

ПИСАРЕНКО F С, ЛЕБЕДЕВ A A. Deformation and strength of materials under complex stress state[M]. JIANG Mingxing, trans. Beijing: Science Press, 1983: 129-140.

[30] KRAJCINOVIC D, SILVA M A G. Statistical aspects of the continuous damage theory[J]. International Journal of Solids & Structures, 1982, 18(7): 551-562.

[31] 秦严. 高温后花岗岩物理力学性质试验研究[D]. 北京:中国地质大学(北京)工程技术学院, 2017: 33-52.

QIN Yan. High-temperature influence on physical and mechanical properties of diorite[D]. Beijing: China University of Geosciences(Beijing). School of Engineering and Technology, 2017: 33-52.

(编辑  陈灿华)

收稿日期：2018-08-08；修回日期：2018-10-12

基金项目(Foundation item)：国家自然科学基金资助项目(41602374, 41674180)；工程地质与岩土防护学术创新基地岩土钻掘与防护教育部工程研究中心开放研究基金资助项目(201703)；中央高校基本科研业务费专项资金资助项目(CUG2170207)；中国地质大学(武汉)教学实验室开放基金资助项目(SKJ2018132) (Projects(41602374, 41674180) supported by the National Natural Science Foundation of China; Project(201703) supported by the Open Fund of Engineering Research Center of Rock-Soil Drilling & Excavation and Protection, Department of Education; Project (CUG2170207) supported by the Fundamental Research Funds for the Central Universities, China University of Geosciences(Wuhan); Project (SKJ2018132) supported by the Teaching Laboratory Open Foundation, China University of Geosciences)

通信作者：田红，博士，副教授，硕士生导师，从事高温岩石力学及数值模拟研究；E-mail:htian@cug.edu.cn

摘要：基于Lemaitre应变等价性理论的岩石损伤模型，假定受热损伤的岩石微元强度服从Normal分布，考虑温度对岩石力学参数的影响，引入热损伤变量，在微元破坏符合Mohr-Coulomb准则条件下，建立高温作用后岩石统计热损伤本构模型。依据岩石屈服的概念，利用极值法确定模型参数，并与不同温度下花岗岩的压缩试验结果进行对比分析。研究结果表明：本模型所得曲线与三轴压缩试验曲线较吻合，并能很好地反映高温作用后岩石软化特征；本模型不包含非常规岩石力学参数，便于工程应用；研究成果能够为高温岩体工程的开发和防护提供参考。