集中容量法模拟器 返回
热分析

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

确认生物数(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
实时冷却/加热动画
当前温度 T°C
经过时间 t0.0s
时间常数 τs
毕渥数 Bi
距平衡%
集总有效?

整个物体颜色均匀=温度均匀(集总假设)。当 Bi>0.1 时内部出现温度梯度,假设失效。

计算结果
生物数 Bi
有效 (Bi < 0.1)
时间常数 τ
s
T at t=τ
°C
99%平衡到达时间
s
几何形状
sphere
形状
有效 Lc
--
mm
不适用集中容量法(Bi > 0.1)— 需进行一维分析

$$T(t) = T_\infty + (T_0 - 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,方法就有效。这个工具会自动计算Bi值。试试增大"对流系数 h",当Bi超过0.1时,会显示"集中容量法需注意"的警告。这就是检查适用条件。
🙋
明白了!那图中的"时间常数τ"是什么?好像跟"99%到达时间"有关系?
🎓
完全正确!时间常数τ是决定"温度变化速率"的最关键参数。τ越小,物体冷却(或加热)越快。例如,比热小、热易传导的铝材τ较小,迅速接近流体温度。相反,绝热材料的τ很大。看右边的图,"无量纲温度"在对数刻度上是一条直线下降——这正是集中容量法的特征。

物理模型与主要公式

集中容量法的核心是假设物体内部的热传导阻抗相对表面对流热传达阻抗可以忽视(Bi < 0.1),用一个方程表示物体的能量平衡。

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

解这个微分方程,可得物体温度随时间变化的曲线(冷却/加热曲线)为指数函数形式。

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

常见问题

一般来说,Bi < 0.1 是应用的标准。在此条件下,物体内部温度分布几乎均匀。本模拟器实时显示Bi值,当超过0.1时会发出警告,可直观确认适用范围。
时间常数τ = ρVcp/(hAs),取决于材料密度、比热、体积。钢的时间常数大,冷却慢;铝的时间常数小,冷却快。在模拟器中切换材料,观察温度曲线的斜率变化,就能理解这一点。
自然对流约5~25 W/(m²·K),强制对流(风扇等)约10~500 W/(m²·K),水冷可达500~1000以上。根据实际情况调整,并确认Bi值小于0.1。
Bi < 0.1 条件下,误差在数%以内,与实测吻合良好。但物体形状复杂或Bi较大时误差会增大。本工具是集中容量法的一次近似,精密设计应补充有限元分析等详细方法。

实际应用

电子器件热设计:发热IC芯片或基板的温升计算中使用。集中容量法确定大致时间常数,评估过渡热失控风险,为散热器和冷却风扇选型提供参考。

金属热处理工艺:炉内加热或冷却后的金属部件,根据部件材质、尺寸和炉内对流条件,估算达到目标温度的时间,是制定热处理工艺制度的基础。

食品和制药的温度管理:灭菌加热、冷藏冷冻过程建模。例如罐头中心温度达到灭菌温度的时间,模拟器给出简化估算,平衡杀菌效果与产品品质。

建筑物节能分析:墙体、窗玻璃等建筑构件对外界气温变化的延迟响应,用集中容量法的一次近似评估。对比不同绝热性能的时间常数差异。

常见误区与注意事项

要充分发挥这个方法的威力,需要了解几个"陷阱"。首先,"Bi<0.1就绝对安全"的误解。虽然这是教科书标准,但实际设计中要考虑安全裕度。例如,有内热源的部件或热应力敏感部件,即使Bi=0.05也应关注内部温度差。模拟器显示"适用"只说明"可用第一近似",最终判断要配合其他详细分析或实验。

其次,特征长度 $L_c$ 的选择陷阱。虽然公式是"体积÷表面积",但细长构件需谨慎。例如冷却细棒时,用整体 $L_c$ 计算会低估冷却速率。此时可采用"分段化"策略——将棒分成多段,分别应用集中容量法。在模拟器中尝试改变形状,观察时间常数变化,能帮助建立直觉。

最后,对流系数 $h$ 的"随意调参"。计算结果与实测不符时,容易冲动地改动 $h$ 来"凑数字"。但 $h$ 由流体状态(自然对流还是强制对流、流速多少)物理决定。静止空气自然对流一般5~25 [W/(m²·K)]。若要调到100才能匹配,说明集中容量法前提(Bi<0.1)可能已破裂。此时应先怀疑内部温度梯度,而非盲目改参数。

使用指南

  1. 在输入框中设置材料密度ρ(kg/m³)、比热cp(J/kg·K)、热导率k(W/m·K)
  2. 根据流体侧热传达系数h(W/m²·K)和代表寸法(m),自动计算生物数Bi=hL/k,确认Bi≤0.1的有效性
  3. 输入初始温度与环境温度差,点击"执行计算",温度衰减曲线T(t)=Tₛ+(T₀-Tₛ)exp(-t/τ)实时图形显示

具体计算例

厚度10mm钢板(ρ=7850kg/m³、cp=490J/kg·K、k=50W/m·K)油冷淬火,h=1000W/m²·K条件下,生物数Bi=1000×0.005/50=0.1,集中容量法适用。时间常数τ=ρcpL/h=7850×490×0.005/1000≈19.2秒,300秒后温度与初始温差降至0.31%。铝合金(ρ=2700kg/m³、cp=900J/kg·K、k=160W/m·K)同条件下τ≈7.6秒,大幅缩短。

工程应用注意