CN110347038A - 一种水泥熟料冷却过程的二自由度Smith预估控制方法 - Google Patents
一种水泥熟料冷却过程的二自由度Smith预估控制方法 Download PDFInfo
- Publication number
- CN110347038A CN110347038A CN201910608000.2A CN201910608000A CN110347038A CN 110347038 A CN110347038 A CN 110347038A CN 201910608000 A CN201910608000 A CN 201910608000A CN 110347038 A CN110347038 A CN 110347038A
- Authority
- CN
- China
- Prior art keywords
- cooling procedure
- follows
- clinker cooling
- control
- controller
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/048—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators using a predictor
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种水泥熟料冷却过程的二自由度Smith预估控制方法。本发明首先法改进了传统的控制结构,基于大林控制算法,选择期望的闭环传递函数,进而得到相应的PI/PID控制器参数形式,使***能够有较好的设定点跟踪。然后,通过将Smith预估控制算法和内模控制算法有机结合,完成干扰抑制控制器的设计。使整个***能够使***同时满足良好的跟踪性能和干扰抑制性能,***的控制要求也能够得到保证。
Description
技术领域
本发明属于自动化技术领域,涉及一种水泥熟料冷却过程的二自由度Smith预估控制方法。
背景技术
在实际工业过程控制中,PID控制器具有结构简单、整定方便的优点,且能够满足大部分的生产控制要求。但是随着生产规模的扩大,控制精度要求的提高,普通的PID控制器已经不足以满足现在的生产要求。对于带有时滞的模型不确定/不匹配过程,设计的控制方法往往比较复杂,又不能同时兼顾设定点跟踪特性和抗干扰特性,因此研究一种二自由度Smith预估控制方法是有必要的。
发明内容
本发明主要是针对传统控制策略在水泥熟料冷却控制过程中所出现的***响应时滞、控制精度不足、模型不确定/不匹配等问题,提出了一种水泥熟料冷却过程的二自由度Smith预估控制方法。
本发明在基于Smith控制结构设计的基础上,设计了一种新的二自由度Smith预估控制结构,对于一阶和二阶时滞过程,通过选择传统的PI/PID形式和期望的闭环传递函数形式,得到相应的设定点跟踪控制器的各个控制参数。然后将这种Smith控制算法和内模控制设计方法有机的结合,选择合适的内部滤波器,得到干扰抑制控制器的设计。
本发明的具体技术方案是:
步骤1、改进二自由度Smith控制结构,具体是:
1.1设***的输出y与设定点输入r之间的闭环传递函数为:
其中,Gc(s)、Q(s)、Gp(s)、Gm(s)、Gm0(s)分别表示***前馈控制器、干扰抑制控制器、被控对象、内部模型和不包含时间延迟项的内部模型传递函数,理想状态下Gp(s)=Gm(s)。
1.2***的输出y与干扰d之间的闭环传递函数为:
步骤2、根据***的一、二阶加时滞过程模型,设计二自由度Smith预估控制器并进行参数的整定,具体是:
2.1首先建立***的一阶加时滞过程模型传递函数:
其中,K是过程增益,T是过程时间常数,θ是延迟时间。
2.2假设***期望的闭环传递函数Gq(s)为:
其中,η表示可调节的控制参数。
2.3根据步骤2.2,选取PI形式设计***的前馈控制器Gc(s):
其中,KP和TI分别表示为比例增益和积分时间常数,且
2.4利用内模控制原理,将过程模型Gm(s)分解为:
Gm(s)=Gm0(s)*G+(s)
其中,G+(s)是过程模型中不可逆的部分。
2.5基于***稳定性要求设计干扰抑制控制器Q(s)为:
其中,F(s)为内模控制低通滤波器,形式如下:
λ是一个调节参数;r是阶次;αi的约束为:这里,p1,...,pn表示过程被控对象中存在的极点。
2.6根据步骤2.5计算得到α值:
α=T[1-(1-λ/T)2e-θ/T]
2.7然后考虑工业设备的不确定性,建立***的二阶加时滞过程模型如下所示:
其中,T1、T2为过程模型时间常数。
2.8根据步骤2.2,并考虑二阶***的振荡特性选取PID形式设计***前馈控制器Gc(s):
其中,可得控制器的各个参数为:
2.9基于二阶过程,选择滤波器的形式为:
2.10根据步骤2.9可以设计二阶过程的干扰抑制控制器Q(s)为:
2.11根据步骤2.5中的αi的约束计算得到α1、α2值(其中p1=-1/T1,p2=-1/T2),然后得到控制器Q(s)并作用于水泥熟料冷却过程的篦冷机。
本发明提出了一种水泥熟料冷却过程的二自由度Smith预估控制方法。首先,该方法改进了传统的控制结构,基于大林控制算法,选择期望的闭环传递函数,进而得到相应的PI/PID控制器参数形式,使***能够有较好的设定点跟踪。然后,通过将Smith预估控制算法和内模控制算法有机结合,完成干扰抑制控制器的设计。使整个***能够使***同时满足良好的跟踪性能和干扰抑制性能,***的控制要求也能够得到保证。
与传统的方法相比,本发明所提出的新的二自由度Smith预估控制方法,能够同时实现良好的设定点跟踪和干扰抑制性能,并且设计简单,具有针对性,控制精度得到了很大的提高。
附图说明
图1为Smith控制框图;
图2为二自由度Smith预估控制框图。
具体实施方式
本发明方法主要是针对水泥生产线上熟料冷却过程的应用,熟料冷却过程是通过水泥熟料在篦冷机篦床上的压力变化来对篦冷机篦速进行调节的,从而保证水泥熟料的充分冷却。当篦床压力过大时控制篦速增大,从而减小水泥熟料厚度(篦压);当篦床压力较小时控制篦速减小,进而增大水泥熟料厚度。篦压为***输出量、篦冷机为***被控对象、篦冷机篦速为***控制输入量,***需要设定适当的篦压大小。
本发明方法的步骤包括:
步骤1、设计Smith控制结构,具体步骤是:
1.1首先,传统的Smith控制结构如图1所示。
1.2根据步骤1.1可知,水泥熟料冷却过程的输出篦压y与设定篦压r之间的闭环传递函数为:
其中,Gc(s)、Gp(s)、Gm(s)、Gm0(s)分别表示前馈控制器、篦冷机被控对象、篦冷机近似数学模型、不包含时间延迟项的篦冷机近似模型传递函数,r、d、y分别表示设定篦压、干扰输入、***输出篦压。
1.3水泥熟料冷却过程的输出篦压y与干扰输入d之间的闭环传递函数为:
1.4如果模型准确,即Gp(s)=Gm(s)时:
可以知道,输出的闭环传递函数的特征方程中并没有过程模型的时间延迟项。
步骤2、改进的二自由度Smith控制结构的设计,具体步骤是:
2.1二自由度Smith控制结构如图2所示。
2.2根据步骤2.1可知,水泥熟料冷却过程的输出篦压y与设定篦压r之间的闭环传递函数为:
2.3水泥熟料冷却过程的输出篦压y与干扰输入d之间的闭环传递函数为:
2.4假设模型准确,即Gp(s)=Gm(s)时:
Gyr(s)=Gc’(s)Gp(s)
Gyd(s)=[1-Gm(s)Q(s)]Gp(s)
可以知道,Gyr(s)和Gc’(s)是相关的,而Gc’(s)是和前馈控制器Gc(s)相关的,Gyd(s)和Q(s)是相关的,而且这两个控制器Gc(s)和Q(s)是可以单独的进行设计。
考虑到水泥熟料冷却过程在实际工业上所采用的篦冷机型号的不确定性,或篦冷机数学模型的不匹配性,接下来我们分别针对一、二阶惯性加时滞过程模型来分析水泥熟料冷却过程的***控制策略。
步骤3、根据水泥熟料冷却过程的一、二阶惯性加时滞过程模型,设计二自由度Smith预估控制器并进行参数的整定,具体是:
3.1首先建立水泥熟料冷却过程的一阶加时滞过程模型传递函数传递函数:
其中,K是冷却过程增益,T是冷却过程时间常数,θ是控制延迟时间。
3.2假设水泥熟料冷却过程期望的闭环传递函数Gq(s)为:
其中,η表示可调节的控制参数。
3.3根据步骤3.2,选取PI形式设计水泥熟料冷却过程的前馈控制器Gc(s):
其中,KP和TI分别表示为前馈控制器的比例增益和积分时间常数,且TI=T。
3.4利用内模控制原理,将水泥熟料冷却过程模型Gm(s)分解为:
Gm(s)=Gm0(s)*G+(s)
其中,G+(s)是水泥熟料冷却过程模型中不可逆的部分。
3.5基于水泥熟料冷却过程稳定性要求设计干扰抑制控制器Q(s)为:
其中,F(s)为内模控制低通滤波器,形式如下:
λ是一个调节参数;r是阶次;αi的约束为:这里,p1,...,pn表示过程被控对象中存在的极点。
3.6根据步骤3.5计算得到α值:
α=T[1-(1-λ/T)2e-θ/T]
3.7然后考虑篦冷机型号的不确定性,或篦冷机数学模型的不匹配性,建立***的二阶加时滞过程模型如下所示:
其中,T1、T2为水泥熟料冷却过程模型时间常数。
3.8根据步骤3.2,并考虑二阶水泥熟料冷却过程***的振荡特性选取PID形式设计***前馈控制器Gc(s):
其中,可得控制器的各个参数为:
3.9基于二阶过程,选择滤波器的形式为:
3.10根据步骤3.9可以设计二阶过程的干扰抑制控制器Q(s)为:
3.11根据步骤3.5中的αi的约束计算得到α1、α2值(其中p1=-1/T1,p2=-1/T2),然后得到控制器Q(s)并作用于水泥熟料冷却过程的篦冷机。
3.12在获得控制器参数后,有时可能需要进一步的微调以获得完美的控制器作用于水泥熟料冷却过程的篦冷机。
Claims (1)
1.一种水泥熟料冷却过程的二自由度Smith预估控制方法,其特征在于该方法包括以下步骤:
步骤1、设计Smith控制结构,具体是:
设水泥熟料冷却过程的输出篦压y与设定篦压r之间的闭环传递函数为:
其中,
Gc(s)、Gp(s)、Gm(s)、Gm0(s)分别表示前馈控制器、篦冷机被控对象、篦冷机近似数学模型、不包含时间延迟项的篦冷机近似模型传递函数,r、d、y分别表示设定篦压、干扰输入、***输出篦压;
水泥熟料冷却过程的输出篦压y与干扰输入d之间的闭环传递函数为:
步骤2、改进的二自由度Smith控制结构的设计,具体步骤是:
设水泥熟料冷却过程的输出篦压y与设定篦压r之间的闭环传递函数为:
水泥熟料冷却过程的输出篦压y与干扰输入d之间的闭环传递函数为:
步骤3、根据水泥熟料冷却过程的一、二阶惯性加时滞过程模型,设计二自由度Smith预估控制器并进行参数的整定,具体是:
3.1 首先建立水泥熟料冷却过程的一阶加时滞过程模型传递函数:
其中,K是冷却过程增益,T是冷却过程时间常数,θ是控制延迟时间;
3.2 假设水泥熟料冷却过程期望的闭环传递函数Gq(s)为:
其中,η表示可调节的控制参数;
3.3 根据步骤3.2,选取PI形式设计水泥熟料冷却过程的前馈控制器Gc(s):
其中,KP和TI分别表示为前馈控制器的比例增益和积分时间常数,且TI=T;
3.4 利用内模控制原理,将水泥熟料冷却过程模型Gm(s)分解为:
Gm(s)=Gm0(s)*G+(s)
其中,G+(s)是水泥熟料冷却过程模型中不可逆的部分;
3.5 基于水泥熟料冷却过程稳定性要求设计干扰抑制控制器Q(s)为:
其中,F(s)为内模控制低通滤波器,形式如下:
λ是一个调节参数;r是阶次;αi的约束为:
这里,p1,...,pn表示过程被控对象中存在的极点;
3.6 根据步骤3.5计算得到α值:
α=T[1-(1-λ/T)2e-θ/T]
3.7 建立***的二阶加时滞过程模型如下所示:
其中,T1、T2为水泥熟料冷却过程模型时间常数;
3.8 根据步骤3.2,并考虑二阶水泥熟料冷却过程***的振荡特性选取PID形式设计***前馈控制器Gc(s):
其中,可得控制器的各个参数为:
3.9 基于二阶过程,选择滤波器的形式为:
3.10 根据步骤3.9设计二阶过程的干扰抑制控制器Q(s)为:
3.11 根据步骤3.5中的αi的约束计算得到α1、α2值,然后得到控制器Q(s)并作用于水泥熟料冷却过程的篦冷机。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910608000.2A CN110347038B (zh) | 2019-07-08 | 2019-07-08 | 一种水泥熟料冷却过程的二自由度Smith预估控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910608000.2A CN110347038B (zh) | 2019-07-08 | 2019-07-08 | 一种水泥熟料冷却过程的二自由度Smith预估控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110347038A true CN110347038A (zh) | 2019-10-18 |
CN110347038B CN110347038B (zh) | 2022-03-11 |
Family
ID=68177956
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910608000.2A Active CN110347038B (zh) | 2019-07-08 | 2019-07-08 | 一种水泥熟料冷却过程的二自由度Smith预估控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110347038B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110879576A (zh) * | 2019-10-15 | 2020-03-13 | 杭州电子科技大学 | 水泥熟料冷却过程的分数阶模型广义二自由度控制方法 |
CN111381494A (zh) * | 2020-05-02 | 2020-07-07 | 苏州科技大学 | 一种基于极点配置的空调温控***及其控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5933345A (en) * | 1996-05-06 | 1999-08-03 | Pavilion Technologies, Inc. | Method and apparatus for dynamic and steady state modeling over a desired path between two end points |
CN101751051A (zh) * | 2008-12-05 | 2010-06-23 | 中国科学院沈阳自动化研究所 | 基于约束史密斯广义预测控制的水泥分解炉温度控制方法 |
CN105865215A (zh) * | 2016-04-18 | 2016-08-17 | 重庆大学 | 水泥窑炉温度多参量控制*** |
CN109116738A (zh) * | 2018-09-27 | 2019-01-01 | 杭州电子科技大学 | 一种工业加热炉的二自由度内模控制分析方法 |
CN109557810A (zh) * | 2018-11-29 | 2019-04-02 | 杭州电子科技大学 | 一种基于新型二自由度内模pid的加热炉温度控制方法 |
-
2019
- 2019-07-08 CN CN201910608000.2A patent/CN110347038B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5933345A (en) * | 1996-05-06 | 1999-08-03 | Pavilion Technologies, Inc. | Method and apparatus for dynamic and steady state modeling over a desired path between two end points |
CN101751051A (zh) * | 2008-12-05 | 2010-06-23 | 中国科学院沈阳自动化研究所 | 基于约束史密斯广义预测控制的水泥分解炉温度控制方法 |
CN105865215A (zh) * | 2016-04-18 | 2016-08-17 | 重庆大学 | 水泥窑炉温度多参量控制*** |
CN109116738A (zh) * | 2018-09-27 | 2019-01-01 | 杭州电子科技大学 | 一种工业加热炉的二自由度内模控制分析方法 |
CN109557810A (zh) * | 2018-11-29 | 2019-04-02 | 杭州电子科技大学 | 一种基于新型二自由度内模pid的加热炉温度控制方法 |
Non-Patent Citations (1)
Title |
---|
邓立广 等: "基于smith预估的水泥分解炉温度MRFAC控制器设计", 《化工自动化及仪表》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110879576A (zh) * | 2019-10-15 | 2020-03-13 | 杭州电子科技大学 | 水泥熟料冷却过程的分数阶模型广义二自由度控制方法 |
CN111381494A (zh) * | 2020-05-02 | 2020-07-07 | 苏州科技大学 | 一种基于极点配置的空调温控***及其控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110347038B (zh) | 2022-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Veenman et al. | Stability analysis with integral quadratic constraints: A dissipativity based proof | |
CN109557810B (zh) | 一种基于新型二自由度内模pid的加热炉温度控制方法 | |
CN110347038A (zh) | 一种水泥熟料冷却过程的二自由度Smith预估控制方法 | |
Weyer | Control of irrigation channels | |
CN105388764B (zh) | 基于动态矩阵前馈预测的电液伺服pid控制方法及*** | |
WO2021027449A1 (zh) | 考虑非定尺订单规格柔性的热轧中厚板组板及板坯设计的方法及模型*** | |
CN109116738B (zh) | 一种工业加热炉的二自由度内模控制分析方法 | |
CN109960149A (zh) | 用于pH中和过程的自抗扰控制参数整定方法 | |
CN108803310A (zh) | 一种pid控制方法、装置及设备 | |
CN104199437A (zh) | 基于区域极点指标的分数阶PIλDμ控制器的参数优化方法 | |
CN109521677A (zh) | 一种基于预测模型的控制方法和控制器 | |
CN104747294A (zh) | 燃气轮机燃料压力控制方法和*** | |
CN105867194A (zh) | 设备控制方法和装置 | |
CN109976405A (zh) | 一种陶瓷窑炉温度复合控制方法、设备及*** | |
CN106647240A (zh) | 基于导前扰动模型的亚临界机组协调预测函数控制算法 | |
CN110308647A (zh) | 含误差积分输入项的无人机三段式模糊pid控制方法 | |
CN105487375B (zh) | 一种离散pid控制器参数整定方法 | |
TWI225977B (en) | Control method and control apparatus | |
CN107704676B (zh) | 金属氧化物界面调控的方法 | |
CN107728481A (zh) | 一种基于模型预测控制的闭环建模方法及装置 | |
CN106026183B (zh) | 水轮机调速器孤网运行控制方法及*** | |
Zhmud et al. | Identification and control of object with time-delay link | |
CN104668493B (zh) | 一种连铸机结晶器的复合振动控制***的控制方法 | |
CN108132596B (zh) | 一种微分超前广义智能内部模型集pid控制器设计方法 | |
CN107870567A (zh) | 一种比例微分超前广义智能内部模型集pid控制器设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |