CN111337255B - 一种滚动轴承剥落位置判定方法 - Google Patents

一种滚动轴承剥落位置判定方法 Download PDF

Info

Publication number
CN111337255B
CN111337255B CN202010212170.1A CN202010212170A CN111337255B CN 111337255 B CN111337255 B CN 111337255B CN 202010212170 A CN202010212170 A CN 202010212170A CN 111337255 B CN111337255 B CN 111337255B
Authority
CN
China
Prior art keywords
bearing
fault
outer ring
inner ring
rolling
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010212170.1A
Other languages
English (en)
Other versions
CN111337255A (zh
Inventor
石怀涛
刘子濛
白晓天
李献文
李思慧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenyang Jianzhu University
Original Assignee
Shenyang Jianzhu University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenyang Jianzhu University filed Critical Shenyang Jianzhu University
Priority to CN202010212170.1A priority Critical patent/CN111337255B/zh
Publication of CN111337255A publication Critical patent/CN111337255A/zh
Application granted granted Critical
Publication of CN111337255B publication Critical patent/CN111337255B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/04Bearings
    • G01M13/045Acoustic or vibration analysis

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明提供了一种滚动轴承剥落位置判定方法,首先建立剥落故障轴承的非线性动力学模型,得到剥落故障轴承的模拟频域振动曲线,然后采集故障或疑似故障轴承上的振动信号,提取实测频域振动曲线,其次计算出故障或疑似故障轴承外圈、内圈的故障特征频率,以模拟频域振动曲线作为参照与实测频域振动曲线进行对比,确定待测轴承是否为故障轴承,并确定出具体故障位置。本发明避免了故障非线性动力学模型对实际故障模拟不准确的问题,同时考虑滚动体运行经过内外圈故障位置过程中的脉冲调制,可以更加精确地模拟出轴承发生早期故障时的运行状态,本方法不仅有效提高检测效率、精度,而且还有利于降低设备维护成本,提高生产效率。

Description

一种滚动轴承剥落位置判定方法
技术领域
本发明涉及轴承故障诊断技术领域,具体地说是一种滚动轴承剥落位置判定方法。
背景技术
剥落故障是造成滚动轴承失效的主要原因之一,通过对滚动轴承剥落故障进行统计分析,并与轴承振动响应的频域信号进行比对,对滚动轴承是否出现剥落故障做出判断,并区分剥落故障的位置处于轴承的内圈或者外圈,以便及时维修或更换,以防因剥落故障过大造成轴承失效,从而引起整机发生故障,造成不必要的经济损失和有可能的人员伤亡,使轴承工作过程始终处于受控状态。现有的轴承剥落故障判别方法无法准确判断出剥落故障发生于轴承的内圈或者外圈,对轴承系统的维护和更换造成了很大的影响。
现有的对轴承故障位置的判别方法大多为通过模拟轴承在运行过程中滚动体与内外圈故障接触时产生的接触变形量,通过Hertz接触定律来模拟轴承内外圈与滚动体的接触刚度,进而建立轴承故障非线性动力学模型,对故障轴承运行时的振动响应进行模拟,将模拟的仿真振动曲线与实际测得的故障轴承运行时的振动信号进行比对,以此判断轴承故障的发生位置。
现有的判断轴承故障发生位置的非线性动力学模型无法准确计算轴承发生故障时滚动体与内外圈的接触刚度,同时也无法准确的模拟出轴承在发生早期故障时的振动信号,因此,现有的判断轴承故障发生位置的方法只能对故障发展到十分严重的程度时的轴承进行诊断,很大程度上会影响设备的正常运行、维护与维修,更有甚者会造成极大的经济损失。
开发并实现可判断滚动轴承剥落故障是否发生及发生位置的监控方法,在故障发生的初期对轴承故障进行识别和定位,配合企业内部设备运转信息,可以实现企业生产过程的智能化和信息化,它的设计及应用对把握轴承运转寿命,提高设备运转效率,降低生产成本,提高加工与生产安全性等各个方面起到积极的辅助作用。
发明内容
本发明首先模拟剥落故障轴承滚动体与内外圈的修正刚度,同时考虑滚动体经过内外圈故障位置处时产生的脉冲调制,然后建立剥落故障轴承的非线性动力学模型,这样避免了以往需要模拟滚动体与轴承内外圈之间的接触变形量来计算接触刚度所造成的不准确性,再者滚动体经过内外圈故障位置处时产生的脉冲调制可以对轴承发生早期故障进行反映,因此,无需等到轴承故障发展到严重影响设备正常运行的程度就可以对故障进行识别。
在将非线性动力学模型模拟出的故障轴承运转仿真曲线与实测运动信号进行比对时,通过计算得出的轴承故障特征频率,将实测信号曲线的故障特征频率所对应的峰值与模拟信号的故障特征频率所对应的峰值进行比对,判断轴承故障的发生位置。
为实现上述目的,本发明提供的技术方案包括如下步骤:
步骤1:建立剥落故障轴承的非线性动力学模型,得到剥落故障轴承的模拟频域振动曲线,包括如下步骤:
步骤1.1:首先采用集中参数模型模拟剥落故障轴承的运行状态,然后根据轴承滚动体与轴承内/外圈的剥落故障位置得到轴承滚动体与轴承内/外圈的接触刚度变化过程,分别计算不同变化过程下的轴承内/外圈修正接触刚度值;
步骤1.2:利用公式(1)计算剥落故障轴承的总接触刚度值K;
Figure BDA0002423206310000021
式中,Ko1为滚动体与轴承外圈接触时的外圈修正接触刚度值,Ki1为滚动体与轴承内圈接触时的内圈修正接触刚度值;
步骤1.3:利用公式(2)计算剥落故障轴承高速旋转时的调制脉冲值ffault
Figure BDA0002423206310000022
式中,d(t)为轴承内圈在轴承载荷方向上的最大位移,Qmax为轴承所受的最大载荷,δmax为最大接触变形量,c为相对间隙,f为轴承转动频率,t为轴承运行时间;
步骤1.4:建立剥落故障轴承的非线性动力学模型如公式(3)所示,
Figure BDA0002423206310000023
式中,m为轴承质量,C为轴承系统阻尼,F为轴承系统所受外载荷;
步骤1.5:根据公式(3)建立剥落故障轴承的仿真模型,并进行仿真运算,得到剥落故障轴承模拟运行时的时域振动信号,根据模拟运行时的时域振动信号得到模拟时域振动曲线,对模拟运行时得到的时域振动信号进行傅里叶变换,得到剥落故障轴承模拟运行时的模拟频域振动曲线;
步骤2:提取实测轴承的实测频域振动曲线,在实测轴承的外圈任意点处设置接触式振动传感器,所述实测轴承包括实际需要检测的故障轴承或者疑似故障轴承,当实测轴承运行时,接触式振动传感器采集实测轴承外圈上的振动信号,采集到的振动信号通过数据采集装置处理后得到实测振动时域信号,对实测振动时域信号进行傅里叶变换,得到实测轴承的实测频域振动曲线;
步骤3:利用公式(4)、(5)分别计算出实测轴承外圈、内圈的故障特征频率,以模拟频域振动曲线作为参照与实测频域振动曲线进行对比,分别查看故障特征频率所对应的峰值,若两组曲线中存在峰值所对应的故障特征频率相差不大于Δf,则定义相差不大于Δf的故障特征频率为相似故障特征频率,Δf的选取方法如下:
Figure BDA0002423206310000031
检查所述相似故障特征频率是否为轴承内圈或者外圈的故障特征频率,若所述相似故障特征频率为内圈的故障特征频率或者外圈的故障特征频率,则判断实测轴承存在剥落故障,若所述相似故障特征频率既不是轴承内圈的故障特征频率也不是轴承外圈的故障特征频率,则判断实测轴承不存在剥落故障,如果所述相似故障特征频率存在于实测轴承的内圈,则故障位置为轴承内圈,如果所述相似故障特征频率存在于实测轴承的外圈,则故障位置为轴承的外圈;
Figure BDA0002423206310000032
Figure BDA0002423206310000033
式中,fbpfo为轴承外圈发生剥落故障时的故障特征频率,fbpfi为轴承内圈发生剥落故障时的故障特征频率,Z为滚动体个数,f为轴承转动频率,d为轴承内径,D为节圆直径,α为轴承接触角。
所述的步骤1.1具体表述为:
步骤1.1.1:根据轴承滚动体与轴承外圈的剥落故障位置将滚动体与轴承外圈的接触刚度变化分为三个外圈故障过程,分别为滚动体进入外圈故障区域的过程Io、滚动体处于外圈故障区域的过程IIo、滚动体离开外圈故障区域的过程IIIo,利用公式(6)计算滚动体分别在所述三个外圈故障过程中的外圈修正接触刚度值Ko1
Figure BDA0002423206310000041
式中,Ko为滚动体与外圈之间的接触刚度,Ro为外圈剥落修正系数,β为滚动体的角位置,
Figure BDA0002423206310000042
为滚动体进入外圈剥落边缘时的角位置,
Figure BDA0002423206310000043
为滚动体离开外圈剥落边缘时的角位置,
Figure BDA0002423206310000044
为外圈剥落中心的角位置,R为剥落的削弱系数,S为形状系数;
所述滚动体进入外圈故障区域的过程Io表述为:滚动体与轴承外圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000045
所述滚动体处于外圈故障区域的过程IIo表述为:滚动体与轴承外圈的剥落故障位置属于双边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000046
所述滚动体离开外圈故障区域的过程IIIo表述为:滚动体与轴承外圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000047
步骤1.1.2:根据轴承滚动体与轴承内圈的剥落故障位置将滚动体与轴承内圈的接触刚度变化分为三个内圈故障过程,分别为滚动体进入内圈故障区域的过程Ii、滚动体处于内圈故障区域的过程IIi、滚动体离开内圈故障区域的过程IIIi,利用公式(7)计算滚动体分别在所述三个内圈故障过程中的内圈修正接触刚度值Ki1
Figure BDA0002423206310000048
式中,Ki为滚动体与内圈之间的接触刚度,Ri为内圈剥落修正系数,
Figure BDA0002423206310000051
为滚动体进入内圈剥落边缘时的角位置,
Figure BDA0002423206310000052
为滚动体离开内圈剥落边缘时的角位置,
Figure BDA0002423206310000053
为内圈剥落中心的角位置;
所述滚动体进入内圈故障区域的过程Ii表述为:滚动体与轴承内圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000054
所述滚动体处于内圈故障区域的过程IIi表述为:滚动体与轴承内圈的剥落故障位置属于双边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000055
所述滚动体离开外圈故障区域的过程IIIi表述为:滚动体与轴承内圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000056
本发明的有益效果是:
本发明提供了一种滚动轴承剥落位置判定方法,1)为避免以往对轴承故障判别方法中,估算滚动体与故障接触时的接触变形量的方法,所造成的故障非线性动力学模型对实际故障模拟不准确的问题,将轴承滚动体经过内外圈故障分为三个过程,分别根据具体情况计算此过程中轴承的修正接触刚度;
2)同时为避免接触变形量的计算,考虑滚动体运行经过内外圈故障位置过程中的脉冲调制,可以更加精确地模拟出轴承发生早期故障时的运行状态,建立故障轴承非线性动力学模型,将所得模拟出的轴承运动曲线与实测故障或者疑似故障轴承的运动曲线进行比对,确定轴承是否发生故障以及发生故障的具体位置;
3)本专利可以更加精确地检测出发生早期剥落故障的轴承,避免了故障的进一步加剧,进而影响轴承甚至整机设备的正常运行,提高了检测效率,增加了检测精度,同时可以对轴承发生故障的位置进行判断,为轴承系统的维修和更换提供了帮助,根据统计轴承故障易发位置使得企业可以有针对性地采购轴承备件,完善设备的设计缺陷,培养工人良好的使用习惯,提高设备的维修与工作效率,降低设备的维护成本,提高了生产效率。
附图说明
图1为本发明中的滚动轴承剥落位置判定方法的流程图。
图2为本发明中的轴承内外圈存在剥落故障时滚动体经过过程示意图,其中图(a)为滚动体进入剥落故障区域时的示意图,与剥落故障位置单边缘接触,图(b)为滚动体完全进入剥落故障区域时的示意图,与剥落故障位置双边缘接触,图(c)为滚动体离开剥落故障区域时的示意图,与剥落故障位置单边缘接触。
图3为本发明中的剥落故障部位的脉冲序列图。
图4为本发明中的某一特定故障下的振动时域信号图,以及经过傅里叶变换后得到的振动频域信号图,其中图(a)表示振动时域信号图;图(b)表示对振动时域信号进行傅里叶变换得到的振动频域信号图。
图5为本发明中的实测频域振动曲线与模拟频域振动曲线比对图。
图6为本发明中的利用MATLAB进行仿真实验的具体编程流程图。
具体实施方式
下面是结合附图对本发明的技术方案进行详细说明。
如图1所示,一种滚动轴承剥落位置判定方法,包括如下步骤:
步骤1:建立剥落故障轴承的非线性动力学模型,得到剥落故障轴承的模拟频域振动曲线,包括如下步骤:
步骤1.1:首先采用集中参数模型模拟剥落故障轴承的运行状态,图2为轴承内外圈存在剥落故障时滚动体经过过程的示意图,然后根据轴承滚动体与轴承内/外圈的剥落故障位置得到轴承滚动体与轴承内/外圈的接触刚度变化过程,分别计算不同变化过程下的轴承内/外圈修正接触刚度值,具体表述为:
步骤1.1.1:根据轴承滚动体与轴承外圈的剥落故障位置将滚动体与轴承外圈的接触刚度变化分为三个外圈故障过程,分别为滚动体进入外圈故障区域的过程Io、滚动体处于外圈故障区域的过程IIo、滚动体离开外圈故障区域的过程IIIo,利用公式(6)计算滚动体在所述三个外圈故障过程中的外圈修正接触刚度值Ko1
Figure BDA0002423206310000061
式中,Ko为滚动体与外圈之间的接触刚度,Ro为外圈剥落修正系数,β为滚动体的角位置,
Figure BDA0002423206310000062
为滚动体进入外圈剥落边缘时的角位置,
Figure BDA0002423206310000063
为滚动体离开外圈剥落边缘时的角位置,
Figure BDA0002423206310000064
为外圈剥落中心的角位置,R为剥落的削弱系数,S为形状系数;
所述滚动体进入外圈故障区域的过程Io表述为:滚动体进入外圈剥落故障区域时,与轴承外圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000065
所述滚动体处于外圈故障区域的过程IIo表述为:滚动体处于外圈剥落故障区域内时,与轴承外圈的剥落故障位置属于双边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000071
所述滚动体离开外圈故障区域的过程IIIo表述为:滚动体离开外圈剥落故障区域时,与轴承外圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000072
步骤1.1.2:根据轴承滚动体与轴承内圈的剥落故障位置将滚动体与轴承内圈的接触刚度变化分为三个内圈故障过程,分别为滚动体进入内圈故障区域的过程Ii、滚动体处于内圈故障区域的过程IIi、滚动体离开内圈故障区域的过程IIIi,利用公式(7)计算滚动体在所述三个内圈故障过程中的内圈修正接触刚度值Ki1
Figure BDA0002423206310000073
式中,Ki为滚动体与内圈之间的接触刚度,Ri为内圈剥落修正系数,
Figure BDA0002423206310000074
为滚动体进入内圈剥落边缘时的角位置,
Figure BDA0002423206310000075
为滚动体离开内圈剥落边缘时的角位置,
Figure BDA0002423206310000076
为内圈剥落中心的角位置;
所述滚动体进入内圈故障区域的过程Ii表述为:滚动体进入内圈剥落故障区域时,与轴承内圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000077
所述滚动体处于内圈故障区域的过程IIi表述为:滚动体处于内圈剥落故障区域内时,与轴承内圈的剥落故障位置属于双边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000078
所述滚动体离开外圈故障区域的过程IIIi表述为:滚动体离开内圈剥落故障区域时,与轴承内圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure BDA0002423206310000079
通过修改不同状态下的内外圈修正接触刚度,可以模拟内外圈的时变刚度,使载荷-变形关系更符合滚动体与内外圈接触的实际情况,也使建模更加真实。
步骤1.2:利用公式(1)计算剥落故障轴承的总接触刚度值K;
Figure BDA00024232063100000710
式中,Ko1为滚动体与轴承外圈接触时的外圈修正接触刚度值,Ki1为滚动体与轴承内圈接触时的内圈修正接触刚度值;
假定在单位载荷作用下,轴承滚动体与故障部位接触产生的冲击力为矩形脉冲,因此,这些伴随着轴承运动的冲击序列可以在一个周期内用无穷大的等幅类矩形函数来表示,如图3所示,当脉冲间隔时间较短时,脉冲信号可视为梯形信号,其中d(t)是轴承内圈在轴承载荷方向上的最大位移,如果轴承外圈固定,内圈转轴旋转不同步,当内外圈产生剥落故障时,Td和Tv分别表示为:
Figure BDA0002423206310000081
Figure BDA0002423206310000082
Figure BDA0002423206310000083
Figure BDA0002423206310000084
式中,n为轴承转速,Td为矩形脉冲序列的周期,Tv为滚动体穿过剥落故障区域的时间,L为剥落宽度,fbpfo为轴承外圈发生剥落故障时的故障特征频率,fbpfi为轴承内圈发生剥落故障时的故障特征频率,Z为滚动体个数,f为轴承转动频率,d为轴承内径,D为节圆直径,α为轴承接触角;
步骤1.3:当轴承高速旋转时,故障位置处的振动信号随轴承的运动而周期性变化,因此滚动体经过载荷区域时产生的冲击力受载荷分布的影响,在脉冲间隔很短的情况下,利用公式(2)计算剥落故障轴承高速旋转时的调制脉冲值ffault
Figure BDA0002423206310000085
式中,d(t)为轴承内圈在轴承载荷方向上的最大位移,Qmax为轴承所受的最大载荷,δmax为最大接触变形量,c为相对间隙,f为轴承转动频率,t为轴承运行时间;
步骤1.4:建立剥落故障轴承的非线性动力学模型如公式(3)所示,
Figure BDA0002423206310000086
式中,m为轴承质量,C为轴承系统阻尼,F为轴承系统所受外载荷;
步骤1.5:根据公式(3)建立剥落故障轴承的仿真模型,利用MATLAB对剥落轴承故障的仿真模型进行仿真运算,得到剥落故障轴承模拟运行时的时域振动信号,根据模拟运行时的时域振动信号得到模拟时域振动曲线,对模拟运行时得到的时域振动信号进行傅里叶变换,得到剥落故障轴承模拟运行时的模拟频域振动曲线,其中当剥落故障位置设置在正下方的六点钟方向,故障大小设置为0.6mm时,利用公式(3)建立该特定剥落故障下的仿真模型,然后利用MATLAB仿真实验,得到的具体模拟时域振动曲线图,以及经过傅里叶变换得到的具体模拟频域振动曲线图,如图4所示;
步骤2:提取实测轴承的实测频域振动曲线,在实测轴承的外圈任意点处设置接触式振动传感器,本实施例中采用的接触式传感器的型号为GB/T 6075,所述实测轴承包括实际需要检测的故障轴承或者疑似故障轴承,当实测轴承运行时,接触式振动传感器采集实测轴承外圈上的振动信号,采集到的振动信号通过数据采集装置处理后得到实测振动时域信号,对实测振动时域信号进行傅里叶变换(简称FFT),得到实测轴承的实测频域振动曲线,本实施例中采用的采集装置具体型号为AIC9000,接触式传感器与信号采集装置连接,信号采集装置与电脑连接,接触式传感器将采集到的信号传输给信号采集装置,经由信号采集装置处理后生成轴承振动的时域数据在电脑上显示;
步骤3:利用公式(4)、(5)分别计算出实测轴承外圈、内圈的故障特征频率,以模拟频域振动曲线作为参照与实测频域振动曲线进行对比,如图5所示,分别查看故障特征频率所对应的峰值,若两组曲线中存在峰值所对应的故障特征频率相差不大于Δf,则定义相差不大于Δf的故障特征频率为相似故障特征频率,Δf的选取方法如下:
Figure BDA0002423206310000091
检查所述相似故障特征频率是否为轴承内圈或者外圈的故障特征频率,若所述相似故障特征频率为内圈的故障特征频率或者外圈的故障特征频率,则判断实测轴承存在剥落故障,若所述相似故障特征频率既不是轴承内圈的故障特征频率也不是轴承外圈的故障特征频率,则判断实测轴承不存在剥落故障,如果所述相似故障特征频率存在于实测轴承的内圈,则故障位置为轴承内圈,如果所述相似故障特征频率存在于实测轴承的外圈,则故障位置为轴承的外圈。
本实施例中利用MATLAB建立仿真模型,并进行仿真运算的具体过程如图6所示,操作示例如下:
S1:登录MATLAB界面,打开文件夹中的baochijia.m文件,M文件中,需要填写的为“输入参数”一栏,其中的“轴承转速r1、陶瓷球半径r、陶瓷球个数Z、外载荷Fx、接触角alpha、节圆直径D_m”为需要填写的轴承相关参数,此外,还需要将实验数据文件至M文件中,以便调用其中的试验数据对保持架磨损尺寸进行验证;
S2:填写“输入参数”,选取7009AC角接触陶瓷球轴承为检测轴承:轴承转速r1=250、陶瓷滚动体半径r=4.75e-3、陶瓷滚动体个数Z=17、外载荷Fx=500、接触角alpha=40/180*pi、节圆直径D_m=60e-3;
S3:导入实验数据,选取三组实测数据作为实验数据进行分析,三组实测数据存入Excel中,命名为“检测数据”,将Excel文件“检测数据”拖拽至MATLAB界面左侧“当前文件夹”一栏中,并在M文件中“导入试验数据”一栏T=xlsread(”,”)中,第一个‘’内填写需要导入的文件名称+.xlsx,即测试数据.xlsx,T=xlsread(”,”)一栏中的第二个‘’内填写需要选取的被检测数据,被检测数据中一共有三组实际测量的数据,分别位于第A、B和C列,每组数据包含300个参数,在MATLAB中调取这三组数据的代码分别为A2:A301、B2:B301、C2:C301,首先选取第一组数据,即将代码A2:A301填写进T=xlsread(”,”)一栏中的第二个‘’内;
S4:运行剥落故障位置判定MATLAB程序,点击MATLAB系统界面“编辑器”一栏中“运行”按键,在MATLAB系统下方“命令行窗口”可以得到对此组检测数据的检验结果,可以看到此组数据检验结果为“未发生剥落故障”。

Claims (1)

1.一种滚动轴承剥落位置判定方法,其特征在于,包括如下步骤:
步骤1:建立剥落故障轴承的非线性动力学模型,得到剥落故障轴承的模拟频域振动曲线,包括如下步骤:
步骤1.1:首先采用集中参数模型模拟剥落故障轴承的运行状态,然后根据轴承滚动体与轴承内/外圈的剥落故障位置得到轴承滚动体与轴承内/外圈的接触刚度变化过程,分别计算不同变化过程下的轴承内/外圈修正接触刚度值;
步骤1.2:利用公式(1)计算剥落故障轴承的总接触刚度值K;
Figure FDA0003315839500000011
式中,Ko1为滚动体与轴承外圈接触时的外圈修正接触刚度值,Ki1为滚动体与轴承内圈接触时的内圈修正接触刚度值;
步骤1.3:利用公式(2)计算剥落故障轴承高速旋转时的调制脉冲值ffault
Figure FDA0003315839500000012
式中,d(t)为轴承内圈在轴承载荷方向上的最大位移,Qmax为轴承所受的最大载荷,δmax为最大接触变形量,c为相对间隙,f为轴承转动频率,t为轴承运行时间;
步骤1.4:建立剥落故障轴承的非线性动力学模型如公式(3)所示,
Figure FDA0003315839500000013
式中,m为轴承质量,C为轴承系统阻尼,F为轴承系统所受外载荷;
步骤1.5:根据公式(3)建立剥落故障轴承的仿真模型,并进行仿真运算,得到剥落故障轴承模拟运行时的时域振动信号,根据模拟运行时的时域振动信号得到模拟时域振动曲线,对模拟运行时得到的时域振动信号进行傅里叶变换,得到剥落故障轴承模拟运行时的模拟频域振动曲线;
步骤2:提取实测轴承的实测频域振动曲线,在实测轴承的外圈任意点处设置接触式振动传感器,所述实测轴承包括实际需要检测的故障轴承或者疑似故障轴承,当实测轴承运行时,接触式振动传感器采集实测轴承外圈上的振动信号,采集到的振动信号通过数据采集装置处理后得到实测振动时域信号,对实测振动时域信号进行傅里叶变换,得到实测轴承的实测频域振动曲线;
步骤3:利用公式(4)、(5)分别计算出实测轴承外圈、内圈的故障特征频率,以模拟频域振动曲线作为参照与实测频域振动曲线进行对比,分别查看故障特征频率所对应的峰值,若两组曲线中存在峰值所对应的故障特征频率相差不大于Δf,则定义相差不大于Δf的故障特征频率为相似故障特征频率,检查所述相似故障特征频率是否为轴承内圈或者外圈的故障特征频率,若所述相似故障特征频率为内圈的故障特征频率或者外圈的故障特征频率,则判断实测轴承存在剥落故障,若所述相似故障特征频率既不是轴承内圈的故障特征频率也不是轴承外圈的故障特征频率,则判断实测轴承不存在剥落故障,如果所述相似故障特征频率存在于实测轴承的内圈,则故障位置为轴承内圈,如果所述相似故障特征频率存在于实测轴承的外圈,则故障位置为轴承的外圈;
Figure FDA0003315839500000021
Figure FDA0003315839500000022
式中,fbpfo为轴承外圈发生剥落故障时的故障特征频率,fbpfi为轴承内圈发生剥落故障时的故障特征频率,Z为滚动体个数,f为轴承转动频率,d为轴承内径,D为节圆直径,α为轴承接触角;
所述的步骤1.1具体表述为:
步骤1.1.1:根据轴承滚动体与轴承外圈的剥落故障位置将滚动体与轴承外圈的接触刚度变化分为三个外圈故障过程,分别为滚动体进入外圈故障区域的过程Io、滚动体处于外圈故障区域的过程IIo、滚动体离开外圈故障区域的过程IIIo,利用公式(6)计算滚动体分别在所述三个外圈故障过程中的外圈修正接触刚度值Ko1
Figure FDA0003315839500000023
式中,Ko为滚动体与外圈之间的接触刚度,Ro为外圈剥落修正系数,β为滚动体的角位置,
Figure FDA0003315839500000031
为滚动体进入外圈剥落边缘时的角位置,
Figure FDA0003315839500000032
为滚动体离开外圈剥落边缘时的角位置,
Figure FDA0003315839500000033
为外圈剥落中心的角位置,R为剥落的削弱系数,S为形状系数;
所述滚动体进入外圈故障区域的过程Io表述为:滚动体与轴承外圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure FDA0003315839500000034
所述滚动体处于外圈故障区域的过程IIo表述为:滚动体与轴承外圈的剥落故障位置属于双边缘接触,滚动体的角位置β满足
Figure FDA0003315839500000035
所述滚动体离开外圈故障区域的过程IIIo表述为:滚动体与轴承外圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure FDA0003315839500000036
步骤1.1.2:根据轴承滚动体与轴承内圈的剥落故障位置将滚动体与轴承内圈的接触刚度变化分为三个内圈故障过程,分别为滚动体进入内圈故障区域的过程Ii、滚动体处于内圈故障区域的过程IIi、滚动体离开内圈故障区域的过程IIIi,利用公式(7)计算滚动体分别在所述三个内圈故障过程中的内圈修正接触刚度值Ki1
Figure FDA0003315839500000037
式中,Ki为滚动体与内圈之间的接触刚度,Ri为内圈剥落修正系数,
Figure FDA0003315839500000038
为滚动体进入内圈剥落边缘时的角位置,
Figure FDA0003315839500000039
为滚动体离开内圈剥落边缘时的角位置,
Figure FDA00033158395000000310
为内圈剥落中心的角位置;
所述滚动体进入内圈故障区域的过程Ii表述为:滚动体与轴承内圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure FDA00033158395000000311
所述滚动体处于内圈故障区域的过程IIi表述为:滚动体与轴承内圈的剥落故障位置属于双边缘接触,滚动体的角位置β满足
Figure FDA00033158395000000312
所述滚动体离开外圈故障区域的过程IIIi表述为:滚动体与轴承内圈的剥落故障位置属于单边缘接触,滚动体的角位置β满足
Figure FDA0003315839500000041
CN202010212170.1A 2020-03-24 2020-03-24 一种滚动轴承剥落位置判定方法 Active CN111337255B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010212170.1A CN111337255B (zh) 2020-03-24 2020-03-24 一种滚动轴承剥落位置判定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010212170.1A CN111337255B (zh) 2020-03-24 2020-03-24 一种滚动轴承剥落位置判定方法

Publications (2)

Publication Number Publication Date
CN111337255A CN111337255A (zh) 2020-06-26
CN111337255B true CN111337255B (zh) 2022-02-01

Family

ID=71180368

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010212170.1A Active CN111337255B (zh) 2020-03-24 2020-03-24 一种滚动轴承剥落位置判定方法

Country Status (1)

Country Link
CN (1) CN111337255B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114777903B (zh) * 2022-04-27 2023-06-27 浙江大学 一种基于信号互相关的光频域反射计多点振动检测方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN207317932U (zh) * 2017-11-10 2018-05-04 孙群丽 一种轴承检测装置
CN108801638A (zh) * 2018-06-19 2018-11-13 武汉理工大学 一种滚动轴承剥落尺寸的估计方法
CN110132597A (zh) * 2019-04-29 2019-08-16 清华大学 一种轴承内滚道剥落宽度的测算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN207317932U (zh) * 2017-11-10 2018-05-04 孙群丽 一种轴承检测装置
CN108801638A (zh) * 2018-06-19 2018-11-13 武汉理工大学 一种滚动轴承剥落尺寸的估计方法
CN110132597A (zh) * 2019-04-29 2019-08-16 清华大学 一种轴承内滚道剥落宽度的测算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A Theoretical Model with the Effect of Cracks in the Local Spalling of Full Ceramic Ball Bearings;Huaitao Shi et al.;《applied sciences》;20191003;第9卷;正文第1-21页 *
滚动轴承外圈剥落故障双冲击特征机理建模;刘倩楠 等;《振动工程学报》;20170831;第30卷(第4期);全文 *
滚动轴承系统局部缺陷位移激励动力学建模及振动响应分析;康建雄;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20180215;正文第1、9-11、16-17、33-41页 *

Also Published As

Publication number Publication date
CN111337255A (zh) 2020-06-26

Similar Documents

Publication Publication Date Title
CN109822399B (zh) 基于并行深度神经网络的数控机床刀具磨损状态预测方法
Cui et al. HVSRMS localization formula and localization law: Localization diagnosis of a ball bearing outer ring fault
CN111947928B (zh) 一种多源信息融合的轴承故障预测系统及方法
CN108181105B (zh) 基于逻辑回归和j散度的滚动轴承故障预诊方法及系统
CN111922095A (zh) 一种冷轧机轧辊扭振故障异常的振动诊断方法
CN111946559B (zh) 一种风机基础和塔架结构检测方法
CN108896299A (zh) 一种齿轮箱故障检测方法
CN105067262A (zh) 一种滚动轴承状态监测方法
CN113027703A (zh) 包括传动系统、齿轮箱和发电机的旋转机械的损害和剩余使用寿命的确定
CN105004523A (zh) 基于加权相似性测度的滚动轴承状态监测方法
KR102040179B1 (ko) 제조 설비의 이상 감지 및 진단 방법
CN113188794B (zh) 一种基于改进pso-bp神经网络齿轮箱故障诊断方法及装置
CN112160761B (zh) 一种基于现场参数的硬岩tbm刀盘滚刀磨损实时评估方法
CN110046409B (zh) 一种基于ResNet的汽轮机部件健康状态评估方法
CN111337255B (zh) 一种滚动轴承剥落位置判定方法
Mufazzal et al. Theoretical and experimental analyses of vibration impulses and their influence on accurate diagnosis of ball bearing with localized outer race defect
CN112131797B (zh) 一种基于应力分析的主轴轴承寿命预测和可靠性评估方法
CN115023598A (zh) 用于确定轴承的剩余可用寿命的方法
CN117629615B (zh) 一种盾构机故障检测方法、设备及介质
CN110646138A (zh) 一种旋转机械无键相无试重动平衡方法和分析装置
JP3108405B2 (ja) 機器の診断方法
EP1896809B1 (de) Verfahren und vorrichtung zur zustandsüberwachung von schiffsanlagen
Hamadache et al. Improving signal-to-noise ratio (SNR) for inchoate fault detection based on principal component analysis (PCA)
CN112621381B (zh) 机床进给系统健康状态智能评估方法及装置
Roemer et al. Machine health monitoring and life management using finite-element-based neural networks

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant