增益预定化学反应器的控制

打开脚本

该示例显示了如何设计和调整从低至高转化率过渡的化学反应器的增益调度控制器。用于背景，见Seborg，D.E.等人，“过程动态和控制”，第2 ED。，2004，Wiley，PP。34-36。

连续搅拌槽式反应器

这里考虑的过程是在从低至高转化率的转变期间连续搅拌釜反应器（CSTR）（高于低残留浓度）。由于化学反应是放热（产生热量），因此必须控制反应器温度以防止热失控。随着过程动态是非线性的，并且随着转换速率的增加，过程动态是从过程动态的非线性和从稳定的转换而转换。反应堆动态在Simulink中建模。金宝app受控变量是残余浓度Cr以及反应堆温度TR.，而被操纵的变量是温度TC.在反应器的冷却夹套中循环的冷却剂的内容。

Open_System（“rct_CSTR_OL”）

我们想从最初的残留浓度8.57 kmol/m^3降到2 kmol/m^3。了解工艺动态如何随着残留浓度的变化而变化Cr，求五个值的平衡条件Cr在8.57和2之间，并在每个平衡周围线性化过程动力学。记录反应堆和冷却剂在每个平衡点的温度。

Creq = Linspace（8.57,2,5）';%浓度TrEQ = 0(5、1);%反应堆温度TcEQ = 0(5、1);％冷却剂温度％指定修剪条件opspec = operspec (“rct_CSTR_OL”5);为k = 1：5％设定所需的残余浓度opspec (k) .Outputs(1)。y = CrEQ (k);opspec (k) .Outputs(1)。知道= true;结束计算平衡条件并记录相应温度(op、报告)= findop (“rct_CSTR_OL”opspec,．..findopOptions (“DisplayReport”，'离开'））;为(2).y =1;TcEQ (k) = op (k) .Inputs.u;结束修剪条件下的％线性化过程动态g =线性化（“rct_CSTR_OL”，“rct_CSTR_OL /装运箱”, op);G.InputName = {Cf的，“助教”，“Tc”};G.OutputName = {“Cr”，“Tr”};

绘制反应堆和冷却剂的平衡温度作为浓度的函数。

次要情节(311),情节(CrEQ“b - *”)、网格、标题(的残余浓度), ylabel (“CrEQ”)次要情节(312)、情节(TrEQ“b - *”)、网格、标题(“反应堆温度”), ylabel ('treq'）子图（313），绘图（TCEQ，“b - *”)、网格、标题(冷却液温度的), ylabel ('tceq'）

一种开环控制策略包括遵循冷却剂温度分布以上平滑过渡之间Cr= 8.57,Cr= 2平衡。然而，这一策略注定要失败，因为反应在中程是不稳定的，必须适当冷却以避免热失控。这是通过考察上述五个均衡点的线性化模型的极点来证实的(五个模型中有三个是不稳定的)。

极(G)

Ans (:，:，3) = 0.5114 + 0.00000 i -0.8229 + 0.00000 i Ans (:，:，4) = 0.0453 + 0.00000 i -0.4991 + 0.00000 i Ans (:，:，5) = -1.1077 + 1.0901i -1.1077 - 1.0901i

BODE PLOT进一步突出显示从转换时植物动力学的显着变化Cr= 8.57Cr= 2。

波德(G (clf:,“Tc”), {0.01, 10})

反馈控制策略

为了防止热失控而逐步降低残留浓度，使用反馈控制来调整冷却剂温度TC.基于残留浓度的测量Cr和反应堆温度TR.．对于这个应用，我们使用一个串级控制架构，其中内环调节反应堆温度，而外环跟踪浓度设置。两个反馈回路都是数字化的，采样周期为0.5分钟。

Open_System（“rct_CSTR”）

目标集中Cref在T = 36的T = 10到2 kmol / m ^ 3以下降到8.57 kmol / m ^ 3（过渡持续26分钟）。相应的简档tr反应器的温度是通过对平衡值的插值得到的TrEQ从分析。控制器计算冷却剂温度调节直接转矩相对于初始平衡值TCEQ（1）= 297.98Cr= 8.57。注意，模型是这样设置的，最初，输出TrSP“浓度控制器”块与反应器温度相匹配，调整直接转矩为零，冷却液温度TC.是否处于均衡值TCEQ（1）．

CLF t = [0 10:36 45];C = interp1([0 10 36 45]，[8.57 8.57 2 2]，t);子图(211)，图(t,C)，网格，集合(gca，'ylim'，[0 10]）标题（“目标剩余浓度”), ylabel ('奶油') subplot(212)， plot(t,interp1(CrEQ,TrEQ,C)) title(“平衡时相应的反应堆温度”), ylabel (“不洁净的”)、网格

控制目标

使用TuningGoal对象来捕获设计需求。首先,Cr应该遵循定位点Cref响应时间约为5分钟。

R1 = TuningGoal。跟踪('奶油'，“Cr”5);

内环(温度)应该以足够的阻尼和足够快的衰减来稳定反应动力学。

MinDecay = 0.2;MinDamping = 0.5;%约束内环的闭环极点，使外环打开R2 = TuningGoal。波兰人(“Tc”，思维，造版）;r2.openings =.“TrSP”；

控制器输出的速率限制块指定冷却液温度TC.变化速度不能超过每分钟10度。这是对控制器权限的严重限制，如果忽略它，可能会导致性能低下或不稳定。要考虑这个利率限制，请注意这一点Cref以0.25 kmol / m ^ 3 / min的速率变化。以确保TC.没有比10度/分钟更大，收益不变Cref来TC.应小于10/0.25=40。

R3 = TuningGoal。获得('奶油'，“Tc”, 40);

最后，在工厂输入中需要至少7 dB的增益余量和45度的相余量TC.．

R4 = TuningGoal。利润(“Tc”7、45);

增长计划的控制器

为了实现这些要求，我们在外环采用PI控制器，在内环采用超前补偿器。由于取样速率缓慢，需要铅补偿器来充分稳定化学反应在中等浓度Cr= 5.28 kmol / m ^ 3 /分钟。因为反应动力学在很大程度上随浓度而变化，我们进一步将控制器增益设定为浓度的函数。这是在Simulink中使用查找表块建金宝app模的，如图1和2所示。

图1：用于集中循环的增益预定PI控制器。

图2：温度循环的增益预定铅补偿器。

调优这个增益调度控制器相当于在一个集中值范围内调优查询表数据。与其调优各个查询表项，不如将控制器增益参数化Kp, Ki, Kt, a, b作为二次多项式Cr例如,

$k_p（c_r）= k_ {p0} + k_ {p1} c_r + k_ {p2} c_r ^ 2。$

除了减少调整的变量数量之外，这种方法确保了平滑的增益转换为Cr各不相同。使用systune，您可以自动调整系数 $K_{p0}， K_{p1}， K_{p2}， K_{i0}， \ldots$ 满足要求R1-R4在上面计算的五个均衡点。这增加了在五个设计点上调了增益预定的控制器Cref弹道。使用tunableSurface对象将每个增益参数化为的二次函数Cr．“调谐网格”设置为五个浓度CrEQ二次参数化的基函数为 C_r美元,C_r ^ 2美元．大多数收益被初始化为相同为零。

TuningGrid =结构(“Cr”, CrEQ);ShapeFcn = @(Cr) [Cr, Cr^2];Kp = tunableSurface (“金伯利进程”， 0, TuningGrid, ShapeFcn);Ki = tunableSurface (“吻”，-2，tuninggrid，shapefcn）;kt = tunballulface（“Kt”， 0, TuningGrid, ShapeFcn);一个= tunableSurface (“一个”， 0, TuningGrid, ShapeFcn);b = tunableSurface (“b”， 0, TuningGrid, ShapeFcn);

控制器调优

由于目标带宽在奈奎斯特频率的10年内，因此更容易在离散域直接调整控制器。将线性化的过程动态离散，采样时间为0.5分钟。使用ZOH方法来反映数字控制器如何与连续时间工厂相互作用。

t = 0.5;Gd =汇集(G, Ts);

创建一个slTuner用于调整上面介绍的二次增益计划的接口。采用块代换法将非线性对象模型替换为5个离散线性模型Gd在设计点获得CrEQ．使用setBlockParam要关联可调谐增益函数Kp，Ki，Kt，一个，b使用同名的查找表块。

BlockSubs =结构(“名字”，“rct_CSTR /装运箱”，“价值”、Gd);ST0 = slTuner (“rct_CSTR”，{“金伯利进程”，“吻”，“Kt”，“一个”，“b”}, BlockSubs);ST0。Ts = Ts;用于调优的样本时间%登记兴趣点ST0.addPoint ({'奶油'，“Cr”，“Tr”，“TrSP”，“Tc”}）参数化查找表块ST0.setBlockParam (“金伯利进程”Kp);ST0.setBlockParam (“吻”、Ki);ST0.setBlockParam (“Kt”、Kt);ST0.setBlockParam (“一个”,);ST0.setBlockParam (“b”, b);

你现在可以使用systune根据要求调整控制器系数R1-R4．使稳定性保证金要求硬限制并优化其余要求。

ST = systune(ST0，[R1 R2 R3]，R4);

Final: Soft = 1.23, Hard = 0.99991, Iterations = 195

由此产生的设计满足了硬约束（硬< 1)，几乎满足其余的要求(软为了验证这一设计，模拟浓度相同的斜坡的响应Cref．每个绘图显示五个设计点的线性响应CrEQ．

T = 0：TS：20;UC = Interp1（[02 5 20]，（ -  0.25）* [0 0 3 3]，t）;子图（211），LSIM（Getiotransfer（St，'奶油'，“Cr”），UC）网格，集（GCA，'ylim'[-1.5 - 0.5])、标题(的残余浓度)次要情节(212),lsim (getIOTransfer(圣,'奶油'，“Tc”)、加州大学)网格,标题(冷却剂温度变化的）

注意，冷却剂温度的变化率保持在物理范围内（每分钟10度或每个样品周期5度）。

控制器验证

检查每个增益如何变化Cr在过渡期间。

%访问调整后的增益计划TGS = GetBlockParam（ST）;绘制增益曲线clf subplot(321)， viewSurf(TGS.Kp)， ylabel(“金伯利进程”）子图（322），ViewsURF（TGS.KI），Ylabel（“吻”）子图（323），ViewsURF（TGS.KT），Ylabel（“Kt”) subplot(324)， viewSurf(TGS.a)， ylabel(“一个”) subplot(325)， viewSurf(TGS.b)， ylabel(“b”）

为了在Simulink中验证增益调度控制器，首先使用金宝appwriteBlockValue将调整结果应用于Simulink模型。金宝app对于每个查找表块，这会评估表断点处的相应二次增益公式，并相应地更新表数据。

writeblockvalue（st）

接下来，按下Play按钮，用调整后的增益计划模拟响应。仿真结果如图3所示。增益调度控制器成功地驱动反应通过过渡充分的响应时间和没有饱和的速率限制(控制器输出杜匹配有效温度变化直接转矩)．反应堆温度保持在它的平衡值附近tr，表示控制器保持反应接近平衡，同时防止热失控。

图3：使用增益预定的级联控制器转换。

控制器在Matlab中调整

或者，您可以在MATLAB中直接调整增益调度而不使用slTuner界面。首先将增益参数化为二次函数Cr正如上面做的。

TuningGrid =结构(“Cr”, CrEQ);ShapeFcn = @(Cr) [Cr, Cr^2];Kp = tunableSurface (“金伯利进程”， 0, TuningGrid, ShapeFcn);Ki = tunableSurface (“吻”，-2，tuninggrid，shapefcn）;kt = tunballulface（“Kt”， 0, TuningGrid, ShapeFcn);一个= tunableSurface (“一个”， 0, TuningGrid, ShapeFcn);b = tunableSurface (“b”， 0, TuningGrid, ShapeFcn);

使用这些增益来构建PI和先导控制器。

π= pid (Kp,吻,“t”Ts,'不时'，'min'）;π。u =ECr的；π。y =“TrSP”；LEAD = Kt * tf([1 -a]，[1 -b]，Ts，'不时'，'min'）;领先。u =“ETr”；领先。y =“Tc”；

使用连接在五个设计点上建立整体控制系统的闭环模型。标记控制器输出TrSP和TC.作为“分析点”，以便可以打开循环，并在这些位置评估稳定性裕度。闭环模型T0线性模型的5乘1数组是否取决于可调系数Kp, Ki, Kt, a, b．每个模型都是离散的，每半分钟采样一次。

gd.timeUnit =.'min'；S1 = sumblk ('ECR = CREF  -  CR'）;s2 = sumblk（'ETr = TrSP - Tr'）;T0 =连接(Gd (:,“Tc”)、铅、π,S1, S2,'奶油'，“Cr”，{“TrSP”，“Tc”}）;

最后，使用systune调整增益计划系数。

T = systune(T0，[R1 R2 R3]，R4);

Final: Soft = 1.23, Hard = 0.99987, Iterations = 193

结果与上述结果相似。通过绘制增益的函数来确定Cr使用调谐系数T．

clf subplot(321)， viewSurf(setBlockValue(Kp,T))， ylabel(“金伯利进程”) subplot(322)， viewSurf(setBlockValue(Ki,T))， ylabel(“吻”）Subplot（323），ViewsURF（SetBlockValue（kt，t）），ylabel（“Kt”) subplot(324)， viewSurf(setBlockValue(a,T))， ylabel(“一个”) subplot(325)， viewSurf(setBlockValue(b,T))， ylabel(“b”）