CN104749959B - 一种基于广义滑模估计器的机组变桨容错控制方法 - Google Patents
一种基于广义滑模估计器的机组变桨容错控制方法 Download PDFInfo
- Publication number
- CN104749959B CN104749959B CN201510205113.XA CN201510205113A CN104749959B CN 104749959 B CN104749959 B CN 104749959B CN 201510205113 A CN201510205113 A CN 201510205113A CN 104749959 B CN104749959 B CN 104749959B
- Authority
- CN
- China
- Prior art keywords
- overbar
- matrix
- fault
- delta
- ranks
- 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.)
- Expired - Fee Related
Links
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于广义滑模估计器的机组变桨容错控制方法,包括如下步骤:1)建立同时存在执行器故障和传感器故障下的变桨系统模型;2)引入增广状态向量,将原来故障下变桨系统模型变换成广义系统模型;3)基于广义系统模型设计广义滑模估计器,实现对桨叶位置状态估计和系统故障重构;4)基于状态估计设计容错控制方法,实现故障系统的状态稳定;5)根据确保估计误差系统和原状态系统稳定的线性矩阵不等式条件,计算参数矩阵,把各参数矩阵配置给所设计的估计器和控制器。本发明建立了估计误差系统和原状态系统稳定的条件,确保估计误差收敛性和状态的有界性,最终实现对执行器故障和传感器故障的重构,以及对故障变桨系统的容错控制。
Description
技术领域
本发明涉及一种基于估计器的容错控制方法,尤其涉及一种基于广义滑模估计器的机组变桨容错控制方法。
背景技术
对于风力发电这类实际工程系统,长期运行时通常不可避免地出现执行器故障或是传感器故障。对于一个发生故障的动态系统,如何维持系统稳定成为控制系统设计的一个重要课题。
目前,大多数容错控制方法都建立在传统的故障检测与诊断(FDD)基础上。但是,现存故障检测与诊断(FDD)技术并不能直接反映故障对系统状态的影响程度。于是在FDD基础上,提出了故障重构的思想:根据已被故障影响的过程变量的测量值来估计出故障的大小和位置,几乎可以将故障信号完整地重现出来,不但能够直观地反映出故障的发生、发展过程,还可通过对信号的分析方便地获得关于故障的类型、性质、严重程度以及可能对系统造成的影响、危害等一系列信息。
目前,针对风力变桨系统,还没有完整地在模型里考虑其非线性、不确定性以及多重故障(包括执行器故障和传感器故障)同时存在的情况。针对变桨系统这样复杂的故障情况,对执行器故障或是传感器故障的重构和还原造成了极大的难度。目前,对于故障重构方法的研究成果主要体现在基于观测器或估计器的故障重构技术上。目前主要有三大类:一是基于广义观测器的估计方法;二是基于滑模观测器的故障重构技术;三是基于鲁棒自适应观测器的方法。这三类方法在不同的应用场合体现着自身的优缺点。但在现有的研究成果中,还没有将故障重构的方法有效地应用到机组变桨的故障系统中去,进而同时实现对执行器故障和传感器故障的重构并有效地获取故障信息,采取更有效的容错控制措施。
发明内容
针对现有技术中风电变桨系统存在的上述不足,本发明充分考虑了系统的非线性、不确定性以及同时存在执行器故障和传感器故障的情况,提供了一种基于广义滑模估计器的机组变桨容错控制方法,不仅能实现执行器故障和传感器故障的重构,还能实现故障系统的状态稳定。
为了解决上述技术问题,本发明采用了如下技术方案:
一种基于广义滑模估计器的机组变桨容错控制方法,其特征在于,该方法包括如下步骤:
1)建立同时存在执行器故障和传感器故障下的变桨系统模型;
首先,考虑存在执行器故障下的叶轮动态方程:
再次,考虑存在传感器故障下的桨距角测量方程为:
y(t)=x(t)+fs(t) (2)
其中, 分别代表实际的桨距角、相应的角速度以及相应的角加速度,是桨距角的实际测量值,是推进器、叶片升力、叶片弯曲、振动及不平衡载荷引起的总力矩,是控制输入,是执行器故障函数,是传感器故障函数,JLB为空气加速引起的等效惯量力矩,JBI是转子叶片沿其纵轴的全部惯性力矩,kDB是阻尼系数,kRL是摩擦系数;
进一步,方程(1)可以写成:
其中,
而测量方程(2)可以写成:
Y(t)=CX(t)+Dsfs(t) (4)
其中为已知的2乘1阶的常值矩阵;
2)引入增广状态向量X,将原来故障下变桨系统模型变换成一种新型的广义系统模型;
定义如下的增广矩阵:
其中,Ip代表维数为p单位矩阵,0k×l代表k行l列的零矩阵,参数δ>0要使得下面条件(6)成立;
其中rank表示矩阵的秩;
最终,可以得到一个增广的广义系统模型(7);
其中,是式(5)中所定义的维数为5的广义系统向量,是维数为2的广义系统测量输出向量,是式(5)中所定义的维数为5的广义非线性向量,是式(5)中所定义的5乘5阶的矩阵,是式(5)中所定义的5乘5阶的参数矩阵,是式(5)中所定义的5乘1阶的参数矩阵,是式(5)中所定义的5乘3阶的矩阵,是式(5)中所定义的2乘5阶的参数矩阵,是式(5)中所定义的维数为3的广义故障向量;
3)基于广义系统模型(7)设计出以下结构的广义滑模估计器(8),实现对桨叶位置状态估计和系统故障重构;
其中,是维数为5的中间变量向量,是维数为5的增广的广义系统状态估计向量,是式(5)中所定义的5乘2阶的矩阵, 分别是所设计的5乘2阶的微分增益矩阵,5乘2阶的比例增益矩阵以及5乘3阶的滑模增益矩阵,是维数为3的连续的输入信号;
3.1)设计微分增益则可证明矩阵是非奇异的;
3.2)设计比例增益
首先,利用如下李雅普诺夫方程(9)求解出矩阵
其中,μ>0满足
进一步,比例增益可通过式子(10)进行计算;
3.3)设计滑模增益
3.4)令和定义如下滑模面:
其中,为正定矩阵,且满足以下条件:
其中,是由所确定的3乘2阶的设计矩阵;
于是,最终设计的滑模面变为:
设计的连续输入信号具有如下的滑模形式:
其中,δ>0满足条件(6),γ>0为设计参数,ο>0是一个小常数,取ο=10-3,α0,α1,β0和d0满足以下条件(15);
4)基于状态估计设计容错控制方法,实现故障系统的状态稳定;
4.1)基于状态估计值,定义一个积分型滑模面(16):
其中,和分别是所设计1乘2阶的参数矩阵;矩阵G要使得GB=1,矩阵K要使得矩阵(A+BK)是霍尔维茨的;
4.2)所设计的控制器结构(17)-(18)如下所示:
且
其中,
5)给出估计误差系统和原状态系统稳定的条件,计算参数矩阵,并把各参数矩阵配置给所设计的估计器和控制器;
稳定性条件:在滑模控制us(t)和控制输入u(t)下,如果存在2乘2阶的正定矩阵2乘2阶的正定矩阵和3乘2阶的矩阵使得下面LMI优化问题(20)成立:
其中,υ是固定常数,
则可以确保估计误差系统和原状态系统的稳定;
通过LMI工具箱的“mincx”求解器,可以方便地求解以上线性矩阵不等式优化问题,计算出所设计的参数矩阵
最后,在风电变桨系统中,将LMI优化问题中求解出的参数矩阵同以上步骤中所设计的微分增益矩阵比例增益矩阵滑模增益矩阵控制器参数矩阵G和K相应配置给所设计的估计器和控制器,实现对桨叶位置的估计和对执行器故障和传感器故障的重构,实现变桨系统稳定。
与现有技术相比,本发明具有如下优点:
1、首次将广义滑模估计器的故障重构方法应用到机组变桨的故障模型中,有效地实现故障的重构和系统的稳定。
2、通过LMI工具箱来求解LMI优化问题,可以容易地得到所设计的估计器增益矩阵,参数设计简单,且易于实现。
附图说明
图1为一种基于广义滑模估计器的机组变桨容错控制方法的流程图;
图2为实际桨距角与估计桨距角的曲线图;
图3为实际的执行器故障曲线和估计的执行器故障曲线对比图;
图4为实际的传感器故障曲线和估计的传感器故障曲线对比图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细说明。
一种基于广义滑模估计器的机组变桨容错控制方法,其流程如图1所示,该方法包括如下步骤:
1)建立同时存在执行器故障和传感器故障下的变桨系统模型;
首先,考虑存在执行器故障下的叶轮动态方程:
再次,考虑存在传感器故障下的桨距角测量方程为:
y(t)=x(t)+fs(t) (2)
其中,M(·)=JLB+JBI>0; 分别代表实际的桨距角、相应的角速度以及相应的角加速度,是桨距角的实际测量值,是推进器、叶片升力、叶片弯曲、振动及不平衡载荷引起的总力矩,是控制输入,是执行器故障函数,是传感器故障函数,JLB为空气加速引起的等效惯量力矩,JBI是转子叶片沿其纵轴的全部惯性力矩,kDB是阻尼系数,kRL是摩擦系数;
进一步,方程(1)可以写成:
其中,
而测量方程(2)可以写成:
Y(t)=CX(t)+Dsfs(t) (4)
其中为已知的2乘1阶的常值矩阵。
2)引入增广状态向量,将原来故障下变桨系统模型变换成一种新型的广义系统模型;
定义如下的增广矩阵:
其中,Ip代表维数为p单位矩阵,0k×l代表k行l列的零矩阵,参数δ>0要使得下面条件(6)成立;
其中rank表示矩阵的秩;
最终,可以得到一个增广的广义系统模型(7);
其中,是式(5)中所定义的维数为5的广义系统向量,是维数为2的广义系统测量输出向量,是式(5)中所定义的维数为5的广义非线性向量,是式(5)中所定义的5乘5阶的矩阵,是式(5)中所定义的5乘5阶的参数矩阵,是式(5)中所定义的5乘1阶的参数矩阵,是式(5)中所定义的5乘3阶的矩阵,是式(5)中所定义的2乘5阶的参数矩阵,是式(5)中所定义的维数为3的广义故障向量。
3)基于广义系统模型(7)设计出以下结构的广义滑模估计器(8),实现对桨叶位置状态估计和系统故障重构;
其中,是维数为5的中间变量向量,是维数为5的增广的广义系统状态估计向量,是式(5)中所定义的5乘2阶的矩阵, 分别是所设计的5乘2阶的微分增益矩阵,5乘2阶的比例增益矩阵以及5乘3阶的滑模增益矩阵,是维数为3的连续的输入信号。
3.1)设计微分增益则可证明矩阵是非奇异的;
3.2)设计比例增益
首先,利用如下李雅普诺夫方程(9)求解出矩阵
其中,μ>0满足
进一步,比例增益可通过式子(10)进行计算;
3.3)设计滑模增益
3.4)令和定义如下滑模面:
其中,为正定矩阵,且满足以下条件:
其中,是由所确定的3乘2阶的设计矩阵;
于是,最终设计的滑模面变为:
设计的连续输入信号具有如下的滑模形式:
其中,δ>0满足条件(6),γ>0为设计参数,ο>0是一个小常数,取ο=10-3,α0,α1,β0和d0满足以下条件(15);
4)基于状态估计设计容错控制方法,实现故障系统的状态稳定;
4.1)基于状态估计值,定义一个积分型滑模面(16):
其中,和分别是所设计1乘2阶的参数矩阵;矩阵G要使得GB=1,矩阵K要使得矩阵(A+BK)是霍尔维茨的;
4.2)所设计的控制器结构(17)-(18)如下所示:
且
其中,
5)给出估计误差系统和原状态系统稳定的条件,计算参数矩阵,并把各参数矩阵配置给所设计的估计器和控制器;
稳定性条件:在滑模控制us(t)和控制输入u(t)下,如果存在2乘2阶的正定矩阵2乘2阶的正定矩阵和3乘2阶的矩阵使得下面LMI优化问题(20)成立:
其中,υ是固定常数,则可以确保估计误差系统和原状态系统的稳定。
通过LMI工具箱的“mincx”求解器,可以方便地求解以上线性矩阵不等式优化问题,计算出所设计的参数矩阵
最后,在风电变桨系统中,将LMI优化问题中求解出的参数矩阵同以上步骤中所设计的微分增益矩阵比例增益矩阵滑模增益矩阵控制器参数矩阵G和K相应配置给所设计的估计器和控制器,实现对桨叶位置的估计和对执行器故障和传感器故障的重构,实现变桨系统稳定。
在机组变桨中采用上述设计的广义滑模估计器结构(8)和容错控制率(17)-(18),按照上述的矩阵参数设计步骤,可以实现对桨距角及桨距角速度的估计,执行器和传感器故障的重构以及桨叶位置的稳定。图2表明通过本发明的广义滑模估计器和容错控制方法,桨距角位置不仅能得到估计,而且能实现状态稳定;图3表明执行器故障的估计情况;图4表明传感器故障的估计情况,均表明了本发明的广义滑模估计器具有较好的故障估计效果。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (1)
1.一种基于广义滑模估计器的机组变桨容错控制方法,其特征在于,该方法包括如下步骤:
1)建立同时存在执行器故障和传感器故障下的变桨系统模型;
首先,考虑存在执行器故障下的叶轮动态方程:
再次,考虑存在传感器故障下的桨距角测量方程为:
y(t)=x(t)+fs(t) (2)
其中,M(·)=JLB+JBI>0; 分别代表实际的桨距角、相应的角速度以及相应的角加速度,是桨距角的实际测量值,是推进器、叶片升力、叶片弯曲、振动及不平衡载荷引起的总力矩,是控制输入,是执行器故障函数,是传感器故障函数,JLB为空气加速引起的等效惯量力矩,JBI是转子叶片沿其纵轴的全部惯性力矩,kDB是阻尼系数,kRL是摩擦系数;
进一步,方程(1)可以写成:
其中,
而测量方程(2)可以写成:
Y(t)=CX(t)+Dsfs(t) (4)
其中为已知的2乘1阶的常值矩阵;
2)引入增广状态向量将原来故障下变桨系统模型变换成一种新型的广义系统模型;
定义如下的增广矩阵:
其中,Ip代表维数为p单位矩阵,0k×l代表k行l列的零矩阵,参数δ>0要使得下面条件(6)成立;
其中rank表示矩阵的秩;
最终,可以得到一个增广的广义系统模型(7);
其中,是式(5)中所定义的维数为5的广义系统向量,是维数为2的广义系统测量输出向量,是式(5)中所定义的维数为5的广义非线性向量,是式(5)中所定义的5乘5阶的矩阵,是式(5)中所定义的5乘5阶的参数矩阵,是式(5)中所定义的5乘1阶的参数矩阵,是式(5)中所定义的5乘3阶的矩阵,是式(5)中所定义的2乘5阶的参数矩阵,是式(5)中所定义的维数为3的广义故障向量;
3)基于广义系统模型(7)设计出以下结构的广义滑模估计器(8),实现对桨叶位置状态估计和系统故障重构;
其中,是维数为5的中间变量向量,是维数为5的增广的广义系统状态估计向量,是式(5)中所定义的5乘2阶的矩阵, 分别是所设计的5乘2阶的微分增益矩阵,5乘2阶的比例增益矩阵以及5乘3阶的滑模增益矩阵,是维数为3的连续的输入信号;
3.1)设计微分增益则可证明矩阵是非奇异的;
3.2)设计比例增益
首先,利用如下李雅普诺夫方程(9)求解出矩阵
其中,μ>0满足
进一步,比例增益可通过式子(10)进行计算;
3.3)设计滑模增益
3.4)令和定义如下滑模面:
其中,为正定矩阵,且满足以下条件:
其中,是由所确定的3乘2阶的设计矩阵;
于是,最终设计的滑模面变为:
设计的连续输入信号具有如下的滑模形式:
其中,δ>0满足条件(6),γ>0为设计参数,ο>0是一个小常数,取ο=10-3,α0,α1,β0和d0满足以下条件(15);
4)基于状态估计设计容错控制方法,实现故障系统的状态稳定;
4.1)基于状态估计值,定义一个积分型滑模面(16):
其中,和分别是所设计1乘2阶的参数矩阵;矩阵G要使得GB=1,矩阵K要使得矩阵(A+BK)是霍尔维茨的;
4.2)所设计的控制器结构(17)-(18)如下所示:
且
其中,
5)给出估计误差系统和原状态系统稳定的条件,计算参数矩阵,并把各参数矩阵配置给所设计的估计器和控制器;
稳定性条件:在滑模控制us(t)和控制输入u(t)下,如果存在2乘2阶的正定矩阵2乘2阶的正定矩阵和3乘2阶的矩阵使得下面LMI优化问题(20)成立:
st.
其中,υ是固定常数,
则可以确保估计误差系统和原状态系统的稳定;
通过LMI工具箱的“mincx”求解器,可以方便地求解以上线性矩阵不等式优化问题,计算出所设计的参数矩阵
最后,在风电变桨系统中,将LMI优化问题中求解出的参数矩阵同以上步骤中所设计的微分增益矩阵比例增益矩阵滑模增益矩阵控制器参数矩阵G和K相应配置给所设计的估计器和控制器,实现对桨叶位置的估计和对执行器故障和传感器故障的重构,实现变桨系统稳定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510205113.XA CN104749959B (zh) | 2015-04-27 | 2015-04-27 | 一种基于广义滑模估计器的机组变桨容错控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510205113.XA CN104749959B (zh) | 2015-04-27 | 2015-04-27 | 一种基于广义滑模估计器的机组变桨容错控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104749959A CN104749959A (zh) | 2015-07-01 |
CN104749959B true CN104749959B (zh) | 2017-05-10 |
Family
ID=53589828
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510205113.XA Expired - Fee Related CN104749959B (zh) | 2015-04-27 | 2015-04-27 | 一种基于广义滑模估计器的机组变桨容错控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104749959B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113969870B (zh) * | 2020-07-23 | 2023-07-25 | 北京金风科创风电设备有限公司 | 用于风力发电机组估计器的监测方法及其装置 |
CN112431717B (zh) * | 2020-11-23 | 2021-11-26 | 山东科技大学 | 一种小型风力机叶片统一变桨系统及变桨方法 |
CN112711190B (zh) * | 2020-12-25 | 2022-03-08 | 四川大学 | 一种自适应容错控制器、控制设备及控制系统 |
CN113998001B (zh) * | 2021-11-05 | 2024-01-16 | 合肥工业大学 | 一种无人驾驶车辆线控转向的容错控制器及其设计方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101900080B (zh) * | 2010-07-21 | 2011-11-23 | 上海电气集团股份有限公司 | 采用变结构pid变桨控制的风机控制系统 |
CN102102630B (zh) * | 2010-10-26 | 2012-09-19 | 重庆大学 | 变速变桨风力发电机组独立变桨控制方法 |
CN103278328B (zh) * | 2013-05-16 | 2015-06-10 | 华北电力大学 | 基于模糊Petri网的风电机组液压变桨系统故障诊断方法 |
-
2015
- 2015-04-27 CN CN201510205113.XA patent/CN104749959B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN104749959A (zh) | 2015-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104749959B (zh) | 一种基于广义滑模估计器的机组变桨容错控制方法 | |
CN104850889B (zh) | 一种飞机旋转作动器驱动装置自适应故障检测、隔离与健康评估的方法 | |
CN104200122B (zh) | 复杂焊接结构随机振动疲劳寿命预测方法 | |
Odgaard et al. | Fault tolerant control of wind turbines using unknown input observers | |
CN107329083A (zh) | 针对高速列车牵引电机非线性系统传感器的故障诊断方法 | |
CN106026822A (zh) | 伺服电机驱动系统的惯量在线辨识方法及系统 | |
Jesussek et al. | Fault detection and isolation for a full-scale railway vehicle suspension with multiple Kalman filters | |
CN105171758A (zh) | 一种机器人的自适应有限时间收敛滑模控制方法 | |
CN104462785A (zh) | 一种两阶段式建筑框架结构损伤检测方法 | |
CN107356523B (zh) | 基于离散切线刚度估计的实时混合模拟试验反馈力修正方法 | |
Shi et al. | A robust adaptive approach to wind turbine pitch actuator component fault estimation | |
Fan et al. | Dynamic state estimation and parameter calibration of a DFIG using the ensemble Kalman filter | |
CN109753744A (zh) | 一种基于简化模型的框架结构模型修正方法 | |
CN106326530A (zh) | 一种基于右矩阵分式模型的时变结构模态参数辨识方法 | |
KR20130094071A (ko) | 부유식 해상풍력발전 하부구조의 동적피로하중산출 시뮬레이션 방법 | |
CN111173688A (zh) | 基于自适应观测器的风力发电机故障诊断与隔离方法 | |
CN107194063A (zh) | 一种高效结构失效概率函数求解的扩展线抽样方法 | |
CN102288159A (zh) | 用于确定塔的偏转的方法和装置 | |
Bila | Power system dynamic state estimation and load modeling | |
Bottasso et al. | Spatial estimation of wind states from the aeroelastic response of a wind turbine | |
Yu et al. | Estimation of unknown large power system dynamics | |
CN113890442B (zh) | 一种粒子群滑模观测器风力发电机转子电流故障检测方法 | |
CN112448410A (zh) | 应用于含光伏电站电力系统低频振荡的模态综合分析方法 | |
CN112925206B (zh) | 一种非线性多倒立摆互联系统分布式鲁棒故障诊断方法 | |
CN105337547A (zh) | 异步电机无速度传感器矢量控制系统的转速辨识方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170510 Termination date: 20180427 |