2. 食品加工设备温度对象特性与建模
要对夹层锅的温度实施可靠的PID控制,第一步不是调参数,而是弄清楚“被控对象到底长什么样”。温度对象特性决定了PID能调多快、能调多稳,也决定了参数整定的基本方向。本节围绕热交换原理、设备热惯性、滞后与非线性三个特性展开,并给出基于阶跃响应辨识一阶加纯滞后(FOPDT)模型的方法,最后以夹层锅为例完成一次建模实战。
一、热交换原理:热量如何进入物料
夹层锅属于间壁式换热设备:蒸汽(或导热油)在锅体夹套内冷凝放热,热量穿过内壁传给锅内物料。以蒸汽夹层锅为例,传热过程可分解为三个串联环节:蒸汽侧冷凝对流换热、不锈钢壁导热、物料侧自然对流或强制对流换热。稳态下传热速率可写作:
Q = K · A · ΔT
K = 1 / (1/α₁ + δ/λ + 1/α₂)
其中 Q 为传热速率(W),A 为有效换热面积(m²),ΔT 为蒸汽温度与物料温度之差,K 为总传热系数,α₁、α₂ 分别为蒸汽侧与物料侧对流换热系数,δ 为壁厚,λ 为壁材导热系数。这个式子说明一个关键事实:物料温度越低,ΔT 越大,升温越快;温度越接近蒸汽温度,升温越慢。这正是温度对象非线性的根源之一。同时,物料侧 α₂ 会随物料黏度、搅拌状态变化,例如熬制高黏度酱料时 α₂ 明显下降,K 随之减小,对象增益发生漂移。
二、热惯性:为什么温度“拉不住”
夹层锅由锅体金属、夹套内蒸汽与冷凝水、锅内物料共同构成一个具有热容的集合体。热量输入后,一部分用于提升物料温度,一部分被锅体和夹套吸收,还有一部分通过锅壁、锅口向环境散失。其能量平衡可近似写为:
C · dT/dt = Q_in − Q_loss
C = m_m·c_m + m_w·c_w (锅体与物料热容之和)
热容 C 越大,同样的热量输入引起的温升越小、越慢,这就是热惯性。工程上表现为:加大蒸汽阀门后温度不会立刻上升,关小阀门后温度还会继续“冲”一段。热惯性直接对应模型中的时间常数 T,T 越大,系统响应越迟缓,PID 的积分作用需要更耐心,微分作用也更容易被噪声干扰。
三、滞后与非线性:控制难度的来源
滞后分为两类。一是纯滞后:蒸汽阀动作后,热量要先改变夹套压力与冷凝状态,再穿过锅壁,才影响到测温点,测温元件本身还有响应时间,这些环节叠加形成一段“死时间” τ。二是容量滞后:物料内部温度分布不均,靠近锅壁处先升温,中心处后升温,测温点位置不同,响应曲线也不同。纯滞后对PID最不友好——控制器在死时间内看不到任何变化,会持续累积输出,等温度终于开始变化时往往已经超调。
非线性则主要体现在三处:升温段与保温段的传热温差不同,导致对象增益随工作点变化;蒸汽阀门小开度时流量特性非线性;物料黏度、装料量变化引起热容与传热系数变化。因此,用一组固定PID参数覆盖整个升温—保温过程,往往在保温段振荡、在升温段偏慢,这也是后续章节要讨论分段控制或增益调度的原因。
四、FOPDT 模型与阶跃响应辨识
食品加工中的温度对象,绝大多数可以用一阶加纯滞后模型近似描述:
G(s) = K · e^(−τs) / (T·s + 1)
三个参数各有明确物理含义:增益 K 表示单位阀门开度变化引起的稳态温升(℃/%);时间常数 T 表示温度完成总变化量 63.2% 所需的时间(s 或 min);纯滞后 τ 表示从阀门动作到温度开始明显变化的时间。辨识方法就是做一次开环阶跃试验,从响应曲线上读出这三个量。
具体做法:设备处于稳定工况(例如物料 60℃保温),手动把蒸汽阀从 30% 阶跃到 40% 并保持不变,同时以固定采样周期(如 5 s)记录温度,直到温度基本稳定。然后按以下步骤处理:
- 计算稳态温升 ΔT_ss = T_终 − T_始,得到 K = ΔT_ss / Δu。示例:Δu = 10%,温升 18℃,则 K = 1.8 ℃/%。
- 在响应曲线拐点处作切线,切线与初始温度水平线的交点对应时刻即 τ;切线与稳态温度水平线的交点对应时刻记为 t₂,则 T = t₂ − τ。
- 也可用“63.2% 法”校核:找到温度达到 ΔT_ss 的 63.2% 的时刻 t₆₃,则 T = t₆₃ − τ。两种方法结果接近即可采用。
五、夹层锅建模实战
假设一台 300 L 蒸汽夹层锅熬制番茄酱,正常装料 200 kg,搅拌 30 r/min。在物料 65℃稳定时,将蒸汽阀由 35% 阶跃至 45%,采样记录得到如下示例数据(合理假设,非实测):
| 时刻 t/s | 温度/℃ | 说明 |
|---|---|---|
| 0 | 65.0 | 阶跃发生 |
| 60 | 65.1 | 死时间内基本无变化 |
| 120 | 66.8 | 开始明显上升 |
| 240 | 72.5 | 快速上升段 |
| 420 | 78.9 | 接近 63.2% 点 |
| 720 | 82.6 | 趋于平缓 |
| 1200 | 83.0 | 基本稳定 |
由表可得 ΔT_ss = 83.0 − 65.0 = 18.0℃,Δu = 10%,故 K = 1.8 ℃/%。从曲线看,前约 90 s 温度几乎不动,取 τ ≈ 90 s。总变化量的 63.2% 为 11.4℃,对应温度 76.4℃,从表中插值约在 t ≈ 380 s,故 T = 380 − 90 = 290 s。于是得到示例模型:
G(s) = 1.8 · e^(−90s) / (290s + 1)
这组参数给出了明确的整定预期:纯滞后 90 s 意味着控制器输出变化后要等一分半才能看到温度响应,积分时间不宜设得过短;时间常数 290 s 说明对象较迟缓,比例带可适当放宽。若后续更换为高黏度物料,可重新做一次阶跃试验,通常会看到 K 下降、T 增大,说明对象变得更“钝”,原PID参数需要相应调整。
需要强调的是,FOPDT 只是工程近似。夹层锅在升温初期与保温后期增益并不相同,严格来说是非线性对象;但在工作点附近小偏差范围内,用 FOPDT 描述已足以支撑PID整定。建模的价值不在于模型多精确,而在于把“感觉温度很难控”转化为 K、T、τ 三个可比较、可复用的量化依据,为后续参数整定和异常诊断提供基准。