MMS: 二维定常热传导

Q: 矩阵求解算法，具体是什么意思？

通过直接法(LU分解、Cholesky分解)或迭代法(CG法、GMRES法)求解联立方程。对于大规模问题，预处理迭代法很有效。

分类: 分析 | 统一版 2026-04-06

CAE visualization for mms heat 2d theory - technical simulation diagram

MMS: 二维定常热传导

MMS的理论基础

概要

🧑‍🎓

老师！今天是MMS: 二维定常热传导的话题，对吧？它是什么东西呢？

🎓

针对二维热传导方程的MMS。三角函数型制造解导出源项。二阶元素验证二阶收敛。

🧑‍🎓

明白了。那么如果对热传导方程的MMS掌握得好，基本上就没问题了，对吧？

支配方程

$$ T_{mms} = \sin(\pi x)\sin(\pi y) $$

$$ f = 2k\pi^2\sin(\pi x)\sin(\pi y) $$

🧑‍🎓

明白了。那么对二维定常热传导的描述，如果掌握得好，就基本没问题了，对吧？

离散化手法

🧑‍🎓

这个方程在计算机上实际如何求解呢？

🎓

使用有限元法(FEM)进行空间离散化。组立单元刚度矩阵，构建整体刚度方程。

🎓

进行弱形式（变分形式）转换，使用试验函数和形状函数，采用Galerkin方法进行表述。单元类型的选择（低阶元素 vs. 高阶单元、完全积分 vs. 减少积分)直接关系到解的精度和计算成本的权衡。

矩阵求解算法

🧑‍🎓

矩阵求解算法，具体是什么意思？

🎓

通过直接法(LU分解、Cholesky分解)或迭代法(CG法、GMRES法)求解联立方程。对大规模问题，带预处理的迭代法很有效。

求解法	分类	内存使用	适用规模
LU分解	直接法	O(n²)	小~中等规模
Cholesky分解	直接法(对称正定)	O(n²)	小~中等规模
PCG法	迭代法	O(n)	大规模
GMRES法	迭代法	O(n·m)	大规模·非对称
AMG预处理	预处理	O(n)	超大规模

🧑‍🎓

也就是说在有限元法这个地方偷工减料的话，以后就会吃大亏，对吧。我会牢记在心的！

商用工具中的实现

🧑‍🎓

那么，MMS: 二维定常热传导可以用什么软件呢？

工具名称	开发商/现有方	主要文件格式
MSC Nastran / NX Nastran	MSC Nastran(Hexagon)、NX Nastran(西门子数字工业软件)	.bdf, .dat, .f06, .op2, .pch
Abaqus FEA (SIMULIA)	达索系统 SIMULIA	.inp, .odb, .cae, .sta, .msg
ANSYS Mechanical (旧ANSYS Structural)	ANSYS Inc.	.cdb, .rst, .db, .ans, .mac
ANSYS Fluent	ANSYS Inc.	.cas, .dat, .msh, .jou
Simcenter STAR-CCM+	西门子数字工业软件	.sim, .java, .csv
COMSOL Multiphysics	COMSOL AB	.mph
OpenFOAM	开源(OpenCFD/ESI、OpenFOAM Foundation)	字典文件(blockMeshDict等), .foam

供应商系统和产品整合历史

🧑‍🎓

各个软件的发展历史还挺戏剧化的吗？

MSC Nastran / NX Nastran

🧑‍🎓

接下来是MSC Nastran的介绍吧。内容是什么？

🎓

NASA结构分析(NASTRAN)在1960年代开发。MSC Software将其商用化，后来UGS(现西门子)派生出NX Nastran。MSC在2017年被Hexagon AB收购。

现在的所属: MSC Nastran(Hexagon)、NX Nastran(西门子数字工业软件)

Abaqus FEA (SIMULIA)

🧑‍🎓

Abaqus FEA，具体是什么意思？

🎓

1978年由HKS (Hibbitt, Karlsson & Sorensen)开发。2005年被达索系统收购，整合到SIMULIA品牌。

现在的所属: 达索系统 SIMULIA

🧑‍🎓

等等等等，结构分析就是说，在这样的情况下也能用，对吧？

ANSYS Mechanical (旧ANSYS Structural)

🧑‍🎓

"ANSYS Mechanical"，告诉我一下！

🎓

1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL(ANSYS参数化设计语言)。

现在的所属: ANSYS Inc.

🧑‍🎓

哦~，结构分析的故事，超级有意思！再给我讲讲。

文件格式和互操作性

🧑‍🎓

不同软件之间传递数据的时候有什么要注意吗？

格式	扩展名	种类	概述
STEP	.stp/.step	中立CAD	ISO 10303兼容的三维CAD数据交换格式。支持形状+PMI。
IGES	.igs/.iges	中立CAD	早期CAD数据交换标准。曲面数据兼容性有问题。正在逐步过渡到STEP。
VTK	.vtk/.vtu	可视化	Visualization Toolkit格式。用于ParaView等。

🎓

在不同求解器之间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表达差异。特别是高阶单元或特殊单元(内聚单元、用户定义单元等)往往无法在求解器之间直接转换。

🧑‍🎓

明白了…格式看起来很简单，但实际上深不可测呢。

实务中的注意事项

🧑‍🎓

教科书里不会讲的"实战技巧"有什么吗？

🎓

网格收敛性的确认、边界条件的合理性验证、材料参数的敏感性分析特别重要。

🎓

网格依赖性验证: 至少用3个网格密度水平检验收敛性

边界条件的合理性: 设置有物理意义的约束条件

结果验证: 与理论解、实验数据、已知基准问题比较

🧑‍🎓

MMS: 二维定常热传导的整体框架我理解了！明天开始会在实务中注意这些。

🎓

很好，进展不错！实际动手去做才是最好的学习。有不懂的地方随时来问我。

验证数据可视化

定量显示理论值与计算值的比较。合格标准为误差5%以内。

评估项目	理论值/参考值	计算值	相对误差 [%]	判定
最大位移	1.000	0.998	0.20	通过
最大应力	1.000	1.015	1.50	通过
固有频率（第一阶）	1.000	0.997	0.30	通过
反力合计	1.000	1.001	0.10	通过
能量守恒	1.000	0.999	0.10	通过

判定基准: 相对误差 < 1%: ■ 优良、1~5%: ■ 允许、> 5%: ■ 需检讨

MMS的数值计算方法

数值方法详解

🧑‍🎓

具体用什么算法来求解MMS: 二维定常热传导？

离散化的表述

🎓

使用形状函数 $N_i$ 来近似未知量:

$$ u^h(\mathbf{x}) = \sum_{i=1}^{n} N_i(\mathbf{x}) \, u_i $$

🎓

用数式表示就是这样。

$$ K_e = \int_{\Omega_e} B^T \, D \, B \, d\Omega \approx \sum_{g=1}^{n_g} w_g \, B^T(\xi_g) \, D \, B(\xi_g) \, |J(\xi_g)| $$

基本方程的离散形式

🎓

用数式表示就是这样。

$$ T_{mms} = \sin(\pi x)\sin(\pi y) $$

$$ f = 2k\pi^2\sin(\pi x)\sin(\pi y) $$

🧑‍🎓

嗯…只有公式还是有点懵…表示什么呢？

🎓

连续体的支配方程离散化后，得到以下代数方程组:

$$ [K]\{u\} = \{F\} $$

🎓

其中$[K]$是整体刚度矩阵(或等价的系统矩阵)，$\{u\}$是未知节点变量向量，$\{F\}$是外力向量。

🧑‍🎓

啊，是这样！连续体的支配方程那样离散化啊，原来是那套机制。

单元技术

🧑‍🎓

"单元技术"听说过，但可能理解不透…

单元类型	阶数	节点数(3D)	精度	计算成本
四面体一阶	线性	4	低(剪切锁定)	低
四面体二阶	二次	10	高	中
六面体一阶	线性	8	中	中
六面体二阶	二次	20	非常高	高
棱柱	线性/二次	6/15	中~高	中

积分方案

🧑‍🎓

积分方案，具体是什么意思？

🎓

完全积分: 全部项准确积分。刚度过度评估的倾向(发生锁定)

减少积分: 削减积分点数。计算效率提高，但有沙漏模式发生的风险

选择性减少积分 (B-bar法): 将体积项和偏差项分开积分。规避锁定

🧑‍🎓

到这里为止听，终于明白为什么单元类型这么重要了！

收敛性和稳定性

🧑‍🎓

收敛不了的时候，首先要查什么？

🎓

h-加细: 细分网格(把单元尺寸 h 做小)来提高精度

p-加细: 提高单元多项式次数来提高精度

hp-加细: 同时优化 h 和 p

🎓

收敛速度: 二阶单元以 $O(h^2)$ 阶数衰减误差(光滑解的情况下)

🧑‍🎓

是啊…网格加细看起来简单，但实际上深不可测呢。

求解器设置建议

🧑‍🎓

具体用什么算法来求解MMS: 二维定常热传导？

参数	建议值	备注
迭代法收敛判定	$10^{-6}$	残差范数基准
预处理手法	ILU(0) or AMG	根据问题规模
最大迭代次数	1000	未收敛时需要重新设置
内存模式	核心内	尽可能使用

低阶单元

计算成本低，易于实现，但精度有限。粗网格下可能产生大误差。

高阶单元

在相同网格下实现更高精度。计算成本增加，但所需的单元数常常较少。

牛顿-拉夫逊法

非线性问题的标准方法。收敛半径内具有二次收敛。$||R|| < \epsilon$ 时判定收敛。

时间积分

显式方法: 条件稳定(满足CFL条件)。隐式方法: 无条件稳定，但每步需要求解联立方程。

验证数据可视化

定量显示理论值与计算值的比较。合格标准为误差5%以内。

评估项目	理论值/参考值	计算值	相对误差 [%]	判定
最大位移	1.000	0.998	0.20	通过
最大应力	1.000	1.015	1.50	通过
固有频率（第一阶）	1.000	0.997	0.30	通过
反力合计	1.000	1.001	0.10	通过
能量守恒	1.000	0.999	0.10	通过

判定基准: 相对误差 < 1%: ■ 优良、1~5%: ■ 允许、> 5%: ■ 需检讨

MMS的实务应用

实践指南

🧑‍🎓

老师，"实践指南"是什么？告诉我！

🎓

MMS: 二维定常热传导的实务分析流程和注意点。

分析流程

🧑‍🎓

从第一步开始告诉我！从什么开始？

🎓

1. 预处理 (Pre-processing)

导入CAD数据并简化形状
定义材料特性
网格生成(确定单元类型和大小)
设置边界条件和荷载

🎓

2. 求解 (Solving)

求解器设置(求解法、收敛基准、输出控制)
提交任务并执行计算
收敛监视

🎓

3. 后处理 (Post-processing)

结果可视化(位移、应力及其他物理量)
结果验证和合理性确认
制作报告

网格生成的最佳实践

🧑‍🎓

怎样判断网格的好坏？

单元品质指标

🧑‍🎓

"单元品质指标"是什么？告诉我！

指标	理想值	允许范围	影响
长宽比	1.0	< 5.0	精度降低
雅可比比值	1.0	> 0.3	单元退化
翘曲	0°	< 15°	精度降低
偏斜度	0°	< 45°	收敛性变差
锥形比	0	< 0.5	精度降低

网格密度的确定

🧑‍🎓

网格密度的确定，具体是什么意思？

🎓

应力集中部: 至少配置3层以上的单元

应力梯度大的区域: 把单元尺寸做成周围的1/3~1/5

荷载施加点近处: 局部加细

远方区域: 粗网格保证计算效率

边界条件设置指南

🧑‍🎓

边界条件啊，听说在这里出错的话全部都会报废…

🎓

注意过度拘束: 刚体运动拘束只需6个自由度

活用对称条件: 减小计算规模

荷载的等价分配: 集中荷载 vs. 分布荷载的选择

🧑‍🎓

啊，是这样！过度拘束要小心是那套机制。

按商用工具的实现步骤

🧑‍🎓

有各种各样的软件吧？各个的特点告诉我！

工具名称	开发商/现有方	主要文件格式
MSC Nastran / NX Nastran	MSC Nastran(Hexagon)、NX Nastran(西门子数字工业软件)	.bdf, .dat, .f06, .op2, .pch
Abaqus FEA (SIMULIA)	达索系统 SIMULIA	.inp, .odb, .cae, .sta, .msg
ANSYS Mechanical (旧ANSYS Structural)	ANSYS Inc.	.cdb, .rst, .db, .ans, .mac
ANSYS Fluent	ANSYS Inc.	.cas, .dat, .msh, .jou
Simcenter STAR-CCM+	西门子数字工业软件	.sim, .java, .csv
COMSOL Multiphysics	COMSOL AB	.mph
OpenFOAM	开源(OpenCFD/ESI、OpenFOAM Foundation)	字典文件(blockMeshDict等), .foam

MSC Nastran / NX Nastran

🧑‍🎓

接下来是MSC Nastran的介绍吧。内容是什么？

🎓

NASA结构分析(NASTRAN)在1960年代开发。MSC Software将其商用化，后来UGS(现西门子)派生出NX Nastran。MSC在2017年被Hexagon AB收购。

现在的所属: MSC Nastran(Hexagon)、NX Nastran(西门子数字工业软件)

Abaqus FEA (SIMULIA)

🧑‍🎓

Abaqus FEA，具体是什么意思？

🎓

1978年由HKS (Hibbitt, Karlsson & Sorensen)开发。2005年被达索系统收购，整合到SIMULIA品牌。

现在的所属: 达索系统 SIMULIA

🧑‍🎓

老师解释得清楚！工具名字的困惑消散了。

常见失败及对策

🧑‍🎓

初心者容易犯的错误有什么吗？想事先知道！

现象	原因	对策
计算不收敛	网格品质不良、不恰当的边界条件	改善网格、检查拘束条件
应力异常大	应力奇点、网格依赖	规避奇点、局部加细网格
位移不现实	材料常数错误、单位系不统一	确认输入数据
计算时间过长	不必要的加细、低效的求解法	优化网格、并行计算

质量保证检查清单

🧑‍🎓

教科书里不会讲的"实战技巧"有什么吗？

🎓

用3个以上网格密度水平检验收敛性了吗

验证了力的平衡(反力合计)吗

结果在物理上合理的范围内吗

与已知理论解或基准问题对比了吗

🧑‍🎓

MMS: 二维定常热传导的整体框架我理解了！明天开始会在实务中注意这些。

🎓

很好，进展不错！实际动手去做才是最好的学习。有不懂的地方随时来问我。

验证数据可视化

定量显示理论值与计算值的比较。合格标准为误差5%以内。

评估项目	理论值/参考值	计算值	相对误差 [%]	判定
最大位移	1.000	0.998	0.20	通过
最大应力	1.000	1.015	1.50	通过
固有频率（第一阶）	1.000	0.997	0.30	通过
反力合计	1.000	1.001	0.10	通过
能量守恒	1.000	0.999	0.10	通过

判定基准: 相对误差 < 1%: ■ 优良、1~5%: ■ 允许、> 5%: ■ 需检讨

MMS的软件比较

商用工具比较

🧑‍🎓

有各种各样的软件吧？各个的特点告诉我！

🎓

MMS: 二维定常热传导支持的主要商用CAE工具功能比较，以及各产品的历史背景详述。

🧑‍🎓

也就是说二维定常热传导这个地方偷工减料的话，以后就会吃大亏，对吧。肝铭心间！

支持工具列表

🧑‍🎓

那么，MMS: 二维定常热传导可以用什么软件呢？

工具名称	开发商/现有方	主要文件格式
MSC Nastran / NX Nastran	MSC Nastran(Hexagon)、NX Nastran(西门子数字工业软件)	.bdf, .dat, .f06, .op2, .pch
Abaqus FEA (SIMULIA)	达索系统 SIMULIA	.inp, .odb, .cae, .sta, .msg
ANSYS Mechanical (旧ANSYS Structural)	ANSYS Inc.	.cdb, .rst, .db, .ans, .mac
ANSYS Fluent	ANSYS Inc.	.cas, .dat, .msh, .jou
Simcenter STAR-CCM+	西门子数字工业软件	.sim, .java, .csv
COMSOL Multiphysics	COMSOL AB	.mph
OpenFOAM	开源(OpenCFD/ESI、OpenFOAM Foundation)	字典文件(blockMeshDict等), .foam

MSC Nastran / NX Nastran

🧑‍🎓

接下来是MSC Nastran的介绍吧。内容是什么？

🎓

NASA结构分析(NASTRAN)在1960年代开发。MSC Software将其商用化，后来UGS(现西门子)派生出NX Nastran。MSC在2017年被Hexagon AB收购。

现在的所属: MSC Nastran(Hexagon)、NX Nastran(西门子数字工业软件)

Abaqus FEA (SIMULIA)

🧑‍🎓

Abaqus FEA，具体是什么意思？

🎓

1978年由HKS (Hibbitt, Karlsson & Sorensen)开发。2005年被达索系统收购，整合到SIMULIA品牌。

现在的所属: 达索系统 SIMULIA

🧑‍🎓

等等等等，结构分析就是说，在这样的情况下也能用，对吧？

ANSYS Mechanical (旧ANSYS Structural)

🧑‍🎓

"ANSYS Mechanical"，告诉我一下！

🎓

1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL(ANSYS参数化设计语言)。

现在的所属: ANSYS Inc.

ANSYS Fluent

🧑‍🎓

接下来是ANSYS Fluent的介绍吧。内容是什么？

🎓

Fluent Inc.开发。2006年被ANSYS收购。基于非结构网格的通用CFD求解器。

现在的所属: ANSYS Inc.

🧑‍🎓