渦励振（VIV）解析
理论与物理
涡激振动现象的物理背景
我听说海洋立管和烟囱在风中摇晃不仅仅是风荷载，具体是什么机制呢？
当流体流过圆柱状结构物时，会在尾流中交替释放出卡门涡街。这种涡旋释放伴随着作用在结构物上的周期性升力波动，导致结构物振动。这就是涡激振动。
涡旋脱落频率 $f_s$ 由斯特劳哈尔数 $St$ 来表征。
其中 $D$ 是圆柱直径，$U$ 是均匀流速。在雷诺数 $300 < Re < 3 \times 10^5$ 范围内，通常使用 $St \approx 0.2$。
当涡旋脱落频率接近结构物的固有频率时，是不是就很危险了？
没错。当 $f_s$ 接近结构物的固有频率 $f_n$ 时，会发生锁定现象。涡旋脱落频率与结构物的振动频率同步，振幅急剧增大。
该区域对应的约化速度 $U^* = U/(f_n D)$ 大致在 $4 < U^* < 8$ 范围内。
控制方程
因为流体和结构相互作用，所以需要联立求解双方的方程吗？
流体侧需要求解纳维-斯托克斯方程。假设不可压缩，则有：
结构侧用运动方程表示。简化的单自由度模型为：
其中 $F_L(t)$ 是流体施加的升力。质量比 $m^* = m / (\rho_f D^2 L)$ 和阻尼比 $\zeta$ 主导响应。
这两个方程通过界面条件连接起来，对吧。
在流固界面上施加速度匹配条件和力平衡条件。
使用ALE方法来处理移动网格是标准的做法。
卡门涡的“规律性”——湍流中的秩序
当流体流过圆柱时，在雷诺数40~200左右的层流区域出现的卡门涡街是自然界“有序的紊乱”的典型例子。涡旋脱落频率由斯特劳哈尔数 St = fD/U ≈ 0.2 惊人地稳定。这意味着无论是1mm还是1m直径的圆柱，只要对流速进行归一化，就会得到相同的无量纲数，体现了“相似律”之美。烟囱、海洋立管、桥梁的吊索，甚至人体的血管，只要有圆柱状结构物，就会出现卡门涡。在工程上，涡激振动虽被视为“麻烦制造者”，但其涡旋的规律性被用于“涡街流量计”（通过涡旋频率测量流量），在化工和石油工厂中广泛应用。
各项的物理含义
- 结构-热耦合项：温度变化引起的热膨胀诱发结构变形，变形又影响温度场。$\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% |
数值解法与实现
ALE法实现的网格移动
结构物移动的话，流体网格也必须跟着变形吧？这个怎么处理呢？
ALE法中引入网格速度 $\mathbf{w}$，修正纳维-斯托克斯方程的对流项。
网格的移动通过拉普拉斯方程或弹簧类比法决定。大变形时需要网格重划分。
重划分的时机怎么决定呢？
监控单元质量指标（长宽比、偏斜度），低于阈值时自动执行重划分。Ansys Fluent中可以通过动态网格功能设置。STAR-CCM+的变形网格也是类似的机制。
耦合算法
按什么顺序求解流体和结构很重要吗？
非常重要。主要分为两大类。
弱耦合（Loose/Weak coupling）: 每个时间步依次求解流体→结构一次。计算速度快，但对于质量比较小（$m^* < 5$ 左右）的系统容易不稳定。
强耦合（Strong coupling）: 在每个时间步内进行子迭代，使界面条件收敛。稳定但计算成本高。
强耦合中，使用Aitken松弛或IQN-ILS（Interface Quasi-Newton Inverse Least Squares）法来加速收敛。
像海洋立管这样质量比小的情况，强耦合是必须的，对吧。
没错。水下结构物的 $m^* \approx 1$ 的情况也存在，弱耦合会因附加质量效应引发数值不稳定。
湍流模型的选择
VIV分析应该用RANS还是LES呢？
要准确捕捉涡旋脱落，LES或DES更理想。RANS的 $k$-$\omega$ SST 模型虽然能捕捉2D涡脱落的时机，但无法表现3D涡结构和展向相关性。
| 湍流模型 | 涡脱落精度 | 计算成本 | 推荐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分析中最重要的概念就是锁定。通常，涡旋脱落频率随流速成比例变化。但当流速接近结构的固有频率时，涡旋的频率会被结构“拖拽”并锁定在固有频率上。这种状态就是锁定，即使流速改变，涡旋也持续以相同频率脱落，振动持续并增大。锁定的速度范围（锁定带宽）取决于结构的质量比和阻尼比。质量轻、阻尼小的结构（如钢制细长烟囱）具有较宽的锁定区域，更加危险。VIV分析的核心在于绘制“在锁定区域内，哪些模态在哪些流速下发生共振”的图谱。
整体式方法
将所有物理场作为一个联立方程组同时求解。对强耦合问题稳定，但实现复杂，内存消耗大。
分区法（分离迭代法）
各物理场独立求解，在界面交换数据。易于实现，可利用现有求解器。适用于弱耦合。
界面数据传递
最近邻法（最简单但精度低）、投影法（守恒性好）、RBF插值（对网格不一致鲁棒性强）。需要在守恒性和精度之间取得平衡。
子迭代
在每个耦合步内进行足够的迭代，确保界面条件的协调性。残差基准基于各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛系数。防止因过度松弛导致发散，是一种加速收敛的自适应方法。
稳定性条件
注意附加质量效应（流固耦合中结构密度≈流体密度时）。不稳定时可应用Robin型界面条件或IQN-ILS法。
Aitken松弛的比喻
Aitken松弛类似于“平衡跷跷板”。一方推得太用力，另一方就会弹起，反弹又导致推得更用力——为了抑制这种振荡，自动调整推力大小的就是Aitken松弛。当耦合迭代振荡不收敛时，根据上一次的修正量自动调整下一次修正量的自适应方法。
实践指南
分析模型的构建步骤
实际开始VIV分析时，应该按什么步骤进行呢？
首先是几何准备。对于圆柱直径 $D$，确保计算域在流入方向至少有 $20D$，尾流方向至少有 $40D$，横向至少有 $20D$。这是《Journal of Fluids and Structures》等期刊推荐的最小尺寸。
网格需要细化到什么程度呢？
圆柱表面的第一层单元高度取决于壁面 $y^+$。URANS需要 $y^+ \approx 1$，LES也需要 $y^+ < 1$。确保圆柱周围至少有200个以上的分割，并使用O型网格解析边界层。
| 网格参数 | URANS推荐 | LES推荐 |
|---|---|---|
| 圆柱周向分割数 | 200以上 | 360以上 |
| 壁面第一层高度 |