发明名称 一种基于变块长频域LMS算法的磁悬浮转子谐波电流抑制方法
摘要 本发明公开了一种基于变块长频域LMS(Least Mean Square,LMS)算法的磁悬浮转子谐波电流抑制方法,该方法针对磁悬浮控制力矩陀螺,首先建立含质量不平衡和传感器谐波的磁悬浮转子动力学模型,然后设计了一种基于变块长频域LMS算法的谐波电流抑制方法,根据误差自适应地调整块长,改善收敛性能。本发明能对磁悬浮转子中磁轴承线圈电流的谐波分量进行抑制,适用于存在质量不平衡和传感器谐波的磁悬浮转子系统谐波电流抑制。
申请公布号 CN105159341B 申请公布日期 2017.04.26
申请号 CN201510502598.9 申请日期 2015.08.14
申请人 北京航空航天大学 发明人 崔培玲;崔健;房建成;张大川;李胜
分类号 G05D19/02(2006.01)I 主分类号 G05D19/02(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 杨学明;顾炜
主权项 一种基于变块长频域LMS算法的磁悬浮转子谐波电流抑制方法,其特征在于:包括以下步骤:步骤(1)建立含质量不平衡和传感器谐波的磁悬浮转子动力学模型主被动磁悬浮控制力矩陀螺转子主动磁轴承控制转子径向两自由度实现悬浮,其余三个自由度由安装在转子和定子上的永磁环实现无源稳定悬浮,设转子质心所在的中心平面为Ι,径向磁轴承的定子中心线与面Ι交于点N,转子几何轴和惯性轴分别交面Ι于O、C两点,在Ι内,以N为原点建立惯性坐标系NXY,以O为原点建立旋转坐标系Oεη,设OC长度为l,OC与Oε坐标轴的夹角为θ,O、C在坐标系NXY中的坐标分别为(x,y)、(X,Y);主被动磁轴承轴承力由主动磁轴承电磁力和被动磁轴承磁力两部分组成,X通道中轴承力f<sub>x</sub>可写为:f<sub>x</sub>=f<sub>ex</sub>+f<sub>px</sub>其中,f<sub>ex</sub>为X通道主动磁轴承电磁力,f<sub>px</sub>为X通道被动磁轴承磁力,被动磁轴承产生的力与位移呈线性关系,表示为:f<sub>px</sub>=K<sub>pr</sub>x其中,K<sub>pr</sub>是被动磁轴承位移刚度;当转子悬浮在磁中心附近时,主动磁轴承电磁力可近似线性化为:f<sub>ex</sub>≈K<sub>er</sub>x+K<sub>i</sub>i<sub>x</sub>其中,K<sub>er</sub>、K<sub>i</sub>分别为主动磁轴承位移刚度、电流刚度,i<sub>x</sub>为功放输出电流;对于含有质量不平衡的转子系统,有:X(t)=x(t)+Θ<sub>x</sub>(t)其中,X(t)为转子质心位移,x(t)为转子几何中心位移,Θ<sub>x</sub>(t)为质量不平衡引起的位移扰动,记为:Θ<sub>x</sub>(t)=l cos(Ωt+θ)其中,l为质量不平衡的幅值,θ为相位,Ω为转子转速;在实际转子系统中,由于机械加工精度和材料的不均匀因素的影响,传感器谐波不可避免,传感器实际测得的位移x<sub>s</sub>(t)可表示为:x<sub>s</sub>(t)=x(t)+x<sub>d</sub>(t)其中,x<sub>d</sub>(t)为传感器谐波,可写为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>x</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>c</mi><mi>i</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mi>&Omega;</mi><mi>t</mi><mo>+</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001189703090000021.GIF" wi="498" he="133" /></maths>其中,c<sub>i</sub>是传感器谐波系数的幅值,θ<sub>i</sub>是传感器谐波系数的相位,k为传感器谐波的最高次数;将i<sub>x</sub>、X、Θ<sub>x</sub>、x<sub>d</sub>依次进行拉普拉斯变换得i<sub>x</sub>(s)、X(s)、Θ<sub>x</sub>(s)、x<sub>d</sub>(s),写出转子动力学方程有:ms<sup>2</sup>X(s)=(K<sub>er</sub>+K<sub>pr</sub>)(X(s)‑Θ<sub>x</sub>(s))+K<sub>i</sub>i<sub>x</sub>(s)其中,i<sub>x</sub>(s)=‑K<sub>s</sub>K<sub>i</sub>G<sub>c</sub>(s)G<sub>w</sub>(s)(X(s)‑Θ<sub>x</sub>(s)+x<sub>d</sub>(s))其中,K<sub>s</sub>为位移传感器环节、G<sub>c</sub>(s)为控制器环节,G<sub>w</sub>(s)为功放环节;从上式可以看出,由于质量不平衡和传感器谐波的存在,不仅导致线圈电流中存在与转速同频的基频成分,还存在倍频成分,线圈电流中的同频和倍频成分合称为谐波电流,谐波电流不仅会增加磁轴承功耗,还会引起谐波振动,并通过磁轴承直接传递给航天器,影响航天器姿态控制的精度,需要加以抑制;步骤(2)设计基于变块长频域LMS算法的谐波电流抑制方法以谐波电流为控制目标,将谐波电流i<sub>x</sub>输入至谐波抑制模块,谐波抑制模块的输出反馈至原控制系统的功放输入端,谐波电流抑制模块的计算过程如下:为了抑制谐波电流,以与谐波电流相关的正弦信号作为参考输入,即由和谐波电流频率相同的单位幅值的正弦信号相加组成,电流作为基本输入,根据基本输入与参考输入,不断自适应地调整权值向量,使谐波抑制模块的输出逐渐趋近于控制信号中的谐波误差,最终使控制信号中的谐波成分趋近于零,进而使功放输出的谐波电流成分趋于零,频域LMS(Least Mean Square,LMS)算法在频域完成滤波器权值向量的自适应,从而减小了计算量,此外,块长是调和高收敛速度和小稳态误差的重要参数,为了提高收敛速度减小稳态误差,采用变块长频域LMS算法,根据误差变化情况设计块长更新算法,从而改善频域LMS算法的收敛性能。
地址 100191 北京市海淀区学院路37号