首页 > 专利 > 杭州电子科技大学 > 一种化工时变工业过程混合控制方法专利详情

一种化工时变工业过程混合控制方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2018-07-11
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2018-12-18
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2021-02-09
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2038-07-11
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201810760549.9 申请日 2018-07-11
公开/公告号 CN108873699B 公开/公告日 2021-02-09
授权日 2021-02-09 预估到期日 2038-07-11
申请年 2018年 公开/公告年 2021年
缴费截止日
分类号 G05B13/04 主分类号 G05B13/04
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 0
权利要求数量 1 非专利引证数量 1
引用专利数量 0 被引证专利数量 0
非专利引证 1、CN 107831662 A,2018.03.23CN 107831662 A,2018.03.23CN 106896800 A,2017.06.27CN 107544255 A,2018.01.05JP 2013069094 A,2013.04.18CN 107991875 A,2018.05.04CN 107991875 A,2018.05.04Limin Wang et al..Robust delaydependent iterative learning fault-tolerant control for batch processes withstate delay and actuator failures. 《Journal of process control》.2012,第22卷(第7期),第1273-1286页. Limin Wang et al..Average dwell time-based optimal iterative learning controlfor multi-phase batch processes《.Journalof Process Control》.2016,第1-12页. Limin Wang et al..Robust design offeedback integrated with iterativelearning control for batch processes withuncertainties and interval time-varyingdelays《.Journal of Process Control》.2011,第21卷(第7期),第987-996页. 王玉中.工业过程的预测控制与模糊PID控制的研究《.中国优秀硕士学位论文全文数据库信息科技辑》.2018,第I140-1189页.;
引用专利 被引证专利
专利权维持 4 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 胡晓敏、李容轩、邹洪波 第一发明人 胡晓敏
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 3
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
杭州浙科专利事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
吴秉中
摘要
本发明公开了一种化工时变批次过程混合控制方法,包括如下步骤:步骤1、建立批次过程时变状态空间模型;步骤2、设计被控对象的批次过程控制器。该方法首先建立批次过程模型,通过引入状态误差和输出误差,将上述模型转化为等效的随机系统模型,根据不同故障发生的概率,将常规的迭代学习控制律设计转化为更灵活的更新律设计。不同于传统的控制策略,本发明所提出的混合控制策略考虑到了执行器出现不同故障的概率,系统全面地分析、处理各类故障,故障处理的灵活性、快速性更好。
  • 摘要附图
    一种化工时变工业过程混合控制方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2021-02-09 授权
2 2018-12-18 实质审查的生效 IPC(主分类): G05B 13/04 专利申请号: 201810760549.9 申请日: 2018.07.11
3 2018-11-23 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.一种化工时变批次过程混合控制方法,包括如下步骤:
步骤1、建立批次过程时变状态空间模型;
步骤2、设计被控对象的批次过程控制器;
步骤1具体如下:
1-
1.建立一个批次过程系统模型,其形式如下:
其中k和t分别表示批次和批次运行时刻,x(t+1,k)、x(t,k)、x(t-d(t),k)分别是k批次t+1时刻、t时刻、t-d(t)时刻的系统状态,d(t)是系统t时刻的状态延迟,dm≤d(t)≤dM,dm、dM分别是状态延迟的下限和上限,y(t,k)∈Rl是k批次t时刻的系统输出,维数为Rl,u(t,k)∈Rm是k批次t时刻的系统输入,维数为Rm,l,m分别是系统输出和输入的阶次,σ(t,k)表示与批次和时刻有关的切换信号,Aσ(t,k),Adσ(t,k),Bσ(t,k),Cσ(t,k)分别表示带有切换信号的适当维度的常数矩阵,ωσ(t,k)(t,k)是k批次t时刻的外部扰动,x(0,k)是k批次系统的初始状态,其初始值设置为x0,k;
1-
2.让批次过程输出跟踪给定的期望轨迹,其定义如下:
其中yr(t)是t时刻系统输出期望轨迹,e(t,k)是k批次t时刻的系统输出误差, 表示‘定义为’;
1-
3.批次过程系统故障的发生概率的定义如下:
0≤P{γ(t+1,k)=1|γ(t,k)=0}=α≤1,0≤P{γ(t+1,k)=0|γ(t,k)=0}=1-α≤1
0≤P{γ(t+1,k)=1|γ(t,k)=1}=1-χ≤1,0≤P{γ(t+1,k)=0|γ(t,k)=1}=χ≤1其中γ(t,k)、γ(t+1,k)分别表示k批次t时刻、t+1时刻系统故障判定函数,取0表示系统正常,取1表示系统故障,α表示系统当前时刻正常运行但下一时刻发生故障的概率,χ表示当前时刻系统故障但下一时刻恢复正常运行的概率;
1-
4.确定判断每批次故障的发生是否与当前时刻有关的概率矩阵:
首先,定义状态传递概率矩阵为
其中p00=1-α,p01=α,p10=χ,p11=1-χ;
进而,得到批次间状态变化的n步传递函数概率矩阵Pn;
1-
5.批次过程系统的随机迭代学习控制律描述如下:
其中Δu(t,k)表示k批次t时刻的系统随机迭代学习输入更新律,u(t,0)表示开始批次t时刻的系统输入,并设置为0;
1-
6.批次过程系统可能发生故障的情况下,系统输入和系统状态误差分别如下:u(t,k)=(1-γ(t,k))u(t,k-1),δ(x(t,k))=x(t,k)-x(t,k-1)
其中x(t,k-1)表示k-1批次t时刻的系统状态,δ(x(t,k))表示k批次t时刻的系统状态误差;
1-
7.根据1-1及1-6,得到下列状态误差和输出误差表达式:
其中 表示拓展的外部扰动,e(t+1,k)、e(t+1,k-1)是k批次、k-1批次t+1时刻的系统输出误差,y(t+1,k)、yr(t+1,k)分别是k批次t+1时刻的系统输出和系统期望输出,δ(x(t+1,k))、δ(x(t-d(t),k))分别是k批次t+1时刻、t-d(t)时刻的系统状态误差,A、Ad、B、C分别是适当维度的常数矩阵;
步骤2具体如下:
2-
1.基于步骤1,进一步得到改写后的系统误差模型如下:
其中 e(t+1-d(t),k-1)是k-1批次t+1-d(t)时刻的系统输出误差,z(t,k)表示
k批次t时刻的系统整体误差,I为适当维度的酉矩阵;
2-
2.进一步得到系统的随机迭代学习更新律如下:
Δu(t,k)=(1-γ(t,k))K0X(t,k)
其中 K0为满足系统要求的增益矩阵;
2-
3.在下一时刻,重复步骤2-1到2-2继续求解新的最优系统随机迭代学习更新律Δu(t,k),得到最优控制量,作用于控制对象,并依次循环。
说明书

技术领域

[0001] 本发明属于自动化工业过程控制领域,涉及到一种化工时变批次过程混合控制方法。

背景技术

[0002] 在工业生产过程中,批次处理过程非常普遍,同时,在复杂的工业生产环境下,长时间运行的生产设备出现故障的情况很普遍,存在的故障,不仅会影响生产效率和产品质量,还会造成财产损失和人员伤亡。综合考虑安全生产和经济效益,当系统出现故障时,系统仍要保持一定的稳定性和可控性。因此,有必要对故障处理方法进行研究。

发明内容

[0003] 本发明目的是为了更好地解决化工批次过程中执行器出现的故障,提出了一种化工时变批次过程混合控制方法。该方法首先建立批次过程模型,通过引入状态误差和输出误差,将上述模型转化为等效的随机系统模型,根据不同故障发生的概率,将常规的迭代学习控制律设计转化为更灵活的更新律设计。不同于传统的控制策略,本发明所提出的混合控制策略考虑到了执行器出现不同故障的概率,系统全面地分析、处理各类故障,故障处理的灵活性、快速性更好。
[0004] 本发明的技术方案是通过模型建立、控制器设计、预测机理、优化等手段,设计了一种化工时变批次过程混合控制方法,利用该方法可以提高系统的安全性和可靠性。其具体技术方案如下:
[0005] 本发明方法的步骤包括:
[0006] 步骤1、建立批次过程时变状态空间模型,具体方法是:
[0007] 1-1.建立一个批次过程系统模型,其形式如下:
[0008]
[0009] 其中k和t分别表示批次和批次运行时刻,x(t+1,k)、x(t,k)、x(t-d(t),k)分别是k批次t+1时刻、t时刻、t-d(t)时刻的系统状态,d(t)是系统t时刻的状态延迟,dm≤d(t)≤dM,dm、dM分别是状态延迟的下限和上限,y(t,k)∈Rl是k批次t时刻的系统输出,维数为Rl,u(t,k)∈Rm是k批次t时刻的系统输入,维数为Rm,l,m分别是系统输出和输入的阶次,σ(t,k)表示与批次和时刻有关的切换信号,Aσ(t,k),Adσ(t,k),Bσ(t,k),Cσ(t,k)分别表示带有切换信号的适当维度的常数矩阵,ωσ(t,k)(t,k)是k批次t时刻的外部扰动,x(0,k)是k批次系统的初始状态,其初始值设置为x0,k。
[0010] 1-2.让批次过程输出跟踪给定的期望轨迹,其定义如下:
[0011] e(t,k)A yr(t)-y(t,k)
[0012] 其中yr(t)是t时刻系统输出期望轨迹,e(t,k)是k批次t时刻的系统输出误差,A表示‘定义为’。
[0013] 1-3.批次过程系统故障的发生概率的定义如下:
[0014] 0≤P{γ(t+1,k)=1|γ(t,k)=0}=α≤1,0≤P{γ(t+1,k)=0|γ(t,k)=0}=1-α≤1
[0015] 0≤P{γ(t+1,k)=1|γ(t,k)=1}=1-χ≤1,0≤P{γ(t+1,k)=0|γ(t,k)=1}=χ≤1
[0016] 其中γ(t,k)、γ(t+1,k)分别表示k批次t时刻、t+1时刻系统故障判定函数,取0表示系统正常,取1表示系统故障,α表示系统当前时刻正常运行但下一时刻发生故障的概率,χ表示当前时刻系统故障但下一时刻恢复正常运行的概率。
[0017] 1-4.确定判断每批次故障的发生是否与当前时刻有关的概率矩阵。
[0018] 首先,定义状态传递概率矩阵为
[0019] 其中p00=1-α,p01=α,p10=χ,p11=1-χ。
[0020] 进而,得到批次间状态变化的n步传递函数概率矩阵Pn。
[0021] 1-5.批次过程系统的随机迭代学习控制律描述如下:
[0022]
[0023] 其中Δu(t,k)表示k批次t时刻的系统随机迭代学习输入更新律,u(t,0)表示开始批次t时刻的系统输入,并设置为0。
[0024] 1-6.批次过程系统可能发生故障的情况下,系统输入和系统状态误差分别如下:
[0025] u(t,k)=(1-γ(t,k))u(t,k),δ(x(t,k))=x(t,k)-x(t,k-1)
[0026] 其中x(t,k-1)表示k-1批次t时刻的系统状态,δ(x(t,k))表示k批次t时刻的系统状态误差。
[0027] 1-7.根据1-1及1-6,得到下列状态误差和输出误差表达式。
[0028]
[0029] 其中 表示拓展的外部扰动,e(t+1,k)、e(t+1,k-1)是k批次、k-1批次t+1时刻的系统输出误差,y(t+1,k)、yr(t+1,k)分别是k批次t+1时刻的系统输出和系统期望输出,δ(x(t+1,k))、δ(x(t-d(t),k))分别是k批次t+1时刻、t-d(t)时刻的系统状态误差,A、、Ad、B C分别是适当维度的常数矩阵。
[0030] 步骤2、设计被控对象的批次过程控制器,具体是:
[0031] 2-1.基于步骤1,进一步得到改写后的系统误差模型如下:
[0032]
[0033] 其中 e(t+1-d(t),k-1)是k-1批次t+1-d(t)时刻的系统输出误差,z(t,k)表示k批次t时刻的系统整体误差,I为适当维度的酉矩阵。
[0034] 2-2.进一步得到系统的随机迭代学习更新律如下:
[0035] Δu(t,k)=(1-γ(t,k))K0X(t,k)
[0036] 其中 K0为满足系统要求的增益矩阵。
[0037] 2-3.在下一时刻,重复步骤2.1到2.2继续求解新的最优系统随机迭代学习更新律Δu(t,k),得到最优控制量,作用于控制对象,并依次循环。

实施方案

[0038] 以注塑成型工艺为例:
[0039] 这里以注塑成型过程中的填料压力加以描述,调节手段是控制比例阀的阀门开度。
[0040] 步骤1、建立注塑成型过程时变状态空间模型,具体步骤是:
[0041] 1-1.建立一个注塑成型过程系统模型,其形式如下:
[0042]
[0043] 其中k和t分别表示注塑成型过程的批次和批次运行时刻,x(t+1,k)、x(t,k)、x(t-d(t),k)分别是注塑成型过程中k批次t+1时刻、t时刻、t-d(t)时刻的系统状态,d(t)是注塑成型过程中t时刻的状态延迟,dm≤d(t)≤dM,dm、dM分别是注塑成型过程中的状态延迟的下限和上限,y(t,k)∈Rl是注塑成型过程中k批次t时刻的填料压力,维数为Rl,u(t,k)∈Rm是注塑成型过程中k批次t时刻的阀门开度,维数为Rm,l,m分别是注塑成型过程中的填料压力和阀门开度的阶次,σ(t,k)表示注塑成型过程中与批次和时刻有关的切换信号,Aσ(t,k),Adσ(t,k),Bσ(t,k),Cσ(t,k)分别表示注塑成型过程中带有切换信号的适当维度的常数矩阵,ωσ(t,k)(t,k)是注塑成型过程中k批次t时刻的外部扰动,x(0,k)是注塑成型过程中k批次注塑成型的初始状态,其初始值设置为x0,k。
[0044] 1-2.让注塑成型过程的填料压力跟踪给定的填料压力轨迹,其定义如下:
[0045] e(t,k)A yr(t)-y(t,k)
[0046] 其中yr(t)是注塑成型过程中t时刻的给定填料压力轨迹,e(t,k)是注塑成型过程中k批次t时刻的填料压力误差,A表示‘定义为’。
[0047] 1-3.注塑成型过程故障的发生概率的定义如下:
[0048] 0≤P{γ(t+1,k)=1|γ(t,k)=0}=α≤1,0≤P{γ(t+1,k)=0|γ(t,k)=0}=1-α≤1
[0049] 0≤P{γ(t+1,k)=1|γ(t,k)=1}=1-χ≤1,0≤P{γ(t+1,k)=0|γ(t,k)=1}=χ≤1
[0050] 其中γ(t,k)、γ(t+1,k)分别表示注塑成型过程中k批次t时刻、t+1时刻注塑成型的系统故障判定函数,取0表示注塑成型过程运行正常,取1表示注塑成型过程运行故障,α表示注塑成型过程当前时刻正常运行但下一时刻发生故障的概率,χ表示注塑成型过程当前时刻故障但下一时刻恢复正常运行的概率。
[0051] 1-4.确定判断注塑成型过程中每批次故障的发生是否与当前时刻有关的概率矩阵。
[0052] 首先,定义注塑成型过程的状态传递概率矩阵为
[0053] 其中p00=1-α,p01=α,p10=χ,p11=1-χ。
[0054] 进而,得到注塑成型过程批次间状态变化的n步传递函数概率矩阵Pn。
[0055] 1-5.注塑成型过程的随机迭代学习控制律描述如下:
[0056]
[0057] 其中Δu(t,k)表示注塑成型过程中k批次t时刻的注塑成型随机迭代学习阀门开度的更新律,u(t,0)表示注塑成型过程中开始批次t时刻的阀门开度。
[0058] 1-6..注塑成型过程可能发生故障的情况下,注塑成型的阀门开度和注塑成型的状态误差分别如下:
[0059] u(t,k)=(1-γ(t,k))u(t,k),δ(x(t,k))=x(t,k)-x(t,k-1)
[0060] 其中x(t,k-1)表示注塑成型过程中k-1批次t时刻的系统状态,δ(x(t,k))表示注塑成型过程中k批次t时刻的系统状态误差。
[0061] 1-7.根据1-1及1-6,得到下列注塑成型过程的状态误差和填料压力误差表达式。
[0062]
[0063] 其中 表示注塑成型过程中拓展的外部扰动,e(t+1,k)、e(t+1,k-1)是注塑成型过程中k批次、k-1批次t+1时刻的填料压力误差,y(t+1,k)、yr(t+1,k)分别是注塑成型过程中k批次t+1时刻的填料压力和期望填料压力,δ(x(t+1,k))、δ(x(t-d(t),k))分别是注塑成型过程中k批次t+1时刻、t-d(t)时刻的系统状态误差,A、、Ad、B C分别是注塑成型过程中的适当维度的常数矩阵。
[0064] 步骤2、设计注塑成型过程的时变混合控制器,具体是:
[0065] 2-1.基于步骤1,进一步得到改写后的注塑成型过程的误差模型如下:
[0066]
[0067] 其中 e(t+1-d(t),k-1)是注塑成型过程中k-1批次t+1-d(t)时刻的填料压力误差,z(t,k)表示注塑成型过程中k批次t时刻的系统整体误差,I为适当维度的酉矩阵。
[0068] 2-2.进一步得到注塑成型过程的随机迭代学习更新律如下:
[0069] Δu(t,k)=(1-γ(t,k))K0X(t,k)
[0070] 其中 K0为满足注塑成型过程的增益矩阵。
[0071] 2-3.在下一时刻,重复步骤2.1到2.2继续求解新的最优注塑成型过程随机迭代学习更新律Δu(t,k),得到最优阀门开度,作用于注塑成型过程,并依次循环。
专利联系人(活跃度排行)
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号