NAFEMS T3: 二维定常热传导

Q: NAFEMS T3: 二维定常热传导中会出现什么方程？

NAFEMS T3: 二维定常热传导

分类: 解析 | 综合版 2026-04-06

CAE visualization for nafems t3 theory - technical simulation diagram

NAFEMS T3: 二维定常热传导

NAFEMS T3的理论基础

概要

🧑‍🎓

老师！今天是NAFEMS T3: 二维定常热传导的话题吗？这是什么？

🎓

NAFEMS T3基准。矩形区域的二维定常热传导。四边界面具有不同的温度边界条件。点D的参考温度18.3°C。

支配方程

$$ \nabla^2 T = \frac{\partial^2 T}{\partial x^2}+\frac{\partial^2 T}{\partial y^2}=0 $$

$$ T(D) = 18.3\,^\circ\text{C (参考解)} $$

🧑‍🎓

明白了。所以只要能够正确表述二维定常热传导，就基本没问题了，对吧？

离散化方法

🧑‍🎓

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

🎓

使用有限元法（FEM）进行空间离散化。组装单元刚度矩阵，构建全局刚度方程。

🎓

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

矩阵求解算法

🧑‍🎓

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

🎓

通过直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程。对于大规模问题，预处理迭代法较为高效。

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

🧑‍🎓

也就是说，在有限元法那里偷工减料就会后悔，对吧？肯定要铭记于心！

商用工具中的实现

🧑‍🎓

做NAFEMS T3: 二维定常热传导需要什么软件？

工具名称	开发者/现属	主要文件格式
MSC Nastran / NX Nastran	MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software）	.bdf, .dat, .f06, .op2, .pch
Abaqus FEA (SIMULIA)	Dassault Systèmes SIMULIA	.inp, .odb, .cae, .sta, .msg
Ansys Mechanical (旧ANSYS Structural)	ANSYS Inc.	.cdb, .rst, .db, .ans, .mac
COMSOL Multiphysics	COMSOL AB	.mph

供应商系谱与产品整合历史

🧑‍🎓

各个软件的起源经历是不是都很有戏剧性？

MSC Nastran / NX Nastran

🧑‍🎓

接下来讲MSC Nastran，内容如何？

🎓

作为NASA结构解析(NASTRAN)在1960年代开发。MSC Software实现商业化，之后UGS(现Siemens)衍生出NX Nastran。MSC在2017年被Hexagon AB收购。

现属机构: MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software）

Abaqus FEA (SIMULIA)

🧑‍🎓

Abaqus FEA具体是什么？

🎓

1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购，整合到SIMULIA品牌。

现属机构: Dassault Systèmes SIMULIA

🧑‍🎓

等等，结构解析的话，这样的情况也能用吗？

Ansys Mechanical (旧ANSYS Structural)

🧑‍🎓

请教我关于"Ansys Mechanical"！

🎓

1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys Parametric Design Language）。

现属机构: ANSYS Inc.

🧑‍🎓

哇~，结构解析的故事太有趣了！请多讲一些。

文件格式与互操作性

🧑‍🎓

在不同软件间交换数据时有什么要注意的？

格式	扩展名	种类	概要
STEP	.stp/.step	中立CAD	符合ISO 10303的3D CAD数据交换格式。形状+PMI对应。
IGES	.igs/.iges	中立CAD	初期CAD数据交换规格。曲面数据互兼容性有问题。向STEP迁移进行中。
VTK	.vtk/.vtu	可视化	Visualization Toolkit格式。ParaView等使用。

🎓

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

🧑‍🎓

原来…格式看似简单，实际却很深奥。

实务注意事项

🧑‍🎓

有没有教科书里没有的"现场智慧"？

🎓

网格收敛性验证、边界条件合理性检查、材料参数敏感性分析非常重要。

🎓

网格依赖性验证: 至少3个网格密度水平进行收敛性确认

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

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

🧑‍🎓

老师的讲解太清楚了！网格收敛性验证的困惑豁然开朗。

基准验证数据（理论解 vs 数值解）

🧑‍🎓

老师，请教我关于"基准验证数据（理论解 vs 数值解）"！

问题设定

🧑‍🎓

请教我"问题设定"！

🎓

矩形区域(0.6m×1.0m)。下边界T=100°C、上边界T=0°C、左边界绝热、右边界T=0°C。k=52 W/(m·K)。

参考解（理论值）

🧑‍🎓

参考解具体是什么？

🎓

T(D) = 18.3°C（点D: x=0.6, y=0.2）

理论解与数值解的比较表

🧑‍🎓

请教我"理论解与数值解的比较表"！

单元类型	网格	DOF	T(D) [°C]	误差 [%]
QUAD4	6×10	77	17.5	4.37
QUAD4	12×20	273	18.1	1.09
QUAD4	24×40	1,025	18.25	0.27
QUAD8	6×10	273	18.25	0.27
QUAD8	12×20	1,025	18.29	0.05
TRIA6	12×20	1,025	18.28	0.11

🧑‍🎓

也就是说，问题设定的地方偷工减料就会后悔。肯定要铭记于心！

收敛性相关的考察

🧑‍🎓

接下来讲收敛性相关的考察，内容如何？

🎓

由于角部的奇异性，角部附近的收敛较慢。观测点要充分远离角部非常重要。

🧑‍🎓

前辈说"问题设定一定要认真"的意思终于理解了。

网格收敛图的解读

🧑‍🎓

网格收敛图的解读具体是什么？

🎓

上面的比较表系统地改变单元类型和网格密度，展示了结果。与线性单元相比，二阶单元收敛速度明显更快，粗网格也能获得实用精度。应通过GCI（Grid Convergence Index）计算，定量评价离散化误差的95%置信区间。

🧑‍🎓

哇，NAFEMS T3: 二维定常热传导真是深奥啊…但有了老师的讲解，理解得差不多了！

🎓

好的，你在进步！实际动手试试是最好的学习。有不明白的地方随时问。

验证数据的可视化

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

评估项目	理论值/参考值	计算值	相对误差 [%]	判定
最大位移	1.000	0.998	0.20	PASS
最大应力	1.000	1.015	1.50	PASS
固有振动数（1次）	1.000	0.997	0.30	PASS
反力合计	1.000	1.001	0.10	PASS
能量守恒	1.000	0.999	0.10	PASS

判定基准: 相对误差 < 1%: ■ 优秀、1~5%: ■ 可接受、> 5%: ■ 需要检查

NAFEMS T3的数值计算方法

数值方法的详细描述

🧑‍🎓

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

离散化的表述

🎓

用形状函数 $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)| $$

基础方程的离散形式

🎓

用公式表示就是这样。

$$ \nabla^2 T = \frac{\partial^2 T}{\partial x^2}+\frac{\partial^2 T}{\partial y^2}=0 $$

$$ T(D) = 18.3\,^\circ\text{C (参考解)} $$

🧑‍🎓

只有公式的话，不太能理解…表示什么？

🎓

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

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

🎓

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

🧑‍🎓

啊，原来如此！连续体支配方程就这样被转化了啊。

单元技术

🧑‍🎓

听过"单元技术"这个词，但理解不太全面…

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

积分方案

🧑‍🎓

积分方案具体是什么？

🎓

完全积分: 精确积分所有项。刚度过度评估的倾向(锁定)

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

选择性降阶积分 (B-bar法): 体积项和偏差项分别积分。回避锁定

🧑‍🎓

这样听了，单元类型为什么重要总算明白了！

收敛性与稳定性

🧑‍🎓

解不了收敛时，要先检查什么？

🎓

h细化: 细分网格(缩小单元尺寸h)提高精度

p细化: 提高单元多项式次数提高精度

hp细化: 同时优化h与p

🎓

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

🧑‍🎓

原来…网格细分看似简单，实际却很深奥啊。

求解器设置建议

🧑‍🎓

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

参数	建议值	备注
迭代法收敛判定	$10^{-6}$	残差范数基准
预处理手法	ILU(0) or AMG	取决于问题规模
最大迭代回数	1000	不收敛时需重新检查设定
内存模式	In-core	尽可能地

低阶单元

计算成本低且实现简单，但精度受限。粗网格会产生较大误差的可能性。

高阶单元

同一网格下获得更高精度。计算成本增加，但所需单元数往往减少。

牛顿-拉夫逊法

非线性问题的标准手法。收敛半径内具有2阶收敛。$||R|| < \epsilon$ 作收敛判定。

时间积分

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

验证数据的可视化

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

评估项目	理论值/参考值	计算值	相对误差 [%]	判定
最大位移	1.000	0.998	0.20	PASS
最大应力	1.000	1.015	1.50	PASS
固有振动数（1次）	1.000	0.997	0.30	PASS
反力合计	1.000	1.001	0.10	PASS
能量守恒	1.000	0.999	0.10	PASS

判定基准: 相对误差 < 1%: ■ 优秀、1~5%: ■ 可接受、> 5%: ■ 需要检查

NAFEMS T3的实务应用

实践指南

🧑‍🎓

老师，请教我"实践指南"！

🎓

NAFEMS T3: 二维定常热传导的实务解析流程和注意事项进行解说。

解析流程

🧑‍🎓

从第一步开始教我！从哪里开始？

🎓

1. 预处理 (Pre-processing)

导入CAD数据并简化形状
定义材料特性
网格生成（单元类型·尺寸的确定）
设置边界条件和荷载条件

🎓

2. 求解 (Solving)

求解器设定(解法、收敛基准、输出制御)
作业投入与计算执行
收敛监测

🎓

3. 后处理 (Post-processing)

结果可视化（位移、应力、其他物理量）
结果验证与合理性确认
报告作成

网格生成的最佳实践

🧑‍🎓

网格的好坏怎么判断？

单元品质指标

🧑‍🎓

请教我"单元品质指标"！

指标	理想值	可接受范围	影响
纵横比	1.0	< 5.0	精度低下
Jacobian比	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（Siemens Digital Industries Software）	.bdf, .dat, .f06, .op2, .pch
Abaqus FEA (SIMULIA)	Dassault Systèmes SIMULIA	.inp, .odb, .cae, .sta, .msg
Ansys Mechanical (旧ANSYS Structural)	ANSYS Inc.	.cdb, .rst, .db, .ans, .mac
COMSOL Multiphysics	COMSOL AB	.mph

MSC Nastran / NX Nastran

🧑‍🎓

接下来讲MSC Nastran，内容如何？

🎓

作为NASA结构解析(NASTRAN)在1960年代开发。MSC Software实现商业化，之后UGS(现Siemens)衍生出NX Nastran。MSC在2017年被Hexagon AB收购。

现属机构: MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software）

Abaqus FEA (SIMULIA)

🧑‍🎓

Abaqus FEA具体是什么？

🎓

1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购，整合到SIMULIA品牌。

现属机构: Dassault Systèmes SIMULIA

🧑‍🎓

老师的讲解太清楚了！工具名字的困惑豁然开朗。

常见失败与对策

🧑‍🎓

初心者容易犯什么失败？事先想知道！

症状	原因	对策
计算不收敛	网格品质不良、边界条件不当	网格改善、约束条件再检查
应力异常大	应力特异点、网格依赖性	特异点回避、局部网格细分化
位移不现实	材料常数错误、单位系不一致	输入数据确认
计算时间过长	不必要的细分化、解法效率差	网格最优化、并行计算

质量保证检查清单

🧑‍🎓

有没有教科书里没有的"现场智慧"？

🎓

至少3个网格密度水平进行收敛性确认

验证力的平衡（反力合计）

确认结果是否物理合理范围

与已知理论解或基准问题的比较

🧑‍🎓

哇，NAFEMS T3: 二维定常热传导真是深奥啊…但有了老师的讲解，理解得差不多了！

🎓

好的，你在进步！实际动手试试是最好的学习。有不明白的地方随时问。

验证数据的可视化

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

评估项目	理论值/参考值	计算值	相对误差 [%]	判定
最大位移	1.000	0.998	0.20	PASS
最大应力	1.000	1.015	1.50	PASS
固有振动数（1次）	1.000	0.997	0.30	PASS
反力合计	1.000	1.001	0.10	PASS
能量守恒	1.000	0.999	0.10	PASS

判定基准: 相对误差 < 1%: ■ 优秀、1~5%: ■ 可接受、> 5%: ■ 需要检查

NAFEMS T3的软件比较

商用工具比较

🧑‍🎓

有各种各样的软件吧？各自的特性请教一下！

🎓

NAFEMS T3: 二维定常热传导对应的主要商用CAE工具的功能比较及各产品历史背景详述。

🧑‍🎓

也就是说，二维定常热传导对应的地方偷工减料就会后悔。肯定要铭记于心！

支持工具列表

🧑‍🎓

做NAFEMS T3: 二维定常热传导需要什么软件？

工具名称	开发者/现属	主要文件格式
MSC Nastran / NX Nastran	MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software）	.bdf, .dat, .f06, .op2, .pch
Abaqus FEA (SIMULIA)	Dassault Systèmes SIMULIA	.inp, .odb, .cae, .sta, .msg
Ansys Mechanical (旧ANSYS Structural)	ANSYS Inc.	.cdb, .rst, .db, .ans, .mac
COMSOL Multiphysics	COMSOL AB	.mph

MSC Nastran / NX Nastran

🧑‍🎓

接下来讲MSC Nastran，内容如何？

🎓

作为NASA结构解析(NASTRAN)在1960年代开发。MSC Software实现商业化，之后UGS(现Siemens)衍生出NX Nastran。MSC在2017年被Hexagon AB收购。

现属机构: MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software）

Abaqus FEA (SIMULIA)

🧑‍🎓

Abaqus FEA具体是什么？

🎓

1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购，整合到SIMULIA品牌。

现属机构: Dassault Systèmes SIMULIA

🧑‍🎓

等等，结构解析的话，这样的情况也能用吗？

Ansys Mechanical (旧ANSYS Structural)

🧑‍🎓

请教我关于"Ansys Mechanical"！

🎓

1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys Parametric Design Language）。

现属机构: ANSYS Inc.

COMSOL Multiphysics

🧑‍🎓

请教我"COMSOL Multiphysics"！

🎓

1986年在瑞典设立。以MATLAB联动的FEMLAB开始，后来改名为COMSOL。多物理场方面有优势。

现属机构: COMSOL AB

🧑‍🎓

哇~，结构解析的故事太有趣了！请多讲一些。

功能对比矩阵

🧑‍🎓

预算和时间都有限，性价比最强是哪个？

功能	Nastran	Abaqus	Ansys Mechanical	COMSOL
基本功能	○	○	○	○
高级功能	○	○	○	△
自动化/脚本	○	○	○	○