欢迎来到知嘟嘟! 联系电话:13095918853 卖家免费入驻,海量在线求购! 卖家免费入驻,海量在线求购!
知嘟嘟
我要发布
联系电话:13095918853
知嘟嘟经纪人
收藏
专利号: 2019103514866
申请人: 闽江学院
专利类型:发明专利
专利状态:已下证
专利领域: 控制;调节
更新日期:2024-02-23
缴费截止日期: 暂无
价格&联系人
年费信息
委托购买

摘要:

权利要求书:

1.一种风机发电偏航控制系统最大功率追踪方法,其特征在于,包括以下步骤:步骤S1:搭建风机发电偏航控制系统100,所述风机发电偏航控制系统包括风机发电机组10和最大功率追踪滑膜输出反馈控制器20;

步骤S2:根据物理学原理,建立具有最大功率追踪的风机偏航非线性系统的数学模型,并表达为奇异II型模糊系统模型;

步骤S3:基于所述数学模型,设计有限时间滑膜输出反馈控制器;

步骤S4:求出闭环控制系统到达滑膜面的时间T*;

*

步骤S5:求出闭环控制系统在T时刻的最大发电功率参考角速度误差值。

2.根据权利要求1所述的一种风机发电偏航控制系统最大功率追踪方法,其特征在于,步骤S2中,建立具有最大功率追踪的风机偏航非线性系统的数学模型并表达为奇异II型模糊系统模型,包括如下步骤:步骤S21:考虑风机转子动态模型如下:式中,Pa表示气动功率,参数ρ、R、Cp、β和v分别表示空气密度、转子半径、功率系数、桨叶角和风速;叶尖速比λ定义如下:式中,参数ωt表示转子速度;

转矩Ta的动力学方程如下:

转子侧的惯性Jr动态方程如下:式中,Br是转子的阻尼,低速转子转矩Tls被定义为:Tls=Kls(θt-θls)+Kls(ωt-ωls)                  (5)式中,Kls表示低速侧转轴的刚度,θt表示转子侧角偏差,θls表示变速箱侧角偏差,ωt表示转子侧角速度,ωls表示变速箱侧角速度,那么有:式中,参数Jg和Tem分别是发电机的惯性和电磁转矩,Ths表示高速侧的转子转矩,Bg表示发电机外部阻尼,ωg表示发电机角速度;

进一步定义减速比ng,得到:根据空间矢量的表达,定子电压在旋转d-q轴的参考框架下可以表达为:式中,ud和uq是d轴和q轴的定子电压,φd和φq是d轴和q轴的定子电流,Rs是定子绕组阻抗,M是级数,ωg是定子的电角速度,Ψq和Ψd是d轴和q轴的磁通,表达为:式中,Ld和Lq是d轴和q轴的的电感,Ψm是磁链;

随着公式(1)-(9),得到:式中,

步骤S22:功率系数曲线Cp(λ,β)具有最大值,其符合如下的情形:式中,

定义 并且随着公式(10)-(13)得到:步骤S23:定义x(t)=[e ωg φd φq Tem]T,并且考虑在状态变量e上存在扰动ω(t);

选择 和 作为模糊前件变量;因此,具有最大功率追踪的风机偏航非线性系统的数学模型,表达为以下的奇异II型模糊系统模型:y(t)=Cx(t)                        (15)式中, 是中立II型的模糊

集,r模糊规则的数量,Al和Bl是对前件变量 线性化后系统矩阵参数,表达如下:

3.根据权利要求1所述的一种风机发电偏航控制系统最大功率追踪方法,其特征在于,步骤S3中,所述有限时间滑膜输出反馈控制器具有以下形式:u(t)=ub(t)+uc(t)                      (17)式中,

Kl是模

糊控制器增益, {l,p}代表模糊集,符号T是预先给定的有限时间周期,αmin(*)代表矩阵*的最小的秩,‖*‖代表矩阵*的范数;

sgn(*)是符号函数,具体的定义如下:式中,s(t)是积分滑膜面函数,定义如下:式中,G是一个给定的矩阵使得GBl是正定的矩阵。

4.根据权利要求1所述的一种风机发电偏航控制系统最大功率追踪方法,其特征在于,步骤S4中,所述闭环控制系统到达滑膜面时间T*,由以下步骤求出:步骤S41:考虑以下函数求导后得到:定义如下方程:

得到:

随着公式,得到:

步骤S42:随着公式(20)-(23)得到:根据公式(23)和(24),得到:随着公式 闭环控制系统到达滑膜面时间T*获得如下:T*≤T             (26)。

5.根据权利要求1所述一种风机发电偏航控制系统最大功率追踪方法,其特征在于,步骤S5中,所述闭环控制系统在T*时刻的最大发电功率参考角速度误差值,由以下步骤求出:步骤S51:将控制器(17)代入(15),获得如下的闭环控制系统:式中,

步骤S52:考虑如下的李雅普诺夫函数:其中ETP=PTE≥0;

沿着系统(27)的轨迹,得到:式中,

引入一个辅助的函数J(t),如下:其中,η是正标量;

随着公式(29)和(30),得到:式中,

以下的不等式成立会使得公式(31)的J(t)<0:Φll<0,1≤l≤r                (33)Φlp+Φpl<0,1≤l

在不等式(35)左乘右乘e-ηt,并且0到t进行积分t∈[0,T*],那么很容易得到:此外,由公式(28)得到:

式中,

步骤S53:将矩阵P分割如下:式中, 并且P3是一个标量,得到ETP=PTE≥0;

进一步分割状态变量空间x(t)如下:式中,

随着公式(38)和(40),有:引入矩阵 并且进一步定义:

基于公式(38)-(42)得到:步骤S54:由于追踪误差的角速度定义如下: 为了使得ε在有限的时间内收敛,将矩阵P1定义如下:式中,P1(1)是一个标量,矩阵 和P1(1)>>‖P1(2)‖随着公式(42)-(44)得到: