CN106340331B - 一种用于核反应堆功率的自抗扰控制方法 - Google Patents
一种用于核反应堆功率的自抗扰控制方法 Download PDFInfo
- Publication number
- CN106340331B CN106340331B CN201610827618.4A CN201610827618A CN106340331B CN 106340331 B CN106340331 B CN 106340331B CN 201610827618 A CN201610827618 A CN 201610827618A CN 106340331 B CN106340331 B CN 106340331B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- mtr
- mfrac
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21D—NUCLEAR POWER PLANT
- G21D3/00—Control of nuclear power plant
- G21D3/001—Computer implemented control
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21D—NUCLEAR POWER PLANT
- G21D3/00—Control of nuclear power plant
- G21D3/001—Computer implemented control
- G21D3/002—Core design; core simulations; core optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Plasma & Fusion (AREA)
- High Energy & Nuclear Physics (AREA)
- Feedback Control In General (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明提供一种用于核反应堆功率的自抗扰控制方法,目的在于针对目前核反应堆对象没有既简单又能取得良好控制效果的控制器,提出将自抗扰控制器应用于核反应堆功率控制中的方法。通过推导将原有模型转化为适合设计自抗扰控制器的非线性模型;并且充分利用模型信息,减轻ESO的负担;还利用系统输出y可测的特点,设计基于RESO的自抗扰控制器;最后通过简单调整带宽、就能容易的得出控制器参数,此控制系统相比预测控制、模糊控制等先进控制算法结构简单,但是有更好的控制效果。
Description
技术领域
本发明属于核电站的控制领域,具体涉及一种用于核反应堆功率的自抗扰控制方法。
背景技术
核电机组是高度复杂的非线性系统,其参数是运行功率、核燃料燃尽程度和控制棒价值的函数,并随时间变化。在负荷跟随条件下,当出现大的功率变动时,就必须特别考虑这些因素的影响。现有的大部分反应堆的控制采用常规调节系统,按照基本负荷工作点参数设计。压水堆通常所采用的三通道非线性控制器实际为带非线性增益补偿的PID控制器。然而,常规控制器的调节性能,在大负荷变动条件下受到挑战。许多其他先进的过程控制方法也在不断出现,比如基于Takagi-Sugeno模糊模型的核反应堆功率积分控制系统,仿真结果表明具有较好的跟踪特性,能够在很小的超调和很少的振荡情况下实现零稳态跟踪,而且满足核反应堆运行安全要求。免疫P-PID串级控制的核反应堆功率调节器,采用基于人工免疫机理的免疫P和PID串级控制策略的核反应堆功率调节器是可行的,而且基于人工免疫的功率调节器比传统PID功率调节器的控制效果要好。广义预测自校正控制算法,在核反应堆功率控制中,结果表明能够较好的控制反应堆功率的输出。另外,自适应控制、鲁棒控制等方法主要是针对处理系统内部不确定性而提出的。这些控制算法都取得了一定的成效,但是由于鲁棒控制的结果相对保守、预测控制方法的算法繁琐、自适应控制难以在大时间延迟系统中得到应用等缺点,它们目前尚未得到广泛工程应用。因此,研究一种既简单又不完全依赖系统模型,且鲁棒性强的控制策略,对于提高现有反应堆功率控制系统的性能有较大的实际意义。
实际控制系统中普遍存在着内部不确定性(参数和未建模动态)和外部不确定性(扰动),因此不确定性系统的控制向来是一个富有挑战性的基本问题。自抗扰控制(ActiveDisturbance Rejection Control,简称ADRC)是中国科学院韩京清教授提出的一种新型的反馈线性化控制策略。自抗扰控制具有较强的抗扰动能力,并且具有较好的性能鲁棒性。此外ADRC还继承了PID控制器结构简单、不依赖数学模型等优点。ADRC是针对同时具有内部和外部不确定性的非线性不确定系统的控制问题而提出的,其核心思想是将系统的内部不确定性(定常或时变,线性或非线性)和外部不确定性(外部扰动)一起作为“总扰动”,通过构造“扩张状态观测器(ESO)”对其进行估计并由控制率来实时补偿,以期获得较强的控制不确定性的能力以及较好的控制精度。ADRC应用非线性控制策略,需整定的参数较多,整定往往依赖于设计者的经验,实践中较难应用。高志强教授等在ADRC控制器的基础上,对其各结构进行线性化设计,提出了线性自抗扰控制(Linear Active Disturbance RejectionControl,简称LADRC),并引入了带宽这个物理意义明确的整定参数,使得控制器参数大量减少,将控制器参数简化为控制器带宽和观测器带宽的函数,提出了一套简单的参数整定方法。线性化和带宽概念的引入给理论研究提供了全新的视角,同时降低了难度。而且LADRC也便于实际应用,并已有大量的应用研究表明LADRC依然对复杂的非线性不确定对象有很强的控制能力。ADRC发展至今,不管在理论还是实践中都走向一个新的台阶。赵志良证明了ADRC在扩张状态有界时的收敛性,Yang等给出了ADRC可以估计的扰动范围。LADRC在许多领域应用的综述和范例,包括精密车床中快速刀具伺服控制、微机电传感器、时滞系统、航天器姿态控制、电力系统中锅炉燃烧和负载频率控制等等,均显示了自抗扰控制技术的巨大潜力。在发电领域,LADRC也在火电单元机组协调系统、火电厂主汽温控制系统、循环流化床锅炉燃烧系统等仿真控制中获得良好的控制效果。但尚无将利用模型信息且基于带宽参数调整的LADRC应用于核反应堆功率控制的研究。
发明内容
为了解决上述问题,本发明提供一种用于核反应堆功率的自抗扰控制方法,所述自抗扰控制方法通过推导核反应堆功率模型,将原来5阶的非线性模型转化为2阶的非线性模型,利用非线性模型中的信息,减轻ESO的负担,设计基于ESO的自抗扰控制器,并且利用系统输出y可测的特点,设计基于RESO的自抗扰控制器,通过自抗扰控制器的设计原理,给定ωc、ωo、nr0,确定该反应堆功率控制系统,完成自抗扰控制;
进一步地,所述自抗扰控制器的设计步骤为:
S1:提供核反应堆功率原始非线性模型;
S2:通过对S1中的原始非线性模型进行模型变换获得用于设计自抗扰控制器的2阶非线性模型;
S3:利用S2中2阶非线性模型信息,减轻ESO的负担,设计基于ESO的自抗扰控制器;
S4:通过降阶状态观测器和变量代换推导,避免RESO的输入是y的导数项,设计基于RESO的自抗扰控制器;
S5:通过自抗扰控制器的设计原理,给定ωc、ωo、nr0,确定该反应堆功率控制系统,自抗扰控制过程;
进一步地,所述步骤S1具体包括:
1-1)假定反应堆内各点的中子密度随时间的变化特性与空间位置不相关的,把核反应堆近似地看做一个没有空间度量的“点”,得到了如下点堆动态方程(1)、(2);
1-2)根据宏观能量守恒定律,得到压水反应堆的关于燃料平均温度和冷却剂出口温度的两个方程(3)、(4);
1-3)控制棒的反应性方程(5);
1-4)根据核反应堆中的温度反馈,所述温度反馈为燃料和冷却剂的反应性温度反馈系数,得到总体的反应性方程(6);
其中,(2)式中,ρ为反应性;β为缓发中子份额;Λ为一代中子寿命,s;λ缓发中子先驱核衰变常数,s-1;nr相对中子密度;cr相对先驱密度;Tf为燃料平均温度,℃;ff为储存在燃料中的反应堆功率份额;p0为初始平衡状态功率,MW;uf为燃料热容量;uc为冷却剂的热容量;Ω为燃料和冷却剂之间传热系数;Te是冷却剂进入反应堆时的温度;Tl冷却剂离开反应堆时的平均温度,℃;M为质量流量与水的热容量的乘积,MW/℃;ρr控制棒引入的反应性;Gr单位长度控制棒反应性价值;Zr控制输入,即控制棒速度;αf燃料温度反应性系数;Tf0反应堆燃料平均温度,℃;αc冷却剂温度反应性系数;Tl0平衡时冷却剂离开反应堆的温度,℃;
进一步地,其中,
M(nr0)=(28nr0+74)MW/℃;
在控制系统设计时,其中的常数取为
β=0.006019,Λ=0.0001s,λ=0.15s-1,ff=0.92;uf=26.3MW.s/℃;Gr=0.0145;
进一步地,所述步骤S2具体为对(1)求导,将(1)(2)代入可得:
设初始条件为:
nr(0)=1,cr(0)=1,ρ(0)=0,ρr(0)=0,Tf(0)=Tf0,Tl(0)=Tl0
设:
nr=nr(0)+δnr (8);
cr=cr(0)+δcr (9);
Tf=Tf(0)+δTf (10);
Tl=Tl(0)+δTl (11);
ρr=ρr(0)+δρr=δρr (12);
δρ=ρ=δρr+δρf+δρl (13);
其中,
δρf=αfδTf (14);
δρl=αlδTl (15);
对(13)求导,并将(5)、(14)、(15)代入可得
由(8)可得:
将(5)、(8)、(13)、(17)代入(7)可得
假设冷却剂入口温度Te不变,由稳态初始条件和(3)可得
将(8)、(10)、(11)、(19)代入(3)可得
其中:
同理,由(4)、(8)、(10)、(11)可得
其中:
将(20)、(21)代入(18)可得
其中:
k1=αfa31+αla41,k2=αfa33+αla43,k3=αfa34+αla44;
进一步地,所述S3具体为:对于二阶对象
y和u分别是输出和输入,g代表模型不确定项,w是外部扰动;所述(24)中参数a1,a2,b为已知;令f=g+w为总扰动,所述(24)的状态方程为:
y=Cx;
其中,x3=f是扩张状态,则线性扩张状态观测器(LESO)为:
其中L是观测器增益,能够用极点配置方法获得
L=[l1 l2 l3]T;
LESO:
A*=A-LC,B*=[B L],C*=I3,D*=(0)3×2;
控制器设计:
PD控制器控制:
u0=kp(r-z1)-kdz2;
r为设定值,闭环传递函数为:
则增益取为:
进一步地,所述S4具体为:
RESO:
令:
V=2ωoy-z1 W=ωo 2y-z2
由(1)得:
由(2)得:
因此:
进一步地,所述S5中具体为当:
nr0=1;
ωc和ωo需要整定,
在实际整定过程中遵循如下规律:
1)ωo越大,ESO的观测能力越强,观测器对噪声的敏感性增加;
ωo从较小的值逐渐增大,直至观测精度满足要求为止;
2)ωc越大,控制作用越强,系统的响应速度越快,但超调和振荡会越严重,同时稳定性下降;
ωo和ωc在保证控制器稳定的前提下通过多次试验和综合比较误差指标加以确定,观测器及控制器的误差上界与其带宽ωo和ωc成反比,带宽越宽,误差越小,稳定性越低;
本发明的有益效果如下:
1)对于本发明专利使用的模型当ωo=16、ωc=28时能够取得良好的控制效果,优于模糊或者预测等一些设计,满足核反应堆控制系统设计的一般要求;
2)在满足运行要求的前提下,控制系统简单可靠;
3)控制系统运行参数的瞬态变化量小,稳态运行参数更接近设计给定值,提高核电厂的输出功率;
4)当控制系统的设计考虑到最坏的工作条件时,在各种条件下,系统都有一定的稳定裕度,不大的超调量和合理的调整时间,控制器有足够强的鲁棒性维持系统稳定;
5)系统处于稳定工作点时,控制器能够有效地抑制系统中的各种扰动。
附图说明
图1为本发明所述方法中基于RESO的自抗扰控制器应用于核反应堆功率控制系统中的结构图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细描述。应当理解,此处所描述的具体实施例仅仅用于解释本发明,并不用于限定本发明。相反,本发明涵盖任何由权利要求定义的在本发明的精髓和范围上做的替代、修改、等效方法以及方案。进一步,为了使公众对本发明有更好的了解,在下文对本发明的细节描述中,详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。
下面结合附图和具体实施例对本发明作进一步说明,但不作为对本发明的限定。下面为本发明的举出最佳实施例:
本发明提供一种用于核反应堆功率的自抗扰控制方法,所述自抗扰控制方法通过推导核反应堆功率模型,将原来5阶的非线性模型转化为2阶的非线性模型,利用非线性模型中的信息,减轻ESO的负担,设计基于ESO的自抗扰控制器,并且利用系统输出y可测的特点,设计基于RESO的自抗扰控制器,通过自抗扰控制器的设计原理,给定ωc、ωo、nr0,确定该反应堆功率控制系统,完成自抗扰控制,如图1所示。
所述自抗扰控制器的设计步骤为:
S1:提供核反应堆功率原始非线性模型;
S2:通过对S1中的原始非线性模型进行模型变换获得用于设计自抗扰控制器的2阶非线性模型;
S3:利用S2中2阶非线性模型信息,减轻ESO的负担,设计基于ESO的自抗扰控制器;
S4:通过降阶状态观测器和变量代换推导,避免RESO的输入是y的导数项,设计基于RESO的自抗扰控制器;
S5:通过自抗扰控制器的设计原理,给定ωc、ωo、nr0,确定该反应堆功率控制系统,自抗扰控制过程。
所述步骤S1具体包括:
1-1)假定反应堆内各点的中子密度随时间的变化特性与空间位置不相关的,把核反应堆近似地看做一个没有空间度量的“点”,得到了如下点堆动态方程(1)、(2);
1-2)根据宏观能量守恒定律,得到压水反应堆的关于燃料平均温度和冷却剂出口温度的两个方程(3)、(4);
1-3)控制棒的反应性方程(5);
1-4)根据核反应堆中的温度反馈,所述温度反馈为燃料和冷却剂的反应性温度反馈系数,得到总体的反应性方程(6);
其中,(2)式中,ρ为反应性;β为缓发中子份额;Λ为一代中子寿命,s;λ缓发中子先驱核衰变常数,s-1;nr相对中子密度;cr相对先驱密度;Tf为燃料平均温度,℃;ff为储存在燃料中的反应堆功率份额;p0为初始平衡状态功率,MW;uf为燃料热容量;uc为冷却剂的热容量;Ω为燃料和冷却剂之间传热系数;Te是冷却剂进入反应堆时的温度;Tl冷却剂离开反应堆时的平均温度,℃;M为质量流量与水的热容量的乘积,MW/℃;ρr控制棒引入的反应性;Gr单位长度控制棒反应性价值;Zr控制输入,即控制棒速度;αf燃料温度反应性系数;Tf0反应堆燃料平均温度,℃;αc冷却剂温度反应性系数;Tl0平衡时冷却剂离开反应堆的温度,℃。
其中,
M(nr0)=(28nr0+74)MW/℃;
在控制系统设计时,其中的常数取为β=0.006019,Λ=0.0001s,λ=0.15s-1,ff=0.92;uf=26.3MW.s/℃;Gr=0.0145。
所述步骤S2具体为对(1)求导,将(1)(2)代入可得
设初始条件为:
nr(0)=1,cr(0)=1,ρ(0)=0,ρr(0)=0,Tf(0)=Tf0,Tl(0)=Tl0
设:
nr=nr(0)+δnr (8);
cr=cr(0)+δcr (9);
Tf=Tf(0)+δTf (10);
Tl=Tl(0)+δTl (11);
ρr=ρr(0)+δρr=δρr (12);
δρ=ρ=δρr+δρf+δρl (13);
其中,
δρf=αfδTf (14);
δρl=αlδTl (15);
对(13)求导,并将(5)、(14)、(15)代入可得
由(8)可得:
将(5)、(8)、(13)、(17)代入(7)可得
假设冷却剂入口温度Te不变,由稳态初始条件和(3)可得
将(8)、(10)、(11)、(19)代入(3)可得
其中:
同理,由(4)、(8)、(10)、(11)可得
其中:
将(20)、(21)代入(18)可得
其中:
k1=αfa31+αla41,k2=αfa33+αla43,k3=αfa34+αla44;
所述S3具体为:对于二阶对象
y和u分别是输出和输入,g代表模型不确定项,w是外部扰动;所述(24)中参数a1,a2,b为已知;令f=g+w为总扰动,所述(24)的状态方程为:
y=Cx;
其中,x3=f是扩张状态,则线性扩张状态观测器(LESO)为:
其中L是观测器增益,能够用极点配置方法获得
L=[l1 l2 l3]T;
LESO:
A*=A-LC,B*=[B L],C*=I3,D*=(0)3×2;
控制器设计:
PD控制器控制:
u0=kp(r-z1)-kdz2;
r为设定值,闭环传递函数为:
则增益取为:
kp=ωc 2,kd=2ωc。
所述S4具体为:
RESO:
令:
V=2ωoy-z1 W=ωo 2y-z2
由(1)得:
由(2)得:
因此:
所述S5中具体为当:
nr0=1;
ωc和ωo需要整定,
在实际整定过程中遵循如下规律:
1)ωo越大,ESO的观测能力越强,观测器对噪声的敏感性增加;ωo从较小的值逐渐增大,直至观测精度满足要求为止;
2)ωc越大,控制作用越强,系统的响应速度越快,但超调和振荡会越严重,同时稳定性下降;
ωo和ωc在保证控制器稳定的前提下通过多次试验和综合比较误差指标加以确定,观测器及控制器的误差上界与其带宽ωo和ωc成反比,带宽越宽,误差越小,稳定性越低。
本发明的目的在于针对目前核反应堆对象没有既简单又能取得良好控制效果的控制器,提出将自抗扰控制器应用于核反应堆功率控制中的方法。通过推导将原有模型转化为适合设计自抗扰控制器的非线性模型;并且充分利用模型信息,减轻ESO的负担;还利用系统输出y可测的特点,设计基于RESO的自抗扰控制器;最后形成给定nr0通过简单调整带宽ωc、ωo就能容易的得出控制器参数,此控制系统相比预测控制、模糊控制等先进控制算法结构简单,但是有更好的控制效果。因此本发明对于将自抗扰控制器将来应用于核电系统具有重要的指导意义。
1.核反应堆功率模型转换
首次将核反应堆功率模型转换成适合设计自抗扰控制器的模型,之前对反应堆功率控制多数将其线性化,而本发明专利是在非线性模型基础上设计线性自抗扰控制器。
1.1原始模型
假定反应堆内各点的中子密度随时间的变化特性与空间位置是不相关的,这样可以把核反应堆近似地看做一个没有空间度量的“点”,就得到了如下点堆动态方程(1)、(2);根据宏观能量守恒定律,可以得到压水反应堆的关于燃料平均温度和冷却剂出口温度的两个方程(3)、(4);控制棒的反应性方程(5);考虑到核反应堆中的温度反馈,即燃料和冷却剂的反应性温度反馈系数,可以得到总体的反应性方程(6)。
(2)式中,ρ为反应性;β为缓发中子份额;Λ为一代中子寿命,s;λ缓发中子先驱核衰变常数,s-1;nr相对中子密度;cr相对先驱密度;Tf为燃料平均温度,℃;ff为储存在燃料中的反应堆功率份额;p0为初始平衡状态功率,MW;uf为燃料热容量;uc为冷却剂的热容量;Ω为燃料和冷却剂之间传热系数;Te是冷却剂进入反应堆时的温度;Tl冷却剂离开反应堆时的平均温度,℃;M为质量流量与水的热容量的乘积,MW/℃;ρr控制棒引入的反应性;Gr单位长度控制棒反应性价值;Zr控制输入,即控制棒速度;αf燃料温度反应性系数;Tf0反应堆燃料平均温度,℃;αc冷却剂温度反应性系数;Tl0平衡时冷却剂离开反应堆的温度,℃。其中,
M(nr0)=(28nr0+74)MW/℃;
在控制系统设计时,其中的常数取为β=0.006019,Λ=0.0001s,λ=0.15s-1,ff=0.92;uf=26.3MW.s/℃;Gr=0.0145。
1.2模型变换
对(1)求导,将(1)(2)代入可得
设初始条件为
nr(0)=1,cr(0)=1,ρ(0)=0,ρr(0)=0,Tf(0)=Tf0,Tl(0)=Tl0
设
nr=nr(0)+δnr (8);
cr=cr(0)+δcr (9);
Tf=Tf(0)+δTf (10);
Tl=Tl(0)+δTl (11);
ρr=ρr(0)+δρr=δρr (12);
δρ=ρ=δρr+δρf+δρl (13);
其中
δρf=αfδTf (14);
δρl=αlδTl (15);
(13)求导,并将(5)、(14)、(15)代入可得
由(8)可得
将(5)、(8)、(13)、(17)代入(7)可得
假设冷却剂入口温度Te不变,由稳态初始条件和(3)可得
将(8)、(10)、(11)、(19)代入(3)可得
其中
同理,由(4)、(8)、(10)、(11)可得
其中
将(20)、(21)代入(18)可得
其中
k1=αfa31+αla41k2=αfa33+αla43k3=αfa34+αla44
2.线性自抗扰控制器的设计
2.1利用模型信息
对于二阶对象
y和u分别是输出和输入,g代表模型不确定项,w是外部扰动。这里参数a1,a2,b是已知的。令f=g+w为总扰动。对象(24)的状态方程为:
这里x3=f是扩张状态,则线性扩张状态观测器(LESO)为:
其中L是观测器增益,能够用极点配置方法获得L=[l1 l2 l3]T
LESO:
A*=A-LC B*=[B L] C*=I3 D*=(0)3×2
控制器按照下面的方法设计:
很容易用PD控制器控制:
u0=kp(r-z1)-kdz2
这里r是设定值,闭环传递函数为:
则增益取为:kp=ωc 2 kd=2ωc。
2.2利用降阶状态观测器(RESO)
由于本系统输出可测不需要估计,所以可使用RESO,取得了更好的控制效果。也是首次将RESO运用到反应堆功率控制系统中,并且经过变量代换推导,避免了RESO的输入是y的导数项。
RESO:
令V=2ωoy-z1 W=ωo 2y-z2
由(1)得
由(2)得
即
因此
3.核反应堆功率控制系统设计
首次将基于带宽调节的LADRC应用于该核反应堆功率控制系统中,并且在设计过程中充分利用了模型信息,效果明显。并且针对本系统输出可测不需要估计,还设计了基于RESO的自抗扰控制器,均取得满意的控制效果。
根据1、2部分的设计原理,给定ωc、ωo、nr0即可确定该反应堆功率控制系统。当nr0=1,可得
至此,ADRC有2个参数需要整定,分别是ωc,ωo。这2个参数在实际整定过程中遵循如下规律:
1)ωo越大,ESO的观测能力越强,但这会增加观测器对噪声的敏感性。因此,ωo应从较小的值逐渐增大,直至观测精度满足要求为止。
2)ωc越大,控制作用越强,系统的响应速度越快,但超调和振荡会越严重,同时稳定性下降.在通常情况下,ωo和ωc则在保证控制器稳定的前提下通过多次试验和综合比较误差指标加以确定.有文献已证明,观测器及控制器的误差上界与其带宽ωo和ωc成反比,即只要带宽足够宽,误差就会足够小,但同时牺牲了稳定性。
对于本发明专利使用的模型当ωo=16、ωc=28时能够取得良好的控制效果,优于模糊或者预测等一些设计,满足核反应堆控制系统设计的一般要求:
1)在满足运行要求的前提下,应尽量使控制系统简单可靠;
2)控制系统的设计应尽量减少运行参数的瞬态变化量,并使稳态运行参数更接近设计给定值,尽可能增加核电厂的输出功率;
3)控制系统的设计应考虑到最坏的工作条件,即在各种条件下,系统仍有一定的稳定裕度,不大的超调量和合理的调整时间,控制器有足够强的鲁棒性维持系统稳定。
4)系统处于稳定工作点时,控制器能够有效地抑制系统中的各种扰动。
正常运行时功率调节的超调量应小于3%FP.冷却剂平均温度的超调量不应大于2.5℃。
以上所述的实施例,只是本发明较优选的具体实施方式的一种,本领域的技术人员在本发明技术方案范围内进行的通常变化和替换都应包含在本发明的保护范围内。
Claims (6)
1.一种用于核反应堆功率的自抗扰控制方法,其特征在于,所述自抗扰控制方法通过推导核反应堆功率模型,将原来5阶的非线性模型转化为2阶的非线性模型,利用非线性模型中的信息,减轻线性扩张状态观测器(LESO)的负担,设计了基于LESO的自抗扰控制器;并且利用系统输出y可测的特点,设计了基于降阶线性扩张状态观测器(RLESO)的自抗扰控制器,进一步地通过变量代换方法使得观测器的输入不是而是y,通过自抗扰控制器的设计原理,给定参数ωc、ωo、nr0,确定该反应堆功率控制系统,完成自抗扰控制;
所述自抗扰控制器的设计步骤为:
S1:提供核反应堆功率原始非线性模型;
S2:通过对S1中的原始非线性模型进行模型变换获得用于设计自抗扰控制器的2阶非线性模型;
S3:利用S2中2阶非线性模型信息,减轻LESO的负担,设计基于LESO的自抗扰控制器;
S4:通过降阶状态观测器和变量代换推导,避免RLESO的输入是y的导数项,设计基于RLESO的自抗扰控制器;
S5:通过自抗扰控制器的设计原理,给定参数ωc、ωo、nr0,确定该反应堆功率控制系统,完成自抗扰控制。
2.根据权利要求1所述的自抗扰控制方法,其特征在于,所述步骤S1具体包括:
1)假定反应堆内各点的中子密度随时间的变化特性与空间位置不相关的,把核反应堆近似地看作一个没有空间度量的“点”,得到了如下点堆动态方程(1)、(2);
2)根据宏观能量守恒定律,得到压水反应堆的关于燃料平均温度和冷却剂出口温度的两个方程(3)、(4);
3)控制棒的反应性方程(5);
4)根据核反应堆中的温度反馈,所述温度反馈为燃料和冷却剂的反应性温度反馈系数,得到总体的反应性方程(6);
<mrow>
<mfrac>
<mrow>
<msub>
<mi>dn</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&rho;</mi>
<mo>-</mo>
<mi>&beta;</mi>
</mrow>
<mi>&Lambda;</mi>
</mfrac>
<msub>
<mi>n</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mi>&Lambda;</mi>
</mfrac>
<msub>
<mi>c</mi>
<mi>r</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>dc</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&lambda;n</mi>
<mi>r</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&lambda;c</mi>
<mi>r</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>dT</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>f</mi>
<mi>f</mi>
</msub>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
</mfrac>
<msub>
<mi>n</mi>
<mi>r</mi>
</msub>
<mo>-</mo>
<mfrac>
<mi>&Omega;</mi>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
</mfrac>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<mfrac>
<mi>&Omega;</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>T</mi>
<mi>l</mi>
</msub>
<mo>+</mo>
<mfrac>
<mi>&Omega;</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>T</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>dT</mi>
<mi>l</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>c</mi>
</msub>
</mfrac>
<msub>
<mi>n</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<mfrac>
<mi>&Omega;</mi>
<msub>
<mi>&mu;</mi>
<mi>c</mi>
</msub>
</mfrac>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>M</mi>
<mo>+</mo>
<mi>&Omega;</mi>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&mu;</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>T</mi>
<mi>l</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>M</mi>
<mo>-</mo>
<mi>&Omega;</mi>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&mu;</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>T</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>d&rho;</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>G</mi>
<mi>r</mi>
</msub>
<msub>
<mi>z</mi>
<mi>r</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&rho;</mi>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>f</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<msub>
<mi>&alpha;</mi>
<mi>c</mi>
</msub>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>l</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>l</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中ρ为反应性;β为缓发中子份额;Λ为一代中子寿命,s;λ缓发中子先驱核衰变时间常数,s-1;nr相对中子密度;cr相对先驱密度;Tf为燃料平均温度,℃;ff为裂变释放能量在燃料中所占份额;p0为反应堆额定功率,MW;uf为燃料热容量;uc为冷却剂的热容量;Ω为燃料和冷却剂之间传热系数;Te是冷却剂入口温度;Tl冷却剂出口温度,℃;M为冷却剂质量流量与其热容量的乘积,MW/℃;ρr控制棒引入的反应性;Gr为控制棒微分价值;Zr控制棒速度;αf燃料反应性温度系数;Tf0为额定功率下燃料平均温度,℃;αc冷却剂反应性温度系数;Tl0额定功率下冷却剂出口温度,℃;
其中αf、αc、μc、Ω和M是nr0的函数:
M(nr0)=(28nr0+74)MW/℃;
在控制系统设计时,其中的常数取为β=0.006019,Λ=0.0001s,λ=0.15s-1,ff=0.92;uf=26.3MW·s/℃;Gr=0.0145δk/k/℃。
3.根据权利要求2所述的自抗扰控制方法,其特征在于,所述步骤S2具体为:
对(1)求导,将(1)、(2)代入可得
<mrow>
<msub>
<mover>
<mi>n</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>=</mo>
<mfrac>
<mover>
<mi>&rho;</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>&Lambda;</mi>
</mfrac>
<msub>
<mi>n</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&rho;</mi>
<mo>-</mo>
<mi>&beta;</mi>
</mrow>
<mi>&Lambda;</mi>
</mfrac>
<msub>
<mover>
<mi>n</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>-</mo>
<mi>&lambda;</mi>
<msub>
<mover>
<mi>n</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>+</mo>
<mfrac>
<mi>&rho;</mi>
<mi>&Lambda;</mi>
</mfrac>
<msub>
<mi>&lambda;n</mi>
<mi>r</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
设初始条件为nr(0)=1,cr(0)=1,ρ(0)=0,ρr(0)=0,Tf(0)=Tf0,Tl(0)=Tl0
设nr=nr(0)+δnr (8)
cr=cr(0)+δcr (9)
Tf=Tf(0)+δTf (10)
Tl=Tl(0)+δTl (11)
ρr=ρr(0)+δρr=δρr (12)
ρ=ρ(0)+δρ=δρ=δρr+δρf+δρl (13)
其中,
δρf=αfδTf (14)
δρl=αlδTl=(αc/2)δTl (15)
对(13)求导,并将(5)、(14)、(15)代入可得
<mrow>
<mover>
<mrow>
<mi>&delta;</mi>
<mi>&rho;</mi>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>G</mi>
<mi>r</mi>
</msub>
<msub>
<mi>z</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mi>f</mi>
</msub>
<mover>
<mrow>
<msub>
<mi>&delta;T</mi>
<mi>f</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mi>l</mi>
</msub>
<mover>
<mrow>
<msub>
<mi>&delta;T</mi>
<mi>l</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
由(8)可得
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<msub>
<mi>n</mi>
<mi>r</mi>
</msub>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mover>
<mrow>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</mtd>
<mtd>
<mrow>
<mover>
<msub>
<mi>n</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>&CenterDot;</mo>
<mo>&CenterDot;</mo>
</mrow>
</mover>
<mo>=</mo>
<mover>
<mrow>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mo>&CenterDot;</mo>
<mo>&CenterDot;</mo>
</mrow>
</mover>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
将(5)、(8)、(13)、(17)代入(7)可得
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mrow>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mo>&CenterDot;</mo>
<mo>&CenterDot;</mo>
</mrow>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>&Lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
<mover>
<mrow>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&beta;</mi>
<mi>&Lambda;</mi>
</mfrac>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mover>
<mrow>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mi>&Lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mi>r</mi>
</msub>
<msub>
<mi>z</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mi>f</mi>
</msub>
<mover>
<mrow>
<msub>
<mi>&delta;T</mi>
<mi>f</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mi>l</mi>
</msub>
<mover>
<mrow>
<msub>
<mi>&delta;T</mi>
<mi>l</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>r</mi>
</msub>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mi>&lambda;</mi>
<mi>&Lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>r</mi>
</msub>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
假设冷却剂入口温度不变,由稳态初始条件和(3)可得
<mrow>
<mfrac>
<mrow>
<msub>
<mi>dT</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>f</mi>
<mi>f</mi>
</msub>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
</mfrac>
<msub>
<mi>n</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mi>&Omega;</mi>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
</mfrac>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mi>&Omega;</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>T</mi>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mi>&Omega;</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>T</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
将(8)、(10)、(11)、(19)代入(3)可得
<mrow>
<mover>
<mrow>
<msub>
<mi>&delta;T</mi>
<mi>f</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>31</mn>
</msub>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>33</mn>
</msub>
<msub>
<mi>&delta;T</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>34</mn>
</msub>
<msub>
<mi>&delta;T</mi>
<mi>l</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
</mrow>
其中
<mrow>
<msub>
<mi>a</mi>
<mn>31</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>f</mi>
<mi>f</mi>
</msub>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
</mfrac>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>33</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mi>&Omega;</mi>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
</mfrac>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>34</mn>
</msub>
<mo>=</mo>
<mfrac>
<mi>&Omega;</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
</mrow>
</mfrac>
</mrow>
同理,由(4)、(8)、(10)、(11)可得
<mrow>
<mover>
<mrow>
<msub>
<mi>&delta;T</mi>
<mi>l</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>41</mn>
</msub>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>43</mn>
</msub>
<msub>
<mi>&delta;T</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>44</mn>
</msub>
<msub>
<mi>&delta;T</mi>
<mi>l</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
其中
<mrow>
<msub>
<mi>a</mi>
<mn>41</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>43</mn>
</msub>
<mo>=</mo>
<mfrac>
<mi>&Omega;</mi>
<msub>
<mi>&mu;</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>44</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>M</mi>
<mo>+</mo>
<mi>&Omega;</mi>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&mu;</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
</mrow>
将(20)、(21)代入(18)可得
<mrow>
<mover>
<mrow>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mo>&CenterDot;</mo>
<mo>&CenterDot;</mo>
</mrow>
</mover>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&beta;</mi>
<mi>&Lambda;</mi>
</mfrac>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mover>
<mrow>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<mfrac>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<mi>&Lambda;</mi>
</mfrac>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<mfrac>
<msub>
<mi>G</mi>
<mi>r</mi>
</msub>
<mi>&Lambda;</mi>
</mfrac>
<msub>
<mi>z</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<mi>g</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>22</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mrow>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mo>&CenterDot;</mo>
<mo>&CenterDot;</mo>
</mrow>
</mover>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mover>
<mrow>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>bz</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<mi>g</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>23</mn>
<mo>)</mo>
</mrow>
</mrow>
其中
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&beta;</mi>
<mi>&Lambda;</mi>
</mfrac>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>f</mi>
</msub>
<mo>*</mo>
<msub>
<mi>f</mi>
<mi>f</mi>
</msub>
<mo>*</mo>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
<mo>/</mo>
<msub>
<mi>&mu;</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mi>c</mi>
</msub>
<mo>*</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
<mo>/</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>*</mo>
<msub>
<mi>&mu;</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>&Lambda;</mi>
</mfrac>
</mrow>
<mrow>
<mi>b</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>G</mi>
<mi>r</mi>
</msub>
<mi>&Lambda;</mi>
</mfrac>
</mrow>
k1=αfa31+αla41 k2=αfa33+αla43 k3=αfa34+αla44
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>&Lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
<mover>
<mrow>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mi>&Lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&delta;T</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&delta;T</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mi>&Lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mi>r</mi>
</msub>
<msub>
<mi>z</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&delta;T</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&delta;T</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mi>&lambda;</mi>
<mi>&Lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&delta;&rho;</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>&delta;n</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
这里g可以看作模型不确定项。
4.根据权利要求3所述的自抗扰控制方法,其特征在于,所述步骤S3具体为:
对于二阶对象
<mrow>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mi>y</mi>
<mo>+</mo>
<mi>b</mi>
<mi>u</mi>
<mo>+</mo>
<mi>g</mi>
<mo>+</mo>
<mi>w</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>24</mn>
<mo>)</mo>
</mrow>
</mrow>
y和u分别是输出和输入,g代表模型不确定项,w是外部扰动;所述(24)中参数a1,a2,b为已知;令f=g+w为总扰动,所述(24)的状态方程为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>A</mi>
<mi>x</mi>
<mo>+</mo>
<mi>B</mi>
<mi>u</mi>
<mo>+</mo>
<mi>E</mi>
<mi>h</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mi>C</mi>
<mi>x</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>b</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>E</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>C</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,x3=f是扩张状态,h为f的导数,则LESO为:
<mrow>
<mover>
<mi>z</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>A</mi>
<mi>z</mi>
<mo>+</mo>
<mi>B</mi>
<mi>u</mi>
<mo>+</mo>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>-</mo>
<mi>L</mi>
<mi>C</mi>
<mo>)</mo>
</mrow>
<mi>z</mi>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>B</mi>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>u</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>y</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>25</mn>
<mo>)</mo>
</mrow>
</mrow>
3
<mrow>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mi>C</mi>
<mi>z</mi>
</mrow>
其中L是观测器增益,能够用极点配置方法获得L=[l1 l2 l3]T,
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>A</mi>
<mo>-</mo>
<mi>L</mi>
<mi>C</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>L</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mn>3</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mn>3</mn>
<msub>
<mi>&omega;</mi>
<mi>o</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>3</mn>
<msup>
<msub>
<mi>&omega;</mi>
<mi>o</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>3</mn>
<msub>
<mi>&omega;</mi>
<mi>o</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<msub>
<mi>&omega;</mi>
<mi>o</mi>
</msub>
<mn>3</mn>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
LESO:
A*=A-LC B*=[B L] C*=I3 D*=(0)3×2,
控制器设计为:
<mrow>
<mi>u</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>u</mi>
<mn>0</mn>
</msub>
</mrow>
<mi>b</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>26</mn>
<mo>)</mo>
</mrow>
</mrow>
PD控制器:
u0=kp(r-z1)-kdz2 (27)
r为设定值,闭环传递函数为:
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>c</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>k</mi>
<mi>p</mi>
</msub>
<mrow>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>k</mi>
<mi>d</mi>
</msub>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>k</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>28</mn>
<mo>)</mo>
</mrow>
</mrow>
则增益取为:kp=ωc 2 kd=2ωc。
5.根据权利要求4所述的自抗扰控制方法,其特征在于,所述步骤S4具体为:
RLESO:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>z</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>A</mi>
<mi>z</mi>
<mo>+</mo>
<mi>B</mi>
<mi>u</mi>
<mo>+</mo>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>-</mo>
<mover>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>-</mo>
<mi>L</mi>
<mi>C</mi>
<mo>)</mo>
</mrow>
<mi>z</mi>
<mo>+</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mi>B</mi>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mi>u</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mover>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>^</mo>
</mover>
<mo>,</mo>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>29</mn>
<mo>)</mo>
</mrow>
</mrow>
这里
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>B</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mi>b</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>C</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>L</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中L是观测器增益,能够用极点配置方法获得,l1=2ωo,l2=ωo 2,
令
V=z1-2ωoy (30)
W=z2-ωo 2y (31)
则由(29)得
<mrow>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mover>
<mi>V</mi>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
<mtr>
<mtd>
<mover>
<mi>W</mi>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&omega;</mi>
<mi>o</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<msub>
<mi>&omega;</mi>
<mi>o</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mi>V</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>W</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mi>b</mi>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>3</mn>
<msup>
<msub>
<mi>&omega;</mi>
<mi>o</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
<msup>
<msub>
<mi>&omega;</mi>
<mi>o</mi>
</msub>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mi>u</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>y</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>32</mn>
<mo>)</mo>
</mrow>
</mrow>
因此
<mrow>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mover>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>^</mo>
</mover>
<mo>=</mo>
<mi>V</mi>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&omega;</mi>
<mi>o</mi>
</msub>
<mi>y</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>33</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mi>W</mi>
<mo>+</mo>
<msup>
<msub>
<mi>&omega;</mi>
<mi>o</mi>
</msub>
<mn>2</mn>
</msup>
<mi>y</mi>
<mo>.</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>34</mn>
<mo>)</mo>
</mrow>
</mrow>
6.根据权利要求5所述的自抗扰控制方法,其特征在于,所述步骤S5具体为:通过自抗扰控制器的设计原理,整定参数ωc、ωo,确定该反应堆功率控制系统;
当nr0=1时,即额定工况下,可以求得此时只有ωc和ωo需要整定;
在实际整定过程中遵循如下规律:
1)ωo越大,LESO的观测能力越强,观测器对噪声的敏感性增加,ωo从较小的值逐渐增大,直至观测精度满足要求为止;
2)ωc越大,控制作用越强,系统的响应速度越快,但超调和振荡会越严重,同时稳定性下降;
ωo和ωc在保证控制器稳定的前提下通过多次试验和综合比较误差指标加以确定,观测器及控制器的误差上界与其带宽成反比,带宽越宽,误差越小,稳定性越低。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610827618.4A CN106340331B (zh) | 2016-09-18 | 2016-09-18 | 一种用于核反应堆功率的自抗扰控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610827618.4A CN106340331B (zh) | 2016-09-18 | 2016-09-18 | 一种用于核反应堆功率的自抗扰控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106340331A CN106340331A (zh) | 2017-01-18 |
CN106340331B true CN106340331B (zh) | 2017-12-05 |
Family
ID=57840026
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610827618.4A Active CN106340331B (zh) | 2016-09-18 | 2016-09-18 | 一种用于核反应堆功率的自抗扰控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106340331B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108983602B (zh) * | 2018-06-26 | 2021-04-27 | 华北电力大学 | 一种用于快堆功率和冷却剂出口温度的自抗扰控制方法 |
CN109116721B (zh) * | 2018-08-23 | 2021-10-19 | 广东工业大学 | 一种将时变系统转变为定常系统的控制方法 |
CN110289114B (zh) * | 2019-06-10 | 2021-03-09 | 清华大学 | 高温气冷堆功率控制方法、控制系统及高温气冷堆核电站 |
CN111665717B (zh) * | 2020-05-28 | 2021-05-14 | 华北电力大学 | 用于大型压水堆轴向功率分布的线性自抗扰控制建模方法 |
CN112764346B (zh) * | 2020-12-24 | 2022-04-01 | 暨南大学 | 分散式自抗扰多变量控制方法 |
CN113031435B (zh) * | 2021-02-03 | 2022-07-12 | 北京航空航天大学 | 数字式双频扩张状态观测器及扰动观测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102354104A (zh) * | 2005-09-19 | 2012-02-15 | 克利夫兰州立大学 | 控制器、观测器及其应用 |
CN105305913A (zh) * | 2015-10-30 | 2016-02-03 | 西安交通大学苏州研究院 | 一种新型的用于滚珠丝杠进给系统的抗扰跟随控制器 |
WO2016118701A1 (en) * | 2015-01-21 | 2016-07-28 | Linestream Technologies | Cascaded active disturbance rejection controllers |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8041436B2 (en) * | 2002-04-18 | 2011-10-18 | Cleveland State University | Scaling and parameterizing a controller |
-
2016
- 2016-09-18 CN CN201610827618.4A patent/CN106340331B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102354104A (zh) * | 2005-09-19 | 2012-02-15 | 克利夫兰州立大学 | 控制器、观测器及其应用 |
WO2016118701A1 (en) * | 2015-01-21 | 2016-07-28 | Linestream Technologies | Cascaded active disturbance rejection controllers |
CN105305913A (zh) * | 2015-10-30 | 2016-02-03 | 西安交通大学苏州研究院 | 一种新型的用于滚珠丝杠进给系统的抗扰跟随控制器 |
Non-Patent Citations (1)
Title |
---|
反应堆功率线性自抗扰控制方法研究;刘玉燕等;《控制工程》;20150930;第2-4节,图1-2 * |
Also Published As
Publication number | Publication date |
---|---|
CN106340331A (zh) | 2017-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106340331B (zh) | 一种用于核反应堆功率的自抗扰控制方法 | |
CN108983602B (zh) | 一种用于快堆功率和冷却剂出口温度的自抗扰控制方法 | |
CN108962410B (zh) | 一种用于铅冷快堆功率的自抗扰控制方法 | |
Kothare et al. | Level control in the steam generator of a nuclear power plant | |
CN106919053A (zh) | 一种基于变结构预测控制算法的火电机组协调控制系统 | |
CN107065518B (zh) | 一种超临界机组协调预测函数控制算法 | |
Bose et al. | An interval approach to nonlinear controller design for load-following operation of a small modular pressurized water reactor | |
US11682497B2 (en) | Method for regulating operating parameters of a nuclear reactor and corresponding nuclear reactor | |
Ansarifar et al. | Gain scheduled dynamic sliding mode control for nuclear steam generators | |
Menon et al. | Gain-scheduled nonlinear control of U-tube steam generator water level | |
Surjagade et al. | Robust optimal integral sliding mode controller for total power control of large PHWRs | |
Dong | Model-free power-level control of MHTGRs against input saturation and dead-zone | |
Hui et al. | Load following control of a pressurized water reactor via finite-time super-twisting sliding mode and extended state observer techniques | |
CN107870564A (zh) | 具有快速、安全、解耦性能的燃料电池抗干扰控制方法 | |
CN106855691B (zh) | 用于超临界火电机组机炉系统负荷跟踪的双层控制系统 | |
Liu et al. | Linear active disturbance rejection control for pressurized water reactor power based on partial feedback linearization | |
Qiu et al. | Research on fuzzy weighted gain scheduling water level control system of U-tube steam generator | |
Hui et al. | Adaptive second-order nonsingular terminal sliding mode power-level control for nuclear power plants | |
Hui | Discrete-time integral terminal sliding mode load following controller coupled with disturbance observer for a modular high-temperature gas-cooled reactor | |
Surjagade et al. | An arbitrary-order continuous sliding mode control technique for nonlinear PWR-type nuclear power plants | |
Reddy et al. | Backstepping based model reference adaptive control for nuclear reactor with matched and unmatched uncertainties | |
Reddy et al. | A robust nonlinear adaptive control for control of nuclear reactor | |
Zhu et al. | Internal model control using LMS filter and its application to superheated steam temperature of power plant | |
Nakamura et al. | Optimal regulation for thermal power plants | |
Andraws et al. | Performance of receding horizon predictive controller for research reactors |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |