机车轴箱轴承冲击故障非平稳振动机理分析_吴冰冰.pdf

基金项目：湖南省自然科学基金资助项目（编号：2020JJ4316）收稿日期：20220620机车轴箱轴承冲击故障非平稳振动机理分析*吴冰冰，王靖（湖南科技大学 信息与电气工程学院，湖南湘潭411201）摘要：滚动轴承是机车车辆走行部的关键部件，其故障概率较大，振动机理复杂，具有非平稳性，故障诊断困难。为解决该问题，从研究机车轴箱滚动轴承故障发生的机理入手，对冲击故障振动信号加以研究。在实际工况下，着重考虑了载荷和工作转速这两个非线性非平稳性因素，并在此基础上建立了机车滚动轴承冲击故障的数学模型。通过所建立的数学模型，对其进行MATLAB仿真信号分析，同时以某大学轴承故障数据集中的外圈故障为例，开展进一步的MATLAB仿真分析，在一定程度上验证了机理分析和所建立数学模型的正确性。关键词：机车滚动轴承；故障振动机理；非平稳性；数学模型；仿真分析中图分类号：U264.7文献标志码：A文章编号：10099492(2023)03008904Analysis on Non-stationary Vibration Mechanism of Rolling Bearing ImpactFailure of LocomotiveWu Bingbing，Wang Jing（School of Information and Electrical Engineering,Hunan University of Science and Technology,Xiangtan,Hunan 411201,China）Abstract:Rolling bearings is the key components of the running part of locomotives Its fault probability is large,its vibration mechanism iscomplex,and it has non-stationarity,so the fault diagnosis is difficult.In order to solve this problem,the impact fault vibration signal wasstudied by studying the mechanism of locomotive axle box rolling bearing failure.In the actual working condition,the load and the workingspeed were considered,and on this basis,the mathematical model of locomotive rolling bearing impact fault was established.Through theestablished mathematical model,the MATLAB simulation signal analysis was carried out.At the same time,the outer ring fault in the bearingfault data set of a university was taken as an example to carry out further MATLAB simulation analysis,which verified the correctness of themechanism analysis and the established mathematical model to a certain extent.Key words:locomotive rolling bearing;failure vibration mechanism;nonstationarity;mathematical model;simulation analysis2023年03月第52卷第03期Mar.2023Vol.52No.03机电工程技术MECHANICAL&ELECTRICAL ENGINEERING TECHNOLOGYDOI:10.3969/j.issn.1009-9492.2023.03.017吴冰冰，王靖.机车轴箱轴承冲击故障非平稳振动机理分析 J.机电工程技术，2023，52（03）：89-92.0引言在机车运行过程中，轴箱轴承作为机车走行部的关键部件，很大程度上减少了轴与轴之间摩擦带来的损失，对零件之间的运转起到支撑作用，同时具备良好的承载能力。轴承长期工作在强振动、负载不均、温差大、多粉尘的恶劣工作环境中，是走行部中故障率较高的部件。轴箱轴承的主要用途是联系构架和轮对，不同于其他静止机械设备上的轴承元件，它不仅要承受机车静态压力，还要受到运行过程中由径向加速度引起的非稳定动态载荷，车速越高，振动越大，动态压力也越大1。轴箱轴承在运行过程中同时受摩擦阻力、负载变化、频繁加减速和装配间隙等非线性因素的影响，因此反映其运行状态的振动信号表现出明显的非线性非平稳性特征，这是机车滚动轴承故障诊断的难点，同样也是当下故障诊断研究的热点。尤其是在实际工况下运行时，轴承故障振动信号的非平稳性愈发明显，对其进行故障诊断十分困难，因此研究其故障发生机理是分析故障原因的重要手段。近年来，许多学者对轴承的故障机理展开了大量研究，本文参考相关文献对其故障机理进行研究2-5，针对载荷和工作转速两个非线性非平稳性因素，分析机车轴箱轴承冲击故障信号，确定故障信号类型；在实际工况下，着重考虑载荷和工作转速这2个非线性非平稳性因素，获得机车滚动轴承冲击故障数学模型，进行仿真验证，并将模型应用于实际工况下的轴承故障分析，验证所建立的故障机理模型。1机车轴箱轴承故障的振动机理分析1.1机车轴承的故障特征频率当机车轴承的某个元件发生故障时，故障点与滚动体相互碰撞所产生的频率是不同的，不同的元件有与之相对应的频率。所以，可以根据故障点与滚动体相互碰撞所产生的频率判断发生故障的元件，此频率称为“故障特征频率”。由轴承的转速和零件尺寸可计算具体元件 89的故障特频率6，如表1所示。表中：fn为轴承的转动频率，fn=N60，N为轴承的转速；Z为滚动体的个数；d为滚子直径；D为轴承节径；为接触角。1.2故障信号振动机理分析当轴箱轴承的某一部件出现损伤故障时，提取故障采样振动信号对其进行频谱分析，可以发现频谱图中通常存在以故障点所在部件的故障特征频率为中心频率，以轴承的转动频率围绕其两侧的调制频率的调制边频带；通过已知采样振动信号中提取到的调制信息，分析其频率成分和强度大小，可以初步判断轴箱轴承的故障位置以及故障程度7。1.2.1调幅分析机车轴箱轴承跟随安装轴转动时，参与滚动的滚动体个数交替变化，其振动信号中一般会存在周期性的瞬态冲击信号。当轴承无故障时，其振动信号的主要成分包含转动频率及其各次谐波频率，正常振动信号x(t)的表达式如下：x()t=m=0amcos()2fnmt+m（1）式中：am为转动频率第m阶谐波分量的幅值；fn为轴承的转动频率；m为转动频率第m阶谐波分量的相位。当轴箱轴承出现损伤类故障时，轴每转一圈，故障点就会和与之相接触的元件表面碰撞一次，故障振动信号的幅值也会产生幅值调制现象，冲击信号的强弱发生周期变化。因此，用am(t)表示幅值调制函数，计算公式为：am()t=n=0bncos()2fknt+n（2）式中：bn为幅值调制函数第n阶分量的幅值；n为幅值调制函数第n阶分量的相位；fk为故障点所在部件的故障特征频率。但是，机车在实际工况下运行时，受曲线、轨道不平顺等因素的影响，轴箱轴承载荷会发生变化，负载越大，振动信号的幅值越大，而轴箱轴承同时承受径向和轴向的联合载荷作用，因此需将实际工作载荷转换为当量负荷8，公式如下：Pr=ha Gr（3）Gr=G2v+G2H（4）式中：Pr为当量负荷；ha为轴向负荷系数；Gr为轴向轴承的径向载荷；Gv为垂直负荷；GH为水平负荷。载荷的变化同时也影响着轴承的接触应力和变形。分析滚动轴承的接触应力及变形的难度在于滚动体与圈体的接触。滚动体与内圈是凸面与凸面相接触，滚动体与外圈是凸面与凹面相接触，所以滚动体与内圈的应力集中程度大于滚动体与外圈的应力集中程度，同时内圈接触变形也会大于外圈接触变形9。在载荷为 0的情况下，滚动体与圈体的接触为点接触，随着载荷的不断增大，点接触变为面接触。当由点接触向面接触过渡时，由于变形较快，接触变形的增量较小；当接触变为面接触时，滚动体与轴承的接触面积基本上不会再发生变化，此时载荷的持续增加会导致接触变形的增量持续增大10。接触区域的位置、形状、大小、接触面压力以及摩擦力分布等接触参数随载荷的变化而变化，是典型的边界非线性问题，尤其是在实际的高速工况下，滚动体接触状态复杂，故障机理分析也复杂。除去载荷的影响，轴箱轴承的工作转速也对振动信号的幅值有一定影响。随转速不断提高，加速度变化极为不稳定，产生越来越大的离心力，导致滚动体与轴承内外圈之间的接触力显著增加，振动信号的幅值也随之增大。记工作转速为n，表达式为：n=h1h2nlim（5）式中：h1为载荷系数；h2为载荷分布系数；nlim为极限转速。因此，综合考虑载荷和工作转速这两个影响因素，幅值bn可表示为bn()t=FPr()t,n()t，则：bn()t=k=0cos()Prcos2kt ncos()2kt（6）故经过幅值调制后的信号模型表示为：ym()t=x()t 1+am()t（7）对调幅后的信号模型ym(t)进行傅里叶变换，得到频谱图并进行分析，可以看出信号频谱是以故障特征频率及其倍频为中心，以调制频率及其倍频为间隔而形成的多对调制边频带。机车轴承振动信号在实际运行过程中产生的载波信号和调制信号均属于多频率成分信号，因此通过傅里叶频谱变换会形成许多围绕故障特征频率及其高次谐波两侧的边频带11。1.2.2调频分析轴箱轴承在振动过程中产生的振动信号受转速波动、故障频率/Hz内圈故障特征频率finner外圈故障特征频率fouter保持架故障特征频率fc滚动体故障特征频率froller计算公式finner=Z2()1+dDcos fnfouter=Z2()1-dDcos fnfc=12()1-dDcos fnfroller=D2d|1-()dD2cos2fn表1机车轴承的故障特征频率公式2023年03月机 电 工 程 技 术第52卷第03期 90负载变化等非线性非平稳性因素影响，使得故障轴承所在轴的转动频率产生一定的变化，进而与故障轴承的故障特征频率发生频率调制现象12。因此，用fm(t)表示频率调制函数。fm()t=fn+Kfm0tx()d（8）式中：Kfm为调频波瞬时频率相对于fn的偏移，为轴频信号随时间变化的函数。对于机车轴承，载荷变化的同时不仅影响振动信号的幅值，也会使轴承的转速发生波动，这种转速波动会造成振动信号的调频现象，从而使振动信号的频率发生偏移。因为任何复杂的运动都可以分解为简单运动的组会，所以将轴承转速波动按简谐规律变化处理，记转速波动为n，表示为：n=n01+Asin()Nn0t（9）式中：n0为平均转速；A为幅值比；N为每周期的波动次数。所以可将频偏函数Kfm展开为：Kfm()t=l=0n01+Asin()Nn0t cos()2lt（10）因此实际轴承调幅调频振动信号的表达式为：y(t)=Fx(t),am(t),fm(t)，将变化的轴频信号公式（8）代入公式（7），替换轴频信号x(t)中的轴频fn，则：y()t=m=0amcos2fm()t mt+m1+am()t（11）为简单明了地分析表达式，对于上述调幅调频振动信号，全部取其一阶信号进行分析，则：y()t=acos2()fn+Kfm()t0tx()d t+1+b()t cos()2fkt+（12）在实际运行过程中，传感器测得的振动信号成分十分复杂，不仅包含故障对振动系统产生的冲击信号，还交杂着环境