连续搅拌罐加热器的输入-输出模型可通过应用非稳态能量衡算获得，将入口温度和加热负荷定义为输入，将罐内温度定义为输出。 对于体积恒定且物性恒定的储罐，其模型为
[ \frac{dT}{dt}
\frac{F}{V}(T_i-T) + \frac{Q}{\rho c_p V}, ]
其中，(T) 为充分混合的罐内温度，(T_i) 为入口温度，(F) 为体积流量，(V) 为液体体积，(Q) 为净传热速率。在稳态附近进行线性化，可得到适用于传递函数分析和反馈控制器设计的一阶输入-输出模型。
核心建模方法是使用偏差变量表示动态能量衡算。以 (T_i) 和 (Q) 为输入、以 (T) 为输出时，恒容CSTH近似表现为一个一阶过程，同时包含入口温度和加热负荷的影响。
定义控制问题
选择被控输出
通常的被控变量是罐内液体温度 (T(t))，因为它决定了离开加热器的工艺物流的热状态。
假设储罐内液体完全混合，因此出口温度等于罐内温度：
[ T_{\text{out}}=T. ]
确定过程输入
主要操纵输入是加热负荷 (Q(t))，可通过蒸汽、电加热器或加热夹套提供。
入口温度 (T_i(t)) 通常被视为扰动输入，但在某些中试装置配置中也可以作为操纵变量。流量 (F(t)) 也可能作为扰动或操纵变量，因为它同时改变停留时间和能量输送量。
建立建模假设
一个有用的一阶模型通常作出以下假设：
- 液体完全混合
- 液体体积恒定
- 密度 (\rho) 恒定
- 热容 (c_p) 恒定
- 入口和出口流量恒定，(F_i=F)
- 忽略热损失，除非对其进行显式建模
- 忽略动能、势能、轴功和膨胀功的影响
这些假设在保留主要热动态特性的同时，形成了紧凑的模型。
建立动态衡算
编写总质量衡算
一般液体体积衡算为
[ \frac{dV}{dt}=F_i-F. ]
对于恒容CSTH，
[ F_i=F, \qquad \frac{dV}{dt}=0. ]
如果液位动态很重要，则必须将储罐体积保留为状态变量，而不能将其视为常数。
编写能量衡算
对于非反应型搅拌罐加热器，非稳态能量衡算为
[ \rho c_p V\frac{dT}{dt}
\rho c_p F(T_i-T)+Q. ]
除以 (\rho c_p V) 后，得到温度状态方程：
[ \frac{dT}{dt}
\frac{F}{V}(T_i-T) + \frac{Q}{\rho c_p V}. ]
等价地，
[ \frac{dT}{dt} + \frac{F}{V}T
\frac{F}{V}T_i + \frac{1}{\rho c_p V}Q. ]
右侧第一项表示随进料进入的对流能量，而 (Q) 表示加入液体的净热量。冷却负荷可以通过取 (Q<0) 来表示。
在需要时加入反应效应
如果CSTH同时作为反应型CSTR运行，能量衡算中可以加入反应热：
[ \frac{dT}{dt}
\frac{F_i}{V}(T_i-T) + \frac{(-\Delta H_r)r_A}{\rho c_p} + \frac{Q}{\rho c_p V}, ]
其中，(r_A) 为反应速率。
在这种情况下，温度通过反应动力学与浓度耦合。对于一级反应，
[ r_A=k_0\exp\left(-\frac{E}{RT}\right)c_A, ]
完整的非线性模型还必须包括 (c_A) 的组分衡算。
将模型转换为偏差变量
定义稳态工作点
设稳态值为 (T_s)、(T_{i,s}) 和 (Q_s)。在稳态下，
[ \frac{dT}{dt}=0, ]
因此工作点满足
[ 0= \frac{F}{V}(T_{i,s}-T_s) + \frac{Q_s}{\rho c_p V}. ]
该工作点是分析动态变化时所使用的参考点。
定义偏差
定义
[ T'=T-T_s, ]
[ T_i'=T_i-T_{i,s}, ]
[ Q'=Q-Q_s. ]
对于恒定的 (F)、(V)、(\rho) 和 (c_p)，从动态方程中减去稳态方程可得
[ \frac{dT'}{dt} + \frac{F}{V}T'
\frac{F}{V}T_i' + \frac{1}{\rho c_p V}Q'. ]
使用停留时间定义
[ \tau=\frac{V}{F}, ]
模型变为
[ \boxed{ \frac{dT'}{dt} + \frac{1}{\tau}T'
\frac{1}{\tau}T_i' + \frac{1}{\rho c_p V}Q' } ]
这是恒容连续搅拌罐加热器的标准线性一阶输入-输出模型。
将模型与通用参数关联
模型也可以写成通用形式
[ \frac{dT'}{dt}+aT'
\frac{1}{\tau}T_i' + K_QQ', ]
其中
[ a=\frac{F}{V}=\frac{1}{\tau}, ]
以及
[ K_Q=\frac{1}{\rho c_p V}. ]
如果 (Q) 表示蒸汽流量信号、阀门开度或加热器指令，而不是实际热量速率，则 (K_Q) 必须包含相应的执行器增益和传热增益。
推导传递函数
应用拉普拉斯变换
对方程进行拉普拉斯变换，并假设初始温度偏差为零，可得
[ sT'(s)+\frac{1}{\tau}T'(s)
\frac{1}{\tau}T_i'(s) + K_QQ'(s). ]
因此，
[ T'(s)
\frac{1/\tau}{s+1/\tau}T_i'(s) + \frac{K_Q}{s+1/\tau}Q'(s). ]
获得入口温度传递函数
入口温度扰动到罐内温度的传递函数为
[ \boxed{ \frac{T'(s)}{T_i'(s)}
\frac{1}{\tau s+1} } ]
因此，入口温度通过时间常数为 (\tau=V/F) 的一阶动态响应影响罐内温度。
获得加热负荷传递函数
加热负荷到罐内温度的传递函数为
[ \boxed{ \frac{T'(s)}{Q'(s)}
\frac{K_Q}{s+1/\tau}
\frac{K_T}{\tau s+1} } ]
其中
[ K_T=\frac{\tau}{\rho c_p V} =\frac{1}{\rho c_p F}. ]
因此，对于直接的热量速率输入，
[ \boxed{ \frac{T'(s)}{Q'(s)}
\frac{1}{\rho c_p F}\frac{1}{\tau s+1} } ]
温度总响应是扰动输入和操纵输入所产生响应之和：
[ T'(s)
\frac{1}{\tau s+1}T_i'(s) + \frac{K_T}{\tau s+1}Q'(s). ]
扩展模型以描述液位动态
将液位作为状态变量加入
如果储罐液位不是恒定的，则定义 (h(t)) 为液体高度，(A) 为储罐横截面积。质量衡算变为
[ A\frac{dh}{dt}=F_i-F. ]
能量衡算可以写为
[ Ah\frac{dT}{dt}
F_i(T_i-T) + \frac{Q}{\rho c_p}. ]
其中因子 (Ah) 是液体体积。因此，液位变化会改变热容，从而改变温度动态特性。
判断是否需要液位耦合
当液位控制使 (V) 基本保持恒定且液位变化较小时，适合使用仅包含温度的模型。
当流量变化导致体积发生显著变化、液位控制属于控制目标的一部分，或运行过程中停留时间发生较大变化时，则需要使用液位-温度耦合模型。
理解模型的控制含义
解释时间常数
时间常数
[ \tau=\frac{V}{F} ]
是液体的名义停留时间。储罐体积越大，温度变化越慢；流量越大，罐内物料更换越快。
发生阶跃变化后，一阶过程在一个时间常数后约达到最终温度变化的63%，并渐近地接近最终值。
解释过程增益
对于加热负荷变化，稳态增益为
[ K_T=\frac{1}{\rho c_p F}. ]
给定热输入所引起的温升会随着流量、密度或热容的增加而降低。
对于入口温度变化，在理想恒物性模型下，稳态增益为1：
[ \lim_{s\to 0}\frac{T'(s)}{T_i'(s)}=1. ]
在稳态下，罐内温度最终会跟随持续的入口温度变化，但还会受到加热负荷控制器作用的影响。
使用模型进行反馈设计
该模型确定了温度控制回路所需的关键要素：
- 被控变量：罐内温度 (T)
- 操纵变量：加热负荷 (Q)
- 主要扰动：入口温度 (T_i)
- 过程时间常数：(\tau=V/F)
- 过程增益：(K_T=1/(\rho c_p F))
- 测量装置：温度变送器
- 最终控制元件：蒸汽阀、加热器或夹套负荷执行器
这些参数支持PID整定、扰动抑制分析、仿真，以及中试装置预测响应与实测响应之间的比较。
理解各种权衡
未经验证，不要将非线性反应器视为一阶过程
对于非反应型或弱耦合热系统，在固定工作点附近，一阶模型是有效的。
对于放热反应器，阿伦尼乌斯反应速率会引入强烈的温度依赖性。由此产生的模型可能具有多个稳态、非线性增益或热不稳定性，而简单的一阶传递函数无法表示这些特性。
热损失显著时不要忽略
如果热量通过储罐壁、保温层或气相空间损失，则可能需要加入额外项：
[ Q_{\text{loss}}=UA(T-T_{\text{amb}}). ]
此时能量衡算中包含的净热输入为
[ Q_{\text{net}}=Q-UA(T-T_{\text{amb}}). ]
忽略显著的热损失会导致稳态增益和控制器整定参数出现误差。
区分热量速率与加热器指令
如果直接使用增益 (1/(\rho c_p V))，变量 (Q) 应表示实际传热速率。
如果测量输入是阀门开度、蒸汽流量、电功率指令或夹套温度，则必须单独建立执行器与传热关系模型，或将其纳入通过实验辨识得到的过程增益中。
考虑测量和执行器动态
实际中试装置通常包括传感器滞后、阀门动态、传输延迟和蒸汽冷凝动态。
如果这些影响很显著，过程模型可以采用以下形式：
[ G(s)= \frac{K_T e^{-\theta s}} {(\tau s+1)(\tau_m s+1)}, ]
其中，(\theta) 为延迟，(\tau_m) 为测量或执行器时间常数。
根据目标做出正确选择
应选择与控制问题和运行条件相匹配的模型复杂度。
- 如果主要关注基础温度控制： 使用恒容一阶模型，将 (T_i) 作为扰动，将 (Q) 作为操纵输入。
- 如果主要关注液位和温度控制： 保留质量衡算和能量衡算，并将 (h) 或 (V) 以及 (T) 作为状态变量。
- 如果主要关注反应器行为： 加入组分衡算和与温度相关的反应动力学，然后仅围绕明确定义的稳态进行线性化。
- 如果主要关注中试装置控制器整定： 当实测执行器、传感器、热损失和延迟效应足以影响观察到的响应时，应将其纳入模型。
- 如果主要关注仿真或故障检测： 在将模型用于操作员培训、优化或诊断之前，应使用阶跃测试数据对其进行验证。
建立良好的CSTH输入-输出模型，可以将守恒定律与可测量的过程行为联系起来，为控制设计和中试装置的安全运行提供可靠依据。
汇总表：
| 模型要素 | 定义 | 控制意义 |
|---|---|---|
| 被控输出 | 罐内温度，T | 决定出口热状态 |
| 操纵输入 | 加热负荷，Q | 调节过程温度 |
| 主要扰动 | 入口温度，Ti | 改变热负荷 |
| 动态衡算 | dT/dt = F/V(Ti - T) + Q/(rho cp V) | 描述温度行为 |
| 时间常数 | tau = V/F | 表示液体停留时间 |
| 加热负荷增益 | KT = 1/(rho cp F) | 关联热输入与温度变化 |
| 传递函数 | T'(s)/Q'(s) = KT/(tau s + 1) | 支持PID整定和反馈设计 |
| 模型扩展 | 液位、反应、热损失、延迟、执行器和传感器动态 | 在需要时提高中试装置的建模准确性 |
通过LABPARK规划更可靠的温度控制实验。我们的教育和职业培训单元操作中试装置帮助大学、研究机构和企业开展化学工程过程的教学、测试与验证。针对化学工程、生物过程与生物技术，或环境与水处理应用，选择具有真实过程动态特性的实用系统，用于能量衡算、传递函数分析、PID整定和操作员培训。联系LABPARK，讨论您所需的规模、配置、仪表和控制目标。
