机床主轴机构模态和振动特性分析与主轴径向回转误测试方法的研究.pdf

Modeling and Simulation 建模与仿真建模与仿真,2023,12(4),3968-3977 Published Online July 2023 in Hans.https:/www.hanspub.org/journal/mos https:/doi.org/10.12677/mos.2023.124362 文章引用文章引用:李希铭.机床主轴机构模态和振动特性分析与主轴径向回转误测试方法的研究J.建模与仿真,2023,12(4):3968-3977.DOI:10.12677/mos.2023.124362 机床主轴机构模态和振动特性分析与主轴机床主轴机构模态和振动特性分析与主轴 径向回转误测试方法的研究径向回转误测试方法的研究 李希铭李希铭 上海理工大学机械工程学院，上海 收稿日期：2023年6月7日；录用日期：2023年7月17日；发布日期：2023年7月24日 摘摘 要要 本文以一特定机床主轴作为研究对象，首先利用传递矩阵法对主轴的参数进行确定、简化轴承等模型，本文以一特定机床主轴作为研究对象，首先利用传递矩阵法对主轴的参数进行确定、简化轴承等模型，使用有限元法对主轴模型进行带预应力的模态仿真。随后，本文使用垂直两点法测试了多组不同转速下使用有限元法对主轴模型进行带预应力的模态仿真。随后，本文使用垂直两点法测试了多组不同转速下的主轴回转误，并使用频域两点法对误差进行分离计算得到同步误差和异步误差。结果显示，由异步误的主轴回转误，并使用频域两点法对误差进行分离计算得到同步误差和异步误差。结果显示，由异步误差异常高导致的回转误差异常高转速区间的转速基频与模态分析结果的第一、四、五阶模态频率相对应。差异常高导致的回转误差异常高转速区间的转速基频与模态分析结果的第一、四、五阶模态频率相对应。为提高机床主轴的精度和稳定性提供理论依据和实验基础。为提高机床主轴的精度和稳定性提供理论依据和实验基础。关键词关键词 振动模态分析，传递矩阵法，主轴回转误差，垂直两点法振动模态分析，传递矩阵法，主轴回转误差，垂直两点法 Study on Analysis of Modal and Vibration Characteristics and Test Method of Radial Rotation Error of Machine Tool Spindle Ximing Li School of Mechanical Engineering,University of Shanghai for Science and Technology,Shanghai Received:Jun.7th,2023;accepted:Jul.17th,2023;published:Jul.24th,2023 Abstract In this article,a specific machine tool spindle is used as the research object.First,the transfer ma-trix method is used to determine the parameters of the spindle,simplify the bearing and other models,and the finite element method is used to perform modal simulation of the spindle model 李希铭 DOI:10.12677/mos.2023.124362 3969 建模与仿真 with prestress.Subsequently,this paper used the vertical two-point method to test the rotation er-rors of multiple sets of spindles at different speeds,and used the frequency domain two-point me-thod to separate the errors and calculate the synchronization error and asynchronous error.The results show that the fundamental frequency of the rotation speed interval of the abnormally high rotation error caused by the abnormally high asynchronous error corresponds to the first,fourth,and fifth-order modal frequencies of the modal analysis results.It provides a theoretical basis and an experimental basis for improving the accuracy and stability of the machine tool spindle.Keywords Spindle Vibration Modal Analysis,Transfer Matrix Method,Spindle Rotation Error,Vertical Two-Point Method Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 随着数控机床的广泛应用，机床主轴的精度和稳定性对加工质量和效率的影响越来越大。其中，主轴径向回转误差是影响主轴精度和稳定性的重要因素之一1。对回转误差的测试和振动模态的研究可用于机床加工的预测和补偿，提高零件的加工精度和表面质量，对于提高机床主轴的精度和稳定性具有重要意义2。目前，国内外学者已经对机床主轴的模态仿真计算方法和回转误差测试方法进行了广泛的研究。在主轴模态计算仿真方面，朱冬等人3通过模态分析得到了电主轴的固有频率、振型和临界转速对数控铣床电主轴故障诊断为高速电主轴结构设计和轴承最佳预紧力的确定提供了分析依据。徐壮4详尽的给出了联轴器传动型直结式铣削主轴的结构分析，并为装配与调试提供了改进方法。刘成颖等人5基于 Hertz理论和最小二乘法理论建立了高速电主轴中主轴机壳振动传递力学模型，利用 Ansys 软件求解出电主轴固定基座的刚度，并进行实验验证与仿真分析。黄康等人6采用多体传递矩阵法对电主轴的某些参数进行简化与确定，在此基础上采用有限元仿真方法建立电主轴仿真模型，并使用仿真软件 Hypermesh 进行 Optistruct 器件的振动仿真，以确定系统的振动特性。王玉伟等人7分别建立主轴转子及主轴部件的有限元模型，对有限元模型中添加的弹簧及电机约束边界条件进行了研究，并分别对以上两个有限元模型进行模态计算。李树森等人8设计了一种精密气悬浮主轴并在 ANSYS Workbench 中建立主轴转子的有限元模型，对设计转速进行模态分析和验证。对于主轴回转误差测试，国际标准 ISO 230-7 9制定了机床旋转轴的几何精度测试标准。陆峰等人10建立误差分离数学模型，对垂直两点法测量电主轴回转精度进行分析，利用垂直两点法分离误差方法对实验数据进行处理，得到电主轴在不同转速下的回转精度，找出影响电主轴回转精度的因素并进行了归纳。美国雄狮精密公司(Lion Precision)11研发了主轴回转误差测量仪，使回转误差的测试可以做到实时检测与快速分析。日本的青木保雄等人12提出了 3 点法圆度测定法，使用三个位移传感器动态监测回转轴的圆度，分离回转误差。韩玉稳等人13对机床主轴精度进行了检测，并用最小二乘法拟合测试结果和机床传感器数据，给出了精度补偿的数学模型。迟玉伦等人14提出了基于多步误差分离法的主轴回转误差检测方法，在不借助标准球或者标准棒的情况下对数控铣床的主轴外轮廓进行直接测量，并通过数据降噪处理后分离得到主轴回转误差。但是，以上研究局限于对振动模态的分析和对回转误差分离方法的研究，而影响主轴加工性能的因Open AccessOpen Access李希铭 DOI:10.12677/mos.2023.124362 3970 建模与仿真 素是多方面的，单一的研究对象无法与主轴系统性的联系起来。本文将从机床主轴的结构和振动特性入手，探讨机床主轴机构模态和振动特性对主轴径向回转误差的影响。为提高机床主轴的精度和稳定性提供理论依据和实验基础。2.直结式机械主轴回转误差分离方法和数学模型的建立直结式机械主轴回转误差分离方法和数学模型的建立 2.1.直结式主轴结构特性及组成直结式主轴结构特性及组成 本文研究对象为沈阳机床 i5m1 立式数控钻铣床的直结式主轴，其结构示意图如图 1 所示。主轴结构包括芯轴、拉杆、外套、轴承、套圈、隔圈、拉爪、隔环以及各种销、垫圈、碟簧等。该主轴使用 BT30 标准刀具，采用气动拉刀装置。在结构上，该主轴的回转轴芯轴由联轴器与主轴电机等驱动机构直接连接。Figure 1.Schematic diagram of spindle structure 图图 1.主轴结构示意图 从结构上看，该直结式主轴主要特点有：1)可实现主轴转速的高速化，有效改善主轴高速综合性能；2)机械结构简单，转动惯量小，能实现高速度和高的角加速度，有利于快速准停；3)减少齿轮传动或皮带传动等的中间传动环节，避免传动链导入的振动和噪声，提高了主轴回转精度。2.2.基于传递矩阵法的主轴动态特性分析基于传递矩阵法的主轴动态特性分析 针对本文这种结构方案的主轴系统，使用传递矩阵法来进行分析建模，设置参数进行振动特性分析来反映电主轴系统的动态特性。传递矩阵法即将系统化分为若干单元，在相邻两个小单元的分界面上，用位移协调和力的平衡条件来相互联系，且每个小单元可运用到牛顿第二定律建立运动方程6。根据运动方程及边界条件得到单元另一侧的力和位移。依此进行下去最后可求得问题的解。通常把分界面上的力和位移变量的列向量记作状态向量。解题从边界开始逐步另一方向，且相邻单元的状态向量用矩阵联系，将这一方法一般称为传递矩阵法。依据传递矩阵法的相关建模原理，针对所研究的电主轴振动系统建立的动力学集总参数模型如图2，在实际计算中，为简化问题，可将圆盘及右边的轴段组合为一构件，如图 3(a)所示，该构件的传递矩阵推导如下。具有弹性支承的圆盘，受力分析见图 3(b)。李希铭 DOI:10.12677/mos.2023.124362 3971 建模与仿真 Figure 2.Simplified model of spindle transfer matrix method 图图 2.主轴传递矩阵法简化模型 Figure 3.Force analysis 图图 3.受力分析 当转子进动角速度=(转子自转角速度)时，系统作同步正向进动，得到圆盘的惯性力为2iimY，惯性力矩为：()2dipiiJJ 由达朗伯原理得到：()22RLiiiriiRLiidipiiiiQQmk YMMJJk=+=+(1)且有：RLiiiRLiiiQQQYYY=(2)综合以上得到具有弹性支承的圆盘的传递矩阵为：()2210000100010001ipidiiiriDJJkmk=+(3)dJ转子直径转动惯量；pJ转子极转动惯量；im节点处总质量(轴段与圆盘质量之和)。同理可以推导无弹性圆盘的传递矩阵iC。李希铭 DOI:10.12677/mos.2023.124362 3972 建模与仿真 ()2210000100010001ipidiiCJJm=(4)故由以上的传递矩阵得到组合构件如图 3(a)所示的传递矩阵U。()()()()()()2322232232222111000260100010102200100100011(1)1622612iiiiiiiiipidiiiiiriiiiiiiriipdiiiiiiiriilllEIEIllUB DJJkEIEImklllllmklJJkEIEIEIEIllmkEI=+=()()()()2222221001iiipdiiiiiiripdiiirillJJkEIEIEIlmkJJklmk+(5)以上得到的传递矩阵,iiiDCBU作为相应的iU，于是，电主轴系统总的传递矩阵：1321nnUU UU U U=(6)其中，总的传递方程：1nZUZ=(7)代入初始边界条件求解，得到总的传递方程，并对应求解出相应的状态向量，即获取系统振动的振型与固有频率。2.3.主轴回转误差分离的垂直两点法主轴回转误差分离的垂直两点法 主轴回转误差是主轴实际回转轴线对其理想回转轴线的位置偏差。可以分解为三种情况：径向误差、轴向误差、倾角摆动三种基本形式。主轴回转径向误差运动是指在回转过程中回转轴线偏离理想轴线位置而出现的附加运动，是评价机床动态性能的一项重要指标，是影响机床工作精度的主要因素。两点测量法是根据日本学者青木保雄和大园成夫在 1966 年提出的三点测量法的基础上，将三个传感器中的两个传感器合并，演化成两点测量法12。三点测量法的测量原理是使用三个同样的传感器按照一定角度垂直布置被测标准球的同一截面上对主轴误差进行测量，并利用频谱技术把圆度误差与回转误差分离。根据ISO 230-7所定义的主轴回转误差测试标准，使用两点法分离主轴回转误差。其原理如图4所示。对处理后的原始信号做傅里叶变换，将转速基频及其倍频部分分离为同步误差，其余部分为异步误差，同步误差还可以进一步分离，其的第一项为偏心误差，其余项为圆度误差。根据三点法误差分离基本方程()()()()()()()()012cossincossinSrxSrxySrxy=+=+=+(8)李希铭 DOI:10.12677/mos.2023.124362 3973 建模与仿真 Figure 4.Schematic diagram of spindle rotation error separation 图图 4.主轴回转误差分离示意图 引入传感器的标定系数012ccc、，则以上三式分别乘标定系数得到：()()()()001122Cc Sc Sc S=+(9)根据(8)式展开，若要分离出圆度误差()r，则需要()x、()y前的系数为 0，则得到：01212coscos0sinsin0ccccc+=+=(10)此线性方程组有无穷多解，取01c=，得到一组特解为：()()0120sinsinsinsinccc=(11)在三点法中，若取2=、=(采样间隔2N=)当 N 较大时，(11)式中的1c、2c可以写成：()()()()()12sinsinsin0sincossin2sinsin21sinsin2cc=+=+李希铭 DOI:10.12677/mos.2023.124362 3974 建模与仿真 当10c=时，传感器1S的输出信号在加权中占有很小的比重，可以忽略不记。那么，那么三点法就演化成了垂直两点法。(9)式就可以写成：()()()02Cc rc r+(12)将(12)式离散化得到：()()0212NC ic r ic r i+(13)对(13)式进行离散傅里叶变换(FFT)，得到：()()()ffRn W nCn (14)()fRn为时域上电主轴的圆度误差()r i，()fCn为时域上圆度误差和回转误差的混合定好；()W n为权函数，如下：()22221e122jnpNW ncNNp+=(15)当()0W n 时，有：()()()ffRnCnW n (16)对(16)式进行离散傅里叶逆变换(IFFT)，得到：()()1ifr nFRn (17)将(17)式带入(8)式可以得到：()()()()()()022xSir iySir ip=+(18)由此得到分离后的同步误差和异步误差。同步误差可以进一步分离为圆度误差和偏心误差。3.电主轴模态分析电主轴模态分析 根据传递矩阵法，选取目标主轴的芯轴进行分析，使用 ANSYS 软件，建立有预应力的模态分析。根据芯轴实际受力情况进行加载和约束，划分网格后计算求解。分析得到的频率和振型图结果如下，表 1 为前 18 阶模态的频率计算结果，取前 6 阶绘制模态振型图如下。Table 1.Modal analysis simulation results 表表 1.模态分析仿真结果 模态 1 2 3 4 5 6 7 8 9 频率(Hz)1219.5 1223.1 2369.6 3764.8 4588.1 4603.4 5059.1 7616.8 7638.2 模态 10 11 12 13 14 15 16 17 18 频率(Hz)7813.2 7891.9 7989.2 8073.4 8445.6 9072.7 9389.8 9832.9 11549.0 前 6 阶模态结果绘图如图 5 所示。李希铭 DOI:10.12677/mos.2023.124362 3975 建模与仿真 (a)一阶振型 (b)二阶振型 (c)三阶振型 (d)四阶振型 (e)五阶振型 (f)六阶振型 Figure 5.Simulation result mode diagram 图图 5.仿真结果振型图 Figure 6.Spindle rotation error test platform 图图 6.主轴回转误差测试平台 李希铭 DOI:10.12677/mos.2023.124362 3976 建模与仿真 可以看到，对于主轴轴端的径向回转误差影响较大的几个模态分别是一阶、四阶、五阶和六阶。4.主轴回转误差测试主轴回转误差测试 使用美国雄狮精密 SEA 主轴回转误差分析仪对该主轴进行测试。测试平台如图 6(a)所示，其中的图6(b)为用来安装传感器的传感器巢和电容式传感器以及安装在主轴上的标准球，电容式传感器插入传感器巢的孔中并由六角螺丝固定，标准球被固定在一定心刀具座(BT40)中，由主轴上的打刀缸和拉爪固定在主轴上。图(c)为传感器驱动模块和数据采集卡，驱动模块包括电源和校零标尺，采集卡使用 NI-6306 采集卡，通过 USB 连接计算机。图(d)为数据采集和分析软件界面，用来设置采样率、平均圈数等参数，并保存测试结果和原始数据。将电容式传感器安装在传感器巢中并校零，设置每 20 圈进行一次平均，间隔 100 rpm 测试主轴从 100 rpm 到 8000 rpm 转速区间的回转误差，计算分离后的同步误差和异步误差值，得到在不同转速下的回转误差结果，部分数据结果如表 2 所示。Table 2.Rotation error analysis results 表表 2.回转误差分析结果 转速(r/min)1000 2000 3000 4000 5000 6000 7000 8000 回转误差(m)0.81 0.76 0.84 0.85 0.92 0.91 1.10 1.42 同步误差(m)0.37 0.27 0.30 0.32 0.40 0.36 0.45 0.50 异步误差(m)0.62 0.55 0.64 0.62 0.71 0.67 0.77 1.11 将所有数据汇总，绘制折线图如图 7 所示。Figure 7.Line chart of error separation result 图图 7.误差分离结果折线图 总的来看，主轴回转误差随转速的提高不断提高，但增加并不是线性的，在 35004000 rpm、54005900 rpm、72008000 这几个转速区间出现了异步误差异常高的情况，导致了整体回转误差偏高。而同步误差随转速的增加不断增加，较为线性。李希铭 DOI:10.12677/mos.2023.124362 3977 建模与仿真 5.结论结论 传递矩阵法主要解决一维结构系统静力学和动力学问题，其建立模型灵活，计算效率高为该方法的优点，对于本文复杂的主轴系统，在采用有限元法的同时，利用传递矩阵法对主轴的参数进行确定、简化轴承等模型，上述工作在有限元分析时起到突出的作用。用垂直两点法进行测量，降低了测量装置的安装难度，装配容易，能很容易地测得主轴的混合误差，并能够分离出同步误差和异步误差，直观地显示回转误差的轨迹图，做出进一步分析。已知主轴电机的极对数 P=6，则三个转速异常高区间的转速基频与模态中的第一、四、五阶模态频率相对应。则可以确定主轴在该转速区间回转精度较差是由于主轴振动导致的。在编制加工工艺参数过程中，应尽量避开该转速区间，使主轴运行更加平稳，进而保证加工质量稳定。参考文献参考文献 1 顾伟德,周亮,万春鹏.主轴回转误差测量技术研究现状J.机械制造,2021,59(9):75-80+84.2 马军旭,赵万华,张根保.国产数控机床精度保持性分析及研究现状J.中国机械工程,2015,26(22):3108-3115.3 朱冬,王玲,殷国富,等.数控铣床电主轴故障诊断及仿真分析J.工具技术,2019,53(1):125-128.4 徐壮.联轴器传动型直结式铣削主轴的装配与调试J.橡塑技术与装备,2016,42(10):46-47.https:/doi.org/10.13520/ki.rpte.2016.10.017 5 刘成颖,郑烽,王立平.高速电主轴中主轴-机壳振动传递力学模型J.清华大学学报(自然科学版),2018,58(7):671-676.https:/doi.org/10.16511/ki.qhdxxb.2018.21.015 6 黄康,张卫霞.基于 HyperMesh 的电主轴系统动态性能建模与仿真J.系统仿真学报,2013,25(6):1393-1398.https:/doi.org/10.16182/ki.joss.2013.06.004 7 王玉伟,吴东旭,索奇,等.气体静压主轴模态仿真分析及试验研究J.制造技术与机床,2023(2):12-20.https:/doi.org/10.19287/j.mtmt.1005-2402.2023.02.002 8 李树森,任毅.一种精密气悬浮主轴的结构设计及其特性研究J.润滑与密封,2017,42(4):80-83+88.9 ISO 230-7(2015)Test Code for Machine Tools.Part 7:Geometric Accuracy of Axes of Rotation,International Stan-dardization Organization.Geneva,Switzerland.10 陆峰,沈鸿波,张啸飞,张丽秀,吴玉厚,张珂.基于垂直两点法的电主轴回转精度分析J.组合机床与自动化加工技术,2014(11):14-18.https:/doi.org/10.13462/ki.mmtamt.2014.11.004 11 Precision,L.(2022)Spindle Measurement RPM and Bandwidth.https:/ 12 青木保雄,大圆成夫.一种基于 3 点法的圆度测量法的展开J.精密机械,1966,32(12):831-836.13 韩玉稳,董应明,王林,等.机床主轴精度检测试验J.机电产品开发与创新,2023,36(2):155-157.14 迟玉伦,王国强.基于多步误差分离法的主轴回转误差检测方法J.轻工机械,2023,41(2):48-56+65.