穿越無充填溶洞時單樁承載特性研究

張乾青 喬勝石 邢宇鋮 張凱 崔偉 王志遠(yuǎn) 陳迪楊

摘要：為解決穿越溶洞型單樁的設(shè)計難題，研究了穿越無充填溶洞時單樁的承載特性. 基于現(xiàn)場物探和室內(nèi)試驗(yàn)結(jié)果明確了溶洞的發(fā)育情況和巖土體的物理力學(xué)參數(shù)，采用有限元軟件建立了14組數(shù)值計算模型，分別對不同溶洞半徑、不同溶洞位置及不同樁徑的工況進(jìn)行了模擬，確定了BoxLucas1荷載傳遞模型的適用性，提出了穿越無充填溶洞單樁的承載特性計算方法，研究結(jié)果表明，樁徑增加和溶洞充填時溶洞對單樁承載力的影響將得到削弱，溶洞的半徑和位置對樁基承載力具有顯著影響.實(shí)際工程中應(yīng)重視前期物探對溶洞發(fā)育情況的準(zhǔn)確探查，降低溶洞對樁基承載力的影響.相關(guān)研究成果可為巖溶區(qū)單樁參數(shù)設(shè)計和承載力預(yù)測提供理論依據(jù).

關(guān)鍵詞：巖溶區(qū)；物探方法；穿越溶洞型單樁；承載特性；數(shù)值模擬

中圖分類號：TU473.1文獻(xiàn)標(biāo)志碼：A

Study on Bearing Characteristics of Single Pile Crossing Unfilled Karst Cave

ZHANG Qianqing1，2，QIAO Shengshi1，2，XING Yucheng1，2，ZHANG Kai3，CUI Wei4，WANG Zhiyuan3，CHEN Diyang2

（1. Geotechnical and Structural Engineering Research Center，Shandong University，Jinan 250061，China；2. Research Institute of New Material and Intelligent Equipment，Shandong University，Dezhou 251114，China；3. Shandong Luqiao Group Co. Ltd.，Jinan 250021，China；4. School of Qilu Transportation，Shandong University，Jinan 250002，China）

Abstract：To solve the design problem of single pile crossing karst caves in practice，the bearing capacity of a single pile crossing unfilled karst cave was analyzed in this paper. Based on the field geophysical and laboratory test results，the distribution of the karst cave and the physical and mechanical parameters of rock and soil were determined. Fourteen groups of numerical calculation models were established by using finite element software to capturethe bearing characteristics of a single pile crossing an unfilled karst cave under different conditions of cave radius，cave position and pile diameter. The applicability of the BoxLucas1 load transfer model was verified by using numerical results. A calculation method was then proposed to predict the response of a single pile crossing unfilled karst caves. The present study showed that the influence of karst caves on the bearing capacity of a single pile was weakened by increasing pile diameter and filling karst caves，and the radius and position of karst caves had a significant influence on the single pile bearing capacity. For practical purposes，the accurate exploration of karst cave distribution should be paid attention to，which was beneficial to reducing the influence of karst caves on single pile bearing capacity. The proposed calculation method can provide a theoretical basis for parameter design and estimation of the bearing capacity of a single pile in the karst area.

Key words：karst area；geophysical prospecting method；single pile crossing karst cave；bearing characteristics；numerical simulation

樁基礎(chǔ)作為一種常見的基礎(chǔ)形式，具有適應(yīng)性廣、承載力強(qiáng)、可減小不均勻沉降等優(yōu)點(diǎn)，在巖溶地區(qū)的施工建設(shè)中被廣泛使用[1-2].嵌巖深度和頂板厚度不滿足設(shè)計要求時樁基需要穿越溶洞的情況屢見不鮮，因此亟待開展穿越溶洞型樁基的相關(guān)研究.

巖溶區(qū)樁承載特性與常規(guī)嵌巖樁不同，Ladanyi[3]、尹平保等[4]和趙明華等[5]通過理論和實(shí)測數(shù)據(jù)分析了嵌巖樁嵌巖深度，給出了不同破壞模式下，嵌巖樁極限承載力的計算方法；劉會球[6]、Poulos等[7]和龔先兵等[8]分析了影響嵌巖樁樁基承載力的因素，得到嵌巖深度、樁的長徑比影響荷載傳遞規(guī)律及樁端下伏溶洞頂板穩(wěn)定性的分析方法；黃明等[9]、鄒新軍等[10]和黎斌等[11]通過有限元、多元線性回歸等手段分析了溶洞頂板的破壞形式、溶洞尺寸和樁基設(shè)計荷載與頂板安全厚度的關(guān)系；李天雨等[12]開展了水平荷載作用于巖溶覆蓋層臨空面室內(nèi)模型試驗(yàn)，揭示了樁前臨空面對樁嵌固端承載力的影響并依據(jù)試驗(yàn)結(jié)果提出了對應(yīng)理論模型；Jiang等[13]、魏鋒等[14]和Yacine等[15]通過巖石力學(xué)試驗(yàn)、數(shù)值模擬和理論分析，發(fā)現(xiàn)溶洞穩(wěn)定性與溶洞大小和相關(guān)地質(zhì)力學(xué)參數(shù)有關(guān)，確定了溶洞頂板破壞的主要影響因素及溶洞頂板安全厚度.目前對于巖溶區(qū)樁基承載特性及穩(wěn)定性評價方法的研究多集中于樁端隱伏溶洞，未能考慮樁基穿越溶洞嵌入基巖的情況，缺乏對于穿越溶洞型單樁沉降計算方法的研究.

本文依托濰煙鐵路昌平特大橋樁基工程，根據(jù)現(xiàn)場物探結(jié)果確定數(shù)值模擬模型形狀，通過現(xiàn)場取樣、室內(nèi)試驗(yàn)確定巖土體的物理力學(xué)參數(shù)，然后采用有限元軟件分析了溶洞半徑、溶洞位置、樁徑等對巖溶區(qū)穿越溶洞型單樁承載特性的影響規(guī)律，并以此確定荷載傳遞模型——BoxLucasl模型，建立穿越無充填溶洞時單樁沉降迭代計算方法，以指導(dǎo)工程實(shí)踐、降低工程風(fēng)險.

1現(xiàn)場鉆探與物探結(jié)果分析

濰煙鐵路昌平特大橋標(biāo)段內(nèi)長28 525 m，需跨越巖溶發(fā)育區(qū)，采用大直徑嵌巖灌注樁基礎(chǔ)類型.本文選定昌平特大橋內(nèi)巖溶較為發(fā)育的445號橋墩（8樁承臺）的4根樁進(jìn)行勘察，通過鉆孔法、管波探測法等組合方法確定巖溶發(fā)育情況，為數(shù)值模型建立提供參考.物探結(jié)果表明（見圖1），樁位區(qū)巖溶發(fā)育，樁體需穿越單溶洞，勘查區(qū)溶洞深度為30～40 m，溶洞半徑為5～10m，其中5號角樁巖溶發(fā)育情況最具代表性.

2現(xiàn)場巖芯物理力學(xué)參數(shù)測定

本文采用TOP INDUSTRIE STAC 600-600多功能巖石三軸試驗(yàn)系統(tǒng)，對現(xiàn)場取出的10塊巖樣（巖樣為直徑5 cm、高10 cm的標(biāo)準(zhǔn)圓柱形試樣）進(jìn)行試驗(yàn).

結(jié)合三軸試驗(yàn)結(jié)果，采用作圖法求解內(nèi)聚力c和內(nèi)摩擦角φ，切點(diǎn)應(yīng)力對應(yīng)破壞的正應(yīng)力σ和剪應(yīng)力τ，近似滿足線性關(guān)系式（1）.

τ=c+σtanφ（1）

以最小主應(yīng)力σ為自變量x，最大主應(yīng)力σ為因變量y，采用最小二乘法y=ax+b確定回歸方程，常數(shù)a、b和c、φ的關(guān)系為：

sinφ=（b-1）/（b+1）（2）

c=a（1-sinφ）/2cosφ（3）

σ-σ關(guān)系曲線如圖2所示.根據(jù)σ-σ應(yīng)力應(yīng)變關(guān)系曲線，利用式（2）和式（3）求得內(nèi)聚力c、內(nèi)摩擦角φ分別為7 250 kPa和53.13°.

利用式（4）求解彈性模量，即：

式中：E為彈性模量；（σ-σ）為主應(yīng)力差，是單軸抗壓強(qiáng)度50%時的應(yīng)力值；ε為最大主應(yīng)力，是單軸抗壓強(qiáng)度50%時的縱向應(yīng)變值.

取最大主應(yīng)力為單軸抗壓強(qiáng)度50%時的縱向應(yīng)變值ε和橫向應(yīng)變值ε計算泊松比v，即：

根據(jù)式（4）（5）求得巖石試樣的彈性模量E和泊松比v.巖樣的物理力學(xué)參數(shù)見表1.

3數(shù)值模型建立

根據(jù)物探結(jié)果及設(shè)計要求，采用有限元數(shù)值模擬軟件建立單樁穿越單溶洞數(shù)值模型，模型主體包含4個部分：土體、巖體、樁和溶洞.初始巖土體尺寸設(shè)計為：20 m×20 m×50 m，其中土體高20 m，巖體高30 m，樁長40 m.采用分級加載方式，共加載6次，每級施加5 MPa.本文共建立14組數(shù)值分析模型，分別設(shè)計了不同溶洞半徑、溶洞位置、樁徑等工況，具體設(shè)計方案見表2.

本文充填溶洞為人工注漿填充，采用片石回填的方法.溶洞充填材料及樁周巖土體采用Mohr- Coulom模型模擬，樁體采用彈性模型.根據(jù)第2節(jié)室內(nèi)試驗(yàn)結(jié)果，數(shù)值模擬材料參數(shù)如表3所示.

將樁土界面接觸定義為各向同性的庫倫摩擦模型，主從界面間允許滑動和分離，摩擦特性在法向上假定為硬接觸，允許分離，在切向上定義為“罰”接觸.摩擦因數(shù)μ可由式（6）和式（7）確定：

μ=tanδ（6）

tanδ=（0.5～1.0）tanφ′（7）

式中：δ為樁與地層摩擦角；φ′為地層內(nèi)摩擦角.本研究中取摩擦因數(shù)與地層系數(shù)的比值為0.75.

樁與地層接觸設(shè)置為面-面接觸，在相互作用中選擇表面-表面的離散方法，計算公式選擇為有限滑移.

4模型結(jié)果分析

4.1不同溶洞半徑

圖3為不同溶洞充填情況及溶洞半徑下地層- 溶洞一樁位移云圖.由圖3可知樁身沉降遠(yuǎn)大于樁周土體沉降，且以樁為對稱中心向土體方向逐漸減小. 無充填情況下，巖土體沉降主要發(fā)生在溶洞底板以上，空溶洞對荷載產(chǎn)生屏蔽作用，使溶洞下方巖土體難以發(fā)揮承載作用，且會產(chǎn)生應(yīng)力集中.由于溶洞充填材料相較于樁身和巖體剛度較小，在樁受荷產(chǎn)生沉降時充填溶洞會產(chǎn)生一定程度的拖空現(xiàn)象，引起局部應(yīng)力增大，使云圖中溶洞段及周圍巖土體位移大于同等深度下巖土體位移.在相同荷載情況下，樁身的沉降隨溶洞半徑的增大而增大，無充填溶洞樁體的位移量明顯大于充填溶洞，溶洞充填可有效改善單樁承載特性.

無充填溶洞情況下溶洞半徑3 m時不同荷載情況下樁側(cè)摩阻力及樁軸力的分布情況，如圖4和圖5所示.

由圖4可知，在荷載較小時（3 927 kN、7 854 kN），樁周土層側(cè)摩阻力尚未完全發(fā)揮；隨著樁頂荷載增大，上覆土層的側(cè)摩阻力逐漸發(fā)揮，側(cè)摩阻力在巖土體分界面處產(chǎn)生突變.對于溶洞段，由于無充填溶洞上頂板巖體較為破碎，導(dǎo)致上部巖體沉降大于樁基，出現(xiàn)負(fù)摩阻力，當(dāng)荷載增大到15 708 kN時，不再有負(fù)摩擦力產(chǎn)生.

由圖5可知，樁身軸力整體沿深度逐漸減小，且隨荷載的增加軸力分布情況逐漸表現(xiàn)為非線性分布.在荷載較小時（3 927 kN、7 854 kN），由于溶洞頂板處存在負(fù)側(cè)摩阻力，位于溶洞上部的樁身軸力局部增大，溶洞底板至樁端區(qū)域軸力下降速度明顯加快，這與摩擦因數(shù)突變，巖層側(cè)摩阻力得以突然發(fā)揮有關(guān).

根據(jù)不同荷載-沉降關(guān)系繪制樁頂及樁端的荷載-沉降曲線，見圖6和圖7.

由圖6可知，在荷載達(dá)到15 000 kN前曲線偏線性，隨著荷載增大，荷載-沉降曲線產(chǎn)生拐點(diǎn).溶洞半徑相同時，無充填溶洞單樁承載力小.而在樁端荷載-沉降曲線（圖7）中，曲線早期斜率較小，但更早到達(dá)拐點(diǎn)，進(jìn)入陡降階段.可見當(dāng)荷載較大時樁端荷載對單樁承載力起決定性作用.同時在相同荷載情況下，溶洞半徑越小單樁沉降越小，無充填溶洞2 m時的荷載沉降曲線位于充填溶洞半徑3 m時荷載- 沉降曲線的下方.根據(jù)溶洞的大小選擇合理的溶洞充填方式和充填材料對于樁體承載力的提高有著顯著作用.

4.2不同樁徑

不同溶洞充填情況及不同樁徑下地層-溶洞-樁位移云圖如圖8所示.

由圖8可知，樁徑增大顯著降低了樁體的沉降；由于巖體強(qiáng)度遠(yuǎn)大于土體，其沉降相對于土體較小，在云圖上表現(xiàn)為位移場以樁為中心對稱向底層擴(kuò)散，且隨著樁徑的增加，擴(kuò)散范圍越大，位移場的連續(xù)性越好；溶洞區(qū)的應(yīng)力集中現(xiàn)象由于溶洞的充填和樁徑的增大明顯減弱，云圖中充填溶洞樁徑1m時溶洞段樁體的變形程度小于無充填溶洞樁徑1.5m時溶洞段樁體的變形程度.結(jié)合上一節(jié)分析，溶洞充填和樁徑增大均可改善單樁承載特性，在穿越溶洞型單樁設(shè)計時應(yīng)協(xié)調(diào)選擇.

不同溶洞充填和樁徑情況下樁頂荷載-沉降曲線和樁端荷載-沉降曲線見圖9和圖10.

由圖9和圖10可知，樁徑越大單樁沉降越小，按照位移控制單樁極限承載力的原則，樁徑越大其承載性能也越強(qiáng)；加載初期，荷載-沉降曲線為線性分布，樁徑越大斜率越小，也越難到達(dá)拐點(diǎn)進(jìn)入陡降階段；無充填溶洞與充填溶洞相比，荷載-沉降曲線變化趨勢相同，無充填溶洞荷載-沉降曲線位于充填溶洞曲線下方.

4.3不同溶洞位置

圖11和圖12分別是溶洞無充填情況下溶洞距離樁頂25 m和30 m時的樁側(cè)摩阻力分布圖.

由圖11、圖12可知，溶洞位置不同，樁側(cè)摩阻力的發(fā)揮形式改變明顯.溶洞距離樁頂30 m時負(fù)側(cè)摩阻力的分布范圍明顯大于溶洞距離樁頂25 m時；溶洞越靠近樁端，溶洞附近產(chǎn)生的應(yīng)力集中場（溶洞上下頂板）與樁端應(yīng)力場相互疊加越明顯，導(dǎo)致樁端側(cè)摩阻力突變加劇.

不同溶洞位置下單樁穿越溶洞時樁頂、樁端荷載-沉降曲線，如圖13與圖14所示.

當(dāng)荷載較小時，相同充填情況不同溶洞位置的樁頂與樁端荷載-沉降曲線幾乎完全重合；相同條件下無充填溶洞沉降大于充填溶洞沉降，整體上溶洞越靠近樁端，單樁承載力越小.溶洞不同位置主要通過影響側(cè)摩阻力和軸力的分布規(guī)律影響樁基承載性狀，體現(xiàn)在荷載-沉降曲線上即樁端沉降量差距較大，根據(jù)圖14可知，不同溶洞位置對樁頂沉降量影響較小.

5樁基穿越無充填溶洞沉降計算方法

5.1荷載傳遞模型

根據(jù)數(shù)值模擬中側(cè)摩阻力的應(yīng)力分布圖，將荷載傳遞模型分為樁-土界面、樁-巖界面、巖-土體分界面、溶洞上頂板、溶洞下底板和樁端等6個部分. 采用荷載傳遞法，并用BoxLucas1模型[16]分別對其在數(shù)值模擬提取不同工況下的側(cè)摩阻力τ與相對位移S數(shù)據(jù)進(jìn)行線性回歸擬合.如圖15-圖20所示，Boxlucas1模型計算值與數(shù)值模型提取值擬合結(jié)果良好，決定系數(shù)R均接近于1，BoxLucas1模型適合作為穿越無填充溶洞時單樁沉降計算的荷載傳遞模型.

荷載傳遞模型的樁側(cè)摩阻力函數(shù)表達(dá)式[17-18]為：

式中：τ為樁側(cè)單位摩阻力；S為樁端位移；m、n為參數(shù).m、n表達(dá)式為：

式中：r為樁身半徑；r為樁的影響半徑；G為樁側(cè)土體剪切模量；R為樁側(cè)土破壞比，表示樁土相對位移達(dá)到極限值對應(yīng)的樁側(cè)摩阻力與極限側(cè)阻差異程度，其值取0.83～0.97[19-23]；τ為樁側(cè)最大靜摩擦力.

樁端荷載傳遞的BoxLucas1模型，具體計算公式[22-23]如下：

式中：P為單位樁端阻力；S為樁端位移；m、n為參數(shù).m、n取值為：

式中：G為樁端土剪切模量；v為樁端土泊松比；r為樁體半徑；P為極限單位樁端阻力；R為樁端土破壞比，表示樁土位移端阻與極限端阻差異，其值取0.85～0.95[24].

5.2單樁沉降計算方法

根據(jù)上文選定的荷載傳遞模型，利用Matlab程序，本文建立了適應(yīng)穿越無充填溶洞時單樁沉降的迭代計算流程，具體步驟如下：

1）將單樁自樁頂?shù)綐抖?，根?jù)地層和溶洞分布情況從1到n分成n段.

2）假定一較小的樁端位移為S.

3）將假定的樁端位移S代入式（14），得到樁端阻力P：

4）假定n樁段中點(diǎn)位移為S，假設(shè)S=S，將S代入式（15）得到樁段n的側(cè)摩阻力τ.

5）根據(jù)P和τ得到樁段n的樁頂荷載P：

P=P+πdLτ（16）

式中：L為第n段樁長；d為樁徑.

6）假定樁段n的樁身軸力線性變化，樁段n中點(diǎn)的彈性壓縮量S表示為：

式中：A為樁身截面面積；E為樁體彈性模量.

樁段n的中點(diǎn)修正位移S′可表示為：

8）樁段n的樁頂位移S和樁頂荷載P可分別表示為：

S=S+S′（19）

P=P+πdLτ′_sn（20）

9）讓P=P，S=S，重復(fù)步驟3）～8），計算其余樁段的樁頂位移和荷載，直至得到溶洞底板下處樁段的樁頂位移S和樁頂荷載P.

10）將溶洞段樁周材料參數(shù)替換為溶洞充填物相關(guān)參數(shù)，假定溶洞底板上處樁段i中點(diǎn)位移為S，假設(shè)S=S，將S代入式（15）得到樁段i的側(cè)摩阻力τ，根據(jù)步驟5）計算i段的樁頂荷載P，根據(jù)步驟6）和式（22）計算i段的中點(diǎn)彈性壓縮量S，根據(jù)步驟7）計算樁段i中點(diǎn)修正位移S′，重復(fù)步驟8）.

式中：L為溶洞段第i段樁段的樁長.

11）樁段i的樁頂位移S和樁頂荷載P可分別表示為：

S=S+S′（24）

P=P+πdLτ′（25）

12）令P=P，S=S，重復(fù)步驟10）、11），計算溶洞段其余樁段的樁頂位移和荷載，直至得到溶洞頂板下處樁段的樁頂位移S和樁頂荷載P.

5.3算例驗(yàn)證

算例來自彭戈[25]關(guān)于某橋墩實(shí)測數(shù)據(jù)，該樁長40 m，樁徑1.5 m，地表5 m內(nèi)為粉質(zhì)黏土與碎石混合物，下部45 m為中風(fēng)化灰?guī)r，溶洞寬度3 m，頂板厚度4 m，實(shí)測數(shù)據(jù)相關(guān)參數(shù)見表4、表5.

穿越無充填溶洞時單樁沉降計算結(jié)果見圖21.由圖21可知，本文計算結(jié)果與數(shù)值模擬結(jié)果較為一致.數(shù)值模擬得到的樁體極限承載力為13 950 kN，實(shí)測得到的豎向極限承載力為16 022 kN，本文迭代方法獲得的單樁極限承載力為16 599 kN（以40 mm 時樁頂荷載作為單樁極限承載力）.由表6中的對比結(jié)果可知，本文建立的計算方法的計算值同算例實(shí)測值對比誤差僅為3.6%，具有較好的一致性.

6結(jié)論

為保證數(shù)值模型結(jié)果的合理性，本文依托實(shí)際工程，根據(jù)物探與室內(nèi)試驗(yàn)結(jié)果構(gòu)建數(shù)值模型，通過數(shù)值模型分析結(jié)果選取合適的荷載傳遞模型，在此基礎(chǔ)上提出穿越無充填溶洞時單樁沉降迭代計算方法，將計算結(jié)果與工程檢測情況對比，獲得了以下結(jié)論：

1）樁基及周圍巖土體整體呈現(xiàn)分層及擴(kuò)散分布，樁徑大小及溶洞充填情況對應(yīng)力和位移分布影響較大，適當(dāng)增大樁徑及充填溶洞可有效改善樁基的承載性能.

2）在外荷載作用下，溶洞區(qū)域上下頂板會出現(xiàn)一定程度的應(yīng)力位移集中現(xiàn)象，溶洞半徑大小及溶洞位置不同，更多的是影響樁側(cè)摩阻力發(fā)揮，在實(shí)際工程中，根據(jù)前期地質(zhì)預(yù)測對巖溶發(fā)育關(guān)鍵部分進(jìn)行預(yù)處理是規(guī)避巖溶病害的有效手段.

3）將樁段分為樁-土界面、樁-巖界面、巖-土體分界面、溶洞上頂板、溶洞下底板、樁端等6個部分，各部分均可較好滿足BoxLucas1模型.本文建立的穿越無充填溶洞時單樁沉降迭代計算方法的計算值與算例值擬合較好，可以有效指導(dǎo)工程實(shí)踐.

參考文獻(xiàn)

[1]吳忠杰，王浩，鄧家喜.巖溶樁基施工地質(zhì)模式的劃分及工程應(yīng)用[J].西部交通科技，2010（11）：33-37.

WU Z J，WANG H，DENG J X. The classification and engineering application of geologic mode of Karst pile foundation [J]. Western China Communications Science & Technology，2010 （11）：33-37. （In Chinese）

[2]張琨，劉建民，鄔歡，等.城市橋梁樁基穿越溶洞預(yù)注漿施工技術(shù)[J].施工技術(shù)，2011，40（1）：77-79.

ZHANG K，LIU J M，WU H，et al. Pre-grouting technology used in pile of urban bridge crossing Karst cave [J]. Construction Technology，2011，40（1）：77-79.（In Chinese）

[3]LADANYI B. Friction and end bearing tests on bedrock for high capacity socket design[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts，1978，15（1）：123-127.

[4]尹平保，楊瑩，雷勇，等.考慮基巖水平破壞的嵌巖樁嵌巖深度計算方法[J].工程科學(xué)與技術(shù)，2019，51（2）：71-77.

YIN P B，YANG Y，LEI Y，et al. Calculation method of rock- socketed depth for rock-socketed piles based on bedrock horizontal failure [J]. Advanced Engineering Sciences，2019，51（2）：7177. （In Chinese）

[5]趙明華，肖堯，徐卓君，等.巖溶區(qū)嵌巖樁樁端承載性能研究[J].巖土工程學(xué)報，2017，39（6）：1123-1129.

ZHAO M H，XIAO Y，XU Z J，et al. Bearing capacity at tip of rock-socketed pile in Karst areas [J]. Chinese Journal of （Geotechnical Engineering，2017，39（6）：1123-1129. （In Chinese）

[6]劉會球.嵌巖樁承載力影響因素數(shù)值分析[J].鐵道科學(xué)與工程學(xué)報，2018，15（10）：2535-2540.

LIU H Q. Numerical analysis of rock socketed pile [J].Journal of Railway Science and Engineering，2018，15（10）：2535-2540. （In Chinese）

[7]POULOS W，PRICE D，TARR K. Jacked piles in London clay a study of load transfer and settlement under working conditions[J]. Geotechnique，1968，5（2）：839-848.

[8]龔先兵，趙明華，張永杰.樁端下伏溶洞頂板穩(wěn)定非概率可靠性分析方法[J].湖南大學(xué)學(xué)報（自然科學(xué)版），2011，38（1）：13-17.

GONG X B，ZHAO M H，ZHANG Y J. Non-probabilistic reliability analysis method of Karst roof stability under pile tip [J]. Journal of Hunan University （Natural Sciences），2011，38（1）：13-17. （In Chinese）

[9]黃明，張冰淇，陳福全，等.串珠狀溶洞地層中樁基荷載傳遞特征的數(shù)值計算[J].工程地質(zhì)學(xué)報，2017，25（6）：1574-1582.

HUANG M，ZHANG B Q，CHEN F Q，et al. Numerical calculation on load transfer progress of pile foundation in beaded Karst cave stratum[J]. Journal of Engineering Geology，2017，25（6）：1574-1582.（In Chinese）

[10]鄒新軍，唐國東，趙明華.串珠狀巖溶區(qū)樁基沉降計算與穩(wěn)定分析[J].建筑結(jié)構(gòu)，2013，43（13）：95-98.

ZOU X J，TANG G D，ZHAO M H.Settlement and stability analysis of pile foundation in string-beaded Karst area[J]. Building Structure，2013，43（13）：95-98.（In Chinese）

[11]黎斌，范秋雁，秦鳳榮.巖溶地區(qū)溶洞頂板穩(wěn)定性分析[J].巖石力學(xué)與工程學(xué)報，2002，21（4）：532-536.

LI B，F(xiàn)AN Q Y，QIN F R.Analysis on roof stability of Karst cave in Karst areas[J].Chinese Journal of Rock Mechanics and Engi- neering，2002，21（4）：532-536.（In Chinese）

[12]李天雨，范秋雁，韓偉，等.覆蓋巖溶臨空面對水平承載樁嵌固端承載力的影響試驗(yàn)研究[J].工程力學(xué)，2020，37 （10）：192-199.

LI T Y，F(xiàn)AN Q Y，HAN W，et al. Experimental study on influence of covered karstic free surface on fixed end bearing capacity of laterally loaded pile[J]. Engineering Mechanics，2020，37 （10）：192-199.（In Chinese）

[13]JIANG C，ZHAO M H，CAO W G. Stability analysis of subgrade cave roofs in Karst region[J].Journal of Central South University of Technology，2008，15（2）：38-44.

[14]魏鋒，陳忠達(dá)，陳峙峰，等.路堤下伏溶洞受力模式和變形破壞的數(shù)值模擬[J].中國公路學(xué)報，2018，31（6）：195-206.

WEI F，CHEN Z D，CHEN Z F，et al.Numerical simulation of themechanical characteristic and failure mode of Karst subgrade[J]. China Journal of Highway and Transport，2018，31（6）：195-206. （In Chinese）

[15]YACINE A，RIDHA M M，LAID H M，et al.Karst sinkholes stability assessment in cheria area，NE Algeria[J]. Geotechnical and Geological Engineering，2014，32（2）：363-374.

[16]WANG Z J，XIE X Y，WANG J C. A new nonlinear method for vertical settlement prediction of a single pile and pile groups in layered soils [J]. Computers and Geotechnics，2012，45：118-126.

[17]VESIC A S. Expansion of cavities in infinite soil mass [J]. Journal of the Soil Mechanics and Foundations Division，1972，98（3）：265-290.

[18]MANANDHAR S，YASUFUKU N. Vertical bearing capacity of tapered piles in sands using cavity expansion theory[J]. Soils and Foundations，2013，53（6）：853-867.

[19]李連祥，李先軍，成曉陽，等.考慮圓孔擴(kuò)張理論的支盤樁荷載傳遞法[J].中國公路學(xué)報，2018，31（8）：20-29.

LI L X，LI X J，CHENG X Y，et al. Load transfer method for squeezed and branch piles considering cavity expansion theory [J]. China Journal of Highway and Transport，2018，31（8）：20- 29.（In Chinese）

[20]張乾青，張忠苗.樁基工程[M].2版.北京：中國建筑工業(yè)出版社，2018.

ZHANG Q Q，ZHANG Z M. Pile foundation engineering [M]. 2nd ed. Beijing：China Architecture & Building Press，2018. （In Chinese）

[21]MINDLIN R D. Influence of rotatory inertia and shear on flexural motions of isotropic，elastic plates[J]. Journal of Applied Mechanics，1951，18（1）：31-38.

[22]ZHANG Q，ZHANG Z，YU F，et al. Field performance of long bored piles within piled rafts[J].Proceedings of the Institution of Civil Engineers-Geotechnical Engineering，2010，163 （6）：293-305.

[23]ZHANG Z M，ZHANG Q Q，YU F. A destructive field study on the behavior of piles under tension and compression[J]. Journal of Zhejiang University：Science A，2011，12（4）：291-300.

[24]李永輝，吳江斌.基于載荷試驗(yàn)的大直徑超長樁承載特性分析[J].地下空間與工程學(xué)報，2011，7（5）：895-902.

LI Y H，WU J B. Bearing behavior study of large-diameter and super-long pile based on full-scale field test[J]. Chinese Journal of Underground Space and Engineering，2011，7（5）：895-902. （In Chinese）

[25]彭戈.穿越多層溶洞橋梁樁基的豎向承載力及沉降變形研究[D].成都：西南交通大學(xué)，2016.

PENG G. Study on vertical bearing capacity and settlement of bridge pile foundation passes through multi-Karst-caves [D]. Chengdu：Southwest Jiaotong University，2016.（In Chinese）