土體剪脹性對地基承載力的影響

胡世敬，姚圣磊，余 苗

(1.貴州交通職業(yè)技術(shù)學(xué)院，貴州 貴陽 550008；2. 重慶交通大學(xué) 交通土建工程材料國家地方聯(lián)合工程實驗室，重慶 400074)

土體由于其顆粒骨料堆積而成的特點而具有兩個基本特性：①壓硬性，即土的強度和剛度隨壓應(yīng)力增大而增大；②剪脹性，即土體的體積應(yīng)變和剪應(yīng)力有關(guān)[1]。筆者通過理論和有限元分析，討論剪脹性對地基承載力的影響，以期為工程建設(shè)和設(shè)計提供一定的參考意義。

1—凝聚分量；2—剪脹分量；3—摩擦分量；4—實測強度

T. W. Lambe[2]將土的抗剪強度分為3個基本分量：凝聚分量、剪脹分量和摩擦分量，如圖1。隨著應(yīng)變的增大，剪脹充分發(fā)揮作用，剪脹分量達(dá)到峰值；到達(dá)某一應(yīng)變后，土體體積不再增加，剪脹分量也逐漸消失。

有關(guān)剪脹角的描述并不統(tǒng)一，剪脹角在滑移線場理論中主要有兩種定義[3]：①實際塑性應(yīng)變ε與剪切應(yīng)變γ之間夾角；②破壞面(應(yīng)力滑移線)與速度矢量方向的夾角[4]。筆者采用后者，如圖2。

陳祖煜[5]分析了具有剪脹性的地基，巖土材料的剪脹角對地基承載力的影響。筆者基于文獻(xiàn)[6-8]，采用ABAQUS有限元軟件，根據(jù)GB 50007—2011《建筑地基基礎(chǔ)設(shè)計規(guī)范》(以下簡稱《規(guī)范》)規(guī)定的地基極限承載力判定標(biāo)準(zhǔn)，給出了地基承載力的數(shù)值計算方法，分析了剪脹角對地基承載力的影響，確定地基承載力時應(yīng)如何考慮土體材料剪脹性的影響，使得工程設(shè)計更加經(jīng)濟合理。

圖2 應(yīng)力比和剪脹角Fig.2 Stress ratio and dilatancy angle

1 Mohr-Coulomb對應(yīng)的剪脹方程

經(jīng)典彈塑性本構(gòu)理論Mohr-Coulomb采用的流動法則為[9-10]：

(1)

對于Mohr-Coulomb相關(guān)聯(lián)的流動法則，塑性勢函數(shù)g與屈服函數(shù)f相同。

由式(1)可以得到：

(2)

(3)

從而可以得到剪脹方程[9]：

(4)

而Mohr-Coulomb理論對應(yīng)的破壞準(zhǔn)則(用p，q，θσ表示)為[11]：

(5)

將式(4)代入式(5)，可以得到：

(6)

式中：φ為內(nèi)摩擦角；ψ為剪脹角。

顯然，式(6)中d<0，只能描述體積膨脹性(剪脹)的。

綜上，Mohr-Coulomb破壞準(zhǔn)則能考慮巖土材料的剪脹性，若要考慮土體材料的剪脹性，需要用Mohr-Coulomb進行計算。筆者在Mohr-Coulomb模型基礎(chǔ)上，借助于有限元軟件ABAQUS對條形地基承載力進行數(shù)值模擬分析。

2 地基承載力有限元分析

2.1 土體模型

土體材料的屈服準(zhǔn)則采用Mohr-Coulomb等面積圓屈服準(zhǔn)則，根據(jù)《規(guī)范》規(guī)定的條形地基承載力的確定方法進行數(shù)值模擬。

取條基寬度B=1 m，即地基承載力作用的部分，并將網(wǎng)格加密，材料的黏聚力c=70 kPa，內(nèi)摩擦角φ=30°或20°，彈性模量E=200 MPa，泊松比υ=0.3。由于求解的是平面應(yīng)變問題，且具有對稱性，故取一半進行分析。幾何圖形和有限元網(wǎng)格如圖3。

圖3 幾何模型及有限元網(wǎng)格Fig.3 Geometric model and finite element grid

荷載施加過程采用位移增量法，以研究地基土破壞的漸進過程?？紤]土體重度和初始應(yīng)力對承載力和變形的影響，將模型底部邊界、左部邊界約束其水平和垂直位移，對稱邊界約束其水平位移。

2.2 計算結(jié)果

圖4為位移云圖(ψ=φ=30° 時)，可見，其與常見的實際地基土破壞模式非常相近，說明有限元分析能真實反映地基土的破壞形態(tài)。

圖4 有限元分析位移云圖Fig.4 Cloud picture of displacement

圖5是等效塑性應(yīng)變云圖(ψ=φ=30° 時)。分析可知，土體的塑性應(yīng)變主要發(fā)生在地基下一定深度范圍內(nèi)，以及基礎(chǔ)側(cè)面一定寬度范圍內(nèi)。

圖5 等效塑性應(yīng)變云圖Fig.5 Cloud picture of equivalent plastic strain

由圖6可知，地基的極限承載力隨剪脹角的增加而增大，曲線可分為3個階段：

1)剪脹角在0～10°(φ=20°)，或0～15°(φ=30°)時為初始增加階段，極限荷載有一定幅度的增加。

2)剪脹角在10～15°(φ=20°)，或15～25°(φ=30°)為繼續(xù)增長階段，極限承載力增幅較大。

3)剪脹角15～20°(φ=20°)，或25～30°(φ=30°)為趨穩(wěn)階段，盡管隨著剪脹角的增加，但極限承載力增加不大，基本上趨于平穩(wěn)。

圖6 不同內(nèi)摩擦角對應(yīng)的地基承載力與剪脹角變化關(guān)系Fig.6 Relation of bearing capacity of ground vs dilatancy angle ψ

對應(yīng)于特定內(nèi)摩擦角，當(dāng)?shù)鼗休d力隨著剪脹角的增加而基本不變時，剪脹角即為合理剪脹角。因此由圖6可知：當(dāng)用有限元分析極限地基承載力時，剪脹角的合理取值范圍：ψ=(3/4～5/6)φ。

當(dāng)采用不同剪脹角計算對應(yīng)于相應(yīng)地基承載力的極限荷載時，由于剪脹角的影響，隨著剪脹角的增加，導(dǎo)致極限地基承載力增加，并且位移矢量圖大幅增加(圖7)，這充分說明剪脹角對承載力的影響較大。

圖7 不同剪脹角對應(yīng)的位移矢量圖Fig.7 The vector of displacement corresponding todifferent dilatancy angle

傳統(tǒng)塑性力學(xué)理論在分析土體極限平衡時，假定剪脹角與內(nèi)摩擦角相等。采用有限元分析極限地基承載力時，剪脹角小于內(nèi)摩擦角，即ψ<φ，并且大量實驗和工程實踐表明[12-15]：有限元計算結(jié)果與巖土材料的基本特性較為相符。

3 結(jié) 論

研究得出以下結(jié)論：

1)摩爾-庫倫模型及劍橋模型均能考慮土體的剪脹性和剪縮性。

2)根據(jù)ABAQUS有限元軟件的摩爾-庫倫模型，采用位移增量加載法模擬地基土的失穩(wěn)破壞過程，并確定條形地基承載力極限荷載。

3)地基承載力隨剪脹角的增加而增加，在計算地基承載力時必須考慮土體剪脹性的影響。

4)當(dāng)采用有限元法計算地基極限承載力時，應(yīng)考慮剪脹角的合理取值，即：ψ= (3/4～5/6)φ，以期得到較為可靠的計算結(jié)果。

[1] 沈珠江.理論土力學(xué)[M].北京:中國水利水電出版社,2000.

Shen Zhujiang.Theory of Soil Mechanics[M].Beijing:China Water Conservancy and Hydropower Press,2000.

[2] Lambe T W.A mechanistic picture of shear strength in clay [C]//Proc of ASCE Research Conference on Shear Strength of Cohesive Soils.Colorado,USA:University of Colorado Boulder ColoPress,1960.

[3] 鄧楚鍵,鄭穎人,王凱,等.有關(guān)巖土材料剪脹的討論[J].巖土工程學(xué)報,2009,31(7):1110-1114.

Deng Chujian,Zheng Yingren,Wang Kai,et al.Some discussion on the dilatancy of geotechnical materials[J].Chinese Journal of Geotechnical Engineering,2009,31(7):1110-1114.

[4] David Muir-Wood.Soil Behavior and Critical State Soil Mechanics [M].New York:Cambridge University Press,1990.

[5] 陳祖煜.土力學(xué)經(jīng)典問題的極限分析上、下限解[J].巖土工程學(xué)報,2002,24(1):1-11.

Chen Zuyu.Limit analysis for the classic problems of soil mechanics[J].Chinese Journal of Geotechnical Engineering,2002,24(1):1-11.

[6] 張國祥.潛在滑移線理論及淺基承載力分析[J].長沙鐵道學(xué)院學(xué)報,1997,15(2):95-99.

Zhang Guoxiang.Potential sliding curve theory and analysis of ground bearing capacity[J].Journal of Changsha Railway University,1997,15(2):95-99.

[7] 張培文,陳祖煜.剪脹角對地基承載力的影響[J].巖土力學(xué),2006,27 (3):463-465.

Zhang Peiwen,Chen Zuyu.Significance of soil dilatancy in bearing capacity of ground[J].Rock and Soil Mechanics,2006,27(3):463-465.

[8] 趙金勇,蔣剛,楊磊,等.剪脹角對地基極限承載力影響的有限元極限分析[J].南京工業(yè)大學(xué)學(xué)報,2008,30(5):28-32.

Zhao Jinyong,Jiang Gang,Yang Lei,et al.Effect of shear dilatancy angle on ultimate bearing capacity of foundation by FEM limit analysis[J].Journal of Nanjing University of Technology,2008,30(5):28-32.

[9] 陳愈炯,俞培基,魏汝龍.論土的剪脹性[J].水利學(xué)報,1964(6):76-80.

Chen Yujiong,Yu Peiji,Wei Rulong.Soil dilatancy theory[J].Hydraulic Engineering,1964(6):76-80.

[10] 孔位學(xué),芮勇勤,董寶弟.巖土材料在非相關(guān)聯(lián)流動法則下剪脹角選取探討[J].巖土力學(xué),2009,30(11):3278-3282.

Kong Weixue,Rui Yongqin,Dong Baodi.Determination of dilatancy angle for geomaterials under non-associated flow rule[J].Rock and Soil Mechanics,2009,30 (11):3278-3282.

[11] 鄭穎人,沈珠江,龔曉南.巖土塑性力學(xué)原理[M].北京:中國建筑工業(yè)出版社,2002.

Zheng Yingren,Shen Zhujiang,Gong Xiaonan.The Principles of Geotechnical Plastic Mechanics[M].Beijing:China Architecture and Building Press,2002.

[12] 費康,張建偉.ABAQUS在巖土工程中的應(yīng)用[M].北京:中國水利水電出版社,2010.

Fei Kang,Zhang Jianwei.The Application of ABAQUS in Geotechnical Engineering[M].Beijing:China Water Conservancy and Hydropower Press,2010.

[13] 莊茁,廖劍暉.基于ABAQUS的有限元分析和應(yīng)用[M].北京:清華大學(xué)出版社,2009.

Zhuang Zhuo,Liao Jianhui.Based on ABAQUS Finite Element Analysis and Application[M].Beijing:Tsinghua University Press,2009.

[14] 夏志高.塑性力學(xué)[M].上海:同濟大學(xué)出版社,1991.

Xia Zhigao.Plastic Mechanics[M].Shanghai:Tongji University Press,1991.

[15] 殷宗澤,鄧肯J M.剪脹土與非剪脹土的應(yīng)力應(yīng)變關(guān)系[J].巖土工程學(xué)報,1984,6(4):24-40.

Yin Zongze,Duncan J M.A stress-strain relationship for dilative and non-dilative soils[J].Chinese Journal of Geotechnical Engineering,1984,6(4):24-40.