胡衛(wèi)東　曹文貴　袁青松

摘要：在現(xiàn)有研究基礎(chǔ)上，結(jié)合臨坡條形基礎(chǔ)雙層粘土地基臨近邊坡的非對(duì)稱性特征和層狀地基特征，構(gòu)建出臨坡條形基礎(chǔ)雙層粘土地基的多滑塊組合單側(cè)滑移破壞模式；根據(jù)速度相容關(guān)系和速度三角形閉合條件，確定出多滑塊離散模式相對(duì)應(yīng)的機(jī)動(dòng)允許速度場(chǎng)；引入上限極限分析理論，導(dǎo)出臨坡條形基礎(chǔ)雙層粘土地基極限承載力計(jì)算模型，在此模型上采用序列二次規(guī)劃算法進(jìn)行優(yōu)化求解，建立了臨坡條形基礎(chǔ)雙層粘土地基極限承載力確定方法.利用Matlab的符號(hào)運(yùn)算功能和優(yōu)化函數(shù)編程計(jì)算進(jìn)行實(shí)例分析，并與有限元數(shù)值計(jì)算結(jié)果進(jìn)行對(duì)比分析，表明了本文研究方法的可行性與合理性.

關(guān)鍵詞：臨坡地基；雙層粘土地基；上限分析；條形基礎(chǔ)；極限承載能力

中圖分類號(hào)： U416.14 文獻(xiàn)標(biāo)識(shí)碼：A

Abstract：Combined with the characteristics of the asymmetry failure mode of the twolayer clay foundation near slope and layered properties， the unilateral sliding failure mode of twolayer foundation under strip footing adjacent to slope was built on the basis of the existing relevant researches. The failure mechanism was formed by many sliders. According to speed compatible relationships and triangle closed condition， the corresponding planar kinematically admissible velocity field was constructed. By introducing the upper limit analysis theory， the ultimate bearing capacity calculation model of strip footings over twolayer foundation soil was derived. A new approach for determining the ultimate bearing capacity of twolayer foundation under strip footing adjacent to slope was put forward by using sequential quadratic programming optimization algorithm. By the use of Matlab symbolic operation function and optimization function， the finite element numerical analysis results were analyzed and compared with programmed calculation. The feasibility and rationality of the proposed approach was shown.

Key words：ground foundation adjacent to slope； twolayer soils； upper limit analysis； strip footing； ultimate bearing capacity

目前，均質(zhì)淺基礎(chǔ)地基的承載力問題已經(jīng)有了較深入研究，但實(shí)際上工程大部分地基都是非均質(zhì)層狀地基，土層之間的力學(xué)指標(biāo)差別較大，使得成層土地基與均質(zhì)土地基在破壞機(jī)理和破壞模式上有著顯著不同，其相應(yīng)地基極限承載力也必然存在較大差別，因此，開展成層土地基承載力確定方法的研究有著重要的工程意義，這也正是本文研究的出發(fā)點(diǎn).

由于受到土層厚度和力學(xué)強(qiáng)度指標(biāo)等的影響，成層土地基的破壞模式更加復(fù)雜，給成層土地基承載力確定方法的研究造成了較大困難，盡管如此，不少學(xué)者通過長(zhǎng)期研究，在成層土地基承載力確定方法方面取得了長(zhǎng)足的進(jìn)步，形成了包括極限平衡法、滑移線法、極限分析法和數(shù)值分析法等多種研究方法[1-11].然而，已有成層土地基承載力確定方法的研究仍然存在一定程度的問題和不足，如Hasen[1]加權(quán)平均法取有效深度范圍內(nèi)不同土層厚度的加權(quán)平均強(qiáng)度后直接用于均質(zhì)土的Hasen公式，其方法由于簡(jiǎn)單實(shí)用得到了較為廣泛的運(yùn)用，但仍然只是一種半經(jīng)驗(yàn)方法；Meyerhof和Hanna[2]提出的剪切理論認(rèn)為上部土層發(fā)生沖剪破壞，下部土層發(fā)生整體剪切破壞，經(jīng)過大量模型試驗(yàn)檢驗(yàn)該理論運(yùn)用于成層土地基極限承載力計(jì)算具有一定合理性，但只能解決硬土層下臥軟土層地基承載力問題；楊永新[3]等基于滑移線場(chǎng)理論，研究了極限平衡區(qū)中上層粉土和下層砂土的地基承載能力，其研究方法的適應(yīng)性仍受到一定限制；Zhu[4]等及Burd和Frydman[5]等采用數(shù)值分析方法能有效解決復(fù)雜邊界問題和非均質(zhì)土層問題，但由于土體強(qiáng)度參數(shù)選擇的困難，也并沒有在實(shí)際工程應(yīng)用中被廣泛接受；Florkiewxicz[6]等基于極限分析上限理論，采用多滑塊組合破壞模式研究雙層和三層非水平分層地基的極限承載力，經(jīng)過與其它理論研究方法結(jié)果和實(shí)驗(yàn)結(jié)果的對(duì)比分析，驗(yàn)證了其正確性和有效性；Michalowski[7]等采用折線間斷面形式構(gòu)建了無粘性土下臥飽和軟粘土雙層地基的機(jī)動(dòng)速度場(chǎng)，對(duì)比分析了多種破壞模式的適用性，但是上述兩者確定的機(jī)動(dòng)速度場(chǎng)實(shí)際應(yīng)用到臨坡地基當(dāng)中仍有待進(jìn)一步改進(jìn).

近年來，臨近斜坡淺基礎(chǔ)下層狀粘土地基廣泛出現(xiàn)在各類巖土工程結(jié)構(gòu)中，但就目前來說，從臨坡地基出發(fā)開展臨坡層狀粘土地基破壞模式和承載力性狀的研究還不夠深入.已有的臨坡均質(zhì)粘土地基破壞模式大多假定滑動(dòng)面為圓弧或?qū)?shù)螺旋線[12]，但現(xiàn)實(shí)情況下，由于斜坡的存在和層狀構(gòu)造土的性質(zhì)各異，臨坡層狀粘土地基應(yīng)力應(yīng)變分布極為復(fù)雜，使得其在地基破壞模式上與臨坡均質(zhì)粘土地基相比已發(fā)生了很大改變，層狀地基越不均質(zhì)，真正的最危險(xiǎn)滑動(dòng)面偏離圓弧滑動(dòng)面越遠(yuǎn)，而且也并非對(duì)數(shù)螺旋線就能逼近真實(shí)滑動(dòng)面，因此，臨坡層狀粘土地基的任意滑動(dòng)面形狀的確定尤為重要.采用多滑塊組合破壞模式的上限極限分析方法，無需事先對(duì)極限狀態(tài)下的連續(xù)滑動(dòng)面形狀做出各種假設(shè)，能有效結(jié)合地基臨近斜坡和層狀構(gòu)造的特征，真實(shí)反映層狀地基粘土重度、粘聚力、土層厚度、埋深和坡角等不同影響因素作用下的破壞模式，而且，上限極限分析方法能避免對(duì)滑塊間作用力過多假設(shè)，不需考慮土體復(fù)雜的應(yīng)力應(yīng)變關(guān)系且理論嚴(yán)謹(jǐn)，這就為臨坡成層土地基承載力確定方法的研究提供了一條有效的研究途徑，這也正是本文研究的重要內(nèi)容之一.

為此，本文將從臨坡條形基礎(chǔ)雙層粘土地基破壞模式研究入手，在考慮基礎(chǔ)埋深和臨坡距離等影響的基礎(chǔ)上，采用多滑塊組合破壞模式，利用上限極限分析方法，構(gòu)建臨坡條形基礎(chǔ)雙層粘土地基機(jī)動(dòng)允許的速度場(chǎng)，建立新的臨坡條形基礎(chǔ)雙層粘土地基極限承載力分析模型，以期完善臨坡地基承載力分析的理論與方法.

1臨坡雙層粘土地基破壞模式

2機(jī)動(dòng)允許速度場(chǎng)的建立

臨坡雙層粘土地基破壞機(jī)構(gòu)機(jī)動(dòng)允許速度場(chǎng)的構(gòu)建是采用極限上限分析方法確定臨坡地基承載力必須解決的關(guān)鍵問題，為了確定允許速度場(chǎng)，首先做如下假定：

1）臨坡地基為水平層狀地基，持力層范圍內(nèi)有兩層飽和粘土，采用不排水剪有效粘聚力，內(nèi)摩擦角φ1=φ2=0.

2）條形基礎(chǔ)基底與地基土間允許產(chǎn)生相對(duì)滑動(dòng)，基礎(chǔ)以速度VP垂直向下運(yùn)動(dòng).

3）臨坡粘土地基服從MohrCoulomb屈服準(zhǔn)則和相關(guān)聯(lián)流動(dòng)法則.

在上述假設(shè)基礎(chǔ)上，根據(jù)速度相容的閉合三角形幾何關(guān)系，可得到如圖 2所示速度關(guān)系和圖3所示破壞機(jī)構(gòu)的機(jī)動(dòng)允許速度場(chǎng)[16]，再由矢量分析原理及速度相容性要求可以推導(dǎo)出各速度間斷面上速度.

本文采用多滑塊組合分析模型，一般來說塊體劃分越小，剛性滑塊的數(shù)量越多，結(jié)果越接近于真實(shí)解，取20個(gè)滑塊時(shí)能夠獲得較好的上限解.根據(jù)本文方法計(jì)算該臨坡條形基礎(chǔ)雙層粘性地基工程算例，得到臨坡地基極限承載力Qu=128.6 kPa，結(jié)果表明本文方法計(jì)算結(jié)果與有限元數(shù)值分析結(jié)果相差較小，本文上限解是較好的上限解，驗(yàn)證了本文上限分析方法的可行性與適用性.同時(shí)繪制本文分析方法對(duì)應(yīng)極限狀態(tài)下的破壞模式如圖6所示，與數(shù)值模擬分析所得潛在滑裂面接近，說明了本文臨坡雙層地基多滑塊組合破壞模型的有效性.

5.2理論計(jì)算與其他上限方法的對(duì)比

為了進(jìn)一步驗(yàn)證本文承載力公式和優(yōu)化算法的正確性與優(yōu)越性，采用以下算例，上下粘土層參數(shù)取相同值，研究問題變?yōu)閱螌油恋呐R坡地基承載力問題，與同樣采用單側(cè)滑移破壞模式，引入上限極限理論計(jì)算臨坡地基極限承載力的文獻(xiàn)[13-15]進(jìn)行比較.某臨坡條形基礎(chǔ)地基，基礎(chǔ)寬度b=2 m，基礎(chǔ)距坡頂距離與基礎(chǔ)寬度比a=0，基礎(chǔ)埋深H=2 m，地基土體重度γ=18 kN/m3，土體內(nèi)摩擦角φ=30°，地基土粘聚力c=20 kPa，邊坡角η=20°.計(jì)算得到臨坡條形基礎(chǔ)地基極限承載力Qu，并與文獻(xiàn)[13-15]的計(jì)算結(jié)果進(jìn)行比較分析，如表1所示.

本文與以上3種文獻(xiàn)都利用了上限極限分析理論和優(yōu)化方法對(duì)臨坡地基極限承載力進(jìn)行求解，通過比較計(jì)算結(jié)果可以看出：本文計(jì)算所得的地基極限承載力Qu數(shù)值略小，是較優(yōu)的上限解，在工程實(shí)際應(yīng)用中偏于安全；同時(shí)幾種方法的結(jié)果較接近，這也進(jìn)一步說明了本文承載力公式和優(yōu)化算法的正確性和合理性.

6結(jié)論

1）本文從臨坡地基出發(fā)，結(jié)合臨坡成層土地基臨近邊坡的非對(duì)稱性特征和層狀地基特征，開展臨坡條形基礎(chǔ)雙層粘土地基承載力性狀的研究，構(gòu)建了臨坡地基多滑塊組合單側(cè)滑移破壞模式，為臨坡成層土地基承載力確定奠定了堅(jiān)實(shí)的理論基礎(chǔ).

2）在臨坡條形基礎(chǔ)雙層粘土地基多滑塊破壞模式的基礎(chǔ)上構(gòu)建對(duì)應(yīng)機(jī)動(dòng)允許速度場(chǎng)，引入上限極限分析理論和優(yōu)化方法建立了可考慮地基破壞模式非對(duì)稱性、地基土成層性狀、基礎(chǔ)埋深和基礎(chǔ)與坡頂距離等影響的臨坡地基極限承載力的確定方法.

3）結(jié)合工程算例，通過與有限元數(shù)值分析方法及其他上限極限分析方法的對(duì)比分析，驗(yàn)證了本文臨坡雙層粘土地基極限承載力上限解是較好的上限解，其滑動(dòng)破壞面接近數(shù)值模擬潛在滑動(dòng)面，表明本文方法的可行性與合理性.

參考文獻(xiàn)

[1]HANSEN J B. A revised and extended formula for bearing capacity[J].Bulletin of the Danish GeoTechnical Institute，1970（28）：5-11.

[3]HANNA A M，MEYERHOF G G.Design charts for ultimate bearing capacity of foundations on sand overlying soft clay[J].Canadian Geotechnical Journal， 1980，17（2）： 300-303.

[4]楊永新，高建紅，高建明，等.雙層地基極限承載力的數(shù)值分析[J].包頭鋼鐵學(xué)院學(xué)報(bào)，2003，22（2）：181-185.

YANG Yongxin， GAO Jianhong， GAO Jianming，et al.Numerical analysis about the ultimate bearing capacity of twolayer foundation[J].Journal of Baotou University of Ironand Steel Technology， 2003， 22（2）：181-185.（In Chinese）

[5]GOURVENCE S， RANDOLPH M.Effect of strength nonhomogeneity on the shape of failure envelopes for combined loading of strip and circular foundations on clay[J].Geotechnique，2003，53（6）：575-586.

[6]BURD H J，F(xiàn)RYDMAN S.Bearing capacity of planestrain footings on layered soils[J].Canadian Geotechnical Journal，1997，34（2）：241-253.

[7]FLORKIEWICZ A.Upper bound to bearing capacity of layered soils[J].Canadian Geotechnical Journal，1989，26（4）：730-736.

[8]MICHALOWSKI R L，SHI L.Bearing capacity of footings over twolayer foundation soils[J].Journal of Geotechnical Engineering， ASCE，1995，121（5）：421-428.

[9]MERIFIELD R S，SLOAN S W，YU H S.Rigorous solutions for the bearing capacity of two layered clay soils[J].Geotechnique，1999，49（4）：471-490.

[10]CHEN W F.Limit analysis and soil plasticity[M]. Amsterdam：Elsevier，1975：170-175.

[11]秦會(huì)來，黃茂松.雙層地基極限承載力的極限分析上限法[J].巖土工程學(xué)報(bào)， 2008，30（4）：611-616.

QIN Huilai， HUANG Maosong.Upperbound method for calculating bearing capacity of strip footings on twolayer soils[J].Chinese Journal of Geotechnical Engineering， 2008，30（4）：611-616. （In Chinese）

[12]陳昌富，曹虹，王純子，等.下臥傾斜硬層對(duì)非均布荷載下地基沉降的影響[J].湖南大學(xué)學(xué)報(bào)：自然科學(xué)版，2014，41（6）：77-84.

CHEN Changfu，CAO Hong，WANG Chunzi，et al. Influnence of inclined bedrock on subgrade settlements under nonuniformly distributed load [J].Journal of Hunan University：Natural Sciences， 2014，41（6）：77-84. （In Chinese）

[13]胡衛(wèi)東，曹文貴.基于非對(duì)稱破壞模式的臨坡地基承載力上限極限分析方法[J].中國公路學(xué)報(bào)，2014，27（6）：1-9.

HU Weidong，CAO Wengui.Upper limit analysis method for ultimate bearing capacity of ground foundation adjacent to slope based on asymmetry failure mode[J].China Journal of Highway and Transport， 2014，27（6）：1-9. （In Chinese）

[14]趙煉恒，羅強(qiáng).臨坡條形基礎(chǔ)極限承載力上限計(jì)算[J].武漢理工大學(xué)學(xué)報(bào)，2010，34（1）：84-87.

ZHAO Lianheng， LUO Qiang.Bearing capacity of strip footing adjacent to slope with upper bound theorem[J]. Journal of Wuhan University of Technology， 2010，34（1）： 84-87. （In Chinese）

[15]楊峰，陽軍生，張學(xué)民，等.斜坡地基單側(cè)滑移破壞模式及承載力上限解[J].工程力學(xué)，2010， 27（6）：162-168.

YANG Feng， YANG Junsheng， ZHANG Xuemin， et al.Oneside slip failure mechanism and upper bound solution for bearing capacity of foundation on slope[J]. Engineering Mechanics，2010， 27（6）： 162-168. （In Chinese）

[16]SOUBRA A H.Upperbound solutions for bearing capacity of foundations[J]. Journal of Geotechnical and Geoenvironmental Engineering，ASCE，1999，125（1）：59-68.

[17]秦會(huì)來，黃茂松，李峰.雙層粘土地基極限承載力的多塊體上限解[J].建筑科學(xué)，2012，28（9）：40-43.

QIN Huilai，HUANG Maosong，LI Feng.Mutiblock upper bound solutions for the twolayer clay foundations[J].Building Science，2012，28（9）：40-43. （In Chinese）

[18]陳飛，練繼建，王海軍.淺埋圓形基礎(chǔ)豎向地基承載力極限分析上限解[J].湖南大學(xué)學(xué)報(bào)：自然科學(xué)版，2014，41（6）：92-98.

CHEN Fei，LIAN Jijian，WANG Haijun.Upperbound limit analysis of the vertical bearing capacity of circular shallow foundations[J].Journal of Hunan University：Natural Sciences， 2014，41（6）：92-98. （In Chinese）

[19]陳祖煜.土力學(xué)經(jīng)典問題的極限分析上、下限解[J].巖土工程學(xué)報(bào)，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. （In Chinese）