化学反应装置内部模型控制器的设计

这个例子展示了如何使用控制系统设计器在串联化学反应器的内模控制结构中设计一个补偿器。在过程控制应用中,基于模型的控制系统通常用于跟踪设定值和拒绝负载干扰。

工厂模型

这个例子中的工厂是一个化学反应器系统,由两个混合良好的容器组成。

反应器为等温反应器,每个反应器中对A组分的反应为一级反应:

$$ R_ {A} = -kC_ {A} $$

物料平衡被应用到系统生成系统的动态模型。罐水平假定留下来,因为溢流喷嘴的恒定,因此没有涉及电平控制。

有关此植物的详细信息,请参见实施例3.3中“过程控制:设计过程和控制系统的动态性能”的第3章由Thomas E.马林。

以下微分方程描述组件结余:

$$ V \压裂{DC_ {A1}} {DT} = F({C_ A0} -C_ {A1}) -  {VkC_ A1} $$

$ $ V \压裂{dC_ {A2}} {dt} = F (C_ {A1} -C_ {A2})——VkC_ {A2} $ $

在稳定状态下,

$$ \压裂{DC_ {A1}} {DT} = 0 $$

$$ \压裂{DC_ {A2}} {DT} = 0 $$

物料平衡如下:

$$˚F^ *(C_ {A0} ^ *  -  C_ {A1} ^ *) -  VkC_ {A1} ^ * = 0 $$

$$F^*(C_{A1}^* - C_{A2}^*) - VkC_{A2}^* = 0$$

哪里$ {C_ A0} ^ * $$ {C_ A1} * $,$ {C_ A2} * $是稳态值。

替代,以下设计规范和反应器参数:

$$˚F^ * $$ = 0.085摩尔$$ /分钟$$

$$C_{A0}^*$$ = 0.925 $$mol/min$$

$$V$$ = 1.05 $$m^3$$

$$k$$ = 0.04 $$min^{-1}$$

在两个反应器将所得的稳态浓度是:

$$ C_ {A1} ^ * = {KC_ A0} ^ * = 0.6191摩尔/立方公尺$$

$$ C_ {A2} ^ * = K ^ 2C_ {A0} ^ * = 0.4144摩尔/立方公尺$$

哪里

$$ K = \压裂{F ^ {*}} {F * + VK} = 0.6693 $$

对于这个例子,设计一个控制器,用于从所述第二反应器维持反应物的出口浓度,美元C_ {A2} ^ * $在任何干扰的在进料浓度的存在下,$ {C_ A0} $。被操作的变量是反应物的摩尔流量,F在进入第一反应器。

线性装置模型

在此控制设计问题中,被控对象模型为

$$ \压裂{C_ {A2}(S)} {F(S)} $$

扰动模型为

$ $ \压裂{C_ {A0} (s)} {C_ {A2} (s)} $ $

该化学过程可以用下面的框图来表示:

哪里

$$ G_ {A1} = \压裂{C_ {A1}(S)} {C_ {A0}(S)} = \压裂{0.6693} {8.2677s + 1} $$

$ $ G_ {F1} = \压裂{C_ {A1} (s)} {F (s)} = \压裂{2.4087}{8.2677 s + 1} $ $

$$ G_ {A2} = \压裂{C_ {A2}(S)} {C_ {A1}(S)} = \压裂{0.6693} {8.2677s + 1} $$

$$ G_ {F2} = \压裂{C_ {A2}(S)} {F(S)} = \压裂{1.6118} {8.2677s + 1} $$

根据框图,得到植物模型和干扰模型如下:

$ $ \压裂{C_ {A2} (s)} {F (s)} = G_ {F1} G_ {A2} + G_ {F2} = \压裂{13.3259 s + 3.2239} {(8.2677 s + 1) ^ 2} $ $

$ $ \压裂{C_ {A2}} {C_ {A0}} = G_ {A1} G_ {A2} = \压裂{0.4480}{(8.2677 s + 1) ^ 2} $ $

创建在命令行的工厂模型:

S = TF('S');G1 =(13.3259 * S + 3.2239)/(8.2677 * S + 1)^ 2;G2 = G1;GD = 0.4480 /(8.2677 * S + 1)^ 2;

G1是在控制器评估中使用的实际设备。G2是在实际设备的近似,它是用来作为IMC结构的预测模型。G2 = G1意味着没有模型失配。是扰动模型。

在控制系统设计中定义内模控制结构

开放式控制系统设计。

controlSystemDesigner

选择IMC控制架构。在控制系统设计中,单击编辑架构。在编辑结构对话框中,选择配置5。

加载系统数据。对于G1G2,,指定一个模型价值

调整补偿器

情节的开环阶跃响应G1

步骤(G1)

右键单击该绘图并选择特点>上升时间子菜单。点击蓝色的上升时间标记。

上升时间是25秒,我们要调整IMC补偿,以实现更快的闭环响应时间。

要调整IMC补偿,在控制系统设计中,单击调优方法,并选择内模控制(IMC)调谐

选择一个占主导地位的闭环时间常数2和一个所需的控制器顺序2

要查看闭环阶跃响应,在控制系统设计,双击IOTransfer_r2y:步骤图选项卡。

模型失配的控制性能

当设计的控制器,我们假设G1等于G2。在实践中,他们往往是不同的,控制器必须足够强大的跟踪设定值和拒绝干扰。

在G1和G2之间创建模型不匹配,并在设置点变化和负载干扰存在的情况下检查MATLAB命令行控制性能。

导出IMC补偿到MATLAB的工作空间。请点击出口。在导出模型对话框中,选择补偿器模型C

请点击出口

转换IMC结构的经典反馈控制结构与所述前馈路径和反馈单元的控制器。

C = ZPK([ -  0.121 -0.121],[ -  0.242,-0.466],2.39);C_new =反馈(C,G2,+ 1)
C_new = 2.39 (s + 0.121) ^ 4——(s - 0.0001594) (s + 0.121) (s + 0.1213) (s + 0.2419)连续时间零/钢管/增益模型。

定义如下工厂模型:

  • 无模型失配:

G1p = (13.3259 * s + 3.2239) / (8.2677 * s + 1) ^ 2;
  • G1时间常数变化5%:

G1t =(13.3259 * S + 3.2239)/(8.7 * S + 1)^ 2;
  • G1增益增加3倍:

G1G = 3 *(13.3259 * S + 3.2239)/(8.2677 * S + 1)^ 2;

评估设定值跟踪性能。

步骤(反馈(G1p * C_new, 1),反馈(G1t * C_new, 1),反馈(G1g * C_new, 1))传说(“没有模型不匹配”“时间常数不匹配”“增益不匹配”

评估干扰抑制性能。

步骤(GD *反馈(1,G1P * C_new),钆*反馈(1,G1t * C_new),钆*反馈(1,G1G * C_new))图例(“没有模型不匹配”“时间常数不匹配”“增益不匹配”

控制器是相当稳健的工厂参数的不确定性。

也可以看看

相关的话题