搅拌罐加热器使用两个状态变量进行建模:液位(h)和液体温度(T)。总质量衡算决定液位如何变化,而能量衡算决定充分混合液体的温度如何变化。这些方程共同预测入口流量、出口流量、入口温度或蒸汽热输入发生变化后的瞬态行为,并在两个时间导数均为零时定义稳态条件。
液位和温度是系统的动态记忆。质量守恒支配(h),能量守恒支配(T),由此得到的微分方程将中试装置的输入和扰动与可测量的过程输出联系起来。
识别状态变量
液位作为状态变量
对于横截面积恒定为(A)的储罐,液体体积为(V=Ah)。因此,在密度视为恒定时,液位(h)足以描述储罐中储存的液体质量。
质量衡算为
[ A\frac{dh}{dt}=F_i-F ]
其中,(F_i)为入口体积流量,(F)为出口体积流量。
该方程表明,当流入量大于流出量时液位上升,而当流出量大于流入量时液位下降。如果(F_i=F),液位保持不变。
温度作为状态变量
由于储罐处于搅拌状态,通常假设液体完全混合。因此,出口温度等于储罐温度(T)。
当控制目标包含加热性能或流出液温度时,温度就成为必需的状态变量。它表示液体中储存的热能,并在能量进入或离开储罐时发生变化。
根据控制目标选择状态变量
液位控制实验可能只需要质量衡算,并将(h)作为状态变量。温度控制实验则需要质量衡算和能量衡算,因为储罐体积会影响其热容。
这是一项实用的建模原则:应包含描述受控输出及其动态储能机制所需的状态,但应避免加入不会影响控制目标的衡算方程。
应用质量守恒
建立动态衡算
一般守恒关系为
[ \text{累积}=\text{输入}-\text{输出} ]
对于液体密度恒定且体积为(V=Ah)的储罐,
[ \frac{dV}{dt}=F_i-F ]
代入(V=Ah)可得
[ A\frac{dh}{dt}=F_i-F ]
这是关于液位状态的一阶微分方程。
解释输入量
入口流量(F_i)通常是操纵变量或扰动,具体取决于装置配置。出口流量(F)可以由泵或阀直接控制,也可以通过重力驱动排放并随液位变化。
质量衡算不会直接决定温度。它决定可用于储存热能的液体量,因此会间接影响温度动态。
确定液位稳态
在稳态下,
[ \frac{dh}{dt}=0 ]
因此质量衡算变为
[ F_i=F ]
该条件有助于检查数学模型和中试装置的运行点。入口流量和出口流量持续不匹配,表明液位必然会继续变化。
应用能量守恒
建立热量衡算
能量衡算考虑进料带入的热能、出口带出的热能,以及蒸汽加热器输入的热能。在密度和热容恒定的假设下,衡算为
[ Ah\frac{dT}{dt}
F_i(T_i-T)+\frac{Q}{\rho c_p} ]
其中:
- (T_i)为入口液体温度,
- (Q)为蒸汽提供的热量,
- (\rho)为液体密度,
- (c_p)为液体热容。
左侧表示与液体体积相关的热储存项。(F_i(T_i-T))表示用入口液体替换储罐液体所产生的对流热效应净值。
理解热输入项
当热负荷(Q)为正时,它会提高储罐温度。除以(\rho c_p)可将输入热功率转换为等效的体积温度变化贡献。
在更详细的装置模型中,蒸汽条件与(Q)之间的关系可以通过传热系数、盘管面积、蒸汽压力或执行器模型表示。守恒方程仍然是过程层面对热量如何改变液体温度的描述。
认识液位耦合
能量方程在储存项(Ah,dT/dt)中包含(h)。液位越高,液体热容越大,因此相同的热输入通常会产生更慢的温度响应。
这种耦合意味着,即使蒸汽输入不变,液位扰动也可能影响温度动态。因此,当两个变量均发生显著变化时,应联立求解液位和温度方程。
确定温度稳态
在稳态下,
[ \frac{dT}{dt}=0 ]
能量衡算变为
[ 0=F_i(T_i-T)+\frac{Q}{\rho c_p} ]
或者,在上述假设并采用一致单位的条件下,
[ T=T_i+\frac{Q}{\rho c_p F_i} ]
该表达式说明稳态温度取决于入口温度、流量和热负荷。
使用模型进行过程动态分析
定义输入、扰动和输出
典型操纵变量包括蒸汽热流量(Q)、入口流量(F_i)或出口流量(F)。常见扰动包括入口温度(T_i)、进料流量、蒸汽压力或液体性质的变化。
常见测量输出包括储罐液位(h)、储罐温度(T)、出口温度,有时还包括流量。状态方程提供了连接这些量的内部动态模型。
仿真瞬态响应
给定初始条件(h(0))和(T(0))后,可以在输入发生阶跃变化或其他变化后,对两个微分方程进行时间积分。
例如,改变(Q)会产生温度瞬态,而改变(F_i)可能同时产生液位瞬态和温度响应。当方程为非线性或时变时,可以使用 MATLAB 的 ode45 或 Python 的 solve_ivp 等数值求解器计算轨迹。
在运行点附近线性化
在控制分析中,变量通常表示为相对于稳态的偏差:
[ h'=h-h_s,\qquad T'=T-T_s,\qquad T_i'=T_i-T_{i,s},\qquad Q'=Q-Q_s ]
其中,下标(s)表示标称运行点。
在液位和流量恒定等附加假设下,温度模型可以简化为一阶形式:
[ \frac{dT'}{dt}+aT'
\frac{1}{\tau}T_i'+KQ' ]
参数(a)、(\tau)和(K)概括了流量、热容和传热行为的影响。该形式适用于推导传递函数、估算响应时间以及设计反馈控制器。
将模型与中试装置数据进行比较
学生可以施加受控扰动,记录液位和温度测量值,并将测量曲线与模型预测结果进行比较。差异可能揭示传感器滞后、混合不完全、热损失、物性变化、执行器限制或热负荷估计不准确等问题。
模型验证应同时包括瞬态和稳态检查。可信的模型应重现初始条件、趋近正确的稳态,并满足控制衡算方程。
理解权衡
简化假设
基本模型通常假设完全混合、密度恒定、热容恒定、热损失可忽略以及储罐温度均匀。它还可能将储罐面积和流量视为已知且恒定。
这些假设使模型适用于教学和初步控制设计,但对于高精度预测或运行条件变化较大的情况可能并不充分。
耦合和非线性行为
能量方程通过(Ah)与液位耦合。如果液位发生较大变化,温度方程将包含时变热储存项,系统的行为可能不再像固定的一阶系统。
假设液位恒定的模型对于温度控制研究可能很方便,但在使用前应检查液位变化相对于运行液位是否足够小。
热损失和执行器动态
实际储罐会通过罐壁、管道和周围环境散失热量。蒸汽阀和加热盘管也具有动态特性,使指令输入与实际传递给液体的热量有所不同。
忽略这些影响可能会使预测的温度响应比测量响应更快或更理想。当模型精度或控制器鲁棒性有要求时,应将这些因素加入模型。
模型过度或不足
仅使用液位衡算不足以实现温度控制,因为它忽略了热储存。相反,当中试装置处理的是物性近似恒定的非反应性液体时,加入组分、反应或详细焓衡算又是不必要的。
适当的模型取决于控制目标、运行范围以及实验所要求的精度。
根据目标做出正确选择
模型应与实验的受控变量、扰动以及所需的预测精度相匹配。
- 如果主要关注液位控制:使用质量衡算(A,dh/dt=F_i-F),并将(h)视为主要动态状态。
- 如果主要关注温度控制:同时加入质量衡算和能量衡算,因为液体体积决定热储存,而(T)是受控热状态。
- 如果主要关注控制器设计:在稳态运行点附近对能量衡算进行线性化,并使用偏差变量估算响应时间和过程增益。
- 如果主要关注模型验证:在已知扰动下,将仿真的液位和温度瞬态与中试装置数据进行比较,并在必要时考虑热损失、传感器滞后和执行器行为。
守恒衡算将搅拌罐加热器从一组设备转变为一个可预测的动态系统,使其液位和温度响应能够被分析、测试和控制。
总结表:
| 模型要素 | 控制关系 | 过程洞见 |
|---|---|---|
| 液位状态,h | A dh/dt = F_i - F | 流入和流出决定液位变化。 |
| 温度状态,T | Ah dT/dt = F_i(T_i - T) + Q/(rho c_p) | 进料条件和加热器热负荷决定温度动态。 |
| 液位稳态 | F_i = F | 当入口和出口流量相等时,储罐液位保持不变。 |
| 温度稳态 | T = T_i + Q/(rho c_p F_i) | 稳态温度取决于进料温度、流量和热负荷。 |
| 典型输入 | F_i, F, Q, T_i | 流量、蒸汽加热和进料温度可作为输入或扰动。 |
| 模型应用 | 仿真、线性化、控制、验证 | 该模型支持瞬态分析以及与中试装置数据的比较。 |
| 关键假设 | 完全混合、物性恒定、热损失可忽略 | 为了提高精度,可能还需要加入热损失、传感器和执行器动态。 |
需要使用可靠设备开展、测试或扩展过程动态教学吗?联系 LABPARK,了解面向化学工程、生物过程与生物技术以及环境与水处理的教育和职业培训单元操作中试装置。LABPARK 帮助大学、研究机构和企业开展实践培训、评估控制策略、验证质量和能量衡算模型,并通过根据其技术和教育目标定制的中试规模系统,展示真实的过程行为。