考慮屈服應力的非飽和土土水特征曲線研究

馮君，巫錫勇，馬悅，李曉寧,，朱寶龍，廖昕

(1.西南交通大學 地質工程系，四川 成都 610031；

2.西南科技大學 土木工程與建筑學院，四川 綿陽 621010)

?

考慮屈服應力的非飽和土土水特征曲線研究

馮君1，2，巫錫勇1，馬悅1，李曉寧1,2，朱寶龍2，廖昕1

(1.西南交通大學 地質工程系，四川 成都 610031；

2.西南科技大學 土木工程與建筑學院，四川 綿陽 621010)

摘要:利用非飽和土固結儀進行控制基質吸力條件下，凈豎向應力逐漸增大的固結試驗，得到不同基質吸力條件下土體的屈服應力；進行凈豎向應力為屈服應力時，不同基質吸力條件下的固結試驗，通過對試驗過程中排出水量的測試，計算出不同基質吸力條件下的穩(wěn)定含水率，進而得到屈服應力條件下土體的土水特征曲線；利用土體孔隙分布的分形模型對土水特征曲線進行擬合。研究結果表明：屈服應力是影響非飽和土土水特征曲線的重要因素，屈服應力下土水特征曲線明顯分為3個區(qū)段；同時，屈服應力下土水特征曲線具有良好的分形特性。

關鍵詞：非飽和土；屈服應力；土水特征曲線；分形特性

土水特征曲線[1-3](soil water content characteristic curve，或 SWCC)是非飽和土理論體系中重要的本構關系之一，它是聯(lián)系基質吸力與體積含水率、質量含水率或飽和度的重要曲線。Barbour[4]曾將土水特征曲線比作飽和土力學中的e-lgp壓縮曲線，可見土水特征曲線在非飽和土力學中占有相當重要的地位。在非飽和土理論體系中，SWCC的作用主要有：1)進行土體中基質吸力和含水率之間的換算；2)可以間接地反映出土體中孔隙大小的分布；3)能反映非飽和土的持水能力；4)分析非飽和土中水分運移規(guī)律。基于其重要性，國內外大量學者從試驗和理論方面對其進行了相應的研究。在試驗研究方面，建立土水特征曲線最重要的是基質吸力和含水率的測定，含水率的測量方法相對簡單，對于基質吸力的測量，已經(jīng)建立了多種測量方法，如Temper 儀法[5]、壓力板儀法[6]、滲析技術[7]、張力計法[8]、熱敏傳感器法[9]和軸平移技術法[10]等。在理論研究方面，主要分為經(jīng)驗模型和分形模型兩大類。常用的經(jīng)驗模型有Broods-Corey模型、Gardner模型、Van-Genuchten模型(簡稱VG模型)、Gardner-Russo模型、Fredlund-Xing模型和Campbell模型等[11]。但是這些模型的共同缺點是它們都是純經(jīng)驗性的，模型中的參數(shù)沒有明確的物理意義且難以測量。近年來，隨著分形理論的發(fā)展，出現(xiàn)了一些關于非飽和土的分形模型，利用分形方法所確定的非飽和土土水特征曲線中的參數(shù)與土體結構性質直接聯(lián)系起來，具有明確的物理意義，是較為理想的方法。目前關于非飽和土土水特征曲線的分形模型主要分為3類[12]：1)土體質量分布的分形模型；2)孔隙表面的分形模型；3)孔隙體積分布的分形模型。XU[13]結合相關試驗研究，理論分析了3種分形模型的優(yōu)缺點，認為孔隙體積分布的分形模型最為適合分析非飽和土土水特征曲線，并在此基礎上提出了利用分形理論建立的非飽和土土水特征曲線的通用表達式。陶高梁等[11]也結合孔隙體積分布的分形模型提出了一種用質量含水率表示的非飽和土土水特征曲線的分形模型。目前對于土水特征曲線的研究，主要側重于天然狀態(tài)下表層非飽和土體的研究，但是實際工程中的土體往往存在一定程度的埋深，故都會受到上覆荷載或者其他附加應力的作用。Charles等[14-15]都曾在土水特征曲線的研究中考慮固結壓力的影響；方祥位等[16]的研究表明應力路徑及凈平均應力對非飽和土的土水特征曲線有著重要的影響。因此進行應力作用下的土水特征曲線的研究具有重要的工程指導意義。本文在進行非飽和土土水特征曲線研究時，考慮凈豎向應力對其的影響，并分析凈豎向應力為土體屈服應力時，土水特征曲線的變化規(guī)律。研究以合肥肥東地區(qū)非飽和黏性土為研究對象，利用非飽和土固結儀對屈服應力作用下的非飽和土土水特征曲線進行試驗研究。

1土樣概況

試驗用土取自合肥肥東地區(qū)。其顆粒級配及初始物性指標見圖1和表1所示。從X衍射試驗結果得出，土中的礦物主要為蒙脫石(30.3%)和伊利石(18.9%)。在地質構造上，該地區(qū)發(fā)育有多道斷層，地下水多為基巖裂隙水并由于其下伏堅硬的砂巖與頁巖斷層的存在而極為匱乏，因此試驗點的的土體長期處于非飽和狀態(tài)。

表1　土樣的初始物理性質指標Table 1 Initial physical properties indicators of the soil sample

圖1　土樣顆粒級配曲線Fig.1 Grain-size distribution curve of the study soil

2試驗儀器與研究方案

試驗儀器為后勤工程學院非飽和土固結儀，參照陳正漢等[17]介紹。試樣制成高度和直徑分別為20 mm和61.8 mm的環(huán)刀樣。加載方案為先對土樣進行控制基質吸力，凈豎向應力增大的試驗，其中吸力值分別控制為50，100，150和200 kPa，凈豎向應力分別為：50，100，200，300，500和800 kPa，以確定研究土樣的屈服應力。然后將凈豎向應力控制為屈服應力時，進行不同基質吸力條件下的固結試驗，通過對試驗過程中排出水量的測試，計算出不同基質吸力條件下的穩(wěn)定體積含水率，進而得到屈服應力條件下土體的土水特征曲線。

3試驗結果

3.1　屈服應力

按照試驗方案進行室內固結試驗，得到土樣體變與凈豎向應力之間的半對數(shù)關系曲線εv-lgp，如圖2所示。

圖2　不同基質吸力條件下的體變~凈豎向應力關系曲線Fig.2 εv-lgp curve with different matric suction

圖2顯示，相同吸力條件下，εv-lgp大致分為2個直線段。0~200 kPa為第1直線段，200~800 kPa為第2直線段，明顯第2直線段斜率大于第1直線段，表明吸力一定時，在較大的凈豎向應力條件下，土樣體積變形增長較快。在同一凈豎向應力條件下，隨著吸力的增大，體積變形變小，表明吸力的存在對土體具有一定的硬化作用，這和文獻[18]的結論一致。參照Alonso等[19-20]確定屈服吸力的方法，以εv-lgp曲線拐點作為土樣豎向應力的屈服點，從圖2可以看出，針對所研究的土樣，其屈服應力值約為200 kPa。

3.2　土水特征曲線

3.2.1土水特性原始數(shù)據(jù)獲取

利用非飽和土固結儀對土樣進行控制凈豎向應力為200 kPa時的吸力增大的固結試驗。利用排水量和初始含水率可以計算出不同吸力條件下，土樣的穩(wěn)定含水率，進而得到屈服應力下非飽和土黏性土的土水特征曲線?？紤]到試驗中排水量量測系統(tǒng)的誤差和水分的蒸發(fā)等因素，由試樣初始和終止含水率的差值，可以算出實際排水量，再根據(jù)實際排水量校核量測值，校正結果如表3所示，土水特征曲線如圖3所示。

表3　試驗含水率校正

圖3　非飽和土土水特征曲線Fig.3 SWCC of the studied unsaturated soil

由圖3可知，屈服應力作用下的土體含水率隨著吸力的增大而減??；曲線以土水特征曲線中的2個特征點(s=50 kPa和s=200 kPa)為界， 將其分為3個區(qū)段。第1區(qū)段(0

3.2.2分形模型

1990年，Tyler和Wheatcraft用Sierpinski地毯模擬土體孔隙分布，導出了具有Sierpinski地毯結構的土體的水分特征曲線。根據(jù)Young-Laplace方程可以得到非飽和土的土水特征曲線[21]：

(1)

式中：θ和θs分別是非飽和土和飽和土的體積含水率；s和se分別是非飽和土的基質吸力和進氣值；D是土體體積的分維，如果是三維空間，式(1)中的D-2變?yōu)镈-3即可。針對所研究的非飽和土黏性土，結合質量含水率和體積含水率的關系，對其屈服應力作用下的土水特征曲線，利用Matlab軟件對試驗數(shù)據(jù)按照式(1)進行擬合，擬合結果如圖4所示。

從圖4可以看出，本文所研究的非飽和黏性土，其屈服應力作用下的土水特征曲線具有很好的分形特性，用土體體積孔隙分布的分形模型能較好地擬合。

圖4　非飽和土土水特征曲線Fig.4 SWCC of the studied unsaturated soil

4結論

1)屈服應力對非飽和土土水特征曲線的影響顯著。屈服應力作用下的土水特征曲線，當200 kPa

2)屈服應力作用下的土水特征曲線具有很好的分形特性，用土體體積孔隙分布的分形模型能較好地擬合。

參考文獻：

[1] 李志清,胡瑞林,王立朝,等.非飽和膨脹土SWCC研究[J].巖土力學,2006(5):49-53.

LI Zhiqing, HU Ruilin, WANG Lichao, et al.Study on SWCC of unsaturated expansive soil[J].Rock and Soil Mechanics,2006(5):49-53.

[2] 張欽喜,陳鵬,楊宇友,等.非飽和土土-水特征曲線試驗及在工程中的應用[J].北京工業(yè)大學學報,2012(8):71-75.

ZHANG Qinxi，CHEN Peng，YANG Yuyou，et al.SWCC test of unsaturated soil and its applications in engineering[J].Journal of Beijing University of Technology,2012(8):71-75.

[3] 李亮,劉興旺,趙煉恒,等.降雨入滲對路基穩(wěn)定性影響因素分析[J].鐵道科學與工程學報，2007,4(2)：19-23.

LI Liang, LIU Xingwang, ZHAO Lianheng, et al.Factors analysis of subgrade stability due to rainfall infiltration[J].Journal of Railway Science and Engineering, 2007,4(2)：19-23.

[4] Barbour S L.Nineteenth canadian geotechnical colloquium: the soil-water characteristic curve: A history cal perspective[J].Can Geotech J, 1998(35):873-894.

[5] 劉曉敏, 趙慧麗, 王連俊, 等.非飽和粉質黏土土-水特性的試驗[J].石家莊鐵道學院學報, 2001, 14(3): 78-82.

LIU Xiaomin ZHAO Huili, WANG Lianjun, et al.Research on soil-water character of unsaturated pulverescent clay by experiment[J].Journal of Shijiazhuang Railly Institute, 2001, 14(3): 78-82.

[6] 黨進謙, 李靖.非飽和黃土含水率與基質吸力的關系[J].水土保持通報, 1995, 15(4): 39-42.

DANG Jinqian, LI Jing.Study on correlation between water content and matric suction of unsaturated loess [J].Bulletin of Soil and Water Conservation, 1995, 15 (4):39-42.

[7] 毛尚之.非飽和膨脹土的土-水特征曲線研究[J].工程地質學報, 2002, 10(2): 129-133.

MAO Shangzhi.A study on soil-water characteristic curve of unsaturated soil[J].Journal of Engineering Geology, 2002, 10(2): 129-133.

[8] 王釗,楊金鑫,況娟娟,等.濾紙法在現(xiàn)場基質吸力量測中的應用[J].巖土工程學報, 2003, 25(4): 405-408.

WANG Zhao, YANG Jinxin, KUANG Juanjuan， et a1.Application of filter paper method in field measurement of matric suction[J].Chinese Journal of Geotechnical Engineering, 2003, 25(4): 405-408.

[9] 薛強,劉勇,劉建軍,等.非飽和土水特征曲線模型參數(shù)反演辨識研究[J].儀器儀表學報,2006,27(6):753-755.

XUE Qiang, LIU Yong, LIU Jianjun, et al.Study on model parameters identification of unsaturated soil water characteristic curve[J].Chinese Journal of Scientific Instrument,2006,27(6):753-755.

[10] 夏達忠,張行南,賈淑彬,等.基于土壤物理特性的土壤水分特征曲線推求方法[J].實驗室研究與探索,2010(10):31-33.

XIA Dazhong， ZHANG Xingnan， JIA Shubing， et al.Study on the determination of SWCC based on soil physical properties[J].Research an Exploration in Laboratory, 2010(10):31-33.

[11] 陶高梁,孔令偉,肖衡林,等.土-水特征曲線的分形特性及其分析擬合[J].巖土力學，2014，35(9)：2443-2447.

TAO Gaoliang, KONG Lingwei, XIAO Henglin，et al.Fractal characteristics and fitting analysis of soil-water characteristic curves[J].Rock and Soil Mechanics，2014，35(9)：2443-2447.

[12] 徐永福,董平.非飽和土的水分特征曲線的分形模型[J].巖土力學, 2002, 23(4): 400-405.

XU Yongfu, DONG Ping.Fractal models forthesoil-water characteristics of unsaturated soils[J].Rock and Soil Mechanics, 2002, 23(4): 400-405.

[13] XU Y F.Calculation of unsaturated hydraulic conductivity using a fractal model for the pore-size distribution[J].Computers and Geotechnics, 2004, 31(7): 549-557.

[14] Charles W W, Pang Y W.Influence of stress state on soil water characteristics and slope stability[J].Journal of Geotechincal and Geoenvironmental Engineering, 2000, 126(2): 157-166.

[15] Fernando A M, Marinho, et al.The Influence of the compaction energy on the SWCC of a residual soil[C]// Proceedings of sessions of Geo-Denver 2000.Colorado: American Society of Civil Engineers, 2000: 125-141.

[16] 方祥位, 陳正漢, 申春妮, 等.剪切對非飽和土土水特征曲線影響的探討[J].巖土力學, 2004, 25(9): 1451-1454.

FANG Xiangwei, CHEN Zhenghan, SHEN Chunni, et al.A study on effect of shear on soil-water charaeteristic curve of an unsaturated soil[J].Rock and Soil Mechanics, 2004, 25(9): 1451-1454.

[17] 陳正漢,扈勝霞,孫樹國,等．非飽和土固結儀和直剪儀的研制及應用[J].巖土工程學報, 2004, 26(2): 161-166.

CHEN Zhenghan, HU Shengxia, SUN Shuguo, et a1.Development and application of consolidation apparatus and direct shear apparatus for unsaturated[J].Chinese Journal of Geotechnical Engineering, 2004, 26(2): 161-166.

[18] 馮君，巫錫勇，朱寶龍,等.考慮土水特性的非飽和土地基沉降特性數(shù)值分析[J].鐵道科學與工程學報，2014,11(2)：54-59.

FENG Jun, WU Xiyong, ZHU Baolong.et al.Numerical analysis of the settlement characteristics for unsaturated soil foundation considering soil water characteristics[J].Journal of Railway Science and Engineering, 2014,11(2)：54-59.

[19] Alonso E E,Gens A, Josa A.A constitutive model for partially saturated soils[J].Geotechnique, 1990, 40(3): 405-430.

[20] 陳正漢.重塑非飽和黃土的變形、強度、屈服和水量變化特性[J].巖土工程學報，1999,21(1)：82-90.

CHEN Zhenghan.Deformation, strength, yield and moisture change of a remolded unsaturated loess[J].Chinese Journal of Geotechnical Engineering, 1999,21(1)：82-90.

[21] Tyler S W，Wheatcraft S W.Fractal process in soil water retention[J].Water Resources Research，1990,26(5)：1047-1054.

(編輯蔣學東)

Study of soil water characteristic curve of unsaturated soil with considering the yielding stress

FENG Jun1,2, WU Xiyong1, MA Yue1, LI Xiaoning1,3, ZHU Baolong3，LIAO Xin1

(1.Department of Geological Engineering, Southwest Jiaotong University, Chengdu 610031, China;

2.School of Civil Engineering and Architecture, Southwest University of Science and Technology,Mianyang 621010, China)

Abstract:The consolidation tests on the condition of controlling matric suction and increasing net vertical stress were carried out by using the unsaturated soil consolidation apparatus and the yielding stress of the studied soil was obtained.Then, the net vertical stress was controlled at the yielding stress, the consolidation tests with different matric suction were carried out.Through the measuring of drainage during the test process, the stable water content under different matric suction could be calculated.Moreover, the soil water characteristic curve under the yielding stress was obtained.Finally, the fractal model of soil pore presure distribution was used for fitting the soil water characteristic curve.The results show that, the yield stress has great influence on the soil water characteristic curve of unsaturated soil.The soil water characteristic curve under the yielding stress can be significantly divided into three sections.At the same time, it has good fractal characteristic.

Key words：unsaturated soil; yielding stress; soil water characteristic curve; fractal characteristic

通訊作者：馮君(1987-)，男，四川綿陽人，博士研究生，從事非飽和土、地基處理和復合地基等方面的研究；E-mail:sckid1987@163.com

基金項目：鐵道部科技研究開發(fā)計劃資助項目(2010G016-B)；西南交通大學2014年優(yōu)秀博士學位論文培育資助項目

收稿日期：2015-05-02

中圖分類號：TU443

文獻標志碼：A

文章編號：1672-7029(2015)06-1331-05