4. 精密温控对象建模

精密温控系统的控制品质，首先取决于对被控对象动态特性的理解程度。冷水机组的温控对象通常由换热器、循环管路、水箱、负载设备以及温度传感器构成，其本质是一个分布参数的热系统。工程上为了控制器设计，需要将其简化为集总参数模型，核心是热容与热阻两个物理量。

热容与热阻：对象建模的两个基本参数

热容 C 表征对象储存热量的能力，单位为 J/K，等于质量与比热容的乘积。对于循环水回路，若水箱容积 20 L、水比热 4.18 kJ/(kg·K)，则水箱热容约为 83.6 kJ/K。热阻 R 表征热量传递的阻碍，单位为 K/W，由导热、对流和辐射共同决定。换热器水侧与制冷剂侧之间的总热阻可写为：

R_total = 1/(h_w·A_w) + δ/(λ·A) + 1/(h_r·A_r)

其中 h 为对流换热系数，A 为换热面积，δ 为壁厚，λ 为导热系数。热容与热阻的乘积 τ = R·C 即为一阶对象的时间常数，直接决定温度响应速度。示例：若某换热器总热阻 0.02 K/W，水箱热容 83.6 kJ/K，则 τ ≈ 1672 s，约 28 min，说明该对象响应缓慢，控制器积分作用需谨慎整定。

负载扰动类型

精密温控对象的扰动可分为三类。第一类是阶跃型扰动，如负载设备突然开启，发热功率从 0 跳变到 2 kW，温度会按指数规律上升。第二类是斜坡型扰动，如激光器功率缓慢爬升，发热量线性增加，稳态会出现恒定的温度偏差。第三类是周期性扰动，如压缩机启停引起的制冷量脉动，频率通常在 0.01～0.1 Hz。不同扰动类型对控制器的要求不同：阶跃扰动考验响应速度，斜坡扰动考验积分能力，周期扰动则要求足够的相位裕度。

传递函数与状态空间建模

对于单容对象，传递函数可写为一阶惯性加纯滞后形式：G(s) = K·e^-Ls/(τs+1)。其中 K 为静态增益，单位 K/kW；L 为纯滞后，主要来自流体传输和传感器响应。若对象包含换热器金属壁与水体两个明显不同的热容，则需用两状态空间描述：

状态变量：x1 = 水温，x2 = 换热器壁温
ẋ1 = [Q_load - (x1 - x2)/R1] / C_w
ẋ2 = [(x1 - x2)/R1 - (x2 - T_amb)/R2] / C_m
输出方程：y = x1

其中 C_w 为水热容，C_m 为金属壁热容，R1 为水与壁之间热阻，R2 为壁与环境之间热阻。状态空间模型能更准确地描述温度响应的双时间尺度特征，适合模型预测控制等高级算法。

实验辨识方法

工程上最常用的辨识方法是阶跃响应法。具体操作：系统稳定在 20°C 后，将制冷量从 3 kW 阶跃增加到 4 kW，记录水温从 20°C 下降到新稳态的全过程。采样周期建议取 τ/10 左右，示例中 τ≈1672 s，则采样周期 150 s 即可。利用两点法求参数：温度变化达到稳态值 28.3% 和 63.2% 对应的时间分别为 t1 和 t2，则 τ = 1.5(t2 - t1)，L = t2 - τ。若实测 t1=420 s、t2=1150 s，则 τ=1095 s，L=55 s。增益 K = ΔT/ΔQ = (20-17.6)/1 = 2.4 K/kW。

模型验证与误差分析

辨识得到的模型必须经过验证。常用方法是将模型预测输出与实际响应对比，计算拟合度：

fit = (1 - ||y_meas - y_sim|| / ||y_meas - mean(y_meas)||) × 100%

拟合度高于 85% 通常认为模型可用于控制器设计。若拟合度不足，需检查：纯滞后是否估计准确、采样周期是否过大、是否存在未建模的动态（如压缩机启停非线性）。残差分析也很重要：若残差呈现周期性，说明存在未建模的周期扰动；若残差有趋势，说明模型增益或时间常数有偏差。示例中若拟合度仅 72%，可尝试将模型改为两状态空间，通常能提升到 90% 以上。

总结而言，精密温控对象建模的完整链条是：从物理结构提取热容与热阻，判断负载扰动类型，建立传递函数或状态空间模型，通过阶跃响应实验辨识参数，最后用拟合度和残差分析验证模型。只有模型精度足够，后续的 PID 整定、前馈补偿和模型预测控制才有可靠基础。