CN107133387A - 转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法 - Google Patents
转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法 Download PDFInfo
- Publication number
- CN107133387A CN107133387A CN201710229607.0A CN201710229607A CN107133387A CN 107133387 A CN107133387 A CN 107133387A CN 201710229607 A CN201710229607 A CN 201710229607A CN 107133387 A CN107133387 A CN 107133387A
- Authority
- CN
- China
- Prior art keywords
- mrow
- rotor
- msub
- search
- unbalance
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02K—DYNAMO-ELECTRIC MACHINES
- H02K15/00—Methods or apparatus specially adapted for manufacturing, assembling, maintaining or repairing of dynamo-electric machines
- H02K15/16—Centering rotors within the stator; Balancing rotors
- H02K15/165—Balancing the rotor
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Manufacturing & Machinery (AREA)
- Power Engineering (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
本发明公开了一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法。包括信号处理、不平衡系数辨识和补偿输出:信号处理是根据位移传感器获取的转子振动信号与转子转速及相位检测装置获取的转子转速和相位信号,计算得到转子振动信号中基频量的幅值;不平衡系数辨识是以转子振动信号中基频量的幅值为判断依据,通过变步长多边形迭代搜寻方法辨识获得不平衡系数;补偿输出是依据不平衡系数产生不平衡补偿电流,输入到电磁线圈中对转子的不平衡振动进行抑制。本发明可应用于带有主动控制单元的主动转子系统,实现对主动转子系统的不平衡进行在线补偿,以对转子系统在整个转速范围内加减速运行过程中的不平衡振动进行抑制。
Description
技术领域
本发明涉及了一种用于主动转子系统中的补偿算法,尤其是涉及了一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法。
背景技术
转子由材质不均匀、加工误差、变形等原因,其质量中心与几何中心之间总是存在着一定的偏差,这个偏差称为残余不平衡距。在转子的旋转过程中,转子的残余不平衡将产生一个与转速同频的不平衡激励力,导致转子产生振动。不平衡振动通过轴承传递给基座,严重时会导致设备损坏。不同于传统转子系统,带有主动控制单元的主动转子系统,如主动电磁轴承转子系统,能够给转子施加一个可控外力,从而能够对转子的不平衡振动进行主动控制,这是主动转子系统最重要也是独一无二的优势之一。
主动转子系统不平衡的补偿,特别是主动电磁轴承-刚性转子系统的不平衡补偿,已经开展了多年的研究。目前主要的方法可以分为两类:一类是采用陷波滤波器提取转子振动信号中的同频量,然而由于陷波滤波器需要设置一个中心频率,并且只能在此频率附近工作,因此不能应用于转子转速发生变化的场合,特别是整个转速区间内转子加减速运行过程中不平衡振动的抑制;另一类是采用自适应控制法,如基于LMS算法的自适应滤波器、神经网络算法等,然而这些算法大多复杂,需要昂贵的先进控制器,增加了系统成本;另外如基于LMS算法的自适应滤波器,其算法中的迭代步长为固定值,不能在保证算法的辨识精度的同时加快算法的收敛速度。因此在系统整个转速区间内加减速运行的场合,应用效果不佳。
因此,现有技术中缺少一种方法,既能在转子不同转速下起作用,算法又简单易实现,同时实时性又能满足要求,能实现转子系统在整个转速区间内加减速运行过程中不平衡振动抑制的同时又降低系统成本。
发明内容
为了克服现有主动转子系统在转子位移高精度控制方面的不足,特别是解决整个转速范围内不平衡振动控制的问题,本发明提供一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,通过在线辨识转子不平衡系数,对主动转子系统的不平衡进行主动补偿,以实现在整个转速范围内不平衡振动的抑制。
如图1所示,本发明的技术解决方案是:
为了克服现有主动转子系统在转子位移高精度控制方面的不足,特别是解决整个转速范围内不平衡振动控制的问题,本发明提供一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,通过在线辨识转子不平衡系数,对主动转子系统的不平衡进行主动补偿,以实现在整个转速范围内不平衡振动的抑制。
本发明的技术解决方案是:
所述的不平衡补偿算法包括实时并依次进行的信号处理、不平衡系数辨识和补偿输出步骤:
信号处理是根据位移传感器获取的转子振动信号e(kT)与转子转速及相位检测装置获取的转子转速ω和相位信号计算得到转子两端的转子振动信号e(kT)中基频量的幅值E(k);
不平衡系数辨识是以转子振动信号e(kT)中基频量的幅值E(k)为判断依据,通过变步长的多边形迭代搜寻方法,辨识计算后输出转子的不平衡系数(α,β)的目标值;
补偿输出是依据不平衡系数(α,β)的目标值产生不平衡补偿电流ic,用不平衡补偿电流ic输入到电磁线圈中对转子的不平衡振动进行抑制。
具体实施中,所述位移传感器采用非接触式位移传感器,转子转速及相位检测装置采用编码盘及光电或磁电开关。
所述的信号处理是依据互相关计算方法,将转子振动信号与依据转子振动信号得到的基频正余弦函数做互相关分析,得到当前转子振动信号中基频成分的两个正交分量,再由两个正交分量相叠加获得基频量的幅值E(k)。
基频正余弦函数中的基频和转子转速ω对应的频率相同。
所述信号处理针对转子其中一端,在离散采样下获得的转子振动信号具体是:
先建立垂直于转子轴向的xy平面坐标系,x方向和y方向与所述位移传感器所采集的转子振动信号中转子振动的方向对应一致;
然后采用以下公式计算获得转子其中一端沿x方向的基频振动信号的幅值:
其中,k表示采样周期数,T表示采样周期,N表示采样个数,ω为转子的旋转角速度,eax表示转子其中一端离散采集到的沿x方向的振动信号,aax表示转子其中一端离散采集到的沿x方向的振动信号eax中基频成分的正弦分量,bax表示转子其中一端离散采集到的沿x方向的振动信号eax中基频成分的余弦分量;
同时采用以下公式计算获得转子其中一端沿y方向的基频振动信号的幅值:
其中,eay表示转子其中一端离散采集到的沿y方向的振动信号,aay表示转子其中一端离散采集到的沿y方向的振动信号eay中基频成分的正弦分量,bay表示转子其中一端离散采集到的沿y方向的振动信号eay中基频成分的余弦分量;
再采用以下公式计算获得转子其中一端振动信号中的基频量的幅值E(k):
其中,E(k)表示转子其中一端振动信号中的第k个采样周期时刻的基频量的幅值。
所述基频振动信号的基频与转子转速ω对应的频率相同。
所述的不平衡系数辨识是对不平衡系数进行辨识,设置辨识值的初始值(α,β),辨识值的初始值任意取值,通常取(0,0),从辨识值的初始值开始,采用多边形迭代搜寻方法实时调整方向进行搜寻,搜寻以幅值E(k)最小为目标,在多边形迭代搜寻进行的同时采用变步长算法实时改变搜寻步长,使辨识值无限接近目标值,最终输出的辨识值作为目标值作为最终的不平衡系数。
所述的多边形迭代搜寻具体是初始时刻用任意搜寻步长沿任意方向开始搜寻,每一步搜寻后对方向进行判断调整再继续进行下一步搜寻:
当判断搜寻方向正确时,不改变方向;
当判断搜寻方向错误时,改变方向;
所述的变步长算法是在多边形迭代搜寻每一步搜寻后对方向进行判断调整的同时,对步长也进行判断调整,再继续进行下一步搜寻:
当判断搜寻方向正确时,减小步长;
当判断搜寻方向错误时,增大步长;通过改变步长能同时保证迭代搜寻算法的辨识精度与收敛速度。
改变步长计算为:
R(k)=R(0)·E(k)/Ec
其中,R(k)为当前步长,R(0)为初始步长,Ec为搜寻启动阈值。
通过上述方式使辨识值无限接近目标值。
所述的多边形迭代搜寻中,改变方向是以固定角度进行改变,所使用的固定角度是0到360度之间的任意值。例如取120度时,搜寻轨迹为三角形;如取90度时,搜寻轨迹为四边形;取60度时,搜寻轨迹为六边形。
所述的多边形迭代搜寻和变步长算法中,搜寻方向判断方式为:
如果第k个采样周期时刻的幅值E(k)≤上一个(第k-1个)采样周期时刻的幅值E(k-1),则判断搜寻方向正确,否则判断搜寻方向错误。
本发明只有当当前转子振动信号中基频量的幅值E(k)大于Ec时,算法才启动,则若要求辨识值与目标值完全相等,Ec理论上取值为零,不过为避免算法频繁启动导致频繁占用控制系统资源,同时又能保证算法的精度,通常将其取为一个很小的正数。
所述补偿输出是将不平衡系数辨识输出转子两端的不平衡系数代入以下公式计算获得补偿电流:
其中,αa表示转子不平衡系数在其中一端α轴方向的分量,βa表示转子不平衡系数在其中一端β轴方向的分量,αb表示转子不平衡系数在另一端α轴方向的分量,βb表示转子不平衡系数在另一端β轴方向的分量,这些分量均由不平衡系数辨识得到;m表示转子质量,l表示两端电磁轴承的距离,kia表示其中一端电磁轴承的电流刚度系数,kib表示另一端电磁轴承的电流刚度系数,表示当前转子相位,表示当前转子转速,表示当前转子转速的变化率。
所述转子振动信号e(kT)中基频量的幅值E(k)对应的基频替换为整数倍转子转速ω对应频率fi下的幅值Ei(k),同样采用上述方法能对频率为fi的振动进行补偿控制。转子振动信号e(kT)中基频量的幅值E(k)对应的基频和转子转速ω对应的频率相同。
本发明的不平衡补偿方法可以应用于带有主动控制单元的主动转子系统,例如主动电磁轴承-转子系统,原理上是利用特殊的不平衡补偿算法,对转子不平衡系数进行在线辨识,并结合转子系统的转速与相位检测装置以及转子系统的动力学模型,实现对主动转子系统的不平衡进行在线补偿,以对转子系统在整个转速范围内加减速运行过程中的不平衡振动进行抑制。
本发明与现有技术相比的优点在于:
(1)本发明采用迭代搜寻方法对转子不平衡进行辨识,克服了闭环系统的稳定性问题。
(2)本发明中的变步长算法既能保证算法的辨识精度,又能加快其收敛速度。
(3)本发明中的多边形迭代搜寻方法,算法简单,易于实现,在现有的常规数字控制器中即可运行,因此能有效降低系统成本。
(4)本发明中的算法启动阈值可以有效避免算法频繁启动的同时保证算法精度,占控制系统资源小。
(5)本发明不仅适用于恒转速运行工况,也适用于转子在加减速运行过程中对不平衡振动的控制。
(6)如果把基频及转子振动信号中基频量的幅值E(k)变为其它频率及对应的幅值,也可以对其它频率的振动实现主动补偿控制。
本发明用输出补偿电流来抑制转子的不平衡振动,特别适用于带有主动控制单元的主动转子系统,如主动电磁轴承-转子系统,在整个工作转速范围内不平衡振动的主动抑制。
附图说明
图1为本发明方法的原理框图。
图2为实施例所用的主动电磁轴承-刚性转子系统不平衡补偿控制系统;
图3为本发明中多边形搜寻法,以三角形为例,以转子某一端为例的搜寻过程示意图。(另一端相同)
图4为本发明算法以转子某一端为例的流程框图。
图5为实施例某主动电磁轴承-刚性转子系统在1337rpm的临界转速附近采用本发明方法前后获得的两个电磁轴承位置转子的运动轨迹。
图6为实施例某主动电磁轴承-刚性转子系统在1337rpm的临界转速附近采用本发明方法的搜寻轨迹。
图7为某主动电磁轴承-刚性转子系统在加速运行过程中采用本发明算法前后转子某端某方向上振动响应曲线。
具体实施方式
下面以主动电磁轴承-刚性转子系统为例对一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法进行具体说明。
本实施例所用的主动电磁轴承-刚性转子系统如图2所示。
转子的轴向由上下一对永磁轴承支承,径向分别由上下电磁轴承支承,转子的位移用位于上下电磁轴承附近的非接触式的位移传感器来进行测量。
如图2,首先可以定义上电磁轴承所在平面以及下电磁轴承所在平面分别为Π1和Π2,上下电磁轴承轴颈几何中心分别位于Π1和Π2内,记为O1、O2点,O1O2两点连线称为定子几何中心线。转子为轴对称转子,因此存在转子几何中心线。
不考虑转子的质量不平衡时,转子质心位于转子几何中心线上,假设位于图中的G点,过G点做平行于Π1、Π2的平面Π,Π(也称为转子的质心横截面)与定子几何中心线交于O点。假设上下电磁轴承所在平面到转子质心横截面的距离分别为la和lb,上下位移传感器所在平面到转子质心横截面的距离分别为lsa和lsb。
为了描述转子的运动状态,以O为原点,建立固定右手坐标系O_xyz,其中Oz轴以沿定子几何中心线朝上的方向为正,Ox轴及Oy轴方向如图2所示。转子的轴向运动由永磁轴承所限定,此处不予考虑。于是,转子在空间的位置可以用G点在Oxy坐标系中的位移(x,y)以及转子几何中心轴绕x和y轴逆时针旋转角度(α,β)唯一确定。转子绕x和y轴旋转角度的正方向如图2所示,因此,转子质心G点的广义坐标为qo=[x β y -α]T,其中[]T符号表示矩阵的转置。
在上下电磁轴承所在平面内,分别以O1及O2点为原点,建立电磁轴承固定坐标系O1_x1y1z1及O2_x2y2z2。电磁轴承固定坐标系各轴的方向与定子固定坐标系O_xyz相同。在上下电磁轴承位置,轴颈几何中心O1及O2分别由坐标(xa,ya)和(xb,yb)来表示,由此构成电磁轴承位置轴颈的坐标矢量为qb=[xa ya xb yb]T。
同理,可在上下传感器所在平面建立对应的固定坐标系,用坐标(xsa,ysa)和(xsb,ysb)表示上下传感器所在平面轴颈几何中心的位置,则传感器位置轴颈的坐标矢量为qs=[xsa ysa xsb ysb]T。
实际中,用传感器测量得到的转子位移是传感器位置轴颈的位置矢量qs,而电磁轴承及转子动力学模型中使用的却是平衡转子质心G点的广义坐标矢量qo及电磁轴承位置轴颈的坐标矢量qb,利用转子的几何关系可以三个坐标之间的关系为:
其中各变量的解释均已在上文对图2的说明中给出。
实际的转子系统是不平衡的,为了便于建立不平衡转子系统的动力学模型,我们把转子的不平衡看作为一个在平衡转子系统质心G以外的附加质量mu造成的。这个不平衡质量与转子的质量相比很小,对转子质心位置的影响可以忽略不计。
假设转子的不平衡质量mu不在转子的质心横截面Π内,而是位于转子质心横截面Π的上端C’点,C’点在平面Π内的投影为C,|CC'|=ez,|GC|=eu。为了描述Π平面内C点的位置,以G为原点建立Gαβ旋转直角坐标系,在初始时刻,旋转坐标系Gαβ的Gα轴与Ox轴平行,在t时刻,Gα轴与Ox轴的夹角为旋转坐标系的旋转速度为ω(t)。
由转子动力学理论,可以得到轴对称电磁轴承-刚性转子系统的动力学方程为:
用矩阵形式表示如下:
其中,M=diag(m,Jr,m,Jr)为转子系统的广义质量矩阵,为转子系统的陀螺矩阵,为坐标变换矩阵,为广义力系数矩阵,fm=[fxa,fya,fxb,fyb]T为上下电磁轴承的电磁力矢量,fu=[fux,fuy,fux,fuy]T为转子的不平衡力矢量。m为转子质量,Jr为转子关于x或y轴的转动惯量,Jz为转子关于z轴的转动惯量,ω为转子的旋转角速度,fxa、fya、fxb、fyb分别为转子两端在x及y方向上受到的电磁力,fux、fuy分别为转子在x及y方向上受到的不平衡力,其它变量的解释已在上文对图2的说明中给出。
如果电磁轴承电磁力采用线性化模型,则有:
fm=Ksqb+Kii (4)
其中,Ks=diag(ksa,ksa,ksb,ksb)为AMB的位移刚度矩阵,Ki=diag(kia,kia,kib,kib)为电磁轴承的电流刚度矩阵;i=[ixa,iya,ixb,iyb]T为电磁轴承线圈中的电流矢量;ksa、kia、ksb、kib分别为对应电磁轴承的位移刚度系数及电流刚度系数,ixa、iya、ixb、iyb分别为两端在x及y方向上电磁线圈中的电流。
变速运行过程中,转子质心在x及y方向上受到的不平衡力可以表示为:
其中,Uu=mueu为转子系统的不平衡量,为不平衡质量mu与不平衡质量到平衡转子几何中心距离eu的乘积。通常可表示为Uu=mueu=mem,这里em为等效转子质量偏心量。θ为C点在Gαβ中位置角。
结合式(1)、式(3)、式(4)与式(5),本实施例所描述的由轴承坐标矢量描述的径向四自由度主动电磁轴承-刚性转子系统的数学模型如式(6)所示,其中,不平衡力fu被视为作用在转子质心的外力。
不平衡补偿的目的是通过补偿器输出补偿电流ic,经过电磁轴承线圈后,产生恰好抵消不平衡力的电磁力。当达到完全补偿时,转子将不受任何外力作用,仅仅围绕其几何中心线自转,并且几何中心线与定子几何中心线重合,即对应式(6)中:qb=0。此时,由式(6)可以得到电磁轴承线圈中的补偿电流为:
ic=-(LmKi)-1Lufu (7)
其具体表达式为:
其中, 和称为转子两端的不平衡系数,统称为转子的不平衡系数,用来表示,其中表示转子不平衡系数在a端α轴方向的分量,表示转子不平衡系数在a端β轴方向的分量,表示转子不平衡系数在b端α轴方向的分量,表示转子不平衡系数在b端β轴方向的分量;转子a端和b端分别指的是转子的两端。
不平衡系数和表示不平衡量在转子中的相对位置,采用本发明特定的方法能够辨识获得不平衡系数和来对转子不平衡进行补偿。
以转子a端为例进行说明,本发明的不平衡补偿算法1包括实时并依次进行的信号处理2、不平衡系数辨识3和补偿输出4步骤,具体实施例及其实施过程如下:
1.振动信号基频幅值获取与算法启动阈值
采用信号处理技术中的相关法,可准确地提取出振动信号中基频的幅值。
对于连续采集的转子振动信号,其中含有多种成份,可以表示为如下形式:
其中,an为n倍转速信号的幅值,nω为n倍转速信号的角频率,为n倍转速信号的初始相位,s(t)为干扰信号。
将转子振动信号Y(t)分别与依据转子振动信号得到的基频正余弦函数做互相关,利用周期函数之间的正交性,就可以得到转子振动信号中基频成分的两个正交分量分别为:
其中,T为积分时间。Ra表示转子振动信号中基频成分的正弦分量,Rb表示转子振动信号中基频成分的余弦分量;
因此,基频信号的幅值为
于是,在本发明离散采集的离散控制系统中,a端电磁轴承x方向基频振动信号幅值可以表示为:
其中,k表示采样周期数,T表示采样周期,N是一个正整数,一般取为光码盘码数的整数倍;eax是a端电磁轴承x方向上的振动信号。
同理,a端电磁轴承y方向上的振动信号基频幅值可以表示为:
其中,eay是a端电磁轴承y方向上的振动信号。
转子a端振动总的基频幅值为以此表征当前采样时刻转子在a端电磁轴承平面内振动的程度。为避免控制器过于频繁执行,当该值超过某一个设定值Eac时,算法才开始启动。因此Eac也称为算法在a端的启动阈值。
2.多边形搜寻法
从辨识值的初始值开始,通常取(0,0),采用多边形迭代搜寻3-1方法实时调整方向进行搜寻,多边形搜寻法的具体过程如下:
(1)初始化:(αa,βa)0=(0,0),r0=R(0)∠0,其中,R(0)为初始步长,r0为初始搜寻矢量。
(2)第k步:(αa,βa)k=(αa,βa)k-1+rk-1;判断Ea(k)≤Ea(k-1)?如满足,则rk=rk-1+R(k)∠0,否则,rk=rk-1+R(k)∠φ。其中,R(k)为当前步长,rk为当前搜寻矢量,φ为多边形搜寻法中搜寻方向角度的改变量,k为搜寻的步数,也是当前采样周期数,k=1,2,3,4,…。
当多边形搜寻法中搜寻方向角度改变量φ=120度时,就为三角形搜寻过程。
下面以图3所示的三角形搜寻过程为例对这个搜寻过程进行详细说明。
●当k=1,(αa,βa)1=(αa,βa)0+r0,搜寻至图中A点,判断Ea(1)≤Ea(0)?如满足,则r1=r0+R(1)∠0,保持该方向;如不满足,则r1=r0+R(1)∠(2π/3),改变方向;
●当k=2,(αa,βa)2=(αa,βa)1+r1,搜寻至图中B点,判断Ea(2)≤Ea(1)?如满足,则r2=r1+R(2)∠0,保持该方向;如不满足,则r2=r1+R(2)∠(2π/3),改变方向;
●……
●当k=4,(αa,βa)4=(αa,βa)3+r3,搜寻至图中D点,判断Ea(4)≤Ea(3)?如满足,则r4=r3+R(4)∠0,保持该方向;如不满足,则r4=r3+R(4)∠(2π/3),改变方向;
●当k=5,(αa,βa)5=(αa,βa)4+r4,搜寻至图中E点,判断Ea(5)≤Ea(4)?如满足,则r5=r4+R(5)∠0,保持该方向;如不满足,则r5=r4+R(5)∠(2π/3),改变方向;
●……
●当k=7,(αa,βa)7=(αa,βa)6+r6,搜寻至图中G点,判断Ea(7)≤Ea(6)?如满足,则r7=r6+R(7)∠0,保持该方向;如不满足,则r7=r6+R(7)∠(2π/3),改变方向;
●当k=8,(αa,βa)8=(αa,βa)7+r7,搜寻至图中H点,判断Ea(8)≤Ea(7)?如满足,则r8=r7+R(8)∠0,保持该方向;如不满足,则r8=r7+R(8)∠(2π/3),改变方向;
●……
最终,辨识值(αa,βa)在目标值附近,稳定在边长为R(k)的三角形内。可见,该算法的精确度为R(k)。
3.变步长算法
由以上的分析可以发现,增大步长可以加快搜寻进程,即增加算法的收敛速度;减小步长可以提高算法的精度。显然,对于固定步长,算法的收敛速度与算法的精度之间是互相矛盾的。因此,在多边形迭代搜寻3-1进行的同时采用变步长算法3-2实时改变搜寻步长,使辨识值无限接近目标值,以达到在不降低精度的前提下,增大算法的收敛速度。
因此,以a端为例,变步长算法为:
R(k)=R(0)·Ea(k)/Eac (13)
其中,R(k)为当前步长,R(0)为初始步长,Ea(k)为转子a端振动信号中基频幅值。
此时,可将R(0)取较小值作为试探,当搜寻方向错误时,改变方向以更大的步长进行搜寻,以增加收敛速度。当搜寻方向正确时,步长逐渐减小,使搜寻保持在一个很小的平面范围内,以提高算法的精度。
如图4所示,首先初始化,令搜寻起点从原点开始((αa,βa)0=(0,0),r0=R(0)∠0),然后进入搜寻循环体;循环体中,首先判断是否满足启动算法的要求,一旦不满足(Ea≤Eac),则不开启算法,退出循环并结束;如果满足(Ea>Eac),则继续搜寻。每迭代一次后,都需要进行判断本次迭代搜寻方向是否正确,如果正确(Ea(k)≤Ea(k-1)),则不改变方向继续搜寻;如果错误(Ea(k)>Ea(k-1)),则改变方向(改变角度为φ)继续搜寻。
本实施例某主动电磁轴承-刚性转子系统在转速1337rpm的临界转速附近采用本发明前后的上下电磁轴承位置转子的运动轨迹如图5所示,搜寻轨迹如图6所示。图中可见,算法对转子振动的抑制效果是非常明显,搜寻轨迹也从原点出发,朝目标值进行。
由于在不同的转速下转子系统不平衡补偿控制器所采用的步长应该是不同的,如果采用了定步长的控制算法,那么就会出现在一些转速区转子的振动得到了良好的控制,在另外一些转速区就难以得到满意的效果,特别是在转子的加减速运行过程中,由于每一个转速下控制时间是非常短的,这就要求算法不仅要有高的收敛精度,还需要有很快的收敛速度。而采用本发明提出的变步长搜寻算法的不平衡补偿控制器既具有高的收敛精度,还具有很快的收敛速度的特点。本实施例某主动电磁轴承-刚性转子系统在加速过程中采用本发明算法前后,转子系统在某端某方向上振动响应曲线对比如图7所示。图中可见,转子在加速过程中,算法对(尤其在临界转速附近)转子同频振动的抑制效果非常明显,能够保证系统在宽转速范围内的稳定运行。
Claims (10)
1.一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,其特征在于:所述的不平衡补偿算法(1)包括实时并依次进行的信号处理(2)、不平衡系数辨识(3)和补偿输出(4)步骤:
信号处理(2)是根据位移传感器获取的转子振动信号e(kT)与转子转速及相位检测装置获取的转子转速ω和相位信号计算得到转子两端的转子振动信号e(kT)中基频量的幅值E(k);
不平衡系数辨识(3)是以转子振动信号e(kT)中基频量的幅值E(k)为判断依据,通过变步长的多边形迭代搜寻方法,辨识计算后输出转子的不平衡
补偿输出(4)是依据不平衡系数产生不平衡补偿电流ic,用不平衡补偿电流ic输入到电磁线圈中对转子的不平衡振动进行抑制。
2.根据权利要求1所述的一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,其特征在于:所述的信号处理(2)是依据互相关计算方法,将转子振动信号与依据转子振动信号得到的基频正余弦函数做互相关分析,得到当前转子振动信号中基频成分的两个正交分量,再由两个正交分量相叠加获得基频量的幅值E(k)。
3.根据权利要求1所述的一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,其特征在于:所述信号处理(2)针对转子其中一端,在离散采样下获得的转子振动信号具体是:
先建立垂直于转子轴向的xy平面坐标系,x方向和y方向与所述位移传感器所采集的转子振动信号中转子振动的方向对应一致;
然后采用以下公式计算获得转子其中一端沿x方向的基频振动信号的幅值:
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>a</mi>
<mrow>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>b</mi>
<mrow>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>e</mi>
<mrow>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mi>k</mi>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</mrow>
<mrow>
<msub>
<mi>b</mi>
<mrow>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>e</mi>
<mrow>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mi>k</mi>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</mrow>
其中,k表示采样周期数,T表示采样周期,N表示采样个数,ω为转子的旋转角速度,eax表示转子其中一端离散采集到的沿x方向的振动信号,aax表示转子其中一端离散采集到的沿x方向的振动信号eax中基频成分的正弦分量,bax表示转子其中一端离散采集到的沿x方向的振动信号eax中基频成分的余弦分量;
同时采用以下公式计算获得转子其中一端沿y方向的基频振动信号的幅值:
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>a</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>a</mi>
<mrow>
<mi>a</mi>
<mi>y</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>b</mi>
<mrow>
<mi>a</mi>
<mi>y</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>a</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>e</mi>
<mrow>
<mi>a</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mi>k</mi>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</mrow>
<mrow>
<msub>
<mi>b</mi>
<mrow>
<mi>a</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>e</mi>
<mrow>
<mi>a</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mi>k</mi>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</mrow>
其中,eay表示转子其中一端离散采集到的沿y方向的振动信号,aay表示转子其中一端离散采集到的沿y方向的振动信号eay中基频成分的正弦分量,bay表示转子其中一端离散采集到的沿y方向的振动信号eay中基频成分的余弦分量;
再采用以下公式计算获得转子其中一端振动信号中的基频量的幅值E(k):
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>A</mi>
<mrow>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>A</mi>
<mrow>
<mi>a</mi>
<mi>y</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
其中,E(k)表示转子其中一端振动信号中的第k个采样周期时刻的基频量的幅值。
4.根据权利要求3所述的一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,其特征在于:所述基频振动信号的基频与转子转速ω对应的频率相同。
5.根据权利要求3所述的一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,其特征在于:所述的不平衡系数辨识(3)是对不平衡系数进行辨识,设置辨识值的初始值(α,β),从辨识值的初始值开始,采用多边形迭代搜寻(3-1)方法实时调整方向进行搜寻,搜寻以幅值E(k)最小为目标,在多边形迭代搜寻(3-1)进行的同时采用变步长算法(3-2)实时改变搜寻步长,使辨识值无限接近目标值,最终输出的辨识值(α,β)作为目标值作为最终的不平衡系数。
6.根据权利要求5所述的一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,其特征在于:
所述的多边形迭代搜寻(3-1)具体是初始时刻用任意搜寻步长沿任意方向开始搜寻,每一步搜寻后对方向进行判断调整再继续进行下一步搜寻:
当判断搜寻方向正确时,不改变方向;
当判断搜寻方向错误时,改变方向;
所述的变步长算法(3-2)是在多边形迭代搜寻(3-1)每一步搜寻后对方向进行判断调整的同时,对步长也进行判断调整,再继续进行下一步搜寻:
当判断搜寻方向正确时,减小步长;
当判断搜寻方向错误时,增大步长;
改变步长计算为:
R(k)=R(0)·E(k)/Ec
其中,R(k)为当前步长,R(0)为初始步长,Ec为搜寻启动阈值。
7.根据权利要求5所述的一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,其特征在于:所述的多边形迭代搜寻(3-1)中,改变方向的角度是固定值,是0到360度之间的任意值。
8.根据权利要求5所述的一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,其特征在于:所述的多边形迭代搜寻(3-1)和变步长算法(3-2)中,搜寻方向判断方式为:如果第k个采样周期时刻的幅值E(k)≤上一个采样周期时刻的幅值E(k-1),则判断搜寻方向正确,否则判断搜寻方向错误。
9.根据权利要求5所述的一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,其特征在于:所述补偿输出(4)是将不平衡系数辨识(3)输出转子两端的不平衡系数代入以下公式计算获得补偿电流:
其中,αa表示转子不平衡系数在其中一端α轴方向的分量,βa表示转子不平衡系数在其中一端β轴方向的分量,αb表示转子不平衡系数在另一端α轴方向的分量,βb表示转子不平衡系数在另一端β轴方向的分量;m表示转子质量,l表示两端电磁轴承的距离,kia表示其中一端电磁轴承的电流刚度系数,kib表示另一端电磁轴承的电流刚度系数,表示当前转子相位,表示当前转子转速,表示当前转子转速的变化率。
10.根据权利要求1所述的一种转子不平衡系数变步长多边形迭代搜寻的不平衡补偿算法,其特征在于:所述转子振动信号e(kT)中基频量的幅值E(k)对应的基频替换为整数倍转子转速ω对应频率fi下的幅值Ei(k),同样采用上述方法能对频率为fi的振动进行补偿控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710229607.0A CN107133387B (zh) | 2017-04-10 | 2017-04-10 | 转子不平衡系数变步长多边形迭代搜寻的不平衡补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710229607.0A CN107133387B (zh) | 2017-04-10 | 2017-04-10 | 转子不平衡系数变步长多边形迭代搜寻的不平衡补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107133387A true CN107133387A (zh) | 2017-09-05 |
CN107133387B CN107133387B (zh) | 2019-10-18 |
Family
ID=59716581
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710229607.0A Active CN107133387B (zh) | 2017-04-10 | 2017-04-10 | 转子不平衡系数变步长多边形迭代搜寻的不平衡补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107133387B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109725533A (zh) * | 2018-12-29 | 2019-05-07 | 南京工业大学 | 一种自适应迭代补偿控制方法和系统 |
CN109765781A (zh) * | 2018-12-29 | 2019-05-17 | 南京工业大学 | 一种同频扰动抑制方法和系统 |
CN110118632A (zh) * | 2018-02-07 | 2019-08-13 | 申克罗泰克有限责任公司 | 借助位移传感器测量轴弹性转子的不平衡度的方法 |
CN110469583A (zh) * | 2019-08-23 | 2019-11-19 | 广东美的暖通设备有限公司 | 磁轴承组件及控制方法、压缩机、空调器和存储介质 |
CN110762120A (zh) * | 2019-11-18 | 2020-02-07 | 南京航空航天大学 | 一种基于磁悬浮轴承转子系统的高回转精度控制方法 |
CN111177943A (zh) * | 2020-01-06 | 2020-05-19 | 浙江大学 | 一种转子不平衡系数变角度迭代搜寻的不平衡补偿方法 |
CN112803852A (zh) * | 2021-02-04 | 2021-05-14 | 江苏大学 | 变步长变角度搜索遗传算法优化转子振动补偿系统及方法 |
CN112886946A (zh) * | 2021-01-26 | 2021-06-01 | 浙江大学 | 电磁轴承主动控制的变相位自适应lms陷波器的设计方法 |
CN112953344A (zh) * | 2021-03-26 | 2021-06-11 | 江苏大学 | 一种无轴承异步电机转子不平衡振动补偿控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331565A (zh) * | 2014-11-10 | 2015-02-04 | 河海大学常州校区 | 轴类磁悬浮刚性转子系统的动力学建模方法及控制方法 |
CN104659990A (zh) * | 2015-03-11 | 2015-05-27 | 河南科技大学 | Lms自适应滤波无轴承电机的不平衡振动位移提取方法 |
CN104660137A (zh) * | 2015-03-11 | 2015-05-27 | 河南科技大学 | Lms自适应滤波无轴承电机的不平衡激振力补偿方法 |
CN105116935A (zh) * | 2015-08-14 | 2015-12-02 | 北京航空航天大学 | 一种基于频域变参数lms算法和陷波器的磁悬浮转子谐波振动抑制方法 |
CN105783898A (zh) * | 2016-04-28 | 2016-07-20 | 北京航空航天大学 | 一种基于频域自适应lms算法的磁悬浮转子谐波振动抑制方法 |
CN105978266A (zh) * | 2016-04-08 | 2016-09-28 | 武义益圣设备有限公司 | 电机转子平衡系统和方法 |
-
2017
- 2017-04-10 CN CN201710229607.0A patent/CN107133387B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331565A (zh) * | 2014-11-10 | 2015-02-04 | 河海大学常州校区 | 轴类磁悬浮刚性转子系统的动力学建模方法及控制方法 |
CN104659990A (zh) * | 2015-03-11 | 2015-05-27 | 河南科技大学 | Lms自适应滤波无轴承电机的不平衡振动位移提取方法 |
CN104660137A (zh) * | 2015-03-11 | 2015-05-27 | 河南科技大学 | Lms自适应滤波无轴承电机的不平衡激振力补偿方法 |
CN105116935A (zh) * | 2015-08-14 | 2015-12-02 | 北京航空航天大学 | 一种基于频域变参数lms算法和陷波器的磁悬浮转子谐波振动抑制方法 |
CN105978266A (zh) * | 2016-04-08 | 2016-09-28 | 武义益圣设备有限公司 | 电机转子平衡系统和方法 |
CN105783898A (zh) * | 2016-04-28 | 2016-07-20 | 北京航空航天大学 | 一种基于频域自适应lms算法的磁悬浮转子谐波振动抑制方法 |
Non-Patent Citations (4)
Title |
---|
JIANG KEJIAN,ZHU CHANGSHENG,CHEN LIANGLIANG: "Unbalance compensation by recursive seeking unbalance mass position in active magnetic bearing-rotor system", 《IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS》 * |
于洁,祝长生: "基于希尔伯特变换的自传感电磁轴承实现", 《浙江大学学报(工学版)》 * |
张德魁,江伟,赵鸿宾: "磁悬浮轴承系统不平衡振动控制的方法", 《清华大学学报(自然科学版)》 * |
蒋科坚, 祝长生: "基于不平衡识别的主动电磁轴承转子系统自动平衡", 《振动工程学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110118632A (zh) * | 2018-02-07 | 2019-08-13 | 申克罗泰克有限责任公司 | 借助位移传感器测量轴弹性转子的不平衡度的方法 |
CN109725533A (zh) * | 2018-12-29 | 2019-05-07 | 南京工业大学 | 一种自适应迭代补偿控制方法和系统 |
CN109765781A (zh) * | 2018-12-29 | 2019-05-17 | 南京工业大学 | 一种同频扰动抑制方法和系统 |
CN110469583A (zh) * | 2019-08-23 | 2019-11-19 | 广东美的暖通设备有限公司 | 磁轴承组件及控制方法、压缩机、空调器和存储介质 |
CN110762120A (zh) * | 2019-11-18 | 2020-02-07 | 南京航空航天大学 | 一种基于磁悬浮轴承转子系统的高回转精度控制方法 |
CN111177943A (zh) * | 2020-01-06 | 2020-05-19 | 浙江大学 | 一种转子不平衡系数变角度迭代搜寻的不平衡补偿方法 |
CN111177943B (zh) * | 2020-01-06 | 2022-03-22 | 浙江大学 | 一种转子不平衡系数变角度迭代搜寻的不平衡补偿方法 |
CN112886946A (zh) * | 2021-01-26 | 2021-06-01 | 浙江大学 | 电磁轴承主动控制的变相位自适应lms陷波器的设计方法 |
CN112803852A (zh) * | 2021-02-04 | 2021-05-14 | 江苏大学 | 变步长变角度搜索遗传算法优化转子振动补偿系统及方法 |
CN112803852B (zh) * | 2021-02-04 | 2022-05-20 | 江苏大学 | 变步长变角度搜索遗传算法优化转子振动补偿系统及方法 |
CN112953344A (zh) * | 2021-03-26 | 2021-06-11 | 江苏大学 | 一种无轴承异步电机转子不平衡振动补偿控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107133387B (zh) | 2019-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107133387B (zh) | 转子不平衡系数变步长多边形迭代搜寻的不平衡补偿方法 | |
Chuan et al. | Unbalance compensation for active magnetic bearing rotor system using a variable step size real-time iterative seeking algorithm | |
CN108716471B (zh) | 一种磁悬浮分子泵转子极小位移主动控制方法 | |
CN111177943B (zh) | 一种转子不平衡系数变角度迭代搜寻的不平衡补偿方法 | |
CN102169046B (zh) | 一种磁悬浮机电设备弹性模态在线测试系统 | |
CN104660137B (zh) | Lms自适应滤波无轴承电机的不平衡激振力补偿方法 | |
Liu et al. | Field dynamic balancing for rigid rotor-AMB system in a magnetically suspended flywheel | |
CN108897218B (zh) | 一种基于混合奇次重复控制器的磁悬浮转子奇次谐波振动抑制方法 | |
CN106647843A (zh) | 一种基于复合分数重复控制器的磁悬浮转子谐波电流抑制方法 | |
CN108279571A (zh) | 一种有限转角机电伺服系统的模型参数辨识方法 | |
CN105783898B (zh) | 一种基于频域自适应lms算法的磁悬浮转子谐波振动抑制方法 | |
US20200217738A1 (en) | Stage-by-stage Measurement, Regulation and Distribution Method for Dynamic Characteristics of Multi-Stage Components of Large-Scale High-Speed Rotary Equipment Based on Multi-Biased Error Synchronous Compensation | |
CN110145541A (zh) | 一种基于相位稳定的磁悬浮轴承转子不平衡运动控制方法 | |
CN110231133B (zh) | 一种磁悬浮轴承电流刚度和位移刚度测量方法 | |
CN104950919B (zh) | 一种设计磁悬浮转子自平衡系统自适应滤波器稳定性参数的方法 | |
CN107992063B (zh) | 基于变参数章动阻尼的变速倾侧动量轮进动控制方法 | |
CN107979305A (zh) | 基于滑模变结构控制的双电机速度同步控制方法及系统 | |
Bu et al. | Adaptive feedforward vibration compensation control strategy of bearingless induction motor | |
Yue et al. | Unbalance identification of speed-variant rotary machinery without phase angle measurement | |
CN108534996A (zh) | 一种基于微振动瀑布图的飞轮刚度参数提取方法 | |
CN108509661A (zh) | 一种轴系动平衡加重影响系数的仿真计算方法和装置 | |
CN109143861A (zh) | 一种基于力矩器的变速倾侧动量轮系统主动振动抑制方法 | |
CN112953344B (zh) | 一种无轴承异步电机转子不平衡振动补偿控制方法 | |
CN109847952B (zh) | 一种基于驱动电流的双轴精密离心机回转台动平衡方法 | |
Zhang et al. | Research on the identification of the moment of inertia of PMSM for industrial robots |
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 |