孫 凱， 李志清， 孔佑興， 周應(yīng)新， 王雙嬌

（1. 中國(guó)科學(xué)院 地質(zhì)與地球物理研究所 頁(yè)巖氣與地質(zhì)工程重點(diǎn)實(shí)驗(yàn)室，北京 100029； 2. 中國(guó)科學(xué)院大學(xué)，北京 100049；3. 中國(guó)科學(xué)院 地球科學(xué)研究院，北京 100029； 4. 云南省交通投資建設(shè)集團(tuán)有限公司，云南 昆明 650100）

0 引言

青藏高原地處高寒高海拔地區(qū)，凍土發(fā)育，地震活動(dòng)頻繁，我國(guó)川藏鐵路、中尼鐵路等工程建設(shè)面臨嚴(yán)峻的挑戰(zhàn)［1］。高原地區(qū)的凍土往往不是單一性質(zhì)的土體，而是含有塊石組成的復(fù)雜巖土體，這種巖土體也廣泛應(yīng)用于路基填筑、土石壩等人造工程中。在地震、車輛荷載等循環(huán)荷載作用下，寒區(qū)土石混合體構(gòu)筑物（如路基和邊坡）的抗震性和穩(wěn)定性逐漸成為研究的重點(diǎn)［2-5］。因此，研究動(dòng)荷載條件下凍結(jié)土石混合體力學(xué)特性具有重要的理論價(jià)值與工程意義。

相較于凍土的研究，凍土石混合體動(dòng)力學(xué)研究成果較少。土石混合體的力學(xué)特性受多種因素影響，包含有含石量、塊石形狀、基質(zhì)組分、礦物成分、粗糙度等因子［6-7］，其復(fù)雜的力學(xué)特性與凍土存在差異，崔凱等［8］在研究地震荷載作用下混合土動(dòng)力學(xué)特性時(shí)，采用Hardin-Drenvich 模型來(lái)擬合試驗(yàn)數(shù)據(jù)，計(jì)算了動(dòng)剪切模量和泊松比，分析得到隨著圍壓的增大動(dòng)剪切模量和泊松比增大。凌沖穎［9］開展了凍結(jié)粗粒土長(zhǎng)期循環(huán)作用下的力學(xué)特性研究，探究了凍結(jié)粗粒土疲勞損傷特性及強(qiáng)度特性。Zhang等［10］在單一溫度下開展了動(dòng)三軸試驗(yàn)，探究了不同粗粒含量下，凍結(jié)土石混合體的動(dòng)剪切模量和滯回環(huán)變化特征。但文中僅在一個(gè)溫度條件下開展實(shí)驗(yàn)，缺乏說(shuō)服力。Payan 等［11］通過(guò)將幾種粒度分布的砂與不同比例的二氧化硅非塑性粉土混合，研究了混合物的小應(yīng)變動(dòng)力特性，討論了非塑性細(xì)粉含量對(duì)動(dòng)彈性模量和阻尼比的影響?？梢?jiàn)，少有學(xué)者研究在不同溫度條件下，凍結(jié)土石混合體的動(dòng)彈性模量和動(dòng)應(yīng)力-應(yīng)變的變化規(guī)律。

厘清巖土體累積塑性應(yīng)變規(guī)律，探究其演化模型是解決循環(huán)荷載作用下路基動(dòng)力穩(wěn)定性和沉降問(wèn)題的關(guān)鍵。關(guān)于循環(huán)荷載對(duì)巖土體累積塑性應(yīng)變的演變模型，學(xué)者多采用經(jīng)驗(yàn)擬合的方法。目前已建立以Monismith 模型為主的經(jīng)驗(yàn)?zāi)Ｐ?，并且因其?jiǎn)潔實(shí)用，而被廣泛應(yīng)用。針對(duì)凍結(jié)巖土體，張斌龍等［12］通過(guò)動(dòng)三軸試驗(yàn)，研究?jī)鼋Y(jié)黏土的軸向累積塑性變化規(guī)律，建立了凍結(jié)黏土塑性安定和塑性蠕變臨界動(dòng)應(yīng)力表達(dá)式。Zhang 等［10］在開展低溫循環(huán)試驗(yàn)的基礎(chǔ)上，提出了修正后的Monismith 累積塑性應(yīng)變方程，能夠更好地?cái)M合試驗(yàn)結(jié)果，但文中未對(duì)模型參數(shù)進(jìn)行探究，試驗(yàn)結(jié)果也僅在-6 ℃條件下開展，合理性有待探討。累積塑性應(yīng)變模型受多種工況影響［13］，目前大部分研究已建立了相應(yīng)的模型，但仍存在試驗(yàn)條件少，試驗(yàn)條件與模型參數(shù)相關(guān)性差的問(wèn)題。

為了更進(jìn)一步了解寒區(qū)凍結(jié)土石混合體的動(dòng)力學(xué)特性，本文以循環(huán)荷載下凍結(jié)土石混合體的應(yīng)力-應(yīng)變機(jī)制為切入點(diǎn)，以累積塑性應(yīng)變和動(dòng)彈性模量為動(dòng)力學(xué)指標(biāo)，研究?jī)鼋Y(jié)土石混合體在不同條件下的動(dòng)應(yīng)力-應(yīng)變關(guān)系、滯回環(huán)特征、動(dòng)彈性模量變化規(guī)律，并給出了循環(huán)荷載作用下累積塑性應(yīng)變模型。

1 材料與方法

1.1 材料與樣品制備

試驗(yàn)采用重塑試樣，基質(zhì)土為粉土，最大干密度為1.5 g·cm-3，土體的塑限和液限分別為17.8%和32.1%。采用的塊石顆粒為石英顆粒，其顆粒尺寸介于6～9 mm 之間、天然密度為2.67 g·cm-3，點(diǎn)荷載強(qiáng)度為10.2 MPa。土顆粒及塊石粒徑級(jí)配曲線和樣品照片如圖1（a）所示。

圖1 土石粒徑Fig. 1 Soil and rock grain size： block photograph （a）， block grain size diagram （b）

試樣制備過(guò)程如下：取烘干后過(guò)2 mm 篩的土樣，加入一定量的蒸餾水，靜置24 h，配制成25%含水量的土樣。將土樣和塊石混合并攪拌均勻，放入直徑50 mm，高100 mm 的圓柱形模具中，分三層進(jìn)行壓實(shí)，配置成含石量為30%、40%、50%的土石混合體試樣。將未脫模的試樣連同模具包裹好保鮮膜后，放入-30 ℃條件的試驗(yàn)箱中，低溫保存24 h以上，而后在-5 ℃、-10 ℃、和-15 ℃條件下保存24 h，進(jìn)行試驗(yàn)。根據(jù)Medley［14］提出的土石閾值為0.05LC［LC為工程特征尺寸，此處取圓柱體試樣的直徑（50 mm）］，可得到本次試驗(yàn)的土石閾值為2.5 mm，故試驗(yàn)選用過(guò)2 mm 篩的粉土作為基質(zhì)土。通過(guò)試驗(yàn)研究表明塊石最大尺寸介于試樣直徑的1/6～1/4，可基本消除尺寸效應(yīng)的影響［15］，因而采用6～10 mm的塊石顆粒。

1.2 試驗(yàn)方法

本文采用KTL-DYN10動(dòng)三軸試驗(yàn)儀，進(jìn)行低溫單軸循環(huán)荷載試驗(yàn)，儀器的最大軸向荷載為10 kN，軸向力分辨率為5 N（軸向傳感器位于壓力室內(nèi)部），允許最大圍壓為2 MPa，軸向加載頻率范圍為0～5 Hz，加載波形包括正弦波和方形波。

為保證試樣在實(shí)驗(yàn)過(guò)程處于恒溫條件，在試驗(yàn)過(guò)程中對(duì)試樣進(jìn)行保溫處理，保溫方法簡(jiǎn)易圖如圖2（a）所示。具體過(guò)程如下，試樣從恒溫箱中取出后，先包裹保鮮膜（減少與空氣進(jìn)行水分交換）再用一層隔熱棉包裹（防止熱輻射和熱傳導(dǎo)），而后在外層包裹冰袋（避免試樣與外界接觸，防止與空氣進(jìn)行熱傳導(dǎo)），放入試樣底座上。試樣處于封閉的環(huán)境，基本無(wú)空氣流動(dòng)，從而降低因?qū)α鞫⑹崃?。該保溫措施參考Li 等［16］研究低溫土石混合體單軸壓縮條件下破壞模式的保溫方法，并對(duì)保溫措施進(jìn)行了改進(jìn)。同時(shí)通過(guò)試驗(yàn)測(cè)試樣表面溫度變化規(guī)律，發(fā)現(xiàn)在25 min 內(nèi)完成試驗(yàn)，樣品表面溫度變化約為1.4 ℃，該條件滿足實(shí)際實(shí)驗(yàn)要求。

圖2 試樣加載及保溫示意圖Fig. 2 diagram of specimen loading and maintaining temperature： thermal insulation （a）， loading （b）

本文在不同含石量（30%、40%、50%）、不同溫度（-5 ℃、-10 ℃、-15 ℃）、不同幅值（1.5 MPa、2.0 MPa、2.5 MPa）條件下開展了正弦振動(dòng)荷載單軸壓縮試驗(yàn)，加載曲線如圖2（b）所示。Rollins 等［17］采用頻率為0.01～2.00 Hz 進(jìn)行的礫石土循環(huán)加載測(cè)試試驗(yàn)，Aghaei 等［18］采用0.01～1.00 Hz 的堆石材料下使用穩(wěn)態(tài)正弦循環(huán)進(jìn)行的動(dòng)力學(xué)試驗(yàn)。國(guó)內(nèi)慣用的粗堆石料振動(dòng)力學(xué)試驗(yàn)加載頻率為1 Hz［10-19］，王慶志等［20］在研究?jī)鐾羺^(qū)青藏鐵路路基的動(dòng)力學(xué)特性時(shí)指出，青藏鐵路在多年凍土區(qū)的平均運(yùn)營(yíng)時(shí)速是120 km·h-1，對(duì)應(yīng)車載動(dòng)荷載頻率為1.3 Hz，綜上所述，本研究結(jié)合前人基礎(chǔ)上選取加載頻率為1 Hz。

按照張獻(xiàn)才［21］研究指出等幅循環(huán)振動(dòng)試驗(yàn)中幅值應(yīng)力取單軸下峰值應(yīng)力的0.5～0.9 倍，每組方案分別取三組不同的幅值，本實(shí)驗(yàn)將采用0.83σmax、0.67σmax、0.50σmax。為對(duì)比不同含石量和溫度下，土石混合體在相同幅值應(yīng)力條件下的動(dòng)力學(xué)特性，σmax取溫度為-5 ℃、含石量40%單軸狀態(tài)下的峰值應(yīng)力（通過(guò)試驗(yàn)得到σmax= 3 MPa）。徐學(xué)燕等［22］指出下限應(yīng)力在0.003～1.525 MPa 間改變時(shí)，對(duì)凍結(jié)試樣變形影響小。本試驗(yàn)選取0.3 MPa為循環(huán)的下限應(yīng)力，即軸向荷載為0.6 kN。綜上，根據(jù)實(shí)際試驗(yàn)條件和最大軸向應(yīng)力，選取的三組軸向幅值為1.5 MPa、2.0 MPa 和2.5 MPa，循環(huán)的下限應(yīng)力為0.3 MPa。具體工況條件見(jiàn)表1。

表1 試驗(yàn)工況Table 1 The test conditions

1.3 振動(dòng)荷載物理參數(shù)

動(dòng)彈性模量和累積塑性應(yīng)變?yōu)樵u(píng)價(jià)巖土材料動(dòng)力學(xué)性能的關(guān)鍵參數(shù)［23］。結(jié)合規(guī)范及實(shí)際試驗(yàn)數(shù)據(jù)，動(dòng)彈性模量是相鄰循環(huán)加載-卸載圍成封閉曲線的割線斜率（圖3），并按照式（1）計(jì)算得到。累積塑性應(yīng)變是指在加載—卸載的往復(fù)循環(huán)過(guò)程中不斷產(chǎn)生的塑性應(yīng)變。循環(huán)加載過(guò)程中由于塑性應(yīng)變持續(xù)累積，相應(yīng)的動(dòng)力學(xué)性能（例如動(dòng)彈性模量）也會(huì)產(chǎn)生變化，因此在達(dá)到不同應(yīng)變后，動(dòng)彈性模量可以認(rèn)為在相應(yīng)累積塑性應(yīng)變下試樣的動(dòng)彈性模量。

圖3 循環(huán)加載試驗(yàn)理論曲線Fig. 3 Theoretical curve of cyclic loading test

動(dòng)彈性模量按照下式計(jì)算：

式中：σd，max為滯回環(huán)上的最大應(yīng)力；σd，max為滯回環(huán)上的最小應(yīng)力；εd，max為滯回環(huán)上最大應(yīng)力對(duì)應(yīng)的應(yīng)變；εd，min為滯回環(huán)上的最小應(yīng)力對(duì)應(yīng)的應(yīng)變。本次試驗(yàn)的終止條件是軸向累積塑性應(yīng)變達(dá)到20%或循環(huán)次數(shù)達(dá)到3 000次。

2 結(jié)果與分析

2.1 動(dòng)應(yīng)力-應(yīng)變特征

土體的動(dòng)應(yīng)力-應(yīng)變關(guān)系是表征土體動(dòng)力學(xué)特性的關(guān)鍵，也是分析土體動(dòng)力破壞的重要基礎(chǔ)。圖4 為不同溫度條件下的動(dòng)應(yīng)力-應(yīng)變曲線試驗(yàn)結(jié)果，可見(jiàn)軸向應(yīng)力在每個(gè)周期內(nèi)波動(dòng)變化。在達(dá)到第一次循環(huán)的幅值前，應(yīng)力-應(yīng)變曲線有明顯的應(yīng)變強(qiáng)化階段。另外，圖4中應(yīng)力-應(yīng)變曲線組成的滯回環(huán)相對(duì)位置呈現(xiàn)出稀疏—密集—稀疏的規(guī)律，圖4（b）尤為明顯（為了清晰地展示試驗(yàn)結(jié)果，筆者在橫軸上設(shè)置了斷點(diǎn)隱去了部分應(yīng)力-應(yīng)變曲線）。以圖4（b）為例，初始8次循環(huán)加載，可逆塑性應(yīng)變逐漸增大（即累積塑性應(yīng)變逐漸增大），滯回環(huán)的相對(duì)位置表現(xiàn)出稀疏的特征，這主要是單個(gè)循環(huán)后產(chǎn)生的應(yīng)變逐漸增大；而后隨著循環(huán)次數(shù)的增大，應(yīng)變達(dá)到7%左右時(shí)，單個(gè)循環(huán)對(duì)應(yīng)的不可逆塑性形變量逐漸變小，即動(dòng)應(yīng)力-應(yīng)變曲線的滯回環(huán)相對(duì)位置表現(xiàn)出密集的規(guī)律；最后，當(dāng)應(yīng)變達(dá)到14%時(shí)，滯回環(huán)的相對(duì)位置又表現(xiàn)出稀疏的特征。在循環(huán)荷載作用下，前后滯回環(huán)的相對(duì)位置不重合，表明凍結(jié)土石混合體具有明顯的應(yīng)變累積性。

圖4 40%含石量條件下不同幅值應(yīng)力的動(dòng)應(yīng)力-應(yīng)變?cè)囼?yàn)結(jié)果Fig. 4 Dynamic stress-strain test results for different amplitude stresses under 40% rock content conditions：15 ℃ （a）， -10 ℃ （b）， -15 ℃ （c）

發(fā)生上述現(xiàn)象的可能原因如下：初期循環(huán)荷載作用下試樣內(nèi)部存在顆粒壓密和微裂紋閉合，導(dǎo)致單個(gè)循環(huán)荷載下對(duì)應(yīng)的塑性應(yīng)變大，滯回環(huán)相對(duì)位置表現(xiàn)出稀疏的特征；隨著微裂紋閉合，試樣抵抗循環(huán)荷載的能力增強(qiáng)，單個(gè)循環(huán)對(duì)應(yīng)的應(yīng)變減小，滯回環(huán)相對(duì)位置表現(xiàn)出密集的性質(zhì)；隨著應(yīng)變的繼續(xù)增大，試樣中裂紋開始發(fā)育和貫通，單個(gè)循環(huán)對(duì)應(yīng)的應(yīng)變?cè)龃?，滯回環(huán)相對(duì)位置表現(xiàn)出稀疏的性質(zhì)。

為了進(jìn)一步明確滯回環(huán)的變化特征，以-10 ℃、幅值應(yīng)力為1.5 MPa 下含石量40%的土石混合體試驗(yàn)結(jié)果為例，進(jìn)一步討論滯回環(huán)分階段的變化規(guī)律。圖5（a）～（d）給出了滯回環(huán)的演化規(guī)律，其中橫坐標(biāo)的長(zhǎng)度均為1%（應(yīng)變）?？偨Y(jié)發(fā)現(xiàn)以下結(jié)論：①單個(gè)滯回環(huán)初期階段循環(huán)的面積最大，但是隨著循環(huán)加載次數(shù)的增加，滯回環(huán)面積先逐漸變小后逐漸增大。此外，在初始階段，加載-卸載曲線重合度小，隨著應(yīng)變的增大，加載-卸載曲線逐漸重合［見(jiàn)圖5（b）中黑色和紅色滯回環(huán)逐漸重合］。而后，隨著當(dāng)循環(huán)次數(shù)達(dá)到183 次時(shí)（以40%含石量，溫度-10 ℃為例），相鄰滯回環(huán)又逐漸分離［圖5（c）和圖5（d）］。②動(dòng)應(yīng)力-應(yīng)變曲線呈現(xiàn)出非線性特征，滯回環(huán)初期為不閉合的橢圓形［如圖5（a）所示］，循環(huán)次數(shù)增大，滯回環(huán)逐漸變化為柳葉形［圖5（b）］，動(dòng)應(yīng)變繼續(xù)增大（循環(huán)加載次數(shù)到達(dá)1 676 時(shí)），最后滯回環(huán)呈現(xiàn)出近似閉合的扁長(zhǎng)橢圓形［圖5（c）和圖5（d）］。③圖5（a）～（d）中單個(gè)滯回環(huán)起點(diǎn)到終點(diǎn)的距離，有先減小后增大的趨勢(shì)，可見(jiàn)，隨著加載循環(huán)的增加，不可逆累積塑性應(yīng)變逐漸增加。

圖5 不同循環(huán)次數(shù)下相鄰循環(huán)處的應(yīng)力-應(yīng)變曲線Fig. 5 Stress-strain curves at adjacent cycles for different number of cycles： N=3， 4（a）；N=318， 319（b）； N=1 237， 1 238（c）； N=1 676， 1 677（d）

2.2 動(dòng)彈性模量的變化規(guī)律

動(dòng)彈性模量作為評(píng)價(jià)巖土體動(dòng)力學(xué)特性的關(guān)鍵指標(biāo)，動(dòng)彈性模量越大，試樣抵抗振動(dòng)荷載破壞的能力越強(qiáng)，反之抵抗循環(huán)荷載破壞的能力弱。動(dòng)彈性模量的影響因素包括材料本身特性和含水量等［23］。而在低溫土石混合體中，溫度是影響土石混合體特性關(guān)鍵因素［24］。本文通過(guò)循環(huán)單軸試驗(yàn)取得了試樣動(dòng)應(yīng)力-應(yīng)變關(guān)系，按照式（1）計(jì)算得到每次循環(huán)條件下的動(dòng)彈性模量。

圖6為不同溫度和含石量條件下試樣累積塑性應(yīng)變與動(dòng)彈性模量之間的關(guān)系，其中每個(gè)數(shù)據(jù)點(diǎn)代表每個(gè)循環(huán)次數(shù)所對(duì)應(yīng)的動(dòng)彈性模量，可見(jiàn)動(dòng)彈性模量存在明顯的波動(dòng)性。其波動(dòng)性可能是由于試樣中塊石咬合導(dǎo)致應(yīng)力集中所致，從而在宏觀上表現(xiàn)出彈性模量的波動(dòng)性。另外，初始變形階段彈性模量存在明顯的增大趨勢(shì)，這與前一節(jié)提到累積塑性應(yīng)變的稀疏階段類似，其原因可能是冰土組成的膠結(jié)物與塊石之間的裂縫逐漸閉合，并且試樣進(jìn)一步壓實(shí)，塊石和土體的咬合和嵌固作用更為明顯。微裂紋的閉合和塊石的嵌固作用增強(qiáng)，能夠在一定程度上提高試樣的動(dòng)彈性模量。

圖6 三種含石量條件下不同溫度動(dòng)彈性模量與累積塑性應(yīng)變的變化情況Fig. 6 Variation of dynamic elastic behavior modulus and cumulative plastic strain at different temperatures for three rock content conditions： 30% （a）， 40% （b）， 50% （c）

通過(guò)圖6可見(jiàn)，相同累積塑性應(yīng)變條件下，溫度越低，凍結(jié)土石混合體的動(dòng)彈性模量越大。另外溫度越低，數(shù)據(jù)點(diǎn)的密度越大，說(shuō)明達(dá)到對(duì)應(yīng)累積塑性應(yīng)變所需要的循環(huán)次數(shù)越多。對(duì)于溫度在-10 ℃（圖6黑色曲線）和-5 ℃（圖6藍(lán)色曲線）條件下的試樣，在5%應(yīng)變條件下，隨著應(yīng)變的增大，動(dòng)彈性模量有增大的趨勢(shì)。但相對(duì)于-15 ℃（圖6 紅色曲線）的試樣，動(dòng)彈性模量先增大到一定峰值而后逐漸減小，與-10 ℃和-5 ℃變化趨勢(shì)不同，這主要是達(dá)到相同應(yīng)變所需要的循環(huán)次數(shù)不同導(dǎo)致的。由于-15 ℃條件下，土石混合體達(dá)到相同應(yīng)變所需的循環(huán)次數(shù)增多，導(dǎo)致試樣內(nèi)部裂紋擴(kuò)展增多，從而降低試樣的動(dòng)彈性模量。

-15 ℃條件下動(dòng)彈性模量隨累積塑性應(yīng)變的增大先增大后減小，初步認(rèn)為可能是由于土石之間裂紋的開展導(dǎo)致。循環(huán)荷載條件下土石之間更易發(fā)生相對(duì)運(yùn)動(dòng)，循環(huán)次數(shù)越多，土石界面更易被破壞而產(chǎn)生微裂縫，土石之間的咬合和嵌固作用下降，這就導(dǎo)致其動(dòng)彈性模量降低。另外，在加載過(guò)程中隨著時(shí)間的推移試樣本身會(huì)產(chǎn)生熱量，內(nèi)部溫度會(huì)有增大趨勢(shì)［25］，導(dǎo)致冰晶的融化，冰的膠結(jié)作用降低，進(jìn)一步降低了土石之間的嵌固作用［26］；其次受設(shè)備保溫效果影響，達(dá)到相同的累積塑性應(yīng)變，-15 ℃所需要的循環(huán)次數(shù)多，加載時(shí)間長(zhǎng)，試樣表面溫度有所升高，導(dǎo)致試樣的彈性模量后期呈現(xiàn)下降的趨勢(shì)。

粗粒含量作為凍結(jié)土石混合體最為關(guān)鍵的參數(shù)之一，下文討論了含石量與凍結(jié)土石混合體彈性模量的規(guī)律。圖7是三種幅值應(yīng)力條件下不同含石量的動(dòng)彈性模量與累積塑性變形之間的關(guān)系，可見(jiàn)當(dāng)累積塑性應(yīng)變小于3%時(shí)，含石量30%試樣的動(dòng)彈性模量小于含石量40%和50%試樣的動(dòng)彈性模量。峰值應(yīng)力為1.5 MPa 和2.0 MPa 時(shí)，40%含石量的試樣大于50%含石量試樣的動(dòng)彈性模量；而當(dāng)峰值應(yīng)力為2.5 MPa 時(shí)，在初始循環(huán)階段含石量40%的試樣明顯小于含石量50%試樣的動(dòng)彈性模量，而隨著累計(jì)塑性變形的增大，含石量40%的試樣明顯大于含石量50%試樣的動(dòng)彈性模量?？梢?jiàn)隨著含石量的增大其動(dòng)彈性模量有增大趨勢(shì)，但在本次試驗(yàn)中極大值在40%與50%之間，受篇幅限制并未進(jìn)一步開展相應(yīng)的試驗(yàn)探究含石量50%含石量的增長(zhǎng)規(guī)律。

圖7 不同幅值應(yīng)力條件下-10 ℃動(dòng)彈性模量與累積塑性應(yīng)變之間的變化情況Fig. 7 Variation between dynamic modulus of elasticity and cumulative plastic strain at -10 °C for different amplitude stress conditions： 1.5 MPa （a）， 2.0 MPa （b）， 2.5 MPa （c）

隨著含石量的增大，動(dòng)彈性模量增大的趨勢(shì)原因可能是塊石的彈性模量明顯大于凍土的彈性模量，塊石含量越大，試樣的動(dòng)彈性模量越大。隨含石量的繼續(xù)增大，動(dòng)彈性模量可能會(huì)出現(xiàn)峰值點(diǎn)。但40%與50%含石量的動(dòng)彈性模量近似，可見(jiàn)彈性模量的增長(zhǎng)并不是線性增大的，這可能是受控于土石的咬合和嵌固作用，由于土與石之間的壓縮系數(shù)不同，這就導(dǎo)致在加載過(guò)程中，凍土的變形量大，而塊石基本不發(fā)生形變，塊石與凍土的膠結(jié)和嵌固作用更易被破壞，進(jìn)而導(dǎo)致其動(dòng)彈性模量降低。

2.3 累積塑性應(yīng)變隨循環(huán)次數(shù)變化的規(guī)律

為了更進(jìn)一步探討累積塑性應(yīng)變的規(guī)律，同時(shí)更直觀地反映低溫土石混合動(dòng)應(yīng)力-應(yīng)變滯回環(huán)從稀疏變密集再到稀疏的過(guò)程，下文探究了累積塑性應(yīng)變與循環(huán)次數(shù)之間的關(guān)系。通過(guò)圖8累積塑性應(yīng)變與循環(huán)次數(shù)之間的關(guān)系，發(fā)現(xiàn)低溫土石混合體的累積塑性應(yīng)變可分為三個(gè)階段：初始變形階段，穩(wěn)定變形階段和破壞階段［圖8（a）］。同時(shí)對(duì)比2.2節(jié)發(fā)現(xiàn)，其初始變形階段對(duì)應(yīng)了滯回環(huán)呈不閉合的橢圓形的階段，穩(wěn)定變形階段對(duì)應(yīng)了柳葉形的階段，破壞階段對(duì)應(yīng)了扁長(zhǎng)橢圓形階段（圖9）。在初始變形階段應(yīng)變隨著循環(huán)次數(shù)增大，累積塑性應(yīng)變快速增加。這是由于軸向荷載的作用，試樣內(nèi)部微裂縫閉合和試樣壓實(shí)過(guò)程中發(fā)生軸向應(yīng)變，其微裂紋是指由于凍脹而引起的裂紋［27］。而后進(jìn)入穩(wěn)定變形階段，隨著循環(huán)加載次數(shù)的增加累積塑性應(yīng)變?cè)鲩L(zhǎng)相對(duì)緩慢，且與循環(huán)次數(shù)線性關(guān)系；線性階段結(jié)束后，進(jìn)入破壞階段，累積塑性應(yīng)變快速增大。該階段試樣已發(fā)生明顯破壞。為了進(jìn)一步闡明累積塑性應(yīng)變的過(guò)程，下文結(jié)合峰值應(yīng)力和溫度，進(jìn)一步探討累積塑性應(yīng)變的規(guī)律。

圖8 不同含石量條件下的累積塑性應(yīng)變與循環(huán)次數(shù)之間的關(guān)系Fig. 8 The relationship between cumulative plastic strain and number of cycles at different stone content：30% （a）， 40% （b）， 50% （c）

圖9 累積塑性應(yīng)變與循環(huán)次數(shù)變化的三階段和滯回環(huán)特征圖Fig. 9 Cycle number change and hysteresis loop feature diagram （a）， and three stages of cumulative plastic deformation （b）， （c）， （d）

圖10 改進(jìn)的Monismith模型40%含石量循環(huán)次數(shù)與累積塑性應(yīng)變的擬合結(jié)果Fig. 10 Results of fitting the improved Monismith model of the number of cycles with 40% stone content to the cumulative plastic strain

相同條件下，幅值應(yīng)力越大，達(dá)到破壞應(yīng)變所需的循環(huán)次數(shù)越多。圖8中在相同循環(huán)加載次數(shù)條件下，幅值應(yīng)力2.5 MPa 條件下對(duì)應(yīng)的累積塑性應(yīng)變最大，其次是幅值應(yīng)力2.5 MPa，最后為幅值應(yīng)力1.5 MPa，即幅值應(yīng)力越大，相同循環(huán)次數(shù)下，累積塑性應(yīng)變?cè)酱蟆?/p>

累積塑性應(yīng)變隨循環(huán)次數(shù)變化的三階段如圖9（a）曲線所示：①在循環(huán)次數(shù)較少時(shí)處于初期變形階段，累積塑性應(yīng)變隨循環(huán)次數(shù)快速增大。②當(dāng)循環(huán)加載次數(shù)達(dá)到318時(shí)，進(jìn)入穩(wěn)定變形階段，累積塑性應(yīng)變的增長(zhǎng)放緩，呈現(xiàn)明顯的線性特征。為了更好地描述第二階段的變形規(guī)律，引入穩(wěn)定變形斜率k，穩(wěn)定變形階段作為試樣抵抗動(dòng)荷載的關(guān)鍵階段，該階段的斜率在一定程度上決定了試樣從穩(wěn)定變形階段向快速變形階段轉(zhuǎn)變所需要的循環(huán)次數(shù)，斜率越低所需的循環(huán)次數(shù)越多，試樣抵抗動(dòng)荷載的能力越強(qiáng)。表2 給出了不同條件下，低溫土石混合體在循環(huán)荷載作用下的穩(wěn)定變形階段斜率k。③當(dāng)循環(huán)次數(shù)大于1 676 時(shí)，進(jìn)入快速變形階段，其塑性變形隨循環(huán)次數(shù)增大而快速增大，曲線的斜率明顯增大。在幅值應(yīng)力為2.5 MPa 條件下，曲線快速越過(guò)穩(wěn)定變形階段進(jìn)入破壞階段，尤其是-5 ℃和-10 ℃條件最為明顯，在破壞階段裂紋擴(kuò)展速度快速增加，100 次循環(huán)加載下試樣應(yīng)變已超過(guò)15%，試樣已破壞。

表2 不同條件下穩(wěn)定變形階段的斜率和改進(jìn)Monismith指數(shù)模型參數(shù)Table 2 Slope of stable deformation stage under different conditions and parameters of improved Monismith exponential model

為了更好地描述凍結(jié)土石混合體累積塑性應(yīng)變的變化規(guī)律，前人提出了Monismith 模型［28］來(lái)預(yù)測(cè)累積塑性應(yīng)變的變化規(guī)律（式2），該模型因其簡(jiǎn)易而被廣泛運(yùn)用。結(jié)合本文試驗(yàn)，運(yùn)用Monismith模型對(duì)-10 ℃含石量40%試驗(yàn)結(jié)果進(jìn)行擬合結(jié)果如圖11（c）所示，可見(jiàn)其擬合度較差，無(wú)法準(zhǔn)確呈現(xiàn)累積塑性應(yīng)變隨循環(huán)次數(shù)增加的三個(gè)階段。

圖11 不同模型下累積塑性應(yīng)變與循環(huán)次數(shù)的擬合結(jié)果Fig. 11 Fitting results for cumulative plastic strain versus number of cycles for different models： results of the improved Monismith model fit （a）， model fitting results from References ［10］ （b）， Monismith exponential model fit results （c）

Monismith指數(shù)模型為

式中：ε為累積塑性應(yīng)變；N為循環(huán)次數(shù)；參數(shù)a和b與應(yīng)力水平和土體性質(zhì)有關(guān)。

Zhang 等［10］在研究低溫條件下碎石土的三軸動(dòng)力學(xué)特性時(shí)提出了修正的Monismith 模型［式（3）］，以期更好地呈現(xiàn)累積塑性應(yīng)變的特征。但在本文試驗(yàn)結(jié)果中擬合度較差，通過(guò)圖11（b）可知，該模型雖在穩(wěn)定變形階段和破壞階段有較好的擬合結(jié)果，但是無(wú)法呈現(xiàn)初期累積塑性變形的快速增長(zhǎng)趨勢(shì)。

式中：ε為累積塑性應(yīng)變；N為循環(huán)次數(shù)；參數(shù)a、b和應(yīng)力水平、土體性質(zhì)有關(guān)。

針對(duì)本文試驗(yàn)的累積塑性應(yīng)變特征，提出了改進(jìn)Monismith 指數(shù)模型［式（4）］，描述凍結(jié)土石混合體的累積塑性應(yīng)變的演化特征。

式中：ε為累積塑性應(yīng)變；N為循環(huán)次數(shù)；參數(shù)a、b和c與應(yīng)力水平、土體性質(zhì)和溫度相關(guān)。當(dāng)a= 0 時(shí)，模型為Monismith 模型，其中N為正整數(shù)，最小值為1。

采用試驗(yàn)結(jié)果來(lái)確定改進(jìn)Monismith 模型參數(shù)，圖10 給出了含石量40%條件下的擬合結(jié)果，可見(jiàn)其能夠較好擬合低溫條件下土石混合體的累積塑性應(yīng)變特征，并能體現(xiàn)累積塑性變形隨循環(huán)次數(shù)變化的三個(gè)階段，其擬合曲線R2均在0.98以上。

2.4 改進(jìn)Monismith 模型參數(shù)與試驗(yàn)條件的相關(guān)性

為了探究不同試驗(yàn)條件下，含石量、溫度、峰值偏應(yīng)力對(duì)改進(jìn)Monismith 指數(shù)模型參數(shù)的影響。本文結(jié)合實(shí)際擬合得到的結(jié)果，通過(guò)統(tǒng)計(jì)學(xué)的方法，利用SPSS 對(duì)表2 中的參數(shù)a、b和c進(jìn)行斯皮爾曼相關(guān)性分析和偏相關(guān)性分析，探究參數(shù)與含石量、溫度和峰值應(yīng)力的相關(guān)性，進(jìn)而明確改進(jìn)Monismith指數(shù)模型中參數(shù)受試驗(yàn)條件影響規(guī)律。斯皮爾曼相關(guān)系數(shù)以Charles Spearman 命名的相關(guān)系數(shù)，是評(píng)估兩個(gè)變量之間依賴性的關(guān)鍵指標(biāo)。斯皮爾曼相關(guān)系數(shù)能夠準(zhǔn)確迅速地確定兩個(gè)變量之間的顯著性關(guān)系，在金融、醫(yī)學(xué)、網(wǎng)絡(luò)等方面廣泛運(yùn)用［29］。且因其不需要變量服從正態(tài)分布，適用范圍廣。偏相關(guān)分析作為一種多變量相關(guān)性分析方法，它可以在控制一個(gè)變量的條件下，研究另外兩個(gè)變量之間的相關(guān)性，即在剔除其他因素影響的條件下，衡量?jī)蓚€(gè)變量之間相關(guān)性的方法［30-31］。以上兩種分析結(jié)果中斯皮爾曼相關(guān)系數(shù)rs和偏相關(guān)系數(shù)rz的絕對(duì)值越接近于1 表明兩個(gè)變量相關(guān)性越高（r＞0 正相關(guān)，r＜0代表負(fù)相關(guān)），P值用于衡量相關(guān)性的顯著水平，P＜0.05 表示兩個(gè)變量之間顯著相關(guān)，P＜0.01 表示2個(gè)變量之間極顯著相關(guān)。

利用斯皮爾曼相關(guān)性分析，對(duì)含石量、溫度和應(yīng)力水平進(jìn)行相關(guān)性分析，得到表三中的結(jié)果。含石量與參數(shù)b兩個(gè)變量之間的斯皮爾曼相關(guān)系數(shù)rs為0.501，P=0.00779（P＜0.01），可以得出兩變量之間的相關(guān)性是極顯著的。結(jié)合斯皮爾曼相關(guān)系數(shù)的定義可知含石量與參數(shù)b的相關(guān)性呈顯著正相關(guān)。說(shuō)明改進(jìn)的Monismith 模型中參數(shù)b受含石量變化的影響。含石量與參數(shù)c兩個(gè)變量之間的斯皮爾曼相關(guān)系數(shù)rs=0.4957，P=0.00855（P＜0.01），可以得出兩變量之間的相關(guān)性是極顯著的。結(jié)合斯皮爾曼相關(guān)系數(shù)的定義可知含石量與參數(shù)c呈極顯著正相關(guān)。說(shuō)明改進(jìn)Monismith 模型中參數(shù)c受含石量變化的影響。但由于含石量與參數(shù)c的rs小于含石量與參數(shù)b的rs，可見(jiàn)含石量與參數(shù)c的相關(guān)性小于與參數(shù)b的相關(guān)性，即含石量的改變對(duì)參數(shù)b影響更為明顯。

另外通過(guò)表3，可見(jiàn)溫度與參數(shù)a、b和c的P值均大于0.05，幅值應(yīng)力與參數(shù)a、b和c的P值同樣大于0.05，無(wú)法滿足合斯皮爾曼相關(guān)系數(shù)中P值的定義，可見(jiàn)溫度和幅值應(yīng)力與參數(shù)a、b和c之間的顯著性弱。考慮到斯皮爾曼相關(guān)性分析會(huì)受到其他變量的影響，不能準(zhǔn)確反映兩個(gè)變量之間的關(guān)系。偏相關(guān)分析通過(guò)控制變量的方法，量化兩個(gè)變量之間的相關(guān)性［30，32］。通過(guò)偏相關(guān)分析（表4）發(fā)現(xiàn)，在控制含石量條件下，溫度與參數(shù)a的偏相關(guān)系數(shù)為rs=-0.45217，P=0.04（P＜0.05），溫度與參數(shù)a呈現(xiàn)顯著負(fù)相關(guān)，可見(jiàn)參數(shù)a受溫度的影響較為明顯；在控制溫度的條件下，含石量與參數(shù)c的偏應(yīng)力系數(shù)rs=0.49086，P=0.01089（P＜0.05），含石量與參數(shù)c呈現(xiàn)顯著正相關(guān)；在控制軸向幅值的條件下，含石量與參數(shù)c的偏應(yīng)力系數(shù)rs=0.49833，P=0.00957（P＜0.01），含石量與參數(shù)c呈現(xiàn)顯著正相關(guān)。綜上，在控制含石量的情況下，Monismith 模型中的參數(shù)a與溫度呈現(xiàn)一定的相關(guān)性，可見(jiàn)參數(shù)a受溫度的影響較為明顯。斯皮爾曼相關(guān)性分析中已得到含石量與參數(shù)b和c呈顯著正相關(guān)性，而在偏相關(guān)分析中得到含石量與c仍存在顯著相關(guān)性。另外，在偏相關(guān)分析的條件下，軸向幅值與a、b、c的P值均大于0.05，由此可見(jiàn)軸向應(yīng)力與參數(shù)a、b和c的相關(guān)性弱，可見(jiàn)該模型在軸向力改變的條件下，對(duì)參數(shù)a、b和c的影響小，說(shuō)明該模型具有更廣的適用性。另外，參數(shù)c與試驗(yàn)條件的相關(guān)性弱，作為擬合測(cè)試的指數(shù)參數(shù)，隨試驗(yàn)條件變化，參數(shù)c的變化幅度小，難以探究其與試驗(yàn)條件的相關(guān)性。

表3 不同實(shí)驗(yàn)條件、不同模型參數(shù)下斯皮爾曼相關(guān)系數(shù)值Table 3 Spearman phase relationship values under different experimental conditions and different model parameters

表4 含石量、溫度和軸向幅值與a、b、c參數(shù)偏相關(guān)分析結(jié)果Table 4 Results of bias correlation analysis of stone content， temperature and axial amplitude with parameters a， b and c

綜上所述，含石量與參數(shù)b和參數(shù)c的相關(guān)性呈顯著正相關(guān)，相較而言含石量對(duì)參數(shù)b的影響更為明顯；在偏相關(guān)分析中，剔除含石量的影響，參數(shù)a受溫度的影響較為明顯。峰值應(yīng)力的改變對(duì)試樣模型參數(shù)a、b、c影響小，可見(jiàn)該模型在不同的峰值應(yīng)力條件下有更好的適用性。

3 結(jié)論

本文通過(guò)開展-5 ℃、-10 ℃和-15 ℃、不同含石量條件下單軸循環(huán)加載試驗(yàn)，討論了凍結(jié)土石混合體在循環(huán)荷載下的動(dòng)力學(xué)特性，比較了不同條件下累積塑性應(yīng)變和動(dòng)彈性模量的變化規(guī)律，得出如下結(jié)論：

（1）循環(huán)荷載作用下動(dòng)應(yīng)力-應(yīng)變曲線表現(xiàn)出疏-密-疏的三個(gè)階段。循環(huán)加載起始階段，滯回環(huán)的相對(duì)位置較為稀疏，這是由于初始循環(huán)階段中土石混合體顆粒間壓密和微裂紋閉合。隨著循環(huán)加載次數(shù)的增大，單個(gè)循環(huán)對(duì)應(yīng)的塑性變形量先逐漸變小，滯回環(huán)的相對(duì)位置表現(xiàn)出從稀疏變轉(zhuǎn)為密集的規(guī)律，其機(jī)理與微裂紋閉合后提升的抵抗動(dòng)荷載能力有關(guān)。而后滯回環(huán)的相對(duì)位置呈現(xiàn)出由密集轉(zhuǎn)變?yōu)橄∈璧囊?guī)律，這與土顆粒冰塊石組成的膠結(jié)體發(fā)生顆?；?，導(dǎo)致整體塑性變形增大有關(guān)。并在此基礎(chǔ)上總結(jié)了滯回環(huán)相對(duì)位置的演化規(guī)律：滯回環(huán)初期為不閉合的橢圓形，隨著循環(huán)加載次數(shù)的增大逐漸接近柳葉形，最后為近似閉合的扁長(zhǎng)橢圓形。

（2）比較了動(dòng)荷載作用下動(dòng)彈性模量的變化情況，隨著含石量的增大，彈性模量呈現(xiàn)增大的趨勢(shì)，而在當(dāng)含石量40%與50%時(shí)其增大趨勢(shì)不明晰。相較于含石量的變化對(duì)動(dòng)彈性模量的影響，溫度對(duì)動(dòng)彈性模量的影響更顯著，隨著溫度的降低，凍結(jié)土石混合體的動(dòng)彈性模量逐漸增大，且溫度越低，動(dòng)彈性模量的增大越明顯，其細(xì)觀機(jī)理與凍土和塊石之間的咬合和嵌固作用有關(guān)。另外，-15 ℃條件下動(dòng)彈性模量隨累積塑性應(yīng)變的增大先增大后減小，這可能與循環(huán)次數(shù)增大后，試樣溫度變化有關(guān)。

（3）比較了不同條件下，循環(huán)次數(shù)與累積塑性應(yīng)變之間的變化規(guī)律，幅值應(yīng)力越大，達(dá)到破壞所需的循環(huán)次數(shù)越少。并提出了累積塑性應(yīng)變與循環(huán)次數(shù)的改進(jìn)Monismith 指數(shù)模型，通過(guò)實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證了該模型的合理性，給出不同溫度、含石量和峰值應(yīng)力條件下的模型參數(shù)值。同時(shí)運(yùn)用統(tǒng)計(jì)學(xué)中斯皮爾曼相關(guān)系數(shù)和偏相關(guān)分析，明確了改進(jìn)Monismith 指數(shù)模型中的參數(shù)b值的受含石量的影響明顯，參數(shù)a受溫度的影響較為明顯，峰值應(yīng)力的改變對(duì)試樣模型參數(shù)a、b、c影響小。