风暴潮、浪作用下中小型跨海桥梁破坏形式与数值模拟研究

王玉琦

摘要：强台风形成的巨浪、风暴潮对跨海桥梁的安全造成了严重威胁，为探讨台风风暴潮、浪共同作用下中小型跨海桥梁上部结构和下部结构的破坏机理及台风风暴潮期间不同水位下的波浪对下部结构的作用规律，开展了台风风暴潮、浪共同作用下中小型跨海桥梁破坏形式调研，发现风暴潮、浪作用下中小型跨海桥梁上部结构的破坏形式主要有梁体离位、梁体脱落和梁体断裂，主要原因为梁体被迫发生无约束运动;下部结构破坏主要表现为结构偏移、倒塌以及盖梁连接失效破坏，主要原因为水动力作用较强及上部结构撞击。通过CFD三维数值模拟研究了桥梁下部结构波浪荷载，给出了考虑风暴潮水位下波浪冲击桥墩的经验公式。根据调研及数值模拟，提出了设计及施工建议。

关键词：风暴潮;跨海桥梁;破坏形式;破坏原因;水动力

中图分类号：U447 文献标志码：A 文章编号：2096-6717（2019）06-0152-06

跨海桥梁在沿海交通建设中具有重要作用，是国家经济大动脉上的关键和重大基础设施。截至2016年底，东海大桥、杭州湾跨海大桥、舟山跨海连岛工程、青岛海湾大桥和港珠澳大桥相继建成。大量中小型跨海桥梁也遍布中国东南沿海地区。中国东南沿海地区普遍受到西北太平洋台风生成区和东亚季风区的影响，台风、寒潮等惡劣天气引发的极端波浪灾害频发。台风伴随的强风和气压骤变，除了引发巨浪之外，还会引起海面的剧烈震荡、海面的异常升高或降低，从而产生严重的风暴潮。由于中国沿海地区桥梁结构常年受台风风暴潮、浪的影响，所以探究台风风暴潮、浪共同作用下跨海桥梁的破坏机理对跨海桥梁的设计具有重要意义。现场灾后调研，是研究灾害机理的重要手段。2005年Katrina飓风造成了美国墨西哥湾沿岸数十座高速公路桥梁损坏，Okeil等现场调研了部分桥梁的桥梁损坏形式，为台风风暴潮、浪共同作用桥梁损坏研究提供了宝贵的资料。徐淑颂通过现场调研，研究了海洋环境因素对沿海公路桥梁的影响，并提出了相应的防治对策。

桥梁结构通常由上部结构、支座和下部结构3部分组成，上部结构与下部结构之间采用支座结构传递力学行为。对于小跨径的预制桥梁，当上部结构由于台风及风暴潮造成损坏时，可进行更换，对灾后救援和重建影响较小;但下部结构发生损坏时，灾后工作将受到严重制约，故跨海桥梁下部结构对结构性损坏的要求较上部结构更为严格。风暴潮是强烈的大气扰乱造成的海面异常升高现象（如强风）。台风、风暴潮期间，风暴潮对波浪作用的桥梁下部结构波浪荷载分布有显著影响。邓江源对现有的跨海大桥下部结构波浪力的计算方法进行了比较研究。房忱等研究了跨海大桥桩基、桥墩等小尺度下部结构波浪荷载的非线性特性，研究发现，对于浅水区的桥墩，非线性波的非线性作用较为显著，与采用Morison方程确定的线性波波浪力有14%的差异。由于中小跨径跨海桥梁大多修建在近海浅水区域，波浪、波浪荷载的非线性显著，采用Morison方程可能低估了波浪荷载，同时，已有研究大多局限于研究不同下部结构类型、不同波浪入射方向角、不同波浪条件等下部结构波浪力的计算方法及规律，较少关注不同风暴潮增水情形对下部结构波浪力的影响。因此，研究台风风暴潮期间不同水位下的波浪对下部结构的作用规律，对跨海桥梁下部结构的设计具有重要意义。

笔者调研了历史台风风暴潮、浪共同作用下中小型跨海桥梁破坏情况，分析了台风风暴潮、浪共同作用下跨海桥梁上部结构和下部结构的破坏原因;为获得台风风暴潮期间不同水位下的波浪对下部结构的作用规律，开展了CFD三维数值模拟，探讨了不同风暴增水对波浪作用桥墩的影响，提出了风暴潮、浪作用下中小型跨海桥梁的设计与施工建议。

1桥梁结构破坏形式调研及成因分析

根据美国沿海桥梁在2005年Katrina飓风期间在台风风暴潮、浪的共同作用下桥梁结构发生的结构损伤特征，可将上部结构的破坏形式分为梁体移位、梁体脱落和梁体断裂等;下部结构的破坏形式可分为结构偏移、倒塌及盖梁连接失效等。每种破坏形式的形成原因有所不同。

1.1上部结构梁体移位

梁体移位主要是由于台风引起的风暴潮冲击梁体，使梁体发生了沿横桥向的位移。梁体移位原因如图1所示。由图1可知，对于沿海桥梁而言，梁体主要受到两个方向的力，一个是沿横桥向的冲击力，另一个是沿竖向的浮力。当风暴潮到达梁位，梁体部分或整体被淹没，梁体将受到竖向的浮力，浮力的存在削弱了梁体与支座之间沿横桥向的抗力，使梁体在横桥向的冲击力下更容易被推动，此时，梁体保持完好，但丧失通行能力。

1.2上部结构梁体脱落

梁体脱落指部分梁体被整体完好地推落至海域里，出现脱落现象。出现该种破坏形式的原因与梁体移位类似，主要是由于竖向浮力高于横向冲击力而先使梁体脱离支座约束，再在横向冲击力的作用下出现脱落破坏。

1.3上部结构梁体断裂

沿海桥梁在台风风暴潮、浪的共同作用下出现的梁体断裂主要出现在跨中位置。该破坏形式主要是由于台风引起的海面异常升高和降低，使海面时高时低。当支座不具备抗拉能力或抗拉能力较弱时，梁体在浮力的作用下时而上升、时而降落，在反复循坏作用下损坏。对于简支梁桥而言，跨中弯矩最大，最容易进入极限状态，出现塑性破坏，如图2所示。

桥梁上部结构还包括桥面构造。桥面构造主要包括桥面铺装、排水防水系统、人行道、路缘石、栏杆、护栏、照明灯和伸缩缝等。栏杆、护栏和照明灯阻水面积较大，水动力效应显著，所以，在风暴潮的作用下较容易发生损坏;护栏结构目前很少进行具体计算设计，其考虑最多的因素是车辆的撞击力，护栏设计通常只布置了抵抗汽车撞击的顶部钢筋，没有设置抵抗来自外部冲击的底部钢筋，当护栏承受风暴潮冲击时，由于阻水面积大，又没有设计抵抗外部冲击力的构造，容易造成护栏与桥面板的连接失效，如图3所示。

从桥梁上部结构的损伤可以看出，梁体在台风风暴潮、浪作用下发生损坏的原因大部分都可以归结于梁体无约束运动。所以，支座的设计对上部结构的安全具有重要意义。对于跨海桥梁，为了适应台风、龙卷风、雷暴等极端灾害性天气，跨海桥梁支座不仅要适应桥梁的变形行为，还要考虑极端灾害的影响。

1.4下部结构破坏

桥梁下部结构是支撑上部结构的建筑物，包括桥墩、桥台和墩台基础。由于跨海桥梁桥墩的墩基础直接立于海中，受到台风风暴潮、浪的影响更加直接。公路桥墩大部分为圆形桥墩、矩形桥墩，铁路桥墩大部分为圆端形桥墩。桥墩的阻水面积较小，对水流还有引导作用，抵抗台风风暴潮、浪的能力较强，但在极端的台风风暴潮和浪作用下也可能发生损伤，其破坏形式如图4所示。

由图4可知，台风风暴潮、浪对桥梁下部结构的破坏极为严重。沿海桥梁下部结构为柱式桥墩时，在台风风暴潮、浪的作用下，柱发生了移位、倒塌、开裂、倾斜及柱与盖梁连接失效等损坏形式。破坏的原因主要有水动力显著、上部结构撞击、下部结构严重冲刷等。同时可以发现，发生严重结构性损伤的位置为盖梁和柱的连接处，这里是抗震设计中布置塑性铰的位置，所以，不能排除台风风暴潮、浪对结构的作用为周期性的作用。从部分柱的倾斜和移位也可以发现，在台风风暴潮、浪的短时间作用下，柱底周围发生了严重的冲刷。

桥梁下部结构施工工期长，对于深海桥梁还具有施工难度大等特点，当桥梁下部结构发生结构性损伤无法使用时，会严重制约线路的抢修。风暴潮引起的桥墩冲刷值得进行深入研究，并需要采用有效的防护措施来保护桥墩，使其免于遭受严重冲刷。同时，需加强桥墩与盖梁之间的连接，使桥墩和盖梁在风暴潮的作用下不至于连接失效。

2风暴潮、浪作用下桥梁下部结构受力数值模拟

下部结构发生损坏时，灾后工作将受到严重制约。采用CFD分析软件Flow-3D，通过选用不同的水位模拟不同的风暴增水情形，研究风暴潮、浪作用下桥梁下部结构的受力特点。假设水位对波浪冲击桥墩的影响即风暴潮对波浪冲击桥墩的影响，不考虑风暴潮和波浪的非线性耦合。在数值计算中，考虑水体为连续不可压缩粘性流体，通过雷诺时均法求解流体连续性方程和动量守恒方程（RANS方程），得到流体的运动状态，同时，采用k-w模型考虑水体的湍流。很多文献都对RANS方程、k-w模型等流体计算方法进行了介绍，本文不再赘述。数值试验设计的桥梁下部结构为半径（D/2）1.8m的近海圆形桥墩，未进行缩尺;波浪选用孤立波，共设计了孤立波波高（H）分别为3、5、7m，水深（d）分别为10、12、15m的9个数值试验工况。

先建立如图5所示的三维数值水槽。在兼顾数值准确性的前提下，采用嵌套网格方式节约计算资源。数值水槽长150m、宽20m、高20m，在人口侧选用源人口造波法进行造波，并在出流侧设置消波块进行消波。数值水槽的边界条件设置为：前侧为波浪边界，后侧为出流边界，两侧为对称边界，上侧为压力边界（大气压），下侧为墙边界。为验证所建立数值模型的准确性，采用波高为5m，水深分别为10、12、15m的数值波形与理论波形进行对比。数值波形取水槽前端3m处水质点的波动轨迹，对比结果如图6所示。由图6可知，数值水槽造波波形与理论波形吻合良好，其中水深10m、波高5m时，造波的最大波高衰减2.2%;水深12m、波高5m时，最大波高衰减1.6%;水深15m、波高5m时，最大波高衰减0.5%，说明造波受破碎指标H/d（波高/水深）的影响。在同一网格下，破碎指标越大，造波的最大波高误差越大，当达到破碎临界时，孤立波将发生破碎。总体而言，数值水槽造波波形与理论波形吻合良好，满足研究需求。在正式开展数值试验时，将模型置于水槽中部，如图5所示。

3不同风暴增水下桥梁下部结构的波浪

各工况数值模拟得到的水平波浪力、纵向弯矩（矩矢垂直波浪前进方向）如图7、图8所示。由图7、图8可知，最大水平波浪力随着水深的增加而减小;除水深10m、波高7m的工况，最大纵向弯矩随着水深的增加而增加，最大水平波浪力与最大纵向弯矩几乎发生在相同时刻。

孤立波的性质由破碎指标H/d决定，不少学者给出的H/d的极限值介于0.714～1.03之间，常用值为0.78。对于水深10m、波高7m的工况，其破碎指标H/d=0.7，临近破碎，所以，水平波浪力和纵向弯矩时程较为陡峭，最大水平波浪力與最大纵向弯矩较破碎指标H/d较小的工况大很多。波浪破碎时与结构物作用，伴随着更剧烈的能量交换，所以，沿岸工程应更关注破碎波与结构物之间的响应;但这不是讨论的重点，本文仅讨论未临近破碎的波浪与结构物的响应，对临近破碎的工况10-7（10-7代表水深为10m、波高为7m的工况）进行了剔除，只关心未临近破碎的孤立波作用圆形桥墩的受力规律。

为得到未破碎（H/d较小）的孤立波作用圆形桥墩的受力规律，采用无量纲分析的方法建立了无因次的最大水平波浪力、最大纵向弯矩力矩与破碎指标（H/d）的关系，如图9、图10所示。无因次的最大水平波浪力采用指数形式拟合，最大纵向弯矩力矩采用线性拟合，拟合结果如图9、图10及式（1）、式（2）所示。采用式（1）、式（2）、式（3）即可进行未临近破碎的孤立波作用圆形桥墩的最大波浪力及最大弯矩计算，这将有助于近海工程结构波浪荷载的初步估计。

4结论与设计施工建议

根据台风风暴潮、浪共同作用下中小型跨海桥梁破坏情况调研，分析了台风风暴潮、浪共同作用下跨海桥梁上部结构和下部结构的破坏原因，并开展了CFD三维数值模拟，探讨了不同风暴增水对波浪作用桥墩的影响，得到以下结论：

1）根据桥梁上部结构破坏形式调研，发现风暴潮、浪作用下跨海桥梁上部结构的破坏形式主要有梁体离位、梁体脱落和梁体断裂，主要原因为风暴潮伴随的海面异常升降、风暴潮的冲击力过大、桥梁结构设计不合理（阻水面积大、上部结构与下部结构连接欠考虑）导致中小跨径跨海桥梁梁体在台风风暴潮、浪共同作用下被迫发生无约束运动。

2）根据桥梁下部结构破坏形式调研，发现下部结构破坏主要表现为结构偏移、倒塌以及盖梁连接失效破坏，主要原因为下部结构严重冲刷、水动力显著、上部结构撞击。

3）桥梁下部结构最大水平波浪力随着水深的增加而减小;最大纵向弯矩在未临近破碎的波浪作用下随着水深的增加而增加。给出了不同风暴增水下近海波浪作用下桥墩的最大波浪力、最大纵向弯矩力矩的经验公式，在确定桥墩直径、桥址水深、来波波高、天文潮+风暴潮总水位后，即可对桥墩受到的波浪力进行简单的初步估计。

在台风风暴潮、浪的共同作用下，桥梁受力是一个复杂的空间问题;浮力、各方向的冲击力都是造成桥梁损坏的因素。为了使跨海桥梁免遭严重破坏，必须制定合理的解决对策，但一味地增大结构抗力不经济。根据台风风暴潮、浪共同作用下的桥梁破坏分析，可采取如下对策：

1）采用阻水面积较小的护栏，采用带有孑L洞的桥面板，释放桥面板下聚积的空气和压力。

2）妥善处理上部结构与下部结构的连接，使支座在竖向、横向都具有耗能能力，并允许具有一定的位移;为了保护下部结构，允许在一定冲击力下发生剪切破坏。

3）采用具有导流作用的桥墩，以减少风暴潮的冲击力。

4）设置防护装置，如在桥墩周围设置防波堤等。

5）加大基础埋深，对桥墩发生的局部冲刷进行防护。