1、武汉理工大学华夏学院课 程 设 计 报 告 书题 目： 工业被控过程建模与控制器参数的工程整定系 名： 信息工程系 专业班级： 自动化 1092 姓 名： 学 号： 指导教师： 2012 年 7 月 5 日课程设计任务书学生姓名： 专业班级： 自动化 1092 指导教师： 工作单位： 信息工程系 题目： 工业被控过程建模与控制器参数的工程整定一初始条件：1. 给出单容储液槽液位机理建模的过程。2. 给出工程建模的方法。3. 给出广义被控对象的单位阶跃响应输出数据。4. 给出经典的工程整定方法。5. 给出使用 MATLAB 软件建模、绘制与拟合曲线、仿真运行和整定调试的方法。二要求完成的主要任务
2、：（一） 设计任务本设计要求完成如下两个大的设计任务，分别是：1.对某一工程对象进行机理建模，应用 MATLAB 软件对给定的工程数据进行工程测试建模。具体要求为： 参考相关资料，参照过程控制系统课程设计指导书给出的单容储液槽液位数学模型的建立方法，应用机理法对某一工业生产过程建立数学模型。要求模型传递函数为： ，写出建模的详细过程，并绘制出示意图。()1sKGseT 按照下表给出的广义被控生产过程的单位阶跃输入下的输出数据，要求应用MATLAB 软件绘制出其响应曲线。t(s) 0 10 20 40 60 80 100 140y 0 0 0.2 0.8 2 3.6 5.4 8.8t(s) 18
3、0 250 300 400 500 600y 11.8 14.4 16.6 18.4 19.2 19.6 写出应用切线法建立数学模型的具体步骤。要求用计算机绘制其切线图形，并在图上做出详细标注。 写出广义对象传递函数。用计算机在同一图形界面下绘制出阶跃响应曲线和拟合后的曲线。2. 对所建立的被控对象（广义）数学模型，应用 MATLAB软件，建立闭环控制系统模型，并进行工程整定的仿真。最终给出仿真结果和结论。具体要求为： 以得到的对象数学模型为广义被控对象，在 MATLAB 中建立闭环控制系统仿真模型。至少用一种方法进行工程整定。要求给出整定过程中的模型及中间得到的图形数据。 写出具体的整定步骤
4、和得到的整定参数。 根据整定参数，进行仿真验证。给出整定后的闭环控制系统仿真模型和阶跃响应曲线。 对研究过程所获得的主要的数据、现象进行定性或定量分析，得出结论。（二） 说明书撰写要求1. 纸张格式：要求统一用 A4 纸打印，页面设置上空 2.5cm，下空 2.0cm，左空2.5cm，右空 2.0cm）：2. 正文层次：正文内容层次序号为：1、1.1、1.1.1，其中正文标题；一级标题 1.（黑体小 2 号加粗） ，二级标题 1.1（黑体小三号） ，三极标题 1.1.1（黑体小四号） 。正文内容格式：宋体五号，1.25 倍行距。3. 正文内容一般包括：选题背景：说明本课题应解决的主要问题及应达
5、到的技术要求；简述本设计的指导思想。 （设计目的中已有阐述）方案论证：说明设计原理并进行方案选择，阐明为什么要选择这个设计方案以及所采用方案的特点。设计内容：对设计工作的详细表述。要求层次分明、表达确切。结果分析：对研究过程总所获得的主要的数据、现象进行定性或定量分析，得出结论和推论。结论或总结：对整个研究工作进行归纳和综合。参考文献：不少于 5 个，并应按文献号、作者、文献题名、出版地：出版社和出版年等顺序书写。如：1 戴军，袁惠新.膜技术在含油废水处理中的应用.膜科学与技术，2002.4. 图纸（或其它）要求图纸要求：图面整洁，布局合理，线条粗细均匀，圆弧连接光滑，尺寸标注规范，要求使用计
6、算机绘图。曲线图表要求：所有曲线、图表、线路图、流程图、程序框图、示意图等不准徒手画，必须按国家规定标准或工程要求绘制（采用计算机辅助绘图） 。课程设计说明书（报告）中图表、公式一律采用采用阿拉伯数字连续编号。图序及图名置于图的下方；表序及表名置于表的上方；说明书（报告）中的公式编号，用括号括起来写在右边行末，其间不加虚线。三时间安排：时 间 设计内容6 月 25-26 日 查阅资料，用机理法为某一工业被控过程建立数 学模型6 月 27 日 查阅 MATLAB 资料，绘制出阶跃响应曲线6 月 28 日 用切线法建立被控过程数学模型6 月 29-30 日 根据所得数学模型，用 MATLAB 建立
7、闭环系统仿真模型7 月 1-3 日 进行闭环系统的工程整定，得出整定结果7 月 4 日 完成设计报告7 月 5-6 日 答辩指导老师签字： 年 月 日目录1 选题背景和设计任务 11.1 选题背景 .12.2 MATLAB 二维绘图简介 23 设计方案 .2.1 建模的一般方法 323.2 控制体统参数的工程整定 .33.3 最佳方案 .54 单容储液槽液位控制系统建模 .64.1 被控对象的解析 .64.2 单容储液槽的建模 .64.3 实验法建模 .75 控制系统的工业整定 85.1 广义对象的阶跃响应曲线 .85.2 数学模型的参数确定 .95.3 拟合后的数学模型与曲线 .105.3.
8、1 拟合模型 105.3.2 拟合曲线 .105.4 控制器参数的工程整定 .115.5 控制器的选择 .136 总结 .137 参考文献 .1311 选题背景和设计任务1.1 选题背景过程控制是自动技术的重要应用领域，它是指对液位、温度、流量等过程变量进行控制，在冶金、机械、化工、电力等方面得到了广泛应用。尤其是液位控制技术在现实生活、生产中发挥了重要作用。锅炉汽包液位的控制，如果锅炉内液位过低，会使锅炉过热，可能发生事故。在这些生产领域里，基本上都是劳动强度大或者操作有一定危险性的工作性质，极容易出现操作失误，引起事故，造成厂家的的损失。可见，在实际生产中，液位控制的准确程度和控制效果直接
9、影响到工厂的生产成本、经济效益甚至设备的安全系数。所以，为了保证安全条件、方便操作，就必须研究开发先进的液位控制方法和策略。由于工业生产过程复杂多样，因此，在设计工业生产过程控制系统时，首先必须花大量的时间和精力去了解该工业生产过程的基本原理、操作过程和过程特性，这是设计和实现一个工业生产过程的首要条件。要实现过程自动控制，还需要对整个工业生产过程的物料流、能源流和生产过程中的有关状态进行准确的测量和计量。根据测量得到的数据和信息，用生产过程工艺和控制理论知识来管理和控制该生产过程。MBATLA 软件是一款进行高效工程计算、数据分析与可视化编程、系统仿真、科学和工程绘图等功能强大的优秀软件。能
10、够用于系统建模和仿真，方便用于系统参数整定。1.2 设计任务工业生产过程的扰动作用使得生产过程操作不稳定，从而影响工厂生产过程的经济效益。过程控制的任务就是使生产过程达到安全、平稳、优质、高效。作为自动化最根本的目标应是使生产过程安全并平稳的运行。本课程设计是工业被控过程建模与控制器参数的工程整定，主要是针对单容储液槽的水位平衡进行机理建模。本设计要求完成如下两个大的设计任务，分别是：1）对某一工程对象进行机理建模，应用 MATLAB 软件对给定的工程数据进行工程测试建模。2）对所建立的被控对象（广义）数学模型，应用 MATLAB 软件，建立闭环控制系统模型，并进行工程整定的仿真。最终给出仿真
11、结果和结论。2 MATLAB 仿真软件的介绍2.1 MATLAB 应用简介MBATLA 是“矩阵实验室(Matrix Laboratory)”的缩写，它是一种以矩阵运算为基础的交互式程2序语言，是专门针对科学和工程中计算和绘图的需求而开发的。MBATLA 软件是一款进行高效工程计算、数据分析与可视化编程、系统仿真、科学和工程绘图等功能强大的优秀软件。其中 simulink 仿真软件包功能强大、操作简单，是进行系统建模、仿真、科学与工程分析计算的强大有力工具。它的基本数据元素是没有维数限制的阵列（矩阵） 。MATLAB 是一种高效的工程计算语言，其典型应用主要包括以下几个方面：数学计算； 算法开
12、发； 数据采集； 系统建模和仿真； 数据分析和可视化； 科学和工程绘图； 应用软件开定义变量并给变量赋值命令可直接在命令窗口中输入：变量名=元素值。例如定义变量 A 为矩阵可直接在命令窗口中输入 A=1,2,3,4;5,6,7,8;9,10,11,12。注意：同一行内不同列元素用空格或逗号分隔开，不同行之间用分号分隔。又如可直接在命令窗口中输入一个行向量 A=1,2,3,4,5。2.2 MATLAB 二维绘图简介使用 plot 命令绘制图形 ()plotx,y其中 x,y 可为向量或矩阵 ,当 x,y 同为 n 维向量时，使用该命令，将绘制出以向量 x 的元素为横坐标，y 向量对应元素为纵坐标
13、的图形。3 设计方案.1 建模的一般方法 3 建立被控过程数学模型的方法一般有：机理建模、试验建模、混合建模。机理建模是根据对象或是生产过程的内部机理，写出各种有关的平衡方程，如物料平衡方程、能量平衡方程、动量平衡方程以及某些设备方程、化学反应定律等，从而得到对象的数学模型。这类模型通常称为机理模型。这种方法建立的模型的最大优点是具有非常明显的物理意义，模型具有很大的1234567890A3适应性，便于模型参数的调整。试验建模是在机理建模难以建立的情况下，可以采用试验建模的方法得到对象的数学模型。试验建模针对所要研究的对象，人为地施加一个输入作用，然后用仪表记录表征对象特性的物理量随着时间变化
14、的规律，得到一系列的试验数据或者是曲线。通过对曲线的分析获得必要的规律信息。混合建模将机理建模和试验建模结合起来就是混合建模。混合建模是一种比较实用的方法，它先由机理分析的方法提出数学结构模式，然后对其中某些未知的或不确定的参数利用试验的方法予以确定。3.2 控制体统参数的工程整定方法一：经验法：根据经验和先验知识确定一组参数，然后根据各参数的影响，调整参数，直至满意为止。由于我们经验和先验知识不足，所以本次设计中不使用此方法。方法二：临界比例度（带）法：比例度（带） ：与比例系数 成反比关系。这种整定方(%)cK法是在闭环情况下进行的。设 ， ，使控制器工作在纯比例，使系统的输出响应情况下，
15、TI0D将比例带由大逐渐变小(对应的比例系数 由小逐渐变大 )呈现等幅振荡。等幅震荡的波形及相关参数cK的获取如图 3.1 所示：图 3.1 等幅振荡波形图临界比例度法整定经验公式如表 1 所示，最后对参数进行微调，直到动态过程满意为止。表 1 临界比例度法整定经验公式 Ti(S) Td(S) P 2 S PI 2.2 S TS/1.2 PID 1.6 S 0.5TS 0.125TS 方法三：衰减曲线法：在闭环系统中，先把控制器设置为纯比例作用，然后把 比例带由大逐渐减小(对应的比例系数 由小逐渐变大 ) ，加阶跃扰动观察输出响应的衰减过程，直至出现 4:1 衰减过cK4程为止，如图 3.2
16、所示。图 3.2 加阶跃扰动的输出曲线这时的比例带称为 4:1 衰减比例带，用 表示之。相邻两波峰间的距离称为衰减周期 ,记录s sT和 。按表 2 所示经验公式整定。最后对参数进行微调，直到动态过程满意为止。sT表 2 衰减曲线法相关参数的整定（%） (S) lT(S) dTP sPI 1.2 s0.5 sTPID 0.8 s0.3 s0.1 sT方法四：响应曲线法：其控制原理图如图 3.3 所示。图 3.3 响应曲线法控制原理图令控制器的输出为幅度为 ,的阶跃信号，则对象经测量变送器后的输出 ， 由该图可确定 、0x )(tyT 和 。相关参数的确定如图 3.4 所示。cK5图 3.4 输
17、出波形及相关参数的确定0)(xyKc通过下式将比例系数转化为比例度：利用表 3 所示的经验公式，就可计算出对应于衰减率为 4:1 时控制器的相关参数。表 3 响应曲线法相关参数整定的经验公式3.3 最佳方案单容储液槽液位控制系统主要是实现对液位的控制，单容储液槽有两个可以控制的变量，一个进水口的流量，另一个是出水口的流量。通过容积与流量的平衡方程很容易建立数学模型。所以被控制过程建模采用机理建模。参数整定可以选用的方法也较多在本设计中采用临界比例度法进行整定。跟其它的方法相比较，临界比例度法简单快捷，所以可以采用临界比例度法进行整定。1%c64 单容储液槽液位控制系统建模4.1 被控对象的解析
18、本设计探讨的是单容储液槽的液位控制问题。为了能更好的选取控制方法和参数，有必要知道被控对象上储液槽的结构和特性。由图 4.1 所示可以知道，单容储液槽的流量特性：储液槽的出水量与水压有关，而水压又与水位高度近乎成正比。这样，当储液槽水位升高时，其出水量也在不断增大。所以，若阀 开度适当，在2V不溢出的情况下，当储液槽的进水量恒定不变时，水位的上升速度将逐渐变慢，最终达到平衡。由此可见，单容储液槽系统是一个自衡系统。图 4.1 单容储液槽结构图4.2 单容储液槽的建模这里研究的被控对象只有一个，那就是单容储液槽（图 4.1） 。要对该对象进行较好的计算机控制，有必要建立被控对象的数学模型。正如前
19、面提到的，单容储液槽是一个自衡系统。根据它的这一特性，我们可以用阶跃响应测试法进行建模。如图 4.1，设储液槽的进水量为 ，出水量为 ，储液槽的液面高度为 h，出水阀 固定于某1Q2 2V一开度值。若 作为被控对象的输入变量，h 为其输出变量，则该被控对象的数学模型就是 h 与 1Q 1Q之间的数学表达式。根据动态物料平衡关系有（4-1）12dhQCt将式（2-1）表示为增量形式（4-2）12dt式中， 、 、 分别为偏离某一平衡状态 、 、 的增量； C储液槽底面积。1Q2h102Q0h7在静态时， ； ；当 发生变化时，液位 h 随之变化，阀 处的静压也随之变21Q0/dth1Q2V化，
20、也必然发生变化。由流体力学可知，流体在紊流情况下，液位 h 与流量之间为非线性关系。但2为简化起见，经线性化处理，则可近似认为 与 成正比，而与阀 的阻力 成反比，即122R或 （4-3）2hQR2hRQ式中， 为阀 的阻力，称为液阻。2V将式（2-3）代入式（2-2）可得（4-4）221dhRCRt在零初始条件下，对上式求拉氏变换，得：（4-5）1)(210 sSTKsQHG式中, 为水箱的时间常数（注意：阀 的开度大小会影响到水箱的时间常数） ， 为过sCRT22V2RK程的放大倍数。令输入流量 ， 为常量，则输出液位的高度为：sR/)(010(4-6)TKTsH/1)()即 （4-7）0
21、)tThtRe当 时， 因而有t)(（4-8）0K输 出 稳 态 值阶 跃 输 入当 t=T 时，则有（4-9）100()().632.()hTReKRh4.3 实验法建模实验法建模是工具工业过程中输入、输出的实测数据进行的某种数学处理后得到数学模型的建模方法。根据给被控过程施加的激励信号和数据分析方法的不同，主要有如下方法：a) 时域方法b) 频域方法8c) 统计相关法本次设计中使用的是时域发中的阶跃响应曲线法：对处于开环、稳态被控过程，使其输入作阶跃变化，测得被控过程的阶跃响应曲线，求出被控过程输入与输出之间的动态数学关系传递函数。要确定被控对象的传递函数，可以分两步来做：1 确定传递函数
22、的形式一般确定传递函数的形式要根据被控过程的先验知识、简历数学模型的目的及对模型的准确性要求拉起选择。在满足精度的情况下，尽量选择低阶形式。本次设计中使用一阶模型即可。2 确定模型参数根据阶跃响应曲线，求出模型参数。对于一阶模型，此次用切线法来确定。实际生产过程的阶跃响应曲线呈现如图 4.2 所示的 S 型曲线是最常见的。用切线法确定一阶惯性加纯滞后的特征参数方法如下。图 4.2 切线法确定一阶惯性加纯滞后的特征参数在此曲线的拐点 D 处做一切线，它与时间轴交于 B 点，与响应稳态值渐近线交于 A 点。图 4.2 中的 OB 即为对象的滞后时间 ，BC 为对象的时间常数 T。5 控制系统的工业
23、整定5.1 广义对象的阶跃响应曲线控制系统的参数整定可以用 MATLAB 仿真之后，进行数据分析实现。按照表 4 给出的广义被控生产过程的单位阶跃输入下的输出数据，要求应用 MATLAB 软件绘制出其响应曲线。表 4 单位阶跃输入下的输出数据t(s) 0 10 20 40 60 80 100 140y 0 0 0.2 0.8 2 3.6 5.4 8.89t(s) 180 250 300 400 500 600y 11.8 14.4 16.6 18.4 19.2 19.6由 MATLAB 绘出的图像如图 5.1 所示图 5.1 单位阶跃输入下的输出响应曲线由式 ，结合图形可得：0)(xyKc 6
24、.190.cK5.2 数学模型的参数确定由上分析数学模型的参数的确定用切线法，由 MATLAN 仿真等到的切线如图 5.2 所示。图 5.2 响应曲线的切线图由图中所标参数值可以近似得到 ,T=270s-37s=233s，又有前所得广义对象的传递函数如下：s3710()1sKGseT将 T 的数值代入式中可以得到广义对象的数学传递函数为：19.6cK ses37126.9)(5.3 拟合后的数学模型与曲线5.3.1 拟合模型根据上面计算出的广义对象的数学传递函数，使用 matlab 画出模型如图 5.3 所示：图 5.3 拟合系统仿真图5.3.2 拟合曲线由拟合系统仿真图得到拟合曲线图，如图
25、5.4 所示：11图 5.4 拟合曲线图5.4 控制器参数的工程整定运行 MATLAB，设计系统结构，通过 MATLAB 的 simulink 模块画出仿真结构图，然后设定仿真时间为 600s，积分微分环节设置为零，输入传函和延迟模块的参数，之后进行调试仿真，仿真结构图如图5.5 所示。图 5.5 系统仿真图本次设计采用临界比例度（带）法。比例度（带） 与比例系数 成反比关系，这种整定(%)cK方法是在闭环情况下进行的。设 ， ，使控制器工作在纯比例情况下，将比例带由大逐TI0D渐变小(对应的比例系数 由小逐渐变大 )，使系统的输出响应呈现等幅振荡，波形如图 5.6 所示。 cK图 5.6 临
26、界比例度法输出响应等幅振荡曲线图由图 5.5 可得 ，由图 5.6 可得 。则按经验公式可得，当采用5.0cKssTs14293712纯比例时，比例系数 可得出 ，将比例器的放大系数改为 0.275，在进行仿12SK275.0.真得到波形如图 5.7 所示。图 5.7 比例控制输出波形当采用 PI 调节时，由经验公式得 , ，则 ，25.0.K120.4iT083.12iK将调节器的参数设定为此值，所得波形如图 5.8 所示。图 5.8 比例积分控制输出波形当采用 PID 控制器时，由经验公式可得出 ， ，348.0615.K72145.0iT， , ，将控制器的参数设定为此0139.72iK
27、42.015.sdTd值得到 PID 调节的波形如图 5.9 所示。13图 5.9 比例积分微分控制输出波形5.5 控制器的选择通过参数整定得出的仿真图形可以看出：比例控制调节速度较快，过渡时间短，系统较达到稳定。比例积分控制调节速度比比例控制稍慢一点，超调较大，系统达到稳定的时间较长。比例积分微分控制调节速度最快，但系统稳定时间也较长。所以对于单容储液槽的液位控制系统可以采用比例控制。比例控制克服干扰能力强、控制及时、过渡时间短。比例控制结构简单，系统稳定性好，能够满足单容水箱的液位控制系统的控制要求。6 总结通过本课程设计加深了对过程控制和自动控制思想的理解，对典型的 PID 调节有更深刻
28、的认识，同时掌握了工业被控过程的建模和控制系统的工程整定的一般方法，提升了对过程控制的兴趣，为今后的进一步学习打下了良好的基础。本课程设计主要是两大部分，一个是过程的建模，另一个是参数的整定。对于每一个板块都有很多的方法实现，在实际的操作中要根据设计的目的、被控的对象，选择合适的方法。在这个过程当中，对于方案的论证就是一个很重要的环节，它决定了整个课程设计的质量和效率。如果犯了方向性的错误，最后的结果是事倍功半。这次的课程设计让我基本学会了了 MATLAB 软件的使用，以前从未学过这个软件，但这次因为要用到所以必需的自己去找资料学习,在实践过程中通过查找资料、分析资料以及请教同学，这使我明白了
29、学习是一个长期积累的过程，在以后的工作、生活中都应该不断的学习，努力提高自己知识和综合素质。147 参考文献1 施仁 刘文江等编.自动化仪表及过程控制.第三版，北京：电子工业出版社,2003.2 林德杰编.过程控制仪表及控制系统.北京:机械工业出版社2009年.3 高志宏主编.过程控制与自动化仪表.浙江，浙江大学出版社，2006年.4 厉玉鸣主编.化工仪表及自动化.化学工业出版社，2004年5 吴勤勤主编.控制仪表及装置. 第二版,北京：化学工业出版社，2002年.6 陈夕松 汪木兰主编.过程控制系统.北京：科学出版社,2005.课程设计评分表评 分 项 目 评分成绩1选题合理、目的明确（10
30、 分）2设计方案论证严谨、正确（30 分）3设计过程阐述清晰（25 分）4设计报告整理规范（10 分）5答辩（25 分）总 分（100 分）答辩记录：1.如何使用 MATLAB 软件绘制被控过程的阶跃响应曲线？答：使用 命令绘制阶跃响应图形，其中 为时间向量， 为阶跃响应输出数据所构()t,yploty成的向量。2.在建立被控过程（对象）数学模型时，对象的比例系数（静态放大倍数）如何确定？答：由阶跃响应曲线测得输入幅值 ，输出量 y 的起始值 和稳态值 ，按公式：0X(0)()进行计算。0()yKx3. 为了得到控制用的实用的数学模型，在建模过程中需要进行什么样的近似处理？答：为了得到控制用的实用的模型，需要突出主要因素，忽略次要因素。为此往往要做许多近视处理，如线性化，模型降阶处理等。实际应用中，被控过程的传递函数或其他动态数学模型的阶次一般不高于三阶，大量采用具有纯滞后的一阶和二阶模型。最常用的是带纯滞后的一阶形式。指导教师综合评语：指导教师（签名） 日 期： 年 月 日
