基于多塊體的非對稱雙層臨坡地基上限承載力

楊博銘

摘 ? 要：為進一步解決條形基礎(chǔ)荷載作用下雙層臨坡地基的承載力和穩(wěn)定性問題，對現(xiàn)有研究成果進行對比分析，并對其假設(shè)的合理性進行判斷，根據(jù)臨坡地基破壞模式的非對稱性的基本特征，考慮到雙層地基分層界面處速度的間斷性，采用一種新的多滑塊模式來構(gòu)建其破壞模型. 根據(jù)速度相容關(guān)系、速度三角形閉合條件確定該模型的相容速度場. 引入極限分析法的上限分析理論，推導(dǎo)條形基礎(chǔ)非對稱雙層臨坡地基上限承載力的計算公式. 采用MATLAB二次規(guī)劃算法的優(yōu)化計算方法對雙層臨坡地基的上限承載力進行優(yōu)化求解. 最后通過工程算例來驗證理論公式的合理性，并與現(xiàn)有研究方法以及軟件計算結(jié)果進行對比，證明了本文方法的可行性.

關(guān)鍵詞：極限承載力;雙層地基;臨坡地基;非對稱破壞模式

中圖分類號：U443.15 ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼：A

Abstract：In order to solve the problems of bearing capacity and stability for the double-layer slope foundation under the load of the strip foundation，this paper compared and analyzed the previous research results，and judges the rationality of its hypothesis. According to the basic characteristics of the asymmetry of the failure mode of the slope foundation and considering the velocity inequality at the stratified interface of the double-layer foundation， a new multi-slider failure mode was used to build its destruction mode. The compatibility velocity field was determined according to the speed-compatible relation and velocity triangle closure condition. At the same time，the upper limit analysis theory of the limit analysis method was introduced to derive the calculation formula of the upper limit bearing capacity of the asymmetric double-layer slope foundation of the strip foundation. The optimization calculation method of the secondary planning algorithm in MATLAB was used to optimize and solve the upper limit bearing capacity of double-layer slope foundation. Finally，the engineering examples was used to verify the rationality of the proposed method， and compared with the existing research methods and the results of the software calculation，proving the feasibility of the method used in this paper.

Key words：ultimate bearing capacity;double layered foundation;slope foundation;asymmetric failure mode

臨坡地基承載力，既可看作地基極限承載力問題，也可看作邊坡整體穩(wěn)定性問題. 臨坡地基與水平地基的區(qū)別在于：臨坡地基的承載力受到坡體的幾何特性以及邊坡土體的物理力學(xué)性質(zhì)的共同影響和制約. 因此推導(dǎo)能綜合考慮臨坡地基影響因素的承載力計算方法是解決工程實際問題的有效途徑.

目前，在進行臨坡地基設(shè)計時，沒有相關(guān)規(guī)范可以參考，絕大多數(shù)工程問題都是采用在水平地基極限承載力公式的基礎(chǔ)上乘以一個修正系數(shù)的方法[1]. 這樣的經(jīng)驗計算方法的可靠性、合理性以及經(jīng)濟性在一定程度上很難實現(xiàn)定量把控[2-3].關(guān)于臨坡地基極限承載力研究，可分為以下三類：第一類是在平地地基的基礎(chǔ)上進行研究[4-6]，考慮坡面傾斜系數(shù)來進行折減修正，其計算簡單，力學(xué)理論清晰明了，但是折減系數(shù)的確定沒有理論公式支持，僅根據(jù)經(jīng)驗取值，缺乏一定的科學(xué)性. 第二類就是直接以臨坡坡體為研究對象，充分考慮其破壞模式的非對稱性，但是其理論基礎(chǔ)還是固體介質(zhì)力學(xué)，這樣導(dǎo)致計算結(jié)果與工程實際存在一定差距[7-8].第三類就是通過實驗來處理工程問題[9-11].雖然實驗的結(jié)果與工程實際更加接近，但是造價以及時間成本是制約其發(fā)展的主要因素. 上述極限承載力計算方法，在工程應(yīng)用中計算值可能過于保守，造成浪費，或者因承載力不足而導(dǎo)致工程事故. 綜上，完善臨坡地基設(shè)計理論、計算方法，在保證建筑物安全的前提下，優(yōu)化工程設(shè)計參數(shù)、降低工程造價具有重要的經(jīng)濟效益和社會效益.

在臨坡地基破壞機理、破壞模式研究方面，Kusakabe等[12]、Mizuno等[13]認為臨坡地基單側(cè)破壞模式的破壞機構(gòu)由3個滑塊組成，即基礎(chǔ)下的三角形主動破壞區(qū)、中間過渡區(qū)和被動滑塊區(qū)，且認為基礎(chǔ)底下的主動區(qū)為對稱的三角形. 但是，對于中間過渡區(qū)，不同學(xué)者又有不同的觀點，Kusakabe認為是對數(shù)螺旋線而Castelli等認為是圓弧滑動面并延伸到坡面[14];對于主動破壞區(qū)，王曉謀等[15]利用有限元法提出基礎(chǔ)底下的彈性區(qū)應(yīng)為不對稱的三角形. 同時魏學(xué)勇等[16]還提出應(yīng)該考慮基礎(chǔ)內(nèi)側(cè)土體應(yīng)力分布對基礎(chǔ)承載力的影響.

對于非均質(zhì)的臨坡地基承載力的研究，部分學(xué)者給出了相應(yīng)的研究方法. 胡衛(wèi)東等[17]考慮了臨坡地基的分層性以及非對稱性，但是沒有考慮到內(nèi)摩擦角對臨坡地基承載力的影響，也沒有考慮基礎(chǔ)與土體之間的相對滑動對承載力的影響;秦會來等[18]考慮了地基的分層特性，但其計算方法是基于水平地基的. 因此，本文在條形基礎(chǔ)荷載作用下，將多塊體理論運用于臨坡地基，考慮臨坡地基破壞模式的非對稱性與地基土體的成層性，同時考慮了土層分層界面處上、下土層的速度間斷性，以及基礎(chǔ)與土體之間的相對滑動，來建立臨坡地基的多塊體破壞模型及其相應(yīng)允許的速度矢量場. 引入極限分析理論，推導(dǎo)了非對稱雙層臨坡地基上限承載力計算方法.

1 ? 臨坡雙層地基破壞模式的確定

1.1 ? 多塊體上限分析法

在本文計算模型中，外力功的功率應(yīng)包括：重力功的功率、基礎(chǔ)傳遞荷載功的功率和考慮埋深的超載功率;對應(yīng)的內(nèi)能耗散的功率應(yīng)包括：塊體間速度間斷面的內(nèi)能耗散率和基礎(chǔ)底面摩擦消耗的內(nèi)能.

1.2 ? 雙層臨坡地基中多塊體上限法的破壞模式

破壞模式?jīng)Q定了相容速度場的分布從而決定了內(nèi)能耗散率. 相比于雙側(cè)破壞模式，單側(cè)破壞模式計算臨坡地基承載力偏于安全，且與工程實際更加接近. 為了能夠充分體現(xiàn)分層界面處速度的間斷性，本文采用如圖1所示的破壞模式，并建立如圖1所示的坐標系：基礎(chǔ)中心點O為坐標原點，水平向右、豎直向下分別為x軸、y軸正方向.

如圖1所示，根據(jù)臨坡地基單側(cè)破壞機理，破壞模式由n個剛性三角形組成. 為便于敘述，分別做如下規(guī)定：滑動面記為A123456…n，沿著滑動面上的節(jié)點1、2、…、n記為外部節(jié)點;外部節(jié)點與坐標原點相連的線段稱為主斜線（如O71）;水平地層分界線與主斜線的交點用n+1、n+2、2n來標記，記為內(nèi)部節(jié)點;條塊標示用①②③…表示，條形基礎(chǔ)正下方的三角形條塊為條塊①;連接主斜線的線段稱為輔斜線（如72）;位于滑裂面上的輔斜線稱為底輔斜線（如12）;滑動面上的n個大三角形（△O12、△O23…）以原點為頂角，且頂角分別設(shè)定為α1、α2、α3…;底輔斜線與豎直線的夾角（順時針）為β1、β2、β3…;破壞面上的底輔斜線長為c1、c2、c3…;輔斜線長為p1、p2、p3;主斜線長為l1、l2、l3.

在圖1中，對于水平層狀土，若主斜線li位于兩層土中則其長度li = lsi + lxi，其中l(wèi)si、lxi分別表示位于上層土的長度和位于下層土的長度. 大三角形被輔斜線分為兩個三角形，三角形編號如圖1所示，三角形的面積分別表示為S1、S2、S3… . 土層分界線將偶數(shù)編號的三角形面積分為上、下層土的面積之和即Si = Ssi + Sxri，其中Ssi為上層土中三角形的面積，Sxri為下層土中的面積. 編號為奇數(shù)的三角形都位于下層土中，Si = Sxli. 下層土中四邊形面積被輔斜線分為左右兩個三角形，分別用Sxli、Sxri表示其面積.

1.3 ? 幾何參數(shù)的求解

已知條件：條形基礎(chǔ)的寬度d、基礎(chǔ)外沿距坡頂距離l=ad（a為距離比系數(shù)）、可變參量αi和βi、上層黏土的厚度h. 根據(jù)圖1所示幾何關(guān)系，可分別求出其他幾何參數(shù).

2 ? 機動允許速度矢量場的確定

2.1 ? 基本假定

根據(jù)條形基礎(chǔ)下的臨坡地基的幾何參數(shù)，以及本文研究的具體情況，在確定臨坡雙層非對稱地基承載力上限值所允許的速度矢量時，做如下假定：

1）土體服從Mohr-Coulomb屈服準則以及相關(guān)流動法則;

2）考慮條形基礎(chǔ)基底和土體的相對滑動，其外摩擦角為δ;

3）基礎(chǔ)在上部荷載作用下向下做垂直運動;

4）土層分界面為水平向.

2.2 ? 速度的確定

根據(jù)假設(shè)以及速度相容原理，求得各滑塊的速度以及間斷面的速度.

1）根據(jù)圖2所示速度關(guān)系，在條形基礎(chǔ)下方的三角形滑塊的速度以及間斷面的速度見式（23）.

要使臨坡地基取得極限承載力Qu的最小值，只有當?shù)鼗幱跇O限狀態(tài)，即式（56）取最小值. 利用Matlab編程計算可獲得承載力的目標函數(shù)表達式，即將問題轉(zhuǎn)化為約束非線性規(guī)劃問題. 對式（56）進行優(yōu)化求解必須滿足式（22）的約束條件，本文采用序列二次規(guī)劃算法（SQP）. SQP非線性優(yōu)化算法是利用目標函數(shù)和約束條件構(gòu)造增廣目標函數(shù)，借此將約束最優(yōu)化問題轉(zhuǎn)化為無約束最優(yōu)化問題，然后利用求解無約束最優(yōu)化問題方法間接求解新目標函數(shù)的局部最優(yōu)解問題，此方法既適合高度非線性問題，也具有較快的運算速度[22-24].

5 ? 工程算例計算與分析

為了驗證本文關(guān)于非對稱雙層臨坡地基上限承載力計算方法的合理性，結(jié)合具體工程案例進行分析，并與相關(guān)的分析方法、軟件計算結(jié)果進行對比. 驗證的具體情況如下.

5.1 ? 算例1

采用文獻[17]中的計算案例：某臨坡條形基礎(chǔ)雙層地基，基礎(chǔ)埋深H = 1.0 m，持力層范圍內(nèi)兩層飽和黏土. 土體相關(guān)計算參數(shù)詳見文獻[17]. 基礎(chǔ)外邊緣到基礎(chǔ)坡頂?shù)木嚯xl = a × b，其中a為距離比例系數(shù)，a = 1，邊坡坡角η = 20°.

5.1.1 ? 計算結(jié)果對比分析

文中還采用有限差分軟件（FLAC-3D）對工程案例1的非對稱臨坡地基進行數(shù)值計算. 在分析過程中按平面應(yīng)變問題處理，材料參數(shù)：E = 45 × 103 ?kPa，ν = 0.3. 土體的本構(gòu)模型服從Mohr-Coulomb屈服準則. 圖5為分級加載的荷載-位移曲線，每一級荷載為10 kPa.

圖5所示為FLAC模擬的分級加載的荷載-位移曲線，可以看到，當加載到110 kPa時位移呈極速下降的狀態(tài)，通過觀察塑性區(qū)可看出，此時塑性區(qū)已經(jīng)貫通. 而文獻[17]中采用Abaqus得到的荷載-位移曲線得出破壞荷載也應(yīng)該在110 kPa左右，此時荷載-位移曲線已經(jīng)處于快速下降階段，在文獻[17]中作者的取值有些偏大.

數(shù)值模擬結(jié)果表明，案例1中的雙層地基加載到Qu = 110 kPa時臨坡地基發(fā)生了破壞. 文獻[17]的理論計算方法所得地基極限承載力Qu = 128.6 kPa，本文所得地基極限承載力Qu = 125.9 kPa，兩者的計算結(jié)果與數(shù)值分析結(jié)果都很接近，很明顯本文結(jié)果更加接近數(shù)值分析結(jié)果. 在多滑塊體分析模型中，塊體劃分得越細其計算值就越精確，本文的塊體劃分比文獻[17]中的塊體數(shù)量少了一半，但計算結(jié)果比文獻[17]中的計算值更加接近于數(shù)值分析結(jié)果，因此本文的分析方法可行且效率更高.

5.1.2 ? 參數(shù)分析與討論

分層地基的承載力一般會受到土層性質(zhì)、土層厚度的影響. 通常，下臥軟弱土層都對基礎(chǔ)承載力有不利影響. 為了探討上、下層土的強度參數(shù)以及幾何尺寸對基礎(chǔ)的極限承載力的影響規(guī)律，下文通過改變案例1中的下臥層土的黏聚力，其他參數(shù)保持不變，來分析土體強度參數(shù)的改變對承載力的影響;在保持強度參數(shù)不變的情況下改變上層土體的厚度來分析土層厚度對承載力的影響.

圖6所示為在不同的上層土體厚度情況下，土體黏聚力對承載力的影響曲線. 由圖6可知：隨著下層土體黏聚力的增加，臨坡條形地基的承載力上限值也逐漸增加且增加速度越來越快. 在上層土厚度不同的情況下，上、下層土體強度差距越小，其地基的承載力也就越接近，這與事實相符. 當上、下層土體強度逐漸接近時，不管上層土的厚度是多少，其地基承載力趨近于同一個值，即接近于均質(zhì)地基承載力.

圖7所示為臨坡條形基礎(chǔ)的上限承載力與上層土厚度的關(guān)系，圖中給出了上、下土層不同黏聚力的承載力隨著上層土的厚度變化的關(guān)系. 當兩層土體強度相差不是很大時，隨著上層土厚度的增加，地基承載力也增加. 當土層厚達4 m左右時，承載力不再受黏聚力比的影響，即轉(zhuǎn)化為單層土體的地基承載力形式. 兩層土體的強度差異很大時，其承載力依然隨著厚度增加而增加.

5.2 ? 算例2

為進一步驗證本文計算理論的正確性和合理性，下面采用不分層的臨坡地基來計算地基上限承載力. 一臨坡地基，基礎(chǔ)寬b = 2 mm，坡頂距比a = 0，基礎(chǔ)埋深2 m，重度γ = 18 kN/m3，內(nèi)摩擦角φ = π/6，黏聚力C = 20 kPa，坡角η = π/9. 計算結(jié)果和文獻[17，25-26]計算結(jié)果對比見表2.

算例2是非成層地基，而本文計算方法主要針對分層地基，但是本文計算方法可以退化為均勻地基來計算臨坡地基的承載力上限值. 計算結(jié)果與其他3種理論計算方法對比表明：本文采用較少的計算單元（10個計算單元）獲得了與其他文獻相近的承載力值. 綜上，本文采用較少的單元劃分，得到較好的計算結(jié)果，證明了本文采用的破壞模式是合理、有效的.

6 ? 結(jié) ? 論

本文綜合考慮了地基的分層性、破壞的非對稱性、基礎(chǔ)與土體的摩擦作用、地基埋深、基礎(chǔ)外緣與坡頂距離等對臨坡地基承載力的影響. 采用一種新的多滑塊模式來構(gòu)建其破壞模式并根據(jù)速度相容的關(guān)系、速度三角形閉合條件確定該模型的相容速度場，同時采用極限分析法的上限分析理論來推導(dǎo)雙層臨坡地基上限承載力的計算公式. 最后，通過兩個算例，分別與數(shù)值分析結(jié)果和理論計算結(jié)果對比，證明了本文計算方法的合理性，并進行參數(shù)分析，得出如下結(jié)論：

1）隨著下層土體黏聚力的增加，臨坡條形地基承載力上限值也逐漸增加且增加得越來越快.

2）當兩層土體強度相差不是很大時，隨著上層土厚度的增加，地基承載力也增加. 當土層厚達4 m左右時，承載力不再受黏聚力比的影響，即轉(zhuǎn)化為單層土體的地基承載力形式. 兩層土體的強度差異很大時，其承載力依然隨著厚度增加而增加.

參考文獻

[1] ? ?GB 50007—2011 建筑地基基礎(chǔ)設(shè)計規(guī)范[S]. 北京：中國建筑工業(yè)出版社，2011：18—33

GB 50007—2011 Code for design of building foundation [S]. Beijing：China Architecture & Building Press，2011：18—33. （In Chinese）

[2] ? ?王紅雨. 臨近基坑既有建筑物地基承載力減損性狀的研究[D].上海：同濟大學(xué)土木工程學(xué)院，2005：1—24.

WANG H Y. Study on the bearing capacity loss of existing building adjacent to foundation pit[D]. Shanghai：College of Civil Engineering，Tongji University，2005：1—24. （In Chinese）

[3] ? ?李海坤. 基坑開挖對臨近建筑物地基極限承載力的影響研究[D]. 上海：同濟大學(xué)土木工程學(xué)院，2003：1—37.

[17] ?胡衛(wèi)東，曹文貴，袁青松. 臨坡雙層黏土地基極限承載力的上限分析[J]. 湖南大學(xué)學(xué)報（自然科學(xué)版），2016，43（1）：110—116.

HU W D，CAO W G，YUAN Q S. Upper bound solution for ultimate bearing capacity of the two-layer clay foundation adjacent to slope [J]. Journal of Hunan University （Natural Sciences），2016，43（1）：110—116. （In Chinese）

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

QING H L，HUANG M S. Limit analysis of ultimate bearing capacity of double-layer foundation by upper bound method [J]. Chinese Journal of Geotechnical Engineering，2008，30（4）：611—616. （In Chinese）

[19] ?MICHLOWSKI R L. Collapse loads over two-layer foundation soils [J]. Soil and Foundations，2002，42（1）：1—7.

[20] ?張學(xué)言. 巖土塑性力學(xué)[M]. 北京：人民交通出版社，1993：232—330.

ZHANG X Y. Plastic mechanics of rock and soil [M]. Beijing：Peoples Communications Press，1993：232—330. （In Chinese）

[21] ?趙澎年. 松散介質(zhì)力學(xué)[M]. 北京：地震出版社，2010：6—36.

ZHAO P N. Mechanics of unconsoildated media [M]. Beijing：Earthquake Press，2010：6—36. （In Chinese）

[22] ?張忠楨.二次規(guī)劃-非線性規(guī)劃與投資組合的計算[M]. 武漢：武漢大學(xué)出版社，2006：152—161.

ZHANG Z Z. Quadratic programming-nonlinear programming and the calculation of portfolio [M]. Wuhan：Wuhan University Press，2006：152—161. （In Chinese）

[23] ?李晶晶，孔令偉，凌賢長. 高鐵塹坡粉質(zhì)黏土原位強度特性與應(yīng)力歷史效應(yīng)[J]. 湖南大學(xué)學(xué)報（自然科學(xué)版），2019，46（3）：99—105.

LI J J，KONG L W，LING X C. In-situ strength characteristics and stress history effect of silty soil on high-speed railway cut slope[J]. Journal of Hunan University （Natural Sciences），2019，46（3）：99—105. （In Chinese）

[24] ?曹文貴，李鵬，許烜，等. 考慮模量變化與三維應(yīng)力影響的地基沉降分析方法[J]. 湖南大學(xué)學(xué)報（自然科學(xué)版），2018，45（3）：142—148.

CAO W G，LI P，XU X，et al. An analytical method of ground foundation settlement influenced by modulus change and three-dimensional stress[J]. Journal of Hunan University （Natural Sciences），2018，45（3）：142—148. （In Chinese）

[25] ?趙煉恒，羅強，楊峰，等. 臨坡條形基礎(chǔ)極限承載力上限計算[J]. 武漢理工大學(xué)學(xué)報（交通科學(xué)與工程版），2010，34（1）：84—87.

ZHAO L H，LUO Q，YANG F，et al. Bearing capacity of strip footing adjacent to slope with upper bound theorem [J]. Journal of Wuhan University of Technology （Transportation Science & Engineering），2010，34（1）：84—87. （In Chinese）

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

YANG F，YANG J S，ZHANG X M，et al. One-side slip failure mechanism and upper bound solution for bearing capacity of foundation on slope [J]. Engineering Mechanics，2010，27（6）：162—168. （In Chinese）