基于IMF能量比的高速铁路钢轨波磨声学诊断方法技术

本发明专利技术涉及一种基于IMF能量比的高速铁路钢轨波磨声学诊断方法，属于高速铁路振动噪声技术领域，其步骤如下：(1)钢轨粗糙度测试；(2)集合经验模态分解；(3)本征模态函数IMF能量比：根据故障特征频率对应IMF信号的能量比进行钢轨波磨的故障识别，筛选得到钢轨波磨对应的IMF分量，通过HHT变换，得到Hilbert边际谱和瞬时频率。基于IMF能量比的高速铁路钢轨波磨声学诊断方法。本发明专利技术在运营线路上采用直接法实测有无波磨区段的钢轨粗糙度特征，对声信号进行EEMD分解后的IMF能量比进行分量筛选，利用IMF能量比畸变特征进行钢轨波磨识别。与直接法实测的钢轨粗糙度对应的理论声学频率进行了对比，提出了有效的高速铁路钢轨波磨的声学诊断策略。

Acoustic diagnosis method of rail corrugation based on IMF energy ratio

The invention relates to an acoustic diagnosis method for rail corrugation of high-speed railway based on IMF energy ratio, which belongs to the technical field of vibration and noise of high-speed railway, and the steps are as follows: (1) rail roughness test; (2) set empirical mode decomposition; (3) eigenmode function IMF energy ratio: according to the energy ratio of IMF signal corresponding to the fault characteristic frequency, the fault identification of rail corrugation is carried out, and the rail is screened Through HHT transformation, the Hilbert marginal spectrum and instantaneous frequency are obtained. Acoustic diagnosis method of rail corrugation based on IMF energy ratio. In the operation line of the invention, the rail roughness characteristics of the section with or without corrugation are measured by direct method, the IMF energy ratio after EEMD decomposition of the acoustic signal is component screened, and the rail corrugation is identified by the distortion characteristics of the IMF energy ratio. Compared with the theoretical acoustic frequency corresponding to the rail roughness measured by the direct method, an effective acoustic diagnosis strategy of rail corrugation is proposed.

【技术实现步骤摘要】
基于IMF能量比的高速铁路钢轨波磨声学诊断方法
本专利技术涉及一种高速铁路钢轨波磨声学诊断方法，特别涉及一种基于IMF能量比的高速铁路钢轨波磨声学诊断方法，属于高速铁路振动噪声

技术介绍
钢轨波磨是一种出现在钢轨表面的周期性波浪形不平顺曲线，高速铁路钢轨波磨导致车辆-轨道系统的中高频振动响应产生的振动噪声直接影响着乘客的舒适度和铁路沿线居民的生活质量，同时恶化系统各部件运行状态，加剧钢轨表面的进一步损伤。关于钢轨波磨的检测，以往采用传统波磨尺运用弦测法进行人工抽样测量，检测效率非常低。近年来，新的检测技术不断被运用，检测精度和效率都得到了提升，如钢轨粗糙度检测手推车、应用振动加速度的惯性基准法和机器视觉法等。对于应用车辆上的振动加速度信号进行波磨诊断，国内外对其故障识别算法逐渐进行了一些研究，Grassie最早提出了应用车辆轴箱加速度信号分析进行轨道的动态监测的构想；Tsunashima等通过车体振动信号的小波包分析进行轨道波磨的识别；曹西宁等通过轴箱加速度信号进行希尔伯特-黄变换对轨道不平顺进行分析和诊断。利用加速度信号等接触式测量方法准备工作相对复杂，且受限于钢轨波磨特征往往与轮轨系统的耦合振动特性相关，易淹没于轮轨动力学系统固有特征内，对信号处理算法也提出了更高的要求。从声学角度进行钢轨的波磨诊断，是一种非接触式的间接测量方法，其以列车运营状态下轮轨振动所产生的声信号作为反映轨道状态的重要信息来源，根据目标结构的声振发生机理和特征，对轨道波磨状态进行诊断，检测效率高，具有明显的早期预警和快速检测优势。高速铁路钢轨波磨初期往往幅值较小，列车在运行状态下所产生的声信号也不可避免的受到噪声的干扰，反映故障信息的脉冲信号很容易被淹没，同时由于钢轨波磨往往与高速铁路车辆-轨道耦合系统或轮轨系统部件的共振频率相关，故无波磨出现时声信号在波磨的特征频率往往也存在峰值频率，所以，运用常用的针对非稳态信号的时频分析技术，也难以在时频谱特征中准确的识别钢轨波磨特征。因此，为了解决高速铁路钢轨波磨的快速检测问题，提供一种基于IMF能量比的高速铁路钢轨波磨声学诊断方法，利用动车组在运营时速下安装于转向架上的传声器采集的声信号进行声学诊断，就成为该
急需解决的技术难题。
技术实现思路
本专利技术的目的是提供一种基于IMF能量比的高速铁路钢轨波磨声学诊断方法，利用动车组在运营时速下安装于转向架上的传声器采集的声信号进行声学诊断。本专利技术的上述目的是通过以下技术方案达到的：一种基于IMF能量比的高速铁路钢轨波磨声学诊断方法，其步骤如下：(1)钢轨粗糙度测试采用小推车直接测量钢轨粗糙度(短波不平顺)，对测试得到的表面不平度进行声学修正；(2)集合经验模态分解(EEMD)对于含有严重噪声的原始信号，根据高速铁路钢轨粗糙度特征频率，按照从高频到低频的顺序，通过重采样和滤波，分解为具有不同振动模态的子信号，获得本征模态函数IMF，有效将噪声成分分离；(3)本征模态函数IMF能量比根据故障特征频率对应IMF信号的能量比进行钢轨波磨的故障识别，筛选得到钢轨波磨对应的IMF分量，通过HHT变换，得到Hilbert边际谱和瞬时频率。优选地，所述步骤(1)中所述声学修正如下：钢轨粗糙度数据处理过程中要进行尖峰去除和曲率修正，其中曲率修正对粗糙度的微观几何特征进行声学角度的处理，以还原钢轨粗糙度对轮轨相互作用的影响；对于每一个粗糙度测试得到的点坐标所构成的实际钢轨粗糙度表面r(x)，接触点位于中心的x0处，相对于理想车轮表面，结合测试得到的车轮半径，通过曲率修正将声学粗糙度修正为r’(xi)-r(xi)。优选地，所述钢轨粗糙度数据表示成圆周长度的函数，其物理含义是不同位置处钢轨表面相对于平均表面的变化值，通常被称为钢轨不平顺幅值，理论研究中常用对数形式的钢轨不平顺等级表示，定义如式1所示(式(1)中的rref指是参考值，k指的是第k个值)，单位为dB；式(1)中，是钢轨粗糙度的均方值在1/3倍频程中进行量化，参考值取1μm，在每个1/3倍频程中将所得的窄带频谱幅值的平方再求和，并除以计算点数即可获得，在声学粗糙度的定义中，10μm粗糙度的有效幅值(均方根值)对应20dB的粗糙度等级，而1μm的粗糙度幅值则对应0dB粗糙度等级。优选地，所述步骤(2)的具体步骤如下：步骤1)：找出信号x(t)(见下式(2))的极大值点和极小值点，用样条插值函数拟合形成上包络线和下包络线，计算上下包络线的均值m1(t)，将原数据序列x(t)减去平均包络m1(t)，得到新的数据序列h1(t)，若h1(t)不满足IMF的条件，则将h1(t)作为原信号重复上面的步骤k次，使得平均包络线趋于零，得到的h1k(t)就是第一个IMF；步骤2)：从原信号减去c1(t)，得到一个新的数据序列，然后重复步骤1)，得到一系列的cn(t)和一个不可再分解的余项序列rn(t)，其中rn(t)表示信号的平均趋势；原信号则可以表示为IMF分量和一个残余项之和。优选地，所述步骤(3)的具体步骤如下：利用车下声信号进行EEMD后，提取的IMF能量比畸变特征进行波磨特征识别，IMF是EEMD方法分解后得到的本征模态分量，反映信号从低频到高频的不同振动模态，IMF的能量熵公式为：其中：pi＝Ei/E(4)pi为第i个IMF能量占总能量的比值，能量的公式为：IMF需要满足条件有两点：一是序列中，极值点数与过0点数必须相等或者最多相差一个；二是在任意时间点上，由信号局部极大值确定的上包络线和由局部极小值确定的下包络线的均值为0。有益效果：本专利技术的基于IMF能量比的高速铁路钢轨波磨声学诊断方法，在运营线路上采用直接法实测了有无波磨区段的钢轨粗糙度(短波不平顺)特征，同时采集了运营动车组通过有无波磨区段的车下声信号，针对车下声信号的特点，对声信号进行EEMD分解后的IMF能量比进行分量筛选，利用IMF能量比畸变特征进行钢轨波磨识别，与直接法实测的钢轨粗糙度对应的理论声学频率进行了对比，提出了有效的高速铁路钢轨波磨的声学诊断策略。本专利技术的基于IMF能量比的高速铁路钢轨波磨声学诊断方法，通过白噪声缓和异常事件的局部干扰，从而有效解决模态混叠问题，有效避免噪声干扰；并且在时频谱特征中准确的识别钢轨波磨特征，特征频率识别较为准确。下面通过附图和具体实施方式对本专利技术做进一步说明，但并不意味着对本专利技术保护范围的限制。附图说明图1为本专利技术基于IMF能量比的高速铁路钢轨波磨声学诊断方法中钢轨粗糙度测试示意图。图2为本专利技术基于IMF能量比的高速铁路钢轨波磨声学诊断方法中钢轨粗糙度声学修正的示意图。图3为本专利技术基于IMF能量比的高速铁路钢轨波磨声学诊断方法中钢轨波磨诊断流程示意图。图4为本专利技术基于IMF能量比的高速铁路钢轨波磨声学诊断方法中高速铁路某区段钢轨波磨测试结果。图5-1为安装于运营车辆车下的传声器监测的声信号进行滤波后得到的时域图(有磨)。图5-2为安装于运营车辆车下的传声器监测的声信号进行滤波后得到的时域图(无磨)。图6-1为声信号通过EEMD分解后钢轨有波磨状态的前4阶本征模态信号(IMF1)。图6-2为声信号通过EEMD分解后钢轨有波磨状态的前4阶本征模态信号(IMF2)。图6-3为声信号通过EE本文档来自技高网...

【技术保护点】
1.一种基于IMF能量比的高速铁路钢轨波磨声学诊断方法，其步骤如下：(1)钢轨粗糙度测试采用小推车直接测量钢轨粗糙度，对测试得到的表面不平度进行声学修正；(2)集合经验模态分解对于含有严重噪声的原始信号，根据高速铁路钢轨粗糙度特征频率，按照从高频到低频的顺序，通过重采样和滤波，分解为具有不同振动模态的子信号，获得本征模态函数IMF，有效将噪声成分分离；(3)本征模态函数IMF能量比根据故障特征频率对应IMF信号的能量比进行钢轨波磨的故障识别，筛选得到钢轨波磨对应的IMF分量，通过HHT变换，得到Hilbert边际谱和瞬时频率。

【技术特征摘要】
1.一种基于IMF能量比的高速铁路钢轨波磨声学诊断方法，其步骤如下：(1)钢轨粗糙度测试采用小推车直接测量钢轨粗糙度，对测试得到的表面不平度进行声学修正；(2)集合经验模态分解对于含有严重噪声的原始信号，根据高速铁路钢轨粗糙度特征频率，按照从高频到低频的顺序，通过重采样和滤波，分解为具有不同振动模态的子信号，获得本征模态函数IMF，有效将噪声成分分离；(3)本征模态函数IMF能量比根据故障特征频率对应IMF信号的能量比进行钢轨波磨的故障识别，筛选得到钢轨波磨对应的IMF分量，通过HHT变换，得到Hilbert边际谱和瞬时频率。2.根据权利要求1所述的基于IMF能量比的高速铁路钢轨波磨声学诊断方法，其特征在于：所述步骤(1)中所述声学修正如下：钢轨粗糙度数据处理过程中要进行尖峰去除和曲率修正，其中，曲率修正对粗糙度的微观几何特征进行声学角度的处理，以还原钢轨粗糙度对轮轨相互作用的影响；对于每一个粗糙度测试得到的点坐标所构成的实际钢轨粗糙度表面r(x)，接触点位于中心的x0处，相对于理想车轮表面，结合测试得到的车轮半径，通过曲率修正将声学粗糙度修正为r’(xi)-r(xi)。3.根据权利要求2所述的基于IMF能量比的高速铁路钢轨波磨声学诊断方法，其特征在于：所述钢轨粗糙度数据表示成圆周长度的函数，其物理含义是不同位置处钢轨表面相对于平均表面的变化值，被称为钢轨不平顺幅值，用对数形式的钢轨不平顺等级表示，定义如式1所示，单位为dB；式(1)中，是钢轨粗糙度的均方值在1/3倍频程中进行量化，参考值取1μm，在每个1/3倍频程中将所得的窄带频谱...

【专利技术属性】
技术研发人员：韩立，伍向阳，刘兰华，李晏良，张毅超，陈迎庆，邵琳，
申请(专利权)人：中国铁道科学研究院集团有限公司节能环保劳卫研究所，
类型：发明
国别省市：北京,11