发明名称 基于非负张量分解振动特征的钢轨伤损探测装置和方法
摘要 基于非负张量分解振动特征的钢轨伤损探测装置和方法,涉及高铁钢轨的伤损探测装置和方法,解决了现有的探伤技术对伤损频率接近的伤损情况只能通过幅值来判断,伤损特征难以区分的问题,基于非负张量分解振动特征的钢轨伤损探测装置,包括振动加速度传感器阵列、信号转换电路、信号特征提取模块、伤损匹配模块和伤损判别模块;基于非负张量分解振动特征的钢轨伤损探测方法,包括:步骤一:采集钢轨的原始振动信号;步骤二:转换原始信号为振动信号;步骤三:将振动信号进行HHT分解;步骤四:利用得到频率幅值之间的幅频信号,建立伤损信号的三维张量;步骤五:利用非负张量分解对构建的三维张量信号分解,提取特征信息;用于高铁钢轨伤损探测。
申请公布号 CN102353717B 申请公布日期 2013.03.13
申请号 CN201110176806.2 申请日期 2011.06.28
申请人 哈尔滨工业大学 发明人 王艳;章欣;冯乃章;沈毅
分类号 G01N29/04(2006.01)I;G01N29/44(2006.01)I 主分类号 G01N29/04(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 牟永林
主权项 1.基于非负张量分解振动特征的钢轨伤损探测方法,其特征是它包括具体步骤如下:步骤一:振动加速度传感器阵列(1)采集到被测钢轨的原始振动信号S(t);步骤二:信号转换电路(2)转换原始振动信号S(t)为振动信号x(t);步骤三:将振动信号x(t)进行希尔伯特-黄变换分解,得到各阶内固模态函数和残差,记为:c<sub>1</sub>,c<sub>2</sub>,...c<sub>n</sub>以及r,然后进行相应的希尔伯特变换,得到频率与幅值之间的幅频信号h(w);步骤四:利用希尔伯特-黄变换分析所得到频率幅值之间的幅频信号h(w),结合相应的测量点和伤损种类建立伤损信号的三维张量;步骤五:利用非负张量分解对构建的三维张量信号x<sup>frequency×point×class</sup>进行分解,提取相应的伤损特征信息,得到每种典型伤损情况相应的特征向量,实测振动信号的特征向量通过与典型伤损情况特征向量的对比,判断出伤损的存在及其位置;步骤五包括的子步骤如下:通过步骤四建立的三维张量信号表示为:X<sup>frequency×point×class</sup>=G×<sub>1</sub>A<sup>(1)</sup>×<sub>2</sub>A<sup>(2)</sup>×<sub>3</sub>A<sup>(3)</sup>其中G,A<sup>(1)</sup>,A<sup>(2)</sup>,A<sup>(3)</sup>≥0,利用非负张量分解对其进行分解;子步骤1、随机初始化A<sup>(n)</sup>,n=1,2,3;算最小化的最小均方误差C<sub>new</sub>=‖X<sub>(n)</sub>-A<sup>(n)</sup>Z<sup>(n)T</sup>‖<sup>2</sup>,其中:X<sub>(n)</sub>为原始信号,<maths num="0001"><![CDATA[<math><mrow><msup><mi>Z</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup><mo>=</mo><msup><mi>A</mi><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></msup><mo>|</mo><mo>&CircleTimes;</mo><mo>|</mo><msup><mi>A</mi><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>|</mo><mo>&CircleTimes;</mo><mo>|</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>|</mo><mo>&CircleTimes;</mo><mo>|</mo><msup><mi>A</mi><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>|</mo><mo>&CircleTimes;</mo><mo>|</mo><msup><mi>A</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>|</mo><mo>&CircleTimes;</mo><mo>|</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>|</mo><mo>&CircleTimes;</mo><mo>|</mo><msup><mi>A</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>,</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>3</mn><mo>;</mo></mrow></math>]]></maths>子步骤2、迭代求解<img file="FDA00002530446900012.GIF" wi="491" he="121" />并计算C<sub>new</sub>=‖X<sub>(n)</sub>-A<sup>(n)</sup>Z<sup>(n)T</sup>‖<sup>2</sup>,n=1,2,3;子步骤3、判断<img file="FDA00002530446900013.GIF" wi="272" he="103" />δ为误差精度要求,成立则C<sub>old</sub>=C<sub>new</sub>,继续执行子步骤2;不成立则执行子步骤4;子步骤4、得到满足要求的A<sup>(1)</sup>,A<sup>(2)</sup>,A<sup>(3)</sup>矩阵;其中:A<sup>(1)</sup>,A<sup>(2)</sup>为新的基底,其物理意义分别代表频率和测量点;A<sup>(3)</sup>为伤损类别的特征系数矩阵,其中每一行是一个特征向量,对应一种伤损情况;通过与典型伤损情况特征向量的对比,判断出伤损的存在及其位置。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号