一种回转支承使用寿命的预测方法技术

一种回转支承使用寿命的预测方法，其特征包括如下步骤：提取回转支承的温度、扭矩和振动原始信号；在时域内选取信号的有量纲和无量纲特征值并选取每个信号对应的敏感特征值；将每个信号的敏感特征值做主成分分析（PCA）降维处理，分别形成回转支承的温度PCA性能衰退指标、扭矩PCA性能衰退指标和振动PCA性能衰退指标；通过自组织特征映射神经网络（SOM）寻求回转支承的寿命周期状态分布，完成输入样本类别的划分。本发明专利技术实现了回转支承的衰退特征指标的建立以及寿命周期状态的分布，从而有效的进行剩余寿命的预测，大大减少因回转支承寿命结束的突发性带来的巨大损失。

【技术实现步骤摘要】
一种回转支承使用寿命的预测方法
本专利技术涉及一种回转支承使用寿命的预测方法，是一种基于智能算法的寿命预测方法。
技术介绍
回转支承是最近几十年发展起来，实现相对回转运动的新型机械部件，由内外圈和滚动体构成，一般需要同时承受轴向力、径向力、倾覆力矩。回转支承的直径通常大于400mm，回转速度运转速度较低，通常在0.1-50r/min，广泛应用于建筑运输设备、风力发电机以及其他机械工程领域，通常安装在大型设备的关键部位。由于回转支承尺寸一般比较大，价格比较昂贵，且储存不便，通常不存备件；载荷作用形式多样、运行环境恶劣，维修难度大，需将几十吨甚至几百吨的东西提升一定高度后才能维修而加大，费用高、停机时间长，因此对回转支承性能衰退评价和寿命预测方法的研究具有重要意义。轴承类部件寿命预测方法主要有基于力学、基于概率统计、基于信息新技术等。传统的基于力学和概率统计的轴承类寿命预测方法已有很多研究及应用。但是集成度和复杂度日益增加的设备的数学模型或失效机理模型难以准确给出，因此基于力学模型的剩余寿命预测方法受到一定限制。虽然基于概率统计剩余寿命预测方法的预测结果更能反映寿命预测的一般规律，但是需要大量试验和数据的积累，然而大型回转支承工况环境复杂、低速重载，样本数量少，基于概率统计剩余寿命预测方法依旧不适用。基于人工智能的寿命预测方法还不完善，是目前寿命预测领域的研究热点。由于回转支承的工作环境复杂，性能衰退特征具有一定的非线性、非平稳性，单一性能参数不能全面有效地评价回转支承性能衰退的全过程。因此本专利选取回转支承的温度、扭矩和振动信号时域内的特征值建立其衰退性能指标。由于大型回转支承样本数量少，为实现较高的预测精度，本专利选取针对小样本寿命预测的支持向量回归的智能算法。近年来，SVR也较多的应用于机械设备状态监测和剩余预测问题，很多研究人员采用SVR算法预测轴承的剩余寿命，取得了较好的预测效果。在寿命预测之前本专利先采用SOM算法对回转支承的寿命状态进行分类，在实际的运用过程中，根据数据所属的类别选择对应的预测模型，这样可以缩短预测时间，同时提高预测精度，另外随着数据信息的不断更新，预测模型也具有一定的应变能力。因此，本专利提出了基于SOM-MSVR的回转支承寿命预测方法。
技术实现思路
为解决上述问题本专利技术提供了一种基于智能算法的回转支承寿命预测方法。目的是实现回转支承性能衰退特征的建立以及寿命周期的分布，从而对回转支承每个状态进行有针对性的寿命预测，不但可以提高预测精度，还可以在实际使用过程中不断完善预测模型。本专利技术为解决以上技术问题，所采用的技术方案是：一种基于智能算法的回转支承寿命预测方法，包括如下步骤：步骤一：提取回转支承的温度、扭矩和振动原始信号，按以下公式计算提取每个信号时域内有量纲和无量纲特征值：Xmax＝max{|xn|}(n＝1,2,...,N)式中：Xmax为数列{xn}的最大值；xn为原始信号第n点的值，n＝1,2,…,N，N为数据采集点的个数。X(p-p)＝max(xn)-min(xn)式中：X(p-p)为数列{xn}的峰峰值。式中：为数列{xn}的方差；为为数列{xn}的均值。式中：Xrms为数列{xn}的均方根值。式中：Xr为数列{xn}的方根幅值。式中：为数列{xn}的绝对平均幅值。式中：β为数列{xn}的峭度。式中：Sf为数列{xn}的波形指标。式中：Kv为数列{xn}的峭度指标。式中：If为数列{xn}的脉冲指标。式中：CLf为数列{xn}的裕度指标。步骤二：对每个特征值按下式做相关度计算分析，选取每个原始信号的敏感特征值。式中r表示相关系数，设置相关系数阈值为0.85，当相关系数大于0.85的特征值作为敏感特征值并被保留，否则剔除。步骤三：将每个信号选取的敏感特征值做PCA降维处理，分别形成回转支承的温度PCA性能衰退指标、扭矩PCA性能衰退指标和振动PCA性能衰退指标。以温度信号为例，具体计算步骤如下所示：对温度信号时域内的敏感特征值提取后组成N×9数据矩阵如下式所示：式中，T为温度信号的特征值矩阵；tnj指第n个时刻的第j个特征值，j＝1,2,...,9，n＝1,2,...,N,N为数据采样点的个数。根据主成分分析原理，通过变换将9个敏感特征值变量组成新的综合变量，用矩阵表示为：式中：Yt1,Yt2,...,Ytp为构造的新变量(p≤9)，且互不相关方差依次减小，依此为第一主成分、第二主成分...第p个主成分；T*为数据矩阵T进行标准归一化处理得到；αpj是主成分系数，αp1+αp2+,...,+αp9＝1,系数αpj是矩阵T*协方差阵E的特征值λp所对应的特征向量αp。选取第一主分量代替原始信号，通过融合温度信号时域内特征值得到可反映原始信息的综合指标，称为温度PCA衰退指标，则第n个点的温度PCA衰退指标为：式中：α1j是温度信号系数矩阵，j＝1,2,...9。同理，确定第n点的扭矩PCA衰退指标为(2-11)所示：Yq(n)＝(Tq1*(n),Tq2*(n),…,Tq6*(n))·β1bT式中：Tq*为Tq标准归一化得到的矩阵，Tq为扭矩信号的特征值矩阵；β1b是扭矩信号系数矩阵，b＝1,2,...6。确定第n点的振动PCA衰退指标为：Ya(n)＝(A1*(n),A2*(n),…,A6*(n))·γ1cT式中：Ac*为Ac标准归一化得到的矩阵，Ac(n)为振动信号的特征值矩阵；γ1c是振动信号系数矩阵，c＝1,2,...,6。步骤四：温度、扭矩、振动PCA性能衰退指标通过SOM聚类寻求回转支承的寿命周期分布状态，完成输入样本类别的划分。SOM算法聚类的具体步骤如图4所示。步骤五：利用PSO算法优化SVR内部核函数的参数g和惩罚因子C，从而优化回转支承的MSVR剩余寿命预测模型，具体步骤如图5所示。步骤六：对每一个子类分别建立MSVR预测模型，本专利选用RBF核函数作为SVR核函数：式中：控制核函数径向作用范围的是变量g，函数中心是Yc。MSVR预测模型如图6所示，u指输入层变量个数，本文有温度、扭矩、振动PCA性能衰退指标三个变量，所以u＝3。K1、K2、Km是SVR核函数层的节点数。W1、W2、Wu为权重值。第n点的温度、扭矩、振动信号的PCA性能衰退指标与剩余寿命的函数关系如下式所示：Z(n)＝F[Yt(n),Yq(n),Ya(n)]式中，Z(n)表示回转支承的剩余寿命值；Yt(n)、Yq(n)和Ya(n)分别表示第n点回转支承的温度PCA值、扭矩PCA值和振动PCA值。回归函数为：F(Yt(n),Yq(n),Ya(n))＝<W·K(Yt(n),Yq(n),Ya(n))>+b式中：ω为权值矢量，b是偏置门限，K(u)是SVR中的核函数。在网络构建中，随机选取M试验数据，选取前m个数据作为训练集，剩余的(M-m)个样本作为测试集，训练集和预测集按照如下公式输入：式中：Xtrain是训练输入集，Ytrain是训练输出集，Xtest是测试输入集，Z(m)表示对应点的实际剩余生命值。由均方根误差RMSE对模型结果进行评定，其公式如下：式中：zi表示真实值，表示评估值。步骤七：在线预测时，实时监测回转支承的温度、扭矩、振动信号数据，其对应的性能衰退指本文档来自技高网...

【技术保护点】
一种回转支承使用寿命的预测方法，其特征包括如下步骤：步骤一：采用温度传感器、扭矩传感器和加速度传感器分别监测回转支承的温度信号、扭矩信号和振动信号；步骤二：分别提取步骤一中测量到的温度信号、扭矩信号和振动信号时域内有量纲和无量纲特征值并做相关度计算分析，选取每个原始信号对应的敏感特征值；步骤三：将每个信号的敏感特征值做PCA降维处理，分别形成回转支承的温度PCA性能衰退指标、扭矩PCA性能衰退指标和振动PCA性能衰退指标；步骤四：温度、扭矩、振动PCA性能衰退指标通过SOM聚类寻求回转支承寿命周期状态的分布，完成输入样本类别的划分；步骤五：采用PSO算法优化SVR内部参数，从而优化回转支承的MSVR剩余寿命预测模型。步骤六：对每一个子类分别建立MSVR预测模型，则第n点的温度、扭矩、振动信号的PCA性能衰退指标与剩余寿命的函数关系如公式所示：Z(n)＝F[Y

【技术特征摘要】
1.一种回转支承使用寿命的预测方法，其特征包括如下步骤：步骤一：采用温度传感器、扭矩传感器和加速度传感器分别监测回转支承的温度信号、扭矩信号和振动信号；步骤二：分别提取步骤一中测量到的温度信号、扭矩信号和振动信号时域内有量纲和无量纲特征值并做相关度计算分析，选取每个原始信号对应的敏感特征值；步骤三：将每个信号的敏感特征值做PCA降维处理，分别形成回转支承的温度PCA性能衰退指标、扭矩PCA性能衰退指标和振动PCA性能衰退指标；步骤四：温度、扭矩、振动PCA性能衰退指标通过SOM聚类寻求回转支承寿命周期状态的分布，完成输入样本类别的划分；步骤五：采用PSO算法优化SVR内部参数，从而优化回转支承的MSVR剩余寿命预测模型。步骤六：对每一个子类分别建立MSVR预测模型，则第n点的温度、扭矩、振动信号的PCA性能衰退指标与剩余寿命的函数关系如公式所示：Z(n)＝F[Yt(n),Yq(n),Ya(n)]式中，Z(n)表示回转支承的剩余寿命值；Yt(n)、Yq(n)和Ya(n)分别表示第n点回转支承的温度PCA值、扭矩PCA值和振动PCA值。步骤七：在线预测时，实时监测回转支承的温度、扭矩、振动信号，通过SOM模型判定回转支承此时寿命状态所属类别，选用相应的预测模型预测k时刻点的剩余寿命T(k)；并不断在线更新数据，实现在线实时预测并不断完善预测模型。2.根据权利要求1所述的回转支承使用寿命的预测方法，其特征所述步骤一中在回转支承的四周每隔90°安装一个温度传感器，另外设置一个温度传感器悬在空中用于测量室温，以便去除室环境温度的影响，提取回转支承的温度信号；所述扭矩传感器安装在马达泵和小齿轮之间，提取回转支承的扭矩信号；所述加速度信号传感器通过磁座固定方式每隔90°安装在靠近滚...

【专利技术属性】
技术研发人员：王华，汤燕，庞碧涛，黄筱调，
申请(专利权)人：南京工业大学，洛阳LYC轴承有限公司，
类型：发明
国别省市：江苏,32