涡激振动（VIV）分析
涡激振动（VIV）理论基础
涡激振动现象的物理背景
海洋立管和烟囱在风中摇摆不仅仅是风荷载的作用吧，听说有其他机制？
当流流过圆形结构物时，下游会交替产生卡门涡列。随着涡脱落，结构物受到周期性升力变化的作用而振动。这就是涡激振动（Vortex-Induced Vibration, VIV）。
涡脱落频率 $f_s$ 由Strouhal数 $St$ 表征。
其中 $D$ 是圆形直径，$U$ 是均匀流速。在 $300 < Re < 3 \times 10^5$ 范围内通常取 $St \approx 0.2$。
涡脱落频率接近结构物固有振动频率会很危险吧？
完全正确。当 $f_s$ 接近结构物固有振动频率 $f_n$ 时会发生锁定（lock-in）现象。涡脱落频率与结构物振动频率同步，振幅会急剧增大。
这个区域对应的换算流速 $U^* = U/(f_n D)$ 约为 $4 < U^* < 8$。
支配方程
流体和结构会相互作用，需要同时求解两套方程吗？
流体侧求解Navier-Stokes方程。在不可压缩流假设下：
结构侧由运动方程表示。简化的单自由度模型为：
其中 $F_L(t)$ 是流体产生的升力。质量比 $m^* = m / (\rho_f D^2 L)$ 和阻尼比 $\zeta$ 控制着响应。
这两个方程通过界面条件连接是吧。
流固界面处施加速度相容条件和力平衡条件。
通常采用ALE（任意拉格朗日-欧拉）法处理移动网格，这是标准方法。
卡门涡的"规律性"——湍流中的秩序
当流体越过圆形时，在Re数40〜200的层流区，出现的卡门涡列是自然界"有秩序的混乱"的典型代表。涡脱落频率Strouhal数 St = fD/U ≈ 0.2 的稳定性令人惊叹。无论圆形直径是1毫米还是1米，归一化流速后得到相同的无量纲数——这就是"相似律"的美妙之处。烟囱、海洋立管、桥梁吊索、甚至人体血管，只要有圆形结构就会出现卡门涡。从工程角度看，涡激振动是"麻烦的家伙"，但利用涡的规律性开发的"涡流量计"（通过涡频率测量流量）在化工和石油工业中广泛应用。
涡激振动（VIV）数值计算方法
ALE法网格移动
结构物移动时流体网格也要变形吧，怎样处理？
ALE法引入网格速度 $\mathbf{w}$，修正Navier-Stokes方程的对流项。
网格移动由拉普拉斯方程或弹簧类比决定。大变形时需要进行网格重新划分。
网格重新划分的时机怎样判断？
通过监控单元质量指标（长宽比或倾斜度），当低于阈值时自动执行重新划分。Ansys Fluent的Dynamic Mesh功能可以设置这个参数。STAR-CCM+的变形网格也有类似机制。
耦合算法
流体和结构的求解顺序很重要吗？
非常重要。大致分为两类。
弱耦合（Loose/Weak coupling）：每个时间步中流体→结构各求解一次。计算速度快，但质量比小（$m^* < 5$ 左右）的系统容易不稳定。
强耦合（Strong coupling）：每个时间步内进行子迭代，直到界面条件收敛。稳定性好但计算成本高。
强耦合中可用Aitken松弛或IQN-ILS（界面拟牛顿逆最小二乘）法加速收敛。
像海洋立管这样质量比很小的情况，强耦合是必须的吗？
是的。水下结构物的 $m^* \approx 1$ 甚至更小，弱耦合会因为付加质量效应产生数值不稳定。
湍流模型选择
VIV分析中应该用RANS还是LES？
要精确捕捉涡脱落过程，LES或DES更佳。RANS的 $k$-$\omega$ SST模型能捕捉二维涡脱落的时序，但无法表现三维涡结构和展向相关性。
| 湍流模型 | 涡脱落精度 | 计算成本 | 推荐Re范围 |
|---|---|---|---|
| URANS (k-omega SST) | 中等 | 低 | Re < 10^4 |
| DES/DDES | 高 | 中等 | 10^4 < Re < 10^6 |
| LES (Wall-Resolved) | 非常高 | 非常高 | Re < 10^5 |
| LES (Wall-Modeled) | 高 | 高 | Re < 10^6 |
"锁定"现象——结构支配涡的时刻
VIV分析中最重要的概念是锁定（Lock-in）。通常涡脱落频率随流速成正比变化。但当流速接近结构固有振动频率时，涡的频率被结构"拖拽"，最终紧紧贴住固有频率。这就是锁定状态，流速变化时涡仍以同样频率脱落，振动持续甚至增强。锁定速度范围（锁定宽度）取决于结构的质量比和阻尼比。轻质低阻尼的结构（如钢制细长烟囱）锁定区域很宽，更危险。VIV分析的核心是绘制"哪个流速下哪次模态产生共振"的映射图，这对设计至关重要。
涡激振动（VIV）实际应用
分析模型构建步骤
实际开始VIV分析时应该如何进行？
首先准备几何模型。相对于圆形直径 $D$，计算域需要：流入方向至少 $20D$，后流方向至少 $40D$，横向至少 $20D$。这是《Fluids and Structures学报》推荐的最小尺寸。
网格应该多细？
圆形表面第一层网格的壁面 $y^+$ 值很关键。URANS取 $y^+ \approx 1$，LES需要 $y^+ < 1$。圆形周向至少分成200份以上，采用O型网格布置以解析边界层。
| 网格参数 | URANS推荐 | LES推荐 |
|---|---|---|
| 圆形周向分割数 | 200以上 | 360以上 |
| 壁面第1层高度 | $y^+ \approx 1$ | $y^+ < 1$ |
| 边界层内层数 | 20以上 | 30以上 |
| 后流区分辨率 | $D/20$ | $D/40$ |
| 展向分割（3D） | - | $\pi D / 20$ 以下 |
时间步长与Courant数
非定常分析的时间步长也很重要吧。
涡脱落一个周期至少需要200个时间步。具体地说：
最大Courant数URANS取 $CFL < 5$，LES目标是 $CFL < 1$。
要运行多少个周期才能得到稳定的结果？
要排除初始瞬态至少20个周期，再运行20个周期以上来获取统计量。通过监控升力系数的时刻历，确认振幅已稳定。
网格收敛验证
如何验证网格的合理性？
用3个及以上的网格密度等级计算Strouhal数、平均阻力系数 $\bar{C}_D$ 和升力系数振幅 $C_{L,rms}$。用Richardson外推法计算网格收敛指数（GCI）是定量的方法。
其中 $F_s = 1.25$（安全系数），$r$ 是网格细化率，$p$ 是收敛阶数。
烟囱的"螺纹"——螺旋凸起对VIV的压制
大型工厂烟囱和海洋桩基表面常见螺旋状凸起（称为Strake），这是VIV对策。有了Strake，后流涡会沿轴向紊乱，脱落频率变得不规则，结构无法共振。效果显著，可将VIV振幅降低80〜90%。但代价是Strake会增加流体阻力20〜25%，海洋结构需要重新计算潮流力。实际海洋平台的立管会采用"Strake覆盖区"和"裸露区"的混合配置，在设计流速、VIV敏感性和成本间寻找最优平衡。这样简单的突起里面也蕴含着深刻的流体物理。
涡激振动（VIV）软件对比
主要工具VIV支持情况
VIV分析可用哪些软件？
主要商用和开源工具如下。
| 工具 | 流体求解器 | 结构求解器 | 耦合方式 | 特点 |
|---|---|---|---|---|
| Ansys Fluent + Mechanical | Fluent（FVM） | Mechanical（FEM） | System Coupling | GUI耦合设置。支持Aitken松弛 |
| STAR-CCM+ | STAR-CCM+（FVM） | 内置FEA / Abaqus耦合 | 强耦合/弱耦合 | 变形网格优秀 |
| COMSOL Multiphysics | 内置CFD | 内置结构 | 完全单体式 | 单一环境完成。中等规模应用 |
| OpenFOAM + preCICE | OpenFOAM（FVM） | CalculiX/deal.II | preCICE库 | 开源。支持IQN-ILS法 |
| Abaqus + Fluent | Fluent | Abaqus/Standard | MpCCI / System Coupling | 非线性结构求解强大 |
Ansys的System Coupling看起来使用方便。
在Ansys Workbench中设置Fluent/Mechanical间的数据转移面，指定耦合控制参数（松弛系数、子迭代次数、收敛判定值）即可。双向FSI通常设最小/最大迭代次数为5〜15次。
专用VIV预测工具
不用CFD也能评估VIV吗？
在海洋工程领域经常用经验方法。
| 工具 | 方法 | 开发机构 | 用途 |
|---|---|---|---|
| SHEAR7 | 模态叠加法 | MIT | 立管VIV疲劳评估 |
| VIVANA | 频率响应法 | MARINTEK/SINTEF | 立管、管线VIV |
| VIV-RISERPRO | 时间域法 | 2H Offshore | 柔性立管 |
这些工具遵循DNV-RP-C205等规范，不是CFD-FSI的替代品，而是在初期设计阶段用于快速筛选。
CFD-FSI和经验方法可以分别在不同设计阶段使用。
概念设计用SHEAR7或VIVANA，详细设计和问题排查用CFD-FSI，这是实务中的标准做法。
VIV分析工具的"精度vs速度"权衡
VIV分析工具可分为3个精度等级。最简的是"半经验式模型"——代入实验系数的Vandiver式或Sklansky式，计算时间以秒计，但精度粗糙。中等级是"用力系数模型的FEM耦合"——Orcina OrcaFlex、DNV Flexcom等属此，是石油气工业的标准工具。最高精度是"LES-CFD+FEM完全耦合"，每个案例需要数百CPU小时，但能精准重现气泡卷入和三维涡结构。实务中采用"先用半经验式全面筛选→用LES精查重点区域"的两步法最有效。各工具的VIV预测精度用公开基准（VIVACE试验数据等）比较，选型时应确认自家分析条件的相关案例实绩。
涡激振动（VIV）先进研究
多圆柱VIV干涉
实际的立管群和烟囱群是多个圆柱近距离排列，干涉怎样处理？
两个间距为 $L/D$ 的圆柱，串联排列时 $L/D < 3.5$ 左右，后方圆柱的涡脱落被上游圆柱的尾流强烈支配。并排排列时 $L/D < 2.5$ 会出现flip-flopping现象。
用CFD-FSI求解需要多体FSI，分别追踪各圆柱运动。OpenFOAM的overset mesh（Chimera法）或STAR-CCM+的overlapping grid很有效。
尾流诱发振动（WIV）
WIV是什么，与VIV不同吗？
置于上游圆柱尾流中的下游圆柱，由不同机制振动。这是尾流诱发振动（Wake-Induced Vibration, WIV）。WIV类似于抖振不稳定性，振幅有时比VIV还大。
与机器学习的融合
最近有用AI做VIV研究的吗？
用物理信息神经网络（PINN）预测VIV响应的研究在进行。Raissi等（2019, Science）的方法被扩展，用少数传感器数据重构整个流场。
另外用LSTM或Transformer做长尺度立管VIV疲劳预测的时间序列模型也被提出。以SHEAR7结果为学习数据，能把计算成本降低几个数量级。
涡抑制装置优化设计
怎样抑制VIV？
主要的涡抑制装置如下。
| 装置 | 抑制效果 | 阻力增加 | 应用例 |
|---|---|---|---|
| 螺旋Strake | 破坏展向涡相关性 | 30〜50% | 海洋立管、烟囱 |
| 导流罩 | 整流后流 | 阻力减小 | 海底管线 |
| 分流板 | 抑制涡相互作用 | 10〜20% | 烟囱、桥墩 |
| 通孔套筒 | 阻碍涡脱落同步 | 20〜40% | 海洋结构物 |
能用形状优化求Strake的最优间距和高度吗？
用CFD-FSI结合伴随法或贝叶斯优化可以优化Strake间距（通常 $5D$〜$17D$）和高度（通常 $0.1D$〜$0.25D$）。计算成本大，所以多用代理模型辅助。
深海立管的VIV——2000米海底的战斗
2000年代之后，石油开发向"超深水"（水深2000米以上）进军时，立管管VIV（连接海底和钻井船的细管）成为严峻问题。2000米深的立管长2km、直径53cm，暴露在深海流中。麻烦在"高次模态"——浅水是1〜2次模态VIV，深水则出现100多次模态同时励振的"多模态VIV"。ExxonMobil、Shell等巨头用OrcaFlex加定制VIV模型求解，但精度始终不尽人意。如今MIT海洋资助研究所开发的VIVANA（VIV分析软件）已成为深海立管设计的实际标准。
涡激振动（VIV）故障排查
耦合计算发散的情况
VIV耦合计算发散了。什么原因？
VIV-FSI发散的主要原因如下。
| 现象 | 原因 | 对策 |
|---|---|---|
| 第1步发散 | 界面映射不匹配 | 检查流固界面网格。验证节点对应关系 |
| 几步后发散 | 付加质量不稳定性 | 切换强耦合。松弛系数降至0.1〜0.3 |
| 振幅增长后发散 | 网格变形过大 | 增加重新划分频率。考虑overset网格 |
| 特定流速发散 | 锁定区共振 | 细化时间步。增加子迭代次数 |
什么是付加质量不稳定性？
流体付加质量与结构质量相当或更大（$m^* \lesssim 1$）时，弱耦合在界面积累人工能量而发散。Causin等（2005）从理论证明了这一点。对策是强耦合+松弛，或采用Robin-Neumann条件。
Strouhal数与文献值不符
计算得到的Strouhal数偏离0.2很大。
检查这些要点。
- 堵塞比：计算域宽度中圆形直径比超过5%会受壁面影响，St会变化。确保 $D/H < 0.03$
- 时间分辨率：FFT采样不足会降低频率同定精度。用Welch法加窗函数，用足够长的时间数据
- 边界条件：出口太近会约束后流。应距离 $30D$ 以上
- 2D vs 3D：Re > 200 时三维效应显著。2D计算易高估St
与疲劳评估的衔接
VIV分析结果如何用于疲劳寿命评估？
将时刻历应力用雨流法计数，结合S-N曲线和Miner线性累积损伤律计算疲劳损伤度。
使用DNV-RP-C203的S-N曲线，满足 $D_{fatigue} < 1/DFF$（DFF：设计疲劳系数，通常3〜10）。应力集中系数（SCF）的设定对结果影响很大。
从CFD-FSI到疲劳评估有一体化流程吗？
"设计流速外出现振动"——VIV的常见隐患
VIV故障调查中常见"应该不出现锁定的设计流速却振动了"的情况。深入调查才发现，流体在结构周围加速，"局部流速"比设计的"平均流速"高40〜50%。贴近地面或壁面会有伯努利效应加速，多个结构群集（并排烟囱）时上游结构会大幅改变下游流场。另一个看漏的是"模态形状变化"——实际支持条件与设计假设（铰接端支持等）不同，导致固有振动频率严重偏离。处理VIV故障时"确认真实局部流速"和"确认真实固有频率"是第一要务。
细节
错误