一维稳态热传导基准验证什么？

平板或杆的稳态热传导，把有精确解的温度分布和热流与求解器结果比较。改变条件（热导率恒定、有内热源、热导率随温度变化）来确认边界条件、物性输入和非线性迭代是否正确。正式的尺寸和条件请以NAFEMS原始文档为准。

一次单元的节点温度与精确解完全一致，精度够了吗？

一维一次单元有一个性质：热导率恒定时，即使有内热源，节点值也严格一致。算例中分2段节点误差就为0，但单元中点温度偏差6.25 K，热流偏差达最大值的25%。不能只比较节点，还要比较单元内部温度和热流。

热导率随温度变化时怎样验证？

用基尔霍夫变换可以构造精确解。热导率为温度一次式时，按单元平均温度计算恰好是精确的，所以要用二次式等条件才能暴露误差。算例中节点误差从2段0.091 K按二阶降到32段0.00049 K；牛顿法5次、逐次代入法8～14次收敛。

NAFEMS T1：一维稳态热传导基准

分类：V&V·NAFEMS基准 | 2026-10-01 改写

CAE visualization for nafems t1 theory - technical simulation diagram

NAFEMS T1: 一维热传导定常分析

理论：有精确解的一维问题

概述

🙋

被告知验证热分析求解器要从一维问题开始。是不是太简单了？

🎓

正因为简单，才能与精确解逐位比较。它能在不混入其他因素的情况下，确认边界条件的施加、热导率和热源的输入、物性随温度变化时的迭代计算是否正确。但有个陷阱：一维一次单元在某些条件下，节点温度与网格粗细无关地严格一致。只看节点就判断“对上了”，可能实际上什么都没验证。另外，NAFEMS热分析基准的正式尺寸和条件请以原始文档为准。

基本方程

$$ \frac{d}{dx}\left(k\,\frac{dT}{dx}\right) + \dot q = 0, \qquad T(0) = T_0,\ \ T(L) = T_L $$

热导率 $k$ 恒定且无热源 $\dot q$ 时温度为直线，有均匀热源时为抛物线。

随温度变化的热导率

$$ \theta(T) = \int_0^T \frac{k(s)}{k_0}\,ds \quad\Rightarrow\quad \frac{d^2\theta}{dx^2} = 0 \ \ (\dot q = 0), \qquad \theta \text{ is linear in } x $$

用基尔霍夫变换 $\theta$，即使热导率随温度变化，$\theta$ 也是 $x$ 的一次式，可通过逆变换求精确解——这是验证非线性热分析的经典方法。

Coffee Break 随笔闲谈

NAFEMS与基准题

NAFEMS于1983年在英国成立，旨在提高有限元分析的质量和可靠性。当时商用有限元软件开始普及，同一问题不同软件给出的答案却可能不同。于是NAFEMS收集了有精确解或可靠参考解的标准问题，按线弹性（LE）、热（T）、固有值（FV）等领域作为基准题发布。至今许多软件的验证手册仍刊载这些问题的结果。

算例：用一次单元求解

厚0.1 m的平板，一侧100 ℃、另一侧0 ℃，热导率50 W/mK，用一次单元求解。

例1：均匀热源（10⁶ W/m³）

分段数	节点误差	单元中点误差	端部热流误差
2	0（精确）	6.25 K	最大值的25%
4	0（精确）	1.56 K	12.5%
8	0（精确）	0.39 K	6.2%
16	0（精确）	0.10 K	3.1%

🙋

节点只分2段就误差为0，单元内部却偏了6 K？

🎓

一维一次单元在热导率恒定时节点值是精确的。但单元内部只是用直线连接，与抛物线温度分布有偏差；热流在每个单元内恒定，所以端部偏差达最大值的25%。而且误差减小方式也不同：单元中点温度为二阶，热流为一阶。只用节点温度验证，这些误差完全看不见。比较热流和单元内部温度很重要。

例2：随温度变化的热导率

热导率取 $k = k_0(1 + \beta T)$（$\beta$ = 0.005 1/K）时，中点精确解为54.95 ℃（恒定时为50 ℃）。但按单元平均温度计算热导率时，任何分段数节点都严格一致，因为对一次式热导率，这种计算在数学上与基尔霍夫变换相同。于是改用 $k = k_0(1 + \beta T)^2$ 比较（中点精确解59.62 ℃）：

分段数	节点最大误差	收敛阶	牛顿法	逐次代入法
2	0.091 K	—	5次	8次
4	0.031 K	1.58	5次	12次
8	0.0078 K	1.97	5次	13次
16	0.0020 K	1.99	5次	14次
32	0.00049 K	2.00	5次	14次

误差按二阶减小，牛顿法与分段数无关地5次收敛。选择验证问题时，关键是选能真正暴露误差的条件。

验证步骤

在热导率恒定、无热源条件下，确认边界条件和单位正确。
加入热源，不仅比较节点，还要把单元内部温度和热流与精确解比较。
引入随温度变化的热导率，与基尔霍夫变换精确解比较，避免恰好精确的条件。
加密网格，确认收敛阶。
记录非线性迭代次数和收敛判据。

Coffee Break 随笔闲谈

“基准通过，实物却不符”

某公司内部求解器验证通过了一维热传导基准，但实际零件分析中热流与测量偏差很大。验证只比较了节点温度，而那个条件下一次单元在节点上恰好精确。实际零件关心的是热源集中处的热流，单元太粗，热流误差大。如果验证阶段也比较热流，就能发现网格不够细。比较什么，要与实际想知道的量相匹配。

常见错误

错误与对策

错误	后果	对策
只比较节点温度	看不到误差	也比较热流和单元内部
使用恰好精确的条件	起不到验证作用	选择会出现误差的条件
只用一套网格判断	不知道收敛阶	多套网格比较
非线性收敛判据太松	残留迭代误差	收紧判据比较
不确认原始条件	与别的问题比较	确认原始尺寸与条件

🙋

我想了解相关内容。

🎓

可以用Richardson外推工具试算收敛阶。相关页面参考T1验证结果、NAFEMS T3、NAFEMS T4、NAFEMS LE1和一维热传导解析解。