一种混合分解与提取的智能旋转机械故障诊断方法技术

本发明专利技术公开了一种混合分解与提取的智能旋转机械故障诊断方法，包括：获取旋转机械的原始振动序列进行去噪处理，得到已去噪处理的振动序列；采用不同消失矩的母小波包函数，对已去噪处理的振动序列进行分解操作，得到各自若干个振动子序列；采用主元分析法，对每个振动子序列进行特征提取，得到旋转机械的特征参数矩阵；获取旋转机械训练样本，且已知训练样本的故障类型，使用训练样本训练支持向量机模型得到旋转机械故障诊断分类器；实时采集待检测旋转机械的原始振动序列，获取其特征参数矩阵，并输入到旋转机械故障诊断分类器，从而判断待检测旋转机械的故障类型。本发明专利技术提高了对旋转机械故障诊断的时效性和精度。

【技术实现步骤摘要】
一种混合分解与提取的智能旋转机械故障诊断方法
本专利技术涉及机械系统故障辨识与预测领域，特别是涉及一种混合分解与提取的智能旋转机械故障诊断方法。
技术介绍
随着科技的不断进步和中国制造2025的提出，各种机械装备广泛使用在各行各业，其中旋转机械在这些机械装备中占据了重要地位。但是旋转机械在长期使用过程中，因为一直存在不同程度的振动工况，极其容易出现不同程度的损耗而产生各种机械故障，如果发现不及时导致机械故障不断积累，轻则出现停机事故导致经济损耗，重则产生安全隐患造成安全事故。目前对旋转机械系统的状态监测与诊断最常用的方法是采集旋转机械系统的随机振动信号，然后对随机振动信号进行时频分析。但由于所测取的振动信号不但有故障信号，而且还包含其它机械零部件的振动信号，并且旋转机械在早期的故障产生阶段，其故障特征信息很微弱，信噪比很小，难以提取。
技术实现思路
本专利技术提出一种混合分解与提取的智能旋转机械故障诊断方法，可以更好更及时的检测出旋转机械的故障类型，即发生故障的部位和故障名称，提高故障诊断的准确性和时效性。为实现上述技术目的，本专利技术的技术方案为：一种混合分解与提取的智能旋转机械故障诊断方法，包括以下步骤：步骤1，获取旋转机械的原始振动序列，对原始振动序列进行去噪处理，得到已去噪处理的振动序列；步骤2，采用若干个具有不同消失矩的母小波包函数，对已去噪处理的振动序列进行分解操作；已去噪处理的振动序列在采用不同消失矩的母小波包函数分解后分别得到各自若干个振动子序列；步骤3，采用主元分析法，对每个振动子序列进行特征提取，得到特征参数矩阵；步骤4，按照步骤1-3对已知故障类型的旋转机械提取特征参数矩阵，作为训练样本；以旋转机械的特征参数矩阵训练样本为输入、故障类型为输出，训练支持向量机模型，得到旋转机械故障诊断分类器；步骤5，实时采集待检测旋转机械的原始振动序列，按照步骤1-3获取待检测旋转机械的特征参数矩阵；将待检测旋转机械的特征参数矩阵输入到旋转机械故障诊断分类器，旋转机械故障诊断分类器输出待检测旋转机械的故障类型。本方案通过采用小波包分解方法，使用多个不同的消失矩对振动序列进行分解，在采用不同消失矩的母小波包函数分解后分别得到多个不同的振动子序列；采用主元分析法，针对每个振动子序列进行特征提取，得到旋转机械的特征参数矩阵；这种小波包分解与主元分析法特征提取相结合的方法，使得到的特征参数矩阵在保证广度和深度的前提下，不但保留了主要特征信息，而且降低了输入支持向量机模型的维度，在后续对待检测旋转机械故障诊断分类器对故障诊断的计算时间，提升了分类器的时效性；而且以训练样本的该特征参数矩阵和故障类型训练的旋转机械故障诊断分类器，对待检测旋转机械的故障诊断精度更高。进一步地，所述支持向量机模型的输入层包括1个节点，输出层包括1个节点，核函数为径向基函数，支持向量机模型的惩罚系数和核函数参数采用狼群-模拟退火混合算法进行优化选择。采用狼群-模拟退火算法对支持向量机的参数进行优化，这样能够有效提高优化参数得到的旋转机械故障诊断分类器在对特征序列进行分类时的精度，增加了故障诊断的准确性和可靠性。进一步地，所述支持向量机模型的惩罚系数和核函数参数采用狼群-模拟退火混合算法进行优化选择的具体过程为：步骤4.1，以个体狼位置作为支持向量机分类模型的惩罚系数、核函数参数，随机初始化狼群中的每只个体狼并设置狼群参数；狼群规模取值范围为：[50，200]，步长因子取值范围为：[50，120]，探狼比例因子取值范围为：[2，6]，最大游走次数取值范围为：[10，40]，距离判定因子取值范围为：[40，100]，最大奔袭次数取值范围为：[4，16]，更新比例因子取值范围为：[2，6]，最大迭代次数取值范围为：[500，1000]，最大搜索精度取值范围为：[0.001，0.005]；设定模拟退火算法的退火初始温度为100、退火速率为退火迭代次数t2＝1，当前温度下最大退火循环次数为Lmax＝6；步骤4.2，设定适应度函数，并确定初始最优头狼位置和迭代次数t，t＝1；将个体狼位置的对应的SVM模型惩罚系数、核函数参数代入支持向量机模型中，将所有训练样本输入支持向量机模型获得旋转机械故障诊断分类器，旋转机械故障诊断分类器输出的预测故障类型正确的数量占所有训练样本数量的比例为诊断准确率，取诊断准确率作为适应度函数；步骤4.3，依次对所有人工狼进行游走行为、奔袭行为、围攻行为，按照个体狼的适应度函数更新狼群，获得更新后的最优头狼位置；步骤4.4，判断是否到达优化精度要求或最大迭代次数，若没有到达，令t＝t+1转至步骤4.5，若到达，转至步骤4.7；步骤4.5，对本次迭代中的最优头狼个体进行模拟退火操作，在得到的头狼位置ai邻域内随机选择新的位置aj并计算两者适应度之差ΔF＝F(ai)-F(aj)，计算选择概率P＝exp(-ΔF/Ti)，Ti为当前温度；如果P＞random[0，1)，则将当前头狼位置由ai替换为aj，并以aj作为下次寻优的开始，否则以ai开始下一次寻优；步骤4.6，令t2＝t2+1，按照进行降温退火，若t2＜Lmax，转至步骤4.5，否则，转至步骤4.3；步骤4.7，输出最优头狼位置对应的支持向量机分类模型所使用的惩罚系数和核函数参数。进一步地，步骤2采用的母小波包函数的消失矩分别为：db5、db7和db9。进一步地，所述步骤3的具体过程为：步骤3.1，所有振动子序列转置组成序列矩阵Q；步骤3.2，将序列矩阵Q中的每个振动子序列进行零均值化处理；步骤3.3，求出零均值化处理后的序列矩阵的协方差矩阵，得到协方差矩阵的特征值和对应的特征向量；步骤3.4，将特征向量按对应特征值大小按序排列，取前k行组成矩阵Q'，该矩阵Q'即为旋转机械的特征参数矩阵。进一步地，所述对原始振动序列进行去噪处理的具体过程为：利用傅里叶变换将原始振动序列分解为低频和高频信号，利用滤波器滤除高频信号，得到的低频信号经傅里叶逆变换后得到已去噪处理的振动序列。进一步地，旋转机械的故障类型包括转子不平衡、转子与联轴器不对中、转轴弯曲、油膜涡动及振荡、连接松动、磨碰喘振。有益效果本方案通过采用小波包分解方法，使用多个不同的消失矩对振动序列进行分解，每个消失矩对应多个不同的振动子序列；采用主元分析法，针对每个振动子序列进行特征提取，得到旋转机械的特征参数矩阵；这种小波包分解与主元分析法特征提取相结合的方法，使得到的特征参数矩阵的广度和深度；从而以训练样本的该特征参数矩阵和故障类型训练的旋转机械故障诊断分类器，对待检测旋转机械的故障诊断精度更高。附图说明图1为本专利技术混合分解与提取的智能旋转机械故障诊断方法的系统框架图。具体实施方式下面结合最佳实施方式与附图对本专利技术的混合分解与提取的智能旋转机械故障诊断方法进行解释说明，如图1所示，具体包括以下步骤。步骤1，获取旋转机械的原始振动序列{Xt}，对原始振动序列进行去噪处理，得到已去噪处理的振动序列。原始振动序列{Xt}的去噪处理过程：利用傅里叶变换将原始振动序列分解为低频和高频信号，利用滤波器滤去高频信号，得到的低频信号经傅里叶逆变换后得到已去噪处理的振动序列{Xt′}，此时的已去噪处理的振动序列更加平稳。步骤2，采用消失矩为本文档来自技高网...

【技术保护点】
1.一种混合分解与提取的智能旋转机械故障诊断方法，其特征在于，包括以下步骤：步骤1，获取旋转机械的原始振动序列，对原始振动序列进行去噪处理，得到已去噪处理的振动序列；步骤2，采用若干个具有不同消失矩的母小波包函数，对已去噪处理的振动序列进行分解操作；已去噪处理的振动序列在采用不同消失矩的母小波包函数分解后分别得到各自若干个振动子序列；步骤3，采用主元分析法，对每个振动子序列进行特征提取，得到特征参数矩阵；步骤4，按照步骤1‑3对已知故障类型的旋转机械提取特征参数矩阵，作为训练样本；以旋转机械的特征参数矩阵训练样本为输入、故障类型为输出，训练支持向量机模型，得到旋转机械故障诊断分类器；步骤5，实时采集待检测旋转机械的原始振动序列，按照步骤1‑3获取待检测旋转机械的特征参数矩阵；将待检测旋转机械的特征参数矩阵输入到旋转机械故障诊断分类器，旋转机械故障诊断分类器输出待检测旋转机械的故障类型。

【技术特征摘要】
1.一种混合分解与提取的智能旋转机械故障诊断方法，其特征在于，包括以下步骤：步骤1，获取旋转机械的原始振动序列，对原始振动序列进行去噪处理，得到已去噪处理的振动序列；步骤2，采用若干个具有不同消失矩的母小波包函数，对已去噪处理的振动序列进行分解操作；已去噪处理的振动序列在采用不同消失矩的母小波包函数分解后分别得到各自若干个振动子序列；步骤3，采用主元分析法，对每个振动子序列进行特征提取，得到特征参数矩阵；步骤4，按照步骤1-3对已知故障类型的旋转机械提取特征参数矩阵，作为训练样本；以旋转机械的特征参数矩阵训练样本为输入、故障类型为输出，训练支持向量机模型，得到旋转机械故障诊断分类器；步骤5，实时采集待检测旋转机械的原始振动序列，按照步骤1-3获取待检测旋转机械的特征参数矩阵；将待检测旋转机械的特征参数矩阵输入到旋转机械故障诊断分类器，旋转机械故障诊断分类器输出待检测旋转机械的故障类型。2.根据权利要求1所述的方法，其特征在于，所述支持向量机模型的输入层包括1个节点，输出层包括1个节点，核函数为径向基函数，支持向量机模型的惩罚系数和核函数参数采用狼群-模拟退火混合算法进行优化选择。3.根据权利要求2所述的方法，其特征在于，所述支持向量机模型的惩罚系数和核函数参数采用狼群-模拟退火混合算法进行优化选择的具体过程为：步骤4.1，以个体狼位置作为支持向量机分类模型的惩罚系数、核函数参数，随机初始化狼群中的每只个体狼并设置狼群参数；狼群规模取值范围为：[50，200]，步长因子取值范围为：[50，120]，探狼比例因子取值范围为：[2，6]，最大游走次数取值范围为：[10，40]，距离判定因子取值范围为：[40，100]，最大奔袭次数取值范围为：[4，16]，更新比例因子取值范围为：[2，6]，最大迭代次数取值范围为：[500，1000]，最大搜索精度取值范围为：[0.001，0.005]；设定模拟退火算法的退火初始温度为100、退火速率为退火迭代次数t2＝1，当前温度下最大退火循环次数为Lmax＝6；步骤4.2，设定适应度函数，并确定初始最优头狼位置和迭代次数t，t＝1；将个体狼位置的对应的...

【专利技术属性】
技术研发人员：刘辉，段铸，李燕飞，施惠鹏，徐一楠，
申请(专利权)人：中南大学，
类型：发明
国别省市：湖南,43