主权项 |
一种MIMO非最小相位CSTR反应器的非线性控制器设计方法,其特征在于,包括以下步骤:(1)根据实际的化工生产过程对CSTR进行状态空间建模,得到系统的非线性数学模型;(2)将系统的非线性数学模型进行状态反馈线性化,得到一般形式的MIMO线性化标准型;(3)将MIMO线性化标准型进行拆分,将其分解为若干线性子系统部分和一个SISO线性化标准型子系统部分,对各子系统单独设计控制器;(4)为了确保所设计的控制器能有效控制具有非最小相位特性的SISO线性化标准型子系统,给出了该控制器的一般形式,控制器中的待定参数通过极点配置和李雅普诺夫稳定性理论来确定;所述步骤(1)中,实际的化工生产过程所选用的是利用CSTR生产环戊烯这一化工过程,建立其状态空间模型如下式所示:状态方程:<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mfrac><mrow><msub><mi>dC</mi><mi>A</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mo>-</mo><msub><mi>k</mi><mn>1</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msub><mi>C</mi><mi>A</mi></msub><mo>-</mo><msub><mi>k</mi><mn>3</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msubsup><mi>C</mi><mi>A</mi><mn>2</mn></msubsup><mo>+</mo><mo>(</mo><msub><mi>C</mi><mrow><mi>A</mi><mn>0</mn></mrow></msub><mo>-</mo><msub><mi>C</mi><mi>A</mi></msub><mo>)</mo><mfrac><mi>F</mi><mi>V</mi></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>dC</mi><mi>B</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><msub><mi>k</mi><mn>1</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msub><mi>C</mi><mi>A</mi></msub><mo>-</mo><msub><mi>k</mi><mn>2</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msub><mi>C</mi><mi>B</mi></msub><mo>-</mo><msub><mi>C</mi><mi>B</mi></msub><mfrac><mi>F</mi><mi>V</mi></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><mi>d</mi><mi>T</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mo>(</mo><mo>-</mo><msub><mi>ΔH</mi><mn>1</mn></msub><mo>)</mo><msub><mi>k</mi><mn>1</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msub><mi>C</mi><mi>A</mi></msub><mo>+</mo><mo>(</mo><mo>-</mo><msub><mi>ΔH</mi><mn>2</mn></msub><mo>)</mo><msub><mi>k</mi><mn>2</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msub><mi>C</mi><mi>B</mi></msub></mrow><mrow><msub><mi>ρC</mi><mi>ρ</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>(</mo><mo>-</mo><msub><mi>ΔH</mi><mn>3</mn></msub><mo>)</mo><msub><mi>k</mi><mn>3</mn></msub><msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mi>A</mi></msub><msubsup><mi>C</mi><mi>A</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>Q</mi><mi>H</mi></msub></mrow><mrow><msub><mi>ρC</mi><mi>ρ</mi></msub></mrow></mfrac><mo>+</mo><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>-</mo><mi>T</mi><mo>)</mo><mfrac><mi>F</mi><mi>V</mi></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001077070860000011.GIF" wi="1574" he="463" /></maths>系统的被控输出为:<maths num="0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mi>T</mi></mtd></mtr><mtr><mtd><msub><mi>y</mi><mn>2</mn></msub><mo>=</mo><msub><mi>C</mi><mi>B</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001077070860000012.GIF" wi="190" he="150" /></maths>其中,C<sub>A0</sub>为物质A'的初始浓度,T<sub>0</sub>为反应初始温度,T是反应温度,C<sub>A</sub>、C<sub>B</sub>分别是物质A'和B'的浓度,F是物质A'的输入流量,V是CSTR的常容积;控制输入u<sub>1</sub>=F/V为稀释率,u<sub>2</sub>=Q<sub>H</sub>为单位体积增加或减少的热量;C<sub>ρ</sub>是热容量,ρ是混合物的密度;ΔH<sub>i</sub>是反应热;速率系数k<sub>i</sub>(T)由下列的Arrhenius方程给出:<img file="FDA0001077070860000021.GIF" wi="598" he="119" />i=1,2,3分别对应CSTR三个反应过程,其中,各过程的反应时间常数为k<sub>i0</sub>,反应激活能量为E<sub>i</sub>,R为气体常数;控制目标是:通过调节被控输入u<sub>1</sub>=F/V和u<sub>2</sub>=Q<sub>H</sub>来使输出y<sub>1</sub>=T和y<sub>2</sub>=C<sub>B</sub>分别稳定在设定的平衡点T<sub>S</sub>、C<sub>BS</sub>上,此时物质A'的浓度相应记为C<sub>AS</sub>,(C<sub>AS</sub>,C<sub>BS</sub>,T<sub>S</sub>)为系统达到稳态时的平衡点;所述步骤(4)中,一般形式的控制器设为v=‑Gx+v<sub>NL</sub>其中,x为SISO线性化标准型子系统的状态向量,G为行增益向量,v<sub>NL</sub>是为使系统内部零动态稳定而引入的非线性补偿项。 |