模型阶数约减方法
模型阶数约减方法的理论基础
概要
老师！今天是讲模型阶数约减方法的内容吧？那是什么东西呢？
大规模耦合模型的高效近似。参数化ROM。系统辨识混合方法。
前辈说\"大规模耦合模型的效果一定要做好\"，现在我明白那句话的意思了。
支配方程
啊，我明白了！模型阶数约减方法就是这样一种机制啊。
离散化方法
这个方程具体怎样在计算机上求解呢？
利用有限元法（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结构） | 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结构）
请给我介绍一下「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等。 |
在不同求解器间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表达差异。特别是高阶单元或特殊单元（内聚单元、用户定义单元等）往往无法在求解器间直接转换。
原来格式看似简单，但实际上隐藏着很深的内容啊。
实务注意事项
教科书上没有的\"现场智慧\"之类的东西有吗？
网格收敛性的确认、边界条件合理性的验证、材料参数的敏感性分析非常重要。
哇，模型阶数约减方法深度真的很深啊…不过在老师的讲解下，我已经整理清楚了！
嗯，学得不错！实际上动手做一下才是最好的学习。有不懂的地方随时来问吧。
POD和SVD——"只提取重要变动模式"的数学
作为模型阶数约减（MOR）基础的固有正交分解（POD）在数学上与奇异值分解（SVD）相同。从众多仿真快照中依次提取\"方差最大的方向（模态）\"。比如在热传导问题中收集100张温度场快照后进行POD，用上位5～10个模态往往能解释超过99%的整体变动。利用这一点，100万自由度的问题可被压缩成10自由度的问题。MOR研究前辈们将\"找到支配性模式\"的想法引入数值解析，由此打开了通向实时仿真的大门。
模型阶数约减方法的数值计算方法
数值方法详解
具体是用什么样的算法来求解模型阶数约减方法的呢？
离散化表述
用形状函数 $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. 前处理（预处理）
- CAD数据导入和形状简化
- 材料特性定义
- 网格生成（单元类型·尺寸的确定）
- 边界条件和荷载条件的设置
2. 求解（Solving）
- 求解器设置（解法、收敛基准、输出控制）
- 作业投入和计算执行
- 收敛监视
3. 后处理（后期处理）
- 结果可视化（位移、应力、其他物理量）
- 结果验证和合理性确认
- 报告编制
网格生成最佳实践
怎样判断网格的好坏呢？
单元质量指标
请给我介绍\"单元质量指标\"！
| 指标 | 理想值 | 许可范围 | 影响 |
|---|---|---|---|
| 纵横比 | 1.0 | < 5.0 | 精度下降 |
| 雅可比比 | 1.0 | > 0.3 | 单元退化 |
| 翘曲 | 0° | < 15° | 精度下降 |
| 斜度 | 0° | < 45° | 收敛性恶化 |
| 锥度比 | 0 | < 0.5 | 精度下降 |
网格密度的确定
网格密度的确定具体是怎么一回事呢？
边界条件设置指南
边界条件要是弄错了，整个就废掉了吧…
啊，我明白了！注意过拘束就是这样一种机制啊。
按商用工具分类的实现步骤
有很多软件吧？请给我讲讲各自的特点！
| 工具名称 | 开发方/现属 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| ANSYS Mechanical（原ANSYS结构） | 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结构）
请给我介绍一下「ANSYS Mechanical」！
1970年由Swanson Analysis Systems Inc.（SASI）开发。以APDL（ANSYS参数化设计语言）为基础。
现属：ANSYS Inc.
老师的说明听得很清楚！对工具名字的困惑消解了。
常见失败与对策
初学者容易犯什么样的错误？提前知道了可以避免！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格质量不良、不适当的边界条件 | 改进网格、重新检视拘束条件 |
| 应力异常大 | 应力奇点、网格依赖 | 避免奇点、局部网格细分 |
| 位移不现实 | 材料常数误差、单位系统不一致 | 确认输入数据 |
| 计算时间过长 | 不必要的细分、解法低效 | 网格优化、并行计算 |
质量保证检查清单
教科书上没有的\"现场智慧\"之类的东西有吗？
哇，模型阶数约减方法深度真的很深啊…不过在老师的讲解下，我已经整理清楚了！
嗯，学得不错！实际上动手做一下才是最好的学习。有不懂的地方随时来问吧。
MOR的\"离线/在线\"阶段——实务导入的流程
MOR的实务导入分为\"离线（事前计算）阶段\"和\"在线（预测）阶段\"两个步骤进行。离线阶段进行多数参数组合的全阶解析，用POD构建基底。这是耗时的投资阶段。在线阶段用构建的缩约模型进行新参数的高速预测。某飞机制造商构建翼的空气动力-热耦合ROM的事例中，离线用时2周，但在线预测每案件0.3秒，年间1000件以上的设计检讨中，总计算时间削减90%以上，这样的报告。
模型阶数约减方法的软件比较
商用工具比较
有很多软件吧？请给我讲讲各自的特点！
说明模型阶数约减方法所对应的主要商用CAE工具的功能比较、以及各产品的历史背景。
支持工具清单
那么，做模型阶数约减方法需要什么样的软件呢？
| 工具名称 | 开发方/现属 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| ANSYS Mechanical（原ANSYS结构） | 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结构）
请给我介绍一下「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开发的开源求解器 |
选择指南
最后到底该选哪个，能教我判断标准吗？
在模型阶数约减方法工具选择时需考虑以下因素：
哇，模型阶数约减方法深度真的很深啊…不过在老师的讲解下，我已经整理清楚了！
嗯，学得不错！实际上动手做一下才是最好的学习。有不懂的地方随时来问吧。
商用MOR工具的现状——ANSYS ROM Builder、Siemens、Modelon
作为商用工具提供MOR（模型阶数约减）的有ANSYS ROM Builder、Siemens Simcenter阶数约减建模、Modelon的Exportlib等。ANSYS ROM Builder可从FEA/CFD结果输出FMU（功能模拟单元）格式的ROM，用于与Simulink等控制系统工具的联动。Siemens版本以与Amesim的整合为特色。相比之下，学术界OpenROMM、FEniCS-ROMs等开源库活跃开发。商用工具优势在导入容易和专门支持，开源版本优势在最新算法试验实现速度快，这样的趋势。
模型阶数约减方法的前沿研究
前沿课题与研究动态
模型阶数约减方法领域今后会怎么发展呢？
来看模型阶数约减方法的最新研究动向和先进方法。
最新数值方法
接下来是最新数值方法的话题吧。是什么内容呢？
嗯，只看公式的话还是有点懵……那是表示什么呢？
高性能计算（HPC）适配
| 并行化方法 | 概述 | 适用求解器 |
|---|---|---|
| MPI（领域分割） | 分布式内存型。大规模问题标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 多数求解器 |
| GPU（CUDA/OpenCL） | GPGPU利用。特别在显式法有效 | LS-DYNA、Fluent等 |
| 混合MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
模型阶数约减方法的故障排除
故障排除
常见错误与对策
老师，您在模型阶数约减方法上也熬过夜进行调试吗？（笑）
1. 收敛失败
收敛失败具体是怎么一回事呢？
症状：求解器在指定迭代次数内未收敛，异常终止
可能的原因：
- 网格质量不足（过度扭曲的单元）
- 材料参数设置不当
- 初始条件不适当
- 非线性性过强（荷载步长不足）
对策：
- 进行网格质量检查（纵横比、雅可比比）
- 确认材料参数的单位系统
- 将荷载分成多个步长（增加子步数）
- 放宽收敛判定基准（但要注意精度）
也就是说，在收敛失败的地方偷工减料的话，后面就会吃大亏了。我牢记在心！
2. 非物理结果
接下来是非物理结果的话题吧。是什么内容呢？
症状：应力/位移/温度等为物理上非现实的值
可能的原因：
- 边界条件误设
- 单位系统混用（SI制与工程制混淆）
- 不适当的单元类型选择
- 应力奇点的存在
对策：
- 确认反力合计（力的平衡）
- 确认单位系统的一致性
- 重新检讨单元类型的适切性
- 消除或通过子建模处理奇点
前辈说过\"收敛失败就一定要好好做\"，现在我明白那句话的意思了。