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]对采用弹性轨枕与道砟垫的有砟轨道进行了疲劳试验,发现该种轨道结构能够保持良好的几何形位及受力状态;刘晓春等[3]对新型钢管混凝土(CFT)双块式轨枕进行了疲劳试验和静载试验,结果表明,双块式轨枕在疲劳试验前后均具有良好的工作性能.
国内外众多学者对优化后轨枕的横向阻力进行了研究:Hayano等[4-6]尝试通过在普通条形混凝土轨枕的侧面或底面布置混凝土条的方法提高轨枕的性能,试验结果表明这类方法能有效提高轨枕横向阻力;郭云龙等[7]为合理设计摩擦型轨枕底部凹槽,提出了3种不同箭头型凹槽并进行了横向阻力试验,结果表明,各类型的摩擦型轨枕均可增大道床横向阻力,最少提升横向阻力7%,最高提升21%. 目前应用钢枕的轨道道床横向阻力明显不足,Jing等[8]和Zakeri等[9]对钢枕进行结构优化并开展横向阻力试验,发现钢枕添加腹板加劲肋或使用枕下橡胶垫可以不同程度地提升道床横向阻力;井国庆等[10-11]对梯形轨枕和双块式轨枕与传统条形轨枕做了横向阻力试验,结果表明,梯形轨枕与双块式轨枕的横向阻力分别比条形轨枕高出55%和4%;Liegner等[12]设计了一种Y型钢枕,并研究了该轨枕在机车运动下的横向位移,结果表明,采用Y型钢枕木施工的轨道运行性能良好,但由于成本较高,在国内外均未大范围应用.
综上所述,目前对轨枕优化的研究主要针对传统条形轨枕进行局部优化,但轨枕整体结构异曲同工. 此外,马春生[13]发现道床纵、横向阻力之间呈现了较强的正相关关系,而道床阻力与支承刚度间的相关性并不明显,建议采用道床横向阻力指标代替道床纵、横向阻力2项指标,以优化有砟道床评价体系. 因此,本文针对自主知识产权的X异形结构轨枕[14-16],开展X形轨枕道床竖向静载和横向阻力试验,并与传统Ⅲ型轨枕结果进行对比分析. 通过离散元仿真模拟研究,从细观上揭示X形轨枕与有砟道床的相互作用机制,为实际工程中新型X形轨枕的应用提供相应的试验支撑.
1. 缩尺道床模型试验
1.1 3D打印轨枕和缩尺道床模型
由于全尺寸有砟轨道试验耗时长、成本和误差高,因此,本研究对X形轨枕和Ⅲ型轨枕道床进行缩尺模型试验. 国内外学者从理论和试验的角度对轨道缩尺模型试验方法进行了详细研究:万轶等[17]设计了有砟轨道道床纵向阻力试验缩尺模型并证明其满足道床纵向阻力的相似要求;Shear等[18]建立1/3缩尺车轮-环轨模型,研究了列车移动荷载作用下有砟轨道的沉降规律以及冲击载荷工况下轮轨的接触响应;邹锦华等[19]通过1∶20的车辆-轨道缩尺模型振动试验进行动力学性能分析,为我国进一步开展轮轨系统的缩尺动力模型试验提供参考. 类似的,朱建华[20]发现,相似模型在静力与动力荷载作用下能够准确反演出原型结果. 以上研究均证明了缩尺轨道模型试验的可行性.
试验所用的Ⅲ型轨枕和X形抗滑动轨枕模型均通过3D打印技术使用聚乳酸(PLA)材料打印,如图1所示. PLA物理力学性能参数如表1所示. 考虑所用3D打印机的最大尺寸限制,本试验所用条形轨枕模型长度为300 mm,弹性模量为
3500 MPa,长度和弹性模量分别相对真实轨枕缩小8.7倍和10.3倍,道砟材料与实际材料密度相同. 基于相似理论确定几何相似常数CL=8.7,弹性模量相似常数CE=10.3,道砟密度相似常数Cρ=1.0,试验相似参数根据量纲矩阵分析法推导如表2所示.根据相似几何常数可以推导出缩尺模型的各部分结构尺寸. 根据《铁路碎石道砟》标准[21],按相似比得到的特级道砟级配曲线如图2所示. 道砟材料为特级花岗岩,在铺设过程中为保证堆砟密实,铺设时采用分层夯实法. 在相同的轨距和轨枕扣件间距情况下,4根Ⅲ型条形轨枕与双根X形轨枕有砟轨道铺设效果图3所示.
表 1 3D打印材料(PLA)参数Table 1. 3D printing material (PLA) parameters参数 线径/mm 密度/(g·cm−3) 熔点/(℃) 拉伸强度/MPa 断裂伸长率/% 弹性模量/MPa 取值 1.75 1.26 176 48 7 3500 表 2 缩尺试验相似数值Table 2. Similarity parameters of scaled test物理量 相似关系 相似常数 长度/位移/m CL 8.7 道砟密度/(kg·m−3) Cρ 1.0 内摩擦角 Cφ 1.0 轨枕弹性模量/MPa CE 10.3 轨枕集中荷载/N CF=CECL2 779.6 1.2 道床竖向静载试验和横向阻力试验方案
道床竖向静载试验采用逐级加载的方式,轨枕的竖向加载示意图和位移计现场布置如图4所示,位移计安装在轨枕两侧,分辨率为0.01 mm. 试验具体加载步骤为逐级放置砝码于钢轨上,对轨枕排施加竖向压力,记录每级荷载下轨枕的沉降情况. 每次试验采用相同夯实方法和道床尺寸,重复3次试验以减少人为误差.
道床横向阻力试验采用侧向逐级加载的方式,采用定滑轮对埋置于道床中的4根Ⅲ型轨枕排和双根X轨枕排施加横向推力. 轨枕的横向加载示意和位移计现场布置如图5所示. 轨枕一端设置位移计记录轨枕排的横向位移,得到道床横向阻力-横向位移曲线. 并用蓝色喷漆在道床肩砟喷涂3条位移参考线,观察试验前后道床的侧向位移情况.
2. 缩尺道床模型试验
2.1 轨枕竖向沉降对比
图6对比了Ⅲ型轨枕和X形轨枕的竖向荷载-沉降以及3次竖向荷载-沉降的平均变化趋势. 由图可见:在加载至约640 N,Ⅲ型轨枕和X形轨枕沉降随着荷载的增加线性增长,且Ⅲ型轨枕沉降增长速率明显高于X形轨枕;当荷载超过该值时,X形轨枕的荷载-沉降曲线出现明显的拐点,沉降速率增大,但仍小于Ⅲ型轨枕;相较Ⅲ型轨枕,X形轨枕最终沉降量分别减小了约0.21、0.18 mm和0.21 mm,平均沉降量减少了约26.3%. 因此,X形轨枕具有降低道床沉降的潜能.
2.2 道床的竖向刚度
静态轨道道砟层的刚度为加载区间内轨枕所受竖向荷载与轨枕沉降增量的比值[22]. 图7对比了3次重复试验的Ⅲ型轨枕和X形轨枕的道床刚度. Ⅲ型轨枕道床刚度分别为1.25、1.17、1.3 MN/m;X形轨枕道床刚度分别为1.81、1.75、1.89 MN/m;X形轨枕的平均道床刚度比Ⅲ型轨枕的平均道床刚度提高了约46.6 %,表明X形轨枕有更好的竖向承载能力.
2.3 道床横向位移
图8为Ⅲ型轨枕和X形轨枕横向荷载-位移曲线. 由图可知,轨枕的横推可以分为3个阶段:在阶段 ① 期间,位移随着荷载的增加缓慢增加,Ⅲ型轨枕和X形轨枕的横向荷载-位移曲线基本上呈线性关系,属于线性变形阶段,相应的Hcr为横向临塑荷载;在阶段 ② 期间,轨枕横向位移增长速率增大,且Ⅲ型轨枕的位移增加量大于X形轨枕,属于局部破坏阶段,相应的Hu为横向极限荷载;在阶段 ③ 期间,轨枕横向位移突增,远大于前期加载时的位移总量,表明肩砟已完全失稳破坏,此时属于整体破坏阶段,相应的Hf为横向破坏荷载. 当Ⅲ型轨枕处于极限状态时,X形轨枕仍处于线性变形阶段,其横向位移较Ⅲ型轨枕分别减小了约1.07、2.64、0.66 mm. 由此可见,X形轨枕横向承载力更大,道床结构更稳定.
此外,对比Ⅲ型轨枕和X形轨枕多次横推试验结果可知:Ⅲ型轨枕横向承载力特性有可能因道床夯实状态的变化而发生显著变化;相较于Ⅲ型轨枕,X形轨枕在3次试验中呈现结果的离散性更小,意味着采用X形轨枕时,轨道的鲁棒性更优,能够有效提升有砟轨道的稳定性.
图9为Ⅲ型轨枕和X形轨枕在不同横向加载阶段时道床的变形情况. 由图9可知:1) 当加载至Ⅲ型轨枕横向极限荷载Hu(Ⅲ)时,肉眼观察Ⅲ型轨枕和X形轨枕道床肩砟标记仍保持完整,X形轨枕道床肩砟部分与加载前无明显变化,可认为X形轨枕道床仍处于线性阶段;但试验中可观察到Ⅲ型轨枕道床肩部少量道砟滑落,Ⅲ型轨枕进入局部破坏阶段. 2) 当加载至Ⅲ型轨枕横向破坏荷载Hf(Ⅲ)时,Ⅲ型轨枕有明显的侧向移动,Ⅲ型轨枕道床肩砟完全被破坏,对应图8(a)中整体破坏阶段;X形轨枕肩砟标记仍保持良好,处于局部破坏阶段. 3) 当加载至Hf(X)时,X形轨枕道床肩砟大范围失稳,肩砟滑落和侧移,3根参考线几乎破坏,X形轨枕横向位移显著增加,对应于图8(b)中的整体破坏阶段. 因此,X形轨枕可有效提升道床横向阻力,提高极限承载力.
2.4 轨枕的横向极限阻力
图10对比了Ⅲ型轨枕和X形轨枕横向荷载-位移的平均变化趋势. 由图可知:在荷载至69.7 N时,Ⅲ型轨枕的横向位移梯度急剧增加,继续增加荷载,轨枕位移增量很大,即处于横向承载极限状态,其道床横向极限阻力约为69.7 N. 同理可得,X形轨枕道床的横向极限阻力约为85.2 N. 对比图10中的突变点可得出,X形轨枕的横向极限阻力比Ⅲ型轨枕提高了约22.4 %.
3. 缩尺与全尺寸道床离散元模拟
采用激光扫描技术生成薄片状、纺锤状和立方体状3种代表性形状道砟颗粒的三维图像. 通过设置球体单元的平滑程度与球体颗粒最小值,在三维模型内表面填充球体单元颗粒簇,生成接近真实道砟颗粒形状的道砟模型. 考虑到全尺寸道床模型计算量过大,采取简化形状道砟颗粒充填道床. 通过采用3种真实形状的道砟颗粒模型进行休止角试验,试验和离散元模拟结果表明,3种道砟颗粒道床休止角与实际休止角接近,能够较好地模拟道砟颗粒间互锁作用,如图11所示.
使用三维STL (Stereolithography)文件在EDEM (engineering discrete element method)环境中创建钢轨、轨枕和道床的几何形状,轨枕尺寸参考陈成等[13-14]研究. 由于本文离散元模拟采用静态加载的方式对轨枕施加横、竖向荷载. 轨枕与钢轨之间相对位移可忽略不计,钢轨与轨枕的连接通过在外接的Recurdyn软件中设置fix命令,将钢轨固定在轨枕上,限制钢轨和轨枕之间的相互移动. 通过颗粒工厂在道床区域内生成道砟颗粒,直到达到所需的道砟层孔隙率0.4. 随后轨枕在重力作用下自由下沉直至稳定. 轨枕与道砟,道砟与道砟的细观接触参数取值参考李朋[23]的研究,如表3所示.
表 3 模型物理参数与接触参数Table 3. Physical and contact parameters of model模型部件 泊松比 剪切模量/
(×1010 Pa)密度/ (kg·m−3) 摩擦系数 道砟 0.30 2.0 2700 0.55 轨枕 0.23 1.5 2800 0.85 钢轨 0.23 1.5 7800 0.85 为了从细观上对比分析Ⅲ型和X形轨枕的竖向及横向承载特性,首先,建立缩尺Ⅲ型和X形轨枕有砟轨道模型(图12(a)),并施加与前文缩尺试验一致的加载方案,以验证模型的可靠性;然后,建立全尺寸型和X形轨枕有砟轨道模型(图12(b)),模拟真实道床刚度和横向阻力试验,以对比真实工况2种轨枕下道床刚度与横向阻力. 模型建立完成后,对轨枕分别进行横、竖向加载,记录轨枕横向位移、竖向位移、道砟与轨枕之间的接触力以及颗粒运动情况.
4. 模拟结果与讨论
4.1 缩尺试验与缩尺模拟结果对比验证
图13为缩尺试验与缩尺模拟的竖向承载力-位移曲线和横向荷载-位移曲线. 结果表明:缩尺模型竖向荷载-沉降关系与缩尺试验结果相似;缩尺模型横向荷载-位移曲线加载前期与缩尺试验变化一致,X形轨枕的极限横向阻力更大. 然而,在横推破坏阶段,模拟中的轨枕横向位移均小于试验. 这可能是由于模拟中颗粒形状相对于真实颗粒的形状多样性低,导致颗粒间咬合降低. 为了进一步研究2种轨枕的工作特性,后续结果将重点进行全尺寸Ⅲ型轨枕和X形轨枕有砟道床模拟,对比真实工况下2种轨枕下道床刚度与横向阻力.
4.2 道床刚度和细观分析
图14为全尺寸Ⅲ型轨枕和X形轨枕道床的竖向承载力-位移曲线. 由图可知:Ⅲ型和X形轨枕的承载力随着位移呈线性增长,属于线性变形阶段,且X形轨枕承载力的增长速率明显大于Ⅲ型轨枕. 参考道床刚度现场试验测试方法,采集7.5 kN 和35.0 kN加载时相对应的位移读数. 根据相应竖向荷载差值与其沉降差值的比值计算道床刚度,计算可得,4根Ⅲ型轨枕排的道床刚度约为141.5 MN/m,符合实际工程中道床刚度要求;双根X形轨枕排的道床刚度约为182.8 MN/m,X形轨枕的道床刚度比Ⅲ型轨枕的道床刚度提高了约29.2%.
图15对比了Ⅲ型轨枕排和X形轨枕排竖向压载至0.5 mm时道砟颗粒间接触力分布情况,不同力链颜色分别代表不同接触力的大小. 由图可见:Ⅲ型轨枕和X形轨枕与道砟接触应力均主要集中在枕底位置下,呈发散状向下扩散,且距离轨枕越远接触力链衰减越明显;相较于Ⅲ型轨枕,X形轨枕与道砟的接触应力明显增加,轨枕与道砟的接触力沿X形4个夹角范围扩散,分布范围更广,使得道砟充分参与分担荷载;相比于Ⅲ型轨枕,X形轨枕的“V”形结构更稳定,具有较高的承载力.
4.3 道床横向阻力和细观分析
图16为4根Ⅲ型轨枕和2根X形轨枕道床的横向阻力-位移曲线. 由图可知:横向阻力在1 mm位移内快速增长,之后增长速率减小,在2 mm时横向阻力趋于稳定,取轨枕位移为2 mm时的轨枕推力为道床横向阻力值. 鉴于X形轨枕在设计初期参照2根条形轨枕中心交叉结合而制,2根X形轨枕等效于4根条形轨枕. 本文道床横向阻力试验模拟的4根Ⅲ型轨枕排的总横向阻力约为35 kN,相应的2根X形轨枕排的总横向阻力约50 kN. 可得出双根X形轨枕的横向阻力相较于4根条形轨枕提升了约31.6%,与缩尺道床横向阻力试验结果接近. 此外,通过Zakeri等[24]和Koike等[25]的研究发现增加枕木的数量并不能成倍地增加横向阻力,这是由于多轨枕排的横推试验中存在枕群效应,减小了单根轨枕作用效果.
轨枕道床横向阻力主要源于枕底道砟、枕间道砟以及肩部道砟的协同作用[26]. 通过提取2种轨枕与各区域道砟接触面的应力,计算得到道床横向阻力分担比,如图17所示. 由于单X形轨枕与双条形轨枕的枕底和端部与道砟接触面积均接近,两者在枕底摩阻力与端部阻力之和上也呈现出较为接近的数值,分别为27.6 kN和29.0 kN. 然而,轨枕间砟对X形轨枕侧边提供的阻力却显著超出条形轨枕,达到3倍. 结果表明X形轨枕主要是通过增强其与轨枕间砟接触,提升道床横向阻力.
图18对比了2种轨枕在横推2 mm时道床内部接力链的分布情况. X形结构轨枕由于其独特的设计,能够更好地嵌固于道床之中,从而显著改善了轨枕侧面与轨枕间砟的接触状态. 当条形轨枕横向移动时,由于间砟与轨枕相对位移较小,几乎处于静止状态,主要是间砟与轨枕侧面的接触在提供侧摩阻力. 然而,当X形轨枕横向移动时,小夹角内侧边间砟受到挤压并产生一定量位移,使得作用在X形轨枕小夹角侧边的压力属于被动压力. 这种被动压力要远大于条形轨枕侧边所承受的静止压力,因此能够更有效地抵抗横向移动. 图18(a)接触力链也很好证明这一结论,X形轨枕小夹角的侧边间砟和相邻X形轨枕的中心间砟受到挤压,导致间砟颗粒的应力以及间砟-轨枕接触力更集中.
选取轨枕A-A剖面,提取道床内部接触力链分布情况,如图18(b)所示. Ⅲ型轨枕的底部接触力分布较为均匀,但是,X形轨枕底部接触力分布主要集中于轨枕两端以及其对应的底部. 在X形轨枕的端部,力链的分布没有条形轨枕明显. 这是因为在小夹角内侧间砟的被动压力分担了X形轨枕端部的作用力,导致与Ⅲ型轨枕相比,X形轨枕的端部接触面积相同,但端阻力有所减小. 根据数值模拟结果可知,X形轨枕可显著提升轨枕的竖向承载力和横向阻力,尽管与缩尺试验结果有一定的差距,但结论趋势大致相同,后续将进行大尺寸模型试验进行验证.
5. 结 论
本文基于3D打印的缩尺试验和全尺寸离散元仿真模拟试验,研究了X形轨枕在竖、横向荷载作用下道床的变形情况和承载特性,并与传统Ⅲ型轨枕进行了对比. 主要结论如下:
1) 在竖向荷载试验中,X形轨枕的平均竖向沉降量比Ⅲ型轨枕降低了约26.3 %,同时X形轨枕的道床刚度比Ⅲ型轨枕提高了约46.6 %. 采用X形轨枕能够有效提升道床的竖向承载力.
2) 在横向阻力的试验中,X形轨枕的道床横向极限阻力提高了约22.4 %. 采用X形轨枕能够提升道床的横向承载力极限.
3) 道床竖向承载力模拟结果表明,双X形轨枕的道床刚度较四根Ⅲ型轨枕提升约29.2 %. 相较于Ⅲ型轨枕,X形轨枕与枕间道砟的接触明显增加,轨枕接触力沿X形4个夹角范围扩散,使得枕间道砟充分参与分担荷载.
4) 道床横向阻力模拟结果表明,X形轨枕的横向阻力较Ⅲ型轨枕提升约31.6 %. 其中轨枕间砟对X形轨枕侧边提供的阻力是条形轨枕的3倍. 因此,X形轨枕主要是通过增强其与轨枕间砟接触,提升道床横向阻力.
-
表 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. -