深圳市地铁集团有限公司运营总部
摘要:地铁车辆走行部作为地铁车辆的重要组成部分，在平常的故障检测中，如车轮踏面的擦伤及剥离等擦伤现象，会影响到车辆轴承损伤、轴温升高及钢轨波磨等问题出现。因此，如何准确识别车辆故障，为维护工作提高效率及质量，是每位地铁工作人员关注的问题之一。本文主要对地铁列车走行部的故障诊断工作中，应用模糊熵得到的准确识别进行分析。
关键词:地铁车辆；轨道振动；模糊熵；故障诊断
前言
地铁车轮的擦伤及剥离问题，对车辆运行产生大量的安全隐患。当前国内外常见的非接触在经的检测技术有：振动分析法、位移检测法、图像检测法、电信号检测法等。而实际应用中，振动分析法技术是较成熟且成本低的一种在线检测方法。在采取振动分析法时，其振动信号在实际应用中容易受到噪声干扰，很难有效提取故障特征信息。而集合经验模态分解(EEMD)通过在原始信号中加入白噪声序列辅助分析，解决了传统EMD存在的模态混叠问题。熵可用于定量描述信号的不确定性和复杂度统计特性，抗噪能力强，稳定性好。由EEMD分解得到若干平稳本征模态函数(IⅧ)并计算相应熵值，可作为车辆平轮故障信息的特征向量。特征提取后需要进行故障识别。支持向量机(SVM)建立在结构风险最小化原则和VC维概念基础上，在解决非线性、较高维、小样本等方面有突出优点，具有很好的泛化能力。
试验采集了地铁车辆在正常情况、擦伤故障、剥离故障等3种运行工况下的轨道振动信号。首先，在轨道两侧各对称安装5个通道传感器来采集列车经过时完整的车轮振动信号，然后通过轨道振动分析实现平轮故障级别的准确判定。
一、特征提取算法理论
（一）EEMD参数设置和主分量筛选
步骤1：设置高斯白噪声的幅值标准差比值系数k和执行EMD的总次数s。则信噪比为：
(1)
由于EEMD分解会产生包含故障信息相对较小的IMF伪分量，若不加筛选直接进行后续故障特征提取，不但增加计算量，而且影响特征向量对故障的识别。本文先采用皮尔逊相关系数法，计算各IMF分量与原始信号的相关系数；然后对相关系数设定阀值，筛选出能保留原信号主要信息的IMF分量，从而有效消除噪声影响。
（二）主分量模糊熵特征提取
近年来，样本熵和模糊熵等熵的概念被应用到机械故障诊断领域。由于熵值具有能反映信号复杂度的特点，故将熵值作为故障信息的特征参数。模糊熵是样本熵的改进算法。二者都是用熵值来表示信号序列的复杂程度，不同之处在于：模糊熵将模糊集合理论引入序列复杂度的计算，利用指数函数将相似性度量模糊化，用模糊隶属度函数替代硬阈值判据，使得模糊熵的值能够随参数稳定变化。与样本熵相比，模糊熵对重构相空间的维数和相似容限度等参数的依赖性更低。因此，本文计算各IMF的模糊熵值来构成特征向量以用于故障识别。模糊熵的计算过程如下：
（5）
二、故障识别算法理论
SVM主要用于模式分类和非线性回归。其主要设计思想是通过核函数的非线性变换把低维数据映射到高维特征空间，并在高维空间中寻找最优分类超平面，将训练样本正确分类并使分类间隔最大。寻找最优超平面的问题等同于解凸二次规划优化问题。优化条件是使两类样本之间的距离最小，根据拉格朗日方程和Karush—Kuhn—Tuker条件，求解得到最优分类函数为
常见的核函数主要有高斯核函数(也称径向基函数(RBF))、Sigmoid线性核函数、多项式核函数等。RBF核函数能将样本非线性地映射到高维空间，平滑性好，分类准确率高，故本文选择RBF核函数。
核函数参数g以及惩罚因子c对SVM的学习能力和泛化能力有着重要影响。相比传统大范围枚举法，启发式遗传算法可不必遍历所有参数点，能快速得到全局最优解。遗传算法(GA)是模拟自然界遗传机制和生物进化论而形成的一种全局并行随机搜索最优化方法，适用于大规模并行寻优计算。
三、组合模型故障诊断流程
地铁车辆轨道振动信号是典型的非线性非平稳信号。平轮故障诊断分为特征提取及故障识别2部分。具体诊断模型的振动信号处理流程如图1所示。
图1 振
图2 信号采集方案
车辆在长期运行和紧急制动过程中，车轮踏面会发生擦伤或剥离进而形成平轮。平轮深度为1mm以下时定义为擦伤故障，平轮深度为1mm以上时定义为剥离故障。轨道振动除了受平轮冲击载荷影响外，受车速和载重干扰影响也较大。因此试验设定在地铁工程车空载且车速为40km／h时，采集轨道两侧共10个通道传感器的振动信号。采样频率为7992 Hz，每个通道取1024个数据长度。如图2所示，磁钢距离第1个传感器0.7m，可通过车速计算该间隔的过车时间，通过采样频率来计算该间隔采样点数，并截断第1个传感器在列车过该间隔期间的采样数据。同理，可由2个传感器间隔为1.4 m来进行每个传感器起始有效数据的对准。采集正常工况、擦伤故障及剥离故障3种工况下的轨道振动信号。每种工况采集32组数据样本，并随机选取其中20组作为训练样本，将其余12组为测试样本。3种状态下训练样本共60组，测试样本共36组。
（一）集合经验模态分解和主分量筛选
地铁车辆平轮故障会引起走行部各部件不同频段的固有振动。振动信号采用EEMD分解，即根据信号自身特点自适应地将不同频段内的固有振动分解到不同IMF分量中。EEMD分解次数s取值70，k取值0.05，剥离故障振动信号分解产生的9个IMF分量如图3所示。IMF分量个数与信号本身特点有关。将IMF分量按频率由高到低排列，则振幅依次减弱。由图3可见，IMF分量基本平稳且没有模态混叠现象。
图3 剥离故障振动信号的原始信号EEMD分解结果
根据分解后的结果，计算原始振动信号和各IMF分量的皮尔逊相关系数，并采用阈值法剔除与原始信号相关性较小的伪分量。经综合考虑，筛选每个样本分解产生的5个相关系数较大的IMF分量以进行后续的特征提取。
（二）模糊熵特征提取
由于模糊熵可分析复杂信号中的确定性成分和随机成分，因此采用模糊熵构造平轮故障特征向量。计算模糊熵时，重构相空间维数m的值越大越能体现信号动态演化过程；而相似容限度r过大会加剧信息的丢失，过小则会增加噪声敏感性并导致熵值不确定性增加。m取值2，r取值0.3。计算每个样本EEMD分解产生的5个主IMF分量的模糊熵值。可视化每个样本点的熵值变化，取前3个主IMF的模糊熵值代表一个样本点，一共96个样本点，可视化结果如图4。由图4可知，模糊熵特征向量代表的样本点拥有良好的类内聚集性，类间边界清晰，能够提高后续故障识别的准确率。为对比模糊熵相比样本熵熵值算法的优越性，本文同时计算样本熵熵值并绘制三维可视图(见图5)。通过图4、图5两图比较表明，模糊熵构成的特征向量比样本熵有更好的分类效果。
图5 不同工况下每个样本取前3个主IMF的样本熵熵值
（三） GA优化SVM故障识别
本文共取96组样本，将其分为训练集和测试集，提取3种工况下的故障特征向量，将训练集特征向量输入到SVM。
根据上述理论，采用GA的SVM参数进行全局寻优。GA设定种群规模为20，最大进化代数为100，利用交叉验证的准确率作为适应度值。在遗传代数为50时取得最优值停止迭代，获得参数优化结果为g=1.701 7，c=0.86527，对训练集的分类准确率为98．333 3％。GA寻优如图6所示。
由表1可见，对轨道原始振动信号进行EEMD分解比EMD分解更有利于后续的特征提取，分类准确率也明显提高。熵能够表征信号的复杂程度，用模糊熵作为特征向量比样本熵有更高的分类准确率。对于小样本情况下的故障分类，用遗传算法优化的SVM比BP神经网络有更高的分类准确率。
五、结束语
针对地铁车辆轨道振动信号复杂、非线性、信噪比低及多源耦合等现象，本文利用EEMD对原始信号进行分解，得到一系列平稳本征模态函数。采用相关系数法筛选主分量，引入熵理论构造模糊熵故障特征向量，通过s VM建立平轮故障识别模型，利用遗传算法优化分类器的核函数和惩罚因子。试验结果表明，本文所提出的特征提取方法和故障识别模型能够有效判别地铁车辆车轮踏面擦伤和剥离等平轮故障形式，识别准确率高。
参考文献:
[1]颜秋，刘永明．基于matlab／simulink的车辆建模与故障分析EJ3．华东交通大学学报，2012．
[2]肖成勇，石博强，冯志鹏．基于EEMD和进化支持向量机的齿轮混合智能诊断方法研究[J]．机械科学与技术，2015(1)：86—89．
论文作者:罗建强
论文发表刊物:《防护工程》2017年第29期
论文发表时间:2018/2/28
标签：样本论文; 故障论文; 信号论文; 模糊论文; 特征论文; 分量论文; 分解论文; 《防护工程》2017年第29期论文;