CN108488036A - 基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法 - Google Patents
基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法 Download PDFInfo
- Publication number
- CN108488036A CN108488036A CN201810416971.2A CN201810416971A CN108488036A CN 108488036 A CN108488036 A CN 108488036A CN 201810416971 A CN201810416971 A CN 201810416971A CN 108488036 A CN108488036 A CN 108488036A
- Authority
- CN
- China
- Prior art keywords
- suspension
- formula
- air gap
- model
- levitating electromagnet
- 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.)
- Granted
Links
- 239000000725 suspension Substances 0.000 title claims abstract description 167
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000005611 electricity Effects 0.000 title claims abstract description 22
- 238000004804 winding Methods 0.000 claims abstract description 73
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000005259 measurement Methods 0.000 claims abstract description 7
- 238000013461 design Methods 0.000 claims abstract description 4
- 238000011217 control strategy Methods 0.000 claims abstract description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 7
- 238000007667 floating Methods 0.000 claims description 7
- 238000013178 mathematical model Methods 0.000 claims description 7
- 230000005484 gravity Effects 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 4
- 230000035699 permeability Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000009897 systematic effect Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000004870 electrical engineering Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 229910052742 iron Inorganic materials 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/0204—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor for orientation in relation to wind direction
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/30—Control parameters, e.g. input parameters
- F05B2270/329—Azimuth or yaw angle
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/60—Control system actuates through
- F05B2270/602—Control system actuates through electrical actuators
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/70—Type of control algorithm
- F05B2270/706—Type of control algorithm proportional-integral-differential
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
- Control Of Vehicles With Linear Motors And Vehicles That Are Magnetically Levitated (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明涉及一种基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法,属电气工程技术领域。该方法通过模型预测控制、PID、模型失配器相结合,解决磁悬浮偏航系统平衡点附近因线性化处理造成的模型失配带来的偏差问题,对平衡点处的悬浮进行实时平稳控制:根据线性化模型采用模型预测控制策略,设计悬浮气隙控制器,实现外环悬浮气隙控制;将悬浮气隙控制器的输出与模型失配器的输出相加,减去悬浮电磁铁绕组的电流测量值,经PID控制器,控制绕组电流,实现内环电流跟踪控制;将悬浮气隙偏差和绕组电流偏差分别乘以各自的调节参数后相加构成模型失配器。本发明控制精度高,能有效抑制外界扰动的影响,确保平衡点附近悬浮快速稳定。
Description
技术领域
本发明涉及一种控制方法,尤其是一种基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法,属于电气工程技术领域。
背景技术
偏航系统是水平轴风力发电机必不可少的重要组成部分。传统的偏航装置均为齿轮驱动,摩擦损耗大、故障率高,对风不够精准,维护不便、停电时间长,降低了发电效率。
风电磁悬浮偏航系统采用磁悬浮驱动技术,取代传统的齿轮驱动技术,具有对风精度高、无需润滑、结构简单、维护方便、停电时间短、运行维护费用低等优势。通过电磁悬浮,一方面机舱处于悬浮状态,可实现精确对风;另一方面,简化了偏航系统结构,使维修简便,可大大缩短停机时间。
风电磁悬浮偏航系统的工作原理是:当风向改变时,首先使其悬浮电磁铁通入直流电,实现悬浮,到达悬浮平衡点后,使其定子通入交流电,悬浮电磁铁开始旋转,直至到达对风位置。在旋转过程中,一方面要实施悬浮控制,使其处于平衡点,另一方面要控制其转速,实现稳定旋转,因此实现平衡点处的稳定悬浮至关重要。
现有文献大多采用在悬浮平衡点附近线性化的方法,忽略了高阶项,这会导致模型失配,进而降低了系统鲁棒性。为此有的文献基于均方差的方法使用卡尔曼滤波状态观测器解决无约束条件下模型预测控制系统中的模型失配问题,有的利用遗传算法解决模型失配对模型预测控制的影响,或者采用无模型自适应控制实现算法的智能推理,以抑制时变参数、模型失配等问题,也有人利用扰动观测器在输出预测中加入干扰项,用来校正由干扰和不确定性因素引起的偏差。但上述方法过程均较为复杂。
发明内容
本发明的主要目的在于:针对现有技术的不足和空白,本发明提供一种基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法,采用结构较为简单的模型失配补偿器,通过模型预测控制(MPC)、PID算法、模型失配器相结合,解决模型失配及不确定干扰给系统带来的偏差问题,实现风电磁悬浮偏航系统在平衡点处的稳定悬浮。因此,本发明的控制方法是一种加模型失配器的MPC-PID方法。
为了达到以上目的,本发明所述风电磁悬浮偏航系统包括悬浮系统、驱动系统,所述悬浮系统由悬浮电磁铁、定子、气隙传感器、悬浮架和悬浮变流器等组成;所述悬浮电磁铁包括铁心和绕组,所述绕组为直流励磁绕组,所述绕组与所述悬浮变流器连接。
本发明基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法,包括以下步骤:
步骤1,建立磁悬浮偏航系统的悬浮动态数学模型,建模过程如下:
所述悬浮电磁铁绕组通电后将产生向上的轴向悬浮力F:
式中,i(t)为所述悬浮电磁铁绕组的电流;δ(t)为所述悬浮电磁铁和所述定子之间的悬浮气隙;k1=μ0N2S/4,其中,μ0为真空磁导率,N为所述悬浮电磁铁绕组的匝数,S为所述悬浮电磁铁铁心的磁极表面有效面积;
根据牛顿第二定律,所述悬浮系统在垂直方向上的力学方程为:
式中,m为悬浮物质量,g为重力加速度;fd(t)为外界扰动力,为悬浮气隙δ(t)对时间t的二阶导数;
所述悬浮电磁铁的电压方程为:
式中,u(t)为所述悬浮电磁铁绕组的电压;R为所述悬浮电磁铁绕组的电阻;ψ(t)为气隙磁链;L为所述悬浮电磁铁绕组的气隙电感,且有L=2k1/δ(t);
综上,可得磁悬浮偏航系统的动态悬浮数学模型:
平衡条件:F(i0,δ0)=k1(i0/δ0)2=mg
式中,i0、δ0分别为稳定悬浮时,即在平衡点处所述悬浮电磁铁绕组的电流和悬浮气隙;
步骤2,建立磁悬浮偏航系统的悬浮系统线性化模型,建模过程如下:
将式(1)在平衡点(i0,δ0)附近进行线性化处理,用泰勒级数公式展开,忽略高阶项,可得:
式中,ΔF=F(i0,δ0)-F;Δδ(t)=δ0-δ(t);Δi(t)=i0-i(t);
类似地,将式(2)在平衡点(i0,δ0)附近进行线性化处理,用泰勒级数公式展开,忽略高阶项,可得:
式中,Δu(t)=u0-u(t),u0为平衡点处稳定悬浮时悬浮电磁铁的电压;
将式(4)、式(5)代入式(3),同时令则可得所述悬浮系统的线性化模型:
式中,kc=2k1i0 2/δ0 3;ki=kcδ0/i0;kf为常数;L0=2k1/δ0;
步骤3,采用模型预测控制策略,设计悬浮气隙控制器,实现外环悬浮气隙控制;具体方法如下:
A)令x1=Δi(t),x2=Δδ(t),代入式(6),整理可得悬浮系统状态空间方程为:
式中,x=[x1,x2,x3]T为状态变量,um=Δu(t)为输入变量,
y1为所述悬浮系统的输出量;
B)根据模型预测理论,在x可测时,式(7)的预测模型为:
x(k+1)=Ax(k)+bum(k) (8)
C)设从k时刻输入发生M步变化,而后保持不变,则由式(8)可以预测出um(k),um(k+1),um(k+2),…,um(k+M-1)作用在未来P(P>=M)个时刻的系统状态,即:
X(k)=Fx(k)+GU(k) (9)
式中,X(k)为状态矩阵,X(k)∈Rn实时可测;n为x(k)的维数;U(k)为系统输入矩阵;F、G分别是系统状态转移矩阵、输入矩阵,且有:
其中:P为预测时域长度;M为控制时域长度;
D)选取所述悬浮系统的优化性能指标为:
式中,Q为误差加权系数;R为控制加权系数;
在不考虑约束时,结合状态预测模型式(8),求出式(10)的最优解为:
U(k)=-(GTQG+R)-1GTQFx(k) (11)
E)由此可得输入变量控制律为:
式中,k2为反馈增益,有:
则式(12)构成了悬浮气隙控制器。
步骤4,将步骤3所述悬浮气隙控制器的输出y1=x1=Δi(t)作为所述悬浮电磁铁绕组的电流变化量参考值Δiref(k),据此求得所述悬浮电磁铁绕组的电流参考值iref(k)=i0+Δiref(k),将此值iref(k)与模型失配器的输出相加,两者之和与所述悬浮变流器输出的电流测量值之差,经PID控制器,送入PWM模块,产生驱动信号,控制所述悬浮变流器的输出电压和电流,实现内环电流跟踪控制;
所述模型失配补偿器按下述方式设计:
令模型失配补偿律为:
y2=C1eδ+C2ei (13)
式中,y2为模型失配补偿器输出的补偿电流;eδ=δ(t)-δl(t);ei=i(t)-il(t),其中,δ(t)、i(t)分别为悬浮气隙测量值、所述悬浮电磁铁绕组电流测量值;δl(t)、il(t)分别为线性化模型式(6)输出的悬浮气隙和所述悬浮电磁铁绕组励磁电流;C1是气隙偏差调节参数;C2是电流偏差调节参数,则式(13)构成了模型失配补偿器。
步骤5,将所述悬浮变流器的输出电压送入线性化模型式(6),线性化模型式(6)输出悬浮气隙δl(t)和所述悬浮电磁铁绕组的电流il(t);所述悬浮气隙δl(t)减去悬浮气隙测量值δ(t),得到悬浮气隙偏差eδ,所述悬浮电磁铁绕组的电流il(t)减去所述悬浮电磁铁绕组的电流测量值i(t),得到绕组电流偏差ei,将所述偏差eδ、ei送入模型失配补偿器式(13),得到补偿电流y2,返回步骤4,实现平衡点附近稳定悬浮。
本发明的有益效果是:本发明采用结构简单的模型失配补偿器,与模型预测控制(MPC)、PID算法相结合,能够使系统在较短时间内达到稳定,且控制精度更高,有效解决了模型失配及不确定干扰给系统带来的偏差问题,确保在平衡点附近整个悬浮过程系统性能实时最优,从而实现风电磁悬浮偏航系统在平衡点附近稳定悬浮。
附图说明
图1为本发明所述风电磁悬浮偏航系统结构示意图。
图2为本发明所述磁悬浮偏航系统悬浮受力分析示意图。
图3为本发明悬浮控制系统结构框图。
图4为不加模型失配器的控制方法的悬浮控制系统结构框图。
图5为本发明与不加模型失配器的控制方法的悬浮气隙对比仿真曲线图。
图6为本发明与不加模型失配器的控制方法的悬浮电磁铁绕组电流对比仿真曲线图。
图中标号:1-悬浮电磁铁,2-定子,3-气隙传感器,4-悬浮架,5-导向轴承,6-悬浮变流器,11-悬浮电磁铁绕组,12-悬浮电磁铁铁心
具体实施方式
下面结合附图,对本发明作进一步详细说明。
如图1、图3所示,本发明用于一种风电磁悬浮偏航系统,包括悬浮系统、驱动系统,悬浮系统由悬浮电磁铁1、定子2、气隙传感器3、悬浮架4、导向轴承5和悬浮变流器6等组成;气隙传感器3用于检测悬浮电磁铁1与定子2之间的悬浮气隙δ(t)。悬浮电磁铁1、气隙传感器3、悬浮架4、机舱等悬浮部分统称为悬浮物。
如图2、图3所示,悬浮电磁铁1包括绕组11和铁心12,悬浮电磁铁绕组11为直流励磁绕组,与悬浮变流器6连接;悬浮变流器6为直流斩波器。
该风电磁悬浮偏航系统的工作原理是:
如图1、图2所示,悬浮电磁铁1与定子2之间的气隙长度为δ(t),当风向变化需要偏航时(或需要解缆时),悬浮变流器6给悬浮电磁铁绕组11通入直流电流,将会有一穿过悬浮电磁铁1与定子2之间的气隙、同时交链定子2与悬浮电磁铁1的主磁场产生,定子2被磁化,产生垂直向上的轴向磁吸力F,悬浮电磁铁1将向上运动,通过悬浮架4使整个机舱向上悬浮,当气隙传感器3检测到悬浮气隙δ(t)达到设定值δ0(如10mm)时,调节悬浮电磁铁绕组11的励磁电流,使磁吸力F与机舱等悬浮物的重力mg达到平衡,机舱悬浮在空中,此处称为悬浮平衡点;然后定子三相绕组通电,则在定子2和悬浮电磁铁1之间的气隙中就会形成一个旋转磁场,旋转磁场和悬浮电磁铁1磁场耦合后,会形成气隙合成磁场。在合成磁场的磁力矩作用下,悬浮电磁铁1通过悬浮架4带动机舱旋转,实现偏航。在偏航旋转过程中,由于导向轴承5的作用,可以保证机舱不会发生偏移,使机舱沿着轨道正常旋转。
由此可见,上述磁悬浮偏航系统是在悬浮平衡点处进行偏航旋转,在此过程中,一方面要控制其转速,满足运动方程,同时更要控制其悬浮气隙保持恒定,使其在平衡点处旋转,实现偏航,因此实现平衡点处的稳定悬浮至关重要。
为此,本发明基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法,包括以下步骤:
步骤1,建立磁悬浮偏航系统的悬浮动态数学模型,建模过程如下:
如图2所示,悬浮电磁铁绕组11通电后将产生向上的轴向悬浮力F:
式中,i(t)为悬浮电磁铁绕组11的电流;δ(t)为悬浮电磁铁1和定子2之间的悬浮气隙;k1=μ0N2S/4,其中,μ0为真空磁导率,N为悬浮电磁铁绕组11的匝数,S为悬浮电磁铁铁心12的磁极表面有效面积。
此外,机舱等悬浮物在垂直方向上还受到向下的重力mg和外界扰动力fd(t),于是由牛顿第二定律,有垂直方向上的力学方程:
式中,m为悬浮物质量,g为重力加速度;为悬浮气隙δ(t)对时间t的二阶导数,即悬浮物的加速度;i(t)为悬浮电磁铁绕组11的电流。
同时悬浮电磁铁1的电压方程为:
式中,u(t)为悬浮电磁铁绕组11的电压;R为悬浮电磁铁绕组11的电阻;ψ(t)为气隙磁链;L为悬浮电磁铁绕组11的气隙电感,且有L=2k1/δ(t)。
综上,可得磁悬浮偏航系统的动态悬浮数学模型,即非线性模型:
平衡条件:F(i0,δ0)=k1(i0/δ0)2=mg
式中,i0、δ0分别为稳定悬浮时,即在平衡点处所述悬浮电磁铁绕组的电流和悬浮气隙。
步骤2,建立磁悬浮偏航系统的悬浮系统线性化模型,建模过程如下:
将式(1)在平衡点(i0,δ0)附近进行线性化处理,用泰勒级数公式展开,忽略高阶项,可得:
式中,ΔF=F(i0,δ0)-F;Δδ(t)=δ0-δ(t);Δi(t)=i0-i(t)。
类似地,将式(2)在平衡点(i0,δ0)附近进行线性化处理,用泰勒级数公式展开,忽略高阶项,可得:
式中,Δu(t)=u0-u(t),u0为平衡点处稳定悬浮时悬浮电磁铁绕组11的电压。
由于外界扰动fd(t)会使悬浮气隙δ(t)偏离δ0,所以可令其中,kf为常数,同时将式(4)、式(5)代入式(3),则可得悬浮系统的线性化模型:
式中,kc=2k1i0 2/δ0 3;ki=kcδ0/i0;kf为常数;L0=2k1/δ0;
步骤3,采用模型预测控制(MPC)策略,设计悬浮气隙控制器,实现外环悬浮气隙控制。具体方法如下:
将平衡点处的悬浮电磁铁绕组11的电流变化量Δi(t)、悬浮气隙变化量Δδ(t)和悬浮气隙变化率作为状态变量,悬浮电磁铁绕组11的电压变化量Δu(t)作为输入变量,即令x1=Δi(t),x2=Δδ(t),代入式(6),整理可得悬浮系统状态空间方程为:
式中,x=[x1,x2,x3]T,um=Δu(t),
y1为悬浮系统的输出量。
根据模型预测理论,在x可测时,式(7)的预测模型为:
x(k+1)=Ax(k)+bum(k) (8)
设从k时刻输入发生M步变化,而后保持不变,则由式(8)可以预测出um(k),um(k+1),um(k+2),…,um(k+M-1)作用在未来P(P>=M)个时刻的系统状态为:
上式可用向量形式表示为:
X(k)=Fx(k)+GU(k) (9)
式中,X(k)为状态矩阵,X(k)∈Rn实时可测;n为x(k)的维数;U(k)为系统输入矩阵;F、G分别是系统状态转移矩阵、输入矩阵,且有:
其中:P为预测时域长度;M为控制时域长度。
选取所述悬浮系统的优化性能指标为:
式中,Q为误差加权系数;R为控制加权系数。
在不考虑约束时,结合状态预测模型式(8),求出式(10)的最优解为:
U(k)=-(GTQG+R)-1GTQFx(k) (11)
由此可得输入变量um(k)控制律为:
式中,k2为反馈增益,有:式(12)即为悬浮气隙控制器。
步骤4,如图3所示,将步骤3悬浮气隙控制器的输出y1=x1=Δi(t)作为悬浮电磁铁绕组11的电流变化量参考值Δiref(k),可求得悬浮电磁铁绕组11的电流参考值iref(k)=i0+Δiref(k),将iref(k)与模型失配器式(13)的输出,即补偿电流y2相加,两者之和减去悬浮变流器6输出的电流测量值(即悬浮电磁铁绕组11的电流i(t)),得到误差e(t),经PID控制器,送入PWM模块,产生驱动信号,控制悬浮变流器6的输出电压和电流(即悬浮电磁铁绕组11的电压u(t)和电流i(t)),从而实现内环电流的跟踪控制。
上述基于平衡点的线性化模型式(6)因忽略了高阶项,使得它与实际非线性模型式(3)存在一定偏差,造成模型失配。为了消除模型失配及外界干扰给悬浮系统带来的影响,必须设计模型失配补偿器,将气隙、绕组电流的偏差送入模型失配补偿器,得到内环电流的补偿值,从而提高整个悬浮系统的控制精度和稳定性。模型失配补偿器按下述方式设计:
令模型失配补偿律为:
y2=C1eδ+C2ei (13)
式中,y2为模型失配补偿器输出的补偿电流;eδ=δ(t)-δl(t);ei=i(t)-il(t),其中,δ(t)、i(t)分别为悬浮气隙测量值(由气隙传感器3测得)、悬浮电磁铁绕组11的电流测量值;δl(t)、il(t)分别为线性化模型式(6)输出的悬浮气隙和悬浮电磁铁绕组11电流的计算值;C1是气隙偏差调节参数,C2是电流偏差调节参数;则式(13)即构成模型失配补偿器。
步骤5,如图3所示,将悬浮变流器6的输出电压u(t)送入悬浮线性化模型式(6),经式(6)计算得悬浮气隙计算值δl(t)和悬浮电磁铁绕组11的电流计算值il(t);将δl(t)减去气隙传感器3测得的悬浮气隙测量值δ(t),得到悬浮气隙偏差eδ,将il(t)减去悬浮电磁铁绕组11的电流测量值i(t),得到绕组电流偏差ei,将偏差eδ、ei送入模型失配补偿器式(13),得到补偿电流值y2,返回步骤4,实现平衡点附近稳定悬浮。
下面用一个优选实施例对本发明做进一步的说明。
由上可知,本发明基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法,实质上是一种加模型失配器的MPC-PID方法。为了验证本发明基于模型失配器的悬浮控制方法的有效性,对一款磁悬浮偏航系统实验样机分别采用本发明的加模型失配器的MPC-PID方法和不加模型失配器的控制方法(以下称为MPC-PID方法)进行对比仿真分析。
如图4所示,不加模型失配器的控制方法主要思路是:由悬浮动态数学模型式(3)得到悬浮线性化模型式(6),其输出Δu(t)作为MPC悬浮气隙控制器式(12)的输入,根据悬浮气隙控制器式(12)的输出,求得悬浮电磁铁绕组11的电流参考值iref(k)=i0+Δiref(k),将iref(k)减去悬浮变流器6输出的电流测量值(即悬浮电磁铁绕组11的电流i(t)),得到误差e(t),经PID控制器,送入PWM模块,产生驱动信号,控制悬浮变流器6的输出电压和电流(即悬浮电磁铁绕组11的电压u(t)和电流i(t)),从而实现内环电流的跟踪控制。
具体仿真参数如下表:
磁悬浮偏航系统的模型参数
名称 | 参数 |
悬浮物质量m(kg) | 250 |
悬浮电磁铁的绕组匝数N | 275 |
悬浮电磁铁的磁极面积S(mm2) | 259050 |
平衡点悬浮气隙δ0(mm) | 10 |
平衡点绕组电流io(A) | 10 |
真空磁导率μ0(N/A2) | 4π×10-7 |
仿真结果如图5、图6所示。
图5为悬浮气隙对比仿真分析曲线图。由图5可知,在启动阶段,两种控制方法都能够快速调节,进入稳态之后,本发明由于加入了模型失配补偿器,悬浮气隙稳定在10mm,而未考虑补偿器的MPC-PID方法,其悬浮气隙稳定在9.4mm。由此可见,加入模型失配补偿器提高了磁悬浮偏航系统的控制精度。
图6为悬浮电磁铁绕组电流对比仿真分析曲线图。由图5可知,本发明输出电流的响应速度更快,同时输出电流稳定在9.92A,接近目标值10A;未加入模型失配补偿器的MPC-PID方法电流稳定在9.84A,误差较大。
从图5、图6可以看出,本发明基于模型失配补偿器的悬浮控制方法能够提高整个系统的快速性、精确性、稳定性。
Claims (1)
1.一种基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法,所述风电磁悬浮偏航系统包括悬浮系统,所述悬浮系统包括悬浮电磁铁、定子、气隙传感器、悬浮架和悬浮变流器,所述悬浮电磁铁包括铁心和绕组,所述绕组为直流励磁绕组,所述绕组与所述悬浮变流器连接;其特征在于,包括以下步骤:
步骤1,建立磁悬浮偏航系统的悬浮动态数学模型,建模过程如下:
所述悬浮电磁铁绕组通电后将产生向上的轴向悬浮力F:
式中,i(t)为所述悬浮电磁铁绕组的电流;δ(t)为所述悬浮电磁铁和所述定子之间的悬浮气隙;k1=μ0N2S/4,其中,μ0为真空磁导率,N为所述悬浮电磁铁绕组的匝数,S为所述悬浮电磁铁铁心的磁极表面有效面积;
根据牛顿第二定律,所述悬浮系统在垂直方向上的力学方程为:
式中,m为悬浮物质量,g为重力加速度;fd(t)为外界扰动力,为悬浮气隙δ(t)对时间t的二阶导数;
所述悬浮电磁铁的电压方程为:
式中,u(t)为所述悬浮电磁铁绕组的电压;R为所述悬浮电磁铁绕组的电阻;ψ(t)为气隙磁链;L为所述悬浮电磁铁绕组的气隙电感,且有L=2k1/δ(t);
综上,可得磁悬浮偏航系统的动态悬浮数学模型:
平衡条件:F(i0,δ0)=k1(i0/δ0)2=mg
式中,i0、δ0分别为稳定悬浮时,即在平衡点处所述悬浮电磁铁绕组的电流和悬浮气隙;
步骤2,建立磁悬浮偏航系统的悬浮系统线性化模型,建模过程如下:
将式(1)在平衡点(i0,δ0)附近进行线性化处理,用泰勒级数公式展开,忽略高阶项,可得:
式中,ΔF=F(i0,δ0)-F;Δδ(t)=δ0-δ(t);Δi(t)=i0-i(t);
类似地,将式(2)在平衡点(i0,δ0)附近进行线性化处理,用泰勒级数公式展开,忽略高阶项,可得:
式中,Δu(t)=u0-u(t),u0为平衡点处稳定悬浮时悬浮电磁铁的电压;
将式(4)、式(5)代入式(3),同时令则可得所述悬浮系统的线性化模型:
式中,kc=2k1i0 2/δ0 3;ki=kcδ0/i0;kf为常数;L0=2k1/δ0;
步骤3,采用模型预测控制策略,设计悬浮气隙控制器,实现外环悬浮气隙控制;具体方法如下:
A)令x1=Δi(t),x2=Δδ(t),代入式(6),整理可得悬浮系统状态空间方程为:
式中,x=[x1,x2,x3]T为状态变量,um=Δu(t)为输入变量,
y1为所述悬浮系统的输出量;
B)根据模型预测理论,在x可测时,式(7)的预测模型为:
x(k+1)=Ax(k)+bum(k) (8)
C)设从k时刻输入发生M步变化,而后保持不变,则由式(8)可以预测出um(k),um(k+1),um(k+2),…,um(k+M-1)作用在未来P(P>=M)个时刻的系统状态,即:
X(k)=Fx(k)+GU(k) (9)
式中,X(k)为状态矩阵,X(k)∈Rn实时可测;n为x(k)的维数;U(k)为系统输入矩阵;F、G分别是系统状态转移矩阵、输入矩阵,且有:
其中:P为预测时域长度;M为控制时域长度;
D)选取所述悬浮系统的优化性能指标为:
式中,Q为误差加权系数;R为控制加权系数;
在不考虑约束时,结合状态预测模型式(8),求出式(10)的最优解为:
U(k)=-(GTQG+R)-1GTQFx(k) (11)
E)由此可得输入变量控制律为:
式中,k2为反馈增益,有:
式(12)构成了悬浮气隙控制器;
步骤4,将步骤3所述悬浮气隙控制器的输出y1=x1=Δi(t)作为所述悬浮电磁铁绕组的电流变化量参考值Δiref(k),据此求得所述悬浮电磁铁绕组的电流参考值iref(k)=i0+Δiref(k),将此值iref(k)与模型失配器的输出相加,两者之和与所述悬浮变流器输出的电流测量值之差,经PID控制器,送入PWM模块,产生驱动信号,控制所述悬浮变流器的输出电压和电流,实现内环电流跟踪控制;
所述模型失配补偿器按下述方式设计:
令模型失配补偿律为:
y2=C1eδ+C2ei (13)
式中,y2为模型失配补偿器输出的补偿电流;eδ=δ(t)-δl(t);ei=i(t)-il(t),其中,δ(t)、i(t)分别为悬浮气隙测量值、所述悬浮电磁铁绕组电流测量值;δl(t)、il(t)分别为线性化模型式(6)输出的悬浮气隙和所述悬浮电磁铁绕组励磁电流;C1是气隙偏差调节参数;C2是电流偏差调节参数;则式(13)构成了模型失配补偿器;
步骤5,将所述悬浮变流器的输出电压送入线性化模型式(6),线性化模型式(6)输出悬浮气隙δl(t)和所述悬浮电磁铁绕组的电流il(t);所述悬浮气隙δl(t)减去悬浮气隙测量值δ(t),得到悬浮气隙偏差eδ,所述悬浮电磁铁绕组的电流il(t)减去所述悬浮电磁铁绕组的电流测量值i(t),得到绕组电流偏差ei;将所述偏差eδ、ei送入模型失配补偿器式(13),得到补偿电流y2,返回步骤4,实现平衡点附近稳定悬浮。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810416971.2A CN108488036B (zh) | 2018-05-04 | 2018-05-04 | 基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810416971.2A CN108488036B (zh) | 2018-05-04 | 2018-05-04 | 基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108488036A true CN108488036A (zh) | 2018-09-04 |
CN108488036B CN108488036B (zh) | 2019-10-25 |
Family
ID=63352948
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810416971.2A Active CN108488036B (zh) | 2018-05-04 | 2018-05-04 | 基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108488036B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110071677A (zh) * | 2019-05-30 | 2019-07-30 | 中国科学院电工研究所 | 高速磁悬浮列车长定子直线同步电机牵引控制方法 |
CN110794680A (zh) * | 2019-11-14 | 2020-02-14 | 浙江工业大学 | 一种基于扩张状态观测器的磁悬浮球系统预测跟踪控制方法 |
CN111682815A (zh) * | 2020-06-17 | 2020-09-18 | 曲阜师范大学 | 一种含高频干扰重构的风力机舱偏航控制方法 |
CN113775474A (zh) * | 2021-08-21 | 2021-12-10 | 曲阜师范大学 | 基于自适应神经网络有限时间控制的垂直轴风电机组悬浮控制方法 |
CN113847202A (zh) * | 2021-09-16 | 2021-12-28 | 国电投河南新能源有限公司 | 一种风电机组变桨稳转速控制系统及控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1669294A (zh) * | 2002-07-16 | 2005-09-14 | 皇家飞利浦电子股份有限公司 | 具有模型失配补偿的回声消除器 |
CN202157900U (zh) * | 2011-06-07 | 2012-03-07 | 华锐风电科技(集团)股份有限公司 | 一种基于内模pid的风电机组变桨控制系统 |
KR20140018497A (ko) * | 2012-08-01 | 2014-02-13 | 한국전력공사 | 단기 풍속 예측을 통한 풍력 발전량 예측 방법과, 그 기능을 이용한 배전선로 전압 예측 방법 |
US9219629B1 (en) * | 2014-08-27 | 2015-12-22 | Zenith Electronics L.L.C. | Adaptive low complexity channel estimation for mobile OFDM systems |
CN106988961A (zh) * | 2017-04-17 | 2017-07-28 | 曲阜师范大学 | 混合悬浮气隙调整型风机偏航系统 |
WO2018037191A1 (fr) * | 2016-08-25 | 2018-03-01 | Institut Polytechnique De Grenoble | Procede et dispositif d'estimation de force |
-
2018
- 2018-05-04 CN CN201810416971.2A patent/CN108488036B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1669294A (zh) * | 2002-07-16 | 2005-09-14 | 皇家飞利浦电子股份有限公司 | 具有模型失配补偿的回声消除器 |
CN202157900U (zh) * | 2011-06-07 | 2012-03-07 | 华锐风电科技(集团)股份有限公司 | 一种基于内模pid的风电机组变桨控制系统 |
KR20140018497A (ko) * | 2012-08-01 | 2014-02-13 | 한국전력공사 | 단기 풍속 예측을 통한 풍력 발전량 예측 방법과, 그 기능을 이용한 배전선로 전압 예측 방법 |
US9219629B1 (en) * | 2014-08-27 | 2015-12-22 | Zenith Electronics L.L.C. | Adaptive low complexity channel estimation for mobile OFDM systems |
WO2018037191A1 (fr) * | 2016-08-25 | 2018-03-01 | Institut Polytechnique De Grenoble | Procede et dispositif d'estimation de force |
CN106988961A (zh) * | 2017-04-17 | 2017-07-28 | 曲阜师范大学 | 混合悬浮气隙调整型风机偏航系统 |
Non-Patent Citations (2)
Title |
---|
李洋: "基于PID控制器的磁悬浮偏航系统悬浮控制", 《电子技术》 * |
李洋: "风电磁悬浮偏航系统悬浮及切换控制研究", 《工程科技II辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110071677A (zh) * | 2019-05-30 | 2019-07-30 | 中国科学院电工研究所 | 高速磁悬浮列车长定子直线同步电机牵引控制方法 |
CN110794680A (zh) * | 2019-11-14 | 2020-02-14 | 浙江工业大学 | 一种基于扩张状态观测器的磁悬浮球系统预测跟踪控制方法 |
CN110794680B (zh) * | 2019-11-14 | 2022-04-05 | 浙江工业大学 | 一种基于扩张状态观测器的磁悬浮球系统预测跟踪控制方法 |
CN111682815A (zh) * | 2020-06-17 | 2020-09-18 | 曲阜师范大学 | 一种含高频干扰重构的风力机舱偏航控制方法 |
CN111682815B (zh) * | 2020-06-17 | 2023-04-14 | 曲阜师范大学 | 一种含高频干扰重构的风力机舱偏航控制方法 |
CN113775474A (zh) * | 2021-08-21 | 2021-12-10 | 曲阜师范大学 | 基于自适应神经网络有限时间控制的垂直轴风电机组悬浮控制方法 |
CN113775474B (zh) * | 2021-08-21 | 2023-09-29 | 曲阜师范大学 | 基于自适应神经网络有限时间控制的垂直轴风电机组悬浮控制方法 |
CN113847202A (zh) * | 2021-09-16 | 2021-12-28 | 国电投河南新能源有限公司 | 一种风电机组变桨稳转速控制系统及控制方法 |
CN113847202B (zh) * | 2021-09-16 | 2023-06-06 | 国电投河南新能源有限公司 | 一种风电机组变桨稳转速控制系统及控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108488036B (zh) | 2019-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108488036B (zh) | 基于模型失配补偿器的风电磁悬浮偏航系统悬浮控制方法 | |
CN108183650A (zh) | 一种基于模型预测控制的风电磁悬浮偏航电机控制方法 | |
Arturo Soriano et al. | Modeling and control of wind turbine | |
Marino et al. | Induction motor control design | |
CN110401378B (zh) | 基于神经网络-模型预测控制的磁悬浮偏航电机控制方法 | |
CN106786677B (zh) | 一种互联电力系统分布式动态矩阵频率控制方法 | |
CN110219774B (zh) | 一种水平轴风力偏航系统机舱悬浮控制方法 | |
CN105048913A (zh) | 基于电流补偿的无轴承异步电机不平衡振动控制系统 | |
CN104660136A (zh) | 无轴承异步电机的不平衡振动控制系统 | |
CN108019316A (zh) | 主被动协同调控的磁悬浮风力偏航系统 | |
CN110345013A (zh) | 基于神经网络模型预测控制的磁悬浮垂直轴风电机组控制方法 | |
CN111835261B (zh) | 基于自适应神经网络的磁悬浮垂直轴风电机组悬浮控制方法 | |
Li et al. | Modeling of maglev yaw system of wind turbines and its robust trajectory tracking control in the levitating and landing process based on NDOB | |
CN110805523B (zh) | 一种基于神经网络的风电磁悬浮偏航系统悬浮控制方法 | |
Sun et al. | Neuron PID control for a BPMSM based on RBF neural network on‐line identification | |
Arribas et al. | Low voltage ride-through in DFIG wind generators by controlling the rotor current without crowbars | |
CN111082443A (zh) | 一种并网调频型飞轮储能系统放电过程中电机侧控制方法 | |
CN107165779A (zh) | 一种垂直轴悬浮永磁风力发电机及其控制方法 | |
Vu et al. | A novel nonlinear observer‐based LQ control system design for wind energy conversion systems with single measurement | |
CN102013870B (zh) | 五自由度无轴承同步磁阻电机逆系统解耦控制器 | |
CN111622899B (zh) | 一种水平轴风力机舱同步悬浮控制方法 | |
CN111173681B (zh) | 基于滑模神经网络的磁悬浮垂直轴风电机组悬浮控制方法 | |
CN113775474B (zh) | 基于自适应神经网络有限时间控制的垂直轴风电机组悬浮控制方法 | |
Wang et al. | Research on suspension control strategy based on finite control set model predictive control with state feedback control‐PID for maglev yaw system of wind turbine | |
Wang et al. | Research on MPC-based suspension control of wind maglev yaw system |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240111 Address after: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province Patentee after: Dragon totem Technology (Hefei) Co.,Ltd. Address before: 273165 Jingxuan West Road, Qufu City, Jining, Shandong Province, No. 57 Patentee before: QUFU NORMAL University |