首页> 中国专利> 一种基于非线性自适应算法的磁轴承系统惯性轴辨识方法

一种基于非线性自适应算法的磁轴承系统惯性轴辨识方法

摘要

本发明涉及一种基于非线性自适应算法的磁轴承系统惯性轴辨识方法。首先建立包含转子不平衡和Sensor Runout的磁轴承系统转子动力学模型;其次提出一种非线性自适应控制律和估计律,保证磁悬浮转子惯性中心位移估计值趋于零的同时,又能估计出Sensor Runout高次谐波分量的傅里叶系数;然后通过改变转子转速的策略,增加系统的可观测度,实现Sensor Runout同频分量和转子不平衡量的辨识,即实现惯性轴的辨识;最后修正自适应算法中同频分量傅里叶系数,准确地抑制多谐波电流和补偿位移刚度力,实现磁悬浮惯性执行机构的多谐波振动抑制。

著录项

  • 公开/公告号CN106289208A

    专利类型发明专利

  • 公开/公告日2017-01-04

    原文格式PDF

  • 申请/专利权人 北京航空航天大学;

    申请/专利号CN201610854370.0

  • 申请日2016-09-27

  • 分类号G01C19/24(20060101);G01C25/00(20060101);

  • 代理机构11251 北京科迪生专利代理有限责任公司;

  • 代理人成金玉;卢纪

  • 地址 100191 北京市海淀区学院路37号

  • 入库时间 2023-06-19 01:17:24

法律信息

  • 法律状态公告日

    法律状态信息

    法律状态

  • 2018-01-19

    授权

    授权

  • 2017-02-15

    实质审查的生效 IPC(主分类):G01C19/24 申请日:20160927

    实质审查的生效

  • 2017-01-04

    公开

    公开

说明书

技术领域

本发明涉及一种基于非线性自适应算法的磁轴承系统惯性轴辨识方法,用于辨识包含转子不平衡和传感器谐波噪声(Sensor Runout)的磁轴承系统惯性轴,实现磁悬浮惯性执行机构的多谐波振动抑制,使磁悬浮惯性执行机构满足未来“超静超稳”卫星平台对极微振动的要求,属于磁轴承系统主动振动控制领域。

背景技术

随着高分辨率对地观测、深空探测、星间激光通信等超高分辨率卫星的发展,超静与敏捷成为衡量卫星平台性能的两项重要指标。超静性能是保证高分辨率有效载荷成像质量的关键因素,敏捷性能是解决高分辨率成像与大范围覆盖这对矛盾的一条技术捷径。越来越高的分辨率指标对卫星平台的指向精度和姿态稳定度的要求越来越高,对星载活动部件所引起的振动越来越敏感。卫星振动主要分为两大类,一类是几到几十Hz的低频高幅振动,这类振动可以通过卫星姿态控制进行抑制;另一类高频低幅振动主要由飞轮、控制力矩陀螺等惯性执行机构引起,这类振动是无法通过姿态控制算法进行抑制的,是影响卫星平台振动水平的技术瓶颈。

对惯性执行机构振动的抑制主要有隔振装置和磁悬浮主动振动控制两种。机械式惯性执行机构通常采用隔振技术来抑制高频振动,但是隔振装置只是将高频低幅振动转化成低频高幅振动,没有从源头上消除振动。磁悬浮惯性执行机构一个重要优点是具有主动振动抑制的能力,其本质是通过调节磁轴承的控制力实现转子绕惯性主轴旋转,从根本上消除高速转子的高频振动。由于加工安装误差、材质不均匀、电子元器件非线性等机械与电气非理想特性,磁悬浮惯性执行机构存在着转子不平衡、Sensor Runout等振动源,从而传递出多谐波振动。目前,磁轴承系统主动振动控制主要集中于转子不平衡振动控制的研究,对包含转子不平衡和Sensor Runout的磁轴承系统主动振动控制研究较少。由于位移传感器输出的同频成分包含Sensor Runout同频分量和转子不平衡,在进行位移刚度补偿时只需要补偿转子不平衡的影响,因此需要进行Sensor Runout同频分量和转子不平衡的辨识,即惯性轴的辨识。现有的辨识方法可分为两类:一类是直接辨识法,直接将磁悬浮转子低速运行时传感器输出的同频成分作为Sensor Runout的同频分量,这类方法辨识误差大;另一类是自适应辨识算法,但现有算法主要实现转子绕几何轴旋转,会对外传递出相当大的振动,只适用于高精度对准磁轴承系统,不适用于磁悬浮惯性执行机构。

发明内容

本发明的技术解决问题是:克服现有技术的不足,发明一种基于非线性自适应算法的磁轴承系统惯性轴辨识算法,通过非线性自适应算法和改变转子转速的策略,提高惯性轴辨识精度,实现磁悬浮转子绕惯性轴旋转,抑制磁轴承系统多谐波振动。此外,本发明中的非线性自适应控制律解决了利用传统线性算法时初始电流过大和参数收敛速度慢等问题。

本发明的技术解决方案是:一种基于非线性自适应算法的磁轴承系统惯性轴辨识算法,首先建立包含转子不平衡和Sensor Runout的磁轴承系统动力学模型,分析振动产生机理及存在形式;其次发明一种非线性自适应控制律和估计律,在保证磁悬浮转子惯性中心位移的估计值收敛于零的同时,实现转子不平衡和Sensor Runout各谐波分量傅里叶系数的估计;然后通过改变转速的策略,增加同频分量的可观测度,实现Sensor Runout同频分量和转子不平衡量的辨识,估计出惯性轴;最后实现磁悬浮惯性执行机构的多谐波主动振动抑制。本发明的具体步骤如下:

(1)建立含转子不平衡和Sensor Runout的磁轴承系统动力学模型

对于两自由度磁轴承系统,x轴和y轴两通道相互解耦。假设x轴和y轴的位移刚度系数和电流刚度系数相同,当磁悬浮转子在平衡位置附近运动时,其线性化的动力学方程为:

式中,m为磁悬浮转子的质量;ki和kh分别为磁轴承系统的电流刚度系数和位移刚度系数;Ic=[icx,icy]T,icx和icy分别为x轴和y轴磁轴承线圈控制电流;χI=[xI,yI]T,xI和yI分别为惯性中心在x轴和y轴方向上的位移,为χI的二阶导数;χg=[xg,yg]T,xg和yg分别为几何中心在x轴和y轴方向上的位移。

由于转子不平衡的影响,使得转子惯性中心与几何中心不重合,则磁悬浮转子几何中心和惯性中心位移之间的关系为:

χI=χg

式中,为转子不平衡量;λ和分别为转子不平衡量的幅值和相位;ω为转子转速。将δ改为矩阵形式有:

δ=PδΦδT

Pδ=-sin(ωt)cos(ωt)cos(ωt)sin(ωt)

其中,Pδ和Φδ分别为转子不平衡量的三角函数矩阵和傅里叶系数。

此外,受位移传感器多谐波噪声Sensor Runout的影响,传感器输出的几何中心位移χs与实际几何中心位移χg存在偏差,两者之间的关系为:

χs=χg+d

式中,为Sensor Runout向量;σi和ξi分别为Sensor>

d=PdΦdT

Pdi=-sin(iωt)cos(iωt)cos(iωt)sin(iωt),i=1,...,k

Φd=[Φd1>dk]

Φdi=[pi>i],pi=σisinξi,qi=σicosξi

其中,Pd和Φd分别为Sensor>di和Φdi分别为Sensor>

磁轴承线圈控制电流Ic为:

Ic=-kadksGw(s)Gc(s)χs

式中,kad为AD采样系数;ks为位移传感器放大倍数;Gc(s)和Gw(s)分别为控制器和功率放大器的传递函数。

则包含转子不平衡和Sensor Runout的磁轴承系统动力学模型为:

mχ··I=-kikadksGw(s)Gc(s)(χI+δ+d)+kh(χI+δ)

(2)非线性自适应控制算法设计

在步骤(1)所述的模型基础上,设计非线性自适应控制算法,该算法主要包括两部分:自适应控制律和估计律。自适应估计律将磁悬浮转子惯性中心位移的估计值作为控制变量,保证收敛于零。自适应估计律能够自适应地估计转子不平衡量和Sensor Runout各谐波分量的傅里叶系数Φδ、Φd,保证各个傅里叶系数估计值收敛。非线性自适应控制算法设计为:

uic=-(khχ^I+mΞχ^·I+ρe+khPδΦ^δT)/(Kamki)

式中,Ξ为正定矩阵;ρ为正常数;为转子不平衡量傅里叶系数的估计值;设计目的是为了补偿转子不平衡量引起的位移刚度力;Kam为功放系统的等效比例系数;e是与和有关的权值函数;为的一阶导数。

Sensor Runout各谐波分量和转子不平衡量傅里叶系数的自适应估计律分别为:

Φ~·d=eTPrdΓd,Φ~·δ=eTPrδΓδ

式中,和分别为三角函数矩阵Pd和Pδ的二阶导数;和分别为Sensor>

(3)磁轴承系统惯性轴辨识

步骤(2)提出的非线性自适应算法能准确地估计出Sensor Runout高次谐波分量的傅里叶系数,以及保证Sensor Runout同频分量和转子不平衡量的傅里叶系数估计值收敛。为了进一步辨识Sensor Runout同频分量和转子不平衡量,需要通过步骤(3)变转速方式实现,即实现惯性轴辨识。惯性轴辨识主要包括三个步骤:a)转子工作在转速ω1下,得到Sensor>2下,得到当前转速下的估计值和c)根据两个转速条件下得到的同频分量傅里叶系数有:

(p1-p^11)+η1(u-u^1)=0(p1-p^12)+η2(u-u^2)=0(q1-q^11)+η1(v-v^1)=0(q1-q^12)+η2(v-v^2)=0p1u=A-1B1q1u=A-1B2

式中,和

最后求解出真实值p1、q1、u和v,则磁轴承系统惯性轴得到辨识。

(4)磁轴承系统多谐波主动振动抑制

将步骤(2)非线性自适应控制律中的Sensor Runout同频分量和转子不平衡傅里叶系数的估计值替换为步骤(3)计算出的真实值,则由转子不平衡引起的位移刚度力得到精确地补偿,转子不平衡量和Sensor Runout引起的多谐波电流得到有效抑制,最终精确地抑制了磁悬浮惯性执行机构的多谐波振动。

本发明的原理是:转子不平衡和Sensor Runout是磁轴承系统的两个主要振动源,两者产生振动的途径以及形式是各不相同的。转子不平衡不仅通过磁轴承系统本身产生位移刚度力,还通过控制器和电流刚度系数产生电流刚度力;而Sensor Runout只产生电流刚度力。因此磁轴承系统多谐波振动的抑制不仅要实现多谐波电流的抑制,即在位移传感器的输出端补偿两种噪声的影响,而且要补偿转子不平衡引起的位移刚度力。然而位移传感器输出的同频分量中既包含转子不平衡分量,又包含Sensor Runout的同频分量。因此在进行位移刚度力精确补偿的前提,要进行转子不平衡和Sensor Runout同频量的辨识,即惯性轴的辨识。从而,实现磁悬浮惯性执行机构高精度多谐波振动抑制。

本发明与现有技术相比的优点在于:一种基于非线性自适应算法的磁轴承系统惯性轴辨识方法,(1)克服了传统直接辨识算法引起辨识误差大的缺点,仅需要一次升速或者降速可实现磁轴承系统惯性轴的辨识;(2)克服了传统线性自适应算法引起的初始控制电流过大和参数辨识收敛速度慢等缺点,利用非线性自适应控制律改善算法性能;(3)避免了传统磁轴承系统多谐波振动控制算法中,Sensor Runout同频分量和转子不平衡量混叠引起的位移刚度力补偿误差,利用辨识后的傅里叶系数进行修正,从而实现磁轴承系统高精度多谐波振动抑制。

附图说明

图1为本发明的实现流程图;

图2为基于非线性自适应控制算法的系统原理框图;

图3为惯性轴辨识流程图。

具体实施方式

下面结合附图以及具体的实施步骤对本发明做进一步说明。

如图1所示,本发明涉及一种基于非线性自适应算法的磁轴承系统惯性轴辨识方法,其实现过程是:首先建立包含转子不平衡和Sensor Runout的磁轴承系统动力学模型,分析磁轴承系统多谐波振动产生原因及存在形式;其次设计非线性自适应控制算法,实现磁悬浮转子的惯性中心位移估计值收敛于零,并且自适应估计转子不平衡量和Sensor Runout各谐波分量的傅里叶系数;然后通过改变转子转速的策略实现Sensor Runout同频分量和转子不平衡的辨识,即实现惯性轴的辨识;最后将非线性自适应算法中Sensor Runout各谐波分量和转子不平衡傅里叶系数进行修正,使磁悬浮转子绕真正惯性轴旋转,从而实现磁轴承系统多谐波振动抑制。图2为非线性自适应控制算法的原理框图。位移传感器检测出转子的位移,并引入Sensor Runout多谐波噪声d,经过AD采样进入控制器,实现闭环控制。非线性自适应算法中的估计律估计出转子不平衡量和Sensor Runout各谐波分量和在控制器输入端χs将两者消除,以实现多谐波电流的抑制;非线性自适应算法利用估计出的转子不平衡量得到进行位移刚度力的补偿,最终得到控制信号uic。磁轴承功放系统根据控制信号驱动磁轴承线圈,产生控制电流Ic,产生相应的磁轴承力F作用于磁悬浮转子,从而改变转子的位置χg。图3为惯性轴辨识流程图,为图1步骤(3)的具体实施流程,分别在转速ω1和ω2下得到稳态时的估计值和通过求解方程得到真实的傅里叶系数p1,q1,u和v,最后将自适应算法中转子不平衡和Sensor>

(1)建立含转子不平衡和Sensor Runout的磁轴承系统动力学模型

对于两自由度磁轴承系统,x轴和y轴两通道相互解耦。假设x轴和y轴的位移刚度系数和电流刚度系数相同,当磁悬浮转子在平衡位置附近运动时,其线性化的动力学方程为:

{mx··I=kiicx+khxgmy··I=kiicy+khyg---(1)

式中,m为磁悬浮转子的质量;ki和kh分别为磁轴承系统的电流刚度系数和位移刚度系数;icx和icy分别为x轴和y轴磁轴承线圈控制电流;xI和yI分别为惯性轴在x轴和y轴方向上的位移;xg和yg分别为几何轴在x轴和y轴方向上的位移。

将式(1)写成矩阵形式为:

mχ··I=kiIc+khχg---(2)

式中,χI=[xI,yI]T,χg=[xg,yg]T,Ic=[icx,icy]T

由于转子不平衡的影响,使得转子惯性中心与几何中心不重合,则转子惯性中心位移χI和几何中心位移χg之间的关系为:

χI=χg-δ>

式中,δ为转子不平衡,表示为:

其中,λ和分别为转子不平衡量的幅值和相位;ω为转子转速。将式(4)写成矩阵形式为:

δ=PδΦδT---(5)

Pδ=-sin(ωt)cos(ωt)cos(ωt)sin(ωt)---(6)

其中,Pδ和Φδ分别为转子不平衡量的三角函数矩阵和傅里叶系数。

令则转子不平衡量的傅里叶系数向量Φδ可表示为:

Φδ=[u>

此外,受位移传感器多谐波噪声Sensor Runout的影响,传感器输出的几何中心位移χs与实际几何中心位移χg存在偏差,两者之间的关系为:

χs=χg+d>

式中,为Sensor Runout向量;σi和ξi分别为Sensor>

d=PdΦdT---(10)

Pdi=-sin(iωt)cos(iωt)cos(iωt)sin(iωt),i=1,...,k---(12)

Φd=[Φd1>dk]>

Φdi=[pi>i],pi=σisinξi,qi=σicosξi>

其中,Pd和Φd分别为Sensor>di和Φdi分别为Sensor>

磁轴承线圈控制电流Ic为:

Ic=-kadksGw(s)Gc(s)χs>

式中,kad为AD采样系数;ks为位移传感器放大倍数;Gc(s)和Gw(s)分别为控制器和功率放大器的传递函数。

则包含转子不平衡和Sensor Runout的磁轴承系统的动力学模型为:

mχ··I=-kikadksGw(s)Gc(s)(χI+δ+d)+kh(χI+δ)---(16)

则振动力F与两种振动源δ和d之间的关系为:

F=Q-1[(khI2×2-kikadksGw(s)Gc(s))δ-kikadksGw(s)Gc(s)d]>

Q=I2×2-khP(s)+kikadksGw(s)Gc(s)P(s)>

式中,I2×2为二阶单位阵;P(s)为磁轴承系统传递函数。

(2)非线性自适应控制算法设计

假设转子不平衡量和Sensor Runout的估计值分别为和则转子惯性中心的位移估计值可表示为:

χ^I=χ^g-δ^=χs-d^-δ^=χI+d~+δ~---(19)

式中,几何中心的位移估计值;和分别为转子不平衡量和Sensor Runout的估计误差,定义为:

δ~=Δδ-δ^=PδΦ~δT---(20)

d~=Δd-d^=PdΦ~dT---(21)

式中,和分别为转子不平衡量和Sensor Runout各谐波分量傅里叶系数的估计误差。

为了保证转子惯性中心位移估计值收敛于零,且能够自适应地估计出转子不平衡和Sensor Runout傅里叶系数,非线性自适应控制算法设计为:

uic=-(khχ^I+mΞχ^·I+ρe+khPδΦ^δT)/(Kamki)---(22)

式中,ρ为正常数;为了补偿转子不平衡量引起的位移刚度力;为的一阶导数;e为与和有关的权值函数。为了克服传统线性权值函数引起的初始控制电流过大,估计参数收敛速度慢等缺点,本发明提出一种非线性权值函数:

e=χ^·I+a>arctan(bχ^I)---(23)

式中,a和b为正常数。则式(22)中正定矩阵Ξ可表示为:

Ξ=ab1+b2x^I200ab1+b2y^I2---(24)

Sensor Runout各谐波分量和转子不平衡量傅里叶系数自适应估计律分别为:

Φ~·d=eTPrdΓd,Φ~·δ=eTPrδΓδ---(25)

Γd=diag(τd1>d1>d2>d2 …>dk>dk)>

Γδ=diag(τδ>δ)>

式中,和分别为三角函数矩阵Pd和Pδ的二阶导数;和分别为Sensor>dii+Δδii)≤1,Δδii=τδ2

(3)磁轴承系统惯性轴辨识

由于系统渐进稳定,当t→∞时,和都将收敛于零,则磁悬浮转子的惯性中心位移估计值收敛于零。根据式(23)和(25)可知,当t→∞时,则转子不平衡量估计值和Sensor Runout傅里叶系数估计值将收敛于定值。根据式(1)、(19)和(22)有:

md~··+mδ~··=khd~---(28)

由于和都趋向于零,则式(28)简化为:

PrdΦ~dT+PrδΦ~δT=0---(29)

根据三角函数正交特性,可得:

(kh+2)Pd1Φ~d1T+PrδΦ~δT=0---(30)

(kh+m(iω)2)PdiΦ~diT=0,i=2,...,k---(31)

由式(30)和(31)可知:Sensor Runout的高次谐波分量傅里叶系数收敛于真实值,而Sensor Runout的同频分量和转子不平衡的傅里叶系数并未收敛到真实值。因此需要通过一种手段增加系统的可观测度,对两者进行辨识,从而实现磁轴承系统惯性轴辨识。

定义则由式(32)可得:

p~1+ηu~=0q~1+ηv~=0(p1-p^1)+η(u-u^)=0(q1-q^1)+η(v-v^)=0---(32)

两个方程四个未知数,上式方程是不可求解的。由η定义可知,改变转子转速以增加方程个数,实现磁轴承系统惯性轴的辨识。

如图3所示,磁轴承系统惯性轴的辨识需要三个步骤:a)转子工作在转速ω1下,得到Sensor>2下,得到当前转速下的估计值和将两个转速条件下的估计值代入式(32),有:

(p1-p^11)+η1(u-u^1)=0(p1-p^12)+η2(u-u^2)=0(q1-q^11)+η1(v-v^1)=0(q1-q^12)+η2(v-v^2)=0p1u=A-1B1q1u=A-1B2---(33)

式中,和

最后求解出真实值p1、q1、u和v,则磁轴承系统惯性轴得到辨识。

(4)磁轴承系统多谐波振动抑制

将步骤(2)中式(22)的Sensor Runout同频分量和转子不平衡傅里叶系数的估计值替换为步骤(3)计算出的真实值,则由转子不平衡引起的位移刚度力得到精确地补偿,转子不平衡量和Sensor Runout引起的多谐波电流得到有效抑制,最终精确地抑制了磁悬浮惯性执行机构的多谐波振动。

由以上可以看出,利用非线性自适应算法进行惯性轴辨识时,只需要在低速条件下进行一次升速或者降速即可,实现简单,辨识精度高于直接辨识法。此外,由于Sensor Runout各谐波分量不随转子转速变化而变化,因此可以将低速条件下辨识出的傅里叶系数值直接用于高速下的磁悬浮惯性执行机构主动振动控制。

本发明说明书中未作详细描述的内容属于本专业领域技术人员公知的现有技术。

去获取专利,查看全文>

相似文献

  • 专利
  • 中文文献
  • 外文文献
获取专利

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号