涡激振动（VIV）
涡激振动（VIV）的理论基础
概要
老师！ 今天是涡激振动（VIV）的话题吧？ 具体是什么呢？
卡门涡产生的结构振动。立管、烟囱、输电线。
等等等等，由卡门涡产生的结构，也就是说这种情况下也能用吗？
控制方程
我明白了…涡激振动看似简单，其实非常深奥呢。
离散化方法
怎样在计算机上实际求解这些方程呢？
使用有限元法（FEM）进行空间离散化。组装要素刚性矩阵，构建整体刚性方程。
矩阵求解算法
矩阵求解算法，具体是怎么回事呢？
通过直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程。对于大规模问题，带预处理的迭代法最有效。
| 求解法 | 分类 | 内存使用量 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小～中规模 |
| Cholesky分解 | 直接法（对称正定） | O(n²) | 小～中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说，如果在有限元法部分掉以轻心，之后会吃苦头呢。我会记住的！
商用工具中的实现
那么，做涡激振动（VIV）用什么软件呢？
| 工具名 | 开发者/现在 | 主要文件格式 |
|---|---|---|
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| Ansys CFX | Ansys Inc. | .cfx, .def, .res, .ccl |
| OpenFOAM | 开源软件（OpenCFD/ESI、OpenFOAM Foundation） | 字典文件（blockMeshDict等）, .foam |
| COMSOL Multiphysics | COMSOL AB | .mph |
供应商系统和产品整合历程
各个软件的起源，是不是很有戏剧性啊？
Ansys Fluent
接下来是Ansys Fluent的话题呢。内容是什么呢？
由Fluent Inc.开发。2006年被Ansys收购。基于非结构网格的通用CFD求解器。
现在的归属: Ansys Inc.
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题呢。内容是什么呢？
由CD-adapco开发。2016年被西门子收购，整合到Simcenter品牌。多面体网格为特色。
现在的归属: Siemens Digital Industries Software
听到这里，我终于理解了开发历史为什么重要！
Ansys CFX
请讲一下「Ansys CFX」！
由AEA Technology (英国) 开发的CFX。2003年被Ansys收购。耦合型求解器方式是特色。
现在的归属: Ansys Inc.
啊～，关于开发的故事，太有意思了！ 请再讲讲。
文件格式与互操作性
在不同软件之间交换数据时有什么要注意的吗？
| 格式 | 扩展名 | 种类 | 概要 |
|---|---|---|---|
| CGNS | .cgns | CFD数据 | CFD通用记号系统。CFD结果的标准交换格式。 |
| VTK | .vtk/.vtu | 可视化 | 可视化工具包格式。用于ParaView等。 |
在不同求解器之间转换模型时，要素类型的对应关系、材料模型的兼容性、荷载·边界条件的表示差异需要特别注意。特别是高阶要素或特殊要素（内聚力要素、用户定义要素等）往往无法在求解器间直接转换。
我明白了…格式看似简单，其实非常深奥呢。
实务注意事项
教科书上没有的「现场智慧」有什么吗？
网格收敛性的确认、边界条件的妥当性验证、材料参数的灵敏度分析非常重要。
啊～，涡激振动（VIV）真的很深奥呢… 但是感谢老师的讲解，我总算理清思路了！
嗯，做得很好！ 实际动手操作才是最好的学习方式。有不明白的地方随时都可以来问我。
海底管道「唱歌」——石油业界为VIV花费数十亿美元
深海石油管道暴露在海流中会因涡激振动（VIV）振动，导致疲劳破损。北海和墨西哥湾的海底管道敷设费用高达数百亿日元，因此VIV抑制是生死攸关的问题。目前石油巨头们还在投入巨额研究费来提高VIV预测精度，光是螺旋状突起（helical strake）的安装对策费就要数十亿日元。看似不起眼的VIV理论实际上是超高成本的问题。
涡激振动（VIV）的数值计算方法
数值方法的详细内容
具体用什么算法来求解涡激振动（VIV）呢？
离散化的定式化
使用形状函数 $N_i$ 逼近未知量：
用公式表示如下。
基本方程的离散形式
用公式表示如下。
仅看公式我不太理解…表示什么呢？
连续体的控制方程离散化后，得到以下代数方程组：
这里 $[K]$ 是整体刚性矩阵（或等效的系统矩阵），$\{u\}$ 是未知节点变量向量，$\{F\}$ 是外力向量。
啊，我明白了！ 连续体的控制方程离散化就是这样一个机制呢。
要素技术
「要素技术」听过，但可能没完全理解…
| 要素类型 | 次数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低（剪切锁定） | 低 |
| 四面体2阶 | 二阶 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二阶 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二阶 | 6/15 | 中～高 | 中 |
积分方案
积分方案，具体是怎么回事呢？
听到这里，我终于理解了要素类型为什么重要！
收敛性和稳定性
收敛不了的时候，首先要检查什么？
收敛速度：二阶要素的误差按 $O(h^2)$ 阶递减（光滑解的情形）
我明白了…细分网格看似简单，其实非常深奥呢。
求解器设置建议
具体用什么算法来求解涡激振动（VIV）呢？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数基准 |
| 预处理手法 | ILU(0) or AMG | 根据问题规模 |
| 最大迭代次数 | 1000 | 未收敛则需要重新设置 |
| 内存模式 | In-core | 尽可能 |
风上差分（Upwind）
1阶风上：数值扩散较大但稳定。2阶风上：精度提高但有振荡风险。对高雷诺数流必须。
中心差分（Central Differencing）
二阶精度，但Pe数 > 2时会产生数值振荡。适用于低雷诺数的扩散主导流。
TVD格式（MUSCL、QUICK等）
通过限制器函数抑制数值振荡同时保持高精度。对冲击波和陡峭梯度的捕捉有效。
有限体积法 vs 有限要素法
FVM：自然满足守恒律。CFD主流。FEM：对复杂形状·多物理场有利。SPH等无网格法也在发展中。
CFL条件（库朗数）
显式法：CFL ≤ 1是稳定条件。隐式法：CFL > 1也稳定，但影响精度和迭代次数。LES：推荐CFL ≈ 1。物理意义：1个时间步内信息不传播超过1个网格。
残差监测
连续方程·动量·能量的各残差下降3～4个数量级即判定收敛。质量守恒的残差尤为重要。
松弛系数
压力：0.2～0.3，速度：0.5～0.7是常见初值。发散时降低松弛系数。收敛后提高以加速。
非定常计算的内部迭代
每个时间步内迭代至定常解。内部迭代数：5～20次为目安。残差在时间步间变动时应调整时间步。
涡激振动（VIV）的实务应用
实践指南
老师，请讲讲「实践指南」！
解说涡激振动（VIV）的实务分析流程和注意事项。
分析流程
从最初的一步开始教我！ 首先要做什么？
1. 预处理 (前处理)
- CAD数据导入和形状简化
- 材料特性定义
- 网格生成（要素类型·尺寸的决定）
- 边界条件和荷载条件的设置
2. 求解 (求解)
- 求解器设置（求解方法、收敛基准、输出控制）
- 作业投入和计算执行
- 收敛监测
3. 后处理 (后处理)
- 结果可视化（位移、应力、其他物理量）
- 结果验证和妥当性确认
- 报告编制
网格生成的最佳实践
网格的好坏怎么判断？
要素质量指标
请讲讲「要素质量指标」！
| 指标 | 理想值 | 允许范围 | 影响 |
|---|---|---|---|
| 长宽比 | 1.0 | < 5.0 | 精度下降 |
| 雅可比比 | 1.0 | > 0.3 | 要素退化 |
| 翘曲 | 0° | < 15° | 精度下降 |
| 斜度 | 0° | < 45° | 收敛性恶化 |
| 锥度比 | 0 | < 0.5 | 精度下降 |
网格密度的决定
网格密度的决定，具体是怎么回事呢？
边界条件设置指南
我听说边界条件这里错误的话全部都会报废…
啊，我明白了！ 过约束注意这样一个机制呢。
各商用工具的实现步骤
有各种软件吧？ 各自的特点教我！
| 工具名 | 开发者/现在 | 主要文件格式 |
|---|---|---|
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| Ansys CFX | Ansys Inc. | .cfx, .def, .res, .ccl |
| OpenFOAM | 开源软件（OpenCFD/ESI、OpenFOAM Foundation） | 字典文件（blockMeshDict等）, .foam |
| COMSOL Multiphysics | COMSOL AB | .mph |
Ansys Fluent
接下来是Ansys Fluent的话题呢。内容是什么呢？
由Fluent Inc.开发。2006年被Ansys收购。基于非结构网格的通用CFD求解器。
现在的归属: Ansys Inc.
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题呢。内容是什么呢？
由CD-adapco开发。2016年被西门子收购，整合到Simcenter品牌。多面体网格为特色。
现在的归属: Siemens Digital Industries Software
先生的说明好懂！ 工具名的困惑解开了。
常见失败及对策
初学者容易犯的失败有哪些？ 想事先知道！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格质量不良、不适当的边界条件 | 网格改善、约束条件重新审视 |
| 应力异常大 | 应力特异点、网格依存 | 特异点回避、局部网格细分 |
| 位移非现实 | 材料常数错误、单位系混乱 | 输入数据确认 |
| 计算时间过长 | 不必要的细分、求解效率低 | 网格优化、并行计算 |
质量保证检查表
教科书上没有的「现场智慧」有什么吗？
啊～，涡激振动（VIV）真的很深奥呢… 但是感谢老师的讲解，我总算理清思路了！
嗯，做得很好！ 实际动手操作才是最好的学习方式。有不明白的地方随时都可以来问我。
螺旋状突起这一毫不起眼的对策的真实力
在全球海底管道中采用的VIV抑制方法是在管表面安装螺旋状突起（helical strake）。外形像「螺旋螺钉」。仅仅安装这个，涡产生的规律就会被打乱，VIV振幅可削减80～90%。在试验施工中与CFD和实测相比较，逐步优化突起形状的工作是地道的，但如果能救一根深海管道，就完全值得。简单的物理工巧与高额CFD分析并肩支撑现场。
涡激振动（VIV）的软件比较
商用工具比较
有各种软件吧？ 各自的特点教我！
支持涡激振动（VIV）的主要商用CAE工具的功能比较和各产品历史背景详述。
支持工具列表
那么，做涡激振动（VIV）用什么软件呢？
| 工具名 | 开发者/现在 | 主要文件格式 |
|---|---|---|
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| Ansys CFX | Ansys Inc. | .cfx, .def, .res, .ccl |
| OpenFOAM | 开源软件（OpenCFD/ESI、OpenFOAM Foundation） | 字典文件（blockMeshDict等）, .foam |
| COMSOL Multiphysics | COMSOL AB | .mph |
Ansys Fluent
接下来是Ansys Fluent的话题呢。内容是什么呢？
由Fluent Inc.开发。2006年被Ansys收购。基于非结构网格的通用CFD求解器。
现在的归属: Ansys Inc.
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题呢。内容是什么呢？
由CD-adapco开发。2016年被西门子收购，整合到Simcenter品牌。多面体网格为特色。
现在的归属: Siemens Digital Industries Software
听到这里，我终于理解了开发历史为什么重要！
Ansys CFX
请讲一下「Ansys CFX」！
由AEA Technology (英国) 开发的CFX。2003年被Ansys收购。耦合型求解器方式是特色。
现在的归属: Ansys Inc.
OpenFOAM
OpenFOAM，具体是怎么回事呢？
Imperial College London的开源CFD。OpenCFD Ltd（ESI Group下）和The OpenFOAM Foundation并行开发。
现在的归属: 开源软件（OpenCFD/ESI、OpenFOAM Foundation）
啊，我明白了！ 开发就是这样一个机制呢。
功能比较矩阵
预算和时间都有限，性价比最强的是哪个？
| 功能 | Fluent | Star-CCM+ | CFX | OpenFOAM |
|---|---|---|---|---|
| 基本功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ | ○ |
| GPU支持 | △ | △ | △ | ○ |
转换时的风险
转换时的风险，具体是怎么回事呢？
啊，我明白了！ 不同工具之间的互换就是这样一个机制呢。
许可证形式
「许可证形式」听过，但可能没完全理解…
| 工具 | 许可证 | 特点 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 高价但有官方支持 |
| OpenFOAM | GPL | 免费但支持有偿 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的OSS求解器 |
选择指南
最后到底该选哪个，判断基准能教我吗？
涡激振动（VIV）的工具选定中要考虑以下因素：
啊～，涡激振动（VIV）真的很深奥呢… 但是感谢老师的讲解，我总算理清思路了！
嗯，做得很好！ 实际动手操作才是最好的学习方式。有不明白的地方随时都可以来问我。
风力发电塔也在与VIV对战
海上风力发电的支柱（单桩）直径4～6米，水深30米以上，暴露在海流中，因此VIV会直击疲劳寿命。随着欧洲海上风力的急速扩大，VIV分析的订单猛增。有意思的是工具选择上，流体侧用OpenFOAM精细处理粘性效应，结构侧用ABAQUS进行疲劳评价的混合运用在实务中很多。「哪个工具最强」不如「哪个组合适合问题」是实务的正解。
涡激振动（VIV）的先进研究
先进课题与研究动向
涡激振动（VIV）领域今后怎么发展呢？
来看看涡激振动（VIV）最新研究动向和先进手法。
最新数值方法
接下来是最新数值方法的话题呢。内容是什么呢？
仅看公式我不太理解…表示什么呢？
高性能计算 (HPC) 的支持
| 并行化手法 | 概要 | 适用求解器 |
|---|---|---|
| MPI (区域分割) | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 许多求解器 |
| GPU (CUDA/OpenCL) | GPGPU活用。特别在显式法有效 | LS-DYNA, Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
涡激振动（VIV）的故障排除
故障排除
常见错误及对策
老师也为涡激振动（VIV）彻夜调试过吗？（笑）
1. 收敛失败
收敛失败，具体是怎么回事呢？
症状：求解器未在指定迭代次数内收敛而异常终止
可能原因：
- 网格质量不足（过度歪斜的要素）
- 材料参数设置不当
- 初始条件不适当
- 非线性过强（荷载步数不足）
对策：
- 进行网格质量检查（