Principle analysis and modeling of rotor imbalance vibration in magnetically suspended control and sensing gyroscope
-
摘要:
磁悬浮控制敏感陀螺(MSCSG)是一种新概念陀螺,采用洛伦兹力磁轴承为力矩器驱动转子径向偏转。针对MSCSG转子旋转过程中产生不平衡振动的问题,分析了不平衡振动产生原理,并建立了解析模型。首先,分析了MSCSG的工作原理。然后,确定了转子不平衡条件下转子几何轴与惯性轴间的几何解析关系;推导了转子不平衡振动力矩数学模型,并对不平衡扰动量的能观性进行了判定;建立了包含振动源的磁轴承-转子控制系统模型,对闭环系统的不平衡振动产生机理进行了分析,并对不同转速下不平衡振动的响应特性进行仿真,仿真结果验证了所提出模型的正确性。最后,根据转子不平衡振动的特点提出了对其进行抑制的要求,为实现MSCSG转子不平衡振动控制奠定了理论基础。
-
关键词:
- 磁悬浮控制敏感陀螺(MSCSG) /
- 洛伦兹力磁轴承 /
- 转子不平衡 /
- 能观性 /
- 振动建模
Abstract:Magnetically suspended control and sensing gyroscope (MSCSG) is a kind of new-concept gyro, which takes Lorentz force magnetic bearing as torquer to drive the rotor to tilt in radial direction. As there is dynamic unbalance in the magnetically suspended rotor system because of the uneven mass distribution, the generation principle of imbalance vibration is analyzed and the analytic model of it is established. First, the working principle of MSCSG is introduced. Then, the geometric analytic relation between geometric and inertial axis of rotor is determined on condition that rotor is unbalanced; the mathematic model of unba-lance vibration torque is established and the observability of imbalance disturbance is demonstrated. The model of bearing-rotor control system containing vibration source is constructed and the vibration generation mechanism in closed-loop system is analyzed. The dynamic response characteristics of unbalance vibration with different rotate speeds are simulated and the simulation result indicates the correctness of the proposed model. Finally, the requirement for suppression of unbalance vibration is put forward according to its vibration characteristics, which lays the theoretical foundation for realizing MSCSG rotor imbalance vibration control.
-
磁悬浮控制敏感陀螺(Magnetically Suspended Control and Sensing Gyroscope, MSCSG)是一种新概念陀螺[1-2],其融合了磁悬浮控制力矩陀螺(Magnetically Suspended Control Momentum Gyroscope, MSCMG)对载体输出控制力矩[3-4]以及转子式速率陀螺仪对载体进行姿态测量[5-6]的功能,将航天器姿态控制系统的执行机构与敏感器合二为一,不仅继承了MSCMG转子系统无接触、无摩擦、寿命长、精度高[7]的优点,还明显降低了航天器姿态控制系统的体积、质量、功耗和成本。MSCSG转子系统采用五自由度全主动控制,具有一定小角度的微框架效应,能够瞬间输出较大的陀螺力矩,实现姿态控制;载体姿态变化时,磁轴承对转子沿径向施加二自由度控制力矩,通过测量偏转磁轴承的控制力矩间接实现姿态敏感。然而,受转子质量分布不平衡因素的影响,磁轴承-转子系统产生的扰动力矩传递到载体上,影响力矩输出精度和姿态测量精度。因此,对MSCSG转子系统进行不平衡振动控制是实现其测控一体化的必然要求,而对不平衡振动机理进行分析及动力学建模是实现不平衡振动控制的前提条件。
转子不平衡是磁悬浮转子系统产生振动的最主要因素[8]。目前,国内外学者对各类磁悬浮转子不平衡振动的原理及建模研究取得了一系列成果。文献[9-12]分析了质量分布不平衡条件下杆状磁悬浮电机转子的振动机理,并在此基础上建立了磁轴承-转子系统动力学模型;文献[13-14]基于模态平衡理论分析了挠性转速下杆状磁悬浮电机转子的运动特性,推导了柔性转子不平衡振动模型;文献[15-17]针对基于磁阻力磁轴承支承的MSCMG转子系统存在不平衡扰动以及转子检测轴与旋转轴不重合2种不同角度,分别建立转子系统动力学模型;文献[18-20]分析了基于混合磁轴承支承的MSCMG转子系统不平衡振动机理,针对混合磁轴承对转子扭转自由度不存在主动可控电流的特殊性,在分析主被动通道磁力非线性的基础上,建立了混合磁轴承支承的磁悬浮转子动力学模型;文献[21-22]在建立包含动不平衡的磁悬浮转子系统动力学模型基础上,分析了模型中不平衡量对转子控制系统的影响。
由于MSCSG是一种新概念陀螺,采用双球形转子结构,通过具有线性特性的洛伦兹力磁轴承(Lorentz Force Magnetic Bearing,LFMB)支承,其支承原理与上述磁悬浮电机、MSCMG有较大区别,因此上述不平衡振动分析方法无法直接应用于MSCSG转子系统。而对这类新型陀螺转子的不平衡振动机理尚未开展系统的理论研究,关于其扰动量能观性的解析性分析也未见报道。因此,本文在对MSCSG转子系统工作机理进行分析的基础上,对转子不平衡问题进行了几何描述,推导了高速转子质量不平衡产生的扰动力矩数学模型,并在此基础上对扰动量的可观测性进行了解析性分析;建立了不平衡振动条件下磁轴承-转子控制系统模型,分析了闭环控制系统不平衡振动的产生机理,并对不平衡振动的响应特性进行仿真,验证了所建立模型的正确性;根据不平衡振动的特点提出了对其进行主动振动控制的要求。
1. MSCSG工作原理
MSCSG结构如图 1所示,主要由陀螺房、陀螺转子、轴向磁轴承、力矩器、径向磁轴承、旋转电机、位移传感器构成。其中,旋转电机驱动转子绕轴向高速旋转,力矩器驱动转子绕径向偏转,径向磁轴承驱动转子沿径向平动,轴向磁轴承驱动转子沿轴向平动。
图 1中,MSCSG的转子组件为双球形包络面结构,由位于转子轴向和径向的2个共球心球缺组成,其中轴向球缺构成轴向磁轴承的转子部分,径向球缺构成径向磁轴承的转子部分。由于径向磁轴承及轴向磁轴承的磁极都呈球状,球面磁极产生的电磁力始终经过磁极球心,无论转子在磁间隙内处于什么位置,都不会产生扭转力矩,从而实现了平动自由度对径向转动自由度的解耦。因此,转子径向偏转只受力矩器控制。
MSCSG采用LFMB为力矩器驱动转子偏转,LFMB产生的电磁力依据安培力定律,即磁感应强度为B的磁场中沿与磁场垂直方向放置长度为L的线圈,当流经线圈中电流为I时,线圈将受到大小为BIL的安培力作用。LFMB组件结构如图 2所示。转子外沿一周的狭长内壁上放置着上、下2层磁钢,2层磁钢的内、外磁钢间充磁方向相反,因此形成了如图 2中虚线所示的闭合磁场。LFMB的定子部分由4组匝数相同的线圈构成,位于内、外磁钢间的狭缝中,沿LFMB周向均匀分布,4组线圈成对使用,正对的2组线圈为一对,用于实现转子的二自由度偏转控制。
以LFMB几何中心O为原点定义定子坐标系O-XYZ,其中X轴与Y轴分别与相对方向2组线圈中心线重合,Z轴方向根据右手定则确定。当与磁场垂直方向放置的线圈通入电流时,线圈的上下两部分将分别产生垂直于线圈及磁场方向的安培力,合力大小为
(1) 式中:n为线圈匝数。
根据如图 2所示的LFMB工作原理图,LFMB提供的X、Y方向偏转力矩分别为
(2) 式中:iX+、iX-分别为X轴正、负方向上的线圈驱动电流;iY+、iY-分别为Y轴正、负方向上线圈驱动电流;lm为LFMB定子半径。令iY+=iα,iX-=iβ,当相对方向线圈通入大小相等、方向相同的电流时,线圈将产生大小相等、方向相反的安培力,形成力偶驱动转子径向偏转,此时,根据式(2)可知偏转力矩pX、pY表达式分别为
(3) 2. MSCSG转子不平衡问题几何分析
MSCSG转子除轴向旋转自由度由电机驱动控制,其余5个由磁轴承悬浮的自由度均需通过传感器来获得转子位置信息。该位置信息既提供给闭环控制器进行反馈控制,也作为转子悬浮工作状态的判定依据。转子偏转自由度为LFMB控制的自由度,需要位移传感器测量转子的位置信息来实现偏转通道的偏转角信息反馈。MSCSG偏转通道4个位移传感器安装在上陀螺房的传感器支架上,配置方式如图 3(a)所示,4个传感器在同一平面上,过定子坐标系原点O对检测面投影,如图 3(b)中黑色虚线所示。
以定子坐标系原点O为原点定义几何坐标系O-xyz,其中x轴与y轴分别与相对方向2组线圈中心线重合,z轴方向根据右手定则确定。由于线圈对向中心线与几何坐标系的x轴、y轴重合,因此定子坐标系与几何坐标系重合。hy+、hy-、hx+、hx-分别为投影在y轴正、负方向以及x轴正、负方向上的轴向传感器测得的转子轴向位移。转子偏转角可以通过同一条直线上2个传感器位置差获得。以转子绕x轴偏转的α通道为例,如图 4所示,转子的偏转角可以通过转子在该通道两端的位置差与对应传感器间跨距比值获取。计算公式为
(4) 同理,转子绕y轴偏转的β通道偏转角计算公式为
(5) 定义转子惯性轴为高速转子旋转自由度所对应的惯性主轴。转子质量分布均匀的情况下,惯性轴与几何轴相一致。然而,转子材料不均匀和加工、装配等误差会使转子质量分布不平衡,导致转子的惯性轴与几何轴不再一致,从而使转子惯性轴相对于几何轴产生如图 5(a)所示的偏转。质量分布不平衡条件下,MSCSG转子惯性轴与几何轴的角位置关系如图 5(b)所示,(αg, βg)为几何轴的角坐标,(αi, βi)为惯性轴的角坐标,g(ξ, ζ)为与转子固连的旋转角坐标系,转子惯性主轴在g(ξ, ζ)中的幅值和相位分别为δ和φ。
由图 5中的转子惯性轴与几何轴间关系可知,惯性坐标系下转子偏转角信息αi、βi表达式为
(6) 式中,α、β分别由式(4)和式(5)得到;Ω为转子轴向转速;Δα、Δβ分别为径向两偏转角在惯性坐标系与几何坐标系下的偏差。
3. 磁轴承-转子控制系统不平衡振动建模
3.1 不平衡扰动项能观性分析
根据欧拉定理写出陀螺技术方程[23]:
(7) 式中:Jx、Jy分别为转子在x、y方向的转动惯量,且Jx、Jy数值上与转子径向转动惯量Jr相等;Jz为转子在z方向的转动惯量。
由于控制系统中各测量值都是在几何坐标系中获取,而陀螺技术方程是在惯性坐标系下定义。因此,将式(6)中转子惯性坐标系与几何坐标系间的关系表达式代入式(7),得到惯性坐标系下转子信息表达式为
(8) 根据式(8)可知,不平衡扰动力矩表达式分别为
(9) 式中:
(10) 根据以上分析可知,当转子质量分布不均匀时,转子的几何轴与惯性轴不再一致,此时转子系统将产生不平衡扰动力矩,作用于与定子固连的基座,引起载体航天器的振动。
将式(3)、式(9)代入式(8),得
(11) 令
,u=[iα iβ]T,可得到式(8)对应的状态方程表达式为(12) 式中:
(13) (14) (15) 对扰动项d列写状态方程,取状态变量T=[Tx Ty]T,则扰动项pd=[pdx pdy]T可用如下模型表述:
(16) 对于式(16)表述的扰动模型,其能观性矩阵N表达式为
(17) 从式(17)表达式可知,转速Ω不等于0时,能观性矩阵的秩rank(N)=2,式(16)能观,因此系统的不平衡扰动可测。
3.2 含振动源的磁轴承-转子控制系统动力学建模
转子不平衡条件下,MSCSG转子偏转控制系统闭环结构如图 6所示。控制通道主要包括传感器、抗混叠滤波器、控制器和功率放大器构成。闭合回路中,传感器检测到转子偏转信息并通过抗混叠滤波器进行降噪处理,控制器根据偏转信息确定偏转指令,随后功率放大器将偏转指令转化为控制电流,驱动LFMB产生控制力矩,使转子偏转至指定参考位置。
由式(7)可知,转子偏转通道间存在耦合,因此为抑制陀螺耦合效应,控制器采用基于交叉滤波的PID控制算法,其中的PID控制器和交叉滤波器传递函数表达式分别为
(18) (19) 式中:kp、ki、kd分别为比例、积分、微分系数;ωD为不完全微分系数;kl和kh分别为低通和高通滤波器增益;a1l、a2l为低通滤波器系数;a1h、a2h为高通滤波器系数。
由如图 4所示的控制系统框图可知,根据转子偏转角可确定控制电流为
(20) 式中:ks为传感器比例增益;gw、ga、gPID和gcr分别为抗混叠滤波器、功率放大器、PID控制器和交叉滤波器的输入-输出传递函数变换算子,对其进行拉普拉斯变换,满足关系式
[gw(d/dt)]=gw(s), [ga(d/dt)]=ga(s), [gPID(d/dt)]=gPID(s), [gcr(d/dt)]=gcr(s)。将式(3)、式(6)、式(20)代入式(7),得到MSCSG转子偏转动力学方程组:
(21) 将式(21)进行拉普拉斯变换得
(22) 对式(22)进行化简,得到
(23) 式(23)的等效闭环控制系统框图如图 7所示,αr、βr为偏转指令信号。对图 7中的转子不平衡条件下MSCSG转子系统模型进行分析可知,转子不平衡量经过磁轴承控制器和功率放大器进入闭环控制系统,通过LFMB产生电流刚度力矩,引起转子产生不平衡振动。
4. 仿真分析
根据本文所建立的MSCSG磁轴承-转子控制系统不平衡振动模型,对惯性轴与几何轴不一致条件下磁轴承-转子控制系统产生的不平衡振动力矩进行仿真分析。仿真参数如表 1所示,其中,转子参数根据MSCSG设计指标确定,控制器相关参数依照文献[23]设定。表 1中:ka为功率放大器增益;ωa为功率放大器截止频率;ωf为抗混叠滤波器截止频率。
表 1 MSCSG系统参数Table 1. Parameters of MSCSG system参数 数值 Jz/(kg·m2) 0.016 6 Jy/(kg·m2) 0.009 7 ki 0.001 ωD 240 kh 0.001 5 n 200 l/m 0.115 8 a2l 2 200 000 a2h 92 100 ls/m 0.078 ka/(V·A-1) 0.22 ωf/Hz 310 Jx/(kg·m2) 0.009 7 kp 15.1 kd 3.2 kl 0.001 5 B/T 0.4 lm/m 0.059 a1l 2 400 a1h 370 ks/(V·m-1) 10 300 δ/(°) 0.009 ωa/Hz 240 φ/(°) 10 转子转速由0升至6 000 r/min时,不平衡振动力矩的响应仿真结果如图 8所示。可知,瀑布图中包含了明显的转速同频成分,不平衡振动频率随转子转速升高而增加,不平衡振动幅值与转子转速呈递增的趋势。因此,不平衡振动仿真结果与动力学模型所描述不平衡因素的频率特性相一致,验证了本文所建立模型的正确性。根据以上分析可知,要想实现MSCSG转子系统的主动振动控制,必须对转子不平衡引起的转速同频扰动量进行抑制。
5. 结论
1) MSCSG转子质量分布不平衡条件下,转子惯性主轴与几何主轴不一致,导致转子系统产生不平衡振动力矩。
2) MSCSG转子转速不为零的条件下,转子质量分布不平衡产生的振动力矩可观测。
3) MSCSG转子不平衡量经过磁轴承-转子控制系统产生电流刚度力矩,引起的振动与转子转速同频,为实现对转子的主动振动控制,必须对转速同频扰动量进行抑制。
-
表 1 MSCSG系统参数
Table 1. Parameters of MSCSG system
参数 数值 Jz/(kg·m2) 0.016 6 Jy/(kg·m2) 0.009 7 ki 0.001 ωD 240 kh 0.001 5 n 200 l/m 0.115 8 a2l 2 200 000 a2h 92 100 ls/m 0.078 ka/(V·A-1) 0.22 ωf/Hz 310 Jx/(kg·m2) 0.009 7 kp 15.1 kd 3.2 kl 0.001 5 B/T 0.4 lm/m 0.059 a1l 2 400 a1h 370 ks/(V·m-1) 10 300 δ/(°) 0.009 ωa/Hz 240 φ/(°) 10 -
[1] 夏长峰, 蔡远文, 任元, 等.磁悬浮控制敏感陀螺转子前馈解耦内模控制[J].北京航空航天大学学报, 2018, 44(3):480-488.XIA C F, CAI Y W, REN Y, et al.Feedforward decoupling and internal model control for rotor of magnetically suspended control and sensing gyroscope[J].Journal of Beijing University of Aeronautics and Astronautics, 2018, 44(3):480-488(in Chinese). [2] 夏长峰, 蔡远文, 任元, 等.MSCSG转子系统的扩展双频Bode图稳定性分析方法[J].宇航学报, 2018, 39(2):168-176.XIA C F, CAI Y W, REN Y, et al.Stability analysis method with extended double-frequency Bode diagram for rotor of MSCSG[J].Journal of Astronautics, 2018, 39(2):168-176(in Chinese). [3] ZHENG S Q, LI H T, HAN B C, et al.Power consumption reduction for magnetic bearing systems during torque output of control moment gyros[J].IEEE Transactions on Power Electronics, 2017, 32(7):5752-5759. doi: 10.1109/TPEL.2016.2608660 [4] HAN B C, ZHENG S Q, LI H T, et al.Weight-reduction design based on integrated radial-axial magnetic bearing of a large scale MSCMG for space station application[J].IEEE Transactions on Industrial Electronics, 2016, 64(3):2205-2214. [5] LIU X K, ZHAO H, YAO Y, et al.Modeling and analysis of micro-spacecraft attitude sensing with gyrowheel[J].Sensors, 2016, 16(8):1321-1328. doi: 10.3390/s16081321 [6] FANG J C, ZHENG S Q, HAN B C.Attitude sensing and dynamic decoupling based on active magnetic bearing of MSDGCMG[J].IEEE Transactions on Instrumentation and Measurement, 2012, 61(2):338-348. doi: 10.1109/TIM.2011.2164289 [7] REN Y, CHEN X C, CAI Y W, et al.Attitude-rate measurement and control integration using magnetically suspended control and sensitive gyroscopes[J].IEEE Transactions on Industrial Electronics, 2018, 65(6):4921-4932. doi: 10.1109/TIE.2017.2772161 [8] ZHENG S Q, HAN B C, FENG R, et al.Vibration suppression control for AMB-supported motor driveline system using synchronous rotating frame transformation[J].IEEE Transactions on Industrial Electronics, 2015, 62(9):5700-5708. doi: 10.1109/TIE.2015.2407857 [9] CHEN Q, LIU G, ZHENG S Q.Suppression of imbalance vibration for AMBs controlled driveline system using double-loop structure[J].Journal of Sound & Vibration, 2015, 337(2):1-13. [10] ZHENG S Q, FENG R.Feedforward compensation control of rotor imbalance for high-speed magnetically suspended centrifugal compressors using a novel adaptive notch filter[J].Journal of Sound & Vibration, 2016, 366(3):1-14. [11] ZHOU J X, ZHENG S Q, HAN B C, et al.Effects of notch filters on imbalance rejection with heteropolar and homopolar magnetic bearings in a 30-kW 60000-rpm motor[J].IEEE Transactions on Industrial Electronics, 2017, 64(10):8033-8039. doi: 10.1109/TIE.2017.2694412 [12] CHEN Q, LIU G, HAN B C.Unbalance vibration suppression for AMBs system using adaptive notch filter[J].Mechanical Systems and Signal Processing, 2017, 93(9):136-150. [13] WANG Y G, FANG J C, ZHENG S Q.A field balancing technique based on virtual trial-weights method for a magnetically levitated flexible rotor[J].Journal of Engineering for Gas Turbines & Power, 2014, 136(9):1-7. [14] TANG E Q, HAN B C, ZHANG Y.Optimum compensator design for the flexible rotor in magnetically suspended motor to pass the first bending critical speed[J].IEEE Transactions on Industrial Electronics, 2016, 63(1):343-354. doi: 10.1109/TIE.2015.2472534 [15] FANG J C, XU X B, TANG J Q, et al.Adaptive complete suppression of imbalance vibration in AMB systems using gain phase modifier[J].Journal of Sound & Vibration, 2013, 332(24):6203-6215. [16] XU X B, FANG J C, LI H T, et al.Active suppression of imbalance vibration in the magnetically suspended control moment gyro[J].Journal of Vibration and Control, 2015, 21(5):990-998. [17] PENG C, FAN Y H, HUANG Z Y, et al.Frequency-varying synchronous micro-vibration suppression for a MSFW with application of small-gain theorem[J].Mechanical Systems and Signal Processing, 2016, 82(6):432-447. [18] CUI P L, WANG Q R, ZHANG G X, et al.Hybrid fractional repetitive control for magnetically suspended rotor systems[J].IEEE Transactions on Industrial Electronics, 2017, 65(4):3491-3495. [19] 崔培玲, 盖玉欢, 李海涛.基于复合前馈补偿的混合磁悬浮CMG转子主动振动控制[J].振动与冲击, 2015, 34(22):162-166.CUI P L, GAI Y H, LI H T.Active vibration control of an active-passive hybrid magnetically suspended rotor based on composite feedforward compensation method[J].Journal of Vibration and Shock, 2015, 34(22):162-166(in Chinese). [20] 刘超, 刘刚, 赵光再.主被动磁悬浮高速转子系统的自动平衡控制[J].光学精密工程, 2015, 23(3):714-722.LIU C, LIU G, ZHAO G Z.Autobalancing control of high-speed rotor suspended by active-passive hybrid magnetic bearings[J].Optics and Precision Engineering, 2015, 23(3):714-722(in Chinese). [21] YU Y J, YANG Z H, HAN B C, et al.Active vibration control of magnetically suspended wheel using active shaft deflection[J].IEEE Transactions on Industrial Electronics, 2017, 64(8):6528-6533. doi: 10.1109/TIE.2017.2682786 [22] HE Y, SHI L, SHI Z G, et al.Unbalance compensation of a full scale test rig designed for HTR-10GT:A frequency-domain approach based on iterative learning control[J].Science and Technology and Nuclear Installations, 2017, 17(1):2-7. [23] 房建成, 任元.磁悬浮控制力矩陀螺技术[M].北京:国防工业出版社, 2014:107-120.FANG J C, REN Y.Magnetically suspended control momentum gyroscope technology[M].Beijing:National Defense Industry Press, 2014:107-120(in Chinese). -