MpCCI 联合仿真
MpCCI 联合的理论基础
概述
老师！今天讲的是 MpCCI 联合仿真吧？这是什么东西呢？
来自 Fraunhofer SCAI 的 MpCCI 代码间联合。异构求解器组合。接口库。
控制方程
哇～联合仿真的话题超有意思！请继续讲！
离散化方法
这个方程实际上怎么用计算机求解呢？
使用有限元法（FEM）进行空间离散化。组装单元刚度矩阵，构建总刚度方程。
矩阵求解算法
矩阵求解算法具体是什么意思？
通过直接法（LU 分解、Cholesky 分解）或迭代法（CG 法、GMRES 法）求解联立方程组。大规模问题中采用预处理迭代法效果显著。
| 求解法 | 分类 | 内存使用量 | 适用规模 |
|---|---|---|---|
| LU 分解 | 直接法 | O(n²) | 小～中规模 |
| Cholesky 分解 | 直接法（对称正定） | O(n²) | 小～中规模 |
| PCG 法 | 迭代法 | O(n) | 大规模 |
| GMRES 法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG 预处理 | 预处理 | O(n) | 超大规模 |
也就是说有限元法的部分马虎不得，否则后续会很痛苦呢。我牢记在心！
商用工具中的实现
进行 MpCCI 联合仿真需要哪些软件呢？
| 工具名称 | 开发方/现属 | 主要文件格式 |
|---|---|---|
| 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。 |
在不同求解器间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载·边界条件的表达差异。特别是高阶单元或特殊单元（粘聚单元、用户定义单元等）在求解器间往往无法直接转换。
原来…格式看似简单，但实际上内涵很深呢。
实务中的注意事项
教科书里没有的「现场知识」是什么呢？
网格收敛性验证、边界条件有效性检查、材料参数的敏感性分析非常重要。
天哪，MpCCI 联合仿真真的很有意思呢……老师的讲解让我理解得更清楚了！
好样的！实际操作最重要。有问题随时问我。
MpCCI 的「代码间通信」——怎样连接不同求解器
MpCCI（基于网格的并行代码耦合接口）是一种实时在不同 CAE 软件间传输数据的中间件。从理论上讲，「参与求解器向 MpCCI 服务器发送数据，服务器进行插值·变换后转送到对方求解器」。Abaqus（结构）和 Fluent（流体）间的流固耦合 FSI 分析是典型例子，各求解器可专注于「自身计算」。开发单位是德国 Fraunhofer 研究所，在欧洲汽车、航空产业中广泛应用。
MpCCI 联合的数值计算方法
数值方法详情
具体用什么算法来求解 MpCCI 联合仿真？
哇～联合仿真的话题超有意思！请继续讲！
离散化的定式化
使用形状函数 $N_i$ 近似未知量：
用公式表示如下。
基本方程的离散形式
用公式表示如下。
嗯…仅看公式还是不太明白…这是什么意思呢？
连续体控制方程离散化后得到下列代数方程系：
其中 $[K]$ 是总体刚度矩阵（或同等系统矩阵），$\{u\}$ 是未知节点变量矢量，$\{F\}$ 是外力矢量。
啊！原来如此！连续体的控制方程就是这种机制呀。
单元技术
听说过「单元技术」，但理解得不够深…
| 单元类型 | 次数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体 1 阶 | 线性 | 4 | 低（剪切锁定） | 低 |
| 四面体 2 阶 | 二次 | 10 | 高 | 中 |
| 六面体 1 阶 | 线性 | 8 | 中 | 中 |
| 六面体 2 阶 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中～高 | 中 |
积分方案
积分方案具体是什么意思？
听到这儿，我终于明白单元类型为什么这么重要了！
收敛性与稳定性
收敛不了时，首先要检查什么？
收敛速度：二阶单元以 $O(h^2)$ 阶数减少误差（光滑解的情况下）
原来…网格细分看似简单，但实际上内涵很深呢。
求解器设置建议
具体用什么算法来求解 MpCCI 联合仿真？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数准则 |
| 预处理方法 | ILU(0) or AMG | 取决于问题规模 |
| 最大迭代次数 | 1000 | 未收敛时需重新调整设置 |
| 内存模式 | In-core | 尽可能采用 |
单体法
将所有物理场作为单个系统联立方程同时求解。对强耦合稳定，但实现复杂且内存消耗大。
分区法（分离迭代法）
各物理场依次求解，在界面处进行数据交换。实现简单，可复用已有求解器，适用于弱耦合。
界面数据转移
最近邻法（最简单但精度低）、投影法（保守性强）、RBF 插值（对网格不一致敏感性强）。需要平衡保守性和精度。
子迭代
在每个耦合时间步内进行充分迭代，确保界面条件的一致性。残差准则应按各物理场的典型值进行缩放。
Aitken 松弛
自动调整耦合迭代的松弛系数。防止过度松弛导致的发散，并加速收敛的自适应技术。
稳定性条件
注意 added mass 效应（流固耦合中结构密度≈流体密度的情况）。不稳定时可采用 Robin 型界面条件或 IQN-ILS 法。
MpCCI 联合的实务应用
实践指南
老师，请讲讲「实践指南」！
解说 MpCCI 联合仿真的实务分析流程和注意事项。
哇～联合仿真的话题超有意思！请继续讲！
分析流程
从第一步开始教我！需要从什么开始呢？
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 |
| 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.
老师的讲解很清楚！工具名称的疑惑解开了。
常见失败与对策
初学者常犯什么错误呢？想事先了解！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不好、边界条件不当 | 改善网格、重新检查约束条件 |
| 应力异常大 | 应力奇点、网格依存 | 避免奇点、局部网格细分 |
| 位移不现实 | 材料常数错误、单位系混乱 | 确认输入数据 |
| 计算时间过长 | 不必要的细分、低效求解 | 网格优化、并行计算 |
质量保证检查表
教科书里没有的「现场知识」是什么呢？
天哪，MpCCI 联合仿真真的很有意思呢……老师的讲解让我理解得更清楚了！
好样的！实际操作最重要。有问题随时问我。
MpCCI 联合的「计时配合」——同步点的设置是关键
在使用 MpCCI 的联合仿真中，何时交换数据（同步点）的设置很重要。如果每时间步都交换，精度高，但通信开销会让计算变慢。如果每 5 步、10 步进行一次，速度快，但精度下降。实务中「物理量变化大的阶段频繁交换，变化缓慢的阶段间隔交换」这种自适应设置很有效。有个案例是在风洞试验的 FSI 分析中，将同步间隔粗化 10 倍后，结果基本不变，计算时间缩短了 40%。
MpCCI 联合的软件比较
商用工具比较
有各种各样的软件吧？各自有什么特点？
MpCCI 联合仿真对应的主要商用 CAE 工具的功能对比和各产品的历史背景详述。
哇～联合仿真的话题超有意思！请继续讲！
支持工具列表
进行 MpCCI 联合仿真需要哪些软件呢？
| 工具名称 | 开发方/现属 | 主要文件格式 |
|---|---|---|
| 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)
现在讲 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 开发的开源求解器 |
选择指南
最后选哪个，判断标准是什么呢？
MpCCI 联合仿真工具选择时需考虑以下因素：
天哪，MpCCI 联合仿真真的很有意思呢……老师的讲解让我理解得更清楚了！
好样的！实际操作最重要。有问题随时问我。
MpCCI 的竞争对手——preCICE 的分工
与 MpCCI 类似的「联合中间件」还有开源软件 preCICE。MpCCI 以商用工具的支持和稳定性为卖点，preCICE 则与 OpenFOAM、FEniCS、SU2 等开源求解器的兼容性好，在学术研究社区快速普及。功能上 preCICE 积极实装最新的插值算法和拟 Newton 加速，在学术研究中表现先进。产业用途选 Abaqus/Fluent 联合时用安心的 MpCCI，学术用途追求自定义性时用 preCICE，这样的分工比较现实。
MpCCI 联合的先进研究
前沿话题与研究动向
MpCCI 联合仿真领域今后会如何发展呢？
MpCCI 联合仿真的最新研究动向和先进手法。
哇～联合仿真的话题超有意思！请继续讲！
最新数值方法
下面讲的是最新数值方法吧。怎样的内容呢？
嗯…仅看公式还是不太明白…这是什么意思呢？
高性能计算 (HPC) 的支持
| 并行化方法 | 概述 | 适用求解器 |
|---|---|---|
| MPI（域分割） | 分布内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 多数求解器 |
| GPU (CUDA/OpenCL) | GPGPU 活用。特别对显式法有效 | LS-DYNA、Fluent 等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模 HPC 环境 |
MpCCI 联合的故障排除
故障排除
哇～联合仿真的话题超有意思！请继续讲！
常见错误与对策
老师也在 MpCCI 联合仿真的调试中通宵过吗？（笑）
1. 收敛失败
收敛失败具体是什么意思？
症状：求解器在指定迭代次数内未收敛，异常终止
可能原因：
- 网格品质不足（过度歪斜的单元）
- 材料参数设置不当
- 初始条件不当
- 非线性性太强（缺乏荷载分步）
对策：
- 进行网格品质检查（纵横比、雅可比）
- 确认材料参数单位系
- 将荷载分为多个分步（增加子步数）
- 放宽收敛判定标准（但注意精度）
也就是说收敛失败的部分马虎不得，否则后续会很痛苦呢。我牢记在心！
2. 非物理结果
下面讲的是非物理结果的话题吧。怎样的内容呢？
症状：应力/位移/温度等在物理上不现实
可能原因：
- 边界条件设置错误
- 单位系混淆（SI 与工程单位混用）
- 不当的单元类型选择
- 应力奇点存在
对策：
- 检查反力合计（力的平衡）
- 确认单位系一致性
- 重新检讨单元类型的适切性
- 奇点消除或子模型分析
前辈说「收敛失败一定要好好处理」的意思我现在明白了。
3. 计算时间超过
计算时间的超过具体是什么意思？
症状：计算耗时是预期的好几倍
对策：
- 网