一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法技术

本发明专利技术公开了一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，包括以下步骤：首先安装并利用加速度传感器获取待测滚珠丝杠副的振动信号；然后构建陷波器对振动信号进行降噪处理；接着对振动信号进行EMD分解获得IMF分量，并计算各IMF分量的FFT；之后重构振动信号并对其进行HHT包络解调处理，获得Hilbert谱及Hilbert边际谱；最后提取滚珠丝杠副疲劳失效特征频率，并将其与理论计算值进行对比分析，从而判断滚珠丝杠副是否发生失效。本发明专利技术的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，可有效解决当前滚珠丝杠副疲劳失效诊断困难、效率低下、准确性差的问题。

【技术实现步骤摘要】
一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法
本专利技术属于滚珠丝杠副故障诊断领域，特别是一种滚珠丝杠副疲劳失效诊断方法。
技术介绍
滚珠丝杠副因其传动效率高、跟随灵敏、磨损小、零间隙高刚度传动等特点，成为高性能数控机床不可或缺的功能部件，始终占据直线运动应用领域的绝大部分市场。然而在实际加工工作中，设备往往存在长时间处于负载过重、润滑不良、维护保养缺乏等现象，使得滚珠丝杠副产生以接触疲劳点蚀为主要形式的疲劳失效。由此产生的周期性冲击振动，使得滚珠丝杠副精度降低，一定程度时致使整个传动系统发生故障而失效，造成巨大经济损失和安全隐患。因此，在实际生产加工中，如何准确高效地对滚珠丝杠副疲劳失效进行诊断对预防设备故障、提高设备工作寿命具有重大意义。迄今为止，鲜有对滚珠丝杠副疲劳失效诊断理论和试验的深入研究，这也使得在实际加工工作中，滚珠丝杠副疲劳失效的判断依据停留在工人师傅肉眼辨别或者加工设备出现故障。因此急需一种简单易行的方法来实现对滚珠丝杠副疲劳失效的高准确性、高效率的诊断。
技术实现思路
本专利技术所解决的技术问题在于提供一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，解决当前滚珠丝杠副疲劳失效诊断困难、效率低下、准确性差的问题。实现本专利技术目的的技术解决方案为：一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，包括以下步骤：步骤1、安装并利用加速度传感器获取待测滚珠丝杠副的振动信号；步骤2、构建陷波器对振动信号进行降噪处理；步骤3、对降噪处理后的振动信号进行EMD分解获得IMF分量，并计算各IMF分量的FFT；步骤4、根据步骤3获得的IMF分量重构振动信号，并对重构信号进行HHT包络解调处理，获得Hilbert谱及Hilbert边际谱；步骤5、从Hilbert谱及Hilbert边际谱中提取滚珠丝杠副疲劳失效特征频率，并将其与所述特征频率的理论推导值进行对比分析，获得滚珠丝杠副疲劳失效诊断结果。本专利技术与现有技术相比，其显著优点为：1)本专利技术的滚珠丝杠副疲劳失效诊断方法，可有效解决当前滚珠丝杠副疲劳失效诊断困难、效率低下、准确性差的问题；2)本专利技术通过分析提取得到的滚珠丝杠副振动信号特征来诊断是否发生疲劳失效，具有更高的准确性和效率；3)本专利技术采用EMD与HHT处理振动信号，弥补了传统方式不能完全意义上进行非线性非平稳信号处理的缺陷；4)本专利技术诊断滚珠丝杠副疲劳失效方法操作简便，适用性好，科学合理。下面结合附图对本专利技术作进一步详细描述。附图说明图1为本专利技术基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法流程图。图2为本专利技术的振动信号经EMD分解得到的有效IMF分量时域图，其中图(a)为有效分量IMF5的时域图，图(b)为有效分量IMF6的时域图，图(c)为有效分量IMF7的时域图。图3为本专利技术的振动信号经EMD分解得到的有效IMF分量经FFT变换后的频域图，其中图(a)为有效分量IMF5经FFT变换后的频域图，图(b)为有效分量IMF6经FFT变换后的频域图，图(c)为有效分量IMF7经FFT变换后的频域图。图4为本专利技术的重构信号的Hilbert边际谱图。具体实施方式结合图1，本专利技术的一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，包括以下步骤：步骤1、安装并利用加速度传感器获取待测滚珠丝杠副的振动信号；步骤2、构建陷波器对振动信号进行降噪处理；步骤3、对降噪处理后的振动信号进行EMD分解获得IMF分量，并计算各IMF分量的FFT；步骤4、根据步骤3获得的IMF分量重构振动信号，并对重构信号进行HHT包络解调处理，获得Hilbert谱及Hilbert边际谱；步骤5、从Hilbert谱及Hilbert边际谱中提取滚珠丝杠副疲劳失效特征频率，并将其与所述特征频率的理论推导值进行对比分析，获得滚珠丝杠副疲劳失效诊断结果。进一步地，步骤1所述安装并利用加速度传感器获取待测滚珠丝杠副的振动信号，包括确定加速度传感器的类型、布置安装在待测滚珠丝杠副上的位置以及采样频率。进一步优选地，步骤1所述安装并利用加速度传感器获取待测滚珠丝杠副的振动信号，具体为：(1)选用的加速度传感器类型为三轴加速度传感器；(2)加速度传感器安装位置为待测滚珠丝杠副螺母的垂直于丝杠轴向的径向处；(3)振动信号的采样频率为fs，取2f0≤fs≤fa，其中f0为振动信号最高频率，fa为数据采集系统上限，单位均为Hz。进一步优选地，步骤2所述构建陷波器对振动信号进行降噪处理，陷波器的公式为：X(z)＝X0(z)H(z)式中，x0(t)为输入陷波器的原始时域信号序列；x(t)为陷波器降噪后输出的时域信号序列即降噪后的振动信号；v(t)为数字滤波器的脉冲响应函数；X(z)和X0(z)分别为输入和输出陷波器的频域信号序列；H(z)为陷波器传递函数；z为Z变换中的复变量；ω0为要滤除的频率，单位为Hz；b0为增益系数；r为与陷波器带宽有关的量，取值为0≤r≤1。示例性优选地，ω0＝50Hz。进一步地，步骤3所述对降噪处理后的振动信号进行EMD分解获得IMF分量，并计算各IMF分量的FFT，具体为：步骤3-1、确定降噪后的振动信号x(t)所有局部极大值点和局部极小值点；步骤3-2、利用三次样条函数分别对所有局部极大值点和局部极小值点进行差值运算，拟合出上包络线和下包络线，并计算上下包络的均值：式中，θ(t)为上包络线，γ(t)为下包络线，m1(t)为均值包络线；步骤3-3、令h1＝x(t)-m1(t)，判断h1是否满足IMF分量条件，若满足条件，则记h1为x(t)的第1个IMF分量；若不满足条件，则将其作为x(t)并重复步骤3-1～步骤3-2直至h1k满足IMF分量条件，则h1k为第1个IMF分量，记为c1；其中h1k为重复k次后的IMF分量，具体为：h1k＝h1(k-1)(t)-m1k(t)式中，h1(k-1)(t)为重复k-1次后的IMF分量，m1k(t)为h1(k-1)的上、下包络线均值；步骤3-4、从x(t)中减去c1，获得第一次提取IMF分量的残余函数r1＝x(t)-c1，并将r1作为x(t)重复步骤3-1～步骤3-3，获得x(t)的第2个IMF分量c2；以此类推，当rn为单调函数不能再提取IMF分量时，循环结束，由此获得n个IMF分量；振动信号x(t)由IMF分量表示为：式中，rn为第n次提取IMF分量的残余函数，代表信号的平均趋势，ci为x(t)的第i个IMF分量；步骤3-5、采用连续两个IMF分量的标准差SD作为再次判断提取IMF分量终止的判据：若SD≤SD0，则结束IMF分量提取；若SD＞SD0，则将残余函数作为x(t)并返回执行步骤3-1继续IMF分量提取；SD0通常取值为0.2～0.3；T表示信号长度；步骤3-6、通过FFT将各IMF分量的时域序列转换为频域序列：ci(t)＝cia(t)+cib(t)式中，ci(t)为时域下各IMF分量，cia(t)为偶数时域序列，cib(t)为奇数时域序列，Ci(e)为频域下各IMF分量；Cia(e)为偶数频域序列；Cib(e)为奇数频域序列；T为IMF分量的个数；为偶数时域序列对应的旋转因子，为奇数时域序列对应的旋转因子，为旋转因子。进一步地，步骤4所述根据步骤3获得的IMF分量重构振动信号，并对重构信本文档来自技高网...

【技术保护点】
1.一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，其特征在于，包括以下步骤：步骤1、安装并利用加速度传感器获取待测滚珠丝杠副的振动信号；步骤2、构建陷波器对振动信号进行降噪处理；步骤3、对降噪处理后的振动信号进行EMD分解获得IMF分量，并计算各IMF分量的FFT；步骤4、根据步骤3获得的IMF分量重构振动信号，并对重构信号进行HHT包络解调处理，获得Hilbert谱及Hilbert边际谱；步骤5、从Hilbert谱及Hilbert边际谱中提取滚珠丝杠副疲劳失效特征频率，并将其与所述特征频率的理论推导值进行对比分析，获得滚珠丝杠副疲劳失效诊断结果。

【技术特征摘要】
1.一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，其特征在于，包括以下步骤：步骤1、安装并利用加速度传感器获取待测滚珠丝杠副的振动信号；步骤2、构建陷波器对振动信号进行降噪处理；步骤3、对降噪处理后的振动信号进行EMD分解获得IMF分量，并计算各IMF分量的FFT；步骤4、根据步骤3获得的IMF分量重构振动信号，并对重构信号进行HHT包络解调处理，获得Hilbert谱及Hilbert边际谱；步骤5、从Hilbert谱及Hilbert边际谱中提取滚珠丝杠副疲劳失效特征频率，并将其与所述特征频率的理论推导值进行对比分析，获得滚珠丝杠副疲劳失效诊断结果。2.根据权利要求1所述的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，其特征在于，步骤1所述安装并利用加速度传感器获取待测滚珠丝杠副的振动信号，包括确定加速度传感器的类型、布置安装在待测滚珠丝杠副上的位置以及采样频率。3.根据权利要求2所述的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，其特征在于，步骤1所述安装并利用加速度传感器获取待测滚珠丝杠副的振动信号，具体为：(1)选用的加速度传感器类型为三轴加速度传感器；(2)加速度传感器安装位置为待测滚珠丝杠副螺母的垂直于丝杠轴向的径向处；(3)振动信号的采样频率为fs，取2f0≤fs≤fa，其中f0为振动信号最高频率，fa为数据采集系统上限，单位均为Hz。4.根据权利要求1所述的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，其特征在于，步骤2所述构建陷波器对振动信号进行降噪处理，陷波器的公式为：X(z)＝X0(z)H(z)式中，x0(t)为输入陷波器的原始时域信号序列；x(t)为陷波器降噪后输出的时域信号序列即降噪后的振动信号；v(t)为数字滤波器的脉冲响应函数；X(z)和X0(z)分别为输入和输出陷波器的频域信号序列；H(z)为陷波器传递函数；z为Z变换中的复变量；ω0为要滤除的频率，单位为Hz；b0为增益系数；r为与陷波器带宽有关的量，取值为0≤r≤1。5.根据权利要求4所述的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，其特征在于，所述ω0＝50Hz。6.根据权利要求4所述的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法，其特征在于，步骤3所述对降噪处理后的振动信号进行EMD分解获得IMF分量，并计算各IMF分量的FFT，具体为：步骤3-1、确定降噪后的振动信号x(t)所有局部极大值点和局部极小值点；步骤3-2、利用三次样条函数分别对所有局部极大值点和局部极小值点进行差值运算，拟合出上包络线和下包络线，并计算上下包络的均值：式中，θ(t)为上包络线，γ(t)为下包络线，m1(t)为均值包络线；步骤3-3、令h1＝x(t)-m1(t)，判断h1是否满足IMF分量条件，若满足条件，则记h1为x(t)的第1个IMF分量；若不满足条件，则将其作为x(t)并重复步骤3-1～步骤3-2直至h1k满足IMF分量条件，则h1k为第1个IMF分量，记为c1；其中h1k为重复k次后的IMF分量，具体为：h1k＝h1(k-1)(t)-m1k(t)式中，h1(k-1)(t)为重复k-1次后的IMF分量，m1k(t)为h1(k-1)的上、下包络线均值；步骤3-4、从x(t)中减去c1，获得第一次提取IMF分量的残余函数r1＝x(t)-c1，并将r1作为x(t)重复步骤3-1～步骤3-3，获得x(t)的第2个IMF分量c2；以此类推，当rn为单调函数不能再提取IMF分量时，循环结束，由此获得n个IMF分量；振动信号x(t)由IMF分量表示为：式中，rn为第n次提取IMF分量的残余函数，代表信号的平均趋势，ci为x(t)的第i个IMF分量；步骤3-5、采用连续两个IMF分量的标准差SD作为再次判断提取IMF分量终止的判据：若SD≤SD0，则结束IMF分量提取；若SD＞SD0，则将残余函数作为x(t)并返回执行步骤3-1继续IMF分量提取；SD0通常取值为0.2～0.3；T表示信号长度；步骤3-6、通过FFT将各IMF分量的时域序列转换为频域序列：ci(t)＝cia(t)+cib(t)式中，ci(t)为时...

【专利技术属性】
技术研发人员：周长光，刘迪一，欧屹，周华西，冯虎田，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏,32