[发明专利]一种微型气体轴承静动态特性及稳定性测试方法

摘要

本发明涉及微型气体轴承测试技术领域，公开了一种微型气体轴承静动态特性及稳定性测试方法，包括以下步骤：1)试验前标定测试轴承转子的初始位置；2)连续采集测试轴承转子的试验数据；3)计算获得测试轴承转子的静动态特性参数和图表，并判断测试轴承转子的稳定性；将所述试验数据代入MATLAB的模型方程中进行计算求解，得到测试轴承转子在不同工作状态下的刚度和阻尼，根据刚度和阻尼构建稳定性特征方程，根据稳定性特征方程及轴承稳定性理论判断测试轴承转子的稳定性，这种微型气体轴承稳定性测试及控制方法操作简单、计算量小，可以连续测量并计算出测试轴承转子在不同工作状态下的静、动态特性参数，并判断测试轴承转子的稳定性。

主权项

1.一种微型气体轴承静动态特性及稳定性测试方法，其特征在于，包括以下步骤：1)试验前标定测试轴承转子的初始位置开启激光位移传感器，此时测试轴承转子与轴承座外圈在最低点接触，Y轴激光位移传感器记录此时Y轴初始最大值y0max,将测试轴承转子沿Y轴提升至与轴承座外圈在最高点接触，Y轴激光位移传感器记录此时Y轴初始的最小值y0min,然后将测试轴承转子Y轴的初始值设置为y0＝y0max‑1/2(y0min+y0max),X、Z轴的初始值设置为0；2)连续采集测试轴承转子的试验数据启动测试轴承转子的压力控制阀，确保测试轴承转子处于悬浮状态时开启电磁驱动系统驱动测试轴承转子，试验机连续采集测试轴承转子在不同工作状态下分别进行空载和激振运行时的试验数据，并将不同工作状态下的试验数据送入信号采集分析系统和MATLAB的模型方程中进行分析计算；3)计算获得测试轴承转子的静动态特性参数和图表，并判断测试轴承转子的稳定性信号采集分析系统将所述试验数据进行分析计算，直接生成测试轴承转子在不同工作状态下静、动态特性参数和图表包括：转速、偏心量、振动波形图、频谱图、轨迹图、振幅‑时间‑频率三维图和振幅‑转速分叉图；将测试轴承转子在不同工作状态下静、动态特性参数和所述试验数据代入MATLAB的模型方程中进行计算求解，得到测试轴承转子在不同工作状态下的刚度和阻尼，根据刚度和阻尼构建稳定性特征方程，根据稳定性特征方程判断测试轴承转子的稳定性；所述步骤3)中的稳定性特征方程是根据MATLAB的模型方程轴承‑转子系统动力学方程和气膜增量与运动参数关系方程获得的；所述轴承‑转子系统动力学方程为：所述气膜增量与运动参数关系方程为：联立方程(1)和(2)得模型方程(3)为：模型方程(3)变换为模型矩阵方程组(4)为：模型矩阵方程组(4)中参数的计算公式(5)为：其中，X0、Y0、Z0为测试轴承转子轴心轨迹的初始位置；Xn、Yn、Zn为tn时刻测试轴承转子轴心轨迹的位置；ΔXn、ΔYn、ΔZn为X、Y、Z轴上激光位移传感器tn时刻测试的振动量；m为测试轴承转子的质量；g为重力加速度；Δt为两次采样的时间间隔；en为tn时刻测试轴承转子的偏心量；Δen为tn时刻测试轴承转子的偏心量变化量；θn为tn时刻测试轴承转子的偏心角；Δθn为tn时刻测试轴承转子的偏心角变化量；为t_n时刻与参数e、θ、Z相关的气膜增量；为t_n时刻与参数e、θ、Z有关的主刚度；为t_n时刻与参数e、θ、Z两两耦合有关的交叉刚度；为t_n时刻与参数e、θ、Z有关的主阻尼；为t_n时刻与参数e、θ、Z两量耦合有关的的交叉阻尼；为t_n时刻测试轴承转子的偏心量变化量一次导数；为t_n时刻测试轴承转子的偏心角变化量一次导数；为t_n时刻测试轴承转子的X轴向位移变化量一次导数；为t_n时刻测试轴承转子的偏心量二次导数；为t_n时刻测试轴承转子的偏心角二次导数；为t_n时刻测试轴承转子的X轴向位移二次导数；将试验数据采集中第n‑2到n+3次采样所对应的测试轴承转子的参数ΔX_n、ΔY_n、ΔZ_n、Δe_n、Δθ_n、代入模型矩阵方程组(4)中进行求解，n≥2，得到第n次采样时刻测试轴承转子的18个刚度和阻尼值矩阵为：将式(6)中的刚度和阻尼值代入下列各式计算测试轴承转子稳定性特征方程的系数：a0、a1、a2、a3、a4、a5和a6a0＝m3a1＝m2(bee+bθθ+bzz)a2＝m2(kee+kθθ+kzz)+m(beebθθ+bθθbzz+bzzbee‑bezbze‑beθbθe‑bθzbzθ)a3＝m[kee(bθθ+bzz)+kθθ(bee+bzz)+kzz(bee+bθθ)]‑m[(kezbze+kzebez)+(keθbθe+kθebeθ)+(kθzbzθ+kzθbθz)]+(bzzbθθbee+bθebezbzθ+bzebeθbθz)‑(bezbzebθθ+bθzbzθbee+bθebeθbzz)a4＝m[(keekθθ+kθθkzz+kzzkee)‑(kezkze+kθzkzθ+kθekeθ)]+kee(bθθbzz‑bzθbθz)+kθθ(bzzbee‑bzebez)+kzz(beebθθ‑beθbθe)+keθ(bθzbze‑bθebzz)+kθz(bzebeθ‑bzθbee)+kze(beθbθz‑bezbθθ)+kθe(bzθbez‑beθbzz)+kzθ(bezbθe‑bθzbee)+kez(bθebzθ‑bzebθθ)a5＝bee(kθθkzz‑kθzkzθ)+bθθ(kzzkee‑kzekez)+bzz(keekθθ‑keθkθe)+bez(kθekzθ‑kzekθθ)+bzθ(kezkθe‑kθzkee)+bθe(kzθkez‑keθkzz)+beθ(kzekθz‑kθekzz)+bθz(keθkze‑kzθkee)+bze(kθzkeθ‑kezkθθ)a6＝(kzzkθθkee+kθzkzekeθ+kezkzθkθe)‑(kezkzekθθ+kθzkzθkee+keθkθekzz)m为测试轴承转子的质量；kee、kθθ、kzz为与偏心量e、偏心角θ和轴向振动z有关的主刚度；keθ、kez、kθe、kθz、kze、kzθ为偏心量e、偏心角θ和轴向振动z两两耦合有关的交叉刚度；bee、bθθ、bzz为与偏心量e、偏心角θ和轴向振动z有关的主阻尼；beθ、bez、bθe、bθz、bze、bzθ为偏心量e、偏心角θ和轴向振动z两两耦合有关的交叉阻尼；根据稳定性特征方程的系数a0、a1、a2、a3、a4、a5和a6的结果构造测试轴承转子的稳定性特征方程(7)为：a0v6+a1v5+a2v4+a3v3+a4v2+a5v1+a6＝0(7)依据轴承稳定性理论，轴承‑转子系统的稳定性取决特征根v在复平面的分布状况：当v<0时，系统将处于稳定状态；当v＝0时，系统处于临界状态；当v>0时，系统处于失稳状态；判断轴承‑转子系统在某一工作状态下是否稳定，只需判定特征根v是否具有负实部。