雨风振动分析
雨风振动的理论基础
现象的概述
什么是雨风振动现象？
斜拉桥的缆绳在降雨时受风的作用会产生大幅度振动。这种现象首次在20世纪80年代的名港西大桥被报告。其特点是在干燥条件下不会发生，缆绳表面的水膜（rivulet）会引发流体动力学不稳定性。
控制方程
怎样对水膜的行为进行建模？
缆绳表面的水膜由薄膜流理论描述。水膜厚度 $h(\theta, t)$ 的时间演化方程为
其中 $R$ 是缆绳半径，$\tau_a$ 是空气的剪应力，$\sigma$ 是水的表面张力。
缆绳本身被视为具有两个自由度（横向和竖向）的振动系统。
$F_L, F_D$ 是随水膜形状变化而变动的升力和阻力。水膜位置的变化会改变升力系数，从而激发不稳定振动。
看来水膜位置是关键因素。
上部水膜（upper rivulet）的周向位置会大幅改变缆绳的空气动力学特性。当上部水膜位于有效分离角附近时，升力梯度变为负值，从而触发跨越式摆振类型的不稳定性。
水滴"制造"振动——雨风振动的神奇发现
雨风振动（RWV）被学术界认识的时间意外地晚。1988年丹麦科学家Hikami和Shet等人首次报告了斜拉桥缆绳的大幅度振动，这才成为起点。在此之前，"大雨天桥梁摇晃"这一现象虽然存在但缺乏理论解释。关键在于：降雨形成的小卷流状水膜为缆绳截面创造了"空气动力学非对称性"。附着在缆绳上的水膜随风的摇摆而改变升力方向，正功输出使振动不断放大。干燥时不会发生，雨量过大时水流走失而消退——正是这个绝妙的"雨量范围"，才是RWV理论的核心所在。
雨风振动的数值计算手法
数值手法
这三者的耦合（风-雨-缆绳）怎样求解？
有全尺度CFD方法和半实验方法两种。
| 方法 | 手法 | 精度 | 计算成本 |
|---|---|---|---|
| 二维 CFD + 水膜模型 | RANS/LES + 薄膜方程 | 中～高 | 中 |
| 三维 CFD-VOF | 多相流CFD | 非常高 | 非常高 |
| 准定常空气动力模型 | 风洞实验数据 + ODE | 低～中 | 低 |
| 实验空气动力系数 + FEM | 风洞数据 + 结构FEM | 中 | 低 |
用VOF法直接求解水膜？
理论上可以，但水膜厚度约0.1～1 mm，缆绳直径100～200 mm，尺度比达1:1000。同时对水膜和外部流场的网格求解非常困难，需要AMR（自适应网格细化）技术。
实用上通常采用半耦合方法：用二维CFD求解圆柱周围的流场，水膜位置和厚度由薄膜方程另外求解。有研究用OpenFOAM的pisoFoam与自定义水膜求解器耦合的例子。
如何处理雨水的表面张力——VOF还是SPH
在RWV仿真中，要再现缆绳表面的水小卷流，气液界面追踪方法是关键。最常用的是VOF（体积流体）法，用流体的体积分数表示界面。OpenFOAM的interFoam求解器被广泛应用，但缆绳表面的接触角（撥水性）处理对结果影响很大。最近受关注的是SPH（光滑粒子流体动力学）法，将水粒子作为离散粒子追踪，能更自然地再现"水下滑""水膜张贴在缆绳上"等行为。但SPH的计算成本比VOF高1～2个数量级，现阶段要用于整个缆绳长度的实际分析还需时日。
雨风振动的实际应用
设计实务中的评估手法
实际桥梁设计中怎样评估雨风振动？
日本遵循道路桥耐风设计便览（日本道路协会）的指导原则。发生条件是风速5～20 m/s、降雨时、缆绳倾斜角20～60°、风向与缆绳的相对角度在特定范围内。
对策手法
有办法抑制雨风振动吗？
多种对策已投入实用。
| 对策 | 原理 | 实绩 |
|---|---|---|
| 缆绳表面凹陷处理 | 改变水膜的稳定位置 | 明石海峡大桥 |
| 螺旋形挡板 | 扰乱水膜的形成模式 | 多摩川天空桥 |
| 安装减振器 | 散逸振动能量 | 多座斜拉桥 |
| 连接式缆绳 | 改变模态 | 鹤见翼桥 |
能用CFD仿真提前评估对策效果吗？
凹陷和挡板的形状优化中CFD得到应用。通过RANS评估不同表面粗糙度带来的空气动力系数变化，确定最优的形状参数。但包括水膜行为的完整FSI仍因计算成本问题而停留在研究阶段。
多摩川桥实测RWV——台风之夜发生了什么
国内斜拉桥缆绳的RWV也确实造成了实际损害。2009年台风18号期间，关东某斜拉桥遭遇强降雨和秒速超20米的风力，缆绳在最大80厘米振幅的振动中摇晃。现场传感器的记录显示，振动在深夜2点突然开始，清晨5点雨势减弱后消退。这一实测数据之后推动了设计标准的修订。如今大型斜拉桥的新设计中RWV验证事实上已成必选项，桥梁设计咨询公司已把"风洞试验＋喷雨"的专用RWV试验当作常规做法。
雨风振动的软件对比
工具对比
雨风振动分析有什么可用的工具？
| 工具 | 用途 | 特点 |
|---|---|---|
| ANSYS Fluent | 二维/三维 CFD | VOF多相流+动态网格FSI |
| OpenFOAM | 二维/三维 CFD | 易于实现自定义水膜模型 |
| STAR-CCM+ | 二维/三维 CFD | VOF+重叠网格 |
| ABAQUS + Fluent | FSI耦合 | 缆绳结构非线性强 |
| MATLAB/Simulink | 准定常模型 | 参数化研究适用 |
实务中多用准定常模型？
设计阶段常采用风洞实验获取的空气动力系数数据库与MATLAB/Python时刻历程应答分析的组合。全尺度CFD-FSI主要在研究机构或特殊项目中使用。
与风洞实验怎样对比？
标准验证方法是将二维CFD结果与风洞试验（旋转圆柱携带人工水膜）的空气动力系数对比。横滨国立大学和丹麦技术大学（DTU）的风洞数据被广泛引用。
桥梁设计咨询中的RWV工具选择实况
RWV分析涉及气液固三相耦合，这是难题，汎用工具无法直接解决。国内桥梁设计咨询常用的组合是"ANSYS Fluent进行气液界面→ANSYS Mechanical进行结构"的单向耦合或通过FLOWizard、Computational Fluidix等专用流程进行验证。欧洲桥梁咨询公司（COWI、Arup等RWV经验丰富的机构）则通过OpenFOAM加自定义脚本积累多年经验，对他们来说实绩数据比工具本身更有价值。RWV分析的国际标准基准至今未统一，"哪个工具最准确"的问题不如"是否有与实测相符的经历"更受重视。
雨风振动的前沿研究
三维效应和跨向相关性
实际缆绳很长。三维效应重要吗？
缆绳长超过100米，跨向风速和水膜状态不均匀。三维LES研究缆绳跨向空气动力相关性。相关长度越短，缆绳整体振动响应越小。
干流入角跨越式摆振
无降雨也会产生振动？
当缆绳对风倾斜或有偏角时，干燥状态下也会发生跨越式摆振型不稳定振动（干流入角跨越式摆振）。Den Hartog稳定性准则
用于判断这类不稳定性。三维CFD评估空气动力系数至关重要。
基于机器学习的振动预测
能实时预测现场的振动吗？
将气象数据（风速、风向、降雨量）和加速度传感器数据输入LSTM或XGBoost进行学习，研究员已开发出振动发生的事前预测系统。将 CFD-FSI结果纳入学习数据，可改善模型对未知气象条件的泛化性能。
用机器学习预测水膜位置——RWV研究前沿
RWV分析的困难在于，缆绳表面水小卷流的位置随风速、降雨量和缆绳倾斜角不断变化。用传统CFD完整求解所有动态边界条件计算成本巨大。东京大学和横滨国立大学的联合研究团队（2022年）开发出通过LES数据学习的PINNs（物理信息神经网络）来预测水膜位置的方法，计算时间削减至1/50以下。RWV涉及"流体、固体、气液界面"三要素，AI与耦合分析的融合最受期待。
雨风振动的故障排除
水膜不稳定
用VOF法求解水膜，总是消散。
水膜极薄（0.1～1 mm），数值扩散容易致其消失。对策如下
- 将网格在水膜区域细化至水膜厚度的1/3以下（推荐AMR）
- 表面张力模型用CSS法而非CSF法
- 时间步CFL < 0.3
- 选用带anti-diffusion项的VOF格式（HRIC、CICSAM）
二维模型无法再现振动怎么办？
列举检查项。
| 检查项 | 详细 |
|---|---|
| 风速范围 | 确认在5～20 m/s范围 |
| 结构阻尼比 | Scruton数 $Sc = 2m\delta / (\rho D^2)$ 越小越易产生振动 |
| 水膜初始位置 | 确认初值中是否给定上部rivulet |
| 空气动力系数更新频率 | 准定常模型更新过粗会遗漏不稳定性 |
| 迎角效应 | 确认正确反映了缆绳倾斜角 |
Scruton数多大时危险？
$Sc < 10$ 时雨风振动风险高。设计基准要求 $Sc > 10$（欧洲）或安装适当的减振器。
"加了减振胶，振动反而增大"——现场反作用案例
RWV对策中常见的是在缆绳上加装制振装置，但设计不当会适得其反。某桥梁在缆绳下端后装粘弹性减振器后，振动反而比安装前更大。调查发现，减振器的固有频率与缆绳-减振器耦合系统频率相同，触发了共振。RWV的发生机制本身对气象条件高度敏感，"某气象条件能抑制但另一条件会恶化"的情况常见。现在已成为标准做法：安装前必须用耦合仿真验证多个气象情景。
更详细
报告