电阻式处理炉温度控制系统设计

1、过程控制大作业 题 目:电阻式处理炉温度控制系统设计 班 级：0602303 学 号：060230328 姓 名： 任课教师： 完成时间：2009年12月10日 目 录 一、 设计任务及要求二、 被控对象数学模型的建立及对象特性分析2.1数学模型计算及仿真验证2.2对象特性分析三、 单回路控制系统设计3.1基本控制方案3.2 控制仪表选择 3.3 控制系统matlab仿真3.4 仿真结果分析四、串级控制系统设计4.1 控制方案选择依据4.2 控制仪表选择4.3 控制系统matlab仿真4.4 仿真结果分析一、 设计任务及要求对一个电阻式处理炉装置作如下实验，在温度控制稳定时，在开环状态下将执行

2、器的输入信号增加，持续后结束，记录炉内温度变化数据如下表，试根据实验数据设计一个超调量的无差温度控制系统。t(min)02468101214161820001.253.757.5010.509.207.256.004.803.70t(min)22241628303234363840422.902.251.851.501.150.850.600.400.300.200.10具体设计要求如下:（1） 根据实验数据选择一定的辨识方法建立对象的模型；（2） 根据辨识结果设计符合要求的控制系统（给出带控制点的控制流程图，控制系统原理图等，选择控制规律）；（3） 根据设计方案选择相应的控制仪表；（4） 对

3、设计系统进行仿真，整定运行参数。二被控对象数学模型的建立及对象特性分析2.1数学模型计算及仿真验证由设计要求与实验数据，可以写出系统阶跃响应数据，如下表：t（min）02468101214161820001.253.757.5010.509.207.256.004.803.70001255.0012.5023.0032.20394545.4550.2553.95t（min）22242628303234363840422.902.251.851.501.150.850.600.400.300.200.1056.8559.1060.9562.4563.6064.4565.0565.4565.756

4、5.9566.05使用matlab制作出阶跃响应和脉冲响应，并作出阶跃响应的归一化曲线。源程序如下：x=0:2:42;y1=0 0 1.25 3.75 7.50 10.50 9.20 7.25 6.00 4.80 3.70 2.90 2.25 1.85 1.50 1.15 0.85 0.60 0.40 0.30 0.20 0.10 ;y2=0 0 1.25 5.00 12.50 23.00 32.20 39.45 45.45 50.25 53.95 56.85 59.10 60.95 62.45 63.60 64.45 65.05 65.45 65.75 65.95 66.05 ;yw=66.

5、15;ys=y2./yw;x0=0:.025:42;y3=spline(x,y1,x0); % y 脉冲响应 y4=spline(x,y2,x0); % y1 阶跃响应 y5=spline(x,ys,x0); % ys 单位阶跃响应plot(x,y1,'r*',x0,y3);gtext('red - 脉冲响应');hold on;plot(x,y2,'g*',x0,y4);gtext('green - 阶跃响应');grid;figure;plot(x0,y5);gtext('单位阶跃响应 ');grid;结果：有

6、上面的归一化图形，可以得到当y=0.4时，t1=10.7min y=0.8时，t2=19.38min对象增益k=yw/i=(66.05+0.1)/3=22.05./mat1-tao/t2-tao=（10.7-2）/（19.38-2）=0.5>0.46,因此选三阶系统近似（n=3）。由公式n*t=(t1+t2)/2.16，3*t=(8.7+17.38)/2.16s,解得t=4.02min=241.2s由于延迟tao=120s，则对象传函为：2.2对象特性分析下面进行仿真验证，源程序如下：num1=22.05;den1=conv(241.2 1,conv(241.2 1,241.2 1);g

7、1=tf(num1,den1);num2,den2=pade(120,3);g=g1*tf(num2,den2);h=feedback(g,1);step(h);gtext('单位阶跃响应');grid;结果：明显的，这个系统是不稳定的系统，需要加入pid控制器进行调节。三单回路控制系统设计3.1基本控制方案3.1基本控制方案电阻炉温度控制系统是闭合的反馈系统，温度系统主要由温度传感器、温度调节仪、执行装置、被控对象四个部分组成。控制量是电阻炉温度，通过调节电流大小控制温度高低。控制系统选择比例积分微分pid调节，比例环节的引入是为了及时成比例地反映控制系统的偏差信号，以最快的

8、速度产生控制作用，使偏差向最小的方向变化。积分环节的引入主要用于消除静差，即当闭环系统处于稳定状态时，则此时控制输出量和控制偏差量都将保持在某一个常值上。微分环节的引入是为了改善系统的稳定性和动态响应速度，它能反映偏差信号的变化趋势，并能在偏差信号值变太大之前，在系统引入一个有效的早期修正信号，从而加速系统的动态速度，减小调节时间。电阻炉tctt系统原理图系统流程图由图，选择被控制量为电阻式处理炉内部温度，调节量为调节阀开度。3.2 控制仪表选择所有仪表均选用ddz-型表。因为大型电阻炉内部温度很高，所以采用热电偶。km>0,从安全角度考虑，调节阀为气开阀，则kv>0。对象kp&g

9、t;0,所以调节器kc>0。温度变送器量程为01000，调节阀范围0100t/min，则gv=（100-0）/(20-4) t/min/ma, gm=(20-4)/(1000-0)ma/.所以go=gv*gp*gm,则广义对象的传函为：3.3 控制系统matlab仿真设调节器传递函数为。应用临界比例带法对pid参数进行整定。源程序如下：num1=2.205;den1=conv(241.2 1,conv(241.2 1,241.2 1);g1=tf(num1,den1);num2,den2=pade(120,3);g=g1*tf(num2,den2);for kcr=0:0.01:500

10、%求临界比例增益gc=kcr; g2=feedback(gc*g,1);num,den=tfdata(g2,'v'); %获得传函的分子和分母p=roots(den); %求闭环特征根pr=real(p);prm=max(pr);pro=find(prm>=-0.00001);n=length(pro);if n>=1 breakend;endstep(g2,0:0.2:5000);grid onkcrkcr = 1.6000从图中可以得出tcr=2160-846=1314s；用临界比例带法确定pid参数kp、ti、tdnum1=2.205;den1=conv(24

11、1.2 1,conv(241.2 1,241.2 1);g1=tf(num1,den1);num2,den2=pade(120,3);g=g1*tf(num2,den2);kcr=1.6;tcr=1314;kp=0.6*kcr;ti=0.55*tcr;td=0.125*tcr;kp,ti,td,s=tf('s');gc=kp*(1+1/(ti*s)+td*s);g=feedback(gc*g,1);step(g)kp = 0.9600ti = 722.7000td = 164.2500kp=0.96 ti=722.7 td=164.253.4 仿真结果分析从上面一节的图形看来，

12、加入pid之后虽然系统最终变为稳定，但是超调量仍然达到21%。要把超调量降低，就要重新调整pid参数。为了减小超调，应该减小比例和积分作用，即减小kp，增大ti。编程在matlab下进行调整：num1=2.205;den1=conv(241.2 1,conv(241.2 1,241.2 1);g1=tf(num1,den1);num2,den2=pade(120,3);g=g1*tf(num2,den2)kp =0.86;ti =800;td =164.25;s=tf('s');gc=kp*(1+1/(ti*s)+td*s);g2=feedback(gc*g,1);step(g

13、2)结果为：从图中可以看出超调量降到了12%，符合要求。下面对系统的抗干扰性进行分析：波形如下：分析：再加入一个脉冲干扰信号后，输出超调60%，说明系统的抗干扰能力比较差。四串级控制系统设计4.1 控制方案选择依据（1）扰动分析：压力，流量，成分，热值（2）控制方案：因为流量信号在加热炉的后侧，有纯滞后环节，不宜采用流量作为副被调量，同样，如果以压力作为副被调量，因为压力检测元件多为细管，易堵塞，所以不采用，而以炉内膛温度作为副被调量。 4.2 控制仪表选择所有仪表均为ddz-型表。炉内温度作为副被调量，拥有克服燃料影响，如温度，成分等，温度点不易找准，实际中应多测几个点温，再去平均温度，。因

14、为炉内温度很高，所以采用热电偶。 串级控制系统原理图 串级控制系统流程图4.3 控制系统matlab仿真先整定副回路，设副对象暂定为，。因此，采用比例控制规律，利用经验法进行整定，编程如下：num=0.1;den=241.2 1;go2=tf(num,den);g2=feedback(2*go2,1);step(g2);经分析，发现kp2越大，ts越小，稳态误差越小，但kp2不能无限制的增加。单位阶跃响应曲线如下：然后整定外环，有刚刚确定的kp2，带入可得外环的广义传递函数为：，编程如下：num=0.1;den=241.2 1;go2=tf(num,den);g2=feedback(2*go2

15、,1);step(g2);num1=22.05;den1=conv(241.2 1,241.2 1);g1=tf(num1,den1);num2,den2=pade(120,3);gp1=g1*tf(num2,den2);go1=g2*gp1*0.016;for kcr=0:0.01:500 %求临界比例增益gc=kcr; g=feedback(gc*go1,1);num,den=tfdata(g,'v'); %获得传函的分子和分母p=roots(den); %求闭环特征根pr=real(p);prm=max(pr);pro=find(prm>=-0.00001);n=l

16、ength(pro);if n>=1 breakend;endstep(g,0:0.2:5000);grid onkcrkcr= 58.6200从图中可以得出tcr=2070-810=1260；kp=0.6*kcr=35.172 , ti=0.5*tcr=630, td=0.125*tcr=157.5编程在matlab下进行：num=0.1;den=241.2 1;go2=tf(num,den);g2=feedback(2*go2,1);step(g2);num1=22.05;den1=conv(241.2 1,241.2 1);g1=tf(num1,den1);num2,den2=pade(120,3);gp1=g1*tf(num2,den2);gp=g2*gp1*0.016;kp =35.172;ti =630;td =157.5;s=tf('s');gc=kp*(1+1/(ti*s)+td*s);g=feedba