多尺度耦合分析
多尺度耦合的理论基础
概述
老师！今天要讲多尺度耦合分析对吧？那是什么呢？
微观结构与宏观结构的尺度间耦合。计算均质化方法(FE²法)。复合材料、多孔材料的有效特性。
前辈说过"微观结构和宏观结构之间的这个部分一定要弄清楚"，现在明白了。
支配方程
等等等等，多尺度耦合分析在这样的情况下也能用吗？
离散化手法
具体怎样在计算机上求解这个方程呢？
使用有限元法(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年在瑞典成立。最初以FEMLAB(MATLAB联动)的名义推出，后改名为COMSOL。在多物理场领域具有优势。
现在的母公司: COMSOL AB
Ansys Mechanical (旧ANSYS Structural)
请为我介绍一下"Ansys Mechanical"！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL(Ansys Parametric Design Language)。
现在的母公司: 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 | 可视化 | Visualization Toolkit格式。ParaView等中使用。 |
在不同求解器间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表示差异。特别是对于高阶单元或特殊单元(内聚单元、用户定义单元等)，求解器间通常无法直接转换。
明白了…格式看似简单，但其实非常深奥呢。
实务中的注意事项
有没有教科书上找不到的那种"现场的智慧"呢？
网格收敛性的确认、边界条件的合理性验证、材料参数的灵敏度分析非常重要。
多尺度耦合分析的全貌掌握了！明天开始在实务中就会着重关注这些。
是的，加油！亲身动手是最好的学习。有问题的时候随时来问我。
均质化理论——通过"平均化"从微观到宏观
多尺度耦合的理论基础之一是均质化(Homogenization)理论。对于微观上不均质的材料(如碳纤维强化塑料的纤维-基体结构)，设置代表体积单元(RVE)，从微观计算得到等效的宏观物性。FE²法中，每个宏观积分点都对应一个RVE，进行微观计算，实现了"设计参数改变时物性自动通过RVE计算更新"的机制。但是当宏观积分点达到1万个，各RVE需要数千个单元时，计算成本会非常巨大，必须依靠高性能HPC环境。
多尺度耦合的数值计算方法
数值手法的详细内容
具体用什么算法来求解多尺度耦合分析呢？
那也就是说，多尺度耦合分析做得好的话，基本上就没问题了是吧？
离散化的定式化
使用形状函数 $N_i$ 进行未知量的逼近:
用公式表示就是这样。
基本方程的离散形式
用公式表示就是这样。
嗯，光看公式的话不太明白…这表示什么呢？
将连续体的控制方程离散化后，得到如下代数方程组:
其中 $[K]$ 是整体刚度矩阵(或等效的系统矩阵)，$\{u\}$ 是未知节点变量矢量，$\{F\}$ 是外力矢量。
啊，原来是这样啊！连续体的控制方程就这样处理啊。
单元技术
"单元技术"听过，但可能没有深入理解…
| 单元类型 | 阶数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低(剪切锁定) | 低 |
| 四面体2阶 | 二次 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二次 | 20 | 非常高 | 高 |
| 棱柱形 | 线性/二次 | 6/15 | 中~高 | 中 |
积分方案
积分方案具体是什么意思？
听到这里，我总算理解单元类型为什么这么重要了！
收敛性和稳定性
不收敛的时候，首先要检查什么呢？
收敛速度: 二阶单元的误差以 $O(h^2)$ 的阶数减小(光滑解的情况)
明白了…网格细化看似简单，但其实非常深奥呢。
求解器设置的建议
具体用什么算法来求解多尺度耦合分析呢？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法的收敛判定 | $10^{-6}$ | 残差范数基准 |
| 预处理方法 | ILU(0) 或 AMG | 根据问题规模 |
| 最大迭代次数 | 1000 | 不收敛时需重新审视设置 |
| 内存模式 | In-core | 尽可能使用 |
单体法
将所有物理场作为1个联立方程组同时求解。对于强耦合能够保证稳定，但实现复杂，内存消耗大。
分割法(分离迭代法)
各物理场独立求解，在界面处进行数据交换。实现简单，可利用现有求解器。适合弱耦合。
界面数据转移
最近邻法(最简单但精度低)、投影法(保守)、RBF插值(对网格不匹配的鲁棒性强)。保守性和精度的平衡很重要。
子迭代
在每个耦合步内进行充分的迭代，确保界面条件的协调性。残差基准应按各物理场的典型值进行缩放。
Aitken加速
自动调整耦合迭代的加速系数。防止过加速引起的发散，加快收敛的自适应方法。
稳定性条件
注意流固耦合中的added mass效应(结构密度≈流体密度的情况)。不稳定时可采用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年在瑞典成立。最初以FEMLAB(MATLAB联动)的名义推出，后改名为COMSOL。在多物理场领域具有优势。
现在的母公司: COMSOL AB
Ansys Mechanical (旧ANSYS Structural)
请为我介绍一下"Ansys Mechanical"！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL(Ansys Parametric Design Language)。
现在的母公司: ANSYS Inc.
老师的讲解很清楚！对工具名称的迷茫消散了。
常见的失败及对策
初学者容易犯什么样的错误呢？事先想了解一下！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、边界条件不当 | 改善网格、重新审视约束条件 |
| 应力异常大 | 应力特异点、网格依赖 | 避免特异点、局部网格细化 |
| 位移非现实 | 材料常数错误、单位系不一致 | 确认输入数据 |
| 计算时间过长 | 不必要的细化、低效的求解法 | 网格优化、并行计算 |
质量保证清单
有没有教科书上找不到的那种"现场的智慧"呢？
多尺度耦合分析的全貌掌握了！明天开始在实务中就会着重关注这些。
是的，加油！亲身动手是最好的学习。有问题的时候随时来问我。
复合材料的多尺度分析——飞行部件的疲劳寿命预测
多尺度耦合最为实用化的领域是碳纤维强化复合材料(CFRP)的疲劳寿命预测。在宏观尺度求解整个翼的应力分布，将疲劳临界位置的应力传送给微观尺度RVE分析。RVE中追踪纤维-基体界面的损伤进展，并将刚度降低反馈到宏观模型。在空客A350和波音787的设计中，这样的多尺度方法被用作材料试验的替代，试样数量削减了传统水平的40~60%。监管机构(FAA等)对仿真认可的信任度也在逐步提高。
多尺度耦合的软件比较
商用工具比较
有很多不同的软件吧？各自有什么特点呢？
详述支持多尺度耦合分析的主要商用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年在瑞典成立。最初以FEMLAB(MATLAB联动)的名义推出，后改名为COMSOL。在多物理场领域具有优势。
现在的母公司: COMSOL AB
Ansys Mechanical (旧ANSYS Structural)
请为我介绍一下"Ansys Mechanical"！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL(Ansys Parametric Design Language)。
现在的母公司: 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开发的开源求解器 |
选型指南
最后要选哪个，有没有判断标准呢？
在选择多尺度耦合分析工具时应考虑:
多尺度耦合分析的全貌掌握了！明天开始在实务中就会着重关注这些。
是的，加油！亲身动手是最好的学习。有问题的时候随时来问我。
Digimat、SwiftComp、MTEX——多尺度分析工具的选法
在复合材料多尺度分析的商用工具中，代表的是e-Xstream公司(Hexagon集团)的Digimat。从纤维配向数据(注塑模流分析结果)自动生成RVE，将物性映射到Ansys/Abaqus的用法最为通用。SwiftComp(普渡大学的工具)在梁-壳单元的均质化上有优势。在开源中，MTEX(MATLAB工具箱)在晶体塑性的多尺度分析中被学界广泛采用。工具选择的关键是"与想要使用的宏观求解器的联动是否建立""是否支持必要的损伤模型"这两点最重要。
多尺度耦合的前沿研究
前沿课题与研究动向
多尺度耦合分析这个领域，以后会怎么发展呢？
观察多尺度耦合分析领域的最新研究动向和先进手法。
那也就是说，多尺度耦合分析做得好的话，基本上就没问题了是吧？
最新的数值方法
接下来讲最新数值方法吧。有什么内容吗？
光看公式的话不太明白…这表示什么呢？
高性能计算 (HPC) 的适应
| 并行化手法 | 概述 | 适用求解器 |
|---|---|---|
| MPI (领域分割) | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 许多求解器 |
| GPU (CUDA/OpenCL) | GPGPU利用。特别对显式法有效 | LS-DYNA, Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
多尺度耦合的故障排查
故障排查
那也就是说，多尺度耦合分析做得好的话，基本上就没问题了是吧？
常见错误及对策
老师也在多尺度耦合分析上搞过通宵调试吗？(笑)
1. 收敛失败
收敛失败具体是什么意思？
症状: 求解器在指定的迭代次数内不收敛，异常终止
可能的原因:
- 网格品质不足(过度歪斜的单元)
- 材料参数设置不当
- 初始条件不当
- 非线性太强(荷载步数不足)
对策:
- 进行网格品质检查(宽高比、雅可比行列式)
- 确认材料参数的单位系
- 将荷载分成多个步(增加子步数)
- 放宽收敛判定基准(但要注意精度)
也就是说收敛失败的地方偷工减料的话，以后就要吃大苦头。牢牢记住！
2. 非物理结果
接下来讲非物理结果吧。有什么内容吗？
症状: 应力/位移/温度等非现实值
可能的原因:
- 边界条件设置错误
- 单位系混淆(SI单位与工程单位混用)
- 单元类型选择不当
- 应力特异点的存在
对策:
- 验证反力合计(力的平衡)
- 确认单位系的一致性
- 重新检