巖石材料蠕變模型中應(yīng)力時間分布規(guī)律的研究

巖石材料蠕變模型中應(yīng)力時間分布規(guī)律的研究

1對非線性流變本構(gòu)關(guān)系的改進(jìn)變形是巖石的一個重要力學(xué)特征之一，包括變形、應(yīng)力松弛、長期強(qiáng)度等。其中,蠕變特指恒定應(yīng)力作用下應(yīng)變隨時間持續(xù)增長的現(xiàn)象,因其與工程實踐聯(lián)系緊密而成為流變研究領(lǐng)域中最重要的內(nèi)容。目前,巖石流變的理論研究多數(shù)基于元件理論,即按照巖石的彈性、塑性和黏滯性設(shè)置基本元件(Hooke體、St.Venant體、Newton體)并建立由其串聯(lián)或并聯(lián)組合而成的模型,目前使用較多的有Kelvin模型、Burgers模型、西原模型等。其中,西原模型(見圖1)因能夠比較清楚且全面地反映流變中的彈-黏彈-黏塑性而受到廣大研究者的重視。元件理論所用的基本元件都是線性的,無論構(gòu)成模型的元件怎么多、組合方式怎么復(fù)雜,最終得到的本構(gòu)關(guān)系還是線性的,不能研究非線性流變,因此,需要對其改進(jìn)以建立非線性流變模型。目前,廣泛使用的方法是用非線性元件來代替模型中的線性元件,比如鄧榮貴等認(rèn)為,由于結(jié)構(gòu)的特殊性,巖石的流變具有非牛頓流體的特性,進(jìn)而提出一種非線性黏滯元件。曹樹剛等認(rèn)為,巖石的黏滯系數(shù)隨裂隙的閉合與擴(kuò)展有先增大再變小的變化趨勢,基于此得到新的黏滯體模型。陳沅江等提出了一種新的非線性變截面黏滯元件。徐衛(wèi)亞等提出了一個新的非線性黏滯元件,對黏性元件的本構(gòu)關(guān)系進(jìn)行改進(jìn)以模擬加速蠕變曲線。用非線性元件(目前多為黏滯元件)代替線性元件的方法能夠體現(xiàn)出巖石流變的非線性特征,尤其是可以模擬出加速蠕變現(xiàn)象,但元件理論中的元件都是一維的,在推導(dǎo)三維本構(gòu)關(guān)系時,目前采用類似于彈性理論中將一維Hooke定律推廣到三維的方法,需要作以下假定:(1)材料的體積變形在受力瞬時完成,且體積變形為彈性,不隨時間產(chǎn)生變化;(2)只有偏應(yīng)力張量引起蠕變,球應(yīng)力張量作用下材料不產(chǎn)生蠕變;(3)流變過程中,泊松比不隨時間發(fā)生變化。以上3個假定比較適用于金屬材料或彈性階段,而對于巖石材料而言,與其蠕變試驗現(xiàn)象矛盾。比如,巖石的體積變形貫穿整個蠕變階段,隨時間不僅有壓縮還會發(fā)生明顯的擴(kuò)容現(xiàn)象,而且不同的體積變形還會引起不同方向的蠕變。試驗表明,蠕變泊松比與瞬時泊松比不同,且泊松比在蠕變過程中并非定值,比如在擴(kuò)容階段,側(cè)向應(yīng)變增加的趨勢要超過軸向,泊松比會增大。因此,如何從一維模型推導(dǎo)符合巖石非線性流變特性的三維模型,仍有待進(jìn)一步研究。除了采用非線性元件求解非線性流變問題外,還可以采用變參數(shù)的方法。孫鈞院士提出:將黏滯系數(shù)η視為非定常的變數(shù)值,由試驗確定,再進(jìn)行非線性流變本構(gòu)關(guān)系的分析計算。這種處理方案比較理想,學(xué)術(shù)理念也較為嚴(yán)格,但計算處理則比較繁復(fù)和困難。采用非定常的流變模型代替?zhèn)鹘y(tǒng)上定常的流變模型能更加準(zhǔn)確地預(yù)測工程巖體的時效非線性變形特征,將是下一步深化研究的一個努力方向。按照此思路,筆者基于西原模型,通過蠕變試驗結(jié)果,對采用變參數(shù)的蠕變方程進(jìn)行了初步探討,希望為巖石非線性流變研究提供新的研究思路與研究方法。2雜材料在地質(zhì)構(gòu)造中的力學(xué)參數(shù)隨時間的變化在目前的流變問題分析中,一般都把流變參數(shù)視為常數(shù),但實際上,巖石這種復(fù)雜材料在地質(zhì)構(gòu)造運(yùn)動、地下水滲流和自然風(fēng)化等諸多因素的作用下,其某些力學(xué)參數(shù)隨時間而變化是十分明顯的。比如,試驗已證實,彈性模量、強(qiáng)度等參數(shù)通常都會隨時間的增長而降低。因此,將流變參數(shù)視為變量符合巖石材料的本質(zhì)屬性。2.1非線性黏滯體模型關(guān)于變參數(shù)流變模型,現(xiàn)有的研究成果主要有:呂學(xué)清基于工程巖體受環(huán)境影響這一事實提出了巖體工程變參數(shù)計算法的思路,假設(shè)彈性模量E(t)與黏滯系數(shù)η(t)按直線規(guī)律或指數(shù)規(guī)律衰減,以Maxwell模型為例介紹了變參數(shù)流變計算的方法;宋德彰、孫鈞在S.S.Vyalov的研究基礎(chǔ)上提出了非線性的Bingham模型,認(rèn)為黏滯系數(shù)η是所施加的荷載及荷載持續(xù)作用時間的函數(shù)η(τ,t);許宏發(fā)根據(jù)軟巖的單軸壓縮蠕變試驗討論了軟巖強(qiáng)度和彈性模量的時間效應(yīng),認(rèn)為軟巖的強(qiáng)度和彈性模量都是時間的函數(shù),其變化規(guī)律均隨時間的延長而降低;丁志坤、呂愛鐘在頁巖蠕變試驗的基礎(chǔ)上,驗證了引入非定常參數(shù)的必要性,認(rèn)為其能比定常模型更準(zhǔn)確地反映巖石的黏彈性變形性能。作者認(rèn)為,H-K蠕變模型中,Kelvin體的彈性模量Ek為非定常參數(shù),而黏滯系數(shù)ηk與Hooke體的彈性模量EH為定值。張貴科、徐衛(wèi)亞在分析巖體與流變模型變形特點的基礎(chǔ)上,認(rèn)為黏滯系數(shù)η在低應(yīng)力作用下隨時間不斷增大,當(dāng)應(yīng)力超過長期強(qiáng)度時η不斷減小。在此認(rèn)識的基礎(chǔ)上,提出了一個非線性黏滯體模型。孫鈞提出在巖石的非線性流變中,黏滯系數(shù)η是應(yīng)力水平σ、應(yīng)力-應(yīng)變狀態(tài)Ω和流變時間t的非線性函數(shù),此時η(σ,Ω,t)是一個變量。2.2西原模型的三維壓縮變量(1)開展非線性流變研究的初衷是解決線性模型無法描述加速蠕變的缺點,但流變的非線性并非僅存在于加速階段,而是貫穿于整個蠕變過程。(2)流變的非線性并非僅體現(xiàn)在黏滯系數(shù)的變化上,模型中各參數(shù)都應(yīng)視為變量。對于西原模型的三維蠕變方程(式(1))而言,由于泊松比v是變量,根據(jù)G=E/[2(1+v)],則瞬時剪切模量G0與黏彈性模量G1也是變量,體積模量K表征著時間的體積變形,顯然也是變量。兩個黏滯系數(shù)η1、η2關(guān)系著非線性流變特性,也應(yīng)是變量。(3)引起蠕變參數(shù)變化的原因包括應(yīng)力、時間、外部條件(滲流、溫度、化學(xué)因素等),這些因素引起巖石微觀層面的變化。蠕變參數(shù)與上述各因素之間的相互關(guān)系可由相應(yīng)的試驗成果得出。(4)目前的變參數(shù)研究多數(shù)為預(yù)先定性地確定黏滯系數(shù)等蠕變參數(shù)的變化規(guī)律再構(gòu)建表達(dá)式,而從試驗結(jié)果得到蠕變參數(shù)的變化規(guī)律后再構(gòu)建蠕變方程則更加客觀。2.3應(yīng)力對各測定儀器膠凝劑的影響根據(jù)分級加載的壓縮蠕變試驗成果,可以得到巖石試件在不同應(yīng)力作用下的蠕變曲線。在同一應(yīng)力作用下,蠕變參數(shù)只與時間相關(guān),取t=t1,t2,…,ti,則在各時刻附近的小區(qū)間內(nèi),可近似認(rèn)為蠕變參數(shù)不變。在此區(qū)間內(nèi)總可以找到n組試驗數(shù)據(jù),從而可通過非線性最小二乘擬合法求出ti時刻的各蠕變參數(shù)。按此方法,可以得到不同時刻的蠕變參數(shù),并分析參數(shù)與時間t之間的關(guān)系。同理,可以得到相同時刻在不同應(yīng)力作用下的各個蠕變參數(shù)。因此,可以分析參數(shù)與應(yīng)力σ之間的關(guān)系。非線性最小二乘法的原理與過程限于篇幅這里不再詳述,請參見文獻(xiàn)。3參數(shù)與時間和電壓之間的關(guān)系3.1流變力學(xué)分析圖2展示了筆者所做雙軸應(yīng)力作用下的巖石滲流-蠕變耦合試驗,試驗方案見表1。試驗得到試件3個方向的蠕變曲線見圖3,其中εz為垂直方向應(yīng)變,εx為水平方向應(yīng)變,εy為滲流方向的應(yīng)變。按照流變力學(xué)研究的慣例,應(yīng)變均以壓縮為正,擴(kuò)大為負(fù)。將3個方向的應(yīng)變相加,可以得到試件的體積應(yīng)變εv,以體積減小為正,擴(kuò)大為負(fù),見圖4。3.2不同參數(shù)試驗研究由于水平壓力與滲透水壓固定,應(yīng)變只隨垂向應(yīng)力與時間變化,因此,可以根據(jù)試驗結(jié)果研究蠕變參數(shù)與σz和t之間的關(guān)系。本文以εz為研究對象,取t=5、15、24、48、…、168h,如圖5所示。對于三維西原蠕變方程的5個參數(shù),體積模量K與瞬時變形模量G0可以通過試驗結(jié)果直接求出,而黏彈性變形模量G1、黏彈性黏滯系數(shù)η1、黏塑性黏滯系數(shù)η2需要由試驗結(jié)果通過非線性最小二乘法求出,根據(jù)試驗結(jié)果σs取為5MPa。限于篇幅,這里僅給出σz為2.5MPa與8.0MPa時不同時刻的蠕變參數(shù),見表2、3。3.3降低氧變速率及時間根據(jù)求得的蠕變參數(shù),繪制了其與時間的相關(guān)曲線,如圖6~12所示。體積模量K直接反映了巖石試件的體積隨時間的變化,K值的大小和變化與試件體積的收縮和膨脹息息相關(guān),不同應(yīng)力作用下的K值很難用統(tǒng)一的方程來描述,見圖6。因此,只能根據(jù)試驗結(jié)果求出。瞬時變形模量G0不隨時間發(fā)生變化,主要決定于瞬時應(yīng)變,而根據(jù)式(1),K值的大小會影響整個蠕變過程中試件的變形。不同應(yīng)力作用下黏彈性變形模量G1隨時間均呈衰減趨勢,見圖7。G1決定減速蠕變階段的變形量,此階段的蠕變曲線為衰減型,應(yīng)變逐漸增大到某一個固定值就基本保持恒定,即黏彈性變形εve在時間t→∞時的極值Sij/2G1,因此,G1隨時間不斷減小到某個值后基本保持穩(wěn)定。黏彈性黏滯系數(shù)η1與時間的關(guān)系見圖8。從蠕變方程分析,對于衰減型蠕變,η1決定蠕變達(dá)到穩(wěn)定階段時間的長短;對于非衰減型蠕變,η1影響蠕變進(jìn)入穩(wěn)定或加速階段的快慢。因此,η1的變化只影響蠕變達(dá)到穩(wěn)定階段的時間,不影響最終總變形量。從圖8可以看到,η1隨時間不斷增大到某一定值即保持穩(wěn)定,說明蠕變已達(dá)到穩(wěn)定階段。由蠕變方程式(1)可知,黏彈性應(yīng)變不僅與η1有關(guān),還應(yīng)同時考慮時間t的變化,因此,嘗試把t/η1作為一個整體變量考慮,其與時間的關(guān)系見圖9,可見二者基本呈直線關(guān)系,規(guī)律性比圖8更加明顯。黏塑性黏滯系數(shù)η2決定非衰減蠕變,即等速蠕變與加速蠕變,同時對減速蠕變也有影響。當(dāng)偏應(yīng)力超過長期強(qiáng)度σs時η2才會起作用,其產(chǎn)生的黏塑性應(yīng)變貫穿于蠕變的全過程。若偏應(yīng)力略大于σs或時間較短時,僅有等速蠕變產(chǎn)生;若偏應(yīng)力大于σs較多或時間足夠長時,等速蠕變則會轉(zhuǎn)變?yōu)榧铀偃渥?。黏塑性黏滯系?shù)η2與時間的關(guān)系見圖10。從圖中看到在減速及等速蠕變階段,η2會隨時間而增長,在加速蠕變階段,η2則會有明顯的降低。由蠕變方程式(1)可知,黏塑性應(yīng)變不僅與η2有關(guān),還應(yīng)同時考慮時間t的變化,因此,嘗試把t/η2作為一個整體變量,其與時間的關(guān)系見圖11,可見其規(guī)律性比圖10更加明顯。3.4擬合曲線分析根據(jù)上面得到的各蠕變參數(shù)與時間的關(guān)系圖,可以用非線性曲線擬合的方法得到二者之間相互關(guān)系的函數(shù)表達(dá)式。在進(jìn)行擬合時,本文是以ti時刻的參數(shù)與初始t0時刻的參數(shù)之間的相對關(guān)系來擬合。這樣做的目的在于,可根據(jù)初始t0時刻的參數(shù)來推求ti時刻的參數(shù),且從損傷力學(xué)的角度來看,建立初始t0時刻與后續(xù)ti時刻蠕變參數(shù)之間的關(guān)系也反映了材料特性隨時間的劣變過程。本文采用的初始t0時刻為t=5h。在尋找參數(shù)的變化規(guī)律時,往往需要根據(jù)所求參數(shù)與其他變量(如時間t)的組合來擬合,這樣得到的擬合曲線精度更高。黏彈性模量G1,以G1(t)/G1(t0)與時間的關(guān)系進(jìn)行擬合,數(shù)據(jù)見表4,數(shù)據(jù)點的分布見圖12。根據(jù)G1(t)/G1(t0)的數(shù)據(jù)分布,用方程y=a+bexp(-x/c)來擬合,結(jié)果如表5所示,擬合曲線見圖12。由擬合結(jié)果可以得到如下關(guān)系式:黏彈性黏滯系數(shù)η1,以[t/η1(t)]/[t0/η1(t0)]進(jìn)行擬合,數(shù)據(jù)見表6,數(shù)據(jù)點分布見圖13。根據(jù)[t/η1(t)]/[t0/η1(t0)]的數(shù)據(jù)分布,用方程y=a+bx來擬合,結(jié)果如表7所示,擬合曲線見圖13。由擬合結(jié)果可以得到如下關(guān)系式:黏彈性黏滯系數(shù)η2,以[t/η2(t)]-[t0/η1(t0)]進(jìn)行擬合,數(shù)據(jù)見表8,數(shù)據(jù)點分布見圖14。根據(jù)[t/η2(t)]-[t0/η2(t0)]的數(shù)據(jù)分布,用方程y=a+exp(x/b)來擬合,結(jié)果如表9所示,擬合曲線見圖14。由擬合結(jié)果可以得到如下關(guān)系式:3.5低應(yīng)力作用下的黏塑性應(yīng)變在三維蠕變方程中,考慮到滲透水壓p的作用,有效側(cè)壓力變?yōu)棣襵-p,則偏應(yīng)力張量Sij=(2σz-σx+p)/3,設(shè)σ=2σz-σx+p。根據(jù)求得的蠕變參數(shù),繪制了其與σ的相關(guān)曲線,如圖15~21所示。瞬時變形模量G0的大小與瞬時應(yīng)變ε0相關(guān),圖15是ε0與G0在不同應(yīng)力σ作用下的關(guān)系曲線?？梢?隨著σ提高,ε0不斷增大,G0則不斷降低。不同應(yīng)力作用下,各時刻黏彈性變形模量G1與σ的關(guān)系均呈先增大后衰減的趨勢,見圖16。在低應(yīng)力作用下試件沒有黏塑性變形,只有黏彈性變形,因此,產(chǎn)生的蠕變均為黏彈性。當(dāng)應(yīng)力提高超過長期強(qiáng)度σs后,既有黏彈性變形又有黏塑性變形,且黏塑性變形較大,則彈黏塑性階段的黏彈性變形量要小于黏彈性階段,因此黏彈性變形模量G1提高。若σ進(jìn)一步加大,則黏彈性變形量增加,G1又會相應(yīng)降低。從蠕變方程式(1)分析,σ/G1是黏彈性應(yīng)變εve的極限值,因此,把σ/G1作為一個整體變量考慮,其與應(yīng)力的關(guān)系見圖17,可見二者之間規(guī)律性比圖16更加明顯。黏彈性黏滯系數(shù)η1與應(yīng)力的關(guān)系見圖18,可以看到,隨著應(yīng)力提高,從只有黏彈性變形到黏彈、黏塑性變形共存,η1有先增大后降低的變化趨勢。這說明,當(dāng)應(yīng)力σ超過長期強(qiáng)度σs后,同一時刻的黏彈性變形達(dá)到穩(wěn)定階段的時間比σ小于長期強(qiáng)度σs時有所增加,但隨著σ進(jìn)一步加大,黏彈性變形會逐漸穩(wěn)定。由蠕變方程式(1)可知,在相同t時刻,黏彈性變形的極限值同時取決于G1/η1。因此,把G1/η1作為一個整體變量考慮,其與應(yīng)力的關(guān)系見圖19,可見二者之間的規(guī)律性比圖18更加明顯。黏塑性黏滯系數(shù)η2與應(yīng)力的關(guān)系見圖20,從圖中看到,η2會隨應(yīng)力的增大而降低,這也反映出黏塑性應(yīng)變在不斷加大。在加速蠕變階段(t=48h),η2會有更加明顯的降低。由蠕變方程式(1)可知,當(dāng)應(yīng)力σ超過長期強(qiáng)度σs后,黏塑性應(yīng)變?yōu)?σ-σs)t/η2。因此,把(σ-σs)/η2作為一個整體變量考慮,其與過應(yīng)力的關(guān)系見圖21。3.6擬合結(jié)果分析根據(jù)上面得到各蠕變參數(shù)與應(yīng)力的關(guān)系圖,用曲線擬合的方法得到二者之間相互關(guān)系的函數(shù)表達(dá)式。下面求解t0=5h時各蠕變參數(shù)與應(yīng)力的變化規(guī)律。瞬時變形模量G0,數(shù)據(jù)見表10,數(shù)據(jù)點分布見圖22。根據(jù)G0的數(shù)據(jù)分布,用方程y=a+bx來擬合,則a=2.07617,b=-0.11792,相關(guān)系數(shù)為0.957,擬合曲線見圖22。由擬合結(jié)果可以得到如下關(guān)系式:黏彈性模量G1,以σ/G1(σ)與應(yīng)力的關(guān)系進(jìn)行擬合,數(shù)據(jù)見表11,數(shù)據(jù)點分布見圖23。根據(jù)σ/G1(σ)的數(shù)據(jù)分布,用方程y=a+bexp(x/c)來擬合,結(jié)果a=1.20223,b=0.12340,c=4.24029,相關(guān)系數(shù)0.99933,擬合曲線見圖23。由擬合結(jié)果可以得到如下關(guān)系式:黏彈性黏滯系數(shù)η1,以G1(σ)/η1(σ)進(jìn)行擬合,數(shù)據(jù)見表12,數(shù)據(jù)點分布見圖24。根據(jù)G1(σ)/η1(σ)的數(shù)據(jù)分布,用方程y=a+bexp(x/c)來擬合,結(jié)果a=0.11071,b=0.00554,c=4.79939,相關(guān)系數(shù)0.96826,擬合曲線見圖24。由擬合結(jié)果可以得到如下關(guān)系式:黏彈性黏滯系數(shù)η2,以(σ-σs)/η2進(jìn)行擬合,數(shù)據(jù)見表13,數(shù)據(jù)點分布見圖25。根據(jù)(σ-σs)/η2的數(shù)據(jù)分布,用方程y=a+bln(x+c)來擬合,結(jié)果a=0.08208,b=0.20552,c=-3.75652,相關(guān)系數(shù)0.99905,擬合曲線見圖25。由擬合結(jié)果可以得到如下關(guān)系式:3.7應(yīng)變參數(shù)西原模量k將t0時刻蠕變參數(shù)與應(yīng)力σ的關(guān)系式,代入相應(yīng)的參數(shù)與時間t的關(guān)系式,可求出蠕變參數(shù)在(σ,t)條件下的表達(dá)式。考慮到蠕變參數(shù)為變量,則變參數(shù)西原蠕變方程可表示為式中:K′、G0′、G1′、η1′、η2′表示可變化的參數(shù)。體積變形模量K′是與試件體積應(yīng)變緊密聯(lián)系的變量,可通過試驗數(shù)據(jù)得到。瞬時變形模量G0′只與應(yīng)力σ有關(guān):黏彈性變形模量G1′與應(yīng)力σ、時間t有關(guān):黏彈性黏滯系數(shù)η1′應(yīng)力σ、時間t有關(guān):黏塑性黏滯系數(shù)η2′與應(yīng)力σ、時間t有關(guān),只在σ>σs時出現(xiàn):上述公式中的a~f