📚 课程目录第2章温度场建模基础

2. 温度场建模基础

多区加热系统的控制难点，往往不在单个加热器本身，而在于“一个区的热量如何影响另一个区”。要理解这种影响，必须先建立温度场的物理描述：热量通过哪些路径传递、在什么条件下可以简化、在什么条件下必须保留空间分布。本节围绕热传导、对流、辐射三种基本传热方式，说明它们在多区系统中的耦合表达，并给出集总参数模型与分布参数模型的适用边界。

一、三种传热方式在多区系统中的物理角色

热传导发生在固体内部或相互接触的固体之间，遵循傅里叶定律：热流密度与温度梯度成正比，方向由高温指向低温。对多区加热平台而言，相邻加热区之间通过基板、隔热层或金属结构直接导热，这是“邻区串扰”的主要来源。若两区之间温差为 ΔT、等效热阻为 R_cond，则传导热流可近似写为 Q_cond = ΔT / R_cond。热阻由材料导热系数、截面积和路径长度决定，因此加宽隔热槽、减小接触面积或选用低导热材料，都会直接削弱传导耦合。

对流发生在固体表面与流体之间，常用牛顿冷却定律描述：Q_conv = h·A·(T_s − T_f)，其中 h 为对流换热系数，A 为换热面积，T_s 与 T_f 分别为表面温度和流体温度。在多区系统中，如果各加热区共用同一风道或同一腔体，热区会先加热空气，热空气再流经冷区表面，形成“热区→空气→冷区”的间接耦合。这种耦合具有方向性和滞后性，是单纯传导模型难以描述的。

辐射通过电磁波传递能量，遵循斯特藩—玻尔兹曼定律：Q_rad = ε·σ·A·(T₁⁴ − T₂⁴)。辐射与绝对温度的四次方相关，因此在高温区（例如 300 ℃ 以上）不可忽略；在低温区（例如 80 ℃ 以下）通常远小于传导和对流。多区系统中，若各区表面相互可见，辐射会形成“面—面”耦合，且随温度升高迅速增强。

二、耦合表达：从单区方程到多区状态空间

把每个加热区视为一个控制体，写出能量平衡：输入电功率 P_i，减去向环境的对流与辐射散热，再减去向相邻区的传导热流，剩余部分用于自身温升。对第 i 区，集总形式可写为：

C_i · dT_i/dt = P_i − (T_i − T_amb)/R_i,amb
                − Σ_j (T_i − T_j)/R_ij
                − ε_i·σ·A_i·(T_i^4 − T_surr^4)

其中 C_i 为第 i 区热容，R_i,amb 为第 i 区到环境的等效热阻，R_ij 为第 i 区与第 j 区之间的等效热阻。把所有区写在一起，就得到多区耦合方程组。若忽略辐射非线性，可线性化为状态空间形式：

dT/dt = A·T + B·u + E·T_amb
y     = C·T

矩阵 A 的对角元表示各区自散热与自加热能力，非对角元表示区与区之间的热耦合强度。非对角元越大，说明一个区的功率变化越容易“拉动”另一个区的温度，解耦控制就越必要。实际工程中，A 的非对角元可通过阶跃实验辨识：给第 j 区一个已知功率阶跃，记录第 i 区的温度响应，用响应幅值和上升时间估计 R_ij 与 C_i。

三、集总参数与分布参数模型的适用边界

集总参数模型假设每个区内部温度均匀，只保留一个温度节点。它的成立条件是毕渥数 Bi = h·L_c/k 远小于 0.1，即内部导热热阻远小于表面换热热阻。对于薄金属加热板、小型加热块，Bi 通常很小，集总模型足够准确。此时模型阶数等于区数，参数少、辨识快、控制器设计简单，适合实时控制。

当区尺寸较大、材料导热系数较低，或加热功率分布不均匀时，区内会出现明显温度梯度，必须使用分布参数模型。分布参数模型用偏微分方程描述：ρ·c_p·∂T/∂t = ∇·(k∇T) + q_gen，边界条件包含对流和辐射。求解方法包括有限差分、有限元或模态展开。分布模型能预测区内热点、边缘效应和热流路径，但计算量大，通常用于离线分析、传感器布置优化或模型降阶的基准。

工程判断可参考以下边界：若区特征尺寸 L 与热扩散长度 √(α·t) 接近或更大，且区内温差超过目标均匀性的 1/3，则应考虑分布模型；若区尺寸小、材料导热好、目标均匀性宽松，集总模型即可。实际做法常是“分布模型离线标定、集总模型在线控制”：用有限元或实验辨识出等效热阻和热容，再嵌入集总状态空间进行解耦控制。

教学示例（合理假设）：某三区铝制加热平台，单区尺寸 100 mm × 100 mm × 10 mm，铝导热系数 k ≈ 237 W/(m·K)，表面对流系数 h ≈ 15 W/(m²·K)。特征长度 L_c ≈ 0.01 m，Bi = 15×0.01/237 ≈ 0.00063 ≪ 0.1，因此单区内部可用集总模型。若换成 10 mm 厚的陶瓷板（k ≈ 2 W/(m·K)），Bi ≈ 0.075，仍接近集总边界；若厚度增至 50 mm，Bi ≈ 0.375，则必须考虑分布效应。

四、耦合强度与解耦需求的量化判断

多区系统是否需要解耦，不能只看“有没有耦合”，而要看耦合强度相对于自调节能力的比例。定义耦合比 γ_ij = |A_ij| / |A_ii|。当 γ_ij 小于约 0.1 时，邻区影响通常可被单区闭环抑制；当 γ_ij 在 0.1～0.3 之间，建议加入前馈或静态解耦；当 γ_ij 大于 0.3 时，强耦合会导致单区 PID 相互“打架”，应优先采用多变量解耦控制或模型预测控制。

耦合比可通过实验估计：稳定运行时，给第 j 区叠加一个幅值为 ΔP 的功率阶跃，记录第 i 区温度变化 ΔT_i 和第 j 区温度变化 ΔT_j，则 γ_ij ≈ (ΔT_i/ΔT_j)·(C_j/C_i)。若两区热容相近，γ_ij 可简化为温升比。这个比值直接告诉工程师：调一个区时，另一个区会被“带跑”多少。

常见误区：把辐射项简单线性化为常数热阻。辐射热流与 T⁴ 相关，在高温区非线性显著；若在 300 ℃ 以上仍用固定 R_rad，模型会在升温段严重低估散热，导致控制器超调。正确做法是在工作点附近线性化，或保留辐射项进行非线性补偿。

五、建模流程与参数获取建议

实际建模可按以下顺序推进：第一，明确各区几何、材料、加热功率和传感器位置；第二，用 Bi 数判断单区是否可集总；第三，通过稳态实验估计各区到环境的等效热阻 R_i,amb；第四，通过邻区阶跃实验估计 R_ij；第五，用降温曲线估计热容 C_i；第六，组装状态空间矩阵并验证预测误差。若预测误差在目标温区内超过 ±2 ℃，应检查辐射项、传感器滞后或分布效应是否被忽略。

参数获取时，建议在多个工作点重复实验：低温段（如 60～100 ℃）以对流和传导为主，高温段（如 200～350 ℃）辐射占比上升。若同一组 R_ij 在低温段拟合良好、高温段偏差大，说明模型缺少辐射非线性项，需要补充或改用分布参数模型。最终模型不必追求物理上完全精确，而应满足控制目标：在目标温区内，模型预测的邻区耦合增益误差小于 20%，即可用于解耦控制器设计。