CN116252820A - 改进频域积分法驱动的高速列车车轮多边形定量检测方法 - Google Patents

改进频域积分法驱动的高速列车车轮多边形定量检测方法 Download PDF

Info

Publication number
CN116252820A
CN116252820A CN202310536835.8A CN202310536835A CN116252820A CN 116252820 A CN116252820 A CN 116252820A CN 202310536835 A CN202310536835 A CN 202310536835A CN 116252820 A CN116252820 A CN 116252820A
Authority
CN
China
Prior art keywords
wheel
signal
aba
polygon
frequency domain
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202310536835.8A
Other languages
English (en)
Other versions
CN116252820B (zh
Inventor
谢清林
陶功权
温泽峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202310536835.8A priority Critical patent/CN116252820B/zh
Publication of CN116252820A publication Critical patent/CN116252820A/zh
Application granted granted Critical
Publication of CN116252820B publication Critical patent/CN116252820B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B61RAILWAYS
    • B61KAUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
    • B61K9/00Railway vehicle profile gauges; Detecting or indicating overheating of components; Apparatus on locomotives or cars to indicate bad track sections; General design of track recording vehicles
    • B61K9/12Measuring or surveying wheel-rims

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明涉及车轮多边形检测技术领域,涉及一种改进频域积分法驱动的高速列车车轮多边形定量检测方法,其包括以下步骤:1)从运行中的列车上采集轴箱加速度ABA和速度信号,并采用滑动窗口技术将数据划分为不同的片段;2)引入稳态试验来计算每个片段的一致相关系数CCC,保存满足阈值要求的片段;3)采用相似性试验从保存的片段中精确提取代表每个车轮旋转周期的准平稳ABA信号;4)对提取的准平稳ABA信号进行改进的频域积分分析以定量检测车轮多边形。本发明能较佳地进行车轮多边形定量检测。

Description

改进频域积分法驱动的高速列车车轮多边形定量检测方法
技术领域
本发明涉及车轮多边形检测技术领域,具体地说,涉及一种改进频域积分法驱动的高速列车车轮多边形定量检测方法。
背景技术
铁道车辆依靠轮轨滚动接触完成载运任务。然而,轮轨之间复杂的相互作用会导致各种车轮不圆度(OOR)缺陷,严重影响铁道车辆的运行安全性与舒适性。
车轮多边形对乘客和车辆-轨道-桥梁/路基系统都有很大的影响,例如,降低乘客舒适度,对沿线居民造成严重的噪声污染,恶化轮轨接触行为并威胁行车安全性和可靠性等。因此,开发一种高效实用的车轮多边形检测方法,对于保障铁路车辆安全、可靠且经济地运行具有重要意义。目前,车轮多边形的测量多采用车轮粗糙度测量仪或车轮镟修镟床。虽然这一类车间内测试方法可以获得较好的精度,但其运行效率较低,无法满足现代铁路的运营维护需求。此外,关于车轮 OOR 线路式检测方法的研究主要集中在车轮扁疤或局部缺陷,只有少数研究涉及到车轮多边形检测,尤其是动态检测方法。已报告的动态车轮多边形测量可分为两类,即轨旁检测和车载检测。
轨旁检测:此类方法是指在一个/多个轨道段上方/附近安装传感器,以收集和分析车轮经过时产生的各种信号。有学者综述总结了轨旁检测式的车轮 OOR 监测方法,并将其分为基于应变、加速度与力学廓形式的3种类型。有学者报道了一种在轨道上布置应变片来获得轮轨接触力以评估车轮 OOR 状态的方法。然而,结果表明该类方法更适用于车轮扁疤和较大尺寸的车轮局部缺陷检测。 有学者对原始的轨旁检测方法进行了改进以检测车轮多边形。此外,有学者在轨道上安装了加速度计来采集轨道振动信号。他们分别采用模糊逻辑和数理统计算法来评估车轮缺陷。有学者基于轮缘顶部不易磨损的原理,利用轮缘的高度来计算车轮多边形的轮廓。除了上述方法,轨旁检测还包括基于其他传感器类型的方法,例如光纤传感技术、机器视觉、声发射设备和超声波技术。然而,由于传感器安装位置和数量的限制,以及轮径、轨枕间距、轮轨接触位置的影响,轨旁检测方法的准确性仍有待提高。此外,轨旁检测方法无法在线监测车轮缺陷。
车载检测:许多学者通过总结现场跟踪测试和高保真数值模拟数据,揭示了车轮多边形与轴箱加速度(ABA)之间的关系。因此,ABA 可以作为检测车轮多边形的良好载体。有学者提出了一种角域同时平均技术的车载检测方法,通过分析垂向ABA来检测铁路车轮多边形。然而,该技术通过构建间接健康指数来反映车轮多边形的严重程度,而不是直接识别车轮多边形的幅值或粗糙度特征。有学者基于ABA与车轮多边形的相关性,定量研究了不同阶次和幅值的车轮多边形激发的ABA特性。该方法受限于车辆恒速运行的要求。有学者开发了一种联合时频脊线检测和信号自适应分解的方法,从齿轮箱振动信号中获取车辆运行速度。然后他们结合重采样和阶次分析技术来处理 ABA信号,以定量诊断车轮多边形阶次和幅值信息。该方法涉及复杂的数学算法,参数设置不当可能会影响检测结果的准确性。鉴于 ABA信号的非线性和非平稳特性,有学者报道了两种动态检测框架,分别基于修正离散傅里叶变换和参数功率谱估计,用于检测车轮多边形。除了传统的信号处理工具外,有学者还建立了一个名为OORNet的基于深度学习的模型,用于高速列车车轮OOR监测,然而该研究是由仿真数据而不是实测信号驱动的。
从上可知,车轮多边形检测已经成为铁路工程中急需解决的关键问题。通过及时并准确地获悉车轮多边形状态,不仅可使研究人员能进一步了解车轮OOR的形成机理和演变趋势,而且铁路运营人员也可以制定更加科学的车轮维护策略。然而,基于ABA的车载检测或轨旁检测方法来定量检测铁路车轮多边形是一项具有挑战性的任务。造成该困境的主要原因是从复杂环境中采集的信号具有强烈的非平稳属性。鉴于上述不足,本技术提出了一种基于ABA信号的高速铁路车轮多边形车载式动态检测方法,以定量检测车轮多边形的阶次与粗糙度特征。
发明内容
本发明的内容是提供一种改进频域积分法驱动的高速列车车轮多边形定量检测方法,其能够定量检测不同运行里程下的车轮多边形阶次与粗糙度特征。
根据本发明的一种改进频域积分法驱动的高速列车车轮多边形定量检测方法,其包括以下步骤:
1)从运行中的列车上采集轴箱加速度ABA和速度信号,并采用滑动窗口技术将数据划分为不同的片段;
2)引入稳态试验来计算每个片段的一致相关系数CCC,保存满足阈值要求的片段;
3)采用相似性试验从保存的片段中精确提取代表每个车轮旋转周期的准平稳ABA信号;
4)对提取的准平稳ABA信号进行改进的频域积分分析以定量检测车轮多边形。
作为优选,1)中,具体包括以下步骤:
1.1)轴箱加速度ABA和速度信号采集;
1.2)进行带通滤波;
1.3)使用滑动窗将ABA信号切割成一系列数据片段;
1.4)使用子滑动窗将切割后的数据片段再次细分为3个数据片段(S m1S m2S m3)。
作为优选,2)中,具体包括以下步骤:
2.1)执行频谱和概率密度分布分析以计算数据片段(S m1S m2S m3)的频谱序列(F 1F 2F 3)和概率密度分布序列(P 1P 2P 3);
2.2)分别计算F 1F 2F 1F 3F 2F 3P 1P 2P 1P 3P 2P 3的CCC值C 1C 2C 3C 4C 5C 6,计算公式如下:
Figure SMS_1
Figure SMS_2
其中C表示一致相关系数CCC,C越接近1,序列AB之间的相关性越高,ρ为皮尔逊相关系数,cov表示求协方差,σ A σ B 分别为序列AB的标准差;
Figure SMS_3
分别为序列AB的向量;/>
Figure SMS_4
和/>
Figure SMS_5
分别为/>
Figure SMS_6
的均值;
通过一致相关系数CCC同时评估两个序列之间的相关性和误差水平,只要(F 1F 2F 3),(P 1P 2P 3)之间的相关性和误差水平满足阈值,就认为数据片段相对稳定;
2.3)计算CCC序列的平均值MCMC m = mean (C 1 C 2 C 3 C 4 C 5 C 6) ;数据片段的MC值被存储在临时组中,用于后续的相似性测试。
作为优选,3)中,具体包括以下步骤:
3.1)将MC值大于阈值T的数据片段保存;
3.2)利用基于迭代算法驱动的相似性试验从保存的数据片段中精确提取代表车轮每个车轮周期的准平稳ABA信号,方法如下:
输入:最小速度V min,最大速度V max,最小半径R min,最大半径R max
采样频率Fs
输出:整旋转周期的信号序列y
Figure SMS_7
其中,T min是车轮旋转一周最短时间,T max 是车轮旋转一周最长时间,j是车轮旋转一周的信号采样点数,SP min是车轮旋转一周最小信号采样点数,SP max 是车轮旋转一周最大信号采样点数;迭代算法的输入是最小速度V min、最大速度V max、最小车轮半径R min、最大车轮半径R max和已知的采样频率Fs;输出y表示预期的整周期ABA信号;x 1x 2x 3是三个端到端的具有不同数据采样点的信号序列,ρ 1 ρ 2ρ 3是对应x 1x 2x 3的三个皮尔逊相关系数,type 表示类型pearson表示皮尔逊相关系数;y 1y 2y 3是对应x 1x 2x 3的三个输出;初始化的功能是基于输入变量获得对应于一个车轮旋转的最小和最大采样点SP minSP max;迭代算法的思想是在平稳性测试后保存的数据片段S m中生成三个端到端的具有不同数据采样点的信号序列(x 1x 2x 3),并计算三个序列之间的平均皮尔逊相关系数ρ,利用周期信号之间的高相似性,最大ρ对应的SP’正是一个车轮转动周期所对应的ABA信号的长度,由此能提取代表每个车轮旋转周期的ABA信号。
作为优选,4)中,改进的频域积分分析采用一种改进的自适应带通滤波的频域积分方法,获得的速度信号ν if(t)和位移信号s if(t),公式为:
Figure SMS_8
Figure SMS_9
Figure SMS_10
Figure SMS_11
Figure SMS_12
其中H(ω)是自适应带通滤波器,
Figure SMS_13
是频率分辨率,A(ω)是加速度信号a(t)经过傅里叶变换后的频域信号,N是采样点数,当分析第k阶车轮多边形时,f lk f uk 是带通滤波器的下截止频率和上截止频率,以f k 为中心的带通滤波器能降低干扰,只提取与第k阶车轮多边形相关的ABA信息;t是时间,v是车辆行驶速度,ω是角频率;R是车轮名义滚动圆处半径;K是总阶数;
基于改进的频域积分方法获得对应于车轮多边形的每一阶的位移信号,就能定量计算每一阶的粗糙度水平;
最后,统计所有保存的数据片段的检测结果以统计其对应的分布特征和数学期望。
本发明提出了一个新颖的数据驱动的检测框架,以便及时并准确地检测车轮多边形。具体地,设计了适当的信号预处理方法、平稳性测试和相似性测试,以精确提取代表每个车轮旋转周期的准静态ABA信号。在提取的ABA信号驱动下,基于惯性基准法和提出的改进频域积分方法,可以获得车轮多边形的定量表达,即阶次和粗糙度特征。采用一个车轮镟修周期内的现场跟踪测试评估了本技术的有效性。本发明提出的高速铁路车轮多边形检测框架能够定量检测不同运行里程下单阶和多阶车轮多边形的主要特征。
附图说明
图1为实施例中一种改进频域积分法驱动的高速列车车轮多边形定量检测方法的流程图;
图2为实施例中ABA信号预处理的示意图;
图3(a)为实施例中频谱曲线的示意图;
图3(b)为实施例中概率密度曲线的示意图;
图4为实施例中迭代算法驱动的相似性试验释义的示意图;
图5为实施例中实测ABA与速度信号 、 提取的准静态ABA信号的示意图;
图6(a)为实施例中车轮镟修后运营里程为50 259 km下的单个主导阶次车轮多边形检测结果的示意图;
图6(b)为实施例中车轮镟修后运营里程为81 694 km下的单个主导阶次车轮多边形检测结果的示意图;
图6(c)为实施例中车轮镟修后运营里程为105 172 km 下的单个主导阶次车轮多边形检测结果的示意图;
图7(a)为实施例中车轮镟修后运营里程为50 259 km下的多个主导阶次车轮多边形检测结果的示意图;
图7(b)为实施例中车轮镟修后运营里程为81 694 km下的多个主导阶次车轮多边形检测结果的示意图;
图7(c)为实施例中车轮镟修后运营里程为105 172 km 下的多个主导阶次车轮多边形检测结果的示意图。
具体实施方式
为进一步了解本发明的内容,结合附图和实施例对本发明作详细描述。应当理解的是,实施例仅仅是对本发明进行解释而并非限定。
如图1所示,本实施例提供了一种改进频域积分法驱动的高速列车车轮多边形定量检测方法,其包括以下步骤:
1)从运行中的列车上采集轴箱加速度ABA和速度信号,并采用滑动窗口技术将数据划分为不同的片段;
2)引入稳态试验来计算每个片段的一致相关系数CCC,保存满足阈值要求的片段;
3)采用相似性试验从保存的片段中精确提取代表每个车轮旋转周期的准平稳ABA信号;
4)对提取的准平稳ABA信号进行改进的频域积分分析以定量检测车轮多边形。
车轮多边形基本特征
车轮多边形定义为周期性偏差,表现为沿车轮圆周方向的不均匀磨损。通常,采用离散傅立叶变换来处理车轮多边形的周向偏差,该周向偏差可以表示为一系列具有不同频率(阶次)、振幅(粗糙度)和相位的空间谐波。车轮多边形由几个具有不同属性的谐波组成,可以表示为
Figure SMS_14
(1)
其中z是车轮的径向偏差,t是时间,k是车轮多边形的阶次,K是最大阶次。A k φ k 分别是第k阶车轮多边形的振幅和相位。R是车轮名义滚动圆处半径,v是车辆行驶速度。当k= 1和2时,分别表示一阶和二阶车轮多边形,称为车轮偏心与椭圆。
车轮多边形动态特性
车轮多边形对车辆的动态影响与轨道不平顺一致。根据广泛使用的惯性基准法原则,通过计算车轮垂直位移(车轮多边形的轮廓)的二阶导数可以获得垂向ABA,a(t)。其数学表达式为
Figure SMS_15
(2)
Figure SMS_16
为/>
Figure SMS_17
的二阶导数,/>
Figure SMS_18
t时刻的车轮的径向偏差,注意到车轮旋转基频f 0=v/2πR,因此对应于第k阶车轮多边形的特征频率为f k =kv/2πR,根据角频率ω和频率f之间的关系,可以获得ω k = 2πf k 。因此,由第k阶车轮多边形(z k )激发的垂向ABA(a k )可以表示为
Figure SMS_19
(3)
Figure SMS_20
z k 的二阶导数;
应该注意的是,除了由车轮多边形引起的激励之外,ABA还包括其他干扰成分,例如由结构振动、系统共振和轨道随机不平顺引起的干扰成分。因此,本技术旨在开发可靠的信号处理算法,从ABA信号中准确检测车轮多边形。
数据采集与信号预处理
在诊断车轮缺陷时,通过分析车轮3 ~ 6个旋转周期的信号,可以减少随机干扰对ABA的影响。通过滑动窗口技术将实测ABA信号分成若干片段(S 1S 2,...,S m , ...,S n ) 以进行针对性分析,如图2所示。然后,采用带通滤波器滤除低频直流(DC)分量和高频噪声。一般地,通常只关注前40阶车轮多边形,本研究涉及的高速列车的运行速度约为200~250 km/h,新车轮的直径为920 mm。因此,建议带通滤波器设置为2~1024 Hz,滑动窗口的长度和步长分别设置为0.25 s和0.05 s,这不仅可以有效地减少随机干扰,确保较高的时间分辨率,还可以提取ABA主要特征。此外,证明了利用子窗口细分每个切片片段并分析这些细分对象(S m1S m2,和S m3)的属性,可以实现更鲁棒的ABA信号处理。在本技术中,子窗口长度和子步长分别设置为0.15 s和0.05 s。
稳态试验
平稳信号的频谱和概率密度分布不应随时间变化,在一定的置信度内,通过一些适当的平稳性测试,可以从ABA测量数据中提取相对平稳的信号。具体地,执行频谱和概率密度分布分析以计算细分对象 (S m1S m2,和S m3)的频谱序列(F 1F 2,和F 3)和概率密度分布序列(P 1P 2,和P 3),所获得的曲线呈现在图3(a)和图3(b)中。只要F 1F 2F 3P 1P 2P 3之间的相关性和误差水平满足合理的阈值,就可以认为片段(S m )相对稳定。通常,皮尔逊相关系数和均方误差分别用于评估两个序列之间的一致性和误差水平。这里使用了一个更全面的一致相关系数(CCC),它可以同时评估两个序列之间的相关性和误差水平。其计算过程为
Figure SMS_21
(4)
Figure SMS_22
(5)
其中C表示一致相关系数CCC,C越接近1,序列AB之间的相关性越高,ρ为皮尔逊相关系数,cov表示求协方差,σ A σ B 分别为序列AB的标准差。
Figure SMS_23
分别为序列AB的向量;/>
Figure SMS_24
和/>
Figure SMS_25
分别为/>
Figure SMS_26
的均值。根据方程(4)和(5),分别计算F 1F 2F 1F 3F 2F 3P 1P 2P 1P 3P 2P 3的CCC值C 1C 2C 3C 4C 5C 6。然后,计算CCC序列的平均值MCMC m = mean(C 1 C 2 C 3 C 4 C 5 C 6)。每个切片片段(S 1S 2,...,S m ,...,S n )的MC值被存储在临时组中,用于后续的相似性测试。/>
通过将MC值与阈值(即MC m 是否大于阈值T)进行比较,可以保存具有较高平稳性的数据片段。注意到列车的实际运行环境非常复杂。因此,阈值T应该是自适应的,而不是恒定值。在本技术中,通过多次实验将阈值T设置为MC序列的99%分位数,它不仅可以根据不同的ABA信号自适应地改变,而且可以确保保存的片段是准平稳信号。对应于保存的片段信号可以被认为是列车以相对平稳的速度在相对平滑的轨道上运行。然而,一段信号通常包含车轮3~6个旋转周期。精确提取对应于每个车轮旋转周期的信号以获得精确的车轮多边形检测结果是非常重要的。否则一个车轮旋转周期内的数据采样点会漏到其他周期,导致车轮多边形与ABA信号的映射关系出现混淆。通过测量车轮半径和车辆速度来估计车轮转动周期是可行的。然而,用这种方法得到的结果只是真实值的近似值。因为轮轨接触位置并不总是保持在名义滚动圆上。当车辆通过小半径曲线时,车轮接触点的横向位置会发生显著变化。此外,车辆的行驶速度在实际运行中总是有轻微的波动,不可能是理想的匀速状态。因此,本研究提出了一种基于迭代算法驱动的相似性试验,以精确提取代表每个车轮旋转周期的准平稳ABA信号,其伪代码如下描述。
输入:最小速度V min,最大速度V max,最小半径R min,最大半径R max
采样频率Fs
输出:整旋转周期的信号序列y
Figure SMS_27
其中,T min是车轮旋转一周最短时间,T max 是车轮旋转一周最长时间,j是车轮旋转一周的信号采样点数,SP min是车轮旋转一周最小信号采样点数,SP max 是车轮旋转一周最大信号采样点数;迭代算法的输入是易获悉的最小速度V min、最大速度V max、最小车轮半径R min、最大车轮半径R max和已知的采样频率Fs。输出y 表示预期的整周期ABA信号。x 1x 2x 3是三个端到端的具有不同数据采样点的信号序列,ρ 1 ρ 2ρ 3是对应x 1x 2x 3的三个皮尔逊相关系数,type 表示类型pearson表示皮尔逊相关系数;y 1y 2y 3是对应x 1x 2x 3的三个输出;初始化的功能是基于输入变量获得对应于一个车轮旋转的最小和最大采样点(SP minSP max)。迭代算法的思想是在S m(平稳性测试后保存的数据片段)中生成三个端到端的具有不同数据采样点 (SP minSP max)的信号序列(x 1x 2x 3),并计算三个序列之间的平均皮尔逊相关系数ρ。利用周期信号之间的高相似性,最大ρ对应的SP’正是一个车轮转动周期所对应的ABA信号的长度。由此可以精确地提取代表每个车轮旋转周期的ABA信号。图4用三个周期的模拟信号说明了上述过程。可以看出,所提出的迭代算法在不需要每个时刻的速度和车轮半径的精确信息下,便可以精确地识别代表每个车轮旋转周期的ABA信号。在迭代算法中有两个问题应予以注意。首先,S m 通常包括车轮的3 ~ 6个旋转周期,因此只从S m 中提取前三个周期以保证算法的通用性。第二,因为生成的序列已经通过了平稳性测试,所以使用皮尔逊相关系数而不是CCC来提高计算效率。
车轮多边形定量检测
基于惯性基准法原理,可以通过适当处理提取的准平稳ABA信号来检测车轮多边形。通常,可以通过在时域或频域中对加速度信号进行积分来产生位移信号。时域积分法主要依赖于离散信号的数值积分,常用的有辛普森法则和梯形公式。然而,时域积分计算不可避免地会引起趋势项和DC分量的干扰。例如,假设a(t)是实测的加速度信号,其数学表达式为
Figure SMS_28
(6)
其中g(t)代表加速度的真实值,通常为谐波信号。ε是测试误差,包括测量过程中的各种噪声。公式(7)和(8)给出了基于时域数值积分的速度信号ν(t)和位移信号s(t)的一般表达式。公式(9)和(10)进一步描述了通过辛普森法则得到的速度信号ν t(t)和位移信号s t(t)。
Figure SMS_29
(7)
Figure SMS_30
(8)
Figure SMS_31
(9)
Figure SMS_32
(10)
其中t= 0,1,...,N-1,N是采样点数。
Figure SMS_33
表示微积分;δη分别是第一次和第二次时域积分计算后产生的常数项。/>
Figure SMS_34
是相邻采样点之间的时间间隔。/>
Figure SMS_35
、/>
Figure SMS_36
分别为t-1、t+1时刻的速度信号,/>
Figure SMS_37
、/>
Figure SMS_38
分别为t-1、t+1时刻的加速度信号,/>
Figure SMS_39
为t-1时刻的位移信号;当t-1小于零时,ν t(t)和s t(t) 的值设置为0。可以看出,由于干扰项εt 2/2,δtη的影响,二次积分计算后的s(t) 可能会偏离真实值。尽管趋势项和DC分量的影响可以通过诸如多项式拟合和计算数学期望的去趋势技术来减少,但是这些干扰不能被完全消除,并且去噪过程可能导致信号能量的损失。
与时域积分相比,频域积分避免了复杂的去趋势计算,可以有效地减小累积误差的影响。频域积分法首先对加速度信号进行傅里叶变换,然后在频域进行积分运算,最后通过傅里叶逆变换将积分信号转换到时域。公式(11)示出了信号f(t)的傅立叶变换的结果
Figure SMS_40
,并且其傅立叶逆变换由公式(12)阐释。
Figure SMS_41
(11)
Figure SMS_42
(12)
其中,ω是频率分量,e是自然常数,j=
Figure SMS_43
。公式(13)和(14)分别列出了通过频域积分获得的速度信号ν f(t)和位移信号s f(t)的表达式。
Figure SMS_44
(13)
Figure SMS_45
(14)
其中,∆f是频率分辨率,A(ω)是加速度信号a(t)经过傅里叶变换后的频域信号。基于上述四个公式和傅里叶变换原理,可以看出传统频域积分有两个重要的缺点。首先,主频信号的能量可能会扩散到其他频率分量,即频谱泄漏问题,这可能会导致通过傅里叶逆变换得到的时域信号失真。其次,实测ABA信号通常是宽频信号,其中低频和高频噪声会导致不必要的误差,尤其是低频趋势项是信号失真的主要原因。幸运的是,本技术所开发的相似性测试可以精确地提取代表每个车轮旋转周期的ABA信号,从而最小化频谱泄漏的影响。针对噪声干扰问题,本技术提出了一种改进的自适应带通滤波的频域积分方法。其旨在构造不同截止频率的带通滤波器对各阶车轮多边形进行分析,这非常有助于对感兴趣阶次的车轮多边形信息进行针对性分析,以获得真实的积分信号。基于改进的频域积分获得的速度信号ν if(t)和位移信号s if(t)由公式(15)至(19)所示
Figure SMS_46
其中H(ω)是自适应带通滤波器。
Figure SMS_47
是频率分辨率,N是采样点数,t是时间,v是车辆行驶速度,ω是角频率;R是车轮名义滚动圆处半径;K是总阶数;A(ω)是加速度信号a(t)经过傅里叶变换后的频域信号;当分析第k阶车轮多边形时,f lk f uk 是带通滤波器的下截止频率和上截止频率。以f k 为中心的带通滤波器可以有效降低干扰,只提取与第k阶车轮多边形相关的ABA信息。带通滤波器的范围设置为10 Hz。注意到SP’代表车轮一个旋转周期对应的采样点数,故可以推导出关系式SP’=Fs×2πR/v。一旦基于所提出的改进频域积分方法获得了对应于车轮多边形的每一阶的位移信号,就可以根据标准EN 15610来定量计算每一阶的粗糙度水平。最后,统计所有保存的数据片段的检测结果以统计其对应的分布特征和数学期望,可更可靠地表达车轮多边形各阶的粗糙度特征。
车轮多边形检测框架试验验证
现场测试
在成都至兴文的一条运营线路上,对一列高速列车进行了长期跟踪试验。详细记录车轮多边形、ABA和速度信号。利用BBM车轮粗糙度测试仪装置测量了沿车轮圆周方向的车轮粗糙度,其空间采样频率为1mm。安装B&K-4520型三轴加速度计于轴箱上,采样频率为5000 Hz。运营速度由GPS模块以5 Hz的采样频率获得。所有测试信号由Somat eDAQ数据采集仪器同步采集,以便于后续的信号处理和分析。
基于本技术的多边形检测结果
图5中的(a)部分示出了实测ABA和速度信号。根据本技术所提出的车轮多边形检测方法,采用2~1024Hz的带通滤波器滤除零漂、DC分量和高频噪声。然后,使用前面介绍的信号预处理方式和平稳性测试来计算每个数据片段的CCC。之后,基于相似性测试提取代表每个车轮旋转周期的准稳态ABA信号,其示例结果在图5中的(b) 部分中示出。注意到在相似性测试的迭代算法中需要确定五个参数。本文研究的高速列车运行速度为200~250 km/h,因此ν minν max分别设为200和250。参考高速列车的维修手册,其初始车轮半径为460 mm,服役的车轮半径应大于425 mm,因此,将R minR max分别设为470和415。上述参数设置可使迭代算法精确地提取整周期ABA信号。
以一个车轮镟修周期内的车轮多边形为研究对象,研究了本技术的有效性。图6(a)、图6(b)、图6(c)展示了不同运行里程下的单阶车轮多边形的预测结果。箱线图显示了各阶车轮多边形的粗糙度分布特征。带方块的蓝线是数学期望值,即利用本技术预测到的粗糙度的平均值。带圆圈的红线表示实测粗糙度值。图6(a)显示,当运行里程为50 259 km时,24阶车轮多边形处于初始阶段,其粗糙度水平相对较低(0.40 dB re 1微米)。随着运营里程增加到81 694 km和105 172km时,24阶车轮多边形逐渐发展,分别达到8.43 dB和14.14 dB。可以观察到,预测值的主要特征与实测值非常接近,这表明所提出的方法可以有效地评估车轮多边形的严重程度。此外,表1给出定量的比较结果。显然,所提出的方法不仅能够以100%的准确度识别车轮多边形的主要阶次,而且能够以非常小的误差检测粗糙度水平。这证明所提出的方法可以高精度地检测初始车轮多边形,其绝对误差仅为1.18 dB。除了单阶车轮多边形之外,图7(a)、图7(b)、图7(c)描述了具有多个主导阶次的车轮多边形检测结果,以进一步说明所提出的方法的通用性。预测值和实测值的趋势非常相似,并且车轮多边形的主要特征几乎相同。这证明了该方法在处理不同行驶里程下的多阶车轮多边形时,仍能保持较强的检测能力,这对该方法在实际工程中的应用具有重要意义。此外,表2表明所提出的方法能够正确地识别车轮多边形的主阶。预测和实测粗糙度值之间的最大和最小绝对误差为-2.33 dB和0.04 dB,分别对应于运行里程为50 259 km和105 172 km的15阶和25阶车轮多边形。结合图6(a)、图6(b)、图6(c)、图7(a)、图7(b)、图7(c)可以看出,预测到的第10至30阶车轮多边形的粗糙度分布特征与实测值大致一致,表明所提出的方法也可以用于非主导车轮多边形特征的定性诊断,这有助于更全面地了解车轮多边形状态。图6(a)、图6(b)、图6(c)、图7(a)、图7(b)、图7(c)中,Dominant characteristics表示主特征,包括主阶次及其对应的粗糙度水平。
表1单个阶次主导的车轮多边形预测值与实测值的定量比较表
Figure SMS_48
表2多个阶次主导的车轮多边形预测值与实测值的定量比较表
Figure SMS_49
结论
本技术提出了一个新颖的数据驱动的检测框架,以便及时并准确地检测车轮多边形。具体地,设计了适当的信号预处理方法、平稳性测试和相似性测试,以精确提取代表每个车轮旋转周期的准静态ABA信号。在提取的ABA信号驱动下,基于惯性基准法和提出的改进频域积分方法,可以获得车轮多边形的定量表达,即阶次和粗糙度特征。采用一个车轮镟修周期内的现场跟踪测试评估了本技术的有效性。结果表明,提出的高速铁路车轮多边形检测框架能够定量检测不同运行里程下单阶和多阶车轮多边形的主要特征,最小和最大绝对误差分别为0.04 dB和-2.33 dB。
以上示意性的对本发明及其实施方式进行了描述,该描述没有限制性,附图中所示的也只是本发明的实施方式之一,实际的结构并不局限于此。所以,如果本领域的普通技术人员受其启示,在不脱离本发明创造宗旨的情况下,不经创造性的设计出与该技术方案相似的结构方式及实施例,均应属于本发明的保护范围。

Claims (5)

1.改进频域积分法驱动的高速列车车轮多边形定量检测方法,其特征在于:包括以下步骤:
1)从运行中的列车上采集轴箱加速度ABA和速度信号,并采用滑动窗口技术将数据划分为不同的片段;
2)引入稳态试验来计算每个片段的一致相关系数CCC,保存满足阈值要求的片段;
3)采用相似性试验从保存的片段中精确提取代表每个车轮旋转周期的准平稳ABA信号;
4)对提取的准平稳ABA信号进行改进的频域积分分析以定量检测车轮多边形。
2.根据权利要求1所述的改进频域积分法驱动的高速列车车轮多边形定量检测方法,其特征在于:1)中,具体包括以下步骤:
1.1)轴箱加速度ABA和速度信号采集;
1.2)进行带通滤波;
1.3)使用滑动窗将ABA信号切割成一系列数据片段(S 1, S 2, ..., S m , ..., S n ),其中m为其中某一数据片段,n为数据片段的总数量;
1.4)使用子滑动窗将切割后的数据片段再次细分为3个数据片段(S m1S m2S m3)。
3.根据权利要求2所述的改进频域积分法驱动的高速列车车轮多边形定量检测方法,其特征在于:2)中,具体包括以下步骤:
2.1)执行频谱和概率密度分布分析以计算数据片段(S m1S m2S m3)的频谱序列(F 1F 2F 3)和概率密度分布序列(P 1P 2P 3);
2.2)分别计算F 1F 2F 1F 3F 2F 3P 1P 2P 1P 3P 2P 3的CCC值C 1C 2C 3C 4C 5C 6,计算公式如下:
Figure QLYQS_1
Figure QLYQS_2
其中C表示一致相关系数CCC,C越接近1,序列AB之间的相关性越高,ρ为皮尔逊相关系数,cov表示求协方差,σ A σ B 分别为序列AB的标准差;
Figure QLYQS_3
分别为序列AB的向量;/>
Figure QLYQS_4
和/>
Figure QLYQS_5
分别为/>
Figure QLYQS_6
的均值;
通过一致相关系数CCC同时评估两个序列之间的相关性和误差水平,只要(F 1F 2F 3),(P 1P 2P 3)之间的相关性和误差水平满足阈值,就认为数据片段相对稳定;
2.3)计算CCC序列的平均值MCMC m = mean (C 1 C 2 C 3 C 4 C 5 C 6) ;数据片段的MC值被存储在临时组中,用于后续的相似性测试。
4.根据权利要求3所述的改进频域积分法驱动的高速列车车轮多边形定量检测方法,其特征在于:3)中,具体包括以下步骤:
3.1)将MC值大于阈值T的数据片段保存;
3.2)利用基于迭代算法驱动的相似性试验从保存的数据片段中精确提取代表车轮每个车轮周期的准平稳ABA信号,方法如下:
输入:最小速度V min,最大速度V max,最小半径R min,最大半径R max
采样频率Fs
输出:整旋转周期的信号序列y
Figure QLYQS_7
其中,T min是车轮旋转一周最短时间,T max 是车轮旋转一周最长时间,j是车轮旋转一周的信号采样点数,SP min是车轮旋转一周最小信号采样点数,SP max 是车轮旋转一周最大信号采样点数;迭代算法的输入是最小速度V min、最大速度V max、最小车轮半径R min、最大车轮半径R max和已知的采样频率Fs;输出y表示预期的整周期ABA信号;x 1x 2x 3是三个端到端的具有不同数据采样点的信号序列,ρ 1 ρ 2ρ 3是对应x 1x 2x 3的三个皮尔逊相关系数,type 表示类型pearson表示皮尔逊相关系数; y 1y 2y 3是对应x 1x 2x 3的三个输出;mean为mean函数,corr为corr函数,end表示结束,max(ρ)→SP’ 表示最大ρ对应SP’ ,初始化的功能是基于输入变量获得对应于一个车轮旋转的最小和最大采样点SP minSP max;迭代算法的思想是在平稳性测试后保存的数据片段S m中生成三个端到端的具有不同数据采样点的信号序列x 1x 2x 3,并计算三个序列之间的平均皮尔逊相关系数ρ,利用周期信号之间的高相似性,最大ρ对应的SP’正是一个车轮转动周期所对应的ABA信号的长度,由此能提取代表每个车轮旋转周期的ABA信号。
5.根据权利要求4所述的改进频域积分法驱动的高速列车车轮多边形定量检测方法,其特征在于:4)中,改进的频域积分分析采用一种改进的自适应带通滤波的频域积分方法,获得的速度信号ν if(t)和位移信号s if(t),公式为:
Figure QLYQS_8
;/>
Figure QLYQS_9
Figure QLYQS_10
Figure QLYQS_11
Figure QLYQS_12
其中H(ω)是自适应带通滤波器,
Figure QLYQS_13
是频率分辨率,A(ω)是加速度信号a(t)经过傅里叶变换后的频域信号,N是采样点数,当分析第k阶车轮多边形时,f lk f uk 是带通滤波器的下截止频率和上截止频率,以f k 为中心的带通滤波器能降低干扰,只提取与第k阶车轮多边形相关的ABA信息;t是时间,v是车辆行驶速度,ω是角频率;R是车轮名义滚动圆处半径;K是总阶数;e是指数函数中的底数;
基于改进的频域积分方法获得对应于车轮多边形的每一阶的位移信号,就能定量计算每一阶的粗糙度水平;
最后,统计所有保存的数据片段的检测结果以统计其对应的分布特征和数学期望。
CN202310536835.8A 2023-05-12 2023-05-12 改进频域积分法驱动的高速列车车轮多边形定量检测方法 Active CN116252820B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310536835.8A CN116252820B (zh) 2023-05-12 2023-05-12 改进频域积分法驱动的高速列车车轮多边形定量检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310536835.8A CN116252820B (zh) 2023-05-12 2023-05-12 改进频域积分法驱动的高速列车车轮多边形定量检测方法

Publications (2)

Publication Number Publication Date
CN116252820A true CN116252820A (zh) 2023-06-13
CN116252820B CN116252820B (zh) 2023-09-05

Family

ID=86688353

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310536835.8A Active CN116252820B (zh) 2023-05-12 2023-05-12 改进频域积分法驱动的高速列车车轮多边形定量检测方法

Country Status (1)

Country Link
CN (1) CN116252820B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116985865A (zh) * 2023-09-25 2023-11-03 成都运达科技股份有限公司 一种轨道交通车轮多边形故障诊断检测方法、装置及系统

Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003166808A (ja) * 2001-11-30 2003-06-13 Central Japan Railway Co 鉄道車両の相対変位測定装置
US20070176763A1 (en) * 2003-12-12 2007-08-02 Bridgestone Corporation Device and method for detecting abnormality of rotating body
US20070208841A1 (en) * 2006-03-01 2007-09-06 L-3 Communications Corporation Self-assembling wireless network, vehicle communications system, railroad wheel and bearing monitoring system and methods therefor
WO2009061291A1 (en) * 2007-11-05 2009-05-14 Square D Company Improvements in hierarchy determination for power monitoring systems
CN101886977A (zh) * 2010-06-12 2010-11-17 苏州大学 一种信号中周期瞬态成分的自适应检测方法
CN102590703A (zh) * 2012-02-18 2012-07-18 西南交通大学 基于零序暂态电荷的谐振接地系统单相接地故障选线方法
CN103471555A (zh) * 2013-09-30 2013-12-25 常州西南交大金路轨道交通科技发展有限公司 一种便携式轨道车辆车轮周向不平顺和直径测量装置
US20140095143A1 (en) * 2012-09-28 2014-04-03 International Business Machines Corporation Transliteration pair matching
CN105092467A (zh) * 2015-09-09 2015-11-25 南车工业研究院有限公司 一种高速铁路车轮踏面擦伤快速监测装置及方法
CN105424388A (zh) * 2015-11-17 2016-03-23 苏州大学 一种基于参数化多普勒瞬态模型的列车轮对轴承故障瞬态特征检测方法
KR20180110783A (ko) * 2017-03-30 2018-10-11 한국철도공사 다수 개의 매개변수에 의한 철도 차량의 진단 시스템 및 그에 의한 철도 차량의 안전 진단 방법
EP3637078A1 (en) * 2018-10-11 2020-04-15 Continental Automotive GmbH Method and system for localizing unbalance of a wheel of a vehicle, a system for localizing said unbalanced wheel, electronic control unit and storage medium
CN111207938A (zh) * 2020-01-13 2020-05-29 西南交通大学 一种铁道车辆故障检测方法
CN111695829A (zh) * 2020-06-17 2020-09-22 南京星邺汇捷网络科技有限公司 一种指标波动周期计算方法、装置、存储介质及电子设备
CN112381027A (zh) * 2020-11-23 2021-02-19 西南交通大学 基于列车轴箱垂向加速度信号的车轮多边形波深估计方法
US20210173045A1 (en) * 2015-07-17 2021-06-10 Yuqian HU Method, apparatus, and system for fall-down detection based on a wireless signal
CN113297439A (zh) * 2021-05-31 2021-08-24 中电福富信息科技有限公司 一种时间序列数据特征提取及分类方法
CN114544206A (zh) * 2022-01-20 2022-05-27 清华大学 一种轨道交通机车车辆轮对多边形故障检测方法及装置
CN114997252A (zh) * 2022-08-05 2022-09-02 西南交通大学 一种基于惯性原理的车轮多边形车载检测方法
CN115130495A (zh) * 2022-05-12 2022-09-30 哈尔滨工程大学 一种滚动轴承故障预测方法及系统
CN115600086A (zh) * 2022-11-15 2023-01-13 西南交通大学(Cn) 一种基于卷积回归的钢轨波磨粗糙度车载式定量检测方法

Patent Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003166808A (ja) * 2001-11-30 2003-06-13 Central Japan Railway Co 鉄道車両の相対変位測定装置
US20070176763A1 (en) * 2003-12-12 2007-08-02 Bridgestone Corporation Device and method for detecting abnormality of rotating body
US20070208841A1 (en) * 2006-03-01 2007-09-06 L-3 Communications Corporation Self-assembling wireless network, vehicle communications system, railroad wheel and bearing monitoring system and methods therefor
WO2009061291A1 (en) * 2007-11-05 2009-05-14 Square D Company Improvements in hierarchy determination for power monitoring systems
CN101886977A (zh) * 2010-06-12 2010-11-17 苏州大学 一种信号中周期瞬态成分的自适应检测方法
CN102590703A (zh) * 2012-02-18 2012-07-18 西南交通大学 基于零序暂态电荷的谐振接地系统单相接地故障选线方法
US20140095143A1 (en) * 2012-09-28 2014-04-03 International Business Machines Corporation Transliteration pair matching
CN103471555A (zh) * 2013-09-30 2013-12-25 常州西南交大金路轨道交通科技发展有限公司 一种便携式轨道车辆车轮周向不平顺和直径测量装置
US20210173045A1 (en) * 2015-07-17 2021-06-10 Yuqian HU Method, apparatus, and system for fall-down detection based on a wireless signal
CN105092467A (zh) * 2015-09-09 2015-11-25 南车工业研究院有限公司 一种高速铁路车轮踏面擦伤快速监测装置及方法
CN105424388A (zh) * 2015-11-17 2016-03-23 苏州大学 一种基于参数化多普勒瞬态模型的列车轮对轴承故障瞬态特征检测方法
KR20180110783A (ko) * 2017-03-30 2018-10-11 한국철도공사 다수 개의 매개변수에 의한 철도 차량의 진단 시스템 및 그에 의한 철도 차량의 안전 진단 방법
EP3637078A1 (en) * 2018-10-11 2020-04-15 Continental Automotive GmbH Method and system for localizing unbalance of a wheel of a vehicle, a system for localizing said unbalanced wheel, electronic control unit and storage medium
CN111207938A (zh) * 2020-01-13 2020-05-29 西南交通大学 一种铁道车辆故障检测方法
CN111695829A (zh) * 2020-06-17 2020-09-22 南京星邺汇捷网络科技有限公司 一种指标波动周期计算方法、装置、存储介质及电子设备
CN112381027A (zh) * 2020-11-23 2021-02-19 西南交通大学 基于列车轴箱垂向加速度信号的车轮多边形波深估计方法
CN113297439A (zh) * 2021-05-31 2021-08-24 中电福富信息科技有限公司 一种时间序列数据特征提取及分类方法
CN114544206A (zh) * 2022-01-20 2022-05-27 清华大学 一种轨道交通机车车辆轮对多边形故障检测方法及装置
CN115130495A (zh) * 2022-05-12 2022-09-30 哈尔滨工程大学 一种滚动轴承故障预测方法及系统
CN114997252A (zh) * 2022-08-05 2022-09-02 西南交通大学 一种基于惯性原理的车轮多边形车载检测方法
CN115600086A (zh) * 2022-11-15 2023-01-13 西南交通大学(Cn) 一种基于卷积回归的钢轨波磨粗糙度车载式定量检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
艾永军: "基于轮轨相对位移的列车脱轨监测算法研究", 中国优秀硕士学位论文全文数据库工程科技Ⅱ辑, no. 4, pages 4 *
郝俊虎;胡毅;: "基于XGBoost和自回归模型的轴承故障诊断和预警方法研究", 组合机床与自动化加工技术, no. 02, pages 145 - 147 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116985865A (zh) * 2023-09-25 2023-11-03 成都运达科技股份有限公司 一种轨道交通车轮多边形故障诊断检测方法、装置及系统
CN116985865B (zh) * 2023-09-25 2023-11-28 成都运达科技股份有限公司 一种轨道交通车轮多边形故障诊断检测方法、装置及系统

Also Published As

Publication number Publication date
CN116252820B (zh) 2023-09-05

Similar Documents

Publication Publication Date Title
Xu et al. High-speed train wheel set bearing fault diagnosis and prognostics: A new prognostic model based on extendable useful life
Alemi et al. Condition monitoring approaches for the detection of railway wheel defects
Molodova et al. Automatic detection of squats in railway infrastructure
Li et al. Fault detection method for railway wheel flat using an adaptive multiscale morphological filter
CN108515984B (zh) 一种车轮伤损检测方法及装置
Cao et al. Wheel-bearing fault diagnosis of trains using empirical wavelet transform
CN106274983B (zh) 一种识别轨道交通车轮轮辋裂纹故障的机理诊断方法
JP5382991B2 (ja) 軌道系交通システムの異常診断方法及び異常診断システム
Mosleh et al. Early wheel flat detection: An automatic data-driven wavelet-based approach for railways
CN105092467B (zh) 一种高速铁路车轮踏面擦伤快速监测装置及方法
EP3219574B1 (en) Method and system for determining a vertical profile of a rail surface
CN106018557A (zh) 一种基于小波分析的轨道表面凹陷检测方法
CN116252820B (zh) 改进频域积分法驱动的高速列车车轮多边形定量检测方法
JP7177028B2 (ja) レール波状摩耗の進展検知方法及び進展検知システム
CN104833535A (zh) 一种铁道车辆车轮踏面擦伤的检测方法
Hoelzl et al. On-board monitoring for smart assessment of railway infrastructure: A systematic review
CN113776760B (zh) 基于整轴振动分析的列车轮对失圆故障监测方法和系统
CN113988326A (zh) 一种地铁设备维修优化方法及系统
Nowakowski et al. Diagnostics of the drive shaft bearing based on vibrations in the high-frequency range as a part of the vehicle's self-diagnostic system
CN114997252A (zh) 一种基于惯性原理的车轮多边形车载检测方法
Nowakowski et al. Tram gearbox condition monitoring method based on trackside acoustic measurement
Zhang et al. Fault diagnosis of axle box bearing with acoustic signal based on chirplet transform and support vector machine
Nacib et al. A comparative study of various methods of gear faults diagnosis
Alten et al. Detecting and classifying rail corrugation based on axle bearing vibration
Xie et al. High-speed railway wheel polygon detection framework using improved frequency domain integration

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant