凍融循環(huán)下含裂隙類巖石凍脹力的力學(xué)模型及實驗研究

王永巖，張金龍，朱思文

(青島科技大學(xué),山東 青島 266601)

寒區(qū)裂隙巖體受到凍融循環(huán)反復(fù)作用影響，在充水裂隙內(nèi)部反復(fù)凍脹、融縮，當(dāng)凍脹荷載超過巖石抗拉強(qiáng)度后裂隙端部發(fā)生宏觀斷裂，并在端部前緣出現(xiàn)局部區(qū)域微裂紋叢集，隨著凍融循環(huán)深入，宏觀斷裂不斷延展，最終導(dǎo)致整體失穩(wěn)。因此，凍融循環(huán)條件下裂隙巖體失穩(wěn)過程實質(zhì)上為反復(fù)凍脹荷載沿裂隙端部的疲勞斷裂過程。

關(guān)于寒區(qū)裂隙巖體損傷機(jī)制研究，近年逐漸引起國內(nèi)外學(xué)者關(guān)注。申艷軍等[1]采用相似材料配制 0°、30°、45°、60°和90°的單裂隙類砂巖試樣，并對裂隙內(nèi)部充分飽水，展開 10、20、30和50 次凍融循環(huán)實驗，而后觀測裂隙端部區(qū)域在不同裂隙傾角、凍融循環(huán)次數(shù)下局部化損傷效應(yīng)差異規(guī)律，對裂隙巖體在凍融循環(huán)作用下的局部化損傷效應(yīng)進(jìn)行分析，并結(jié)合斷裂力學(xué)應(yīng)力疊加理論，驗證了因局部化損傷效應(yīng)造成的裂隙端部斷裂特性及擴(kuò)展路徑規(guī)律。劉泉聲等[2]針對巖體在凍融循環(huán)下裂隙中會經(jīng)歷凍脹力的萌生、發(fā)展與消散，裂隙凍脹擴(kuò)展和巖體凍脹損傷程度受凍脹力控制，基于熱力學(xué)、滲流理論、界面力學(xué)和彈性理論建立了柱形封閉裂隙中凍脹力演化模型，對考慮水分遷移和不遷移兩種情況下的凍脹力量值進(jìn)行了研究。裴向軍等[3]對干燥及飽水狀態(tài)下裂隙巖石凍融的特征展開研究,研究發(fā)現(xiàn)裂隙巖石的脹縮率k都隨凍融次數(shù)的增加而增大，脹縮率的變化程度跟不同巖樣的孔隙率特征有關(guān)，飽水裂隙巖石在凍融過程中對溫度具有較強(qiáng)的敏感性，損傷表現(xiàn)在溫度損傷和凍脹損傷。路亞妮等[4]針對不同幾何特征的閉合裂隙類砂巖模型試樣，進(jìn)行凍融后的單軸壓縮實驗，分析裂隙巖體經(jīng)不同凍融循環(huán)次數(shù)后的物理力學(xué)特性，以及不同裂隙傾角、不同裂隙長度的巖樣對凍融巖體強(qiáng)度及破壞形態(tài)的影響，雖然對裂隙類巖石試件進(jìn)行了研究，但僅僅通過實驗現(xiàn)象總結(jié)出不同裂隙擴(kuò)展現(xiàn)象[5-11]，未通過理論公式得出凍融循環(huán)對含不同角度裂隙試件進(jìn)行深入分析、以及凍融循環(huán)強(qiáng)度、凍融循環(huán)次數(shù)之間的關(guān)系。本文的水泥砂漿類巖石材料為研究對象，結(jié)合彈性力學(xué)、斷裂力學(xué)和流體力學(xué)相關(guān)理論，通過實驗驗證了凍脹力模型的可行性。

1 實驗準(zhǔn)備

1.1 含單裂隙類巖石試件制作

圖1 試件模具和類巖石材料Fig.1 Specimen mold and similar rock material

圖2 實驗儀器Fig.2 Experimental apparatus

1.2 類巖石試件的加載儀器

圖2加載儀器為TAW-200電子式多功能材料力學(xué)試驗機(jī)(長春朝陽試驗儀器有限公司與青島科技大學(xué)聯(lián)合研制)。

1.3 類巖石試件凍融循環(huán)儀器

實驗的凍融循環(huán)是在恒溫條件下進(jìn)行的，為了提供恒溫環(huán)境，實驗中運(yùn)用RPH-80型恒溫恒濕試驗箱(東莞市環(huán)瑞環(huán)境測試設(shè)備有限公司)，凍融循環(huán)實驗的方法：把恒溫恒濕凍融循環(huán)試驗機(jī)設(shè)置為首先在溫度為-30 ℃條件下凍結(jié)6 h，然后在30 ℃的條件下融化6 h，如此反復(fù)，即每個凍融循環(huán)周期為12 h。該凍融循環(huán)試驗機(jī)的最低溫度控制在-30 ℃，濕度設(shè)置在100%。

2 類巖石試件凍結(jié)過程凍脹力

當(dāng)水的溫度低于其冰點時，水將發(fā)生相變，一般來講孔隙的初始飽和度超過91%，由于孔隙密封，結(jié)冰過程中水的體積變化就不能通過水的流動來補(bǔ)償，水-冰體系體積增加受到巖石骨架的約束，引起孔隙中未凍水的壓力及骨架中應(yīng)力的增加，導(dǎo)致孔隙體積的增大(彈性變形)，當(dāng)孔隙中應(yīng)力集中處的應(yīng)力達(dá)到巖石的抗壓強(qiáng)度時，裂隙開始擴(kuò)展，在其擴(kuò)展前，巖石骨架變形為彈性變形。由于凍脹作用下裂隙的擴(kuò)展本質(zhì)上是由冰的生長驅(qū)動的(雖然并不一定是由冰直接破裂，但是孔隙以及裂隙中水結(jié)冰過程中將會產(chǎn)生一定的凍脹力)，將水與巖石骨架作為冰的外部環(huán)境，計算中各參數(shù)的符號規(guī)定如下：

(1)孔隙的體積為Vp，孔隙中冰的體積為Vi，未凍水的體積為Vw，巖石骨架的體積為Vr，在計算中規(guī)定體積減小為正(或壓縮)；

(2)水的密度為ρw，冰的密度為ρi；

(3)水在無壓條件下結(jié)冰前后體積的膨脹系數(shù)為λa=9%，而在承壓條件下結(jié)冰前后體積的膨脹系數(shù)為λ，且λ<λa；

(4)水、冰及巖石骨架內(nèi)的應(yīng)力分別為Pw、Pi、Pr；

下面將分孔隙(裂隙)擴(kuò)展和不擴(kuò)展兩種情況，討論孔隙中水結(jié)冰過程中凍脹力的變化以及巖石骨架的響應(yīng)。

2.1 裂隙不擴(kuò)展條件下凍脹力計算模型

2.1.1 基本假設(shè)

(1)不考慮熱力耦合；

(2)孔隙的初始飽和度為100%，即在結(jié)冰過程中孔隙內(nèi)部為水-冰兩相系統(tǒng)；

(3)物體在P0的壓力下體積為V0。若壓強(qiáng)增加(P0→P0+dP)，則體積減小dV。則有K=dP/(-dV/V0)，K被稱為該物體的體積模量。水、冰均視作線彈性體，其體積模量分別為Kw(水的彈性模量)和Ki(冰的彈性模量)，且在整個過程中保持恒定；巖石對水和冰的約束簡化為彈性模型，忽略巖石骨架中的剪應(yīng)力和剪切變形，其等效體積模量為Kr(巖石骨架的彈性模量)且在整個過程中保持恒定；

(4)在整個過程中冰始終被水包圍，凍結(jié)從中心開始，整個過程中水、冰及巖石骨架的應(yīng)力均相等，即Pw=Pi=Pr。

2.1.2 凍脹力計算模型的建立

在水冰微小的相變過程中，水和冰的體積變化可表示為：

(1)

(2)

巖石、水和冰體系滿足以下條件：

(1)質(zhì)量守恒(冰和水之間相互轉(zhuǎn)化,質(zhì)量保持不變)

(3)

(4)

(2)彈性本構(gòu)關(guān)系

由于假設(shè)水、冰和巖石骨架均為彈性體，則三者增量形式的本構(gòu)關(guān)系可表示為：

(5)

(6)

(7)

其中，Kr由孔隙的形狀和巖石的彈性參數(shù)決定；dvr為孔隙體積的變化，在不考慮裂隙擴(kuò)展的條件下僅為巖石骨架的彈性變形。本質(zhì)上λ是由凍脹力(水相變?yōu)楸w積變大產(chǎn)生的壓力)決定的，首先體積為1的水在無約束條件下結(jié)冰，冰的體積為1+λa，此時冰壓力為零，之后冰的體積被壓縮至1+λ，此時冰的壓力為pi，冰被壓縮的體積為λa-λ。冰的壓力可表示為：

(8)

整理可得：

(9)

(3)變形相容條件

水在封閉孔隙中結(jié)冰引起體系壓力的增加，冰、水和巖石骨架中壓力的增加均來源于自身體積的壓縮。三者體積的壓縮均來源于相變膨脹，因此可得以下相容條件：

(10)

(4)能量轉(zhuǎn)換條件

相變膨脹所做的功=水內(nèi)力變化+冰內(nèi)力變化+巖石內(nèi)力變化。

聯(lián)合式(2)、(5)、(6)、(7)、(9)、(10)得(11)：

(11)

其中，Vr=Vr0×(P-P0)/Kr=Vr0×Pw/K；上式僅是dVw與dPw之間的關(guān)系式。

2.1.3 孔隙等效體積模量的計算方法

預(yù)制的矩形裂隙由于裂隙水結(jié)冰引起的裂紋擴(kuò)展視為平面應(yīng)變問題，裂隙受凍脹力影響，形狀發(fā)生變化，將變化后的預(yù)制的矩形裂隙簡化為橢圓柱狀孔隙結(jié)果。如此簡化也不無依據(jù)，可參考制作的類巖石試件，其預(yù)制的裂紋也并非真正的矩形，其周邊由于制作的誤差類似于橢圓，經(jīng)過應(yīng)力變化后更加近似為橢圓。裂隙的截面變化為橢圓，橢圓的長軸為a,短軸為b。在內(nèi)部壓力p的作用下孔隙發(fā)生變形，利用彈性力學(xué)理論，結(jié)合復(fù)變函數(shù)可求得p作用下長短軸的位移如下：

(12)

(13)

其中l(wèi)a,lb分別為長、短軸的位移，G=E/2(1+μ)為巖石的剪切模量；μ為巖石泊松比；m=(a-b)/(a+b)，R=(a+b)/2。裂隙的體積應(yīng)變?yōu)?/p>

(14)

其中，V是裂隙的體積。

將位移公式(12)、(13)帶入(14)可得(15)：

(15)

由于裂隙的變形為彈性，其變形應(yīng)滿足胡克定律，即：

(16)

對比式(15)、(16)可得出:

(17)

上式即為含矩形裂隙巖石等效體積模量的計算公式。

2.1.4 裂隙不擴(kuò)展條件下巖石彈性模量對凍脹力發(fā)展影響

由式(11)可以看出，孔隙水結(jié)冰過程中的凍脹力的大小是由孔隙的等效體積模量Kr決定的，而Kr又決定于巖石彈性模量和孔隙的形狀決定。其中計算的參數(shù)如表1所示。

表1 力學(xué)參數(shù)Table 1 Mechanical parameters

根據(jù)表1的力學(xué)參數(shù)，假設(shè)孔隙水結(jié)冰比例為1，根據(jù)公式(11)，Ki、Kw為已知量，可以利用MATLAB得到巖石體積模量與體系壓力之間的關(guān)系如圖3所示。

圖3 凍脹力和體積模量曲線Fig.3 Frost heaving force and bulk modulus curve

巖石體積模量決定著凍脹力的大小，體積模量越大相同凍結(jié)階段孔隙內(nèi)凍脹力越大。如圖3所示，巖石體積模量和孔隙水完全凍結(jié)時最大凍脹力呈非線性的正相關(guān)關(guān)系，且隨著模量的增加最大凍脹力趨向于一特定值，體積模量K和拉伸模量(或稱彈性模量)E、泊松比μ之間關(guān)系：E=3K(1-2μ)。假設(shè)泊松比為固定值，體積模量K和彈性模量E之間為正比例關(guān)系，所以凍脹力和彈性模量呈非線性的正相關(guān)關(guān)系。

3 含裂隙單軸壓縮實驗

3.1 常規(guī)單軸壓縮實驗

利用TAW-200電子式多功能材料力學(xué)實驗機(jī)對類巖石試件進(jìn)行單軸壓縮實驗，采取負(fù)荷加載方式，以50 N/s的加載速率進(jìn)行加載，實驗?zāi)康臏y量類巖石試件的彈性模量、泊松比、峰值強(qiáng)度，為類巖石的凍融實驗提供必要的信息。

因為本文采取的是類巖石相似材料做成的試件，而且類巖石試件內(nèi)的裂隙厚度0.8 mm比前人的0.3 mm更寬，所以由于材料和厚度的影響，這里的實驗結(jié)果與前人的研究結(jié)果不同。根據(jù)圖4～6可知，單軸抗壓強(qiáng)度隨角度的增加顯著增加，含裂隙試件的最低抗壓強(qiáng)度大約是完整試件的30%左右，含裂隙試件的最高抗壓強(qiáng)度大約是完整試件的67%左右，含貫通單裂隙對試件強(qiáng)度具有較大影響；彈性模量隨角度增加逐漸增加，泊松比隨角度的增加反而減小。因此，一方面可以驗證凍脹力和彈性模量之間的關(guān)系，另一方面可以研究不同角度裂隙試件凍融循環(huán)后試件強(qiáng)度性能變化及破壞方式。

圖4 不同裂隙傾角單軸抗壓強(qiáng)度Fig.4 Uniaxial compressive strength of specimens with different crack angle

圖5 不同裂隙傾角彈性模量Fig.5 Elastic modulus of specimens with different crack angle

圖6 不同裂隙傾角泊松比Fig.6 Poisson's ratio of specimens with different crack angle

3.2 凍融循環(huán)后單軸壓縮實驗

利用TAW-200電子式多功能材料力學(xué)實驗機(jī)對類巖石試件進(jìn)行單軸壓縮實驗，采取負(fù)荷加載方式，以50 N/s的加載速率進(jìn)行加載，針對不同循環(huán)次數(shù)(5次、10次、15次、20次)進(jìn)行單軸實驗，分析凍融循環(huán)后不同角度類巖石試件峰值強(qiáng)度，以及相同凍融循環(huán)次數(shù)對試件峰值強(qiáng)度的影響。

如圖7所示,含不同角度的裂隙試件的峰值強(qiáng)度隨著凍融循環(huán)次數(shù)的增加，不同角度裂隙試件的峰值強(qiáng)度逐漸減少，在20次凍融循環(huán)階段含不同角度的裂隙試件的峰值強(qiáng)度出現(xiàn)大幅度下降，不同角度的裂隙試件隨角度增加峰值強(qiáng)度仍然保持與常規(guī)單軸實驗相同趨勢，峰值強(qiáng)度隨角度增加而顯著增加。

由于不同角度裂隙試件的彈性模量不同，所以凍脹力隨彈性模量增加呈現(xiàn)增加的趨勢，也就是說，彈性模量大的可以承受更大的凍脹力，所以對于不同角度裂隙，0°裂隙更容易凍融循環(huán)損傷破壞，90°裂隙試件比0°裂隙試件更能承受凍脹力(圖8)，驗證了上文提出的凍脹力和彈性模量呈非線性的正相關(guān)關(guān)系。

圖7 峰值強(qiáng)度與凍融循環(huán)次數(shù)曲線Fig.7 Curves of peak strength and freeze-thaw cycle frequency

圖8 不同裂隙傾角強(qiáng)度比曲線Fig.8 Strength ratio curve

4 結(jié)論

(1) 巖石體積模量決定著凍脹力的大小，體積模量越大相同凍結(jié)階段孔隙內(nèi)凍脹力越大。如圖8所示，巖石體積模量和孔隙水完全凍結(jié)時最大凍脹力呈非線性的正相關(guān)關(guān)系，且隨著模量的增加最大凍脹力趨向于一特定值，體積模量K和拉伸模量(或稱彈性模量)E、泊松比μ之間有關(guān)系：E=3K(1-2μ)。假設(shè)泊松比為固定值，體積模量K和彈性模量E之間為正比例關(guān)系，所以凍脹力和彈性模量呈非線性的正相關(guān)關(guān)系。

(2) 由于不同角度裂隙試件的彈性模量不同，所以凍脹力隨彈性模量增加呈現(xiàn)增加的趨勢，也就是說，彈性模量大的可以承受更大的凍脹力，所以對于不同角度裂隙，0°裂隙更容易被凍融循環(huán)損傷破壞，90°裂隙試件比0°裂隙試件更能承受凍脹力。

參考文獻(xiàn)：

[1]申艷軍，楊更社，榮騰龍，等.凍融循環(huán)作用下單裂隙類砂巖局部化損傷效應(yīng)及端部斷裂特性分析[J].巖石力學(xué)與工程學(xué)報,2017,36(3): 562-570.

[2]劉泉聲,康永水,黃興,等.裂隙巖體凍融損傷關(guān)鍵問題及研究狀況[J].巖土力學(xué),2012,33(4): 971-978.

[3]裴向軍，蒙明輝，袁進(jìn)科，等.干燥及飽水狀態(tài)下裂隙巖石凍融特征研究[J].巖石力學(xué)，2017,38(7)：1999-2006.

[4]路亞妮，李新平，肖家雙.單裂隙巖體凍融力學(xué)特性試驗分析[J].地下空間與工程學(xué)報， 2014，10(6)：593-598.

[5]崔宏環(huán)， 劉建坤， 張立群，等.寒區(qū)路基改良土凍融循環(huán)與荷載耦合作用下?lián)p傷力學(xué)研究[J].冰川凍土，2016，38(8):1183-1188.

[6]耿珂.凍融循環(huán)對寒區(qū)隧道結(jié)構(gòu)凍脹力的影響[J].冰川凍土，2013,35(8):913-919.

[7]劉泉聲，黃詩冰，康永水，等.低溫凍結(jié)巖體單裂隙凍脹力與數(shù)值計算研究[J].巖土工程學(xué)報，2015,37(9)：1572-1580.

[8]路亞妮，李新平，吳興宏.三軸壓縮條件下凍融單裂隙巖樣裂縫貫通機(jī)制[J].巖土力學(xué)，2014,35(6)：1579-1584.

[9]蘇海健，靖洪文，趙洪輝.高溫后砂巖單軸壓縮加載速率效應(yīng)的試驗研究.2014,36(6):1064-1071.

[10]王永巖,朱羽萌,范夕燕，等.孔隙率對頁巖相似材料破壞影響的試驗研究[J].應(yīng)用力學(xué)學(xué)報，2016,33(4)：695-700.

[11]朱羽萌.孔隙率對頁巖相似材料破壞影響的試驗研究[D].青島：青島科技大學(xué)，2014.