[发明专利]一种基于多特征和相空间的滚动轴承寿命动态预测方法

摘要

一种基于多特征和相空间的滚动轴承寿命动态预测方法，先获取滚动轴承的振动信号，提取其时域和频域特征指标，分别进行标准化和滑移处理得到相对特征指标；然后利用PCA技术对多个特征指标进行融合，获得一个综合衰退性能指标序列进行寿命预测；再将历史退化过程的综合衰退性能指标进行相空间重构，在相空间中与当前退化过程进行对比，得到一个预测的失效时间；最后将这个预测的失效时间和历史失效时间匹配组合估计其概率密度分布，计算平均寿命；且通过分析样本的不断积累扩充，得到不同观测时刻的平均寿命，本发明能有效反映滚动轴承退化过程中多因素的作用，不受假设条件的约束，同时减小了预测误差，实现寿命动态和准确地预测。

主权项

1.一种基于多特征和相空间的滚动轴承寿命动态预测方法，其特征在于，包括以下步骤：第一步，获取滚动轴承的振动信号，提取其多个时域和频域特征指标；第二步，将时域和频域特征指标分别进行标准化和滑移处理得到相对特征指标，计算步骤如下：步骤1：任意选取正常期内一段趋势平稳的特征值，计算该段特征值的平均数；步骤2：计算原始特征值与平均数之比，得到一个相对特征指标；步骤3：对相对特征指标进行M点滑移平均处理，M＝5～20，xRRX为原始特征与平均数的比值；xMA为经过滑移平均后的序列，即最终的相对特征指标，计算公式如下：第三步，利用PCA技术对多个特征指标进行融合，获得一个综合衰退性能指标序列进行寿命预测；第四步，将历史退化过程的综合衰退性能指标进行相空间重构，在相空间中与当前退化过程进行对比，得到一个预测的失效时间；第五步，将这个预测的失效时间和历史失效时间匹配组合估计其概率密度分布，计算平均寿命；第六步，通过伸缩时间窗不断积累扩充分析样本，得到不同观测时刻的平均寿命；所述的第四步的具体过程为：步骤1：假设有M组历史退化过程的综合衰退性能指标序列，对这些指标序列进行相空间重构，实现滚动轴承退化过程的动力学轨道在相空间中完全展开；步骤2：相空间中退化轨迹函数的学习，利用径向基(RBF)神经网络对非线性函数优良的逼近能力获得非线性轨迹函数，通过RBF神经网络实现相空间中退化轨迹函数的学习；假设第i组历史退化过程的综合衰退性能指标序列为Xi＝{xi(t1),xi(t2),…,xi(tn)}，其中i＝1,2,…,M、n为第i组指标序列的长度、k＝1,2,…L.L＝n‑(mi‑1)τi经过相空间重构后矢量为：xxik(t)＝{xi(tk),xi(tk+τi),…,xi(tk+(mi‑1)τi)}m和τ是该序列的嵌入维数和延迟时间，将重构后的矢量{xxik(t)|k＝1,2,…,L}和对应的服役时间{tk|k＝1,2,…,L}分别作为输入和输出进行RBF神经网络的训练，得到第i组历史退化过程的相空间轨迹函数，其他各组历史退化过程都经历相同训练得到不同的轨迹函数；步骤3：预测失效时间，假设当前需要预测的退化过程为第M+1次退化过程，在这一观测时刻进行预测对应的运行时间为TM+1，建立的综合衰退性能指标序列为{xM+1(t1),xM+1(t2),…,xM+1(tn)}，对该序列进行相空间重构，把重构矢量输入到训练完成后的RBF经网络中进行学习，获得一组估计的运行时间向量[T1,T2,…TM]，计算该运行时间与实际运行时间的偏差：ei＝(Ti‑TM+1)2i＝1,2,…M则在时刻TM+1处得到一个误差向量e＝[e1,e2,…,eM]，从而得到该时刻退化轨迹与历史退化轨迹的相似程度：tM+1为当前时刻预测的失效时间，ti为滚动轴承历史服役寿命，根据相似程度和历史服役寿命可预测滚动轴承在当前时刻的失效时间：在下一次预测时，通过伸缩窗扩大性能指标数据，重复以上的步骤得到一个新误差向量，将新误差向量与上次计算的误差向量相加重新赋值给新误差向量；所述的第五步的具体过程为：步骤1：输入样本数据，将历史失效时间和当前预测的失效时间进行匹配作为输入样本，假设当前预测是在第K次观测时，得到的样本集为其中m为样本点数；步骤2：扩充样本，原始输入样本Xk通过RBF神经网络扩充后得到新的样本Yk＝{y1,y2,…,yn}，其中n为扩充后样本点数，且n＞m；步骤3：概率估计，σ为平滑因子，对数据归一化后使用核密度法估计这一观测时刻失效时间的概率密度分布：步骤4：用变量t代替上式中的变量y，则得到平均寿命为：