CN109406147B - 一种变速工况下的列车轴承轨边声学诊断方法 - Google Patents
一种变速工况下的列车轴承轨边声学诊断方法 Download PDFInfo
- Publication number
- CN109406147B CN109406147B CN201811270155.1A CN201811270155A CN109406147B CN 109406147 B CN109406147 B CN 109406147B CN 201811270155 A CN201811270155 A CN 201811270155A CN 109406147 B CN109406147 B CN 109406147B
- Authority
- CN
- China
- Prior art keywords
- train
- working condition
- fault
- speed
- neural network
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic or vibration analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M17/00—Testing of vehicles
- G01M17/08—Railway vehicles
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种变速工况下的列车轴承轨边声学诊断方法,包括:1)通过安装在铁轨两侧的麦克风采集列车在不同车速通过时列车轴承发出的声学信号并对其进行预处理;2)提取峰值因子、方差等7个特征并和车速Vs一起制作成具有8个特征量的历史样本集;3)利用历史样本集训练得到4层BP神经网络;4)采集当前列车轴承发出的轨边声学信号和列车车速,通过步骤3)的BP神经网络进行故障诊断。本发明消除了列车变速工况下多普勒声学信号畸变对故障特征的影响,实现了在变速工况下列车轴承故障的精准辨识;只需测量车速值,所需传感器数量更少;实现“不矫正前提下的精确诊断”,计算量更少;充分利用历史数据信息,诊断结果可靠性得到有效提升。
Description
技术领域
本发明属于高速列车轴承轨边声学故障诊断领域,具体涉及一种变速工况下的列车轴承轨边声学诊断方法。
背景技术
铁路作为重要的交通运输方式之一,由于运行成本低、运载量大等特点,其在货物运输及中短途旅客运输中承担着越来越重要的角色。而列车轴承作为列车的关键部件,其健康状况对列车的整体运行和乘客的生命安全有着重大的意义和影响。因此,对轴承进行在线故障检测具有非常重要的现实意义。轨边声学诊断技术使用安装在铁轨两侧的麦克风采集列车通过时轴承发出的声音信号,通过信号处理技术识别轴承的健康状况。由于列车的高速运行,轨边信号具有的多普勒畸变问题一直是困扰列车轴承轨边声学诊断的核心技术难题。现有方法是采用多普勒畸变矫正方法来消除多普勒畸变,但矫正处理需要预先获取车速、麦克风与轴承声源运动直线的垂直距离、横向距离等多个参数,并且需要进行插值处理,计算量较大。本发明旨在另辟蹊径,提出一种“不矫正前提下的精确诊断”,直接从采集信号中提取特征,使用神经网络来构建故障特征、车速与诊断结果之间的映射关系。
发明内容
本发明提供了一种变速工况下的列车轴承轨边声学诊断方法。本发明创新性地使用神经网络构建故障特征、车速Vs与诊断结果之间的关系,消除了多普勒畸变对传统故障特征的影响,实现了在变速工况下对轴承故障类型的精准判断。
本发明采用的技术方案为:一种变速工况下的列车轴承轨边声学诊断方法,包括以下5个步骤:
步骤1:历史轨边信号采集与预处理
通过麦克风阵列采集M种不同故障类型的轴承声源在N种不同速度下的声音信号,并将采集到的信号进行带通滤波处理,除去低频噪声,得到去噪后的共M*N个历史轨边信号;
步骤2:制作历史样本集
提取滤波后的轨边信号的传统特征,并将与之对应的列车车速Vs作为一个新的特征和传统特征一起制作成历史样本集。
步骤3:构建BP神经网络
创建一个4层的BP神经网络,包括一个输入层、两个隐藏层和一个输出层。根据经验和试探法搭建合适的隐层神经元个数,按实际需求给定学习率和迭代步数、选取合适激励函数、拟定最优优化算法与迭代步数;
步骤4:建立BP神经网络模型
用步骤2中的历史样本集训练步骤3中构建的网络,得到BP神经网络模型
步骤5:利用模型对当前类别未知的样本进行诊断
重新采集待检测的当前样本集,利用步骤4中训练的网络对当前样本集进行故障类别诊断。
其中,所述步骤2中历史数据集的制作,其具体为:
③将列车车速Vs作为一个新的特征加入,则矩阵X中的每个向量都可以提取到包含类车车速Vs在内的共8种特征参数;
④将8种特征进行标准化处理,使每种特征的特征值按比例放缩至0附近;再将故障类型标签作标签二值化处理:第1种故障类型对应输出为第i种故障类型输出为即第i个值为1,其余均为零。对特征值和故障类型标签进行数据预处理后,即完成了历史数据集的制作过程;
其中,所述步骤3中神经网络具体参数设置如下:
①输入层神经元的个数与输入特征种类数相同,两个隐层神经元个数分别为100、100,输出层神经元个数与故障类型数目相同;
②神经网络激励函数为S型函数,其表达式为:
其中,x为神经元接收到的总输入值,f(x)为通过激励函数处理后神经元的输出;
其中,所述步骤4中训练BP神经网络所用的方法为带有动量项的Adam算法,其算法公式为:
其中,μ为动量因子,0≤μ<1。gt为梯度,mt、nt分别是对梯度的一阶矩估计和二阶矩估计,是对mt、nt的校正。动量法降低了网络对于误差曲面局部细节的敏感性,有效地抑制网络陷入局部极小,促进超参数动态调整。
本发明与现有技术相比的优点在于:
(1)本发明消除了列车变速工况下多普勒声学信号畸变对故障特征的影响,实现了在变速工况下对列车轴承故障的精准辨识。
(2)本发明只需测量车速值,需要的传感器数量更少。
(3)本发明实现“不矫正前提下的精确诊断”,计算量更少。
(4)本发明充分利用历史数据信息,诊断结果可靠性得到有效提升。
附图说明
图1为本发明方法流程图;
图2(a)为峭度Ku随速度变化图;
图2(b)为峰值因子C随速度变化图;
图3为本发明提供的方法和传统方法对比实施流程图;
图4为传统特征作为输入时BP神经网络结构示意图;
图5为提出的将车速Vs和传统特征同时作为输入时BP神经网络结构示意图;
图6为输入特征有无车速精度对比图;
图7(a)为传统特征作为输入时网络输出融合特征的类内间效果图;
图7(b)为车速Vs和传统特征同时作为输入时网络输出融合特征的类内间效果图。
具体实施方式
下面结合附图和实际的实验数据对本发明作进一步的详细说明。
本实施案例包括正常、内圈故障程度0.014英寸、外圈故障程度0.014英寸、滚子故障程度0.014英寸4种故障类型。将这4种故障类型对应神经网络输出类型标签分别标记为y1,y2,y3,y4。
如图1所示,本发明的实施包括以下步骤:
步骤1:历史轨边信号采集与预处理
首先通过麦克风采集正常、内圈故障程度0.014英寸、外圈故障程度0.014英寸、滚子故障程度0.014英寸4种故障类型的轴承声源在7种不同速度下的声音信号,采样频率为20kHz,并将采集到的信号用巴特沃斯(Butterworth)滤波器滤波,带通频率范围设置为[1000,4000],除去低频噪声,最后得到每种类型147组、共588组历史畸变信号,且每组信号都通过光电传感器测得相应的列车车速Vs。
步骤2:制作历史样本集
对上述588组历史轨边信号进行特征提取,提取峰值因子C、方差Va、均方根值Rm、峭度Ku、波形因子S、脉冲因子I、裕度因子L等7种特征。如图2(a)和2(b)所示,同一故障类型在不同速度下峭度Ku、峰值因子C等特征均会发生不同程度改变,导致不同故障间的特征值交融到一起,给故障诊断带来一定的误诊率。
如图3所示,为了验证本发明提出的加入列车车速Vs来消除多普勒畸变对特征分类影响的有效性,本次制作了A和B两个历史样本集。其中历史样本集A中的输入特征只含有7种传统特征,历史样本集B中的输入特征除传统特征外,还加入与每次采集信号时相对应的列车车速Vs。
将提取的峰值因子C、方差Va、均方根值Rm、峭度Ku、波形因子S、脉冲因子I、裕度因子L等7种传统特征分别标记为x1~x7,对7种特征的特征值进行标准化处理,使每种特征特征值按比例放缩至0附近,将特征对应的故障类型标签y1,y2,y3,y4作标签二值化处理,将预处理后的特征值和标签一起作为历史样本集A;将每种特征值下对应的列车车速Vs标记为x8作为新的特征和7种传统特征一起作为神经网络的输入特征,并对这8种特征进行标准化处理,并和经过标签二值化处理后的标签一起作为历史样本集B。
步骤3:构建BP神经网络
创建一个包含两个隐藏层的4层BP神经网络,根据经验和试探法最后确定两个隐层神经元的个数分别为50和50,输入层神经元个数与输入特征种类数目相同,输出层神经元个数为4;激励函数选定为S型函数,其表达式为:
最后确定学习率为0.0001,优化算法为带有动量项的Adam算法,具体过程如下:
其中,μ为动量因子,0≤μ<1。gt为梯度,mt、nt分别是对梯度的一阶矩估计和二阶矩估计,是对mt、nt的校正。动量法降低了网络对于误差曲面局部细节的敏感性,有效地抑制网络陷入局部极小,促进超参数动态调整。
步骤4:建立BP神经网络模型
历史样本集A和历史样本集B训练步骤3中构建的神经网络,其过程分别如图4和图5所示。经过8万次迭代后,网络的均方误差收敛趋于平缓,得到训练好的BP神经网络模型A和模型B。
步骤5:利用模型对当前类别未知的样本进行诊断
采集测试信号样本集,采样频率为20kHz,并将采集到的信号用巴特沃斯(Butterworth)滤波器滤波,带通频率范围设置为[1000,4000],除去低频噪声,得到252组当前测试信号,且每组信号都通过光电传感器测得相应的列车车速Vs。
分别利用模型A和模型B对当前样本集进行类别诊断,诊断准确率对比如图6所示,可以看出加入列车车速Vs后,模型B相比于模型A诊断准确率得到了显著提升。
为了检验输入网络的特征是否具有代表性,实验最后建立LDA模型进行可视化分析,将n_components数值设为2,即把第二个隐藏层输出的高维融合特征降维至2维,通过得到的类内间图可以直观观测到输入特征是否具有分类代表性;图7(a)和图7(b)分别为模型A和模型B对当前样本诊断时得到的隐层融合特征的类内间距离图。由图可明显看出,加入列车车速Vs后,类内间融合特征更加集中,类与类之间融合特征更加分散。
综上所述,本发明提出的方法有效消除了多普勒畸变给故障特征带来的影响,使网络的输入特征更具有分类代表性,从而大大提升了诊断的准确率。
Claims (1)
1.一种变速工况下的列车轴承轨边声学诊断方法,其特征在于:包括以下步骤:
步骤1:历史轨边信号采集与预处理
通过麦克风阵列采集M种不同故障类型的轴承声源在N种不同速度下的声音信号,并将采集到的信号进行带通滤波处理,除去低频噪声,得到去噪后的共M*N个历史轨边信号;
步骤2:制作历史样本集
③将列车车速Vs作为一个新的特征加入,则矩阵X中的每个向量都可以提取到包含列车车速Vs在内的共8种特征参数;
④将8种特征值进行标准化处理,使每种特征的特征值按比例放缩至0附近;再将故障类型标签作标签二值化处理,则第1种故障类型对应输出为第i种故障类型输出为即第i个值为1,其余均为零,对特征值和故障类型标签进行数据预处理后,即完成了历史数据集的制作过程;
步骤3:构建BP神经网络
创建一个4层的BP神经网络,包括一个输入层、两个隐藏层和一个输出层,根据经验和试探法搭建合适的隐藏层神经元个数,按实际需求给定学习率和迭代步数、选取合适激励函数、拟定最优优化算法与迭代步数;神经网络具体参数设置如下:
①输入层神经元的个数与输入特征种类数相同,两个隐藏层的神经元个数均为100,输出层神经元个数与故障类型数目相同;
②神经网络激励函数为S型函数,其表达式为:
其中,x为神经元接收到的总输入值,f(x)为通过激励函数处理后神经元的输出;
步骤4:建立BP神经网络模型
用步骤2中的历史样本集训练步骤3中构建的网络,得到BP神经网络模型;
训练BP神经网络所用的方法为带有动量项的Adam算法,其算法公式为:
其中,μ,ν分别为一阶矩衰减系数、二阶矩衰减系数,且0≤μ<1,0<ν<1,μt是μ的t次方,表示第t次迭代时一阶矩的衰减系数,νt是ν的t次方,表示第t次迭代时二阶矩的衰减系数,gt为梯度,mt、nt分别是对梯度的一阶矩估计和二阶矩估计,是对mt、nt的校正,η为学习率,ε是常量,ε取10-8,Δθt是参数更新前后的差值,动量法降低了网络对于误差曲面局部细节的敏感性,有效地抑制网络陷入局部极小,促进超参数动态调整;
步骤5:利用模型对当前待检轴承进行故障诊断
采集待检测的轴承的轨边声学信号和当前车速,利用步骤4中训练的网络进行故障诊断;
该变速工况下的列车轴承轨边声学诊断方法消除了列车变速工况下多普勒声学信号畸变对故障特征的影响,实现了在变速工况下对列车轴承故障的精准辨识;
该变速工况下的列车轴承轨边声学诊断方法只需测量车速值,需要的传感器数量更少;
该变速工况下的列车轴承轨边声学诊断方法实现“不矫正前提下的精确诊断”,计算量更少;
该变速工况下的列车轴承轨边声学诊断方法充分利用历史数据信息,诊断结果可靠性得到有效提升。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811270155.1A CN109406147B (zh) | 2018-10-29 | 2018-10-29 | 一种变速工况下的列车轴承轨边声学诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811270155.1A CN109406147B (zh) | 2018-10-29 | 2018-10-29 | 一种变速工况下的列车轴承轨边声学诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109406147A CN109406147A (zh) | 2019-03-01 |
CN109406147B true CN109406147B (zh) | 2020-11-13 |
Family
ID=65469761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811270155.1A Active CN109406147B (zh) | 2018-10-29 | 2018-10-29 | 一种变速工况下的列车轴承轨边声学诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109406147B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110580915B (zh) * | 2019-09-17 | 2022-03-25 | 中北大学 | 基于可穿戴式设备的声源目标识别系统 |
CN110444225B (zh) * | 2019-09-17 | 2022-03-25 | 中北大学 | 基于特征融合网络的声源目标识别方法 |
CN113049251A (zh) * | 2021-03-16 | 2021-06-29 | 哈工大机器人(合肥)国际创新研究院 | 一种基于噪声的轴承故障诊断方法 |
CN113436645A (zh) * | 2021-06-09 | 2021-09-24 | 杭州电子科技大学 | 复杂噪声环境下的机电系统故障在线监测声学处理方法 |
CN114167837B (zh) * | 2021-12-02 | 2023-09-15 | 中国路桥工程有限责任公司 | 一种铁路信号系统的智能故障诊断方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN204556257U (zh) * | 2015-04-22 | 2015-08-12 | 北京康拓红外技术股份有限公司 | 一种铁路车辆滚动轴承声学诊断装置 |
CN106226078A (zh) * | 2016-07-01 | 2016-12-14 | 中国科学技术大学 | 一种基于麦克风阵列的多普勒畸变声学信号的校正方法 |
CN106872171A (zh) * | 2017-04-10 | 2017-06-20 | 中国科学技术大学 | 一种多普勒声学信号的自适应学习校正方法 |
CN108398267A (zh) * | 2018-02-27 | 2018-08-14 | 安徽大学 | 一种高速列车轨边运动参数自适应识别方法 |
CN108709744A (zh) * | 2018-03-05 | 2018-10-26 | 安徽大学 | 一种变负载工况下电机轴承故障诊断方法 |
-
2018
- 2018-10-29 CN CN201811270155.1A patent/CN109406147B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN204556257U (zh) * | 2015-04-22 | 2015-08-12 | 北京康拓红外技术股份有限公司 | 一种铁路车辆滚动轴承声学诊断装置 |
CN106226078A (zh) * | 2016-07-01 | 2016-12-14 | 中国科学技术大学 | 一种基于麦克风阵列的多普勒畸变声学信号的校正方法 |
CN106872171A (zh) * | 2017-04-10 | 2017-06-20 | 中国科学技术大学 | 一种多普勒声学信号的自适应学习校正方法 |
CN108398267A (zh) * | 2018-02-27 | 2018-08-14 | 安徽大学 | 一种高速列车轨边运动参数自适应识别方法 |
CN108709744A (zh) * | 2018-03-05 | 2018-10-26 | 安徽大学 | 一种变负载工况下电机轴承故障诊断方法 |
Non-Patent Citations (1)
Title |
---|
ADAM: A METHOD FOR STOCHASTIC OPTIMIZATION;Diederik P. Kingma et al;《ICLR 2015》;20151231;1-15 * |
Also Published As
Publication number | Publication date |
---|---|
CN109406147A (zh) | 2019-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109406147B (zh) | 一种变速工况下的列车轴承轨边声学诊断方法 | |
Shao et al. | Modified stacked autoencoder using adaptive Morlet wavelet for intelligent fault diagnosis of rotating machinery | |
CN108830127A (zh) | 一种基于深度卷积神经网络结构的旋转机械故障特征智能诊断方法 | |
CN112651167A (zh) | 一种基于图神经网络的半监督滚动轴承故障诊断方法 | |
CN110879254B (zh) | 一种基于改进最小二乘生成式对抗网络的钢轨裂纹声发射信号检测方法 | |
CN108256556A (zh) | 基于深度信念网络的风力发电机组齿轮箱故障诊断方法 | |
CN104833534A (zh) | 一种基于多源信息融合的列车走行系故障诊断装置及方法 | |
CN106314438A (zh) | 一种司机驾驶轨迹中异常轨迹的检测方法和系统 | |
CN108197731B (zh) | 一种基于共同训练的机动车尾气遥测与车检结果一致性方法 | |
CN104568444A (zh) | 变转速火车滚动轴承故障特征频率提取方法 | |
CN106198020A (zh) | 基于子空间和模糊c均值聚类的风电机组轴承故障诊断法 | |
CN106447040A (zh) | 基于异类多传感器数据融合的机械设备健康状态评估方法 | |
CN109446669B (zh) | 一种矿浆浓度的软测量方法 | |
CN110232415B (zh) | 一种基于生物信息特征的列车转向架故障识别方法 | |
CN104634872A (zh) | 高铁钢轨伤损在线监测方法 | |
CN112985830B (zh) | 一种abs结果自动判定算法 | |
CN112364706A (zh) | 一种基于类不平衡的小样本轴承故障诊断方法 | |
CN108267312A (zh) | 一种基于快速搜索算法的地铁列车轴承智能诊断方法 | |
CN113203464A (zh) | 一种动态车载称重系统的传感器故障检测方法 | |
CN105241665A (zh) | 一种基于IRBFNN-AdaBoost分类器的滚动轴承故障诊断方法 | |
CN113537152A (zh) | 一种基于深度神经网络的流场状态故障检测方法 | |
CN111461183A (zh) | 改进人工鱼群算法优化bp神经网络的平轮故障检测方法 | |
CN115586023A (zh) | 一种轨道车辆传动系统故障诊断方法及系统 | |
CN115371724A (zh) | 轨道状态评价方法及装置 | |
CN116644278A (zh) | 基于样本筛选和特征提取的双模型轴承异常检测算法 |
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 |