1. 反应釜温度单回路PID参数整定与阶跃响应验证

对象建模：从夹套加热实验到 FOPDT 近似

夹套加热反应釜的温度对象具有明显的热惯性与传输滞后。工程整定前，先用一阶加纯滞后（FOPDT）模型近似：

G(s) = K·e^(-τs) / (T·s + 1)

以某 5 m³ 夹套反应釜的教学示例为例：夹套蒸汽阀开度阶跃 +10%，釜内温度由 80.0 ℃ 稳态升至 88.0 ℃ 新稳态，纯滞后约 45 s，时间常数约 300 s。由此得到：

• 增益 K = ΔT / Δu = 8.0 ℃ / 10% = 0.8 ℃/%
• 滞后 τ = 45 s，时间常数 T = 300 s，τ/T = 0.15（属可控性较好的对象）

该组参数为教学示例假设值，实际装置应通过阶跃测试辨识，并注意测试期间进料、搅拌与室温的稳定。

PID 参数整定：IMC 与 Z-N 对照

对 FOPDT 对象，采用内模控制（IMC）整定公式，取闭环时间常数 λ = max(τ, 0.1T)：

Kc = T / (K·(λ + τ))
Ti = T
Td = 0（PI）或 τ/2（PID）

代入 K=0.8、T=300、τ=45、λ=45：

工程上先以 IMC-PID 作为初值，再按超调量微调 Kc 与 Td。注意：温度回路执行器为蒸汽调节阀，Kc 过大易引起阀门频繁动作，需配合输出限幅与速率限制。

阶跃仿真验证

用离散 PID 对 FOPDT 模型做设定值阶跃仿真（设定值由 80 ℃ 阶跃至 85 ℃），采样周期 Ts = 5 s，仿真时长 1800 s。核心代码如下（Python，可直接运行）：

import numpy as np

K, T, tau = 0.8, 300.0, 45.0
Kc, Ti, Td = 4.17, 300.0, 22.5
Ts, Tsim = 5.0, 1800.0
sp = 85.0
N = int(Tsim / Ts)

# 对象离散化（一阶惯性 + 纯滞后用队列实现）
a = np.exp(-Ts / T)
y = 80.0
u = 0.0
buf = [0.0] * int(tau / Ts)
integ = 0.0
prev_e = 0.0
ys, us = [], []

for k in range(N):
    e = sp - y
    integ += e * Ts
    deriv = (e - prev_e) / Ts
    u = Kc * (e + integ / Ti + Td * deriv)
    u = min(max(u, 0.0), 100.0)          # 阀门开度限幅
    prev_e = e
    buf.append(u)
    u_delay = buf.pop(0)
    y = a * y + (1 - a) * (80.0 + K * u_delay)
    ys.append(y); us.append(u)

ys = np.array(ys)
overshoot = (ys.max() - sp) / (sp - 80.0) * 100
# 调节时间：进入 ±2% 带且不再越出
band = 0.02 * (sp - 80.0)
idx = np.where(np.abs(ys - sp) > band)[0]
ts = (idx[-1] + 1) * Ts if len(idx) else 0
print(f"超调量 = {overshoot:.1f}%")
print(f"调节时间 = {ts:.0f} s")

示例运行预期（由给定模型与参数直接推出，非实测）：超调量约 3%~6%，调节时间约 600~750 s。若工艺指标为“超调 ≤ 10%、调节时间 ≤ 900 s”，则该组参数合格。

整定与验证关系图

现场落地要点

• 仿真参数下装前，确认控制器采样周期与仿真一致（本例 5 s），否则积分与微分作用会偏离预期。
• 蒸汽阀存在死区与摩擦，实际 Kc 可比仿真值略小 10%~20%，先保守再逐步放开。
• 阶跃幅值不宜过大，建议 3~5 ℃，避免触发工艺联锁或物料分解。
• 记录每次阶跃的设定值、实际温度曲线、阀门开度曲线，作为后续参数微调依据。