群桩水平承载特性数值模拟研究_宋修广

铁道建筑ＲailwayEngineering

97

文章编号：1003-1995（2015）08-0097-04

群桩水平承载特性数值模拟研究

1，21，21，21，21，2

宋修广，岳红亚，周志东，庄培芝，张宏博

（1．山东大学土建与水利学院，山东济南250061；2．山东省路基安全工程技术研究中心，山东济南250061）

摘要：针对目前群桩基础水平承载特性研究落后于工程应用的现状，采用数值分析的方法研究了软土地

基管桩群桩基础在水平荷载作用下的承载特性。计算结果表明：群桩承受水平荷载时，承台前土体压缩变形范围明显小于承台后受拉变形范围；随着桩间距的增加，群桩侧向影响范围逐渐减小，而深度影响范围和桩顶位置处弯矩逐渐增大；相同水平荷载作用下前桩最大弯矩值明显大于后桩，前排桩桩身剪力零点位置高于后桩；当桩距达到6～8倍桩径时，群桩效应明显降低。关键词：群桩水平荷载管桩承载特性影响范围

+

中图分类号：TU473.11文献标识码：ADOI：10．3969/j．issn．1003-1995．2015．08．28软土地基建筑物尤其是海洋平台、输电塔、岸防工程和大跨度桥梁等，水平荷载成为设计控制因素。软土地基中建筑物的基础一般采用桩基形式，除少数采

［1-4］

。群桩用大直径单桩基础外，一般采用群桩基础基础受桩—土—承台的共同作用，在承受水平荷载时，

土体会产生应力和位移重叠区域，其水平承载特性较单桩复杂，群桩水平承载力并不等于各单桩承载力之［5-6］

。由于桩基属于地下工程，和现场试验受限于试验条件、环境、人为等因素的影响，往往只能得出群桩

难以揭示土体和桩体工作机理，因此水平承载力水平，

很多学者采用数值模拟的手段研究桩基受荷特性。课题组在文献中通过现场试验和数值计算比

较分析，确定了数值计算在管桩单桩水平计算中的合理性和适用性，并分析了管桩单桩水平承载特性的影

表1

编号②③④④-1⑤⑥⑦⑧⑨-1⑨土质名称粉土粉质黏土粉土粉质黏土粉土粉质黏土粉土粉质黏土粉土粉砂土层厚度/重度/

m（kN/m3）319.3419.31019.6318.1519.20.1520.1520.0319.7618.5

泊松比0.3

0.30.30.30.30.30.30.30.30.3

［7］

响因素，明确了各因素的主次关系。基于上述研究，本

文通过现场试验和数值模拟的对比分析，明确了数值模拟在群桩水平承载性能研究中的适用性和合理性。并采用数值模拟的手段，主要针对群桩中桩体的桩身位移及弯矩分布情况进行对比分析，以初步揭示群桩中桩—土—承台的荷载传递规律及协调变形特性，可为工程中群桩布置形式的优化提供参考。

1

1.1

群桩水平承载特性验证性分析

现场概况

试验区域位于鲁东地区黄河下游冲积平原，勘测揭

al

及第四系上露地层主要为第四系全新统冲击层（Q4），

al

更新统冲积层（Q4），岩性为粉土、粉质黏土和粉砂，局部地段分布有少量人工填土。物理性能指标见表1。

地基土物理力学性质指标压缩模量/MPa16.88.320.88.619.910.128.912.316.222.4

变形模量/MPa12.486.1715.456.3914.787.5021.479.1412.0316.

体积模量/MPa10.405.1412.885.3212.326.2517.7.6110.0313.87

切变模量/MPa4.802.375.942.465.692.8.263.514.636.40

黏聚力/kPa241423923222666843

内摩擦角/（°）22.913.226.216.028.217.530.411.224.826.8

收稿日期：2014-12-09；修回日期：2015-06-03

51208284）；山东省科技基金项目：国家自然科学基金项目（41372300，

发展计划项目（2013GSF11603）

作者简介：宋修广（1966—），男，山东威海人，教授，博士。

预应力管桩类别为AB桩，外径为600mm、壁厚

130mm，C80混凝土。群桩布置形式为桩长为25m，1×2型，桩间距为4D。采用慢速维持荷载法，荷载等级为50kN。由于试验桩为工程桩，为避免本次试验

98

铁道建筑August，2015

影响桩体以后工作性能，试验控制水平位移为10mm。

1.2数值计算模型的选取

采用有限差分程序FLAC3D进行计算分析，建立与试验桩相同尺寸的计算模型。桩体和承台采用弹性Coulomb弹塑性模型，模型；土体采用Mohr-模型的长、20和50m；桩—土界面采用Interface宽、高分别为40，命令，允许滑移和分离；用FLAC3D软件自带History命令对不同荷载作用下桩身水平位移进行跟踪监

［8-9］

。测1.3

试验结果对比分析

图1

荷载—位移曲线

2数值计算模型建立

受现场试验条件的，现场试验仅测量桩顶位移曲线，见图1。由图1可知现场试验和数值计算吻合度较高，同一水平荷载作用下数值计算对应桩顶水平位移较小，数值计算6mm位移对应水平承载力为370kN，为现场试验结果的1.06倍。通过比较分析，数值计算结果与现场试验结果具有较高的吻合度，数值计算模型选取和计算参数较合理，能够满足群桩效应综合系数的分析精度要求。

表2

桩径/

mm600

壁厚/mm130

桩长/m29

数值计算过程中假定承台与桩顶固结，并通过提高承台模量的方法予以反映。承台边缘与桩体最小距离等于0.5倍桩径，沿荷载方向模型长度为40m；竖向深度50m；垂直荷载方向保证承台边缘至模型边界AB-距离大于10倍桩径。选用桩型为PHC-600（130）-Lb型（桩径600mm、壁厚130mm、桩长29.0m），计算参数取值见表2。采用分级加载的方式，每级荷载取

［10］

1×2型群桩单级预估水平承载力的1/10～1/15，荷载为60kN。

群桩计算标准工况参数取值

土体密度/

（kg/m3）1700

土体黏聚力/

kPa

25

土体内摩擦角/

（°）

15

桩体弹性模量/承台弹性模量/土体压缩模量/

MPaMPaMPa3.8×104

2.85×105

7.0

3

3.1

计算结果及分析

土体影响范围

本文选取3D和9D桩间距计算工况，研究分析加

～2.5倍；③随着桩间距的增大，前后桩深度影响范围

逐渐增大。前桩的深度影响范围小于后桩，但两者相600kN荷载作用下位移影响深度约为10倍差不大，桩径。

表3

工况3D9D

桩土影响范围汇总

m

载到10级荷载即600kN时土体水平方向和深度方向影响范围，加载示意如图2。

水平方向影响范围侧向4..2

前桩4.37.2

后桩10.911.1

深度方向影响范围前桩5.65.8

后桩5.76.5

3.2

前、后桩承载力不均匀性分析

600kN水平荷载作用下，1×2桩型由图3可知，

图2加载示意

时，前桩最大弯矩值明显大于后桩，根据材料力学的基

本理论可知前排桩所受剪力（即桩前土体抗力）值大于后桩。桩距为3倍桩径时前桩桩顶弯矩（39.1kN·m）为后桩（5.63kN·m）的6.9倍，而桩距为9倍桩径·m）仅为后桩（29.5kN·m）时前桩桩顶弯矩（69.7kN

的2.4倍。前桩最大弯矩截面高于后桩，即前排桩桩身剪力零点位置高于后桩。

表3为桩土影响范围结果汇总表，可以看出：①相同水平荷载时随着桩间距的增大，群桩侧向影响范围逐渐减小；②沿着荷载方向，承台前土体压缩变形范围明显小于承台后受拉变形范围，后者约为前者的1.3

2015年第8期宋修广等：群桩水平承载特性数值模拟研究

99

图3前、后桩桩身弯矩对比

3.3

桩间距对群桩承载特性的影响分析

图4为不同桩间距对应桩身位移对比图。分析图为3倍桩径时（5.38mm）的76%，表明随着桩间距的增大，群桩效应逐渐减弱，群桩的水平承载力逐渐提高。本文在数值计算的过程中，由于不考虑承台变形，故相同水平荷载作用下前、后桩桩顶处位移一致。

4可知：相同水平荷载作用下，桩顶位移随桩距的增大

不断减小；桩距为9倍桩径时桩顶位移（4.12mm）仅

图41×2型群桩不同桩距下桩身位移对比

图5为不同桩距下桩身弯矩对比图。由图5可知：①相同水平荷载作用下，随着桩间距的增加，桩顶9D桩间距下前、位置处弯矩逐渐增大，后桩桩顶0.5m处弯矩分别为3D桩间距的1.78和5.24倍，桩身最大弯矩随桩间距的增加逐渐减小；②由图5（c）可知，相同水平荷载作用下前桩弯矩明显大于后桩弯矩；前、

后桩上部桩身弯矩变化规律基本一致，首先随桩距的增大而明显增大，变化速率逐渐趋向平缓，当桩距达到6～8倍桩径时，5（b）中弯矩基本保持不变。对比图5（a），7D和9D桩身弯矩曲线，两者基本重合，表明两者承载特性基本一致。随着桩间距的增大，群桩效应影响逐渐减弱，当桩间距大于7D时，可不考虑群桩效应的影响。

图5不同桩距下桩身弯矩对比

4结论

本文通过数值计算对水平荷载作用下1×2型群

揭示了群桩前后桩影桩基础受力特性进行分析研究，

响范围和受力不均匀性以及桩间距对承载特性的影响规律，得出主要结论如下：

100

铁道建筑437．

August，2015

1）群桩承受水平荷载时，承台前土体压缩变形范

围明显小于承台后受拉变形范围；随着桩间距的增加，1×2型群桩侧向影响范围逐渐减小，而前后桩深度影

前桩的深度影响范围小于后桩，但两响范围逐渐增大，

者相差不大。

2）相同水平荷载作用下前、后桩荷载分布不均匀，前桩最大弯矩值明显大于后桩，且前桩最大弯矩截面高于后桩，即前排桩桩身剪力零点位置高于后桩。3）相同水平荷载作用下，随着桩间距的增加，桩

9D桩间距下前、顶处弯矩逐渐增大，后桩桩顶0.5m位置处弯矩分别为3D桩间距的1.78和5.24倍；0.5

m处桩身弯矩随桩间距增加变化速率逐渐趋向平缓，当桩间距达到6～8倍桩径时，弯矩基本保持不变，群桩效应影响减弱。

参

考

文

献

［2］陶学俊，丁选明．水平荷载作用下PCC桩群桩效应数值分

J］．防震减灾工程，2010（30）：287-292．析［

［3］．岩土力学，2005汤斌，陈晓平．群桩效应有限元分析［J］

（26）：299-302．

［4］．水运工蔡亮，鲁子爱，韩洁．横向荷载群桩效应分析［J］

2003（1）：4-7．程，

［5］周洪波，杨敏，茜平一．水平荷载作用下群桩相互作用的弹

J］．水文地质工程地质，2003（3）：29-34．塑性数值分析［

［6］汤斌，王月香，雷学文．承台底土阻力群桩效应系数有限元

J］．岩土力学与工程学报，2004（5）：804-808．分析［

［7］庄培芝，亓乐，周志东．水平荷载作用下单桩承载力影响因

J］．公路，2013（8）：12-18．素分析［

［8］J］．山东科马志涛．混凝土管桩水平受力性状有限元分析［

2009（5）：30-34．技大学学报（自然科学版），

［9］．扬张华华．水平荷载桩的现场试验与数值模拟分析［D］

2007．州：扬州大学，

［10］中华人民共和国交通部．JGJ94—2008

S］．北京：人民交通出版社，2008．范［

建筑桩基技术规

［1］孟凡丽，卢成原，汤丽君．水平荷载作用下群桩—土—结构

．浙江工业大学学报，2005（33）：434-共同作用分析［J］

Numericalsimulationstudyonhorizontalload-bearing

characteristicsofgroupofpiles

22222

SONGXiuguang1，，YUEHongya1，，ZHOUZhidong1，，ZHUANGPeizhi1，，ZHANGHongbo1，

（1．SchoolofCivilEngineering，ShandongUniversity，JinanShandong250061，China；2．ShandongEngineering＆

TechnologyＲesearchCenterforSubgradeSafety，JinanShandong250061，China）

Abstract：Ascurrentstudiesonhorizontalload-bearingcharacteristicsofgroupofpilescannotyetmeettheengineeringreality，thepapertakesthefoundationofpipepilegroupthatlocatesonsoftsoilbaseasthestudyobjectandappliesnumericalsimulationinitsresearch．Theresultsindicatethatashorizontalloadingimposes，soilinfrontofthebearingplatformdeformswithinanevidentlymorelimitedscopethanthatbehindtheplatform．Asthepileintervalgoesup，thelateralscopeofinfluencefallsgradually，whiletheverticalscopeofinfluenceandthebendingthemaximumbendingmomentofthefrontpilelargelymomentatpiletiprise．Underequalhorizontalloading，

exceedsthatofthebackpile，asthenullpointofshearforcelocateshigheronthefrontpilethanthatofthebackpile．Asthepileintervalraisesto6～8timesofthepilediameter，groupofpileseffectisevidentlyweakened．Keywords：Groupofpiles；Horizontalloading；Pipepile；Load-bearingcharacteristics；Influencescope

（责任审编

赵其文）