现代控制理论大作业概要

现代控制理论(主汽温对象模型)班级:学号:姓名:目录一.背景及模型建立1.火电厂主汽温研究背景及意义2.主汽温对象的特性3.主汽温对象的数学模型二.分析1.状态空间表达2.化为约当标准型状态空间表达式并进行分析3.系统状态空间表达式的求解4.系统的能控性和能观性5.系统的输入输出传递函数6.分析系统的开环稳定性7.闭环系统的极点配置8.全维状态观测器的设计9.带状态观测器的状态反馈控制系统的状态变量图10.带状态观测器的闭环状态反馈控制系统的分析三.结束语1.主要内容2.问题及分析3.评价一.背景及模型建立1.火电厂主汽温研究背景及意义火电厂锅炉主汽温控制决定着机组生产的经济性和安全性。由于锅炉的蒸汽容量非常大、过热汽管道很长，主汽温调节对象往往具有大惯性和大延迟，导致锅炉主汽温控制存在很多方面的问题，影响机组的整个工作效率。主汽温系统是表征锅炉特性的重要指标之一，主汽温的稳定对于机组的安全运行至关重要。其重要性主要表现在以下几个方面：(1)汽温过高会加速锅炉受热面以及蒸汽管道金属的蠕变，缩短其使用寿命。例如，12CrMoV钢在585℃环境下可保证其应用强度的时间约为10万小时，而在595℃时，其保证应用强度的时间可能仅仅是3万小时。而且一旦受热面严重超温，管道材料的强度将会急剧下降，最终可能会导致爆管。再者，汽温过高也会严重影响汽轮机的汽缸、汽门、前几级喷嘴和叶片、高压缸前轴承等部件的机械强度，从而导致设备损坏或者使用年限缩短。(2)汽温过低，会使得机组循环热效率降低，增大煤耗。根据理论估计可知：过热汽温每降低10℃，会使得煤耗平均增加0.2%。同时，汽温降低还会造成汽轮机尾部的蒸汽湿度增大，其后果是，不仅汽轮机内部热效率降低，而且会加速汽轮机末几级叶片的侵蚀。此外，汽温过低会增大汽轮机所受的轴向推力，不利于汽轮机的安全运行。(3)汽温变化过大会使得管材及有关部件产生疲劳，此外还将引起汽轮机汽缸的转子与汽缸的胀差变化，甚至产生剧烈振动，危及机组安全运行。据以上所述，工艺上对汽温控制系统的质量要求非常严格，一般控制误差范围在±5℃。主汽温太高会缩短管道的使用寿命，太低又会降低机组效率。所以必须实现汽温系统的良好控制。而汽温被控对象往往具有大惯性、大延时、非线性，时变一系列的特性，造成对象的复杂性，增加了控制的难度。现代控制系统中有很多关于主汽温的控制方案，本文我们着重研究带状态观测器的状态反馈控制对主汽温的控制[1]。2.主汽温对象的特性2.1主汽温对象的静态特性主汽温被控对象的静态特性是指汽温随锅炉负荷变化的静态关系。过热器的传热形式、结构和布置将直接影响过热器的静态特性。现代大容量锅炉多采用对流过热器、辐射过热器和屏式过热器。对流过热器布置在450℃～1000℃烟气温度的烟道中，受烟气的横向和纵向冲刷，烟气以对流方式将热量传给管道。而辐射过热器则是直接吸收火焰和高温烟气的辐射能。屏式过热器布置在炉膛内上部或出口处，属于辐射或半辐射过热器。对于辐射过热器，当锅炉负荷增加时，必须增加燃料量和风量。随着炉膛温度的提高，辐射传热量也将增加。但炉膛温度提高得并不多，而蒸汽量与燃料量是成正比的，这样就造成辐射传热量的增加赶不上蒸发量的增加。再者来说，当负荷增加时，强化燃烧后炉膛出口烟气温度将升高，这表明每公斤燃料产生的烟气带出炉膛的热量增多，这也说明了炉膛辐射吸热量的相对减少。所以，辐射受热面的汽温是随着锅炉负荷的增加而降低的。对于对流式受热面，当锅炉负荷增加和燃料消耗量增大时，流过对流受热面的烟气流速也增加，从而使对流的放热系数增大。另外，因为炉膛出口烟温升高，即进入对流受热面的吸热量增加值，超过了流过对流受热面的蒸汽流量的增加值。所以，对流受热面的汽温是随锅炉负荷的增加而升高的。2.2主汽温对象的动态特性[3]2.2.1蒸汽流量扰动下主汽温的动态特性在蒸汽流量D产生阶跃扰动时，主汽温变化的响应曲线如图2-2-1所示：图2-2-1蒸汽流量扰动下主汽温响应曲线蒸汽流量D扰动下，主汽温调节对象动态特性的特点是:有延迟、有惯性、有自衡能力。当锅炉负荷增加时，通过对流式过热器的烟气温度和流速都增加，因此对流式过热器出口汽温升高；但对于辐射式过热器，炉膛内烟温升高增加的辐射传热量小于蒸汽流量增加所需的吸热量，因此辐射式过热器出口汽温下降。2.2.2烟气热量扰动下汽温的动态特性在烟气热量yQ产生阶跃扰动下，汽温θ变化的响应曲线如图2-2-2所示：图2-2-2烟气热量扰动下汽温响应曲线烟气热量yQ扰动下，汽温调节对象动态特性的特点是:有迟延、有惯性、有自平能力。由于烟气热量变化时，沿过热器长度使烟气和过热蒸汽之间的传热量同时变化，因此，汽温θ反应较快，其时间常数Tc和迟延时间τ均较小。2.2.3减温水量扰动下汽温的动态特性减温水量BW扰动下，汽温调节对象动态特性的特点是:有迟延、有惯性、有自衡能力。由于现代大型锅炉过热器管路很长，因此减温水量BW变化时汽温θ反应较慢，其时间常数Tc和迟延时间τ均较大。在减温水量BW产生阶跃扰动下，汽温θ变化的响应曲线如图2-2-3所示：图2-2-3减温水量扰动下汽温响应曲线3.主汽温对象的数学模型3.1主汽温对象模型的建立分析和设计自动控制系统的一个首要任务是建立系统的数学模型，因为不论要了解的是简单系统的特性还是复杂系统的特性，都必须掌握系统中各变量之间的相互动态关系。尽管，电厂汽温控制系统的动态特性复杂，具有多变量、非线性和分布复杂的特点，难以建立精确的数学模型，但是运用多种知识，建立能相对反映汽温控制系统动态性能的数学模型，还是对研究相应的自动控制方法大有益处。锅炉过热器是由辐射过热器、对流过热器和减温器等组成，其任务是将从汽包出来的饱和蒸汽加热到一定的数值，然后送往汽轮机去做功。通常称减温器前的过热器为前级过热器，减温器后的过热器为后级过热器。图2-3-1过热蒸汽喷水减温系统示意图过热器布置在高温烟道中，大型锅炉的过热器往往分为若干段，在各段之间设置喷水减温器，即采用汽温的分段控制，温度调节用减温水由锅炉的给水系统提供。其示意图如图2-3-1。图中，θs为过热器出口蒸汽温度，它是控制系统的被调量；Ws1、Ws2是减温器的喷水量，它是控制系统的调节量。3.2主汽温对象的模型目前电厂中广泛使用的汽温调节系统主要是基于常规控制器构成的串级控制系统。因为常规PID控制器的参数是根据被控对象的数学模型来整定的。而汽温调节对象的动态特性随机组负荷的变化而有较大变化，因此，当机组负荷变化时，该控制系统的工作点建立多固定模型，设计相应的控制器，在运行过程中根据系统工况变化选择相应的控制器输出以获得满意的控制效果，通过控制器在不同工况下的切换使系统达到全局稳定性的要求。某机组导前区模型：11()1.1()()121ysGsyss，惰性区模型：122()2.8()()(301)ysGsuss，总体模型：2()1.12.8()()(121)30s+1)ysGsuss（。对应的一种状态空间模型为1/301/300001/301/300001/121.1*2.8/12xxuyx100二.分析1.状态空间表达式1/301/300001/301/300001/121.1*2.8/12xxuyx1002.化为约当标准型状态空间表达式并进行分析程序:clear;A=[-1/301/300;0-1/301/30;00-1/12];B=[001.1*2.8/12]';C=[100];D=0;[P,J]=jordan(A);A1=inv(P)*A*PB1=inv(P)*BC1=C*P结果:变换矩阵000000.16667.006667.04444.00222.04444.0P特征值0333.00000333.00000833.0J0333.0000000.10333.00000833.0~AB~2567.002567.0C~4444.00222.04444.0因此约当标准型状态空间表达式为uxx2567.002567.0~0333.0000000.10333.00000833.0~xy~4444.00222.04444.0系统矩阵A的特征根有重根,因此,经线性变换,将A化为约当标准型.3.系统状态空间表达式的求解线性定常非齐次状态方程的解可表示为tdButxttx0)()()0(Atet矩阵A有重特征根时,10333.00333.00333.00833.000000PeteeePettttAt0000.15000.1000000.300000.450000.10000000000000.16667.006667.04444.00222.04444.00333.00333.00333.00833.0tttteteeetttttttttteeeeeteeetee0833.00333.00833.00333.00333.00333.00833.00333.00333.00333.0006667.06667.000005.104444.02222.04444.00006.00333.0999.0当输入脉冲函数)()(ttu时,则系统状态方程的解Btxttx)()0()()(因此tttttteeeetee0833.00333.00833.00333.00333.00833.02567.01711.01711.01141.00570.01141.04.系统的能控性和能观性(1)能控性①代数判据:线性定常系统状态完全能控的充要条件为:矩阵BAABBUnC1...是满秩的,或表示成nBAABBrankrankUnC1...②模态判据1:若系统矩阵A为对角型,则系统能控的充要条件是输入矩阵B没有任何一行的元素全部为零③模态判据2:若系统矩阵为约当型,则系统能控的充要条件是:输入矩阵B中对应于互异的特征值的各行,没有一行的元素全为零;输入矩阵B中与每个约当块最后一行相对应的各行,没有一行的元素全为零.Qc=ctrb(A,B)rank(Qc)0018.00214.02567.00010.00086.000003.000Qc3ans因此系统能控(2)能观性①代数判据:系统定常系统),(CA状态能观测的充要条件为:TnOCACACV1...的秩为n②模态判据1:若系统A为对角型,则系统能观测的充要条件是输出矩阵C没有任何一列的元素全部为零③模态判据2:若系统A为约当型,则系统能观测的充要条件是:)0(006667.06667.000005.104444.02222.04444.00006.00333.0999.0)(0833.00333.00833.00333.00333.00333.00833.00333.00333.00333.0xeeeeeteeeteetxtttttttttt输出矩阵C中对应于互异的特征值的各列,没有一列的元素全为零;输出矩阵C中与每个约当块的第一列相对应的各列,没有一列的元素全为零Qo=obsv(A,C)rank(Qo)0011.00022.00011.000333.00333.0000000.1Qo3ans因此系统能观5.系统的输入输出传递函数G=ss(A,B,C,D);Gc=tf(G)Tran