当滚动轴承的振动信号包含强干扰噪声,振动信号的频谱划分严重噪音干扰。传统的经验小波变换(易)将信号分解成大量的组件,并且很难选择合适的组件,包含故障信息。为了解决上述问题,在本文中,我们提出了改进的经验小波变换(IEWT)方法。仿真实验证明,IEWT可以解决大量的问题易组件和单独的组件的影响有效地从噪声包含轴承故障信息。IEWT方法结合支持向量机(SVM)诊断滚动轴承的故障。排列熵(PE)是用于构造特征向量的感应能力强动态变化的非平稳、非线性信号。惩罚因子的关键参数<我t一个l我c>
C我t一个l我c>和内核参数<我nl我ne-formula>
滚动轴承是旋转机械的重要组成部分。其工作环境恶劣,需要复杂的负载。它是旋转机械的重要组成部分,容易断裂。研究状态监测和故障诊断是一个必要的措施确保旋转机械的正常运行。
滚动轴承的振动信号是非线性和非平稳的。最常用的方法来处理这些信号经验模态分解法(EMD)提出( 对EMD缺乏完整的理论推导,mode-mixing等问题和端点效应,实证小波变换(易),提出了基于小波变换的 习( 轴承信号采集的过程中,它是很难获得大量的轴承数据和大量的数据也很难分析和过程。SVM是一个智能故障诊断方法基于结构风险最小化的原则,突出优点为小样本数据( 本文提出了一种基于MI谱重新分区方法来解决这个问题,易被噪声严重影响并产生许多谱划分点,导致太多的分解组件。成功提取出故障影响组件从一个复杂的时域系列,使信号分解的准确性。
本文提出了一种新的方法来解决这个问题所引起的过度易组件。和一个基于IEWT和支持向量机的滚动轴承故障诊断方法。应用组件的PEs作为SVM的输入特征向量,以及惩罚因子<我t一个l我c>
C我t一个l我c>和内核参数<我nl我ne-formula>
易是一种自适应小波变换理论的框架下分解方法。首先,振动信号频谱分为自适应,然后,使用的正交小波滤波器组分解获得紧支持特色调幅-调频信号( 计算振动信号的频谱、频谱定义的范围内(0<我t一个l我c>
π我t一个l我c>),频谱分为<我t一个l我c>
N我t一个l我c>段,<我t一个l我c>
ω我t一个l我c>0年代ub>= 0和<我t一个l我c>
ω我t一个l我c>
根据Meyer和Littlewood-Paley方程,实证构造小波。相应的尺度函数和小波函数定义如下: 实证根据经典的小波变换的小波构造方法。细节系数和近似系数描述如下: 根据小波理论,经验模态函数得到以下方程: 根据方程(信号重建 根据方程(
支持向量机是一种机器学习方法基于VC维理论和经验风险最小化原则。支持向量机在分类和回归显示了强大的优势的小样本数据( 考虑到数据{<我t一个l我c>
x我t一个l我c>我年代ub>,<我t一个l我c>
y我t一个l我c>我年代ub>;<我t一个l我c>
我我t一个l我c>= 1,2,…<我t一个l我c>
n我t一个l我c>},<我t一个l我c>
x我t一个l我c>我年代ub>样本数据和吗<我t一个l我c>
y我t一个l我c>
当松弛因子<我t一个l我c>
ξ我t一个l我c>和惩罚因子<我t一个l我c>
C我t一个l我c>引入解决方案过程中,优化问题表示为 内核函数<我t一个l我c>
K我t一个l我c>(<我t一个l我c>
x我t一个l我c>
本文采用RBF核函数和方程
支持向量机分类过程中,惩罚因子<我t一个l我c>
C我t一个l我c>和内核参数<我nl我ne-formula>
尽管遗传算法的结构和算法简单,容易实现,和快速收敛,很容易落入过早成熟,很难选择最优参数( 相结合的关键是构造适当的量子门和更新量子比特量子门旋转。具体操作如下所示: 假设有一个包含人口<我t一个l我c>
N我t一个l我c>个人,<我nl我ne-formula>
量子位编码相结合的优化是实现量子旋转门和更新。建立一个新的人口、初始化种群和测量量子位二进制值都是由量子位编码完成。为了进一步优化模型参数,计算人口健身价值和最优解保存,通过量子旋转搜索和更新当前量子位,毛圈和输出最优解来满足终止条件(
排列熵是衡量时间序列的复杂性,具有良好的能力分析的复杂性和意义的非平稳、非线性信号的微小变化和良好的抗噪音。它显示了巨大的优势代表了旋转机械故障振动信号的状态( 定义一个时间序列{<我t一个l我c>
X我t一个l我c>(<我t一个l我c>
我我t一个l我c>),<我t一个l我c>
我我t一个l我c>= 1,2,…<我t一个l我c>
k我t一个l我c>},根据香农熵。然后,排列熵可以被定义为 的范围<我t一个l我c>
H我t一个l我c>
嵌入维度<我t一个l我c>
米我t一个l我c>和时间延迟<我t一个l我c>
τ我t一个l我c>有一个很大的影响排列熵的计算结果。在[
在尺度空间平面上,每个最小值对应于一个尺度空间曲线。Ostu阈值确定Th的方法,和截止点的阈值大于保留。然而,Ostu法获得的阈值相对较小,导致太多的截止点选择,导致过度易组件( 互信息(MI)是由信息论中熵的概念。总是用来衡量两个随机变量和的关联度比相关系数(更准确 易解决过度组件的问题,提出一种改进的易的方法。的最大阈值忽略多余的截止点,只有有效的频谱截止点依然存在。根据原始信号的相关性,进一步结合冗余组件。振动信号最终分解成有意义的组件。该方法包括四个步骤,写如下: 步骤1:FFT进行振动信号获取频谱。
步骤2:Ostu法用于确定阈值尺度空间曲线的Th,和尺度空间方法用于确定初始边界点获取初始组件。
步骤3:计算初始的MI组件,并将根据MI组件。如果组件的MI大于均值和相邻的MI组件大于均值,这两个组件合并。如果相邻的MI组件小于均值,组件是独立的。如果这个组件的MI小于均值和相邻的MI组件小于平均值,那么这两个组件合并。如果相邻的MI组件大于均值,保持独立的组件。根据组件的组合,选择在初始边界点获得保留划分点。
第四步:根据保留边界点,信号的傅里叶谱改变。
第五步:傅里叶谱的基础上,构造一个新的正交滤波器组和一种新的分解组件。IEWT如图的实现过程
模拟信号是用方程(
步骤1后,FFT进行模拟信号获取频谱。
步骤2后,阈值是通过如图所示的Ostu法 获得边界点和分裂。(一)尺度空间方法获得初始边界点。(b)划分频谱。
IEWT把频谱。(一)初始组件的MI。(b)的重新分配频谱。
步骤3后,只保留3个边界点,如图 步骤4后,频谱分为4段根据保留截止点如图 步骤4后,IEWT组件,如图 时域图IEWT组件。
从IEWT组件如图 相关系数。
峰度IEWT组件。
根据每个组件之间的相关系数模拟信号和IEWT组件表 为了验证IEWT方法的分解效果,EEMD方法被用于比较。图 EEMD组件。
从图可以看出 EEMD组件的相关系数。
峰度EEMD组件。
根据表 根据IEWT和EEMD组件的数量,优化组件峰度值之间的相关系数组成组件,和相应的组件在上面的表中,分解两种方法综合比较的影响。
根据表 比较IEWT EEMD。
总之,模拟信号的对比分析证明了IEWT解决传统易产生的缺陷冗余组件,大大减少了组件的数量,并提供方便后续的分析和处理。EEMD通过比较,证明比EEMD IEWT有更好的分解作用,不仅将每个组件分解为单个组件也将更好的故障影响的组件,同时,有一个更快的分解效率。
组件
IEWT1
IEWT2
IEWT3
IEWT4
8.2056<我t一个l我c>
e我t一个l我c>−04
0.1855
0.1445
−2.7467<我t一个l我c>
e我t一个l我c>−05
8.3583<我t一个l我c>
e我t一个l我c>−07
−5.3453<我t一个l我c>
e我t一个l我c>−06
0.0125
0.4286
0.2312
组件
IEWT1
IEWT2
IEWT3
IEWT4
峰度
1.5891
4.3010
3.1950
组件
IMF1
IMF2
IMF3
IMF4
IMF5
IMF6
IMF7
IMF8
IMF9
IMF10
IMF11
IMF12
IMF13
0.4385
0.4337
0.0869
0.0614
0.0431
0.0372
0.0051
0.0010
0.0036
0.0002
0.0000
0.0029
0.5654
0.0029
0.0060
0.0023
0.0026
0.0103
0.0055
0.8514
0.0279
0.0171
0.0217
0.0661
0.6487
0.3905
0.3491
0.2510
0.1743
0.1064
0.0009
0.0189
0.0059
0.0083
0.0040
0.0008
IMF1
IMF2
IMF3
IMF4
IMF5
IMF6
IMF7
IMF8
IMF9
IMF10
IMF11
IMF12
IMF13
3.6889
4.6746
3.3461
3.4055
3.0461
3.0961
3.1959
1.5658
2.1009
1.5896
1.4678
1.8313
方法
相关系数
峰度的最佳组件
数量的组件
运行时间/(年代)
IEWT
0.7501
8.2330
4
4.285
0.9928
0.7710
EEMD
0.4565
4.7932
13
20.074
0.9799
0.6636
IEWT方法提出了改进传统的易。的PEs IEWT组件作为SVM的输入特征向量计算。相结合用于支持向量机的参数优化模型,用于预测模型和最优参数。写具体实现步骤如下:
步骤1:IEWT方法用于不同状态的振动信号分解和保留前5组件。
步骤2:计算所有组件的PEs,规范化,作为支持向量机的特征向量输入。
步骤3:选择10组样本在每个州训练模型。采用相结合的优化算法。适应度函数的预测精度作为交叉验证(简历)。最大进化一代是200,人口是20,的范围<我t一个l我c>
C我t一个l我c>(0.1,100),的范围<我nl我ne-formula>
第四步:输入到训练SVM模型最优参数,并选择15组样本在每个州的错分类预测。图
故障诊断流程图。
为了验证本文提出的方法的有效性,轴承的数据采用凯斯西储大学的实验室实验( 轴承参数。
一群随机选择的时域图绘制在图状态信号 每个状态信号的时域图。
特征向量。
它可以观察到从表 十组训练样本在不同状态的特征向量输入到支持向量机,分别实现采用交叉验证与分类精度适应度函数找到最优的解决方案。图 实现优化的结果。
根据实现迭代优化结果如图 IEWT-SVM的分类结果。
IEWT-SVM的分类结果表明,该方法可以有效地、准确地诊断轴承故障。为了验证IEWT-SVM的优势,传统的易和VMD方法是用来比较。传统的易方法的分类结果图所示 EWT-SVM的分类结果。
当实现优化参数<我nl我ne-formula>
VMD-SVM的分类结果。
组件通过VMD也发送到支持向量机分类。当相结合优化结果<我nl我ne-formula>
为了证明SVM的分类预测效果,不同的分类器被用来分类特征向量。再(资讯)算法是一种基于先验知识的监督学习算法,具有结构简单和强大的分类性能,因此被广泛应用于故障诊断领域( 预测精度在不同<我t一个l我c>
K我t一个l我c>值。
从图可以看出 IEWT-KNN的分类结果<我t一个l我c>
K我t一个l我c>= 1。
极端学习机(ELM)是一种监督学习算法基于一个单隐层前馈神经网络(SLFN)。与传统的梯度下降算法相比,榆树有较短的训练时间和更好的泛化性能,被广泛用于机械故障诊断领域的( IEWT-ELM的分类结果。
故障诊断结果的比较。
通过比较在表5的结果的方法
内座圈直径(毫米)
外环直径(毫米)
节圆直径(毫米)
滚动体直径(毫米)
球的数量
接触角(°)
25
52
39.04
7.94
9
90年
轴承状态
样品标签
排列熵
E1
E2
E3
E4
E5
正常的
1
1.163
1.988
2.351
2.844
2.676
1.205
1.519
1.957
2.119
2.368
1.190
1.523
1.960
2.125
2.358
外环的错1
2
3.092
3.229
3.493
2.888
2.092
3.109
3.244
3.409
2.991
2.414
3.144
3.188
3.468
2.966
2.069
外环故障2
3
1.245
2.018
3.877
3.549
3.895
1.140
1.784
3.952
3.366
3.928
1.158
1.984
3.950
3.589
3.894
内套错1
4
1.205
1.519
1.957
2.119
2.368
1.190
1.524
1.960
2.125
2.358
1.199
1.492
1.961
2.133
2.346
内套错2
5
3.008
3.474
2.859
3.410
3.259
3.012
3.497
2.764
3.678
2.892
2.982
3.571
2.855
3.390
3.299
滚动体故障1
6
1.526
2.197
4.148
3.148
3.165
1.486
1.995
4.118
3.171
3.111
1.484
2.020
4.114
3.239
2.982
滚动体故障2
7
1.140
2.164
4.406
2.802
3.088
1.158
2.120
4.323
2.880
2.972
1.161
2.143
4.325
2.835
2.972
方法
分类错误的样本
Correct-classified样品
样本总数
准确率(%)
运行时间(s)
IEWT-SVM
3
102年
105年
97.14
129.29
EWT-SVM
24
81年
105年
77.14
126.40
VMD-SVM
9
96年
105年
91.43
379.10
IEWT-KNN
9
96年
105年
91.43
121.11
IEWT-ELM
11
94年
105年
89.52
111.04
振动信号分析是一种常见的和可用的方法来诊断轴承故障。易有很大的优势在信号分解mode-mixing,跑得很快。然而,会有过度易导致组件组件分解的选择困难。要克服的问题,易产生许多冗余组件,IEWT方法提出。模拟信号分析验证IEWT方法不仅可以有效地减少组件的数量也从复杂的信号中提取故障影响组件。摘要IEWT方法结合支持向量机和PEs IEWT组件构造的特征向量。相结合是用来优化惩罚因子<我t一个l我c>
C我t一个l我c>和内核参数<我nl我ne-formula>
实验数据支持分析由凯斯西储大学,和它的免费下载的网站已经被引用。
作者宣称没有利益冲突有关的出版。
这项工作是由中国国家自然科学基金(11790282;11572206;U1534204;11802184),河北省自然科学基金(A2016210099),河北省人才项目的科研项目训练(A2016002036)和项目在河北大学的高级人才(GCC2014021)。