不同桩头约束下微倾单桩纵横向受荷响应计算的三参数法

张磊+龚晓南

摘 要：为研究桩头转动约束及桩身初始微倾斜对纵横向组合荷载作用下桩身侧向响应的影响，基于三参数形式的地基水平抗力系数，通过矩阵运算提出了桩身变形和内力的半解析解，并与模型试验结果及已有解计算结果进行对比以验证其可靠性。计算结果表明：桩头转动刚度增加时，桩顶位移和地表以下桩身最大弯矩减小，桩顶弯矩和地表以下桩身最大弯矩距离桩顶的距离增大。桩身初始倾角增加时，桩身最大位移和最大弯矩均线性增大，且随纵向荷载的增加其变化速率逐渐增大；纵向荷载增加时，桩身最大位移和最大弯矩均增大，且随纵向荷载和桩身初始倾角的增加其变化速率逐渐增大，而地表以下桩身最大弯矩距离地表的距离呈线性减小。

关键词：桩基础；转动刚度；倾斜；横向荷载；半解析解

中图分类号：TU473.1

文献标志码：A 文章编号：1674-4764（2017）05-0023-08

Abstract：In order to study the influences of the rotational restraint at pile head and the slight initial inclination of pile shaft on lateral pile responses under simultaneous vertical and lateral loads， semi-analytical solutions were put forward by matrix calculation， based on the coefficient of subgrade reaction expressed in the three-parameter form. The results were validated by the comparison with model test results and the existing solutions. The results indicated that when the coefficient of rotational restraint at pile head increases， the displacement at pile head and the maximum bending moment under ground decrease， while the bending moment at pile head and the distance between the maximum bending moment under ground and pile head increased. Both the maximum displacement and bending moment increased linearly with the increase of the initial inclination angle of pile shaft， and their changing rates increased when the vertical load increased. When the vertical load increased， the maximum displacement and bending moment increased， and their changing rates increased with increases of the vertical load and the initial inclination angle of pile shaft. The distance between the maximum bending moment under ground and ground decreased linearly when the vertical load increases.

Keywords：pile foundations； coefficient of rotational restraint； inclination； lateral loads； semi-analytical solutions

樁头转动约束对桩基水平承载性状有显著影响[1]：相同水平荷载下，桩顶固定时桩头水平位移约为桩顶自由时的25%；桩顶自由时桩身最大弯矩位于地表以下某点处，桩顶固定时位于桩顶。现有的桩基水平承载性状研究一般把桩顶的边界条件取为自由和固定两种理想化的情形[2-4]。对于存在上部承台的桩基，当桩头嵌入承台的深度不足或承台配筋量较少时，桩顶附近承台中的混凝土变形或开裂，导致桩头转动。即使桩与承台锚固牢靠，荷载作用下承台的旋转仍会带动桩头转动。因此，承台对桩头的约束弱于桩顶固定并强于桩顶自由。Mokwa等[1]和Prakash等[5]分别采用桩头转动刚度（桩头弯矩与桩头转角的比值）和桩头约束度（桩头弯矩与桩顶固定时桩头弯矩的比值）定量表示桩头转动约束。孔令刚等[6]讨论了几种桩头约束定量表征参数之间的关系，分析了桩身抗弯刚度、水平地基比例系数、自由段及地表以下桩身长度等因素对各表征参数的影响。姜丽红等[7]基于经验公式提出了桩顶位移、桩顶弯矩、地表以下桩身最大弯矩及其所在位置的简化计算方法。不仅横向荷载可以使桩身产生较大的侧向变形和内力，由上部结构及桩身自重等引起的纵向荷载也将因桩身的挠曲变形而产生附加的桩身弯矩和变形[8]。张玲等[9]、Zhang等[10]、Hirai[11]等针对不同工况，分别提出了纵横向组合荷载作用下桩身侧向变形和内力的计算方法，并分析了纵向荷载对桩身响应的影响。然而，已有的考虑桩头转动约束的桩基水平承载性状研究没有考虑纵向荷载的影响。

在实际工程中，桩身初始微倾斜普遍存在并可由多种原因引起[12-13]：压入或打入土体的预制桩，垂直度控制不好导致桩身倾斜；灌注桩施工时，桩孔倾斜导致桩身倾斜。对于纵横向受荷的桩基，由纵向荷载引起的桩身附加弯矩和附加变形可因桩身初始微倾斜而进一步加大。基于室内模型试验，赵明华等[14]研究了纵横向组合荷载作用下微倾木桩的水平承载性状；李微哲等[15]提出了层状地基中地基水平抗力系数随深度线性增加时桩身侧向响应的幂级数解。endprint

本文基于更具一般性的三参数形式的地基水平抗力系数，采用桩头转动刚度表征桩头转动约束的程度，推导出纵横向组合荷载作用下微倾单桩侧向变形和内力的半解析解。编制了计算程序，通过与模型试验结果及已有解计算结果的对比，验证了所得解的可靠性。分析了桩头转动约束、桩身初始微倾斜、纵向荷载等对桩身响应的影响。

1 控制方程与解

1.1 建立方程

桩顶弯矩等于桩头转动刚度与桩顶转角的乘积。至此，桩顶位移、转角、弯矩和剪切力全部得到。由式（6）～（8）及各分界点处桩身响应的连续性，可得任意点处桩身变形和内力。

2 算例验证

采用Fortran语言编制了桩身响应计算程序，现验证解及程序的可靠性。通过室内模型试验，赵明华等[14]研究了微倾木桩在纵横向组合荷载作用下的水平承载性状。第2组第5号试桩主要参数：地表以上桩长为66.0 cm，总桩长为135.0 cm，桩身直径为31.57 mm，抗弯刚度为407.686 6 N·m2，桩身初始倾角θ=0.012 9 rad，桩顶自由，即Km=0 kN·m/rad，桩底自由。上覆黏土层厚度为30 cm，下卧砂土层厚度为70 cm。地基水平抗力系数中的参数：m=2.1×105 kN/m3.8，z0=0.1 m，n=0.8。桩、土其余参数详见文献[14]。地表以上桩身无分布荷载作用，水平力Qp=26.95 N。地表处桩身水平位移与纵向荷载的实测关系及不计桩侧摩阻力和桩身重度影响时，本文解和文献[15]解的计算结果如图3所示。

由图3可见，本文解计算值与模型试验实测结果吻合度较高，所得解及程序是可靠的。另外，本文地基水平抗力系数比文献[15]中的更具一般性，计算值也比文献[15]解计算值更接近实测结果，说明本文解更先进。

3 影响因素分析

为探讨桩头转动约束、桩身初始微倾斜及纵向荷载等对纵横向组合荷载作用下桩身侧向响应的影响，以某桥梁桩基为例进行分析。地表以上及以下樁长分别为15、30 m，桩径为1.8 m，桩底自由，桩身抗弯刚度为9.275×106 kN·m2；比例系数m=6×103 kN/m3.8，z0=0.4 m，n=0.8。地表以上桩身无分布荷载作用，且该段桩身轴力增长系数f0=62.345 kN/m；为方便计算，假定桩身轴力自地表至桩底线性减小到零。

3.1 桩头转动约束的影响

分别利用桩顶自由时桩顶位移ufree、桩顶固定时桩顶弯矩Mfixed及总桩长H对桩顶位移up、桩顶弯矩Mp、地表以下桩身最大弯矩Mmax及其距离桩顶的距离zmax归一化，得到无量纲的量：桩顶位移比up/ufree、桩顶弯矩比Mp/Mfixed、地表以下桩身最大弯矩比Mmax/Mfixed及地表以下桩身最大弯矩距离桩顶的距离比zmax/H。纵向荷载Np=10 MN、桩身初始倾角θ=5×10-3rad，不同的水平荷载下，up/ufree、Mp/Mfixed、Mmax/Mfixed及zmax/H与桩头转动刚度Km的关系如图4所示。

由图4可见，水平荷载的大小对up/ufree、Mp/Mfixed、Mmax/Mfixed及zmax/H的影响很小，可忽略不计。当Km<1×104 kN·m/rad时，近于桩顶自由状态，桩顶弯矩接近零；桩身最大弯矩位于地表以下某点处，其值约为桩顶固定时桩顶弯矩的1.8倍。当Km的值在1×104 kN·m/rad和1×108 kN·m/rad之间时，随着桩头转动刚度的增加，桩顶位移和地表以下桩身最大弯矩迅速减小，而桩顶弯矩和地表以下桩身最大弯矩距离桩顶的距离迅速增大。桩头转动刚度在1×106 kN·m/rad附近存在一临界值，在该点处桩顶弯矩等于地表以下桩身最大弯矩，且桩身最大弯矩达到最小值；当桩头转动刚度大于该临界值时，桩顶弯矩超过地表以下桩身最大弯矩而成为最大弯矩。当Km>1×108 kN·m/rad时，近于桩顶固定状态；桩顶位移约为桩顶自由时的22%，这与文献[1]的结论相近。

3.2 桩身初始微倾斜的影响

分别利用θ=0 rad时的桩顶位移uvp和桩身最大弯矩Mvmax及总桩长H对桩顶位移up、桩身最大弯矩Mmax及其距离桩顶的距离zmax归一化，得到无量纲的量：桩顶位移比up/uvp、最大弯矩比Mmax/Mvmax及桩身最大弯矩距离桩顶的距离比zmax/H。纵向荷载Np=10 MN、水平荷载Qp=250 kN，不同的桩头转动刚度下，up/uvp、Mmax/Mvmax及zmax/H与桩身初始倾角的关系如图5所示。

由图5（a）、（b）可见，随着桩身初始倾角的增加，桩身最大位移和最大弯矩均线性增大，且在桩头转动刚度增加时其变化速率保持不变。桩身初始倾角从0 rad增至1×10-2 rad时，桩身最大位移和最大弯矩均增加约42%。由图5（c）可见， Km=1×103 kN·m/rad及Km=5×105 kN·m/rad时，桩身最大弯矩位于地表以下某点处，且其位置随桩身初始倾角的增加而保持不变。另外，计算结果表明当Km=1×108 kN·m/rad时，桩身最大弯矩位于桩顶。因此，纵横向组合荷载作用下桩身侧向响应受桩身初始微倾斜的影响较大，不容忽视。

3.3 纵向荷载的影响

桩头转动刚度Km=5×104 kN·m/rad、水平力Qp=350 kN；桩身最大位移、最大弯矩以及最大弯矩距离地表的距离在不同的桩身初始倾角下与纵向荷载的关系如图6所示。

由图6（a）、（b）可见，随着纵向荷载的增加，桩身最大位移和桩身最大弯矩均增大，当纵向荷载较小时其变化速率较小。随着纵向荷载和桩身初始倾角的增加，其变化速率逐渐增大。在纵向荷载施加的后期，随着纵向荷载的增加，桩身最大位移和桩身最大弯矩急剧增大，此时，桩基处于失稳状态。桩身最大位移和桩身最大弯矩受桩身初始微倾斜的影响随纵向荷载的增加而增大。由图5（c）可见，随着纵向荷载的增加，桩身最大弯矩距离地表的距离呈线性减小。endprint

4 结论

基于三参数形式的地基水平抗力系数，采用桩头转动刚度表征桩头转动约束的程度，推导出纵横向组合荷载作用下微倾单桩侧向变形和内力的半解析解。计算分析后得到以下结论：

1）随着桩头转动刚度的增加，桩顶位移和地表以下桩身最大弯矩均减小，桩顶弯矩和地表以下桩身最大弯矩距离桩顶的距离均增大。桩头转动刚度存在一临界值，在该点处桩顶弯矩等于地表以下桩身最大弯矩；当桩头转动刚度大于该临界值时，桩顶弯矩超过地表以下桩身最大弯矩而成为最大弯矩。

2）随着桩身初始倾角的增加，桩身最大位移和最大弯矩均线性增大，且在桩头转动刚度增加时其变化速率保持不变。桩身最大位移和最大弯矩受桩身初始微倾斜的影响随纵向荷载的增加而增大。桩身最大弯矩所在位置随桩身初始倾角的增加而保持不变。

3）随着纵向荷载的增加，桩身最大位移和桩身最大弯矩均增大，其变化速率也随纵向荷载及桩身初始倾角的增加而逐渐增大，而地表以下桩身最大弯矩距离地表的距离呈线性减小。

参考文献：

[1] MOKWA R L，DUNCAN J M.Rotational restraint of pile caps during lateral loading [J].Journal of Geotechnical and Geoenvironmental Engineering， ASCE，2003，129（9）：829-837.

[2] FAYYAZI M S，TAIEBAT M，FINN W D L.Group reduction factors for analysis of laterally loaded pile groups [J].Canadian Geotechnical Journal，2014，51（7）：758-769.

[3] LAORA R D，ROVITHIS E.Kinematic bending of fixed-head piles in nonhomogeneous soil [J].Journal of Geotechnical and Geoenvironmental Engineering，ASCE，2015，141（4）：04014126.

[4] TEHRANI F S，PREZZI M，SALGADO R.A multidirectional semi-analytical method for analysis of laterally loaded pile groups in multi-layered elastic strata [J].International Journal for Numerical and Analytical Methods in Geomechanics，2016，40（12）：1730-1757.

[5] PRAKASH S，SHARMA H D.Pile foundations in engineering practice [M].New York：Wiley，1990：252-256.

[6] 孔令刚，姜丽红，陈仁朋，等.考虑桩头约束的水平受荷单桩响应分析[J].应用力学学报，2014，31（3）：393-399.

KONG L G，JIANG L H，CHEN R P，et al.Response of single pile with partial pile head fixity subjected to lateral loading [J].Chinese Journal of Applied Mechanics，2014，31（3）：393-399. （in Chinese）

[7] 姜丽红，孔令刚，蔺港，等.不同桩头约束下水平受荷单桩简化计算方法[J].工程力学，2014，31（11）：63-69.

JIANG L H，KONG L G，LIN G，et al.Simplified method of laterally loaded single piles with partial pile head fixity [J].Engineering Mechanics，2014，31（11）：63-69. （in Chinese）

[8] LIANG F Y，CHEN H B，CHEN S L.Influences of axial load on the lateral response of single pile with integral equation method [J].International Journal for Numerical and Analytical Methods in Geomechanics，2012，36（16）：1831-1845.

[9] 張玲，赵明华，赵衡.倾斜荷载下桩柱式桥墩受力变形分析传递矩阵法[J].中国公路学报，2015，28（2）：69-76.

ZHANG L，ZHAO M H，ZHAO H.Transfer matrix method for deformation of pile type bridge pier under axial transverse load [J].China Journal of Highway and Transport，2015，28（2）：69-76. （in Chinese）

[10] ZHANG L，GONG X N，YU J L.Elastic solutions for partially embedded single piles subjected to simultaneous axial and lateral loading [J].Journal of Central South University，2014，21（11）：4330-4337.endprint

[11] HIRAI H.A winkler model approach for vertically and laterally loaded piles in nonhomogeneous soil [J].International Journal for Numerical and Analytical Methods in Geomechanics，2012，36（17）：1869-1897.

[12] 鄭刚，李帅，杜一鸣，等.竖向荷载作用下倾斜桩的承载力特性[J].天津大学学报，2012，45（7）：567-576.

ZHENG G，LI S，DU Y M，et al.Bearing capacity behaviors of inclined pile under vertical load [J].Journal of Tianjin University，2012，45（7）：567-576. （in Chinese）

[13] NADEEM M，CHAKRABORTY T，MATSAGAR V.Nonlinear buckling analysis of slender piles with geometric imperfections [J].Journal of Geotechnical and Geoenvironmental Engineering，ASCE，2015，141（1）：06014014.

[14] 赵明华，李微哲，杨明辉，等.成层地基中倾斜偏心荷载下基桩位移特性室内模型试验研究[J].土木工程学报，2006，39（12）：95-99.

ZHAO M H，LI W Z，YANG M H，et al.A model test study on displacement of piles under inclined and eccentric loads in layered soils [J].China Civil Engineering Journal，2006，39（12）：95-99. （in Chinese）

[15] 李微哲，赵明华，单远铭，等.倾斜偏心荷载下基桩内力位移分析[J].中南公路工程，2005，30（3）：53-57.

LI W Z，ZHAO M H，SHAN Y M，et al.Analysis of single pile under eccentric and inclined loading [J].Central South Highway Engineering，2005，30（3）：53-57. （in Chinese）

[16] 张磊，龚晓南，俞建霖.纵横荷载下单桩地基反力法的半解析解[J].哈尔滨工业大学学报，2011，43（6）：96-100.

ZHANG L，GONG X N，YU J L.Semi-analytical solutions for single piles subjected to combined vertical and lateral loads by subgrade reaction method [J].Journal of Harbin Institute of Technology，2011，43（6）：96-100. （in Chinese）

（编辑 王秀玲）endprint