趙明華 馬思齊 肖堯 楊超煒

摘? ?要：為解決現(xiàn)有方法難以考慮加筋材料實際受力狀態(tài)而低估筋箍碎石樁復合地基沉降的問題，假設樁土等應變且均為線彈性材料，選取單樁有效加固單元整體作為分析對象，其中碎石樁同時受到加筋材料和土環(huán)的約束作用，而土環(huán)則可以考慮為同時受到單樁有效加固范圍外土體的靜止土壓力和內部碎石樁鼓脹壓力共同作用的厚壁圓筒，再結合廣義胡克定律得到應力應變關系，進而導出了筋箍碎石樁復合地基沉降計算的新方法.采用工程實例驗證并與已有方法進行比較，同時分析了外荷載水平、置換率對復合地基總沉降和樁土應力比的影響.與已有方法相比，該方法可以使上部荷載和筋箍碎石樁側向受力變形聯(lián)動進而調整加筋材料內力，從而更符合筋箍碎石樁實際受力變形情況，計算值與實測值相對誤差為+5.70%，與現(xiàn)有保守計算方法相比誤差最小.參數(shù)分析表明：置換率對于控制復合地基總沉降具有重要作用，外荷載一定時復合地基總沉降與樁土應力比均隨著置換率提高而減小.

關鍵詞：復合地基;筋箍碎石樁;徑向位移;沉降

中圖分類號：TU443? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻標志碼：A

文章編號：1674—2974（2019）01—0101—08

Abstract：The existing method is difficult to consider the actual state of geogrid ribs under loading， which could underestimate the settlement of the geogrid-encased stone column （GESC） composite foundation. In order to solve this problem， in this paper， the unit cell concept was employed as the objects in the analysis， while it was assumed that GESC together with soil was elastic materials and evenly deformed. Stone column was constrained by geogrid and earth ring at the same time， while the earth ring can be considered as a ring with earth pressure at rest provided by soil outside and bulging pressure caused by the stone column installed inside. Hook's law was then introduced to gain the stress-strain formula of the unit cell. Then， a different method of settlement calculation of GESC composite foundation was established. The rationality of the results of this method was verified by an engineering example and also compared with the existing method. The result demonstrates that the method agrees well with the engineering example. Through the comparison with the existing method， it is found that the relative error between the calculated value and measured value is +5.70%， which demonstrates the method presented in this paper is safer and more accurate as a result of the combined top loading and lateral deformation of GESC. The parametric study shows that the replacement ratio has a significant influence on the controlling settlement of composite foundation， and the increase of the replacement ratio can effectively reduce the pile - soil ratio and settlement of composite foundation when a certain level load is applied.

Key words：composite foundation;geogrid-encased column;radial displacement;settlement

作為取材便利、造價經濟且工程性能良好的軟土地基處理方式，碎石樁在道路工程上得到了大量應用.然而碎石樁復合地基存在易發(fā)生鼓脹破壞的缺陷，使得承載力提升和沉降控制均受到原有軟土基礎的極大制約.在該背景下，筋箍碎石樁應運而生，有效解決了碎石樁的鼓脹破壞問題[1-2]，因此筋箍碎石樁復合地基的承載力和沉降計算方法也成為關注重點.其中，找到有效的沉降計算方法，是使筋箍碎石樁復合地基更加適用于工程實際的關鍵.

目前對于普通碎石樁復合地基的沉降計算，得到廣泛應用的有等效模量法、應力修正法和樁身壓縮量法[3-4].然而筋箍碎石樁的筋箍段不再完全呈現(xiàn)散體材料樁的性質，相較而言其受力機理更為復雜，散體材料樁的沉降計算方法不再適用，一些學者便開始尋求適用于計算筋箍碎石樁沉降的方法，在試驗、數(shù)值模擬和理論計算方面均做出了研究和探討.在試驗方面，Afshar等[5]對碎石樁分別做了豎向和側向加筋試驗，發(fā)現(xiàn)增加筋材量總能提高復合地基承載力，并且在所有的試驗組中，筋箍碎石樁的鼓脹破壞是主要破壞形式，起著決定承載力大小的作用;Gu等[6]依托大比例室內試驗，研究了荷載增長對筋箍碎石樁樁土應力比的影響，發(fā)現(xiàn)筋箍套筒的存在能有效減少碎石樁體鼓脹破壞的產生，從而達到控制復合地基沉降的目的;Almeida等[7]進行足尺試驗，結果表明箍筋的存在對于減小碎石樁徑向變形進而減小復合地基沉降有很大貢獻;Hong等[8]通過縮尺砂樁試驗，發(fā)現(xiàn)筋箍套筒的存在不僅能夠明顯減少徑向應變并且加筋材料在破壞后仍能對復合地基承載力提供一定貢獻.

數(shù)值分析最大程度上還原了工程實際中筋箍碎石樁復合地基的受力狀態(tài)和變形過程，是計算筋箍碎石樁復合地基沉降的有效方法[9-10].在數(shù)值模擬方面，Yoo[11]通過數(shù)值建模分析認為碎石樁的有效加固深度只存在于加固區(qū)軟土深度大于10 m的時候，而在其他情況下，為了使沉降量最小，有必要對碎石樁進行通長加筋;Khabbazian等[12]應用有效加固單元概念，建立了3種不同數(shù)值分析模型，發(fā)現(xiàn)筋箍碎石樁三維全樁模型能更準確地反應其受力過程中土工格柵內部產生的拉應力，同時結果顯示格柵縱橫兩方向的拉力并不相等;Hosseinpour等[13]基于有效加固單元建立了二維有限元分析模型，分析認為能提供徑向約束的土工合成材料對碎石樁的承載力提高和地基沉降量的減少的作用更明顯，并且徑向約束之間存在最優(yōu)間距.數(shù)值模擬中的模型建立和參數(shù)選取極為復雜，需要扎實的理論基礎和實際經驗，且計算過程費時費力，不具備易得性和普適性，而理論計算法則更具可操作性.在理論計算方面，曹文貴等[14-16]將筋箍碎石樁劃分為筋箍段、非筋箍段和下臥層三部分，在筋箍段考慮樁土間的相對滑移，在非筋箍段考慮樁土的鼓脹變形，分別建立了不同的沉降計算式，取得了較好的成果;另外，曹文貴等[17]利用孔隙介質分析方法，建立了散體材料復合地基的孔隙介質力學模型，并引入分級加載思想計算復合地基總沉降，具有很好的精度;Zhang等[18]基于彈性理論和有效加固單元概念，忽略圓周方向接觸面上的摩擦力，并將樁周土提供的側向限制近似取為被動土壓力而推導出了一種計算筋箍碎石樁復合地基沉降變形的計算方法;張玲等[19]基于荷載傳遞，并綜合考慮了散體材料豎向及側向變形的特性，得到了計算散體材料樁復合地基沉降的一種新方法.

本文對筋箍碎石樁沉降進行更進一步探討，旨在理論計算方面提供一些新思路.鑒于Alamgir等[20]創(chuàng)新性地提出了“有效加固單元”這一概念，本文即以此為依托，主要分析了筋箍碎石樁復合地基的沉降.基于合理的前提假設，主要運用厚壁圓筒模型考慮單元體中裹敷筋箍碎石樁的土體受力狀態(tài)，從而進一步分析筋材的實際受力狀態(tài)，在考慮樁、土實際分擔的荷載后，提出了基于厚壁圓筒模型的筋箍碎石樁沉降計算新方法.

1? ?計算模型與基本假定

筋箍碎石樁的沉降由樁體加固區(qū)沉降和下臥層沉降兩部分組成：

下臥層沉降目前已有較好的計算方法，關鍵在于樁體加固區(qū)的沉降[21-22]，其中加固區(qū)的沉降又以樁體的沉降計算最為復雜，下面將重點對加固區(qū)樁體的沉降進行分析和探討.

1.1? ?基本假定

碎石樁復合地基多以群樁形式出現(xiàn)，樁、土和墊層形成一個整體來承受上部荷載，在樁、土的相互作用下，其受力情況較為復雜.為簡化分析，現(xiàn)取筋箍碎石樁復合地基中的單根樁加固區(qū)域作為分析單元，如圖1所示.

圖1中：rp為樁體半徑;re為樁體影響半徑（即分析單元半徑）;sd為樁中心距.

根據布樁方式的不同，分析單元半徑可以由單樁加固范圍直徑公式得到：

對于復合地基而言，剛性樁和柔性樁在樁土界面上均會產生相對位移.在理論分析中，由于樁土界面處的互鎖作用，可假設樁土在界面處完全粘結[23]，且目前在實驗與數(shù)值模型分析的文獻中，也多見剛性板加載測試[24-25]，忽略樁土界面相對滑動也可以達到較好的精度.故為了簡化計算，本文計算復合地基最終沉降過程中亦忽略樁側摩阻力的影響.在計算復合地基沉降時，一般認為樁體所受荷載在承載范圍內，未達到極限狀態(tài)，可在計算中只考慮樁土的彈性變形.為簡化計算，根據圖1所示分析單元，作如下假定：

1）單元體之間相互不影響，即忽略群樁效應對單元體受力特性的影響;

2）筋箍碎石樁與樁間土變形相同，樁土界面不產生摩阻力，即樁土等應變;

3）碎石樁體和樁間土為彈性體，沉降過程中只產生彈性變形.

1.2? ?樁單元應力、應變分析

在樁土彈性變形階段，已知單元體的應力狀態(tài)，便可求解出單元體應變.因此在求解樁體沉降時，可從樁體的受力情況入手，將樁體沿深度方向等分為N段微單元，求解出每段微單元的沉降后再進行疊加，便可得到樁體總沉降.

首先將樁體等分成N段：

如圖2所示，取第i段樁身微單元進行分析，并用微單元上表面的應力來計算其變形.

2? ?基于厚壁圓筒模型的沉降計算方法

通過對微單元的應力應變分析可知，求解樁單元的沉降關鍵在于求解樁體的應力狀態(tài).接下來即從樁單元軸力和所受圍限力兩方面進行分析.

2.1? ?樁單元軸力

根據假定2），在不考慮樁土界面摩阻力的情況下，并忽略樁身荷載傳遞損失部分，樁身微單元上表面的豎向應力只由樁頂荷載和樁體自重兩部分組成，其表達式為：

2.2? ?樁單元圍限力

對于筋箍碎石樁加筋部分，樁體所受圍限力由土體和筋箍套筒兩方面提供，即

對于非加箍段，令式（6）中筋箍套筒提供的圍限力σgr，i為零，則樁體所受圍限力只由土體σsr，i提供，即

接下來分別對土體圍限力和筋材圍限力進行求解.

2.2.1? ?土體圍限力

如圖2所示樁體微單元，樁周土體受到樁體鼓脹產生的壓力時，可將其受力變形過程視為厚壁圓筒模型內壁受力變形過程，土體在樁土界面處受到樁體擴張壓力作用.取微單元上表面所在平面進行分析，根據假定3），當土體尚處于彈性變形階段時，利用彈性力學的厚壁圓筒公式[26]來描述其位移場為：

式中：σr0，i為半徑re處土體徑向應力;ur，i 為任意半徑r處位移;vs為樁周土泊松比;Es為樁周土彈性模量.

在該平面上，可合理假設單根樁影響半徑re處的側向土壓力為靜止土壓力，土環(huán)內部壓力與土體對碎石樁圍限力數(shù)值上相等，如圖3所示.

式中：σsz，i為土體微單元上表面豎向應力;K0，s 為樁間土體靜止土壓力系數(shù)，完全彈性條件下，K0，s = vs /（1- vs）;γs為樁周土體重度;Ps為樁間土體上部荷載;n為樁土應力比;Pp為樁體上部荷載.

同時，對于樁土界面處有：

令式（8）中變量r為rp，并把式（10）（12）代入式（8），為使公式更加簡化，令t = re2/rp2，即可解得土體所提供的圍限力為：

2.2.2? ?筋材圍限力

當鼓脹變形發(fā)生時，筋材與樁、土發(fā)生協(xié)調變形，其受力示意圖如圖4所示.

■

從圖4中可知，將樁周土體圍限力σsr，i和碎石樁鼓脹壓力σpr，i沿著Y軸積分，引入筋材拉力Tg，可由受力平衡得到：

根據以往研究成果，當樁體處于彈性變形階段時，由胡克定律可知，樁身微單元的應力應變關

系為：

式中：σpz，i 為微單元豎向應力;σpr，i為微單元徑向應力;Ep為碎石樁彈性模量;vp為碎石樁泊松比;εpz，i為微單元豎向應變;εpr，i為微單元徑向應變.

因此，已知樁身微單元的豎向應力σpz，i和徑向應力σpr，i，就可計算出對應的豎向應變εpz，i，即樁體微單元沉降，從而求解出樁體的總沉降.

根據式（17），微單元徑向應變表達式為：

聯(lián)立式（6）（15）和式（16），便可算出筋材圍限力的表達式為：

土體圍限力σsr，i和微單元上表面豎向應力σpz，i可分別由式（13）和式（5）計算得到，再將筋材圍限力和土體圍限力代入式（6），便得到樁單元的徑向圍限力，則樁體微單元徑向應力得解.

2.3? ?樁、土沉降計算

根據上述分析，碎石樁與樁間土豎向變形可分別通過下式計算：

式中：εpz，i和εsz，i分別為第i段碎石樁與第i段樁間土的豎向應變.

將式（12）（18）代入式（6），得到σpr，i后再與式（5）一并代入式（20），即可得到第i段加筋段碎石樁的豎向應變.而對于非筋箍段，只需令式（18）中筋材抗拉剛度Jg的值等于零，所得結果即為第i段非筋箍段樁體微單元的豎向應變.將各段應變求和即得到樁體沉降.

由于本文將復合地基分為N段進行求解，故可利用編程語言對問題進行計算，思路如下：

給定樁土應力比范圍與計算步長，利用循環(huán)嵌套語句，根據式 （20）分別計算樁體壓縮量ΔSp和土體壓縮量ΔSs，若兩壓縮量滿足ΔSp - ΔSs≤ δ（δ為給定誤差限），則認為在該給定樁土應力比情況下，滿足本文的樁土等應變假設，取此時的樁土應力比和復合地基沉降量為符合要求的計算值.再根據《建筑地基基礎設計規(guī)范》[22]方法可得到下臥層土體的沉降量，從而計算出整個復合地基的沉降.在實際工程中，采用散體材料加固處理軟土地基時，一般會設置褥墊層調整樁土應力比.該方法是將褥墊層在復合地基中的作用簡化為樁土應力比n，在計算復合地基沉降時以樁土等沉為收斂條件，通過循環(huán)迭代，可以得到合理的樁土應力比.同時該方法簡便快捷，避免了沉降公式過于復雜冗長的問題，契合本文計算思路，并且可以較快達到精度要求.計算流程如圖5所示.

3? ?工程實例驗證與分析

3.1? ?算例驗證

四川某繞城高速公路互通式立交橋A匝道，采用土工格柵加筋碎石樁-砂墊層對路基進行加固處理[27]，樁徑d為0.6 m，樁長L為6 m，筋箍段長Lg為2 m，樁中心距sd為2.0 m，正方形布樁;碎石樁體壓縮模量Esp為100 MPa，泊松比vp為0.4，重度γp為20 kN/m3;筋材抗拉剛度Jg為120 kN/m，總外荷載P0為132 kPa，置換率8.16%.根據樁土應力比、置換率和樁體荷載Pp、土體荷載Ps之間的關系可知樁體荷載Pp、土體荷載Ps分別為350.84 kPa和112.45 kPa[10];路基實測沉降ws為110 mm[16];土層分布情況如圖6所示.土層各參數(shù)信息如表1所示.

將表1中各參數(shù)值代入前文所示的沉降步驟計算公式中進行計算，當誤差限δ取0.1 mm時，得到復合地基沉降計算值為116.3 mm，從工程實例中測得的沉降實測值為110 mm，二者的相對誤差為+5.70%，可見本文沉降計算式所得結果與實測值吻合較好，能夠滿足工程要求，證明本文方法能夠在一定程度上正確反映工程中的實際情況.

將以往文獻對該實例的沉降分析結果與本文方法所得結果進行比較，如表2所示.

一般來說，理論沉降計算結果通常為最終沉降值，由于實測值通常也是在一段時間內的觀測值，而非最終沉降值，計算結果較實測值偏大是比較合理的，由表2可知，本文計算方法在偏于保守的結果中相對誤差最小，可以從側面證明本方法的合理性.

3.2? ?影響因素分析

筋箍碎石樁復合地基主要參數(shù)有土體性質參數(shù)、筋材性質參數(shù)、樁徑以及面積置換率等.這些因素都在一定程度上影響著復合地基受力變形性質.以工程實例參數(shù)為基礎，利用本文方法進一步探討置換率對筋箍碎石樁復合地基的樁土應力比和總沉降的影響.分別改變置換率及復合地基頂面總外荷載，其他非討論參數(shù)因素不變，得到筋箍碎石樁復合地基樁土應力比以及復合地基沉降量的變化曲線圖，計算結果分別如圖7和圖8所示.

3.2.1? ?置換率對樁土應力比的影響

圖7為筋箍碎石樁復合地基樁土應力比隨置換率變化的變化趨勢.可以看出，置換率一定時，隨著外荷載增大樁土應力比減小，各級荷載下的減小程度相近，故增大置換率可以使樁土應力分布更加均勻，從而可以降低樁土應力比.相比較而言，樁土應力比在外荷載水平較低時偏大.

3.2.2? ?置換率對總沉降的影響

圖8為筋箍碎石樁復合地基沉降隨著復合地基置換率增加時兩者的變化關系曲線.可以看出復合地基總沉降量隨著置換率增加而減小，隨著荷載增加而增加.而且在外荷載較大時，曲線斜率較大，故此時若提高置換率可以比較明顯地減小復合地基的沉降.

4? ?結? ?論

本文在深入研究筋箍碎石樁復合地基受力變形的基礎上，以單樁有效加固范圍內的復合地基作為研究對象，結合現(xiàn)有學者在現(xiàn)場試驗和理論研究等各方面的成果，得到了筋箍碎石樁復合地基的沉降計算新方法，并利用工程實例對該模型進行驗證以及相應的參數(shù)分析.

1）基于厚壁圓筒理論，將單樁有效加固范圍內的土體考慮為厚壁圓筒，充分考慮了樁、土、筋材三者之間的變形協(xié)調;

2）從加筋段實際受力狀態(tài)出發(fā)，考慮筋材在任意應力條件下的應力，進而更加準確地反映筋箍碎石樁復合地基的沉降;

3）參數(shù)分析表明：提高置換率可以優(yōu)化樁-土應力分布，減小樁土應力比，從而減小復合地基總沉降;

4）本文筋箍碎石樁沉降計算結果與工程實際測量結果吻合良好，證明了本文沉降計算方法能較準確反映筋箍碎石樁復合地基的受力變形情況.

參考文獻

[1]? ? 唐承鐵，劉猛. 加筋碎石樁格柵強度設計初探[J].公路工程，2012，37（6）：25—28.

TANG C T，LIU M. A tentative exploration on the strength design of the geogrid in geosynthetic-encased stone columns[J]. Highway Engineering，2012，37（6）：25—28.（In Chinese）

[2]? ? 陳慶. 格柵套筒加筋碎石樁復合地基承載力計算方法研究[D].長沙：湖南大學土木工程學院，2011：11—12.

CHEN Q. The research on the calculation method for eosynthetic-encased gravel pile composite foundation[D]. Changsha： College of Civil Engineering，Hunan University，2011：11—12.（In Chinese）

[3]? ? 龔曉南. 復合地基理論及工程應用[M].北京：中國建筑工業(yè)出版社，2007：130—131.

GONG X N. Theory and engineering application of composite foundation[M]. Beijing：China Architecture and Building Press，2007：130—131.（ In Chinese）

[4]? ? OMINE K，OHNO S. Deformation analysis of composite ground by homogenization method[C]// Proceedings of the Fourteenth?International Conference on Soil Mechanics and Foundation Engineering.Rotterdam： Balkema A A，1997：719—722.

[5]? AFSHAR J N，GHAZAVI M. Experimental studies on bearing capacity of geosynthetic reinforced stone columns[J].Arabian Journal for Science & Engineering，2014，39（3）：1559—1571.

[6]? ? GU M，ZHAO M，ZHANG L，et al. Effects of geogrid encasement on lateral and vertical deformations of stone columns in model tests[J].Geosynthetics International，2016，23（2）：100—112.

[7]? ? ALMEIDA M S S，HOSSEINPOUR I，RICCIO M，et al. Behavior of geotextile-encased granular columns supporting test embankment on soft deposit[J]. Journal of Geotechnical & Geoenvironmental?Engineering，2015，141（3）：04014116.

[8]? ? HONG Y S，WU C S，YU Y S. Model tests on geotextile-encased granular columns under 1-g and undrained conditions[J].Geotextiles & Geomembranes，2016，44 （1）：13—27.

[9]? ? 李小青，周偉.考慮樁土接觸特性的復合地基承載機理的數(shù)值分析研究[J]. 巖土工程學報，2006，28（4）： 529—532.

LI X Q，ZHOU W. Numerical analysis on bearing capacity of composite foundation in consideration of contacting characteristics of pile-soil interface[J]. Chinese Journal of Geotechnical Engineering，2006，28（4）： 529—532.（In Chinese）

[10]? 崔，張志耕，閆澍旺. 碎石樁聯(lián)合土工格柵復合地基處理濕地軟基的機制研究[J]. 巖土力學，2009，30（6）： 1764—1768.

CUI W，ZHANG Z G，YAN S W. Mechanism of composite ground with stone columns and geogrid in wetland[J]. Rock and Soil Mechanics，2009，30（6）：1764—1768.（In Chinese）

[11] YOO C. Settlement behavior of embankment on geosynthetic-encased stone column installed soft ground-A numerical investigation[J]. Geotextiles & Geomembranes，2015，43（6）：484—492.

[12] KHABBAZIAN M，KALIAKIN V N，MEEHAN C L. Column supported embankments with geosynthetic encased columns： validity of the unit cell concept[J]. Geotechnical and Geological Engineering，2015，33（3）：425—442.

[13] HOSSEINPOUR I，RICCIO M，ALMEIDA M S S. Numerical evaluation of a granular column reinforced by geosynthetics using encasement and laminated disks[J].Geotextiles & Geomembranes，2014，42（4）：363—373.

[14]? 曹文貴，王江營，羅忠濤，等.考慮樁-土滑移的柔性基礎下復合地基沉降分析方法[J].中國公路學報，2012，25（6）： 9—16.

CAO W G，WANG J Y，LUO Z T，et al. Settlement calculation method for composite ground considering pile-soil slip under flexible foundation[J]. China Journal of Highway and Transport，2012，25（6）： 9—16.（In Chinese）

[15]? 曹文貴，楊澤華.柔性基礎下碎石樁加筋復合地基沉降分析方法[J].巖土工程學報，2012，34（11）：1997—2004.

CAO W G，YANG Z H. Method for analyzing settlement of reinforced-gravel-pile composite ground under flexible foundation[J]. Chinese Journal of Geotechnical Engineering，2012，34（11）： 1997—2004.（In Chinese）

[16]? 曹文貴，趙聚才，賀敏，等.柔性基礎下筋箍碎石樁復合地基變形機理及其沉降分析方法[J].巖土工程學報，2014，36（5）： 818—826.

CAO W G，ZHAO J C，HE M，et al. Deformation mechanism and?settlement analysis method of reinforced-hoop-gravel-pile? composite ground under flexible foundation[J]. Chinese Journal of Geotechnical Engineering，2014，36（5）：818—826.（In Chinese）

[17]? 曹文貴，唐旖旎，王江營.基于孔隙介質模型的散體材料樁復合地基沉降分層總和分析方法[J].湖南大學學報（自然科學版），2014，41（11）： 87—95.

CAO W G，TANU Y N，WAND J Y. A layer-wise summation analysis method for the settlement of composite ground with friable material piles based on porous medium model[J].Journal of Hunan University（Natural Sciences），2014，41（11）： 87—95.（In Chinese）

[18]? ZHANG L，ZHAO M. Deformation analysis of geotextile-encased stone columns[J].International Journal of Geomechanics，2015，15（3）：04014053.

[19]? 張玲，趙明華.考慮鼓脹變形的散體材料樁復合地基沉降計算[J].湖南大學學報（自然科學版），2016，43（5）：120—124.

ZHANG L，ZHAO M H. Settlement calculation of granular material pile composite foundation with consideration of lateral bulging[J].Journal of Hunan University（Natural Sciences），2016，43（5）：120—124.（In Chinese）

[20] ALAMGIR M，MIURA N，POOROOSHASB H B. Deformation analysis of soft ground reinforced by columnar inclusions[J].Computers and Geotechnics，1996，18（4）：267—290.

[21]? 章定文，劉松玉.路堤荷載下柔性樁復合地基沉降實用計算方法[J].巖土力學，2007，28（6）：1133—1138.

ZHANG D W，LIU S Y. Practical method for settlement calculation of flexible columns composite foundation under embankment[J].Rock and Soil Mechanics，2007，28（6）：1133—1138.（In Chinese）

[22]? GB 50007—2011? ?建筑地基基礎設計規(guī)范[S].北京：中國建筑工業(yè)出版，2011：28—30.

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

[23] CASTRO J. An analytical solution for the settlement of stone columns beneath rigid footings[J].Acta Geotechnica，2014，11（2）：1—16.

[24]? DEBNATH P，DEY A K. Bearing capacity of geogrid reinforced sand over encased stone columns in soft clay[J]. Geotextiles & Geomembranes，2017，45（6）：653—664.

[25]? CASTRO J.Groups of encased stone columns： Influence of column length and arrangement[J]. Geotextiles & Geomembranes，2017， 45（2）：68—80.

[26]? 徐秉業(yè)，劉信聲.應用彈塑性力學[M].北京： 清華大學出版社，1995：188—191.

XU B Y，LIU X S. Applied elastic and plastic mechanics[M]. Beijing：Tsinghua University Press，1995：188—191.（In Chinese）

[27]? 黃兵.GSPS軟基加固系統(tǒng)機理及其在成都繞城高速公路的應用研究[D].成都：西南交通大學土木工程學院，2002：41—51.

HUANG B. The mechanism of GSPS soft ground treatment system and its application in Chengdu City Highway[D].Chengdu：School of Civil Engineering，Southwest Jiaotong University，2002：41—51.（In Chinese）

[28]? 張贊威.筋箍碎石樁復合地基承載機理及其變形分析[D].長沙：湖南大學土木工程學院，2014：52—53

ZHANG Z W. The analysis of geotextile encased stone columns composite foundation bear deformation mechanism[D].Changsha： College of Civil Engineering，Hunan University，2014：52—53.（In Chinese）