CN108181105A - 基于逻辑回归和j散度的滚动轴承故障预诊方法及*** - Google Patents
基于逻辑回归和j散度的滚动轴承故障预诊方法及*** Download PDFInfo
- Publication number
- CN108181105A CN108181105A CN201711216940.4A CN201711216940A CN108181105A CN 108181105 A CN108181105 A CN 108181105A CN 201711216940 A CN201711216940 A CN 201711216940A CN 108181105 A CN108181105 A CN 108181105A
- Authority
- CN
- China
- Prior art keywords
- bearing
- divergences
- value
- sample
- divergence
- 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.)
- Granted
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于逻辑回归和J散度的滚动轴承故障预诊方法,包括以下步骤:采集轴承运行时不同故障位置的故障传感数据和正常状态轴承运行的正常传感数据,分别对故障传感数据和正常传感数据进行预处理和特征提取,建立故障位置特征样本和正常状态特征样本;通过建立好的故障位置的样本和正常状态的样本对逻辑回归模型进行训练,得到逻辑模型参数,建立逻辑回归模型;本发明依据轴承的健康衰退程度对轴承进行有针对性的故障诊断,同时根据不同工况不同对象,改变训练样本类型、特征值类型、健康阈值等参数,训练后的模型可以调整,具有实时性强、数据处理精度高、核心算法鲁棒性好、状态评估准确性高、诊断结果准确度高等优点。
Description
技术领域
本发明涉及机械产品质量可靠性评估与故障诊断技术领域,尤其涉及了一种基于逻辑回归和J散度的滚动轴承故障预诊方法及***。
背景技术
滚动轴承是旋转机械的重要支撑件,其实时运行状态直接决定了机械设备的可靠性,因此,对滚动轴承在线故障诊断研究非常重要,实时评估滚动轴承当前运行状态、定位故障部位和量化故障严重程度,对机械设备的维护和设计有极其重要的指导意义。本专利针对滚动轴承实现了一种基于逻辑回归和J散度的在线故障预诊办法,能够实时评估滚动轴承的运行状态,基于当前运行状态地定位故障部位。
现有在线故障诊断方法的研究有很多,目前成熟的方法有时域有效值峰值判断法、振幅概率密度分析法,冲击脉冲法等。这些方法只能给出故障或正常两种状态,事实上,轴承大部分时间的运行状态介于故障和正常之间,既非完全故障,亦非完全正常,可以将这种中间态称为健康衰退状态。轴承健康衰退状态不同于常见的轴承故障,是一个隐性的、渐变的、漫长的过程。处于此状态的轴承,在稳定衰退期,轴承的性能呈现稳定和缓慢的下降趋势并形成早期故障,从早期故障到最终故障状态,轴承性能会有较快速的下降,通过监测技术和信号处理技术将整个衰退过程量化和显性化,有助于更好地管理轴承的健康状态。显然,上述方法只给出故障或正常两种状态数据而丢失大部分健康衰退状态数据,无法准确评估轴承运行状态,有较高的误判率和漏判率。
发明内容
本发明针对现有技术中无法准确地评估轴承运行状态的缺点,提供了一种基于逻辑回归和J散度的滚动轴承故障预诊方法及***。
为了解决上述技术问题,本发明通过下述技术方案得以解决:
一种基于逻辑回归和J散度的滚动轴承故障预诊方法,包括以下步骤:
采集轴承运行时不同故障位置的故障传感数据和正常状态轴承运行的正常传感数据,分别对故障传感数据和正常传感数据进行预处理和特征提取,建立故障位置特征样本和正常状态特征样本;
通过建立好的故障位置的样本和正常状态的样本对逻辑回归模型进行训练,得到逻辑模型参数,建立逻辑回归模型;
采集待测轴承的实时传感数据,对所述实时传感数据进行预处理和特征提取,得到待测轴承的实时状态特征指标,将所述特征指标代入到建立好的所述逻辑回归模型中,经过计算得到待测轴承当前时刻的健康度;
将计算得到的待测轴承当前时刻的健康度和预先设定好的健康度阈值作比较,若待测轴承当前时刻的健康度低于设定的阈值,则通过J散度的故障诊断方法分别计算待测轴承的实时状态特征指标与故障位置特征样本的J散度、正常状态特征样本的J散度,根据故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,实现对待测轴承的故障诊断。
作为一种可实施方式,所述逻辑回归模型获取方法的具体步骤包括,
假设N个独立变量的向量,XT=(x1,x2,…,xN),yi∈{0,1},y为0或1,逻辑回归模型的数学表达式为:
yi=1表示事件发生;yi=0表示事件未发生;
pi(yi=1/xi)表示观测量相对于第i个事件发生的概率(取值在0到1之间);
表示回归截距;β1,β2,…βN表示回归系数;回归截距及回归系数通过极大似然估计法来求取。
作为一种可实施方式,所述轴承运行时不同故障位置的故障传感数据包括轴承的内圈故障、外圈故障和滚动体故障中的一种或者几种故障传感数据。
作为一种可实施方式,所述故障位置的特征样本包括内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本中的一种或者几种故障特征样本。
作为一种可实施方式,当故障位置的特征样本为内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本时,所述通过J散度的故障诊断方法分别计算待测轴承的实时状态特征指标与故障位置特征样本的J散度、正常状态特征样本的J散度,根据故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,实现对待测轴承的故障诊断,过程如下,
通过公式J散度公式:分别计算实时轴承状态的特征指标与内圈故障特征样本、外圈故障特征样本、滚动体故障特征样本、正常状态特征样本之间的J散度,得到四个J散度值,分别为:散度值J1,散度值J2,散度值J3,散度值J4,J散度公式中,S为正常状态信号的特征值;τ为未知状态信号的特征值;J(s,τ)为两者之间的J散度,N为信号特征值的个数,i为信号特征值的序列;
对比散度值J1、散度值J2、散度值J3和散度值J4,找出四个J散度值中的最小值,散度值越小,则轴承当前时刻状态与预设的预先设定好的健康度阈值对应的状态越相近,通过散度值的大小来确定轴承的故障类型。
一种基于逻辑回归和J散度的滚动轴承故障预诊***,包括:
采集预处理模块,用于采集轴承运行时不同故障位置的故障传感数据和正常状态轴承运行的正常传感数据,分别对故障传感数据和正常传感数据进行预处理和特征提取,建立故障位置特征样本和正常状态特征样本;
模型建立模块,用于通过建立好的故障位置的样本和正常状态的样本对逻辑回归模型进行训练,得到逻辑模型参数,建立逻辑回归模型;
采集计算模块,用于采集待测轴承的实时传感数据,对所述实时传感数据进行预处理和特征提取,得到待测轴承的实时状态特征指标,将所述特征指标代入到建立好的所述逻辑回归模型中,经过计算得到待测轴承当前时刻的健康度;
判断模块,用于将计算得到的待测轴承当前时刻的健康度和预先设定好的健康度阈值作比较,若待测轴承当前时刻的健康度低于设定的阈值,则通过J散度的故障诊断方法分别计算待测轴承的实时状态特征指标与故障位置特征样本的J散度、正常状态特征样本的J散度,根据故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,实现对待测轴承的故障诊断。
作为一种可实施方式,当故障位置的特征样本为内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本时,模型建立模块被设置为,
假设N个独立变量的向量,XT=(x1,x2,…,xN),yi∈{0,1},y为0或1,逻辑回归模型的数学表达式为:
yi=1表示事件发生;yi=0表示事件未发生;
pi(yi=1/xi)表示观测量相对于第i个事件发生的概率(取值在0到1之间);
表示回归截距;β1,β2,…βN表示回归系数;回归截距及回归系数通过极大似然估计法来求取。
作为一种可实施方式,所述采集预处理模块被设置为,所述轴承运行时不同故障位置的故障传感数据包括轴承的内圈故障、外圈故障和滚动体故障中的一种或者几种故障传感数据。
作为一种可实施方式,所述采集预处理模块被设置为,所述故障位置的特征样本包括内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本中的一种或者几种故障特征样本。
作为一种可实施方式,所述判断模块被设置为,
通过公式J散度公式:分别计算实时轴承状态的特征指标与内圈故障特征样本、外圈故障特征样本、滚动体故障特征样本、正常状态特征样本之间的J散度,得到四个J散度值,分别为:散度值J1,散度值J2,散度值J3,散度值J4,J散度公式中,S为正常状态信号的特征值;τ为未知状态信号的特征值;J(s,τ)为两者之间的J散度,N为信号特征值的个数,i为信号特征值的序列;
对比散度值J1、散度值J2、散度值J3和散度值J4,找出四个J散度值中的最小值,散度值越小,则轴承当前时刻状态与预设的预先设定好的健康度阈值对应的状态越相近,通过散度值的大小来确定轴承的故障类型。
本发明由于采用了以上技术方案,具有显著的技术效果:
本发明依据轴承的健康衰退程度对轴承进行有针对性的故障诊断,同时根据不同工况不同对象,改变训练样本类型、特征值类型、健康阈值等参数,训练后的模型可以调整,具有实时性强、数据处理精度高、核心算法鲁棒性好、状态评估准确性高、诊断结果准确度高等优点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明的示例性方法的流程示意图;
图2是本发明的具体流程示意图;
图3是本发明的实验装置示意图;
图4是本发明中绘制的全寿命周期的评估指数示意图;
图5是本发明的***结构示意图。
具体实施方式
下面结合实施例对本发明做进一步的详细说明,以下实施例是对本发明的解释而本发明并不局限于以下实施例。
示例性方法:
一种基于逻辑回归和J散度的滚动轴承故障预诊方法,如图1所示,包括以下步骤:
S1:采集轴承运行时不同故障位置的故障传感数据和正常状态轴承运行的正常传感数据,分别对故障传感数据和正常传感数据进行预处理和特征提取,建立故障位置特征样本和正常状态特征样本;
S2:通过建立好的故障位置的样本和正常状态的样本对逻辑回归模型进行训练,得到逻辑模型参数,建立逻辑回归模型;
S3:采集待测轴承的实时传感数据,对所述实时传感数据进行预处理和特征提取,得到待测轴承的实时状态特征指标,将所述特征指标代入到建立好的所述逻辑回归模型中,经过计算得到待测轴承当前时刻的健康度;
S4:将计算得到的待测轴承当前时刻的健康度和预先设定好的健康度阈值作比较,若待测轴承当前时刻的健康度低于设定的阈值,则通过J散度的故障诊断方法分别计算待测轴承的实时状态特征指标与故障位置特征样本的J散度、正常状态特征样本的J散度,根据故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,实现对待测轴承的故障诊断。
为了能更好的说明通过J散度的故障诊断方法数据,所述轴承运行时不同故障位置的故障传感数据包括轴承的内圈故障、外圈故障和滚动体故障中的一种或者几种故障传感数据,通过获得的不同故障位置的故障传感数据经过预处理之后得到故障位置的特征样本,与一种或者几种故障传感数据对应的故障位置的特征样本包括内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本,也就是说,如果选取一种或者几种故障传感数据,那么,就选择一种或者几种对应的故障位置的特征样本,在此,可以将其解释为:为适应不同类型轴承、不同工况,可选择提取不同的特征类型,以适应并满足不同的分析需求,同时针对实际应用的情况,可以对故障位置进行自适应调整,例如在某些情况下保持架的故障发生频率较高,可以增加保持架故障的数据样本,进行特征提取;或者在某些工况下,外圈故障几乎不会发生,可以删除外圈故障的数据样本。
如果是轴承运行时不同故障位置的故障传感数据同时包括轴承的内圈故障、外圈故障和滚动体故障这三种故障传感数据时,所述通过J散度的故障诊断方法分别计算待测轴承的实时状态特征指标与故障位置特征样本的J散度、正常状态特征样本的J散度,根据故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,实现对待测轴承的故障诊断,过程如下,
通过公式J散度公式:分别计算实时轴承状态的特征指标与内圈故障特征样本、外圈故障特征样本、滚动体故障特征样本、正常状态特征样本之间的J散度,得到四个J散度值,分别为:散度值J1,散度值J2,散度值J3,散度值J4,J散度公式中,S为正常状态信号的特征值;τ为未知状态信号的特征值;J(s,τ)为两者之间的J散度,N为信号特征值的个数,i为信号特征值的序列;
对比散度值J1、散度值J2、散度值J3和散度值J4,找出四个J散度值中的最小值,散度值越小,则轴承当前时刻状态与预设的预先设定好的健康度阈值对应的状态越相近,通过散度值的大小来确定轴承的故障类型。
实施例2:
结合具体数据和附图来解释本发明的方法,在本实施例中,是以三种故障传感数据为例来说明的,具体装置图如附图2所示,一种基于逻辑回归和J散度的滚动轴承故障预诊方法,包括以下步骤:
步骤1:采集轴承运行时不同故障位置的三种故障传感数据和正常状态轴承运行的正常传感数据,分别对三种故障传感数据和正常传感数据进行预处理和特征提取,建立三种故障位置特征样本和正常状态特征样本;
步骤2:通过建立好的三种故障位置的样本和正常状态的样本对逻辑回归模型进行训练,得到逻辑模型参数,建立逻辑回归模型;
步骤3:采集待测轴承的实时传感数据,对所述实时传感数据进行预处理和特征提取,得到待测轴承的实时状态特征指标,将所述特征指标代入到建立好的所述逻辑回归模型中,经过计算得到待测轴承当前时刻的健康度;
步骤4:将计算得到的待测轴承当前时刻的健康度和预先设定好的健康度阈值作比较,若待测轴承当前时刻的健康度低于设定的阈值,则通过J散度的故障诊断方法分别计算待测轴承的实时状态特征指标与三种故障位置特征样本的J散度、正常状态特征样本的J散度,根据三种故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,实现对待测轴承的故障诊断。
该方法只有在轴承当前时刻的健康度低于设定的阈值,才通过J散度的故障诊断方法分别计算当前时刻的轴承状态与三种故障位置特征样本的J散度、正常状态特征样本的J散度,根据故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,也就是说,如果轴承当前时刻的健康度不低于设定的阈值时,则说明轴承不存在故障,不需要进行故障诊断。
在本实施例中,三种故障传感数据为轴承的内圈故障、外圈故障和滚动体故障的故障传感数据,通过获得的不同故障位置的故障传感数据经过预处理之后得到故障位置的特征样本,与三种对应的故障位置的特征样本包括内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本。
在本实施例中,逻辑回归模型的建立依据以下理论:假设N个独立变量的向量,XT=(x1,x2,…,xN),yi∈{0,1},y为0或1,逻辑回归模型的数学表达式为:
yi=1表示事件发生;yi=0表示事件未发生;
pi(yi=1/xi)表示观测量相对于第i个事件发生的概率(取值在0到1之间);
表示回归截距;β1,β2,…βN表示回归系数;回归截距及回归系数通过极大似然估计法来求取。
当故障位置的特征样本为内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本时,所述通过J散度的故障诊断方法分别计算待测轴承的实时状态特征指标与故障位置特征样本的J散度、正常状态特征样本的J散度,根据故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,实现对待测轴承的故障诊断,过程如下,
通过公式J散度公式:分别计算实时轴承状态的特征指标与内圈故障特征样本、外圈故障特征样本、滚动体故障特征样本、正常状态特征样本之间的J散度,得到四个J散度值,分别为:散度值J1,散度值J2,散度值J3,散度值J4,J散度公式中,S为正常状态信号的特征值;τ为未知状态信号的特征值;J(s,τ)为两者之间的J散度,N为信号特征值的个数,i为信号特征值的序列;
对比散度值J1、散度值J2、散度值J3和散度值J4,找出四个J散度值中的最小值,散度值越小,则轴承当前时刻状态与预设的预先设定好的健康度阈值对应的状态越相近,通过散度值的大小来确定轴承的故障类型。
结合实际情况来讲,在本发明中以Rexnord ZA-2115轴承为例,将4个该型号的轴承安装在轴上,交流电机转速为2000rpm,施加27kN径向载荷,利用PCB353B33振动传感器实时采集轴承振动信号,实验装置如附图3所示。
计算轴承的健康度,直至其中的某一轴承完全失效,绘制失效轴承的全生命健康评估曲线,其健康度变化如附图4所示:将健康度阈值设置为0.3(预设的健康度阈值),从附图3中可以看出,轴承状态评估指标CV值在0-1之间变化,在此,CV值的物理意义为轴承当前健康状态的健康度。CV值为1则说明滚动轴承完全属于正常性能状态,CV值为0则说明滚动轴承处于性能最差状态,为适应不同的工况、轴承型号,上述两种状态及其CV值是可以做自定义设置的。例如,若应用工程对轴承性能要求比较严格,可以选取健康衰退的中期阶段的样本的CV值设置为0,并采用逻辑回归模型进行训练,得到严格性能要求下的健康度曲线,CV值在第一天到第二十五天基本处于[0.7,1]的范围内,数值变化不大,定义为轴承健康期。随着轴承继续工作,当轴承工作至第25天时,轴承的CV值开始逐渐变小,并且低于预先设置的健康期的健康度阈值,处于[0.3,0.7]的范围内,定义为健康衰退期。当轴承工作至第27天时,轴承的健康度急剧下降,CV值小于0.3,轴承的健康度已达到轴承故障阈值,则开始进行故障诊断并对故障类型进行判断。在本示例中,启动故障诊断的健康度阈值被设定成0.3,根据轴承实际使用情况,该值可以进行调整,以适应不同的分析需求,在高可靠性要求的应用场景下,可以将阈值提高;在追求高经济性,且意外停机损失不大的应用场景下,可以将阈值降低。
当轴承的健康度已达到轴承故障阈值,则开始进行故障诊断并对故障类型进行判断,此时,计算实时轴承状态的特征指标与正常状态,内圈故障、外圈故障、滚动体故障的特征指标之间的J散度,共得到四个J散度值:散度值J1,散度值J2,散度值J3,散度值J4,如表1所示。
表1实时状态与标准故障样本之间的J散度
由表1可知,通过J散度可以清楚的区分轴承的状态,实时轴承状态与三种标准状态之间的J散度值差异明显,其中实时轴承状态与内圈故障状态的J散度值J2越小,代表轴承当前状态与内圈故障状态最相近。经过多次实验计算验证,本方法可以清楚的区分轴承不同的状态,有效的对轴承故障进行诊断,准确率和灵敏度较高。
实施例2:
一种基于逻辑回归和J散度的滚动轴承故障预诊***,如图5所示,包括:
采集预处理模块1,用于采集轴承运行时不同故障位置的故障传感数据和正常状态轴承运行的正常传感数据,分别对故障传感数据和正常传感数据进行预处理和特征提取,建立故障位置特征样本和正常状态特征样本;
模型建立模块2,用于通过建立好的故障位置的样本和正常状态的样本对逻辑回归模型进行训练,得到逻辑模型参数,建立逻辑回归模型;
采集计算模块3,用于采集待测轴承的实时传感数据,对所述实时传感数据进行预处理和特征提取,得到待测轴承的实时状态特征指标,将所述特征指标代入到建立好的所述逻辑回归模型中,经过计算得到待测轴承当前时刻的健康度;
判断模块4,用于将计算得到的待测轴承当前时刻的健康度和预先设定好的健康度阈值作比较,若待测轴承当前时刻的健康度低于设定的阈值,则通过J散度的故障诊断方法分别计算待测轴承的实时状态特征指标与故障位置特征样本的J散度、正常状态特征样本的J散度,根据故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,实现对待测轴承的故障诊断。
具体地,模型建立模块2被设置为,
假设N个独立变量的向量,XT=(x1,x2,…,xN),yi∈{0,1},y为0或1,逻辑回归模型的数学表达式为:
yi=1表示事件发生;yi=0表示事件未发生;
pi(yi=1/xi)表示观测量相对于第i个事件发生的概率(取值在0到1之间);
表示回归截距;β1,β2,…βN表示回归系数;回归截距及回归系数通过极大似然估计法来求取。
更进一步地,所述采集预处理模块1被设置为,所述轴承运行时不同故障位置的故障传感数据包括轴承的内圈故障、外圈故障和滚动体故障中的一种或者几种故障传感数据。所述采集预处理模块被设置为,所述故障位置的特征样本包括内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本中的一种或者几种。
更进一步来讲,当故障位置的特征样本为内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本时,所述判断模块被设置为,
通过公式J散度公式:分别计算实时轴承状态的特征指标与内圈故障特征样本、外圈故障特征样本、滚动体故障特征样本、正常状态特征样本之间的J散度,得到四个J散度值,分别为:散度值J1,散度值J2,散度值J3,散度值J4,J散度公式中,S为正常状态信号的特征值;τ为未知状态信号的特征值;J(s,τ)为两者之间的J散度,N为信号特征值的个数,i为信号特征值的序列;
对比散度值J1、散度值J2、散度值J3和散度值J4,找出四个J散度值中的最小值,散度值越小,则轴承当前时刻状态与预设的预先设定好的健康度阈值对应的状态越相近,通过散度值的大小来确定轴承的故障类型。
此外,需要说明的是,本说明书中所描述的具体实施例,其零、部件的形状、所取名称等可以不同。凡依本发明专利构思所述的构造、特征及原理所做的等效或简单变化,均包括于本发明专利的保护范围内。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,只要不偏离本发明的结构或者超越本权利要求书所定义的范围,均应属于本发明的保护范围。
Claims (10)
1.一种基于逻辑回归和J散度的滚动轴承故障预诊方法,其特征在于包括以下步骤:
采集轴承运行时不同故障位置的故障传感数据和正常状态轴承运行的正常传感数据,分别对故障传感数据和正常传感数据进行预处理和特征提取,建立故障位置特征样本和正常状态特征样本;
通过建立好的故障位置的样本和正常状态的样本对逻辑回归模型进行训练,得到逻辑模型参数,建立逻辑回归模型;
采集待测轴承的实时传感数据,对所述实时传感数据进行预处理和特征提取,得到待测轴承的实时状态特征指标,将所述特征指标代入到建立好的所述逻辑回归模型中,经过计算得到待测轴承当前时刻的健康度;
将计算得到的待测轴承当前时刻的健康度和预先设定好的健康度阈值作比较,若待测轴承当前时刻的健康度低于设定的阈值,则通过J散度的故障诊断方法分别计算待测轴承的实时状态特征指标与故障位置特征样本的J散度、正常状态特征样本的J散度,根据故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,实现对待测轴承的故障诊断。
2.根据权利要求1所述的基于逻辑回归和J散度的滚动轴承故障预诊方法,其特征在于:所述逻辑回归模型获取方法的具体步骤包括,
假设N个独立变量的向量,XT=(x1,x2,…,xN),yi∈{0,1},y为0或1,逻辑回归模型的数学表达式为:
yi=1表示事件发生;yi=0表示事件未发生;
pi(yi=1/xi)表示观测量相对于第i个事件发生的概率(取值在0到1之间);
表示回归截距;β1,β2,…βN表示回归系数;回归截距及回归系数通过极大似然估计法来求取。
3.根据权利要求1所述的基于逻辑回归和J散度的滚动轴承故障预诊方法,其特征在于:所述轴承运行时不同故障位置的故障传感数据包括轴承的内圈故障、外圈故障和滚动体故障中的一种或者几种故障传感数据。
4.根据权利要求3所述的基于逻辑回归和J散度的滚动轴承故障预诊方法,其特征在于:所述故障位置的特征样本包括内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本中的一种或者几种故障特征样本。
5.根据权利要求4所述的基于逻辑回归和J散度的滚动轴承故障预诊方法,其特征在于:当故障位置的特征样本为内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本时,所述通过J散度的故障诊断方法分别计算待测轴承的实时状态特征指标与故障位置特征样本的J散度、正常状态特征样本的J散度,根据故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,实现对待测轴承的故障诊断,过程如下,
通过公式J散度公式:分别计算实时轴承状态的特征指标与内圈故障特征样本、外圈故障特征样本、滚动体故障特征样本、正常状态特征样本之间的J散度,得到四个J散度值,分别为:散度值J1,散度值J2,散度值J3,散度值J4,J散度公式中,S为正常状态信号的特征值;τ为未知状态信号的特征值;J(s,τ)为两者之间的J散度,N为信号特征值的个数,i为信号特征值的序列;
对比散度值J1、散度值J2、散度值J3和散度值J4,找出四个J散度值中的最小值,散度值越小,则轴承当前时刻状态与预设的预先设定好的健康度阈值对应的状态越相近,通过散度值的大小来确定轴承的故障类型。
6.一种基于逻辑回归和J散度的滚动轴承故障预诊***,其特征在于,包括:
采集预处理模块,用于采集轴承运行时不同故障位置的故障传感数据和正常状态轴承运行的正常传感数据,分别对故障传感数据和正常传感数据进行预处理和特征提取,建立故障位置特征样本和正常状态特征样本;
模型建立模块,用于通过建立好的故障位置的样本和正常状态的样本对逻辑回归模型进行训练,得到逻辑模型参数,建立逻辑回归模型;
采集计算模块,用于采集待测轴承的实时传感数据,对所述实时传感数据进行预处理和特征提取,得到待测轴承的实时状态特征指标,将所述特征指标代入到建立好的所述逻辑回归模型中,经过计算得到待测轴承当前时刻的健康度;
判断模块,用于将计算得到的待测轴承当前时刻的健康度和预先设定好的健康度阈值作比较,若待测轴承当前时刻的健康度低于设定的阈值,则通过J散度的故障诊断方法分别计算待测轴承的实时状态特征指标与故障位置特征样本的J散度、正常状态特征样本的J散度,根据故障位置的样本的J散度和正常状态的样本的J散度的大小判断轴承的故障位置,实现对待测轴承的故障诊断。
7.根据权利要求6所述的基于逻辑回归和J散度的滚动轴承故障预诊***,其特征在于:模型建立模块被设置为,
假设N个独立变量的向量,XT=(x1,x2,…,xN),yi∈{0,1},y为0或1,逻辑回归模型的数学表达式为:
yi=1表示事件发生;yi=0表示事件未发生;
pi(yi=1/xi)表示观测量相对于第i个事件发生的概率(取值在0到1之间);
表示回归截距;β1,β2,…βN表示回归系数;回归截距及回归系数通过极大似然估计法来求取。
8.根据权利要求6所述的基于逻辑回归和J散度的滚动轴承故障预诊***,其特征在于:所述采集预处理模块被设置为,所述轴承运行时不同故障位置的故障传感数据包括轴承的内圈故障、外圈故障和滚动体故障中的一种或者几种故障传感数据。
9.根据权利要求8所述的基于逻辑回归和J散度的滚动轴承故障预诊***,其特征在于:所述采集预处理模块被设置为,所述故障位置的特征样本包括内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本中的一种或者几种故障特征样本。
10.根据权利要求9所述的基于逻辑回归和J散度的滚动轴承故障预诊***,其特征在于:当故障位置的特征样本为内圈故障特征样本、外圈故障特征样本和滚动体故障特征样本时,所述判断模块被设置为,
通过公式J散度公式:分别计算实时轴承状态的特征指标与内圈故障特征样本、外圈故障特征样本、滚动体故障特征样本、正常状态特征样本之间的J散度,得到四个J散度值,分别为:散度值J1,散度值J2,散度值J3,散度值J4,J散度公式中,S为正常状态信号的特征值;τ为未知状态信号的特征值;J(s,τ)为两者之间的J散度,N为信号特征值的个数,i为信号特征值的序列;
对比散度值J1、散度值J2、散度值J3和散度值J4,找出四个J散度值中的最小值,散度值越小,则轴承当前时刻状态与预设的预先设定好的健康度阈值对应的状态越相近,通过散度值的大小来确定轴承的故障类型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711216940.4A CN108181105B (zh) | 2017-11-28 | 2017-11-28 | 基于逻辑回归和j散度的滚动轴承故障预诊方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711216940.4A CN108181105B (zh) | 2017-11-28 | 2017-11-28 | 基于逻辑回归和j散度的滚动轴承故障预诊方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108181105A true CN108181105A (zh) | 2018-06-19 |
CN108181105B CN108181105B (zh) | 2020-08-07 |
Family
ID=62545456
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711216940.4A Active CN108181105B (zh) | 2017-11-28 | 2017-11-28 | 基于逻辑回归和j散度的滚动轴承故障预诊方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108181105B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110987433A (zh) * | 2019-12-13 | 2020-04-10 | 华中科技大学 | 一种基于高频信号特征幅值的轴承故障预警方法 |
CN111076962A (zh) * | 2020-01-03 | 2020-04-28 | 中国电建集团华东勘测设计研究院有限公司 | 一种面向智慧水电厂的机电设备故障诊断方法 |
CN111260822A (zh) * | 2019-12-31 | 2020-06-09 | 杭州安脉盛智能技术有限公司 | 一种基于大数据的轨道交通车辆健康状态分析方法及终端 |
CN111651933A (zh) * | 2020-05-22 | 2020-09-11 | 宁波诺丁汉新材料研究院有限公司 | 基于统计推断的工业锅炉故障预警方法及*** |
CN112396344A (zh) * | 2020-11-30 | 2021-02-23 | 天津大学 | 一种基于产品质量的化工过程可靠性在线评估方法 |
WO2021042748A1 (zh) * | 2019-09-02 | 2021-03-11 | 佛山科学技术学院 | 一种基于概率分布的轴承工况稳定性评估方法及装置 |
CN112614342A (zh) * | 2020-12-10 | 2021-04-06 | 大唐高鸿数据网络技术股份有限公司 | 一种道路异常事件的预警方法、车载设备及路侧设备 |
CN112948223A (zh) * | 2019-11-26 | 2021-06-11 | 北京沃东天骏信息技术有限公司 | 一种监测运行情况的方法和装置 |
US20220351554A1 (en) * | 2021-04-28 | 2022-11-03 | Dana Automotive Systems Group, Llc | Systems and methods for prediction of component degradation |
JP7521461B2 (ja) | 2021-03-09 | 2024-07-24 | 株式会社明電舎 | 設備の異常診断装置及び異常診断方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102682180A (zh) * | 2011-03-07 | 2012-09-19 | 上海宝信软件股份有限公司 | 旋转式机械设备性能退化的评估方法 |
CN103679378A (zh) * | 2013-12-20 | 2014-03-26 | 北京航天测控技术有限公司 | 基于遥测数据评估航天器健康状态的方法及装置 |
CN103868692A (zh) * | 2014-03-18 | 2014-06-18 | 电子科技大学 | 基于核密度估计和k-l散度的旋转机械故障诊断方法 |
CN104392082A (zh) * | 2014-07-10 | 2015-03-04 | 中山火炬职业技术学院 | 一种基于振动监测的风力发电机组齿轮箱早期故障诊断方法 |
CN105241680A (zh) * | 2015-08-26 | 2016-01-13 | 电子科技大学 | 一种基于概率密度函数的旋转机械健康状态评估方法 |
CN105510023A (zh) * | 2015-11-24 | 2016-04-20 | 国网内蒙古东部电力有限公司电力科学研究院 | 基于散度指标的变工况风电行星齿轮箱故障诊断方法 |
CN106053066A (zh) * | 2016-05-23 | 2016-10-26 | 华东交通大学 | 基于经验模态分解和逻辑回归的滚动轴承性能退化评估方法 |
-
2017
- 2017-11-28 CN CN201711216940.4A patent/CN108181105B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102682180A (zh) * | 2011-03-07 | 2012-09-19 | 上海宝信软件股份有限公司 | 旋转式机械设备性能退化的评估方法 |
CN103679378A (zh) * | 2013-12-20 | 2014-03-26 | 北京航天测控技术有限公司 | 基于遥测数据评估航天器健康状态的方法及装置 |
CN103868692A (zh) * | 2014-03-18 | 2014-06-18 | 电子科技大学 | 基于核密度估计和k-l散度的旋转机械故障诊断方法 |
CN104392082A (zh) * | 2014-07-10 | 2015-03-04 | 中山火炬职业技术学院 | 一种基于振动监测的风力发电机组齿轮箱早期故障诊断方法 |
CN105241680A (zh) * | 2015-08-26 | 2016-01-13 | 电子科技大学 | 一种基于概率密度函数的旋转机械健康状态评估方法 |
CN105510023A (zh) * | 2015-11-24 | 2016-04-20 | 国网内蒙古东部电力有限公司电力科学研究院 | 基于散度指标的变工况风电行星齿轮箱故障诊断方法 |
CN106053066A (zh) * | 2016-05-23 | 2016-10-26 | 华东交通大学 | 基于经验模态分解和逻辑回归的滚动轴承性能退化评估方法 |
Non-Patent Citations (2)
Title |
---|
郭艳平 等: ""基于经验模态分解和散度指标的风力发电机滚动轴承故障诊断方法"", 《电力***保护与控制》 * |
钟鑫 等: ""基于逻辑回归的滚动轴承性能退化评估"", 《科技信息》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021042748A1 (zh) * | 2019-09-02 | 2021-03-11 | 佛山科学技术学院 | 一种基于概率分布的轴承工况稳定性评估方法及装置 |
CN112948223A (zh) * | 2019-11-26 | 2021-06-11 | 北京沃东天骏信息技术有限公司 | 一种监测运行情况的方法和装置 |
CN110987433A (zh) * | 2019-12-13 | 2020-04-10 | 华中科技大学 | 一种基于高频信号特征幅值的轴承故障预警方法 |
CN110987433B (zh) * | 2019-12-13 | 2021-05-18 | 华中科技大学 | 一种基于高频信号特征幅值的轴承故障预警方法 |
CN111260822A (zh) * | 2019-12-31 | 2020-06-09 | 杭州安脉盛智能技术有限公司 | 一种基于大数据的轨道交通车辆健康状态分析方法及终端 |
CN111260822B (zh) * | 2019-12-31 | 2022-07-26 | 杭州安脉盛智能技术有限公司 | 一种基于大数据的轨道交通车辆健康状态分析方法及终端 |
CN111076962A (zh) * | 2020-01-03 | 2020-04-28 | 中国电建集团华东勘测设计研究院有限公司 | 一种面向智慧水电厂的机电设备故障诊断方法 |
CN111076962B (zh) * | 2020-01-03 | 2022-03-29 | 中国电建集团华东勘测设计研究院有限公司 | 一种面向智慧水电厂的机电设备故障诊断方法 |
CN111651933A (zh) * | 2020-05-22 | 2020-09-11 | 宁波诺丁汉新材料研究院有限公司 | 基于统计推断的工业锅炉故障预警方法及*** |
CN111651933B (zh) * | 2020-05-22 | 2023-09-26 | 宁波诺丁汉新材料研究院有限公司 | 基于统计推断的工业锅炉故障预警方法及*** |
CN112396344A (zh) * | 2020-11-30 | 2021-02-23 | 天津大学 | 一种基于产品质量的化工过程可靠性在线评估方法 |
CN112614342A (zh) * | 2020-12-10 | 2021-04-06 | 大唐高鸿数据网络技术股份有限公司 | 一种道路异常事件的预警方法、车载设备及路侧设备 |
JP7521461B2 (ja) | 2021-03-09 | 2024-07-24 | 株式会社明電舎 | 設備の異常診断装置及び異常診断方法 |
US20220351554A1 (en) * | 2021-04-28 | 2022-11-03 | Dana Automotive Systems Group, Llc | Systems and methods for prediction of component degradation |
Also Published As
Publication number | Publication date |
---|---|
CN108181105B (zh) | 2020-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108181105A (zh) | 基于逻辑回归和j散度的滚动轴承故障预诊方法及*** | |
CN110823576B (zh) | 基于生成对抗网络的机械异常检测方法 | |
CN108267312B (zh) | 一种基于快速搜索算法的地铁列车轴承智能诊断方法 | |
CN110687896A (zh) | 故障诊断方法、装置、设备及可读介质 | |
CN106429689B (zh) | 基于物联网大数据支撑的电梯维保*** | |
CN109213121A (zh) | 一种风机制动***卡缸故障诊断方法 | |
CN109858104A (zh) | 一种滚动轴承健康评估与故障诊断方法及监测*** | |
CN111665066B (zh) | 基于卷积神经网络的设备故障自适应上下预警界生成方法 | |
CN108896299A (zh) | 一种齿轮箱故障检测方法 | |
CN110060368A (zh) | 基于潜在特征编码的机械异常检测方法 | |
CN111611294B (zh) | 星敏感器数据异常检测方法 | |
CN112729834A (zh) | 一种轴承故障的诊断方法、装置和*** | |
CN109855873A (zh) | 汽轮机大轴的故障诊断方法及装置 | |
CN114414245A (zh) | 轴承故障检测方法与计算机可读存储介质 | |
CN109894476B (zh) | 一种冷轧硅钢生产线液压设备故障诊断方法和装置 | |
CN109990803A (zh) | 检测***异常的方法、装置及传感器处理的方法、装置 | |
CN116611184A (zh) | 一种齿轮箱的故障检测方法、装置及介质 | |
CN113487141B (zh) | 一种多源信息聚类融合的轴承状态评估方法 | |
CN114611247A (zh) | 基于时频特征序列的轴承剩余寿命预测方法及装置 | |
Kolar et al. | Identification of inability states of rotating machinery subsystems using industrial IoT and convolutional neural network–initial research | |
Zhang et al. | A health indicator construction method based on the information entropy of specific frequency band and GGRBM | |
CN117851954B (zh) | 基于数据分析的轴承加工设备运行质量检测***和方法 | |
CN111368379A (zh) | 一种改进hwpso-wphm模型的滚动轴承可靠度评估方法和*** | |
Lee et al. | Fuel Consumption Prediction and Life Cycle History Management System Using Historical Data of Agricultural Machinery | |
Yamashita et al. | Defect Detection in Multiple Product Variants Using Hammering Test with Machine Learning |
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 |