一种基于表面磁记忆信号的无缝铁轨内部温度应力早期诊断方法技术

本发明专利技术属于高速铁路无缝线路轨道稳定性的检测和评估领域，具体涉及一种基于表面磁记忆信号的无缝铁轨内部温度应力早期诊断方法。包括提出一种铁路现场测量方案对无缝铁轨内部在不同程度温度应力积累下表面金属磁记忆信号进行采集；对各帧磁记忆信号进行时域特征分析并进行主元分析PCA特征选优；利用傅立叶变换以及功率谱分析确定小波分解层数；最大小波能量‑香农熵比值确定小波基；对信号进行小波包分解并提取其子频带能量信息作为磁记忆信号频域特征；结合时频域特征，利用机器学习算法建立数据模型并对铁路内部温度应力积累程度进行诊断。本发明专利技术能够快速判断无缝线路内部应力累积程度，实现铁轨损伤的早期预防，减少事故的发生，具有实际的工程应用价值。

An early diagnosis method for internal temperature stress of seamless rail based on surface magnetic memory signal

【技术实现步骤摘要】
一种基于表面磁记忆信号的无缝铁轨内部温度应力早期诊断方法
本专利技术属于高速铁路无缝线路轨道稳定性的检测和评估领域,具体涉及一种基于表面磁记忆信号的无缝铁轨内部温度应力早期诊断方法。
技术介绍
我国铁路事业飞速发展，无缝铁轨被广泛的应用在高速铁路建设中，其特点是每段轨之间以焊接形式连接，因此因气温变化时铁轨的热胀冷缩被限制，其内部往往会积累巨大的温度应力，加上高速行驶的列车碾压造成的残余应力，如果不能及时处理，当应力累积到一定程度时，铁轨易发生扭曲甚至断裂，因此，对无缝线路内部应力积累情况的快速早期评估具有现实意义。金属磁记忆技术是新兴的唯一被动无损检测方法，其原理可以简明阐述为铁磁性工件受工作载荷和地磁场的共同作用，在应力和变形集中区域会发生磁致伸缩性质的磁畴组织定向、不可逆的重新取向，形成漏磁场，其区别于其他无损检测手段最主要的方面是其不仅能检测出铁磁构件宏观缺陷，同时也能反映工件内部应力集中及早期损伤。自1999年俄国学者DOUBV将此技术在中国提出后，国内大量学者做了很多磁记忆机理及应用方面的研究，但大都针对工件或工业设备焊接裂缝及宏观裂纹，其在交通尤其新兴无缝线路的隐性损伤诊断甚少。在铁路检测方面，超声波检测技术应用广泛，但其只能检测已存在的宏观缺陷，无法通过检测内部应力对其寿命进行估计。
技术实现思路
针对目前磁记忆检测技术在无缝线路检测方面应用甚少，且在铁磁构件应力检测中仅能定性检测,达不到定量检测的，本专利技术的目的是提供一种铁磁构件应力定量检测方法,通过提取无缝铁轨内部不同程度温度应力积累情况下表面金属磁记忆信号特征的方法，将无缝线路内部温度应力状况与其表面金属磁记忆信号特征联系起来，并建立预测模型，可用于无缝线路沿线检测对铁路内部应力积累状况的评估及早期预警。为达到上述目的,本专利技术采用的技术方案是:对不同温度应力下法线方向上的磁记忆信号进行采集；对所得的磁记忆信号进行时域分析并进行主元分析PCA,进行优选，得到时域特征向量；对所得的磁记忆信号进行频域分析，得到频域特征向量；将得到的时频域特征向量汇总得到最终特征向量；用提取的各特征向量建立训练集，并训练预测模型；将需要诊断的检测数据输入预测模型，得到定量检测结果。进一步,所述磁记忆信号的采集.需要在测量的无缝线路段每隔一段距离进行标记，采集各标记点轨腰处内部温度应力值(单位Mpa)及该点前后0.5m范围内法线方向磁记忆信号为一帧数据(单位：A/mm)，如图2。进一步,所述得到的不同应力下各帧数据x＝{x1,x2,...,xN},时域统计量包括：“均方根”、“极差”、“偏度”、“峰度”、“峭度”、“整流平均值”、“方差”及“标准差”，其计算公式分别为：均方根：极差：pd＝xmax-xmin(2)偏度：峰度：峭度：整流平均值：方差及标准差：初步得到时域特征向量：T＝[rms,pd,ske,ku,kur,av,va,st]。PCA降维处理将上述时域特征向量T由8维降至k'维：t＝[pc(m)],(m＝1,2,…,k')，k'由主元分析PCA,分析后贡献较大的成分个数决定。进一步,所述得到的不同应力下各帧数据x＝{x1,x2,...,xN}进行离散傅立叶变换并求出其功率谱重心，观察其主要幅值所处频率段fd。；得到的频谱分析结果，确定小波分解层数J；计算在不同小波基下磁记忆信号的能量-香农熵比，选择能量-香农熵比比值最大的为合适小波基。根据所得小波基与所确定的小波分解层数，对得到的磁记忆信号进行J层分解及重构，得到2n个子频带并求出各子频带能量，取前四个子频带信号能量作为磁记忆信号频域特征::k取值由实际情况决定。进一步,将所述得到的时域特征向量与频域特征向量合并：为最终提取的信号特征。计算公式如下:离散傅里叶变换公式为：其中N为采样序列长度。功率谱重心公式为：小波分解层数J确定公式为：其中fs为采样频率，fd为需要识别的频率，在实际应用中，可以在此基础上增加1层提高分辨率。小波能量-香农熵比值公式为：其中为小波分解第i级的能量，di(m)为信号小波分解第m个小波系数，M为信号小波分解第i级小波系数的总个数，为2n个，n为小波分解层数，且有：其中是小波系数的能量概率密度，且若pm＝0,则相应的pmlogapm＝0。其中为小波包树底层各子频带重构信号，xjk为离散点的幅值且有总信号进一步,现有大部分机器学习算法均可用于训练数据模型，如SVM，ELM等。与现有技术相比，本专利技术的有益效果是：1、为磁记忆技术在无缝线路的内部应力定量研究方向打下基础，为无缝线路沿线检测及维护提供了新的方法思路。2、本专利技术在随机信号处理基础上，从时域与频域两方面分析不同应力下的磁记忆信号特征，建立数据模型，能够快速判断无缝线路内部应力累积程度，实现铁轨损伤的早期预防，减少事故的发生，具有实际的工程应用价值。附图说明图1是本专利技术实施例的一种基于表面磁记忆信号的无缝铁轨内部温度应力早期诊断方法流程图。图2是本专利技术各帧数据采集示意图。具体实施方式以下结合具体实施例对本专利技术作进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本专利技术，并不用于限定本专利技术。本专利技术中具体实验数据来源为中国河北保定工务段保满线下行，长度为1100m无缝线路，所测量的固定区长度为900m。磁记忆信号由俄国动力诊断公司的TSC-1M-4型仪器采集，各测量点内部应力值由本实验室自主开发且拥有自主产权的金属磁巴克豪森噪声应力仪测量。分析软件环境为Matlab2016a的编程脚本下进行。本专利技术方法具体实施步骤如下所述：步骤1，现场获取数据：在上述公务路段每隔一段距离进行标记，采集各标记点轨腰处内部温度应力值(单位Mpa)及该点前后0.5m范围内法线方向磁记忆信号为一帧数据(单位：A/mm)。步骤2，将获取的磁记忆信号进行时域分析，根据式(1)～(8)计算各统计量，得到原始时域特征向量T，通过PCA主成分分析，时域特征量中标准差与极差对主成分贡献最大，故去掉其他时域特征得到最终时域特征向量：t＝[Ran,St]；步骤3.1,对步骤1获取的磁记忆信号按式(9)，(10)进行傅里叶变换并求出其功率谱重心，得到信号能量集中在在低频段0～10Hz。步骤3.2，通过式(11)，确定小波分解级数为3层。步骤3.3，通过式(12)，(13)，选定能量-香农熵比值最大的小波基为Demy小波步骤3.4，利用Matlab对步骤1得到的磁记忆信号进行J层分解及重构，得到2n个子频带并按式(14)求出各子频带能量，取前4个子频带信号能量作为磁记忆信号频域特征:步骤4，将步骤2与步骤3中得到的时频域特征向量汇总得到最终特征向量：步骤5，本次实验使用K-ELM训练预测模本文档来自技高网...

【技术保护点】
1.一种基于表面磁记忆信号的无缝铁轨内部温度应力早期诊断方法，其特征是，包括以下步骤：/n步骤1：对不同温度应力下法线方向上的磁记忆信号进行采集；/n步骤2：对步骤1所得的磁记忆信号进行时域分析并进行主元分析(PCA)优选，得到时域特征向量；/n步骤3：对步骤1所得的磁记忆信号进行频域分析，得到频域特征向量；/n步骤4：将步骤2与步骤3中得到的时域和频域特征向量汇总得到最终特征向量；/n步骤5：将步骤4中提取的各最终特征向量建立训练集，并训练预测模型；/n步骤6：将需要诊断的检测数据输入步骤5中得到的预测模型中，生成检测结果。/n

【技术特征摘要】
1.一种基于表面磁记忆信号的无缝铁轨内部温度应力早期诊断方法，其特征是，包括以下步骤：
步骤1：对不同温度应力下法线方向上的磁记忆信号进行采集；
步骤2：对步骤1所得的磁记忆信号进行时域分析并进行主元分析(PCA)优选，得到时域特征向量；
步骤3：对步骤1所得的磁记忆信号进行频域分析，得到频域特征向量；
步骤4：将步骤2与步骤3中得到的时域和频域特征向量汇总得到最终特征向量；
步骤5：将步骤4中提取的各最终特征向量建立训练集，并训练预测模型；
步骤6：将需要诊断的检测数据输入步骤5中得到的预测模型中，生成检测结果。


2.根据权利要求1所述的方法，其特征在于，所述步骤1包括：
在需要测量的无缝线路段每隔一段距离进行标记，采集各标记点轨腰处内部温度应力值(单位Mpa)及该点前后0.5m范围内法线方向磁记忆信号为一帧数据(单位：A/mm)。


3.根据权利要求1所述的方法，其特征在于，所述步骤2包括：
求出步骤1得到的不同应力下,各帧数据x＝{x1,x2,...,xN}时域统计量，包括：“均方根RMS”、“极差pd”、“偏度ske”、“峰度Ku”、“峭度kur”、“整流平均值av”、“方差va”及“标准差st”，初步得到时域特征向量：T＝[RMS,pd,sk...

【专利技术属性】
技术研发人员：陈娟，黄顺昊，杨巧宁，
申请(专利权)人：北京化工大学，
类型：发明
国别省市：北京;11