首页 / 专利库 / 电脑图像 / 霍夫变换 / 广义霍夫变换 / 一种基于广义滑模估计器的机组变桨容错控制方法

一种基于广义滑模估计器的机组变桨容错控制方法

阅读:691发布:2020-05-26

专利汇可以提供一种基于广义滑模估计器的机组变桨容错控制方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种基于广义滑模估计器的机组变桨容错控制方法,包括如下步骤:1)建立同时存在执行器故障和 传感器 故障下的变桨系统模型;2)引入增广状态向量,将原来故障下变桨系统模型变换成广义系统模型;3)基于广义系统模型设计广义滑模估计器,实现对桨叶 位置 状态估计和系统故障重构;4)基于状态估计设计容错控制方法,实现故障系统的状态稳定;5)根据确保估计误差系统和原状态系统稳定的线性矩阵不等式条件,计算参数矩阵,把各参数矩阵配置给所设计的估计器和 控制器 。本发明建立了估计误差系统和原状态系统稳定的条件,确保估计误差收敛性和状态的有界性,最终实现对执行器故障和传感器故障的重构,以及对故障变桨系统的容错控制。,下面是一种基于广义滑模估计器的机组变桨容错控制方法专利的具体信息内容。

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相应配置给所设计的估计器和控制器,实现对桨叶位置的估计和对执行器故障和传感器故障的重构,实现变桨系统稳定。

说明书全文

一种基于广义滑模估计器的机组变桨容错控制方法

技术领域

[0001] 本发明涉及一种基于估计器的容错控制方法,尤其涉及一种基于广义滑模估计器的机组变桨容错控制方法。

背景技术

[0002] 对于发电这类实际工程系统,长期运行时通常不可避免地出现执行器故障或是传感器故障。对于一个发生故障的动态系统,如何维持系统稳定成为控制系统设计的一个重要课题。
[0003] 目前,大多数容错控制方法都建立在传统的故障检测与诊断(FDD)基础上。但是,现存故障检测与诊断(FDD)技术并不能直接反映故障对系统状态的影响程度。于是在FDD基础上,提出了故障重构的思想:根据已被故障影响的过程变量的测量值来估计出故障的大小和位置,几乎可以将故障信号完整地重现出来,不但能够直观地反映出故障的发生、发展过程,还可通过对信号的分析方便地获得关于故障的类型、性质、严重程度以及可能对系统造成的影响、危害等一系列信息。
[0004] 目前,针对风力变桨系统,还没有完整地在模型里考虑其非线性、不确定性以及多重故障(包括执行器故障和传感器故障)同时存在的情况。针对变桨系统这样复杂的故障情况,对执行器故障或是传感器故障的重构和还原造成了极大的难度。目前,对于故障重构方法的研究成果主要体现在基于观测器或估计器的故障重构技术上。目前主要有三大类:一是基于广义观测器的估计方法;二是基于滑模观测器的故障重构技术;三是基于鲁棒自适应观测器的方法。这三类方法在不同的应用场合体现着自身的优缺点。但在现有的研究成果中,还没有将故障重构的方法有效地应用到机组变桨的故障系统中去,进而同时实现对执行器故障和传感器故障的重构并有效地获取故障信息,采取更有效的容错控制措施。

发明内容

[0005] 针对现有技术中风电变桨系统存在的上述不足,本发明充分考虑了系统的非线性、不确定性以及同时存在执行器故障和传感器故障的情况,提供了一种基于广义滑模估计器的机组变桨容错控制方法,不仅能实现执行器故障和传感器故障的重构,还能实现故障系统的状态稳定。
[0006] 为了解决上述技术问题,本发明采用了如下技术方案:
[0007] 一种基于广义滑模估计器的机组变桨容错控制方法,其特征在于,该方法包括如下步骤:
[0008] 1)建立同时存在执行器故障和传感器故障下的变桨系统模型;
[0009] 首先,考虑存在执行器故障下的叶轮动态方程:
[0010]
[0011] 再次,考虑存在传感器故障下的桨距测量方程为:
[0012] y(t)=x(t)+fs(t)  (2)
[0013] 其中,分别代表实际的桨距角、相应的角速度以及相应的角加速度,
是桨距角的实际测量值, 是推进器叶片升力、叶片弯曲、振动及不平衡
载荷引起的总力矩, 是控制输入, 是执行器故障函数, 是传感器故
障函数,JLB为空气加速引起的等效惯量力矩,JBI是转子叶片沿其纵轴的全部惯性力矩,kDB是阻尼系数,kRL是摩擦系数
[0014] 进一步,方程(1)可以写成:
[0015]
[0016] 其中,
[0017] 而测量方程(2)可以写成:
[0018] Y(t)=CX(t)+Dsfs(t)  (4)
[0019] 其中 为已知的2乘1阶的常值矩阵;
[0020] 2)引入增广状态向量X,将原来故障下变桨系统模型变换成一种新型的广义系统模型;
[0021] 定义如下的增广矩阵:
[0022]
[0023] 其中,Ip代表维数为p单位矩阵,0k×l代表k行l列的零矩阵,参数δ>0要使得下面条件(6)成立;
[0024]
[0025] 其中rank表示矩阵的秩;
[0026] 最终,可以得到一个增广的广义系统模型(7);
[0027]
[0028] 其中, 是式(5)中所定义的维数为5的广义系统向量, 是维数为2的广义系统测量输出向量, 是式(5)中所定义的维数为5的广义非线性向量,
是式(5)中所定义的5乘5阶的矩阵, 是式(5)中所定义的5乘5阶的参数矩
阵, 是式(5)中所定义的5乘1阶的参数矩阵, 是式(5)中所定义的5乘3阶的
矩阵, 是式(5)中所定义的2乘5阶的参数矩阵, 是式(5)中所定义的维数为3
的广义故障向量;
[0029] 3)基于广义系统模型(7)设计出以下结构的广义滑模估计器(8),实现对桨叶位置状态估计和系统故障重构;
[0030]
[0031] 其中, 是维数为5的中间变量向量, 是维数为5的增广的广义系统状态估计向量, 是式(5)中所定义的5乘2阶的矩阵,
分别是所设计的5乘2阶的微分增益矩阵,5乘2阶的比例增益矩阵以及5乘3阶的滑模增益矩阵, 是维数为3的连续的输入信号
[0032] 3.1)设计微分增益 则可证明矩阵 是非奇异的;
[0033] 3.2)设计比例增益
[0034] 首先,利用如下李雅普诺夫方程(9)求解出矩阵
[0035]
[0036] 其中,μ>0满足
[0037] 进一步,比例增益 可通过式子(10)进行计算;
[0038]
[0039] 3.3)设计滑模增益
[0040] 3.4)令 和 定义如下滑模面:
[0041]
[0042] 其中, 为正定矩阵,且满足以下条件:
[0043]
[0044] 其中, 是由 所确定的3乘2阶的设计矩阵;
[0045] 于是,最终设计的滑模面变为:
[0046]
[0047] 设计的连续输入信号具有如下的滑模形式:
[0048]
[0049] 其中,δ>0满足条件(6),γ>0为设计参数,ο>0是一个小常数,取ο=10-3,α0,α1,β0和d0满足以下条件(15);
[0050]
[0051] 4)基于状态估计设计容错控制方法,实现故障系统的状态稳定;
[0052] 4.1)基于状态估计值,定义一个积分型滑模面(16):
[0053]
[0054] 其中, 和 分别是所设计1乘2阶的参数矩阵;矩阵G要使得GB=1,矩阵K要使得矩阵(A+BK)是霍尔维茨的;
[0055] 4.2)所设计的控制器结构(17)-(18)如下所示:
[0056]
[0057] 且
[0058]
[0059] 其中,
[0060]
[0061] 5)给出估计误差系统和原状态系统稳定的条件,计算参数矩阵,并把各参数矩阵配置给所设计的估计器和控制器;
[0062] 稳定性条件:在滑模控制us(t)和控制输入u(t)下,如果存在2乘2阶的正定矩阵2乘2阶的正定矩阵 和3乘2阶的矩阵 使得下面LMI优化问题(20)成
立:
[0063]
[0064] 其中,υ是固定常数,
[0065] 则可以确保估计误差系统和原状态系统的稳定;
[0066] 通过LMI工具箱的“mincx”求解器,可以方便地求解以上线性矩阵不等式优化问题,计算出所设计的参数矩阵
[0067] 最后,在风电变桨系统中,将LMI优化问题中求解出的参数矩阵 同以上步骤中所设计的微分增益矩阵 比例增益矩阵 滑模增益矩阵 控制器参数矩阵G和K相应配置给所设计的估计器和控制器,实现对桨叶位置的估计和对执行器故障和传感器故障的重构,实现变桨系统稳定。
[0068] 与现有技术相比,本发明具有如下优点:
[0069] 1、首次将广义滑模估计器的故障重构方法应用到机组变桨的故障模型中,有效地实现故障的重构和系统的稳定。
[0070] 2、通过LMI工具箱来求解LMI优化问题,可以容易地得到所设计的估计器增益矩阵,参数设计简单,且易于实现。附图说明
[0071] 图1为一种基于广义滑模估计器的机组变桨容错控制方法的流程图
[0072] 图2为实际桨距角与估计桨距角的曲线图;
[0073] 图3为实际的执行器故障曲线和估计的执行器故障曲线对比图;
[0074] 图4为实际的传感器故障曲线和估计的传感器故障曲线对比图。

具体实施方式

[0075] 下面结合附图和具体实施方式对本发明作进一步详细说明。
[0076] 一种基于广义滑模估计器的机组变桨容错控制方法,其流程如图1所示,该方法包括如下步骤:
[0077] 1)建立同时存在执行器故障和传感器故障下的变桨系统模型;
[0078] 首先,考虑存在执行器故障下的叶轮动态方程:
[0079]
[0080] 再次,考虑存在传感器故障下的桨距角测量方程为:
[0081] y(t)=x(t)+fs(t)  (2)
[0082] 其中,M(·)=JLB+JBI>0;分别代表实际的桨距角、相应的角速度以及相应的角加速
度, 是桨距角的实际测量值, 是推进器、叶片升力、叶片弯曲、振动及不
平衡载荷引起的总力矩, 是控制输入, 是执行器故障函数, 是传感
器故障函数,JLB为空气加速引起的等效惯量力矩,JBI是转子叶片沿其纵轴的全部惯性力矩,kDB是阻尼系数,kRL是摩擦系数;
[0083] 进一步,方程(1)可以写成:
[0084]
[0085] 其中,
[0086] 而测量方程(2)可以写成:
[0087] Y(t)=CX(t)+Dsfs(t)  (4)
[0088] 其中 为已知的2乘1阶的常值矩阵。
[0089] 2)引入增广状态向量 ,将原来故障下变桨系统模型变换成一种新型的广义系统模型;
[0090] 定义如下的增广矩阵:
[0091]
[0092] 其中,Ip代表维数为p单位矩阵,0k×l代表k行l列的零矩阵,参数δ>0要使得下面条件(6)成立;
[0093]
[0094] 其中rank表示矩阵的秩;
[0095] 最终,可以得到一个增广的广义系统模型(7);
[0096]
[0097] 其中, 是式(5)中所定义的维数为5的广义系统向量, 是维数为2的广义系统测量输出向量, 是式(5)中所定义的维数为5的广义非线性向量,
是式(5)中所定义的5乘5阶的矩阵, 是式(5)中所定义的5乘5阶的参数矩
阵, 是式(5)中所定义的5乘1阶的参数矩阵, 是式(5)中所定义的5乘3阶的
矩阵, 是式(5)中所定义的2乘5阶的参数矩阵, 是式(5)中所定义的维数为3
的广义故障向量。
[0098] 3)基于广义系统模型(7)设计出以下结构的广义滑模估计器(8),实现对桨叶位置状态估计和系统故障重构;
[0099]
[0100] 其中, 是维数为5的中间变量向量, 是维数为5的增广的广义系统状态估计向量, 是式(5)中所定义的5乘2阶的矩阵,
分别是所设计的5乘2阶的微分增益矩阵,5乘2阶的比例增益矩阵以及5乘3阶的滑模增益矩阵, 是维数为3的连续的输入信号。
[0101] 3.1)设计微分增益 则可证明矩阵 是非奇异的;
[0102] 3.2)设计比例增益
[0103] 首先,利用如下李雅普诺夫方程(9)求解出矩阵
[0104]
[0105] 其中,μ>0满足
[0106] 进一步,比例增益 可通过式子(10)进行计算;
[0107]
[0108] 3.3)设计滑模增益
[0109] 3.4)令 和 定义如下滑模面:
[0110]
[0111] 其中, 为正定矩阵,且满足以下条件:
[0112]
[0113] 其中, 是由 所确定的3乘2阶的设计矩阵;
[0114] 于是,最终设计的滑模面变为:
[0115]
[0116] 设计的连续输入信号具有如下的滑模形式:
[0117]
[0118] 其中,δ>0满足条件(6),γ>0为设计参数,ο>0是一个小常数,取ο=10-3,α0,α1,β0和d0满足以下条件(15);
[0119]
[0120] 4)基于状态估计设计容错控制方法,实现故障系统的状态稳定;
[0121] 4.1)基于状态估计值,定义一个积分型滑模面(16):
[0122]
[0123] 其中, 和 分别是所设计1乘2阶的参数矩阵;矩阵G要使得GB=1,矩阵K要使得矩阵(A+BK)是霍尔维茨的;
[0124] 4.2)所设计的控制器结构(17)-(18)如下所示:
[0125]
[0126] 且
[0127]
[0128] 其中,
[0129]
[0130] 5)给出估计误差系统和原状态系统稳定的条件,计算参数矩阵,并把各参数矩阵配置给所设计的估计器和控制器;
[0131] 稳定性条件:在滑模控制us(t)和控制输入u(t)下,如果存在2乘2阶的正定矩阵2乘2阶的正定矩阵 和3乘2阶的矩阵 使得下面LMI优化问题(20)
成立:
[0132]
[0133] 其中,υ是固定常数,则可以确保估计误差系统和原状态系统的稳定。
[0134] 通过LMI工具箱的“mincx”求解器,可以方便地求解以上线性矩阵不等式优化问题,计算出所设计的参数矩阵
[0135] 最后,在风电变桨系统中,将LMI优化问题中求解出的参数矩阵 同以上步骤中所设计的微分增益矩阵 比例增益矩阵 滑模增益矩阵 控制器参数矩阵G和K相应配置给所设计的估计器和控制器,实现对桨叶位置的估计和对执行器故障和传感器故障的重构,实现变桨系统稳定。
[0136] 在机组变桨中采用上述设计的广义滑模估计器结构(8)和容错控制率(17)-(18),按照上述的矩阵参数设计步骤,可以实现对桨距角及桨距角速度的估计,执行器和传感器故障的重构以及桨叶位置的稳定。图2表明通过本发明的广义滑模估计器和容错控制方法,桨距角位置不仅能得到估计,而且能实现状态稳定;图3表明执行器故障的估计情况;图4表明传感器故障的估计情况,均表明了本发明的广义滑模估计器具有较好的故障估计效果。
[0137] 最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈