王 唯，劉遠明

(貴州大學(xué)土木工程學(xué)院，貴州貴陽550025)

0 引 言

在諸多巖土工程實際中，廣泛涉及巖石蠕變問題。因此，巖石蠕變方面的研究具有非常重要的理論意義。近年來，隨著技術(shù)手段和理論的創(chuàng)新，巖石流變領(lǐng)域的發(fā)展取得了豐碩的成果。許多研究學(xué)者通過將簡單的彈性、粘性及塑性元件串并聯(lián)組合成復(fù)雜的模型來描述流變現(xiàn)象。由于受到加載水平的限制，能描述巖石蠕變?nèi)^程的模型并不多，尤其是加速蠕變階段。學(xué)者們對模型參數(shù)進行了改進，取得良好的研究成果[1-2]。夏才初等[3]基于基本流變力學(xué)模型進行串聯(lián)組合，形成15個流變力學(xué)模型，并提出辨識模型的方法；朱昌星等[4]將能夠描述加速流變特性的時效損傷變量引入巖石非線性蠕變損傷模型中，并結(jié)合板巖剪切流變試驗驗證了模型的合理性；楊圣奇等[5]對龍灘水電站泥板巖進行剪切流變試驗，將非線性流變元件與西原模型串聯(lián)起來，建立了新的能夠描述加速流變特性的巖石非線性流變模型；趙延林等[6]對橄欖巖試樣進行流變試驗，提出由瞬彈性Hooke體、塑性體、改進Bingham體串聯(lián)形成的非線性粘塑性流變模型，能反應(yīng)不同應(yīng)力水平下巖石流變規(guī)律；孫鈞院士[7]對非線性流變問題的求解提供了2種方案：一種是由試驗擬合的經(jīng)驗本構(gòu)關(guān)系式，另一種是在粘塑性部分內(nèi)加上一項非線性的經(jīng)驗粘性元件作為對線性流變模型的修正；唐佳[8]等對白云石英片巖和綠泥鈉長片巖進行蠕變試驗，基于應(yīng)力-應(yīng)變曲線，提出在Mohr-Coulomb準(zhǔn)則基礎(chǔ)上引入新的M-C塑性元件，構(gòu)成了改進的Burgers蠕變模型。

綜上，巖石蠕變本構(gòu)模型建立過程中，通過將模型中的線性元件替換為非線性元件，即可得到能夠描述巖石的加速蠕變過程。本文在已有的流變理論基礎(chǔ)上，建立了一種新的對非線性巖石流變本構(gòu)模型。

1 非線性流變模型的建立

1.1 巖石的蠕變曲線

巖石的蠕變曲線圖1[9]。圖中，3條蠕變曲線是在不同的應(yīng)力條件下得到的。其中，σa>σb>σc。當(dāng)巖石在低應(yīng)力水平的恒定荷載持續(xù)作用下時，應(yīng)變量雖然會隨時間有所增加，但是蠕變變形的速率會逐漸減小，最后會趨于一個穩(wěn)定的數(shù)值。當(dāng)應(yīng)力水平較大時，曲線不在穩(wěn)定于一個定值，而是無限增加直到破壞。根據(jù)應(yīng)變速率的不同，將蠕變過程劃分為3個階段：第一蠕變階段(曲線中ab段)，應(yīng)變速率隨時間的增加而減小，稱為減速蠕變階段；第二蠕變階段(曲線中bc段)應(yīng)變速率保持不變，稱為等速蠕變階段；第三蠕變階段(曲線中cd段)，應(yīng)變速率迅速增加直到巖石破壞，稱為加速蠕變破壞。

圖1 巖石蠕變曲線

由于傳統(tǒng)的元件基本都是線性的，無論通過何種串并聯(lián)組合方式，只能描述其減速和等速蠕變。因此，通過引入非線性粘滯系數(shù)的牛頓體，與應(yīng)力閥值并聯(lián)組合新元件來描述加速蠕變的過程。

1.2 非線性Bingham模型的改進

傳統(tǒng)的Bingham模型由1個塑性元件與牛頓粘壺并聯(lián)后再與1個Hooke體串聯(lián)。蠕變方程為

(1)

式中，ε為蠕變量；E和η為蠕變系數(shù)；t為蠕變時間；σs為應(yīng)力閥值。傳統(tǒng)的Bingham模型見圖2。

圖2 傳統(tǒng)的Bingham模型

當(dāng)作用應(yīng)力水平時，傳統(tǒng)的Bingham的蠕變量隨時間呈線性關(guān)系。然而在工程實際中，巖石蠕變通常表現(xiàn)為非線性的函數(shù)關(guān)系；而傳統(tǒng)的Bingham模型只能描述穩(wěn)態(tài)蠕變，并不能描述整個巖石的蠕變過程。因此，通過將Binham模型中的某些線性元件替換為非線元件，可建立描述巖石蠕變?nèi)^程的改進的Bingham模型。

為描述巖石的加速蠕變與衰減蠕變過程，曹樹剛等[10]提出粘滯系數(shù)是時間的函數(shù)關(guān)系；宋飛等[11]提出粘滯系數(shù)是應(yīng)力水平與時間的函數(shù)關(guān)系；周家文等[12]基于流變伺服儀得到綠片巖三軸試驗數(shù)據(jù)，構(gòu)造出能夠模擬巖石衰減蠕變和加速蠕變特征的非線性函數(shù)；馬林建等[13]對鹽巖的三軸壓縮試驗表明，鹽巖的彈性模量隨加載的時間呈指數(shù)遞減趨勢。本文基于流變力學(xué)理論，假定在靜力荷載作用下，該模型的彈性元件系數(shù)與時間t呈衰減趨勢，衰減規(guī)律為

E2(α)=E2/[1-exp(-αt)]

(2)

式中，E2為Bingham模型彈性元件初值，該參數(shù)主要反映巖石的瞬時彈性特性；α為控制彈性系數(shù)衰減的無量綱參數(shù)，可由試驗測定。

則Hooke體的一維方程為

(3)

當(dāng)σ=σ0時，蠕變方程為

(4)

對于描述巖石加速蠕變階段，不少文獻中組合模型的粘滯體不再是常值，而是用非線性牛頓體來代替。因此，本文從粘滯系數(shù)入手，用冪函數(shù)形式的粘塑性體來描述巖石加速蠕變隨時間的變化關(guān)系。非線性牛頓體的一維方程為

(5)

(6)

本文基于傳統(tǒng)的Binham模型，將模型中的常值牛頓粘體替換為非線性的牛頓體，并串聯(lián)1個描述巖石彈性衰減的Hooke體，建立了如圖3所示的改進的Binham模型。圖中，通過引入開關(guān)應(yīng)力閥值σs來表示巖石材料加速蠕變的初始應(yīng)力，只有當(dāng)σ>σs時，并聯(lián)的粘塑性體才會發(fā)生作用。其中，原模型彈性體描述巖石的瞬時彈性變形，非線性彈性體描述巖石的衰減蠕變，粘塑性體描述巖石的加速蠕變。

圖3 改進的Bingham流變模型

根據(jù)流變力學(xué)相關(guān)理論，串聯(lián)組合各元件上應(yīng)力相等，總應(yīng)變等于各元件應(yīng)變和。可得基于圖3的狀態(tài)方程如下

σ=σ1=σ2=σ3

(7)

ε=ε1+ε2+ε3

(8)

σ1=E1ε1

(9)

(10)

(11)

(1) 當(dāng)σ≤σs時，該模型退化為2個彈性體串聯(lián)，本構(gòu)方程為

(12)

(2)當(dāng)σ>σs時，有

(13)

1.3 一維蠕變方程

令σ=σ0H(t)，代入式(12)、(13)本構(gòu)方程中，方程兩邊經(jīng)過Laplace變換以及逆變換可得到蠕變方程。

(1)當(dāng)σ≤σs時，有

(14)

(2)當(dāng)σ>σs時，有

(15)

1.4 三維蠕變方程

在三維應(yīng)力狀態(tài)下，流變方程很難用形象具體的元件去描述。參照彈性力學(xué)方法，將巖石內(nèi)部1點的應(yīng)力張量σij分解為球應(yīng)力張量σm和偏應(yīng)力張量Sij這2個部分，巖石的三軸應(yīng)力狀態(tài)可由3個不變的主應(yīng)力σ1、σ2、σ3表示，即

(16)

(17)

同理，也可將應(yīng)變張量εij分解為球應(yīng)變張量εm和偏應(yīng)變張量eij，即

εij=eij+δijεm

(18)

對于大多數(shù)各向同性材料而言，一般條件下均存在以下規(guī)律[14]：球應(yīng)變張量只與球應(yīng)力張量有關(guān)，而偏應(yīng)變張量只與偏應(yīng)力張量有關(guān)，這給建立三維本構(gòu)方程帶來很大方便。對于一般的線彈性材料而言，球應(yīng)變張量與偏應(yīng)變張量分別與球應(yīng)力張量與偏應(yīng)力張量成正比。對于彈塑性材料，在塑性范圍內(nèi)，球應(yīng)變張量與球應(yīng)力張量成正比，偏應(yīng)力張量與偏應(yīng)變張量的關(guān)系比彈性材料要復(fù)雜。但是為了尋找6個獨立應(yīng)變分量與6個獨立應(yīng)力分量的關(guān)系，可簡化為偏應(yīng)變張量與偏應(yīng)力張量之間的關(guān)系。根據(jù)彈性力學(xué)理論，可直接寫出三維條件下流微分本構(gòu)形式為

(19)

式中，eij為偏應(yīng)變張量；σii為應(yīng)力；δii為應(yīng)變；K為體積模量；G為剪切模量。

綜上分析，對于簡單的元件組合模型，可采用單軸蠕變試驗參數(shù)確定其參數(shù)。孫鈞[15]認(rèn)為，在三維應(yīng)力狀態(tài)下，蠕變模型很難用形象具體的元件去描述，但三維蠕變本構(gòu)方程可以采用類比法通過一維蠕變本構(gòu)方程推導(dǎo)得出，從一維模型到三維模型，只需要將下列的符號對應(yīng)：σ→Sij、ε→eij、E→G、β→η，則非線性Bingham模型的三維形式為

(20)

exp(tλ)(Sij>σs)

(21)

(22)

2 模型的適用性驗證

2.1 巖石蠕變試驗曲線

為驗證模型的適用性，楊彩虹等[16]采用遼寧工程技術(shù)大學(xué)能源與新技術(shù)研究所自主研制的三軸滲透儀進行了巖石蠕變試驗。天然頁巖單軸壓縮蠕變試驗曲線見圖4。從圖4可以看出，應(yīng)力水平較低且小于38 MPa時，該模型只發(fā)生衰減蠕變和等速蠕變；應(yīng)力水平大于38 MPa且作用的時間足夠長時，該模型會發(fā)生衰減蠕變、等速蠕變以及加速蠕變。因此，該試驗?zāi)軌蚰M巖石完整蠕變過程，且當(dāng)應(yīng)力水平越大時，進入到加速蠕變階段所需時間越短。

圖4 不同應(yīng)力水平下天然頁巖蠕變

2.2 本構(gòu)模型的參數(shù)確定

1stOpt為目前世界領(lǐng)先的非線性曲線擬合及綜合優(yōu)化分析計算軟件，該軟件在曲線擬合、非線性復(fù)雜模型參數(shù)估算求解等領(lǐng)域居世界領(lǐng)先地位。通用全局優(yōu)化算法(UGO)是1stOpt內(nèi)置的一種線性優(yōu)化控制算法，該算法最大特點是克服了當(dāng)今世界上在優(yōu)化計算領(lǐng)域中使用迭代法必須給出合適初始值的難題，即用戶勿需給出參數(shù)初始值，而由1stOpt隨機給出，通過其獨特的全局優(yōu)化算法最終找出最優(yōu)解。本以天然頁巖蠕變試驗為依據(jù)，給出了求解粘彈塑性模型參數(shù)的方法[17]。

(1)在低應(yīng)力狀態(tài)下，即σ≤σs，此時的本構(gòu)模型退化為2個Hooke體串聯(lián)，蠕變方程為式(14)。具體的求解方法是：先求解瞬時彈性模量E1，當(dāng)t=0時，E1=σ0/ε，其中，ε為天然頁巖的瞬時形變量，由圖4蠕變曲線得出，由于σ0已知即可求得；再通過1stopt軟件，通過對試驗數(shù)據(jù)的擬合求得彈性模量E2以及蠕變參數(shù)α。

(2)在高應(yīng)力狀態(tài)下，即σ>σs，單軸壓縮擬合蠕變方程為式(15)，三軸壓縮擬合蠕變方程為式(21)。通過輸入數(shù)據(jù)并繪制散點圖，自定義函數(shù)進行非線性擬合。運用全局優(yōu)化算法對提出改進后模型蠕變方程的參數(shù)進行求取以及模型的適用性驗證。程序運行完成即可得到粘滯系數(shù)η1與蠕變指數(shù)λ，5個流變參數(shù)全部求得。各流變參數(shù)值見表1。從表1中可以看出，在不同的應(yīng)力水平下試驗數(shù)據(jù)與擬合曲線的相關(guān)系數(shù)均達到0.9以上。因此，該模型能夠較好地描述巖石蠕變的3個階段，證明了該模型的適用性。

表1 非線性Bingham模型巖石蠕變參數(shù)

蔣海飛等采用RLW-2000巖石三軸流變試驗系統(tǒng)，對取自重慶某深基坑細砂巖進行三軸蠕變試驗。該巖石的單軸抗壓強度為63.40 MPa。蠕變試驗設(shè)計有2種加載工況，見表2。選取1號試件在高于屈服應(yīng)力水平下的相關(guān)試驗和數(shù)據(jù)，對本文提出的本構(gòu)方程式(21)進行擬合和模型參數(shù)的求取。三軸條件下的砂巖應(yīng)變曲線的擬合參數(shù)見表3。擬合曲線見圖5。

表2 三軸蠕變試驗工況

注：1號試件自然干燥。

2.3 非線性參數(shù)影響分析

從非線性模型的擬合曲線選取3種蠕變曲線代表曲線(見圖6)以及表1的參數(shù)反演可得，在應(yīng)力水平未達到開關(guān)閥值，且當(dāng)α越小時，在相同時間內(nèi)，蠕變變形量和蠕變速率越大，而且進入到穩(wěn)態(tài)蠕變階段的時間延長，同時衰減蠕變的曲率半徑越大；當(dāng)應(yīng)力水平達到開關(guān)閥值時，通過分析41、45、55 MPa與65 MPa數(shù)據(jù)曲線可知，隨著λ的增大，巖石經(jīng)歷加速蠕變階段的時間縮短，且在相同時間內(nèi)，蠕變變形量和蠕變速率也越大。λ的值與巖石非線性加速有很大關(guān)系，這也說明了本文考慮這2個系數(shù)對巖石蠕變過程的必要性與合理性。

表3 砂巖應(yīng)變曲線的擬合參數(shù)

圖5 擬合曲線

圖6 不同應(yīng)力狀態(tài)下擬合曲線

比較頁巖和砂巖在加速蠕變階段的擬合值發(fā)現(xiàn)，砂巖的蠕變指數(shù)λ較頁巖高1個數(shù)量級，可能參數(shù)λ的大小與巖石種類以及受到的荷載有關(guān)，具體的變化規(guī)律需進一步研究確定。

從圖5、6可知，試驗數(shù)據(jù)與擬合曲線兩者吻合效果良好，特別在高應(yīng)力狀態(tài)下擬合的相關(guān)系數(shù)R2接近于1，因此基于改進的Binghan模型所建立的本構(gòu)方程能較好地描述靜力荷載下巖石蠕變規(guī)律。

3 結(jié) 語

本文基于傳統(tǒng)的Bingham模型，將其模型中的常值牛頓粘體替換為與時間相關(guān)的非常值元件并串聯(lián)一個非線性彈性元件，建立了基于靜力荷載下非線性Bingham模型，得出下以結(jié)論 ：

(1)與傳統(tǒng)的Bingham模型相比，改進后的非線性模型能夠較好地模擬巖石的衰減蠕變與加速蠕變階段，且該模型參數(shù)較少，模型簡單方便用于工程實際。

(2)基于改進后的Bingham模型所推導(dǎo)的流變本構(gòu)方程對頁巖與砂巖的試驗數(shù)據(jù)擬合表明，所建立的本構(gòu)方程能較好地描述巖石蠕變整個階段，尤其是加速蠕變階段，擬合曲線與試驗數(shù)據(jù)的相關(guān)系數(shù)R2均在0.9以上。因此，與傳統(tǒng)的流變模型對比，本文模型更適合描述巖石材料的流變性質(zhì)。

(3)通過對控制衰減蠕變系數(shù)α與加速蠕變系數(shù)λ進行參數(shù)分析，α隨應(yīng)力水平增大而減小，且α越小，衰減蠕變的時間越長；λ隨應(yīng)力水平增大而增大，且λ越大，巖石在加速蠕變階段的蠕變速度越大。