一种基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法
【专利摘要】一种基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法，包括以下步骤：1)采集轴承全寿命周期振动信号；2)利用振动信号计算有效值，基于有效值构建出反映轴承健康状态的指数；3)计算各时刻健康指数变化量，利用矩形窗函数和K?S检验截取出轴承耗损期内的健康指数数据；4)对耗损期的健康指数数据进行拟合分析，构建表征轴承退化过程的状态空间模型，利用当前观测得到的健康指数数据和无迹卡尔曼滤波算法更新模型参数并预测轴承剩余寿命。本发明提供了一种预测精度较高的基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法。
【专利说明】
-种基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法
技术领域
[0001] 本发明属于轴承故障诊断及预测领域，尤其设及一种基于无迹卡尔曼滤波算法的 轴承剩余寿命预测方法。
【背景技术】
[0002] 轴承作为一个关键零部件，在电力、石化、冶金、机械、航空航天等旋转机械中被广 泛的应用，但不幸的是，它也是旋转机械中故障的多发部件。据统计，航空发动机主轴承服 役寿命仅为数百小时，数控机床高速轴承精度寿命为数千小时，一旦运行时间超出服役寿 命，轴承的运行精度会急剧下降，进而导致航空发动机、数控机床等无法正常工作。因此轴 承状态监测、故障诊断W及故障预测一直是近年来的研究重点。考虑到轴承从早期故障发 生，发展直至失效是一个非线性的过程，由此利用贝叶斯理论的非线性滤波算法，如扩展卡 尔曼滤波、粒子滤波等在轴承的故障预测方面得到了快速发展。然而当非线性系统中的非 线性函数化ylor展开式的高阶项无法忽略时，利用扩展卡尔曼滤波进行线性化会使系统产 生较大误差;利用粒子滤波要得到较高精度的预测，需要较多数目的粒子，由此而产生了很 大的计算量，并且粒子经过迭代后会产生粒子退化的问题。
【发明内容】
[0003] 为了克服现有的非线性贝叶斯滤波算法在解决轴承故障预测时剩余寿命预测精 度较低的不足，本发明提供了一种预测精度较高的基于无迹卡尔曼滤波算法的轴承剩余寿 命预测方法。
[0004] 为了解决上述技术问题提供的技术方案为：
[0005] -种基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法，所述方法包括W下步 骤：
[0006] Sl.采集轴承的全寿命周期振动信号；
[0007] S2.利用振动信号计算有效值，基于有效值构建出反映轴承健康状态的指数；
[000引S3.计算各时刻轴承健康指数的变化量，并利用矩形窗函数进行截断，运用K-S检 验检验矩形窗内健康指数的变化量是否符合正态分布，随时间推移不断推进矩形窗，K-S检 验结果将轴承寿命周期分为磨合期、有效寿命期和耗损期=个时期，截取出轴承耗损期的 健康指数数据；
[0009] S4.拟合分析轴承耗损期的健康指数数据，构建表征其退化过程的非线性状态空 间模型，利用当前观测到的健康指数数据和无迹卡尔曼滤波算法更新模型参数，并预测剩 余寿命，过程如下：
[0010] 对耗损期的健康指数数据，拟合分析构建如下的退化模型：
[0011] HI(k) =a ? exp(b ? k)+c ? exp(d ? k) (I)
[001^ 上式中，HKk)为轴承在k时刻的健康指数，k为时间参数，a，b，c，d为退化模型参数 初始值，基于该退化模型构建状态方程：
[0013]
(2)
[0014] 上式中，ak，bk，Ck，dk，和祉-1，bk-i，Ck-I，dk-i分别为在k时刻和k-1时刻状态变量a，b， c，d的值，N(0,〇2)为高斯白噪声；
[0015] 同时构建观测方程：
[0016] HIk =化? exp(bk ? k)+ck ? e邱（dk ? k)+vk (3)
[0017]上式中，HIk为在k时亥讚康指数的观测值，Vk为在k时刻的观测噪声；
[0018] 利用无迹卡尔曼滤波算法更新状态方程和观测方程参数至k时刻，按公式（1)计算 k+1时刻的健康指数HI化+1):
[0019] HI(k+l) =ak ? e邱（bk ?化+1))+ck ? exp(dk ?化+1)) (4)
[0020] 上式中，1 = 1，2，...，-;计算使得不等式(5)成立的1值，并记录1的最小值为加寸 刻预测的轴承剩余寿命；
[0021] HI 化+1)〉故障阔值 （5)。
[0022] 进一步，所述S2中，对Sl所得的轴承全寿命周期振动信号，构建健康指数，过程如 下：
[0023] 设第k时刻振动信号Xi，其包含N个采样点，则该时刻对应的有效值x_rms化)为：
[0024]
(6)
[0025] 健康指数包含水平和垂直两个方向的信息，其由下式计算得到：
[0026]
(7)
[0027] 上式中，x_rms和y_rms分别为在水平振动信号和垂直振动信号上计算得到的有效 值。
[0028] 再进一步，所述S3中，对S2所得的健康指数，截取轴承损耗期的健康指数数据；
[0029] 计算各时刻有效值的变化量A HI，其由下式计算得到：
[0030] AHI(k)=HI 化）-HI(k-l) (8)
[0031 ]上式中，HK k)和HI化-1)分别为k时刻和k-1时刻健康指数；
[0032] 用矩形窗函数截取健康指数变化量，用Fn(Z)表示截取出健康指数变化量的累计 分布函数，F(Z)表示高斯分布的累计概率分布函数，贝化-S检验统计量是F(Z)和Fn(Z)之间 的最大偏差Dn:
[0033] Dn=max{ IF(Z)-Fn(Z) I } (9)
[0034] 若对于每一个Z来说，Fn(Z)和F(Z)都十分接近，则表明矩形框所截取的健康指数 变化量分布与高斯分布的拟合度很高，不断推进矩形窗，根据K-S检验的结果，将轴承寿命 周期分为磨合期、有效寿命期和耗损期=个时期，从而截取出轴承在耗损期的数据。
[0035] 本发明的技术构思为:通过采集轴承振动信号，基于对振动信号有效值的计算构 建健康指数，利用矩形窗函数和K-S检验截取出轴承耗损期的健康指数数据。通过拟合分析 轴承在耗损期内的健康指数数据，构建用于描述轴承退化过程的状态空间模型，利用当前 观测到的健康指数和无迹卡尔曼滤波算法更新模型参数，并预测剩余寿命。
[0036] 本发明的有益效果为:准确截取出轴承耗损期的数据，剩余寿命预测精度较高。
【附图说明】
[0037] 图1为基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法流程图；
[0038] 图2为轴承全寿命周期健康指数示意图；
[0039] 图3为轴承健康指数变化量图；
[0040] 图4为矩形窗截取数据示意图；
[0041 ]图5为轴承健康指数变化量经矩形窗框截取后的K-S检验结果图；
[0042] 图6为轴承在耗损期的健康指数数据；
[0043] 图7为轴承在耗损期阶段的剩余寿命预测。
【具体实施方式】
[0044] 下面结合附图对本发明作进一步描述。
[0045] 参照图1~图7,一种基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法，所述方 法包括W下步骤：
[0046] Sl.采集轴承的全寿命周期振动信号；
[0047] S2.利用振动信号计算有效值，基于有效值构建出反映轴承健康状态的指数，方便 后续步骤利用该指数进行剩余寿命预测；
[0048] S3.计算各时刻轴承健康指数的变化量，并利用矩形窗函数进行截断，运用K-S检 验检验矩形窗内健康指数的变化量是否符合正态分布，随时间推移不断推进矩形窗，K-S检 验结果将轴承寿命周期分为磨合期、有效寿命期和耗损期=个时期，截取出轴承耗损期的 健康指数数据；
[0049] S4.拟合分析轴承耗损期的健康指数数据，构建表征其退化过程的非线性状态空 间模型，利用当前观测到的健康指数数据和无迹卡尔曼滤波算法更新模型参数，并预测剩 余寿命。
[0050] 所述S2中，对Sl所得的轴承全寿命周期振动信号，建立健康指数；
[0051] 设第k时刻振动信号XI，其包含N个采样点，则该时刻对应的有效值x_rms化)为：
[0化2]
(6)
[0053] 健康指数包含水平和垂直两个方向的信息，其由下式计算得到：
[0054]
巧）
[0055] 上式中，x_rms和y_rms分别为在水平振动信号和垂直振动信号上计算得到的有效 值；
[0056] 所述S3中，对S2所得的健康指数，截取轴承损耗期的健康指数数据；
[0057] 计算各时刻有效值的变化量A HI，其由下式计算得到：
[005引 AHKlO=HI 化）-HKk-I) (8)
[0059] 上式中，HKk)和HKk-1)分别为k时刻和k-1时刻健康指数；
[0060] 用矩形窗函数截取健康指数变化量，用Fn(Z)表示截取出健康指数变化量的累计 分布函数，F (Z)表示高斯分布的累计概率分布函数，贝化-S检验统计量是F (Z)和Fn (Z)之间 的最大偏差Dn:
[0061] Dn=max{ IF(Z)-Fn(Z) I } (9)
[0062] 若对于每一个Z来说，Fn(Z)和F(Z)都十分接近，则表明矩形框所截取的健康指数 变化量分布与高斯分布的拟合度很高，不断推进矩形窗，如附图4所示，根据K-S检验的结 果，将轴承寿命周期分为磨合期、有效寿命期和耗损期=个时期，如附图5所示，从而截取出 轴承在耗损期的数据；
[0063] 所述S4中，拟合分析轴承耗损期的健康指数数据，构建退化模型并建立状态空间 模型，利用当前观测到的健康指数数据和无迹卡尔曼滤波算法更新模型参数，并预测剩余 寿命，过程如下：
[0064] 对耗损期的健康指数数据，拟合分析构建如下的退化模型：
[00化]HI(k) =a ? exp(b ? k)+c ? exp(d ? k) (I)
[0066] 上式中，HKk)为轴承在k时刻的健康指数，k为时间参数，a，b，c，d为耗损期轴承健 康指数拟合出的退化模型参数，基于该退化模型构建状态方程：
[0067]
(2)
[006引上式中，ak，bk，Ck，dk，和祉-1，bk-1，Ck-I，dk-1分别为在k时刻和k-1时刻状态变量a，b， c，d的值，N(0,02)为高斯白噪声；
[0069] 同时构建观测方程：
[0070] HIk =化? exp(bk ? k)+ck ? e邱（dk ? k)+vk (3)
[0071] 上式中，HIk为在k时刻健康指数的观测值，Vk为在k时刻的观测噪声；
[0072] 利用无迹卡尔曼滤波算法更新状态方程和观测方程参数至k时刻，按公式（1)计算 k+1时刻的健康指数HI化+1):
[0073] HI(k+l) =ak ? e邱（bk ?化+1))+ck ? exp(dk ?化+1)) (4)
[0074] 上式中，1 = 1，2，...，-;计算使得不等式(5)成立的1值，并记录1的最小值为加寸 刻预测的轴承剩余寿命；
[00巧]HI化+1)〉故障阔值 （5)。
[0076] 本实施例利用PR0N0STIA平台轴承全周期寿命数据对基于无迹卡尔曼滤波的轴承 剩余寿命预测方法进行验证。具体过程如下：
[0077] (1)采集轴承的振动信号。通过加速度传感器采集水平方向和垂直方向的振动信 号，信号每IOs采集一次，每一次采集时长为0.1 s，数据采样频率为25.6曲Z;
[0078] (2)利用振动信号计算有效值，基于有效值构建出反映轴承健康状态的指数，方便 后续步骤利用该指数进行剩余寿命预测，构建出轴承健康指数反应其健康状态如附图2所 示；
[0079] (3)计算各时刻健康指数变化量，如附图3所示，并利用矩形窗函数进行截断，运用 K-S检验检验矩形窗内健康指数的变化量是否符合正态分布，随时间推移不断推进矩形框， 根据K-S检验的结果，将轴承寿命周期分为磨合期、有效寿命期和耗损期=个时期，从而截 取出耗损期的寿命进行剩余寿命预测，截取出轴承耗损期的健康指数数据，如附图6所示；
[0080] (4)利用双指数函数拟合耗损期轴承数据，构建退化模型并建立状态空间模型，利 用当前观测到的健康指数数据和无迹卡尔曼滤波算法更新模型参数，并预测剩余寿命。利 用无迹卡尔曼滤波算法更新模型参数和预测剩余寿命，建立剩余寿命预测模型为：
[0081] HI(k+l)=￡ik ? e邱（bk ?化+l))+ck ? exp(dk ?化+1)) (4)
[0082] 上式中，1 = 1，2，...，-;计算使得不等式(5)成立的1值，并记录1的最小值为加寸 刻预测的轴承剩余寿命；
[0083] HI 化+1)〉故障阔值 （5)。
[0084] 附图7表示轴承数据的预测曲线，从曲线中可W看出，一开始由于数据不足，预测 曲线与实际寿命曲线的误差较大，随着数据量的不断增大，最终的预测曲线与实际曲线相 吻合。有效的验证了无迹卡尔曼滤波在轴承故障检测和故障预测中的可行性。
【主权项】
1. 一种基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法，其特征在于:所述方法包 括以下步骤：51. 采集轴承的全寿命周期振动信号；52. 利用振动信号计算有效值，基于有效值构建出反映轴承健康状态的指数；53. 计算各时刻轴承健康指数的变化量，并利用矩形窗函数进行截断，运用K-S检验矩 形窗内健康指数的变化量是否符合正态分布，随时间推移不断推进矩形窗，K-S检验结果将 轴承寿命周期分为磨合期、有效寿命期和耗损期三个时期，截取出轴承耗损期的健康指数 数据；54. 拟合分析轴承耗损期的健康指数数据，构建表征其退化过程的非线性状态空间模 型，利用当前观测到的健康指数数据和无迹卡尔曼滤波算法更新模型参数，并预测剩余寿 命，过程如下： 对耗损期内的健康指数数据，通过拟合分析构建如下的退化模型： HI(k)=a · exp(b · k)+c · exp(d · k) (I) 上式中，HI (k)为轴承在k时刻的健康指数，k为时间参数，a，b，c，d为退化模型参数初始 值，基于该退化模型构建状态方程：(2) 上式中，ak，bk，ck，dk，和ak-i，bk-i，ck-i，dk-i分别为在k时刻和k-Ι时刻状态变量a，b，c，d的 值，Ν(0, σ2)为高斯白噪声； 同时构建观测方程： HIk = ak · exp(bk · k)+ck · exp(dk · k)+vk (3) 上式中，HIk为在k时刻健康指数的观测值，Vk为在k时刻的观测噪声； 利用无迹卡尔曼滤波算法更新状态方程和观测方程参数至k时刻，按公式（1)计算k+1 时刻的健康指数HI( k+1): HI(k+l)=ak · exp(bk · (k+l))+ck · exp(dk · (k+1)) (4) 上式中，1 = 1,2, . . .;计算使得不等式(5)成立的I值，并记录I的最小值为k时刻预 测的轴承剩余寿命； HI (k+Ι )>故障阈值 （5)。2. 如权利要求1所述的一种基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法，其特 征在于:所述S2中，对Sl所得的轴承全寿命周期振动信号，构建健康指数，过程如下： 设第k时刻振动信号Xi，其包含N个采样点，则该时刻对应的有效值x_rms(k)为：?(,) 健康指数HI包含水平和垂直两个方向的信息，其由下式计算得到： 0) 上式中，X_rmi^Py_rmS分别为在水平振动信号和垂直振动信号上计算得到的有效值。3.如权利要求1或2所述的一种基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法，其 特征在于:所述S3中，对S2所得的健康指数，截取轴承耗损期的健康指数数据； 计算各时刻有效值的变化量A HI，其由下式计算得到： AHI(k)=HI(k)-HI(k-l) (8) 上式中，HI (k)和HKk-I)分别为k时刻和k-1时刻健康指数； 用矩形窗函数截取健康指数变化量，用Fn(Z)表示截取出健康指数变化量的累计分布函 数，F(Z)表示高斯分布的累计概率分布函数，则K-S检验统计量是F(zWPFn(z)之间的最大 偏差D n: Dn=max{ IF(z)-Fn(z) I} (9) 若对于每一个z来说，Fn(Z)和F(Z)都十分接近，则表明矩形框所截取的健康指数变化量 分布与高斯分布的拟合度很高，不断推进矩形窗，根据K-S检验的结果，将轴承寿命周期分 为磨合期、有效寿命期和耗损期三个时期，从而截取出轴承在耗损期的数据。
【文档编号】G06F17/50GK106021719SQ201610334109
【公开日】2016年10月12日
【申请日】2016年5月19日
【发明人】金晓航, 阙子俊, 孙毅
【申请人】浙江工业大学