基于wav文件的柴油机齿轮异响识别方法与流程

1.本发明涉及一种柴油机故障检测方法，尤其涉及一种基于wav(wave form，即波形)文件的柴油机齿轮异响识别方法。


背景技术：

2.异响是柴油机的常见故障，例如：齿轮啮合异响、柴油机气门间隙过大或过小导致的异响、气缸垫的窜气异响等。其中，齿轮啮合异响的形成原因较复杂，各齿轮孔的形状和位置误差、齿轮偏心误差、齿形误差、齿距误差等均会造成齿轮啮合的周期性冲击和非周期性振动，从而产生周期性和非周期性的齿轮啮合异响。
3.现有技术中识别齿轮异响的主要方法是：采用西门子公司、bk公司、朗德公司等研发的专业测试与分析软件，但这些软件需要基于nvh的理论知识、通过复杂的fft、阶次、colormap等技术和经验来完成，得到的信息量较多，需通过专业人员进行分析和识别，无法得出直观的柴油机好或坏的故障排查结论，因此在线识别方面很少有运用。
4.中国发明专利zl201710343372.8公开了一种柴油机发电机组齿轮异响的故障排查方法，并具体公开了：步骤一，对发电机机组负载进行故障排查，计算运行参数的平稳波动；步骤二，对调控系统进行故障排查，在额定负载下调节pi参数，获取噪声不均匀度；步骤三，对增压系统和配气机构进行故障排查，测量进排气温度；步骤四，对喷油器进行故障排查，拆检喷油器。该专利是用于柴油发电机组的齿轮异响故障排查，不是单独的柴油机，且该故障排查方法的前提是柴油机发电机组的初始状态无故障，且在使用后出现了较明显的故障，该方法的应用场景特殊，不具有通用性，也无法用于处于生产阶段的单独的柴油机齿轮异响故障排查。


技术实现要素：

5.本发明的目的在于提供一种基于wav文件的柴油机齿轮异响识别方法，能在柴油机生产时实现在线产品合格情况的实时排查，确保加柴油机的生产节拍不受影响。
6.本发明是这样实现的：
7.一种基于wav文件的柴油机齿轮异响识别方法，包括以下步骤：
8.步骤1：在柴油机的上方设置噪声传感器，噪声传感器采集柴油机的噪声信号并生成wav噪声音频文件；
9.步骤2：噪声传感器将wav噪声音频文件传输至信号记录仪，信号记录仪将wav噪声音频文件传输至计算机；
10.步骤3：采集n个柴油机样品的wav噪声音频文件，并记录n个wav噪声音频文件对应的齿轮异响故障状态；
11.步骤4：基于n个wav噪声音频文件建立齿轮异响识别模型，齿轮异响识别模型的输入为wav噪声音频文件，齿轮异响识别模型的输出为判断值y；
12.步骤5：将柴油机的任意wav噪声音频文件输入齿轮异响识别模型，根据判断值y判
断柴油机是否合格。
13.所述的步骤4包括：
14.步骤4.1：分别对每一个wav噪声音频文件进行快速傅里叶变换，得到模型参数xi1、xi5、xi6；
15.其中，xi1为fft变换后的频率峰度值，i∈[1,n]；
[0016]
xi5为fft变换后正态分布的频域偏度值，i∈[1,n]；
[0017]
xi6为fft变换后正态分布的频域峰度值，i∈[1,n]；
[0018]
步骤4.2：对快速傅里叶变换后的信号进行归一化处理，得到模型参数xi2、xi3、xi4。
[0019]
其中，xi2为信号归一化后正态分布的时域偏度值，i∈[1,n]；
[0020]
xi3为信号归一化后正态分布的时域峰度值，i∈[1,n]；
[0021]
xi4为信号归一化后柯西分布的cdf，i∈[1,n]；
[0022]
x1＝xi1，x2＝xi2，x3＝xi3，x4＝xi4，x5＝xi5，x6＝xi6，i∈[1,n]
ꢀꢀꢀꢀ
方程(1)
[0023]
步骤4.3：建立参数方程：
[0024]
y＝xβ ε
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
方程(2)
[0025][0026]
柴油机合格对应的yi＝1，柴油机不合格对应的yi＝0，i∈[1,n]
ꢀꢀꢀꢀꢀ
方程(4)
[0027]
xi7＝xi1^2，xi8＝xi1*xi2，xi9＝xi1*xi3，
…
，xi12＝xi1*xi6，xi13＝xi2*xi2，
…
，xi17＝xi2*xi6，xi18＝xi3*xi3，
…
，xi27＝xi6*xi6
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
方程(5)；
[0028]
步骤4.4：将方程(1)、(5)和(4)带入方程(3)中，计算β和ε的值，并根据方程(2)得到模型关系式：
[0029]
y＝(1937786915459869*x1)/9007199254740992 (3676422554321869*x2)/17592186044416 (301004670484087*x3)/35184372088832 (5453649764968379*x4)/70368744177664
‑
(4810986767505429*x5)/1125899906842624 (7599815182024087*x6)/18014398509481984 (267752683928105*x1*x2)/70368744177664
‑
(5436930222024705*x1*x3)/576460752303423488
‑
(2406166496145273*x1*x4)/4503599627370496
‑
(8620055050281217*x2*x3)/562949953421312 (4149610817155857*x1*x5)/72057594037927936
‑
(5061598959295695*x2*x4)/35184372088832 (872003768366459*x1*x6)/36028797018963968
‑
(7373407150549555*x2*x5)/140737488355328
‑
(5016021620551169*x3*x4)/1125899906842624
‑
(3674189446593231*x2*x6)/562949953421312
‑
(5036524658204117*x3*x5)/4503599627370496
‑
(6101276554733531*x3*x6)/72057594037927936
‑
(1257559061559625*x4*x5)/70368744177664 (349969757861721*x4*x6)/562949953421312
‑
(7914723802440157*x5*x6)/18014398509481984
‑
(8150881690082041*x1*x1)/1152921504606846976
(5330672480157831*x2*x2)/140737488355328
‑
(6840852579128251*x3*x3)/18014398509481984
‑
(2315181900228411*x4*x4)/281474976710656 (4109815306072147*x5*x5)/562949953421312
‑
(7168538800631285*x6*x6)/1152921504606846976
‑
3080342811025019/70368744177664。
[0030]
所述的步骤5中，柴油机是否合格的判断方法是：当
‑
0.582≤y≤0.513时，柴油机合格；当y＜
‑
0.582或y＞0.513时，柴油机不合格。
[0031]
所述的柴油机合格是指柴油机齿轮啮合异响故障概率小于50％，柴油机不合格是指柴油机齿轮啮合异响故障概率大于等于50％。
[0032]
在计算机的软件上设置红色灯和绿色灯，当判定为所述的柴油机合格时，绿色灯亮，当判定为所述的柴油机不合格时，红色灯亮。
[0033]
所述的步骤3中，柴油机样品的数量n≥28，且n个wav噪声音频文件的长度均为30s。
[0034]
本发明与现有技术相比，具有如下有益效果：
[0035]
1、本发明由于基于wav噪声音频文件建立模型，通过模型参数xi1～xi6提高模型关系式的准确性和适应性，能在柴油机生产过程中进行在线噪声数据采集并通过模型学习完成在线质量实时排查，减少在柴油机生产下线后的质检步骤；同时，wav噪声音频文件可通过设置在柴油机自由端上方的噪声传感器直接采集，数据采集方便、快捷且针对性强。
[0036]
2、本发明由于建立了输入为wav噪声音频文件、输出为判断值的齿轮异响识别模型，能根据在线采集的噪声音频文件直接输出判断值用于判断柴油机是否存在齿轮异响故障，判断准确性高，不影响柴油机的生产节拍，确保生产效率和质量。
[0037]
本发明通过对wav噪声音频文件的模型学习，能在柴油机生产时实现在线产品合格情况的实时排查，数据采集和处理高效、准确，识别结果准确性高，确保柴油机的生产节拍不受影响。
附图说明
[0038]
图1是本发明基于wav文件的柴油机齿轮异响识别方法的流程图；
[0039]
图2是本发明基于wav文件的柴油机齿轮异响识别方法采用的硬件结构图。
[0040]
图中，1计算机，2信号记录仪，3柴油机，4噪声传感器，5信号线，6lan数据传输线。
具体实施方式
[0041]
下面结合附图和具体实施例对本发明作进一步说明。
[0042]
请参见附图1，一种基于wav文件的柴油机齿轮异响识别方法，包括以下步骤：
[0043]
步骤1：请参见附图2，在柴油机3的上方设置噪声传感器4，噪声传感器4采集柴油机3的噪声信号并生成wav噪声音频文件。优选的，噪声传感器4的采样频率为2khz，噪声传感器4的采样时长由柴油机试验台架的控制程序确定。
[0044]
所述的噪声传感器4设置在柴油机3的自由端的上方，优选的，噪声传感器4距离柴油机3的自由端顶面20cm。由于柴油机3不同部位的噪声信号不同，柴油机3自由端处的噪声信号主要来自柴油机发电机组的齿轮，有利于对齿轮啮合异响进行故障的精确诊断。
[0045]
步骤2：噪声传感器4将wav噪声音频文件通过信号线5传输至信号记录仪2，信号记
录仪2将wav噪声音频文件通过lan数据传输线6传输至计算机1。信号记录仪2用于识别并记录wav格式的文件，确保传输至计算机1的文件均为wav噪声音频文件，避免其他格式信号的干扰。
[0046]
步骤3：采集n个柴油机样品的wav噪声音频文件，并记录n个wav噪声音频文件对应的齿轮异响故障状态。柴油机样品的数量n≥28；n个wav噪声音频文件的音频长度均为30s左右，确保建立的模型具有较高的准确性，且兼顾了数据处理量。
[0047]
步骤4：基于n个wav噪声音频文件建立齿轮异响识别模型，齿轮异响识别模型的输入为wav噪声音频文件，齿轮异响识别模型的输出为判断值y。
[0048]
步骤4.1：分别对每一个wav噪声音频文件进行快速傅里叶变换(fft)，得到模型参数xi1、xi5、xi6，快速傅里叶变换公式为：
[0049][0050]
其中，k＝0，1，
…
，n
‑
1，i∈[1,n]。
[0051]
其中，xi1为fft变换后的频率峰度值，i∈[1,n]。峰度值用于研究数据分布陡峭或者平滑的统计量，峰度值是噪声运行点的主要特征之一，频率峰度值不同说明噪声频率的平滑程度不同。
[0052]
xi5为fft变换后正态分布的频域偏度值，i∈[1,n]。偏度值用于衡量统计数据的不对称性，频域偏度值是噪声运行点的主要特征之一，相同情况下频域偏度值不同说明噪声不对称性不同。
[0053]
xi6为fft变换后正态分布的频域峰度值，i∈[1,n]。峰度值用于研究数据分布陡峭或者平滑的统计量，频域峰度值是噪声运行点的主要特征之一，相同情况下频域峰度值不同说明噪声频率分布的平滑程度不同。
[0054]
步骤4.2：分别对每一个wav噪声音频文件进行归一化处理，得到模型参数：xi2、xi3、xi4。归一化转换函数为：x*＝(x
‑
min)/(max
‑
min)，用于对wav噪声音频文件进行线性变换，使其映射到[0,1]区间。
[0055]
其中，xi2为信号归一化后正态分布的时域偏度值，i∈[1,n]。偏度用于衡量统计数据的不对称性，时域偏度值是噪声运行点的主要特征之一，相同情况下时域偏度值不同说明噪声不对称性状态不同。
[0056]
xi3为信号归一化后正态分布的时域峰度值，i∈[1,n]。峰度值用于研究数据分布陡峭或者平滑的统计量，时域峰度值是噪声运行点的主要特征之一，相同情况下时域峰度值不同说明噪声随转速变化的平滑影响不同。
[0057]
xi4为信号归一化后柯西分布的cdf(cumulative distribution function，累计分布函数)，i∈[1,n]。柯西分布是一个数学期望不存在的连续型概率分布，方差和高阶矩均不存在，可以描述受迫共振的微分方程的解，用于表征噪声的变化率。
[0058]
快速傅里叶变换和归一化处理是信号处理的常用处理方式，此处不再赘述。
[0059]
x1＝xi1，x2＝xi2，x3＝xi3，x4＝xi4，x5＝xi5，x6＝xi6，i∈[1,n]
ꢀꢀꢀꢀ
方程(1)
[0060]
1与β0相乘还是β0，表示线性方程组中的常数部分。
[0061]
步骤4.3：建立参数方程。
[0062]
y＝xβ ε
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
方程(2)
[0063][0064]
柴油机合格对应的yi＝1，柴油机不合格对应的yi＝0，i∈[1,n]
ꢀꢀꢀꢀꢀꢀ
方程(4)。
[0065]
通过近似算法拟合，得到：
[0066]
xi7＝xi1^2，xi8＝xi1*xi2，xi9＝xi1*xi3，
…
，xi12＝xi1*xi6，xi13＝xi2*xi2，
…
，xi17＝xi2*xi6，xi18＝xi3*xi3，
…
，xi27＝xi6*xi6
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
方程(5)。
[0067]
步骤4.4：将方程(1)、(5)和(4)带入方程(3)中，计算β和ε的值，并根据方程(2)得到模型关系式：
[0068]
y＝(1937786915459869*x1)/9007199254740992 (3676422554321869*x2)/17592186044416 (301004670484087*x3)/35184372088832 (5453649764968379*x4)/70368744177664
‑
(4810986767505429*x5)/1125899906842624 (7599815182024087*x6)/18014398509481984 (267752683928105*x1*x2)/70368744177664
‑
(5436930222024705*x1*x3)/576460752303423488
‑
(2406166496145273*x1*x4)/4503599627370496
‑
(8620055050281217*x2*x3)/562949953421312 (4149610817155857*x1*x5)/72057594037927936
‑
(5061598959295695*x2*x4)/35184372088832 (872003768366459*x1*x6)/36028797018963968
‑
(7373407150549555*x2*x5)/140737488355328
‑
(5016021620551169*x3*x4)/1125899906842624
‑
(3674189446593231*x2*x6)/562949953421312
‑
(5036524658204117*x3*x5)/4503599627370496
‑
(6101276554733531*x3*x6)/72057594037927936
‑
(1257559061559625*x4*x5)/70368744177664 (349969757861721*x4*x6)/562949953421312
‑
(7914723802440157*x5*x6)/18014398509481984
‑
(8150881690082041*x1*x1)/1152921504606846976 (5330672480157831*x2*x2)/140737488355328
‑
(6840852579128251*x3*x3)/18014398509481984
‑
(2315181900228411*x4*x4)/281474976710656 (4109815306072147*x5*x5)/562949953421312
‑
(7168538800631285*x6*x6)/1152921504606846976
‑
3080342811025019/70368744177664。步骤5：将柴油机3的任意wav噪声音频文件输入齿轮异响识别模型，根据判断值y判断柴油机是否合格。
[0069]
所述的柴油机是否合格的判断方法是：当
‑
0.582≤y≤0.513时，柴油机合格；当y＜
‑
0.582或y＞0.513时，柴油机不合格。
[0070]
所述的柴油机合格是指柴油机齿轮啮合异响故障概率小于50％，柴油机不合格是指柴油机齿轮啮合异响故障概率大于等于50％。
[0071]
在计算机1的软件上设置红色灯和绿色灯，当判定为柴油机合格时，绿色灯亮，当判定为柴油机不合格时，红色灯亮。
[0072]
本发明以xi1、xi2、xi3、xi4、xi5、xi6为模型参数，基于28个柴油机样品建立模型
关系式，该模型参数对柴油机齿轮异响识别的在线排查准确率能达到89％以上。可通过增加柴油机样品的数量进一步修正模型关系式，提高模型关系式的适应性和准确性，从而提高柴油机齿轮异响识别的在线排查准确率。
[0073]
实施例1：
[0074]
对型号为d20的柴油机是否合格进行在线识别，本实施例中采用的噪声传感器4型号为bk 1/2”microphone type 4189(serial no:3023196)，噪声传感器4安装在柴油机3的自由端顶面上方20cm处，信号记录仪2的型号为zsy2048。
[0075]
通过噪声传感器4分别采集n＝28个柴油机3(标记为1#
‑
28#)的wav噪声音频文件，并分别记录该28个柴油机3是否合格。
[0076]
对28个wav噪声音频文件分别进行fft变换，得到模型参数xi1、xi5、xi6，对快速傅里叶变换后的信号进行归一化处理，得到模型参数xi2、xi3、xi4，i∈[1，28]。
[0077]
根据方程(1)至(5)计算得到：
[0078]
x1＝43.4219155160640，x2＝0.0372748012972，x3＝4.8225703199080，x4＝0.5265644078846，x5＝1.491275626，x6＝15.83873625。
[0079]
β0＝
‑
3080342811025019/70368744177664＝
‑
43.7743041604935。
[0080]
β1＝1937786915459869/9007199254740992＝0.215137565036091。
[0081]
β2＝3676422554321869/17592186044416＝208.98042716464。
[0082]
……
[0083]
β
27
＝
‑
7168538800631285/1152921504606846976＝
‑
0.00621771627295286。
[0084]
ε＝0.02937962032596890。
[0085]
从而得到模型关系式：
[0086]
y＝(1937786915459869*x1)/9007199254740992 (3676422554321869*x2)/17592186044416 (301004670484087*x3)/35184372088832 (5453649764968379*x4)/70368744177664
‑
(4810986767505429*x5)/1125899906842624 (7599815182024087*x6)/18014398509481984 (267752683928105*x1*x2)/70368744177664
‑
(5436930222024705*x1*x3)/576460752303423488
‑
(2406166496145273*x1*x4)/4503599627370496
‑
(8620055050281217*x2*x3)/562949953421312 (4149610817155857*x1*x5)/72057594037927936
‑
(5061598959295695*x2*x4)/35184372088832 (872003768366459*x1*x6)/36028797018963968
‑
(7373407150549555*x2*x5)/140737488355328
‑
(5016021620551169*x3*x4)/1125899906842624
‑
(3674189446593231*x2*x6)/562949953421312
‑
(5036524658204117*x3*x5)/4503599627370496
‑
(6101276554733531*x3*x6)/72057594037927936
‑
(1257559061559625*x4*x5)/70368744177664 (349969757861721*x4*x6)/562949953421312
‑
(7914723802440157*x5*x6)/18014398509481984
‑
(8150881690082041*x1*x1)/1152921504606846976 (5330672480157831*x2*x2)/140737488355328
‑
(6840852579128251*x3*x3)/18014398509481984
‑
(2315181900228411*x4*x4)/281474976710656 (4109815306072147*x5*x5)/562949953421312
‑
(7168538800631285*x6*x6)/1152921504606846976
‑
3080342811025019/70368744177664。
[0087]
在线生产时，将29#柴油机3的wav噪声音频文件输入齿轮异响识别模型，齿轮异响
识别模型输出判断值y＝0.970621＞0.513，29#柴油机3判定为不合格，即29#柴油机3存在齿轮异响故障。采用本发明的柴油机齿轮异响识别方法对后续生产的柴油机进行在线故障识别，故障排查准确率达到89％。
[0088]
以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围，因此，凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。