基于经验模态分解与倒频谱的摆线锥齿轮故障诊断方法
【专利摘要】基于经验模态分解与倒频谱的摆线锥齿轮故障诊断方法，属于故障诊断【技术领域】，该方法包括：1、利用加速度传感器对摆线锥齿轮副进行测量，采集加速度振动信号作为待分析信号；2、将采集的信号导入Matlab中，得到原始信号，利用经验模态分解(EMD)方法将原始信号分解为一系列固有模态函数(IMF)分量；3、对前几阶固有模态函数分量进行倒频谱分析，得到其幅值倒频谱；4、采用Matlab软件的绘图工具绘制出幅值倒频谱图，根据倒频谱图中幅值的分布，提取故障特征信息。本发明方能够运用于摆线锥齿轮的故障诊断中，并能准确地提取出故障特征频率。本发明为摆线锥齿轮的故障诊断提供了新方法，且为其他旋转机械故障诊断技术提供了有效借鉴。
【专利说明】基于经验模态分解与倒频谱的摆线锥齿轮故障诊断方法
【技术领域】
[0001]本发明属于齿轮故障诊断【技术领域】，特别涉及一种基于经验模态分解与倒频谱的摆线锥齿轮故障诊断方法。
【背景技术】
[0002]摆线锥齿轮作为重型车铣复合数控机床重载摆角铣头的核心传动部件，由于重载摆角铣头的加工对象为大型复杂曲面，因此摆线锥齿轮经常在振动、冲击、变载荷、频繁换向等非稳态工况下工作，会产生点蚀、磨损、崩裂甚至断齿等故障。故障状态又直接影响了摆角铣头的性能，导致数控机床加工出的零件不合格或整台机床停运，甚至产生连锁反应，严重影响企业的安全生产和经济效益。因此，研究摆线锥齿轮的故障诊断技术及方法，尽早发现并及时消除故障对于保障机床设备安全、避免事故和巨额经济损失及提高设备使用性能，具有非常重要的现实意义。
[0003]为有效提高重载摆角铣头的加工性能及摆线锥齿轮的寿命和可靠性，有必要深入研究摆线锥齿轮的故障诊断技术及方法。本发明通过分析摆线锥齿轮故障振动信号的时域及频域特征，根据经验模态分解(EMD)方法及倒频谱分析方法，提出了一种基于经验模态分解与倒频谱的摆线锥齿轮故障诊断方法，该发明方法能够准确提取出摆线锥齿轮的故障信息，实现摆线锥齿轮的故障诊断。本发明方法不仅能达到摆线锥齿轮状态监测和故障诊断的目的，且为其他旋转机械复合故障诊断技术提供了有效借鉴。
【发明内容】
[0004]本发明的目的是提供一种基于经验模态分解与倒频谱的摆线锥齿轮故障诊断方法，通过对摆线锥齿轮振动信号进行分析和处理，准确地提取出摆线锥齿轮故障特征。本发明为摆线锥齿轮的故障诊断和特征提取提供了新方法，且为其他旋转机械复合故障诊断技术提供了有效借鉴。
[0005]本发明是采用以下技术手段实现的:
[0006]基于经验模态分解与倒频谱的摆线锥齿轮故障诊断方法，其特征在于:具体步骤如下:
[0007]SI利用加速度传感器对摆线锥齿轮副进行测量，采集该种齿轮的加速度振动信号作为待分析信号；
[0008]S2将采集的待分析信号导入Matlab软件中，得到原始信号x (t)，采用经验模态分解法将原始信号进行分解，得到一系列固有模态函数IMF分量；
[0009]S2.1.确定原始信号x(t)的所有局部极值点，用三次样条曲线将所有的极大值点连接起来形成X (t)的上包络线，同样连接所有的极小值点形成X (t)的下包络线，上下包络线应该包络所有的数据点；记上下包络线的平均值为Hl1U),并计算X(t)和!111(0的差值Mt),得到:
[0010]Ii1 (t) = X (t) -Hi1 (t)；[0011]S2.2.对于不同的原始信号x(t)，hjt)可能满足固有模态函数MF的条件，也可能不满足；若不满足IMF条件，此时将Ill (t)作为原始信号，重复步骤S2.1得到Ill (t)的分解步骤，得到:
[0012]hn(t) = K1 (t) -mn (t)
[0013]其中，mn(t)为匕⑴的上、下包络线均值；
[0014]S2.3.如果hn (t)不满足MF的条件，则分解继续，重复上述步骤2.1中的分解k次’得..[0015]hlk(t) = hi Q5^1) (t)-mlk(t)
[0016]S2.4.判定hlk(t)是否为一个MF分量，必须要有分解终止准则，它可以定义为连续两个分解结果之间的标准差SD值:
【权利要求】
1.基于经验模态分解与倒频谱的摆线锥齿轮故障诊断方法，其特征在于:具体步骤如下: Si利用加速度传感器对摆线锥齿轮副进行测量，采集该种齿轮的加速度振动信号作为待分析信号； S2将采集的待分析信号导入Matlab软件中，得到原始信号x (t)，采用经验模态分解法将原始信号进行分解，得到一系列固有模态函数IMF分量； S2.1.确定原始信号x(t)的所有局部极值点，用三次样条曲线将所有的极大值点连接起来形成x(t)的上包络线，同样连接所有的极小值点形成x(t)的下包络线，上下包络线应该包络所有的数据点；记上下包络线的平均值为Hl1U),并计算X(t)和!111(0的差值Ii1 (t)，得到:
hi (t) = X (t) -1ii1 (t)； S2.2.对于不同的原始信号XUhh1(t)可能满足固有模态函数IMF的条件，也可能不满足；若不满足MF条件，此时将Ill (t)作为原始信号，重复步骤S2.1得到Ill (t)的分解步骤，得到:
hn(t) = hi (t)-mn(t) 其中，mn(t) Sh1U)的上、下包络线均值； S2.3.如果hn(t)不满足MF的条件，则分解继续，重复上述步骤2.1中的分解k次，得:
hik(t) = Ii1(H) (t)-mlk(t) S2.4.判定hlk(t)是否为一个MF分量，必须要有分解终止准则，它可以定义为连续两个分解结果之间的标准差SD值:
【文档编号】G01M13/02GK104006961SQ201410181013
【公开日】2014年8月27日 申请日期:2014年4月29日 优先权日:2014年4月29日
【发明者】杨文通, 罗兵, 刘志峰, 湛承鹏 申请人:北京工业大学