多因素耦合下深水钢管桩偏位简化计算方法

张刚领，陈永亮，崔福刚

（1. 重庆大学 土木工程学院,重庆 400045； 2. 中铁十一局集团第五工程有限公司,重庆 400037）

为方便跨江、跨河桥梁主体结构施工作业，栈桥已成为不可缺少的临时性结构物，其多采用贝雷梁体系作为上部结构，钢管桩作为下部结构[1-4]。近年来施工栈桥质量、安全事故频发，例如2018年新建福州至平潭铁路站前工程因栈桥的下部结构施工延期，造成大桥整体工期延长6个月；2017年益娄衡高速大桥栈桥下部结构垮塌，造成较大的经济损失。这些事故不但严重影响施工进度，更造成恶劣社会影响和巨大生命财产损失。施工栈桥是跨江桥梁施工中的主要辅助措施，而施工栈桥的关键在于下部钢管桩结构的安全性。

目前，对于桩长较小的钢管桩主要通过有限元分析整体受力情况，例如文献[5-8]结合桥位地质、水文条件，在探究钢栈桥设计思路及施工工艺的基础上，考虑栈桥自重、风荷载、流水压力及其他可变荷载等工况，对钢管桩各杆件承载力、刚度和稳定性进行验算。从受力性质角度来看,栈桥钢管桩属于典型的压杆构件，当压杆长度较长时，较多研究表明轴力将使构件对初始几何缺陷的敏感性显著增大[9-12]；同时，较为精确的二阶分析方法在考虑几何非线性后，计算所得构件内力及位移比传统一阶分析结果更大[13-15]。因此，当钢管桩较长时，应综合考虑施工竖向荷载、水流水平向荷载、几何初始缺陷和几何非线性的耦合作用。

虽然目前可通过精细化的有限元模型计算得到耦合作用下钢管桩的内力及位移，但在实际工程中不便于广大工程人员应用，故往往采用增大安全系数的方法来保障设计安全。针对这一问题，基于弹性稳定理论，在同时考虑施工荷载、水流荷载以及几何非线性、几何初始缺陷耦合作用下，通过理论推导得出钢管桩桩顶偏位的近似计算公式，为栈桥钢管桩的偏位计算提供新方法。

1 桩顶偏位规范简化计算方法

1.1 栈桥钢管桩简化分析模型

栈桥钢管桩往往单排设置2根桩，桩底部通常嵌入基岩并浇筑混凝土；2根钢管桩之间采用横向连系连接成整体；桩顶设置横梁，在横梁上再布置贝雷梁。因此，单排钢管柱可简化为高度为l的单根悬臂柱。

施工荷载可视化为作用在墩顶的集中力N和弯矩M。几何初始缺陷选用沿桩高具有连续几何分布的曲线y0表示，墩顶处的偏位值为Δ。记钢管桩顶部为A，底部为C，水流高度处为B，如图1所示。

图1 水流荷载简化模型Fig. 1 Simplified model of flow load

对常见等截面圆形钢管桩，水流荷载可按《公路桥涵通用设计规范》（JTG D60—2015）及《铁路桥涵设计规范》（J460/TB10002—2017）规定，水深h处的水流速度V及流水压力标准值Fw可分别取为

式中，Vmax为水面上的流速，H为总水深；n为常数，由水流性质、河床性质等多种因素确定，通常取为1/4~1/10；γ为水重力密度；A为阻水面积；g为重力加速度；K为桥墩形状系数。

当n取1/2作为一般河流水流流速的分布曲线已偏于安全，同时将式（1）带入式（2）可得：

由式（3）可以看出，水流荷载沿河流垂直面线性增大，为准确计算水流荷载的作用，可将水流荷载沿钢管桩的分布简化为底部为0、上部为q的倒三角形。倒三角水流荷载作用下钢管桩偏位模型以及倒三角水流荷载的计算模型沿钢管桩的分布如图1所示。

1.2 规范计算公式

按现行规范规定，假定水流压力的合力点在设计水位线以下0.3倍水深处，F为流水压力标准值，为集中力，Δ为桩顶偏位值，基于力法原理，可得到水流荷载及几何初始缺陷共同作用下钢管桩桩顶的偏位w1为：

1.3 耦合作用下桩顶偏位计算方法

为方便计算，根据水位高度将钢管桩分为上、下两部分，直接承受水流荷载的下部记为第I部分；其余部分记为第Ⅱ部分，如图1所示。

当考虑钢管桩几何非线性时，水流荷载及钢管桩几何初始缺陷产生的侧移无法直接叠加，须同时考虑两者耦合作用。为方便计算，假设钢管桩顶部整体初始几何缺陷位移函数为：

由式（5）可得钢管桩高度为l1时的钢管桩几何初始缺陷值y0(l1)，记为Δ1。则第I部分和第Ⅱ部分钢管桩的初始几何缺陷y0I和y0Ⅱ可分别表示为：

记整根钢管桩桩顶几何初始缺陷值与高度为l1处几何初始缺陷的差值为e，由力平移定理可知，竖向集中力在钢管桩l1高度处形成的附加弯矩MB为：

取钢管桩第I部分的上部l1-x1节段为隔离体，建立力平衡微分方程：

式中：wI为同时考虑钢管桩几何非线性及水流荷载下桩顶侧向位移；yI为l1-x1截面处侧向位移。式(8)的通解为：

式中，C1、C2为常数。利用带入式(9)中，并令β=N/NE，α=l1/l，可得第I部分钢管桩顶部偏位wI为：

为简化计算，采用三角函数简化式(10)，可得：

对式（10）求导并带入边界条件，可得第I部分钢管桩顶部转角θI：

采用三角函数简化式(12)，可得：

取钢管桩上部l2-x2节段为隔离体，并建立力平衡微分方程：

式中：wⅡ为考虑钢管桩几何非线性及水流荷载下桩顶侧向位移；yⅡ为l2-x2截面处侧向位移；MA=N×wI，为第I部分钢管桩偏位引起第Ⅱ部分钢管桩的横向弯矩。

与第I部分钢管桩偏位求解过程相同，通过方程求解、三角函数化简可得wⅡ的表达式为：

2 钢管桩垂直度偏位计算公式简化

通过上述分析可以发现，考虑几何非线性的钢管桩垂直度偏位计算结果较不考虑几何非线性的钢管桩垂直度偏位较大，更符合实际情况，但是倒三角水流荷载计算公式较为复杂，在实际工程中应用不太方便，因此应对多种因素作用下倒三角水流荷载的钢管桩垂直度偏位计算公式进行简化。

钢栈桥在设计时为满足施工要求，钢管桩高度要始终大于水位高度，为简化计算，取钢管桩上部不受水流荷载作用部分，即第Ⅱ部分，由式（15）可得：

由式（17）可以看出，由几何初始缺陷引起的钢管桩顶偏位随α值变化而变化，因此可取f，如式（18）所示：

通常情况下，水位高度大于0.5倍钢管桩高度，因此α的取值范围为0.5~1，增量为0.01，可得f的取值如图2所示。

图2 不同α取值下f值Fig. 2 f at different values of α

由图2可知，f的取值与α成反比关系，当α>0.65即水位高度与钢管桩高度之比大于0.65时，f的取值小于0.05，同时当α>0.85，f的取值趋近于零，可忽略不计，所以为简化计算，当α>0.65时，式（17）可简化为：

根据式（19）可得，当α>0.65时，可将式（16）简化为：

3 简化公式的准确性

由第2节可知，式（20）主要针对式（17）进行了简化，由式（17）可以看出，f的取值只对几何初始缺陷引起的偏位有影响。为探究简化式（20）的准确性，考虑具有普适性的栈桥深水钢管桩特征参数，选取高度为60 m的钢管桩，水流速度取0.5 m/s，轴压比取0.05~0.2，增量为0.05；设水位高度与钢管桩高度比值α为0.65~0.95，增量为0.1。钢管桩外径为1 m，壁厚为12 mm，采用Q235钢材，容重为7 850 kN/m3，弹性模量取为2×105MPa。钢管桩在运输和插打过程中的初始缺陷根据《建筑施工临时支撑结构技术规范》取值，即考虑钢管桩初始垂直度偏差范围为20~100 mm，增量为20 mm。由式（17）可以看出，wⅡ中几何初始缺陷引起的钢管桩顶偏位与α、几何初始缺陷取值以及轴压比有关，通过对比不同α取值及不同几何初始缺陷在不同轴压比下简化公式（20）与原公式（16）钢管桩垂直度偏位差值如图3~4所示。

图3 不同α取值下简化公式垂直度误差Fig. 3 The error of perpendicularity of the simplified formula under different values of α

由图3可以看出，随着α取值的提高，简化公式（20）与原公式（16）垂直度偏位计算值的差值略有提高，但仍控制在5%以内，可以忽略不计。由图4可知，随着钢管桩几何初始缺陷取值的增加，简化公式（20）与原公式（16）垂直度偏位计算值的差值略有提高，但控制在3%以内，可以忽略不计。

图4 不同几何初始缺陷下简化公式垂直度误差Fig. 4 Deviation value of perpendicularity of steel pipe piles of different heights

由上述计算结果可知，简化计算公式（20）和原计算公式（16）垂直度偏位计算值差值随着轴压比提高而略有增加，但基本一致，因此，在正常施工应用中，可按照简化计算公式（20）进行钢管桩垂直度偏位值计算。

4 结 论

通过倒三角水流荷载简化模型，根据力平衡微分方程推导钢管桩在考虑几何非线性条件下垂直度偏位计算公式，并根据规范提供的方法进行优化，主要结论如下：

1）钢管桩垂直度偏位受钢管桩高度和水流荷载作用影响很大，随着钢管桩高度和水流速度的增加，钢管桩偏位呈非线性上升趋势。

2）结果表明，钢管桩受几何非线性影响显著，且随着轴压比的增加，几何非线性的影响更加明显，因此在钢管桩垂直度偏位计算时必须考虑P-Δ效应的影响。

3）在实际工程中，考虑几何非线性的钢管桩垂直度偏位计算公式更精确、更符合实际情况，但是倒三角水流荷载计算公式较为复杂，在实际工程中应用不太方便，因此针对多种因素作用下倒三角水流荷载的钢管桩垂直度偏位计算公式进行优化。

4）上述分析可以看出轴压比对钢管桩垂直度偏位值的影响很大，因此施工过程中，在超深水期以及三峡水库泄洪时应严格控制施工荷载以及监测钢管桩垂直度偏位的变化，降低几何非线性对钢管桩垂直度偏位值的影响程度，保障施工安全。