CN109976403B - 多陷波器并联的磁悬浮转子振动谐波抑制方法及系统 - Google Patents

多陷波器并联的磁悬浮转子振动谐波抑制方法及系统 Download PDF

Info

Publication number
CN109976403B
CN109976403B CN201910254380.4A CN201910254380A CN109976403B CN 109976403 B CN109976403 B CN 109976403B CN 201910254380 A CN201910254380 A CN 201910254380A CN 109976403 B CN109976403 B CN 109976403B
Authority
CN
China
Prior art keywords
parallel
phase shift
rotor
harmonic
trap
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
CN201910254380.4A
Other languages
English (en)
Other versions
CN109976403A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201910254380.4A priority Critical patent/CN109976403B/zh
Publication of CN109976403A publication Critical patent/CN109976403A/zh
Application granted granted Critical
Publication of CN109976403B publication Critical patent/CN109976403B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D19/00Control of mechanical oscillations, e.g. of amplitude, of frequency, of phase
    • G05D19/02Control of mechanical oscillations, e.g. of amplitude, of frequency, of phase characterised by the use of electric means

Abstract

本发明公开了一种多陷波器并联的磁悬浮转子振动谐波抑制方法及系统,包括依次连接在一起的并联相移陷波器、磁悬浮控制器、磁轴承转子,其中,并联相移陷波器的信号提取点放在磁悬浮控制器的输入端,并联相移陷波器的信号插入点放在位移误差信号处,发明针对磁悬浮转子,考虑转子质量不平衡和传感器振动的影响,建立包含谐波扰动的磁悬浮转子动力学模型,通过并联不同相移的多个陷波器有效抑制变转速频率下谐波电流的不同频率扰动分量,进而有效抑制谐波振动,同时保证系统的闭环稳定性,并对磁悬浮‑转子系统变速情况下的谐波抑制问题具有参考意义。

Description

多陷波器并联的磁悬浮转子振动谐波抑制方法及系统
技术领域
本发明涉及一种不同相位的多陷波器并联的磁悬浮转子控制系统,从而获得全转速范围内优越的谐波电流抑制效果,进而有效抑制谐波振动。本发明属于磁悬浮-转子系统谐波振动抑制领域。
背景技术
磁悬浮飞轮转子由于具有非接触、低功耗、长寿命等优点,成为卫星稳定性研究的热点。为了获得超高精度的卫星姿态控制性能,非常有必要尽可能地减少磁悬浮转子不需要的振动力。电磁轴承-转子系统中不需要的振动力主要是由谐波电流引起的,这是磁悬浮转子系统中普遍存在的问题,它对系统稳定性和性能下降构成了潜在的威胁。谐波电流包含同步电流和多频电流。质量不平衡是主要的振源,并将在控制系统中引入同步电流。考虑到敏感表面的圆度误差和非理想的电磁特性,传感器振动会产生同步的多频电流。
为了控制振动,谐波电流的抑制的方法也在被不断提出。然而方法大多针对固定频率的多频振动控制,由于磁悬浮转子的转速变化范围很大,不能满足实际运行的需要。在电磁轴承振动主动控制的应用中也发现了一些现代控制算法,但是大多数现代控制方法都是基于模型的。数学模型的准确性是保证系统性能良好的关键。然而,获取被控对象模型并不是一件容易的事情,在实践中的应用受到了一定的限制。
在闭环系统中加入陷波器的方法对于振动抑制问题的解决有很大意义,例如之前提出了一种相移陷波器来并联基线磁悬浮控制器,该陷波器在抑制同步电流的同时,保证了系统在整个速度范围内的稳定性。通常都采用多陷波器抑制谐波电流,但是,多陷波器具有相同的结构和参数,在整个转速范围内不能保证闭环稳定性。由于不同的谐波在较宽的速度范围内交叉,当陷波器插入到基线系统中时,不可避免地会影响闭环稳定性。在全转速范围内有效抑制谐波振动同时保证系统的全局稳定性,是磁轴承转子振动控制系统中的很有挑战性的一个问题。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种多陷波器并联的磁悬浮转子振动谐波抑制方法及系统,通过并联不同相移的多个陷波器有效抑制变转速频率下谐波电流的不同频率扰动分量,进而获得高效的抑制谐波振动性能,同时保证系统的闭环稳定性。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种多陷波器并联的磁悬浮转子振动谐波抑制方法,通过引入不同的移相的多个陷波器并联来抑制磁悬浮-转子系统在时变转速频率下谐波电流的不同频率扰动分量,从而在全转速范围内有效抑制谐波振动。
具体包括以下骤:
步骤1),采集转子转速、谐波的直流分量,考虑转子质量不平衡和传感器振动,建立包含谐波扰动的磁悬浮转子动力学模型;
步骤2),搭建带有以并行模式连接的多个相移陷波器的闭环控制系统,明确相移陷波器的插入方式;
步骤21),设计用于消除时变转速频率下控制电流中同频扰动分量的相移陷波器。
传统的陷波器通常具有如下的传递函数形式:
Figure GDA0003068569840000021
其中,N(s)为传统陷波器的传递函数,s为拉普拉斯算子,Ω表示转子转速,ε是相移陷波器的增益系数,它影响着算法的收敛速度。
步骤22),在传统陷波器的基础上引入调节相位φ,构成相移陷波器,表示如下:
Figure GDA0003068569840000022
多相移陷波器并联表示为:
Figure GDA0003068569840000023
等价形式为:
Figure GDA0003068569840000024
其中,Nms(s)表示多相移陷波器并联的传递函数,n表示谐波次数,调节相位φ在时变转速频率下进行自我调节以保证整个闭环系统的稳定性,从而达到消除时变转速频率下控制电流中同频扰动分量的效果;
步骤3),根据步骤2得到的相移陷波器设计具有相同相移的多个并联相移陷波器;
步骤4),根据步骤2)搭建的闭环控制系统,对谐波电流抑制闭环系统进行稳定性分析,确定不同相位并联的多陷波器下整个系统的稳定条件;
步骤5),根据步骤2)、步骤3)和步骤4),利用有多个不同相位并联的相移滤波器对谐波电流进行有效抑制:
并行模式下不同相位并联的相移滤波器的抑制性能显示如下:
Figure GDA0003068569840000025
Figure GDA0003068569840000031
其中,Np(s)表示灵敏度函数,s表示拉普拉斯算子,ε表示相移陷波器的增益系数,Nmd(s)表示并联相移陷波器的传递函数,Ω表示转子转速,n表示谐波次数,j表示虚部单位,
Figure GDA0003068569840000032
具有并行模式不同相位并联的相移滤波器的闭环系统中从e(s)到xd(s)的传递函数表示为:
Figure GDA0003068569840000033
将整个系统的特征方程表示为:
Figure GDA0003068569840000034
其中,GNd(s)表示从e(s)到xd(s)的传递函数,Gc(s)表示磁悬浮控制器,P(s)表示磁悬浮转子,q(s)表示整个系统的特征方程,
在ε=0的情况下,根是s=±jΩ,±j2Ω,...,±jnΩ,在虚轴附近的根应该配置在左半复平面,对ε求导,满足下列式中所示的所有条件,以使整个系统是稳定的:
Figure GDA0003068569840000035
根据上式求导可推导出在全转速范围内不同频率下系统稳定的条件为:
Figure GDA0003068569840000036
根据基本磁悬浮控制系统的灵敏度函数的相角变化,对不同频率的转速进行分段处理,从而在各转速段选取满足上述稳定性条件的相移陷波器的调节相位值。
优选的:所述步骤1)中的建立包含谐波扰动的磁悬浮转子动力学模型过程如下:
位移传感器检测到的真实信号表示为:
xs(t)=x(t)+xd(t)
ys(t)=y(t)+yd(t)
其中,xs(t)和ys(t)表示x和y方向位移传感器检测到的真实信号,x(t)和y(t)是x和y方向几何中心的位移信号,xd(t)和yd(t)是x和y方向谐波扰动,用傅里叶分解的方法描述为:
Figure GDA0003068569840000041
Figure GDA0003068569840000042
其中,xd0和yd0是x和y方向谐波的直流分量,n是谐波次数,axi,bxi,ayi和byi是x和y方向第i倍频谐波的傅里叶系数,Ω是转子转速,t表示时间变量;
根据牛顿第二定律,径向磁力轴承-转子系统的磁悬浮转子动力学模型为:
Figure GDA0003068569840000043
Figure GDA0003068569840000044
其中,
Figure GDA0003068569840000045
Figure GDA0003068569840000046
表示磁悬浮转子x和y方向的所受的力,kx和ky为x和y方向的位移刚度参数,kix和kiy为x和y方向的当前刚度参数,m为转子质量,η和
Figure GDA0003068569840000047
为不平衡质量的偏心和相位,ix(t)和iy(t)为x和y方向的磁轴承电流,t为时间变量。
优选的:步骤3)中通过具有相同相移的多个并联相移陷波器稳定闭环系统的方法:
将闭环系统特征方程写为:
1+Gc(s)P(s)=0
通过步骤二明确多个并联具有相同相移的陷波器在闭环系统中的表达形式,因此进一步推导出插入陷波器后系统的特征方程为:
Figure GDA0003068569840000048
其中,s表示拉普拉斯算子,n表示谐波次数,j表示虚部单位。
对于移相来说稳定条件为:
Figure GDA0003068569840000049
也就是说,移相同时满足以下条件:
-90°≤arg[S(jΩ)]+φ≤90°
优选的:步骤4)中不同相位并联的多陷波器,其传递函数表达如下:
Figure GDA0003068569840000051
它的等价形式为:
Figure GDA0003068569840000052
其中,Nmd(s)表示不同相位并联的多陷波器的传递函数,n表示第几次谐波,φi是第i倍频的相位。
一种多陷波器并联的磁悬浮转子振动谐波抑制系统,包括依次连接在一起的并联相移陷波器Nmd(s)、磁悬浮控制器Gc(s)、磁轴承转子P(s),其中,并联相移陷波器Nmd(s)的信号提取点放在磁悬浮控制器Gc(s)的输入端,并联相移陷波器Nmd(s)的信号插入点放在位移误差信号处,参考位移信号x*(t)和位移传感器检测到的位移信号xs(t)得到位移误差信号,参考位移信号x*(t)、位移误差信号、以及并联相移陷波器Nmd(s)的信号进入到磁悬浮控制器Gc(s),磁悬浮控制器Gc(s)产生控制信号对磁轴承转子P(s)、扰动信号xd(s)进行控制,位移传感器检测到转子的位移信号xs(t);
并联相移陷波器Nmd(s)中:
不同相位并联的多陷波器传递函数如下:
Figure GDA0003068569840000053
它的等价形式为:
Figure GDA0003068569840000054
并行模式下不同陷波器的抑制性能显示如下:
Figure GDA0003068569840000055
其中,
Figure GDA0003068569840000056
具有并行模式不同相位的多陷波器的闭环系统中从e(s)到xd(s)的传递函数表示为:
Figure GDA0003068569840000061
将整个系统的特征方程表示为:
Figure GDA0003068569840000062
在ε=0的情况下,根是s=±jΩ,±j2Ω,...,±jnΩ,在虚轴附近的根应该配置在左半复平面上,对ε求导,满足下列式中所示的所有条件,以使整个系统是稳定的:
Figure GDA0003068569840000063
根据上式求导可推导出在全转速范围内不同频率下系统稳定的条件为:
Figure GDA0003068569840000064
其中:Np(s)表示灵敏度函数,s表示拉普拉斯算子,ε表示相移陷波器的增益系数,Nmd(s)表示并联相移陷波器的传递函数,Ω表示转子转速,n表示谐波次数,j表示虚部单位。
优选的:位移传感器检测到的真实信号表示为:
xs(t)=x(t)+xd(t)
ys(t)=y(t)+yd(t)
其中,xs(t)和ys(t)表示x和y方向位移传感器检测到的真实信号,x(t)和y(t)是x和y方向几何中心的位移信号,xd(t)和yd(t)是x和y方向谐波扰动,用傅里叶分解的方法描述为:
Figure GDA0003068569840000065
Figure GDA0003068569840000066
其中,xd0和yd0是x和y方向谐波的直流分量,n是谐波次数,axi,bxi,ayi和byi是x和y方向第i倍频谐波的傅里叶系数,Ω是转子转速,t表示时间变量。
优选的:磁轴承转子P(s)的磁悬浮转子动力学模型为:
Figure GDA0003068569840000071
Figure GDA0003068569840000072
其中,
Figure GDA0003068569840000073
Figure GDA0003068569840000074
表示磁悬浮转子x和y方向的所受的力,kx和ky为x和y方向的位移刚度参数,
Figure GDA0003068569840000075
Figure GDA0003068569840000076
为x和y方向的当前刚度参数,m为转子质量,η和
Figure GDA0003068569840000077
为不平衡质量的偏心和相位,ix(t)和iy(t)为x和y方向的磁轴承电流,t为时间变量。
本发明相比现有技术,具有以下有益效果:
1.本发明提出一种通过并联不同相位多个陷波器的磁悬浮转子振动谐波抑制方法,其主要控制思想是对于不同频率的扰动分量设计不同的相位补偿,在达到谐波电流有效抑制的同时,对于不同的相位能够确保整个系统的稳定性,从而实现全转速范围内的谐波电流抑制,来有效抑制谐波振动。
2.本发明提出的谐波振动抑制主要有两个优势。首先,通过在磁悬浮控制系统中插入多个并联的不同相位陷波器,以获得全转速范围内优越的谐波振动抑制效果。其次,在不同的频率分段采用不同的相角补偿以满足系统稳定条件,从而保证整个系统的全局稳定性。
3.为了保证系统的稳定性,需要多个不同的换挡阶段,以满足文权利要求5中所述的相应条件。这样,即使灵敏度函数过了一个严重的斜率,并且需要消除高倍频分量,每个补偿相角也可以被调整到一个适当的值,并保持系统的稳定。
附图说明
图1为基本磁悬浮控制系统框图;
图2为带有相移陷波器的磁悬浮控制系统框图;
图3为带有不同相位陷波器并联的磁悬浮控制系统框图;
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
一种多陷波器并联的磁悬浮转子振动谐波抑制方法,研究对象磁悬浮转子是由两个径向磁力轴承和一个3自由度轴向磁力轴承悬浮在一起的。前者用来控制转子的两个径向平移运动。三自由度磁力轴承用于两种径向旋转运动。由于平动运动和转动运动是解耦的,因此主要研究径向平动运动的多重振动力抑制问题。其主要控制思想包含两个方面,一方面利用不同相位并联的陷波器在全转速范围内抑制同频电流扰动;另一方面在不同的频率分段采用不同的相角补偿保证系统的全局稳定性。通过引入不同的移相的多个陷波器并联来抑制磁悬浮-转子系统在时变转速频率下谐波电流的不同频率扰动分量,从而在全转速范围内有效抑制谐波振动。针对磁悬浮转子,考虑转子质量不平衡和传感器振动的影响,建立包含谐波扰动的磁悬浮转子动力学模型,通过并联不同相移的多个陷波器有效抑制变转速频率下谐波电流的不同频率扰动分量,进而有效抑制谐波振动,同时保证系统的闭环稳定性,并对磁悬浮-转子系统变速情况下的谐波抑制问题具有参考意义。
步骤一、考虑转子质量不平衡和传感器振动的影响,建立包含谐波扰动的磁悬浮转子动力学模型。
考虑到质量不平衡和传感器振动对转子动力学的影响。
位移传感器检测到的真实信号xs(t)和ys(t)可以表示为:
xs(t)=x(t)+xd(t)
ys(t)=y(t)+yd(t)
其中x(t)和y(t)是x和y方向几何中心的位移信号,xd(t)和yd(t)是谐波扰动,可以用傅里叶分解的方法更详细地描述为:
Figure GDA0003068569840000081
Figure GDA0003068569840000082
其中Ω是转子的转速,xd0和yd0是谐波的直流分量,n是谐波次数,axi,bxi,ayi和byi是第i倍频谐波的傅里叶系数。从式中可以看出,位移信号含有相当大的谐波。
根据牛顿第二定律,径向磁力轴承-转子系统可以建模为:
Figure GDA0003068569840000083
Figure GDA0003068569840000084
其中kx和ky为x和y方向的位移刚度参数,kix和kiy为当前刚度参数,m为转子质量,η和
Figure GDA0003068569840000085
为不平衡质量的偏心和相位。
结合图1,质量不平衡和传感器振动会在控制电流中引入一些谐波分量。因此,为了减小电磁轴承-转子系统中不需要的振动力,必须考虑谐波电流的抑制。
步骤二、搭建并联的多个相移陷波器的闭环控制系统,明确相移陷波器的插入方式。
首先,设计用于消除时变转速频率下控制电流中同频扰动分量的相移陷波器。
传统的陷波器通常具有如下的传递函数形式:
Figure GDA0003068569840000091
其中,s为拉普拉斯算子,Ω表示磁悬浮飞轮转子的转速,ε是相移陷波器的增益系数,它影响着算法的收敛速度。
然而,考虑磁悬浮转子需要在较宽的转速范围内工作,不平衡振动抑制算法的稳定性范围很宽,并且与转子转速自适应变化。而传统陷波器将会在虚轴部分引入极点,因此传统陷波器仅仅可以在一段转速范围内消除某个频率的扰动分量。为了克服变转速情况下的潜在失稳特性,让陷波器可以在全转速范围内维持整个系统的稳定性,本发明在传统陷波器的基础上引入调节相位φ,构成相移陷波器,表示如下:
Figure GDA0003068569840000092
多陷波器并联可表示为:
Figure GDA0003068569840000093
等价形式为:
Figure GDA0003068569840000094
其中,调节相位φ在时变转速频率下进行自我调节以保证整个闭环系统的稳定性,从而达到消除时变转速频率下控制电流中同频扰动分量的效果。
其次,在基本磁悬浮控制系统中插入相移陷波器,得到带有相移陷波器的控制系统框图。
实施例中,基本磁悬浮控制系统主要包括磁悬浮控制器,功率放大器和磁轴承转子。
步骤三、明确插入多陷波器后系统的稳定条件。
通过图2磁轴承控制框图,可以将闭环系统特征方程写为:
1+Gc(s)P(s)=0
通过步骤二可以明确多个并联具有相同相移的陷波器在闭环系统中的表达形式,因此进一步推导出插入陷波器后系统的特征方程为:
Figure GDA0003068569840000095
对于一些移相来说稳定条件为:
Figure GDA0003068569840000101
也就是说,移相同时满足以下条件:
-90°≤arg[S(jΩ)]+φ≤90°
基本磁悬浮控制系统如图1所示,x*(t)是参考位移信号,xs(t)为位移传感器检测到的位移信号,xd(s)为扰动信号,Gc(s)为磁悬浮控制器,相移陷波器Nms(s)的信号提取点放在磁悬浮控制器Gc(s)的输入端,相移陷波器Nms(s)的信号插入点放在位移误差信号处,得到带有多个相移陷波器并联的悬浮控制系统框图如图2所示。也就是说,信号提取点位于信号插入点之后,相移陷波器与磁悬浮控制器串联,在这种信号处理顺序下,控制电流中的谐波扰动能够被有效消除。
步骤四、设计一种基于不同相位并联的多陷波器的闭环控制系统。
从步骤三闭环系统稳定条件可以看出,相同相位的多个陷波器并联并不能实现全转速范围内的稳定从而不能在全转速范围进行有效的谐波电流抑制,因此设计不同相位并联的多陷波器,其传递函数表达如下:
Figure GDA0003068569840000102
可以推导出它的等价形式为:
Figure GDA0003068569840000103
其中n表示第几次谐波,φi是第i倍频的相位,可以看出在不同频率下采取了不同的相位补偿。
步骤五、对插入不同相位陷波器并联的控制系统进行抑制性能分析及稳定性分析,确定不同相位的稳定性条件。
并行模式下不同陷波器的抑制性能显示如下:
Figure GDA0003068569840000104
其中,
Figure GDA0003068569840000105
具有并行模式不同相位的多陷波器的闭环系统中从e(s)到xd(s)的传递函数可以表示为:
Figure GDA0003068569840000111
将整个系统的特征方程表示为:
Figure GDA0003068569840000112
从上述特征方程中可以看出,在ε=0的情况下,根可以是s=±jΩ,±j2Ω,...,±jnΩ。为了确保所有的闭环极点都有负实部,从而使整个系统稳定,在虚轴附近的根应该配置在左半复平面上。因此,让我们对ε求导,满足下列式中所示的所有条件,以使整个系统是稳定的:
Figure GDA0003068569840000113
根据上式求导可推导出在全转速范围内不同频率下系统稳定的条件为:
Figure GDA0003068569840000114
根据基本磁悬浮控制系统的灵敏度函数的相角变化,对不同频率的转速进行分段处理,从而在各转速段选取满足上述稳定性条件的相移陷波器的调节相位值。
一种多陷波器并联的磁悬浮转子振动谐波抑制系统,如图3所示,包括依次连接在一起的并联相移陷波器Nmd(s)、磁悬浮控制器Gc(s)、磁轴承转子P(s),其中,并联相移陷波器Nmd(s)的信号提取点放在磁悬浮控制器Gc(s)的输入端,并联相移陷波器Nmd(s)的信号插入点放在位移误差信号处,参考位移信号x*(t)和位移传感器检测到的位移信号xs(t)得到位移误差信号,参考位移信号x*(t)、位移误差信号、以及并联相移陷波器Nmd(s)的信号进入到磁悬浮控制器Gc(s),磁悬浮控制器Gc(s)产生控制信号对磁轴承转子P(s)、扰动信号xd(s)进行控制,位移传感器检测到转子的位移信号xs(t)。其中,x*(t)是参考位移信号,xs(t)为位移传感器检测到的位移信号,xd(s)为扰动信号,Gc(s)为磁悬浮控制器,相移陷波器Nmd(s)的信号提取点放在磁悬浮控制器Gc(s)的输入端,相移陷波器Nmd(s)的信号插入点放在位移误差信号处。通过图3所示的磁悬浮转子谐波振动控制系统,能够以较高精度获得全转速范围内谐波振动的抑制效果同时保证系统的稳定性。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (3)

1.一种多陷波器并联的磁悬浮转子振动谐波抑制方法,其特征在于,包括以下步骤:
步骤1),采集转子转速、谐波的直流分量,考虑转子质量不平衡和传感器振动,建立包含谐波扰动的磁悬浮转子动力学模型;
建立包含谐波扰动的磁悬浮转子动力学模型过程如下:
位移传感器检测到的真实信号表示为:
xs(t)=x(t)+xd(t)
ys(t)=y(t)+yd(t)
其中,xs(t)表示x方向位移传感器检测到的真实信号,ys(t)表示y方向位移传感器检测到的真实信号,x(t)是x方向几何中心的位移信号,y(t)是y方向几何中心的位移信号,xd(t)是x方向谐波扰动,yd(t)是y方向谐波扰动,用傅里叶分解的方法描述为:
Figure FDA0003068569830000011
Figure FDA0003068569830000012
其中,xd0是x方向谐波的直流分量,yd0是y方向谐波的直流分量,n是谐波次数,axi是x方向第i倍频谐波的傅里叶正弦系数,bxi是x方向第i倍频谐波的傅里叶余弦系数,ayi是y方向第i倍频谐波的傅里叶正弦系数,byi是y方向第i倍频谐波的傅里叶余弦系数,Ω是转子转速,t表示时间变量;
根据牛顿第二定律,径向磁力轴承-转子系统的磁悬浮转子动力学模型为:
Figure FDA0003068569830000013
Figure FDA0003068569830000014
其中,
Figure FDA0003068569830000015
表示磁悬浮转子x方向的所受的力,
Figure FDA0003068569830000016
表示磁悬浮转子y方向的所受的力,kx为x方向的位移刚度参数,ky为y方向的位移刚度参数,kix为x方向的当前刚度参数,kiy为y方向的当前刚度参数,m为转子质量,η为不平衡质量的偏心,ξ为不平衡质量的相位,ix(t)为x方向的磁轴承电流,iy(t)为y方向的磁轴承电流,t为时间变量;
步骤2),搭建带有以并行模式连接的多个相移陷波器的闭环控制系统,明确相移陷波器的插入方式;
步骤21),设计用于消除时变转速频率下控制电流中同频扰动分量的相移陷波器;
传统的陷波器具有如下的传递函数形式:
Figure FDA0003068569830000021
其中,N(s)为传统陷波器的传递函数,s为拉普拉斯算子,Ω表示转子转速,ε是相移陷波器的增益系数;
步骤22),在传统陷波器的基础上引入调节相位φ,构成相移陷波器,表示如下:
Figure FDA0003068569830000022
多相移陷波器并联表示为:
Figure FDA0003068569830000023
等价形式为:
Figure FDA0003068569830000024
其中,Nms(s)表示并联的多相移陷波器,n表示谐波次数,调节相位φ在时变转速频率下进行自我调节以保证整个闭环系统的稳定性,从而达到消除时变转速频率下控制电流中同频扰动分量的效果;
步骤3),根据步骤2)得到的相移陷波器设计具有相同相移的多个并联相移陷波器,稳定闭环系统;
步骤4),根据步骤2)搭建的闭环控制系统,对谐波电流抑制闭环系统进行稳定性分析,确定不同相位并联的多陷波器下整个系统的稳定条件;
相同相位的多个陷波器并联并不能实现全转速范围内的稳定从而不能在全转速范围进行有效的谐波电流抑制,因此设计不同相位并联的多陷波器,其传递函数表达如下:
Figure FDA0003068569830000025
它的等价形式为:
Figure FDA0003068569830000026
其中,Nmd(s)表示并联相移陷波器的传递函数,φi是第i倍频的相位;
步骤5),根据步骤2)、步骤3)和步骤4),利用有多个不同相位并联的相移滤波器对谐波电流进行有效抑制:
并行模式下不同相位并联的相移滤波器的抑制性能显示如下:
Figure FDA0003068569830000027
Figure FDA0003068569830000031
其中,k=1,2,...,n-1,n,Np(s)表示灵敏度函数,s表示拉普拉斯算子,ε表示相移陷波器的增益系数,Nmd(s)表示并联相移陷波器的传递函数,Ω表示转子转速,j表示虚数单位,
Figure FDA0003068569830000032
具有并行模式不同相位并联的相移滤波器的闭环系统中从e(s)到xd(s)的传递函数表示为:
Figure FDA0003068569830000033
将整个系统的特征方程表示为:
Figure FDA0003068569830000034
其中,GNd(s)表示闭环系统中从e(s)到xd(s)的传递函数,Gc(s)表示磁悬浮控制器,P(s)表示磁轴承转子,q(s)表示整个系统的特征方程,
在ε=0的情况下,根是s=±jΩ,±j2Ω,...,±jnΩ,在虚轴附近的根应该配置在左半复平面,对ε求导,满足下列式中所示的所有条件,以使整个系统是稳定的:
Figure FDA0003068569830000035
根据上式求导可推导出在全转速范围内不同频率下系统稳定的条件为:
Figure FDA0003068569830000036
根据基本磁悬浮控制系统的灵敏度函数的相角变化,对不同频率的转速进行分段处理,从而在各转速段选取满足上述稳定性条件的相移陷波器的调节相位值。
2.根据权利要求1所述多陷波器并联的磁悬浮转子振动谐波抑制方法,其特征在于:步骤3) 中通过具有相同相移的多个并联相移陷波器稳定闭环系统的方法:
将闭环系统特征方程写为:
1+Gc(s)P(s)=0
通过步骤2)明确多个并联具有相同相移的陷波器在闭环系统中的表达形式,因此进一步推导出插入陷波器后系统的特征方程为:
Figure FDA0003068569830000041
其中,s表示拉普拉斯算子,n表示谐波次数;
对于移相来说稳定条件为:
Figure FDA0003068569830000042
也就是说,移相同时满足以下条件:
-90°≤arg[S(jΩ)]+φ≤90°。
3.一种基于权利要求1所述的多陷波器并联的磁悬浮转子振动谐波抑制方法的抑制系统,其特征在于:包括依次连接在一起的并联相移陷波器Nmd(s)、磁悬浮控制器Gc(s)、磁轴承转子P(s),其中,并联相移陷波器Nmd(s)的信号提取点放在磁悬浮控制器Gc(s)的输入端,并联相移陷波器Nmd(s)的信号插入点放在位移误差信号处,参考位移信号x*(t)和位移传感器检测到的位移信号xs(t)得到位移误差信号,参考位移信号x*(t)、位移误差信号、以及并联相移陷波器Nmd(s)的信号进入到磁悬浮控制器Gc(s),磁悬浮控制器Gc(s)产生控制信号对磁轴承转子P(s)、扰动信号xd(s)进行控制,位移传感器检测到转子的位移信号xs(t);
并联相移陷波器Nmd(s)中:
不同相位并联的多陷波器传递函数如下:
Figure FDA0003068569830000043
它的等价形式为:
Figure FDA0003068569830000044
并行模式下不同陷波器的抑制性能显示如下:
Figure FDA0003068569830000045
Figure FDA0003068569830000051
其中,
Figure FDA0003068569830000052
具有并行模式不同相位的多陷波器的闭环系统中从e(s)到xd(s)的传递函数表示为:
Figure FDA0003068569830000053
将整个系统的特征方程表示为:
Figure FDA0003068569830000054
在ε=0的情况下,根是s=±jΩ,±j2Ω,...,±jnΩ,在虚轴附近的根应该配置在左半复平面上,对ε求导,满足下列式中所示的所有条件,以使整个系统是稳定的:
Figure FDA0003068569830000055
根据上式求导可推导出在全转速范围内不同频率下系统稳定的条件为:
Figure FDA0003068569830000056
其中:Np(s)表示灵敏度函数,s表示拉普拉斯算子,ε表示相移陷波器的增益系数,Nmd(s)表示并联相移陷波器的传递函数,Ω表示转子转速,j表示虚部单位。
CN201910254380.4A 2019-03-31 2019-03-31 多陷波器并联的磁悬浮转子振动谐波抑制方法及系统 Active CN109976403B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910254380.4A CN109976403B (zh) 2019-03-31 2019-03-31 多陷波器并联的磁悬浮转子振动谐波抑制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910254380.4A CN109976403B (zh) 2019-03-31 2019-03-31 多陷波器并联的磁悬浮转子振动谐波抑制方法及系统

Publications (2)

Publication Number Publication Date
CN109976403A CN109976403A (zh) 2019-07-05
CN109976403B true CN109976403B (zh) 2021-08-06

Family

ID=67081904

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910254380.4A Active CN109976403B (zh) 2019-03-31 2019-03-31 多陷波器并联的磁悬浮转子振动谐波抑制方法及系统

Country Status (1)

Country Link
CN (1) CN109976403B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112432634B (zh) * 2020-11-30 2023-03-14 北京航空航天大学 一种基于多同步旋转坐标变换的谐波振动力抑制方法
CN113741181A (zh) * 2021-01-18 2021-12-03 南京航空航天大学 一种转速自适应的磁悬浮转子系统奇次谐波电流抑制方法
CN113029440B (zh) * 2021-03-23 2022-09-23 华驰动能(北京)科技有限公司 一种磁悬浮转子动不平衡振动的抑制方法
CN114826406B (zh) * 2021-05-27 2024-01-09 中航光电科技股份有限公司 一种基于磁悬浮技术的超宽带射频光纤稳相传输系统
CN113485472B (zh) * 2021-08-04 2022-02-11 北京航空航天大学 一种基于双通道陷波器的磁悬浮转子同频振动力矩抑制方法
CN115514281B (zh) * 2022-09-27 2023-08-29 华驰动能(北京)科技有限公司 磁悬浮转子电流抑制方法、装置、电子设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6239772B1 (en) * 1998-02-17 2001-05-29 U.S. Philips Corporation Video moire reduction
CN102305886A (zh) * 2011-05-31 2012-01-04 浙江大学 电网电压谐波畸变及不平衡时基波电压同步信号检测方法
CN103322923A (zh) * 2013-06-09 2013-09-25 中国科学院长春光学精密机械与物理研究所 基于方波均值法消除非线性误差的光外差干涉法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104503238B (zh) * 2014-12-15 2017-03-22 北京航空航天大学 一种基于自适应重复控制器的磁悬浮转子系统电流谐波抑制方法
CN105116935B (zh) * 2015-08-14 2017-03-22 北京航空航天大学 一种基于频域变参数lms算法和陷波器的磁悬浮转子谐波振动抑制方法
WO2017040223A1 (en) * 2015-08-28 2017-03-09 Skyworks Solutions, Inc. Tunable notch filter and contour tuning circuit
CN106773673A (zh) * 2016-11-23 2017-05-31 北京航空航天大学 一种基于频率自适应的分数补偿重复控制器的磁悬浮转子谐波电流抑制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6239772B1 (en) * 1998-02-17 2001-05-29 U.S. Philips Corporation Video moire reduction
CN102305886A (zh) * 2011-05-31 2012-01-04 浙江大学 电网电压谐波畸变及不平衡时基波电压同步信号检测方法
CN103322923A (zh) * 2013-06-09 2013-09-25 中国科学院长春光学精密机械与物理研究所 基于方波均值法消除非线性误差的光外差干涉法

Also Published As

Publication number Publication date
CN109976403A (zh) 2019-07-05

Similar Documents

Publication Publication Date Title
CN109976403B (zh) 多陷波器并联的磁悬浮转子振动谐波抑制方法及系统
CN106647843B (zh) 一种基于复合分数重复控制器的磁悬浮转子谐波电流抑制方法
Chen et al. Unbalance vibration suppression for AMBs system using adaptive notch filter
CN109058292B (zh) 一种新型磁悬浮轴承不平衡振动力直接抑制方法
CN112432634B (zh) 一种基于多同步旋转坐标变换的谐波振动力抑制方法
CN104503238B (zh) 一种基于自适应重复控制器的磁悬浮转子系统电流谐波抑制方法
CN108227499B (zh) 一种谐波振荡器并联分数阶重复控制器的磁悬浮转子谐波电流抑制方法
CN108897218B (zh) 一种基于混合奇次重复控制器的磁悬浮转子奇次谐波振动抑制方法
JPH0572177B2 (zh)
JP2005502004A (ja) 同期擾乱の自動的補償のために使用される装置および方法
CN111624888B (zh) 基于重复控制器和补偿的磁悬浮转子振动力抑制方法及系统
CN113741181A (zh) 一种转速自适应的磁悬浮转子系统奇次谐波电流抑制方法
CN114637202A (zh) 基于扰动补偿的电磁轴承柔性转子系统同频振动抑制方法
Bu et al. Unbalanced displacement LMS extraction algorithm and vibration control of a bearingless induction motor
CN114326409B (zh) 基于双通道谐波重构的磁悬浮转子直接振动力抑制方法
CN115001336A (zh) 一种全转速磁悬浮高速电机转子振动控制方法
Okubo et al. Unbalance vibration control for active magnetic bearing using automatic balancing system and peak-of-gain control
Nakamura et al. A phase stabilization method for unbalance vibration control of five-axes active magnetic bearing systems
Kezhen et al. Harmonic vibration suppression of maglev rotor system under variable rotational speed without speed measurement
CN114371622B (zh) 基于多谐波逆Park变换的磁悬浮转子谐波振动力抑制方法
CN109687799B (zh) 基于假定旋转坐标的磁轴承转子分层振动补偿方法
Ren et al. Rotation modes stability analysis and phase compensation for magnetically suspended flywheel systems with cross feedback controller and time delay
CN112525179B (zh) 一种基于虚拟采样二阶重复控制的谐波电流抑制方法
Nakamura et al. An interpretation of unbalance vibration compensator for five-axes active magnetic bearing systems based on internal model principle
CN114962450B (zh) 磁悬浮转子系统同步振动抑制方法、系统、存储介质及终端

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