汪婧

基于上限分析原理的巖溶樁基破壞模式與極限承載力計算

汪婧

(中鐵十一局集團有限公司，湖北 武漢 430061)

針對下伏溶洞頂板極限承載力問題，提出等截面樁與階梯型變截面樁的頂板極限承載力計算方法?；跇O限上限法引入摩爾庫倫準則，構建圓臺狀頂板破壞滑裂面，提出巖溶等截面樁的樁頂荷載與頂板厚度的關系式并分析式中參數(shù)對樁頂荷載與頂板厚度曲線的影響規(guī)律，同時結合算例對比數(shù)值解和解析算法的結果，論證該算法的可靠性。在此基礎上，結合極限分析法推導階梯型變截面樁各變階處滑裂面的上限解，并尋找最先滿足上限解的滑裂面，提出該溶洞頂板破壞模式并推導出樁頂荷載與頂板厚度的關系式，通過模型試驗對比驗證計算方法的有效性。該計算方法對于上述2類截面樁的極限承載力計算均較適用，對巖溶區(qū)樁基設計有一定參考價值。

巖溶樁基；極限上限定理；極限承載力；階梯型變截面樁

隨著我國西部大開發(fā)戰(zhàn)略的實施，西部地區(qū)的交通事業(yè)快速發(fā)展。但西部部分地區(qū)地質地貌的特殊性，溶洞、采空區(qū)等特殊地質的存在給工程設計和施工帶來了很大的影響。由于工程選址時不可能完全避開巖溶區(qū)，所以對巖溶的形成和承載力的研究顯得十分必要[1?3]。極限分析理論最早由Fellenius等人在20世紀20年代提出，從50年代以來得到了較大的發(fā)展。近年來，部分學者將極限分析理論應用于下伏空洞的樁基承載力計算研究中，劉輝等[4]利用極限分析上限法推導出了圓形空洞上方條形基礎的地基極限承載力解析公式。趙明華等[5]還利用上限有限元法計算出多種工況下的地基極限承載力上限解。雷勇等[6]引入格里菲斯非線性巖石強度準則，采用極限分析方法對頂板沖切破壞模式極限承載力計算方法進行了分析。目前對于巖溶樁基的承載穩(wěn)定性問題研究主要以數(shù)值模擬[7]以及模型試驗為主[8]，理論解析方面的研究相對較少。基于格里菲斯準則的極限分析方法在進行巖溶樁基極限承載力的確定時主要考慮脆性完整巖石材料，但實際工程中的巖層并不完全發(fā)生脆性破壞，該準則存在一定的局限性，而且目前的研究對象主要以等截面樁為主，而針對階梯型變截面樁的研究則少有涉及。隨著巖溶地層階梯型變截面樁應用的逐漸興起[9]，對其極限承載力的研究就顯得極為重要?；谏鲜鲅芯康牟蛔?，本文將極限分析法與摩爾?庫倫屈服準則相結合，首先考慮巖溶地層等截面樁基破壞模式，對其樁基破壞問題進行上限分析，得到巖溶等截面樁的溶洞頂板破壞模式及頂板厚度的關系式。在此基礎上，將該計算方法推廣到階梯型變截面樁的溶洞破壞模式的上限分析之中，從而得到巖溶階梯型變截面樁的溶洞頂板破壞模式及頂板厚度的關系式。

1 等截面樁下伏溶洞頂板破壞模式上限分析

1.1 基于摩爾庫倫準則的極限分析上限定理

將式(2)和(3)代入式(1)有：

由式(4)可以說明，在任意一個滿足運動容許的速度場內，物體所處的真實荷載小于等于由虛功原理所解出的荷載。

1.2 巖溶頂板破壞模型

溶洞頂板破壞近似形成由圓臺體母線發(fā)展的破壞面[10]。當溶洞較小時，頂板范圍內滑裂面無法自由發(fā)展，破壞邊界從樁端延伸至溶洞邊界，如圖1(a)?？紤]下伏溶洞極為復雜且在地層中分布范圍大，樁徑往往會遠小于溶洞大小，因此本文統(tǒng)一假定為相對樁徑大得多的空溶洞，將模型進行簡化，忽略溶洞邊界對破壞的影響，簡化成如圖1(b)所示。

(a) 模型1；(b) 模型2

根據(jù)極限上限法假定破壞模式，并考慮對稱性構建一種由剛性塊體構成的破壞機制，考慮一簡單均質溶洞頂板的巖溶樁基破壞模式，取對稱的半樁進行分析，如圖2(a)所示。樁端上部荷載為/2，樁端距溶洞頂板高度為，樁端下巖溶頂板臨界滑裂面與豎直方向的夾角為，巖層內摩擦角為，黏聚力為，并做出以下假設：

1) 巖溶樁基頂板發(fā)生沖切破壞，分析屬于平面問題；

2) 樁為嵌巖樁，樁端下方為均質巖體，且允許向下移動，破壞范圍在溶洞頂板范圍以內；

3) 巖層破壞服從Mohr-Coulomb屈服準則，且假定溶洞相對于樁徑大得多，忽略溶洞邊界對破壞的影響。

圖2 極限平衡破壞模式和速度模式

如圖2(a)所示，頂板破壞區(qū)分為2塊，一塊是樁端下方矩形巖層，跟隨樁基一起沉降，另一塊是滑裂面邊三角形區(qū)域。將2塊區(qū)域看做剛體，那么相互之間就會產(chǎn)生相對速度及塑性流動區(qū)。假定樁體的初速度為0，方向豎直向下，則樁端下方巖層速度也為0，1表示三角形區(qū)域巖層的向溶洞移動的速度，其相對速度可用0, 1表示，在相鄰剛性塊之間，存在著速度躍變線，表示躍變線與速度方向之間的夾角，幾個速度之間的關系用矢量三角形表示，如圖2(b)所示。

由速度矢量圖可以得到1和0,1與樁端下巖溶頂板臨界滑裂面與豎直方向的夾角，巖層內摩擦角為以及初速度0之間的關系：

1.3 巖溶頂板破壞模式

考慮溶洞頂板自重，外功率包括外部荷載/2所做的功率、矩形巖層以及靠近滑裂面三角形巖層重力所做的功率。

外部荷載/2所做的功率，為樁端上部荷載/2與初速度0的乘積：

巖層的外功率為其重力與沿重力方向的分速度的乘積，矩形巖層與靠近滑裂面三角形巖層產(chǎn)生的外功率分別為：

總的外功率為三者之和，其表達式如下：

內部損耗功率為2剛性塊之間滑裂面的損耗功率，滑裂面上的能量損耗表示為：

滑裂面上的能量損耗可通過建立如圖3所示的坐標系，并通過對進行積分得到：

其中()為的直線方程：

則總的內部損耗功率為滑裂面和上的能量損耗之和，表達式如下：

根據(jù)極限上限定理，外部功率等于內部損耗功率，即

對上式進行整理化簡，可得上部荷載關于溶洞頂板厚度的目標函數(shù)：

當樁端的上部荷載遠大于巖層重力時，可忽略巖層自重對溶洞頂板破壞的影響。在不考慮土的自重時，外功率只有外部荷載/2所做的功率，推理過程與考慮荷載的情況相似，此處不再贅述，其相應的目標函數(shù)為：

1.4 模型參數(shù)的影響分析

采用MATLAB軟件編寫相應的計算程序進行分析。選取巖溶區(qū)常見的石灰?guī)r為研究基礎，其巖層容重=25 kN/m3，泊松比=0.2，保持二者在研究參數(shù)過程中不變。圖4為樁徑對上部荷載與頂板厚度曲線的影響。保持研究參數(shù)和不變，改變樁徑(分別取為1，2，3和4 m)。如圖所示，相同溶洞頂板安全厚度的條件下，隨著樁徑的增大，承載能力增大，能承受的上部極限荷載也隨之增大；相同的上部荷載的條件下，隨著樁徑的增大，承載能力也相應增大，所需要的溶洞頂板安全厚度隨之減小。

圖4 不同樁徑下溶洞頂板安全厚度與上部極限荷載的關系曲線

圖5為摩擦因數(shù)對上部荷載與頂板厚度曲線的影響。保持研究參數(shù)和不變，改變摩擦因數(shù) (分別取為10°，20°，30°，40°和50°)。由圖可知，相同溶洞頂板安全厚度的條件下，隨著內摩擦角的增大，極限承載能力也隨之提高；相同的上部荷載的條件下，隨著內摩擦角的增大，所需要的溶洞頂板安全厚度隨之減小。并且隨著溶洞頂板安全厚度增大，上部荷載的極限值增加幅度也越高，內摩擦角越小其曲線的增幅趨勢越平緩，可見，隨著巖層的增大，破壞時的剪切破壞面也會隨著增大，內摩擦角的抗破壞作用也會隨之展現(xiàn)出來。

圖5 不同內摩擦角下溶洞頂板安全厚度與上部極限荷載的關系曲線

圖6 不同黏聚力下溶洞頂板安全厚度與上部極限荷載的關系曲線

圖6為黏聚力對上部荷載與頂板厚度曲線的影響。保持研究參數(shù)和不變，改變黏聚力(分別取為10，20，30和40 MPa)。由圖可知，相同溶洞頂板安全厚度的條件下，隨著巖體的黏聚力的增大，承載能力增大，能承受的上部極限荷載也隨之增大；相同的上部荷載的條件下，隨著巖層的黏聚力的增大，承載能力也相應增大，所需要的溶洞頂板安全厚度隨之減小。

表1 計算參數(shù)

1.5 算例驗證

利用式(16)和(17)分別求出各對應溶洞頂板安全厚度下的上部荷載，同時引入文獻[8]中的計算方法，計算圍巖及樁體參數(shù)取值如表1所示。各對比曲線如圖7所示。由圖可見，溶洞頂板安全厚度在0~2 m范圍內，本文采用的方法與文獻[8]中的方法有較好的相似性。當增大時，與文獻[8]方法之間就有了較為明顯的偏差，相較于文獻[8]，由于考慮了嵌巖深度的因素，在較小時，相同厚度的溶洞頂板安全下，能夠承受的上部荷載極限值會較大，而隨著的增大，其計算結果卻相較于其他幾種方法較低，這是因為文獻[8]中的安全系數(shù)考慮較高，安全儲備預留較多。然而通過數(shù)值模擬的結果顯示，本文的方法較為接近，相較于文獻中的方法既不會偏于保守，也不會過分考慮溶洞邊界對其的影響。

圖7 不同計算方法的溶洞頂板安全厚度H與上部極限荷載P的關系曲線

2 階梯型截面樁下伏溶洞頂板破壞模式上限分析

2.1 階梯型截面樁下伏溶洞頂板的破壞模式

參考等截面巖溶樁基的溶洞頂板破壞模式，構建如圖8所示的階梯型變截面樁溶洞頂板破壞 模式。

圖8 階梯型變截面樁溶洞頂板潛在破壞模式

類比普通等截面樁在樁端樁周處發(fā)展滑裂面，假定存在3個潛在滑裂面Ⅰ，Ⅱ和Ⅲ，分別是從每級變階邊緣以及樁端邊緣開始破壞并發(fā)展至溶洞頂板的臨空面，則變截面空心樁的巖溶頂板破壞模式可能由滑裂面Ⅰ，Ⅱ和Ⅲ中其中1條構成。

根據(jù)式(18)可得到各滑裂面的上限分析方程，以滑裂面Ⅱ和Ⅲ為例：

為研究階梯型變截面樁最有可能的溶洞頂板滑裂面，將式(19)和(20)相除，得

則將式(20)改寫為

同理可以獲得任意2滑裂面上限解的比值，根據(jù)實際工程中的相應位置的上部荷載比值來推斷先發(fā)生破壞的滑裂面，最終尋找出最先發(fā)生破壞的滑裂面，即為階梯型變截面樁溶洞頂板的破壞 模式。

2.2 階梯型變截面樁下伏溶洞頂板破壞模型

當階梯型變截面空心樁的巖溶頂板破壞發(fā)生在樁端底面，最后一階的底端，即可構建如圖9所示的破壞模式，虛線表示溶洞頂板因上部荷載產(chǎn)生的滑裂面，計算假定同等截面樁。

圖9 階梯型變截面溶洞頂板破壞模式

根據(jù)極限上限法假定破壞模式，取對稱的半樁進行分析，如圖10(a)所示。樁端上部荷載為/2，樁端距溶洞頂板高度為，樁端下巖溶頂板臨界滑裂面與豎直方向的夾角為，土體內摩擦角為，黏聚力為。

由于存在變階面的關系，部分荷載從變階處傳遞到周邊土體，實際引起溶洞頂板破壞的荷載為傳至樁端的豎向荷載。結合如圖10(b)所示的速度矢量分析圖，通過極限上限定理，使得外部功率等于內部損耗功率，可得

2.3 樁端荷載的計算方法

從而樁端荷載可表示為：

圖10 極限平衡破壞模式與速度模式

圖11 2次變截面樁計算模型圖

圖11所示為2次階梯型變截面樁的細部尺寸圖，則根據(jù)樁基極限承載力的計算公式，樁側摩阻力可表示為：

其中：q為樁側土層的極限側阻力標準值。

變階處阻力可表示為：

其中：q極限端阻力標準值。

2.4 n階變截面樁P(H)函數(shù)構建

通過上述推導，將其推廣到階變截面樁的情況，如圖12所示，由于變截面樁的特殊性和下伏溶洞地層的影響下，樁周土體會存在松動現(xiàn)象，無法完全發(fā)揮極限側摩阻力及變階處端阻力，因此引入折減系數(shù)1和2，其值可由數(shù)值模擬或模型試驗確定。

根據(jù)上述的推導過程，樁側摩阻力和變階處阻力可表示為：

圖12 n階變截面樁計算模型圖

2.5 算例驗證

針對階梯型變截面樁的溶洞頂板破壞模式的研究，開展了多組模型試驗，將試驗樁的幾何參數(shù)及試驗參數(shù)代入式(32)進行極限荷載的計算。圖13為試驗結果與理論計算曲線。由圖中曲線可知，理論與試驗結果較為相近，總體趨勢較為吻合。

圖13 試驗結果與理論計算對比曲線

3 結論

1) 針對基于極限上限法推導的上部荷載關于溶洞頂板厚度的目標函數(shù)進行參數(shù)分析，溶洞巖層黏聚力，內摩擦角和樁的直徑對上部極限荷載均為正相關，有利于提高巖溶樁基的極限承載力。

2) 本文提出的方法所需參數(shù)較少，并且忽略溶洞邊界對破壞的影響，提供了一定的安全儲備，對工程設計有很大的便利之處，能對初步設計提供較好的參考。

3) 通過對巖溶等截面樁基溶洞頂板破壞的研究分析，類比推理階梯型變截面空心樁溶洞頂板潛在滑裂面，提出以樁端處延伸至溶洞的階梯型變截面樁的溶洞頂板破壞模式，并得到了巖溶等截面樁以及階梯型截面樁的頂板極限承載力計算方法。

[1] Waltham A C, Swift G M. Bearing capacity of rock over mined cavities in Nottingham[J]. Engineering Geology, 2004, 75(1):15?31.

[2] Hatzor Y H, Wainshtein I, Mazor D B. Stability of shallow karstic caverns in blocky rock masses[J]. International Journal of Rock Mechanics & Mining Sciences, 2010, 47(8):1289?1303.

[3] JIANG C, LIU L, WU J P. A new method determining safe thickness of karst cave roof under pile tip[J]. Journal of Central South University, 2014, 21(3):1190?1196.

[4] 劉輝, 楊峰, 陽軍生. 空洞上方淺基礎地基破壞模式與極限承載力分析[J]. 巖土力學, 2010, 31(11): 3373? 3378.LIU Hui, YANG Feng, YANG Junsheng. Failure mechanism and ultimate bearing capacity of strip footing above void[J]. Rock and Soil Mechanics, 2010, 31(11): 3373?3378.

[5] 趙明華, 盧曉明, 張銳. 下伏空洞地基極限承載力與破壞模式上限有限元法[J]. 巖土力學, 2017, 38(1): 229?236, 246.ZHAO Minghua, LU Xingming, ZHANG Rui. Upper bound finite element method for ultimate bearing capacity and failure mechanism of subgrade above void[J]. Rock and Soil Mechanics, 2017, 38(1): 229?236, 246.

[6] 雷勇, 尹君凡, 陳秋南, 等. 基于極限分析法的溶洞頂板極限承載力研究[J]. 巖土力學, 2017, 38(7): 1926? 1932.LEI Yong, YIN Junfan, CHEN Qiunan, et al. Study on the ultimate bearing capacity of karst roof based on limit analysis[J]. Geotechnical Mechanics, 2013, 38(7): 1926? 1932.

[7] 王華牢, 張鵬, 李寧. 巖溶洞穴對嵌巖單樁承載力的影響研究[J]. 西安理工大學學報, 2010, 26(1): 31?36.WANG Hualao, ZHANG Peng, LI Ning. Research on the single rock-socketed pile bearing capacity influenced by karst caves[J]. Journal of Xi’an University of Technology, 2010, 26(1): 31?36.

[8] 黃明, 付俊杰, 陳福全, 等. 樁端巖溶頂板的破壞特征試驗與理論計算模型研究[J]. 工程力學, 2018, 35(10): 172?182.HUANG Ming, FU Junjie, CHEN Fuquan, et al. Theoretical calculation model and model test on the failure characteristic of karst roof under rock-socketed pile[J]. Rock and Soil Mechanics, 2018, 35(10): 172? 182.

[9] 黃明, 江松, 許德祥, 等. 超大直徑變截面空心樁的荷載傳遞特征與理論模型[J]. 巖石力學與工程學報, 2018, 37(10): 2370?2383.HUANG Ming, JIANG Song, XU Dexiang, et al. Load transfer mechanism and theoretical model of step tapered hollow pile with huge diameter[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(10): 2370? 2383.

[10] 雷勇, 劉一新, 鄧加政, 等. 沖切破壞模式下溶洞頂板極限承載力計算[J]. 巖石力學與工程學報, 2018, 37(9): 2162?2169. LEI Yong, LIU Yixin, DENG Jiazheng, et al. Calculation of ultimate bearing capacity of cave roof under punching failure mode[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(9): 2162?2169.

Failure mode and ultimate bearing capacity of karst pile foundation based on upper-bound theorem

WANG Jing

(China Railway 11 Bureau Group Co. Ltd., Wuhan 430061, China)

In view of the problem of the ultimate bearing capacity of the roof of the underlying karst cave, the calculation method of the ultimate bearing capacity of the roof of the equal section pile and the stepped variable section pile was put forward. Mohr coulomb criterion, introduced based on the upper-bound theorem, build the frustum of a cone shaped roof damage critical slip surface, put forward such as karst section pile top load and the relation between the thickness of roof and analytical parameters on pile top load and the influence law of roof thickness curve, at the same time, combined with the numerical solution of examples and the result of the parsing algorithm, proved the reliability of the algorithm. On this basis, combined with the limit analysis method, the upper limit solution of sliding surface at each variable order of stepped type variable section pile was deduced, and the first sliding surface that meets the upper limit solution was found, the failure mode of roof of karst cave was put forward, and the relationship between pile top load and roof thickness was derived. The calculation method was suitable for the calculation of the ultimate bearing capacity of the above two kinds of cross-section piles and has certain reference value for the design of karst area pile foundation.

karst pile foundation; upper-bound theorem; ultimate bearing capacity; stepped variable section pile

TU473

A

1672 ? 7029(2019)09? 2207 ? 08

10.19713/j.cnki.43?1423/u.2019.09.011

2018?11?14

國家重點研發(fā)計劃項目(2018YFC0808704)

汪婧(1982?)，女，江西景德鎮(zhèn)人，高級工程師，從事工程技術與管理相關工作；E?mail：8295161@qq.com

(編輯 涂鵬)