强耦合与弱耦合
强耦合与弱耦合的理论基础
概述
老师！今天是关于强耦合和弱耦合的话题，对吧？这是什么东西呢？
按耦合强弱进行分类。当附加质量效应较大时强耦合是必需的。计算成本与精度的权衡。
控制方程
离散化方法
这个方程在计算机上实际是怎样解的呢？
使用有限元法（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 |
| 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。 |
在不同求解器之间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表达差异。特别是高阶单元或特殊单元（内聚单元、用户定义单元等）在求解器间通常无法直接转换。
我明白了…格式看似简单，实际上非常深奥。
实务注意事项
教科书里没有的"现场经验"之类的东西存在吗？
网格收敛性确认、边界条件合理性验证、材料参数敏感性分析都非常重要。
哇，强耦合和弱耦合的内容确实很深。不过在老师的讲解下已经理清了不少！
很好的样子！实际上亲手操作是最好的学习。有不懂的地方随时来问。
"耦合强度"的数学指标——算子分裂误差论
强耦合是否必需、是否弱耦合足够的量化定义其实并不简单。从理论上讲，应该用"算子分裂误差（splitting error）"来评估。弱耦合（算子分裂）中，"分别求解A和B→连接"的操作会产生O(Δt)的分裂误差。而强耦合则无此误差，耦合影响通过反复迭代（或同时求解）消除。哪个更重要取决于耦合系数的大小。"一个场的改变会导致另一个场产生巨大变化的问题"需要强耦合。热应变（温度→应变单向）可用弱耦合，但热生成依赖于变形的问题（如摩擦生热）需要强耦合。
强耦合与弱耦合的数值计算方法
数值方法详解
具体用什么算法来求解强耦合和弱耦合呢？
离散化公式化
用形状函数 $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 |
| 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.
老师的讲解很容易理解！工具名字的困惑都解开了。
常见失败与对策
初学者容易犯的错误有哪些？想提前知道！
| 现象 | 原因 | 对策 |
|---|---|---|
| 计算无法收敛 | 网格品质不良、边界条件不当 | 网格改善、约束条件重审 |
| 应力异常大 | 应力奇异点、网格依赖 | 避免奇异点、局部网格细化 |
| 位移非现实 | 材料常数错误、单位系不一致 | 确认输入数据 |
| 计算时间过长 | 不必要的细化、低效的求解方法 | 网格优化、并行计算 |
质量保证检查清单
教科书没有的"现场经验"之类的东西存在吗？
哇，强耦合和弱耦合的内容确实很深。不过在老师的讲解下已经理清了不少！
很好的样子！实际上亲手操作是最好的学习。有不懂的地方随时来问。
"弱耦合够不够"的判断基准——耦合敏感性分析的实务方法
判断是否需要强耦合还是弱耦合足够的实务方法是"耦合敏感性分析"。步骤是（1）用弱耦合计算，（2）固定或改变耦合变量（如温度），看相对场（如变形）变化多少，（3）相互影响小（如温度变化10%时应力仅变化1%）则弱耦合足够。实际上有个发动机热管理项目，弱耦合和强耦合的最大应力差仅0.3%，所以决定继续用弱耦合，计算时间减为1/3。这样的"正当偷工减料"是计算战略上重要的。
强耦合与弱耦合的软件比较
商用工具比较
有各种各样的软件吧？请告诉我各自的特点！
详细阐述支持强耦合和弱耦合的主要商用CAE工具的功能比较及各产品的历史背景。
支持工具列表
那么，要做强耦合和弱耦合分析用什么软件呢？
| 工具名称 | 开发商/现有公司 | 主要文件格式 |
|---|---|---|
| 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年被西门子收购，整合为Simcenter品牌。多面体网格为特色。
现在的所属：Siemens Digital Industries Software
那么…在瑞典成立看似简单，实际上特别深奥。
功能比较矩阵
预算和时间都有限的，性价比最强是哪个？
| 功能 | COMSOL | Ansys Mechanical | Abaqus | Star-CCM+ |
|---|---|---|---|---|
| 基础功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ | ○ |
| GPU支持 | △ | △ | △ | ○ |
转换时的风险
转换时的风险具体是什么意思呢？
啊，我明白了！不同工具之间的转就是这么一回事。
许可证形式
"许可证形式"听说过，但可能没有完全理解……
| 工具 | 许可证 | 特征 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 高价但有官方支持 |
| OpenFOAM | GPL | 免费但支持收费 |
| COMSOL | 节点锁定/浮动 | 按模块单独购买 |
| Code_Aster | GPL | EDF开发的OSS求解器 |
选择指南
最后到底该选哪个，请给我讲讲判断标准。
强耦合和弱耦合工具选择时应考虑以下因素：
哇，强耦合和弱耦合的内容确实很深。不过在老师的讲解下已经理清了不少！
很好的样子！实际上亲手操作是最好的学习。有不懂的地方随时来问。
COMSOL的"弱贡献（Weak Contribution）"功能——在方程中描述边界耦合
COMSOL有"弱贡献（Weak Contribution）"功能，可直接将界面物理场耦合写入偏微分方程弱形式。比如，流体-结构界面的力平衡可作为"对试函数的弱形式残差项"描述，这就直接作为耦合项加入方程。通过这样可用统一的方程体系实现单晶（强耦合）多物理场计算。对初学者有难度，但一旦理解，就可实现"COMSOL中几乎没有实现不了的耦合"这样的强大灵活性。适合自定义实现热-电-结构3场耦合等高度定制化需求。
强耦合与弱耦合的前沿研究
前沿话题与研究动向
强耦合和弱耦合这个领域以后会怎么发展呢？
看强耦合和弱耦合的最新研究动向和先进手法。
等等等等，强耦合和弱耦合在…就是说这样的情形下也能用吗？
最新数值方法
接下来是最新数值方法的话题吧。内容是什么呢？
只有公式的话我还是有点不太明白…这代表什么意思啊？
高性能计算（HPC）适配
| 并行化方法 | 概述 | 适用求解器 |
|---|---|---|
| MPI（域分割） | 分布式内存型。大规模问题的标准 | 全部主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 许多求解器 |
| GPU（CUDA/OpenCL） | GPGPU活用。特别对显式法有效 | LS-DYNA、Fluent等 |
| 混合MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
强耦合与弱耦合的故障排除
故障排除
常见错误与对策
老师也为强耦合和弱耦合的调试而通宵过吗？（笑）
1. 收敛失败
收敛失败具体是什么意思呢？
现象：求解器在指定迭代次数内无法收敛并异常结束
可能原因：
- 网格品质不足（过度畸变的单元）
- 材料参数设置不当
- 初始条件不合理
- 非线性过强（荷载步不足）
对策：
- 进行网格品质检查（纵横比、雅可比）
- 确认材料参数的单位系
- 分步加载（增加子步数）
- 放松收敛判定基准（但注意精度）
也就是说，在收敛失败这里马虎了，后来要吃苦头，对吧？铭记在心！