船行波激励下施工栈桥动力响应及控制研究

2022-06-21陈永亮

重庆大学学报 2022年5期

陈永亮，刘纲

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

施工栈桥作为跨江、跨海桥梁主体施工期间重要的临时辅助设施，其结构安全直接影响桥梁主体的正常施工和安全。实践表明，往来栈桥附近的船舶引起的波浪在一定条件下会迫使栈桥产生较大振动，影响施工安全甚至造成栈桥损伤，埋下一定安全隐患。为保证栈桥安全，亟待展开船行波下栈桥动力响应特性分析，制定具有针对性的控制措施。在船行波理论及模拟方面，王新宇等[1]基于三维势流理论，针对无限水深情况，通过Rankine源高阶面元法求解非线性船行波的波高；Abbasnia等[2]基于势能数值拖曳槽和非均匀有理B样条曲线进行全非线性和线性船舶波浪的建模[2]；在船行波对船只及结构作用求解方面，主要围绕CFD数值模拟或试验展开研究[3-9]，例如，石博文等[3]采用CFD方法分析了船舶在骑浪状态下的稳性损失；赵文宾等[4]借助FLOW-3D三维数值模拟软件模拟渡槽建成后的船行波、流场分布情况，对渡槽内船行波分布规律及消能设施效果展开研究；徐博等[5]基于RANS方程，模拟分析波浪对跨海桥梁高桩承台的冲击作用；李志松等[6]基于RANS方程、RNGk-ε湍流模型建立数值波浪水池，用VOF法追踪自由面的位置，并模拟内河航道中船行波的生成、传播、爬高和岸壁反射过程，对数值模型的有效性进行验证；姜辉等[8]通过实测船行波的特征参数，结合数值模拟波浪水槽仿真桥墩所受波浪力，分析了桥墩在波浪下的动力特性。上述方法具有较高的计算精度，但计算较复杂且耗时较多，特别是在水位大幅波动情况下，需根据不同工况建立不同模型，不便于工程实际应用。目前，专门针对船行波作用下施工栈桥动力响应分析的研究较少，文中以郑万高铁某特大桥施工栈桥为例，按照线性波理论简化波浪荷载，采用附加质量与附加阻尼比简化考虑流固耦合效应，在分析不同船速、船型、水深情况下，船行波对施工栈桥动力响应的影响，根据响应特性制定相应的处治措施，为保障栈桥结构安全提供技术参考。

1 工程概况

郑万高铁某特大桥的施工栈桥主体部分跨径组合为6 m+3 m×9 m+6 m，其支桥跨径组合为6×6.3 m，桥面宽9.0 m。栈桥下部结构为φ1 000 mm×12 mm钢管桩，钢管桩最长为65.0 m，入岩最大深度为5 m。单排设置2根钢管桩，钢管桩横向间距6.3 m，并采用φ630 mm×10 mm钢管连成整体，桩顶设置双拼I40b作为横梁。上部结构采用贝雷梁，布于横梁上，横向间距90.0 cm、45.0 cm，梁上铺设I20b分配梁，间距352.5 mm，其上铺设10 mm厚花纹钢板作为桥面。栈桥布置如图1所示。

2 船行波荷载确定

为研究船行波对施工栈桥动力特性的影响，需明确船行波荷载大小，进而分析栈桥的动力响应。通过调研及施工现场的长期观察，发现该栈桥所处河流段内的船型相差不大，主要为客船和驳船，故仅针对客船和驳船两种船型作用下船行波荷载进行分析，为简化计算，采用线性波理论确定波浪荷载。

2.1 船行波波高

船行波荷载的大小与船行波波高有较大关系。根据文献[10-11]的研究以及文献[14]的实际测量，在内河深水条件下，库斯科公式所得船行波参数与实际测量船行波参数有较好的吻合度，即船行波波高Hm为

(1)

式中：v为船舶与水的相对速度；B′为船舶离岸距离；g为重力加速度；ρ为水的密度；R为船舶兴波阻尼；V为船舶排水量；L为船舶长度；d为航道水深。

2.2 船行波波长与周期

船行波的波长和周期对船行波荷载有较大影响，在确定船行波荷载之前需对波长和周期进行计算。船行波分为散波和横波，在深水条件下，影响结构动力特性的主要是散波，其扩散方向不受航速影响，在长距扩散过程中，能量衰减较少，且移动方向总是垂直于波峰线[3]。因依托栈桥工程的桩长达60 m，河流水位较深，故主要考虑散波对栈桥结构的影响。

根据OCDI经验公式[12]，在水位较深情况下可得到横波的波长L0：

(2)

式中:T0为横波周期。在此基础上，可得到对应散波的波长Ld和周期Td为

(3)

式中:θ为凯尔文角，文献[14]建议一般情况下取19.467°。

2.3 船行波荷载

根据《港口与航道水文规范》(JTS145—2015)，采用Morison公式计算栈桥钢管桩波浪压力分布[13]。由于钢管桩相互之间的中心距大于4倍桩径，可忽略波浪力之间的相互干扰，无需考虑群桩效应。作用于某一高度柱体全断面上的等效波浪力分布p为

p=pd+pI,

(4)

其中：

(5)

式中:pd为波浪力速度分力；pI为波浪力惯性分力；γ为水重度；DZ为桩体直径；A为桩体断面面积；u为水质点轨道运动水平速度；H为钢管桩处船行波波高；z为计算点在水底面以上的高度；ω为波浪运动圆频率；CD为速度力系数;CM为惯性力系数。以上各系数具体取值可由《港口与航道水文规范》(JTS145—2015)确定。

2.4 客船、驳船船行波荷载计算

经现场统计观测和调研，以栈桥所在河段内行驶客船、驳船最不利尺寸考虑船行波对钢管桩的荷载。客船长度L取16.3 m，满载排水量取11.60 t，最大船速取30 km/h；驳船长度L取65.0 m，满载排水量取987 t，最高船速取12 km/h。鉴于船行波作用影响时间较短，假设所有船只仅从河中央附近通行，每次只有1艘且只有1种船型通过航道。由式(1)～式(3)可得不同船型、不同船速下船行波的计算参数，如表1、表2所示。

表1 客船在不同船速下船行波参数表

表2 驳船在不同船速下船行波参数表

从表1、表2可知，船速越高，所产生的波高越大，波长越长，频率越小；驳船的排水量较大，所产生的波浪作用总体较大。采用式(4)～式(5)，分别计算客船、驳船在不同行驶速度下的船行波荷载峰值如图2、图3所示。以30 km/h速度下客船为例，钢管桩波浪力的速度分力、惯性力分力,如表3所示。

图2 客船在不同船速下行波荷载沿深度变化Fig.2 Variation of Wave Load along Depth for Ships with Different Velocities

图3 驳船在不同船速下船行波荷载沿深度变化Fig.3 Variation of Wave Load along Depth for Barges with Different Velocities

表3 不同水深的船行波荷载值

由图可知，船行波荷载峰值随水深增大而迅速减小；船速越小，荷载峰值变化越剧烈，其荷载峰值及辐射深度减小得越快，荷载总体趋势表现为随着船速的降低而减小。由于排水量大，驳船所产生的船行波的荷载峰值较大。

由表3可知，船行波波浪力的惯性分力pI远大于速度分力pd，说明惯性力分力pI起主要影响作用，根据公式(5)可知，由于船行波的水质点轨道运动的水平速度u对pd的影响更大，而船行波所引起的u却较小，从而导致速度分力pd较小。

3 栈桥结构动力特性分析

3.1 有限元模型的建立

船行波往往仅作用在栈桥下部的钢管桩上，故仅建立钢管桩体系的模型，栈桥上部的横梁、贝雷梁和桥面板简化为质量直接作用在钢管桩顶部。采用梁单元，根据设计图纸建立钢管桩体系有限元模型，模型划分34 239个节点，17 268个单元。因钢管桩嵌入岩层较深，且混凝土锚固桩的混凝土浇筑质量较好，将钢管桩底部近似看作固接；因上部有贝雷梁的限制作用，在钢管桩顶端限制其发生转角位移；钢管桩与横向支撑之间采用固接。建立的有限元模型如图4所示。

图4 栈桥有限元模型Fig. 4 Finite element model of trestle

根据附加质量原理可知，水中结构做变速运动时，受到水体质量的影响，其作用相当于在结构上附加一部分质量，如公式(6)所示：

F=(m+ma)a，

(6)

式中：F为外力；m为结构质量；ma为附加质量；a为结构加速度。

水体对结构的作用将导致结构自振频率下降。而当前流固耦合的数值模拟较为复杂，结构频率会随着水深的变化而变化，工程在栈桥使用期间，河水的水位发生较大变动，流固耦合方法不便于在该类工程中应用。因此，将水对钢管桩的影响考虑为附加于浸水部分桩体上的分布质量。根据文献[16]的研究，对于圆筒结构，其附加质量的经验公式为

(7)

式中：Dc为钢管桩截面外径；Lc为圆筒的长度。系数Cam为1.9，由实验数据得到[16]。

阻尼对结构动力响应的影响较大，因钢管桩入水的深度较大，需同时考虑钢管桩结构阻尼和水对钢管桩的阻尼作用。文献[17]表明，水体引起的阻尼作用可视为等效外阻尼比，与结构内阻尼比线性叠加，从而确定结构在水中的总阻尼比。在计算中，根据舒适度验算要求，栈桥结构内阻尼比取为0.015，水体引起的阻尼比根据文献[17]所给曲线，针对不同水深进行选取，对于圆筒钢管桩，水体引起的阻尼比取值宜为0.003 5～0.007 5。栈桥使用期间，河水水位将发生较大变动，需考虑不同水位深度波浪力作用下栈桥结构动力响应的变化。根据栈桥所处河道情况，主要考虑3种典型水位条件，如表4所示，不同水位高度与栈桥高度的关系如图1所示。为简化计算，将栈桥全桥进行节段划分，并根据第2节推导的荷载计算方法，计算每个单元的波浪荷载。假定相邻节点间的波浪荷载为线性分布，将各单元的波浪荷载分别施加在模型上，再对栈桥动力响应进行分析。结合舒适度与实际施工需求，在船行波作用下栈桥响应最大位移不应超过1 cm，最大加速度不应超过0.1 m/s2。

表4 不同水位工况

3.2 客船动力响应分析

根据所建立的有限元模型及船行波荷载，采用谐响应分析，提取对应频率下栈桥的动力响应进行讨论。分析不同施工水位状态下，不同客船行驶速度产生船行波对栈桥的动力响应，结果表明，栈桥最大位移和最大加速度均发生在桩长最大且靠近河道的钢管柱顶端，如图4中A点所示。该点的动力响应结果如表5所示。

表5 客船船行波作用下施工栈桥动力响应

由表5可知，栈桥顺桥向最大位移与最大加速均出现在船行速度为20 km/h时，对应的船行波频率为0.316 Hz，与结构第一自振频率0.327 Hz较接近，易发生共振现象，故栈桥结构将产生较大的位移或加速度；当船速加快后，船行波频率将远离结构自振频率，结构的动力响应也相应减小。为此，需对船行速度进行限制，避免栈桥结构出现共振。因船速高于15 km/h时，栈桥结构的最大位移超过了1 cm的限值，仅针对船行速度为15 km/h时，不同水位条件下栈桥结构的动力响应进行分析，其结果如表6所示。

表6 不同水位下船行波作用栈桥动力响应

由表可知在最低水位时，结构自振频率为0.383 Hz，波浪力频率为0.421 Hz，两者的频率较接近，此时的响应大于其他水位条件下的动力响应。在最高水位时，结构频率降低，与波浪力频率相差更远，此时，荷载作用位置上移以及受力面积增大，但结构响应增加程度却不大。因此，不同船速导致船行波频率改变对栈桥产生的作用远大于不同水位高度产生的作用。

通过上述分析可得，为保证栈桥结构在主桥施工期间的安全，综合考虑对位移与加速度的限制，应限制客船等小型船只的船行速度在15 km/h以下。

3.3 驳船动力响应分析

分析常规施工水位条件下，不同驳船行驶速度产生的船行波对栈桥动力响应的影响。与客船规律相同，驳船行驶下最大位移和加速度发生在图4中的A点处，其位移和加速度如表7所示。

表7 驳船船行波作用下施工栈桥动力响应

由表可知，当船速增加时，栈桥结构最大位移和最大加速度均增加，且增加速度较快。因此，对此类大型船只船行速度的控制应更加严格。由图2、图3可知，相对于客船，驳船所产生的船行波荷载较大，因此在船速较低的时候也有较大响应。船速在15 km/h以下时，船行波所能影响水位深度较浅，并且船速降低时，对应荷载峰值将大幅减小，影响深度也将进一步减小，在降低船速的时候栈桥的响应迅速减小。结合表7结果，为保证施工安全，驳船的航行速度应限制在10 km/h以下。

4 参数影响分析

4.1 附加质量的影响

为验证文中所采用的附件质量简化方法的准确性，针对文献[16]中的圆柱桩模型试验条件建立考虑附加质量的模型，进行结构频率分析，并与文献[17]所给FSI边界流固耦合模型的计算结果及文献[16]实验所得数据进行对比，如表8所示。

表8 不同方案结构频率对比

由表8计算结果可知，水深越小，附加质量简化方法和流固耦合计算方法所得结果与试验结果越接近。在所有水深条件下，附加质量简化算方法计算结果与试验结果较为接近，误差在10%以内，与精细化的数值耦合计算方法结果的误差小于3.5%，基本满足工程应用精度的实际要求。

采用附加质量法计算不同水位条件下栈桥自振频率如表9所示。

表9 不同工况下结构的自振频率

由表9可知，水位高度对结构自振频率影响较大，随着水位的上升，结构自振频率下降；最高水位与不考虑水体影响情况下自振频率将降低47.2%，必须重视水位高度对结构频率的影响。

4.2 船行波频率的影响

因客船这类中小型船产生波浪力的频率可能与结构自振频率相近，此时结构的动力响应将急剧增加，通过谐响应分析，得到不同频率下栈桥结构最大位移响应如图5所示。

图5 不同波浪频率下栈桥最大位移Fig.5 Maximum Displacement of Trestle under Different Wave Frequencies

由图5可知，当波浪力频率接近结构自振频率时，结构的位移响应急剧增大；而波浪力频率远离结构自振频率时，位移响应均较小。根据结构动力学原理，取栈桥基频的±15%为共振区，当波浪力频率在共振区以外时，由于客船船行波荷载值较小，由船只造成的结构位移响应的变化不大。若船只速度限制在15km/h以内，可保证栈桥结构不发生共振，从而预防栈桥发生过大变形。