付宏淵 戚雙星 史振寧 高乾豐 匡波 曾鈴

DOI：?10.11835/j.issn.2096-6717.2022.028

基金項目：國家自然科學(xué)基金（51838001、51908069、52078066、52078067）；湖湘青年科技創(chuàng)新人才項目（2020RC306）；湖南省重點領(lǐng)域研發(fā)項目（2019SK2171）；長沙市杰出創(chuàng)新青年培養(yǎng)計劃（kq1905043）；國家級大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練項目（202120536003）；長沙理工大學(xué)公路養(yǎng)護技術(shù)國家工程實驗室開放基金（kfj170103）

作者簡介：付宏淵（1965-?），男，教授，博士生導(dǎo)師，主要從事巖土工程防災(zāi)減災(zāi)研究，E-mail： fhy001@126.com。通信作者：曾鈴（通信作者），男，教授，博士生導(dǎo)師，E-mail：?zl001@csust.edu.cn。

Received： 2021?12?01

Foundation items： National Natural Science Foundation of China （Nos. 51838001， 51908069， 52078066， 52078067）， Youth Scientifc and Technological Innovation Talents of Hunan Province （No. 2020RC306）， Research and Development Projects in Key Fields of Hunan Province （No. 2019SK2171）， Outstanding Innovative Youth Training Program of Changsha City （No. kq1905043）， College Students Innovation and Entrepreneurship Program of China （No. 202120536003）， Open Fund of National Engineering Laboratory of Highway Maintenance Technology， Changsha University of Science & Technology （No. kfj170103）

Author brief： FU Hongyuan （1965-?）， professor， doctorial supervisor， main research interest： disaster prevention and mitigation in geotechnical engineering， E-mail： fhy001@126.com.

corresponding author：ZENG Ling （corresponding author），?professor，?doctorial supervisor，?E-mail：?zl001@csust.edu.cn.

（1. 長沙理工大學(xué)?a.土木工程學(xué)院；?b. 公路養(yǎng)護技術(shù)國家工程實驗室；?c. 交通運輸工程學(xué)院，長沙?410114；?2. 廣西交通設(shè)計集團有限公司，?南寧530029）

摘要：針對邊坡、隧道工程中粉砂質(zhì)泥巖長期穩(wěn)定性問題，采用DZSZ-150型多場耦合巖石三軸壓力試驗機，利用粉砂質(zhì)泥巖相似材料試樣，開展不同圍壓下粉砂質(zhì)泥巖蠕變試驗。根據(jù)蠕變試驗結(jié)果建立可描述粉砂質(zhì)泥巖蠕變?nèi)^程的非線性黏彈塑性本構(gòu)模型，以7 MPa圍壓下蠕變試驗結(jié)果為例，對所提出的蠕變本構(gòu)模型進行合理性和可行性驗證。結(jié)果表明：試樣的穩(wěn)態(tài)蠕變速率與偏應(yīng)力大小有關(guān)，隨著偏應(yīng)力的增加，試樣穩(wěn)態(tài)蠕變速率與偏應(yīng)力之間符合指數(shù)函數(shù)關(guān)系；在加速蠕變階段的蠕變速率相對低偏應(yīng)力水平下減速蠕變和穩(wěn)態(tài)蠕變速率明顯要增大很多，且非常短時間內(nèi)就發(fā)生了蠕變破壞；提出一種描述粉砂質(zhì)泥巖全蠕變過程的非線性黏彈塑性本構(gòu)模型，該模型能描述粉砂質(zhì)泥巖蠕變-破壞全過程的蠕變變形規(guī)律，反映了粉砂質(zhì)泥巖在不同偏應(yīng)力水平作用下的蠕變特性，克服了經(jīng)典模型中無法準確描述加速蠕變過程的難題。

關(guān)鍵詞：粉砂質(zhì)泥巖；蠕變試驗；加速蠕變；蠕變特性；本構(gòu)模型

中圖分類號：TU411.6 ????文獻標志碼：A ????文章編號：2096-6717（2024）03-0001-08

Creep characteristics of silty mudstone and nonlinear viscoelastic-plasticity constitutive model

FU Hongyuan^1a，1b，?QI Shuangxing^1a，?SHI Zhenning^1b，1c，?GAO Qianfeng^1b，1c，?KUANG Bo²，?ZENG Ling^1a

（1a. School of Civil Engineering；?1b. National Engineering Laboratory of Highway Maintenance Technology；?1c. School of Traffic & Transportation Engineering， Changsha University of Science & Technology， Changsha 410114， P. R. China；?2. Guangxi Communications Design Group Co.， Ltd.， Nanning 530029， P. R. China）

Abstract： In order to guarantee the long-term stability of silty mudstone in slope and tunneling projects， the DZSZ-150 multi-field coupling rock triaxial pressure testing machine was used to conduct creep tests on silty mudstone specimens of similar materials under different surrounding pressures. Based on the testing results， a nonlinear viscoelastic-plastic constitutive model describing the whole creep process of silty mudstone was established， and the proposed creep model was verified for its rationality and feasibility by taking the results at a confining pressure of 7 MPa as an example. It show that the steady-state creep rate of the specimen is related to the magnitude of the deviatoric ?stress， and with increase of the deviatoric stress， the steady-state creep rate and the deviatoric stress conform to an exponential function relationship; the creep rate in the accelerated creep stage is significantly greater than that in the decelerated creep and steady-state creep at low stress levels， and the creep damage occurs in a very short time; a nonlinear viscoelastic-plastic constitutive model is proposed to describe the whole creep process of the silty mudstone， which can describe the creep-deformation law of the whole process of creep-damage of silty mudstone， reflecting the creep characteristics of silty mudstone under different levels of deviatoric ?stress， and avoiding the problem that the classical model cannot accurately describe the accelerated creep process.

Keywords： silty mudstone；?creep test；?accelerated creep；?creep characteristics；?constitutive model

隨著基礎(chǔ)設(shè)施建設(shè)的快速發(fā)展，越來越多的粉砂質(zhì)泥巖工程問題不斷涌現(xiàn)，尤其是切割粉砂質(zhì)泥巖巖體形成的高陡邊坡、地下洞室等工程巖體的蠕變性強，長期服役過程中會產(chǎn)生蠕變變形，導(dǎo)致支護結(jié)構(gòu)破壞，防護體系失穩(wěn)等災(zāi)害頻發(fā)^[1-2]。粉砂質(zhì)泥巖蠕變問題已成為威脅大型基礎(chǔ)設(shè)施安全運營的關(guān)鍵難題。故迫切需要開展粉砂質(zhì)泥巖蠕變試驗，研究粉砂質(zhì)泥巖長期蠕變規(guī)律，建立合理的蠕變本構(gòu)模型。

學(xué)者們在對不同類型巖石的蠕變行為及蠕變本構(gòu)模型方面開展了大量研究^[3-4]。在研究巖石的蠕變行為方面，辛亞軍等^[5]采用紅砂巖開展了不同加載速率下三軸蠕變試驗，分析了蠕變變形特征，結(jié)果表明，隨著加載速率增大，減速蠕變應(yīng)變增量越大，等速蠕變應(yīng)變增量越小。胡波等^[6]在進行粉砂巖三軸蠕變試驗時發(fā)現(xiàn)，粉砂巖的起始蠕變速率會隨著應(yīng)力的增加而增大，蠕變進入穩(wěn)態(tài)階段后蠕變速率也會隨時間增加。Wang等^[7]研究了周期循環(huán)荷載下鹽巖的蠕變行為，發(fā)現(xiàn)當(dāng)超出循環(huán)周期數(shù)臨界值后，周期越短，蠕變發(fā)生不可逆變形越大。在研究巖石的蠕變本構(gòu)模型方面，王永巖等^[8]以西原模型為基礎(chǔ)，引入了可考慮溫度和圍壓的軟巖蠕變參數(shù)，同時將黏塑性元件替換成非線性黏塑性元件，得到了可描述溫度和圍壓共同影響下的非線性蠕變本構(gòu)模型，并用該模型對三軸蠕變試驗結(jié)果進行了擬合驗證。劉家順等^[9]針對含水條件下弱膠結(jié)軟巖的蠕變問題，通過引入分數(shù)階微積分理論，建立并驗證了可用于描述全蠕變過程的分數(shù)階弱膠結(jié)軟巖蠕變本構(gòu)模型。Liu等^[10]在對鹽巖蠕變試驗數(shù)據(jù)分析的基礎(chǔ)上，引入了一個損傷系數(shù)來描述加速蠕變階段的非線性應(yīng)變，得到的新模型與實驗數(shù)據(jù)擬合效果非常好，可為鹽巖工程的蠕變性行為研究提供理論依據(jù)。

綜上所述，現(xiàn)有對巖石蠕變的研究已經(jīng)取得了較豐碩的成果，對認識不同類型巖石蠕變特性有著重要意義。但對泥巖蠕變特性，尤其對于粉砂質(zhì)泥巖蠕變特性及本構(gòu)模型的系統(tǒng)研究還較少。一直以來，粉砂質(zhì)泥巖強度低、易破碎、制樣困難，并且開展長期蠕變試驗對設(shè)備要求高，導(dǎo)致粉砂質(zhì)泥巖的長期蠕變力學(xué)特性的試驗研究鮮有涉及，尚未建立可描述粉砂質(zhì)泥巖非線性蠕變特征的本構(gòu)模型。現(xiàn)有研究成果不能準確表述粉砂質(zhì)泥巖工程巖體的蠕變規(guī)律。為此，開展蠕變試驗，分析粉砂質(zhì)泥巖應(yīng)變與時間關(guān)系，建立基于加速蠕變改進的非線性黏彈塑性本構(gòu)模型，并驗證該模型的合理性和可行性，以期為粉砂質(zhì)泥巖工程的長期穩(wěn)定性分析研究奠定基礎(chǔ)。

1 巖石試樣蠕變試驗

1.1 試樣制備

鑒于粉砂質(zhì)泥巖內(nèi)部裂隙較多且分布不均，在取樣過程中難以獲取均一的粉砂質(zhì)泥巖原樣試樣，為降低試驗結(jié)果的離散性，采用相似材料代替粉砂質(zhì)泥巖開展室內(nèi)試驗。已有研究結(jié)果已經(jīng)證明，所配制相似材料試樣的強度及變形規(guī)律可與粉砂質(zhì)泥巖試樣的規(guī)律相一致^[11-12]。試驗所采用的相似材料制作方法見文獻[13]，基本物理參數(shù)如表1所示。

1.2 試驗設(shè)備

試驗采用DZSZ-150型多場耦合巖石三軸壓力試驗機（圖1）。該試驗機系統(tǒng)主要包括圍壓加載系統(tǒng)、軸壓加載系統(tǒng)、變形測量系統(tǒng)和控制系統(tǒng)，可實現(xiàn)蠕變試驗過程中對試樣的應(yīng)力、應(yīng)變及時間數(shù)據(jù)的連續(xù)采集。圍壓和軸向荷載控制范圍分別為0～150 MPa、0～50 kN，精度分別為±0.000 1 MPa、±0.001 kN。變形測量極值為25 mm，精度為±0.001 mm。

1.3 試驗方案

根據(jù)國際巖石力學(xué)學(xué)會（IRSM）推薦的試驗規(guī)程^[14]，制作試樣尺寸均為標準50 mm×h100 mm的圓柱體。在開展三軸蠕變試驗時，采用加圍壓到設(shè)定值后保持圍壓不變，再施加每一級荷載，直到發(fā)生蠕變破壞，具體應(yīng)力路徑和加載方式見圖2。根據(jù)現(xiàn)有研究成果^[15-17]，第1級荷載可采用單軸抗壓強度的30%～50%，然后再依次增加每級荷載，共設(shè)置6～8級軸向應(yīng)力，直到試樣發(fā)生蠕變破壞后停止。根據(jù)試樣的單軸抗壓強度^[13]，蠕變試驗的第一級采用5 MPa保持偏應(yīng)力不變，持續(xù)時間為24 h，再依次施加下一級偏應(yīng)力荷載，每一級增加2.5 MPa，持續(xù)時間為24 h，直至試樣發(fā)生加速蠕變破壞；在此過程中圍壓分別為3、5、7、9 MPa，始終保持穩(wěn)定不變。

2 試驗結(jié)果分析

2.1 蠕變變形分析

不同圍壓下的三軸蠕變特征曲線如圖3所示。從整個蠕變特征曲線可看出，除最終發(fā)生失穩(wěn)破壞的階段外，每一級偏應(yīng)力荷載下試樣都會經(jīng)歷一段減速蠕變后逐漸趨于穩(wěn)定。隨著加載級數(shù)的不斷增加，試樣蠕變變形進入非穩(wěn)態(tài)蠕變區(qū)，最終發(fā)生了加速蠕變。此時施加的偏應(yīng)力已超過試樣的屈服強度，試樣很快從減速蠕變過渡到穩(wěn)態(tài)蠕變階段，應(yīng)變將持續(xù)增加，進入加速蠕變階段，最終發(fā)生了失穩(wěn)破壞。由圖3可知，除圍壓3 MPa外，加載偏應(yīng)力等級共經(jīng)歷了7級，前6級每級持續(xù)時間為1 440 min，第7個等級持續(xù)時間存在一定差異，最終發(fā)生加速蠕變破壞。具體有如下特點：1）當(dāng)試樣處在應(yīng)力低于屈服強度時，減速蠕變現(xiàn)象不明顯，蠕變曲線都經(jīng)歷減速蠕變過渡到穩(wěn)態(tài)蠕變，最終蠕變趨于穩(wěn)定。2）當(dāng)偏應(yīng)力荷載超過屈服強度后，即蠕變試驗進入最后一級，試樣先減速蠕變，然后穩(wěn)態(tài)蠕變，最后快速進入加速蠕變，最終發(fā)生了蠕變破壞。

2.2 穩(wěn)態(tài)蠕變速率特性分析

在蠕變試驗中，蠕變速率是蠕變特征的重要指標之一^[18]。根據(jù)蠕變試驗數(shù)據(jù)，重點分析穩(wěn)態(tài)蠕變速率和偏應(yīng)力的關(guān)系。圖4為穩(wěn)態(tài)蠕變速率與偏應(yīng)力擬合的關(guān)系曲線，由圖4可知：試樣的穩(wěn)態(tài)蠕變速率與偏應(yīng)力有關(guān)，在不同圍壓下，隨著偏應(yīng)力的增加，穩(wěn)態(tài)蠕變速率與偏應(yīng)力之間符合指數(shù)函數(shù)關(guān)系

，R²≥0.922 8 （1）

式中：為穩(wěn)態(tài)蠕變速率；A和B都是大于零的參數(shù)；為偏應(yīng)力。其中，不同偏應(yīng)力等級下穩(wěn)態(tài)蠕變速率的曲線所對應(yīng)的參數(shù)A和B變化關(guān)系如表2所示。

2.3 軸向應(yīng)變速率分析

為進一步分析蠕變試驗過程中軸向蠕變速率隨時間的變化規(guī)律，將軸向應(yīng)變對加載時間求導(dǎo)，所得結(jié)果為試驗實測的軸向應(yīng)變速率。在7 MPa圍壓下，每級荷載下試樣軸向應(yīng)變速率-時間曲線如圖5所示。由圖5可知：OA段，軸向蠕變速率隨時間增加而減小，并很快降低，即初期蠕變階段；AB段，軸向蠕變速率保持不變，并趨于恒定值，即穩(wěn)定蠕變階段；此外，初期蠕變階段，在相同時刻下，不同等級應(yīng)力水平下軸向初始蠕變速率隨應(yīng)力水平等級的增加而逐漸增大，且表現(xiàn)出良好的一致性。

圖6為7 MPa圍壓下，在第7級荷載（偏應(yīng)力為20 MPa）下試樣進入非穩(wěn)態(tài)蠕變階段過程中軸向應(yīng)變和軸向應(yīng)變率隨時間的變化曲線。由圖6可知，隨時間的增加，軸向應(yīng)變速率在不斷下降，先經(jīng)歷了t₁=100 min的減速蠕變階段，蠕變速率降到非常??；然后過渡到穩(wěn)態(tài)蠕變階段，這一過程持續(xù)到t₂=360 min，軸向應(yīng)變速率接近定值，軸向應(yīng)變接近平緩的直線；最后進入加速蠕變階段，持續(xù)約1 min，試樣迅速發(fā)生加速蠕變破壞。

3 粉砂質(zhì)泥巖試樣蠕變模型研究

結(jié)合圖6的分析可知，在不高于屈服強度的加載等級（）時，對蠕變試驗曲線特征進行分析，可得出試樣的蠕變表現(xiàn)出明顯的黏彈性特征；在高于屈服強度的加載等級（）時，試樣等速蠕變階段會持續(xù)很短時間后進入加速蠕變破壞階段，表現(xiàn)出明顯的非線性黏彈塑性特征。

目前，描述巖石黏彈性蠕變特性的元件模型有許多種，最常用的有三元件的廣義Kelvin模型、四元件的伯格斯模型等^[19-20]。根據(jù)上述分析可知，現(xiàn)有模型中只有賓漢姆模型和西原模型符合，考慮到賓漢姆模型在低于破壞應(yīng)力水平下，理想黏塑性體沒有變形，只有彈性變形，沒有蠕變變形，不符合本文粉砂質(zhì)泥巖試樣的蠕變特點^[21]。而現(xiàn)有的經(jīng)典西原模型，僅能夠描述試樣除加速蠕變階段的總體變化趨勢。

針對上述本構(gòu)模型存在的缺陷，在經(jīng)典西原模型的基礎(chǔ)上，采用現(xiàn)有線性黏彈性模型（Modified Kelvin模型，簡稱MK模型）與非線性黏塑性模型^[17]（簡稱NVPB模型）串聯(lián)組合，建立非線性黏彈塑性本構(gòu)模型來表征粉砂質(zhì)泥巖的全蠕變過程，在該模型中假定非線性黏塑性模型與時間呈指數(shù)函數(shù)關(guān)系^[22]，其蠕變力學(xué)模型示意圖如圖7所示。

如圖7所示，構(gòu)建的新蠕變模型可根據(jù)疊加原理得到一維蠕變方程，考慮到應(yīng)力水平對試驗蠕變特性的影響，故蠕變方程要分成兩種情況。

1）當(dāng)時，模型退化為廣義Kelvin三元件模型，塑性模型不發(fā)揮其作用，該一維蠕變方程為 （2）

式中：為應(yīng)變；為應(yīng)力；為巖石屈服強度；彈性體模量；為Kelvin彈性模量；為Kelvin黏滯系數(shù)；為時間。

2）當(dāng)時，試樣發(fā)生加速蠕變，此時該模型的黏塑性部分將發(fā)揮作用，則該一維蠕變方程為 （3）

式中：為非線性損傷黏塑性元件黏滯系數(shù)；為蠕變參數(shù)。

大量蠕變研究成果表明^[23]：在蠕變過程中，當(dāng)巖石所受荷載超過屈服應(yīng)力的閾值時，應(yīng)力長時間作用下巖石內(nèi)部產(chǎn)生蠕變損傷。因此，引入損傷變量D隨時間累積呈負指數(shù)函數(shù)形式來描述黏性系數(shù)的蠕變損傷劣化，構(gòu)建了NVPB模型。根據(jù)以往研究成果^[22]，軟巖發(fā)生蠕變過程中由于損傷變量D隨時間t累積之間滿足式（4）所示關(guān)系，即損傷變量D的表達式為 （4）

式中：D為損傷變量；為與巖石材料性質(zhì)相關(guān)受加載應(yīng)力的影響的系數(shù)，其主要影響加速蠕變發(fā)生的時間；為時間。

引入開關(guān)函數(shù)將式（2）和式（3）統(tǒng)一為 （5）

其中，開關(guān)函數(shù)。

4 本構(gòu)模型的驗證和分析

為了驗證新蠕變本構(gòu)模型的合理性和可行性，根據(jù)試驗數(shù)據(jù)，采用Levenberg-Marquardt（L-M）算法對7 MPa圍壓下非線性黏彈塑性本構(gòu)模型進行參數(shù)識別，結(jié)果見表3。

圖8為時，擬合的理論曲線和試驗曲線的對比結(jié)果。由圖8對比曲線可知，試驗值和擬合值吻合效果比較理想。表明蠕變本構(gòu)模型采用新建的非線性黏彈塑性本構(gòu)模型可很好地反映粉砂質(zhì)泥巖的減速蠕變階段和穩(wěn)態(tài)蠕變階段的蠕變特征。

圖9為時，加速蠕變階段擬合的理論值曲線和試驗值曲線的對比結(jié)果。由圖9通過基于三軸蠕變試驗值與西原模型、伯格斯模型和本文提出的模型擬合值進行對比分析發(fā)現(xiàn)，加速蠕變階段下建立的蠕變本構(gòu)模型精度明顯高于西原模型和伯格斯模型。即模型擬合的相關(guān)性系數(shù)R²為0.96，明顯優(yōu)于西原模型（R²=0.33）和伯格斯模型（R²=0.95）的結(jié)果，與試驗值曲線最為接近。

為進一步分析不同蠕變模型與蠕變試驗結(jié)果的精度，以試驗值為依據(jù)，對上述3種蠕變模型每個小時的蠕變應(yīng)變與試驗值進行對比，得到的誤差見表4、誤差絕對值如圖10所示，需要說明的是，0 h的數(shù)據(jù)為第2 s采集的數(shù)據(jù)，也就是在試驗開始后第1次采集的數(shù)據(jù)，即時間為第2 s。由表4和圖10可知，相比西原模型和伯格斯模型得到的蠕變應(yīng)變值，建立的粉砂質(zhì)泥巖非線性黏彈塑性本構(gòu)模型每小時的應(yīng)變值與試驗值的誤差最小。說明建立的粉砂質(zhì)泥巖非線性黏彈塑性本構(gòu)模型誤差明顯更小，更適用于對粉砂質(zhì)泥巖蠕變特性中的加速蠕變階段規(guī)律的描述。

綜合圖8～圖10可知，對以7 MPa圍壓為例的粉砂質(zhì)泥巖試樣全過程蠕變試驗結(jié)果進行擬合驗證后可知，本文建立的非線性黏彈塑性蠕變本構(gòu)模型克服了經(jīng)典西原模型和伯格斯模型不能精確地對加速蠕變過程進行描述的難題，更適用于準確描述粉砂質(zhì)泥巖蠕變?nèi)^程。

為進一步揭示粉砂質(zhì)泥巖蠕變破壞機理，圖11給出了粉砂質(zhì)泥巖試樣蠕變破壞過程的機理示意圖。圖11（a）所示為初始狀態(tài)下的結(jié)構(gòu)示意圖，內(nèi)部巖石骨架顆粒周圍包裹著一定量的膠結(jié)物，但是膠結(jié)物與其膠結(jié)程度較低；同時巖石骨架內(nèi)會存在相對均勻分布的初始原生閉合孔隙。圖11（b）所示為等速蠕變階段的示意圖，在不高于屈服強度的加載等級（）時，試樣內(nèi)部主要以巖石骨架顆粒彈性變形為主，原生閉合孔隙被不斷擠壓，密度增大，此時的作用力不足以造成巖石內(nèi)部產(chǎn)生新缺陷，在該級偏應(yīng)力荷載下巖石以結(jié)構(gòu)強度為主。圖11（c）所示為加速蠕變階段的示意圖，隨著偏應(yīng)力荷載的不斷增加，彈性應(yīng)變引起的總應(yīng)變減少，造成巖石內(nèi)部結(jié)構(gòu)以非線性變形為主，試樣內(nèi)部部分孔隙結(jié)構(gòu)會發(fā)生破壞貫穿，匯集產(chǎn)生微裂紋，甚至可能產(chǎn)生裂隙簇，骨架結(jié)構(gòu)重新分布，損傷不斷累積，承載力會不斷下降，但是巖石仍具有一定的承載力，且應(yīng)變增量非常緩慢，以塑性變形為主。圖11（d）所示為破壞階段的示意圖，偏應(yīng)力荷載繼續(xù)增加，當(dāng)高于屈服強度的加載等級（）時，上述現(xiàn)象會不斷加劇，最終導(dǎo)致巖石內(nèi)部骨架顆粒發(fā)生崩潰，裂紋連接產(chǎn)生斷裂面，喪失承載能力，發(fā)生該現(xiàn)象的根本原因是，累積損傷致使粉砂質(zhì)泥巖承載能力下降到臨界破壞荷載，迅速發(fā)生加速蠕變。綜上可知，粉砂質(zhì)泥巖發(fā)生蠕變破壞的整個過程是從初始蠕變階段到相對緩慢的穩(wěn)態(tài)蠕變階段，再到突增的加速蠕變階段，這與文獻[24-25]中提到深部巖體的蠕變損傷機理相吻合。

對于粉砂質(zhì)泥巖工程而言，可采用FLAC^3D軟件提供的二次開發(fā)程序接口，實現(xiàn)提出的非線性黏彈塑性本構(gòu)模型的二次開發(fā)，采用開發(fā)后的本構(gòu)模型進行數(shù)值模擬，并結(jié)合實際巖體工程的長期監(jiān)測數(shù)據(jù)，對粉砂質(zhì)泥巖相關(guān)工程的長期變形進行預(yù)測，為該類巖體長期穩(wěn)定性研究及防護提供理論指導(dǎo)。

5 結(jié)論

1）三軸蠕變試驗研究表明，試樣發(fā)生臨界破壞時表現(xiàn)出明顯的加速蠕變。

2）在粉砂質(zhì)泥巖蠕變試驗過程中，試樣的穩(wěn)態(tài)蠕變速率與偏應(yīng)力大小有關(guān)，在不同圍壓下，隨著偏應(yīng)力的增加，穩(wěn)態(tài)蠕變速率與偏應(yīng)力之間符合指數(shù)函數(shù)關(guān)系。

3）在蠕變進入破壞階段時，在加速蠕變階段的蠕變速率相對低應(yīng)力水平下減速蠕變和穩(wěn)態(tài)蠕變速率明顯要增大很多，且非常短時間內(nèi)就發(fā)生了蠕變破壞。

4）結(jié)合三軸蠕變試驗結(jié)果，基于西原模型和伯格斯模型的不足，提出了一種可描述粉砂質(zhì)泥巖全蠕變過程的非線性黏彈塑性本構(gòu)模型。

5）對非線性黏彈塑性本構(gòu)模型的合理性和可行性進行驗證表明，該模型不僅可描述低應(yīng)力水平過程中的減速蠕變階段和穩(wěn)態(tài)蠕變階段，還能很好地描述高應(yīng)力下的加速蠕變階段，有效克服了經(jīng)典蠕變模型存在的不足。

參考文獻

[1] ?付宏淵， 陳鏡丞， 曾鈴， 等. 溫、濕度對粉砂質(zhì)泥巖單軸力學(xué)性能的影響試驗[J]. 土木工程學(xué)報， 2019， 52（1）： 89-98.FU H Y， CHEN J C， ZENG L， et al. Experiment on the effects of temperature and humidity on uniaxial mechanical properties of silty mudstone [J]. China Civil Engineering Journal， 2019， 52（1）： 89-98. （in Chinese）

[2] ?ZHANG S， XU Q， HU Z M. Effects of rainwater softening on red mudstone of deep-seated landslide， Southwest China [J]. Engineering Geology， 2016， 204： 1-13.

[3] ?YANG S Q， HU B. Creep and long-term permeability of a red sandstone subjected to cyclic loading after thermal treatments [J]. Rock Mechanics and Rock Engineering， 2018， 51（10）： 2981-3004.

[4] ?王游， 趙翔. 改進的Bingham巖石蠕變模型研究[J]. 力學(xué)季刊， 2021， 42（1）： 168-177.WANG Y， ZHAO X. Studies on improved Bingham rock creep model [J]. Chinese Quarterly of Mechanics， 2021， 42（1）： 168-177. （in Chinese）

[5] ?辛亞軍， 郝海春， 呂鑫， 等. 高應(yīng)力區(qū)峰后紅砂巖蠕變特性及模型分析[J]. 實驗力學(xué)， 2019， 34（5）： 809-818.XIN Y J， HAO H C， LV X， et al. Creep characteristics and model analysis of red sandstone in post-peak high stress zone [J]. Journal of Experimental Mechanics， 2019， 34（5）： 809-818. （in Chinese）

[6] ?胡波， 王宗林， 梁冰， 等. 巖石蠕變特性試驗研究[J]. 實驗力學(xué)， 2015， 30（4）： 438-446.HU B， WANG Z L， LIANG B， et al. Experimental study of rock creep properties [J]. Journal of Experimental Mechanics， 2015， 30（4）： 438-446. （in Chinese）

[7] ?WANG J B， ZHANG Q， SONG Z P， et al. Experimental study on creep properties of salt rock under long-period cyclic loading [J]. International Journal of Fatigue， 2021， 143： 106009.

[8] ?王永巖， 王鴻偉， 崔立樁， 等. 溫度-圍壓共同作用的軟巖蠕變模型及試驗驗證[J]. 應(yīng)用力學(xué)學(xué)報， 2020， 37（6）： 2561-2566， 2705.WANG Y Y， WANG H W， CUI L Z， et al. Investigation and experimental verification of a soft rock creep model at the effect of temperature and confining pressure [J]. Chinese Journal of Applied Mechanics， 2020， 37（6）： 2561-2566， 2705. （in Chinese）

[9] ?劉家順， 靖洪文， 孟波， 等. 含水條件下弱膠結(jié)軟巖蠕變特性及分數(shù)階蠕變模型研究[J]. 巖土力學(xué)， 2020， 41（8）： 2609-2618.LIU J S， JING H W， MENG B， et al. Research on the effect of moisture content on the creep behavior of weakly cemented soft rock and its fractional-order model [J]. Rock and Soil Mechanics， 2020， 41（8）： 2609-2618. （in Chinese）

[10] ?LIU J， WU F， ZOU Q L， et al. A variable-order fractional derivative creep constitutive model of salt rock based on the damage effect [J]. Geomechanics and Geophysics for Geo-Energy and Geo-Resources， 2021， 7（2）： 1-16.

[11] ?FU H Y， QI S X， SHI Z N， et al. Effect of nano-CaCO₃on the physical and mechanical properties of analogue to silty mudstone materials [J]. Arabian Journal of Geosciences， 2021， 14（23）： 1-11.

[12] ?WANG J J， SHI Z N， ZENG L， et al. The effects of different nanoadditives on the physical and mechanical properties of similar silty mudstone materials [J]. Advances in Civil Engineering， 2020， 2020： 8850436.

[13] ?FU H Y， QI S X， SHI Z N， et al. Mixing ratios and cementing mechanism of similar silty mudstone materials for model tests [J]. Advances in Civil Engineering， 2021， 2021： 2426130.

[14] ?HOEK E，WOOD D，SHAH S.A modified Hoek-Brown failure criterion for jointed rock masses[C]//Hudson J A.Proceedings of the Rock Characterization， Symposium of ISRM. London： Thomas Telford Publishing， 1992：209-214.

[15] ?劉東燕. 巖石力學(xué)[M]. 重慶： 重慶大學(xué)出版社， 2014.LIU D Y. Rock mechanics [M]. Chongqing： Chongqing University Press， 2014. （in Chinese）

[16] ?張亮亮， 王曉健， 周瑞鶴. 一種新的巖石非線性黏彈塑性蠕變模型研究[J]. 力學(xué)季刊， 2020， 41（1）： 116-124.ZHANG L L， WANG X J， ZHOU R H. A new nonlinear viscoelastic-plastic creep model for rocks [J]. Chinese Quarterly of Mechanics， 2020， 41（1）： 116-124. （in Chinese）

[17] ?魏堯， 楊更社， 申艷軍， 等. 白堊系飽和凍結(jié)砂巖蠕變試驗及本構(gòu)模型研究[J]. 巖土力學(xué)， 2020， 41（8）： 2636-2646.WEI Y， YANG G S， SHEN Y J， et al. Creep test and constitutive model of Cretaceous saturated frozen sandstone [J]. Rock and Soil Mechanics， 2020， 41（8）： 2636-2646. （in Chinese）

[18] ?李安潤， 鄧輝， 王小雪， 等. 飽水-失水循環(huán)條件下紅層泥巖蠕變特性及本構(gòu)模型研究[J]. 工程地質(zhì)學(xué)報， 2021， 29（3）： 843-850.LI A R， DENG H， WANG X X， et al. Research on creep characteristics and constitutive model of red bed mudstone under saturated-dehydrated cycle [J]. Journal of Engineering Geology， 2021， 29（3）： 843-850. （in Chinese）

[19] ?PU S Y， ZHU Z D， SONG L， et al. Fractional-order visco-elastoplastic constitutive model for rock under cyclic loading [J]. Arabian Journal of Geosciences， 2020， 13（9）： 1-11.

[20] ?杜超， 楊春和， 馬洪嶺， 等. 深部鹽巖蠕變特性研究[J]. 巖土力學(xué)， 2012， 33（8）： 2451-2456， 2520.DU C， YANG C H， MA H L， et al. Study of creep characteristics of deep rock salt [J]. Rock and Soil Mechanics， 2012， 33（8）： 2451-2456， 2520. （in Chinese）

[21] ?蔡美峰， 何滿潮， 劉東燕. 巖石力學(xué)與工程[M]. 北京： 科學(xué)出版社， 2002.CAI M F， HE M C， LIU D Y. Rock mechanics and engineering [M]. Beijing： Science Press， 2002. （in Chinese）

[22] ?王春萍， 陳亮， 梁家瑋， 等. 考慮溫度影響的花崗巖蠕變?nèi)^程本構(gòu)模型研究[J]. 巖土力學(xué)， 2014， 35（9）： 2493-2500， 2506.WANG C P， CHEN L， LIANG J W， et al. Creep constitutive model for full creep process of granite considering thermal effect [J]. Rock and Soil Mechanics， 2014， 35（9）： 2493-2500， 2506. （in Chinese）

[23] ?陳陸望， 李圣杰， 陳逸飛， 等. 巖石含水蠕變損傷模型的開發(fā)與應(yīng)用[J]. 固體力學(xué)學(xué)報， 2018， 39（6）： 642-651.CHEN L W， LI S J， CHEN Y F， et al. Further development and application of a creep damage model for water-bearing rocks [J]. Chinese Journal of Solid Mechanics， 2018， 39（6）： 642-651. （in Chinese）

[24] ?張晉勛， 宋永威， 楊昊， 等. 凍結(jié)飽水砂卵石蠕變損傷特性試驗研究[J]. 科學(xué)技術(shù)與工程， 2021， 21（29）： 12665-12671.ZHANG J X， SONG Y W， YANG H， et al. Experimental research on the creep damage characteristics of frozen saturated sandy gravel [J]. Science Technology and Engineering， 2021， 21（29）： 12665-12671. （in Chinese）

[25] ?邵珠山， 靳冬冬， 陳浩哲， 等. 含水狀態(tài)對石英砂巖單軸分級蠕變性能影響研究[J]. 應(yīng)用力學(xué)學(xué)報， 2021， 38（5）： 1839-1845.SHAO Z S， JIN D D， CHEN H Z， et al. Investigation on creep properties of different moisture state quartz sandstone by multi-stage test [J]. Chinese Journal of Applied Mechanics， 2021， 38（5）： 1839-1845. （in Chinese）

（編輯??胡英奎）