CN107741324B - 一种滚动轴承外圈故障定位诊断方法 - Google Patents
一种滚动轴承外圈故障定位诊断方法 Download PDFInfo
- Publication number
- CN107741324B CN107741324B CN201710954465.4A CN201710954465A CN107741324B CN 107741324 B CN107741324 B CN 107741324B CN 201710954465 A CN201710954465 A CN 201710954465A CN 107741324 B CN107741324 B CN 107741324B
- Authority
- CN
- China
- Prior art keywords
- outer ring
- angle position
- failure
- vibration acceleration
- bearing
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic 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
本发明涉及一种滚动轴承外圈故障定位诊断方法。该方法首先建立滚动轴承***非线性动力学模型,并由第二类拉格朗日方程推导出轴承***的动力学微分方程组。基于***动力学微分方程组,仿真出不同故障角位置时的轴承***水平方向和垂直方向振动加速度信号。然后由仿真数据,采用最小二乘法拟合出外圈故障角位置与平垂同步均方根值的线性关系式。其次测取待诊断轴承***的定位诊断基准数据,并计算定位诊断基准数据的平垂同步均方根值,从而确定待诊断轴承***的外圈故障角位置θ与平垂同步均方根值的关系式。最终由外圈故障角位置定位法则确诊出被诊断外圈故障的角位置θd。本发明实现轴承外圈故障的精确定位诊断,具有重要的理论意义和应用价值。
Description
技术领域
本发明属于滚动轴承故障诊断领域,具体涉及一种滚动轴承外圈故障定位诊断方法。
背景技术
滚动轴承是旋转机械等设备中使用最多的零部件之一,也是最容易损坏的零部件。当其产生故障时,将诱发整个设备产生故障,因此,滚动轴承故障诊断研究对于维护设备的正常运转非常关键。目前轴承故障诊断研究主要集中在轴承故障定量诊断、轴承故障模式识别等方面。对于轴承外圈故障的定位诊断研究几乎没有,然而不同角位置的故障可能对应着不同的产生原因,比如制造、装配以及工作状态等。另外,在其它因素相同的情况下,不同角位置的故障可能对应着不同的剩余寿命,比如缺陷角位置越靠近承载中心的轴承的剩余寿命往往会越短。再有,实现轴承外圈故障定位诊断将有助于拆机后更快速的排查出故障位置,同时有助于轴承的维护以及维修工作,提高工作效率。因此解决如何快速、准确地预测出轴承外圈故障角位置的问题对轴承故障诊断以及轴承剩余寿命估计具有重要的理论意义和工程应用价值。
VR是指平-垂同步均方根值,是由发明人为实现轴承外圈的定位诊断而提出的一种新指标,其详细的算法和应用技术将在本发明内容和实施例中具体说明。
发明内容
为了快速、准确地预测出滚动轴承外圈故障的角位置,本发明基于机理研究所挖掘出的VR指标提出了一种滚动轴承外圈故障定位诊断方法,同时为轴承故障技术开拓了新的发展方向。
为实现上述目的,本发明的技术方案如下:
一种滚动轴承外圈故障定位诊断方法,包括以下具体步骤:
步骤1建立滚动轴承***非线性动力学模型及其微分方程组,
建立轴承***全局绝对直角坐标系,以面向动力输入端为参照,将水平方向设定为x轴,且水平向右为正方向,将垂直方向设定为y轴,且垂直向下为正方向,轴承内圈表示为质点m1,轴承外圈表示为质点m2,轴承内圈的水平方向位移和垂直方向位移分别表示为x1和y1,轴承外圈的水平方向位移和垂直方向位移分别表示为x2和y2,轴承外圈由水平方向的弹簧阻尼器1和垂直方向的弹簧阻尼器2支撑,弹簧阻尼器1的刚度和阻尼系数分别表示为k1和c1,弹簧阻尼器2的刚度和阻尼系数分别表示为k2和c2,将外圈故障角位置表示为θ,基于上述所建模型,由第二类拉格朗日方程推导出轴承***的动力学微分方程组;
步骤2建立平垂同步均方根值与外圈故障角位置之间的函数关系,设定外圈故障角位置的变化区间为[240°,300°],将该变化区间进行7等分,则得到7个由小到大的外圈故障角位置θ1,θ2…θj…θ7,j表示外圈故障角位置的编号,取值1~7之间的整数;由步骤1中推导的轴承***动力学微分方程组求解每个故障角位置θj对应的水平方向仿真振动加速度信号ax,j和垂直方向仿真振动加速度信号ay,j,再计算每个故障角位置对应的平垂同步均方根值:
上式中,N为仿真数据点数,其中,以θj为自变量,VRj为因变量,由最小二乘法线性拟合出外圈故障角位置θ与VR值的关系式:
θ=270±k·(VR-VR4) (2)
步骤3测取待诊断轴承***的定位诊断基准数据,
以面向动力输入端为参照,设定水平向右为0°,将振动加速度传感器1安装在待诊断轴承***的轴承座的90°位置,将振动加速度传感器2安装在待诊断轴承***的轴承座的180°位置,由振动加速度传感器1和振动加速度传感器2分别测取外圈故障角位置在270°的垂直方向振动加速度信号a′y和水平方向振动加速度信号a′x,再计算待诊断轴承***的定位诊断基准VR值:
将VR'替换公式(2)中的VR4,得到待诊断轴承***的外圈故障角位置θ与VR值的关系式:
θ=270±k·(VR-VR′) (4)
步骤3外圈故障角位置定位法则,
以面向动力输入端为参照,设定水平向右为0°,将振动加速度传感器1安装在待诊断轴承***的轴承座的90°位置,将振动加速度传感器2安装在待诊断轴承***的轴承座的180°位置,由振动加速度传感器1和振动加速度传感器2分别测取被诊断故障轴承的垂直方向振动加速度信号a″y和水平方向振动加速度信号a″x,并计算被诊断故障轴承的VR值:
分别绘制a″x和a″y的时域波形图,观察波形图中的故障冲击起始方向,当a″x波形图中的故障冲击起始方向为正且a″y波形图中的故障冲击起始方向为负时,则该外圈故障可判定为位于270°的左侧,即被诊断的故障角位置可确定为:
θd=270-k·(VR″-VR′) (6)
当a″x波形图中的故障冲击起始方向为负且a″y波形图中的故障冲击起始方向为负时,则该外圈故障可判定为位于270°的右侧,即被诊断的故障角位置可确定为:
θd=270+k·(VR″-VR′) (7)
本发明的有益效果是:提出了一种滚动轴承外圈故障定位诊断方法,在建立待诊断轴承***的外圈故障角位置θ与VR值的函数关系式之后,由外圈故障角位置定位法则就能确诊出被诊断外圈故障的角位置,有重要的实际应用价值。
下面具体结合附图与实例对本发明作进一步的说明。
附图说明
图1本发明的工作流程图
图2滚动轴承***非线性动力学模型示意图
图3外圈故障角位置θ与平垂同步均方根值的变化关系图
图4为待诊断故障轴承***的水平方向和垂直方向振动加速度信号
具体实施方式
下面具体结合附图与实施例对本发明的诊断方法进行详细说明,但是本发明的保护范围不局限于所述实施例。
如图1所示,是本发明的一种滚动轴承外圈故障定位诊断方法的工作流程图。具体实施过程如下:
步骤1建立滚动轴承***非线性动力学模型和其微分方程组,
建立轴承***全局绝对直角坐标系,以面向动力输入端为参照,将水平方向设定为x轴,且水平向右为正方向,将垂直方向设定为y轴,且垂直向下为正方向,轴承内圈表示为质点m1,轴承外圈表示为质点m2,轴承内圈的水平方向位移和垂直方向位移分别表示为x1和y1,轴承外圈的水平方向位移和垂直方向位移分别表示为x2和y2,轴承外圈由水平方向的弹簧阻尼器1和垂直方向的弹簧阻尼器2支撑,弹簧阻尼器1的刚度和阻尼系数分别表示为k1和c1,弹簧阻尼器2的刚度和阻尼系数分别表示为k2和c2,将外圈故障角位置表示为θ,完成后的轴承***非线性动力学模型如图2所示;基于上述所建模型,由第二类拉格朗日方程推导出轴承***的动力学微分方程组;
步骤2建立平垂同步均方根值与外圈故障角位置之间的函数关系,本发明采用NSK6308轴承作为具体实施例,轴承滚珠数量Nb=8,滚珠半径rb=7.54mm,外滚道半径Ro=40.04mm,选取的故障角跨度为Δφf=1°,故障深度h=0.04mm,采样频率为Fs=65536Hz,设定外圈故障角位置的变化区间为[240°,300°],将该变化区间进行7等分,则得到7个由小到大的外圈故障角位置θ1,θ2…θj…θ7,j表示外圈故障角位置的编号,取值1~7之间的整数,由步骤1中推导的轴承***动力学微分方程组求解每个故障角位置θj对应的水平方向仿真振动加速度信号ax,j和垂直方向仿真振动加速度信号ay,j,再计算每个故障角位置对应的平垂同步均方根值:
根据公式(1)和已知参数,可计算出每个故障角位置θj对应的VR值分别为VR1=0.5649,VR2=0.3384,VR3=0.1648,VR4=0.01165,VR5=0.1638,VR6=0.3615,VR7=0.5943。以θj为自变量,VRj为因变量,由最小二乘法线性拟合出外圈故障角位置θ与VR值的关系式:
由公式2确定的本实施例中外圈故障角位置θ与VR值的变化曲线如图3所示。
步骤3测取待诊断轴承***的定位诊断基准数据,
以面向动力输入端为参照,设定水平向右为0°,将振动加速度传感器1安装在待诊断轴承***的轴承座的90°位置,将振动加速度传感器2安装在待诊断轴承***的轴承座的180°位置,由振动加速度传感器1和振动加速度传感器2分别测取外圈故障角位置在270°的垂直方向振动加速度信号a′y和水平方向振动加速度信号a′x,再计算待诊断轴承***的定位诊断基准VR值:
将VR'=0.497替换公式(2)中的VR4,得到待诊断轴承***的外圈故障角位置θ与VR值的关系式为:
θ=270±50·(VR-0.497) (4)
步骤3外圈故障角位置定位法则,
以面向动力输入端为参照,设定水平向右为0°,将振动加速度传感器1安装在待诊断轴承***的轴承座的90°位置,将振动加速度传感器2安装在待诊断轴承***的轴承座的180°位置,由振动加速度传感器1和振动加速度传感器2分别测取被诊断故障轴承的垂直方向振动加速度信号a″y和水平方向振动加速度信号a″x,并计算被诊断故障轴承的VR值:
分别绘制a″x和a″y的时域波形图,如图4所示,观察波形图中的故障冲击起始方向可知,a″x中的故障冲击起始方向为负且a″y中的故障冲击起始方向也为负,所以本实施例中的外圈故障可判定为位于270°的右侧,即被诊断的故障角位置可确诊为:
由此实施例可知,在轴承***的物理和几何参数给定之后,由本发明方法提出的公式只需VR这一参量就可预测出对应的故障角位置。充分说明了本发明方法的重要意义和应用价值。
Claims (1)
1.一种滚动轴承外圈故障定位诊断方法,其特征在于:该方法包括以下具体步骤,
步骤1建立滚动轴承***非线性动力学模型及其微分方程组,
建立轴承***全局绝对直角坐标系,以面向动力输入端为参照,将水平方向设定为x轴,且水平向右为正方向,将垂直方向设定为y轴,且垂直向下为正方向,轴承内圈表示为质点m1,轴承外圈表示为质点m2,轴承内圈的水平方向位移和垂直方向位移分别表示为x1和y1,轴承外圈的水平方向位移和垂直方向位移分别表示为x2和y2,轴承外圈由水平方向的第一弹簧阻尼器和垂直方向的第二弹簧阻尼器支撑,第一弹簧阻尼器的刚度和阻尼系数分别表示为k1和c1,第二弹簧阻尼器的刚度和阻尼系数分别表示为k2和c2,将外圈故障角位置表示为θ,基于上述所建模型,由第二类拉格朗日方程推导出轴承***的动力学微分方程组;
步骤2建立平垂同步均方根值与外圈故障角位置之间的函数关系,
设定外圈故障角位置的变化区间为[240°,300°],将该变化区间进行7等分,则得到7个由小到大的外圈故障角位置θ1,θ2…θj…θ7,j表示外圈故障角位置的编号,取值1~7之间的整数;由步骤1中推导的轴承***动力学微分方程组求解每个故障角位置θj对应的水平方向仿真振动加速度信号ax,j和垂直方向仿真振动加速度信号ay,j,再计算每个故障角位置对应的平垂同步均方根值:
上式中,N为仿真数据点数,其中,以θj为自变量,VRj为因变量,由最小二乘法线性拟合出外圈故障角位置θ与VR值的关系式:
θ=270±k·(VR-VR4) (2)
步骤3测取待诊断轴承***的定位诊断基准数据,
以面向动力输入端为参照,设定水平向右为0°,将第一振动加速度传感器安装在待诊断轴承***的轴承座的90°位置,将第二振动加速度传感器安装在待诊断轴承***的轴承座的180°位置,由第一振动加速度传感器和第二振动加速度传感器分别测取外圈故障角位置在270°的垂直方向振动加速度信号a′y和水平方向振动加速度信号a′x,再计算待诊断轴承***的定位诊断基准VR值:
将VR'替换公式(2)中的VR4,得到待诊断轴承***的外圈故障角位置θ与VR值的关系式:
θ=270±k·(VR-VR′) (4)
步骤4外圈故障角位置定位法则,
以面向动力输入端为参照,设定水平向右为0°,将第一振动加速度传感器安装在待诊断轴承***的轴承座的90°位置,将第二振动加速度传感器安装在待诊断轴承***的轴承座的180°位置,由第一振动加速度传感器和第二振动加速度传感器分别测取被诊断故障轴承的垂直方向振动加速度信号a″y和水平方向振动加速度信号a″x,并计算被诊断故障轴承的VR值:
分别绘制a″x和a″y的时域波形图,观察波形图中的故障冲击起始方向,当a″x波形图中的故障冲击起始方向为正且a″y波形图中的故障冲击起始方向为负时,则该外圈故障可判定为位于270°的左侧,即被诊断的故障角位置可确定为:
θd=270-k·(VR″-VR′) (6)
当a″x波形图中的故障冲击起始方向为负且a″y波形图中的故障冲击起始方向为负时,则该外圈故障可判定为位于270°的右侧,即被诊断的故障角位置可确定为:
θd=270+k·(VR″-VR′) (7)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710954465.4A CN107741324B (zh) | 2017-10-13 | 2017-10-13 | 一种滚动轴承外圈故障定位诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710954465.4A CN107741324B (zh) | 2017-10-13 | 2017-10-13 | 一种滚动轴承外圈故障定位诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107741324A CN107741324A (zh) | 2018-02-27 |
CN107741324B true CN107741324B (zh) | 2019-05-24 |
Family
ID=61236602
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710954465.4A Expired - Fee Related CN107741324B (zh) | 2017-10-13 | 2017-10-13 | 一种滚动轴承外圈故障定位诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107741324B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109211566B (zh) * | 2018-08-21 | 2020-03-13 | 北京工业大学 | 一种滚动轴承外圈缺陷二维量化诊断方法 |
CN109992926B (zh) * | 2019-04-23 | 2021-02-05 | 清华大学 | 轴承外圈缺陷角位置量化评估方法 |
CN110146290B (zh) * | 2019-05-29 | 2020-08-18 | 清华大学 | 一种基于定位普遍法则的滚动轴承外圈缺陷定位诊断方法 |
CN111046502B (zh) * | 2019-11-13 | 2023-08-25 | 长江大学 | 一种管道穿越断层的土弹簧刚度计算方法及装置 |
CN112267979B (zh) * | 2020-10-26 | 2021-07-23 | 积成电子股份有限公司 | 一种判断偏航轴承失效的预警方法及*** |
CN113029570B (zh) * | 2021-04-01 | 2022-04-26 | 温州大学 | 一种谐波轴承故障样本生成模型与诊断方法 |
CN114235414B (zh) * | 2021-12-28 | 2023-08-04 | 频率探索智能科技江苏有限公司 | 适用于外圈缺陷定位诊断的信号处理方法 |
CN114486261B (zh) * | 2022-01-21 | 2022-09-23 | 沈阳科网通信息技术有限公司 | 基于同步均方根之比的轴承外圈故障定位方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102009033137A1 (de) * | 2009-07-15 | 2011-02-03 | Schaeffler Technologies Gmbh & Co. Kg | Linearführung und Verfahren zum Betreiben eines Mess- und Auswertesystems |
CN102183366A (zh) * | 2011-03-08 | 2011-09-14 | 上海大学 | 滚动轴承振动测量和故障分析装置及方法 |
CN205067090U (zh) * | 2015-11-11 | 2016-03-02 | 内蒙古科技大学 | 滚动轴承故障检测与诊断*** |
CN106198008A (zh) * | 2016-06-16 | 2016-12-07 | 温州大学 | 一种电机轴承外圈故障辨识方法 |
CN106295072A (zh) * | 2016-08-30 | 2017-01-04 | 北京工业大学 | 一种针对轴承内外圈故障的定量趋势诊断方法 |
CN106769030A (zh) * | 2016-11-10 | 2017-05-31 | 浙江工业大学 | 一种基于mea‑bp神经网络算法的轴承状态跟踪与预测方法 |
-
2017
- 2017-10-13 CN CN201710954465.4A patent/CN107741324B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102009033137A1 (de) * | 2009-07-15 | 2011-02-03 | Schaeffler Technologies Gmbh & Co. Kg | Linearführung und Verfahren zum Betreiben eines Mess- und Auswertesystems |
CN102183366A (zh) * | 2011-03-08 | 2011-09-14 | 上海大学 | 滚动轴承振动测量和故障分析装置及方法 |
CN205067090U (zh) * | 2015-11-11 | 2016-03-02 | 内蒙古科技大学 | 滚动轴承故障检测与诊断*** |
CN106198008A (zh) * | 2016-06-16 | 2016-12-07 | 温州大学 | 一种电机轴承外圈故障辨识方法 |
CN106295072A (zh) * | 2016-08-30 | 2017-01-04 | 北京工业大学 | 一种针对轴承内外圈故障的定量趋势诊断方法 |
CN106769030A (zh) * | 2016-11-10 | 2017-05-31 | 浙江工业大学 | 一种基于mea‑bp神经网络算法的轴承状态跟踪与预测方法 |
Non-Patent Citations (2)
Title |
---|
基于振动相应机理的轴承故障定量诊断及量化分析;崔玲丽;《北京工业大学学报》;20151130;第41卷(第11期);1681-1687 |
滚动轴承缺陷频谱特征及故障诊断实例;李丽红;《通用机械》;20141231;48-51 |
Also Published As
Publication number | Publication date |
---|---|
CN107741324A (zh) | 2018-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107741324B (zh) | 一种滚动轴承外圈故障定位诊断方法 | |
CN107704695B (zh) | 一种滚动轴承外圈缺陷全尺寸定量诊断方法 | |
CN105928707B (zh) | 一种滚动轴承-转子***动态耦合建模方法 | |
CN110376882A (zh) | 基于有限时间扩张状态观测器的预定性能控制方法 | |
CN105716603B (zh) | 根据视频中物体的太阳影子变化对物体进行定位的方法 | |
CN107505841A (zh) | 一种基于干扰估计器的机械臂姿态鲁棒控制方法 | |
CN109765026A (zh) | 一种低速动态试验洞壁干扰修正方法 | |
CN107843428A (zh) | 轴承载荷在线获取方法及装置、轴承寿命评估方法及装置 | |
CN110542439A (zh) | 基于三维气浮的惯性器件残余力矩测量装置及方法 | |
CN107870371A (zh) | 一种动基座重力梯度仪自梯度补偿方法 | |
CN205175667U (zh) | 用于研究密封动力特性的转子试验台 | |
CN108843654A (zh) | 一种基于子空间辨识的阀控液压缸泄漏判断装置及方法 | |
CN110146290B (zh) | 一种基于定位普遍法则的滚动轴承外圈缺陷定位诊断方法 | |
CN102853848A (zh) | 基于捷联惯导***定位精度的惯性器件误差仿真方法 | |
CN112115575A (zh) | 一种设备润滑油状态评价***和方法 | |
CN109145501B (zh) | 一种航空发动机中介轴承局部损伤故障振动仿真方法 | |
CN206195966U (zh) | 一种钢结构网架监测*** | |
Johansen | On developing a digital twin for fault detection in drivetrains of offshore wind turbines | |
Ruan et al. | Modeling research of MEMS gyro drift Based on Kalman filter | |
Choux et al. | Leakage-detection in blade pitch control systems for wind turbines | |
CN105571784A (zh) | 一种发动机曲轴平衡率的检测方法 | |
CN109472055B (zh) | 基于最大似然法的密封动力特性系数识别方法 | |
CN206556890U (zh) | 一种转臂定位悬挂***综合性能试验台 | |
CN114486261B (zh) | 基于同步均方根之比的轴承外圈故障定位方法 | |
CN109388885A (zh) | 一种基于矩估计法的密封动力特性系数数值获取方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190524 Termination date: 20211013 |