集总热容法模拟器 返回
Interactive Tool — Transient Heat Transfer

集总热容法·冷却/加热曲线模拟器

在确认毕奥数(Bi)有效性的同时,实时可视化非稳态导热的温度历程。可对比钢、铝、陶瓷等各种材料的时间常数。

参数设置
物体形状
材料预设
密度 ρ [kg/m³]
kg/m³
比热容 c_p [J/kgK]
J/kg·K
导热系数 k [W/mK]
W/m·K
对流系数 h [W/m²K]
W/m²K
特征长度 L_c [mm]
mm
初始温度 T_i [°C]
°C
流体温度 T_∞ [°C]
°C
计算时间 t_max [s]
s
计算结果
毕奥数 Bi
有效 (Bi < 0.1)
时间常数 τ
s
T at t=τ
°C
99%平衡到达时间
s
集总热容法不适用(Bi > 0.1)— 需要一维解析
温度 T(t) — 冷却/加热曲线(τ标记)
无量纲温度 θ*(t) — 半对数坐标
毕奥数指示器
理论与关键公式

$$T(t) = T_\infty + (T_i - T_\infty)\,e^{-t/\tau}$$

集总热容法解:\(\tau = \rho c_p V / (hA_s)\),表示热响应时间常数 [s]

$$\mathrm{Bi} = \frac{hL_c}{k} \leq 0.1$$

适用条件:\(L_c = V/A_s\) 为特征长度,\(k\) 为固体导热系数 [W/(m·K)]

什么是集总热容法

🙋
「集总热容法」是什么?听起来好复杂。
🎓
简单来说,就是把整个物体当成一个点,认为它内部温度瞬间就变得一样了。这就像你刚泡好一杯茶,觉得整个茶杯都烫手,而不会去想茶杯壁和茶杯把手的温度是不是不同。在实际工程中,当物体导热很快、表面散热相对较慢时,这个假设就成立。你可以在模拟器里选个「铝」试试,它的导热快,就挺符合这个条件。
🙋
诶,真的吗?那怎么判断什么时候能用这个“当成一个点”的方法呢?
🎓
这就靠一个叫「毕奥数(Bi)」的无量纲数来判断了。简单说,如果物体内部的导热阻力远小于表面的对流散热阻力,Bi数就很小(通常小于0.1),就能用集总热容法。你试着在模拟器里把「对流系数 h」的滑块拖到很大,或者把「特征长度 L_c」拖长,看看上面的Bi数会不会变红超过0.1,那就是在提醒你“这个方法可能不准了哦”。
🙋
原来那个Bi数是干这个用的!那温度具体怎么算呢?比如一块烧红的钢块扔进水里,多久会凉?
🎓
问得好!这就是模拟器核心要算的「冷却曲线」。温度随时间是指数下降(或上升)的,公式就在下面。关键是一个叫「时间常数 τ」的参数,它代表了物体热响应的快慢。τ越大,凉得越慢。你可以在模拟器里对比一下「钢」和「铝」,把其他参数设成一样,你会发现铝的τ小得多,曲线掉得飞快,这就是为什么铝锅热得快凉得也快!

物理模型与关键公式

集总热容法的核心控制方程,来源于物体的能量平衡:物体内能减少的速率等于表面通过对流散失热量的速率。

$$ \rho c_p V \frac{dT}{dt}= -h A_s (T - T_{\infty}) $$
ρ 密度 kg/m³ V 体积 c_p 比热容 J/(kg·K) T 物体温度 °C 或 K t 时间 s h 对流换热系数 W/(m²·K) A_s 表面积 T_∞ 流体温度 °C 或 K

对上述方程进行求解,得到预测物体温度随时间变化的解析解,以及定义系统响应速度的关键参数——时间常数。

$$ T(t) = T_{\infty}+ (T_i - T_{\infty}) \cdot \exp\left(-\frac{t}{\tau}\right), \quad \tau = \frac{\rho c_p L_c}{h}$$
T_i 初始温度 °C 或 K τ 时间常数 s L_c 特征长度 V/A_s t=τ 剩余温差 约36.8%

现实世界中的应用

电子元件散热设计:工程师用集总热容法快速估算芯片或散热片的温升过程。比如评估一个功率MOSFET在启动后多久会达到稳定温度,这有助于防止过热损坏,是设计冷却风扇和散热器的基础。

热处理工艺:在金属淬火过程中,将炽热的工件浸入冷却液(油或水)。通过估算时间常数τ,可以预测工件心部温度下降到马氏体转变温度所需的时间,从而控制材料的硬度和金相组织。

食品工程与储存:估算罐头食品在杀菌锅中的加热时间,或者预测一块黄油从冰箱取出后,在室温下软化到可涂抹状态需要多久。这关系到食品安全和产品体验。

建筑围护结构响应:虽然墙壁通常不符合Bi<0.1的条件,但对于某些轻薄的内饰板或玻璃,在估算其对室内温度骤变的响应时间时,集总热容法能提供一个非常快速的初步判断。

常见误解与注意事项

要熟练掌握此方法,了解几个“陷阱”至关重要。首先,“Bi<0.1即绝对安全”的误解。这虽是教科书标准,但实际设计中需考虑安全裕度。例如,对于发热元件或材料特性对温度变化敏感(热应力成问题)的部件,即使Bi数在0.05左右,有时也需关注内部温差。即使模拟器显示“适用”,其含义是“可作为第一近似使用”,专业做法是通过其他详细分析或实验来辅助最终判断。

其次,特征长度 $L_c$ 的选取方法。这是易错点。虽然记住“体积÷表面积”,但需根据形状注意。例如冷却细长棒时,若将整体作为单一模块计算 $L_c$,会预测出比实际更慢的冷却速度。此时可采用“分段”技术,将棒体分割为短段并对各段应用集中容量法。在模拟器中尝试改变形状观察时间常数如何变化,有助于掌握其规律。

最后,对流系数 $h$ 的“魔术式调整”。当计算结果与实测不符时,常会想调整 $h$ 值来拟合。但 $h$ 是由流体状态(自然对流或强制对流、流速大小)物理决定的参数。例如静止空气中的自然对流通常在 $5\sim25$ [W/(m²·K)] 量级。若为“拟合”而将其改为100等数值,恰恰说明模型本身已脱离实际。此时应首先怀疑集中容量法的前提(Bi<0.1)是否已不成立。