Ansys System Coupling 战略
Ansys System Coupling战略的理论基础
概要
老师！今天讲的是Ansys System Coupling战略，这是什么东西呢？
Ansys Workbench的系统耦合。Fluent-Mechanical间的双向连成。数据转送和收敛控制。
前辈说"系统耦合一定要好好做"，现在我理解了。
支配方程
离散化方法
这个方程具体怎样用计算机来求解呢？
利用有限元法（FEM）的空间离散化。组装单元刚度矩阵，构筑全体刚度方程。
矩阵求解算法
矩阵求解算法具体是什么意思呢？
用直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程。大规模问题中带预处理的迭代法效果显著。
| 求解法 | 分类 | 内存使用量 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小〜中规模 |
| Cholesky分解 | 直接法（对称正定） | O(n²) | 小〜中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说，有限元法环节偷工减料的话，以后就吃亏！牢记于心！
商用工具中的实现
做Ansys System Coupling战略用什么软件呢？
| 工具名 | 开发商/现属 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Mechanical (原ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
供应商谱系与产品整合经过
各个软件的历史背景很有趣吗？
COMSOL Multiphysics
请给我讲讲"COMSOL Multiphysics"!
1986年成立于瑞典。以MATLAB联动的FEMLAB起始，后改名为COMSOL。多物理场领域优势明显。
现属: COMSOL AB
Ansys Mechanical (原ANSYS Structural)
请给我讲讲"Ansys Mechanical"!
1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys参数化设计语言）。
现属: Ansys Inc.
Abaqus FEA (SIMULIA)
Abaqus FEA具体是什么意思呢？
1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购，整合到SIMULIA品牌。
现属: Dassault Systèmes SIMULIA
啊，这样啊！"年在瑞典成立"原来是这样的机制。
文件格式与互操作性
不同软件间传递数据时有什么需要注意的吗？
| 格式 | 扩展名 | 种类 | 概要 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | 符合ISO 10303的3D CAD数据交换格式。形状+PMI对应。 |
| IGES | .igs/.iges | 中立CAD | 初期CAD数据交换规范。曲面数据互相兼容有课题。逐步向STEP迁移。 |
| VTK | .vtk/.vtu | 可视化 | 可视化工具包格式。ParaView等使用。 |
| MED | .med | 网格/结果 | EDF/CEA开发。Code_Aster等使用。HDF5基础。 |
在不同求解器间转换模型时，必须注意单元类型的对应关系、材料模型的兼容性、荷载边界条件的表示差异。特别是高阶单元或特殊单元（内聚单元、用户定义单元等）往往无法在求解器间直接转换。
原来格式看似简单，实际上内涵深刻呢。
实务上的注意事项
教科书里没有的"现场智慧"之类的有吗？
网格收敛性确认、边界条件妥当性验证、材料参数敏感性分析非常重要。
好的，进展顺利！实际动手操作是最好的学习，有不懂的随时问我。
Ansys System Coupling的辅助求解器架构——为何需要"中介角色"
Ansys System Coupling能连接不同求解器（Fluent与Mechanical等），是因为有"System Coupling服务"这个中介进程。各求解器连接到System Coupling服务，进行数据收发。采用间接连接而非直接连接的设计有其道理：（1）数据补间可在服务侧集中管理，（2）时间步同步由中央控制，（3）某个求解器崩溃也能记录和重启。这是分布式计算的"代理模式"这一经典设计思想，与MpCCI采用相同哲学。
Ansys System Coupling战略的数值计算方法
数值方法的详细
具体用什么算法来求解Ansys System Coupling战略呢？
离散化的定式化
用形状函数 $N_i$ 近似未知量:
用式子表示就是这样。
基础方程的离散形式
用式子表示就是这样。
只有式子的话，我理解不了……这表示什么呢？
连续体的支配方程离散化后，得到如下代数方程组:
其中$[K]$是全体刚度矩阵（或等价系统矩阵），$\{u\}$是未知节点变量向量，$\{F\}$是外力向量。
啊，这样啊！连续体的支配方程经过什么处理变成这样的机制。
单元技术
听说过"单元技术"，但可能没完全理解……
| 单元类型 | 阶 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低（剪切锁定） | 低 |
| 四面体2阶 | 二次 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中〜高 | 中 |
积分方案
积分方案具体是什么意思呢？
听到这里，终于明白为什么单元类型那么重要了！
收敛性与稳定性
收敛不了的时候，首先检查什么？
收敛速度: 二次单元中，误差以$O(h^2)$阶减少（光滑解的情况）
原来网格细分看似简单，实际上内涵深刻呢。
求解器设置的建议
具体用什么算法来求解Ansys System Coupling战略呢？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数基准 |
| 预处理方法 | ILU(0) 或 AMG | 依问题规模而定 |
| 最大迭代次数 | 1000 | 不收敛时需重新设置 |
| 内存模式 | In-core | 尽可能使用 |
单片法
将全物理场作为1个联立方程系统同时求解。对强耦合问题稳定，但实现复杂、内存消耗大。
分割法（分离迭代法）
各物理场独立求解，在界面处交换数据。实现容易、可活用现有求解器。适于弱耦合。
界面数据转写
最近邻法（最简单但精度低）、投影法（保守）、RBF补间（对网格不匹配强）。保守性与精度的平衡很重要。
子迭代
各耦合步内进行充分迭代，确保界面条件的整合性。残差基准按各物理场的典型值缩放。
Aitken缓和
自动调整耦合迭代的缓和系数。防止过缓和导致的发散，加快收敛的自适应方法。
稳定性条件
注意added mass效应（流体-结构耦合中结构密度≈流体密度时）。不稳定情况下使用Robin型界面条件或IQN-ILS法。
Ansys System Coupling战略的实务应用
实践指南
老师，请给我讲讲"实践指南"！
解说Ansys System Coupling战略实务上的分析流程与注意事项。
也就是说在战略实务上的分析流程偷工减料，以后就吃亏！牢记于心！
分析流程
从最初的一步开始教我！应该从什么开始？
1. 前处理 (Pre-processing)
- CAD数据导入与形状简化
- 材料特性定义
- 网格生成（单元类型尺寸决定）
- 边界条件与荷载条件设置
2. 求解 (Solving)
- 求解器设置（求解法、收敛基准、输出控制）
- 作业投入与计算执行
- 收敛监控
3. 后处理 (Post-processing)
- 结果可视化（位移、应力、其他物理量）
- 结果验证与妥当性确认
- 报告制作
网格生成的最佳实践
怎样判断网格好坏？
单元质量指标
请给我讲讲"单元质量指标"！
| 指标 | 理想值 | 允许范围 | 影响 |
|---|---|---|---|
| 纵横比 | 1.0 | < 5.0 | 精度低下 |
| Jacobian比 | 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 |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
COMSOL Multiphysics
请给我讲讲"COMSOL Multiphysics"!
1986年成立于瑞典。以MATLAB联动的FEMLAB起始，后改名为COMSOL。多物理场领域优势明显。
现属: COMSOL AB
Ansys Mechanical (原ANSYS Structural)
请给我讲讲"Ansys Mechanical"!
1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys参数化设计语言）。
现属: Ansys Inc.
老师的说明清楚！工具名的疑惑消散了。
常见失败与对策
初学者容易犯的失败有哪些？想事先知道！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格质量不良、不适切边界条件 | 网格改进、拘束条件重新考虑 |
| 应力异常巨大 | 应力奇点、网格依赖 | 奇点回避、局部网格细分 |
| 位移非现实 | 材料常数错误、单位系不一致 | 输入数据确认 |
| 计算时间过长 | 不必要细分、求解法低效 | 网格优化、并行计算 |
品质保证检查表
教科书没有的"现场智慧"有吗？
好的，进展顺利！实际动手操作是最好的学习，有不懂的随时问我。
System Coupling的调试——日志文件看什么
Ansys System Coupling连成上不顺利时，调试的关键在于日志文件的读法。System Coupling执行中输出"scService.log"等日志，记录各时间步的转送量、收敛残差、各求解器的响应时间。要看的指标有三点：（1）转送变量的最大值最小值物理上妥当吗，（2）各时间步的收敛迭代次数是否有递增倾向（递增是不稳定的征兆），（3）特定求解器是否长期等待其他求解器（等待时间不均衡）。日志分析工具Ansys官方没提供，但社内工具中常用Python脚本对日志进行分析并可视化。
Ansys System Coupling战略的软件比较
商用工具比较
有各种软件呢？请分别讲讲特点！
对应Ansys System Coupling战略的主要商用CAE工具的功能比较与各产品的历史背景详述。
前辈说"对应战略的主要工具一定要好好做"，现在我理解了。
对应工具清单
做Ansys System Coupling战略用什么软件呢？
| 工具名 | 开发商/现属 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Mechanical (原ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
COMSOL Multiphysics
请给我讲讲"COMSOL Multiphysics"!
1986年成立于瑞典。以MATLAB联动的FEMLAB起始，后改名为COMSOL。多物理场领域优势明显。
现属: COMSOL AB
Ansys Mechanical (原ANSYS Structural)
请给我讲讲"Ansys Mechanical"!
1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys参数化设计语言）。
现属: Ansys Inc.
Abaqus FEA (SIMULIA)
Abaqus FEA具体是什么意思？
1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购，整合到SIMULIA品牌。
现属: Dassault Systèmes SIMULIA
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题。什么内容呢？
CD-adapco开发。2016年被Siemens收购整合到Simcenter品牌。多面体网格是特点。
现属: Siemens Digital Industries Software
原来"年在瑞典成立"看似简单，实际内涵深刻呢。
功能比较矩阵
预算时间都有限，成本效益最高是哪一个？
| 功能 | COMSOL | Ansys Mechanical | Abaqus | Star-CCM+ |
|---|---|---|---|---|
| 基本功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ | ○ |
| GPU对应 | △ | △ | △ | ○ |
转换时的风险
转换时的风险具体是什么意思？
啊，这样啊！不同工具间的模式转移原来是这样的机制。
许可证形式
听说过"许可证形式"，但可能没完全理解……
| 工具 | 许可证 | 特点 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 高额但官方支持 |
| OpenFOAM | GPL | 免费但支持有偿 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的OSS求解器 |
选择指南
最后应该选哪个，请给判断标准！
Ansys System Coupling战略工具选择应考虑以下因素:
好的，进展顺利！实际动手操作是最好的学习，有不懂的随时问我。
System Coupling与COMSOL——设计思想的根本差异
Ansys System Coupling与COMSOL多物理场，实现相同"连成分析"但设计思想截然相反。System Coupling是"既有专用求解器（Fluent、Mechanical等）通过消息传递连接的分布式架构"。COMSOL是"1个统一环境中多个物理场方程共存的单片式架构"。前者强点是"保持习惯求解器的精度功能"，后者强点是"物理场间的边界条件在方程级统一，理论一致性高"。哪个更优取决于问题，产业界通常"大规模FSI用System Coupling，科研教育用多物理场用COMSOL"这样分工。
Ansys System Coupling战略的前沿研究
前沿课题与研究动向
Ansys System Coupling战略这个领域今后怎样发展？
看Ansys System Coupling战略中的最新研究动向与先进方法。
等等，这样说的话也能用这样的情况吗？
最新的数值方法
接下来讲最新数值方法。什么内容？
只有式子的话，我理解不了……这表示什么呢？
高性能计算(HPC)的应对
| 并行化方法 | 概要 | 适用求解器 |
|---|---|---|
| MPI (域分割) | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 多数求解器 |
| GPU (CUDA/OpenCL) | GPGPU活用。特别对显式解法有效 | LS-DYNA, Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
Ansys System Coupling战略的故障排除
故障排除
常见错误与对策
老师也被Ansys System Coupling战略折腾过一整夜吗？（笑）
1. 收敛失败
收敛失败具体是什么意思？
症状: 求解器在指定迭代次数内不收敛，异常结束
考虑的原因:
- 网格质量不足（过度变形单元）
- 材料参数设置不恰当
- 初始条件不恰当
- 非线性性太强（荷载步不足）
对策:
- 实施网格质量检查（纵横比、Jacobian）
- 确认材料参数的单位系
- 将荷载分为多个步（增加子步数）
- 缓和收敛判定基准（但要注意精度）
也就是说收敛失败环节偷工减料，以后就吃亏！牢记于心！
2. 非物理结果
接下来讲非物理结果。什么内容？
症状: 应力/位移/温度等出现物理上非现实的值
考虑的原因:
- 边界条件设置错误
- 单位系混入（SI单位与工程单位混同）
- 不适切单元类型的选择
- 应力奇点的存在
对策:
- 确认反力合计（力的平衡）
- 确认单位系的一致性
- 重新检讨单元类型的适切性
- 奇点除去或子建模
前辈说"收敛失败一定要好好做"，现在我理解了。
3. 计算时间超过
计算时间超过具体是什么意思？
症状: 计算耗时预想时间数倍以上
对策:
- 网格粗密