CN114204870B - 一种基于区间ⅱ型t-s模糊模型的风力发电系统积分滑模控制方法 - Google Patents

一种基于区间ⅱ型t-s模糊模型的风力发电系统积分滑模控制方法 Download PDF

Info

Publication number
CN114204870B
CN114204870B CN202111343818.XA CN202111343818A CN114204870B CN 114204870 B CN114204870 B CN 114204870B CN 202111343818 A CN202111343818 A CN 202111343818A CN 114204870 B CN114204870 B CN 114204870B
Authority
CN
China
Prior art keywords
sliding mode
permanent magnet
synchronous motor
magnet synchronous
power generation
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
Application number
CN202111343818.XA
Other languages
English (en)
Other versions
CN114204870A (zh
Inventor
赵孟皕
高忠飞
李猛
陈晨
张茜茹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jinzhou Electric Power Supply Co Of State Grid Liaoning Electric Power Supply Co ltd
State Grid Corp of China SGCC
Original Assignee
Jinzhou Electric Power Supply Co Of State Grid Liaoning Electric Power Supply Co ltd
State Grid Corp of China SGCC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jinzhou Electric Power Supply Co Of State Grid Liaoning Electric Power Supply Co ltd, State Grid Corp of China SGCC filed Critical Jinzhou Electric Power Supply Co Of State Grid Liaoning Electric Power Supply Co ltd
Priority to CN202111343818.XA priority Critical patent/CN114204870B/zh
Publication of CN114204870A publication Critical patent/CN114204870A/zh
Application granted granted Critical
Publication of CN114204870B publication Critical patent/CN114204870B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/22Current control, e.g. using a current control loop
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/0003Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/0003Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • H02P21/0007Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control using sliding mode control
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/0003Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • H02P21/001Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control using fuzzy control
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/28The renewable source being wind energy
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Fuzzy Systems (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

一种基于区间Ⅱ型T‑S模糊模型的风力发电系统积分滑模控制方法,属于风力发电系统稳定控制领域包括如下步骤:步骤1:分析风力发电系统中永磁同步电机工作原理,建立非线性动态控制系统,基于区间Ⅱ型T‑S模糊原理,建立T‑S模糊模型;步骤2:设计积分滑模面,依据李雅普诺夫稳定理论设计系统线性矩阵不等式,求解反馈矩阵,并设计积分滑模控制律;步骤3:根据风力发电系统运行参数,进行仿真分析。本发明利用区间Ⅱ型T‑S模糊原理和积分滑模控制原理相结合来设计控制律,利用区间Ⅱ型T‑S模糊原理实现了用线性系统贴切描述风力发电系统实际工作状态,积分滑模控制原理确保了控制系统的渐近稳定性,实现风力发电系统的稳定控制。

Description

一种基于区间Ⅱ型T-S模糊模型的风力发电系统积分滑模控制方法
技术领域
本发明属于风力发电系统稳定控制领域,特别涉及一种基于区间Ⅱ型T-S模糊模型的风力发电系统积分滑模控制方法。
背景技术
随着世界各国对环保意识的增强,各行各业都在推动低碳行动,大力发展新能源已经成为未来趋势所在。众所周知,我国海上风能资源及其丰富,到目前为止,我国风电装机容量稳居世界之首。随着我国风力发电装机规模逐步扩大,风力发电有望代替火力发电占据我国整体供电的大部分市场。与火力发电相比,风力发电虽然拥有环保、可再生的优点,但由于其受天气因素影响以及并网时依赖逆变器等大量电力电子设备,风力发电具有较大的间歇性与不确定性等问题。由于缺少传统火力发电机组的可提供较大系统惯性的旋转设备,因此,研究风力发电系统的稳定运行控制方法具有重要意义。
风力发电系统的稳定运行是保证电力系统安全稳定的基本要求。随着风力发电装机容量大幅度提升,对风力发电系统的稳定控制提出了更高的要求。目前,国内外科研领域,关于风力发电系统稳定运行技术的研究,主要涉及风力发电系统数学模型的搭建、控制系统控制参数优化以及风力发电系统分布式控制等方面。应用的主要技术包括T-S模糊逼近、神经网络建模、自适应控制参数优化、分布式下垂控制以及PID控制。
传统的T-S模糊逼近方法并不能精确描述风力发电系统的实际运行状态;而神经网络建模往往伴随着复杂的算法,间接提高了稳定控制的复杂程度,在实际应用中并未展现出良好的应用效果;自适应控制参数优化以及PID控制策略局限于参数优化问题,稳定性较差;而分布式下垂控制是通过设定风力发电系统中电气参数之间的下垂关系来模拟传统的发电机组电气参数之间的影响特性,但目前应用最多的先行下垂关系并不能贴切表达电气参数之间的运行特性,因此,分布式下垂控制策略并不能满足风力发电系统稳定运行的更高要求。
发明内容
本发明要解决的技术问题是提供一种基于区间Ⅱ型T-S模糊模型的风力发电系统积分滑模控制方法,该方法利用区间Ⅱ型T-S模糊技术和积分滑模控制技术相结合来设计控制律,利用区间Ⅱ型T-S模糊原理实现了用线性系统贴切描述风力发电系统实际工作状态,积分滑模控制原理确保了控制系统的渐近稳定性,实现风力发电系统的稳定控制。
本发明采用技术方案如下:
一种基于区间Ⅱ型T-S模糊模型的风力发电系统积分滑模控制方法,包括如下步骤:
步骤1:对风力发电系统的工作原理进行分析,根据风力发电系统中的永磁同步电机的实际工作状态建立非线性控制系统,考虑永磁同步电机中电机动摩擦系数的变化,基于区间Ⅱ型T-S模糊原理,选择前件变量并设计隶属度函数,对非线性模型进行线性化处理,建立区间Ⅱ型T-S模糊模型;
步骤1.1:考虑风力发电系统的组成以及永磁同步电机的工作特性,抽象风力发电系统的动态模型,风力发电系统包括风轮机、传动轴、永磁同步电机、逆变器和变压器五部分,其动态模型表达式为:
Figure BDA0003353224870000021
Figure BDA0003353224870000022
Figure BDA0003353224870000023
其中:ωg(t)表示永磁同步电机转速;id(t)和iq(t)表示永磁同步电机的直轴电流、交轴电流;J表示永磁同步电机惯性系数;f(t)表示永磁同步电机摩擦系数;Pn表示永磁同步电机磁极对数;ψf表示永磁同步电机永磁体的磁通量;Rs表示永磁同步电机定子绕组;Ld和Lq分别表示永磁同步电机的直轴电感、交轴电感;Vd(t)和Vq(t)分别表示永磁同步电机的机端电压直轴、交轴分量;t代表时间;
Figure BDA0003353224870000024
T0表示风力发电系统转矩;Cpmax表示风力发电系统最大功率因数;ρa表示空气密度;R表示风轮机叶片长度;λopt表示对应风力发电系统最大功率因数的叶尖速比;
步骤1.2:以永磁同步电机转速ωg(t)、永磁同步电机的直轴电流id(t)、永磁同步电机的交轴电流iq(t)为系统状态变量x(t),以永磁同步电机机端电压直轴、交轴分量Vd(t)、Vq(t)为系统控制输入u(t)建立永磁同步电机非线性动态模型:
Figure BDA0003353224870000025
式中:x(t)=[x1(t),x2(t),x3(t)]T=[ωg(t),iq(t),id(t)]T,u(t)=[Vq(t),Vd(t)]T
Figure BDA0003353224870000026
代表系统状态变量变化函数的微分。
Figure BDA0003353224870000031
其中A∈R3×3、B∈R3×2都是根据风力发电系统、永磁同步电机实际运行状况选取的实际参数归纳的实常数矩阵;
步骤1.3:根据永磁同步电机的非线性动态模型,选取前件变量,规定模糊规则并设计隶属度函数:
在实际工作中永磁同步电机摩擦系数是时刻变化的,因此,令f(t)=T0μ(t),系统状态变量x1(t)∈[μ-d,μ+d],其中μ(t)∈[1,7],d=10,μ(t)代表永磁同步电机摩擦系数变化的函数,d代表永磁同步电机转速的变化范围上限,则可建立以下T-S模糊模型:
模糊规则i:若x1(t)隶属于模糊集合
Figure BDA0003353224870000032
Figure BDA0003353224870000033
式中
Figure BDA0003353224870000034
Figure BDA0003353224870000035
则全局模糊模型可表示为:
Figure BDA0003353224870000036
式中hi(x1(t))代表隶属度函数,h1(x1(t))=0.5-(μ(t)-x1(t))/(2d),h2(x1(t))=0.5+(μ(t)-x1(t))/(2d)。由于在隶属度函数中存在不确定参数μ(t),因此区间Ⅰ型T-S模糊模型并不能贴切描述永磁同步永磁同步电机的原非线性模型。
步骤1.4:利用上下界隶属度函数,引入权重系数,建立区间Ⅱ型T-S模糊模型:
基于区间Ⅰ型T-S模糊模型的隶属度函数,设计以下上下界隶属度函数;
Figure BDA0003353224870000041
Figure BDA0003353224870000042
引入权重系数v1(x1(t)),其上下界分别为
Figure BDA0003353224870000043
Figure BDA0003353224870000044
则对应模糊规则1的隶属度函数可表示为:
Figure BDA0003353224870000045
对应模糊规则2的隶属度函数h2(x1(t))=1-h1(x1(t)),依据全局区间Ⅱ型T-S模糊模型的隶属度函数可实现由线性系统向原非线性系统的高精度逼近。
步骤2:为保证永磁同步电机运行状态稳定,针对系统状态变量设计积分滑模控制律,利用线性矩阵不等式求解积分滑模函数的反馈矩阵,并通过李雅普诺夫理论对积分滑模控制律进行分析,保证当永磁同步电机运行状态偏离预设的滑模函数时,电机运行状态能以较快的速度恢复到滑模函数上并保持稳定;步骤2的具体方法为:
步骤2.1:根据步骤1中建立的永磁同步电机区间Ⅱ型T-S模糊模型,利用积分滑模控制的方法,基于模糊模型状态变量设计积分滑模函数;
定义以下全局积分滑模函数s(t)为:
Figure BDA0003353224870000046
式中,H∈R2×3是给定的常数矩阵,满足HBi(i=1,2)为非奇异矩阵。Ki∈R2×3是积分滑模函数中的反馈矩阵,可对以下步骤中的线性矩阵不等式求解获得。
步骤2.2:基于T-S模糊系统状态变量设计线性矩阵不等式,求解线性矩阵不等式获得积分滑模函数中反馈矩阵Ki,运用等效控制原理确保系统状态运行到滑模函数后的稳定性;
对步骤2.1中的积分滑模函数进行求导得:
Figure BDA0003353224870000051
当系统状态到达滑模函数后应保持稳定,即
Figure BDA0003353224870000052
则系统等效控制律为:
Figure BDA0003353224870000053
将系统等效控制律ueqi(t)带入区间Ⅱ型T-S模糊模型中,原模型可表达为:
Figure BDA0003353224870000054
考虑系统状态稳定性,建立下列李雅普诺夫函数代表系统状态变化趋势:
Figure BDA0003353224870000055
式中,Pi是正定矩阵,对上述李雅普诺夫函数进行求导:
Figure BDA0003353224870000056
Figure BDA0003353224870000057
则表示系统状态在达到滑模函数后是稳定的,即确保
Figure BDA0003353224870000058
令QiPi -1,在上述不等式两边同时乘以矩阵Qi可得
Figure BDA00033532248700000510
令Mi=KiQi,则可获得以下线性矩阵不等式
Figure BDA00033532248700000511
对上述线性矩阵不等式进行求解可得到积分滑模函数中反馈矩阵Ki
步骤2.3:根据滑模控制原理,基于上述步骤所设计的积分滑模函数,设计切换控制项使永磁同步电机状态在偏离滑模函数后可在有限时间内达到滑模函数上,从而保证整个系统渐进稳定性;
设计切换控制律为:
uswi=-β1s(t)-β2sgn(s(t))
其中β1、β2是大于零的常数,β1主要影响系统状态在偏离滑模函数到达滑模函数的时间,β2主要影响系统状态达到滑模函数后的抖振情况,为确保控制性能,控制律参数0≤β1≤50、0≤β2≤1。
系统模糊滑模控制律由等效控制律和切换控制律构成,可表示为:
Figure BDA0003353224870000061
选取李雅普诺夫函数,此函数代表系统状态运行趋势,并证明其导数小于零,使永磁同步电机系统渐进稳定;
取李雅普诺夫函数
Figure BDA0003353224870000062
Figure BDA0003353224870000063
将系统控制律带入系统中,并对积分滑模函数进行求导得:
Figure BDA0003353224870000064
将上式代入
Figure BDA0003353224870000065
Figure BDA0003353224870000066
步骤3:基于步骤1、2的控制律设计过程,对线性矩阵不等式进行求解,获得反馈矩阵Ki,进行仿真分析;
步骤3.1:根据风力发电系统的历史运行数据,将系统参数带入控制系统中,并定义系统状态变量x1(t)变化范围以及电机滑动摩擦系数变化范围;
Pn=8,Rs=1.13Ω,Ld=Lq=0.0027H,ψf=0.16Wb,J=0.005N·m,f(t)∈[0.002,0.014],Cpmax=0.41H,ρa=1.2Kg/m2,R=0.5m,λopt=8.1。
依据上述系统参数可求得相应系统参数矩阵,将系统参数矩阵带入步骤2.2设计的线性矩阵不等式中可求得积分滑模函数中反馈矩阵,
Figure BDA0003353224870000071
选取切换控制律中参数β1=10、β2=0.15,并对风力发电系统进行仿真分析。
本发明的有益效果:
本发明依据风力发电系统中永磁同步电机的工作特性建立非线性控制模型,考虑永磁同步电机中电机动摩擦系数变化,运用区间Ⅱ型T-S模糊原理实现用线性控制系统对非线性控制系统的精准逼近,逼近效果相比于传统的T-S模糊方法更加精确;并采用积分滑模控制策略实现对风力发电系统的稳定控制,相比于传统的滑模控制方法,所采用的积分滑模控制策略具有更好的鲁棒性,提高了控制系统的控制精度,运用线性矩阵不等式求解相应的反馈矩阵,相比于极点配置等方法可得到更可靠的更精确的反馈矩阵,随着风力发电的快速发展,本发明可以确保风力发电系统的稳定运行,对维持风电场安全稳定运行、提高供电质量具有指导意义和推广价值。
附图说明
图1为本发明风力发电系统结构示意图;
图2为本发明电机转速变化仿真图;
图3为本发明电机直轴电流变化仿真图;
图4为本发明电机交轴电流变化仿真图;
图5为本发明风力发电系统系统状态变量变化三维仿真图。
具体实施方式
下面结合附图及实施例对本发明做进一步解释。
本发明为一种基于区间Ⅱ型T-S模糊模型的风力发电系统积分滑模控制方法,针对风力发电系统运行状态,考虑了永磁同步电机中动摩擦系数的变化情况,运用区间Ⅱ型T-S模糊原理实现对非线性系统的精准逼近,设计积分滑模面,并基于李雅普诺夫稳定理论设计系统控制律,实现风力发电系统的稳定控制,包括以下步骤:
步骤1:对风力发电系统的工作原理进行分析,根据风力发电系统中的永磁同步电机的实际工作状态建立非线性控制系统,考虑永磁同步电机中电机动摩擦系数的变化,基于区间Ⅱ型T-S模糊原理,选择前件变量并设计隶属度函数,对非线性模型进行线性化处理,建立区间Ⅱ型T-S模糊模型;
步骤1.1:考虑风力发电系统的组成以及永磁同步电机的工作特性,抽象风力发电系统的动态模型,如图1所示,风力发电系统包括风轮机、传动轴、永磁同步电机、逆变器和变压器五部分,其动态模型表达式为:
Figure BDA0003353224870000081
Figure BDA0003353224870000082
Figure BDA0003353224870000083
其中:ωg(t)表示永磁同步电机转速;id(t)和iq(t)表示永磁同步电机的直轴电流、交轴电流;J表示永磁同步电机惯性系数;f(t)表示永磁同步电机摩擦系数;Pn表示永磁同步电机磁极对数;ψf表示永磁同步电机永磁体的磁通量;Rs表示永磁同步电机定子绕组;Ld和Lq分别表示永磁同步电机的直轴电感、交轴电感;Vd(t)和Vq(t)分别表示永磁同步电机的机端电压直轴、交轴分量;t代表时间;
Figure BDA0003353224870000084
T0表示风力发电系统转矩;Cpmax表示风力发电系统最大功率因数;ρa表示空气密度;R表示风轮机叶片长度;λopt表示对应风力发电系统最大功率因数的叶尖速比;
步骤1.2:以永磁同步电机转速ωg(t)、永磁同步电机的直轴电流id(t)、永磁同步电机的交轴电流iq(t)为系统状态变量x(t),以永磁同步电机机端电压直轴、交轴分量Vd(t)、Vq(t)为系统控制输入u(t)建立永磁同步电机非线性动态模型:
Figure BDA0003353224870000085
式中:x(t)=[x1(t),x2(t),x3(t)]T=[ωg(t),iq(t),id(t)]T,u(t)=[Vq(t),Vd(t)]T
Figure BDA0003353224870000086
代表系统状态变量变化函数的微分。
Figure BDA0003353224870000091
其中
Figure BDA0003353224870000092
都是根据风力发电系统、永磁同步电机实际运行状况选取的实际参数归纳的实常数矩阵;
步骤1.3:根据永磁同步电机的非线性动态模型,选取前件变量,规定模糊规则并设计隶属度函数:
在实际工作中永磁同步电机摩擦系数是时刻变化的,因此,令f(t)=T0μ(t),系统状态变量x1(t)∈[μ-d,μ+d],其中μ(t)∈[1,7],d=10,μ(t)代表永磁同步电机摩擦系数变化的函数,d代表永磁同步电机转速的变化范围上限,则可建立以下T-S模糊模型:
模糊规则i:若x1(t)隶属于模糊集合
Figure BDA0003353224870000093
Figure BDA0003353224870000094
式中
Figure BDA0003353224870000095
Figure BDA0003353224870000096
则全局模糊模型可表示为:
Figure BDA0003353224870000097
式中hi(x1(t))代表隶属度函数,h1(x1(t))=0.5-(μ(t)-x1(t))/(2d),h2(x1(t))=0.5+(μ(t)-x1(t))/(2d)。由于在隶属度函数中存在不确定参数μ(t),因此区间Ⅰ型T-S模糊模型并不能贴切描述永磁同步永磁同步电机的原非线性模型。
步骤1.4:利用上下界隶属度函数,引入权重系数,建立区间Ⅱ型T-S模糊模型:
基于区间Ⅰ型T-S模糊模型的隶属度函数,设计以下上下界隶属度函数;
Figure BDA0003353224870000101
Figure BDA0003353224870000102
引入权重系数v1(x1(t)),其上下界分别为
Figure BDA0003353224870000103
Figure BDA0003353224870000104
则对应模糊规则1的隶属度函数可表示为:
Figure BDA0003353224870000105
对应模糊规则2的隶属度函数h2(x1(t))=1-h1(x1(t)),依据全局区间Ⅱ型T-S模糊模型的隶属度函数可实现由线性系统向原非线性系统的高精度逼近。
步骤2:为保证永磁同步电机运行状态稳定,针对系统状态变量设计积分滑模控制律,利用线性矩阵不等式求解积分滑模函数的反馈矩阵,并通过李雅普诺夫理论对积分滑模控制律进行分析,保证当永磁同步电机运行状态偏离预设的滑模函数时,电机运行状态能以较快的速度恢复到滑模函数上并保持稳定;步骤2的具体方法为:
步骤2.1:根据步骤1中建立的永磁同步电机区间Ⅱ型T-S模糊模型,利用积分滑模控制的方法,基于模糊模型状态变量设计积分滑模函数;
定义以下全局积分滑模函数s(t)为:
Figure BDA0003353224870000106
式中,H∈R2×3是给定的常数矩阵,满足HBi(i=1,2)为非奇异矩阵。Ki∈R2×3是积分滑模函数中的反馈矩阵,可对以下步骤中的线性矩阵不等式求解获得。
步骤2.2:基于T-S模糊系统状态变量设计线性矩阵不等式,求解线性矩阵不等式获得积分滑模函数中反馈矩阵Ki,运用等效控制原理确保系统状态运行到滑模函数后的稳定性;
对步骤2.1中的积分滑模函数进行求导得:
Figure BDA0003353224870000111
当系统状态到达滑模函数后应保持稳定,即
Figure BDA0003353224870000112
则系统等效控制律为:
Figure BDA0003353224870000113
将系统等效控制律ueqi(t)带入区间Ⅱ型T-S模糊模型中,原模型可表达为:
Figure BDA0003353224870000114
考虑系统状态稳定性,建立下列李雅普诺夫函数代表系统状态变化趋势:
Figure BDA0003353224870000115
式中,Pi是正定矩阵,对上述李雅普诺夫函数进行求导:
Figure BDA0003353224870000116
Figure BDA0003353224870000117
则表示系统状态在达到滑模函数后是稳定的,即确保
Figure BDA0003353224870000118
令Qi=Pi -1,在上述不等式两边同时乘以矩阵Qi可得
Figure BDA00033532248700001110
令Mi=KiQi,则可获得以下线性矩阵不等式
Figure BDA00033532248700001111
对上述线性矩阵不等式进行求解可得到积分滑模函数中反馈矩阵Ki
步骤2.3:根据滑模控制原理,基于上述步骤所设计的积分滑模函数,设计切换控制项使永磁同步电机状态在偏离滑模函数后可在有限时间内达到滑模函数上,从而保证整个系统渐进稳定性;
设计切换控制律为:
uswi=-β1s(t)-β2sgn(s(t))
其中β1、β2是大于零的常数,β1主要影响系统状态在偏离滑模函数到达滑模函数的时间,β2主要影响系统状态达到滑模函数后的抖振情况,为确保控制性能,控制律参数0≤β1≤50、0≤β2≤1。
系统模糊滑模控制律由等效控制律和切换控制律构成,可表示为:
Figure BDA0003353224870000121
选取李雅普诺夫函数,此函数代表系统状态运行趋势,并证明其导数小于零,使永磁同步电机系统渐进稳定;
取李雅普诺夫函数
Figure BDA0003353224870000122
Figure BDA0003353224870000123
将系统控制律带入系统中,并对积分滑模函数进行求导得:
Figure BDA0003353224870000124
将上式代入
Figure BDA0003353224870000125
Figure BDA0003353224870000126
步骤3:基于步骤1、2的控制律设计过程,对线性矩阵不等式进行求解,获得反馈矩阵Ki,进行仿真分析;
步骤3.1:根据风力发电系统的历史运行数据,将系统参数带入控制系统中,并定义系统状态变量x1(t)变化范围以及电机滑动摩擦系数变化范围;
Pn=8,Rs=1.13Ω,Ld=Lq=0.0027H,ψf=0.16Wb,J=0.005N·m,f(t)∈[0.002,0.014],Cpmax=0.41H,ρa=1.2Kg/m2,R=0.5m,λopt=8.1。
依据上述系统参数可求得相应系统参数矩阵,将系统参数矩阵带入步骤2.2设计的线性矩阵不等式中可求得积分滑模函数中反馈矩阵,
Figure BDA0003353224870000131
选取切换控制律中参数β1=10、β2=0.15,并对风力发电系统进行仿真分析。
首先,输入风力发电系统的系统参数,并构成相应的系统矩阵;其次利用线性矩阵不等式工具箱对线性矩阵不等式进行求解,解得相应的积分滑模函数中的反馈矩阵;最后,基于系统矩阵、反馈矩阵描述系统控制律和系统状态微分方程,并对微分方程进行求解,得到风力发电系统状态变量变化曲线。
图2为永磁同步电机转速变化仿真图如图2所示,永磁同步电机直轴电流变化仿真图如图3所示,永磁同步电机交轴电流变化仿真图如图4所示,风力发电系统系统状态变量变化三维仿真图如图5所示;根据图2中容易看出,永磁同步电机转速变化曲线在10秒内收敛至零,充分体现了积分滑模控制策略良好的稳定性;在图3与图4中,永磁同步电机电流的直轴分量、交轴分量呈现出了不同的变化趋势,电流的直轴分量相比于交轴分量变化幅度更小,收敛的速度更快,永磁同步电机电流的交轴分量也可以在10秒内收敛至零;分析图5,系统状态变量的运动轨迹都呈现出收敛至滑模函数的趋势,并都在10秒内运动到滑模函数上,这充分体现了积分滑模控制策略的良好控制性能以及区间Ⅱ型T-S模糊原理的有效性,对于永磁同步电机的稳定运行具有重要意义。
以上仅为本发明的具体实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种基于区间Ⅱ型T-S模糊模型的风力发电系统积分滑模控制方法,其特征在于:包括如下步骤:
步骤1:对风力发电系统的工作原理进行分析,根据风力发电系统中的永磁同步电机的实际工作状态建立非线性控制系统,考虑永磁同步电机中电机动摩擦系数的变化,基于区间Ⅱ型T-S模糊原理,选择前件变量并设计隶属度函数,对非线性模型进行线性化处理,建立区间Ⅱ型T-S模糊模型;
步骤1.1:考虑风力发电系统的组成以及永磁同步电机的工作特性,抽象风力发电系统的动态模型,风力发电系统包括风轮机、传动轴、永磁同步电机、逆变器和变压器五部分,其动态模型表达式为:
Figure FDA0003353224860000011
Figure FDA0003353224860000012
Figure FDA0003353224860000013
式中:ωg(t)表示永磁同步电机转速;id(t)和iq(t)表示永磁同步电机的直轴电流、交轴电流;J表示永磁同步电机惯性系数;f(t)表示永磁同步电机摩擦系数;Pn表示永磁同步电机磁极对数;ψf表示永磁同步电机永磁体的磁通量;Rs表示永磁同步电机定子绕组;Ld和Lq分别表示永磁同步电机的直轴电感、交轴电感;Vd(t)和Vq(t)分别表示永磁同步电机的机端电压直轴、交轴分量;t代表时间;
其中,风力发电系统转矩T0的表达式:
Figure FDA0003353224860000014
式中:Cpmax表示风力发电系统最大功率因数;ρa表示空气密度;R表示风轮机叶片长度;λopt表示对应风力发电系统最大功率因数的叶尖速比;
步骤1.2:以永磁同步电机转速ωg(t)、永磁同步电机的直轴电流id(t)、永磁同步电机的交轴电流iq(t)为系统状态变量x(t),以永磁同步电机机端电压直轴、交轴分量Vd(t)、Vq(t)为系统控制输入u(t),通过系统状态变量变化函数的微分
Figure FDA0003353224860000015
形式,建立永磁同步电机非线性动态模型:
Figure FDA0003353224860000021
式中:x(t)=[x1(t),x2(t),x3(t)]T=[ωg(t),iq(t),id(t)]T,u(t)=[Vq(t),Vd(t)]T
Figure FDA0003353224860000022
其中A∈R3×3、B∈R3×2都是根据风力发电系统、永磁同步电机实际运行状况选取的实际参数归纳的实常数矩阵;
步骤1.3:根据永磁同步电机的非线性动态模型,选取前件变量,规定模糊规则并设计隶属度函数:
在实际工作中永磁同步电机摩擦系数是时刻变化的,因此,令f(t)=T0μ(t),系统状态变量x1(t)∈[μ-d,μ+d],其中μ(t)∈[1,7],d=10,μ(t)代表永磁同步电机摩擦系数变化的函数,d代表永磁同步电机转速的变化范围上限,则可建立以下T-S模糊模型:
模糊规则i:若x1(t)隶属于模糊集合
Figure FDA0003353224860000023
Figure FDA0003353224860000024
式中
Figure FDA0003353224860000025
Figure FDA0003353224860000026
则全局模糊模型可表示为:
Figure FDA0003353224860000031
式中hi(x1(t))代表隶属度函数,h1(x1(t))=0.5-(μ(t)-x1(t))/(2d),h2(x1(t))=0.5+(μ(t)-x1(t))/(2d);由于在隶属度函数中存在不确定参数μ(t),因此区间Ⅰ型T-S模糊模型并不能贴切描述永磁同步永磁同步电机的原非线性模型;
步骤1.4:利用上下界隶属度函数,引入权重系数,建立区间Ⅱ型T-S模糊模型:
基于区间Ⅰ型T-S模糊模型的隶属度函数,设计以下上、下界隶属度函数;
Figure FDA0003353224860000032
Figure FDA0003353224860000033
引入权重系数v1(x1(t)),其上下界分别为
Figure FDA0003353224860000034
Figure FDA0003353224860000035
则对应模糊规则1的隶属度函数可表示为:
Figure FDA0003353224860000036
对应模糊规则2的隶属度函数h2(x1(t))=1-h1(x1(t)),依据全局区间Ⅱ型T-S模糊模型的隶属度函数可实现由线性系统向原非线性系统的高精度逼近;
步骤2:为保证永磁同步电机运行状态稳定,针对系统状态变量设计积分滑模控制律,利用线性矩阵不等式求解积分滑模函数的反馈矩阵,并通过李雅普诺夫理论对积分滑模控制律进行分析,保证当永磁同步电机运行状态偏离预设的滑模函数时,电机运行状态能以较快的速度恢复到滑模函数上并保持稳定;步骤2的具体方法为:
步骤2.1:根据步骤1中建立的永磁同步电机区间Ⅱ型T-S模糊模型,利用积分滑模控制的方法,基于模糊模型状态变量设计积分滑模函数;
定义以下全局积分滑模函数s(t)为:
Figure FDA0003353224860000037
式中,H∈R2×3是给定的常数矩阵,满足HBi(i=1,2)为非奇异矩阵;Ki∈R2×3是积分滑模函数中的反馈矩阵,可对以下步骤中的线性矩阵不等式求解获得;
步骤2.2:基于T-S模糊系统状态变量设计线性矩阵不等式,求解线性矩阵不等式获得积分滑模函数中反馈矩阵Ki,运用等效控制原理确保系统状态运行到滑模函数后的稳定性;
对步骤2.1中的积分滑模函数进行求导得:
Figure FDA0003353224860000041
当系统状态到达滑模函数后应保持稳定,即
Figure FDA0003353224860000042
则系统等效控制律为:
Figure FDA0003353224860000043
将系统等效控制律ueqi(t)带入区间Ⅱ型T-S模糊模型中,原模型可表达为:
Figure FDA0003353224860000044
考虑系统状态稳定性,建立下列李雅普诺夫函数代表系统状态变化趋势:
Figure FDA0003353224860000045
式中,Pi是正定矩阵,对上述李雅普诺夫函数进行求导:
Figure FDA0003353224860000046
Figure FDA0003353224860000047
则表示系统状态在达到滑模函数后是稳定的,即确保
Figure FDA0003353224860000048
令Qi=Pi -1,在上述不等式两边同时乘以矩阵Qi可得
Figure FDA0003353224860000049
令Mi=KiQi,则可获得以下线性矩阵不等式
Figure FDA00033532248600000410
对上述线性矩阵不等式进行求解可得到积分滑模函数中反馈矩阵Ki
步骤2.3:根据滑模控制原理,基于上述步骤所设计的积分滑模函数,设计切换控制项使永磁同步电机状态在偏离滑模函数后可在有限时间内达到滑模函数上,从而保证整个系统渐进稳定性;
设计切换控制律uswi为:
uswi=-β1s(t)-β2sgn(s(t))
其中β1、β2是大于零的常数,β1主要影响系统状态在偏离滑模函数到达滑模函数的时间,β2主要影响系统状态达到滑模函数后的抖振情况,为确保控制性能,控制律参数0≤β1≤50、0≤β2≤1;
系统模糊滑模控制律由等效控制律和切换控制律构成,可表示为:
Figure FDA0003353224860000051
选取李雅普诺夫函数,此函数代表系统状态运行趋势,并证明其导数小于零,使永磁同步电机系统渐进稳定;
取李雅普诺夫函数
Figure FDA0003353224860000052
Figure FDA0003353224860000053
将系统控制律带入系统中,并对积分滑模函数进行求导得:
Figure FDA0003353224860000054
将上式代入
Figure FDA0003353224860000055
Figure FDA0003353224860000056
步骤3:基于步骤1、2的控制律设计过程,对线性矩阵不等式进行求解,获得反馈矩阵Ki,进行仿真分析。
2.根据权利要求1所述的基于区间Ⅱ型T-S模糊模型的风力发电系统积分滑模控制方法,其特征是:
进行仿真分析具体过程如下:
根据风力发电系统的历史运行数据,将系统参数带入控制系统中,并定义系统状态变量x1(t)变化范围以及电机滑动摩擦系数变化范围;
依据上述系统参数可求得相应系统参数矩阵,将系统参数矩阵带入步骤2.2设计的线性矩阵不等式中可求得积分滑模函数中反馈矩阵Ki
选取切换控制律中参数βi,并对风力发电系统进行仿真分析。
3.根据权利要求2所述的基于区间Ⅱ型T-S模糊模型的风力发电系统积分滑模控制方法,其特征是:所述切换控制律中参数为β12其中β1=10、β2=0.15。
CN202111343818.XA 2021-11-14 2021-11-14 一种基于区间ⅱ型t-s模糊模型的风力发电系统积分滑模控制方法 Active CN114204870B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111343818.XA CN114204870B (zh) 2021-11-14 2021-11-14 一种基于区间ⅱ型t-s模糊模型的风力发电系统积分滑模控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111343818.XA CN114204870B (zh) 2021-11-14 2021-11-14 一种基于区间ⅱ型t-s模糊模型的风力发电系统积分滑模控制方法

Publications (2)

Publication Number Publication Date
CN114204870A CN114204870A (zh) 2022-03-18
CN114204870B true CN114204870B (zh) 2023-05-09

Family

ID=80647658

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111343818.XA Active CN114204870B (zh) 2021-11-14 2021-11-14 一种基于区间ⅱ型t-s模糊模型的风力发电系统积分滑模控制方法

Country Status (1)

Country Link
CN (1) CN114204870B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110566403A (zh) * 2019-08-08 2019-12-13 天津科技大学 一种风力发电t-s模糊鲁棒调度容错控制方法
CN112186789A (zh) * 2020-09-24 2021-01-05 沈阳工程学院 一种电动汽车参与微电网负荷调频的滑模控制方法
CN113534665A (zh) * 2021-07-23 2021-10-22 杭州电子科技大学 基于区间二型t-s模型的有限时间稳定的滑模控制方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103032265B (zh) * 2012-12-12 2014-11-05 天津市电力公司 一种基于极值搜索的风电机组最大出力追踪控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110566403A (zh) * 2019-08-08 2019-12-13 天津科技大学 一种风力发电t-s模糊鲁棒调度容错控制方法
CN112186789A (zh) * 2020-09-24 2021-01-05 沈阳工程学院 一种电动汽车参与微电网负荷调频的滑模控制方法
CN113534665A (zh) * 2021-07-23 2021-10-22 杭州电子科技大学 基于区间二型t-s模型的有限时间稳定的滑模控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
谭川 ; 刘欣然 ; .基于滑模模糊方法的变速风电系统的最大风能捕捉控制器设计.电力科学与工程.2010,(第02期),全文. *

Also Published As

Publication number Publication date
CN114204870A (zh) 2022-03-18

Similar Documents

Publication Publication Date Title
Mani et al. Adaptive fractional fuzzy integral sliding mode control for PMSM model
Naik et al. Design and implementation of interval type-2 fuzzy logic-PI based adaptive controller for DFIG based wind energy system
Bouderbala et al. Direct and indirect vector control of a doubly fed induction generator based in a wind energy conversion system
CN106788052B (zh) 基于观测器的异步电机命令滤波误差补偿模糊控制方法
Mansour et al. A grid-connected variable-speed wind generator driving a fuzzy-controlled PMSG and associated to a flywheel energy storage system
CN105179164A (zh) 基于t-s模糊模型的风能转换系统滑模控制方法及装置
Laghridat et al. Comparative Analysis between PI and Linear‐ADRC Control of a Grid Connected Variable Speed Wind Energy Conversion System Based on a Squirrel Cage Induction Generator
Barambones Robust wind speed estimation and control of variable speed wind turbines
CN107276471B (zh) 一种基于状态受限的异步电机模糊位置跟踪控制方法
Jabal Laafou et al. Dynamic Modeling and Improved Control of a Grid‐Connected DFIG Used in Wind Energy Conversion Systems
González-Hernández et al. A novel MPPT PI discrete reverse-acting controller for a wind energy conversion system
Abolvafaei et al. Two novel approaches to capture the maximum power from variable speed wind turbines using optimal fractional high-order fast terminal sliding mode control
Zhang et al. Enhanced wind turbine maximum wind-energy capture based on the inverse-system method
CN114204870B (zh) 一种基于区间ⅱ型t-s模糊模型的风力发电系统积分滑模控制方法
Yu et al. Anti-saturation coordination control of permanent magnet synchronous wind power system
Wang et al. Research on speed control of permanent magnet synchronous motor based on RBF neural network tuning PID
Huang et al. Error-based active disturbance rejection control for wind turbine output power regulation
Kumari et al. A comprehensive review of traditional and smart MPPT techniques in PMSG based Wind energy conversion system
Doumi et al. Robust fuzzy gains scheduling of RST controller for a WECS based on a doubly-fed induction generator
Kumar et al. ANN Control for Improved Performance of Wind Energy System Connected to Grid
Laafou et al. Proposed robust ADRC control of a DFIG used in wind power production
BARRA et al. Sensorless speed and reactive power control of a DFIG-wind turbine
Jiang et al. No-load cutting-in control of the doubly fed induction generator based on grey prediction PI control
Abderrahim et al. Power optimisation of wind energy conversion system using TS fuzzy approach
Liu et al. Current control with active disturbance rejection control of immersion and invariance for the grid-connected inverter

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
GR01 Patent grant
GR01 Patent grant