一种基于阶次边带乘积谱的齿轮健康监测指标构造方法技术

一种基于阶次边带乘积谱的齿轮健康监测指标构造方法，先使用计算阶次跟踪技术对测得振动加速度信号进行角度域重采样，并对重采样信号进行傅里叶变换得到阶次谱；再根据阶次分辨率和齿数信息，使用阶次谱计算阶次边带乘积谱，并以轮齿故障特征阶次为中心设计一个矩形窗，在阶次边带乘积谱中抓取对应轮齿故障信息；最后进行健康监测指标计算，并依托基于数据驱动的统计学习方法建立报警阈值自学习方法；本发明专利技术通过阶次边带乘积谱技术有效地抑制了带内噪声对监测指标构造的影响，所提出的报警阈值自学习方法简洁通用，避免了对专业技术人员的依赖，实现了齿轮传动实时健康监测及其阈值自设定功能。阈值自设定功能。阈值自设定功能。

【技术实现步骤摘要】
一种基于阶次边带乘积谱的齿轮健康监测指标构造方法


[0001]本专利技术属于机械故障诊断
，具体涉及一种基于阶次边带乘积谱的齿轮健康监测指标构造方法。

技术介绍

[0002]齿轮传动发生齿面失效时，其故障特征信号在时域中主要表现为幅值调制信号，在频域中则表现为以各阶啮合频率为中心，失效轮齿故障特征频率为边频带带宽的中心—边频形式，监测齿轮传动各阶啮合频率及其边带幅值变化是齿轮运行健康状态监测的基本思路。实际测得齿轮传动振动加速度信号中，总是不可避免的存在背景噪声干扰的问题，常用的窄带滤波方式虽能消除频带范围外的背景噪声，提纯故障信息，但无法去除带内噪声对啮合频率及其边带幅值的影响。故而在背景噪声干扰的情况下，对利用测得振动加速度信号监测齿轮运行健康状态提出了较高的要求。
[0003]计算阶次跟踪技术通过对测得振动加速度信号进行角度域等角度间隔重采样，使原先时域中因转速波动而产生的非平稳信号，转化为角度域平稳或循环平稳信号。故障设备响应在振动加速度信号频域中产生的故障特征频率，在角度域信号阶次谱中不会因转速波动或变化而产生频率模糊或者频率偏移，其始终以固定的阶次出现在阶次谱当中，这一特点使其在监测某一特定设备故障响应变化方面具备独有的优势。
[0004]乘积谱技术最早源于语音信号处理领域，因自然界中的语音信号通常都是由基础频率及其各阶谐波分量组成，乘积谱技术通过对频带内各分量按整数倍作乘，从而起到增强周期性谐波分量并抑制非周期性噪声分量的作用。机械系统故障响应与语音信号类似，其诱发的故障响应在振动加速度信号中同样表现为故障特征频率及其各阶谐波分量。通常使用窄带滤波技术，便能够十分清晰的将故障特征频率及其各阶谐波分量展示在信号频谱或包络谱中，但当滤波信号中存在带内噪声时，故障特征频率及其各阶谐波分量极易被带内噪声抑制，乘积谱技术增强周期性分量抑制非周期性分量的特点，使其在应对带内噪声时具备广阔的应用前景。
[0005]机械设备运行健康监测指标阈值的设定，一直以来都是机械故障诊断领域的关键问题，一个恰当的报警阈值设定是对所构造健康监测指标的有力补充。传统的阈值设定方法包括基于人为经验的设定方法、基于国家标准的设定方法以及数学推导方法三大类；基于人为经验的设定方法过于依赖工程技术人员的专业素养，而基于国家标准的设定方法和数学推导方法仅对特定指标有效，因此非常有必要开发一种不依赖人工经验并且通用性强的健康监测指标阈值设定方法。

技术实现思路

[0006]为了克服上述现有技术的缺点，本专利技术的目的在于提供了一种基于阶次边带乘积谱的齿轮健康监测指标构造方法，可在存在带内噪声的情况下监测齿轮传动的运行健康状态，并基于数据驱动的统计学习方法，实现了构造指标报警阈值的自学习自设定功能。
[0007]为了达上述目的，本专利技术采取的技术方案为：
[0008]一种基于阶次边带乘积谱的齿轮健康监测指标构造方法，包括以下步骤：
[0009]步骤一，采集得到齿轮传动振动加速度信号序列v(n)以及齿轮传动轴的转速脉冲信号序列p(n)；
[0010]步骤二，通过设定脉冲阈值T，提取转速脉冲信号序列中的脉冲沿所对应的时刻，脉冲沿计存方式包括上升沿和下降沿两种，具体计算公式如下：
[0011][0012][0013]上述公式中，tp
upper
(n)和tp
lower
(n)分别代表计算得到的上升沿和下降沿脉冲时刻序列，Fs表示信号采样频率，N为转速脉冲信号序列的总点数，并利用计存的脉冲时刻序列对振动加速度信号序列v(n)进行键相截齐，得到经键相截齐后的振动加速度信号序列v
c
(n)；
[0014]步骤三，根据脉冲时刻序列的总点数计算脉冲时刻转轴转过的弧度，Ap＝L
×
2π/ppr，其中Ap代表脉冲时刻内转轴转过的弧度，L为脉冲时刻序列总点数，ppr代表每转一圈所产生的脉冲个数；在[0,Ap]区间内，按照Δθ进行等角度间隔划分，生成重采样弧度序列AngArr(n)，其中Δθ＝π/order，order＝2
×
km
×
om，order代表最大可分辨阶次，km和om分别代表待计算啮合阶次阶数和啮合阶次；
[0015]步骤四，对脉冲时刻序列进行归零化处理，即对脉冲时刻序列中所有时刻对应的值减去序列中第一个时刻对应的值，并根据脉冲时刻序列总点数L计算脉冲时刻对应的弧度序列，即脉冲弧度序列AngP(n)，AngP(n)＝(0：1：L
‑
1)
×
2π/ppr；
[0016]步骤五，使用线性内插技术，根据已知的脉冲时刻序列tp(n)、脉冲弧度序列AngP(n)和重采样弧度序列AngArr(n)，按下述原理公式计算重采样时刻序列tang(n)，
[0017][0018]式中Num表示插值后重采样时刻序列的总点数，tang(1)＝tp(1),tang(Num)＝tp(Num)；
[0019]步骤六，根据步骤二得到的振动加速度信号序列vc(n)，利用采样频率Fs计算对应的时刻序列tc(n),tc(n)＝(0：1：length(vc(n))
‑
1)/Fs，式中length(
·
)表示求长度运算符，进而利用vc(n)、tc(n)和tang(n)根据步骤五中的线性内插技术原理求得等角度间隔重采样之后的重采样信号序列va(n)；
[0020]步骤七，对步骤六中得到的重采样信号序列va(n)，根据下述公式执行离散傅里叶变换，
[0021]式中O(k)即为所求阶次谱序列；
[0022]步骤八，从阶次谱序列当中提取各阶啮合阶次及两侧边带，进行加权求和并翻转；
[0023]步骤九，利用步骤八得到的单边加权和序列进行乘积谱运算，即阶次边带乘积谱技术，计算方法如下：
[0024]式中OSPS(k)即为所求阶次边带乘积谱；
[0025]步骤十，使用步骤九求得的阶次边带乘积谱构造齿轮健康监测指标，其计算公式如下：
[0026]CI＝log2(1+sum(OSPS
win
))，
[0027][0028]上述公式中，CI表示构造的齿轮健康监测指标，OSPS
win
表示阶次边带乘积谱中使用矩形窗截取的一段，α表示阶次系数，O
f
表示待监测故障阶次，sum(
·
)表示累加求和运算符；
[0029]步骤十一，使用基于数据驱动的统计学习方法，设定齿轮健康监测指标报警阈值。
[0030]所述的步骤八的具体过程为：
[0031]8.1)根据乘积谱乘积阶次和阶次分辨率计算单侧边带长度Ls，
[0032][0033]上式中，K表示乘积谱的乘积阶次，O
max
为最大阶次即角度采样率；
[0034]8.2)以各阶啮合阶次为中心，单侧边带长度为上下限，截取阶次谱序列中各阶啮合本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于阶次边带乘积谱的齿轮健康监测指标构造方法，其特征在于，包括以下步骤：步骤一，采集得到齿轮传动振动加速度信号序列v(n)以及齿轮传动轴的转速脉冲信号序列p(n)；步骤二，通过设定脉冲阈值T，提取转速脉冲信号序列中的脉冲沿所对应的时刻，脉冲沿计存方式包括上升沿和下降沿两种，具体计算公式如下：tp
upper
(n)＝(i+1)/Fss.t.p(i)＜T&&p(i+1)＞＝T，i＝1，2，3，
…
，N
’
tp
lower
(n)＝(i+1)/Fss.t.p(i)＞T&&p(i+1)＜＝T，i＝1，2，3，
…
，N
°
上述公式中，tp
upper
(n)和tp
lower
(n)分别代表计算得到的上升沿和下降沿脉冲时刻序列，Fs表示信号采样频率，N为转速脉冲信号序列的总点数，并利用计存的脉冲时刻序列对振动加速度信号序列v(n)进行键相截齐，得到经键相截齐后的振动加速度信号序列v
c
(n)；步骤三，根据脉冲时刻序列的总点数计算脉冲时刻转轴转过的弧度，Ap＝L
×
2π/ppr，其中Ap代表脉冲时刻内转轴转过的弧度，L为脉冲时刻序列总点数，ppr代表每转一圈所产生的脉冲个数；在[0，Ap]区间内，按照Δθ进行等角度间隔划分，生成重采样弧度序列AngArr(n)，其中Δθ＝π/order，order＝2
×
km
×
om，order代表最大可分辨阶次，km和om分别代表待计算啮合阶次阶数和啮合阶次；步骤四，对脉冲时刻序列进行归零化处理，即对脉冲时刻序列中所有时刻对应的值减去序列中第一个时刻对应的值，并根据脉冲时刻序列总点数L计算脉冲时刻对应的弧度序列，即脉冲弧度序列AngP(n)，AngP(n)＝(0∶1：L
‑
1)
×
2π/ppr；步骤五，使用线性内插技术，根据已知的脉冲时刻序列tp(n)、脉冲弧度序列AngP(n)和重采样弧度序列AngArr(n)，按下述原理公式计算重采样时刻序列tang(n)，式中Num表示插值后重采样时刻序列的总点数，tang(1)＝tp(1)，tang(Num)＝tp(Num)；步骤六，根据步骤二得到的振动加速度信号序列vc(n)，利用采样频率Fs计算对应的时刻序列tc(n)，tc(n)＝(0：1：length(vc(n))
‑
1)/Fs，式中length(
·
)表示求长度运算符，进而利用vc(n)、tc(n)和tang(n)根据步...

【专利技术属性】
技术研发人员：区瑞坚，孙鹏，郭美娜，
申请(专利权)人：苏州微著设备诊断技术有限公司，
类型：发明
国别省市：