海斯勒图模拟器 返回
传热学

海斯勒图模拟器 — 瞬态导热

用海斯勒图背后的一项近似,计算初始温度均匀的物体突然置于对流环境后的瞬态导热。调整毕渥数与傅里叶数,即可实时看到无限大平板与球体的中心温度以及任意位置的无量纲温度。

参数设置
物体形状
特征值方程与空间因子随形状而变
毕渥数 Bi
表面对流与内部导热之比 Bi = hL/k
傅里叶数 Fo
无量纲时间 Fo = αt/L²
无量纲位置
平板为 x/L,球为 r/r0(0=中心,1=表面)
计算结果
第一特征值 ζ₁
系数 C₁
中心无量纲温度 θ₀*
指定位置无量纲温度 θ*
毕渥数 Bi
傅里叶数 Fo
物体截面与温度分布 — Fo 扫描

在物体截面上叠加无量纲温度分布 θ*(位置)。中心变化最慢,对流表面(边缘)最快。颜色渐变表示温度高低,竖线标出当前位置。

中心温度 θ₀* 随傅里叶数 Fo 变化
温度分布 θ* 随位置变化
理论与主要公式

$$\theta_0^{ }=C_1\,e^{-\zeta_1^{2} Fo},\qquad \theta^{ }=\theta_0^{*}\cdot f(\zeta_1,\,\text{position})$$

中心无量纲温度 θ₀* 与指定位置的无量纲温度 θ*。C₁ 为系数,ζ₁ 为第一特征值,Fo 为傅里叶数。空间因子 f 在平板为 cos(ζ₁·x/L),在球为 sin(ζ₁·r/r0)/(ζ₁·r/r0)。

$$\zeta_1\tan\zeta_1=Bi \quad (\text{平板}),\qquad 1-\zeta_1\cot\zeta_1=Bi \quad (\text{球})$$

确定第一特征值 ζ₁ 的超越方程。毕渥数 Bi 表示表面对流与内部导热之比。一项近似仅用该第一项近似瞬态解,从而取代经典的海斯勒图。

什么是海斯勒图模拟器

🙋
传热课上出现了"海斯勒图",那种很难读的图到底是干什么用的?
🎓
简单说,它是用一张图就能读出"把热的东西突然冷却时,中心会是多少度"的工具。比如把淬火的钢棒丢进水里:表面立刻冷下来,中心却迟迟不降。海斯勒图就是让你不必计算,直接从图上读出"时间与中心温度的关系"。这是海斯勒在 1947 年整理出来的。
🙋
读图就是拿尺子描线对吧……说实话我总读不准。
🎓
对,这正是它的弱点。所以现在我们像这个工具一样用"一项近似"直接计算。瞬态导热的精确解其实是无穷级数,但当傅里叶数 Fo 超过 0.2 时,第二项及以后几乎消失。于是只保留第一项即可,这就是一项近似。海斯勒图不过是把第一项画成图,所以直接计算反而比读图更准,答案是一样的。
🙋
明白了。可是左边的"毕渥数"和"傅里叶数"有什么区别?看起来都是热的事。
🎓
好问题。可以把毕渥数 Bi 看作"空间之比",傅里叶数 Fo 看作"时间"。Bi = hL/k 是表面散热的难易(对流热阻)与内部传热的难易(导热热阻)之比。Bi 小则内部温度近乎均匀,Bi 大则中心与表面温差很大。Fo = αt/L² 是无量纲时间,Fo 增大,整个物体就逐渐接近环境温度。同时拖动两个滑块,就能看到上方图表的曲线形状随之变化。
🙋
还出现了"第一特征值 ζ₁",这又是什么?
🎓
ζ₁ 是决定温度分布"形状与衰减快慢"的数。对平板,它是超越方程 ζ₁·tan(ζ₁) = Bi 的解,手算解不出来,所以本工具用二分法数值求解。ζ₁ 越大,exp(−ζ₁²·Fo) 衰减越快,即冷得越快。增大 Bi 会使 ζ₁ 增大:表面散热越强,中心也冷得越快,这一物理规律已经写进了方程里。
🙋
那在实际工作里,用它能得到什么有用的信息呢?
🎓
比如"一块肉要加热几分钟中心才能达到安全温度""钢件淬火时中心达到规定温度要多少秒"这类问题。中心滞后最严重,所以只要中心满足条件,整个物体都满足。反过来,如果只想热处理表面,就看表面的 Fo。本工具同时输出中心 θ₀* 和任意位置 θ*,可以直接用于这类设计判断。

常见问题

海斯勒图是一种经典的列线图,用于读取初始温度均匀的物体突然置于对流环境后的瞬态导热。横轴为傅里叶数 Fo,纵轴为中心无量纲温度 θ0*,参数为毕渥数 Bi 的倒数。本工具直接计算该图背后的一项近似(仅保留级数第一项的解),用数值结果替代读图。
瞬态导热的精确解是无穷级数,但当傅里叶数 Fo 大于约 0.2 时,高阶模态迅速衰减,仅第一项即可近乎精确地表示温度,这就是一项近似。当 Fo 小于 0.2,即加热或冷却刚开始的初期阶段,高阶项的贡献不可忽略,因此本工具在 Fo<0.2 时给出提示。海斯勒图本身也以 Fo≳0.2 为使用前提。
毕渥数 Bi = hL/k 是物体内部导热热阻与表面对流热阻之比。Bi 较小(≲0.1)时内部温度近乎均匀,可用集总热容法;Bi 较大时内部会产生明显的温度梯度。傅里叶数 Fo = αt/L² 是无量纲时间,表示热量已渗透到物体内部的程度。Fo 越大,物体越接近稳态(与环境达到平衡)。
θ0* 是物体中心(平板的中心面、球的中心)的无量纲温度,由 θ0* = C1·exp(−ζ1²·Fo) 求得。θ* 是偏离中心的指定位置的无量纲温度,等于 θ0* 乘以与位置相关的空间因子。平板的空间因子为 cos(ζ1·x/L),球为 sin(ζ1·r/r0)/(ζ1·r/r0)。中心的温度变化最慢,对流表面最快接近环境温度。

实际应用

热处理与淬火:在钢材的淬火、退火、回火中,质量取决于部件中心是否达到规定温度,或冷却时中心是否越过马氏体转变温度。中心响应最慢,因此用海斯勒图(一项近似)估算中心温度的时间历程,从而确定保温时间与冷却速率。Bi 越大的急冷使中心跟随越快,但表面与中心的温差也越大,会引起热应力与开裂。

食品加热与杀菌:在罐头杀菌或烤制中,升温最慢的"冷点(中心附近)"必须达到安全温度或规定的 F 值。把食材近似为球或平板,再把傅里叶数换算为时间,即可预测中心达到条件的时刻。"只要中心满足条件,整体就满足"这一思路,是食品工程中常用的安全裕度取法。

电子元件与蓄热体的过渡响应:本工具可快速估算突然受到发热或环境温度变化的部件,需多少秒才能使内部温度达到平衡。Bi 较小(薄、高导热)的部件内部近乎均匀,接近集总热容法;Bi 较大的厚部件则需要考虑内部梯度的一项近似。

CAE 过渡传热分析的前期评估与验证:在运行有限元或有限体积的过渡导热分析之前,用一项近似可估出中心温度的数量级。若详细分析结果与该估算相差一个数量级,便可作为排查初始温度、对流边界条件或物性参数输入错误的合理性检查。反之,若与近似一致,则可对网格与时间步长的合理性更有信心。

常见误解与注意事项

最大的陷阱是"在 Fo 很小的初期阶段使用一项近似"。本工具的计算仅保留级数的第一项,仅当傅里叶数 Fo 大致在 0.2 及以上时才准确。当 Fo 更小,即冷却或加热刚开始时,第二、第三项等高阶模态仍然很大,仅用第一项会低估(或高估)温度。本工具在 Fo<0.2 时给出提示,但在该区域使用时应理解误差可能达到百分之几到十几个百分点。

其次,"弄错毕渥数的特征长度 L"。平板的 Bi 与 Fo 使用的 L 是"半厚"(板厚的一半),而不是板厚本身;球使用半径 r0。而集总热容法(Bi≲0.1 判据)使用的特征长度是"体积/表面积",对平板为半厚、对球为 r0/3,与海斯勒图所用的 L 定义不同。即便同名"毕渥数",特征长度取法不同就会使数值偏差,因此务必确认是哪种定义的 Bi。

最后,"θ* 不是实际温度本身,而是无量纲温度"。本工具输出的 θ0* 与 θ* 是无量纲量 θ* = (T−T∞)/(Ti−T∞):越接近 0 表示已到达环境温度 T∞,越接近 1 表示仍处于初始温度 Ti。要换回实际温度需用 T = T∞ + θ*·(Ti−T∞)。此外本工具假设物性(导热系数、热扩散率)恒定、初始温度均匀、对流换热系数恒定等理想条件。若存在相变、随温度变化的物性或非均匀初始条件,则超出适用范围,需要数值分析。

使用指南

  1. 输入毕渥数(Bi):根据物料导热系数λ、表面传热系数h和特征长度Lc计算Bi=hLc/λ,典型值范围0.01~100
  2. 输入傅里叶数(Fo):由时间t、物料导热系数α和特征长度平方Lc²确定,Fo=αt/Lc²,反映热扩散程度
  3. 设定位置系数(0~1):0表示物体中心,1表示表面;模拟器自动查表计算第一特征值ζ₁和系数C₁
  4. 输出中心温度θ₀*和指定位置温度θ*,即(T-T∞)/(T₀-T∞)的无量纲值

具体计算示例

混凝土厚板(Lc=0.1m,λ=1.4W/m·K,α=8.6×10⁻⁷m²/s)表面突然冷却,h=25W/m²·K。计算:Bi=25×0.1/1.4≈1.79;加热60分钟(t=3600s)后,Fo=8.6×10⁻⁷×3600/0.01≈0.31。在Bi=1.79、Fo=0.31处,查表得ζ₁≈1.52、C₁≈1.28,中心无量纲温度θ₀*≈0.68,距表面0.05m处(x/Lc=0.5)的θ*≈0.75,表示中心尚有68%的温差未消散。

实务注意事项

  1. 毕渥数判别物体内温度梯度:Bi<0.1时可采用集总热容法(温度均匀),Bi>40时主要受表面传热阻力控制
  2. 傅里叶数≥0.2表示瞬态阶段,≥2时逐渐接近准稳态;钢铁类物料(α≈2×10⁻⁵m²/s)升温快,陶土类(α≈6×10⁻⁷m²/s)升温慢
  3. 无限大平板模型适用于厚度远小于其他尺寸的零部件;球体模型用于球形铸件或球形传感器探头的温度响应预测
  4. 实际工程中应考虑物性随温度变化、接触热阻等因素修正结果