毕业论文火电厂锅炉燃烧过程控制系统仿真研究

1、 本 科 毕 业 论 文火电厂锅炉燃烧过程控制系统仿真研究 Simulation Study on CombustionControl System of Boiler in Thermal Power Plant 学院名称：电子信息与电气工程学院 专业班级： 自动化2010级2班 学生姓名： 学 号： 指导教师名称： 指导老师职称：2014年 5月目 录摘要IAbstractII引言1第一章 控制系统及MATLAB语言的应用基础31.1控制系统性能要求3 1.2 控制系统的时域分析3 1.3 控制系统的频域分析4 频域性能指标4 频域分析的三种分析方法4 1.4 控制系统的根轨迹分析5 1.

2、5MATLAB软件认识5 的特点5在控制系统分析中的应用5 1.5.3 根轨迹绘制6 控制系统的频域分析6 1.6 MATLAB环境下的Simulink简介9第二章 火电厂锅炉燃烧过程控制系统设计10 2.1 蒸汽压力控制系统和燃料空气比值控制系统12 2.2 炉膛负压控制系统13第三章 火电厂锅炉燃烧过程控制系统数学模型16 3.1 蒸汽压力控制系统和燃料空气比值控制系统数学模型16 3.2 炉膛负压控制系统数学模型16第四章 系统稳定性分析17 4.1 燃料控制系统17 4.2 空气流量控制系统18 4.3 炉膛负压控制系统18第五章 火电厂锅炉燃烧过程控制各子系统仿真20 5.1 燃料控

3、制系统20 5.2 蒸汽压力控制系统23 5.3 空气流量控制系统25 5.4 炉膛负压控制系统28第六章 火电厂锅炉燃烧过程控制系统Simulink仿真31结论35致谢36参考文献37火电厂锅炉燃烧过程控制系统仿真研究摘要：在现代化生产过程中，为保证生产安全顺利进行，达到优质高产，提高经济效益和劳动生产率，必须对生产过程的各种参数进行自动控制，这种控制就是过程控制。过程控制是自动化生产的重要方向和组成部分，已经被广泛应用在石油、化工、轻工、冶金、机械、电力等许多国民经济的重要领域。锅炉是工业生产中的动力设备，因而对锅炉进行有效的控制也就显得至关重要。燃烧是锅炉控制中必要的一环，从燃烧角度来看

4、，有燃气、燃油、燃煤的分别。燃烧过程的控制包括燃烧程度控制、温度控制、压力控制、节能控制、安全性控制等。本文根据火电厂锅炉燃烧过程生产工艺和流程，利用串级控制、比值控制等控制系统的特点，分别设计火电厂锅炉燃烧过程控制系统的三个子系统：蒸汽压力控制系统，燃料与空气比值控制系统，炉膛负压控制系统，并建立数学模型。在仿真软件MATLAB/Simulink中，根据控制系统的原理方框图，作出仿真模型图，分别进行相应的仿真。经PID参数整定后，得出仿真结果，并进行分析总结。关键词：燃烧过程控制；MATLAB仿真；SIMULINKSimulation Study on CombustionControl S

5、ystem of Boiler in Thermal Power PlantAbstract: in the modern production process, to ensure safety in production smoothly,achieve the high yield and quality, enhance economic efficiency and labor productivity, must carry out the automatic control of the various parameters of production process, this

6、 control is the process control. Process control is an importantdirection of automated production and component, has been widely used inpetroleum, chemical industry, light industry, an important field, metallurgy, machinery,electricity and many other national economy. Boiler is an dynamic equipmenti

7、n industrial production, so the effective control of boiler is crucial.Combustion is a necessary link in the boiler control, from the burning point of view,gas, fuel, coal respectively. Control of burning process including combustion control,temperature control, pressure control, energy saving contr

8、ol, safety control.According to theproduction processand the process of combustion of boiler in thermal power plant, the cascade control, ratio control system, the three subsystems were designed thermal power plant boiler combustion process control system: the steam pressure control system, ratio co

9、ntrol system of fuel and air, hearth negative pressure control system and the establishment of mathematical model. In thesimulation software MATLAB/Simulink, according to the principle of control systemblock diagram, make a simulation model, respectively, corresponding simulation. ThePID parameter t

10、uning, simulation results obtained, analyzing and summarizing.Keywords: combustion control; MATLAB simulation; SIMULINK引言过程控制系统是工业中控制系统的主要表现形式，一般指工业生产过程中自动控制系统的被控变量为温度、压力、流量、液位、成分等变量的系统。由于被控过程的多样性，因此过程控制系统的形式也多样，相应的控制方案也丰富多彩。由于负荷变化从本质上说是非线性和时变的，锅炉侧又存在很大的延迟，负荷适应能力和主蒸汽压力稳定的矛盾一直是锅炉燃烧自动控制中有待于进一步解决的问题。而且

11、锅炉燃烧控制系统不同于大多数生产过程控制系统，它不仅存在动态品质调节的问题还要考虑到锅炉的经济燃烧。如何使主蒸汽压力既具有良好的动态特性，又能使入炉燃料得以充分燃烧(即高效燃烧问题)，是燃烧优化控制的真正内涵，也是燃烧过程控制的关键。目前，主汽压调节系统的控制器采用的控制规律以常规的基于数学模型的PID控制为多。在实际工程中，控制系统的结构往往很复杂，如果不借助专用的系统建模软件，很难准确地把一个控制系统的复杂模型输入计算机，对其进行进一步的分析与仿真，Simulink的出现给控制系统分析与设计带来了福音。它有两个主要功能：Simu（仿真）和Link（连接），即该软件可以利用鼠标在模型窗口上搭

12、建出所需要的控制系统模型，然后对系统进行仿真和分析。Simulink仿真的模型是具体的函数模型，因此，过程的建模就显的十分的重要。过程建模方法主要有这么几种：1机理法建模用机理法建模就是根据过程的内在机理，写出各种有关的平衡方程，例如物质平衡方程、能量平衡方程、动量平衡方程、相平衡方程，以及反映流体流动、传热、传质、化学反应等基本规律的运动方程、物性参数方程和某些设备的特性方程等，从中获得所需的数学模型。2测试法建模测试法建模通常只用于建立输入输出模型。它是根据过程的输入和输出的实测数据进行某种数学运算后得到的模型，其主要特点是把被研究的过程视为一个黑匣子，完全从外特性上描述它的动态性质，也称

13、为“黑箱模型”。复杂过程一般都采用测试法建模。测试建模法又可分为经典辨识法和系统辨识法两大类：（1）经典辨识法不考虑测试数据中偶然性误差的影响，只需对少量的测试数据进行比较简单的数学处理，计算工作量一般较小。经典辨识法包括时域法、频域法和相关分析法。（2）系统辨识法其特点是可以清除测试数据中的偶然性误差即噪声的影响，为此就需要处理大量的测试数据，计算机是不可缺少的工具。本文通过搜集资料得到系统辨识函数并以此确立数学模型，搭建仿真框图，调整PID参数，最终使系统在较好的指标内达到稳定。第一章 控制系统及MATLAB语言的应用基础1.1 控制系统性能要求控制的任务是形成控制作用的变化规律，使得不管

14、是否存在扰动对象都能得到所期望的行为。自动控制系统是在无人直接操作或干预的条件下，通过控制器使控制对象自动的按照给定的规律运行，使被控制量能够按照给定的规律变化。对于一个控制系统的性能要求可以概括为三个方面：稳定性、快速性和准确性。稳定性是指动态过程的平稳性，就是动态过程的振荡倾向于系统受干扰偏离平衡状态后恢复平衡状态的能力。控制系统的稳定性是系统使用的基础，不稳定的系统没有实际意义。快速性是指动态过程的快速性，即动态过程的时间长短，时间越短，快速性能越好。准确性是指动态过程的最终精度，即控制系统输出量的期望值与实际终止值之差（即稳态误差），其差值越小精度越高。受控对象不同，对稳、准、快的技术

15、要求也各有所侧重。对于一个实际系统，不能要求三项性能指标都很优良，那样系统的成本会很高。而且，同一个系统的稳、准、快各个指标往往相互制约。1.2控制系统的时域分析时域分析法是根据自动控制系统微分方程，用拉普拉斯变换来求解动态响应的过程曲线。时域分析的目的是求解响应的性能指标。通常将控制系统跟踪或复现阶跃输入信号响应的指标作为系统控制性能的指标。阶跃响应的一般性能指标有：上升时间tr、峰值时间tp、超调量sigma(%)、调节时间ts及稳态误差ess。超调量%是指阶跃响应曲线h（t）中对稳态值的最大超出量与稳态值之比。即： （1-1）峰值时间tp：是指从0到阶跃响应曲线h(t)进入稳态值附近的误

16、差带而不再超出的最小时间。恢复时间：从阶跃扰动作用开始到输出量基本上恢复稳态的过程中，输出量与新的稳态值之差进入某基准量的范围之内所需的时间。在MATLAB中实现时域分析有两种方法：一种是MATLAB的函数指令方式，如求取连续系统的单位阶跃响应函数step()，单位冲激响应函数impulse()，零输入响应函数initial()等；另一种是在Simulink窗口菜单操作方式下进行时域仿真。1.3控制系统的频域分析频域性能指标峰值：幅频特性A()的最大值。频带：幅频特性A()的数值衰减到0.707A(0)时对应的频率。相频宽：相频特性(j)等于/2时对应的频率。频域分析的三种分析方法Bode图（

17、幅频/相频特征曲线）、Nyquist曲线、Nichols图。（1）Bode图：Bode图即对数频率特性曲线。（2）Nyquist曲线：Nyquist曲线是根据开环频率特性在复平面上会出的幅相轨迹。根据开环Nyquist曲线，可判断闭环系统的稳定性。反馈控制系统稳定的充分条件是，Nyquist曲线按逆时针包围临界点（-1，j0）的圈数p等于开环传递函数位于右半S平面的极点数，否则闭环系统不稳定。当开环传递函数包含虚轴上的极点时，闭合曲线应以0的半圆从右侧绕过该极点。（3）Nichols图：Nichols图又称为对数幅频率特性图，它以开环频率特性函数的对数幅值为纵坐标，以相角值为横坐标，以角频率为

18、参变量绘制的频率特性图。通常是先画出Bode图，再根据Bode图绘制Nichols图。频域分析法是经典控制理论中重要的系统分析方法。传统的频域分析需要计算数据，绘制Bode图，并求出频域性能指标。这种方法要耗费大量时间、精力，而且绘图、计算精确度不能得到保证。应用MATLAB软件，运行它提供的频率分析函数，就能够方便、简单、快捷的绘制三种曲线，并计算出频域性能指标。1.4 控制系统的根轨迹分析根轨迹法是一种表示特征根与某一参数的全部数值关系的图解方法。与该参数的某一特定数值相应的根，可在根轨迹图上找到。上述可变参数可以是开环传递函数中任何可变参数，但通常取开环增益作为可变参数。即当系统某个可变

19、参数有0时，系统闭环极点在s平面上所描绘出来的轨迹，根轨迹能直观的反映系统参数变化对根分布位置的影响。系统闭环特征方程的根满足1+G(s)H(s)=0则系统的根轨迹方程为：（1-2）在MATLAB中，系统专门提供了绘制根轨迹有关的函数，这些函数能简单快捷的绘制根轨迹或者进行相关的计算。1.5MATLAB软件认识MATLAB 的特点(1)超强的数值运算功能。在MATLAB里，有超过500种的数学、统计、科学及工程方面的函数可供使用，而且使用简单快捷。(2)语法限制不严格，程序设计自由度大。例如，在MATLAB里，用户不用对矩阵预定义就可以使用。(3)程序的可移植性好，基本上不做修改就可以在各种型

20、号的计算机和操作系统上运行。(4)强大的数据可视化功能。在FORTRAN和C语言里，绘图都很不容易，但在MATLAB里，数据可视化非常简单。MATLAB还有较强的编辑图形界面的能力。(5)丰富的工具箱。由各学科领域水平很高的专家编写的功能强劲的工具箱，用户无需编写自己学科范围内的基础程序，而直接进行高、精、尖的研究。MATLAB在控制系统分析中的应用利用MATLAB可对控制系统进行如下有效分析：(1)分析系统的稳定性；(2)求取系统在典型和任意输入信号作用下的时域响应；(3)绘制系统的根轨迹，在根轨迹上可确定任意点的根轨迹增益K值，从而得到系统稳定的根轨迹增益K值范围；(4)利用MATLAB绘

21、制系统的Bode图、Nichols图和Nyquist图等，并求取系统的幅值裕量和相位裕量；(5)分析具有时间延迟系统的频率特性；(6)求取频率响应数据，且根据频率响应数据辩识系统的模型参数；(7)分析系统的能控性和能观测性，并能对不完全能控或不完全能观测的系统进行结构分解。这里重点介绍一下利用MATLAB进行系统根轨迹分析和频域分析的方法。根轨迹绘制(1)绘制系统的零极点图pzmap( )函数可绘制系统的零极点图，其调用格式为：p，z=pzmap(A，B，C，D)或p，z=pzmap(p，z)；p，z=pzmap(num，den)其中列向量P为系统极点位置，列向量z为系统的零点位置。可通过pz

22、map(p，z)绘制出零极点图，图中的极点用 “×”表示，零点用 “o”表示。(2)根轨迹绘制的命令闭环系统的特征根在s平面变化的轨迹，该函数的调用格式为：r，K=rlocus(num，den)；r，K=rlocus(num，den，K)其中返回值r为系统的闭环极点，K为相应的增益。在系统分析过程中，常常希望确定根轨迹上某一点处的增益值K，这时可利用MATLAB中的rlocfind()函数，在使用此函数前要首先得到系统的根轨迹，然后再执行如下命令：K，poles=rlocfind(num，den)或K，poles=rlocfind(A，B，C，D)执行上述命令后，将在屏幕上的图形中生

23、成一个十字光标，使用鼠标移动它至所希望的位置，然后敲击鼠标左键即可得到该极点的位置坐标值poles以及它所对应的增益K值。控制系统的频域分析MATLAB的控制系统工具箱提供了多种求取线性系统频率响应曲线的函数。（1）频率向量 频率向量可由logspace()函数来构成。此函数的调用格式为：=logspace(m，n，npts)（2）系统的Bode图连续系统的伯德图可利用bode()函数来绘制，其调用格式为：mag，phase，=bode(num，den)mag，phase，=bode(num，den，)mag，phase，=bode(A，B，C，D)mag，phase，=bode(A，B，C，

24、D，iu)mag，phase，=bode(A，B，C，D，iu，)式中num，den和A，B，C，D分别为系统的开环传递函数和状态方程的参数，而为频率点构成的向量。有了这些数据就可以使用下面的MATLAB命令：subplot(2，1，1)；semilogx(w，20*log10(mag)subplot(2，1，2)；semilogx(w，phase) 此时便在一个窗口上同时显示出系统的Bode图，其中前一条命令中对幅值向量mag求取分贝(dB)值。如果只想绘制出系统的Bode图，而对获得幅值和相位的具体数值并不感兴趣，则可以采用简单的调用格式：bode(num，den，)；bode(A，B，C

25、，D，iu，）或bode(num，den)；bode(A，B，C，D，iu)（3）幅值裕量和相位裕量在判断系统稳定性时，常常需要求出系统的幅值裕量和相位裕量。利用MATLAB控制系统工具箱提供的margin()函数可以求出系统的幅值裕量与相位裕量，该函数的调用格式为：Gm，Pm，Wcg，Wcp=margin(num，den) 或Gm，Pm，Wcg，Wcp=margin(A，B，C，D)式中Gm和Pm分别为系统的幅值裕量和相位裕量，而Wcg和Wcp分别为幅值裕量和相位裕量处相应的频率值。除了根据系统模型直接求取幅值和相位裕量之外，MATLAB的控制系统工具箱中还提供了由幅值和相位相应数据来求取裕

26、量的方法，这时函数的调用格式为：Gm，Pm，Wcg，Wcp=margin(mag，phase，)式中频率响应可以是由bode()函数获得的幅值和相位向量，也可以是系统的实测幅值与相位向量，为相应的频率点向量。（4）频率响应值MATLAB也提供了直接的求取频率响应数据的函数freqresp()，其调用格式为：F=freqresp(num，den，sqrt(-1)*) 或F=freqresp(A，B，C，D，iu，sqrt(-1)*)式中F为频率响应，为给定的频率范围向量。（5）系统的奈奎斯特图（Nyquist图）nyquist()函数的调用格式为：Re，Im，=nyquist(num，den)R

27、e，Im，=nyquist(num，den，)Re，Im，=nyquist(A，B，C，D)Re，Im，=nyquist(A，B，C，D，iu)Re，Im，=nyquist(A，B，C，D，iu，)其中返回值Re，Im和分别为频率特性的实部向量、虚部向量和对应的频率向量，有了这些值就可利用命令plot(Re，Im)来直接绘出系统的奈奎斯特图。当然也可使用下面的简单命令来直接绘出系统的奈奎斯特图：nyquist(num，den，)或nyquist(A，B，C，D，iu)（6）尼柯尔斯图(Nichols图)函数调用格式为：mag，phase，=nichols(num，den，)或mag，phase

28、，=nichols(A，B，C，D，iu，)可见该函数的调用格式以及返回的值与bode()函数完全一致，事实上虽然它们使用的算法不同， 但这两个函数得出的结果还是基本一致的。但Nichols图的绘制方式和Bode图是不同的，它可由以下命令绘制plot(phase，20*log10(mag)1.6MATLAB环境下的Simulink简介MATLAB环境下的Simulink是当前众多仿真软件中功能最强大、最优秀、最容易使用的一个系统建模、仿真和分析的动态仿真集成环境工具箱。Simulink为用户提供了一个图形化的用户界面（GUI）。它使用方框图表示一个系统，通过图形界面，利用鼠标点击和拖拉方式建立

29、系统模型。它与用微分方程和差分方程建模的传统仿真软件包相比，具有更直观、更方便、更灵活的优点。它不但实现了可视化的动态仿真，也实现了与MATLAB、C或者FORTRAN甚至和硬件之间的数据传递，大大扩展了它的功能。Simulink的主要功能如下：(1)交互建模Simulink提供了大量的功能块以方便用户快速的建立动态系统模型。建模时只需要使用鼠标拖放库中的功能模块并将它们接连起来。(2)交互仿真Simulink提供了交互性很强的非线性仿真环境，可以通过下拉菜单执行仿真，或使用命令进行批处理。仿真结果可以在运行的同时通过示波器或图形窗口查看。(3)扩充和制定Simulink的开放结构允许扩展仿真

30、环境的功能：用户可以用MATLAB、Fortran和C代码生成自定义模块图，并拥有自己的图标和界面；用户还可以将原来的Fortran和C语言代码连接起来。(4)工具箱集成由于Simulink可以直接利用MATLAB的数学图形和编程功能，用户可以直接在Simulink下完成诸如数据分析、过程自动化、优化参数等工作。工具箱提供的高级设计和分析能力可以通过Simulink的屏蔽手段在仿真过程中执行。第二章火电厂锅炉燃烧过程控制系统设计 锅炉燃烧过程是一个将燃料的化学能转化成热能，以蒸汽形式向负荷设备提供热能的能量转换过程。煤由皮带输送到锅炉车间的煤斗，进入磨煤机磨成煤粉，然后通过粉仓与阀门进入炉内，

31、同时由鼓风机送入空气，使煤粉充分燃烧，将煤的化学能转换成热能，烟气经除尘器清除灰分后，由引风机抽出，经烟囱排入大气。以中储式燃煤锅炉燃烧控制系统为例，工艺流程图如图2.1所示。 图2.1燃煤锅炉燃烧控制系统的工艺流程图燃烧过程控制系统有三大任务： 第一个任务是维持汽压恒定。汽压变化表示锅炉蒸汽量和负荷的耗汽量不相适应，必须相应地改变燃料量，以改变锅炉的蒸汽量； 第二个任务是保证燃烧过程的经济性。当燃料量改变时，必须相应地调节送风量，使它与燃料量相配合，保证燃烧过程有较高的经济性； 第三个任务是调节引风量与送风量相配合，以保证炉膛压力不变。 因此，蒸汽锅炉的燃烧过程主要由三个子系统构成：蒸汽压力

32、控制系统、燃料空气比值控制系统和炉膛负压控制系统。图2.2是燃烧过程控制系统示意图，图2.3是原理方框图。图2.2 燃烧过程控制系统示意图图2.3 原理方框图2.1 蒸汽压力控制系统和燃料空气比值控制系统锅炉燃烧的目的是为后续的生产环节提供稳定的压力。一般生产过程中蒸汽的控制是通过调节压力实现的，随着后续环节的蒸汽用量不同，会造成蒸汽压力的波动，蒸汽压力波动较大时，燃料用量和送风量波动也会很大,这样就会给后续的生产造成不良的影响，因此，维持蒸汽压力恒定是保证生产正常进行的首要条件。保证蒸汽压力恒定的主要手段是随着蒸汽压力波动及时调节燃料产生的热量，而燃烧产生热量的调节是通过控制所供应的燃料量以

33、及适当比例的助燃空气实现的。因此，各个控制环节的关系如下：蒸汽压力是最终被控量，根据蒸汽的生成量确定；燃料量根据蒸汽压力确定；送风量根据空气量与燃料量的比值确定。控制量示意图如图2.4所示。图2.5为蒸汽压力控制系统和燃料空气比值控制系统流程图。图2.6为蒸汽压力控制系统和燃料空气比值控制系统框图。 图2.4 控制量示意图图2.5 蒸汽压力控制系统和燃料空气比值控制系统流程图图2.6 蒸汽压力控制系统和燃料空气比值控制系统框图2.2 炉膛负压控制系统 炉膛负压是指炉膛顶部的烟气压力。炉膛负压是反映燃烧工况稳定与否的重要参数，是运行中需要监视和控制的重要参数之一。炉内燃烧工况一旦发生变化，炉膛负

34、压随即发生相应变化。当锅炉的燃烧系统发生故障或异常时，最先在炉膛负压上反映出来，而后才是火焰的变化，其次才是蒸汽参数的变化。因此，监视和控制炉膛负压对于保证炉内燃烧工况的稳定、分析炉内燃烧工况、烟道运行工况、分析某些事故的原因均有极其重要的意义。大多数锅炉采用平衡通风方式，使炉内烟气压力低于外界大气压力，即炉内烟气负压，炉膛内烟气压力最高的部位是炉膛顶部。当炉膛负压过大时，引风量增大，不完全燃烧损失、排烟热损失均增大。甚至使燃烧不稳定，炉膛压力变为正压，火焰及飞灰从炉膛不严处冒出，恶化工作环境，危及人身及设备安全，故应保持炉膛负压在正常范围内。保证炉膛负压的措施是引风量和送风量的平衡。如果负压

35、波动不大，调节引风量即可以实现负压控制；当蒸汽压力波动较大时，燃料用量和送风量波动也会很大，此时，经常采用的控制方法为动态前馈-反馈控制，如图2.7所示。前馈控制的基本概念是测取进入控制过程的干扰信号，在炉膛负压控制系统中，当蒸汽压力波动较大时，燃料用量和送风量的波动会较大，所以通过测取送风量，就可以得到干扰信号，利用反应较快的动态前馈控制，就可以很好的减小干扰信号对系统的影响。将前馈与反馈有效的结合，运用前馈控制在扰动发生后，及时抑制由主要扰动引起的被控量所产生的偏差，同时运用反馈控制消除其余的扰动对负压的影响。图2.7为炉膛负压控制系统流程图，图2.8为炉膛负压控制系统框图。前馈控制系统主

36、要特点如下：（1）属于开环控制 只要系统中各环节是稳定的，则控制系统必然稳定。但若系统中有一个环节不稳定，或局部不稳定，系统就不稳定。另外，系统的控制精度取决于构成控制系统的每一部分的精度，所以对系统各环节精度要求较高。（2）很强的补偿局限性 前馈控制实际是利用同一干扰源经过干扰通道和前馈通道对系统的作用的叠加来消除干扰的影响。因此，固定的前馈控制只对相应的干扰源起作用，而对其他干扰没有影响。而且，在工程实际中，影响生产过程的原因多种多样，系统随时间、工作状态、环境等情况的变化，也会发生变化甚至表现出非线性，这些都导致不可能精确确定某一干扰对系统影响的程度或数学描述关系式。因此，前馈控制即使对

37、单一干扰也难以完全补偿。（3）前馈控制反应迅速 在前馈控制系统中，信息流只向前运行，没有反馈问题，因此相应提高了系统反应的速度。当扰动发生后，前馈控制器及时动作，对抑制被控制量由于扰动引起的动静态偏差比较有效。这非常有利于大迟滞系统的控制。（4）只能用于可测的干扰 对不可测干扰，由于无法构造前馈控制器而不能使用。 该方案中以负压为控制目标，用引风量做成控制闭环，利用前馈控制消除送风量变动对负压的影响。图2.7 炉膛负压控制系统流程图图2.8 炉膛负压控制系统框图第三章火电厂锅炉燃烧过程控制系统数学模型 本文通过查找资料，得到中储式燃煤锅炉的系统辨识函数，并以此为数学模型搭建仿真平台。3.1 蒸

38、汽压力控制系统和燃料空气比值控制系统数学模型燃烧流量被控过程的数学模型为：（3-1）燃料流量与蒸汽压力关系为： （3-2）蒸汽压力与燃料流量关系为: （3-3）蒸汽压力检测变换系统数学模型为： （3-4）燃料流量检测变换系统数学模型为： （3-5）燃料流量与空气流量比值为： （3-6）空气流量被控过程的数学模型为:（3-7）3.2 炉膛负压控制系统数学模型引风量与炉膛负压关系： （3-8）送风量对炉膛负压的干扰为：（3-9）第四章 系统稳定性分析4.1 燃料控制系统燃料流量被控过程的数学模型为：（4-1）MATLAB程序：G_Fuel =tf(2,13 1,'inputdelay

39、9;,3);margin(G_Fuel);Gm,Pm,Wg,Wp=margin(G_Fuel)图4.1 燃料被控对象Bode图经软件仿真结果如图4.1所示，由图可知，无调节器时，燃料控制系统是稳定的。由程序运行结果可知其幅值稳定裕量为11.4dB,相位稳定裕量为97.1dB，对应增益为3.7283。4.2 空气流量控制系统空气流量被控过程的数学模型为： （4-2）MATLAB程序：G_Air =tf(3,11 1,'inputdelay',2);margin(G_Air);Gm,Pm,Wg,Wp=margin(G_Air)图4.2 空气流量被控对象Bode图经软件仿真结果如图4

40、.2所示，由图可知，无调节器时，空气流量控制系统是稳定的。由程序运行结果可知其幅值稳定裕量为9.81dB,相位稳定裕量为80deg，对应增益为3.0956。4.3 炉膛负压控制系统引风量与炉膛负压关系： （4-3）MATLAB程序：G_LeAir_NegPres =tf(10,7 1,'inputdelay',1);margin(G_LeAir_NegPres);Gm,Pm,Wg,Wp=margin(G_LeAir_NegPres)图4.3 引风量与炉膛负压关系的Bode图经软件仿真结果如图4.3所示，由图可知，无调节器时，炉膛负压控制系统是稳定的。由程序运行结果可知幅值稳定裕

41、量为1.32dB,相位稳定裕量为14.3deg，对应增益为1.1641。第五章 火电厂锅炉燃烧过程控制各子系统仿真5.1 燃料控制系统为使系统无静差，燃料流量调节器采用PI形式，即：（5-1）其中，参数Kp和Ki采用稳定边界法整定。先让Ki=0，调整Kp使系统等幅振荡，即系统临界稳定状态。稳定边界法整定的计算方法如下：调节规律KpKiPK/2PIK/2.2Kp/(0.85T) 系统临界振荡仿真框图如图5.1所示，其振荡响应如图5.2所示：图5.1 系统临界振荡仿真框图图5.2 系统临界振荡响应曲线记录此时的振荡周期T=11s和比例参数K=3.8,则Kp=K/2.2=1.73，Ki=Kp/(0.

42、85T)=0.18。在Kp=1.73，Ki=0.18的基础上，对PI参数进一步整定，燃料流量控制系统的仿真框图如5.4所示，其中PI模块的结构如图5.3所示。调节Kp=1.1，Ki=0.1，系统响应如图5.5所示，可见系统有约10%的超调量。图5.3 PI模块结构图5.4 燃料流量控制系统的仿真框图图5.5 燃料流量控制系统单位阶跃响应曲线5.2 蒸汽压力控制系统在燃料流量控制系统整定的基础上，采用试误法整定压力控制系统参数。系统整定仿真框图如图5.6所示。当Ki=0，Kp=1时(此时相当于无调节器，因此系统最简单)，仿真结果如图所示，图5.6为系统仿真框图，图5.7为阶跃输出。图5.6 蒸汽

43、压力控制系统仿真框图图5.7 蒸汽压力控制系统单位阶跃响应曲线由仿真结果可以看出，系统响应超调量约为50%。此时系统调节器最简单，工程上系统响应速度和稳定程度一般。5.3 空气流量控制系统空气流量控制系统的整定方法和燃料流量控制参数整定方法类似，参数Kp和Ki也采用稳定边界法整定。先让Ki=0，调整Kp使系统等幅振荡，即系统临界稳定状态。系统临界振荡仿真框图如图5.8所示，其振荡响应如图5.9所示：图5.8 系统临界振荡仿真框图图5.9 系统临界振荡响应曲线记录此时的振荡周期T=7s和比例参数K=3.3,则Kp=K/2.2=1.5，Ki=Kp/(0.85T)=0.25。在Kp=1.5，Ki=0

44、.25的基础上，对PI参数进一步整定，空气流量控制系统单位阶跃输入的仿真框图如图5.10所示。调节Kp=1.0，Ki=0.1，系统响应如图5.11所示，可见系统有约10%的超调量。图5.10 空气流量控制系统仿真框图图5.11 空气流量控制系统单位阶跃响应曲线5.4 炉膛负压控制系统由系统概述可知，利用动态前馈与反馈控制能使负压控制系统达到比较好的效果，反馈控制在燃料控制系统和空气流量控制系统中已经详细的论述与应用过了，在此就前馈控制做一下讨论。图5.12 前馈控制原理图在图5.12中实线部分为未引入控制措施前的控制框图，虚线部分为所引入的控制部分。Gf(s)G3(s)为干扰源至系统输出的干扰

45、通道传递函数；Gd(s)为前馈调节器函数；G2(s)G3(s)为干扰源至系统输出的控制通道传递函数；G1(s)为给定环节传递函数。系统输出为： （5-2） 由上式可知，消除可测干扰M的影响，可令M的系数为0，即可得到: （5-3）因此，可以得到炉膛负压控制系统前馈补偿函数为： （5-4）图5.13为炉膛负压前馈控制系统参数整定仿真框图，图5.14为未加前馈控制系统的单位阶跃响应曲线，图5.15为加入了前馈控制系统的单位阶跃响应曲线。图5.13 炉膛负压前馈控制系统参数整定仿真框图图5.14 未加前馈控制系统的单位阶跃响应曲线图5.15 加入前馈控制系统的单位阶跃响应曲线 两图对比可知，加入的前

46、馈控制起到了很好的作用。第六章火电厂锅炉燃烧过程控制系统Simulink仿真搭建完整的系统Simulink仿真框图，图6.1为不加扰动时系统的仿真框图，图6.3为不加扰动的响应曲线图。但是在实际的工作环境中往往会出现各种各样的扰动信号，为了模拟真正的工业控制环境，在特定的环节，我们添加了幅值为-1至+1的噪声干扰信号，图6.2为仿真框图，图6.4为相应的响应曲线图。图6.1 未加扰动的燃烧控制系统仿真框图图6.2 加入扰动的燃烧控制系统仿真框图炉膛负压蒸汽压力空气流量图6.3 未加扰动的燃烧控制系统响应曲线随机扰动炉膛负压蒸汽压力空气流量图6.4 加入扰动的燃烧控制系统响应曲线引入干扰信号之后

47、，系统较好的抑制了干扰，曲线达到了稳定。结 论论文从控制系统基础理论入手，简要的了解MATLAB在控制系统中应用的基础后，着重分析了火电厂锅炉燃烧过程控制系统各个环节，其次对各个控制环节建立数学模型，通过MATLAB软件对建立的数学模型进行化简，分析，得出不同参数对系统稳定性和动态性能的影响，最终确定系统各个环节的最佳参数，用于指导实际工作中参数的选择，对控制系统进行仿真确实有很大的实际意义。MATLAB语言易学易懂，对使用者不要求有高深的数学、计算机语言知识以及对编程技巧的灵活运用，并且其数值计算功能强大，研究控制系统便利，图形界面丰富，十分适合用于控制系统的仿真研究。在应用过程中，MATLAB的优越性十分突出，鉴于其开放、扩充能力强，随着其版本的升级，它的功能将越来越强，其应用也必将越来越广。在仿真阶段，本文采用的控制方法是应用最广泛