Modeling and Dynamics Analysis of High-Temperature Magnetic Bearing-Rotor System
-
摘要:
在多电航空发动机中,主动磁悬浮轴承因其耐高温、非接触等特性可以突破温度对支承部位的限制,使支承部位能够更靠近燃烧室. 为探究温度对磁悬浮轴承转子系统动态特性的影响规律,提出一种高温磁悬浮轴承转子系统动力学建模方法. 通过仿真得到转子在不同温度下的温度分布,并使用多项式拟合转子轴向温度分布;基于有限元方法推导柔性转子单元的动力学模型,引入温度影响,建立考虑温度影响的磁悬浮轴承转子系统整体动力学模型,并通过模态试验验证模型的准确性;基于理论动力学模型分析系统的动态特性. 结果表明:温度升高会导致转子的前三阶支承模态频率下降,增大各阶幅频响应幅值;当温度从常温升至450 ℃时,转子的前三阶弯曲支承模态频率分别降低了3.818%、5.670%、3.183%,前三阶弯曲模态幅频响应幅值分别升高了83.4%、34.4%、24.1%.
Abstract:In the multi-electric aircraft engine, an active magnetic bearing can break through the limitation of temperature on the support part due to its high temperature resistance and non-contact characteristics, which enables its support part to be closer to the combustion chamber. In order to investigate the influence of temperature on the dynamic characteristics of the magnetic bearing-rotor system, a dynamics modeling method for a high-temperature magnetic bearing-rotor system was proposed. The temperature distributions of the rotor at different temperatures were obtained through simulation, and the axial temperature distribution of the rotor was fitted using polynomials. Based on the finite element method, the dynamics model of the flexible rotor unit was derived. The temperature influence was introduced, and the overall dynamics model of the magnetic bearing- rotor system considering the temperature influence was established. The accuracy of the model was verified by a modal test. The dynamic characteristics of the system were analyzed based on the theoretical dynamics model, and the results show that an increase in temperature leads to a decrease in the first three orders of the support modal frequency of the rotor and an increase in the amplitude of the amplitude frequency response of each order. When the temperature increases from room temperature to 450 ℃, the first three orders of the bending support modal frequency of the rotor decrease by 3.818%, 5.670%, and 3.183%, respectively, and the amplitudes of the first three orders of the bending modal amplitude frequency response increase by 83.4%, 34.4%, and 24.1%, respectively.
-
在长期服役过程中,杂质颗粒会逐渐侵入道砟骨架空隙,造成道床脏污. 道砟脏污不仅会产生高昂的维护成本,还会增加铁路运行的安全风险. 研究表明,脏污会降低道床稳定性、承载力和排水能力,并引发道床板结和翻浆冒泥等病害[1-2]. 脏污质主要包括:1) 基床表层及底砟层渗入的细小土石颗粒;2) 上部侵入的煤渣、矿渣、砂石颗粒以及枕木风化后的碎屑;3) 发生磨损与破碎的道砟颗粒[3]. 目前,国内外关于脏污道砟的力学行为特性已有大量研究. Koohmishi等[4]研究了脏污道砟颗粒级配及尺寸对其导水性能的影响,结果显示:脏污道砟的排水性能受道砟级配、污染物级配及污染程度的影响;受细砂污染道砟的渗透系数比受粗砂污染道砟的渗透系数低10倍左右. Kashiani等[5]针对磨损的脏污道砟,开展了一系列大型三轴压缩试验,研究了颗粒含水率对其力学特性的影响. Indraratna等[6-7]利用室内三轴试验,并结合相应离散元(DEM)数值模型,通过分析颗粒运动及颗粒间接触特征,提出了污染物对道砟的“润滑剂”效应,认为脏污质填充于道砟空隙,弱化了颗粒间的咬合,从而降低了集料的抗剪强度.
为提升有砟轨道道床性能,在铁道工程领域常常采用格栅、土工织物等材料. 研究表明,格栅可显著提高道砟集料的抗剪强度,并有效减小道床沉降[8–10]. 土工格栅主要通过约束道砟集料的侧向变形,增加道砟颗粒间的咬合,从而实现其加固机理[11-12]. 然而,由于脏污质填充于集料骨架空隙间(如煤渣、土石颗粒等),或附着在道砟颗粒表面(如湿润的黏土质等),道砟颗粒间以及道砟颗粒与格栅纵横肋之间的嵌固咬合作用将会发生改变.
目前,国内外对脏污道砟或格栅加筋道砟的力学性能已进行了充分研究,但针对不同脏污质对格栅加固性能影响进行比较研究的文献较少. Indraratna等[13-14]分析了不同脏污程度下格栅加筋道砟的强度及变形特性,但脏污质仅局限于煤渣颗粒,不同脏污情况下格栅加筋道砟的力学性能演化规律仍不明晰. 根据Huang等[15-16]关于不同污染物对道砟的强度及变形特性的研究,污染物的种类及含量对道砟力学行为具有明显影响. 因我国现有的有砟铁路列车运载性质多样、所处地理环境复杂等原因,道床脏污来源不尽相同. 如处于华北地区的大秦线、蒙华铁路等多条货运线路主要受到上部的煤渣污染;对于南方山区、丘陵地带的铁路轨道,由于地基层不充分排水,产生的黏土污染最为严重;煤渣及黏土的材料物理性质、颗粒级配及强度特性等具有明显差异. 因此,本文采用煤渣、黏土作为脏污质,开展了不同脏污程度下格栅加筋道砟集料的大型直剪试验,分析不同法向压力条件下各试样强度及形变发展规律,探讨了其抗剪强度、峰值摩擦角、剪胀角等力学参数的变化特性,分析比较了不同脏污质对格栅加筋道砟性能影响的差异,以期为实际铁路运维过程中格栅加固的方法设计及效果评判提供理论依据.
1. 试验设计
1.1 试验仪器
由于道砟颗粒尺寸较大,传统直剪设备无法满足试验要求,因此本文自主设计一套大型道砟集料直剪设备,如图1所示. 该设备由4部分构成:1) 反力系统,包括反力架(构件)及反力底座(构件);2) 加载系统,包括水平千斤顶及竖向千斤顶,该设备使用的千斤顶可提供15 t推力,满足该直剪试验所需的要求,在剪切过程中竖向千斤顶施加的法向压力恒定,而水平千斤顶由数控油压设备控制,实现匀速水平推力;3) 装样系统,包括上、下剪切盒及盖板;4) 数据采集系统,包括位移传感器、静态应变仪.
需要特别指出的是,已有的大型直剪设备多采用推动上剪切盒方式,而研发的设备采用水平千斤顶(构件)推动下剪切盒(构件)方式施加水平推力,该方式可避免在剪切过程中由竖向千斤顶(构件)施加的法向压力产生偏心. 除此之外,为保证剪切过程中剪切面积恒定,下剪切盒剪切方向尺寸较上剪切盒大100 mm.
1.2 试样制备
该试验使用的道砟级配满足我国《铁路碎石道砟》(TB/T 2140—2008)[17]标准,相关参数见图2.
从图2可知:道砟最大粒径Dmax = 64 mm,有效粒径D10 = 24 mm,中值粒径D30 = 34 mm,控制粒径D60 = 42 mm,平均粒径D50 = 39 mm,曲率系数Cc = 1.75 mm,不均匀系数Cu = 1.15 mm. 分别采用煤渣、黏土作为污染杂质,其相关力学指标见表1,所用格栅如图3所示. 根据Mcdowell等[18]的研究,最优格栅孔径与道砟尺寸之比为1.4,即格栅孔径为1.4倍道砟尺寸可提供最优约束. 此试验所用道砟D50 = 39 mm,故孔径尺寸选为55 mm,格栅相关参数见表2. 根据Indraratna等[19]对道砟脏污程度的定义,采用VCI (void contamination index, VVCI)作为该试验脏污质污染指标,见式(1).
VVCI=1+efebGs,bGs,fMfMb×100%, (1) 式中:ef、 Gs,f、 Mf和eb、 Gs,b、 Mb分别为脏污质和道砟集料的孔隙比、比重、干重.
Indraratna等[14]对煤渣脏污下格栅加筋道砟力学特性的研究显示,当VCI大于40%时,脏污集料的归一化峰值剪应力相较于干净集料降低50%以上. 在实际服役中,当道砟集料的抗剪强度由于外部污染或颗粒劣化而明显降低时,需采用重新夯实或材料换新等措施,故本文仅对VCI小于40%的脏污情况进行了研究. 该试验分别制备煤渣、黏土脏污条件下VCI为0、20%、40%的直剪试样各4组,为反映脏污质存在条件下格栅加固性能的变化,本研究设置了相同条件下不含格栅的各组试样进行对比,试样具体设置情况见表3.
表 1 道砟及脏污质相关力学指标Table 1. Mechanical indexes of ballast and fouling substances名称 比重 堆积密度/
(kg•m−3)含水
率/%孔隙
比 e液限ωL/% 塑限ωP/% 道砟 2.66 1432 0.858 黏土 2.70 1178 22.2 1.207 42.10 22.40 煤渣 1.22 715 0.707 表 2 格栅参数Table 2. Physical and technical properties of geogrid格栅参数 数值 孔径形式 双向孔径 55 mm × 55 mm 材料 聚丙烯 2% 应变时抗拉强度/(kN•m−1) 11 5% 应变时抗拉强度/(kN•m−1) 15 峰值抗拉强度/(kN•m−1) 30 屈服点伸长率/% 13 表 3 试样设置情况Table 3. Details of experimental specimens类型 试样编号 脏污质种类 VCI/% 法向压力 σn/kPa 含格栅 GN1 黏土 0 15、35、55、75 GN2 20 GN3 40 GM1 煤渣 0 15、35、55、75 GM2 20 GM3 40 不含格栅 N1 黏土 0 15、35、55、75 N2 20 N3 40 M1 煤渣 0 15、35、55、75 M2 20 M3 40 1.3 试验步骤
为减小试验误差,试验前开展了预试验并以此确定直剪试样中道砟集料的堆积密度(压实后为1432 kg/m3),从而计算各组试验中道砟的合理用量. 为减小剪切过程中下剪切盒与底板以及上下剪切盒之间的摩擦力,在底板上布置10根直径为10 mm的圆棒(如图1),并在上下盒剪切面及内部均涂抹润滑油. 制样时,每组试样称取等量道砟集料275 kg,与式(1)计算得到相应质量的脏污杂质充分拌合. 装样过程中,每150 mm高度采用低锤重压的方式对试样进行压实. 格栅布置在剪切面,即距下剪切盒底面高250 mm处.
每组试样加载前,对试样施加100 kPa的法向压力进行过载压实,该压力持续60 s,待试样稳定后缓慢卸载至0,随后缓慢施加法向压力至设定值,即15、35、55、75 kPa. 剪切过程中,水平千斤顶以3 mm/min的速度缓慢施加水平推力;法向压力通过伺服控制系统维持压力恒定. 通过预试验发现,在剪切应变达到10%之前,剪切力已逐渐趋于稳定,剪切面出现明显破坏. 因此本试验在加载过程中待剪切应变达到10%,即剪切位移达到60 mm时停止.
2. 试验结果
2.1 强度特性分析
黏土及煤渣分别作为脏污质时,在不同污染程度下道砟试样的剪切应力τs随剪切位移 Δs发展情况如图4. 从图可以看出:随着法向压力σn的提高,各试样的抗剪强度均逐渐提高;相同试验环境下(即污染物种类、VCI值及法向压力均维持一致),不论脏污质材料及脏污质含量,受格栅加固道砟集料的抗剪强度明显高于未受格栅加固道砟集料的抗剪强度,格栅的存在可显著提高道砟集料的抗剪强度;在道床受到外界杂质污染时,布置格栅可作为减轻道砟抗剪强度降低的手段,这是由于格栅的纵横肋限制了道砟颗粒的平移与旋转,从而加强了道砟颗粒之间的咬合作用[11-12].
为探究不同脏污质对格栅加固效果的影响,图5展示了不同污染程度下,格栅加固道砟集料分别受黏土和煤渣污染时,峰值抗剪强度τsmax与σn关系. 从图5可以看出:干净的受格栅加固道砟集料表现出最大峰值抗剪强度;同一法向压力下,随着VCI从20%增加至40%,格栅加固道砟的峰值抗剪强度逐渐降低;受黏土污染或受煤渣污染的格栅加固道砟在不同法向压力条件下具有不同的响应,在σn=15 kPa时,VCI从20%增加至40%,受黏土污染的试样峰值强度降低20.3%,随着σn增加至75 kPa时,该试样峰值强度降低9.1%;对于受煤渣污染的试样,峰值强度降低比受σn的影响较小,始终维持在7.3%左右.
值得注意的是,同一污染程度下,受煤渣污染的格栅加固道砟展现出更低的抗剪强度. 这说明相较于黏土脏污,煤渣会对格栅加固道砟的剪切性能产生更不利的影响,这是由黏土和煤渣不同的颗粒粒径所导致,如图6所示. 由图可知:当污染程度较低时,由于黏土的粒径较小,其附着在道砟颗粒表面,道砟颗粒仍充分接触,而煤渣颗粒则填充于部分道砟骨架中间,阻碍了道砟颗粒的接触,导致道砟颗粒间力链的断裂;当污染程度较高时,黏土完全填充道砟骨架孔隙,或格栅与道砟颗粒之间,限制了道砟的运动;相反而言,由于煤渣颗粒粒径较大,其之间的相对运动导致道砟颗粒更易移动与旋转. 当VCI值达到40%时,受黏土污染和受煤渣污染的格栅加固道砟峰值剪切强度最低,且两者抗剪强度几乎相等,此时已达到道床需要维护的临界点.
图7为各污染情况下道砟集料的峰值摩擦角Φmax与σn关系曲线. 对于同种污染物,在相同法向压力条件下,随着污染程度的增加,格栅加固道砟的峰值摩擦角会逐渐减小. 而对比不同脏污质影响下的道砟集料,可以发现煤渣污染使格栅加固道砟展现出更低的峰值摩擦角. 这说明,相较于黏土脏污,煤渣脏污对格栅加固集料抗剪性能具有更不利的影响.
2.2 变形特性分析
为研究脏污情况下道砟集料的变形特性,图8展示了受黏土或煤渣污染时,道砟试样的法向位移Δn随Δs变化曲线. 从图中可以看出:所有试样在剪切前期均经历了短暂的剪缩阶段,且随着法向压力的增加,剪缩现象更明显,剪缩阶段持续的时间越长;随后试样体积逐渐增加,出现剪胀现象直到整个剪切过程结束. 相同压力条件下,受黏土或煤渣污染的道砟试样法向位移随剪切位移变化如图9所示. 为节约篇幅,此处仅展示了法向压力为55 kPa时的各试样法向位移变化图. 对比VCI为0 (即未受污染的格栅加固道砟试样)与VCI为20%和40%的格栅加固道砟试样,可以发现:干净道砟的体积变形最大;随着道砟脏污程度增加,试样的体积变形会减小. 这是因为当污染程度逐渐增加时,更多的道砟骨架孔隙被脏污质颗粒所填充,道砟颗粒的相对运动更困难. 尽管格栅可限制道砟的侧向变形,但剪切过程中道砟颗粒的重排布受到了脏污质颗粒的限制,故造成试样体积变化更不明显.
对比同一污染程度下受黏土和煤渣污染的格栅加固道砟试样,受煤渣污染试样的体积变形比受黏土污染试样的大30.5%. 这是由于脏污杂质材料性质不同导致:相较于煤渣颗粒粗糙的表面,黏土颗粒的表面相对较光滑,故黏土颗粒与道砟颗粒之间的摩擦力较小,道砟颗粒更容易重新排布,导致体积变形更小.
剪胀角常用来表示材料在剪切过程中体积变化率. 图10为不同压力条件下,分别受煤渣、黏土脏污质污染的格栅加固道砟峰值剪胀角ϕmax与σn关系. 由图可以看出:对于所有试样,法向压力的增加会导致峰值剪胀角的降低;不受格栅加固的干净道砟ϕmax值最大,格栅加固干净道砟次之;随着污染程度的增加,试样的ϕmax值也逐渐减小;同时,受煤渣污染的格栅加固道砟ϕmax值比受黏土污染的格栅加固道砟大. 这说明,在相同法向压力条件下,发生相同剪切位移时,前者比后者更容易发生体积膨胀.
3. 结 论
为研究不同脏污杂质对格栅加固性能的影响,本文选取煤渣和黏土两种道床常见脏污材料,对格栅加固道砟以及未受格栅加固道砟开展了一系列不同污染程度下的直剪试验,研究了脏污质对格栅加固机理的影响,并分析了不同脏污质的影响差异机理,可得到以下结论:
1) 随着污染程度的增加,煤渣或黏土质污染均会导致道砟集料剪应力、峰值摩擦角降低,法向剪胀位移、峰值剪胀角减小.
2) 同一脏污程度下,煤渣脏污引起的格栅加固道砟剪切强度下降现象更显著. 当污染指标VCI从20%增加至40%时,在低围压下黏土污染的格栅加固道砟强度降低20.3%,在高围压下该值降低9.1%;而煤渣污染的格栅加固道砟受法向压力变化影响较小,始终为7.3%左右. 这是由于煤渣颗粒粒径较黏土颗粒粒径更大所导致.
3) 对于脏污质对格栅加固道砟变形特性的影响,煤渣污染的格栅加固道砟具有更明显的剪缩现象,峰值剪胀角也更大. 同一污染程度下,煤渣污染的格栅加固道砟体积变形比黏土污染的试样大30.5%. 这是由于煤渣颗粒表面更粗糙,导致道砟颗粒在受到剪切时更难发生重排布,从而引起更大的体积变形.
4) 相较于黏土脏污,煤渣脏污会对格栅加固道砟的剪切力学行为产生更不利的影响. 在实际工程中应该更注意道床上部列车产生的煤渣污染.
-
表 1 机械系统所使用的材料
Table 1. Materials used in mechanical systems
序号 所用材料 1,2,3,4,6,8 1Cr11Ni2W2MoV 5 GH4169 7,10 高温软磁合金(1J22) 9,11 38 黄铜 表 2 数学模型拟合的R2值
Table 2. R2 by mathematical model fitting
m 150 ℃ 300 ℃ 450 ℃ 1 0.799 123 0.815 718 0.824 486 2 0.993 046 0.993 757 0.993 782 3 0.996 804 0.995 234 0.994 565 4 0.998 708 0.998 45 0.998 344 5 0.999 406 0.999 188 0.999 120 6 0.999 495 0.999 374 0.999 328 表 3 磁悬浮轴承的结构参数和电控参数
Table 3. Structural parameters and electrical control parameters of magnetic bearing
参数 取值 参数 取值 转子半径/mm 89.4 比例系数 kp 1 磁极面积 As/mm2 200 积分系数 ki 1 磁极夹角 θ/(°) 22.5 微分系数 kd 0.0 005 单边气隙 sa/mm 0.3 偏置电流 I/A 2 定子齿数 8 线圈匝数 N/匝 240 表 4 支承模态频率对比
Table 4. Comparison of support modal frequency
温度/℃ 阶次 计算频率/Hz 试验频率/Hz 误差/% 25 1 994.490 994.375 0.0 115 2 1506.200 1506.250 0.0 033 3 2437.090 2494.380 2.2 960 150 1 990.606 983.750 0.6 900 2 1471.570 1496.880 0.0 169 3 2421.910 2477.610 2.3 000 300 1 982.417 976.250 0.6 300 2 1450.280 1450.990 0.0 493 3 2404.800 2460.590 2.2 320 450 1 968.743 958.750 1.0 000 2 1412.160 1413.460 0.0 922 3 2374.860 2391.080 2.6 830 -
[1] 吴志琨,李军,时瑞军. 多电航空发动机研究现况及关键技术[J]. 航空工程进展,2012,3(4): 463-467. doi: 10.3969/j.issn.1674-8190.2012.04.015WU Zhikun, LI Jun, SHI Ruijun. Current research status and key technologies of more-electric aeroengine[J]. Advances in Aeronautical Science and Engineering, 2012, 3(4): 463-467. doi: 10.3969/j.issn.1674-8190.2012.04.015 [2] 刘程子,湛江,杨艳,等. 主动磁悬浮轴承–柔性转子的研究和发展综述[J]. 中国电机工程学报,2020,40(14): 4602-4614,4739.LIU Chengzi, ZHAN Jiang, YANG Yan, et al. Review of research status and development of flexible rotor-magnetic bearing[J]. Proceedings of the CSEE, 2020, 40(14): 4602-4614,4739. [3] 金超武. 高温磁悬浮轴承若干关键技术研究[D]. 南京:南京航空航天大学,2011. [4] ZIAEI-RAD S, KOUCHAKI E, IMREGUN M. Thermoelastic mdeling of rotor response with shaft rub[J]. Journal of Applied Mechanics, 2010, 77(6): 061010.1-061010.1. [5] 朱向哲,贺威,袁惠群. 稳态温度场对转子系统振动特性的影响[J]. 东北大学学报(自然科学版),2008,29(1): 113-116.ZHU Xiangzhe, HE Wei, YUAN Huiqun. Effects of steady temperature field on vibrational characteristics of a rotor system[J]. Journal of Northeastern University (Natural Science), 2008, 29(1): 113-116. [6] 何鹏,刘占生,刘镇星. 考虑杨氏模量随轴向温度分布变化的转子有限元建模方法研究[J]. 振动与冲击,2012,31(14): 22-26,55. doi: 10.3969/j.issn.1000-3835.2012.14.005HE Peng, LIU Zhansheng, LIU Zhenxing. Finite element modelling of rotor considering the variation of Yang’s modulus with axial temperature distribution[J]. Journal of Vibration and Shock, 2012, 31(14): 22-26,55. doi: 10.3969/j.issn.1000-3835.2012.14.005 [7] GOU X F, ZHU L Y, QI C J. Nonlinear dynamic model of a gear-rotor-bearing system considering the flash temperature[J]. Journal of Sound and Vibration, 2017, 410: 187-208. [8] XIANG F G, CHEN X, ZHANG B, et al. Simulation method of rotor dynamic characteristics considering temperature distribution and aerodynamic load[C]//Proceedings of ASME Turbo Expo 2022:Turbomachinery Technical Conference and Exposition. Rotterdam:ASME, 2022:13-17. [9] DU T C, SUN Y H, GENG H P, et al. Dynamic Analysis on Rotor System Supported by Active Magnetic Bearings based on Sliding Mode Control[C]//2018 IEEE International Conference on Mechatronics and Automation (ICMA). Changchun:IEEE, 2018:1960-1965. [10] MUSHI S E, LIN Z L, ALLAIRE P E. Design, construction, and modeling of a flexible rotor active magnetic bearing test rig[J]. ASME Transactions on Mechatronics, 2012, 17(6): 1170-1182. [11] XU Y P, ZHOU J, LIN Z L, et al. Identification of dynamic parameters of active magnetic bearings in a flexible rotor system considering residual unbalances[J]. Mechatronics, 2018, 49: 46-55. [12] 沈权,周瑾,马彦超,等. 基础激励下磁悬浮转子系统动力学建模与分析[J]. 振动与冲击,2022,41(17): 35-47,72.SHEN Quan, ZHOU Jin, MA Yanchao, et al. Dynamic modeling and analysis of maglev rotor system under base excitation[J]. Journal of Vibration and Shock, 2022, 41(17): 35-47,72. [13] WU C, SU Z Z, WANG D, et al. Dynamic modeling method for active magnetic bearings-rotor system of steam turbines[J]. Journal of Mechanical Science and Technology, 2023, 37(4): 1665-1673. [14] 王戈一. 磁悬浮多电发动机的研究[J]. 燃气涡轮试验与研究,2007,20(4): 15-18,35. doi: 10.3969/j.issn.1672-2620.2007.04.003WANG Geyi. Study of a more-electric engine with active magnetic bearings[J]. Gas Turbine Experiment and Research, 2007, 20(4): 15-18,35. doi: 10.3969/j.issn.1672-2620.2007.04.003 [15] BURDET L. Active magnetic bearing design and characterization for high temperature applications[R]. Swiss:Federal Institute of Technology in Lausanne, 2006. [16] MONTAGUE G T, JANSEN M J, PROVENZA A, et al. Experimental high temperature characterization of a magnetic bearing for turbomachinery[C]//59th Annual Forum and Technology Display. Phoenix:Vertical Flight Society, 2003:1-18. [17] SHI Z, SUN X D, LEI G, et al. Multiobjective optimization of a five-phase bearingless permanent magnet motor considering winding area[J]. IEEE/ASME Transactions on Mechatronics, 2022, 27(5): 2657-2666. [18] 周扬,周瑾,张越,等. 基于RBF近似模型的磁悬浮轴承结构优化设计[J]. 西南交通大学学报,2022,57(3): 682-692. doi: 10.3969/j.issn.0258-2724.20210766ZHOU Yang, ZHOU Jin, ZHANG Yue, et al. Optimum structural design of active magnetic bearing based on RBF approximation model[J]. Journal of Southwest Jiaotong University, 2022, 57(3): 682-692. doi: 10.3969/j.issn.0258-2724.20210766 [19] 徐园平. 柔性转子磁悬浮轴承支承特性辨识[D]. 南京:南京航空航天大学,2018. -