管束流力振动解析
管束流力振动的理论基础
概要
老师！今天是管束流力振动解析的课题呢？ 是什么样的内容呢？
热交换器管束的流力弹性不安定现象。Connors的不安定速度基准。螺距比和减衰的影响。
明白了。那么如果热交换器管束的流力弹性已经满足条件的话，基本上就没问题了吧？
支配方程
离散化手法
这个方程，在计算机上实际怎样求解呢？
使用有限元法（FEM）进行空间离散化。组装单元刚性矩阵，构建全局刚性方程。
矩阵求解算法
矩阵求解算法，具体是什么意思呢？
通过直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程组。对于大规模问题，预处理迭代法很有效。
| 求解法 | 分类 | 内存使用量 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小～中规模 |
| Cholesky分解 | 直接法（对称正定） | O(n²) | 小～中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说在有限元法的地方马虎的话，后面就会吃亏啊。铭记于心！
商用工具中的实现
那么，管束流力振动解析可以用什么样的软件呢？
| 工具名 | 开发方/现所属 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Mechanical（旧ANSYS Structural） | ANSYS Inc. | .cdb, .rst, .db, .ans, .mac |
| Ansys Fluent | ANSYS Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
供应商系谱与产品整合历史
各个软件的来历，好像都有些戏剧性的故事呢？
COMSOL Multiphysics
请给我讲讲"COMSOL Multiphysics"！
1986年在瑞典创立。最初以MATLAB联动的FEMLAB开始，后来改名为COMSOL。多物理场领域有强项。
现所属: COMSOL AB
Ansys Mechanical（旧ANSYS Structural）
请给我讲讲"Ansys Mechanical"！
1970年由Swanson Analysis Systems Inc.（SASI）开发。基于APDL（Ansys Parametric Design Language）。
现所属: ANSYS Inc.
Ansys Fluent
接下来是Ansys Fluent的话题。什么内容呢？
由Fluent Inc.开发。2006年被ANSYS收购。基于非结构网格的通用CFD求解器。
现所属: ANSYS Inc.
哦，这样啊！在瑞典创立是这个意思啊。
文件格式与互操作性
在不同软件之间传递数据时有什么注意事项吗？
| 格式 | 扩展名 | 类型 | 概要 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | 符合ISO 10303的3D CAD数据交换格式。支持形状+PMI。 |
| IGES | .igs/.iges | 中立CAD | 最初的CAD数据交换规范。曲面数据互操作性有问题。正向STEP迁移。 |
| CGNS | .cgns | CFD数据 | CFD通用符号系统。CFD结果的标准交换格式。 |
| VTK | .vtk/.vtu | 可视化 | 可视化工具包格式。可用ParaView等工具。 |
在不同求解器之间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载与边界条件的表示差异。特别是高阶单元或特殊单元（内聚单元、用户定义单元等）通常无法在求解器间直接转换。
明白了…格式看起来简单，实际上很深奥呢。
实务上的注意事项
教科书上没有的"现场智慧"之类的东西有吗？
网格收敛性的确认、边界条件的妥当性验证、材料参数的敏感性分析非常重要。
管束流力振动解析的整体面貌掌握了！ 明天开始在实务中多加留意。
好的，干得不错！ 实际动手操作是最好的学习方法。有不懂的地方随时来问。
核反应堆蒸汽发生器的管束振动——切尔诺贝利之后世界学到的教训
管束流力振动对最造成严重影响的案例之一是核电站蒸汽发生器。压水堆（PWR）蒸汽发生器内有数千根细传热管组成的管束，在高压高温的一级冷却水与二级侧蒸汽之间进行热交换。这些传热管在流力弹性不安定性（Fluidelastic Instability）作用下振动，管与管之间碰撞、磨损而破损，会导致一级冷却水泄漏。美国在1980~1990年代的San Onofre核电站发生了大规模管束破损事故，电站最终被廃炉。该事故促使流力弹性不安定的设计基准（Connors式）进一步严格化，管束间隔和支撑结构的设计被从根本上重新审视。
管束流力振动的数值计算手法
数值手法的详细
具体用什么算法来求解管束流力振动解析呢？
离散化的定式化
使用形状函数 $N_i$ 来近似未知量:
用数式表达如下。
基础方程的离散形
用数式表达如下。
嗯，只看式子的话感觉模糊…到底表示什么呢？
将连续体的控制方程离散化后，得到如下代数方程组:
这里 $[K]$ 是全局刚性矩阵（或等效的系统矩阵），$\{u\}$ 是未知节点变量向量，$\{F\}$ 是外力向量。
哦，原来如此！将连续体的控制方程离散化这样就清楚了。
单元技术
"单元技术"这个名词听过，但可能理解得不够深…
| 单元类型 | 阶次 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低（剪切锁定） | 低 |
| 四面体2阶 | 二次 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中～高 | 中 |
积分方案
积分方案，具体是什么意思呢？
听到这里，终于明白为什么单元类型这么重要了！
收敛性和稳定性
不收敛的话，首先应该检查什么？
收敛速率: 二次单元时 $O(h^2)$ 阶的误差减小（光滑解的情况）
明白了…细分网格看起来简单，其实特别深奥呢。
求解器设置的推荐
具体用什么算法来求解管束流力振动解析呢？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数基准 |
| 预处理手法 | ILU(0) or AMG | 按问题规模 |
| 最大迭代次数 | 1000 | 不收敛时重新检查设置 |
| 内存模式 | In-core | 尽可能 |
单体耦合法
将所有物理场作为一个联立方程组同时求解。对于强耦合问题稳定，但实现复杂、内存消耗大。
分割法（分离迭代法）
各物理场独立求解，在界面处进行数据交换。实现容易、可重用现有求解器。适合弱耦合问题。
界面数据转移
最近邻法（最简单但精度低）、射影法（守恒）、RBF插值（对非匹配网格有利）。守恒性与精度的平衡很重要。
子迭代
在各耦合步骤内进行充分的迭代，确保界面条件的一致性。残差基准应根据各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛系数。防止过度松弛导致的发散，加快收敛的自适应手法。
稳定性条件
注意附加质量效应（流体-结构耦合中结构密度≈流体密度的情况）。如果不稳定，采用Robin型界面条件或IQN-ILS法。
管束流力振动的实务应用
实践指南
老师，请给我讲讲"实践指南"！
讲解管束流力振动解析的实务解析流程和注意事项。
老师的讲解容易理解！管束流力振动解析的困惑消散了。
解析流程
从最初的第一步教我！从什么开始比较好呢？
1. 预处理 (Pre-processing)
- CAD数据的导入与形状简化
- 材料特性的定义
- 网格生成（单元类型·尺寸的决定）
- 边界条件和荷载条件的设定
2. 求解 (Solving)
- 求解器设置（求解方法、收敛基准、输出控制）
- 作业提交和计算实行
- 收敛监视
3. 后处理 (Post-processing)
- 结果的可视化（位移、应力、其他物理量）
- 结果的验证与妥当性确认
- 报告作成
网格生成的最佳实践
网格的好坏怎样判断呢？
单元品质指标
请给我讲讲"单元品质指标"！
| 指标 | 理想值 | 容许范围 | 影响 |
|---|---|---|---|
| 纵横比 | 1.0 | < 5.0 | 精度下降 |
| 雅可比比 | 1.0 | > 0.3 | 单元退化 |
| 扭曲度 | 0° | < 15° | 精度下降 |
| 倾斜度 | 0° | < 45° | 收敛性恶化 |
| 锥形比 | 0 | < 0.5 | 精度下降 |
网格密度的决定
网格密度的决定，具体是什么意思呢？
边界条件的设置指南
听说边界条件这里出错的话全部都完蛋了…
哦，这样啊！过度约束注意这个意思了。
按商用工具分类的实现步骤
有各种各样的软件呢？ 各自的特点教一下！
| 工具名 | 开发方/现所属 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Mechanical（旧ANSYS Structural） | ANSYS Inc. | .cdb, .rst, .db, .ans, .mac |
| Ansys Fluent | ANSYS Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
COMSOL Multiphysics
请给我讲讲"COMSOL Multiphysics"！
1986年在瑞典创立。最初以MATLAB联动的FEMLAB开始，后来改名为COMSOL。多物理场领域有强项。
现所属: COMSOL AB
Ansys Mechanical（旧ANSYS Structural）
请给我讲讲"Ansys Mechanical"！
1970年由Swanson Analysis Systems Inc.（SASI）开发。基于APDL（Ansys Parametric Design Language）。
现所属: ANSYS Inc.
老师的讲解容易理解！工具名的困惑消散了。
常见失败与对策
初学者容易犯的错误有什么？想事先知道！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质差、边界条件不适当 | 网格改进、约束条件重新检查 |
| 应力异常大 | 应力特异点、网格相关 | 特异点回避、局部网格细化 |
| 位移不现实 | 材料常数错误、单位系不统一 | 确认输入数据 |
| 计算时间过长 | 不必要的细化、低效求解 | 网格最优化、并列计算 |
品质保证检查表
教科书上没有的"现场智慧"之类的东西有吗？
管束流力振动解析的整体面貌掌握了！ 明天开始在实务中多加留意。
好的，干得不错！ 实际动手操作是最好的学习方法。有不懂的地方随时来问。
化工厂热交换器更换——"以为是腐蚀，结果是振动"
实务中常见的管束振动漏诊的案例。在定期检查中发现传热管有伤痕时，通常最初判断为"腐蚀或冲蚀吧"。但当伤痕的分布集中在特定支撑位置，或管外周有规则的磨耗时，流力振动引起的"摩擦磨损"的可能性很高。如果发电厂运行中"为什么每5年同一地方的管子都损伤"这样的情况持续发生，首先用Connors式从流速、管配置、支撑间隔检查一下比较快。没有在设计阶段进行流力振动评价的旧式设备中，这类"谜样的重复损伤"并不少见。
管束流力振动的软件比较
商用工具比较
有各种各样的软件呢？ 各自的特点教一下！
对管束流力振动解析有对应的主要商用CAE工具的功能比较，以及各产品的历史背景进行详述。
支持工具一览
那么，管束流力振动解析可以用什么样的软件呢？
| 工具名 | 开发方/现所属 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Mechanical（旧ANSYS Structural） | ANSYS Inc. | .cdb, .rst, .db, .ans, .mac |
| Ansys Fluent | ANSYS Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
COMSOL Multiphysics
请给我讲讲"COMSOL Multiphysics"！
1986年在瑞典创立。最初以MATLAB联动的FEMLAB开始，后来改名为COMSOL。多物理场领域有强项。
现所属: COMSOL AB
Ansys Mechanical（旧ANSYS Structural）
请给我讲讲"Ansys Mechanical"！
1970年由Swanson Analysis Systems Inc.（SASI）开发。基于APDL（Ansys Parametric Design Language）。
现所属: ANSYS Inc.
Ansys Fluent
接下来是Ansys Fluent的话题。什么内容呢？
由Fluent Inc.开发。2006年被ANSYS收购。基于非结构网格的通用CFD求解器。
现所属: ANSYS Inc.
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题。什么内容呢？
由CD-adapco开发。2016年被Siemens收购，整合到Simcenter品牌。多面体网格是特征。
现所属: Siemens Digital Industries Software
哦，这样啊！在瑞典创立那个意思啊。
功能比较矩阵
预算和时间都有限，成本效益最佳是哪个呢？
| 功能 | COMSOL | Ansys Mechanical | Fluent | Star-CCM+ |
|---|---|---|---|---|
| 基本功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并列计算 | ○ | ○ | ○ | ○ |
| GPU对应 | △ | △ | △ | ○ |
转换时的风险
转换时的风险，具体是什么意思呢？
哦，这样啊！不同工具间的模型是这个意思啊。
许可形式
"许可形式"这个名词听过，但理解可能不够深…
| 工具 | 许可 | 特点 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 高额但有官方支持 |
| OpenFOAM | GPL | 免费但支持需付费 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的OSS求解器 |
选择指南
结果怎样才能判断该选哪个，判断标准教一下？
管束流力振动解析工具选择时要考虑以下因素:
管束流力振动解析的整体面貌掌握了！ 明天开始在实务中多加留意。
好的，干得不错！ 实际动手操作是最好的学习方法。有不懂的地方随时来问。
管束振动解析软件——原子能领域独特的生态系统
管束流力振动的解析工具有与通用CFD不同的生态系统。原子能领域中，法国EDF（电力公司）开发·维护的代码"Cast3m"和"Code_Aster"被用作欧洲标准。这些虽然开源，但比商用软件需要更多专门知识。日本方面，三菱重工、东芝能源、日立各自拥有独立的管束振动评估代码，并将验证实绩用于向规制当局（原子力规制委员会）的申请。另一方面，石化工厂中则使用HTRI（Heat Transfer Research Inc.）的Xist等工艺热交换器专用代码，内置有流力振动评估功能。有趣的是，不同行业的工具系谱完全不同。
管束流力振动的前沿研究
前沿主题与研究动向
管束流力振动解析领域今后要怎样演变呢？
讨论管束流力振动解析领域的最新研究动向和先进手法。
老师的讲解容易理解！管束流力振动解析的困惑消散了。
最新的数值手法
接下来是最新数值手法的话题。什么内容呢？
嗯，只看式子的话感觉模糊…到底表示什么呢？
高性能计算 (HPC) 的应对
| 并列化手法 | 概要 | 适用求解器 |
|---|---|---|
| MPI（领域分割） | 分布内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并列 | 多数求解器 |
| GPU（CUDA/OpenCL） | GPGPU活用。特别对显式法有效 | LS-DYNA、Fluent等 |
| 混合MPI+OpenMP | 节点间+节点内并列 | 大规模HPC环境 |
管束流力振动的故障排除
故障排除
也就是说在有限元法的地方马虎的话，后面就会吃亏啊。铭记于心！
常见错误与对策
老师也在管束流力振动解析里通宵调试过呢？（笑）
1. 收敛失败
收敛失败，具体是什么意思呢？
症状: 求解器在指定迭代次数内不收敛并异常结束
可能的原因:
- 网格品质不足（过度歪斜的单元）
- 材料参数设置不适当
- 初始条件不适当
- 非线性性太强（荷载步骤不足）