非等樁徑葷素樁組合咬合樁受力計算模型研究

洪源，彭正坤，任晨灝，曾習凡，梅東，李繁

非等樁徑葷素樁組合咬合樁受力計算模型研究

洪源1，彭正坤1，任晨灝1，曾習凡1，梅東1，李繁2

(1. 中國水利水電第八工程局有限公司，湖南 長沙 410004；2. 中南大學 土木工程學院，湖南 長沙 410075)

針對咬合樁等效成地下連續墻的理論計算模型中的不足，本研究考慮鉆孔咬合樁中葷素樁之間樁徑相對大小的影響及素混凝土樁的剛度折減效應，推導出修正后的等效剛度的咬合樁受力與形變理論計算模型。依托深圳市城市軌道交通12號線太子灣站主體圍護結構咬合樁工程，采用Midas/GTS NX有限元軟件建立了數值模擬，分析了基坑開挖全過程圍護結構的受力機理與變形規律，驗證理論計算模型。研究結果表明：隨著基坑開挖，圍護結構最大水平位移最終收斂于10 mm左右。開挖工況下，圍護結構的內力與變形增長較快，尤其第三層土開挖完成后結構內力迅速增長。驗證了該理論模型能夠有效預測基坑施工過程中鉆孔咬合樁的內力與位移發展規律，可為圍護結構的合理設計與安全施工提供理論指導。

鉆孔咬合樁；等效剛度法；數值模擬；受力機理；變形規律

鉆孔咬合樁基坑圍護結構主要包括素混凝土樁與鋼筋混凝土樁，與排樁等復合樁基的不同處是咬合樁中的葷素樁相互咬合，在土體中形成了致密的止水帷幕，同時具有良好的圍護性能[1?2]。鉆孔咬合樁結構從空間形態上比較接近地下連續墻的形式。因此，在受力特點上，類似地連墻的受力模 式[3?4]。

目前，國內外研究學者多采用等效剛度的分析理論對咬合樁支護結構進行受力機理分析[5?7]。任紅林[8]等人采用等效剛度法對鉆孔咬合樁結構進行內力求解及配筋設計，證明了該方法可以滿足咬合樁設計要求。廖中林[9]等人將鉆孔咬合樁簡化成相應厚度的地連墻，再進行咬合樁的設計與驗算分析。張前[10]等人根據抗彎等效剛度原則并結合工程實例，分析咬合樁的變形與穩定時采取等效地下連續墻的力學模式。張楨[11]運用有限元數值模擬的手段，建立了鉆孔咬合樁等效剛度法的有限元模型，并與實測數據進行對比，驗證了模型的有效性。閆文進[12]以直線型咬合樁的等效剛度分析方法為基礎，調整葷素樁的咬合角度，將其設計為折線形咬合樁結構，結合理論計算法與數值模擬方法，綜合分析了樁體抗彎承載性能。本研究基于剛度等效原則，擬將鉆孔咬合樁等效為地下連續墻的計算模型，并根據實際工程中葷素樁的樁體特征差異對現有的等效地連墻的計算模型進行參數修正，建立適應拋填石層超深基坑鉆孔咬合樁的結構受力變形分析模型，并結合數值模擬[13?14]的方法，驗證其合理性。

1 咬合樁結構理論計算模型

不同樁徑的葷素樁組成的鉆孔咬合樁結構計算模型如圖1所示，其修正后的等效剛度的計算公式為：

式中：I1為鋼筋混凝土樁的截面慣性矩；I2為素混凝土樁的有效截面慣性矩；E1為鋼筋混凝土樁的混凝土彈性模量；E2為素樁的混凝土彈性模量；E3為等效地連墻的混凝土彈性模量；d1為鋼筋混凝土樁的樁徑；d2為素樁的樁徑；a為樁體咬合量；b為等效的地下連續墻厚度；k為素樁的剛度折減系數，參考文獻[8]的折減系數取值范圍，考慮該工程實際取0.8。

鋼筋混凝土樁樁體截面慣性矩為：

素混凝土樁的有效截面慣性矩為：

3、4及5的計算簡圖如圖2所示，計算式為：

式中：3為素混凝土樁的圓形截面慣性矩；4為素混凝土樁截面弓形部分的截面慣性矩；5為鋼筋混凝土樁截面弓形部分的截面慣性矩；為截面上任意點的橫坐標；為截面上任意點的縱坐標；1為咬合截面高度的一半；1為鋼筋混凝土樁的截面半徑；2為素混凝土樁的截面半徑。

將公式(4)～(6)代入公式(3)，可得到素樁有效截面慣性矩：

如圖3所示，咬合樁截面等效計算模型通過分析，確定鉆孔咬合樁等效剛度法的計算單元長度為(1+2?2)。

圖2咬合樁截面慣性矩I2計算簡圖

圖3 咬合樁截面等效計算簡圖

根據式(1)，可變換成：

根據式(8)，可求得等效墻體的等效厚度：

再完成等效墻體彎矩、剪力及位移的計算，并將其換算成咬合樁的彎矩、剪力及位移。換算的計算式為：

式中：為素樁與葷樁的剛度比；1為葷樁樁體的彎矩；1為葷樁樁體的剪力；1為葷樁樁體的位移；2為葷樁樁體的彎矩；2為葷樁樁體的剪力；2為葷樁樁體的位移；為等效墻體彎矩；為等效墻體剪力；為等效墻體位移。

2 咬合樁數值模擬

2.1 工程地質概況

以深圳市城市軌道交通12號線工程土建一工區太子灣站主體圍護結構咬合樁工程為背景?？紤]選取主體圍護結構A-A橫截面進行建模分析。此基坑深約30 m，寬度約26 m，場地地下水位埋深取0.88 m。自上而下分布的地層依次為：填碎石、素填土、淤泥質粉質黏土、粉質黏土、強風化片麻狀混合花崗巖、中風化片麻狀混合花崗巖。計算水土壓力時，采取水土合算的方法，模擬區間基坑所在地層條件如圖4所示。

基坑開挖分5層進行，土層厚度由上往下依次為6.9、3.1、6.3、5、5 m?；觾蓉Q向設5道支撐，第一、二、三、五道支撐采用800×1 000砼支撐，第四支撐采用1 000×1 200砼支撐。所有內支撐選用C30混凝土。咬合樁素混凝土樁樁徑為1 200 mm，采用C20混凝土。配筋樁樁徑1 500 mm，采用C35混凝土。葷樁與素樁中心距為1 000 mm，兩樁的咬合量為350 mm。等效地下連續墻的有效厚度為1 250 mm，采用C30混凝土。從圖4中可以看出，第一層碎填石與第二層素填土的交界線并不水平，因現場實際中碎填石層并不是均勻等厚，所以給出的地質剖面圖也只能代表某一區間段的地層情況。因此，考慮到模型的準確性與簡易性，本研究對于碎填石層和其他土層采取等效均勻厚度進行模擬分析。

圖4 模擬區間所在地層地質縱剖面

2.2 數值分析模型

采用Midas/GTS NX有限元分析軟件建立圍護結構模型進行數值分析?；迂Q直方向為Y軸，沿基坑水平方向為X軸，建立計算坐標系。為提高計算精度，利用混合四面體網格，在基坑內土體及圍護結構施工位置處，網格劃分較密；向外部土體網格劃分，逐漸變疏。在不影響計算精度的前提下，減少網格數量，節省計算時間，網格單元尺寸在1.0～1.5 m，共劃分6 626個單元，6 582個節點。

二維模型邊界的約束為：對模型側面垂直邊界施加水平向約束，底部水平邊界施加水平向約束與垂直向約束，地層表面為自由面。根據現場基坑施工的工藝流程，通過Midas/GTS NX軟件，模擬圍護結構施作和基坑開挖。在模型中，依據實際施工數據模擬地面施工超載及土體自重，依據相關資料的數據模擬圍護結構與內支撐。其中，地面施工荷載為線壓力形式，施加于基坑兩側一定范圍內。其主要步驟為：

1) 初始地應力平衡。將模型內初始狀態參與計算的各土層、邊界條件激活，同時施加模型整體水位線，使后續計算時考慮滲流，在重力作用下達到平衡狀態，并將其位移清零。

2) 圍護結構施作模擬。激活該位置處的等效地下連續墻、施工荷載。

3) 基坑內土體開挖。進行該處的內支撐施作與土體開挖過程的模擬，激活該處內支撐，并鈍化該處的開挖土體。

4) 待分層開挖的第一層土開挖過程計算平衡后，重復步驟三，開挖至預定的基坑底部，完成整個模型的計算。有限元數值模型如圖5所示。

本模型地層及主要結構材料參數見表1。

圖5 圍護結構數值模型示意

表1 模型地層參數表

2.3 理論模型與數值模型對比驗證

該研究區體圍護結構A-A橫斷面與數值模擬分析中選擇的斷面一致?；诓煌瑯稄降娜澦貥督M成的鉆孔咬合樁結構，推導出修正后的等效剛度的計算模型。將理論模型預測結果與數值模擬分析結果進行相互印證，分析兩者異同，探明拋填石層超深基坑鉆孔咬合樁的受力機理。在同一施工工況下，理論分析計算結果與數值模擬分析結果相互對比，給出圍護結構的彎矩、剪力及水平位移變化情況，如圖6～11所示。

圖6 各工況下圍護結構剪力模擬結果

圖7 各工況下圍護結構彎矩模擬結果

圖8 理論分析與數值模擬最大彎矩對比

圖6～11中的1～55種開挖工況依次代表第一層、第二層、第三層、第四層、第五層土開挖的施工工況。

圖9 理論分析與數值模擬最大剪力對比

從圖6中可以看出，對應開挖工況1～5，圍護結構所承受的最大剪力分別為214.234、255.177、395.866、337.91、324.601 kN，最大剪力一般發生在開挖深度16～20 m的范圍內。隨著開挖的進行，最大剪力的變化規律表現為先增大后降低，且剪力的增大幅度大于降低幅度。同時，圍護結構在開挖完成第三層土體時，承受到最大的剪力395.866 kN。

從圖7中可以看出，對應開挖工況1～5，圍護結構所承受的最大彎矩分別為563.924、629.942、763.608、832.241、722.982 kN·m，最大彎矩一般發生在開挖深度18～20 m的范圍內，這一變化規律與圍護結構的剪力變化情況相同。最大彎矩的變化規律表現為開挖至第五層土之前持續增大，開挖完第五層土后略有降低，且降低的幅度相對較小。圍護結構在開挖完成第四層土體時，承受到最大的彎矩832.241 kN·m，表明：施工到第五層土時結構所受到的彎矩最大。從每層土開挖時圍護結構彎矩的增量來看，第三層土開挖完成時彎矩增量最大，為113.666 kN·m，這表明在第四層土開挖的過程中，結構的內力迅速增長。因此，為了保證圍護結構的穩定性與施工安全，考慮對第四道支撐進行加強處理，如：增大其截面尺寸、增加其他預施加荷載等措施。

從圖10中可以看出，對應開挖工況1～5，圍護結構產生的最大水平位移分別為8.29、8.95、9.75、10、10 mm，最大水平位移發生位置分別在開挖深度8、9、11、11、11 m。隨著基坑開挖的進行，圍護結構的水平位移持續增大，位移最終收斂于10 mm左右。最大水平位移主要發生在11 m左右的位置，正好處于素填土的地層中，這表明該層素填土的土層性質相對較差。因此，在設計施工中，應采取相應措施進行加固。

圖10 各工況下圍護結構水平位移模擬結果

圖11 理論分析與數值模擬最大水平位移對比

根據圖8、圖9及圖11的對比情況，綜合分析出基坑圍護結構的內力與變形的變化趨勢與數值模擬基本一致，在中間土體的開挖工況下具有良好的符合性，但在第一、二、四、五層土體的開挖工況下，呈現出一定的差異性。一方面，兩者彎矩的最大偏差約為12%，剪力的最大偏差約為5%，水平位移的最大偏差約為20%，但誤差均在允許范圍之內，這表明預測結果與數值模擬結果相符合。另一方面，相對于數值模型，理論分析模型關于內力的計算結果偏差較小，而關于計算咬合樁在施工過程中的最大水平位移值高出一定范圍，這表明采用基于不同樁徑葷素樁組成的鉆孔咬合樁結構，推導出修正后的等效剛度計算模型的結果是偏安全、保守的。當涉及到具體的鉆孔咬合樁的設計時，可根據實際選擇適應的計算方法。

3 結論與建議

基于剛度等效原則，由不同樁徑葷素樁組成的鉆孔咬合樁結構，推導出修正后的等效剛度的計算模型。采用Midas/GTS NX有限元軟件，模擬分析了拋填石層超深基坑施工過程中圍護結構的內力變化與位移演變規律，并與理論計算結果進行對比分析，得到結論為：

1) 從理論計算模型分析來看，鋼筋混凝土樁與素樁承受彎矩之比約為7:3，設計過程若能合理利用素樁的承載受力性能，有利于發揮工程經濟性，降低施工成本。

2) 隨著基坑開挖的進行，圍護結構的水平位移持續增大，最大水平位移保持10～13 mm的范圍，一般發生在埋深11 m處素填土的地層中。因此，在設計施工中，應對該層土進行土體加固，可以采取預注漿加固等施工措施。

3) 在基坑開挖過程中，開挖工況下圍護結構的內力與變形增長較快，及時增加內支撐，有利于降低圍護結構水平位移，防止內力變形過大，出現不穩定性。

4) 從理論模型分析結果來看，第三層土開挖完成后，結構內力迅速增長，有必要對這一階段的內支撐進行加強處理，可以適當增大第四道內支撐截面尺寸，增加預施加軸力，改善內撐混凝土強度。

5) 本研究理論模型結果是偏安全、保守的，但分析圍護結構的內力與變形的變化趨勢與數值模擬一致，相互印證后較符合，能夠有效預測鉆孔咬合樁內力與位移的發展規律，進而指導圍護結構的合理設計與安全施工。

[1] Cervia A R D. Construction of the deep cut-off at the walter F. George dam[C]//GeoSupport Conference 2004. Orlando, Florida, USA. Reston, VA, USA: American Society of Civil Engineers, 2004: 530?544.

[2] Anderson T C. Secant piles support access shafts for tunnel crossing in difficult geologic conditions[C]// GeoSupport Conference 2004. Orlando, Florida, USA. Reston, VA, USA: American Society of Civil Engineers, 2004: 299?308.

[3] Finno R J, Bryson L S. Response of building adjacent to stiff excavation support system in soft clay[J]. Journal of Performance of Constructed Facilities,2002,16(1):10?20.

[4] 陳海兵,梁發云,何招智.咬合樁在鄰近高填土基坑中的工程應用與實測分析[J].土木建筑與環境工程,2014, 36(3): 1?5. (CHEN Hai-bing, LIANG Fa-yun, HE Zhao-zhi. Application of secant piles in excavation pit adjacent to the high-filling surcharge loads and in situ measurements[J]. Journal of Civil, Architectural & Environmental Engineering,2014,36(3):1?5.(in Chinese))

[5] 胡琦,陳彧,柯瀚,等.深基坑工程中的咬合樁受力變形分析[J].巖土力學,2008,29(8):2144?2148,2176.(HU Qi, CHEN Yu, KE Han, et al. Analysis of bearing and deforming behavior of secant piles of deep foundation pit excavation[J]. Rock and Soil Mechanics, 2008, 29(8): 2144?2148,2176.(in Chinese))

[6] 楊虹衛,楊新偉.鉆孔咬合樁的配筋計算方法[J].地下空間與工程學報, 2008, 4(3): 402?405, 414. (YANG Hong-wei, YANG Xin-wei. Computational method of reinforcement of the borehole cast-in-place pile[J]. Chinese Journal of Underground Space and Engineering, 2008,4(3):402?405,414.(in Chinese))

[7] 廖少明,周學領,宋博,等.咬合樁支護結構的抗彎承載特性研究[J].巖土工程學報,2008,30(1):72?78.(LIAO Shao-ming, ZHOU Xue-ling, SONG Bo, et al. Study on the bending performance of secant pile retaining wall[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(1):72?78.(in Chinese))

[8] 任紅林,吳建忠,朱正誼.鉆孔咬合樁基坑圍護結構設計計算方法探討[J].城市軌道交通研究,2010,13(9):29?32, 37.(REN Hong-lin, WU Jian-zhong, ZHU Zheng-yi. Calculation method of secant pile as excavation retaining structure[J].Urban Mass Transit,2010,13(9):29?32,37.(in Chinese))

[9] 廖中林,徐傳鵬.地鐵車站基坑鉆孔咬合樁施工與設計[J].交通科技與經濟,2013,15(2):24?27.(LIAO Zhong-lin, XU Chuan-peng. Design and construction of borehole cast-in-place piles in the subway station pit[J]. Technology & Economy in Areas of Communications, 2013,15(2):24?27.(in Chinese))

[10] 張前,楊霜,楊振甲.咬合樁等效剛度法在基坑工程中應用[J].低溫建筑技術,2015,37(8):140?141.(ZHANG Qian, YANG Shuang, YANG Zhen-jia. The application of the equivalent stiffness method of occlusal piles on foundation pit engineering[J]. Low Temperature Architecture Technology, 2015, 37(8): 140?141. (in Chinese))

[11] 張楨.蘭州地鐵鉆孔咬合樁設計優化研究[D].蘭州:蘭州大學,2016.(ZHANG Zhen.Design optimization of the secant pile wall at Lanzhou subway[D].Lanzhou: Lanzhou University,2016.(in Chinese))

[12] 閆文進.折線型咬合樁圍護結構承載特性分析及應用研究[D].蘭州:蘭州理工大學,2018.(YAN Wen-jin. Analysis for the bearing characteristics and application of folded linear secant pile[D].Lanzhou:Lanzhou University of Technology,2018.(in Chinese))

[13] 張健,周中,邢凱.盾構隧道近接施工對地層變形的影響[J].交通科學與工程,2019,35(4):63?70.(ZHANG Jian, ZHOU Zhong, XING Kai. The influence of shield tunnel construction adjacent to the existing structure on deep soil displacement[J].Journal of Transport Science and Engineering,2019,35(4):63?70.(in Chinese))

[14] 周廣友,李聰,胡勇,等.箱涵頂進施工過程中路面沉降的數值模擬分析[J].公路工程,2020,45(3):29?37.(ZHOU Guang-you, LI Cong, HU Yong, et al. Numerical simulation analysis of pavement settlement during the construction process of box culvert-jacking[J]. Highway Engineering,2020,45(3):29?37.(in Chinese))

Research on the force calculation model of the bored occluding pile composed of unequal diameter reinforced concrete piles and plain piles

HONG Yuan1, PENG Zheng-kun1, REN Chen-hao1, ZENG Xi-fan1, MEI Dong1, LI Fan2

(1. The Eighth Sinohydro Bureau Co., Ltd., Changsha 410004, China;2. School of Civil Engineering, Central South University, Changsha 410075, China)

Aiming at the shortcomings in the current theoretical calculation model of occluding piles equivalent to underground continuous walls, this paper considered the influence of the relative size of the pile diameter between the reinforced concrete piles and the plain concrete piles in the bored occluding piles and the stiffness reduction effect of the plain concrete piles. The corrected equivalent rigidity occluding pile force and deformation theoretical calculation model was derived. At the same time, relying on the occluding pile project of the main enclosure structure of the Taiziwan Station of Shenzhen Urban Rail Transit Line 12, the finite element software (Midas/GTS NX) was used to establish a numerical simulation test to analyze the force mechanism and the deformation law of the enclosure structure during the whole process of excavation of the foundation pit. The deformation law verifies the accuracy of the theoretical calculation model. The analysis results of the theoretical model and the numerical model were combined to point out the structural weaknesses in the bored occluding pile engineering example, and suitable measures were given. The research results show that with the excavation of the foundation pit, the maximum horizontal displacement of the enclosure structure finally converges to about 10 mm. Under the excavation conditions, the internal force and deformation of the enclosure structure increase rapidly, especially after the third layer of soil is excavated. The theoretical model in this paper can effectively predict the internal force and displacement development law of the bored occluding pile during the construction of the foundation pit, and then guide the rational design and safe construction of the enclosure structure.

bored occluding pile; equivalent stiffness method; numerical simulation; stress mechanism; deformation law

TU753.3

A

1674 ? 599X(2021)01 ? 0075 ? 07

2020?08?22

國家自然科學基金資助項目(50908234)；湖南省自然科學基金項目(2020JJ4743)

洪源(1991?)，男，中國水利水電第八工程局有限公司工程師。