网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

轮对轴承安全域增量建模方法

  • 刘志亮 1
  • 孙文君 1
  • 康金龙 1
  • 秦勇 2
1. 电子科技大学 机械与电气工程学院 成都 611731; 2. 北京交通大学 轨道交通控制与安全国家重点实验室 北京 100044

发布日期:2020-06-29

DOI:10.6052/1672-6553-2020-027

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
目录contents

摘要

安全域是一种从域的角度描述轮对轴承可安全运行区域的定性模型.然而,高速列车运行状态的样本数据通常是正常样本,故障样本的缺乏意味着只能使用正常样本进行安全域建模.支持向量数据描述是一种单分类方法,可以只利用一类样本进行建模.现有的滚动轴承安全域建模方法基本采用的是批量学习算法,而在样本数据不断增加的情况下,批量算法的效率低下.因此,本文将增量支持向量数据描述方法引入轮对轴承的故障诊断领域,提出了一种基于样本增量学习的高速列车轮对轴承安全域的建模方法.实验结果表明,相比于传统批量学习算法,本文提出的方法在轮对轴承安全域的建模中更加高效,并且分类准确率基本不受影响.

* 国家重点研发计划项目(2016YFB1200401-106),国家自然科学基金资助项目(6183302)

引言

高速列车是铁路运输中迅速发展起来的一种新型交通工具,具有快速性、稳定性、安全性及舒适性等优点.滚动轴承作为高速列车的关键部件之一,其工作性能的好坏直接影响列车的运行安全.然而由于滚动轴承在列车运行中起着承受和传递载荷的作用,因此,它是易损零件,故障率高.可见,准确高效的滚动轴承故障诊断对于高速列车安全运行具有重要意义.

安全域是一种从域的角度描述系统整体可安全稳定运行区域的定性模型.利用这种模型,可以通过比较安全域边界与系统运行点的相对关系分析得到系统的运行状[

1].最早在20世纪80年代,在电力系统安全领域,Wu[2]就提出了安全域分析及估计的相关理论,之后安全域理论扩展到信息安全以及轨道交通等领域.在轨道交通领域,秦勇[3]将安全域理论应用到列车关键设备的运行安全状态监测及评价中,提出了关于滚动轴承的安全域估计和状态辨识方法,该方法需要正常及故障样本进行安全域边界的估计.然而高速列车滚动轴承运行状态数据主要是正常的,并且获取故障样本的成本很高.因此需要研究仅使用正常运行状态的样本数据进行安全域建模的方法.

SVDD是Tax和Duin[

4]在支持向量[5](support vector machine,简称SVM)的基础上提出的一种单分类算法.与SVM的最优超平面不同,SVDD旨在求解一个尽可能包容所有目标样本数据的最优超球体.考虑到SVDD能仅利用一类样本进行建模,而高速列车轴承又缺乏故障样本,因此本文利用SVDD算法对正常样本进行学习获得轴承的安全域模型,从而在故障诊断时就能对故障数据进行异常检测.

SVDD发展至今,除了传统的批量学习算法外,还有不少学者对增量学习的SVDD算法进行了研究和改进.例如,刘健[

6]提出一种基于有效向量集的SVDD增量学习算法—EISVDD.段修生[7]基于准确挖掘出样本中所有的支持向量,提出一种新的超球SVDD增量学习算法.花小朋[8]通过对原样本和新增样本进行筛选,舍弃无用样本,提出一种改进的SVDD增量学习算法—IISVDD.冯国瑜[9]针对在线学习问题,总结出一种SVDD增量学习方法—ISVDD.

在故障分析与诊断领域中,蒋勉等[

10]针对转子-滑动轴承系统支承松动-碰摩故障动力学行为进行了分析,并提出基于动力学行为非线性度量的转子-轴承系统支承松动状态评估方法,韩清凯等[11]指出转子系统故障诊断的未来研究主要聚焦在基于大数据的智能诊断方面.而在轴承故障诊断领域,现在依然采用批量学习算法对高速列车滚动轴承进行安全域建模.在实际应用中发现,要在训练初期就收集到一个非常完整的训练集是很难实现的,更多情况下,训练样本是不断加入的.而批量学习算法在面对新加入的训练样本时,必须将新样本和原有样本组合重新进行训练,这个过程将占用大量的时间和内存,使得建模效率低下.因此用增量学习算法进行高速列车滚动轴承安全域的建模具有重要意义.

有鉴于此,本文将增量学习算法引入到轴承的故障诊断中,提出一种基于样本增量学习的高速列车滚动轴承安全域进行建模方法.该方法使用增量学习的SVDD进行轴承安全域建模,对于新加入的训练样本,可以不用舍弃之前的训练好的安全域模型,而是渐进地更新安全域边界,同时实现了仅使用正常样本建模及增量学习的目的.并且设计实验比较了批量学习的SVDD与本文方法在构造的高速列车轴承安全域模型上的差异,为高速列车滚动轴承安全域的建模方法提供了一种参考.

本文的组织结构如下.首先,第1节介绍了滚动轴承安全域的基本概念;第2节提出了基于样本增量学习的滚动轴承安全域建模方法;第3节通过实验验证本文提出方法的有效性,并与基于批量学习SVDD安全域建模方法进行比较;最后得出结论.

1 滚动轴承安全域的基本理论

在滚动轴承的故障诊断中,安全域是一种从域的角度描述滚动轴承可安全稳定运行区域的定性模型.安全域建模是指在轴承运行状态的状态特征变量所确定的空间内,通过分类方法求解得到一个安全域边界,该边界将空间分为两个区域,安全域与非安全域.一个二维空间下的安全域示意图如图1

图1 安全域概念示意图

Fig.1 Schematic diagram of safety region concept

示,其中xy表示两个轴承的两个状态特征变量,阴影部分表示安全域区域,阴影部分外的区域即表示非安全域区域.若轴承的状态特征变量落在安全域内,则表示其运行状态是正常的;若落在安全域以外的区域,则表示其处于故障的运行状态.因此可以用安全域辨识当前轴承的运行状态是否正常.

接下来本文将增量学习SVDD应用于高速列车轴承安全域模型,详细介绍基于样本增量学习的安全域模型建立过程以及如何将建立好的安全域模型应用到轴承的故障诊断,辨识轴承运行状态的安全性.

2 基于样本增量学习的安全域建模方法

如同名称所示,SVDD是一种描述数据域的专业方法.因此,它可以用来描述轴承正常运行状态所对应的安全域.利用轴承正常数据建立的SVDD模型便可等同于安全域模[

4,12,13].

基于样本增量学习SVDD的安全域建模方法介绍如下:

设训练样本集为{xi|xiRn, i =1,2,⋯,N},其中xi是由反映轴承运行状态的特征变量所构成的向量,n表示向量的维度,N是训练样本个数.SVDD旨在求解包含所有训练样本的最小超球体的球心a和半径R,其数学表达式为

min  R2+Ciξis. t. Φxi-a2R2+ξi, ξi>0, i=1,2,,N (1)

其中,Φ(xi)是xi经核函数在核空间中的映射,C是惩罚参数,用以平衡R与松弛变量ξi.

利用拉格朗日乘子法和Wolfe对偶定[

14]可以将等价为无约束二次凸优化问题.

minL=i,j=1nαiαjκ(xi,xj)-i=1nαiκ(xi,xi)s.t.0αiC,i=1nαi-1=0 (2)

其中,κ(xi, xj)=<Φ(xi), Φ(xj)>,是核函数.一般取高斯径向基函数作为核函数,其表达式为κ(xi,xj) = exp(-xi-xj2/(2σ2)).求解即可得到αi,并且超球体半径R与中心a可由αi表示,如所示.

a=iαiΦ(xi) (3)
R2=Φ(xsv)-a2=κ(xsv,xsv)-2iαiκ(xsv,xi)+i,jαiαjκ(xi,xj) (4)

其中,xsv是指0<α<C对应的支持向量.

在轴承故障诊断中,随着新样本的不断到来,可以通过渐进地更新αi的值实现增量学习,不需要重新训练SVDD模型.接下来将介绍在新样本加进来时,如何更新αi的值.

首先将转换为无约束凸规划问题,即:

minL=i,j=1nαiαjκ(xi,xj)-i=1nαiκ(xi,xi)+            b(i=1nαi-1) (5)

的最优解需要满足Kuhn-Tucker(KT)条件,即:

gi=Lαi=2j=1nαjκ(xi,xj)-κ(xi,xi)+b            0,αi=0,(xiR)=0,0αiC,(xiS)0,αi=C,(xiE) (6)
Lb=i=1nαi-1 (7)

其中,giLαi的偏导,表示使L取到最小值的αi*值有以下三种情况,如图2所示.

图2 αi取值的三种情况

Fig.2 Three cases of αi values

(1)若αi的最优值αi*=0,即区间[0,C]的左端点,αi*R,则L在该区间单调递增,如图2 (a);

(2)若αi的最优值0<αi*<C,即区间[0,C]内,则L在该区间先减后增,如图2 (b);

(3)若αi的最优值αi*=C,即区间[0,C]内,则L在该区间单调递减,如图2 (c).

若通过已有样本求出最优解为α=(α1,α2,⋯,αN),当有一个新样本xc加入时,设新的最优解为α*=(α1*,α2*,⋯,αN*,αc),也必须满足KT条件,通过调整αc,并在满足KT条件下相应地调整其他αi值,使αc也满足KT条件,最后得到满足KT条件的最优解α*.

根据这个原理,假设∆αc足够小以致于不引起任意一集合RSE中的点变为另一集合中的点.由可知集合RE中的点所对应的αi分别固定为0和C,其值不能改变,否则将不满足KT条件,因此可以改变值的αi只能对应集合S中的点,并且此类αi对应的gi必须固定为0,即其增量∆gi也应为0.将的等式两端写成增量形式,得到下面两式:

Δgi=2jRSEΔαjκ(xi,xj)+2Δαcκ(xi,xc)        +Δb,iRSE (8)
iRSEΔαi+Δαc=0 (9)

除去不可变的αi(αiRE),只剩下αiS,且对应gi=0,则可以简化为:

0=2j=s1smΔαjκ(xi,xj)+2Δαcκ(xi,xc)+Δb,iS (10)
0=i=s1smΔαi+Δαc (11)

其中,s1为集合S的第一个元素的下标,sm为集合S的最后一个元素的下标.

结合,得到如下矩阵形式:

01112κs1s12κs1sm12κsms12κsmsmΔbΔαs1Δαsm+12κs1c2κsmcΔαc=0000 (12)

其中,κs1sm表示κ(xs1,xsm),其他以此类推.记矩阵Q为:

Q=01112κs1s12κs1sm12κsms12κsmsm (13)

做简单移项,并令矩阵U=Q-1可得:

ΔbΔαs1Δαsm=-U12κs1c2κsmcΔαc (14)

记:

ββs1βsm=-U12κs1c2κsmc (15)

联立得:

Δb=βΔαc (16)
Δαi=βiΔαc,iS (17)

对于iRE,设βi=0,则对所有的iRSE均满足如下关系:

Δαi=βiΔαc,iRSE (18)

ββi都称为系数敏感[

12](coefficient sensitivities),将带入,可得

Δgi=2κic+2jSκijβij+βΔαc,iREc0,iS (19)

于是,总结出增量学习方法为:对于新训练样本的xc,设其对应的αc为0,检验αcgc的关系是否满足KT条件(即).若满足,则原始的αi不变,αc取0;若不满足,则给αc一个增量∆αc,并利用给出的∆αc,∆αi和∆gi之间的关系,求出∆αi和∆gi,并更新αcαigi,再次检验新的αcgc是否满足KT条件,若不满足,则反复迭代此增量过程直至满足.

上述结论基于这个假设:∆αc足够小以致于不引起任一集合RSE中的点变为另一集合中的点.但在寻找最优值的过程中,一次迭代必然要引起点从一个集合变到另一集合,那么就需要相应地更新集合RSE以及βi的值,以新的βi去计算增量∆αi和∆gi.因此,每次迭代需校对如下条件以实现集合中样本点的更新.

(1) gc≤0,若gc取0,则xc进入集合S

(2) αcC,若αcC,则xc进入集合E

(3) 0<αi<C,i∈S,若αi取0,则xi从集合S转移R;若αiC,则xi从集合S转移E

(4) gi0iE,若gi取0,则xi从集合E转移S

(5) gi≥0,iR,若gi取0,则xi从集合R转移S.

其算法流程图如图3所示:

图3 ISVDD算法流程图

Fig.3 ISVDD algorithm flow chart

在新样本加入的情况下,可以通过上面的方法增量学习更新SVDD模型的参数,进而得到最小超球体的球心a和半径R,从而建立轴承的安全域模型.

设轴承状态特征变量构成的测试样本z,可以用这个安全域模型判断轴承的运行状态,其决策函数如下:

f(z)=(Φ(z)-a)T(Φ(z)-a)=κ(z,z)-      2i=1nαiκ(z,xi)+i,j=1nαiαjκ(xi,xj) (20)

如果f(z) ≤ R2,表示z处于安全域内,即此时轴承处于安全状态;如果f(z) > R2,表示z处于安全域外,即此时轴承处于故障状态.

3 本文方法在实验数据上的验证

为了验证本文方法的有效性,本小节采用高速列车轮对轴承的实验数据.该实验数据通过青岛四方所的轮对轴承实验台采集,如图4所示,该实验台具有75 kW驱动电机、传动比为1.4的带传动系统、最大载荷为250 kN的垂直载荷加载台、最大载荷为100 kN的横向载荷加载台、最大转速为1480的两个风扇马达和控制系统.垂直和横向载荷加载台用来模拟实际列车运行中的二维载荷,风扇电机可产生与列车运行方向相反的风,实验轴承安装在轴箱中.对于振动数据采集系统,加速度传感器安装在轴箱箱体的十二点位置,采样频率为5120 Hz.

图4 轮对轴承实验台

Fig.4 Wheelset bearing test bench

实验采集了一种工况下,滚动轴承正常状态和内圈故障、外圈故障及滚动体故障三种故障状态下的振动信号.然后,将每个状态下的原始振动信号进行分割,得到50个信号片段,则总计共有50个正常样本和150个故障样本(内圈故障、外圈故障和滚动体故障各50个).当旋转机械发生故障时,其状态监测信号会有不同的响应,典型的检测方法是计算振动信号在时域或者频域中的特征.Antoni[

15]受热力学领域的启发,验证了振动信号的平方包络谱的负熵(negentropy of squared envelope spectrum,简称NESES)和平方包络的负熵(negentropy of squared envelope,简称NESE)是反映滚动轴承故障的有效特征值.因此,本文即选择这两个状态量作为增量安全域模型的输入变量进行实验.该实验使用Windows 10操作系统,主频为4.2GHz的i7-7700K处理器,内存为16GB.

3.1 批量算法和增量算法在滚动轴承安全域建模中的差异

为了验证本文方法的可行性,分别用批量学习的SVDD算法和本文方法对上述轴承正常数据进行安全域建模.惩罚参数根据C=1/(Nv)公式计算,其中v为误检率,本文固定为0.1.选用高斯径向基函数,参数σ = 0.2.分别用30个正常样本和50个正常样本得到的安全域模型如图5图6所示.

图5 使用30个样本训练得到的安全域模型

Fig.5 Safety region model with 30 training samples

图6 使用50个样本训练得到的安全域模型

Fig.6 Safety region models with 50 training samples

通过对比图5图6所示的安全域边界,可以看出本文方法建立的安全域模型与基于批量学习SVDD建立的安全域基本无差异,因此可以认为两种算法得到的优化结果是一致的,说明本文方法建立安全域模型是有效的.

上面的实验结果验证了本文基于样本增量学习的安全域建模方法是可行的.接下来,为了进一步验证本文方法的高效性与有效性,设计实验如下:将正常样本随机分为5份,其中4份作为训练样本,另外1份再加上150个故障样本作为测试样本.从30个样本开始依次增加训练集的样本数量,用批量算法和本文方法分别进行训练,并用测试集进行测试,记录安全域建模时间和分类正确率.其中,批量安全域建模时间为将样本批量导入之后进行训练最终得到模型的时间,增量安全域建模时间为新样本加入前已有样本建立模型所需时间与新样本加入时样本集合更新所需时间之和.

将五折交叉验证实验重复进行100次,得到实验结果如图7和8所示. 实验结果表明:随着样本的不断加入,本文方法的安全域建模时间曲线上升缓慢,而批量学习的SVDD算法的建模时间曲线上升很迅速,这是因为批量学习方法在新样本加入的时候,需要舍弃之前的模型,重新训练安全域模型,而本文方法可以增量地更新安全域模型,因此本文方法的效率比批量算法的效率高;本文方法的分类正确率曲线与批量学习的SVDD算法的分类正确率相差不多,其平均相差不到1.5%,并且在此实验条件下与样本量并无明显关联,综合以上分析可以说明本文方法在故障识别上是有效的.需要特别说明的是,在本次实验中,分类准确率和样本数量并无特定的相关关系.

图7 模型建立时间对比

Fig.7 Model establishment time comparison

图8 分类正确率对比

Fig.8 Comparison of classification accuracy

4 结论

本文针对高速列车滚动轴承故障样本少和批量学习算法效率低问题,提出了一种基于样本增量学习的SVDD安全域建模方法.在高速列车轴承实验数据上的验证结果表明,本文方法在保证准确率的同时能极大地缩短模型训练时间.

参 考 文 献

1

张媛. 基于安全域的列车关键设备服役状态辨识与预测方法研究[博士学位论文].北京:北京交通大学,2014 (Zhang Y.Service status identification and prediction method research based on safety region estimation for the key experiments in rail vehicles[Ph.D Thesis].Beijing: Beijing Jiaotong University, 2014 (in Chinese))

2

Wu F, Kumagai S. Steady-state security regions of power systems. IEEE Transactions on Circuits & Systems, 1982, 29(11): 703~711

3

张媛, 秦勇, 邢宗义,.基于LMD-PCA-LSSVM的滚动轴承安全域估计和状态辨识方法.振动与冲击,2013, 32(20):172~178

Zhang Y, Qin Y, Xing Z Y, et al. Safety region estimation and state identification of rolling bearings based on LMD-PCA-LSSVM.Journal of Vibration and Shock, 2013, 32(20):172~178 (in Chinese)

4

Tax D M J, Duin R P W. Support vector data description. Machine Learning, 2004, 54(1): 45~66

5

Vapnik V N.The nature of statistical learning theory. New York: Springer Verlag, 1999

6

刘建, 万坚,黄渊凌.一种基于有效向量集的SVDD增量学习算法.电信技术研究, 2013215~20 (Liu J, Wan J, Huang Y L. An incremental learning algorithm for SVDD based on effective vector set. Research on Telecommunication Technology,2013,2:15~20(in chinese))

7

段修生, 曹健, 孙世宇, .一种新的超球SVDD增量学习方法. 计算机应用与软件, 2015, 32(2):163~166

Duan X S, Cao J, Sun S Y, et al. A new incremental learning method for hypersphere SVDD. Computer Applications and Software, 2015,32(2):163~166 (in Chinese)

8

花小朋, 皋军, 田明, .改进的SVDD增量学习算法. 计算机工程, 2009, 35(22):210~211

Hua X P, Gao J, Tian M, et al. Improved incremental learning algorithm for support vector data description. Computer Engineering, 2009, 35(22):210~211 (in Chinese)

9

冯国瑜, 肖怀铁, 付强, .一种适于在线学习的增量支持向量数据描述方法.信号处理,2012, 28(2):186~192

Feng G Y, Xiao H T, Fu Q, et al. An incremental support vector data description method for online learning. Signal Processing, 2012, 28(2): 186~192 (in Chinese)

10

蒋勉, 伍济钢, 彭鑫胜, .转子-滑动轴承系统支承松动-碰摩故障动力学行为及评估方法.动力学与控制学报, 2017, 15(6): 550~557

Jiang M, Wu J G, Peng X S, et al. Dynamic behaviors and assessment of loose-rubbing faults in rotor-sliding bearing system.Journal of Dynamics and Control, 2017, 15(6): 550~557 (in Chinese)

11

韩清凯, 马辉.转子动力学研究进展.动力学与控制学报, 2018,16(6): 481~482

Han Q K, Ma H. The research progress of rotor dynamics.Journal of Dynamics and Control, 2018, 16(6): 481~482 (in Chinese)

12

Cauwenberghs G, Poggio T.Incremental and decremental support vector machine learning// Dietterich T G, Becker S, Ghahramani Z.Advances in neural information processing systems: Cambridge MIT, London MIT, 2001: 409~415

13

Tax D M J, Laskov P. Online SVM learning: from classification to data description and back. In:IEEE XIII Workshop on Neural Networks for Signal Processing, 2003: 499~508

14

邓乃扬.数据挖掘中的新方法-支持向量机.北京:科学出版社,2004(Deng N Y. A new method in data mining-support vector machine.Beijing: Science Press, 2004 (in Chinese))

15

Antoni J.The infogram:Entropic evidence of the signature of repetitive transients. Mechanical Systems & Signal Processing, 2016, 74:73~94

微信公众号二维码

手机版网站二维码