冲击激励下轴线失准转子-磁轴承系统不对中定量研究

肖玲; 李园超; 赵晨曦; 程文杰; 冯圣

doi:10.3969/j.issn.0258-2724.20230454

冲击激励下轴线失准转子-磁轴承系统不对中定量研究

doi: 10.3969/j.issn.0258-2724.20230454

1.
西安科技大学理学院，陕西西安 710054
2.
西安交通大学机械工程学院，陕西西安 710049

基金项目: 国家自然科学基金（52275271）；陕西省自然科学基金（2022JM-194）

详细信息

作者简介:
肖玲（1983—），女，教授，博士，研究方向为电磁轴承、高速电机转子、软磁复合材料，E-mail：xiaoling@xust.edu.cn

中图分类号: TH133.3
计量
- 文章访问数: 373
- HTML全文浏览量: 105
- PDF下载量: 92
- 被引次数: 0
出版历程
- 收稿日期: 2023-09-04
- 修回日期: 2024-03-22
- 网络出版日期: 2024-04-29
- 刊出日期: 2024-03-28

Quantitative Research on Misalignment Magnitude of Rotor-Magnetic Bearing System with Axis Misalignment Under Shock Excitation

1.
College of Science, Xi’an University of Science and Technology, Xi’an 710054, China
2.
School of Mechanical Engineering, Xi’an Jiaotong University, Xi’an 710049, China

摘要

摘要:
为研究和识别转子系统在轴承处发生的平行角度混合不对中，提出一种频谱辨识转子-磁轴承系统固有不对中量大小的方法. 采用动量矩定理将圆盘不平衡力对转轴的影响等效到转子轴向力上，建立考虑轴向径向耦合效应的刚性双偏置圆盘转子-磁轴承系统的动力学模型；通过SIMULINK仿真得到系统时域下的位移和电流响应，分析不对中条件下转子系统动力学特性，并利用快速傅里叶变换将时域响应转换为频域响应，基于频域下最小二乘算法得到转子系统不对中量大小. 结果表明：在冲击激励影响条件下，采用该方法计算的不对中量大小误差均在5.0%以内，当转子受到外界扰动力时，该算法能够准确定量识别转子的不对中量，可为不对中转子-磁轴承系统故障诊断及自修复提供理论参考.
- 冲击载荷 /
- 轴承-转子系统 /
- 不对中量 /
- 快速傅里叶变换 /
- 轴向径向耦合
Abstract:
A spectral identification method for calculating the magnitude of inherent misalignment in a rotor-magnetic bearing system was proposed to study and identify mixed parallel misalignment of the rotor system occurring at the bearing. The momentum moment theorem was used to equate the effect of the disc unbalance force on the rotating shaft to the axial force of the rotor and establish a dynamics model of the rigid double offset disc rotor-magnetic bearing system considering the axial and radial coupling effects. The SIMULINK simulation was used to calculate the displacement and current response of the system in the time domain, and the dynamics characteristics of the rotor system under the misalignment condition were analyzed. Furthermore, the fast Fourier transform was utilized to convert the response in the time domain into that in the frequency domain. The magnitude of the misalignment of the rotor system was then calculated based on the least squares algorithm in the frequency domain. The results show that the error in the calculated magnitude of misalignment by using this method is within 5.0% under the influence of shock excitation, indicating that even if the rotor is affected by external disturbance forces, the algorithm can accurately quantify the rotor’s misalignment. This provides a theoretical reference for fault diagnosis and self-repair of misaligned rotor-magnetic bearing systems.
- shock excitation /
- bearing-rotor system /
- misalignment magnitude /
- fast Fourier transform /
- axial and radial coupling

HTML全文

随着现代工业的不断发展，磁轴承因高转速和无摩擦等优点被广泛应用于航空航天、高速机床等领域. 在实际运行中，转子不平衡、弯曲和转子不对中等最为常见，是机器振幅增加的主要来源，严重影响系统的安全稳定运转. 为使高速转子能够平稳高效运行，研究和识别轴承转子系统不对中非常重要，以便于纠正和消除转子轴承系统中不对中故障.

旋转机械转子不对中可分为2种情况，即2个耦合转子的轴线彼此不对中和转子轴线和支撑转子的轴承不对中，这种轴承不对中故障的原因可能是由于安装或组装不当，传感器测量误差以及涡流效应导致膨胀和热变形引起的轴承和转子轴线失准. 陈宏等^[1]提出一种联轴器处不对中定量分析方法，该方法通过安装在转子系统上的电涡流传感器获得转子运动的轴心位置信息，进而得到转子耦合处平行不对中量以及角度不对中量. 甄满等^[2]研究了在弹性联轴器连接下含有不对中故障的转子系统，发现不对中会使系统的非线性特征更加复杂. 当膜盘联轴器径向刚度增大后会出现偶数倍频分量. Shari等^[3]通过多通道振动数据采集系统监测系统振动，并将快速傅里叶变换（FFT）技术和独立分量分析（ICA）相结合来检测旋转机械中发生的不平衡和不对中故障. Kuppa等^[4]提出一种基于最小二乘的识别算法，应用动态降阶方法估计转子-磁轴承-联轴器系统中不平衡、联轴器不对中和磁轴承的动态参数，并针对测量误差测试了所提出方法的有效性和鲁棒性. Lees等^[5]提出一种利用转子静态和动态响应特性变化的诊断方法检测联轴器不对中故障，并发现，随着不对中程度的增加，转子轨道随不对中量的增加而趋于平坦.

对于轴承不对中故障，Kumar等^[6]发现磁轴承不对中会导致转子振动幅值增大，但对于转子系统仅考虑了平行不对中，实际中每个支撑转子的磁轴承不对中量并不是完全相同的. Kärkkäinen等^[7]研究磁轴承-转子系统轴承不对中产生的动力学特性，并给出了水平和垂直方向上的不对中量，发现水平方向的不对中严重影响转子系统稳定性. Kumar等^[8]研究了磁轴承-转子系统不同偏心距和不对中量大小对转子系统的影响，发现磁轴承不对中和圆盘不平衡会产生更高的电流和位移幅值. Kumar等^[9]还提出了一种利用FFT将磁轴承处的位移和电流时域响应转变为频域来估计不平衡故障参数和磁轴承不对中量，发现该算法在多个转子转速下，对噪声信号和建模误差具有鲁棒性和有效性，但未考虑陀螺效应影响下转子的轴向振动. Tiwari等^[10]分析磁轴承不对中以及传感器不对中所引起的振动特性，并使用陀螺动态降阶方法简化柔性转子-磁轴承不对中系统的数学模型.

随着转子转速的增加，磁轴承系统的陀螺效应更加明显. 实际工程中发生过多起汽轮发电机组轴向振动偏大的问题，影响机组的安全运行. 因此，在高速转动条件下，轴径向耦合问题不容忽视. 汤华涛等^[11]分析了不平衡转子会发生轴向振动的来源，利用多体系统传递矩阵法得到转子轴向力的推导方法，并搭载实验平台验证了转子的横向运动效应. 为防止轴向振动过高，朱熀秋等^[12]提出了三自由度混合磁轴承来抑制转子轴向传动. 汪希平等^[13]对五自由度电磁轴承转子系统进行分析，并搭建试验平台验证所提出转子系统动力学分析方法的可行性. Darpe等^[14]研究周期冲击载荷作用下裂纹转子引起的轴径向耦合振动，从径向振动的频谱中发现，裂纹会引起轴径向耦合振动. 但在工程实际中，转子-轴承系统往往会存在随机扰动. 为保证转子受到扰动情况下可以稳定运行，邹博等^[15]分析不同冲击激励下磁悬浮转子系统的响应，发现冲击时间和转速对转子振动具有显著影响. 此外，在相同冲击幅值下，矩形冲击相对于半正弦冲击，会使转子产生更明显的幅值响应. 陈伟等^[16]研究航空发动机转子系统受到冲击激励下的振动特性，并进行试验研究，发现转子系统瞬态振动响应加剧并会激起转子的正反进动和横向振动模态. 因此，需要感知机器的动态特性并识别故障，以便安全、平稳和有效地运行机器.

基于以上分析，冲击扰动下磁轴承固有不对中故障对转子系统的影响不可忽视. 本文研究对象为径向磁轴承和混合磁轴承支撑的五自由度转子系统，考虑冲击激励下和轴径向耦合作用下平行角度混合不对中对转子系统的影响，并利用频域下的最小二乘法计算转子-磁轴承系统的不对中量，使用冲击扰动和噪声验证该算法的鲁棒性，为该类转子的故障诊断和处理提供了参考.

1. 五自由度轴承-转子系统模型

1.1 系统模型

本文以刚性双偏置圆盘转子-磁轴承系统为分析对象，该旋转部件由转轴和圆盘组成，左右两端分别由径向磁轴承和混合磁轴承所支承，具体结构和参数如所示. 其中， ${A_1}{A_1}$ ， ${A_2}{A_2}$ 分别为径向磁轴承和混合磁轴承径向中心轴； ${\delta _{x1}}$ 、 ${\delta _{x2}}$ 分别为两磁轴承与转子旋转轴 ${\textit{z}}$ 的不对中量； ${A_3}{A_3}$ 为混合磁轴承轴向方向中心轴，其不对中量为 ${\delta _{{\textit{z}}2}}$ ； $G$ 为转子重心； ${a_1}$ 、 ${a_2}$ 分别为径向磁轴承和混合磁轴承距转子重心的距离； ${b_1}$ 、 ${b_2}$ 分别为圆盘1、2距转子重心的距离.

图 1 刚性双偏置圆盘转子-磁轴承系统模型

Figure 1. Model of a rigid double offset disc rotor-magnetic bearing system

下载: 全尺寸图片幻灯片

根据几何关系有

$\left\{\begin{gathered} u_x=u_{x 2}+\left(u_{x 1}-u_{x 2}\right) a_2 / l=\left(a_2 u_{x 1}+a_1 u_{x 2}\right) / l, \\ u_y=u_{y 2}+\left(u_{y 1}-u_{y 2}\right) a_2 / l=\left(a_2 u_{y 1}+a_1 u_{y 2}\right) / l, \end{gathered}\right.$

(1)

$\left\{\begin{gathered} \varphi_y=\left(u_{x 2}-u_{x 1}\right) / l, \\ \varphi_x=\left(u_{y 1}-u_{y 2}\right) / l, \end{gathered}\right.$

(2)

式中： $l$ 为转子长度， ${u_x}$ 、 ${u_y}$ 分别为转子重心处在x、y向的平动位移， ${\varphi _y}$ 、 ${\varphi _x}$ 分别为转轴在xOz平面和yOz平面内的横向角位移， ${u_{x1}}$ （ ${u_{y1}}$ ）、 ${u_{x2}}$ （ ${u_{y2}}$ ）分别为转子在径向磁轴承和混合磁轴承处x（y）方向的平动位移.

1.2 转子受力分析

假设转子在圆盘1处x、y方向分别受到加速度幅值为 $R{{g}}$ ，持续时间为 $t$ 的半正弦冲击激励和矩形脉冲激励^[15]. 其中，R为幅值系数. 以半正弦冲击为例，转子在x、y方向受到的冲击力为

$\left\{\begin{gathered} {f_x} = m{\rho _1}\ddot u\left( t \right) , \\ {f_y} = m{\rho _2}\ddot u\left( t \right) , \end{gathered}\right.$

(3)

式中： ${\rho _1}$ 、 ${\rho _2}$ 为衰减系数， ${\rho _1} = 0.6$ ， ${\rho _2} = 0.3$ ； $m$ 为转子的质量； $\ddot u(t)$ 为半正弦加速度激励信号，如式（4）.

$\ddot{u}(t)=\left\{\begin{array}{l} Rg \sin\;{\varOmega }_{1}t,\;\;\;\text{ }{t}_{0}\leqslant t\leqslant {t}_{1}, \\ 0, \; \; \; \text{ }t < {t}_{0}\;或\;t > {t}_{1},\end{array} \right. \\$

(4)

式中： ${\varOmega _1}$ 为冲击激励信号的角频率， ${t_0}$ 及 ${t_1}$ 分别为冲击开始、冲击结束时间.

当转子以角速度 $\omega$ 旋转时，x、y、z方向产生的不平衡力组成的矩阵为

${{\boldsymbol{f}}_{{\mathrm{unb}}}}=\left[ {\begin{array}{*{20}{c}} {{f_{{\mathrm{unb}}xk}}} \\ {{f_{{\mathrm{unb}}yk}}} \\ {{f_{{\mathrm{unb}}{\textit{z}}k}}} \end{array}} \right] = {m_{{\mathrm{d}}k}}{\omega ^2}\left[ {\begin{array}{*{20}{c}} {{e_k}\cos(\omega t + {\beta _k})} \\ {{e_k}\sin(\omega t + {\beta _k})} \\ {{e_{{\text{z}}k}}} \end{array}} \right] ，$

(5)

式中： ${m_{{\mathrm{d}k}}}$ 为圆盘质量， $\;{\beta _k}$ 为圆盘不平衡对应相位， ${e_k}$ 和 ${e_{{\text{z}}k}}$ 分别为径向和轴向偏心距，下标 $k = 1, 2$ 分别为圆盘1和圆盘2.

1.3 轴承电磁力计算

当转子受到外界干扰发生微小振动时，导致一边磁极与转子的间隙为 ${s_0} + {u_{nq}}$ ，另一边磁极与转子的间隙为 ${s_0} - {u_{nq}}$ （其中： ${s_0}$ 为气隙； $q = 1, 2$ ，分别代表径向磁轴承和混合磁轴承；n=x, y, z）. 为使转子系统平稳运行，磁轴承控制系统会输出动态控制电流 ${i_{nq}}$ ，并与偏置电流 ${i_0}$ 共同作用. 对转子产生的电磁力为

${F_{nq}} = {k_q}\left\{ {\dfrac{{{{\left( {{i_0} + {i_{nq}}} \right)}^2}}}{{{{\left( {{s_0} - {u_{nq}}} \right)}^2}}} - \dfrac{{{{\left( {{i_0} - {i_{nq}}} \right)}^2}}}{{{{\left( {{s_0} + {u_{nq}}} \right)}^2}}}} \right\}，$

(6)

式中： ${{{k_q} = {\mu _0}N_q^2{S_q}\cos \;\alpha } \mathord{\left/ {\vphantom {{{k_q} = {\mu _0}N_q^2{S_q}\cos \alpha } 4}} \right. } 4}$ ， ${\mu _0}$ 为磁导率， $N_q$ 为线圈匝数， ${S_q}$ 为磁极的面积， $\alpha$ 为相邻磁极的夹角.

由于 ${u_{nq}} \ll {s_0}$ ，对式（6）进行一阶泰勒展开：

${F_{nq}} = {K_{{\mathrm{a}}nq}}{u_{nq}} + {K_{{\mathrm{i}}nq}}{i_{nq}}，$

(7)

式中： ${K_{{\rm{a}}nq}} = 4{k_q}{{i_0^2} \mathord{\left/ {\vphantom {{i_0^2} {s_0^3}}} \right. } {s_0^3}}$ ； $K_{{\mathrm{i}}nq}$ 为力位移刚度， ${K_{{\text{i}}nq}} =$ $4{k_q}{{{i_0}} \mathord{\left/ {\vphantom {{{i_0}} {s_0^2}}} \right. } {s_0^2}}$ ，为力电流刚度.

转子和磁轴承轴线不对中侧视图如所示. 径向和混合磁轴承中轴线偏离刚性转子的不对中量为 ${\delta _{nq}}$ . 导致磁轴承与转子的初始气隙变为 ${s_0} + {\delta _{nq}}$ 和 ${s_0} - {\delta _{nq}}$ ，则

${F_{nq}} = {k_q}\left( {\dfrac{{{{\left( {{i_0} + {i_{nq}}} \right)}^2}}}{{{{\left( {{s_0} - {\delta _{nq}} + {u_{nq}}} \right)}^2}}} - \dfrac{{{{\left( {{i_0} - {i_{nq}}} \right)}^2}}}{{{{\left( {{s_0} + {\delta _{nq}} - {u_{nq}}} \right)}^2}}}} \right) .$

(8)

由于 ${u_{nq}} \ll ({s_0} - {\delta _{nq}}) < ({s_0} + {\delta _{nq}})$ ，忽略高阶小项 ${u_{nq}}$ 、 $i_{nq}^2$ 、 ${i_{nq}}{u_{nq}}$ 、 $i_{nq}^2{u_{nq}}$ ，电磁力可写成式（9）形式.

${F_{nq}} = {k_q}\left( {\dfrac{{4{s_0}i_0^2{u_{nq}} + 4{i_0}{i_{nq}}\left( {s_0^2 + \delta _{nq}^2} \right) + 4{s_0}{\delta _{nq}}i_0^2}}{{{{\left( {{s_0} - {\delta _{nq}}} \right)}^2}{{\left( {{s_0} + {\delta _{nq}}} \right)}^2}}}{\text{ }}} \right).$

(9)

所以，当转子不对中时，磁轴承对刚性转子的电磁力可表示为

${F_{{\mathrm{m}}1nq}} = {k_{{\mathrm{m1a}}nq}}{u_{nq}} + {k_{{\mathrm{m1i}}nq}}{i_{nq}} + {f_{{\mathrm{m1c}}nq}}，$

(10)

式中： ${k_{{\mathrm{m1a}}nq}} = {{{K_{{\mathrm{a}}nq}}}}/{{{{\left( {1 - \zeta _{nq}^2} \right)}^2}}}$ ，为位移刚度； ${k_{{\mathrm{m1i}}nq}} = {{{K_{{\mathrm{i}}nq}}\left( {1 + \zeta _{nq}^2} \right)}}/{{{{\left( {1 - \zeta _{nq}^2} \right)}^2}}}$ ，为电流刚度； ${\zeta _{nq}} = {{{\delta _{nq}}}}/{{{s_0}}}$ ； ${f_{{\mathrm{m1c}}nq}}$ $={{4{k_q}i_0^2{\zeta _{nq}}}}/{\left({s_0^2{{\left( {1 - \zeta _{nq}^2} \right)}^2}}\right)}$ ，为径向常力.

因此，转子左右磁轴承的电磁力可以表示为

$\begin{split} &\left[\begin{gathered} {F_{n1}} \\ {F_{n2}} \\ \end{gathered} \right]= \left[ {\begin{array}{*{20}{c}} {{k_{{\mathrm{m1a}}n1}}}&&0 \\ 0&&{{k_{{\mathrm{m1a}}n2}}} \end{array}} \right]\left[ \begin{gathered} {u_{n1}} \\ {u_{n2}} \\ \end{gathered} \right] + \\ &\quad{\text{ }}\left[ {\begin{array}{*{20}{c}} {{k_{{\mathrm{m1i}}n1}}}&&0 \\ 0&&{{k_{{\mathrm{m1i}}n2}}} \end{array}} \right]\left[ \begin{gathered} {i_{n1}} \\ {i_{n2}} \\ \end{gathered} \right] + \left[ \begin{gathered} {f_{{\mathrm{m1c}}n1}} \\ {f_{{\mathrm{m1c}}n2}} \\ \end{gathered} \right]. \\ \end{split}$

(11)

为能够准确定量计算系统固有不对中量大小，需要引入一个额外偏差条件下的不对中量 ${{\varDelta}}$ ，可以通过改变电磁轴承作动器的高度等方法实现. 引入该量后，转子的轴线与磁轴承之间不对中量变为 ${\delta _{nq}} + {\varDelta _{nq}}$ ， ${\varDelta _{nq}}$ 为径向磁轴承和混合磁轴承各个方向的额外不对中量. 那么，额外偏差条件下轴承在刚性转子上的电磁力可以表示为

${F_{{\mathrm{m2}}nq}} = {k_{{\mathrm{m2a}}nq}}{u_{nq}} + {k_{{\mathrm{m2i}}nq}}{i_{nq}} + {f_{{\mathrm{m2c}}nq}},$

(12)

式中： ${k_{{\mathrm{m2a}}nq}} = {{{K_{{\mathrm{a}}nq}}}}/{{{{\left( {1 - p_{nq}^2} \right)}^2}}}$ ； ${k_{{\mathrm{m2i}}nq}} = {{{K_{{\mathrm{i}}nq}}\left( {1 + p_{nq}^2} \right)}}/ {{{{\left( {1 - p_{nq}^2} \right)}^2}}}$ ， ${p_{nq}} = {\zeta _{nq}} + {\varDelta _{{\mathrm{a}}nq}}$ ， ${\varDelta _{{\mathrm{a}}nq}} = \dfrac{{{\varDelta _{nq}}}}{{{s_0}}}$ ； ${f_{{\mathrm{m2c}}nq}} = {{4{k_q}i_0^2{p_{nq}}}}/ {\left({s_0^2{{\left( {1 - p_{nq}^2} \right)}^2}}\right)}$ .

1.4 考虑轴向径向耦合的转子轴向力

由于偏置圆盘的作用会使得转子旋转时产生轴向运动，当转速上升，即使刚度很大的转子也会由于离心力的作用发生微弹性变形， ${\rho _{\mathrm{w}}} = e{\mu ^2}/\left( {1 - } \right.$ $\left. {{\mu ^2}} \right)$ ，其中： $\mu = {\omega / {{\omega _0}}}$ ， ${\omega _0}$ 为系统固有频率， $e$ 为偏心距（e₁、e₂）. 如所示，图中：( ${x_w}$ ， ${y_w}$ )、( ${x_s}$ ， ${y_s}$ )分别为形心w和重心s坐标. 以圆盘1为例，由于圆盘的自重影响，会产生 $x$ 方向的静挠度， ${\rho _{\mathrm{G}}} = {m_{{\mathrm{d}}1}}g({a_1} -$ ${b_1})({b_1} + {a_2}) [{l^2} - {({a_1} - {b_1})^2} - {({b_1} + {a_2})^2}]$ $/(6EIl)$ ，其中： $E$ 为弹性模量， $I$ 为惯性矩.

图 2 圆盘1形心和重心位置关系

Figure 2. Relationship between shape center and gravity center for disc 1

下载: 全尺寸图片幻灯片

建立轴截面形心运动微分方程为

$\left\{ \begin{gathered} {{\ddot x}_w} + \omega _0^2{x_w} = {e_1}{\omega ^2}\cos(\omega t + {\beta _1}) ，\\ {{\ddot y}_w} + \omega _0^2{y_w} = {e_1}{\omega ^2}\sin(\omega t + {\beta _1}). \\ \end{gathered}\right.$

(13)

求解微分方程，并根据形心和重心坐标关系可得圆盘重心的运动方程为

$\left\{\begin{array}{l} {x_s} = {{{e_1}} / {(1 - {\mu ^2})}}\cos(\omega t + {\beta _1}) - {\rho _{\mathrm{G}}} = \\ \quad{\rho _{\mathrm{s}}}\cos(\omega t + {\beta _1}) - {\rho _{\mathrm{G}}}， \\ {y_s} = {{{e_1}} / {(1 - {\mu ^2})}}\sin(\omega t + {\beta _1}) = \\ \quad{\rho _{\mathrm{s}}}\sin(\omega t + {\beta _1}), \end{array}\right.$

(14)

式中： ${\rho _{\mathrm{s}}} = {{{e_1}} / {(1 - {\mu ^2})}}$ .

根据动量矩定理得到空间坐标系中x、y方向的分力为

$\left\{ \begin{gathered} {F_{{\mathrm{l}}x}} = \dfrac{{{M_{{\mathrm{l}}x}}}}{{{a_1} - {b_1}}} = \dfrac{{{J_{{\mathrm{l}}1}}{{\ddot \theta }_x}}}{{{a_1} - {b_1}}} ， \\ {F_{{\mathrm{l}}y}} = \dfrac{{{M_{{\mathrm{l}}y}}}}{{{a_1} - {b_1}}} = \dfrac{{{J_{{\mathrm{l}}1}}{{\ddot \theta }_y}}}{{{a_1} - {b_1}}} ， \\ \end{gathered} \right.$

(15)

式中：M_lx、M_ly分别为x、y方向的力矩； ${\theta _x}$ 、 ${\theta _y}$ 分别为圆盘重心和左支撑点的连线与轴向方向的夹角在zOx平面和zOy平面的投影, ${{{\theta _x} = {y_s}} / {\left( {{a_1} - {b_1}} \right)}}$ ， ${\theta _y} = {{x_s}} / {\left( {{a_1} - {b_1}} \right)}$ ； ${J_{{\mathrm{l}}1}}$ 为圆盘1绕左支撑点的转动惯量； ${\ddot \theta _x}$ 、 ${\ddot \theta _y}$ 分别为 ${\theta _x}$ 和 ${\theta _y}$ 对时间的二次导数.

将各方向的分力投影到z轴可得圆盘1绕左支撑点转动产生的轴向力为

$\begin{split}&{F_{{\mathrm{l}}{\textit{z}}1}} = \dfrac{{{J_{{\mathrm{l}}1}}{\rho _s}{\rho _{\mathrm{G}}}{\omega ^2}}}{{{{({a_1} - {b_1})}^3}}}\cos(\omega t + {\beta _1}) -\\ & \quad\dfrac{{{J_{{\mathrm{l}}1}}\rho _s^2{\omega ^2}}}{{{{({a_1} - {b_1})}^3}}}\cos(2\omega t + {\beta _1}).\end{split}$

(16)

同理可得圆盘1绕右支撑点转动产生的轴向力为

$\begin{split}&{F_{{\mathrm{r}}{\textit{z}}1}} = \dfrac{{{J_{{\mathrm{r}}1}}{\rho _s}{\rho _{\mathrm{G}}}{\omega ^2}}}{{{{({a_2} + {b_1})}^3}}} \cos (\omega t + {\beta _1}) - \\&\quad\dfrac{{{J_{{\mathrm{r}}1}}\rho _s^2{\omega ^2}}}{{{{({a_2} + {b_1})}^3}}}\cos (2\omega t + {\beta _1})， \end{split}$

(17)

式中： ${J_{r1}}$ 为圆盘1绕右支撑点的转动惯量.

双偏置圆盘转动产生的轴向合力为

${F_{\textit{z}}} = ({F_{{\mathrm{l}}{\textit{z}}1}} - {F_{{\mathrm{r}}{\textit{z}}1}}) - ({F_{{\mathrm{l}}{\textit{z}}2}} - {F_{{\mathrm{r}}{\textit{z}}2}})，$

(18)

式中： ${F_{{\mathrm{l}}{\textit{z}}2}}$ 、 ${F_{{\mathrm{r}}{\textit{z}}2}}$ 分别为圆盘2绕左、右支撑点转动产生的轴向力.

2. 耦合效应下磁轴承转子系统动力学建模及固有不对中量推导

2.1 动力学建模

以刚性双偏置圆盘转子-磁轴承系统模型为分析对象，建立考虑旋转机械系统的惯性力、不平衡力、陀螺力矩、冲击力等因素下不对中转子-磁轴承系统平衡方程为

$\qquad\quad\;\; \left\{ \begin{array}{l} m{a_2}{{\ddot u}_x} = {I_{\rm{d}}}{{\ddot \varphi }_y} - ({I_{{\rm{p}}1}} + {I_{{\rm{p}}2}}){{\ddot \varphi }_x}\omega + \left\{ {{f_{{\rm{unb}}x1}}({a_{\text{2}}} + {b_1}) + {f_{{\rm{unb}}x2}}({a_{\text{2}}} - {b_2})} \right\} + {F_{{\rm{m1}}x1}}l + {f_x}({a_{\text{2}}} + {b_1}), \\ m{a_2}{{\ddot u}_y} = - {I_{\rm{d}}}{{\ddot \varphi }_x} + ({I_{{\rm{p}}1}} + {I_{{\rm{p}}2}}){{\ddot \varphi }_y}\omega + \left\{ {{f_{{\rm{unb}}y1}}({a_{\text{2}}} + {b_1}) + {f_{{\rm{unb}}y2}}({a_{\text{2}}} - {b_2})} \right\} + {F_{{\rm{m1}}y1}}l + {f_y}({a_{\text{2}}} + {b_1}), \\ m{a_1}{{\ddot u}_x} = - {I_{\rm{d}}}{{\ddot \varphi }_y} + ({I_{{\rm{p}}1}} + {I_{{\rm{p}}2}}){{\ddot \varphi }_x}\omega + \left\{ {{f_{{\rm{unb}}x1}}({a_{\text{1}}} - {b_1}) + {f_{{\rm{unb}}x2}}({a_{\text{1}}} + {b_2})} \right\} + {F_{{\rm{m1}}x2}}l + {f_x}({a_{\text{1}}} - {b_1}), \\ m{a_1}{{\ddot u}_y} = {I_{\rm{d}}}{{\ddot \varphi }_x} - ({I_{{\rm{p}}1}} + {I_{{\rm{p}}2}}){{\ddot \varphi }_y}\omega + \left\{ {{f_{{\rm{unb}}y1}}({a_{\text{1}}} - {b_1}) + {f_{{\rm{unb}}y2}}({a_{\text{1}}} + {b_2})} \right\} + {F_{{\rm{m1}}y2}}l + {f_y}({a_{\text{1}}} - {b_1}), \\ m{{\ddot u}_{{\textit{z}}{\text{2}}}} = {m_{{\mathrm{d}}1}}{\omega ^2}{e_{{\textit{z}}1}} + {m_{{\mathrm{d}}2}}{\omega ^2}{e_{\textit{z}}}_2 + {F_{{\rm{m1}}{\textit{z}}2}} + {F_{\textit{z}}}， \end{array} \right.$

(19)

式中： ${I_{{\mathrm{p}}1}}$ 、 ${I_{{\mathrm{p}}2}}$ 分别为圆盘1、2的惯性矩， ${I_{\mathrm{d}}}$ 为转子惯性矩， $u_{{\textit{z}}2}$ 为轴向位移， $F_{{\mathrm{m}}1{\textit{z}}2}$ 为混合磁铁z方向的电磁，同样可通过式（10）得到.

将平衡方程写成矩阵形式，如式（20）.

${\boldsymbol{M}}{\ddot {\boldsymbol{u}}_{nq}} - \omega {\boldsymbol{G}}{\dot {\boldsymbol{u}}_{nq}} = {{\boldsymbol{f}}_{{\mathrm{unb}}}} + {{\boldsymbol{F}}_{{\mathrm{m1}}nq}}，$

(20)

式中： ${\boldsymbol{M}}$ 、 ${\boldsymbol{G}}$ 为转子系统的质量矩阵、陀螺矩阵， ${{\boldsymbol{f}}_{{\mathrm{unb}}}}$ 为圆盘不平衡力向量， ${{\boldsymbol{F}}_{{\mathrm{m1}}nq}}$ 为不对中条件下磁轴承电磁力向量， ${{\boldsymbol{u}}_{nq}}$ 为刚性转子振动位移向量.

2.2 固有不对中量推导

本文采用FFT算法将时域响应转换为频域响应，该算法可以高效分析信号频域特征，频域信号中的主要频率成分是机械设备故障特征提取的主要依据. 由于对时域信号的随机采样而导致相位会发生改变，因此，使用正交参考相位信号^[17]来获得位移和电流响应的正确相位.

转子的时域位移响应和磁轴承的时域电流响应以及圆盘的离心力可以在频域中表示，如式（21）^[6].

$\left\{ \begin{array}{l} {u_{xq}}\left( t \right) + {\mathrm{j}}{u_{yq}}\left( t \right) = \left( {{R_{{\mathrm{i}}xq{\text{,}}{\mathrm{Re}}}}\left( \omega \right) + {\mathrm{j}}{R_{{\mathrm{i}}yq{\text{,}}{\mathrm{Im}}}}\left( \omega \right)} \right){{\mathrm{e}}^{{\mathrm{j}}i\omega t}}, \\ {i_{xq}}\left( t \right) + {\mathrm{j}}{{{i}}_{yq}}\left( t \right) = \left( {{I_{{\mathrm{i}}xq{\text{,}}{\mathrm{Re}}}}\left( \omega \right) + {\mathrm{j}}{I_{{\text{i}}yq{\text{,}}{\mathrm{Im}}}}\left( \omega \right)} \right){{\mathrm{e}}^{{\mathrm{j}}i\omega t}},\\ {m_{\mathrm{d}}}_k{e_k}{\omega ^2}{{\mathrm{e}}^{{\mathrm{j}}{\beta _k}}} = {m_{\mathrm{d}}}_k{\omega ^2}\left( {{e_k}\cos\; {\beta _k} + {\mathrm{j}}{e_k}\sin \;{\beta _k}} \right) =\\ \quad{m_{\mathrm{d}}}_k{\omega ^2}\left( {{e_{k,{\mathrm{Re}}}} + {\mathrm{j}}{e_{k,{\mathrm{Im}}}}} \right) \text{} ，\end{array}\right.$

(21)

式中： ${R_{{\mathrm{i}}xq{\text{,}}{\mathrm{Re}}}}\left( \omega \right)$ 、 ${R_{{\mathrm{i}}yq{\text{,}}{\mathrm{Im}}}}\left( \omega \right)$ 、 ${I_{{\mathrm{i}}xq{\text{,}}{\mathrm{Re}}}}\left( \omega \right)$ 、 ${I_{{\text{i}}yq{\text{,}}{\mathrm{Im}}}}\left( \omega \right)$ 分别为时域位移响应 ${u_{xq}}\left( t \right)$ 、 ${u_{yq}}\left( t \right)$ 和电流响应 ${i_{xq}}\left( t \right)$ 、 ${i_{yq}}\left( t \right)$ 在x、y方向下的频域表示形式； ${e_{k,{\mathrm{Re}}}}$ 、 ${e_{k,{\mathrm{Im}}}}$ 为圆盘偏心距在频域下的表示形式； $i = 0$ 和 $i = 1$ ，即对应系统方程的直流形式和基波形式，方程的实、虚部分别对应.

2.2.1 径向方向不对中量δ_x₁、δ_y₁、δ_x₂、δ_y₂求解推导

为求解系统不对中量的大小，将式（21）代入式（20）中. 频域下的式（21）可以重新排列成矩阵形式，令圆盘偏心距为 ${e}_{1,{\rm{Re}}}、{e}_{2,{\rm{Re}}}、{e}_{1,{\rm{Im}}}、\text{ }{e}_{2,{\rm{Im}}}$ ，不对中情况下力位移刚度 ${{k_{{\rm{m1a}}nq}}}$ 、力电流刚度 ${{k_{{\rm{m1i}}nq}}}$ 、磁轴承的径向常力 ${{f_{{\rm{m1c}}nq}}}$ 作为未知量，其他作为已知量，建立矩阵方程：

${\boldsymbol{Kx_1 = J }},$

(22)

式中：K、J为矩阵的系数矩阵与常数矩阵， ${\boldsymbol{x}}_1 = ({{e_{1,{\mathrm{Re}}}}} ,$ ${e_{2,{\mathrm{Re}}}},\;{e_{1,{\mathrm{Im}}}},\;{e_{2,{\mathrm{Im}}}},\;{k_{{\text{m}}{\mathrm{1a}}x1}},\;{k_{{\mathrm{m1a}}y1}},\;{k_{{\mathrm{m1a}}x2}},\;{k_{{\mathrm{m1a}}y2}}{k_{{\mathrm{m1i}}x1}}, \;{k_{{\mathrm{m1i}}y1}},$ ${k_{{\mathrm{m1i}}x2}},\;{k_{{\mathrm{m1i}}y2}},\;{f_{{\mathrm{m1c}}x1}},\;{f_{{\mathrm{m1c}}y1}},\;{f_{{\mathrm{m1c}}x2}},\;{ {{f_{{\mathrm{m1c}}y2}}} )^{\mathrm{T}}}$ .

未知数有16个，但是方程仅有8个，选择远离临界转速的2种转速 ${\omega _1}$ 、 ${\omega _2}$ 运行. 不对中条件下和额外偏差条件下的回归矩阵均可以表示为

$\left[ {\begin{array}{*{20}{c}} {{\boldsymbol{K}}({\omega _1})} \\ {{\boldsymbol{K}}({\omega _2})} \end{array}} \right]{\boldsymbol{x}} = \left[ {\begin{array}{*{20}{c}} {{\boldsymbol{J}}({\omega _1})} \\ {{\boldsymbol{J}}({\omega _2})} \end{array}} \right] .$

(23)

利用最小二乘法求解式（23），得到

${\boldsymbol{x}} = {({{\boldsymbol{K}}^{\rm{T}}}{\boldsymbol{K}})^{ - 1}}{{\boldsymbol{K}}^{\rm{T}}}{\boldsymbol{J}}.$

(24)

对于回归矩阵（23）的求解，至少需要2种转速. 可以通过采取多种转速以实现对未知量x更全面估计. 得到未知量矩阵 ${\boldsymbol{x}}$ 后，提取不同情形下力位移刚度为

$\dfrac{{{k_{{\rm{m1a}}nq}}}}{{{k_{{\mathrm{m2a}}nq}}}} = \dfrac{{{{\left( {1 - p_{nq}^2} \right)}^2}}}{{{{\left( {1 - \zeta _{nq}^2} \right)}^2}}} .$

(25)

通过以上推导，可得到

$\begin{split} {\delta _{nq}} = \dfrac{{2{\varDelta _{nq}} + \sqrt {4\varDelta _{nq}^2 - 4\left( {\sqrt {\dfrac{{{k_{{\rm{m1a}}nq}}}}{{{k_{{\rm{m2a}}nq}}}}} - 1} \right)\left( {1 - \sqrt {\dfrac{{{k_{{\rm{m1a}}nq}}}}{{{k_{{\rm{m2a}}nq}}}}} - \varDelta _{nq}^2} \right)} }}{{2\left( {\sqrt {\dfrac{{{k_{{\rm{m1a}}nq}}}}{{{k_{{\rm{m2a}}nq}}}}} - 1} \right)}}{s_0}. \end{split}$

(26)

2.2.2 轴向方向不对中量 ${\delta _{{\textit{z}}2}}$ 求解推导

同理，为求解轴向方向的不对中量，将轴向时域位移和电流响应与径向 $x$ 方向的位移和电流响应分别组成复数形式，并在频域中表示为

$\left\{ \begin{gathered} {u_{xq}}\left( t \right) + {\rm{j}}{u_{{\textit{z}}q}}\left( t \right) = \left( {{R_{{\rm{i}}xq{\text{,}}{\rm{Re}}}}\left( \omega \right) + {\mathrm{j}}{R_{{\text{i}}{\textit{z}}q{\text{,}}{\rm{Im}}}}\left( \omega \right)} \right){{\mathrm{e}}^{{\rm{j}}i\omega t}}，\\ {i_{xq}}\left( t \right) + {\rm{j}}{i_{{\textit{z}}q}}\left( t \right) = \left( {{I_{{\rm{i}}xq{\text{,}}{\rm{Re}}}}\left( \omega \right) + {\mathrm{j}}{I_{i{\textit{z}}q{\text{,}}{\rm{Im}}}}\left( \omega \right)} \right){{\mathrm{e}}^{{\rm{j}}i\omega t}}， \end{gathered}\right.$

(27)

式中： ${R_{{\text{i}}{\textit{z}}q{\text{,}}{\rm{Im}}}}\left( \omega \right)$ 、 ${I_{{\rm{i}}{\textit{z}}q{\text{,}}{\rm{Im}}}}\left( \omega \right)$ 分别为轴向时域位移响应 ${u_{{\textit{z}}q}}\left( t \right)$ 和电流响应 ${i_{{\textit{z}}q}}\left( t \right)$ 在频域下的表示形式.

同样，将圆盘偏心距 $\left( {{e_{{{{\textit{z}}}}1}},{e_{{{{\textit{z}}}}2}}} \right)$ ，固有不对中情况下的力位移刚度 ${k_{{\rm{m1a}}{\textit{z}}2}}$ 、力电流刚度 ${k_{{\rm{m1i}}{\textit{z}}2}}$ 、电磁力常力项 ${f_{{\rm{m1c}}{\textit{z}}2}}$ 以及额外偏差条件下的力位移刚度 ${k_{{\rm{m2a}}{\textit{z}}2}}$ 、力电流刚度 ${k_{{\rm{m2i}}{\textit{z}}2}}$ 、电磁力常力项 ${f_{{\rm{m2c}}{\textit{z}}2}}$ 作为未知量，其他作为已知量，建立矩阵方程：

${\boldsymbol{K}}{\boldsymbol{x}}_2 = {\boldsymbol{J}} ，$

(28)

式中： ${\boldsymbol{x }}_2 = {\left( {{e_{{{\textit{z}}}1}},{e_{{{{\textit{z}}}}2}},{k_{{\rm{m1a}}{\textit{z}}2}},{k_{{\rm{m1i}}{\textit{z}}2}},{f_{{\rm{m1c}}{\textit{z}}2}},{k_{{\rm{m2a}}{\textit{z}}2}},{k_{{\rm{m2i}}{\textit{z}}2}},{f_{{\rm{m2c}}{\textit{z}}2}}} \right)^{\rm{T}}}$ .

为求解该矩阵未知量，以2种转速运行，并利用最小二乘法求得轴向力位移刚度，进一步得到轴向不对中量大小.

3. 系统动力学特性分析

3.1 SIMULINK仿真模型

为获得系统的时域位移和电流响应，使用比例-积分-微分（PID）方法控制轴承线圈电流，比例因子K_P=5 500 A/m，积分因子K_I=8 000 A/（m·s），微分因子K_D=3 A·s/m. 建立不对中条件下SIMULINK模型，如图3所示. 图中：u₁、u₂为simulink中函数常量，f_c为磁轴承常力，F₁、F₂分别为不平衡力和冲击力，K_i、K_a分别为力电流、力位移刚度矩阵. 表1为系统结构参数与仿真参数. 使用定步长0.0001的四阶龙格库塔（ode4）微分方程求解器求解动力学方程.

图 3 动力学模型仿真

Figure 3. Dynamics model simulation

下载: 全尺寸图片幻灯片

表 1 系统结构参数与仿真参数

Table 1. System structure parameters and simulation parameters

参数	数值	参数	数值
m/kg l/m	3.910 0.400	s₀/mm	0.400
m_d1/kg	1.065	δ_x1，δ_x2/mm	0.140，0.150
m_d2/kg	2.081	δ_y1，δ_y2/mm	0.160，0.145
e₁，e₂/μm	80，100	δ_z2，Δ_x1/mm	0.130，0.100
${e_{{\textit{z}}1}}，{e_{{\textit{z}}2}}$ /μm	8，10	Δ_y1，Δ_x2/mm	0.130，0.120
β₁，β₂/（°）	20，30	Δ_y2，Δ_z2/mm	0.110，0.150
I_d/（kg·m²）	0.0455	K_P/（A·m⁻¹）	5500
a₁/mm	0.226	K_I/（A·（m·s）⁻¹）	8000
a₂/mm	0.174	K_D/（A·s·m⁻¹）	3
b₁/mm	0.126	R t/ms	20 6
I_p1/（kg·m²）	0.0019	b₂/mm	0.074
I_p2/（kg·m²）	0.0059

下载: 导出CSV

| 显示表格

3.2 转子时域位移和电流响应

转子在工作过程中，难免会受到各种扰动冲击作用. 因此，考虑冲击作用下轴承转子系统不对中量的精确求解非常必要. 图4为转子受半正弦冲击激励和矩形脉冲激励后位移响应，将矩形脉冲激励幅值设置为半正弦冲击幅值的一半，用于判断冲击幅值对不对中量识别算法准确性的影响. 可以看出，转子受到冲击后位移幅值非常大，对矩形冲击，即使冲击幅值减小了一半，但转子运动状态仍然会发生较大变化，在不同冲击激励下，转子最终达到稳态后的轴心轨迹几乎一致.

图 4 转子受冲击后轴心轨迹

Figure 4. Axial trajectory of rotor subjected to impact

下载: 全尺寸图片幻灯片

根据表1所示的转子系统结构参数，计算得到系统前两阶临界转速分别为250.08、508.91 rad/s，假设考虑转子为刚性，故选择低于临界转速的2种转速进行不对中量计算，分别为188.40、219.80 rad/s，其对应转速频率为30、35 Hz.

图5为30 Hz和35 Hz转速时，径向磁轴承处转子位移和控制电流系统响应. 可以看出：1）转子由于存在不对中偏差，径向磁轴承处转子位移幅值高于对中时的转子位移幅值，且会产生更高的控制电流. 2）随着系统转速的增加，转子不对中引起的高位移幅值和高电流消耗特征更加明显. 3）完全对中情况下，转子位移和电流响应均在平衡位置处，不对中情况下，磁轴承电流偏离平衡位置处，这是由于电磁力常力项导致. 4）这些常力项也会使位移响应偏离平衡位置，但由于PID控制的补偿作用，通过输出控制电流，使转子恢复至平衡位置.

图 5 不同转速下径向磁轴承处位移、电流稳态响应

Figure 5. Steady-state response of displacement and current at radial magnetic bearing at different rotational speeds

下载: 全尺寸图片幻灯片

图6为30 Hz和35 Hz转速时，混合磁轴承处转子空间振动以及控制电流信号. 可以发现：不对中条件下，径向位移和控制电流均高于对中状态，随着转速的增加，转子位移和磁轴承控制电流均有增加；对于轴向方向，尽管轴向振动相比径向振动幅值较小，但会产生更高的控制电流，随着转速增加，这种轴向振动会愈加明显，因此，轴向振动同样不容忽视.

图 6 不同转速下混合磁轴承处位移、电流稳态响应

Figure 6. Steady-state response of displacement and current at hybrid magnetic bearing at different rotational speeds

下载: 全尺寸图片幻灯片

3.3 转子频域位移和电流响应

利用FFT将SIMULINK模型中求解得到的时域下稳态位移响应和电流响应变换成频域位移响应和频域电流响应，这些响应的实部和虚部可由复位移和复相位计算得到，并用于回归矩阵式（23）的求解，其部分数值如表2所示.

表 2 频域下的径向磁轴承处转子位移和控制电流响应

Table 2. Rotor displacement and control current response at radial magnetic bearing in frequency domain

频率/ Hz	i	位移		电流
频率/ Hz	i	幅值/A	相位/（°）	幅值/m	相位/（°）
30	0	3.69 × 10⁻⁹	−16.09	0.929	−132.73
30	1	1.36 × 10⁻⁴	−32.43	0.754	152.99
35	0	1.13 × 10⁻⁹	112.79	0.930	−133.20
35	1	1.83 × 10⁻⁴	−61.97	1.016	124.49

下载: 导出CSV

| 显示表格

在30 Hz时，转子位移和电流响应的幅值分别为1.36 × 10⁻⁴ m和0.929 A. 同样，35 Hz的转速下，位移和电流响应的幅值分别为1.83 × 10⁻⁴ m和1.016 A. 可以发现：随着系统转速的增加，位移和电流的幅值均有不同程度增加. 并且从电流幅值可以看出，系统具有很高的电流直流谐波. 这是由于转子系统初始不对中导致电磁轴承对静态转子产生控制电流. 利用这些信号特征能够实时监测系统是否存在不对中故障，为系统故障诊断和安全高效运行提供了参考.

3.4 系统不对中量的计算

在实际工程中，转子系统的响应容易受到系统本身噪声和附近环境各种噪声影响. 为在仿真过程中模拟真实情况，将1% ～ 5%的噪声添加到位移和电流响应中，利用最小二乘法求解回归矩阵得到磁轴承动态参数，进而得到系统的不对中量误差，如图7所示. 可以看到，该算法在计算δ_x1时容易受到噪声的影响，最大误差为18.3%，其余不对量的计算误差随着噪声信号的增加均有不同程度的增加. 在无噪声影响条件下，半正弦冲击扰动和矩形冲击扰动下计算的不对中量大小均如表3所示，这是由于不同冲击扰动后转子稳态轨迹几乎一致，导致2种扰动下不对中量估计值接近相同. 对比实际设定的不对中量和估计值大小，最大误差均在5%以内，可能由频谱泄漏等原因造成，具有较好的准确性. 综合来看，识别系统不对中量大小将会为在线矫正转子轴线失准提供依据，并能指导实验中转子的精准组装.

图 7 不同噪声影响下计算的不对中量误差百分比

Figure 7. Percentage of misalignment error calculated under different noise interferences

下载: 全尺寸图片幻灯片

表 3 无噪声影响下不对中量实际值和估计值比较

Table 3. Comparison of actual and estimated values of misalignment without noise interference

不对中量	实际值/mm	估计值/mm	误差/%
δ_x1	0.140	0.1388	−0.857
δ_y1	0.160	0.1575	−1.563
δ_x2	0.150	0.1439	−4.067
δ_y2	0.145	0.1456	0.414
δ_z2	0.130	0.1285	−1.154

下载: 导出CSV

| 显示表格

4. 结　论

本文对冲击激励下转子-磁轴承系统不对中条件下的系统动力学特性进行分析，同时利用FFT和最小二乘法对转子系统不对中量的大小进行计算. 主要结论如下：

1）利用动量矩定理，将圆盘不平衡力对转轴的影响等效到转子轴向力上，推导转子横向运动产生的轴向力，并建立考虑轴径向耦合效应的刚性双偏置圆盘转子-磁轴承系统动力学模型.

2）研究不同转速条件下磁轴承不对中对转子系统的影响，发现不对中会使得转子位移幅值变大，同时产生更高的控制电流，随着转速增加，不对中对转子系统的影响逐渐变高.

3）利用FFT将系统时域下的位移和电流响应转换成频域形式，并将其应用于系统不对中参数的识别算法中，通过频域下最小二乘法来计算不对中量的大小，在无噪声影响条件下，对比实际设定的不对中量，计算误差均在5.0%以内，在转子受到外界冲击力的影响时，该算法对不对中量的估计非常有效.

图 1 刚性双偏置圆盘转子-磁轴承系统模型

Figure 1. Model of a rigid double offset disc rotor-magnetic bearing system

下载: 全尺寸图片幻灯片

图 2 圆盘1形心和重心位置关系

Figure 2. Relationship between shape center and gravity center for disc 1

下载: 全尺寸图片幻灯片

图 3 动力学模型仿真

Figure 3. Dynamics model simulation

下载: 全尺寸图片幻灯片

图 4 转子受冲击后轴心轨迹

Figure 4. Axial trajectory of rotor subjected to impact

下载: 全尺寸图片幻灯片

图 5 不同转速下径向磁轴承处位移、电流稳态响应

Figure 5. Steady-state response of displacement and current at radial magnetic bearing at different rotational speeds

下载: 全尺寸图片幻灯片

图 6 不同转速下混合磁轴承处位移、电流稳态响应

Figure 6. Steady-state response of displacement and current at hybrid magnetic bearing at different rotational speeds

下载: 全尺寸图片幻灯片

图 7 不同噪声影响下计算的不对中量误差百分比

Figure 7. Percentage of misalignment error calculated under different noise interferences

下载: 全尺寸图片幻灯片

表 1 系统结构参数与仿真参数

Table 1. System structure parameters and simulation parameters

参数	数值	参数	数值
m/kg l/m	3.910 0.400	s₀/mm	0.400
m_d1/kg	1.065	δ_x1，δ_x2/mm	0.140，0.150
m_d2/kg	2.081	δ_y1，δ_y2/mm	0.160，0.145
e₁，e₂/μm	80，100	δ_z2，Δ_x1/mm	0.130，0.100
${e_{{\textit{z}}1}}，{e_{{\textit{z}}2}}$ /μm	8，10	Δ_y1，Δ_x2/mm	0.130，0.120
β₁，β₂/（°）	20，30	Δ_y2，Δ_z2/mm	0.110，0.150
I_d/（kg·m²）	0.0455	K_P/（A·m⁻¹）	5500
a₁/mm	0.226	K_I/（A·（m·s）⁻¹）	8000
a₂/mm	0.174	K_D/（A·s·m⁻¹）	3
b₁/mm	0.126	R t/ms	20 6
I_p1/（kg·m²）	0.0019	b₂/mm	0.074
I_p2/（kg·m²）	0.0059

下载: 导出CSV

表 2 频域下的径向磁轴承处转子位移和控制电流响应

Table 2. Rotor displacement and control current response at radial magnetic bearing in frequency domain

频率/ Hz	i	位移		电流
频率/ Hz	i	幅值/A	相位/（°）	幅值/m	相位/（°）
30	0	3.69 × 10⁻⁹	−16.09	0.929	−132.73
30	1	1.36 × 10⁻⁴	−32.43	0.754	152.99
35	0	1.13 × 10⁻⁹	112.79	0.930	−133.20
35	1	1.83 × 10⁻⁴	−61.97	1.016	124.49

下载: 导出CSV

表 3 无噪声影响下不对中量实际值和估计值比较

Table 3. Comparison of actual and estimated values of misalignment without noise interference

不对中量	实际值/mm	估计值/mm	误差/%
δ_x1	0.140	0.1388	−0.857
δ_y1	0.160	0.1575	−1.563
δ_x2	0.150	0.1439	−4.067
δ_y2	0.145	0.1456	0.414
δ_z2	0.130	0.1285	−1.154

下载: 导出CSV

参考文献(17)

[1]	陈宏,雷文平,陈磊,等. 一种转子动态不对中量计算方法[J]. 中国机械工程,2016,27(17): 2379-2383. doi: 10.3969/j.issn.1004-132X.2016.17.018 CHEN Hong, LEI Wenping, CHEN Lei, et al. A novel quantitive calculation method of dynamic misalignment in rotor systems[J]. China Mechanical Engineering, 2016, 27(17): 2379-2383. doi: 10.3969/j.issn.1004-132X.2016.17.018
[2]	甄满,田拥胜,孙涛,等. 具有不对中故障的双跨转子系统非线性动力学特性[J]. 机械工程学报,2020,56(16): 109-117. doi: 10.3901/JME.2020.16.109 ZHEN Man, TIAN Yongsheng, SUN Tao, et al. Nonlinear dynamics of two-span rotor-bearing system with flexible coupling misalignment[J]. Journal of Mechanical Engineering, 2020, 56(16): 109-117. doi: 10.3901/JME.2020.16.109
[3]	SHARI A, ALI A A, ALMUDHAFFER M. Combination of FFT & ICA methods for faults analysis of rotating machine[C]//Proceedings of the Information Conference on Information and Communication Technology. New York:[s.n.],2019:196-202.
[4]	KUPPA S K, LAL M. Dual flexible rotor system with active magnetic bearings for unbalance and coupling misalignment faults analysis[J]. Sādhanā,2019, 44:188.1-188.16.
[5]	LEES A W, SINHA J K, FRISWELL M I. Estimating rotor unbalance and misalignment from a single Run-down[J]. Materials Science Forum, 2003, 440/441: 229-236. doi: 10.4028/www.scientific.net/MSF.440-441.229
[6]	KUMAR P, TIWARI R. Development of a novel approach for quantitative estimation of rotor unbalance and misalignment in a rotor system levitated by active magnetic bearings[J]. Iranian Journal of Science and Technology, Transactions of Mechanical Engineering, 2021, 45: 769-786. doi: 10.1007/s40997-020-00364-7
[7]	KÄRKKÄINEN A, HELFERT M, AESCHLIMANN B, et al. Dynamic analysis of rotor system with misaligned retainer bearings[J]. Journal of Tribology-transactions of The Asme, 2008, 130: 319-320.
[8]	KUMAR P, TIWARI R. A numerical study on the effect of unbalance and misalignment fault parameters in a rigid rotor levitated by active magnetic bearings. [C]//ASME 2019 GasTurbine India Conference. [S.L.]:American Society of Mechanical Engineers,2019:1-7.
[9]	KUMAR P, TIWARI R. Dynamic analysis and identification of unbalance and misalignment in a rigid rotor with two offset discs levitated by active magnetic bearings: a novel trial misalignment approach[J]. Propulsion and Power Research, 2021, 10(1): 58-82. doi: 10.1016/j.jppr.2020.06.003
[10]	TIWARI R, KUMAR P. An innovative virtual trial misalignment approach for identification of unbalance, sensor and active magnetic bearing misalignment along with its stiffness parameters in a magnetically levitated flexible rotor system[J]. Mechanical Systems and Signal Processing, 2022, 167: 108540.1-108540.30.
[11]	汤华涛,吴新跃,刘海涛. 转子横向运动对轴向振动的影响分析[J]. 华中科技大学学报(自然科学版),2013,41(10): 74-78. TANG Huatao, WU Xinyue, LIU Haitao. Influence of a rotor′s lateral motion on its axial vibration[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2013, 41(10): 74-78.
[12]	朱熀秋,王绍帅. 六极径向–轴向主动磁轴承电磁特性分析及实验研究[J]. 中国电机工程学报,2020,40(5): 1653-1663. ZHU Huangqiu, WANG Shaoshuai. Electromagnetic characteristics analysis and experiment study of six-pole radial-axial active magnetic bearing[J]. Proceedings of the CSEE, 2020, 40(5): 1653-1663.
[13]	汪希平,朱礼进,于良,等. 主动磁轴承转子系统动力学特性的研究[J]. 机械工程学报,2001,37(11): 7-12. doi: 10.3321/j.issn:0577-6686.2001.11.002 WANG Xiping, ZHU Lijin, YU Liang, et al. Investigation on dynamic performance of active magnetic bearing rotor system[J]. Chinese Journal of Mechanical Engineering, 2001, 37(11): 7-12. doi: 10.3321/j.issn:0577-6686.2001.11.002
[14]	DARPE A, CHAWLA A, GUPTA K. Analysis of the response of a crachked Jeffcott rotor to axial excitation[J]. Jouranl of Sound and Vibration, 2002, 249(3): 429-445. doi: 10.1006/jsvi.2001.3870
[15]	邹博,徐园平,周瑾. 冲击激励下磁悬浮转子系统的响应特性分析[J]. 机械制造与自动化,2020,49(5): 41-43. ZOU Bo, XU Yuanping, ZHOU Jin. Analysis of response characteristics of magnetic suspension rotor system under shock excitation[J]. Machine Building & Automation, 2020, 49(5): 41-43.
[16]	陈伟,吴泽宇,韩佳奇,等. 突加基础冲击激励下转子系统振动特性试验[J]. 航空动力学报,2023,38(4): 878-888. CHEN Wei, WU Zeyu, HAN Jiaqi, et al. Test on vibration characteristics of rotor system under sudden base shock excitation[J]. Journal of Aerospace Power, 2023, 38(4): 878-888.
[17]	SINGH S, TIWARI R. Model-based fatigue crack identification in rotors integrated with active magnetic bearings[J]. Journal of Vibration and Control, 2017, 23(6): 980-1000. doi: 10.1177/1077546315587146