穿越无充填溶洞时单桩承载特性研究

张乾青 乔胜石 邢宇铖 张凯 崔伟 王志远 陈迪杨

摘要：为解决穿越溶洞型单桩的设计难题，研究了穿越无充填溶洞时单桩的承载特性. 基于现场物探和室内试验结果明确了溶洞的发育情况和岩土体的物理力学参数，采用有限元软件建立了14组数值计算模型，分别对不同溶洞半径、不同溶洞位置及不同桩径的工况进行了模拟，确定了BoxLucas1荷载传递模型的适用性，提出了穿越无充填溶洞单桩的承载特性计算方法，研究结果表明，桩径增加和溶洞充填时溶洞对单桩承载力的影响将得到削弱，溶洞的半径和位置对桩基承载力具有显著影响.实际工程中应重视前期物探对溶洞发育情况的准确探查，降低溶洞对桩基承载力的影响.相关研究成果可为岩溶区单桩参数设计和承载力预测提供理论依据.

关键词：岩溶区；物探方法；穿越溶洞型单桩；承载特性；数值模拟

中图分类号：TU473.1文献标志码：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

桩基础作为一种常见的基础形式，具有适应性广、承载力强、可减小不均匀沉降等优点，在岩溶地区的施工建设中被广泛使用[1-2].嵌岩深度和顶板厚度不满足设计要求时桩基需要穿越溶洞的情况屡见不鲜，因此亟待开展穿越溶洞型桩基的相關研究.

岩溶区桩承载特性与常规嵌岩桩不同，Ladanyi[3]、尹平保等[4]和赵明华等[5]通过理论和实测数据分析了嵌岩桩嵌岩深度，给出了不同破坏模式下，嵌岩桩极限承载力的计算方法；刘会球[6]、Poulos等[7]和龚先兵等[8]分析了影响嵌岩桩桩基承载力的因素，得到嵌岩深度、桩的长径比影响荷载传递规律及桩端下伏溶洞顶板稳定性的分析方法；黄明等[9]、邹新军等[10]和黎斌等[11]通过有限元、多元线性回归等手段分析了溶洞顶板的破坏形式、溶洞尺寸和桩基设计荷载与顶板安全厚度的关系；李天雨等[12]开展了水平荷载作用于岩溶覆盖层临空面室内模型试验，揭示了桩前临空面对桩嵌固端承载力的影响并依据试验结果提出了对应理论模型；Jiang等[13]、魏锋等[14]和Yacine等[15]通过岩石力学试验、数值模拟和理论分析，发现溶洞稳定性与溶洞大小和相关地质力学参数有关，确定了溶洞顶板破坏的主要影响因素及溶洞顶板安全厚度.目前对于岩溶区桩基承载特性及稳定性评价方法的研究多集中于桩端隐伏溶洞，未能考虑桩基穿越溶洞嵌入基岩的情况，缺乏对于穿越溶洞型单桩沉降计算方法的研究.

本文依托潍烟铁路昌平特大桥桩基工程，根据现场物探结果确定数值模拟模型形状，通过现场取样、室内试验确定岩土体的物理力学参数，然后采用有限元软件分析了溶洞半径、溶洞位置、桩径等对岩溶区穿越溶洞型单桩承载特性的影响规律，并以此确定荷载传递模型——BoxLucasl模型，建立穿越无充填溶洞时单桩沉降迭代计算方法，以指导工程实践、降低工程风险.

1现场钻探与物探结果分析

潍烟铁路昌平特大桥标段内长28 525 m，需跨越岩溶发育区，采用大直径嵌岩灌注桩基础类型.本文选定昌平特大桥内岩溶较为发育的445号桥墩（8桩承台）的4根桩进行勘察，通过钻孔法、管波探测法等组合方法确定岩溶发育情况，为数值模型建立提供参考.物探结果表明（见图1），桩位区岩溶发育，桩体需穿越单溶洞，勘查区溶洞深度为30～40 m，溶洞半径为5～10m，其中5号角桩岩溶发育情况最具代表性.

2现场岩芯物理力学參数测定

本文采用TOP INDUSTRIE STAC 600-600多功能岩石三轴试验系统，对现场取出的10块岩样（岩样为直径5 cm、高10 cm的标准圆柱形试样）进行试验.

结合三轴试验结果，采用作图法求解内聚力c和内摩擦角φ，切点应力对应破坏的正应力σ和剪应力τ，近似满足线性关系式（1）.

τ=c+σtanφ（1）

以最小主应力σ为自变量x，最大主应力σ为因变量y，采用最小二乘法y=ax+b确定回归方程，常数a、b和c、φ的关系为：

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

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

σ-σ关系曲线如图2所示.根据σ-σ应力应变关系曲线，利用式（2）和式（3）求得内聚力c、内摩擦角φ分别为7 250 kPa和53.13°.

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

式中：E为弹性模量；（σ-σ）为主应力差，是单轴抗压强度50%时的应力值；ε为最大主应力，是单轴抗压强度50%时的纵向应变值.

取最大主应力为单轴抗压强度50%时的纵向应变值ε和横向应变值ε计算泊松比v，即：

根据式（4）（5）求得岩石试样的弹性模量E和泊松比v.岩样的物理力学参数见表1.

3数值模型建立

根据物探结果及设计要求，采用有限元数值模拟软件建立单桩穿越单溶洞数值模型，模型主体包含4个部分：土体、岩体、桩和溶洞.初始岩土体尺寸设计为：20 m×20 m×50 m，其中土体高20 m，岩体高30 m，桩长40 m.采用分级加载方式，共加载6次，每级施加5 MPa.本文共建立14组数值分析模型，分别设计了不同溶洞半径、溶洞位置、桩径等工况，具体设计方案见表2.

本文充填溶洞为人工注浆填充，采用片石回填的方法.溶洞充填材料及桩周岩土体采用Mohr- Coulom模型模拟，桩体采用弹性模型.根据第2节室内试验结果，数值模拟材料参数如表3所示.

将桩土界面接触定义为各向同性的库伦摩擦模型，主从界面间允许滑动和分离，摩擦特性在法向上假定为硬接触，允许分离，在切向上定义为“罚”接触.摩擦因数μ可由式（6）和式（7）确定：

μ=tanδ（6）

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

式中：δ为桩与地层摩擦角；φ′为地层内摩擦角.本研究中取摩擦因数与地层系数的比值为0.75.

桩与地层接触设置为面-面接触，在相互作用中选择表面-表面的离散方法，计算公式选择为有限滑移.

4模型结果分析

4.1不同溶洞半径

图3为不同溶洞充填情况及溶洞半径下地层- 溶洞一桩位移云图.由图3可知桩身沉降远大于桩周土体沉降，且以桩为对称中心向土体方向逐渐减小. 无充填情况下，岩土体沉降主要发生在溶洞底板以上，空溶洞对荷载产生屏蔽作用，使溶洞下方岩土体难以发挥承载作用，且会产生应力集中.由于溶洞充填材料相较于桩身和岩体刚度较小，在桩受荷产生沉降时充填溶洞会产生一定程度的拖空现象，引起局部应力增大，使云图中溶洞段及周围岩土体位移大于同等深度下岩土体位移.在相同荷载情况下，桩身的沉降随溶洞半径的增大而增大，无充填溶洞桩体的位移量明显大于充填溶洞，溶洞充填可有效改善单桩承载特性.

无充填溶洞情况下溶洞半径3 m时不同荷载情况下桩侧摩阻力及桩轴力的分布情况，如图4和图5所示.

由图4可知，在荷载较小时（3 927 kN、7 854 kN），桩周土层侧摩阻力尚未完全发挥；随着桩顶荷载增大，上覆土层的侧摩阻力逐渐发挥，侧摩阻力在岩土体分界面处产生突变.对于溶洞段，由于无充填溶洞上顶板岩体较为破碎，导致上部岩体沉降大于桩基，出现负摩阻力，当荷载增大到15 708 kN时，不再有负摩擦力产生.

由图5可知，桩身轴力整体沿深度逐渐减小，且随荷载的增加轴力分布情况逐渐表现为非线性分布.在荷载较小时（3 927 kN、7 854 kN），由于溶洞顶板处存在负侧摩阻力，位于溶洞上部的桩身轴力局部增大，溶洞底板至桩端区域轴力下降速度明显加快，这与摩擦因数突变，岩层侧摩阻力得以突然发挥有关.

根据不同荷载-沉降关系绘制桩顶及桩端的荷载-沉降曲线，见图6和图7.

由图6可知，在荷载达到15 000 kN前曲线偏线性，随着荷载增大，荷载-沉降曲线产生拐点.溶洞半径相同时，无充填溶洞单桩承载力小.而在桩端荷载-沉降曲线（图7）中，曲线早期斜率较小，但更早到达拐点，进入陡降阶段.可见当荷载较大时桩端荷载对单桩承载力起决定性作用.同时在相同荷载情况下，溶洞半径越小单桩沉降越小，无充填溶洞2 m时的荷载沉降曲线位于充填溶洞半径3 m时荷载- 沉降曲线的下方.根据溶洞的大小选择合理的溶洞充填方式和充填材料对于桩体承载力的提高有着显著作用.

4.2不同桩径

不同溶洞充填情况及不同桩径下地层-溶洞-桩位移云图如图8所示.

由图8可知，桩径增大显著降低了桩体的沉降；由于岩体强度远大于土体，其沉降相对于土体较小，在云图上表现为位移场以桩为中心对称向底层扩散，且随着桩径的增加，扩散范围越大，位移场的连续性越好；溶洞区的应力集中现象由于溶洞的充填和桩径的增大明显减弱，云图中充填溶洞桩径1m时溶洞段桩体的变形程度小于无充填溶洞桩径1.5m时溶洞段桩体的变形程度.結合上一节分析，溶洞充填和桩径增大均可改善单桩承载特性，在穿越溶洞型单桩设计时应协调选择.

不同溶洞充填和桩径情况下桩顶荷载-沉降曲线和桩端荷载-沉降曲线见图9和图10.

由图9和图10可知，桩径越大单桩沉降越小，按照位移控制单桩极限承载力的原则，桩径越大其承载性能也越强；加载初期，荷载-沉降曲线为线性分布，桩径越大斜率越小，也越难到达拐点进入陡降阶段；无充填溶洞与充填溶洞相比，荷载-沉降曲线变化趋势相同，无充填溶洞荷载-沉降曲线位于充填溶洞曲线下方.

4.3不同溶洞位置

图11和图12分别是溶洞无充填情况下溶洞距离桩顶25 m和30 m时的桩侧摩阻力分布图.

由图11、图12可知，溶洞位置不同，桩侧摩阻力的发挥形式改变明显.溶洞距离桩顶30 m时负侧摩阻力的分布范围明显大于溶洞距离桩顶25 m时；溶洞越靠近桩端，溶洞附近产生的应力集中场（溶洞上下顶板）与桩端应力场相互叠加越明显，导致桩端侧摩阻力突变加剧.

不同溶洞位置下单桩穿越溶洞时桩顶、桩端荷载-沉降曲线，如图13与图14所示.

当荷载较小时，相同充填情况不同溶洞位置的桩顶与桩端荷载-沉降曲线几乎完全重合；相同条件下无充填溶洞沉降大于充填溶洞沉降，整体上溶洞越靠近桩端，单桩承载力越小.溶洞不同位置主要通过影响侧摩阻力和轴力的分布规律影响桩基承载性状，体现在荷载-沉降曲线上即桩端沉降量差距较大，根据图14可知，不同溶洞位置对桩顶沉降量影响较小.

5桩基穿越无充填溶洞沉降计算方法

5.1荷载传递模型

根据数值模拟中侧摩阻力的应力分布图，将荷载传递模型分为桩-土界面、桩-岩界面、岩-土体分界面、溶洞上顶板、溶洞下底板和桩端等6个部分. 采用荷载传递法，并用BoxLucas1模型[16]分别对其在数值模拟提取不同工况下的侧摩阻力τ与相对位移S数据进行线性回归拟合.如图15-图20所示，Boxlucas1模型计算值与数值模型提取值拟合结果良好，决定系数R均接近于1，BoxLucas1模型适合作为穿越无填充溶洞时单桩沉降计算的荷载传递模型.

荷载传递模型的桩侧摩阻力函数表达式[17-18]为：

式中：τ为桩侧单位摩阻力；S为桩端位移；m、n为参数.m、n表达式为：

式中：r为桩身半径；r为桩的影响半径；G为桩侧土体剪切模量；R为桩侧土破坏比，表示桩土相对位移达到极限值对应的桩侧摩阻力与极限侧阻差异程度，其值取0.83～0.97[19-23]；τ为桩侧最大静摩擦力.

桩端荷载传递的BoxLucas1模型，具体计算公式[22-23]如下：

式中：P为单位桩端阻力；S为桩端位移；m、n为参数.m、n取值为：

式中：G为桩端土剪切模量；v为桩端土泊松比；r为桩体半径；P为极限单位桩端阻力；R为桩端土破坏比，表示桩土位移端阻与极限端阻差异，其值取0.85～0.95[24].

5.2单桩沉降计算方法

根据上文选定的荷载传递模型，利用Matlab程序，本文建立了适应穿越无充填溶洞时单桩沉降的迭代计算流程，具体步骤如下：

1）将单桩自桩顶到桩端，根据地层和溶洞分布情况从1到n分成n段.

2）假定一较小的桩端位移为S.

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

4）假定n桩段中点位移为S，假设S=S，将S代入式（15）得到桩段n的侧摩阻力τ.

5）根据P和τ得到桩段n的桩顶荷载P：

P=P+πdLτ（16）

式中：L為第n段桩长；d为桩径.

6）假定桩段n的桩身轴力线性变化，桩段n中点的弹性压缩量S表示为：

式中：A为桩身截面面积；E为桩体弹性模量.

桩段n的中点修正位移S′可表示为：

8）桩段n的桩顶位移S和桩顶荷载P可分别表示为：

S=S+S′（19）

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

9）让P=P，S=S，重复步骤3）～8），计算其余桩段的桩顶位移和荷载，直至得到溶洞底板下处桩段的桩顶位移S和桩顶荷载P.

10）将溶洞段桩周材料参数替换为溶洞充填物相关参数，假定溶洞底板上处桩段i中点位移为S，假设S=S，将S代入式（15）得到桩段i的侧摩阻力τ，根据步骤5）计算i段的桩顶荷载P，根据步骤6）和式（22）计算i段的中点弹性压缩量S，根据步骤7）计算桩段i中点修正位移S′，重复步骤8）.

式中：L为溶洞段第i段桩段的桩长.

11）桩段i的桩顶位移S和桩顶荷载P可分别表示为：

S=S+S′（24）

P=P+πdLτ′（25）

12）令P=P，S=S，重复步骤10）、11），计算溶洞段其余桩段的桩顶位移和荷载，直至得到溶洞顶板下处桩段的桩顶位移S和桩顶荷载P.

5.3算例验证

算例来自彭戈[25]关于某桥墩实测数据，该桩长40 m，桩径1.5 m，地表5 m内为粉质黏土与碎石混合物，下部45 m为中风化灰岩，溶洞宽度3 m，顶板厚度4 m，实测数据相关参数见表4、表5.

穿越无充填溶洞时单桩沉降计算结果见图21.由图21可知，本文计算结果与数值模拟结果较为一致.数值模拟得到的桩体极限承载力为13 950 kN，实测得到的竖向极限承载力为16 022 kN，本文迭代方法获得的单桩极限承载力为16 599 kN（以40 mm 时桩顶荷载作为单桩极限承载力）.由表6中的对比结果可知，本文建立的计算方法的计算值同算例实测值对比误差仅为3.6%，具有较好的一致性.

6结论

为保证数值模型结果的合理性，本文依托实际工程，根据物探与室内试验结果构建数值模型，通过数值模型分析结果选取合适的荷载传递模型，在此基础上提出穿越无充填溶洞时单桩沉降迭代计算方法，将计算结果与工程检测情况对比，获得了以下结论：

1）桩基及周围岩土体整体呈现分层及扩散分布，桩径大小及溶洞充填情况对应力和位移分布影响较大，适当增大桩径及充填溶洞可有效改善桩基的承载性能.

2）在外荷载作用下，溶洞区域上下顶板会出现一定程度的应力位移集中现象，溶洞半径大小及溶洞位置不同，更多的是影响桩侧摩阻力发挥，在实际工程中，根据前期地质预测对岩溶发育关键部分进行预处理是规避岩溶病害的有效手段.

3）将桩段分为桩-土界面、桩-岩界面、岩-土体分界面、溶洞上顶板、溶洞下底板、桩端等6个部分，各部分均可较好满足BoxLucas1模型.本文建立的穿越无充填溶洞时单桩沉降迭代计算方法的计算值与算例值拟合较好，可以有效指导工程实践.

参考文献

[1]吴忠杰，王浩，邓家喜.岩溶桩基施工地质模式的划分及工程应用[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]张琨，刘建民，邬欢，等.城市桥梁桩基穿越溶洞预注浆施工技术[J].施工技术，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].工程科学与技术，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]赵明华，肖尧，徐卓君，等.岩溶区嵌岩桩桩端承载性能研究[J].岩土工程学报，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]刘会球.嵌岩桩承载力影响因素数值分析[J].铁道科学与工程学报，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]龚先兵，赵明华，张永杰.桩端下伏溶洞顶板稳定非概率可靠性分析方法[J].湖南大学学报（自然科学版），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]黄明，张冰淇，陈福全，等.串珠状溶洞地层中桩基荷载传递特征的数值计算[J].工程地质学报，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]邹新军，唐国东，赵明华.串珠状岩溶区桩基沉降计算与稳定分析[J].建筑结构，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]黎斌，范秋雁，秦凤荣.岩溶地区溶洞顶板稳定性分析[J].岩石力学与工程学报，2002，21（4）：532-536.

LI B，FAN 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]李天雨，范秋雁，韩伟，等.覆盖岩溶临空面对水平承载桩嵌固端承载力的影响试验研究[J].工程力学，2020，37 （10）：192-199.

LI T Y，FAN 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]魏锋，陈忠达，陈峙峰，等.路堤下伏溶洞受力模式和变形破坏的数值模拟[J].中国公路学报，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]李连祥，李先军，成晓阳，等.考虑圆孔扩张理论的支盘桩荷载传递法[J].中国公路学报，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版.北京：中国建筑工业出版社，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]李永辉，吴江斌.基于载荷试验的大直径超长桩承载特性分析[J].地下空间与工程学报，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].成都：西南交通大学，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）