一种基于线性回归和TCN的滚动轴承退化趋势预测方法技术

本发明专利技术公开了一种基于线性回归和TCN的滚动轴承退化趋势预测方法，其特征在于，首先通过加速度传感器采集滚动轴承的全生命周期振动数据，计算获得数据的样本熵作为轴承的性能退化指标，对样本熵进行平滑预处理，然后根据线性回归模型和μ+3δ确定轴承的SPT点，分析不同的窗口长度对确定SPT的影响，最后将预处理后的样本熵值输入到训练好的TCN网络中进行退化趋势预测。本发明专利技术能有效且及时地确定轴承的SPT点，TCN能够很好地拟合轴承的性能退化趋势，避免由于发现轴承故障时间过晚而带来的设备不能安全运行的问题；本发明专利技术达到对设备健康管理的方法相比采用循环神经网络预测的方法准确率明显提高，为轴承的退化趋势预测方法提供了新思路。供了新思路。供了新思路。

【技术实现步骤摘要】
一种基于线性回归和TCN的滚动轴承退化趋势预测方法


[0001]本专利技术属于轴承性能退化评估
，具体涉及一种基于线性回归和TCN的滚动轴承退化趋势预测方法。

技术介绍

[0002]随着控制科学和各学科的发展，机械设备朝着复杂、高效和智能化的趋势发展。滚动轴承作为大多数机械设备的关键零部件，特别是像卫星、航天飞机这样的高精尖设备对轴承的振动性能要求特别高，对轴承振动时间序列进行分析为设备的安全运转提供了重要保障。据统计，在含有滚动轴承的机械设备装置中，因滚动轴承故障而导致的机械设备整体故障的概率为50％；如果能对滚动轴承的性能退化作出评估，则可通过预测性维护来避免灾难的发生，保障了生命安全的同时也最大限度的降低设备的维修成本。及时地发现轴承的初始退化时间是退化趋势预测的关键，在轴承的退化趋势预测中，数据驱动的方法已成为预测方法的主流，其中，循环神经网络(RNN)在轴承退化趋势预测中得到广泛应用，但是模型内部都含有大量参数，这些可训练的参数将直接影响模型的性能，且参数没有通用的选择标准。时间卷积网络(TCN)模型的卷积架构在各种任务和数据集上都优于典型的循环网络(例如LSTM和GRU)，且TCN对参数选择的敏感性低于RNN等模型。这使得TCN有着更稳定的性能。
[0003]因此，本文提出一种基于线性回归和TCN的滚动轴承退化趋势预测方法。

技术实现思路

[0004]为了解决上述技术问题，本专利技术设计了一种基于线性回归和TCN的滚动轴承退化趋势预测方法，提取滚动轴承的样本熵作为性能退化指标，将退化指标输入到TCN中实现对滚动轴承的退化趋势预测，提高预测精度。
[0005]为了达到上述技术效果，本专利技术是通过以下技术方案实现的：一种基于线性回归和TCN的滚动轴承退化趋势预测方法，其特征在于：
[0006]Step1：通过加速度传感器采集滚动轴承的全生命周期振动数据，根据采样时间和采样频率将数据进行划分；
[0007]Step2：计算划分后各数据组的样本熵，作为滚动轴承性能退化指标；
[0008]Step3：对样本熵进行平滑预处理；
[0009]Step4：利用线性回归模型确定SPT，若连续多次超出退化阈值则可判定滚动轴承发生了初始退化，并分析不同的窗口长度大小对确定SPT的影响；
[0010]Step5：利用VMD分解确定发生初始退化的数据组数，并对数据进行包络谱分析，判断确定的SPT是否正确；
[0011]Step6：构建TCN模型；
[0012]Step7：根据Step4中确定的SPT将数据划分为训练集和测试集，SPT点之前的数据作为训练集，SPT点之后的数据作为测试集输入到Step6中的模型，实现滚动轴承的退化趋
势预测，并利用评价函数对预测结果进行评估。
[0013]进一步的，所述Step1中的采样时间为1s，采样频率为20000HZ，将轴承全生命周期数据划分为984组，每组含20480个数据。
[0014]进一步的，所述Step2中，划分后的矩阵为其中m＝20480，l＝984。
[0015]进一步的，所述Step3中，平滑方法采用局部加权回归进行平滑，在局部加权回归中，使用给定跨度内的相邻数据点确定平滑值，计算每个数据点的回归权重W
i
，计算方法为其中s
i
为s的相邻元素，d(s)表示跨度中最远值之间的水平距离。
[0016]进一步的，所述Step4中，线性回归模型可表示为：其中，g对应样本熵梯度值，d
i
为对应时间t
i
时的样本熵的值，w为滑动窗口的大小。退化阈值由μ+3δ确定，对于一个近似正态分布的数据，设其均值为μ，方差为σ2，则其数值分布在(μ
‑
3σ,μ+3σ)中的概率为99.73％，多次超出退化阈值次数为5次；将窗口长度取为10、20、30、40、50。
[0017]进一步的，所述Step5中，VMD的分解模态数为7，根据峭度最大原则重构信号后进行包络谱分析。
[0018]进一步的，所述Step6中，构建的TCN模型中包含因果卷积、空洞卷积、残差连接和一维全卷积；所述因果卷积当前时刻预测值y
t
的输出只与当前时刻s
t
和之前的输入相关，而与未来时刻的值s
t+1
无关，可以表示为：y
t
＝f(s1,s2,...,s
t
),t＝1,2,...,T，其中，f(
·
)表示卷积操作；对于一维输入s1,s2,...,s
T
和过滤器f，序列的扩张卷积可以表示为：
[0019]其中，d为膨胀因子，n为滤波器大小，s
x
‑
d
·
i
表示序列元素乘以卷积核中的元素；膨胀因子d越大，输入的范围也就越大，从而加大了卷积网络的感受野；为避免由于深层网络导致的梯度消失问题，在TCN的输出层中引入了残差模块，同时将输入s融合到卷积网络的输出y中，可表示为：y＝G(s,{W
i
})+s，其中，函数G(s,{W
i
})+s表示学习残差映射；通过一维全卷积保证了TCN的输入序列和输出序列的长度保持一致。
[0020]进一步的，所述Step7中，使用TCN网络中完全连接的层来预测输出；预测评价指标
为拟合决定系数R2，其取值范围越接近于1则说明预测效果越好。
[0021]本专利技术的有益效果是：
[0022]1、本专利技术通过线性回归方法和μ+3δ退化阈值自适应法准确且及时地确定了滚动轴承地退化初始时间；
[0023]2、分析了不同窗口长度对确定SPT的影响，避免了由于人工经验选取窗口长度对确定SPT的主观性；
[0024]3、本专利技术通过将样本熵作为性能退化指标输入到TCN网络模型对滚动轴承的退化趋势进行了预测，与目前主流的循环神经网络相比，大大提高了预测的准确性。
附图说明
[0025]为了更清楚地说明本专利技术实施例的技术方案，下面将对实施例描述所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0026]图1为本专利技术的流程图；
[0027]图2为样本熵平滑后反映轴承性能退化趋势图；
[0028]图3为不同窗口长度下确定的SPT；
[0029]图4为对第525组数据进行包络谱分析结果图；
[0030]图5为TCN模型结构图；
[0031]图6为不同预测模型预测退化趋势对比图；
[0032]图7为不同模型的预测准确率对比图。
具体实施方式
[0033]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于线性回归和TCN的滚动轴承退化趋势预测方法，其特征在于：Step1：通过加速度传感器采集滚动轴承的全生命周期振动数据，根据采样时间和采样频率将数据进行划分；Step2：计算划分后各数据组的样本熵，作为滚动轴承性能退化指标；Step3：对样本熵进行平滑预处理；Step4：利用线性回归模型确定SPT，若连续多次超出退化阈值则可判定滚动轴承发生了初始退化，并分析不同的窗口长度大小对确定SPT的影响；Step5：利用VMD分解确定发生初始退化的数据组数，并对数据进行包络谱分析，判断确定的SPT是否正确；Step6：构建TCN模型；Step7：根据Step4中确定的SPT将数据划分为训练集和测试集，SPT点之前的数据作为训练集，SPT点之后的数据作为测试集输入到Step6中的模型，实现滚动轴承的退化趋势预测，并利用评价函数对预测结果进行评估。2.根据权利要求1所述的一种基于线性回归和TCN的滚动轴承退化趋势预测方法，其特征在于：所述Step1中的采样时间为1s，采样频率为20000HZ，将轴承全生命周期数据划分为984组，每组含20480个数据。3.根据权利要求1所述的一种基于线性回归和TCN的滚动轴承退化趋势预测方法，其特征在于：所述Step2中，划分后的矩阵为其中m＝20480，l＝984。4.根据权利要求1所述的一种基于线性回归和TCN的滚动轴承退化趋势预测方法，其特征在于：所述Step3中，平滑方法采用局部加权回归进行平滑，在局部加权回归中，使用给定跨度内的相邻数据点确定平滑值，计算每个数据点的回归权重W
i
，计算方法为其中s
i
为s的相邻元素，d(s)表示跨度中最远值之间的水平距离。5.根据权利要求1所述的一种基于线性回归和TCN的滚动轴承退化趋势预测方法，其特征在于：所述Step4中，线性回归模型可表示为：其中，g对应样本熵梯度值，d
i
为对应时间t
i
时的...

【专利技术属性】
技术研发人员：严帅，邓为权，
申请(专利权)人：昆明理工大学，
类型：发明
国别省市：