Modal simulation analysis and experimental study of aerostatic spindle
-
摘要: 为准确获取气体静压主轴的动力学模态特性,分别建立主轴转子及主轴部件的有限元模型,对有限元模型中添加的弹簧及电机约束边界条件进行了研究,并分别对以上两个有限元模型进行模态计算。比较两种情况下的结构固有频率、振型趋势及各阶模态下的有效参与质量,分析主轴部件中的外壳对主轴转子模态参数的影响,其中仅保留主轴转子的前六阶固有频率为1.5 Hz 、277 Hz、279 Hz 、316 Hz、384 Hz、385 Hz,而主轴部件的前六阶固有频率为1.5 Hz 、259 Hz、270 Hz 、281 Hz、352 Hz、363 Hz,以上结果中均包含刚体模态及重根模态。两者比较结果表明主轴部件的固有频率因其质量的增加而有所降低,且主轴部件不再满足对称结构形式,因而相对于仅保留主轴转子重根模态的情况有所减弱。基于有限元仿真分析结果,进一步对主轴部件进行模态试验验证,通过仿真与试验对比,结果表明主轴转子的有限元仿真结果与试验测量结果更加贴近,其前三阶误差分别为0.6%、10.6%和2.6%。上述仿真边界条件的加载与计算过程及试验测试方法对进一步提高气体静压主轴模态参数仿真精度及修正动力学仿真模型提供了有益参考。Abstract: In order to accurately obtain the dynamic modal characteristics of aerostatic spindle, the finite element models of the spindle rotor and the spindle components are established respectively, the boundary conditions of spring and motor constraints added in the finite element model are studied, the modal calculation of the above two finite element models is carried out respectively. This paper makes a comparative analysis of the structural natural frequency, mode of vibration and effective participation mass of each mode under different conditions, and investigates the influence of the housing in the spindle components on the modal parameters of the spindle rotor. The first six natural frequencies of spindle rotor are 1.5 Hz 、277 Hz、279 Hz 、316 Hz、384 Hz、385 Hz. Correspondingly, the first six natural frequencies of the spindle components are 1.5 Hz、259 Hz、270 Hz、281 Hz、352 Hz、363 Hz. The above results include rigid body mode and doublet mode. The comparison results show that the natural frequency of the spindle component decreases due to the increase of its mass. Moreover, the spindle components do not meet the requirements of the symmetrical structure. Therefore, the result of the doublet mode of the spindle components is weakened. To verify the results of finite element simulation analysis, the modal test for spindle components is conducted. Through the comparison between simulation and test, the results show that the simulation results of the spindle rotor are closer to the experimental measurement results, and the first three order deviations are 0.6%, 10.6% and 2.6% respectively. The loading and calculation process of the above simulation boundary conditions and testing methods provide significant reference for further improving the calculation accuracy of modal parameters of aerostatic spindle and revising the dynamic simulation model.
-
Key words:
- aerostatic spindle /
- modal analysis /
- doublet mode /
- modal test
-
表 1 不同弹簧刚度及受力下的总体计算刚度
弹簧总数 单个刚度
/(N/μm)力方向 变形
/μm计算刚度
/(N/μm)8 180 轴向 0.19 526 径向 0.38 263 90 轴向 0.38 263 径向 0.77 130 表 2 主轴转子仿真结果模态参数
阶数 固有频率
/Hz参与系数 有效质量 X Y Z RX RY RZ X Y Z RX RY RZ 1 1.5 0 0 0 0 1 0 0 0 0 0 0.99 0 2 277 0.02 0 1 0.3 0 0 0.003 0 0.91 0.07 0 0 3 279 1 0 0.02 0 0 0.31 0.91 0 0.003 0 0 0.07 4 316 0 1 0 0 0 0 0 0.99 0 0 0 0 5 384 0.3 0 0.12 0.44 0 1 0.07 0 0.02 0.15 0 0.76 6 385 0.12 0 0.3 1 0 0.44 0.02 0 0.07 0.76 0 0.15 表 3 主轴部件仿真结果模态参数
阶数 固有频率
/Hz参与系数 有效质量 X Y Z RX RY RZ X Y Z RX RY RZ 1 1.5 0 0 0 0 1 0 0 0 0 0 0.05 0 2 259 0 0.9 1 0.24 0 0 0 0.25 0.15 0.02 0 0 3 270 1 0 0 0 0.05 0.2 0.29 0 0 0 0 0.01 4 281 0 1 0.9 0.48 0 0 0 0.29 0.12 0.07 0 0 5 352 0 0.16 0.35 1 0 0.01 0 0 0.02 0.32 0 0 6 363 0.3 0 0 0.01 0.01 1 0.03 0 0 0 0 0.4 表 4 传感器及激励点的位置与方向分布表
位置 传感器 方向 激励 方向 主轴部件 1 +Y 1 −Y 2 −Z 2 +Z 3 +X 3 +X 轴瓦 4 +Y 4 −Y 5 -Z 5 +Z 6 +X 6 +X 外壳左端
及中间7 +Y 7 −Y 8 −Z 8 +Z 9 +X 9 +X 10 +X 10 −X 11 +X 电机外套 11 +X 12 +Z 13 +X 外壳右端 12 −Y 14 +X 13 −Z 14 +X 15 +Y 表 5 仿真与试验结果对比
仿真结果/Hz 试验结果/Hz 与主轴转子的平均误差 与主轴部件
的平均误差阶次 主轴转子 主轴部件 激励点11 激励点13 激励点14 平均值 1 277~279 259~270 272.9 278.2 277.9 276.3 0.6% 4.5% 2 316 281 349.8 349.4 349.7 349.6 10.6% 24.4% 3 384~385 352~363 374.5 374.9 374 374.4 2.6% 4.7% -
[1] 胡秋, 夏仰球, 米良, 等. 超精密静压支承部件及其在超精密加工机床的应用[J]. 机床与液压, 2015, 43(22): 4-6,24. doi: 10.3969/j.issn.1001-3881.2015.22.002 [2] 潘高星. 液体静压功能部件发展及应用[J]. 机床与液压, 2013, 41(20): 137-139. [3] 管峰, 李圣怡, 关朝亮, 等. 气体静压轴承设计及关键零件加工工艺研究[J]. 机械科学与技术, 2016, 35(7): 1083-1088. doi: 10.13433/j.cnki.1003-8728.2016.0716 [4] 侯国安. 流体静压支承对超精密金刚石车床动态特性影响的研究[D]. 哈尔滨: 哈尔滨工业大学, 2013. [5] 熊万里, 原帅, 胡灿, 等. 液体静压主轴的回转精度规律及其极限预测[J]. 机械工程学报, 2021, 57(13): 70-82. [6] 安晨辉, 许乔, 张飞虎. 空气静压主轴气膜的压强分布与动态特性分析[J]. 纳米技术与精密工程, 2009, 7(5): 459-468. [7] 蔡力钢, 马仕明, 赵永胜, 等. 多约束状态下重载机械式主轴有限元建模及模态分析[J]. 机械工程学报, 2012, 48(3): 165-173. [8] 李茂源. 高速永磁电机机械强度与转子模态分析[D]. 青岛: 山东科技大学, 2019. [9] 郑晋, 韩明, 黄文祥, 等. 车用永磁同步电机转子受力浅析[J]. 上海汽车, 2015(5): 18-21. [10] 赖文海, 黄开胜, 周游. 永磁同步电机转子偏心电磁力和挠度的分析与计算[J]. 电机与控制应用, 2022, 49(1): 74-79,95. doi: 10.12177/emca.2021.160 [11] 武新峰, 雷勇军, 李家文. 基于有限元分析的复杂结构弹性振动传递函数建模[J]. 振动与冲击, 2012, 31(11): 57-61,77. doi: 10.3969/j.issn.1000-3835.2012.11.013 [12] 曹志浩. 矩阵特征值问题[M]. 上海: 上海科学技术出版社, 1980 [13] 徐洪, 王扬渝, 文东辉. 基于多点脉冲激励的主轴系统运行模态分析方法[J]. 机械工程学报, 2014, 50(19): 175-181. [14] Avitabile P.Modal testing: a practitioner's guide[M]. NY:John Wiley&Sons, 2017. [15] 谭详军. 从这里学NVH [M]. 北京: 机械工业出版社, 2021. [16] 刘进明, 应怀樵, 谭祥军, 等. 振动和声的互易性原理在结构模态试验中的应用[C]. 现代振动与噪声技术(第九卷), 2011: 156-165. [17] 谭祥军, 张伟. 某汽车制动盘MIMO模态测试分析[C]. 第二十三届全国振动与噪声控制学术会议论文集, 2010: 486-494. [18] 王荣辉, 王战铎. 计算机控制系统中采样周期的选择方法[J]. 陕西工学院学报, 1997(2): 41-46. [19] 陈茉莉, 李舜酩. 信号源识别的相干函数法[J]. 中国机械工程, 2007(1): 95-100. doi: 10.3321/j.issn:1004-132X.2007.01.024