CN109870908A - 一种用于连续搅拌釜式反应器系统的采样控制方法 - Google Patents
一种用于连续搅拌釜式反应器系统的采样控制方法 Download PDFInfo
- Publication number
- CN109870908A CN109870908A CN201910146505.1A CN201910146505A CN109870908A CN 109870908 A CN109870908 A CN 109870908A CN 201910146505 A CN201910146505 A CN 201910146505A CN 109870908 A CN109870908 A CN 109870908A
- Authority
- CN
- China
- Prior art keywords
- formula
- stirred tank
- tank reactor
- continuous stirred
- reactor system
- 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.)
- Pending
Links
Abstract
本发明公开了一种用于连续搅拌釜式反应器系统的采样控制方法,该方法包括以下步骤:首先建立连续搅拌釜式反应器系统的状态空间模型;然后对于一个采样周期内系统不发生切换时与一个采样周期内系统发生切换时的情况,分别构造相应的李雅普诺夫函数,并设置相应的控制率;最后设置采样周期与平均驻留时间以使构造的李雅普诺夫函数有界,证明闭环系统的稳定性,从而完成连续搅拌釜式反应器系统的采样控制方法。本发明的采样控制方法使连续搅拌釜式反应器系统在不同运行条件下具有更好的性能与鲁棒性,且方法简单,使控制成本显著降低,且更加便于用微处理器实现,具有较高的工程实用价值。
Description
技术领域
本发明属于非线性系统控制技术领域,特别涉及一种用于连续搅拌釜式反应器系统的采样控制方法。
背景技术
连续搅拌釜式反应器,简称为一种常见的非线性化学反应器,由于其具有热交换能力强、产品质量稳定和低成本等优点,在化工生产的核心设备中占有重要地位,被广泛应用于染料、医药试剂、食品及合成材料等工业生产中。系统的控制变量主要包括反应物的温度、浓度等,对这些变量的控制效果将直接影响到化工产品的生产质量。同时,连续搅拌釜式反应器系统在模型和控制方面具有较高的研究价值。
随着数字电路技术的发展和进步,数字控制技术和微型处理器在工程实践和科学研究等领域被广泛应用。采用离散时间采样控制器来控制连续时间对象,这样的控制系统被称为计算机控制系统或采样控制系统。目前,计算机控制技术和数字技术已广泛应用于工业控制系统中,针对采样控制系统的研究具有重要的理论价值。
目前,国内已基本上实现了连续搅拌釜式反应器系统的理论研究,但对连续搅拌釜式反应器系统的控制策略研究不够深入,尤其是采样控制非线性控制策略的研究还存在空白。
发明内容
本发明所要解决的技术问题在于提供一种方法较为简单的连续搅拌釜式反应器系统的采样控制方法,使系统在不同运行条件下具有更好的性能与鲁棒性。
实现本发明目的的技术解决方案为:一种用于连续搅拌釜式反应器系统的采样控制方法,包括以下步骤:
步骤1、根据连续搅拌釜式反应器系统结构图,建立连续搅拌釜式反应器系统的模型,然后建立连续搅拌釜式反应器系统的状态空间模型;
步骤2、根据步骤1建立的连续搅拌釜式反应器系统的状态空间模型,构造当一个采样周期内系统不发生切换时的李雅普诺夫函数,并设置相应的控制率;
步骤3、根据步骤1建立的连续搅拌釜式反应器系统的状态空间模型,构造当一个采样周期内系统发生切换时的李雅普诺夫函数,并设置相应的控制率;
步骤4、在步骤2和步骤3的基础上,设置采样周期与平均驻留时间以使构造的李雅普诺夫函数有界,即完成连续搅拌釜式反应器系统的采样控制。
本发明与现有技术相比,其显著优点在于:1)本发明的方法通过设置控制增益与采样周期,设计过程灵活,可根据不同情况构建出更合适的控制器;2)本发明中采样控制作为一种新型的控制策略,控制效果很显著;3)本发明控制方法可以使控制成本显著降低,且更加便于用微处理器实现,具有较高的工程实用价值。
附图说明
图1为本发明用于连续搅拌釜式反应器系统的采样控制方法的流程图。
图2为本发明连续搅拌釜式反应器系统结构图。
具体实施方式
结合图1、图2,本发明一种用于连续搅拌釜式反应器系统的采样控制方法,包括以下步骤:
步骤1、根据连续搅拌釜式反应器系统结构图,建立连续搅拌釜式反应器系统的模型,然后建立连续搅拌釜式反应器系统的状态空间模型;
步骤2、根据步骤1建立的连续搅拌釜式反应器系统的状态空间模型,构造当一个采样周期内系统不发生切换时的李雅普诺夫函数,并设置相应的控制率;
步骤3、根据步骤1建立的连续搅拌釜式反应器系统的状态空间模型,构造当一个采样周期内系统发生切换时的李雅普诺夫函数,并设置相应的控制率;
步骤4、在步骤2和步骤3的基础上,设置采样周期与平均驻留时间以使构造的李雅普诺夫函数有界,即完成连续搅拌釜式反应器系统的采样控制。
进一步地,步骤1所述根据连续搅拌釜式反应器系统结构图,建立连续搅拌釜式反应器系统的模型,然后建立连续搅拌釜式反应器系统的状态空间模型,具体为:
步骤1-1、根据连续搅拌釜式反应器系统结构图,建立连续搅拌釜式反应器系统的模型为:
式中,CA为反应物A浓度,为反应物入口浓度,为反应物入口温度,T1为反应器温度,V为反应器的体积,Vρ为流体的体积,CP为流体的热容量,R为气体常数,E为活化能,qσ为进料流速,UA为传热常数,Tc为冷却液温度,a0为反应速率常数,ΔH为反应焓,f为饲料流指数,ρ为流体密度,为“*”的导数;
步骤1-2、由步骤1-1的模型建立连续搅拌釜式反应器系统的状态空间模型为:
式中,x1,x2均代表系统状态,h1,σ(t),h2,σ(t)代表控制增益,f1,σ(t),f2,σ(t)代表未知连续非线性函数,σ(t)为切换信号,u代表控制器输入,表示为如下形式:
式中,tk为采样点,k为正整数。
进一步地,步骤2所述根据步骤1建立的连续搅拌釜式反应器系统的状态空间模型,构造一个当采样周期内系统不发生切换时的李雅普诺夫函数,并设置相应的控制率,具体为:
当一个采样周期内系统不发生切换时,即当切换信号满足如下形式时:
σ(tk)=σ(tk+1)=l (4)
式中,l为切换的第l个子系统;
步骤2-1、引入坐标变换:
式中,z1,z2,l代表坐标变换后的系统状态,α1,l为虚拟控制器;
步骤2-2、根据步骤2-1的公式(5)构造李雅普诺夫函数V1为:
步骤2-3、对步骤2-2的公式(6)进行求导放缩得为:
式中, h 1,l分别为h1,l的最大值、最小值,c1,l,c2,l均为正常数,h1,l、h2,l为控制增益;
步骤2-4、设置相应的控制率u(t):
式中,
通过杨氏不等式变换与式(8),式(7)变换为:
式中,d,α均为正常数;
令则式(9)变换为:
由此可得:
Vn,l(tk+1)≤δVn,l(tk) (11)
式中,δ=(e-T+d4(1-e-T))2,其中d4为可调节变量以使连续搅拌釜式反应器系统稳定,Vn,l(*)代表李雅普诺夫函数。
进一步地,步骤3所述根据步骤1建立的连续搅拌釜式反应器系统的状态空间模型,构造当一个采样周期内系统发生切换时的李雅普诺夫函数,并设置相应的控制率,具体为:
当一个采样周期内系统发生切换时,即当切换信号满足如下形式时:
σ(tk)=l,σ(tk+1)=m (12)
步骤3-1、结合步骤1中构建的连续搅拌釜式反应器系统的状态空间模型与步骤2-1 中的坐标变换,构造一个保证连续搅拌釜式反应器系统稳定性的李雅普诺夫函数V2,m:
步骤3-2、对步骤3-1中的公式(13)进行求导放缩获得为:
式中,z1,m,z2,m为系统状态;
步骤3-3、设置相应的控制率
式中,γ为正常数;
步骤3-4、通过杨氏不等式变换与式(15),式(14)变换为:
式中,T为采样周期;
令则式(16)变换为:
式中,d3为可调节变量以使连续搅拌釜式反应器系统稳定;
令V2,m=ZTPmZ,其中Z=[x1,x2]T,Pm为正定矩阵;
由此可得:
V2,m(tk+1)<υV2,l(tk) (18)
式中, maxm∈2λ(Pm),minm∈2λ(Pm)分别为矩阵Pm的最大特征值和最小特征值,V2,m(*)、V2,l(*)均代表李雅普诺夫函数。
进一步地,步骤4所述在步骤2和步骤3的基础上,设置采样周期与平均驻留时间以使构造的李雅普诺夫函数有界,即完成连续搅拌釜式反应器系统的采样控制,具体为:
步骤4-1、设置采样周期T为:
设置平均滞留时间Ta为:
设置李雅普诺夫函数Vσ(t)为:
式中,zi,σ(t)为系统状态;
步骤4-2、由步骤2中的式(11)和步骤3中的式(18),获得:
式中,N0为抖振界,t0为时间初值;
由式(22)可知,构造的李雅普诺夫函数有界。
本发明的采样控制方法使连续搅拌釜式反应器系统在不同运行条件下具有更好的性能与鲁棒性;此外,本发明的方法简单,可以使控制成本显著降低,且更加便于用微处理器实现,具有较高的工程实用价值。
Claims (5)
1.一种用于连续搅拌釜式反应器系统的采样控制方法,其特征在于,包括以下步骤:
步骤1、根据连续搅拌釜式反应器系统结构图,建立连续搅拌釜式反应器系统的模型,然后建立连续搅拌釜式反应器系统的状态空间模型;
步骤2、根据步骤1建立的连续搅拌釜式反应器系统的状态空间模型,构造当一个采样周期内系统不发生切换时的李雅普诺夫函数,并设置相应的控制率;
步骤3、根据步骤1建立的连续搅拌釜式反应器系统的状态空间模型,构造当一个采样周期内系统发生切换时的李雅普诺夫函数,并设置相应的控制率;
步骤4、在步骤2和步骤3的基础上,设置采样周期与平均驻留时间以使构造的李雅普诺夫函数有界,即完成连续搅拌釜式反应器系统的采样控制。
2.根据权利要求1所述的用于连续搅拌釜式反应器系统的采样控制方法,其特征在于,步骤1所述根据连续搅拌釜式反应器系统结构图,建立连续搅拌釜式反应器系统的模型,然后建立连续搅拌釜式反应器系统的状态空间模型,具体为:
步骤1-1、根据连续搅拌釜式反应器系统结构图,建立连续搅拌釜式反应器系统的模型为:
式中,CA为反应物A浓度,为反应物入口浓度,为反应物入口温度,T1为反应器温度,V为反应器的体积,Vρ为流体的体积,CP为流体的热容量,R为气体常数,E为活化能,qσ为进料流速,UA为传热常数,Tc为冷却液温度,a0为反应速率常数,ΔH为反应焓,f为饲料流指数,ρ为流体密度,为“*”的导数;
步骤1-2、由步骤1-1的模型建立连续搅拌釜式反应器系统的状态空间模型为:
式中,x1,x2均代表系统状态,h1,σ(t),h2,σ(t)代表控制增益,f1,σ(t),f2,σ(t)代表未知连续非线性函数,σ(t)为切换信号,u代表控制器输入,表示为如下形式:
式中,tk为采样点,k为正整数。
3.根据权利要求2所述的用于连续搅拌釜式反应器系统的采样控制方法,其特征在于,步骤2所述根据步骤1建立的连续搅拌釜式反应器系统的状态空间模型,构造一个当采样周期内系统不发生切换时的李雅普诺夫函数,并设置相应的控制率,具体为:
当一个采样周期内系统不发生切换时,即当切换信号满足如下形式时:
σ(tk)=σ(tk+1)=l (4)
式中,l为切换的第l个子系统;
步骤2-1、引入坐标变换:
式中,z1,z2,l代表坐标变换后的系统状态,α1,l为虚拟控制器;
步骤2-2、根据步骤2-1的公式(5)构造李雅普诺夫函数V1为:
步骤2-3、对步骤2-2的公式(6)进行求导放缩得为:
式中, h 1,l分别为h1,l的最大值、最小值,c1,l,c2,l均为正常数,h1,l、h2,l为控制增益;
步骤2-4、设置相应的控制率u(t):
式中,
通过杨氏不等式变换与式(8),式(7)变换为:
式中,d,α均为正常数;
令则式(9)变换为:
由此可得:
Vn,l(tk+1)≤δVn,l(tk) (11)
式中,δ=(e-T+d4(1-e-T))2,其中d4为可调节变量以使连续搅拌釜式反应器系统稳定,Vn,l(*)代表李雅普诺夫函数。
4.根据权利要求3所述的用于连续搅拌釜式反应器系统的采样控制方法,其特征在于,步骤3所述根据步骤1建立的连续搅拌釜式反应器系统的状态空间模型,构造当一个采样周期内系统发生切换时的李雅普诺夫函数,并设置相应的控制率,具体为:
当一个采样周期内系统发生切换时,即当切换信号满足如下形式时:
σ(tk)=l,σ(tk+1)=m (12)
步骤3-1、结合步骤1中构建的连续搅拌釜式反应器系统的状态空间模型与步骤2-1中的坐标变换,构造一个保证连续搅拌釜式反应器系统稳定性的李雅普诺夫函数V2,m:
步骤3-2、对步骤3-1中的公式(13)进行求导放缩获得为:
式中,z1,m,z2,m为系统状态;
步骤3-3、设置相应的控制率
式中,γ为正常数;
步骤3-4、通过杨氏不等式变换与式(15),式(14)变换为:
式中,T为采样周期;
令则式(16)变换为:
式中,d3为可调节变量以使连续搅拌釜式反应器系统稳定;
令V2,m=ZTPmZ,其中Z=[x1,x2]T,Pm为正定矩阵;
由此可得:
V2,m(tk+1)<υV2,l(tk) (18)
式中,maxm∈2λ(Pm),minm∈2λ(Pm)分别为矩阵Pm的最大特征值和最小特征值,V2,m(*)、V2,l(*)均代表李雅普诺夫函数。
5.根据权利要求4所述的用于连续搅拌釜式反应器系统的采样控制方法,其特征在于,步骤4所述在步骤2和步骤3的基础上,设置采样周期与平均驻留时间以使构造的李雅普诺夫函数有界,即完成连续搅拌釜式反应器系统的采样控制,具体为:
步骤4-1、设置采样周期T为:
设置平均滞留时间Ta为:
设置李雅普诺夫函数Vσ(t)为:
式中,zi,σ(t)为系统状态;
步骤4-2、由步骤2中的式(11)和步骤3中的式(18),获得:
式中,N0为抖振界,t0为时间初值;
由式(22)可知,构造的李雅普诺夫函数有界。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910146505.1A CN109870908A (zh) | 2019-02-27 | 2019-02-27 | 一种用于连续搅拌釜式反应器系统的采样控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910146505.1A CN109870908A (zh) | 2019-02-27 | 2019-02-27 | 一种用于连续搅拌釜式反应器系统的采样控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109870908A true CN109870908A (zh) | 2019-06-11 |
Family
ID=66919324
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910146505.1A Pending CN109870908A (zh) | 2019-02-27 | 2019-02-27 | 一种用于连续搅拌釜式反应器系统的采样控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109870908A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112731811A (zh) * | 2020-12-28 | 2021-04-30 | 南京理工大学 | 一种针对双级化学反应器的采样控制方法 |
CN114019802A (zh) * | 2021-11-08 | 2022-02-08 | 河海大学 | 一种用于两级化学反应器系统的采样控制方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090143871A1 (en) * | 2002-04-18 | 2009-06-04 | Cleveland State University | Controllers, observers, and applications thereof |
CN102023574A (zh) * | 2010-12-31 | 2011-04-20 | 浙江大学 | 一阶反应连续搅拌釜式反应器的混杂模型优化控制方法 |
CN103217904A (zh) * | 2013-04-15 | 2013-07-24 | 浙江工业大学 | 一种具有通信受限的广域测量电力系统控制方法 |
CN105242540A (zh) * | 2015-10-27 | 2016-01-13 | 西安石油大学 | 基于平均驻留时间的网络控制系统动态切换控制方法 |
CN105487515A (zh) * | 2015-12-29 | 2016-04-13 | 浙江工业大学 | 一种连续搅拌釜式反应过程工艺设计与控制的集成优化方法 |
-
2019
- 2019-02-27 CN CN201910146505.1A patent/CN109870908A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090143871A1 (en) * | 2002-04-18 | 2009-06-04 | Cleveland State University | Controllers, observers, and applications thereof |
CN102023574A (zh) * | 2010-12-31 | 2011-04-20 | 浙江大学 | 一阶反应连续搅拌釜式反应器的混杂模型优化控制方法 |
CN103217904A (zh) * | 2013-04-15 | 2013-07-24 | 浙江工业大学 | 一种具有通信受限的广域测量电力系统控制方法 |
CN105242540A (zh) * | 2015-10-27 | 2016-01-13 | 西安石油大学 | 基于平均驻留时间的网络控制系统动态切换控制方法 |
CN105487515A (zh) * | 2015-12-29 | 2016-04-13 | 浙江工业大学 | 一种连续搅拌釜式反应过程工艺设计与控制的集成优化方法 |
Non-Patent Citations (1)
Title |
---|
SHI LI等: "Global Stabilization of a Class of Switched Nonlinear Systems Under Sampled-Data Control", 《IEEE TRANSACTIONS ON SYSTEMS, MAN, AND CYBERNETICS》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112731811A (zh) * | 2020-12-28 | 2021-04-30 | 南京理工大学 | 一种针对双级化学反应器的采样控制方法 |
CN114019802A (zh) * | 2021-11-08 | 2022-02-08 | 河海大学 | 一种用于两级化学反应器系统的采样控制方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109870908A (zh) | 一种用于连续搅拌釜式反应器系统的采样控制方法 | |
CN104650104B (zh) | 锌离子-卟啉纳米络合物的制备方法 | |
CN102875493B (zh) | 在连续流微通道反应器中制备环氧脂肪酸甲酯的方法 | |
CN203890063U (zh) | 磷酸铁锂连续制备装置 | |
Moroi | Distribution of solubilizates among micelles and kinetics of micelle-catalyzed reactions | |
CN202928343U (zh) | 钟罩炉负压烧结装置 | |
Sue et al. | Continuous hydrothermal synthesis of Fe2O3 nanoparticles using a central collision-type micromixer for rapid and homogeneous nucleation at 673 K and 30 MPa | |
CN105273113A (zh) | 一种聚丙烯酸合成的反应装置及聚丙烯酸合成方法 | |
Niu et al. | Optimal control of nonlinear switched system in an uncoupled microbial fed-batch fermentation process | |
CN102176178A (zh) | 超临界水处理系统中电加热器出口温度的控制方法 | |
CN102070332A (zh) | 一种兼有高起始磁导率和高热磁灵敏度的纳米晶锰锌铁氧体磁性材料及其制备方法 | |
Shen et al. | Nonlinear state-dependent impulsive system and its parameter identification in microbial fed-batch culture | |
CN108023087A (zh) | 一种SiO2包覆的镍锰酸锂正极材料的制备方法 | |
CN104261477A (zh) | 一种Mn3O4八面体结构的制备方法 | |
CN114019802A (zh) | 一种用于两级化学反应器系统的采样控制方法及系统 | |
CN201470381U (zh) | 聚合釜升温恒温反应控制装置 | |
CN105511269B (zh) | 一种基因调节非线性动力系统控制模型 | |
Liu et al. | Modeling and optimal control of a nonlinear dynamical system in microbial fed-batch fermentation | |
Chang | The snowball effect in fed‐batch bioreactions | |
CN105600814B (zh) | 一种制备花状结构Cu2O光电材料的方法 | |
CN110407515B (zh) | 一种耐热压力传感器材料及其制备方法 | |
CN104788090A (zh) | 一种高性能LiZn铁氧体材料 | |
CN111704157A (zh) | 一种纳米氧化锌的微通道制备方法 | |
CN104829525A (zh) | 一种连续化反应生产2-氯-3-氰基吡啶的方法 | |
Olsen et al. | Modeling and simulation of single cell protein production |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190611 |