过控课设精馏塔温度控制系统设计与仿真.doc

1、 自动化工程与科学学院自动化工程与科学学院过程控制系统与仪表课程设计精馏塔温度控制系统设计与仿真班级：11 级自动化 2班姓名： 董文杰学号：9指导老师： 哀薇日期： 2过程控制系统课程设计精馏塔温度控制系统设计与仿真一、 研究对象FB , xBVQHPDGLLDLBFmD , xDmmLQVVRTTSu1u2Pp图 1 精馏塔温度控制问题某精馏塔的工艺流程如图 1 所示，现要求对精馏段温度 TR 和提馏段温度 TS都进行有效的控制，以确保塔顶和塔底产品的质量。图 1 中，F 为进料量，它受上游流程控制，为精馏塔温度的主要干扰之一，其它干扰包括进料组成与温度变化、塔底蒸汽量变化、塔顶回流冷凝后

2、温度变化等；L 为塔顶冷回流量，拟作为精馏段温度 TR 的控制手段；塔底蒸汽量 QH 拟作为提馏段温度 TS的控制手段。u 1为调节阀 VL 的相对输入信号，u 2为调节阀 VQ的相对输入信号（以 DDZ III 型为例，当输入电流为 4 mA 时，对应相对输入信号为 0 %；当输入电流为 20 mA 时，对应相对输入信号为 100 %） ，P 为精馏塔顶压力，其变化可基本忽略，P p 为泵出口压力， Pp受塔顶产品调节阀 VD 开度的影响，变化范围较大。图 1 中Lm、V m、F m 分别为 L、V、F 的测量值。为便于控制方案研究，假设如下：（1） 该精馏塔的静态工作点为 T0 = 140

3、 ，F 0 = 60 T/hr（吨/小时） ，L 0 = 20 T/hr，V 0 = 15 T/hr，u 10 = 25 %，f L0 = 75 %， u20 =25 %，f Q0 = 25 %，P p0 = 0.9 MPa，P t0 = 0.86MPa。这里， fQ 为调节阀 VQ 相对流通面积，f L 为调节阀 L 相对流通面积。（2） 精馏段和提馏段温度的测量范围都为 0 200 ，进料量 F 的测量范围为 0 100 T/hr，塔顶冷回流量 L 的测量范围为 0 50 T/hr，塔底回流量 V 的测量范围为 0 25 T/hr。L 、V 、F 的测量值：L m、V m、F m 均用%来

4、表示，即Lm、V m、F m 的最小值为 0，最大值为 100。（3） 流量测量仪表的动态滞后忽略不计；而温度测量环节可用带纯滞后的一阶环节来近似，温度测量环节的一阶时间常数 ， ，纯滞后时间 ，7.RT6.s 3.0R，单位为分。4.s（4） 考虑到精馏塔操作的安全性，控制阀 VL 选用气关阀，控制阀 VQ 选用气开阀，3假设控制阀都为线性阀，其动态滞后忽略不计，动态特性可表示为， 。1)()(1sufGLL 1)()(2sufGQ（5） 对于塔顶冷回流对象，假设控制通道与扰动通道的动态特性可表示为：， 。)()(2RTKfRLp 2R2)(dpRdKPL其中 为控制阀 VL 相对流通面积的

5、变化量，% ；T R2 基本不变，这里设sfL分； KR2、K Rd2 在一定范围内变化，这里设 KR2、K Rd2 的变化范围分4.02RT别为 (T/hr)/%； (T/hr)/MPa。.,0,12Rd（6） 对于塔底回流对象，假设控制通道与扰动通道的动态特性可表示为:， 。)()(22sTfQsGSVHSp 2S2)()(ddsVG其中 为控制阀 VR 相对流通面积的变化量， %；T S2 基本不变，这里设fV分；K S2、K Sd2在一定范围内变化，这里设 KS2、K Sd2的变化范围分别5.02ST为 (T/hr)/%； T/hr。.,15,102Sd（7） 对于温度对象，假设控制通

6、道与扰动通道的动态特性可表示为；1221212 sRssSkkeTGe ；sTKFs RdRdRd exp)()(1sGSSSS 其中对象特性参数均可能在以下范围内变化： /(T/hr)，0.4,51/(T/hr) ， /(T/hr)， /(T/hr) ，0.2,512K0.3,21K,32K分， 分， 分，3T.4,12T.,T分， 分， 分，.,41分； /(T/hr)， 分，2 .,5Rd 0.4,1RdT分， /(T/hr)， 分，0.,Rd 061S 53S分。S二、 研究任务对于上述被控过程，假设被控变量 T 所受的主要扰动为进料量 F、泵出口压力 Pp的变化、蒸汽压力变化 Pf，

7、而且变化范围为： T/hr， MPa, 206F02.9pP；另外，被控变量 T 的设定范围为 。试应用单回路、205fPKpa 154,SRT串级、前馈、比值、选择、Smith 预估、解耦等控制方法，设计至少 2 套控制系统，达到控制精馏塔温度的目的。对于每一套控制方案，具体要求：1、 说明所采用的控制方案以及采用该方案的原因，并在工艺流程上表明该控制系统。42、 确定所用控制器的正反作用，画出控制系统完整的方框图（需注明方框图各环节的输入输出信号） ，并选择合适的 PID 控制规律。3、 在 SIMULINK 仿真环境下，对所采用的控制系统进行仿真研究。具体步骤包括：（1） 在对象特性参数

8、的变化范围内，确定各环节对象的传递函数模型，并构造SIMULINK 对象模型；（2） 引入手动/自动切换环节，在手动状态下对控制通道、干扰通道分别进行阶跃响应试验，以获得“广义对象”开环阶跃响应曲线；（3） 依据 PID 参数整定方法，确定各控制器的参数；（4） 在控制系统处于“闭环”状态下，进行温度设定值跟踪响应试验、干扰塔底回流量、P p与 F 对系统输出的扰动响应试验，并获得相应的响应曲线；（5） 在各控制器参数均保持不变的前提下，当对象特性在其变化范围内发生变化时，重新进行温度设定值跟踪试验与扰动响应试验，并获得相应的响应曲线。4、 根据不同控制方案的闭环响应曲线，比较控制性能（包括是

9、否稳定、衰减比、超调量、过渡过程时间等） 。三、 控制方案精馏塔温度控制的开环回路如下图所示：温度控制开环回路由图我们可知，被控变量 TS、T R控制回路都存在两种主要干扰。被控变量 TR主要受到干扰量是进料量 F和泵出口压力 PP；被控变量 TS主要受到干扰量为进料量 F和塔底蒸汽量V。且被控变量 TS、T R相互耦合。为此，我设计了两种方案来实现控制精馏塔的温度。第一种，忽略耦合的影响，在每条控制回路上加上反馈控制，构成单回路反馈控制；第二种，对耦合系统先进行解耦，分别采用单回路反馈控制回路。（1）第一种方案，构造类似单回路反馈电路。虽然被控变量之间有耦合，但是两者之间的耦合强度不强，所以

10、可以不做考虑。从传递函数可以看出，干扰量对于被控变量的影响不是很大，可以采取直接反馈回路消除。51、控制方案的工艺控制流程图图 1-1 控制工艺流程图2、控制系统方框图系统方框图6图 1-2 系统仿真方框图其中，被控变量 TR 测量环节的传递函数模块 为：图 1-3 被控变量 TR 测量环节的传递函数模块(Gmr)被控变量 TS 测量环节的传递函数模块为：图 1-4 被控变量 TS 测量环节的传递函数模块(Gms)控制通道与扰动通道的动态特性传递函数模块：图 1-5 GR 传递函数模块( 包含 G11、G 21)图 1-6 GS 传递函数模块(包含 G12、G 22)进料量 F 的干扰通道的传

11、递函数模块：7图 1-7 进料量 F 的干扰通道的传递函数模块控制器 PID1、PID2 都是正作用，都选择 PID 控制规律。3、确定各环节对象的传递函数模型（a）对于塔顶冷回流对象，控制通道与扰动通道的动态特性传递函数模型为：，14.03)(L)(RsfsGp 2()()150RdpLsGP（b）对于塔底回流对象，假设控制通道与扰动通道的动态特性传递函数模型为:，15.0)()(2sfQsVHSp 2()10dTsV（c）对于温度对象，假设控制通道与扰动通道的动态特性传递函数模型为；3221.04sRssS eGe；FTsRRd 0.3exp1.3)()( ssGSS .20.45（e）温

12、度测量环节的传递函数模型：，0.31().7sRTse 0.4.1()6STGse4、切换到手动状态，得出各开环阶跃响应曲线(所有控制变量及干扰量初始状态都为 0)（a）开始 10min 后，仅给 TR 控制信号+50 的阶跃信号；再过 40min 后，给 TR 控制信号为-50 的阶跃信号（横轴表示时间，纵轴表示幅度）：8TR 阶跃响应曲线 TS 阶跃响应曲线（b）开始 10min 后，仅给 TS 控制信号+50 的阶跃信号；再过 40min 后，给 TS 控制信号为-50 的阶跃信号（横轴表示时间，纵轴表示幅度）：TR 阶跃响应曲线 TS 阶跃响应曲线（c）开始 10min 后，仅给干扰量

13、 F+20 的阶跃信号；再过 40min 后，给干扰量 F 为-20 的阶跃信号（横轴表示时间，纵轴表示幅度）：进料量阶跃作用的 TR 响应曲线 进料量阶跃作用的 TS 阶跃响应曲线（d）开始 10min 后，仅给干扰量 PP+0.02 的阶跃信号；再过 40min 后，给干扰量 PP 为-0.02 的阶跃信号（横轴表示时间，纵轴表示幅度）：9泵出口压力阶跃作用的 TR 响应曲线 泵出口压力阶跃作用的 TS 响应曲线（e）开始 10min 后，仅给干扰量 V+15 的阶跃信号；再过 40min 后，给干扰量 V 为-15的阶跃信号（横轴表示时间，纵轴表示幅度）：塔底回流量阶跃作用的 TR 响应

14、曲线 塔底回流量阶跃作用的 TS 响应曲线由图（a） 、 （b）可知，被控变量 TR、T S 的控制通道是相互耦合的，T R 的控制变量对被控变量 TS 的影响大于 TS 的控制变量对被控变量 TR 的影响；由图（c）可知，干扰量 F 对于对于被控变量 TS 的影响较大点；有图（d）可知，虽然干扰量 PP 在被控变量 TR 控制回路，但因为耦合作用，也影响到被控变量 TS；有图（e ）可知，干扰量 V 只影响被控变量TS，因为干扰量 V 作用在被控变量回路 TS 的末端。5、控制器参数的整定（a）控制 TR 的控制器 PID1 的参数为：P=1.8，I=0.7，D=0.5（b）控制 TS 的控

15、制器 PID2 的参数为：10P=3.2，I=1.2，D=0.56、控制系统处于“闭环”状态下的性能测试（a）温度设定值跟踪响应试验（温度设定值为 140）：温度设定值 TR 响应曲线 温度设定值 TS 响应曲线（b）泵出口压力 PP 对系统输出的扰动响应试验（开始 40min 后，给阶跃信号+0.02MPa）：泵出口压力干扰 TR 响应曲线 泵出口压力干扰 TS 响应曲线（c）进料量 F 对系统输出的扰动响应试验（开始 40min 后，给阶跃信号+20T/hr）：进料量干扰 TR 响应曲线 进料量干扰 TS 响应曲线11（d）塔底回流量 V 对系统输出的扰动响应试验（开始 40min 后，给

16、阶跃信号+15T/hr）：塔底回流量干扰 TR 响应曲线 塔底回流量干扰 TS 响应曲线（f）三种干扰量同时对系统输出的扰动响应试验（开始 40min 后，三种干扰量同时给阶跃信号）：三种干扰同时作用 TR 响应曲线 三种干扰同时作用 TS 响应曲线由以上图形比较可以得出结论：控制系统在“闭环”状态下，温度设定值跟踪响应性能好，响应时间短；干扰量 PP 与 F 可以通过反馈控制迅速消除，保持系统的稳定；干扰量V 是主要干扰量，仅靠反馈控制可以消除，但是超调量大，调节时间长，不利于系统安全。7、变换对象特性传递函数，测定系统“闭环”动态性能，RL()0.241psGf2()()20RdpLsGP

17、，2().()05HSpVQsfs2()15dTsV12；41.534.302sRssS eGe；()2()xp4RRdTsF0.e3()1SSdGss（a）温度设定值跟踪响应试验（温度设定值为 140）：温度设定值 TR 响应曲线 温度设定值 TS 响应曲线（b）泵出口压力 PP 对系统输出的扰动响应试验（开始 40min 后，给阶跃信号+0.02MPa）：泵出口压力干扰 TR 响应曲线 泵出口压力干扰 TS 响应曲线（c）进料量 F 对系统输出的扰动响应试验（开始 40min 后，给阶跃信号+20T/hr）：13进料量干扰 TR 响应曲线 进料量干扰 TS 响应曲线（d）塔底回流量 V 对

18、系统输出的扰动响应试验（开始 40min 后，给阶跃信号+15T/hr）：塔底回流量干扰 TR 响应曲线 塔底回流量干扰 TS 响应曲线通过以上图形，可以看出：无论是温度设定值跟踪响应曲线还是扰动响应曲线，都发生强烈的振荡。所以在各控制器参数均保持不变的前提下，当对象特性在发生变化时，系统将不稳定，即每当对象特性发生变化时，我们都要进行 PID 整定，才能维持系统的稳定。（2）第二种方案，先对控制通道进行解耦，塔顶塔底均采用单回路反馈控制方法。被控变量 TR、T S的控制通道是相互耦合的，当 TS受到大干扰时，由于耦合作用，同样会严重影响被控变量 TR。一般而已，我们对于塔顶精馏段要求比较高，

19、所以增加解耦环节，使被控变量 TR不再受到塔底回流量 V的干扰，有利于被控变量 TR的稳定快速控制。对于塔顶回路，干扰量影响不大，采用单回路反馈控制，即可以达到控制要求。对于塔底回路，由于控制过程由两条通道并联，且滞后时间都不一样，所以无法采用 Smith预估补偿控制，故仍然采用单回路反馈控制。141、控制方案的工艺控制流程图图 2-1 控制工艺流程图2、控制系统方框图系统方框图15图 2-2 系统仿真方框图其中，被控变量 TR 测量环节的传递函数模块为：图 2-3 被控变量 TR 测量环节的传递函数模块(Gmr)被控变量 TS 测量环节的传递函数模块为：图 2-4 被控变量 TS 测量环节的

20、传递函数模块(Gms)控制通道与扰动通道的动态特性传递函数模块：图 2-5 GR 传递函数模块( 包含 G11、G 21)图 2-6 GS 传递函数模块(包含 G12、G 22)进料量 F 的干扰通道的传递函数模块：16图 2-7 进料量 F 的干扰通道的传递函数模块解耦环节的传递函数模块：图 2-8 解耦环节传递函数模块 N(包含 N21、N 12)控制器 PID1、PID2 都是正作用，都选择 PID 控制规律。3、确定各环节对象的传递函数模型（a）对于塔顶冷回流对象，控制通道与扰动通道的动态特性传递函数模型为：，14.03)(L)(RsfsGp 2()()150RdpLsGP（b）对于塔

21、底回流对象，假设控制通道与扰动通道的动态特性传递函数模型为:，15.0)()(2sfQsVHSp 2()10dTsV（c）对于温度对象，假设控制通道与扰动通道的动态特性传递函数模型为；3221.04sRssS eGe；FTsRRd 0.3exp1.3)()( 17ssFTsGSSd 0.2exp10.45)()( （e）温度测量环节的传递函数模型：，0.31().7sRTse 0.4.1()6STGse4、切换到手动状态，得出各开环阶跃响应曲线(所有控制变量及干扰量初始状态都为 0)（a）开始 10min 后，仅给 TR 控制信号+50 的阶跃信号；再过 40min 后，给 TR 控制信号为-

22、50 的阶跃信号（横轴表示时间，纵轴表示幅度）：TR 阶跃响应曲线 TS 阶跃响应曲线（b）开始 10min 后，仅给 TS 控制信号+50 的阶跃信号；再过 40min 后，给 TS 控制信号为-50 的阶跃信号（横轴表示时间，纵轴表示幅度）：TR 阶跃响应曲线 TS 阶跃响应曲线（c）开始 10min 后，仅给干扰量 F+20 的阶跃信号；再过 40min 后，给干扰量 F 为-20 的阶跃信号（横轴表示时间，纵轴表示幅度）：18进料量阶跃作用的 TR 响应曲线 进料量阶跃作用的 TS 阶跃响应曲线（d）开始 10min 后，仅给干扰量 PP+0.02 的阶跃信号；再过 40min 后，给

23、干扰量 PP 为-0.02 的阶跃信号（横轴表示时间，纵轴表示幅度）：泵出口压力阶跃作用的 TR 响应曲线 泵出口压力阶跃作用的 TS 响应曲线（e）开始 10min 后，仅给干扰量 V+15 的阶跃信号；再过 40min 后，给干扰量 V 为-15的阶跃信号（横轴表示时间，纵轴表示幅度）：塔底回流量阶跃作用的 TR 响应曲线 塔底回流量阶跃作用的 TS 响应曲线由以上图形相互比较，可以得出：通过 N21(s)、N 12(s)前馈解耦环节解耦以后，被控变量 TR、T S 的控制回路都转变为单回路控制通道。控制信号、干扰量都只对其控制通道有作19用。但是，由于每条回路都是由并联的过程特性函数组成

24、，由于过程特性函数滞后时间不一样，存在过渡时间，使得阶跃响应曲线有略微的振荡。5、控制器参数的整定（a）控制 TR 的控制器 PID1 的参数为：P=0.8，I=0.7，D=0.6（b）控制 TS 的控制器 PID2 的参数为：P=4.6，I=0.9，D=0.56、控制系统处于“闭环”状态下的性能测试（a）温度设定值跟踪响应试验（温度设定值为 140）：温度设定值 TR 响应曲线 温度设定值 TS 响应曲线（b）泵出口压力 PP 对系统输出的扰动响应试验（开始 40min 后，给阶跃信号+0.02MPa）：泵出口压力干扰 TR 响应曲线 泵出口压力干扰 TS 响应曲线20（c）进料量 F 对系

25、统输出的扰动响应试验（开始 40min 后，给阶跃信号+20T/hr）：进料量干扰 TR 响应曲线 进料量干扰 TS 响应曲线（d）塔底回流量 V 对系统输出的扰动响应试验（开始 40min 后，给阶跃信号+15T/hr）：塔底回流量干扰 TR 响应曲线 塔底回流量干扰 TS 响应曲线（f）三种干扰量同时对系统输出的扰动响应试验（开始 40min 后，三种干扰量同时给阶跃信号）：三种干扰同时作用 TR 响应曲线 三种干扰同时作用 TS 响应曲线由以上图形比较可以得出结论：控制系统在“闭环”状态下，温度设定值跟踪响应性能好，响应时间短；干扰量 PP 与 F 可以通过反馈控制迅速消除，保持系统的稳

26、定；通过解21耦环节以后，干扰量 V 不再对被控变量 TR 造成影响。但干扰量 V 仍让是被控变量 TS 的主要干扰量，仅靠反馈控制可以消除，但是超调量大，调节时间长，不利于系统安全。7、变换对象特性传递函数，测定系统“闭环”动态性能，RL()0.241psGf2()()20RdpLsGP，2().()05HSpVQsfs2()15dTsV；42.52.541.3sRssS eGe；()0()xp41RRdTsF.e3()SSdGss（a）温度设定值跟踪响应试验（温度设定值为 140）：温度设定值 TR 响应曲线 温度设定值 TS 响应曲线（b）泵出口压力 PP 对系统输出的扰动响应试验（开始

27、 40min 后，给阶跃信号+0.02MPa）：22泵出口压力干扰 TR 响应曲线 泵出口压力干扰 TS 响应曲线（c）进料量 F 对系统输出的扰动响应试验（开始 40min 后，给阶跃信号+20T/hr）：进料量干扰 TR 响应曲线 进料量干扰 TS 响应曲线（d）塔底回流量 V 对系统输出的扰动响应试验（开始 40min 后，给阶跃信号+15T/hr）：塔底回流量干扰 TR 响应曲线 塔底回流量干扰 TS 响应曲线通过以上图形，可以看出：塔顶回路可以维持稳定，但是对于塔底回路来说，无论是温度设定值跟踪响应曲线还是扰动响应曲线，都发生强烈的振荡。所以在各控制器参数均保持不变的前提下，当对象特

28、性在发生变化时，系统将可能出现不稳定状况，即每当对象特性发生变化时，我们都要进行 PID 整定，才能维持系统的稳定，才能使系统获得最佳23状态。四、 控制性能比较表 1 不同方案的控制性能比较表格方案一 方案二TR TS TR TS温度设定值跟踪衰减比 2.1：1 6：1 3.5:1 3:1超调量/ 63 36 35 9调节时间/min 14 19.5 11.5 11.5是否稳定 是 是 是 是泵压力干扰影响最大动态偏差/ 7 3.5 6.6 0恢复时间/min 3 4 3.4 0是否稳定 是 是 是 是进料量干扰影响最大动态偏差/ 2.5 6 2.8 6恢复时间/min 5.3 9 6.3

29、12是否稳定 是 是 是 是塔底回流量干扰影响最大动态偏差/ 16 150 0 150恢复时间/min 25 25.6 0 20.5是否稳定 是 是 是 是由表格中的数据比较可以得出：控制系统未解耦时，温度设定值跟踪时超调量大，调节时间长，相对于解耦的系统，稳定性要差一些。而且，干扰量会同时对两个被控变量造成影响。解耦之后，控制通道的干扰量只会对自己通道的控制变量造成影响，便于控制以及调整；另外，解耦后相当于单回路控制回路，PID 整定简单。24五、 心得体会过程控制是自动化专业的一门重要分支，也是自动化本科学生学习的重点。过程控制的课程设计是精馏塔温度控制系统设计与仿真，是一个比较实际，切合

30、实际工业生产的课题。课程设计综合了过程控制系统的各种理论，涉及到各种总控制系统的，像简单的 pid 控制，解耦，smith 预估，前馈控制，串级控制等，是检验过程控制课程学习的很好的试金石。在整个过孔课程设计中，最难的就是控制方法的确定，究竟应该选择怎样的控制方法。控制方法的选择不单单只是考虑控制效果，还应该考虑实际生产的成本控制，是否简单易行等各种因素。因此一开始我选择一个简单的 pid 控制系统，看看简单系统的控制是否能达到设计的要求。再者，本设计中涉及两个被控变量，两个输入，很明显需要检验是否存在强耦合关系，所以解耦的设计也是十分必要的。结果发现，控制系统未解耦时，温度设定值跟踪时超调量大，调节时间长，相对于解耦的系统，稳定性要差一些。而且，干扰量会同时对两个被控变量造成影响。解耦之后，控制通道的干扰量只会对自己通道的控制变量造成影响，便于控制以及调整；另外，解耦后相当于单回路控制回路，PID 整定简单。此题中还是存在强耦合关系，实际应用这应该考虑解耦。本次课程设计，时间不是很长，也比较赶，中间还穿插着计控课程设计，plc 创新实践，生产实习（对于我来说还要考虑实习等）能按时完成任务也为自己点个赞，希望以后的设计都要认真做好。继续努力。