桥梁桁 颤振空气动力不稳定 返回
桥梁空气动力·抗风设计

桥梁桁 颤振空气动力不稳定 — 塔科马纳罗斯桥教训

长大桥梁桁断面发生的空气动力不稳定振动(颤振)和涡致振动的临界风速实时计算工具。通过改变桥梁类型、跨径、桁宽、桁高、固有频率等参数,评估设计风速的裕度和涡放出频率的锁定情况,学习塔科马纳罗斯桥1940年崩塌事件的教训,掌握现代抗风设计的要点。

参数设置
桥梁类型
自动设置断面形状系数sF(流线形越稳定)
中央跨径 L
m
桁宽 B
m
直接影响颤振临界风速的关键参数
桁高 D
m
影响涡放出频率 f_v = St·U/D
单位长度质量 m
kg/m
竖向弯曲固有频率 f_b
Hz
扭转固有频率 f_t
Hz
f_t/f_b的比值越大越抗颤振
设计风速 V
m/s
百年复现期望值的10分钟平均风速
计算结果
质量比 μ
颤振临界风速 U_cr (m/s)
设计风速裕度 U_cr/V
还原风速 U_red
涡放出频率 f_v (Hz)
锁定风速 U_lock (m/s)
桁断面·风流·扭转振动·涡放出

蓝色箭头表示风,黄色S形表示桁的交替涡(Kármán涡街),桁断面因扭转振动而倾斜。当裕度低于1.2时,振幅增大,颤振发散模式启动。

颤振临界风速 vs 频率比 f_t/f_b
桥梁类型别 U_cr 比较(相同条件)
理论·主要公式

$$\mu = \frac{m}{\pi\,\rho\,(B/2)^{2}}, \qquad U_{cr} = 2.5\,B\,f_{t}\,\frac{\sqrt{\mu}}{s_{F}}$$

质量比μ与Selberg公式的颤振临界风速U_cr。m:单位长度质量,ρ:空气密度(1.225 kg/m³),B:桁宽,f_t:扭转固有频率,s_F:断面形状系数(流线形越小)。

$$U_{red} = \frac{V}{f_{t}\,B}, \qquad f_{v} = \frac{St\cdot V}{D}, \qquad U_{lock} = \frac{f_{b}\,D}{St}$$

还原风速U_red、涡放出频率f_v、涡致振动锁定风速U_lock。St=0.12(钝头桥梁桁的典型值),V:设计风速,D:桁高,f_b:竖向弯曲固有频率。当f_v接近f_b时进入锁定状态。

桥梁桁 颤振·空气动力不稳定 — 塔科马纳罗斯桥教训

🙋
塔科马纳罗斯桥就是那个著名的"摇晃崩塌"视频吧?据说当时的风速其实不太强,这是真的吗?
🎓
是的,1940年11月7日发生崩塌时,风速约19 m/s(约38英里/小时),在日本相当于强风注意报水平,远非台风级别。设计者当时只考虑了静风压(横向推力),完全没想到桁会"自己开始扭转"的颤振现象。19 m/s在通常的风压计算中不足以造成破坏,但由于桁断面是薄H形板梁,气动力与扭转振动耦合,就超过了临界风速。教科书的风压计算方法完全无法预测这种破坏模式。
🙋
我在左边把"桥梁类型"改成"桁架桁",马上裕度就降到1以下,判定为危险。但桁架桁看起来比"悬索桥"更坚固啊,为什么反而变弱了?
🎓
很好的问题!"坚固"分两个维度。桁架桁的弯曲强度和刚度确实更高,但从空气动力学角度它很不利。桁架的空隙和尖角会产生大量涡,这些涡的力与桁的扭转振动容易耦合。本工具用断面形状系数s_F表示:桁架=1.3(不利),扁平箱桁=0.7(有利)。公式U_cr = 2.5·B·f_t·√μ/s_F中,s_F在分母,所以s_F越大U_cr越小。这正是为什么明石海峡大桥选择了流线形补刚桁而不是桁架。
🙋
那么要提高U_cr,就是让桁"重一点、宽一点、扭转刚度高一点",对吧?但还有个"还原风速 U_red"也在显示,这是什么意思?
🎓
U_red = V/(f_t·B)是一个无量纲风速,表示"桁完成一次扭转的时间内,风推进了桁宽的多少倍"。这是风洞试验的关键相似性参数——只要U_red相匹配,模型和实桥的流场就是相似的。经验上,扁平断面颤振发生在U_red≈8~12的范围,本工具的默认值(U_red=10)就在警戒区。当你增加风速时,U_red也会增大,会进入涡励振锁定区。桥梁设计者必须计算设计风速下的U_red,并对比风洞试验数据。
🙋
还有"涡励振"和"颤振"两个词都出现了,但它们有什么不同?都是风吹导致的摇晃啊?
🎓
关键差别在"是否自己增长振幅"。涡励振(VIV)是由Kármán涡的放出频率f_v = St·U/D等于桥的竖向弯曲固有频率f_b时产生的强制振动,这只在风速进入特定"锁定区域"时发生,离开该区域后自然衰减。颤振是自励振动,超过临界风速就会指数级发散,无法停止。塔科马桥是扭转颤振,振幅膨胀到±35°左右,桁断裂。VIV是疲劳和乘坐舒适的问题,颤振是整桥致命倒塌的问题。
🙋
最后一个问题,现在的长大桥梁怎样防止颤振?单纯加重是不够的吧?
🎓
现代长大桥采用四层递进策略。(1)断面形状流线化——明石海峡大桥(中央跨1991m)和本州四国连络桥采用扁平箱桁,从源头抑制涡放出。(2)中央切口——在桁中央开竖缝,连通上下表面气压,消除非对称性。(3)整流罩和小翼——装设各种空气动力装置。(4)制振装置——调谐质量阻尼器(TMD)或粘性阻尼器吸收振动能量。从设计流程看,先用像本工具这样的概算估算U_cr,再用部分模型风洞试验实测,最后用巨大全桥模型(10~30米级)确认,三个阶段都要过。设计风速的1.2倍以上余度是业界标准,本州四国连络桥等甚至保持2倍以上的余度。

常见问题

颤振是由风产生的气动力与桥梁自身的竖向弯曲振动、扭转振动耦合,导致振幅自励发散的现象。超过临界风速U_cr时,气流持续向振动系统输入净能量,阻尼无法阻止振幅增长。1940年塔科马纳罗斯桥在风速约19m/s的条件下发生扭转颤振,不到一小时桁完全破坏。即使风速恒定,桥梁也会以"自身想振动"的频率继续振动,与共振不同,只要有优势风向就会导致破坏。
Selberg概算公式用于扁平板桁的颤振临界风速评估,利用桁宽B、扭转固有频率f_t、质量比μ进行计算。质量比μ = m/(π·ρ·(B/2)²)表示桁的单位长质量m与空气惯性的比值,μ越大(重量大、宽度窄)越不易被风扰动,U_cr越高。实际桥梁需要根据断面形状进行大幅补正,桁架桁流动扰乱使U_cr降低,箱桁流线形使U_cr升高。本工具用断面形状系数进行经验补正。
涡致振动(VIV)是桁上下表面交替放出的涡频率f_v = St·U/D与桥梁竖向弯曲固有频率f_b相等时的强制振动。仅在风速进入"锁定区域"时振幅增大,风速离开该区域后自然衰减。颤振是自励振动,超过临界风速后振幅指数级发散并无法停止。VIV的问题是疲劳和舒适性,颤振的问题是整个桁的致命崩塌。
现代长大桥梁设计采用四步骤:(1)通过风洞试验对部分模型和全桥模型进行测试,实测U_cr;(2)采用流线形箱桁(明石海峡大桥、本州四国连络桥)抑制涡放出;(3)装设中央切口、整流罩等空气动力装置;(4)通过调谐质量阻尼器(TMD)或粘性阻尼器增加阻尼。通常将设计风速的1.2~1.5倍以上作为临界风速裕度,本工具的margin也采用这一标准。塔科马桥之后,设计标准从桁架补刚桁转变为扁平箱桁,正是出于空气动力稳定性的考虑。

实际应用

长大悬索桥·斜拉桥基本设计:中央跨径超过500米的悬索桥、斜拉桥,颤振临界风速的确保成为"选择何种桁断面"的第一约束条件。明石海峡大桥(中央跨1991m、扁平箱桁+中央切口)、本州四国连络桥各桥、英国亨伯桥(1410m)等,20世纪后期以后的长大桥基本上都采用本工具涉及的Selberg公式及其延伸版作为出发点,并通过风洞试验确认U_cr的方法进行设计。

既有桥梁抗风补强:1980年代以前建成的悬索桥、桁架桥中,部分桥梁对当代设计风速的裕度不足。典型案例如加拿大Leonard Moissan桥、英国旧式桥等,通过后期装设TMD、增加整流罩、对桁进行中央切口改造等抗风补强工程来提升颤振余度。用本工具把桥梁类型设为"桁架桁"并输入既有条件,就能直观理解为何需要补强。

人行桥·中小桥梁涡励振对策:径间100米左右的人行桥和中等规模桥梁,如果桁高D小、竖向弯曲频率f_b低,锁定风速U_lock会落在日常风速范围(5~15m/s),导致每天都有涡励振晃动的事故。伦敦千禧桥(人行桥、横向同步步行加振)、东京湾横断道路付属桥等多个案例都是设计后追加TMD才解决问题。本工具中锁定风速低于设计风速时的显示,就是设计阶段需要采取对策的信号。

CFD·风洞试验的前期检讨:详细的CFD(计算流体力学)解析和风洞试验成本巨大(每个工况数百万至数千万元)。用本工具的概算对各桁方案进行筛选,判断"临界风速裕度是否低于1.5",可以筛出需要进行本格试验的方案。反过来,CFD结果与这个概算差异过大时,说明模型化(特别是振动模态和气动导数的设定)可能有问题,可以作为理智检查工具。

常见误解和注意事项

最大的陷阱是把颤振与共振混淆。共振是外力的频率等于固有振动频率时振幅增大的强制振动,外力停止后摇晃也停止。颤振是自励振动,风速超过临界U_cr后,即使风的具体频率不定,桥梁也会自选频率导致幅度膨胀。看塔科马桥视频很容易误解成"当日风速恰好等于桥的共振频率",其实应该理解为"风速超过了临界值,所以桥自己选择了一个频率发散"。增加阻尼对VIV有效,对颤振效果不大——颤振的根本对策是改善断面形状。

其次是盲目接受Selberg公式的数值。本工具的公式U_cr = 2.5·B·f_t·√μ/s_F是对扁平平板的古典近似,再加上桥梁类型的经验补正系数s_F。实际设计用的是从风洞试验实测的气动导数,代入Scanlan方程求解U_cr。本工具用于建立定性认识,数值不能直接用于实设计。特别是s_F=1.3(桁架)与s_F=0.7(箱桁)的系数差,实桥往往差2~3倍,远大于本工具的差值。

最后是为了裕度1.2就认为安全。设计风速V是"百年复现期的10分钟平均风速",实际阵风是平均风的1.3~1.6倍,考虑突风还要更大。另外V与地表粗度(城市vs海上)、桥梁高度(地上100m与10m差很大)有关,本工具假设的是"桁高处的代表风速"。必须单独检讨现地的风观测数据和地形效应(峡谷风、吹上气流)。业界标准是裕度1.5以上,本州四国连络桥等甚至2.0以上,这才是真正的安全。

使用指南

  1. 输入桥梁桁的跨径(典型值:100~500m)、桁宽(12~40m)、桁厚(1~3m)
  2. 设定桁的单位长质量(钢铁:8~12吨/m,混凝土:15~25吨/m),并根据纵横比确定空气动力阻尼系数
  3. 按计算按钮,输出颤振临界风速U_cr、还原风速U_red、涡致振动锁定风速U_lock。用U_cr/V与设计基准风速(例:50年复现期40 m/s)的比值评估安全余度

具体计算示例

塔科马纳罗斯桥复现工况:跨径520m、桁宽12m、桁高2.4m、单位质量8.5吨/m(钢制)。对于纵横比H/B=0.2的非流线形桁,颤振临界风速U_cr≈18 m/s,还原风速U_red=U_cr×B/ω_n≈3.5为低值。设计基准风速40 m/s时,U_cr/V=0.45,说明安全余度严重不足。涡致振动锁定风速U_lock≈12 m/s,共振发生风险显著。

实务注意事项