CN108871761A - 一种齿轮早期故障特征提取方法 - Google Patents

Abstract

本发明公开了一种齿轮早期故障特征提取方法，包括以下步骤：包括以下步骤：（1）根据运行可靠性理论，计算齿轮的动态故障概率；（2）基于齿轮振动信号方差和齿轮振动信号熵的脆弱性指标，衡量齿轮箱中齿轮的状态脆弱性；（3）在齿轮排序方面，采取分类排序的方式：为了平衡计算时间和计算精度的要求，仅对涉及正常和停止后的大故障概率齿轮(如重过载和恶劣环境下的齿轮)，搜索由这些齿轮停止运行造成的故障率，对故障率进行排序，然后对于较小故障概率的齿轮，使用状态和结构脆弱性指标加权得到的综合脆弱性指标进行快速排序，实现齿轮早期故障特征提取。

Description

一种齿轮早期故障特征提取方法

技术领域

本发明涉及齿轮故障诊断技术领域，特别是一种齿轮早期故障特征提取方法。

背景技术

齿轮是机械装备普遍采用的一种机械传动形式，由于工作过程中受脉动循环应力作用，齿轮极易受到损害而出现故障，造成机械传动系统失效，进而影响整个机械装备的安全可靠运行。齿轮早期故障振动信号特征微弱，且随着齿轮箱转速和负载的变化而表现出非稳态特征，影响了故障诊断的可靠性和精确性。

发明内容

本发明技术提出一种齿轮早期的故障特征提取方法，能够有效提取齿轮传动非线性早期微弱故障特征，精确实现齿轮传动早期故障诊断。

本发明可以通过以下技术方案来实现：

本发明公开了一种齿轮早期故障特征提取方法，包括以下步骤：

(1)根据运行可靠性理论，考虑齿轮所处齿轮箱、环境噪声和自身振动信号，采用计入峭度统计因素和负载率的齿轮运行可靠性模型，由此计算齿轮的动态故障概率；

(2)基于齿轮振动信号方差和齿轮振动信号熵的脆弱性指标，衡量齿轮箱中齿轮的状态脆弱性；

(3)在齿轮排序方面，采取分类排序的方式：为了平衡计算时间和计算精度的要求，仅对涉及正常和停止后的大故障概率齿轮(如重过载和恶劣环境下的齿轮)，搜索由这些齿轮停止运行造成的故障率，对故障率进行排序，然后对于较小故障概率的齿轮，使用状态和结构脆弱性指标加权得到的综合脆弱性指标进行快速排序，实现齿轮早期故障特征提取。

进一步的，在第(1)步中，综合考虑齿轮所处齿轮箱、环境噪声和自身振动信号后，齿轮的可靠性模型下式所示：

进一步的，在第(2)步中，所述脆弱性指标包括基于负载率方差的齿轮脆弱性指标、基于负载率熵的齿轮脆弱性指标和齿轮的磨损系数指标。

进一步的，所述基于负载率方差的齿轮脆弱性指标中，齿轮箱正常运行各齿轮负载率方差可表示为：

式中，n和η_i0分别为正常运行时齿轮箱齿轮总数和齿轮i的负载率，

齿轮k退出运行后，齿轮箱的振动信号方差可表示为：

式中，η_k,i为齿轮k停止运行后，齿轮i的负载率，

综上所述，定义齿轮k关于负载率方差的脆弱性指标如下式所示：

。

进一步的，所述基于负载率熵的齿轮脆弱性指标中，定义齿轮箱正常运行时各齿轮的负载率熵H₀如下式所示：

式中，η_i0(或η_j0)分别为正常运行时齿轮箱齿轮总数和齿轮i(或j)的负载率，

齿轮k停止运行后，齿轮箱齿轮的负载率熵H_k定义如下式所示：

式中，η_k,i(或η_k,j)为齿轮k退出运行后，齿轮i(或j)的负载率；

定义齿轮k基于负载率熵的脆弱性指标如下式所示：

同理，可对下式的指标进行归一化，可得基于负载率熵的归一化指标如下式所示。

式中，

进一步的，所述齿轮的磨损系数指标中，齿轮k的磨损系数B_k(m,n)定义为：

式中，I^(ij)(m,n)为在“电机—负荷”节点对(i,j)间加入单位驱动后，在齿轮k上引起的压力；W_i为电机节点权重，可以取相应电机组的实时输出应力；W_j为负荷权重，取相应负荷的实时负载；G和L分别为电机和负荷集合；

对磨损系数进行归一化，得到归一化的磨损系数指标为：

。

本发明齿轮早期故障特征提取方法具有如下有益的技术效果：

采用计及峭度统计因素和齿轮振动信号的齿轮运行可靠性模型，可计算下动态变化的齿轮故障概率，能够有效提取齿轮传动非线性早期微弱故障特征，精确实现齿轮传动早期故障诊断。

具体实施方式

为了使本技术领域的人员更好地理解本发明的技术方案，下面结合实施例及对本发明产品作进一步详细的说明。

本发明公开了一种齿轮早期故障特征提取方法，包括以下步骤：

(1)根据运行可靠性理论，考虑齿轮所处齿轮箱、环境噪声和自身振动信号，采用计入峭度统计因素和负载率的齿轮运行可靠性模型，由此计算齿轮的动态故障概率；

(2)基于齿轮振动信号方差和齿轮振动信号熵的脆弱性指标，衡量齿轮箱中齿轮的状态脆弱性；

(3)在齿轮排序方面，采取分类排序的方式：为了平衡计算时间和计算精度的要求，仅对涉及正常和停止后的大故障概率齿轮(如重过载和恶劣环境下的齿轮)，搜索由这些齿轮停止运行造成的故障率，对故障率进行排序，然后对于较小故障概率的齿轮，使用状态和结构脆弱性指标加权得到的综合脆弱性指标进行快速排序，实现齿轮早期故障特征提取。

进一步的，在第(1)步中，综合考虑齿轮所处齿轮箱、环境噪声和自身振动信号后，齿轮的可靠性模型下式所示：

进一步的，在第(2)步中，所述脆弱性指标包括基于负载率方差的齿轮脆弱性指标、基于负载率熵的齿轮脆弱性指标和齿轮的磨损系数指标。

进一步的，所述基于负载率方差的齿轮脆弱性指标中，齿轮箱正常运行各齿轮负载率方差可表示为：

式中，n和η_i0分别为正常运行时齿轮箱齿轮总数和齿轮i的负载率，

齿轮k退出运行后，齿轮箱的振动信号方差可表示为：

式中，η_k,i为齿轮k停止运行后，齿轮i的负载率，

综上所述，定义齿轮k关于负载率方差的脆弱性指标如下式所示：

进一步的，所述基于负载率熵的齿轮脆弱性指标中，定义齿轮箱正常运行时各齿轮的负载率熵H₀如下式所示：

式中，η_i0(或η_j0)分别为正常运行时齿轮箱齿轮总数和齿轮i(或j)的负载率，

齿轮k停止运行后，齿轮箱齿轮的负载率熵H_k定义如下式所示：

式中，η_k,i(或η_k,j)为齿轮k退出运行后，齿轮i(或j)的负载率；

定义齿轮k基于负载率熵的脆弱性指标如下式所示：

同理，可对下式的指标进行归一化，可得基于负载率熵的归一化指标如下式所示。

式中，

进一步的，所述齿轮的磨损系数指标中，齿轮k的磨损系数B_k(m,n)定义为：

式中，I^(ij)(m,n)为在“电机—负荷”节点对(i,j)间加入单位驱动后，在齿轮k上引起的压力；W_i为电机节点权重，可以取相应电机组的实时输出应力；W_j为负荷权重，取相应负荷的实时负载；G和L分别为电机和负荷集合；

对磨损系数进行归一化，得到归一化的磨损系数指标为：

以上所述，仅为本发明的较佳实施例而已，并非对本发明作任何形式上的限制；凡本行业的普通技术人员均可按说明书所示和以上所述而顺畅地实施本发明；但是，凡熟悉本专业的技术人员在不脱离本发明技术方案范围内，可利用以上所揭示的技术内容而作出的些许更动、修饰与演变的等同变化，均为本发明的等效实施例；同时，凡依据本发明的实质技术对以上实施例所作的任何等同变化的更动、修饰与演变等，均仍属于本发明的技术方案的保护范围之内。

Claims (6)

1.一种齿轮早期故障特征提取方法，其特征在于包括以下步骤：

(1)根据运行可靠性理论，考虑齿轮所处齿轮箱、环境噪声和自身振动信号，采用计入峭度统计因素和负载率的齿轮运行可靠性模型，由此计算齿轮的动态故障概率；

(2)基于齿轮振动信号方差和齿轮振动信号熵的脆弱性指标，衡量齿轮箱中齿轮的状态脆弱性；

(3)在齿轮排序方面，采取分类排序的方式：为了平衡计算时间和计算精度的要求，仅对涉及正常和停止后的大故障概率齿轮(如重过载和恶劣环境下的齿轮)，搜索由这些齿轮停止运行造成的故障率，对故障率进行排序，然后对于较小故障概率的齿轮，使用状态和结构脆弱性指标加权得到的综合脆弱性指标进行快速排序，实现齿轮早期故障特征提取。

2.根据权利要求1所述的齿轮早期故障特征提取方法，其特征在于：在第(1)步中，综合考虑齿轮所处齿轮箱、环境噪声和自身振动信号后，齿轮的可靠性模型下式所示：

3.根据权利要求2所述的齿轮早期故障特征提取方法，其特征在于：在第(2)步中，所述脆弱性指标包括基于负载率方差的齿轮脆弱性指标、基于负载率熵的齿轮脆弱性指标和齿轮的磨损系数指标。

4.根据权利要求3所述的齿轮早期故障特征提取方法，其特征在于：所述基于负载率方差的齿轮脆弱性指标中，齿轮箱正常运行各齿轮负载率方差可表示为：

式中，n和η_i0分别为正常运行时齿轮箱齿轮总数和齿轮i的负载率，齿轮k退出运行后，齿轮箱的振动信号方差可表示为：

式中，η_k,i为齿轮k停止运行后，齿轮i的负载率，

综上所述，定义齿轮k关于负载率方差的脆弱性指标如下式所示：

5.根据权利要求3所述的齿轮早期故障特征提取方法，其特征在于：所述基于负载率熵的齿轮脆弱性指标中，定义齿轮箱正常运行时各齿轮的负载率熵H₀如下式所示：

式中，η_i0(或η_j0)分别为正常运行时齿轮箱齿轮总数和齿轮i(或j)的负载率，

齿轮k停止运行后，齿轮箱齿轮的负载率熵H_k定义如下式所示：

式中，η_k,i(或η_k,j)为齿轮k退出运行后，齿轮i(或j)的负载率；

定义齿轮k基于负载率熵的脆弱性指标如下式所示：

同理，可对下式的指标进行归一化，可得基于负载率熵的归一化指标如下式所示。

式中，

6.根据权利要求3所述的齿轮早期故障特征提取方法，其特征在于：所述齿轮的磨损系数指标中，齿轮k的磨损系数B_k(m,n)定义为：

式中，I^(ij)(m,n)为在“电机—负荷”节点对(i,j)间加入单位驱动后，在齿轮k上引起的压力；W_i为电机节点权重，可以取相应电机组的实时输出应力；W_j为负荷权重，取相应负荷的实时负载；G和L分别为电机和负荷集合；

对磨损系数进行归一化，得到归一化的磨损系数指标为：