何忠明，劉雅欣，曾新發(fā)，劉登生，職孟林

?

非飽和粗粒土基質(zhì)吸力與含水率及級配關(guān)系試驗研究

何忠明1, 2，劉雅欣2，曾新發(fā)3，劉登生4，職孟林5

(1. 長沙理工大學(xué) 公路養(yǎng)護(hù)技術(shù)國家工程實驗室，湖南 長沙，410014；2. 長沙理工大學(xué) 交通運輸工程學(xué)院，湖南 長沙，410014；3. 湖南城市學(xué)院 土木工程學(xué)院，湖南 益陽，413000；4. 中交(天津)生態(tài)環(huán)保設(shè)計研究院有限公司，天津，300450；5．河南省交通規(guī)劃勘察設(shè)計院，河南 鄭州，450000)

為了研究非飽和粗粒土基質(zhì)吸力與含水率及級配之間的關(guān)系，通過控制曲率系數(shù)和不均勻系數(shù)配制5種不同級配的粗粒土土樣，采用自行設(shè)計的模型裝置進(jìn)行降雨入滲試驗，并根據(jù)試驗結(jié)果繪制基質(zhì)吸力與含水率關(guān)系擬合曲線，得到滿足試驗條件的函數(shù)關(guān)系式，從而得到土?水特征曲線方程。研究結(jié)果表明：不均勻系數(shù)越大，土樣初始含水率越小，達(dá)到飽和狀態(tài)時含水率越小，且土樣達(dá)到飽和所需的時間越長；隨著脫水過程的進(jìn)行，土樣不均勻系數(shù)越大，其含水率的下降速度越慢，殘余含水率也越大；初始基質(zhì)吸力隨著土樣不均勻系數(shù)增大而不斷提高，同時，土樣的殘余基質(zhì)吸力隨著不均勻系數(shù)增大而增大；利用Van Genuchten模型對實測土?水特征曲線進(jìn)行擬合可獲得擬合參數(shù)，從而得到可應(yīng)用于數(shù)值分析的土?水特征曲線方程。

粗粒土；級配；基質(zhì)吸力；含水率；土水特征曲線

土?水特征曲線主要研究土體基質(zhì)吸力與含水率之間的關(guān)系，是非飽和土力學(xué)中主要研究內(nèi)容[1?2]。隨著粗粒土在工程中的廣泛應(yīng)用，研究者對粗粒土土?水特征曲線的研究也逐步深入。對非飽和土而言，土體的級配對其基質(zhì)吸力與含水率的分布特征有重要影響[3?4]。土體內(nèi)部水分子的定向移動總是從能量高的地方向能量低的地方流動，當(dāng)土體孔隙內(nèi)的水所含能量比孔隙之外的水能量高時，孔隙內(nèi)的水會流向孔隙外部，反之，土體就會將孔隙外部的水吸入其內(nèi)部[5?6]。土?水特征曲線反映非飽和土中含水率與基質(zhì)吸力之間的關(guān)系[7?8]，它對研究非飽和土的相關(guān)性質(zhì)非常重要，因為土體的諸多參數(shù)如抗剪強(qiáng)度、滲透性能等均能在土水特征曲線中得到體現(xiàn)[9?14]。VANAPALLI 等[15]研究了非飽和土的級配與土?水特征曲線之間的關(guān)系。李志清等[16]通過室內(nèi)外降雨試驗，推算了土體在非飽和情況下的抗剪強(qiáng)度、滲透性、土體礦物組成以及應(yīng)力等多種因素對土?水特征曲線的影響。毛尚之[17]基于特定膨脹土在不同應(yīng)力條件下的降雨試驗，得到應(yīng)力對土?水特征曲線的影響規(guī)律。SILLERS 等[18]通過對不同土體的土?水特征曲線模型進(jìn)行統(tǒng)計分析，得出現(xiàn)有土?水特征曲線模型的均值、中值和均方差均存在較大偏差。目前，人們針對級配變化對基質(zhì)吸力與含水率關(guān)系的影響研究較少，為此，本文作者擬通過控制曲率系數(shù)c和不均勻系數(shù)u來配制5種不同級配的粗粒土土樣；然后，采用自行研制的模型裝置進(jìn)行室內(nèi)降雨試驗，對試驗結(jié)果進(jìn)行分析后得到粗粒土基質(zhì)吸力、含水率隨級配的變化規(guī)律，以便為工程應(yīng)用提供參考。

1 土?水特征曲線基本理論

1.1 典型土?水特征曲線

典型土?水特征曲線一般分為3個階段：邊界效應(yīng)段、過渡段和非飽和殘余段，見圖1。在邊界效應(yīng)段內(nèi)，基質(zhì)吸力變化對含水率無顯著影響；當(dāng)基質(zhì)吸力增大到某個臨界值時，土體的基質(zhì)吸力大于孔隙的抵抗力，土體孔隙中的水將被排出，土體含水率會相應(yīng)降低，這個臨界值稱為進(jìn)氣值[11]。隨后，土?水特征曲線進(jìn)入過渡段，此時，土體基質(zhì)吸力逐漸增大，造成土中較大的孔隙開始向外排水，較小的孔隙也開始排水。隨著基質(zhì)吸力繼續(xù)增大，曲線進(jìn)入非飽和殘余段，在此階段，無論土體的基質(zhì)吸力如何變化，其含水率并不會產(chǎn)生明顯變化，而是趨向于某一個穩(wěn)定值，土體中除了細(xì)小的孔隙外，其他孔隙均不含自由水，此時，土體的含水率稱為殘余含水率[12]。

圖1 典型土?水特征曲線

1.2 土?水特征曲線數(shù)學(xué)模型

理論上，直接建立土體基質(zhì)吸力與含水率的關(guān)系非常困難，目前，國內(nèi)外學(xué)者通過試驗結(jié)果并結(jié)合相應(yīng)的數(shù)學(xué)模型來獲得土?水特征曲線，其中較經(jīng)典的數(shù)學(xué)模型有 Mualem模型、Fredlund和Xing模型、Van Genuchten模型、Gardner模型和Brooks-corey模型等，其中最常用的Van Genuchten模型公式為[19?21]

式中：為土壤含水率；r為殘余含水率；θ為飽和含水；，和為經(jīng)驗系數(shù)；為壓力水頭。

2 不同級配粗粒土室內(nèi)降雨模型試驗

2.1 土樣級配的確定

固定曲率系數(shù)c不變，改變不均勻系數(shù)u，制備試驗所需的5種級配土樣，具體級配見表1，試驗土樣粒徑見表2，試驗土樣的不均勻系數(shù)和曲率系數(shù)見表3。根據(jù)不同不均勻系數(shù)u，將土樣依次編號為土樣1、土樣2、土樣3、土樣4和土樣5，級配曲線按照橫坐標(biāo)為對數(shù)坐標(biāo)，如圖2所示。

表1 不同試驗土樣顆粒級配

表2 不同試驗土樣粒徑

注：10，30和60分別為分布曲線上小于該粒徑的試樣質(zhì)量占總質(zhì)量的10%，30%和60%的粒徑。

表3 試驗土樣的不均勻系數(shù)Cu和曲率系數(shù)Cc

1—土樣1；2—土樣2；3—土樣3；4—土樣4；5—土樣5。

2.2 試驗設(shè)備的制作及試驗用土的制備

本文使用的降雨入滲試驗設(shè)備主要由有機(jī)玻璃筒、土工篩、電子秤、含水率測試儀、張力計等組成。

1) 本試驗采用定制的有機(jī)玻璃筒如圖3所示，筒高為800 mm，直徑為400 mm，壁厚為5 mm。筒壁鉆有4個直徑為25 mm的圓孔。為了方便控制試驗用土的壓實效果，在有機(jī)玻璃筒的頂部設(shè)有壓實控制線，該線距離有機(jī)玻璃筒的頂部50 mm。

圖3 有機(jī)玻璃筒

2) 本試驗采用的含水率測試儀(型號B?10A，電源電壓為3.3 V，輸出電流為4~20 mA)測量范圍為0~100%，測量精度為1%；張力計(北京奧特斯達(dá)科技有限公司)測量范圍為0~80 kPa，測量精度為1 kPa。

試驗時，首先將從高速公路路堤施工現(xiàn)場取回的粗粒土晾曬，隨后將土樣分別采用圓孔直徑為0.10，0.25，1.00，2.00，5.00，10.00和20.00 mm的土工篩進(jìn)行篩分，將篩分后的土體按照試驗所需的級配曲線進(jìn)行稱量配比并攪拌均勻，制備符合級配要求的粗粒土。根據(jù)現(xiàn)場的擊實試驗得到土體的最大干密度為1.5 g/cm3。由于玻璃筒能承受土樣的最大體積為 94 200 cm3，因此，計算出試驗所需土體的總質(zhì)量為141.3 kg。

2.3 降雨入滲試驗過程

1) 首先，將蒸餾水沿著張力計內(nèi)壁緩慢注入張力計中，直至水完全充滿張力計，然后塞上橡皮塞，擰緊玻璃帽；通過紙巾擦拭張力計的陶瓷頭將張力計表盤讀數(shù)歸零。

2) 將試驗所需的粗粒土進(jìn)行分層壓實，當(dāng)土體到達(dá)鉆孔位置時，將調(diào)試好的張力計通過鉆孔插入，將含水率傳感器探頭埋設(shè)在與張力計相同的高度，并利用玻璃膠將張力計與鉆孔處進(jìn)行密封，以防漏水。安裝完成后，靜至24 h，記錄3個張力計和3個含水率測試儀的讀數(shù)，分別求其平均值，將平均值分別作為土樣的初始基質(zhì)吸力和初始含水率。

3) 本次試驗設(shè)置的降雨強(qiáng)度為1.16×10?6m/s，該降雨強(qiáng)度為工程所在地近10年來單日最大降雨強(qiáng)度，對降雨裝置進(jìn)行調(diào)節(jié)，使降雨均勻覆蓋整個土樣表面。

4) 開始降雨之后，每隔1 h測1次數(shù)據(jù)，直至土樣達(dá)到飽和狀態(tài)，隨后，打開玻璃筒底部的排水閥排水，同樣每隔1 h讀1次數(shù)據(jù)，直至數(shù)據(jù)不再發(fā)生變化為止。為了盡量減少誤差，每次測量讀3次，并將3次試驗結(jié)果求平均值，作為測量的基質(zhì)吸力與含水率的最終值。

2.4 試驗結(jié)果

2.4.1 含水率變化規(guī)律

圖4所示為不同試驗土樣含水率隨降雨時間的變化曲線。

1—土樣1；2—土樣2；3—土樣3；4—土樣4；5—土樣5。

從圖4可以看出：土樣1至土樣5的初始含水率分別為7.6%，7.3%，6.9%，6.5%和6.3%，呈現(xiàn)下降趨勢；飽和含水率分別為37.6%，36.7%，35.3%，34.6%和33.4%，同樣呈現(xiàn)下降趨勢；殘余含水率分別為17.2%，18.7%，19.6%，20.2%和21.6%，呈現(xiàn)上升趨勢。經(jīng)過分析可知：不均勻系數(shù)越大，土樣初始含水率越小，降雨開始后，土樣的體積含水率不斷升高，在達(dá)到飽和之前，各個土樣含水率變化趨勢幾乎一致；同時，隨著不均勻系數(shù)增大，土樣達(dá)到飽和所需時間越長，且土樣飽和含水率越小。降雨停止時，打開模型底部的排水閥，土樣進(jìn)入脫水過程，在這個過程中，土樣不均勻系數(shù)越大，其含水率下降速度越慢，并最終保持不變，此時，含水率稱為土體殘余含水率，隨著土樣不均勻系數(shù)增大而增大。

產(chǎn)生上述現(xiàn)象的原因是：土樣不均勻系數(shù)越大，土樣顆粒的粒徑分布越均勻，土體級配越好而土顆粒之間的孔隙越小，導(dǎo)致土體飽和滲透系數(shù)越小，進(jìn)而造成降雨條件下雨水下滲速度緩慢，表現(xiàn)出在相同降雨強(qiáng)度下，土樣不均勻系數(shù)越大，土體達(dá)到飽和狀態(tài)所需要的時間越長。因為土樣的不均勻系數(shù)大，造成土顆粒之間的孔隙小，雨水可填充的孔隙就少，因此，在降雨條件下，土樣的飽和含水率隨著不均勻系數(shù)增大而不斷減小，造成土樣在脫水過程中含水率降低速度慢，土樣的孔隙小，增強(qiáng)了土樣的持水能力。

2.4.2 基質(zhì)吸力變化規(guī)律

圖5所示為不同試驗土樣基質(zhì)吸力隨降雨時間的變化規(guī)律。

1—土樣1；2—土樣2；3—土樣3；4—土樣4；5—土樣5。

分析圖5可以得出：降雨開始前土樣1至土樣5初始基質(zhì)吸力分別為48，52，57，61和65 kPa，初始基質(zhì)吸力隨著土樣不均勻系數(shù)的增大而不斷升高；降雨開始后，土樣的基質(zhì)吸力隨著降雨時間持續(xù)而不斷降低，且最終基質(zhì)吸力為0 kPa，此時，土樣達(dá)到飽和狀態(tài)；降雨停止并打開模型底部的排水閥時，土樣進(jìn)入脫水過程，在這個過程中，土樣的基質(zhì)吸力隨著脫水時間的持續(xù)而不斷升高，經(jīng)過一段時間后，土樣的基質(zhì)吸力保持穩(wěn)定，此時，土樣1至土樣5的殘余基質(zhì)吸力分別為22，24，26，28和29 kPa，證明土樣的殘余基質(zhì)吸力隨著不均勻系數(shù)的增大而增大。

土樣1至土樣5的初始基質(zhì)吸力和最終殘余基質(zhì)吸力均隨著土樣不均勻系數(shù)增大而增大，這是因為土樣不均勻系數(shù)越大，土體的級配越好，而土顆粒之間的孔隙越小，造成土樣對水的吸附能力越強(qiáng)。土樣體積含水率的變化規(guī)律與基質(zhì)吸力的變化規(guī)律具有一定的相對性，這是因為土樣基質(zhì)吸力隨著土體體積含水率的增大而減小，如當(dāng)土樣的體積含水率達(dá)到最大時，其基質(zhì)吸力基本為0 kPa；當(dāng)土樣的體積含水率降到最低值時，則基質(zhì)吸力達(dá)到最大值。

3 土?水特征曲線繪制及擬合

3.1 土?水特征曲線繪制

土水特征曲線是研究非飽和土性質(zhì)中非常重要的曲線，其含義為土的含水率與吸力之間的關(guān)系曲線，它反映了土的持水能力和特性。本文選用體積含水率和基質(zhì)吸力等參數(shù)對試驗土樣進(jìn)行研究，通過模型試驗繪制基質(zhì)吸力與含水率關(guān)系擬合曲線。不同試驗土樣的基質(zhì)吸力與含水率關(guān)系曲線的擬合參數(shù)見表4，不同試驗土樣基質(zhì)吸力與含水率的關(guān)系見圖6。

從表4可知：土樣的不均勻系數(shù)對土體的基質(zhì)吸力及含水率均有重要影響，對于含水率而言，不均勻系數(shù)越大，土體含水率下降速度越小。從圖6可知：土樣1至土樣5這的初始含水率分別為37.6%，36.7%，35.3%，34.6%和33.4%，它們隨著土樣不均勻系數(shù)增大而減少；土樣1至土樣5的殘余含水率分別為17.2%，18.7%，19.6%，20.2%和21.6%，它們隨著土樣不均勻系數(shù)的增大而增大?？傮w上，土樣不均勻系數(shù)越大，曲線初始含水率越小，且含水率隨著基質(zhì)吸力的增大而不斷減?。划?dāng)基質(zhì)吸力基本保持不變時，不均勻系數(shù)大的土樣，其殘余含水率反而比不均勻系數(shù)小的土樣大。造成上述現(xiàn)象的原因是：土樣不均勻系數(shù)增大導(dǎo)致孔隙比減小和土體自身滲透系數(shù)減小，且空氣進(jìn)出更加困難，進(jìn)而導(dǎo)致土樣的排水能力較弱；同時，孔隙比小的土樣在脫水過程中，由于其孔隙間的持水能力較強(qiáng)，造成不均勻系數(shù)大的土樣殘余含水率比不均勻系數(shù)小的土樣大。

3.2 土?水特征曲線擬合

根據(jù)試驗結(jié)果對土?水特征曲線的形態(tài)和變化規(guī)律進(jìn)行描述，利用數(shù)學(xué)模型得到土?水特征曲線的經(jīng)驗方程[16]。本文以基質(zhì)吸力為自變量，體積含水率為因變量，選用Van Genuchten模型對土?水特征曲線進(jìn)行擬合分析。由于土樣1和土樣5的級配差距最大，為了便于對比，選取土樣1和土樣5進(jìn)行擬合，擬合結(jié)果分別如圖7和圖8所示。

從表5可以看出：土樣1和土樣5的土?水特征曲線擬合相關(guān)系數(shù)R均在0.99以上，說明擬合效果良好。模型經(jīng)驗系數(shù)和先通過擬合得到，代入Van Genuchten模型公式即可進(jìn)一步得到土?水特征曲線經(jīng)驗方程，用于后續(xù)粗粒土路堤穩(wěn)定性計算或者降雨入滲量計算。應(yīng)用本文得到的土水特征曲線，可對粗粒土數(shù)值分析時所需相關(guān)參數(shù)進(jìn)行合理賦值，使數(shù)值分析結(jié)果更加符合工程實際。

表4 不同試驗土樣基質(zhì)吸力與含水率關(guān)系曲線擬合參數(shù)

1—土樣1；2—土樣2；3—土樣3；4—土樣4；5—土樣5。

圖7 土樣1的基質(zhì)吸力與含水率關(guān)系擬合曲線

圖8 土樣5的基質(zhì)吸力與含水率關(guān)系擬合曲線

表5 土樣1和土樣5基質(zhì)吸力與含水率關(guān)系曲線的擬合參數(shù)

Table 5 Moisture content and matrix suction curve fitting parameter of soil sample 1 and 5

4 結(jié)論

1) 不均勻系數(shù)越大，土樣初始含水率越??；隨著不均勻系數(shù)增大，土樣達(dá)到飽和所需的時間越長，且土樣飽和含水率越小。

2) 土樣在脫水過程中，不均勻系數(shù)越大，其含水率下降速度越小，且最終殘余含水率隨著土樣不均勻系數(shù)的增大而增大。

3) 初始基質(zhì)吸力隨著土樣不均勻系數(shù)增大而不斷升高；土樣在脫水過程中，其基質(zhì)吸力隨著脫水時間的持續(xù)而不斷升高，且最終土樣的殘余基質(zhì)吸力隨著不均勻系數(shù)增大而增大。

4) 利用Van Genuchten模型對實測土?水特征曲線進(jìn)行擬合，獲得擬合參數(shù)，從而可得到應(yīng)用于數(shù)值分析過程中的土?水特征曲線經(jīng)驗方程。在后續(xù)粗粒土路堤穩(wěn)定性計算或者降雨入滲計算過程中，應(yīng)用本文方法得到的土水特征曲線對數(shù)值分析所需相關(guān)參數(shù)合理地賦值，可使數(shù)值分析結(jié)果更加符合工程實際。

[1] 黃文熙. 土的工程性質(zhì)[M]. 北京: 水利電力出版社, 1983: 68?91. HUANG Wenxi. Engineering properties of soil[M]. Beijing: China Water Power Press, 1983: 68?91.

[2] 馬吉倩, 付宏淵, 王桂堯, 等. 降雨條件下成層土質(zhì)邊坡的滲流特征[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2018, 49(2): 464?471. MA Jiqian, FU Hongyuan, WANG Guiyao, et al. Seepage characteristics of layered soil slope under rainfall conditions[J]. Journal of Central South University(Science and Technology), 2018, 49(2): 464?471.

[3] 王協(xié)群, 鄒維列, 駱以道, 等. 壓實度與級配對路基重塑黏土土?水特征曲線的影響[J]. 巖土力學(xué), 2011, 32(S1): 181?186.WANG Xiequn, ZOU Weilie, LUO Yidao, et al. Effects of compaction degree and gradation on soil-water characteristic curve of roadbed remolded clay[J]. Rock and Soil Mechanics, 2011, 32(S1): 181?186.

[4] 鄭俊杰, 郭震山, 崔嵐, 等. 考慮非飽和滲流與增濕膨脹下的膨脹土隧道穩(wěn)定性分析[J]. 巖土力學(xué), 2017, 38(11): 3271?3277. ZHENG Junjie, GUO Zhenshan, CUI Lan, et al. Stability analysis of expansive soil tunnel considering unsaturated seepage and moistening swelling deformation[J]. Rock and Soil Mechanics, 2017, 38(11): 3271?3277.

[5] 張鵬程, 湯連生, 姜力群, 等. 基質(zhì)吸力與含水量及干密度定量關(guān)系研究[J], 巖石力學(xué)與工程學(xué)報, 2013, 23(5): 2792?2797. ZHANG Pengcheng, TANG liansheng, JIANG Liqun, et al. Quantitative relationship between soil suction and water content and dry density[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 23(5): 2792?2797.

[6] 文海家, 張巖巖, 付紅梅, 等. 降雨型滑坡失穩(wěn)機(jī)理及穩(wěn)定性評價方法研究進(jìn)展[J]. 中國公路學(xué)報, 2018, 31(2): 15?29, 96. WEN Haijia, ZHANG Yanyan, FU Hongmei, et al. Research status of instability mechanism of rainfall-induced landslide and stability evaluation methods[J]. China Journal of Highway Transporation, 2018, 31(2): 15?29, 96.

[7] 簡文星, 許強(qiáng), 童龍云. 三峽庫區(qū)黃土滑坡降雨入滲模型研究[J]. 巖土力學(xué), 2013, 34(12): 3527?3533. JIAN Wenxing, XU Qiang, TONG Longyun. Rainfall infiltration model of Huangtupo landslide in Three Gorges Reservoir area[J]. Rock and Soil Mechanics, 2013, 34(12): 3527?3533.

[8] 張建, 李江騰, 林杭, 等. 降雨觸發(fā)淺層坡體失穩(wěn)的遲滯現(xiàn)象及其與土質(zhì)參數(shù)的關(guān)聯(lián)性[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2018, 49(1): 150?157. ZHANG Jian, LI Jiangteng, LIN Hang, et al. Delay phenomenon of shallow slope failure triggered by rainfall and its correlation with soil parameters[J]. Journal of Central South University(Science and Technology), 2018, 49(1): 150?157.

[9] 申存科, 遲世春, 賈宇峰. 考慮顆粒破碎影響的粗粒土本構(gòu)關(guān)系[J]. 巖土力學(xué), 2010, 31(7): 2111?2115. SHEN Cunke, CHI Shichun, JIA Yufeng. Constitutive relations of coarse grained soil with particle crushing[J]. Rock and Soil Mechanics, 2010, 31(7): 2111?2115.

[10] 劉磊, 殷坤龍, 徐勇, 等. 考慮降雨及庫水位變動的區(qū)域滑坡災(zāi)害穩(wěn)定性評價研究[J]. 巖石力學(xué)與工程學(xué)報, 2018, 37(2): 403?414. LIU Lei, YIN Kunlong, XU Yong, et al. Evaluation of regional landslide stability considering rainfall and variation of water level of reservoir[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(2): 403?414.

[11] 張磊, 張璐璐, 程演, 等. 考慮潛蝕影響的降雨入滲邊坡穩(wěn)定性分析[J]. 巖土工程學(xué)報, 2014, 36(9): 1680?1687. ZHANG Lei, ZHANG Lulu, CHENG Yan, et al. Study on stability of rainfall infiltration slope considering the influence of potential erosion[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(9): 1680?1687.

[12] 盧應(yīng)發(fā), 陳高峰, 羅先啟, 等. 土?水特征曲線及其相關(guān)性研究[J]. 巖土力學(xué), 2014, 29(9): 2481?2486. LU Yingfa, CHEN Gaofeng, LUO Xianqi, et al. Study on soil-water characteristic curve and its correlation[J]. Rock and Soil Mechanics, 2014, 29 (9): 2481?2486.

[13] 左威, 鄭明新, 孫斌. 基于PF土?水特征曲線的非飽和土滲透特性預(yù)測[J]. 中外公路, 2011, 38(2): 26?28. ZUO Wei, ZHENG Mingxin, SUN Bin. Performance prediction of unsaturated soil based on PF soil-water characteristic curve[J]. Journal of China & Foreign Highway, 2011, 38(2): 26?28.

[14] FREDLUND D G, XING A, HUANG S. Predicting the permeability function for unsaturated soils using the soil-water characteristic curve[J]. Canadian Geotechnical Journal, 1994, 31(4): 533?546.

[15] VANAPALLI S K, FREDLUND D G, PUFAHL D E. The relationship between the soil-water characteristic curve and the unsaturated shear strength of a compacted glacial till[J]. Geotechnical Testing Journal, 1996, 19(3): 259?268.

[16] 李志清, 胡瑞林, 王立朝, 等. 非飽和膨脹土SWCC研究[J]. 巖土力學(xué), 2016, 27(5): 730?737. LI Zhiqing, HU Ruilin, WANG Lichao, et al. SWCC study of unsaturated expansive soils[J]. Rock and Soil Mechanics, 2016, 27(5): 730?737.

[17] 毛尚之. 非飽和膨脹土的土?水特征曲線研究[J]. 工程地質(zhì)學(xué)報, 2012, 10(2): 129?133. MAO Shangzhi. Study on soil-water characteristic curves of unsaturated expansive soils[J]. Journal of Engineering Geology, 2012, 10(2): 129?133.

[18] SILLERS W S, FREDLUND D G, ZAKERZADEH N. Mathematical attributes of some soil-water characteristic curve models[J]. Geotechnical and Geological Engineering, 2011, 19(3/4): 243?283.

[19] GENUCHTEN V. A closed-from equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Science Society of America Journal, 1980, 18(14): 892?898.

[20] NIE Wen, HUANG Runqiu, CHEN Yulong, et al. Optimization of groundwater 1D multi-tank model based on physical experiments[J]. Electronic Journal of Geotechnical Engineering, 2014, 19(6): 8827?8841.

[21] 張昭, 劉奉銀, 張國平. 土在全含水率范圍內(nèi)持水及非飽和滲透特性的模型描述[J]. 巖土工程學(xué)報, 2014, 36(11): 2069?2077. ZHANG Zhao, LIU Fengyin, ZHANG Guoping. Models for water retention and unsaturated permeability in full range ofwater content[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(11): 2069?2077.

Experimental study on relationship among matric suction and moisture content and gradation of unsaturated coarse grained soil

HE Zhongming1, 2, LIU Yaxin2, ZENG Xinfa3, LIU Dengsheng4, ZHI Menglin5

(1. State Engineering Laboratory of Highway Maintenance Technology, Changsha University of Science & Technology, Changsha 410014, China;2. School of Traffic and Transportation Engineering, Changsha University of Science & Technology, Changsha 410014, China;3. College of Civil Engineering, Hunan City University, Yiyang 413000, China;4. China Communications(Tianjin) Environmental Protection Design and Research Institute Co. Ltd., Tianjin 300450, China;5. Communications Planning and Design Institute Co. Ltd of Henan Province, Zhengzhou 450000, China)

In order to study the relationship among matric suction and moisture content and gradation of unsaturated coarse grained soil, five kinds of coarse-grained soil samples with different gradations were prepared by controlling curvature coefficient and non-uniform coefficient. The rainfall infiltration test was conducted using a self-designed model device. The fit curves of moisture content with matrix suction were plotted based on the experimental results, and the functional relationship under the experimental conditions was obtained. The results show that the bigger the non-uniform coefficient, the smaller the initial moisture content of the soil sample, and the smaller the water content, the longer it takes for the soil sample to reach saturation. With the progress of the dehydration process, the larger the non-uniform coefficient of the soil sample, the slower the rate of decline of its moisture content, and the greater the final residual moisture content. The initial matrix suction increases with the increase of the non-uniform coefficient of the soil sample. The residual matrix suction of the soil sample increases with the increase of the non-uniform coefficient. Using the Van Genuchten model to fit the measured soil-water characteristic curve, the fitting parameters can be obtained, so that the empirical equation of the soil-water characteristic curve applied in the numerical analysis process can be obtained.

coarse grained soil; gradation; matric suction; water content; soil water characteristic curve

U416

A

1672?7207(2019)03?0712?07

10.11817/j.issn.1672-7207.2019.03.026

2018?06?08；

2018?08?18

國家自然科學(xué)基金資助項目(51678073，51508042，51838001) (Projects(51678073, 51508042, 51838001) supported by the National Natural Science Foundation of China)

曾新發(fā)，博士，副教授，從事巖土工程、結(jié)構(gòu)振動控制研究；E-mail：zengxinfa126@126.com

(編輯 陳燦華)