DOI： 10.11817/j.issn.1672-7207.2021.05.026

考虑黏弹塑性应变分离的岩石复合蠕变模型研究

张亮亮¹，王晓健^{1, 2}

(1. 安徽理工大学 土木建筑学院，安徽 淮南，232001；

2. 安徽理工大学 省部共建深部煤矿采动响应与灾害防控国家重点实验室，安徽 淮南，232001)

摘要：依据不同偏应力下石灰岩循环加卸载蠕变试验数据提出和建立瞬时弹性体(IEB)、瞬时塑性体(IPB)、非线性黏弹性体(NVEB)和黏塑性体(VPB)分别用于描述瞬时弹性应变、瞬时塑性应变、非线性黏弹性应变和黏塑性应变，并将这四者串联，建立能够同时反映岩石加、卸载过程的复合蠕变模型。根据非线性流变理论推导岩石在一维、三维应力状态下的复合蠕变模型加卸载蠕变本构方程，最后采用试验数据对模型合理性进行验证。研究结果表明：当复合蠕变模型综合考虑瞬时弹性应变、瞬时塑性应变、非线性黏弹性应变和黏塑性应变时，模型理论曲线与试验曲线高度吻合；所建模型不仅能够描述石灰岩加载过程中非线性蠕变特征，还能描述卸载过程中应变随时间的变化规律，不同偏应力下试验曲线与理论曲线相关系数均大于0.968，对比结果验证复合蠕变模型的合理性和可行性。

关键词：岩石力学；蠕变；分级加卸载；黏弹塑性；非线性

中图分类号：TU45            文献标志码：A               开放科学(资源服务)标识码(OSID)

文章编号：1672-7207（2021）05-1655-11

Study on composite creep model of rock considering viscoelastic plastic strain separation

ZHANG　Liangliang¹, WANG　Xiaojian^{1, 2}

(1. School of Civil Engineering and Architecture, Anhui University of Science and Technology, Huainan 232001, China;

2. State Key Laboratory of Mining Response and Disaster Prevention and Control in Deep Coal Mine, Anhui University of Science and Technology, Huainan 232001, China)

Abstract: Based on the cyclic loading and unloading creep test data at different deviator stress levels of limestone, the instantaneous elastic body(IEB), the instantaneous plastic body(IPB), the nonlinear viscoelastic body(NVEB) and the viscoplastic body(VPB) were proposed and established to describe the instantaneous elastic strain, the instantaneous plastic strain, the nonlinear viscoelastic strain and the viscoplastic strain respectively. And they were connected in series to establish a composite creep model which could simultaneously reflect the loading and unloading creep process of rock. According to the nonlinear rheological theory, the loading and unloading constitutive equations of the composite creep model of rock in one-dimensional and three-dimensional stress state were deduced. Finally, the rationality of the model was verified by the test data. The results show that the theoretical curve of the model is highly consistent with the test curve, when the composite creep model considers the instantaneous elastic strain, the instantaneous plastic strain, the nonlinear viscoelastic strain and the viscoplastic strain. The model can not only describe the nonlinear creep characteristics in the loading process of limestone, but also describe the change rule of strain with time in the unloading process, and the correlation coefficients between the test curve and the theoretical curve are greater than 0.968 under different deviator stresses. The comparison results verify the rationality and feasibility of the composite creep model.

Key words: rock mechanics; creep; step loading and unloading; viscoelastic plasticity; nonlinearity

岩石蠕变力学特性是岩石固有属性之一，尤其是深部赋存条件复杂的软岩，受高地应力、高温和高水压影响，其节理裂隙往往发育较好，宏观上表现为显著的蠕变变形，该变形特性对服务期限较长结构的安全性和长期稳定性至关重要^[1-2]，因而，在进行深部地下结构(地下硐室群、核废料储存室及水电站地基等)设计时，不仅要考虑荷载作用下的瞬时变形，还要着重考虑结构随时间的蠕变变形。目前，关于岩石蠕变特性研究主要包括室内单、三轴分级加载压缩蠕变试验研究以及根据试验结果建立相应的元件(弹性元件、黏性元件和塑性元件)组合模型理论研究。如韩阳等^[3]提出了一种非线性黏性元件，将其替换Burgers模型中的2个线性黏性元件，建立了一种可描述岩石三阶段蠕变特性的非定常参数Burgers模型；陈卫忠等^[4]基于泥岩现场大型真三轴蠕变试验结果，提出了一种能够真实反映深部软岩不同应力下的流变特征的非线性经验幂函数型蠕变模型；HU等^[5]基于分数阶理论和损伤理论，提出了一种非线性损伤蠕变模型，该模型可描述硬岩加速蠕变阶段；刘开云等^[6]在Maxwell模型基础上串联1个非线性黏塑性体，建立了可以描述岩石非稳态蠕变的参数非定常非线性黏弹塑性蠕变本构模型；刘东燕等^[7]基于不同围压作用下重庆地区砂岩蠕变结果，将黏滞系数随时间按幂函数衰减的经验模型与Kelvin模型串联，建立了可描述不同围压作用下加速蠕变特征的改进黏塑性模型；ROBERTS等^[8]针对盐岩进行了三轴压缩、拉伸和循环拉压蠕变试验，研究了不同加载方式及应力路径对盐岩蠕变及扩容特性的影响；SINGH等^[9]基于Maxwell模型采用声发射技术获得了盐岩单轴压缩过程中的弹性模量和黏性系数，并根据该参数预测盐岩的单轴蠕变行为；BOUKHAROV等^[10]分析了脆性岩石流变曲线中的3个典型过程，并将应力触发的虎克体和应变触发的非线性黏壶元件串联用于描述非稳态蠕变过程。通过对上述研究成果及目前岩石蠕变力学发展过程的分析总结，本文作者认为仍存在以下2个方面不足：

1) 缺乏岩石循环加卸载蠕变试验研究。在地下硐室开挖过程中，围岩原始应力平衡状态被打破，围岩由于卸载应力得到释放产生较大变形，随硐室进一步开挖，扰动围岩继续受开挖过程中爆破及机械振动等荷载，此过程相当于加载过程；再如大型水电站基础，在蓄水过程中，水压通过大坝将部分压力传递给坝基，而在放水过程中这部分压力又会随水位的降低逐渐释放。实际岩体工程多数受循环加卸载作用，而目前针对不同类型岩石的单、三轴分级加载试验研究成果较多，但关于循环加卸载作用下岩石变形响应的试验报道却相对较少。李娜等^[11-13]进行了该方面的研究，但由于岩石种类繁多、地质条件复杂，所取得的成果应用受限，因此，在不同应力、不同围压及其他复杂应力路径下岩石循环加卸载蠕变试验尚需进一步研究。

2) 岩石蠕变模型建立尚需完善。目前岩石蠕变模型理论最成熟、应用最广泛的是元件组合模型，该模型将弹性元件、黏性元件和塑性元件通过不同的串并联方式进行组合来反映岩石在不同应力下的蠕变特性。但研究发现经典元件组合模型(Bingham，Burgers和Nishihara模型等)只能描述岩石在低应力状态下发生的稳态蠕变，而不能描述岩石在高应力状态下发生的非稳态蠕变，尤其是非线性特征尤为明显的加速蠕变阶段。因此，目前大多数研究重点是基于非线性损伤理论，采用非线性黏性元件代替线性元件来对传统模型进行改进和完善，且多数研究都是建立加载条件下岩石的蠕变模型，针对岩石卸载条件下的蠕变特征研究及相应条件下的蠕变本构方程研究较少。同时，不少研究都采用弹性体来描述岩石在加载过程中产生的瞬时应变，认为该应变都为瞬时弹性应变，忽略了瞬时塑性应变，赵延林等^[14]的研究表明瞬时塑性应变约占瞬时应变的一半且占比随应力的提高而提高，因此，仅考虑瞬时弹性应变忽视瞬时塑性应变所建立的蠕变模型不够严谨。

本文作者以文献[15]中石灰岩循环加卸载蠕变试验数据为基础，采用瞬时弹性体(IEB)、瞬时塑性体(IPB)、非线性黏弹性体(NVEB)和黏塑性体(VPB)分别描述瞬时弹性应变、瞬时塑性应变、非线性黏弹性应变和黏塑性应变，并将四者串联建立复合蠕变模型。最后，根据非线性岩石流变理论推导岩石单、三轴蠕变及卸载方程，结合试验曲线与模型理论曲线对比结果验证模型的合理性。

1  石灰岩循环加卸载蠕变试验数据分析

赵延林等^[14]采用MTS815岩石试验系统对石灰岩进行了不同围压下的循环加卸载蠕变试验，石灰岩基本力学参数及蠕变试验中偏应力分别见表1和表2。

表1　石灰岩基本力学参数

Table 1　Basic mechanical parameters of limestone

表2　石灰岩蠕变试验各级偏应力

Table 2　Different deviator stresses in creep test of limestone

围压分别为4 MPa和8 MPa的L2和L7试样的石灰试样的循环加卸载蠕变曲线如图1所示，图1中，σ₃为第3主应力。由图1可见：在前2级应力下，蠕变曲线由瞬时应变通过衰减蠕变进入等速蠕变后，蠕变速率基本为0，卸载后瞬时恢复应变比加载瞬时应变小，说明加载过程中不仅有瞬时弹性应变，还有瞬时塑性应变；后3级应力下等速蠕变曲线斜率明显大于0，且卸载时的永久残余应变大于加载瞬时塑性应变，这表明在蠕变过程中，不仅有黏弹性应变还存在黏塑性应变。

图1　不同围压下石灰岩试样的循环加卸载蠕变曲线

Fig. 1　Cyclic loading unloading creep curves of limestone samples at different confining pressures

典型的岩石加卸载蠕变曲线示意图^[16]如图2所示。

图2　岩石加卸载应变分离示意图

Fig. 2　Schematic diagram of strain separation of rock loading and unloading

根据图2及上述试验结果分析，可以将加载过程中总应变分为与时间无关的瞬时应变和与时间相关的蠕变应变，其中瞬时应变又可分为瞬时弹性应变和瞬时塑性应变，蠕变应变又可分为可随时间逐渐恢复的黏弹性应变和不可恢复的黏塑性应变，即

(1)

式中：为残余塑性应变。

基于岩石蠕变与卸载时应变恢复特征，将蠕变试验曲线根据式(1)进行分解，得到不同应力下各应变，如表3所示。

表3　岩石瞬时弹、塑性应变和黏弹、塑性应变试验结果

Table 3　Testing results of instantaneous elasto-plastic strain and viscoelasto-plastic strain

由表3可以看出：瞬时弹性应变、瞬时塑性应变、黏弹性应变和黏塑性应变都随应力的提高而提高。以试样L2为例，5种应力水平下瞬时塑性应变占瞬时应变的比例分别为33.11%，32.16%，31.28%，33.61%和34.50%，由于所占比例较大，因此，在建立描述瞬时应变的蠕变模型时，不能为了简化模型而忽略瞬时塑性应变，应将瞬时弹性应变和瞬时塑性应变分开考虑，分别建立相应的模型。此外，随偏应力逐级增大，5种应力下黏塑性应变占蠕变应变的比例分别为21.67%，28.70%，38.73%，42.74%和43.64%，这表明随偏应力增加，岩石不可恢复的黏塑性应变不断增加且逐渐累积储存于试样内部，当试验条件不能有效扼制住试样的黏塑性应变发展时就会发生破坏。在低应力下，岩石黏塑性应变所占比例较小，岩石的蠕变变形更大程度上表现为黏弹性变形；而在高应力下，岩石的蠕变变形表现为黏弹性和黏塑性共存的状态。

2  复合蠕变模型的建立

2.1　瞬时弹性体与瞬时塑性体

根据表1试验处理结果，绘制偏应力与瞬时弹性应变和瞬时塑性应变的关系曲线，结果如图3所示。

图3　瞬时弹、塑性应变与偏应力的关系曲线

Fig. 3　Relation curves among instantaneous elastic, plastic strain and deviator stress

由于应力加载阶段时间比后期应力维持加载时间短，因此，可认为瞬时弹性应变和瞬时塑性应变瞬间完成，两者量值只与应力有关而与时间无关。由图3可知：瞬时弹性应变、塑性应变与应力近似为线性增长关系，因此，可将建立瞬时弹性体与瞬时塑性体串联来描述，力学模型见图4。

图4　瞬时弹性体和塑性体模型

Fig. 4　Models of instantaneous elastic and plastic body

蠕变方程如下式所示：

(2)

式中：为应力；为瞬时弹性体和塑性体的瞬时弹性模量；为瞬时塑性体的瞬时弹性模量。

2.2　非线性黏弹性体

蠕变曲线由瞬时应变进入等速蠕变之前会经历衰减蠕变阶段，该阶段蠕变曲线斜率会逐渐减小，表现出非常明显的非线性特征，且衰减过程、特征及岩石进入等速蠕变的时间都与岩石应力密切相关。随蠕变试验进行，在常应力条件下岩石内部弱胶结物会逐渐发生破坏，宏观表现为蠕变应变随加载时间逐渐增加，在该阶段，黏性参数会随时间逐渐变大并在岩石进入等速蠕变阶段后逐渐趋于稳定值^[17]，说明蠕变已达稳定状态。传统Kelvin模型虽然能够较好地描述特定应力条件下蠕变衰减特征，但由于其黏性参数始终为不随时间变化的常数，因此，难以精确地反映不同应力下蠕变非线性特征。本文考虑Kelvin模型中黏性系数的时间相关性，假定该系数在蠕变过程中与时间满足幂函数关系^[18]，从而建立非线性黏弹性模型，如图5所示。

图5　非线性黏弹性体模型

Fig. 5　Model of nonlinear viscoelastic body

该模型微分型本构方程为

(3)

式中：为黏弹性应变；和分别为非线性黏弹性体的弹性模量和黏性系数；为常数；为时间。

解式(3)微分方程得到非线性黏弹性体的蠕变方程为

(4)

图6所示为不同对应的非线性黏弹性体蠕变曲线。由图6可知：越大，岩石由衰减蠕变进入等速蠕变所需的时间越短，衰减蠕变阶段蠕变斜率衰减的速度越快，说明该模型能够反映不同蠕变时间和不同应力状态下岩石衰减蠕变曲线特征，适用性较传统Kelvin模型更广泛。

图6　不同λ岩石蠕变曲线

Fig. 6　Creep curves of rock at different λ

假设在时刻卸载，式(3)可表示为

(5)

解式(5)得

(6)

式中：A为积分常数。

当时，，代入式(6)得到积分常数A为

(7)

把式(7)代入式(6)得到卸载阶段蠕变为

   (8)

2.3　黏塑性体

当试验应力介于某一临界应力和长期强度之间时，岩石等速蠕变阶段蠕变曲线表现出斜率大于0且随时间呈近似线性增长关系。该阶段岩石试样受长期蠕变特性及试样非均质性影响，在应力长时间作用下试样内部一些软弱结构处发生了局部破坏，此破坏由卸载后永久残余应变大于加载瞬时塑性应变所致。根据岩石流变理论，采用黏塑性体可描述该蠕变阶段岩石黏塑性应变随时间的变化关系，其力学模型如图7所示。

图7　黏塑性体模型

Fig. 7　Model of viscoplastic body

黏塑性体的蠕变方程为

(9)

式中：为黏塑性体的黏性系数。

2.4　复合蠕变模型

至此，建立了瞬时弹性体、瞬时塑性体、非线性黏弹性体和黏塑性体，分别描述岩石瞬时弹性应变、瞬时塑性应变、非线性黏弹性应变和黏塑性应变。根据叠加原理，将上述4个模型串联，建立可以同时描述加载过程中的上述4种应变、完全卸载时瞬时恢复的瞬时弹性应变和随时间逐渐恢复的黏弹性应变的复合蠕变模型，该复合蠕变力学模型如图8所示。

图8　复合蠕变模型

Fig. 8　Model of composite creep

当时，岩石总应变由瞬时弹性应变、瞬时塑性应变和非线性黏弹性应变组成，其本构方程为

(10)

当时，岩石总应变由瞬时弹性应变、瞬时塑性应变、非线性黏弹性应变和黏塑性应变组成，其本构方程为

(11)

根据拉普拉斯变换和逆变换，得到复合蠕变模型的蠕变方程为

(12)

由于卸载过程中瞬时塑性应变和黏塑性应变不可恢复，故复合蠕变模型的卸载方程为

(13)

3  三维蠕变本构方程的建立

由于实际岩体往往处于复杂的三维应力状态下，且本文引用数据为三轴蠕变数据，因此，为反映实际岩石受力情况以及更好地将模型理论结果与试验结果进行对比，需建立岩石在三维应力状态下复合蠕变模型的三维蠕变本构关系。假设岩石在三维应力状态下总应变为，瞬时弹、塑性应变分别为和，非线性黏弹性应变和黏塑性应变分别为和，根据叠加原理，有

(14)

根据岩石弹塑性理论得到岩石三维应力状态下瞬时弹性应变和瞬时塑性应变可表示为^[19]

(15)

式中：和分别为瞬时弹性体的剪切模量和体积模量；为瞬时塑性体的剪切模量；为偏应力张量；为平均应力；为Kronecker张量。

主应力空间下偏应力张量可表示为

(16)

非线性黏弹性体的三维蠕变方程可以参考一维应力状态下的蠕变方程：

(17)

式中：为非线性黏弹性体的剪切模量。

对于黏塑性体，其三维本构方程不仅与岩石屈服函数F和塑性势函数Q相关，而且与流动法则相关，并不是简单地将偏应力张量替换一维蠕变方程中的应力^[20]。根据流变力学理论得到黏塑性体三维蠕变方程为

(18)

式中：为开关函数；为屈服函数初始值，为简便计算一般取1^[21]；一般为指数函数或者幂函数，本文假设为幂函数且幂指数为1^[22]，即。此外，采用相关联流动法则，即^[23]，且可采用Mises屈服准则作为判定岩石屈服条件^[24]，因此，有

(19)

式中：为应力偏量第2不变量；为第1主应力。

在常规三轴压缩蠕变试验中，第2主应力，则

(20)

结合式(14)~(20)，得到复合蠕变模型三维蠕变方程为

(21)

卸载过程中的三维蠕变方程为

  (22)

4  模型验证

本文提出的复合蠕变模型参数包括瞬时弹性体的剪切模量和体积模量，该参数可以根据同等围压条件下常规三轴试验得到的弹性模量和泊松比确定；瞬时塑性体的剪切模量根据各偏应力下瞬时塑性应变结合式(15)确定。剩余参数，，和结合试验曲线根据最小二乘法确定。

以试样L7为例，根据式(21)得到复合蠕变模型各级偏应力下理论曲线与试验曲线的对比结果如图9所示。由图9可知：模型理论曲线与试验曲线吻合程度非常高，相关系数均在0.96以上，尤其在卸载过程，两者基本一致，由此说明本文建立的复合蠕变模型不仅能够描述石灰岩加载过程中非线性蠕变特征，而且能描述卸载过程中应变随时间的变化规律，证明了该模型的合理性。

图9　试样L7在不同偏应力水平下模型理论曲线与试验曲线的对比结果

Fig. 9　Comparisons of theoretical curves and experimental curves of L7 specimen at different deviator stresses

为进一步验证复合蠕变模型的适用性及可行性，以试样L2试验曲线为基础，采用复合蠕变模型理论曲线和广义Kelvin蠕变模型理论曲线分别与试验曲线进行对比，结果如图10所示。广义Kelvin蠕变模型三维轴向蠕变方程为

图10　试样L2在不同偏应力水平下试验曲线、复合蠕变模型拟合曲线和广义Kelvin模型拟合曲线对比结果

Fig. 10　Comparison results of test curve, composite creep model fitting curve and generalized Kelvin model fitting curve of sample L2 at different deviator stresses

(23)

式中：和分别为广义Kelvin模型中弹性体的剪切模量和体积模量；和分别为广义Kelvin模型中黏弹性体的剪切模量和黏性系数。

相应的三维卸载蠕变方程为

(24)

由图10可知，广义Kelvin蠕变模型和复合蠕变模型都能很好地描述加载过程中石灰岩的蠕变特征，但在衰减蠕变阶段，复合蠕变模型描述效果比广义Kelvin蠕变模型的效果好。同时，复合蠕变模型卸载蠕变曲线基本上与试验曲线保持一致，而广义Kelvin模型曲线明显比试验曲线低，这是由于广义Kelvin蠕变模型不能将瞬时塑性应变从瞬时应变中分离，也不能将黏塑性应变从蠕变应变中分离，而是认为加载过程中产生的瞬时应变在卸载的一瞬间完全恢复，从而使得模型理论曲线与试验曲线产生较大误差。而复合蠕变模型综合考虑了瞬时弹性应变、瞬时塑性应变、非线性黏弹性应变和黏塑性应变，认为卸载过程中瞬时弹性应变瞬间恢复，非线性黏弹性应变随蠕变时间逐渐恢复，而瞬时塑性应变和黏塑性应变为永久残余应变，因此，复合蠕变模型卸载过程中模型曲线与试验曲线高度吻合，相关系数均大于0.968，对比结果进一步验证本文建立的能够同时描述石灰岩分级加卸载蠕变特征的复合蠕变模型是合理、可行的。

石灰岩试样L2和L7在不同偏应力下复合蠕变模型参数见表4。

表4　各级偏应力下复合蠕变模型参数

Table 4　Parameters of composite creep model under different deviator stresses

5  结论

1) 本文根据石灰岩循环加卸载蠕变曲线特征，将瞬时弹性应变、瞬时塑性应变、非线性黏弹性应变和黏塑性应变从总应变中分离，采用瞬时弹性体、瞬时塑性体、非线性黏弹性体和黏塑性体对上述4种应变分别进行描述，再将四者串联，建立了一种能够同时描述石灰岩加、卸载过程蠕变特征的复合蠕变模型。

2) 根据叠加原理和非线性流变理论推导了岩石在一维、三维应力状态下的加卸载蠕变本构方程。石灰岩不同偏应力下加卸载蠕变试验曲线和复合蠕变模型理论曲线的对比结果表明两者高度吻合，该模型不仅能够描述石灰岩加载过程中非线性蠕变特征，而且能描述卸载过程中应变随时间的变化规律。

3) 用复合蠕变模型理论曲线和广义Kelvin蠕变模型理论曲线分别与试验曲线进行对比，结果表明广义Kelvin蠕变模型和复合蠕变模型都能很好地描述加载过程中石灰岩的蠕变特征，但在衰减蠕变阶段及卸载蠕变阶段，复合蠕变模型描述效果要优于广义Kelvin蠕变模型。复合蠕变模型综合考虑了瞬时弹性应变、瞬时塑性应变、非线性黏弹性应变和黏塑性应变，在卸载过程中模型曲线与试验曲线高度吻合，相关系数均大于0.968，对比结果验证了本文建立的复合蠕变模型是合理、可行的。

参考文献：

[1] 刘新喜, 李盛南, 周炎明, 等. 高应力泥质粉砂岩蠕变特性及长期强度研究[J]. 岩石力学与工程学报, 2020, 39(1): 138-146.

LIU Xinxi, LI Shengnan, ZHOU Yanming, et al. Study on creep behavior and long-term strength of argillaceous siltstone under high stresses[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(1): 138-146.

[2] 张亮亮, 王晓健, 周瑞鹤. 一种新的岩石非线性黏弹塑性蠕变模型研究[J]. 力学季刊, 2020, 41(1): 116-124.

ZHANG Liangliang, WANG Xiaojian, ZHOU Ruihe. A new nonlinear viscoelastic-plastic creep model for rocks[J]. Chinese Quarterly of Mechanics, 2020, 41(1): 116-124.

[3] 韩阳, 谭跃虎, 李二兵, 等. 岩石非定常Burgers蠕变模型及其参数识别[J].工程力学, 2018, 35(3): 210-217.

HAN Yang, TAN Yuehu, LI Erbing, et al. Non-stationary Burgers creep model of rock and its parameter identification[J]. Engineering Mechanics, 2018, 35(3): 210-217.

[4] 陈卫忠, 谭贤君, 吕森鹏, 等. 深部软岩大型三轴压缩流变试验及本构模型研究[J]. 岩石力学与工程学报, 2009, 28(9): 1735-1744.

CHEN Weizhong, TAN Xianjun, LU Senpeng, et al. Research on large-scale triaxial compressive rheological test of soft rock in depth and its constitutive model[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(9): 1735-1744.

[5] HU Bo, YANG Shengqi, XU Peng. A nonlinear rheological damage model of hard rock[J]. Journal of Central South University, 2018, 25(7): 1665-1677.

[6] 刘开云, 薛永涛, 周辉. 参数非定常的软岩非线性黏弹塑性蠕变模型[J]. 中国矿业大学学报, 2018, 47(4): 921-928.

LIU Kaiyun, XUE Yongtao, ZHOU Hui. A nonlinear viscoelastic-plastic creep model of soft rock with unsteady parameters[J]. Journal of China University of Mining & Technology, 2018, 47(4): 921-928.

[7] 刘东燕, 谢林杰, 庹晓峰, 等. 不同围压作用下砂岩蠕变特性及非线性黏弹塑性模型研究[J]. 岩石力学与工程学报, 2017, 36(S2): 3705-3712.

LIU Dongyan, XIE Linjie, TUO Xiaofeng, et al. Creep properties of sandstone under different confining pressures and research on a nonlinear viscoelasto-plastic creep model[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(S2): 3705-3712.

[8] ROBERTS L A, BUCHHOLZ S A, MELLEGARD K D, et al. Cyclic loading effects on the creep and dilation of salt rock[J]. Rock Mechanics and Rock Engineering, 2015, 48(6): 2581-2590.

[9] SINGH A, KUMAR C, KANNAN L G, et al. Estimation of creep parameters of rock salt from uniaxial compression tests[J]. International Journal of Rock Mechanics and Mining Sciences, 2018, 107: 243-248.

[10] BOUKHAROV G N, CHANDA M W, BOUKHAROV N G. The three processes of brittle crystalline rock creep[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1995, 32(4): 325-335.

[11] 李娜, 曹平, 衣永亮, 等. 分级加卸载下深部岩石流变实验及模型[J]. 中南大学学报(自然科学版), 2011, 42(11): 3465-3471.

LI Na, CAO Ping, YI Yongliang, et al. Creep properties experiment and model of deep rock with step loading and unloading[J]. Journal of Central South University(Science and Technology), 2011, 42(11): 3465-3471.

[12] 张治亮, 徐卫亚, 王伟. 向家坝水电站坝基挤压带岩石三轴蠕变试验及非线性黏弹塑性蠕变模型研究[J]. 岩石力学与工程学报, 2011, 30(1): 132-140.

ZHANG Zhiliang, XU Weiya, WANG Wei. Study of triaxial creep tests and its nonlinear visco-elastoplastic creep model of rock from compressive zone of dam foundation in xiangjiaba hydropower station[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(1): 132-140.

[13] 赵延林, 曹平, 文有道, 等. 岩石弹黏塑性流变试验和非线性流变模型研究[J]. 岩石力学与工程学报, 2008, 27(3): 477-486.

ZHAO Yanlin, CAO Ping, WEN Youdao, et al. Elastovisco-plastic rheological experiment and nonlinear rheological model of rocks[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(3): 477-486.

[14] 赵延林, 唐劲舟, 付成成, 等. 岩石黏弹塑性应变分离的流变试验与蠕变损伤模型[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.

[15] ZHAO Yanlin, ZHANG Lianyang, WANG Weijun, et al. Separation of elastoviscoplastic strains of rock and a nonlinear creep model[J]. International Journal of Geomechanics, 2018, 18(1): 04017129.

[16] 徐鹏, 杨圣奇. 循环加卸载下煤的黏弹塑性蠕变本构关系研究[J]. 岩石力学与工程学报, 2015, 34(3): 537-545.

XU Peng, YANG Shengqi. Study of visco-elasto-plastic constitutive model of coal under cyclic loading[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(3): 537-545.

[17] 阎岩, 王思敬, 王恩志. 基于西原模型的变参数蠕变方程[J]. 岩土力学, 2010, 31(10): 3025-3035.

YAN Yan, WANG Sijing, WANG Enzhi. Creep equation of variable parameters based on Nishihara model[J]. Rock and Soil Mechanics, 2010, 31(10): 3025-3035.

[18] 张亮亮, 王晓健. 岩石黏弹塑性损伤蠕变模型研究[J]. 岩土工程学报, 2020, 42(6): 1085-1092.

ZHANG Liangliang, WANG Xiaojian. Viscoelastic-plastic damage creep model for rock[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(6): 1085-1092.

[19] 张亮亮, 王晓健. 一种新的岩石非线性损伤蠕变模型[J]. 应用力学学报, 2020, 37(1): 372-377.

ZHANG Liangliang, WANG Xiaojian. A new nonlinear damage creep model of rocks[J]. Chinese Journal of Applied Mechanics, 2020, 37(1): 372-377.

[20] 齐亚静, 姜清辉, 王志俭, 等. 改进西原模型的三维蠕变本构方程及其参数辨识[J]. 岩石力学与工程学报, 2012, 31(2): 347-355.

QI Yajing, JIANG Qinghui, WANG Zhijian, et al. 3D creep constitutive equation of modified Nishihara model and its parameters identification[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(2): 347-355.

[21] 姜鹏, 潘鹏志, 赵善坤, 等. 基于应变能的岩石黏弹塑性损伤耦合蠕变本构模型及应用[J]. 煤炭学报, 2018, 43(11): 2967-2979.

JIANG Peng, PAN Pengzhi, ZHAO Shankun, et al. A coupled elasto-viscoplastic damage model based on strain energy theory of rock and application[J]. Journal of China Coal Society, 2018, 43(11): 2967-2979.

[22] ABU AL-RUB R K, DARABI M K, KIM S M, et al. Mechanistic-based constitutive modeling of oxidative aging in aging-susceptible materials and its effect on the damage potential of asphalt concrete[J]. Construction and Building Materials, 2013, 41(4): 439-454.

[23] 潘晓明, 杨钊, 许建聪. 非定常西原黏弹塑性流变模型的应用研究[J]. 岩石力学与工程学报, 2011, 30(S1): 2640-2646.

PAN Xiaoming, YANG Zhao, XU Jiancong. Application study of nonstationary nishihara viscoelasto-plastic rheological model[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(S1): 2640-2646.

[24] 单仁亮, 宋立伟, 李东阳, 等. 冻结红砂岩非线性蠕变模型的研究[J]. 岩土力学, 2014, 35(6): 1541-1546.

SHAN Renliang, SONG Liwei, LI Dongyang, et al. Study of nonlinear creep model of frozen red sandstone[J]. Rock and Soil Mechanics, 2014, 35(6): 1541-1546.

(编辑  刘锦伟)

收稿日期： 2020 -07 -06; 修回日期： 2020 -09-15

基金项目(Foundation item)：国家自然科学基金资助项目(51874005) (Project(51874005) supported by the National Natural Science Foundation of China)

通信作者：张亮亮，博士研究生，从事地下结构设计与支护技术研究；E-mail: zllaust@163.com

引用格式： 张亮亮, 王晓健. 考虑黏弹塑性应变分离的岩石复合蠕变模型研究[J]. 中南大学学报(自然科学版), 2021, 52(5): 1655-1665.

Citation: ZHANG Liangliang, WANG Xiaojian. Study on composite creep model of rock considering viscoelastic plastic strain separation[J]. Journal of Central South University(Science and Technology), 2021, 52(5): 1655-1665.