考慮土水吸附效應(yīng)的黏土收縮曲線方程

陳 盼，向 銳，魏小棋,3，韋昌富,3，王吉利

(1. 巖土力學(xué)與工程國家重點實驗室 中國科學(xué)院武漢巖土力學(xué)研究所，武漢 430071； 2. 中交城鄉(xiāng)建設(shè) 規(guī)劃設(shè)計研究院有限公司，武漢 430052； 3. 桂林理工大學(xué) 土木與建筑工程學(xué)院，廣西 桂林 541004)

在干燥脫濕條件下，隨著含水量的降低，土體具有體積收縮的特性.土體收縮曲線是定量描述土體收縮特性的數(shù)學(xué)方程，描述土體孔隙率與含水量之間本構(gòu)關(guān)系[1].土體收縮曲線是土體水力學(xué)行為描述中的一個十分重要的本構(gòu)關(guān)系，研究發(fā)現(xiàn)，土體收縮曲線在非飽和土力學(xué)中與土水特征關(guān)系曲線具有同等重要的地位.因此，對于土體收縮曲線的研究逐漸得到研究者們的重視.已有研究表明，土體的收縮變形特性對土體的力學(xué)與水力學(xué)特性均具有重要的影響[2-4].因此，土體收縮曲線的研究對土體體積收縮開裂、不均勻沉降、地基承載力與邊坡穩(wěn)定性的評估均具有重要的意義.

不少學(xué)者已經(jīng)針對土體的收縮行為開展了大量的研究工作，主要包括試驗測試[5-7]、機理分析[8-9]及數(shù)學(xué)模型等[10-11].一些學(xué)者還開展了鹽分對黏性土收縮特性的影響[12].當(dāng)前文獻中土體的收縮曲線按照收縮變形的發(fā)展特征通常分為4個階段：結(jié)構(gòu)收縮、比例收縮、殘余收縮及0收縮階段[10].不少研究者根據(jù)這一分段特性建立了土體收縮曲線的分段模型，通過確定收縮變形的特征點，每個區(qū)段分別采用一個函數(shù)進行描述[13].這類模型的缺點是采用的模型參數(shù)較多，不便于推廣應(yīng)用.因此，學(xué)者們進一步發(fā)展了用1個方程來描述土體的整個收縮過程.Fredlund等[14]提出了1個可用于描述重塑土體收縮的曲線方程.Peng等[15]提出了1個能夠用于結(jié)構(gòu)性土收縮特性描述的收縮曲線方程，該方程得到了廣泛的應(yīng)用[16-17].

上述方程均假定土體在收縮過程中一定的含水量范圍內(nèi)必定滿足0收縮階段或0收縮斜率狀態(tài).然而上述假定并不總是有效，試驗已經(jīng)證實低含水量條件下干燥的土體仍然能夠發(fā)生持續(xù)的變形.例如，持續(xù)干燥條件下黏土仍然會繼續(xù)收縮，低含水量條件下的持續(xù)收縮體變在總體變中占的比較較大，尤其對于膨脹土，低含水量條件下的體變不能忽略[16,18-19].因此，現(xiàn)有的土體收縮曲線方程在描述這類土體的體變收縮時存在缺陷，并且不能準(zhǔn)確預(yù)測全含水量變化下黏性土的體積變形.研究表明，土體收縮行為與土體的土水特性關(guān)系密切，然而文獻中較少考慮土體收縮與土水特性的關(guān)聯(lián)機制，建立的數(shù)學(xué)模型多屬于經(jīng)驗方程，缺乏一定的物理機制[16,20].此外，當(dāng)前文獻中采用分段預(yù)測收縮曲線的模型應(yīng)用較為廣泛，能夠僅利用1個方程統(tǒng)一描述結(jié)構(gòu)性與非結(jié)構(gòu)性土體的收縮曲線的研究甚少.

本文提出1個描述從飽和到完全干燥狀態(tài)下全含水量范圍的土體收縮曲線方程.該方程不僅能描述高含水率條件下結(jié)構(gòu)性與非結(jié)構(gòu)性黏土的收縮行為，而且能夠描述低含水量條件下黏性土的持續(xù)收縮行為.隨后，利用文獻中測得的土體收縮曲線數(shù)據(jù)，通過預(yù)測曲線與實測結(jié)果的對比驗證本文提出收縮曲線的有效性與可靠性.在此基礎(chǔ)上，與文獻中現(xiàn)有的方程進行對比，證實本文提出方程的優(yōu)勢所在.本文的研究成果將為干旱與半干旱地區(qū)水分蒸發(fā)引起的地表土體變形開裂及核廢料緩沖層在干燥過程中的變形老化等方面的研究及評估提供理論支持.

1 土體收縮曲線的基本概念

土體收縮曲線是指土體在干燥失水條件下土體積變化與含水量的關(guān)系，見圖1，圖中IL為液限指數(shù).土體從泥漿(或完全飽和)狀態(tài)開始失水收縮，其體積的收縮等于排出水的體積.隨著干燥過程的進行，由于氣體進入土體孔隙，土體開始達到非飽和狀態(tài)，土體的收縮體積大于排出水的體積.隨著含水量進一步降低，土體體積持續(xù)收縮，直至達到0含水量狀態(tài).按照土體形態(tài)與界限含水量可將干燥收縮過程中泥漿經(jīng)歷的性態(tài)的變化分為4個狀態(tài)：流態(tài)(液態(tài))、可塑態(tài)、半固態(tài)和固態(tài).從圖1中可以看出，流態(tài)與可塑態(tài)的界限含水量是土體的液限指標(biāo)，塑限是可塑態(tài)與半固態(tài)的分界，而縮限是半固態(tài)與固態(tài)的分界.值得一提的是，土力學(xué)中縮限的定義對黏土仍然值得探討.

對于結(jié)構(gòu)性黏土，從飽和狀態(tài)開始失水收縮，土體的收縮體積小于排出水的體積，導(dǎo)致收縮曲線偏離飽和線，并不會出現(xiàn)像泥漿等非結(jié)構(gòu)性黏土沿飽和線的收縮階段[10-11].此外，黏土收縮曲線的4階段劃分法中，通常最后1個階段均定為0收縮階段.然而試驗結(jié)果表明，并不是所有的土類均具備0收縮階段或0收縮斜率的變形特征，失水收縮過程與土體的礦物組成與含量密切相關(guān).對于黏性土(尤其是膨脹土)，即使在較低的含水量條件下，由于黏土礦物團聚體及層間的水分在較為干燥的大氣環(huán)境下，土中吸附水也能夠遷移出來，土體仍會發(fā)生持續(xù)的收縮[16, 18-19].這些研究結(jié)果表明被廣泛采用的4階段土體收縮曲線的劃分方法存在明顯的局限性.

文獻中土體收縮曲線的描述主要有以下3種：質(zhì)量含水率w-孔隙率e(圖2(a))、體積含水率θ-孔隙率e(圖2(b))，水土體積比δ-孔隙率e(圖2(c))[11, 14, 16, 20].其中e=Vv/Vs，水土體積比δ=Vw/Vs，Vv為孔隙的體積，Vs為固體顆粒的體積，Vw為孔隙水的體積.在δ-e描述框架下，土體處于完全飽和狀態(tài)時，δs=es，δs為飽和時水土體積比，es為飽和時孔隙率.土體開始干燥，仍然處于飽和狀態(tài)，將沿著對角線即飽和線δ=e收縮，直到氣體進入土體，出現(xiàn)非飽和狀態(tài)，收縮曲線偏離飽和線，見圖3(c).采用δ-e定義土體收縮曲線可以清晰地描述土體收縮過程中從飽和到非飽和狀態(tài)的轉(zhuǎn)變.因此，本文采用δ-e來表達土體收縮曲線.

2 考慮土水吸附效應(yīng)的土體收縮曲線方程

試驗表明，低含水量條件下土體的收縮變形及收縮斜率并不為0，土體在低含水量條件下的收縮由吸附水控制，而在高含水量條件下，土體的收縮行為主要受毛細水影響.因此，土體的收縮可以分為毛細水變化引起的收縮與吸附水變化引起的收縮兩部分[19-20].土體從飽和狀態(tài)開始脫濕，收縮主要由毛細水控制，脫濕到一定含水量狀態(tài)，毛細水排干，對土體收縮的貢獻可忽略不計.由于Fredlund等的模型不能很好地描述結(jié)構(gòu)性黏土的收縮行為，而Peng等[15]建立的S型的土體收縮曲線方程能很好地描述結(jié)構(gòu)性與非結(jié)構(gòu)性黏土由毛細水主導(dǎo)的土體收縮過程的孔隙率演化：

(1)

式中：er為毛細水排干土體收縮達到的孔隙率；α、β及γ為模型參數(shù).

考慮吸附含水量對毛細收縮過程中孔隙率變化的影響，基于式(1)，提出毛細水變化引起的土體孔隙率e1演化方程：

(2)

式中：x為考慮吸附含水量對毛細水控制的孔隙率變化的影響量.

(3)

因此，脫濕過程中土體收縮曲線的方程為

e=e1+Δe

(4)

(5)

當(dāng)δ=0時土體的孔隙率為

(6)

對式(5)求導(dǎo)可得土體的收縮率：

(7)

當(dāng)δ=0時，吸附收縮斜率為

(8)

將式(6)代入式(5)，得到本文建立的收縮曲線的簡單表達形式為

(9)

式(9)共有5個參數(shù)e0,α,β,γ及k.與分段收縮曲線方程相比[11, 13]，本文僅用1個方程即可描述考慮土水吸附效應(yīng)的全含水量變化范圍內(nèi)土體的收縮特性.與文獻中建立的單一收縮曲線方程相比[14, 15]，本文僅引入?yún)?shù)吸附收縮斜率k，即能夠描述由于土水吸附效應(yīng)控制的土體持續(xù)收縮及非0收縮斜率的行為.以下將對上述方程的準(zhǔn)確性與適用性進行驗證，并對土體的收縮特性進行分析與討論.

3 模型的驗證及結(jié)果分析

3.1 采用的試驗土樣及基本物性參數(shù)

本文用于驗證模型的試驗數(shù)據(jù)來自文獻[18，21-22].采用了黏土、膨脹土與高嶺土共7種土樣，其中倫敦黏土、倫敦黏土與粉土混合土及藍黏土3種土樣來自文獻[21]，3種膨脹土來自文獻[22]，Norwegian clay來自文獻[18].這7種土樣的基本物理特性參數(shù)列于表1.

表1 土樣物理特性參數(shù)Tab.1 Physical property index for used soils

3.2 模型驗證及結(jié)果討論

利用式(9)對實測結(jié)果進行預(yù)測，利用Excel中的Solver通過最小二乘法獲得式(9)的各個參數(shù).獲得的方程模型參數(shù)及相關(guān)系數(shù)見表2.不同黏土收縮體變的預(yù)測曲線與實測數(shù)據(jù)對比見圖3，其中飽和線表示e=δ曲線.從圖3(a)可以看出初始處于飽和狀態(tài)的倫敦黏土不具有結(jié)構(gòu)性，預(yù)測曲線與實測結(jié)果具有很高的相關(guān)性，表明式(9)能很好地預(yù)測倫敦黏土從飽和到完全干燥過程中土體孔隙率的變化.盡管低含水量條件下土體收縮變化較緩，但倫敦黏土在低含水量條件下仍繼續(xù)收縮.因此，擬合獲得的斜率k值較小.而對于倫敦黏土與粉土混合土樣，初始狀態(tài)稍有結(jié)構(gòu)性(起始狀態(tài)偏離飽和線).在干燥過程中土體持續(xù)收縮，特別是低含水量條件下土體收縮較倫敦黏土更為明顯，建立的方程能很好地定量描述這種具有初始結(jié)構(gòu)的混合土樣從飽和到干燥過程中土體的體積收縮.同時捕捉到了黏性土干燥過程中的非0收縮過程及非零收縮斜率，見圖3(b).同時預(yù)測的k也比倫敦黏土的收縮斜率大，表明k是1個反映低含水量條件下黏土收縮強烈程度的表征量.

隨后采用文獻[22]中的3種膨脹土對提出的收縮曲線方程進行了進一步的驗證，預(yù)測曲線與實測結(jié)果見圖4，方程的模型參數(shù)與相關(guān)系數(shù)見表2，表中R2為決定系數(shù).從圖4中可以看出3種膨脹土均具有一定的結(jié)構(gòu)性.初期脫濕干燥過程中，土體孔隙率的降低小于土體排出水體積的變化.從表2的相關(guān)性系數(shù)可以看出，式(9)很好地預(yù)測了3種具有初始結(jié)構(gòu)性的膨脹土在干燥條件下的完整收縮行為.從低含水量條件下的收縮體變來看，膨脹土2收縮程度較大，膨脹土1與3收縮程度較為接近，這種變化趨勢也體現(xiàn)在方程的k上.

3.3 與現(xiàn)有收縮曲線方程的比較

本節(jié)將式(9)與文獻中較為常見的兩個收縮曲線方程進行對比.Peng等提出的方程已在式(1)中給出，而Fredlund等提出的方程以δ為自變量，表示如下：

表2 式(9)的模型參數(shù)與模型預(yù)測的相關(guān)系數(shù)

(10)

式中：a、b是方程的模型參數(shù).

3個收縮曲線方程預(yù)測曲線與實測結(jié)果對比見圖5，模型參數(shù)及相關(guān)系數(shù)見表3.圖5(a)給出的藍黏土屬非結(jié)構(gòu)性黏土，含水比低于0.2條件下仍有明顯的收縮.從預(yù)測結(jié)果與實測數(shù)據(jù)的對比可以看出，本文提出的方程很好地預(yù)測了藍黏土從飽和到完全干燥過程中的體積收縮，特別是低含水量范圍內(nèi)的持續(xù)收縮與非0收縮斜率的行為.Fredlund等與Peng等提出的方程在一定含水量范圍內(nèi)均預(yù)測到0收縮階段，并不能較好地描述低含水量條件下黏土的持續(xù)收縮行為.圖5(b)所示為1種結(jié)構(gòu)性黏土的收縮曲線.從飽和到干燥條件下也觀察到了挪威黏土的持續(xù)收縮過程.預(yù)測結(jié)果表明，本文提出的方程預(yù)測曲線與實測數(shù)據(jù)相比，具有最高的相關(guān)性系數(shù)，主要原因在于式(9)更好地預(yù)測了結(jié)構(gòu)性挪威黏土在干燥過程中展示的非零收縮段及非0收縮斜率的行為.Fredlund等提出的方程的預(yù)測效果最差，相關(guān)性系數(shù)最低.主要原因在于Fredlund等提出的方程并不適合用來描述結(jié)構(gòu)性黏土的收縮行為.同樣Fredlund等與Peng等提出的方程也不能較好地描述低含水量條件下黏土的非0收縮段及非0收縮斜率的行為.上述結(jié)果表明本文建立的方程在定量描述結(jié)構(gòu)性與非結(jié)構(gòu)性黏土收縮行為上更為合理與適用.

表3 不同方程的模型參數(shù)與相關(guān)系數(shù)Tab.3 Model parameters and correlation coefficients for different equations

4 結(jié)論

通過總結(jié)試驗測定的土體收縮曲線數(shù)據(jù)，克服現(xiàn)有土體收縮曲線方程的局限性，本文提出了1個能夠描述從完全飽和到干燥條件下土體收縮行為的數(shù)學(xué)方程.通過與文獻中實測數(shù)據(jù)及現(xiàn)有模型的對比分析，驗證了本文提出方程的可靠性與有效性.主要結(jié)論如下：

(1) 在低含水量條件下，土水吸附效應(yīng)主導(dǎo)土體的收縮變形行為，使得黏性土仍能夠發(fā)生持續(xù)體變，具有非0收縮段及非0收縮斜率的典型特征.

(2) 本文建立的方程不僅能夠描述非結(jié)構(gòu)性黏土的收縮曲線，同時也能夠描述結(jié)構(gòu)性黏土的收縮變形特性.

(3) 由于考慮了土水吸附效應(yīng)對收縮變形的影響，本文建立的方程能夠很好地描述低含水量條件下黏土的非0收縮段及非0收縮斜率行為.