雨风振动分析
理论与物理
现象概述
雨振现象具体是指什么？
这是斜拉桥缆索在降雨时受风引发大振幅振动的现象。于1980年代在名港西大桥首次被报告。其特征是干燥时不会发生，由缆索表面的水膜（rivulet）引发流体力学不稳定性所致。
控制方程
如何对水膜行为进行建模？
缆索表面的水膜用薄膜流理论描述。水膜厚度 $h(\theta, t)$ 的时间演化方程为：
$R$ 是缆索半径，$\tau_a$ 是来自空气的剪切应力，$\sigma$ 是水的表面张力。
缆索本身作为2自由度（面内方向和升力方向）振动系统处理。
$F_L, F_D$ 是随水膜形状变化而波动的升力和阻力。水膜位置改变升力系数，从而激发不稳定振动。
水膜的位置是关键呢。
上部水膜（upper rivulet）的周向位置会大幅改变缆索的空气动力学特性。当上部水膜位于有效分离角附近时，升力梯度变为负值，从而产生驰振型不稳定性。
水滴“制造”振动——雨振现象的神奇发现
雨振现象（RWV）在学术上被认知的时间意外地近，起点是1988年丹麦的Hikami和Shet等人报告了斜拉桥缆索的大振幅振动。在此之前，“大雨天桥梁摇晃”的现象虽为人知，却缺乏理论解释。关键在于，雨水形成的溪流状水膜创造了缆索截面的“空气动力学不对称性”。附着在缆索上的水膜随风移动时，升力方向随之改变，做正功从而放大振动。干燥时不发生，雨太大时水流淌下而停止——存在这种微妙的“雨量范围”是RWV理论的核心。
各项的物理意义
- 结构-热耦合项：温度变化引起的热膨胀诱发结构变形，变形又影响温度场。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常例子】夏天铁轨伸长导致缝隙变窄——温度上升→热膨胀→产生应力的典型例子。电子基板在焊接后翘曲也是由于不同材料热膨胀率差异所致。发动机的气缸体因高温部和低温部的温差产生热应力，最坏情况下会导致裂纹。
- 流体-结构耦合（FSI）项：流体压力·剪切力使结构变形，结构变形又改变流体区域，是双向的相互作用。【日常例子】强风下吊桥缆索振动（涡激振动）——风力使结构摇晃，摇晃的结构改变风流，进一步放大振动。心脏血流与血管壁的弹性变形、飞机机翼的颤振（气动弹性不稳定性）也是典型的FSI问题。有时仅单向耦合即可，但变形大时双向耦合是必须的。
- 电磁-热耦合项：焦耳发热 $Q = J^2/\sigma$ 引起温度上升，温度变化又改变电阻，形成反馈回路。【日常例子】电暖炉的镍铬丝通电后发热（焦耳热）变红——温度上升电阻改变，电流分布也变化。IH电磁炉的涡流发热、输电线路温度上升导致的垂度增加也是此耦合的例子。
- 数据传递项：通过插值解决不同物理场间网格不匹配的问题。【日常例子】天气预报中结合“气温数据”和“风速数据”计算体感温度时，若各自的观测地点不同就需要插值——CAE的耦合分析中，结构网格和CFD网格通常也不一致，因此界面处的数据传递（插值）精度直接关系到结果的可信度。
假设条件与适用范围
- 弱耦合假设（单向耦合）：当一方物理场影响另一方而反向影响可忽略时有效
- 需要强耦合的情况：FSI中的大变形、电磁-热耦合中温度依赖性强的场合
- 时间尺度的分离：各物理场特征时间差异大时，可通过子循环提高效率
- 界面条件的协调性：需确认耦合界面处的能量·动量守恒在数值上得到满足
- 不适用的场合：三个以上物理场同时强耦合时，有时需要整体式方法
量纲分析与单位制
| 变量 | SI单位 | 注意点·换算备忘 |
|---|---|---|
| 热膨胀系数 $\alpha$ | 1/K | 钢: 约12×10⁻⁶、铝: 约23×10⁻⁶ |
| 耦合界面力 | N/m²（压力）或N（集中力） | 确认流体侧与结构侧的力平衡 |
| 数据传递误差 | 无量纲（%） | 插值精度依赖于网格密度比。5%以下为目标 |
数值解法与实现
数值方法
这个三者耦合（风-雨-缆索）如何求解？
有全CFD方法和半经验方法。
| 方法 | 技术 | 精度 | 计算成本 |
|---|---|---|---|
| 2D CFD + 水膜模型 | RANS/LES + 薄膜方程 | 中~高 | 中 |
| 3D CFD-VOF | 多相流CFD | 非常高 | 非常高 |
| 准定常空气动力模型 | 风洞实验数据 + ODE | 低~中 | 低 |
| 实验空气动力系数 + FEM | 风洞数据 + 结构FEM | 中 | 低 |
用VOF法直接解析水膜吗？
可以，但水膜厚度约0.1~1 mm，而缆索直径约100~200 mm，尺度比达到1:1000。要同时解析水膜的网格和外部流动的网格非常困难，需要AMR（自适应网格细化）。
实用中多采用2D CFD求解圆柱周围流动，水膜位置和厚度用薄膜方程另行求解的半耦合方法。有研究案例将OpenFOAM的pisoFoam与定制的水膜求解器耦合。
如何处理雨水的表面张力——VOF还是SPH
在RWV分析中，要数值再现缆索表面的水膜，气液界面的追踪方法是关键。代表性的是VOF（流体体积）法，用流体的体积分数表示界面。OpenFOAM的interFoam求解器被广泛使用，但缆索表面接触角（疏水性）的处理方式会极大影响结果。最近备受关注的是SPH（光滑粒子流体动力学），它将水粒子作为粒子追踪，因此能更自然地再现“水滴落”、“桥缆索上形成水膜”的行为。但SPH计算成本比VOF高1~2个数量级，因此用于实际的缆索全长分析尚需时日。
整体式法
将所有物理场作为一个联立方程组同时求解。对强耦合问题稳定，但实现复杂，内存消耗大。
分区法（分离迭代法）
各物理场独立求解，在界面交换数据。易于实现，可利用现有求解器。适用于弱耦合。
界面数据传递
最近邻法（最简单但精度低）、投影法（具有守恒性）、RBF插值（对网格不一致鲁棒性强）。守恒性与精度的平衡很重要。
子迭代
在每个耦合步内进行充分迭代，确保界面条件的协调性。残差基准基于各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛系数。防止因过度松弛导致发散，是加速收敛的自适应方法。
稳定性条件
注意附加质量效应（流体-结构耦合中结构密度≈流体密度时）。不稳定时可应用Robin型界面条件或IQN-ILS法。
Aitken松弛的比喻
Aitken松弛类似于“平衡跷跷板”。一方推得太用力，对面就会弹起，反弹又导致推得太用力——为了抑制这种振荡，自动调整推力大小的就是Aitken松弛。当耦合迭代振荡不收敛时，根据上次的修正量自动调整下次修正量的自适应方法。
实践指南
设计实务中的评估方法
在实际桥梁设计中如何评估雨振现象？
在日本遵循《道路桥抗风设计指南》（日本道路协会）的指导方针。发生条件是风速5~20 m/s、降雨时，缆索倾斜角20~60°，风向与缆索的相对角度在特定范围内。
对策方法
有抑制雨振现象的方法吗？
有多种对策已实用化。
| 对策 | 原理 | 实绩 |
|---|---|---|
| 缆索表面凹坑加工 | 改变水膜的稳定位置 | 明石海峡大桥 |
| 螺旋肋条 | 扰乱水膜的形成模式 | 多摩川天空桥 |
| 设置阻尼器 | 耗散振动能量 | 众多斜拉桥 |
| 连接缆索 | 连接缆索间以改变模态 | 鹤见翼桥 |
能用CFD仿真事前评估对策效果吗？
CFD被用于凹坑和肋条的形状优化。用RANS评估表面粗糙度差异引起的空气动力系数变化，决定最优形状参数。不过，包含水膜行为的完整FSI因计算成本问题仍停留在研究层面。
多摩川桥的RWV实测——台风之夜发生了什么
在日本国内，斜拉桥缆索的RWV也作为实际危害发生过。2009年台风18号期间，关东地区某斜拉桥遭遇强降雨和超过20m/s的风，缆索以最大振幅约80cm振动。现场传感器记录显示，振动于深夜2点左右突然开始，在雨势减弱的早晨5点收敛。此实测数据促成了后续设计标准的修订。如今，大型斜拉桥的新设计事实上必须进行RWV验证，桥梁设计咨询中“风洞试验+喷雨”的RWV专用试验已成为常态。
解析流程的比喻
吹过气球吗？那个瞬间，实际上发生了高级的流体-结构耦合。内部气压（流体）推开橡胶壁（结构）→扩大的壁改变内部压力分布→改变的压力进一步使壁变形……在计算步中重复这种“投接球”的就是FSI分析。
初学者易犯的错误
“单向耦合就够了吧？”——这个判断失误在耦合分析中最危险。结构变形微小的话确实单向就够了，但像心脏瓣膜开闭那样变形大幅改变流路的情况，单向就完全不行了。判断基准是“变形量是否超过特征长度的1%”。超过的话双向耦合是必须的。若误用单向耦合，结果会“看似合理实则大错特错”——这是最可怕的模式。
边界条件的思考方式
耦合界面的数据交换如同“国境的出入境管理”。各国（物理场）有独自的法律（控制方程），但若不在国境（界面）准确管理人和物（力、温度、位移）的往来，两国的经济（能量平衡）就会崩溃。网格不一致时的插值如同“翻译”——误译（插值误差）越小，结果越好。
软件比较
工具比较
有哪些可用于雨振现象分析的工具？
| 工具 | 用途 | 特点 |
|---|---|---|
| Ansys Fluent | 2D/3D CFD | VOF多相流+动态网格FSI |
| OpenFOAM | 2D/3D CFD | 易于实现定制水膜模型 |
| Star-CCM+ | 2D/3D CFD | 内置FSI耦合器，易于设置 |
| Ansys Mechanical | 结构FEM | 与Fluent的双向FSI耦合 |
| Abaqus | 结构FEM | 通过协同仿真与CFD耦合 |
| COMSOL Multiphysics | 多物理场 | 内置FSI模块，易于设置 |
哪个工具最适合初学者？
COMSOL或Star-CCM+的向导式界面容易上手。但若想深入定制水膜模型，OpenFOAM的灵活性更高。研究用途常选OpenFOAM，工业设计则多用Ansys或Star-CCM+。
有开源的FSI耦合库吗？
preCICE库（precice.org）是用于分区耦合的开源库，支持多种求解器（OpenFOAM、CalculiX、SU2等）间的数据交换。它处理界面插值和耦合算法，让用户能专注于各物理场的求解器。
从FORTRAN到Python——耦合分析工具的进化
早期的耦合分析是“手工作业”。用FORTRAN写结构求解器，用另一个FORTRAN写CFD求解器，通过文件交换数据，用shell脚本控制流程——这是2000年代常见的景象。随着多物理场需求的增加，出现了像MpCCI（现为ANSYS的一部分）这样的专用耦合器。2010年代后，preCICE等开源库降低了门槛。如今，在Jupyter Notebook中调用OpenFOAM和CalculiX进行FSI参数研究也成为可能。工具的进化让耦合分析从“专家专属”逐渐走向“普通工程师也能触及”。
工具选择的要点
| 判断基准 | 推奨ツール | 理由 |
|---|---|---|
| 快速验证概念 | COMSOL | GUI设置简单，内置物理模型丰富 |
| 大规模工业分析 | Ansys/Star-CCM+ | 并行效率高，技术支持完善 |
| 定制模型开发 | OpenFOAM | 源代码开放，可自由修改 |
| 多软件协同仿真 | preCICE | 中立库，支持多种求解器组合 |
| 教育/研究 | OpenFOAM + preCICE | 零许可成本，透明度高 |
许可成本比较
| ツール | ライセンス形態 | 年間コスト目安 |
|---|---|---|
| Ansys套件 | 商业许可 | 数万~数十万美元 |
| Star-CCM+ | 商业许可 | 类似Ansys |
| COMSOL | 商业许可 | 数千~数万美元 |
| OpenFOAM | 开源（LGPL） | 免费 |
| preCICE | 开源（LGPL） | 免费 |
注：商业工具的成本因模块、核心数、维护服务而异。开源工具虽无许可费，但需自担开发/支持成本。
学习曲线
| ツール | 初学者到达基本操作 | 到达定制开发 |
|---|---|---|
| COMSOL | 1~2周 | 1~2个月（通过GUI/API） |
| Ansys Workbench | 2~4周 | 3~6个月（需学习APDL/UDF） |
| Star-CCM+ | 2~4周 | 2~4个月（通过宏/Java） |
| OpenFOAM | 1~2个月 | 3~6个月（需C++技能） |
未来展望
未来的耦合分析将走向“智能化”和“云化”。AI用于自动调整耦合参数（松弛系数、时间步长），云平台提供按需的HPC资源。更重要的是，随着数字孪生的普及，实时耦合仿真将成为可能——例如，根据实际桥梁的风雨传感器数据，实时预测RWV风险并调整阻尼器设置。耦合分析正从“设计工具”演变为“运营支持系统”。
なった
詳しく
報告